the astrophysicist subrahmanyan chandrasekhar once referred to black holes as `` the most perfect macroscopic objects there are in the universe : the only elements in their construction are our concepts of space and time.''@xcite this special issue of journal of physics a is devoted to a completely different topic , the quantum - mechanical amplitudes for the scattering of relativistic particles , in which a revolution in our understanding has been taking place recently .
however , scattering amplitudes represent a kind of short - distance complement to the perfection chandrasekhar saw in black holes .
we might consider them to be the most perfect microscopic structures in the universe .
black holes are gigantic and ( almost ) eternal .
scattering amplitudes , in contrast , portray events that wink in and out of existence much faster than the blink of an eye .
they describe processes at the shortest distance scales that can be probed in the laboratory . when gravitons scatter , interacting via einstein s field equations , again in some sense the only elements in the construction are our concepts of space and time  although the boundary conditions at infinity can be more detailed than for black holes , involving multiple plane - wave ripples ( the gravitons ) propagating in various directions through minkowski space - time . when other relativistic particles scatter , as in non - abelian gauge theory , more information has to be provided , such as the gauge group , the coupling constant , the matter content of the theory , the particle masses , and the precise set of particle types and spins being scattered
. however , in the ultra - relativistic limit , if there are only gauge interactions , much of the group - theoretical information can be stripped from the amplitudes , and certain color - ordered `` primitive '' amplitudes can be defined which are quite universal .
these amplitudes have intricate analytic properties , which have been teased out over the decades , but especially in the past few years . a variety of new , efficient techniques have been found to construct amplitudes , and hidden symmetries have been identified which illustrate the amplitudes `` perfection . ''    in the real universe , black holes are not found in an absolute vacuum . typically they are surrounded by accretion disks of gas and dust , which muddy their perfection somewhat , but also enrich the range of phenomena involving them , and make them much easier to observe , albeit indirectly .
similarly , scattering amplitudes for quantum chromodynamics , or qcd , the most relativistic non - abelian gauge theory we can study in the laboratory , are shrouded from direct view by confinement .
the quarks and gluons that theorists scatter mentally are never seen experimentally . in the initial state of a hadronic collision
they are bound into a proton or other hadron . in the final state
they emerge as collimated jets of particles . despite this fact , quantities that are sensitive only to very short distances , so - called infrared - safe observables such as jet production cross sections ,
can be computed reliably in terms of quarks and gluons , in a systematic expansion in the strong coupling @xmath1 .
if we could go to asymptotically high energies , so that @xmath1 were infinitesimal , it would be enough to consider just the leading order in this expansion , which is generated by the born approximation , or tree - level amplitudes .
the tree - level amplitudes of qcd have even more `` perfection '' than the generic , loop - level amplitudes , because they coincide with the tree amplitudes of a much more symmetric theory , @xmath0 supersymmetric yang - mills theory ( @xmath0 sym ) .
this theory , which is the subject of many of the articles in this special issue , has the maximum amount of supersymmetry possible in a gauge theory .
its four supercharges can be used to transform the helicity @xmath2 gluon state , by @xmath3 unit of helicity at a time , all the way to the helicity @xmath4 cpt conjugate gluon state .
all states in the theory belong to a single supermultiplet , transforming in the adjoint representation of the gauge group
. therefore all interactions are related by supersymmetry to the triple - gluon vertex , and the theory has a single dimensionless coupling @xmath5 . the @xmath6 function for @xmath0 sym vanishes exactly for any gauge group , @xmath7 , so that the theory is exactly conformally invariant at the quantum level , as well as classically .
( the classical conformal invariance of qcd with massless quarks is of course spoiled by its nonvanishing one - loop @xmath6 function , which leads to asymptotic freedom . )    in investigations of scattering amplitudes in @xmath0 sym , the gauge group under consideration is usually su@xmath8 . quite often , the limit @xmath9 is also taken , at a fixed value of the t  hooft coupling , @xmath10 . for @xmath11 ,
perturbation theory can be applied , and for large @xmath12 only planar feynman diagrams contribute . for @xmath13 ,
perturbation theory breaks down , but the ads / cft correspondence  @xcite can be used . in the large-@xmath12 or planar limit of @xmath0 sym , even more remarkable symmetries and relations emerge for scattering amplitudes .
as we will discuss further below , these properties include dual conformal ( super)symmetry , which is part of a larger yangian invariance ; a description of amplitudes in terms of spaces of complex planes ( grassmannians ) ; and in terms of various types of twistors ; and a relation between scattering amplitudes and the expectation values of wilson loops for closed polygons bounded by light - like edges .
many of these properties are intimately connected with a strong - coupling picture of gluon scattering in terms of strings moving in anti - de - sitter space  @xcite .
the exceptional simplicity and numerous hidden symmetries of @xmath0 sym have made it a playground for theorists interested in scattering amplitudes for a couple of decades .
this interest has accelerated rapidly in the past few years . because the tree - level amplitudes of qcd and @xmath0 sym coincide , the early discoveries about qcd helicity amplitudes , motivated by the physics of jet production , were also discoveries about @xmath0 sym .
for example , the parke - taylor formula  @xcite for the maximal - helicity - violating ( mhv ) sequence of @xmath14-gluon amplitudes , which is just a single - term expression , even for an arbitrarily large value of @xmath14 , is also consistent with all the symmetries of @xmath0 sym  @xcite , including some that were only unveiled decades later .
many general properties of tree amplitudes were understood early on , such as their factorization properties  @xcite and the ( @xmath15 or @xmath16 ) supersymmetry ward identities that they obey  @xcite . however , a fuller understanding and exploitation of these properties have come only in recent years .
for instance , the complete solutions of the long - known @xmath0 supersymmetry ward identities , and the corresponding @xmath17 identities in supergravity , have been worked out recently by elvang , freedman and kiermaier  @xcite , and are reviewed here  @xcite .
several years ago , witten  @xcite discovered that the twistor space developed decades earlier by penrose  @xcite gave a natural description of tree - level amplitudes in massless gauge theory .
as discussed in , such amplitudes have a natural description in terms of two - component weyl spinors , a left - handed and a right - handed spinor for each external state .
to go to twistor space , one performs a fourier transform on half the variables , namely the left - handed spinors , leaving the right - handed ones alone .
witten found that gauge theory could be reformulated in terms of a topological string moving in twistor space .
gauge amplitudes turn out to be localized on particular curves in twistor space . in one approach , the curves are intersecting lines : a single line for the simplest mhv amplitudes , a pair of lines for the next - to - maximally - helicity - violating ( nmhv ) amplitudes , and so on .
this version gave rise to the csw or mhv rules for gauge theory  @xcite , whose mhv vertices are a particular off - shell continuation of the parke - taylor formulae . these developments are reviewed here by brandhuber , spence and travaglini  @xcite , and by adamo , bullimore , mason and skinner  @xcite .
the latter reference also discusses another class of twistors , namely momentum twistors , which also provides a natural set of kinematic variables for amplitudes , in that they automatically satisfy the kinematic constraints of momentum conservation and the mass - shell conditions .
finally , the correspondence between wilson loops and scattering amplitudes is described in this article  @xcite from the perspective of twistor space .
britto , cachazo , feng and witten ( bcfw )  @xcite recognized that factorization properties are powerful enough to allow the derivation of recursion relations  @xcite for tree amplitudes , in terms of on - shell lower - point amplitudes , which are evaluated at particular complex momenta .
this observation , reviewed here by brandhuber , spence and travaglini  @xcite , makes essential use of the analyticity , or _
plasticity _ , of amplitudes ; that is , how they vary under smooth deformations of the kinematics .
bcfw embedded an amplitude @xmath18 into a one - complex - parameter family of on - shell amplitudes @xmath19 with shifted complex momenta , and associated the poles of @xmath19 in the @xmath20 plane with factorization of the amplitude onto simpler lower - point amplitudes .
one can also construct a supersymmetrized version of the bcfw recursion relation , in which grassmann variables @xmath21 associated with @xmath0 supersymmetry are shifted along with the momenta  @xcite .
this recursion relation can be solved explicitly to yield all tree amplitudes in @xmath0 sym  @xcite ; the solution and its many symmetries are reviewed here by drummond  @xcite .
scattering amplitudes in massless gauge theories all have quite similar structure at tree level , once the color factors have been stripped off , as discussed in .
on the other hand , the structure of loop amplitudes depends critically on the theory .
the simplicity of one - loop multi - gluon amplitudes in @xmath0 sym played a key role in the development of the unitarity method  @xcite .
this method reconstructs loop amplitudes from their unitarity cuts , again exploiting analyticity . in general at one loop , as explained in this issue by britto  @xcite and ita  @xcite , unitarity cuts are matched against a decomposition of the amplitude in terms of a set of scalar integrals  boxes , triangles , bubbles and ( sometimes ) tadpoles  in order to determine the coefficients of the integrals . in @xmath0 sym ,
the high degree of supersymmetry implies that only boxes have non - vanishing coefficients . using unitarity , an infinite sequence of one - loop amplitudes ( the mhv amplitudes )
could be determined in @xmath0 sym from just the product of two tree - level mhv ( parke - taylor ) amplitudes  @xcite .
later , generalized unitarity was applied to the problem of computing one - loop amplitudes , by extracting additional information from products of three and four tree amplitudes , the so - called triple  @xcite and quadruple  @xcite cuts .
quadruple cuts have four on - shell conditions .
these four equations determine the four components of the ( four - dimensional ) loop momenta completely , up to a ( possible ) two - fold degeneracy .
this realization allows the box coefficients to be computed very simply in terms of the product of four tree amplitudes , glued together as the four corners of a box , and evaluated at each of the two solutions for the loop momenta  @xcite . for @xmath0 sym , there are no other coefficients to determine , so the one - loop problem has been `` reduced to quadrature . ''
( it turns out that maximal @xmath17 supergravity can be proven to have the same `` no - triangle '' property  @xcite , so it too is solved at one loop by the quadruple cuts . )
generalized unitarity can also be applied very effectively at the multi - loop level . in principle
it can be used for any gauge ( or gravitational ) theory . as an example , the two - loop four - gluon scattering amplitudes in qcd have been computed in this way  @xcite . in practice
the method has been pushed the furthest in @xmath0 sym .
the basic techniques of multi - loop generalized unitarity are reviewed here by bern and huang  @xcite , and by carrasco and johansson  @xcite .
returning now to one - loop amplitudes for generic gauge theories , including qcd , the triangle and bubble coefficients have to be determined , as well as certain _
rational parts _ , which have no unitarity cuts in four dimensions .
the triple cuts contain information about the triangle coefficients , but they also receive contributions from boxes . similarly , the ordinary double cuts determine the bubble coefficients , once the contributions of triangles and boxes are removed .
this separation can be done analytically , by making use of the different analytic behavior of the different types of cut integrals , as reviewed by britto  @xcite .
it can also be done numerically , by a suitable sampling of the continuum of loop momenta solving the triple or double cut on - shell conditions  @xcite , as described by ita  @xcite .
these two articles also review methods for computing the rational parts of one - loop amplitudes .
although the rational parts have no unitarity cuts in four dimensions ( by definition ) , they can be constructed via their cuts in dimensional regularization in @xmath22 dimensions .
alternatively , it is possible to construct them recursively in the number of legs  @xcite .    on - shell loop amplitudes in massless non - abelian
gauge theories always contain infrared divergences , due to the exchange of soft gluons or virtual collinear splittings . in qcd
it is conventional to regulate these divergences , as well as the ultraviolet ones , using dimensional regularization . in supersymmetric theories , dimensional reduction  @xcite , or the related four - dimensional helicity scheme  @xcite ,
can be used to keep the number of bosonic and fermionic states the same , and preserve ward identities for ordinary supersymmetry .
the general structure of the infrared divergences is well understood from decades of work in qed as well as qcd  @xcite .
it has been worked out in the most detail in the context of dimensional regularization  @xcite .
the basic picture is shown in .
soft divergences and collinear divergences associated with the amplitude @xmath23 each have a universal form .
they can be factorized from each other and from a hard , short - distance part of the amplitude .
soft divergences , denoted by the blob marked @xmath24 , come from exchange of long - wavelength gluons .
these gluons do not have the resolving power to probe the internal structure of the `` jets '' @xmath25 of virtual collinear particles which capture the collinear divergences .
soft gluons can only see the overall color charges of the jets .
an individual jet function @xmath26 depends on the type of particle @xmath27 , but not on the full amplitude kinematics .
the soft function @xmath24 does not depend on the particle types , but only on their momenta and color quantum numbers .
in general these quantum numbers can be mixed by gluon exchange , so @xmath24 is a matrix in color space .
the hard function @xmath28 has no infrared singularities , but generically depends on the particle types , colors and kinematics .
4.5 truein    in the planar or large-@xmath12 limit , the picture simplifies considerably , to that shown in .
now @xmath23 represents the coefficient of a particular color structure , such as @xmath29 ( assuming that all external states are in the adjoint representation ; see ) . in the planar limit
, individual soft gluons can only connect color - adjacent external partons .
there is no mixing of different color structures at large @xmath12 .
one can absorb the entire soft function @xmath24 into jet functions , which corresponds to breaking up the right - hand side of into @xmath14 wedges .
each wedge is bounded by two hard lines , and is composed of `` half '' of each of the two jet functions , as well as the soft gluons exchanged between them . up to nonsingular terms
, the wedge controlling the infrared divergences represents the square root of the sudakov form factor , which is defined as the amplitude for a color - singlet state to decay into a pair of ( adjoint ) gluons .
4.5 truein    in dimensional regularization , the sudakov form factor obeys a particular differential equation  @xcite , whose source term is called the cusp anomalous dimension @xmath30  @xcite .
there is an additional constant of integration , called @xmath31 .
these two functions ( along with the @xmath6 function in a non - conformal theory ) control the infrared divergences to all orders for any amplitude in a planar massless gauge theory . in planar
@xmath0 sym , the cusp anomalous dimension has been determined to all orders in the coupling using integrability  @xcite .
the bds ansatz  @xcite was built up from this description of the infrared singularities of planar amplitudes , plus the observation for the two- and three - loop four - point amplitudes that the hard function @xmath28 essentially reduces to a constant , independent of the kinematics .
besides the cusp anomalous dimension , another place that integrability certainly enters multi - loop scattering amplitudes is in multi - regge - kinematics , a particular class of high - energy or small - angle limits of @xmath32 particle scatterings .
these kinematics gave the first indication that the bds ansatz for mhv amplitudes has to be corrected beginning at the six - point level  @xcite .
( it was previously argued  @xcite using the properties of wilson loops , that the ansatz should be corrected at two loops for a large value of @xmath14 , unless the correspondence between mhv amplitudes and wilson loops were to break down .
the hexagon wilson loop was then computed and found to differ from the ansatz prediction  @xcite . )
the dynamics of gluons in the transverse plane is given , in the leading - logarithmic approximation , by the hamiltonian for a integrable open spin chain , as reviewed here by bartels , lipatov and prygarin  @xcite .
in planar @xmath0 sym , another infrared regulator is more convenient for many purposes , in particular for exploring the loop - level consequences of dual conformal symmetry , which is not preserved by dimensional regularization .
the features of a recently - developed `` higgs regulator ''  @xcite are reviewed in this issue by henn  @xcite .
for this regulator , vacuum expectation values are given to some of the adjoint scalar fields in the theory , breaking the gauge symmetry in such a way that ( in the planar limit ) the amplitudes are regulated by massive particles circulating around the outside of the loop diagrams .
dual conformal symmetry remains exact in a certain sense .
when the particle masses @xmath33 are taken to be much less than the momentum - invariants for the scattering process , logarithmic divergences develop , the analogs of the @xmath34 infrared poles in dimensional regularization .
there is a hierarchy of simplicity in the scattering amplitudes for various types of gauge theory , as sketched in .
this hierarchy begins to be revealed at one loop .
the outer region of the diagram stands for a generic gauge theory with massive matter fields , and perhaps massive gauge bosons , if the gauge symmetry is spontaneously broken , as in electroweak theory .
one - loop amplitudes in such a theory generically contain tadpole integrals .
one - particle cuts are nontrivial , and are particularly delicate because of external - leg contributions  @xcite .
the cut structure of loop integrals containing massive propagators in the loop is generically somewhat more complicated than the purely massless case .
massive particles in the loop can be unstable , which usually necessitates complex masses . when one enters the `` massless '' ring in , corresponding to massless gauge bosons and matter fields ,
most of these complications vanish , although there are still generically rational parts to compute .
the ring `` sym '' stands for supersymmetric gauge theories .
their one - loop amplitudes can be constructed from four - dimensional unitarity cuts alone , _
i.e. _  there are no non - trivial rational parts  @xcite .
5.2 truein    moving further inward in , we arrive at @xmath0 sym .
as mentioned earlier , at one loop the coefficients of bubble and triangle integrals now vanish , as well as the independent rational parts .
( there are other gauge theories with vanishing bubble and triangle coefficients , at least for their @xmath14-gluon amplitudes  @xcite . )
the theory becomes conformally invariant .
it has been conjectured that the _ leading singularities _  the multi - loop analogs of the quadruple cuts  are sufficient to determine the amplitudes at any loop order  @xcite .
in addition , scattering amplitudes have empirically a predictable , uniform transcendental weight  @xcite .
this weight refers to their construction out of polylogarithms , logarithms , and riemann @xmath35 values .
for example , the finite ( @xmath36 ) terms in one - loop @xmath0 sym amplitudes are of weight two : they contain some terms proportional to the polylogarithm @xmath37 , and others which are products of two logarithms , or proportional to @xmath38 , but they do not contain any terms of lower transcendentality .
( at one loop this result just follows from the absence of bubble and tadpole integrals . ) wherever analytic results are available , this uniform transcendentality property holds . for example , at two loops ( weight four for the finite terms ) it holds for the non - planar , subleading - color terms as well as the planar ones  @xcite .
finally , for planar @xmath0 sym , as mentioned earlier , additional symmetries appear , dual superconformal invariance and an associated yangian symmetry .
the extra symmetries are related to integrability .
they give rise to the prospect of adapting methods developed for other integrable systems to determine the @xmath24-matrix exactly , at least in the large @xmath12 limit .
indeed , as reviewed by drummond  @xcite , an integrand exhibiting manifest yangian invariance has been proposed recently  @xcite for the general @xmath39-loop @xmath14-point amplitude , based on the bcfw recursion relations ; see also ref .
@xcite .    at the level of integrated amplitudes ,
dual conformal invariance has an anomaly due to the need for an infrared regulator .
the anomaly was first understood in terms of wilson loops rather than amplitudes  @xcite ( for which divergences are ultraviolet in nature , rather than infrared ) .
nevertheless , dual conformal symmetry completely fixes the form of the four- and five - point amplitudes .
the form to which they are fixed turns out to be precisely the abdk / bds ansatz  @xcite , which was developed earlier based on the structure of infrared divergences of amplitudes , plus patterns observed at two and three loops . as discussed in , for six external legs one can build three different combinations of kinematic variables that are invariant under all dual conformal transformations , the cross - ratios @xmath40 , @xmath41 and @xmath42 .
at the six - point level , the abdk / bds ansatz fails at the first nontrivial order , two loops  @xcite , and dual conformal symmetry allows for an additional `` remainder function '' , @xmath43 , which depends on the three cross - ratios .    at strong coupling ,
large @xmath44 and large-@xmath12 , alday and maldacena used the ads / cft correspondence to map the problem of gluon scattering to that of strings moving in anti - de - sitter space  @xcite .
the ads background is weakly curved in this limit , relative to the string scale .
another way of saying this is that the two - dimensional string world - sheet stretches over a large area , relative to the scale on which strings fluctuate , allowing a semi - classical expansion to be used at large @xmath44 . in the leading term in the expansion
, different helicity configurations are not distinguished .
the amplitude is given simply by @xmath45 , where @xmath46 is the classical action that minimizes the area ( at least for scattering configurations with a euclidean interpretation ) .
using a @xmath47-duality transformation in string theory , alday and maldacena showed that the boundary conditions for the world - sheet are closed polygons , where each edge is a light - like segment given by the momentum @xmath48 of the @xmath49 gluon .
the solution to the minimal area problem was found explicitly for four - gluon scattering .
more recently , in order to handle more complicated kinematical configurations than the four - point case , the integrability of the string sigma - model action  @xcite has been exploited .
the minimal area problem has been solved by mapping it to a system of equations identical to those of the thermodynamical bethe ansatz  @xcite .
the first computation at strong coupling was for the four - gluon amplitude , and matched precisely the prediction of the abdk / bds ansatz  @xcite .
but it also suggested a weak - coupling correspondence between scattering amplitudes and wilson loops for closed light - like polygons , which was rapidly established for a variety of cases , initially for the simplest case of mhv scattering amplitudes  @xcite .
the wilson - loop correspondence has been very important both conceptually and technically .
conceptually , it is still not entirely clear why it happens at all at weak coupling .
technically , given its existence , it allows amplitudes to be computed in terms of wilson line integrals , which are generally somewhat more tractable than feynman loop integrals .
some recent conceptual advances have come from considering certain correlation functions with close - to - light - like separations  @xcite .
the relative simplicity of the wilson line integrals for the six - point ( hexagon ) case allowed the remainder function @xmath43 to be evaluated analytically in terms of goncharov polylogarithms  @xcite .
the expression in ref .
@xcite was then simplified considerably , to a few lines involving classical polylogarithms @xmath50 , using properties of the symbol operation which captures all of the essential analytic behavior of a multi - variable function  @xcite .
these results open the door to the possibility of finding simple analytic results for more general kinematic configurations .
for some configurations in which the external momenta are restricted to two space - time dimensions , two - loop analytic results are already available beyond six external legs  @xcite .
there is certainly the prospect of more multi - loop analytic results on the horizon , for both wilson loops and scattering amplitudes in planar @xmath0 super - yang - mills theory .
gravity seems to be a completely separate force from gauge theory .
it has a spin two force carrier instead of spin one , no color degrees of freedom , and a dimensionful coupling constant .
nevertheless , the two theories are intimately connected .
the ads / cft correspondence
@xcite relates the two , holographically , as a _ weak - strong _ duality .
as mentioned earlier , strong - coupling scattering of gluons has an alternate description in terms of strings moving in a weakly - curved five - dimensional gravitational background . on the other hand ,
there is also a _ weak - weak _ duality of some kind between gravity and gauge theory , directly in four dimensions : it is possible to write perturbative scattering amplitudes for gravitons in terms of `` double copies '' of gluon amplitudes .
the original examples of such relations are due to kawai , lewellen and tye ( klt )  @xcite , who found them by first deriving relations between closed and open string theory tree amplitudes , whose low - energy limits give tree - level relations in field - theory .
graviton amplitudes are represented in terms of sums of products of pairs of gluon amplitudes with different color orderings .
the klt relations hold not only for pure - graviton scattering , in terms of pure - gauge theory , but also for @xmath17 supergravity ( or any subsector of it ) in terms of amplitudes for @xmath0 sym .
more recently it has been recognized that other double - copy formulas exist  @xcite .
the new formulas are simpler in some sense than the klt relations ; they involve the direct squares of gauge - theory components , whereas the klt relations employ multiple permutations of full gauge - theory amplitudes .
as reviewed in this issue by carrasco and johansson  @xcite , these relations , combined with generalized unitarity , have important consequences for multi - loop amplitudes in @xmath17 supergravity , because they map the problem to a much simpler one , the corresponding multi - loop amplitudes for ( non - planar ) @xmath0 sym .    in the next two sections we discuss a few of the basic tools used to organize and illuminate the color and kinematic structure of amplitudes , which underlay the discussions in many of the other articles in this special issue .
in this section we describe some common conventions for organizing the color structure of su@xmath8 gauge theory amplitudes , which are used elsewhere in this special issue .
most of the articles implicitly discuss color - ordered partial amplitudes , which are particularly convenient in the large-@xmath12 limit in which planar diagrams dominate .
color - ordered amplitudes emerge from a `` trace - based '' color decomposition ( as also reviewed in _
e.g. _  refs .
@xcite ) . however , in some cases , in particular for describing subleading - color terms in amplitudes , and color - kinematic duality relations  @xcite , decompositions based on the su@xmath8 structure constants are more useful .    in general
, we consider two different su@xmath8 representations for the external states :    * the adjoint representation , for the gluon and any of its superpartners ( _ i.e. _  for all states in @xmath0 sym ) .
adjoint color indices are denoted by @xmath51 . * the fundamental ( defining ) representation @xmath12 , and its conjugate representation @xmath52 , for quarks and anti - quarks respectively .
fundamental color indices are denoted by @xmath53 , and anti - fundamental @xmath52 indices by @xmath54 .
the generators of su@xmath8 in the fundamental representation are traceless hermitian @xmath55 matrices , denoted by @xmath56 .
it is conventional when discussing helicity amplitudes to normalize the generators by @xmath57 , in order to avoid a proliferation of @xmath58 s .    in qcd ,
the group theory ( color ) factors in feynman diagrams are of two types : @xmath56 for the gluon - quark - antiquark vertex , and the su@xmath8 structure constants @xmath59 ( or products of @xmath59 s ) for all other vertices .
( it is also convenient to normalize the structure constants differently , using @xmath60 , where @xmath59 is the standard , textbook normalization . )
the color factors can be represented diagrammatically  @xcite using feynman - diagram notation , as in .
lines with arrows ( quarks ) carry fundamental representation indices ; curly lines ( gluons ) carry adjoint indices .
the adjoint representation is also a bi - fundamental ( @xmath61 ) representation from which the singlet , or trace component , has been projected out , as in the second equation for @xmath62 in the figure .
this relation also gives rise to the su@xmath8 fierz identity , ( t^a)_i_1^  _ 1 ( t^a)_i_2^  _ 2  =   _ i_1^  _ 2 _ i_2^  _ 1 - 1n_c _
i_1^  _ 1 _ i_2^  _ 2 , [ sunfierz ] which allows one to simplify products of traces of @xmath63 s .    4.8 truein      color - ordering  @xcite is related to the t  hooft double - line formalism  @xcite .
we begin by rewriting the color factors entirely in terms of the @xmath63 generators , using the relation ^abc if^abc= ( [ t^a , t^b ] t^c ) , [ removefabc ] which follows from the definition of the structure constants , and is depicted graphically in . after inserting repeatedly into the color factor for a typical feynman diagram , one obtains a large number of traces of the generic form @xmath64 . if the amplitude has external quark legs , then there will also be strings of @xmath63 s terminated by fundamental indices , of the form @xmath65 , one for each external quark - antiquark pair .
the number of traces can be reduced considerably by repeated use of the su@xmath8 fierz identity  ( [ sunfierz ] ) . in the case of tree - level amplitudes with external states in the adjoint representation , such as @xmath14-gluon amplitudes , sketches
how the color factors all may be reduced to a single trace , @xmath66 , for some permutation @xmath67 of the @xmath14 gluons .
this reduction leads to the _ trace - based color decomposition _ for @xmath14-gluon tree amplitudes , -20 mm a^_n(\{k_i , h_i , a_i } ) = g^n-2 -1.3 mm _ s_n / z_n -1.3 mm ( t^a_(1)t^a_(n ) )   a_n^((1^h_1),
,(n^h_n ) ) . [ treegluecolor ] here
@xmath68 is the full amplitude , with dependence on the external gluon momenta @xmath48 , @xmath69 , helicities @xmath70 , and adjoint indices @xmath71 . in qcd
, the gauge coupling @xmath5 is related to the strong coupling by @xmath72 .
the _ partial amplitudes _ ( or _ primitive _ or _ color - stripped amplitudes _ ) @xmath73 have had all the color factors removed , but contain all the kinematic information .
cyclic permutations of the arguments of a partial amplitude , denoted by @xmath74 , leave it invariant , because the associated trace is invariant under these operations .
however , all @xmath75 non - cyclic permutations , or orderings , of the partial amplitude appear in .
these permutations are denoted by @xmath76 .
5.0 truein    looking again at the way the different trace factors arise in , one sees that the partial amplitude @xmath77 only receives contributions from tree - level feynman diagrams that can be drawn on a plane , in which the cyclic ordering of the external legs , @xmath78 , matches the ordering of the arguments in @xmath79 .
therefore each partial amplitude can only have singularities in momentum invariants formed by squaring color - adjacent sums of momenta , such as @xmath80 , @xmath81 , _ etc .
_ , where all indices are defined modulo @xmath14 . in this way , the color decomposition  ( [ treegluecolor ] ) disentangles the kinematic complexity of the full amplitude @xmath82 .
similarly , tree amplitudes with two external quarks and @xmath83 gluons can be reduced to single strings of @xmath63 matrices , ^_n(q_1,g_3,
,g_n-1,|q_n ) & = & g^n-2 -1.3 mm _ s_n-2 -1.3 mm ( t^a_(2)t^a_(n-1))_i_1^  _
n + & & 1.3 cm a_n^(1_q,(2),  ,(n-1),n_|q ) , [ treequarkgluecolor ] where we have omitted the helicity labels , and numbers without subscripts in the argument of @xmath79 refer to gluons . in this case
there are @xmath84 terms , corresponding to all possible gluon orderings between the quarks .
color decompositions for amplitudes containing additional external quark pairs have also been described  @xcite .
we note that the partial amplitudes appearing in could also be used to describe the scattering of two fermions with different color quantum numbers , and @xmath83 gluons .
for example , if the fermions are gluinos in the adjoint representation , then the color decomposition has the form of , but the partial amplitudes , with two external fermions , are of the type given in .
this is easy to see if the two fermions are color - adjacent ; essentially all one needs to do to go between the two cases is remove one color line running between the two gluinos .
however , the other cyclic orderings can also be obtained from @xmath85 , using the kleiss - kuijf relations discussed in the next subsection .
this property illustrates the universality of the color - ordered primitive amplitudes alluded to in the introduction .    many articles in this issue consider loop - level amplitudes in which all of the external states are gluons in the adjoint representation .
the general color decomposition here is similar to , except that multiple color traces are now possible . roughly speaking , a loop of gluons carries a fundamental and an anti - fundamental index around the loop with it , and an external gluon can attach to either index .
the number of traces at @xmath39 loops is equal to @xmath86 .    at one loop ,
the full color decomposition for @xmath14 gluons is  @xcite -24 mm a^_n(\{k_i , h_i , a_i } ) & = & g^n .
[ loopgluecolor ] here @xmath87 are the partial amplitudes , @xmath74 and @xmath88 are the subsets of @xmath89 that leave the corresponding single and double trace structures invariant , and @xmath90 is the greatest integer less than or equal to @xmath91 .
the formula is more complicated than at tree level , but only by the need to keep track of more distinct trace structures and their various symmetries .    if there are @xmath92 flavors of quarks in the loop , it is easy to see that they only contribute to the single - trace coefficient @xmath93 in , and with a weight @xmath94 relative to the adjoint gluons .
( similar color decompositions are available for one - loop amplitudes with an external quark pair  @xcite . )
when one constructs the color - summed cross section , as is usually required for qcd applications , the contribution of the double - trace coefficients , @xmath87 for @xmath95 , is suppressed by a power of @xmath96 with respect to that of @xmath93 .
( it also turns out that the @xmath87 are not really independent , but can be computed as a sum over permutations of the @xmath93  @xcite . )
thus simplifies a lot in the large-@xmath12 limit .
the factor of @xmath97 simply comes from those terms in which all external gluons attach to the same fundamental index line , leaving the trace of the identity matrix for the untouched line .
correspondingly , if we only want the leading terms in the large-@xmath12 limit at @xmath39 loops , we can write the compact color decomposition , & & -1.5 cm a^_n(\{k_i , h_i , a_i})|_leadingcolor + & & -0.5 cm = g^n-2 ( g^2 n_c)^l _ s_n / z_n ( t^a_(1)t^a_(n ) )   a_n^(l)((1^h_1),
,(n^h_n ) ) , [ lloopleadingcolor ] where we have dropped the `` ; 1 '' index on the leading - color partial amplitude .
the t  hooft coupling @xmath10 emerges naturally in this limit .
the articles in this issue about multi - loop amplitudes in planar @xmath0  sym are generally concerned with the @xmath98 in . because each loop integration typically brings a factor of @xmath99 ( or @xmath100 when using dimensional regularization in @xmath22 )
, in some cases the precise normalizations of the @xmath98 will differ by this factor , as well as possibly other conventional factors .
there is another type of color decomposition that is actually more useful than the above trace - based decomposition for addressing certain issues , and particularly for working beyond the large-@xmath12 approximation .
we could call this approach the @xmath101-based decomposition , because it goes back to the original color factors built out of the structure constants @xmath59 .
( now we assume that all states are in the adjoint representation . )
the diagrammatic representation of a general @xmath101-based color factor at tree level has the structure of a cubic tree graph , with vertices given by @xmath102 s and propagators given by @xmath62 s .
however , these color factors are not all independent , due to the color jacobi identity .
this identity can be written in different ways .
one way , ^dae ^bce - ^dbe ^ace = ^abe ^cde , [ colorjacobiidentity ] which is also depicted in , corresponds to the fact that the structure constants are also the su@xmath8 generators in the adjoint representation , @xmath103 , so that their commutator gives @xmath104 contracted with @xmath105 .
4.0 truein    at tree level , one can use the identity  ( [ colorjacobiidentity ] ) repeatedly to transform all @xmath101-based color factors into a `` multi - peripheral '' form in which two selected gluons , say 1 and @xmath14 , are always at the end of a long chain of structure constants , and the other @xmath83 gluons are emitted from along the ladder  @xcite , -20 mm a_n^(\{a_i } ) & = & g^n-2 _
s_n-2 ^a_1 a_(2 ) x_1 ^x_1 a_(3 ) x_2 ^x_n-3
a_(n-1 ) a_n + & & a_n^(1,(2),
 ,(n-1),n ) + & = & g^n-2 _ s_n-2 ( f^a_(2 ) f^a_(n-1 ) ) _ a_1
