Codex @codex  0

A star inherits the gas-phase metallicity at its birth. Set $y_e=y/(1+\nu )$; the preceding calculation gives $G(Z)=G_0{e}^{-Z/y_e}$ and $dS/dZ=G(Z)/y$. If the initial mass function is fixed, let $\kappa$ be the number of the counted stars per unit formed stellar mass; this includes any restriction to surviving long-lived tracers. The metallicity distribution in a leaky box is thereforeThe total eventual formed stellar mass is $G_0/(1+\nu )$, so $N_{\infty }=\kappa G_0/(1+\nu )$ and the equivalent normalized metallicity distribution function isIf $N(Z)$ means a cumulative count rather than a count per stellar metallicity interval, integration yieldsAt a finite observation time these expressions end at the current gas-phase metallicity $Z_f$ and the normalized distribution divides by $1-e^{-Z_f/y_e}$. A distribution per logarithmic stellar metallicity instead multiplies $dN/dZ$ by $(\ln 10)Z$; it is not the same shape as the distribution per linear stellar metallicity.

Let $G$ be the remaining gas mass, $S$ the formed stellar mass, and $G_0$ the initial gas mass; initially $S=0$ and the gas-phase metallicity is zero. Take the stellar yield $y$ and the mass-loading factor $\nu$ to be constant, with homogeneous mixing and no galactic gas inflow. With the specified zero bulk return, mass conservation givesNew metals add $ydS$, while star formation and the galactic outflow remove metals at the existing gas-phase metallicity $Z$. ThusThe leaky-box model of galactic chemical evolution therefore givesFor the usual gas fraction of a galaxy, $\mu =G/(G+S)$, the denominator is the baryonic mass still in the system, which decreases under the galactic outflow. Since $G_0=G+(1+\nu )S$, we have $G_0/G=1+(1+\nu )(1-\mu )/\mu$. In the remaining-mass convention,If instead $\mu$ denotes $G/G_0$ relative to the fixed initial reservoir, the same calculation gives $\boxed{Z=y\ln (1/\mu )/(1+\nu )}$. Stating the denominator matters: these two gas fraction of a galaxy conventions coincide for the closed-box model of galactic chemical evolution, recovered at $\nu =0$, but differ when gas escapes. Neglecting bulk return while retaining a nonzero stellar yield is the intended trace-metal approximation.

The escape-speed scaling of wind mass loading makes the proposed trend dynamically plausible: the shallower Newtonian gravitational potential of a dwarf galaxy lets stellar feedback expel more gas per unit stellar mass. If coupled feedback supplies specific energy $e_{\ast }$, escape requires $\nu v_{\mathrm{esc}}^2/2\lesssim e_{\ast }$; fixed coupling suggests $\nu \propto v_{\mathrm{esc}}^{-2}$ for an energy-driven outflow. Fixed specific momentum $p_{\ast }$ instead gives $\nu v_{\mathrm{esc}}\lesssim p_{\ast }$, or $\nu \propto v_{\mathrm{esc}}^{-1}$ for a momentum-driven outflow.

Since escape velocity has $v_{\mathrm{esc}}^2\sim GM/R$, galaxies whose characteristic radius grows more slowly than their gravitating mass have increasing binding energy per unit mass with increasing mass. For $R\propto M^a$ with $a<1$, the energy-driven outflow estimate gives $\nu \propto M^{a-1}$. Equivalently, greater velocity dispersion in massive spheroids implies greater escape energy and more effective retention of metals.

For scale, one core-collapse supernova of $10^{51}\mathrm{erg}$ per $100M_{\odot }$ formed supplies about $5\times 10^{15}\mathrm{c}{\mathrm{m}}^2{\mathrm{s}}^{-2}$ before coupling losses. At ten percent coupling the energetic upper estimate $2e_{\ast }/v_{\mathrm{esc}}^2$ is about $110$ for $v_{\mathrm{esc}}=30\mathrm{km}{\mathrm{s}}^{-1}$, but about $0.4$ for $500\mathrm{km}{\mathrm{s}}^{-1}$. Efficient loss from shallow potentials and retention in deep potentials can span the required range. These are illustrative energetics rather than a unique prediction: cooling, entrainment, fallback and the actual initial mass function determine the escaping mass-loading factor.

