来源论文: https://arxiv.org/abs/2605.22522v1 生成时间: Jun 07, 2026 18:21

无公度势、强关联与超导的纠缠：混杂过渡金属二硫属化物 4Hb-TaS2 的三维莫尔物理学深度解析

0. 执行摘要

在二维范德华（vdW）材料中，通过人工转角或层间堆叠构建莫尔超晶格（Moiré Superlattices）已成为调控电子结构、诱导强关联相变（如 Mott 绝缘态、非常规超导、拓扑磁性）的黄金法则。然而，这种“莫尔工程”长期以来被局限于二维单层或双层体系中。对于三维块体材料，是否存在天然的、具有类似物理效应的“体莫尔势（Bulk Moiré Potentials）”，以及它如何与材料内部固有的多体关联（电子-电子相互作用、电荷密度波、超导）发生耦合，一直是凝聚态物理和量子化学领域悬而未决的前沿课题。

近期，来自魏茨曼科学研究所的 Jean C. Souza、Haim Beidenkopf 等人组成的联合研究团队，在一项具有里程碑意义的工作中，聚焦于混杂过渡金属二硫属化物（TMD）多型体 4Hb-TaS2。该材料由强关联的 1T-TaS2 层与超导金属 1H-TaS2 层交替堆叠而成。研究团队结合极低温扫描隧道显微镜/光谱（STM/STS）、密度泛函理论（DFT）自洽计算以及动力学平均场理论（DMFT）多体模拟，系统性地证明了 4Hb-TaS2 中存在由 1T 和 1H 层之间约 1.1% 的晶格失配（Lattice Mismatch）所引起的涌现无公度势（Emergent Incommensurate Potential）。

这一无公度势通过空间局部调控 1T-1H 层间距离，精细调控了层间杂化强度和电荷转移量，使原本因强电荷转移而退化的 1T 层 Mott 态局部“复活”，进入一个掺杂 Mott 机制调控的自屏蔽（Self-Screening）区间，展现出特征性的零偏压电导峰（ZBCP）及不对称的 Hubbard 卫星带。更为惊人的是，研究发现该无公度局域势与材料的体超导转变（$T_c \approx 2.6\text{ K}$）存在强烈的非平凡耦合——在超导态下，无公度势的晶格“呼吸”相变被显着调制。这一成果不仅颠覆了此前“硒杂质导致电荷转移抑制”的传统认知，更将莫尔工程的设计理念成功拓宽至三维块体强关联量子材料中。

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

1.1 核心科学问题

混杂过渡金属二硫属化物 4Hb-TaS2（空间群 $P6_3/mmc$）是一种极为独特的天然异质结。在单层极限下，1T-TaS2 是典型的 Mott 绝缘体，其 $\sqrt{13} \times \sqrt{13}$ 的电荷密度波（CDW）重构锁定了单粒子能带中的平带，在半填充时产生强关联 Mott 间隙；而 1H-TaS2 则表现为金属性并伴随超导电性（$T_c \approx 0.8\text{ K}$，在块体 4Hb 中提升至 $2.7\text{ K}$）。

当这两种极性截然不同的层交替堆叠形成 4Hb-TaS2 时，1T 层的电子会大量向金属性的 1H 层转移（平均电荷转移量约为 $92\% \pm 11\%$ 电子/CDW 单元），导致 1T 层的 Mott 状态由于过度空穴掺杂而退化（Depleted Mott State）。然而，实验上在 1T 终止面上却总是能观测到约 $11\% \pm 1\%$ 的 CDW 位点保留了完整的下 Hubbard 带（LHB）、上 Hubbard 带（UHB）以及费米能级处的零偏压电导峰（ZBCP，常被理解为 Kondo 共振或掺杂 Mott 准粒子峰）。

长期以来，学术界普遍将这大约 11% 的“未退化（Undepleted）”位点归因于材料生长过程中不可避免引入的硒（Se）代位缺陷。本工作核心要解决的科学问题是：这些未退化的强关联位点究竟是杂质诱导的平庸缺陷效应，还是由 1T/1H 层间本征无公度晶格失配所产生的涌现体莫尔势（Incommensurate Moiré Potential）所主导的本征量子物态？

1.2 理论基础

