来源论文: https://arxiv.org/abs/2605.21585v1 生成时间: Jun 14, 2026 10:31

强关联体系中交织电荷条纹相的微观图景：基于Emery模型的DMRG与DQMC深度解析

0. 执行摘要

探索铜氧化物（Cuprates）高温超导体中的交织序（Intertwined Orders）是凝聚态物理学数十年来的核心物理前沿之一。长期以来，学术界普遍采用简化的单带Hubbard模型来近似描述高温超导体的低能物理行为。然而，单带模型是否能够完全等效地捕获包含氧轨道自由度的多带物理本质（如张-Rice单态的局域化及其退禁闭机制），一直存在巨大的理论争议。

本文针对最新发表的关于三带Emery模型中涨落交织电荷条纹相的理论工作进行深度学术剖析。该研究联合应用了最前沿的数值关联方法——基态密度矩阵重整化群（DMRG）和有限温行列式量子蒙特卡洛（DQMC），首次在包含显式铜 $d$ 轨道和氧 $p$ 轨道的Emery模型中，系统性地构建出了强涨落电荷条纹相的微观图景。研究证实，在广泛使用的物理参数区间内，氧轨道上存在着鲁棒的、呈现强轨道选择性的电荷条纹，这些条纹与铜轨道上的自旋畴壁（Spin Domain Walls）高度锁定。更重要的是，在极低温DQMC计算下，研究揭示了系统从高温下的 $B_{2g}$ 向列向易感性（Nematic Susceptibility）向低温下 $B_{1g}$ 键电荷向列易感性占优的交叉行为，且这一向列涨落的增强与实空间单向电荷条纹的形成以及动能各向异性的非单调掺杂依赖性表现出强烈的微观耦合。本解析将详尽拆解其理论基础、算法细节、基准数据，并评述该工作的深远物理意义与技术局限性。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题：超越单带Hubbard模型的必要性

强关联电子体系（Strongly Correlated Electron Systems）的本质难题在于，极小的能量尺度差异决定了截然不同的量子基态。在铜氧化物高温超导体中，超导相（Superconductivity）与反铁磁自旋密度波（SDW）、电荷密度波（CDW）及向列相（Nematicity）在相图上交织共存，其基态能量差异往往仅在每个格点毫电子伏特（meV）量级（相当于跃迁能 $t$ 的千分之一）。

传统的单带Hubbard模型通过张-Rice单态（Zhang-Rice Singlet, ZRS）近似，将氧 $p$ 轨道上的掺杂空穴与铜 $d$ 轨道上的局域空穴束缚在一起，作为单一的有效带电载流子。然而，真实的铜氧化物是典型的电荷转移型绝缘体（Charge-Transfer Insulators）。实验表征（如STM、RIXS）表明，电荷不均匀调制和向列性对称破缺在氧轨道上有着极其显著的表现。这就带来了一个根本性的科学问题：在显式包含氧原子的三带Emery模型中，电荷条纹是如何在空间和特定轨道上显现的？它与自旋条纹的锁定关系是否与单带模型一致？向列涨落的微观动力学起源是什么？

1.2 理论基础：三带Emery模型及其关联算符

三带Emery模型构建在二维 Lieb 晶格（Lieb Lattice）之上，该晶格由铜格点（位于胞中心）和氧格点（位于 $x$ 和 $y$ 桥部）组成，分别承载 Cu $d_{x^2-y^2}$ 轨道、O $p_x$ 轨道和 O $p_y$ 轨道。在空穴语言下，其系统哈密顿量可表示为：

$$H = K + \Delta_{pd}\sum_{i,\sigma} n_{i,\sigma} + \sum_i U_i n_{i,\uparrow} n_{i,\downarrow} - \mu \sum_{i,\sigma} n_{i,\sigma}$$

其中，第一项 $K$ 为动能项（跃迁项），包含铜-氧之间的最近邻跃迁 $t_{pd}$ 以及氧-氧之间的次近邻直接跃迁 $t_{pp}$：

$$K = -\sum_i t_{pd} \left( \xi^a_i + \xi^b_i + \xi^c_i + \xi^d_i \right) - \sum_i t_{pp} \left( \xi^e_i + \xi^f_i + \xi^g_i + \xi^h_i \right)$$

