@bookBudker:2013,
title=Optical Magnetometry,
author=edited by D. Budker and D. F. Jackson Kimball,
year=2013,
publisher=Cambridge University Press, Cambridge, England
@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 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.,
@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
,
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 () () 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.
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 D 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 -axis and having the frequency . We assume that the atoms are exposed to a combination of (reference) static and (test) oscillating magnetic fields. The static is applied along , while oscillates along the -axis, as shown in Fig. 1. The resulting magnetic field has the strength
(1)
oscillating between and , while its direction is characterized by the angle
(2)
varying in the range . A typical variation of and with is shown in Fig. 2. We take fairly close values of and so that the direction of the resulting magnetic field oscillates with a sufficiently large amplitude.
By design, the field is always perpendicular to the light propagation direction and oscillates around . 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 , 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 , 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 . Both choices of the quantization axis lead to the same result.
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:
(3)
where is a constant amplitude, is its frequency, , , and are cylindrical coordinates, and are longitudinal and transverse components of the linear momentum, respectively, and 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
(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 and the projection of the total angular momentum upon the propagation direction. Moreover, its longitudinal momentum and absolute value of the transverse momentum are also fixed, see Ref. Matula et al. (2013) for further details. The vector potential for the Bessel beam can be written as
(5)
where and is a weight function given by:
(6)
It follows from these expressions that the Bessel beam can be seen as a superposition of plane waves whose wave vectors are uniformly distributed upon the surface of a cone with a polar opening angle .
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 . For arbitrary values of the opening angle , 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, Matula et al. (2013). In this regime is small and the vector potential (4) can be considerably simplified to:
(7)
Let us remark here that the vector beam in the paraxial approximation (7) has azimuthally varying linear polarization that lies in the transverse - 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 , 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 and excited atomic states whose properties can be traced back to the first-order matrix element
(8)
with , where and are the nuclear and electron angular momenta, respectively, is the projection of on the quantization axis, and denotes all additional quantum numbers required to specify the state uniquely. Moreover, runs over all electrons in a target atom and denotes the vector of Dirac matrices for the th particle Johnson (2007). In Eq. (8), we have introduced the impact parameter to specify the position of the atom within the wave front Knyazev and Serbo (2018). In particular, 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 can be expressed in terms of its Bessel counterparts as
(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
(10)
where we have used the notation to denote the reduced matrix element for magnetic () and electric () transitions. Furthermore, the Euler angles as the arguments of the Wigner -functions and in Eq. (10) characterize the rotation from the atomic frame with the quantization axis along the magnetic field to the photon frame with the quantization axis along the wave vector Rose (1957). The time-dependent angle 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 satisfying the Liouville-von Neumann equation:
(11)
Here is the total Hamiltonian of the atom in the presence of external fields, and is introduced to take into account phenomenologically spontaneous decay Auzinsh et al. (2010). We consider transitions between the Zeeman sublevels of the ground to those of the excited state , described by a density matrix of size . In this basis, we can write the elements of the density matrix as:
(12a)
(12b)
(12c)
(12d)
In Eqs. (12), the diagonal elements and are the probabilities of finding an atom in the Zeeman substates and , 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 , , , and . To solve these equations, we introduce , , , and employ the rotating-wave approximation, which is valid when 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 , so that we can rewrite the Liouville-von Neumann equation as
(13a)
(13b)
(13c)
(13d)
where denotes the light frequency detuning from resonance, is the atomic transition frequency in the absence of external fields, is the Larmor frequency, and is the transition matrix element, which is proportional to the Rabi frequency. The contribution of spontaneous decay to the terms, obtained from the rate for emission summed over polarizations and integrated over angles, is
(14a)
(14b)
(14c)
(14d)
where and is the decay rate of the upper level Tremblay and Jacques (1990); Schmidt et al. (2024). Moreover, corresponds to a dipole transition, to a quadrupole, 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 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 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 which depend on the position of the target atom through the transition amplitudes , as well as on the properties of the magnetic field through the Larmor frequencies and the angle .
III Results and Discussion
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 () () electric dipole (E1) transition in 87Rb at zero detuning (i.e. 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 and the spontaneous decay rate have been calculated using the JAC code Fritzsche (2019). Here the theoretically obtained value s-1 is relatively close to the measured decay rate s-1 of the excited state Volz and Schmoranzer (1996). In what follows we shall only deal with the vector potential (4) with the total angular momentum projections . We have chosen the parameters and such that the Bessel solution (4) reproduces the experimentally realistic Laguerre-Gaussian mode of waist m and total power 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 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 is much smaller than the Larmor frequency , 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.
III.1 Time evolution of the excited-state population
We start our discussion by considering the time evolution of the excited-state population
(15)
for several selected positions 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 and the radial distance is varied (, , , , and m). In the right column, the radial distance is fixed rather at m and the azimuthal angle is varied (, , , , and ). The subfigures (a) and (b) indicate the population of the exited state for the atoms in the external magnetic field with G and . For this static magnetic field, reaches a steady state after several tens of microseconds irrespective of the atomic position. The explicit value of the steady-state population is, however, very sensitive to and . Here the dependence of on is mainly due to the radial distance dependence of the light intensity. For example, is greater at m than at m and m, since the light intensity is higher at this point. The variation of the excited-state population with the azimuthal angle can in turn be understood with the help of Fermi’s golden rule for the electric dipole transition rate
(16)
where is the dipole moment operator of the atom Bransden and Joachain (2003). As seen from Eq. (16), the absorption depends on the direction 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 , is different for different . For instance, is much lower at than at or . A more detailed discussion of this angular dependence with explicit expressions for can be found in the work Castellucci et al. (2021).
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 kHz and G, i.e. . Here the excited-state population 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 s. As seen from Fig. 6, the shape of the oscillations of changes also with the frequency of the applied magnetic field, and this is most noticeable when increases from kHz to kHz.
III.2 Time evolution of the light absorption profile
Let us next discuss the excited-state population in the entire - 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, , , and . One should notice that is depending on , 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 and frequency of an unknown oscillating component of the magnetic field if the static component 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
(17)
can be considered a more realistic characteristic to measure. It is worth noting that since oscillations of repeat themselves at each cycle of magnetic field oscillations, averaging over one period and averaging over many periods give the same results. The time-averaged light absorption profile for G, G, and kHz as well as the polar plot of absorption intensity at m are displayed in Fig. 8 (a) and (b), respectively. Fig. 8 (c) similarly shows polar plots of at different strengths of the oscillating component of the magnetic field, assuming that kHz and G. The petals adjacent to each other in the time-averaged profile begin to merge as varies from to G. With a further increase of the magnetic field strength up to G, the pattern rotates by . This effect can be explained by the time evolution of the magnetic field direction . For G the time evolution of resembles a rectangle function where most of the time the magnetic field is almost parallel to the -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 and , as follows from Eq. (16). In contrast, for G the direction oscillates sinusoidal around the -axis. Here stronger absorption of (4) occurs nearby . This is especially visible in the polar plots in Fig. 8 (c). An even larger enhancement of the magnetic field component , however, will not significantly affect the absorption profile . Such loss of sensitivity is explained by the fact that the magnetic field direction ceases to depend on at high , see Fig. 2.
We are also interested in how the magnetic field frequency affects the time-averaged absorption profile . It can be seen from Fig. 9 that the minima of become less pronounced as increases from zero in the case of G. We also see that the absorption profile ceases to change noticeably at kHz. Thus the determination of the magnetic field based on the analysis of is difficult at such high magnetic field frequencies.
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 G and G 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 .
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 . This sensitivity grows with the angle between both fields and reaches its maximum when . 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 () () 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 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 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).
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).
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).
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).
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).