Skip to content

Linewidth and Coherence

A laser can oscillate in one spatial and longitudinal mode without being perfectly monochromatic. The selected mode fixes a carrier-frequency neighborhood and field pattern; fluctuations of its phase make the emitted line finite. A laser linewidth is therefore not merely a property of a passive cavity resonance. It is an operational summary of temporal coherence, stochastic dynamics, observation time, and measurement method.

The shortest reliable statement is

phase fluctuations⟷decay of g(1)(τ)⟷finite optical spectrum.\text{phase fluctuations} \longleftrightarrow \text{decay of }g^{(1)}(\tau) \longleftrightarrow \text{finite optical spectrum}.

For ideal white frequency noise, this chain gives an exponential coherence envelope and a Lorentzian line. Real lasers also have drift, flicker noise, mechanical resonances, modulation sidebands, mode hops, and servo features. Then one full width at half maximum cannot encode the complete behavior.

Required background. Optical Cavities fixes the passive decay convention; Rate-Equation Lasers supplies the mean operating point; Correlation Functions supplies field correlations and spectra. This page treats the fluctuating phase of a selected laser mode.

Optical linewidth and frequency-noise conventions

Section titled “Optical linewidth and frequency-noise conventions”

Factors of two and 2π2\pi cause more linewidth disagreements than new physics. Unless stated otherwise:

  • the positive-frequency field is written E(+)(t)=A(t)e−i[2πν0t+ϕ(t)]\mathcal E^{(+)}(t)=A(t)e^{-i[2\pi\nu_0t+\phi(t)]};
  • ν\nu and Fourier offset ff are ordinary frequencies in hertz;
  • ω=2πν\omega=2\pi\nu is an angular frequency in radians per second;
  • Δν\Delta\nu means an optical power-spectrum FWHM in hertz;
  • κ\kappa is the intracavity photon-number decay rate, defined by n˙=−κn\dot n=-\kappa n when gain is absent;
  • Sν(2)(f)S_\nu^{(2)}(f) is a two-sided classical frequency-noise PSD;
  • Sν(1)(f)=2Sν(2)(f)S_\nu^{(1)}(f)=2S_\nu^{(2)}(f) for f>0f>0 is its one-sided version;
  • PSD formulas assume stationary noise unless a finite record or explicit low-frequency cutoff is declared.

The instantaneous frequency fluctuation associated with the phase convention above is

δν(t)=12πdϕdt.\delta\nu(t) = \frac{1}{2\pi} \frac{d\phi}{dt}.

Changing the sign in the field convention changes the sign of δν\delta\nu, but not any autocorrelation or PSD used below.

For one selected polarization and detected spatial mode, define

G(1)(τ)=⟨E(−)(t)E(+)(t+τ)⟩.G^{(1)}(\tau) = \left\langle \mathcal E^{(-)}(t) \mathcal E^{(+)}(t+\tau) \right\rangle.

For a stationary beam, this depends only on the delay τ\tau. Its normalized form is

g(1)(τ)=G(1)(τ)G(1)(0).g^{(1)}(\tau) = \frac{ G^{(1)}(\tau) }{ G^{(1)}(0) }.

The normalization gives g(1)(0)=1g^{(1)}(0)=1 for nonzero mean intensity and ∣g(1)(τ)∣≤1|g^{(1)}(\tau)|\le 1. In a balanced delayed interferometer, the fringe visibility is ∣g(1)(τ)∣|g^{(1)}(\tau)| when polarization, spatial mode, and intensities are matched. Coherence is therefore a relation between two field samples, not a statement that an absolute optical phase is known.

If amplitude noise is negligible on the timescale of interest,

E(+)(t)≃E0e−i[2πν0t+ϕ(t)],\mathcal E^{(+)}(t) \simeq \mathcal E_0 e^{-i[2\pi\nu_0t+\phi(t)]},

and

g(1)(τ)=e−i2πν0τ⟨e−i[ϕ(t+τ)−ϕ(t)]⟩.g^{(1)}(\tau) = e^{-i2\pi\nu_0\tau} \left\langle e^{-i[\phi(t+\tau)-\phi(t)]} \right\rangle.

Only phase increments appear. The absolute phase may diffuse without bound while the increment statistics remain stationary.

Use the ordinary-frequency Wiener–Khinchin pair

SE(ν)=∫−∞∞dτ ei2πντG(1)(τ),G(1)(τ)=∫−∞∞dν e−i2πντSE(ν).\begin{aligned} S_E(\nu) &= \int_{-\infty}^{\infty} d\tau\, e^{i2\pi\nu\tau} G^{(1)}(\tau), \\ G^{(1)}(\tau) &= \int_{-\infty}^{\infty} d\nu\, e^{-i2\pi\nu\tau} S_E(\nu). \end{aligned}

An overall normalization depends on whether SES_E is a power, photon-flux, or field spectrum. The line shape and widths do not. For an isolated stationary line, the normalized spectrum