Let $m_{\chi }$ be the dark matter particle mass, so its number density is $n_{\chi }=\rho /m_{\chi }$. Interpret $ϵ$ as total emitted power per particle, and assume isotropic optically thin dark-matter decay radiation; volume emissivity is $j_d=ϵ\rho /m_{\chi }$. At projected radius $R$, a line-of-sight coordinate $z$ has $r^2=R^2+z^2$. Writing $A=V^2/(2\pi G)$, the surface brightness isHere $z=R\tan \theta$ turns the integral into $\pi /R$. For the isothermal outer halo,If emissivity was defined per unit solid angle, the factor $1/(4\pi )$ is already included in $ϵ$; this only changes the normalization. The radial law holds between any core and outer truncation. For example, a phenomenological cored profile $\rho =A/(r^2+r_c^2)$ gives $I_d=ϵA/[4m_{\chi }\sqrt{R^2+r_c^2}]$: finite central surface brightness, but the same $R^{-1}$ outer law. This cored interpolation is not being asserted to be the exact isotropic constant-dispersion equilibrium.

Use a steady, isotropic, self-gravitating dark matter halo, with $V$ its constant one-dimensional velocity dispersion. The Spherical Jeans equation and spherical mass conservation readFor the scale-free outer region, put $\rho =A{r}^{-p}$. Its enclosed mass scales as $r^{3-p}$, so $GM(r)/r$ is constant only if $p=2$; the Spherical Jeans equation then gives $2V^2=GM(r)/r=4\pi GA$. Thus the singular isothermal sphere hasIts circular speed is $v_c=\sqrt{2}V$. If a speed denoted $V$ were instead the circular speed, the coefficient would be $V^2/(4\pi G)$, but the printed variable is a velocity dispersion.

The assumptions matter. A constant velocity dispersion and spherical shape alone do not specify a unique central mass density or orbital distribution. Even with isotropy, the Spherical Jeans equation integrates to $\rho ={\rho }_0\exp [-(\Phi -{\Phi }_0)/V^2]$, and the self-gravitating Poisson equation for Newtonian gravity givesThere are finite-central-density isothermal solutions as well as the singular one; the $r^{-2}$ expression is the intended outer scale-free halo approximation. A finite halo must also end or steepen at large radius, since the untruncated singular isothermal sphere has infinite total mass.

Let the extrapolated lower limit be $a=m_{\min }/M_{\odot }$ and keep the upper surviving limit at $x=1$. If the same cubic mass-luminosity relation is formally continued, the initial mass function givesMatching $1000$ therefore requires a planetary-scale cutoffThe mass per logarithmic interval is proportional to $m^{-1}$ for this initial mass function, so most of the added mass is near the cutoff. Substellar objects do not actually obey the unevolved cubic mass-luminosity relation: they cool, and their luminosity depends strongly on age. If they are approximated as entirely dark while only the original $0.1$--$1M_{\odot }$ stars supply light, then ${{\rm Y}}_{\ast }=(1/a-1)/0.9$ instead, giving $a=1/901\simeq 1.11\times 10^{-3}$. Both conventions lead to the same order-of-magnitude conclusion; the exact $10^{-3}$ result uses the question's formal power-law extrapolation.

The result is in the broad range expected for an old, relatively faint stellar population, but is high compared with many observed stellar mass-to-light ratios of a few solar units in optical bands. The comparison must specify the band and whether the measured mass includes dark matter: for example, early-type galaxy dynamical measurements find an $I$-band normalization near $3.8$ solar units at velocity dispersion $200\mathrm{km}{\mathrm{s}}^{-1}$, including a contribution from dark matter. That is an illustrative observational comparison, not a direct equality with a bolometric stellar calculation.

For an old stellar population, an upper surviving main sequence mass near $M_{\odot }$ is reasonable because more massive stars have shorter lifetimes; low-mass main sequence stars also contain much of the surviving stellar mass. However, red giants and other evolved stars can contribute a large fraction of the light despite their small numbers, so omitting them generally makes this calculated mass-to-light ratio too large. Omitting stellar remnants has the opposite effect on the inferred mass, and a single cubic mass-luminosity relation is only approximate across the stated range. Finally, a single $m^{-3}$ initial mass function down to $0.1M_{\odot }$ is steeper and more low-mass-heavy than commonly adopted distributions, which flatten at low mass. Ten solar units is a useful toy-model estimate, not a precision stellar-population prediction.

