Spin Chain Shocks: Exploring Quantum Dynamics at the Edge of Integrability

Author: Denis Avetisyan


New research delves into the complex behavior of spin-1 chains as they are driven from ordered to chaotic states, revealing surprising insights into their non-equilibrium dynamics.

The evolution of quantum states within a phantom helix - specifically those with a torsional angle of <span class="katex-eq" data-katex-display="false">\pi/2</span> and winding parameters of <span class="katex-eq" data-katex-display="false">\pi/4</span>, <span class="katex-eq" data-katex-display="false">\pi/2</span>, and <span class="katex-eq" data-katex-display="false">3\pi/4</span> - demonstrates a complex interplay between entanglement, fidelity, and correlations, governed by the ratio <span class="katex-eq" data-katex-display="false">J_{z}/J = \cos(Q_{p})</span> and observed across lattices of size 30, suggesting that even seemingly constrained systems harbor rich, dynamic behaviors shaped by subtle shifts in their fundamental parameters.
The evolution of quantum states within a phantom helix – specifically those with a torsional angle of \pi/2 and winding parameters of \pi/4, \pi/2, and 3\pi/4 – demonstrates a complex interplay between entanglement, fidelity, and correlations, governed by the ratio J_{z}/J = \cos(Q_{p}) and observed across lattices of size 30, suggesting that even seemingly constrained systems harbor rich, dynamic behaviors shaped by subtle shifts in their fundamental parameters.

This study investigates quantum quenches in a tunable SU(3)-symmetric spin-1 chain, focusing on entanglement entropy and conserved quantities to characterize the resulting non-equilibrium behavior.

Understanding the dynamics of many-body quantum systems far from equilibrium remains a central challenge in modern physics. This work, ‘Quantum quenches in a spin-1 chain with tunable symmetry’, investigates non-equilibrium phenomena in an anisotropic spin-1 Heisenberg chain, leveraging the time-evolving block decimation (TEBD) method to explore a parameter space tunable between integrable and non-integrable regimes. By modulating the quadrupolar interaction strength, we demonstrate the emergence of a novel conserved quantity at the SU(3) symmetric point and characterize its influence on local observables, entanglement entropy, and spin correlations. How might these findings inform the design of quantum simulators capable of realizing and controlling complex, out-of-equilibrium states in spin-1 lattice models?


Beyond Simplification: The Echo of True Spin Interaction

Conventional spin models, while foundational in understanding magnetic phenomena, frequently operate under simplifying assumptions that neglect the nuanced reality of interacting quantum spins. These models often treat spins as purely dipolar, overlooking the significant contribution of quadrupolar interactions arising from the non-spherical charge distribution of certain atomic nuclei. This simplification proves inadequate when studying materials exhibiting strong spin-orbit coupling or containing nuclei with significant quadrupolar moments. Consequently, these traditional approaches struggle to accurately describe phenomena such as anisotropic magnetic behavior, complex magnetic ordering, and the influence of electric field gradients on spin dynamics. The failure to incorporate these quadrupolar effects limits the predictive power of standard models, necessitating more sophisticated frameworks capable of capturing the full complexity of spin interactions within a material.

The Anisotropic Heisenberg Hamiltonian represents a significant advancement in the study of magnetic materials, offering a more complete description than traditional isotropic models. While simpler models assume spin interactions are uniform in all directions, this Hamiltonian accounts for the inherent directionality often present due to crystal structures and electronic configurations. This capability allows researchers to explore phenomena-like magnetic anisotropy and complex ordering patterns-that remain hidden when using conventional approaches. By incorporating parameters that define the strength of interactions along specific axes, the Hamiltonian accurately simulates the behavior of materials where spin alignment is not equivalent in all directions, ultimately providing a more realistic and nuanced understanding of their magnetic properties and potential applications.

The Anisotropic Heisenberg Hamiltonian distinguishes itself by directly incorporating quadrupolar operators, a critical feature for accurately modeling systems where interactions aren’t uniform in all directions. These operators, representing the electric quadrupole moment of a nucleus or electron, account for the distribution of charge that deviates from spherical symmetry. Consequently, the Hamiltonian can describe anisotropic exchange interactions – forces that depend on the relative orientation of spins – which are prevalent in materials exhibiting complex magnetic behavior. Unlike traditional models that assume isotropic interactions, this approach allows for a nuanced understanding of how the shape of the electronic charge distribution influences magnetic ordering and excitations, offering a more complete picture of phenomena observed in materials like certain rare-earth compounds and strongly correlated electron systems. The inclusion of these operators, mathematically represented by tensor quantities, fundamentally alters the energy landscape of the system, leading to unique magnetic phases and dynamic responses.

The model, described by <span class="katex-eq" data-katex-display="false">Eq.II.1</span>, transitions between non-integrable SU(2) symmetry and integrable SU(3) symmetry based on the strength of quadrupolar interactions, parameterized by <span class="katex-eq" data-katex-display="false">J_q/J</span>.
The model, described by Eq.II.1, transitions between non-integrable SU(2) symmetry and integrable SU(3) symmetry based on the strength of quadrupolar interactions, parameterized by J_q/J.

