Source author record

Marek M. Rams

Marek M. Rams 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

19works
8topics
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

19 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.

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

Breaking the entanglement barrier: Tensor network simulation of quantum transport

The recognition that large classes of quantum many-body systems have limited entanglement in the ground and low-lying excited states led to dramatic advances in their numerical simulation via so-called tensor networks. However, global dynamics elevates many particles into excited states, and can lead to macroscopic entanglement and the failure of tensor networks. Here, we show that for quantum transport -- one of the most important cases of this failure -- the fundamental issue is the canonical basis in which the scenario is cast: When particles flow through an interface, they scatter, generating a "bit" of entanglement between spatial regions with each event. The frequency basis naturally captures that -- in the long-time limit and in the absence of inelastic scattering -- particles tend to flow from a state with one frequency to a state of identical frequency. Recognizing this natural structure yields a striking -- potentially exponential in some cases -- increase in simulation efficiency, greatly extending the attainable spatial- and time-scales, and broadening the scope of tensor network simulation to hitherto inaccessible classes of non-equilibrium many-body problems.

preprint2020arXiv

Open System Tensor Networks and Kramers' Crossover for Quantum Transport

Tensor networks are a powerful tool for many-body ground states with limited entanglement. These methods can nonetheless fail for certain time-dependent processes - such as quantum transport or quenches - where entanglement growth is linear in time. Matrix-product-state decompositions of the resulting out-of-equilibrium states require a bond dimension that grows exponentially, imposing a hard limit on simulation timescales. However, in the case of transport, if the reservoir modes of a closed system are arranged according to their scattering structure, the entanglement growth can be made logarithmic. Here, we apply this ansatz to open systems via extended reservoirs that have explicit relaxation. This enables transport calculations that can access steady states, time dynamics and noise, and periodic driving (e.g., Floquet states). We demonstrate the approach by calculating the transport characteristics of an open, interacting system. These results open a path to scalable and numerically systematic many-body transport calculations with tensor networks.

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.

preprint2016arXiv

Inhomogeneous quasi-adiabatic driving of quantum critical dynamics in weakly disordered spin chains

We introduce an inhomogeneous protocol to drive a weakly disordered quantum spin chain quasi-adiabatically across a quantum phase transition and minimize the residual energy of the final state. The number of spins that simultaneously reach the critical point is controlled by the length scale in which the magnetic field is modulated, introducing an effective size that favors adiabatic dynamics. The dependence of the residual energy on this length scale and the velocity at which the magnetic field sweeps out the chain is shown to be nonmonotonic. We determine the conditions for an optimal suppression of the residual energy of the final state and show that inhomogeneous driving can outperform conventional adiabatic schemes based on homogeneous control fields by several orders of magnitude.

preprint2016arXiv

Matrix product state renormalization

The truncation or compression of the spectrum of Schmidt values is inherent to the matrix product state (MPS) approximation of one-dimensional quantum ground states. We provide a renormalization group picture by interpreting this compression as an application of Wilson's numerical renormalization group along the imaginary time direction appearing in the path integral representation of the state. The location of the physical index is considered as an impurity in the transfer matrix and static MPS correlation functions are reinterpreted as dynamical impurity correlations. Coarse-graining the transfer matrix is performed using a hybrid variational ansatz based on matrix product operators, combining ideas of MPS and the multi-scale entanglement renormalization ansatz. Through numerical comparison with conventional MPS algorithms, we explicitly verify the impurity interpretation of MPS compression, as put forward by [V. Zauner et al., New J. Phys. 17, 053002 (2015)] for the transverse-field Ising model. Additionally, we motivate the conceptual usefulness of endowing MPS with an internal layered structure by studying restricted variational subspaces to describe elementary excitations on top of the ground state, which serves to elucidate a transparent renormalization group structure ingrained in MPS descriptions of ground states.

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.

preprint2015arXiv

Truncating an exact Matrix Product State for the XY model: transfer matrix and its renormalisation

