Distributed circuit model. A transmission line must be treated as a distributed system when its propagation time is appreciable compared with the signal timescale. Divide a uniform two-conductor line into elements of length $dx$. Each has series resistance $R^{\prime }dx$ and inductance $L^{\prime }dx$, and shunt electrical conductance $G^{\prime }dx$ and capacitance $C^{\prime }dx$. The primes denote quantities per unit length, with units $\mathrm{\Omega }{\mathrm{m}}^{-1}$, $\mathrm{H}{\mathrm{m}}^{-1}$, $\mathrm{S}{\mathrm{m}}^{-1}$ and $\mathrm{F}{\mathrm{m}}^{-1}$, respectively.

Kirchhoff's voltage and current laws give the Telegrapher's equations:For a lossless line, $R^{\prime }=G^{\prime }=0$, differentiation gives $V_{xx}=L^{\prime }{C}^{\prime }{V}_{tt}$ and the same wave equation for $I$. The propagation speed and characteristic impedance areA general lossless signal has $V=V^{+}(t-x/u)+V^{-}(t+x/u)$ and $I=[V^{+}(t-x/u)-V^{-}(t+x/u)]/Z_0$. Thus characteristic impedance is the voltage/current ratio of one travelling wave, not simply the resistance of the wires. It can be real even when the line dissipates no energy, because energy is carried away by propagation.

For sinusoidal quantities with convention $e^{i\omega t}$, the general lossy-line formulas areChoose the square-root branches for forward attenuation and positive forward power. The real part of ${\gamma }_{\mathrm{prop}}$ is attenuation per unit length; its imaginary part is the phase wavenumber. For small losses, attenuation is approximately $R^{\prime }/(2Z_0)+G^{\prime }{Z}_0/2$, using the lossless real $Z_0$ in this estimate.

Fields and geometry. Two separated conductors support a transverse electromagnetic mode, with electric and magnetic fields transverse to the propagation axis. For a homogeneous dielectric and negligible conductor penetration, both examples below have $L^{\prime }{C}^{\prime }=\mu ϵ$ and hence $u=1/\sqrt{\mu ϵ}$. The distributed circuit model follows from the field energy and charges in each slice.

For a coaxial cable with inner radius $a$ and inner surface of the outer conductor at $b$, Gauss's law and Ampère's law giveThe electric field is radial and the magnetic field circles the axis. Integration of voltage and field energy givesThe equal and opposite return current on the outer conductor confines the ideal external field to the annulus.

For wide parallel plates of width $w$ separated by $d\ll w$, neglect fringing. The electric field crosses the gap, $E\simeq V/d$, while the magnetic field lies in the gap parallel to the plates and perpendicular to propagation, $H\simeq I/w$. ThereforeFor a forward wave, $\mathbf{E}\times \mathbf{H}$ points along the line. Integrating this Poynting vector over either cross-section gives instantaneous power $VI$. For complex peak amplitudes, average power is $\mathrm{Re}(V{I}^{\ast })/2$, or $|V^{+}{|}^2/(2Z_0)$ for a single forward wave on a lossless line. The conductors guide an energy flux predominantly in the dielectric between them.

Terminations and reflection. At a load $Z_L$, require $V_L=Z_L{I}_L$. If $V^{+}$ and $V^{-}$ are the incident and reflected voltage amplitudes there, the total current is $(V^{+}-V^{-})/Z_0$. ThusA matched load, $Z_L=Z_0$, absorbs the incident power with no reflection. An open circuit has ${\Gamma }_V=+1$, a current node and a doubled voltage at the end. A short circuit has ${\Gamma }_V=-1$, a voltage node and a doubled current there. A purely reactive load on a lossless line has $|{\Gamma }_V|=1$ and accepts no average power. For real $Z_0$, the reflected power fraction is $|{\Gamma }_V{|}^2$ and the accepted fraction is $1-|{\Gamma }_V{|}^2$.

