Spin Correlation Functions Research Articles

Holographic Analysis of Boundary Correlation Functions for the Hyperbolic-Lattice Ising Model

Abstract We analyze boundary spin correlation functions of the hyperbolic-lattice Ising model from the holographic point of view. Using the corner-transfer-matrix renormalization-group (CTMRG) method, we demonstrate that the boundary correlation function exhibits power-law decay with quasiperiodic oscillation, while the bulk correlation function always decays exponentially. On the basis of the geometric relation between the bulk correlation path and distance along the outer edge boundary, we find that scaling dimensions for the boundary correlation function can be well explained by a combination of the bulk correlation length and background curvatures inherent to the hyperbolic lattice. We also investigate the cutoff effect of the bond dimension in CTMRG, revealing that the long-distance behavior of the boundary spin correlation is accurately described even with a small bond dimension. In contrast, the short-distance behavior rapidly loses its accuracy.

Four interacting spins: Addition of angular momenta, spin–spin correlation functions, and entanglement

We study four spins on a ring coupled through competing Heisenberg interactions between nearest neighbors, J, and next-nearest neighbors, J2≡αJ>0. From the pedagogical point of view, dealing with few spins illustrates how to add more than two angular momenta in a systematic way. The spectrum is obtained by using the rules for addition of four angular momenta, which allows us to follow the evolution of the ground state with α, characterized by level crossings and by spin–spin correlation functions. The reduced number of spins also allows us to illustrate how to quantify bipartite entanglement.

Twisted graphene superlattices: resonating valence bond states and magnetic properties

Resonating valence bond (RVB) states are fundamental for understanding quantum spin liquids in two-dimensional (2D) systems. The RVB state is a collective phenomenon in which spins are uncoupled. 2D lattices such as triangular, honeycomb, and dice lattices were investigated using the Hubbard model and exact diagonalization method. We analyzed the total spin, spin–spin correlation functions, local magnetic moments, and spin and charge gaps as a function of on-site Coulomb repulsion, electron concentration, and electronic hopping parameters. Phase diagrams showed that RVB states can live in half-filled and hole-doped anisotropic triangular lattices. We found two types of RVB states: one in the honeycomb sublattice and the other in the centered hexagons in the triangular lattices. Owing to the novel discovery of exotic magnetic ordering in triangular moiré patterns in twisted bilayer graphene and transition metal dichalcogenide systems, our results provide physical insights into the onset of magnetism and possible spin liquid states in these layered materials.

Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core

In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin–spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on a graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. We show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS which can improve the performance by 26% ∼ 33% compared to the implementation based on cuBLAS. Furthermore, we perform the ‘on-the-fly’ calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes large-scale atomistic spin dynamics simulations feasible and affordable.

Magnetic phases of XY model with three-spin terms: interplay of topology and entanglement

Magnetic and topological properties along with quantum correlations in terms of several entanglement measures have been investigated for an antiferromagnetic (AFM) spin-1/2 XY model in the presence of transverse magnetic field and XZX−YZY type of three-spin interactions. Symmetries of the spin Hamiltonian have been identified. Under the Jordan–Wigner transformation, the spin Hamiltonian converted into spinless superconducting model with nearest neighbor (NN) hopping and Cooper pairing terms in addition to next NN Cooper pairing potential. Long range AFM order has been studied in terms of staggered spin–spin correlation functions, while the topological orders have been characterized by winding numbers. Magnetic and topological phase diagrams have been prepared. Faithful coexistence of magnetic and topological superconducting phases is found in the entire parameter regime. Boundaries of various quantum phases have been marked and positions of bicritical points have been identified.

Coherent THz spin dynamics in antiferromagnets beyond the approximation of the Néel vector

Controlled generation of coherent spin waves with highest possible frequencies and shortest possible wavelengths is a cornerstone of spintronics and magnonics. Here, using Heisenberg antiferromagnet RbMnF3, we demonstrate that laser-induced THz spin dynamics corresponding to pairs of mutually coherent counter-propagating spin waves with the wavevectors up to the edge of the Brillouin zone cannot be understood in terms of magnetization and antiferromagnetic (Néel) vectors, conventionally used to describe spin waves. Instead, we propose to model such spin dynamics using the spin correlation function. We derive a quantum-mechanical equation of motion for the latter and emphasize that unlike the magnetization and antiferromagnetic vectors the spin correlations in antiferromagnets do not exhibit inertia.

Raman study of the spin–lattice excitations of the layered antiferromagnets MPSe3 (M = Fe, Mn)

