Source author record

D. Blume

D. Blume 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

48works
5topics
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

48 published item(s)

preprint2024arXiv

Driven generalized quantum Rayleigh-van der Pol oscillators: Phase localization and spectral response

Driven classical self-sustained oscillators have been studied extensively in the context of synchronization. Using the master equation, this work considers the classically driven generalized quantum Rayleigh-van der Pol oscillator, which is characterized by linear dissipative gain and loss terms as well as three non-linear dissipative terms. Since two of the non-linear terms break the rotational phase space symmetry, the Wigner distribution of the quantum mechanical limit cycle state of the undriven system is, in general, not rotationally symmetric. The impact of the symmetry-breaking dissipators on the long-time dynamics of the driven system are analyzed as functions of the drive strength and detuning, covering the deep quantum to near-classical regimes. Phase localization and frequency entrainment, which are required for synchronization, are discussed in detail. We identify a large parameter space where the oscillators exhibit appreciable phase localization but only weak or no entrainment, indicating the absence of synchronization. Several observables are found to exhibit the analog of the celebrated classical Arnold tongue; in some cases, the Arnold tongue is found to be asymmetric with respect to vanishing detuning between the external drive and the natural oscillator frequency.

preprint2023arXiv

Finite SSH chains coupled to a two-level emitter: Hybridization of edge and emitter states

The Hamiltonian for the one-dimensional SSH chain is one of the simplest Hamiltonians that supports topological states. This work considers between one and three finite SSH chains with open boundary conditions that either share a lattice site (or cavity), which -- in turn -- is coupled to a two-level emitter, or are coupled to the same two-level emitter. We investigate the system properties as functions of the emitter-cavity coupling strength $g$ and the detuning between the emitter energy and the center of the band gap. It is found that the energy scale introduced by the edge states that are supported by the uncoupled finite SSH chains leads to a $g$-dependent hybridization of the emitter and edge states that is unique to finite-chain systems. A highly accurate analytical three-state model that captures the band gap physics of $k$-chain ($k \ge 1$) systems is developed. To quantify the robustness of the topological system characteristics, the inverse participation ratio for the cavity-shared and emitter-shared systems consisting of $k$ chains is analyzed as a function of the onsite disorder strength. The $g$-dependent hybridization of the emitter and uncoupled edge states can be probed dynamically.

preprint2023arXiv

Photon-induced droplet-like bound states in one-dimensional qubit array

We consider an array of $N_e$ non-interacting qubits or emitters that are coupled to a one-dimensional cavity array with tunneling energy $J$ and non-linearity of strength $U$. The number of cavities is assumed to be larger than the number of qubits. Working in the two-excitation manifold, we focus on the bandgap regime where the energy of two excited qubits is off-resonant with the two-photon bound state band. A two-step adiabatic elimination of the photonic degrees of freedom gives rise to a one-dimensional spin Hamiltonian with effective interactions; specifically, the Hamiltonian features constrained single-qubit hopping and pair hopping interactions not only between nearest neighbors but also between next-to-nearest and next-to-next-to-nearest spins. For a regularly arranged qubit array, we identify parameter combinations for which the system supports novel droplet-like bound states whose characteristics depend critically on the pair hopping. The droplet-like states can be probed dynamically. The bound states identified in our work for off-resonance conditions are distinct from localized hybridized states that emerge for on-resonance conditions.

preprint2022arXiv

Structural properties of bosonic He clusters with N=2-10 for different potential models at the physical point and at unitarity

Since the $^4$He dimer supports only one weakly bound state with an average interatomic distance much larger than the van der Waals length and no deeply bound states, $^4$He$_N$ clusters with $N>2$ are a paradigmatic model system with which to explore foundational concepts such as large $s$-wave scattering length universality, van der Waals universality, Efimov physics, and effective field theories. This work presents structural properties such as the pair and triple distribution functions, the hyperradial density, the probability to find the $N$th particle at a given distance from the center of mass of the other $N-1$ atoms, and selected contacts. The kinetic energy release, which can be measured via Coulomb explosion in dedicated size-selected molecular beam experiments -- at least for small $N$ -- , is also presented. The structural properties are determined for three different realistic $^4$He-$^4$He interaction potentials and contrasted with those for an effective low-energy potential model from the literature that reproduces the energies of $^4$He$_N$ clusters in the ground state for $N=2$ to $N=\infty$ at the $\gtrsim 95$~percent level with just four input parameters. The study is extended to unitarity (infinite $s$-wave scattering length) by artificially weakening the interaction potentials. In addition to contributing to the characterization of small bosonic helium quantum droplets, our study provides insights into the effective low-energy theory's predictability of various structural properties.

preprint2022arXiv

Two emitters coupled to a bath with Kerr-like non-linearity: Exponential decay, fractional populations, and Rabi oscillations

We consider two non-interacting two-level emitters that are coupled weakly to a one-dimensional non-linear wave guide. Due to the Kerr-like non-linearity, the wave guide considered supports -- in addition to the scattering continuum -- a two-body bound state. As such, the wave guide models a bath with non-trivial mode structure. Solving the time-dependent Schrödinger equation, the radiation dynamics of the two emitters, initially prepared in their excited states, is presented. Changing the emitter frequency such that the two-emitter energy is in resonance with one of the two-body bound states, radiation dynamics ranging from exponential decay to fractional populations to Rabi oscillations is observed. Along with the detuning, the dependence on the separation of the two emitters is investigated. Approximate reduced Hilbert space formulations, which result in effective emitter-separation and momentum dependent interactions, elucidate the underlying physical mechanisms and provide an avenue to showcase the features that would be absent if the one-dimensional wave guide did not contain a non-linearity. Our theoretical findings apply to a number of experimental platforms and the predictions can be tested with state-of-the-art technology. In addition, the weak-coupling Schrödinger equation based results provide critical guidance for the development of master equation approaches.

preprint2020arXiv

Mean-field spin-oscillation dynamics beyond the single-mode approximation for a harmonically trapped spin-1 Bose-Einstein condensate

