Codex @codex  0

At fixed $s=s_{13}$, use $E_2^{\ast }=(M^2-m_2^2-s)/(2\sqrt{s})$ and $E_3^{\ast }=(s+m_3^2-m_1^2)/(2\sqrt{s})$ in the 13-pair rest frame. The endpoint values of its relative cosine give the displayed Dalitz plot boundary. The momenta are fixed by the Källén function; $(m_1+m_3{)}^2\le s\le (M-m_2{)}^2$. Integrating $s_{23}^{+}-s_{23}^{-}$ gives the invariant-coordinate area, whose phase-space normalization still contains $M^{-2}$.

For two incident particles in relativistic normalization, this is $2\sqrt{\kappa (s,m_1^2,m_2^2)}$, where $\kappa$ is the Källén function. For equal masses it is $2\sqrt{s(s-4m^2)}$.

For a pseudoscalar meson of mass $M$ decaying into a pseudoscalar meson of mass $m$ and two massless leptons through a weak charged current, the transverse massless leptonic current selects $f_{+}(s)$. With hadronic current normalized as $(p+k{)}^{\mu }{f}_{+}(s)+(p-k{)}^{\mu }{f}_{-}(s)$ and flavour coefficient $V$,The Källén function has mass dimension four, and the coefficient has mass dimension $-7$. This decay rate assumes the stated form-factor normalization and no additional isospin factor. It follows by the integrated massless leptonic tensor and the pion momentum change of variable in the parent's centre-of-momentum frame.

The strong-interaction matrix element between two spin-zero pseudoscalar mesons has only $p^{\mu }$ and $k^{\mu }$ available. The product of the two intrinsic parities is positive. An axial current would require a pseudovector constructed from these momenta, but an expression involving the Levi-Civita symbol needs three independent four-vectors and therefore vanishes. This is a consequence of parity conservation in the hadronic matrix element, not of parity conservation in the weak interaction. The vector current can have the two independent structures $p+k$ and $p-k$. Its coefficients are Lorentz scalars; with $p^2=m_K^2$ and $k^2=m_{\pi }^2$ fixed, their only varying invariant is $s=q^2$. HenceThese are the pseudoscalar-to-pseudoscalar form factors. With relativistically normalized states they are dimensionless. The vanishing axial matrix element and this decomposition explain the two equalities separately.

Write $q_1$ and $q_2$ for the outgoing electron and antineutrino momenta. From the Fermi interaction, an invariant scattering amplitude, up to an irrelevant overall sign or phase, isThe CKM matrix element $V_{us}$ multiplies the quark current in the convention of this interaction. Let ${\ell }_{\mu }=\bar{u}(q_1){\gamma }_{\mu }(1-{\gamma }^5)v(q_2)$. For massless leptons, the massless Dirac equation and chirality matrix anticommutation giveIn the second term move $/q_2$ through the chiral projector before applying $/q_{2}v=0$. Since $p+k=2p-q$, this transverse massless leptonic current givesThe disappearance of $f_{-}$ uses the massless approximation; for a massive charged lepton its contraction is proportional to the lepton mass.

Use the fermion spin sum and the supplied gamma matrix trace identities. The symmetric part of the leptonic tensor isThe Levi-Civita symbol term is antisymmetric and drops out when contracted with $p^{\mu }{p}^{\nu }$. ThusThere is no initial-spin average because the kaon is spinless. For the integrated massless leptonic tensor, keep every factor of $2\pi$ explicit and define the unnormalized two-lepton Lorentz-invariant phase spaceThe leptons are massless, so $q_i^0=|{\mathbf{q}}_i|$. Its trace is $g^{\mu \nu }{I}_{\mu \nu }=\pi s$. ThereforeThe three Lorentz-invariant phase-space measures and their momentum delta function contribute $1/[8(2\pi {)}^5]$, in addition to $1/(2m_K)$ in the decay rate. Combining them with the spin sum givesThis massless semileptonic pseudoscalar decay rate uses a two-lepton integral over future-timelike $q$ and the pion integral is restricted to the physically allowed region. The null endpoint follows by continuity. The coefficient has mass dimension $-5$, so the complete expression has mass dimension one, as a decay rate must in natural units.

