Quantum Spin Chains Research Articles

Spectral crossovers and universality in quantum spin chains coupled to random fields

We study the spectral properties of and spectral crossovers between different random matrix ensembles [Poissonian, Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE)] in correlated spin-chain systems, in the presence of random magnetic fields, and the scalar spin-chirality term, competing with the usual isotropic and time-reversal invariant Heisenberg term. We have investigated these crossovers in the context of the level-spacing distribution and the level-spacing ratio distribution. We use random matrix theory (RMT) analytical results to fit the observed Poissonian-to-GOE and GOE-to-GUE crossovers, and examine the relationship between the RMT crossover parameter $\ensuremath{\lambda}$ and scaled physical parameters of the spin-chain systems in terms of a scaling exponent. We find that the crossover behavior exhibits universality, in the sense that, it becomes independent of lattice size in the large Hamiltonian matrix dimension limit. Moreover, the scaling exponent obtained from such a finite size scaling analysis, seems to be quite robust and independent of the type of crossover considered or the specific spectral correlation measure used.

A Lieb–Robinson bound for quantum spin chains with strong on-site impurities

We consider a quantum spin chain with nearest neighbor interactions and sparsely distributed on-site impurities. We prove commutator bounds for its Heisenberg dynamics which incorporate the coupling strengths of the impurities. The impurities are assumed to satisfy a minimum spacing, and each impurity has a non-degenerate spectrum. Our results are proven in a broadly applicable setting, both in finite volume and in thermodynamic limit. We apply our results to improve Lieb–Robinson bounds for the Heisenberg spin chain with a random, sparse transverse field drawn from a heavy-tailed distribution.

Exact ground state and elementary excitations of a competing spin chain with twisted boundary condition

A novel Bethe ansatz scheme is proposed to investigate the exact physical properties of an integrable anisotropic quantum spin chain with competing interactions among the nearest, next nearest neighbor and chiral three spin couplings, where the boundary condition is the twisted one. The eigenvalue of the transfer matrix is characterized by its zero roots instead of the traditional Bethe roots. Based on the exact solution, the conserved momentum and charge operators of this U(1)-symmetry broken system are obtained. The ground state energy and density of rapidities are calculated. It is found that there exist three kinds of elementary excitations and the corresponding dispersion relations are obtained, which gives a different picture from that with periodic boundary condition.

Dynamical quantum phase transition in XY chains with the Dzyaloshinskii-Moriya and XZY–YZX three-site interactions

We study the dynamical quantum phase transitions (DQPTs) in the XY chains with the Dzyaloshinskii–Moriya interaction and the XZY–YZX type of three-site interaction after a sudden quench. Both the models can be mapped to the spinless free fermion models by the Jordan–Wigner and Bogoliubov transformations with the form , where the quasiparticle excitation spectra εk may be smaller than 0 for some k and are asymmetrical (εk ≠ ε –k ). It is found that the factors of Loschmidt echo equal 1 for some k corresponding to the quasiparticle excitation spectra of the pre-quench Hamiltonian satisfying εk ⋅ ε –k < 0, when the quench is from the gapless phase. By considering the quench from different ground states, we obtain the conditions for the occurrence of DQPTs for the general XY chains with gapless phase, and find that the DQPTs may not occur in the quench across the quantum phase transitions regardless of whether the quench is from the gapless phase to gapped phase or from the gapped phase to gapless phase. This is different from the DQPTs in the case of quench from the gapped phase to gapped phase, in which the DQPTs will always appear. Moreover, we analyze the different reasons for the absence of DQPTs in the quench from the gapless phase and the gapped phase. The conclusion can also be extended to the general quantum spin chains.

Localization crossover and subdiffusive transport in a classical facilitated network model of a disordered interacting quantum spin chain

We consider the random-field Heisenberg model, a paradigmatic model for many-body localization (MBL), and add a Markovian dephasing bath coupled to the Anderson orbitals of the model's non-interacting limit. We map this system to a classical facilitated hopping model that is computationally tractable for large system sizes, and investigate its dynamics. The classical model exhibits a robust crossover between an ergodic (thermal) phase and a frozen (localized) phase. The frozen phase is destabilized by thermal subregions (bubbles), which thermalize surrounding sites by providing a fluctuating interaction energy and so enable off-resonance particle transport. Investigating steady state transport, we observe that the interplay between thermal and frozen bubbles leads to a clear transition between diffusive and subdiffusive regimes. This phenomenology both describes the MBL system coupled to a bath, and provides a classical analogue for the many-body localization transition in the corresponding quantum model, in that the classical model displays long local memory times. It also highlights the importance of the details of the bath coupling in studies of MBL systems coupled to thermal environments.

Engineered dissipation induced entanglement transition in quantum spin chains: From logarithmic growth to area law