Compared to single-component Bose-Einstein condensates, spinor Bose-Einstein condensates display much richer dynamics. In addition to density oscillations, spinor Bose-Einstein condensates exhibit intriguing spin dynamics that is associated with population transfer between different hyperfine components. This work analyzes the validity of the widely employed single-mode approximation when describing the spin dynamics in response to a quench of the system Hamiltonian. The single-mode approximation assumes that the different hyperfine states all share the same time-independent spatial mode. This implies that the resulting spin Hamiltonian only depends on the spin interaction strength and not on the density interaction strength. Taking the spinor sodium Bose-Einstein condensate in the $f=1$ hyperfine manifold as an example and working within the mean-field theory framework, it is found numerically that the single-mode approximation misses, in some parameter regimes, intricate details of the spin and spatial dynamics. We develop a physical picture that explains the observed phenomenon. Moreover, using that the population oscillations described by the single-mode approximation enter into the effective potential felt by the mean-field spinor, we derive a semi-quantitative condition for when dynamical mean-field induced corrections to the single-mode approximation are relevant. Our mean-field results have implications for a variety of published and planned experimental studies.

preprint2020arXiv

Rabi oscillations and Ramsey-type pulses in ultracold bosons: Role of interactions

Double-well systems loaded with one, two, or many quantum particles give rise to intriguing dynamics, ranging from Josephson oscillation to self-trapping. This work presents theoretical and experimental results for two distinct double-well systems, both created using dilute rubidium Bose-Einstein condensates with particular emphasis placed on the role of interaction in the systems. The first is realized by creating an effective two-level system through Raman coupling of hyperfine states. The second is an effective two-level system in momentum space generated through the coupling by an optical lattice. Even though the non-interacting systems can, for a wide parameter range, be described by the same model Hamiltonian, the dynamics for these two realizations differ in the presence of interactions. The difference is attributed to scattering diagrams that contribute in the lattice coupled system but vanish in the Raman coupled system. The internal dynamics of the Bose-Einstein condensates for both coupling scenarios is probed through a Ramsey-type interference pulse sequence, which constitutes a key building block of atom interferometers. These results have important implications in a variety of contexts including lattice calibration experiments and momentum space lattices used for quantum analog simulations.

preprint2016arXiv

Path integral Monte Carlo determination of the fourth-order virial coefficient for unitary two-component Fermi gas with zero-range interactions

The unitary equal-mass Fermi gas with zero-range interactions constitutes a paradigmatic model system that is relevant to atomic, condensed matter, nuclear, particle, and astro physics. This work determines the fourth-order virial coefficient $b_4$ of such a strongly-interacting Fermi gas using a customized \textit{ab initio} path integral Monte Carlo (PIMC) algorithm. In contrast to earlier theoretical results, which disagreed on the sign and magnitude of $b_4$, our $b_4$ agrees within error bars with the experimentally determined value, thereby resolving an ongoing literature debate. Utilizing a trap regulator, our PIMC approach determines the fourth-order virial coefficient by directly sampling the partition function. An on-the-fly anti-symmetrization avoids the Thomas collapse and, combined with the use of the exact two-body zero-range propagator, establishes an efficient general means to treat small Fermi systems with zero-range interactions.

preprint2016arXiv

Scattering framework for two particles with isotropic spin-orbit coupling applicable to all energies

Previous work developed a K-matrix formalism applicable to positive energies for the scattering between two $s$-wave interacting particles with two internal states, isotropic spin-orbit coupling and vanishing center-of-mass momentum [H. Duan, L. You and B. Gao, Phys. Rev A {\bf{87}}, 052708 (2013)]. This work extends the formalism to the entire energy regime. Explicit solutions are obtained for the total angular momentum $J=0$ and $1$ channels. The behavior of the partial cross sections in the negative energy regime is analyzed in detail. We find that the leading contributions to the partial cross sections at the negative energy thresholds are governed by the spin-orbit coupling strength $k_{\text{so}}$ and the mass ratio. The fact that these contributions are independent of the two-body scattering length $a_s$ is a direct consequence of the effective reduction of the dimensionality, and hence of the density of states, near the scattering thresholds due to the single-particle spin-orbit coupling terms. The results are analytically continued to the energy regime where bound states exist. It is shown that our results are consistent with results obtained by alternative approaches. Our formulation, which can be regarded as an extension of the standard textbook partial wave decomposition, can be generalized to two-body systems with other types of spin-orbit coupling, including cases where the center-of-mass momentum does not vanish.

preprint2015arXiv

Efimov physics and the three-body parameter for shallow van der Waals potentials

Extremely weakly-bound three-boson systems are predicted to exhibit intriguing universal properties such as discrete scale invariance. Motivated by recent experimental studies of the ground and excited helium trimers, this work analyzes the three-body parameter and the structural properties of three helium atoms as the s-wave scattering length is tuned artificially. Connections with theoretical and experimental studies of the Efimov scenario as it pertains to cold atom systems are made.

preprint2015arXiv

Energy and structural properties of $N$-boson clusters attached to three-body Efimov states: Two-body zero-range interactions and the role of the three-body regulator

The low-energy spectrum of $N$-boson clusters with pairwise zero-range interactions is believed to be governed by a three-body parameter. We study the ground state of $N$-boson clusters with infinite two-body $s$-wave scattering length by performing {\em{ab initio}} Monte Carlo simulations. To prevent Thomas collapse, different finite-range three-body regulators are used. The energy and structural properties for the three-body Hamiltonian with two-body zero-range interactions and three-body regulator are in much better agreement with the "ideal zero-range Efimov theory" results than those for Hamiltonian with two-body finite-range interactions. For larger clusters we find that the ground state energy and structural properties of the Hamiltonian with two-body zero-range interactions and finite-range three-body regulators are not universally determined by the three-body parameter, i.e., dependences on the specific form of the three-body regulator are observed. For comparison, we consider Hamiltonian with two-body van der Waals interactions and no three-body regulator. For the interactions considered, the ground state energy of the $N$-body clusters is---if scaled by the three-body ground state energy---fairly universal, i.e., the dependence on the short-range details of the two-body van der Waals potentials is small. Our results are compared with the literature.

preprint2015arXiv

Incorporating exact two-body propagators for zero-range interactions into $N$-body Monte Carlo simulations

