pith. sign in

arxiv: 2508.11417 · v2 · submitted 2025-08-15 · ⚛️ physics.atom-ph · physics.optics

Wavelength dependence of laser pulse filamentation in the close spectral vicinity of atomic resonances

Gabor Demeter This is my paper

Pith reviewed 2026-05-18 23:29 UTC · model grok-4.3

classification ⚛️ physics.atom-ph physics.optics
keywords laser filamentationrubidium vaporatomic resonancesplasma channelsself-focusinganomalous dispersionmultiphoton ionizationwavelength dependence
0
0 comments X

The pith

Laser pulses in rubidium vapor self-focus sharply below the D2 resonance but weaken above it.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper investigates how the central wavelength of intense laser pulses affects their nonlinear propagation and plasma channel formation in rubidium vapor near atomic resonances. Simulations of pulse propagation show that wavelengths below the 780 nm D2 resonance lead to strong self-focusing and sharp plasma boundaries, much like the resonant case, whereas wavelengths above 780 nm produce weaker self-focusing and more diffuse boundaries. The authors attribute this to an interplay of anomalous dispersion near the resonance, additional resonant transitions between excited rubidium states in the 740-780 nm range, and wavelength-dependent multiphoton ionization rates. A reader would care because it explains observed differences in plasma channel quality and suggests ways to control filamentation through wavelength choice.

Core claim

Numerical solutions of the propagation equations reveal that pulse wavelengths below the rubidium D2 resonance exhibit strong self-focusing and sharp plasma boundaries similar to the resonant 780 nm case, while wavelengths above 780 nm yield gradually weaker self-focusing and increasingly diffuse plasma boundaries. This behavior stems from the combined influence of anomalous dispersion around atomic resonances, resonant transitions between excited states of Rb in the 740-780 nm range, and wavelength-dependent multiphoton ionization rates.

What carries the argument

The numerical solutions of the propagation equations incorporating anomalous dispersion, resonant excited-state transitions, and multiphoton ionization rates.

If this is right

  • Plasma channels are more homogeneous with sharper boundaries for sub-resonant wavelengths.
  • Weaker self-focusing occurs for super-resonant wavelengths leading to diffuse boundaries.
  • The model accounts for the distinct behaviors seen at 780 nm versus 810 nm in experiments.
  • The interplay of dispersion, resonances, and ionization determines the filamentation outcome near atomic lines.

Where Pith is reading between the lines

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

  • This could allow wavelength tuning to optimize plasma channel homogeneity for applications like laser guiding.
  • Similar wavelength dependencies might occur in other atomic vapors with nearby resonances, affecting filamentation there.
  • Direct measurements of ionization rates at these wavelengths could confirm the model's key factors.

Load-bearing premise

The numerical solutions of the propagation equations capture all relevant nonlinear and resonant effects without missing important contributions from the atomic medium or requiring unstated adjustments to model parameters.

What would settle it

Observing no significant change in plasma boundary sharpness or self-focusing strength when tuning across the resonance in rubidium, or finding that multiphoton ionization rates are independent of wavelength in this range, would falsify the proposed interplay.

Figures

Figures reproduced from arXiv: 2508.11417 by Gabor Demeter.

