Codex @codex  0

For noninteracting particles, a many-particle configuration has $E={\sum }_j{n}_j{ϵ}_j$ and $N={\sum }_j{n}_j$. Summing independently over allowed occupations therefore factorizes the grand partition function:At fixed $T,V$, differentiating the original state sum gives ${\partial }_{\mu }\ln \Theta =\beta ⟨N⟩$ and ${\partial }_{\mu }^2\ln \Theta ={\beta }^2\mathrm{Var}N$, soFor one level, write $q=e^{-\beta (ϵ-\mu )}$. Fermions allow $n=0,1$, giving ${\Theta }_F=1+q$ and $\bar{n}=q/(1+q)$. Bosons allow all nonnegative integers, giving ${\Theta }_B=1/(1-q)$ and $\bar{n}=q/(1-q)$, with $q<1$. ThusThe grand-canonical single-level occupation variance follows from the same derivative identity applied to one level:In the high-temperature nondegenerate regime, $\bar{n}\ll 1$ and both approach classical Poisson-like fluctuations, $\mathrm{Var}n\simeq \bar{n}$. At low temperature, nearly filled or empty fermion levels have strongly suppressed fluctuations; only the Fermi-surface window changes occupation. Highly occupied low-energy boson levels can instead have $\mathrm{Var}n\sim {\bar{n}}^2$, showing bosonic enhancement. This last statement depends on occupation and reservoir chemical potential, not temperature alone; grand-canonical condensate fluctuations must not be confused with a fixed-total-number ensemble.

For an infinite cubic potential well of side $L$, the one-particle energy eigenvalues are $E_{n_x{n}_y{n}_z}=E_0(n_x^2+n_y^2+n_z^2)$, where $E_0={\pi }^2{\hslash }^2/(2m{L}^2)$ and each quantum number is a positive integer. The lowest spatial orbital is $(1,1,1)$, with energy $3E_0$ and orbital degeneracy one. The next energy is $6E_0$, with three orbitals given by permutations of $(2,1,1)$.

Spin-zero identical particles are bosons. Both can occupy the lowest orbital, so the two-boson ground state is nondegenerate. The first excited total energy is $9E_0$: one particle occupies the lowest orbital and one an excited orbital. Symmetrizing the two-particle wavefunction gives one state for each of the three excited orbitals, so the first-excited degeneracy is three.

Spin-half identical particles are fermions, so the total state must be antisymmetric. Two particles in the same lowest orbital have a symmetric spatial wavefunction and must form the unique spin singlet state. Thus the fermion ground state is also nondegenerate. At the first excited energy, each choice of excited spatial orbital permits a symmetric spatial state with a singlet state, or an antisymmetric spatial state with a three-dimensional triplet state. ThereforeThe spin multiplicity is essential: counting only the three orbital choices would omit allowed triplet states. The result assumes ordinary spin-statistics, an unsplit spin degeneracy and the hard-wall cube model.

For conventional superconductivity, an attractive interaction pairs opposite-spin fermions. A complex pairing Hubbard–Stratonovich transformation introduces a Cooper pair field $\Delta$. The fermionic fields then appear quadratically, giving the pair-field effective action of a superconductorThe trace convention counts each independent spin-pair block once. The uniform saddle is the BCS gap equation; the fermionic excitations have $E_k=\sqrt{{\xi }_k^2+|\Delta {|}^2}$ in the fully gapped case. A slow-field expansion near the transition yields Ginzburg-Landau theory of superconductivity. Normalize its pair field $\Psi$ so $|\Psi {|}^2$ is an effective pair density, with pair mass $m_{\ast }$ and signed charge $q=-2e$:For $a<0$, the uniform minimum has ${\rho }_0=|\Psi {|}^2=-a/b$, with arbitrary constant phase. Writing $\Psi =(\sqrt{{\rho }_0}+h)e^{i\theta }$ separates the amplitude mode $h$ from the phase. The local potential has curvature term $2|a|h^2$ and coherence scale $\xi =\hslash /\sqrt{2m_{\ast }|a|}$. A neutral condensate has a gapless phase Goldstone mode when suitable dynamics is included.

