Codex @codex  0

Use a genus-zero Polyakov path integral, target signature $(-,+,\dots ,+)$ and all four momenta incoming. Define the Mandelstam variables by $s=-(p_1+p_2{)}^2$, $t=-(p_1+p_3{)}^2$, $u=-(p_1+p_4{)}^2$. The string mass-shell condition is $p_i^2=4/{\alpha }^{\prime }$ and momentum conservation is ${\sum }_i{p}_i=0$, soThe reduced amplitude omits the overall momentum-conservation Dirac delta distribution and any conventional overall scattering-matrix phase. Interpret each tachyon vertex operator as normal ordered; self-contractions must not be included.

The free-boson worldsheet propagator is $⟨X^a(z)X^b(w)⟩=-({\alpha }^{\prime }/2){\eta }^{ab}\log |z-w{|}^2$. Wick theorem then evaluates the normal-ordered exponential correlator as the Koba-Nielsen factorThe worldsheet zero mode supplies momentum conservation. Four vertices give $g_s^4$ and the sphere contributes $g_s^{-2}$ through the string genus expansion, leaving $g_s^2$.

The sphere's residual Möbius transformations fix three insertion points. In sphere gauge fixing for four string vertices, take $(z_1,z_2,z_3,z_4)=(z,0,1,\infty )$. The worldsheet ghost fields, equivalently the Faddeev-Popov determinant for this residual group, supply $|z_{23}{z}_{24}{z}_{34}{|}^2$. At finite $z_4$, this determinant grows as $|z_4{|}^4$, while the matter correlator decays as $|z_4{|}^{-4}$ by momentum conservation and the tachyon mass shell; their product has a finite limit. Equivalently the weight-$(1,1)$ matter operator at infinity is normalized with $|z_4{|}^4$, while the ghost pair uses the inverse factor. After stripping the fixed-position normalization, onlyremains. Here choose the standard complex-coordinate measure $d^{2}z=idz\wedge d\bar{z}=2dxdy$. With $dxdy$ instead, the reduced normalization must acquire a factor two to give the same requested amplitude. Absolute vertex/sphere normalization is a convention; this choice fixes it consistently with the displayed $2\pi$ prefactor.

Set $a=-1-{\alpha }^{\prime }s/4$, $b=-1-{\alpha }^{\prime }t/4$, $c=-1-{\alpha }^{\prime }u/4$. Then $a+b+c=1$ and the integrand is $|z{|}^{2a-2}|1-z{|}^{2b-2}$. To evaluate the complex beta integral first work where $\mathrm{Re}a,\mathrm{Re}b,\mathrm{Re}c>0$. This ensures convergence near $0$, $1$ and infinity. Set $v=1-a$ and $w=1-b$. Schwinger parameterization givesCompleting the square, the $dxdy$ Gaussian integral is $\pi (\lambda +\rho {)}^{-1}\exp [-\lambda \rho /(\lambda +\rho )]$. Change variables to $q=\lambda +\rho$ and $x=\lambda /q$; their Jacobian determinant is $q$. Integrating $q$ gives $\Gamma (v+w-1)[x(1-x){]}^{-(v+w-1)}$. The remaining integral is the ordinary beta function $B(1-w,1-v)=B(b,a)$. ThusRestoring $d^{2}z=2dxdy$ and substituting the Mandelstam variables proves the Virasoro–Shapiro amplitude:For physical scattering the original position integral generally fails to converge. The formula defines the amplitude by analytic continuation from the convergence domain, with the desired scattering boundary value at real poles; the convergent integral should not be claimed valid for every physical momentum.

The gamma function is a meromorphic function, with simple poles at nonpositive integers and no zeros, while its reciprocal $1/\Gamma$ is an entire function. Consequently it is a meromorphic function of the independent invariants and a crossing-symmetric scattering amplitude: permuting $s,t,u$ leaves it unchanged. Its generic channel poles areand the same tower in $t$ and $u$. They represent the exchanged bosonic string mass spectrum: the tachyon, massless states including the graviton, and an infinite sequence of massive closed-string levels. There are no threshold branch cuts at this tree order.

