Search

The collective interaction of electronic excitations with the electromagnetic field in atomic arrays can lead to significantly reduced decay rates, resulting in the formation of subradiant states with promising applications in quantum information and quantum memories. In many previous studies, atoms are assumed to be perfectly positioned and fixed. In this work, we investigate the effects of finite trap strength on the long-lived subradiant states for two, three, and many atoms in a 1D waveguide or free space. We demonstrate that the decay rate is lower bounded by the spatial spread of the motional wavefunction, which is non-zero even for the vibrational ground state. Additionally, in relatively tight traps, the decay of a shared electronic excitation imparts recoil energy to the atoms, which is proportional to the lifetime of the excited state. In a 1D waveguide, this energy can reach a vibrational quantum, even in the Lamb-Dicke regime. Furthermore, for atoms in weaker traps, the induced dipole forces become significant, causing the atoms to undergo motion prior to decay. This motion results in a time-dependent decay rate, an additional energy transfer to the atoms following decay, and mixing between the different electronic eigenstates. Consequently, these effects lead to decoherence and a reduction in the fidelity of a stored electronic eigenstate. Infidelity can be mitigated by either decreasing the induced forces or increasing the trap strength. For quantum information storage, these considerations suggest optimal array configurations, both in terms of geometric arrangements and polarization directions. Our findings provide insights for the behavior of quantum memories and atom array experiments.

We theoretically and computationally investigate the role that the spatial spread of atoms plays in the transmission and reflection of weak light from atom arrays. In particular, we investigate whether coherent wave functions for the atoms' positions leads to different results from a thermal distribution with the same spatial spread. We find that the coherence is not relevant when the light is weak and the electronic states evolve on time scales shorter than the oscillation period of the atoms in their traps. Full numerical calculations and derivations using the sudden approximation show that reflection and transmission agree with the simple averaging over atom positions for these conditions. For parameters outside these restrictions, the simple spatial averaging may lead to inaccurate results.

We theoretically and computationally investigate the cooling of antihydrogen, $\bar{H}$, using optical molasses cooling. This updates the results in Ref. [1] to the current capabilities of the ALPHA experiment. Through Monte Carlo simulation, we show that $\bar{H}$s do not give the standard cooling even in an ideal optical molasses because of their small mass and large transition frequency. For optical molasses cooling in the ALPHA trap, the photons are constrained to travel in one direction only. It is only through the phase space mixing in the trap that cooling in all directions can be achieved. We explore the nontrivial role that laser intensity plays in the cooling. We also investigate the possibility for simultaneously cooling atoms in either of the trapped ground states.

The driven Dicke model, with interesting quantum phases induced by parameterized driving, has been intensively studied in cavities, where permutation symmetry applies due to the atoms' equal coupling to the field and identical interaction. As a result, the system, with proper initialization, can remain in a highly symmetric subset of the state space, where the photon emission of each atom constructively interferes with each other, leading to superradiance at steady state. However, because of the degeneracy of steady states for the driven Dicke model, the steady state can be qualitatively altered by an infinitesimal perturbation. In this work, we simulate superconducting qubits coupled to a 1D waveguide as the extended system and theoretically investigate four kinds of perturbations: local dephasing, individual driving phases, the separation between adjacent qubits, and individual detunings. Using an angular momentum basis, we predict the dimension of the degenerate subspace and study the transition within the subspace due to the perturbation.

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of point-like particles, typically the dominant process for reaching thermal equilibrium in higher dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential $V_{ij} \propto 1/|x_i-x_j|^n$ or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature $T$ toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low density limit. For the inverse power-law potential, the rate of thermalization is proportional to $\rho^2 T^{\frac{1}{2}-\frac{1}{n}}$ where $\rho$ is the number density. The corresponding proportionality constant decreases with increasing $n$., Comment: 13 pages, 12 figures