Ultracold atomic gases are, to a very good approximation, described by pairwise zero-range interactions. This paper demonstrates that $N$-body systems with two-body zero-range interactions can be treated reliably and efficiently by the finite temperature and ground state path integral Monte Carlo approaches, using the exact two-body propagator for zero-range interactions in the pair product approximation. Harmonically trapped one- and three-dimensional systems are considered. A new propagator for the harmonically trapped two-body system with infinitely strong zero-range interaction, which may also have applications in real time evolution schemes, is presented.

preprint2015arXiv

One-dimensional Fermi gas with a single impurity in a harmonic trap: Perturbative description of the upper branch

The transition from "few to many" has recently been probed experimentally in an ultra cold harmonically confined one-dimensional lithium gas, in which a single impurity atom interacts with a background gas consisting of one, two, or more identical fermions [A. N. Wenz {\em{et al.}}, Science {\bf{342}}, 457 (2013)]. For repulsive interactions between the background or majority atoms and the impurity, the interaction energy for relatively moderate system sizes was analyzed and found to converge toward the corresponding expression for an infinitely large Fermi gas. Motivated by these experimental results, we apply perturbative techniques to determine the interaction energy for weak and strong coupling strengths and derive approximate descriptions for the interaction energy for repulsive interactions with varying strength between the impurity and the majority atoms and any number of majority atoms.

preprint2015arXiv

Spin structure of harmonically trapped one-dimensional atoms with spin-orbit coupling

We introduce a theoretical approach to determine the spin structure of harmonically trapped atoms with two-body zero-range interactions subject to an equal mixture of Rashba and Dresselhaus spin-orbit coupling created through Raman coupling of atomic hyperfine states. The spin structure of bosonic and fermionic two-particle systems with finite and infinite two-body interaction strength $g$ is calculated. Taking advantage of the fact that the $N$-boson and $N$-fermion systems with infinitely large coupling strength $g$ are analytically solvable for vanishing spin-orbit coupling strength $k_{so}$ and vanishing Raman coupling strength $Ω$, we develop an effective spin model that is accurate to second-order in $Ω$ for any $k_{so}$ and infinite $g$. The three- and four-particle systems are considered explicitly. It is shown that the effective spin Hamiltonian, which contains a Heisenberg exchange term and an anisotropic Dzyaloshinskii-Moriya exchange term, describes the transitions that these systems undergo with the change of $k_{so}$ as a competition between independent spin dynamics and nearest-neighbor spin interactions.

preprint2015arXiv

Trapped unitary two-component Fermi gases with up to ten particles

The properties of two-component Fermi gases with zero-range interactions are universal. We use an explicitly correlated Gaussian basis set expansion approach to investigate small equal-mass two-component Fermi gases under spherically symmetric external harmonic confinement. At unitarity, we determine the ground state energy for systems with up to ten particles interacting through finite-range two-body potentials for both even and odd number of particles. We extrapolate the energies to the zero-range limit using a novel scheme that removes the linear and, in some cases, also the quadratic dependence of the ground state energies on the two-body range. Our extrapolated zero-range energies are compared with results from the literature. We also calculate the two-body Tan contact and structural properties.

preprint2015arXiv

Tunneling dynamics of two interacting one-dimensional particles

We present one-dimensional simulation results for the cold atom tunneling experiments by the Heidelberg group [G. Zürn {\em{et al.}}, Phys. Rev. Lett. {\bf{108}}, 075303 (2012) and G. Zürn {\em{et al.}}, Phys. Rev. Lett. {\bf{111}}, 175302 (2013)] on one or two $^6$Li atoms confined by a potential that consists of an approximately harmonic optical trap plus a linear magnetic field gradient. At the non-interacting particle level, we find that the WKB (Wentzel-Kramers-Brillouin) approximation may not be used as a reliable tool to extract the trapping potential parameters from the experimentally measured tunneling data. We use our numerical calculations along with the experimental tunneling rates for the non-interacting system to reparameterize the trapping potential. The reparameterized trapping potentials serve as input for our simulations of two interacting particles. For two interacting (distinguishable) atoms on the upper branch, we reproduce the experimentally measured tunneling rates, which vary over several orders of magnitude, fairly well. For infinitely strong interaction strength, we compare the time dynamics with that of two identical fermions and discuss the implications of fermionization on the dynamics. For two attractively-interacting atoms on the molecular branch, we find that single-particle tunneling dominates for weakly-attractive interactions while pair tunneling dominates for strongly-attractive interactions. Our first set of calculations yields qualitative but not quantitative agreement with the experimentally measured tunneling rates. We obtain quantitative agreement with the experimentally measured tunneling rates if we allow for a weakened radial confinement.

preprint2014arXiv

Energy spectrum of a harmonically trapped two-atom system with spin-orbit coupling

Ultracold atomic gases provide a novel platform with which to study spin-orbit coupling, a mechanism that plays a central role in the nuclear shell model, atomic fine structure and two-dimensional electron gases. This paper introduces a theoretical framework that allows for the efficient determination of the eigenenergies and eigenstates of a harmonically trapped two-atom system with short-range interaction subject to an equal mixture of Rashba and Dresselhaus spin-orbit coupling created through Raman coupling of atomic hyperfine states. Energy spectra for experimentally relevant parameter combinations are presented and future extensions of the approach are discussed.

preprint2014arXiv

Generalized Efimov scenario for heavy-light mixtures

Motivated by recent experimental investigations of Cs-Cs-Li Efimov resonances, this work theoretically investigates the few-body properties of $N-1$ non-interacting identical heavy bosons, which interact with a light impurity through a large $s$-wave scattering length. For Cs-Cs-Cs-Li, we predict the existence of universal four-body states with energies $E_4^{(n,1)}$ and $E_4^{(n,2)}$, which are universally linked to the energy $E_3^{(n)}$ of the $n$th Efimov trimer. For infinitely large $^{133}$Cs-$^6$Li and vanishing $^{133}$Cs-$^{133}$Cs scattering lengths, we find $(E_4^{(1,1)}/E_3^{(1)})^{1/2} \approx 1.51$ and $(E_4^{(1,2)}/E_3^{(1)})^{1/2} \approx 1.01$. The $^{133}$Cs-$^6$Li scattering lengths at which these states merge with the four-atom threshold, the dependence of these energy ratios on the mass ratio between the heavy and light atoms, and selected aspects of the generalized Efimov scenario for $N>4$ are also discussed. Possible implications of our results for ongoing cold atom experiments are presented.