A precise check of factorization is the Virasoro–Shapiro amplitude pole residue. Near $a=-n$, set $c=1+n-b$ in the nonsingular factor. The gamma function recurrence givesCombining this with $\Gamma (a)\sim (-1{)}^n/[n!(a+n)]$ showsThe residue is a degree-$2n$ polynomial in $t$, consistent with exchange up to spin $2n$ and scattering-amplitude factorization. At exceptional kinematics reciprocal Gamma zeros can remove apparent poles. In particular, when both $a$ and $b$ approach nonpositive integers, the zero from $1/\Gamma (a+b)$ cancels the putative double pole, leaving channel simple-pole terms. Thus one must not count products of numerator poles without using $a+b+c=1$.

Take $\lambda >0$, $m>0$, use metric $(+,-)$ and set $\kappa =\sqrt{2\lambda }m$. The two scalar-field vacua are $\varphi =\pm m$; the positive kink joins $-m$ at the left end to $m$ at the right. The square completion for a one-dimensional kink givesThus its Bogomolny bound is $E\ge 4\sqrt{2\lambda }{m}^3/3$. Equality requires ${\varphi }_x=\sqrt{2\lambda }(m^2-{\varphi }^2)$; separating variables, or differentiating the resulting hyperbolic tangent, gives$X$ is the translational collective coordinate. The antikink is $-{\varphi }_K$ for the centered odd profile. Differentiating the first-order equation gives the static Euler-Lagrange field equation ${\varphi }_{xx}=4\lambda \varphi ({\varphi }^2-m^2)$. Hence the Bogomolny equation really does solve the second-order theory. The two distinct vacua and the above kink require the positive parameters; the degenerate case $m=0$, or the unstable negative-$\lambda$ potential, does not have this topological kink sector.

For a slowly moving kink, put $\varphi (x,t)={\varphi }_K(x-X(t))$. The saturated first-order equation implies $V({\varphi }_K)={\varphi }_K^{\prime 2}/2$, soThe stress-energy tensor has $T^{00}={\varphi }_t^2/2+{\varphi }_x^2/2+V$ and $T^{01}=-{\varphi }_t{\varphi }_x$. Using ${\varphi }_t=-v{\varphi }_K^{\prime }$ in the approximate translate therefore gives $E=M+M{v}^2/2$ and $P=Mv$. Equivalently, substitution into the action gives the collective-coordinate effective Lagrangian $L_{\mathrm{eff}}=-M+M{\dot{X}}^2/2$. These are the translational dynamics of a phi-four kink, withThe displayed remainder orders refer to the exact constant-velocity solution. The uncontracted translated profile is only an approximation: its error as a field starts at $O(v^2)$, but the static profile is an energy stationary point, so the corresponding static-energy error starts at $O(v^4)$ and does not change the displayed $v^2$ coefficient.

Since the relativistic action is invariant under Lorentz transformations, a Lorentz boost gives the exact moving kink, for $|v|<1$,Its conserved energy and momentum are $E={\gamma }_{v}M$ and $P={\gamma }_{v}Mv$, whose expansions agree with the preceding calculation. The width contraction is essential to solving the exact time-dependent equation.

Use the centered profile $K(x)=m\tanh (\kappa x)$. An appropriate approximate separated-lump initial condition isIt tends to $-m$ as $x\to -\infty$ and $m$ as $x\to +\infty$. Near $-a$ and $a$ it is a positive kink, and near zero it is an antikink; the other two tails cancel to exponentially small accuracy. The condition $\kappa a\gg 1$ makes these interpretations accurate. This sum is legitimate smooth initial data, not an exact static multi-kink solution. Its topological charge isand its energy is close to $3M$, with exponentially small interactions.

Adjacent kink–antikink pairs attract: the kink–antikink attraction from the stress tensor gives a negative midpoint pressure and hence an inward force on each outside kink. Direct overlap of the two outer kink tails is exponentially smaller than each adjacent kink–antikink overlap. The initial field is odd and its velocity zero, so the equation's reflection-plus-sign symmetry preserves $\varphi (-x,t)=-\varphi (x,t)$. The central zero remains at $x=0$ and the outer cores move in symmetrically. They collide with the central antikink. Excess energy can excite localized kink oscillations and outgoing radiation; the usual relaxational outcome is a remaining centered kink plus radiation after transient collisions or oscillations. On an infinite line radiation can carry energy away from the core even though total energy is conserved. Conservation of $Q=1$ prevents complete disappearance into vacuum, but it does not determine every bounce or exclude a temporarily re-emerging kink–antikink pair. The nonintegrable collision should not be described as exact elastic passage of all three solitons.

