Propagation of angular momentum in charged pion decay and related processes

There are confusions about angular momentum propagation in scattering or decay processes involving the transition between particle systems that appear to transform differently under Lorentz transformations. This paper provides an analysis of the transformation properties of the states and interactions for a few typical processes within the standard model of particle physics, and performs explicit calculations showing how angular momentum transfers in these processes. We shall show that angular momentum conservation is completely consistent with the calculation in quantum field theory.


I. Introduction
Within the standard model (SM), the decay of a charged pion 1 (take π − to be definite) into a lepton l − and an anti-neutrino νl occurs through weak interactions mediated by a W boson.This process is described in Fig. 1.The corresponding matrix element The decay of a charged pion into a lepton and an anti-neutrino: π − → l − νl .The light blob represents purely perturbative interactions in the graph, while the dark blob also includes non-perturbative contributions. is where α is the fine structure constant, V ud the "ud" element of the quark mixing matrix, and θ W the weak mixing angle.The definition of the amplitude M follows the convention of Ref. [1].The pion and leptons are in "in" and "out" states respectively, and the fields between the states are all renormalized.The polarizations of the fermions 2 are not labeled explicitly.The V − A quark current is defined by Note that the decay occurs inevitably via the interaction term W + µ (x)J µ+ q (x), which is explicitly extracted in Eq. 1.The insertion of the vacuum states indicates then that the fields can only be contracted in such a way as to maintain a factorized form for the transition amplitude represented by the graph in Fig. 1.The first factor of the integrand in Eq. 1 gives the vertex function of two external fermions and an off-shell propagator [2,3] with all order corrections including the self-energy graphs. 3These are denoted by the blob in light grey in Fig. 1.The second factor has a simple tensor structure determined by the properties of the quark current and of the states under Lorentz transformations 4 and translations: where f π is the pion decay constant (in the convention of the particle data group [4]).It receives contributions from perturbative and non-perturbative interactions within the π − induced by the quark current J µ+ q (x).These interactions result in the decay of the pion and are represented by the blob in dark grey in Fig. 1.
Perturbative calculations can be done for the branching ratio R π e/µ ≡ Γ(π − → e − νe )/Γ(π − → µ − νµ ), which is independent of f π .The agreement between the SM prediction of R π e/µ (with radiative corrections) and measurements has reached a level within 0.5%, and is limited by the accuracy of the experiment [5].This provides a good precision test of the SM and perturbative calculations in quantum field theory.Furthermore, f π can be calculated in lattice QCD (see Ref. [4], Table 72.1).It can also be extracted from experiments once the W decay subprocess (i.e., the first factor of the integrand on the r.h.s. of Eq. 1) is computed (with parameters such as V ud taken from other measurements).Comparison of results from both approaches shows an agreement at the percent level [4].This is also the level of uncertainty of the Lattice QCD result.
On the other hand, measurements of the angular distributions of pion decay products in the pion rest frame were done in many early experiments.While most of them obtained isotropic results (e.g., Ref. [6][7][8]), small asymmetry were found in some of the studies (e.g., Ref. [9] ).There are no reports of significant deviation from a uniform angular distribution in more recent experiments of pion decays [10,11].The limits for the sizes of possible Lorentz violating terms are set to the level of 10 −5 ∼ 10 −4 relative to the SM terms [12][13][14].
It is a consensus that no anisotropic angular distributions can be deduced from the charged pion's SM decays, since the pion is a pseudoscalar and the interactions involved are all Lorentz invariant.Studies of possible Lorentz violating effects all seek to find the sources for Lorentz violation from extensions of the SM [12][13][14][15][16].It may be conceptually confusing, however, that in the SM a spin-0 initial state propagates via a spin-1 vector boson, and remains isotropic before being measured.This is usually explained by noting that there is a longitudinal polarization mode of the massive vector boson that behaves like a scalar with 0 angular momentum [1,17,18].In Ref. [19] the decay is considered as occuring via a spin-0 goldstone boson to avoid non-conservation of angular momentum.Nevertheless, it is still helpful to analyze the decay and other related processes using symmetry principles to show explicitly how angular momentum conservation manifests at the amplitude level.
In the following we shall first explore the rotational and Lorentz symmetries in production processes of a vector boson, and then illustrate the results based on symmetry considerations by simple tree level calculations.Consequences of the results in experiments will then be discussed.A summary and discussion of the results will be given at the end.

