@bookBudker:2013, title=Optical Magnetometry, author=edited by D. Budker and D. F. Jackson Kimball, year=2013, publisher=Cambridge University Press, Cambridge, England

@articleBudker/N:2007, title=Optical magnetometry, author=Budker, Dmitry and Romalis, Michael, journal=Nat. Phys., volume=3, pages=227–234, year=2007, doi=10.1038/nphys566

@articleJohnson/PMB:2013, doi = 10.1088/0031-9155/58/17/6065, url = https://dx.doi.org/10.1088/0031-9155/58/17/6065, year = 2013, month = aug, publisher = IOP Publishing, volume = 58, number = 17, pages = 6065, author = Cort N Johnson and P D D Schwindt and M Weisend, title = Multi-sensor magnetoencephalography with atomic magnetometers, journal = Phys. Med. Biol.,

@articleVasilakis/PRL:2009, title = Limits on New Long Range Nuclear Spin-Dependent Forces Set with a 𝐊3Hesuperscript3𝐊He\mathbf{K}\mathrm{\text{$-$}}^{3}\mathrm{He}bold_K - start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_He Comagnetometer, author = Vasilakis, G. and Brown, J. M. and Kornack, T. W. and Romalis, M. V., journal = Phys. Rev. Lett., volume = 103, issue = 26, pages = 261801, numpages = 4, year = 2009, month = Dec, publisher = American Physical Society, doi = 10.1103/PhysRevLett.103.261801, url = https://link.aps.org/doi/10.1103/PhysRevLett.103.261801

@articleDang/APL:2010, author = Dang, H. B. and Maloof, A. C. and Romalis, M. V., title = ”Ultrahigh sensitivity magnetic field and magnetization measurements with an atomic magnetometer”, journal = Appl. Phys. Lett., volume = 97, number = 15, pages = 151110, year = 2010, month = 10, issn = 0003-6951, doi = 10.1063/1.3491215, url = https://doi.org/10.1063/1.3491215

@articleShah/NP:2007, title=Subpicotesla atomic magnetometry with a microfabricated vapour cell, author=Shah, V. and Knappe, S. and Schwindt, P. D. D. and Kitching, J., journal=Nat. Photonics, volume=1, pages=649, year=2007, doi=10.1038/nphoton.2007.201

@articleRubinsztein-Dunlop/JO:2017, doi = 10.1088/2040-8978/19/1/013001, url = https://dx.doi.org/10.1088/2040-8978/19/1/013001, year = 2017, month = nov, publisher = IOP Publishing, volume = 19, number = 1, pages = 013001, author = Halina Rubinsztein-Dunlop and Andrew Forbes and M V Berry and M R Dennis and David L Andrews and Masud Mansuripur and Cornelia Denz and Christina Alpmann and Peter Banzer and Thomas Bauer and Ebrahim Karimi and Lorenzo Marrucci and Miles Padgett and Monika Ritsch-Marte and Natalia M Litchinitser and Nicholas P Bigelow and C Rosales-Guzmán and A Belmonte and J P Torres and Tyler W Neely and Mark Baker and Reuven Gordon and Alexander B Stilgoe and Jacquiline Romero and Andrew G White and Robert Fickler and Alan E Willner and Guodong Xie and Benjamin McMorran and Andrew M Weiner, title = Roadmap on structured light, journal = J. Opt.,

@bookGbur:2017, title=Singular optics, author=Gbur, Gregory J, year=2017, publisher=CRC press

@articleCastellucci/PRL:2021, title = Atomic Compass: Detecting 3D Magnetic Field Alignment with Vector Vortex Light, author = Castellucci, Francesco and Clark, Thomas W. and Selyem, Adam and Wang, Jinwen and Franke-Arnold, Sonja, journal = Phys. Rev. Lett., volume = 127, issue = 23, pages = 233202, numpages = 6, year = 2021, month = Nov, publisher = American Physical Society, doi = 10.1103/PhysRevLett.127.233202, url = https://link.aps.org/doi/10.1103/PhysRevLett.127.233202

@articleQiu/PR:2021, author = Shuwei Qiu and Jinwen Wang and Francesco Castellucci and Mingtao Cao and Shougang Zhang and Thomas W. Clark and Sonja Franke-Arnold and Hong Gao and Fuli Li, journal = Photon. Res., keywords = CCD cameras; Cylindrical vector beams; Nitrogen vacancy centers; Optical fields; Tunable diode lasers; Vector beams, number = 12, pages = 2325–2331, publisher = Optica Publishing Group, title = Visualization of magnetic fields with cylindrical vector beams in a warm atomic vapor, volume = 9, month = Dec, year = 2021, url = https://opg.optica.org/prj/abstract.cfm?URI=prj-9-12-2325, doi = 10.1364/PRJ.418522,

@articleSavukov/PRL:2005, title = Tunable Atomic Magnetometer for Detection of Radio-Frequency Magnetic Fields, author = Savukov, I. M. and Seltzer, S. J. and Romalis, M. V. and Sauer, K. L., journal = Phys. Rev. Lett., volume = 95, issue = 6, pages = 063004, numpages = 4, year = 2005, month = Aug, publisher = American Physical Society, doi = 10.1103/PhysRevLett.95.063004, url = https://link.aps.org/doi/10.1103/PhysRevLett.95.063004

@articleLedbetter/PRA:2007, title = Detection of radio-frequency magnetic fields using nonlinear magneto-optical rotation, author = Ledbetter, M. P. and Acosta, V. M. and Rochester, S. M. and Budker, D. and Pustelny, S. and Yashchuk, V. V., journal = Phys. Rev. A, volume = 75, issue = 2, pages = 023405, numpages = 6, year = 2007, month = Feb, publisher = American Physical Society, doi = 10.1103/PhysRevA.75.023405, url = https://link.aps.org/doi/10.1103/PhysRevA.75.023405

@bookBransden:2003, title=Physics of Atoms and Molecules, author=Bransden, Brian Harold and Joachain, Charles Jean, year=2003, publisher=Prentice Hall, Harlow, England

@bookBlum:2012, title=Density Matrix Theory and Applications, author=Blum, Karl, year=2012, publisher=Springer, Berlin

@bookAuzinsh:2010, title=Optically Polarized Atoms: Understanding Light-Atom Interactions, author=Auzinsh, Marcis and Budker, Dmitry and Rochester, Simon M, year=2010, publisher=Oxford University, Oxford

@articleWense:2020, title=The theory of direct laser excitation of nuclear transitions, author=von der Wense, Lars and Bilous, Pavlo V and Seiferle, Benedict and Stellmer, Simon and Weitenberg, Johannes and Thirolf, Peter G and Pálffy, Adriana and Kazakov, Georgy, journal=Eur. Phys. J. A, volume=56, pages=176, year=2020, publisher=Springer, url = https://doi.org/10.1140/epja/s10050-020-00177-x

@articleTremblay/PRA:1990, title = Optical pumping with two finite linewidth lasers, author = Tremblay, P. and Jacques, C., journal = Phys. Rev. A, volume = 41, issue = 9, pages = 4989–4999, numpages = 0, year = 1990, month = May, publisher = American Physical Society, doi = 10.1103/PhysRevA.41.4989, url = https://link.aps.org/doi/10.1103/PhysRevA.41.4989

@articleSchmidt/PRA:2024, title = Atomic photoexcitation as a tool for probing purity of twisted light modes, author = Schmidt, R. P. and Ramakrishna, S. and Peshkov, A. A. and Huntemann, N. and Peik, E. and Fritzsche, S. and Surzhykov, A., journal = Phys. Rev. A, volume = 109, issue = 3, pages = 033103, numpages = 11, year = 2024, month = Mar, publisher = American Physical Society, doi = 10.1103/PhysRevA.109.033103, url = https://link.aps.org/doi/10.1103/PhysRevA.109.033103

@articleMatula/JPB:2013, doi = 10.1088/0953-4075/46/20/205002, url = https://dx.doi.org/10.1088/0953-4075/46/20/205002, year = 2013, month = oct, publisher = IOP Publishing, volume = 46, number = 20, pages = 205002, author = O Matula and A G Hayrapetyan and V G Serbo and A Surzhykov and S Fritzsche, title = Atomic ionization of hydrogen-like ions by twisted photons: angular distribution of emitted electrons, journal = J. Phys. B,

@articleKnyazev/PU:2018, doi = 10.3367/UFNe.2018.02.038306, url = https://dx.doi.org/10.3367/UFNe.2018.02.038306, year = 2018, month = may, publisher = Uspekhi Fizicheskikh Nauk, Russian Academy of Sciences and IOP Publishing, volume = 61, pages = 449, author = B A Knyazev and V G Serbo, title = Beams of photons with nonzero projections of orbital angular momenta: new results, journal = Phys.-Usp.,

@articleSchulz/PRA:2020, title = Generalized excitation of atomic multipole transitions by twisted light modes, author = Schulz, S. A.-L. and Peshkov, A. A. and Müller, R. A. and Lange, R. and Huntemann, N. and Tamm, Chr. and Peik, E. and Surzhykov, A., journal = Phys. Rev. A, volume = 102, issue = 1, pages = 012812, numpages = 10, year = 2020, month = Jul, publisher = American Physical Society, doi = 10.1103/PhysRevA.102.012812, url = https://link.aps.org/doi/10.1103/PhysRevA.102.012812

@bookJohnson:2007, title=Atomic Structure Theory, author=Johnson, Walter R, year=2007, publisher=Springer, New York

@articleSolyanik-Gorgone/JOSAB:2019, author = Maria Solyanik-Gorgone and Andrei Afanasev and Carl E. Carlson and Christian T. Schmiegelow and Ferdinand Schmidt-Kaler, journal = J. Opt. Soc. Am. B, keywords = Circular polarization; Laser beams; Light beams; Light matter interactions; Optical vortices; Orbital angular momentum multiplexing, number = 3, pages = 565–574, publisher = Optica Publishing Group, title = Excitation of E1-forbidden atomic transitions with electric, magnetic, or mixed multipolarity in light fields carrying orbital and spin angular momentum

Invited𝐼𝑛𝑣𝑖𝑡𝑒𝑑Inviteditalic_I italic_n italic_v italic_i italic_t italic_e italic_d

, volume = 36, month = Mar, year = 2019, url = https://opg.optica.org/josab/abstract.cfm?URI=josab-36-3-565, doi = 10.1364/JOSAB.36.000565

@articlePeshkov/AdP:2023, author = Peshkov, Anton A. and Jordan, Elena and Kromrey, Markus and Mehta, Karan K. and Mehlstäubler, Tanja E. and Surzhykov, Andrey, title = Excitation of Forbidden Electronic Transitions in Atoms by Hermite–Gaussian Modes, journal = Ann. Phys., volume = 535, number = 9, pages = 2300204, doi = https://doi.org/10.1002/andp.202300204, url = https://onlinelibrary.wiley.com/doi/abs/10.1002/andp.202300204, year = 2023

@bookRose:1957, title=Elementary Theory of Angular Momentum, author=Rose, Morris Edgar, year=1957, publisher=John Wiley & Sons, New York

@articleFritzsche/CPC:2019, title=A fresh computational approach to atomic structures, processes and cascades, author=Stephan Fritzsche, journal=Comput. Phys. Commun., volume=240, number=1, pages=1-14, year=2019, doi = 10.1016/j.cpc.2019.01.012,

@articleMaldonado/OE:2024, title=Sensitivity of a vector atomic magnetometer based on electromagnetically induced transparency, author=Mario Gonzalez Maldonado and Owen Rollins and Alex Toyryla and James A. McKelvy and Andrey Matsko and Isaac Fan and Yang Li and Ying-Ju Wang and John Kitching and Irina Novikova and Eugeniy E. Mikhailov, journal=Opt. Express, volume=32, pages=25062-25073, year=2024, doi = 10.1364/OE.529276,

@articleWang/AVSQS:2020, title=Vectorial light–matter interaction: Exploring spatially structured complex light fields, author=Jinwen Wang and Francesco Castellucci and Sonja Franke-Arnold, journal=AVS Quantum Sci., volume=2, pages=031702, year=2020, doi = 10.1116/5.0016007,

@articleWang/PRL:2024, title=Measuring the Optical Concurrence of Vector Beams with an Atomic-State Interferometer, author=Jinwen Wang and Sphinx J. Svensson and Thomas W. Clark and Yun Chen and Mustafa A. Al Khafaji and Hong Gao and Niclas Westerberg and Sonja Franke-Arnold, journal=Phys. Rev. Lett., volume=132, pages=193803, year=2024, doi = 10.1103/PhysRevLett.132.193803,

@articleVolz/PHysScripta:1996, title=Precision lifetime measurements on alkali atoms and on helium by beam–gas–laser spectroscopy, author=U. Volz and H. Schmoranzer, journal=Physica Scripta, volume=T65, pages=48-56,, year=1996, doi = 10.1088/0031-8949/1996/T65/007,

thanks: These two authors contributed equally to this work.thanks: These two authors contributed equally to this work.

Interaction of vector light beams with atoms
exposed to a time-dependent magnetic field