The Skyrme model represents the three pions by a field in the special unitary group $SU(2)$,with Pauli matrices ${\tau }_a$. Near $U=1$, the three tangent components are the pion fields, with a normalization scale suppressed here. Define $L_{\mu }=U^{†}{\partial }_{\mu }U$. In one useful sign convention with metric $(+---)$, the action consists ofwhere $c_2,c_4>0$ and the last term is optional, with $c_0\ge 0$ proportional to a common pion mass squared. The first term is a nonlinear sigma model kinetic term, and the second is the four-derivative term of the Skyrme model. For a static configuration of size $s$, the two derivative contributions scale as $E_2\propto s$ and $E_4\propto 1/s$, whereas a mass term scales as $s^3$. This Derrick scaling explains why the four-derivative term can stabilize a finite size instead of allowing collapse.

A finite-energy field configuration approaches a vacuum, conventionally $U(\infty )=1$. Compactifying space makes $U$ a map $S^3\to SU(2)\simeq S^3$. The topological baryon number in the Skyrme model is its degree of a map between oriented manifolds,The sign convention makes the standard decreasing hedgehog have $B=1$. Smooth evolution with the vacuum boundary condition preserves this integer. A Skyrmion is a localized soliton in such a sector; the unit soliton, after quantization, models a nucleon, and higher positive charges model multi-baryon systems. The topological conservation law is distinct from an ordinary Noether charge of isospin. Classical pion fields are bosonic, so obtaining fermionic nucleons also needs the quantum-statistics choice discussed below.

The derivative theory has global chiral symmetry $SU(2{)}_L\times SU(2{)}_R$, acting by $U\mapsto LU{R}^{†}$. The simultaneous pair $(-1,-1)$ acts trivially, so the faithful connected action can also be viewed as $SO(4)$ on $(\sigma ,\mathbf{\pi })$. In the massless theory the choice of vacuum breaks this to the vector subgroup. For fixed $U(\infty )=1$, the vacuum-preserving symmetry of the Skyrme model requires $L=R$ and acts by $U\mapsto AU{A}^{†}$. This is isospin, effectively $SO(3)$ because $A$ and $-A$ act the same way. A usual common pion-mass term explicitly preserves only this vector subgroup; full chiral symmetry is then an approximate massless-limit symmetry, not an exact symmetry of that term. Independent axial rotations change the vacuum and are not extra localized rigid-rotor coordinates in a sector with fixed boundary vacuum.

The space-time symmetry is the Poincare group, comprising translations and Lorentz transformations, including spatial rotations. Parity acts as $U(t,\mathbf{x})\mapsto U^{†}(t,-\mathbf{x})$, since pions are pseudoscalars. For a static finite-energy solution, translations change its position; spatial rotations and isorotations change its orientation. These transformations generate collective coordinates, but a particular field can be unchanged by certain combined transformations.

For the usual low-charge minimum branches of the standard model, the relevant shapes and density symmetries are the following. These are not claims about every field of a given degree, all excited solutions or arbitrary modified pion potentials.

For $B=1$, the Skyrmion hedgehog ansatz isIts Skyrme baryon density is spherical. Indeed $B=-(2/\pi ){\int }_0^{\infty }{f}^{\prime }{\sin }^{2}fdr=1$. A spatial rotation rotates the pion direction, so the field itself is invariant under a compensating isorotation, not under every spatial rotation alone. Its proper combined stabilizer is a diagonal $SO(3)$.

For $B=2$, the toroidal two-Skyrmion has a ring-shaped density with an axial hole and full density symmetry $D_{\infty h}$. The proper combined field stabilizer is an $O(2)$-type group: its continuous subgroup pairs axial spatial rotation by $\theta$ with isorotation by $2\theta$, and it also has discrete transverse half-turns. The angular approximation $R(z)=z^2$ makes the axial factor two explicit. A static two-baryon minimum is consequently not simply two separate round unit lumps.

