Codex @codex  0

Choose the Fourier series conventionFor a real potential energy, $V_{-1}=V_1^{\ast }$. The symbol $V(G_1)$ means the Fourier coefficient $V_1$, with units of energy, rather than the pointwise value of $V(x)$ at a position $G_1$.

Retain the two plane waves $|1⟩=L^{-1/2}{e}^{ikx}$ and $|2⟩=L^{-1/2}{e}^{i(k-G_1)x}$ in the nearly-free electron model. Their kinetic energies are ${ϵ}_1={\hslash }^2{k}^2/(2m)$ and ${ϵ}_2={\hslash }^2(k-G_1{)}^2/(2m)$. The matrix representation of a linear map for the Quantum Hamiltonian projected onto this two-dimensional basis isThus $\det (H-EI)=0$ givesSolving this quadratic equation yieldsChoosing the zero of potential energy so that its spatial average is $V_0=0$ reproduces the requested expression. If that average has not been removed, both energy bands are shifted by $V_0$; no band gap, group velocity or effective mass is changed by this constant shift.

The two kinetic energies coincide at $k=G_1/2$. Their coupling gives an avoided crossing with band gap $2|V_1|$. The result is exact within the projected two-wave Quantum Hamiltonian. It is an approximation to the full crystal problem, appropriate near this crossing for a weak periodic potential whose other Fourier coefficients do not cause strong mixing with remote plane waves. Small $|V_1|$ alone would not justify ignoring arbitrarily large other Fourier coefficients.

Let the infinite outer walls be at $x=0,D$, the finite steps at $x=a,b$, and the three potential energy values be $V_L,V_B,V_R$. The original drawing specifies the orderingThe energy eigenstates can be chosen real. The Dirichlet boundary conditions are $\psi (0)=\psi (D)=0$; wavefunction matching at a finite potential step requires both $\psi$ and ${\psi }^{\prime }$ to be continuous at $a,b$. In particular, neither finite step acts as an additional hard wall.

The time-independent Schrödinger equation givesThus a wavefunction is sinusoidal in a classically allowed region, with wavenumber $k=\sqrt{2m(E-V)}/\hslash$, and a linear combination of exponential functions in a classically forbidden region, with inverse decay length $q=\sqrt{2m(V-E)}/\hslash$. Its curvature changes sign relative to $\psi$ at a step crossing $E$; its slope does not jump.

By the Sturm oscillation theorem, the three energy eigenstates have exactly zero, one, and two interior nodes, respectively. Choose their global phases so that their slopes at the left wall are positive:

For $E_0,E_1$, the right-hand piece is proportional to $\sinh [q_R(D-x)]$, rather than a pure decaying exponential: this enforces the finite-distance hard wall exactly. The barrier pieces likewise contain both exponential solutions. Their small tails are an instance of quantum tunnelling.

This is an original illustrative calculation, with vertical lines marking only the region boundaries. It uses $D=1$, $a=0.4$, $b=0.6$, $V_L=0$, $V_R=0.62$, $V_B=0.82$, and ${\hslash }^2/(2m{D}^2)=0.002$ in arbitrary energy units. Its lowest energy eigenvalues are approximately $0.098,0.384,0.696$, so it has exactly the required ordering. Here both nodes of the third wavefunction happen to lie towards the left; its dominant lobe is on the right. These precise positions and relative amplitudes are illustrative: the original drawing supplies no numerical lengths, barrier heights, or mass, so it does not determine a unique quantitative wavefunction. The required sketch features are the matching, curvature, hard-wall zeros, and ordered node counts.

For a particle in a ring, the arc-length coordinate is $s=R\varphi$. The tangential kinetic energy isThe constraint removes radial motion; the constant potential energy along the ring can be chosen as zero. The tangential momentum operator is ${\hat{p}}_s=-i\hslash {\partial }_s$, or equivalently the angular momentum operator is ${\hat{L}}_z=-i\hslash {\partial }_{\varphi }$. ThusIts physical domain has periodic boundary conditions, $\psi (\varphi +2\pi )=\psi (\varphi )$ and ${\psi }^{\prime }(\varphi +2\pi )={\psi }^{\prime }(\varphi )$. These make the Hamiltonian self-adjoint and reflect that $\varphi$ and $\varphi +2\pi$ describe the same point. No magnetic flux or additional angle-dependent potential energy is present.

