pith. sign in

arxiv: 2505.17930 · v2 · submitted 2025-05-23 · ⚛️ physics.atom-ph · nlin.PS

Self-induced transparency and optical transients in atomic vapors

Pith reviewed 2026-05-19 13:36 UTC · model grok-4.3

classification ⚛️ physics.atom-ph nlin.PS
keywords self-induced transparencyoptical transientsMaxwell-Bloch equationsatomic vaporsrubidiumsolitonssimultonsDoppler broadening
0
0 comments X

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.

The paper examines what happens when a strong continuous-wave laser is switched on quickly while resonant with an atomic transition in a vapor. It finds that the propagating light field develops oscillations resembling a series of damped solitons that eventually decay, leaving the system in the stationary state predicted by self-induced transparency theory. The authors reach this conclusion by integrating the Maxwell-Bloch equations for rubidium, incorporating homogeneous and Doppler broadening together with the full hyperfine structure. They also treat doubly resonant V-systems, where the transients instead appear as trains of damped simultons. The same behavior is expected in any atomic vapor provided the laser turn-on is fast compared with the relaxation times.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2505.17930 by B. S. Cartwright, R. M. Potvliege, S. A. Wrathmall.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) The applied field in the 2-state model discussed [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Solid curve: The intensity of the field at a distance [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The variation of the field in the 2-state model dis [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Coupling field with [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The fields and the eigenenergies of the rotating-wave [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. The coupling field and the probe field as obtained [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The ratio Ω [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The claim rests on the applicability of the Maxwell-Bloch framework to the transient regime and the assumption that included broadening mechanisms dominate; no free parameters or new entities are introduced in the abstract.

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.
    Invoked as the computational foundation for all results in the abstract.
  • domain assumption Homogeneous broadening, Doppler broadening, and full hyperfine structure are the relevant physical effects to include for rubidium vapor.
    Explicitly stated as accounted for in the models.

pith-pipeline@v0.9.0 · 5682 in / 1529 out tokens · 64451 ms · 2026-05-19T13:36:40.162875+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages

  1. [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. [2]

    II, calculating the density operator ˆρ(z, t) involves calculating the velocity-specific density operators ˆρv(x, t, uz) for all the relevant velocity classes

    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. [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. [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. [5]

    S. L. McCall and E. L. Hahn, Phys. Rev. Lett. 18, 908 (1967); Phys. Rev. 183, 457 (1969)

  6. [6]

    M. D. Crisp, Phys. Rev. Lett. 22, 820 (1969)

  7. [7]

    J. H. Eberly, Phys. Rev. Lett. 22, 760 (1969)

  8. [8]

    G. L. Lamb, Jr., Rev. Mod. Phys. 43, 99 (1971)

  9. [9]

    L. A. Bol’shov and V. V. Likhanski˘ ı, Kvantovaya Elek- tron. 12, 1339 (1985) [Sov. J. Quantum Electron.15, 889 (1985)]

  10. [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)

  11. [11]

    Allen and J

    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)

  12. [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)

  13. [13]

    M. V. Arkhipov, A. A. Shimko, N. N. Rosanov, I. Babushkin, and R. M. Arkhipov, Phys. Rev. A 101, 013803 (2020)

  14. [14]

    Z. Bai, C. S. Adams, G. Huang, and W. Li, Phys. Rev. Lett. 125, 263605 (2020)

  15. [15]

    Jirauschek, M

    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)

  16. [16]

    Hughes, Phys

    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)

  17. [17]

    A. M. Alhasan, J. Fiutak and W. Miklaszewski, Z. Phys. B 88, 349 (1992)

  18. [18]

    Miklaszewski and J

    W. Miklaszewski and J. Fiutak, Z. Phys. B 93, 491 (1994)

  19. [19]

    M. J. Konopnicki and J. H. Eberly, Phys. Rev. A 24, 2567 (1981)

  20. [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)

  21. [21]

    Huang and C

    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)

  22. [22]

    F. T. Hioe and R. Grobe, Phys. Rev. Lett. 73, 2559 (1994)

  23. [23]

    M. D. Crisp, Phys. Rev. A 5, 1365 (1972)

  24. [24]

    S´ egard, B

    B. S´ egard, B. Macke, J. Zemmouri, and W. Sergent, Ann. Phys. (Paris) 15, 167 (1990)

  25. [25]

    de Lamare, Ph

    J. de Lamare, Ph. Kupecek, and M. Comte, Opt. Com- mun. 95, 305 (1993)

  26. [26]

    M. D. Crisp, Phys. Rev. A 1, 1604 (1970); ibid. 2, 2172 15 (1970)

  27. [27]

    J. E. Rothenberg, D. Grischkowsky, and A. C. Balant, Phys. Rev. Lett. 53, 552 (1984)

  28. [28]

    Avenel, E

    O. Avenel, E. Varoquaux, and G. A. Williams, Phys. Rev. Lett. 53, 2058 (1984)

  29. [29]

    E. M. Pessina, B. S´ egard, and B. Macke, Opt. Comm. 81, 397 (1991)

  30. [30]

    W. R. LeFew, S. Venakides, and D. J. Gauthier, Phys. Rev. A 79, 063842 (2009)

  31. [31]

    D. Wei, J. F. Chen, M. M. T. Loy, G. K. L. Wong, and S. Du, Phys. Rev. Lett. 103, 093602 (2009)

  32. [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)

  33. [33]

    Horovitz and N

    B. Horovitz and N. Rosenberg, Phys. Rev. A 26, 2799 (1982)

  34. [34]

    Macke and B

    B. Macke and B. S´ egard, Phys. Rev. A81, 015803 (2010)

  35. [35]

    D. J. Kaup, Phys. Rev. A 16, 704 (1977)

  36. [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)

  37. [37]

    B. S. Cartwright, Optical transients in atomic vapours, Durham theses, Durham University (2022), http://etheses.dur.ac.uk/14532/

  38. [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. [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. [40]

    V. V. Kozlov and E. B. Kozlova, Opt. Spektrosk. 107, 139 (2009) [Opt. Spectrosc. 107, 129 (2009)]

  41. [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)]

  42. [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]

  43. [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)]

  44. [44]

    V. V. Kozlov and E. E. Fradkin, Pis’ma Zh. Eksp. Teor. Fiz. 68, 359 (1998) [JETP Lett. 68, 383 (1998)]

  45. [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)]

  46. [46]

    V. V. Kozlov, P. G. Polynkin, and M. O. Scully, Phys. Rev. A 59, 3060 (1999)

  47. [47]

    Paspalakis, N

    E. Paspalakis, N. J. Kylstra, and P. L. Knight, Phys. Rev. A 61, 045802 (2000)

  48. [48]

    V. V. Kozlov and E. B. Kozlova, Opt. Spektrosk. 108, 824 (2010) [Opt. Spectrosc. 108, 780 (2010)]

  49. [49]

    O. M. Fedotova, O. K. Khasanov, G. A. Rusetsky, J. Degert, and E. Freysz, Phys. Rev. A 90, 053843 (2014)

  50. [50]

    R. M. Potvliege and S. A. Wrathmall, Comput. Phys. Commun. 306, 109374 (2025)

  51. [51]

    D. A. Steck, Rubidium 85 D line data, available online at http://steck.us/alkalidata

  52. [52]

    H. P. Grieneisen, J. Goldhar, N. A. Kurnit, and A. Javan, Appl. Phys. Lett. 21, 559 (1972)

  53. [53]

    S. M. Hamadani, J. Goldhar, N. A. Kurnit, and A. Javan, Appl. Phys. Lett. 25, 160 (1974)

  54. [54]

    Matusovsky, B

    M. Matusovsky, B. Vaynberg, and M. Rosenbluh, J. Opt. Soc. Am. B 13, 1994 (1996)

  55. [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. [56]

    Rahman, Phys

    A. Rahman, Phys. Rev. A 60, 4187 (1999); A. Rahman and J. H. Eberly, Opt. Express 4, 133 (1999)

  57. [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. [58]

    E.g., C. R. Higgins and I. G. Hughes, J. Phys. B: At. Mol. Opt. Phys. 54, 165403 (2021)

  59. [59]

    T. Y. Abi-Salloum, Phys. Rev. A 81, 053836 (2010)

  60. [60]

    S. Khan, V. Bharti, and V. Natarajan, Phys. Lett. A 380, 4100 (2016)

  61. [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)

  62. [62]

    Siddons, J

    P. Siddons, J. Phys. B: At. Mol. Opt. Phys. 47, 093001 (2014)

  63. [63]

    Chang, W.-C

    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