Codex @codex  0

For the elastic pendulum, set $r=l+x>0$ and measure $\theta$ from the downward vertical. In the vertical plane, Cartesian coordinates are $X=r\sin \theta$ and $Z=-r\cos \theta$, with $Z$ upward. Differentiating gives ${\dot{X}}^2+{\dot{Z}}^2={\dot{x}}^2+r^2{\dot{\theta }}^2$. Thus the Lagrangian isFor $x$, the Euler-Lagrange equation givesFor $\theta$, the canonical momentum is $p_{\theta }=m(l+x{)}^2\dot{\theta }$. Its Euler-Lagrange equation isor, for positive length,At an equilibrium, all time derivatives vanish, so $kx=mg\cos \theta$ and $\sin \theta =0$. Modulo full rotations, the two formal possibilities areThe lower configuration always has positive length. The upper one is a physical straight-spring configuration only if $l>mg/k$. For example, $l=mg/(2k)$ would give the purported upper state a negative length $-mg/(2k)$; at $l=mg/k$ its zero length makes the angular coordinate singular. Also, the straight compressed spring is an idealization; a real unsupported compressed spring can buckle.

The potential energy is $V=k{x}^2/2-mg(l+x)\cos \theta$. At either formal stationary point its mixed second derivative vanishes, whileBoth curvatures are positive at the lower state, making it a stable equilibrium; the upper state, when it has positive length, has negative angular curvature and is unstable.

To linearize about the stable state, define the small displacement $\xi =x-mg/k$, not $x$ itself, and $r_0=l+mg/k$. Products such as ${\dot{\theta }}^2$, $\dot{\xi }\dot{\theta }$ and $\xi \ddot{\theta }$ are second order, while $\cos \theta =1+O({\theta }^2)$ and $\sin \theta =\theta +O({\theta }^3)$. The first-order equations of motion therefore decouple:There are two independent normal modes. The radial one is vertical simple harmonic motion with angular frequency ${\omega }_r=\sqrt{k/m}$. The transverse one is a simple pendulum oscillation at fixed, statically stretched length $r_0$, with angular frequency ${\omega }_{\theta }=\sqrt{g/r_0}$. Their decoupling is a first-order small-oscillation result, not an assertion about the nonlinear motion of the elastic pendulum.

First take $a\ne 0$, put $A=|a|$, and suppose $k>0$. The force is the negative derivative of the potential energy:The potential energy is an odd function, passes through $U(0)=0$, and approaches zero from above as $x\to -\infty$ and from below as $x\to +\infty$. Its stationary points are $x=\pm A$:Thus the sole stable equilibrium is the local minimum at $x=A$; the local maximum at $x=-A$ is an unstable equilibrium. The requested values at the stable equilibrium areThe force is an even function, positive for $|x|<A$, negative for $|x|>A$, and zero at $x=\pm A$. Its central maximum is $F(0)=k/A^2$. Sinceits two minima occur at $x=\pm \sqrt{3}A$, with $F=-k/(8A^2)$, and $F\to 0^{-}$ at both ends. These values fix the shapes of both requested sketches.

In one dimension, attraction towards the origin means $xF(x)<0$, while repulsion from the origin means $xF(x)>0$. With that explicitly stated convention,A positive force alone means rightward, not necessarily outward. At the origin $F(0)\ne 0$, so it is not an equilibrium; the inward/outward distinction is undefined exactly there.

The sign assumptions matter because the paper calls $a,k$ merely constants. For $k<0$ both curves reverse sign, the stable equilibrium is $x=-A$, its potential energy is $-|k|/(2A)$, and the attraction/repulsion intervals interchange. For $k=0$ the force vanishes everywhere on the nonsingular domain and there is no isolated potential energy minimum. If $a=0$ and $k\ne 0$, then $U=-k/x$ for $x\ne 0$ and $F=-k/x^2$, with no finite equilibrium.

