Source author record

Anatoli Polkovnikov

Anatoli Polkovnikov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

54works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2022arXiv

Counter-diabatic driving in the classical $β$-Fermi-Pasta-Ulam-Tsingou chain

Shortcuts to adiabaticity (STAs) have been used to make rapid changes to a system while eliminating or minimizing excitations in the system's state. In quantum systems, these shortcuts allow us to minimize inefficiencies and heating in experiments and quantum computing protocols, but the theory of STAs can also be generalized to classical systems. We focus on one such STA, approximate counter-diabatic (ACD) driving, and numerically compare its performance in two classical systems: a quartic anharmonic oscillator and the $β$ Fermi-Pasta-Ulam-Tsingou (FPUT) lattice. In particular, we modify an existing variational technique to optimize the approximate driving and then develop classical figures of merit to quantify the performance of the driving. We find that relatively simple forms for the ACD driving can dramatically suppress excitations regardless of system size. ACD driving in classical nonlinear oscillators could have many applications, from minimizing heating in bosonic gases to finding optimal local dressing protocols in interacting field theories.

preprint2021arXiv

Shortcuts to Dynamic Polarization

Dynamic polarization protocols aim to hyperpolarize a spin bath by transferring spin polarization from a well-controlled qubit such as a quantum dot or a color defect. Building on techniques from shortcuts to adiabaticity, we design fast and efficient dynamic polarization protocols in central spin models that apply to dipolarly interacting systems. The protocols maximize the transfer of polarization via bright states at a nearby integrable point, exploit the integrability-breaking terms to reduce the statistical weight on dark states that do not transfer polarization, and realize experimentally accessible local counterdiabatic driving through Floquet-engineering. A master equation treatment suggests that the protocol duration scales linearly with the number of bath spins with a prefactor that can be orders of magnitude smaller than that of unassisted protocols. This work opens new pathways to cool spin baths and extend qubit coherence times for applications in quantum information processing and metrology.

preprint2020arXiv

Adiabatic landscape and optimal paths in ergodic systems

Whether one is interested in quantum state preparation or in the design of efficient heat engines, adiabatic (reversible) transformations play a pivotal role in minimizing computational complexity and energy losses. Understanding the structure of these transformations and identifying the systems for which such transformations can be performed efficiently and quickly is therefore of primary importance. In this paper we focus on finding optimal paths in the space of couplings controlling the system's Hamiltonian. More specifically, starting from a local Hamiltonian we analyze directions in the space of couplings along which adiabatic transformations can be accurately generated by local operators, which are both realizable in experiments and easy to simulate numerically. We consider a non-integrable 1D Ising model parametrized by two independent couplings, corresponding to longitudinal and transverse magnetic fields. We find regions in the space of couplings characterized by a very strong anisotropy of the variational adiabatic gauge potential (AGP), generating the adiabatic transformations, which allows us to define optimal adiabatic paths. We find that these paths generally terminate at singular points characterized by extensive degeneracies in the energy spectrum, splitting the parameter space into adiabatically disconnected regions. The anisotropy follows from singularities in the AGP, and we identify special robust weakly-thermalizing and non-absorbing many-body "dark" states which are annihilated by the singular part of the AGP and show that their existence extends deep into the ergodic regime.

preprint2020arXiv

Emergent conservation laws and nonthermal states in the mixed-field Ising model

This paper presents a method of computing approximate conservation laws and eigenstates of integrability-broken models using the concept of adiabatic continuation. Given some Hamiltonian, eigenstates and conserved operators may be computed by using those of a simple Hamiltonian close by in parameter space, dressed by some unitary rotation. However, most adiabatic continuation analyses only use this unitary implicitly. In this work, approximate adiabatic gauge potentials are used to construct a state dressing using variational methods, to compute eigenstates via a rotated truncated spectrum approximation. These methods allow construction of both low and high-energy approximate nonthermal eigenstates, as well as quasi-local almost-conserved operators, in models where integrability may be non-perturbatively broken. These concepts will be demonstrated in the mixed-field Ising model.

preprint2020arXiv

Quantum echo dynamics in the Sherrington-Kirkpatrick model

Understanding the footprints of chaos in quantum-many-body systems has been under debate for a long time. In this work, we study the echo dynamics of the Sherrington-Kirkpatrick (SK) model with transverse field under effective time reversal. We investigate numerically its quantum and semiclassical dynamics. We explore how chaotic many-body quantum physics can lead to exponential divergence of the echo of observables and we show that it is a result of three requirements: i) the collective nature of the observable, ii) a properly chosen initial state and iii) the existence of a well-defined chaotic semi-classical (large-$N$) limit. Under these conditions, the echo grows exponentially up to the Ehrenfest time, which scales logarithmically with the number of spins $N$. In this regime, the echo is well described by the semiclassical (truncated Wigner) approximation. We also discuss a short-range version of the SK model, where the Ehrenfest time does not depend on $N$ and the quantum echo shows only polynomial growth. Our findings provide new insights on scrambling and echo dynamics and how to observe it experimentally.

preprint2020arXiv

Quantum eigenstates from classical Gibbs distributions

We discuss how the language of wave functions (state vectors) and associated non-commuting Hermitian operators naturally emerges from classical mechanics by applying the inverse Wigner-Weyl transform to the phase space probability distribution and observables. In this language, the Schrödinger equation follows from the Liouville equation, with $\hbar$ now a free parameter. Classical stationary distributions can be represented as sums over stationary states with discrete (quantized) energies, where these states directly correspond to quantum eigenstates. Interestingly, it is now classical mechanics which allows for apparent negative probabilities to occupy eigenstates, dual to the negative probabilities in Wigner's quasiprobability distribution. These negative probabilities are shown to disappear when allowing sufficient uncertainty in the classical distributions. We show that this correspondence is particularly pronounced for canonical Gibbs ensembles, where classical eigenstates satisfy an integral eigenvalue equation that reduces to the Schrödinger equation in a saddle-point approximation controlled by the inverse temperature. We illustrate this correspondence by showing that some paradigmatic examples such as tunneling, band structures, Berry phases, Landau levels, level statistics and quantum eigenstates in chaotic potentials can be reproduced to a surprising precision from a classical Gibbs ensemble, without any reference to quantum mechanics and with all parameters (including $\hbar$) on the order of unity.