Simulations and theory are presented for the power spectral density functions (PSDs) of particles in time dependent and anharmonic potentials including the effects of a thermal environment leading to damping and fluctuating forces. We investigate three one dimensional perturbations to the harmonic oscillator of which two are time dependent changes in the natural frequency of the oscillator, while the other is a time independent extension of the quadratic potential to include a quartic term. We investigate the effect of these perturbations on two PSDs of the motion that are used in experiments on trapped nano-oscillators. We also derive and numerically test the PSDs for the motion of a spherical nanoparticle in a Paul trap. We found that the simple harmonic Langevin oscillator's PSDs are good approximations for the $x$-and $y$-coordinates' PSDs for small values of the parameter $q$ of the Mathieu equation, but the difference can be more than a factor of two as '$q$' increases. We also numerically showed that the presence of a permanent electric dipole on the nanosphere can significantly affect the PSDs in the $x$-and $y$-coordinates., Comment: 15 pages, 15 figures

We theoretically and computationally investigate the behavior of infinite atom arrays when illuminated by nearly resonant light. We use higher order mean field equations to investigate the coherent reflection and transmission and incoherent scattering of photons from a single array and from a pair of arrays as a function of detuning for different values of the Rabi frequency. For the single array case, we show how increasing the light intensity changes the probabilities for these different processes. For example, the incoherent scattering probability initially increases with light intensity before decreasing at higher values. For a pair of parallel arrays at near resonant separation, the effects from increasing light intensity can become apparent with incredibly low intensity light. In addition, we derive the higher order mean field equations for these infinite arrays giving a representation that can be evaluated with a finite number of equations.

Previous work on imaging wave packet dynamics with x-ray scattering revealed that the scattering patterns deviate substantially from the notion of instantaneous momentum density of the wave packet. Here we show that scattering patterns can provide clear insights into the electron wave packet dynamics if the final state of the scattered electron and the scattered photon momentum are determined simultaneously. The scattering probability is shown to be proportional to the modulus square of the Fourier transform of the instantaneous spatial wave function weighted by the final state of the electron. Several cases for the choice of final state of the electron are explored. First, the case where the final state can be measured up to a given principal quantum number $n$ and orbital angular momentum $l$ are presented. Next, the case where the final states can only be determined up to a given energy is discussed. Finally, the case of an initial wave packet consisting of a large amount of a known stationary state and a small amount of an unknown stationary state is examined. The scattering profile is used to determine the properties of the unknown state in the wave packet., Comment: 12 pages, 5 figures

The recoil of atoms in dense ensembles during light matter interactions is studied using quantized vibrational states for the atomic motion. The recoil resulting from the forces due to the near-field collective dipole interactions and far-field laser and decay interactions are explored. The contributions to the recoil and the dependence on the trap frequency of the different terms of the Hamiltonian and Lindbladian are studied. These calculations are compared with previous results using the impulse model in the slow oscillation approximation. Calculations in highly subradiant systems show enhanced recoil indicating that recoil effects cannot be ignored in such cases., Comment: 9 pages, 5 figures

Photoionization of atoms immersed in an environment such as an ultracold gas is investigated. We show that the interference of two ionization pathways, one passing directly to the continuum and one accounting for scattering processes between the photoelectron and a neighboring atom, produces a fine structure in the photoionization cross-section over an energy range less than 1 eV above threshold. This fine structure includes all the details of the corresponding three-body system, e.g. the interatomic distance or the scattering information of the electron-atom subsystem; therefore, photoelectrons produced in a multi-particle environment can be utilized as structural probes. As an illustration, for experimentally relevant parameters, we propose a scheme based on the photoionization of a Rydberg molecule where the low-energy electron-atom phase shifts are extracted from the fine structure spectra using neural networks., Comment: 9 pages (main text and supplemental material), 4 figures

We explore conditions for Dicke superradiance in a cloud of atoms by examining the Taylor series expansion of the photon emission rate at $t= 0$. By defining superradiance as an increasing photon emission rate for $t\sim 0$, we have calculated the conditions for superradiance for a variety of cases. We investigate superradiance as defined for photon emission into all angles as well as directional superradiance where the photon emission is only detected in a particular direction. Although all of the examples are for two level atoms that are fully inverted at $t=0$, we also give equations for partially inverted two level atoms and for fully inverted multilevel atoms. We give an algorithm for efficiently evaluating these equations for atom arrays and determine superradiance conditions for large atom number.

