Codex @codex  0

False. Adiabatic invariance of an orbital action preserves actions under sufficiently slow regular evolution, not the ratios of their frequencies. With action-angle variables, ${\Omega }_i=\partial H/\partial J_i$ can change even when $J_i$ does not. A closed bound orbit requires commensurate radial and angular frequencies; rational commensurability is not an adiabatic invariant.

For an explicit counterexample, slowly change the physical potential from Kepler form towith a small, slowly varying positive $\alpha$. The radial epicyclic frequency and circular angular frequency at a nearly circular orbit areInitially their ratio is one. For generic nonzero $\alpha$, $\Omega /\kappa$ is irrational, producing the Kepler-harmonic epicyclic precession. An initially slightly eccentric ellipse with nonzero radial action therefore becomes a precessing nonclosed orbit, while its radial action remains approximately fixed and its angular momentum remains exactly fixed by spherical symmetry.

Adiabatic actions do not preserve orbital closure. Exceptional circular or resonantly commensurate final orbits do not rescue a claim about every initially closed orbit.

For a normal mode, ${\partial }_t\mapsto -i\omega$, ${\partial }_x\mapsto ik$, and ${\partial }_x^2\mapsto -k^2$. Substitution into the linear complex Ginzburg-Landau equation givesTherefore the dispersion relation isFor real $k$, its temporal growth rate is ${\omega }_i=\mu -k^2$, with unique maximum $\mu$ at $k=0$. Its real angular frequency is $Uk+c_d{k}^2$, displaying the separate advection and wave dispersion terms.

Use the geometrical optics conventionBecause ${\partial }_t=\epsilon {\partial }_T$ and ${\partial }_x=\epsilon {\partial }_X$, differentiating $e^{i\theta /\epsilon }$ gives an order-one physical angular frequency and order-one physical wavenumbers. Absorb a constant leading amplitude phase into $\theta$ and take $A_0$ real; for a complex amplitude the corresponding quadratic density is $|A_0{|}^2$.

Substitution into the wave equation gives the eikonal equation for a variable-speed wave at leading order:Choose the positive-frequency branch $\omega =\Omega =cK$, $K=\sqrt{k^2+{\ell }^2}$; the negative branch can be treated in exactly the same way. The imaginary first-order terms giveMultiplying by $A_0$ puts this wave-action transport for a variable-speed wave in conservation form:Its density is $\omega A_0^2$ and its flux is that density times the group velocity. Independently, equality of mixed wave phase derivatives givesThis phase compatibility in geometrical optics is essential when converting the differentiated dispersion relation into equations along rays.

For completeness let $\omega =\Omega (\mathbf{k};\mathbf{X},T)$ be a general local dispersion relation. Define ray tracing by $d\mathbf{X}/dT={\mathbf{c}}_g={\nabla }_{\mathbf{k}}\Omega$. Using the chain rule in ${\nabla }_{\mathbf{X}}\omega$ and the symmetry ${\partial }_{X_i}{k}_j={\partial }_{X_j}{k}_i$ givesLikewise ${\omega }_T={\Omega }_T+{\mathbf{c}}_g\cdot {\mathbf{k}}_T$ cancels its advective term because ${\mathbf{k}}_T=-\nabla \omega$. Hence the Hamiltonian ray equations for a local dispersion relation areThe derivatives on the right hold the wavevector fixed. For this wave equation they reduce to

When $c=c(Y)$, the ray invariants are $k$ and $\omega$, while $d\ell /dT=-K{c}_Y<0$. For a monochromatic component with these fixed invariants write $\theta =kX-\omega T+S(Y)$, soUsing the stated independence of $A_0$ from $X,T$, the wave-action transport for a variable-speed wave becomes $(c^2\ell A_0^2{)}_Y=0$. ThereforeThis is the stationary wave amplitude in a stratified wavespeed; $\mu$ is the constant for that component or ray family.

