The canonical anticommutation relations arewith the mixed and equal-type anticommutators zero. Insert the mode expansion of a Dirac field. The phases cancel at equal time, and a givesSince $\bar{\psi }{\gamma }^0={\psi }^{†}$, the requested componentwise answer isThis is the Equal-time canonical anticommutator of a Dirac field: the canonical momentum conjugate to ${\psi }_{\alpha }$ is ${\pi }_{\alpha }=i{\psi }_{\alpha }^{†}$. It is the fermionic counterpart of the canonical field canonical commutation relation, ensuring the correct local canonical quantization. For this first-order Lagrangian density, the conjugate variables are constrained; the displayed relation is the standard result after handling that constraint. It does not assert a nonzero anticommutator between the fields at separated spatial points.

Let $A=\mathbf{p}\cdot \mathbf{\sigma }$ and $E=\sqrt{{\mathbf{p}}^2+m^2}$. The Pauli multiplication identity gives $A^2={\mathbf{p}}^2{I}_2$, and the orthonormal two-component spinors obey ${\sum }_s{\xi }^s{\xi }^{s†}=I_2$. Taking the Dirac adjoint of the supplied Dirac spinor therefore gives, directly,For the antiparticle deduction, specify the standard normalization $v^{r†}(-\mathbf{p})v^s(-\mathbf{p})=2E{\delta }^{rs}$, just as $u^{r†}{u}^s=2E{\delta }^{rs}$. Spanning alone does not specify a normalization: replacing every $v$ by $2v$ preserves the spanning property but multiplies its proposed fermion spin sum by four. The standard normalization is the convention implicit in the mode expansion of a Dirac field with the standard canonical anticommutation relations.

To use completeness correctly, introduce the Hermitian Dirac HamiltonianThe Dirac equation makes $u^s(\mathbf{p})$ and $v^s(-\mathbf{p})$ its eigenvectors with eigenvalues $E$ and $-E$. Their mutual orthogonality and spanning property give complementary orthogonal projectorsMultiply the second equality by $2E{\gamma }^0$ on the right and replace $\mathbf{p}$ by $-\mathbf{p}$. For the positive-energy four-momentum $p=(E,\mathbf{p})$,

Put $Q=M-3m$. For $Q<0$ the decay width is zero by conservation of energy. Just above the three-particle decay threshold, all three daughters are nonrelativistic. The decay amplitude of iv remains finite: the internal scalar propagator denominator approaches $(2m{)}^2-m^2=3m^2$, and the Dirac scalar bilinear does not vanish at rest. Thus, for nonzero fixed $g$ and $\lambda$, the leading suppression comes entirely from three-body threshold phase space.

Use the spatial Dirac delta function to eliminate one daughter momentum. The remaining integral has six independent momentum coordinates and one conservation of energy constraint,Rescaling both momenta by $\sqrt{mQ}$ gives $Q^3$ from the measure and $Q^{-1}$ from the Dirac delta function. The smooth factors $1/(2E)$ do not change the power. Consequentlywith a finite positive leading coefficient when $g\lambda \ne 0$; $\Gamma =0$ for $m\ge M/3$. The $1/2!$ for the identical bosons changes only this coefficient. Equivalently, at fixed $M$ the opening is quadratic in $M/3-m$. If $g\lambda =0$, this leading tree-level Feynman diagram vanishes, so the stated leading coefficient is zero.

The printed integration labels need correction. The parent's four-momentum is fixed, so the two Lorentz-invariant phase-space measures must integrate $p_2$ and $p_3$, rather than the printed $p_1$ and $p_2$. Write $P=p_1=(M,0)$ and $p=p_2$, $k=p_3$. The physical decay width isThe spatial Dirac delta function sets $\mathbf{k}=-\mathbf{p}$. Performing the angular integral in spherical coordinates, with $E=\sqrt{p^2+m^2}$, gives the integrated two-body decay phase spaceThe positive root is $p_{\ast }=\sqrt{M^2/4-m^2}$. Since $d(2E)/dp=2p/E$, the Dirac delta function integration yields ${\Phi }_2=p_{\ast }/(4\pi M)$. ThusThe two final particles have different species, so there is no identical final-state symmetry factor. The result tends to $g^{2}M/(16\pi )$ as $m\to 0$ and vanishes at the two-particle decay threshold.

