Search

We identify persistent oscillations in a nonintegrable quantum Ising chain left behind by a rapid transition into a ferromagnetic phase. In the integrable chain with nearest-neighbor (NN) interactions, the nature, origin, and decay of post-transition oscillations are tied to the Kibble-Zurek mechanism (KZM). However, when coupling to the next nearest neighbor (NNN) is added, the resulting nonintegrable Ising chain (still in the quantum Ising chain universality class) supports persistent post-transition oscillation: KZM-like oscillations turn into persistent oscillations of transverse magnetization. Their longevity in our simulations is likely limited only by the numerical accuracy. Their period differs from the decaying KZM oscillation but their amplitude depends on quench rate. Moreover, they can be excited by driving in resonance with the excitations' energy gap. Thus, while one might have expected that the integrability-breaking NNN coupling would facilitate relaxation, the oscillations we identify are persistent. At low to medium transverse fields, they are associated with Cooper pairs of Bogoliubov quasiparticles -- kinks. This oscillation of the pair condensate is a manifestation of quantum coherence., Comment: 6 pages, 5 figures, plus Supplementary Material

Adiabatic preparation of a critical ground state is hampered by the closing of its energy gap as the system size increases. However, this gap is directly relevant only for a uniform ramp, where a control parameter in the Hamiltonian is tuned uniformly in space towards the quantum critical point. Here, we consider inhomogeneous ramps in two dimensions: initially, the parameter is made critical at the center of a lattice, from where the critical region expands at a fixed velocity. In the 1D and 2D quantum Ising models, which have a well-defined speed of sound at the critical point, the ramp becomes adiabatic with a subsonic velocity. This subsonic ramp can prepare the critical state faster than a uniform one. Moreover, in both a model of $p$-wave paired 2D fermions and the Kitaev model, the critical dispersion is anisotropic -- linear with a nonzero velocity in one direction and quadratic in the other -- but the gap is still inversely proportional to the linear size of the critical region, with a coefficient proportional to the nonzero velocity. This suffices to make the inhomogeneous ramp adiabatic below a finite crossover velocity and superior to the homogeneous one., Comment: 12 pages, 9 figures; Updated with changes made for publication in Physical Review B

We present an open-source tensor network Python library for quantum many-body simulations. At its core is an abelian-symmetric tensor, implemented as a sparse block structure managed by logical layer on top of dense multi-dimensional array backend. This serves as the basis for higher-level tensor networks algorithms, operating on matrix product states and projected entangled pair states, implemented here. Using appropriate backend, such as PyTorch, gives direct access to automatic differentiation (AD) for cost-function gradient calculations and execution on GPUs or other supported accelerators. We show the library performance in simulations with infinite projected entangled-pair states, such as finding the ground states with AD, or simulating thermal states of the Hubbard model via imaginary time evolution. We quantify sources of performance gains in those challenging examples allowed by utilizing symmetries., Comment: 10 pages, 6 figures, source code at https://github.com/yastn/yastn

Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations demand resources that grow exponentially with system size. Here we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schr\"odinger equation. We demonstrate area-law scaling of entanglement in the model quench in two-, three- and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for classical approaches. We assess approximate methods based on tensor networks and neural networks and conclude that no known approach can achieve the same accuracy as the quantum annealer within a reasonable timeframe. Thus quantum annealers can answer questions of practical importance that classical computers cannot.

The Minimally Entangled Typical Thermal States (METTS) are an ensemble of pure states, equivalent to the Gibbs thermal state, that can be efficiently represented by tensor networks. In this article, we use the Projected Entangled Pair States (PEPS) ansatz as to represent METTS on a two-dimensional (2D) lattice. While Matrix Product States (MPS) are less efficient for 2D systems due to their complexity growing exponentially with the lattice size, PEPS provide a more tractable approach. To substantiate the prowess of PEPS in modeling METTS (dubbed as PEPS-METTS), we benchmark it against the purification method for the 2D quantum Ising model at its critical temperature. Our analysis reveals that PEPS-METTS achieves accurate long-range correlations with significantly lower bond dimensions. We further corroborate this finding in the 2D Fermi Hubbard model at half-filling. At a technical level, we introduce an efficient \textit{zipper} method to obtain PEPS boundary matrix product states needed to compute expectation values. The imaginary time evolution is performed with the neighbourhood tensor update., Comment: 14 pages, 14 figures

Resonant transport occurs when there is a matching of frequencies across some spatial medium, increasing the efficiency of shuttling particles from one reservoir to another. We demonstrate that in a periodically driven, many--body titled lattice, there are sets of spatially fractured resonances. These ``emanate'' from two essential resonances due to scattering off internal surfaces created when the driving frequency and many--body interaction strength vary, a scattering reminiscent of lens flare. The confluence of these fractured resonances dramatically enhances transport. At one confluence, the interaction strength is finite and the essential resonance arises due to the interplay of interaction with the counter--rotating terms of the periodic drive. We discuss the origin and structure of the fractured resonances, as well as the scaling of the conductance with system parameters. These results furnish a new example of the richness of open, driven, many--body systems., Comment: significantly revised version! comments welcome