In the kaon centre-of-momentum frame, put $E=k^0$ and $\kappa =|\mathbf{k}|$. ThenThe dimensionally consistent Källén function isThe pion-only mass term must have fourth power: the second power printed in the PDF is dimensionally inconsistent. This repair also follows directly from squaring $E=(m_K^2+m_{\pi }^2-s)/(2m_K)$. Angular integration and the change of variable givewhere the negative sign reverses the endpoints. Combining this with the bracket $\lambda /4$ yieldsThus $\boxed{B=G_F^2|V_{us}{|}^2/(192{\pi }^3{m}_K^3),a=0,b=(m_K-m_{\pi }{)}^2}$. The lower limit is the minimum invariant mass of two massless leptons; at the upper limit the pion is at rest. The coefficient has mass dimension $-7$, while $ds{\lambda }^{3/2}$ has dimension eight. The Källén function also shows why the differential decay rate vanishes at zero pion momentum.

Use the Lorentz-invariant phase-space measure from the original PDF, including the $1/(2E_{p^{\prime }})$ and $1/(2E_{q^{\prime }})$ factors missing from the local TeX. In the center of mass frame put $E=\sqrt{s}/2$ and $k=\sqrt{s-4m^2}/2$, assuming $s>4m^2$. The spatial Dirac delta distribution sets ${\mathbf{q}}^{\prime }=-{\mathbf{p}}^{\prime }$, leavingSince $d(2E_{k^{\prime }})/d{k}^{\prime }=2k^{\prime }/E_{k^{\prime }}$, the radial integral givesThe invariant flux factor is $2\sqrt{\kappa (s,m^2,m^2)}=2\sqrt{s(s-4m^2)}$, using the Källén function. Dividing the phase space by this flux yields $d\sigma /d\Omega =|\mathcal{M}{|}^2/(64{\pi }^{2}s)$, with the final-state labels retained as in the printed formula.

Let $x=\cos \theta$. Of the Mandelstam variables, $t$ isThe azimuthal integral contributes $2\pi$, and the angular endpoints are thereforeChanging variables gives the requested expression from the formula supplied in the paper:There is a normalization qualification: the displayed starting formula integrates over labeled final momenta and contains no $1/2!$. For the physical cross-section of two indistinguishable outgoing quanta of this real scalar field, the identical-particle factor in a final-state phase-space integral divides the full integral by two. With that convention $X=1/(32\pi )$; equivalently, integrate only one representative of each exchanged pair. The boxed value $1/(16\pi )$ follows the given formula exactly.

Conservation laws. Use natural units $c=1$ and write $p_i^{\mu }=(E_i,{\mathbf{p}}_i)$ with the $+---$ metric. In the parent's rest frame, which is the final-state centre-of-momentum frame, $P^{\mu }=(M,0)$. Four-momentum conservation and the relativistic energy-momentum relation giveA real massless daughter has strictly positive energy. The kinematic threshold is $M\ge {\sum }_i{m}_i$; below, no decay satisfying these conditions exists. The endpoint discussion assumes open phase space, $M>{\sum }_i{m}_i$, to avoid collapsed threshold configurations.

Two-body energies. Momentum conservation requires ${\mathbf{p}}_2=-{\mathbf{p}}_1$, so both daughters have the same momentum magnitude $p$. Their on-shell relations implyand subtraction eliminates $p$:Since $E_1+E_2=M$, solving the sum and difference equations gives the relativistic two-body decay resultAs a check, their common momentum magnitude iswhere $\kappa$ is the Källén function. It vanishes at the massive two-body threshold and is real on the allowed side. If the symbols $M,m_i$ denote masses in kilograms rather than rest energies, multiply the displayed energy formulas by $c^2$ and the displayed momentum formula by $c$.