Recent theoretical work has shown that the competition between coherent unitary dynamics and stochastic measurements, performed by the environment, along wavefunction trajectories can give rise to transitions in the entanglement scaling. In this work, complementary to these previous studies, we analyze a situation where the role of Hamiltonian and dissipative dynamics is reversed. We consider an engineered dissipation, which stabilizes an entangled phase of a quantum spin$-1/2$ chain, while competing single-particle or interacting Hamiltonian dynamics induce a disentangled phase. Focusing on the single-particle unitary dynamics, we find that the system undergoes an entanglement transition from a logarithmic growth to an area law when the competition ratio between the unitary evolution and the non-unitary dynamics increases. We evidence that the transition manifests itself in state-dependent observables at a finite competition ratio for Hamiltonian and measurement dynamics. On the other hand, it is absent in trajectory-averaged steady-state dynamics, governed by a Lindblad master equation: although purely dissipative dynamics stabilizes an entangled state, for any non-vanishing Hamiltonian contribution the system ends up irremediably in a disordered phase. In addition, a single trajectory analysis reveals that the distribution of the entanglement entropy constitutes an efficient indicator of the transition. Complementarily, we explore the competition of the dissipation with coherent dynamics generated by an interacting Hamiltonian, and demonstrate that the entanglement transition also occurs in this second model. Our results suggest that this type of transition takes place for a broader class of Hamiltonians, underlining its robustness in monitored open quantum many-body systems.

False vacuum decay in quantum spin chains

The false vacuum decay has been a central theme in physics for half a century with applications to cosmology and to the theory of fundamental interactions. This fascinating phenomenon is even more intriguing when combined with the confinement of elementary particles. Due to the astronomical time scales involved, the research has so far focused on theoretical aspects of this decay. The purpose of this Letter is to show that the false vacuum decay is accessible to current optical experiments as quantum analog simulators of spin chains with confinement of the elementary excitations, which mimic the high energy phenomenology but in one spatial dimension. We study the non-equilibrium dynamics of the false vacuum in a quantum Ising chain and in an XXZ ladder. The false vacuum is the metastable state that arises in the ferromagnetic phase of the model when the symmetry is explicitly broken by a longitudinal field. This state decays through the formation of "bubbles" of true vacuum. Using iTEBD simulations, we are able to study the real-time evolution in the thermodynamic limit and measure the decay rate of local observables. We find that the numerical results agree with the theoretical prediction that the decay rate is exponentially small in the inverse of the longitudinal field.

Long Time Behaviour of a Local Perturbation in the Isotropic XY Chain Under Periodic Forcing

We study the isotropic XY quantum spin chain with a time-periodic transverse magnetic field acting on a single site. The asymptotic dynamics is described by a highly resonant Floquet–Schrödinger equation, for which we show the existence of a periodic solution if the forcing frequency is away from a discrete set of resonances. This in turn implies the state of the quantum spin chain to be asymptotically a periodic function synchronised with the forcing, also at arbitrarily low non-resonant frequencies. The behaviour at the resonances remains a challenging open problem.

Electric field- and strain-induced quantum phase transitions in a spin chain

The effect of the external electric field and/or the external strain on the low-temperature behavior of the quantum spin chain model is studied. The external electric field or the strain can cause the quantum magnetic structural phase transition between two magnetically ordered phases with different order parameters. Such a quantum critical point can be observed in the special behavior of the low-temperature specific heat: At the critical value the low-temperature specific heat manifests linear in temperature behavior instead of the exponentially small one for other values of the field. The magnetic susceptibility also manifests the special behavior, in which quantum phase transitions are revealed.

Weak three-dimensional coupling of Heisenberg quantum spin chains in SrCuTe2O6

The magnetic Hamiltonian of the Heisenberg quantum antiferromagnet SrCuTe$_{2}$O$_{6}$ is studied by inelastic neutron scattering technique on powder and single crystalline samples above and below the magnetic transition temperatures at 8 K and 2 K. The high temperature spectra reveal a characteristic diffuse scattering corresponding to a multi-spinon continuum confirming the dominant quantum spin-chain behavior due to the third neighbour interaction J$_{intra}$ = 4.22 meV (49 K). The low temperature spectra exhibits sharper excitations at energies below 1.25 meV which can be explained by considering a combination of weak antiferromagnetic first nearest neighbour interchain coupling J$_1$ = 0.17 meV (1.9 K) and even weaker ferromagnetic second nearest neighbour J$_2$ = -0.037 meV (-0.4 K) or a weak ferromagnetic J$_2$ = -0.11 meV (-1.3 K) and antiferromagnetic J$_6$ = 0.16 meV (1.85 K) giving rise to the long-range magnetic order and spin-wave excitations at low energies. These results suggest that SrCuTe$_{2}$O$_{6}$ is a highly one-dimensional Heisenberg system with three mutually perpendicular spin-chains coupled by a weak ferromagnetic J$_2$ in addition to the antiferromagnetic J$_1$ or J$_6$ presenting a contrasting scenario from the highly frustrated hyper-hyperkagome lattice (equally strong antiferromagnetic J$_1$ and J$_2$) found in the iso-structural PbCuTe$_{2}$O$_{6}$.

