Source author record

Jacek Dziarmaga

Jacek Dziarmaga 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

42works
12topics
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

42 published item(s)

preprint2022arXiv

Kink correlations, domain size distribution, and the emptiness formation probability after the Kibble-Zurek quench in the quantum Ising chain

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.

preprint2022arXiv

Time evolution of an infinite projected entangled pair state: a gradient tensor update in the tangent space

Time evolution of an infinite 2D many body quantum lattice system can be described by the Suzuki-Trotter decomposition applied to the infinite projected entangled pair state (iPEPS). Each Trotter gate increases the bond dimension of the tensor network, $D$, that has to be truncated back in a way that minimizes a suitable error measure. This paper goes beyond simplified error measures -- like the one used in the full update (FU), the simple update (SU), and their intermediate neighborhood tensor update (NTU) -- and directly maximizes an overlap between the exact iPEPS with the increased bond dimension and the new iPEPS with the truncated one. The optimization is performed in a tangent space of the iPEPS variational manifold. This gradient tensor update (GTU) is benchmarked by a simulation of a sudden quench of a transverse field in the 2D quantum Ising model and the quantum Kibble-Zurek mechanism in the same 2D system.

preprint2021arXiv

Simulation of many body localization and time crystals in two dimensions with the neighborhood tensor update

The Heisenberg antiferromagnet with discrete disorder on an infinite square lattice is evolved in time from an initial Néel state. The simulation is performed with the neighborhood tensor update (NTU) algorithm for an infinite projected entangled pair state (iPEPS) [Phys. Rev. B 104, 094411 (2021)]. Ancillary spins are used to average over $2$ or $5$ discrete values of disorder. With a bond dimension up to $20$, evolution times are long enough to identify a many body localized regime for a strong enough disorder. Furthermore, the same Hamiltonian is subject to periodic spin flips. Simulations of the Floquet dynamics show that it can sustain a time crystalline stage for a strong enough disorder.

preprint2021arXiv

Tensor network study of the $m=1/2$ magnetization plateau in the Shastry-Sutherland model at finite temperature

The two-dimensional infinite projected entangled pair states tensor network is evolved in imaginary time with the full update (FU) algorithm to simulate the Shastry-Sutherland model in a magnetic field at finite temperature directly in the thermodynamic limit. We focus on the phase transition into the $m=1/2$ magnetization plateau, which was observed in experiments on SrCu$_2$(BO$_3$)$_2$. For the largest simulated bond dimension, the early evolution in the high-temperature regime is simulated with the simple update (SU) scheme and then, as the correlation length increases, continued with the FU scheme towards the critical regime. We apply a small symmetry-breaking bias field and then extrapolate towards zero bias using a simple scaling theory in the bias field. The combined SU + FU scheme provides an accurate estimate of the critical temperature, even though the results could not be fully converged in the bond dimension in the vicinity of the transition. The critical temperature estimate is improved with a generalized scaling theory that combines two divergent length scales: One due to the bias, and the other due to the finite bond dimension. The obtained results are consistent with the transition being in the universality class of the two-dimensional classical Ising model. The estimated critical temperature is $3.5(2)$ K, which is well above the temperature $2.1$ K used in the experiments.

preprint2020arXiv

Determining non-Abelian topological order from infinite projected entangled pair states

We generalize the method introduced in Phys. Rev. B 101, 041108 (2020) of extracting information about topological order from the ground state of a strongly correlated two-dimensional system represented by an infinite projected entangled pair state (iPEPS) to non-Abelian topological order. When wrapped on a torus the unique iPEPS becomes a superposition of degenerate and locally indistinguishable ground states. We find numerically symmetries of the iPEPS, represented by infinite matrix product operators (MPO), and their fusion rules. The rules tell us how to combine the symmetries into projectors onto states with well defined anyon flux. A linear structure of the MPO projectors allows for efficient determination for each state its second Renyi topological entanglement entropy on an infinitely long cylinder directly in the limit of infinite cylinder's width. The same projectors are used to compute topological $S$ and $T$ matrices encoding mutual- and self-statistics of emergent anyons. The algorithm is illustrated by examples of Fibonacci and Ising non-Abelian string net models.

