Write . In the uniform equilibrium, the linearized ideal magnetohydrodynamic equations reduce to
Integrating from the undisplaced reference state, and using the uniform background magnetic field, gives
These are the perturbations induced by the fluid displacement; independent time-independent changes to the reference state are excluded. The linear Lorentz force density is . Its magnetic pressure and magnetic tension parts give
For a Fourier mode, put , and . Then and . Substitution gives
If , the fluid displacement is perpendicular to both the wave vector and the magnetic field. Only magnetic tension restores it, and
This is the Alfvén wave. The Alfvén velocity is the field-directed propagation vector: . The signed phase velocity normal to the wavefront is ; the group velocity is . The distinction matters for oblique propagation. For generic directions the Alfvén wave has one transverse polarization, with .
To obtain the other magnetohydrodynamic waves, take the scalar products with and :
The determinant condition for a nonzero pair is
Thus the fast magnetosonic wave and slow magnetosonic wave have
Their fluid displacements lie in the plane of and and are generally compressive. In the fast magnetosonic wave, gas and magnetic pressure provide the stronger restoring combination; in the slow magnetosonic wave, their perturbations oppose one another. The phase speeds depend on direction. They satisfy and . For parallel propagation the two speeds are and , with a degeneracy between a transverse branch and the Alfvén wave. For perpendicular propagation and ; the latter is a nonpropagating limiting disturbance. These special directions require interpreting the polarizations by continuity rather than assuming three distinct nonzero frequencies.
In a steady flow the ideal magnetohydrodynamic induction equation gives . Since both vectors have only poloidal magnetic field components, write . Axisymmetry implies
For a field regular on the axis, ; equivalently one may impose zero toroidal electromotive force. Under this physical condition, is parallel to . Write . Continuity and then give . The poloidal magnetic flux function has , so locally on each connected magnetic surface,
This regular-axis alignment of steady poloidal ideal flow needs its zero-circulation condition: axisymmetry alone does not imply this alignment on a domain excluding the axis. An explicit counterexample is available in a cylindrical annulus. Take nonzero constants , positive , and
Here , and has zero curl. The flow satisfies continuity, , and . Its acceleration and Lorentz force density vanish, so it solves the steady equations while and are perpendicular. The nonzero toroidal circulation and the axial singularity explain why this counterexample is excluded from a regular polar accretion column. The remaining derivation uses the regular, aligned branch.
For a narrow stream tube, integrating continuity gives constant . On a nonzero-flow tube the magnetohydrodynamic mass loading is constant, so . Therefore is constant and . This also follows directly from conserved magnetic flux through a flux tube.
For constant isothermal sound speed, , with arbitrary reference density . Project momentum along . The Lorentz force density has no component in that direction, and alignment implies . Thus the isothermal magnetic Bernoulli integral is
Changing only shifts the Bernoulli function by a constant. Since is constant along the magnetic field,
Differentiating the isothermal magnetic Bernoulli integral therefore gives
In the slender polar accretion column, the dipolar flux-tube area is . The decreasing axial field is a leading approximation, not an exactly solenoidal field throughout a cylinder. Indeed for that literal field. Near the axis a small radial component supplies the required radial divergence. Its magnitude is smaller by , even though its divergence is leading order. Keeping this expanding flux tube while neglecting transverse forces yields the intended one-dimensional model.
Use the positive inward speed ; the signed sonic velocity is . Conservation of mass and the Bernoulli equation give
Eliminating yields the isothermal dipolar accretion equation,
A smooth transonic branch requires both sides to vanish at its sonic point, so
This is the critical point of the flow equation; a generic subsonic solution need not cross it. A crossing in the exterior column requires , and a sonic point strictly outside the star requires . Differentiating the equation at the crossing gives . For inward accretion that accelerates toward the star, .
The reservoir boundary condition fixes . At the sonic point, and , giving
Finally use the sonic mass flux through the dipolar flux-tube area:
This rate belongs to the smooth transonic branch of the idealized column, with the supplied total loaded area. The reservoir condition alone does not force every steady solution onto this branch. A real narrow dipolar column must match an outer flow where the slender approximation ceases to apply.
Use the ideal long-solenoid approximation: negligible pitch and end effects, vacuum permeability , and negligible displacement current. By cylindrical symmetry the magnetic field is axial. An Ampère's law rectangular loop with one axial side inside and the other outside gives . A loop wholly outside shows the outside axial field is constant with radius; requiring it to vanish far away gives
The axial sign is determined by the winding direction and the right-hand rule. The small axial component of current from the spiral pitch is neglected in this idealization.
There are turns, and the magnetic flux through each is . The flux linkage is , so the self-inductance is
The induced electromotive force opposes increasing current by Faraday's law. Together with Ohm's law, the circuit equation is . Solving this first-order ordinary differential equation with gives
The time constant is , and the limiting current is .
The magnetic field satisfies . Expanding the divergence of its Maxwell stress tensor gives
For the Newtonian gravitational field , Poisson equation for Newtonian gravity gives , and the gradient representation gives . Consequently
The negative of the Newtonian gravitational stress tensor, in this force-stress convention, therefore supplies . Adding the pressure stress supplies , so
Each summand is symmetric. Absence of external gravitational sources is needed to represent the full gravitational force by this self-gravitating stress.
Integrate Poisson equation for Newtonian gravity through a narrow slab around the disk. The horizontal derivative contributions vanish as its thickness tends to zero, leaving the normal-derivative jump
The even function symmetry of makes the derivatives opposite, so
In the current-free simply connected upper half-space, Ampère's circuital law gives and permits a magnetic scalar potential. Rescale it so that . The divergence-free condition makes satisfy Laplace's equation, with
Compare with the gravitational jump condition. Subject to the same isolated-field boundary condition at infinity, is the harmonic potential of the effective surface density
This gravity-equivalent magnetic surface density can have either sign; it is a mathematical representation of the exterior magnetic field, not physical negative mass. An imposed nondecaying field would require additional boundary data and would not be fixed by the disk surface density alone.
The Lorentz force in the uniform magnetic field gives
where the last step is the vector triple product identity. Hence
On the other hand, , so
The derivatives cancel, proving the magnetic angular-momentum invariant
The magnetic Lorentz force is perpendicular to velocity, so its power is zero:
Thus the kinetic energy is constant.
For the alternative kinetic energy expression, work on an interval where so that the radial unit vector is defined. Differentiating once and twice gives
Since , it follows that and . Therefore the unit-direction kinetic-energy identity is
This identity holds for any twice differentiable trajectory away from the origin, not only for motion in a magnetic field.
The solenoidal constraint on the magnetic field reduces to . Locally, or globally in a simply connected cross-section, it permits a Cartesian magnetic flux function satisfying
The transverse magnetic field lines are contours of , because . The independent component supplies the twist of the flux tube; it is not restricted by the solenoidal constraint.
For the Lorentz force, take the curl explicitly:
Its cross product with has transverse components and component . Consequently
The term is the transverse gradient of the magnetic pressure associated with the axial field, while the other terms include magnetic tension.
Use , , and a prime for . For the plane-parallel magnetohydrodynamic flow with imposed shear, , while
The solenoidal magnetic field constraint gives . The component of the MHD induction equation then gives . Thus is constant in both space and time.
Use the magnetic tension and magnetic pressure decomposition
The continuity equation and momentum equation now reduce to
The transverse components of the MHD induction equation are
The terms involving respectively accelerate the flow across the imposed shear flow and wind the transverse magnetic field. The closure is the isothermal equation of state , as printed in the PDF.
The logarithmic density term is the barotropic fluid energy for , rather than the thermodynamic heat content of an isolated gas. Define
For any scalar , the continuity equation implies . Therefore , , and dotting the momentum equations with gives
The MHD induction equation gives the complementary magnetic energy balance
Adding these two identities, and collecting the product derivatives, yields the barotropic magnetic energy equation
Here . The two stress terms in describe magnetohydrodynamic shear work: the maintained background shear flow can supply energy to, or remove energy from, the perturbation flow and magnetic field. The sign depends on the off-diagonal total stress; this is why the energy excluding the background shear flow is not generally conserved. Changing merely adds a constant multiple of the conserved density to and the corresponding mass flux to , leaving unchanged.
For a steady flow, mass conservation gives and . Define the Alfvén velocity and . Combining the horizontal momentum and MHD induction equations without dividing by gives
Dot these identities with and , respectively:
The vertical momentum equation and the isothermal equation of state give
Multiply by and eliminate the magnetic derivative to obtain
No division by was needed in deriving this necessary relation.
The magnetosonic critical speeds in the direction are
The plus sign gives the fast magnetosonic wave speed and the minus sign the slow magnetosonic wave speed. Thus the differential coefficient is . A smooth outflow proceeding from below both speeds to above both must normally pass through both magnetosonic critical speeds. At each crossing, the right-hand side must also vanish; this is the regularity at a magnetosonic point condition. The derivative coefficient changes sign at each nondegenerate crossing. For and , the driving term must be positive below the slow point, negative between the points and positive above the fast point.
The Alfvén speed component satisfies . With downward gravity , the gravitational term is consequently nonnegative at the slow point and nonpositive at the fast point. The magnetohydrodynamic shear work contribution must balance it at each point, and can provide the upward driving needed to pass the fast point. For and strictly positive , generic separated slow and fast points cannot satisfy the required zero numerator; special vanishing-gravity or coincident-speed cases need separate treatment.
There is also Alfvén-point compatibility in a plane-parallel sheared flow. At , the original transverse equations require
These restrictions are not generally visible as a zero of the scalar differential coefficient, which there equals . In particular, a strictly accelerating regular solution must have at that point. The scalar relation is therefore a necessary wind equation, not a substitute for regularity of all the original ideal magnetohydrodynamic equations. Degenerate cases such as a purely longitudinal magnetic field can merge characteristic speeds and reduce the number of distinct critical conditions.
Use the central-field identity and substitute the prescribed relation between the toroidal magnetic field and axial magnetic field:
The flux-tube axial field virial identity bounds the last average by , while . Therefore the uniform-twist central axial field bound is
The stronger exact expression is .
In SI units, magnetostatic equilibrium without gravity is . For a purely axial magnetic field, the magnetic tension term vanishes. Thus the radial balance is
Matching the magnetohydrodynamic total pressure to the field-free exterior gives the gas pressure and axial field bound:
The inequality follows directly from nonnegative gas pressure. It is the magnetic pressure limit for an untwisted flux tube, without inward hoop tension.
The azimuthal derivative of a cylindrical vector Fourier mode must include basis rotation: the azimuthal dependence refers to the components in the rotating cylindrical basis. Since and , a vector amplitude obeys
Also in Cartesian coordinates, so . Consequently
and the basis-rotation terms cancel in the ideal magnetohydrodynamic induction equation:
Here the printed is the coefficient in the specified magnetic field: . If denotes SI electric current density, then . The field and all subsequent printed coefficients are mutually consistent with this normalization.
For a normal mode with , induction gives . Substitution into the momentum equation, including the Coriolis force, gives the reduced velocity equation
Thus the coefficients are the printed and . Let . Inserting the allowed relation and multiplying by yields
Hence the uniform-current rotating magnetohydrodynamic wave dispersion is
The radicand is real. For , because , so both roots are purely imaginary. For the roots of the quadratic are also imaginary; any zero-frequency case must be checked in the original equations, since the elimination divided by . A positive real part requires
For integer this is possible only when and . Taking the conventional representative gives the printed exception , with and sufficiently large . Literally allowing negative integers also gives , ; complex conjugation maps to , so these describe the conjugate real disturbance. Purely imaginary roots describe neutral oscillatory modes; a repeated root at threshold does not itself prove boundedness of arbitrary initial perturbations.
Use perturbations for velocity, magnetic field, magnetohydrodynamic total pressure and buoyancy displacement variable. Set
This is advection by the background Keplerian shearing sheet. Subtracting the equilibrium before retaining first-order terms gives the nine linearized equations:
The azimuthal coefficient combines the Coriolis acceleration with perturbation advection of the background shear. The term is the winding of radial magnetic field by that shear. The tidal force cancels when the equations at a fixed position are subtracted; it does not add a separate Eulerian perturbation force. Magnetic pressure is already included in , leaving only linear magnetic tension. In this normalization has dimensions of length because .
When the imposed magnetic field vanishes, the dynamical dispersion relation becomes
The nonstationary hydrodynamic branch therefore has . The radial Solberg–Høiland instability criterion in this local model is
The radial buoyancy frequency supplies either restoration or a driving force, while the radial epicyclic frequency supplies rotational restoration. Zero sum is marginal, and positive sum gives stable oscillations. The neutral roots do not themselves signify exponential growth.
Because both equilibrium pressure and dimensionless specific entropy decrease outward, their radial gradients have the same sign. The leading minus sign in the expression for therefore makes : the radial stratification is adverse. For smooth profiles varying over radius , and , so
For a thin disk, the magnitude of this negative radial buoyancy frequency squared is much less than . Rotation stabilizes the hydrodynamic mode despite the adverse entropy gradient. This estimate assumes gradients on the global radial scale; a sharp thermal feature may be different. It addresses the ideal nondiffusive equations here. Processes such as convective overstability require additional thermal relaxation and are not ruled out by this particular criterion.
A cataclysmic variable is a close semidetached binary in which a white dwarf accretes from a cool, usually low-mass donor star filling its Roche lobe. In an ordinary hydrogen-rich system the donor star is often near the lower main sequence. In a nonmagnetic cataclysmic variable (a disk-fed CV), the white dwarf's magnetic field is too weak to control the flow over the disk. This leaves a characteristic disk-fed geometry, summarized in the original schematic below. NASA's introduction to cataclysmic variables describes the basic components.
Roche-lobe overflow passes through the inner Lagrange point . A nearly ballistic stream bends in the rotating frame and strikes the outer accretion disk, producing an accretion-disk stream-impact hot spot. Its retained angular momentum prevents direct radial infall. Viscous evolution of an accretion disk transports angular momentum outward while gas moves inward through a nearly Keplerian accretion disk. Close to the white dwarf, gas slows from orbital rotation toward the stellar rotation in an accretion-disk boundary layer. The disk, accretion-disk stream-impact hot spot, accretion-disk boundary layer, white dwarf and donor star all contribute to the spectrum and, where the orbital inclination permits, to the eclipses of an eclipsing binary. In particular, “nonmagnetic” describes the accretor's control of the flow; it does not imply that the donor star cannot sustain a magnetic field.
The gravitational power available from accretion is approximately . For a slowly rotating white dwarf and a thin steady Keplerian accretion disk, the specific energy changes from approximately zero far out to at the inner disk. Roughly half the available power is radiated by the accretion disk; the rest is released in the accretion-disk boundary layer as the orbital kinetic energy is dissipated. Stellar rotation and departures from a steady thin disk change this partition.
A classical nova has a different energy source. Transferred hydrogen accumulates on the white dwarf; compression heats the base of its envelope until temperature-sensitive hydrogen burning accelerates. Electron degeneracy pressure initially weakens the expansion response to heating, helping a thermonuclear runaway develop. The envelope subsequently expands and ejects material, producing a large optical outburst followed by a decline as the ejecta expand and residual burning ends. The white dwarf usually survives, so continued accretion can build another fuel layer. A rough recurrence scale is the ignition-envelope mass divided by the mean accretion rate; both this mass and the rate vary strongly among systems. The event is an envelope eruption, and the retained fraction is not automatically unity. Starrfield, Iliadis and Hix's nova calculations explains this nuclear mechanism. A classical nova is powered by unstable nuclear burning on the white dwarf.
A dwarf nova undergoes recurrent, shorter brightenings powered principally by enhanced gravitational accretion. The hydrogen-ionization disk instability creates cold, mostly neutral and hot, ionized branches of the disk's accretion-disk thermal S-curve, separated by unstable equilibria. In quiescence the cool disk stores matter because inward transport is slow. Once a critical surface density is reached, a heating transition puts the disk into a hotter, more state with higher effective viscosity: the inward accretion rate and luminosity rise and the disk drains. A cooling transition returns it to quiescence, completing the cycle. A persistent increase in the donor star's transfer rate is not required. Sufficiently high transfer rates can keep the disk on its hot stable branch, giving a nova-like variable rather than ordinary disk cycles. Lasota's disk-instability analysis and Hameury's disk-instability review develop this picture. A dwarf-nova outburst is a disk instability, not a white-dwarf thermonuclear explosion. The names classify mechanisms and need not identify permanently distinct binaries: a nova-producing binary can also possess an unstable disk between nuclear eruptions.
The cataclysmic-variable orbital-period distribution is not smooth. For ordinary hydrogen-rich cataclysmic variables, prominent features are a cataclysmic-variable period gap around two to three hours, a cataclysmic-variable period minimum near eighty minutes, and an accumulation near that minimum. These are population features rather than absolute exclusions. Selection effects matter: luminous high-accretion rate systems are easier to find than faint evolved systems. Helium-transferring binaries have a different period range and are not described by the hydrogen-rich minimum. Gänsicke and collaborators' period-minimum study documents the observed accumulation.
The Roche-lobe-filling period-density relation makes the orbital period a measure of donor structure. Combining with Kepler's third law gives
As the donor star loses mass, its stellar radius response exponent implies
A donor with evolves toward shorter orbital periods. When its effective response falls below , continued mass loss instead lengthens the orbital period: this is the cataclysmic-variable period bounce. A very low-mass donor may be substellar and increasingly affected by electron degeneracy pressure; the ideal degenerate scaling illustrates the reversal. The precise cataclysmic-variable period minimum depends on thermal disequilibrium and the strength of orbital angular momentum loss. Near a turning point is small; for an approximately steady evolutionary flow of systems, the number per period interval scales as , explaining the accumulation. Knigge, Baraffe and Patterson's donor-based evolutionary study relates the donor star's response to these period features.
Long-term transfer is driven mainly by losses of orbital angular momentum, rather than by disk outbursts. Magnetic braking of a binary star removes the cool donor star's spin through a magnetized stellar wind. Tidal synchronization makes the orbit replenish that spin, so the wind extracts orbital angular momentum. This is usually the dominant standard driving mechanism above the cataclysmic-variable period gap. Gravitational-wave emission from a binary system supplies a baseline loss, especially important below the gap. For a weak-field, slowly moving circular orbit, the circular gravitational-wave inspiral gives
Because ordinary transfer from the lighter donor star tends to expand its Roche lobe if orbital angular momentum is conserved, an external loss is needed to sustain contact. In the conservative contact approximation, the binary mass-transfer contact equation reads
where the stellar radius response exponent must match the evolutionary timescale and additional donor expansion has been neglected. The positive denominator on the stable branch makes drive . Nova ejecta can introduce additional nonconservative binary mass transfer. Knigge's evolutionary discussion describes the standard loss mechanisms and their limitations.
In the disrupted magnetic braking model, relatively rapid mass loss above the cataclysmic-variable period gap keeps the donor star inflated relative to stellar thermal equilibrium. When the donor approaches the fully convective star transition, the model postulates a substantial reduction in magnetic braking of a binary star. The donor star can contract within its Roche lobe, suppressing Roche-lobe overflow near the upper edge of the gap. Gravitational-wave emission from a binary system continues to shrink the detached binary; near the lower edge the Roche lobe again reaches the donor radius and transfer resumes. No mass transfer is needed during the detached crossing. The density relation predicts a radius ratio between contact at three and two hours if the masses remain nearly fixed, illustrating the required inflation before detachment. The torque reduction is a model ingredient, not a claim that all fully convective stars lose their magnetic fields. Zorotovic and collaborators' detached-binary study tests the predicted detached population in the gap.
The standard formation channel starts with an initially wider binary star containing two main sequence stars. The initially more massive component evolves first and becomes a giant. Unstable Roche-lobe overflow can engulf the companion in a common envelope. Drag causes inward orbital motion, releasing two-body orbital energy and transferring angular momentum to the envelope. If the envelope is expelled before merger, a close detached binary survives, containing the exposed core, which becomes a white dwarf, and the lower-mass companion. This explains how a binary becomes much tighter than the giant progenitor's radius would have allowed. The outcome depends on envelope binding and the efficiency of energy deposition, summarized approximately by the common-envelope energy formalism; ejection is not guaranteed. Ivanova and collaborators' common-envelope analysis discusses the relevant physics and uncertainties.
Subsequent magnetic braking of a binary star and gravitational-wave emission from a binary system shrink the detached binary until its companion fills its Roche lobe. Stable Roche-lobe overflow, for a suitable binary mass ratio and stellar radius response exponent, then creates a cataclysmic variable. Its later nuclear eruptions, disk cycles and secular orbital period evolution occur on different timescales. The formation sequence is a wide binary, envelope ejection, a close detached white-dwarf binary, and angular-momentum-driven contact.
Use the Minkowski metric with signature . The source point is the intersection of the worldline with the observation event's past light cone: and . In components this is with . For a timelike worldline the retarded solution is unique in the usual complete-source setting.
Writing , and gives . Thus the Lienard-Wiechert potentials are and . At leading order in and first order in , set and expand the denominator:
Now use and let . Derivatives of and give near-field terms of order ; the radiation terms come from and . The electric field and magnetic field are therefore
They are transverse to the observation direction. The Poynting vector is
Integrating over a sphere uses , yielding the nonrelativistic Larmor formula . For simple harmonic motion, , whose period average is . Multiplication by the period gives
This is the leading radiation energy in the nonrelativistic and far-field approximation.
For the magnetic part of the Lorentz force, . Since velocity has physical dimension , the magnetic field has dimension
For the field generated by the infinitely long, axially symmetric shell currents, axisymmetry and the right-hand rule give . A circular loop has . Its enclosed electric current is zero inside the inner shell, between the shells and outside both. Ampère's law therefore gives
The positive direction is counterclockwise viewed from positive . The ideal current sheets make the tangential magnetic field discontinuous at each shell; these formulas specify the three open regions, with no additional imposed background field.
For a moving closed circuit bounding an oriented surface , Faraday's law is
Here is the electric field, the magnetic field, the local velocity of a circuit element, and the electromagnetic electromotive force per unit charge around the circuit. The line orientation and normal obey the right-hand rule. The flux derivative is a total derivative and includes both changing field and changing surface. The term is the motional electromotive force. A battery contributes an additional non-electromagnetic source electromotive force.