1. 无公度结构与超空间群理论：当两个具有不同晶格常数的层叠加时，若其晶格矢量的线性组合无法用有理数表示，系统便失去平移对称性。在数学上，可以通过将低维无公度结构投影到高维（如 4D 或 6D）周期性超空间（Superspace）中进行严格描述。在本体系中，1T 和 1H 的晶格常数分别为 $3.38\text{ \AA}$ 和 $3.342\text{ \AA}$，存在 1.1% 的晶格失配，形成一个跨越约 550 个 CDW 位点的巨型莫尔超晶格。
2. 层间电荷转移与局部能带杂化：电荷转移量（CT）强烈依赖于层间阴离子（S 2p 轨道）与阳离子（Ta 5d 轨道）的重叠程度。晶格失配导致沿面内空间不同的局部堆叠位置（Registries），其局域层间距离 $d_{\perp}$ 发生周期性起伏。根据紧束缚模型，杂化矩阵元 $V_0$ 随 $d_{\perp}$ 呈指数衰减： $$V_0(d_{\perp}) \propto e^{-\beta d_{\perp}}$$ 因此，层间局域晶格重合度的微小变化将导致电荷转移量和局域电子自能（Self-energy）的剧烈空间调制。
3. 周期性安德森模型（PAM）与多体关联：为了描述 1T 层（局部化 d 轨道，具强库仑排斥 $U$）与 1H 层（非局域传导 c 电子）的耦合，建立周期性安德森模型： $$H = H_{1H} + H_{1T} + H_V + H_U$$ 通过动力学平均场理论（DMFT），将格点自能自洽映射为单杂质安德森模型（SIAM），从而在动力学层次精确捕获零温附近的相干准粒子峰（Kondo 尺度）和高能 Hubbard 分支。

1.3 技术难点

• 计算尺度的灾难：一个包含 1.1% 晶格失配的完整 4Hb-TaS2 莫尔胞包含有数千个原子，且包含 $\sqrt{13} \times \sqrt{13}$ 的 CDW 面内重构。直接对其进行常规的 DFT+U 计算，乃至 DMFT 求解，在计算资源上是完全不可行的。
• 谱学特征的微细分辨：实验上 ZBCP 处于零偏压附近，半峰全宽（FWHM）仅为几个毫电子伏特（meV），这对 STM 谱学测量的信噪比以及极低温（1.1 K）下的热展宽控制提出了极其苛刻的要求。
• 动态自洽双重计数（Double Counting）修正：在 PAM 的 DMFT 求解中，由于电荷转移是由 DFT 预先给出并需要在 DMFT 中自洽维持，如何动态调整双重计数势 $\mu_{dc}$ 以精确匹配局域占据数 $n_d = 1 - \text{CT}$，是多体算法收敛的难点。

1.4 方法细节

为了克服上述计算瓶颈，研究团队巧妙地采用了局域晶格重合近似（Local-Registry Approximation）。具体步骤如下：

1. 结构采样与第一性原理计算（DFT）： 使用 VASP 软件包，构建包含单层 1T-TaS2 和单层 1H-TaS2 的双层异质结模型（引入面内 $\sqrt{13} \times \sqrt{13}$ 的 CDW 面内晶格调制）。通过人为施加微小的面内平移矢量，在莫尔胞内采样四个具有高对称性或代表性的局域堆叠构型：

   • 稳定堆叠（Square, $\square$）：系统能量最低的构型，层间距最小，杂化最强，对应最大电荷转移（CT $\approx 0.92$）。
   • 中间堆叠（Star, $\star$）：介于稳定与亚稳态之间。
   • 亚稳堆叠 I（Pyramid, $\triangle$）：层间距中等。
   • 亚稳堆叠 II（Inverted Pyramid, $\nabla$）：层间距最大（达 $6.9\text{ \AA}$），杂化最弱，对应电荷转移受阻最严重的构型（CT $\approx 0.57$）。

   使用 PBE 泛函和 DFT-D3 色散修正进行结构弛豫，在 $Ta\ 5d$ 轨道上施加局域库仑作用（1T 层 $U=1.76\text{ eV}$，1H 层 $U=2.82\text{ eV}$，采用 Dudarev 方法）。

2. 紧束缚参数化（Tight-Binding Extraction）： 从 DFT 计算的局域构型能带结构中，提取出 $H_{1H}$（包含 13 个金属性能带）和 $H_{1T}$（单能带平带模型）的层内跃迁积分，并确定层间杂化项 $H_V$ 的耦合标量 $V_0$。