The charged theory is invariant underThe gauge-invariant condensate momentum is $\hslash \nabla \theta -q\mathbf{A}$. Only $\hslash \nabla \theta -q\mathbf{A}$ and $\hslash \dot{\theta }+q\varphi$ enter physical phase dynamics. At nearly constant amplitude, varying $F$ givesFor a smooth phase, taking the curl and using Ampère's law gives ${\nabla }^2\mathbf{B}=\mathbf{B}/{\lambda }^2$: the Meissner effect screens a static bulk magnetic field over the London penetration depth $\lambda$. In a gauge with $\theta =0$, the phase stiffness is a quadratic term in $\mathbf{A}$, the superconducting version of the Higgs mechanism.

To identify the dynamical mode, a simple lossless long-wave superconducting phase-density action isHere ${\chi }_p$ is the pair compressibility coefficient. Eliminating the longitudinal electric potential givesThus the neutral sound mode becomes a gapped longitudinal plasma oscillation; transverse electromagnetic waves acquire ${\omega }^2=c_{\mathrm{vac}}^2{k}^2+c_{\mathrm{vac}}^2/{\lambda }^2$ in this simple bulk model. This redistribution of the phase degree of freedom into electromagnetic collective modes is the Anderson-Higgs mechanism in a superconductor. Local gauge symmetry is a redundancy, not an independently observable symmetry whose breaking creates another massless particle. ^[ref]

The static functional alone does not determine a frequency or damping rate. The displayed mode dispersions require the added dynamical assumptions; real screening, normal carriers, boundaries and dissipation modify them. The amplitude fluctuation is distinct from the phase mode and is not removed by the gauge choice.

The spin commutation relations give $[S_{\mathrm{tot}}^z,S_n^{\pm }]=\pm S_n^{\pm }$. Using the product rule for a commutator,and the conjugate term also vanishes. The longitudinal field commutes with every $S_n^z$, so $[H,S_{\mathrm{tot}}^z]=0$. The Hamiltonian is invariant under the global spin rotation $U(\theta )=e^{-i\theta S_{\mathrm{tot}}^z}$: its continuous symmetry group is $U(1)$. Simultaneous energy and total-magnetization eigenstates can be chosen, and evolution preserves each fixed-magnetization sector. Local magnetization can still flow between sites; generally $[H,S_n^z]\ne 0$. In the fermion representation the conserved number will be $N_f=N/2-S_{\mathrm{tot}}^z$.

Use the canonical anticommutation relations $\{c_n,c_m^{†}\}={\delta }_{nm}$ and $\{c_n,c_m\}=0$. In particular, $n_m=c_m^{†}{c}_m$ is a projector, so $(1-2n_m{)}^2=1$. Define $K_n={\prod }_{m<n}(1-2n_m)$. Then $K_n^{†}=K_n$, $K_n^2=1$, and $K_n$ commutes with $c_n$ and $c_n^{†}$. Taking the adjoint of the given raising operator gives $S_n^{-}=c_n^{†}{K}_n$. HenceAlso $[n_n,c_n]=-c_n$ and $[n_n,c_n^{†}]=c_n^{†}$ imply $[S_n^z,S_n^{\pm }]=\pm S_n^{\pm }$. For different sites, exchanging two fermion operators contributes one minus sign and commuting an operator through the other site's parity string contributes another; their product is positive. Thus distinct spin algebras commute. The string supplies precisely the signs required for the Jordan–Wigner transformation. The fermion vacuum corresponds to spin up, while one occupied site corresponds to spin down; this convention fixes subsequent signs.

For a bulk bond, $K_{n+1}=K_n(1-2n_n)$ and $c_n(1-2n_n)=-c_n$. Reordering the remaining different-site operators then givesThe periodic spin boundary condition requires separate attention. Set $P={\prod }_{n=1}^N(1-2n_n)=(-1{)}^{N_f}$. The end bond becomes $S_N^{+}{S}_1^{-}=-P{c}_1^{†}{c}_N$, with the conjugate relation for its adjoint. Consequently the exact finite-ring Hamiltonian isfermion parity is conserved, so diagonalize in a sector $P=\pm 1$. There the end bond can be written like every other bond provided $c_{N+1}=-P{c}_1$. This parity-dependent boundary condition of the Jordan-Wigner chain is antiperiodic for even fermion number and periodic for odd fermion number. Simply assuming periodic fermions in both sectors changes the finite-size spectrum.

