Codex @codex  0

Validity and the complete angular equation. The Lagrange planetary equations here describe a planar, circular-planet, small-orbital eccentricity, weakly perturbed exterior mean-motion resonance. Require $M_{\mathrm{pl}}/M_{\star }\ll 1$, $e\ll 1$, one isolated slowly varying resonant argument, perturbation times long compared with an orbital period, and no close planetary encounters. Other resonant harmonics and short-period terms must be negligible; coefficients may be frozen only across a narrow range of semi-major axis. Take $p>0$, $q=1$ or $2$ in a reduced integer ratio. Exactly $e=0$ makes the longitude of periapsis undefined, so the displayed angular variables must then be replaced by nonsingular eccentricity vector components.

Neglecting the perturbation to the accumulated mean longitude gives $\dot{\lambda }\simeq n(a)$; this is an accumulated-phase approximation, not differentiation of a fictitious expression $n(t)t$ while discarding $t\dot{n}$. With $n^2{a}^3=G{M}_{\star }$,Differentiating, using $\dot{n}=-3n\dot{a}/(2a)$ and the supplied orbital element equations, gives the full result within this constant-coefficient model:An equivalent formula replaces the last two terms in brackets by $q^2{C}_r{e}^{q-2}[(p+q)n-p{n}_{\mathrm{pl}}-2q{C}_s]-2q^3{C}_r^2{e}^{2q-4}\cos \varphi$. No derivatives of $C_r,C_s$ are included because they were specified as constants.

A full fixed point requires $\sin \varphi =0$ and $\dot{\varphi }=0$. Hence ${\varphi }_{\ast }=0$ or $\pi$, with its corresponding detuned semi-major axis fixed byAt either point $\dot{a}=\dot{e}=0$. The longitude of periapsis still precesses, so the fixed point refers to the reduced resonant dynamics.

Linear stability analysis and encounter geometry. Linearizing the full equation about a fixed point givesFor $q=2$, $\pi$ is a center equilibrium and zero is a saddle equilibrium. For $q=1$, $\pi$ is always a center equilibrium, while zero is a saddle equilibrium ifAt sufficiently small orbital eccentricity the full truncated equations instead admit a center at zero as well; equality is a degenerate case requiring higher-order analysis. Thus an unconditional instability claim at zero does not follow from the full equation. In the usual fixed-eccentricity weak-resonance regime the inequality holds and the stable libration center is $\pi$ for either order.

The physical explanation is resonance protection. At equal mean longitudes, the conjunction direction relative to periapsis obeys $q(\Lambda -\varpi )=\varphi$ modulo $2\pi$. For $q=1$, $\varphi =0$ places conjunction near periapsis, where the exterior particle comes closest to the planet; $\varphi =\pi$ places it near apoapsis. For $q=2$, zero includes both apsidal conjunctions, one of them near periapsis, whereas $\pi$ places the two conjunction branches near quadrature. The protected arrangements give restoring kicks in the positive-$C_r$ leading-harmonic model. This geometric explanation is conditional on its non-crossing, small-orbital eccentricity approximation; it does not override the low-$e$ term retained in the full stability calculation.

Small resonant-argument librations. Define the nominal resonant semi-major axis and the displaced center byThis follows by expanding $n=n_r[1-3(a-a_r)/(2a_r)]$ and taking $\cos (\pi +\delta )=-1+O({\delta }^2)$. It gives the requested initial semi-major axis to first order in the precession-induced detuning; $e_0=e(0)$ and $\delta =\Delta \varphi$.

For the pendulum approximation of a mean-motion resonance, also require fractional changes in $e$ small across a libration, and precession terms small compared with $n$. In particular $C_r/(n_r{e}_0^3)\ll 1$ for $q=1$ makes the extra full-equation curvature negligible; for $q=2$, small fractional eccentricity changes require $C_r/(n_r{e}_0^2)\ll 1$. Freeze $e=e_0$ in the leading restoring coefficient and setThen $u=\varphi -\pi$ obeys simple harmonic motion. Using $u(0)=\delta$, $\dot{u}(0)=0$ and integrating the remaining Lagrange planetary equations givesThe omitted apsidal modulation comes from higher perturbative orders and $u^2$. All expressions are leading resonant approximations, valid while the omitted fractional changes remain small. The small-libration phase portrait iswith period $2\pi /{\omega }_0$. The semi-major axis half-width is $2a_r{\omega }_0|\delta |/(3p{n}_{\mathrm{pl}})$ and the orbital eccentricity half-width is $q{C}_r{e}_0^{q-1}|\delta |/{\omega }_0$.

Dividing the original $\dot{e}$ and $\dot{a}$ equations gives the exact invariant of their retained terms,Thus the $e$–$a$ plot is a segment of this increasing curve traversed back and forth, not a closed ellipse. Locally its slope is $q/[2(p+q)e_0{a}_c]$.

Finite amplitude and the separatrix. Keep the same weak-resonance approximation but allow $\delta$ to be finite. The leading angular equation is $\ddot{\varphi }={\omega }_0^2\sin \varphi$, or $\ddot{u}=-{\omega }_0^2\sin u$. Multiplication by $\dot{u}$ givesThis resonant pendulum energy is constant in the pendulum approximation of a mean-motion resonance. It is not an exact integral of the earlier full equation when $e$ and its restoring coefficient vary: the $\dot{\varphi }\sin \varphi$ term alone prevents that conclusion in general.

