Squeezed Light Bench

Bow-tie OPO (PPKTP, 532 nm pump → 1064 nm squeezed vacuum) + balanced homodyne detection · 1200 × 900 mm table

In short: Light always carries a tiny, unavoidable flicker called quantum noise. This bench shows how a special crystal inside a ring of mirrors "squeezes" that noise below its normal floor in one property of the light, and how a sensitive detector reads it out. The same trick makes the LIGO gravitational-wave detectors more sensitive.

© 2026 Fabian Magaña-Sandoval, PhD. All rights reserved.
fabianoptics@gmail.com
drag a mirror knob (PM1, PM2, SQ1, SQ2, LM1, LM2, LM3) to turn it · drag elsewhere to orbit · scroll to zoom · pump 532 · LO 1064 · squeezed vacuum

OPO and pump

threshold P_th
—
pump mode overlap
—
pump parameter x = √(P/P_th)
—
OPO linewidth
—
OPO waist in the crystal
—

Pump alignment: knobs on PM1/PM2. Misaligning the pump lowers its overlap with the OPO mode, which acts like lower pump power.

Homodyne detection

no knob selected
homodyne visibility
—
total efficiency η
—
beam w at PD A: squeezed / LO
—
knobs (yaw, pitch µrad)
—

Squeezing

squeezing / anti-squeezing
—
measured now
—
Noise ellipse (purple) vs vacuum (dashed circle); yellow line = LO phase. X₁ is the squeezed quadrature at θ = 0.

Model: below-threshold OPO, V±(f) = 1 ± η·4x / ((1 ∓ x)² + (f/γ)²), with LO phase jitter and detector dark noise. Parameters marked ASSUMED in the README (nonlinearity, crystal temperature bandwidth) are typical values.

Spectrum analyzer, zero span at the analysis frequency: noise vs LO phase (0 dB = shot noise)
squeezed and anti-squeezed quadrature vs frequency; yellow = analysis frequency

Why these knobs

Pump power sets the squeezing strength and bandwidth. OC transmission and loss set the escape efficiency, the main limit. LS position and the LO mirrors set the homodyne overlap; the loss scales as visibility². LO phase picks which quadrature you measure. Phase jitter mixes anti-squeezing in near threshold.

Detailed abstract

Purpose

Answers "what must be on the table to make and measure squeezed light, and how good can it be?" It joins an exact Gaussian-beam layout (bow-tie OPO, pump delivery, homodyne overlap) with the standard below-threshold OPO noise formulas, so alignment, mode matching and losses show up directly as dB of squeezing.

What it can do

Limitations

How realistic is it?

The modeling: core physics and how it is handled

Core physics and how it is handled

1. Bow-tie cavity eigenmode. The stable mode satisfies the round-trip ABCD condition q = (Aq + B)/(Cq + D); stability needs |A + D|/2 < 1. Handled: the engine traces a general astigmatic Gaussian (Arnaud X/Θ form) once round the real ring: crystal faces (PPKTP n_z = 1.8302 at 1064 nm), M4 and M3 concave R = 100 mm, flat M1 and M2, on a 120 × 50 mm rectangle with crossing diagonals. The round-trip matrices give tangential and sagittal waists 53.7 × 55.1 µm, optical round trip 508.4 mm, FSR = c/L = 589.6 MHz, angle of incidence 11.3°.

2. Gaussian-beam propagation through real optics. Handled: exact 3D chief ray plus astigmatic q-parameter through each surface (refraction at spherical and plane surfaces, reflection at tilted mirrors). Includes M3's concave face as a weak lens for the pump, the 6.35 mm N-BK7 output coupler at 11.3° (≈ 0.43 mm lateral shift), the 3 mm BS plate, and plano-convex lenses LP (R = 51.5 mm, f ≈ 100 mm) and LS (R = 103 mm, f ≈ 200 mm). Mirror knob tilts change the chief ray.

3. Mode overlap. Power coupling η = |∫E₁E₂ dA|² / (∫|E₁|² dA · ∫|E₂|² dA). Handled:* numerically on a 72 × 72 grid spanning ±4.5 w plus offset, including offset, tilt (phase ramp) and astigmatism. Used twice: pump vs the OPO's 532 nm mode (w/√2 at the crystal centre, same Rayleigh range) and LO vs squeezed beam at PD A. Defaults: pump 99.1 %, homodyne 99.4 % (visibility 99.7 %).

4. Threshold and phase matching. P_th = (T + L)² / (4 E_NL · s²), with s = sinc(2.7832 · ΔT / 1.0 °C), so the conversion efficiency s² has a 1.0 °C FWHM. Handled: E_NL = 0.017 W⁻¹ (ASSUMED). Pump parameter x = √(η_pump · P / P_th). Default P_th = 153 mW, x = 0.805 at 100 mW.

