密度矩阵重整化群

密度矩阵重整化群 (Density Matrix Renormalization Group)，简称DMRG，是一种数值算法，于公元1992年由美国物理学家史提芬·怀特提出。密度矩阵重整化群是用来计算量子多体系统（例如：Hubbard model、t-J模型、海森堡模型，等等）的一个非常精准的数值算法，在一维或准一维的系统可以得到系统尺寸很大且很准确的计算结果，但是在二维的量子多体系统中却很难达到所需要的精确度。目前此算法仍无法计算三维的量子系统。

从数值计算的角度来看，量子多体物理主要的困难之处就在于系统的希尔伯特空间维度随着系统的尺寸呈指数成长，例如，一个由${\displaystyle N}$个自旋1/2的粒子所组成的一维晶格系统其希尔伯特空间维度大小为 ${\displaystyle 2^N}$。 传统的解决方法有两种：

史提芬·怀特最先意识到，NRG在演算Hubbard模型中的失败，是由于在NRG的迭代过程中忽略了环境对系统的影响。换句话说，NRG的重整化方法——只保留低能量本征态——并不能正确得出下一次迭代时的低能状态。
DMRG的重整化方法不同于NRG。DMRG在重整化前，把整个系统视为两个部分，一部分为系统，一部分为环境，而系统和环境的整体称为超块。接着，计算超块的基态，有了基态之后便计算约化密度矩阵，然后对角化这个约化密度矩阵，选出拥有较大的本征值的本征态。这些拥有较大的本征值的本征态正是基态性质最重要的态，然后根据此标准对系统部分做重整化。

实际实行DMRG是一个很冗长的工作，一些主要常用的计算手段如下：

如缺少上述的一些计算手段，DMRG可能难以完成对实际物理模型的演算。

DMRG 已经成功的在许多不同的一维模型上计算低能态的一些性质，如易辛模型，海森保模型等自旋模型，费米子系统如 Hubbard 模型 ，杂质系统如近藤效应，玻色子系统，混合玻色子与费米子的系统。随着现代电脑硬件技术的进步，DMRG应用在二维系统上可行性愈来愈高，目前一般的作法是将二维系统视为一个多腿的梯子，再将梯子的长度拉长。2011年发表在《Science》封面的一篇文章中，利用 DMRG 探讨二维Kagome晶格中的自旋-1/2系统的基态。由这篇文章来看， DMRG 可能仍是对付二维系统最强大的武器。

DMRG之所以在一维系统中如此成功，背后的理论可以用矩阵积态来加以解释。有限尺度的DMRG中，扫荡的过程等同于将此系统的波函数写在矩阵积态空间做变分法。以自旋-1/2的系统为例，矩阵积态如以下形式：

${\displaystyle |\mathrm{\Phi <!--\; \Phi \; -->}⟩<!--\; ⟩\; -->=\underset{{\mathrm{\sigma <!--\; \sigma \; -->}}_1\cdots <!--\; \cdots \; -->{\mathrm{\sigma <!--\; \sigma \; -->}}_N}{\sum <!--\; \sum \; -->}(A_1^{}{A}_2^{}\cdots <!--\; \cdots \; -->A_n^{}\cdots <!--\; \cdots \; -->A_{N-<!--\; -\; -->1}^{}{A}_N^{})|{\mathrm{\sigma <!--\; \sigma \; -->}}_1{\mathrm{\sigma <!--\; \sigma \; -->}}_2\cdots <!--\; \cdots \; -->{\mathrm{\sigma <!--\; \sigma \; -->}}_N⟩<!--\; ⟩\; -->}$

其中${\displaystyle {\mathrm{\sigma <!--\; \sigma \; -->}}_1{\mathrm{\sigma <!--\; \sigma \; -->}}_2\cdots <!--\; \cdots \; -->{\mathrm{\sigma <!--\; \sigma \; -->}}_N}$表示每一个格点上自旋${\displaystyle z}$方向的分量，${\displaystyle A_i^{}}$表示第${\displaystyle i}$格点、自旋${\displaystyle z}$方向的分量为${\displaystyle {\mathrm{\sigma <!--\; \sigma \; -->}}_i}$的矩阵。${\displaystyle A_1^{}}$矩阵大小是1×d、${\displaystyle A_2^{}}$矩阵大小是d×d2、${\displaystyle A_3^{}}$矩阵大小是d2×d3、……直到第${\displaystyle n}$格点时，dn≥m，${\displaystyle A_n^{}}$矩阵大小是dn-1×m、${\displaystyle A_{n+1}^{}}$矩阵大小是m×m、……，${\displaystyle A_{N-<!--\; -\; -->1}^{}}$矩阵大小是d2×d、${\displaystyle A_N^{}}$矩阵大小是d×1。当m趋近无穷大时，所有的波函数皆可写成矩阵积态的形式。

DMRG的巨大成功带给人们许多冲击与启示，可惜的是由于波函数被表示成矩阵积态（Matrix Product State），造成DMRG在处理二维量子晶格系统时特别困难，更别说是三维的量子系统。继承DMRG的知识和技术，许多物理学家着手发展适合研究二维甚至三维系统中的数值方法，例如：TEBD（Time-evolving block decimation）、PEPS（Projected Entangled Pair States）、MERA（multi-scale entanglement renormalization ansatz），等等。另一方面，也有许多物理学家在原有的DMRG方法上加以改良，让科学家可以处理更多有趣的一维量子晶格系统的问题，例如：时间演化、有限温度，等等。