The marginal likelihood integrates the likelihood function against the prior density:It is the normalizing constant in Bayes' theorem. For the subsequent posterior distribution and importance sampling calculations, assume a normalized prior density and $0<Z<\infty$.

Let $m$ count the parallel neighboring pairs, or frustrated bonds in an antiferromagnetic Ising ring. The periodic boundary conditions implyThere are $N-m$ antiparallel pairs, so $(-1{)}^{N-m}=1$: the number of antiparallel pairs must be even. For odd $N$, $m$ must consequently be odd. Each frustrated bond in an antiferromagnetic Ising ring raises the energy by $2J$ relative to a satisfied pair, givingThe minimum is attained at $m=1$. Choose its position among the $N$ labeled bonds, choose either orientation of one Ising spin, then alternate across every other bond. This constructs exactly two spin configurations per position. ThusThis is geometric frustration: every local interaction prefers opposite Ising spins, but the closed odd ring prevents all preferences being satisfied at once. No particular location is preferred for the single frustrated bond in an antiferromagnetic Ising ring.

Each neighboring pair of Ising spins contributes at least $-J$ to the energy, with equality exactly when the two Ising spins are opposite. For even $N$, alternation closes consistently under the periodic boundary conditions. The ground states are thereforeHere $g_0$ is the ground-state degeneracy. Choosing $S_1$ fixes every other Ising spin, so these are all the ground states; translations do not give additional spin configurations. The two ground states are exchanged by spin inversion symmetry.

Take $L\ge 2$, so neighbouring operators act at distinct sites. The identityexpresses the Hamiltonian operator of this quantum XY chain in spin ladder operators. Substituting the Schwinger boson representation gives the requested quartic expression:Operators at different sites commute, allowing normal ordering without a correction here. For the first quartic term, the inverse discrete Fourier transforms giveHere ${\delta }^{(L)}$ is one when its two momenta agree modulo $2\pi$ and zero otherwise; it expresses lattice momentum conservation. The other quartic term has the same ordered operators with phase $e^{i(p-r)}$. On the support of this Kronecker delta, $p-r=t-q$ modulo $2\pi$, so the two phases sum to $2\cos (q-t)$. Thus the Fourier vertex of a Schwinger boson chain isAn equivalent three-sum expression sets $t=p+q-r$ modulo $2\pi$, so its coefficient is $\cos (r-p)$. The factor $1/L$ comes from four factors $L^{-1/2}$ and one site sum, while the factor of two from the two exchange terms cancels the original $1/2$. The formula is Hermitian: conjugation exchanges $p$ with $t$ and $q$ with $r$, and the coefficient is unchanged after imposing lattice momentum conservation. It also preserves each local total boson number in its real-space form, hence preserves the physical constrained subspace.

For the specified positive-exponent discrete Fourier transform, the inverse discrete Fourier transforms have negative exponents:Every sum runs over $q=2\pi j/L$, $j=0,\dots ,L-1$. Indeed, orthogonality of roots of unity givesso substituting the forward discrete Fourier transform into either proposed inverse recovers the original operator. Taking Hermitian conjugates supplies the inverse formulas for the creation operators.

Multiplying the inverse discrete Fourier transforms gives the local number operator combinationThis local expression includes off-diagonal momentum terms. Only after summing over sites does orthogonality of roots of unity giveIn a physical fixed-$s$ state the first boxed operator acts as $s$ at every site, by the local spin constraint in the Schwinger boson representation.

Let $w$ count the unequal neighbouring pairs, or domain walls. The energy is $E=-JL+Jw$. A single domain wall is impossible with periodic boundary conditions: following all the other, equal-state bonds around the ring would force its endpoints to be equal too. Two domain walls are possible, splitting the ring into two nonempty uniform arcs carrying different states. Thus the first excited state hasChoose the two wall bonds in $\left(\genfrac{}{}{0ex}{}{L}{2}\right)$ ways. For that choice, identify one of the resulting arcs by its smallest site label. Its state has $Q$ choices and the other arc has $Q-1$. This counts each configuration exactly once, giving the two-domain-wall degeneracy of a Potts ringFor example, $Q=2,L=4$ gives $d_0=2$ and $d_1=12$; the remaining two configurations are alternating and have four domain walls. The formula also covers $L=2$, where the specified ring energy counts the two directed neighbouring bonds, and an unequal pair breaks both. If $Q=1$ or $L=1$, every allowed configuration is a ground state and there is no first excited state. For $Q\ge 3$, the number of domain walls need not be even: three distinct states on a three-site ring give three, although one remains impossible.