5. Linewidth. FWHM = FSR · (T + L) / (2π), half-width γ = FWHM/2. Handled: closed form (low-loss ring approximation); default 9.57 MHz FWHM, γ = 4.79 MHz.

6. Squeezing spectra. V∓(f) = 1 ∓ η · 4x / ((1 ± x)² + (f/γ)²), relative to shot noise (1 = vacuum). Handled: closed form, x capped at 0.999. Total efficiency η = η_esc · (1 − L_prop) · QE · η_hd, with η_esc = T/(T + L). Default η = 0.980 × 0.97 × 0.95 × 0.994 = 0.898.

7. LO phase, jitter and dark noise. V(θ) = V₋⟨cos²θ⟩ + V₊⟨sin²θ⟩ with ⟨cos²θ⟩ = (1 + cos 2θ · e^(−2σ²))/2 for Gaussian jitter σ. Dark noise: V_meas = (V + N_e)/(1 + N_e), N_e = 10^(−clearance/10). Handled: closed form; σ = 10 mrad, clearance 15 dB. Result at 1 MHz: −8.12 dB squeezing, +15.48 dB anti-squeezing.

8. Phase-space picture. Handled: the noise ellipse is drawn with semi-axes √V₊ and √V₋ against the vacuum circle; the LO phase is a line through it (display only).

How the simulation is built

How to recreate it

  1. Build the ring: four mirrors at the rectangle corners, two R = 100 mm curved, a 10 mm crystal of index 1.8302 between them; compute the round-trip ABCD (tangential and sagittal, with effective R·cos θ and R/cos θ at 11.3°) and solve for q.
  2. Compute L_opt, FSR and FWHM = FSR(T + L)/2π.
  3. Propagate the pump (532 nm) and LO (1064 nm) Gaussian beams through lenses and mirrors to the crystal and PD A; compute overlaps (analytic for aligned Gaussians, or numerically with offsets/tilts).
  4. Compute P_th, x, η and the spectra V∓(f) from items 4–6.
  5. Add LO-phase averaging and dark noise (item 7); plot dB = 10 log₁₀ V vs θ and vs f.
  6. Check defaults: P_th ≈ 153 mW, −8.1 / +15.5 dB at 1 MHz.
  7. To rerun: node sq_design.js (optional, re-solves positions), then python build.py and open squeezed_light_bench.html.

Key parameters

Parameter Value Source
Wavelengths 1064 nm (squeezed, LO), 532 nm (pump) design choice
Cavity geometry 120 × 50 mm bow-tie, R = 100 mm curved mirrors Chua et al. 2011 (style); design choice
Optical round trip / FSR 508.4 mm / 589.6 MHz computed
Waist in crystal 53.7 × 55.1 µm computed
PPKTP length, n_z 10 mm; 1.8302 (1064), 1.8887 (532) datasheet (approximate)
E_NL 0.017 W⁻¹ ASSUMED
Phase-matching width 1.0 °C FWHM (of sinc²) ASSUMED
OC transmission T 10 % design choice
Intracavity loss L 0.2 % ASSUMED
Pump power 100 mW design choice
Detector QE 95 % ASSUMED
Propagation loss 3 % ASSUMED
LO phase jitter 10 mrad rms ASSUMED
Dark-noise clearance 15 dB ASSUMED
Analysis frequency 1 MHz design choice
The shared 3D bench kit (how the table, mounts, beams and parts list are drawn)

What the shared bench engine is

Every interactive page loads three shared files. engine.js is the optics math, kit.js is the three.js scene, stand-in hardware, beam drawing and parts table, and clamps.js places the posts and forks. Units are mm in a table frame (x, y, z up). Each experiment's physics lives in its own model file. The engine only does the optics that a model explicitly asks for.

How beams are traced

  • Chief ray: exact 3D tracing at infinite planes, spheres and cylinders. Refraction uses vector Snell's law and fails on total internal reflection. Mirrors use d' = d − 2(d·N)N. A 'stop' surface only intersects. Air is n = 1.000278 and N-BK7 comes from its Sellmeier formula. Each wavelength gets its own instance via make(λ).
  • Gaussian beam: a general astigmatic beam in Arnaud form. Complex 2 × 2 matrices X and Θ give W = Θ·X⁻¹ = n·Q⁻¹. Free space is X += (L/n)·Θ, the matrix form of the q-parameter (ABCD) law. At a surface, X' = A·X and Θ' = C·X + D·Θ. A and D come from projecting the beam basis onto the surface tangents, and C = −Δ·P₂⁻ᵀ·K·P₁⁻¹ with Δ = n₂(d₂·N) − n₁(d₁·N) and K the curvature tensor. Tilted and cylindrical surfaces therefore produce correct astigmatism. Gouy phase is accumulated per axis.
  • Sampling: samples are taken about every 6 mm and refined wherever the radius changes by more than 12 % between samples. The 1/e² ellipse at each sample comes from Im(W).
  • Splits: there is no automatic splitting. A model clones the beam (cloneB) at a beamsplitter and traces each branch with its own surface list. Surfaces have no apertures, so the model must handle clipping.
  • Extras: Hermite and Laguerre polynomials, and an FFT angular-spectrum propagator for mode-shape pages.

