Resonator Stability and the g-Parameter Diagram

A laser will not oscillate just because it has mirrors and gain between them. The feedback has to close on itself: some transverse field pattern must leave a mirror, cross the cavity, reflect, return, and arrive as a faithful scaled copy of what set out — a self-reproducing mode. For the cleanest such pattern that copy is a Gaussian beam, and whether one exists with a real, finite spot size is settled entirely by geometry: the mirror radii R1R_1 and R2R_2 and their separation LL. Move the mirrors a centimetre and a cavity that held a tidy sub-millimetre spot can stop supporting any confined mode at all. This synthesis leans on two results it does not re-derive. From aq-optics-abcd-matrix-beam-propagation it takes the ray-transfer matrix, which compresses a whole round trip of gaps and mirror reflections into one 2×22\times2 matrix of unit determinant. From aq-optics-gaussian-beams it takes the complex beam parameter qq, the single number carrying a Gaussian's spot size and wavefront curvature together. Demanding that qq reproduce itself after one pass turns "does a mode exist?" into "does this matrix rotate or magnify?" — and the answer collapses into two dimensionless numbers, gi=1−L/Rig_i = 1 - L/R_i, and one inequality, 0≤g1g2≤10 \le g_1 g_2 \le 1. Plotted in the g1g_1–g2g_2 plane, that inequality carves out the stable region and gives every resonator a home address. It is the single most useful design chart in laser physics, telling you not just whether a cavity works but how close it sits to the cliff edge of failure.

Key takeaways
  • A two-mirror cavity confines a stable Gaussian mode only when 0≤g1g2≤10 \le g_1 g_2 \le 1, with the g-parameter gi=1−L/Rig_i = 1 - L/R_i folding each mirror's curvature and the spacing into one number.
  • The criterion is the ray-and-Gaussian statement that the round-trip 2×22\times2 matrix must rotate rather than magnify: its half-trace (A+D)/2=2g1g2−1(A+D)/2 = 2g_1g_2 - 1 must have magnitude ≤1\le 1.
  • In the g1g_1–g2g_2 plane the stable set is the two leaf-shaped regions hugging the origin between the coordinate axes and the hyperbola branches g1g2=1g_1 g_2 = 1.
  • Named cavities are landmarks: confocal sits at the origin (0,0)(0,0), plane-parallel at the corner (1,1)(1,1), concentric at (−1,−1)(-1,-1), hemispherical on an axis at (1,0)(1,0).
  • Being inside the region is not the same as being good: near a boundary the mode radius blows up as 1/1−g1g21/\sqrt{1-g_1g_2} and alignment tolerance vanishes; confocal is the most forgiving point.
  • Instability is not failure — deliberately unstable resonators (g1g2>1g_1g_2 > 1 or <0< 0) fill a large gain volume and are the standard choice for high-gain, high-power lasers.
  • Key idea: fold R1R_1, R2R_2, LL into g1=1−L/R1g_1 = 1-L/R_1 and g2=1−L/R2g_2 = 1-L/R_2; the cavity holds a confined Gaussian mode if and only if 0≤g1g2≤10 \le g_1 g_2 \le 1, the stable region of the g-parameter plane.
The core insight

Fold the three geometric quantities R1R_1, R2R_2, LL into g1=1−L/R1g_1 = 1-L/R_1 and g2=1−L/R2g_2 = 1-L/R_2: the cavity supports a confined Gaussian mode if and only if 0≤g1g2≤10 \le g_1 g_2 \le 1, the stable region of the g-parameter plane.

Where this lives today

The g-diagram is the first sketch on the whiteboard whenever someone designs a laser. A titanium-sapphire mode-locked oscillator, the workhorse of ultrafast science, is a folded cavity whose curved mirrors place a tight waist inside the crystal for strong Kerr-lensing while keeping the resonator just inside the stable region — and Kerr-lens mode-locking works by exploiting how the operating point shifts on the diagram between the continuous-wave and pulsed states, so a soft aperture favours the pulse. High-power solid-state lasers live at the opposite extreme: their gain rods develop a thermal lens that changes RR, hence gg, as pump power rises, and the design brief is to keep the operating point from wandering across a boundary over the full power range — a trajectory drawn directly on the g-plane. At the precision frontier the same criterion governs the kilometre-scale Fabry–Pérot arms of gravitational-wave detectors such as LIGO and Virgo, where mirror curvatures are engineered to hundreds of parts per million so the resonator sits at a safe point with a large, well-defined spot on each test mass. And where a large mode beats a clean Gaussian — high-energy pulsed CO2_2 and solid-state systems — engineers step outside the stable region on purpose, using an unstable resonator whose magnifying round trip fills the whole gain aperture. In every case the first question is the one this article answers: where does the cavity sit in the g1g_1–g2g_2 plane, and how much margin does it have?

