来源论文: https://arxiv.org/abs/2607.07067v1 生成时间: Jul 09, 2026 18:45

从第一性原理热力学统计建模到强关联金属-绝缘体转变：SrTi1-xVxO3 固溶体的 GQCA+DFT+DMFT 深度解析

0. 执行摘要

在过渡金属氧化物（$\text{ABO}_3$ 钙钛矿）研究中，控制电子关联强度与化学无序的竞争是凝聚态物理的核心前沿之一。固溶体 $\text{SrTi}_{1-x}\text{V}_x\text{O}_3$ 完美地架起了从带状绝缘体 $\text{SrTiO}_3$（$3d^0$ 电子排布）到中等关联 Mott 金属 $\text{SrVO}_3$（$3d^1$ 电子排布）之间的桥梁。然而，由于该体系同时存在着强电子-电子关联（Mott 局域化）与局域化学无序（Anderson 局域化），传统的第一性原理方法在处理该体系时面临着巨大的理论瓶颈。单独的超胞（Supercell）计算无法捕捉热力学无序的固溶体特征，而大规模特殊准随机结构（SQS）与动力学平均场理论（DMFT）的结合在计算成本上又是无法承受的。

为了克服这一瓶颈，奥地利格拉茨工业大学（TU Graz）的 Luka Wibmer 及其合作者提出了一种全新的无参数、自洽热力学统计建模框架。该框架将**广义准化学近似（Generalized Quasichemical Approximation, GQCA）与第一性原理密度泛函理论（DFT）及动力学平均场理论（DMFT）**深度耦合，实现了对 $\text{SrTi}_{1-x}\text{V}_x\text{O}_3$ 合金全组分范围的热力学混溶间隙（Miscibility Gap）和金属-绝缘体转变（Metal-Insulator Transition, MIT）的统一描述。

本研究的核心成果包括：

  1. 热力学相图预测：首次在无经验参数修正的前提下，通过统计配分函数预测出该固溶体在高温下存在一个临界温度 $T_c = 1443\text{ K}$ 的混溶间隙。该预测与淬火实验中固溶体在 $1100^\circ\text{C}$ 下发生相分离的实验事实高度吻合，大幅修正了以往基于 $\text{GGA}+U$ 团簇膨胀预测的 $\sim 500\text{ K}$ 偏低结果。
  2. 金属-绝缘体相变重现：单纯的 DFT 计算错误地预测该体系在所有 $x > 0$ 组分下均为金属，而耦合了局域关联效应的 $\text{DFT}+\text{DMFT}$ 团簇系综平均方法成功重现了实验观测到的 MIT，其在 $x \approx 0.3$ 以下呈现绝缘态，而在 $x \to 1$ 时演化为相干金属态。
  3. 渗流理论微观图像：将不同团簇分类为金属或绝缘体，并在简单立方（SC）格子上进行位置渗流（Site Percolation）模拟，揭示出在 $x \approx 0.4$ 附近存在一个尖锐的系统尺度导电通路建立过程，完美解释了实验测得的电导率突变。

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

1.1 核心科学问题：Mott-Anderson 局域化的共存与竞争

在 $\text{SrTi}_{1-x}\text{V}_x\text{O}_3$ 固溶体中，MIT 的物理机制长期存在 Mott 型局域化（电子库仑排斥导致能隙打开）与 Anderson 型局域化（无序势场导致波函数局域化）的争议。实验上，由于制备工艺、薄膜厚度和衬底约束的差异，测得的临界组分 $x_c$ 在 $0.4 \sim 0.7$ 之间宽幅分布。在微观尺度，不同阳离子（$\text{Ti}^{4+}$ 与 $\text{V}^{4+}$）的局域空间排布不仅直接决定了静电势无序度，更通过调控 $\text{V}-t_{2g}$ 轨道的局域重叠积分（即带宽 $W$）强力支配着局部关联参数 $U/W$ 的大小。因此，如何在一个统一的、热力学自洽的理论框架下，同时且对等地处理局域化学无序和强电子关联，是该领域亟待解决的瓶颈问题。

