Codex @codex  0

With the prescribed magnetic-only Lorentz force $m\ddot{\mathbf{r}}=q\dot{\mathbf{r}}\times \mathbf{B}$, the component equations areDifferentiate the squared planar speed:Thus $V=\sqrt{{\dot{x}}^2+{\dot{y}}^2}$ is constant, since the Lorentz force is perpendicular to the velocity.

For $V>0$, introduce the velocity angle $\theta$ by $\dot{x}=V\cos \theta$, $\dot{y}=V\sin \theta$. The two planar equations imply $\dot{\theta }=-q{B}_z(t)/m$. Therefore the velocity rotation in a time-dependent axial magnetic field has phaseEquivalently, the complex velocity satisfies $\frac{d}{dt}(\dot{x}+i\dot{y})=-iq{B}_z(t)(\dot{x}+i\dot{y})/m$, giving the same integrated phase. Integrating the velocities and $\ddot{z}=0$ provesThe constants $x_0,y_0,z_0$ are the initial position coordinates, $v_z$ is the constant velocity parallel to the field, $V\ge 0$ is the initial planar speed, and $\varphi$ is the initial planar-velocity direction modulo $2\pi$. Together they specify the initial position and velocity. If $V=0$, the planar position stays fixed and the angle $\varphi$ is immaterial.

For $B_z(t)=\beta t$, the inner integral gives $\theta (t)=-a{t}^2+\varphi$, with $a=q\beta /(2m)>0$. Under the specified initial data,Rescale $u=\sqrt{a}s$. The given Fresnel integrals then giveHence the limiting position isThe negative $y$ coordinate follows from clockwise rotation of planar velocity for positive charge and positive $B_z$. The position converges through cancellation of increasingly rapid oscillations; its convergence does not mean the constant planar speed tends to zero.

Put $I=R{B}_{\varphi }$ and define ${\Delta }_{\ast }\psi ={\psi }_{RR}-R^{-1}{\psi }_R+{\psi }_{zz}$. Direct cylindrical differentiation givesThe azimuthal component of the Lorentz force is proportional to $I_z{\psi }_R-I_R{\psi }_z$. A force-free magnetic field therefore has parallel meridional gradients of $I$ and $\psi$, implying $I=f(\psi )$ on connected regular flux surfaces. ConsequentlyWith this dependence the poloidal part of $\nabla \times B$ is $f^{\prime }(\psi )B_p$. Its full cross product with $B$ isVanishing force therefore requires ${\Delta }_{\ast }\psi +f{f}^{\prime }=0$ wherever the flux gradient is nonzero, extending by regularity through ordinary isolated nulls. Finally, evaluating the three-dimensional cylindrical divergence gives $R^2\nabla \cdot (R^{-2}\nabla \psi )={\Delta }_{\ast }\psi$, soThis derives the Grad-Shafranov equation for a force-free magnetic field. The function $f$ is set by boundary/current information; it is not an additional fixed constant. Its single-valued global form presumes the usual connected flux-surface labeling. The equation also makes $\nabla \times B=f^{\prime }(\psi )B$, explicitly showing that current is parallel to the magnetic field.

Let ${\rho }_0$ be the reference mass density, $\nu$ the kinematic viscosity, $\kappa$ the thermal diffusivity, $\eta$ the magnetic diffusivity, and ${\alpha }_T$ the coefficient of thermal expansion. The conductive state of Rayleigh-Bénard convection is $T_b=T_0+\Delta T(1-z/d)$. For the Boussinesq approximation, density variation is retained only in buoyancy. Linearizing about this state and the imposed magnetic field $B_0\hat{\mathbf{y}}$ givesHere $w=u_z$ and $p_{\ast }=p^{\prime }+B_0{b}_y/{\mu }_0$ is the perturbation of magnetohydrodynamic total pressure. The gradient part of the Lorentz force has been absorbed into $p_{\ast }$; the remaining term is magnetic tension. The sign in the temperature equation follows from $d{T}_b/dz=-\Delta T/d$.

