compute the analytic structure of the Feynman Diagrams of a model - Maple Programming Help

Online Help

All Products    Maple    MapleSim

Home : Support : Online Help : Physics : Physics/FeynmanDiagrams

Physics[FeynmanDiagrams] - compute the analytic structure of the Feynman Diagrams of a model

Calling Sequence

FeynmanDiagrams(InteractionLagrangian, vertices = ..., numberoflegs = ..., numberofloops = ..., externallegs = ..., excludepropagators = ..., crossedpropagators = ..., includetadpoles = ..., fields = ...)




the interaction terms of a Lagrangian depending on quantum fields

vertices = ...


(optional) the same as points = ..., list of points at which the Lagrangian is evaluated; these are the vertices of the corresponding Feynman diagrams; default value is 1 vertex

numberoflegs = ...


(optional) the number of external legs of the Feynman diagrams, these legs are represented in the output as _NP(...) normal ordered products; default value is all possible cases

numberofloops = ...


(optional) integer or set of integers; the number of loops of the Feynman diagrams: each line joining two vertices of the loop is a propagator, represented in the output as a list of two fields

externallegs = ...


(optional) the same as normalproducts = ..., the right-hand side is a function _NP(...), or a set of them, where each _NP(...) contains fields or field names as arguments, which are the normal ordered products of unpaired fields expected in the output; to these unpaired fields are associated the external legs of the related Feynman diagrams

excludepropagators = ...


(optional) a list of two fields representing a propagator, or a set of these lists, indicating that terms involving these propagators are to be excluded from the output; the right-hand side can also be one of the keywords free or cross

crossedpropagators = ...


(optional) it can be true or false (default) to indicate that crossed propagators (related to crossed terms in the quadratic part of the Lagrangian) are not zero and so must be taken into account

includetadpoles = ...


(optional) it can be the false (default) or true to include or discard terms containing tadpoles, that is, lines that start and end in the same vertex.

fields = ...


(optional) either a function or its name, or a set or list of them; indicates the fields of the InteractionLagrangian;



A scattering matrix S relates the initial and final states of an interacting system. In an N-dimensional spacetime with coordinates X, S can be written as:



where I is the imaginary unit, L is the interaction Lagrangian, written in terms of quantum fields depending on the spacetime coordinates X. The T symbol means time-ordered. For the terminology used in this page, see chapter IV of Bogoliubov, N.N., and Shirkov, D.V. Quantum Fields. Benjamin Cummings, 1982. This exponential can be expanded as a sum of terms.



Each Sn contains integrals of time-ordered products of the interaction Lagrangians evaluated at different points (vertices of the Feynman Diagram). For example, the second order term is of the form:



where X and Y represent different points in spacetime.


Using Wick's theorems, the FeynmanDiagrams command computes the expanded form of these time-ordered products of interaction Lagrangians in the integrands of each Sn.


The expansion of each Snconsists of a sum of terms, each of which is a product of certain pairings of the operator field functions entering L with the normal product of the remaining unpaired operators of free fields. The pairings between fields are the propagators, the Green functions of the associated free fields.


To each of the terms of the expansion of Sn corresponds a Feynman diagram with n vertices (the number of interaction Lagrangians in Sn). Within the diagram, each of the pairings of fields has a corresponding internal line, and each field entering the normal products of the unpaired field operators has a corresponding external leg. The output of the FeynmanDiagrams command is thus a sum of terms that has a one-to-one correspondence to the sequence of Feynman diagrams for the given interaction Lagrangian.


NOTE: the current implementation of FeynmanDiagrams can expand up to S3 (products of up to three interaction Lagrangians), and therefore can produce the analytic representation of Feynman diagrams with no more than three vertices.


In the Maple Standard GUI, to have the output of this help page use textbook notation, open it as a worksheet (see icon in the toolbar) and execute the input lines below after setting Physics[Setup](mathematicalnotation = true);

Load the package and set three coordinate systems to work in.






Default differentiation variables for d_, D_ and dAlembertian are: X=x1,x2,x3,x4

Systems of spacetime Coordinates are: X=x1,x2,x3,x4,Y=y1,y2,y3,y4,Z=z1,z2,z3,z4






The expressions entering S1 (only one vertex and so one evaluation point), representing the connected Feynman graphs for this interaction Lagrangian and discarding terms with tadpoles, is:




To include terms with tadpoles use the option includetadpoles = true or just includetadpoles




You can filter these result to obtain only the tree-level terms (see the Options subsection of the Description).




After using FeynmanDiagrams once, the fields phi and eta get automatically set as quantum operators, so you do not need to set them as such again. Alternatively, you can always set them as quantumoperators (before or after calling FeynmanDiagrams) using the Setup command.


* Partial match of 'quantumop' against keyword 'quantumoperators'



The structure of S2 (that is, the expanded form of the time-ordered product of interaction Lagrangians entering the integrand of the second term in the expansion of the scattering matrix) is:




In the result above, the tree-level terms involving propagators are associated with connected, one-particle-reducible Feynman graphs, and can be obtained by specifying numberofloops=0.




NOTE: the current implementation of FeynmanDiagrams can expand products of up to three interaction Lagrangians, so a maximum of three spacetime points are allowed, and therefore produces the analytic representation of Feynman diagrams with no more than three vertices.

Compute only the terms corresponding to Feynman diagrams with two external legs.




Compute only the terms corresponding to Feynman diagrams with two external legs corresponding to the field phi, regardless of the vertex to which they are attached (it could be X or Y).




The same computation but including tadpoles




Compute only the terms with two external legs corresponding to the normal products etax1,x2,x3,x4,etay1,y2,y3,y4.




Compute the terms leading to Feynman graphs with two loops.




In S2, among the terms resulting from expanding the time-ordered product of two interaction Lagrangians (Feynman diagrams with two vertices), there are no terms involving three loops.




Compute the terms entering the integrand of S3, having four external legs (normal ordered products involving four fields); in this example, these are of the form phi,eta, where each field is applied at X, Y, or Z.




Following is an example of an interaction Lagrangian with a spinor-electromagnetic vertex. First, define the vector field Amu.


Defined objects with tensor properties



Check the type of letter used to represent spinor indices (or, change this setting according to your preference using Setup) and set A as a quantum operator.


* Partial match of 'op' against keyword 'quantumoperators'



So, to represent spinor indices, use lowercase Latin letters, and the spinor fields will be represented by any name prefixed by Q, and the matrix indices related to the matrix product involving Dirac matrices are written explicitly in the interaction Lagrangian.




Compute the expanded form of the integrand entering S2, only the terms with 4 external legs, so all the corresponding graphs have no loops.




The terms entering the integrand of S2 and that only have 2 external legs, so all the corresponding graphs have one loop.




Compute the expanded form of the integrand entering S3 for this model.