3. 多体关联模拟（DMFT）： 在提取出的四个局域紧束缚哈密顿量基础上，加入局域 Hubbard 作用项（$U_{1T} = 100\text{ meV}$，对应于重构 CDW 单元的屏弊后有效相互作用强度）。 采用连续时间杂化膨胀量子蒙特卡洛（CT-HYB）算法作为杂质求解器（使用开源多体计算包 w2dynamics）。模拟温度跨越 $230\text{ K}$ 到极低温 $8\text{ K}$。在每一步 DMFT 迭代中，动态调节双重计数势，强制使 $d$ 轨道的局域占据数严格匹配 DFT 给出的各向异性电荷转移量，最终自洽输出虚轴自能 $\Sigma(i\omega_n)$，并通过最大熵方法（MaxEnt）进行解析延拓，得到实轴上的局域态密度（LDOS）。

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

2.1 关键 Benchmark 体系：Se 掺杂浓度依赖性试验

为了彻底澄清“未退化 Mott 态由 Se 缺陷导致”的既往假说，研究团队设计了一组极具说服力的对照实验体系：$4\text{Hb-TaS}_{2-x}\text{Se}_x$ 晶体，其中 $x = 0\%, 1\%, 2\%$。这一设计的物理精妙之处在于：引入 Se 杂质不仅会产生局域杂质势，更会系统性地增大 $c$ 轴层间距（$c/a$ 比值增加），从而在全局上削弱层间杂化。

实验观测数据：

• ** topography（形貌图）**：在 $x=0\%, 1\%, 2\%$ 的所有样品中，1T 终止面均展现出完美的 $\sqrt{13} \times \sqrt{13}$ CDW 超结构。通过偏置电压成像，可以清晰区分出“亮位点（Bright Sites，对应未退化强关联态）”和“暗位点（Dim Sites，对应电荷退化 Mott 态）”。
• 未退化位点比例的定量分析（Fig. 1f）： 对三种样品的独立解理面进行大规模统计，亮位点的占比 $n_{occ}(x)$ 随 Se 掺杂量 $x$ 呈线性增长关系： $$n_{occ}(x) \approx 0.08 + 3x$$
  • 关键结论：当 $x \to 0$ 时，公式的外推截距不为零，而是指向约 $8\%$ 的本征占比。在完全不含 Se 的纯 $4\text{Hb-TaS}_2$ 样品中，实验依然稳定测得含有约 $8\%$ 的未退化强关联位点（Fig. 1e）。这极其有力地证明了未退化 Mott 态存在非平庸的、不依赖于任何外在杂质的本征起源——即由 1T/1H 晶格失配产生的天然无公度势。

2.2 DFT 与 DMFT 关键计算数据分析

研究团队通过对四个局域堆叠构型的系统计算，给出了层间距离、电荷转移以及多体关联之间的定量链条关系。

理论计算的核心发现：

1. 能带演化与晶格呼吸（Fig. 2e）：当系统从基态 $\square$ 移动到亚稳态 $\nabla$ 时，层间距离由于原子排斥被迫增大了约 $1.2\text{ \AA}$。这种局域膨胀极大地削弱了层间 $Ta-S-Ta$ 轨道的重叠，导致电荷转移量从近乎饱和的 $0.92\ e^-$ 锐减到 $0.57\ e^-$。这表明，1.1% 的晶格失配在空间中形成了一个强烈的、波动的层间势垒。
2. 多体谱线的不对称性——自屏蔽的确凿证据（Fig. 3a, 3c）：
   • 在经典 Kondo 屏蔽模型中，若 ZBCP 源自 1T 局域自旋与 1H 传导电子的杂化，在无电荷转移（即半填充）时，其 LHB 和 UHB 应该相对于费米能级对称分布。
   • 然而，DMFT 计算所得的光谱（Fig. 3a）以及实验谱学统计（Fig. 3c）均显示出显着的能带不对称性：LHB 距离费米面较近，而 UHB 被推向更高的能量（能量重心差值达 $\approx 80\text{ meV}$）。
   • 物理阐释：DMFT 完美复现了这一不对称性，证明其本质上是由于电荷转移受阻后，1T 层仍处于非半填充（空穴掺杂）的 Mott 状态。在这种状态下，局域磁矩的屏蔽主要是通过 1T 层内自有的空穴流动进行相干自屏蔽（Self-Screening），而非平庸的层间 Kondo 耦合。这一结论澄清了长期以来关于该体系是“Kondo 晶格”还是“掺杂 Mott 绝缘体”的理论争端。