Incident and reflected waves produce a standing-wave pattern with voltage standing wave ratio $(1+|{\Gamma }_V|)/(1-|{\Gamma }_V|)$. The input impedance of a lossless length $\ell$ isThis gives $i{Z}_0\tan (k\ell )$ for a short and $-i{Z}_0\cot (k\ell )$ for an open. A quarter-wave section gives $Z_{\mathrm{in}}=Z_0^2/Z_L$, enabling impedance transformation. After a source launches a pulse, a distant termination affects it only after the reflected pulse returns; the return time is $2\ell /u$. Real conductor resistance, dielectric loss and departures from a TEM field profile limit the ideal model at high frequency.

Assume all rings have the same circulation $\Gamma$ and orientation, and that there is no imposed background flow. Vortex lines describe the direction of vorticity, not the direction of velocity. For an incompressible fluid, the velocity induced by a thin vortex filament has the same Biot-Savart integral form as the magnetic field of a current filament:A closely spaced coaxial array of vortex rings is therefore analogous to a solenoid. The velocity field passes along the common axis through the rings and returns outside a finite array. Its orientation reverses if every circulation is reversed. In a long array, far from its ends and vortex cores, the velocity field of a dense coaxial vortex-ring array is approximatelyThe exterior field is exactly zero only in the ideal continuous infinitely long limit, with the background velocity set to zero. In a finite array the streamlines spread out at one end, run back outside, and converge into the other. A large number of closely spaced rings alone does not guarantee that the total array length greatly exceeds $r$; if it does not, the end and return flows remain substantial. Different unspecified ring circulations would require superposing their actual signed contributions.

The symmetry is electromagnetic reciprocity applied to mutual inductance. For thin closed circuits in vacuum, the magnetic vector potential generated by circuit B isStokes theorem converts the magnetic flux through circuit A into its magnetic vector potential circulation. Dividing by $I_B$ gives Neumann's mutual-inductance formula,The dot product and separation are unchanged by interchanging A and B, soEach has units $\mathrm{Wb}{\mathrm{A}}^{-1}=\mathrm{H}$. The winding orientations must be used consistently to determine the common sign. The same reciprocity holds in ordinary linear reciprocal magnetic media; it is not an unrestricted statement about biased nonreciprocal materials. Equivalently, a reversible magnetic energy function has equal mixed derivatives with respect to the two currents.

Let $s(t)$ be the length of the free straight segment from the current contact point to the mass, and let $\theta$ measure the contact point's rotation from first contact. For an inextensible taut rope, winding through angle $\theta$ uses length $R\theta$:The point mass follows an involute of the pole's circular cross-section. Its instantaneous velocity is perpendicular to the free rope. To verify the geometry explicitly, use planar polar unit vectors at the contact point:The rope's tension therefore does no work, and the speed stays at its initial value $v$. The geometry gives $v=s\dot{\theta }$. ConsequentlyThe distance from the pole's axis is $\sqrt{R^2+s^2}$, so a point mass first reaches the pole when $s=0$. The winding time around a circular pole is thusThis has units of time. The model assumes a light inextensible rope, no energy loss at contact and a point mass. The required tension $m{v}^2/s$ diverges at the ideal final contact, so finite strength and size matter very close to impact.

Treat the lake as a large thermal reservoir at $T_f=283.15\mathrm{K}$, and take the copper's heat capacity $C=50\mathrm{J}{\mathrm{K}}^{-1}$ constant over the cooling interval from $T_i=783.15\mathrm{K}$. Although the actual cooling is irreversible, entropy is a state function, so its change can be calculated along a hypothetical reversible path:The lake receives $Q=C(T_i-T_f)=2.50\times 10^4\mathrm{J}$ at constant temperature, givingTherefore, neglecting other changes,The positive total is the entropy generated by transferring heat across a finite temperature difference. The copper's entropy decreases, but the colder lake gains a greater amount.

