Source author record

Charles W. Clark

Charles W. Clark appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

63works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

63 published item(s)

preprint2022arXiv

Experimental Realization of Neutron Helical Waves

Methods of preparation and analysis of structured waves of light, electrons, and atoms have been advancing rapidly. Despite the proven power of neutrons for material characterization and studies of fundamental physics, neutron science has not been able to fully integrate such techniques due to small transverse coherence lengths, the relatively poor resolution of spatial detectors, and low fluence rates. Here, we demonstrate methods that are practical with the existing technologies, and show the experimental achievement of neutron helical wavefronts that carry well-defined orbital angular momentum (OAM) values. We discuss possible applications and extensions to spin-orbit correlations and material characterization techniques.

preprint2020arXiv

Producing flow in "racetrack" atom circuits by stirring at zero and non-zero temperature

We present a study of how macroscopic flow can be produced in Bose-Einstein condensate confined in a "racetrack" potential by stirring with a wide rectangular barrier. This potential consists of two half-circle channels separated by straight channels of length $L$ and is a ring potential if $L=0$. We present the results of a large set of simulations where racetrack condensates stirred with a barrier under varying conditions of barrier height, stir speed, racetrack geometry, and temperature. We found that stirring was readily able to produce circulation in ring and non-ring geometries but that the exact amount of flow produced was complicated. We therefore also studied the mechanism by which flow was produced in the stirring process. We found that circulation was induced by the swap of a vortex/anti-vortex pair that was initially created by backflow of the condensate in the region of depressed density by the barrier. When the barrier strength reached a critical value a number of these vortex-swap events occurred in rapid succession so that flow speed of the circulation produced was enough to exceed the stir speed of the barrier. Flow that was initially localized in the vortices involved in the vortex swap was converted into macroscopic flow around the racetrack by pairs of disturbances each generated during the vortex swap. Each pair consisted of a vortex/anti-vortex pair moving in the anti-stir direction and a compression wave moving in the stir direction. This picture of the mechanism for making flow will enable the design of stirring schedules that create a desired amount of flow.

preprint2019arXiv

Generation and detection of spin-orbit coupled neutron beams

Spin-orbit coupling of light has come to the fore in nano-optics and plasmonics, and is a key ingredient of topological photonics and chiral quantum optics. We demonstrate a basic tool for incorporating analogous effects into neutron optics: the generation and detection of neutron beams with coupled spin and orbital angular momentum. $^3$He neutron spin-filters are used in conjunction with specifically oriented triangular coils to prepare neutron beams with lattices of spin-orbit correlations, as demonstrated by their spin-dependant intensity profiles. These correlations can be tailored to particular applications, such as neutron studies of topological materials.

preprint2016arXiv

Entanglement and spin-squeezing without infinite-range interactions

Infinite-range interactions are known to facilitate the production of highly entangled states with applications in quantum information and metrology. However, many experimental systems have interactions that decay with distance, and the achievable benefits in this context are much less clear. Combining recent exact solutions with a controlled expansion in the system size, we analyze quench dynamics in Ising models with power-law ($1/r^α$) interactions in $D$ dimensions, thereby expanding the understanding of spin squeezing into a broad and experimentally relevant context. In spatially homogeneous systems, we show that for small $α$ the scaling of squeezing with system size is identical to the infinite-range ($α=0$) case. This indifference to the interaction range persists up to a critical value $α=2D/3$, above which squeezing degrades continuously. Boundary-induced inhomogeneities present in most experimental systems modify this picture, but it nevertheless remains qualitatively correct for finite-sized systems.

preprint2016arXiv

Magic wavelengths, matrix elements, polarizabilities, and lifetimes of Cs

Motivated by recent interest in their applications, we report a systematic study of Cs atomic properties calculated by a high-precision relativistic all-order method. Excitation energies, reduced matrix elements, transition rates, and lifetimes are determined for levels with principal quantum numbers $n \leq 12$ and orbital angular momentum quantum numbers $l \leq 3$. Recommended values and estimates of uncertainties are provided for a number of electric-dipole transitions and the electric dipole polarizabilities of the $ns$, $np$, and $nd$ states. We also report a calculation of the electric quadrupole polarizability of the ground state. We display the dynamic polarizabilities of the $6s$ and $7p$ states for optical wavelengths between 1160 nm and 1800 nm and identify corresponding magic wavelengths for the $6s-7p_{1/2}$, $6s-7p_{3/2}$ transitions. The values of relevant matrix elements needed for polarizability calculations at other wavelengths are provided.

preprint2016arXiv

Spin-Orbit States of Neutron Wavepackets

We propose a method to prepare an entangled spin-orbit state between the spin and the orbital angular momenta of a neutron wavepacket. This spin-orbit state is created by passing neutrons through the center of a quadrupole magnetic field, which provides a coupling between the spin and orbital degrees of freedom. A Ramsey fringe type measurement is suggested as a means of verifying the spin-orbit correlations.

preprint2016arXiv

Superfluid transport dynamics in a capacitive atomtronic circuit

We simulate transport in an atomtronic circuit of a Bose-Einstein condensate that flows from a source region into a drain through a gate channel. The time-dependent Gross-Pitaevskii equation (GPE) solution matches the data of a recent experiment. The atomtronic circuit is found to be similar to a variable-resistance RLC circuit, which is critically damped at early times and shows LC oscillations later. The GPE also predicts atom loss from the drain. Studies of the dependence of condensate transport upon gate parameters suggest the utility of the GPE for investigation of atomtronic circuits.

preprint2015arXiv

Resonant wavepackets and shock waves in an atomtronic SQUID

The fundamental dynamics of ultracold atomtronic devices are reflected in their phonon modes of excitation. We probe such a spectrum by applying a harmonically driven potential barrier to a $^{23}$Na Bose-Einstein condensate in a ring-shaped trap. This perturbation excites phonon wavepackets. When excited resonantly, these wavepackets display a regular periodic structure. The resonant frequencies depend upon the particular configuration of the barrier, but are commensurate with the orbital frequency of a Bogoliubov sound wave traveling around the ring. Energy transfer to the condensate over many cycles of the periodic wavepacket motion causes enhanced atom loss from the trap at resonant frequencies. Solutions of the time-dependent Gross-Pitaevskii equation exhibit quantitative agreement with the experimental data. We also observe the generation of supersonic shock waves under conditions of strong excitation, and collisions of two shock wavepackets.

preprint2015arXiv

Spinor Bose-Einstein Condensates of Positronium