preprint2020arXiv

Determining topological order from infinite projected entangled pair states

We present a method of extracting information about the topological order from the ground state of a strongly correlated two-dimensional system computed with the infinite projected entangled pair state (iPEPS). For topologically ordered systems, the iPEPS wrapped on a torus becomes a superposition of degenerate, locally indistinguishable ground states. Projectors in the form of infinite matrix product operators (iMPO) onto states with well-defined anyon flux are used to compute topological $S$ and $T$ matrices (encoding mutual- and self-statistics of emergent anyons). The algorithm is shown to be robust against a perturbation driving string-net toric code across a phase transition to a ferromagnetic phase. Our approach provides accurate results near quantum phase transition, where the correlation length is prohibitively large for other numerical methods. Moreover, we used numerically optimized iPEPS describing the ground state of the Kitaev honeycomb model in the toric code phase and obtained topological data in excellent agreement with theoretical prediction.

preprint2020arXiv

Kibble-Zurek Dynamics in a Trapped Ultracold Bose Gas

The dynamical evolution of an inhomogeneous ultracold atomic gas quenched at different controllable rates through the Bose-Einstein condensation phase transition is studied numerically in the premise of a recent experiment in an anisotropic harmonic trap. Our findings based on the stochastic (projected) Gross-Pitaevskii equation are shown to be consistent at early times with the predictions of the homogeneous Kibble-Zurek mechanism. This is demonstrated by collapsing the early dynamical evolution of densities, spectral functions and correlation lengths for different quench rates, based on an appropriate characterization of the distance to criticality felt by the quenched system. The subsequent long-time evolution, beyond the identified dynamical critical region, is also investigated by looking at the behaviour of the density wavefront evolution and the corresponding phase ordering dynamics.

preprint2020arXiv

Sonic horizons and causality in the phase transition dynamics

A system gradually driven through a symmetry-breaking phase transition is subject to the Kibble-Zurek mechanism (KZM). As a consequence of the critical slowing down, its state cannot follow local equilibrium, and its evolution becomes non-adiabatic near the critical point. In the simplest approximation, that stage can be regarded as "impulse" where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such "freeze-out" might suggest that the coherence length of the nascent order parameter remains unchanged as the critical region is traversed. By contrast, the original causality-based discussion emphasized the role of the {\it sonic horizon}: domains of the broken symmetry phase can expand with a velocity limited by the speed of the relevant sound. This effect was demonstrated in the quantum Ising chain where the dynamical exponent $z=1$ and quasiparticles excited by the transition have a fixed speed of sound. To elucidate the role of the sonic horizon, in this paper we study two systems with $z>1$ where the speed of sound is no longer fixed, and the fastest excited quasiparticles set the size of the sonic horizon. Their effective speed decays with the increasing transition time. In the extreme case, the dynamical exponent $z$ can diverge such as in the Griffiths region of the random Ising chain where localization of excited quasiparticles freezes the growth of the correlation range when the critical region is traversed. Of particular interest is an example with $z<1$ -- the long-range extended Ising chain, where there is no upper limit to the velocity of excited quasiparticles with small momenta. Initially, the power-law tail of the correlation function grows adiabatically, but in the non-adiabatic stage it lags behind the adiabatic evolution -- in accord with a generalized Lieb-Robinson bound.

preprint2020arXiv

Tensor-Network Approach for Quantum Metrology in Many-Body Quantum Systems

Identification of the optimal quantum metrological protocols in realistic many particle quantum models is in general a challenge that cannot be efficiently addressed by the state-of-the-art numerical and analytical methods. Here we provide a comprehensive framework exploiting matrix product operators (MPO) type tensor networks for quantum metrological problems. Thanks to the fact that the MPO formalism allows for an efficient description of short-range spatial and temporal noise correlations, the maximal achievable estimation precision in such models, as well as the optimal probe states in previously inaccessible regimes can be identified. Moreover, the application of infinite MPO (iMPO) techniques allows for a direct and efficient determination of the asymptotic precision of optimal protocols in the limit of infinite particle numbers. We illustrate the potential of our framework in terms of an atomic clock stabilization (temporal noise correlation) example as well as for magnetic field sensing in the presence of locally correlated magnetic field fluctuations (spatial noise correlations). As a byproduct, the developed methods for calculating the quantum Fisher information via MPOs may be used to calculate the fidelity susceptibility - a parameter widely used in many-body physics to study phase transitions.