Intersecting defects in gauge theory, quantum spin chains, and Knizhnik-Zamolodchikov equations

We propose an interesting BPS/CFT correspondence playground: the correlation function of two intersecting half-BPS surface defects in four-dimensional mathcal{N} = 2 supersymmetric SU(N) gauge theory with 2N fundamental hypermultiplets. We show it satisfies a difference equation, the fractional quantum T-Q relation. Its Fourier transform is the 5-point conformal block of the {hat{mathfrak{sl}}}_N current algebra with one of the vertex operators corresponding to the N-dimensional {mathfrak{sl}}_N representation, which we demonstrate with the help of the Knizhnik-Zamolodchikov equation. We also identify the correlator with a state of the {XXX}_{{mathfrak{sl}}_2} spin chain of N Heisenberg-Weyl modules over Y ( {mathfrak{sl}}_2 ). We discuss the associated quantum Lax operators, and connections to isomonodromic deformations.

Simple mitigation of global depolarizing errors in quantum simulations.

To get the best possible results from current quantum devices error mitigation is essential. In this work we present a simple but effective error mitigation technique based on the assumption that noise in a deep quantum circuit is well described by global depolarizing error channels. By measuring the errors directly on the device, we use an error model ansatz to infer error-free results from noisy data. We highlight the effectiveness of our mitigation via two examples of recent interest in quantum many-body physics: entanglement measurements and real-time dynamics of confinement in quantum spin chains. Our technique enables us to get quantitative results from the IBM quantum computers showing signatures of confinement, i.e., we are able to extract the meson masses of the confined excitations which were previously out of reach. Additionally, we show the applicability of this mitigation protocol in a wider setting with numerical simulations of more general tasks using a realistic error model. Our protocol is device-independent, simply implementable, and leads to large improvements in results if the global errors are well described by depolarization.

Slow dynamics of the Fredkin spin chain

The dynamical behavior of a quantum many-particle system is characterized by the lifetime of its excitations. When the system is perturbed, observables of any non-conserved quantity decay exponentially, but those of a conserved quantity relax to equilibrium with a power law. Such processes are associated with a dynamical exponent $z$ that relates the spread of correlations in space and time. We present numerical results for the Fredkin model, a quantum spin chain with a three-body interaction term, which exhibits an unusually large dynamical exponent. We discuss our efforts to produce a reliable estimate $z$=3.16(1) through direct simulation of the quantum evolution and to explain the slow dynamics in terms of an excited bond that executes a constrained random walk in Monte Carlo time.

Efficient matrix product state methods for extracting spectral information on rings and cylinders

Based on the MPS formalism, we introduce an ansatz for capturing excited states in finite systems with open boundary conditions, providing a very efficient method for computing, e.g., the spectral gap of quantum spin chains. This method can be straightforwardly implemented on top of an existing DMRG or MPS ground-state code. Although this approach is built on open-boundary MPS, we also apply it to systems with periodic boundary conditions. Despite the explicit breaking of translation symmetry by the MPS representation, we show that momentum emerges as a good quantum number, and can be exploited for labeling excitations on top of MPS ground states. We apply our method to the critical Ising chain on a ring and the classical Potts model on a cylinder. Finally, we apply the same idea to compute excitation spectra for 2-D quantum systems on infinite cylinders. Again, despite the explicit breaking of translation symmetry in the periodic direction, we recover momentum as a good quantum number for labeling excitations. We apply this method to the 2-D transverse-field Ising model and the half-filled Hubbard model; for the latter, we obtain accurate results for, e.g., the hole dispersion for cylinder circumferences up to eight sites.

Fourth-neighbour two-point functions of the XXZ chain and the fermionic basis approach

We give a descriptive review of the fermionic basis approach to the theory of correlation functions of the XXZ quantum spin chain. The emphasis is on explicit formulae for short-range correlation functions which will be presented in a way that allows for their direct implementation on a computer. Within the fermionic basis approach a huge class of stationary reduced density matrices, compatible with the integrable structure of the model, assumes a factorized form. This means that all expectation values of local operators and all two-point functions, in particular, can be represented as multivariate polynomials in only two functions ρ and ω and their derivatives with coefficients that are rational in the deformation parameter q of the model. These coefficients are of ‘algebraic origin’. They do not depend on the choice of the density matrix, which only impacts the form of ρ and ω. As an example we work out in detail the case of the grand canonical ensemble at temperature T and magnetic field h for q in the critical regime. We compare our exact results for the fourth-neighbour two-point functions with asymptotic formulae for h, T = 0 and for finite h and T.

