Source author record

Kaden R. A. Hazzard

Kaden R. A. Hazzard 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

35works
10topics
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

35 published item(s)

preprint2026arXiv

Quantum science with arrays of metastable helium-3 atoms

The motion of atoms in programmable optical tweezer arrays offers many new opportunities for neutral atom quantum science. These include inter- and intra-site atom motion for resource-efficient implementations of fermionic and bosonic modes, respectively, as well as tweezer transport for efficient compilation of arbitrary circuits. However, the exploitation of atomic motion for all three purposes and others is limited by the inertia of the atoms. We present a comprehensive architectural blueprint for the use of fermionic metastable helium-3 ($^3$He$^*$) atoms -- the lightest trappable atomic species -- in programmable optical tweezer arrays. This includes a concrete analysis of atomic structure considerations as well as Rydberg-mediated interactions. We show that inter-tweezer hopping of $^3$He$^*$ atoms can be $\gtrsim3\times$ faster than previous demonstrations with lithium-6. We also demonstrate a new toolbox for encoding and manipulating qubits directly in the tweezer trap potential, uniquely enabled by the light mass of $^3$He$^*$. Finally, we provide several examples of new opportunities for fermionic quantum simulation and computation that leverage the transport and inter-tweezer hopping of $^3$He$^*$ atom arrays. These tools present new methods to improve the resource efficiency of neutral atom quantum science that may also enable quantum simulations of lattice gauge theories and quantum chemistry outside the Born-Oppenheimer approximation

preprint2022arXiv

A tensor network discriminator architecture for classification of quantum data on quantum computers

We demonstrate the use of matrix product state (MPS) models for discriminating quantum data on quantum computers using holographic algorithms, focusing on classifying a translationally invariant quantum state based on $L$ qubits of quantum data extracted from it. We detail a process in which data from single-shot experimental measurements are used to optimize an isometric tensor network, the tensors are compiled into unitary quantum operations using greedy compilation heuristics, parameter optimization on the resulting quantum circuit model removes the post-selection requirements of the isometric tensor model, and the resulting quantum model is inferenced on either product state or entangled quantum data. We demonstrate our training and inference architecture on a synthetic dataset of six-site single-shot measurements from the bulk of a one-dimensional transverse field Ising model (TFIM) deep in its antiferromagnetic and paramagnetic phases. We experimentally evaluate models on Quantinuum's H1-2 trapped ion quantum computer using entangled input data modeled as translationally invariant, bond dimension 4 MPSs across the known quantum phase transition of the TFIM. Using linear regression on the experimental data near the transition point, we find predictions for the critical transverse field of $h=0.962$ and $0.994$ for tensor network discriminators of bond dimension $χ=2$ and $χ=4$, respectively. These predictions compare favorably with the known transition location of $h=1$ despite training on data far from the transition point. Our techniques identify families of short-depth variational quantum circuits in a data-driven and hardware-aware fashion and robust classical techniques to precondition the model parameters, and can be adapted beyond machine learning to myriad applications of tensor networks on quantum computers, such as quantum simulation and error correction.

preprint2022arXiv

Quantum Membrane Phases in Synthetic Lattices of Cold Molecules or Rydberg Atoms

