In charge-free, current-free vacuum, Maxwell's equations are
Taking the curl of the third equation and using the fourth, the curl of the curl identity, and the zero divergence of the electric field gives
The magnetic flux density obeys the same wave equation. For the plane wave factor , the spatial Laplacian contributes and the second time derivative contributes . Thus
for propagation in the positive direction. The zero-divergence conditions require , while the curl equations require . These conditions show that the proposed plane waves satisfy all four Maxwell's equations, rather than only their derived wave equations. Physical electric fields and magnetic fields are the real parts of the complex expressions.
A pure stationary state of a time-independent Hamiltonian operator is . Its probability density and all expectations of time-independent observables are constant. This does not require zero probability current.
For real scalar potential energy and no vector potential, substitute and its complex conjugate into . The potential terms cancel:
Thus the probability current is . In a stationary state : flow may circulate without changing the density.
Here only the azimuthal phase contributes an imaginary logarithmic derivative. Since ,
There is no radial or polar flow; the azimuthal current is independent of , so its divergence vanishes. The apparent coordinate singularity at cancels in the final expression.
The z-component of orbital angular momentum follows from the moment of the mass current:
Indeed . This agrees with acting on the wavefunction.
displacement current completes Ampère's law:
In vacuum . Taking the divergence and using gives , exactly the continuity equation. Without the added term a time-dependent charge distribution would contradict the curl equation.
A charging parallel-plate capacitor illustrates the physics. For plate area , , so the gap displacement current is , equal to the conduction current in its leads. An Ampère contour therefore gives the same magnetic circulation whether its spanning surface crosses the wire or the gap. No electrons have to cross an ideal vacuum gap. In matter, includes the polarization response as well as the vacuum term. Combined with Faraday induction, the displacement term permits source-free electromagnetic waves.
stress describes internal force per area. A cut with unit normal has traction ; its normal component is and its remainder is shear traction. force balance constrains the stress divergence, while ordinary local angular-momentum balance makes the Cauchy stress tensor symmetric in a medium without body couples.
strain describes deformation rather than force. For small displacement gradients,
The antisymmetric displacement-gradient part is a rigid rotation to first order and contributes no infinitesimal strain. Engineering shear strain is , so it must not be confused with the tensor component. stress has units of pressure; strain is dimensionless.
An isotropic linear elastic solid obeys
Two independent moduli are required. The shear modulus and bulk modulus must be positive for stability. Young's modulus and Poisson ratio satisfy
For a tensile specimen with free side faces, gives and . Pure shear gives ; hydrostatic loading gives . The incompressible limit has at finite shear stiffness, not vanishing shear resistance. Crystals can be anisotropic and need a fourth-rank elastic stiffness tensor instead of the isotropic two-constant law.
elastic deformation is recoverable on unloading; it need not remain linear at large strain. Beyond yield, irreversible rearrangements such as dislocation motion cause plastic strain, often followed by hardening. Ductile samples can neck before fracture; brittle failure has much less plastic deformation. Engineering stress uses initial area, while true stress uses the changing instantaneous area. A decreasing engineering stress during necking therefore need not mean that local true stress decreases. stress concentrations, flaws and loading geometry can determine failure before a uniform-material strength is reached. Viscoelastic materials also show rate dependence, creep and stress relaxation; their losses cannot be described by a static Hooke law alone.
A tensile test should measure force and gauge-length extension independently, correct machine compliance, align the grips, and control strain rate and temperature. Unloading distinguishes recoverable strain from a residual plastic extension. Separate shear or hydrostatic tests determine moduli that one axial measurement does not fix. Compression specimens need friction and buckling checks. Work per volume is ; the recoverable part is elastic energy, while a dissipative cycle produces a loss rather than creating energy. The original loading/unloading curve below is a chosen constitutive sketch, not measured material data.
Figure 1. . Small-strain geometry and illustrative unloading. Details.
At fixed positive pressure, the extrapolated term tends to as . That prediction is not physically realistic for an equilibrium sample. The classical ideal-gas description with constant heat capacity cannot extend unchanged to absolute zero. A real substance generally undergoes phase changes, and quantum statistics and thermal freeze-out change its heat capacity and entropy. For a nondegenerate equilibrium ground state, the third law gives entropy approaching zero; residual ground-state degeneracy can give a finite residual value rather than this negative divergence. This is a breakdown of the extrapolation, not a violation of the third law by a real gas.
At fixed positive temperature, tends to as . This is a physically meaningful unbounded-volume limit for a fixed mole number: , so the number of accessible translational configurations grows without bound. For an ideal gas the entropy difference is at fixed temperature. No finite vessel reaches an infinite entropy, and at finite temperature is not a finite-volume one-mole state. The divergence is approached by dilution; it does not conflict with the third law, which concerns a different temperature limit.
The microscopic magnetic response becomes a macroscopic field problem through the magnetization , the magnetic dipole moment per unit volume. Define
in magnetostatics. Bound currents are represented by and surface current . Normal is continuous across an interface, while . It is generally incorrect to impose continuity of every component of both fields. A linear isotropic material has and , but a ferromagnet usually requires a nonlinear, history-dependent constitutive relation.
Internal field differs from externally applied field. For uniform magnetization along a principal axis of an ellipsoid, the demagnetizing field is and . A linear specimen therefore has apparent susceptibility relative to the applied field. Shape is especially important for large susceptibility. A long closed magnetic circuit reduces stray fields; a bar magnet still has demagnetizing fields despite the absence of magnetic monopoles. Equivalent magnetic-pole descriptions are bookkeeping for magnetization divergence, not free magnetic charge.[ref]
In paramagnetism, permanent localized moments are preferentially aligned by the field, against thermal disorder. For independent spin-one-half moments with projections , the two Zeeman energies give
in the dilute weak-field limit with . This is a Curie law; the response saturates at when alignment dominates. More generally, spin gives . Interactions can produce a Curie–Weiss response over a suitable range. Conduction-electron Pauli paramagnetism need not show the localized-moment law.
In diamagnetism, the field induces a magnetic moment opposing it, so the linear susceptibility is negative. The microscopic cause is the change in orbital electronic motion, and a closed-shell insulating material has no permanent atomic moment to orient. This weak response is often only mildly temperature dependent away from other transitions. It also exists as a background in materials with a stronger positive response. Perfect flux expulsion in a superconductor is a distinct collective phenomenon and should not be identified with ordinary weak atomic diamagnetism.
In ferromagnetism, exchange interactions stabilize spontaneous magnetic order below the Curie temperature. A specimen can have zero net moment while its magnetic domains remain locally magnetized in different directions. Domain-wall motion and moment rotation produce a large response; anisotropy and pinning lead to remanence, coercivity and hysteresis. Increasing field eventually approaches saturation. Heating above the ordering transition removes spontaneous order, although an applied field still induces a response. Soft magnetic cores favor easy reversal and small cyclic loss; hard permanent-magnet materials favor retained magnetization and resistance to reversal. These properties concern collective order and barriers, not merely independent moments all following a weak field.[ref]
A practical distinction requires a controlled field sweep and a separate magnetic readout. A toroidal core with drive winding , path length and current has approximately . A pickup winding measures ; integrating after calibration gives , and . Reversing the sweep reveals a loop whose positive loss area per volume is . An air-core/empty-holder reference, known geometry and turn count, offset correction, temperature control and slower sweeps test background, heating and rate-dependent eddy-current losses. A single force or susceptibility reading cannot establish spontaneous order or hysteresis. The original curves below are illustrative responses, not measured material data.
Figure 1.
Magnetic response and hysteresis
. Details.
Choose the guide axis as , with and . The electric field is along , parallel to the walls' dimension. Take a homogeneous lossless nonmagnetic dielectric with relative permittivity , ideal conducting walls, and time dependence . The lowest transverse electric mode of a rectangular waveguide is
At , this tangential electric field vanishes. At the electric field is normal to the wall; a nonzero value is supported by surface electric charge and violates no conducting boundary condition. The field has zero divergence. The associated normal magnetic field also vanishes at each wall, as verified below. The fundamental sine across the wider dimension minimizes the nonzero transverse wave number; is not a physical propagating mode of a hollow simply connected guide.
Below the waveguide cutoff frequency , put with to select decay for . The evanescent transverse electric mode in a rectangular waveguide then has
For the fixed frequency in the question,
The decay length increases monotonically and diverges as ; its formal limit at is . At the mode is at cutoff, with zero axial wave number. For , is real and the mode propagates without axial attenuation in this ideal lossless model. There is no real finite evanescent decay length to continue beyond cutoff. The paper's graph notation denotes the same as its preceding formula.
Figure 1.
Evanescent decay length and propagating axial wave number on opposite sides of cutoff
. Details.
From Faraday's law of induction, . Writing gives
Thus vanishes on the walls and vanishes on the walls, meeting the normal magnetic boundary conditions. The fields also satisfy . For real the physical fields are , , and .
The time-averaged Poynting vector is
All products are imaginary because the electric field and magnetic field are in quadrature. A single decaying evanescent wave stores reactive energy without transmitting average power into a semi-infinite lossless guide. A finite below-cutoff section can transmit power when both independent evanescent solutions are present; that is a different boundary-value problem.
The skin effect also gives exponential field decay, but in a good metal has both real and imaginary parts and there is a real inward mean Poynting vector. That flux supplies Joule heating . Here attenuation comes from geometric cutoff in a lossless region, with no conductivity and no dissipative power sink; exponential decay alone does not establish absorption.
Maxwell's displacement current is the electric-field contribution needed to make Ampère's law consistent with charge conservation. In macroscopic form,
The displacement current density has units . Since , it contains the vacuum field term and, in matter, the polarization-current contribution. It need not represent free charge moving through a material.
Taking the divergence gives . Gauss's law, , then yields the continuity equation . Without the displacement term, Ampère's original law would incorrectly require every current distribution to be divergence-free even while charge accumulates.
A charging parallel-plate capacitor makes this concrete. Ignoring fringing, , so the displacement current across its dielectric gap is , equal to the conduction electric current in its wire. The integral law
therefore gives the same answer whether the spanning surface cuts the wire or passes through the gap. Together with Faraday's law, this field term also permits self-consistent electromagnetic waves in vacuum, with speed .
At every , ,
There is no constant-pressure coexistence interval, unstable loop or finite-temperature critical inflection. The isothermal compressibility is positive:
Thus this equation of state predicts no liquid-gas phase transition in its physical domain. Its divergence at the excluded-volume boundary is a repulsive constraint, not a liquid-gas coexistence mechanism. This conclusion concerns the stated model and cannot rule out an ordering transition requiring additional physics or an additional thermodynamic variable. Extrapolation to absolute zero is also outside the classical gas model's reliable regime.
Integral conservation. Electric charge cannot be created by ordinary transport. For a fixed volume , outward current reduces the enclosed charge:
Here has units and has units . Applying the divergence theorem to an arbitrary volume gives the local continuity equation
It is local charge conservation: a net outflow must be accompanied by depletion. In a steady state, . Integrating over a small circuit junction then gives Kirchhoff's current law. During a transient the junction can accumulate charge and the incoming and outgoing conduction currents need not yet match.
Consistency with Maxwell's equations. In vacuum, take the divergence of the Ampère-Maxwell equation:
Gauss's law, , reproduces the continuity equation. The displacement current term makes these equations consistent for time-dependent charge distributions, such as a charging capacitor. It is not a stream of free charges across the insulating gap. In matter, and give the analogous conservation law for free charge; total charge also includes bound charge.
Origin of diffusion. In a material near local thermal equilibrium, microscopic carriers transfer energy from hotter regions to colder regions. On distances long compared with their mean free path and times long compared with relaxation times, the local constitutive law is Fourier's law,
Combining this temperature-gradient-driven flux with local conservation of energy produces a diffusion equation. This description assumes diffusive thermal conduction; it is not generally valid for ballistic transport, and flow would add advective transport.
Energy balance with heating. Choose any fixed volume , let be its uniform volumetric heat capacity and its heat generation rate per unit volume. The rate of increase of internal energy equals heat generation minus outward heat flux:
Since is arbitrary, for constant thermal conductivity
The ratio is the thermal diffusivity in . Here is in , in and in , so both right-hand terms have units . If thermal conductivity varies with position or temperature, the divergence must be retained rather than replacing it by .
In a homogeneous linear medium, write and . Without free charge or current, Maxwell's equations become
Taking the curl of Faraday's law, using , gives
Their plane-wave solutions have and speed ; in vacuum this is . For a wave , the divergence equations give , while Faraday's law gives . Thus the fields and propagation direction form an orthogonal triad, with .
The two independent transverse electric components allow linear, circular and elliptical polarization of an electromagnetic wave. Their relative amplitudes and phase determine the polarization state. In a lossless nondispersive medium the energy density and energy flux are
A travelling plane wave has equal electric and magnetic energy contributions and time-averaged intensity , where is its wave impedance and denotes peak electric-field amplitude. Reflection and refraction follow by imposing the electromagnetic boundary conditions. Frequency-dependent constitutive response introduces dispersion; conductivity introduces attenuation, so the simple real constant-speed result has a specified domain of validity.
Use the divergence identity for a cross product, followed by Faraday's law of induction and the Ampère-Maxwell equation:
With the electromagnetic energy density of part (b), these give
Integrate over a fixed volume and apply the divergence theorem, choosing outward:
This is Poynting's theorem. The outward electromagnetic energy flux equals minus the power transferred to matter and minus the rate at which electromagnetic energy accumulates inside. For the dissipative electrical conductor in this question, the matter term is Ohmic heating.
In SI units, with microscopic charge density and current density , Maxwell's equations can be written
The first relates electric flux to charge. The second states that the net magnetic flux through a closed surface is zero. Faraday's law relates electric circulation to changing magnetic flux. The Ampère-Maxwell equation includes both conduction current and displacement current; the latter is essential to time-dependent fields and local charge conservation.
For fixed contours and spanning surfaces, the integral curl equations are and . For a moving material contour the first must include the motional force term; its total emf is . This prevents confusing a changing circuit area with a partial time derivative of a static magnetic field.
Taking the divergence of the Ampère-Maxwell equation and using electric Gauss's law gives . In a source-free vacuum, combining the two curl equations gives , and similarly for , predicting electromagnetic waves of speed . The Lorentz force law then specifies how these fields act on a charge. Macroscopic matter descriptions may instead use and with free sources and constitutive relations; those should not be silently mixed with the microscopic-source equations above.
In the northern hemisphere the coast is on the wave's right when looking in its direction of propagation. In the ocean- convention above it travels northward, with land east of the flow; for the opposite ocean side it travels southward. This follows from combining with a sea-level bulge that decreases offshore, rather than choosing propagation independently of trapping.
For a positive northward bulge, is northward and strongest at its crest, and there is no cross-shore velocity in the ideal mode. Ahead of the crest the current decreases northwards, so : horizontal flow convergence increases sea level, . Behind the crest : divergence lowers sea level. At the crest itself the alongshore derivative and instantaneous horizontal flow convergence vanish. These signs reverse appropriately for a trough or the opposite coast, while the trapping rule remains the same.
The first panel below sketches a positive bulge, alongshore velocity and the leading convergence/trailing divergence regions. Its cross-shore decay is the pressure-gradient/Coriolis balance, not a current directed into the coast.
Figure 1.
Coastal Kelvin-wave trapping, convergence and amplitude under conserved energy transport
. Original normalized model. The ocean lies left of the coast, the bulge travels north, and the amplitude curves assume slow changes with negligible dissipation and reflection. Details.
Use a shallow, homogeneous active upper layer of mean thickness and mass density , over a much deeper, effectively stationary layer with mass density . The Boussinesq approximation keeps the small density contrast only in buoyancy, giving reduced gravity . The hydrostatic approximation and shallow-water approximation make horizontal velocity nearly depth independent and turn the interface slope into the horizontal force .
Linearize about rest and uniform , with , neglecting products of perturbations and convective acceleration. The f-plane approximation sets the Coriolis parameter to constant . Neglect wind forcing, viscosity, entrainment and bottom drag. The two horizontal momentum equations then balance local acceleration, the Coriolis force and the pressure gradient. The signs and follow from eastward and northward coordinates.
The thickness equation is the linearized conservation of mass law: horizontal divergence decreases , while horizontal convergence increases it. Here is active-layer thickness, not total ocean depth, and rather than gives the restoring force because the motion principally displaces the thermocline.
Take the curl of the Euler equations for an inviscid fluid. The gradient on the right has zero curl. Use the stated vector identity with , :
incompressible flow gives , while a curl always has zero divergence, so . The remaining terms give
For a two-dimensional velocity independent of , and . Thus the scalar vorticity is conserved following each fluid particle. It can still vary from particle to particle; this is advection, not a claim of spatially uniform vorticity.
Use the two-dimensional streamfunction convention , . Here
The divergence is , so the flow is incompressible.
In vacuum, the source-free Maxwell equations are
Put and denote the constant electric amplitude by . A corresponding traveling-wave magnetic field is
Both divergences vanish because both the electric field and the magnetic field depend only on and have zero component. Direct differentiation gives and . A time-independent background magnetic field satisfying the homogeneous equations could be added; the formula is the wave component.
Introduce the magnetic vector potential and electric potential by and . These identities automatically give the divergence equation for and Faraday's law. For this wave choose and ; the derivatives just calculated verify both potential equations.
The electromagnetic force is the Lorentz force, . On the smooth constraint plane , the constraint force is normal to the plane. Both and are tangential, so is normal and is balanced by that constraint force; it does not affect the tangential motion. Thus
A particular motion is , a circle of radius
The general planar motion also includes a constant drift velocity; the particle moves on a circle when that drift is zero. The constraint force can do no work because it is perpendicular to the velocity.