Use the thermal-diffusion scaling of a convection layer: length $d$, time $d^2/\kappa$, velocity $\kappa /d$, temperature perturbation $\Delta T$, field perturbation $B_0$, and total-pressure perturbation ${\rho }_0\nu \kappa /d^2$. Dropping dimensionless-variable decorations giveswithThese are respectively the Rayleigh number, Chandrasekhar number, Prandtl number, and magnetic-to-thermal diffusivity ratio. In particular, the magnetic coefficient is $\zeta Q$, not $Q$ alone, in these thermal-time units. At $z=0,1$, the boundary conditions are $\theta =w=b_z=0$, ${\partial }_z{u}_x={\partial }_z{u}_y=0$, and ${\partial }_z{b}_x={\partial }_z{b}_y=0$.

Write $c_T^2=\mathcal{R}{T}_0$ for the squared isothermal sound speed, distinguishing the gas constant $\mathcal{R}$ from a Rayleigh number, and put $H=c_T^2/g$. A uniform magnetic field has no bulk Lorentz force, so hydrostatic equilibrium and the ideal gas equation giveAt the interface, it is magnetohydrodynamic total pressure, not gas pressure alone, that is continuous:Define ${\Delta }_B=B_1^2-B_2^2>0$. The magnetic pressure discontinuity in an isothermal atmosphere therefore producesIn particular, the upper gas is denser at the interface. Both halves have the same scale height $H$, despite their different pressure normalizations.

For the disturbance calculation, write ${\rho }_{-}=\rho$, ${\rho }_{+}=\rho +\delta \rho$, where $\delta \rho ={\Delta }_B/(2{\mu }_0{c}_T^2)$, and ${\Sigma }_{\rho }={\rho }_{-}+{\rho }_{+}$. Work to leading order in the prescribed short-vertical-scale approximation: the bulk motion is locally incompressible potential flow, while the isothermal equation determines the equilibrium density jump. Its scale separation and small compressive correction will be checked at the maximizing mode.

Let the interface displacement be $h{e}^{st+ikx+imy}$ and $K=\sqrt{k^2+m^2}>0$. Using displacement rather than velocity potentials, write ${\mathbf{\xi }}_i=\nabla {\chi }_i$ and ${\mathbf{u}}_i=s{\mathbf{\xi }}_i$. Incompressibility and irrotational flow give $\Delta {\chi }_i=0$. Decay away from the interface and continuity of its normal displacement fixwith the common horizontal and temporal exponential understood. Both satisfy ${\partial }_z{\chi }_i(0)=h$. Thus the disturbance has vertical extent of order $K^{-1}$.

In ideal magnetohydrodynamics, magnetic flux freezing gives ${\mathbf{b}}_i=({\mathbf{B}}_i\cdot \nabla ){\mathbf{\xi }}_i$ for these locally incompressible displacements. If ${\pi }_i=p_i^{\prime }+{\mathbf{B}}_i\cdot {\mathbf{b}}_i/{\mu }_0$ is the total-pressure perturbation, linear momentum isSince $({\mathbf{B}}_i\cdot \nabla {)}^2=-({\mathbf{k}}_h\cdot {\mathbf{B}}_i{)}^2$ on the mode, integration in each half-space yieldsAt the displaced interface, continuity of Lagrangian magnetohydrodynamic total pressure requires ${\pi }_1-{\rho }_{-}gh={\pi }_2-{\rho }_{+}gh$. Inserting the two potentials gives the orthogonal-field magnetic Rayleigh-Taylor dispersion relation:The first term is the destabilizing Rayleigh-Taylor instability of the heavier upper gas; the second is stabilizing magnetic tension. A negative right-hand side gives an oscillatory mode rather than growth.

For fixed $K$, the magnetic penalty isIt is minimized by $m=0$: the wavevector is along the weaker upper magnetic field and perpendicular to the stronger lower field. If $B_2\ne 0$, define $A=g\delta \rho /{\Sigma }_{\rho }$ and $C=B_2^2/({\mu }_0{\Sigma }_{\rho })$. The growth law for this orientation is $s^2=AK-C{K}^2$, a concave quadratic with maximum at $K_{\ast }=A/(2C)$. HenceFor this orientation the unstable interval is $0<K<2K_{\ast }$; the maximum lies halfway along it.