跃迁算符 $\xi^m_i$ 的具体形式与轨道位相约定（Sign Convention）密切相关。为了保证物理波函数的正确对称性，定义如下：

$$\xi^a_i = \sum_\sigma \left( p^{x\dagger}_{i,\sigma} d_{i,\sigma} + \text{h.c.} \right)$$

$$\xi^b_i = -\sum_\sigma \left( p^{y\dagger}_{i,\sigma} d_{i,\sigma} + \text{h.c.} \right)$$

$$\xi^c_i = -\sum_\sigma \left( p^{x\dagger}_{i-\hat{x},\sigma} d_{i,\sigma} + \text{h.c.} \right)$$

$$\xi^d_i = \sum_\sigma \left( p^{y\dagger}_{i-\hat{y},\sigma} d_{i,\sigma} + \text{h.c.} \right)$$

$$\xi^e_i = -\sum_\sigma \left( p^{x\dagger}_{i,\sigma} p^y_{i+\hat{x}-\hat{y},\sigma} + \text{h.c.} \right)$$

$$\xi^f_i = \sum_\sigma \left( p^{x\dagger}_{i,\sigma} p^y_{i-\hat{y},\sigma} + \text{h.c.} \right)$$

$$\xi^g_i = -\sum_\sigma \left( p^{x\dagger}_{i,\sigma} p^y_{i,\sigma} + \text{h.c.} \right)$$

$$\xi^h_i = \sum_\sigma \left( p^{x\dagger}_{i,\sigma} p^y_{i+\hat{x},\sigma} + \text{h.c.} \right)$$

在哈密顿量中，$\Delta_{pd} = \epsilon_p - \epsilon_d$ 为电荷转移能，$U_i$ 为格点排斥，其中铜格点上 $U_d = U_{dd}$，氧格点上 $U_p = U_{pp}$。为了探究向列性与条纹相的本质，研究引入了静态电荷易感率与自旋易感率算符：

$$\chi_{c,ij}(\omega = 0) = \int_0^\beta \left( \langle n_i(\tau) n_j(0) \rangle - \langle n_i(\tau) \rangle \langle n_j(0) \rangle \right) d\tau$$

$$\chi_{s,ij}(\omega = 0) = \int_0^\beta \langle S^z_i(\tau) S^z_j(0) \rangle d\tau$$

针对格点体系在 $C_4$ 对称性下的自发破缺，特别定义了键电荷向列算符。在 $B_{1g}$ 键电荷向列对称性下，其物理响应由静态 $B_{1g}$ 键电荷向列易感率（$\chi_{B_{1g}}$）表征，对应算符形式为：

$$\rho_{B_{1g}} = \sum_i \left( \xi^a_i - \xi^b_i + \xi^c_i - \xi^d_i \right)$$

相应地，对于 $B_{2g}$ 对称性（即沿对角线方向的向列畸变），其算符定义为：

$$\rho_{B_{2g}} = \xi^e_i - \xi^f_i + \xi^g_i - \xi^h_i$$

1.3 技术难点：多带体系的指数爆炸与负号问题

多带Emery模型的数值求解面临双重极其严苛的物理瓶颈：

  1. 费米子负号问题（Fermion Sign Problem）： 在辅助场量子蒙特卡洛（DQMC）模拟中，由于非半满掺杂以及强关联导致的费米子配分函数项出现负权重，平均负号 $\langle S \rangle$ 随系统尺寸 $N$、交互强度 $U$ 和反温 $\beta = 1/T$ 呈指数衰减。这限制了 DQMC 在真实物理参数下的低温探查能力。
  2. 希尔伯特空间的轨道维数灾难： 在 DMRG 计算中，由于每个 Lie b 原胞包含三个物理格点，其局部希尔伯特空间维度大幅增加。要在宽尺寸的二维圆柱格点上保持极高的能量精度并精确捕获微弱的电荷条纹振幅，要求极高且收敛极其困难的矩阵键维度（Bond Dimension）。

1.4 方法细节：互补数值技术的联合出击