The physics that determines the line shape of the 1S--2S transition in magnetically trapped $\bar{\rm H}$ is explored. Besides obtaining an understanding of the line shape, one goal is to replace the dependence on large scale simulations of $\bar{\rm H}$ with a simpler integration over well defined functions. For limiting cases, analytic formulas are obtained. Example calculations are performed to illustrate the limits of simplifying assumptions. We also describe a $\chi^2$ method for choosing experimental parameters that can lead to the most accurate determination of the transition frequency.

We report the experimental observation of superradiant emission emanating from an elongated dense ensemble of laser cooled two-level atoms, with a radial extent smaller than the transition wavelength. In the presence of a strong driving laser, we observe that the system is superradiant along its symmmetry axis. This occurs even though the driving laser is orthogonal to the superradiance direction. This superradiance modifies the spontaneous emission, and, resultantly, the Rabi oscillations. We also investigate Dicke superradiance in the emission of an almost fully-inverted system as a function of atom numnber. The experimental results are in qualitative agreement with ab-initio, beyond-mean-field calculations.

The ALPHA collaboration has successfully demonstrated the production and the confinement of cold antihydrogen, $\overline{\mathrm{H}}$. An analysis of trapping data allowed a stringent limit to be placed on the electric charge of the simplest antiatom. Charge neutrality of matter is known to a very high precision, hence a neutrality limit of $\overline{\mathrm{H}}$ provides a test of CPT invariance. The experimental technique is based on the measurement of the deflection of putatively charged $\overline{\mathrm{H}}$ in an electric field. The tendency for trapped $\overline{\mathrm{H}}$ atoms to be displaced by electrostatic fields is measured and compared to the results of a detailed simulation of $\overline{\mathrm{H}}$ dynamics in the trap. An extensive survey of the systematic errors is performed, with particular attention to those due to the silicon vertex detector, which is the device used to determine the $\overline{\mathrm{H}}$ annihilation position. The limit obtained on the charge of the $\overline{\mathrm{H}}$ atom is \mbox{$ Q = (-1.3\pm1.8\pm0.4)\times10^{-8}$}, representing the first precision measurement with $\overline{\mathrm{H}}$., Comment: 5 pages, 3 figures

Results from higher order mean field calculations of light interacting with atom arrays are presented for calculations of one- and two-time expectation values. The atoms are approximated as two-levels and are fixed in space. Calculations were performed for mean field approximations that include the expectation value of one operator (mean field), two operators (mean field-2), and three operators (mean field-3). For the one-time expectation values, we examined three different situations to understand the convergence with increasing order of mean field and some limitations of higher order mean field approximations. As a representation of a two-time expectation value, we calculated the $g^{(2)}(\tau )$ for a line of atoms illuminated by a perpendicular plane wave at several emission angles and two different intensities. For many cases, the mean field-2 will be sufficiently accurate to quantitatively predict the response of the atoms as measured by one-time expectation values. However, the mean field-3 approximation will often be needed for two-time expectation values.

We present a study of the Rydberg spectrum in \ts{166}Er for series connected to the $4f^{12} (^3H_6) 6s$, $J_c=13/2 $ and $J_c=11/2 $ ionic core states using an all-optical detection based on electromagnetically induced transparency in an effusive atomic beam. Identifying approximately 550 individual states, we find good agreement with a multi-channel quantum defect theory (MQDT) which allows assignment of most states to $ns$ or $nd$ Rydberg series. We provide an improved accuracy for the lowest two ionization thresholds to $E_{\textrm{IP}, J_c = 13/2 } = 49260.750(1)\,$cm$^{-1}$ and $E_{\textrm{IP}, J_c = 11/2 } = 49701.184(1)\,$cm$^{-1}$ as well as the corresponding quantum defects for all observed series. We identify Rydberg states in five different isotopes, and states between the two lowest ionization thresholds. Our results open the way for future applications of Rydberg states for quantum simulation using erbium and exploiting its special open-shell structure., Comment: 17 pages, 9 figures