1.2 理论基础：广义准化学近似 (GQCA)

传统合金建模常用“团簇膨胀”（Cluster Expansion）方法，但其对于强关联电子结构的耦合极其困难。GQCA 提供了一种替代路径：它将无序合金视为由一系列在统计和能量上相互独立的团簇（Clusters）组成的系综。对于拟二元合金 $\text{A}_{1-x}\text{B}_x\text{C}$，体系被划分为 $J$ 个对称非等价的团簇类别，每个团簇 $j$ 具有特定的总能量 $E_j$、简并度 $g_j$、可替代位点数 $n$ 以及包含 $\text{B}$ 原子个数 $n_j$。在该系综中，某特定团簇 $j$ 的出现概率（即发生概率）记为 $p_j$，定义为：

$$p_j = \frac{M_j}{M}$$

其中 $M_j$ 为第 $j$ 类团簇的数量，$M$ 为总团簇数。

GQCA 的核心在于通过最小化吉布斯混合自由能 $\Delta G$ 来确定这些概率 $p_j(x, T)$：

$$\Delta G = \Delta H - T\Delta S$$

其中，混合焓 $\Delta H(x, T)$ 表达为各团簇还原超额焓（Reduced Excess Enthalpies）$\Delta_j$ 的加权和：

$$\Delta H(x, T) = M \sum_{j=1}^J p_j \Delta_j$$$$\Delta_j = H_j - \frac{n - n_j}{n} H_A - \frac{n_j}{n} H_B$$

这里 $H_j$ 是第 $j$ 个团簇通过第一性原理结构弛豫后得到的总能量，$H_A$ 和 $H_B$ 分别代表端点组分 $\text{SrTiO}_3$ 和 $\text{SrVO}_3$ 单胞能量。配置熵 $\Delta S(x, T)$ 利用玻尔兹曼熵公式在微正则系综下导出：

$$\Delta S(x, T) = k_B \ln W$$

合金配置数 $W$ 可以通过将 $N$ 个位置上放置 $N_A$ 个 $\text{A}$ 原子和 $N_B$ 个 $\text{B}$ 原子的排布数，乘以随机配置下产生相应团簇系综的多元概率来近似表达：

$$W = \frac{N!}{N_A! N_B!} \frac{M!}{\prod_{j=1}^J M_j!} \prod_{j=1}^J \left(p_j^0\right)^{M_j}$$

其中 $p_j^0$ 表示在完全随机混排下，第 $j$ 类团簇出现的几率：

$$p_j^0 = g_j x^{n_j} (1-x)^{n - n_j}$$

应用斯特林公式（$\ln \alpha! \approx \alpha \ln \alpha - \alpha$），混合自由能 $\Delta G$ 最终写为：

$$\Delta G(x, T) = M \sum_{j=1}^J p_j \Delta_j + N k_B T \left( x \ln x + (1-x) \ln (1-x) \right) + M k_B T \sum_{j=1}^J p_j \ln \left(\frac{p_j}{p_j^0}\right)$$

在约束条件 $\sum_{j=1}^J p_j = 1$ 和 $\sum_{j=1}^J n_j p_j = n x$ 下，利用拉格朗日乘子法求自由能极值，可导出平衡状态下的团簇概率 $p_j$：

$$p_j = \frac{g_j \eta^{n_j} e^{-\Delta_j \beta}}{\sum_{i=1}^J g_i \eta^{n_i} e^{-\Delta_i \beta}}$$

其中 $\beta = (k_B T)^{-1}$，而参数 $\eta$ 由拉格朗日乘子 $\lambda_L$ 决定：$\eta = \frac{x e^{\lambda_L \beta}}{1-x}$。通过求解以下 $n$ 阶多项式方程的唯一实数根，可以唯一确定 $\eta$：

$$\sum_{j=1}^J (nx - n_j) g_j e^{-\Delta_j \beta} \eta^{n_j} = 0$$

任何组分与温度依赖的材料宏观性质 $\mathcal{P}(x, T)$ 均可通过对各独立团簇的相应性质 $\mathcal{P}_j$ 进行系综平均得到：