Write $x=m/M_{\odot }$ and absorb units into the initial mass function normalization, $dN=A{x}^{-3}dx$. The specified mass-luminosity relation is $L(x)=L_{\odot }{x}^3$, normalized to the Sun. Therefore the total stellar mass and luminosity areThe unknown normalization cancels. The mass-to-light ratio isLow-mass main sequence stars dominate the stellar mass, while equal intervals of $x$ contribute equal luminosity under these particular exponents. More generally the mass-to-light ratio of a power-law stellar population with these exponents on $a<x<b$ is $(M_{\odot }/L_{\odot })/(ab)$.

Objects near the derived cutoff are too faint for their absence in ordinary star counts alone to exclude a massive dark matter halo. The decisive test is gravitational microlensing: compact objects passing in front of background stars produce transient magnification even if they emit no detectable light. For a lens at distance $D_l$ and a source at $D_s$, the physical Einstein radius and event timescale areFor $D_l=10\mathrm{kpc}$, $D_s=50\mathrm{kpc}$ and $v_{\perp }=200\mathrm{km}{\mathrm{s}}^{-1}$, a $10^{-3}{M}_{\odot }$ lens gives $t_E\simeq 2.2$ days, within the regime tested by suitably sampled gravitational microlensing surveys. The microlensing optical depth isso at fixed compact-object mass density it is independent of the individual lens mass: reducing individual mass changes the event durations rather than hiding the total lensing probability.

EROS observations of the Magellanic Clouds found far fewer candidate events than a halo dominated by compact objects would predict and excluded such objects as the principal Galactic halo component over $6\times 10^{-8}<m/M_{\odot }<15$, under their tested halo models. Their often-quoted eight-percent bound is specifically near $0.4M_{\odot }$, not a uniform bound at every lens mass. A halo dominated by the approximately Jupiter-mass objects needed here is incompatible with those microlensing constraints. Deep optical and infrared counts also restrict a large population of luminous low-mass stars and brown dwarfs, while the cosmic baryon fraction provides a separate constraint on explaining all cosmological dark matter with baryonic objects.

First specify the dynamical estimate: assume an initially isolated virialized system, stars of different masses well mixed with the same velocity distribution, and ejecta escaping rapidly compared with a stellar crossing time, without imparting kicks to surviving stars. If $f$ is the retained stellar mass fraction, the original virial theorem gives $2T+W=0$, where $T$ is kinetic energy and $W<0$ is Newtonian gravitational potential energy. Immediately after loss, positions and velocities are unchanged, soThus impulsive disruption by stellar mass loss makes the energy of the whole remaining system nonnegative when at least half the initial mass is removed.

For the power-law initial mass function, the fraction removed is a mass fraction, not a number fraction:For integer slopes, $F(0)\simeq 0.990$, $F(1)\simeq 0.901$, $F(2)=0.333$, and $F(3)\simeq 0.009$. Interpolating between $\alpha =1$ and $2$ puts the threshold around $1.7$--$1.8$. Solving $F=1/2$ more accurately givesThe rearranged equation has a spurious root at $\alpha =2$ introduced when the vanishing denominator was cleared; the logarithmic limit $F(2)=1/3$ excludes it. Increasing $\alpha$ shifts the mass distribution towards lower masses, so the top-heavy side, $\alpha <{\alpha }_{\mathrm{crit}}$, has nonnegative remnant energy in the impulsive model; equality is marginal. One can see strict monotonicity from the covariance obtained by differentiating $F$: $F^{\prime }=-\mathrm{Cov}(1_{x>10},\ln x)<0$ under the normalized mass-weighted distribution.

The initial mass function alone does not guarantee complete disruption. If losses occur slowly, the adiabatic mass-loss expansion law gives expansion rather than the instantaneous half-mass threshold, and positive total energy can still leave a bound core after selective escape. The stated value is therefore the conventional global impulsive estimate, with no external confining Newtonian gravitational potential.

Massive stars rapidly return oxygen, magnesium, neon and other alpha elements through core-collapse supernovae, whereas Type Ia supernovae supply a substantial delayed iron contribution. If the galactic bulge and stellar halo formed much of their stellar mass rapidly, their stars locked in enhanced alpha-element to iron ratios before much delayed iron enrichment, while prolonged star formation in the local galactic disk included more iron-rich gas. The abundance contrast chiefly records different enrichment timescales, with the initial mass function and gas exchange also affecting the ratios.

The differential gravitational field of the Milky Way removes loosely bound outer material by tidal stripping, giving globular clusters and satellite dwarf galaxies a characteristic tidal radius. For a low-mass satellite of mass $m$ on a circular orbit of radius $R$ in a spherical host, the Jacobi tidal radius is $r_t\simeq R\{m/[(3-d\ln M_h/d\ln R)M_h(<R)]{\}}^{1/3}$; eccentric orbits often make the strongest restriction near pericentre. Tidal limitation explains the sharp outer falloff, although escaping stars can form extended tidal tails, so a fitted luminosity cutoff is not an exact boundary of all associated material.