The Tuning Parameter: A Gateway to Symmetry

The parameter J_q/J governs the relative strength of isotropic and anisotropic exchange interactions within the system. J represents the strength of the isotropic Heisenberg interaction, promoting spin alignment in all directions, while J_q quantifies the quadrupolar interaction, which favors specific orientations of spins due to the electronic structure of the ions. A higher ratio of J_q/J increases the influence of the anisotropic quadrupolar interaction, while a lower ratio emphasizes the isotropic Heisenberg exchange. This control over the balance between these competing interactions allows for precise tuning of the system’s magnetic properties and dynamical behavior.

When the quadrupolar interaction is strong, indicated by a high value of J_q/J, the system’s behavior converges towards that of the integrable SU3 Heisenberg Model. This convergence results in the presence of conserved quantities, most notably M^2, which represents the squared eigenvalue of the quadratic Casimir operator for the SU3 symmetry. The existence of such conserved quantities significantly constrains the system’s dynamics, allowing for exact solutions and predictable long-term behavior. This integrability arises from the enhanced symmetry present at higher values of J_q/J, effectively limiting the available phase space for the system’s evolution.

Decreasing the J_q/J parameter effectively weakens quadrupolar interactions, pushing the system’s behavior closer to that of the non-integrable SU2 Heisenberg model. This transition is significant because the SU2 model lacks the conserved quantities present in the integrable SU3 model, preventing solutions via traditional analytical methods. Consequently, reducing J_q/J introduces chaotic dynamics and allows for the observation of tunable system behavior, effectively creating a pathway to explore the transition from an integrable to a non-integrable regime. This tunability enables controlled investigations into the emergence of complex dynamics and the breakdown of established analytical techniques.

The time evolution of spin and quadrupole correlations between adjacent sites reveals consistent behavior across different initial states, system sizes, and values of <span class="katex-eq" data-katex-display="false">J_q/J</span>, as demonstrated by the patterns in the out-of-plane spin, in-plane spin, out-of-plane quadrupole, and in-plane quadrupole correlations.
The time evolution of spin and quadrupole correlations between adjacent sites reveals consistent behavior across different initial states, system sizes, and values of J_q/J, as demonstrated by the patterns in the out-of-plane spin, in-plane spin, out-of-plane quadrupole, and in-plane quadrupole correlations.

Probing the System: Initial States and Validation

Initialization of the system in Nematic (NM) and Domain Wall (DW) states allows for the direct observation of quadrupolar interaction effects on spin ordering. Nematic states, characterized by long-range orientational order, provide a platform to study how quadrupolar interactions promote alignment of spins along a preferred direction. Domain Wall states, representing boundaries between regions of differing spin orientation, facilitate the investigation of how these interactions mediate correlations across spatial domains. By analyzing the resulting spin configurations and correlation functions from these initial conditions, we can quantify the strength and range of the quadrupolar interactions and their influence on the system’s ground state properties and low-energy excitations. These states serve as controlled starting points to probe the mechanisms driving spin alignment and the emergence of spatial correlations due to the quadrupolar coupling.

Phantom Helix States are utilized to analyze the system’s behavior when subjected to a global momentum input, enabling the study of its non-equilibrium dynamics. These states introduce a defined, macroscopic momentum which drives the system away from equilibrium, allowing observation of the subsequent relaxation processes and the emergence of spatial correlations. Analysis of the system’s response to this induced momentum provides information regarding the transport properties and the mechanisms governing energy dissipation within the system, complementing the static analysis provided by the Nematic and Domain Wall initial states.

Exact Diagonalization (ED) was utilized as a benchmark to validate the accuracy of Time-Evolving Block Decimation (TEBD) simulations performed on the system. This validation process is critical for ensuring the reliability of results obtained from the TEBD method, which approximates the time evolution of quantum states. TEBD simulations were optimized by systematically increasing the bond dimension, χ, until convergence was achieved. Specifically, a bond dimension of \chi = 600 was found to be sufficient, as results remained consistent for values exceeding 400, indicating a negligible impact of truncation error at this level.

TEBD simulations of spin chains reveal that the absolute difference in local magnetization, entanglement entropy, and fidelity remain stable across varying time steps (<span class="katex-eq" data-katex-display="false">dt=5\times 10^{-2}, 10^{-2}, 5\times 10^{-3}</span>) for initial states including DW I, NPI, and phantom helix with <span class="katex-eq" data-katex-display="false">Q_p = \pi/4</span> and <span class="katex-eq" data-katex-display="false">\pi/2</span>, indicating the robustness of the results given <span class="katex-eq" data-katex-display="false">J_{xy} = J_q</span> and differing <span class="katex-eq" data-katex-display="false">J_z</span> values.
TEBD simulations of spin chains reveal that the absolute difference in local magnetization, entanglement entropy, and fidelity remain stable across varying time steps (dt=5\times 10^{-2}, 10^{-2}, 5\times 10^{-3}) for initial states including DW I, NPI, and phantom helix with Q_p = \pi/4 and \pi/2, indicating the robustness of the results given J_{xy} = J_q and differing J_z values.