Assume a hollow square electromagnetic waveguide with perfectly conducting walls, filled with vacuum or air. For side $a$, the waveguide cutoff frequency isA transverse electric mode can have one zero index, but not both. A transverse magnetic mode requires both indices positive. Thus the lowest cutoff belongs to the two degenerate transverse electric modes ${\mathrm{TE}}_{10}$ and ${\mathrm{TE}}_{01}$, withThe corresponding cutoff vacuum wavelength is $2a=2.0\mathrm{cm}$. Strictly, a travelling mode carrying power along the guide needs $f>f_c$: at the cutoff itself, the axial wavenumber and group velocity are zero. The number above is therefore the lower frequency threshold, not a propagating mode with nonzero longitudinal energy flux exactly at that frequency. A dielectric filling would replace $c$ by the wave speed in that medium.

Constitutive law and Poisson ratio. Use positive stress for tension. In an isotropic linearly elastic medium, Hooke's law givesThe analogous equations follow by cyclic permutation. Under uniaxial tension in direction 1, $e_1={\tau }_1/E$ and $e_2=e_3=-\sigma {\tau }_1/E$. Thus Poisson's ratio $\sigma =-e_2/e_1$ measures transverse contraction relative to longitudinal extension when the transverse stresses are zero. It is dimensionless; it is not the ratio of transverse stress to longitudinal stress.

Radial and hoop strains. Let a material point originally at radius $r$ move to $r+R(r)$, with small axisymmetric displacement. Two neighbouring radii separated by $dr$ move apart by $dr+dR$, so the radial fractional extension is $dR/dr$. A circle's circumference changes from $2\pi r$ to $2\pi (r+R)$, so its fractional extension is $R/r$. ConsequentlyThe absence of longitudinal stress is a plane stress condition, not a claim of zero longitudinal strain.

Equilibrium of an annular wedge. Consider a slice of unit axial length, radial thickness $dr$ and small angular width $d\varphi$. The difference of the outward and inward radial face forces is $(r{\tau }_r{)}^{\prime }drd\varphi$ to leading order. The two tangential face tensions have inward radial resultant ${\tau }_{\varphi }drd\varphi$: each contributes ${\tau }_{\varphi }dr\sin (d\varphi /2)$. With no body force, their balance gives the radial equilibrium of a thick elastic cylinder,This hoop contribution is why the radial stress need not be spatially constant in a thick wall.

Displacement equation. With ${\tau }_z=0$, the two Hooke's law equations readInvert them:Substitute in ${\tau }_{\varphi }=(r{\tau }_r{)}^{\prime }$. After cancelling the common factor, the right-hand side is $r{R}^{\prime \prime }+(1+\sigma )R^{\prime }$, soThe cancellation explains why this displacement equation does not depend on the value of Poisson's ratio.

Power laws and wall pressure. Try $R\propto r^n$. The indicial equation is $n(n-1)+n-1=n^2-1=0$, soThe resulting Lamé stress distribution in a pressurized cylindrical wall can be writtenThe pressure boundary tractions are ${\tau }_r(a)=-p_1$ and ${\tau }_r(b)=-p_0$. Solving them yieldsHence the positive compressive radial pressure in the wall isIt equals $p_1$ at $a$ and $p_0$ at $b$. This directional radial pressure is not an isotropic fluid pressure: radial and tangential stresses differ. The associated displacement is $R=[(1-\sigma )Pr+(1+\sigma )Q/r]/E$. Also $e_z=-\sigma ({\tau }_r+{\tau }_{\varphi })/E=-2\sigma P/E$, generally nonzero even though ${\tau }_z=0$.

Maximum absolute tangential stress. For $p_1>p_0$, $Q>0$ and ${\tau }_{\varphi }^{\prime }=-2Q/r^3<0$. Therefore the maximum of its absolute value must be at one of the two surfaces. Keeping the specified external pressure gives the fully general answerFor the usual high internal pressure with negligible external pressure, the largest magnitude is the tensile stress on the inner wall:For finite $p_0$, the inner-wall expression is $[p_1(a^2+b^2)-2p_0{b}^2]/(b^2-a^2)$, but its absolute value is not always the larger endpoint value. For example, nearly equal internal and external pressures can leave the outer surface more compressive. With $p_1=p_0=p$, both stresses correctly reduce to $-p$ everywhere. These results use a homogeneous linear elastic wall without axial loading; a closed-ended pipe would require a different longitudinal-stress condition.