Exact Steady State of the Open XX -Spin Chain: Entanglement and Transport Properties

We study the reduced dynamics of open quantum spin chains of arbitrary length $N$ with nearest neighbour $XX$ interactions, immersed within an external constant magnetic field along the $z$ direction, whose end spins are weakly coupled to heat baths at different temperatures, via energy preserving couplings. We find the analytic expression of the unique stationary state of the master equation obtained in the so-called global approach based on the spectralization of the full chain Hamiltonian. Hinging upon the explicit stationary state, we reveal the presence of sink and source terms in the spin-flow continuity equation and compare their behaviour with that of the stationary heat flow. Moreover, we also obtain analytic expressions for the steady state two-spin reduced density matrices and for their concurrence. We then set up an algorithm suited to compute the stationary bipartite entanglement along the chain and to study its dependence on the Hamiltonian parameters and on the bath temperatures.

Constraint-Induced Delocalization.

We study the impact of quenched disorder on the dynamics of locally constrained quantum spin chains, that describe 1D arrays of Rydberg atoms in both the frozen (Ising-type) and dressed (XY-type) regime. Performing large-scale numerical experiments, we observe no trace of many-body localization even at large disorder. Analyzing the role of quenched disorder terms in constrained systems we show that they act in two, distinct and competing ways: as an on-site disorder term for the basic excitations of the system, and as an interaction between excitations. The two contributions are of the same order, and as they compete (one towards localization, the other against it), one does never enter a truly strong disorder, weak interaction limit, where many-body localization occurs. Such a mechanism is further clarified in the case of XY-type constrained models: there, a term which would represent a bona fide local quenched disorder term acting on the excitations of the clean model must be written as a series of nonlocal terms in the unconstrained variables. Our observations provide a simple picture to interpret the role of quenched disorder that could be immediately extended to other constrained models or quenched gauge theories.

Thermal Transistor Effect in Quantum Systems

We study a quantum system composed of three interacting qubits, each coupled to a different thermal reservoir. We show how to engineer it in order to build a quantum device that is analogous to an electronic bipolar transistor. We outline how the interaction among the qubits plays a crucial role for the appearance of the effect, also linking it to the characteristics of system-bath interactions that govern the decoherence and dissipation mechanism of the system. By comparing with previous proposals, the model considered here extends the regime of parameters where the transistor effect shows up and its robustness with respect to small variations of the coupling parameters. Moreover, our model appears to be more realistic and directly connected in terms of potential implementations to feasible setups in the domain of quantum spin chains and molecular nanomagnets.

Worldline algorithm by oracle-guided variational autoregressive network

The variational autoregressive network is extended to the Euclidean path integral representation of quantum partition function. An essential challenge is adapting the sequential process of sample generation by an autoregressive network to a nonlocal constraint due to the periodic boundary condition in path integral. An add-on oracle is devised for this purpose, which accurately identifies and stalls unviable configurations as soon as they occur. The oracle enables rejection-free sampling conforming to the periodic boundary condition. As a demonstration, the oracle-guided autoregressive network is applied to obtain variational solutions of quantum spin chains at finite temperatures with relatively large system sizes and numbers of time slicing, and to efficiently compute thermodynamic quantities.

Driven quantum spin chain in the presence of noise: Anti-Kibble-Zurek behavior

We study defect generation in a quantum XY-spin chain arising due to the linear drive of the many-body Hamiltonian in the presence of a time-dependent fast Gaussian noise. The main objective of this work is to quantify analytically the effects of noise on the defect density production. In the absence of noise, it is well known that in the slow sweep regime, the defect density follows the Kibble-Zurek (KZ) scaling behavior with respect to the sweep speed. We consider time-dependent fast Gaussian noise in the anisotropy of the spin-coupling term $[{\ensuremath{\gamma}}_{0}=({J}_{1}\ensuremath{-}{J}_{2})/({J}_{1}+{J}_{2})]$ and show via analytical calculations that the defect density exhibits anti-Kibble-Zurek (AKZ) scaling behavior in the slow sweep regime. In the limit of large chain length and long time, we calculate the entropy and magnetization density of the final decohered state and show that their scaling behavior is consistent with the AKZ picture in the slow sweep regime. We have also numerically calculated the sublattice spin correlators for finite separation by evaluating the Toeplitz determinants and find results consistent with the KZ picture in the absence of noise, while in the presence of noise and slow sweep speeds the correlators exhibit the AKZ behavior. Furthermore, by considering the large $n$-separation asymptotes of the Toeplitz determinants, we further quantify the effect of the noise on the spin-spin correlators in the final decohered state. We show that while the correlation length of the sublattice correlator scales according to the AKZ behavior, we obtain different scaling for the magnetization correlators.