preprint2016arXiv

Space-time renormalization in phase transition dynamics

When a system is driven across a quantum critical point at a constant rate its evolution must become non-adiabatic as the relaxation time $τ$ diverges at the critical point. According to the Kibble-Zurek mechanism (KZM), the emerging post-transition excited state is characterized by a finite correlation length $\hatξ$ set at the time $\hat t=\hat τ$ when the critical slowing down makes it impossible for the system to relax to the equilibrium defined by changing parameters. This observation naturally suggests a dynamical scaling similar to renormalization familiar from the equilibrium critical phenomena. We provide evidence for such KZM-inspired spatiotemporal scaling by investigating an exact solution of the transverse field quantum Ising chain in the thermodynamic limit.

preprint2016arXiv

Variational tensor network renormalization in imaginary time: benchmark results in the Hubbard model at finite temperature

A Gibbs operator $e^{-βH}$ for a 2D lattice system with a Hamiltonian $H$ can be represented by a 3D tensor network, the third dimension being the imaginary time (inverse temperature) $β$. Coarse-graining the network along $β$ results in an accurate 2D projected entangled-pair operator (PEPO) with a finite bond dimension. The coarse-graining is performed by a tree tensor network of isometries that are optimized variationally to maximize the accuracy of the PEPO. The algorithm is applied to the two-dimensional Hubbard model on an infinite square lattice. Benchmark results are obtained that are consistent with the best cluster dynamical mean-field theory and power series expansion in the regime of parameters where they yield mutually consistent results.

preprint2016arXiv

Variational tensor network renormalization in imaginary time: Two-dimensional quantum compass model at finite temperature

Progress in describing thermodynamic phase transitions in quantum systems is obtained by noticing that the Gibbs operator $e^{-βH}$ for a two-dimensional (2D) lattice system with a Hamiltonian $H$ can be represented by a three-dimensional tensor network, the third dimension being the imaginary time (inverse temperature) $β$. Coarse-graining the network along $β$ results in a 2D projected entangled-pair operator (PEPO) with a finite bond dimension $D$. The coarse-graining is performed by a tree tensor network of isometries. The isometries are optimized variationally --- taking into account full tensor environment --- to maximize the accuracy of the PEPO. The algorithm is applied to the isotropic quantum compass model on an infinite square lattice near a symmetry-breaking phase transition at finite temperature. From the linear susceptibility in the symmetric phase and the order parameter in the symmetry-broken phase the critical temperature is estimated at ${\cal T}_c=0.0606(4)J$, where $J$ is the isotropic coupling constant between $S=1/2$ pseudospins.

preprint2015arXiv

Projected Entangled Pair States at Finite Temperature: Iterative Self-Consistent Bond Renormalization for Exact Imaginary Time Evolution

A projected entangled pair state (PEPS) with ancillas can be evolved in imaginary time to obtain thermal states of a strongly correlated quantum system on a 2D lattice. Every application of a Suzuki-Trotter gate multiplies the PEPS bond dimension $D$ by a factor $k$. It has to be renormalized back to the original $D$. In order to preserve the accuracy of the Suzuki-Trotter (S-T) decomposition, the renormalization has in principle to take into account full environment made of the new tensors with the bond dimension $k\times D$. Here we propose a self-consistent renormalization procedure operating with the original bond dimension $D$, but without compromising the accuracy of the S-T decomposition. The iterative procedure renormalizes the bond using full environment made of renormalized tensors with the bond dimension $D$. After every renormalization, the new renormalized tensors are used to update the environment, and then the renormalization is repeated again and again until convergence. As a benchmark application, we obtain thermal states of the transverse field quantum Ising model on a square lattice - both infinite and finite - evolving the system across a second-order phase transition at finite temperature.

