UltrafastEngine

15 -- Spatially resolved gain and the thermal lens

A 2 ps, 100 µJ pulse at 10 kHz crosses a 4 mm Yb:YAG rod once, end-pumped by 60 W at 940 nm (super-Gaussian, 0.6 mm, double pass). The rod is modelled on a radial grid, and then the beam travels 0.4 m in free space.

15 -- Spatially resolved gain and the thermal lens

Top row 60 W pump, bottom row 5 W. Left to right: the beam in the rod (the input is a collimated Gaussian), the wavefront it leaves with (a parabola: the thermal lens), the beam 0.4 m behind the rod, and the Zernike content of the wavefront (the blue bar is defocus, Noll 4). At 60 W the lens converges the beam, so the peak fluence 0.4 m behind the rod is higher than in it (colour bars); at 5 W the lens is 10 times weaker. python make_figures.py runs the commands below and redraws it.

ultrafast run rod.ufs run.json -o results-60W
ultrafast run rod.ufs run-lowpump.json -o results-5W     # same system, 5 W pump (a run-spec override)
ultrafast-plot results-60W

On the grid:

  • every transverse point has its own inversion, so the gain and depletion follow the pump and beam profiles;
  • the beam diffracts inside the rod;
  • the heat load (quantum defect × absorbed pump) gives ΔT(r, z) and Δn = dn/dT·ΔT, applied as a position-dependent index (docs/spatial-properties.md).
60 W pump 5 W pump
absorbed pump / heat 44.3 W / 3.87 W 4.05 W / 0.35 W
on-axis inversion 0.20 0.02
energy gain 1.42 0.92 (reabsorption: Yb:YAG is quasi-three-level)
thermal lens f_th 0.54 m 5.9 m
wavefront after the rod (R) −0.51 m (converging) −3.7 m
beam radius after the rod / after 0.4 m 471 µm / 319 µm 497 µm / 524 µm
M² 1.034 (aberrated lens: the pump is not Gaussian) 1.002

gain.png shows the inversion across the beam (flat inside the pump, zero outside). The pulse train converges in ≈ 70 pulses (τ = 0.95 ms ≫ 100 µs period). The far-field camera's Strehl ratio (0.41 at 60 W) is computed against a flat-phase beam, so it includes the defocus of the thermal lens as well as its aberration.