$$\mathcal{P}(x, T) = \sum_{j=1}^J p_j(x, T) \mathcal{P}_j$$

其统计标准差（表征无序引起的局域涨落）计算公式为：

$$\Delta \mathcal{P}(x, T) = \sqrt{\sum_{j=1}^J p_j(x, T) \mathcal{P}_j^2 - \left( \sum_{j=1}^J p_j(x, T) \mathcal{P}_j \right)^2}$$

1.3 技术难点与应对策略

  1. 多维配置空间的离散化：要在全组分尺度上捕捉局部无序，超胞的选择必须在“能充分描述局部配位环境多样性”与“计算成本”之间妥协。作者采用了 $2 \times 2 \times 2$ 钙钛矿超胞（包含 8 个 $\text{B}$ 位阳离子位点，公式单元数为 8）。利用对称性群论分析，将 $2^8 = 256$ 种可能的占据态压缩为 22 个对称非等价的团簇类别，极大降低了计算成本。
  2. 强关联体系的无经验参数能谱描述：常规的 $\text{DFT}+U$ 存在强烈的经验依赖。为了解决这一问题，作者在每个团簇弛豫后，提取其 $\text{V}-t_{2g}$ 轨道作为关联空间，构建低能有效 Kanamori 汉密尔顿量，采用单基点（Single-site）连续时间量子蒙特卡洛（CT-QMC）杂化膨胀算法求解 $\text{DFT}+\text{DMFT}$ 方程。此方法在完全无偏的能谱自能修正下，保证了电子结构的客观性和准确性。

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

为了验证 GQCA+DFT+DMFT 联合框架的可靠性，作者对 $\text{SrTi}_{1-x}\text{V}_x\text{O}_3$ 体系的结构性质、热力学稳定性和能谱特征进行了全方位的 Benchmark 评估。

2.1 晶格参数与 Vegard 定律的偏离分析

在 $T = 1800\text{ K}$（该温度与实验上采用 $\text{SrCO}_3$、$\text{TiO}_2$ 和 $\text{V}_2\text{O}_5$ 进行固相反应合成固溶体的温度一致）下，通过 GQCA 系综平均计算得到的晶格常数 $a(x)$ 如图 1 所示。

  • 与实验数据的对比：实验测得端点组分 $\text{SrTiO}_3$ 的晶格常数为 $3.90\text{ \AA}$，$\text{SrVO}_3$ 为 $3.84\text{ \AA}$。计算弛豫结果与实验值的相对误差低于 1%。
  • 线性关系（Vegard's Law）：系综平均的晶格常数 $a(x)$ 随钒组分 $x$ 呈现出几乎完美的线性下降趋势。这是由于 $\text{V}^{4+}$ 的离子半径（约 $0.58\text{ \AA}$）小于 $\text{Ti}^{4+}$（约 $0.605\text{ \AA}$）。
  • 实验中细微弯曲（Bowing）的物理阐释：实验在富钛端（$x$ 较小）观察到晶格常数偏离线性的微小向下弯曲。作者指出，理论模型未能在富钛端完全重现这一细节，原因在于：
    1. 结构弛豫是在 DFT 级别完成的，未考虑 DMFT 自能修正对局域晶格受力的自洽反馈；
    2. 实际样品中由于退火和氧分压控制，可能存在复杂的多价态（如 $\text{V}^{3+}$/$\text{V}^{5+}$）共存以及氧空位缺陷，改变了局域 B 位离子的有效半径。

2.2 混合热力学与混溶间隙 (Miscibility Gap)

这是本项工作最为惊艳的成果。图 2 详细拆解了该体系在 $400\text{ K} \sim 2000\text{ K}$ 温度区间内的各项热力学本征量：