The recoil of atoms in arrays due to the emission or absorption of photons is studied for sub-wavelength interatomic spacing. The atoms in the array interact with each other through collective dipole-dipole interactions and with the incident laser field in the low intensity limit. Shining uniform light on the array gives rise to patterns of excitation and recoil in the array. These arise due to the interference of different eigenmodes of excitation. The relation between the recoil and the decay dynamics is studied when the array is in its excitation eigenstates. The recoil experienced by a subradiant collective decay is substantially larger than from independent atom decay. A method to calculate the rate of recoil when steady state has been achieved with a constant influx of photons is also described., Comment: 13 pages, 6 figures

Limits to Rydberg gate fidelity that arise from the entanglement of internal states of neutral atoms with the motional degrees of freedom due to the momentum kick from photon absorption and re-emission is quantified. This occurs when the atom is in a superposition of internal states but only one of these states is manipulated by visible or UV photons. The Schr\"odinger equation that describes this situation is presented and two cases are explored. In the first case, the entanglement arises because the spatial wave function shifts due to the separation in time between excitation and stimulated emission. For neutral atoms in a harmonic trap, the decoherence can be expressed within a sudden approximation when the duration of the laser pulses are shorter than the harmonic oscillator period. In this limit, the decoherence is given by simple analytic formulas that account for the momentum of the photon, the temperature of the atoms, the harmonic oscillator frequency, and atomic mass. In the second case, there is a reduction in gate fidelity because the photons causing absorption and stimulated emission are in focused beam modes. This leads to a dependence of the optically induced changes in the internal states on the center of mass atomic position. In the limit where the time between pulses is short, the decoherence can be expressed as a simple analytic formula involving the laser waist, temperature of the atoms, the trap frequency and the atomic mass. These limits on gate fidelity are studied for the standard $\pi-2\pi-\pi$ Rydberg gate and a new protocol based on a single adiabatic pulse with Gaussian envelope., Comment: 8 figures

The interference between Compton scattering and nonlinear Compton scattering from a two-color field in the X-ray regime is theoretically examined for bound electrons. The underlying phase shifts are analysed using a perturbative approach in the incoming classical field. The perturbative approach is bench marked with a non-perturbative approach in the classical field. The interference for different combinations of linear polarization of the two fields is examined when the Compton and the nonlinear Compton scattered waves have the same wave vector and polarization. Only two cases exhibit interference. When there is interference, the calculations reveal an intrinsic phase difference between the Compton scattered wave function and the nonlinear Compton scattered wave function of either 0 or $\pi$ depending on the scattering angle., Comment: 11 pages, 6 figures

Forces and torques exerted on dielectric disks trapped in a Gaussian standing wave are analyzed theoretically for disks of radius $2~\mu\text{m}$ with index of refraction $n=1.45$ and $n=2.0$ as well as disks of radius 200 nm with $n=1.45$. Calculations of the forces and torques were conducted both analytically and numerically using a discrete-dipole approximation method. Besides harmonic terms, third order ro-translational coupling terms in the potential energy can be significant and a necessary consideration when describing the dynamics of disks outside of the Rayleigh limit. The coupling terms are a result of the finite extension of the disk coupling to both the Gaussian and standing wave geometry of the beam. The resulting dynamics of the degrees of freedom most affected by the coupling terms exhibit several sidebands as evidenced in the power spectral densities. Simulations show that for Gaussian beam waists of $1-4~\mu\text{m}$ the disk remains stably trapped., Comment: 9 pages, 5 figures

We predict the existence of a universal class of ultralong-range Rydberg molecular states whose vibrational spectra form trimmed Rydberg series. A dressed ion-pair model captures the physical origin of these exotic molecules, accurately predicts their properties, and reveals features of ultralong-range Rydberg molecules and heavy Rydberg states with a surprisingly small Rydberg constant. The latter is determined by the small effective charge of the dressed anion, which outweighs the contribution of the molecule's large reduced mass. This renders these molecules the only known few-body systems to have a Rydberg constant smaller than $R_\infty/2$., Comment: 6 pages, 3 figures and supplemental material (4 pages and 4 figures)

