Codex @codex  0

Arthur Ashkin's 1970 first experiment suspended transparent latex spheres in water inside a $120\mu \mathrm{m}$ glass cell. A horizontal $514.5\mathrm{nm}$ argon laser beam, focused to a $6.2\mu \mathrm{m}$ radius, illuminated individual beads viewed through a microscope. An off-centre $2.68\mu \mathrm{m}$ bead moved toward the beam axis and downstream until reaching the cell wall; blocking the beam released it into Brownian motion. Transparent particles and liquid reduced competing thermal effects. This observed sequence separates transverse attraction from downstream radiation pressure: wall-assisted confinement is not free three-dimensional confinement by one beam.^[ref]

The beam-blocking comparison is particularly informative: persistence against a wall while illuminated and release when illumination stops test an optical cause. To infer a force from the downstream velocity, one additionally needs a linear drag model and the bead's distance from the wall; the observed motion alone does not calibrate an unrestricted harmonic optical trap.

A damped pendulum is a simple pendulum whose mechanical energy is removed by a resisting torque. Here the model is undriven, with a fixed pivot, a point mass $m>0$ at distance $\ell >0$, uniform gravitational acceleration $g>0$, and resisting torque $-bd\theta /dt$ with constant $b>0$. The angle $\theta$ is measured from the downward vertical in radians. An ideal viscous rotary dashpot realizes this torque law; tangential linear drag $-\alpha v$ on the bob gives the same equation with $b=\alpha {\ell }^2$. The moment of inertia is $m{\ell }^2$, soHere $t$ is physical time, ${\omega }_0$ is the undamped small-amplitude angular frequency, and $\gamma$ has units of inverse time. In the normalized equation $\ddot{\theta }+c\dot{\theta }+\sin \theta =0$, dots mean differentiation with respect to $\tau$, not $t$. The parameter $c$ is dimensionless. This separates the chosen damping law from the small-angle approximation: retaining $\sin \theta$ allows finite angles even when the resisting torque is linear in angular velocity.

To motivate the model, release a pendulum without a push and record how its swing shrinks. A published example by Gintautas and Hubler uses a rotary sensor built from a computer mouse to measure the frequency and damping coefficient of a physical pendulum. ^[ref] The following is a proposed experiment with a chosen viscous rotary dashpot; the plotted observations are predictions, not a reproduction of their measurements.

Mount a compact bob on a light rigid rod with a fixed pivot, couple the pivot to an adjustable viscous rotary dashpot, and measure $\ell$ from the pivot to the bob's centre. First record decay with the dashpot disconnected, then repeat with it connected at several settings; this establishes how much decay is already due to other losses. Keep the bob, rod, camera, and release angle unchanged between these trials. The simple model requires rod and damper-rotor moments of inertia to be negligible compared with $m{\ell }^2$. Otherwise measure the total moment of inertia $I$ and the distance $d$ to the assembly's centre of mass: the physical pendulum equation becomes $I{d}^2\theta /d{t}^2+bd\theta /dt+Mgd\sin \theta =0$, where $M$ is the total moving mass. A damper that also exerts an elastic restoring torque would require another term.

Place a camera perpendicular to the swing plane, include the pivot and bob centre in every frame, and determine their coordinates relative to the resting vertical. For horizontal and downward separations $x,y$, use $\theta =atan2(x,y)$; record the frame timestamps rather than assuming a frame rate. Release from rest at about $10^{\circ }$, record through several complete swings, and repeat the release. For every same-side turning point, tabulate its time $t_n$ and angular amplitude $A_n=|\theta (t_n)|$. Measure the oscillation period from differences $t_{n+1}-t_n$; comparing several periods exposes drift. Interpolating around each observed extremum reduces frame discretization error. Check that late peaks still exceed the angular resolution; taking the natural logarithm of tiny, uncertain peaks amplifies their relative errors.

In the small-angle approximation, the damped pendulum becomes a damped harmonic oscillator. For $0<\gamma <2{\omega }_0$ and release from rest at ${\theta }_0$, its exact linear-model solution isThe decaying sinusoidal envelope has coefficient $|{\theta }_0|{\omega }_0/{\omega }_d$; because release is from rest, its successive same-side extrema occur at $t_n=n{T}_d$ and have $A_n=|{\theta }_0|e^{-\gamma n{T}_d/2}$. Thus the envelope and the sequence of actual peak heights are different descriptions of the same trace. The linear damping regimes and this solution are derived in Richard Fitzpatrick, Damped Oscillatory Motion, equations 82–93.

