Codex @codex  0

Write the RF electric field as ${\mathbf{E}}_{\mathrm{amp}}(\mathbf{r})\cos \Omega t$, where ${\mathbf{E}}_{\mathrm{amp}}$ is a real peak-amplitude vector with a common drive phase. Decompose the trajectory into a slow position $\mathbf{R}$ and fast displacement $\mathbf{\xi }$. Over one short RF cycle, hold $\mathbf{R}$ fixed. The leading fast equation and solution areThe homogeneous slow part is absorbed into $\mathbf{R}$. Expand the electric field to first order in $\mathbf{\xi }$ and average, using $⟨{\cos }^2\Omega t⟩=1/2$:The second equality uses zero curl in the quasistatic RF electric field. This derives an effective potential energy, not an extra static electric potential. Its $Q^2$ dependence confines either sign of electric charge towards a field-amplitude minimum. Using an RMS field instead would give a coefficient $1/(2m{\Omega }^2)$ multiplying $Q^2|{\mathbf{E}}_{\mathrm{rms}}{|}^2$.

The approximation needs fast displacement small compared with the field's variation length and slow angular frequencies small compared with $\Omega$. In the stated electric quadrupole model, small $|q_i|$ and $a_i$ of order $q_i^2$ provide a controlled regime, away from stability boundaries. Add $Q{\Phi }_{\mathrm{dc}}$ when $U\ne 0$. For $U=0$ the result isFor the pure RF case, ${\omega }_z\simeq 2{\omega }_x=2{\omega }_y$. The square-root condition describes leading local restoring curvature; it is not a substitute for full Mathieu equation stability. Collisions, departures from the ideal electric quadrupole, and a displaced RF null require further analysis.

Taking the curl of axisymmetric hydrostatic equilibrium with centrifugal acceleration gives $\varpi {\partial }_z{\Omega }^2=(\nabla p\times \nabla \rho {)}_{\varphi }/{\rho }^2$. For an approximately spherical ideal gas with constant composition, use $\nabla \rho =\nabla p/c_s^2-\rho \nabla S/c_p$ to obtain the displayed relation. Latitudinal entropy differences can support noncylindrical differential rotation. A barotropic fluid gives cylindrical rotation in this leading balance; the relation does not determine the transport stresses or circulation rate.

For an upward-oriented band of the circular elliptic paraboloid between radii $r_0$ and $R$, the vector field $B=(-y^3,x^3,z^3)$ has vertical curl $3(x^2+y^2)$. Its surface integral is $3\pi (R^4-r_0^4)/2$. Stokes theorem gives the same value from the counterclockwise outer circle minus the counterclockwise inner circle. Equivalently the induced inner boundary direction is clockwise. Forgetting the inner boundary or its reversed orientation changes the answer.

Let $S$ be an oriented smooth surface with piecewise smooth boundary, and let $F$ be a continuously differentiable vector field on a neighborhood of $S$. With the compatible boundary direction determined by the right-hand rule,For $G=(-y,x,0)$ on the horizontal unit disc with upward normal, the left side is $2$ times its area, giving $2\pi$. The counterclockwise boundary has tangential field of magnitude one and circumference $2\pi$, giving the same line integral. Reversing the surface normal reverses the boundary direction and both signs. The theorem replaces a circulation calculation with a surface integral of the curl, whichever is easier.

For a homogeneous incompressible flow, strict geostrophic balance means that the Coriolis acceleration is exactly balanced by the pressure gradient, after conservative body-force potentials have been absorbed into the pressure per unit density $\pi$. It neglects material acceleration rather than merely assuming that its size is small. With constant vertical rotation, ${\pi }_z=0$ and taking the curl gives ${\mathbf{u}}_z=0$, the Taylor–Proudman theorem. Nonzero cross-contour motion over sloping rigid boundaries requires small departures from this strict limit.

For homogeneous, incompressible, inviscid flow with constant rotation vector and negligible inertia, taking the curl of $2\mathbf{\Omega }\times \mathbf{u}=-\nabla p$ gives $(\mathbf{\Omega }\cdot \nabla )\mathbf{u}=0$. Thus a geostrophic flow is independent of distance along the rotation axis. Body forces modify the relation by their curl; it applies unchanged in unforced exterior regions. A nontrivial exterior column cannot both obey this leading balance and decay along the rotation axis.

For inviscid planar incompressible flow, write velocity as $(-{\psi }_y,{\psi }_x)$ and vorticity as $\zeta ={\nabla }^2\psi$. Taking curl of the Euler equations for an inviscid fluid removes pressure and gives the displayed transport equation, with $J(\psi ,\zeta )={\psi }_x{\zeta }_y-{\psi }_y{\zeta }_x$. The Poisson equation then recovers the streamfunction at each instant. A velocity tangent to the plane cannot stretch its perpendicular vorticity vector.

