arXiv:1612.05529 · Computer Physics Communications

ARC: An open-source library for calculating properties of alkali Rydberg atoms

Nikola Šibalić, Jonathan D. Pritchard, Charles S. Adams, Kevin J. Weatherill

Published Aug 5, 2017• 466 citations

Abstract

We present an object-oriented Python library for the computation of properties of highly-excited Rydberg states of alkali atoms. These include single-body effects such as dipole matrix elements, excited-state lifetimes (radiative and black-body limited) and Stark maps of atoms in external electric fields, as well as two-atom interaction potentials accounting for dipole and quadrupole coupling effects valid at both long and short range for arbitrary placement of the atomic dipoles. The package is cross-referenced to precise measurements of atomic energy levels and features extensive documentation to facilitate rapid upgrade or expansion by users. This library has direct application in the field of quantum information and quantum optics which exploit the strong Rydberg dipolar interactions for two-qubit gates, robust atom-light interfaces and simulating quantum many-body physics, as well as the field of metrology using Rydberg atoms as precise microwave electrometers. Program Title: ARC: Alkali Rydberg Calculator Program Files doi: http://dx.doi.org/10.17632/hm5n8w628c.1 Licensing provisions: BSD-3-Clause Programming language: Python 2.7 or 3.5, with C extension External Routines: NumPy [1], SciPy [1], Matplotlib [2] Nature of problem: Calculating atomic properties of alkali atoms including lifetimes, energies, Stark shifts and dipole–dipole interaction strengths using matrix elements evaluated from radial wavefunctions. Solution method: Numerical integration of radial Schrödinger equation to obtain atomic wavefunctions, which are then used to evaluate dipole matrix elements. Properties are calculated using second order perturbation theory or exact diagonalisation of the interaction Hamiltonian, yielding results valid even at large external fields or small interatomic separation. Restrictions: External electric field fixed to be parallel to quantisation axis. Supplementary material: Detailed documentation (.html), and Jupyter notebook with examples and benchmarking runs (.html and .ipynb). [1] T.E. Oliphant, Comput. Sci. Eng. 9, 10 (2007). http://www.scipy.org/. [2] J.D. Hunter, Comput. Sci. Eng. 9, 90 (2007). http://matplotlib.org/.

DOI: 10.1016/j.cpc.2017.06.015

Related papers

ARC 3.0: An expanded Python toolbox for atomic physics calculations

Elizabeth Robertson, Nikola Šibalić, R. M. Potvliege +1 more

ARC 3.0 is a modular, object-oriented Python library combining data and algorithms to enable the calculation of a range of properties of alkali and divalent atoms. Building on the initial version of the ARC library (Šibalić et al., 2017), which focused on Rydberg states of alkali atoms, this major upgrade introduces support for divalent atoms. It also adds new methods for working with atom–surface interactions, for modelling ultracold atoms in optical lattices and for calculating valence electron wave functions and dynamic polarisabilities. Such calculations have applications in a variety of fields, e.g., in the quantum simulation of many-body physics, in atom-based sensing of DC and AC fields (including in microwave and THz metrology) and in the development of quantum gate protocols. ARC 3.0 comes with an extensive documentation including numerous examples. Its modular structure facilitates its application to a wide range of problems in atom-based quantum technologies. Program Title: ARC 3.0 CPC Library link to program files: https://doi.org/10.17632/c4z4n2cdf7.1 Licencing provisions: BSD-3-Clause Programming language: Python External Routines: NumPy [1], SciPy [1], Matplotlib [2], SymPy [3], LmFit [4] Nature of problem: The calculation of atomic properties of alkali and divalent atoms including energies, Stark shifts and dipole–dipole interaction strengths using matrix elements evaluated through a variety of means. Solution method: Dipole matrix elements are calculated using an analytical semi-classical approximation or wave functions obtained by numerical integration of the radial Schrödinger equation for a one-electron model potential. Interaction energies and shifts due to external fields are calculated using second order degenerate perturbation theory or exact diagonalisation of the interaction Hamiltonian, yielding results valid even at large external fields or small interatomic separation. Additional comments including restrictions and unusual features: External electric and magnetic field must be parallel to the quantisation axis. The accuracy of short range (≲1μm) atom - atom interaction potentials is limited by the truncation of the basis. Only weak magnetic fields are supported as only linear Zeeman shifts are taken into account. Calculations for divalent atoms use a single-electron approximation and calculation of their wave functions is not supported. References [1] T.E. Oliphant, Comput. Sci. Eng. 9 (2007) 10. http://www.scipy.org/. [2] J.D. Hunter, Comput. Sci. Eng. 9 (2007) 90. http://matplotlib.org/. [3] A. Meurer et al., PeerJ Comput. Sci. 3 (2017) e103. https://doi.org/10.7717/peerj-cs.103 [4] M. Newville et al., lmfit/lmfit-py 1.0.0 (Version 1.0.0), Zenodo (2019). https://doi.org/10.5281/zenodo.3588521