preprint2014arXiv

Harmonically trapped two-atom systems: Interplay of short-range s-wave interaction and spin-orbit coupling

The coupling between the spin degrees of freedom and the orbital angular momentum has a profound effect on the properties of nuclei, atoms and condensed matter systems. Recently, synthetic gauge fields have been realized experimentally in neutral cold atom systems, giving rise to a spin-orbit coupling term with "strength" $k_{\text{so}}$. This paper investigates the interplay between the single-particle spin-orbit coupling term of Rashba type and the short-range two-body $s$-wave interaction for cold atoms under external confinement. Specifically, we consider two different harmonically trapped two-atom systems. The first system consists of an atom with spin-orbit coupling that interacts with a structureless particle through a short-range two-body potential. The second system consists of two atoms that both feel the spin-orbit coupling term and that interact through a short-range two-body potential. Treating the spin-orbit term perturbatively, we determine the correction to the ground state energy for various generic parameter combinations. Selected excited states are also treated. An important aspect of our study is that the perturbative treatment is not limited to small $s$-wave scattering lengths but provides insights into the system behavior over a wide range of scattering lengths, including the strongly-interacting unitary regime. Our perturbative results are confirmed by numerical calculations that expand the eigenfunctions of the two-particle Hamiltonian in terms of basis functions that contain explicitly correlated Gaussians.

preprint2014arXiv

Molecular branch of a small highly-elongated Fermi gas with an impurity: Full three-dimensional versus effective one-dimensional description

We consider an impurity immersed in a small Fermi gas under highly-elongated harmonic confinement. The impurity interacts with the atoms of the Fermi gas through an isotropic short-range potential with three-dimensional free-space s-wave scattering length $a_{3\text{d}}$. We investigate the energies of the molecular branch, i.e., the energies of the state that corresponds to a gas consisting of a weakly-bound diatomic molecule and "unpaired" atoms, as a function of the s-wave scattering length $a_{3\text{d}}$ and the ratio $η$ between the angular trapping frequencies in the tight and weak confinement directions. The energies obtained from our three-dimensional description that accounts for the dynamics in the weak and tight confinement directions are compared with those obtained within an effective one-dimensional framework, which accounts for the dynamics in the tight confinement direction via a renormalized one-dimensional coupling constant. Our theoretical results are related to recent experimental measurements.

preprint2014arXiv

Temperature-dependence of small harmonically trapped atom systems with Bose, Fermi and Boltzmann statistics

While the zero-temperature properties of harmonically trapped cold few-atom systems have been discussed fairly extensively over the past decade, much less is known about the finite-temperature properties. Working in the canonical ensemble, we characterize small harmonically trapped atomic systems as a function of the temperature using analytical and numerical techniques. We present results for the energetics, structural properties, condensate fraction, superfluid fraction, and superfluid density. Our calculations for the two-body system underline that the condensate and superfluid fractions are distinctly different quantities. Our work demonstrates that the path integral Monte Carlo method yields reliable results for bosonic and fermionic systems over a wide temperature range, including the regime where the de Broglie wave length is large, i.e., where the statistics plays an important role. The regime where the Fermi sign problem leads to reasonably large signal to noise ratios is mapped out for selected parameter combinations. Our calculations for bosons focus on the unitary regime, where the physics is expected to be governed by the three-body parameter. If the three-body parameter is large compared to the inverse of the harmonic oscillator length, we find that the bosons form a droplet at low temperature and behave approximately like a non-interacting Bose and eventually Boltzmann gas at high temperature. The change of the behavior occurs over a fairly narrow temperature range. A simple model that reproduces the key aspects of the phase transition like feature, which can potentially be observed in cold atom Bose gas experiments, is presented.

preprint2014arXiv

Three-body bound states in a harmonic waveguide with cylindrical symmetry

Highly-elongated quasi-one-dimensional cold atom samples have been studied extensively over the past years experimentally and theoretically. This work determines the energy spectrum of two identical fermions and a third distinguishable particle as functions of the mass ratio $κ$ and the free-space $s$-wave scattering length $a_{3\text{D}}$ between the identical fermions and the distinguishable third particle in a cylindrically symmetric waveguide whose symmetry axis is chosen to be along the $z$-axis. We focus on the regime where the mass of the identical fermions is equal to or larger than that of the third distinguishable particle. Our theoretical framework accounts explicitly for the motion along the transverse confinement direction. In the regime where excitations in the transverse direction are absent (i.e., for states with projection quantum number $M_{\text{rel}}=0$), we determine the binding energies for states with odd parity in $z$. These full three-dimensional energies deviate significantly from those obtained within a strictly one-dimensional framework when the $s$-wave scattering length is of the order of or smaller than the oscillator length in the confinement direction. If transverse excitations are present, we predict the existence of a new class of universal three-body bound states with $|M_{\text{rel}}|=1$ and positive parity in $z$. These bound states arise on the positive $s$-wave scattering length side if the mass ratio $κ$ is sufficiently large. Implications of our results for ongoing cold atom experiments are discussed.

preprint2014arXiv

Universal and non-universal effective $N$-body interactions for ultracold harmonically-trapped few-atom systems

We derive the ground-state energy for a small number of ultracold atoms in an isotropic harmonic trap using effective quantum field theory (EFT). Atoms are assumed to interact through pairwise energy-independent and energy-dependent delta-function potentials with strengths proportional to the scattering length $a$ and effective range volume $V$, respectively. The calculations are performed systematically up to order $l^{-4}$, where $l$ denotes the harmonic oscillator length. The effective three-body interaction contains a logarithmic divergence in the cutoff energy, giving rise to a non-universal three-body interaction in the EFT. Our EFT results are confirmed by nonperturbative numerical calculations for a Hamiltonian with finite-range two-body Gaussian interactions. For this model Hamiltonian, we explicitly calculate the non-universal effective three-body contribution to the energy.