preprint2020arXiv

Semiclassical dynamics of a disordered two-dimensional Hubbard model with long-range interactions

Quench dynamics in a two-dimensional system of interacting fermions is analyzed within the semiclassical truncated Wigner approximation (TWA). The models with short-range and long-range interactions are considered. We show that in the latter case, the TWA is very accurate, becoming asymptotically exact in the infinite-range limit, provided that the semiclassical Hamiltonian is correctly identified. Within the TWA, different dynamical timescales of charges and spins can be clearly distinguished. Interestingly, for a weak and moderate disorder strength, we observe subdiffusive behavior of charges, while spins exhibit diffusive dynamics. At strong disorder, the quantum Fisher information shows logarithmic growth in time with a slower increase for charges than for spins. It is shown that in contrast to the short-range model, strong inhomogeneities such as domain walls in the initial state can significantly slow down thermalization dynamics, especially at weak disorder. This behavior can put additional challenges in designing cold-atom experimental protocols aimed to analyze possible many-body localization in such systems. While within this approach we cannot make any definite statements about the existence of a many-body localized phase, we see a very fast crossover as a function of disorder strength from rapidly thermalizing to a slow glassy like regime both for the short-range and long-range models.

preprint2019arXiv

Swift heat transfer by fast-forward driving in open quantum systems

Typically, time-dependent thermodynamic protocols need to run asymptotically slowly in order to avoid dissipative losses. By adapting ideas from counter-diabatic driving and Floquet engineering to open systems, we develop fast-forward protocols for swiftly thermalizing a system oscillator locally coupled to an optical phonon bath. These protocols control the system frequency and the system-bath coupling to induce a resonant state exchange between the system and the bath. We apply the fast-forward protocols to realize a fast approximate Otto engine operating at high power near the Carnot Efficiency. Our results suggest design principles for swift cooling protocols in coupled many-body systems.

preprint2019arXiv

Variational Schrieffer-Wolff Transformations for Quantum Many-Body Dynamics

Building on recent results for adiabatic gauge potentials, we propose a variational approach for computing the generator of Schrieffer-Wolff transformations. These transformations consist of block diagonalizing a Hamiltonian through a unitary rotation, which leads to effective dynamics in a computationally tractable reduced Hilbert space. The generator of these rotations are computed variationally and thus go beyond standard perturbative methods; the error is controlled by the locality of the variational ansatz. The method is demonstrated on two models. First, in the attractive Fermi-Hubbard model with on-site disorder, we find indications of a lack of observable many-body localization in the thermodynamic limit due to the inevitable mixture of different spinon sectors. Second, in the low-energy sector of the XY spin model with a broken U(1)-symmetry, we analyze ground state response functions by combining the variational SW transformation with the truncated spectrum approach.

preprint2018arXiv

Quantum Hydrodynamics in Spin Chains with Phase Space Methods

Connecting short time microscopic dynamics with long time hydrodynamics in strongly correlated quantum systems is one of the outstanding questions. In particular, it is very difficult to determine various hydrodynamic coefficients like the diffusion constant or viscosity starting from a microscopic model: exact quantum simulations are limited to either small system sizes or to short times, which are insufficient to reach asymptotic behavior. In this Letter, we show that these difficulties, at least for a particular model, can be circumvented by using the cluster truncated Wigner approximation (CTWA), which maps quantum Hamiltonian dynamics into classical Hamiltonian dynamics in auxiliary high-dimensional phase space. We apply CTWA to a XXZ next-nearest-neighbor spin 1/2 chain and find behavior consisting of short time spin relaxation which gradually crosses over to emergent diffusive behavior at long times. For a random initial state we show that CTWA correctly reproduces the whole spin spectral function. Necessary in this construction is sampling from properly fluctuating initial conditions: the Dirac mean-field (variational) ansatz, which neglects such fluctuations, leads to incorrect predictions.

preprint2016arXiv

Enabling Adiabatic Passages Between Disjoint Regions in Parameter Space through Topological Transitions

We explore topological transitions in parameter space in order to enable adiabatic passages between regions adiabatically disconnected within a given parameter manifold. To this end, we study the Hamiltonian of two coupled qubits interacting with external magnetic fields, and make use of the analogy between the Berry curvature and magnetic fields in parameter space, with spectrum degeneracies associated to magnetic charges. Symmetry-breaking terms induce sharp topological transitions on these charge distributions, and we show how one can exploit this effect to bypass crossing degeneracies. We also investigate the curl of the Berry curvature, an interesting but as of yet not fully explored object, which together with its divergence uniquely defines this field. Finally, we suggest a simple method for measuring the Berry curvature, thereby showing how one can experimentally verify our results.

preprint2016arXiv

From Quantum Chaos and Eigenstate Thermalization to Statistical Mechanics and Thermodynamics

This review gives a pedagogical introduction to the eigenstate thermalization hypothesis (ETH), its basis, and its implications to statistical mechanics and thermodynamics. In the first part, ETH is introduced as a natural extension of ideas from quantum chaos and random matrix theory (RMT). To this end, we present a brief overview of classical and quantum chaos, as well as RMT and some of its most important predictions. The latter include the statistics of energy levels, eigenstate components, and matrix elements of observables. Building on these, we introduce the ETH and show that it allows one to describe thermalization in isolated chaotic systems without invoking the notion of an external bath. We examine numerical evidence of eigenstate thermalization from studies of many-body lattice systems. We also introduce the concept of a quench as a means of taking isolated systems out of equilibrium, and discuss results of numerical experiments on quantum quenches. The second part of the review explores the implications of quantum chaos and ETH to thermodynamics. Basic thermodynamic relations are derived, including the second law of thermodynamics, the fundamental thermodynamic relation, fluctuation theorems, and the Einstein and Onsager relations. In particular, it is shown that quantum chaos allows one to prove these relations for individual Hamiltonian eigenstates and thus extend them to arbitrary stationary statistical ensembles. We then show how one can use these relations to obtain nontrivial universal energy distributions in continuously driven systems. At the end of the review, we briefly discuss the relaxation dynamics and description after relaxation of integrable quantum systems, for which ETH is violated. We introduce the concept of the generalized Gibbs ensemble, and discuss its connection with ideas of prethermalization in weakly interacting systems.