In a fixed fermion-parity sector use lattice spacing one and the unitary discrete Fourier transformThe allowed momenta satisfy $e^{ikN}=-P$: they are $k=2\pi (\ell +1/2)/N$ for $P=+1$ and $k=2\pi \ell /N$ for $P=-1$, with $N$ distinct values modulo $2\pi$. The hopping and number sums become ${\sum }_{k}2\cos k{c}_k^{†}{c}_k$ and ${\sum }_k{c}_k^{†}{c}_k$. ThusHere $E$ is the energy of the all-up spin configuration, not generally the ground-state energy. The dispersion has minimum $h-J$ at $k=0$, maximum $h+J$ at $k=\pi$, and positive curvature at its minimum. Since the fermions do not interact, each mode can be occupied zero or one time. A negative-energy occupation lowers the many-body energy, so in the thermodynamic limitThis is a filled Fermi sea containing many fermions; it is not obtained by putting a single particle into the band minimum. For $h>J$ every mode has positive energy and the vacuum is the ground state. Its spin interpretation is complete longitudinal polarization. At $h=J$ the zero-momentum one-particle excitation first becomes gapless.

For an exact finite ring use parity-consistent filling of a free-fermion band. In each momentum grid fill negative-energy levels, check that their particle number has the chosen parity, and if necessary toggle a mode with the smallest $|{ϵ}_k|$. Compare the two sector energies. A zero-energy mode allows degenerate occupations. This restores the physical spin Hilbert space without changing the stated band.

For one coordinate with Hamiltonian $H={\hat{p}}^2/(2m)+V(\hat{x})$, the propagator is $K(x_f,t_f;x_i,t_i)=⟨x_f|e^{-iHT/\hslash }|x_i⟩$, where $T=t_f-t_i>0$. Divide $T$ into $M$ intervals $\delta t$, use the Trotter product formula to separate kinetic and potential evolution, and insert complete position and momentum bases at each slice. The momentum integral is Gaussian:Multiplying the short-time kernels gives the time-sliced configuration-space path integral,with $x_0=x_i$, $x_M=x_f$. Continuum notation abbreviates this limit as ${\int }_{x_i}^{x_f}\mathcal{D}x{e}^{iS[x]/\hslash }$, where $S=\int dt[m{\dot{x}}^2/2-V(x)]$. The measure includes the displayed normalization; it is not an unweighted count of paths. A small damping or analytic-continuation prescription defines the oscillatory integrals. The paths contribute complex amplitudes, so alternative histories interfere before probabilities are calculated. This is the central difference from a classical statistical mixture of trajectories.^[ref]

When actions vary rapidly relative to $\hslash$, stationary phase selects classical paths. Put $x=x_{\mathrm{cl}}+\eta$, with $\eta$ zero at the endpoints. Vanishing first variation gives $m{\ddot{x}}_{\mathrm{cl}}+V^{\prime }(x_{\mathrm{cl}})=0$. The quadratic term isThe resulting Gaussian fluctuation determinant is the semiclassical propagator prefactor. In one dimension, away from caustics,The branch convention makes $\nu =0$ for a direct free-particle path before any caustic. The Maslov index tracks the phase changes as fluctuations pass through a focal singularity. In several dimensions replace the endpoint derivative by the Van Vleck determinant. The approximation sums all relevant classical paths; real-time action need only be stationary, not a minimum. It fails in this simple form at a caustic or coalescing saddles, where a uniform approximation is needed. Quadratic actions have no higher fluctuation terms, making their Gaussian evaluation exact. For example $S_{\mathrm{cl}}=m(x_f-x_i{)}^2/(2T)$ gives the exact free-particle kernel $\sqrt{m/(2\pi i\hslash T)}{e}^{i{S}_{\mathrm{cl}}/\hslash }$. Tunnelling involves classically forbidden regions and can be described by imaginary-time saddles rather than real classical trajectories.

Quantum statistical mechanics follows by Wick rotation $t=-i\tau$. The imaginary-time evolution operator is $e^{-\tau H/\hslash }$, with Euclidean actionTaking a thermal trace identifies and integrates the endpoint positions:This thermal trace as a closed Euclidean path integral explains the periodic imaginary-time direction. A density-matrix element instead retains open endpoints. For many identical particles the trace also sums over endpoint permutations, with a plus sign for bosons and permutation parity for fermions. In coherent-state language this becomes periodic bosonic fields and antiperiodic Grassmann fields, with corresponding Matsubara frequencies.