preprint2013arXiv

Abnormal Superfluid Fraction of Harmonically Trapped Few-Fermion Systems

Superfluidity is a fascinating phenomenon that, at the macroscopic scale, leads to dissipationless flow and the emergence of vortices. While these macroscopic manifestations of superfluidity are well described by theories that have their origin in Landau's two-fluid model, our microscopic understanding of superfluidity is far from complete. Using analytical and numerical \textit{ab initio} approaches, this paper determines the superfluid fraction and local superfluid density of small harmonically trapped two-component Fermi gases as a function of the interaction strength and temperature. At low temperature, we find that the superfluid fraction is, in certain regions of the parameter space, negative. This counterintuitive finding is traced back to the symmetry of the system's ground state wave function, which gives rise to a diverging quantum moment of inertia $I_{\text{q}}$. Analogous abnormal behavior of $I_{\text{q}}$ has been observed in even-odd nuclei at low temperature. Our predictions can be tested in modern cold atom experiments.

preprint2013arXiv

Correlations of the upper branch of 1d harmonically trapped two-component Fermi gases

We present highly-accurate energy spectra and eigen functions of small 1d harmonically trapped two-component Fermi gases with interspecies $δ$-function interactions, and analyze the correlations of the so-called upper branch (i.e., the branch that describes a repulsive Fermi gas consisting of atoms but no molecules) for positive and negative coupling constants. Changes of the two-body correlations as a function of the interspecies coupling strength reflect the competition of the interspecies interaction and the effective repulsion due to the Pauli exclusion principle, and are interpreted as a few-body analog of a transition from a non-magnetic to a magnetic phase. Moreover, we show that the eigenstate $ψ_{\rm{adia}}$ of the infinitely strongly-interacting system with $|n_1+n_2|>2$ and $|n_1-n_2| < n$ ($n_1$ and $n_2$ denote the number of fermions of components 1 and 2, respectively), which is reached experimentally by adiabatically changing the system parameters, does not, as previously proposed, coincide with the wave function $ψ_{\rm{G}}$ obtained by applying a generalized Fermi-Fermi mapping function to the eigen function of the non-interacting single-component Fermi gas.

preprint2013arXiv

Harmonically trapped Fermi gas: Temperature dependence of the Tan contact

Ultracold atomic gases with short-range interactions are characterized by a number of universal species-independent relations. Many of these relations involve the two-body Tan contact. Employing the canonical ensemble, we determine the Tan contact for small harmonically trapped two-component Fermi gases at unitarity over a wide range of temperatures, including the zero and high temperature regimes. A cluster expansion that describes the properties of the N-particle system in terms of those of smaller subsystems is introduced and shown to provide an accurate description of the contact in the high temperature regime. Finite-range corrections are quantified and the role of the Fermi statistics is elucidated by comparing results for Fermi, Bose and Boltzmann statistics.

preprint2013arXiv

Non-universal bound states of two identical heavy fermions and one light particle

We study the behavior of the bound state energy of a system consisting of two identical heavy fermions of mass M and a light particle of mass m. The heavy fermions interact with the light particle through a short-range two-body potential with positive s-wave scattering length a_s. We impose a short-range boundary condition on the logarithmic derivative of the hyperradial wavefunction and show that, in the regime where Efimov states are absent, a non-universal three-body state "cuts through" the universal three-body states previously described by Kartavtsev and Malykh [O. I. Kartavtsev and A. V. Malykh, J. Phys. B 40, 1429 (2007)]. The presence of the non-universal state alters the behavior of the universal states in certain regions of the parameter space. We show that the existence of the non-universal state is predicted accurately by a simple quantum defect theory model that utilizes hyperspherical coordinates. An empirical two-state model is employed to quantify the coupling of the non-universal state to the universal states.

preprint2013arXiv

Small two-component Fermi gases in a cubic box with periodic boundary conditions

The properties of two-component Fermi gases become universal if the interspecies s-wave scattering length $a_s$ and the average interparticle spacing are much larger than the range of the underlying two-body potential. Using an explicitly correlated Gaussian basis set expansion approach, we determine the eigen energies of two-component Fermi gases in a cubic box with periodic boundary conditions as functions of the interspecies s-wave scattering length and the effective range of the two-body potential. The universal properties of systems consisting of up to four particles are determined by extrapolating the finite-range energies to the zero-range limit. We determine the eigen energies of states with vanishing and finite momentum. In the weakly-attractive BCS regime, we analyze the energy spectra and degeneracies using first-order degenerate perturbation theory. Excellent agreement between the perturbative energy shifts and the numerically determined energies is obtained. For the infinitely large scattering length case, we compare our results - where available - with those presented in the literature.

preprint2013arXiv

Tunable high-temperature thermodynamics of weakly-interacting dipolar gases

We consider dilute gases of dipolar bosons or fermions in the high-temperature limit in a spherically symmetric harmonic trapping potential. We examine the system using a virial expansion up to second order in the fugacity. Using the Born approximation and assuming purely dipolar interactions, we find that the second-order virial coefficient for both bosons and fermions depends quadratically on the dipole length and is negative at high temperature, indicating that to lowest order in the dipole-dipole interactions the dipolar single-component quantum gases are repulsive. If the $s$-wave scattering length for the bosonic system is tunable and its absolute value is made small, then the $s$-wave interactions dominate and the dipolar as behaves like a weakly-interacting Bose gas with isotropic $s$-wave interactions. If the generalized scattering lengths for the fermionic system are tunable, then the dipole length can enter linearly in the virial equation of state, enhancing the dipole-dipole effects in the thermodynamic observables.

preprint2012arXiv

Degeneracies in trapped two-component Fermi gases

We report on previously unobserved inter-system degeneracies in two-component equal-mass Fermi gases with interspecies zero-range interactions under isotropic harmonic confinement. Over the past 10 years, two-component Fermi gases consisting of $n_1$ spin-up and $n_2$ spin-down atoms with interspecies zero-range interactions have become a paradigm for modeling condensed matter systems, nuclear matter and neutron matter. We show that the eigen energies of the $(n_1+1,n_2-1)$ system are degenerate with the eigen energies of the $(n_1,n_2)$ system for any s-wave scattering length $a_s$, including infinitely large, positive and negative $a_s$. The existence of the inter-system degeneracies is demonstrated explicitly for few-body systems with $n_1+n_2=4, 5$ and 6. The degeneracies and associated symmetries are explained within a group theoretical framework.