Bose-Einstein condensates (BECs) of positronium (Ps) have been of experimental and theoretical interest due to their potential application as the gain medium of a $γ$-ray laser. Ps BECs are intrinsically spinor due to the presence of ortho-positronium (o-Ps) and para-positronium (p-Ps), whose annihilation lifetimes differ by three orders of magnitude. In this paper, we study the spinor dynamics and annihilation processes in the p-Ps/o-Ps system using both solutions of the time-dependent Gross-Pitaevskii equations and a semiclassical rate-equation approach. The spinor interactions have an $O(4)$ symmetry which is broken to $SO(3)$ by an internal energy difference between o-Ps and p-Ps. For an initially unpolarized condensate, there is a threshold density of $\approx 10^{19}$ cm$^{-3}$ at which spin mixing between o-Ps and p-Ps occurs. Beyond this threshold, there are unstable spatial modes accompanied by spin mixing. To ensure a high production yield above the critical density, a careful choice of external field must be made to avoid the spin mixing instability.

preprint2015arXiv

The magic road to precision

We predict a sequence of magic-zero wavelengths for the Sr excited $5s5p~ ^3P_0$ state, and provide a general roadmap for extracting transition matrix elements using precise frequency measurements. We demonstrate that such measurements can serve as a best global benchmark of the spectroscopic accuracy that is required for the development of high-precision predictive methods. These magic-zero wavelengths are also needed for state-selective atom manipulation for implementation of quantum logic operations. We also identify five magic wavelengths of the $5s^2\ ^1S_0 - 5s5p\ ^3P_0$ Sr clock transition between 350 nm and 500 nm which can also serve as precision benchmarks.

preprint2014arXiv

Decay of a superfluid current of ultra-cold atoms in a toroidal trap

Using a numerical implementation of the truncated Wigner approximation, we simulate the experiment reported by Ramanathan et al. in Phys. Rev. Lett. 106, 130401 (2011), in which a Bose-Einstein condensate is created in a toroidal trap and set into rotation via a phase imprinting technique. A potential barrier is then placed in the trap to study the decay of the superflow. We find that the current decays via thermally activated phase slips, which can also be visualized as vortices crossing the barrier region in the radial direction. Adopting the notion of critical velocity used in the experiment, we determine it to be lower than the local speed of sound at the barrier, in contradiction to the predictions of the zero-temperature Gross-Pitaevskii equation. We map out the superfluid decay rate and critical velocity as a function of temperature and observe a strong dependence. Thermal fluctuations offer a partial explanation of the experimentally observed reduction of the critical velocity from the phonon velocity.

preprint2014arXiv

Half-Quantum Vortex Molecules in a Binary Dipolar Bose Gas

We study the ground state phases of a rotating two-component, or binary Bose-Einstein condensate, wherein one component possesses a large magnetic dipole moment. A variety of non-trivial phases emerge in this system, including a half-quantum vortex (HQV) chain phase and a HQV molecule phase, where HQVs of opposite charge bind at short distances. We attribute the emergence of these phases to the development of a minimum in the adiabatic HQV interaction potential, which we calculate explicitly. We thus show that the presence of dipolar interactions in this system leads to a rich phase diagram, and the formation of HQV molecules.

preprint2014arXiv

High-resolution, vacuum-ultraviolet absorption spectrum of boron trifluoride

In the course of investigations of thermal neutron detection based on mixtures of $^{10}$BF$_3$ with other gases, knowledge was required of the photoabsorption cross sections of $^{10}$BF$_3$ for wavelengths between 135 and 205 nm. Large discrepancies in the values reported in existing literature led to the absolute measurements reported in this communication. The measurements were made at the SURF III synchrotron radiation facility at the National Institute of Standards and Technology. The measured absorption cross sections vary from 10$^{-20}$ cm$^2$ at 135 nm to less than 10$^{-21}$ cm$^2$ in the region from 165 to 205 nm. Three previously unreported absorption features with resolvable structure were found in the regions 135 to 145 nm, 150 to 165 nm and 190 to 205 nm. Quantum mechanical calculations, using the TD-B3LYP/aug-cc-pVDZ variant of time-dependent density functional theory implemented in Gaussian 09, suggest that the observed absorption features arise from symmetry-changing adiabatic transitions.

preprint2014arXiv

Hysteresis in a quantized, superfluid atomtronic circuit

Atomtronics is an emerging interdisciplinary field that seeks new functionality by creating devices and circuits where ultra-cold atoms, often superfluids, play a role analogous to the electrons in electronics. Hysteresis is widely used in electronic circuits, e.g., it is routinely observed in superconducting circuits and is essential in rf-superconducting quantum interference devices [SQUIDs]. Furthermore, hysteresis is as fundamental to superfluidity (and superconductivity) as quantized persistent currents, critical velocity, and Josephson effects. Nevertheless, in spite of multiple theoretical predictions, hysteresis has not been previously observed in any superfluid, atomic-gas Bose-Einstein condensate (BEC). Here we demonstrate hysteresis in a quantized atomtronic circuit: a ring of superfluid BEC obstructed by a rotating weak link. We directly detect hysteresis between quantized circulation states, in contrast to superfluid liquid helium experiments that observed hysteresis directly in systems where the quantization of flow could not be observed and indirectly in systems that showed quantized flow. Our techniques allow us to tune the size of the hysteresis loop and to consider the fundamental excitations that accompany hysteresis. The results suggest that the relevant excitations involved in hysteresis are vortices and indicate that dissipation plays an important role in the dynamics. Controlled hysteresis in atomtronic circuits may prove to be a crucial feature for the development of practical devices, just as it has in electronic circuits like memory, digital noise filters (e.g., Schmitt triggers), and magnetometers (e.g., SQUIDs).

preprint2014arXiv

Nearly-linear light cones in long-range interacting quantum systems

In non-relativistic quantum theories with short-range Hamiltonians, a velocity $v$ can be chosen such that the influence of any local perturbation is approximately confined to within a distance $r$ until a time $t \sim r/v$, thereby defining a linear light cone and giving rise to an emergent notion of locality. In systems with power-law ($1/r^α$) interactions, when $α$ exceeds the dimension $D$, an analogous bound confines influences to within a distance $r$ only until a time $t\sim(α/v)\log r$, suggesting that the velocity, as calculated from the slope of the light cone, may grow exponentially in time. We rule out this possibility; light cones of power-law interacting systems are algebraic for $α>2D$, becoming linear as $α\rightarrow\infty$. Our results impose strong new constraints on the growth of correlations and the production of entangled states in a variety of rapidly emerging, long-range interacting atomic, molecular, and optical systems.