The Euclidean weight resembles a classical Boltzmann weight. A time-sliced local quantum field in $d$ spatial dimensions becomes a classical-looking field model in $d+1$ dimensions, with anisotropic couplings and extra extent $\beta \hslash$. Thus zero-temperature quantum criticality can be treated using an infinitely extended imaginary-time direction, while nonzero temperature imposes a finite thickness. The dynamical scaling of time need not equal the spatial scaling: the dynamical critical exponent specifies their relationship. This quantum-to-classical imaginary-time mapping requires qualification when weights have signs or imaginary terms; such a functional integral is not automatically an ordinary positive classical probability distribution. Even with a positive weight, the mapping concerns equilibrium statistics, not literal dissipative real-time trajectories. For a free particle the Euclidean kernel is a diffusion kernel with diffusion coefficient $\hslash /(2m)$, illustrating the mathematical correspondence.

The classical high-temperature limit is a further approximation, rather than simply another name for imaginary time. If quantum level spacings and exchange effects are negligible, a one-particle partition function in $d$ dimensions becomesFor indistinguishable particles there is a factor $1/N!$ in the dilute classical gas, and exchange is negligible when the number density times the thermal de Broglie wavelength to power $d$ is small. For a harmonic oscillator, the exact quantum result $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}$ tends to $1/(\beta \hslash \omega )$ as $\beta \hslash \omega \to 0$, recovering the classical phase-space integral. stationary phase explains classical mechanics; the short imaginary-time, weak-exchange limit explains classical statistical mechanics. These are related but distinct limits of the same quantum formulation.

Consider a conventional spin-singlet, weak-coupling, three-dimensional superconductor. Set $\hslash =k_B=1$ in the microscopic action and restore electromagnetic units below. An attraction $\lambda >0$ in an energy shell about the Fermi surface gives a Grassmann actionThe attractive sign permits pairing even when the single-particle Fermi gas is otherwise stable. A complex Hubbard–Stratonovich transformation introduces a pair field $\Delta$ proportional at the saddle to $\lambda ⟨{\psi }_{\downarrow }{\psi }_{\uparrow }⟩$. The action is then quadratic in fermions. Integrating them out gives the BCS pairing-field effective actionHere the trace uses one independent Nambu spinor pairing block; a fully doubled particle-hole representation instead carries its corresponding one-half factor. At a uniform saddle the Bogoliubov--de Gennes Hamiltonian can be written asThe negative branch describes the paired vacuum; excitations are Bogoliubov quasiparticles mixing particles and holes. Extremizing the effective action givesThe physical cutoff or a renormalized interaction is part of the microscopic model; an unregulated contact integral must not be used indiscriminately. With constant density of states $N(0)$ per spin and shell half-width ${\Omega }_c$, the zero-temperature weak-coupling solution is ${\Delta }_0\simeq 2{\Omega }_c{e}^{-1/[N(0)\lambda ]}$. The exponentially small scale follows from the logarithmic pairing susceptibility of states near the Fermi surface. For this ideal weak-coupling model, $2{\Delta }_0/(k_B{T}_c)\simeq 3.53$. The microscopic theory relates the gap, pairing and thermodynamic transition.^[ref]

Close to a continuous transition, for fields varying slowly compared with microscopic scales, expand the effective action in powers and gradients of $\Delta$. Its static part is the Ginzburg-Landau theory of superconductivity:The normalization of $\Delta$ fixes the coefficients; here $a(T)$ changes sign at $T_c$. The uniform minimum has $|{\Delta }_0{|}^2=-a/b$ for $a<0$, condensation free-energy density $-a^2/(2b)$, and characteristic coherence-length scale $\xi =\sqrt{K/|a|}$. Above the transition the minimum is at zero. The quartic stabilization and gradient cost are precisely the ingredients of Landau theory for a complex order parameter. This controlled near-transition expansion supplies a bridge between microscopic BCS theory and macroscopic magnetic response; its coefficients and validity depend on the material and fluctuation regime.^[ref]