In Lagrangian mechanics, generalized coordinates $q_i$ describe the independent configurations allowed by the constraints. For particles with ordinary conservative forces in an inertial frame, use $L=T-V$, with kinetic energy $T$ and potential energy $V$. Holonomic constraints can often be built into the coordinates, avoiding explicit constraint forces.

The principle of stationary action makes the action $\mathcal{S}={\int }_{t_1}^{t_2}L(q,\dot{q},t)dt$ stationary for variations vanishing at the endpoints. Integrating the velocity variation by parts givesSince each independent variation is arbitrary, the Euler-Lagrange equations areIf nonconservative generalized forces are present, their components appear on the right. The formalism is invariant under suitable coordinate changes and adding a total time derivative to $L$ leaves the fixed-endpoint equations unchanged.

The generalized momentum is $p_i=\partial L/\partial {\dot{q}}_i$. If $q_i$ is cyclic, $\partial L/\partial q_i=0$ and $p_i$ is conserved. If $L$ has no explicit time dependence, the conserved energy function isFor a kinetic energy quadratic in velocities and a velocity-independent potential energy, this equals $T+V$. This connects symmetry to conservation laws.

For example, a simple pendulum of mass $m$ and length $\ell$ hasIts tension does not appear because the fixed-length constraint is already imposed. Expanding for small $\theta$ gives simple harmonic motion with $\omega =\sqrt{g/\ell }$. More generally, expansion of $T$ and $V$ to quadratic order about a stable equilibrium leads directly to the normal-mode eigenvalue problem.

Expand the cloud's potential energy about a stable equilibrium: $U(z)=U(0)+U^{\prime }(0)z+U^{\prime \prime }(0)z^2/2+\cdots$. Equilibrium gives $U^{\prime }(0)=0$, stability gives $U^{\prime \prime }(0)=k>0$, and reflection symmetry also removes odd terms. For small displacements, $U-U(0)\simeq k{z}^2/2$ and $F=-dU/dz\simeq -kz$.

Including the electric and damping forces, Newton's second law is $m\ddot{z}=qE(t)-b\dot{z}-kz$. Dividing by $m$ gives the Lorentz oscillator modelThe specified convention uses $2\gamma$ as the damping coefficient per unit mass.

For the steady response, represent the field by $E_0{e}^{i\omega t}$ and the displacement by $\tilde{z}{e}^{i\omega t}$, taking real parts for physical quantities. Substitution givesDefine the nonnegative amplitude $z_0=|\tilde{z}|$ and $\delta =atan2(2\gamma \omega ,{\omega }_0^2-{\omega }^2)$. ThenFor positive $q$, $\varphi =-\delta$. For electrons $q<0$, the phase of a negatively charged driven oscillator is $\varphi =\pi -\delta$: the low-frequency displacement is opposite to the field, rather than in phase with it. A signed amplitude could absorb this extra $\pi$, but the phase convention must then be stated. The induced dipole $p=qz$ has phase $-\delta$ for either charge sign. Transient free oscillations have been omitted because the question asks for steady state.

For the electron cloud, $q/m=-e/m_e$ andThe visible range corresponds to $\omega =2\pi c/\lambda$ between $2.35\times 10^{15}$ and $4.71\times 10^{15}\mathrm{rad}{\mathrm{s}}^{-1}$, below this resonance. The resonant wavelength is about $180\mathrm{nm}$, in the ultraviolet.

The acceleration amplitude is ${\omega }^2{z}_0$. At fixed incident field amplitude, the scattered intensity is therefore proportional toWell below resonance and outside strong damping corrections, $z_0$ is nearly frequency-independent and this becomes the Rayleigh scattering wavelength dependence $I_{\mathrm{sc}}\propto {\omega }^4\propto {\lambda }^{-4}$. The shortest visible wavelengths are consequently scattered much more strongly into directions away from the direct solar beam than the longest ones, producing a blue sky. The simple fourth-power approximation gives a factor 16 between 400 and 800 nm at equal incident intensity. Retaining the undamped denominator with the given resonance raises that factor to about 22.7, since the upper end is not extremely far below resonance. The exact ratio needs the unspecified damping and incident solar spectrum; the stronger scattering of the blue end is the robust conclusion here.