In the convention of iii, the invariant decay amplitude is ${\mathcal{M}}_{s_1{s}_2}=-g{\bar{u}}_{\chi }^{s_2}(p_2)u_{\Psi }^{s_1}(p_1)$. Including the common overall $i$ in the definition changes no squared decay amplitude. An unpolarized ensemble is an incoherent equal mixture of the two initial spin states; hence take a spin average over $s_1$ and a spin sum over the unobserved final $s_2$. Different final spin states are orthogonal outcomes, so their probabilities, rather than their decay amplitudes, add.

Using the Dirac adjoint and the fermion spin sum givesThe gamma-matrix trace identities used here are $\mathrm{tr}{\gamma }^{\mu }=0$ and $\mathrm{tr}({\gamma }^{\mu }{\gamma }^{\nu })=4{\eta }^{\mu \nu }$. The mass-shell conditions $p_1^2=M^2$ and $p_2^2=p_3^2=m^2$, together with four-momentum conservation, implyThereforeIn the parent's rest frame this also says $E_2=E_3=M/2$.

The momentum space Feynman rules use the scalar propagator and two Dirac propagatorsThe cubic scalar field interaction vertex is $-i\lambda$. Each of the two flavour-changing Yukawa coupling interaction vertices, $\varphi \bar{\Psi }\chi$ and $\varphi \bar{\chi }\Psi$, has factor $-ig{\delta }_{ab}$ in Dirac spinor indices. A fermion flow arrow continues through each interaction vertex, but its species changes. There is no $\varphi \bar{\Psi }\Psi$ or $\varphi \bar{\chi }\chi$ interaction vertex.

Every interaction vertex carries $(2\pi {)}^4{\delta }^{(4)}(\sum p)$ with inward four-momenta. Integrate each independent loop momentum with $d^4\ell /(2\pi {)}^4$, divide by the Feynman-diagram symmetry factor, and include a minus sign and a gamma-matrix trace for each closed fermion loop. The $i0$ is the Feynman i-epsilon prescription.

For an amputated scattering amplitude, incoming fermions supply $u^s(p)$ and outgoing fermions supply ${\bar{u}}^s(p)$; incoming antifermions supply ${\bar{v}}^s(p)$ and outgoing antifermions supply $v^s(p)$. An external real scalar field particle supplies $1$. Use the standard relativistic normalization of a one-particle state; the external factors $1/\sqrt{2E}$ belong to that normalization and the Lorentz-invariant phase-space measure, rather than an additional factor in the invariant scattering amplitude.

Hold the external endpoints fixed and distinguish the three Feynman diagrams by their left-to-right order in the source. Write $S$ for the divisor and $w=1/S$ for the multiplying combinatorial weight. Each line below includes external Feynman propagators; amputation leaves the same $w$.

Left: two-point bubble. There are two cubic interaction vertices and two parallel internal lines. Interchanging those lines preserves the labelled external endpoints, giving $S=2$. Direct Wick contraction counting gives two assignments of the external endpoints to the vertices, $3^2$ choices of their field slots, and $2!$ pairings of the remaining slots:With $\Delta ={\Delta }_F$ and external points $x_1,x_2$, its contribution is

Middle: three-point triangle. The three labelled external endpoints fix their three adjacent cubic interaction vertices, so no nonidentity graph automorphism remains. Counting gives $3!$ assignments of endpoints to vertices, $3^3$ choices of external field slots, and $2^3$ assignments of remaining slots to the two specified neighbours:Its contribution isThe apparent rotational symmetry of an unlabelled triangle does not survive fixing three distinct external endpoints.

Right: two-point graph with a chord. There are four cubic interaction vertices: the left and right attachment vertices $A,B$ and the two chord ends $C,D$. Its internal edges are $AC,AD,BC,BD,CD$. Reflection interchanges $C,D$, fixes $A,B$, and preserves the two external labels. This is the only nontrivial graph automorphism, so $S=2$.

