Dynamical Correlation Functions Research Articles

Space-like asymptotics of the thermal two-point functions of the XXZ spin-1/2 chain

This work proposes a closed formula for the leading term of the large-distance and long-time asymptotics in a cone of the space-like regime for the transverse dynamical two-point functions of the XXZ spin 1/2 chain at finite temperatures. The result follows from a simple analysis of the thermal form factor series for dynamical correlation functions. The obtained leading asymptotics are driven by the Bethe Ansatz data associated with the first sub-leading eigenvalue of the quantum transfer matrix.

Robust Measurements of n-Point Correlation Functions of Driven-Dissipative Quantum Systems on a Digital Quantum Computer.

We propose and demonstrate a unified hierarchical method to measure n-point correlation functions that can be applied to driven, dissipative, or otherwise open or nonequilibrium quantum systems. In this method, the time evolution of the system is repeatedly interrupted by interacting an ancilla qubit with the system through a controlled operation, and measuring the ancilla immediately afterward. We discuss the robustness of this method as compared to other ancilla-based interferometric techniques (such as the Hadamard test), and highlight its advantages for near-term quantum simulations of open quantum systems. We implement the method on a quantum computer in order to measure single-particle Green's functions of a driven-dissipative fermionic system. This Letter shows that dynamical correlation functions for driven-dissipative systems can be robustly measured with near-term quantum computers.

Dual symplectic classical circuits: An exactly solvable model of many-body chaos

We propose a general exact method of calculating dynamical correlation functions in dual symplectic brick-wall circuits in one dimension. These are deterministic classical many-body dynamical systems which can be interpreted in terms of symplectic dynamics in two orthogonal (time and space) directions. In close analogy with quantum dual-unitary circuits, we prove that two-point dynamical correlation functions are non-vanishing only along the edges of the light cones. The dynamical correlations are exactly computable in terms of a one-site Markov transfer operator, which is generally of infinite dimensionality. We test our theory in a specific family of dual-symplectic circuits, describing the dynamics of a classical Floquet spin chain. Remarkably, expressing these models in the form of a composition of rotations leads to a transfer operator with a block diagonal form in the basis of spherical harmonics. This allows us to obtain analytical predictions for simple local observables. We demonstrate the validity of our theory by comparison with Monte Carlo simulations, displaying excellent agreement with the latter for a choice of observables.

Scale-Free Chaos in the 2D Harmonically Confined Vicsek Model.

Animal motion and flocking are ubiquitous nonequilibrium phenomena that are often studied within active matter. In examples such as insect swarms, macroscopic quantities exhibit power laws with measurable critical exponents and ideas from phase transitions and statistical mechanics have been explored to explain them. The widely used Vicsek model with periodic boundary conditions has an ordering phase transition but the corresponding homogeneous ordered or disordered phases are different from observations of natural swarms. If a harmonic potential (instead of a periodic box) is used to confine particles, then the numerical simulations of the Vicsek model display periodic, quasiperiodic, and chaotic attractors. The latter are scale-free on critical curves that produce power laws and critical exponents. Here, we investigate the scale-free chaos phase transition in two space dimensions. We show that the shape of the chaotic swarm on the critical curve reflects the split between the core and the vapor of insects observed in midge swarms and that the dynamic correlation function collapses only for a finite interval of small scaled times. We explain the algorithms used to calculate the largest Lyapunov exponents, the static and dynamic critical exponents, and compare them to those of the three-dimensional model.

Superdiffusion from Nonabelian Symmetries in Nearly Integrable Systems