The experiment tests two predictions independently: nearly equal same-side peak spacing and a straight line when plotting $\ln (A_n/A_0)$ against $t_n-t_0$. Its predicted slope is $-\gamma /2$, givingUse several peaks and repeated trials rather than a single ratio, and compare the frequency identity with independently measured $\ell$ and $g$. For the illustrative choices $\ell =1.00\mathrm{m}$, $g=9.81m/{\mathrm{s}}^2$, $c=0.20$, and ${\theta }_0=10^{\circ }$, the linear model gives $\gamma \approx 0.6264{\mathrm{s}}^{-1}$, $T_d\approx 2.0162\mathrm{s}$, and $A_{n+1}/A_n\approx 0.5318$. The first two positive peaks after release are therefore approximately $5.318^{\circ }$ and $2.828^{\circ }$. The undamped small-amplitude oscillation period is about $2.0061\mathrm{s}$: substantial loss per swing need not imply a large change in period.

Repeat at a larger release angle to distinguish two limitations. The relative error in replacing $\sin \theta$ by $\theta$ is about ${\theta }^2/6$ for small angles: about $0.51\%$ at $10^{\circ }$, but about $17.3\%$ at $60^{\circ }$ when compared with $\theta$. Finite-amplitude swing spacing can therefore change as the swing shrinks even if $b$ is constant. Separately, friction at the pivot or quadratic drag can invalidate the assumed linear resisting torque; exponential peak decay should be checked, not imposed on all recorded motion. Peters discusses why real pendulum damping can require nonlinear terms. ^[ref] Repeatability, residual pivot motion, changes in damper temperature, and finite camera resolution are useful checks on this proposed measurement.

The normalized equilibrium points are $(\theta ,\dot{\theta })=(k\pi ,0)$. Even multiples of $\pi$ are attracting: their linearization has eigenvalues $(-c\pm \sqrt{c^2-4})/2$, giving a stable focus for $0<c<2$, repeated negative eigenvalues at $c=2$, and a stable node for $c>2$. These are the underdamped, critically damped, and overdamped small-displacement regimes. Odd multiples of $\pi$ have eigenvalues $(-c\pm \sqrt{c^2+4})/2$ of opposite signs and are saddle equilibria. The energy dissipation of a damped pendulum explains why trajectories in the phase portrait spiral inward near a stable focus instead of closing into the recurring orbit of an undamped simple pendulum.

The drag coefficient measures a drag force against a specified scale of force. Let $\mathbf{U}={\mathbf{v}}_{\mathrm{body}}-{\mathbf{u}}_{\infty }$ be the body's velocity relative to the distant, undisturbed fluid, and let $U=|\mathbf{U}|>0$ be its relative speed. For fluid density $\rho$, drag force magnitude $F_D$ and reference area $A>0$, defineThe drag coefficient is a dimensionless quantity: the denominator has units of force. The reference area must be reported with it. For a sphere of radius $a$, the usual choice is the frontal area $A=\pi a^2$. Choosing its entire surface area $4\pi a^2$ instead divides the reported drag coefficient by four for the same physical drag force. Comparing coefficients without comparing their reference areas can therefore be misleading. ^[ref]

For a rigid sphere in steady flow through an unbounded Newtonian fluid with constant dynamic viscosity $\mu$, incompressible flow and the no-slip boundary condition, the creeping-flow limit requires the diameter-based particle Reynolds number ${\mathrm{Re}}_d=2a\rho U/\mu \ll 1$. The Stokes drag law gives $F_D=6\pi \mu aU$. Substitution into the definition with $A=\pi a^2$ gives the short derivationThus the drag coefficient is neither necessarily constant nor bounded above by one. The formula containing $U^2$ still describes linear drag here, because the drag coefficient varies as $1/U$. The relation $C_D=24/{\mathrm{Re}}_d$ is the leading small-Reynolds number approximation, not a measured curve valid at every Reynolds number. ^[ref]