物理量变化特征与物理意义
团簇超额焓 $\Delta_j$22 个非等价配置的 $\Delta_j$ 全部为正值（分布在 $5 \sim 15\text{ meV/atom}$，见图 2(a)）。这表明任何局域混排状态在焓上都是不利的，体系本征具有极强的自发相分离（Phase Separation）趋势。
混合焓 $\Delta H$在全组分范围内恒为正值，在 $x=0.5$ 附近达到最大值 $\sim 11.5\text{ meV/atom}$（见图 2(b)）。
混合自由能 $\Delta G$在 $T < 1443\text{ K}$ 时，$\Delta G(x)$ 曲线在中间组分出现下凹（凹陷区），即双轴点（Binodal Points）之间可做公切线。这意味着在此温度以下，均匀固溶体在热力学上是不稳定的（见图 2(d)）。
  • 临界温度预测：计算预测出混溶间隙的最高临界温度 $T_c = 1443\text{ K}$（约 $1170^\circ\text{C}$，见图 2(e)）。
  • 与实验的一致性验证：Patino 等人的淬火实验报道指出，在 $1100^\circ\text{C}$（$1373\text{ K}$，低于理论 $T_c$）长期退火后，$\text{SrTi}_{1-x}\text{V}_x\text{O}_3$ 固溶体彻底相分离为纯的 $\text{SrTiO}_3$ 和 $\text{SrVO}_3$。而通过快速淬火可以将其保留在亚稳的均匀单相。这一实验现象与本工作预测的混溶间隙边界高度吻合。
  • 与以往理论模型的巨大修正：Liu 等人此前使用 $\text{GGA}+U$ 结合传统 Monte Carlo 模拟预测的临界温度仅为 $\sim 500\text{ K}$，严重低估了相分离驱动力。这证明了 GQCA 框架在无需人工调控参数的前提下，展现出了极高的定量预测精度。

2.3 电子能谱特征与金属-绝缘体转变 (MIT)

作者通过对比单纯 DFT 密度泛函计算与 $\text{DFT}+\text{DMFT}$ 能谱结果，揭示了关联效应在 MIT 中的支配地位（见图 3、图 4）：

  1. 单纯 DFT 的失效：如图 3(a) 所示，对于所有 $x > 0$ 的组分，DFT 预测的费米能级态密度 $N(\epsilon_F)$ 均大于零，表现为金属。这与实验中观察到的大范围电绝缘性完全违背。这是因为传统的局域密度近似（LDA/GGA）高估了电子去局域化能，无法描述 Mott 绝缘体的能隙打开。
  2. $\text{DFT}+\text{DMFT}$ 零频谱权重 $A(\omega=0)$：引入局域库仑排斥（$U=4.0\text{ eV}$）后，当 $x < 0.3$ 时，所有团簇的费米能级谱权重 $A(\omega=0)$ 严格归零，表明带隙完全打开，体系进入 Mott 绝缘态。当 $x \sim 0.5$ 时，系综中绝缘团簇与金属团簇并存，表现出剧烈的局域配置涨落（标准差达到最大值）。在 $x \to 1$ 时，谱权重平滑过渡到强关联金属 $\text{SrVO}_3$ 的特征值（见图 3(b)）。
  3. 谱函数演化（MaxEnt 动力学自能解析延拓）：通过最大熵方法，作者重构了实频谱函数 $A(\omega)$（见图 4）。富钛端（如 $x=0.1$）能谱在费米能级处表现为明晰的 Mott 能隙。随着 $x$ 增加，V 浓度提升，局域 V 团簇相互连通，导致 $t_{2g}$ 有效带宽 $W$ 增大。当其克服库仑排斥时，能隙在 $x \approx 0.35$ 处闭合，并在费米能级处激发出相干准粒子峰（Coherent Quasiparticle Peak），同时伴随着 $\pm 2\text{ eV}$ 附近的 Hubbard 子带（Hubbard Bands）演化。这一典型的 Mott 转变图像在多物理场系综平均下得到了极为精细的呈现。

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

为便于广大计算化学与固体物理同行复现本项工作，以下详细拆解其多尺度模拟管线（Pipeline），并提供关键配置指南。

3.1 总体多尺度模拟流程

+------------------+
| 1. Supercell     | --> 生成 2x2x2 钙钛矿超胞，利用群论对称性分析
|    Code          |     提取 22 个非等价团簇类别 (含有不同 Ti/V 占据)
+------------------+
         |
         v