At the proposed simple turning point of a variable-speed wave, take $c_s=c(Y_s)=\omega /k$ and $c_s^{\prime }=c_Y(Y_s)>0$. On the chosen positive branch this entails $k>0$. DefineAs the ray approaches from below, $\ell \sim [q(Y_s-Y){]}^{1/2}$ andThe amplitude diverges and the vertical group velocity $c^2\ell /\omega$ vanishes. The WKB approximation also fails its slow-variation condition $\epsilon |{\ell }_Y|/{\ell }^2\ll 1$ when $Y_s-Y=O({\epsilon }^{2/3})$. The turning-point enhancement of wave amplitude shows that this divergence is a failure of the outer ray approximation, rather than an actual infinite solution.

To derive the inner wavefield directly, separate the conserved horizontal wave phase and frequency in the original differential equation:Dividing by $c^2$ gives ${\epsilon }^2{\Psi }_{YY}+2{\epsilon }^2(c_Y/c){\Psi }_Y+Q(Y)\Psi =0$. Set $Y-Y_s=\delta \eta$. The second derivative is of size ${\epsilon }^2/{\delta }^2$, whereas $Q\sim -q\delta \eta$, so the Airy scaling at a variable-speed wave turning point isThe first-derivative term is smaller by $O(\delta )$, and higher coefficients in the expansion of $Q$ give the same relative order. Thus the general leading inner profile is a linear combination of the two Airy functions $\mathrm{Ai}(\eta )$ and $\mathrm{Bi}(\eta )$. For decaying Airy continuation and ray reflection, if the physical continuation is bounded and decays into $Y>Y_s$, its growing $\mathrm{Bi}$ coefficient is zero. In that usual evanescent continuation the wavefield has the formWithout a condition in the evanescent region the corresponding $\mathrm{Bi}$ term remains admissible locally; the given existence of a wavefield below the turning point alone does not select its coefficient.

Finally, the outer amplitude at distance $\delta$ is $O({\delta }^{-1/4})=O({\epsilon }^{-1/6})$ for $\mu =O(1)$ and fixed nonzero $q$. The oscillatory Airy turning-point connection formula has exactly the same $(-\eta {)}^{-1/4}$ envelope on the allowed side. Therefore the inner wavefield has magnitude $O({\epsilon }^{-1/6})$. The decaying Airy continuation matches an incident and reflected pair of rays, so a single upward $\ell >0$ branch must be supplemented by its reflected branch near the turning point; no formal coefficient matching is needed to identify this magnitude.

With no electrical resistance or applied electric potential difference, divide the equation by the inductance $L$ to get a simple harmonic oscillator with angular frequency ${\omega }_0=(LC{)}^{-1/2}$. Its general real solution isEach nonzero solution has period $2\pi \sqrt{LC}$. Here ${\omega }_0$ is an angular frequency; the ordinary frequency in cycles per unit time is ${\omega }_0/(2\pi )$.

For nonzero $\omega$, the equilibrium points satisfy $\sin 2\theta =0$. The equilibrium point $\theta =0$ is stable: it puts the bead at the point farthest from the rotation axis. Writing $\theta =\eta$ small in the tangential ordinary differential equation givesThus the small oscillation is simple harmonic motion, with angular frequency $|\omega |$ and ordinary frequency $|\omega |/(2\pi )$.

The endpoints $\theta =\pm \pi /2$ represent the same rim point $O$. A perturbation there instead satisfies $\ddot{\eta }=+{\omega }^2\eta$ and is unstable. Equivalently, the centrifugal potential along the hoop is $-\frac{1}{2}{\omega }^2|\mathbf{r}{|}^2=-2a^2{\omega }^2{\cos }^2\theta$, whose minimum is at zero and maximum at $O$. If $\omega =0$, every bead position with zero relative velocity is an equilibrium point, but perturbing that velocity produces unrestored drift; there is no nonzero restoring frequency.

Consider a horizontal fluid layer under the Boussinesq approximation heated from below, rotating uniformly about its vertical axis. To make the linear formulas definite, take depth one, stress-free impermeable fixed-temperature plates, a horizontally infinite or sufficiently large periodic domain, and negligible centrifugal-buoyancy modifications. Rigid plates change the vertical eigenfunctions and thresholds and require their own boundary value problem; the explicit formulas below use the free-slip model. Let $P=\nu /{\kappa }_T$ be the Prandtl number, $\mathrm{Ta}=(2\Omega d^2/\nu {)}^2$ the Taylor number, and $R$ the Rayleigh number. Time is again measured on the thermal diffusion time.