Fragmentation and the Limits of Thermalization

The Anisotropic Heisenberg Hamiltonian, a cornerstone in modeling interacting quantum spins, surprisingly displays a phenomenon known as Hilbert Space Fragmentation. This isn’t a traditional breakdown of the system, but rather a partitioning of its entire state space – the Hilbert space – into isolated, disconnected sectors. Each sector behaves as a self-contained quantum system, unable to directly interact with states in other sectors, even with strong interactions between the constituent spins. This fragmentation fundamentally alters the system’s dynamics, preventing it from exploring its full range of possible states and leading to a suppression of thermalization – the expected equilibration to a uniform temperature. Consequently, the system retains ‘memory’ of its initial conditions, defying the typical ergodic behavior where all accessible states are eventually visited with equal probability. This discovery challenges conventional understandings of how quantum systems evolve and opens new avenues for exploring many-body localization and the preservation of quantum information.

The emergence of Hilbert space fragmentation fundamentally alters the expected behavior of quantum systems, giving rise to localization phenomena that actively impede thermalization. Typically, closed quantum systems are expected to eventually reach a state of thermal equilibrium, where energy is evenly distributed across all available degrees of freedom – a concept rooted in ergodicity. However, fragmentation creates disconnected sectors within the system’s Hilbert space, effectively trapping excitations and preventing the free flow of energy. This localization means that energy remains confined to specific regions, inhibiting the system from exploring all possible states and thus suppressing the approach to thermal equilibrium. The resulting non-ergodic behavior challenges long-held assumptions about how quantum systems evolve and suggests the potential for novel states of matter where energy transport is fundamentally limited.

The degree of Hilbert space fragmentation and its subsequent impact on quantum information were carefully quantified through a suite of observables. Measurements of Local Magnetization revealed spatially localized excitations, while calculations of Entanglement Entropy consistently yielded values around S \approx 6.4, indicating a significant departure from thermalized states. Fidelity measurements further confirmed the persistence of these non-ergodic behaviors. These calculations were enabled by Time-Evolving Block Decimation (TEBD) simulations optimized with a time step of 5 \times 10^{-3}; this refinement reduced accumulated error by up to four orders of magnitude when compared to simulations employing a larger time step of 5 \times 10^{-2}, ensuring the reliability and precision of the reported fragmentation effects.

Relative frequencies of local bond alignments within Hilbert subspaces reveal distinctions between symmetry sectors with <span class="katex-eq" data-katex-display="false">J_q/J = 1</span> (red) and <span class="katex-eq" data-katex-display="false">J_q/J \neq 1</span> (blue) for initial states including DW I, DW II, AFM I, AFM II, NM, and NPI, categorized by alignment types like FM-Z, AFM-Z, FQ-Z, and AFQ-dZ.
Relative frequencies of local bond alignments within Hilbert subspaces reveal distinctions between symmetry sectors with J_q/J = 1 (red) and J_q/J \neq 1 (blue) for initial states including DW I, DW II, AFM I, AFM II, NM, and NPI, categorized by alignment types like FM-Z, AFM-Z, FQ-Z, and AFQ-dZ.

The study of this spin-1 chain, oscillating between order and chaos, echoes a fundamental truth about all complex systems. It isn’t about building a static structure, but rather observing the unfolding of inherent dynamics. The researchers demonstrate how subtle parameter tuning dramatically alters the system’s response – a delicate balance between integrable and non-integrable regimes. As Isaac Newton observed, “An object at rest stays at rest and an object in motion stays in motion with the same speed and in the same direction unless acted upon by a force.” This principle extends metaphorically to these quantum systems; the ‘force’ being the manipulation of symmetry and initial conditions, revealing conserved quantities and intricate entanglement behaviors. The system isn’t silent, but constantly revealing itself through non-equilibrium dynamics, and its ultimate state is less a destination and more a continuous trajectory.

What Lies Ahead?

The exploration of non-equilibrium dynamics, even within the seemingly controlled environment of a spin-1 chain, reveals a truth often obscured by the pursuit of precision: systems do not solve problems; they redistribute them. The identification of novel, emergent conserved quantities is not a destination, but a map pointing to previously unseen regions of instability. Long-lived stability, particularly in regimes approaching integrability, should not be interpreted as robustness, but as the prolonged accumulation of potential for unexpected evolution.

Future work will undoubtedly focus on expanding the scope of these investigations. However, a more fruitful path may lie in accepting the inherent limitations of current analytical tools. The temptation to engineer ‘solutions’ through precise parameter tuning should be resisted. A truly revealing study would embrace disorder, examining how these systems fail to maintain coherence under realistic conditions – how they transform into something other than what was initially intended.

The observed sensitivity to initial states suggests that the precise definition of ‘equilibrium’ is itself a transient illusion. The question is not whether these systems can be driven to a desired state, but how they will become something else entirely. This chain, and others like it, are not tools to be wielded, but gardens to be observed – and the most interesting blooms will always be the ones that defy expectation.


Original article: https://arxiv.org/pdf/2604.18425.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/

See also:

2026-04-21 15:55