2.3 莫尔势与超导的纠缠数据 (Fig. 4)

实验上，研究团队在极低温下探测了 CDW 位点占有率随温度以及磁场的变化，揭示了无公度莫尔势与块体超导性之间的强耦合：

1. 温度驱动的晶格呼吸（Fig. 4a, 4b）： 在零磁场下，当样品从 $4.2\text{ K}$（正常态）冷却至 $1.1\text{ K}$（超导态）时，约有 $40\%$ 的 CDW 位点改变了其占有构型。这种变化并不是杂乱无章的，而是沿特定的莫尔域界（Moiré Domains）集中发生。其中，$68\%$ 的相变仅涉及最邻近的高对称性构型跃迁（$\Delta n_{occ} = \pm 1/3$），高度契合无公度势中相邻稳定/亚稳格点之间的热力学竞争。
2. 超导态对晶格呼吸的抑制（Fig. 4c, 4d）：
   • 当施加垂直 $c$ 轴的外加磁场 $\mu_0 H_{\parallel c} = 1.7\text{ T}$（该强度已超过 $4\text{Hb-TaS}_2$ 垂直方向的超导临界磁场 $H_{c2} \approx 1.1\text{ T}$，系统退回正常金属态）时，这种温度驱动的局域构型跃迁被近乎完全抑制（Fig. 4c）。
   • 相反，当施加面内磁场 $\mu_0 H_{\perp c} = 1.7\text{ T}$（远低于该系统庞大的面内超导临界磁场 $H_{c2}^{\parallel} \approx 12\text{ T}$，系统仍保持强超导相干性）时，构型的空间跃迁行为得以完全恢复（Fig. 4d）。
   • 深层物理机制：超导态与无公度莫尔势之间存在相互竞争和调制。在超导相干态下，由于超导能隙（$\Delta \approx 0.6\text{ meV}$）的打开，局域电子态密度在费米能级附近发生重构，这导致超导凝聚能（Condensation Energy）与无公度势的静电能发生微妙的竞争。超导波函数的空间调制（如潜在的对密度波 Pair Density Wave, PDW）与莫尔超晶格的空间不均匀性产生共振，促使了局域堆叠构型在低温下的重新排布。这一发现提供了三维莫尔超导材料中“电子结构-晶格失配-配对对称性”三者纠缠的直接谱学证据。

3. 代码实现细节、复现指南与开源链接

为了方便量子化学与凝聚态计算领域的科研人员复现本工作的核心理论模拟，以下给出基于 DFT+U 参数提取以及使用 DMFT 求解周期性安德森模型（PAM）的完整技术路线。

3.1 所需软件包及开源仓链接

1. 第一性原理计算：VASP (Vienna Ab initio Simulation Package) (商业软件，需授权)
2. 紧束缚能带投影：Wannier90 (开源，提取最大局域化 Wannier 函数) GitHub Repo
3. 多体动力学平均场计算：w2dynamics (开源，基于连续时间量子蒙特卡洛算法的强关联杂质求解器) GitHub Repo
4. 数据后处理与解析延拓：TRIQS (Toolbox on Research in Interacting Quantum Systems) 及其 MaxEnt 模块 GitHub Repo

3.2 理论模型参数化与复现工作流

步骤 1：第一性原理局域构型构建与弛豫 (VASP)

针对无公度体系的四个代表性局域堆叠（如基态 Square $\square$ 和亚稳态 Inv-Pyramid $\nabla$），构建双层 $1\text{T}/1\text{H}-\text{TaS}_2$ 晶胞。由于要考虑 1T 层的 $\sqrt{13} \times \sqrt{13}$ CDW 畸变，晶胞包含 26 个 Ta 原子和 52 个 S 原子。

关键 VASP 输入参数示例 (INCAR):

# 基础电子结构与自洽收敛控制
ENCUT = 400.0       # 平面波截断能
ALGO = Normal
NELM = 100
EDIFF = 1E-7        # 极高的电子能量收敛阈值

