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 (ソルボンヌ大学、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 ハイゼンベルグモデルを研究した。少数の空孔が基底状態に与える影響(非共面構成を有利にする)の分析後、低温カイラル相の涌現と共線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次元系で離散秩序を破壊する方法を明確にしたカイラル涌現機構 : 3つの副格子が空孔存在下で相互に垂直な配向を示す傾向があり、局所カイラル性を決定する方法を説明した交換パラメータ : 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 χ :秩序強度を定量化各パラメータ点について十分な熱化と統計処理を実施 複数の無序実現に対して平均化 有限サイズスケーリング分析を使用して臨界指数を抽出 完全な温度-希釈率相図を確立し、3つの主要領域を識別:
低希釈領域 (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₂相当の炭素排出に相当し、大規模数値研究の環境コスト考慮を反映している。