A detailed theoretical framework for highly excited Rydberg molecules is developed based on the generalized local frame transformation. Our approach avoids the use of pseudopotentials and yields analytical expressions for the body-frame reaction matrix. The latter is used to obtain the molecular potential energy curves, but equally it can be employed for photodissociation, photoionization, or other processes. To illustrate the reliability and accuracy of our treatment we consider the Rb$^*-$Rb Rydberg molecule and compare our treatment with state-of-the-art alternative approaches. As a second application, the present formalism is used to re-analyze the vibrational spectra of Sr$^*-$Sr molecules, providing additional physical insight into their properties and a comparison of our results with corresponding measurements., Comment: 13 pages, 8 figures and 1 table

A proposal for cooling the translational motion of optically levitated magnetic nanoparticles is presented. The theoretical cooling scheme involves the sympathetic cooling of a ferromagnetic YIG nanosphere with a spin-polarized atomic gas. Particle-atom cloud coupling is mediated through the magnetic dipole-dipole interaction. When the particle and atom oscillations are small compared to their separation, the interaction potential becomes dominantly linear which allows the particle to exchange energy with the $N$ atoms. While the atoms are continuously Doppler cooled, energy is able to be removed from the nanoparticle's motion as it exchanges energy with the atoms. The rate at which energy is removed from the nanoparticle's motion was studied for three species of atoms (Dy, Cr, Rb) by simulating the full $N+1$ equations of motion and was found to depend on system parameters with scalings that are consistent with a simplified model. The nanoparticle's damping rate due to sympathetic cooling is competitive with and has the potential to exceed commonly employed cooling methods., Comment: 8 pages, 2 figures

The results of simulations of cooling based on Sawtooth-Wave Adiabatic Passage (SWAP) are presented including the possibility of population leaking to states outside of the cycling transition. The amount of population leaking can be substantially suppressed compared to Doppler cooling, which could be useful for systems that are difficult to repump back to the cycling transition. The suppression of the leaked population was more effective when simulating the slowing of a beam than in cooling a thermal distribution. As expected, calculations of the leaked population versus branching ratio of spontaneous emission show that the suppression is more effective for narrow linewidth transitions. In this limit, using SWAP to slow a beam may be worth pursuing even when the branching ratio out of the cycling transition is greater than 10%.

This paper quantifies the rate at which laser shot noise energy is delivered to a nanoparticle for the various scenarios commonly encountered in levitated optomechanics. While previous articles have the same form and dependencies, the proportionality constants often differ in the literature. This paper resolves these discrepancies. The rate at which energy is delivered to an optically trapped particle's respective degrees of freedom depends on the radiation pattern of scattered light as well as the direction of laser propagation. For a traveling plane wave with linearly polarized light, in the Rayleigh regime this leads the translational shot noise heating rate to be proportional to $1/10$ of the total rate in the laser polarization direction, $7/10$ in the laser propagation direction, and $2/10$ in the direction perpendicular to both. Analytical expressions for the shot noise heating rate are provided in the Rayleigh limit as well as numerical calculations for particles in the Mie regime for silica and diamond. For completeness, numerical calculations of the shot noise heating for silica Mie particles at the focus of a strongly focused laser beam are calculated for varying numerical aperture and common laser wavelengths. Both numerical calculations show that the Rayleigh expression generally gives an overestimate of the shot noise heating especially for larger radii, but is still a good approximation even for incident focal fields. The exception to the relative decrease is when a Mie resonance is reached which was found for diamond. Lastly, Rayleigh expressions for the rotational shot noise heating for a symmetric top-like particle for linear, elliptically, and unpolarized light are also provided., Comment: 11 pages, 3 figures

We study the effect of a weakly driven atomic cloud's polarization distribution on its photon scattering lineshape. In doing this, we find three distinct polarization regimes. First, for dilute clouds, the polarization magnitude is relatively constant. Second, for denser clouds, polarization builds at the front of the cloud for near-resonant light. Third, when the cloud condenses to the point where its dimensions become comparable to the wavelength, light refocuses towards the back of the cloud for red detuning. For these regimes, we show which `mean-field' frameworks accurately describe the differing photon scattering lineshapes. Finally, for even denser clouds, mean field models become inaccurate and necessitate the full point dipole model that includes atom-atom correlations.

