[论文解读] Phase transition in the three-dimensional $O(N)\otimes O(M)$ model: a Monte Carlo study

A. O. Sorokin|arXiv (Cornell University)|May 15, 2022

Theoretical and Computational Physics被引用 4

一句话总结

本蒙特卡洛研究调查了三维 $O(N)\otimes O(M)$ 模型在 $2 \leq M \leq 4$ 且 $M \leq N \leq 8$ 条件下的相变行为，发现当 $N \geq M+4$ 时为连续相变。该工作通过表明临界维度 $N_c^+(M,3)$ 位于一阶与二阶大 $N$ 展开估计之间，解决了先前理论预测中的不一致之处，且在 $M=2$ 时各类方法得到的临界指数一致，并在拓扑缺陷缺失的情况下识别出二阶相变。

ABSTRACT

Using Monte Carlo simulations, we consider the lattice version of the $O(N)\otimes O(M)$ sigma model for $2\leq M\leq4$ and $M\leq N \leq8$. We find a continuous transition for $N\geq M+4$. Estimates of the critical exponents for cases of second-order and weak first-order transitions are found. For $M=2$ our estimates of the exponents and marginal dimensionality $N_c^+(M)$ are in good agreement with the results of the non-perturbative renormalization group approach. For $M\geq2$ we find estimates of the exponents and marginal dimensionality between the values obtained in the first and second orders of the large-N expansion. To complete the picture, we also consider the usual $O(N)$ model ($M=1$).

研究动机与目标

• 解决三维 $O(N)\otimes O(M)$ 模型临界行为理论预测中长期存在的不一致问题。
• 通过有限尺寸标度法，确定不同 $N$ 和 $M$ 值下的相变阶数（二阶、弱一阶或一阶）。
• 估算临界指数和临界维度 $N_c^+(M,3)$，特别是针对 $M=2,3,4$，并与大 $N$ 展开、非微扰重整化群（RG）和共形 bootstrap 方法的结果进行比较。
• 通过识别拓扑缺陷缺失的情况，阐明拓扑缺陷在决定相变阶数中的作用，此时相变仍保持连续。

提出的方法

• 在立方格点上对格点 $O(N)\otimes O(M)$ 胞射模型进行数值蒙特卡洛模拟，格点尺寸最大达 $L=60$。
• 通过内能分布的有限尺寸标度分析，检测可能指示一阶相变的双峰结构。
• 通过能量和比热等可观测量的标度塌缩法估算临界指数。
• 计算费舍尔指数 $\eta$，以区分弱一阶相变（$\eta < 0$）与二阶相变（$\eta \geq 0$）。
• 将临界温度的倒数 $K_c$ 作为 $N$ 的函数，对每个 $M$ 拟合经验标度关系。
• 将结果与大 $N$ 展开、非微扰 RG 和共形 bootstrap 的理论预测进行比较，尤其针对 $M=2$ 的情况。

实验结果

研究问题

• RQ1 三维 $O(N)\otimes O(M)$ 模型是否在 $N \geq M+4$ 时表现出连续相变，如某些 RG 方法所预测？
• RQ2 临界维度 $N_c^+(M,3)$ 的精确值是多少？其与 $4-\varepsilon$、大 $N$ 和非微扰 RG 方法的估计值相比如何？
• RQ3 为何微扰 RG 方法对 $M=2$ 且 $N=2,3,4$ 预测存在焦点型固定点，而其他方法未发现？
• RQ4 数值模拟能否解决 $O(2)\otimes O(2)$ 和 $O(3)\otimes O(2)$ 模型中相变性质的争议？
• RQ5 在何种条件下——尤其是关于拓扑缺陷——尽管存在多个耦合常数，相变仍保持为二阶？

主要发现

• 对于 $M=2$，估算的临界指数和临界维度 $N_c^+(2,3)$ 与非微扰重整化群（NPRG）结果高度一致。
• 对于 $M \geq 2$，估算的临界指数和 $N_c^+(M,3)$ 位于一阶与二阶大 $N$ 展开估计之间，表明重整化后收敛性得到改善。
• 当 $N \geq M+4$ 时观察到二阶相变，且在 $L=60$ 范围内内能分布无双峰结构。
• 对于 $M=2$，在 $N=4$ 和 $N=5$ 时为弱一阶相变（$\eta < 0$），在 $N=6,7,8$ 时为二阶相变；对于 $M=3$，在 $N=7,8$ 时为二阶相变；对于 $M=4$，在 $N=8$ 时为二阶相变。
• 临界温度倒数 $K_c$ 对每个 $M$ 均与 $N$ 呈线性关系，拟合结果为：$K_c \approx 0.244812N - 0.082677$（$M=2$），$K_c \approx 0.245994N - 0.135670$（$M=3$），$K_c \approx 0.250345N - 0.169116$（$M=4$）。
• 本研究未发现 $M=2$ 且 $N=2,3,4$ 时存在焦点型固定点，与微扰 RG 预测相矛盾，反而支持 $N \geq M+4$ 时为连续相变。