为了克服上述物理瓶颈，本工作设计了 DMRG 与 DQMC 互补的联合架构：

  • 基态DMRG（密度矩阵重整化群）： 在 $T=0$ 基态下，采用具有严格 $SU(2)_{\text{spin}} \otimes U(1)_{\text{charge}}$ 对称性的多轨道 DMRG 算法。在 $24 \times 4$ 和 $32 \times 4$ 圆柱体几何下（边界端接有 $d$ 和 $p_y$ 轨道，无悬空的 $p_x$ 轨道），将矩阵乘态（MPS）的键维度 $D$ 提升至前所未有的 $16,000$（等效于无非阿贝尔对称性下的 $m \approx 47,000$ 态），从而在控制截断误差 $\epsilon \le 3 \times 10^{-7}$ 的高精度下，获取氧轨道电荷分布的基态定量基准。
  • 有限温DQMC（行列式量子蒙特卡洛）： 在有限温度下，在具有完全周期性边界条件（PBC）的 $16 \times 4$ 和 $8 \times 8$ 格点上，执行精确且无偏差的蒙特卡洛抽样。通过 Jackknife 重抽样方法解决误差相关性，Imaginary Time 间隔设为 $\delta\tau \le 0.05 \text{ eV}^{-1}$（对于 $U_{dd}=8.5$ eV）或 $\delta\tau \le 0.08 \text{ eV}^{-1}$（对于 $U_{dd}=6$ eV）。

2. 关键 Benchmark 体系、计算所得数据与性能分析

2.1 Benchmark 参数集设计

研究精细设计了三套具有不同物理侧重的 Emery 模型参数体系（单位均为 eV）：

参数集标识$t_{pd}$$t_{pp}$$\Delta_{pd}$$U_{dd}$$U_{pp}$物理侧重及数值计算特点
参数集 A (Standard)$1.13$$0.49$$3.24$$8.5$$4.1$真实铜氧化物超导体的物理模拟。费米负号最严重，最低可达温度受限（$\beta_{max} = 6.5$）。
参数集 B (Standard, $U_{pp}=0$)$1.13$$0.49$$3.24$$8.5$$0.0$探究氧格点库仑排斥对电荷条纹稳定性的直接贡献。负号问题略有缓解（$\beta_{max} = 10.0$）。
参数集 C (Large $J$)$1.13$$0.49$$3.0$$6.0$$0.0$有效反铁磁超交换作用 $J$ 较大，且负号问题最轻微。允许 DQMC 达到极低温（$\beta_{max} = 12.0$, $T/J \approx 0.29$）。

2.2 基态 DMRG 的精准定量基准

在 $32 \times 4$（空穴掺杂 $\delta = 1/8$）以及 $24 \times 4$（$\delta = 1/12, 1/6$）圆柱体上，DMRG 清晰地给出了空间电荷密度分布与自旋关联的基态解（如图1所示）：

  • 自旋与电荷周期的锁定关系： 当掺杂浓度为 $\delta = 1/8$ 时，反铁磁自旋条纹的周期为 $\lambda_s = 8$。电荷密度涨落的周期则完美锁定在 $\lambda_c = 4$（即 $\lambda_c = 1/(2\delta)$，对应半满畴壁条纹）。自旋条纹在通过电荷河流时发生 $\pi$ 位相移动。
  • 极强的轨道选择性（Orbital-Selective Modulations）： DMRG 数据表明，电荷密度调制的幅值呈现出强烈的各向异性： $$\text{Amplitude}(p_y) \gg \text{Amplitude}(p_x) \gg \text{Amplitude}(d)$$ 具体而言，在纵向条纹构型下，指向“电荷河流”方向（即垂直方向 $y$）的氧 $p_y$ 轨道上的电荷调制振幅最大。这微观地阐明了，电荷积累并非均匀分布在铜-氧面上，而是以特定的 $p$-对称性局部集中在畴壁之上。

2.3 DQMC 有限温静态易感率分析

利用 DQMC 静态电荷易感率的傅里叶变换 $\chi^{\mu\nu}_c(\mathbf{q}, \omega=0)$，可以在动量空间对电荷条纹的短程序进行精确定量。拟合采用周期性双洛伦兹函数（Periodic Double Lorentzian）：

$$\chi(\mathbf{q}) = \frac{A}{\left[(\mathbf{q} - \mathbf{q}^*)^2 + \xi^{-2}\right]} + \frac{A}{\left[(\mathbf{q} + \mathbf{q}^*)^2 + \xi^{-2}\right]} + \text{BG}$$