preprint2012arXiv

Effective renormalized multi-body interactions of harmonically confined ultracold neutral bosons

We calculate the renormalized effective 2-, 3-, and 4-body interactions for N neutral ultracold bosons in the ground state of an isotropic harmonic trap, assuming 2-body interactions modeled with the combination of a zero-range and energy-dependent pseudopotential. We work to third-order in the scattering length a defined at zero collision energy, which is necessary to obtain both the leading-order effective 4-body interaction and consistently include finite-range corrections for realistic 2-body interactions. The leading-order, effective 3- and 4-body interaction energies are U3 = -(0.85576...)(a/l)^2 + 2.7921(1)(a/l)^3 + O[(a/l)^4] and U4 = +(2.43317...)(a/l)^3 + O[(ał)^4], where w and l are the harmonic oscillator frequency and length, respectively, and energies are in units of hbar*w. The one-standard deviation error 0.0001 for the third-order coefficient in U3 is due to numerical uncertainty in estimating a slowly converging sum; the other two coefficients are either analytically or numerically exact. The effective 3- and 4-body interactions can play an important role in the dynamics of tightly confined and strongly correlated systems. We also performed numerical simulations for a finite-range boson-boson potential, and it was comparison to the zero-range predictions which revealed that finite-range effects must be taken into account for a realistic third-order treatment. In particular, we show that the energy-dependent pseudopotential accurately captures, through third order, the finite-range physics, and in combination with the multi-body effective interactions gives excellent agreement with the numerical simulations, validating our theoretical analysis and predictions.

preprint2012arXiv

Hyperspherical explicitly correlated Gaussian approach for few-body systems with finite angular momentum

Within the hyperspherical framework, the solution of the time-independent Schroedinger equation for a n-particle system is divided into two steps, the solution of a Schroedinger like equation in the hyperangular degrees of freedom and the solution of a set of coupled Schroedinger like hyperradial equations. The solutions to the former provide effective potentials and coupling matrix elements that enter into the latter set of equations. This paper develops a theoretical framework to determine the effective potentials, as well as the associated coupling matrix elements, for few-body systems with finite angular momentum L=1 and negative and positive parity. The hyperangular channel functions are expanded in terms of explicitly correlated Gaussian basis functions and relatively compact expressions for the matrix elements are derived. The developed formalism is applicable to any n; however, for n greater or equal to 6, the computational demands are likely beyond present-day computational capabilities. A number of calculations relevant to cold atom physics are presented, demonstrating that the developed approach provides a computationally efficient means to solving four-body bound and scattering problems with finite angular momentum on powerful desktop computers. Details regarding the implementation are discussed.

preprint2012arXiv

Natural and unnatural parity states of small trapped equal-mass two-component Fermi gases at unitarity and fourth-order virial coefficient

Equal-mass two-component Fermi gases under spherically symmetric external harmonic confinement with large s-wave scattering length are considered. Using the stochastic variational approach, we determine the lowest 286 and 164 relative eigenenergies of the (2,2) and (3,1) systems at unitarity as a function of the range $r_0$ of the underlying two-body potential and extrapolate to the $r_0 \rightarrow 0$ limit. Our calculations include all states with vanishing and finite angular momentum $L$ (and natural and unnatural parity $Π$) with relative energy up to $10.5 \hbar Ω$, where $Ω$ denotes the angular trapping frequency of the external confinement. Our extrapolated zero-range energies are estimated to have uncertainties of 0.1% or smaller. The (2,2) and (3,1) energies are used to determine the fourth-order virial coefficient of the trapped unitary two-component Fermi gas in the low-temperature regime. Our results are compared with recent predictions for the fourth-order virial coefficient of the homogeneous system. We also calculate small portions of the energy spectra of the (3,2) and (4,1) systems at unitarity.

preprint2012arXiv

Occupation numbers of the harmonically trapped few-boson system

We consider a harmonically trapped dilute $N$-boson system described by a low-energy Hamiltonian with pairwise interactions. We determine the condensate fraction, defined in terms of the largest occupation number, of the weakly-interacting $N$-boson system ($N \ge 2$) by employing a perturbative treatment within the framework of second quantization. The one-body density matrix and the corresponding occupation numbers are compared with those obtained by solving the two-body problem with zero-range interactions exactly. Our expressions are also compared with high precision {\em{ab initio}} calculations for Bose gases with $N=2-4$ that interact through finite-range two-body model potentials. Non-universal corrections are identified to enter at subleading order, confirming that different low-energy Hamiltonians, constructed to yield the same energy, may yield different occupation numbers. Lastly, we consider the strongly-interacting three-boson system under spherically symmetric harmonic confinement and determine its occupation numbers as a function of the three-body "Efimov parameter".

preprint2012arXiv

Thermodynamics of the two-component Fermi gas with unequal masses at unitarity

We consider mass-imbalanced two-component Fermi gases for which the unequal-mass atoms interact via a zero-range model potential with a diverging s-wave scattering length $a_s$, i.e., with $1/a_s=0$. The high temperature thermodynamics of the harmonically trapped and homogeneous systems are examined using a virial expansion approach up to third order in the fugacity. We find that the universal part of the third-order virial coefficient associated with two light atoms and one heavy atom is negative, while that associated with two heavy and one light atom changes sign from negative to positive as the mass ratio $κ$ increases, and diverges when Efimov physics sets in at $κ=13.61$. By examining the Helmholtz free energy, we find that the equilibrium polarization of the trapped and homogeneous systems is 0 for $κ=1$, but finite for $κ\ne 1$ (with a majority of heavy particles). Compared to the equilibrium polarization of the non-interacting system, the equilibrium polarization at unitarity is increased for the trapped system and decreased for the homogeneous system. We find that unequal-mass Fermi gases are stable for all polarizations.

preprint2012arXiv

