2025-11-25T06:07:17.212139

Chiral order emergence driven by quenched disorder

Letouzé, Viot, Messio
Quenched disorder can destroy magnetic order, for example when a random field is applied in a 2-dimensional Ising model. Even when an order exists in the presence of quenched disorder, it is usually only the survival of the order of the clean model. We present here a surprising phenomenon where an order emerges, driven by quenched disorder. This order has nothing in common with the order present in the clean model. This type of \textit{order by disorder} differs from the usual thermal or quantum one. The classical $J_1-J_3$ Heisenberg model on the kagome lattice is studied by parallel tempering Monte Carlo simulations, with site dilution. After analyzing the effect of a few vacancies on the ground state, favoring non-coplanar configurations, we show the emergence of a low-temperature chiral phase and the progressive destruction of the collinear $q=4$ Potts order, the only order present in the absence of vacancies.
academic

Emergenza dell'ordine chirale guidata dal disordine congelato

Informazioni Fondamentali

  • ID Articolo: 2503.24139
  • Titolo: Emergenza dell'ordine chirale guidata dal disordine congelato
  • Autori: Coraline Letouzé, Pascal Viot, Laura Messio (Sorbonne Université, CNRS, LPTMC)
  • Classificazione: cond-mat.stat-mech, cond-mat.str-el
  • Data di Pubblicazione: 24 settembre 2025
  • Link Articolo: https://arxiv.org/abs/2503.24139v3

Riassunto

Il disordine congelato generalmente distrugge l'ordine magnetico, ad esempio quando si applica un campo casuale al modello di Ising bidimensionale. Anche quando l'ordine persiste in presenza di disordine congelato, di solito rappresenta solo un residuo dell'ordine del modello puro. Questo articolo dimostra un fenomeno sorprendente: l'emergenza di ordine guidata dal disordine congelato. Questo ordine non ha alcuna relazione con l'ordine presente nel modello puro. Questo meccanismo di "disordine che induce ordine" differisce dal consueto disordine termico o quantistico che induce ordine. Attraverso simulazioni Monte Carlo con tempering parallelo, è stato studiato il modello di Heisenberg classico J1J3J_1-J_3 sul reticolo kagome con diluizione di siti. Dopo aver analizzato l'influenza di poche vacanze sullo stato fondamentale (che favorisce configurazioni non complanari), è stata dimostrata l'emergenza della fase chirale a bassa temperatura e la progressiva distruzione dell'ordine di Potts colineare q=4q=4.

Contesto di Ricerca e Motivazione

Importanza del Problema

  1. Sfida alla conoscenza tradizionale: Tradizionalmente, il disordine congelato è considerato come inibitore o distruttore di vari stati ordinati. Questa ricerca sfida questo punto di vista convenzionale, scoprendo che il disordine congelato può guidare l'emergenza di stati ordinati completamente nuovi.
  2. Significato Teorico: Nei magneti frustrati bidimensionali con gradi di libertà di spin continui, l'ordine di spin non colineare o non complanare può rompere simmetrie discrete (come l'inversione temporale o le simmetrie reticolari), il che, secondo il teorema di Mermin-Wagner, consente transizioni di fase a temperatura finita.
  3. Limitazioni Esistenti:
    • I meccanismi esistenti di "disordine che induce ordine" includono principalmente il disordine termico che induce ordine (ObTD) e il disordine quantistico che induce ordine (ObQD), che generalmente selezionano le configurazioni di spin più colineari
    • La ricerca sul disordine congelato si è concentrata principalmente sui suoi effetti distruttivi, mancando di una comprensione profonda dei suoi effetti costruttivi
  4. Motivazione della Ricerca: Esplorare se il disordine congelato può indurre l'emergenza di nuovi stati ordinati, in particolare l'ordine chirale nei modelli con simmetria O(3), che ha un significato teorico importante.

Contributi Principali

  1. Scoperta di un nuovo meccanismo di disordine che induce ordine: Prima osservazione dell'emergenza di ordine chirale discreto guidato dal disordine congelato in modelli con simmetria di rotazione di spin O(3)
  2. Costruzione di un diagramma di fase completo: Attraverso simulazioni Monte Carlo su larga scala, è stato costruito il diagramma di fase temperatura-tasso di diluizione del modello di Heisenberg kagome J1J3J_1-J_3 sotto diluizione di siti
  3. Chiarimento del meccanismo teorico: Spiegazione dal punto di vista microscopico di come le vacanze distruggono l'ordine K4 attraverso la produzione di campi casuali efficaci, promuovendo contemporaneamente la formazione dell'ordine chirale
  4. Verifica della classe di universalità: Attraverso l'analisi di scaling a dimensione finita, è stato confermato che la transizione di fase chirale appartiene alla classe di universalità di Ising, mentre la transizione di fase K4 appartiene alla classe di universalità di Potts a 4 stati

Descrizione Dettagliata dei Metodi

Definizione del Compito

Studio del comportamento delle transizioni di fase del modello di Heisenberg J1J3J_1-J_3 sul reticolo kagome in condizioni di diluizione di siti, con particolare attenzione al meccanismo di emergenza dell'ordine chirale.