# 结构弛豫控制
IBRION = 2
NSW = 200
EDIFFG = -1E-3      # 严格的受力收敛标准：原子受力 < 1 meV/Å
ISIF = 2            # 仅弛豫内部原子坐标，保持晶胞大小（对应局域刚性重合近似）

# 范德华力色散修正 (关键！用于准确捕获层间距变化)
IVDW = 11           # 激活 DFT-D3 修正

# DFT+U 强关联修正 (Dudarev 方案)
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 2  2        # 对 1T-Ta 和 1H-Ta 的 d 轨道施加 U
LDAUU = 1.76  2.82  # 根据文献 [48] 参数化
LDAUJ = 0.00  0.00
LMAXMIX = 4         # 对 d 轨道元素，此项必设为 4，保证电荷密度混叠的稳定性

步骤 2：投影最大局域化 Wannier 函数 (Wannier90)

利用 Wannier90 投影费米能级附近的能带，提取高精度紧束缚哈密顿量。其中：

• 1H 轨道：选择 Ta $5d_{z^2}$ 轨道投影（共 13 个能带，代表非关联导带 c）。
• 1T 轨道：选择处于费米能级处的孤立平带进行投影（共 1 个能带，代表关联轨道 d）。

通过该投影，可直接获取能谱矩阵元，作为 DMFT 输入的 local_hamiltonian。

步骤 3：多体强关联计算 (w2dynamics)

在 w2dynamics 中配置周期性安德森模型。利用自洽调整的 $\mu_{dc}$ 保证各构型对应的目标占据数 $n_d = 1 - \text{CT}$（例如对于 $\nabla$ 构型，设定 $n_d = 1 - 0.57 = 0.43$）。

w2dynamics 输入配置文件片段 (input.ini):

[General] 
# 模拟物理参数
beta = 150.0             # 对应极低温极限 (k_B T ≈ 0.67 meV)
U = 0.1                  # 1T 层的有效局域排斥相互作用 U_d = 100 meV
mu_dc = -0.05            # 初始双重计数势势能

[Hamiltonian]
# 导入 Wannier90 输出的紧束缚 Hamiltonian
model = "periodic_anderson"
tb_file = "tas2_tb_parameters.json" 

[Solver]
# 量子蒙特卡洛抽样参数
solver = "cthyb"         # 杂化膨胀算法
measure_gtau = .true.
measure_gleg = .false.
warmups = 10000
qmc_steps = 1000000      # 确保高能 Hubbard 卫星带自能项的高频精度

[SelfConsistency]
max_loops = 25           # 自洽迭代上限
conv_threshold = 1E-5    # 收敛判定阈值
target_occupation = 0.43 # 动态调节 mu_dc 以锁定目标电荷转移量

步骤 4：最大熵解析延拓 (TRIQS/MaxEnt)

量子蒙特卡洛计算是在虚时间/虚频轴 $G(i\omega_n)$ 上进行的。为了与 STM 测量的微分电导谱 $dI/dV$ 进行对比，需要将其延拓至实频轴上得到谱函数 $A(\omega)$。

Python 脚本片段 (maxent_postprocess.py):

from triqs.gf import *
from triqs_maxent import *
import numpy as np

# 1. 载入 w2dynamics 自洽收敛的虚轴格林函数
g_iw = GfImFreq(indices=[0], beta=150.0)
g_iw.from_file("tas2_dmft_results.h5", "g_iw_last_loop")

# 2. 初始化最大熵求解器
tm = TauMaxEnt(chi2_target=1.0)

# 3. 设置实轴能量网格 (从 -300 meV 到 +300 meV)
omega_grid = np.linspace(-0.3, 0.3, 1000)

# 4. 执行解析延拓，得到谱函数 A(omega)
result = tm.run(g_iw, omega_grid=omega_grid)
spectral_function = result.get_spectral_function()

# 5. 保存并绘图
np.savetxt("tas2_spectral_function.dat", spectral_function)

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

4.1 关键引用文献深度关联分析

本研究立足于深厚的物理学脉络。以下五篇文献对理解本工作起到了基石作用：

