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发表时间 : September 24, 2025论文链接 : https://arxiv.org/abs/2503.24139v3 淬火无序通常会破坏磁有序,例如在二维Ising模型中施加随机场时。即使在淬火无序存在下仍有有序存在,通常也只是干净模型有序的残存。本文展示了一个令人惊讶的现象:由淬火无序驱动的有序涌现。这种有序与干净模型中存在的有序毫无共同之处。这种"无序致序"机制不同于通常的热力学或量子无序致序。通过并行回火蒙特卡罗模拟,研究了kagome晶格上具有位点稀释的经典J 1 − J 3 J_1-J_3 J 1 − J 3 海森堡模型。在分析少数空位对基态的影响(有利于非共面构型)后,展示了低温手性相的涌现和共线q = 4 q=4 q = 4 Potts有序的渐进破坏。
传统认知的挑战 : 传统上,淬火无序被认为会抑制或破坏各种有序态。本研究挑战了这一传统观点,发现淬火无序可以驱动全新有序态的涌现。理论意义 : 在具有连续自旋自由度的二维受挫磁体中,非共线或非共面自旋有序可能破坏离散对称性(如时间反演或晶格对称性),根据Mermin-Wagner定理,这允许有限温度相变。现有局限性 :现有的"无序致序"机制主要包括热力学无序致序(ObTD)和量子无序致序(ObQD),它们通常选择最共线的自旋构型 淬火无序的研究主要集中在其破坏性作用上,缺乏对其建构性作用的深入理解 研究动机 : 探索淬火无序是否能够诱导新有序态的涌现,特别是在O(3)对称模型中的手性有序,这在理论上具有重要意义。发现新的无序致序机制 : 首次在具有O(3)自旋旋转对称性的模型中观察到由淬火无序驱动的离散手性有序涌现建立完整相图 : 通过大规模蒙特卡罗模拟,构建了J 1 − J 3 J_1-J_3 J 1 − J 3 海森堡kagome模型在位点稀释下的温度-稀释率相图理论机制阐释 : 从微观角度解释了空位如何通过产生有效随机场破坏K4有序,同时促进手性有序的形成普适性类验证 : 通过有限尺寸标度分析,确认了手性相变属于Ising普适性类,而K4相变属于4态Potts普适性类研究kagome晶格上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 为经典三分量单位向量η 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论证在二维系统中破坏离散有序手性涌现机制 : 解释了三个子晶格如何在空位存在下趋向相互垂直取向,从而确定局域手性交换参数 : 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 kagome模型尺寸限制 : 由于热化困难,大稀释区域的系统尺寸限制在L ≤ 80 L \leq 80 L ≤ 80 温度范围 : 低温区域的数值精度受到限制普适性验证 : 在其他受挫磁体系统中验证该机制实验验证 : 寻找实验系统验证理论预测量子效应 : 研究量子自旋系统中的类似现象理论创新 : 发现了全新的无序致序机制,挑战了传统认知数值严谨 : 大规模蒙特卡罗模拟,充分的有限尺寸标度分析机制清晰 : 从微观到宏观的完整理论解释实验相关 : 与真实材料系统的联系系统特殊性 : 结果的普适性需要进一步验证数值限制 : 某些参数区域的数值精度有待提高实验验证 : 缺乏直接的实验验证理论贡献 : 为受挫磁体理论提供新视角实用价值 : 为材料设计提供指导原则方法学 : 为研究淬火无序系统提供范例受挫磁体系统 : 特别是具有多子晶格结构的系统材料设计 : 通过控制缺陷实现特定磁有序理论研究 : 为其他无序致序现象提供参考本文引用了48篇相关文献,涵盖了无序致序、受挫磁体、相变理论等多个领域的重要工作,为研究提供了坚实的理论基础。
计算环境影响 : 该项目使用约40,000核心日的计算资源,相当于2,800kg CO₂当量的碳排放,体现了大规模数值研究的环境成本考虑。