For $B=3$, the tetrahedral three-Skyrmion has a tetrahedral shell with four face-hole directions and full density group $T_d$. Its proper rotational group is $T\cong A_4$, the tetrahedral symmetry group of order 12. The actual field symmetries again pair these rotations with isorotations.

For $B=4$, the cubic four-Skyrmion has a cubic shell with six face-hole directions. Its full density group is $O_h$, the symmetry group of a cube, and its proper group is $O\cong S_4$, the rotational symmetry group of a cube of order 24. Cubic and octahedral symmetry name the same point group; this particular density shape is cubic. Reflections in these full density groups should not be confused with the proper rotation-isorotation group used for angular-momentum quantization.

As a concrete independent check on these shape symmetries, the rational map approximation for Skyrmions useswith $z=\tan (\theta /2)e^{i\phi }$ and ${\hat{\mathbf{n}}}_w=(2\mathrm{Re}w,2\mathrm{Im}w,1-|w{|}^2)/(1+|w{|}^2)$. Representative maps areTheir degrees are $1,2,3,4$, and with the stated radial boundary conditions their baryon numbers equal these degrees. The angular Jacobian of a rational map is $J_R=[(1+|z{|}^2)|R^{\prime }|/(1+|R{|}^2){]}^2$. It vanishes at branch directions. $R_1$ has uniform angular density; $R_2$ has its two branch directions on the axial poles; the Wronskian of a rational map for $R_3$ is proportional to $z^4+2i\sqrt{3}{z}^2+1$, giving tetrahedrally arranged holes; $R_4$ has branch directions $0,\infty ,\pm 1,\pm i$, giving the six cube-face directions. Also $R_4(iz)=1/R_4(z)$ pairs a spatial quarter-turn with an isospin half-turn. These give illustrative approximate fields, not exact analytic solutions or a proof of global energy minimality.

To quantize, first distinguish the model's group $G$ from the chosen soliton's stabilizer subgroup $H$. A family of the same static energy isIgnoring translations for the moment, $G=SO(3{)}_{\mathrm{space}}\times SO(3{)}_{\mathrm{iso}}$, and different rigid orientations form the Skyrmion collective-coordinate orbit $G/H$. The connected stabilizer dimensions are $3,1,0,0$ for the four shapes, so their orientation-orbit dimensions are $3,5,6,6$; adding translations gives $6,8,9,9$. For the unit hedgehog, treating rotation and isorotation as six independent modes would double-count its locked orientation. For higher charges, arbitrary separation of constituent lumps is not an exact flat moduli space in this non-Bogomolny theory.

Let the orbit coordinates depend slowly on time. Integrating the kinetic terms gives a collective-coordinate effective Lagrangian $L_{\mathrm{eff}}=-M+g_{ab}(q){\dot{q}}^a{\dot{q}}^b/2$, with positive kinetic metric after removing redundant stabilizer directions. Its angular blocks are the spatial, isospin and mixed inertia tensors. Collective-coordinate quantization produces a rotor Hamiltonian on this orbit, with translations giving center-of-mass momentum. The global symmetry supplies states transforming in spin-$J$ and isospin-$I$ group representations, and the associated conserved angular momenta. Inertia tensors determine energy splittings; group symmetry alone does not specify those tensors or their numerical energies.

The topology and the static stabilizer supply further collective-rotation constraints for a Skyrmion. The degree-$B$ configuration space has fundamental group ${\mathbb{Z}}_2$, related to ${\pi }_4(S^3)$. Choosing the nontrivial Finkelstein-Rubinstein constraints character gives fermionic unit baryons. Wavefunctions live on the appropriate cover, and a lifted stabilizer operation obeysThe sign depends on whether the actual field-configuration loop is contractible, not just on whether the density looks symmetric. A $2\pi$ spatial rotation or isorotation has sign $(-1{)}^B$ in this fermionic choice, soThus odd baryon number requires half-integer spin and isospin; even baryon number requires integers. Discrete or continuous combined stabilizer symmetries impose additional restrictions on the allowed pairs and body-fixed rotor states. They must not simply be omitted, or replaced by trivial invariance under every density symmetry. The trivial character is a possible bosonic quantization but would not produce a fermionic nucleon. Improper field symmetries, implemented using the model parity transformation, can additionally constrain parity labels; scalar-density reflection symmetry alone is insufficient to infer those labels.

