Chiral order emergence driven by quenched disorder

Basic Information

• Paper ID: 2503.24139
• Title: Chiral order emergence driven by quenched disorder
• Authors: Coraline Letouzé, Pascal Viot, Laura Messio (Sorbonne Université, CNRS, LPTMC)
• Classification: cond-mat.stat-mech, cond-mat.str-el
• Publication Date: September 24, 2025
• Paper Link: https://arxiv.org/abs/2503.24139v3

Abstract

Quenched disorder typically destroys magnetic order, such as when random fields are applied in the two-dimensional Ising model. Even when order persists in the presence of quenched disorder, it usually represents merely a remnant of the order in the clean model. This paper demonstrates a surprising phenomenon: order emergence driven by quenched disorder. This order bears no resemblance to the order present in the clean model. This "disorder-induced order" mechanism differs from conventional thermodynamic or quantum disorder-induced order. Using parallel tempering Monte Carlo simulations, we investigate the classical $J_1-J_3$ Heisenberg model with site dilution on the kagome lattice. After analyzing the effects of a small number of vacancies on the ground state (which favor non-coplanar configurations), we demonstrate the emergence of a low-temperature chiral phase and the progressive destruction of collinear $q=4$ Potts order.

Research Background and Motivation

Significance of the Problem

1. Challenge to Conventional Wisdom: Traditionally, quenched disorder is believed to suppress or destroy various ordered states. This research challenges this conventional view by discovering that quenched disorder can drive the emergence of entirely new ordered phases.
2. Theoretical Significance: In two-dimensional frustrated magnets with continuous spin degrees of freedom, non-collinear or non-coplanar spin order may break discrete symmetries (such as time-reversal or lattice symmetries), which according to the Mermin-Wagner theorem allows finite-temperature phase transitions.
3. Existing Limitations:
   • Existing "disorder-induced order" mechanisms primarily include thermodynamic disorder-induced order (ObTD) and quantum disorder-induced order (ObQD), which typically select the most collinear spin configurations
   • Research on quenched disorder has primarily focused on its destructive effects, lacking deep understanding of its constructive role
4. Research Motivation: To explore whether quenched disorder can induce the emergence of new ordered phases, particularly chiral order in O(3)-symmetric models, which has important theoretical significance.

Core Contributions

1. Discovery of a New Disorder-Induced Order Mechanism: First observation of discrete chiral order emergence driven by quenched disorder in models with O(3) spin rotational symmetry
2. Establishment of Complete Phase Diagram: Through large-scale Monte Carlo simulations, construction of the temperature-dilution rate phase diagram for the $J_1-J_3$ Heisenberg kagome model under site dilution
3. Elucidation of Theoretical Mechanism: Microscopic explanation of how vacancies destroy K4 order through effective random fields while promoting chiral order formation
4. Verification of Universality Classes: Through finite-size scaling analysis, confirmation that the chiral phase transition belongs to the Ising universality class, while the K4 phase transition belongs to the 4-state Potts universality class

Methodology Details

Task Definition

Investigation of phase transition behavior in the $J_1-J_3$ Heisenberg model on the kagome lattice under site dilution conditions, with particular focus on the emergence mechanism of chiral order.

Model Architecture

Hamiltonian

The Hamiltonian under site dilution conditions is: $H = J_1 \sum_{\langle i,j\rangle} \eta_i\eta_j \mathbf{S}_i \cdot \mathbf{S}_j + J_3 \sum_{\langle i,j\rangle_3} \eta_i\eta_j \mathbf{S}_i \cdot \mathbf{S}_j$

where:

• $\mathbf{S}_i$ is a classical three-component unit vector
• $\eta_i = 1$ or $0$ indicates the presence or absence of magnetic spin at site $i$
• The dilution rate is defined as $p = 1 - \sum_i \eta_i/N$

Order Parameter Definitions

1. K4 Order Parameter $\Sigma$: Based on local variables on triangles $\sigma_{ABC} = \begin{pmatrix} \mathbf{S}_B \cdot \mathbf{S}_C \\ \mathbf{S}_C \cdot \mathbf{S}_A \\ \mathbf{S}_A \cdot \mathbf{S}_B \end{pmatrix}$
2. Chiral Order Parameter $\chi$: $\chi_{ABC} = \mathbf{S}_A \cdot (\mathbf{S}_B \times \mathbf{S}_C)$

Numerical Methods

Monte Carlo Simulations

• Algorithm: Parallel tempering method combined with heat bath local updates
• System Sizes: Linear size $L$ ranging from 24 to 128 ($N = 3L^2$ lattice sites)
• Disorder Averaging: Average over at least 10 disorder realizations for each $L$
• Parameter Selection: $J_1 = -1$, $J_3 = 1$

Finite-Size Scaling Analysis

The following relationships are used to determine critical exponents of phase transitions: $C_V^{\max}(L) = \frac{\langle H^2\rangle - \langle H\rangle^2}{NT^2} \simeq aL^{\alpha/\nu} + b$$S_X^{\max}(L) = \frac{\langle X^2\rangle - \langle X\rangle^2}{NT} \simeq cL^{\gamma/\nu} + d$$T_{Y^{\max}}(L) \simeq eL^{-1/\nu} + T_c$

Technical Innovations

1. Microscopic Mechanism Analysis: Revelation of how vacancies create effective magnetic fields through single-vacancy analysis, leading to non-collinear orientation of adjacent spins
2. Random Field Effects: Clarification of how site disorder produces an effective random field on K4 order, which according to the Imry-Ma argument destroys discrete order in two-dimensional systems
3. Chiral Emergence Mechanism: Explanation of how the three sublattices tend toward mutually perpendicular orientations in the presence of vacancies, thereby determining local chirality