preprint2016arXiv

Geodesic Paths for Quantum Many-Body Systems

We propose a method to obtain optimal protocols for adiabatic ground-state preparation near the adiabatic limit, extending earlier ideas from [D. A. Sivak and G. E. Crooks, Phys. Rev. Lett. 108, 190602 (2012)] to quantum non-dissipative systems. The space of controllable parameters of isolated quantum many-body systems is endowed with a Riemannian quantum metric structure, which can be exploited when such systems are driven adiabatically. Here, we use this metric structure to construct optimal protocols in order to accomplish the task of adiabatic ground-state preparation in a fixed amount of time. Such optimal protocols are shown to be geodesics on the parameter manifold, maximizing the local fidelity. Physically, such protocols minimize the average energy fluctuations along the path. Our findings are illustrated on the Landau-Zener model and the anisotropic XY spin chain. In both cases we show that geodesic protocols drastically improve the final fidelity. Moreover, this happens even if one crosses a critical point, where the adiabatic perturbation theory fails.

preprint2016arXiv

Heating and many-body resonances in a periodically driven two-band system

We study the dynamics and stability in a strongly interacting resonantly driven two-band model. Using exact numerical simulations, we find a stable regime at large driving frequencies where the time evolution is governed by a local Floquet Hamiltonian that is approximately conserved out to very long times. For slow driving, on the other hand, the system becomes unstable and heats up to infinite temperature. While thermalization is relatively fast in these two regimes (but to different "temperatures"), in the crossover between them we find slow nonthermalizing time evolution: temporal fluctuations become strong and temporal correlations long lived. Microscopically, we trace back the origin of this nonthermalizing time evolution to the properties of rare Floquet many-body resonances, whose proliferation at lower driving frequency removes the approximate energy conservation, and thus produces thermalization to infinite temperature.

preprint2016arXiv

Schrieffer-Wolff Transformation for Periodically Driven Systems: Strongly Correlated Systems with Artificial Gauge Fields

We generalize the Schrieffer-Wolff transformation to periodically driven systems using Floquet theory. The method is applied to the periodically driven, strongly interacting Fermi-Hubbard model, for which we identify two regimes resulting in different effective low-energy Hamiltonians. In the nonresonant regime, we realize an interacting spin model coupled to a static gauge field with a nonzero flux per plaquette. In the resonant regime, where the Hubbard interaction is a multiple of the driving frequency, we derive an effective Hamiltonian featuring doublon association and dissociation processes. The ground state of this Hamiltonian undergoes a phase transition between an ordered phase and a gapless Luttinger liquid phase. One can tune the system between different phases by changing the amplitude of the periodic drive.

preprint2016arXiv

Slow quenches in a quantum Ising chain; dynamical phase transitions and topology