S. Ramakrishna [email protected] Helmholtz-Institut Jena, D-07743 Jena, Germany GSI Helmholtzzentrum für Schwerionenforschung GmbH, D-64291 Darmstadt, Germany Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universität Jena, D-07743 Jena, Germany    R. P. Schmidt [email protected] Physikalisch-Technische Bundesanstalt, D-38116 Braunschweig, Germany Institut für Mathematische Physik, Technische Universität Braunschweig, D-38106 Braunschweig, Germany    A. A. Peshkov Physikalisch-Technische Bundesanstalt, D-38116 Braunschweig, Germany Institut für Mathematische Physik, Technische Universität Braunschweig, D-38106 Braunschweig, Germany    S. Franke-Arnold School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, United Kingdom    A. Surzhykov Physikalisch-Technische Bundesanstalt, D-38116 Braunschweig, Germany Institut für Mathematische Physik, Technische Universität Braunschweig, D-38106 Braunschweig, Germany Laboratory for Emerging Nanometrology Braunschweig, D-38106 Braunschweig, Germany    S. Fritzsche Helmholtz-Institut Jena, D-07743 Jena, Germany GSI Helmholtzzentrum für Schwerionenforschung GmbH, D-64291 Darmstadt, Germany Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universität Jena, D-07743 Jena, Germany
(July 25, 2024)
Abstract

During recent years interest has been rising for applications of vector light beams towards magnetic field sensing. In particular, a series of experiments were performed to extract information about properties of static magnetic fields from absorption profiles of light passing through an atomic gas target. In the present work, we propose an extension to this method for oscillating magnetic fields. To investigate this scenario, we carried out theoretical analysis based on the time-dependent density matrix theory. We found that absorption profiles, even when averaged over typical observation times, are indeed sensitive to both strength and frequency of the time-dependent field, thus opening the prospect for a powerful diagnostic technique. To illustrate this sensitivity, we performed detailed calculations for the 5sS1/225𝑠superscriptsubscript𝑆1225s\;\,{}^{2}S_{1/2}5 italic_s start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT (F=1𝐹1F=1italic_F = 1) -- 5pP3/225𝑝superscriptsubscript𝑃3225p\;\,{}^{2}P_{3/2}5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT (F=0𝐹0F=0italic_F = 0) transition in rubidium atoms, subject to a superposition of an oscillating (test) and a static (reference) magnetic field.

I Introduction

In optical magnetometry, magnetic field properties are measured by observing changes in the optical properties of an atomic medium immersed in the field edited by D. Budker and Kimball (2013). The great advantage of this technology compared to superconducting magnetic field sensors is that it offers high sensitivity without requiring cryogenic temperatures Budker and Romalis (2007). Optical magnetometers are finding applications in a wide variety of fields including medicine Johnson et al. (2013), fundamental physics Vasilakis et al. (2009), and geophysics Dang et al. (2010). Significant progress has been made in miniaturizing these devices and improving their operating characteristics Shah et al. (2007).

In most optical magnetometry experiments, the polarization across the light beam is approximately uniform. Meanwhile, recent advances in optics made it possible to generate light fields with space-varying polarization Rubinsztein-Dunlop et al. (2017). The best known examples of such vector beams include radially and azimuthally polarized beams with an azimuthally varying linear polarization surrounding an optical vortex Gbur (2017). These spatially varying light polarizations can excite locally varying magnetization profiles in atoms Wang et al. (2020). It has been demonstrated recently in cold Castellucci et al. (2021) and warm Qiu et al. (2021); Maldonado et al. (2024) atomic vapors to test the magnetic field components transverse and along the optical axis simultaneously. These experiments determined static magnetic fields from absorption profiles of a vector beam after its passage through the rubidium vapor.

Beyond static magnetic fields, the detection of time-dependent magnetic fields is important, especially in the radio-frequency domain. Detection of fields in the kilohertz to gigahertz frequency range finds many applications, from radio communication to detection of nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) signals Savukov et al. (2005). An example of an atomic magnetometer for detection of radio-frequency magnetic fields is magnetometer Ledbetter et al. (2007) based on a nonlinear magneto-optical rotation. In this paper, we explore the possibility of detecting oscillating (radio-frequency) magnetic fields based on another effect, namely the dependence of the absorption profile of a vector beam, propagating through an atomic vapor, on the strength and frequency of the magnetic field. In addition to the light field and the magnetic field to be measured, here we need an additional reference static magnetic field with well-known direction and strength. As in previous works Savukov et al. (2005); Ledbetter et al. (2007), the static magnetic field is perpendicular to the oscillating one. While, however, the static field in Refs. Savukov et al. (2005); Ledbetter et al. (2007) was added to achieve the Zeeman resonance of atoms, it is here applied in order to trigger oscillations in the orientation of the resulting magnetic field relative to the light polarization.

Refer to caption
Figure 1: Proposed geometry of the experimental setup: The atomic target is subjected to the combination of a reference DC magnetic field 𝑩constsubscript𝑩const\bm{B}_{\mathrm{const}}bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT and an AC test field 𝑩AC=𝑩osccos(ωosct)subscript𝑩ACsubscript𝑩oscsubscript𝜔osc𝑡\bm{B}_{\mathrm{AC}}=\bm{B}_{\mathrm{osc}}\cos(\omega_{\mathrm{osc}}t)bold_italic_B start_POSTSUBSCRIPT roman_AC end_POSTSUBSCRIPT = bold_italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT roman_cos ( italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT italic_t ), and illuminated by a vector light beam propagating along z𝑧zitalic_z. The quantization axis is chosen along the time varying direction of the total magnetic field in the x𝑥xitalic_x-y𝑦yitalic_y plane.

To determine the effect of the applied magnetic field on the light absorption profile, we employ density matrix theory, whose basic formulas are briefly reviewed in Sec. II. We show, in particular, how to calculate the time evolution of the atomic density matrix from the transition amplitudes for vector beams and Zeeman shifts caused by the superposition of the reference and test magnetic fields. From this density matrix we obtain in Sec. III the light absorption profile illustrated for the case of the 87Rb D2222 line. We first investigate the time evolution of the absorption profiles and show that the petal-like absorption patterns rotate about the beam axis at a rate that depends on the applied AC magnetic field. In general, this effect can be used for magnetometry, but it requires a high time resolution of the detector. An easier characteristic to observe is the absorption profile averaged over the measurement time. We found that even this profile can be sensitive to both the strength and frequency of the test magnetic field, and this sensitivity is most pronounced for frequencies up to a hundred kHz and strengths up to several Gauss for the reference field of about one Gauss. Finally, Sec. IV provides a brief summary and outlook.

II Theory

II.1 Geometry of the process

In the present work we consider the interaction of atoms with a structured light beam propagating along the z𝑧zitalic_z-axis and having the frequency ω𝜔\omegaitalic_ω. We assume that the atoms are exposed to a combination of (reference) static and (test) oscillating magnetic fields. The static 𝑩const=Bconst𝒆xsubscript𝑩constsubscript𝐵constsubscript𝒆𝑥\bm{B}_{\mathrm{const}}=B_{\mathrm{const}}\bm{e}_{x}bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT bold_italic_e start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT is applied along x𝑥xitalic_x, while 𝑩AC=𝑩osccos(ωosct)=Bosccos(ωosct)𝒆ysubscript𝑩ACsubscript𝑩osccossubscript𝜔osc𝑡subscript𝐵osccossubscript𝜔osc𝑡subscript𝒆𝑦\bm{B}_{\mathrm{AC}}=\bm{B}_{\mathrm{osc}}\mathrm{cos}(\omega_{\mathrm{osc}}t)% =B_{\mathrm{osc}}\mathrm{cos}(\omega_{\mathrm{osc}}t)\bm{e}_{y}bold_italic_B start_POSTSUBSCRIPT roman_AC end_POSTSUBSCRIPT = bold_italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT roman_cos ( italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT italic_t ) = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT roman_cos ( italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT italic_t ) bold_italic_e start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT oscillates along the y𝑦yitalic_y-axis, as shown in Fig. 1. The resulting magnetic field 𝑩total(t)=𝑩const+𝑩osccos(ωosct)subscript𝑩total𝑡subscript𝑩constsubscript𝑩osccossubscript𝜔osc𝑡\bm{B}_{\mathrm{total}}(t)=\bm{B}_{\mathrm{const}}+\bm{B}_{\mathrm{osc}}% \mathrm{cos}(\omega_{\mathrm{osc}}t)bold_italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ) = bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT + bold_italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT roman_cos ( italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT italic_t ) has the strength

Btotal(t)=Bconst2+Bosc2cos2(ωosct)subscript𝐵total𝑡superscriptsubscript𝐵const2superscriptsubscript𝐵osc2superscriptcos2subscript𝜔osc𝑡\displaystyle B_{\mathrm{total}}(t)=\sqrt{B_{\mathrm{const}}^{2}+B_{\mathrm{% osc}}^{2}\mathrm{cos}^{2}(\omega_{\mathrm{osc}}t)}italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ) = square-root start_ARG italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT italic_t ) end_ARG (1)

oscillating between Bconstsubscript𝐵constB_{\mathrm{const}}italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT and Bconst2+Bosc2superscriptsubscript𝐵const2superscriptsubscript𝐵osc2\sqrt{B_{\mathrm{const}}^{2}+B_{\mathrm{osc}}^{2}}square-root start_ARG italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG, while its direction is characterized by the angle

ϕB(t)=arctan(Bosccos(ωosct)/Bconst)subscriptitalic-ϕ𝐵𝑡subscript𝐵osccossubscript𝜔osc𝑡subscript𝐵const\displaystyle\phi_{B}(t)=\arctan(B_{\mathrm{osc}}\mathrm{cos}(\omega_{\mathrm{% osc}}t)/B_{\mathrm{const}})italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) = roman_arctan ( italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT roman_cos ( italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT italic_t ) / italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT ) (2)

varying in the range arctan(Bosc/Bconst)ϕB+arctan(Bosc/Bconst)arctansubscript𝐵oscsubscript𝐵constsubscriptitalic-ϕ𝐵arctansubscript𝐵oscsubscript𝐵const-\mathrm{arctan}(B_{\mathrm{osc}}/B_{\mathrm{const}})\leq\phi_{B}\leq+\mathrm{% arctan}(B_{\mathrm{osc}}/B_{\mathrm{const}})- roman_arctan ( italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT / italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT ) ≤ italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ≤ + roman_arctan ( italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT / italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT ). A typical variation of Btotal(t)subscript𝐵total𝑡B_{\mathrm{total}}(t)italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ) and ϕB(t)subscriptitalic-ϕ𝐵𝑡\phi_{B}(t)italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) with t𝑡titalic_t is shown in Fig. 2. We take fairly close values of Bconstsubscript𝐵constB_{\mathrm{const}}italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT and Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT so that the direction ϕB(t)subscriptitalic-ϕ𝐵𝑡\phi_{B}(t)italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) of the resulting magnetic field oscillates with a sufficiently large amplitude.

By design, the field 𝑩total(t)subscript𝑩total𝑡\bm{B}_{\mathrm{total}}(t)bold_italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ) is always perpendicular to the light propagation direction and oscillates around 𝑩constsubscript𝑩const\bm{B}_{\mathrm{const}}bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT. For such a complex geometry, particular attention should be paid to the choice of the quantization axis of the entire system. In our work we take the quantization axis to be along the resulting magnetic field 𝑩total(t)subscript𝑩total𝑡\bm{B}_{\mathrm{total}}(t)bold_italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ), because this choice simplifies the calculation of the Zeeman splitting Bransden and Joachain (2003). On the other hand, the description of the coupling between atoms and photons becomes somewhat more complicated because the direction of 𝑩total(t)subscript𝑩total𝑡\bm{B}_{\mathrm{total}}(t)bold_italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ), and hence of the quantization axis, is not stationary in the reference frame of the light beam. Of course, the observables should not depend on the specific choice of the quantization axis. To see this, we additionally performed calculations with the quantization axis along the static magnetic field 𝑩constsubscript𝑩const\bm{B}_{\mathrm{const}}bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT. Both choices of the quantization axis lead to the same result.

Refer to caption
Figure 2: Time dependence of the strength (a) and orientiation (b) of the magnetic field 𝑩total(t)=𝑩const+𝑩osccos(ωosct)subscript𝑩total𝑡subscript𝑩constsubscript𝑩osccossubscript𝜔osc𝑡\bm{B}_{\mathrm{total}}(t)=\bm{B}_{\mathrm{const}}+\bm{B}_{\mathrm{osc}}% \mathrm{cos}(\omega_{\mathrm{osc}}t)bold_italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ) = bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT + bold_italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT roman_cos ( italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT italic_t ). Results are shown for Bosc=subscript𝐵oscabsentB_{\mathrm{osc}}=italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G (black solid lines), Bosc=subscript𝐵oscabsentB_{\mathrm{osc}}=italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 G (blue dotted lines), Bosc=subscript𝐵oscabsentB_{\mathrm{osc}}=italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 5 G (red dashed lines), and Bosc=subscript𝐵oscabsentB_{\mathrm{osc}}=italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 10 G (green dash-dotted lines). For all curves, ωosc=2π×\omega_{\mathrm{osc}}=2\pi\timesitalic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 100 kHz and Bconst=subscript𝐵constabsentB_{\mathrm{const}}=italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G.

