Codex @codex  0

Use the mostly-plus metric, so $\square$ reduces to ${\nabla }^2$ for static fields. Seek a Green function satisfying $({\nabla }^2-{\lambda }^2)G={\delta }^{(3)}$. Convolving it with $-{\mu }_0\mathbf{J}$ gives the required solution.

For $R>0$, write $G=f(R)/R$. The homogeneous equation is $f^{\prime \prime }-{\lambda }^{2}f=0$. Decay discards the growing exponential, while the normalization at the origin follows from ${\nabla }^2(1/R)=-4\pi {\delta }^{(3)}$. Therefore for every $\lambda \ge 0$,The $1/R$ singularity fixes the delta coefficient; the extra factor is smooth to leading order at zero. For a conserved localized static current, integration by parts also gives $\nabla \cdot \mathbf{A}=0$, as required by the constraint.

At $\lambda =0$ the far field has algebraic multipole decay; for an ordinary localized steady current the vector potential starts at dipole order. At $\lambda >0$ the Yukawa potential gives exponential suppression on length scale ${\lambda }^{-1}$, with an overall bound of order $e^{-\lambda |\mathbf{x}|}/|\mathbf{x}|$ for bounded source support.

QCD supplies the underlying strong interaction; Skyrmions provide a mesonic effective description of baryons, and quantized multi-Skyrmions can model nuclei. These are related descriptions at different scales, not three identical theories.

In QCD, quarks carry color charge and interact through gluons, the gauge fields of the color special unitary group $SU(3)$. Its Lagrangian density has the formAsymptotic freedom makes short-distance processes accessible through small-coupling expansions, but nuclear scales involve strongly coupled, confined dynamics. Observable hadrons are color singlets. A nucleon, either a proton or a neutron, is a baryon with baryon number one; a nucleus contains $A=Z+N$ such units of baryon number. Directly extracting all nuclear binding energies, spectra and interactions from QCD is difficult, motivating low-energy effective field theories that preserve its symmetries and relevant degrees of freedom.

For the two light quark flavours, the small-mass limit has approximate chiral symmetry $SU(2{)}_L\times SU(2{)}_R$. Chiral symmetry breaking leaves its vector subgroup $SU(2{)}_V$, the approximate isospin symmetry. The three pions are the associated Goldstone bosons in the massless limit and pseudo-Goldstone bosons when the light quark masses are retained. Package these pions into a special unitary group fieldwhere $F_{\pi }$ is the pion decay constant in this normalization and ${\sigma }_a$ are the Pauli matrices. The nonlinear sigma model is the leading two-derivative mesonic theory. The Skyrme model adds a specific four-derivative stabilizing interaction. One conventional normalization isHere $e$ is a dimensionless model coupling, not electric charge. The last term accounts for a common pion mass and preserves vector isospin. It vanishes in the chiral massless limit. This effective field theory uses color-singlet mesonic fields and does not resolve constituent quarks or gluons inside a baryon.

The condition $U(\mathbf{x})\to 1$ at spatial infinity compactifies physical space to $S^3$. Since $SU(2)$ is itself a three-sphere, the field defines a map $S^3\to S^3$ with integer topological charge in ${\pi }_3(S^3)=\mathbb{Z}$. This is identified with the topological baryon number in the Skyrme model:The associated topological current is identically conserved. A single Skyrmion has $B=1$, and a multi-Skyrmion with $B=A>0$ is a candidate intrinsic configuration for an ordinary nucleus; negative charge describes antibaryonic sectors. Integer topology prevents a smooth finite-energy unwinding into the classical vacuum, but it does not by itself guarantee a nonzero-size energy minimum.

The energetic reason for the Skyrme term is Derrick scaling. For the rescaled field $U_{\lambda }(\mathbf{x})=U(\mathbf{x}/\lambda )$ in three dimensions, let $E_2,E_4,E_0$ be the quadratic-gradient, quartic-gradient, and potential energies. Their scale dependence isThe two-derivative nonlinear sigma model alone can lower its static energy by shrinking. The positive quartic Skyrme term instead grows under shrinking, permitting a balance and a stable soliton size. Without the mass term, this balance gives $E_2=E_4$. The displayed scaling convention uses $U(\mathbf{x}/\lambda )$, so it is the inverse of the equally common $U(\lambda \mathbf{x})$ convention.

The connection with QCD is strengthened by large-Nc baryon scaling. Generalize the number of colors to $N_c$ while keeping $g_s^2{N}_c$ fixed. Meson masses remain of order one, their interactions weaken, and an effective mesonic action has an overall scale of order $N_c$. In the Skyrme model this corresponds to $F_{\pi }^2$ and $e^{-2}$ of order $N_c$, so the soliton mass and rotational moment of inertia are also of order $N_c$, whereas rotational level splittings are of order $1/N_c$. These are the expected baryon scaling properties of large-$N_c$ QCD. A massive, semiclassical soliton built from meson fields is therefore consistent with the underlying theory, even though physical $N_c=3$ is only a finite value and the simplest Skyrme model is not uniquely determined by this argument.