Time dependence in quantum mechanics describes how a prepared quantum state changes and thereby changes the probabilities of later measurements. For an isolated system, the Schrödinger picture assigns this evolution to a state vector $|\psi (t)⟩$, governed by the Time-dependent Schrödinger equationThe Hamiltonian is a self-adjoint operator; it has units of energy, while $\hslash$ has units of action, so $\hat{H}/\hslash$ has units of inverse time. For a particle with potential energy $V(\mathbf{r},t)$, this becomesGiven a suitable initial wavefunction and the physical boundary conditions, this first-order equation determines its subsequent propagation. Its linearity means that the evolution of a quantum superposition is the corresponding quantum superposition of evolved states.

For a real potential energy, this equation and its complex conjugate imply the local continuity equationHere $\rho$ is the probability density and $\mathbf{j}$ the probability current. Vanishing net boundary flux preserves wavefunction normalization, as required by the Born rule.

The time-evolution operator packages the solution for every initial state:Because ${\hat{H}}^{†}=\hat{H}$, differentiation givesTherefore $U^{†}U=I$: this is unitary time evolution, preserving inner products and all state norms. It is reversible, with $U(t,t_0{)}^{-1}=U(t_0,t)=U(t,t_0{)}^{†}$, and obeys the composition law $U(t_2,t_0)=U(t_2,t_1)U(t_1,t_0)$.

For a time-independent Hamiltonian, the solution is the operator exponentialIf the Hamiltonian varies with time, the analogous simple exponential of its integral is valid when $[\hat{H}(t),\hat{H}(t^{\prime })]=0$ for all relevant times. Otherwise time ordering is essential:The meaning is explicit in the first terms of the ordered series,the later-time operator stands on the left. In the position representation, the same propagation is expressed by the quantum-mechanical propagatorIts initial value is the Dirac delta function ${\delta }^3(\mathbf{r}-{\mathbf{r}}^{\prime })$.

A time-independent Hamiltonian has energy eigenstates satisfying $\hat{H}|n⟩=E_n|n⟩$. A single such initial state evolves asIt is a stationary state: its global phase changes, but its probability density and the distributions of time-independent observables do not. Its energy is sharp. Any quantum superposition entirely within one degenerate energy eigenspace also acquires one common global phase and is stationary. Stationarity does not require zero probability current; for instance, $u_n(\varphi )=e^{in\varphi }/\sqrt{2\pi }$ on the ring has a constant density and persistent circulation when $n\ne 0$.

An arbitrary normalized state can instead contain different energy eigenvalues:A continuous spectrum requires corresponding integrals. For the discrete case, an energy measurement gives $E_n$ with fixed probability $|c_n{|}^2$ (summed over degeneracies), but relative quantum phases change with time. ConsequentlyThe interference terms can oscillate at angular frequencies $(E_m-E_n)/\hslash$. For two populated energy levels, this supplies a characteristic period $2\pi \hslash /|E_m-E_n|$ when the interference term is nonzero. Thus fixed energy populations do not imply a stationary spatial distribution.

For a normalized state and an observable $\hat{A}(t)$, its expectation value isThis is the statistical mean over repeated preparations and measurements. In an eigenbasis of energy eigenstates for a time-independent Hamiltonian and a time-independent $\hat{A}$,Only matrix elements connecting different energies can generate these oscillations. If $[\hat{A},\hat{H}]=0$, its expectation value is constant for every initial state, provided it has no explicit time dependence.

Differentiating the expectation value and substituting the Schrödinger equation for both the bra and ket proves the general Ehrenfest theorem:Equivalently,These manipulations require the state to lie in the relevant operator domains. Taking $\hat{A}=\hat{H}$ gives $d⟨H⟩/dt=⟨{\partial }_t\hat{H}⟩$, so energy is conserved under a time-independent Hamiltonian.