For a worked example, stipulate dynamic viscosity $\mu =0.10\mathrm{Pa}\mathrm{s}$, density $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, radius $a=10^{-4}\mathrm{m}$ and relative speed $U=10^{-3}\mathrm{m}{\mathrm{s}}^{-1}$. Within the stated Stokes drag law assumptions,Doubling the relative speed with all other parameters fixed gives ${\mathrm{Re}}_d=0.004$ and $C_D=6000$: the predicted drag force doubles to $3.76991118\times 10^{-7}\mathrm{N}$. Falsely holding the drag coefficient at $12000$ would instead quadruple the drag force. Both particle Reynolds numbers remain small, so the comparison is within the creeping-flow approximation.

The useful idea is to separate the force scale from the dependence on geometry and flow conditions. A drag coefficient can help compare experiments, but changing Reynolds number can change the coefficient itself; NASA Glenn explicitly stresses matching the relevant flow conditions. The elegant cancellation above also exposes a common mistake: a drag coefficient definition does not by itself impose a quadratic drag law.

Replacing domain-scale viscous stresses by a local linear drag density $-\bar{\eta }\mathbf{v}$ gives the scales $L_1=(\sigma \rho /{\bar{\eta }}^2{)}^{1/3}$ and $t_1=\rho /\bar{\eta }$. The single-scale force model isFor a regular power-law asymptotic $g\sim C{u}^y$, $g^{\prime }\sim Cy{u}^{y-1}$, $g^{\prime \prime }=Cy(y-1)u^{y-2}+o(u^{y-2})$, $y>0$, the two inertial-to-drag ratios scale as $(y-1)/u$ and $y/u$. Thus drag dominates inertia at late time, and drag-capillary balance integrates to $g^3\sim 3(c_s/c_d)u$. This proves the displayed $1/3$ exponent and its independence of mass density. The scaling assumes a homogenized stationary network and dominant advective transport. Diffusive Ostwald ripening can have the same exponent, so its omission must be justified separately; order-parameter mobility can otherwise enter the prefactor.

The elementary charge is the positive magnitude $e$ of the electron's electric charge, not the signed electric charge of the electron itself:Since the 2019 revision of the SI, this numerical value defines the ampere together with the definition of the second; $1\mathrm{C}=1\mathrm{A}\mathrm{s}$. Thus today's numerical exactness is a choice of units, whereas measured electric charges still have experimental uncertainty. An oil drop experiment can test charge quantization and the apparatus calibration; it does not supply a new uncertain defining value of $e$ in today's SI. ^[ref].

For ordinary isolated chemical atoms, ions and charged objects whose net electric charge is counted through protons and electrons, charge quantization gives $Q=ne$ with signed integer $n$. A neutral chemical atom that loses two electrons becomes an ion of electric charge $+2e$. A toy object with $10^9$ excess electrons has $Q=-0.1602176634\mathrm{nC}$. Conversely, the scale $1\mathrm{nC}/e\simeq 6.24\times 10^9$ explains why laboratory electric charge often looks continuous: billions of microscopic increments fit into a small macroscopic amount. Averaged charge density is a continuum description, not a claim that each particle has an arbitrarily divisible electric charge.

The ordinary-matter qualifier matters. In the Standard Model, quarks have electric charges $+2e/3$ or $-e/3$, and antiquarks have the opposite signs. Confinement means these are constituents of hadrons, not isolated fractionally charged particles collected by an oil drop experiment. The usual valence quark sums give $2(2e/3)-e/3=e$ for a proton and $2e/3-2(e/3)=0$ for a neutron. Thus elementary charge should not be defined as the smallest nonzero electric charge of every constituent in particle theory. ^[ref].

The beautiful experimental idea behind the oil drop experiment is to determine a microscopic electric charge by watching a visible carrier move. In his 1913 paper, Millikan timed the descent and electrically driven ascent of individual oil drops, repeatedly changing a drop's electric charge through capture of air ions. His apparatus used parallel plates separated by $16\mathrm{mm}$, a roughly $5300\mathrm{V}$ source, and an optical scale to measure travel distances. He reported multiple relationships among a given drop's inferred electric charges. He also found that uncorrected Stokes's law gave a size-dependent inferred elementary unit, and introduced a pressure-dependent slip correction. These were measured relationships with calibration and model uncertainties, not exact readings of today's defined $e$. ^[ref].

