On 11 December 1924, Edward Victor Appleton and Miles Aylmer Fulton Barnett perform an experiment: the Bournemouth broadcasting radio transmitter slowly changes wavelength; an Oxford radio receiver records alternating maxima and minima from ground wave-skywave wave interference. The inferred equivalent height is about . A separate comparison of simultaneous loop and vertical antenna signals tests the downward arrival direction. The E layer is established first; the higher F layer is found in winter 1926–27. These dates describe different findings.[ref]
Why change the wavelength deliberately? Ordinary fading alone could arise from changing emission, receiver drift or a moving propagation path. A controlled sweep predicts a specific progression of phase difference between two arrivals. The directional antenna comparison supplies an independent kind of evidence for an overhead contribution. Agreement supports an atmospheric return; fitting a distance alone would not establish a unique electron density profile. In the model below, comparable received amplitudes make the interference fringes easier to see. A practical repetition would also calibrate the frequency sweep, track radio transmitter output and receiver gain, and compare repeated sweeps with fixed-frequency intervals to detect time drift. These are controls inferred from the measurement model, rather than additional procedures attributed to the historical account.
The radio receiver does not need to resolve two pulses to detect two paths. If one path takes longer, the two contributions reaching the receiver during a slow frequency sweep carry slightly different emission frequencies and can produce a slow beat in the detected amplitude. At strictly fixed frequency with stationary paths, the arrivals instead have a fixed phase difference: there is no continuing beat merely because one path is longer. This distinction connects the historical observation to ionospheric interference fringe spacing.
Releasing the central bob of three identical spring-coupled simple pendulums from displacement excites only the symmetric normal modes. In displacement ratios, . Thus the outer coordinates are and the central coordinate is , where and . The antisymmetric mode is absent by reflection symmetry. The central modulation envelope is , so it never vanishes completely. Slow beats require a small frequency difference; the decomposition itself is exact in the linear model.
Comb self-referencing 2026-10-18
Comb self-referencing compares two parts of one optical frequency comb. Coherent broadening supplies an octave, a factor of two in frequency. Second-harmonic generation doubles a red tooth; a photodiode detects its beat with an original blue tooth.[ref] Applying the tooth formula gives
The doubled branch carries , while the original high-frequency branch carries : their difference retains exactly one offset.
Figure 1.
Read the tooth spacing and offset through separate electronic signals
. Original functional comb self-referencing scheme. Solid optical paths bring doubled red light and original blue light to one detector. Dashed electronic paths connect both readouts to a common reference and feedback controls. Actuators can affect both parameters. Details.
The displayed difference uses the order doubled tooth minus original tooth. An ordinary electronic frequency counter returns a nonnegative beat magnitude; spectral selection and known electronic shifts must identify whether it represents , , or a shifted version. Coherent broadening, overlapping photon polarization, spatial alignment and temporal overlap are required. Merely producing a broad patch of colored light does not establish a usable offset signal.
One beat magnitude does not tell us the integer . In the example, every proposed optical frequency produces the same signed beat with its corresponding tooth. A coarse independent optical measurement can select the integer: once the beat sign is known, an error bound below suffices if the candidate is selected by nearest distance.
There is also a useful two-spacing check. Keep the unknown laser and fixed, track the same tooth continuously, and increase by a stipulated . Its frequency rises by , so the signed beat falls from to . Hence
This inference fails if the unknown drifts appreciably, the offset changes unnoticed, the beat sign is lost, or the tracker jumps to another tooth. The assumed tooth motion stays well below half the spacing in this example. If also changes by a measured amount, use .
For a fixed verified tooth number, the first-order frequency error is . Errors in these readings can be correlated because they share a reference. A small beat uncertainty therefore does not imply the same absolute optical accuracy: the reference's calibration and the factor remain in the error budget. A wrong tooth number is a separate discrete error, not a small correction. The beauty of the comb is that two controlled parameters organize millions of optical reference points, while the unknown's position is read electronically.
Coupled oscillator Created 2026-10-20 Updated 2026-11-05
Coupled oscillators influence one another's motion through shared forces or other interactions. Two pendulums connected by a light spring provide a visible example: moving either bob changes the force on the other. For equal masses and lengths, the symmetric motion leaves the connecting spring unstretched, while the antisymmetric motion stretches it and has a higher frequency.
For a conservative linear system, the equations of motion take the form , where is the mass matrix and the elastic stiffness matrix. A normal mode has a fixed displacement pattern satisfying . Expanding the initial motion in normal modes separates the evolution. Displacing just one pendulum generally excites several normal modes, producing beats and visible transfer of energy. A driving force excites only the modes onto which its spatial pattern has a nonzero projection.
For two equal masses with equal wall springs released with displacements , their second displacement is , where is the mean of the two positive normal mode frequencies and their difference. The first nonzero time with is a fast turning point, not the slow maximum of its beat envelope. With , writing gives . Thus , using the first root of in . The envelope instead first reaches its maximum at .
Consider an original two-path model with a locally flat ground baseline and a symmetric equivalent return height . Replace the skywave by two straight segments traveling at the speed of light , and replace the ground wave by a straight segment of length . Their lengths and positive excess distance are
This is an equivalent geometric construction, not a metallic layer overhead. It neglects the curvature of Earth, refraction within the return region, ground-induced wave phase shifts and optical dispersion.
For received complex amplitudes and , absorb their constant relative wave phase into and choose their magnitudes nonnegative. An ideal detector reading proportional to squared received field has
Assume and are independent of frequency across the sweep, with slowly changing amplitudes. Consecutive maxima differ by in phase difference, so
The often-used approximation therefore carries a frequency-independent-path assumption. Counting complete fringe periods between and gives ; equivalently, when the vacuum wavelength is . Count whole periods, not maxima plus minima as though each were a whole period.
Choose and , with . Then , and . A sweep covers four full fringe periods. Conversely, the inferred equivalent height is
These chosen values illustrate the inversion; they are not a transcription of the 1924 measurements. With and , the normalized reading ranges from to : incomplete cancellation does not remove the distance information.
Figure 1.
Two arrival paths and the interference fringes produced by a controlled frequency sweep
. An original functional diagram distinguishes a curved skywave from a dashed equivalent straight-ray construction. The neighboring graph uses the stipulated nondispersive 150 km baseline and 100 km equivalent height, giving 3 kHz fringe spacing. The receiver detects the superposed arrivals. Details.
For a sweep with linearly changing frequency and two constant delays , the received instantaneous frequency difference has magnitude . In this model , giving a beat frequency . Equivalently, traversing each of the sweep crosses one fringe. This calculation requires a sufficiently slow sweep that propagation, tuning and amplitude changes can be treated consistently; it is not a beat between two independent transmitters.
More generally, measure the total phase difference , including propagation and any reflection or apparatus contribution. Define the differential group delay
For a locally monotone wave phase with slowly varying derivative, the local spacing is . If , this becomes
Here is a phase-equivalent optical path difference, not necessarily the geometric excess distance . A drifting ionosphere, several simultaneous skywaves, changing wave polarization, a frequency-dependent reflection wave phase or significant loss can invalidate the simple constant-spacing inversion.
Ammonia gas entered a vacuum chamber through fine channels; the selected beam entered a cavity with input and output waveguides.[ref] With external microwaves swept across a molecular resonance, the 1954 experiment read an increase in transmitted power; with sufficient beam emission, it sustained oscillation.[ref]
The 1955 readout mixed two maser signals with a local oscillator in a crystal detector, amplified them at an intermediate frequency, and combined them in a diode detector. An oscilloscope showed their beat frequency. Changing cavity tuning, source pressure or focuser field tested sensitivity; excessive detuning stopped oscillation.[ref] The beat tests relative frequency stability; it does not by itself establish an absolute frequency or an intrinsic linewidth independent of detector noise. Detuning also changes the field–molecule interaction, so the threshold model's constant cannot describe arbitrary tuning scans.
The 2000 experiment prepared optical pulses with a mode-locked laser and broadened them in an air-silica microstructure optical fiber. A wavelength-selective beam splitter separated red and blue light. A nonlinear crystal doubled the red branch; the blue branch received a known acousto-optic shift. A filter selected recombined light for an avalanche photodiode, supplying an electronic beat. Mirror translation controlled repetition frequency; mirror tilt controlled offset with coupled effects. Separate delayed-pulse interference tested programmed phase changes. Alternative envelope fits and unmatched optical reflections were examined as checks and limitations.[ref]
These readouts answer different questions. The electronic beat closes a feedback loop; changing its setpoint and checking the delayed-pulse phase tests whether that control reaches the optical waveform. Neither readout alone establishes the accuracy of an unknown optical frequency. That additionally requires a calibrated time reference, the correct tooth number and beat sign, and checks for missed cycles and drift. The ideal unshifted diagram above omits the experiment's known acousto-optic shift; its actual raw detector frequency must be corrected for that shift before identifying .
Mix an unknown single-frequency laser with a nearby comb tooth on a photodiode. Define the signed beat , so
The beat arises from the cross term in light intensity:
An electronic detector can follow this slow difference even though it cannot resolve either optical oscillation. The square-law expression also shows why a real-valued intensity trace by itself loses the sign of . A small known upward shift of the chosen tooth decreases when the unknown lies above it and increases when it lies below, provided the tooth assignment remains unchanged.
For an original example, stipulate
Then and . Choosing the opposite sign would give instead. The stipulated optical frequency is countable because the large optical value is expressed using the integer and three electronic frequencies; the diagram's two controls establish the comb, and this additional beat places the unknown against it.
Both comb readouts use one calibrated microwave reference. Feedback fixes and , determining each tooth once is identified. A frequency counter measures the slow beat, transferring the calibration to the optical frequency.[ref]
Write . The initial conditions on the normal mode coordinates are and . Solving the two harmonic oscillator equations and recovering the original angles gives the exact solution within the harmonic approximation:
Their angular velocities are and , which verify all four initial conditions.
To expose the beat, define and . Replacing by changes each angle by uniformly in time. The sum-to-product identities then give
The displayed angular velocities are exact in the linear system. The absolute and signed beat envelopes have leading amplitudes and . Because , at we have and . Hence
The first simple pendulum has negligible oscillation amplitude. Its kinetic energy is , while the second has . Thus, to the retained order, all instantaneous kinetic energy belongs to the second simple pendulum. The fast phase can make that kinetic energy small: the whole initial energy is transferred as the second simple pendulum's combined kinetic energy and gravitational potential energy, since
This is energy transfer after an impulse in two weakly coupled oscillators. It does not require .
At , , so the beat has returned the oscillation to the first simple pendulum:
The energy localization returns, although the fast oscillation phase need not equal its initial phase. The exact first minimum of the leading beat envelope is at , whose leading value is . Keeping the slow phase is essential: expanding at fixed and then using would lose the energy exchange. Although angles are retained only to order , their leading energies are of order .
Figure 1. . The exact linear solution at , with each pendulum's kinetic plus uncoupled gravitational potential energy divided by the initial total energy. The small remaining fraction is spring potential energy. Vertical lines mark and . Details.
In second-harmonic generation, the quadratic nonlinear polarization contains products of pairs of comb components. Their angular frequencies are
including when . Original comb components have angular frequencies , hence
When the doubled low-frequency region overlaps the original high-frequency region, choosing gives the beat at . More generally the nearest positive beat has angular frequency
which is for the stipulated small offset. For it is , provided the corresponding teeth are in the overlap.
In comb self-referencing, overlap the doubled red part and the original blue part on a photodiode. The beat measures the carrier-envelope offset frequency. Feedback on an intracavity dispersion control or pump setting can set modulo , giving . A phase-sensitive lock or calibrated approach to zero resolves the sign ambiguity of an unsigned beat measurement.[ref]
The coherent optical spectrum must span an octave, , with sufficient overlapping spectral strength after second-harmonic generation. This calls for very short broadband optical pulses, or coherent external spectral broadening. The two branches must preserve optical coherence, overlap in time and polarization, and have adequate phase matching; bandwidth alone does not guarantee a detectable beat.
Release from rest gives . At , and , so . Consequently
The frequency ratio is exactly two, so both motions repeat after . These are commensurate normal modes, rather than a pair of nearly equal frequencies producing slow beats.
For a central displacement in three freely ended coupled pendulums, the initial vector decomposes as . The antisymmetric second mode is absent. Hence
Let and . The outer displacement is , explicitly exhibiting an amplitude modulation. The middle displacement is , whose modulation envelope is
It ranges from to , so it does not vanish completely. These are beats between the two excited frequencies. A slow envelope well separated from the carrier requires ; in that limit . The exact two-mode formulas remain valid without that additional weak-coupling assumption.
Figure 1.
Normal-mode displacement ratios and beating after releasing the central pendulum
. Original eigenvectors and dimensionless time-domain solutions. The central initial displacement excites the uniform and highest-frequency modes, leaving the antisymmetric mode absent. Details.
Let and . The sum and difference identities give
For weak coupling, , so the slow sine and cosine envelopes produce beats. The first maximum of the second envelope is at , when the first envelope vanishes.
A literal first turning point of an initially stationary coupled oscillator is different from that slow envelope maximum. The first nonzero time at which the second mass is instantaneously at rest satisfies
Set and . Dividing the derivative by and expanding at of order one gives . The first positive root lies in , where , so
There is no earlier root: for , both terms in the exact derivative are positive when is sufficiently small. At the root the second displacement reaches its first local maximum. If the question instead intends the first envelope maximum, its time is . That is a slow amplitude-transfer time, and is not generally a time of zero instantaneous velocity. The distinction matters because both masses keep oscillating while the envelopes change.
Figure 1.
The first turning point and the much slower envelope transfer of coupled oscillators
. Original dimensionless solutions at weak nonzero coupling. An early-time panel marks the first zero of the second velocity, while long-time curves show the separate beat envelopes. Details.
As approaches , the driving motion couples resonantly to the pendulum. The particular-solution amplitude has denominator , so a steady harmonic-response formula alone is singular. For a specified finite initial state, the near-resonant motion instead develops beats, and its exact resonant limit has a growing envelope.
For example, with the vertical zero-angular-velocity initial condition above, taking the limit gives
The ideal undamped linear theory predicts amplitude growing linearly with time, not an infinite displacement at the first instant. In reality the motion leaves the small-angle regime, and damping and nonlinear changes of the pendulum's oscillation period alter the response. Below and above resonance, the sign of the harmonic response coefficient changes, corresponding to a phase reversal.
In SI units, the Maxwell equations in vacuum with prescribed charge density and current density are
When and , take the curl of the third equation and substitute the fourth. Using the curl of the curl identity and gives the electromagnetic wave equation:
For example, , is a plane electromagnetic wave provided . The electric field, magnetic field and propagation direction are mutually perpendicular.
The Maxwell equations are linear, so the superposition principle applies to fields. However, the Poynting vector is quadratic:
One may add the average wave intensities only when the cross terms vanish under the actual detector average. For co-propagating plane electromagnetic waves of parallel linear polarization, averaging over optical cycles but potentially resolving the difference frequency gives
If the detector integrates for , the last term acquires , with sinc function convention , evaluated at the midpoint phase. Thus suppresses beats, and the infinite-time average adds the wave intensities exactly for distinct fixed frequencies. Random mutually uncorrelated wave phases also remove the cross terms in an ensemble average; orthogonal linear polarizations remove them for co-propagating waves in polarization-insensitive detection. Merely having different frequencies does not remove a resolved beat signal: two nearly equal frequencies with fixed relative wave phase are a concrete counterexample. For different propagation directions, the relevant addition is of average Poynting vectors, or their components through the chosen detector surface.
A normal mode is a motion in which all coordinates have a fixed relative pattern and oscillate at one frequency. Near a stable equilibrium, let the quadratic kinetic and potential energies be and . With positive-definite mass matrix , the equations and harmonic ansatz give
The stiffness and mass matrices are symmetric. For distinct frequencies, their symmetry gives , so the mode vectors are orthogonal in the mass inner product. A basis can also be chosen in a degenerate eigenspace. Expanding decouples the equations into ; initial displacements and velocities determine the coefficients. Positive stiffness gives real frequencies, a zero eigenvalue can describe a free symmetry motion, and a negative one signals instability rather than an oscillatory mode.[ref]
For two equal masses between fixed walls with three equal springs, the stiffness matrix is and the mass matrix is . The patterns and have frequencies and . An initial displacement of one mass excites both modes and can produce beats.
For a uniform taut string fixed at both ends, the wave equation , with , has patterns and angular frequencies
The boundary conditions select the spectrum. Fourier expansion of arbitrary admissible initial data provides the mode amplitudes, illustrating the extension from a finite set of coordinates to a field.
For the normal modes of inductively coupled LC circuits, take two identical circuits with self-inductance , capacitance and mutual inductance , where . The capacitor charges obey and the exchanged equation. In-phase and opposite-phase charge patterns have
The positive inductance matrix for two identical coils is the physical counterpart of positive mass. Here inductance plays the role of the mass matrix. Mechanical, wave and electrical examples therefore share the same eigenvalue structure. Weak forcing excites modes according to its spatial overlap with each pattern; damping broadens the resonances, while strong nonlinearity can invalidate a fixed linear-mode decomposition.
Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form is
The vector describes the mode shape, while is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.
Let be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:
The mass matrix defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillations
A practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix . This guarantees real and a complete orthogonal basis in mass-weighted coordinates.
For distinct eigenvalues, the fact that and are real symmetric matrices implies
Thus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing decouples the motion into , with total energy . A positive-definite matrix gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.
Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.
For two equal masses , each attached to a fixed support by a spring of stiffness and coupled to the other by a spring of stiffness ,
The symmetric mode leaves the coupling spring unstretched and has . The antisymmetric mode stretches it and has . Superposing the two gives beats when their frequencies are close.
Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length , tension and mass per unit length , the wave equation and fixed endpoints give
Orthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.
A uniform chain of equal masses joined by springs , with spacing , has travelling normal modes with
Their quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion . These examples show how restoring forces and geometry determine the frequency spectrum.
Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to or can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.
For sinusoidal steady state with convention , complex impedance is the voltage phasor divided by the current phasor. Its modulus is the amplitude ratio and its argument is the voltage-current phase difference. Since ,
The original circuit is a series RLC circuit with output across the inductor. Its transfer function is
For the original capacitance, the outputs are at and at , a ratio . The current resonance is , so the circuit favors the first transmitter.
The resistor voltage is in phase with the current; the inductor voltage leads it by and the capacitor voltage lags by . Large reactive voltages therefore nearly cancel in their phasor sum. A large inductor voltage is consistent with a small supply voltage: it is inductor voltage magnification, not power gain in a passive circuit.
Reducing the capacitance raises to , between the two transmitters. Direct evaluation gives and , respectively. Their similar detunings give approximately equal amplitudes; they are not exactly equal, partly because the numerator of also depends on frequency.
The two-tone output of a linear resonant circuit, for coherent monochromatic sources of equal input amplitude and arbitrary fixed input phases, is
With , this is approximately
Thus a rapid oscillation has a beat envelope repeating in magnitude at . The signed envelope cosine has frequency . Unequal amplitudes leave a small nonzero envelope minimum; independently drifting transmitter phases make the beating drift. No new frequency is generated by this ideal linear circuit.
Figure 1.
Original calculated inductor-voltage gains and illustrative two-tone output, including the circuit's frequency-dependent phase shifts
. Details.