Three s-wave interacting fermions under anisotropic harmonic confinement: Dimensional crossover of energetics and virial coefficients

We present essentially exact solutions of the Schroedinger equation for three fermions in two different spin states with zero-range s-wave interactions under harmonic confinement. Our approach covers spherically symmetric, strictly two-dimensional, strictly one-dimensional, cigar-shaped, and pancake-shaped traps. In particular, we discuss the transition from quasi-one-dimensional to strictly one-dimensional and from quasi-two-dimensional to strictly two-dimensional geometries. We determine and interpret the eigenenergies of the system as a function of the trap geometry and the strength of the zero-range interactions. The eigenenergies are used to investigate the dependence of the second- and third-order virial coefficients, which play an important role in the virial expansion of the thermodynamic potential, on the geometry of the trap. We show that the second- and third-order virial coefficients for anisotropic confinement geometries are, for experimentally relevant temperatures, very well approximated by those for the spherically symmetric confinement for all s-wave scattering lengths.

preprint2012arXiv

Universal four-body states in heavy-light mixtures with positive scattering length

The number of four-body states known to behave universally is small. This work adds a new class of four-body states to this relatively short list. We predict the existence of a universal four-body bound state for heavy-light mixtures consisting of three identical heavy fermions and a fourth distinguishable lighter particle with mass ratio $κ\gtrsim 9.5$ and short-range interspecies interaction characterized by a positive s-wave scattering length. The structural properties of these universal states are discussed and finite-range effects are analyzed. The bound states can be experimentally realized and probed utilizing ultracold atom mixtures.

preprint2011arXiv

Few-body physics with ultracold atomic and molecular systems in traps

Few-body physics has played a prominent role in atomic, molecular and nuclear physics since the early days of quantum mechanics. It is now possible---thanks to tremendous progress in cooling, trapping, and manipulating ultracold samples---to experimentally study few-body phenomena in trapped atomic and molecular systems with unprecedented control. This review summarizes recent studies of few-body phenomena in trapped atomic and molecular gases, with an emphasis on small trapped systems. We start by introducing the free-space scattering properties and then investigate what happens when two particles, bosons or fermions, are placed in an external confinement. Next, various three-body systems are treated analytically in limiting cases. Our current understanding of larger two-component Fermi systems and Bose systems is reviewed, and connections with the corresponding bulk systems are established. Lastly, future prospects and challenges are discussed. Throughout this review, commonalities with other systems such as nuclei or quantum dots are highlighted.

preprint2011arXiv

Modification of roton instability due to the presence of a second dipolar Bose-Einstein condensate

We study the behavior of two coupled purely dipolar Bose-Einstein condensates, each located in a cylindrically symmetric pancake-shaped external confining potential, as the separation b between the traps along the tight confining direction is varied. The solutions of the coupled Gross-Pitaevskii and Bogoliubov-de Gennes equations, which account for the full dynamics, show that the system behavior is modified by the presence of the second dipolar BEC. For sufficiently small b, the presence of the second dipolar BEC destabilizes the system dramatically. In this regime, the coupled system collapses through a mode that is notably different from the radial roton mode that induces the collapse of the uncoupled system. Finally, we comment on the shortcomings of an approach that neglects the dynamics in the z-direction, which is assumed to be a good approximation for highly pancake-shaped dipolar BECs in the literature.

preprint2010arXiv

Breakdown of universality for unequal-mass Fermi gases with infinite scattering length

We treat small trapped unequal-mass two-component Fermi gases at unitarity within a non-perturbative microscopic framework and investigate the system properties as functions of the mass ratio $κ$, and the numbers $N_{1}$ and $N_2$ of heavy and light fermions. While equal-mass Fermi gases with infinitely large interspecies s-wave scattering length $a_s$ are universal, we find that unequal mass Fermi gases are, for sufficiently large $κ$ and in the regime where Efimov physics is absent, not universal. In particular, the $(N_1,N_2)=(2,1)$ and $(3,1)$ systems exhibit three-body (3b) and four-body (4b) resonances at $κ= 12.314(2)$ and $10.4(2)$, respectively, as well as surprisingly large finite-range (FR) effects. These findings have profound implications for ongoing experimental efforts and quantum simulation proposals that utilize unequal-mass atomic Fermi gases.

preprint2010arXiv

Energy spectrum of harmonically trapped two-component Fermi gases: Three- and Four-Particle Problem

Trapped two-component Fermi gases allow for the investigation of the so-called BCS-BEC crossover by tuning the interspecies atom-atom $s$-wave scattering length scattering $a^{(aa)}$ from attractive to repulsive, including vanishing and infinitely large values. Here, we numerically determine the energy spectrum of the equal-mass spin-balanced four-fermion system---the smallest few-particle system that exhibits BCS-BEC crossover-like behavior---as a function of $a^{(aa)}$ using the stochastic variational approach. For comparative purposes, we also treat the two- and three-particle systems. States with vanishing and finite total angular momentum as well as with natural and unnatural parity are considered. In addition, the energy spectrum of weakly-attractive and weakly-repulsive gases is characterized by employing a perturbative framework that utilizes hyperspherical coordinates. The hyperspherical coordinate approach allows for the straightforward assignment of quantum numbers and furthermore provides great insights into the strongly-interacting unitary regime.

preprint2010arXiv

Few-body resonances of unequal-mass systems with infinite interspecies two-body s-wave scattering length

Two-component Fermi and Bose gases with infinitely large interspecies s-wave scattering length $a_s$ exhibit a variety of intriguing properties. Among these are the scale invariance of two-component Fermi gases with equal masses, and the favorable scaling of Efimov features for two-component Bose gases and Bose-Fermi mixtures with unequal masses. This paper builds on our earlier work [D. Blume and K. M. Daily, arXiv:1006.5002] and presents a detailed discussion of our studies of small unequal-mass two-component systems with infinite $a_s$ in the regime where three-body Efimov physics is absent. We report on non-universal few-body resonances. Just like with two-body systems on resonance, few-body systems have a zero-energy bound state in free space and a diverging generalized scattering length. Our calculations are performed within a non-perturbative microscopic framework and investigate the energetics and structural properties of small unequal-mass two-component systems as functions of the mass ratio $κ$, and the numbers $N_{1}$ and $N_2$ of heavy and light atoms. For purely attractive Gaussian two-body interactions, we find that the $(N_1,N_2)=(2,1)$ and $(3,1)$ systems exhibit three-body and four-body resonances at mass ratios $κ= 12.314(2)$ and 10.4(2), respectively. The three- and four-particle systems on resonance are found to be large. This suggests that the corresponding wave function has relatively small overlap with deeply-bound dimers, trimers or larger clusters and that the three- and four-body systems on resonance have a comparatively long lifetime. Thus, it seems feasible that the features discussed in this paper can be probed experimentally with present-day technology.

