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

Chiral order emergence driven by quenched disorder

基本信息

  • 论文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晶格上具有位点稀释的经典J1J3J_1-J_3海森堡模型。在分析少数空位对基态的影响(有利于非共面构型)后,展示了低温手性相的涌现和共线q=4q=4 Potts有序的渐进破坏。

研究背景与动机

问题的重要性

  1. 传统认知的挑战: 传统上,淬火无序被认为会抑制或破坏各种有序态。本研究挑战了这一传统观点,发现淬火无序可以驱动全新有序态的涌现。
  2. 理论意义: 在具有连续自旋自由度的二维受挫磁体中,非共线或非共面自旋有序可能破坏离散对称性(如时间反演或晶格对称性),根据Mermin-Wagner定理,这允许有限温度相变。
  3. 现有局限性:
    • 现有的"无序致序"机制主要包括热力学无序致序(ObTD)和量子无序致序(ObQD),它们通常选择最共线的自旋构型
    • 淬火无序的研究主要集中在其破坏性作用上,缺乏对其建构性作用的深入理解
  4. 研究动机: 探索淬火无序是否能够诱导新有序态的涌现,特别是在O(3)对称模型中的手性有序,这在理论上具有重要意义。

核心贡献

  1. 发现新的无序致序机制: 首次在具有O(3)自旋旋转对称性的模型中观察到由淬火无序驱动的离散手性有序涌现
  2. 建立完整相图: 通过大规模蒙特卡罗模拟,构建了J1J3J_1-J_3海森堡kagome模型在位点稀释下的温度-稀释率相图
  3. 理论机制阐释: 从微观角度解释了空位如何通过产生有效随机场破坏K4有序,同时促进手性有序的形成
  4. 普适性类验证: 通过有限尺寸标度分析,确认了手性相变属于Ising普适性类,而K4相变属于4态Potts普适性类

方法详解

任务定义

研究kagome晶格上J1J3J_1-J_3海森堡模型在位点稀释条件下的相变行为,特别关注手性有序的涌现机制。

模型架构

哈密顿量

位点稀释条件下的哈密顿量为: 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

其中:

  • Si\mathbf{S}_i为经典三分量单位向量
  • ηi=1\eta_i = 100表示位点ii上磁性自旋的存在或缺失
  • 稀释率定义为p=1iηi/Np = 1 - \sum_i \eta_i/N

序参量定义

  1. K4序参量 Σ\Sigma:基于三角形上的局域变量 σ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. 手性序参量 χ\chiχABC=SA(SB×SC)\chi_{ABC} = \mathbf{S}_A \cdot (\mathbf{S}_B \times \mathbf{S}_C)

数值方法

蒙特卡罗模拟

  • 算法: 并行回火方法结合heat bath局域更新
  • 系统尺寸: 线性尺寸LL从24到128 (N=3L2N = 3L^2个格点)
  • 无序平均: 每个LL至少10个无序实现的平均
  • 参数选择: J1=1J_1 = -1, J3=1J_3 = 1

有限尺寸标度分析

使用以下关系确定相变的临界指数: 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

技术创新点

  1. 微观机制分析: 通过单空位分析揭示了空位如何产生有效磁场,导致相邻自旋的非共线取向
  2. 随机场效应: 阐明了位点无序如何对K4有序产生有效随机场效应,根据Imry-Ma论证在二维系统中破坏离散有序
  3. 手性涌现机制: 解释了三个子晶格如何在空位存在下趋向相互垂直取向,从而确定局域手性

实验设置

模拟参数

  • 交换参数: J1=1J_1 = -1 (铁磁), J3=1J_3 = 1 (反铁磁)
  • 稀释范围: 0p0.20 \leq p \leq 0.2
  • 温度范围: 根据相图调整,覆盖所有相变区域
  • 系统尺寸: L=24,36,48,64,80,96,128L = 24, 36, 48, 64, 80, 96, 128

评价指标

  1. 比热容 CVC_V:标识相变位置
  2. 序参量磁化率 SΣS_\Sigma, SχS_\chi:确定相变性质
  3. 序参量 Σ\Sigma, χ\chi:量化有序强度

数据处理

  • 每个参数点进行充分的热化和统计
  • 对多个无序实现进行平均
  • 使用有限尺寸标度分析提取临界指数