Write $\Delta =(|{\Delta }_0|+\eta )e^{i\varphi }$. The radial fluctuation has positive local restoring curvature, whereas the phase has no local potential in a neutral theory with global $U(1)$ symmetry. Freezing the amplitude givesThe neutral case supports a gapless collective phase mode. In the charged case the electromagnetic field couples to this phase rigidity. Under $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$ and $\varphi \mapsto \varphi +q\chi /\hslash$, the displayed combination is invariant. Locally choosing a gauge with constant phase leaves a quadratic term in $\mathbf{A}$. Varying the free energy gives the transverse London current and magnetic screening scale,Combining this relation with static Maxwell equations gives exponential field penetration, the Meissner effect. Single-valuedness of the pair field around a loop with negligible current gives flux in units ${\Phi }_0=2\pi \hslash /|q|=\pi \hslash /e$. These are gauge-invariant physical consequences of the condensate stiffness.

For the dynamics, the neutral phase and density mode is reorganized by the Coulomb interaction and the dynamical electromagnetic field. In a bulk three-dimensional system the long-wave longitudinal mode has a plasma gap, rather than remaining the neutral sound mode. This is the Anderson-Higgs mechanism in a superconductor: a phase coordinate becomes part of the electromagnetic collective response instead of surviving as an independent massless physical particle. Local gauge transformations are redundancies of the description, so a literal breaking of a local redundancy is not an observable. Phase coherence, current response and magnetic screening provide the physical statements.^[ref]

A static Landau functional alone cannot determine that full dynamical spectrum. A two-dimensional film, for example, need not have the same long-wave plasma dispersion as a bulk charged system. Nor does a positive radial curvature automatically guarantee a sharp Higgs amplitude resonance: coupling to quasiparticles and other modes can damp it. The chosen mean-field phase is also a thermodynamic symmetry-breaking description; an isolated finite system can be formulated with definite particle number. The microscopic pairing saddle, the near-transition Landau expansion and the electromagnetic Anderson-Higgs response describe successive scales of the same superconducting state, with distinct assumptions at each step.

For one bosonic mode, define the unnormalized bosonic coherent state $|z⟩=e^{z{a}^{†}}|0⟩$. The commutator $[a,a^{†}]=1$ gives $a|z⟩=z|z⟩$ and $⟨z^{\prime }|z⟩=e^{{\bar{z}}^{\prime }z}$. The normalized state is $e^{-|z{|}^2/2}|z⟩$; it is not orthogonal to a different coherent state. The coherent-state resolution of identity follows from Gaussian moments,For many modes take products of these states and measures. The complex amplitudes have the natural scale of square root occupation; coherent states turn normal-ordered matrix elements into classical symbols without replacing the quantum trace by a single classical configuration.

For fermions, $\{c_i,c_j^{†}\}={\delta }_{ij}$ and a nilpotent annihilation operator cannot have an arbitrary nonzero ordinary complex eigenvalue. Introduce independent Grassmann variables ${\eta }_i,{\bar{\eta }}_i$, anticommuting among themselves and with the fermionic operators, and define $|\eta ⟩=\exp (-{\sum }_i{\eta }_i{c}_i^{†})|0⟩$. Then $c_i|\eta ⟩={\eta }_i|\eta ⟩$ and $⟨{\bar{\eta }}^{\prime }|\eta ⟩=e^{{\sum }_i{\bar{\eta }}_i^{\prime }{\eta }_i}$. With an ordered Berezin integral measure satisfying $\int d\bar{\eta }d\eta \bar{\eta }\eta =-1$ for one pair,for an even operator $O$; this is the fermionic coherent-state trace. In one mode the trace integrand is $e^{-\bar{\eta }\eta }(O_{00}-\bar{\eta }\eta O_{11})$, whose integral is $O_{00}+O_{11}$. Without the sign reversal it would be the supertrace, not the thermal trace. These formal Grassmann-valued coherent states do not represent ordinary complex-valued fermion amplitudes.