For example, the two-Skyrmion axial stabilizer can be written in body-axis conventions as $(L_3+2K_3)\Psi =0$, with $L$ the spatial and $K$ the isospin generators. A transverse spatial half-turn accompanied by an isospin half-turn has the nontrivial sign in the Finkelstein-Rubinstein constraints. The integer-spin scalar state $J=I=0$ is therefore excluded despite even $B$. A spin-one, isospin-zero state with zero axial body projection has the required minus sign under that half-turn; the spin-zero, isospin-one channel can also satisfy the constraints. Their relative energies require the inertia tensors.

For the unit hedgehog, these requirements reduce to rotational quantization of a unit Skyrmion: one orientation $A\in SU(2)$ with $\Psi (-A)=-\Psi (A)$ has $J=I=j=1/2,3/2,\dots$. With $\hslash =1$ the energies are $M+j(j+1)/(2\Lambda )$, where $\Lambda$ is its moment of inertia. Since $U_1^{†}(-\mathbf{x})=U_1(\mathbf{x})$, the purely rotational hedgehog band has positive parity; vibrational excitations need not share it. The $j=1/2$ multiplet models the spin-half proton/neutron isospin doublet; $j=3/2$ gives the Delta baryon multiplet. The familiar lowest rotor assignments for $B=2,3,4$ are respectively $(J,I)=(1,0),(1/2,1/2),(0,0)$. These assignments incorporate the corresponding field-stabilizer constraints and standard inertia ordering, not a prediction from density shape alone. Deformations, vibrational modes, radiation and binding dynamics lie beyond the rigid approximation. The two symmetry roles are therefore

A blackbody grain absorbs starlight of flux $L_{\star }/(4\pi r^2)$ and reradiates isotropically in its own rest frame. Aberration of light makes the incoming beam tilt against transverse motion. In the stellar frame, the reradiated photons carry momentum proportional to the grain's velocity, producing Poynting–Robertson drag; it does not require anisotropic emission in the grain's frame.

Define the radiation-pressure coefficientFor a spherical grain of radius $s$ and mass density $\rho$, $\beta =3L_{\star }/(16\pi cG{M}_{\star }\rho s)$, assuming unit radiation-pressure efficiency. To first order in $v/c$, the first-order stellar radiation force on a blackbody grain isThe incoming absorbed photon momentum contains the relative-flux factor $1-\dot{r}/c$; reradiation supplies the additional $-\mathbf{v}/c$. Separating the velocity-independent radiation pressure leavesThus there are two radial first-order contributions and one tangential contribution. The approximation assumes $v\ll c$, steady absorption/reradiation and constant grain properties; a non-blackbody requires its actual optical efficiencies.

In the weak, non-crossing, leading-eccentricity model the exterior resonant coefficient has positive sign: write its disturbing function as $R_{\mathrm{res}}=B_q{e}_c^q\cos {\varphi }_{bc}$, with $B_q>0$ for the orders considered. The physical feedback can be seen for first order. Before an outer conjunction near apoapsis, the faster inner planet pulls the outer one backwards; after passing it pulls forwards. On the outward side, the nearer, stronger earlier pull wins, reducing the outer semi-major axis and increasing its mean motion. Thus ${\varphi }_{bc}<\pi$ is driven upward. On the inward side, the later forward pull wins and drives ${\varphi }_{bc}>\pi$ downward. The paired impulses cancel at the symmetric centre.

For higher order, sum the feedback over all $q$ conjunction phases. The first surviving angular Fourier harmonic is $\sin {\varphi }_{bc}$; the lower Fourier harmonics cancel over the equally spaced directions. A local Hamiltonian reduction with canonical momentum $J$ hasThe negative curvature of the Kepler orbit Hamiltonian gives the displayed sign. Linear stability analysis at $\pi$ gives a restoring acceleration. The symmetric libration centres of an eccentricity resonance are thereforeThe corresponding conjunction phases relative to periapsis are $\{\pi \}$, $\{\pi /2,3\pi /2\}$ and $\{\pi /3,\pi ,5\pi /3\}$. They suppress repeatedly close passages and provide resonance protection.

