We investigate an extension of the quantum Ising (QI) model in one spatial dimension including long-range interactions in its statics and dynamics with possible applications from heteronuclear polar molecules in optical lattices to trapped ions described by two-state spin systems. We introduce the statics of the system via both numerical techniques with finite size and infinite size matrix product states (MPSs) and a theoretical approaches using a truncated Jordan–Wigner transformation for the ferromagnetic and antiferromagnetic case and show that finite size effects have a crucial role shifting the quantum critical point of the external field by fifteen percent between thirty-two and around five-hundred spins. We numerically study the Kibble–Zurek hypothesis in the long-range QI model with MPSs. A linear quench of the external field through the quantum critical point yields a power-law scaling of the defect density as a function of the total quench time. For example, the increase of the defect density is slower for longer-range models and the critical exponent changes by twenty-five percent. Our study emphasizes the importance of such long-range interactions in statics and dynamics that could point to similar phenomena in a different setup of dynamical systems or for other models.
We investigate an extension of the quantum Ising (QI) model in one spatial dimension including long-range interactions in its statics and dynamics with possible applications from heteronuclear polar molecules in optical lattices to trapped ions described by two-state spin systems. We introduce the statics of the system via both numerical techniques with finite size and infinite size matrix product states (MPSs) and a theoretical approaches using a truncated Jordan–Wigner transformation for the ferromagnetic and antiferromagnetic case and show that finite size effects have a crucial role shifting the quantum critical point of the external field by fifteen percent between thirty-two and around five-hundred spins. We numerically study the Kibble–Zurek hypothesis in the long-range QI model with MPSs. A linear quench of the external field through the quantum critical point yields a power-law scaling of the defect density as a function of the total quench time. For example, the increase of the defect density is slower for longer-range models and the critical exponent changes by twenty-five percent. Our study emphasizes the importance of such long-range interactions in statics and dynamics that could point to similar phenomena in a different setup of dynamical systems or for other models.
We analyze the quantum phase transitions taking place in a one-dimensional transverse-field Ising model with long-range couplings that decay algebraically with distance. We are interested in the Kibble-Zurek universal scaling laws emerging in nonequilibrium dynamics and in the potential for the unambiguous observation of such behavior in a realistic experimental setup based on trapped ions. To this end, we determine the phase diagram of the model and the critical exponents characterizing its quantum phase transitions by means of density-matrix renormalization group calculations and finite-size scaling theory, which allows us to obtain good estimates for different range of ferromagnetic and antiferromagnetic interactions. Beyond critical equilibrium properties, we tackle a nonequilibrium scenario in which quantum Kibble-Zurek scaling laws may be retrieved. Here, it is found that the predicted nonequilibrium universal behavior, i.e., the scaling laws as a function of the quench rate and critical exponents, can be observed in systems comprising an experimentally feasible number of spins. Finally, a scheme is introduced to simulate the algebraically decaying couplings accurately by means of a digital quantum simulation with trapped ions. Our results suggest that quantum Kibble-Zurek physics can be explored and observed in state-of-the-art experiments with trapped ions realizing long-range Ising models.
The trapped-ion quantum simulator has demonstrated qualitative properties of different physical models for up to tens of ions.In particular, a linear ion chain naturally hosts long-range Ising interactions under the laser driving, which has been used for various phenomena, such as quantum phase transition, localization, thermalization, and information propagation.For near-term practical usage, a central task is to find more quantitative applications of the noisy quantum simulators that are robust to small errors in the parameters.Here we report the quantum simulation of a long-range transverse-field Ising model using up to 61 ions and probe the critical behavior of its quantum phase transition through the Kibble-Zurek mechanism.By calibrating and verifying the coupling coefficients, we realize the same model for increasing ion numbers, so as to extract a critical exponent free of the finite-size effect.For ferromagnetic interaction, our experimental result agrees with the previous numerical prediction.As for the antiferromagnetic case, signals are too weak to fit a critical exponent due to the frustration in the interaction, but still consistent with the theory.
When a system is swept through a quantum critical point, the quantum Kibble-Zurek mechanism makes universal predictions for quantities such as the number and energy of excitations produced. This mechanism is now being used to obtain critical exponents on emerging quantum computers and emulators, which in some cases can be compared to matrix product state (MPS) numerical studies. However, the mechanism is modified when the divergence of entanglement entropy required for a faithful description of many quantum critical points is not fully captured by the experiment or classical calculation. In this Letter, we study how low-energy dynamics of quantum systems near criticality are modified by finite entanglement, using conformally invariant critical points described approximately by a MPS as an example. We derive that the effect of finite entanglement on a Kibble-Zurek process is captured by a dimensionless scaling function of the ratio of two length scales, one determined dynamically and one by the entanglement restriction. Numerically we confirm first that dynamics at finite bond dimension χ is independent of the algorithm chosen, then obtain scaling collapses for sweeps in the transverse field Ising model and the three-state Potts model. Our result establishes the precise role played by entanglement in time-dependent critical phenomena and has direct implications for quantum state preparation and classical simulation of quantum states.
Slow quenches of the magnetic field across the paramagnetic-ferromagnetic phase transition of spin systems produce heat. In systems with short-range interactions the heat exhibits universal power-law scaling as a function of the quench rate, known as Kibble-Zurek scaling. In this work we analyze slow quenches of the magnetic field in the Lipkin-Meshkov-Glick (LMG) model, which describes fully connected quantum spins. We analytically determine the quantum contribution to the residual heat as a function of the quench rate δ by means of a Holstein-Primakoff expansion about the mean-field value. Unlike in the case of short-range interactions, scaling laws in the LMG model are only found for a ramp starting or ending at the critical point. If instead the ramp is symmetric, as in the typical Kibble-Zurek scenario, then the number of excitations exhibits a crossover behavior as a function of δ and tends to a constant in the thermodynamic limit. Previous, and seemingly contradictory, theoretical studies are identified as specific limits of this dynamics. Our results can be tested on several experimental platforms, including quantum gases and trapped ions.
State-of-the-art tensor networks are employed to simulate the Hamiltonian ramp in the analog-digital quantum simulation of the quantum phase transition to the quasi-long-range ordered phase of the two-dimensional square-lattice XX model [T. I. Andersen et al. Nature (London) 638, 79 (2025)]. We focus on the quantum Kibble-Zurek (KZ) mechanism near the quantum critical point. Using the infinite projected entangled pair state, we simulate an infinite lattice and demonstrate the KZ scaling hypothesis for the XX correlations across a wide range of ramp times. We use the time-dependent variational principle algorithm to simulate a finite $8\ifmmode\times\else\texttimes\fi{}8$ lattice, similar to the one in the quantum simulation, and find that adiabatic finite-size effects dominate for longer ramp times, where the correlation length's growth with increasing ramp time saturates and the excitation energy's dependence on the ramp time crosses over to a power-law decay characteristic of adiabatic transitions. This finding contradicts the quantum simulation data where the correlation length seems to obey KZ-like power laws, although with modified exponents.
When a quantum system exhibiting a second order phase transition is quenched across the critical point in large but finite time, the dynamics are not adiabatic in the critical region and the Kibble-Zurek (KZ) mechanism provides a framework to determine local observables such as the mean defect density. However, to find higher-point functions, one has to go beyond the KZ paradigm asshown in recent works on one-dimensional transverse field Ising model (TFIM) following a linear quench. It has been found that (i) besides the KZ scale, the quench dynamics depend on another length scale that arises due to the finite phase difference between the low energy modes, and (ii) contrary to the expectations based on the KZ mechanism, in general, the correlation functions do not decay exponentially with distance. Motivated by these results for the linear quench, we are interested in understanding if these properties are universal, and consider the 1D TFIM when the transverse field varies algebraically in the vicinity of the critical field. We focus on the equal-time,longitudinal kink-kink correlation function at the end of the quench from the paramagnetic to the ferromagnetic phase, and find that (i) the correlator depends only on the KZ length for superlinear quenches, otherwise an additional dephasing length is required to describe it, and (ii) the dephased correlator decays as a compressed exponential with an exponent that changes continuously with the quench exponent. Our results are obtained using an adiabatic perturbation theory, analytical arguments and exact numerical integration of the relevant equations.