For a separate Wick contraction count, choose the ordered external attachment vertices in $4\cdot 3$ ways and their field slots in $3^2$ ways. Assign their two remaining slots to $C,D$ in $2^2$ ways. At each of $C,D$, choose the slot joining $A$ and then the slot joining $B$ in $3\cdot 2$ ways; the last slots join each other. ThusIts contribution isConsequently the three multiplying weights are $\boxed{(1/2,1,1/2)}$, or divisors $\boxed{(2,1,2)}$. The central chord is a single internal line joining two vertices, not a vertex at the centre of the circle.

For the interacting vacuum state defined perturbatively by the same Feynman i-epsilon prescription, the normalized time-ordered product isEach functional derivative inserts $i\varphi$, accounting for $i^{-m}$. The denominator implements cancellation of vacuum bubbles, while disconnected pieces carrying external insertions still contribute. The lower limit $i=i$ in the printed product is a typographical error; the product here runs from 1 to $m$.

At order $V$ in the cubic scalar field interaction term, the expansion of the generating functional supplies $1/[V!(3!{)}^V]$. The $V!$ counts permutations of indistinguishable integrated interaction vertices, and $3!$ counts permutations of the three field slots at each interaction vertex. Applying the Wick theorem pairs fields into Feynman propagators. Several pairings produce the same Feynman diagram; their number $N_G$ multiplies the expansion coefficient. For fixed labelled external points, the net combinatorial weight isHere $S_G$ is the order of the graph's automorphism group when the external labels are fixed, with permutations of identical parallel lines and reversal of a self-loop included as applicable. This is the Feynman-diagram symmetry factor convention in which a Feynman diagram is divided by $S_G$. Each cubic interaction vertex then contributes $-i\lambda$, each line contributes ${\Delta }_F$, and each internal vertex is integrated over its position. External labelled points are not divided by $m!$.

Keep $\lambda$ for the cubic coupling and $j$ for the external source. Expanding the interaction part of the scalar field path integral yieldsA functional derivative of the source exponential inserts $i\varphi (y)$. In particular, ${\delta }_j(y{)}^3{e}^{i\int j\varphi }=-i\varphi (y{)}^3{e}^{i\int j\varphi }$, so $-i\lambda \varphi (y{)}^3$ can be replaced by $\lambda {\delta }_j(y{)}^3$. The functional derivatives commute and can be taken outside the free path integral. Resumming the exponential series givesThis is a formal perturbation series; a real cubic scalar potential is unbounded and the expression does not assert a stable nonperturbative vacuum state.

To evaluate the free Gaussian path integral, use integration by parts to writeThe inverse kernel isCompleting the square with $\varphi =\chi -K^{-1}j$ changes the exponent to $i\int \chi K\chi /2-i\int j{K}^{-1}j/2$. The translation-invariant functional measure makes the remaining Gaussian path integral precisely $Z[0;0]$, henceFor $ϵ>0$, the exponent includes $-ϵ\int {\varphi }^2/2$, which damps the regulated finite-dimensional integrals. The Feynman i-epsilon prescription displaces the positive-energy pole below the real $p^0$ axis and the negative-energy pole above it, selecting the time-ordered vacuum expectation value. Indeed the Feynman propagator is ${\Delta }_F=iD$. Thus $iϵ$ specifies both convergence and the vacuum state boundary conditions, with $ϵ\to 0^{+}$ taken after integration.

Take parity symmetry in quantum field theory to preserve free time evolution, so the conjugate initial data reflect along with the field. Substitute the mode expansion of a free field into the defining parity operator relation and change the integration variable $\mathbf{p}\mapsto -\mathbf{p}$. Since $E_{-\mathbf{p}}=E_{\mathbf{p}}$, coefficient comparison, or comparison of both the field and its conjugate initial data, gives the intrinsic-even conventionHere parity symmetry in quantum field theory reverses spatial momentum and preserves time. The creation operator formulas follow by taking Hermitian adjoints. The given field relation has fixed the intrinsic phase to $+1$.