preprint2014arXiv

Noble gas excimer scintillation following neutron capture in boron thin films

Far-ultraviolet (FUV) scintillation signals have been measured in heavy noble gases (argon, krypton, xenon) following boron-neutron capture ($^{10}$B($n,α$)$^7$Li) in $^{10}$B thin films. The observed scintillation yields are comparable to the yields from some liquid and solid neutron scintillators. At noble gas pressures of 107 kPa, the number of photons produced per neutron absorbed following irradiation of a 1200 nm thick $^{10}$B film was 14,000 for xenon, 11,000 for krypton, and 6000 for argon. The absolute scintillation yields from the experimental configuration were calculated using data from (1) experimental irradiations, (2) thin-film characterizations, (3) photomultiplier tube calibrations, and (4) photon collection modeling. Both the boron films and the photomultiplier tube were characterized at the National Institute of Standards and Technology. Monte Carlo modeling of the reaction cell provided estimates of the photon collection efficiency and the transport behavior of $^{10}$B($n,α$)$^7$Li reaction products escaping the thin films. Scintillation yields increased with gas pressure due to increased ionization and excitation densities of the gases from the $^{10}$B($n,α$)$^7$Li reaction products, increased frequency of three-body, excimer-forming collisions, and reduced photon emission volumes (i.e., larger solid angle) at higher pressures. Yields decreased for thicker $^{10}$B thin films due to higher average energy loss of the $^{10}$B($n,α$)$^7$Li reaction products escaping the films. The relative standard uncertainties in the measurements were determined to lie between 14 % and 16 %. The observed scintillation signal demonstrates that noble gas excimer scintillation is promising for use in practical neutron detectors.

preprint2014arXiv

Relativistic all-order calculations of Th, Th$^{+}$ and Th$^{2+}$ atomic properties

Excitation energies, term designations, and $g$-factors of Th, Th$^{+}$ and Th$^{2+}$ are determined using a relativistic hybrid configuration interaction (CI) + all-order approach that combines configuration interaction and linearized coupled-cluster methods. The results are compared with other theory and experiment where available. We find some "vanishing" $g$-factors, similar to those known in lanthanide spectra. Reduced matrix elements, oscillator strengths, transition rates, and lifetimes are determined for Th$^{2+}$. To estimate the uncertainties of our results, we compared our values with the available experimental lifetimes for higher $5f7p\ ^3G_{4}$, $7s7p\ ^3P_{0}$, $7s7p\ ^3P_{1}$, and $6d7p\ ^3F_{4}$ levels of Th$^{2+}$. These calculations provide a benchmark test of the CI+all-order method for heavy systems with several valence electrons and yield recommended values for transition rates and lifetimes of Th$^{2+}$.

preprint2014arXiv

Relativistic calculations of $C_6$ and $C_8$ coefficients for strontium dimers

The electric dipole and quadrupole polarizabilities of the $5s5p~^3\!P_1^o$ state and the $C_6$ and $C_8$ coefficients for the $^1\!S_0 +\, ^1\!S_0$ and $^1\!S_0 +\, ^3\!P_1^o$ dimers of strontium are calculated using a high-precision relativistic approach that combines configuration interaction and linearized coupled-cluster methods. Our recommended values of the long range dispersion coefficients for the $0_u$ and $1_u$ energy levels are $C_6(0_u)=3771(32)$ a.u. and $C_6(1_u)= 4001(33)$ a.u., respectively. They are in good agreement with recent results from experimental photoassociation data. We also calculate $C_8$ coefficients for Sr dimers, which are needed for precise determination of long-range interaction potential. We confirm the experimental value for the magic wavelength, where the Stark shift on the $^1\!S_0$-$^3\!P_1^o$ transition vanishes. The accuracy of calculations is analyzed and uncertainties are assigned to all quantities reported in this work.

preprint2014arXiv

Relativistic many-body calculations of van der Waals coefficients for Yb-Li and Yb-Rb dimers

We derive the relativistic formulas for the van der Waals coefficients of Yb-alkali dimers that correlate to ground and excited separated-atom limits. We calculate $C_6$ and $C_8$ coefficients of particular experimental interest. We also derive a semi-empirical formula that expresses the $C_8$ coefficient of heteronuclear $A+B$ dimers in terms of the $C_6$ and $C_8$ coefficients of homonuclear dimers and the static dipole and quadrupole polarizabilities of the atomic states $A$ and $B$. We report results of calculation of the $C_6$ coefficients for the Yb-Rb $^3/!P_1^o+5s\, ^2/!S_{1/2}$ and $^1/!S_0+5p\, ^2/!P^o_{1/2}$ dimers, and the $C_8$ coefficients for the Yb-Li $^1/!S_0+2s\, ^2/!S_{1/2}$ and Yb-Rb $^1/!S_0+5s\, ^2/!S_{1/2}$ dimers. Uncertainties are estimated for all predicted properties.

preprint2014arXiv

Robust finite-temperature disordered Mott insulating phases in inhomogeneous Fermi-Fermi mixtures with density and mass imbalance

Ultracold mixtures of different atomic species have great promise for realizing novel many-body phenomena. In a binary mixture of femions with a large mass difference and repulsive interspecies interactions, a disordered Mott insulator phase can occur. This phase displays an incompressible total density, although the relative density remain compressible. We use strong-coupling and Monte Carlo calculations to show that this phase exists for a broad parameter region for ultracold gases confined in a harmonic trap on a three-dimensional optical lattice, for experimentally accessible values of the trap parameters.

preprint2013arXiv

Laser cooling and trapping of potassium at magic wavelengths

We carry out a systematic study of the static and dynamic polarizabilities of the potassium atom using a first-principles high-precision relativistic all-order method in which all single, double, and partial triple excitations of the Dirac-Fock wave functions are included to all orders of perturbation theory. Recommended values are provided for a large number of electric-dipole matrix elements. Static polarizabilities of the 4s, 4p_j, 5s, 5p_j, and 3d_j states are compared with other theory and experiment where available. We use the results of the polarizability calculations to identify magic wavelengths for the 4s-np transitions for $n = 4, 5$, i.e. those wavelengths for which the two levels have the same ac Stark shifts. These facilitate state-insensitive optical cooling and trapping. The magic wavelengths for the $4s-5p$ transitions are of particular interest for attaining a quantum gas of potassium at high phase-space density. We find 20 such wavelengths in the technically interest region of 1050-1130 nm. Uncertainties of all recommended values are estimated.