Assume $Q\ge 2$ and $L\ge 2$; the trivial exceptions are described below. In this periodic Potts chain, every satisfied bond contributes $-J$ to the energy, and no bond can contribute less. All $L$ bonds can be satisfied simultaneously by choosing the same state at every site. Conversely, satisfying every bond forces this uniform configuration. Therefore the ground state energy and degeneracy factor areThe $Q$ choices are distinct configurations on the labelled lattice: a global relabelling of the states is a symmetry, not an identification of microstates.

For commuting real scalar fields, remove the closed fermion loop sign and replace the Dirac propagators by scalar propagators. There are no spinor indices, hence no matrix trace. With two independent real scalar fields ${\chi }_1,{\chi }_2$ and interaction $g\varphi {\chi }_1{\chi }_2$, the same normalization givesThe overall sign has reversed relative to part (i) because $(ig{)}^2$ is now unaccompanied by a closed fermion loop sign. The same formula holds in the source's kernel convention, since there are again four propagators.

There is an additional convention to specify if the two symbols are intended to denote the same real scalar field. For the literal interaction $g\varphi {\chi }^2$, there are two identical connected Wick contractions between ${\chi }^2(x)$ and ${\chi }^2(y)$, so the coefficient becomes $-2g^2$. If the conventional interaction is instead $(g/2)\varphi {\chi }^2$, each vertex is $ig$ and the loop has symmetry factor $1/2$, giving coefficient $-g^2/2$. In both interpretations the closed fermion loop sign and spinor matrix trace disappear; the differing numerical factor comes from the interaction normalization and field identity.

Normalize the generating functional by its value with all sources zero, and take the connected part of the time-ordered correlation function. Source differentiation of the interaction exponential is equivalent to the Dyson series followed by Wick's theorem. Write $D(x-y)=⟨T\varphi (x)\varphi (y){⟩}_0$ and $F_{ij}(x-y)=⟨T{\psi }_i(x){\bar{\psi }}_j(y){⟩}_0$ temporarily, so that every contraction has an unambiguous meaning.

The second-order Dyson series contributesThere are two connected external scalar field attachments:They become equal after exchanging the integration variables $x,y$, cancelling $2!$. The connected fermion contractions have the closed fermion loop sign:The minus sign can also be checked directly: the pairing of the first with the fourth field and the second with the third has positive pairing parity, but $⟨T{\bar{\psi }}_i(x){\psi }_j(y)⟩=-F_{ji}(y-x)$. Multiplying it by $(ig{)}^2=-g^2$ leaves a positive result:A contraction of the two external fields with each other times a vacuum bubble cancels against the normalization. Separate fermion tadpole integrals at $x$ and $y$ belong to a product of one-point functions, removed by taking the connected part.

To match the source's generating functional convention, distinguish its displayed kernels from contractions. The positive $iJ\Delta J/2$ exponent gives $D=-i\Delta$. With the source insertion order in equation (5), the positive $i\bar{\eta }{S}_F\eta$ exponent similarly gives $F=-i{S}_F$. The two scalar propagators and two Dirac propagators supply a total phase $(-i{)}^4=1$. Thus, writing $S=S_F$ as in the requested formula,The matrix trace closes the spinor indices; no extra factor of two or additional symmetry factor remains for a Dirac field with this Yukawa interaction.

This cannot be the sole symmetry-breaking sector of the Standard Model. Its gauge-boson mass matrix has two zero eigenvalues, so besides the electromagnetic photon it predicts an extra massless neutral gauge boson. Every linear combination of $W^3$ and $B$, including the combination normally called $Z$, is massless. The observed massive Z boson is a concrete contradiction. Choosing $g,g^{\prime },m,\lambda$ differently cannot change these two zero eigenvalues.

A real scalar triplet can occur as an additional sector, but a further vacuum expectation value, for example that of a Higgs doublet, is then necessary to obtain the observed neutral gauge boson mass pattern. This does not rescue the triplet-only model posed here.