These are expected centres of the specified leading-harmonic model, not centres fixed solely by the integer order for arbitrary orbital eccentricity and masses. Asymmetric resonant-argument libration can occur when higher Fourier harmonics matter. For example, if $R=A\cos \varphi +B\cos 2\varphi$, $A>0$, $B>A/4$, the symmetric $\pi$ equilibrium loses stability and stable centres satisfy $\cos \varphi =-A/(4B)$. Thus further dynamical information is needed beyond the PDF's unrestricted eccentric orbit to assert unique centres. Exterior asymmetric branches are documented by Winter and Murray, “Resonance and chaos. II. Exterior resonances and asymmetric libration”.

Let $f(\mathbf{p})d^{3}p$ be particle number per unit volume in a momentum element, including the spin-state count in $f$. This is the local momentum part of a phase-space distribution function, normalized by the number density $n=\int f(\mathbf{p})d^{3}p$. For an isotropic distribution with particle energy $E(p)=\sqrt{m^2{c}^4+p^2{c}^2}$, the kinetic pressure of an isotropic gas is the average momentum flux:The kinetic energy density, explicitly excluding rest mass, isFor $p\ll mc$, $E-m{c}^2=p^2/(2m)$ and $\mathbf{v}=\mathbf{p}/m$, giving $P=2u_{\mathrm{kin}}/3$. For $p\gg mc$, $E\simeq pc$ and $v\simeq c$, giving $P\simeq u_{\mathrm{kin}}/3$. ThusNeither relation assumes a Maxwellian distribution. At intermediate momenta neither constant ratio is exact, and anisotropic distributions require a pressure tensor rather than this scalar pressure.

For a classical Maxwell-Boltzmann distribution, $f=n(2\pi m{k}_{B}T{)}^{-3/2}\exp [-p^2/(2m{k}_{B}T)]$. Three Gaussian component integrals give $⟨p^2⟩=3m{k}_{B}T$, so the ideal gas hasIn a fully ionized mixture, sum over independent species to obtain $P_{\mathrm{gas}}=\rho k_{B}T/({\mu }_{\mathrm{mol}}{m}_u)$, where the mean molecular weight counts ions and free Electrons. This excludes partial-ionization and interaction corrections. A dilute classical ultrarelativistic gas still has $P=n{k}_{B}T$ but $u_{\mathrm{kin}}\simeq 3n{k}_{B}T$; its momentum distribution is proportional to $e^{-pc/(k_{B}T)}$ instead of the nonrelativistic Gaussian.

For a thermal photon gas, the two polarization states and zero photon chemical potential give the Planck photon distributionUse $x=pc/(k_{B}T)$ and ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$. The radiation constant and photon equation of state are thereforeLikewise $n_{\gamma }=2\zeta (3)(k_{B}T/\hslash c{)}^3/{\pi }^2$. Unlike a gas with a fixed particle number, photons do not have pressure proportional to baryonic mass density.

For fully degenerate Electrons, the Fermi-Dirac distribution becomes a filled momentum sphere with two spin states. Counting them gives the Fermi momentumHere ${\mu }_e$ is the mean molecular weight per electron. The equation of state of a cold electron gas follows by integrating momentum flux up to $p_F$:Its two limits areThe corresponding kinetic energy densities are $3n_e[p_F^2/(2m_e)]/5$ and $3n_e{p}_{F}c/4$, respectively. Thus in the high-density relativistic limit the Electron pressure scales as ${\rho }^{4/3}$ and is nearly independent of temperature. This is the ideal noninteracting, fixed-composition Electron result, not a universal equation of state at nuclear densities where captures, interactions and the composition change.

The boundaries on a stellar equation-of-state regime diagram concern the Electron component. Define the electron relativistic density threshold and thermal electron relativistic threshold byA degenerate Electron gas is nonrelativistic well below ${\rho }_{\ast }$ and ultrarelativistic well above it: this is an approximately vertical division on a log-density plot. A nondegenerate Electron gas instead becomes thermally relativistic near $T_{\ast }$: this is an approximately horizontal division. The regimes have broad crossovers rather than a discontinuity at either line.