For $\hat{H}={\hat{p}}^2/(2m)+V(\hat{x})$, the canonical commutation relation givesHence the particularly useful position-and-momentum form of Ehrenfest theorem isIt resembles Newton's second law, but the mean force generally differs from the force at the mean position. If a narrow wave packet samples a region where the force varies little, $⟨V^{\prime }(x)⟩\simeq V^{\prime }(⟨x⟩)$ gives approximate classical motion. For a quadratic potential energy this replacement is exact. In particular the quantum harmonic oscillator obeysThe centroid follows the classical harmonic oscillator, although this does not specify the full evolving wavefunction or its width.

The Heisenberg picture expresses the same physics by assigning evolution to observables instead:The expectation value is unchanged, and differentiation yields $d{\hat{A}}_H/dt=(i/\hslash )[{\hat{H}}_H,{\hat{A}}_H]+U^{†}({\partial }_t{\hat{A}}_S)U$. Thus state propagation, evolving expectation values, and the classical correspondence in Ehrenfest theorem are mutually consistent descriptions of the same unitary time evolution.

An Euler strut is a slender initially straight elastic column under axial compression. For a uniform bending stiffness $EI$ and a small transverse displacement $y(x)$, the second-order potential energy under a conservative compressive load $P$ isThe bending term is positive; shortening releases load potential. Stationarity yields $EI{y}^{(4)}+P{y}^{\prime \prime }=0$. Pinned ends impose $y(0)=y(L)=0$ and zero end bending moments, $y^{\prime \prime }(0)=y^{\prime \prime }(L)=0$. Substitution of $y=A\sin (n\pi x/L)$ givesThe straight configuration loses stability first in the longest admissible bending mode. The second moment of area $I=\int z^{2}dA$ is geometric, not a mass moment of inertia. Boundary conditions matter: a fixed-free column has effective length $2L$, while fixed-fixed ends have effective length $L/2$, changing the threshold by factors $1/4$ and $4$. Slenderness, elastic material response, alignment and the weak bending direction must be checked. Imperfections produce bending before the ideal threshold, and a short column can yield or crush first. Linear buckling theory does not predict an unlimited deflection or a complete postbuckling curve.

A normal mode is a motion in which all displacement coordinates oscillate at one angular frequency with fixed amplitude ratios. For small conservative oscillations, the symmetric stiffness and mass matrices give the generalized eigenvalue problem for small oscillations, $\mathsf{K}\mathbf{u}={\omega }^2\mathsf{M}\mathbf{u}$. For positive mass and a stable equilibrium, its mass-orthogonal eigenvectors span the motion. A superposition of these modes, with phases and amplitudes fixed by initial displacements and velocities, gives the general solution; different modes need not share a phase.

Use upward end displacements. The spring compressions are $x_0-x_1$ and $x_0-x_2$, while the beam's centre rises by $q=(x_1+x_2)/2$. The potential energy isStatic force balance gives $2k{x}_0=mg$, cancelling the linear terms and leavingFor a uniform beam, its central moment of inertia in the vertical plane is $m{l}^2/12$. To first order in the tilt, $\theta =(x_2-x_1)/l$. Its translational and rotational kinetic energies areThus $\mathsf{M}=(m/6)\left(\begin{array}{cc}2& 1\\ 1& 2\end{array}\right)$ and $\mathsf{K}=k1$. The normal modes of a rigid beam on two end springs have eigenvectorsThe first translates without tilt; the second rocks about a stationary centre. Ordinary frequencies are ${\omega }_s/(2\pi )$ and ${\omega }_a/(2\pi )$.

Assume the system had settled before the end load leaves, and put that load at end 1. Its extra static spring load is $Mg$ at that end and zero at the other. Therefore $x_1(0)=-Mg/k$, $x_2(0)=0$, and both initial velocities vanish because no impulse is imparted. Decomposing into the two normal modes givesOnly the symmetric mode contributes to the centre position:The datum $z_{\mathrm{eq}}$ is the unloaded beam's equilibrium height; no absolute height is supplied. The rocking mode is present even though it cancels out of this coordinate. These results use small tilt, $Mg/(kl)\ll 1$, ideal light vertical springs and negligible damping.

