BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Peter D. Drummond

Peter D. Drummond contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.quant-gasquant-phcond-mat.str-elgr-qcastro-ph.COcond-mat.stat-mechcond-mat.supr-conmath-phmath.MPmath.NANumerical Analysisphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

24 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
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
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
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