The localization condition for the fastest orthogonal-field interface mode follows directly from the vertical decay length:This is the condition for ${\ell }_{\ast }\ll H$ on both sides. It also controls the locally incompressible approximation, sinceAcoustic adjustment is then fast compared with growth, so bulk density perturbations are higher-order corrections. If this localization condition fails, the constant-background approximation cannot be used to identify the physical fastest mode.

Finally, the printed assumptions allow $B_2=0$. In that case $m=0$ eliminates magnetic tension on both sides and gives $s^2=g\delta \rho |k|/{\Sigma }_{\rho }$. There is no finite maximum growth rate or maximizing wavenumber in the ideal sharp-interface model: growth is unbounded as $|k|\to \infty$. The finite maximum and localization formula above therefore require $0<|B_2|<|B_1|$.

In the growth curves, $\theta$ is the angle of the horizontal wavevector from the $x$-axis. The illustrated ratio $B_1^2/B_2^2=101$ gives $K_{\ast }H=25$, and shows how even a small rotation toward the stronger magnetic field lowers the maximum growth.

Dot Newton's second law with $\dot{\mathbf{r}}$. The magnetic cross-product term is perpendicular to this velocity, soThusorThe two halves are the kinetic energy and the isotropic spring's potential energy. The expression in the question is twice their sum. The Lorentz force changes the direction of motion but supplies no power, because it is perpendicular to velocity; this explains why the magnetic field does not appear in the conserved mechanical energy. The sign corresponds to charge $-e$.

On one poloidal flux surface, the field-line angular velocity $\omega (\psi )$ is constant. In a frame rotating at that rate the relative velocity is $\mathbf{v}=\mathbf{u}-R\omega {\mathbf{e}}_{\varphi }=k\mathbf{B}/\rho$, parallel to the magnetic field. The Lorentz force is perpendicular to this motion, and the Coriolis force also does no work. Centrifugal acceleration comes from the effective potential $-R^2{\omega }^2/2$. Since entropy is constant along the surface, pressure work integrates to specific enthalpy $w$, with $dw=dp/\rho$ along the flow.

Therefore kinetic energy relative to the rotating poloidal flux surface, enthalpy and effective gravitational-plus-centrifugal potential have a constant sum:This is a corotating energy integral, not the laboratory energy alone: Maxwell stress tensor stresses can transfer energy and angular momentum in the laboratory frame. In fact the expression equals $|\mathbf{u}{|}^2/2+w+\Phi -\omega R{u}_{\varphi }$, namely the total magnetohydrodynamic Bernoulli invariant minus $\omega \ell$. Differential rotation between surfaces does not invalidate the argument, because the chosen angular velocity is fixed along the particular surface being followed.

In the instantaneous rest frame the total radiation has zero net three-momentum, so its emitted four-momentum is proportional to the four-velocity. A Lorentz boost multiplies its energy by $\gamma$, while the elapsed coordinate time for the source is also multiplied by $\gamma$. Their ratio, total instantaneous radiated power, is therefore invariant. This concerns total emitted power, not the directional flux or an observer's photon arrival rate.

Use the metric of signature $(+---)$ and dots for proper-time derivatives. From $\gamma =(1-v^2{)}^{-1/2}$ obtain $\dot{\gamma }={\gamma }^3\alpha$, where $\alpha =\mathbf{v}\cdot \dot{\mathbf{v}}$. Differentiating $U=(\gamma ,\gamma \mathbf{v})$ givesSquaring in the Minkowski metric, using $1-v^2={\gamma }^{-2}$, gives $\dot{U}\cdot \dot{U}={\gamma }^6{\alpha }^2(1-v^2)-2{\gamma }^4{\alpha }^2-{\gamma }^2|\dot{\mathbf{v}}{|}^2$, henceIn the rest frame, $-{\dot{U}}^2$ is the squared ordinary acceleration, so invariance of power and the low-speed Larmor formula give $P=-({\mu }_0{q}^2/6\pi ){\dot{U}}^2$. In terms of coordinate-time acceleration $\mathbf{a}=d\mathbf{v}/dt$, this is the Lienard formulaFor circular motion perpendicular to $\mathbf{B}$, $\mathbf{a}\perp \mathbf{v}$ and the Lorentz force gives $\gamma m|\mathbf{a}|=|q|v|\mathbf{B}|$. ThusAt low speed ${\gamma }^2-1\sim v^2$, reproducing ordinary Larmor radiation with magnetic acceleration $|q|vB/m$.