Let $g$ and $g^{\prime }$ be the two gauge couplings. In unitary gauge, $\Phi =(0,0,v+h{)}^T$ and the weak hypercharge term in the gauge covariant derivative vanishes. Using $(T^a{)}_{bc}=-i{ϵ}_{abc}$, the kinetic term becomesThus, in the basis $(W^1,W^2,W^3,B)$, the gauge-boson mass matrix and its eigenvalues areHere $W^{\pm }=(W^1\mp i{W}^2)/\sqrt{2}$. The two massless neutral gauge bosons may be rotated to any orthonormal basis; no neutral mass term selects a unique mixing angle.

Expanding the same kinetic term to first order in $h$ givesFunctional differentiation gives the Feynman rulesand, in the real basis, $h{W}_{\mu }^a{W}_{\nu }^b:2i{g}^{2}v{\delta }^{ab}{g}_{\mu \nu }$ for $a,b\in \{1,2\}$. All such three-point Feynman vertices involving $W^3$ or $B$ vanish. The factor two in the real basis comes from the two identical gauge boson fields; in the charged basis it is already explicit in ${\mathcal{L}}_{hVV}$.

Put $r^2={\Phi }^T\Phi$. The potential has its minimum atChoose the vacuum expectation value ${\Phi }_0=(0,0,v{)}^T$. Since the weak hypercharge is zero, $Q=T^3$ on this real scalar triplet, and $T^3{\Phi }_0=0$. Neither $T^1$ nor $T^2$ annihilates ${\Phi }_0$. The stabilizer subgroup is consequentlyIn particular, spontaneous symmetry breaking has not reduced the two neutral gauge generators to a single one: $Y$ annihilates every component and remains unbroken independently of $T^3$.

The scalar mass matrix isWith canonical kinetic normalization $\frac{1}{2}(\partial h{)}^2$, the radial Higgs boson therefore hasThe two tangential zero eigenvalues are Goldstone modes, absorbed by the two massive gauge bosons through the Higgs mechanism.

There are two tree-level Feynman diagrams, with intermediate photon and gluon respectively. Their color charge tensors are ${\delta }_{ji}{\delta }_{kl}$ and ${\sum }_A{T}_{ji}^A{T}_{kl}^A$. The pure photon factor isTheir interference vanishes after the color charge sum because of photon-gluon colour orthogonality:Thus there is no term proportional to $e^2{g}_s^2$, and with the $F_u$ of part (iii),Here $Q_u=Q_t=2/3$. The absence of interference concerns the specified averaged and summed result; a fixed choice of external color charge need not have zero interference.

The exchanged particle is now a gluon, and the incoming up quark mass replaces $m_e$. DefineThe Dirac trace and spin average are unchanged. Normalize the fundamental gauge generators by $\mathrm{tr}(T^A{T}^B)={\delta }^{AB}/2$. Averaging over the $N_c^2$ initial color charge choices and summing the final color charge choices givesReplace the factor $N_c{e}^4{Q}_e^2{Q}_t^2$ from part (ii) by $g_s^4{C}_g$. ThereforeKeeping $m_u$ is allowed; if the usual negligible up quark mass approximation is wanted, simply set $m_u=0$. The factor $1/N_c^2$ is essential: initial quark color charge states are averaged, whereas final ones are summed.

Let $\Sigma =m_e^2+m_t^2$, $Q_e=-1$, $Q_t=2/3$ and $N_c=3$. The single tree-level Feynman diagram is annihilation through a virtual photon. Up to an irrelevant common phase, its scattering amplitude isThere are two initial spin angular momentum states for each particle, so the spin average is $1/4$. The final color charge sum gives ${\sum }_{c,d}|{\delta }_{cd}{|}^2=N_c$. There is no additional photon polarization average for a virtual photon and no separate polarization label beyond the Dirac spinor spin sums.

The spin sums and the gamma-matrix trace identities giveHere $p_1\cdot p_2+m_e^2=p_3\cdot p_4+m_t^2=s/2$. Contracting givesSince $2p_1\cdot p_3=2p_2\cdot p_4=\Sigma -t$ and $2p_1\cdot p_4=2p_2\cdot p_3=\Sigma -u$, this is the massive vector-current annihilation kernelThus the absolutely normalized result isAfter taking out the displayed prefactor, $\boxed{A=1,B=4,C=-2}$. Because the source uses a proportionality sign, only their ratio is fixed until that prefactor is specified. Equivalently, after taking out just $e^4/s^2$, their values are $8/3,32/3,-16/3$. No mass has been neglected. At $m_e=0$ and $s=4m_t^2$, the bracket yields the finite isotropic threshold value of the scattering amplitude; the vanishing threshold scattering cross-section comes from relativistic two-body phase space.

