Two-level Bloch-vector solver written from scratch. The plotted HHG spectrum is the dipole-power spectrum: I(ω) = |FFT[W(t)P(t)]|², using a super-Gaussian window with near-zero wings.
Ready.
Diagnostics will appear after calculation.
Time-domain field and dipole
Dipole spectrum: |FFT[P(t)]|²
This is intentionally minimal and explicit: no phase cycling, no subtraction-based nonlinear residual in the spectrum, no Hann/Hanning window, and no spectral smoothing. The perturbative odd harmonics are marked as H1, H3, H5, ... . The plotted spectrum is the dipole power spectrum without ω⁴ multiplication. For a simple two-level system, visible higher harmonics require sufficiently strong coupling dE₀ and/or resonance enhancement.