The reaction sign and the meaning of recoil translation matter here. In Lighthill elongated-body theory, put $\mathcal{D}={\partial }_t+U{\partial }_x$. The transverse momentum imparted to the fluid per unit length is $m\mathcal{D}h$, so the lateral force on the fish is $-\mathcal{D}(m\mathcal{D}h)$. Thus body mass and positive added mass enter effective inertia with the same sign. The periodic trailing-edge thrust, to quadratic order in slope, isThe reaction sign and this small-amplitude thrust formula can also be checked in the primary derivation of the elongated-body force, equations (2.6), (3.2) and the discussion after (3.5). For the long-body limit below, the recoiling leading edge contributes only at relative order ${\alpha }^{-2}$.

The requested recoil and thrust formulas cannot all hold for the literal fish described. In particular, the printed thrust would remain $-m{U}^2/4$ when $h_1=0$, although a straight, unperturbed fish has $h_t=h_x=0$ and zero reactive thrust. Its negative term is missing an amplitude-squared slope factor. The stated effective translational coefficient $M^{\prime }-1$ also has the wrong sign for positive added mass, and the literal centre of mass translation cannot be equated to the translation of an arbitrarily chosen undeformed body reference line. The PDF defines $I^{\prime }=I/[m(\alpha L{)}^3]$; the converted TeX incorrectly uses the square.

Here are the force balance and angular momentum calculations for a consistent interpretation of the stated uniform fish. They give a corrected recoil solution, as well as a direct test of the printed equations. Defineand let $R$ be actual centre of mass displacement. The prescribed deformation has mean $a(t)=h_{1}f(t)/[2(\alpha +1)]$. Its mass-weighted mean must be removed in the recoil correction in elongated-body theory:Indeed $\int hdx=\ell R$, since $\int \xi dx=0$ and $\int h_{0}dx=\ell a$. Introduce the two geometric integralsFor uniform mass per unit length, $I={\mu }_{b}J$ at leading small-amplitude order. Retaining the symbol $I$ makes the rotational momentum balance transparent.

The needed integrals of transverse acceleration and convective derivatives areThe second spatial derivative includes the hinge's slope jump, as a distributional derivative; discarding it would incorrectly remove the $U^2$ terms. Newton's lateral equation is $M\ddot{R}=-m\int {\mathcal{D}}^{2}hdx$. The angular momentum of the deforming fish is $I\dot{\theta }+{\mu }_{b}C\dot{f}$, so its time derivative equals $-m\int \xi {\mathcal{D}}^{2}hdx$. Consequently the consistent uniform-added-mass recoil equations areThese follow directly from the local Lighthill elongated-body theory model with the physical reaction sign. An alternative reference-line translation replaces $R$ by $R-a$, and changes the shape forcing, but cannot change $M+m\ell$ to a difference of two positive inertias. Leading-edge details of a uniform-depth idealization can require an end model; they do not alter this reaction-sign test or the zero-stroke counterexample.

For $f=\cos \omega t$, write $R=\mathrm{Re}(r{e}^{i\omega t})$ and $\theta =\mathrm{Re}(\tau e^{i\omega t})$, choosing the periodic response and omitting free rigid-motion constants. The two balances give the explicit recoil amplitudesAt the tail the displacement and slope amplitudes are $H_L=h_1+r+\ell \tau /2-h_1/[2(\alpha +1)]$ and $S_L=h_1/L+\tau$. Thus a corrected, directly calculable mean thrust iswithin the trailing-edge approximation, with any retained leading-edge contribution treated consistently. This expression has the correct dimensions and vanishes quadratically with stroke amplitude.

A particularly simple independent check is the limit $U=0$. For the uniform fish, $I={\mu }_{b}J$, the balances give $\ddot{R}=0$ and $\theta =-Cf/J$ for the periodic response. The tail displacement amplitude is then exactlyHence, to the requested relative order,The leading-edge contribution is only $O({\alpha }^{-2})$ here. This recoil reduction contradicts the printed positive correction, for example for $M^{\prime }=2$, $I^{\prime }=1/6$. It also shows explicitly why treating the prescribed shape's moving mass as stationary would give a different model.