preprint2015arXiv

Second-order Peierls transition in the spin-orbital Kumar-Heisenberg model

We add a Heisenberg interaction term $\proptoλ$ in the one-dimensional SU(2)$\otimes$XY spin-orbital model introduced by B. Kumar. At $λ=0$ the spin and orbital degrees of freedom can be separated by a unitary transformation leading to an exact solution of the model. We show that a finite $λ>0$ leads to spontaneous dimerization of the system which in the thermodynamic limit becomes a smooth phase transition at $λ\to 0$, whereas it remains discontinuous within the first order perturbation approach. We present the behavior of the entanglement entropy, energy gap and dimerization order parameter in the limit of $λ\to 0$ confirming the critical behavior. Finally, we show the evidence of another phase transition in the Heisenberg limit, $λ\to\infty$, and give a qualitative analytical explanation of the observed dimerized states both in the limit of small and large $λ$.

preprint2015arXiv

Variational Approach to Projected Entangled Pair States at Finite Temperature

The projected entangled pair state (PEPS) ansatz can represent a thermal state in a strongly correlated system. We introduce a novel variational algorithm to optimize this tensor network. Since full tensor environment is taken into account, then with increasing bond dimension the optimized PEPS becomes the exact Gibbs state. Our presentation opens with a 1D version for a matrix product state (MPS) and then generalizes to PEPS in 2D. Benchmark results in the quantum Ising model are presented.

preprint2014arXiv

Exotic Spin Order due to Orbital Fluctuations

We investigate the phase diagrams of the spin-orbital $d^9$ Kugel-Khomskii model for increasing system dimensionality: from the square lattice monolayer, via the bilayer to the cubic lattice. In each case we find strong competition between different types of spin and orbital order, with entangled spin-orbital phases at the crossover from antiferromagnetic to ferromagnetic correlations in the intermediate regime of Hund's exchange. These phases have various types of exotic spin order and are stabilized by effective interactions of longer range which follow from enhanced spin-orbital fluctuations. We find that orbital order is in general more robust and spin order melts first under increasing temperature, as observed in several experiments for spin-orbital systems.

preprint2014arXiv

Fermionic Projected Entangled Pair States at Finite Temperature

An algorithm for imaginary time evolution of a fermionic projected entangled pair state (PEPS) with ancillas from infinite temperature down to a finite temperature state is presented. As a benchmark application, it is applied to spinless fermions hopping on a square lattice subject to $p$-wave pairing interactions. With a tiny bias it allows to evolve the system across a high-temperature continuous symmetry-breaking phase transition.

preprint2014arXiv

Quench in 1D Bose-Hubbard Model: Topological Defects and Excitations from Kosterlitz-Thouless Phase Transition Dynamics

Kibble-Zurek mechanism (KZM) uses critical scaling to predict density of topological defects and other excitations created in second order phase transitions. We point out that simply inserting asymptotic critical exponents deduced from the immediate vicinity of the critical point to obtain KZM predictions can lead to results that are inconsistent with a more careful KZM-like analysis based on causality -- on the comparison of the relaxation time of the order parameter with the time "distance" from the critical point. As a result, scaling of quench-generated excitations with quench rates can exhibit behavior that is locally (i.e., in the neighborhood of any given quench rate) well approximated by the power law, but with exponents that depend on that rate, and that are quite different from the naive prediction based on the critical exponents relevant for asymptotically long quench times. Kosterlitz-Thouless scaling (that governs e.g. Mott insulator to superfluid transition in the Bose-Hubbard model in one dimension) is investigated as an example of this phenomenon.

preprint2014arXiv

Striped critical spin liquid in a spin-orbital entangled RVB state in a projected entangled-pair state representation