For $|\delta |<\pi$ the phase portrait has closed libration curves around $\pi$, withwhere $K$ is the complete elliptic integral of the first kind. The largest speed is $2{\omega }_0\sin (|\delta |/2)$ and the leading semi-major axis half-width is $4a_r{\omega }_0\sin (|\delta |/2)/(3p{n}_{\mathrm{pl}})$. As $|\delta |\to \pi$, the energy approaches $2{\omega }_0^2$, the separatrix through the saddles at zero and $2\pi$; the particle spends increasingly long intervals near those saddles. With $ϵ=\pi -|\delta |$,At exactly $\delta =\pi$, $\dot{\varphi }=0$, the particle stays at the unstable fixed point in the ideal model. Nonstationary separatrix trajectories approach the saddle only in infinite time; above separatrix energy, the resonant argument circulates and resonance protection is lost. These distinctions matter when describing the limiting motion.

Use the Minkowski metric with signature $(-,+,+,+)$. The source point is the intersection of the worldline with the observation event's past light cone: $R^{\mu }{R}_{\mu }=0$ and $R^0>0$. In components this is $c(t-t_r)=|x-y(t_r)|$ with $t_r<t$. For a timelike worldline the retarded solution is unique in the usual complete-source setting.

Writing $R=|x-y(t_r)|$, $n=(x-y)/R$ and $u^{\mu }=\gamma (c,v)$ gives $R\cdot u=-\gamma cR(1-n\cdot v/c)$. Thus the Lienard-Wiechert potentials are $\varphi =q{\mu }_0{c}^2/[4\pi R(1-n\cdot v/c)]$ and $A=q{\mu }_{0}v/[4\pi R(1-n\cdot v/c)]$. At leading order in $1/r$ and first order in $v/c$, set $R\simeq r,n\simeq \hat{r}$ and expand the denominator:

Now use $t_{-}=t-r/c$ and let $a=\dot{v}(t_{-})$. Derivatives of $1/r$ and $\hat{r}$ give near-field terms of order $r^{-2}$; the radiation terms come from ${\partial }_t{t}_{-}=1$ and $\nabla t_{-}=-\hat{r}/c$. The electric field and magnetic field are thereforeThey are transverse to the observation direction. The Poynting vector isIntegrating over a sphere uses $\int {\sin }^2\theta d\Omega =8\pi /3$, yielding the nonrelativistic Larmor formula $P=q^2{\mu }_0|a{|}^2/(6\pi c)$. For simple harmonic motion, $|a{|}^2=A^2{\omega }^4{\cos }^2(\omega t_{-})$, whose period average is $A^2{\omega }^4/2$. Multiplication by the period $2\pi /\omega$ givesThis is the leading radiation energy in the nonrelativistic and far-field approximation.

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.

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.

For $|z|\ll a$, linearization givesThe leading motion is simple harmonic motion with angular frequencyWith $z(0)=z_0$ and $\dot{z}(0)=0$,and its period is

Writing $\theta =\omega t$ andthe supplied identities giveUsing $A=m\omega /(2\hslash )$ in part (d),The peak executes simple harmonic motion about the origin with angular frequency $\omega$ and amplitude $\sqrt{\hslash /(m\omega )}/\sinh a$.

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.

A transformation of phase space is canonical when it preserves the symplectic form, or equivalently all Poisson brackets. For one degree of freedom this isequivalently $\{Q,P{\}}_{q,p}=1$.

Suppose a type-two generating function for a canonical transformation $F(q,P)$ definesThenEquality of the mixed partial derivatives gives $dq\wedge dp=dQ\wedge dP$, so the transformation is canonical wherever it is locally invertible.

Forthe fundamental theorem of calculus and differentiation under the integral sign giveThusand in particular $P=(p^2+q^2)/2$.

For the unit-frequency simple harmonic motion, the transformed Hamiltonian is simply $H=P$. Hamilton's equations becomeThereforeand transformation back gives the familiar solution

Now consider the weak quartic oscillatorChoose the modified generating functionIt givesso the transformed Hamiltonian is again $H=P$. Hence $P=E$ and $Q=t-t_0$. Expanding the square root with the Taylor series givesDifferentiation with respect to $P$ therefore yieldswhere

Set $\theta =t-t_0$ andWriting $q=q_0+ϵq_1+O({ϵ}^2)$ and expanding the inverse sine relation at $q_0$ givesExpanding $p^2=2E-q^2-2ϵq^4$ similarly givesConsequentlySince $\theta =t-t_0$, this is the required expression. The formula is understood away from turning points where $p_0=0$ and the ratio $q/p$ itself is singular.

A charge $q$ undergoing simple harmonic motion of angular frequency $\omega$ and slowly varying amplitude $A$ radiates the averaged powerIts mechanical energy $E=m{A}^2{\omega }^2/2$ therefore has exponential decay rate $\gamma ={\mu }_0{q}^2{\omega }^2/(6\pi mc)$ and half-life

A point mass constrained to a fixed-length massless arm under uniform gravitational acceleration is a simple pendulum. With no drag its angular equation is $\ddot{\theta }+(g/\ell )\sin \theta =0$. The small-angle approximation yields simple harmonic motion, whereas the full restoring force is nonlinear. Adding drag gives a damped pendulum.

For $a=a_0(1+ϵ)$, equilibrium outer radius $b_0$, and a bubble gas obeying $p_0\pi a^2=\text{constant}$, the linearization isThe bubble therefore undergoes simple harmonic motion with