Use the Mandelstam variables $s=(p_1+p_2{)}^2$, $t=(p_1-p_3{)}^2$ and $u=(p_1-p_4{)}^2$. The external on shell relations are $p_1^2=p_2^2=m_e^2$ and $p_3^2=p_4^2=m_t^2$. Expanding the three Mandelstam variables and using four-momentum conservation, $p_3+p_4=p_1+p_2$, gives

There are four tree-level Feynman diagrams: three photon exchange graphs with pairings $12|34$, $13|24$, $14|23$, and the four-photon contact interaction. They are drawn in the bottom row of the figure. With no external fermions or Faddeev-Popov ghost fields, an internal matter or ghost line would have to close into a loop; it cannot add a tree-level Feynman diagram here.

Use all-incoming null four-momenta $k_i$ with ${\sum }_i{k}_i=0$. Let ${\epsilon }_i$ be the external photon polarization vectors, including the complex conjugation for physical outgoing photons, and impose $k_i\cdot {\epsilon }_i=0$. For one pairing $ij|kl$, put $K=k_i+k_j=-(k_k+k_l)$. At the first cubic Feynman vertex, the terms proportional to $k_i\cdot {\epsilon }_i$ and $k_j\cdot {\epsilon }_j$ vanish, leavingAt the other end, the incoming internal four-momentum is $+K$, giving $+2\lambda ({\epsilon }_k\cdot {\epsilon }_l)K^{\sigma }$. The Feynman gauge photon propagator then givesThe last expression is evaluated at generic kinematics, with singular channel limits taken after cancellation. Summing the three pairings produces $+4i{\lambda }^{2}P$, whereThe quartic Feynman vertex contributes $i{\mathcal{M}}_{\mathrm{contact}}=-4i{\lambda }^{2}P$. ThusThis is the tree photon-scattering cancellation in a quadratic Abelian gauge: the individual gauge fixing interactions cancel in the physical scattering amplitude. The argument uses transverse external photon polarization vectors; it does not assert that an arbitrary off shell four-point correlation function vanishes. Physical light-by-light scattering first occurs through charged-particle loop Feynman diagrams.

Expand the gauge fixing square:The second term is a derivative cubic photon interaction; the third is a quartic photon contact interaction. Both arise from nonlinear Abelian gauge fixing, rather than an additional physical electromagnetic force.

For $\xi =1$, differentiating the quartic term with respect to four identical gauge fields gives eight copies for each pairing, henceWith the Fourier convention $A(x)=\int d^{4}p(2\pi {)}^{-4}{e}^{ipx}A(p)$ and all four-momenta incoming, the differentiated cubic term givesmatching the supplied Feynman rules. Reversing the Fourier convention reverses this cubic rule's sign; each exchange graph contains two such Feynman vertices, so its contribution is unchanged.

For fixed $A$, the delta functional sets $\sigma (x)=F[A](x)$. Therefore its weighted functional integral isSince the action enters as $e^{iS}$, the gauge fixing contribution isThe Faddeev-Popov determinant does not depend on $\sigma$ and remains in the path integral. The limit $\xi \to 0$ enforces the original condition sharply; for nonzero $\xi$ this is a weighted average of the gauge slices.

An ordinary bosonic Gaussian functional integral produces an inverse determinant. An integral over independent Grassmann fields $c,\bar{c}$ instead produces a determinant, allowing a local action to represent the Faddeev-Popov determinant. With one convenient antighost sign convention,Multiplying the operator by the field-independent factor $-i$ changes only the normalization of its determinant. After integration by parts, an equivalent density isThus the Faddeev-Popov ghost fields have a massless kinetic term and a photon-ghost-antighost interaction proportional to $\lambda$. They are auxiliary anticommuting scalar fields, not physical external particles. Closed ghost loops carry the minus sign of Grassmann variables. Replacing $\bar{c}$ by $-\bar{c}$ reverses both displayed signs without changing any observable; both the kinetic and interaction conventions must be changed together.