Normal-mode idea. Near a stable equilibrium, expand the potential energy to quadratic order and the kinetic energy as a quadratic form in generalized velocities. The Euler-Lagrange equations then take the linear form $\mathsf{M}\ddot{\mathbf{x}}+\mathsf{K}\mathbf{x}=0$. A normal mode has fixed displacement ratios, $\mathbf{x}=\mathbf{q}\cos (\omega t+\phi )$, with $\mathsf{K}\mathbf{q}={\omega }^2\mathsf{M}\mathbf{q}$. Real systems admit this description for sufficiently small motion, before higher-order nonlinearities or strong damping matter. Superposition permits a general motion to be assembled from the modes. A free symmetry direction may instead be a zero mode, with constant or linearly changing coordinate.

Lagrangian and bead equations. Take each $x_i$ as an arc-length displacement positive in the same sense around the fixed circular wire. Spring length changes are $x_2-x_1$, $x_3-x_2$ and $x_1-x_3$. Equal equilibrium spring tensions give no net tangential force; their linear energy terms sum to zero. Dropping an irrelevant constant, the Lagrangian isFor example, differentiating the potential with respect to $x_1$ gives $k(2x_1-x_2-x_3)$. The three Euler-Lagrange equations therefore areThe wire supplies normal constraint forces but does no tangential work. These are the equal-mass case of the normal modes of three beads connected by springs on a circular hoop.

Frequencies. The stiffness matrix divided by $k$ isIt annihilates $1$. On the two-dimensional subspace $q_1+q_2+q_3=0$, its action is multiplication by three. Thus the angular frequencies of the normal modes areThe repeated frequency is a degeneracy of normal modes, while the zero frequency expresses free rigid rotation of all beads along the wire.

General solution. Two independent vectors spanning the vibrating subspace are ${\mathbf{q}}_A=(1,-1,0{)}^{\mathsf{T}}$ and ${\mathbf{q}}_B=(1,0,-1{)}^{\mathsf{T}}$. Each amplitude satisfies $\ddot{Q}+{\omega }^{2}Q=0$, with $\omega =\sqrt{3k/m}$. The uniform coordinate satisfies ${\ddot{Q}}_0=0$. HenceThe six freely specifiable real constants $A,\varphi ,B,\psi ,C,D$ give the six initial data; $\omega$ is fixed by $k,m$ and is not an additional independent integration constant. The linear zero-mode term must be included to allow nonzero total tangential momentum.

Shapes and uniqueness. In the first term beads 1 and 2 have opposite equal displacements and bead 3 remains at its equilibrium position. In the second term beads 1 and 3 do this while bead 2 stays fixed. The last term moves all three identically, leaving spring extensions unchanged: $D$ is a rigid phase shift and $C$ a uniform drift around the wire. The sketches show positive modal amplitudes; all arrows reverse half a vibration period later.

The two nonzero-frequency eigenvectors are not unique, because any independent linear combinations in their degenerate two-dimensional eigenspace form an equally valid basis. In particular the displayed pair is independent but not orthogonal; one could use $(1,-1,0)$ and $(1,1,-2)$ instead. The zero eigenspace is one-dimensional, so its vector is fixed up to scale and sign.

Specified release. Zero initial velocities permit $\varphi =\psi =0$ and require $C=0$. At $t=0$, the three displacement conditions become $A+B+D=a$, $-A+D=0$, $-B+D=0$. Therefore one convenient set of constants isEquivalently,Their sum remains $a$, as required by zero initial total tangential momentum. Alternative amplitude signs with compensating phase shifts represent the same motion.