preprint2013arXiv

Long range interaction coefficients for ytterbium dimers

We evaluate the electric-dipole and electric-quadrupole static and dynamic polarizabilities for the 6s^2 ^1S_0, 6s6p ^3P_0, and 6s6p ^3P_1 states and estimate their uncertainties. A methodology is developed for an accurate evaluation of the van der Waals coefficients of dimers involving excited state atoms with strong decay channel to the ground state. This method is used for evaluation of the long range interaction coefficients of particular experimental interest, including the C_6 coefficients for the Yb-Yb ^1S_0+^3P_{0,1} and ^3P_0+^3P_0 dimers and C_8 coefficients for the ^1S_0+^1S_0 and ^1S_0+^3P_1 dimers.

preprint2013arXiv

Meron Ground State of Rashba Spin-Orbit-Coupled Dipolar Bosons

We study the effects of dipolar interactions on a Bose-Einstein condensate with synthetically generated Rashba spin-orbit coupling. The dipolar interaction we consider includes terms that couple spin and orbital angular momentum in a way perfectly congruent with the single-particle Rashba coupling. We show that this internal spin-orbit coupling plays a crucial role in the rich ground-state phase diagram of the trapped condensate. In particular, we predict the emergence of a thermodynamically stable ground state with a meron spin configuration.

preprint2013arXiv

Probing the circulation of ring-shaped Bose-Einstein condensates

This paper reports the results of a theoretical and experimental study of how the initial circulation of ring-shaped Bose-Einstein condensates (BECs) can be probed by time-of-flight (TOF) images. We have studied theoretically the dynamics of a BEC after release from a toroidal trap potential by solving the 3D Gross-Pitaevskii (GP) equation. The trap and condensate characteristics matched those of a recent experiment. The circulation, experimentally imparted to the condensate by stirring, was simulated theoretically by imprinting a linear azimuthal phase on the initial condensate wave function. The theoretical TOF images were in good agreement with the experimental data. We find that upon release the dynamics of the ring--shaped condensate proceeds in two distinct phases. First, the condensate expands rapidly inward, filling in the initial hole until it reaches a minimum radius that depends on the initial circulation. In the second phase, the density at the inner radius increases to a maximum after which the hole radius begins slowly to expand. During this second phase a series of concentric rings appears due to the interference of ingoing and outgoing matter waves from the inner radius. The results of the GP equation predict that the hole area is a quadratic function of the initial circulation when the condensate is released directly from the trap in which it was stirred and is a linear function of the circulation if the trap is relaxed before release. These scalings matched the data. Thus, hole size after TOF can be used as a reliable probe of initial condensate circulation. This connection between circulation and hole size after TOF will facilitate future studies of atomtronic systems that are implemented in ultracold quantum gases.

preprint2013arXiv

Spatial correlations of one dimensional driven-dissipative systems of Rydberg atoms

We consider a one-dimensional lattice of atoms with laser excitation to a Rydberg state and spontaneous emission. The atoms are coupled due to the dipole-dipole interaction of the Rydberg states. This driven-dissipative system has a broad range of non-equilibrium phases, such as antiferromagnetic ordering and bistability. Using the quantum trajectory method, we calculate the spatial correlation function throughout the parameter space for up to 20 lattice sites. We show that bistability significantly strengthens the spatial correlations and entanglement.

preprint2013arXiv

Spin Waves and Dielectric Softening of Polar Molecule Condensates

We consider an oblate Bose-Einstein condensate of heteronuclear polar molecules in a weak applied electric field. This system supports a rich quasiparticle spectrum that plays a critical role in determining its bulk dielectric properties. In particular, in sufficiently weak fields the system undergoes a polarization wave rotonization, leading to the development of textured electronic structure and a dielectric instability that is characteristic of the onset of a negative static dielectric function.

preprint2012arXiv

Approximate mean-field equations of motion for quasi-2D Bose-Einstein condensate systems

We present a method for approximating the solution of the three-dimensional, time-dependent Gross-Pitaevskii equation (GPE) for Bose-Einstein condensate systems where the confinement in one dimension is much tighter than in the other two. This method employs a hybrid Lagrangian variational technique whose trial wave function is the product of a completely unspecified function of the coordinates in the plane of weak confinement and a gaussian in the strongly confined direction having a time-dependent width and quadratic phase. The hybrid Lagrangian variational method produces equations of motion that consist of (1) a two-dimensional, effective GPE whose nonlinear coefficient contains the width of the gaussian and (2) an equation of motion for the width that depends on the integral of the fourth power of the solution of the 2D effective GPE. We apply this method to the dynamics of Bose-Einstein condensates confined in ring-shaped potentials and compare the approximate solution to the numerical solution of the full 3D GPE.

preprint2012arXiv

Blackbody radiation shift in the Sr optical atomic clock

We evaluated the static and dynamic polarizabilities of the 5s^2 ^1S_0 and 5s5p ^3P_0^o states of Sr using the high-precision relativistic configuration interaction + all-order method. Our calculation explains the discrepancy between the recent experimental 5s^2 ^1S_0 - 5s5p ^3P_0^o dc Stark shift measurement Δα= 247.374(7) a.u. [Middelmann et. al, arXiv:1208.2848 (2012)] and the earlier theoretical result of 261(4) a.u. [Porsev and Derevianko, Phys. Rev. A 74, 020502R (2006)]. Our present value of 247.5 a.u. is in excellent agreement with the experimental result. We also evaluated the dynamic correction to the BBR shift with 1 % uncertainty; -0.1492(16) Hz. The dynamic correction to the BBR shift is unusually large in the case of Sr (7 %) and it enters significantly into the uncertainty budget of the Sr optical lattice clock. We suggest future experiments that could further reduce the present uncertainties.

preprint2012arXiv

Bond order solid of two-dimensional dipolar fermions

Recent experimental realization of dipolar Fermi gases near or below quantum degeneracy provides opportunity to engineer Hubbard-like models with long range interactions. Motivated by these experiments, we chart out the theoretical phase diagram of interacting dipolar fermions on the square lattice at zero temperature and half filling. We show that in addition to p-wave superfluid and charge density wave order, two new and exotic types of bond order emerge generically in dipolar fermion systems. These phases feature homogeneous density but periodic modulations of the kinetic hopping energy between nearest or next-nearest neighbors. Similar, but manifestly different, phases of two-dimensional correlated electrons have previously only been hypothesized and termed "density waves of nonzero angular momentum". Our results suggest that these phases can be constructed flexibly with dipolar fermions, using currently available experimental techniques.