arXiv:2007.12016Computer Physics CommunicationsJan 1, 202196 citations

Calculation of Rydberg interaction potentials

Sebastian Weber, Christoph Tresp, Henri Menke +4 more

Abstract The strong interaction between individual Rydberg atoms provides a powerful tool exploited in an ever-growing range of applications in quantum information science, quantum simulation and ultracold chemistry. One hallmark of the Rydberg interaction is that both its strength and angular dependence can be fine-tuned with great flexibility by choosing appropriate Rydberg states and applying external electric and magnetic fields. More and more experiments are probing this interaction at short atomic distances or with such high precision that perturbative calculations as well as restrictions to the leading dipole–dipole interaction term are no longer sufficient. In this tutorial, we review all relevant aspects of the full calculation of Rydberg interaction potentials. We discuss the derivation of the interaction Hamiltonian from the electrostatic multipole expansion, numerical and analytical methods for calculating the required electric multipole moments and the inclusion of electromagnetic fields with arbitrary direction. We focus specifically on symmetry arguments and selection rules, which greatly reduce the size of the Hamiltonian matrix, enabling the direct diagonalization of the Hamiltonian up to higher multipole orders on a desktop computer. Finally, we present example calculations showing the relevance of the full interaction calculation to current experiments. Our software for calculating Rydberg potentials including all features discussed in this tutorial is available as open source.

arXiv:1612.08053Journal of Physics B Atomic Molecular and Optical PhysicsMay 19, 2017174 citations

A Python program to draw energy level diagram and to calculate quantum defects, the high lying Rydberg level

Hasan Nafees, Saba Javaid, Misha Shafi +4 more

Rydberg atoms have strong and adjustable atomic interactions that can be tweaked by choosing different states with various principal quantum numbers or orbital angular momentum. They are highly used in quantum information applications. Furthermore, these atoms have long lifetimes and many possible energy levels, and their separations enable coupling to electromagnetic fields with frequencies ranging over six orders of magnitude. We studied some properties of these atoms and have developed a python program; the program gives two options. Firstly, it calculates energies of elements in s, p, d, and f orbitals and draws an energy level diagram of the element listed in the program. Secondly, it gives an option for calculating energy and quantum defects of atoms and ions, either listed or not listed in the program and generates their graphs accordingly. There are 21 elements in the database whose energy level diagrams could be drawn. The database can be extended by adding more elements.

Revista Mexicana de Física EJan 13, 20231 citation

Creation and manipulation of quantum states of light and cold atoms using Rydberg states

T Laycock

Rydberg atoms are often proposed as the basis of quantum computing and quantum information protocols. One of the central reasons for this is that they provide a strong and long-ranged interaction that can be coherently switched on and off. This thesis details two techniques which use the exaggerated properties of Rydberg atoms to manipulate both the quantum state of the atom itself and that of the external light field.\n\nThe first proposal initially focuses on the creation of many-body quantum states from two-level atoms trapped in a two-dimensional lattice. This approach uses the van der Waals interaction present between alkali metal atoms in highly excited Rydberg s-states. The approximate solution of the corresponding Hamiltonian is detailed in the regime where the laser driving is the largest energy scale of the system. The states which are most likely achieved using an oscillating laser detuning are then determined. These states are then taken as the basis for the creation of deterministic single-photons, whose properties are shown to rely on the interplay between interatomic spacing and the geometry of the lattice.\n\nThe second technique described uses the coupling between a Rydberg atom and a moving electron to manipulate the atomic state. In this system, the atom is initially excited to a Rydberg s-state and trapped at a finite distance from an electron waveguide. Two analytical methods are used to show that the final state of the atom depends strongly on the direction and modulus of the electron momentum. A complementary numerical simulation shows that the atoms may be left in a polarised state, suggesting the possibility of using this setup to ‘switch on’ permanent electric dipoles in the atoms. This investigation leads naturally to a system where multiple interacting atoms are trapped close to the waveguide, allowing various many-body states to be accessed.