We study the slow quenching dynamics (characterized by an inverse rate, $τ^{-1}$) of a one-dimensional transverse Ising chain with nearest neighbor ferromagentic interactions across the quantum critical point (QCP) and analyze the Loschmidt overlap {measured using the subsequent temporal evolution of the final wave function (reached at the end of the quenching) with the final time-independent Hamiltonian}. Studying the Fisher zeros of the corresponding generalized "partition function", we probe non-analyticities manifested in the rate function of the return probability known as dynamical phase transitions (DPTs). In contrast to the sudden quenching case, we show that DPTs survive {in the subsequent temporal evolution following the quenching across two critical points of the model for a sufficiently slow rate; furthermore, an interesting "lobe" structure of Fisher zeros emerge.} We have also made a connection to topological aspects studying the dynamical topological order parameter ($ν_D(t)$), as a function of time ($t$) {measured from the instant when the quenching is complete. Remarkably, the time evolution of $ν_D(t)$ exhibits drastically different behavior following quenches across a single QCP and two QCPs. } {In the former case, $ν_D (t)$ increases step-wise by unity at every DPT (i.e., $Δν_D =1$). In the latter case, on the other hand, $ν_D(t)$ essentially oscillates between 0 and 1 (i.e., successive DPTs occur with $Δν_D =1$ and $Δν_D =-1$, respectively), except for instants where it shows a sudden jump by a factor of unity when two successive DPTs carry a topological charge of same sign.

preprint2015arXiv

Dynamic trapping near a quantum critical point

The study of dynamics in closed quantum systems has recently been revitalized by the emergence of experimental systems that are well-isolated from their environment. In this paper, we consider the closed-system dynamics of an archetypal model: spins near a second order quantum critical point, which are traditionally described by the Kibble-Zurek mechanism. Imbuing the driving field with Newtonian dynamics, we find that the full closed system exhibits a robust new phenomenon -- dynamic critical trapping -- in which the system is self-trapped near the critical point due to efficient absorption of field kinetic energy by heating the quantum spins. We quantify limits in which this phenomenon can be observed and generalize these results by developing a Kibble-Zurek scaling theory that incorporates the dynamic field. Our findings can potentially be interesting in the context of early universe physics, where the role of the driving field is played by the inflaton or a modulus.

preprint2015arXiv

Dynamical Stability of a Many-body Kapitza Pendulum

We consider a many-body generalization of the Kapitza pendulum: the periodically-driven sine-Gordon model. We show that this interacting system is dynamically stable to periodic drives with finite frequency and amplitude. This finding is in contrast to the common belief that periodically-driven unbounded interacting systems should always tend to an absorbing infinite-temperature state. The transition to an unstable absorbing state is described by a change in the sign of the kinetic term in the effective Floquet Hamiltonian and controlled by the short-wavelength degrees of freedom. We investigate the stability phase diagram through an analytic high-frequency expansion, a self-consistent variational approach, and a numeric semiclassical calculations. Classical and quantum experiments are proposed to verify the validity of our results.

preprint2015arXiv

Negative mass corrections in a dissipative stochastic environment

We study the dynamics of a macroscopic object interacting with a dissipative stochastic environment using an adiabatic perturbation theory. The perturbation theory reproduces known expressions for the friction coefficient and, surprisingly, gives an additional negative mass correction. The effect of the negative mass correction is illustrated by studying a harmonic oscillator interacting with a dissipative stochastic environment. While it is well known that the friction coefficient causes a reduction of the oscillation frequency we show that the negative mass correction can lead to its enhancement. By studying an exactly solvable model of a magnet coupled to a spin environment evolving under standard non-conserving dynamics we show that the effect is present even beyond the validity of the adiabatic perturbation theory.

preprint2015arXiv

Quantum versus classical annealing: insights from scaling theory and results for spin glasses on 3-regular graphs

We discuss an Ising spin glass where each $S=1/2$ spin is coupled antiferromagnetically to three other spins (3-regular graphs). Inducing quantum fluctuations by a time-dependent transverse field, we use out-of-equilibrium quantum Monte Carlo simulations to study dynamic scaling at the quantum glass transition. Comparing the dynamic exponent and other critical exponents with those of the classical (temperature-driven) transition, we conclude that quantum annealing is less efficient than classical simulated annealing in bringing the system into the glass phase. Quantum computing based on the quantum annealing paradigm is therefore inferior to classical simulated annealing for this class of problems. We also comment on previous simulations where a parameter is changed with the simulation time, which is very different from the true Hamiltonian dynamics simulated here.

preprint2015arXiv

SU(3) semiclassical representation of quantum dynamics of interacting spins

We present a formalism for simulating quantum dynamics of lattice spin-one systems by first introducing local hidden variables and then doing semiclassical (truncated Wigner) approximation in the extended phase space. In this way we exactly take into account the local on-site Hamiltonian and approximately treat spin-spin interactions. In particular, we represent each spin with eight classical SU(3) variables. Three of them represent the usual spin components and five others are hidden variables representing local spin-spin correlations. We compare our formalism with exact quantum dynamics of fully connected spin systems and find very good agreement. As an application we discuss quench dynamics of a Bose-Hubbard model near the superfluid-insulator transition for a 3D lattice system consisting of 1000 sites.

preprint2015arXiv

Universal High-Frequency Behavior of Periodically Driven Systems: from Dynamical Stabilization to Floquet Engineering

We give a general overview of the high-frequency regime in periodically driven systems and identify three distinct classes of driving protocols in which the infinite-frequency Floquet Hamiltonian is not equal to the time-averaged Hamiltonian. These classes cover systems, such as the Kapitza pendulum, the Harper-Hofstadter model of neutral atoms in a magnetic field, the Haldane Floquet Chern insulator and others. In all setups considered, we discuss both the infinite-frequency limit and the leading finite-frequency corrections to the Floquet Hamiltonian. We provide a short overview of Floquet theory focusing on the gauge structure associated with the choice of stroboscopic frame and the differences between stroboscopic and non-stroboscopic dynamics. In the latter case one has to work with dressed operators representing observables and a dressed density matrix. We also comment on the application of Floquet Theory to systems described by static Hamiltonians with well-separated energy scales and, in particular, discuss parallels between the inverse-frequency expansion and the Schrieffer-Wolff transformation extending the latter to driven systems.

preprint2014arXiv

Dynamic scaling at classical phase transitions approached through non-equilibrium quenching

We use Monte Carlo simulations to demonstrate generic scaling aspects of classical phase transitions approached through a quench (or annealing) protocol where the temperature changes as a function of time with velocity $v$. Using a generalized Kibble-Zurek ansatz, we demonstrate dynamic scaling for different types of stochastic dynamics (Metropolis, Swendsen-Wang, and Wolff) on Ising models in two and higher dimensions. We show that there are dual scaling functions governing the dynamic scaling, which together describe the scaling behavior in the entire velocity range $v \in [0, \infty)$. These functions have asymptotics corresponding to the adiabatic and diabatic limit, and close to these limits they are perturbative in $v$ and $1/v$, respectively. Away from their perturbative domains, both functions cross over into the same universal power-law scaling form governed by the static and dynamic critical exponents (as well as an exponent characterizing the quench protocol). As a by-product of the scaling studies, we obtain high-precision estimates of the dynamic exponent $z$ for the two-dimensional Ising model subject to the three variants of Monte Carlo dynamics; for single-spin Metropolis updates $z_M=2.1767(5)$, for Swendsen-Wang multi-cluster updates $z_{\rm SW}=0.297(3)$, and for Wolff single-cluster updates $z_{W}=0.30(2)$. For Wolff dynamics, we find an interesting behavior with a non-analytic break-down of the quasi-adiabatic and diabatic scaling, instead of the generic smooth cross-over described by a power law. We interpret this disconnect between the two scaling regimes as a dynamic phase transition of the Wolff algorithm, caused by an effective sudden loss of ergodicity at high velocity.

preprint2014arXiv

Emergent Newtonian dynamics and the geometric origin of mass

We consider a set of macroscopic (classical) degrees of freedom coupled to an arbitrary many-particle Hamiltonian system, quantum or classical. These degrees of freedom can represent positions of objects in space, their angles, shape distortions, magnetization, currents and so on. Expanding their dynamics near the adiabatic limit we find the emergent Newton's second law (force is equal to the mass times acceleration) with an extra dissipative term. In systems with broken time reversal symmetry there is an additional Coriolis type force proportional to the Berry curvature. We give the microscopic definition of the mass tensor relating it to the non-equal time correlation functions in equilibrium or alternatively expressing it through dressing by virtual excitations in the system. In the classical (high-temperature) limit the mass tensor is given by the product of the inverse temperature and the Fubini-Study metric tensor determining the natural distance between the eigenstates of the Hamiltonian. For free particles this result reduces to the conventional definition of mass. This finding shows that any mass, at least in the classical limit, emerges from the distortions of the Hilbert space highlighting deep connections between any motion (not necessarily in space) and geometry. We illustrate our findings with four simple examples.

preprint2014arXiv

Measuring a topological transition in an artificial spin 1/2 system

We present measurements of a topological property, the Chern number ($C_\mathrm{1}$), of a closed manifold in the space of two-level system Hamiltonians, where the two-level system is formed from a superconducting qubit. We manipulate the parameters of the Hamiltonian of the superconducting qubit along paths in the manifold and extract $C_\mathrm{1}$ from the nonadiabitic response of the qubit. By adjusting the manifold such that a degeneracy in the Hamiltonian passes from inside to outside the manifold, we observe a topological transition $C_\mathrm{1} = 1 \rightarrow 0$. Our measurement of $C_\mathrm{1}$ is quantized to within 2 percent on either side of the transition.

preprint2014arXiv

Stroboscopic versus non-stroboscopic dynamics in the Floquet realization of the Harper-Hofstadter Hamiltonian

We study the stroboscopic and non-stroboscopic dynamics in the Floquet realization of the Harper-Hofstadter Hamiltonian. We show that the former produces the evolution expected in the high-frequency limit only for observables which commute with the operator to which the driving protocol couples. On the contrary, non-stroboscopic dynamics is capable of capturing the evolution governed by the Floquet Hamiltonian of any observable associated with the effective high-frequency model. We provide exact numerical simulations for the dynamics of the density operator following a quantum cyclotron orbit on a $2\times 2$ plaquette, as well as the chiral current operator flowing along the legs of a $2\times 20$ ladder. The exact evolution is compared with its stroboscopic and non-stroboscopic counterparts, including finite-frequency corrections.

preprint2014arXiv

The Second Law of Thermodynamics under Unitary Evolution and External Operations

A microscopic definition of the thermodynamic entropy in an isolated quantum system must satisfy (i) additivity, (ii) extensivity and (iii) the second law of thermodynamics. We show that the diagonal entropy, which is the Shannon entropy in the energy eigenbasis at each instant of time, meets the first two requirements and that the third requirement is satisfied if an arbitrary external operation is performed at typical times. In terms of the diagonal entropy, thermodynamic irreversibility follows from the facts that the Hamiltonian dynamics restricts quantum trajectories under unitary evolution and that the external operation is performed without referring to any particular information about the microscopic state of the system.

preprint2013arXiv

Classifying and measuring the geometry of the quantum ground state manifold

From the Aharonov-Bohm effect to general relativity, geometry plays a central role in modern physics. In quantum mechanics many physical processes depend on the Berry curvature. However, recent advances in quantum information theory have highlighted the role of its symmetric counterpart, the quantum metric tensor. In this paper, we perform a detailed analysis of the ground state Riemannian geometry induced by the metric tensor, using the quantum XY chain in a transverse field as our primary example. We focus on a particular geometric invariant -- the Gaussian curvature -- and show how both integrals of the curvature within a given phase and singularities of the curvature near phase transitions are protected by critical scaling theory. For cases where the curvature is integrable, we show that the integrated curvature provides a new geometric invariant, which like the Chern number characterizes individual phases of matter. For cases where the curvature is singular, we classify three types -- integrable, conical, and curvature singularities -- and detail situations where each type of singularity should arise. Finally, to connect this abstract geometry to experiment, we discuss three different methods for measuring the metric tensor, namely via integrating a properly weighted noise spectral function and by using leading order responses of the work distribution to ramps and quenches in quantum many-body systems.

preprint2013arXiv

Dynamical Quantum Phase Transitions in the Transverse Field Ising Model

A phase transition indicates a sudden change in the properties of a large system. For temperature-driven phase transitions this is related to non-analytic behavior of the free energy density at the critical temperature: The knowledge of the free energy density in one phase is insufficient to predict the properties of the other phase. In this paper we show that a close analogue of this behavior can occur in the real time evolution of quantum systems, namely non-analytic behavior at a critical time. We denote such behavior a dynamical phase transition and explore its properties in the transverse field Ising model. Specifically, we show that the equilibrium quantum phase transition and the dynamical phase transition in this model are intimately related.

preprint2013arXiv

Efficiency bounds for nonequilibrium heat engines

We analyze the efficiency of thermal engines (either quantum or classical) working with a single heat reservoir like atmosphere. The engine first gets an energy intake, which can be done in arbitrary non-equilibrium way e.g. combustion of fuel. Then the engine performs the work and returns to the initial state. We distinguish two general classes of engines where the working body first equilibrates within itself and then performs the work (ergodic engine) or when it performs the work before equilibrating (non-ergodic engine). We show that in both cases the second law of thermodynamics limits their efficiency. For ergodic engines we find a rigorous upper bound for the efficiency, which is strictly smaller than the equivalent Carnot efficiency. I.e. the Carnot efficiency can be never achieved in single reservoir heat engines. For non-ergodic engines the efficiency can be higher and can exceed the equilibrium Carnot bound. By extending the fundamental thermodynamic relation to nonequilibrium processes, we find a rigorous thermodynamic bound for the efficiency of both ergodic and non-ergodic engines and show that it is given by the relative entropy of the non-equilibrium and initial equilibrium distributions.These results suggest a new general strategy for designing more efficient engines. We illustrate our ideas by using simple examples.

preprint2013arXiv

Many-body energy localization transition in periodically driven systems

According to the second law of thermodynamics the total entropy of a system is increased during almost any dynamical process. The positivity of the specific heat implies that the entropy increase is associated with heating. This is generally true both at the single particle level, like in the Fermi acceleration mechanism of charged particles reflected by magnetic mirrors, and for complex systems in everyday devices. Notable exceptions are known in noninteracting systems of particles moving in periodic potentials. Here the phenomenon of dynamical localization can prevent heating beyond certain threshold. The dynamical localization is known to occur both at classical (Fermi-Ulam model) and at quantum levels (kicked rotor). However, it was believed that driven ergodic systems will always heat without bound. Here, on the contrary, we report strong evidence of dynamical localization transition in periodically driven ergodic systems in the thermodynamic limit. This phenomenon is reminiscent of many-body localization in energy space.

preprint2013arXiv

Quasi-adiabatic quantum Monte Carlo algorithm for quantum evolution in imaginary time

We propose a quantum Monte Carlo (QMC) algorithm for non-equilibrium dynamics in a system with a parameter varying as a function of time. The method is based on successive applications of an evolving Hamiltonian to an initial state and delivers results for a whole range of the tuning parameter in a single run, allowing for access to both static and dynamic properties of the system. This approach reduces to the standard Schrödinger dynamics in imaginary time for quasi-adiabatic evolutions, i.e., including the leading non-adiabatic correction to the adiabatic limit. We here demonstrate this quasi-adiabatic QMC (QAQMC) method for linear ramps of the transverse-field Ising model across its quantum-critical point in one and two dimensions. The critical behavior can be described by generalized dynamic scaling. For the two-dimensional square-lattice system we use the method to obtain a high-precision estimate of the quantum-critical point $(h/J)_c=3.04463(12)$, where $h$ is the transverse magnetic field and $J$ the nearest-neighbor Ising coupling. The QAQMC method can also be used to extract the Berry curvature and the metric tensor.

preprint2013arXiv

Transport-induced correlations in weakly interacting systems

We study spatial correlations in the transport of energy between two baths at different temperatures. To do this, we introduce a minimal model in which energy flows from one bath to another through two subsystems. We show that the transport-induced energy correlations between the two subsystems are of the same order as the energy fluctuations within each subsystem. The correlations can be either positive or negative and we give bounds on their values which are associated with a dynamic energy scale. The different signs originate as a competition between fluctuations generated near the baths, and fluctuations of the current between the two subsystems. This interpretation sheds light on known results for spatially-dependent heat and particle conduction models.

preprint2013arXiv

Universal Rephasing Dynamics after a Quantum Quench via Sudden Coupling of Two Initially Independent Condensates

We consider a quantum quench in which two initially independent condensates are suddenly coupled, and study the subsequent "rephasing" dynamics. For weak couplings, the time-evolution of physical observables is predicted to follow universal scaling laws, connecting the short-time dynamics to the long-time non-perturbative regime. We first present a two-mode model valid in two and three dimensions and then move to one dimension, where the problem is described by a gapped Sine-Gordon theory. Combining analytical and numerical methods, we compute universal time-dependent expectation values, allowing a quantitative comparison with future experiments.

preprint2012arXiv

Entropy of isolated quantum systems after a quench

A diagonal entropy, which depends only on the diagonal elements of the system's density matrix in the energy representation, has been recently introduced as the proper definition of thermodynamic entropy in out-of-equilibrium quantum systems. We study this quantity after an interaction quench in lattice hard-core bosons and spinless fermions, and after a local chemical potential quench in a system of hard-core bosons in a superlattice potential. The former systems have a chaotic regime, where the diagonal entropy becomes equivalent to the equilibrium microcanonical entropy, coinciding with the onset of thermalization. The latter system is integrable. We show that its diagonal entropy is additive and different from the entropy of a generalized Gibbs ensemble, which has been introduced to account for the effects of conserved quantities at integrability.

preprint2012arXiv

Geometric phase contribution to quantum non-equilibrium many-body dynamics

We study the influence of geometry of quantum systems underlying space of states on its quantum many-body dynamics. We observe an interplay between dynamical and topological ingredients of quantum non-equilibrium dynamics revealed by the geometrical structure of the quantum space of states. As a primary example we use the anisotropic XY ring in a transverse magnetic field with an additional time-dependent flux. In particular, if the flux insertion is slow, non-adiabatic transitions in the dynamics are dominated by the dynamical phase. In the opposite limit geometric phase strongly affects transition probabilities. We show that this interplay can lead to a non-equilibrium phase transition between these two regimes. We also analyze the effect of geometric phase on defect generation during crossing a quantum critical point.

preprint2012arXiv

Quantum phase slips in one-dimensional superfluids in a periodic potential

We study the decay of superflow of a one-dimensional (1D) superfluid in the presence of a periodic potential. In 1D, superflow at zero temperature can decay via quantum nucleation of phase slips even when the flow velocity is much smaller than the critical velocity predicted by mean-field theories. Applying the instanton method to the O(2) quantum rotor model, we calculate the nucleation rate of quantum phase slips $Γ$. When the flow momentum $p$ is small, we find that the nucleation rate per unit length increases algebraically with $p$ as $Γ/L \propto p^{2K-2}$, where $L$ is the system size and $K$ is the Tomonaga-Luttinger parameter. Based on the relation between the nucleation rate and the quantum superfluid-insulator transition, we present a unified explanation on the scaling formulae of the nucleation rate for periodic, disorder, and single-barrier potentials. Using the time-evolving block decimation method, we compute the exact quantum dynamics of the superflow decay in the 1D Bose-Hubbard model at unit filling. From the numerical analyses, we show that the scaling formula is valid for the case of the Bose-Hubbard model, which can quantitatively describe Bose gases in optical lattices.

preprint2012arXiv

Weak and strong typicality in quantum systems

We study the properties of mixed states obtained from eigenstates of many-body lattice Hamiltonians after tracing out part of the lattice. Two scenarios emerge for generic systems: (i) the diagonal entropy becomes equivalent to the thermodynamic entropy when a few sites are traced out (weak typicality); and (ii) the von Neumann (entanglement) entropy becomes equivalent to the thermodynamic entropy when a large fraction of the lattice is traced out (strong typicality). Remarkably, the results for few-body observables obtained with the reduced, diagonal, and canonical density matrices are very similar to each other, no matter which fraction of the lattice is traced out. Hence, for all physical quantities studied here, the results in the diagonal ensemble match the thermal predictions.

preprint2011arXiv

Localized phase structures growing out of quantum fluctuations in a quench of tunnel-coupled atomic condensates

We investigate the relative phase between two weakly interacting 1D condensates of bosonic atoms after suddenly switching on the tunnel-coupling. The following phase dynamics is governed by the quantum sine-Gordon equation. In the semiclassical limit of weak interactions, we observe the parametric amplification of quantum fluctuations leading to the formation of breathers with a finite lifetime. The typical lifetime and density of the these 'quasibreathers' are derived employing exact solutions of the classical sine-Gordon equation. Both depend on the initial relative phase between the condensates, which is considered as a tunable parameter.

preprint2011arXiv

Nonequilibrium dynamics of closed interacting quantum systems

This colloquium gives an overview of recent theoretical and experimental progress in the area of nonequilibrium dynamics of isolated quantum systems. We particularly focus on quantum quenches: the temporal evolution following a sudden or slow change of the coupling constants of the system Hamiltonian. We discuss several aspects of the slow dynamics in driven systems and emphasize the universality of such dynamics in gapless systems with specific focus on dynamics near continuous quantum phase transitions. We also review recent progress on understanding thermalization in closed systems through the eigenstate thermalization hypothesis and discuss relaxation in integrable systems. Finally we overview key experiments probing quantum dynamics in cold atom systems and put them in the context of our current theoretical understanding.

preprint2011arXiv

Oscillating fidelity susceptibility near a quantum multicritical point

We study scaling behavior of the geometric tensor $χ_{α,β}(λ_1,λ_2)$ and the fidelity susceptibility $(χ_{\rm F})$ in the vicinity of a quantum multicritical point (MCP) using the example of a transverse XY model. We show that the behavior of the geometric tensor (and thus of $χ_{\rm F}$) is drastically different from that seen near a critical point. In particular, we find that is highly non-monotonic function of $λ$ along the generic direction $λ_1\simλ_2 = λ$ when the system size $L$ is bounded between the shorter and longer correlation lengths characterizing the MCP: $1/|λ|^{ν_1}\ll L\ll 1/|λ|^{ν_2}$, where $ν_1<ν_2$ are the two correlation length exponents characterizing the system. We find that the scaling of the maxima of the components of $χ_{αβ}$ is associated with emergence of quasi-critical points at $λ\sim 1/L^{1/ν_1}$, related to the proximity to the critical line of finite momentum anisotropic transition. This scaling is different from that in the thermodynamic limit $L\gg 1/|λ|^{ν_2}$, which is determined by the conventional critical exponents. We use our results to calculate the defect density following a rapid quench starting from the MCP and show that it exerts a step-like behavior for small quench amplitudes. Study of heat density and diagonal entropy density also show signatures of quasi-critical points.

preprint2011arXiv

Prethermalization in quenched spinor condensates

Motivated by recent experiments, we consider the dynamics of spin-one spinor condensates after a quantum quench from the polar to ferromagnetic state from varying the quadratic Zeeman field q. We apply the Truncated Wigner Approximation (TWA) to the spinor system, including all spatial and spin degrees of freedom. For short times, we find full agreement with the linearized Bogoliubov analysis. For longer times, where the Bogoliubov theory fails, we find that the system reaches a quasi-steady prethermalized state. We compute the Bogoliubov spectrum about the ferromagnetic state with general q and show that the resulting finite temperature correlation functions grossly disagree with the full TWA results, thus indicating that the system does not thermalize even over very long time scales. Finally we show that the absence of thermalization over realistic time scales is consistent with calculations of Landau damping rates of excitations in the finite-temperature condensate.

preprint2011arXiv

Superfluid to Mott insulator transition in the one-dimensional Bose-Hubbard model for arbitrary integer filling factors

We study the quantum phase transition between the superfluid and the Mott insulator in the one-dimensional (1D) Bose-Hubbard model. Using the time-evolving block decimation method, we numerically calculate the tunneling splitting of two macroscopically distinct states with different winding numbers. From the scaling of the tunneling splitting with respect to the system size, we determine the critical point of the superfluid to Mott insulator transition for arbitrary integer filling factors. We find that the critical values versus the filling factor in 1D, 2D, and 3D are well approximated by a simple analytical function. We also discuss the condition for determining the transition point from a perspective of the instanton method.

preprint2011arXiv

Universal energy fluctuations in thermally isolated driven systems

When an isolated system is brought in contact with a heat bath its final energy is random and follows the Gibbs distribution -- a cornerstone of statistical physics. The system's energy can also be changed by performing non-adiabatic work using a cyclic process. Almost nothing is known about the resulting energy distribution in this setup, which is especially relevant to recent experimental progress in cold atoms, ions traps, superconducting qubits and other systems. Here we show that when the non-adiabatic process comprises of many repeated cyclic processes the resulting energy distribution is universal and different from the Gibbs ensemble. We predict the existence of two qualitatively different regimes with a continuous second order like transition between them. We illustrate our approach performing explicit calculations for both interacting and non-interacting systems.

preprint2010arXiv

Microscopic diagonal entropy and its connection to basic thermodynamic relations

We define a diagonal entropy (d-entropy) for an arbitrary Hamiltonian system as $S_d=-\sum_n ρ_{nn}\ln ρ_{nn}$ with the sum taken over the basis of instantaneous energy states. In equilibrium this entropy coincides with the conventional von Neumann entropy $S_n=-{\rm Tr}\, ρ\lnρ$. However, in contrast to $S_n$, the d-entropy is not conserved in time in closed Hamiltonian systems. If the system is initially in stationary state then in accord with the second law of thermodynamics the d-entropy can only increase or stay the same. We also show that the d-entropy can be expressed through the energy distribution function and thus it is measurable, at least in principle. Under very generic assumptions of the locality of the Hamiltonian and non-integrability the d-entropy becomes a unique function of the average energy in large systems and automatically satisfies the fundamental thermodynamic relation. This relation reduces to the first law of thermodynamics for quasi-static processes. The d-entropy is also automatically conserved for adiabatic processes. We illustrate our results with explicit examples and show that $S_d$ behaves consistently with expectations from thermodynamics.

preprint2010arXiv

Phase space representation of quantum dynamics

We discuss a phase space representation of quantum dynamics of systems with many degrees of freedom. This representation is based on a perturbative expansion in quantum fluctuations around one of the classical limits. We explicitly analyze expansions around three such limits: (i) corpuscular or Newtonian limit in the coordinate-momentum representation, (ii) wave or Gross-Pitaevskii limit for interacting bosons in the coherent state representation, and (iii) Bloch limit for the spin systems. We discuss both the semiclassical (truncated Wigner) approximation and further quantum corrections appearing in the form of either stochastic quantum jumps along the classical trajectories or the nonlinear response to such jumps. We also discuss how quantum jumps naturally emerge in the analysis of non-equal time correlation functions. This representation of quantum dynamics is closely related to the phase space methods based on the Wigner-Weyl quantization and to the Keldysh technique. We show how such concepts as the Wigner function, Weyl symbol, Moyal product, Bopp operators, and others automatically emerge from the Feynmann's path integral representation of the evolution in the Heisenberg representation. We illustrate the applicability of this expansion with various examples mostly in the context of cold atom systems including sine-Gordon model, one- and two-dimensional Bose Hubbard model, Dicke model and others.

preprint2010arXiv

Quantum damping of Fermi-Pasta-Ulam revivals in ultracold Bose gases

We propose an experimental scheme for studying the Fermi-Pasta-Ulam (FPU) phenomenon in a quantum mechanical regime using ultracold atoms. Specifically, we suggest and analyze a setup of one-dimensional Bose gases confined into an optical lattice. The strength of quantum fluctuations is controlled by tuning the number of atoms per lattice sites (filling factor). By simulating the real-time dynamics of the Bose-Hubbard model by means of the exact numerical method of time-evolving block decimation, we investigate the effects of quantum fluctuations on the FPU recurrence and show that strong quantum fluctuations cause significant damping of the FPU oscillation.

preprint2010arXiv

Universal Dynamics Near Quantum Critical Points

We give an overview of the scaling of density of quasi-particles and excess energy (heat) for nearly adiabatic dynamics near quantum critical points (QCPs). In particular we discuss both sudden quenches of small amplitude and slow sweeps across the QCP. We show close connection between universal scaling of these quantities with the scaling behavior of the fidelity susceptibility and its generalizations. In particular we argue that the Kibble-Zurek scaling can be easily understood using this concept. We discuss how these scalings can be derived within the adiabatic perturbation theory and how using this approach slow and fast quenches can be treated within the same framework. We also describe modifications of these scalings for finite temperature quenches and emphasize the important role of statistics of low-energy excitations. In the end we mention some connections between adiabatic dynamics near critical points with dynamics associated with space-time singularities in the metrics, which naturally emerges in such areas as cosmology and string theory.

preprint2009arXiv

Accurate numerical verification of the instanton method for macroscopic quantum tunneling: dynamics of phase slips

Instanton methods, in which imaginary-time evolution gives the tunneling rate, have been widely used for studying quantum tunneling in various contexts. Nevertheless, how accurate instanton methods are for the problems of macroscopic quantum tunneling (MQT) still remains unclear because of lack of their direct comparison with exact time evolution of the many-body Schroedinger equation. Here, we verify instanton methods applied to coherent MQT. Specifically applying the quasi-exact numerical method of time-evolving block decimation to the system of bosons in a ring lattice, we directly simulate the real-time quantum dynamics of supercurrents, where a coherent oscillation between two macroscopically distinct current states occurs due to MQT. The tunneling rate extracted from the coherent oscillation is compared with that given by the instanton method. We show that the error is within 10% when the effective Planck's constant is sufficiently small. We also discuss phase slip dynamics associated with the coherent oscillations.

preprint2009arXiv

Breakdown of macroscopic quantum self-trapping in coupled mesoscopic one dimensional Bose gases

Two coupled BECs with a large population imbalance exhibit macroscopic quantum self-trapping (MQST) if the ratio of interaction energy to tunneling energy is above a critical value. Here we investigate effect of quantum fluctuations on MQST. In particular, we analyze the dynamics of a system of two elongated Bose gases prepared with a large population imbalance, where due to the quasi one dimensional character of the gas we no longer have true long range order, and the effect of quantum fluctuations is much more important. We show that MQST is possible in this system, but even when it is achieved it is not always dynamically stable. Using this instability one can construct states with sharply peaked momentum distributions around characteristic momenta dependent on system parameters. Our other finding is the nonmonotonic oscillating dependence of the decay rate of the MQST on the length of the wires. We also address the interesting question of thermalization in this system and show that it occurs only in sufficiently long wires.

preprint2009arXiv

Superfluid-insulator transition of disordered bosons in one-dimension

We study the superfluid-insulator transition in a one dimensional system of interacting bosons, modeled as a disordered Josephson array, using a strong randomness real space renormalization group technique. Unlike perturbative methods, this approach does not suffer from run-away flows and allows us to study the complete phase diagram. We show that the superfluid insulator transition is always Kosterlitz- Thouless like in the way that length and time scales diverge at the critical point. Interestingly however, we find that the transition at strong disorder occurs at a non universal value of the Luttinger parameter, which depends on the disorder strength. This result places the transition in a universality class different from the weak disorder transition first analyzed by Giamarchi and Schulz [Europhys. Lett. {\bf 3}, 1287 (1987)]. While the details of the disorder potential are unimportant at the critical point, the type of disorder does influence the properties of the insulating phases. We find three classes of insulators which arise for different classes of disorder potential. For disorder only in the charging energies and Josephson coupling constants, at integer filling we find an incompressible but gapless Mott glass phase. If both integer and half integer filling factors are allowed then the corresponding phase is a random singlet insulator, which has a divergent compressibility. Finally in a generic disorder potential the insulator is a Bose glass with a finite compressibility.

preprint2007arXiv

Spectroscopy of Collective Excitations in Interacting Low-Dimensional Many-Body Systems Using Quench Dynamics

We study the problem of rapid change of the interaction parameter (quench) in many-body low-dimensional system. It is shown that, measuring correlation functions after the quench the information about a spectrum of collective excitations in a system can be obtained. This observation is supported by analysis of several integrable models and we argue that it is valid for non-integrable models as well. Our conclusions are supplemented by performing exact numerical simulations on finite systems. We propose that measuring power spectrum in dynamically split 1D Bose-Einsten condensate into two coupled condensates can be used as experimental test of our predictions.

preprint2007arXiv

Unzipping Vortices in Type-II Superconductors

The unzipping of vortex lines using magnetic-force microscopy from extended defects is studied theoretically. We study both the unzipping isolated vortex from common defects, such as columnar pins and twin-planes, and the unzipping of a vortex from a plane in the presence of other vortices. We show, using analytic and numerical methods, that the universal properties of the unzipping transition of a single vortex depend only on the dimensionality of the defect in the presence and absence of disorder. For the unzipping of a vortex from a plane populated with many vortices is shown to be very sensitive to the properties of the vortices in the two-dimensional plane. In particular such unzipping experiments can be used to measure the ``Luttinger liquid parameter'' of the vortices in the plane. In addition we suggest a method for measuring the line tension of the vortex directly using the experiments.

preprint2006arXiv

Full quantum distribution of contrast in interference experiments between interacting one dimensional Bose liquids

We analyze interference experiments for a pair of independent one dimensional condensates of interacting bosonic atoms at zero temperature. We show that the distribution function of fringe amplitudes contains non-trivial information about non-local correlations within individual condensates and can be calculated explicitly using methods of conformal field theory. We point out interesting relations between these distribution functions, the partition function for a quantum impurity in a one-dimensional Luttinger liquid, and transfer matrices of conformal field theories. We demonstrate the connection between interference experiments in cold atoms and a variety of statistical models ranging from stochastic growth models to two dimensional quantum gravity. Such connection can be used to design a quantum simulator of unusual two-dimensional models described by nonunitary conformal field theories with negative central charges.