For the scattering calculation, take the heavy target as fixed and define $\kappa =2Z{e}^2/(4\pi {ϵ}_0)>0$. The incoming kinetic energy and angular momentum are $E=m{v}^2/2$ and $J=mvb$. In the original sketch, the horizontal dashed incoming asymptote has perpendicular distance $b$ from the nucleus, while $d$ is the radius from the nucleus to the point where $\dot{r}=0$.

At closest approach, the velocity is tangential and $v_t=J/(md)=vb/d$. conservation of energy givesChoose the positive root and put $s=\kappa /(m{v}^2)$:Head-on incidence gives $d_{\mathrm{rep}}=2s$, and repulsion always gives $d_{\mathrm{rep}}>b$. With an attractive Coulomb force, change the sign of the potential term:For $b>0$, the attractive and repulsive inverse-square closest approaches obey $d_{\mathrm{att}}{d}_{\mathrm{rep}}=b^2$. For $b=0$ attraction causes a collision with the point centre, so $d=0$ is a singular endpoint, not a regular scattering turning point. Finite nuclear size is excluded by the stated point-force model.

For the repulsive inverse-cube central force $\mathbf{F}={\kappa }_3\hat{\mathbf{r}}/r^3$, the torque vanishes, $\mathbf{r}\times \mathbf{F}=0$, and angular momentum is still conserved. Its potential energy is $+{\kappa }_3/(2r^2)$, since $-dV/dr=+{\kappa }_3/r^3$. ThusThe new coefficient is denoted ${\kappa }_3$ because its units are $\mathrm{N}{\mathrm{m}}^3$, unlike the $\mathrm{N}{\mathrm{m}}^2$ units of the inverse-square coefficient. Conserving angular momentum depends on centrality, not on a particular radial power law.

A velocity field $\mathbf{u}(\mathbf{x},t)$ assigns a velocity to each point. Following a fluid parcel rather than holding position fixed gives the material or convective derivativeThe convective term can accelerate fluid even in a steady flow, for example as a stream narrows. Conservation of mass is ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$; constant-density incompressible flow has $\nabla \cdot \mathbf{u}=0$.

For an incompressible Newtonian fluid of constant density and viscosity, with conservative body-force potential per mass $\Phi$, momentum balance isNeglect viscosity and assume steady flow. Dotting with an element along a streamline, or using $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (u^2/2)-\mathbf{u}\times (\nabla \times \mathbf{u})$, givesThis is Bernoulli's equation: pressure work, kinetic energy and gravitational potential energy interconvert. In a connected irrotational region under the same assumptions, $B$ is the same for all streamlines. In a rotational flow it can vary between them. Viscous dissipation, pumps or unsteady motion require the corresponding modified balance; the steady inviscid expression cannot be applied indiscriminately.

Vorticity $\mathbf{\zeta }=\nabla \times \mathbf{u}$ measures local fluid rotation. For constant-density incompressible flow it evolves byshowing vortex stretching and viscous diffusion. A potential flow has $\mathbf{u}=\nabla \varphi$ locally, so irrotational incompressible motion satisfies ${\nabla }^2\varphi =0$. Boundary conditions then determine flows around simple shapes. A multiply connected region can have circulation even when its local vorticity vanishes; ideal potential flow also misses no-slip boundary layers and viscous drag.

For a Newtonian fluid, shear stress is proportional to strain rate. No-slip at solid walls produces velocity gradients and dissipation. Fully developed laminar pipe flow, for example, hasThe Reynolds number $\mathrm{Re}=UL/\nu$ compares inertial and viscous terms. Small values favour creeping-flow approximations; large values can permit thin boundary layers and turbulence, without guaranteeing turbulence in every geometry. The appropriate flow model follows from conservation laws, constitutive properties, boundaries and scale comparisons.

For a moving loop, Faraday's law of induction is $\mathcal{E}=\oint (\mathbf{E}+\mathbf{u}\times \mathbf{B})\cdot d\mathbf{l}=-d{\Phi }_B/dt$, where $\mathbf{u}$ is the local wire velocity, ${\Phi }_B=\int \mathbf{B}\cdot d\mathbf{S}$, and the oriented boundary and surface follow the right-hand rule. For a fixed loop the motional term vanishes. Lenz's law is the minus sign: the induced response opposes the change of linked flux, not necessarily the applied field itself.