We discuss how to analytically obtain an -- essentially infinite -- Matrix Product State (MPS) representation of the ground state of the XY model. On the one hand this allows to illustrate how the Ornstein-Zernike form of the correlation function emerges in the exact case using standard MPS language. On the other hand we study the consequences of truncating the bond dimension of the exact MPS, which is also part of many tensor network algorithms, and analyze how the truncated MPS transfer matrix is representing the dominant part of the exact quantum transfer matrix. In the gapped phase we observe that the correlation length obtained from a truncated MPS approaches the exact value following a power law in effective bond dimension. In the gapless phase we find a good match between a state obtained numerically from standard MPS techniques with finite bond dimension, and a state obtained by effective finite imaginary time evolution in our framework. This provides a direct hint for a geometric interpretation of Finite Entanglement Scaling at the critical point in this case. Finally, by analyzing the spectra of transfer matrices, we support the interpretation put forward by [V. Zauner at. al., New J. Phys. 17, 053002 (2015)] that the MPS transfer matrix emerges from the quantum transfer matrix though the application of Wilson's Numerical Renormalisation Group along the imaginary-time direction.

preprint2013arXiv

Exact results for fidelity susceptibility of the quantum Ising model: The interplay between parity, system size, and magnetic field

We derive an exact closed-form expression for fidelity susceptibility of even- and odd-sized quantum Ising chains in the transverse field. To this aim, we diagonalize the Ising Hamiltonian and study the gap between its positive and negative parity subspaces. We derive an exact closed-form expression for the gap and use it to identify the parity of the ground state. We point out misunderstanding in some of the former studies of fidelity susceptibility and discuss its consequences. Last but not least, we rigorously analyze the properties of the gap. For example, we derive analytical expressions showing its exponential dependence on the ratio between the system size and the correlation length.

preprint2012arXiv

A quantum phase transition in a quantum external field: Superposing two magnetic phases

We study an Ising chain undergoing a quantum phase transition in a quantum magnetic field. Such a field can be emulated by coupling the chain to a central spin initially in a superposition state. We show that - by adiabatically driving such a system - one can prepare a quantum superposition of any two ground states of the Ising chain. In particular, one can end up with the Ising chain in a superposition of ferromagnetic and paramagnetic phases -- a scenario with no analogue in prior studies of quantum phase transitions. Remarkably, the resulting magnetization of the chain encodes the position of the critical point and universal critical exponents, as well as the ground state fidelity.

preprint2012arXiv

Assisted finite-rate adiabatic passage across a quantum critical point: Exact solution for the quantum Ising model

The dynamics of a quantum phase transition is inextricably woven with the formation of excitations, as a result of the critical slowing down in the neighborhood of the critical point. We design a transitionless quantum driving through a quantum critical point that allows one to access the ground state of the broken-symmetry phase by a finite-rate quench of the control parameter. The method is illustrated in the one-dimensional quantum Ising model in a transverse field. Driving through the critical point is assisted by an auxiliary Hamiltonian, for which the interplay between the range of the interaction and the modes where excitations are suppressed is elucidated.

preprint2011arXiv

Scaling of ground state fidelity in the thermodynamic limit: XY model and beyond

We study ground state fidelity defined as the overlap between two ground states of the same quantum system obtained for slightly different values of the parameters of its Hamiltonian. We focus on the thermodynamic regime of the XY model and the neighborhood of its critical points. We describe in detail cases when fidelity is dominated by the universal contribution reflecting quantum criticality of the phase transition. We show that proximity to the multicritical point leads to anomalous scaling of fidelity. We also discuss fidelity in a regime characterized by pronounced oscillations resulting from the change of either the system size or the parameters of the Hamiltonian. Moreover, we show when fidelity is dominated by non-universal contributions, study fidelity in the extended Ising model, and illustrate how our results provide additional insight into dynamics of quantum phase transitions. Special attention is put to studies of fidelity from the momentum space perspective. All our main results are obtained analytically. They are in excellent agreement with numerics.

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 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.

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.