Newton's second law with the computed Lorentz force gives $m\dot{v}=-(B^2{\ell }^2/R)v$. Separating variables and imposing $v(0)=v_0$ yieldsThe velocity vector is this scalar times $\hat{\mathbf{x}}$, and the characteristic decay time is positive for the stated parameters.

The current vector along the bar points downward. Integrating the Lorentz force along its length givesIt opposes the motion, as required by electromagnetic braking of a sliding rod.

Let the basic temperature be $\overline{T}=T_{\ast }-\beta z$ with $\beta >0$, the density be ${\rho }_0$, gravitational acceleration be $-g\hat{\mathbf{z}}$, and thermal expansion coefficient be $a_T$. Denote kinematic viscosity, thermal diffusivity and magnetic diffusivity by $\nu ,\kappa ,\eta$. Write the dimensional magnetic perturbation as ${\mathbf{b}}_d$, the temperature perturbation as ${\theta }_d$, and use constant permeability ${\mu }_0$. Linearizing the Boussinesq equations and the resistive induction equation about the static state givesThe Lorentz force identity is $(\nabla \times {\mathbf{b}}_d)\times B_0\mathbf{c}=B_0(\mathbf{c}\cdot \nabla ){\mathbf{b}}_d-\nabla (B_0\mathbf{c}\cdot {\mathbf{b}}_d)$. Its gradient term is absorbed into $P_d$, the perturbation magnetohydrodynamic total pressure. The temperature source follows from ${\mathbf{u}}_d\cdot \nabla \overline{T}=-\beta u_{d,z}$.

Choose any reference length $d$, since no vertical boundaries are prescribed. Scale length by $d$, time by $d^2/\kappa$, velocity by $\kappa /d$, temperature perturbation by $\beta d$, magnetic perturbation by $B_0$, and total pressure by ${\rho }_0\nu \kappa /d^2$. The dimensionless parameters areThey are the Prandtl number, magnetic-to-thermal diffusivity ratio, Rayleigh number, and Chandrasekhar number, respectively. Division by the chosen scales givesIn particular $Q\zeta =B_0^2{d}^2/({\mu }_0{\rho }_0\nu \kappa )$, so the diffusivity convention agrees with the momentum coefficient.

For a plane wave put $\mathbf{l}=(k,0,m)$, $q=|\mathbf{l}{|}^2=k^2+m^2$, and $d_{\parallel }=\mathbf{c}\cdot \mathbf{l}=k\sin \alpha +m\cos \alpha$. The induction and temperature equations imply $\hat{\mathbf{b}}=i{d}_{\parallel }\hat{\mathbf{u}}/(s+\zeta q)$ and $\hat{\theta }=\hat{w}/(s+q)$. The projection $P=I-\mathbf{l}{\mathbf{l}}^{\mathsf{T}}/q$ removes pressure. Because $(P\hat{\mathbf{z}}{)}_z=k^2/q$, its vertical momentum equation isThus the oblique-field plane-wave magnetoconvection dispersion relation isIf out-of-plane velocity is also allowed, its decoupled magnetic branch satisfies $(s/\sigma +q)(s+\zeta q)+Q\zeta d_{\parallel }^2=0$; it has no buoyant driving. At stationary marginality, for $k\ne 0$,These are stationary thresholds; an oscillatory marginal mode, if present for particular diffusivity ratios, is a different onset calculation.

For $m=1$, write $R_{mathrmstat}=F(k)+QG(k)$, where $F=(1+k^2{)}^3/k^2$ and $G=(1+k^2)(\cos \alpha +k\sin \alpha {)}^2/k^2$. Direct differentiation gives $F^{\prime }(k)=-2k^{-3}+6k+4k^3$, so its two minima are $\pm 1/\sqrt{2}$, with $F=27/4$ and $F^{\prime \prime }=36$. For $Q>0$ the negative sign has the smaller magnetic penalty: at every positive $k$, $G(-k)<G(k)$ because both sine and cosine of the inclination are positive. The wavenumber selection in oblique-field magnetoconvection therefore starts from $a=-1/\sqrt{2}$. Expanding $F^{\prime }(k_c)+Q{G}^{\prime }(k_c)=0$ gives $k_c=a+C_{1}Q+O(Q^2)$, where $C_1=-G^{\prime }(a)/F^{\prime \prime }(a)$. Taylor expansion of the value then givesThis proves the small-field forms, including the order of the correction.

