BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Peter D. Drummond

Peter D. Drummond appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.quant-gasquant-phcond-mat.str-elcond-mat.supr-congr-qcmath.NAastro-ph.COcond-mat.stat-mechmath-phmath.MPNumerical Analysisphysics.comp-ph

Catalog footprint

What is connected

31works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2022arXiv

A midpoint projection algorithm for stochastic differential equations on manifolds

Stochastic differential equations projected onto manifolds occur in physics, chemistry, biology, engineering, nanotechnology and optimization, with interdisciplinary applications. Intrinsic coordinate stochastic equations on the manifold are often computationally impractical, and numerical projections are useful in many cases. We show that the Stratonovich interpretation of the stochastic calculus is obtained using adiabatic elimination with a constraint potential. We derive intrinsic stochastic equations for spheroidal and hyperboloidal surfaces for comparison purposes, and review some earlier projection algorithms. In this paper, a combined midpoint projection algorithm is proposed that uses a midpoint projection onto a tangent space, combined with a subsequent normal projection to satisfy the constraints. Numerical examples are given for a range of manifolds, including circular, spheroidal, hyperboloidal, and catenoidal cases, as well as higher-order polynomial constraints and a ten-dimensional hypersphere. We show that in all cases the combined midpoint method has greatly reduced errors compared to methods using a combined Euler projection approach or purely tangential projection. Our technique can handle multiple constraints. This allows manifolds that embody several conserved quantities. The algorithm is accurate, simple and efficient. An order of magnitude error reduction in diffusion distance is typically found compared to the other methods, with reductions of several orders of magnitude in constraint errors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Massive particle interferometry with lattice solitons: robustness against ionization

We revisit the proposal of Castin and Weiss [Phys. Rev. Lett. vol. 102, 010403 (2009)] for using the scattering of a quantum matter-wave soliton on a barrier in order to create a coherent superposition state of the soliton being entirely to the left of the barrier and being entirely to the right of the barrier. In that proposal, is was assumed that the scattering is perfectly elastic, i.e. that the center-of-mass kinetic energy of the soliton is lower than the chemical potential of the soliton. Here we relax this assumption. Also, we introduce an interferometric scheme, which uses interference of soltions, that can be used to detect the degree of coherence between the reflected and transmitted part of the soliton. Using exact diagonalization, we numerically simulate a complete interferometric cycle for a soliton consisting of six atoms. We find that the interferometric fringes persist even when the center-of-mass kinetic energy of the soliton is above the energy needed for complete dissociation of the soliton into constituent atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulating complex networks in phase space: Gaussian boson sampling

We show how phase-space simulations of Gaussian quantum states in a photonic network permit verification of measurable correlations of Gaussian boson sampling (GBS) quantum computers. Our results agree with experiments for up to 100-th order correlations, provided decoherence is included. We extend this to more than 16,000 modes, and describe how to simulate genuine multipartite entanglement.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Simulating macroscopic quantum correlations in linear networks

Many developing quantum technologies make use of quantum networks of different types. Even linear quantum networks are nontrivial, as the output photon distributions can be exponentially complex. Despite this, they can still be computationally simulated. The methods used are transformations into equivalent phase-space representations, which can then be treated probabilistically. This provides an exceptionally useful tool for the prediction and validation of experimental results, including decoherence. As well as experiments in Gaussian boson sampling, which are intended to demonstrate quantum computational advantage, these methods are applicable to other types of entangled linear quantum networks as well. This paper provides a tutorial and review of work in this area, to explain quantum phase-space techniques using the positive-P and Wigner distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Initial states and apodisation for quantum field simulations in phase-space