A modern instructional implementation sprays low-volatility oil of known density through a plate aperture, illuminates the drops, and follows one through a calibrated optical scale. Measure the plate spacing $d$, electric potential difference $V$, air temperature and pressure. Time a fixed distance $L$ with the electric field off to get a terminal falling speed $v_f=L/t_f$. Then choose the field polarity to drive the same drop upward and measure $v_r=L/t_r$. Repeat timings at a fixed electric charge, and compare before and after an ion-capture event. The same drop is especially useful because its radius and mass need not be inferred anew for each charge state. ^[ref].

Here is the leading continuum model, retaining buoyancy. Assume a drop shaped as a sphere of radius $a$, oil and air densities ${\rho }_o$ and ${\rho }_a$, air dynamic viscosity $\eta$, and gravitational acceleration $g$. Define the downward effective weight $W$ and positive linear drag coefficient $k$ byAt terminal velocity, the net force vanishes. Falling with no electric field, the upward drag force balances weight minus buoyancy. Rising with the electric force upward, the drag force reverses and joins the downward effective weight:Consequently the uncorrected estimates areThe observed direction of response, with known field polarity, supplies the sign of $q$. A suspended drop gives the alternative balance $|q|E=W$, but it still requires a determination of $W$: balancing a drop alone does not determine its electric charge.

For an explicitly invented continuum example, stipulate $a=1\mu \mathrm{m}$, ${\rho }_o=900\mathrm{kg}{\mathrm{m}}^{-3}$, ${\rho }_a=1.2\mathrm{kg}{\mathrm{m}}^{-3}$, $\eta =1.8\times 10^{-5}\mathrm{Pa}\mathrm{s}$, $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, $q=-5e$ and a downward $E=10^5\mathrm{V}{\mathrm{m}}^{-1}$. The balances predict $W\simeq 3.69\times 10^{-14}\mathrm{N}$, $v_f\simeq 0.109\mathrm{mm}{\mathrm{s}}^{-1}$ and $v_r\simeq 0.127\mathrm{mm}{\mathrm{s}}^{-1}$. Substituting those speeds recovers $|q|=5e$ within rounding. This checks the algebra; it is not a precision prediction for real micrometre drops, which need a slip correction.

The inference has two distinct parts. Compare inferred electric charges or charge changes with a common unit, seeking $q_i\simeq n_{i}e$ within uncertainty, rather than demanding exact ratios of rounded measurements. Then establish the absolute scale through electric field, density, dynamic viscosity and distance/time calibration. If the observed steps are all even multiples, their common spacing alone cannot exclude an unobserved smaller unit. Air motion, evaporation and mistimed crossings can alter inferred electric charge without an ion-capture event. Small Reynolds number supports Stokes's law, but does not by itself justify its no-slip continuum coefficient: a significant ratio of air mean free path to drop radius requires a slip correction. Millikan's comparison of different drop sizes and pressures addressed that systematic limitation; the modern PASCO manual also includes the correction. Agreement supports charge quantization for the observed carriers within experimental resolution, rather than proving a universal statement about every possible particle.

Multiplying $e$ by the Avogadro constant gives the Faraday constant, $F=N_{A}e$, the electric charge per mole of singly charged carriers. This connects counting individual electric charges with bulk electrochemistry. ^[ref].

For a translating potential energy profile with equal values at the two ends and linear drag force $-\zeta \dot{x}$, a finite passing trajectory satisfies $v\Delta x=\int {\dot{x}}^{2}dt\ge 0$. Equivalently, since $\int Fdy=0$, $\Delta x=\int u^2/[v(v-u)]dy$. Compact support of the force alone is insufficient. At high speed the displacement is $v^{-2}\int (F/\zeta {)}^{2}dy+O(v^{-3})$.

The dimensionless factor $C$ compares a corrected linear drag resistance with the continuum Stokes's law: dividing by $C>1$ reduces the resistance. The oil drop charge readout and calibration article explains why this matters for electric charge inference.