II. Angular momentum transfer in the production of a vector boson
It is helpful to first work out the angular momentum of a vector boson (vector current) created from a collision of two spin-1/2 fermions.This is somewhat more complicated than the case of charged pion decay, but it provides a richer picture as to how angular momentum is transferred between different parts of an amplitude.To make the discussion general, we start by considering the matrix element where J µ denotes the SM current that tranforms like a vector, axial-vector, or a combination of them, e.g., J µ = Ψ e γ µ Ψ e , Ψ e γ µ 1−γ 5 2 Ψ ν , etc.The momenta and helicities of the fermions are indicated in the state ket.This matrix element represents the all order contribution to the fermion-vector vertex shown in Fig. 2 (a) (with the vector boson part excluded, as indicated by the vertical bar).
The initial state is made up of two fermions with definite momenta and helicities, as indicated in the state ket in Eq. 4. This is not an eigenstate of the angular momentum.To analyse the angular structure of the process, we expand the initial state as a combination of states with definite J, M , λ a , λ b , P + ≡ p a + p b , and 20,21].J and M are respectively the quantum numbers for the total angular momentum of the fermions, and for its projection on the 3-axis in the center-of-mass (CM) frame of the fermions, in which the 3-direction is defined by the direction of the 3-momentum of a: p ≡ p a = −p b .It is convenient to work in this CM frame and express the expansion as [20,21]  FIG.2: Vertices that produce a vector boson in the processes (a): fermion-fermion scattering, and (b): decay (of e.g., π ± , Ψ, etc.).The light blob represents purely perturbative interactions in the graph, while the dark blob also includes non-perturbative contributions.
where P + = 0 is suppressed, and E + is replaced by |p| in the states on the r.h.s. of Eq. 5. J runs over the integers 0, 1, . ... There is no sum over M , because the 3-component of the angular momentum coincides with the difference of the helicities.For an anti-fermion, the physical spin in the direction of its momentum is opposite to its helicity, so an extra minus sign 5 is to be included in its helicity λ a or λ b when obtaining M = λ a − λ b .It is shown [20,21] that C J = (2J + 1)/(4π) in order to be consistent with the normalization of the momentum and angular momentum eigenstates: Note that Eq. 7 is expressed in the CM frame of P + , in which P +′ may be non-zero.For a non-zero P +′ , Eq. 7 gives 0 irrespective of the 3-direction chosen in the CM frame of P +′ to define M ′ .Whereas if P +′ = P + = 0, both M and M ′ are defined with respect to the direction of p.However, in Eq. 7, M and M ′ are not restricted to M = λ a − λ b and M ′ = λ ′ a − λ ′ b , since the states are not restricted to those in the expansion in Eq. 5.The matrix element in Eq. 4 can also be expanded as or, for a full process with the final state |f , We are particularly interested in analysing the contribution to the matrix element in Eq. 8 from each value of J.A vanishing contribution from a specific J signifies that the propagation of this angular momentum eigenvalue is forbidden.To determine the allowed J values we should consider the transformation properties of both the operator J µ and the state |J, M = λ a − λ b , λ a , λ b , |p| .The angular momentum eigenstates with fixed J furnish an irreducible representation under the spatial rotation R(α, β, γ), which is a 3 by 3 matrix that rotates the spatial coordinates and is parameterized by the Euler angles α, β and γ: Note that λ a , λ b , and |p| are invariant under rotations.M (on the l.h.s.) is fixed to λ a − λ b , but is changed by a rotation to M ′ (on the r.h.s.).The subscripts of the Wigner D-matrix D j m ′ ,m (The explicit form can be found in e.g., [20,22]) take the values −J, −J + 1, . . ., J. (All helicity values λ a and λ b are allowed as long as −J ≤ M = λ a − λ b ≤ J).The current operator J µ (0) is a 4-vector that transforms according to under an arbitrary Lorentz transformation Λ (The x dependence of J µ is removed by translation in Eq. 4 so that the transformation of the current is simpler).On the r.h.s. of Eq. 11 µ and ν label respectively the row and column of the matrix Λ −1 .This is equivalent to the self-representation of the Lorentz group, which can be seen from the property for an arbitrary Lorentz transformation matrix Λ Therefore J µ (0) furnish a 4-dimensional irreducible representation of the Lorentz group.However, upon restriction to its spatial rotation subgroup SO(3), the current operators are no longer irreducible.Clearly, J 0 is invariant under rotations, while J 1 , J 2 , and J 3 are mixed.This reduces the 4 components to an SO(3) singlet labeled by J = 0 and a triplet in the SO(3) self-representation labeled by J = 1: As representations of the SO(3) group with the same dimension are equivalent, the D J=0 and D J=1 matrices above are equivalent to the representations D 0 and D 1 in Eq. 10, respectively.Hence the only non-vanishing contribution in Eq. 8 is from J = 0 and J = 1 according to Wigner-Eckart theorem. 6This is not exactly right because they could still vanish due to constraint from further symmetries (See section IV for an example), but at least all final states with J ≥ 2 cannot be produced via a vector boson coupling to the vertex in Fig. 2 (a).
One may want to also investigate the transformation of J µ (0) and |J, M = λ − λ ′ , λ, λ ′ , |p| under the full Lorentz group, rather than its subgroup.However, it is not straightforward to see the transformation law of the angular momentum eigenstates in this case, except that |J = 0, M = λ − λ ′ = 0, λ, λ ′ , |p| will no longer be invariant (a boost changes the eigenvalues P + and E + ).
The analysis above can be carried over immediately to the case of the charged pion decay (or more generally, charged pseudoscalar meson decay) via production of a W boson, as in Fig. 2 (b).The initial state pion is a pseudoscalar with J = M = 0, and the amplitude 0|J µ+ q (x)|π − (k) in Eq. 3 is non-zero from the same reasoning as given above of the J = 0 contribution being non-zero in two fermion scatterings.