Bosonic quantum fields can be simulated with `quantum software' in phase-space. The positive-P, Wigner and Q-function phase-space methods are reviewed. Initial quantum states and boundaries for infinite domains are considered in detail. The quantum initial conditions treated include both multi-mode Gaussian states and number states, which are sampled using phase-shifted weighting techniques. This is more efficient than direct probabilistic sampling. Algorithms for treating periodic boundaries within infinite domains are developed via apodisation, or absorbers near the boundary. Similar truncation or anti-aliasing issues arise in the frequency domain with Fourier transform methods. Complex apodisation, in which there are phase-shifts as well as absorption, improves accuracy by reducing unwanted errors, and this is illustrated numerically. The complex power law apodisation techniques described here can be applied to numerical calculations or even directly in optics. Phase-space simulations, used for both quantum optical and Bose-Einstein condensate quantum dynamical calculations, require quantum noise terms to be included. A method is described to analyze number conservation in apodised calculations. Some techniques developed here are potentially applicable to all digital simulations, whether using classical or quantum logic.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Phase-space representations of thermal Bose-Einstein condensates

Phase-space methods allow one to go beyond the mean-field approximation to simulate the quantum dynamics of interacting fields. Here, we obtain a technique for initializing either Wigner or positive-P phase-space simulations of Bose-Einstein condensates with quantum states at a finite temperature. As a means to calculate the initial states, we introduce the idea of a nonlinear chemical potential, which removes the zero-momentum phase-noise divergences of Bogoliubov theory to give a diagonal Hamiltonian. The resulting steady-state quantum theory is then directly applicable to the calculations of initial conditions for quantum simulations of BEC dynamics using phase-space techniques. These methods allow efficient and scalable simulation of large Bose-Einstein condensates. We suggest that nonlinear chemical potentials may have a general applicability to cases of broken symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comment on: `Single-shot simulations of dynamic quantum many-body systems' [arXiv:1501.03224]