We introduce a spin-orbital entangled (SOE) resonating valence bond (RVB) state on a square lattice of spins-$\frac12$ and orbitals represented by pseudospins-$\frac12$. Like the standard RVB state, it is a superposition of nearest-neighbor hard-core coverings of the lattice by spin singlets, but adjacent singlets are favoured to have perpendicular orientations and, more importantly, an orientation of each singlet is entangled with orbitals' state on its two lattice sites. The SOE-RVB state can be represented by a projected entangled pair state (PEPS) with a bond dimension $D=4$. This representation helps to reveal that the state is a superposition of striped coverings conserving a topological quantum number. The stripes are a critical quantum spin liquid. We propose a spin-orbital Hamiltonian supporting a SOE-RVB ground state.

preprint2013arXiv

Dynamics of the modified Kibble-Żurek mechanism in antiferromagnetic spin-1 condensates

We investigate the dynamics and outcome of a quantum phase transition from an antiferromagnetic to phase separated ground state in a spin-1 Bose-Einstein condensate of ultracold atoms. We explicitly demonstrate double universality in dynamics within experiments with various quench time. Furthermore, we show that spin domains created in the nonequilibrium transition constitute a set of mutually incoherent quasicondensates. The quasicondensates appear to be positioned in a semi-regular fashion, which is a result of the conservation of local magnetization during the post-selection dynamics.

preprint2013arXiv

Dynamics of the Mott Insulator to Superfluid quantum phase transition in the truncated Wigner approximation

The quantum phase transition from the Mott insulator state to the superfluid in the Bose-Hubbard model is investigated. We research one, two and three dimensional lattices in the truncated Wigner approximation. We compute both kinetic and potential energy and they turn out to have a power law behaviour as a function of the transition rate, with the power equal to 1/3. The same applies to the total energy in a system with a harmonic trap, which is usually present in the experimental set-up. These observations are in agreement with the experiment of [8], where such scalings were also observed and the power of the decay was numerically close to 1/3. The results confirm the Kibble-Zurek (adiabatic-impulse-adiabatic approximation) scenario for this transition.

preprint2013arXiv

Exotic spin orders driven by orbital fluctuations in the Kugel-Khomskii model

We study zero temperature phase diagram of the three-dimensional Kugel-Khomskii model on a cubic lattice using the cluster mean field theory and different perturbative expansions in the orbital sector. The phase diagram is rich, goes beyond the single-site mean field theory due to spin-orbital entanglement. In addition to the antiferromagnetic (AF) and ferromagnetic (FM) phases, one finds also a plaquette valence-bond phase with singlets ordered either on horizontal or vertical bonds. More importantly, for increasing Hund's exchange we identify three phases with exotic magnetic order stabilized by orbital fluctuations in between the AF and FM order: (i) an AF phase with two mutually orthogonal antiferromagnets on two sublattices in each $ab$ plane and AF order along the c axis (ortho-$G$-type phase), (ii) a canted-$A$-type AF phase with a non-trivial canting angle between nearest neighbor FM layers along the c axis, and (iii) a striped-AF phase with anisotropic AF order in the $ab$ planes. We elucidate the mechanism responsible for each of the above phases by deriving effective spin models which involve second and third neighbor Heisenberg interactions as well as four-site spin interactions going beyond Heisenberg physics, and explain how the entangled nearest neighbor spin-orbital superexchange generates spin interactions between more distant spins.

preprint2013arXiv

Topological Order in an Entangled SU(2)$\otimes$XY Spin-Orbital Ring

We present rigorous topological order which emerges in a one-dimensional spin-orbital model due to the ring topology. Although an exact solution of a spin-orbital ring with SU(2) spin and XY orbital interactions separates spins from orbitals by means of a unitary transformation, the spins are not independent when the ring is closed, but form two half-rings carrying opposite pseudomomenta. We show that an inverse transformation back to the physical degrees of freedom entangles the spin half-rings with the orbitals once again. This surprising correlation arises on changing the topology from an open to a closed chain, which reduces the degeneracy of the ground-state manifold, leaving in it only the states in which pseudomomenta compensate each other.

preprint2012arXiv

Double universality of a quantum phase transition in spinor condensates: the Kibble-Żurek mechanism and a conservation law