II.2 Vector light beams

In optical physics, light beams are usually described in terms of the electric field. Such choice allows for an easier treatment of light-matter interaction in the dipole approximation. In this work we assume that the electric field of the incident vector light beam in the paraxial regime has the form:

𝐄(vec)(𝐫,t)=superscript𝐄vec𝐫𝑡absent\displaystyle\mathbf{E}^{\mathrm{(vec)}}(\mathbf{r},t)=bold_E start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT ( bold_r , italic_t ) = E0[cos(2ϕ)𝒆x+sin(2ϕ)𝒆y]subscript𝐸0delimited-[]cos2italic-ϕsubscript𝒆𝑥sin2italic-ϕsubscript𝒆𝑦\displaystyle\,E_{0}\left[\mathrm{cos}\left(2\phi\right)\bm{e}_{x}+\mathrm{sin% }\left(2\phi\right)\bm{e}_{y}\right]italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ roman_cos ( 2 italic_ϕ ) bold_italic_e start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + roman_sin ( 2 italic_ϕ ) bold_italic_e start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ]
×J2(ϰr)eikzzeiωt,absentsubscript𝐽2italic-ϰsubscript𝑟bottomsuperscript𝑒𝑖subscript𝑘𝑧𝑧superscript𝑒𝑖𝜔𝑡\displaystyle\times J_{2}(\varkappa r_{\bot})e^{ik_{z}z}e^{-i\omega t},× italic_J start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_ϰ italic_r start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_z end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT , (3)

where E0subscript𝐸0E_{0}italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is a constant amplitude, ω𝜔\omegaitalic_ω is its frequency, rsubscript𝑟perpendicular-tor_{\perp}italic_r start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT, ϕitalic-ϕ\phiitalic_ϕ, and z𝑧zitalic_z are cylindrical coordinates, kzsubscript𝑘𝑧k_{z}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT and ϰitalic-ϰ\varkappaitalic_ϰ are longitudinal and transverse components of the linear momentum, respectively, and Jn(ϰr)subscript𝐽𝑛italic-ϰsubscript𝑟perpendicular-toJ_{n}(\varkappa r_{\perp})italic_J start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_ϰ italic_r start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT ) is the Bessel function. In order to develop a general and relativistic theory as well as to simplify the discussion, it is more convenient to describe the incident radiation in terms of the vector potential. For example, in order to obtain electric field (3), one has to start from the vector potential

𝑨(vec)(𝒓,t)=12superscript𝑨vec𝒓𝑡12\displaystyle\bm{A}^{\mathrm{(vec)}}(\bm{r},t)=\frac{1}{\sqrt{2}}bold_italic_A start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT ( bold_italic_r , italic_t ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG [𝑨mγ,λ=+1(B)(𝒓,t)𝑨mγ,λ=1(B)(𝒓,t)]delimited-[]subscriptsuperscript𝑨Bsubscript𝑚𝛾𝜆1𝒓𝑡subscriptsuperscript𝑨Bsubscript𝑚𝛾𝜆1𝒓𝑡\displaystyle\left[\bm{A}^{\mathrm{(B)}}_{-m_{\gamma},\,\lambda=+1}(\bm{r},t)-% \bm{A}^{\mathrm{(B)}}_{m_{\gamma},\,\lambda=-1}(\bm{r},t)\right]\,[ bold_italic_A start_POSTSUPERSCRIPT ( roman_B ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT , italic_λ = + 1 end_POSTSUBSCRIPT ( bold_italic_r , italic_t ) - bold_italic_A start_POSTSUPERSCRIPT ( roman_B ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT , italic_λ = - 1 end_POSTSUBSCRIPT ( bold_italic_r , italic_t ) ] (4)

which is a linear combination of two Bessel beams. Since these wave solutions with an annular intensity structure have been frequently discussed in the past Knyazev and Serbo (2018); Matula et al. (2013); Schulz et al. (2020), we may restrict ourselves to a rather short account of basic formulas. The Bessel beam is characterized by the well-defined helicity λ𝜆\lambdaitalic_λ and the projection mγsubscript𝑚𝛾m_{\gamma}italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT of the total angular momentum upon the propagation direction. Moreover, its longitudinal momentum kzsubscript𝑘𝑧k_{z}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT and absolute value of the transverse momentum ϰ=|𝒌|italic-ϰsubscript𝒌bottom\varkappa=|\bm{k}_{\bot}|italic_ϰ = | bold_italic_k start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT | are also fixed, see Ref. Matula et al. (2013) for further details. The vector potential for the Bessel beam can be written as

𝑨mγ,λ(B)(𝒓)=A0aϰmγ(𝒌)𝒆𝒌λei𝒌𝒓eiωtd2𝒌(2π)2,subscriptsuperscript𝑨Bsubscript𝑚𝛾𝜆𝒓subscript𝐴0subscript𝑎italic-ϰsubscript𝑚𝛾subscript𝒌bottomsubscript𝒆𝒌𝜆superscript𝑒𝑖𝒌𝒓superscript𝑒𝑖𝜔𝑡superscript𝑑2subscript𝒌bottomsuperscript2𝜋2\displaystyle\bm{A}^{\mathrm{(B)}}_{m_{\gamma},\lambda}(\bm{r})=A_{0}\int a_{% \varkappa m_{\gamma}}(\bm{k}_{\bot})\;\bm{e}_{\bm{k}\lambda}e^{i\bm{k}\cdot\bm% {r}}e^{-i\omega t}\;\frac{d^{2}\bm{k}_{\bot}}{(2\pi)^{2}},bold_italic_A start_POSTSUPERSCRIPT ( roman_B ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT , italic_λ end_POSTSUBSCRIPT ( bold_italic_r ) = italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∫ italic_a start_POSTSUBSCRIPT italic_ϰ italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_italic_k start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT ) bold_italic_e start_POSTSUBSCRIPT bold_italic_k italic_λ end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i bold_italic_k ⋅ bold_italic_r end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT bold_italic_k start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (5)

where A0=E0ωsubscript𝐴0subscript𝐸0𝜔A_{0}=-E_{0}\omegaitalic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ω and aϰmγ(𝒌)subscript𝑎italic-ϰsubscript𝑚𝛾subscript𝒌bottoma_{\varkappa m_{\gamma}}(\bm{k}_{\bot})italic_a start_POSTSUBSCRIPT italic_ϰ italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_italic_k start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT ) is a weight function given by:

aϰmγ(𝒌)=2πϰ(i)mγeimγϕkδ(kϰ).subscript𝑎italic-ϰsubscript𝑚𝛾subscript𝒌bottom2𝜋italic-ϰsuperscript𝑖subscript𝑚𝛾superscript𝑒𝑖subscript𝑚𝛾subscriptitalic-ϕ𝑘𝛿subscript𝑘bottomitalic-ϰ\displaystyle a_{\varkappa m_{\gamma}}(\bm{k}_{\bot})=\frac{2\pi}{\varkappa}(-% i)^{m_{\gamma}}e^{im_{\gamma}\phi_{k}}\delta(k_{\bot}-\varkappa).italic_a start_POSTSUBSCRIPT italic_ϰ italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_italic_k start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT ) = divide start_ARG 2 italic_π end_ARG start_ARG italic_ϰ end_ARG ( - italic_i ) start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_δ ( italic_k start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT - italic_ϰ ) . (6)

It follows from these expressions that the Bessel beam can be seen as a superposition of plane waves whose wave vectors 𝒌=(𝒌,kz)𝒌subscript𝒌bottomsubscript𝑘𝑧\bm{k}=(\bm{k}_{\bot},k_{z})bold_italic_k = ( bold_italic_k start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) are uniformly distributed upon the surface of a cone with a polar opening angle θk=arctan(ϰ/kz)subscript𝜃𝑘arctanitalic-ϰsubscript𝑘𝑧\theta_{k}=\mathrm{arctan}(\varkappa/k_{z})italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = roman_arctan ( italic_ϰ / italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ).

Refer to caption
Figure 3: The intensity profile (displayed as a density plot) and polarization pattern (black arrows) of the vector beam (7) in the paraxial approximation with θk=0.11°subscript𝜃𝑘0.11°\theta_{k}=0.11\degreeitalic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = 0.11 ° and ω=2π×384𝜔2𝜋384\omega=2\pi\times 384italic_ω = 2 italic_π × 384 THz.

The light field (4) is often called the vector beam because its polarization pattern is varying across the profile Rubinsztein-Dunlop et al. (2017); Gbur (2017). In this work we consider the vector beam consisting of Bessel modes (5) with the total angular momentum projections mγ=±1subscript𝑚𝛾plus-or-minus1m_{\gamma}=\pm 1italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT = ± 1. For arbitrary values of the opening angle θksubscript𝜃𝑘\theta_{k}italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT, the vector beam exhibits spatially dependent polarization along all three axes. Atomic physics experiments, however, commonly use vector beams produced in the paraxial regime where the transverse momentum of the photon is much smaller than the longitudinal momentum, ϰkzmuch-less-thanitalic-ϰsubscript𝑘𝑧\varkappa\ll k_{z}italic_ϰ ≪ italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT Matula et al. (2013). In this regime θksubscript𝜃𝑘\theta_{k}italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT is small and the vector potential (4) can be considerably simplified to:

𝐀(vec)(𝐫,t)superscript𝐀vec𝐫𝑡absent\displaystyle\mathbf{A}^{\mathrm{(vec)}}(\mathbf{r},t)\approxbold_A start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT ( bold_r , italic_t ) ≈ A0(i)[cos(2ϕ)𝒆x+sin(2ϕ)𝒆y]subscript𝐴0𝑖delimited-[]cos2italic-ϕsubscript𝒆𝑥sin2italic-ϕsubscript𝒆𝑦\displaystyle A_{0}(-i)\left[\mathrm{cos}\left(2\phi\right)\bm{e}_{x}+\mathrm{% sin}\left(2\phi\right)\bm{e}_{y}\right]italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( - italic_i ) [ roman_cos ( 2 italic_ϕ ) bold_italic_e start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + roman_sin ( 2 italic_ϕ ) bold_italic_e start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ]
×J2(ϰr)eikzzeiωt.absentsubscript𝐽2italic-ϰsubscript𝑟bottomsuperscript𝑒𝑖subscript𝑘𝑧𝑧superscript𝑒𝑖𝜔𝑡\displaystyle\times J_{2}(\varkappa r_{\bot})e^{ik_{z}z}e^{-i\omega t}.× italic_J start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_ϰ italic_r start_POSTSUBSCRIPT ⊥ end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_z end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT . (7)

Let us remark here that the vector beam in the paraxial approximation (7) has azimuthally varying linear polarization that lies in the transverse x𝑥xitalic_x-y𝑦yitalic_y plane Qiu et al. (2021). The corresponding spatial field distribution is depicted in Fig. 3. From the vector potential (7) and by using the standard relation 𝑬=t𝑨𝑬subscript𝑡𝑨\bm{E}=-\partial_{t}\bm{A}bold_italic_E = - ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT bold_italic_A, we can finally obtain the electric field (3).

II.3 Transition amplitudes

Having discussed the vector potential for the incident light, we are ready now to examine its interaction with atoms. In particular, we will question the laser-induced transition between ground |αgFgMgketsubscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔\left|\alpha_{g}F_{g}M_{g}\right>| italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ and excited |αeFeMeketsubscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒\left|\alpha_{e}F_{e}M_{e}\right>| italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩ atomic states whose properties can be traced back to the first-order matrix element

Veg(vec)subscriptsuperscript𝑉vec𝑒𝑔\displaystyle\,V^{\mathrm{(vec)}}_{eg}italic_V start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT
=\displaystyle== ecαeFeMe|q𝜶q𝑨(vec)(𝒓q𝒃)|αgFgMg𝑒𝑐quantum-operator-productsubscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒subscript𝑞subscript𝜶𝑞superscript𝑨vecsubscript𝒓𝑞𝒃subscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔\displaystyle\,ec\;\left\langle\alpha_{e}F_{e}M_{e}\left|\sum_{q}\bm{\alpha}_{% q}\cdot\bm{A}^{\mathrm{(vec)}}(\bm{r}_{q}-\bm{b})\right|\alpha_{g}F_{g}M_{g}\right\rangleitalic_e italic_c ⟨ italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT | ∑ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT bold_italic_α start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ⋅ bold_italic_A start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT ( bold_italic_r start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT - bold_italic_b ) | italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ (8)

with 𝑭=𝑰+𝑱𝑭𝑰𝑱\bm{F}=\bm{I}+\bm{J}bold_italic_F = bold_italic_I + bold_italic_J, where 𝑰𝑰\bm{I}bold_italic_I and 𝑱𝑱\bm{J}bold_italic_J are the nuclear and electron angular momenta, respectively, M𝑀Mitalic_M is the projection of 𝑭𝑭\bm{F}bold_italic_F on the quantization axis, and α𝛼\alphaitalic_α denotes all additional quantum numbers required to specify the state uniquely. Moreover, q𝑞qitalic_q runs over all electrons in a target atom and 𝜶qsubscript𝜶𝑞\bm{\alpha}_{q}bold_italic_α start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT denotes the vector of Dirac matrices for the q𝑞qitalic_qth particle Johnson (2007). In Eq. (8), we have introduced the impact parameter 𝒃=(bcosϕb,bsinϕb,0)𝒃𝑏subscriptitalic-ϕ𝑏𝑏subscriptitalic-ϕ𝑏0\bm{b}=(b\cos\phi_{b},b\sin\phi_{b},0)bold_italic_b = ( italic_b roman_cos italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , italic_b roman_sin italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , 0 ) to specify the position of the atom within the wave front Knyazev and Serbo (2018). In particular, b=0𝑏0b=0italic_b = 0 corresponds to an atom located on the vortex line of a light beam. The impact parameter plays an important role in the analysis due to the complex spatial structure of the vector beam, see Fig. 3.

Similar to the vector potential (4), the transition amplitude Veg(vec)subscriptsuperscript𝑉vec𝑒𝑔V^{\mathrm{(vec)}}_{eg}italic_V start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT can be expressed in terms of its Bessel counterparts as

Veg(vec)superscriptsubscript𝑉𝑒𝑔vec\displaystyle\,V_{eg}^{\mathrm{(vec)}}italic_V start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT
=\displaystyle== 12[Veg(B)(mγ,λ=+1)Veg(B)(mγ,λ=1)].12delimited-[]superscriptsubscript𝑉𝑒𝑔Bsubscript𝑚𝛾𝜆1subscriptsuperscript𝑉B𝑒𝑔subscript𝑚𝛾𝜆1\displaystyle\,\frac{1}{\sqrt{2}}\,\left[V_{eg}^{(\mathrm{B})}(-m_{\gamma},% \lambda=+1)-V^{(\mathrm{B})}_{eg}(m_{\gamma},\lambda=-1)\right].divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG [ italic_V start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_B ) end_POSTSUPERSCRIPT ( - italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT , italic_λ = + 1 ) - italic_V start_POSTSUPERSCRIPT ( roman_B ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT , italic_λ = - 1 ) ] . (9)

