Away from the origin, the electric field is , with here denoting the position vector in the numerator and its Euclidean norm in the denominator. The spherical-coordinate formula for divergence gives
At the origin the electric field is singular, so the ordinary divergence is undefined there. The equation above must not be extended to the whole space as an ordinary smooth-field identity. Its correct distributional divergence, derived in solution, includes the point source.
Use the real magnetic field in the coupled equations; its leading periodic approximation will be . In the conducting incompressible flow, the magnetohydrodynamic momentum equation, resistive induction equation and solenoidal magnetic-field constraint are
Here is pressure, is electric current density, and and are magnetic diffusivity and kinematic viscosity. The moving-conductor Ohm law gives .
In the insulating exterior, neglecting displacement current, the magnetic field is curl-free away from the wire and solenoidal. The imposed electric current density supplies
The exterior electric field also satisfies Faraday's law. At the rigid wall, impose the no-slip boundary condition . At finite conductivity, the interface conditions for electromagnetic fields give continuous normal magnetic field, continuous tangential magnetic field for equal permeability and no prescribed singular surface current, and continuous tangential electric field. The normal electric current density is zero at an insulating wall; it is automatic for the present two-dimensional fields with current only along . The induced fields decay far from the wire and wall, and initial conditions, or selection of the periodic state after transients, complete the specification. In the limiting perfect conductor approximation a surface current can emerge, permitting a tangential-field jump.
The fully coupled solution need not remain monochromatic: a harmonic magnetic field produces both a mean and a twice-frequency Lorentz force density, and the resulting velocity can generate further harmonics. The following parts consistently use the leading magnetic response with fluid motion neglected.
The wind integrals require a steady axisymmetric magnetohydrodynamic wind in ideal magnetohydrodynamics, with no viscous azimuthal torque or imposed toroidal loop voltage. Axisymmetry by itself is insufficient. Assume smooth connected poloidal flux surfaces with nonzero poloidal magnetic field; functions of the flux label below are understood on each such branch. In cylindrical coordinates, the given poloidal magnetic flux function obeys
For a steady ideal wind, the electric field is and . Its toroidal component satisfies and , so . Regularity at the axis, or equivalently the absence of an imposed loop voltage in an annular domain, gives . Thus and . The continuity equation and now yield
Consequently is the magnetohydrodynamic mass loading. Including the toroidal magnetic field, rewrite the full velocity as
Since , steady ideal induction gives . Hence , establishing
The invariant is the field-line angular velocity of an axisymmetric wind, not generally the fluid's angular velocity . When an ideal magnetic surface is anchored to a rigidly rotating star, equals the stellar angular velocity. With no poloidal outflow this reduces to Ferraro's law of isorotation.
To calculate the magnetic torque, use the azimuthal momentum equation. Axisymmetric pressure and Newtonian gravity have no azimuthal components; the magnetic pressure has none either. The azimuthal acceleration and magnetic tension contain the curvature terms characteristic of cylindrical coordinates:
Multiplication by gives the Maxwell torque conservation in an axisymmetric wind equation
Using , and proves
This is the total angular momentum transport per unit poloidal magnetic flux. The first term is the material contribution and the second is the contribution of the Maxwell stress tensor. The specific magnetohydrodynamic angular-momentum invariant is when ; its normalization differs from the paper's .
Define the poloidal Alfvén Mach number by
Substituting into gives
At a smooth crossing of the Alfvén surface, and the toroidal magnetic field is finite. The Alfvén-surface regularity condition for an axisymmetric wind is therefore
One can also expose the apparent singularity explicitly:
Finite passage through the Alfvén surface requires the numerator to vanish with the denominator. Equality of the numerator zeros is necessary; a globally smooth wind must also satisfy its other dynamical critical conditions. Here is the cylindrical distance from the rotation axis, not automatically the spherical Alfvén radius.
For a tube carrying outward mass flux , the outward angular momentum flux is
Thus the stellar torque is , with the integral over the escaping wind. Magnetic stresses give each mass element an effective angular-momentum lever arm equal to the cylindrical Alfvén crossing radius. In a strongly sub-Alfvénic region the fluid approximately corotates with the magnetic surface; beyond that region it can carry the accumulated angular momentum away. For an outward field and positive , a super-Alfvénic expanding tube with has , so the magnetic stress transports positive angular momentum outwards. If greatly exceeds the stellar surface lever arm, a modest mass-loss rate can cause substantial wind-driven magnetic braking of a solar-type star. An unmagnetized outflow would instead carry approximately the material angular momentum fixed at its launch radius.
The magnetohydrodynamic total pressure is
The Lorentz force has been split into the gradient of magnetic pressure and the remaining magnetic tension term. An appropriate boundary model is impermeable, perfectly conducting stationary cylindrical walls with no initial normal magnetic flux: at . In an inviscid fluid this is a slip condition, not a no-slip condition on or . It also makes the tangential ideal electric field vanish at the walls. Take perturbations periodic in and periodic or Fourier-resolved in ; the wall reaction supplies the normal pressure force without an extra prescribed pressure value.
Write the steady velocity and magnetic field as and . They satisfy both divergence constraints. Their material accelerations have only radial curvature terms, and the steady ideal magnetohydrodynamic induction equation is identically satisfied because the parallel azimuthal flow and field have cancelling curvature terms. The radial momentum equation is
Hence the most general such equilibrium is
Its gas pressure is ; on a bounded annulus the additive constant can be chosen to keep it positive. Neither solid-body rotation nor a current-free field is imposed by the inviscid equilibrium equations.
For the linear stability calculation, put , , and define . Use axisymmetric normal modes proportional to . Denote the velocity amplitudes by and the magnetic amplitudes in velocity units by . Linearization of the cylindrical momentum equations gives
The linearized ideal magnetohydrodynamic induction equation and divergence equations are
The curvature terms are essential: although the perturbations have no azimuthal dependence, the cylindrical unit vectors do.
For , . Introduce the radial Lagrangian fluid displacement by . The azimuthal amplitudes reduce to
Consequently the radial and vertical equations, combined with the incompressible flow constraint, give
where the Michael criterion for axisymmetric toroidal-field interchange coefficient is
For , the required single displacement equation is therefore
The same differential equation holds for , which is a constant multiple of .
Set . Multiplication by and integration by parts with the wall conditions gives the global variational form of the Michael criterion:
The denominator is strictly positive for a nonzero displacement, so every mode's is real. If everywhere, no exponentially growing mode exists. Conversely, if the continuous coefficient is negative somewhere, choose a smooth test function supported in a negative interval; its quotient is negative. For sufficiency, let be the positive Sturm-Liouville operator with Dirichlet boundary conditions and let multiply by . The generalized mode equation is equivalent to the compact operator and self-adjoint operator problem . The negative test quotient guarantees a negative eigenvalue and therefore a growing solution . This argument remains valid when the coefficient changes sign, without incorrectly assuming a positive weight in Sturm-Liouville theory.
Thus, for any fixed nonzero vertical wavenumber,
Equality is marginal. When , the incompressible flow constraint and both rigid-wall conditions force the radial velocity to vanish, so that special case has no radial interchange mode. When the field vanishes, the coefficient reduces to , recovering the Rayleigh centrifugal stability criterion.
For a Keplerian disk, . If , the toroidal interchange field threshold in a thin Keplerian disk is
For ordinary radial gradients on the scale , the Alfvén speed must therefore be comparable to the orbital speed. A field on the gas pressure scale has in a thin disk, and cannot normally meet this condition. Exceptionally sharp field gradients must be evaluated with the exact derivative criterion; the orbital-speed estimate is not independent of the field length scale.
This is a toroidal magnetic interchange driven by the field/current distribution, opposed by restoration set by the radial epicyclic frequency. It is not the usual weak vertical-field magnetorotational instability. In the vertical-field problem, a nonzero bends magnetic field lines and couples displaced fluid elements, allowing decreasing to supply rotational kinetic energy even when . Here the axisymmetric perturbation has no scalar variation along the purely toroidal background field, and no such weak-field coupling: a sufficiently adverse toroidal-field gradient, rather than merely , is required. Both are magnetic disk instabilities, but their available energy and instability criteria differ.
Write horizontal components with a subscript . The steady ideal magnetohydrodynamic induction equation gives . Since , its horizontal components imply
Thus is a spatial constant. A vanishing axial electric field, , is needed to infer horizontal alignment. It is satisfied, for example, if the constant electric field vanishes at a boundary or at infinity. With this additional condition and , define the magnetohydrodynamic mass loading by , and put . Then
The steady continuity equation and the solenoidal magnetic-field constraint give
Also . Using the vector identity for the curl of a cross product, together with and , its curl is . Hence
These are steady planar ideal-MHD field-line invariants: both quantities are constant along a magnetic field line.
The zero-electric-field condition does not follow from the printed stationarity assumptions alone. Constant , , and , with , satisfy every displayed ideal magnetohydrodynamics equation but have perpendicular horizontal velocity and magnetic field. Their is nonzero, so no give the stated decomposition. The field-line decomposition is valid on the branch; as an unrestricted implication, the first claim has this counterexample. The subsequent calculations use that branch.
Let denote geodesic distance from the charge; is the fixed curvature radius, not the radial variable. In round three-dimensional hyperbolic space, a sphere has area
Using the convention that electric flux is , Gauss's law gives the radial orthonormal electric field
Its magnitude is , and a negative charge reverses the direction. In rationalized natural units set ; Gaussian conventions instead use flux .
For this approaches . For , . This hyperbolic-space Coulomb field falls rapidly because the sphere area grows exponentially; conserved electric flux is unchanged. A potential vanishing at infinity is .
The original PDF actually prints rather than in the angular metric tensor. The round angular factor used above is required by its description as hyperbolic space. If the printed angular metric tensor is taken literally, the sphere-area coefficient is instead
and a radial Gauss-law solution has . That angular geometry is not a smooth round sphere, so the literal metric tensor does not describe the asserted homogeneous hyperbolic space. Both readings retain the same radial area-growth factor; the geometric repair is explicit.
Fix the electric potential to vanish at infinity and assume a localized, nonsingular charge density of finite electrostatic energy. Assemble the charge by scaling from to . Linearity gives , so the work is
The factor prevents double counting of pair interactions. Using Gauss law, , and integration by parts,
The boundary integral at infinity vanishes for a localized distribution since and . Singular point-charge self-energies require a separate regularization and are not finite instances of this identity.
For the uniform solid sphere, spherical symmetry and Gauss law give
Integrating inward from infinity and matching the continuous electric potential at gives
The value of the electric field at the centre is the zero vector.
Writing , the electrostatic energy of a uniformly charged solid sphere follows either from the charge-potential integral or from the field integral including the exterior:
For the nuclear scaling, while constant volume per proton gives . Therefore the electric contribution scales as in this uniform-density model.
The magnetic vector potential gives the upper poloidal magnetic field , and is its toroidal magnetic field; play the same roles below. This representation automatically satisfies . Diffusion acts in both regions. Above the interface the alpha effect converts toroidal magnetic field into poloidal magnetic field, while below it a mean shear converts the poloidal component into toroidal magnetic field. Separating these source regions still allows an alpha-Omega dynamo through diffusion across the interface.
For electromagnetic matching, keep the closure actually represented by these equations. Its retained turbulent mean-field electromotive force is above and zero below: only toroidal-to-poloidal alpha generation is retained. Let the mean velocity be , continuous at the interface, and write the mean electric field as
This yields the specified scalar equations with above and below. In particular, the upper tangential electric component is
where is the scalar toroidal-field amplitude. The other tangential component is .
Assume equal permeability, no imposed singular surface current and ordinary electromagnetic matching. Continuity of the normal magnetic field gives , and tangential magnetic continuity gives and . For the nonzero horizontal Fourier wavenumber used below, normal-field continuity is equivalent to ; the spatially constant potential difference is a gauge freedom. Continuity of , with continuous and , gives the final condition. Thus electromagnetic matching for an alpha-Omega interface is
Continuity of follows from and their time evolution, rather than supplying a fifth independent amplitude condition. A full isotropic mean-field electromotive force would add toroidal alpha induction and an alpha contribution to ; that would be a different bulk and interface problem. Retaining it only in a boundary condition would be inconsistent with the equations being solved here.
Set , without loss of generality for a real wave and its conjugate, and write . Define the two complex vertical wavenumbers and so the upper fields decay as and the lower fields as . The original PDF has a lowercase in the lower phase. The homogeneous upper and lower equations give
Both and have positive real part. Since their squares agree, they must be equal, rather than opposite. Put
so and follow from the equations rather than being assumed in advance.
Suppress the common factor and write
The identity gives the upper forcing equation, while the same operator on gives . Hence
The first and third interface dynamo matching give , . The derivative conditions then give
Thus and . Eliminating the nonzero amplitudes produces the dispersion relation of an equal-diffusivity interface dynamo:
In terms of ,
Each branch has a square root with positive real part and is therefore compatible with vertical localization. The branch with the plus sign is the potentially growing one:
The minus branch has . At the first branch is neutrally oscillatory. The frequency sign follows from the stated convention: positive give propagation in the negative direction. Reversing the coupling sign reverses propagation, while the threshold for the conjugate sign depends on the magnitude of the dynamo number.
Figure 1.
Separated interface-dynamo source layers and the growth threshold at dynamo number 32
. Details.
The equal-diffusivity interface dynamo idealizes a solar configuration in which helical convection regenerates poloidal magnetic field above the tachocline and differential rotation generates strong toroidal magnetic field in or below it. Diffusion communicates the fields between the two regions, closing the feedback loop. This separation is the mechanism emphasized by Parker's interface-wave model. An oscillatory traveling mode gives migrating magnetic patterns and polarity reversals; buoyant toroidal flux can appear as sunspot belts. If is a local latitude coordinate, the direction of their migration depends on the alpha–shear sign and the coordinate orientation.
The magnetic period of this branch is ; polarity reverses after half a period, and activity that depends on field strength rather than sign can repeat on that half-period. This is the qualitative relation between the magnetic solar cycle and its activity cycle. The model does not itself predict the observed period or guarantee equatorward migration: those require actual geometry, transport coefficients and shear signs. Real stars also have diffusivity contrasts, meridional transport, spatially varying regeneration and nonlinear feedback. The kinematic exponential mode describes onset, while saturation and a maintained finite-amplitude cycle require further physics.
Decompose the electric field into its orthogonal parts , where . The uniaxial conductivity tensor acts on them as
If , its perpendicular and parallel components each vanish. Since and , this forces both field components to be zero. The contrapositive gives
Equivalently, the tensor has nonzero eigenvalues and determinant .
Writing the uniaxial conductivity tensor in vector form gives the electric current density
Choose the plane through the origin perpendicular to :
For an electric field in this plane, the anisotropic term vanishes and . Thus every nonzero such field is parallel or antiparallel to its current according to the sign of the nonzero . The zero field gives zero current. This establishes collinearity throughout the plane without assuming that the tensor is isotropic.
Take zero electric potential at infinity and put . By spherical symmetry and Gauss's law, the electric field is radial, with magnitude determined by the total enclosed electric charge. Thus, away from the idealized shell surfaces,
Integrating from infinity, or adding the three spherical-shell potentials, gives
The expressions agree at each interface, so the potential is continuous everywhere, including the origin. The electric field has different interior and exterior limits at each infinitely thin charged surface, with jump equal to surface charge density divided by . An ideal sheet has no single classical field value on the sheet itself; the displayed adjacent limits specify the field there. In conducting material of nonzero thickness the internal field is zero.
The electrostatic energy is the integral of . There is no contribution from , and spherical integration gives
The dependence on cancels because the magnitude of enclosed electric charge is on both sides of the middle shell, even though the field direction changes.
Take the divergence of the Ampere-Maxwell law and use Gauss's law . Since the divergence of a curl is zero, this gives the charge continuity equation . Uniform conductivity and Ohm's law imply . Therefore
The spatially uniform initial density stays uniform in the interior; is the charge relaxation time.
Denote the cylinder radius by . Cylindrical symmetry makes the interior electric field radial. Gauss's law on a coaxial cylinder of radius and length gives . Thus
The current does not disappear at the interface: the vacuum conducts no current, so charge accumulates on the cylindrical surface. Initially taking the charge to be the stated volume charge alone, its surface charge density satisfies , giving the charge relaxation in an isolated conducting cylinder law
The total enclosed charge per unit length is . Consequently the exterior electric field is
independent of time. The boundary jump obeys , as it must. Applying Gauss's law to the decaying volume charge alone would incorrectly predict a decaying exterior electric field and violate electric charge conservation. The vector formulas also specify the inward direction if .
In SI units, the source-free Maxwell equations in vacuum are
Use complex amplitudes with phase and take real parts at the end. The electric divergence equation gives . Faraday's equation gives , so
The magnetic divergence equation follows because . Ampere's equation becomes . Since , a nonzero wave satisfies
Thus both fields are transverse, mutually perpendicular for linear polarization, and obey the vacuum plane-wave dispersion relation.
For the interface let point from the minus side to the plus side and write , similarly for . Applying the integral Maxwell equations to a thin pillbox and a thin loop gives the electromagnetic boundary conditions
Here is the tangential surface current. The jump notation fixes the signs, including the direction of the normal.
For the perfect conductor, take the time-varying interior fields to vanish. Choose positive with . The incident electric field has complex form . To cancel its tangential value at the surface, use the reflected field , whose wavevector points upward. The total vacuum fields are therefore
The sign of each magnetic contribution follows from . At , the electric field is zero and the magnetic field is tangential. The jump conditions yield
For a real incident amplitude this is . It supplies the reflected wave while satisfying the normal-field and tangential-electric conditions. Any independently imposed static interior field is separate from this oscillatory response.
Use rationalized electromagnetic units with and electrostatic constant . The four-potential transforms by the same matrix as the coordinates:
Write and . Gauss's law for the line at rest gives , . A suitable potential is , . The inverse Lorentz transformation gives , . The fields from and are therefore
The sign agrees with a line charge moving in the positive direction: its density is and its current . In SI units insert in the electric expressions and replace by .
The two invariants check explicitly:
For , the electric field differs from the static field only at order , while the magnetic field is of order and is the field of the moving line current. These fields are evaluated off the charged axis.
Neglect end effects, so the electrostatic potential depends only on the cylindrical radius . The radial Laplace equation is , with general solution . Imposing the conductor boundary conditions gives
The radial electric field is . Just inside the middle cylinder it is , directed towards the inner cylinder; just outside it is , directed towards the outer cylinder. By Gauss's law, its charge per unit length is
Hence the capacitance per unit length is
The two coaxial cylindrical capacitors contribute in parallel because both connect the middle cylinder to the same grounded potential. This is the coaxial capacitor with grounded inner and outer cylinders construction.
For with , while . Therefore
The inner gap has width and area per unit length approximately . Its parallel-plate capacitor approximation gives exactly the same leading result; the outer gap contributes only a bounded correction.
For the extremum, set , , so . Then
Thus the unique stationary point is . Strict convexity and divergence at both endpoints show it is a global minimum:
There is no maximum on .
For a moving circuit, Faraday's law takes the form
Here is the instantaneous wire contour, its local material velocity, and any spanning surface. The contour orientation and surface normal vector obey the right-hand convention. The integral on the right is the magnetic flux, differentiated with both the field and circuit position allowed to vary. The electric field term is the electric force per unit electric charge around the wire, and supplies the motional electromotive force. Their sum is the electromotive force driving charge around the moving circuit. The condition makes the magnetic flux independent of the spanning surface.
Choose the initial loop normal vector to be , so rotation about gives . The spatially uniform magnetic field gives
Using Ohm's law for the resistive loop,
Positive electric current is measured in the contour orientation associated with . Reversing that convention reverses the signed answer. When , the magnetic flux is constant and the induced electric current vanishes. This uses the usual thin resistive-circuit approximation with self-inductance neglected.
Use a steady axisymmetric flow governed by ideal magnetohydrodynamics and take the Newtonian gravitational potential to be axisymmetric. Work locally on connected regular magnetic flux surfaces with nonzero poloidal flow. As usual for a wind without an imposed toroidal loop voltage, ; regularity on an included rotation axis guarantees this condition. It is needed to exclude steady cross-field flow in an annular domain.
Since , the poloidal magnetic flux function gives
The poloidal magnetic field is divergence-free. The zero toroidal electric field, with , gives . Define by . Steady mass conservation then gives
Thus is the magnetohydrodynamic mass loading, the mass flux per unit poloidal magnetic flux. Subtracting from the full velocity leaves only an azimuthal component. The cross product is
The steady ideal magnetohydrodynamic induction equation requires . Consequently the field-line angular velocity is and
Neither gas pressure nor axisymmetric gravity exerts an azimuthal torque. Multiplying the azimuthal momentum equation by combines the cylindrical curvature terms into angular-momentum transport:
Insert and use . It follows that
so the magnetohydrodynamic angular-momentum invariant is
The two terms account respectively for matter and Maxwell stress tensor transport of angular momentum.
Finally, subtracting times the logarithmic mass density equation from the logarithmic pressure equation gives . The specific entropy of a fixed-composition ideal gas is a function of . Steadiness and axisymmetry reduce its advection to . Hence
These flux-label functions are local on connected regular surfaces; disconnected components with the same numerical label need not share constants without an additional matching condition.
For the charge-generated electric field, use the electric potential
The vector potential is time independent, so its usual additional contribution is zero. The continuity equation now gives
Thus the second term in part (iii) is exactly the displacement current contribution , and
This is the Ampère-Maxwell equation. The constructed field also has zero curl, consistent with the time-independent magnetic field and Faraday's law. A static addition to the electric field does not change this conclusion; the integral construction selects the fields sourced by the given densities rather than an independently added radiation field.
Take the electric dipole moment along the polar axis, . The electric field is minus the gradient of the electric potential; in spherical polar coordinates,
Since and , the radial and angular components combine as
The origin is excluded because the point-dipole potential is singular there.
Spherical symmetry and the integral Gauss law give for . Together with the stipulated zero interior field,
The electric field jumps across the ideal charged surface; its value exactly on that surface is not needed for the volume integral. The electrostatic energy of a uniformly charged spherical shell is
As a check, the shell has constant electric potential , with zero reference at infinity, so gives the same result.