The linearized rotating Rayleigh-Bénard convection equations about conduction areThe Coriolis force does no work on the fluid, but it couples vertical motion to vertical vorticity and thereby changes the balance of buoyancy and viscous damping. On the plates require $w=\theta =0$ and ${\partial }_z{u}_x={\partial }_z{u}_y=0$. Write $\zeta =(\nabla \times \mathbf{u}{)}_z$. Taking a curl and a double curl removes pressure and givesA normal mode with horizontal wavenumber $k>0$, vertical index $n\ge 1$ and complex growth rate $s$ has $w=W\sin (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$, $\theta =\Theta \sin (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$ and $\zeta =Z\cos (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$. Put $a=k^2+n^2{\pi }^2$. The three equations become $a(s/P+a)W+\sqrt{\mathrm{Ta}}n\pi Z=R{k}^2\Theta$, $(s/P+a)Z=\sqrt{\mathrm{Ta}}n\pi W$, and $(s+a)\Theta =W$. Their determinant gives the rotating-convection growth-rate polynomialThis cubic includes viscous and thermal modes as well as the convective instability. It avoids excluding a root by division by $s/P+a$. Linear stability means every root has negative real part; stationary onset has $s=0$, while overstability has $s=\pm i\omega$ with nonzero angular frequency.

Setting $s=0$ gives the stationary neutral curve of rotating convectionFor oscillatory convection, set $s=i\omega$ and equate real and imaginary parts. Eliminating $R$ yieldsThe oscillatory neutral curve of rotating convection is physically admissible only when ${\omega }^2>0$. Thus $P<1$ and sufficiently rapid rotation are necessary; a formal minimum with ${\omega }^2\le 0$ is not a Hopf bifurcation. Convection begins at the smaller of the stationary minimum and the admissible oscillatory minimum. The angular frequency approaches zero where a fixed-wave-number stationary and oscillatory threshold meet; that is a double-zero limit of the cubic, requiring a different slow-time reduction.

For wavenumber selection in rotating convection, the lowest vertical index is $n=1$. To compare indices, set $k^2=n^2{\pi }^{2}y$ at fixed $y>0$. Both neutral thresholds have a positive term proportional to $n^4$ and a rotation term independent of $n$, while the admissible angular frequency squared decreases with $n$. Their minima therefore cannot improve on $n=1$. Set $t=k^2$, $m={\pi }^2$. Minimizing each neutral curve givesFor the latter equation also check ${\omega }^2>0$, or minimize over the admissible set instead. Without rotation, $k_c=\pi /\sqrt{2}$ and $R_c=27{\pi }^4/4$. At rapid rotation, $k_s\sim (\mathrm{Ta}{\pi }^2/2{)}^{1/6}$ and $R_{s,c}\sim 3(\mathrm{Ta}{\pi }^2/2{)}^{2/3}$: rotation narrows the horizontal convection rolls and raises their threshold. In a finite box only discrete wavenumbers are permitted, so the minimum must be taken over those allowed modes.

The rapid-rotation comparison has $R_{o,c}/R_{s,c}\sim 2P^{4/3}/(1+P{)}^{1/3}$. Its equality gives $8P^4-P-1=0$, with positive root about $0.67660$. Below this value, sufficiently rapid rotation can make oscillatory convection the primary instability; low Prandtl number allows inertial motions to interact with the thermal field before viscosity damps them. This estimate describes the selected neutral-curve comparison in the stated free-slip asymptotic model, rather than claiming every $P<1$ has oscillatory primary onset at every rotation rate.

Linear selection does not determine the nonlinear planform or its stability. A weakly nonlinear expansion with a solvability condition gives amplitude equations. A stationary convection roll amplitude has the Landau amplitude equation $A_T=rA-g|A{|}^{2}A$; a supercritical branch requires $g>0$ and has $|A{|}^2=r/g$, whereas $g<0$ requires higher-order saturation and permits subcritical behaviour. The coefficient depends on the physical parameters and boundary conditions, and should not be assumed positive merely because the linear threshold is known.

A nonzero-frequency Hopf bifurcation supplies oppositely travelling convection roll amplitudes $Z_{+},Z_{-}$. At a fixed orientation, a cubic normal form iswith complex $g,h$ and the understood symmetry exchanging propagation directions. For $r>0$ and $g_r=\mathrm{Re}g>0$, a travelling-roll branch has one nonzero amplitude, of squared modulus $r/g_r$; it is stable to the opposite travelling amplitude when $h_r=\mathrm{Re}h>g_r$. A standing-roll branch has equal squared moduli $r/(g_r+h_r)$, and is stable within this pair when $g_r>h_r$ and $g_r+h_r>0$. The imaginary parts shift angular frequencies. These conditions concern amplitude perturbations in that reduced subspace; orientation and modulation modes still have to be tested. This is the competition of travelling and standing convection rolls.

Near simultaneous stationary and oscillatory thresholds, a codimension-two bifurcation requires retaining both types of mode. After selecting an orientation and fixing a steady spatial phase, a schematic nonresonant steady–Hopf mode interaction isHere $a$ is real, $z$ is the chosen complex oscillatory amplitude, $g_s,h_s$ are real, and $g_o,h_o$ may be complex. Assume $r_s>0$, $g_s>0$ for a radially stable pure steady branch, and $r_o>0$, $g_{o,r}=\mathrm{Re}(g_o)>0$ for a radially stable pure oscillatory branch. Steady convection rolls suppress the oscillatory mode if $r_o-h_{o,r}{r}_s/g_s<0$, where $h_{o,r}=\mathrm{Re}(h_o)$; oscillatory convection rolls suppress the steady mode if $r_s-h_s{r}_o/g_{o,r}<0$. These are transverse tests for branches whose existence and radial stability have already been checked. Mixed states have positive intensities $u=a^2$, $v=|z{|}^2$ solving $g_{s}u+h_{s}v=r_s$, $h_{o,r}u+g_{o,r}v=r_o$. If $D=g_s{g}_{o,r}-h_s{h}_{o,r}\ne 0$, then $u=(r_s{g}_{o,r}-h_s{r}_o)/D$, $v=(g_s{r}_o-h_{o,r}{r}_s)/D$. For positive $u,v$ and the positive self-saturation coefficients above, the intensity Jacobian matrix has negative trace and determinant $4uvD$, so the mixed state is attracting in intensities when $D>0$ and is a saddle equilibrium when $D<0$; phase and other-mode perturbations remain separate tests. Depending on the cross-couplings, one obtains coexistence or competition and bistability; degeneracies or resonances require extra terms. A full travelling/standing-wave competition must retain both $Z_{\pm }$, not just this illustrative one-mode $z$. If the Hopf angular frequency tends to zero, averaging over fast oscillations fails and a double-zero normal form must keep the corresponding two slow variables. A crossing of global minima at distinct wavenumbers is instead an interaction of distinct modes, not automatically the same double-zero problem.

Weakly nonlinear convection rolls can also lose stability spatially. The real Ginzburg–Landau equation $A_T=rA+\xi A_{XX}-g|A{|}^{2}A$, with $r,\xi ,g>0$, has detuned convection rolls with $|A{|}^2=(r-\xi Q^2)/g$. Linear phase modulation gives diffusion coefficient $\xi (r-3\xi Q^2)/(r-\xi Q^2)$, so the Eckhaus instability excludes $Q^2>r/(3\xi )$ even though convection rolls exist up to $Q^2<r/\xi$. Transverse bending and mean-flow couplings impose additional restrictions in the actual rotating layer; the scalar equation is a local longitudinal example, not its complete stability theory.

Most distinctively, rotation breaks mirror symmetry between competing convection roll orientations. For two sets of convection rolls at relative angle $\vartheta$, write ${\dot{A}}_2=r{A}_2-g_0|A_2{|}^2{A}_2-g(\vartheta )|A_1{|}^2{A}_2$ and interchange indices with $\vartheta \mapsto -\vartheta$ for the other equation. The established convection roll $|A_1{|}^2=r/g_0$ is unstable to the new orientation when its linear growth rate $r[1-g(\vartheta )/g_0]$ is positive. The inequality need not be symmetric under $\vartheta \mapsto -\vartheta$ because the imposed rotation supplies handedness. This is the Küppers–Lortz instability: sufficiently rapid rotation can destabilize steady convection rolls immediately above their stationary onset to oblique convection rolls of another orientation. Successive replacements can produce time-dependent convection roll switching or heteroclinic cycles, instead of a stable single convection roll pattern. It requires three-dimensional perturbations even when the original convection roll solution is independent of its axial coordinate.

At finite Prandtl number and stress-free plates, nearly parallel convection rolls can couple resonantly to a slowly damped large-scale Eulerian mean flow. This small-angle instability of rotating convection rolls makes a regular two-roll expansion nonuniform as the angle tends to zero; one must retain the mean-flow mode. It should not be identified automatically with the finite-angle Küppers–Lortz instability or assigned a universal threshold from a one-amplitude equation. Consequently the onset type and selected scale follow from the admissible neutral curves, while persistent convection roll patterns require a separate nonlinear and sideband stability calculation.

Measure angles from periapsis in the direction of motion. The true anomaly $f$ is the angle at the stellar focus between the pericentre ray and the star-particle radius vector. Draw the auxiliary circle of radius $a$ around the ellipse's center; the point on that circle with the same projected major-axis coordinate as the particle defines the eccentric anomaly $E$ at the center.

The mean anomaly $M=n(t-\tau )$ is a uniformly advancing time phase, with $M=0$ at pericentre passage $t=\tau$. It is the angle swept by a fictitious uniformly moving circular phase marker, not generally the actual position angle. The mean motion is the angular frequencyThe second formula is Kepler's equation, which converts the uniform time phase into geometric position.

Let a slowly varying wave have phase $S(\mathbf{x},t)$, wavevector $\mathbf{k}=\nabla S$, and angular frequency $\omega =-{\partial }_{t}S$. The local dispersion relation is the Hamilton-Jacobi equation $\omega =\Omega (\mathbf{k};\mathbf{x},t)$. Equality of mixed derivatives, followed along a curve with velocity $\dot{\mathbf{x}}={\nabla }_{\mathbf{k}}\Omega$, gives the Hamiltonian ray-tracing equationsHere $d/dt={\partial }_t+\dot{\mathbf{x}}\cdot \nabla$ is the total derivative evaluated along a ray, and $\dot{\mathbf{x}}$ is the group velocity.

For $\Omega =kc(z)$, $\omega$, $k_x$, and $k_y$ are conserved. If $\psi$ is the angle between the ray and the $z$-axis, then $k\sin \psi =(k_x^2+k_y^2{)}^{1/2}$ and $k=\omega /c$. Hence Snell's law is

For the ray launched with $\mathbf{k}=k(\cos \varphi ,0,\sin \varphi )$ at $z=0$, this constant is $|\cos \varphi |/c_0$. Thereforewhich proves

For $0<\varphi <\pi /2$, write $a=\cos \varphi$. Since the ray and wavevector are parallel,Near the upper turning point, put $z=z_m-\zeta$. Expansion givesIf $x=x_m$ at the turn, integration yields the local ray pathThe ray is tangent to the line $z=z_m$, turns smoothly, crosses $z=0$, and repeats symmetrically between $-z_m$ and $z_m$. Increasing $\varphi$ decreases $z_m$, while $\varphi \downarrow 0$ sends the turning level outward.

Taking the gradient of the quadratic potential energy found in part (ii) resolves the linearized force into identical $x$ and $y$ components:Each component therefore performs simple harmonic motion withThus the angular frequency is

The electric-dipole approximation for radiation requires the source size to be much smaller than the radiation wavelength:The displayed $1/r$ fields also describe the radiation zone, so the observation point must satisfy $r\gg a$ and $r\gg c/\omega$.

Write $t_r=t-r/c$. The given vector potential and the time derivative of the electric dipole moment giveWhen taking the curl, differentiating $1/r$ produces a near-field term of order $r^{-2}$, while differentiating the retarded time produces the leading radiation term because $\nabla t_r=-\hat{\mathbf{x}}/c$. Hence

The assumed electric field $\mathbf{E}=-c\hat{\mathbf{x}}\times \mathbf{B}$ and the transversality $\hat{\mathbf{x}}\cdot \mathbf{B}=0$ make the Poynting vectorOn a sphere of radius $R$, the power through surface element $R^{2}d\Omega$ is thereforeChoose the polar axis along $\ddot{\mathbf{p}}$ and use the angular distribution of electric dipole radiation, for which $|\hat{\mathbf{x}}\times \ddot{\mathbf{p}}{|}^2=|\ddot{\mathbf{p}}{|}^2{\sin }^2\theta$. Since $\int {\sin }^2\theta d\Omega =8\pi /3$,

For the charged particle in simple harmonic motion,ThusHere the angular frequency obeys $\omega =2\pi /T$, and the average uses $⟨{\sin }^2⟩=1/2$. The dipole condition is

Allow the amplitude to vary slowly as radiation removes energy. The mechanical energy isso the radiation damping of a harmonically oscillating charge givesThis is exponential decay with rate $\gamma ={\mu }_0{q}^2{\omega }^2/(6\pi mc)$. Its half-life isFinally, the largest particle speed is $v_{\max }=A\omega$. Therefore the dipole condition $A\omega /c\ll 1$ itself implies $v_{\max }\ll c$, precisely the nonrelativistic limit required by the mechanical energy formula.

The physical rapid phase isDefine the local wavevector, angular frequency, and group velocity byThe equality of mixed partial derivatives gives the phase-compatibility equationsA ray is defined to move with the local group velocity:

Differentiate the local dispersion relationwith respect to $x_i$. Along a ray,The last two terms cancel because $\mathbf{k}$ is a gradient. Similarly,Finally,Thus the Hamiltonian ray-tracing equations are

A stationary state has the formFor a time-independent quantum observable $Q$, the global phases cancel in the expectation value:It is therefore independent of time.

For the infinite square well,The possible energy measurement results arewith respective Born rule probabilities $|c_1{|}^2$ and $|c_2{|}^2$. Normalization gives $|c_1{|}^2+|c_2{|}^2=1$.

At time $t$,The required matrix elements areConsequentlywhere the angular frequency isFinally, $2|c_1{c}_2|\le |c_1{|}^2+|c_2{|}^2=1$, so

At time $t$ the state isDirect evaluation of the radial integral gives the matrix elementsandIt follows thatThus $R(t)$ oscillates sinusoidally about $39a/8$. Its angular frequency is

In the Toomre stability criterionthe disk's gravitational potential supplies the destabilizing attraction, gas pressure represented by $c_s$ suppresses short wavelengths, and the Keplerian shear represented by the angular frequency $\Omega$ prevents coherent collapse on long scales through epicyclic motion. Axisymmetric perturbations are stable for $Q>1$ and gravitationally unstable for $Q<1$.

The state $Q\sim 1$ is a self-regulating fixed point. If cooling lowers $c_s$ until $Q<1$, gravitational instability produces shocks and turbulence that heat the disk and raise its effective velocity dispersion. If $Q$ becomes appreciably larger than one, the instability and its heating switch off, while turbulent dissipation and radiative cooling reduce $c_s$. A sustained gravito-turbulent state therefore remains close to marginal stability.

Equating the two linearized expressions for $p(a,t)-p_{\infty }$ from parts (iii) and (iv) givesThis is the equation of simple harmonic motion. Thus the bubble performs small oscillations with angular frequencyThis is the small oscillation of a planar gas bubble in an annular liquid.

For this metric, $\sqrt{-g}=a^2$, while $g^{\eta \eta }=-a^{-2}$, $g^{zz}=a^{-2}$ and $g^{ij}={\delta }^{ij}$. The massless Klein-Gordon equation is consequentlyInsert the spatial Fourier transform given in the question. Each wavenumber mode then satisfiesThus every field mode is a simple harmonic oscillator with a time-dependent angular frequency.