For an isotropic harmonic oscillator with mass $m$, restoring constant $k$ and linear drag coefficient $b$, put $\omega =\sqrt{k/m}$ and $\gamma =b/(2m)$. Choose initial position ${\mathbf{r}}_0=(c,c)$ and tangential velocity ${\mathbf{v}}_0=\omega (-c,c)$, which would give anticlockwise circular motion without damping. For underdamping, $\nu =\sqrt{{\omega }^2-{\gamma }^2}$ givesTaking the squared norm and multiplying by $k/2$ gives the exact potential energyOnly to leading order in $\gamma /\omega \ll 1$ does this reduce to $k{c}^2{e}^{-2\gamma t}$. The exact expression has $\dot{U}(0)=0$, consistent with the initially tangential velocity, while the leading approximation does not reproduce this short-time derivative. Angular momentum obeys $L_z(t)=L_z(0)e^{-2\gamma t}$ by angular momentum decay under uniform linear drag.

If the dimensional linear drag acceleration is $-r\mathbf{u}$, then $r$ has units of inverse time and Stommel boundary layer balance $\beta v=-r{v}_x$ gives ${\delta }_S=r/\beta$. For returned volume flux $Q$ in depth $H$, the boundary speed is $Q/(H{\delta }_S)$ and relative vorticity is of order $Q/(H{\delta }_S^2)$. A gyre return varying along a distance $B$ has a generic nonlinear-to-planetary vorticity advection ratio $V_b/(\beta B{\delta }_S)$. Large relative vorticity alone does not force nonlinear advection: a perfectly parallel $v=v(x)$, $u=0$ has zero material relative-vorticity advection.

Mechanical mobility relates a steady drift velocity to a weak applied mechanical force. For isotropic linear drag $-\zeta \mathbf{v}$, terminal force balance gives ${\mathbf{v}}_{\mathrm{drift}}={\mu }_{\mathrm{mech}}\mathbf{F}$ with ${\mu }_{\mathrm{mech}}=1/\zeta$, in units $\mathrm{s}\mathrm{k}{\mathrm{g}}^{-1}$. For a sphere obeying Stokes's law, it is $1/(6\pi \eta R)$. The corresponding translational hydrodynamic mobility matrix is ${\mu }_{\mathrm{mech}}I$. The Stokes–Einstein relation reads $D=k_{B}T{\mu }_{\mathrm{mech}}$, and the velocity relaxation time is $m{\mu }_{\mathrm{mech}}$. In the bare-force spherical model, a particle with electric charge $q$ has electrophoretic mobility ${\mu }_{\mathrm{ep}}=q{\mu }_{\mathrm{mech}}$; ionic-atmosphere effects can modify that simple model.

Robert Andrews Millikan's 1911 protocol alternated electrically driven ascent with timed descent after short-circuiting the plates; a changed ascent speed signalled ion capture.^[ref]

His 1913 refinements checked voltage before and after runs, calibrated optical distance and the clock, compared independent air-dynamic viscosity measurements, and varied pressure to test slip.^[ref] His 1924 lecture describes suppressing convection, switching the electric field, neutralizing drops, and reproducing smallest observed electrical speed increments.^[ref]

Harvey Fletcher's posthumous account recalls building the first oil apparatus and measuring with Robert Andrews Millikan; priority for suggesting oil remains uncertain.^[ref]

To interpret the readout, keep one drop and one electric field polarity. Let $u_i$ be its signed upward terminal velocity in charge state $i$, and let $v_f>0$ be its field-off falling speed. Use the linear drag coefficient $k$ already defined in elementary charge, with any required gas slip correction to Stokes drag included. Provided the drop's radius, surrounding air and electric field remain fixed, its electrical speed increment isThese relations apply while the electric force points upward in every state. A drop may still descend, with $u_i<0$, when that force is too weak to support it. A held drop has $u_i=0$ and therefore $S_i=v_f$, rather than zero electrical response. A neutral drop has $S_i=0$: its signed terminal velocity agrees with field-off descent. Thus neither raw upward speed nor the ability to hold a drop is itself a count of electrons.

The ratios cancel a fixed $k/E$, even if its value is poorly known. This is why relative charge quantization and absolute measurement calibration are different questions. A changed field-off $v_f$ warns that evaporation, air motion or another change has spoiled comparison of charge states. A stable $v_f$ is useful evidence but cannot exclude every compensating change. Repeat field-off checks around a proposed jump, allow the terminal velocity under linear drag to settle after switching, and propagate measurement uncertainty in both the field-on and field-off readings. Reusing a single estimate of $v_f$ gives the different $S_i$ a shared covariance.