Radial energy and effective potential. For a central force with potential energy $U(r)$, angular momentum conservation fixes $J=m{r}^2\dot{\varphi }$ and confines motion to a plane. The kinetic energy separates into radial and angular parts, givingThe second term is the centrifugal barrier. It converts the radial motion to a one-dimensional energy problem: allowed radii have $E\ge U_{\mathrm{eff}}$, turning points have equality, and a circular orbit occurs at a stationary point of $U_{\mathrm{eff}}$. Stability of that radial orbit requires a local minimum.

The first-order orbit equation. Write $A=GMm$. Since $f=-dU/dr=-A/r^2-A\gamma /r^3$, choose zero potential at infinity to getAlso $\dot{r}=(dr/d\varphi )J/(m{r}^2)$. Substitute this into the energy equation and multiply by $2m{r}^4/J^2$:Here $\gamma$ has units of length. A regular bound orbit with a centrifugal barrier requires $J^2>mA\gamma$ as well as the appropriate negative energy.

Solving in inverse radius. Let $u=1/r$ and defineDividing the first-order equation by $r^4$ givesDifferentiating, or directly using the Binet equation, gives $u^{\prime \prime }+qu=b$. The continuation through turning points follows from this smooth second-order equation. Its noncircular bound solution, with angular origin at a pericentre, isSubstitution into the first-order equation fixes the amplitude through this $e$. Therefore the requested representation holds withsince ${\alpha }_1/r-{\alpha }_2=\cos ({\alpha }_3\varphi )$. A general angular origin replaces $\varphi$ by $\varphi -{\varphi }_0$. This is the orbit equation for combined inverse-square and inverse-cube attraction. The circular limit $e=0$ is regular in $u=b/q$ but is not represented by finite ${\alpha }_1,{\alpha }_2$ in the source's normalized cosine form.

Precession and the inferred length. The last expression immediately givesSuccessive pericentres occur after angular advance $2\pi /{\alpha }_3$, so the apsidal precession per orbit isFor ${\alpha }_3$ near one, $1-{\alpha }_3^2\simeq \Delta \varpi /\pi$. Inferring an inverse-cube perturbation from apsidal precession consequently gives $\gamma \simeq [J^2/(mA)]\Delta \varpi /\pi$. For the stipulated nearly circular approximation, $J^2/(mA)\simeq r_0$, and Kepler's third law givesusing $P=88(86400)\mathrm{s}$, $M=2.0\times 10^{30}\mathrm{kg}$ and $G=6.6743\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$. A century contains about $36525/88=415.1$ orbits. Convert the observed advance to radians per orbit before using it:HenceThe positive sign strengthens the attraction and gives prograde precession. The fractional perturbation $\gamma /r_0\simeq 1.6\times 10^{-7}$ validates the small-deviation expansion. This is the parameter fitting the proposed phenomenological force law under the specified small-eccentricity approximation, not a derivation of that law as the physical explanation of Mercury's anomaly.

Changing basis and acceleration. Attach orthonormal unit vectors ${\mathbf{e}}_i(t)$ to the rotating frame, with a common origin fixed in the inertial frame. An infinitesimal rotation $d\mathbf{\theta }=\mathbf{\omega }dt$ changes a vector fixed in those axes by $d\mathbf{A}=d\mathbf{\theta }\times \mathbf{A}$: its tip moves tangentially to the circle about the rotation axis. More generally, write $\mathbf{A}={\sum }_i{A}_i{\mathbf{e}}_i$. Differentiating the components and the rotating basis separately proves the rotating-frame derivative formula,

Apply the formula first to position: ${\dot{\mathbf{r}}}_0=\mathbf{v}+\mathbf{\omega }\times \mathbf{r}$, where $\mathbf{v}=(d\mathbf{r}/dt{)}_S$. Apply it again to both terms. For constant $\mathbf{\omega }$,One cross term comes from differentiating $\mathbf{v}$ in the inertial frame and another from differentiating $\mathbf{r}$ in $\mathbf{\omega }\times \mathbf{r}$. The displayed source formula implicitly requires constant angular velocity and an unaccelerated common origin. For changing rotation rate there is also $\dot{\mathbf{\omega }}\times \mathbf{r}$; an accelerating origin supplies its own translational acceleration. Moving the rotation terms to the force side gives the apparent Coriolis force and centrifugal force with negative signs.

