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.
논문 ID : 2503.24139제목 : Chiral order emergence driven by quenched disorder저자 : Coraline Letouzé, Pascal Viot, Laura Messio (Sorbonne Université, CNRS, LPTMC)분류 : cond-mat.stat-mech, cond-mat.str-el발표 시간 : 2025년 9월 24일논문 링크 : https://arxiv.org/abs/2503.24139v3 담금질 무질서는 일반적으로 2차원 Ising 모델에 무작위 장을 가할 때와 같이 자기 질서를 파괴한다. 담금질 무질서가 존재하더라도 질서가 유지되는 경우, 이는 보통 순수 모델 질서의 잔존일 뿐이다. 본 논문은 놀라운 현상을 보여준다: 담금질 무질서에 의해 구동되는 질서의 출현. 이러한 질서는 순수 모델에 존재하는 질서와 전혀 공통점이 없다. 이러한 "무질서에 의한 질서화" 메커니즘은 통상적인 열역학적 또는 양자 무질서에 의한 질서화와 다르다. 병렬 담금질 몬테카를로 시뮬레이션을 통해 위치 희석을 갖는 카고메 격자 위의 고전적 J 1 − J 3 J_1-J_3 J 1 − J 3 하이젠베르그 모델을 연구했다. 소수의 공공(vacancy)이 기저 상태에 미치는 영향(비공면 구성을 선호함)을 분석한 후, 저온 카이랄 상의 출현과 공선 q = 4 q=4 q = 4 Potts 질서의 점진적 파괴를 보여준다.
전통적 인식의 도전 : 전통적으로 담금질 무질서는 다양한 질서 상태를 억제하거나 파괴한다고 여겨져 왔다. 본 연구는 이러한 전통적 관점에 도전하여 담금질 무질서가 완전히 새로운 질서 상태의 출현을 구동할 수 있음을 발견했다.이론적 의의 : 연속 스핀 자유도를 갖는 2차원 좌절 자기체에서, 비공선 또는 비공면 스핀 질서는 이산 대칭성(예: 시간 반전 또는 격자 대칭성)을 파괴할 수 있으며, Mermin-Wagner 정리에 따르면 이는 유한 온도 상전이를 허용한다.기존 한계 :기존의 "무질서에 의한 질서화" 메커니즘은 주로 열역학적 무질서에 의한 질서화(ObTD)와 양자 무질서에 의한 질서화(ObQD)를 포함하며, 이들은 보통 가장 공선적인 스핀 구성을 선택한다 담금질 무질서에 대한 연구는 주로 그 파괴적 작용에 집중되어 있으며, 그 건설적 작용에 대한 깊이 있는 이해가 부족하다 연구 동기 : 담금질 무질서가 새로운 질서 상태의 출현을 유도할 수 있는지, 특히 O(3) 대칭 모델에서의 카이랄 질서를 탐구하는 것은 이론적으로 중요한 의미를 갖는다.새로운 무질서에 의한 질서화 메커니즘 발견 : O(3) 스핀 회전 대칭성을 갖는 모델에서 담금질 무질서에 의해 구동되는 이산 카이랄 질서의 출현을 처음으로 관찰완전한 상도 구축 : 대규모 몬테카를로 시뮬레이션을 통해 위치 희석 조건 하에서 J 1 − J 3 J_1-J_3 J 1 − J 3 하이젠베르그 카고메 모델의 온도-희석율 상도 구성이론적 메커니즘 규명 : 공공이 유효 무작위 장을 생성하여 K4 질서를 파괴하면서 동시에 카이랄 질서 형성을 촉진하는 방식을 미시적 관점에서 설명보편성 클래스 검증 : 유한 크기 스케일링 분석을 통해 카이랄 상전이는 Ising 보편성 클래스에, K4 상전이는 4상태 Potts 보편성 클래스에 속함을 확인위치 희석 조건 하에서 카고메 격자 위의 J 1 − J 3 J_1-J_3 J 1 − J 3 하이젠베르그 모델의 상전이 거동, 특히 카이랄 질서의 출현 메커니즘을 연구한다.
위치 희석 조건 하의 해밀토니안은 다음과 같다:
H = J 1 ∑ ⟨ i , j ⟩ η i η j S i ⋅ S j + J 3 ∑ ⟨ i , j ⟩ 3 η i η j S i ⋅ S j 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 H = J 1 ∑ ⟨ i , j ⟩ η i η j S i ⋅ S j + J 3 ∑ ⟨ i , j ⟩ 3 η i η j S i ⋅ S j
여기서:
S i \mathbf{S}_i S i 는 고전적 3성분 단위 벡터η i = 1 \eta_i = 1 η i = 1 또는 0 0 0 은 위치 i i i 에서 자기 스핀의 존재 또는 부재를 나타냄희석율은 p = 1 − ∑ i η i / N p = 1 - \sum_i \eta_i/N p = 1 − ∑ i η i / N 으로 정의됨 K4 질서 매개변수 Σ \Sigma Σ : 삼각형 위의 국소 변수에 기반
σ A B C = ( S B ⋅ S C S C ⋅ S A S A ⋅ S B ) \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} σ A BC = S B ⋅ S C S C ⋅ S A S A ⋅ S B 카이랄 질서 매개변수 χ \chi χ :
χ A B C = S A ⋅ ( S B × S C ) \chi_{ABC} = \mathbf{S}_A \cdot (\mathbf{S}_B \times \mathbf{S}_C) χ A BC = S A ⋅ ( S B × S C ) 알고리즘 : 병렬 담금질 방법과 heat bath 국소 업데이트 결합시스템 크기 : 선형 크기 L L L 은 24에서 128까지 (N = 3 L 2 N = 3L^2 N = 3 L 2 격자점)무질서 평균 : 각 L L L 에 대해 최소 10개의 무질서 실현의 평균매개변수 선택 : J 1 = − 1 J_1 = -1 J 1 = − 1 , J 3 = 1 J_3 = 1 J 3 = 1 다음 관계식을 사용하여 상전이의 임계 지수를 결정:
C V max ( L ) = ⟨ H 2 ⟩ − ⟨ H ⟩ 2 N T 2 ≃ a L α / ν + b C_V^{\max}(L) = \frac{\langle H^2\rangle - \langle H\rangle^2}{NT^2} \simeq aL^{\alpha/\nu} + b C V m a x ( L ) = N T 2 ⟨ H 2 ⟩ − ⟨ H ⟩ 2 ≃ a L α / ν + b S X max ( L ) = ⟨ X 2 ⟩ − ⟨ X ⟩ 2 N T ≃ c L γ / ν + d S_X^{\max}(L) = \frac{\langle X^2\rangle - \langle X\rangle^2}{NT} \simeq cL^{\gamma/\nu} + d S X m a x ( L ) = NT ⟨ X 2 ⟩ − ⟨ X ⟩ 2 ≃ c L γ / ν + d T Y max ( L ) ≃ e L − 1 / ν + T c T_{Y^{\max}}(L) \simeq eL^{-1/\nu} + T_c T Y m a x ( L ) ≃ e L − 1/ ν + T c
미시적 메커니즘 분석 : 단일 공공 분석을 통해 공공이 어떻게 유효 자기장을 생성하여 인접 스핀의 비공선 배향을 초래하는지 규명무작위 장 효과 : 위치 무질서가 K4 질서에 어떻게 유효 무작위 장 효과를 생성하는지, 그리고 Imry-Ma 논증에 따라 2차원 시스템에서 이산 질서를 파괴하는 방식을 규명카이랄 출현 메커니즘 : 세 개의 부분 격자가 공공 존재 하에서 어떻게 상호 수직 배향으로 향하여 국소 카이랄성을 결정하는지 설명교환 매개변수 : J 1 = − 1 J_1 = -1 J 1 = − 1 (강자성), J 3 = 1 J_3 = 1 J 3 = 1 (반강자성)희석 범위 : 0 ≤ p ≤ 0.2 0 \leq p \leq 0.2 0 ≤ p ≤ 0.2 온도 범위 : 상도에 따라 조정, 모든 상전이 영역 포함시스템 크기 : L = 24 , 36 , 48 , 64 , 80 , 96 , 128 L = 24, 36, 48, 64, 80, 96, 128 L = 24 , 36 , 48 , 64 , 80 , 96 , 128 비열 C V C_V C V : 상전이 위치 식별질서 매개변수 자기 감수율 S Σ S_\Sigma S Σ , S χ S_\chi S χ : 상전이 성질 결정질서 매개변수 Σ \Sigma Σ , χ \chi χ : 질서 강도 정량화각 매개변수 지점에서 충분한 열화 및 통계 수행 여러 무질서 실현에 대한 평균 유한 크기 스케일링 분석을 사용하여 임계 지수 추출 온도-희석율 완전 상도를 구축하여 세 가지 주요 영역을 식별:
낮은 희석 영역 (p < 0.05 p < 0.05 p < 0.05 ):고온 K4 질서 상전이, 임계 온도는 p p p 에 따라 선형 감소 저온 카이랄 상전이 출현, T c ≈ 4 p T_c \approx 4p T c ≈ 4 p 중간 희석 영역 (0.05 ≤ p ≤ 0.1 0.05 \leq p \leq 0.1 0.05 ≤ p ≤ 0.1 ):K4 질서 소멸 카이랄 질서 지속, 임계 온도는 p p p 증가에 따라 증가 높은 희석 영역 (0.1 < p < 0.2 0.1 < p < 0.2 0.1 < p < 0.2 ):카이랄 질서만 존재 p = 0.2 p = 0.2 p = 0.2 일 때 카이랄 질서도 소멸K4 상전이 : 4상태 Potts 보편성 클래스에 속함 (α / ν = 1 \alpha/\nu = 1 α / ν = 1 , γ / ν = 7 / 4 \gamma/\nu = 7/4 γ / ν = 7/4 , 1 / ν = 3 / 2 1/\nu = 3/2 1/ ν = 3/2 )카이랄 상전이 : Ising 보편성 클래스에 속함 (α / ν = 0 \alpha/\nu = 0 α / ν = 0 , γ / ν = 7 / 4 \gamma/\nu = 7/4 γ / ν = 7/4 , 1 / ν = 1 1/\nu = 1 1/ ν = 1 )공공 효과 분석 : B 부분 격자 위의 단일 공공이 생성하는 유효 자기장 강도는 J 1 J_1 J 1 이며, A, C 부분 격자 스핀이 결실 스핀 방향에 수직으로 향하도록 함무작위 장 메커니즘 : 위치 무질서가 K4 질서에 유효 무작위 장을 생성하며, 에너지 차이는 공공 위치에 따라 달라짐: 같은 색 위치는 8 J 3 + J 1 8J_3 + J_1 8 J 3 + J 1 , 다른 색은 8 J 3 − J 1 8J_3 - J_1 8 J 3 − J 1 침투 임계값 : p = 0.2 p = 0.2 p = 0.2 일 때, 완전한 삼각형의 확률 ( 1 − p ) 3 = 0.512 (1-p)^3 = 0.512 ( 1 − p ) 3 = 0.512 는 벌집 격자의 위치 침투 임계값 0.697에 가까움열역학적 무질서에 의한 질서화 : Villain 등의 개척적 연구, 엔트로피 기여가 높은 밀도 낮은 에너지 여기의 기저 상태를 선택양자 무질서에 의한 질서화 : 양자 변동이 가장 공선적인 스핀 구성을 선택담금질 무질서 효과 : 전통적으로 질서를 파괴한다고 여겨졌으나, 본 논문에서 처음으로 그 건설적 작용 발견순수 모델 : 풍부한 영온 상도 및 ObTD 현상이 알려짐실험적 실현 : Ba-Vesignieite 화합물 BaCu₃V₂O₈(OH)₂가 이 모델의 실험적 대응물O(3) 대칭 하이젠베르그 모델에서 담금질 무질서에 의해 구동되는 카이랄 질서의 출현을 처음으로 관찰 완전한 상도를 구축하여 K4 질서에서 카이랄 질서로의 전환 과정 규명 전통적인 열역학적 또는 양자 무질서에 의한 질서화와 구별되는 "담금질 무질서에 의한 질서화"의 새로운 메커니즘 제시 모델 특수성 : 연구는 특정 J 1 − J 3 J_1-J_3 J 1 − J 3 카고메 모델로 제한됨크기 제한 : 열화의 어려움으로 인해 높은 희석 영역의 시스템 크기는 L ≤ 80 L \leq 80 L ≤ 80 으로 제한됨온도 범위 : 저온 영역의 수치 정확도가 제한됨보편성 검증 : 다른 좌절 자기체 시스템에서 이 메커니즘 검증실험적 검증 : 이론 예측을 검증할 실험 시스템 탐색양자 효과 : 양자 스핀 시스템에서의 유사 현상 연구이론적 혁신 : 전통적 인식에 도전하는 완전히 새로운 무질서에 의한 질서화 메커니즘 발견수치적 엄밀성 : 대규모 몬테카를로 시뮬레이션, 충분한 유한 크기 스케일링 분석메커니즘 명확성 : 미시적에서 거시적 수준의 완전한 이론적 설명실험 관련성 : 실제 물질 시스템과의 연결시스템 특수성 : 결과의 보편성은 추가 검증 필요수치 제한 : 특정 매개변수 영역의 수치 정확도 개선 필요실험적 검증 : 직접적인 실험적 검증 부재이론적 기여 : 좌절 자기체 이론에 새로운 관점 제공실용적 가치 : 재료 설계에 지침 원리 제공방법론 : 담금질 무질서 시스템 연구를 위한 범례 제시좌절 자기체 시스템 : 특히 다중 부분 격자 구조를 갖는 시스템재료 설계 : 결함 제어를 통한 특정 자기 질서 실현이론 연구 : 다른 무질서에 의한 질서화 현상의 참고 자료본 논문은 48편의 관련 문헌을 인용하며, 무질서에 의한 질서화, 좌절 자기체, 상전이 이론 등 여러 분야의 중요한 연구를 포함하여 연구의 견고한 이론적 기초를 제공한다.
계산 환경 영향 : 본 프로젝트는 약 40,000 코어일의 계산 자원을 사용하였으며, 이는 약 2,800kg CO₂ 상당의 탄소 배출에 해당하여 대규모 수치 연구의 환경 비용 고려를 반영한다.