For many optimization algorithms the time-to-solution depends not only on the problem size but also on the specific problem instance and may vary by many orders of magnitude. It is then necessary to investigate the full distribution and especially its tail. Here we analyze the distributions of annealing times for simulated annealing and simulated quantum annealing (by path integral quantum Monte Carlo) for random Ising spin glass instances. We find power-law distributions with very heavy tails, corresponding to extremely hard instances, but far broader distributions - and thus worse performance for hard instances - for simulated quantum annealing than for simulated annealing. Fast, non-adiabatic, annealing schedules can improve the performance of simulated quantum annealing for very hard instances by many orders of magnitude.
We study the out-of-equilibrium dynamics of bosonic atoms in a 1D optical lattice, after the ground-state is excited by a single spontaneous emission event, i.e. after an absorption and re-emission of a lattice photon. This is an important fundamental source of decoherence for current experiments, and understanding the resulting dynamics and changes in the many-body state is important for controlling heating in quantum simulators. Previously it was found that in the superfluid regime, simple observables relax to values that can be described by a thermal distribution on experimental time-scales, and that this breaks down for strong interactions (in the Mott insulator regime). Here we expand on this result, investigating the relaxation of the momentum distribution as a function of time, and discussing the relationship to eigenstate thermalization. For the strongly interacting limit, we provide an analytical analysis for the behavior of the system, based on an effective low-energy Hamiltonian in which the dynamics can be understood based on correlated doublon-holon pairs.
We explore two complementary modifications of the hybridization-expansion continuous-time Monte Carlo method, aiming at large multi-orbital quantum impurity problems. One idea is to compute the imaginary-time propagation using a matrix product states representation. We show that bond dimensions considerably smaller than the dimension of the Hilbert space are sufficient to obtain accurate results, and that this approach scales polynomially, rather than exponentially with the number of orbitals. Based on scaling analyses, we conclude that a matrix product state implementation will outperform the exact-diagonalization based method for quantum impurity problems with more than 12 orbitals. The second idea is an improved Monte Carlo sampling scheme which is applicable to all variants of the hybridization expansion method. We show that this so-called sliding window sampling scheme speeds up the simulation by at least an order of magnitude for a broad range of model parameters, with the largest improvements at low temperature.
Variational studies of the t-J model on the square lattice based on infinite projected-entangled pair states (iPEPS) confirm an extremely close competition between a uniform d-wave superconducting state and different stripe states. The site-centered stripe with an in-phase d-wave order has an equal or only slightly lower energy than the stripe with anti-phase d-wave order. The optimal stripe filling is not constant but increases with J/t. A nematic anisotropy reduces the pairing amplitude and the energies of stripe phases are lowered relative to the uniform state with increasing nematicity.
We study the thermalization of excitations generated by spontaneous emission events for cold bosons in an optical lattice. Computing the dynamics described by the many-body master equation, we characterize equilibration timescales in different parameter regimes. For simple observables, we find regimes in which the system relaxes rapidly to values in agreement with a thermal distribution, and others where thermalization does not occur on typical experimental timescales. Because spontaneous emissions lead effectively to a local quantum quench, this behavior is strongly dependent on the low-energy spectrum of the Hamiltonian, and undergoes a qualitative change at the Mott Insulator-superfluid transition point. These results have important implications for the understanding of thermalization after localized quenches in isolated quantum gases, as well as the characterization of heating in experiments.
We investigate the electronic structure of cobalt atoms on a copper surface and in a copper host by combining density functional calculations with a numerically exact continuous-time quantum Monte Carlo treatment of the five-orbital impurity problem. In both cases we find low energy resonances in the density of states of all five Co $d$-orbitals. The corresponding self-energies indicate the formation of a Fermi liquid state at low temperatures. Our calculations yield the characteristic energy scale -- the Kondo temperature -- for both systems in good agreement with experiments. We quantify the charge fluctuations in both geometries and suggest that Co in Cu must be described by an Anderson impurity model rather than by a model assuming frozen impurity valency at low energies. We show that fluctuations of the orbital degrees of freedom are crucial for explaining the Kondo temperatures obtained in our calculations and measured in experiments.
In this chapter we discuss aspects of the quantum critical behavior that occurs at a quantum phase transition separating a topological phase from a conventionally ordered one. We concentrate on a family of quantum lattice models, namely certain deformations of the toric code model, that exhibit continuous quantum phase transitions. One such deformation leads to a Lorentz-invariant transition in the 3D Ising universality class. An alternative deformation gives rise to a so-called conformal quantum critical point where equal-time correlations become conformally invariant and can be related to those of the 2D Ising model. We study the behavior of several physical observables, such as non-local operators and entanglement entropies, that can be used to characterize these quantum phase transitions. Finally, we briefly consider the role of thermal fluctuations and related phase transitions, before closing with a short overview of field theoretical descriptions of these quantum critical points.
Using the Ehrenfest urn model we illustrate the subtleties of error estimation in Monte Carlo simulations. We discuss how the smooth results of correlated sampling in Markov chains can fool ones perception of the accuracy of the data, and show (via numerical and analytical methods) how to obtain reliable error estimates from correlated samples.
We discuss how to construct models of interacting anyons by generalizing quantum spin Hamiltonians to anyonic degrees of freedom. The simplest interactions energetically favor pairs of anyons to fuse into the trivial (identity) channel, similar to the quantum Heisenberg model favoring pairs of spins to form spin singlets. We present an introduction to the theory of anyons and discuss in detail how basis sets and matrix representations of the interaction terms can be obtained, using non-Abelian Fibonacci anyons as example. Besides discussing the golden chain, a one-dimensional system of anyons with nearest neighbor interactions, we also present the derivation of more complicated interaction terms, such as three-anyon interactions in the spirit of the Majumdar-Ghosh spin chain, longer range interactions and two-leg ladders. We also discuss generalizations to anyons with general non-Abelian su(2)_k statistics. The k to infinity limit of the latter yields ordinary SU(2) spin chains.
Two-dimensional quantum loop gases are elementary examples of topological ground states with Abelian or non-Abelian anyonic excitations. While Abelian loop gases appear as ground states of local, gapped Hamiltonians such as the toric code, we show that gapped non-Abelian loop gases require non-local interactions (or non-trivial inner products). Perturbing a local, gapless Hamiltonian with an anticipated ``non-Abelian ground-state wavefunction immediately drives the system into the Abelian phase, as can be seen by measuring the Hausdorff dimension of loops. Local quantum critical behavior is found in a loop gas in which all equal-time correlations of local operators decay exponentially.
