Self-induced transparency and optical transients in atomic vapors
Pith reviewed 2026-05-19 13:36 UTC · model grok-4.3
The pith
Rapid turn-on of a strong resonant laser field can produce transient trains of damped solitons in atomic vapors before the system settles to a steady state.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The rapid turn-on of a strong, resonant, continuous wave laser field may trigger the formation of a transient oscillation akin to a train of damped solitons, before the vapor-field system relaxes into a stationary state. We study this transient dynamic on theoretical models of a rubidium vapor. We also consider doubly resonant V-systems, for which the transients take the form of trains of damped simultons. We compute the propagating field(s) by solving the Maxwell-Bloch equations, taking homogeneous broadening, Doppler broadening and the full hyperfine structure of the atoms into account. We also compare the actual fields to the stationary dnoidal fields predicted by the Maxwell-Bloch in the
What carries the argument
Numerical integration of the Maxwell-Bloch equations that include homogeneous and Doppler broadening plus the full hyperfine structure, used to generate the time-dependent propagating fields and compare them with the stationary dnoidal solutions of self-induced transparency.
If this is right
- The transients appear as trains of damped solitons for a two-level transition.
- In doubly resonant V-systems the transients take the form of trains of damped simultons.
- After the transient the propagating fields approach the stationary dnoidal fields of self-induced transparency.
- The same transient dynamics should occur in any atomic vapor when the turn-on is sufficiently rapid relative to relaxation times.
Where Pith is reading between the lines
- The effect may appear whenever a strong resonant laser is gated on or off in laboratory vapor cells, affecting the early-time behavior of optical signals.
- Extensions that include velocity-changing collisions or spatial inhomogeneities could be tested by adding those terms to the same Maxwell-Bloch integration.
- The duration and amplitude of the transient train should scale with the ratio of the Rabi frequency to the relaxation rates.
Load-bearing premise
The laser turn-on must be fast compared with atomic relaxation rates, and the Maxwell-Bloch equations with only homogeneous and Doppler broadening plus hyperfine structure must capture the dominant physics.
What would settle it
Time-resolved measurement of the transmitted intensity immediately after a fast resonant turn-on of a strong continuous-wave laser through a rubidium vapor cell, checking for the predicted train of damped intensity oscillations before the field reaches its stationary value.
Figures
read the original abstract
The rapid turn-on of a strong, resonant, continuous wave laser field may trigger the formation of a transient oscillation akin to a train of damped solitons, before the vapor-field system relaxes into a stationary state. We study this transient dynamic on theoretical models of a rubidium vapor. We also consider doubly resonant V-systems, for which the transients take the form of trains of damped simultons. We compute the propagating field(s) by solving the Maxwell-Bloch equations, taking homogeneous broadening, Doppler broadening and the full hyperfine structure of the atoms into account. We also compare the actual fields to the stationary dnoidal fields predicted by the Maxwell-Bloch equations in conditions of self-induced transparency. A similar dynamics is expected to occur in any atomic vapor at the turn-on of a strong resonant continuous wave field provided the turn-on is sufficiently fast compared to relaxation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims that the rapid turn-on of a strong resonant continuous-wave laser field in atomic vapors (modeled on rubidium) can trigger transient oscillations resembling trains of damped solitons before relaxation to a stationary state. These dynamics are obtained by direct numerical integration of the Maxwell-Bloch equations that incorporate homogeneous broadening, Doppler broadening, and the full hyperfine structure. For doubly resonant V-systems the transients appear as trains of damped simultons. The computed propagating fields are compared to the stationary dnoidal solutions of the Maxwell-Bloch equations under self-induced transparency conditions. The effect is asserted to occur in any atomic vapor provided the turn-on is sufficiently fast relative to relaxation times.
Significance. If the central numerical results hold, the work identifies a previously under-examined transient regime at the onset of self-induced transparency. The direct integration of the Maxwell-Bloch equations without fitted parameters or circular definitions supplies a non-circular foundation for the claim. Inclusion of realistic Doppler broadening and hyperfine structure improves experimental relevance, while the explicit comparison to stationary dnoidal fields usefully demarcates the transient window from the long-time attractor.
major comments (2)
- [Abstract] Abstract and conclusion: the central claim requires that the laser turn-on be 'sufficiently fast' compared to atomic relaxation times, yet no specific rise time, functional form of the turn-on, or scan over rise times is supplied to show when the damped-soliton train disappears. Without this threshold the applicability to laboratory lasers remains unquantified and the claim cannot be assessed for realistic conditions.
- [Numerical methods] Numerical methods / results sections: the manuscript provides no information on numerical stability, time-step or grid-size convergence tests, or quantitative error measures for the transient trains. These checks are load-bearing for establishing that the reported damped-soliton-like oscillations are physical rather than numerical artifacts.
minor comments (1)
- Figure captions and axis labels should explicitly indicate the time window over which the transient train is visible versus the stationary regime to aid reader interpretation.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the positive assessment of its significance. We respond to each major comment below and outline the revisions we will make to address the concerns.
read point-by-point responses
-
Referee: [Abstract] Abstract and conclusion: the central claim requires that the laser turn-on be 'sufficiently fast' compared to atomic relaxation times, yet no specific rise time, functional form of the turn-on, or scan over rise times is supplied to show when the damped-soliton train disappears. Without this threshold the applicability to laboratory lasers remains unquantified and the claim cannot be assessed for realistic conditions.
Authors: We agree that a more quantitative characterization of the required turn-on speed would improve the manuscript's relevance to laboratory conditions. In the revised version we will specify the functional form of the turn-on (a smooth, monotonic ramp) used in the simulations, report the rise times employed, and add a brief parameter scan showing the disappearance of the damped-soliton train as the rise time approaches the atomic relaxation timescales. These results will be summarized in the abstract and discussed in the conclusion. revision: yes
-
Referee: [Numerical methods] Numerical methods / results sections: the manuscript provides no information on numerical stability, time-step or grid-size convergence tests, or quantitative error measures for the transient trains. These checks are load-bearing for establishing that the reported damped-soliton-like oscillations are physical rather than numerical artifacts.
Authors: We acknowledge that explicit documentation of numerical convergence is necessary to substantiate the physical origin of the transients. We will expand the numerical methods section to describe the integration algorithm, the spatial and temporal grid parameters, and the results of systematic convergence tests. These tests will demonstrate that the amplitude, frequency, and damping of the reported oscillations remain unchanged under successive refinement of the time step and grid spacing, together with quantitative error estimates (relative L2 norm of the field) for the transient regime. revision: yes
Circularity Check
No significant circularity in numerical Maxwell-Bloch simulations
full rationale
The paper obtains its results on transient damped-soliton-like oscillations by direct numerical integration of the Maxwell-Bloch equations that incorporate homogeneous broadening, Doppler broadening, and full hyperfine structure for rubidium and V-systems. No parameters are fitted to data subsets and then presented as independent predictions, no self-citations supply load-bearing uniqueness theorems or ansatzes, and no known empirical patterns are merely renamed. The transient dynamics and comparisons to stationary dnoidal fields follow from the stated initial conditions and the equations themselves, making the derivation self-contained without reduction to its own inputs.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math The Maxwell-Bloch equations accurately describe the coherent interaction between the propagating laser field and the atomic medium in this regime.
- domain assumption Homogeneous broadening, Doppler broadening, and full hyperfine structure are the relevant physical effects to include for rubidium vapor.
Reference graph
Works this paper leans on
-
[1]
II A to the case considered in Sec
Two-color fields This appendix extends Sec. II A to the case considered in Sec. IV, namely that of two fields addressing a 3-level system in a V configuration. As usual in this context, we refer to one of these fields as the probe field and to the other as the coupling field. The former and the latter couple the ground state (state 1) to two different exc...
-
[2]
The Maxwell-Bloch equations without homogeneous broadening As noted in Sec. II, calculating the density operator ˆρ(z, t) involves calculating the velocity-specific density operators ˆρv(x, t, uz) for all the relevant velocity classes. Since we neglect homogeneous broadening here, we set ˆρv(z, t, uz) ≡ |Ψ(z, t)⟩⟨Ψ(z, t)| (A17) with |Ψ(z, t)⟩ = c1(z, t)|1...
-
[3]
Z τ K(k) −τ K(k) Ωp(z, t) dt #2 +
Equal propagation coefficients The above equations can be solved in closed form when the propagation coefficients µp and µc are equal. In par- ticular, setting µp = µc = µ. (A30) makes it possible to seek solutions of either Eqs. (A20a)– (A20e) or Eqs. (A22a)–(A22e) for which the two fields are in a constant ratio, i.e., solutions for which Ωp(z, t) ≡ rΩc...
-
[4]
Exact results applicable to the case where µp ̸= µc are scarce
Unequal propagation coefficients Whether inhomogeneous broadening is taken into ac- count or not, Ω p(z, t) and Ω c(z, t) cannot be obtained in 13 closed form without further approximation if the propa- gation coefficients µp and µc are unequal. Exact results applicable to the case where µp ̸= µc are scarce. To the best of our knowledge, they are limited ...
-
[5]
S. L. McCall and E. L. Hahn, Phys. Rev. Lett. 18, 908 (1967); Phys. Rev. 183, 457 (1969)
work page 1967
-
[6]
M. D. Crisp, Phys. Rev. Lett. 22, 820 (1969)
work page 1969
-
[7]
J. H. Eberly, Phys. Rev. Lett. 22, 760 (1969)
work page 1969
-
[8]
G. L. Lamb, Jr., Rev. Mod. Phys. 43, 99 (1971)
work page 1971
-
[9]
L. A. Bol’shov and V. V. Likhanski˘ ı, Kvantovaya Elek- tron. 12, 1339 (1985) [Sov. J. Quantum Electron.15, 889 (1985)]
work page 1985
-
[10]
A. I. Maimistov, A. M. Basharov, S. O. Elyutin, and Yu. M. Sklyarov, Phys. Rep. 191, 1 (1990); A. I. Maimistov and A. M. Basharov, Nonlinear Optical Waves , Kluwer, Dordrecht (1999)
work page 1990
-
[11]
For an introduction to the early work on SIT, see, e.g., L. Allen and J. H. Eberly,Optical Resonances and Two-level Atoms, Wiley, New York (1975)
work page 1975
-
[12]
T. P. Ogden, K. A. Whittaker, J. Keaveney, S. A. Wrath- mall, C. S. Adams, and R. M. Potvliege, Phys. Rev. Lett. 123, 243604 (2019)
work page 2019
-
[13]
M. V. Arkhipov, A. A. Shimko, N. N. Rosanov, I. Babushkin, and R. M. Arkhipov, Phys. Rev. A 101, 013803 (2020)
work page 2020
-
[14]
Z. Bai, C. S. Adams, G. Huang, and W. Li, Phys. Rev. Lett. 125, 263605 (2020)
work page 2020
-
[15]
See, e.g., C. Jirauschek, M. Riesch, and P. Tzenov, Adv. Theory Simul. 2, 1900018 (2019); G. T. Adamashvili, Eur. Phys. J. D 74, 41 (2020); S. K. Hazra, P. K. Pathak, and T. N. Dey, Phys. Rev. B 107 235409 (2023); M. S. Najafabadi, L. L. S´ anchez-Soto, J. F. Corney, N. Kalinin, A. A. Sorokin, and G. Leuchs, Phys. Rev. Res. 6, 023142 (2024)
work page 2019
-
[16]
See, e.g., S. Hughes, Phys. Rev. Lett. 81, 3363 (1998); D. V. Novitsky, Phys. Rev. A 84, 013817 (2011); H. Wu, J. Tang, M. Chen, M. Xiao, Y. Lu, K. Xia, and F. Nori, Opt. Express 32, 11010 (2024); A. Pakhomov, Phys. Rev. A 111, 013502 (2025)
work page 1998
-
[17]
A. M. Alhasan, J. Fiutak and W. Miklaszewski, Z. Phys. B 88, 349 (1992)
work page 1992
- [18]
-
[19]
M. J. Konopnicki and J. H. Eberly, Phys. Rev. A 24, 2567 (1981)
work page 1981
-
[20]
See also M. J. Konopnicki, P. H. Drummond, and J. H. Eberly, Opt. Commun. 36, 313 (1981); C. R. Stroud, Jr. and D. A. Cardimona, Opt. Commun. 37, 221 (1981)
work page 1981
-
[21]
Many-state V-systems may also support many-color si- multons — see, e.g., G. Huang and C. Hang, Phys. Lett. A 354, 406 (2006); G. Huang, C. Hang and L. Deng, Eur. Phys. J. D 40, 437 (2006)
work page 2006
-
[22]
F. T. Hioe and R. Grobe, Phys. Rev. Lett. 73, 2559 (1994)
work page 1994
-
[23]
M. D. Crisp, Phys. Rev. A 5, 1365 (1972)
work page 1972
-
[24]
B. S´ egard, B. Macke, J. Zemmouri, and W. Sergent, Ann. Phys. (Paris) 15, 167 (1990)
work page 1990
- [25]
-
[26]
M. D. Crisp, Phys. Rev. A 1, 1604 (1970); ibid. 2, 2172 15 (1970)
work page 1970
-
[27]
J. E. Rothenberg, D. Grischkowsky, and A. C. Balant, Phys. Rev. Lett. 53, 552 (1984)
work page 1984
- [28]
-
[29]
E. M. Pessina, B. S´ egard, and B. Macke, Opt. Comm. 81, 397 (1991)
work page 1991
-
[30]
W. R. LeFew, S. Venakides, and D. J. Gauthier, Phys. Rev. A 79, 063842 (2009)
work page 2009
-
[31]
D. Wei, J. F. Chen, M. M. T. Loy, G. K. L. Wong, and S. Du, Phys. Rev. Lett. 103, 093602 (2009)
work page 2009
-
[32]
K. E. Oughstun, N. A. Cartwright, D. J. Gauthier, and H. Jeong, J. Opt. Soc. Am. B 27, 1664 (2010); B. Macke and B. S´ egard,ibid. 28, 450 (2011); K. E. Oughstun, N. A. Cartwright, D. J. Gauthier, and H. Jeong, ibid. 28, 468 (2011)
work page 2010
- [33]
- [34]
-
[35]
D. J. Kaup, Phys. Rev. A 16, 704 (1977)
work page 1977
-
[36]
M. A. Newbold and G. J. Salamo, Phys. Rev. Lett. 42, 887; J. L. Shultz and G. J. Salamo, Phys. Rev. Lett. 78, 855 (1997); M. O. Scully, G. S. Agarwal, O. Kocharovskaya, V. V. Kozlov, and A. B. Matsko, Opt. Express 8, 66 (2001); S. M. Saadeh, J. L. Shultz, and G. J. Salamo, Opt. Express 8, 153 (2001)
work page 1997
-
[37]
B. S. Cartwright, Optical transients in atomic vapours, Durham theses, Durham University (2022), http://etheses.dur.ac.uk/14532/
work page 2022
-
[38]
We follow the NIST Digital Library of Mathematical Functions (https://dlmf.nist.gov/) in the definition of this function: K(z) = Z π/2 0 dθp 1 − z2 sin2 θ
-
[39]
They need to be supple- mented by a dispersion relation in more general cases [1]
These results, as stated, only apply to the case where g(∆) is an even function of ∆. They need to be supple- mented by a dispersion relation in more general cases [1]
-
[40]
V. V. Kozlov and E. B. Kozlova, Opt. Spektrosk. 107, 139 (2009) [Opt. Spectrosc. 107, 129 (2009)]
work page 2009
-
[41]
L. A. Bol’shov, N. N. Elkin, T. K. Kirichenko, V. V. Likhanski˘ ı, and A. P. Napartovich, Kvantovaya Elek- tron. 9, 1476 (1982) [Sov. J. Quantum Electron. 12, 941 (1982)]
work page 1982
-
[42]
L. A. Bol’shov, N. N. Elkin, T. K. Kirichenko, V. V. Likhanski˘ ı, and M. I. Persiantsev, Preprint IAE-3732/16, Atomic Energy Inst., Moscow, 1983 [in Russian]
work page 1983
-
[43]
L. A. Bol’shov, N. N. Yelkin, V. V. Likhanski˘ ı, and M. I. Persiantsev, Zh. Eks. Teor. Fiz. 94, 101 (1988) [Sov. Phys. JETP 67, 2013 (1988)]
work page 1988
-
[44]
V. V. Kozlov and E. E. Fradkin, Pis’ma Zh. Eksp. Teor. Fiz. 68, 359 (1998) [JETP Lett. 68, 383 (1998)]
work page 1998
-
[45]
N. V. Denisova, V. S. Egorov, V. V. Kozlov, N. M. Reutova, P. Yu. Serdobintsev, and E. E. Fradkin, Zh. Eksp. Teor. Fiz. 113, 71 (1998) [J. Exp. Theor. Phys. 86, 39 (1998)]
work page 1998
-
[46]
V. V. Kozlov, P. G. Polynkin, and M. O. Scully, Phys. Rev. A 59, 3060 (1999)
work page 1999
-
[47]
E. Paspalakis, N. J. Kylstra, and P. L. Knight, Phys. Rev. A 61, 045802 (2000)
work page 2000
-
[48]
V. V. Kozlov and E. B. Kozlova, Opt. Spektrosk. 108, 824 (2010) [Opt. Spectrosc. 108, 780 (2010)]
work page 2010
-
[49]
O. M. Fedotova, O. K. Khasanov, G. A. Rusetsky, J. Degert, and E. Freysz, Phys. Rev. A 90, 053843 (2014)
work page 2014
-
[50]
R. M. Potvliege and S. A. Wrathmall, Comput. Phys. Commun. 306, 109374 (2025)
work page 2025
-
[51]
D. A. Steck, Rubidium 85 D line data, available online at http://steck.us/alkalidata
-
[52]
H. P. Grieneisen, J. Goldhar, N. A. Kurnit, and A. Javan, Appl. Phys. Lett. 21, 559 (1972)
work page 1972
-
[53]
S. M. Hamadani, J. Goldhar, N. A. Kurnit, and A. Javan, Appl. Phys. Lett. 25, 160 (1974)
work page 1974
-
[54]
M. Matusovsky, B. Vaynberg, and M. Rosenbluh, J. Opt. Soc. Am. B 13, 1994 (1996)
work page 1994
-
[55]
4(b) and α = (1.2 mm)−1 in the case of Fig
According to the Beer-Lambert law, the intensity of the probe field would decrease like I in p exp(−αz) with α = (0.42 mm)−1 in the case of Fig. 4(b) and α = (1.2 mm)−1 in the case of Fig. 7, in the absence of the coupling field. See, e.g., Ref. [46] for the calculation of the absorption coefficient α with Doppler broadening
-
[56]
A. Rahman, Phys. Rev. A 60, 4187 (1999); A. Rahman and J. H. Eberly, Opt. Express 4, 133 (1999)
work page 1999
-
[57]
[2] denotes the negative of the detuning ∆ defined in the present work
These equations take into the fact that the symbol ∆ used in Ref. [2] denotes the negative of the detuning ∆ defined in the present work
-
[58]
E.g., C. R. Higgins and I. G. Hughes, J. Phys. B: At. Mol. Opt. Phys. 54, 165403 (2021)
work page 2021
-
[59]
T. Y. Abi-Salloum, Phys. Rev. A 81, 053836 (2010)
work page 2010
-
[60]
S. Khan, V. Bharti, and V. Natarajan, Phys. Lett. A 380, 4100 (2016)
work page 2016
-
[61]
Quantum interference plays a role for weak coupling fields, though. See C. Zhu, C. Tan and G. Huang, Phys. Rev. A 87, 043813 (2013)
work page 2013
- [62]
-
[63]
Doubly dressed states have been considered previously for ladder systems [see, e.g., R.-Y. Chang, W.-C. Fang, Z.-S. He, B.-C. Ke, P.-N. Chen, and C.-C. Tsai, Phys. Rev. A 76, 053420 (2007)] but, to our knowledge, not for V-systems
work page 2007
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.