For one real scalar field, let $R$ act on one-particle wave functions by $(Rf)(\mathbf{p})=f(-\mathbf{p})$. It is a Hermitian operator with $R^2=1$, so $P_{-}=(1-R)/2$ is the orthogonal projector onto odd wave functions. Its second quantization is the odd-mode number operatorThe full momentum space is integrated in both expressions; the factors prevent double counting the pairs $\{\mathbf{p},-\mathbf{p}\}$. A suitable parity operator isFor verification, $[N_{\mathrm{odd}},a_{\mathbf{p}}]=-(a_{\mathbf{p}}-a_{-\mathbf{p}})/2$. Thus the even combination $a_{\mathbf{p}}+a_{-\mathbf{p}}$ is unchanged, and the odd combination $a_{\mathbf{p}}-a_{-\mathbf{p}}$ changes sign under conjugation by $e^{i\pi N_{\mathrm{odd}}}$, exactly reproducing $a_{\mathbf{p}}\mapsto a_{-\mathbf{p}}$. The same holds for creation operators. The odd-mode number operator has integer eigenvalues, so ${\mathcal{P}}^2=1$, and $\mathcal{P}|0⟩=|0⟩$. Multiplying by an overall constant phase would not change its conjugation action. This is parity as odd-mode number parity.

Use the combinations $a_{\pm }=a\pm b$ at each momentum. Part a gives $[Q,a_{+}]=0$, $[Q,a_{-}]=-2a_{-}$, while their creation operators obey $[Q,a_{+}^{†}]=0$, $[Q,a_{-}^{†}]=2a_{-}^{†}$. Applying b givesThe creation operators transform by the Hermitian adjoint formulas with $e^{+2i\lambda }$. At $\lambda =\pi /2$, both exponentials are $-1$, soSubstituting into the fixed-time mode expansion of a free field of the complex scalar field givesThus $C$ exchanges particles and antiparticles, as required for charge conjugation. To see the charge reversal explicitly, the conserved charge $Q_{\mathrm{em}}={\int }_{\mathbf{p}}(a_{\mathbf{p}}^{†}{a}_{\mathbf{p}}-b_{\mathbf{p}}^{†}{b}_{\mathbf{p}})$ obeys $C{Q}_{\mathrm{em}}{C}^{†}=-Q_{\mathrm{em}}$. The generator $Q$ in the question is distinct from this charge.

The two occupation terms in the Hamiltonian are exchanged, with the same relativistic energy $E_{\mathbf{p}}$, and thereforeIn fact $[Q,H]=0$ because its mixing terms exchange equal-energy modes. With normal ordering, $:Q:=2{\int }_{\mathbf{p}}{c}_{-,\mathbf{p}}^{†}{c}_{-,\mathbf{p}}$ for $c_{-}=(a-b)/\sqrt{2}$. Then $C=e^{i\pi N_{-}}$ fixes the vacuum state and satisfies $C^2=1$. The construction is a bosonic species-swap generator.

For a Hermitian operator $Q$, the operator exponential $C_{\lambda }=e^{i\lambda Q}$ is a unitary operator, with $C_{\lambda }^{†}=e^{-i\lambda Q}$. Differentiation and $[Q,C_{\lambda }]=0$ giveWhen $[Q,A]=cA$, this is the first-order ordinary differential equation $A^{\prime }(\lambda )=icA(\lambda )$ with $A(0)=A$. Multiplication by $e^{-ic\lambda }$ makes its derivative zero, provingEquivalently, the commutator expansion for exponential conjugation has successive nested commutators $c^{r}A$, whose exponential series sums to this expression.

Write ${\int }_{\mathbf{p}}=\int d^{3}p/(2\pi {)}^3$ and $d_{\mathbf{p}}=a_{\mathbf{p}}-b_{\mathbf{p}}$. The nonzero canonical commutation relations arewith cross-species brackets zero. The commutator product rule givesLinearity of the commutator now yields, for the real coefficient $s$,Strictly, the printed $Q$ contains an infinite constant in continuum Fock space. Work first with finitely many modes, or replace it by its normal ordering $:Q:={\int }_{\mathbf{p}}{d}_{\mathbf{p}}^{†}{d}_{\mathbf{p}}$. Their difference is a constant, so every displayed commutator and every conjugation below is unchanged. The latter convention fixes the vacuum state phase without an infinite phase ambiguity.