Let $s$ increase downhill and choose positive current $I$ to produce an uphill magnetic force. The normal component of the imposed field is $B\cos \theta$, so the changing external flux is $Bl\cos \theta s$. With initially zero current and no resistance, conservation of flux linkage gives $LI=Bl\cos \theta s$. The rod feels gravity vertically downward, a rail reaction normal to the slope, and magnetic force of magnitude $IlB$ horizontally uphill. Its downhill component is $-IlB\cos \theta$; its normal component presses the rod into the rails, giving $N=mg\cos \theta +IlB\sin \theta$.

The inductive oscillation of a sliding rod on inclined rails follows fromThe magnetic force is not itself parallel to the rails. Defining ${\omega }_B=Bl\cos \theta /\sqrt{mL}$ and $s_e=mg\sin \theta /(m{\omega }_B^2)$, release from rest givesThe first stop occurs at $t=\pi /{\omega }_B$; the rod then returns and oscillates between $0$ and $2s_e$. At $s_e$ the force is zero but the speed is maximal. Mechanical potential energy is reversibly exchanged with kinetic energy and $L{I}^2/2$; a perfectly conducting circuit provides no damping.

For the resistive long-time drift of an inductive sliding rod, with loop resistance $R>0$, useJoule heating damps the transient oscillation. On sufficiently long rails the asymptotic motion is steady sliding, not permanent rest:The zero-resistance model has persistent oscillations; its late-time limit cannot be interchanged with taking $R\to 0$. These results assume $B\cos \theta \ne 0$, smooth contacts, negligible rod self-inductance and constant geometry.

The shell theorem gives enclosed halo mass $4\pi \rho r^3/3$. Define $\alpha =4\pi \rho /3$. Then $F_r=-GMm/r^2-\alpha Gmr$. Interior to the uniform halo its potential energy may be writtenOuter shells alter the additive constant, not this interior force. An actually infinite uniform mass distribution has no finite potential referenced to infinity; only local potential differences are needed here.

The central force has zero torque, so $\mathbf{J}=\mathbf{r}\times m\dot{\mathbf{r}}$ is conserved and the motion is planar. With $J=m{r}^2\dot{\varphi }$, the effective radial potential and equation areFor circular radius $r_0$, $J=m{r}_0^2\Omega$ givesSet $r=r_0+\xi$ at fixed $J$. The derivative of the right-hand side at equilibrium is $-3{\Omega }^2+2GM/r_0^3-\alpha G$, soBoth frequencies are positive, but radial and azimuthal motion no longer close with the Kepler ratio. At first order in eccentricity $e$, one may write $r=r_0(1-e\cos \omega t)$ and $\varphi =\Omega t+(2\Omega e/\omega )\sin \omega t$. For a weak halo this is a nearly Keplerian ellipse whose apsis rotates at angular rate $\Omega -\omega$.

Put ${\Omega }_K=\sqrt{GM/r_0^3}$ and $\epsilon =\alpha r_0^3/M\ll 1$. Expansion yieldsThe precession is retrograde. In cycles per unit time the precession frequency is $\dot{\varpi }/(2\pi )$; per radial period the apsidal shift is $2\pi (\Omega /\omega -1)\simeq -3\pi \epsilon$. This weak-halo interpretation avoids confusing a focus-centred Kepler ellipse with the centred ellipse of a purely harmonic potential.

In electrostatics, the Maxwell equations in matter areHere ${\rho }_f$ is free charge density and $\mathbf{P}$ is polarization density. For a homogeneous isotropic linear dielectric, the constitutive equation is $\mathbf{D}={ϵ}_0{ϵ}_r\mathbf{E}$. Integrating the first equation over a thin pillbox straddling an interface, with unit normal $\mathbf{n}$ from medium 1 to medium 2, gives $\mathbf{n}\cdot ({\mathbf{D}}_2-{\mathbf{D}}_1)={\sigma }_f$. Integrating the second around a vanishingly narrow rectangular loop gives $\mathbf{n}\times ({\mathbf{E}}_2-{\mathbf{E}}_1)=0$. These are the electrostatic interface boundary conditions. When there is no free surface electric charge, the normal electric displacement field and tangential electric field are continuous; the normal electric field generally is not.