For completeness, the positive correction quoted in the question has a simple conditional origin if its two displayed recoil equations are accepted as formal equations rather than derived physical balances. Put $J^{\prime }=I^{\prime }+1/6$ and keep $U/(\omega L)$ fixed as $\alpha \to \infty$, away from $M^{\prime }=1$. Adding those equations and then using the first gives their leading periodic responseUnder the additional reference-line interpretation $h=h_0+R+\xi \theta$, the in-phase tail amplitude is $h_1[1+1/(2\alpha (M^{\prime }-1))+1/(8\alpha J^{\prime })]$. Squaring it in the trailing-edge thrust gives the printed frequency correction. Even this formal route, however, gives the amplitude-consistent resultnot a term $-m{U}^2/4$ independent of $h_1$. The quadrature part of $R$ first changes the squared displacement at relative order ${\alpha }^{-2}$, while recoil changes the tail slope at that order as well. This conditional calculation explains the intended algebra without presenting the flawed starting equations as a physical theorem.

Reactive thrust is positive when the tail's lateral-motion term exceeds its streamwise slope penalty; for an unrecoiling hinged tail this requires $\omega L>U$. A real steady swimming speed additionally balances the thrust against resistive drag force, absent from this inertial idealization. Recoil can alter either term and is not universally thrust-enhancing. A pole at $M^{\prime }=1$ is not a physical negative-added-mass resonance: the genuine combined inertias remain positive.

Work initially in units $G=c=\hslash =k_B=1$ and consider gravitational collapse settling to a nonextremal Schwarzschild black hole. Classically, outgoing rays sufficiently close to the forming event horizon suffer an arbitrarily large redshift. Quantum mechanically, the meaning of positive frequency at past null infinity differs from that at future null infinity; a Bogoliubov transformation between these two mode descriptions can therefore mix creation operators and annihilation operators. Hawking's argument computes that mixing and obtains a thermal occupation, even when the incoming state has no particles.

To see why the frequency mixing is universal, use final Schwarzschild time $t$, the Schwarzschild tortoise coordinate $r_{\ast }$, and retarded time $u=t-r_{\ast }$. Near the final event horizon, $r_{\ast }\sim (2\kappa {)}^{-1}\log (r-2M)$, with Schwarzschild surface gravity $\kappa =1/(4M)$. On a fixed advanced-time section, $u=v-2r_{\ast }$, so $r-2M$ is proportional to $e^{-\kappa u}$. Tracing a ray back through the regular collapsing region maps this small separation smoothly to its separation from the last escaping incoming ray. If $U$ is an affine incoming null coordinate and $U_H$ labels that limiting ray, the Hawking exponential ray map isThe logarithm at a simple horizon root fixes the exponent. Smooth propagation through the earlier collapse changes $A$ and phases but not the leading exponential redshift.

For a massless bosonic test scalar field, an outgoing mode of frequency $\omega >0$ at future null infinity behaves as $e^{-i\omega u}$. Trace it backward in the high-frequency geometric optics approximation. Its relevant incoming dependence is, up to normalization and a phase,Here $\Theta$ is the Heaviside step function. Incoming rays with $U>U_H$ do not escape. This truncated logarithmic phase is not purely positive frequency with respect to $U$. Decomposing it into incoming positive-frequency solutions $e^{-i{\omega }^{\prime }U}$ and negative modes $e^{+i{\omega }^{\prime }U}$ gives the Bogoliubov transformation. With $x=U_H-U$, the two relevant Fourier transform integrals areThe former is the positive-frequency coefficient and the latter is the negative-frequency coefficient; the Klein-Gordon inner product introduces the same frequency normalization into their modulus ratio. By the defining integral of the gamma function and analytic continuation in its Laplace parameter,For a complex number $z$ in the right half-plane, $|z^{-1-ia}|=|z{|}^{-1}{e}^{a\arg z}$. As $ϵ\downarrow 0$, the two arguments tend to $+\pi /2$ and $-\pi /2$. This gives the thermal ratio of Hawking Bogoliubov coefficients:The regulator fixes the branches and is essential to the sign of the thermal exponent.