实验结果

主要结果

相图特征

建立了完整的温度-稀释率相图,识别出三个主要区域:

  1. 小稀释区域 (p<0.05p < 0.05):
    • 高温K4有序相变,临界温度随pp线性下降
    • 低温手性相变涌现,Tc4pT_c \approx 4p
  2. 中等稀释区域 (0.05p0.10.05 \leq p \leq 0.1):
    • K4有序消失
    • 手性有序持续存在,临界温度随pp增加
  3. 大稀释区域 (0.1<p<0.20.1 < p < 0.2):
    • 仅存在手性有序
    • p=0.2p = 0.2时手性有序也消失

临界指数验证

  • K4相变: 属于4态Potts普适性类 (α/ν=1\alpha/\nu = 1, γ/ν=7/4\gamma/\nu = 7/4, 1/ν=3/21/\nu = 3/2)
  • 手性相变: 属于Ising普适性类 (α/ν=0\alpha/\nu = 0, γ/ν=7/4\gamma/\nu = 7/4, 1/ν=11/\nu = 1)

关键发现

  1. 空位效应分析: 单个空位在B子晶格上产生的有效磁场强度为J1J_1,导致A、C子晶格自旋趋向垂直于缺失自旋方向
  2. 随机场机制: 位点无序对K4有序产生有效随机场,能量差依赖于空位位置:相同颜色位点为8J3+J18J_3 + J_1,不同颜色为8J3J18J_3 - J_1
  3. 渗流阈值: 当p=0.2p = 0.2时,完整三角形的概率(1p)3=0.512(1-p)^3 = 0.512接近蜂窝晶格的位点渗流阈值0.697

相关工作

无序致序研究

  1. 热力学无序致序: Villain等人的开创性工作,熵贡献选择高密度低能激发的基态
  2. 量子无序致序: 量子涨落选择最共线自旋构型
  3. 淬火无序效应: 传统上认为破坏有序,本文首次发现其建构性作用

kagome晶格研究

  1. 纯净模型: 已知存在丰富的零温相图和ObTD现象
  2. 实验实现: Ba-Vesignieite化合物BaCu₃V₂O₈(OH)₂为该模型的实验对应物

结论与讨论

主要结论

  1. 首次在O(3)对称海森堡模型中观察到由淬火无序驱动的手性有序涌现
  2. 建立了完整的相图,揭示了K4有序到手性有序的转变过程
  3. 提出了"淬火无序致序"的新机制,区别于传统的热力学或量子无序致序

局限性

  1. 模型特殊性: 研究局限于特定的J1J3J_1-J_3kagome模型
  2. 尺寸限制: 由于热化困难,大稀释区域的系统尺寸限制在L80L \leq 80
  3. 温度范围: 低温区域的数值精度受到限制

未来方向

  1. 普适性验证: 在其他受挫磁体系统中验证该机制
  2. 实验验证: 寻找实验系统验证理论预测
  3. 量子效应: 研究量子自旋系统中的类似现象

深度评价

优点

  1. 理论创新: 发现了全新的无序致序机制,挑战了传统认知
  2. 数值严谨: 大规模蒙特卡罗模拟,充分的有限尺寸标度分析
  3. 机制清晰: 从微观到宏观的完整理论解释
  4. 实验相关: 与真实材料系统的联系

不足

  1. 系统特殊性: 结果的普适性需要进一步验证
  2. 数值限制: 某些参数区域的数值精度有待提高
  3. 实验验证: 缺乏直接的实验验证

影响力

  1. 理论贡献: 为受挫磁体理论提供新视角
  2. 实用价值: 为材料设计提供指导原则
  3. 方法学: 为研究淬火无序系统提供范例

适用场景

  1. 受挫磁体系统: 特别是具有多子晶格结构的系统
  2. 材料设计: 通过控制缺陷实现特定磁有序
  3. 理论研究: 为其他无序致序现象提供参考

参考文献

本文引用了48篇相关文献,涵盖了无序致序、受挫磁体、相变理论等多个领域的重要工作,为研究提供了坚实的理论基础。


计算环境影响: 该项目使用约40,000核心日的计算资源,相当于2,800kg CO₂当量的碳排放,体现了大规模数值研究的环境成本考虑。