Neglect end fringing field of a capacitor. In either axial segment of the coaxial cylindrical capacitor, a cylindrical Gaussian surface gives $D_r={\lambda }_f/(2\pi r)$, where ${\lambda }_f$ is free electric charge per length. Integrating the radial electric field between the conductors givesThe liquid and air segments share the same conductor electric potential difference, so their capacitances add in parallel. Thus the partially filled coaxial capacitor hasHere the question's $ϵ$ is relative permittivity, and air's relative permittivity is approximated by one.

Write $x$ for the rise above the external level. The wide reservoir keeps that level approximately fixed. At constant applied voltage, one must include work done by the source: the appropriate electrical contribution to the mechanical potential is $-C{V}^2/2$. The upward dielectric insertion force at fixed voltage is consequentlyThe excess hydrostatic pressure $\rho gx$ acts over the annular area $A=\pi (b^2-a^2)$. Alternatively its gravitational potential energy is $\rho gA{x}^2/2$, since liquid initially at the reservoir level must be raised through all heights from zero to $x$. Balancing the derivatives givesThis assumes $h+x<L$, a communicating liquid reservoir, negligible capillarity, and a sufficiently long capacitor that end effects are small. Differentiating gives the dielectric-rise voltmeter's sensitivity at a bias $V$:To obtain its numerical value in $\mu \mathrm{m}{\mathrm{V}}^{-1}$, multiply the right-hand side in SI units by $10^6$. A convenient parametric scale, with $\beta =b/a$, isFor a concrete estimate, choose $a=1\mathrm{mm}$, $b=2\mathrm{mm}$, mass density $\rho =10^3\mathrm{kg}{\mathrm{m}}^{-3}$, relative permittivity $ϵ=2$ and operating voltage $V=100\mathrm{V}$. These are illustrative choices, not supplied apparatus data. They giveThus this millimetre-scale example requires fine optical displacement measurement for volt-level changes. The parametric formula states exactly how the estimate changes with the chosen geometry, liquid and bias.

The displacement is quadratic in voltage, cannot distinguish its sign, and has zero first-order sensitivity at zero bias. Smaller gaps and larger relative permittivity increase sensitivity, but capillarity, leakage, dielectric breakdown, temperature-dependent liquid properties, and finite available rise then limit performance. A nonconducting liquid avoids a sustained conduction load, although charging still draws electric charge and leakage is possible. A high voltage can give a readable displacement; a low-voltage instrument would require fine displacement measurement and careful calibration.

For the attractive inverse-square force, choose zero potential energy at infinity, so $U=-A/r$. At either radial turning point $r_p,r_a$, the kinetic energy is purely tangential, with conserved angular momentum $J$. ThusIts two roots are the periapsis and apoapsis distances. Their sum is $r_p+r_a=-A/E$, whereas the ellipse has $r_p+r_a=2a$. Therefore $a=-A/(2E)$, with $E<0$ for a nondegenerate bound ellipse. The circular case follows by the double-root limit.

For the two-body problem, define $\mathbf{r}={\mathbf{r}}_B-{\mathbf{r}}_A$, total mass $M=m_A+m_B$, reduced mass $\mu =m_A{m}_B/M$, and center of mass $\mathbf{R}$. The relative equation is $\mu \ddot{\mathbf{r}}=\mathbf{F}$, while the isolated center of mass has $\ddot{\mathbf{R}}=0$. Once the relative orbit is known,The relative energy uses $\mu$ in its kinetic energy and $A=G{m}_A{m}_B$ in the gravitational potential energy.

In the equal-mass case the initial relative speed is $u$, $\mu =m/2$, and the center-of-mass velocity is half B's initial velocity. ThusNegative relative energy requires $u^2<4Gm/a_1$. The escape threshold is $u_E=2\sqrt{Gm/a_1}$; equality is marginal parabolic escape, so a strictly bound pair has $u<u_E$.