研究得到了极具物理冲击力的发现（如图3所示）：

  • Yamada 关系的严格成立： 随着掺杂浓度 $p$ 从 $0.05$ 逐渐增加到 $0.25$，提取的自旋调制波矢 $\delta_{SDW}$ 与一半的电荷调制波矢 $\delta_{CDW}/2$ 表现出高度线性相关的轨迹，完美符合经典的 $\delta_{CDW} = 2p$ 以及 $\delta_{SDW} = p$ 的偏离半满条纹相 Yamada 关联（如图3(c)）。这表明，哪怕在温度高于基态的强涨落区间，Emery 模型中多轨道的电荷-自旋交织关联仍然是动力学上最主导的低能物理激发。
  • 局域各向异性响应： 在 $16 \times 4$ 格点上，由于几何对称性破缺，$\chi^{xx}_c$（$x$ 方向氧格点关联）和 $\chi^{yy}_c$（$y$ 方向氧格点关联）表现出不同的物理行为。静态电荷关联显示出，在氧格点上的局部化学势扰动倾向于在与其平行的轨道方向上诱导出强烈的、各向异性的单向电荷调制（如图2）。

2.4 向列易感性与温度演化的交叉行为

在具有完整 $C_4$ 旋转对称性的 $8 \times 8$ 格点上，$\chi_{B_{1g}}$ 和 $\chi_{B_{2g}}$ 的温度依赖性（图5）揭示了系统非平庸的物理图像：

  • 向列对称破缺的主导通道切换： 在高温区（$T > 0.4$ eV），$\chi_{B_{2g}}$ 显著大于 $\chi_{B_{1g}}$。这主要来源于铜-氧电荷转移能所诱导的局域高频电荷波动。然而，随着温度下降，在 $T \approx 0.2$ eV（对应自旋和电荷条纹相开始强烈涨落的特征温度），$\chi_{B_{1g}}$ 呈现出指数级的急剧发散增长；而 $\chi_{B_{2g}}$ 则发生断崖式下跌。这微观地指明，系统的自发对称破缺倾向从高温的 $B_{2g}$ 通道向低温的 $B_{1g}$ 键电荷向列通道进行交叉（Crossover），且这一交叉行为直接由一维单向条纹的物理成核所驱动。

3. 代码实现细节、复现指南及开源软件库

3.1 代码实现架构与算法核心

  • DMRG部分： 基于张量网络态（Matrix Product States）框架，核心逻辑在于将二维 Lieb 晶格在一维一维扫描路径（Snake Path）上进行紧凑的一维化映射。对于 $L_x \times L_y$ 胞体系，每一个原胞对应 3 个 MPS 物理格点。为保证大矩阵键维度下计算的进行，哈密顿量矩阵元显式实现了 $SU(2)_{spin} \otimes U(1)_{charge}$ 守恒。非阿贝尔 $SU(2)$ 守恒极大缩减了奇异值分解（SVD）时对局域哈密顿子空间的维数要求。
  • DQMC部分： 采用格林函数（Green’s Function）蒙特卡洛抽样。通过引入 Hubbard-Stratonovich 变换，将强关联排斥项 $U_i$ 转换为辅助磁场，通过行列式更新求得配分函数。静态关联函数通过沿虚时间方向的数值积分 $\int_0^\path d\tau$ 实现精确测量，为压制误差，计算中利用了晶格的平移对称性。

3.2 复现指南与步骤

为了复现本工作的有限温 DQMC 模拟数据，可基于 Edwin Huang 开发并公开维护的开源高性能 DQMC 软件包。其详细步骤如下：

第一步：克隆开源代码库并进行本地编译

该研究的有限温 DQMC 测量引擎完全基于 C++/Fortran 架构，可从开源 Repo 链接获取：

git clone https://github.com/edwnh/dqmc.git
cd dqmc
make all

注：编译需要本地已配置高性能数学库（如 Intel MKL、LAPACK、BLAS）以及具备 MPI 编译环境。

第二步：配置物理参数输入文件（Input File）

在配置文件中，定义三带 Emery 模型的跃迁参数、排斥力以及 Lieb 格点几何结构。以下为复现参数集 C（$U_{dd}=6$ eV，$\beta = 12$ eV$^{-1}$，格点尺寸 $8 \times 8$）的示例输入控制段落：

# 格点几何定义
LATTICE = Lieb
Nx = 8
Ny = 8