L(ν)=SE(ν)∫SE(ν′) dν′\mathcal L(\nu) = \frac{S_E(\nu)} {\int S_E(\nu')\,d\nu'}

has units of inverse hertz and integrates to one.

A Lorentzian with FWHM Δν\Delta\nu is

LL(ν)=1πΔν/2(ν−ν0)2+(Δν/2)2.\mathcal L_{\mathrm L}(\nu) = \frac{1}{\pi} \frac{\Delta\nu/2}{ (\nu-\nu_0)^2+(\Delta\nu/2)^2 }.

Its normalized field correlation is

gL(1)(τ)=e−i2πν0τe−πΔν∣τ∣.g_{\mathrm L}^{(1)}(\tau) = e^{-i2\pi\nu_0\tau} e^{-\pi\Delta\nu|\tau|}.

Thus the 1/e1/e decay time of the field-correlation amplitude is

τ1/e=1πΔν.\tau_{1/e} = \frac{1}{\pi\Delta\nu}.

The integral convention used on the general correlation page,

τc(2)=∫−∞∞dτ ∣g(1)(τ)∣2,\tau_c^{(2)} = \int_{-\infty}^{\infty} d\tau\, |g^{(1)}(\tau)|^2,

happens to give the same numerical value for a Lorentzian:

τc(2)=1πΔν.\tau_c^{(2)} = \frac{1}{\pi\Delta\nu}.

That equality is special to the exponential envelope and these integration limits. A one-sided integral, an intensity-correlation width, or a different threshold gives another number.

A normalized Gaussian line with rms width σν\sigma_\nu is

LG(ν)=12πσνexp⁡[−(ν−ν0)22σν2].\mathcal L_{\mathrm G}(\nu) = \frac{1}{ \sqrt{2\pi}\sigma_\nu } \exp \left[ -\frac{(\nu-\nu_0)^2}{2\sigma_\nu^2} \right].

Then

gG(1)(τ)=e−i2πν0τe−2π2σν2τ2,g_{\mathrm G}^{(1)}(\tau) = e^{-i2\pi\nu_0\tau} e^{-2\pi^2\sigma_\nu^2\tau^2},

and

ΔνG=22ln⁡2 σν.\Delta\nu_{\mathrm G} = 2\sqrt{2\ln2}\,\sigma_\nu.

Slow Gaussian-distributed frequency excursions often create a Gaussian core, while white frequency noise creates Lorentzian wings or a Lorentzian core. Their convolution is Voigt only when the corresponding broadening processes are independent and have the required stationary distributions.

A coherence time becomes a distance only after the propagation geometry is specified. For group velocity vgv_g,

ℓc=vgτc.\ell_c = v_g\tau_c.

In a Michelson interferometer, a mirror displacement xx changes the optical path by approximately 2x2x, so the mechanical displacement associated with a delay τ\tau is x=vgτ/2x=v_g\tau/2. Quoting c/Δνc/\Delta\nu as “the coherence length” without a time definition and path convention can therefore differ from another correct convention by factors of π\pi or two.

Three linewidths that must not be conflated

Section titled “Three linewidths that must not be conflated”
Object What has a width? Primary mechanism
passive cavity driven cavity power response photon escape and internal loss
material transition absorption or emission response radiative, collisional, Doppler, power, and other broadening
laser output oscillator field spectrum phase noise plus technical fluctuations

A laser output can be far narrower than both its gain bandwidth and its passive-cavity resonance. Gain replenishes the decaying field amplitude, but it does not restore a unique absolute phase after random phase perturbations.

Above threshold, gain saturation supplies a restoring tendency for amplitude: if the field grows, it depletes inversion and reduces net gain. A free-running autonomous laser has continuous phase symmetry. Rotating the steady complex amplitude by a constant phase produces an equally valid solution, so there is no analogous restoring force for absolute phase.

A small random quadrature kick has a radial part and a tangential part in complex-amplitude space. The radial part is damped; the tangential part shifts phase. Repeated tangential kicks accumulate into diffusion.

Write the phase stochastic differential equation as

dϕ(t)=2Dϕ dWt,d\phi(t) = \sqrt{2D_\phi}\,dW_t,

where dWtdW_t is a Wiener increment with

⟨dWt⟩=0,⟨dWt2⟩=dt.\langle dW_t\rangle=0, \qquad \langle dW_t^2\rangle=dt.

The phase increment is Gaussian with

⟨[ϕ(t+τ)−ϕ(t)]2⟩=2Dϕ∣τ∣.\left\langle [\phi(t+\tau)-\phi(t)]^2 \right\rangle = 2D_\phi|\tau|.

For any zero-mean Gaussian variable XX,

⟨e−iX⟩=e−⟨X2⟩/2.\langle e^{-iX}\rangle = e^{-\langle X^2\rangle/2}.

Therefore

g(1)(τ)=e−i2πν0τe−Dϕ∣τ∣.g^{(1)}(\tau) = e^{-i2\pi\nu_0\tau} e^{-D_\phi|\tau|}.

Comparison with the Lorentzian Fourier pair gives

Dϕ=πΔν,Δν=Dϕπ.D_\phi = \pi\Delta\nu, \qquad \Delta\nu = \frac{D_\phi}{\pi}.

Some authors define their diffusion constant by ⟨(Δϕ)2⟩=D∣τ∣\langle(\Delta\phi)^2\rangle=D|\tau| rather than 2D∣τ∣2D|\tau|. The physics is unchanged; the symbol differs by two.

Define the two-sided frequency-noise PSD by

Sν(2)(f)=∫−∞∞dτ e−i2πfτ⟨δν(t+τ)δν(t)⟩.S_\nu^{(2)}(f) = \int_{-\infty}^{\infty} d\tau\, e^{-i2\pi f\tau} \langle \delta\nu(t+\tau)\delta\nu(t) \rangle.

Since δν=ϕ˙/(2π)\delta\nu=\dot\phi/(2\pi),

Sν(2)(f)=f2Sϕ(2)(f),S_\nu^{(2)}(f) = f^2S_\phi^{(2)}(f),

away from singular contributions at f=0f=0. Equivalently,

Sϕ(2)(f)=Sν(2)(f)f2.S_\phi^{(2)}(f) = \frac{S_\nu^{(2)}(f)}{f^2}.

Suppose the two-sided frequency-noise PSD is white:

Sν(2)(f)=h0.S_\nu^{(2)}(f)=h_0.

Then its time correlation is h0δ(τ)h_0\delta(\tau), and integration over a delay gives

⟨[ϕ(t+τ)−ϕ(t)]2⟩=4π2h0∣τ∣.\left\langle [\phi(t+\tau)-\phi(t)]^2 \right\rangle = 4\pi^2h_0|\tau|.

Hence

Dϕ=2π2h0,Δν=2πh0.D_\phi = 2\pi^2h_0, \qquad \Delta\nu = 2\pi h_0.

For the one-sided white plateau H0=2h0H_0=2h_0, the same result is

Δν=πH0.\Delta\nu = \pi H_0.

This useful formula is valid only for the declared one-sided PSD in Hz2/Hz\mathrm{Hz^2/Hz}. Applying it to an angular-frequency PSD or a two-sided PSD without conversion gives a wrong factor.

For stationary, zero-mean Gaussian frequency noise,

Δϕ(τ)=2π∫tt+τδν(t′) dt′.\Delta\phi(\tau) = 2\pi \int_t^{t+\tau} \delta\nu(t')\,dt'.

Using a one-sided PSD gives

12⟨[Δϕ(τ)]2⟩=2∫0∞df Sν(1)(f)sin⁡2(πfτ)f2.\frac12 \left\langle [\Delta\phi(\tau)]^2 \right\rangle = 2 \int_0^\infty df\, S_\nu^{(1)}(f) \frac{ \sin^2(\pi f\tau) }{ f^2 }.

Therefore

g(1)(τ)=e−i2πν0τ×exp⁡[−2∫0∞df Sν(1)(f)sin⁡2(πfτ)f2].g^{(1)}(\tau) ={} e^{-i2\pi\nu_0\tau} \times \exp \left[ -2 \int_0^\infty df\, S_\nu^{(1)}(f) \frac{\sin^2(\pi f\tau)}{f^2} \right].

The optical line follows by Fourier transformation. This relation is more fundamental than assigning each PSD a single “integrated linewidth.” It also shows why low-frequency noise is dangerous: the phase response contains 1/f21/f^2.

For practical one-sided frequency-noise data, the β-separation line is

Sν,β(f)=8ln⁡2π2 f.S_{\nu,\beta}(f) = \frac{8\ln2}{\pi^2}\,f.

Noise above this line has a high modulation index and contributes strongly to an approximate FWHM. With record duration T0T_0, define

A(T0)=∫1/T0∞df Sν(1)(f) Θ[Sν(1)(f)−Sν,β(f)].\begin{aligned} A(T_0) = \int_{1/T_0}^{\infty} df\, S_\nu^{(1)}(f) \,\Theta \left[ S_\nu^{(1)}(f) -S_{\nu,\beta}(f) \right]. \end{aligned}

The estimate is

Δνβ≃8ln⁡2 A.\Delta\nu_\beta \simeq \sqrt{8\ln2\,A}.

This is a useful diagnostic, not an exact replacement for Fourier transforming g(1)g^{(1)}. Its explicit lower limit exposes the observation-time dependence caused by flicker or random-walk frequency noise. Noise below the β line can still make spectral wings, sidebands, and phase errors important to an application.

A white one-sided frequency-noise plateau, its exponential field coherence, and the resulting Lorentzian power spectrum with a correctly marked half-maximum width.

A white one-sided frequency-noise plateau H0H_0 gives phase diffusion Dφ=π2H0D_φ=π^2H_0 and optical FWHM Δν=πH0Δν=π H_0. The middle panel uses delay in units of 1/Dφ1/D_φ. The lower panel uses offset in units of the FWHM, so its half-maximum points lie at ±1/2±1/2. These exact curves describe the ideal Gaussian phase-diffusion model; colored noise requires the PSD integral above.

A deterministic sinusoidal phase

ϕ(t)=βsin⁡(2πfmt)\phi(t) = \beta\sin(2\pi f_mt)

produces discrete sidebands with Bessel-function weights. It does not by itself produce irreversible phase diffusion. Random modulation, unresolved sidebands, or a finite instrument resolution may make the measured trace look broadened, so line broadening must be distinguished from coherent modulation.

Stimulated emission adds field in phase with the occupied mode, whereas spontaneous emission into that mode carries an unpredictable quadrature. For a large coherent amplitude of radius nˉ\sqrt{\bar n}, a transverse fluctuation of fixed scale produces a phase kick proportional to 1/nˉ1/\sqrt{\bar n}. The phase variance per kick is therefore proportional to 1/nˉ1/\bar n.

This gives the robust scaling

quantum-limited linewidth∝noise injection rateintracavity photon number.\text{quantum-limited linewidth} \propto \frac{ \text{noise injection rate} }{ \text{intracavity photon number} }.

Since output power is proportional to photon number at fixed output coupling, the linewidth narrows approximately as 1/Pout1/P_{\mathrm{out}}. The original Schawlow–Townes argument captured this central idea. This random-phasor picture is heuristic. In the conventional quantum model used below, gain and loss make equal contributions to the phase-diffusion limit, as shown by Wiseman; attributing the entire result to spontaneous gain can give an incorrect factor. The division also depends on the stochastic representation, while the measured linewidth does not.

Consider an ideal single-mode, far-above-threshold, Markovian laser with:

  • an adiabatically eliminated gain medium;
  • no amplitude–phase coupling;
  • an orthogonal cavity mode;
  • complete inversion in the effective transition;
  • mean intracavity photon number nˉ≫1\bar n\gg1;
  • photon-number loss rate κ\kappa.

In a widely used conventional linear-gain model, the angular-frequency power-spectrum FWHM is

ΔωSQL=κ2nˉ.\Delta\omega_{\mathrm{SQL}} = \frac{\kappa}{2\bar n}.

Because Δω=2πΔν\Delta\omega=2\pi\Delta\nu, this page’s ordinary-frequency FWHM is

ΔνSQL=κ4πnˉ.\Delta\nu_{\mathrm{SQL}} = \frac{\kappa}{4\pi\bar n}.

Equivalently, the phase-diffusion coefficient is

Dϕ=κ4nˉ.D_\phi = \frac{\kappa}{4\bar n}.

The label “standard quantum limit” here names this conventional laser model, not a theorem that no engineered gain process can do better. Quantum gain models with different noise properties require their own derivation.

Let κout\kappa_{\mathrm{out}} be the photon escape rate through the monitored output port. Then

Pout=ℏω0κoutnˉ,P_{\mathrm{out}} = \hbar\omega_0 \kappa_{\mathrm{out}} \bar n,

so

ΔνSQL=ℏω0κκout4πPout.\Delta\nu_{\mathrm{SQL}} = \frac{ \hbar\omega_0 \kappa\kappa_{\mathrm{out}} }{ 4\pi P_{\mathrm{out}} }.

If output coupling dominates all loss, κout≃κ\kappa_{\mathrm{out}}\simeq\kappa. Since the passive cavity FWHM is

Δνc=κ2π,\Delta\nu_c = \frac{\kappa}{2\pi},

the result becomes

ΔνSQL≃ℏω0κ24πPout=πhν0(Δνc)2Pout.\Delta\nu_{\mathrm{SQL}} \simeq \frac{ \hbar\omega_0\kappa^2 }{ 4\pi P_{\mathrm{out}} } = \frac{ \pi h\nu_0(\Delta\nu_c)^2 }{ P_{\mathrm{out}} }.

Every equality in this chain uses a declared definition. Formulas in the literature can differ because the symbol called “cavity linewidth” may be a field HWHM, power HWHM, power FWHM, or angular-frequency width, and because some derivations use total dissipated power rather than one output port.

Real devices can depart from the ideal scale in several conceptually different ways.

Factor or effect Physical origin Typical role
nspn_{\mathrm{sp}} incomplete inversion and excess spontaneous emission multiplies the ideal diffusion scale in common laser models
1+αH21+\alpha_H^2 carrier-induced gain fluctuations coupled to refractive-index fluctuations Henry broadening in semiconductor lasers
KK nonorthogonal left and right cavity modes Petermann excess-noise factor
dispersive pulling gain or intracavity dispersion changes phase response and group delay can alter the conversion from noise to frequency
multiple modes shared gain, mode beating, hopping, and non-normal coupling can invalidate a one-phase description
technical noise pumps, mechanics, electronics, temperature, feedback usually adds a frequency-dependent PSD rather than one constant factor

Under a restricted model, one often encounters the schematic estimate

Δν≃ΔνSQL nsp(1+αH2)K.\Delta\nu \simeq \Delta\nu_{\mathrm{SQL}}\, n_{\mathrm{sp}} (1+\alpha_H^2) K.

This is a correction map, not a universal product law. The factors can be frequency dependent, correlated, or already included in a more microscopic noise model.

In a simple gain medium with upper- and lower-state populations N2N_2 and N1N_1, a conventionally normalized spontaneous-emission factor has the form

nsp∼N2N2−N1,n_{\mathrm{sp}} \sim \frac{N_2}{N_2-N_1},

when emission and absorption cross sections are treated symmetrically. Complete inversion gives nsp→1n_{\mathrm{sp}}\to1. Near transparency, N2−N1N_2-N_1 becomes small and this simplified factor grows. Actual multilevel and semiconductor expressions must use the relevant gain and occupation model.

In a semiconductor, a carrier fluctuation changes both gain and refractive index. The corresponding phase and amplitude quadratures are coupled. Henry’s linewidth-enhancement parameter αH\alpha_H measures that coupling in a declared susceptibility or gain convention. In the standard approximation it broadens the intrinsic line by

1+αH2.1+\alpha_H^2.

The same coupling also produces chirp under current modulation. Semiconductor Lasers Overview defines αH\alpha_H in an intensity-gain convention and develops its device context. The important lesson here is that amplitude restoration does not prevent amplitude noise from being converted into phase noise.

Open, gain-guided, lossy, or otherwise non-Hermitian resonators can have nonorthogonal eigenmodes. For one isolated resonance, the Petermann factor K≥1K\ge1 quantifies enhanced spontaneous-emission noise associated with the left–right mode overlap. It is not simply a statement that two transverse intensity patterns overlap. Near exceptional points or strongly overlapping resonances, a scalar isolated-mode factor may itself become inadequate.

The far-above-threshold phase-diffusion picture can fail:

  • near threshold, where amplitude and phase fluctuations are strongly coupled and the line may be non-Lorentzian;
  • in few-emitter or strong-coupling lasers, where adiabatic elimination and Gaussian noise fail;
  • in bad-cavity or superradiant lasers, where the active medium stores much of the phase;
  • in multimode or mode-hopping operation;
  • under delayed optical feedback or injection locking;
  • when technical noise exceeds the quantum scale.

A smaller fitted line than a naïve Schawlow–Townes estimate is not by itself a violation of quantum mechanics. First verify that the same definitions, ports, loss rates, operating regime, and gain model were used.

For a simple cavity resonance,

νq≃qc2nL,\nu_q \simeq \frac{q c}{2nL},

so small optical-path changes give

δνqνq≃−(δLL+δnn).\frac{\delta\nu_q}{\nu_q} \simeq - \left( \frac{\delta L}{L} + \frac{\delta n}{n} \right).

Mirror displacement, acoustic vibration, thermal expansion, thermorefractive noise, air-pressure changes, coating fluctuations, and pump-induced thermal lensing can therefore become laser frequency noise. A fractional optical path change of 10−1510^{-15} shifts a 300 THz300\ \mathrm{THz} carrier by about 0.3 Hz0.3\ \mathrm{Hz}.

The cold-cavity relation is only a starting point. Gain dispersion and frequency pulling determine how strongly the lasing frequency follows the cavity and material resonances.

Pump noise changes inversion, gain saturation, temperature, and refractive index. It can produce:

  • direct amplitude noise;
  • amplitude-to-phase conversion;
  • carrier-density chirp in semiconductors;
  • thermal frequency drift;
  • relaxation-oscillation peaks;
  • changes of spatial mode or polarization.

Relative intensity noise is commonly defined as

SRIN(f)=SP(1)(f)Pˉ2,S_{\mathrm{RIN}}(f) = \frac{ S_P^{(1)}(f) }{ \bar P^2 },

with units Hz−1\mathrm{Hz^{-1}}. Pure amplitude noise does not automatically broaden the optical carrier, but nonlinear index changes, detector conversion, saturation, and semiconductor αH\alpha_H coupling can turn it into phase noise.

Unintended optical feedback returns a delayed field to the laser. Depending on delay, phase, and strength, it can narrow a line, pull the carrier, create external-cavity modes, produce sidebands, or drive coherence collapse. An isolator reduces this pathway but does not prove that all feedback is absent.

Electronic control loops also reshape noise. Inside the useful loop bandwidth they may suppress free-running noise; near the unity-gain frequency they can create a servo bump; outside the bandwidth the free-running laser remains. A linewidth measured only through the in-loop error signal can be artificially optimistic.

Frequency behavior Typical phase behavior Optical-spectrum consequence
white SνS_\nu Brownian phase Lorentzian component
flicker Sν∝1/fS_\nu\propto1/f increasingly large slow wander observation-time-dependent core
random-walk Sν∝1/f2S_\nu\propto1/f^2 nonstationary frequency drift no unique long-time stationary linewidth
narrow mechanical or electronic peak quasi-periodic phase modulation sidebands or shoulders
discrete mode hop discontinuous carrier change multiple peaks or intermittent spectrum
slow deterministic drift moving center frequency broadened time-averaged trace

Calling all of these effects “linewidth” discards useful diagnostic information.

A direct Fourier spectrum from a record of duration TT has resolution of order 1/T1/T. Noise below that frequency is not physically absent: it can appear as a record-to-record offset or an unresolved drift, depending on detrending and the estimator. Using 1/T1/T as a hard PSD cutoff is a declared approximation. Increasing TT can therefore make a real laser’s fitted line broader by admitting slower noise, even while the instrument’s transform resolution improves.

For white frequency noise, the Lorentzian linewidth approaches a stationary value. For flicker and random-walk noise, the fitted width can depend persistently on:

  • acquisition duration;
  • detrending and windowing;
  • dead time;
  • sweep speed;
  • low-frequency servo cutoff;
  • whether mode hops are excluded;
  • the estimator used.

Frequency-noise PSD, phase-noise PSD, Allan deviation, and optical line shape answer related but different questions. Long-term clock stability should not be inferred from a short-time Lorentzian core alone.

An optical spectrum analyzer, scanning Fabry–Pérot interferometer, or spectrometer can reveal:

  • longitudinal and side modes;
  • broad pedestals;
  • amplified spontaneous emission;
  • mode hops;
  • modulation sidebands.

Its resolution bandwidth and instrument line shape must be narrower and better characterized than the feature being inferred. A trace at the instrument resolution is an upper bound, not a resolved laser linewidth. Deconvolution is model dependent and unstable when the signal-to-noise ratio is poor.

Two fields on a square-law detector produce a beat phase

Φb(t)=2π(ν1−ν2)t+ϕ1(t)−ϕ2(t).\Phi_b(t) = 2\pi(\nu_1-\nu_2)t + \phi_1(t)-\phi_2(t).

For independent lasers, phase variances and frequency-noise PSDs add:

Sν,b(f)=Sν,1(f)+Sν,2(f).S_{\nu,b}(f) = S_{\nu,1}(f) + S_{\nu,2}(f).

The beat spectrum is the convolution of one optical line shape with the frequency-reversed other line shape, because the beat measures a difference frequency. For symmetric lines this reduces to the usual width-addition rules:

  • independent Lorentzian FWHMs add, Δνb=Δν1+Δν2\Delta\nu_b=\Delta\nu_1+\Delta\nu_2;
  • independent Gaussian variances add, σb2=σ12+σ22\sigma_b^2=\sigma_1^2+\sigma_2^2.

Two nominally identical Lorentzians each have half the beat FWHM. Two identical Gaussians each have the beat FWHM divided by 2\sqrt2. Applying the wrong rule is a common factor error.

The measurement is relative. Correlated environmental noise can cancel in the beat, while reference-laser noise can make the device under test appear worse.

In delayed self-heterodyne detection, the laser is split into two arms. One arm is delayed by τd\tau_d, and one is shifted by an acousto-optic or electro-optic frequency fsf_s. Recombination produces a beat near fsf_s with stochastic phase

ϕ(t)−ϕ(t−τd).\phi(t)-\phi(t-\tau_d).

The delay gives the phase-noise transfer magnitude

∣Hd(f)∣2=4sin⁡2(πfτd).|H_d(f)|^2 = 4\sin^2(\pi f\tau_d).

If τd\tau_d greatly exceeds the coherence time, the two samples are nearly uncorrelated. For an ideal Lorentzian laser, the long-delay beat then has

Δνbeat=2Δνlaser.\Delta\nu_{\mathrm{beat}} = 2\Delta\nu_{\mathrm{laser}}.

At finite delay, the spectrum contains correlation fringes and must be fit with the finite-delay model. The delay fiber itself adds thermomechanical and acoustic phase noise, polarization mismatch reduces contrast, and the frequency shifter and detector add their own noise.

An unbalanced interferometer converts small frequency fluctuations into phase and then voltage. Its calibrated response is periodic and frequency-dependent. A stable optical cavity converts detuning into transmitted or reflected intensity or phase. Both methods measure frequency relative to a physical delay or resonance:

δV(f)=Kν(f) δν(f)+δVreadout(f).\delta V(f) = K_\nu(f)\,\delta\nu(f) + \delta V_{\mathrm{readout}}(f).

Recovering SνS_\nu requires the complex discriminator calibration Kν(f)K_\nu(f), detector-noise subtraction, and a check that fluctuations remain inside the linear range. The discriminator’s own length and index noise are indistinguishable from laser frequency noise unless measured independently.

A heterodyne beat can be digitized as

z(t)=I(t)+iQ(t)=A(t)eiΦb(t).z(t) = I(t)+iQ(t) = A(t)e^{i\Phi_b(t)}.

After phase unwrapping,

δνb(t)=12πdΦbdt−ν‾b.\delta\nu_b(t) = \frac{1}{2\pi} \frac{d\Phi_b}{dt}-\overline\nu_b.

Here ν‾b\overline\nu_b is the nominal beat carrier (or zero after demodulation). Differentiation amplifies high-frequency detector noise, while cycle slips produce impulsive errors. A trustworthy analysis reports sampling rate, anti-alias filtering, estimator bandwidth, phase-unwrapping criteria, and the treatment of missing or low-amplitude samples.

Method Strongest use Main limitation
optical spectrum analyzer side modes and broad spectral structure finite resolution and instrument line shape
reference heterodyne relative line shape and phase noise requires a quieter or independently known reference
delayed self-heterodyne no second laser required long stable delay and finite-delay modeling
delay-line discriminator broadband frequency-noise PSD periodic response and delay noise
cavity discriminator high sensitivity near resonance reference-cavity noise and limited linear range
frequency comb or counter traceable frequency comparisons transfer-oscillator and dead-time details

A reproducible linewidth claim should state:

  1. center wavelength or frequency and operating power;
  2. line-shape model and width convention;
  3. measurement method and reference;
  4. resolution bandwidth, sampling rate, and record duration;
  5. Fourier-frequency range for a PSD-derived result;
  6. one-sided or two-sided PSD and hertz or angular-frequency units;
  7. treatment of drift, detrending, sidebands, and mode hops;
  8. uncertainty and instrument-noise floor;
  9. whether the result is in loop or out of loop;
  10. the laser operating point and stabilization state.

Frequency stabilization compares the laser with a discriminator, filters the error, and drives one or more actuators. References include:

  • a passive optical cavity;
  • an atomic or molecular transition;
  • another laser;
  • an optical frequency comb;
  • an interferometric delay.

Actuators can change injection current, cavity length, temperature, intracavity loss, or an external acousto-optic or electro-optic frequency. Fast and slow actuators are often combined because no one actuator has unlimited range and bandwidth.

For scalar open-loop gain G(f)G(f), a schematic uncorrelated noise budget is

Sν,out(f)≃∣11+G(f)∣2Sν,free(f)+∣G(f)1+G(f)∣2[Sν,ref(f)+Sν,sens(f)].S_{\nu,\mathrm{out}}(f) \simeq{} \left| \frac{1}{1+G(f)} \right|^2 S_{\nu,\mathrm{free}}(f) + \left| \frac{G(f)}{1+G(f)} \right|^2 \left[ S_{\nu,\mathrm{ref}}(f) + S_{\nu,\mathrm{sens}}(f) \right].

High gain suppresses free-running laser noise but transfers reference and sensing noise to the output. Delay limits stable bandwidth. Actuator resonances, saturation, and gain peaking can add noise.

A laser tightly locked to one cavity resonance can have excellent short-term coherence while the cavity drifts slowly. An atomic reference can provide better long-term accuracy but lower short-term signal-to-noise ratio. Hybrid systems use a quiet cavity for short times and steer it toward an atomic or comb reference at long times.

The appropriate metric depends on the experiment:

Requirement Useful metric
resolve a narrow transition optical line shape over the interrogation time
preserve a Rabi phase integrated phase noise weighted by the pulse sequence
Ramsey interrogation relative phase variance at the dark time
optical clock fractional frequency instability and systematic accuracy
coherent communication phase-error variance in the receiver bandwidth
frequency-comb transfer residual phase noise of relevant beat locks

Laser Stabilization owns detailed discriminator design, Pound–Drever–Hall modulation, loop shaping, actuator allocation, and stability margins. The central principle here is that a lock does not erase noise; it redistributes noise among the laser, reference, sensor, actuator, and frequency bands.

For a Lorentzian laser with

Δν=10 kHz,\Delta\nu = 10\ \mathrm{kHz},

the field-correlation 1/e1/e time is

τ1/e=1πΔν≃31.8 μs.\tau_{1/e} = \frac{1}{\pi\Delta\nu} \simeq 31.8\ \mu\mathrm{s}.

In vacuum, cτ1/e≃9.55 kmc\tau_{1/e}\simeq9.55\ \mathrm{km}. That is a propagation distance, not a Michelson mirror displacement; the latter is approximately half as large for the same round-trip delay.

Suppose a measured one-sided white plateau is

Sν(1)=100 Hz2/Hz.S_\nu^{(1)} = 100\ \mathrm{Hz^2/Hz}.

The associated Lorentzian FWHM is

Δν=πSν(1)≃314 Hz.\Delta\nu = \pi S_\nu^{(1)} \simeq 314\ \mathrm{Hz}.

Low-frequency noise above that plateau can still make a much broader long-record line.

A target laser is beaten with a reference. The beat is well fit by a Lorentzian of FWHM 180 kHz180\ \mathrm{kHz}, while an independent measurement gives the reference FWHM as 30 kHz30\ \mathrm{kHz}. If the lasers are independent,

Δνtarget=180 kHz−30 kHz=150 kHz.\Delta\nu_{\mathrm{target}} = 180\ \mathrm{kHz} - 30\ \mathrm{kHz} = 150\ \mathrm{kHz}.

If both lines were Gaussian, widths would not subtract linearly:

Δνtarget=(180 kHz)2−(30 kHz)2≃177.5 kHz.\Delta\nu_{\mathrm{target}} = \sqrt{ (180\ \mathrm{kHz})^2 - (30\ \mathrm{kHz})^2 } \simeq 177.5\ \mathrm{kHz}.

Equating passive-cavity and laser linewidth

Section titled “Equating passive-cavity and laser linewidth”

The passive cavity width is set by field storage and loss. The laser output width is set by oscillator phase fluctuations. The former enters the quantum-limited scale but is not the latter.

A Lorentzian FWHM, Gaussian FWHM, β-separation estimate, integrated phase noise, and finite-record spectral width are not interchangeable.

State field versus photon decay, HWHM versus FWHM, angular versus ordinary frequency, and one-sided versus two-sided PSD before substituting numbers.

In diffusion, the absolute phase variance grows without bound. Its increments can be stationary, and those increments determine g(1)g^{(1)}.

White frequency noise supports a stationary Lorentzian width. Flicker, random-walk, resonant, and deterministic components generally require the full line shape, a finite observation time, or an application-specific phase error.

That rule applies to two independent, equal Lorentzians. Equal Gaussian FWHMs combine by 2\sqrt2, and unequal or correlated sources need an explicit model.

An in-loop sensor can be squashed below its own physical noise through correlation with feedback. Use an independent out-of-loop reference to validate the stabilized output.

Coherent modulation redistributes power into discrete spectral components. A width fit that ignores those components can hide the actual modulation physics.

1. Transform the phase-diffusion correlation

Section titled “1. Transform the phase-diffusion correlation”

Starting from

g(1)(τ)=e−i2πν0τe−Dϕ∣τ∣,g^{(1)}(\tau) = e^{-i2\pi\nu_0\tau} e^{-D_\phi|\tau|},

compute the normalized spectrum and identify its HWHM and FWHM in both angular and ordinary frequency.

Solution

Let δν=ν−ν0\delta\nu=\nu-\nu_0. The relevant Fourier transform is

∫−∞∞dτ ei2πδντe−Dϕ∣τ∣=2∫0∞dτ e−Dϕτcos⁡(2πδντ)=2DϕDϕ2+(2πδν)2.\begin{aligned} \int_{-\infty}^{\infty} d\tau\, e^{i2\pi\delta\nu\tau} e^{-D_\phi|\tau|} &= 2 \int_0^\infty d\tau\, e^{-D_\phi\tau} \cos(2\pi\delta\nu\tau) \\ &= \frac{ 2D_\phi }{ D_\phi^2+(2\pi\delta\nu)^2 }. \end{aligned}

After normalization,

L(ν)=1πDϕ/(2π)(ν−ν0)2+[Dϕ/(2π)]2.\mathcal L(\nu) = \frac{1}{\pi} \frac{ D_\phi/(2\pi) }{ (\nu-\nu_0)^2+[D_\phi/(2\pi)]^2 }.

The ordinary-frequency HWHM is Dϕ/(2π)D_\phi/(2\pi) and the FWHM is

Δν=Dϕπ.\Delta\nu = \frac{D_\phi}{\pi}.

The angular-frequency HWHM is DϕD_\phi and the angular-frequency FWHM is 2Dϕ2D_\phi.

A Gaussian optical line has FWHM ΔνG\Delta\nu_{\mathrm G}. Find the 1/e1/e delay of ∣g(1)∣|g^{(1)}| and the two-sided integral ∫∣g(1)(τ)∣2dτ\int|g^{(1)}(\tau)|^2d\tau. Show that they are not equal.

Solution

The rms width is

σν=ΔνG22ln⁡2.\sigma_\nu = \frac{ \Delta\nu_{\mathrm G} }{ 2\sqrt{2\ln2} }.

Since

∣g(1)(τ)∣=e−2π2σν2τ2,|g^{(1)}(\tau)| = e^{-2\pi^2\sigma_\nu^2\tau^2},

the 1/e1/e delay is

τ1/e=12πσν=2ln⁡2πΔνG.\tau_{1/e} = \frac{ 1 }{ \sqrt2\pi\sigma_\nu } = \frac{ 2\sqrt{\ln2} }{ \pi\Delta\nu_{\mathrm G} }.

For the squared magnitude,

∣g(1)(τ)∣2=e−4π2σν2τ2.|g^{(1)}(\tau)|^2 = e^{-4\pi^2\sigma_\nu^2\tau^2}.

Using the Gaussian integral,

∫−∞∞dτ ∣g(1)(τ)∣2=12πσν=2ln⁡2πΔνG.\int_{-\infty}^{\infty} d\tau\, |g^{(1)}(\tau)|^2 = \frac{ 1 }{ 2\sqrt{\pi}\sigma_\nu } = \frac{ \sqrt{2\ln2} }{ \sqrt{\pi}\Delta\nu_{\mathrm G} }.

The two conventions have different numerical coefficients. A coherence time must therefore name its definition.

Let the one-sided frequency-noise PSD be a constant H0H_0. Derive the phase increment variance and Lorentzian FWHM.

Solution

The corresponding two-sided PSD is H0/2H_0/2, so

⟨δν(t)δν(t′)⟩=H02δ(t−t′).\langle \delta\nu(t)\delta\nu(t') \rangle = \frac{H_0}{2} \delta(t-t').

Because

Δϕ=2π∫tt+τδν(t′) dt′,\Delta\phi = 2\pi \int_t^{t+\tau} \delta\nu(t')\,dt',

its variance is

⟨(Δϕ)2⟩=(2π)2H02∣τ∣=2π2H0∣τ∣.\begin{aligned} \langle(\Delta\phi)^2\rangle &= (2\pi)^2 \frac{H_0}{2} |\tau| \\ &= 2\pi^2H_0|\tau|. \end{aligned}

Comparison with 2Dϕ∣τ∣2D_\phi|\tau| gives

Dϕ=π2H0.D_\phi = \pi^2H_0.

Therefore

Δν=Dϕπ=πH0.\Delta\nu = \frac{D_\phi}{\pi} = \pi H_0.

4. Estimate an ideal quantum-limited width

Section titled “4. Estimate an ideal quantum-limited width”

An ideal laser has photon-number decay rate

κ2π=1.0 MHz\frac{\kappa}{2\pi} = 1.0\ \mathrm{MHz}

and mean intracavity photon number nˉ=108\bar n=10^8. Estimate the conventional standard-quantum-limit FWHM.

Solution

Using

ΔνSQL=κ4πnˉ,\Delta\nu_{\mathrm{SQL}} = \frac{\kappa}{4\pi\bar n},

gives

ΔνSQL=2π(106 s−1)4π(108)=5.0×10−3 Hz.\begin{aligned} \Delta\nu_{\mathrm{SQL}} &= \frac{ 2\pi(10^6\ \mathrm{s^{-1}}) }{ 4\pi(10^8) } \\ &= 5.0\times10^{-3}\ \mathrm{Hz}. \end{aligned}

This millihertz-scale number is an ideal model result. It does not include technical noise, incomplete inversion, amplitude–phase coupling, or mode-nonorthogonality.

For the laser in Exercise 4, suppose a restricted model gives nsp=1.2n_{\mathrm{sp}}=1.2, αH=3\alpha_H=3, and K=1.5K=1.5. Estimate the corrected intrinsic width using the schematic multiplicative formula. Explain why the result is not a universal prediction.

Solution

The correction factor is

nsp(1+αH2)K=1.2(1+9)(1.5)=18.n_{\mathrm{sp}} (1+\alpha_H^2) K = 1.2(1+9)(1.5) = 18.

Hence

Δν≃18(5.0 mHz)=90 mHz.\Delta\nu \simeq 18(5.0\ \mathrm{mHz}) = 90\ \mathrm{mHz}.

The product assumes that each factor was defined relative to the same single-mode Markov model and that correlations or frequency dependence can be neglected. A microscopic semiconductor or non-Hermitian multimode model may combine the effects differently. Technical noise must be added through its PSD rather than hidden in this product.

6. Find the coherent peak at finite self-heterodyne delay

Section titled “6. Find the coherent peak at finite self-heterodyne delay”

Let ϕ\phi be a Wiener phase with increment variance 2Dϕ∣τ∣2D_\phi|\tau|, and let ψ(t)=ϕ(t)−ϕ(t−τd)\psi(t)=\phi(t)-\phi(t-\tau_d) be the delayed beat phase. Neglect fiber and detector noise. Find the normalized beat correlation Cb(τ)=⟨e−i[ψ(t+τ)−ψ(t)]⟩C_b(\tau)=\langle e^{-i[\psi(t+\tau)-\psi(t)]}\rangle and the fraction of beat power in its coherent delta peak.

Solution

The two length-τd\tau_d Wiener increments overlap for ∣τ∣<τd|\tau|<\tau_d. Their covariance is 2Dϕ(τd−∣τ∣)2D_\phi(\tau_d-|\tau|) in that interval and zero outside it. Therefore

Var⁡[ψ(t+τ)−ψ(t)]=4Dϕmin⁡(∣τ∣,τd),\operatorname{Var}[\psi(t+\tau)-\psi(t)] =4D_\phi\min(|\tau|,\tau_d),

and

Cb(τ)=e−2Dϕmin⁡(∣τ∣,τd).C_b(\tau)=e^{-2D_\phi\min(|\tau|,\tau_d)}.

Its nonzero long-delay plateau produces a delta peak at the frequency shifter’s beat frequency, with integrated fraction e−2Dϕτd=∣g(1)(τd)∣2e^{-2D_\phi\tau_d}=|g^{(1)}(\tau_d)|^2. The remaining spectrum is a continuous pedestal with finite-delay fringes. For Dϕτd≫1D_\phi\tau_d\gg1, the delta weight is negligible and the pedestal approaches a Lorentzian with twice the original optical FWHM. A finite instrument broadens the delta peak; it does not remove its correlated origin.

7. Design a long-delay self-heterodyne measurement

Section titled “7. Design a long-delay self-heterodyne measurement”

A Lorentzian laser has FWHM 50 kHz50\ \mathrm{kHz}. Estimate its field 1/e1/e coherence time. How much fiber with group velocity 2.0×108 m s−12.0\times10^8\ \mathrm{m\,s^{-1}} gives a delay of five coherence times? What ideal long-delay beat FWHM is expected?

Solution

The coherence time is

τ1/e=1π(50×103 Hz)≃6.37 μs.\tau_{1/e} = \frac{ 1 }{ \pi(50\times10^3\ \mathrm{Hz}) } \simeq 6.37\ \mu\mathrm{s}.

Five coherence times are 31.8 μs31.8\ \mu\mathrm{s}. The required fiber length is

L=vg(5τ1/e)≃(2.0×108)(31.8×10−6)≃6.37 km.L = v_g(5\tau_{1/e}) \simeq (2.0\times10^8) (31.8\times10^{-6}) \simeq 6.37\ \mathrm{km}.

For an ideal Lorentzian in the uncorrelated long-delay limit,

Δνbeat=2Δνlaser=100 kHz.\Delta\nu_{\mathrm{beat}} = 2\Delta\nu_{\mathrm{laser}} = 100\ \mathrm{kHz}.

Five coherence times is a design estimate, not a guarantee. The actual fit should include finite-delay fringes and fiber-added phase noise.

At one Fourier frequency, let G=100G=100 be real. Suppose

Sν,free=106 Hz2/Hz,S_{\nu,\mathrm{free}} = 10^6\ \mathrm{Hz^2/Hz},

while the reference and sensing PSDs are 44 and 9 Hz2/Hz9\ \mathrm{Hz^2/Hz}, respectively. Use the schematic closed-loop equation to estimate the output PSD.

Solution

The free-running contribution is

106∣1+100∣2=10610201≃98.0 Hz2/Hz.\frac{ 10^6 }{ |1+100|^2 } = \frac{10^6}{10201} \simeq 98.0\ \mathrm{Hz^2/Hz}.

The transferred reference-plus-sensor contribution is

10021012(4+9)≃12.7 Hz2/Hz.\frac{ 100^2 }{ 101^2 } (4+9) \simeq 12.7\ \mathrm{Hz^2/Hz}.

Thus

Sν,out≃111 Hz2/Hz.S_{\nu,\mathrm{out}} \simeq 111\ \mathrm{Hz^2/Hz}.

The loop suppresses the large free-running noise, but at high gain it nearly copies the 13 Hz2/Hz13\ \mathrm{Hz^2/Hz} reference and sensing floor. Correlations, complex phase, actuator noise, and plant dynamics would require the full control model.

For a compact separation of passive-cavity width, laser-output linewidth, coherence length, QQ, finesse, Rabi frequency, and saturation language, see Laser Nomenclature.