The model describes frequency-dependent scattering, not the detailed perceived color after solar-spectrum weighting and human visual response. It also explains why treating the electrons as free at visible frequencies would miss the bound-charge fourth-power limit.

A normal mode is a motion in which all generalized coordinates oscillate at one frequency with fixed amplitude ratios. Small oscillations about stable equilibrium in a time-independent conservative system give $\mathsf{M}\ddot{\mathbf{\theta }}+\mathsf{K}\mathbf{\theta }=0$, with positive mass matrix and symmetric stiffness matrix, allowing a generalized normal-mode eigenproblem.

For the three-to-one mass double-pendulum modes, the angles are absolute displacements from vertical, as printed. To quadratic order,For $\mathbf{\theta }=\mathbf{A}\cos \omega t$ and $q={\omega }^2\ell /g$, the determinant isThusThe low mode has the strings leaning together, with the lower angle twice as large; the high mode has opposite angular signs. Both reverse together after half their respective periods.

For the specified release, decompose $(\alpha ,\alpha )=(3\alpha /4)(1,2)+(\alpha /4)(1,-2)$. Thereforeand the lower mass's horizontal displacement isThis is the supremum of two incommensurate harmonic modes, conventionally described as the maximum-amplitude estimate. Strictly, ${\omega }_{+}/{\omega }_{-}=\sqrt{3}$ is irrational, so the opposite cosine extrema are approached arbitrarily closely but never coincide at a finite time. The quoted value is a supremum of the ideal linearized motion.

Consider an idealized planet that is a sphere of radius $R$, illuminated by parallel sunlight of irradiance $S$, measured on a plane perpendicular to the rays. Its intercepting projected area is $\pi R^2$, whereas its total surface area is $4\pi R^2$. Therefore the incoming power averaged per unit total surface area iswhere $\alpha$ is the broadband Bond albedo. The factor four is a sphere geometry result, not an atmospheric transmission factor. Local sunlight still varies over the illuminated hemisphere and vanishes on the dark side.

Now stipulate a single spatially uniform temperature $T$, no atmosphere, no internal heating, and emissivity one for outgoing thermal radiation. This is a thermal blackbody model using the Stefan–Boltzmann law; the prescribed solar-band reflection is treated separately, so the body is not a perfect absorber at all wavelengths. Let $C>0$ be the effective heat capacity per unit surface area, in $\mathrm{J}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$. Conservation of energy gives the ordinary differential equationHere $F$ is net heating per unit surface area, in $\mathrm{W}{\mathrm{m}}^{-2}$, and $\sigma$ is the Stefan-Boltzmann constant. Use $\sigma =5.670374419\times 10^{-8}\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-4}$, rounded from the exact constant reported by NIST. ^[ref] The uniform temperature is a model assumption: generally the surface average of $T^4$ differs from the fourth power of the surface-average temperature.

For fixed Bond albedo $\alpha$, radiative equilibrium requires $F(T_{\ast })=0$, givingWith the illustrative solar constant $S=1360\mathrm{W}{\mathrm{m}}^{-2}$, the global incoming average is $340\mathrm{W}{\mathrm{m}}^{-2}$:

Doubling albedo here lowers the equilibrium temperature by $33.231\mathrm{K}$, rather than halving it: the absorbed fraction changes from $0.7$ to $0.4$, and thermal radiation scales as $T^4$. These stipulated Bond albedos are whole-planet examples, not measurements of bare water and ice. The model predicts neither Earth's modern surface temperature nor its response to actual climate change; it omits the greenhouse effect, spatial heat transfer and seasons.

Each fixed-albedo equilibrium is a stable equilibrium: for $T_{\ast }>0$, the derivative is $F^{\prime }(T_{\ast })=-4\sigma T_{\ast }^3<0$. A small temperature increase produces net cooling, and a small decrease produces net heating. This is the simplest expression of the stabilizing Planck feedback. Allowing albedo to fall with temperature adds a competing term, as in ice-ocean albedo feedback.