III. Amplitudes with angular momentum eigenstates
Now we derive formulae for practical calculations of the amplitudes on the r.h.s. of Eq. 9.The spatial part of p a in Eq. 5 can be rotated by R(φ, θ, 0) to a direction determined by the azimuthal angle φ and polar angle θ, with respect to the original axes.Denoting the initial state under this rotation by we obtain the expansion of the momentum eigenstates in an arbitrary direction by applying the rotation on Eq. 5: The inverse expansion can be obtained using the properties of the D-matrices: This can be used to compute perturbatively the amplitude for the production of a vector boson from a state with a particular J: or, for a full process with the final state |f , The Feynman rules for the momentum eigenstates on the r.h.s. of Eqs. 20, 21 are well known.In the following we shall compute a few simple processes to show the propagation of J values.
IV. Example processes FIG. 3: Example 2 to 2 processes at tree level with the production and decay of a vector boson, where the fermion polarizations are defined by their helicities given in the main text.
We consider 2 to 2 tree level fermion scattering processes in s-channel mediated by a vector boson.The two processes in Fig. 3 will be treated in sequence.In these calculations, we would like to see explicitly the contribution from the J = 0 and J = 1 eigenstates, in accord with the generic result obtained in section II.Therefore the initial and final state fermions are taken to have definite helicities such that the angular structure of the processes are easily revealed.Specifically, we require that the two fermions in the initial state have the opposite spin components along the direction of the incoming beams.This is necessary to ensure that the J = 0 partial wave is non-zero in the expansion of the momentum eigenstates in Eq. 5.In contrast, if the fermions have the same spin along the beam, then the total spin of the system is S = 1, which cannot be canceled by the orbital angular momentum to give J = 0 because the orbital angular momentum is perpendicular to the beam direction.In this case, the partial wave expansion only contains terms with J ≥ 1.
Here the subscript "L" for all the four fermions refers to the helicity configuration λ a = −λ b = λ c = −λ d = −1, rather than definite chiralities.In other words, we use "L" and "R" to describe the relation between the direction of a particle's physical spin and the direction in which it moves.This choice of polarizations is not contradicting to the fact that the anti-neutrino can only be right-handed in chirality, as long as it is a massive particle.Recall that when obtaining M = λ a − λ b we give an extra minus sign to λ b for the anti-neutrino, so that in this case M = 0 instead of −2.In the massless limit, the helicity and chirality eigenstates coincide, but this amplitude vanishes by helicity conservation.However, we shall keep the fermion masses so that the J = 0 mode is allowed, as will be seen below.Recall that we are examining the propagation of each J-mode in principle rather than evaluating the practical importance of the contribution from each J (which will be discussed in section V).
The initial momenta are where m a ≡ m e , and m b ≡ m νe .The tree level vertex for the subprocess e − L (p a )ν Le (p b ) → W is (excluding the W propagator and an overall factor −i √ 2πα/ sin θ W at the vertex) where J µ (0) = Ψ νe γ µ 1−γ5 2 Ψ e , and the particle species are omitted in the states and spinors.Here we follow the method in Ref. [1] to compute the amplitude for polarized states and employ the Weyl representation for the Dirac γ-matrices: in which each element is a 2 by 2 matrix and σ i are the Pauli matrices.The corresponding spinors for the spin-up electron and spin-down anti-neutrino (both with respect to the 3-direction) are Inserting these expressions into Eq.24 gives The vertex for The term in the square brackets is of the same form as in Eq. 24, except the momenta p a = −p b are rotated from the 3-direction to p c = −p d .This term must be a 4-vector that transforms the same way as p a and p b under the rotation Λ(R) that rotates them to p c and p d .This simply follows from the transformation of the vertex function which gives 7 immediately at tree level where E a = E c and E b = E d .Now restore the couplings (and other overall factors) at the vertices and contract with the W propagator in unitary gauge8 to get the amplitude for the full process where M W is the mass of the W boson.The dependence of the result on the external momenta is only through the angle θ between the initial and final fermion beams.
It is straightforward to compute the amplitude with an angular momentum eigenstate using Eq.21, in the form The matrix element on the r.h.s. of Eq 34 is computed with the initial and final states where on the r.h.s. of Eq. 37 (and below), we suppress the factor i(2π) 4 and the delta function that implements the overall momentum conservation.N ab represents the coefficient in Eq. 33 multiplying the square bracket, and θ ac is the angle between p ′ a and p c : There is no contradiction between Eq. 37 and Eq. 6, because Eq. 37 is a product of "in" and "out" states that makes a transition amplitude in the presence of interactions, while Eq. 6 is a normalization condition when both states are "in", "out", or free states.
We do not label these states with "in", "out", or "free" since they can be distinguished by context.Inserting the explicit expressions for the D-matrices [20,22], the angular integration of Eq. 34 can be trivially performed to give Hence, both J = 0 and J = 1 partial waves contribute to the amplitude Eq. 33.9 Indeed, plugging the result of Eq. 39 into the expansion Eq. 9 reproduces Eq. 33 immediately.Note the J ≥ 2 contribution vanishes.This can be seen from the form of the terms on the r.h.s. of Eq. 37.Those terms proportional to 1 and cos θ ′ can be expressed by the matrix elements D J=0 00 and D J=1 00 , which are orthogonal to the D J * 00 with J ≥ 2 in the integration in Eq. 34.The other terms are proportional to eigher cos φ ′ or sin φ ′ , which vanish after the integral over φ ′ .This is consistent with our general argument that states with J ≥ 2 transform differently from the transformation of the current under rotations.Note that the argument was made for the vertex in Fig. 2 (a), with all order contributions.However, the vanishing of the vertex must be order by order since the perturbation series is zero for arbitrary values of the coupling constants involved in the interactions.
A calculation similar to the one carried out above gives where m a = m e and m c = m µ .In place of Eq. 39 we have where now N ab stands for the factor multiplying − cos θ in Eq. 40.This result shows that only the J = 1 mode contributes to the angular distribution of the final state, which by itself reproduces Eq. 40.The amplitudes with J ≥ 2 are zero by the similar arguments given above.The new feature compared to the e − -ν e scattering is the vanishing of the J = 0 result, which is not a coincidence and has a deeper origin.
To see this we note that a positronium state with definite spin S and orbital angular momentum L is an eigenstate of the charge conjugation operator C, with the eigenvalue η C = (−1) L+S [17,23].This property can be carried over to analyzing our process with an unbounded electron-positron pair.The state |J, M, λ a , λ b , p with J = M = 0 has a definite J rather than L or S.However, the rule of adding two angular momenta requires that J = 0 can only be obtained by combining states with definite L and S that satisfy L = S.For all these states, (−1) L+S = 1, and |J = 0, M = 0 is therefore a state that is even under charge conjugation: where the quantum numbers other than J and M are suppressed in the state.Since the electromagnetic current is C-odd, we have which merely reflects the fact that J µ and |J = 0, M = 0 belong to different representations of the charge conjugation operation (one C-even and the other C-odd).So the J = 0 mode does vanish. 10Apparently, this argument cannot be applied to the case of e − -ν e scattering or charged pion decay, because the initial state of either process has a net electric charge that is reversed by C, and is therefore not an eigenstate of C. The initial state must then have a non-zero projection onto the C-odd space of the states, such that the amplitude with the current is not restricted to zero by charge conjugation.