Figure 1
Figure 1. Figure 1: FIG. 1. a) Electronic levels of the rubidium atom that are [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Transmitted pulse properties at [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Contour plots of the ionization probability [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. a) Plasma channel radius [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Contour plots of the ionization probability [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Contour plots of the radiant fluence [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Contour plots of a) the plasma radius [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Plots of the real part of the atomic polarization [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
read the original abstract

We investigate the propagation and nonlinear self-focusing of TW power laser pulses that create 10-m-scale, highly homogeneous plasma channels in rubidium vapor. Using computational solutions of the relevant propagation equations, we study the effects of the ionizing pulse central wavelength in relation to the resonance frequencies of atomic rubidium. Recent experiments show that pulse propagation and plasma channel creation is distinctly different for 780 nm laser pulses (resonant with the rubidium $D_2$ line) and 810 nm laser pulses. We study pulse propagation in a $\pm$30 nm range around the $D_2$ resonance and find that the results are distinctly different when tuning to sub-resonant wavelengths from those obtained for super-resonant wavelengths. For pulse wavelengths below the resonance the behavior is similar to the resonant case, characterized by strong self-focusing and a sharp plasma boundary. Pulse wavelengths above 780 nm on the other hand yield gradually weaker self-focusing and an increasingly diffuse plasma boundary. Our results suggest that the observed behavior can be attributed to an interplay between multiple factors: anomalous dispersion around atomic resonances, resonant transitions between excited states of Rb lying in the 740-780 nm range and wavelength-dependent multiphoton ionization rates.

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 / 2 minor

Summary. The paper investigates the wavelength dependence of TW-power laser pulse propagation and filamentation in rubidium vapor near the D2 resonance at 780 nm. Using numerical solutions of the propagation equations over a ±30 nm range, it reports that sub-resonant wavelengths produce strong self-focusing and sharp plasma boundaries similar to the resonant case at 780 nm, while super-resonant wavelengths yield weaker self-focusing and increasingly diffuse plasma boundaries. The authors attribute the sub- versus super-resonant asymmetry to the combined effects of anomalous dispersion, resonant transitions between excited Rb states in the 740–780 nm range, and wavelength-dependent multiphoton ionization rates.

Significance. If the numerical results are robust, the work offers a plausible mechanistic account for experimentally observed differences in plasma-channel formation near atomic resonances. This could inform control of long, homogeneous plasma channels in high-intensity laser–atom interactions, with relevance to filamentation studies and applications requiring precise plasma boundaries.

major comments (2)
  1. [Computational model / propagation equations] The central attribution to resonant excited-state transitions (740–780 nm) and wavelength-dependent MPI rates is load-bearing, yet the manuscript provides no explicit description of how these are incorporated into the propagation equations or susceptibility model. If an effective two-level or generic Keldysh treatment is used without multi-level tracking of 5P–5D transitions, the reported asymmetry could arise from incomplete physics rather than the claimed interplay.
  2. [Numerical methods] No information is given on the numerical scheme, spatial/temporal grid resolution, convergence tests, or direct validation against the referenced experiments. Without these, it is impossible to assess whether the solutions reliably capture the claimed wavelength dependence of self-focusing and plasma-boundary sharpness.
minor comments (2)
  1. [Abstract] The abstract states that results are 'distinctly different' for sub- versus super-resonant tuning but does not quantify the wavelength range or the magnitude of the differences in self-focusing or boundary sharpness.
  2. [References] Ensure that all cited experimental works on Rb filamentation are referenced with complete bibliographic details.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments, which will help strengthen the presentation of our results. We address each major comment below and will revise the manuscript to provide the requested details.

read point-by-point responses

  1. Referee: [Computational model / propagation equations] The central attribution to resonant excited-state transitions (740–780 nm) and wavelength-dependent MPI rates is load-bearing, yet the manuscript provides no explicit description of how these are incorporated into the propagation equations or susceptibility model. If an effective two-level or generic Keldysh treatment is used without multi-level tracking of 5P–5D transitions, the reported asymmetry could arise from incomplete physics rather than the claimed interplay.

    Authors: We agree that an explicit description of the model is essential for supporting the central claims. In the revised manuscript we will add a dedicated subsection that details the multi-level susceptibility model. This will specify how the resonant transitions between excited Rb states (including 5P–5D lines in the 740–780 nm window) enter the nonlinear polarization term, how the anomalous dispersion is computed from the full set of nearby resonances, and how the wavelength-dependent multiphoton ionization rate is obtained from a resonance-adjusted Keldysh-type expression. These additions will make clear that the reported sub- versus super-resonant asymmetry arises from the combined action of dispersion, excited-state resonances, and ionization rather than from an oversimplified two-level treatment. revision: yes

  2. Referee: [Numerical methods] No information is given on the numerical scheme, spatial/temporal grid resolution, convergence tests, or direct validation against the referenced experiments. Without these, it is impossible to assess whether the solutions reliably capture the claimed wavelength dependence of self-focusing and plasma-boundary sharpness.

    Authors: We acknowledge that the current version omits these technical details. The revised manuscript will include a new subsection describing the numerical scheme (split-step Fourier propagation with adaptive step-size control), the spatial and temporal grid resolutions employed, and the results of systematic convergence tests performed by successively refining the grids until the self-focusing distance and plasma-boundary sharpness stabilize to within a few percent. We will also add a brief comparison of simulated plasma-channel lengths and boundary profiles with the qualitative features reported in the referenced experiments, noting the level of agreement that can be expected given the available experimental data. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation relies on standard numerical propagation modeling

full rationale

The paper computes solutions to propagation equations for TW laser pulses in Rb vapor and attributes sub- vs super-resonant differences in self-focusing and plasma boundary sharpness to the interplay of anomalous dispersion, excited-state resonances (740-780 nm), and wavelength-dependent MPI rates. These attributions rest on known atomic physics and nonlinear optics rather than any self-definitional loop, fitted parameter renamed as prediction, or load-bearing self-citation. No equations or results in the manuscript reduce by construction to the paper's own inputs; the numerical model is externally falsifiable against independent measurements of dispersion, transition strengths, and ionization rates. The central claim therefore remains independent of the present work's fitted values or definitions.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based on abstract only; no explicit free parameters, axioms, or invented entities are stated. The work relies on standard nonlinear propagation models and known atomic physics effects.

pith-pipeline@v0.9.0 · 5742 in / 1001 out tokens · 47032 ms · 2026-05-18T23:29:27.604419+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

39 extracted references · 39 canonical work pages

  1. [1]

    Couairon and A

    A. Couairon and A. Mysyrowicz, Physics Reports441, 47– (2007)

    work page 2007

  2. [2]

    Berg´ e, S

    L. Berg´ e, S. Skupin, R. Nuter, J. Kasparian, and J.-P. Wolf, Reports on Progress in Physics70, 1633– (2007)

    work page 2007

  3. [3]

    V. P. Kandidov, S A Shlenov, and O. G. Kosareva, Quan- tum Electronics39, 205– (2009)

    work page 2009

  4. [4]

    P. B. Corkum and F. Krausz, Nature Physics3, 381 (2007)

    work page 2007

  5. [5]

    Kasparian and J.-P

    J. Kasparian and J.-P. Wolf, Opt. Express16, 466 (2008)

    work page 2008

  6. [6]

    P. Chen, J. Dawson, R. Huff, and T. Katsouleas, Phys- ical Review Letters54, 693 (1985)

    work page 1985

  7. [7]

    Blumenfeld, C

    I. Blumenfeld, C. Clayton, F. Decker, M. Hogan, C. Huang, C. Joshi, T. Katsouleas, W. Lu, W. Mori, P. Muggli,et al., Nature445, 741 (2007)

    work page 2007

  8. [8]

    Joshi and T

    C. Joshi and T. Katsouleas, Physics Today56, 47 (2003)

    work page 2003

  9. [9]

    Esarey, C

    E. Esarey, C. B. Schroeder, and W. P. Leemans, Rev. Mod. Phys.81, 1229 (2009)

    work page 2009

  10. [10]

    Leemans and E

    W. Leemans and E. Esarey, Physics Today62, 44 (2009)

    work page 2009

  11. [11]

    Tajima, X

    T. Tajima, X. Q. Yan, and T. Ebisuzaki, Reviews of Modern Plasma Physics4, 7 (2020)

    work page 2020

  12. [12]

    Caldwell, K

    A. Caldwell, K. Lotov, A. Pukhov, and F. Simon, Nature Physics5, 363– (2009)

    work page 2009

  13. [13]

    E. Gschwendtneret al., Nuclear Instruments and Meth- ods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 829, 76– (2016), 2nd European Advanced Accelerator Concepts Workshop - EAAC 2015

    work page 2016

  14. [14]

    Adliet al.(AWAKE Collaboration), Nature561, 363 (2018)

    E. Adliet al.(AWAKE Collaboration), Nature561, 363 (2018)

    work page 2018

  15. [15]

    K. V. Lotov, A. Pukhov, and A. Caldwell, Physics of Plasmas20, 013102 (2013)

    work page 2013

  16. [16]

    ¨Oz and P

    E. ¨Oz and P. Muggli, Nuclear Instruments and Meth- ods in Physics Research Section A: Accelerators, Spec- trometers, Detectors and Associated Equipment740, 197 (2014), proceedings of the first European Advanced Ac- celerator Concepts Workshop 2013

    work page 2014

  17. [17]

    Collaboration,AWAKE Design Report: A Proton- Driven Plasma Wakefield Acceleration Experiment at CERN, Tech

    A. Collaboration,AWAKE Design Report: A Proton- Driven Plasma Wakefield Acceleration Experiment at CERN, Tech. Rep. CERN-SPSC-2013-013, SPSC-TDR- 003 (CERN, 2013)

    work page 2013

  18. [18]

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

    work page 1971

  19. [19]

    de Lamare, M

    J. de Lamare, M. Comte, and P. Kupecek, Phys. Rev. A 50, 3366– (1994)

    work page 1994

  20. [20]

    Delagnes and M

    J. Delagnes and M. Bouchene, Optics Communications 281, 5824– (2008)

    work page 2008

  21. [21]

    Plyushchev, R

    G. Plyushchev, R. Kersevan, A. Petrenko, and P. Muggli, Journal of Physics D: Applied Physics51, 025203 (2017)

    work page 2017

  22. [22]

    Demeter, Phys

    G. Demeter, Phys. Rev. A99, 063423 (2019)

    work page 2019

  23. [23]

    Demeter, J

    G. Demeter, J. T. Moody, M. A. Kedves, B. R´ aczkevi, M. Aladi, A.-M. Bachmann, F. Batsch, F. Braunm¨ uller, G. P. Djotyan, V. Fedosseev, F. Friebel, S. Gess- ner, E. Granados, E. Guran, M. H¨ uther, V. Lee, M. Martyanov, P. Muggli, E. ¨Oz, H. Panuganti, L. Verra, and G. Z. Della Porta, Phys. Rev. A104, 033506 (2021)

    work page 2021

  24. [24]

    Demeter, J

    G. Demeter, J. Moody, M. Kedves, F. Batsch, M. Berga- maschi, V. Fedosseev, E. Granados, P. Muggli, H. Panu- ganti, and G. Zevi Della Porta, Optics & Laser Technol- ogy168, 109921 (2024)

    work page 2024

  25. [25]

    Gschwendtner, K

    E. Gschwendtner, K. Lotov, P. Muggli, M. Wing, et al.(AWAKE Collaboration), Symmetry14(2022), 10.3390/sym14081680

    work page doi:10.3390/sym14081680 2022

  26. [26]

    Houard, P

    A. Houard, P. Walsch, T. Produit,et al., Nature Photon- ics (2023), https://doi.org/10.1038/s41566-022-01139-z

    work page doi:10.1038/s41566-022-01139-z 2023

  27. [27]

    Brabec and F

    T. Brabec and F. Krausz, Phys. Rev. Lett.78, 3282– (1997)

    work page 1997

  28. [28]

    Couairon, E

    A. Couairon, E. Brambilla, T. Corti, D. Majus, O. de J. Ram´ ırez-G´ ongora, and M. Kolesik, The Eu- ropean Physical Journal Special Topics199, 5– (2011)

    work page 2011

  29. [29]

    Boyd,Nonlinear Optics(Elsevier Science, 2003)

    R. Boyd,Nonlinear Optics(Elsevier Science, 2003)

    work page 2003

  30. [30]

    A. M. Perelomov, V. S. Popov, and M. V. Terent’ev, Soviet Physics JETP23, 924 (1966)

    work page 1966

  31. [31]

    A. M. Perelomov, V. S. Popov, and M. V. Terent’ev, Soviet Physics JETP24, 207 (1967)

    work page 1967

  32. [32]

    A. M. Perelomov and V. S. Popov, Soviet Physics JETP 25, 336 (1967)

    work page 1967

  33. [33]

    B´ ejot,Theoretical and experimental investigations of ultrashort laser filamentation in gases, Ph.D

    P. B´ ejot,Theoretical and experimental investigations of ultrashort laser filamentation in gases, Ph.D. the- sis (2008), , Universit´ e de Gen´ eve 2008 https://archive- ouverte.unige.ch/unige:4568

    work page 2008

  34. [34]

    B. C. Duncan, V. Sanchez-Villicana, P. L. Gould, and H. R. Sadeghpour, Phys. Rev. A63, 043411 (2001)

    work page 2001

  35. [35]

    D. A. Steck,Rubidium 85 D Line Data(available online at http://steck.us/alkalidata, (revision 2.1.2, 12 August 2009))

    work page 2009

  36. [36]

    M. S. Safronova, C. J. Williams, and C. W. Clark, Phys. Rev. A69, 022509 (2004)

    work page 2004

  37. [37]

    Kramida, Y

    A. Kramida, Y. Ralchenko, J. Reader, and N. A. Team, NIST Atomic Spectra Database(National Institute of Standards and Technology, Gaithersburg, MD., 2018)

    work page 2018

  38. [38]

    B. E. A. Saleh and M. C. Teich,Fundamentals of Pho- tonics(John Wiley and Sons, Inc., 2019)

    work page 2019

  39. [39]

    H´ eder, E

    M. H´ eder, E. Rig´ o, D. Medgyesi, R. Lovas, S. Tenczer, F. T¨ or¨ ok, A. Farkas, M. Em˝ odi, J. Kadlecsik, G. Mez˝ o, ´A. Pint´ er, and P. Kacsuk, Inform´ aci´ os T´ arsadalom22, 128 (2022)

    work page 2022