Expand the outgoing annihilator as $b=\int d{\omega }^{\prime }({\alpha }_{\omega {\omega }^{\prime }}^{\ast }{a}_{{\omega }^{\prime }}-{\beta }_{\omega {\omega }^{\prime }}^{\ast }{a}_{{\omega }^{\prime }}^{†})$, including a complete set of incoming channels. The canonical identities for a bosonic Bogoliubov transformation express $[b,b^{†}]=1$. For normalized outgoing wave packets, they give $\int (|\alpha {|}^2-|\beta {|}^2)d{\omega }^{\prime }=1$. In the incoming Fock vacuum, the outgoing mean of the particle number operator is $N_{\omega }=\int |\beta {|}^{2}d{\omega }^{\prime }$. Combining the two identities with the thermal ratio, in the late-time narrow-frequency packet limit, yieldsWave packets avoid interpreting the divergent normalization of infinitely long continuum modes as a finite particle count. The temperature is the Hawking temperature, and the occupation is the zero-chemical-potential Bose-Einstein distribution. Restoring constants, the Schwarzschild black hole temperature is $T_H=\hslash c^3/(8\pi GM{k}_B)$ for physical mass $M$.

The preceding thermal occupation describes the horizon-originating channel before exterior scattering. The exterior angular and curvature potential partly reflects it. For a scalar field, the number flux at future null infinity is consequentlywhere the greybody factor ${\Gamma }_{\ell }$ is the transmission probability of a partial wave. Multiplication by $\omega$ gives its energy flux. Thus the asymptotic spectrum has the thermal denominator but is not an exact blackbody spectrum. For fermionic fields the canonical anticommutation relations replace the minus sign in the occupation denominator by a plus sign, at the same Hawking temperature.

The relevant state is the collapse vacuum for Hawking radiation: no incoming thermal bath at past null infinity, and regular short-distance behavior in freely falling coordinates at the forming event horizon. It predicts an outgoing flux, unlike equilibrium with a thermal bath on an eternal black hole. The derivation uses quantum field theory in curved spacetime on a prescribed classical geometry; it does not require outgoing particles to follow forbidden classical paths out of the interior. Positive energy carried to infinity is accompanied, in the semiclassical description, by negative Killing-energy flux into the black hole. Including that flux in slow backreaction decreases its mass, giving black-hole evaporation.

The approximation applies to the late radiation of a large nonextremal black hole while its surface gravity changes slowly on a time scale ${\kappa }^{-1}$. Backward propagation involves very high locally measured frequencies, so the calculation assumes the usual regular short-distance quantum-field state. It determines the leading radiation and temperature, not the Planck-scale endpoint of black-hole evaporation. The exponential horizon redshift forces positive/negative frequency mixing with a thermal ratio, giving outgoing Hawking radiation at $T_H=\kappa /(2\pi )$.

Use mostly-plus Lorentzian signature and $\{{\gamma }^{\mu },{\gamma }^{\nu }\}=2{\eta }^{\mu \nu }$, with antisymmetrized gamma products of unit weight. A Majorana spinor obeys a reality condition under charge conjugation:A Weyl spinor instead has definite chirality,A four-dimensional Majorana spinor has four real components, while a Weyl spinor has two complex components. Charge conjugation reverses chirality in this signature, so a nonzero spinor cannot satisfy both conditions simultaneously.

For the Dirac equation $(/\partial +m)\psi =0$, the derivative term reverses chirality because ${\gamma }_5$ anticommutes with every gamma matrix, whereas the mass term preserves it. Projecting a purely chiral spinor equation onto its two chiralities therefore gives both $/\partial \psi =0$ and $m\psi =0$. Thus a nonzero four-component Weyl spinor satisfying this uncoupled Dirac equation must be massless. A Majorana spinor contains conjugate opposite-chirality components and can have a nonzero mass. Its charge-conjugation reality condition identifies particle and antiparticle operators, so a Majorana particle is its own antiparticle. A Weyl field does not by itself impose that identification.

This does not forbid writing a Majorana mass term for a neutral two-component Weyl field: such a term couples the field to its conjugate, and the resulting four-component massive field is Majorana, not a purely chiral solution of the uncoupled Dirac equation.

For the Rarita-Schwinger field, set $\chi ={\gamma }^{\mu }{\psi }_{\mu }$ and $d={\partial }^{\mu }{\psi }_{\mu }$. The antisymmetric kinetic operator is invariant underbecause commuting derivatives contract to zero against ${\gamma }^{\mu \nu \rho }$. Expanding its gamma product givesGamma contraction gives ${\gamma }_{\mu }{E}^{\mu }=2(/\partial \chi -d)$, so the equation implies $d=/\partial \chi$ and $/\partial {\psi }^{\mu }={\partial }^{\mu }\chi$.