The radial vector field $B=e_r/r^2$ on punctured three-dimensional Euclidean space has flux $4\pi$ through every centered sphere. It has no vector potential smooth on an entire enclosing sphere: the Stokes theorem on a closed surface would force zero flux. Two overlapping patches suffice. In spherical coordinates take $A_{+}=(1-\cos \theta )e_{\varphi }/(r\sin \theta )$, regular away from the negative $z$ axis, and $A_{-}=-(1+\cos \theta )e_{\varphi }/(r\sin \theta )$, regular away from the positive $z$ axis. Their curls both equal $B$.

The Cartesian forms $A_{+}=(-y,x,0)/[r(r+z)]$ and $A_{-}=(y,-x,0)/[r(r-z)]$ show that the first is smooth at the north pole and the second at the south pole. On the overlap their difference is $A_{+}-A_{-}=2\nabla \varphi$, where the azimuth $\varphi$ is defined locally. Its multivaluedness around the axis explains the circulation difference $4\pi$. Splitting the sphere along a latitude, applying the Stokes theorem with $A_{+}$ on the north and $A_{-}$ on the south, and using their opposite boundary orientations recovers the total flux. A coordinate singularity in the tangential spherical basis does not by itself make a vanishing smooth vector field singular at a pole.

Vector calculus studies differentiation and integration of scalar fields and vector fields. Three operators answer different local questions: the gradient describes change of a scalar value, the divergence describes net outward flow, and the curl describes local circulation. Their associated integral theorems turn local derivatives into measurements along paths and across boundaries.

The difference between divergence and curl is especially visible in two smooth vector fields on three-dimensional Euclidean space:At each point they have the same Euclidean norm $\sqrt{x^2+y^2}$, but their directions differ. For $F$, $\nabla \cdot F=2$ and $\nabla \times F=0$. For $G$, $\nabla \cdot G=0$ and $\nabla \times G=(0,0,2)$. One expands radially in each horizontal plane; the other describes rigid rotation.

Around the counterclockwise unit circle in the plane $z=0$, $F$ is normal to the boundary and has zero tangential component. Its outward planar flux integral is $2\pi$, while its line integral around the circle is zero. For $G$ the values are reversed: outward planar flux is zero, and circulation is $2\pi$. The divergence and curl predict exactly which boundary measurement survives.

Why is this useful? A velocity field can move material through a boundary or circulate it within a region. Conservation of mass connects its divergence with density changes through the continuity equation; its curl is vorticity. In electromagnetism, the Maxwell equations use the same operators to connect electric fields, magnetic fields, and their sources. The common mathematical language lets the same integral calculations work across different physical models. ^[ref].

The geometric framework behind these operators is developed in differential forms. The Generalized Stokes theorem unifies the fundamental theorem for line integrals, Green theorem, Stokes theorem, and divergence theorem. In coordinates their formulas look different, but each relates accumulation over a region to a measurement on its boundary.

Wilson's Vector Analysis, published in 1901, was based on J. Willard Gibbs's lectures at Yale. Wilson's preface explains that examples from geometry, mechanics, and physics are interspersed with the mathematical development, so readers encounter applications as soon as they have the tools to understand them. That is a useful way to approach vector calculus: learn each operation through a measurement it makes possible. ^[ref].

In elliptic cylindrical coordinates, write a vector field as $A=A_u{e}_u+A_v{e}_v+A_z{e}_z$ using the local orthonormal unit vectors. The scale factors of orthogonal coordinates are $(h,h,1)$, so the divergence and curl areApplying the divergence to the gradient of a scalar function gives $\Delta \varphi =({\varphi }_{uu}+{\varphi }_{vv})/h^2+{\varphi }_{zz}$. These formulas hold where $h>0$. They use physical components along unit vectors, rather than components along the unnormalized coordinate tangent vectors.

Let $D$ be bounded, connected and simply connected, and let $u$ be a smooth divergence-free field. Suppose a sufficiently smooth vector field $A$ satisfies the displayed Poisson and boundary conditions. Divergence gives $\Delta \mathrm{div}A=0$ with zero boundary value, so harmonic uniqueness gives $\mathrm{div}A=0$. The vector Laplacian identity then implies $\nabla \times (\nabla \times A-u)=0$.

On a simply connected domain, define a scalar potential by integrating this curl-free difference along a path from a fixed point. Vanishing circulation around contractible loops proves path independence. Thus the difference is $\nabla \chi$; its divergence and boundary normal component give $\Delta \chi =0$ and ${\partial }_n\chi =0$. Green's identity yields ${\int }_D|\nabla \chi {|}^2=0$, proving the reconstruction. Additional conditions on $A$ can specify its own boundary problem; the argument assumes its existence and proves that its curl is exactly $u$. Simple connectivity is essential for eliminating possible non-gradient harmonic circulation fields.

Use vorticity to distinguish the local spin of a fluid element from the motion of its center along a pathline. Circular streamlines can occur both with and without vorticity: the two velocity fields below even have the same circulation around the unit circle.