Southward motion on Earth. Use local east, north and up axes, with latitude $\lambda$ positive northward. Then $\mathbf{\Omega }=\Omega (0,\cos \lambda ,\sin \lambda )$ and $\mathbf{v}=(0,-v,0)$. The Coriolis acceleration in local Earth coordinates givesIt is westward in the northern hemisphere, eastward in the southern hemisphere and zero for exactly southward motion at the equator. There is no vertical Coriolis component for this particular velocity.

Correct steering. On the desired meridional track, drag has no eastward component because it is antiparallel to the velocity. If the thrust makes signed angle $\varphi$ east of south, its eastward component is $T\sin \varphi$. The cross-track force balance is thereforesoThe argument must have magnitude at most one for this simple thrust-only correction to be feasible. The signed angle changes sign at the equator: correction is east of south in the north and west of south in the south. Thus the phrase “east of south” cannot describe a positive correction throughout the entire flight.

Wrong-sign correction and the miss distance. A literally constant velocity vector cannot develop a cross-track displacement. The physically useful interpretation of the last task is constant speed, with an erroneous lateral steering correction. In the usual leading-track model, thrust and drag balance along the actual motion, and the applied transverse correction has the correct magnitude but the wrong sign. It therefore reinforces rather than cancels the Coriolis force. Starting at the equator, use the unperturbed southward track $\lambda \simeq -vt/R$, and let $x_E$ be the eastward departure from that track. To first order in the departure,There is a factor two relative to uncompensated Coriolis deflection because the wrong correction adds an equal lateral force. Integration gives the wrong-sign Coriolis steering error,For southern angular distance $\Lambda =16\pi /180=0.27925$, $R=6.4\times 10^6\mathrm{m}$ and $v=800/3.6=222.2\mathrm{m}{\mathrm{s}}^{-1}$, the nominal flight time after the equator isTaking $\Omega \simeq 2\pi /(86400\mathrm{s})$ givesThe small-latitude approximation $\sin (vt/R)\simeq vt/R$ gives ${\ddot{x}}_E\simeq 4\Omega v^{2}t/R$ and $x_E\simeq 2\Omega v^2{t}^3/(3R)$, or $1.95\times 10^5\mathrm{m}$, confirming the scale and direction.

This is an estimate under an explicit steering model, not an exact trajectory forced by the wording. For example, retaining spherical transport of the local axes while still linearizing in the cross-track motion gives, with $s=-\lambda$, $\dot{u}-(v/R)u\tan s=4\Omega v\sin s$. Consequently $u\cos s=2\Omega R{\sin }^{2}s$ and the longitude departure is $\Delta \ell =2\Omega R(\tan \Lambda -\Lambda )/v$. The distance along the destination parallel is then $R\cos \Lambda \Delta \ell \simeq 1.93\times 10^5\mathrm{m}$. Corrections from the changed track and a more detailed control law are not fixed by the problem. If “keeps the velocity constant” instead means the full vector, the pilot would need an additional compensating force and there would be no miss; if thrust is kept at a specified geographical bearing despite drift, its along-track component and lateral drag must also be modelled.

Field equations. In a homogeneous source-free nonmagnetic medium, $\mathbf{D}=ϵ{\epsilon }_0\mathbf{E}$ and $\mathbf{B}={\mu }_0\mathbf{H}$. Maxwell's equations becomeThe first equation follows from $\nabla \cdot \mathbf{D}=0$ for spatially constant nonzero $ϵ$. When the effective relative permittivity is negative below, its charge response is already included in $\mathbf{D}$; it must not also be counted as an extra free source.