The Heisenberg spin chain is a canonical integrable model. As such, it features stable ballistically propagating quasiparticles, but spin transport is subballistic at any nonzero temperature: An initially localized spin fluctuation spreads in time t to a width t2/3. This exponent as well as the functional form of the dynamical spin correlation function suggest that spin transport is in the Kardar–Parisi–Zhang (KPZ) universality class. However, the full counting statistics of magnetization is manifestly incompatible with KPZ scaling. A simple two-mode hydrodynamic description, derivable from microscopic principles, captures both the KPZ scaling of the correlation function and the coarse features of the full counting statistics, but remains to be numerically validated. These results generalize to any integrable spin chain invariant under a continuous nonabelian symmetry and are surprisingly robust against moderately strong integrability-breaking perturbations that respect the nonabelian symmetry.

Dynamical correlations in simple disorder and complex disorder liquids

Liquids in equilibrium exhibit two types of disorder, simple and complex. Typical simple disorder liquid is liquid nitrogen, or weakly polar liquids. Complex liquids concern those who can form long lived local assemblies, and cover a large range from water to some soft matter and biological liquids. The existence of such structures leaves characteric features upon the atom-atom correlation functions, concerning both atoms which directly participate to these structures and those who do not. The question we ask here is: do these features have also characteristic dynamical aspects, which could be tracked through dynamical correlation functions? Herein, we compare the van Hove function, intermediate scattering function and the dynamical structure factor, for both types of liquids, using force field models and computer simulations. The calculations reveal the paradoxical fact that neighbouring atom correlations for simple disorder liquids relax slower than that for complex disorder liquids, while prepeak features typical of complex disorder liquids relax even slower. This is an indication of the existence of fast kinetic self-assembly processes in complex disorder liquids, while the lifetime of such assemblies itself is quite slow. This is further confirmed by the existence of a very low-k dynamical pre-peak uncovered in the case of water and ethanol.

Extending the thermodynamic form factor bootstrap program: multiple particle-hole excitations, crossing symmetry, and reparameterization invariance

In this study, we further the thermodynamic bootstrap program which involves a set of recently developed ideas used to determine thermodynamic form factors of local operators in integrable quantum field theories. These form factors are essential building blocks for dynamic correlation functions at finite temperatures or non-equilibrium stationary states. In this work we extend this program in three ways. Firstly, we demonstrate that the conjectured annihilation pole axiom is valid in the low energy particle-hole excitations. Secondly, we introduce a crossing relation, which establishes a connection between form factors with different excitation content. Typically, the crossing relation is a consequence of Lorentz invariance, but due to the finite energy density of the considered states, Lorentz invariance is broken. Nonetheless a crossing relation involving excitations with both particles and holes can established using the finite volume representation of the thermodynamic form factors. Finally, we demonstrate that the thermodynamic form factors satisfy a reparameterization invariance, an invariance which encompasses crossing. Reparameterization invariance exploits the fact that the details of the representation of the thermodynamic state are unimportant. In the course of developing these results, we demonstrate the internal consistency of the thermodynamic form factor bootstrap program in a number of ways. Finally, we provide explicit computations of form factors of conserved charges and densities with crossed excitations and show our results can be used to infer information about thermodynamic form factors in the Lieb-Liniger model.

Integrable crosscap states: from spin chains to 1D Bose gas

The notion of a crosscap state, a special conformal boundary state first defined in 2d CFT, was recently generalized to 2d massive integrable quantum field theories and integrable spin chains. It has been shown that the crosscap states preserve integrability. In this work, we first generalize this notion to the Lieb-Liniger model, which is a prototype of integrable non-relativistic many-body systems. We then show that the defined crosscap state preserves integrability. We derive the exact overlap formula of the crosscap state and the on-shell Bethe states. As a byproduct, we prove the conjectured overlap formula for integrable spin chains rigorously by coordinate Bethe ansatz. It turns out that the overlap formula for both models take the same form as a ratio of Gaudin-like determinants with a trivial prefactor. Finally we study quench dynamics of the crosscap state, which turns out to be surprisingly simple. The stationary density distribution is simply a constant. We also derive the analytic formula for dynamical correlation functions in the Tonks-Girardeau limit.

Exact results of dynamical structure factor of Lieb–Liniger model