The evaluation of the transition amplitude for Bessel beams (4) has already been discussed in detail in Refs. Schulz et al. (2020); Solyanik-Gorgone et al. (2019); Peshkov et al. (2023). For the geometry shown in Fig. 1, the final form of this amplitude is

Veg(B)(mγ,λ)=A0ec2πMiL+M[L,Fg]1/2(iλ)psuperscriptsubscript𝑉𝑒𝑔Bsubscript𝑚𝛾𝜆subscript𝐴0𝑒𝑐2𝜋subscript𝑀superscript𝑖𝐿𝑀superscript𝐿subscript𝐹𝑔12superscript𝑖𝜆𝑝\displaystyle V_{eg}^{(\mathrm{B})}(m_{\gamma},\lambda)=A_{0}ec\sqrt{2\pi}\;% \sum_{M}i^{L+M}\;[L,F_{g}]^{1/2}\;(i\lambda)^{p}\,italic_V start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_B ) end_POSTSUPERSCRIPT ( italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT , italic_λ ) = italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e italic_c square-root start_ARG 2 italic_π end_ARG ∑ start_POSTSUBSCRIPT italic_M end_POSTSUBSCRIPT italic_i start_POSTSUPERSCRIPT italic_L + italic_M end_POSTSUPERSCRIPT [ italic_L , italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT ( italic_i italic_λ ) start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT
×(1)mγei(mγM)ϕbJmγM(ϰb)dM,λL(θk)absentsuperscript1subscript𝑚𝛾superscript𝑒𝑖subscript𝑚𝛾𝑀subscriptitalic-ϕ𝑏subscript𝐽subscript𝑚𝛾𝑀italic-ϰ𝑏subscriptsuperscript𝑑𝐿𝑀𝜆subscript𝜃𝑘\displaystyle\;\;\;\;\;\times(-1)^{m_{\gamma}}\;e^{i(m_{\gamma}-M)\phi_{b}}\;J% _{m_{\gamma}-M}(\varkappa b)\,d^{L}_{M,\lambda}(\theta_{k})\,× ( - 1 ) start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT - italic_M ) italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT - italic_M end_POSTSUBSCRIPT ( italic_ϰ italic_b ) italic_d start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_M , italic_λ end_POSTSUBSCRIPT ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT )
×DMeMg,ML(π,π/2,πϕB(t))absentsubscriptsuperscript𝐷𝐿subscript𝑀𝑒subscript𝑀𝑔𝑀𝜋𝜋2𝜋subscriptitalic-ϕ𝐵𝑡\displaystyle\;\;\;\;\;\times D^{L}_{M_{e}-M_{g},M}(\pi,\pi/2,\pi-\phi_{B}(t))\,× italic_D start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_M end_POSTSUBSCRIPT ( italic_π , italic_π / 2 , italic_π - italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) )
×FgMgLMeMg|FeMe(1)I+Fg+L+Jeabsentinner-productsubscript𝐹𝑔subscript𝑀𝑔𝐿subscript𝑀𝑒subscript𝑀𝑔subscript𝐹𝑒subscript𝑀𝑒superscript1𝐼subscript𝐹𝑔𝐿subscript𝐽𝑒\displaystyle\;\;\;\;\;\times\langle F_{g}M_{g}LM_{e}-M_{g}|F_{e}M_{e}\rangle% \;(-1)^{I+F_{g}+L+J_{e}}\,× ⟨ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_L italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT | italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩ ( - 1 ) start_POSTSUPERSCRIPT italic_I + italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT + italic_L + italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUPERSCRIPT
×{FeFgLJgJeI}αeJeHγ(pL)αgJg,absentmatrixsubscript𝐹𝑒subscript𝐹𝑔𝐿subscript𝐽𝑔subscript𝐽𝑒𝐼delimited-⟨⟩subscript𝛼𝑒subscript𝐽𝑒normsubscript𝐻𝛾𝑝𝐿subscript𝛼𝑔subscript𝐽𝑔\displaystyle\;\;\;\;\;\times\begin{Bmatrix}F_{e}&F_{g}&L\\ J_{g}&J_{e}&I\end{Bmatrix}\langle\alpha_{e}J_{e}||H_{\gamma}(pL)||\alpha_{g}J_% {g}\rangle,× { start_ARG start_ROW start_CELL italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_L end_CELL end_ROW start_ROW start_CELL italic_J start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_I end_CELL end_ROW end_ARG } ⟨ italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT | | italic_H start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_p italic_L ) | | italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ , (10)

where we have used the notation αeJeHγ(pL)αgJgdelimited-⟨⟩subscript𝛼𝑒subscript𝐽𝑒normsubscript𝐻𝛾𝑝𝐿subscript𝛼𝑔subscript𝐽𝑔\langle\alpha_{e}J_{e}||H_{\gamma}(pL)||\alpha_{g}J_{g}\rangle⟨ italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT | | italic_H start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_p italic_L ) | | italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ to denote the reduced matrix element for magnetic (p=0𝑝0p=0italic_p = 0) and electric (p=1𝑝1p=1italic_p = 1) transitions. Furthermore, the Euler angles as the arguments of the Wigner D𝐷Ditalic_D-functions dM,λL(θk)subscriptsuperscript𝑑𝐿𝑀𝜆subscript𝜃𝑘d^{L}_{M,\lambda}(\theta_{k})italic_d start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_M , italic_λ end_POSTSUBSCRIPT ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) and DMeMg,ML(π,π/2,πϕB(t))subscriptsuperscript𝐷𝐿subscript𝑀𝑒subscript𝑀𝑔𝑀𝜋𝜋2𝜋subscriptitalic-ϕ𝐵𝑡D^{L}_{M_{e}-M_{g},M}(\pi,\pi/2,\pi-\phi_{B}(t))italic_D start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_M end_POSTSUBSCRIPT ( italic_π , italic_π / 2 , italic_π - italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) ) in Eq. (10) characterize the rotation from the atomic frame with the quantization axis along the magnetic field 𝑩total(t)subscript𝑩total𝑡\bm{B}_{\mathrm{total}}(t)bold_italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ) to the photon frame with the quantization axis along the wave vector 𝒌𝒌\bm{k}bold_italic_k Rose (1957). The time-dependent angle πϕB(t)𝜋subscriptitalic-ϕ𝐵𝑡\pi-\phi_{B}(t)italic_π - italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) refers to the oscillation of the polarization vector of the incident light at the position of the atom in its (time-dependent) reference frame.

II.4 Density-matrix formalism

Due to the presence of both the incident radiation and the time-dependent magnetic field, the populations of atomic ground and excited states can vary with time. To investigate time dependence of atomic level populations, it is convenient to use the time-dependent density matrix theory Blum (2012). In this approach, the state of a system is represented by the density operator ρ^(t)^𝜌𝑡\hat{\rho}(t)over^ start_ARG italic_ρ end_ARG ( italic_t ) satisfying the Liouville-von Neumann equation:

ddtρ^(t)=i[H^(t),ρ^(t)]+R^(t).𝑑𝑑𝑡^𝜌𝑡𝑖Planck-constant-over-2-pi^𝐻𝑡^𝜌𝑡^𝑅𝑡\displaystyle\frac{d}{dt}\hat{\rho}(t)=-\frac{i}{\hbar}[\hat{H}(t),\hat{\rho}(% t)]+\hat{R}(t).divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG over^ start_ARG italic_ρ end_ARG ( italic_t ) = - divide start_ARG italic_i end_ARG start_ARG roman_ℏ end_ARG [ over^ start_ARG italic_H end_ARG ( italic_t ) , over^ start_ARG italic_ρ end_ARG ( italic_t ) ] + over^ start_ARG italic_R end_ARG ( italic_t ) . (11)

Here H^(t)^𝐻𝑡\hat{H}(t)over^ start_ARG italic_H end_ARG ( italic_t ) is the total Hamiltonian of the atom in the presence of external fields, and R^(t)^𝑅𝑡\hat{R}(t)over^ start_ARG italic_R end_ARG ( italic_t ) is introduced to take into account phenomenologically spontaneous decay Auzinsh et al. (2010). We consider transitions between the Zeeman sublevels of the ground |αgFgMgketsubscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔|\alpha_{g}F_{g}M_{g}\rangle| italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ to those of the excited state |αeFeMeketsubscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒|\alpha_{e}F_{e}M_{e}\rangle| italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩, described by a density matrix of size (2Fg+2Fe+2)×(2Fg+2Fe+2)2subscript𝐹𝑔2subscript𝐹𝑒22subscript𝐹𝑔2subscript𝐹𝑒2(2F_{g}+2F_{e}+2)\times(2F_{g}+2F_{e}+2)( 2 italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT + 2 italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT + 2 ) × ( 2 italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT + 2 italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT + 2 ). In this basis, we can write the elements of the density matrix as:

ρgg(t)=subscript𝜌𝑔superscript𝑔𝑡absent\displaystyle\rho_{gg^{\prime}}(t)=italic_ρ start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = αgFgMg|ρ^(t)|αgFgMg,quantum-operator-productsubscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔^𝜌𝑡subscript𝛼𝑔subscript𝐹𝑔subscriptsuperscript𝑀𝑔\displaystyle\langle\alpha_{g}F_{g}M_{g}|\hat{\rho}(t)|\alpha_{g}F_{g}M^{% \prime}_{g}\rangle,⟨ italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT | over^ start_ARG italic_ρ end_ARG ( italic_t ) | italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ , (12a)
ρee(t)=subscript𝜌𝑒superscript𝑒𝑡absent\displaystyle\rho_{ee^{\prime}}(t)=italic_ρ start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = αeFeMe|ρ^(t)|αeFeMe,quantum-operator-productsubscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒^𝜌𝑡subscript𝛼𝑒subscript𝐹𝑒subscriptsuperscript𝑀𝑒\displaystyle\langle\alpha_{e}F_{e}M_{e}|\hat{\rho}(t)|\alpha_{e}F_{e}M^{% \prime}_{e}\rangle,⟨ italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT | over^ start_ARG italic_ρ end_ARG ( italic_t ) | italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩ , (12b)
ρge(t)=subscript𝜌𝑔𝑒𝑡absent\displaystyle\rho_{ge}(t)=italic_ρ start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) = αgFgMg|ρ^(t)|αeFeMe,quantum-operator-productsubscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔^𝜌𝑡subscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒\displaystyle\langle\alpha_{g}F_{g}M_{g}|\hat{\rho}(t)|\alpha_{e}F_{e}M_{e}\rangle,⟨ italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT | over^ start_ARG italic_ρ end_ARG ( italic_t ) | italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩ , (12c)
ρeg(t)=subscript𝜌𝑒𝑔𝑡absent\displaystyle\rho_{eg}(t)=italic_ρ start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) = αeFeMe|ρ^(t)|αgFgMg.quantum-operator-productsubscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒^𝜌𝑡subscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔\displaystyle\langle\alpha_{e}F_{e}M_{e}|\hat{\rho}(t)|\alpha_{g}F_{g}M_{g}\rangle.⟨ italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT | over^ start_ARG italic_ρ end_ARG ( italic_t ) | italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ . (12d)