The physical picture

Treat one round trip as a single lens-like operation that keeps re-imaging the beam onto itself. A stable resonator is one in which a ray launched near the axis bounces back and forth forever, oscillating about the centre but never running away; an unstable one magnifies the ray's transverse excursion on every pass until it walks off the edge of a mirror. The g-parameters measure this focusing balance. Too little net focusing — mirrors too flat or too far apart — and the beam diffracts wider each pass and leaks out the sides. Too much — mirrors so strongly curved that the geometry drives past concentric — and the beam over-focuses and diverges. Between the two lies a corridor where the round-trip focusing is just restoring enough to hold a bounded, self-reproducing beam. A homely analogy is a marble rolling between two curved bowls: if the combined curvature nudges it back toward the centre each pass its path stays bounded, but if the geometry defocuses on net, any small nudge grows without limit. The product g1g2g_1 g_2 is the single parameter deciding between oscillatory (bounded) and hyperbolic (runaway) ray dynamics — and the same product decides whether a Gaussian mode has a real spot size or an imaginary one.

Figure
A ray bouncing between two cavity mirrors stays bounded when 0 ≤ g₁g₂ ≤ 1 (left) but its transverse excursion grows every pass and walks off the mirror edge when g₁g₂ lies outside that range (right)
Stable: bounded ray · 0 ≤ g₁g₂ ≤ 1 Unstable: ray walks off · g₁g₂ outside 0…1 mirror 1 mirror 2 excursion stays within the mirrors → mode survives excursion grows each pass → beam misses the mirror

The self-consistency picture is complementary and actually fixes the mode. Insist that the Gaussian's parameter qq come back to itself after a full round trip: that is an eigenvalue problem for the round-trip matrix, and a confined mode with a real waist exists only when the eigenvalues are complex conjugates on the unit circle, e±iθe^{\pm i\theta}, rather than a real magnifying pair. That condition is ∣A+D∣≤2|A+D| \le 2, nothing other than 0≤g1g2≤10 \le g_1 g_2 \le 1. Distance to the boundary carries real physics: near g1g2→1g_1g_2 \to 1 or →0\to 0 the spot swells to overfill the mirrors, diffraction loss spikes, and the slightest tilt throws the mode off. The chart below turns three lengths into one point and shows at a glance how much margin that point has.

Figure
The stability diagram in the g₁–g₂ plane: stable when 0 ≤ g₁g₂ ≤ 1 (green), bounded by the axes and the hyperbola g₁g₂ = 1 — confocal (0,0), plane-parallel (1,1), concentric (−1,−1), hemispherical (1,0)
Every two-mirror cavity is one point (g₁, g₂); green is stable, 0 ≤ g₁g₂ ≤ 1 g₁ g₂ 1 −1 1 −1 g₁g₂ = 1 STABLE confocal (0,0) plane-parallel (1,1) concentric (−1,−1) hemispherical (1,0) (−0.5,−0.5), g₁g₂ = 0.25 HeNe (1, 0.17)

The core equations

The whole apparatus rests on one definition. Each mirror gets a dimensionless g-parameter measuring how strongly its curvature focuses relative to the cavity length: a flat mirror gives g=1g = 1, one with radius equal to the spacing gives g=0g = 0. The sign convention is fixed here: R>0R > 0 for a mirror concave toward the cavity.

gi=1−LRi,i=1,2(1)g_i = 1 - \frac{L}{R_i}, \qquad i = 1,2 \qquad (1)

where: gig_i — dimensionless stability parameter of mirror ii; LL — on-axis mirror separation; RiR_i — radius of curvature of mirror ii, positive if concave toward the cavity. Physically: as LL grows toward RiR_i the focusing uses up the cavity length and gi→0g_i \to 0; push past RiR_i and gig_i goes negative. Everything geometric about stability funnels into the product of the two g-parameters — the single discriminant. Not g1g_1 and g2g_2 separately, only their product decides whether a confined mode exists. It is the master result, the boundary of the diagram.

0≤g1g2≤1(2)0 \le g_1 g_2 \le 1 \qquad (2)

