BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jacek Dziarmaga

Jacek Dziarmaga contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
quant-phcond-mat.quant-gascond-mat.str-elcond-mat.mes-hallcond-mat.stat-mechcond-mat.dis-nn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#34;impulse&#34; where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such &#34;freeze-out&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;, proportional to inhomogeneity of the transition, where the excitations are suppressed exponentially. Interestingly, the adiabatic threshold v&#39; is nonzero despite vanishing minimal group velocity of low energy quasiparticles. In the adiabatic regime below v&#39; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access