For a Fourier mode with momentum $k$, if $k^2\ne 0$ the latter identity gives $k^2{\psi }^{\mu }=k^{\mu }/k\chi$. The mode is proportional to $k^{\mu }$ and hence pure gauge. Every physical nonzero-momentum mode therefore has $k^2=0$. For a null mode, choose an index with $k^{\mu }\ne 0$. The same equation makes $\chi$ an element of the image of $/k$, allowing a gauge choice with $\chi =0$. The field equations then reduce toResidual gauge transformations obey $/\partial \eta =0$.

Take a null momentum in the third spatial direction. Residual gauge freedom and transversality remove the temporal and longitudinal vector components, leaving two transverse vector polarizations. Their helicities are $\pm 1$, while on-shell spinor helicities are $\pm 1/2$. The gamma-trace constraint removes the two combinations of total helicity $\pm 1/2$ and retains the aligned combinations of helicity $\pm 3/2$. Explicitly, in a chiral representation the transverse gamma contractions use ${\sigma }_1\pm i{\sigma }_2$, which annihilate the correspondingly aligned spin states and act nontrivially on the opposite states. Thus a Majorana Rarita-Schwinger field has two physical massless polarizations. The Majorana assumption is the one relevant to simple supergravity; a complex vector-spinor has these helicities for its particle and independent antiparticle modes.

A consistent free massive deformation is the algebraic termadded to a kinetic term $-\frac{1}{2}{\overline{\psi }}_{\mu }{\gamma }^{\mu \nu \rho }{\partial }_{\nu }{\psi }_{\rho }$. Its equation isFor $m\ne 0$, divergence and gamma contraction give, respectively,These imply the massive Rarita-Schwinger constraints $\chi =0$ and $d=0$, and the remaining equation is $(/\partial +m){\psi }_{\mu }=0$. Squaring gives $(\square -m^2){\psi }_{\mu }=0$ in the chosen signature. The gauge symmetry is lost. In the rest frame, transversality sets the temporal component to zero; three spatial vector components times two on-shell spinor states give six states, and gamma trace removes the two spin-one-half states. The result is four massive spin-$3/2$ states, with spin projections $\pm 3/2,\pm 1/2$. In supergravity the additional longitudinal states can be supplied by the super-Higgs mechanism.

Use the mostly-plus Clifford algebra and Riemann curvature tensor conventions of the preceding solutions, with $e=\det (e_{\mu }{}^a)$, ${\kappa }^2=8\pi G$, and a Majorana gravitino. An action for minimal four-dimensional supergravity in the exact transformation normalization given isThe kinetic coefficient is tied to the normalization of the supersymmetry transformations: if ${\Psi }_{\mu }=2{\psi }_{\mu }$ and $\eta =2ϵ$, this becomes the canonical coefficient $-1/2$, with $\delta e_{\mu }{}^a=(\kappa /2)\bar{\eta }{\gamma }^a{\Psi }_{\mu }$ and $\delta {\Psi }_{\mu }={\nabla }_{\mu }\eta /\kappa$. Mixing those two normalizations would spoil the cancellation.

To the retained order, use the torsion-free connection. More generally the 1.5-order formalism allows its variation to be omitted after imposing its own algebraic equation; the resulting contorsion is quadratic in fermions and first affects the displayed action at quartic order. The metric tensor variation induced by the vierbein is $\delta g_{\mu \nu }=4\kappa \overline{ϵ}{\gamma }_{(\mu }{\psi }_{\nu )}$, so the Einstein-Hilbert action varies asFor the fermionic action, varying both Majorana factors gives the same contribution after integration by parts and the Grassmann bilinear interchange rule. ThusThe vector-index Riemann curvature tensor contribution vanishes by the first Bianchi identity. The spin-index contribution is determined byFor the last equality, expand the antisymmetric three-gamma/two-gamma product: the totally antisymmetric Riemann curvature tensor terms vanish by Bianchi, and the surviving terms are the Ricci contraction minus half its trace. Hencewhich cancels the gravitational variation. Variations of vierbeins, gamma matrices and the connection in the fermionic kinetic term carry three dynamical fermions. They lie beyond this linear-fermion cancellation and are cancelled by the omitted four-fermion completion. Thus the action is invariant up to boundary terms at the order specified, rather than claiming that a truncated action is exactly invariant at every fermion order.