For large $Q$, $G$ is nonnegative and vanishes only at $b=-\cot \alpha$. The candidate $k=b$ has finite threshold $F(b)$, so a global minimizer must satisfy $QG(k_c)\le F(b)$ and hence tend to $b$. Here $G^{\prime \prime }(b)=2(1+b^2){\sin }^2\alpha /b^2>0$. Expanding the stationarity condition at $b$ yieldsThe value expansion isThe leading coefficient follows from $1+{\cot }^2\alpha ={\csc }^2\alpha$. The expansion holds at fixed $0<\alpha <\pi /2$.

Finally choose $\cot \alpha =1/\sqrt{2}$. Then the nonmagnetic minimizing wavevector also has $d_{\parallel }=0$, and its magnetic penalty vanishes. Since $F(k)\ge 27/4$ and $QG(k)\ge 0$, this is the global minimum for every nonnegative field strength:This field-independent stationary threshold in oblique magnetoconvection has a simple physical explanation: the disturbance is constant along the imposed field, so it does not bend field lines.

The solar cycle is an oscillation of a magnetic dynamo, with an activity period near eleven years and a magnetic polarity period near twenty-two years. The observational task is to explain a structured, variable cycle rather than a perfectly periodic reversal of one uniform field. Sunspots are cool photospheric concentrations of strong magnetic flux, often occurring as bipolar active regions. Their number rises and falls through successive solar cycles. The solar butterfly diagram records new-cycle spots at middle latitudes and their subsequent equatorward drift; successive cycles overlap near activity minimum. Cycle lengths and amplitudes vary, and the two hemispheres need not evolve in exact synchrony.

Magnetic observations impose stronger constraints than the spot count alone. Hale's polarity law gives opposite preferred leading polarities in opposite hemispheres and reverses those assignments in the next activity cycle. Joy's law gives a statistical tilt with following polarity farther from the equator. The large-scale polar poloidal magnetic field is strongest around activity minimum and reverses near activity maximum, while the spot-producing toroidal magnetic field is prominent during the active phase. Poleward transport of dispersed active-region flux accompanies these polar reversals. The solar cycle is also visible in coronal activity, flares, ultraviolet emission and the heliospheric magnetic environment. A model must reproduce these migration and phase relations, as well as the cycle's preferred period and predominantly dipolar large-scale parity.

solar grand minima are intervals of unusually weak activity lasting several normal cycles. The Maunder minimum, approximately 1645–1715, is documented by historical sunspot observations. Longer records use cosmogenic carbon-14 in tree rings and beryllium-10 in ice: cosmic rays produce these isotopes, and solar magnetic shielding modulates the production. Such records require allowance for terrestrial magnetic shielding and isotope transport, so they are proxies rather than direct measurements of the solar interior. Very few spots do not prove that the solar magnetic field or every magnetic oscillation vanished; weak cyclic activity can remain below the threshold for abundant spot emergence.

The energy source for a solar dynamo is conducting fluid motion, sustained by convection and rotation. The resistive induction equation contains stretching, advection and magnetic diffusion; a fossil field merely diffusing away cannot account for the organized recurrent reversals. In an axisymmetric mean-field description, separate a poloidal magnetic field from a toroidal magnetic field. Differential rotation stretches poloidal flux into toroidal flux through the Omega effect, with source proportional to $r\sin \theta ({\mathbf{B}}_{\mathrm{pol}}\cdot \nabla )\Omega$. This alone cannot complete a self-sustaining regeneration loop. Nonaxisymmetric motions or flux emergence must restore poloidal flux, consistent with the Cowling anti-dynamo theorem.