In Eqs. (12), the diagonal elements ρgg(t)subscript𝜌𝑔𝑔𝑡\rho_{gg}(t)italic_ρ start_POSTSUBSCRIPT italic_g italic_g end_POSTSUBSCRIPT ( italic_t ) and ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) are the probabilities of finding an atom in the Zeeman substates |αgFgMgketsubscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔|\alpha_{g}F_{g}M_{g}\rangle| italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ and |αeFeMeketsubscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒|\alpha_{e}F_{e}M_{e}\rangle| italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩, whereas the off-diagonal elements describe the coherence between them.

In its matrix form the Liouville-von Neumann equation (11) represents a system of coupled differential equations for the evolution of the density matrix elements ρgg(t)subscript𝜌𝑔superscript𝑔𝑡\rho_{gg^{\prime}}(t)italic_ρ start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ), ρee(t)subscript𝜌𝑒superscript𝑒𝑡\rho_{ee^{\prime}}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ), ρge(t)subscript𝜌𝑔𝑒𝑡\rho_{ge}(t)italic_ρ start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ), and ρeg(t)subscript𝜌𝑒𝑔𝑡\rho_{eg}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ). To solve these equations, we introduce ρ~gg(t)=ρgg(t)subscript~𝜌𝑔superscript𝑔𝑡subscript𝜌𝑔superscript𝑔𝑡\tilde{\rho}_{gg^{\prime}}(t)=\rho_{gg^{\prime}}(t)over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = italic_ρ start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ), ρ~ee(t)=ρee(t)subscript~𝜌𝑒superscript𝑒𝑡subscript𝜌𝑒superscript𝑒𝑡\tilde{\rho}_{ee^{\prime}}(t)=\rho_{ee^{\prime}}(t)over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = italic_ρ start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ), ρ~ge(t)=ρge(t)eiωtsubscript~𝜌𝑔𝑒𝑡subscript𝜌𝑔𝑒𝑡superscript𝑒𝑖𝜔𝑡\tilde{\rho}_{ge}(t)=\rho_{ge}(t)e^{-i\omega t}over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) = italic_ρ start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) italic_e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT, ρ~eg(t)=ρeg(t)eiωtsubscript~𝜌𝑒𝑔𝑡subscript𝜌𝑒𝑔𝑡superscript𝑒𝑖𝜔𝑡\tilde{\rho}_{eg}(t)=\rho_{eg}(t)e^{i\omega t}over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) = italic_ρ start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) italic_e start_POSTSUPERSCRIPT italic_i italic_ω italic_t end_POSTSUPERSCRIPT and employ the rotating-wave approximation, which is valid when ω𝜔\omegaitalic_ω is sufficiently close to resonance Auzinsh et al. (2010); von der Wense et al. (2020). This approximation allows us to eliminate the fast-oscillating terms proportional to e±2iωtsuperscript𝑒plus-or-minus2𝑖𝜔𝑡e^{\pm 2i\omega t}italic_e start_POSTSUPERSCRIPT ± 2 italic_i italic_ω italic_t end_POSTSUPERSCRIPT, so that we can rewrite the Liouville-von Neumann equation as

ddtρ~gg(t)=𝑑𝑑𝑡subscript~𝜌𝑔superscript𝑔𝑡absent\displaystyle\allowdisplaybreaks\frac{d}{dt}\tilde{\rho}_{gg^{\prime}}(t)=divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = iΩg(L)(t)[MgMg]ρ~gg(t)i2[MeVegρ~eg(t)MeVegρ~ge(t)]+Rgg(t),𝑖superscriptsubscriptΩ𝑔L𝑡delimited-[]subscript𝑀𝑔subscriptsuperscript𝑀𝑔subscript~𝜌𝑔superscript𝑔𝑡𝑖2Planck-constant-over-2-pidelimited-[]subscriptsubscript𝑀𝑒superscriptsubscript𝑉𝑒𝑔subscript~𝜌𝑒superscript𝑔𝑡subscriptsubscript𝑀𝑒subscript𝑉𝑒superscript𝑔subscript~𝜌𝑔𝑒𝑡subscript𝑅𝑔superscript𝑔𝑡\displaystyle-i\Omega_{g}^{\mathrm{(L)}}(t)\left[M_{g}-M^{\prime}_{g}\right]% \tilde{\rho}_{gg^{\prime}}(t)\;-\frac{i}{2\hbar}\left[\sum_{M_{e}}V_{eg}^{*}\,% \tilde{\rho}_{eg^{\prime}}(t)-\sum_{M_{e}}V_{eg^{\prime}}\,\tilde{\rho}_{ge}(t% )\right]+R_{gg^{\prime}}(t),- italic_i roman_Ω start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ) [ italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT - italic_M start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ] over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) - divide start_ARG italic_i end_ARG start_ARG 2 roman_ℏ end_ARG [ ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) - ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) ] + italic_R start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) , (13a)
ddtρ~ee(t)=𝑑𝑑𝑡subscript~𝜌𝑒superscript𝑒𝑡absent\displaystyle\frac{d}{dt}\tilde{\rho}_{ee^{\prime}}(t)=divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = iΩe(L)(t)[MeMe]ρ~ee(t)i2[MgVegρ~ge(t)MgVegρ~eg(t)]+Ree(t),𝑖superscriptsubscriptΩ𝑒L𝑡delimited-[]subscript𝑀𝑒subscriptsuperscript𝑀𝑒subscript~𝜌𝑒superscript𝑒𝑡𝑖2Planck-constant-over-2-pidelimited-[]subscriptsubscript𝑀𝑔subscript𝑉𝑒𝑔subscript~𝜌𝑔superscript𝑒𝑡subscriptsubscript𝑀𝑔superscriptsubscript𝑉superscript𝑒𝑔subscript~𝜌𝑒𝑔𝑡subscript𝑅𝑒superscript𝑒𝑡\displaystyle\,-i\Omega_{e}^{\mathrm{(L)}}(t)\left[M_{e}-M^{\prime}_{e}\right]% \tilde{\rho}_{ee^{\prime}}(t)\;-\frac{i}{2\hbar}\left[\sum_{M_{g}}V_{eg}\,% \tilde{\rho}_{ge^{\prime}}(t)-\sum_{M_{g}}V_{e^{\prime}g}^{*}\,\tilde{\rho}_{% eg}(t)\right]+R_{ee^{\prime}}(t),- italic_i roman_Ω start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ) [ italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - italic_M start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ] over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) - divide start_ARG italic_i end_ARG start_ARG 2 roman_ℏ end_ARG [ ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) - ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) ] + italic_R start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) , (13b)
ddtρ~ge(t)=𝑑𝑑𝑡subscript~𝜌𝑔𝑒𝑡absent\displaystyle\frac{d}{dt}\tilde{\rho}_{ge}(t)=divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) = iΔρ~ge(t)+i[Ωe(L)(t)MeΩg(L)(t)Mg]ρ~ge(t)i2[MeVegρ~ee(t)MgVegρ~gg(t)]+Rge(t),𝑖Δsubscript~𝜌𝑔𝑒𝑡𝑖delimited-[]superscriptsubscriptΩ𝑒L𝑡subscript𝑀𝑒superscriptsubscriptΩ𝑔L𝑡subscript𝑀𝑔subscript~𝜌𝑔𝑒𝑡𝑖2Planck-constant-over-2-pidelimited-[]subscriptsuperscriptsubscript𝑀𝑒superscriptsubscript𝑉superscript𝑒𝑔subscript~𝜌superscript𝑒𝑒𝑡subscriptsuperscriptsubscript𝑀𝑔superscriptsubscript𝑉𝑒superscript𝑔subscript~𝜌𝑔superscript𝑔𝑡subscript𝑅𝑔𝑒𝑡\displaystyle\;-i\Delta\tilde{\rho}_{ge}(t)+i\left[\Omega_{e}^{\mathrm{(L)}}(t% )M_{e}-\Omega_{g}^{\mathrm{(L)}}(t)M_{g}\right]\tilde{\rho}_{ge}(t)\;-\frac{i}% {2\hbar}\left[\sum_{M_{e}^{\prime}}V_{e^{\prime}g}^{*}\,\tilde{\rho}_{e^{% \prime}e}(t)-\sum_{M_{g}^{\prime}}V_{eg^{\prime}}^{*}\,\tilde{\rho}_{gg^{% \prime}}(t)\right]+R_{ge}(t),- italic_i roman_Δ over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) + italic_i [ roman_Ω start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ) italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - roman_Ω start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ) italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ] over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) - divide start_ARG italic_i end_ARG start_ARG 2 roman_ℏ end_ARG [ ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_e end_POSTSUBSCRIPT ( italic_t ) - ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) ] + italic_R start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) , (13c)
ddtρ~eg(t)=𝑑𝑑𝑡subscript~𝜌𝑒𝑔𝑡absent\displaystyle\frac{d}{dt}\tilde{\rho}_{eg}(t)=divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) = iΔρ~eg(t)i[Ωe(L)(t)MeΩg(L)(t)Mg]ρ~eg(t)i2[MgVegρ~gg(t)MeVegρ~ee(t)]+Reg(t),𝑖Δsubscript~𝜌𝑒𝑔𝑡𝑖delimited-[]superscriptsubscriptΩ𝑒L𝑡subscript𝑀𝑒superscriptsubscriptΩ𝑔L𝑡subscript𝑀𝑔subscript~𝜌𝑒𝑔𝑡𝑖2Planck-constant-over-2-pidelimited-[]subscriptsuperscriptsubscript𝑀𝑔subscript𝑉𝑒superscript𝑔subscript~𝜌superscript𝑔𝑔𝑡subscriptsuperscriptsubscript𝑀𝑒subscript𝑉superscript𝑒𝑔subscript~𝜌𝑒superscript𝑒𝑡subscript𝑅𝑒𝑔𝑡\displaystyle\;i\Delta\tilde{\rho}_{eg}(t)-i\left[\Omega_{e}^{\mathrm{(L)}}(t)% M_{e}-\Omega_{g}^{\mathrm{(L)}}(t)M_{g}\right]\tilde{\rho}_{eg}(t)\;-\frac{i}{% 2\hbar}\left[\sum_{M_{g}^{\prime}}V_{eg^{\prime}}\;\tilde{\rho}_{g^{\prime}g}(% t)-\sum_{M_{e}^{\prime}}V_{e^{\prime}g}\;\tilde{\rho}_{ee^{\prime}}(t)\right]+% R_{eg}(t),italic_i roman_Δ over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) - italic_i [ roman_Ω start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ) italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - roman_Ω start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ) italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ] over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) - divide start_ARG italic_i end_ARG start_ARG 2 roman_ℏ end_ARG [ ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g end_POSTSUBSCRIPT ( italic_t ) - ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g end_POSTSUBSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) ] + italic_R start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) , (13d)

where Δ=ωω0Δ𝜔subscript𝜔0\Delta=\omega-\omega_{0}roman_Δ = italic_ω - italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT denotes the light frequency detuning from resonance, ω0subscript𝜔0\omega_{0}italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the atomic transition frequency in the absence of external fields, Ω(L)(t)=gFμBBtotal(t)/superscriptΩL𝑡subscript𝑔𝐹subscript𝜇𝐵subscript𝐵total𝑡Planck-constant-over-2-pi\Omega^{\mathrm{(L)}}(t)=g_{F}\mu_{B}B_{\mathrm{total}}(t)/\hbarroman_Ω start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ) = italic_g start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT roman_total end_POSTSUBSCRIPT ( italic_t ) / roman_ℏ is the Larmor frequency, and Vegsubscript𝑉𝑒𝑔V_{eg}italic_V start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT is the transition matrix element, which is proportional to the Rabi frequency. The contribution of spontaneous decay to the R(t)𝑅𝑡R(t)italic_R ( italic_t ) terms, obtained from the rate for emission summed over polarizations and integrated over angles, is

