Hopf bifurcation Created 2026-09-24 Updated 2026-10-13
A Hopf bifurcation creates or destroys a small periodic orbit near an equilibrium point as a pair of complex conjugate eigenvalues crosses the imaginary axis. It explains how a resting state can give way to sustained oscillations: the linearization of a dynamical system supplies rotation and growth, while a nonlinear term can saturate that growth at a nonzero amplitude.
The simplest supercritical Hopf bifurcation is the exact Hopf normal form
Here is a complex number, is the bifurcation parameter, and is the angular frequency. In polar coordinates , the equations become and for . At , each turn shrinks towards the equilibrium point. At , small radii grow and large radii shrink towards , producing a limit cycle. For example, and give radius , period , and frequency cycles per unit time. The limiting angular frequency stays nonzero while the amplitude tends to zero: a tiny newborn oscillation can still repeat at a finite rate.
Figure 1.
Computed phase portraits across a supercritical Hopf bifurcation
. Exact orbits of the displayed Hopf normal form with , at . Arrows follow increasing time; the two starting radii are and . The green circle at has radius . At , nonlinear decay still attracts orbits, despite purely imaginary linear eigenvalues. These are dimensionless model calculations, not measured data. Details.
For a precise generic statement, consider a smooth function defining an autonomous differential equation on , with a smooth branch of equilibrium points and critical parameter . Require:
Then a locally unique branch of small periodic orbits, up to phase, emerges on the side . On the two-dimensional centre manifold, the branch attracts for and repels for . Attraction in the full phase space additionally requires all remaining eigenvalues to have negative real part; any transverse unstable direction makes the periodic orbit unstable there.[ref]
An eigenvalue-crossing bifurcation test alone cannot establish these nonlinear conclusions. For example, the purely linear equation has the same crossing but no isolated periodic orbit: at every nonzero radius is constant, and at none is. If the cubic coefficient vanishes, higher-order terms must be examined; a generalized Hopf bifurcation is one possible degeneracy. The usefulness of the generic theorem is that a large system can have its small oscillatory onset understood through one Hopf normal form, subject to these hypotheses and local limits.
Linear stability of principal-axis rotation Created 2026-10-03 Updated 2026-10-14
For , the transverse linear stability rates about steady rotation with signed component are
Axes one and three have transverse eigenvalues and , producing bounded oscillations. Axis two has transverse eigenvalues , giving growing and decaying exponential functions. A neutral longitudinal perturbation corresponds to a small change of . Generic intermediate-axis disturbances contain a growing component; a precisely chosen decaying direction is exceptional. For the normalized intermediate-axis model with , the values are , and .
Proof outline: the linearization of the Euler equations for a torque-free rigid body about the third principal axis gives transverse perturbations obeying
and hence
Cyclic permutation gives the other rates. For nonlinear Lyapunov stability of the extreme rotations, use the torque-free rigid-body invariants: the conserved quadratic forms and have strictly positive transverse coefficients about axes one and three respectively. They bound transverse angular velocity, while bounds its longitudinal change near nonzero rotation. Imaginary transverse eigenvalues alone would not establish this nonlinear conclusion.[ref]
The practical lesson is to choose a stable principal axis when using rotation to maintain an instrument's pointing direction. Both extreme principal moments of inertia are stable in the ideal torque-free model; if internal damping lowers rotational kinetic energy while approximately preserving angular momentum, the minimum energy instead favors the largest principal moment of inertia. A control system must account for torque, changing mass distributions and flexibility before applying this ideal result.
Magnetic dipole torque 2026-10-14
An externally imposed magnetic field that is nearly uniform over a rigid current loop produces the torque
The angle is measured between the magnetic dipole moment and the external magnetic field, not between the plane of the current loop and the magnetic field. The magnetic dipole potential energy uses a dot product, whereas the torque uses a cross product. The torque is largest when the magnetic dipole moment is perpendicular to the magnetic field; parallel alignment minimizes the potential energy.[ref]
For the sign convention and , positive rotation increases and
Thus at the torque turns the magnetic dipole moment toward the external magnetic field. Both parallel and antiparallel alignment have zero torque, but they are respectively a stable equilibrium and an unstable equilibrium. With no damping, a freely released current loop can undergo oscillation around the stable equilibrium; settling requires damping.
Figure 1.
Angular potential energy and signed torque of a magnetic dipole in an imposed field
. The ideal model has and . The diagram fixes the external magnetic field along positive , so the torque is clockwise for the illustrated positive angle. These curves are calculated, not measured data. Details.
In the twelve-turn example, take and . The predicted torque is and the magnetic dipole potential energy is . The unit check is ; torque and energy have the same physical dimensions but different physical meanings.
A short derivation starts with the magnetic force on a current-carrying wire, . In an exactly uniform magnetic field, the net force is , but its torque need not vanish: opposite sides produce a turning couple. Integrating gives . Integrating at fixed and gives the orientation potential energy. This is the mechanical orientation potential energy for a maintained electric current, with the external source held fixed; the energy accounting for field and power supply is a separate problem. A varying external magnetic field requires corrections across a finite current loop, and can also give a net force.
Malus's law 2026-10-13
For incident light with linear polarization, an ideal linear polarizer transmits
where is light intensity and is the angle between the incident electric field's line of oscillation and the polarizer's pass axis. The input and output are compared in the same medium. The classical explanation is an orthogonal projection: for a unit vector along the pass axis,
The dot product reduces the transmitted wave amplitude by . For a plane electromagnetic wave in the same medium, light intensity is proportional to squared electric field wave amplitude, giving Malus's law. For example, transmits one quarter of the incident light intensity, although the transmitted wave amplitude is one half. [ref].
Start with light already in linear polarization along the horizontal axis and let its light intensity be . Two ideal absorbing linear polarizers with pass axes and give . Insert a polarizer between them: the first leaves unchanged, the middle transmits with linear polarization along , and the last transmits with linear polarization along . Each angle in Malus's law is relative to the preceding output, so the final step uses :
Figure 1.
Two crossed polarizers transmit zero; inserting a middle polarizer transmits for initially horizontal linear polarization
. Pass axes are measured from horizontal. Each blue segment represents transmitted light intensity, and each inset shows the pass axis in the transverse plane. The two successive changes each halve light intensity by Malus's law. Details.
For these ideal absorbing polarizers, the rejected electric field component transfers energy to the filter. In the three-filter arrangement the middle absorbs the fraction , and the last absorbs , leaving transmitted. Conservation of energy holds at each step. Inserting the middle polarizer changes the linear polarization arriving at the last one; it does not create light. The effect illustrates why successive orthogonal projections cannot generally be replaced by a single comparison of the first and last axes.
The initially horizontal assumption matters. For unpolarized light of light intensity , the first ideal linear polarizer transmits : averaging over uniformly distributed incident orientations gives . The same sequence then transmits . This extra first-stage factor does not apply to the already horizontally polarized input in the figure.
Molecular geometry 2026-10-13
Why do methane, ammonia, and water have different shapes when each central chemical atom is surrounded by eight valence electrons? The answer separates the positions of atomic nuclei from the electron domains used in a VSEPR model.
Molecular geometry describes the three-dimensional arrangement of the atomic nuclei in a molecule, including bond lengths and bond angles. For a larger molecule, a shape name often describes just the neighbours of one central chemical atom. A Lewis structure records connections and electron bookkeeping; its flat layout does not fix molecular geometry.
Figure 1. . Original geometric constructions. The upper row is an ideal VSEPR model counting scaffold: orange LP markers identify lone pairs, not additional chemical atoms or measured electron density. The lower row shows only chemical atoms; ammonia and water use the sourced approximate bond angles discussed below. Sizes and bond lengths are arbitrary. Three-dimensional projections can change the apparent bond angles on the page. Dotted outlines are geometric guides, not extra chemical bonds. Details.
A static drawing also needs a physical interpretation. The atomic nuclei undergo oscillations; an arrangement at a minimum of potential energy and an experimentally averaged arrangement need not have identical bond lengths or bond angles. This distinction matters when comparing measurements with quantum mechanics calculations; see the IUPAC Gold Book: experimental bond-length conventions. Gas electron diffraction gives a concrete route from a laboratory signal to molecular geometry rather than inferring shape from a drawing alone.
Moment of inertia Created 2026-09-24 Updated 2026-10-14
The moment of inertia about a specified axis measures how a mass distribution resists changes of angular velocity about that axis. For point masses at perpendicular distances , or a continuous distribution with mass density ,
Its SI unit is . The distance is to the entire axis, not to one point on it: moving a point mass along the axis leaves its contribution unchanged. The moment of inertia depends on both the distribution and the chosen axis; it is not determined by total mass alone. About an axis fixed in an inertial frame, a rigid body with larger moment of inertia has smaller angular acceleration under the same net axial torque, and greater rotational kinetic energy at the same angular speed.[ref]
For a simple model, put four equal point masses of total mass on light rigid spokes at radius . Their moment of inertia is . Doubling quadruples the moment of inertia without changing . The illustrated placements and give and . Actual spokes and finite-sized weights contribute their own moments of inertia; the point mass model neglects them. Repositioning weights during rotation also breaks the fixed-distribution assumption, so it requires a separate work and angular momentum calculation.
Figure 1.
Same total mass, different radial placement and moment of inertia
. Original model diagrams, not experimental data. The upper pair has four equal point masses on negligible-mass spokes, at and from a perpendicular central axis. The lower pair compares a uniform disc with an ideal thin ring of the same mass and outer radius . The displayed moments of inertia determine rotational kinetic energy at a shared angular speed. The visible thickness of the ring and size of the point markers are only drawing conventions. Details.
The moment of inertia of a uniform disc follows by dividing it into concentric annuli: and . For an ideal thin ring every mass element is at , so . A ring of finite radial thickness instead uses the moment of inertia of a uniform annulus, , with inner radius . A moment of inertia is thus a compact summary of radial placement rather than a complete description of shape. For several joined pieces, their moments of inertia add only when each is calculated about the same axis.[ref]
A physical pendulum makes the moment of inertia measurable through timing. Martín Monteiro, Cecilia Stari, Cecilia Cabeza and Arturo C. Martí used a metal bar with suspension holes apart and an attached smartphone. They report bar length , bar mass and phone mass . Changing the suspension hole changes the pivot axis. Their readouts were phone gyroscope angular velocity records, analyzed by autocorrelation or a fit to a sinusoid, and an independent photogate; the logged records for the sinusoid fits lasted . They observed agreement between the three timing methods and the geometry-based model, including a minimum of oscillation period as the suspension point varied. These are reported apparatus details and observations, not measurements made for this article.[ref]
To infer a moment of inertia from such a measurement, determine the total mass including every attached part, locate the assembly's centre of mass , and measure its distance below the horizontal pivot . Record successive crossings in the same direction, or the elapsed time for complete oscillations, giving oscillation period . The physical pendulum model with a fixed pivot, a uniform gravitational field, negligible friction and the small-angle approximation gives
The crossing times, masses and geometric distances are the measured inputs; is an inference conditional on this model, while an integral over the assumed shape supplies an independent model prediction. The parallel axis theorem converts the inferred pivot moment of inertia to one about a parallel axis through . The controls and interpretation of that comparison are described there. An axial measurement determines one scalar moment of inertia, not the complete inertia tensor needed for general rotation.
Negative feedback 2026-10-13
Negative feedback closes a causal loop whose net effect opposes an initiating change. In signal transduction, activation followed by a downstream inhibition of the earlier activating stage is one example of negative feedback. Loop sign alone does not establish asymptotic stability or complete signaling adaptation; see the worked signaling feedback example.
A useful biological example is chemotaxis feedback regulation: an Escherichia coli biological cell adjusts its response to a sustained attractant instead of remaining at its immediate response. Increased CheA activity promotes CheB phosphorylation; CheB reverses receptor methylation, opposing their activation. After attractant addition, reduced activity instead favors receptor methylation by CheR, restoring activity. This is a return from receptor activity to the machinery regulating that activity. The response need not return exactly at every dose. [ref]
A motivating measurement is the chemotaxis FRET assay. CheY and CheZ carrying yellow and cyan fluorescent tags were observed in immobilized Escherichia coli biological cells on a coverslip in a flow chamber. Adding and removing attractant changed the yellow-to-cyan emission ratio, reporting the association of phosphorylated CheY with CheZ. This measures a molecular readout upstream of motor behavior. Comparisons with biological cells lacking CheR, CheB, or both tested the dependence on receptor modification; these perturbations also change starting activity and sensitivity, so they are not pure changes of one feedback gain. [ref]
The later chemotaxis feedback regulation experiment used only the Tar receptor and nonmetabolizable methylaspartate, supplied and removed by flow. Long recordings separated slow recovery from incomplete signaling adaptation; correction for bleaching addressed drift in the optical readout. Activity recovered incompletely after sufficiently large additions, with detectable imprecision already at in that preparation. Restricting the receptor type and preventing consumption of the stimulus help interpret the perturbation, but limit transfer to other receptors and conditions. The readout infers CheA activity through CheY-CheZ association, rather than directly measuring every molecular step. [ref]
Negative feedback can also overshoot or sustain oscillations if its return arrives too late. The worked signaling feedback example gives the contrasting instantaneous case; the delayed negative feedback model below isolates the effect of timing. Loop sign describes causation, while asymptotic stability describes what happens to displacements over time.
Figure 1.
Negative feedback with an instantaneous or delayed return
. Original synthetic delay differential equation modes, not measured chemotaxis responses. All variables and time are normalized. The same negative return coefficient is used throughout; only its delay changes. Each curve has its own matching history. A growing oscillation signals instability of this linear mathematical model, not a bounded biological rhythm. Details.
Normal mode Created 2026-09-24 Updated 2026-10-14
A normal mode, also called an eigenmode, is a pattern of free oscillation in a linear system whose relative displacements stay fixed while a single temporal factor changes their size and sign. For an undamped stable vibrational mode,
Here is the mode shape, sets the amplitude, is the eigenfrequency expressed as an angular frequency, and is the oscillation phase. Every snapshot is a scalar multiple of the same mode shape; neighbouring regions may move in opposite directions. Normal modes of a double pendulum show that normal modes also describe coupled discrete systems, beyond strings. Their shapes are eigenvectors, rather than spatial eigenfunctions. [ref]
For the ideal wave equation on a string, take a uniform string of length with constant tension and linear mass density . Assume small transverse slope, negligible bending stiffness, no damping, and fixed endpoints. Its wave speed is , and the nonzero normal modes are
The angular frequency is in radians per second; the frequency is in hertz. The wavelength is . The mode shape has loops separated by nodes of a standing wave, including both fixed endpoints. [ref]
As a worked example, choose , , and . Then and , , . Their oscillation periods are , , and . These are deliberately simple model parameters, rather than measurements of a musical instrument. The plot divides displacement by so that the three mode shapes can be compared, and uses each mode's own oscillation period for its snapshots.
Figure 1.
First three ideal fixed-string normal modes at fractions of their own periods; black dots mark permanent nodes
. Details.
For , at every point has zero displacement, but most points have nonzero velocity. Only zeros of the mode shape are permanent nodes of a standing wave. At the original mode shape has reversed sign, and at it returns. The synchronized conversion between kinetic energy and potential energy is shown in normal-mode energy exchange.
The short derivation uses separation of variables: substituting into the wave equation gives and . The fixed-end boundary conditions remove the cosine spatial term and require , hence and . Different geometries and boundary conditions give different eigenfunctions; a sine mode shape is a property of this uniform fixed string, not a universal form of an eigenmode.
The usefulness is that a complicated vibration becomes a collection of harmonic oscillators. For sufficiently regular fixed-end initial data, the Fourier sine series gives
The initial displacement and velocity determine and . This is the superposition principle in a Fourier sine basis, and modal energy distribution measures how strongly each normal mode is excited. The striking simplification is that each component keeps its identity while their sum can continually change shape. [ref]
These independent normal modes belong to the linear model. At larger slopes, changing tension introduces nonlinear effects and can couple the modes. A linearly damped string still admits a Fourier sine basis, but its temporal factors decay and need not oscillate; its energy dissipation identity for a linearly damped string replaces conservation of energy. More general damping need not leave undamped mode shapes independent.
Parallel axis theorem Created 2026-09-24 Updated 2026-10-14
The parallel axis theorem transfers a moment of inertia from an axis through the centre of mass to a parallel axis through :
Here is the total mass and is the shortest perpendicular distance between the two axes. Both moments of inertia concern the same mass distribution and the same axis direction. The axis for must pass through the centre of mass; may refer to an axis outside the body. No uniform-density or planar-shape assumption is needed. Within any family of parallel axes, the centre of mass axis minimizes the moment of inertia.[ref]
A short proof isolates why the centre of mass matters. Choose parallel axes perpendicular to an plane, put at the origin, and place at . Then
because in centre of mass coordinates. Individual mass elements can move closer to the new axis, but their weighted cross terms cancel. For arbitrary parallel axes and , use , where each distance is measured from the parallel centre of mass axis; adding times the squared separation of and is generally incorrect. A displacement along the axis has . Changing the axis direction instead requires the inertia tensor or, in its applicable planar setting, the perpendicular axis theorem.
For a uniform disc, mass and radius , the central perpendicular moment of inertia is . Shifting that axis to the rim gives and , three times . In the computed , example, , and . At angular speed the rim-pivoted disc has rotational kinetic energy . Its centre of mass now moves, so already includes that translation; adding a separate centre of mass translation term to it would count that contribution twice.
Figure 1.
Parallel axis theorem: a uniform disc's central and rim axes
. Original geometric and numerical model, not measured data. Both axes are perpendicular to the disc; the separation is . The bar chart separates from the additional in , normalized by . All parts of the disc remain in the calculation, including those nearer to the shifted axis. Details.
The physical pendulum measurement described under moment of inertia provides a practical check: infer at several pivot positions, then subtract to obtain . Keep all pivot axes parallel and measure to the centre of mass of the entire moving assembly. The same assembly should yield the same within measurement precision. If a smartphone is attached, its mass and its own moment of inertia belong in both the geometric prediction and the centre of mass calculation. The published bar-and-phone study explicitly adds both contributions using the parallel axis theorem.[ref]
For a careful replication, release small swings without a push, repeat timings at several small amplitudes, and check that the extrapolated oscillation period is stable as amplitude decreases. Time many complete oscillations and repeat the run; two crossings of a photogate per cycle must not be mistaken for two cycles. Inspect decay of amplitude and compare early and late timings to detect significant pivot friction or air drag. Mark the bar and phone mount to detect relative slip, tighten the attachment, and check that the pivot support does not rock or migrate between holes. Repeat with an independent photogate or video readout to check the phone's timing. These are proposed controls for the replication, not additional results reported by the authors. The ideal model excludes appreciable damping, large swings and flexing; subtraction of also magnifies relative errors when is much smaller than .
The oscillation period need not increase monotonically with pivot distance: the shift raises the moment of inertia but also increases the gravitational restoring torque. For a computed slender uniform bar model of length , and
This expression has its minimum at . With , and stipulated , the prediction is ; it is not a recorded timing or a prediction for the bar with a phone attached. The competing terms explain why a timing experiment can reveal a mass distribution rather than merely its weight.
Phase difference 2026-10-13
The difference of two wave phases at the same place and time, considered modulo . It determines whether two same-frequency oscillations reinforce, cancel, or are in phase quadrature.
Polarization (waves) 2026-10-13
Polarization of waves describes the direction and time-dependent geometry of a transverse wave's vector during its oscillation at a fixed position. In three spatial dimensions, a transverse wave can have two independent directions of oscillation, both perpendicular to its direction of propagation. For polarization of an electromagnetic wave, the relevant vector is the transverse electric field; for a mechanical transverse wave, it can be the transverse displacement. A purely longitudinal wave has no choice between these transverse directions, and the orientation of a static electric field is not polarization of waves.
Choose transverse axes in a Cartesian coordinate system. A monochromatic wave with a fixed phase difference has the real electric field
Here are wave amplitudes, is the angular frequency, and is the phase difference of the oscillation relative to the oscillation. The vector's tip traces its path during the wave period . Both plotted coordinates are real electric field components; neither axis represents the imaginary part of a complex amplitude.
For nonzero wave amplitudes, or modulo gives linear polarization: the vector moves back and forth along a fixed line. Equal nonzero wave amplitudes and modulo give circular polarization: the tip traces a circle of constant radius. Unequal nonzero wave amplitudes with give a noncircular ellipse whose axes coincide with . More generally, nonzero wave amplitudes and give elliptical polarization, with circular polarization as a special case. If exactly one wave amplitude vanishes, the nonzero wave has linear polarization for every ; if both vanish, there is no nonzero wave to assign a polarization of waves. [ref].
For example, take , so the wave period is . The choices , and , with electric field wave amplitudes in , produce respectively the line segment , the circle and the ellipse , when coordinates are measured in . At , their tips are respectively , and . The last two share a point without having the same path.
Figure 1. . All panels use the same coordinate scale. Orange arrows follow increasing wave phase ; the dot marks . The blue vector is the electric field at . These are paths at one fixed position, rather than spatial drawings of a propagating wave. Details.
Resolving the electric field into two transverse components makes polarization of waves useful: the same frequency and wave amplitudes can describe different paths when their phase difference changes. An ideal linear polarizer selects the component along its pass axis, letting its orientation probe the direction of linear polarization through Malus's law. A fixed pair of wave amplitudes and a fixed phase difference describe one stable path; unpolarized light requires a statistical description rather than one fixed ellipse.
Reduced mass Created 2026-09-24 Updated 2026-10-14
The reduced mass is the effective mass for the two-particle relative coordinate of two bodies with constant positive masses :
Both bodies can respond to their mutual force, so their relative acceleration is greater than if one were fixed. Consequently . This explains why using either individual mass, or their sum, generally gives the wrong internal oscillation rate. The same reduced mass describes a two-body problem, a free two-mass spring, and the approximate molecular vibration of a diatomic molecule.[ref]
For equal masses , . Two bodies therefore have : in their stationary centre of mass frame a relative speed means individual speeds and total kinetic energy . With and , instead . At the same relative speed, the individual speeds are and , and the total kinetic energy is . These are calculations for stipulated masses and velocities.
For fixed and , . The heavy-partner limit gives from below; the relative error incurred by replacing with is . At this is . Treating the heavy body as fixed is thus a controlled approximation to internal motion, rather than an exact consequence of its being heavier.
RLC circuit Created 2026-10-05 Updated 2026-10-14
An RLC circuit combines electrical resistance , inductance , and capacitance . The capacitor stores energy in its electric field; the inductor stores energy in its magnetic field. Their exchange can produce an oscillation, while Joule heating removes energy. A series RLC circuit gives a particularly clear example: changing turns ringing into a direct return to equilibrium, and changing or tunes a resonance.
These circuits connect two useful viewpoints. In time, they describe how stored energy is released after switching. In frequency, they select a range of signals: the resonant series RLC response explains radio tuning and simple filters. The same mathematics describes a damped harmonic oscillator, so insight about mechanical ringing transfers directly to an electrical circuit. [ref]; [ref].
Rüchardt experiment 2026-10-14
The Rüchardt experiment infers heat capacity ratio from the period of an oscillation of a piston supported by trapped gas. Use a vertical precision cylinder with a close-fitting movable piston, known moving mass , face area , trapped equilibrium volume , and local atmospheric pressure . With negligible static friction, force balance gives the absolute gas pressure , where is gravitational acceleration. Displace the piston slightly and release it: compression raises pressure, and expansion lowers it, producing a restoring force. The oscillation is a mechanical probe of thermal response, without separately measuring two heat inputs. [ref]
Let be upward displacement from equilibrium, so . For constant heat capacity ratio and a reversible adiabatic process, the net upward force, including the ambient load and weight, is
The linear approximation requires . Applying Newton's second law gives simple harmonic motion with
Here is the undamped angular frequency and the undamped oscillation period, rather than temperature. Record repeated cycles, measure independently, and include all connected trapped volume. An illustrative model with , , , , and gives and . These stipulated values are not measurements.
Figure 1. . Original schematic and analytical ideal gas curves. Upward piston displacement increases volume; the blue curve uses constant , its dashed tangent gives the small-amplitude model, and the orange curve is the isothermal process limit. The shaded central interval illustrates a small fractional volume change, without specifying an experimental acceptance threshold. No measured waveform is shown. Details.
In Shearwood and Sloan's implementation, piston-mounted magnets induce an electromotive force in a stationary coil through Faraday's law of induction; an oscilloscope records the signal. They observed faster oscillations at smaller trapped volume values and stronger damping at small cylinder volume values. The signal tracks velocity locally, so its amplitude is not directly piston displacement. Calibration also identified inlet dead space beyond the graduated cylinder volume. [ref]
Repeat at smaller initial amplitudes to test the linear approximation, and compare several values. Heat transfer weakens the restoring force toward the isothermal process limit, biasing the inference toward . Check leakage through resting-piston drift and before-and-after trapped volume readings. For viscous damping with envelope , the measured angular frequency obeys ; this correction does not remove thermal leakage. Internal equilibration must be faster than an oscillation period, with heat transfer and leakage slow over a run. Large amplitude or sticking friction require a richer model. [ref]