Results obtained by numerically integrating the Maxwell Bloch equations for a flat top pulse
============================================================================================

B S Cartwright, S A Wrathmall and R M Potvliege, "Self-induced transparency and optical transients in atomic vapors", accepted for publication in Physical Review A (2026).


Atomic parameters, pulse shape, atomic density, temperature, etc.: As described in this article.

Numerical results for the 2- and 3-state models
-----------------------------------------------

Data set 1: 2-state model, Iin = 5 W/cm2, 780 nm field only, no broadening.

Data set 2: 2-state model, Iin = 5 W/cm2, 780 nm field only, with Doppler broadening, no homogeneous broadening.

Data set 3: 2-state model, Iin = 15 W/cm2, 780 nm field only, with Doppler broadening, no homogeneous broadening.

Data set 4: 2-state model, Iin = 5 W/cm2, 780 nm field only, with Doppler broadening and homogeneous broadening.

Data set 5: 2-state moedl, Iin = 5 W/cm2, 795 nm field only, with Doppler broadening and homogeneous broadening.

Data set 6: 3-state model, superposition of a 780 nm field with Iin = 5 W/cm2 and a 795 nm field with Iin = 1 muW/cm2, with Doppler broadening and homogeneous broadening.

Data set 7: 3-state model, superposition of a 780 nm field with Iin = 5 W/cm2 and a 795 nm field with Iin = 5 W/cm2, with Doppler broadening and homogeneous broadening.

Data set 8: 2-state model, Iin = 0.5 W/cm2, 780 nm field only, with Doppler broadening and homogeneous broadening.

prop_k?.dat and prop_c?.dat files: Input files for the driveall program of the CoOMBE suite of codes, used for generating these numerical results. See R M Potvliege and S A Wrathmall, CoOMBE: A suite of open-source programs for the integration of the optical Bloch questions and Maxwell-Bloch equations, Comput. Phys. Commun. 106, 109374 (2025).

outamplitudes?.dat:  Electric field amplitudes, E.

Structure of these files for single-field calculations:
First column: time in mus; second column: z in mum; third and fourth columns: Re E and Im E in V/m.

Structure of these files for two-field calculations:
First column: time in mus; second column: z in mum; third and fourth columns: Re E and Im E for the 780 nm field in V/m; fifth and sixth columns: Re E and Im E for the 795 nm field, in V/m.


Numerical results from the calculations taking into the hyperfine structure of the states (Fig. 7 of the article)
-----------------------------------------------------------------------------------------------------------------
 
r2_dnoidal_80deg_all_10muW_truetime: data for Fig. 7b
r3_dnoidal_80deg_all_10muW_truetime: data for Fig. 7a
r2_dnoidal_80deg_all_both50: data for Fig. 7d
r3_dnoidal_80deg_all_both50: data for Fig. 7c

Structure of these files:
First column: t in mus; second column: z in mum; third column: normalised |E|.

The values of |E| given in these four files can be transformed into the corresponding intensity by dividing |E(z,t)|**2 by |E(z = 0,t = 4 ns)|**2 and multiplying by the relevant value of Iin, as stated in the article.
 