# 跃迁矩阵元定义（依据文中位相约定）
tpd = 1.13
tpp = 0.49
delta_pd = 3.0

# 局域相互作用设置
U_d = 6.00
U_p = 0.00

# 虚时切片设置
beta = 12.0
L = 240         # 虚时步数，对应 dtau = 0.05

# 蒙特卡洛采样控制
warmup = 10000
sweeps = 37500
num_chains = 100  # 使用MPI多链并行运行，总链数建议不少于100

第三步：启动计算与数据收集

在计算集群上，使用 MPI 启动大规模并行采样计算：

mpirun -np 100 ./dqmc_exec -i config_input.txt -o results_output.dat

第四步：执行 Jackknife 误差分析与静态物理量提取

计算完毕后，运行代码仓中附带的 Python 数据分析脚本，对虚时间关联函数进行格林函数傅里叶变换、实空间积分以及误差传播评估，最终绘制出如 $\chi_{B_{1g}}(T)$ 与 $\chi_{B_{2g}}(T)$ 等特征演化曲线。

4. 关键引用文献与局限性批判

4.1 关键引用文献

  1. Emery (1987) [3]: 理论物理基石。首次提出三带模型来描述高温超导体的 $CuO_2$ 平面，指出掺杂空穴主要定域在氧 $p$ 轨道，明确了氧自由度在超导配对及关联效应中的核心作用。
  2. Zhang & Rice (1988) [4]: 著名的 Zhang-Rice 单态假说。认为铜氧化物低能激发可映射为单带 Hubbard 模型。本工作通过展示氧轨道上强烈的条纹不均匀性与轨道选择响应，在更高能级上对该映射的严格性提出了挑战。
  3. Huang et al. (Science 2017) [9]: 首次使用数值 DQMC 技术证实了 Emery 模型中自旋条纹的存在。本工作是其在电荷条纹、氧格点选择性及向列易感性领域的深层物理拓展。
  4. Mai et al. (PNAS 2024) [22]: 探讨了多带模型中弱电荷密度波的可能特征。本工作直接对其进行了精确的基态 DMRG 标定，解决了多带参数集敏感性这一科学疑虑。
  5. Huang et al. (PRB 2023) [23]: 系统化探究了单带 Hubbard 模型中的涨落交织条纹。为本工作提供了直接的物理对照标准。

4.2 局限性深度评论

尽管本工作在多带关联数值模拟领域取得了突破，但从严苛的量子物理研究视角来看，仍存在以下技术局限性与理论不确定性：

  1. 极度严苛的低温极限费米负号瓶颈： 对于最贴近真实铜氧化物超导材料的参数集 A（$U_{dd}=8.5$ eV, $U_{pp}=4.1$ eV），在掺杂区域，DQMC 模拟所面临的负号问题极其严重，导致实际计算无法超越反温 $\beta = 6.5$ eV$^{-1}$（物理温度约 $1800$ 开尔文）。尽管本工作通过调小 $U_{dd}$ 至 $6$ eV（参数集 C）将计算推进到 $\beta = 12$ eV$^{-1}$，但这并非真实铜氧化物的物理区间。在最受关注的高温超导配对区（几百开尔文至零温），电荷条纹与超导态的直接微观竞争仍然超出了目前 DQMC 的直接处理能力。
  2. 准一维圆柱几何的强人工条纹偏置效应： DMRG 采用的 $L_x \times 4$（宽度 $W_y = 4$）圆柱格点，由于其准一维极限制约和边界截断，天然对单向电荷条纹（Unidirectional Stripes）具有极强的束缚和强化作用。在真正的二维热力学极限（Thermodynamic Limit, $2D$ 无边界偏置）下，这种单向电荷条纹是否会被强热涨落或量子涨落熔化，退化为双向棋盘状（Checkerboard）电荷密度波或是无序的液晶相，DMRG 无法给出毫无争议的最终定论。
  3. 长程库仑相互作用的完全缺失： Emery 模型仅包含局域的 Hubbard 斥力 $U_{dd}$ 和 $U_{pp}$。然而在真实材料中，由于不完美的屏蔽效应，铜-氧、氧-氧之间存在非零的长程排斥（如次近邻长程库仑作用 $V$）。长程排斥会带来极其强烈的静电能惩罚，倾向于剧烈摧毁任何电荷聚集（即电荷条纹的去稳定化）。忽视长程库仑作用可能导致模型系统性地高估了电荷条纹在相图中的鲁棒性。

