基于纯化的变分方法已经被广泛应用于研究多体系统热平衡态的理论计算中。但是，在处理低温问题，特别是系统热态存在自发对称性破缺时，基于纯化的方法存在着一些理论问题，并不能给出体系正确的热态。

　　为了克服这一困难，中国科学院理论物理研究所副研究员石弢与美国哈佛大学、德国马普量子光学研究所研究人员合作，发展了一套系统研究多体系统热态的变分方法。在这个方法中，他们提出了一套全新的虚实演化方程，使得自由能单调递减，从而得到变分流形上的不动点，对应于系统近似的热平衡态。基于对非高斯态的理解，他们提出了10种非高斯变换，并在解析上给出了所有变分参数满足的运动方程。通过数值求解运动方程，可以高效地得到系统热态。研究结果发表在Physical Review Letters上。