Rgg(t)=Γ[Fg,Je]{FeFgLJgJeI}2subscript𝑅𝑔superscript𝑔𝑡Γsubscript𝐹𝑔subscript𝐽𝑒superscriptmatrixsubscript𝐹𝑒subscript𝐹𝑔𝐿subscript𝐽𝑔subscript𝐽𝑒𝐼2\displaystyle R_{gg^{\prime}}(t)=\;\Gamma[F_{g},J_{e}]\begin{Bmatrix}F_{e}&F_{% g}&L\\ J_{g}&J_{e}&I\end{Bmatrix}^{2}italic_R start_POSTSUBSCRIPT italic_g italic_g start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = roman_Γ [ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ] { start_ARG start_ROW start_CELL italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_L end_CELL end_ROW start_ROW start_CELL italic_J start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_I end_CELL end_ROW end_ARG } start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
×Me,Me,MFgMgLM|FeMeρ~ee(t)FgMgLM|FeMe,\displaystyle\times\sum\limits_{M_{e},M^{\prime}_{e},M}\langle F_{g}M_{g}LM|F_% {e}M_{e}\rangle\,\tilde{\rho}_{ee^{\prime}}(t)\,\langle F_{g}M^{\prime}_{g}LM|% F_{e}M^{\prime}_{e}\rangle\,,× ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , italic_M start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , italic_M end_POSTSUBSCRIPT ⟨ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_L italic_M | italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩ over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) ⟨ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_L italic_M | italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩ , (14a)
Ree(t)=Γ[Fg,Je]{FeFgLJgJeI}2ρ~ee(t),subscript𝑅𝑒superscript𝑒𝑡Γsubscript𝐹𝑔subscript𝐽𝑒superscriptmatrixsubscript𝐹𝑒subscript𝐹𝑔𝐿subscript𝐽𝑔subscript𝐽𝑒𝐼2subscript~𝜌𝑒superscript𝑒𝑡\displaystyle R_{ee^{\prime}}(t)=\,-\Gamma[F_{g},J_{e}]\begin{Bmatrix}F_{e}&F_% {g}&L\\ J_{g}&J_{e}&I\end{Bmatrix}^{2}\tilde{\rho}_{ee^{\prime}}(t),italic_R start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) = - roman_Γ [ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ] { start_ARG start_ROW start_CELL italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_L end_CELL end_ROW start_ROW start_CELL italic_J start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_I end_CELL end_ROW end_ARG } start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_t ) , (14b)
Rge(t)=12Γ[Fg,Je]{FeFgLJgJeI}2ρ~ge(t),subscript𝑅𝑔𝑒𝑡12Γsubscript𝐹𝑔subscript𝐽𝑒superscriptmatrixsubscript𝐹𝑒subscript𝐹𝑔𝐿subscript𝐽𝑔subscript𝐽𝑒𝐼2subscript~𝜌𝑔𝑒𝑡\displaystyle R_{ge}(t)=\,-\frac{1}{2}\Gamma[F_{g},J_{e}]\begin{Bmatrix}F_{e}&% F_{g}&L\\ J_{g}&J_{e}&I\end{Bmatrix}^{2}\tilde{\rho}_{ge}(t)\,,italic_R start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG roman_Γ [ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ] { start_ARG start_ROW start_CELL italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_L end_CELL end_ROW start_ROW start_CELL italic_J start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_I end_CELL end_ROW end_ARG } start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_g italic_e end_POSTSUBSCRIPT ( italic_t ) , (14c)
Reg(t)=12Γ[Fg,Je]{FeFgLJgJeI}2ρ~eg(t),subscript𝑅𝑒𝑔𝑡12Γsubscript𝐹𝑔subscript𝐽𝑒superscriptmatrixsubscript𝐹𝑒subscript𝐹𝑔𝐿subscript𝐽𝑔subscript𝐽𝑒𝐼2subscript~𝜌𝑒𝑔𝑡\displaystyle R_{eg}(t)=\,-\frac{1}{2}\Gamma[F_{g},J_{e}]\begin{Bmatrix}F_{e}&% F_{g}&L\\ J_{g}&J_{e}&I\end{Bmatrix}^{2}\tilde{\rho}_{eg}(t),italic_R start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG roman_Γ [ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ] { start_ARG start_ROW start_CELL italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_L end_CELL end_ROW start_ROW start_CELL italic_J start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_CELL start_CELL italic_I end_CELL end_ROW end_ARG } start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT ( italic_t ) , (14d)

where [Fg,Je]=(2Fg+1)(2Je+1)subscript𝐹𝑔subscript𝐽𝑒2subscript𝐹𝑔12subscript𝐽𝑒1[F_{g},J_{e}]=(2F_{g}+1)(2J_{e}+1)[ italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ] = ( 2 italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT + 1 ) ( 2 italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT + 1 ) and ΓΓ\Gammaroman_Γ is the decay rate of the upper level |αeJeketsubscript𝛼𝑒subscript𝐽𝑒|\alpha_{e}J_{e}\rangle| italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ⟩ Tremblay and Jacques (1990); Schmidt et al. (2024). Moreover, L=1𝐿1L=1italic_L = 1 corresponds to a dipole transition, L=2𝐿2L=2italic_L = 2 to a quadrupole, L=3𝐿3L=3italic_L = 3 to an octupole, etc.

II.5 Light absorption profile

Solving the Liouville-von Neumann equation (13) numerically allows the determination of the atomic density matrix at any instant of time. The elements of this matrix are directly related to physical observables. In optical magnetometry experiments with vector beams and atoms, the absorption profile of the light is most commonly observed Castellucci et al. (2021); Qiu et al. (2021). Different approaches can be used to analyze such profiles. In Ref. Castellucci et al. (2021), for example, the approach based on Fermi’s golden rule and spatially dependent partially dressed states has been successfully applied to explain the experimental findings. We take a different approach here in which we focus on the diagonal density matrix elements ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) representing the population of photoexcited atomic states Blum (2012). The method relies on a simple assumption that atoms excited to the upper state must decay back to the ground state by the emission of photons in all directions. As a result, regions with many excited atoms appear darker than those with less excitations as the detector measures the intensity of light in the direction of the incoming beam, see Fig. 1. In other words, high values of ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT imply a large imaginary part of the refractive index of the medium. Thus the analysis of the light absorption profile may be reduced to the analysis of the density matrix elements ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT which depend on the position 𝒃𝒃\bm{b}bold_italic_b of the target atom through the transition amplitudes Veg(vec)subscriptsuperscript𝑉vec𝑒𝑔V^{\mathrm{(vec)}}_{eg}italic_V start_POSTSUPERSCRIPT ( roman_vec ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT, as well as on the properties of the magnetic field through the Larmor frequencies Ω(L)superscriptΩL\Omega^{\mathrm{(L)}}roman_Ω start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT and the angle ϕBsubscriptitalic-ϕ𝐵\phi_{B}italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT.

III Results and Discussion

Refer to caption
Figure 4: The 5sS1/225𝑠superscriptsubscript𝑆1225s\;\,{}^{2}S_{1/2}5 italic_s start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT (F=1𝐹1F=1italic_F = 1) -- 5pP3/225𝑝superscriptsubscript𝑃3225p\;\,{}^{2}P_{3/2}5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT (F=0𝐹0F=0italic_F = 0) transition in 87Rb. The lower sublevels are split by the energy Ωg(L)Planck-constant-over-2-pisuperscriptsubscriptΩ𝑔L\hbar\Omega_{g}^{\mathrm{(L)}}roman_ℏ roman_Ω start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT as given by the Larmor frequency of the atom in the magnetic field. The arrows indicate the interaction with light at zero detuning.
Refer to caption
Figure 5: Population (15) of the 5pP3/225𝑝superscriptsubscript𝑃3225p\;\,{}^{2}P_{3/2}5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT (F=0𝐹0F=0italic_F = 0) excited state of 87Rb as a function of time and impact parameter 𝒃𝒃\bm{b}bold_italic_b. The atom interacts with the vector beam (4) in the presence of either the static magnetic field with Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G and Bosc=0subscript𝐵osc0B_{\mathrm{osc}}=0italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 0 (a),(b) or the AC magnetic field with Bconst=Bosc=1subscript𝐵constsubscript𝐵osc1B_{\mathrm{const}}=B_{\mathrm{osc}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G and ωosc=2π×10subscript𝜔osc2𝜋10\omega_{\mathrm{osc}}=2\pi\times 10italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 10 kHz (c),(d). Calculations were performed for different positions 𝒃𝒃\bm{b}bold_italic_b of the atom relative to the beam center (see the top row of the figure). (a),(c) Dependence on b𝑏bitalic_b for ϕb=60°subscriptitalic-ϕ𝑏60°\phi_{b}=60\degreeitalic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 60 °; (b), (d) Dependence on ϕbsubscriptitalic-ϕ𝑏\phi_{b}italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT for b=200𝑏200b=200italic_b = 200 μ𝜇\muitalic_μm. For all curves, Δ=0Δ0\Delta=0roman_Δ = 0, mγ=±1subscript𝑚𝛾plus-or-minus1m_{\gamma}=\pm 1italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT = ± 1, θk=0.11°subscript𝜃𝑘0.11°\theta_{k}=0.11\degreeitalic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = 0.11 °, and A0=3.07×1014subscript𝐴03.07superscript1014A_{0}=3.07\times 10^{-14}italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 3.07 × 10 start_POSTSUPERSCRIPT - 14 end_POSTSUPERSCRIPT.

In the previous section we have outlined the necessary theory for describing the interaction of vector beams with atoms in the presence of an external magnetic field. This formalism can be applied to analyze transitions between two arbitrary hyperfine-structure levels over a wide frequency range. Here we focus on the 5sS1/225𝑠superscriptsubscript𝑆1225s\;\,{}^{2}S_{1/2}5 italic_s start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT (F=1𝐹1F=1italic_F = 1) -- 5pP3/225𝑝superscriptsubscript𝑃3225p\;\,{}^{2}P_{3/2}5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT (F=0𝐹0F=0italic_F = 0) electric dipole (E1) transition in 87Rb at zero detuning (i.e. ω=ω0=2π×384𝜔subscript𝜔02𝜋384\omega=\omega_{0}=2\pi\times 384italic_ω = italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 italic_π × 384 THz, see Fig. 4). This transition has already been utilized in the atomic magnetometer based on vector beams Castellucci et al. (2021). The atom is assumed to be initially unpolarized. The required reduced matrix element 5pP3/22Hγ(E1)5sS1/22delimited-⟨⟩5𝑝superscriptsubscript𝑃322normsubscript𝐻𝛾𝐸15𝑠superscriptsubscript𝑆122\left<5p\;\,{}^{2}P_{3/2}||H_{\gamma}(E1)||5s\;\,{}^{2}S_{1/2}\right>⟨ 5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT | | italic_H start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ( italic_E 1 ) | | 5 italic_s start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ⟩ and the spontaneous decay rate ΓΓ\Gammaroman_Γ have been calculated using the JAC code Fritzsche (2019). Here the theoretically obtained value Γtheo=4.042107subscriptΓtheo4.042superscript107\Gamma_{\mathrm{theo}}=4.042\cdot 10^{7}roman_Γ start_POSTSUBSCRIPT roman_theo end_POSTSUBSCRIPT = 4.042 ⋅ 10 start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT s-1 is relatively close to the measured decay rate Γexp=3.811107subscriptΓexp3.811superscript107\Gamma_{\mathrm{exp}}=3.811\cdot 10^{7}roman_Γ start_POSTSUBSCRIPT roman_exp end_POSTSUBSCRIPT = 3.811 ⋅ 10 start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT s-1 of the P3/22superscriptsubscript𝑃322{}^{2}P_{3/2}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT excited state Volz and Schmoranzer (1996). In what follows we shall only deal with the vector potential (4) with the total angular momentum projections mγ=±1subscript𝑚𝛾plus-or-minus1m_{\gamma}=\pm 1italic_m start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT = ± 1. We have chosen the parameters A0=3.07×1014subscript𝐴03.07superscript1014A_{0}=3.07\times 10^{-14}italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 3.07 × 10 start_POSTSUPERSCRIPT - 14 end_POSTSUPERSCRIPT and θk=0.11°subscript𝜃𝑘0.11°\theta_{k}=0.11\degreeitalic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = 0.11 ° such that the Bessel solution (4) reproduces the experimentally realistic Laguerre-Gaussian mode of waist 200200200200 μ𝜇\muitalic_μm and total power 0.40.40.40.4 μ𝜇\muitalic_μW in the vicinity of the beam center. Such choice of parameters produce Rabi frequencies in the range of MHz in the regions of high beam intensity. We assume that the magnetic field frequency ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT is much smaller than these Rabi frequencies and the decay rate so that the atom-light interaction can be considered adiabatic. Moreover, we also suppose that ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT is much smaller than the Larmor frequency Ω(L)(t)superscriptΩL𝑡\Omega^{\mathrm{(L)}}(t)roman_Ω start_POSTSUPERSCRIPT ( roman_L ) end_POSTSUPERSCRIPT ( italic_t ), which is in the range of MHz as well. For this reason we neglect magnetic field induced transitions, that can occur between Zeeman substates within the same level.

Refer to caption
Figure 6: Excited-state populations for different frequencies of the magnetic field: (a) ωosc=2π×10subscript𝜔osc2𝜋10\omega_{\mathrm{osc}}=2\pi\times 10italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 10 kHz, (b) ωosc=2π×100subscript𝜔osc2𝜋100\omega_{\mathrm{osc}}=2\pi\times 100italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 100 kHz, and (c) ωosc=2π×200subscript𝜔osc2𝜋200\omega_{\mathrm{osc}}=2\pi\times 200italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 200 kHz. It is assumed that b=200𝑏200b=200italic_b = 200 μ𝜇\muitalic_μm, ϕb=60°subscriptitalic-ϕ𝑏60°\phi_{b}=60\degreeitalic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 60 °, Bconst=Bosc=1subscript𝐵constsubscript𝐵osc1B_{\mathrm{const}}=B_{\mathrm{osc}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G (solid lines), Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G and Bosc=0subscript𝐵osc0B_{\mathrm{osc}}=0italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 0 (dashed lines). Vertical bars indicate period of magnetic field oscillations TBsubscript𝑇𝐵T_{B}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT. All other parameters are the same as in Fig. 5.

III.1 Time evolution of the excited-state population

We start our discussion by considering the time evolution of the excited-state population

ρee(t)=5pP3/22(F=0)|ρ^(t)|5pP3/22(F=0)subscript𝜌𝑒𝑒𝑡quantum-operator-product5𝑝superscriptsubscript𝑃322𝐹0^𝜌𝑡5𝑝superscriptsubscript𝑃322𝐹0\displaystyle\rho_{ee}(t)=\left\langle 5p\;\,{}^{2}P_{3/2}\,(F=0)\left|\hat{% \rho}(t)\right|5p\;\,{}^{2}P_{3/2}\,(F=0)\right\rangle\,italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) = ⟨ 5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ( italic_F = 0 ) | over^ start_ARG italic_ρ end_ARG ( italic_t ) | 5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ( italic_F = 0 ) ⟩ (15)

for several selected positions 𝒃=(bcosϕb,bsinϕb,0)𝒃𝑏subscriptitalic-ϕ𝑏𝑏subscriptitalic-ϕ𝑏0\bm{b}=(b\cos\phi_{b},b\sin\phi_{b},0)bold_italic_b = ( italic_b roman_cos italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , italic_b roman_sin italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , 0 ) of the target atom with respect to the zero-intensity center of the vector beam. The results of the calculation are shown in Fig. 5. In the left column, the azimuthal angle is fixed at ϕb=60°subscriptitalic-ϕ𝑏60°\phi_{b}=60\degreeitalic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 60 ° and the radial distance is varied (b=100𝑏100b=100italic_b = 100, 150150150150, 200200200200, 250250250250, and 300300300300 μ𝜇\muitalic_μm). In the right column, the radial distance is fixed rather at b=200𝑏200b=200italic_b = 200 μ𝜇\muitalic_μm and the azimuthal angle is varied (ϕb=0°subscriptitalic-ϕ𝑏0°\phi_{b}=0\degreeitalic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 0 °, 30°30°30\degree30 °, 45°45°45\degree45 °, 60°60°60\degree60 °, and 90°90°90\degree90 °). The subfigures (a) and (b) indicate the population ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) of the exited state for the atoms in the external magnetic field with Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G and Bosc=0subscript𝐵osc0B_{\mathrm{osc}}=0italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 0. For this static magnetic field, ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) reaches a steady state after several tens of microseconds irrespective of the atomic position. The explicit value of the steady-state population ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT is, however, very sensitive to b𝑏bitalic_b and ϕbsubscriptitalic-ϕ𝑏\phi_{b}italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT. Here the dependence of ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT on b𝑏bitalic_b is mainly due to the radial distance dependence of the light intensity. For example, ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT is greater at b=200𝑏200b=200italic_b = 200 μ𝜇\muitalic_μm than at b=100𝑏100b=100italic_b = 100 μ𝜇\muitalic_μm and b=300𝑏300b=300italic_b = 300 μ𝜇\muitalic_μm, since the light intensity is higher at this point. The variation of the excited-state population with the azimuthal angle ϕbsubscriptitalic-ϕ𝑏\phi_{b}italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT can in turn be understood with the help of Fermi’s golden rule for the electric dipole transition rate