For a connected amputated Feynman diagram, write $L$ for the number of independent loop momenta. Each loop momentum integration scales as ${\lambda }^4$, each internal scalar propagator as ${\lambda }^{-2}$, and each internal Dirac propagator as ${\lambda }^{-1}$ when both four-momenta and mass are scaled. The Yukawa coupling has mass dimension zero. ThusCounting the two fermion ends and one scalar end at each interaction vertex, and using connectedness, givesConsequently $n_f=n_v$, $L=n_s+1$, andThe requested exponent requires amputation and removal of the overall momentum delta function. An ordinary unamputated connected Green function has an additional $n$ external scalar propagators, giving mass dimension $4-3n$ after stripping that delta function, or $-3n$ if it is retained. For a concrete check, the free two-point scalar propagator scales as ${\lambda }^{-2}$, not ${\lambda }^2$. The exponent $4-n$ is the engineering result for the amputated object. A fixed renormalization scale introduces logarithms; the displayed homogeneity is a naive power counting statement rather than an exact quantum scaling law at fixed subtraction scale.

More generally, with $E_s$ external scalar ends and $E_f$ external fermion ends, the same counting gives the superficial degree of divergenceThis is an upper bound on the overall ultraviolet divergence. Numerator cancellations and Lorentz invariance may lower it; an apparently convergent whole Feynman diagram can still have ultraviolet subdivergences. Local subtraction of those ultraviolet subdivergences is required first.

For pure scalar ends, $E_s=0,1,2,3,4$ can require counterterms. Lorentz invariance reduces the two-point terms to a mass counterterm and wave-function renormalization; derivatives of a single scalar field are total derivatives, and a one-derivative cubic scalar term cannot be a Lorentz scalar. With two external fermions, $E_s=0$ requires their kinetic term and mass counterterm, while $E_s=1$ requires Yukawa coupling renormalization. Four external fermions give $D=-2$, so no independent four-fermion counterterm is required.

A complete counterterm Lagrangian consistent with the species symmetry can therefore be writtenThe vacuum energy density term matters for an unnormalized vacuum functional and cancels from normalized Green functions. The tadpole counterterm can impose $⟨\varphi ⟩=0$. There is no $\varphi \mapsto -\varphi$ symmetry of the massive Yukawa interaction, so linear, cubic and quartic scalar terms cannot be excluded. For instance, in a constant background the divergent fermion loop contribution contains $(m+g\varphi {)}^4$, exhibiting all these powers. This is counterterm closure of a massive Yukawa theory.

The two species have identical couplings, so their common unitary flavour rotations preserve the Lagrangian density and their counterterms can be chosen species independent. Equality of the scalar and fermion masses, however, is not protected: $\delta m_s^2$ and $\delta m_f$ are independent. Including the displayed finite set of operators yields a renormalizable quantum field theory. Restricting the bare Lagrangian density to the terms printed in the question is not closed under renormalization.

Define $\mathcal{M}$ by removing the factor $i(2\pi {)}^4{\delta }^{(4)}(p_1+p_2-p_3-p_4)$ from the connected transition matrix element. The single tree-level Feynman diagram in i givesThere is no identical-particle exchange term because the species differ. Each incoming fermion has two equally populated spin angular momentum states in an unpolarized beam. Their probabilities multiply to $1/4$, giving a spin average over the four initial choices. The unobserved final spin angular momentum states are distinct outcomes, so their probabilities are added by a spin sum, with no final averaging:Taking the complex conjugate of ${\bar{u}}_3{u}_1$ gives ${\bar{u}}_1{u}_3$. The fermion spin sum therefore converts each squared bilinear into a gamma-matrix trace:Using $\mathrm{tr}1=4$, $\mathrm{tr}{\gamma }^{\mu }=0$, $\mathrm{tr}({\gamma }^{\mu }{\gamma }^{\nu })=4{\eta }^{\mu \nu }$, and the mass shell conditions $p_j^2=m^2$, each gamma-matrix trace isThe second equality follows from $q^2=2m^2-2p_1\cdot p_3$; momentum conservation gives the identical result for $p_2,p_4$. HenceThe printed prefactor $g^2$ is incorrect for the stated coupling. Each of the two Yukawa interaction vertices contributes one $g$ to the scattering amplitude, so its square is proportional to $g^4$. As a concrete counterexample, take $g=2$, $m>0$, and forward elastic kinematics $p_3=p_1$, $p_4=p_2$. Then $q^2=0$ and the spin average is $16g^4=256$, whereas the printed expression gives $16g^2=64$.