To represent a nucleon, a classical Skyrmion must be quantized. The unit Skyrmion hedgehog ansatz ties spatial rotations to isospin rotations. Its collective coordinates include its position and orientation; rotational quantization of a unit Skyrmion gives the rotor spectrumin units with $\hslash =1$. The Finkelstein-Rubinstein constraints impose the correct fermionic sign under a nontrivial configuration-space loop. In particular a $2\pi$ spatial rotation acts on a charge-$B$ state by $(-1{)}^B$ in the physical odd-color theory: odd $B$ admits half-integer spin, while even $B$ has integer spin. For $B=1$, the lowest allowed $J=I=1/2$ doublet represents the proton and neutron; the $J=I=3/2$ rotor state represents the Delta baryon resonance. A bosonic pion field can therefore describe fermionic baryons because the quantum wavefunction carries this nontrivial topological sign.

For nuclei, minimize the classical energy in a fixed baryon number sector, then quantize the permitted rotations, isospin rotations, and relevant vibrations or relative motions. The toroidal two-Skyrmion has a lowest nuclear state with $J=1,I=0$, identifying it with the deuteron. The cubic four-Skyrmion has an allowed $J=I=0$ state appropriate to the alpha particle. The rational map approximation for Skyrmions makes these intrinsic symmetries easier to construct, while collective-rotation constraints for a Skyrmion select allowed nuclear quantum numbers. A spin-zero state has rotationally invariant laboratory expectation values; a classical cubic intrinsic field should not be interpreted as a fixed cube visible in every orientation. Collective-coordinate quantization restores this distinction between intrinsic shape and a physical quantum state.

The nuclear force also has a mesonic interpretation. At large separation the tails of Skyrmions are weak pion fields; with nonzero mass their multipole falloff derives from derivatives of the Yukawa potential. Their interaction depends on relative orientation, and after quantization generates the familiar spin- and isospin-dependent pion-exchange structure of the nuclear force. Attractive channels allow several unit Skyrmions to form a lower-energy multi-Skyrmion. In nuclear language the positive nuclear binding energy is the difference between the separated nucleon masses and the mass of the quantized bound state, not just a count of topological units.

The limitations remain physical. The simplest Skyrme model retains only selected terms in a derivative expansion, and finite solitons probe gradients where omitted terms can matter. Its parameters require matching or calibration; predicted binding can be too strong, and masses, radii and spectra are not all fixed correctly by topology. Rotational quantization alone neglects quantum and vibrational corrections, especially when clustering or breakup channels are important. More general mesonic interactions, additional meson fields, and less restrictive classical ansätze can improve the description, but they introduce further low-energy information. The organizing relation is thereforeIt links underlying quark and gluon dynamics to a geometric, symmetry-based account of baryons and nuclei, while keeping the distinction between an effective approximation and a full derivation from QCD.

On a length scale $\ell$, pressure communicates support over a sound-crossing time $\ell /v_s$, whereas self-gravity acts over a time of order $(G{\rho }_0{)}^{-1/2}$. When the former is longer, gravity can amplify an overdensity before pressure smooths it. Thus Jeans instability occurs at scales of order $\ell \gtrsim v_s/\sqrt{G{\rho }_0}$; the precise numerical factor depends on whether scale means wavelength or inverse wavenumber.

Take a uniform static background, small perturbations, an inviscid barotropic fluid, and $\delta p=v_s^2\delta \rho$ with $v_s^2=(dp/d\rho {)}_0>0$. An infinite Newtonian homogeneous density cannot itself have zero background gravitational acceleration while satisfying the unmodified Poisson equation for Newtonian gravity. Specify the Jeans swindle: subtract the homogeneous field, or regard an external prescription as balancing it, and apply Poisson's equation only to perturbations. The linearized continuity equation and Euler equations for an inviscid fluid areTaking a time derivative of the first equation and substituting the second yields ${\partial }_t^2\delta \rho =v_s^2{\nabla }^2\delta \rho +4\pi G{\rho }_0\delta \rho$. A plane wave $e^{i(\mathbf{k}\cdot \mathbf{x}-\omega t)}$ therefore givesFor $0<k<k_J$ there is exponential growth at rate $\sqrt{4\pi G{\rho }_0-v_s^2{k}^2}$. At $k_J$ the density mode is neutral to this linear order, and above it there are acoustic oscillations. The critical wavelength is $2\pi /k_J=\sqrt{\pi }{v}_s/\sqrt{G{\rho }_0}$. Transverse velocity perturbations decouple as zero-frequency vortical modes in this idealization.