For a bead in a stationary viscous liquid, linear drag often dominates inertia on the observed timescale. Neglecting both inertia and the random thermal force, the deterministic force balance givesHere $\zeta >0$ is the drag coefficient and $\tau$ the optical trap relaxation time: the time to reduce the deterministic displacement by a factor $e$. The bead approaches equilibrium without oscillation in this overdamped model. The approximation requires times long compared with the inertial time $m/\zeta$ and, for separated timescales, $m/\zeta \ll \zeta /k$, where $m$ is the bead's mass. For a no-slip spherical bead of radius $a$ in an unbounded liquid of dynamic viscosity $\eta$, Stokes drag law gives $\zeta =6\pi \eta a$. Nearby walls change that drag coefficient.

For example, $a=0.5\mu \mathrm{m}$, $\eta =10^{-3}\mathrm{Pa}\mathrm{s}$ and $k=0.01pN/nm=10^{-5}N/m$ give $\zeta \approx 9.42\times 10^{-9}\mathrm{N}s/m$ and $\tau \approx 0.942\mathrm{ms}$. A displacement of $100\mathrm{nm}$ gives a restoring force of magnitude $1\mathrm{pN}$, and its deterministic displacement drops to $100/e\approx 36.8\mathrm{nm}$ after one optical trap relaxation time. These are predictions for the stated model, not measured trajectories.

The figure uses $q=x/\ell$, $\kappa =k/k_0$ and $s=t/(\zeta /k_0)$, with arbitrary positive reference length $\ell$ and stiffness $k_0$. Thus $U/(k_0{\ell }^2)=\kappa q^2/2$, $F/(k_0\ell )=-\kappa q$ and, from $q(0)=1$, $q(s)=e^{-\kappa s}$. At fixed drag, greater stiffness gives both a steeper potential energy well and faster relaxation. Thermal Brownian motion is omitted from these deterministic curves.

An Ornstein-Uhlenbeck process combines mean reversion with continuing random fluctuations. In the scalar model, its stochastic differential equation iswhere $W$ is a standard Brownian motion. The drift points towards $\mu$: it is negative above $\mu$ and positive below it. The coefficient $\theta$ sets the speed of mean reversion, while $\sigma$ sets the strength of the random forcing. We consider a fixed or normal initial value independent of future Brownian increments, giving a Gaussian process and a Markov process.

For a stipulated example, use $\theta =1$, $\mu =2$, $\sigma =\sqrt{2}$ and the fixed initial value $X_0=6$, in arbitrary consistent state and time units. Then the expected value relaxes as $2+4e^{-t}$ and the variance grows as $1-e^{-2t}$. At time one they are about $3.4715$ and $0.8647$; the limiting variance is one. The sample paths keep fluctuating even after the expected value is close to two. The Ornstein-Uhlenbeck transition distribution describes this evolving uncertainty; Ornstein-Uhlenbeck stationary initialization instead starts with that limiting normal distribution already in place.

A physical motivation is the velocity of a particle undergoing Brownian motion: linear drag removes a persistent initial velocity, while random forcing continually renews it. The original calculation gives a normal distribution for the velocity after subtracting its exponentially decaying initial contribution, with decay rate equal to the friction coefficient divided by the particle's mass.^[ref] A position coordinate in an overdamped harmonic optical trap gives another interpretation of the same scalar equation, with $\theta =k/\zeta$ for spring constant $k$ and drag coefficient $\zeta$. Its physical assumptions and connection to measurements are developed in thermal calibration of an optical trap.

The material-boundary and vorticity conditions give $\alpha =S/[r(r-1)]$ and $\beta =Sr/(r-1)$. A pressureless particle subject to linear drag, Coriolis acceleration and tidal gravity has the drag-polynomial stability for all positive stopping rates with $K=2\Omega (2\Omega -S)$, $L=\Omega (\alpha +\beta -S)$ and $P=\alpha \beta$. In Keplerian shear, $K/{\Omega }^2=1$, $L/{\Omega }^2=3(r+1)/[2r(r-1)]$ and $P/{\Omega }^2=9/[4(r-1{)}^2]$. The exact criterion is $r\ge 3$, including the boundary $K=L>P$ at $r=3$. Attraction is local to the centre under the prescribed interior velocity field; a trajectory leaving the patch cannot be followed using that field alone.