a_n a_n^(1,(2),  ,(n-1),n ) . [ gluondecompnew ] notice that this representation is identical in form to the decomposition for two quarks and @xmath83 gluons in , except for the representation used for the su@xmath8 generator matrices .
although correspond to the cases of fundamental and adjoint representations , respectively , the same type of decomposition clearly holds for two external matter fields in an arbitrary su@xmath8 representation . the fact that the coefficients of the @xmath101-based decomposition  ( [ gluondecompnew ] ) are identical to the color - ordered partial amplitudes @xmath79 appearing in the trace - based decomposition can be demonstrated by contracting both representations with @xmath29 and keeping only the leading terms at large @xmath12  @xcite .
note that there are only @xmath84 terms in the @xmath101-based color decomposition , in contrast to the @xmath75 terms in the trace - based decomposition  ( [ treegluecolor ] ) .
therefore there are identities between the @xmath79 , which are group - theoretical in nature  generalizations of certain @xmath106 decoupling identities  @xcite .
these identities can be used to put any two legs next to each other , for example legs 1 and @xmath14 : a_n^(1,\{},n,\ { } ) = ( -1)^n _ _ \{}\{^t } a_n^(1,(\{}\{^t}),n ) , [ kkequation ] where @xmath107 and @xmath108 are two sets of external gluons , whose union is @xmath109 .
also , @xmath110 is the number of elements in the set @xmath108 , the set @xmath111 is @xmath108 with the ordering reversed , and @xmath112 is the set of ordered permutations , or `` mergings '' , of the two sets @xmath107 and @xmath111 that preserve the ordering of the @xmath113 within @xmath107 and of the @xmath114 within @xmath111 , while allowing for all possible relative orderings of the @xmath113 with respect to the @xmath114 .
these relations are known as the kleiss - kuijf relations  @xcite .
they can be derived from by expanding out the @xmath102 factors using , identifying the coefficients of the trace structures , and comparing them with the original trace - based decomposition  ( [ treegluecolor ] )  @xcite .
in fact , it has been realized recently that there are not even @xmath84 independent color - ordered amplitudes @xmath79 , but only @xmath115  @xcite , as reviewed in this issue by carrasco and johansson  @xcite .
the additional linear relations between @xmath79 s are not purely group - theoretical in nature , but also involve kinematical factors .
they allow one to put any _ three _ external gluon legs next to each other .
for example , for legs 1 , 2 and @xmath14 , the identities take the form , -1.0 cm a_n^((1),(2),  ,(n ) ) = _ s_n-3 k^ ( ) _ a_n^(1,2,(3),  ,(n-1),n ) , [ bcj3 ] where @xmath116 and @xmath117 are permutations and @xmath118 are kinematic - dependent ( but not state - dependent ) coefficients  @xcite .
these relations were first proved using open string theory  @xcite , using a contour deformation argument  @xcite similar to that used in deriving the klt relations  @xcite between gravity and gauge theory amplitudes .
they have also been proven directly in field theory , using a recursive argument  @xcite .
although expresses the consequences of certain color - kinematic identities in the trace basis , the identities themselves are more naturally stated in an @xmath101-based approach . for this purpose
it is useful to _ not _ transform the @xmath101-based color factors into multi - peripheral form , but leave the color decomposition free , in terms of the set of all cubic graphs @xmath119 , writing _
n^= g^n-2 _ i_3 .
[ bcjtreerep ] the color factors @xmath120 are products of structure constants for each vertex @xmath121 of the @xmath49 graph , contracted together along the propagators , @xmath122 .
the factors of @xmath123 are scalar propagators , one for each internal line in the graph , while the @xmath124 are kinematical numerator factors . consider a triplet of graphs in @xmath119 , which are identical except for a region from which four lines emanate . within this region
they differ according to the three ways of joining four legs with two cubic vertices , shown in .
call the three graphs @xmath125 , @xmath126 and @xmath127 , and let the color factors @xmath128 be normalized ( signed ) so that they obey the jacobi identity  ( [ colorjacobiidentity ] ) as @xmath129
. then the statement of color - kinematic duality  @xcite is that the associated numerator factors should be related by @xmath130 , for every such triplet of graphs .
a simple and intuitive argument for these relations has been provided based on the heterotic string  @xcite .
the existence of such a representation with _
@xmath124 ( _ i.e. _  polynomials in the momenta ) has been checked in various examples , not only at tree level , but also through three loops for @xmath0 sym amplitudes  @xcite . as reviewed here by carrasco and johansson  @xcite ,
the existence of a representation for gauge theory amplitudes satisfying @xmath130 for each triplet of jacobi - related cubic graphs has important implications .
large classes of numerator factors are related to each other , and gravitational amplitudes are easily constructed from gauge theory ones .
now that we have described two different ways to organize the color quantum numbers for scattering amplitudes , let us turn to kinematical issues , including convenient choices of external states and kinematical variables .      in qcd ,
the helicities of massless quarks are conserved by their chirality - preserving interactions with gluons .
hence it is natural to use a helicity basis for the quark amplitudes . in four - component notation ,
the spinors for the external states are taken to be @xmath131 for a quark with momentum @xmath132 , and @xmath133 for an anti - quark . however , in the massless limit these spinors can be chosen to be equal to each other , @xmath134 , and we can use two - component weyl spinors as well .
we will see that helicity amplitudes for external gluons can be built from the same objects .
we want to consider amplitudes with @xmath14 different momenta @xmath48 , @xmath69 , and to do so in a way that respects crossing symmetry .
thus we take all the momenta to be outgoing , so that momentum conservation reads @xmath135 . in a realistic process with some incoming particles ,
the physical momentum for each of the incoming particles is simply the _ negative _ of the corresponding @xmath48 . also , the physical helicity @xmath136 is the negative of the one by which we label the amplitude .
( the spin @xmath137 does not reverse under crossing , but @xmath138 does . )
the two - component spinor for the @xmath49 state , if it is an outgoing fermion with helicity @xmath139 , has several notations  @xcite : ( _ i ) _  & &  |i^+  |k_i^+  [ u_+(k_i ) ] _
=  [ v_-(k_i ) ] _ , + ( _ i ) _  & &  |i^-  |k_i^-  [ u_-(k_i ) ] _  =  [ v_+(k_i ) ] _ . similarly ,
the conjugate spinors are ( _ i)^  & &  i^+|  k_i^+|  [ ] ^=  [ ] ^ , + ( _ i)^  & &  i^-|  k_i^-|  [ ] ^  =  [ ] ^. two - component indices are raised and lowered with the two - dimensional antisymmetric tensors @xmath140 , @xmath141 , _ etc . _
the spinor inner products come from contracting spinors for different momenta with these tensors , .l  & = &   ^ ( _ j)_(_l ) _  =  j^- | l^+  =   u_+(k_l ) , [ spadef ] + .l  & = &   ^ ( _ j ) _
_  =  j^+ | l^-  =   u_-(k_l ) .
[ spbdef ] spinor products are antisymmetric under exchange of labels , @xmath142 , @xmath143 , and @xmath144 .
lorentz vectors can be written as bi - spinors , or @xmath145 matrices , by contracting them with the pauli matrices .
a massless vector written in this way factorizes into the product of two massless spinors , as ( _ i ) _
k_i^ ( _ ) _  =  u_+(k_i )  =  ( _ i)_(_i ) _ .
[ masslesskfact ] this factorization plays a role in the bcfw complex - momentum shift , which is best described in terms of the @xmath44 and @xmath146 variables  @xcite , as is done elsewhere in this issue  @xcite .
amplitudes with external gluons can also be described in terms of the @xmath147 and @xmath148 variables , thanks to the spinor - helicity formalism  @xcite .
the polarization vector @xmath149 for an outgoing massless vector particle with momentum @xmath132 and helicity @xmath150 , is required to be transverse to @xmath132 , @xmath151 . in the spinor - helicity formalism
, @xmath152 is also chosen to be transverse to another massless momentum @xmath153 , called the _ reference momentum _
( which should not be parallel to @xmath132 but is otherwise arbitrary ) .
note that @xmath132 and @xmath153 span a two - dimensional subspace of four - dimensional space - time . because @xmath154 , the two orthogonal directions to this subspace are spanned by @xmath155 and @xmath156 .
the spinor - helicity polarization vectors live in this subspace ; indeed , they are proportional to these two vectors , and are normalized by the condition that @xmath157 : _
^(k , q ) = . [ spinorheldef ] as a bi - spinor , @xmath158 is given by -1 cm [ ^+(k_i , q_i ) ] _
= , _ = - . [ spinorheldef2 ] if each @xmath159 is chosen to be another momentum in the process , say @xmath160 , then it is clear from the general form of the feynman rules ( using also ) that the full amplitude can be built entirely out of the spinor products @xmath161 and @xmath162 , for @xmath163 .    what are the basic properties of the spinor products which make them so convenient variables for helicity amplitudes ?
first of all , for real momenta , the two types of spinor products are complex conjugates of each other , @xmath164 . also , for either real or complex momenta they `` square '' to give the ordinary dot products of momenta ,
.j.l  =  =  2 k_jk_l = s_jl .
[ sfromspinors ] ( is written in the convention usually used in the qcd literature , but the reader should be aware that another convention exists , in which the sign of @xmath162 is reversed so that @xmath165 . )    for real momenta , the conjugation property and imply that .l  =   e^i_jl , .l  =   e^-i_jl , [ realmomsqrt ] where @xmath166 is a phase
this phase shifts as @xmath160 and @xmath167 are rotated in azimuthal angle around the axis corresponding to their sum , @xmath168 .
one way of seeing the utility of the spinor products for helicity amplitudes is to examine the limits in which two particles become collinear , @xmath169 .
( a ) shows that the typical behavior of amplitudes in a massless scalar field theory is @xmath170 , coming just from the scalar propagator .
however , in massless gauge theory , there are numerator factors in the feynman diagrams which lessen this divergence .
more physically , the factorization in the collinear limit should be onto a physical gluon state with helicity @xmath171 .
however , the sum of the two external helicities must be @xmath172 or 0 .
therefore there is at least a @xmath171 mismatch of the spin angular momentum along the @xmath173 axis , as illustrated in ( b ) .
the spin angular - momentum mismatch requires some orbital angular momentum , which in turn causes a suppression of the magnitude of the amplitude , from @xmath174 to @xmath175 .
it also dictates a phase shift under azimuthal rotation of @xmath160 and @xmath167 about the @xmath173 axis .
both of these properties are captured by the spinor products , making them ideal for describing helicity amplitudes .
4.5 truein    the simplest nonvanishing tree - level gluon amplitudes are the parke - taylor or mhv amplitudes  @xcite , a_n^(1^+,  ,l^-,  ,m^-,  ,n^+ ) = i .m^4 1.22.3.1 .
[ mhvpt ] here exactly two gluons , @xmath176 and @xmath177 , have negative helicity , and the remaining @xmath83 gluons have positive helicity .
the amplitudes with zero or one negative gluon helicity vanish by supersymmetry ward identities  @xcite .
the only denominator factors in these amplitudes are the spinor products for color - adjacent pairs of momenta , @xmath178 .
they capture such universal collinear limits  @xcite as , & & -1 cm a_n^(1^+,  ,l^-,  ,m^-,  ,(n-1)^+,n^+)|_n-1n [ collimitexample ] + & ~ & a_n-1^(1^+,
,l^-,  ,m^-,
,p^+ ) , where legs @xmath179 and @xmath14 are becoming parallel , with @xmath180 , @xmath181 and @xmath182 .      another way to see why the spinor products are so useful is to consider the special case @xmath183 for the mhv amplitudes , a_3^(1 ^ -,2 ^ -,3^+ ) = i .2 ^ 4 1.22.33.1 .
[ mmp ] one might think that the scattering of three massless particles is singular , because all of the @xmath184 vanish since they are equal to the squared momentum of the remaining particle . indeed , for real momenta
, implies that all the spinor products vanish too , and there is no way to build a nonvanishing and nonsingular amplitude .
however , for complex momenta , the complex - conjugation relation @xmath164 no longer holds , and it is possible for half the spinor products to be nonvanishing , while the other half vanish .    specifically , we can choose the three negative - helicity two - component spinors to be proportional  @xcite , _
1^  _ 2^  _ 3 ^ .
[ spbvanish ] then according to we have @xmath185 . however , the other three spinor products , @xmath186 , @xmath187 and @xmath188 , are allowed to be nonzero .
( this is consistent with the vanishing of , which still holds for complex momenta . ) for this choice of complex kinematics , the mhv three - point amplitude  ( [ mmp ] ) is nonvanishing and nonsingular , while the conjugate `` @xmath189 '' amplitude , a_3^(1^+,2^+,3 ^ - ) = -i 1.2 ^ 41.22.33.1 , [ ppm ] vanishes .
in contrast , for the conjugate kinematics with _ 1^  _
2^  _ 3 ^ , [ spavanish ] vanishes while is nonvanishing and nonsingular .    the three - point amplitudes  ( [ mmp ] ) and ( [ ppm ] ) , as well as related amplitudes with two matter particles and one gluon , are very important in massless gauge theory .
the bcfw relations build all tree - level scattering amplitudes recursively from them , using factorization onto multi - particle poles , as well as onto complex - momentum versions of the collinear poles shown in .
of course feynman diagrams can also be used to build amplitudes from three- and four - point vertices .
however , feynman vertices are typically evaluated with off - shell lines emanating from them , which makes them gauge - dependent .
in contrast , are on shell , albeit for complex momenta , and therefore they are fully gauge invariant . in the recursive construction of @xmath14-gluon amplitudes ,
the four - gluon vertex is never needed , essentially because it is related to the three - vertex by gauge transformations .
another reason why the three - point amplitudes are important is for building up loop amplitudes through generalized unitarity , as reviewed in this issue by ita  @xcite and britto  @xcite at one loop , and by bern and huang  @xcite and carrasco and johansson  @xcite at the multi - loop level .
typically , generalized unitarity cuts place many propagators on shell , so they very often contain at least one three - point tree amplitude
. finally , the relations between gravity and gauge theory tree amplitudes are the simplest of all for the three - point case , where they can be written as an exact square , m_3^(1,2,3 ) = [ a_3^(1,2,3 ) ] ^2 , [ klt3 ] ignoring overall coupling - constant factors .      for planar ( large @xmath12 ) gauge theory amplitudes  which includes all tree - level amplitudes  new sets of kinematic variables have proved very useful for identifying new structures and symmetries .    because the external legs have a definite cyclic order in the planar case , they carve the plane into sectors . to each such sector
, one can assign a vertex @xmath190 for a dual graph , as shown in ( a ) .
lines connecting vertices of the dual graph cross lines of the original graph .
the latter lines carry momenta .
differences between @xmath191 s are computed according to the net momentum flowing through the lines of the original graph .
thus in ( a ) , adjacent sectors are separated by @xmath192 , or more generally , x_i , j^x_i^- x_j^= k_j^+ k_j+1^+ + k_i-1^. [ sectorvariable ] the @xmath191 are often called _ dual coordinates _ ( or sometimes region or sector variables ) , but they are essentially momenta , not coordinates . in a planar loop graph
we can assign additional @xmath191 variables to the faces inside each loop , corresponding to the independent loop momenta .
( b ) shows a two - loop example . here
the thick lines of the dual graph that are shown each cross ( and therefore correspond to ) one propagator of the planar double box integral .
4.5 truein    more precisely , differences of the @xmath191 are momenta .
the @xmath191 themselves are not constrained by momentum conservation .
they are constrained by on - shell conditions , @xmath193 . however , unlike momenta , it is possible to invert the @xmath191 s , according to x_i^ , x_i , j^2 .
[ dualinversion ] remarkably , this inversion is a symmetry of both integrands and amplitudes in planar @xmath0 sym .
the integrands are functions of the combinations @xmath194 , which are manifestly invariant under lorentz transformations and translations of the @xmath191 s .
( translation invariance is guaranteed simply because the amplitude depends on momenta , which are differences of @xmath191 s . )
therefore the inversion  ( [ dualinversion ] ) combines with poincar invariance to generate the conformal group so(4,2 ) .
this group acts on the dual variables @xmath191 , not on space - time coordinates ; hence the symmetry is referred to as dual conformal invariance .    at the loop level
, dual conformal invariance can be spoiled by the infrared regulator needed for on - shell amplitudes . a loop - integration measure in four dimensions , @xmath195 ,
transforms under an inversion  ( [ dualinversion ] ) as , @xmath196 . in @xmath197
this factor precisely cancels terms from inverting the associated propagators in graphs such as the one in ( b ) .
however , in dimensional regularization , we have @xmath198 , and the invariance is lost due to extra factors in the transformation of the loop integration measure . as discussed in detail in this issue by henn
@xcite , for planar @xmath0 sym it is possible to retain dual conformal invariance at the loop level by using instead a higgs regulator  @xcite , which generates particle masses that can also be thought of as extra components of the @xmath190 .
although amplitude computations with the higgs regulator have not been pushed quite as far yet as in dimensional regularization , there are certain advantages to this approach .
for example , in dimensional regularization , terms in lower - loop amplitudes that vanish as @xmath199 often have to be computed , because they can multiply pole terms , @xmath200 , in other amplitudes ( or in the squares of amplitudes needed to construct differential cross sections ) .
the article by schabinger in this issue discusses such @xmath201 contributions  @xcite .
in contrast , with the higgs regulator , the singularities have the form of powers of logarithms , @xmath202 , and power - suppressed terms of the form @xmath203 can always be dropped , because they vanish faster than any power of a logarithm .
how constraining is dual conformal symmetry ?
in general it can only determine amplitudes up to arbitrary functions of the cross - ratios , u_ijkl   , [ genu ] because these variables are invariant under the inversion  ( [ dualinversion ] ) .
however , the on - shell constraints @xmath204 imply that there are no non - trivial cross ratios for four- and five - particle scattering . at the six - point level
there are three non - trivial cross ratios , utilizing the allowed @xmath205 and @xmath206 kinematic invariants : [ defu ] -1.8 cm u_1 = = , u_2 = = , u_3 = = . as mentioned earlier , the remainder function @xmath207 characterizing the two - loop mhv six - point amplitude in planar @xmath0 sym is a function of these three variables , and this function is now known analytically in terms of polylogarithms  @xcite .    in planar
@xmath0 sym , dual conformal invariance combines with dual @xmath0 supersymmetry to generate dual superconformal symmetry  @xcite .
it is natural to consider super - amplitudes , which package sets of amplitudes together by using an @xmath0 on - shell superfield  @xcite , [ super - wave ] -2.3 cm ( ) = g^+ + ^a g_a + ^a ^b _ ab + ^a^b^c _ abcd |g^d + ^a^b^c ^d _ abcd g^- . here
@xmath208 are the @xmath209 helicity gluons , @xmath210 and @xmath211 the four flavors of @xmath212 helicity gluinos , and @xmath213 the six real @xmath214 helicity scalar states
. integrations over the grassmann variables @xmath21 , @xmath215 , can be used to pick off the desired component amplitudes from the super - amplitudes , defined by [ super - amplitude ] a_n(_i )  a ( _ 1(_1 ) ,  , _
just as the @xmath191 bosonic variables automatically satisfy momentum conservation , @xmath216 , one can construct their superpartners @xmath217 , whose differences satisfy _ i^a - _
i+1^a  =  _ i^_i^a , [ dualsuper ] for @xmath69 .
super - momentum conservation constitutes eight constraints : q^a  _ i=1^n _ i^_i^a  =  0 . [ supercons ] like ordinary momentum conservation , it is satisfied by virtue of the periodicity of the dual coordinates , @xmath218 , @xmath219 .
the mhv super - amplitude generalizing is _
n^,mhv = i = i . [ mhvsuper ] it transforms covariantly under dual superconformal transformations  @xcite which extend so(4,2 ) to the superalgebra psu(2,2@xmath2204 ) .
furthermore , the solutions to the bcfw super - recursion relation for amplitudes with more negative helicities ( nmhv , nnmhv , _ etc .
_ ) are given by the product of this super - amplitude with collections of dual superconformal invariants , generically denoted by @xmath221  @xcite .
superconformal and dual superconformal symmetry are separate symmetries , which partly overlap .
they close into a very large symmetry group called the yangian  @xcite , as discussed in this issue by bargheer , beisert and loebbert  @xcite and by drummond  @xcite .
the generic generator of the yangian is not local in the sense of being a sum of single - particle terms ( like @xmath132 and @xmath153 ) ; rather it is a sum of multi - particle operators .
( for the dual superconformal generators , the nonlocality is mild ; only the sum of two - particle operators appears . )
there are certain `` anomalies '' in the action of the yangian on amplitudes . at tree level
, the anomalies only act at the boundary of the @xmath14-particle phase space , where two particles become collinear , as in  @xcite . at loop - level
they act in the full phase - space , due to singularities in the loop integration  @xcite . however , it is possible to deform either the symmetry algebra or the representation in a suitable way so as to preserve the invariance  @xcite
. it will be very interesting to see whether the constraints from the yangian are sufficient to determine multi - loop amplitudes directly , without having to pass through the computation of loop integrals .
another very recent development , in the same vein of using symmetries or dynamical principles to bypass standard loop ( or wilson line ) integrals , concerns polygonal wilson loops .
one considers an operator - product - like expansion for the wilson loop , associated with a limit in which multiple momenta become collinear  @xcite .
this expansion can be carried out at both weak and strong coupling . at weak coupling
it has been used to compute the discontinuity of the integrated two - loop amplitude , directly in terms of one - loop quantities  @xcite .
there has been a remarkable array of breakthroughs in our understanding of scattering amplitudes over the past few years , as described in the review articles in this special issue of journal of physics a. this article has represented an overview of many of these developments , as well as providing some basic material as an introduction to the remaining contributions . however , it certainly did not do justice to all of the directions being pursued currently .
one of the other remarkable features of this field has been the interplay and cross - fertilization between formal developments and phenomenology .
the parke - taylor tree amplitudes were found in the course of trying to understand patterns arising in the structure of five- and six - gluon tree - level amplitudes . a similar analysis of the structure of the five - gluon amplitude at one loop in qcd led to the development of the unitarity method , which has been indispensable for computing high - loop - order amplitudes in maximally supersymmetric gauge theory and gravity , but which has also produced new nlo qcd results for colliders  @xcite .
a `` radiation zero '' found over thirty years ago in the electroweak process @xmath222
@xcite was studied for a general gauge theory soon thereafter  @xcite ; the relations found there were recognized much later as the four - point versions of a more general color - kinematics duality  @xcite .
wilson loops , used to characterize the effects of soft gluons in qcd , also turned out to be in exact correspondence to full mhv amplitudes in planar @xmath0 super - yang - mills theory .
this special issue should be of interest , not only to experienced practitioners in the field , but also to newcomers who want to get started . as the history of the field shows , new developments can and will come from totally unanticipated angles .
many of the new ideas in the future may well come to physicists whose interest in the remarkable aspects of scattering amplitudes was sparked , at least in part , by reading the articles assembled in this special issue .
* acknowledgments *    i am grateful to radu roiban , mark spradlin and anastasia volovich for inviting me to contribute this introductory article , and i thank radu roiban also for valuable comments on the manuscript .
this work was supported by the us department of energy under contract de  ac0276sf00515 .
the figures were generated using jaxodraw  @xcite , based on axodraw  @xcite .
10      j.  m.  maldacena , adv .
phys .  * 2 * , 231 ( 1998 ) [ int .
j.  theor .
phys .  * 38 * , 1113 ( 1999 ) ] [ hep - th/9711200 ] ; + s.  s.  gubser , i.  r.  klebanov and a.  m.  polyakov , phys .
b * 428 * , 105 ( 1998 ) [ hep - th/9802109 ] ; + o.  aharony , s.  s.  gubser , j.  m.  maldacena , h.  ooguri and y.  oz , phys .  rept .   *
323 * , 183 ( 2000 ) [ hep - th/9905111 ] .
l.  f.  alday and j.  m.  maldacena , jhep * 0706 * , 064 ( 2007 ) [ 0705.0303 [ hep - th ] ] . s.  j.  parke and t.  r.  taylor , phys .
lett .   * 56 * , 2459 ( 1986 ) .
v.  p.  nair , phys .
b * 214 * , 215 ( 1988 ) .
m.  l.  mangano and s.  j.  parke , phys .
* 200 * , 301 ( 1991 ) [ hep - th/0509223 ] .
m.  t.  grisaru , h.  n.  pendleton and p.  van nieuwenhuizen , phys .
d * 15 * , 996 ( 1977 ) .
m.  t.  grisaru and h.  n.  pendleton , nucl .
b * 124 * , 81 ( 1977 ) .
s.  j.  parke and t.  r.  taylor , phys .
b * 157 * , 81 ( 1985 ) [ erratum - ibid .   * 174b * , 465 ( 1986 ) ] ; + z.  kunszt , nucl .  phys .
b * 271 * , 333 ( 1986 ) .
h.  elvang , d.  z.  freedman and m.  kiermaier , jhep * 1010 * , 103 ( 2010 ) [ 0911.3169 [ hep - th ] ] .
h.  elvang , d.  z.  freedman and m.  kiermaier , to appear in j.  phys .
a [ 1012.3401 [ hep - th ] ] .
e.  witten , commun .
phys .   * 252 * , 189 ( 2004 ) [ hep - th/0312171 ] .
r.  penrose , j.  math .
* 8 * , 345 ( 1967 ) .
f.  cachazo , p.  svrek and e.  witten , jhep 0409:006 ( 2004 ) [ hep - th/0403047 ] .
a.  brandhuber , b.  spence and g.  travaglini , to appear in j.  phys .
a [ 1103.3477 [ hep - th ] ] .
t.  adamo , m.  bullimore , l.  mason and d.  skinner , to appear in j.  phys .
a [ 1104.2890 [ hep - th ] ] .
r.  britto , f.  cachazo , b.  feng and e.  witten , phys .
lett .   * 94 * , 181602 ( 2005 ) [ hep - th/0501052 ] .
r.  britto , f.  cachazo and b.  feng , nucl .
b * 715 * , 499 ( 2005 ) [ hep - th/0412308 ] .
n.  arkani - hamed , f.  cachazo and j.  kaplan , jhep * 1009 * , 016 ( 2010 ) [ 0808.1446 [ hep - th ] ] .
a.  brandhuber , p.  heslop and g.  travaglini , phys .
d * 78 * , 125005 ( 2008 ) [ 0807.4097 [ hep - th ] ] .
j.  m.  drummond and j.  m.  henn , jhep * 0904 * , 018 ( 2009 ) [ 0808.2475 [ hep - th ] ]
. j.  drummond , to appear in j.  phys .
a.    z.  bern , l.  j.  dixon , d.  c.  dunbar and d.  a.  kosower , nucl .
b * 425 * , 217 ( 1994 ) [ hep - ph/9403226 ] .
z.  bern , l.  j.  dixon , d.  c.  dunbar and d.  a.  kosower , nucl .
b * 435 * , 59 ( 1995 ) [ hep - ph/9409265 ] .
r.  britto , to appear in j.  phys .
a [ 1012.4493 [ hep - th ] ] .
z.  bern , l.  j.  dixon and d.  a.  kosower , nucl .
b * 513 * , 3 ( 1998 ) [ hep - ph/9708239 ] .
r.  britto , f.  cachazo and b.  feng , nucl .
b * 725 * , 275 ( 2005 ) [ hep - th/0412103 ] .
n.  e.  j.  bjerrum - bohr and p.  vanhove , jhep * 0804 * , 065 ( 2008 ) [ 0802.0868 [ hep - th ] ] ; jhep * 0810 * , 006 ( 2008 ) [ 0805.3682 [ hep - th ] ] ; fortsch .
phys .   * 56 * , 824 ( 2008 ) [ 0806.1726 [ hep - th ] ] .
z.  bern , l.  j.  dixon and d.  a.  kosower , jhep * 0001 * , 027 ( 2000 ) [ hep - ph/0001001 ] ; + z.  bern , a.  de freitas and l.  j.  dixon , jhep * 0203 * , 018 ( 2002 ) [ hep - ph/0201161 ] .
z.  bern and y .- t .
huang , to appear in j.  phys .
a [ 1103.1869 [ hep - th ] ] .
g.  ossola , c.  g.  papadopoulos and r.  pittau , nucl .
b * 763 * , 147 ( 2007 ) [ hep - ph/0609007 ] ; jhep * 0803 * , 042 ( 2008 ) [ 0711.3596 [ hep - ph ] ] .
r.  k.  ellis , w.  t.  giele and z.  kunszt , jhep * 0803 * , 003 ( 2008 ) [ 0708.2398 [ hep - ph ] ] . c.  f.  berger _
_ , phys .
d * 78 * , 036003 ( 2008 ) [ 0803.4180 [ hep - ph ] ] . w.  siegel , phys .  lett .
b * 84 * , 193 ( 1979 ) ; + d.  m.  capper , d.  r.  t.  jones and p.  van nieuwenhuizen , nucl .
b * 167 * , 479 ( 1980 ) ; + l.  v.  avdeev and a.  a.  vladimirov , nucl .
b * 219 * , 262 ( 1983 ) .
z.  bern and d.  a.  kosower , nucl .
b * 379 * , 451 ( 1992 ) ; + z.  bern , a.  de freitas , l.  dixon and h.  l.  wong , phys .  rev .
d * 66 * , 085002 ( 2002 ) [ hep - ph/0202271 ] .
r.  akhoury , phys .
d * 19 * , 1250 ( 1979 ) ; + a.  h.  mueller , phys .
d * 20 * , 2037 ( 1979 ) ; + j.  c.  collins , phys .
d * 22 * , 1478 ( 1980 ) ; in _ perturbative qcd _ , ed .
a.  h. mueller , advanced series on directions in high energy physics , vol .
5 ( world scientific , singapore , 1989 ) [ hep - ph/0312336 ] ; + a.  sen , phys .  rev .
d * 24 * , 3281 ( 1981 ) ; phys .
d * 28 * , 860 ( 1983 ) .
l.  magnea and g.  sterman , phys .
d * 42 * , 4222 ( 1990 ) .
s.  catani , phys .
b * 427 * , 161 ( 1998 ) [ hep - ph/9802439 ] .
g.  sterman and m.  e.  tejeda - yeomans , phys .
b * 552 * , 48 ( 2003 ) [ hep - ph/0210130 ] .
g.  p.  korchemsky , mod .
lett .  a * 4 * , 1257 ( 1989 ) ; + g.  p.  korchemsky and g.  marchesini , nucl .
b * 406 * , 225 ( 1993 ) [ hep - ph/9210281 ] .
n.  beisert , b.  eden and m.  staudacher , j.  stat .
mech .   * 0701 * , p01021 ( 2007 ) [ hep - th/0610251 ] .
z.  bern , l.  j.  dixon and v.  a.  smirnov , phys .
d * 72 * , 085001 ( 2005 ) [ hep - th/0505205 ] .
j.  bartels , l.  n.  lipatov and a.  sabio  vera , phys .
d * 80 * , 045002 ( 2009 ) [ 0802.2065 [ hep - th ] ] ; eur .
j.   c * 65 * , 587 ( 2010 ) [ 0807.0894 [ hep - th ] ] .
j.  bartels , l.  n.  lipatov and a.  prygarin , to appear in j.  phys .  a [ 1104.0816 [ hep - th ] ]
. l.  f.  alday , j.  m.  henn , j.  plefka and t.  schuster , jhep * 1001 * , 077 ( 2010 ) [ 0908.0684 [ hep - th ] ] . j.  m.  henn , to appear in j.  phys .
a [ 1103.1016 [ hep - th ] ] .
r.  k.  ellis , w.  t.  giele , z.  kunszt and k.  melnikov , nucl .
b * 822 * , 270 ( 2009 ) [ 0806.3467 [ hep - ph ] ] .
s.  lal and s.  raju , phys .
d * 81 * , 105002 ( 2010 ) [ 0910.0930 [ hep - th ] ] .
a.  v.  kotikov , l.  n.  lipatov , a.  i.  onishchenko and v.  n.  velizhanin , phys .
b * 595 * , 521 ( 2004 ) [ erratum - ibid .
b * 632 * , 754 ( 2006 ) ] [ hep - th/0404092 ] .
z.  bern , m.  czakon , l.  j.  dixon , d.  a.  kosower and v.  a.  smirnov , phys .
d * 75 * , 085010 ( 2007 ) [ hep - th/0610248 ] .
s.  g.  naculich , h.  nastase and h.  j.  schnitzer , jhep * 0811 * , 018 ( 2008 ) [ 0809.0376 [ hep - th ] ] . n.  arkani - hamed , j.  l.  bourjaily , f.  cachazo , s.  caron - huot and j.  trnka , jhep * 1101 * , 041 ( 2011 ) [ 1008.2958 [ hep - th ] ] .
r.  h.  boels , jhep * 1011 * , 113 ( 2010 ) [ 1008.3101 [ hep - th ] ] . j.  m.  drummond , j.  henn , g.  p.  korchemsky and e.  sokatchev , nucl .
b * 826 * , 337 ( 2010 ) [ 0712.1223 [ hep - th ] ] . c.  anastasiou , z.  bern , l.  j.  dixon and d.  a.  kosower , phys .
lett .   * 91 * , 251602 ( 2003 ) [ hep - th/0309040 ] .
j.  m.  drummond , j.  henn , g.  p.  korchemsky and e.  sokatchev , nucl .
b * 815 * , 142 ( 2009 ) [ 0803.1466 [ hep - th ] ] ; + z.  bern , l.  j.  dixon , d.  a.  kosower , r.  roiban , m.  spradlin , c.  vergu and a.  volovich , phys .
d * 78 * , 045007 ( 2008 ) [ 0803.1465 [ hep - th ] ] .
i.  bena , j.  polchinski and r.  roiban , phys .
d * 69 * , 046002 ( 2004 ) [ hep - th/0305116 ] .
l.  f.  alday , d.  gaiotto and j.  maldacena , 0911.4708 [ hep - th ] ; + l.  f.  alday , j.  maldacena , a.  sever and p.  vieira , j.  phys .
a * 43 * , 485401 ( 2010 ) [ 1002.2459 [ hep - th ] ] . j.  m.  drummond , g.  p.  korchemsky and e.  sokatchev , nucl .
b * 795 * , 385 ( 2008 ) [ 0707.0243 [ hep - th ] ] ; + a.  brandhuber , p.  heslop and g.  travaglini , nucl .
b * 794 * , 231 ( 2008 ) [ 0707.1153 [ hep - th ] ] ; + j.  m.  drummond , j.  henn , g.  p.  korchemsky and e.  sokatchev , nucl .
b * 795 * , 52 ( 2008 ) [ 0709.2368 [ hep - th ] ] .
l.  f.  alday , b.  eden , g.  p.  korchemsky , j.  maldacena and e.  sokatchev , 1007.3243 [ hep - th ] ; + b.  eden , g.  p.  korchemsky and e.  sokatchev , 1007.3246 [ hep - th ] ; + l.  j.  mason and d.  skinner , jhep * 1012 * , 018 ( 2010 ) [ 1009.2225 [ hep - th ] ] ; + s.  caron - huot , 1010.1167 [ hep - th ] ; + a.  v.  belitsky , g.  p.  korchemsky and e.  sokatchev , 1103.3008 [ hep - th ] . v.  del duca , c.  duhr and v.  a.  smirnov , jhep * 1003 * , 099 ( 2010 ) [ 0911.5332 [ hep - ph ] ] ; jhep * 1005 * , 084 ( 2010 ) [ 1003.1702 [ hep - th ] ] .
a.  b.  goncharov , m.  spradlin , c.  vergu and a.  volovich , phys .
lett .   * 105 * , 151605 ( 2010 ) [ 1006.5703 [ hep - th ] ] . v.  del duca , c.  duhr and v.  a.  smirnov , jhep * 1009 * , 015 ( 2010 ) [ 1006.4127 [ hep - th ] ] ; + p.  heslop and v.  v.  khoze , jhep * 1011 * , 035 ( 2010 ) [ 1007.1805 [ hep - th ] ]
. h.  kawai , d.  c.  lewellen and s .- h .
h.  tye , nucl .
b * 269 * , 1 ( 1986 ) .
z.  bern , j.  j.  m.  carrasco and h.  johansson , phys .
d * 78 * , 085011 ( 2008 ) [ 0805.3993 [ hep - ph ] ] .
z.  bern , j.  j.  m.  carrasco and h.  johansson , phys .
* 105 * , 061602 ( 2010 ) [ 1004.0476 [ hep - th ] ] . l.  j.  dixon , in _ qcd & beyond : proceedings of tasi 95 _ , ed .
d.  e.  soper ( world scientific , 1996 ) [ hep - ph/9601359 ] .
p.  cvitanovi , _ group theory _
( nordita , 1984 ) .
f.  a.  berends and w.  giele , nucl .
b * 294 * , 700 ( 1987 ) ; + m.  l.  mangano , s.  j.  parke and z.  xu , nucl .  phys .
b * 298 * , 653 ( 1988 ) ; + m.  l.  mangano , nucl .
b * 309 * , 461 ( 1988 ) .
g.  t hooft , nucl .
b * 72 * , 461 ( 1974 ) ; nucl .
b * 75 * , 461 ( 1974 ) .
z.  bern and d.  a.  kosower , nucl .
b * 362 * , 389 ( 1991 ) .
z.  bern , l.  j.  dixon and d.  a.  kosower , nucl .
b * 437 * , 259 ( 1995 ) [ hep - ph/9409393 ] . v.  del duca , a.  frizzo and f.  maltoni , nucl .
b * 568 * , 211 ( 2000 ) [ hep - ph/9909464 ] .
v.  del duca , l.  j.  dixon and f.  maltoni , nucl .
b * 571 * , 51 ( 2000 ) [ hep - ph/9910563 ] .
r.  kleiss and h.  kuijf , nucl .
b * 312 * , 616 ( 1989 ) .
n.  e.  j.  bjerrum - bohr , p.  h.  damgaard and p.  vanhove , phys .
lett .   * 103 * , 161602 ( 2009 ) [ 0907.1425 [ hep - th ] ] ; + s.  stieberger , 0907.2211 [ hep - th ] .
e.  plahte , nuovo cim .
a * 66 * , 713 ( 1970 ) .
b.  feng , r.  huang and y.  jia , phys .
b * 695 * , 350 ( 2011 ) [ 1004.3417 [ hep - th ] ] ; + y.  x.  chen , y.  j.  du and b.  feng , jhep * 1102 * , 112 ( 2011 ) [ 1101.0009 [ hep - th ] ] . s .- h .  h. tye and y.  zhang , jhep * 1006 * , 071 ( 2010 ) [ 1003.1732 [ hep - th ] ] .
f.  a.  berends , r.  kleiss , p.  de causmaecker , r.  gastmans and t.  t.  wu , phys .
b * 103 * , 124 ( 1981 ) ; + p.  de causmaecker , r.  gastmans , w.  troost and t.  t.  wu , nucl .
b * 206 * , 53 ( 1982 ) ; + r.  kleiss and w.  j.  stirling , nucl .
b * 262 * , 235 ( 1985 ) ; + j.  f.  gunion and z.  kunszt , phys .
b * 161 * , 333 ( 1985 ) ; + z.  xu , d.  h.  zhang and l.  chang , nucl .
b * 291 * , 392 ( 1987 ) ; + r.  gastmans and t.  t.  wu , _ the ubiquitous photon : helicity method for qed and qcd _ ( clarendon press , 1990 ) .
n.  berkovits and j.  maldacena , jhep * 0809 * , 062 ( 2008 ) [ 0807.3196 [ hep - th ] ] . n.  beisert , r.  ricci , a.  a.  tseytlin and m.  wolf , phys .
d * 78 * , 126004 ( 2008 ) [ 0807.3228 [ hep - th ] ] . j.  m.  drummond , j.  m.  henn and j.  plefka , jhep * 0905 * , 046 ( 2009 ) [ 0902.2987 [ hep - th ] ]
. t.  bargheer , n.  beisert and f.  loebbert , to appear in j.  phys .  a [ 1104.0700 [ hep - th ] ] . t.  bargheer , n.  beisert , w.  galleas , f.  loebbert and t.  mcloughlin , jhep * 0911 * , 056 ( 2009 ) [ 0905.3738 [ hep - th ] ] .
n.  beisert , j.  henn , t.  mcloughlin and j.  plefka , jhep * 1004 * , 085 ( 2010 ) [ 1002.1733 [ hep - th ] ] . l.  f.  alday , d.  gaiotto , j.  maldacena , a.  sever and p.  vieira , jhep * 1104 * , 088 ( 2011 ) [ 1006.2788 [ hep - th ] ] .
d.  gaiotto , j.  maldacena , a.  sever and p.  vieira , 1102.0062 [ hep - th ] . k.  o.  mikaelian , m.  a.  samuel and d.  sahdev , phys .  rev .
lett .   *
43 * , 746 ( 1979 ) .
d.  zhu , phys .
d * 22 * , 2266 ( 1980 ) ; + c.  j.  goebel , f.  halzen and j.  p.  leveille , phys .  rev .
d * 23 * , 2682 ( 1981 ) .
d.  binosi and l.  theussl , comput .
commun .  * 161 * , 76 ( 2004 ) [ hep - ph/0309015 ] ; + d.  binosi , j.  collins , c.  kaufhold and l.  theussl , comput .
commun .   * 180 * , 1709 ( 2009 ) [ 0811.4113 [ hep - ph ] ] .
j.  a.  m.  vermaseren , comput .
commun .  * 83 * , 45 ( 1994 ) .