For the cosmological extension, add a gravitino mass parameter $m$ and deform the derivative:Indeed the fermion terms combine as $-2{\bar{\psi }}_{\mu }{\gamma }^{\mu \nu \rho }{\mathcal{D}}_{\nu }{\psi }_{\rho }$ because ${\gamma }^{\mu \nu \rho }{\gamma }_{\nu }=-2{\gamma }^{\mu \rho }$. The modified spinor connection curvature isand contraction givesThe same action-variation cancellation therefore works with $G$ replaced by $G+\Lambda g$ precisely forThis is the Anti-de Sitter deformation of minimal supergravity. A cosmological term alone would leave its variation uncancelled; the mass term and the gamma shift of the gravitino transformation are both necessary. A positive cosmological constant cannot arise from this real unbroken minimal deformation.

For electromagnetism, retain $\mathcal{N}=1$ and introduce an Abelian vector multiplet $(A_{\mu },\lambda )$, with Majorana gaugino $\lambda$. A complete quadratic-fermion extension isThe accompanying leading transformations areThey preserve the rigid Maxwell-gaugino action for constant $ϵ$. For a local parameter, the variation containsTo obtain this, vary $F^2$, vary both gaugino factors, and use ${\gamma }^{ab}{\gamma }^{\mu }={\gamma }^{ab\mu }+g^{b\mu }{\gamma }^a-g^{a\mu }{\gamma }^b$. The terms with $\nabla \lambda$ cancel between the two actions; integrating the three-gamma term uses ${\nabla }_{[\mu }{F}_{\nu \rho ]}=0$, leaving exactly the displayed derivative of the parameter. The interaction is $-\kappa e{\bar{\psi }}_{\mu }{j}^{\mu }$, whose variation under $\delta {\bar{\psi }}_{\mu }={\nabla }_{\mu }\overline{ϵ}/\kappa$ cancels that term.

The other linear-fermion variation is also fixed. Put $\mathcal{F}=F_{ab}{\gamma }^{ab}$ and $T_{\mathrm{EM}}^{\mu \nu }=F^{\mu \rho }{F}^{\nu }{}_{\rho }-g^{\mu \nu }{F}^2/4$. Clifford multiplication givesVarying $\lambda$ in the gravitino-current interaction then gives $2\kappa e{T}_{\mathrm{EM}}^{\mu \nu }{\bar{\psi }}_{\mu }{\gamma }_{\nu }ϵ=-2\kappa e{T}_{\mathrm{EM}}^{\mu \nu }\overline{ϵ}{\gamma }_{\nu }{\psi }_{\mu }$, which cancels the Maxwell vierbein variation. This proves the required supergravity coupling of an Abelian vector multiplet through the retained order.

Beyond that order one uses a supercovariant gauge field strength, here ${\hat{F}}_{\mu \nu }=F_{\mu \nu }-2\kappa {\bar{\psi }}_{[\mu }{\gamma }_{\nu ]}\lambda$, and the associated higher-fermion transformation terms and four-fermion interactions. An auxiliary scalar may be retained for off-shell closure, or eliminated in this on-shell formulation. A Maxwell field without a fermionic partner cannot be coupled by merely adding $-F^2/4$ while leaving all transformations unchanged. An alternative is an $\mathcal{N}=2$ gravity multiplet, where the vector is a graviphoton and a second gravitino supplies its extra fermionic states; the explicit construction above is the minimal $\mathcal{N}=1$ matter extension.

Use $s=p+1$ and the exterior mean-motion resonance condition $s{n}_r=p{n}_{\mathrm{pl}}$. Then the resonant semi-major axis is $a_r=a_{\mathrm{pl}}(s/p{)}^{2/3}$, neglecting the small planet-to-star mass correction. For the geometric construction, take $a_{\mathrm{pl}}=n_{\mathrm{pl}}=1$, $p=4$, $e=0.2$, $\varpi =\pi /4$ and start the particle at pericentre. Solve Kepler's equation $M=E-e\sin E$ with $M=(4/5)t$. Its coordinates in the rotating reference frame areFour particle orbital periods correspond to five planetary orbital periods, so the curve closes. At a pericentre, the mean longitude and the true longitude both equal $\varpi$, and the resonant argument becomes $\varphi =4(\varpi -{\lambda }_{\mathrm{pl}})$. Thus the four pericentre directions are $\varphi /4$ modulo $\pi /2$. The figure chooses $\varphi =\pi$ and marks the corresponding $\pi /4$ angle to the planet; other phases rotate this pattern. Its second panel shows an instantaneous conjunction at apocentre, where the angle from pericentre to the common longitude is directly $\varphi =\pi$: here both true and mean anomaly equal $\pi$. The loops near pericentre reflect the fact that the particle can temporarily overtake the planet in angular speed, despite its lower mean motion.