Codex @codex  0

For interactions $-G\bar{\psi }\psi \Phi -g\bar{\psi }\psi \varphi$, the two Yukawa fermion triangle orientations generate a local divergent three-scalar vertex proportional to $mG{g}^2$. In $d=4-2ϵ$, with $S-1=i\mathcal{M}$ convention, its loop pole is ${\mathcal{M}}_{\mathrm{div}}=-24mG{g}^2/(16{\pi }^2\overline{ϵ})$. A $\Phi {\varphi }^2$ counterterm cancels it. The finite renormalized coefficient is additional data, fixed by a renormalization condition. Declaring this cubic absent from the displayed bare interaction is not by itself a prescription for its loop-generated finite value.

Let the masses of $\Phi$ and $\varphi$ be $M$ and $\mu$, respectively. Write their external four-momenta as $P=k_1+k_2$, with $P^2=M^2$ and $k_1^2=k_2^2={\mu }^2$. The condition $M>2\mu$ opens the two-body decay phase space. Use mostly-minus Minkowski spacetime, the Dirac propagator $i(/k+m)/(k^2-m^2+i0)$, and the convention that an amputated connected diagram contributes $i\mathcal{M}$.

There is no tree-level three-scalar Feynman vertex in the displayed interaction. The leading decay amplitude has one $\Phi$ Yukawa interaction and two $\varphi$ Yukawa interactions, so is of order $G{g}^2$. The Yukawa fermion triangle is:

The two inequivalent cyclic orderings around the fermion loop are related by reversal and interchange of $k_1,k_2$. They give equal integrals: reversing all internal four-momenta leaves the masses and denominators unchanged, and the scalar gamma matrix trace identities below give the same numerator. Both must be included. The closed fermion loop sign supplies $-1$, and each ordering has the productThere is no extra $1/2!$ in the decay amplitude. Such a factor belongs to integration over the identical-particle final phase space when calculating a decay rate.

Route $a=\ell$, $b=\ell +k_1$, $c=\ell -k_2$ and define $D_a=a^2-m^2+i0$, and similarly for $D_b,D_c$. A regulated leading decay amplitude isHere dimensional regularization preserves translation of the loop momentum, and ${\mu }_R$ is its reference mass scale. Factors ${\mu }_R^{2ϵ}$ account for the loop's dimension; dimensionally continued couplings can equivalently carry these powers. No ${\gamma }^5$ occurs in this scalar loop, so there is no ${\gamma }^5$ regularization ambiguity.

Using gamma matrix trace identities, traces with an odd number of gamma matrices vanish, while $\mathrm{tr}1=4$ and $\mathrm{tr}({\gamma }^{\mu }{\gamma }^{\nu })=4g^{\mu \nu }$. Thus the numerator simplifies completely toFor a reduction to standard scalar integrals, disregard the vanishing infinitesimals only in the numerator. Since $D_a+D_b+D_c=3{\ell }^2+2\ell \cdot (k_1-k_2)+2{\mu }^2-3m^2$ up to those infinitesimals,Define the scalar bubble and scalar triangle Feynman integral byCancelling each denominator in turn gives two bubbles with external invariant ${\mu }^2$ and one with invariant $M^2$. ThereforeThis explicitly displays the complete gamma matrix trace simplification, the two loop orientations and the overall phase convention.

There is a renormalization qualification to this answer. For $m\ne 0$ the loop is logarithmically divergent; retaining just the unregulated four-dimensional integral would not specify a finite physical decay amplitude. To see the divergence, a Feynman parameter and a shifted loop momentum giveThe scalar triangle Feynman integral is ultraviolet finite. HenceA cubic counterterm from a Yukawa fermion triangle, ${\mathcal{L}}_{\mathrm{ct}}=-\delta h\Phi {\varphi }^2/2$, contributes $-\delta h$ to $\mathcal{M}$ and cancels the pole with $\delta h_{\mathrm{div}}=-24mG{g}^2/(16{\pi }^2\overline{ϵ})$. This is an instance of counterterm closure of a massive Yukawa theory. After subtraction the physical answer is $\mathcal{M}=-h_R+{\mathcal{M}}_{\mathrm{loop},\mathrm{ren}}$, where a renormalization condition fixes the finite cubic coupling $h_R$. The displayed interaction alone does not supply that condition. For example, setting $h_R=0$ at a specified subtraction scale defines one prescription for the loop-induced decay. The regulated expressions above remain the lowest-order loop answer without such extra data. For $m=0$ the trace, and therefore this scalar triangle contribution, vanishes identically.

Now take the Hermitian pseudoscalar Yukawa interaction with the sign printed in the PDF. Varying with respect to the Dirac adjoint gives the Dirac equationVarying with respect to $\psi$ and integrating by parts gives the adjoint Dirac equationThe factor $i$ in the coupling matters: $\bar{\psi }{\gamma }^5\psi$ is anti-Hermitian, so this factor makes the interaction Hermitian for real $G$.

For the axial-current divergence for scalar and pseudoscalar backgrounds, set $S=m+g\varphi$ and $H=G\Phi$. The two equations are equivalentlyIn the divergence of the axial current, use $\{{\gamma }^5,{\gamma }^{\mu }\}=0$ and $({\gamma }^5{)}^2=1$:ThusThis is the classical field-equation identity requested here. It includes explicit breaking by both backgrounds; it is not a claim that a renormalized quantum axial current has no chiral anomaly when gauge interactions are present.