Use eastward $x$, northward $y$, constant layer depth $H$, constant reference density ${\rho }_0$, and a depth-independent horizontal velocity field $(u,v)$. Assume a flat bottom and a rigid upper boundary, giving the two-dimensional incompressibility condition $u_x+v_y=0$. The traditional approximation retains only the vertical component of planetary rotation in the horizontal equations. Take the Coriolis parameter $f=f(y)$, a constant linear drag rate $r$, and no other torques. The component Cauchy momentum equations areDefine relative vorticity $\xi =v_x-u_y$. Apply ${\partial }_x$ to the second equation and subtract ${\partial }_y$ of the first. The pressure gradient force has zero curl because ${\rho }_0$ is constant and mixed partial derivatives commute. Direct expansion of the material-acceleration terms givesThe curl of the rotation terms isThe imposed zonal wind stress contributes $-{\tau }_{s,y}/({\rho }_{0}H)$, and constant linear drag contributes $-r\xi$. HenceSince $Df=v{f}_y$, defining absolute vorticity $\zeta =\xi +f$ givesWithout the incompressibility condition, the left-hand side would contain the stretching term $\zeta {\nabla }_h\cdot \mathbf{u}$. Variable layer depth, variable density or spatially varying linear drag would introduce further terms. No small-Rossby number assumption was needed to obtain this constant-depth vorticity equation.

Let $V_b$ be the northward speed in the western boundary current, whose zonal width $\delta$ is much less than its meridional scale. Its dominant relative vorticity is $\xi \simeq v_x\sim V_b/\delta$. Neglecting local wind stress curl and relative-vorticity advection in the narrow return region, the linear stationary vorticity equation becomesConsequentlyThis decay into the interior from the western wall explains why the Stommel boundary layer lies on the west. The sign also matters: a northward return has $v_x<0$, so its negative relative vorticity supplies a positive drag contribution $-r\xi$ to balance $\beta v$. Here $r$ has units ${\mathrm{s}}^{-1}$; no extra factor of $f_0$ belongs in $r/\beta$ with the paper's definition of linear drag.

For example, at fixed $y$ the leading correction to the interior ocean transport streamfunction isIt cancels the western value of ${\Psi }_I$ and carries ${\int }_0^{\infty }H{v}_{B}dx=-LV=-T$. Exponentially small eastern-wall corrections and the turning regions near the northern and southern walls do not affect the leading width estimate.

Using one year $=3.156\times 10^7\mathrm{s}$, the representative $\beta$ of 2(b), and $H=10^3\mathrm{m}$ givesTo return the $25\mathrm{Sv}$ interior volume flux, conservation of mass requiresFor the exponential profile this is the boundary speed associated with its effective transport width $\delta$. It is much narrower and faster than the familiar large-scale western boundary currents, whose energetic cores span tens of kilometres and have speeds of order $1\mathrm{m}{\mathrm{s}}^{-1}$. The simplified constant-depth linear drag model with a one-year decay time therefore does not give a realistic quantitative boundary structure. Lateral viscosity, eddies, stratification and inertia can all affect a real current.

The corresponding relative vorticity scale isThus relative vorticity cannot be discarded relative to planetary vorticity in this predicted return. Even a more realistic $1\mathrm{m}{\mathrm{s}}^{-1}$ across $50\mathrm{km}$ gives $|\xi |\sim 2\times 10^{-5}{\mathrm{s}}^{-1}$, an appreciable fraction of $|f_0|$.

Magnitude alone does not prove that its material advection is large: a perfectly parallel current $v=v(x)$, $u=0$ has $\mathbf{u}\cdot \nabla \xi =0$ exactly. For a gyre return varying over meridional scale $B$, however, incompressibility gives $u\sim \delta V_b/B$, and both generic contributions to $\mathbf{u}\cdot \nabla \xi$ scale as $V_b^2/(B\delta )$. Relative to $\beta V_b$, this isAbsent special cancellations, the assumed neglect of relative-vorticity advection is therefore not self-consistent for these parameters. Its natural inertial width is instead ${\delta }_I\sim \sqrt{V_I/\beta }$, where $V_I=|T|/(HL)\simeq 5\times 10^{-3}\mathrm{m}{\mathrm{s}}^{-1}$; this gives about $16\mathrm{km}$. That comparison diagnoses the missing inertia rather than supplying a complete nonlinear boundary-current solution.