We consider a phase transition from antiferromagnetic to phase separated ground state in a spin-1 Bose-Einstein condensate of ultracold atoms. We demonstrate the occurrence of two scaling laws, for the number of spin fluctuations just after the phase transition, and for the number of spin domains in the final, stable configuration. Only the first scaling can be explained by the standard Kibble-Żurek mechanism. We explain the occurrence of two scaling laws by a model including post-selection of spin domains due to the conservation of condensate magnetization.

preprint2012arXiv

Noncollinear Magnetic Order Stabilized by Entangled Spin-Orbital Fluctuations

Quantum phase transitions in the two-dimensional Kugel-Khomski model on a square lattice are studied using the plaquette mean field theory and the entanglement renormalization ansatz. When $3z^2-r^2$ orbitals are favored by the crystal field and Hund's exchange is finite, both methods give a noncollinear exotic magnetic order which consists of four sublattices with mutually orthogonal nearest neighbor and antiferromagnetic second neighbor spins. We derive effective frustrated spin model with second and third neighbor spin interactions which stabilize this phase and follow from spin-orbital quantum fluctuations involving spin singlets entangled with orbital excitations.

preprint2012arXiv

Quench from Mott Insulator to Superfluid

We study a linear ramp of the nearest-neighbor tunneling rate in the Bose-Hubbard model driving the system from the Mott insulator state into the superfluid phase. We employ the truncated Wigner approximation to simulate linear quenches of a uniform system in 1,2, and 3 dimensions, and in a harmonic trap in 3 dimensions. In all these setups the excitation energy decays like one over third root of the quench time. The -1/3 scaling arises from an impulse-adiabatic approximation - a variant of the Kibble-Zurek mechanism - describing a crossover from non-adiabatic to adiabatic evolution when the system begins to keep pace with the increasing tunneling rate.

preprint2012arXiv

Topological Schrödinger cats: Non-local quantum superpositions of topological defects

Topological defects (such as monopoles, vortex lines, or domain walls) mark locations where disparate choices of a broken symmetry vacuum elsewhere in the system lead to irreconcilable differences. They are energetically costly (the energy density in their core reaches that of the prior symmetric vacuum) but topologically stable (the whole manifold would have to be rearranged to get rid of the defect). We show how, in a paradigmatic model of a quantum phase transition, a topological defect can be put in a non-local superposition, so that - in a region large compared to the size of its core - the order parameter of the system is "undecided" by being in a quantum superposition of conflicting choices of the broken symmetry. We demonstrate how to exhibit such a "Schrödinger kink" by devising a version of a double-slit experiment suitable for topological defects. Coherence detectable in such experiments will be suppressed as a consequence of interaction with the environment. We analyze environment-induced decoherence and discuss its role in symmetry breaking.

preprint2011arXiv

A ring of BEC pools as a trap for persistent flow

Mott insulator - superfluid transition in a periodic lattice of Josephson junctions can be driven by tunneling rate increase. Resulting winding numbers $W$ of the condensate wavefunction decrease with increasing quench time in accord with the Kibble-Zurek mechanism (KZM). However, in very slow quenches Bose-Hubbard dynamics rearranges wavefunction phase so that its random walk cools, $\bar{W^2}$ decreases and eventually the wavefunction becomes too cold to overcome potential barriers separating different $W$. Thus, in contrast with KZM, in very slow quenches $\bar{W^2}$ is set by random walk with "critical" step size, independently of $τ_Q$.

preprint2011arXiv

Excitation energy after a smooth quench in a Luttinger liquid

Low energy physics of quasi-one-dimensional ultracold atomic gases is often described by a gapless Luttinger liquid (LL). It is nowadays routine to manipulate these systems by changing their parameters in time but, no matter how slow the manipulation is, it must excite a gapless system. We study a smooth change of parameters of the LL (a smooth "quench") with a variable quench time and find that the excitation energy decays with an inverse power of the quench time. This universal exponent is -2 at zero temperature, and -1 for slow enough quenches at finite temperature. The smooth quench does not excite beyond the range of validity of the low energy LL description.