Architettura del Modello

Hamiltoniana

L'Hamiltoniana in condizioni di diluizione di siti è: H=J1i,jηiηjSiSj+J3i,j3ηiηjSiSjH = 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

dove:

  • Si\mathbf{S}_i è un vettore unitario classico a tre componenti
  • ηi=1\eta_i = 1 o 00 indica la presenza o l'assenza di uno spin magnetico nel sito ii
  • Il tasso di diluizione è definito come p=1iηi/Np = 1 - \sum_i \eta_i/N

Definizione dei Parametri d'Ordine

  1. Parametro d'ordine K4 Σ\Sigma: basato su variabili locali sui triangoli σABC=(SBSCSCSASASB)\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. Parametro d'ordine chirale χ\chi: χABC=SA(SB×SC)\chi_{ABC} = \mathbf{S}_A \cdot (\mathbf{S}_B \times \mathbf{S}_C)

Metodi Numerici

Simulazioni Monte Carlo

  • Algoritmo: Metodo di tempering parallelo combinato con aggiornamenti locali heat bath
  • Dimensioni del Sistema: Dimensione lineare LL da 24 a 128 (N=3L2N = 3L^2 siti reticolari)
  • Media sul Disordine: Media su almeno 10 realizzazioni di disordine per ogni LL
  • Scelta dei Parametri: J1=1J_1 = -1, J3=1J_3 = 1

Analisi di Scaling a Dimensione Finita

Utilizzo delle seguenti relazioni per determinare gli esponenti critici della transizione di fase: CVmax(L)=H2H2NT2aLα/ν+bC_V^{\max}(L) = \frac{\langle H^2\rangle - \langle H\rangle^2}{NT^2} \simeq aL^{\alpha/\nu} + bSXmax(L)=X2X2NTcLγ/ν+dS_X^{\max}(L) = \frac{\langle X^2\rangle - \langle X\rangle^2}{NT} \simeq cL^{\gamma/\nu} + dTYmax(L)eL1/ν+TcT_{Y^{\max}}(L) \simeq eL^{-1/\nu} + T_c

Punti di Innovazione Tecnica

  1. Analisi del Meccanismo Microscopico: Rivelazione attraverso l'analisi di singole vacanze di come le vacanze producono campi magnetici efficaci, portando a orientamenti non colineari degli spin adiacenti
  2. Effetto di Campo Casuale: Chiarimento di come il disordine di sito produce un effetto di campo casuale efficace sull'ordine K4, che secondo l'argomento di Imry-Ma distrugge l'ordine discreto nei sistemi bidimensionali
  3. Meccanismo di Emergenza Chirale: Spiegazione di come i tre sottoreticoli tendono ad orientarsi mutuamente perpendicolari in presenza di vacanze, determinando così la chiralità locale

Configurazione Sperimentale

Parametri di Simulazione

  • Parametri di Scambio: J1=1J_1 = -1 (ferromagnetico), J3=1J_3 = 1 (antiferromagnetico)
  • Intervallo di Diluizione: 0p0.20 \leq p \leq 0.2
  • Intervallo di Temperatura: Regolato secondo il diagramma di fase, coprendo tutte le regioni di transizione di fase
  • Dimensioni del Sistema: L=24,36,48,64,80,96,128L = 24, 36, 48, 64, 80, 96, 128

Indicatori di Valutazione

  1. Calore Specifico CVC_V: Identifica la posizione della transizione di fase
  2. Suscettibilità del Parametro d'Ordine SΣS_\Sigma, SχS_\chi: Determina la natura della transizione di fase
  3. Parametri d'Ordine Σ\Sigma, χ\chi: Quantifica l'intensità dell'ordine

Elaborazione dei Dati

  • Termalizzazione adeguata e raccolta di statistiche per ogni punto di parametro
  • Media su più realizzazioni di disordine
  • Utilizzo dell'analisi di scaling a dimensione finita per estrarre gli esponenti critici

Risultati Sperimentali

Risultati Principali

Caratteristiche del Diagramma di Fase

È stato stabilito un diagramma di fase completo temperatura-tasso di diluizione, identificando tre regioni principali:

  1. Regione di Piccola Diluizione (p<0.05p < 0.05):
    • Transizione di fase K4 ad alta temperatura, con temperatura critica che diminuisce linearmente con pp
    • Emergenza della transizione di fase chirale a bassa temperatura, Tc4pT_c \approx 4p
  2. Regione di Diluizione Moderata (0.05p0.10.05 \leq p \leq 0.1):
    • Scomparsa dell'ordine K4
    • Persistenza dell'ordine chirale, con temperatura critica che aumenta con pp
  3. Regione di Grande Diluizione (0.1<p<0.20.1 < p < 0.2):
    • Presenza solo dell'ordine chirale
    • Scomparsa dell'ordine chirale anche a p=0.2p = 0.2

Verifica degli Esponenti Critici

  • Transizione di Fase K4: Appartiene alla classe di universalità di Potts a 4 stati (α/ν=1\alpha/\nu = 1, γ/ν=7/4\gamma/\nu = 7/4, 1/ν=3/21/\nu = 3/2)
  • Transizione di Fase Chirale: Appartiene alla classe di universalità di Ising (α/ν=0\alpha/\nu = 0, γ/ν=7/4\gamma/\nu = 7/4, 1/ν=11/\nu = 1)