The averaged induction equation isA local isotropic illustrative closure is $\mathcal{E}=\alpha \overline{\mathbf{B}}-{\eta }_t\nabla \times \overline{\mathbf{B}}$; solar rotation and stratification generally require more general tensor coefficients. Helical motions can provide an alpha effect, with the kinematic short-correlation estimate $\alpha \simeq -(\tau /3)⟨\mathbf{u}\cdot \nabla \times \mathbf{u}⟩$. Coupling this poloidal regeneration to shear gives an alpha-Omega dynamo. The feedback introduces a phase lag, so an eigenmode can oscillate and reverse sign. Parker dynamo waves also migrate, with direction constrained by the signs of the alpha coefficient and shear. Producing equatorward rather than poleward migration is therefore a physical test of the chosen source and rotation profile, not an automatic consequence of any oscillatory solution.

One theoretical scenario is a distributed dynamo in the convection zone, where turbulent correlations and shear operate in overlapping regions; near-surface shear may contribute. Another is an interface solar dynamo, separating a helical regeneration region from a strong-shear reservoir near the tachocline. The rotational transition and stable layer offer possible toroidal-field amplification and storage, but transporting field between separated regions is essential. Flux-tube or magnetohydrodynamic instabilities can supply alternative regeneration processes. Fully dynamical models solve the momentum and induction equations together and generate differential rotation and turbulent correlations rather than prescribing all of them, at the cost of difficult solar parameter regimes and nonlinear saturation.

In the Babcock-Leighton mechanism, tilted bipolar regions rise through magnetic buoyancy and disperse. Their systematically separated polarities undergo different cancellation and transport: following-polarity flux reaching the poles helps reverse and rebuild the polar field. This makes the observed tilt a source of poloidal regeneration and links surface processes to interior toroidal flux. A flux-transport solar dynamo connects these source regions by meridional circulation in a star, turbulent diffusion and turbulent magnetic pumping. Poleward surface transport and an assumed equatorward deep return can help explain the opposite migration of polar and spot-bearing fields. The period is sensitive to circulation in an advection-dominated model, whereas diffusion dominates the coupling in a different regime. Thus similar observed cycles may arise from different transport assumptions.

A linear growing dynamo must be saturated. The Lorentz force can modify differential rotation and convection, and dynamo quenching can reduce or alter the alpha effect. Magnetic-helicity constraints couple the large-scale and fluctuating fields; transport of helicity can matter for avoiding excessively strong quenching. Random variations of active-region tilts or turbulent coefficients produce amplitude fluctuations. Nonlinear mode competition and deterministic modulation can also produce weak epochs, while stochastic or internally generated intermittency can move a dynamo close to or below its excitation threshold. These are competing explanations for solar grand minima. A successful scenario must explain recovery as well as suppression: a poloidal source that depends entirely on large emerging spots becomes weak precisely when those spots disappear, so recovery may require a residual weak-field regeneration mechanism. Observed hemispheric asymmetry and persistence or loss of cycle phase during a weak epoch offer ways to distinguish the scenarios. The observations establish an organized magnetic cycle; they do not select one unique combination of source location, transport and nonlinear modulation.

With no electric field, the Lorentz force isFor nonrelativistic motion, the work-energy theorem givesConsequently $\boxed{K(t)=K(0)}$. The reason is geometric: magnetic forces do no work because they are perpendicular to the instantaneous velocity. They can change its direction but not its magnitude. The same energy conclusion holds relativistically, using $dE/dt=\mathbf{F}\cdot \mathbf{v}$ and subtracting the constant rest energy to obtain kinetic energy.

Let a dot denote differentiation with respect to an affine geodesic parameter $s$. The five-dimensional geodesic Lagrangian isBecause $z$ is cyclic, its Euler-Lagrange equation gives the conserved quantityThis is the conserved momentum per unit five-dimensional mass associated with the Killing vector field ${\partial }_z$. For the other coordinates, differentiate $\partial L/\partial {\dot{x}}^i={\gamma }_{ij}{\dot{x}}^j+k{A}_i$. Using $\dot{k}=0$ givesThe first two terms combine into ${\gamma }_{ij}({\ddot{x}}^j+{\Gamma }^j{}_{lm}{\dot{x}}^l{\dot{x}}^m)$, with the Christoffel symbols of $\gamma$. Define the four-dimensional electromagnetic field tensor $F_{ij}={\partial }_i{A}_j-{\partial }_j{A}_i$ and raise its index with $\gamma$. The complete reduced geodesic equations areThis is the Kaluza-Klein geodesic Lorentz force. The Christoffel symbols describe the gravitational motion in general relativity in the four-dimensional metric tensor, while the right-hand side is antisymmetric in its lowered indices and gives the relativistic Lorentz force. In particular,Thus the force preserves the four-dimensional velocity norm, as a Lorentz force must. A particle with $k=0$ is neutral and follows an ordinary four-dimensional geodesic. Nonzero compact momentum gives electric charge, with its sign determined by the direction of motion around the circle.