The dynamical structure factor (DSF) represents a measure of dynamical density–density correlations in a quantum many-body system. Due to the complexity of many-body correlations and quantum fluctuations in a system of an infinitely large Hilbert space, such kind of dynamical correlations often impose a big theoretical challenge. For one-dimensional (1D) quantum many-body systems, qualitative predictions of dynamical response functions are usually carried out by using the Tomonaga–Luttinger liquid (TLL) theory. In this scenario, a precise evaluation of the DSF for a 1D quantum system with arbitrary interaction strength remains a formidable task. In this paper, we use the form factor approach based on algebraic Bethe ansatz theory to calculate precisely the DSF of Lieb–Liniger model with an arbitrary interaction strength at a large scale of particle number. We find that the DSF for a system as large as 2000 particles enables us to depict precisely its line-shape from which the power-law singularity with corresponding exponents in the vicinities of spectral thresholds naturally emerge. It should be noted that, the advantage of our algorithm promises an access to the threshold behavior of dynamical correlation functions, further confirming the validity of nonlinear TLL theory besides Kitanine et al (2012 J. Stat. Mech. P09001). Finally we discuss a comparison of results with the results from the ABACUS method by J-S Caux (2009 J. Math. Phys. 50 095214) as well as from the strongly coupling expansion by Brand and Cherny (2005 Phys. Rev. A 72 033619).

A high-performance implementation of atomistic spin dynamics simulations on x86 CPUs

Atomistic spin dynamics simulations provide valuable information about the energy spectrum of magnetic materials in different phases, allowing one to identify instabilities and the nature of their excitations. However, the time cost of evaluating the spin-spin correlation functions in real space increases quadratically as the number of spins N, leading to significant computational effort, making the simulation of large spin systems very challenging. In this work, we propose to use a highly optimized general matrix multiply (GEMM) subroutine to calculate the dynamical spin-spin correlation function that can achieve near-optimal hardware utilization. Furthermore, we fuse the element-wise operations in the calculation of S(q,t) into the in-house GEMM kernel, which results in further performance improvements of 44% - 71% on several relatively large lattice sizes when compared to the implementation that uses the GEMM subroutine in OpenBLAS, which is the state-of-the-art open source library for Basic Linear Algebra Subroutine (BLAS).

Self-similarity of spectral response functions for fractional quantum Hall states

Spectral response functions are central quantities in the analysis of quantum many-body states, since they describe the response of many-body systems to external perturbations and hence directly correspond to observables in experiments. In this paper, we evaluate a momentum-averaged dynamical density structure factor for the fermionic ν = 1 / 3 fractional quantum Hall (FQH) state on a torus, using the continued fraction method to compute the dynamical correlation function. We establish the scaling behaviour of the screened Coulomb structure factor with respect to interaction range, and expose an inherent self-similarity of structure factors in the frequency domain. These results highlight the statistical properties of spectral response functions for FQH states and show how they can be efficiently approximated in numerical models.

Lotka-Volterra predator-prey model with periodically varying carrying capacity.

We study the stochastic spatial Lotka-Volterra model for predator-prey interaction subject to a periodically varying carrying capacity. The Lotka-Volterra model with on-site lattice occupation restrictions (i.e., finite local carrying capacity) that represent finite food resources for the prey population exhibits a continuous active-to-absorbing phase transition. The active phase is sustained by the existence of spatiotemporal patterns in the form of pursuit and evasion waves. Monte Carlo simulations on a two-dimensional lattice are utilized to investigate the effect of seasonal variations of the environment on species coexistence. The results of our simulations are also compared to a mean-field analysis in order to specifically delineate the impact of stochastic fluctuations and spatial correlations. We find that the parameter region of predator and prey coexistence is enlarged relative to the stationary situation when the carrying capacity varies periodically. The (quasi-)stationary regime of our periodically varying Lotka-Volterra predator-prey system shows qualitative agreement between the stochastic model and the mean-field approximation. However, under periodic carrying capacity-switching environments, the mean-field rate equationspredict period-doubling scenarios that are washed out by internal reaction noise in the stochastic lattice model. Utilizing visual representations of the lattice simulations and dynamical correlation functions, we study how the pursuit and evasion waves are affected by ensuing resonance effects. Correlation function measurements indicate a time delay in the response of the system to sudden changes in the environment. Resonance features are observed in our simulations that cause prolonged persistent spatial correlations. Different effective static environments are explored in the extreme limits of fast and slow periodic switching. The analysis of the mean-field equationsin the fast-switching regime enables a semiquantitative description of the (quasi-)stationary state.