where: g1g2g_1 g_2 — product of the two g-parameters; a value in [0,1][0,1] is necessary and sufficient for a confined Gaussian mode. At g1g2=0g_1g_2 = 0 or 11 the cavity is marginally stable; values outside make it unstable. The bridge from the matrix to this inequality is the round-trip half-trace. Any lossless round-trip matrix is unimodular (AD−BC=1AD - BC = 1), so its eigenvalues multiply to one: either a conjugate pair on the unit circle (stable) or a reciprocal real pair (unstable), and the half-trace decides which. Multiplying out the four matrices of a two-mirror trip gives the compact result below.

A+D2=2 g1g2−1  ≡  cos⁡θ(3)\frac{A+D}{2} = 2\,g_1 g_2 - 1 \;\equiv\; \cos\theta \qquad (3)

where: A,DA, D — diagonal elements of the round-trip matrix; θ\theta — transverse round-trip (Gouy) phase per pass. Stability is ∣cos⁡θ∣≤1|\cos\theta| \le 1, i.e. −1≤2g1g2−1≤1-1 \le 2g_1g_2 - 1 \le 1, exactly equation (2); the identity θ=2arccos⁡g1g2\theta = 2\arccos\sqrt{g_1 g_2} ties the geometry to the transverse-mode spacing. To connect the matrix to an actual beam we need the object a round trip acts on. A Gaussian beam is fully described at any plane by one complex number, the beam parameter qq, whose real and imaginary parts package the wavefront curvature and the spot size.

1q=1R−i λπw2(4)\frac{1}{q} = \frac{1}{R} - i\,\frac{\lambda}{\pi w^2} \qquad (4)

where: qq — complex beam parameter; RR — wavefront radius of curvature; ww — beam radius (1/e21/e^2-intensity spot size); λ\lambda — wavelength in the medium. The imaginary part makes the spot real and finite; a mode with Im(1/q)=0\mathrm{Im}(1/q) = 0 is a plane wave, not a confined beam. A round trip transforms qq by the fractional-linear map governing any ABCD system. A self-reproducing mode is a fixed point of this map — the beam that comes back identical to itself.

q2=A q1+BC q1+D(5)q_2 = \frac{A\,q_1 + B}{C\,q_1 + D} \qquad (5)

where: A,B,C,DA, B, C, D — round-trip matrix elements; q1,q2q_1, q_2 — beam parameters before and after the trip. Setting q2=q1q_2 = q_1 gives a quadratic with complex roots — a physical, finite-waist mode — precisely when ∣A+D∣≤2|A+D| \le 2, recovering equation (2). Solving that fixed-point condition gives the mode radii on the two mirrors — how the diagram earns its keep, turning a (g1,g2)(g_1, g_2) point into the actual spot a designer must fit inside a mirror and match an external beam to.

w12=Lλπg2g1(1−g1g2),w22=Lλπg1g2(1−g1g2)(6)w_1^2 = \frac{L\lambda}{\pi}\sqrt{\frac{g_2}{g_1(1-g_1 g_2)}}, \qquad w_2^2 = \frac{L\lambda}{\pi}\sqrt{\frac{g_1}{g_2(1-g_1 g_2)}} \qquad (6)

where: w1,w2w_1, w_2 — mode radii at mirrors 1 and 2; L,λL, \lambda — spacing and wavelength. The factor 1/1−g1g21/\sqrt{1-g_1g_2} is the tell-tale of the outer boundary: as g1g2→1g_1g_2 \to 1 both spots diverge and the mode overfills the mirrors. The internal beam waist rounds out the toolkit, collapsing to a memorable form for the symmetric confocal cavity where the Rayleigh range equals half the cavity length.

w02=Lλπg1g2(1−g1g2)(g1+g2−2g1g2)2,w0,conf2=Lλ2π(7)w_0^2 = \frac{L\lambda}{\pi}\sqrt{\frac{g_1 g_2 (1 - g_1 g_2)}{(g_1 + g_2 - 2 g_1 g_2)^2}}, \qquad w_{0,\text{conf}}^2 = \frac{L\lambda}{2\pi} \qquad (7)

where: w0w_0 — intracavity waist radius; w0,confw_{0,\text{conf}} — waist of the symmetric confocal resonator (g1=g2=0g_1 = g_2 = 0, R1=R2=LR_1=R_2=L); L,λL, \lambda — length and wavelength. The confocal waist is the smallest mode any symmetric stable cavity of length LL can hold — the natural mode-matching reference.