To construct the many-body thermal path integral, write $\mathcal{Z}=\mathrm{Tr}{e}^{-\beta K}$ with $K=H-\mu \hat{N}$ and divide the imaginary-time interval $\mathcal{T}=\beta \hslash$ into slices of width $\delta \tau$. Insert the appropriate coherent-state resolution of identity between factors. For a normal-ordered Hamiltonian, the adjacent-slice kernel iscoherent-state time slicing makes the overlap and Gaussian measure combine into the first-order time derivative. The grand-canonical coherent-state action in the regulated continuum limit isFor bosons the integration variables are complex, with the conjugate labels on the original integration contour; for fermions the barred and unbarred Grassmann variables are independent. The continuum notation retains the time-sliced ordering prescription. For example, $n(n-1)=a^{†2}{a}^2$ has normal symbol ${\bar{z}}^2{z}^2$, whereas $n^2$ has that symbol plus $\bar{z}z$. Dropping that distinction changes a theory. The coherent-state thermal boundary conditions follow because the bosonic trace closes paths periodically and the fermionic sign in the trace closes them antiperiodically:The barred fields obey the same respective boundary conditions. Hence Matsubara frequencies are $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. The derivative term encodes the canonical structure; omitting it would remove the quantum dynamics.^[ref]

As a check, a free diagonal mode with $K=\xi a^{†}a$ gives $Z_B=(1-e^{-\beta \xi }{)}^{-1}$ for $\xi >0$, whereas $K=\xi c^{†}c$ gives $Z_F=1+e^{-\beta \xi }$. Gaussian integration produces an inverse determinant for bosons and a determinant for independent complex fermions, with their respective temporal boundary conditions. Differentiating $\log Z$ recovers Bose-Einstein and Fermi-Dirac occupations. Interactions become higher powers in the normal-ordered symbols, and sources generate time-ordered correlations. The coherent-state construction preserves statistics, ordering and boundary conditions together; a formal continuum expression without its measure and temporal prescription is incomplete.

For four-momentum $k$ flowing along the fermion flow, the momentum space Feynman rules areThe scalar propagator and Dirac propagator invert their respective quadratic kernels. The interaction vertex follows by expanding $\exp (i\int {\mathcal{L}}_{\mathrm{int}})$, fixing its sign. Species is conserved along each open fermion line; there is no vertex changing species 1 into species 2.

An interaction vertex carries $(2\pi {)}^4{\delta }^{(4)}(\sum k)$ with inward four-momenta. Impose momentum conservation and integrate each independent loop momentum with $d^{4}k/(2\pi {)}^4$. Include the closed fermion loop sign $-1$ for every fermion loop, the relevant fermionic signs when rearranging external states, and the Feynman-diagram symmetry factor. External incoming and outgoing fermions carry $u_s(p)$ and ${\bar{u}}_s(p)$; incoming and outgoing antifermions carry ${\bar{v}}_s(p)$ and $v_s(p)$. An external real scalar field leg has factor 1 in an amputated scattering amplitude. The Feynman i-epsilon prescription specifies the poles.

Let $p_1,p_3$ label the incoming and outgoing fermion, and $p_2,p_4$ the incoming and outgoing antifermion. The bottom row of the figure gives the two tree-level Feynman diagrams:

Thus the Mandelstam variables in the two denominators areThe two scattering amplitude terms have a relative fermionic sign of minus when written with spinor products $({\bar{u}}_3{u}_1)({\bar{v}}_2{v}_4)$ and $({\bar{v}}_2{u}_1)({\bar{u}}_3{v}_4)$. The overall phase depends on the chosen ordering of the external states. No $u$ channel joins an incoming fermion to an outgoing antifermion through a species-preserving Yukawa interaction.

With all four external fermions of species 1, both final attachments are possible. The middle row of the figure gives the two tree-level Feynman diagrams:

For a fixed ordering of external fermion states, the fermionic sign from exchanging the final identical fermions givesHere $u_j=u_{s_j}(p_j)$ are Dirac spinors. In particular, the scattering amplitude changes sign when the complete final labels $3$ and $4$ are exchanged, as required for identical fermions.

Label the initial species-1 and species-2 fermions by $(p_1,s_1)$ and $(p_2,s_2)$, and the final species-1 and species-2 fermions by $(p_3,s_3)$ and $(p_4,s_4)$. There is exactly one tree-level Feynman diagram: the $t$-channel scalar propagator joins the continuous species-1 line $1\to 3$ to the continuous species-2 line $2\to 4$. Its exchanged four-momentum $q_t$ determines $t$, in the notation of Mandelstam variables:The top row of the figure displays this Feynman diagram. A crossed attachment $1\to 4$ would change species at an interaction vertex and is forbidden. There is no initial antifermion for an annihilation channel.

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$.

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.

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.