preprint2010arXiv

Trapped two-component Fermi gases with up to six particles: Energetics, structural properties, and molecular condensate fraction

We investigate small equal-mass two-component Fermi gases under external spherically symmetric confinement in which atoms with opposite spins interact through a short-range two-body model potential. We employ a non-perturbative microscopic framework, the stochastic variational approach, and determine the system properties as functions of the interspecies s-wave scattering length a, the orbital angular momentum L of the system, and the numbers N1 and N2 of spin-up and spin-down atoms (with N1-N2 =0 or 1 and N < 7, where N=N1+N2). At unitarity, we determine the energies of the five- and six-particle systems for various ranges r0 of the underlying two-body model potential and extrapolate to the zero-range limit. These energies serve as benchmark results that can be used to validate and assess other numerical approaches. We also present structural properties such as the pair distribution function and the radial density. Furthermore, we analyze the one-body and two-body density matrices. A measure for the molecular condensate fraction is proposed and applied. Our calculations show explicitly that the natural orbitals and the momentum distributions of atomic Fermi gases approach those characteristic for a molecular Bose gas if the s-wave scattering length a, a>0, is sufficiently small.

preprint2009arXiv

Aligned dipolar Bose-Einstein condensate in a double-well potential: From cigar-shaped to pancake-shaped

We consider a Bose-Einstein condensate (BEC), which is characterized by long-range and anisotropic dipole-dipole interactions and vanishing s-wave scattering length, in a double-well potential. The properties of this system are investigated as functions of the height of the barrier that splits the harmonic trap into two halves, the number of particles (or dipole-dipole strength) and the aspect ratio $λ$, which is defined as the ratio between the axial and longitudinal trapping frequencies $ω_z$ and $ω_ρ$. The phase diagram is determined by analyzing the stationary mean-field solutions. Three distinct regions are found: a region where the energetically lowest lying stationary solution is symmetric, a region where the energetically lowest lying stationary solution is located asymmetrically in one of the wells, and a region where the system is mechanically unstable. For sufficiently large aspect ratio $λ$ and sufficiently high barrier height, the system consists of two connected quasi-two-dimensional sheets with density profiles whose maxima are located either at $ρ=0$ or away from $ρ=0$. The stability of the stationary solutions is investigated by analyzing the Bogoliubov de Gennes excitation spectrum and the dynamical response to small perturbations. These studies reveal unique oscillation frequencies and distinct collapse mechanisms. The results derived within the mean-field framework are complemented by an analysis based on a two-mode model.

preprint2009arXiv

Excitation spectrum and effective interactions of highly-elongated Fermi gas

Full 3D calculations of small two-component Fermi gases under highly-elongated confinement, in which unlike fermions interact through short-range potentials with variable atom-atom s-wave scattering length, are performed using the correlated Gaussian approach. In addition, microscopic 1D calculations are performed for effective "atomic" and "molecular" 1D model Hamiltonian. Comparisons of the 3D and 1D energies and excitation frequencies establish the validity regimes of the effective 1D Hamiltonian. Our numerical results for three- and four-particle systems suggest that the effective 1D atom-dimer and dimer-dimer interactions are to a good approximation determined by simple analytical expressions. Implications for the description of quasi-1D Fermi gases within strict 1D frameworks are discussed.

preprint2009arXiv

Illustration of universal relations for trapped four-fermion system with arbitrary s-wave scattering length

A two-component four-fermion system with equal masses, interspecies s-wave scattering length a and vanishing intraspecies interactions under external spherically symmetric harmonic confinement is considered. Using a correlated Gaussian basis set expansion approach, we determine the energies and various structural properties of the energetically lowest-lying gas-like state throughout the crossover for various ranges of the underlying two-body potential. Extrapolating to the zero-range limit, our numerical results show explicitly that the total energy, the trap energy as well as certain aspects of the pair distribution function and of the momentum distribution are related through the so-called integrated contact intensity I(a). Furthermore, it is shown explicitly that the total energy and the trap energy are related through a generalized virial theorem that accounts for a non-zero range.

preprint2009arXiv

Tuning the structural and dynamical properties of a dipolar Bose-Einstein condensate: Ripples and instability islands

It is now well established that the stability of aligned dipolar Bose gases can be tuned by varying the aspect ratio of the external harmonic confinement. This paper extends this idea and demonstrates that a Gaussian barrier along the strong confinement direction can be employed to tune both the structural properties and the dynamical stability of an oblate dipolar Bose gas aligned along the strong confinement direction. In particular, our theoretical mean-field analysis predicts the existence of instability islands immersed in otherwise stable regions of the phase diagram. Dynamical studies indicate that these instability islands, which can be probed experimentally with present-day technology, are associated with the going soft of a Bogoliubov--de Gennes excitation frequency with radial breathing mode character. Furthermore, we find dynamically stable ground state densities with ripple-like oscillations along the radial direction. These structured ground states exist in the vicinity of a dynamical radial roton-like instability.

preprint2006arXiv

Dipolar Bose gases: Many-body versus mean-field description

We characterize zero-temperature dipolar Bose gases under external spherical confinement as a function of the dipole strength using the essentially exact many-body diffusion Monte Carlo (DMC) technique. We show that the DMC energies are reproduced accurately within a mean-field framework if the variation of the s-wave scattering length with the dipole strength is accounted for properly. Our calculations suggest stability diagrams and collapse mechanisms of dipolar Bose gases that differ significantly from those previously proposed in the literature.