Systems passing through quantum critical points at finite rates have a finite probability of undergoing transitions between different eigenstates of the instantaneous Hamiltonian. This mechanism was proposed by Kibble as the underlying mechanism for
the formation of topological defects in the early universe and by Zurek for condensed matter systems. Here, we use a system of nuclear spins as an experimental quantum simulator undergoing a non-equilibrium quantum phase transition. The experimental data confirm the validity of the Kibble-Zurek mechanism of defect formation.
The Trotter-Suzuki approximation leads to an efficient algorithm for solving the time-dependent Schrodinger equation. Using existing highly optimized CPU and GPU kernels, we developed a distributed version of the algorithm that runs efficiently on a
cluster. Our implementation also improves single node performance, and is able to use multiple GPUs within a node. The scaling is close to linear using the CPU kernels, whereas the efficiency of GPU kernels improve with larger matrices. We also introduce a hybrid kernel that simultaneously uses multicore CPUs and GPUs in a distributed system. This kernel is shown to be efficient when the matrix size would not fit in the GPU memory. Larger quantum systems scale especially well with a high number nodes. The code is available under an open source license.
Various fundamental phenomena of strongly-correlated quantum systems such as high-$T_c$ superconductivity, the fractional quantum-Hall effect, and quark confinement are still awaiting a universally accepted explanation. The main obstacle is the compu
tational complexity of solving even the most simplified theoretical models that are designed to capture the relevant quantum correlations of the many-body system of interest. In his seminal 1982 paper [Int. J. Theor. Phys. 21, 467], Richard Feynman suggested that such models might be solved by simulation with a new type of computer whose constituent parts are effectively governed by a desired quantum many-body dynamics. Measurements on this engineered machine, now known as a quantum simulator, would reveal some unknown or difficult to compute properties of a model of interest. We argue that a useful quantum simulator must satisfy four conditions: relevance, controllability, reliability, and efficiency. We review the current state of the art of digital and analog quantum simulators. Whereas so far the majority of the focus, both theoretically and experimentally, has been on controllability of relevant models, we emphasize here the need for a careful analysis of reliability and efficiency in the presence of imperfections. We discuss how disorder and noise can impact these conditions, and illustrate our concerns with novel numerical simulations of a paradigmatic example: a disordered quantum spin chain governed by the Ising model in a transverse magnetic field. We find that disorder can decrease the reliability of an analog quantum simulator of this model, although large errors in local observables are introduced only for strong levels of disorder. We conclude that the answer to the question Can we trust quantum simulators? is... to some extent.
Systems with long-range interactions show a variety of intriguing properties: they typically accommodate many meta-stable states, they can give rise to spontaneous formation of supersolids, and they can lead to counterintuitive thermodynamic behavior
. However, the increased complexity that comes with long-range interactions strongly hinders theoretical studies. This makes a quantum simulator for long-range models highly desirable. Here, we show that a chain of trapped ions can be used to quantum simulate a one-dimensional model of hard-core bosons with dipolar off-site interaction and tunneling, equivalent to a dipolar XXZ spin-1/2 chain. We explore the rich phase diagram of this model in detail, employing perturbative mean-field theory, exact diagonalization, and quasiexact numerical techniques (density-matrix renormalization group and infinite time evolving block decimation). We find that the complete devils staircase -- an infinite sequence of crystal states existing at vanishing tunneling -- spreads to a succession of lobes similar to the Mott-lobes found in Bose--Hubbard models. Investigating the melting of these crystal states at increased tunneling, we do not find (contrary to similar two-dimensional models) clear indications of supersolid behavior in the region around the melting transition. However, we find that inside the insulating lobes there are quasi-long range (algebraic) correlations, opposed to models with nearest-neighbor tunneling which show exponential decay of correlations.
We use NMR quantum simulators to study antiferromagnetic Ising spin chains undergoing quantum phase transitions. Taking advantage of the sensitivity of the systems near criticality, we detect the critical points of the transitions using a direct meas
urement of the Loschmidt echo. We test our simulators for spin chains of even and odd numbers of spins, and compare the experimental results to theoretical predictions.
I show how to perform a Loschmidt echo (time reversal) in the Bose-Hubbard model implemented with cold bosonic atoms in an optical lattice. The echo is obtained by applying a linear phase imprint on the lattice and a change in magnetic field to tune
the boson-boson scattering length through a Feshbach resonance. I discuss how the echo can measure the fidelity of the quantum simulation, the intensity of an external potential (e.g. gravity), or the critical point of the superfluid-insulator quantum phase transition.
We study the decay rate of the Loschmidt echo or fidelity in a chaotic system under a time-dependent perturbation $V(q,t)$ with typical strength $hbar/tau_{V}$. The perturbation represents the action of an uncontrolled environment interacting with th
e system, and is characterized by a correlation length $xi_0$ and a correlation time $tau_0$. For small perturbation strengths or rapid fluctuating perturbations, the Loschmidt echo decays exponentially with a rate predicted by the Fermi Golden Rule, $1/tilde{tau}= tau_{c}/tau_{V}^2$, where typically $tau_{c} sim min[tau_{0},xi_0/v]$ with $v$ the particle velocity. Whenever the rate $1/tilde{tau}$ is larger than the Lyapunov exponent of the system, a perturbation independent Lyapunov decay regime arises. We also find that by speeding up the fluctuations (while keeping the perturbation strength fixed) the fidelity decay becomes slower, and hence, one can protect the system against decoherence.
A neutral impurity atom immersed in a dilute Bose-Einstein condensate (BEC) can have a bound ground state in which the impurity is self-localized. In this small polaron-like state, the impurity distorts the density of the surrounding BEC, thereby cre
ating the self-trapping potential minimum. We describe the self-localization in a strong coupling approach.
The Loschmidt echo (LE) is a measure of the sensitivity of quantum mechanics to perturbations in the evolution operator. It is defined as the overlap of two wave functions evolved from the same initial state but with slightly different Hamiltonians.
Thus, it also serves as a quantification of irreversibility in quantum mechanics. In this thesis the LE is studied in systems that have a classical counterpart with dynamical instability, that is, classically chaotic. An analytical treatment that makes use of the semiclassical approximation is presented. It is shown that, under certain regime of the parameters, the LE decays exponentially. Furthermore, for strong enough perturbations, the decay rate is given by the Lyapunov exponent of the classical system. Some particularly interesting examples are given. The analytical results are supported by thorough numerical studies. In addition, some regimes not accessible to the theory are explored, showing that the LE and its Lyapunov regime present the same form of universality ascribed to classical chaos. In a sense, this is evidence that the LE is a robust temporal signature of chaos in the quantum realm. Finally, the relation between the LE and the quantum to classical transition is explored, in particular with the theory of decoherence. Using two different approaches, a semiclassical approximation to Wigner functions and a master equation for the LE, it is shown that the decoherence rate and the decay rate of the LE are equal. The relationship between these quantities results mutually beneficial, in terms of the broader resources of decoherence theory and of the possible experimental realization of the LE.