+------------------+
| 2. VASP (DFT)    | --> 对 22 个超胞进行晶格常数与内部原子位置全弛豫
|    Relaxation    |     获得每个团簇的平衡总能量 E_j 和晶格常数 a_j
+------------------+
         |
         v
+------------------+
| 3. PLOVasp &     | --> 投影 V-3d 轨道，获取局部规范化 Wannier 投影算符
|    DFTTools      |     确定各超胞自适应能量区间 (Energy Window)
+------------------+
         |
         v
+------------------+
| 4. TRIQS/CTSEG   | --> 在 Kanamori-Hubbard 相互作用下自洽求解 DMFT 方程
|    (DMFT Solver) |     提取虚频格林函数 G(tau)，并获取零频谱权重 A(w=0)
+------------------+
         |
         v
+------------------+
| 5. EGQCA         | --> 输入弛豫能量和能谱数据，最小化混合自由能 DG
|    Code          |     计算团簇概率 p_j(x, T)，输出宏观性质系综平均值
+------------------+
         |
         v
+------------------+
| 6. Percolation   | --> 基于简单立方格子进行 site-percolation 模拟
|    Simulation    |     计算电学穿透概率 p_perc (判断宏观导电通道是否建立)
+------------------+

3.2 步骤 1：对称非等价超胞生成

使用开源工具 Supercell 进行随机位置替换与对称性分类：

  • 开源 Repo 链接：https://github.com/gloseira/supercell
  • 操作指令： 在 $\text{SrTiO}_3$ cubic $2 \times 2 \times 2$ 晶胞基础上，对 B 位阳离子位点进行所有可能的二元 $\text{Ti}_{1-x}\text{V}_x$ 替换。程序将自动输出 22 个非等价空间群对称结构，并给出各自的微观简并度 $g_j$（对应配置权重）。

3.3 步骤 2：DFT 结构弛豫（VASP）

对于 22 个团簇独立超胞，进行无对称性约束的完全弛豫（ISIF = 3）：

  • VASP 6.3 推荐 INCAR 关键参数配置：
    PREC   = Accurate
ENCUT  = 550.0  # 高切断能，避免晶格体积变化引起的 Pulay 应力
EDIFF  = 1.0E-7
EDIFFG = -1.0E-4 # 极严的力收敛阈值
ISMEAR = 0       # 高斯展宽，避免金属-绝缘体转变区域电荷振荡
SIGMA  = 0.05
LCHARG = .FALSE.
LWAVE  = .TRUE.  # 保留波函数，以便后续投影
ADDGRID= .TRUE.
  • 赝势选择：PAW_PBE 数据库中的 Sr_sv（带半芯态）、Ti_sv、V、O。

3.4 步骤 3 & 4：强关联能谱解算（TRIQS 平台）

使用基于 Python/C++ 的多体物理计算库 TRIQS：

关键配置脚本片段逻辑：

from triqs.applications.dft.converters.vasp import VaspConverter
from triqs.applications.dft.sumk_dft import SumkDFT
from triqs_ctseg import Solver
import numpy as np

# 1. 转换 VASP 输出的 Wannier 投影矩阵 (PROJCAR)
Converter = VaspConverter(archive_name='srtivo3_cluster_j.h5')
Converter.convert_dft_input()

# 2. 构建 DFT+DMFT 关联环路
SK = SumkDFT(archive_name='srtivo3_cluster_j.h5', use_dft_blocks=False)

# 3. 设置 Kanamori 参数 (在富钒端和过渡端对等适用)
U = 4.0   # eV
J = 0.65  # eV
beta = 40.0 # eV^-1 (对应约 290 K)

# 构建双体排斥相互作用矩阵
# (利用常规 Kanamori 形式)
# ... [设置相互作用张量 U_matrix] ...

# 4. 初始化 CT-QMC 杂化膨胀 solver
S = Solver(beta=beta, gf_struct=[('l_t2g', 3)]) # 考虑 t2g 的 3 个简并轨道