Antihydrogen production has reached such a level that precision spectroscopic measurements of its properties are within reach. In particular, the ground-state level population is of central interest for experiments aiming at antihydrogen spectroscopy. The positron density and temperature dependence of the ground-state yield is a result of the interplay between recombination, collisional, and radiative processes. Considering the fact that antihydrogen atoms with the principal quantum number n=15 or lower quickly cascade down to the ground state within 1ms, the number of such states are adopted as a measure of useful antihydrogen atoms. It has been found that the scaling behaviour of the useful antihydrogen yield is different depending on the positron density and positron temperature., Comment: 6 pages, 6 figures

The properties of antihydrogen ($\bar{\rm H}$) have, thus far, been probed at magnetic fields of $\sim 1$ T. It may be fruitful to perform some of these measurements at magnetic fields approaching 0 T. In this case, there could occur zeros in the magnitude of the $B$-field. The number and properties of the magnetic field zeros are investigated. For typical magnetic field geometries in $\bar{\rm H}$ traps, the zeros will occur as two groups of 5 closely spaced points instead of as a single point. Except in special cases, results from calculations show that these 10 zeros can be treated as independent sources of spin flip probability. Although the behavior of Majorana spin flip near higher order zeros should not be important in the $\bar{\rm H}$ traps, the probability for spin flip is calculated for the case of a quadratic zero. Finally, results are presented for a simple model of how magnetic field zeros would affect the trapped population of $\bar{\rm H}$.

We theoretically investigate the rigid body dynamics of an optically levitated nanodumbbell under parametric feedback cooling and provide a simplified model for describing the motion. Differing from previous studies, the spin of the nanoparticle about its symmetry axis is considered non-negligible. Simulations reveal that standard parametric feedback cooling can extract energy from two of the five rotational degrees of freedom when the nanoparticle is levitated using a linearly polarized laser beam. The dynamics after feedback cooling are characterized by a normal mode describing precession about the laser polarization axis together with spin about the nanoparticle's symmetry axis. Cooling the remaining mode requires an asymmetry in the two librational frequencies associated with motion about the polarization axis as well as information about the two frequencies of rotation about the polarization axis. Introducing an asymmetric potential allows full cooling of the librational coordinates if the frequencies of both are used in the feedback modulation and is an avenue for entering the librational quantum regime. The asymmetry in the potential needs to be large enough for practical cooling times as the cooling rate of the system depends non-linearly on the degree of asymmetry, a condition that is easily achieved experimentally., Comment: 11 pages, 5 figures