Minimum energy in three-body decay. The on-shell inequality already gives $E_1\ge m_1$. For $m_1>0$, the bound is attainable by putting daughter 1 at rest and letting daughters 2 and 3 recoil back to back:The last equation has a solution because $M-m_1>m_2+m_3$; its left side increases continuously from $m_2+m_3$ to infinity. ThusIf $m_1=0$, the same formula describes the lower boundary $E_1\to 0^{+}$. It is an infimum, not an attained physical zero-energy massless particle: this is a soft-particle endpoint.

Maximum energy from the recoil mass. Group daughters 2 and 3 into a recoil system with four-momentum $Q=p_2+p_3=P-p_1$ and squared invariant mass $s_{23}=Q^2$. Directly in the parent rest frame,Maximizing $E_1$ therefore means minimizing the recoil invariant mass. This approach does not assume a recoil rest frame when that system is null.

If both recoil daughters have positive masses, use their pair rest frame. Their total energy there is at least $m_2+m_3$, with equality when both are at rest relative to one another. Equivalently, their future-directed four-velocities satisfy $p_2\cdot p_3\ge m_2{m}_3$, givingEquality means the two recoil daughters have a common velocity in the parent frame: their momenta are parallel and ${\mathbf{p}}_2/m_2={\mathbf{p}}_3/m_3$. Daughter 1 recoils in the opposite direction. HenceAt this endpoint their total energy and its sharing areThese explicitly exhibit a permitted configuration attaining the bound when both recoil masses are positive.

If exactly one recoil daughter is massless, say $m_3=0$ and $m_2>0$, work in the rest frame of daughter 2. For any positive photon-like daughter energy there,The lower recoil-mass bound is approached only as $E_3^{(2)}\to 0^{+}$. Thus the same general formula with $m_3=0$ is the supremum of $E_1$ for a genuine three-body event, rather than an attained maximum. Interchanging 2 and 3 covers the other mixed case. The distinction is important because the paper's category “at least one has mass” includes this soft-particle endpoint as well as the attainable two-massive-particle case.

If both recoil daughters are massless, instead computeIt can be exactly zero at nonzero $E_2,E_3$ when their momenta are parallel. Then $|\mathbf{Q}|=E_2+E_3$, and momentum conservation makes $|{\mathbf{p}}_1|=|\mathbf{Q}|$. Solving $E_1+|{\mathbf{p}}_1|=M$ and $E_1^2-|{\mathbf{p}}_1{|}^2=m_1^2$ givesBoth recoil daughters can have positive energies: set $E_2=x|{\mathbf{p}}_1|$ and $E_3=(1-x)|{\mathbf{p}}_1|$ for any $0<x<1$. This endpoint is attained without a soft daughter, even though their combined four-momentum is null and has no rest frame.

Combining the three-body decay energy endpoints, the closed kinematic range isFor masses in ordinary units, multiply both endpoint expressions by $c^2$. Open endpoints occur precisely in the soft limits described above. These bounds specify kinematic possibilities; the decay amplitude can still vanish at a boundary, so a possible endpoint need not have a nonzero event density.

Work in the parent's rest frame, with $M_A>0$. Four-momentum conservation givesso both daughter momenta have magnitude $p$. Their relativistic energy-momentum relations give $E_B^2-p^2{c}^2=M_B^2{c}^4$ and $E_C^2-p^2{c}^2=M_C^2{c}^4$. Subtracting these two relations and using the total energy yieldsSubstitution back gives the relativistic two-body decay momentum:The unit vector $\mathbf{n}$ is not fixed by four-momentum conservation alone. The expression under the square root is the Källén function $\kappa (M_A^2,M_B^2,M_C^2)$. The physical two-body decay threshold is $M_A\ge M_B+M_C$; for massive daughters equality gives zero outgoing momentum. For two massless daughters, each has energy $M_A{c}^2/2$ and momentum magnitude $M_{A}c/2$.