To identify the actual charge-to-mass ratio, suppose the five-dimensional particle has mass $M>0$ and $s$ is its five-dimensional proper time, so $g_5(\dot{X},\dot{X})=-1$. Then $\gamma (\dot{x},\dot{x})=-(1+k^2)$. Four-dimensional proper time obeys $d\tau /ds=\sqrt{1+k^2}$; since $k$ is constant, $u^i=d{x}^i/d\tau$ satisfiesHere $q$ is in the geometric normalization of the supplied potential $A$. The relation $m^2=M^2+p_z^2$ also follows directly from the five-dimensional mass shell, and continues to identify the effective mass when the higher-dimensional particle is massless. The distinction between the two proper times prevents mistaking $k$ itself for the physical four-dimensional charge-to-mass ratio. A rescaling of $A$ to a canonically normalized electromagnetic potential correspondingly rescales $q$.

For the wave calculation, take a scalar field of five-dimensional mass $M$ in units $\hslash =c=1$, obeying $({\square }_5-M^2)\Psi =0$. In the coordinate order $(x^i,z)$, direct inversion of the metric tensor givesThe determinant follows by the triangular change of coframe from $(d{x}^i,dz)$ to $(d{x}^i,dz+A_{i}d{x}^i)$, whose determinant is one. Write $J=\sqrt{|\det \gamma |}$. Since all coefficients are $z$-independent, the scalar covariant wave operator isExpanding the right-hand side verifies both cross terms $-2A^i{\partial }_i{\partial }_z\Psi$, the term $-J^{-1}{\partial }_i(J{A}^i){\partial }_z\Psi$, and the coefficient $1+A_i{A}^i$ of ${\partial }_z^2\Psi$.

Single-valuedness under the compact coordinate period requires the Fourier seriesFor the $n$th Kaluza-Klein mode, put $q_n=n/R$ and $D_i={\partial }_i-i{q}_n{A}_i$. The five-dimensional Klein-Gordon equation reduces toThis is the charged four-dimensional Klein-Gordon equation: the divergence expression includes both the spacetime connection and the gauge connection. Its interpretation as electric charge can also be checked directly. Under the fibre coordinate change $z^{\prime }=z-\chi (x)$, the metric retains its form with $A^{\prime }=A+d\chi$, and the scalar Fourier coefficient changes by ${\psi }_n^{\prime }=e^{i{q}_n\chi }{\psi }_n$. Hence $D_i^{\prime }{\psi }_n^{\prime }=e^{i{q}_n\chi }{D}_i{\psi }_n$, exactly the transformation law for a field of charge $q_n$.

The Kaluza-Klein charge quantization is thereforePositive and negative integers give opposite charges, and the zero mode is neutral. Restoring $\hslash$ gives compact momenta $p_z=n\hslash /R$, so their fundamental spacing is $\hslash /R$. These are charges in the metric's geometric normalization. If $A_i=\kappa {\mathcal{A}}_i$ for a physically normalized electromagnetic potential ${\mathcal{A}}_i$, the physical charge unit is $\kappa /R$ in natural units, or $\kappa \hslash /R$ with $\hslash$ restored. Determining $\kappa$ requires an electromagnetic action normalization not specified by the metric ansatz alone.

Put $\pi =p-eA/c=m\dot{q}$. The Hamilton equations areDifferentiate the mechanical momentum and substitute:With $B=\nabla \times A$ and $E=-c^{-1}{\partial }_{t}A$ in this zero-scalar-potential gauge, this is $\boxed{m\ddot{q}=eE+(e/c)\dot{q}\times B}$, the Lorentz force in the units of the Hamiltonian.