Scoperte Chiave

  1. Analisi dell'Effetto di Vacanza: Una singola vacanza sul sottoreticolo B produce un campo magnetico efficace di intensità J1J_1, portando gli spin dei sottoreticoli A e C a tendere verso un'orientazione perpendicolare alla direzione dello spin mancante
  2. Meccanismo di Campo Casuale: Il disordine di sito produce un effetto di campo casuale efficace sull'ordine K4, con differenza di energia dipendente dalla posizione della vacanza: 8J3+J18J_3 + J_1 per siti dello stesso colore, 8J3J18J_3 - J_1 per siti di colore diverso
  3. Soglia di Percolazione: Quando p=0.2p = 0.2, la probabilità di triangoli completi (1p)3=0.512(1-p)^3 = 0.512 è prossima alla soglia di percolazione di sito del reticolo a nido d'ape 0.697

Lavori Correlati

Ricerca su Disordine che Induce Ordine

  1. Disordine Termico che Induce Ordine: Lavoro pionieristico di Villain e altri, dove il contributo entropico seleziona lo stato fondamentale con alta densità di eccitazioni a bassa energia
  2. Disordine Quantistico che Induce Ordine: Le fluttuazioni quantistiche selezionano le configurazioni di spin più colineari
  3. Effetti del Disordine Congelato: Tradizionalmente considerato come distruttore di ordine, questo articolo rivela per la prima volta i suoi effetti costruttivi

Ricerca sul Reticolo Kagome

  1. Modello Puro: Noto per possedere un ricco diagramma di fase a temperatura zero e fenomeni di ObTD
  2. Realizzazione Sperimentale: Il composto Ba-Vesignieite BaCu₃V₂O₈(OH)₂ rappresenta la controparte sperimentale di questo modello

Conclusioni e Discussione

Conclusioni Principali

  1. Prima osservazione dell'emergenza di ordine chirale guidato dal disordine congelato in un modello di Heisenberg con simmetria O(3)
  2. Costruzione di un diagramma di fase completo che rivela il processo di transizione dall'ordine K4 all'ordine chirale
  3. Proposizione di un nuovo meccanismo di "disordine congelato che induce ordine", distinto dai tradizionali meccanismi di disordine termico o quantistico che inducono ordine

Limitazioni

  1. Specificità del Modello: La ricerca è limitata al modello kagome J1J3J_1-J_3 specifico
  2. Limitazioni di Dimensione: A causa delle difficoltà di termalizzazione, le dimensioni del sistema nella regione di grande diluizione sono limitate a L80L \leq 80
  3. Intervallo di Temperatura: La precisione numerica nella regione a bassa temperatura è limitata

Direzioni Future

  1. Verifica di Universalità: Verifica del meccanismo in altri sistemi di magneti frustrati
  2. Verifica Sperimentale: Ricerca di sistemi sperimentali per verificare le previsioni teoriche
  3. Effetti Quantistici: Studio di fenomeni analoghi in sistemi di spin quantistici

Valutazione Approfondita

Vantaggi

  1. Innovazione Teorica: Scoperta di un nuovo meccanismo di disordine che induce ordine, sfidando la conoscenza convenzionale
  2. Rigore Numerico: Simulazioni Monte Carlo su larga scala, analisi di scaling a dimensione finita sufficientemente completa
  3. Chiarezza del Meccanismo: Spiegazione teorica completa dal microscopico al macroscopico
  4. Rilevanza Sperimentale: Connessione con sistemi di materiali reali

Insufficienze

  1. Specificità del Sistema: La universalità dei risultati richiede ulteriore verifica
  2. Limitazioni Numeriche: La precisione numerica in alcune regioni di parametri necessita miglioramento
  3. Verifica Sperimentale: Mancanza di verifica sperimentale diretta

Impatto

  1. Contributo Teorico: Fornisce una nuova prospettiva alla teoria dei magneti frustrati
  2. Valore Pratico: Fornisce principi guida per la progettazione di materiali
  3. Metodologia: Fornisce un esempio per lo studio di sistemi con disordine congelato

Scenari Applicabili

  1. Sistemi di Magneti Frustrati: In particolare sistemi con struttura multisottoreticolo
  2. Progettazione di Materiali: Realizzazione di ordine magnetico specifico attraverso il controllo dei difetti
  3. Ricerca Teorica: Fornisce riferimento per altri fenomeni di disordine che induce ordine

Bibliografia

Questo articolo cita 48 articoli correlati, coprendo importanti lavori in molteplici campi inclusi disordine che induce ordine, magneti frustrati, teoria delle transizioni di fase, fornendo una base teorica solida per la ricerca.


Impatto Ambientale Computazionale: Questo progetto ha utilizzato circa 40.000 giorni di core computazionali, equivalenti a circa 2.800 kg di CO₂, riflettendo la considerazione dei costi ambientali della ricerca numerica su larga scala.