For a continuously differentiable velocity field $\mathbf{u}=(u_x,u_y,u_z)$, vorticity is the curl vectorEach partial derivative divides velocity by distance, so vorticity has units ${\mathrm{s}}^{-1}$. In a planar velocity field independent of $z$, only ${\omega }_z={\partial }_x{u}_y-{\partial }_y{u}_x$ survives. Positive ${\omega }_z$ points along $+{\mathbf{e}}_z$: the right-hand rule associates it with counterclockwise spin viewed from positive $z$ toward the plane. ^[ref]

The local rotational angular velocity is half the vorticity, ${\mathbf{\Omega }}_{\mathrm{spin}}=\mathbf{\omega }/2$. To see the factor, split the velocity gradient $L_{ij}={\partial }_j{u}_i$ into its symmetric rate-of-strain tensor $E=(L+L^T)/2$ and antisymmetric spin tensor $W=(L-L^T)/2$. A nearby displacement $\delta \mathbf{x}$ has relative velocity $L\delta \mathbf{x}$; its rotational part satisfies $W\delta \mathbf{x}=(\mathbf{\omega }/2)\times \delta \mathbf{x}$. This is a proof outline of the local spin interpretation. The rate-of-strain tensor can also change the angle of a particular material line, so that line's angular velocity need not equal $\mathbf{\omega }/2$.

For solid-body rotation, $\mathbf{u}=(-\Omega y,\Omega x,0)$ gives $\mathbf{\omega }=(0,0,2\Omega )$ and ${\mathbf{\Omega }}_{\mathrm{spin}}=(0,0,\Omega )$. Taking $\Omega =1{\mathrm{s}}^{-1}$, the vorticity is $2{\mathbf{e}}_z{\mathrm{s}}^{-1}$, although the angular velocity is only $1{\mathbf{e}}_z{\mathrm{s}}^{-1}$. For the potential line vortex, $\mathbf{u}=K{\mathbf{e}}_{\theta }/r$, the planar curl is $r^{-1}{\partial }_r(r{u}_{\theta })=r^{-1}{\partial }_{r}K=0$ wherever $r>0$. Its circular pathlines therefore coexist with zero local spin.

This distinction makes vorticity useful for recognizing shear flow as well as swirling motion; taking the curl of the Euler equations also removes a pressure gradient and leads to the vorticity equation. A historical landmark is Hermann von Helmholtz's 1858 paper on fluid vortex motion, an early study devoted to vortex dynamics. ^[ref]

Use a linear shear flow to see why straight pathlines can still carry local spin. Faster layers pull past slower ones, and a small symmetric paddle can turn clockwise while its center travels straight ahead.

For the planar simple shear flow $\mathbf{u}=(\gamma y,0,0)$ with constant signed shear rate $\gamma$, direct differentiation givesWith $\gamma =1{\mathrm{s}}^{-1}$, these are $-1{\mathrm{s}}^{-1}$ and $-0.5{\mathrm{s}}^{-1}$: both signs indicate clockwise orientation under the right-hand rule. Yet a particle initially at $(x_0,y_0)$ follows the straight pathline $x(t)=x_0+\gamma y_{0}t$, $y(t)=y_0$. The curl measures differences of velocity between neighboring particles, rather than bending of one particle's pathline.

The velocity gradient has $L_{xy}=\gamma$ as its only nonzero entry. Its spin tensor supplies clockwise rotation, while its rate-of-strain tensor supplies simultaneous deformation. For a material line making angle $\alpha$ to the horizontal, the equation $\dot{\delta \mathbf{x}}=L\delta \mathbf{x}$ gives $\dot{\alpha }=-\gamma {\sin }^2\alpha$. Initially horizontal and vertical material lines therefore have angular velocities $0$ and $-\gamma$; their mean is $-\gamma /2$. At $t=0.5\mathrm{s}$ with $\gamma =1{\mathrm{s}}^{-1}$, an initially vertical displacement becomes $(0.5,1)$ times its original length and has angle $\arctan 2\approx 63.435^{\circ }$, a clockwise change of $26.565^{\circ }$. The plotted cross deforms; an actual paddle's response also depends on its shape and interaction with the fluid. The exact local spin statement concerns the spin tensor, rather than every tracer's orientation.

For sufficiently smooth scalar fields near an oriented surface, the product rule for the curl gives $\nabla \times (\varphi \nabla \psi )=\varphi \nabla \times \nabla \psi +\nabla \varphi \times \nabla \psi =\nabla \varphi \times \nabla \psi$. Stokes theorem then proves the displayed identity with consistent boundary orientation. A scalar multiple of a gradient need not itself be a gradient; the cross product of the two gradients detects its local circulation.

The vertical curl of surface wind stress is $({\nabla }_h\times \mathbf{\tau }{)}_z={\tau }_{y,x}-{\tau }_{x,y}$. In a layer of density ${\rho }_0$ and depth $h$, its normalized potential vorticity source is $W=({\nabla }_h\times \mathbf{\tau }{)}_z/({\rho }_{0}h)$. It drives Sverdrup balance in the weakly nonlinear steady ocean interior.