Amplitude relations and bulk dispersion. Work with complex phasors proportional to $e^{ikx-Kz-i\omega t}$. For $\mathbf{H}=(0,H_y,0)$, the Ampere-Maxwell law givesso the amplitudes must obeyThey are not arbitrary independent constants: the dispersion relation alone would not make arbitrary amplitudes satisfy the equations. These relations give $ik{E}_x-K{E}_z=0$, so the electric divergence vanishes; the magnetic divergence vanishes because there is no $y$ dependence. The only nonzero component of Faraday's law of induction isSubstitution then givesThis proves that all four field equations hold for the specified transverse magnetic polarization when both the amplitude relations and the dispersion relation are imposed.

Interface matching. Let region 1 be vacuum above the interface, with ${ϵ}_1=1$, $K={\beta }_1>0$; region 2 has ${ϵ}_2=ϵ$, $K={\beta }_2<0$. These signs make the evanescent fields decay as $z\to \pm \infty$. In the absence of free surface current, tangential $H_y$ is continuous, giving a common interface amplitude. Tangential $E_x$ is also continuous, henceNormal $D_z$ is continuous as well: both sides have $D_z=-k{H}_y/\omega$. This does not require $E_z$ itself to be continuous; the material can have polarization surface charge. The normal magnetic field is zero on both sides.

Existence and surface dispersion. The opposite signs of ${\beta }_1$ and ${\beta }_2$ require $ϵ<0$. The two bulk equations areUse ${\beta }_2^2={ϵ}^2{\beta }_1^2$ to eliminate the decay constants. For the nontrivial interface branch,A positive ${\beta }_1^2$ together with real propagating $k$ gives $ϵ<-1$. A negative permittivity between $-1$ and zero is insufficient. This lossless planar surface-plasmon dispersion describes a surface plasmon polariton; $ϵ=1$ would instead be identical media and supplies no decaying bound mode.

Frequency range and long wavelengths. With $ϵ=1-{\omega }_p^2/{\omega }^2$, the existence condition becomesThe upper value is a limiting frequency approached as $k\to \infty$, not an attained maximum at finite $k$. Substitution into the dispersion relation givesFor $\omega \ll {\omega }_p$, choosing propagation in the positive $x$ direction,It approaches the vacuum light line from larger $k$. The vacuum decay constant is then ${\beta }_1\simeq {\omega }^2/(c{\omega }_p)$, so confinement becomes weak at low frequency. This calculation assumes the stated real, lossless effective permittivity; dissipation would require complex propagation and decay constants.

Field law and centre value. The Biot-Savart law for a steady thin wire isHere $d{\mathbf{\ell }}^{\prime }$ follows the conventional electric current and the displacement points from source to observation point. For a circular loop of radius $a$, every element at the centre has distance $a$ and contributes in the same axial direction. Choose $+z$ by the right-hand rule. Then

On-axis integration. At axial distance $z$, every wire element has distance $\sqrt{a^2+z^2}$. Azimuthal symmetry cancels transverse components. Parameterizing ${\mathbf{r}}^{\prime }=(a\cos \phi ,a\sin \phi ,0)$ gives axial cross-product element $a^{2}d\phi$, so the on-axis magnetic field of a circular current loop isIt points the same way on the two halves of the axis; its strength is even in $z$.

Directions near the axial point. At a point on the positive axis the field is axial, but nearby lines spread outward because $B_0$ decreases as $z$ increases. On the negative-axis side they converge toward the axis while still progressing in the positive axial direction. There is no azimuthal field. The near-axis sketch shows this local spreading; the complete magnetic field lines return around the outside of the loop and form closed curves.

Radial component from zero divergence. Smooth axial symmetry gives $B_z(\rho ,z)=B_0(z)+O({\rho }^2)$ and $B_{\rho }=O(\rho )$. By Gauss's law for magnetism,Integrate radially from zero, excluding a singular $1/\rho$ term, to obtain the near-axis radial magnetic field, $B_{\rho }=-\rho B_0^{\prime }(z)/2+O({\rho }^3)$. Differentiating the axial field therefore givesFor $z>0$ it is outward; for $z<0$ it is inward. This supplies the direction reasoning used in the sketch.