Dec 10, 20130 citations

Controlling Rydberg atoms in external electric fields

Jens Grimmel

In this work, numerical calculations were implemented based on the diagonalisation of a matrix representation of the Stark Hamiltonian on a subset of the basis given by the zero-field Hamiltonian. Dipole matrix elements between the involved states were calculated, which allow for the prediction of transition strengths to Stark shifted Rydberg states. The applicability of the numerical methods were confirmed experimentally by precision spectroscopy in a rubidium vapour cell with electrodes using electromagnetically induced transparency (EIT). The numerical methods were then extended to predict ionisation rates of Rydberg states in the highly Stark shifted regime beyond the classical ionisation threshold. For this purpose, a complex absorbing potential (CAP) was specifically adjusted to the Stark Hamiltonian. An experiment was conducted, in which cold atoms from a magneto-optical trap (MOT) were excited to Stark shifted Rydberg states and the arising ions were detected. The obtained ionisation spectra agree well with the results from the numerical calculations. Another extension of the numerical methods was implemented for crossed electric and magnetic fields. This was used for approximate calculations of the motional Stark effect (MSE), which acts on atoms moving through static magnetic fields and thereby experiencing a Lorentz electric field. In a precision spectroscopy, once more using EIT in a vapour cell, the MSE was observed in low magnetic fields for the first time in rubidium. The spectral shift, which is caused by the MSE, matches the prediction from the numerical calculations. In a final experiment, the practical use of the knowledge of several properties of Stark shifted Rydberg states was illustrated. The numerical methods including the CAP were used to specifically select states for their optical accessibility, permanent electric dipole moments and ionisation rates. In the experiment, atoms from a MOT were excited to these states to observe different characteristics of the dipole blockade effect. The detection scheme for this experiment was realised by controlled ionisation using only small changes in the external electric field.

Jan 1, 20180 citations

ElecSus: Extension to arbitrary geometry magneto-optics

James Keaveney, Charles S. Adams, Ifan G. Hughes

We present a major update to ElecSus, a computer program and underlying model to calculate the electric susceptibility of an alkali-metal atomic vapour. Knowledge of the electric susceptibility of a medium is essential to predict its absorptive and dispersive properties. In this version we implement several changes which significantly extend the range of applications of ElecSus, the most important of which is support for non-axial magnetic fields (i.e. fields which are not aligned with the light propagation axis). Supporting this change requires a much more general approach to light propagation in the system, which we have now implemented. We exemplify many of these new applications by comparing ElecSus to experimental data. In addition, we have developed a graphical user interface front-end which makes the program much more accessible, and have improved on several other minor areas of the program structure. Program Title: ElecSus Program Files doi: http://dx.doi.org/10.17632/h7cj8bz4bd.1 Licensing provisions: Apache License, Version 2.0 Programming language: Python External routines/libraries: SciPy library [1] 0.15.0 or later, NumPy [1], matplotlib [2], sympy [3], lmfit 0.9.5 or later [4], wxpython (required for GUI only) Nature of problem: Calculating the weak-probe electric susceptibility of an alkali-metal vapour. The electric susceptibility can be used to calculate spectra such as transmission and Stokes parameters. Measurements of experimental parameters can be made by fitting the theory to data. Solution method: The transition frequencies and wavelengths are calculated using a matrix representation of the Hamiltonian in the completely uncoupled basis. A suite of fitting methods are provided in order to allow user supplied experimental data to be fit to the theory, thereby allowing experimental parameters to be extracted. Restrictions: Results are only valid in the weak-probe regime. [1] T. E. Oliphant, Comput. Sci. Eng. 9, 10 (2007). http://www.scipy.org/ [2] J. D. Hunter, Comput. Sci. Eng. 9, 10 (2007). http://matplotlib.org/ [3] A. Meurer et. al, PeerJ Comp. Sci. 3, e103 (2017) http://www.sympy.org/ [4] M. Newville et al., LMFIT: Non-Linear Least-Square Minimization and Curve-Fitting for Python, Zenodo (2014). DOI:10.5281/zenodo.11813 https://lmfit.github.io/lmfit-py/

arXiv:1708.05305Computer Physics CommunicationsDec 11, 201758 citations