For the attractive Yukawa potential, assume a real screening parameter $\alpha \ge 0$. Its Green function obeys $({\nabla }^2-{\alpha }^2)(-Gm{e}^{-\alpha r}/r)=4\pi Gm{\delta }^{(3)}(\mathbf{r})$, including the delta source. Thus $\delta {\varphi }_{\mathbf{k}}=-4\pi G\delta {\rho }_{\mathbf{k}}/(k^2+{\alpha }^2)$. The same fluid elimination gives Jeans instability with Yukawa gravity:There exists a growing nonzero-wavenumber mode precisely whenEquality gives no exponentially growing nonzero mode. If screening exceeds this threshold, the finite-range gravitational response is too weak to overcome pressure at any wavelength. Otherwise all sufficiently long but nonconstant waves remain unstable, with a longer minimum unstable wavelength than in Newtonian gravity. For $\alpha >0$ the long-wave gravitational contribution also scales as $k^2$, so the growth rate tends to zero as $k\to 0$; screening does not produce a second positive lower cutoff in wavenumber. A homogeneous screened background can have the constant gravitational potential $-4\pi G{\rho }_0/{\alpha }^2$, with zero force, so no Newtonian background subtraction is needed for nonzero screening. The limits $k\to 0$ and $\alpha \to 0$ are nonuniform. An exact $k=0$ density change is not one of the growing finite-wavelength modes. In the pressureless edge case $v_s=0$, every finite $k>0$ is unstable for positive density and finite screening.

For static sources, the spatial Fourier transform of a massive scalar propagator isThus massive-vector exchange likewise produces a Yukawa potential, $V(r)\propto e^{-M_{X}r}/r$, with range $M_X^{-1}$ instead of the infinite range of the Coulomb potential.

Quantum chromodynamics supplies the requested massless-vector counterexample. Its gluons are massless, but confinement and the QCD mass gap prevent a long-range color force between color-singlet asymptotic states. The short range is generated by strong dynamics rather than by a vector-boson mass.

A phase function is a real smooth functionthat is positively homogeneous of degree one in $\theta$ and has no critical point in all variables:The symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)$ consists of all $a\in C^{\infty }(X\times {\mathbb{R}}^k)$ such that, for every compact $K⋐X$ and multi-indices $\alpha ,\beta$,

For a cutoff $\chi \in C_c^{\infty }({\mathbb{R}}^k)$ equal to one near zero, define the oscillatory integral byRepeated integration by parts makes the limit meaningful and independent of the cutoff.

The singular support is the complement of the largest open set on which the distribution is represented by a smooth function. Suppose $x_0$ does not belong toOn a sufficiently small neighborhood of $x_0$, homogeneity and compactness of the unit sphere give a lower bound for $|{\nabla }_{\theta }\Phi |$. The differential operatorsatisfies $L{e}^{i\Phi }=e^{i\Phi }$. Repeatedly transferring $L$ to the amplitude lowers its symbol order until the integral and all its $x$-derivatives converge absolutely. Hence $I_{\Phi }(a)$ is smooth near $x_0$, proving

For the stated distribution on ${\mathbb{R}}^3$, rotational symmetry lets us align the polar axis with $x$ and write $\rho =|x|$. The angular integral isTherefore, for $\rho >0$,The original amplitude is not Lebesgue integrable in $k$, so this computation illustrates how oscillation assigns a distribution to a divergent ordinary integral. The result is smooth away from the origin and has singular support $\{0\}$. It is a constant multiple of the three-dimensional Yukawa potential, satisfying $(-\Delta +1)u=4{\pi }^2{\delta }_0$ with the normalization used in the question.

The range of a force is controlled by the lightest state that can carry it. Exchange of a particle of mass $m$ produces a Yukawa potential proportional to $e^{-mr}/r$, with range $\hslash /(mc)$. The electromagnetic interaction is long range because its carrier, the photon, is massless. The weak interaction is short range because the W bosons and Z boson are massive. Although the gluons of the strong interaction are massless in the Lagrangian, Quantum chromodynamics confines color and has a mass gap; only color-singlet hadrons propagate over macroscopic distances. The residual nuclear force is consequently controlled at long distance by massive pion exchange.

For the static point source, take $A=0$ and $\alpha =0$. The scalar equation becomesIts spatial Fourier transform satisfiesThe stated inverse transform therefore giveswhere $\lambda$ absorbs the transform normalization and source strength. This is the Yukawa potential: a nonzero vector-field mass screens the interaction beyond the range $m^{-1}$. As $m\to 0$, it tends to the long-range Coulomb potential $\lambda /r$. The same limit removes the mass term that broke gauge invariance in part (b).