The exact electron Fermi temperature, subtracting Electron rest energy, isFor Electrons with negligible thermal pairs, strong degeneracy requires $T\ll T_F$, whereas $T\gg T_F$ gives a nondegenerate gas. In the pair-rich regime the actual Electron and Positron distributions must instead be treated with their chemical potentials, as discussed below. The two limiting degeneracy boundaries areThese slopes $2/3$ and $1/3$ explain the bent degeneracy boundary on the logarithmic graph. A thermal-wavelength test $n_e[2\pi {\hslash }^2/(m_e{k}_{B}T){]}^{3/2}\sim 2$ gives the same nonrelativistic density/temperature scaling, with an order-one definition of the crossover.

The radiation-to-gas pressure boundary in the nondegenerate fully ionized regime isRadiation dominates above this line. Once Electrons are strongly degenerate, compare radiation pressure with $P_e(\rho )$ instead of continuing the ideal-gas comparison into that regime. The radiation-to-degeneracy pressure boundary is $T=[3P_e(\rho )/a_{\mathrm{rad}}{]}^{1/4}$, with logarithmic slopes $5/12$ and $1/3$ in the nonrelativistic and ultrarelativistic limits. Degeneracy of the Electrons and dominance of their pressure are distinct criteria.

For the pair curve, distinguish baryonic net Electrons from thermally created Electrons and Positrons. Chemical equilibrium with photons requires opposite Electron/Positron chemical potentials when the one-particle energy includes rest energy, as in the distributions below. At low temperature and low degeneracy the zero-chemical-potential density per charge species isCharge neutrality gives $n_{-}-n_{+}=n_b=\rho /({\mu }_e{m}_u)$, while the nondegenerate equilibrium product is $n_{-}{n}_{+}=n_0^2$. HencePairs become important when $n_0$ is comparable with or exceeds $n_b$, approximately the electron-positron thermal pair abundance curve ${\rho }_{\mathrm{pair}}(T)={\mu }_e{m}_u{n}_0(T)$. Below its density at a fixed temperature, pairs dominate over the net charge Electrons. The exponential makes the low-temperature portion steep in a log-log plot. For ${\mu }_e=2$, the full zero-potential integral gives approximate pair-marker densities $1.80\times 10^3\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ at $10^9\mathrm{K}$ and $7.52\times 10^5\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ at $3\times 10^9\mathrm{K}$. At relativistic temperature use the full zero-chemical-potential Fermi-Dirac distribution instead; it gives $n_0\to [3\zeta (3)/(2{\pi }^2)](k_{B}T/\hslash c{)}^3$, so the high-temperature pair curve approaches slope three in log-density versus log-temperature. Strong net-electron degeneracy suppresses Positrons and requires the full chemical-potential-dependent distribution.

The original diagram uses fully ionized helium, ${\mu }_e=2$ and ${\mu }_{\mathrm{mol}}=4/3$, only to set numerical locations. It plots the full $T_F$ curve and the zero-chemical-potential pair integral, rather than extending their asymptotes into the crossover. The pair boundary is deliberately an approximate abundance marker, not a phase transition; at high temperature the nonrelativistic ion approximation and fixed-composition model also have limits. For $T\gg T_{\ast }$ with abundant nondegenerate pairs, the combined Electron/Positron energy density is $7a_{\mathrm{rad}}{T}^4/4$, so photons plus pairs have $u=11a_{\mathrm{rad}}{T}^4/4$ and $P=u/3$. Real stellar matter adds Coulomb effects, partial ionization, nuclear reactions and, at sufficiently high density, nuclear-matter physics beyond the ideal regime map.

Suppress tildes in the normalized equation and let $X=x-vt$. A Gross–Pitaevskii solitary wave stationary in these translating coordinates obeysThe boundary fixes the bulk complex argument and describes a localized disturbance with zero net vortex winding number; a vortex pair is possible. No extra plane-wave factor is needed when only the coordinates are changed. A full Galilean transformation of the field would additionally change the bulk flow and its complex argument convention.

