Non-destructive selective probing of phononic excitations in a cold Bose gas using impurities

We introduce a detector that selectively probes the phononic excitations of a cold Bose gas. The detector is composed of a single impurity atom confined by a double-well potential, where the two lowest eigenstates of the impurity form an effective probe qubit that is coupled to the phonons via density-density interactions with the bosons. The system is analogous to a two-level atom coupled to photons of the radiation field. We demonstrate that tracking the evolution of the qubit populations allows probing both thermal and coherent excitations in targeted phonon modes. The targeted modes are selected in both energy and momentum by adjusting the impuritys potential. We show how to use the detector to observe coherent density waves and to measure temperatures of the Bose gas down to the nano-Kelvin regime. We analyze how our scheme could be realized experimentally, including the possibility of using an array of multiple impurities to achieve greater precision from a single experimental run.

Capturing exponential variance using polynomial resources: applying tensor networks to non-equilibrium stochastic processes

Estimating the expected value of an observable appearing in a non-equilibrium stochastic process usually involves sampling. If the observables variance is high, many samples are required. In contrast, we show that performing the same task without sampling, using tensor network compression, efficiently captures high variances in systems of various geometries and dimensions. We provide examples for which matching the accuracy of our efficient method would require a sample size scaling exponentially with system size. In particular, the high variance observable $mathrm{e}^{-beta W}$, motivated by Jarzynskis equality, with $W$ the work done quenching from equilibrium at inverse temperature $beta$, is exactly and efficiently captured by tensor networks.

What is a quantum simulator?

Quantum simulators are devices that actively use quantum effects to answer questions about model systems and, through them, real systems. Here we expand on this definition by answering several fundamental questions about the nature and use of quantum simulators. Our answers address two important areas. First, the difference between an operation termed simulation and another termed computation. This distinction is related to the purpose of an operation, as well as our confidence in and expectation of its accuracy. Second, the threshold between quantum and classical simulations. Throughout, we provide a perspective on the achievements and directions of the field of quantum simulation.

Solving search problems by strongly simulating quantum circuits

Simulating quantum circuits using classical computers lets us analyse the inner workings of quantum algorithms. The most complete type of simulation, strong simulation, is believed to be generally inefficient. Nevertheless, several efficient strong simulation techniques are known for restricted families of quantum circuits and we develop an additional technique in this article. Further, we show that strong simulation algorithms perform another fundamental task: solving search problems. Efficient strong simulation techniques allow solutions to a class of search problems to be counted and found efficiently. This enhances the utility of strong simulation methods, known or yet to be discovered, and extends the class of search problems known to be efficiently simulable. Relating strong simulation to search problems also bounds the computational power of efficiently strongly simulable circuits; if they could solve all problems in $mathrm{P}$ this would imply the collapse of the complexity hierarchy $mathrm{P} subseteq mathrm{NP} subseteq # mathrm{P}$.

Breathing oscillations of a trapped impurity in a Bose gas

Motivated by a recent experiment [J. Catani et al., arXiv:1106.0828v1 preprint, 2011], we study breathing oscillations in the width of a harmonically trapped impurity interacting with a separately trapped Bose gas. We provide an intuitive physical picture of such dynamics at zero temperature, using a time-dependent variational approach. In the Gross-Pitaevskii regime we obtain breathing oscillations whose amplitudes are suppressed by self trapping, due to interactions with the Bose gas. Introducing phonons in the Bose gas leads to the damping of breathing oscillations and non-Markovian dynamics of the width of the impurity, the degree of which can be engineered through controllable parameters. Our results reproduce the main features of the impurity dynamics observed by Catani et al. despite experimental thermal effects, and are supported by simulations of the system in the Gross-Pitaevskii regime. Moreover, we predict novel effects at lower temperatures due to self-trapping and the inhomogeneity of the trapped Bose gas.

Impurity transport through a strongly interacting bosonic quantum gas

Using near-exact numerical simulations we study the propagation of an impurity through a one-dimensional Bose lattice gas for varying bosonic interaction strengths and filling factors at zero temperature. The impurity is coupled to the Bose gas and confined to a separate tilted lattice. The precise nature of the transport of the impurity is specific to the excitation spectrum of the Bose gas which allows one to measure properties of the Bose gas non-destructively, in principle, by observing the impurity; here we focus on the spatial and momentum distributions of the impurity as well as its reduced density matrix. For instance we show it is possible to determine whether the Bose gas is commensurately filled as well as the bandwidth and gap in its excitation spectrum. Moreover, we show that the impurity acts as a witness to the cross-over of its environment from the weakly to the strongly interacting regime, i.e., from a superfluid to a Mott insulator or Tonks-Girardeau lattice gas and the effects on the impurity in both of these strongly-interacting regimes are clearly distinguishable. Finally, we find that the spatial coherence of the impurity is related to its propagation through the Bose gas, giving an experimentally controllable example of noise-enhanced quantum transport.

Phonon resonances in atomic currents through Bose-Fermi mixtures in optical lattices

We present an analysis of Bose-Fermi mixtures in optical lattices for the case where the lattice potential of the fermions is tilted and the bosons (in the superfluid phase) are described by Bogoliubov phonons. It is shown that the Bogoliubov phonons enable hopping transitions between fermionic Wannier-Stark states; these transitions are accompanied by energy dissipation into the superfluid and result in a net atomic current along the lattice. We derive a general expression for the drift velocity of the fermions and find that the dependence of the atomic current on the lattice tilt exhibits negative differential conductance and phonon resonances. Numerical simulations of the full dynamics of the system based on the time-evolving block decimation algorithm reveal that the phonon resonances should be observable under the conditions of a realistic measuring procedure.

Dynamical simulations of classical stochastic systems using matrix product states

We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of 1D quantum systems, to simulate the time-evolution of non-equilibrium stochastic systems. We describe this method in detail; a systems probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time-evolution is efficiently simulated by repeatedly updating and approximately re-factorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible, even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size.