For $u=u_E/2=\sqrt{Gm/a_1}$, $E_{\mathrm{rel}}=-3G{m}^2/(4a_1)$ and the relative semimajor axis is $a=2a_1/3$. The initial transverse speed is below the relative circular value $\sqrt{2Gm/a_1}$, so the initial turning point is apoapsis. HenceChoose the initial separation along $+x$ from A to B, with B's initial velocity $+u\hat{\mathbf{y}}$. Conservation of angular momentum gives relative periapsis speed $3u$. At the opposite end of the ellipse the direction is $-\hat{\mathbf{y}}$, soThus B has the same speed as initially but the opposite velocity. A has $2u\hat{\mathbf{y}}$ at that instant. A value $3u/2$ for B would be its center-of-mass-frame speed, not the requested initial laboratory frame result.

The relative orbital period follows from Kepler's third law with total gravitating mass $2m$:In the center-of-mass frame each star follows an ellipse of semimajor axis $a_1/3$ and eccentricity $1/2$, with the center of mass at a focus; they remain opposite on their simultaneous position vectors. In the laboratory frame add the uniform drift $ut/2$ along $y$. The laboratory paths of an equal-mass binary are therefore not closed ellipses, even though their relative motion repeats.

In Lagrangian mechanics, choose independent generalized coordinates $q_i$ compatible with the mechanical constraints, write $L=T-U$, and apply the Euler-Lagrange equationEquivalently, the physical path makes the action stationary under variations with fixed endpoints. For ideal holonomic constraints, their unknown reaction forces need not appear in the reduced equations. Direct use of Newton's laws of motion instead resolves vector forces and reactions and enforces the constraints separately. Both approaches describe the same mechanics; the scalar Lagrangian is often shorter in curvilinear coordinates. Nonconservative generalized forces must be added separately. Here a drive maintains the prescribed azimuthal angular velocity $\Omega$.

The mass has polar speed $a\dot{\theta }$ and azimuthal speed $a\Omega \sin \theta$, so its kinetic energy is $m{a}^2({\dot{\theta }}^2+{\Omega }^2{\sin }^2\theta )/2$. The elastic potential energy is $\mu {\theta }^2/2$, and the height is $a\cos \theta$, giving gravitational potential energy $mga\cos \theta$. Thereforeand the Euler-Lagrange equation givesFor $\mu =2mga$, define $q={\Omega }^{2}a/g$ and $F(\theta )=q\sin \theta \cos \theta -2\theta +\sin \theta$, so $\ddot{\theta }=(g/a)F(\theta )$. Near zero,When $q<1$ the upright equilibrium is linearly stable and the small oscillations have ${\omega }_1^2=g/a-{\Omega }^2$.

For $q>1$, the graph of $F$ initially rises, whereas $F(\pi /2)=1-\pi <0$. There must therefore be a zero inside the interval. More strongly, for $0<\theta <\pi /2$,is strictly decreasing: $\cos \theta$ decreases and $\theta /\sin \theta$ increases because $\sin \theta -\theta \cos \theta >0$. The ratio goes from $q-1>0$ to $1-\pi <0$, giving exactly one zero ${\theta }_0$ in that interval. At the zero, $F^{\prime }({\theta }_0)<0$, so it is a stable equilibrium. The rotating torsion-stabilized inverted pendulum also has a symmetric stable angle $-{\theta }_0$ if both signs of tilt are permitted.

At ${\theta }_0=\pi /4$, setting $F({\theta }_0)=0$ gives $q/2-\pi /2+1/\sqrt{2}=0$, soFor $\theta ={\theta }_0+\xi$, the linearized equation is $\ddot{\xi }=(g/a)F^{\prime }({\theta }_0)\xi$. Since $F^{\prime }=q\cos 2\theta -2+\cos \theta$,At the threshold $q=1$ the linear restoring term vanishes, but the cubic term remains:The upright state remains nonlinearly stable, with no nonzero infinitesimal harmonic frequency. Small finite disturbances produce quartic oscillator motion whose frequency is proportional to its small angular amplitude; its period tends to infinity as that amplitude tends to zero. Zero linear stiffness is not the same as an absent restoring torque.