Role of Charge Ordering in the Dynamics of Cluster Formation in Associated Liquids.

Liquids are archetypes of disordered systems, yet liquids of polar molecules are locally more ordered than nonpolar molecules, due to the Coulomb interaction based charge ordering phenomenon. Hydrogen bonded liquids, such as water or alcohols, for example, represent a special type of polar liquids, in that they form labile clustered local structures. For water, in particular, hydrogen bonding and the related local tetrahedrality, play an important role in the various attempts to understand this liquid. However, labile structures imply dynamics, and it is not clear how it affects the understanding of this type of liquids from purely static point of view. Herein, we propose to reconsider hydrogen bonding as a charge ordering process. This concept allows us to demonstrate the insufficiency of the analysis of the microscopic structure based solely on static pair correlation functions, and the need for dynamical correlation functions, both in real and reciprocal space. The subsequent analysis allows to recover several aspects of our understanding of hydrogen bonded liquids, but from the charge order perspective. For water, it confirms the jump rotation picture found recently, and it allows to rationalize the contradicting pictures that arise when following the interpretations based on hydrogen bonding. For alcohols, it allows to understand the dynamical origin of the scattering prepeak, which does not exist for water, despite the fact that both these liquids have very similar hydroxyl group chain clusters. The concept of charge ordering complemented by the analysis of dynamical correlation functions appear as a promising way to understand microheterogeneity in complex liquids and mixtures from kinetics point of view.

Dynamics of active particles with translational and rotational inertia

Inertial effects affecting both the translational and rotational dynamics are inherent to a broad range of active systems at the macroscopic scale. Thus, there is a pivotal need for proper models in the framework of active matter to correctly reproduce experimental results, hopefully achieving theoretical insights. For this purpose, we propose an inertial version of the active Ornstein–Uhlenbeck particle (AOUP) model accounting for particle mass (translational inertia) as well as its moment of inertia (rotational inertia) and derive the full expression for its steady-state properties. The inertial AOUP dynamics introduced in this paper is designed to capture the basic features of the well-established inertial active Brownian particle model, i.e. the persistence time of the active motion and the long-time diffusion coefficient. For a small or moderate rotational inertia, these two models predict similar dynamics at all timescales and, in general, our inertial AOUP model consistently yields the same trend upon changing the moment of inertia for various dynamical correlation functions.

Ternary Unitary Quantum Lattice Models and Circuits in 2+1 Dimensions.

We extend the concept of dual unitary quantum gates introduced in Phys. Rev. Lett. 123, 210601 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.210601 to quantum lattice models in 2+1 dimensions, by introducing and studying ternary unitary four-particle gates, which are unitary in time and both spatial dimensions. When used as building blocks of lattice models with periodic boundary conditions in time and space (corresponding to infinite temperature states), dynamical correlation functions exhibit a light ray structure. We also generalize solvable matrix product states introduced in Phys. Rev. B 101, 094304 (2020)PRBMDO2469-995010.1103/PhysRevB.101.094304 to two spatial dimensions with cylindrical boundary conditions, by showing that the analogous solvable projected entangled pair states can be identified with matrix product unitaries. In the resulting tensor network for evaluating equal-time correlation functions, the bulk ternary unitary gates cancel out. We delineate and implement a numerical algorithm for computing such correlations by contracting the remaining tensors.

