matlab费曼编码输入,多点格林函数数值积分(费曼参数积分)的程序分析及应用

0引言在粒子物理中研究物质基本结构和基本相互作用中所涉及的物理量是由计算相应过程发生的跃迁几率振幅,一般称为散射矩阵或S矩阵来得到的.跃迁矩阵元的计算是在微扰论的基础上,根据费曼规则一阶一阶的计算的.从第二阶开始,费曼图包括了复杂的圈图,其中要对某个(或某几个)独立的动量进行积分,这些积分可能是发散的,需要用重整化处理.重整化方案已经在理论上比较完美地解决了,但另一方面,还面临着非常复杂的费曼参数积分问题.圈图积分的一般表达式可写为:TN1…N=(2)4-Di2dDqq1…qN[q2-m21][(q+k1)2-m22]…[(q+kN-1)2-m2N]对动量的费曼积分展开式见文献[3].可见,计算跃迁矩阵元的核心问题就是计算不同形式的四动量标量积分,统称格林函数.目前面临的复杂的多重积分,随着一圈图顶角的增多,即使有基本公式,手工推导的工作量也是极大的,如增加到二圈图,计算的复杂度更是不可想象的,更可怕的是出现奇点的情况,这使计算跃迁矩阵变得十分困难.解决上述问题只有借助计算机程序来避免高阶计算时手工推导的繁琐,这种程序设计工作的工作量也是极大的.目前处理此问题已经取得了一定的成绩,可以一定程度地避免大部分情况下的奇点的影响,得出四点以下格林函数有限部分的数值结果.1软件包LoopTools1LoopTools简介LoopTools是一个用于数值计算高能物理一圈图费曼参数积分的软件包,由ThomasHahn于2001年在FF软件包基础上开发完成.它提供完整的一到四点格林函数的计算程序,并兼容Fortran,C和Mathematica多种语言.1.2LoopTools的获取与安装LoopTools以压缩的tar文件形式出现,目前最新版本为LoopTools21.tar.gz,可在http://www.feynarts.de/looptools获得.LoopTools的安装使命完整写在自带的makefile文件中,只需将压缩文件解压到指定目录下后直接运行make命令编译即可.1.3应用LoopTools的具体操作方法见文献[3].2LoopTools的工作原理LoopTools的工作原理比较简单,在仔细研究了./ff,./fortran,./check几个目录中的子过程文件后可以发现LoopTools的fortran源程序的运行原则,具体内容见文献[3].3具体应用实例和结论为了验证LoopTools软件包的正确性,下面将考虑一具体的散射一圈图(见图1),代入数据后将LoopTools计算的结果与手工推导的结果做比较说明.图1散射一圈图Fig1TheboxdiagramwhichcontributestoP0-P0mixingM=ig44i,jijd4q(2)4l(p3)v(1-5)h(p1)l(p2)(1-5)h(p4)q(q-p1)v1(q2-M2W)(q2-m2i)[(q-p1)2-M2W][(q-p1)2-m2j]公式推导动量积分部分I=d4qq(q-p1)v(q2-M2W)(q2-m2i)[(q-p1)2-M2W][(q-p1)2-m2j]=1(M2W-m2i)(M2W-m2j)d4qq(q-p1)v{1(q2-M2W)[(q-p1)2-M2W]+1(q2-m2i)[(q-p1)2-m2j]-1(q2-M2W)[(q-p1)2-m2j]-1(q2-m2i)[(q-p1)2-M2W]}设:Iij=d4qq(q-p1)v(q2-m2i)[(q-p1)2-m2j]=d4qq(q-p1)v10dx1[(q-p1x)2+p21(x-x2)-m2i(1-