preprint2010arXiv

Adiabatic dynamics of an inhomogeneous quantum phase transition: the case of z > 1 dynamical exponent

We consider an inhomogeneous quantum phase transition across a multicritical point of the XY quantum spin chain. This is an example of a Lifshitz transition with a dynamical exponent z = 2. Just like in the case z = 1 considered in New J. Phys. 12, 055007 (2010) when a critical front propagates much faster than the maximal group velocity of quasiparticles vq, then the transition is effectively homogeneous: density of excitations obeys a generalized Kibble-Zurek mechanism and scales with the sixth root of the transition rate. However, unlike for z = 1, the inhomogeneous transition becomes adiabatic not below vq but a lower threshold velocity v', proportional to inhomogeneity of the transition, where the excitations are suppressed exponentially. Interestingly, the adiabatic threshold v' is nonzero despite vanishing minimal group velocity of low energy quasiparticles. In the adiabatic regime below v' the inhomogeneous transition can be used for efficient adiabatic quantum state preparation in a quantum simulator: the time required for the critical front to sweep across a chain of N spins adiabatically is merely linear in N, while the corresponding time for a homogeneous transition across the multicritical point scales with the sixth power of N. What is more, excitations after the adiabatic inhomogeneous transition, if any, are brushed away by the critical front to the end of the spin chain.

preprint2010arXiv

Disorder-Induced Order in Quantum XY Chains

We observe signatures of disorder-induced order in 1D XY spin chains with an external, site-dependent uni-axial random field within the XY plane. We numerically investigate signatures of a quantum phase transition at T=0, in particular an upsurge of the magnetization in the direction orthogonal to the external magnetic field, and the scaling of the block-entropy with the amplitude of this field. Also, we discuss possible realizations of this effect in ultra-cold atom experiments.

preprint2010arXiv

Dynamics of a Quantum Phase Transition and Relaxation to a Steady State

We review recent theoretical work on two closely related issues: excitation of an isolated quantum condensed matter system driven adiabatically across a continuous quantum phase transition or a gapless phase, and apparent relaxation of an excited system after a sudden quench of a parameter in its Hamiltonian. Accordingly the review is divided into two parts. The first part revolves around a quantum version of the Kibble-Zurek mechanism including also phenomena that go beyond this simple paradigm. What they have in common is that excitation of a gapless many-body system scales with a power of the driving rate. The second part attempts a systematic presentation of recent results and conjectures on apparent relaxation of a pure state of an isolated quantum many-body system after its excitation by a sudden quench. This research is motivated in part by recent experimental developments in the physics of ultracold atoms with potential applications in the adiabatic quantum state preparation and quantum computation.

preprint2010arXiv

Dynamics of an inhomogeneous quantum phase transition

We argue that in a second order quantum phase transition driven by an inhomogeneous quench density of quasiparticle excitations is suppressed when velocity at which a critical point propagates across a system falls below a threshold velocity equal to the Kibble-Zurek correlation length times the energy gap at freeze-out divided by $\hbar$. This general prediction is supported by an analytic solution in the quantum Ising chain. Our results suggest, in particular, that adiabatic quantum computers can be made more adiabatic when operated in an "inhomogeneous" way.

preprint2010arXiv

Spontaneous symmetry breaking in a generalized orbital compass model

We introduce a generalized two-dimensional orbital compass model, which interpolates continuously from the classical Ising model to the orbital compass model with frustrated quantum interactions, and investigate it using the multiscale entanglement renormalization ansatz (MERA). The results demonstrate that increasing frustration of exchange interactions triggers a second order quantum phase transition to a degenerate symmetry broken state which minimizes one of the interactions in the orbital compass model. Using boson expansion within the spin-wave theory we unravel the physical mechanism of the symmetry breaking transition as promoted by weak quantum fluctuations and explain why this transition occurs only surprisingly close to the maximally frustrated interactions of the orbital compass model. The spin waves remain gapful at the critical point, and both the boson expansion and MERA do not find any algebraically decaying spin-spin correlations in the critical ground state.