1. McMillan, W. L. Phys. Rev. B 14, 1496 (1976) [[1]]：电荷密度波（CDW）无公度-公度相变理论的开山之作，奠定了体系中 $\sqrt{13} \times \sqrt{13}$ CDW 调制与点阵畸变相互作用的经典唯象框架。
2. Wang, Y. et al. Nat. Commun. 11, 4215 (2020) [[21]]：首次系统实验论证了单层 1T-TaS2 随着温度降解展现出从能带绝缘体到强关联 Mott 绝缘体的转变，提供了 1T 层孤立平带的标准单粒子参数。
3. Nayak, A. et al. Proc. Natl. Acad. Sci. USA 120, e2304274120 (2023) [[38]]：首次在 4Hb-TaS2 中报道了一阶量子相变现象，本工作正是在此基础上，将该相变的物理根源从平庸的“Se 缺陷阻碍”修正为本征的“无公度势晶格呼吸”。
4. Crippa, L. et al. Nat. Commun. 15, 1357 (2024) [[46]]：构建了 1T/1H 异质结的理论大厦，论证了多体自能中 Kondo 物理与掺杂 Mott 物理的理论边界，为本工作 DMFT 模型的设计直接提供了 PAM 参数模板。
5. Wallerberger, M. et al. Comput. Phys. Commun. 235, 388 (2019) [[70]]：详细阐述了本工作所用强关联求解器 w2dynamics 的算法架构、自洽算法细节及蒙特卡洛误差控制方法。

4.2 本工作局限性之独家批判性评论

尽管本工作在实验技术和多体建模的结合上堪称典范，但站在量子化学与凝聚态计算的严苛视角，其理论框架仍存在以下不容忽视的局限性与改进空间：

1. 局域晶格重合近似（Local-Registry Approximation）的局限性

• 批判：该工作在 DFT 计算中假设四个堆叠状态是彼此完全隔离的刚性结构。然而，在真实的天然无公度莫尔超胞中，不同构型之间存在着连续的剪切应变（Shear Strain）和晶格重构（Lattice Reconstruction）。
• 物理缺失：这种“刚性切片”式的处理忽略了晶格滑移界面的拓扑重构（如莫尔孤子，Moiré Solitons）以及层间剪切力学对电荷转移梯度的平滑效应。真实系统的电荷转移并不是在四个离散值之间跳跃，而是呈空间连续调制分布。这一近似会导致计算出的亚稳态能量偏高，且高估了相变的突变性。

2. 双层模型（Bilayer Model）对三维体效应（Bulk Effects）的过度简化

• 批判：4Hb-TaS2 是一个极其复杂的三维散装（Bulk）系统，沿 $c$ 轴方向 1T 和 1H 是无限交替重叠的。然而，本工作中的所有 DFT+U 计算以及紧束缚参数投影，均是在**简化的 1T/1H 双层异质结（Bilayer）**上进行的。
• 物理缺失：在三维晶体中，1T 层夹在上下两个 1H 层之间，其所受到的静电屏蔽和电荷转移是双向的。双层模型会显着削弱电荷转移的层间屏蔽通道，高估 1T 层内部的局域库仑排斥强度。这或许正是计算中需要人为调小 $U$ 值（$U_{eff} = 100\text{ meV}$，远低于单层 1T-TaS2 真实的局域库仑作用）的深层根源。

3. 有效局域排斥作用参数 $U$ 的唯象定位

• 批判：在 DMFT 计算中，将 CDW 单元（包含 13 个 Ta 原子）等效收缩为一个单一的活性格点，并赋予其 $U = 100\text{ meV}$。这种处理虽然在凝聚态唯象模型中十分流行，但在量子化学视角下，它极度粗糙地抹杀了 CDW 单元内部 13 个 Ta 原子之间的局域电荷涨落和轨道杂化细节。如果能在未来结合约束性随机相位近似（cRPA），从头算确定这一有效 $U$ 值，理论说服力将大幅跃升。

4. 对超导调制的唯象性解读

• 批判：文章指出超导态会调制莫尔势和 CDW 占有率，但未能给出微观的定量计算。DMFT 模型中完全不包含超导配对自能项，超导仅作为实验现象被外加讨论。未能建立起一个统一包含超导序参量、电荷密度波和多体自能的空间自洽多体模型，是本工作理论部分的终极遗憾。

5. 补充探讨、领域展望与前沿碰撞

5.1 Se 代位掺杂物理机制的再审视

