The intriguing phenomenon of many-body localization (MBL) has attracted significant interest recently, but a complete characterization is still lacking. In this work, we introduce the total correlations, a concept from quantum information theory capt
uring multi-partite correlations, to the study of this phenomenon. We demonstrate that the total correlations of the diagonal ensemble provides a meaningful diagnostic tool to pin-down, probe, and better understand the MBL transition and ergodicity breaking in quantum systems. In particular, we show that the total correlations has sub-linear dependence on the system size in delocalized, ergodic phases, whereas we find that it scales extensively in the localized phase developing a pronounced peak at the transition. We exemplify the power of our approach by means of an exact diagonalization study of a Heisenberg spin chain in a disordered field.
We propose and analyze a scheme for parametrically cooling bilayer cuprates based on the selective driving of a $c$-axis vibrational mode. The scheme exploits the vibration as a transducer making the Josephson plasma frequencies time-dependent. We sh
ow how modulation at the difference frequency between the intra- and interbilayer plasmon substantially suppresses interbilayer phase fluctuations, responsible for switching $c$-axis transport from a superconducting to resistive state. Our calculations indicate that this may provide a viable mechanism for stabilizing non-equilibrium superconductivity even above $T_c$, provided a finite pair density survives between the bilayers out of equilibrium.
In this work we analyze the simultaneous emergence of diffusive energy transport and local thermalization in a nonequilibrium one-dimensional quantum system, as a result of integrability breaking. Specifically, we discuss the local properties of the
steady state induced by thermal boundary driving in a XXZ spin chain with staggered magnetic field. By means of efficient large-scale matrix product simulations of the equation of motion of the system, we calculate its steady state in the long-time limit.We start by discussing the energy transport supported by the system, finding it to be ballistic in the integrable limit and diffusive when the staggered field is finite. Subsequently, we examine the reduced density operators of neighboring sites and find that for large systems they are well approximated by local thermal states of the underlying Hamiltonian in the nonintegrable regime, even for weak staggered fields. In the integrable limit, on the other hand, this behavior is lost, and the identification of local temperatures is no longer possible. Our results agree with the intuitive connection between energy diffusion and thermalization.
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 sam
pling, 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.
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.
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 s
imulation 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}$.
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 c
onfined 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.
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.
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 probab
ility 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.
The competition between electron localization and de-localization in Mott insulators underpins the physics of strongly-correlated electron systems. Photo-excitation, which re-distributes charge between sites, can control this many-body process on the
ultrafast timescale. To date, time-resolved studies have been performed in solids in which other degrees of freedom, such as lattice, spin, or orbital excitations come into play. However, the underlying quantum dynamics of bare electronic excitations has remained out of reach. Quantum many-body dynamics have only been detected in the controlled environment of optical lattices where the dynamics are slower and lattice excitations are absent. By using nearly-single-cycle near-IR pulses, we have measured coherent electronic excitations in the organic salt ET-F2TCNQ, a prototypical one-dimensional Mott Insulator. After photo-excitation, a new resonance appears on the low-energy side of the Mott gap, which oscillates at 25 THz. Time-dependent simulations of the Mott-Hubbard Hamiltonian reproduce the oscillations, showing that electronic delocalization occurs through quantum interference between bound and ionized holon-doublon pairs.
