Cavity Layout β€” Folded Beam Path
Tangential w Sagittal w
Beam width shown schematically, exaggerated for visibility (not to length scale). Drag the small handles along any arm to change its length.
Round-Trip Matrix
Beam Parameters
Astigmatism
Stability Diagram
Beam Profile w(z)
Sweep & Optimize
Alignment Sensitivity
Mode Patterns
Q-Switch
MΒ² & Coupling
πŸ”₯ Thermal Lens
🎲 Monte Carlo
⌨ Scripting
Round-trip ABCD Matrix
M = [
β€”β€” β€”β€”
]
Det(M) = ADβˆ’BCβ€”
(A+D)/2 (trace/2)β€”
Stability |trace/2| < 1β€”
Both Planes at a Glance

Tangential (yz)

Sagittal (xz)

Koechner eq. 5.69: a fold mirror used off-axis has ftangential=(R/2)cosΞΈ and fsagittal=(R/2)/cosΞΈ β€” the two planes generally see different stability.

Individual Element Matrices (selected plane)
Equivalent g-Parameters & Cavity Facts
Round-trip lengthβ€”
Free spectral rangeβ€”
Gouy phase / round tripβ€”
Transverse Mode Frequency Degeneracy

Transverse modes (m,n) sit at Ξ½ = FSRΓ—[q + (m+Β½)ψt/2Ο€ + (n+Β½)ψs/2Ο€] relative to the longitudinal comb. When ψ/2Ο€ is a simple fraction p/N, mode families from different q overlap in frequency β€” the well-known transverse-mode-degeneracy effect (strongest at confocal, ψ=90Β°).

Eigenmode at Reference Plane (M1)

Tangential

Sagittal

Beam Ellipticity

Ratio wsagittal/wtangential at the reference plane. A perfectly round beam gives 1.00; folded cavities with fold mirrors are astigmatic and deviate from 1.00 unless compensated.

Mode Size at Each Optical Element
Elementz (mm) wt (mm)ws (mm) Area (mmΒ²)
Crystal Face Mode Sizes (input / center / output)
Fold Mirror Astigmatism (Koechner eq. 5.69)
Brewster Element Path Asymmetry (eq. 5.70)
Astigmatic Compensation Check (eq. 5.71)
g₁–gβ‚‚ Diagram (both planes)

✦ tangential Β· ✚ sagittal. Shaded = stable (0<g₁gβ‚‚<1 or βˆ’1<g₁gβ‚‚<0). This point uses only the two terminal mirrors β€” it does NOT move for fold-mirror curvature or thermal lensing. For the real stability of the full cavity, use the "effective g₁gβ‚‚" numbers below and the header badge, not this point's position.

Sensitivity β€” Scan Total Cavity Length
1.000Γ—
Beam Radius w(z) Along Unfolded Path β€” Both Planes
Rayleigh Range & Divergence
Diffraction Limit
Parameter Sweep (RP-Resonator-style dioptric-power / wavelength scan)

Mode Radius vs Swept Parameter
wβ‚€ tangentialwβ‚€ sagittalβ–¨ unstable region
Local Optimizer (coordinate-descent figure-of-merit minimization)

Mark parameters to optimize with the 🎯 button next to any slider in the element list (left panel), then run.

Misalignment Setup (ABCDEF extended matrices)
mrad
mm

Solves the self-consistent axis offset (Iβˆ’M)⁻¹·(E,F) induced by a fixed per-round-trip perturbation, then propagates it through the cavity with ordinary ABCD matrices. Tilt and decenter of the selected mirror are applied together (2Γ— doubling on reflection for each, matching the existing tilt convention); set either to 0 to isolate the other's effect.

Induced Beam-Axis Offset per Element
Higher-Order Transverse Mode TEMmn

Pattern uses I_mn(x,y) = [H_m(√2x/w_t)]Β² [H_n(√2y/w_s)]Β² exp(βˆ’2xΒ²/w_tΒ²βˆ’2yΒ²/w_sΒ²), evaluated at the cavity waist (M1 reference plane) using the current wβ‚€β‚œ, wβ‚€β‚›. MΒ² of a pure TEMmn mode is exactly (2m+1) in the tangential direction and (2n+1) in the sagittal direction β€” a standard, closed-form result of Gaussian-beam theory (not a wave-optics simulation).

Q-Switch Design Calculator (Koechner eq. 8.6–8.13, optimum-coupling analysis)
%
mm
Γ—10⁻¹⁹ cmΒ²
mm

z is found from the desired output when in target mode by solving E_out/E_sc = zβˆ’1βˆ’ln z numerically. Defaults reproduce Koechner's own Nd:YAG worked example (Β§8.1): expect R_optβ‰ˆ0.65, t_pβ‰ˆ11ns, P_pβ‰ˆ9MW.

Results
Photon Lifetime & Cavity Q (Koechner eq. 3.8–3.9, using this cavity's actual OC)
MΒ² & External Beam Propagation
mm (from waist)

Scales the cavity's ideal TEMβ‚€β‚€ waist by MΒ² using the standard real-beam relations ΞΈ_MΒ²=MΒ²Ξ»/(Ο€wβ‚€) and w(z)=wβ‚€βˆš[1+(MΒ²Ξ»z/Ο€wβ‚€Β²)Β²]. Use this to predict the actual output spot size downstream when the real beam isn't diffraction-limited.

Pump/Mode Overlap (simplified, co-located Gaussians)
ΞΌm

Ξ· = [2wβ‚€w_p/(wβ‚€Β²+w_pΒ²)]Β² β€” assumes the pump and cavity mode share the same waist location with matched wavefronts. This is a fast estimate, not Koechner's full 3-D overlap integral (eq. 6.71), which additionally needs the pump's MΒ², absorption profile, and focus offset solved numerically.

Fresnel Number (for any Aperture elements in the cavity)

N_F = aΒ²/(Ξ»L). Large N_F (≫1) β†’ the aperture is many Fresnel zones across, diffraction loss is negligible and the paraxial ABCD treatment is trustworthy. N_F ≲ 1 β†’ significant diffraction, the simple clipping estimate elsewhere in this tool becomes unreliable.

Monte Carlo Tolerance Analysis
%
%
Β°

Independently perturbs every element's ROC, length/thickness, and (for fold mirrors) AOI with Gaussian noise at the Οƒ set above, then reruns the same round-trip stability/eigenmode solve used everywhere else in this tool for each trial. Angular tilt and lateral decenter of a single mirror are handled separately in the Alignment Sensitivity tab β€” this tab is for dimensional manufacturing/fabrication tolerances across the whole cavity.

Yield & Waist Statistics

Run the analysis to see results.

Macro Script (runs in your browser, against the current design)
Output

        
Available Helpers

setParam(label, key, value) β€” set a numeric parameter on the element named "label"
getParam(label, key) β€” read a parameter back
metrics() β€” {stableT, stableS, w0t, w0s (mm), g1, g2} for the current design
log(...args) β€” print to this Output panel
elementLabels() β€” array of current element labels, in order