Experimental Setup

Simulation Parameters

• Exchange Parameters: $J_1 = -1$ (ferromagnetic), $J_3 = 1$ (antiferromagnetic)
• Dilution Range: $0 \leq p \leq 0.2$
• Temperature Range: Adjusted according to phase diagram, covering all phase transition regions
• System Sizes: $L = 24, 36, 48, 64, 80, 96, 128$

Evaluation Metrics

1. Specific Heat $C_V$: Identifies phase transition locations
2. Order Parameter Susceptibilities $S_\Sigma$, $S_\chi$: Determines phase transition nature
3. Order Parameters $\Sigma$, $\chi$: Quantifies order strength

Data Processing

• Sufficient thermalization and statistics for each parameter point
• Averaging over multiple disorder realizations
• Extraction of critical exponents using finite-size scaling analysis

Experimental Results

Main Results

Phase Diagram Characteristics

A complete temperature-dilution rate phase diagram has been established, identifying three main regions:

1. Small Dilution Region ($p < 0.05$):
   • High-temperature K4 ordered phase transition with critical temperature decreasing linearly with $p$
   • Low-temperature chiral phase transition emergence with $T_c \approx 4p$
2. Intermediate Dilution Region ($0.05 \leq p \leq 0.1$):
   • Disappearance of K4 order
   • Persistence of chiral order with critical temperature increasing with $p$
3. Large Dilution Region ($0.1 < p < 0.2$):
   • Existence of chiral order only
   • Disappearance of chiral order at $p = 0.2$

Critical Exponent Verification

• K4 Phase Transition: Belongs to 4-state Potts universality class ($\alpha/\nu = 1$, $\gamma/\nu = 7/4$, $1/\nu = 3/2$)
• Chiral Phase Transition: Belongs to Ising universality class ($\alpha/\nu = 0$, $\gamma/\nu = 7/4$, $1/\nu = 1$)

Key Findings

1. Vacancy Effect Analysis: A single vacancy on the B sublattice produces an effective magnetic field of strength $J_1$, causing A and C sublattice spins to tend toward perpendicular orientation relative to the missing spin direction
2. Random Field Mechanism: Site disorder produces an effective random field on K4 order, with energy differences depending on vacancy location: $8J_3 + J_1$ for same-color sites and $8J_3 - J_1$ for different-color sites
3. Percolation Threshold: When $p = 0.2$, the probability of complete triangles $(1-p)^3 = 0.512$ approaches the site percolation threshold of the honeycomb lattice 0.697

Related Work

Disorder-Induced Order Research

1. Thermodynamic Disorder-Induced Order: Pioneering work by Villain et al., where entropy contributions select ground states with high density of low-energy excitations
2. Quantum Disorder-Induced Order: Quantum fluctuations select the most collinear spin configurations
3. Quenched Disorder Effects: Traditionally believed to destroy order; this work first reveals its constructive role

Kagome Lattice Research

1. Pure Model: Known to exhibit rich zero-temperature phase diagrams and ObTD phenomena
2. Experimental Realization: The Ba-Vesignieite compound BaCu₃V₂O₈(OH)₂ serves as an experimental counterpart to this model

Conclusions and Discussion

Main Conclusions

1. First observation of chiral order emergence driven by quenched disorder in O(3)-symmetric Heisenberg models
2. Establishment of a complete phase diagram revealing the transition process from K4 order to chiral order
3. Proposal of a new "quenched disorder-induced order" mechanism, distinct from conventional thermodynamic or quantum disorder-induced order

Limitations

1. Model Specificity: Research limited to the specific $J_1-J_3$ kagome model
2. Size Constraints: System sizes limited to $L \leq 80$ in the large dilution region due to thermalization difficulties
3. Temperature Range: Numerical precision in the low-temperature region is limited

Future Directions

1. Universality Verification: Verification of this mechanism in other frustrated magnetic systems
2. Experimental Verification: Search for experimental systems to verify theoretical predictions
3. Quantum Effects: Investigation of analogous phenomena in quantum spin systems

In-Depth Evaluation

Strengths

1. Theoretical Innovation: Discovery of a novel disorder-induced order mechanism challenging conventional wisdom
2. Numerical Rigor: Large-scale Monte Carlo simulations with sufficient finite-size scaling analysis
3. Clear Mechanism: Complete theoretical explanation from microscopic to macroscopic scales
4. Experimental Relevance: Connection to real material systems

Weaknesses

1. System Specificity: Universality of results requires further verification
2. Numerical Limitations: Numerical precision in certain parameter regions needs improvement
3. Experimental Verification: Lack of direct experimental verification

Impact

1. Theoretical Contribution: Provides new perspective for frustrated magnet theory
2. Practical Value: Provides guiding principles for material design
3. Methodological: Serves as an example for studying quenched disorder systems

Applicable Scenarios

1. Frustrated Magnetic Systems: Particularly systems with multi-sublattice structures
2. Material Design: Realization of specific magnetic order through defect control
3. Theoretical Research: Reference for other disorder-induced order phenomena

References

This paper cites 48 relevant references, covering important works in multiple fields including disorder-induced order, frustrated magnets, and phase transition theory, providing a solid theoretical foundation for the research.

Computational Environment Impact: This project utilized approximately 40,000 core-days of computational resources, equivalent to 2,800 kg CO₂ emissions, reflecting consideration of environmental costs in large-scale numerical research.