preprint2012arXiv

Magic wavelengths for optical cooling and trapping of lithium

Using first-principles calculations, we identify magic wavelengths for the 2s-2p and 2s-3p transitions in lithium. The ns and np atomic levels have the same ac Stark shifts at the corresponding magic wavelength, which facilitates state-insensitive optical cooling and trapping. Tune-out wavelengths for which the ground-state frequency-dependent polarizability vanishes are also calculated. Differences of these wavelengths between 6Li and 7Li are reported. Our approach uses high-precision, relativistic all-order methods in which all single, double, and partial triple excitations of the Dirac-Fock wave functions are included to all orders of perturbation theory. Recommended values are provided for a large number of Li electric-dipole matrix elements. Static polarizabilities for the 2s, 2p, 3s, 3p, and 3d levels are compared with other theory and experiment where available. Uncertainties of all recommended values are estimated. The magic wavelengths for the uv 2s-3p transition are of particular interest for the production of a quantum gas of lithium [Duarte et al., Phys. Rev. A 84, 061406R (2011)].

preprint2012arXiv

Noise Correlation Scalings: Revisiting the Quantum Phase Transition in Incommensurate Lattices with Hard-Core Bosons

Finite size scalings of the momentum distribution and noise correlations are performed to study Mott insulator, Bose glass, and superfluid quantum phases in hard-core bosons (HCBs) subjected to quasi-periodic disorder. The exponents of the correlation functions at the Superfluid to Bose glass (SF-BG) transition are found to be approximately one half of the ones that characterizes the superfluid phase. The derivatives of the peak intensities of the correlation functions with respect to quasiperiodic disorder are shown to diverge at the SF-BG critical point. This behavior does not occur in the corresponding free fermion system, which also exhibits an Anderson-like transition at the same critical point, and thus provides a unique experimental tool to locate the phase transition in interacting bosonic systems. We also report on the absence of primary sublattice peaks in the momentum distribution of the superfluid phase for special fillings.

preprint2012arXiv

Polarizabilities of Si^{2+}: a benchmark test of theory and experiment

We have calculated electric-dipole polarizabilities of the 3s^2 ^1S_0, 3s3p ^3P_0, and 3s3p ^1P_1 states of the Si^{2+} ion using recently developed configuration interaction + all-order method. Detailed evaluation of the uncertainties of the final results is carried out. Our value for the ground state electric-dipole polarizability 11.670(13) a.u. is in excellent agreement with the resonant excitation Stark ionization spectroscopy value 11.669(9) a.u. [Komara et al., J. Phys. B 38, 87 (2005); Mitroy, Phys. Rev. A 78, 052515 (2008)]. This work represents the most precise benchmark test to date of theory and experiment in divalent atoms. The near cancellation of the ns^2 ^1S_0 ground state and the lowest nsnp ^3P_0 polarizabilities previously observed in B+, Al+, In+, Tl+, and Pb^{2+} is also found in Si^{2+} ion.

preprint2012arXiv

Three-Dimensional Spin-Orbit Coupling in a Trap

We investigate the properties of an atom under the influence of a synthetic three-dimensional spin-orbit coupling (Weyl coupling) in the presence of a harmonic trap. The conservation of total angular momentum provides a numerically efficient scheme for finding the spectrum and eigenfunctions of the system. We show that at large spin-orbit coupling the system undergoes dimensional reduction from three to one dimension at low energies, and the spectrum is approximately Landau level-like. At high energies, the spectrum is approximately given by the three-dimensional isotropic harmonic oscillator. We explore the properties of the ground state in both position and momentum space. We find the ground state has spin textures with oscillations set by the spin-orbit length scale.

preprint2012arXiv

Unconventional Spin Density Waves in Dipolar Fermi Gases

The conventional spin density wave (SDW) phase (Overhauser, 1962), as found in antiferromagnetic metal for example (Fawcett 1988), can be described as a condensate of particle-hole pairs with zero angular momentum, $\ell=0$, analogous to a condensate of particle-particle pairs in conventional superconductors. While many unconventional superconductors with Cooper pairs of finite $\ell$ have been discovered, their counterparts, density waves with non-zero angular momenta, have only been hypothesized in two-dimensional electron systems (Nayak, 2000). Using an unbiased functional renormalization group analysis, we here show that spin-triplet particle-hole condensates with $\ell=1$ emerge generically in dipolar Fermi gases of atoms (Lu, Burdick, and Lev, 2012) or molecules (Ospelkaus et al., 2008; Wu et al.) on optical lattice. The order parameter of these exotic SDWs is a vector quantity in spin space, and, moreover, is defined on lattice bonds rather than on lattice sites. We determine the rich quantum phase diagram of dipolar fermions at half-filling as a function of the dipolar orientation, and discuss how these SDWs arise amidst competition with superfluid and charge density wave phases.

preprint2012arXiv

Ytterbium in quantum gases and atomic clocks: van der Waals interactions and blackbody shifts

We evaluated the C_6 coefficients of Yb-Yb and Yb -alkali/group II van der Waals interactions with 2% uncertainty. The only existing results for such quantities are for the Yb-Yb dimer. Our value, C_6=1929(39) a.u., is in excellent agreement with the recent experimental determination of 1932(35) a.u. We have also developed a new approach for the calculation of the dynamic correction to the blackbody radiation shift. We have calculated this quantity for the Yb 6s^2 ^1S_0 - 6s6p ^3P_0 clock transition with 3.5% uncertainty. This reduces the fractional uncertainty due to the blackbody radiation shift in the Yb optical clock at 300 K to 10^{-18} level.

preprint2011arXiv

Bright and Dark Solitons and Breathers in Strongly Repulsive Bose-Einstein Condensates

Collisional dynamics of solitary matter waves of superfluid hard core bosons, consisting of dark and bright non-linear waves as well as supersonic periodic trains, reveals remarkable richness and coherence, with the phase of the condensate playing a key role. Depending upon the condensate density and their relative velocity, two distinctive collisional types emerge: the intuitively expected repulsive collisions due to the hard core boson constraint; and, also collisions in which they "pass through" one another. In addition to confirming the soliton status of both bright and dark solitary waves, our studies reveal a variety of multi-solitons including multiple families of breathers, that can be produced and precisely controlled via quantum phase engineering.

preprint2011arXiv

Detecting paired and counterflow superfluidity via dipole oscillations