5. 补充物理机制：微观微扰图像与实验对照

5.1 强关联极限下的条纹形成机制（Perturbative Picture）

为了直观理解氧 $p$ 轨道上电荷条纹形成的微观机理，我们可以诉诸强关联极限（$U_{dd} \to \infty$）下的有效微扰图像：

在未掺杂（半满）基态下，每个铜格点由于强关联局域化承载一个空穴，通过超交换作用 $J = 4t_{pd}^4 / \Delta_{pd}^3$ 形成长程反铁磁背景。当我们在氧轨道上注入一个额外空穴时，该空穴由于强反铁磁交换作用 $J_{pd} \sim t_{pd}^2 / \Delta_{pd}$，倾向于与相邻的两个铜格点自旋组合成张-Rice单态。

从能量优化角度来看，在反铁磁背景中，空穴的运动会受到自旋阻挫的惩罚（即空穴跃迁会破坏其路径上的反铁磁自旋对齐，产生高能量惩罚）。为了减少这一自旋磁能惩罚，多个注入的空穴会自发地聚集在空间的一维狭长区域上。这一区域充当了自旋畴壁的角色（畴壁两侧的反铁磁畴位相相差 $\pi$）。

在该畴壁上，自旋秩序被完全熔化，形成自旋液体般的高度起伏区。因为此处的反铁磁秩序已经破缺，空穴跃迁不再受自旋阻挫惩罚，从而大幅释放了动能。这就构成了“电荷河流”（Rivers of Charge）。因为纵向畴壁分布在 $y$ 方向上，为了最大化空穴在畴壁内部一维流动的跃迁概率，空穴主要倾向于占据指向条纹方向的氧 $p_y$ 轨道（如图4(e, f)所示的实空间局部积聚），这就是氧轨道选择性电荷条纹的微观起源。

5.2 前沿光谱实验的交相呼应

本工作所得数值图像在近年多项高温超导前沿实验中找到了清晰的物理对应：

  • STM（扫描隧道显微镜）： STM 在 BSCCO（铋系铜氧化物）中观测到的非均匀局域局域态密度波动，其空间分布呈现出强烈的内胞向列性（Intra-unit cell Nematicity），且该向列自发破缺主要定域在氧原子轨道上，与本工作 DMRG 揭示的氧轨道振幅极大性完美契合。
  • RIXS（共振非弹性X射线散射）： 在 YBCO 等多种铜氧化物中，RIXS 探测到了在特定温度和能量窗口下极其活跃的电荷密度波涨落。这一涨落即使在超导转变温度以上依然以短程强起伏形式存在，本工作 DQMC 拟合的静态电荷关联函数及其与 Yamada 关系的相符性，为解释 RIXS 的有限温电荷涨落提供了第一手、微观非微扰的计算依据。
  • 向列阻抗实验： 超导转变附近的输运测量探测到了剧烈的 $B_{1g}$ 对称性电阻率各向异性（电阻向列性）。这与本工作发现的随着温度下降，B1g 键电荷向列易感性陡然发散的特征高度一致。

5.3 展望：未来计算方法演进的方向

要彻底解决铜氧化物相图这一世纪难题，未来理论计算方法需在以下维度取得突破：

  1. 张量网络态的2D扩展（2D Tensor Networks）： 引入不依赖边界的无限PEPS（iPEPS）方法，直接在 $2D$ 热力学极限下搜索三带 Emery 模型的基态，确认单向条纹相在真实 2D 环境下的稳定边界。
  2. 机器学习辅助的量子蒙特卡洛（ML-assisted QMC）： 引入人工神经网络对配分函数进行正则化映射，大幅缓解甚至消除强关联下的费米子负号问题，从而使得多带模型的有限温 DQMC 模拟能够真正下探至极低温物理区间（$T \le 100$ K），直接揭示超导态与条纹相的竞争相图。
  3. 非局域关联动力学方法（DCA / CDMFT）： 将动态平均场理论（DMFT）拓展至多带团簇体系，直接计算多轨道空间分辨的自能函数与动态谱函数，为实验光谱学提供全方位的理论支撑。