✨ 要点🔬 技术摘要
在广袤的物理学领域中，存在着一类特殊的时刻，此时物质会从根本上改变其本质。这些是相变，就像水变成冰或铁变成磁铁那样的转变。在量子世界中，粒子表现得更像波而非固体物体，这些转变可能极其微妙。当一个系统被缓慢地推过这种变化时，它通常会遵循阻力最小的路径，保持在其最低能量状态。但如果推动得太快，系统就会掉队，无法跟上不断变化的情况。这种滞后产生了一种普遍的行为模式，一种仅取决于转变类型而非具体材料的宇宙指纹。科学家们称之为基布尔-祖雷克机制（Kibble-Zurek mechanism）。这是一个强大的概念，因为它表明，通过观察一个系统如何无法适应环境，我们可以了解它在完全静止时是如何表现的深层真相。然而，实时观察这些量子系统是极其困难的，尤其是当它们规模庞大且粒子处于“受挫”状态时——即粒子受到来自邻居的相互矛盾方向的拉扯。
来自洛桑联邦理工学院和弗拉蒂伦研究所的一个研究小组发现了一种观察这场戏剧上演的新方法。他们开发了一种计算方法，就像是量子系统的高速摄像机，能够模拟材料跨越临界阈值的确切时刻。他们并没有试图计算每一个单个粒子同时运动的、不可能实现的复杂性，而是构建了一个“基础”模型。想象一张精确描述系统在广泛条件下平静、低能状态的主图。研究人员训练了一个基于神经网络的高级计算机程序来创建这张图。一旦地图准备就绪，他们便从中选择了一组经过精心挑选的小规模状态作为固定的舞台。然后，他们观察系统在被驱动通过相变时是如何在这个舞台上移动的。因为舞台已经搭建完毕，他们可以利用不同的速度和系统规模运行数千次模拟，而无需每次都重新计算底层的物理规律。这种方法绕过了此前阻碍科学家研究二维系统动力学的常规瓶颈。
该团队首先在著名的横场伊辛模型（transverse-field Ising model）上测试了他们的方法，该模型描述了一个磁自旋网格。他们将系统推过临界点，并测量了剩余能量以及自旋之间的相关性。结果完美符合预期的普遍模式，证实了他们的方法能够高精度地捕捉基布尔-祖雷克机制。随后，他们将同样的技术应用于一个更困难的问题：受挫的正方晶格 J1-J2 海森堡模型（frustrated square-lattice J1-J2 Heisenberg model）。在这个系统中，磁相互作用相互竞争，形成了一场拉锯战，这种拉锯战会融化通常的磁有序，并可能产生一种被称为量子自旋液体（quantum spin liquid）的神秘状态。多年来，科学家们一直在争论从有序磁态到这种液体态的转变是平滑的连续变化，还是突然的跳跃。
通过在高达 16 x 16 个格点的晶格上模拟这个受挫系统，研究人员观察到了清晰的基布尔-祖雷克机制信号。系统对驱动速度的响应方式揭示了该转变确实是连续的。他们提取了描述临界行为的具体数值，发现这些数值与之前的静态计算高度吻合，但现在是通过动态视角得到了确认。这提供了强有力的、独立的证据，证明该转变是平滑的，从而解决了一个关于该量子相性质的长期争论。这项研究表明，曾经被认为只能在平衡态下测量的普遍标度律，也可以在实时转变的混沌奔流中被可靠地测量。通过将一个复杂的多元体问题转化为预建舞台上几个简单数字的移动，研究人员为探索此前难以触及的量子临界动力学开辟了一条实用的路径。
技术摘要：基于神经基础态子空间方法的二维受挫系统中的 Kibble-Zurek 动力学
问题陈述 J 1 J_1 J 1 J 2 J_2 J 2
方法论：神经基础态子空间 (NFS)
基础神经网络量子态 (FNQS)： 一个单一的神经网络（基于 Transformer 架构）被训练一次，用于表示涵盖整个淬火范围（从初始 γ 0 \gamma_0 γ 0 γ f \gamma_f γ f ∣ Φ θ ( γ ) ⟩ |\Phi_\theta(\gamma)\rangle ∣ Φ θ ( γ )⟩ 子空间构建： 从 FNQS 中提取一小组经过保真度筛选的基组 M M M ∣ ϕ k ⟩ = ∣ Φ θ ( γ k ) ⟩ |\phi_k\rangle = |\Phi_\theta(\gamma_k)\rangle ∣ ϕ k ⟩ = ∣ Φ θ ( γ k )⟩ 线性变分拟设： 时间演化的态可以表示为线性叠加：∣ Ψ ( t ) ⟩ = ∑ k = 1 M c k ( t ) ∣ ϕ k ⟩ |\Psi(t)\rangle = \sum_{k=1}^M c_k(t) |\phi_k\rangle ∣Ψ ( t )⟩ = ∑ k = 1 M c k ( t ) ∣ ϕ k ⟩ 投影动力学： 将随时间变化的薛定谔方程投影到这个固定的子空间上，从而得到关于系数 c ( t ) c(t) c ( t ) S c ˙ ( t ) = − i H ( γ ( t ) ) c ( t ) S \dot{c}(t) = -i H(\gamma(t)) c(t) S c ˙ ( t ) = − i H ( γ ( t )) c ( t )
S S S H ( γ ( t ) ) H(\gamma(t)) H ( γ ( t )) ⟨ ϕ i ∣ h ^ μ ∣ ϕ j ⟩ \langle \phi_i | \hat{h}_\mu | \phi_j \rangle ⟨ ϕ i ∣ h ^ μ ∣ ϕ j ⟩
效率： 由于基底是固定的且算符已预先计算，该方法在实时演化过程中避免了神经网络优化和采样。一旦建立子空间，不同的斜率时间（τ Q \tau_Q τ Q
关键结果
二维横场伊辛模型：
与 6 × 6 6 \times 6 6 × 6
有限尺寸缩放分析（6 ≤ L ≤ 12 6 \le L \le 12 6 ≤ L ≤ 12
提取的临界指数为：α = 0.362 ( 7 ) \alpha = 0.362(7) α = 0.362 ( 7 ) ν = 0.567 ( 18 ) \nu = 0.567(18) ν = 0.567 ( 18 ) η = 0.043 ( 29 ) \eta = 0.043(29) η = 0.043 ( 29 ) z = 1 z=1 z = 1
受挫方格点 J 1 J_1 J 1 J 2 J_2 J 2
应用于高达 16 × 16 16 \times 16 16 × 16 M = 28 M=28 M = 28 J 2 / J 1 = 0.49 J_2/J_1 = 0.49 J 2 / J 1 = 0.49
该方法为跨越 Néel 到自旋液体转变的 KZ 缩放提供了强有力的数值证据。
提取的指数为：ν = 1.23 ( 15 ) \nu = 1.23(15) ν = 1.23 ( 15 ) η = 0.409 ( 19 ) \eta = 0.409(19) η = 0.409 ( 19 ) z = 1 z=1 z = 1 ν ≈ 1.21 , η ≈ 0.410 \nu \approx 1.21, \eta \approx 0.410 ν ≈ 1.21 , η ≈ 0.410
意义与主张 J 1 J_1 J 1 J 2 J_2 J 2
作者强调，NFS 方法将计算负担转移到了单次基态优化上，避免了 tVMC 在时间演化过程中的不稳定性。他们指出，虽然该方法天然适用于慢速斜率和近临界动力学（即系统保持在基态流形附近时），但对于快速淬火或长时后临界演化（即状态探索了瞬时基态未捕捉到的方向时），其精度会下降。尽管如此，这项工作为在目前直接实时方法受到限制的二维受挫系统中进行动力学检验开辟了一条切实可行的路径。
每周获取最佳 condensed matter 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题，再试一次？
无垃圾邮件，随时退订。