We suggest an experimentally feasible procedure to observe paired and counterflow superfluidity in ultra-cold atom systems. We study the time evolution of one-dimensional mixtures of bosonic atoms in an optical lattice following an abrupt displacement of an additional weak confining potential. We find that the dynamic responses of the paired superfluid phase for attractive inter-species interactions and the counterflow superfluid phase for repulsive interactions are qualitatively distinct and reflect the quasi long-range order that characterizes these states. These findings suggest a clear experimental procedure to detect these phases, and give an intuitive insight into their dynamics.

preprint2011arXiv

Generalized Thermalization in an Integrable Lattice System

After a quench, observables in an integrable system may not relax to the standard thermal values, but can relax to the ones predicted by the generalized Gibbs ensemble (GGE) [M. Rigol et al., Phys. Rev. Lett. 98, 050405 (2007)]. The GGE has been shown to accurately describe observables in various one-dimensional integrable systems, but the origin of its success is not fully understood. Here we introduce a microcanonical version of the GGE and provide a justification of the GGE based on a generalized interpretation of the eigenstate thermalization hypothesis, which was previously introduced to explain thermalization of nonintegrable systems. We study relaxation after a quench of one-dimensional hard-core bosons in an optical lattice. Exact numerical calculations for up to 10 particles on 50 lattice sites (~10^10 eigenstates) validate our approach.

preprint2011arXiv

Magic-zero wavelengths of alkali-metal atoms and their applications

Using first-principles calculations, we identify "magic-zero" optical wavelengths, λ_zero, for which the ground-state frequency-dependent polarizabilities of alkali-metal atoms vanish. Our approach uses high-precision, relativistic all-order methods in which all single, double, and partial triple excitations of the Dirac-Fock wave functions are included to all orders of perturbation theory. We discuss the use of magic-zero wavelengths for sympathetic cooling in two-species mixtures of alkalis with group-II and other elements of interest. Special cases in which these wavelengths coincide with strong resonance transitions in a target system are identified.

preprint2011arXiv

Precision calculation of blackbody radiation shifts for optical frequency metrology

We show that three group IIIB divalent ions, B+, Al+, and In+, have anomalously small blackbody radiation (BBR) shifts of the ns^2 1S0 - nsnp 3P0 clock transitions. The fractional BBR shifts for these ions are at least 10 times smaller than those of any other present or proposed optical frequency standards at the same temperature, and are less than 0.3% of the Sr clock shift. We have developed a hybrid configuration interaction + coupled-cluster method that provides accurate treatment of correlation corrections in such ions, considers all relevant states in the same systematic way, and yields a rigorous upper bound on the uncertainty of the final results. We reduce the BBR contribution to the fractional frequency uncertainty of the Al+ clock to 4 \times 10^{-19} at T=300K. We also reduce the uncertainties due to this effect at room temperature to 10^{-18} level for B+ and In+ to facilitate further development of these systems for metrology and quantum sensing. These uncertainties approach recent estimates of the feasible precision of currently proposed optical atomic clocks.

preprint2011arXiv

Prototyping method for Bragg-type atom interferometers

We present a method for rapid prototyping of new Bragg ultra-cold atom interferometer (AI) designs useful for assessing the performance of such interferometers. The method simulates the overall effect on the condensate wave function in a given AI design using two separate elements. These are (1) modeling the effect of a Bragg pulse on the wave function and (2) approximating the evolution of the wave function during the intervals between the pulses. The actual sequence of these pulses and intervals is then followed to determine the approximate final wave function from which the interference pattern can be calculated. The exact evolution between pulses is assumed to be governed by the Gross-Pitaevskii (GP) equation whose solution is approximated using a Lagrangian Variational Method to facilitate rapid prototyping. The method presented here is an extension of an earlier one that was used to analyze the results of an experiment [J.E. Simsarian, et al., Phys. Rev. Lett. 83, 2040 (2000)], where the phase of a Bose-Einstein condensate was measured using a Mach- Zehnder-type Bragg AI. We have developed both 1D and 3D versions of this method and we have determined their validity by comparing their predicted interference patterns with those obtained by numerical integration of the 1D GP equation and with the results of the above experiment. We find excellent agreement between the 1D interference patterns predicted by this method and those found by the GP equation. We show that we can reproduce all of the results of that experiment without recourse to an ad hoc velocity-kick correction needed by the earlier method, including some experimental results that the earlier model did not predict. We also found that this method provides estimates of 1D interference patterns at least four orders-of-magnitude faster than direct numerical solution of the 1D GP equation.

preprint2010arXiv

Exploring complex phenomena using ultracold atoms in bichromatic lattices

With an underlying common theme of competing length scales, we study the many-body Schrödinger equation in a quasiperiodic potential and discuss its connection with the Kolmogorov-Arnold-Moser (KAM) problem of classical mechanics. We propose a possible visualization of such connection in experimentally accessible many-body observables. Those observables are useful probes for the three characteristic phases of the problem: the metallic, Anderson and band insulator phases. In addition, they exhibit fingerprints of non-linear phenomena such as Arnold tongues, bifurcations and devil's staircases. Our numerical treatment is complemented with a perturbative analysis which provides insight on the underlying physics. The perturbation theory approach is particularly useful in illuminating the distinction between the Anderson insulator and the band insulator phases in terms of paired sets of dimerized states.

preprint2010arXiv

Far-ultraviolet signatures of the 3He(n,tp) reaction in noble gas mixtures

Previous work showed that the 3He(n,tp) reaction in a cell of 3He at atmospheric pressure generated tens of far-ultraviolet photons per reacted neutron. Here we report amplification of that signal by factors of 1000 and more when noble gases are added to the cell. Calibrated filter-detector measurements show that this large signal is due to noble-gas excimer emissions, and that the nuclear reaction energy is converted to far-ultraviolet radiation with efficiencies of up to 30%. The results have been placed on an absolute scale through calibrations at the NIST SURF III synchrotron. They suggest possibilities for high-efficiency neutron detectors as an alternative to existing proportional counters.

preprint2010arXiv

Momentum-space engineering of gaseous Bose-Einstein condensates

We show how the momentum distribution of gaseous Bose--Einstein condensates can be shaped by applying a sequence of standing-wave laser pulses. We present a theory, whose validity for was demonstrated in an earlier experiment [L.\ Deng, et al., \prl {\bf 83}, 5407 (1999)], of the effect of a two-pulse sequence on the condensate wavefunction in momentum space. We generalize the previous result to the case of $N$ pulses of arbitrary intensity separated by arbitrary intervals and show how these parameters can be engineered to produce a desired final momentum distribution. We find that several momentum distributions, important in atom-interferometry applications, can be engineered with high fidelity with two or three pulses.