# 5. 自洽循环解算
for dmft_step in range(40):
    SK.symm_deg_gf() # 对称化格林函数
    S.G0_iw << SK.G0_iw # 更新杂化函数
    S.solve(h_int=h_int, n_cycles=1000000) # QMC 抽样
    SK.set_Sigma([S.Sigma_iw]) # 传回自能
    SK.dft_quantum_step() # 结合双重计数修正 (FLL)

  • 零频谱函数计算（公式 13 的直接实现）： 为了避免数值解析延拓（如 MaxEnt）带来的不确定性，在 QMC 收敛后，直接读取虚频（虚时间）格林函数在 $\tau = \beta / 2$ 处的数值：

    $$A(\omega=0) \approx -\frac{eta}{\pi} G\left(\tau = \frac{\beta}{2}\right)$$

    该公式在低温极限下是严格精确的，避免了人工调能级展宽的偏差。

3.5 步骤 5：GQCA 统计平均计算

  • 自研热力学核心解算器：作者使用了名为 EGQCA 的程序（可向通讯作者索取，或参考类似开源 Egqca 规范：https://github.com/EGQCA）。
  • 该脚本通过数值求解公式 10，在设定的 $T = 1800\text{ K}$ 下解出拉格朗日因子的等效实根 $\eta$，进而快速获取每个团簇的统计概率 $p_j(x, T)$，并对 $a_j$ 及 $A_j(\omega=0)$ 进行加权平均（公式 11）。

4. 关键引用文献与本工作局限性评述

4.1 关键文献追溯

  1. GQCA 理论奠基：
    • A. Sher, M. van Schilfgaarde, A.-B. Chen, and W. Chen, “Quasichemical approximation in binary alloys”, Phys. Rev. B 36, 4279 (1987).
    • 贡献：奠定了拟二元合金热力学配分函数和多元团簇统计分布的数学物理基础。
  2. $\text{SrTi}_{1-x}\text{V}_x\text{O}_3$ 局域超胞关联探究：
    • A. D. N. James, M. Aichhorn, and J. Laverock, “Composition-driven mott transition within $\text{SrTi}_{1-x}\text{V}_x\text{O}_3$”, Electronic Structure 6, 015010 (2024).
    • 贡献：证实了 MIT 对局域阳离子排布（$\text{Ti}/\text{V}$ 的短程有序）极其敏感，指出单一超胞无法等效热力学固溶体，启发了本工作系综平均框架的建立。
  3. 实验热力学关键证据：
    • J. Macías, A. A. Yaremchenko, et al., “Compromising between phase stability and electrical performance: $\text{SrVO}_3\text{-SrTiO}_3$ solid solutions…”, ChemSusChem 12, 240 (2019).
    • 贡献：详实报道了均匀固溶体在 $1100^\circ\text{C}$ 下彻底相分离为两相的物理化学事实，直接充当了本工作相图预测的关键 Benchmark 对象。

4.2 局限性深度评述

尽管该工作代表了强关联合金计算领域的顶尖水平，但从极度苛刻的学术视角来看，仍存在以下值得深入改进的局限性：

1. 静态团簇近似（Static Cluster Assumption）

GQCA 假设超胞团簇之间是“空间和能量相互独立”的，这意味着它完全忽略了长程静电交互作用（Madelung Potential）和弹性长程应变场（Long-range Elastic Interactions）。在实际钙钛矿合金中，$\text{V}$ 与 $\text{Ti}$ 的尺寸和局域电荷分布差异会导致长程应变场的自发协同调节。这种长程弹性各向异性无法通过简单的、孤立的 $2\times 2\times 2$ 超胞边界完全消解，可能导致对团簇概率 $p_j$ 的高估或低估。

2. 能谱-结构弛豫的非自洽性（Non-self-consistent Mott-elastic Coupling）