Schwinger-Dyson Equations-of-Motion Method for Quantum Spin-1/2

The essential difficulty of a large portion of the strongly interacting systems can be traced down to the “noncanonicality” of the underlying operators which are used to formulate the theory, like t-J model, Hubbard-type models in the large-U limit, quantum spin, etc. We proposed a new Schwinger-Dyson equations-of-motion method to compute the dynamical correlation functions of quantum spins at zero temperature. This method can also compute all the static observables from the dynamical correlation functions. We propose that this method can formulate dynamical perturbation calculations for quantum spin-1/2 models.

Scale-free chaos in the confined Vicsek flocking model.

The Vicsek model encompasses the paradigm of active dry matter. Motivated by collective behavior of insects in swarms, we have studied finite-size effects and criticality in the three-dimensional, harmonically confined Vicsek model. We have discovered a phase transition that exists for appropriate noise and small confinement strength. On the critical line of confinement versus noise, swarms are in a state of scale-free chaos characterized by minimal correlation time, correlation length proportional to swarm size and topological data analysis. The critical line separates dispersed single clusters from confined multicluster swarms. Scale-free chaotic swarms occupy a compact region of space and comprise a recognizable "condensed" nucleus and particles leaving and entering it. Susceptibility, correlation length, dynamic correlation function, and largest Lyapunov exponent obey power laws. The critical line and a narrow criticality region close to it move simultaneously to zero confinement strength for infinitely many particles. At the end of the first chaotic window of confinement, there is another phase transition to infinitely dense clusters of finite size that may be termed flocking black holes.

Nearly Optimal Quantum Algorithm for Estimating Multiple Expectation Values.

Many quantum algorithms involve the evaluation of expectation values. Optimal strategies for estimating a single expectation value are known, requiring a number of state preparations that scales with the target error ϵ as O(1/ϵ). In this Letter, we address the task of estimating the expectation values of M different observables, each to within additive error ϵ, with the same 1/ϵ dependence. We describe an approach that leverages Gilyén etal.'s quantum gradient estimation algorithm to achieve O(sqrt[M]/ϵ) scaling up to logarithmic factors, regardless of the commutation properties of the M observables. We prove that this scaling is worst-case optimal in the high-precision regime if the state preparation is treated as a black box, even when the operators are mutually commuting. We highlight the flexibility of our approach by presenting several generalizations, including a strategy for accelerating the estimation of a collection of dynamic correlation functions.

Efficient simulation of dynamics in two-dimensional quantum spin systems with isometric tensor networks

Simulating the dynamics of finite two-dimensional (2D) quantum many-body systems poses challenges to numerical methods due to the exponential growth of the many-body Hilbert space. In this work, the authors explore the variational power of recently introduced isometric tensor network states (isoTNS) --- a subset of general tensor-network states that allow an efficient numerical simulation. As a efficient tool to find ground states, the variational density matrix renormalization group (DMRG${}_{2}$) is introduced. It is then demonstrated that dynamical correlation functions for the 2D transverse field Ising model and the Kitaev Honeycomb model can be obtained.

Finite-temperature dynamics in gapped one-dimensional models in the sine-Gordon family

The sine-Gordon model appears as the low-energy effective field theory of various one-dimensional gapped quantum systems. Here we investigate the dynamics of generic, nonintegrable systems belonging to the sine-Gordon family at finite temperature within the semiclassical approach. Focusing on timescales where the effect of nontrivial quasiparticle scatterings becomes relevant, we obtain universal results for the long-time behavior of dynamical correlation functions. We find that correlation functions of vertex operators behave neither ballistically nor diffusively but follow a stretched exponential decay in time. We also study the full counting statistics of the topological current and find that distribution of the transferred charge is non-Gaussian with its cumulants scaling nonuniformly in time.