preprint2010arXiv

Noise correlations of one-dimensional Bose mixtures in optical lattices

We study the noise correlations of one-dimensional binary Bose mixtures, as a probe of their quantum phases. In previous work, we found a rich structure of many-body phases in such mixtures, such as paired and counterflow superfluidity. Here we investigate the signature of these phases in the noise correlations of the atomic cloud after time-of-flight expansion, using both Luttinger liquid theory and the time-evolving block decimation (TEBD) method. We find that paired and counterflow superfluidity exhibit distinctive features in the noise spectra. We treat both extended and inhomogeneous systems, and our numerical work shows that the essential physics of the extended systems is present in the trapped-atom systems of current experimental interest. For paired and counterflow superfluid phases, we suggest methods for extracting Luttinger parameters from noise correlation spectroscopy.

preprint2010arXiv

Phase fluctuations in anisotropic Bose condensates: from cigars to rings

We study the phase-fluctuating condensate regime of ultra-cold atoms trapped in a ring-shaped trap geometry, which has been realized in recent experiments. We first consider a simplified box geometry, in which we identify the conditions to create a state that is dominated by thermal phase-fluctuations, and then explore the experimental ring geometry. In both cases we demonstrate that the requirement for strong phase fluctuations can be expressed in terms of the total number of atoms and the geometric length scales of the trap only. For the ring-shaped trap we discuss the zero temperature limit in which a condensate is realized where the phase is fluctuating due to interactions and quantum fluctuations. We also address possible ways of detecting the phase fluctuating regime in ring condensates.

preprint2010arXiv

State-insensitive bichromatic optical trapping

We propose a scheme for state-insensitive trapping of neutral atoms by using light with two independent wavelengths. In particular, we describe the use of trapping and control lasers to minimize the variance of the potential experienced by a trapped Rb atom in ground and excited states. We present calculated values of wavelength pairs for which the 5s and 5p_{3/2} levels have the same ac Stark shifts in the presence of two laser fields.

preprint2010arXiv

Theory and applications of atomic and ionic polarizabilities

Atomic polarization phenomena impinge upon a number of areas and processes in physics. The dielectric constant and refractive index of any gas are examples of macroscopic properties that are largely determined by the dipole polarizability. When it comes to microscopic phenomena, the existence of alkaline-earth anions and the recently discovered ability of positrons to bind to many atoms are predominantly due to the polarization interaction. An imperfect knowledge of atomic polarizabilities is presently looming as the largest source of uncertainty in the new generation of optical frequency standards. Accurate polarizabilities for the group I and II atoms and ions of the periodic table have recently become available by a variety of techniques. These include refined many-body perturbation theory and coupled-cluster calculations sometimes combined with precise experimental data for selected transitions, microwave spectroscopy of Rydberg atoms and ions, refractive index measurements in microwave cavities, ab initio calculations of atomic structures using explicitly correlated wave functions, interferometry with atom beams, and velocity changes of laser cooled atoms induced by an electric field. This review examines existing theoretical methods of determining atomic and ionic polarizabilities, and discusses their relevance to various applications with particular emphasis on cold-atom physics and the metrology of atomic frequency standards.

preprint2009arXiv

Blackbody-radiation shift in a 88Sr+ ion optical frequency standard

The blackbody radiation (BBR) shift of the 5s - 4d_{5/2} clock transition in 88Sr+ is calculated to be 0.250(9) Hz at room temperature, T=300K, using the relativistic all-order method where all single and double excitations of the Dirac-Fock wave function are included to all orders of perturbation theory. The BBR shift is a major component in the uncertainty budget of the optical frequency standard based on the 88Sr+ trapped ion. The scalar polarizabilities of the 5s and 4d_{5/2} levels, as well as the tensor polarizability of the 4d_{5/2} level, are presented together with the evaluation of their uncertainties. The lifetimes of the 4d_{3/2}, 4d_{5/2}, 5p_{1/2}, and 5p_{3/2} states are calculated and compared with experimental values.

preprint2009arXiv

Collisional cooling of ultra-cold atom ensembles using Feshbach resonances

We propose a new type of cooling mechanism for ultra-cold fermionic atom ensembles, which capitalizes on the energy dependence of inelastic collisions in the presence of a Feshbach resonance. We first discuss the case of a single magnetic resonance, and find that the final temperature and the cooling rate is limited by the width of the resonance. A concrete example, based on a p-wave resonance of $^{40}$K, is given. We then improve upon this setup by using both a very sharp optical or radio-frequency induced resonance and a very broad magnetic resonance and show that one can improve upon temperatures reached with current technologies.

preprint2009arXiv

Counterflow and paired superfluidity in one-dimensional Bose mixtures in optical lattices

We study the quantum phases of mixtures of ultra-cold bosonic atoms held in an optical lattice that confines motion or hopping to one spatial dimension. The phases are found by using Tomonaga-Luttinger liquid theory as well as the numerical method of time evolving block decimation (TEBD). We consider a binary mixture with repulsive intra-species interactions, and either repulsive or attractive inter-species interaction. For a homogeneous system, we find paired- and counterflow-superfluid phases at different filling and hopping energies. We also predict parameter regions in which these types of superfluid order coexist with charge density wave order. We show that the Tomonaga-Luttinger liquid theory and TEBD qualitatively agree on the location of the phase boundary to superfluidity. We then describe how these phases are modified and can be detected when an additional harmonic trap is present. In particular, we show how experimentally measurable quantities, such as time-of-flight images and the structure factor, can be used to distinguish the quantum phases. Finally, we suggest applying a Feshbach ramp to detect the paired superfluid state, and a $π/2$ pulse followed by Bragg spectroscopy to detect the counterflow superfluid phase.

preprint2009arXiv

Mesoscopic Effects in Quantum Phases of Ultracold Quantum Gases in Optical Lattices

We present a wide array of quantum measures on numerical solutions of 1D Bose- and Fermi-Hubbard Hamiltonians for finite-size systems with open boundary conditions. Finite size effects are highly relevant to ultracold quantum gases in optical lattices, where an external trap creates smaller effective regions in the form of the celebrated "wedding cake" structure and the local density approximation is often not applicable. Specifically, for the Bose-Hubbard Hamiltonian we calculate number, quantum depletion, local von-Neumann entropy, generalized entanglement or Q-measure, fidelity, and fidelity susceptibility; for the Fermi-Hubbard Hamiltonian we also calculate the pairing correlations, magnetization, charge-density correlations, and antiferromagnetic structure factor. Our numerical method is imaginary time propagation via time-evolving block decimation. As part of our study we provide a careful comparison of canonical vs. grand canonical ensembles and Gutzwiller vs. entangled simulations. The most striking effect of finite size occurs for bosons: we observe a strong blurring of the tips of the Mott lobes accompanied by higher depletion, and show how the location of the first Mott lobe tip approaches the thermodynamic value as a function of system size.