Extended reservoirs provide a framework for capturing macroscopic, continuum environments, such as metallic electrodes driving a current through a nanoscale contact, impurity, or material. We examine the application of this approach to periodically driven systems, specifically in the context of quantum transport. As with non--equilibrium steady states in time--independent scenarios, the current displays a Kramers' turnover including the formation of a plateau region that captures the physical, continuum limit response. We demonstrate that a simple stability criteria identifies an appropriate relaxation rate to target this physical plateau. Using this approach, we study quantum transport through a periodically driven tilted lattice coupled to two metallic reservoirs held at a finite bias and temperature. We use this model to benchmark the extended reservoir approach and assess the stability criteria. The approach recovers well--understood physical behavior in the limit of weak system--reservoir coupling. Extended reservoirs enable addressing strong coupling and non--linear response as well, where we analyze how transport responds to the dynamics inside the driven lattice. These results set the foundations for the use of extended reservoir approach for periodically driven, quantum systems, such as many--body Floquet states., Comment: 14 pages, 8 figures

The simulation of open many-body quantum systems is challenging, requiring methods to both handle exponentially large Hilbert spaces and represent the influence of (infinite) particle and energy reservoirs. These two requirements are at odds with each other: Larger collections of modes can increase the fidelity of the reservoir representation but come at a substantial computational cost when included in numerical many-body techniques. An increasingly utilized and natural approach to control the growth of the reservoir is to cast a finite set of reservoir modes themselves as an open quantum system. There are, though, many routes to do so. Here, we introduce an accumulative reservoir construction -- an ARC -- that employs a series of partial refreshes of the extended reservoirs. Through this series, the representation accumulates the character of an infinite reservoir. This provides a unified framework for both continuous (Lindblad) relaxation and a recently introduced periodically refresh approach (i.e., discrete resets of the reservoir modes to equilibrium). In the context of quantum transport, we show that the phase space for physical behavior separates into discrete and continuous relaxation regimes with the boundary between them set by natural, physical timescales. Both of these regimes "turnover" into regions of over- and under-damped coherence in a way reminiscent of Kramers' crossover. We examine how the range of behavior impacts errors and the computational cost, including within tensor networks. These results provide the first comparison of distinct extended reservoir approaches, showing that they have different scaling of error versus cost (with a bridging ARC regime decaying fastest). Exploiting the enhanced scaling, though, will be challenging, as it comes with a substantial increase in (operator space) entanglement entropy., Comment: 27 pages, 18 figures

The Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network, an infinite projected entangled pair state, evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With U(1)xU(1) symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators, measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit, provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime., Comment: 18 pages; minor changes; accepted in PRB

Linear quench of the transverse field drives the quantum Ising chain across a quantum critical point from the paramagnetic to the ferromagnetic phase. We focus on normal and anomalous quadratic correlators between fermionic kink creation and annihilation operators. They depend not only on the Kibble-Zurek (KZ) correlation length but also on a dephasing length scale, which differs from the KZ length by a logarithmic correction. Additional slowing down of the ramp in the ferromagnetic phase further increases the dephasing length and suppresses the anomalous correlator. The quadratic correlators enter Pfaffians that yield experimentally relevant kink correlation functions, the probability distribution of ferromagnetic domain sizes, and, closely related, emptiness formation probability. The latter takes the form of a Pfaffian of a block Toeplitz matrix that allows for some analytic asymptotes. Finally, we obtain further insight into the structure of the state at the end of the ramp by interpreting it as a paired state of fermionic kinks characterized by its pair wave function. All those quantities are sensitive to quantum coherence between eigenstates with different numbers of kinks, thus making them a convenient probe of the quantumness of a quantum simulator platform., Comment: 14 pages; minor improvements in presentation; to appear in Phys. Rev. B

It is now widely accepted that quenches through the critical region of quantum phase transitions result in post-transition states populated with topological defects -- analogs of the classical topological defects. However, consequences of the very non-classical fact that the state after a quench is a {\it superposition} of distinct, broken-symmetry vacua with different numbers and locations of defects have remained largely unexplored. We identify coherent quantum oscillations induced by such superpositions in observables complementary to the one involved in symmetry breaking. These oscillations satisfy Kibble-Zurek dynamical scaling laws with the quench rate, with an instantaneous oscillation frequency set primarily by the gap of the system. In addition to the obvious fundamental significance of a superposition of different broken symmetry states, quantum coherent oscillations can be used to verify unitarity and test for imperfections of the experimental implementations of quantum simulators., Comment: version to appear in Phys. Rev. Lett. with a new example of non-integrable model

Sampling a diverse set of high-quality solutions for hard optimization problems is of great practical relevance in many scientific disciplines and applications, such as artificial intelligence and operations research. One of the main open problems is the lack of ergodicity, or mode collapse, for typical stochastic solvers based on Monte Carlo techniques leading to poor generalization or lack of robustness to uncertainties. Currently, there is no universal metric to quantify such performance deficiencies across various solvers. Here, we introduce a new diversity measure for quantifying the number of independent approximate solutions for NP-hard optimization problems. Among others, it allows benchmarking solver performance by a required time-to-diversity (TTD), a generalization of often used time-to-solution (TTS). We illustrate this metric by comparing the sampling power of various quantum annealing strategies. In particular, we show that the inhomogeneous quantum annealing schedules can redistribute and suppress the emergence of topological defects by controlling space-time separated critical fronts, leading to an advantage over standard quantum annealing schedules with respect to both TTS and TTD for finding rare solutions. Using path-integral Monte Carlo simulations for up to 1600 qubits, we demonstrate that nonequilibrium driving of quantum fluctuations, guided by efficient approximate tensor network contractions, can significantly reduce the fraction of hard instances for random frustrated 2D spin-glasses with local fields. Specifically, we observe that by creating a class of algorithmic quantum phase transitions, the diversity of solutions can be enhanced by up to 40% with the fraction of hard-to-sample instances reducing by more than 25%., Comment: 13 pages, 6 figures, close to published version