Use a top-hat plume model per unit span. Define volume and kinematic momentum fluxes and , and set the pure-ice reduced gravity . The suspension buoyancy is . The supplied equations express vertical momentum conservation, conservation of relative enthalpy including latent heat, and conservation of excess salt, respectively. At the release height, prescribe compatible source volume/momentum, salt and enthalpy fluxes, or equivalent source width, velocity, composition and ice fraction; impose ambient entrained water at with no vertical momentum. A pure far-field similarity solution may instead be described using a virtual origin, but its singular near-origin values are not physical inlet conditions.
The missing balance is total mass, equivalently volume under the Boussinesq approximation. In a wall plume, only the exposed outer edge entrains, so the Batchelor entrainment hypothesis gives
If two free edges are intended, replace by everywhere; the numerical coefficient is a geometry convention. The glacier boundary itself has no normal flow. Local phase conversion does not add volume at retained Boussinesq order, and crystals are assumed to share the liquid velocity.
For a box of height , the upward momentum flux changes by . The ambient hydrostatic pressure has been subtracted; the remaining force is the suspension buoyancy . Entrained ambient water carries no vertical momentum, and drag is neglected, giving . Salt entering through the exposed edge is , so and therefore
These control-volume derivations account for entrainment rather than treating as constant.
Thermodynamic elimination. Let the conserved excess enthalpy and salt fluxes be and :
The ambient satisfies , while equilibrium in the ice-bearing plume gives
Substitute this temperature in the enthalpy invariant. Exact rearrangement gives
Thus, at dilute-crystal order,
The ice-bearing wall-plume thermodynamic invariant fixes the leading advected ice volume and buoyancy flux from the source deficits. A rising ice-dominated branch requires . It is not determined from reduced gravity alone. The small- expansion also presumes finite source fluxes and no cancellation making the discarded terms comparable to .
The three remaining leading equations are , , and . Seek a similarity form in with , . Continuity gives , while momentum gives . Therefore . The coefficients obey , giving the self-similar ice-bearing wall plume
The width grows linearly, the speed approaches a constant, and the crystals dilute inversely with height. This is a buoyancy-conserving wall line plume, with buoyancy supplied mainly by the ice fraction rather than a temperature anomaly. The dilute-crystal regime requires and lies away from the formal virtual-origin singularity. If all source flux deficits vanish, and this nontrivial rising branch does not exist. The exact thermodynamic invariants above, rather than their truncated version, are needed near a concentrated source.
For a massive particle of rest mass , its four-momentum is , where the four-velocity is , and . With metric signature , . A photon has four-momentum , a future-pointing null vector. Four-momentum conservation states that the total incoming and outgoing four-momenta agree for an isolated interaction; this includes both energy and momentum conservation in every inertial frame.
Take the incoming photon direction as . Its energy is , so the total initial four-momentum is
The invariant mass available after photon absorption is therefore . In the centre-of-momentum frame, the two equal daughters have opposite momenta. Their total energy is at least their combined rest energy , with equality precisely when both daughters have zero momentum in that frame. Thus the mass threshold after photon absorption is
The bound is attainable kinematically by taking each daughter four-momentum equal to ; each then has the required rest mass. It immediately gives as .
At this threshold both daughters move with the centre-of-momentum frame. The common laboratory velocity follows from :
For , its continuous limiting value is zero. At smaller allowed daughter rest masses, some of the centre-of-momentum energy appears as their relative kinetic energy, so their four-momenta need not be .
First apply momentum conservation to the rocket and the infinitesimal amount of expelled gas. During , the exhaust mass is and its inertial velocity, to first order, is . With no dust, the momentum balance is
It gives the rocket equation , while .
For a dust-collecting rocket, the incoming dust mass is with zero initial momentum. The same balance has rocket mass after the interval, so
Consequently
The second term in the acceleration equation accounts for the momentum required to bring collected dust up to rocket speed.
Set and , so . As long as the engines operate and , the velocity increases and time can be eliminated:
Integrating from , , gives
Thus the required mass–velocity relation is
If , this formula would give as . In this interval, . A real rocket retains positive body mass even after its fuel is gone, and collected dust only increases that residual mass. It therefore cannot follow this branch all the way to zero mass: its fuel must run out while . Indeed the hypothetical time integral is integrable at when , so even the formal zero-mass endpoint would occur at finite time.
If , the same relation instead gives as . The mass first decreases, reaches its minimum at , and then increases because dust arrives faster than fuel is expelled. The time integral now diverges at : under indefinitely maintained thrust, approaches asymptotically while stored dust grows without bound. At that limiting speed, thrust balances the dust-loading term , and . This is the physical reason the total rocket mass no longer signals fuel exhaustion. A rocket with a finite initial fuel supply still runs out at finite time and before reaching ; increasing total mass does not create new fuel.
Take an oriented closed string in flat, critical bosonic string theory, with and a mostly-plus target Minkowski metric. After continuation of the worldsheet to Euclidean signature, the Polyakov path integral sums over embeddings and worldsheet metrics, divided by worldsheet diffeomorphisms and Weyl transformations. In a flat target its kinetic action is
Target-time continuation or analytic continuation of external momenta defines the Lorentzian scattering amplitude; a naive real Euclidean Gaussian for timelike would not be convergent.
An external tachyon is represented by the tachyon vertex operator , with conformal weights . The physical integrated vertex has , so . Schematically its tree amplitude is
The Faddeev-Popov determinant from conformal gauge is represented by the worldsheet ghost fields. At tree level the worldsheet is a Riemann sphere, whose unpunctured complex structure has no moduli. Its residual conformal automorphisms are the Möbius transformations, . Fix three insertion points, accompanying their unintegrated vertices by the required ghost factors. The remaining complex insertion positions are integrated over the sphere; equivalently one integrates all positions and divides by the residual conformal group.
The embedding fields are free, with
Their zero-mode integral gives momentum conservation, and their nonzero-mode Gaussian integral gives the Koba-Nielsen factor
This explains the sphere tachyon position integral without needing its evaluation. For four tachyons the resulting Virasoro–Shapiro amplitude displays the exchanged string spectrum directly.
Use all-incoming external momenta and introduce alongside the two printed dimensionless Mandelstam variables. Since and , one obtains . The physical center-of-mass energy squared in the channel is . The channel measures the analogous crossed momentum transfer, with sign determined by the mostly-plus convention. Rewriting the Gamma factors in a symmetric form gives
The Gamma function poles imply, at generic fixed values of the other invariant,
The denominator Gamma factors can remove residues at special intersecting channel kinematics; the statement concerns a generic single-channel limit. An -channel pole occurs when the intermediate momentum satisfies the mass-shell condition for a closed-string state:
The pole at exchanges the ground-state tachyon, the pole at exchanges massless states, and the positive integer poles exchange the infinite massive tower. The -channel interpretation is the crossed version. Factorization means that each residue is a sum of products of couplings to intermediate physical states that couple to the chosen external particles. It need not expose every representation at that mass.
For example, the Gamma function recurrence and Gamma function residue at a nonpositive integer give the dimensionless Virasoro–Shapiro amplitude pole residue
The residue polynomial has degree , consistent with maximum spin in the exchanged level. The amplitude also has the corresponding -channel poles by crossing symmetry.
The massless fields can be treated as target backgrounds rather than separate asymptotic insertions. Write the target metric as , introduce a Kalb–Ramond field , and a dilaton . In conventional Euclidean signs their string nonlinear sigma model action is
Here is the antisymmetric tensor density. The three backgrounds correspond to the graviton, antisymmetric tensor and scalar states at closed-string level one. Expanding the vacuum functional in , then Fourier expanding the backgrounds, produces exactly their integrated string vertex operators. Its functional derivatives therefore generate the amplitudes with massless external strings. The connected vacuum functional organizes connected amplitudes; the spacetime effective action organizes the corresponding vertices after treating massless propagation consistently.
At momenta small compared with , massive string propagators can be expanded in powers of momenta over their masses. The analytic part of the amplitudes consequently determines local higher-derivative interactions, ordered by powers of . Massless exchange poles are retained through propagation of the massless fields, rather than expanded into local contact terms. Up to field redefinitions, the leading massless-sector action in the string-frame metric is
The terms denoted contain additional derivatives, including curvature-squared terms in the bosonic theory. The expansion concerns the massless sector around the perturbative bosonic background; the tachyon instability remains and is not cured by omitting its field from this displayed action. Thus this is a formal perturbative effective description, not a claim of a stable bosonic vacuum.
Finally split the dilaton into a constant and its variation, . By the Gauss-Bonnet theorem, its dilaton Euler-characteristic weighting on a connected closed oriented surface follows from
Define the string coupling by . A genus- path integral is weighted by
The sphere carries , the torus carries , and each extra handle adds . At higher genus one integrates over complex-structure moduli as well as insertion points, with the associated antighost insertions supplying the correct moduli measure. With canonically normalized external vertices an -point genus- amplitude scales as .
Summing connected worldsheets of every genus of a surface yields the string-loop effective action expansion
Each coefficient has its own low-energy expansion. The two parameters have different roles: resolves finite string size through higher derivatives, while counts additional string loops.
All dots in this solution denote differentiation with respect to conformal time . This is required by comoving , the streaming term , and . The PDF's reference to proper-time derivatives is inconsistent with these coefficients. Proper-time equations follow by dividing the conformal-time equations by , replacing streaming by and the collision rate by .
The monopole equation of the photon Boltzmann hierarchy is ; its collision term vanishes. Consequently the photon continuity equation is
For the dipole the hierarchy gives
Using the relations from the preceding part, the photon Euler equation becomes
Because , the collision term drives the photon velocity towards the baryon velocity, as Thomson scattering should.
For pressureless baryons, stress-energy conservation without energy exchange gives . At this order Thomson scattering changes momentum but not the local photon energy: energy transfer by the relative bulk motion is second order, and recoil/thermal energy exchange is neglected. Total momentum conservation fixes the baryon force to be the negative of the photon force. The momentum densities are proportional to their enthalpy densities, and . Thus
Including the pressureless expansion and gravitational terms yields the baryon Euler equation with Thomson drag
The expansion coefficient for photons cancels against the changing background pressure, while pressureless baryons retain the Hubble drag term. The opposite drag signs and the enthalpy ratio ensure that the combined photon-baryon fluid conserves momentum.
Use flavor hypercharge , where is baryon number and is strangeness. The isospin coordinate is . The baryon octet has coordinates
Here the nucleons have , the Sigma baryons and Lambda baryon have , and the Xi baryons have . The central Sigma baryon belongs to an isospin triplet while the central Lambda baryon is an isospin singlet; equal coordinates do not identify the states.
The pseudoscalar meson octet is
The pions have , the upper kaons have , and their lower antiparticles have . The two central states are the neutral pion and the Eta octet state. This is the octet basis of flavor symmetry; the physical eta can also mix with the flavor-singlet state.
Figure 1.
Flavor SU(3) baryon and pseudoscalar meson octets in isospin and strong hypercharge coordinates
. Details.
For the flavor SU(3) Cartan generators, choose the Hermitian physics convention and the inner product . Then an orthonormal diagonal basis is
Strictly, are elements of the anti-Hermitian SU(3) Lie algebra; are the corresponding Hermitian observables. In the quark basis , the up and down quarks form an isospin doublet and the strange quark is a singlet. Their values are . Each quark has baryon number , and their strangeness values are . It follows that
These flavor hypercharge conventions differ from the electroweak hypercharge convention. If instead the inner product is , the orthonormal basis is , , and the same operators are , .
In the ordinary quark model, the proton has valence content and charge , while the neutron has content and charge zero. Additivity of electric charge gives and , hence , . The Sigma baryon has content and charge , giving . Thus the quark triplet's electric charges, in units of the positive elementary charge, are
The Gell-Mann--Nishijima formula is consequently
It also reproduces every baryon octet charge from the first diagram. On antiquarks the additive quantum numbers reverse sign, and combining a quark with an antiquark reproduces the meson octet charges.
Because the down and strange quarks have identical electric charge, commutes with the U-spin generators
They satisfy ; equivalently the anti-Hermitian matrices span an subalgebra. The entries on the block of are equal, so for all three U-spin generators. The electric charge is therefore constant within each irreducible U-spin multiplet. For example, U-spin relates and , and relates and , without changing their charge. It does not imply exact mass degeneracy: unequal down- and strange-quark masses break U-spin.
For pion-nucleon octet channels, assume the collision is governed by the strong interaction. The initial baryon number is one and strangeness is zero, so an outgoing meson-baryon pair must preserve , , and electric charge. Thus its total flavor hypercharge is . The allowed types are
A kaon of can accompany a Lambda baryon or Sigma baryon of . An antikaon cannot balance the nonpositive strangeness of an octet baryon. A Xi baryon would require a meson of , which the meson octet does not contain.
Resolving these types by electric charge gives all possible pairs:
In the isospin-symmetric approximation, total isospin is conserved as well: the incoming contains . The and channels contain both values, while and contain only . The extreme-charge initial states are pure , consistently excluding and . Clebsch-Gordan coefficients relate amplitudes in different charge channels; the table establishes permission, not equal probabilities. Electromagnetism and unequal up- and down-quark masses introduce small violations of isospin symmetry.
Finally, energy and momentum conservation require for a particular pair, where is the squared total four-momentum. Only channels above their own threshold can occur. Total angular momentum and parity symmetry in quantum field theory constrain the partial waves of the meson-baryon scattering: a pseudoscalar meson and a positive-parity spin-one-half baryon have pair parity and total angular momentum (only for ). Initial and final partial waves must have matching and parity. Sufficient energy can open a channel, but cannot remove the electric charge conservation, baryon number, or strangeness constraints.
Let be longitudinal displacements from equilibrium. The original diagram has masses , and no attachments to fixed walls. The Lagrangian is
The Euler-Lagrange equations give
The forces sum to zero, consistent with momentum conservation for an isolated molecule.
Work in one fixed inertial frame, and take positive velocity in the rocket's direction. During , the rocket ejects the positive amount of mass . Its exhaust has inertial velocity , by the Galilean transformation of velocities. Once expelled, that material retains its velocity because it experiences no external force. The system's momentum is therefore the rocket momentum plus the sum of all these exhaust momenta:
The minus sign comes from the positive expelled mass being .
Differentiate, applying the product rule to the rocket momentum and the fundamental theorem of calculus to the integral:
The momentum conservation law thus gives the nonrelativistic rocket equation
The same bookkeeping for kinetic energy gives
Its derivative is
where the rocket equation cancels the first two terms. Thus kinetic energy increases during mass ejection, with when and . This does not violate conservation of energy: the increase comes from stored internal energy.
Now denotes the instantaneous rest mass of the rocket, including its fuel and stored internal energy, and is measured in a fixed inertial frame. Assume the exhaust is expelled backwards with speed in the rocket's instantaneous rest frame. By the relativistic velocity-addition formula, material emitted at time has inertial velocity
It coasts at this velocity after emission.
The crucial change is that the expelled material's rest mass is not simply : some of the rocket's rest energy becomes exhaust motion. Instead, use conservation of relativistic energy. Put , so the rocket's energy is . The exhaust emitted during carries energy . For material moving at velocity , the relativistic energy-momentum relation gives momentum , since and . Hence the modified total momentum is
Meanwhile the total relativistic energy is . These expressions consistently account for both energy and momentum of the expelled material.
Differentiate the momentum formula. The momentum conservation law gives
To simplify it, use
Multiplying by , the equation becomes
The two terms containing cancel. Dividing by proves the relativistic rocket equation
Define the open-string endpoint momentum flux along the open string by . Variation of the Nambu-Goto phase-space action gives
The two last equations are the Virasoro constraints. Integration of the momentum equation gives
Thus the necessary and sufficient condition for total momentum conservation is equality of the two endpoint fluxes, component by component. Vanishing of both fluxes is a sufficient local boundary condition.
The spatial boundary term in the action variation is . For free ends the endpoint variations are arbitrary, so
In a gauge with and nonzero this reduces to the usual Neumann boundary condition . The conclusion is conservation of total momentum, not that total momentum must vanish. The PDF contains the derivative in this conclusion; the TeX transcription drops it.
For endpoints confined to , the Dirichlet boundary condition fixes . It therefore imposes no requirement that vanish. Preserving the fixed position requires at each endpoint, but that is a different condition. For example, in gauge the endpoint momentum is zero while can be nonzero. Consequently is generally not conserved: the normal endpoint forces transfer momentum to the hyperplane. A dynamical D-brane recoils and carries the missing momentum. If the hyperplane is idealized as fixed, its external support absorbs the momentum. The combined string-and-support system conserves momentum; freely varying directions tangent to the hyperplane retain their vanishing-flux condition.
Use the Fourier transform convention . Statistical homogeneity and isotropy define the primordial bispectrum by
The Dirac delta function enforces momentum conservation; the three magnitudes describe a closed triangle. Set . For the centered quadratic Gaussian transformation, its quadratic Fourier term is
At first order in , choose the quadratic field at one of the three external positions. For the first position, Wick contractions pair the two internal fields with the two remaining external fields in two ways. The self-pairing is precisely removed by the subtracted mean. The surviving contribution is times the momentum delta. Adding the other positions gives the tree-level local-type primordial non-Gaussianity result
This expression is at leading order in . The exact quadratic map also has a loop correction to the local primordial bispectrum, from three quadratic vertices:
There is no quadratic-in- bispectrum term, because it would involve five centered Gaussian fields. A regulator may be needed for this higher-order integral. The displayed leading formula uses the Gaussian power , consistently neglecting its order- correction.
Write the total remaining rocket mass as , with . Over an interval , an exhaust mass is expelled backwards with inertial velocity to first order. Applying momentum conservation to the rocket and this exhaust gives
Cancellation and division by give the rocket equation
With , , and , integration gives
At burnout,
The logarithm depends on the mass ratio, not on the burn rate. This follows from momentum conservation of the complete rocket–exhaust system; applying a constant-mass equation directly to the rocket alone would omit the momentum carried away.
Use momentum conservation and conservation of energy in the laboratory inertial frame. The photon has momentum magnitude , while the outgoing Electron satisfies the relativistic energy-momentum relation. Its total energy and momentum obey
Substitute these into . Cancelling and gives . Both photon energies are positive, so division yields
This is Compton scattering. The energies include the Electron's rest energy; using only its kinetic energy in the mass-shell equation would give the wrong result. The outgoing photon has unchanged energy at and smaller energy for a nonzero scattering angle.
Take Euclidean gamma matrices with , and independent Grassmann fields . The minimally coupled Dirac action is
Its mass term is the Dirac mass term, and the interaction is . With Fourier convention , the free kinetic operator is . The Euclidean adjoint Dirac Feynman rules are therefore
The vertex sign follows by expanding , not . Momentum conservation supplies the usual Dirac delta function at each vertex. There is only a one-gluon fermion vertex, with no two-gluon seagull because the fermionic covariant derivative is first order. A closed fermion loop has the usual fermionic sign . The identity checks the Dirac propagator.
Use a two-sided top-hat line plume with full width , upward speed , and reduced gravity uniform across the plume. Fluxes are measured per unit span along the line source: the volume flux is , the kinematic momentum flux is , and the buoyancy flux is . If the inward edge speed is , where is the entrainment coefficient, the two exposed edges give
These are the volume conservation, momentum conservation, and buoyancy flux balances for a Boussinesq approximation plume in an unstratified ambient. Entrained ambient fluid supplies neither vertical momentum nor reference buoyancy.
A pure plume has no persistent source length scale. Since a line-source has dimensions , dimensional analysis gives constant , , and . Substituting constant into and the momentum balance determines the coefficients:
Thus and when width is full width and velocity and buoyancy have top-hat profiles. If width means half-width, instead. Other prescribed profile shapes change these numerical factors; the linear growth laws remain the same. The point-source in Question 1 has different dimensions, so its law must not be used here.
Let be local contaminant concentration and let be horizontally integrated plume concentration. Its conserved amount per unit span is . The dilute contaminant is a passive scalar; the maintained plume is unaffected by its impulsive release. A one-dimensional effective transport closure takes the integrated scalar flux to be
The constant advective speed and the eddy diffusivity scaling follow from the line-plume scales. The entrainment coefficient alone does not determine the scalar dispersion coefficient: this Fickian closure for the integrated variable is an additional modelling assumption. Mass conservation then gives the displayed transport equation in the PDF. With this effective closure and advective speed chosen as , ; if , then , with independent dimensionless mixing coefficient .
There is a second possible closure convention. If one instead applies local diffusive flux to horizontally averaged plume concentration , its integrated diffusive flux is . For , the total scalar flux is then . The same printed equation is obtained with , rather than . Thus the printed and should be regarded as effective coefficients unless the averaging and closure conventions are specified; no equality between and the velocity prefactor is universal.
Assume , , an initial impulse at the origin, no further contaminant input, and zero endpoint scalar flux for . For a similarity solution, put
Substituting into the advection-diffusion equation gives
Decay at infinity makes the integrated constant zero. Hence , and normalization with the Gamma function gives . The resulting gamma impulse solution for linearly increasing diffusivity is
Its integral over is . Its flux is , which vanishes at both endpoints for . The scale shrinks to zero as , so the normalized solution has weak convergence of probability measures to the required unit source impulse. The mass at the boundary is a full unit impulse on the half-line, not half of a whole-line impulse.
For the integrated quantity, , and consequently
This verifies the printed location for horizontally integrated plume concentration. It is a maximum because the derivative changes from positive to negative there. The normalized profile is a gamma distribution with shape : its expected value is and its variance is , so neither the mean position nor the spread should be mistaken for the modal position.
The PDF changes from “integral” to “averaged” concentration in its final request. A genuine horizontally averaged plume concentration divides by the growing width , and is therefore
It is normalized by , not by . When , its interior maximum is
When it decreases from a finite boundary supremum; when it is singular at the ideal point source and has no positive interior maximum. A finite source regularizes that singularity. Thus the final printed maximum is correct for the integrated variable , or for an “average” using a fixed reference width, but is not generally the maximum of the local mean . Both quantities have been given rather than silently identifying them. If , the smooth similarity formula is replaced by the advected impulse ; positive eddy diffusivity is essential to the gamma profile.
Use and the physical real part of the printed plane wave. Its null wavevector obeys , so each component satisfies . Direct differentiation of the defined stress tensor gives
The trace is proportional to for a single plane wave and hence vanishes. Therefore
This is local conservation of energy and momentum conservation; spatial integration gives conserved total charges when the boundary flux vanishes. For a finite-amplitude complex phasor, one takes the real field before quadratic products and may time-average: in transverse-traceless gauge. Its energy density is nonnegative and its flux travels in the positive direction.
Let the charge parameter depend on position temporarily. Using the gauge-covariant derivative of a charged scalar field, the variation of the matter action is . The resulting Noether current is
Here . The potential does not contribute, and . Because has not been rescaled to a canonical kinetic term, there is no factor in this Noether current or in the matter interaction vertices.
A local change of variables in the path integral gives . Integration by parts in the term containing produces the two opposite contact terms of the Ward identity.
For clarity about phases, transform the current correlator with and suppress its overall momentum conservation Dirac delta function. Denote the result by . The coordinate Ward identity becomes
Amputation of the two external quantum field theory propagators then fixes the longitudinal part of the current interaction vertex. Since the action couples to , the gauge interaction vertex has the opposite sign. The phase convention for the paper's gauge interaction vertices is , where is the corresponding derivative of the Euclidean quantum effective action.
One can derive the exact one-particle-irreducible correlation function directly, without assuming the entire connected current correlator is one-particle irreducible. Local gauge invariance of the quantum effective action implies
A linear covariant gauge fixing adds a breaking term independent of the scalars, which disappears after the scalar differentiations below. The identity requires a regularization in quantum field theory and counterterms preserving the Abelian gauge theory Ward identity; scalar quantum electrodynamics has no gauge anomaly.
Differentiate with respect to and , then set all background fields to zero. The two-point derivative is the inverse exact quantum field theory propagator. Use the Fourier transform convention , , with the photon momentum incoming. This gives
At tree level , so , and the identity reduces to .
To obtain the two-photon scalar Ward identity, differentiate the same functional Ward identity also with respect to . In terms of Euclidean action derivatives it reads
This construction automatically includes the seagull vertex: the Noether current depends on , with . It cannot be discarded when differentiating a current insertion.
Take the first photon momentum to be incoming , and the second incoming . On the first scalar leg the contact term shifts to ; on the other it shifts to . Consequently
With , the result is
There is a momentum-routing sign error in the printed second identity. With the scalar momentum convention of the first identity, its first shifted argument must be , not . This is already forced at tree level: the seagull vertex has and hence . The corrected right-hand side is , matching the left-hand side, whereas the printed right-hand side is . Changing to an all-incoming scalar convention would also change the first identity, so it does not fix both printed formulas simultaneously.
Initially every particle has , so its Jacobi energy in a shearing sheet is . During an inelastic collision, positions are fixed at the instant of impact, while momentum conservation preserves the sum of the tangential velocities. Hence the sum of , and therefore the sum of the epicyclic guiding center positions , is unchanged. Between collisions, these quantities are individually conserved.
The collision dissipates kinetic energy without changing the instantaneous tidal potential. Consequently the total Jacobi energy in a shearing sheet decreases. At any later time, writing for the particles' current epicyclic guiding center positions gives
If is the accumulated energy dissipated in the inelastic collisions, comparison with the initial circular orbits gives
Thus the mean guiding-center position stays fixed, while its variance grows. The initial ensemble has no preference for positive or negative , and the local equations and collision law preserve the symmetry . The spreading is therefore symmetric in the ensemble average: angular momentum transport moves some particles inward and others outward. A particular finite random realization need not be exactly symmetric. This is the microscopic energy argument for dissipative spreading of a planetary ring.
Let and , so mass conservation for the two-dimensional incompressible flow is automatic. Write for buoyancy and introduce the diffusive operators
The Linearized Boussinesq equations are
Taking the curl of the momentum conservation equations eliminates the pressure and gives
Apply and use . Since the constant-coefficient linear partial differential operators commute,
This is the viscous-diffusive internal gravity wave equation.