Where readings come from

The kit's photodiode and camera objects are bodies on posts and compute nothing. Photodiode voltages, quadrant read-outs, spectra and noise are all computed in each model file. These calculations are often analytic and use no engine at all, as in the Faraday, EOM and direction-finder pages. Labels are CSS2D text set by the app.

Parts table

Each stand-in part and post stack registers a row: reference, category, part number, description. A stack registers a TR post, a PH40 or PH50 holder (chosen for the beam height), a BE1/M base, and a CF125C/M fork unless the post sits directly on a hole. partsTable() groups and counts identical rows, sorts them by category, and writes the HTML table and the CSV.

Clamps and obstacles

Each fork's M6 screw must land on a 25 mm-grid hole 26.4–49.6 mm from the post. Forks perpendicular to the beam are preferred. A sliding part keeps its hole while that hole still fits, then re-clamps. Forks may not overlap other forks, pedestals or registered footprints ("obstacles"). Conflicts turn red, and a guard stops posts from colliding. All of this is mechanical layout only; it has no optical effect.

Purely visual

  • Beam tubes draw the ellipse samples, optionally magnified, with opacity on a log scale of power.
  • Some pages draw hand-made tube samples instead of traced ones.
  • Meshes are parametric stand-ins, not CAD.

Parts list (generated from the 3D model)

Thorlabs-style suggestions; verify part numbers before ordering. "custom" means a coating or crystal vendor.

CSV (select and copy)

References and related literature

A standard continuous-wave squeezer: a below-threshold degenerate OPO in a bow-tie cavity, read out by balanced homodyne detection. The bow-tie geometry follows the gravitational-wave squeezer designs (Chua et al.); the noise spectra follow Collett & Gardiner.

  1. based on: bow-tie OPOS. S. Y. Chua, M. S. Stefszky, C. M. Mow-Lowry, B. C. Buchler, S. Dwyer, D. A. Shaddock, P. K. Lam and D. E. McClelland, "Backscatter tolerant squeezed light source for advanced gravitational-wave detectors," Opt. Lett. 36, 4680 (2011).
  2. classic experimentL.-A. Wu, H. J. Kimble, J. L. Hall and H. Wu, "Generation of squeezed states by parametric down conversion," Phys. Rev. Lett. 57, 2520 (1986).
  3. M. J. Collett and C. W. Gardiner, "Squeezing of intracavity and traveling-wave light fields produced in parametric amplification," Phys. Rev. A 30, 1386 (1984): the OPO spectra V±(f).
  4. H. Vahlbruch, M. Mehmet, K. Danzmann and R. Schnabel, "Detection of 15 dB squeezed states of light and their application for the absolute calibration of photoelectric quantum efficiency," Phys. Rev. Lett. 117, 110801 (2016).
  5. R. Schnabel, "Squeezed states of light and their applications in laser interferometers," Phys. Rep. 684, 1 (2017).
  6. LIGO Scientific Collaboration, "A gravitational wave observatory operating beyond the quantum shot-noise limit," Nat. Phys. 7, 962 (2011).
  7. M. Tse et al., "Quantum-enhanced Advanced LIGO detectors in the era of gravitational-wave astronomy," Phys. Rev. Lett. 123, 231107 (2019).
  8. G. D. Boyd and D. A. Kleinman, "Parametric interaction of focused Gaussian light beams," J. Appl. Phys. 39, 3597 (1968): optimum focusing in the crystal.
  9. H.-A. Bachor and T. C. Ralph, A Guide to Experiments in Quantum Optics, 2nd ed. (Wiley-VCH, 2004): homodyne detection, OPO practice.
  10. H. Kogelnik and T. Li, "Laser beams and resonators," Appl. Opt. 5, 1550 (1966).
  11. J. A. Arnaud, "Nonorthogonal optical waveguides and resonators," Bell Syst. Tech. J. 49, 2311 (1970): general astigmatic Gaussian beams (the X/Θ form used by the optics engine).