Derivation AStability from the round-trip ray matrix

The goal is to derive 0≤g1g2≤10 \le g_1 g_2 \le 1 from nothing but the demand that a ray launched near the axis stay bounded after arbitrarily many round trips. We model one trip as an ordered product of ray-transfer matrices, use its unimodularity to classify the eigenvalues, and show that boundedness is the single statement that the half-trace lies between −1-1 and +1+1. Multiplying the four matrices then delivers a strikingly clean half-trace, and the inequality falls straight out.

Step 1. Build the round-trip matrix as gaps and reflections in sequence, M=MR1MLMR2MLM = M_{R_1} M_L M_{R_2} M_L, with ML=(1L01)M_L = \left(\begin{smallmatrix} 1 & L \\ 0 & 1 \end{smallmatrix}\right) for a gap of length LL and MR=(10−2/R1)M_R = \left(\begin{smallmatrix} 1 & 0 \\ -2/R & 1 \end{smallmatrix}\right) for reflection off a mirror of radius RR. Why: a mirror of radius RR focuses like a lens of focal length R/2R/2, hence the entry −2/R-2/R; stringing the operations in the order the ray meets them is the whole content of one round trip.

Step 2. Each factor has determinant 11, so det⁡M=1\det M = 1 and the eigenvalues satisfy μ+μ−=1\mu_+\mu_- = 1 with μ±=e±iθ\mu_\pm = e^{\pm i\theta}, cos⁡θ=12(A+D)\cos\theta = \tfrac12(A+D). A ray after NN trips is MN(y0,y0′)TM^N(y_0, y_0')^{\mathrm T}; it stays bounded iff θ\theta is real, i.e. ∣12(A+D)∣≤1|\tfrac12(A+D)| \le 1. Why: a unimodular 2×22\times2 matrix either rotates (complex eigenvalues on the unit circle, bounded orbits) or magnifies (real reciprocal eigenvalues, runaway) — the half-trace alone decides.

Step 3. Multiply the four matrices and simplify using gi=1−L/Rig_i = 1 - L/R_i. The off-diagonal clutter cancels and the half-trace reduces to 12(A+D)=2(1−L/R1)(1−L/R2)−1=2g1g2−1\tfrac12(A+D) = 2(1-L/R_1)(1-L/R_2) - 1 = 2g_1 g_2 - 1. Why: this is the algebraic miracle of the topic — three lengths and two curvatures collapse into a single product, which is why a two-dimensional chart suffices to classify every two-mirror cavity.

Step 4. Impose boundedness: −1≤2g1g2−1≤1-1 \le 2g_1 g_2 - 1 \le 1, so g1g2≥0g_1 g_2 \ge 0 and g1g2≤1g_1 g_2 \le 1, together 0≤g1g2≤10 \le g_1 g_2 \le 1. Why: the two halves are physically distinct — the lower bound rules out defocusing runaway (product negative), the upper bound rules out over-focusing runaway (product past unity) — and only their intersection is a genuine resonator.

Go deeper · Derivation B · 1 minStability from Gaussian self-consistency
Derivation BStability from Gaussian self-consistency

The ray argument proves stability but never mentions the beam. The goal here is to reach the same inequality from the Gaussian side, demanding that qq return to itself after one round trip with a real, finite spot size. This route not only reconfirms 0≤g1g2≤10 \le g_1 g_2 \le 1 but hands us the mode: the qq that solves the fixed-point equation is the physical Gaussian the cavity supports, its imaginary part the waist.

Step 1. Write the self-consistency condition q=(Aq+B)/(Cq+D)q = (Aq+B)/(Cq+D), the statement that the round-trip map leaves the mode unchanged. Cross-multiplying gives B(1/q)2+(D−A)(1/q)−C=0B(1/q)^2 + (D-A)(1/q) - C = 0. Why: a resonator mode is by definition the field that reproduces after one trip, and for a Gaussian that field is entirely captured by qq — so the mode is just a fixed point of a fractional-linear map.

Step 2. Solve for 1/q1/q:   1/q=D−A2B±12B(A−D)2+4BC\;1/q = \dfrac{D-A}{2B} \pm \dfrac{1}{2B}\sqrt{(A-D)^2 + 4BC}. Using unimodularity AD−BC=1AD - BC = 1, the discriminant becomes (A+D)2−4(A+D)^2 - 4. Why: the determinant constraint turns a messy discriminant into the clean (A+D)2−4(A+D)^2 - 4, so its sign is governed by the same half-trace that controlled the ray picture.

Step 3. Demand a physical mode. The spot is real and finite only if 1/q1/q has a nonzero imaginary part, requiring a negative discriminant: (A+D)2−4<0(A+D)^2 - 4 < 0, i.e. ∣A+D∣≤2|A+D| \le 2. Why: Im(1/q)=−λ/πw2\mathrm{Im}(1/q) = -\lambda/\pi w^2 must be nonzero for a finite ww; a real 1/q1/q is a plane wave of infinite spot, so "reality of the waist" is literally the stability boundary.

Step 4. Substitute 12(A+D)=2g1g2−1\tfrac12(A+D) = 2g_1g_2 - 1: ∣2g1g2−1∣≤1⇔0≤g1g2≤1|2g_1g_2 - 1| \le 1 \Leftrightarrow 0 \le g_1 g_2 \le 1. Why: ray and Gaussian pictures land on the identical inequality — "rays stay bounded" and "a Gaussian mode has a real waist" are two faces of one geometric fact.

Go deeper · Derivation C · 1 minThe confocal cavity waist
Derivation CThe confocal cavity waist

The goal is to specialize the mode formulas to the symmetric confocal resonator — the cavity at the origin of the diagram — and extract its waist in closed form. Confocal geometry sets R1=R2=LR_1 = R_2 = L, hence g1=g2=0g_1 = g_2 = 0, and the fastest route sidesteps the 0/00/0 the general waist formula throws up by using the defining fact that each mirror sits one Rayleigh range from the central waist. This yields the most-quoted mode size in resonator design.

Step 1. Place the waist at the cavity centre by symmetry; the confocal condition puts the mirrors at ±zR\pm z_R, so b=2zR=Lb = 2z_R = L and zR=L/2z_R = L/2. Why: "confocal" means the two mirror focal points coincide at the centre, forcing the Rayleigh range — a Gaussian's natural depth of focus — to equal half the cavity length, which fixes zRz_R.

Step 2. Invert zR=πw02/λz_R = \pi w_0^2/\lambda to get w02=λzR/π=Lλ/2πw_0^2 = \lambda z_R/\pi = L\lambda/2\pi, so w0=Lλ/2πw_0 = \sqrt{L\lambda/2\pi}. Why: zRz_R couples waist size to wavelength, so once geometry pins zRz_R the waist follows with no further physics.

Step 3. Read off the spot on each mirror. At z=±zRz = \pm z_R a Gaussian has expanded by 2\sqrt 2, so wmirror=2 w0w_{\text{mirror}} = \sqrt2\,w_0 and wmirror2=Lλ/πw_{\text{mirror}}^2 = L\lambda/\pi. Why: this is the smallest mirror spot any symmetric stable cavity of length LL can carry, which is why confocal is the deepest, most alignment-tolerant point on the diagram and the default mode-matching target.

Limiting cases

The named resonators are landmarks — (g1,g2)(g_1, g_2) points where the geometry has a familiar name and the mode behaviour is worth memorizing.

CaseResultMeaning
Plane-parallel (Fabry–Pérot)R1=R2=∞⇒g1=g2=1R_1=R_2=\infty \Rightarrow g_1=g_2=1Corner (1,1)(1,1) on the boundary g1g2=1g_1g_2=1; marginally stable, extreme tilt sensitivity, largest mode volume
Symmetric confocalR1=R2=L⇒g1=g2=0R_1=R_2=L \Rightarrow g_1=g_2=0Origin of the diagram; deepest inside the stable region, most forgiving alignment, smallest symmetric mode
Symmetric concentricR1=R2=L/2⇒g1=g2=−1R_1=R_2=L/2 \Rightarrow g_1=g_2=-1Corner (−1,−1)(-1,-1) on g1g2=1g_1g_2=1; foci meet at centre, tiny waist but again only marginally stable
HemisphericalR1=∞, R2=L⇒g1=1, g2=0R_1=\infty,\ R_2=L \Rightarrow g_1=1,\ g_2=0On the g2=0g_2=0 axis; one flat and one curved mirror, a common and tunable near-stable design
Boundary (g1g2→0g_1g_2 \to 0 or 11)mode radius ∝1/1−g1g2\propto 1/\sqrt{1-g_1g_2} divergesMarginal stability; spot overfills the mirrors, diffraction loss and alignment sensitivity blow up
Unstableg1g2>1g_1g_2>1 or g1g2<0g_1g_2<0No confined Gaussian mode; used on purpose for high-gain, high-power lasers with geometric output coupling

Worked examples

1 · Mode sizes of a hemispherical HeNe cavity (confidence builder)

A helium–neon laser at λ=632.8\lambda = 632.8 nm uses a flat output coupler (R1=∞R_1 = \infty) and a concave back mirror R2=0.30R_2 = 0.30 m, separated by L=0.25L = 0.25 m. Is it stable, and how big is the mode on each mirror?

Step 1. Compute the g-parameters. The flat mirror gives g1=1−0=1g_1 = 1 - 0 = 1; the curved mirror gives g2=1−0.25/0.30=0.1667g_2 = 1 - 0.25/0.30 = 0.1667. The product g1g2=0.1667g_1 g_2 = 0.1667 sits inside [0,1][0,1], so the cavity is stable — it plots on the g1=1g_1 = 1 line just inside the boundary, the red HeNe point above.

Step 2. Because g1=1g_1 = 1, the wavefront is flat there, so the waist sits on the flat mirror and w0=w1w_0 = w_1. With Lλ/π=(0.25)(632.8×10−9)/π=5.036×10−8 m2L\lambda/\pi = (0.25)(632.8\times10^{-9})/\pi = 5.036\times10^{-8}\ \mathrm{m}^2, equation (6) gives 0.1667/0.8333=0.4472\sqrt{0.1667/0.8333} = 0.4472, so w12=5.036×10−8×0.4472=2.252×10−8 m2w_1^2 = 5.036\times10^{-8}\times0.4472 = 2.252\times10^{-8}\ \mathrm{m}^2 and w1=w0=0.150 mmw_1 = w_0 = \mathbf{0.150\ mm}.

Step 3. The spot on the curved mirror uses the companion formula: g1/[g2(1−g1g2)]=1/(0.1667×0.8333)=7.20=2.683\sqrt{g_1/[g_2(1-g_1g_2)]} = \sqrt{1/(0.1667\times0.8333)} = \sqrt{7.20} = 2.683, so w22=5.036×10−8×2.683=1.351×10−7 m2w_2^2 = 5.036\times10^{-8}\times2.683 = 1.351\times10^{-7}\ \mathrm{m}^2 and w2=0.368 mmw_2 = \mathbf{0.368\ mm}. The mode fans from a 0.150.15 mm waist on the flat mirror to a 0.370.37 mm spot on the curved one. Push LL toward 0.300.30 m and g2→0g_2 \to 0: the cavity nears the hemispherical boundary, the mode swells, and alignment gets touchy — the tuning knob every HeNe builder knows. Stable, with w0=0.150w_0 = 0.150 mm at the flat mirror and 0.3680.368 mm at the curved mirror.

2 · An asymmetric cavity that fails the sign test

A Nd:YAG design (λ=1064\lambda = 1064 nm) proposes mirrors R1=1.0R_1 = 1.0 m and R2=0.5R_2 = 0.5 m at spacing L=0.8L = 0.8 m. Will it lase on a clean mode? If not, fix it.

Step 1. Evaluate the g-parameters. g1=1−0.8/1.0=0.20g_1 = 1 - 0.8/1.0 = 0.20; g2=1−0.8/0.5=−0.60g_2 = 1 - 0.8/0.5 = -0.60. The product g1g2=−0.12g_1 g_2 = -0.12 is negative — outside [0,1][0,1]. The cavity is unstable: no confined Gaussian mode, so no clean TEM00_{00} beam. The point lands in a product-negative quadrant, where the green leaves never reach.

Step 2. Diagnose the culprit. Mirror 2 has R2=0.5R_2 = 0.5 m but L=0.8L = 0.8 m >R2> R_2, driving g2g_2 negative; one negative and one positive gg forces a negative product. Shorten the cavity so both g's share a sign. Try L=0.4L = 0.4 m: g1=0.60g_1 = 0.60, g2=1−0.4/0.5=0.20g_2 = 1 - 0.4/0.5 = 0.20, product g1g2=0.12∈[0,1]g_1 g_2 = 0.12 \in [0,1] — now stable.

Step 3. Size the repaired mode at λ=1064\lambda = 1064 nm, L=0.4L = 0.4 m, with Lλ/π=1.355×10−7 m2L\lambda/\pi = 1.355\times10^{-7}\ \mathrm{m}^2 and 1−g1g2=0.881 - g_1g_2 = 0.88. Mirror 1: 0.20/(0.60×0.88)=0.6155\sqrt{0.20/(0.60\times0.88)} = 0.6155, so w1=0.289 mmw_1 = \mathbf{0.289\ mm}. Mirror 2: 0.60/(0.20×0.88)=1.846\sqrt{0.60/(0.20\times0.88)} = 1.846, so w2=0.500 mmw_2 = \mathbf{0.500\ mm}. A thermal lens in the pumped rod effectively shortens RR and shifts both g's, so a high-power design keeps margin here against drifting across the boundary as pump power ramps. Original cavity unstable (g1g2=−0.12g_1g_2 = -0.12); shortening to L=0.4L = 0.4 m restores stability (g1g2=0.12g_1g_2 = 0.12), giving mode radii 0.2890.289 mm and 0.5000.500 mm.

Physics of the portrait

The stability diagram looks timeless, but it is barely older than the laser. After 1958 the open resonator — two mirrors with the sides thrown wide open — was a genuine puzzle: unlike a closed microwave cavity it has a continuum of diffraction loss, and it was not obvious that any low-loss mode should survive. The people who answered that built the entire language this article speaks.

Origins · 1 minHerwig Kogelnik and how this was discovered
HK
Herwig Kogelnik
born 1932, Austrian-American physicist at Bell Labs · laser resonators and integrated optics

Herwig Kogelnik trained in Vienna and Oxford and joined Bell Telephone Laboratories in 1961, exactly as the laser became a technology in want of a theory of its own cavities. The open resonator was the pressing problem: Fox and Li had just shown numerically (1961) that iterating the diffraction integral converges to a self-reproducing field with finite loss, and Boyd and Gordon had solved the confocal case analytically, but a general design theory was missing. Kogelnik supplied it. With Boyd he generalized the confocal analysis to arbitrary curvatures and spacings, and in the landmark 1966 review "Laser Beams and Resonators", written with Tingye Li, he distilled the subject into the tools physicists still use: the complex beam parameter qq, the ABCD law that propagates it, and the compact criterion 0≤g1g2≤10 \le g_1 g_2 \le 1 read straight off the g-chart. That review turned resonator design from solving diffraction integrals into multiplying a few 2×22\times2 matrices. Kogelnik went on to help found integrated optics, but the g-diagram remains the piece of his work every laser course teaches — a result compressed until its statement, two numbers and one inequality, fits on a single line.

History
1958
— Schawlow and Townes propose the optical maser, raising the open-resonator mode question
1961
— Fox and Li compute self-reproducing cavity modes by iterating the diffraction integral
1961
— Boyd and Gordon solve the confocal resonator modes analytically
1962
— Boyd and Kogelnik extend the theory to general stable resonators
1966
— Kogelnik and Li publish "Laser Beams and Resonators", codifying the ABCD method and the g-diagram

Going further

Texts. Saleh and Teich, Fundamentals of Photonics, 3rd ed., Ch. 10–11 (Resonator Optics) — the cleanest modern treatment of the g-parameters · Siegman, Lasers, Ch. 19–21 (stable and unstable resonator theory) — the definitive reference · Hecht, Optics, 5th ed., §13.1 (laser resonator geometry) · Milonni and Eberly, Laser Physics, Ch. 7 (cavity modes) · Demtröder, Laser Spectroscopy, Vol. 1, §5.2 (open optical resonators) · Svelto, Principles of Lasers, Ch. 5 (passive optical resonators).

Provenance. The definition gi=1−L/Rig_i = 1 - L/R_i and the inequality 0≤g1g2≤10 \le g_1 g_2 \le 1 follow Kogelnik and Li's 1966 review and are reproduced in Saleh and Teich (Ch. 11) and Siegman (Ch. 19); the ABCD self-consistency route is standard in both. The named-cavity landmarks and the confocal waist w0=Lλ/2πw_0 = \sqrt{L\lambda/2\pi} trace to Boyd and Gordon (1961) and Boyd and Kogelnik (1962). A verifier should confirm: (1) multiplying the four round-trip ray matrices symbolically gives half-trace 12(A+D)=2g1g2−1\tfrac12(A+D) = 2g_1g_2 - 1 (checked here numerically for symmetric and asymmetric cases); (2) the Gaussian self-consistency discriminant (A+D)2−4<0(A+D)^2 - 4 < 0 yields the same inequality; (3) the confocal waist from zR=L/2z_R = L/2; (4) the worked mode sizes w1=0.150w_1 = 0.150 mm, w2=0.368w_2 = 0.368 mm for the HeNe cavity and the stability sign flip in the Nd:YAG example. No numerical constants beyond standard textbook geometry are asserted.

Prerequisites. Gaussian Beams · ABCD Matrix Beam Propagation · Optical Resonators and Cavity Modes. Leads to. Transverse Cavity Modes (Hermite–Gaussian) · Laser Threshold Condition · Cavity Q Factor and Photon Lifetime.

Ask the tutor: Give me your mirror radii and spacing and I will plot your cavity on the g-diagram, tell you how much stability margin it has, and size the mode on each mirror.

What exactly does the g-parameter measure?

It measures how much a mirror focuses relative to the cavity length. In gi=1−L/Rig_i = 1 - L/R_i, a flat mirror (R=∞R=\infty) gives g=1g=1 and no focusing; a mirror with R=LR=L gives g=0g=0; push RR below LL and gg goes negative, meaning it over-focuses relative to the gap. The beauty is that only the product g1g2g_1 g_2 enters the stability test, so the two mirrors' focusing combines into one number deciding whether the round trip nets out as a rotation (stable) or a magnification (unstable).

Why does only the product g₁g₂ decide stability, not the individual g's?

Because stability is a property of the round trip, and the round-trip half-trace collapses to 2g1g2−12g_1g_2 - 1 — a function of the product alone. Physically, what matters is the total focusing accumulated in one circuit, and the two mirrors' contributions multiply rather than add in that trace, which is why a two-dimensional chart suffices. Be careful, though: the product classifies survival, but the individual g1g_1 and g2g_2 set the mode's shape — the spot on each mirror and where the waist sits. Two cavities with equal g1g2g_1g_2 can hold very differently shaped modes.

Is a cavity deep inside the stable region always the best choice?

No — "deepest inside" (near confocal) buys the smallest mode and the greatest alignment tolerance, but a small mode may fill the gain medium poorly and cap output. Closer to a boundary the mode grows, overlapping more gain and raising power, at the price of steeper diffraction loss and knife-edge alignment. There is no universal optimum; the g-diagram gives the trade-off, not the answer. A HeNe alignment laser lives near confocal for robustness; a high-power rod laser edges toward a boundary for mode volume while leaving margin for its thermal lens.

What happens right on a boundary, like the plane-parallel cavity?

On a boundary — g1g2=0g_1 g_2 = 0 or 11 — the cavity is marginally stable: the round-trip eigenvalues sit exactly at ±1\pm 1, poised between rotation and magnification. The mathematics still admits a mode, but a limiting, infinitely sensitive one. The plane-parallel cavity at (1,1)(1,1) is the classic case: its flat mirrors impose no restoring focusing, so the "mode" is essentially a plane wave that only finite mirror size and diffraction keep marginally bounded. This is why plane-parallel resonators are notoriously twitchy to align, and why real designs back off a boundary into the interior.

Why would anyone build an unstable resonator on purpose?

Because for a high-gain, high-energy laser a small stable Gaussian mode is a liability: it samples only a thin thread of the gain medium and clamps the output. An unstable resonator (g1g2>1g_1 g_2 > 1 or <0< 0) has a magnifying round trip, so the beam expands to fill the whole gain aperture each pass. The light that magnifies past the small feedback mirror spills around it — that geometric leakage is the output coupling, no partial mirror needed. Pulsed CO2_2, excimer, and large solid-state systems use this routinely; the trade-off is a doughnut-shaped near field rather than a clean Gaussian, acceptable when raw energy beats beam quality.

How does the g-parameter connect to the transverse-mode frequencies?

Through the round-trip Gouy phase θ\theta, with cos⁡θ=2g1g2−1\cos\theta = 2g_1g_2 - 1, i.e. θ=2arccos⁡g1g2\theta = 2\arccos\sqrt{g_1g_2}. This is the extra advance a Gaussian picks up per trip beyond the plane-wave part, and it separates transverse modes (TEM00_{00}, TEM01_{01}, …) in frequency: higher-order modes accumulate integer multiples of θ\theta, so their spacing is θ/2π\theta/2\pi times the free spectral range. At confocal (g1g2=0g_1g_2 = 0) θ=π\theta = \pi and the transverse modes fall exactly halfway between longitudinal ones — the confocal degeneracy. Outside the stable range θ\theta turns imaginary, the signature of runaway. So the geometry that sets stability also sets the transverse spectrum.

Before this

Where this leads