本项研究中的结构原子弛豫全部是在 DFT 级别完成的。在发生金属-绝缘体转变（MIT）时，局域电子排布发生剧烈突变（相干准粒子消亡，局域磁矩形成，库仑排斥能急剧攀升）。这种局域电子态的重构会反过来施加一个强大的“Mott 弹性能”（Mott-elastic driving force），彻底改变局域 $\text{V}-\text{O}$ 键长和八面体旋转角。由于未采用自洽的 $\text{DFT}+\text{DMFT}$ 全力弛豫（Force Relaxation），计算得到的局域团簇能量 $H_j$ 存在一定程度的偏差。

3. 渗流模型中金属/绝缘体二分法的粗糙性

在图 5 所示的渗流模拟中，作者将光谱权重 $A(\omega=0) < A_c = 0.1\text{ (eV}\cdot\text{f.u.)}^{-1}$ 的团簇强行定义为“绝缘体”，其余为“金属”。在真实的多相共存微观图像中，绝缘与金属边界是极其平滑且存在量子隧道效应（Quantum Tunneling）的。将复杂强关联体系简化为一维、无相互作用的经典 Site Percolation，虽然物理直觉清晰，但定量上丢失了弱局域化区间和介观尺度电荷输运的细节。

5. 补充：Mott-Anderson 协同物理与未来展望

5.1 Mott 物理与 Anderson 无序的微观协同机制

在这项工作中，GQCA 团簇平均框架揭示出了一个极其精妙的微观协同机制。 当 $x$ 较小时，$\text{V}$ 原子处于孤立状态，被 $\text{Ti}$ 原子包围。由于 $\text{Ti}^{4+}$ 的 $3d$ 轨道完全空置且能量较高，$\text{V}$ 轨道上的单个 $d$ 电子无法轻易向外跃迁。在微观能谱上，这表现为 $\text{V}-t_{2g}$ 局部子能带带宽 $W$ 的急剧收缩（在 $x=0.25$ 时收缩至仅有 $0.67\text{ eV}$，而在 $x=1$ 的纯 $\text{SrVO}_3$ 中为 $2.52\text{ eV}$）。 随着带宽 $W$ 的严重变窄，局域库仑排斥作用阻碍电荷跃迁的相对能力（即有效关联强度 $U/W$）被成倍放大，自发触发了 Mott 局域化，打开能隙（见图 6）。

与此同时，$\text{V}$ 原子的随机空间排布产生了显著的 Anderson 势能涨落，使不同局部区域的电学活性迥异。只有当 $x$ 达到一定临界值（$\sim 0.4$），这些局部已经克服 Mott 局域化并展现相干准粒子特征的金属化团簇相互接触、连通，并在空间上跨越宏观样品的渗流阈值时（见图 5 中的紫色实线），电导率才会发生骤增。这一图像完美地将 “局部 Mott 能隙闭合” 与 “空间 Anderson 渗流连通” 统一在一个极具说服力的物理框架下。

5.2 未来展望与方法推广

本工作建立的 $\text{GQCA}+\text{DFT}+\text{DMFT}$ 无参数预测方法，其重大价值不仅在于解决了 $\text{SrTi}_{1-x}\text{V}_x\text{O}_3$ 一个体系的争议，更为宽温区、强无序、强关联多组分体系的探索开辟了通用通路：

  1. 高熵氧化物（High-Entropy Oxides, HEOs）：HEOs 通常由 5 种及以上的过渡金属阳离子无序混配而成。其局部晶格无序和多体电荷关联极其显著。利用本工作发展的对称团簇分类配合极高统计精度的 GQCA，可在可承受的计算开销内预测 HEOs 独特的磁性转变、莫特转变及超导特性。
  2. 超导合金及铁基超导体中的等价掺杂：例如 $\text{BaFe}_2(\text{As}_{1-x}\text{P}_x)_2$ 等体系，掺杂直接作用于费米面本身。通过基于 GQCA 框架的自能和费米面平均重构，可自洽评估无序散射与超导电性配对振幅的演化趋势。
  3. 高通量无序物相筛选：将此技术封装入自动化高通量计算工作流（如 AiiDA 或 Fireworks 平台），能够大规模、快速筛选并设计具备特定 Mott 转变温度的智能外场响应窗口材料、关联透明导电氧化物（Correlated TCOs）以及低能耗 Mott 忆阻器件。