Ising模型暂态定义及transient_results函数功能与模拟影响问询

背景说明 #

我下载了一个基于C#的简单Ising模型实现，大致理解了全部代码，但对以下函数逻辑存在疑问：

private static void transient_results(double T)
{
    for (int a = 1; a <= transient; a++)
    {
        array_to_list();
        for (int b = 1; b <= L * L; b++)
        {
            choose = choose_random_site("i", 0);
            posx = choose_random_site("x", choose);
            posy = choose_random_site("y", choose);

            if (test_flip(posx,posy,T))
            {
                flip(posx,posy);
            }

            list.RemoveAt(choose);
        }
    }
}

• transient_results()接收温度值T作为参数。
• transient是从控制台读取的整数，代表暂态迭代次数。
• array_to_list()会清空字符串列表，并以"i , j"格式的新字符串初始化该列表，用作位点定位器，跟踪已处理/未处理的位点位置。
• 代码块choose = choose_random_site("i", 0); posx = choose_random_site("x", choose); posy = choose_random_site("y", choose);用于随机选择一个位点及其对应的(x, y)坐标。
• test_flip()用于检查位点状态是否可翻转，返回布尔值；若位点可翻转，则执行flip(posx,posy)完成翻转。
• 无论位点是否可翻转，其对应的位点定位器字符串都会从列表中移除，标记该位点已处理。

技术问题解答 #

1. Ising模型中暂态（transient state）或暂态阶段（transient phase）的定义是什么？ #

在Ising模型模拟里，暂态指的是系统从初始状态（比如随机自旋分布）出发，还没达到热力学平衡态的阶段。这个阶段里，系统的宏观物理量（比如磁化强度、能量）会随着模拟步数变化，还没稳定在对应温度下的平衡值。简单说就是系统"热身"的过程，还没进入稳定的可观测状态。

2. 上述transient_results函数整体实现了什么功能？ #

这个函数就是实现Ising模型的暂态演化过程：

• 外层循环执行transient次完整的"遍历所有位点"操作；
• 每次外层循环开始时，重新初始化所有位点的待处理列表；
• 内层循环会随机遍历系统的每一个位点（每个位点恰好被处理一次），对每个位点尝试用Metropolis准则（test_flip应该是实现这个准则）判断是否翻转自旋，无论是否翻转都标记该位点已处理；
• 整个过程就是让系统在指定温度T下，经过transient轮的全位点随机更新，远离初始状态，向平衡态靠近。

3. 该函数会对Ising模型的模拟过程及结果产生怎样的影响？ #

• 模拟过程：它会占用一部分模拟时间，让系统先完成"热身"，之后再进行平衡态的采样（比如统计磁化强度、能量等），避免初始状态的随机性干扰平衡态数据；
• 结果准确性：如果跳过这个过程，直接从初始状态开始采样，得到的物理量会包含初始态的偏差，无法反映对应温度下的真实平衡性质；
• 参数影响：transient的取值很关键——取值太小，系统没达到平衡，采样数据不准；取值太大则会浪费计算资源。一般需要根据系统大小、温度（比如接近相变温度时需要更长暂态）调整这个参数。

内容的提问来源于stack exchange，提问作者user366312