Magnetism in van der Waals (vdW) materials has recently garnered significant attention due to their potential for advanced applications. In this study, we investigate the magnetic transition in layered MPSe3 (M = Fe, Mn) bulk flakes utilizing temperature-dependent Raman spectroscopy. Our findings reveal unique features in the Raman spectrum that are strongly correlated with the antiferromagnetic ordering in MPSe3. We analyze the spin-dependent Raman scattering modes in MPSe3 by applying the Suzuki–Kamimura formalism. Consequently, our comparative analyses of the spin correlation function and the reduced magnetization provide compelling evidence of the 2D Ising characteristic especially in FePSe3. Our results clearly show the potential of MPSe3 as a prominent platform for exploring fundamental aspects of low-dimensional magnetism.

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.

Monte Carlo simulation with Wolff algorithm for scaling behavior of two dimensional XY model with KT phase transition

The property of XY model with its KT phase transition has been discussed in many aspects such as specific heat and magnetic susceptibility along with correlation function. Also, this is one of the original models for topological phase transition which gives a new type of phase transition without symmetry breaking. Since the calculation of these properties above can be achieved by Monte Carlo Method with Wolff algorithm, the visualization of vortex and antivortex is also an interesting task which is accomplished in this article. By successfully modeling the system and getting the step-by-step configuration, the relation between the generation of vortex and the variation of temperature can be concluded in a qualitative way. The calculation of spin correlation function shows the phase transition between disorder and quasi long-range order which represents the topological phase transition. As a reliable method, the Monte Carlo method with Wolff algorithm can be a choice of generating the training set for neural network which can calculated much more sophisticated case in a relatively short time and the transition between classical simulation and quantum simulation can be approached as mentioned in the end of this article.

Yu-Shiba-Rusinov states assisted by asymmetric Coulomb repulsion in a bipartite molecular device

Interfacing magnetism with superconducting condensates are promising candidates holding Majorana bound states with which fault-tolerant quantum computation could be implemented. Within this field, understanding the detailed dynamics is important both for fundamental reasons and for the development of innovative quantum technologies. Herein, motivated by a molecular magnet Tb2Pc3 interacting with a superconducting Pb(111) substrate, which results in spin-orbital Yu-Shiba-Rusinov (YSR) states, as is affirmed by a theoretical simulation with the aid of the numerical renormalization group technique (see Xia et al 2022 Nat. Commun. 13 6388), we study the YSR states and quantum phase transitions (QPTs) in a bipartite molecular device adsorbed on an wave superconducting substrate. We highlight the effect of asymmetric Coulomb repulsion by computing the spectrum function and spin correlation function in various parameter regimes. We demonstrate that if one impurity is non-interacting, there are no YSR states in both impurities with any repulsion value in the other impurity. Whereas if the repulsion in one impurity is strong, the YSR states are observed in both impurities, and a QPT arises as the repulsion in the other impurity sweeps, assisted by the competition between the superconducting singlet (Cooper pair) and the Kondo singlet. The evolution of YSR states distinguishes from the single impurity case and can be well interpreted by the energy scales of the isotropic superconducting gap parameter, as well as the two Kondo temperatures. Our findings provide theoretical insights into the phase diagram of two magnetic impurities on a superconducting host and shine light on the effects induced by asymmetric Coulomb repulsion on many-body interactions.

Evidence of Kardar-Parisi-Zhang scaling on a digital quantum simulator

Understanding how hydrodynamic behaviour emerges from the unitary evolution of the many-particle Schrödinger equation is a central goal of non-equilibrium statistical mechanics. In this work we implement a digital simulation of the discrete time quantum dynamics of a spin-frac{1}{2} XXZ spin chain on a noisy near-term quantum device, and we extract the high temperature transport exponent at the isotropic point. We simulate the temporal decay of the relevant spin correlation function at high temperature using a pseudo-random state generated by a random circuit that is specifically tailored to the ibmq-montreal 27 qubit device. The resulting output is a spin excitation on a homogeneous background on a 21 qubit chain on the device. From the subsequent discrete time dynamics on the device we are able to extract an anomalous super-diffusive exponent consistent with the conjectured Kardar-Parisi-Zhang (KPZ) scaling at the isotropic point. Furthermore we simulate the restoration of spin diffusion with the application of an integrability breaking potential.

Interplay of dissipation and memory in the quantum Langevin dynamics of a spin in a magnetic field