We have experimentally demonstrated both heating and cooling of electrons in an ultracold neutral plasma (UNP) by embedding Rydberg atoms into the plasma soon after its creation. We have determined the relationship between the initial electron temperature, $T_{e,i}$, and the binding energy of the added Rydberg atoms, $E_b$, at the crossover between heating and cooling behaviors (that is, the binding energy of the atoms, which, when they are added to the plasma, neither accelerate or slow down the plasma expansion). Specifically, this condition is $|E_b| \approx 2.7 \times k_B T_{e,i}$ when the diagnostic used is the effect of the Rydberg atoms on the plasma asymptotic expansion velocity. Additionally, we have obtained experimental estimates for the amount of heating or cooling which occurs when the Rydberg binding energy does not satisfy the crossover condition. The experimental results for the crossover condition, and the degree of heating or cooling away from the crossover, are in agreement with predictions obtained from numerical modeling of the interactions between Rydberg atoms and the plasma. We have also developed a simple intuitive picture of how the Rydberg atoms affect the plasma which supports the concept of a `bottleneck' in the Rydberg state distribution of atoms in equilibrium with a co-existing plasma., Comment: 15 pages, 5 figures

The theory is developed for one and two atom interactions when the atom has a Rydberg electron attached to a hyperfine split core state. This situation is relevant for some of the rare earth and alkaline earth atoms that have been proposed for experiments on Rydberg-Rydberg interactions. For the rare earth atoms, the core electrons can have a very substantial total angular momentum, $J$, and a non-zero nuclear spin, $I$. In the alkaline earth atoms there is a single, $s$, core electron whose spin can couple to a non-zero nuclear spin for odd isotopes. The resulting hyperfine splitting of the core state can lead to substantial mixing between the Rydberg series attached to different thresholds. Compared to the unperturbed Rydberg series of the alkali atoms, the series perturbations and near degeneracies from the different parity states could lead to qualitatively different behavior for single atom Rydberg properties (polarizability, Zeeman mixing and splitting, etc) as well as Rydberg-Rydberg interactions ($C_5$ and $C_6$ matrices).

In this work, we demonstrate the importance of considering correlations between degenerate Zeeman sublevels that develop in dense atomic ensembles. In order to do this, we develop a set of equations capable of simulating large numbers of atoms while still incorporating correlations between degenerate Zeeman sublevels. This set of equations is exact in the single-photon limit, and may be interpreted as a generalization of the coupled harmonic oscillator equations typically used the literature. Using these equations, we demonstrate that in sufficiently dense systems, correlations between Zeeman sublevels can cause non-trivial differences in the photon scattering lineshape in arrays and clouds of atoms.

In this paper, results of quantum calculations are presented for feedback cooling of an optically trapped nanoparticle in the laser-shot-noise-dominant regime. We numerically investigate the system using both parametric and force feedback cooling schemes. For the same measurement efficiency, the cooling limit from the force feedback is lower than that from the parametric feedback. We also develop a set of semi-classical equations for feedback cooling that accurately match the quantum results. It is demonstrated, by rescaling the semi-classical equations, that the cooling dynamics is uniquely determined by the parameter set: the feedback strength, the measurement efficiency and the change of occupation number over one oscillation period due to the shot noise. The minimum occupation number is determined by the measurement efficiency and the change of occupation number over one period.

Results on the heating and the parametric feedback cooling of an optically trapped anisotropic nanoparticle in the laser shot noise dominant regime are presented. The related dynamical parameters, such as the oscillating frequency and shot noise heating rate, depend on the shape of the trapped particle. For an ellipsoidal particle, the ratio of the axis lengths and the overall size controls the shot noise heating rate relative to the frequency. For a particle with smaller ellipticity or bigger size, the relative heating rate for rotation tends to be smaller than that for translation indicating a better rotational cooling. For one feedback scheme, we also present results on the lowest occupation number that can be achieved as a function of the heating rate and the amount of classical uncertainty in the position measurement.

We observe that high-Q electromagnetic cavity resonances increase the cyclotron cooling rate of pure electron plasmas held in a Penning-Malmberg trap when the electron cyclotron frequency, controlled by tuning the magnetic field, matches the frequency of standing wave modes in the cavity. For certain modes and trapping configurations, this can increase the cooling rate by factors of ten or more. In this paper, we investigate the variation of the cooling rate and equilibrium plasma temperatures over a wide range of parameters, including the plasma density, plasma position, electron number, and magnetic field., Comment: 5 pages, 5 figures

In this comment, we agree with the formulas derived in Refs. [1,2] but show that the results are not due to "interference between spatially separated states of the center of mass of a mesoscopic harmonic oscillator".

An optically levitated nanoparticle in vacuum is a paradigm optomechanical system for sensing and studying macroscopic quantum mechanics. While its center-of-mass motion has been investigated intensively, its torsional vibration has only been studied theoretically in limited cases. Here we report the first experimental observation of the torsional vibration of an optically levitated nonspherical nanoparticle in vacuum. We achieve this by utilizing the coupling between the spin angular momentum of photons and the torsional vibration of a nonspherical nanoparticle whose polarizability is a tensor. The torsional vibration frequency can be one order of magnitude higher than its center-of-mass motion frequency, which is promising for ground state cooling. We propose a simple yet novel scheme to achieve ground state cooling of its torsional vibration with a linearly-polarized Gaussian cavity mode. A levitated nonspherical nanoparticle in vacuum will also be an ultrasensitive nanoscale torsion balance with a torque detection sensitivity on the order of $10^{-29} ~\mathrm{N}\cdot \mathrm{m}/\sqrt{\mathrm{ Hz}}$ under realistic conditions.

The mechanism of decoherence for a quantum system with rotational degrees of freedom is studied. From a simple model of elastic scattering, we show that the non-diagonal density matrix elements of the system exponentially decay. The decay rate depends on the difference of scattering amplitudes for different rotational configurations, leading to the gradual loss of quantum coherence between the pointer states in the system orientational space. For a dielectric ellipsoid immersed in a photon-gas environment (assuming no absorption), the decay rate is found to be proportional to the seventh power of the temperature. For an ellipsoidal object interacting with mass particles, the decay rate is proportional to the 5/2 power of the temperature. Both are different from the case of translational decoherence induced by the same environment scattering. For photon scattering, the coherence time in the rotational degrees of freedom is shown to be much shorter than that in the translational degrees of freedom.

The coherent dipole-dipole interactions of atoms in an atomic array are studied. It is found that the excitation probability of an atom in an array parallel to the direction of laser propagation ($\boldsymbol{\hat{k}}$) will either grow or decay logarithmically along $\boldsymbol{\hat{k}}$, depending on the detuning of the laser. The symmetry of the system for atomic separations of $\delta r = j\lambda/2$, where $j$ is an integer, causes the excitation distribution and scattered radiation to abruptly become symmetric about the center of the array. For atomic separations of $\delta r < \lambda/2$, the appearance of a collection of extremely subradiant states ($\Gamma\sim 0$), disrupts the described trend. In order to interpret the results from a finite array of atoms, a band structure calculation in the $N\rightarrow \infty$ limit is conducted where the decay rates and the Collective Lamb Shifts of the eigenmodes along the Brillouin zone are shown. Finally, the band structure of an array strongly affects its scattered radiation, allowing one to manipulate the Collective Lamb Shift as well as the decay rate (from superradiant to subradiant) by changing the angle of the driving laser.

The standard closed-orbit theory is extended for the photodetachment of negative ions in a time-dependent electric field. The time-dependent photodetachment rate is specifically studied in the presence of a single-cycle terahertz pulse, based on exact quantum simulations and semiclassical analysis. We find that the photodetachment rate is unaffected by a weak terahertz field, but oscillates complicatedly when the terahertz pulse gets strong enough. Three types of closed classical orbits are identified for the photoelectron motion in a strong single-cycle terahertz pulse, and their connections with the oscillatory photodetachment rate are established quantitatively by generalizing the standard closed-orbit theory to a time-dependent form. By comparing the negative hydrogen and fluorine ions, both the in-phase and antiphase oscillations can be observed, depending on a simple geometry of the contributed closed classical orbits. On account of its generality, the presented theory provides an intuitive understanding from a time-dependent viewpoint for the photodetachment dynamics driven by an external electric field oscillating at low frequency.

A rigorous theoretical framework is developed for a generalized local frame transformation theory (GLFT). The GLFT is applicable to the following systems: to Rydberg atoms or molecules in an electric field, or to negative ions in any combination of electric and/or magnetic fields. A first test application to the photoionization spectra of Rydberg atoms in an external electric field demonstrates dramatic improvement over the first version of the local frame transformation theory developed initially by Fano and Harmin. This revised GLFT theory yields non-trivial corrections because it now includes the full on-shell Hilbert space without adopting the truncations in the original theory. Comparisons of the semi-analytical GLFT Stark spectra with {\it ab initio} numerical simulations yields errors in the range of a few tens of MHz, an improvement over the original Fano-Harmin theory whose errors are 10-100 times larger. Our analysis provides a systematic pathway to precisely describe the corresponding photoabsorption spectra that should be accurate enough to meet most modern experimental standards., Comment: 9 pages, 3 figure

Sufficiently high densities in Bose-Einstein condensates provide favorable conditions for the production of ultralong-range polyatomic molecules consisting of one Rydberg atom and a number of neutral ground state atoms. The chemical binding properties and electronic wave functions of these exotic molecules are investigated analytically via hybridized diatomic states. The effects of the molecular geometry on the system's properties are studied through comparisons of the adiabatic potential curves and electronic structures for both symmetric and randomly configured molecular geometries. General properties of these molecules with increasing numbers of constituent atoms and in different geometries are presented. These polyatomic states have spectral signatures that lead to non-Lorentzian line-profiles., Comment: 13 pages, 8 figures