An Extended Version of the Simple and Fast Feynman Diagram Generator EasyFeynDiag.
Bo Xiao,Hao Wang
DOI: https://doi.org/10.1016/j.cpc.2014.08.025
IF: 4.717
2014-01-01
Computer Physics Communications
Abstract:The main improvement in this version of EasyFeynDiag is to prohibit equivalent Feynman diagrams to two-loop order (while the previous version is to one-loop). The symmetry factor calculation with equivalent diagrams prohibition to two-loop is a nontrivial problem, which is well solved in the new program. The code of the program is also reorganized, hopefully in a way more convenient for people to read.New version program summaryProgram title: EasyFeynDiagCatalogue identifier: AEPD_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEPD_v2_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 51148No. of bytes in distributed program, including test data, etc.: 55,295Distribution format: tar.gzProgramming language: C.Computer: Platforms on which an ANSI C compiler is available.Operating system: Operating Systems on which an ANSI C compiler is available.RAM: 236,000 bytes for a typical two-loop order uū→tt̄ process in the QCD model.Classification: 4.4.Does the new version supersede the previous version?: The new version generates Feynman diagrams better, but is more time consuming, than the previous version. Overall, the new version supersedes the previous version.Nature of problem:Automatic generation of Feynman diagrams from given interactions, scattering process, and number of loops.Solution method:The algorithm of EasyFeynDiag is translated from the perturbative expansion formula of the S-matrix [1]. It follows an ordered iterative traversing procedure to find all the possible combinations of the fields and creation/annihilation operators.Two selection rules are adopted in turn, to suppress equivalent Feynman diagrams originating from vertices relabeling, and from multiple identical particles of an interaction term, respectively. This simple algorithm prohibits equivalent diagrams to one-loop order.To further prohibit equivalent diagrams to higher order, more careful treating of the identical particles is required. By introducing one additional “order rule” for the identical particles and one additional “selection rule” for identical vertices, the equivalent diagrams are prohibited to two-loop order. Those two rules are explained briefly below:1. order rule for identical particlesWhen identical particles from the same vertex are to connect to other legs, those connections must be ordered: the second identical particle must connect to a leg equal or behind the leg connected by the first identical particle, and so on the Nth to the (N−1)th identical particle. At the same time, proper symmetry factors need to be generated to count for this ordering.2. selection rule for identical verticesWhen m identical particles from the same vertex have connected to m new vertices that are of the same interaction type, these m vertices are defined as identical vertices. If later another particle is to connect to one of these identical vertices, it is only permitted to connect to the most front one, the other (m−1) connections are represented by a factor m.Reasons for new version:High energy calculations to two-loop order are becoming more and more common, e.g. [2]. To introduce, into EasyFeynDiag, the ability to prohibit equivalent diagrams to two-loop order would make it more attractive to researchers in High Energy Physics.Summary of revisions:1.the code is organized in a totally new way, more convenient for people to read and modify;2.the algorithm for Feynman diagram generation in the program is extended to a version that prohibits equivalent diagrams to two-loop order (the previous version is to one-loop order).Restrictions:1.There are default values for the maximum external particles and internal vertices for a scattering process, and maximum legs or leg types for an interaction term. Users need to modify the corresponding values in the code to get larger limits.2.Equivalent Feynman diagrams arise at three-loop order or higher, although this fact does not hurt the correctness of the Feynman amplitudes to be obtained from those diagrams.Unusual features:Allow for arbitrary physical model and arbitrary number of loops. Symmetry factors are naturally generated. Small and fast.Running time:About 0.139 s to generate all the Feynman diagrams of a typical two-loop order uū→tt̄ process in the QCD model.References:[1]Bo Xiao, Hao Wang, Shou-hua Zhu, A simple algorithm for automatic Feynman diagram generation, Comput. Phys. Commn. 184 (2013) 1966.[2]Peter Barnreuther, Michal Czakon, Alexander Mitov, Percent level precision physics at the tevatron: first genuine NNLO QCD corrections to qq̄→tt̄+X, Phys. Rev. Lett. 109, 132001 (2012).