We derive a quantum Langevin equation for a quantum spin in the presence of a magnetic field and coupled to a bath. We have considered two models for the bath: an Ohmic bath model and a non-Markovian Drude bath model with finite memory. The effect of finite memory on the dynamics of the quantum spin is studied in detail. In particular, we have studied the effect of memory on the expectation values of the magnetic moment and the spin correlation functions. The time evolution study of the spin correlation functions shows that the correlation functions exhibit a damped oscillatory behavior with the randomization time being determined by the damping rate. We have analytically shown that the damping gets affected by the finite memory time when the bath retains memory of the interaction with the system. We also analyze the spin response function of the system for the Ohmic bath. Our results qualitatively agree with the findings in spin noise spectroscopy experiments and our predictions can be tested in future ultracold atom experiments.

Distinction between transport and Rényi entropy growth in kinetically constrained models

Conservation laws and the associated hydrodynamic modes have important consequences on the growth of higher R\'enyi entropies in isolated quantum systems. It has been shown in various random unitary circuits and Hamiltonian systems that the dynamics of the R\'enyi entropies in the presence of a U(1) symmetry obey ${S}^{(n\ensuremath{\ge}2)}(t)\ensuremath{\propto}{t}^{1/z}$, where $z$ is identified as the dynamical exponent characterizing transport of the conserved charges. Here, however, we demonstrate that this simple identification may not hold in certain quantum systems with kinetic constraints. In particular, we study two types of U(1)-symmetric quantum automaton circuits with XNOR (exclusive NOR) and Fredkin constraints, respectively. We find numerically that while spin transport in both models is subdiffusive, the second R\'enyi entropy grows diffusively in the XNOR model, and superdiffusively in the Fredkin model. For systems with the XNOR constraint, this distinction arises since the spin correlation function can be attributed to an emergent tracer dynamics of tagged particles, whereas the R\'enyi entropies are constrained by collective transport of the particles. Our results suggest that care must be taken when relating transport and entanglement entropy dynamics in generic quantum systems with conservation laws.

Spatiotemporal dynamics of classical and quantum density profiles in low-dimensional spin systems

We provide a detailed comparison between the dynamics of high-temperature spatiotemporal correlation functions in quantum and classical spin models. In the quantum case, our large-scale numerics are based on the concept of quantum typicality, which exploits the fact that random pure quantum states can faithfully approximate ensemble averages, allowing the simulation of spin-$1/2$ systems with up to $40$ lattice sites. Due to the exponentially growing Hilbert space, we find that for such system sizes even a single random state is sufficient to yield results with extremely low noise that is negligible for most practical purposes. In contrast, a classical analog of typicality is missing. In particular, we demonstrate that, in order to obtain data with a similar level of noise in the classical case, extensive averaging over classical trajectories is required, no matter how large the system size. Focusing on (quasi-)one-dimensional spin chains and ladders, we find a remarkably good agreement between quantum and classical dynamics. This applies not only to cases where both the quantum and classical model are nonintegrable, but also to cases where the quantum spin-$1/2$ model is integrable and the corresponding classical $s\to\infty$ model is not. Our analysis is based on the comparison of space-time profiles of the spin and energy correlation functions, where the agreement is found to hold not only in the bulk but also in the tails of the resulting density distribution. The mean-squared displacement of the density profiles reflects the nature of emerging hydrodynamics and is found to exhibit similar scaling for quantum and classical models.

Spin Excitation Spectra of Anisotropic Spin-1/2 Triangular Lattice Heisenberg Antiferromagnets.

Investigation of dynamical excitations is difficult but crucial to the understanding of many exotic quantum phenomena discovered in quantum materials. This is particularly true for highly frustrated quantum antiferromagnets whose dynamical properties deviate strongly from theoretical predictions made based on the spin-wave or other approximations. Here, we present a large-scale numerical calculation on the dynamical correlation functions of spin-1/2 triangular Heisenberg model using a state-of-the-art tensor network renormalization group method. The calculated results allow us to gain for the first time a comprehensive picture on the nature of spin excitation spectra in this highly frustrated quantum system. It provides a quantitative account for all the key features of the dynamical spectra disclosed by inelastic neutron scattering measurements for Ba_{3}CoSb_{2}O_{9}, revealing the importance of the interplay between low- and high-energy excitations and its renormalization effect to the low-energy magnon bands and high-energy continuums. We identify the longitudinal Higgs modes in the intermediate-energy scale and predict the energy and momentum dependence of spectral functions along the three principal axes that can be verified by polarized neutron scattering experiments. Furthermore, we find that the spin excitation spectra weakly depend on the anisotropic ratio of the antiferromagnetic interaction.