本工作极具学术冲击力的一点在于重新阐释了 Se 掺杂的作用。过去的研究（如文献 [[38, 39]]）简单地认为，Se 杂质作为“化学障碍”直接阻碍了局域的层间电荷转移。而 Souza 等人通过精细数据指出：

• 真实的物理图像：Se 原子的共价半径（$1.16\text{ \AA}$）显着大于 S 原子（$1.02\text{ \AA}$）。当 Se 替代 S 后，主要的效应是由于空间位阻效应导致晶格在 $c$ 轴方向全局发生膨胀。
• 力学-电学耦合：这种膨胀系统性地拉大了 1T 与 1H 层的平均间距，使整个材料内部的层间杂化整体削弱。其结果是，原本属于稳定构型（$\square$）的区域，其层间距也被迫拉大到接近亚稳态的临界值，从而使本征莫尔超晶格中原本能量偏高的“未退化强关联亚稳态”（如 $\nabla$ 构型）的能量被相对拉低，在统计学上表现为亮位点比例 $n_{occ}$ 随 Se 浓度线性激增。 这是一种优雅的**“全局形变调制局域关联物态”**的物理机制，而非简单的“局域散射局域缺陷”机制。

5.2 掺杂 Mott 自屏蔽与经典 Kondo 效应的本质区别

在凝聚态理论中，4Hb-TaS2 长期被视为研究“重费米子 Kondo 晶格”的天然实验室：1T 层的定域磁矩与 1H 层的传导电子发生反铁磁杂化，并在费米面产生 Kondo 共振。然而，本工作通过 DMFT 给出的光谱特征，强有力地支持了掺杂 Mott 绝缘体自屏蔽图像。这两者存在本质的量子化学区别：

【传统 Kondo 屏蔽图像】
  1T 定域自旋 <======[ 层间强杂化 V_0 ]======> 1H 传导电子 (c)
  物理后果：能带高度对称，ZBCP 的形成极其依赖层间电荷的动力学共振交换。

【本工作：掺杂 Mott 自屏蔽图像】
  1T 定域自旋 <======[ 1T 层内空穴流 ]======> 1T 掺杂定域空穴
  物理后果：由于层间高度脱耦（d_⊥ 增大至 6.9 Å），1T 层局部退化为独立单层。
            由于 CT=0.57，1T 层处于非半填充状态，利用其内部空穴多体关联形成准粒子重构。
            能带呈现极强的不对称性（UHB 与 LHB 高度失衡）。

这一物理图像的修正，对于指导如何在该类材料中寻找非费米液体、分数量子态等非常规物态具有决定性的指导意义。

5.3 展望：涌现“体莫尔物理学”（Bulk Moiré Physics）的诞生

长期以来，“莫尔物理”与“二维材料转角异质结”几乎划等号，这意味着极其低下的样品产出率（依赖微纳撕裂堆叠工艺，无法大规模量产）。

4Hb-TaS2 这一典型案例向我们昭示：在特定的天然块体化合物中，只要存在层间晶格失配（或不同的面内自旋/电荷调制周期），就能本征地涌现出跨越数千原子的体莫尔势。

这开辟了全新的**“三维体莫尔材料体系”**。未来的研究方向将极富吸引力：

1. 设计新型 misfit 混合晶体：例如通过气相输运法（CVT）生长具有无公度堆叠的 $(\text{BiS})_{1.2}\text{TaS}_2$、$(\text{LaS})_{1.14}\text{NbS}_2$ 等异质结构，在三维体材料中实现超高密度的莫尔平带调制。
2. 利用超声或高压精细调控体莫尔势：由于无公度势的能量景观极其平坦（堆叠能差仅在几十 meV 数量级），使用微弱的外部物理场（如吉帕级的流体静力压，或外加压电应变）即可戏剧性地改变层间距，从而在全球范围内实现 Mott 态的“淬灭”与“复活”的可逆调控。
3. 莫尔超导相干性的相干调控：本工作展示了磁场可以抑制或恢复莫尔势的局域跃迁。在未来，若能结合极低温微区约瑟夫森光谱（Josephson Scanning Tunneling Spectroscopy, JSTS），将可以直接在实空间绘制出超导库珀对密度波在无公度莫尔超晶格中的空间波动图谱，彻底揭开强关联、无公度势与超导电性三者纠缠的世纪之谜。