V. Phenomenology
Given a non-zero J = 0 mode carried by a massive gauge boson in processes from previous sections, the immediate question about the role of this mode arises when determining the spin of the vector boson from experiments.For instance, the spin of the W boson was determined by measurement of the angular distribution of the electrons or positrons from the decay W → eν e in the experiment carried out by the UA1 collaboration [24][25][26], in which W bosons were produced via Drell-Yan like processes in proton-antiproton collisions.The measured unpolarized 11 distribution is consistent with the shape (1 + cos θ) 2 , where θ is the angle between the directions of the electron (positron) and the proton (antiproton) in the rest frame of the W boson.The shape of the distribution can be traced back to the Wigner D-matrix D 1 * 1,1 (R(φ, θ, 0)) = e −iφ (1 + cos θ)/2.It shows that the W boson is produced in a state with J = M = 1.
The scalar mode 12 of the W boson contributes an isotropic term to the distribution of the decay product (see Eq. 39).This would in principle distort the distribution deduced from the vector mode of the W.However, the V − A form only allows left-handed fermions and right-handed antifermions to participate in the interaction in the massless limit, in which case the propagation of the scalar mode is forbidden.In fact, it is straightforward to show that a fermion (antifermion) with a righthanded (left-handed) spin would introduce a factor E − |p| = m/ E + |p| to the scattering amplitude, where E, p, and m are the energy, 3-momentum, and mass of the fermion (antifermion) (See Eqs.28 and 31 for instance).If m ≪ E, this gives a suppression factor m/(2E) (known as helicity suppression) relative to the amplitude that is not suppressed (e.g., the one proportional to D 1  1,1 (R(φ, θ, 0)) above).Here E is of the order of the parton level CM energy for the scattering.For an estimate with the UA1 experiment we could take E ∼ 100 GeV, and consider the scalar mode that propagates via the process d R ūR → W → µ − R νµR (The process with left handed spins is much smaller since the final state suppression factor is from the mass of the neutrino rather than from that of the muon.Also we consider the muon instead of the electron in the final state since it has a larger mass).The suppression factor at the cross section level is ∼ O(( 16 .This is far too small to be observed.The statistical uncertainty of the measured angular distribution in the UA1 experiment is as large as 50%.For total cross section and forward-backward asymmetry observables, the accuracy improves to ∼ 10%.More recent measurements at the Fermilab Tevatron [27,28] and the LHC [29,30] reached a percent or even per-mille level accuracy, but are still far from being sensitive to the mass-suppressed modes.
To further reduce the mass-suppression effect on the scalar mode, we might consider the s-channel single top production, so that the final state suppression factor comes from the top quark mass.However, not much could be done about the initial state.
Even if we consider the "sea-sea" contribution to the production of the W boson and take the charm quark from an incoming beam, we get a factor ∼ O(( mcmt 4E 2 ) 2 ) ∼ 10 −6 , where E ∼ 300 GeV for the Tevatron experiment.This will further be suppressed by the sea quark distributions.Also the identification of the top final state is much more difficult than the leptonic decay product of the W [31][32][33][34].
At this point it appears hopeless to probe the J = 0 state of a massive vector boson produced via Drell-Yan (like) processes. 13t is out of the scope of this study to perform a scan of processes, but here we leave the possibility open that this scalar state could manifest in some other production modes of the vector bosons, in addition to the decay processes of e.g., the charged pions.

VI. Summary and discussion
In this study, we analyze the propagation of angular momentum via a vector boson in several processes by expanding the initial states into a superposition of angular momentum eigenstates of the system.Through a general discussion based on symmetry principles and two practical calculations, we are able to show that both the singlet (J = 0) and triplet (J = 1) representations of the rotation group can be transferred to the final state via an intermediate vector boson, unless constraints from further symmetries are present to forbid it.However, we have shown that the scalar state of a massive vector boson is produced with an extremely small rate that it has no impact on existing experiments at hadron colliders.The decays of several spin-zero particles, such as the charged pions, are the only obvious phenomena we have found related with the scalar state of gauge bosons.
The crucial point to understand the J = 0 case is that an initial state that is invariant under spatial rotations is not a singlet under a generic Lorentz transformation.The J = 0 state can be changed by boosts in any direction and gives a non-zero matrix element for producing a Lorentz vector current to which a massive gauge boson couples.Hence there is no Lorentz violation in the production of an off-shell gauge boson from such a state via interactions that are Lorentz invariant.The spin of a vector boson as a virtual particle has to be measured via its decay product and is not necessarily equal to its spin as a free field.
On the other hand, an on-shell vector boson should only be in a state with spin 1, and this should also be reflected at the amplitude level.We have shown that a state with J = M = 0 can propagate though the µ = 0 component of a vector current, as appears in Eq. 8.However, this fact alone does not ensure the full amplitude with a J = M = 0 initial state is none zero, since the initial vertex function will be contracted with the rest part of the amplitude via the propagator of a gauge boson.For instance, we may think of the decay amplitude for π − → l − νl at leading order.After accounting for Eq. 3, the amplitude takes the form where the tensor structure of the propagator in unitary gauge is obtained on the W mass shell by summing over three physical polarization vectors of the W boson These polarization vectors are purely spatial in the W rest frame and make an SO(3) triplet under rotations.Therefore the contraction of the propagator with k µ (whose spatial components are all zero) in Eq. 44 gives Of course, experimentally the pion mass m π − is much smaller than M W .But as M W is a free parameter of the SM, we allow k 2 to get on the W mass shell (i.e., m π − = M W ) for the purpose of checking the consistency of the calculation.The infinitesimal imaginary part of the denominator in Eq. 46 then ensures that the contraction is zero.This structure of contraction persists in loop corrections of the decay amplitude to all orders.Therefore an on-shell W boson indeed forbids the propagation of a J = 0 state.In the off-shell case, however, the state of the W boson cannot be described by the polarization vectors in Eq. 45 that furnish a D 1 representation of the rotation group, and hence allows for a non-zero decay rate.This discontinuous transition of the decay amplitude to zero at the W mass shell shows clearly the consistency of SM dynamics with angular momentum conservation 14 .

VII. Acknowledgement
This work was supported by National Natural Science Foundation of China (12105068) and Hangzhou Normal University Start-up Funds.