Room-temperature antiferromagnetic order in monolayer Fe2C

We have studied the room-temperature long-range antiferromagnetic order in the ${\mathrm{Fe}}_{2}\mathrm{C}$ monolayer using the combination of first-principles calculations and the Green's function analysis within the anisotropic Heisenberg model. The ${\mathrm{Fe}}_{2}\mathrm{C}$ monolayer is a semimetal with the out-of-plane antiferromagnetic order between two $\mathrm{Fe}$ planes. The Green's function approach within the random phase approximation is formulated to calculate the temperature-dependent antiferromagnetic magnon energy and the spin correlation function in a two-dimensional antiferromagnetic monolayer. The correlation function is used to evaluate the sublattice magnetization and study the magnetic phase transition in monolayer ${\mathrm{Fe}}_{2}\mathrm{C}$. Moreover, the spin Hamiltonian and the Green's function formalism are developed to investigate the antiferromagnetic ${\mathrm{Fe}}_{2}\mathrm{C}$ monolayer in the presence of an external magnetic field applied along the easy axis. The H-T phase diagram shows that the antiferromagnetic to spin-flop and the spin-flop to paramagnetic phase transitions occur in low temperatures. Finally, we estimate the N\'eel temperature and the critical values of the magnetic field strength for these two phase transitions.

On general-n coefficients in series expansions for row spin–spin correlation functions in the two-dimensional Ising model

We consider spin–spin correlation functions for spins along a row, R n = ⟨σ 0,0 σ n,0⟩, in the two-dimensional Ising model. We discuss a method for calculating general-n expressions for coefficients in high-temperature and low-temperature series expansions of R n and apply it to obtain such expressions for several higher-order coefficients. In addition to their intrinsic interest, these results could be useful in the continuing quest for a nonlinear ordinary differential equation whose solution would determine R n , analogous to the known nonlinear ordinary differential equation whose solution determines the diagonal correlation function ⟨σ 0,0 σ n,n ⟩ in this model.

Analytic average magnetization expression for the body-centered cubic Ising lattice

Recently we have performed large-scale average magnetization studies of the Ising model for the square, honeycomb, triangular, and simple cubic lattice. We want to complement those studies with the structurally more complicated body-centered cubic lattice Ising model in this paper. We have relevantly calculated the order parameter or average magnetization expression for the body-centered cubic Ising lattice based on the recently developed interrelation, given by $\langle\sigma_{0, i}\rangle= \langle\tanh[K(\sigma_{1,i}+\sigma_{2,i}+\dots +\sigma_{z,i})+H]\rangle$. We will use the same conjectures for the odd spins correlation functions as in our previous works. The mathematical forms of the odd spin correlation functions will be eventually elucidated by some physical discussions. We have seen that the obtained average magnetization expression with the numerical data for the body-centered cubic lattice Ising model is in good agreement.

Boundary logarithmic corrections to the dynamical correlation functions of one-dimensional spin- 12 chains

The asymptotic dynamical correlation functions in one-dimensional spin chains are described by power-laws. The corresponding exponents characterize different bulk and boundary critical behavior. We present novel results for the logarithmic contribution to the boundary correlations of an isotropic Heisenberg chain. The exponent of the logarithm, $\lambda=1$, is derived using a renormalization group technique. We confirm our analytical results by comparing with numerical quantum Monte Carlo data.

Correlation functions for open XXZ spin 1/2 quantum chains with unparallel boundary magnetic fields

In this paper we continue our derivation of the correlation functions of open quantum spin 1/2 chains with unparallel magnetic fields on the edges; this time for the more involved case of the XXZ spin 1/2 chains. We develop our study in the framework of the quantum separation of variables, which gives us both the complete spectrum characterization and simple scalar product formulae for separate states, including transfer matrix eigenstates. Here, we leave the boundary magnetic field in the first site of the chain completely arbitrary, and we fix the boundary field in the last site N of the chain to be a specific value along the z-direction. This is a natural first choice for the unparallel boundary magnetic fields. We prove that under these special boundary conditions, on the one side, we have a simple enough complete spectrum description in terms of homogeneous Baxter like TQ-equation. On the other side, we prove a simple enough description of the action of a basis of local operators on transfer matrix eigenstates as linear combinations of separate states. Thanks to these results, we achieve our main goal to derive correlation functions for a set of local operators both for the finite and half-infinite chains, with multiple integral formulae in this last case.