preprint2009arXiv

Particle-Hole Asymmetry and Brightening of Solitons in A Strongly Repulsive BEC

We study solitary wave propagation in the condensate of a system of hard-core bosons with nearest-neighbor interactions. For this strongly repulsive system, the evolution equation for the condensate order parameter of the system, obtained using spin coherent state averages is different from the usual Gross-Pitaevskii equation (GPE). The system is found to support two kinds of solitons when there is a particle-hole imbalance: a dark soliton that dies out as the velocity approaches the sound velocity, and a new type of soliton which brightens and persists all the way up to the sound velocity, transforming into a periodic wave train at supersonic speed. Analogous to the GPE soliton, the energy-momentum dispersion for both solitons is characterized by Lieb II modes.

preprint2009arXiv

Quantum Many-Body Dynamics of Dark Solitons in Optical Lattices

We present a fully quantum many-body treatment of dark solitons formed by ultracold bosonic atoms in one-dimensional optical lattices. Using time-evolving block decimation to simulate the single-band Bose-Hubbard Hamiltonian, we consider the quantum dynamics of density and phase engineered dark solitons as well as the quantum evolution of mean-field dark solitons injected into the quantum model. The former approach directly models how one may create quantum entangled dark solitons in experiment. While we have already presented results regarding the latter approach elsewhere [Phys. Rev. Lett. {\bf 103}, 140403 (2009)], we expand upon those results in this work. In both cases, quantum fluctuations cause the dark soliton to fill in and may induce an inelasticity in soliton-soliton collisions. Comparisons are made to the Bogoliubov theory which predicts depletion into an anomalous mode that fills in the soliton. Our many-body treatment allows us to go beyond the Bogoliubov approximation and calculate explicitly the dynamics of the system's natural orbitals.

preprint2008arXiv

Phase transitions, entanglement and quantum noise interferometry in cold atoms

We show that entanglement monotones can characterize the pronounced enhancement of entanglement at a quantum phase transition if they are sensitive to long-range high order correlations. These monotones are found to develop a sharp peak at the critical point and to exhibit universal scaling. We demonstrate that similar features are shared by noise correlations and verify that these experimentally accessible quantities indeed encode entanglement information and probe separability.

preprint2006arXiv

Lifshitz-like transition and enhancement of correlations in a rotating bosonic ring lattice

We study the effects of rotation on one-dimensional ultra-cold bosons confined to a ring lattice. For commensurate systems, at a critical value of the rotation frequency, an infinitesimal interatomic interaction energy opens a gap in the excitation spectrum, fragments the ground state into a macroscopic superposition of two states with different circulation and generates a sudden change in the topology of the momentum distribution. These features are reminiscent of the topological changes in the Fermi surface that occurs in the Lifshitz transition in fermionic systems. The entangled nature of the ground state induces a strong enhancement of quantum correlations and decreases the threshold for the Mott insulator transition. In contrast to the commensurate case, the incommensurate lattice is rather insensitive to rotation. Our studies demonstrate the utility of noise correlations as a tool for identifying new physics in strongly correlated systems.

preprint2006arXiv

Vortices in attractive Bose-Einstein condensates in two dimensions

The form and stability of quantum vortices in Bose-Einstein condensates with attractive atomic interactions is elucidated. They appear as ring bright solitons, and are a generalization of the Townes soliton to nonzero winding number $m$. An infinite sequence of radially excited stationary states appear for each value of $m$, which are characterized by concentric matter-wave rings separated by nodes, in contrast to repulsive condensates, where no such set of states exists. It is shown that robustly stable as well as unstable regimes may be achieved in confined geometries, thereby suggesting that vortices and their radial excited states can be observed in experiments on attractive condensates in two dimensions.

preprint2004arXiv

Quantum Kinetic Theory of BEC Lattice Gas:Boltzmann Equations from 2PI-CTP Effective Action

We continue our earlier work [Ana Maria Rey, B. L. Hu, Esteban Calzetta, Albert Roura and Charles W. Clark, Phys. Rev. A 69, 033610 (2004)] on the nonequilibrium dynamics of a Bose Einstein condensate (BEC) selectively loaded into every third site of a one-dimensional optical lattice. From the two-particle irreducible (2PI) closed-time-path (CTP) effective action for the Bose- Hubbard Hamiltonian, we show how to obtain the Kadanoff-Baym equations of quantum kinetic theory. Using the quasiparticle approximation, we show that the local equilibrium solutions of these equations reproduce the second- order corrections to the self-energy originally derived by Beliaev. This work paves the way for the use of effective action methods in the derivation of quantum kinetic theory of many atom systems.

preprint2004arXiv

Scalable quantum computation in systems with Bose-Hubbard dynamics

Several proposals for quantum computation utilize a lattice type architecture with qubits trapped by a periodic potential. For systems undergoing many body interactions described by the Bose-Hubbard Hamiltonian, the ground state of the system carries number fluctuations that scale with the number of qubits. This process degrades the initialization of the quantum computer register and can introduce errors during error correction. In an earlier manuscript we proposed a solution to this problem tailored to the loading of cold atoms into an optical lattice via the Mott Insulator phase transition. It was shown that by adding an inhomogeneity to the lattice and performing a continuous measurement, the unit filled state suitable for a quantum computer register can be maintained. Here, we give a more rigorous derivation of the register fidelity in homogeneous and inhomogeneous lattices and provide evidence that the protocol is effective in the finite temperature regime.

preprint2002arXiv

Dynamics of a period-three pattern loaded Bose-Einstein condensate in an optical lattice

We discuss the dynamics of a Bose-Einstein condensate initially loaded into every third site of an optical lattice using a description based upon the discrete nonlinear Schrodinger equation. An analytic solution is developed for the case of a periodic initial condition and is compared with numerical simulations for more general initial configurations. We show that mean field effects in this system can cause macroscopic quantum self-trapping, a phenomenon already predicted for double well systems. In the presence of a uniform external potential, the atoms exhibit generalized Bloch oscillations which can be interpreted in terms of the interference of three different Bloch states. We also discuss how the momentum distribution of the system can be used as experimental signature of the macroscopic self trapping effect.