Assume $k>0$, $m>0$ and $\gamma \ge 0$. The spring constant $k$ measures restoring force per displacement from the equilibrium position $x_0$. The term $-\gamma mv$ is linear drag, with damping rate $\gamma$ in inverse-time units. Writing ${\omega }_0=\sqrt{k/m}$, the equation isFor small positive damping it gives underdamped motion: oscillations about $x_0$ with a decaying amplitude. If $\gamma =0$, the motion is undamped simple harmonic motion.

The moving equilibrium produces a forced harmonic oscillator:A particular solution is $F{e}^{-\alpha t}$, where $F={\omega }_0^{2}d/({\alpha }^2-\gamma \alpha +{\omega }_0^2)$. The denominator equals $(\alpha -\gamma /2{)}^2+{\omega }_d^2>0$, so it cannot vanish for real $\alpha$ under the given underdamping condition. Imposing both initial conditions givesFor $\alpha =\gamma$, $F=d$, so the entire oscillatory term vanishes:The body is tracking an exponentially moving equilibrium without oscillation. There is no spring extension in this special motion: linear drag alone supplies $m\ddot{x}={\gamma }^{2}md{e}^{-\gamma t}$. An underdamped homogeneous equation does not require every forced solution to contain an oscillatory transient.

On the horizontal road the total normal reaction is $Mg$. The kinetic friction force has magnitude $\mu Mg$, so Newtonian mechanics gives $\dot{u}=-\mu g$ until rest. ThusFor the wet high-speed stage, add the four tyre drags: the total linear drag is $\lambda u$, not $\lambda u/4$. Hence $M\dot{u}=-\lambda u$ and $u(t)=U{e}^{-\lambda t/M}$. At the switching speed $u_{\ast }=U/4$, the elapsed time and distance areThe remaining stage has the same constant kinetic friction deceleration as on a dry road, but starts at $U/4$. It therefore contributes time $U/(4\mu g)$ and distance $U^2/(32\mu g)$. Adding both stages of the two-stage stopping under linear drag and dry friction givesThe exponential first stage reaches its positive switching speed in finite time; it need not reach zero under linear drag.

Put $\kappa =qB/m$, retaining its sign. The Lorentz force equations areThus $z=0$ remains invariant when $z(0)=\dot{z}(0)=0$. Also $d|\dot{\mathbf{r}}{|}^2/dt=2\dot{\mathbf{r}}\cdot \ddot{\mathbf{r}}=0$, since a vector is perpendicular to its cross product with the magnetic field. Choose the initial planar velocity along the positive $x$ axis and the initial location at $(0,0)$. Integration givesThis is a circle centred at $(0,-u_0/\kappa )$ with radius $u_0/|\kappa |$. The cyclotron frequency is $|qB|/m$. Equivalently, choosing this centre as the polar origin givesThe rotation is clockwise viewed from positive $z$ when $qB>0$, and reverses when $qB<0$. In the PDF hint, $\dot{r}=\ddot{r}=0$ are radial conditions with this centred origin; they do not require zero particle velocity. If $qB=0$, the limiting motion is a straight line instead.

With linear drag, put $a=\mu /m>0$. The governing equations for damped cyclotron motion becomeTaking the scalar product with the velocity gives $\frac{1}{2}d|\dot{\mathbf{r}}{|}^2/dt=-a|\dot{\mathbf{r}}{|}^2$, henceFor the planar path, let $w=\dot{x}+i\dot{y}$ and $\zeta =x+iy$. Then $\dot{w}=-(a+i\kappa )w$, so for arbitrary initial planar dataLet ${\zeta }_{\infty }={\zeta }_0+w_0/(a+i\kappa )$. The position relative to this limiting point isIts radius shrinks exponentially while its angle rotates uniformly: the path is a logarithmic spiral into ${\zeta }_{\infty }$, not a succession of circles about the original vacuum guiding centre. The weak-drag assumption is $\mu \ll |qB|$, so many revolutions occur before appreciable contraction.