In their recent paper [Nature Physics 15, 451 (2006)], Sakmann and Kasevich study the formation of fringe patterns in ultra-cold Bose gases and claim: `Here, we show how single shots can be simulated from numerical solutions of the time-dependent many-body Schrödinger equation.' It would be remarkable if they had solved this exponentially complex equation. Instead they solve nonlinear equations with the aim to approximate the solution of the Schrödinger equation. The authors proceed to criticize phase-space approaches to simulating quantum dynamics and claim the impossibility of interpreting single trajectories of the truncated Wigner (tW) method as single-shot experimental outcomes. Here we aim to provide relevant context and elaborate why we disagree with the authors' claims.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Critical fluctuations in an optical parametric oscillator: when light behaves like magnetism

We study the nondegenerate optical parametric oscillator in a planar interferometer near threshold, where critical phenomena are expected. These phenomena are associated with nonequilibrium quantum dynamics that are known to lead to quadrature entanglement and squeezing in the oscillator field modes. We obtain a universal form for the equation describing this system, which allows a comparison with other phase transitions. We find that the unsqueezed quadratures of this system correspond to a two-dimensional XY-type model with a tricritical Lifshitz point. This leaves open the possibility of a controlled experimental investigation into this unusual class of statistical models. We evaluate the correlations of the unsqueezed quadrature using both an exact numerical simulation and a Gaussian approximation, and obtain an accurate numerical calculation of the non-Gaussian correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Parallel optimized sampling for stochastic equations

Stochastic equations play an important role in computational science, due to their ability to treat a wide variety of complex statistical problems. However, current algorithms are strongly limited by their sampling variance, which scales proportionate to 1/N_S for N_S samples. In this paper, we obtain a new class of variance reduction methods for treating stochastic equations, called parallel optimized sampling. The objective of parallel optimized sampling is to reduce the sampling variance in the observables of an ensemble of stochastic trajectories. This is achieved through calculating a finite set of observables - typically statistical moments - in parallel, and minimizing the errors compared to known values. The algorithm is both numerically efficient and unbiased. Importantly, it does not increase the errors in higher order moments, and generally reduces such errors as well. The same procedure is applied both to initial ensembles and to changes in a finite time-step. Results of these methods show that errors in initially optimized moments can be reduced to the machine precision level, typically around 10^(-16) in current hardware. For nonlinear stochastic equations, sampled moment errors during time-evolution are larger than this, due to error propagation effects. Even so, we provide evidence for error reductions of up to two orders of magnitude in a nonlinear equation example, for low order moments, which is a large practical benefit. The sampling variance typically scales as 1/N_S, but with the advantage of a very much smaller prefactor than for standard, non-optimized methods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The scaling of boson sampling experiments

Boson sampling is the problem of generating a quantum bit stream whose average is the permanent of a $n\times n$ matrix. The bitstream is created as the output of a prototype quantum computing device with $n$ input photons. It is a fundamental challenge to verify boson sampling, and the question of how output count rates scale with matrix size $n$ is crucial. Here we apply results from random matrix theory to establish scaling laws for average count rates in boson sampling experiments with arbitrary inputs and losses. The results show that, even with losses included, verification of nonclassical behaviour at large $n$ values is indeed possible.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probabilistic quantum phase-space simulation of Bell violations and their dynamical evolution

Quantum simulations of Bell inequality violations are numerically obtained using probabilistic phase space methods, namely the positive P-representation. In this approach the moments of quantum observables are evaluated as moments of variables that have values outside the normal eigenvalue range. There is thus a parallel with quantum weak measurements and weak values. Nevertheless, the representation is exactly equivalent to quantum mechanics. A number of states violating Bell inequalities are sampled, demonstrating that these quantum paradoxes can be treated with probabilistic methods. We treat quantum dynamics by simulating the time evolution of the Bell state formed via parametric down-conversion, and discuss multi-mode generalizations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probabilistic simulation of mesoscopic "Schrödinger cat" states

We carry out probabilistic phase-space sampling of mesoscopic Schrödinger cat quantum states, demonstrating multipartite Bell violations for up to 60 qubits. We use states similar to those generated in photonic and ion-trap experiments. These results show that mesoscopic quantum superpositions are directly accessible to probabilistic sampling, and we analyze the properties of sampling errors. We also demonstrate dynamical simulation of super-decoherence in ion traps. Our computer simulations can be either exponentially faster or slower than experiment, depending on the correlations measured.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum probabilistic sampling of multipartite 60-qubit Bell inequality violations

We show that violation of genuine multipartite Bell inequalities can be obtained with sampled, probabilistic phase space methods. These genuine Bell violations cannot be replicated if any part of the system is described by a local hidden variable theory. The Bell violations are simulated probabilistically using quantum phase-space representations. We treat mesoscopically large Greenberger-Horne-Zeilinger (GHZ) states having up to 60 qubits, using both a multipartite SU(2) Q-representation and the positive P-representation. Surprisingly, we find that sampling with phase-space distributions can be exponentially faster than experiment. This is due to the classical parallelism inherent in the simulation of quantum measurements using phase-space methods. Our probabilistic sampling method predicts a contradiction with local realism of "Schrödinger-cat" states that can be realized as a GHZ spin state, either in ion traps or with photonic qubits. We also present a quantum simulation of the observed super-decoherence of the ion-trap "cat" state, using a phenomenological noise model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Functional Wigner representation of BEC quantum dynamics

We develop a method of simulating the full quantum field dynamics of multi-mode multi-component Bose-Einstein condensates in a trap. We use the truncated Wigner representation to obtain a probabilistic theory that can be sampled. This method produces c-number stochastic equations which may be solved using conventional stochastic methods. The technique is valid for large mode occupation numbers. We give a detailed derivation of methods of functional Wigner representation appropriate for quantum fields. Our approach describes spatial evolution of spinor components and properly accounts for nonlinear losses. Such techniques are applicable to calculating the leading quantum corrections, including effects like quantum squeezing, entanglement, EPR correlations and interactions with engineered nonlinear reservoirs. By using a consistent expansion in the inverse density, we are able to explain an inconsistency in the nonlinear loss equations found by earlier authors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum simulations of the early universe

A procedure is described whereby a linearly coupled spinor Bose condensate can be used as a physically accessible quantum simulator of the early universe. In particular, an experiment to generate an analog of an unstable vacuum in a relativistic scalar field theory is proposed. This is related to quantum theories of the inflationary phase of the early universe. There is an unstable vacuum sector whose dynamics correspond to the quantum sine-Gordon equations in one, two or three space dimensions. Numerical simulations of the expected behavior are reported using a truncated Wigner phase-space method, giving evidence for the dynamical formation of complex spatial clusters. Preliminary results showing the dependence on coupling strength, condensate size and dimensionality are obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating Bell violations without quantum computers

We demonstrate that it is possible to simulate Bell violations using probabilistic methods. A quantum state corresponding to optical experiments that violate the Bell inequality is generated, demonstrating that these quantum paradoxes can indeed be simulated probabilistically. This provides an explicit counter-example to Feynman's claim that such classical simulations could not be carried out.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Half-quantum vortex state in a spin-orbit coupled Bose-Einstein condensate

We investigate theoretically the condensate state and collective excitations of a two-component Bose gas in two-dimensional harmonic traps subject to isotropic Rashba spin-orbit coupling. In the weakly interacting regime when the inter-species interaction is larger than the intra-species interaction ($g_{\uparrow\downarrow}>g$), we find that the condensate ground state has a half-quantum-angular-momentum vortex configuration with spatial rotational symmetry and skyrmion-type spin texture. Upon increasing the interatomic interaction beyond a threshold $g_{c}$, the ground state starts to involve higher-order angular momentum components and thus breaks the rotational symmetry. In the case of $g_{\uparrow\downarrow}<g$, the condensate becomes unstable towards the superposition of two degenerate half-quantum vortex states. Both instabilities (at $g>g_{c}$ and $g_{\uparrow\downarrow}<g$) can be determined by solving the Bogoliubov equations for collective density oscillations of the half-quantum vortex state, and by analyzing the softening of mode frequencies. We present the phase diagram as functions of the interatomic interactions and the spin-orbit coupling. In addition, we directly simulate the time-dependent Gross-Pitaevskii equation to examine the dynamical properties of the system. Finally, we investigate the stability of the half-quantum vortex state against both the trap anisotropy and anisotropy in the spin-orbit coupling term.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Confinement induced resonances in anharmonic waveguides

We develop the theory of anharmonic confinement-induced resonances (ACIR). These are caused by anharmonic excitation of the transverse motion of the center of mass (COM) of two bound atoms in a waveguide. As the transverse confinement becomes anisotropic, we find that the COM resonant solutions split for a quasi-1D system, in agreement with recent experiments. This is not found in harmonic confinement theories. A new resonance appears for repulsive couplings ($a_{3D}>0$) for a quasi-2D system, which is also not seen with harmonic confinement. After inclusion of anharmonic energy corrections within perturbation theory, we find that these ACIR resonances agree extremely well with anomalous 1D and 2D confinement induced resonance positions observed in recent experiments. Multiple even and odd order transverse ACIR resonances are identified in experimental data, including up to N=4 transverse COM quantum numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Universal contact of strongly interacting fermions at finite temperatures

The recently discovered universal thermodynamic behaviour of dilute, strongly interacting Fermi gases also implies a universal structure in the many-body pair-correlation function at short distances, as quantified by the contact ${\cal I}$. This quantity is an excellent indicator of the presence of strong correlations in these systems, which provide a highly accessible physical model for other strongly correlated quantum fluids. Here we theoretically calculate the temperature dependence of this universal contact for a Fermi gas in free space and in a harmonic trap. At high temperatures above the Fermi degeneracy temperature, $T\gtrsim T_{F}$, we obtain a reliable non-perturbative quantum virial expansion up to third order. At low temperatures we compare different approximate strong coupling theories. These make different predictions, which need to be tested either by future experiments or advanced quantum Monte Carlo simulations. We conjecture that in the universal unitarity limit, the contact or correlation decreases monotonically with increasing temperature, unless the temperature is significantly lower than the critical temperature, $T\ll T_{c}\sim0.2T_{F}$. We also discuss briefly how to measure the universal contact either in homogeneous or harmonically trapped Fermi gases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Confinement-induced resonance in quasi-one-dimensional systems under transversely anisotropic confinement

We theoretically investigate the confinement-induced resonance for quasi-one-dimensional quan- tum systems under transversely anisotropic confinement, using a two-body s-wave scattering model in the zero-energy collision limit. We predict a single resonance for any transverse anisotropy, whose position shows a slight downshift with increasing anisotropy. We compare our prediction with the recent experimental result by Haller et al. [Phys. Rev. Lett. 104, 153203 (2010)], in which two resonances are observed in the presence of transverse anisotropy. The discrepancy between theory and experiment remains to be resolved.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamic response of strongly correlated Fermi gases in the quantum virial expansion

By developing a quantum virial expansion theory, we quantitatively calculate the dynamic density response function of a trapped strongly interacting Fermi gas at high temperatures near unitarity. A clear transition from atomic to molecular responses is identified in the spectra when crossing from the BCS to BEC regimes, in qualitative agreement with recent Bragg spectroscopy observations. Our virial expansion method provides a promising way to solve the challenging strong-coupling problems and is applicable to other dynamical properties of strongly correlated Fermi gases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exact few-body results for strongly correlated quantum gases in two dimensions

The study of strongly correlated quantum gases in two dimensions has important ramifications for understanding many intriguing pheomena in solid materials, such as high-$T_{c}$ superconductivity and the fractional quantum Hall effect. However, theoretical methods are plagued by the existence of significant quantum fluctuations. Here, we present two- and three-body exact solutions for both fermions and bosons trapped in a two-dimensional harmonic potential, with an arbitrary $s$-wave scattering length. These few-particle solutions link in a natural way to the high-temperature properties of many-particle systems via a quantum virial expansion. As a concrete example, using the energy spectrum of few fermions, we calculate the second and third virial coefficients of a strongly interacting Fermi gas in two dimensions, and consequently investigate its high-temperature thermodynamics. Our thermodynamic results may be useful for ongoing experiments on two-dimensional Fermi gases. These exact results also provide an unbiased benchmark for quantum Monte Carlo simulations of two-dimensional Fermi gases at high temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mean-field study of itinerant ferromagnetism in trapped ultracold Fermi gases: Beyond the local density approximation

We theoretically investigate the itinerant ferromagnetic transition of a spherically trapped ultracold Fermi gas with spin imbalance under strongly repulsive interatomic interactions. Our study is based on a self-consistent solution of the Hartree-Fock mean-field equations beyond the widely used local density approximation. We demonstrate that, while the local density approximation holds in the paramagnetic phase, after the ferromagnetic transition it leads to a quantitative discrepancy in various thermodynamic quantities even with large atom numbers. We determine the position of the phase transition by monitoring the shape change of the free energy curve with increasing the polarization at various interaction strengths.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pseudo-gap pairing in ultracold Fermi atoms

The crossover from a BEC (Bose-Einstein condensation) to a BCS (Bardeen-Cooper-Schrieffer) superfluid in dilute gases of ultracold Fermi atoms creates an ideal environment to enrich our knowledge of strongly correlated many-body systems. These experiments are relevant to a wide range of fields from condensed matter to astrophysics. The nature of pairing in strongly interacting Fermi gases can be readily studied, thus aiding our understanding of related problems in high-T_{c} superconductors, whose mechanism is still under debate. These are not well-understood due to the large interaction parameter. Here, we calculate the dynamical properties of a normal, trapped, and strongly correlated Fermi gas, by developing a quantum cluster expansion. In ultra-cold atomic physics one can measure the elementary excitations, using rf or Bragg spectroscopy. Our calculations for the single-particle spectral function agree with the recent measurements, and clearly demonstrate pseudogap pairing in the strongly interacting regime.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Static structure factor of a strongly correlated Fermi gas at large momenta

We theoretically investigate the static structure factor of an interacting Fermi gas near the BEC-BCS crossover at large momenta. Due to short-range two-body interactions, we predict that the structure factor of unlike spin correlations $S_{\uparrow\downarrow}(q)$ falls off as $1/q$ in a universal scaling region with large momentum $\hbar q$ and large scattering length. The scaling coefficient is determined by the celebrated Tan's contact parameter, which links the short-range behavior of many-body systems to their universal thermodynamic properties. By implementing this new Tan relation together with the random-phase approximation and the virial expansion theory in various limiting cases, we show how to calculate $S_{\uparrow\downarrow}(q)$ at zero and finite temperatures for arbitrary interaction strengths, at momentum transfer higher than the Fermi momentum. Our results provide a way to experimentally confirm a new Tan relation and to accurately measure the value of contact parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Three attractively interacting fermions in a harmonic trap: Exact solution, ferromagnetism, and high-temperature thermodynamics

Three fermions with strongly repulsive interactions in a spherical harmonic trap, constitute the simplest nontrivial system that can exhibit the onset of itinerant ferromagnetism. Here, we present exact solutions for three trapped, attractively interacting fermions near a Feshbach resonance. We analyze energy levels on the upper branch of the resonance where the atomic interaction is effectively repulsive. When the s-wave scattering length a is sufficiently positive, three fully polarized fermions are energetically stable against a single spin-flip, indicating the possibility of itinerant ferromagnetism, as inferred in the recent experiment. We also investigate the high-temperature thermodynamics of a strongly repulsive or attractive Fermi gas using a quantum virial expansion. The second and third virial coefficients are calculated. The resulting equations of state can be tested in future quantitative experimental measurements at high temperatures and can provide a useful benchmark for quantum Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Universal thermodynamics of a strongly interacting Fermi gas: theory versus experiment

Strongly interacting, dilute Fermi gases exhibit a scale-invariant, universal thermodynamic behaviour. This is notoriously difficult to understand theoretically because of the absence of a small interaction parameter. Here we present a systematic comparison of theoretical predictions from different quantum many-body theories with recent experimental data of Nascimbene et. al. (arXiv:0911.0747v1). Our comparisons have no adjustable parameters, either theoretically or experimentally. A simple Gaussian pair fluctuation theory is shown to give the best quantitative agreement in the superfluid state below threshold. In the normal state, we also calculate the equation of state by using a quantum cluster expansion theory and explore in detail its applicability to low temperatures. Using the accurate experimental result for the thermodynamic function $S(T)$, we determine the temperature $T$ of a trapped Fermi gas at unitarity as a function of a non-interacting temperature $T_{i}$ which can be obtained by an adiabatic sweep to the free gas limit. By analyzing the recent experimental data, we find a normal-superfluid transition temperature $(T/T_{F})_{c}=0.19\pm0.02$ or $(T_{i}/T_{F})_{c}=0.16\pm0.02$ in a harmonic trap, where $T_{F}$ is the Fermi temperature for a trapped ideal, non-interacting Fermi gas.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Fulde-Ferrell-Larkin-Ovchinnikov states in one-dimensional spin-polarized ultracold atomic Fermi gases

We present a systematic study of quantum phases in a one-dimensional spin-polarized Fermi gas. Three comparative theoretical methods are used to explore the phase diagram at zero temperature: the mean-field theory with either an order parameter in a single-plane-wave form or a self-consistently determined order parameter using the Bogoliubov-de Gennes equations, as well as the exact soluble Bethe ansatz method. We find that a spatially inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov phase, which lies between the fully paired BCS state and the fully polarized normal state, dominates most of the phase diagram of a uniform gas. The phase transition from the BCS state to the Fulde-Ferrell-Larkin-Ovchinnikov phase is of second order, and therefore there are no phase separation states in one-dimensional homogeneous polarized gases. This is in sharp contrast to the three-dimensional situation, where a phase separation regime is predicted to occupy a very large space in the phase diagram. We conjecture that the prediction of the dominance of the phase separation phases in three dimension could be an artifact of the non-self-consistent mean-field approximation, which is heavily used in the study of three-dimensional polarized Fermi gases. We consider also the effect of a harmonic trapping potential on the phase diagram, and find that in this case the trap generally leads to phase separation, in accord with the experimental observations for a trapped gas in three dimension. We finally investigate the local fermionic density of states of the Fulde-Ferrell-Larkin-Ovchinnikov ansatz. A two-energy-gap structure is shown up, which could be used as an experimental probe of the Fulde-Ferrell-Larkin-Ovchinnikov states.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Universal thermodynamics of strongly interacting Fermi gases

Strongly interacting Fermi gases are of great current interest. Not only are fermions the most common particles in the universe, but they are also thought to have a universal thermodynamic behavior for strong interactions \cite{heiselberg,carlson,ho}. Recent experiments on ultra-cold Fermi gases provide an unprecedented opportunity to test universality in the laboratory \cite{hara,thomas,partridge,stewart,luo}. In principle this allows - for example - the interior properties of hot, dense neutron stars to be investigated on earth. Here we carry out a detailed test of this prediction. We analyze results from three ultra-cold fermion experiments involving two completely distinct atomic species in different kinds of atomic trap environments \cite{partridge,stewart,luo}. The data is compared with the predictions of a recent strong interaction theory \cite{hldepl,hldpra}. Excellent agreement is obtained, with no adjustable parameters. By extrapolating to zero temperature, we show that the experimental measurements yield a many-body parameter} $\mathbf{β\simeq-0.59\pm.07}$, \textbf{describing the universal energy of strongly interacting Fermi gases.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Monte Carlo techniques for real-time quantum dynamics

The stochastic-gauge representation is a method of mapping the equation of motion for the quantum mechanical density operator onto a set of equivalent stochastic differential equations. One of the stochastic variables is termed the "weight", and its magnitude is related to the importance of the stochastic trajectory. We investigate the use of Monte Carlo algorithms to improve the sampling of the weighted trajectories and thus reduce sampling error in a simulation of quantum dynamics. The method can be applied to calculations in real time, as well as imaginary time for which Monte Carlo algorithms are more-commonly used. The method is applicable when the weight is guaranteed to be real, and we demonstrate how to ensure this is the case. Examples are given for the anharmonic oscillator, where large improvements over stochastic sampling are observed.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The time-reversal test for stochastic quantum dynamics

The calculation of quantum dynamics is currently a central issue in theoretical physics, with diverse applications ranging from ultra-cold atomic Bose-Einstein condensates (BEC) to condensed matter, biology, and even astrophysics. Here we demonstrate a conceptually simple method of determining the regime of validity of stochastic simulations of unitary quantum dynamics by employing a time-reversal test. We apply this test to a simulation of the evolution of a quantum anharmonic oscillator with up to $6.022\times10^{23}$ (Avogadro's number) of particles. This system is realisable as a Bose-Einstein condensate in an optical lattice, for which the time-reversal procedure could be implemented experimentally.

0 citations0 reviewsNo code linkedOpen access