W|Veg|2|αeFeMe|𝒆𝒅|αgFgMg|2,proportional-to𝑊superscriptsubscript𝑉𝑒𝑔2superscriptquantum-operator-productsubscript𝛼𝑒subscript𝐹𝑒subscript𝑀𝑒𝒆𝒅subscript𝛼𝑔subscript𝐹𝑔subscript𝑀𝑔2\displaystyle W\propto\left|V_{eg}\right|^{2}\approx\left|\left\langle\alpha_{% e}F_{e}M_{e}\left|\bm{e}\cdot\bm{d}\right|\alpha_{g}F_{g}M_{g}\right\rangle% \right|^{2}\,,italic_W ∝ | italic_V start_POSTSUBSCRIPT italic_e italic_g end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≈ | ⟨ italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT | bold_italic_e ⋅ bold_italic_d | italic_α start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⟩ | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (16)

where 𝒅𝒅\bm{d}bold_italic_d is the dipole moment operator of the atom Bransden and Joachain (2003). As seen from Eq. (16), the absorption depends on the direction 𝒆𝒆\bm{e}bold_italic_e of local polarization of light with respect to the quantization axis given by the magnetic field. Since the polarization of the vector beam (4) varies with the azimuthal angle, the transition rate, and hence ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT, is different for different ϕbsubscriptitalic-ϕ𝑏\phi_{b}italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT. For instance, ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT is much lower at ϕb=90°subscriptitalic-ϕ𝑏90°\phi_{b}=90\degreeitalic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 90 ° than at ϕb=30°subscriptitalic-ϕ𝑏30°\phi_{b}=30\degreeitalic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 30 ° or 60°60°60\degree60 °. A more detailed discussion of this angular dependence with explicit expressions for W𝑊Witalic_W can be found in the work Castellucci et al. (2021).

Refer to caption
Figure 7: Absorption profile of the vector beam at different times t=TB/4𝑡subscript𝑇𝐵4t=T_{B}/4italic_t = italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / 4, TB/2subscript𝑇𝐵2T_{B}/2italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / 2, and TBsubscript𝑇𝐵T_{B}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT for rubidium atoms immersed in (a) the static magnetic field with Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G, Bosc=0subscript𝐵osc0B_{\mathrm{osc}}=0italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 0 and the time-dependent magnetic field with (b) Bconst=Bosc=1subscript𝐵constsubscript𝐵osc1B_{\mathrm{const}}=B_{\mathrm{osc}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G, ωosc=2π×10subscript𝜔osc2𝜋10\omega_{\mathrm{osc}}=2\pi\times 10italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 10 kHz, and (c) Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G, Bosc=10subscript𝐵osc10B_{\mathrm{osc}}=10italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 10 G, ωosc=2π×18subscript𝜔osc2𝜋18\omega_{\mathrm{osc}}=2\pi\times 18italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 18 kHz. The darker shading indicates the maximum absorption and the brighter shading represents the minimum absorption. Parameters are the same as in Fig. 5.

In contrast to the static magnetic field, the density matrix never reaches a steady state when the atoms are immersed in a time-dependent magnetic field. This is illustrated in Fig. 5 (c) and (d) for ωosc=2π×10subscript𝜔osc2𝜋10\omega_{\mathrm{osc}}=2\pi\times 10italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 10 kHz and Bconst=Bosc=1subscript𝐵constsubscript𝐵osc1B_{\mathrm{const}}=B_{\mathrm{osc}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G, i.e. 45°ϕB(t)+45°45°subscriptitalic-ϕ𝐵𝑡45°-45\degree\leq\phi_{B}(t)\leq+45\degree- 45 ° ≤ italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) ≤ + 45 °. Here the excited-state population ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) undergoes oscillations whose character strongly depends on the atomic position. Moreover, the period of these oscillations is shorter than the period of magnetic field oscillations TB=2π/ωosc=100subscript𝑇𝐵2𝜋subscript𝜔osc100T_{B}=2\pi/\omega_{\mathrm{osc}}=100italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT = 2 italic_π / italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 100 μ𝜇\muitalic_μs. As seen from Fig. 6, the shape of the oscillations of ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) changes also with the frequency of the applied magnetic field, and this is most noticeable when ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT increases from 2π×102𝜋102\pi\times 102 italic_π × 10 kHz to 2π×1002𝜋1002\pi\times 1002 italic_π × 100 kHz.

III.2 Time evolution of the light absorption profile

Let us next discuss the excited-state population ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) in the entire x𝑥xitalic_x-y𝑦yitalic_y plane. According to the discussion in Sec. II.5, this distribution can be related to the light absorption profile. Fig. 7 shows such profiles at different times, t=TB/4𝑡subscript𝑇𝐵4t=T_{B}/4italic_t = italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / 4, TB/2subscript𝑇𝐵2T_{B}/2italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / 2, and TBsubscript𝑇𝐵T_{B}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT. One should notice that TBsubscript𝑇𝐵T_{B}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT is depending on ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT, and hence each subfigure features different moments in time. In the figure, the rate of absorption is indicated by shading: lighter shades are minimum values (weak absorption) and darker shades are maximum values (strong absorption). For the case of an external static magnetic field, we observe a stationary petal-like pattern, see Fig. 7 (a). This result agrees well with experiment and previous theoretical predictions Castellucci et al. (2021). If the atoms are additionally exposed to an oscillating magnetic field, the petal-like pattern begins to rotate about the beam axis. As seen from Figs. 7 (b) and (c), the rotation itself depends on the field parameters. This dependence can be used, for example, to determine the strength Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT and frequency ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT of an unknown oscillating component of the magnetic field if the static component Bconstsubscript𝐵constB_{\mathrm{const}}italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT is known.

III.3 Time-averaged light absorption profile

Direct observation of the above-mentioned rotation of the petal-like absorption pattern can be difficult because of the limited time resolution of typical detectors and the requirement that the experiment would need to be executed repeatedly to obtain the time traces. The light absorption profile averaged over a cycle of oscillation

ρ¯ee=1TBTBρee(t)𝑑tsubscript¯𝜌𝑒𝑒1subscript𝑇𝐵subscriptsubscript𝑇𝐵subscript𝜌𝑒𝑒𝑡differential-d𝑡\displaystyle\overline{\rho}_{ee}=\frac{1}{T_{B}}\int_{T_{B}}\rho_{ee}(t)\,dtover¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT end_ARG ∫ start_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) italic_d italic_t (17)

can be considered a more realistic characteristic to measure. It is worth noting that since oscillations of ρee(t)subscript𝜌𝑒𝑒𝑡\rho_{ee}(t)italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT ( italic_t ) repeat themselves at each cycle of magnetic field oscillations, averaging over one period TBsubscript𝑇𝐵T_{B}italic_T start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT and averaging over many periods give the same results. The time-averaged light absorption profile for Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G, Bosc=0.3subscript𝐵osc0.3B_{\mathrm{osc}}=0.3italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 0.3 G, and ωosc=2π×100subscript𝜔osc2𝜋100\omega_{\mathrm{osc}}=2\pi\times 100italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 100 kHz as well as the polar plot of absorption intensity at b=200𝑏200b=200italic_b = 200 μ𝜇\muitalic_μm are displayed in Fig. 8 (a) and (b), respectively. Fig. 8 (c) similarly shows polar plots of ρ¯eesubscript¯𝜌𝑒𝑒\overline{\rho}_{ee}over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT at different strengths Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT of the oscillating component of the magnetic field, assuming that ωosc=2π×100subscript𝜔osc2𝜋100\omega_{\mathrm{osc}}=2\pi\times 100italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 100 kHz and Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G. The petals adjacent to each other in the time-averaged profile begin to merge as Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT varies from 0.30.30.30.3 to 1111 G. With a further increase of the magnetic field strength up to Bosc=5subscript𝐵osc5B_{\mathrm{osc}}=5italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 5 G, the pattern rotates by 45°45°45\degree45 °. This effect can be explained by the time evolution of the magnetic field direction ϕBsubscriptitalic-ϕ𝐵\phi_{B}italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT. For Bosc>2subscript𝐵osc2B_{\mathrm{osc}}>2italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT > 2 G the time evolution of ϕBsubscriptitalic-ϕ𝐵\phi_{B}italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT resembles a rectangle function where most of the time the magnetic field is almost parallel to the y𝑦yitalic_y-axis, see Fig. 2. In the case of such a magnetic field, absorption of the incident radiation (4) is most pronounced for atoms at the angles ϕb=0subscriptitalic-ϕ𝑏superscript0\phi_{b}=0^{\circ}italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 0 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT and 90superscript9090^{\circ}90 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT, as follows from Eq. (16). In contrast, for Bosc2subscript𝐵osc2B_{\mathrm{osc}}\leq 2italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT ≤ 2 G the direction ϕB(t)subscriptitalic-ϕ𝐵𝑡\phi_{B}(t)italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) oscillates sinusoidal around the x𝑥xitalic_x-axis. Here stronger absorption of (4) occurs nearby ϕb=45subscriptitalic-ϕ𝑏superscript45\phi_{b}=45^{\circ}italic_ϕ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 45 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT. This is especially visible in the polar plots in Fig. 8 (c). An even larger enhancement of the magnetic field component Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT, however, will not significantly affect the absorption profile ρ¯eesubscript¯𝜌𝑒𝑒\overline{\rho}_{ee}over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT. Such loss of sensitivity is explained by the fact that the magnetic field direction ϕB(t)subscriptitalic-ϕ𝐵𝑡\phi_{B}(t)italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_t ) ceases to depend on Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT at high Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT, see Fig. 2.