We calculate properties of dipolar interacting ultracold molecules or Rydberg atoms in a semi-synthetic three-dimensional configuration -- one synthetic dimension plus a two-dimensional real space optical lattice or periodic microtrap array -- using the stochastic Green function Quantum Monte Carlo method. Through a calculation of thermodynamic quantities and appropriate correlation functions, along with their finite size scalings, we show that there is a second order transition to a low temperature phase in which two-dimensional `sheets' form in the synthetic dimension of internal rotational or electronic states of the molecules or Rydberg atoms, respectively. Simulations for different values of the interaction $V$, which acts between atoms or molecules that are adjacent both in real and synthetic space, allow us to compute a phase diagram. We find a finite-temperature transition at sufficiently large $V$, as well as a quantum phase transition -- a critical value $V_c$ below which the transition temperature vanishes.

preprint2022arXiv

Topological correlations in three dimensional classical Ising models: an exact solution with a continuous phase transition

We study a three-dimensional (3D) classical Ising model that is exactly solvable when some coupling constants take certain imaginary values. The solution combines and generalizes the Onsager-Kaufman solution of the 2D Ising model and the solution of Kitaev's honeycomb model, leading to a three-parameter phase diagram with a third order phase transition between two distinct phases. Interestingly, the phases of this model are distinguished by topological features: the expectation value of a certain family of loop observables depend only on the topology of the loop (whether the loop is contractible), and are quantized at rational values that differ in the two phases. We show that a related exactly solvable 3D classical statistical model with real coupling constants also shows the topological features of one of these phases. Furthermore, even in the model with complex parameters, the partition function has some physical relevance, as it can be interpreted as the transition amplitude of a quantum dynamical process and may shed light on dynamical quantum phase transitions.

preprint2021arXiv

Universal thermodynamics of an SU($N$) Fermi-Hubbard Model

The SU(2) symmetric Fermi-Hubbard model (FHM) plays an essential role in strongly correlated fermionic many-body systems. In the one particle per site and strongly interacting limit ${U/t \gg 1}$, it is effectively described by the Heisenberg Hamiltonian. In this limit, enlarging the spin and extending the typical SU(2) symmetry to SU($N$) has been predicted to give exotic phases of matter in the ground state, with a complicated dependence on $N$. This raises the question of what -- if any -- are the finite-temperature signatures of these phases, especially in the currently experimentally relevant regime near or above the superexchange energy. We explore this question for thermodynamic observables by numerically calculating the thermodynamics of the SU($N$) FHM in the two-dimensional square lattice near densities of one particle per site, using determinant Quantum Monte Carlo and Numerical Linked Cluster Expansion. Interestingly, we find that for temperatures above the superexchange energy, where the correlation length is short, the energy, number of on-site pairs, and kinetic energy are universal functions of $N$. Although the physics in the regime studied is well beyond what can be captured by low-order high-temperature series, we show that an analytic description of the scaling is possible in terms of only one- and two-site calculations.

preprint2020arXiv

Collective modes of ultracold fermionic alkaline-earth gases with SU(N) symmetry

We calculate the collective modes of ultracold trapped alkaline-earth fermionic atoms, which possess an SU($N$) symmetry of the nuclear spin degree of freedom, and a controllable $N$, with $N$ as large as $10$. We calculate the breathing and quadrupole modes of two-dimensional and three-dimensional harmonically trapped gases in the normal phase. We particularly concentrate on two-dimensional gases, where the shift is more accessible experimentally, and the physics has special features. We present results as a function of temperature, interaction strength, density, and $N$. We include calculations across the collisionless to hydrodynamic crossover. We assume the gas is interacting weakly, such that it can be described by a Boltzmann-Vlasov equation that includes both mean-field terms and the collision integral. We solve this with an approximate scaling ansatz, taking care in two-dimensions to preserve the scaling symmetry of the system. We predict the collective mode frequency shifts and damping, showing that these are measurable in experimentally relevant regimes. We expect these results to furnish powerful tools to characterize interactions and the state of alkaline-earth gases, as well as to lay the foundation for future work, for example on strongly interacting gases and SU($N$) spin modes.

preprint2020arXiv

Numerical linked cluster expansions for inhomogeneous systems

We develop a numerical linked cluster expansion (NLCE) method that can be applied directly to inhomogeneous systems, for example Hamiltonians with disorder and dynamics initiated from inhomogeneous initial states. We demonstrate the method by calculating dynamics for single-spin expectations and spin correlations in two-dimensional spin models on a square lattice, starting from a checkerboard state. We show that NLCE can give moderate to dramatic improvement over an exact diagonalization of comparable computational cost, and that the advantage in computational resources grows exponentially as the size of the clusters included grows. Although the method applies to any type of NLCE, our explicit benchmarks use the rectangle expansion. Besides showing the capability to treat inhomogeneous systems, these benchmarks demonstrate the rectangle expansion's utility out of equilibrium.

preprint2020arXiv

Observation of antiferromagnetic correlations in an ultracold SU($N$) Hubbard model

Mott insulators are paradigms of strongly correlated physics, giving rise to phases of matter with novel and hard-to-explain properties. Extending the typical SU(2) symmetry of Mott insulators to SU($N$) is predicted to give exotic quantum magnetism at low temperatures, but understanding the effect of strong quantum fluctuations for large $N$ remains an open challenge. In this work, we experimentally observe nearest-neighbor spin correlations in the SU(6) Hubbard model realized by ytterbium atoms in optical lattices. We study one-dimensional, two-dimensional square, and three-dimensional cubic lattice geometries. The measured SU(6) spin correlations are dramatically enhanced compared to the SU(2) correlations, due to strong Pomeranchuk cooling. We also present numerical calculations based on exact diagonalization and determinantal quantum Monte Carlo. The experimental data for a one-dimensional lattice agree with theory, without any fitting parameters. The detailed comparison between theory and experiment allows us to infer from the measured correlations a lowest temperature of $\left[{0.096 \pm 0.054 \, \rm{(theory)} \pm 0.030 \, \rm{(experiment)}}\right]/k_{\rm B}$ times the tunneling amplitude. For two- and three-dimensional lattices, experiments reach entropies below where our calculations converge, highlighting the experiments as quantum simulations. These results open the door for the study of long-sought SU($N$) quantum magnetism.

preprint2020arXiv

Spin-imbalanced ultracold Fermi gases in a two-dimensional array of tubes

Motivated by a recent experiment [Revelle et al. Phys. Rev. Lett. 117, 235301 (2016)] that characterized the one- to three-dimensional crossover in a spin-imbalanced ultracold gas of $^6$Li atoms trapped in a two-dimensional array of tunnel-coupled tubes, we calculate the phase diagram for this system using Hartree-Fock Bogoliubov-de Gennes mean-field theory, and compare the results with experimental data. Mean-field theory predicts fully spin-polarized normal, partially spin-polarized normal, spin-polarized superfluid, and spin-balanced superfluid phases in a homogeneous system. We use the local density approximation to obtain density profiles of the gas in a harmonic trap. We compare these calculations with experimental measurements in Revelle {\em et al.} as well as previously unpublished data. Our calculations qualitatively agree with experimentally-measured densities and coordinates of the phase boundaries in the trap, and quantitatively agree with experimental measurements at moderate-to-large polarizations. Our calculations also reproduce the experimentally-observed universal scaling of the phase boundaries for different scattering lengths at a fixed value of scaled inter-tube tunneling. However, our calculations have quantitative differences with experimental measurements at low polarization, and fail to capture important features of the one- to three-dimensional crossover observed in experiments. These suggest the important role of physics beyond-mean-field theory in the experiments. We expect that our numerical results will aid future experiments in narrowing the search for the FFLO phase.

preprint2020arXiv

Tightening the Lieb-Robinson Bound in Locally-Interacting Systems

The Lieb-Robinson (LR) bound rigorously shows that in quantum systems with short-range interactions, the maximum amount of information that travels beyond an effective "light cone" decays exponentially with distance from the light-cone front, which expands at finite velocity. Despite being a fundamental result, existing bounds are often extremely loose, limiting their applications. We introduce a method that dramatically and qualitatively improves LR bounds in models with finite-range interactions. Most prominently, in systems with a large local Hilbert space dimension $D$, our method gives an LR velocity that grows much slower than previous bounds with $D$ as $D\to \infty$. For example, in the Heisenberg model with spin $S$, we find $v\leq$ const. compared to the previous $v\propto S$ which diverges at large $S$, and in multiorbital Hubbard models with $N$ orbitals, we find $v\propto \sqrt{N}$ instead of previous $v\propto N$, and similarly in the $N$-state truncated Bose-Hubbard model and Wen's quantum rotor model. Our bounds also scale qualitatively better in some systems when the spatial dimension or certain model parameters become large, for example in the $d$-dimensional quantum Ising model and perturbed toric code models. Even in spin-1/2 Ising and Fermi-Hubbard models, our method improves the LR velocity by an order of magnitude with typical model parameters, and significantly improves the LR bound at large distance and early time.

preprint2019arXiv

Quantum Simulators: Architectures and Opportunities

Quantum simulators are a promising technology on the spectrum of quantum devices from specialized quantum experiments to universal quantum computers. These quantum devices utilize entanglement and many-particle behaviors to explore and solve hard scientific, engineering, and computational problems. Rapid development over the last two decades has produced more than 300 quantum simulators in operation worldwide using a wide variety of experimental platforms. Recent advances in several physical architectures promise a golden age of quantum simulators ranging from highly optimized special purpose simulators to flexible programmable devices. These developments have enabled a convergence of ideas drawn from fundamental physics, computer science, and device engineering. They have strong potential to address problems of societal importance, ranging from understanding vital chemical processes, to enabling the design of new materials with enhanced performance, to solving complex computational problems. It is the position of the community, as represented by participants of the NSF workshop on "Programmable Quantum Simulators," that investment in a national quantum simulator program is a high priority in order to accelerate the progress in this field and to result in the first practical applications of quantum machines. Such a program should address two areas of emphasis: (1) support for creating quantum simulator prototypes usable by the broader scientific community, complementary to the present universal quantum computer effort in industry; and (2) support for fundamental research carried out by a blend of multi-investigator, multi-disciplinary collaborations with resources for quantum simulator software, hardware, and education.

preprint2016arXiv

Accessing Rydberg-dressed interactions using many-body Ramsey dynamics

We demonstrate that Ramsey spectroscopy can be used to observe Rydberg-dressed interactions. In contrast to many prior proposals, our scheme operates comfortably within experimentally measured lifetimes, and accesses a regime where quantum superpositions are crucial. The key idea is to build a spin-1/2 from one level that is Rydberg-dressed and another that is not. These levels may be hyperfine or long-lived electronic states. An Ising spin model governs the Ramsey dynamics, for which we derive an exact solution. Due to the structure of Rydberg interactions, the dynamics differs significantly from that in other spin systems. As one example, spin echo can increase the rate at which coherence decays. The results also apply to bare (undressed) Rydberg states as a special case, for which we quantitatively reproduce recent ultrafast experiments without fitting.

preprint2015arXiv

Effective many-body parameters for atoms in non-separable Gaussian optical potentials

We analyze the properties of particles trapped in three-dimensional potentials formed from superimposed Gaussian beams, fully taking into account effects of potential anharmonicity and non-separability. Although these effects are negligible in more conventional optical lattice experiments, they are essential for emerging ultracold atom developments. We focus in particular on two potentials utilized in current ultracold atom experiments: arrays of tightly focused optical tweezers and a one-dimensional optical lattice with transverse Gaussian confinement and highly excited transverse modes. Our main numerical tools are discrete variable representations (DVRs), which combine many favorable features of spectral and grid-based methods, such as the computational advantage of exponential convergence and the convenience of an analytical representation of Hamiltonian matrix elements. Optimizations, such as symmetry adaptations and variational methods built on top of DVR methods, are presented and their convergence properties discussed. We also present a quantitative analysis of the degree of non-separability of eigenstates, borrowing ideas from the theory of matrix product states (MPSs), leading to both conceptual and computational gains. Beyond developing numerical methodologies, we present results for construction of optimally localized Wannier functions and tunneling and interaction matrix elements in optical lattices and tweezers relevant for constructing effective models for many-body physics.

preprint2015arXiv

Synthetic gauge fields stabilize a chiral spin liquid phase

We calculate the phase diagram of the SU($N$) Hubbard model describing fermionic alkaline earth atoms in a square optical lattice with on-average one atom per site, using a slave-rotor mean-field approximation. We find that the chiral spin liquid predicted for $N\ge5$ and large interactions passes through a fractionalized state with a spinon Fermi surface as interactions are decreased before transitioning to a weakly interacting metal. We also show that by adding an artificial uniform magnetic field with flux per plaquette $2π/N$, the chiral spin liquid becomes the ground state for all $N\ge 3$ at large interactions, persists to weaker interactions, and its spin gap increases, suggesting that the spin liquid physics will persist to higher temperatures. We discuss potential methods to realize the artificial gauge fields and detect the predicted phases.

preprint2015arXiv

Ultracold nonreactive molecules in an optical lattice: connecting chemistry to many-body physics

We derive effective lattice models for ultracold bosonic or fermionic nonreactive molecules (NRMs) in an optical lattice, analogous to the Hubbard model that describes ultracold atoms in a lattice. In stark contrast to the Hubbard model, which is commonly assumed to accurately describe NRMs, we find that the single on-site interaction parameter $U$ is replaced by a multi-channel interaction, whose properties we elucidate. The complex, multi-channel collisional physics is unrelated to dipolar interactions, and so occurs even in the absence of an electric field or for homonuclear molecules. We find a crossover between coherent few-channel models and fully incoherent single-channel models as the lattice depth is increased. We show that the effective model parameters can be determined in lattice modulation experiments, which consequently measure molecular collision dynamics with a vastly sharper energy resolution than experiments in an ultracold gas.

preprint2014arXiv

Many-body dynamics of dipolar molecules in an optical lattice

Understanding the many-body dynamics of isolated quantum systems is one of the central challenges in modern physics. To this end, the direct experimental realization of strongly correlated quantum systems allows one to gain insights into the emergence of complex phenomena. Such insights enable the development of theoretical tools that broaden our understanding. Here, we theoretically model and experimentally probe with Ramsey spectroscopy the quantum dynamics of disordered, dipolar-interacting, ultracold molecules in a partially filled optical lattice. We report the capability to control the dipolar interaction strength, and we demonstrate that the many-body dynamics extends well beyond a nearest-neighbor or mean-field picture, and cannot be quantitatively described using previously available theoretical tools. We develop a novel cluster expansion technique and demonstrate that our theoretical method accurately captures the measured dependence of the spin dynamics on molecule number and on the dipolar interaction strength. In the spirit of quantum simulation, this agreement simultaneously benchmarks the new theoretical method and verifies our microscopic understanding of the experiment. Our findings pave the way for numerous applications in quantum information science, metrology, and condensed matter physics.

preprint2014arXiv

Quantum correlations and entanglement in far-from-equilibrium spin systems

By applying complementary analytic and numerical methods, we investigate the dynamics of spin-$1/2$ XXZ models with variable-range interactions in arbitrary dimensions. The dynamics we consider is initiated from uncorrelated states that are easily prepared in experiments, and can be equivalently viewed as either Ramsey spectroscopy or a quantum quench. Our primary focus is the dynamical emergence of correlations and entanglement in these far-from-equilibrium interacting quantum systems: we characterize these correlations by the entanglement entropy, concurrence, and squeezing, which are inequivalent measures of entanglement corresponding to different quantum resources. In one spatial dimension, we show that the time evolution of correlation functions manifests a non-perturbative dynamic singularity. This singularity is characterized by a universal power-law exponent that is insensitive to small perturbations. Explicit realizations of these models in current experiments using polar molecules, trapped ions, Rydberg atoms, magnetic atoms, and alkaline-earth and alkali atoms in optical lattices, along with the relative merits and limitations of these different systems, are discussed.

preprint2013arXiv

Kitaev honeycomb and other exotic spin models with polar molecules

We show that ultracold polar molecules pinned in an optical lattice can be used to access a variety of exotic spin models, including the Kitaev honeycomb model. Treating each molecule as a rigid rotor, we use DC electric and microwave fields to define superpositions of rotational levels as effective spin degrees of freedom, while dipole-dipole interactions give rise to interactions between the spins. In particular, we show that, with sufficient microwave control, the interaction between two spins can be written as a sum of five independently controllable Hamiltonian terms proportional to the five rank-2 spherical harmonics Y_{2,q}(theta,phi), where (theta,phi) are the spherical coordinates of the vector connecting the two molecules. To demonstrate the potential of this approach beyond the simplest examples studied in [S. R. Manmana et al., arXiv:1210.5518v2], we focus on the realization of the Kitaev honeycomb model, which can support exotic non-Abelian anyonic excitations. We also discuss the possibility of generating spin Hamiltonians with arbitrary spin S, including those exhibiting SU(N=2S+1) symmetry.

preprint2013arXiv

Realizing a lattice spin model with polar molecules

With the recent production of polar molecules in the quantum regime, long-range dipolar interactions are expected to facilitate the understanding of strongly interacting many-body quantum systems and to realize lattice spin models for exploring quantum magnetism. In atomic systems, where interactions require wave function overlap, effective spin interactions on a lattice can be realized via superexchange; however, the coupling is weak and limited to nearest-neighbor interactions. In contrast, dipolar interactions exist in the absence of tunneling and extend beyond nearest neighbors. This allows coherent spin dynamics to persist even at high entropy and low lattice filling. Effects of dipolar interactions in ultracold molecular gases have so far been limited to the modification of chemical reactions. We now report the observation of dipolar interactions of polar molecules pinned in a 3D optical lattice. We realize a lattice spin model with spin encoded in rotational states, prepared and probed by microwaves. This spin-exchange interaction arises from the resonant exchange of rotational angular momentum between two molecules. We observe clear oscillations in the evolution of the spin coherence in addition to an overall decay. The frequency of these oscillations, the strong dependence of the spin coherence time on the lattice filling, and the effect of a multi-pulse sequence designed to reverse dynamics due to two-body exchange interactions all provide evidence of dipolar interactions. We also demonstrate suppression of loss in weak lattices due to a quantum Zeno mechanism. Measurements of these tunneling-induced losses allow us to independently determine the lattice filling factor. These results comprise an initial exploration of the behavior of many-body spin models with direct, long-range spin interactions and lay the groundwork for future studies of many-body dynamics in spin lattices.

preprint2013arXiv

Suppressing the loss of ultracold molecules via the continuous quantum Zeno effect

We investigate theoretically the suppression of two-body losses when the on-site loss rate is larger than all other energy scales in a lattice. This work quantitatively explains the recently observed suppression of chemical reactions between two rotational states of fermionic KRb molecules confined in one-dimensional tubes with a weak lattice along the tubes [Yan et al., Nature 501, 521-525 (2013)]. New loss rate measurements performed for different lattice parameters but under controlled initial conditions allow us to show that the loss suppression is a consequence of the combined effects of lattice confinement and the continuous quantum Zeno effect. A key finding, relevant for generic strongly reactive systems, is that while a single-band theory can qualitatively describe the data, a quantitative analysis must include multiband effects. Accounting for these effects reduces the inferred molecule filling fraction by a factor of five. A rate equation can describe much of the data, but to properly reproduce the loss dynamics with a fixed filling fraction for all lattice parameters we develop a mean-field model and benchmark it with numerically exact time-dependent density matrix renormalization group calculations.

preprint2012arXiv

Adiabatic loading of one-dimensional SU(N) alkaline earth fermions in optical lattices

Ultracold fermionic alkaline earth atoms confined in optical lattices realize Hubbard models with internal SU(N) symmetries, where N can be as large as ten. Such systems are expected to harbor exotic magnetic physics at temperatures below the superexchange energy scale. Employing quantum Monte Carlo simulations to access the low-temperature regime, we show that after adiabatically loading a weakly interacting gas into the strongly interacting regime of an optical lattice, the final temperature decreases with increasing N. Furthermore, we estimate the temperature scale required to probe correlations associated with low-temperature SU(N) magnetism. Our findings are encouraging for the exploration of exotic large-N magnetic states in ongoing experiments.

preprint2012arXiv

Far from equilibrium quantum magnetism with ultracold polar molecules

Recent theory has indicated how to emulate tunable models of quantum magnetism with ultracold polar molecules. Here we show that present molecule optical lattice experiments can accomplish three crucial goals for quantum emulation, despite currently being well below unit filling and not quantum degenerate. The first is to verify and benchmark the models proposed to describe these systems. The second is to prepare correlated and possibly useful states in well-understood regimes. The third is to explore many-body physics inaccessible to existing theoretical techniques. Our proposal relies on a non-equilibrium protocol that can be viewed either as Ramsey spectroscopy or an interaction quench. It uses only routine experimental tools available in any ultracold molecule experiment.

preprint2012arXiv

Non-equilibrium dynamics of Ising models with decoherence: an exact solution

The interplay between interactions and decoherence in many-body systems is of fundamental importance in quantum physics: Decoherence can degrade correlations, but can also give rise to a variety of rich dynamical and steady-state behaviors. We obtain an exact analytic solution for the non-equilibrium dynamics of Ising models with arbitrary interactions and subject to the most general form of local Markovian decoherence. Our solution shows that decoherence affects the relaxation of observables more than predicted by single-particle considerations. It also reveals a dynamical phase transition, specifically a Hopf bifurcation, which is absent at the single-particle level. These calculations are applicable to ongoing quantum information and emulation efforts using a variety of atomic, molecular, optical, and solid-state systems.

preprint2012arXiv

Topological phases in ultracold polar-molecule quantum magnets

We show how to use polar molecules in an optical lattice to engineer quantum spin models with arbitrary spin S >= 1/2 and with interactions featuring a direction-dependent spin anisotropy. This is achieved by encoding the effective spin degrees of freedom in microwave-dressed rotational states of the molecules and by coupling the spins through dipolar interactions. We demonstrate how one of the experimentally most accessible anisotropies stabilizes symmetry protected topological phases in spin ladders. Using the numerically exact density matrix renormalization group method, we find that these interacting phases -- previously studied only in the nearest-neighbor case -- survive in the presence of long-range dipolar interactions. We also show how to use our approach to realize the bilinear-biquadratic spin-1 and the Kitaev honeycomb models. Experimental detection schemes and imperfections are discussed.

preprint2011arXiv

Spectroscopy of dipolar fermions in 2D pancakes and 3D lattices

Motivated by ongoing measurements at JILA, we calculate the recoil-free spectra of dipolar interacting fermions, for example ultracold heteronuclear molecules, in a one-dimensional lattice of two-dimensional pancakes, spectroscopically probing transitions between different internal (e.g., rotational) states. We additionally incorporate p-wave interactions and losses, which are important for reactive molecules such as KRb. Moreover, we consider other sources of spectral broadening: interaction-induced quasiparticle lifetimes and the different polarizabilities of the different rotational states used for the spectroscopy. Although our main focus is molecules, some of the calculations are also useful for optical lattice atomic clocks. For example, understanding the p-wave shifts between identical fermions and small dipolar interactions coming from the excited clock state are necessary to reach future precision goals. Finally, we consider the spectra in a deep 3D lattice and show how they give a great deal of information about static correlation functions, including \textit{all} the moments of the density correlations between nearby sites. The range of correlations measurable depends on spectroscopic resolution and the dipole moment.

preprint2011arXiv

SU(N) magnetism in chains of ultracold alkaline-earth-metal atoms: Mott transitions and quantum correlations

We investigate one dimensional SU$(N)$ Hubbard chains at zero temperature, which can be emulated with ultracold alkaline earth atoms, by using the density matrix renormalization group (DMRG), Bethe ansatz (BA), and bosonization. We compute experimental observables and use the DMRG to benchmark the accuracy of the Bethe ansatz for $N>2$ where the BA is only approximate. In the worst case, we find a relative error $ε\lesssim 4%$ in the BA ground state energy for $N \leq 4$ at filling 1/N, which is due to the fact that BA improperly treats the triply and higher occupied states. Using the DMRG for $N \leq 4$ and the BA for large $N$, we determine the regimes of validity of strong- and weak-coupling perturbation theory for all values of $N$ and in particular, the parameter range in which the system is well described by a SU$(N)$ Heisenberg model at filling 1/N. We find this depends only weakly on $N$. We investigate the Berezinskii-Kosterlitz-Thouless phase transition from a Luttinger liquid to a Mott-insulator by computing the fidelity susceptibility and the Luttinger parameter $K_ρ$ at 1/N filling. The numerical findings give strong evidence that the fidelity susceptibility develops a minimum at a critical interaction strength which is found to occur at a finite positive value for $N>2$.

preprint2011arXiv

Universality class of quantum criticality in the two-dimensional Hubbard model at intermediate temperatures ($t^2/U\ll T\ll t$)

We show that the dilute Fermi gas quantum critical universality class quantitatively describes the Mott/metal crossover of the two-dimensional Hubbard model for temperatures somewhat less than (roughly half) the tunneling but much greater than (roughly twice) the superexchange energy. We calculate the observables expected to be universal near the transition --- density and compressibility --- with numerically exact determinantal quantum Monte Carlo. We find they are universal functions of the chemical potential. Despite arising from the strongly correlated regime of the Hubbard model, these functions are given by the weakly interacting, dilute Fermi gas model. These observables and their derivatives are the only expected universal static observables of this universality class, which we also confirm by verifying there is no scaling collapse of the kinetic energy, fraction of doubly occupied sites, and nearest neighbor spin correlations. Our work resolves the universality class of the intermediate temperature Mott/metal crossover, which had alternatively been proposed to be described by more exotic theories. However, in the presence of a Zeeman magnetic field, we find that interplay of spin with itinerant charge can lead to physics beyond the dilute Fermi gas universality class.

preprint2010arXiv

Candidate theories to explain the anomalous spectroscopic signatures of atomic H in molecular H$_2$ crystals

We analyze a number of proposed explanations for spectroscopic anomalies observed in atomic hydrogen defects embedded in a solid molecular hydrogen matrix. In particular, we critically evaluate the possibility that these anomalies are related to Bose-Einstein condensation (both global and local). For each proposed mechanism we discuss which aspects of the experiment can be explained, and make predictions for future experiments.

preprint2010arXiv

High temperature thermodynamics of fermionic alkaline earth atoms in optical lattices

We calculate experimentally relevant properties of trapped fermionic alkaline earth atoms in an optical lattice, modeled by the SU(N) Hubbard model. Our calculation is accurate when the temperature is much larger than the tunneling rate, similar to current regimes in ultracold atom experiments. In addition to exploring the Mott insulator-metal crossover, we calculate final temperatures achieved by the standard experimental protocol of adiabatically ramping from a non-interacting gas, as a function of initial gas temperature and final state lattice parameters. Of particular experimental interest, we find that increasing $N$ gives substantially \textit{colder} Mott insulators, up to more than a factor of five for relevant parameters. This cooling happens for all $N$, fixing the initial entropy, or for all $N \lsim 20$ (the exact value depends on dimensionality), fixing the initial temperature.

preprint2010arXiv

Local versus global equilibration near the bosonic Mott-superfluid transition

We study the response of trapped two dimensional cold bosons to time dependent lattices. We find that in lattice ramps from 11 (superfluid, $\hbar/U_{\text{i}} = 3$ms, $\hbar/J_{\text{i}} = 45$ms) to 16 recoils (Mott, $\hbar/U_{\text{f}} = 2$ms, $\hbar/J_{\text{f}} = 130$ms) the local number fluctuations remains at their equilibrium values if ramps are slower than 3 ms. Global transport, however, is much slower (1s), especially in the presence of Mott shells. This separation of timescales has practical implications for cold atom experiments and cooling protocols.

preprint2010arXiv

Techniques to measure quantum criticality in cold atoms

Attempts to understand zero temperature phase transitions have forced physicists to consider a regime where the standard paradigms of condensed matter physics break down [1-4]. These quantum critical systems lack a simple description in terms of weakly interacting quasiparticles, but over the past 20 years physicists have gained deep insights into their properties. Most dramatically, theory predicts that universal scaling relationships describe their finite temperature thermodynamics up to remarkably high temperatures. Unfortunately, these universal functions are hard to calculate: for example there are no reliable general techniques [4,5] to calculate the scaling functions for dynamics. Viewing a cold atom experiment as a quantum simulator [6], we show how to extract universal scaling functions from (non-universal) atomic density profiles or spectroscopic measurements. Such experiments can resolve important open questions about the Mott-Metal crossover [7,8] and the dynamics of the finite density O(2) rotor model [1,9], with direct impact on theories of, for example, high temperature superconducting cuprates [10,11], heavy fermion materials [12], and graphene [13].

preprint2009arXiv

Many-body physics in the radio frequency spectrum of lattice bosons

We calculate the radio-frequency spectrum of a trapped cloud of cold bosonic atoms in an optical lattice. Using random phase and local density approximations we produce both trap averaged and spatially resolved spectra, identifying simple features in the spectra that reveal information about both superfluidity and correlations. Our approach is exact in the deep Mott limit and in the deep superfluid when the hopping rates for the two internal spin states are equal. It contains final state interactions, obeys the Ward identities (and the associated conservation laws), and satisfies the $f$-sum rule. Motivated by earlier work by Sun, Lannert, and Vishveshwara [Phys. Rev. A \textbf{79}, 043422 (2009)], we also discuss the features which arise in a spin-dependent optical lattice.

preprint2009arXiv

On-site correlations in optical lattices: band mixing to coupled quantum Hall puddles

We extend the standard Bose-Hubbard model to capture arbitrarily strong on-site correlations. In addition to being important for quantitatively modeling experiments, for example, with Rubidium atoms, these correlations must be included to describe more exotic situations. Two such examples are when the interactions are made large via a Feshbach resonance, or when each site rotates rapidly, making a coupled array of quantum Hall puddles. Remarkably, even the mean field approximation to our model includes all on-site correlations. We describe how these on-site correlations manifest themselves in the system's global properties: modifying the phase diagram and depleting the condensate.

preprint2008arXiv

Explanation of 100-fold reduction of spectral shifts for hydrogen on helium films

We show that helium film-mediated hydrogen-hydrogen interactions account for a two orders of magnitude discrepancy between previous theory and recent experiments on cold collision shifts in spin-polarized hydrogen adsorbed on a helium film. These attractive interactions also explain the anomalous dependence of the cold collision frequency shifts on the $^3$He covering of the film. Our findings suggest that the gas will become mechanically unstable before reaching the Kosterlitz-Thouless transition unless the experiment is performed in a drastically different regime, for example with a much different helium film geometry.

preprint2007arXiv

Hyperfine spectra of trapped Bosons in optical lattices

We calculate the interaction induced inhomogeneous broadening of spectral lines in a trapped Bose gas as a function of the depth of a three-dimensional cubic optical lattice. As observed in recent experiments, we find that the terraced "wedding-cake" structure of Mott plateaus splits the spectrum into a series of discrete peaks. The spectra are extremely sensitive to density corrugations and trap anharmonicities. For example, even when the majority of the cloud is superfluid the spectrum displays discrete peaks.