preprint2008arXiv

Dynamics of a quantum phase transition with decoherence: the quantum Ising chain in a static spin environment

We consider a linear quench from the paramagnetic to ferromagnetic phase in the quantum Ising chain interacting with a static spin environment. Both decoherence from the environment and non-adiabaticity of the evolution near a critical point excite the system from the final ferromagnetic ground state. For weak decoherence and relatively fast quenches the excitation energy, proportional to the number of kinks in the final state, decays like an inverse square root of a quench time, but slow transitions or strong decoherence make it decay in a much slower logarithmic way. We also find that fidelity between the final ferromagnetic ground state and a final state after a quench decays exponentially with a size of a chain, with a decay rate proportional to average density of excited kinks, and a proportionality factor evolving from 1.3 for weak decoherence and fast quenches to approximately 1 for slow transitions or strong decoherence. Simultaneously, correlations between kinks randomly distributed along the chain evolve from a near-crystalline anti-bunching to a Poissonian distribution of kinks in a number of isolated Anderson localization centers randomly scattered along the chain.

preprint2008arXiv

Multi-scale Entanglement Renormalization Ansatz in Two Dimensions: Quantum Ising Model

We propose a symmetric version of the multi-scale entanglement renormalization Ansatz (MERA) in two spatial dimensions (2D) and use this Ansatz to find an unknown ground state of a 2D quantum system. Results in the simple 2D quantum Ising model on the $8\times8$ square lattice are found to be very accurate even with the smallest non-trivial truncation parameter.

preprint2006arXiv

Dynamics of a quantum phase transition in the random Ising model

A quantum phase transition from paramagnetic to ferromagnetic phase is driven by a time-dependent external magnetic field. For any rate of the transition the evolution is non-adiabatic and finite density of defects is excited in the ferromagnetic state. The density of excitations has only logarithmic dependence on the transition rate. This is much weaker than any usual power law scaling predicted for pure systems by the Kibble-Zurek mechanism.

preprint2005arXiv

Images of a Bose-Einstein condensate: diagonal dynamical Bogoliubov vacuum

Evolution of a Bose-Einstein condensate subject to a time-dependent external perturbation can be described by a time-dependent Bogoliubov theory: a condensate initially in its ground state Bogoliubov vacuum evolves into a time-dependent excited state which can be formally written as a time-dependent Bogoliubov vacuum annihilated by time-dependent quasiparticle annihilation operators. We prove that any Bogoliubov vacuum can be brought to a diagonal form in a time-dependent orthonormal basis. This diagonal form is tailored for simulations of quantum measurements on an excited condensate. As an example we work out phase imprinting of a dark soliton followed by a density measurement.

preprint2004arXiv

Quantum fluctuations of a vortex in a dilute Bose-Einstein condensate

A vortex in a quasi two-dimensional Bose-Einstein condensate is subject to the Magnus force and can be effectively described as a planar particle in a uniform magnetic field. Quantization of this effective particle leads to the lowest Landau level where the most localized wave function is a gaussian. In this gaussian state vortex position seems to fluctuate with an average magnitude set by the magnetic width of the gaussian. We readdress this problem using the number-conserving version of the Bogoliubov theory. We find that the Bogoliubov mode that might be interpreted as a fluctuation of vortex position actually does not contribute to position fluctuation at all. The only non-zero contribution comes from phonons but it is an order of magnitude less than the simple estimate, based on the magnetic width of the effective gaussian wave packet.

preprint1994arXiv

Vortices and hierarchy of states in double-layer fractional Hall effect

Vortex solutions in $U(1)\times U(1)$ Chern-Simons theory coupled to a pair of hard-core bosons representing two layers of electrons are analysed. It is shown that there is such a range of parameters $(αβ<γ^{2})$ in which finite charge $(1,0)$ or $(0,1)$ vortices can not exist. Instead the minimal flux quanta are $(1,1)$ hybrids. Their statistical dominance is shown to have an influence on Haldane-Halderin hierarchy of states. The form of Coulomb potential and its influence on the range of physically acceptable states is discussed.