Subtract the uniform background from the grand potential. A convenient dimensionless energy and the renormalized momentum of a condensate areThe subtraction in $\mathbf{p}$ removes the troublesome background complex argument gradient. In particular, for the usual subsonic localized far field $\psi -1=O(r^{-1})$, $\nabla \psi =O(r^{-2})$, and $1-|\psi {|}^2=O(r^{-2})$, the energy density is $O(r^{-4})$ and the momentum integrand is $O(r^{-3})$, both integrable in two dimensions. Merely writing $\psi \to 1$ without sufficient decay is not by itself a convergence proof. The common solitary-wave branch is subsonic: linearizing this normalization gives ${\omega }^2=k^2/2+k^4/4$ and speed of sound $1/\sqrt{2}$.

Integration by parts, with decaying variations and their boundary terms vanishing, givesThus the translating equation is precisely $\delta (E-v{p}_x)=0$. Along a differentiable family ${\psi }_{\lambda }$ with fixed bulk normalization, evaluate this stationarity on ${\partial }_{\lambda }\psi$ to obtainwhere $d{p}_x/d\lambda \ne 0$. This energy–momentum slope of a solitary wave does not require differentiating the Lagrange multiplier $v$ inside the field variation. At a turning point the parametrized identity remains true, while $E$ need not be a single-valued function of $p_x$ locally. The chosen factors make $i{\psi }_t=\delta E/\delta {\psi }^{\ast }$, consistent with the normalized equation; changing the normalization of $E$ alone would change the claimed slope.

The Maxwell equations with vacuum permittivity and permeability give $\rho ={ϵ}_0\nabla \cdot E$ and $J={\mu }_0^{-1}\nabla \times B-{ϵ}_0{\partial }_{t}E$. The current identity is the Ampère-Maxwell equation, and Faraday's law gives ${\partial }_{t}B=-\nabla \times E$. With $g={ϵ}_{0}E\times B$,The identity $(\nabla \times B)\times B=(B\cdot \nabla )B-\frac{1}{2}\nabla B^2$, and its electric analogue, turn this into the divergence of the Maxwell stress tensor $T$. Thus the tensor in the question is its negative:The vector $g$ is electromagnetic momentum density, and ${\sigma }_{ij}$ is the outward flux of its $i$th component across a surface normal in the $j$ direction. In static conditions the total force on matter is $F_i=-\oint {\sigma }_{ij}{n}_{j}dS$; generally subtract the time derivative of the enclosed electromagnetic field momentum as well.

For the two infinite sheets, superposition and the electromagnetic field jumps give $E=(\sigma /{ϵ}_0)e_x$, $B={\mu }_{0}K{e}_y$ between them and zero outside. Substitution yields exactlyA thin pillbox enclosing the $x=0$ sheet has only its $+x$ face in the interior electromagnetic field. ThereforeDirectly, omit the sheet's self-electromagnetic field: the other sheet supplies $E_{\mathrm{other}}=\sigma e_x/(2{ϵ}_0)$ and $B_{\mathrm{other}}={\mu }_{0}K{e}_y/2$. The Lorentz force per area is $\sigma E_{\mathrm{other}}+K{e}_z\times B_{\mathrm{other}}$, giving the same result. The electric contribution is attractive and the magnetic contribution repulsive.

If $K=\sigma v$ for a single moving charge population, use ${\mu }_0{ϵ}_0=c^{-2}$ to obtain $F/A={\sigma }^2(1-v^2/c^2)e_x/(2{ϵ}_0)$. For physical $|v|<c$ this is attractive, so motion of these charges cannot make the force repulsive. Independently imposed currents need not satisfy that restriction on $K/\sigma$.

For covariant quantization of the bosonic string,All other commutators between independent canonical variables vanish. The covariant Fock vacuum obeys ${\alpha }_k^m|0;p⟩=0$ for every $k>0$ and every $m$. Normal ordering moves these annihilating string oscillators to the right and givesThen ${\hat{N}}_{\mathrm{cov}}|0;p⟩=0$. Its commutator $[{\hat{N}}_{\mathrm{cov}},{\alpha }_{-k}^m]=k{\alpha }_{-k}^m$ grades a Fock state basis by $N=\sum k{n}_{km}$, although its inner product is indefinite because of the timelike oscillator.