We are also interested in how the magnetic field frequency ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT affects the time-averaged absorption profile ρ¯eesubscript¯𝜌𝑒𝑒\overline{\rho}_{ee}over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT. It can be seen from Fig. 9 that the minima of ρ¯eesubscript¯𝜌𝑒𝑒\overline{\rho}_{ee}over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT become less pronounced as ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT increases from zero in the case of Bconst=Bosc=1subscript𝐵constsubscript𝐵osc1B_{\mathrm{const}}=B_{\mathrm{osc}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G. We also see that the absorption profile ceases to change noticeably at ωosc>2π×100subscript𝜔osc2𝜋100\omega_{\mathrm{osc}}>2\pi\times 100italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT > 2 italic_π × 100 kHz. Thus the determination of the magnetic field based on the analysis of ρ¯eesubscript¯𝜌𝑒𝑒\overline{\rho}_{ee}over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT is difficult at such high magnetic field frequencies.

Refer to caption
Figure 8: (a) Time-averaged absorption profile of the vector beam in the presence of the time-dependent magnetic field with Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G, Bosc=0.3subscript𝐵osc0.3B_{\mathrm{osc}}=0.3italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 0.3 G, and ωosc=2π×100subscript𝜔osc2𝜋100\omega_{\mathrm{osc}}=2\pi\times 100italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 100 kHz. All other parameters are the same as in Fig. 5. (b) Polar plot of the absorption intensity which refers to the radius of the circle in (a) at b=200𝑏200b=200italic_b = 200 μ𝜇\muitalic_μm. (c) Polar plots of absorption intensity as a function of Boscsubscript𝐵oscB_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT. Here Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G and ωosc=2π×100subscript𝜔osc2𝜋100\omega_{\mathrm{osc}}=2\pi\times 100italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 100 kHz.
Refer to caption
Figure 9: (a) Time-averaged absorption profile of the vector beam in the presence of the time-dependent magnetic field with Bconst=Bosc=1subscript𝐵constsubscript𝐵osc1B_{\mathrm{const}}=B_{\mathrm{osc}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G and ωosc=2π×10subscript𝜔osc2𝜋10\omega_{\mathrm{osc}}=2\pi\times 10italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 2 italic_π × 10 kHz. All other parameters are the same as in Fig. 5. (b) Polar plot of the absorption intensity which refers to the radius of the circle in (a) at b=200𝑏200b=200italic_b = 200 μ𝜇\muitalic_μm. (c) Polar plots of absorption intensity as a function of ωoscsubscript𝜔osc\omega_{\mathrm{osc}}italic_ω start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT. Here Bconst=Bosc=1subscript𝐵constsubscript𝐵osc1B_{\mathrm{const}}=B_{\mathrm{osc}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT = 1 G.

In the calculations above we have investigated the sensitivity of absorption profile to the strength and frequency of the oscillating magnetic field for a particular choice of computational parameters. For instance, in Figs. 8 and 9 we have chosen Bconst=1subscript𝐵const1B_{\mathrm{const}}=1italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT = 1 G and 0.30.30.30.3 G Boscabsentsubscript𝐵oscabsent\leq B_{\mathrm{osc}}\leq≤ italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT ≤ 10101010 G. The conclusion based on these results, however, can be generalized on a wider parameter range. In particular, we performed additional calculations which indicated that sensitivity to the AC magnetic field strength and frequency is most pronounced for the case when strengths of reference and test field are comparable to each other, i.e. when 0.1Bconst/Bosc10less-than-or-approximately-equals0.1subscript𝐵constsubscript𝐵oscless-than-or-approximately-equals100.1\lessapprox B_{\mathrm{const}}/B_{\mathrm{osc}}\lessapprox 100.1 ⪅ italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT / italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT ⪅ 10.

In the present study, moreover, we have mainly focused on the scenario where the static (reference) and the oscillating (test) magnetic field are perpendicular to each other. As mentioned above, our theoretical approach is general and can be used for any angle between both fields. For the sake of brevity, we will not discuss here results for different field directions in detail and just mention briefly the most important findings. Namely, we found that time-averaged absorption profile is not sensitive to the strength and frequency of the oscillating magnetic field when 𝑩const𝑩ACconditionalsubscript𝑩constsubscript𝑩AC\bm{B}_{\mathrm{const}}\parallel\bm{B}_{\mathrm{AC}}bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT ∥ bold_italic_B start_POSTSUBSCRIPT roman_AC end_POSTSUBSCRIPT. This sensitivity grows with the angle between both fields and reaches its maximum when 𝑩const𝑩ACperpendicular-tosubscript𝑩constsubscript𝑩AC\bm{B}_{\mathrm{const}}\perp\bm{B}_{\mathrm{AC}}bold_italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT ⟂ bold_italic_B start_POSTSUBSCRIPT roman_AC end_POSTSUBSCRIPT. Such dependence on the angle between both magnetic fields may allow to determine the direction and properties of the AC test field, but this would require application of the DC reference field at several directions and measurements of the absorption profile.

IV Summary and outlook

In this paper we have studied the propagation of a vector light beam through an atomic target, exposed to an external magnetic field having static and oscillating components. Special attention was paid to the absorption profile of transmitted light, as it is often measured in experiments. To compute the absorption profile, we used the density matrix theory. The resulting expressions are general and applicable to any atomic system. As an example, calculations have been performed for the 5sS1/225𝑠superscriptsubscript𝑆1225s\;\,{}^{2}S_{1/2}5 italic_s start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT (F=1𝐹1F=1italic_F = 1) -- 5pP3/225𝑝superscriptsubscript𝑃3225p\;\,{}^{2}P_{3/2}5 italic_p start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_P start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT (F=0𝐹0F=0italic_F = 0) transition in 87Rb. These calculations have shown that the absorption profile of a vector beam varies with time, and its temporal evolution depends on the parameters of the oscillating magnetic field. A measurable signature of the temporal changes to the population of ρeesubscript𝜌𝑒𝑒\rho_{ee}italic_ρ start_POSTSUBSCRIPT italic_e italic_e end_POSTSUBSCRIPT and hence the absorption can be found in the time-averaged absorption profile. We find that the averaged absorption profile is sensitive to both the strength and frequency of the oscillating magnetic field. For the static field of about one Gauss, the highest sensitivity was observed at field strengths in the range from zero to several Gauss and at frequencies in the range from zero to a hundred kHz.

Our study indicates that measurements of the absorption profile of vector light beams can be utilized to diagnose oscillating magnetic fields. The combination with a DC (reference) field enables one to extract information about frequency and strength of an AC (test) field. Based on our calculations, we found that sensitivity to both these parameters is most pronounced when the ratio Bconst/Boscsubscript𝐵constsubscript𝐵oscB_{\mathrm{const}}/B_{\mathrm{osc}}italic_B start_POSTSUBSCRIPT roman_const end_POSTSUBSCRIPT / italic_B start_POSTSUBSCRIPT roman_osc end_POSTSUBSCRIPT is close to unity.

In the present publication several assumptions have been made to simplify our theoretical treatment. In particular, we restricted our work to (i) the 87Rb D2 line induced by (ii) vector mode with particular polarization pattern given by Eq. (4), interacting with (iii) cold atoms whose center-of-mass motion was neglected. While these assumptions are feasible for analysis of current experiments, they have to be questioned for optimizing future measurement setups. In a forthcoming study, we therefore plan to investigate coupling of hot atomic gas with vector light mode exhibiting richer polarization structure and pay special attention to operation of hyperfine transitions.

Acknowledgments

We acknowledge support from the Research School of Advanced Photon Science of the Helmholtz Institute Jena, HPC cluster DRACO of FSU Jena, QuantERA II Programme with funding received via the EU H2020 research and innovation programme under Grant No. 101017733, EPSRC under Grant No. EP/Z000513/1 (V-MAG), and Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy- EXC-2123 QuantumFrontiers-390837967. The authors are grateful for fruitful discussions with K. Essink, N. Huntemann, and E. Peik.

References

  • edited by D. Budker and Kimball (2013) edited by D. Budker and D. F. J. Kimball, Optical Magnetometry (Cambridge University Press, Cambridge, England, 2013).
  • Budker and Romalis (2007) D. Budker and M. Romalis, Nat. Phys. 3, 227 (2007).
  • Johnson et al. (2013) C. N. Johnson, P. D. D. Schwindt,  and M. Weisend, Phys. Med. Biol. 58, 6065 (2013).
  • Vasilakis et al. (2009) G. Vasilakis, J. M. Brown, T. W. Kornack,  and M. V. Romalis, Phys. Rev. Lett. 103, 261801 (2009).
  • Dang et al. (2010) H. B. Dang, A. C. Maloof,  and M. V. Romalis, Appl. Phys. Lett. 97, 151110 (2010).
  • Shah et al. (2007) V. Shah, S. Knappe, P. D. D. Schwindt,  and J. Kitching, Nat. Photonics 1, 649 (2007).
  • Rubinsztein-Dunlop et al. (2017) H. Rubinsztein-Dunlop, A. Forbes, M. V. Berry, M. R. Dennis, D. L. Andrews, M. Mansuripur, C. Denz, C. Alpmann, P. Banzer, T. Bauer, E. Karimi, L. Marrucci, M. Padgett, M. Ritsch-Marte, N. M. Litchinitser, N. P. Bigelow, C. Rosales-Guzmán, A. Belmonte, J. P. Torres, T. W. Neely, M. Baker, R. Gordon, A. B. Stilgoe, J. Romero, A. G. White, R. Fickler, A. E. Willner, G. Xie, B. McMorran,  and A. M. Weiner, J. Opt. 19, 013001 (2017).
  • Gbur (2017) G. J. Gbur, Singular optics (CRC press, 2017).
  • Wang et al. (2020) J. Wang, F. Castellucci,  and S. Franke-Arnold, AVS Quantum Sci. 2, 031702 (2020).
  • Castellucci et al. (2021) F. Castellucci, T. W. Clark, A. Selyem, J. Wang,  and S. Franke-Arnold, Phys. Rev. Lett. 127, 233202 (2021).
  • Qiu et al. (2021) S. Qiu, J. Wang, F. Castellucci, M. Cao, S. Zhang, T. W. Clark, S. Franke-Arnold, H. Gao,  and F. Li, Photon. Res. 9, 2325 (2021).
  • Maldonado et al. (2024) M. G. Maldonado, O. Rollins, A. Toyryla, J. A. McKelvy, A. Matsko, I. Fan, Y. Li, Y.-J. Wang, J. Kitching, I. Novikova,  and E. E. Mikhailov, Opt. Express 32, 25062 (2024).
  • Savukov et al. (2005) I. M. Savukov, S. J. Seltzer, M. V. Romalis,  and K. L. Sauer, Phys. Rev. Lett. 95, 063004 (2005).
  • Ledbetter et al. (2007) M. P. Ledbetter, V. M. Acosta, S. M. Rochester, D. Budker, S. Pustelny,  and V. V. Yashchuk, Phys. Rev. A 75, 023405 (2007).
  • Bransden and Joachain (2003) B. H. Bransden and C. J. Joachain, Physics of Atoms and Molecules (Prentice Hall, Harlow, England, 2003).
  • Knyazev and Serbo (2018) B. A. Knyazev and V. G. Serbo, Phys.-Usp. 61, 449 (2018).
  • Matula et al. (2013) O. Matula, A. G. Hayrapetyan, V. G. Serbo, A. Surzhykov,  and S. Fritzsche, J. Phys. B 46, 205002 (2013).
  • Schulz et al. (2020) S. A.-L. Schulz, A. A. Peshkov, R. A. Müller, R. Lange, N. Huntemann, C. Tamm, E. Peik,  and A. Surzhykov, Phys. Rev. A 102, 012812 (2020).
  • Johnson (2007) W. R. Johnson, Atomic Structure Theory (Springer, New York, 2007).
  • Solyanik-Gorgone et al. (2019) M. Solyanik-Gorgone, A. Afanasev, C. E. Carlson, C. T. Schmiegelow,  and F. Schmidt-Kaler, J. Opt. Soc. Am. B 36, 565 (2019).
  • Peshkov et al. (2023) A. A. Peshkov, E. Jordan, M. Kromrey, K. K. Mehta, T. E. Mehlstäubler,  and A. Surzhykov, Ann. Phys. 535, 2300204 (2023).
  • Rose (1957) M. E. Rose, Elementary Theory of Angular Momentum (John Wiley & Sons, New York, 1957).
  • Blum (2012) K. Blum, Density Matrix Theory and Applications (Springer, Berlin, 2012).
  • Auzinsh et al. (2010) M. Auzinsh, D. Budker,  and S. M. Rochester, Optically Polarized Atoms: Understanding Light-Atom Interactions (Oxford University, Oxford, 2010).
  • von der Wense et al. (2020) L. von der Wense, P. V. Bilous, B. Seiferle, S. Stellmer, J. Weitenberg, P. G. Thirolf, A. Pálffy,  and G. Kazakov, Eur. Phys. J. A 56, 176 (2020).
  • Tremblay and Jacques (1990) P. Tremblay and C. Jacques, Phys. Rev. A 41, 4989 (1990).
  • Schmidt et al. (2024) R. P. Schmidt, S. Ramakrishna, A. A. Peshkov, N. Huntemann, E. Peik, S. Fritzsche,  and A. Surzhykov, Phys. Rev. A 109, 033103 (2024).
  • Fritzsche (2019) S. Fritzsche, Comput. Phys. Commun. 240, 1 (2019).
  • Volz and Schmoranzer (1996) U. Volz and H. Schmoranzer, Physica Scripta T65, 48 (1996).