14-kohn-sham 问答笔记¶

源文件：14-kohn-sham.md

Q1：能量泛函 $E[\rho] = T + V_{ext} + J + E_{xc}$ 每一项的物理含义¶

原文："$E[\rho] = T[\rho] + V_{ext}[\rho] + J[\rho] + E_{xc}[\rho]$"

项	数学形式	物理含义	单位/量级（Si 原胞参考）
$T[\rho]$ 动能泛函	难写成 $\rho$ 的显式形式；KS 中替换成 $T_s = \sum_i \int \phi_i^* (-\tfrac{1}{2}\nabla^2) \phi_i$	N 个电子运动总动能的期望值（量子力学意义上）	+100 eV/原胞（正值，"动"总是正能量）
$V_{ext}[\rho]$ 外势能	$\int V_{ext}(\vec{r}) \rho(\vec{r}) d^3r$（简单积分）	电子云与原子核之间的静电吸引能——电子被核"吸引"获得能量	−300 eV/原胞（负值，吸引性主导）
$J[\rho]$ 经典库仑能 / Hartree 能	$\frac{1}{2} \iint \frac{\rho(\vec{r})\rho(\vec{r}')}{\\|\vec{r}-\vec{r}'\\|} d^3r d^3r'$	把电子密度当经典电荷分布算它的自库仑排斥（忽略了量子性）	+60 eV/原胞（正值，电子互斥）
$E_{xc}[\rho]$ 交换-关联能	真实多体贡献与 $T_s, J$ 等的差	两类量子效应：交换（同自旋费米子回避）+ 关联（所有电子动态躲避）+ $T - T_s$ 动能修正	−90 eV/原胞（负值，降低能量）

相加：$T + V_{ext} + J + E_{xc} \approx 100 - 300 + 60 - 90 = -230$ eV/原胞，和 Si 实际总能量符合。

关键观察： - $|V_{ext}|$ 最大（电-核吸引主导） - $T$ 和 $|E_{xc}|$ 也是 100 eV 级 - $J$（经典库仑）次之 - 相互之间大量抵消，最终 $E_{total}$ 只有 -230 eV

这正是 DFT 需要高精度的原因——几百 eV 的大数相减，要求每项误差 < 几 meV 才能给出 meV/atom 的总能量精度。

Q2：为什么 Thomas-Fermi 误差巨大¶

原文："早期 Thomas-Fermi 用 $T^{TF}[\rho] \propto \int \rho^{5/3} d^3r$ 硬写，但误差巨大——甚至不能束缚原子成键"

Thomas-Fermi 模型（1927）：

\[T^{TF}[\rho] = \frac{3}{10}(3\pi^2)^{2/3} \int \rho^{5/3}(\vec{r}) d^3r\]

来自均匀电子气（uniform electron gas）的动能密度——假设电子像自由电子气一样填充费米海。

为什么这个假设对真实原子 / 分子 / 固体误差巨大¶

局域近似太粗
TF 的动能"每个点只看当地密度"，完全忽略密度的梯度和曲率
真实原子：核附近密度剧变（内层电子抱核），远处平滑（价电子延伸）
动能对"密度变化快不快"极其敏感——梯度信息关键，TF 丢掉了
不能描述壳层结构
TF 算出的原子密度是光滑球对称——没有 1s / 2s / 2p 壳层
现实：电子层层排布（Bohr 模型 → 量子力学轨道）
没有壳层就没有化学键、没有周期表规律
Teller 不稳定定理（1962）
严格证明：Thomas-Fermi 理论不能让任何分子稳定
即：H₂ 分子在 TF 下总能量高于两个分离 H 原子的能量之和
→ TF 认为 H₂ 会自动解离成两个 H 原子
这直接否定了 TF 作为化学方法的可能性
数值精度
对重原子（如 Hg）动能偏差 ~10%，绝对值 ~1 keV
DFT 精度需求：~1 meV/atom
差 $10^6$ 倍——TF 完全不可用

根本问题¶

动能对密度的依赖是强非局域的。"相同密度但不同节点/振荡模式"的波函数，动能可以相差几 eV/电子。TF 只看局域密度，等于把这些量子结构信息全丢了。

加上梯度修正（TF + $|\nabla\rho|^2$ 项）也只能小幅改善——核心问题是纯密度泛函捕不到动能的精细结构。

Kohn-Sham 的妙招¶

既然 $T[\rho]$ 不好直接写，就保留 N 个单粒子轨道 $\phi_i^{KS}$ 用来算动能：

\[T_s = \sum_i \int \phi_i^{KS*} \left(-\tfrac{1}{2}\nabla^2\right) \phi_i^{KS} d^3r\]

这个 $T_s$ 可以精确算（对给定 KS 轨道只是个积分）。真实 $T$ 与 $T_s$ 的差是个小量（~1 eV 级），塞到 $E_{xc}$ 里。

KS 的核心洞察：宁可多保留 N 个 3D 辅助对象（轨道），也要换取动能的显式可算性。

Q3：$E_{xc}$ 是什么？和其他项的量级关系¶

原文："$E_{xc}[\rho]$ 精确形式未知"

$E_{xc}$ 的组成¶

$E_{xc}$ 按物理来源分成三部分：

交换能 $E_x$（exchange）
费米子反对称约束导致同自旋电子互相回避（Pauli hole）
对 HF 这一项可以精确算（$-\frac{1}{4}\sum_{ij} \iint \frac{\phi_i^*(\vec{r})\phi_j^*(\vec{r}')\phi_i(\vec{r}')\phi_j(\vec{r})}{|\vec{r}-\vec{r}'|}$）
量级：~总能量 10–30%
关联能 $E_c$（correlation）
真实电子的库仑斥力导致动态躲避（correlation hole）
HF 缺失这一项
量级：~总能量 1–5%
动能关联修正 $T_c = T - T_s$
真实相互作用动能 vs 虚构非相互作用动能的差
量级：~几 eV/atom

量级对比表（He 原子最清晰）¶

He 原子是最简单的多体系统（2 电子），常作为 benchmark：

能量项	值（eV）	占精确总能量
$T$（真实动能）	+78.87	99.8% 正
$V_{ext}$（核吸引）	−199.45	252% 负（主导）
$J$（经典库仑）	+34.03	43% 正
$E_x$（交换）	−27.89	35%（大！）
$E_c$（关联）	−1.14	1.4%（小但化学决定）
$E_{xc}$ 总	−29.03	37%
$E_{total}$（精确）	−79.01	100%

几个重要观察¶

$E_x$ 本身大（约 1/3 总能量）——但 HF 能精确算，不是 DFT 的痛点
$E_c$ 虽小（1–5%）——但化学键能、反应活化能都在这个量级，决定性
DFT 的挑战：$E_{xc} = E_x + E_c$ 作为密度泛函（不用轨道）没精确形式，必须近似
XC 近似泛函的差异：LDA vs PBE vs SCAN 给出的 $E_{xc}$ 本身可能差 10 eV/atom，但总能量差通常只几百 meV/atom——大量相互抵消

在 DFT 实际计算的角色¶

总能量 E = T_s[{φ}] + V_ext[ρ] + J[ρ] + E_xc[ρ]
          ↑         ↑           ↑       ↑
          精确（轨道） 精确（积分）  精确（库仑） 近似（LDA/PBE/...）

精度瓶颈全在 $E_{xc}$——前三项都能算到 $10^{-8}$ 精度，$E_{xc}$ 近似误差决定了最终精度。

对 benchmark 的影响¶

磁性氧化物（NiO, MnO）：DFT 失败就是因为这类体系 $E_{xc}$ 近似严重偏离（关联太强）
半导体带隙：$E_{xc}$ 的 "derivative discontinuity" 缺失导致带隙低估 30–50%
范德华力：标准 $E_{xc}$ 抓不住长程关联，需要加色散修正（DFT-D3）

这些都是化简 5 的范畴。

Q4：为什么 KS 方程写成特征值等式的形式¶

原文："$\left[-\tfrac{1}{2}\nabla^2 + V_{eff}(\vec{r})\right] \phi_i^{KS}(\vec{r}) = \varepsilon_i \phi_i^{KS}(\vec{r})$"

因为 KS 方程在数学结构上就是一个单粒子薛定谔方程——标准形式 $\hat{H}\psi = E\psi$ 的特征值问题。

从变分到特征值¶

KS 方程的来源：最小化能量 + 轨道正交约束。

最小化目标（KS 形式能量）： $$E[\{\phi_i\}] = T_s[\{\phi_i\}] + \int V_{eff}(\vec{r}) \rho(\vec{r}) d^3r + \text{（与 } \phi_i \text{ 无关的项）}$$

约束：$\langle \phi_i | \phi_j \rangle = \delta_{ij}$

拉格朗日乘子法： $$\mathcal{L}[\{\phi_i\}; \{\varepsilon_{ij}\}] = E - \sum_{ij} \varepsilon_{ij} (\langle \phi_i | \phi_j \rangle - \delta_{ij})$$

对 $\phi_i^*$ 变分，设为 0： $$\hat{h}_{KS} \phi_i = \sum_j \varepsilon_{ji} \phi_j$$

其中 $\hat{h}_{KS} = -\tfrac{1}{2}\nabla^2 + V_{eff}$。

做幺正变换让矩阵 $\varepsilon_{ij}$ 对角化（不改变密度），得到： $$\hat{h}_{KS} \phi_i = \varepsilon_i \phi_i$$

拉格朗日乘子 $\varepsilon_i$ 自然变成特征值——这是变分带约束问题的一般结果。

在基组下变成矩阵特征值¶

离散化（展开到平面波基组 $\phi_i = \sum_k c_{ik} \chi_k$）后：

\[\sum_l H_{kl} c_{il} = \varepsilon_i \sum_l S_{kl} c_{il}\]

矩阵形式： $$\mathbf{H} \vec{c}_i = \varepsilon_i \mathbf{S} \vec{c}_i$$

对平面波基组 $\mathbf{S} = \mathbf{I}$，变成标准特征值问题： $$\mathbf{H} \vec{c}_i = \varepsilon_i \vec{c}_i$$

这就是 QE 每次 SCF 迭代里做的事——用 Davidson / CG / PPCG 算法对角化，求前 $N_{nbnd}$ 个最小特征值和对应特征向量。

物理意义¶

$\phi_i^{KS}$ = 虚构非相互作用电子系统的单粒子态（"单电子轨道"）
$\varepsilon_i$ = 该状态的轨道能量
填充顺序：从最低 $\varepsilon_i$ 开始，填 N 个电子（Aufbau 原理）
填满的态 $\rho = \sum_{i \leq N} |\phi_i^{KS}|^2$ 给出基态密度

ML 类比¶

KS 方程在基组下 = PCA 求协方差矩阵特征向量的结构完全同构：

KS-DFT	PCA
哈密顿矩阵 $\mathbf{H}$	协方差矩阵 $\mathbf{\Sigma}$
KS 轨道（特征向量）$\vec{c}_i$	主成分向量 $\vec{v}_i$
轨道能 $\varepsilon_i$	方差（特征值）$\lambda_i$
占据低能 $N$ 条	保留最大 $N$ 个主成分

只是 KS 求"最小"特征值（能量低者填充），PCA 求"最大"特征值（方差大者保留）。底层数学是同一件事。

为什么不是一般的 PDE / ODE 而是特征值问题¶

因为 KS 算子 $\hat{h}_{KS}$ 是厄米（Hermitian）算符——任何厄米算符的变分带正交约束问题都归结到特征值方程。这是泛函分析 / 算子理论的标准结论。

Q5：有效势 $V_{eff}$ 是什么？为什么突然出现？符号似乎与 $E$ 项不重叠¶

原文："$V_{eff}(\vec{r}) = V_{ext}(\vec{r}) + V_H(\vec{r}) + V_{xc}(\vec{r})$"

$V_{eff}$ 是什么¶

KS 系统是虚构的非相互作用电子系统——电子之间没有直接相互作用，每个电子只感受一个统一的有效势 $V_{eff}(\vec{r})$。

$V_{eff}$ 就是这个统一势场：它把真实系统中"电子互相推搡"的所有复杂效应，平均、打包成一个等效的外部势场，让虚构电子"像自由电子一样"各自在其中运动。

为什么需要 $V_{eff}$——因为 KS 方程是单粒子方程¶

真实多体薛定谔方程： $$\hat{H} \Psi = \left[\sum_i \hat{T}_i + \sum_i V_{ext}(\vec{r}_i) + \sum_{i<j}\frac{1}{|\vec{r}_i-\vec{r}_j|}\right] \Psi = E \Psi$$

有相互作用项 $\sum_{i<j} 1/|\vec{r}_i - \vec{r}_j|$——每个电子感受的"势场"依赖其他 N-1 个电子的位置。不是单粒子问题。

KS 方程： $$\left[-\tfrac{1}{2}\nabla^2 + V_{eff}(\vec{r})\right]\phi_i = \varepsilon_i \phi_i$$

只有单粒子算符。每个电子只感受 $V_{eff}$，不看其他电子。是 N 个独立的单粒子问题。

→ $V_{eff}$ 的作用就是把"多体"的信息编码成"单粒子势场"。这是 KS 妙招的关键。

为什么 $V_{eff} = V_{ext} + V_H + V_{xc}$¶

从能量泛函对密度的泛函导数推出：

\[\frac{\delta E[\rho]}{\delta \rho(\vec{r})} = \mu \quad \text{（化学势约束）}\]

展开 $E[\rho] = T_s + V_{ext}[\rho] + J[\rho] + E_{xc}[\rho]$：

\[\frac{\delta T_s}{\delta \rho(\vec{r})} + V_{ext}(\vec{r}) + \frac{\delta J}{\delta \rho(\vec{r})} + \frac{\delta E_{xc}}{\delta \rho(\vec{r})} = \mu\]

后三项合起来就是 $V_{eff}(\vec{r})$： - $\delta V_{ext}[\rho] / \delta \rho = V_{ext}(\vec{r})$（因为 $V_{ext}[\rho] = \int V_{ext} \rho$） - $\delta J[\rho] / \delta \rho = V_H(\vec{r}) = \int \rho(\vec{r}')/|\vec{r}-\vec{r}'| d^3r'$ - $\delta E_{xc}[\rho] / \delta \rho = V_{xc}(\vec{r})$

所以 $V_{eff}$ 是 $E[\rho]$ 对 $\rho$ 求变分后的自然结果——不是无中生有，而是数学推导的必然。

"符号似乎不重叠"的问题¶

你可能觉得：$E[\rho]$ 里有 $E_{xc}$，KS 方程里却是 $V_{xc}$ 且同时有 $V_H$ 和 $V_{ext}$——感觉符号错位。

澄清：$V_{xc}$ 和 $E_{xc}$ 不是同一量，它们是"能量"和"能量对密度的导数"的关系：

符号	类型	关系
$E_{xc}[\rho]$	能量泛函（输入密度，输出数）	$E_{xc} = \int \epsilon_{xc}(\vec{r}) \rho(\vec{r}) d^3r$
$V_{xc}(\vec{r})$	势场函数（输入位置，输出数）	$V_{xc} = \frac{\delta E_{xc}}{\delta \rho(\vec{r})}$

类比： - $E(x) = x^2$ vs $\frac{dE}{dx} = 2x$ - 同一物理实体的"能量形式"vs"其导数（驱动力/势场）"

能量泛函的各项都有对应的"势"： | 能量项 | 对应的势 | |--------|---------| | $V_{ext}[\rho]$ | $V_{ext}(\vec{r})$ （本身就是势，外势）| | $J[\rho]$ | $V_H(\vec{r})$（Hartree 势）| | $E_{xc}[\rho]$ | $V_{xc}(\vec{r})$（XC 势）|

$T_s$ 不出现在 $V_{eff}$ 里，因为它通过 KS 轨道的 $-\tfrac{1}{2}\nabla^2$ 作用直接给出（不需要作为势场）。

物理图像¶

想象一个 Si 原子： - $V_{ext}$：两个 Si 核的吸引力场（两个负 $-Z/r$ 漏斗） - $V_H$：电子云产生的经典排斥场（平滑的正值分布） - $V_{xc}$：量子效应的"修正场"（负值为主，反映交换-关联洞） - $V_{eff}$ = 三者叠加：每个虚构电子在这个"整体势场"里独立运动，它们的密度和加起来等于真实 Si 的密度

Q6：KS 轨道没有物理意义？为什么叫"轨道"¶

原文："KS 轨道和真实电子波函数没有直接物理对应，它们只是'用来造出正确密度的辅助数学对象'"

这是个微妙但重要的点。答案是：KS 轨道的物理意义是"软"的——严格说是辅助构造，但实际用起来效果接近真实物理。

严格意义：只是数学辅助对象¶

真实 N 电子系统的基态是多体波函数 $\Psi(\vec{r}_1, ..., \vec{r}_N)$，不能分解成 N 个独立单电子函数
KS 轨道 $\phi_i^{KS}$ 属于一个虚构的非相互作用系统
它们唯一的"任务"是：$\sum_i |\phi_i^{KS}|^2 = \rho_{\text{真实}}$
单个 $\phi_i^{KS}$ 不对应真实电子态，$\varepsilon_i$ 也不是真实电子能级

软物理意义：实际效果很好¶

虽然严格说是"辅助"，但 KS 轨道在实际应用里：

画能带图：$\varepsilon_n(\vec{k})$ 的能带结构和实验 + 高精度方法定性一致（除带隙系统性低估）
化学键分析：KS 轨道的空间分布和节点结构反映键合/反键（用于 HOMO、d 带中心理论、Bader 电荷分析）
HOMO 严格对应负电离能（Janak 定理）——这点 HF 轨道也有（Koopmans 定理），但 KS 还带一些关联
工业上无人怀疑：QE 输出的 band structure、DOS、PDOS 全是基于 KS 轨道，被材料学界广泛使用

为什么叫"轨道"¶

有三个原因：

历史传承：Hartree 方程、Hartree-Fock 方程的单粒子函数都叫"轨道"（orbital），KS 延续了习惯
数学形式相同：$\hat{h}_{eff} \phi = \varepsilon \phi$ 的形式和原子轨道 $\hat{h}_{atom} \phi_{1s} = E_{1s} \phi_{1s}$ 完全一样
直觉类比：化学家/物理学家习惯从"轨道"角度理解电子结构（哪些原子轨道组合成分子/能带）

KS 轨道 vs HF 轨道的对比¶

	HF 轨道	KS 轨道
来源	单 Slater 行列式对真实波函数的近似	虚构非相互作用系统
物理意义	更"硬"（有 Koopmans 定理）	更"软"（只有 HOMO 严格）
精度	不含关联，能量差	含关联（经 $V_{xc}$），能量准
能带图定性	类似	类似

实际经验：KS 轨道在"画能带 / 看轨道分布"这些定性任务上，结果通常比 HF 更接近"真实"。

ML 类比¶

KS 轨道很像神经网络的 latent representation / hidden features：

KS 轨道	NN latent representation
没有固有物理意义	没有固有语义
组合后给正确密度	组合后给正确预测
用于分析定性（能带、化学键）很有用	用于可视化、特征重要性分析很有用
设计上是"辅助"	设计上是"中间变量"

你可以把 KS 轨道理解为"一组有数学结构的隐变量，恰好对应化学/物理直觉中的'电子轨道'"。

小结¶

严格数学角度：KS 轨道是辅助构造，不对应真实单电子态
实际应用角度：画能带、分析化学键、算 DOS 都用 KS 轨道，效果好
名字 "轨道"：历史 + 形式同构 + 直觉保留——大家都这么叫

所以不是"没有物理意义"，而是物理意义比 HF 更宽松但更实用。

Q7：化简 4（KS）有 LLM / agent 需要处理的部分吗¶

原文："一次 SCF 迭代 = 解一次 KS 方程 + 更新密度 + 混合"

有，但比前几步少——大部分是和化简 2、5、11 交叉的参数。

KS 本身对应的 QE 参数（LLM 基本不调）¶

参数	作用	LLM 是否需要调
`diagonalization`	对角化算法（`david`/`cg`/`ppcg`/`paro`）	❌ 基本默认 `david` 够用
`diago_full_acc`	强制所有空带也精确对角化	❌ 默认 `.false.` 够用（为速度）
`diago_david_ndim`	Davidson 子空间大小	❌ 默认 4 够用；收敛慢时调大
`diago_thr_init`	初始对角化阈值	❌ 默认就好

这些是"算法细节"，QE 默认值对绝大多数体系都合适。LLM 除非出了问题，否则不用碰。

和其他化简交叉的参数（这是 LLM 真正的任务）¶

参数	出处	LLM 任务
`nbnd`	化简 4 直接（KS 轨道数）	带隙/DOS 任务需加空带 buffer；默认 $N/2$（绝缘体）或 auto（金属）
`nspin` + `starting_magnetization`	化简 2	磁性体系必设
`occupations` + `smearing` + `degauss`	化简 9（金属 BZ）+ 化简 11（SCF）	金属必加 smearing
`input_dft` / 赝势选择	化简 5、6	选 LDA/PBE/PBEsol
`mixing_beta` + `electron_maxstep`	化简 11（SCF）	难收敛体系调整

所以 KS 方程"怎么解"的具体行为，由这一族交叉参数共同决定。

KS 层面 LLM 的元能力¶

识别"需要多少条 KS 轨道 nbnd"
绝缘体：$N_{val}/2$ 默认够
金属：需要多几条（smearing）
带隙/DOS：需要加 buffer（典型 +10–20 空带）
磁性：up/down 各半，算 $N_{val}$
错了：结果不完整（算不出带隙 / DOS 缺段 / 磁矩错）
识别"对角化算法是否合适"
默认 Davidson 对大多数问题够
超大体系（>1000 原子）可能需要 PPCG 或换 CG
金属 + Davidson 偶尔震荡 → 换 CG
识别"SCF 收敛是否受 KS 细节影响"
diago_full_acc=.true. 在某些难收敛体系有帮助
mixing_ndim 大一点（Broyden 历史）加速收敛
后处理任务的轨道需求
band_gap 任务：需要 nscf + 更密 k + 更多 nbnd
dos 任务：nscf + 非常密 k
LLM 要知道"scf → nscf 的链式调用"

对 benchmark 的直接影响¶

TritonDFT benchmark 的任务：

task	KS 层面要求
`vc_relax`（多数任务）	默认 nbnd 够，不用调 KS 细节
`scf_energy`	同上
`band_gap`	nbnd 加 buffer（+10 空带），必需！
`dos`	nbnd 加更大 buffer（+20–50），k 点更密

所以 LLM 在 band_gap 和 dos 任务上，如果不加 nbnd 空带 buffer，算不出带隙（只看到占据态）——这是 benchmark 的一个考查点。

与其他化简的任务量对比¶

化简	LLM 任务量
1 BO	低（只给核位置）
2 Slater	高（nspin/magnetization/occupations）
3 HK	间接（元能力）
4 KS	中低（主要是 nbnd + 元判断）
5 XC	高（选泛函）
6 赝势	高（选赝势）
8 平面波	高（ecutwfc）
9 k 点	高（k-mesh）
11 SCF	中（mixing_beta）
12 BFGS	中（conv_thr）

KS 的主要 LLM 任务不在 KS 算法本身，而在"配合化简 2 + 5 + 11 的交叉参数"上。单独 KS 的直接调参基本交给 QE 默认。

本轮 7 个问题速览¶

#	主题	核心
Q1	$E[\rho]$ 每项含义与量级	$T$ / $V_{ext}$ / $J$ / $E_{xc}$；Si 原胞分别 +100 / −300 / +60 / −90 eV，大量抵消成 −230 eV（决定 DFT 精度要求）
Q2	Thomas-Fermi 为什么不行	局域近似太粗；无壳层；Teller 定理证明 TF 不能束缚分子；~10% 动能误差远超化学精度；KS 用轨道算动能避开这个问题
Q3	$E_{xc}$ 组成和量级	交换 $E_x$（~10–30%）+ 关联 $E_c$（~1–5%）+ 动能修正；He 原子 $E_x = -28$ eV, $E_c = -1.1$ eV；化学性质由 $E_c$ 决定
Q4	为什么 KS 是特征值等式	变分 + 正交约束的 Lagrange 乘子就是本征值；基组下变成矩阵特征值 $\mathbf{H}\vec{c} = \varepsilon\vec{c}$，和 PCA 同构
Q5	$V_{eff}$ 是什么 / 符号错位	把多体相互作用打包成单粒子势场；$V_{eff} = V_{ext} + V_H + V_{xc}$；$V_{xc} = \delta E_{xc}/\delta\rho$（能量导数 vs 能量本身）
Q6	KS 轨道的物理意义	严格是辅助对象，实际用于能带/DOS/化学键分析效果好；类比 NN 的 latent representation
Q7	LLM 任务	KS 本身算法细节基本默认；真正任务是 nbnd buffer（带隙/DOS）+ 和化简 2/5/11 交叉参数

KS-DFT	PCA
哈密顿矩阵 \(\mathbf{H}\)	协方差矩阵 \(\mathbf{\Sigma}\)
KS 轨道（特征向量）\(\vec{c}_i\)	主成分向量 \(\vec{v}_i\)
轨道能 \(\varepsilon_i\)	方差（特征值）\(\lambda_i\)
占据低能 \(N\) 条	保留最大 \(N\) 个主成分

14-kohn-sham 问答笔记¶

Q1：能量泛函 \(E[\rho] = T + V_{ext} + J + E_{xc}\) 每一项的物理含义¶

Q2：为什么 Thomas-Fermi 误差巨大¶

为什么这个假设对真实原子 / 分子 / 固体误差巨大¶

根本问题¶

Kohn-Sham 的妙招¶

Q3：\(E_{xc}\) 是什么？和其他项的量级关系¶

\(E_{xc}\) 的组成¶

量级对比表（He 原子最清晰）¶

几个重要观察¶

在 DFT 实际计算的角色¶

对 benchmark 的影响¶

Q4：为什么 KS 方程写成特征值等式的形式¶

从变分到特征值¶

在基组下变成矩阵特征值¶

物理意义¶

ML 类比¶

为什么不是一般的 PDE / ODE 而是特征值问题¶

Q5：有效势 \(V_{eff}\) 是什么？为什么突然出现？符号似乎与 \(E\) 项不重叠¶

\(V_{eff}\) 是什么¶

为什么需要 \(V_{eff}\)——因为 KS 方程是单粒子方程¶

为什么 \(V_{eff} = V_{ext} + V_H + V_{xc}\)¶

"符号似乎不重叠"的问题¶

物理图像¶

Q6：KS 轨道没有物理意义？为什么叫"轨道"¶

严格意义：只是数学辅助对象¶

软物理意义：实际效果很好¶

为什么叫"轨道"¶

KS 轨道 vs HF 轨道的对比¶

ML 类比¶

小结¶

Q7：化简 4（KS）有 LLM / agent 需要处理的部分吗¶

KS 本身对应的 QE 参数（LLM 基本不调）¶

和其他化简交叉的参数（这是 LLM 真正的任务）¶

KS 层面 LLM 的元能力¶

对 benchmark 的直接影响¶

与其他化简的任务量对比¶

本轮 7 个问题速览¶

项	数学形式	物理含义	单位/量级（Si 原胞参考）
\(T[\rho]\) 动能泛函	难写成 \(\rho\) 的显式形式；KS 中替换成 \(T_s = \sum_i \int \phi_i^* (-\tfrac{1}{2}\nabla^2) \phi_i\)	N 个电子运动总动能的期望值（量子力学意义上）	+100 eV/原胞（正值，"动"总是正能量）
\(V_{ext}[\rho]\) 外势能	\(\int V_{ext}(\vec{r}) \rho(\vec{r}) d^3r\)（简单积分）	电子云与原子核之间的静电吸引能——电子被核"吸引"获得能量	−300 eV/原胞（负值，吸引性主导）
\(J[\rho]\) 经典库仑能 / Hartree 能	\(\frac{1}{2} \iint \frac{\rho(\vec{r})\rho(\vec{r}')}{\\|\vec{r}-\vec{r}'\\|} d^3r d^3r'\)	把电子密度当经典电荷分布算它的自库仑排斥（忽略了量子性）	+60 eV/原胞（正值，电子互斥）
\(E_{xc}[\rho]\) 交换-关联能	真实多体贡献与 \(T_s, J\) 等的差	两类量子效应：交换（同自旋费米子回避）+ 关联（所有电子动态躲避）+ \(T - T_s\) 动能修正	−90 eV/原胞（负值，降低能量）

能量项	值（eV）	占精确总能量
\(T\)（真实动能）	+78.87	99.8% 正
\(V_{ext}\)（核吸引）	−199.45	252% 负（主导）
\(J\)（经典库仑）	+34.03	43% 正
\(E_x\)（交换）	−27.89	35%（大！）
\(E_c\)（关联）	−1.14	1.4%（小但化学决定）
\(E_{xc}\) 总	−29.03	37%
\(E_{total}\)（精确）	−79.01	100%

符号	类型	关系
\(E_{xc}[\rho]\)	能量泛函（输入密度，输出数）	\(E_{xc} = \int \epsilon_{xc}(\vec{r}) \rho(\vec{r}) d^3r\)
\(V_{xc}(\vec{r})\)	势场函数（输入位置，输出数）	\(V_{xc} = \frac{\delta E_{xc}}{\delta \rho(\vec{r})}\)

参数	出处	LLM 任务
`nbnd`	化简 4 直接（KS 轨道数）	带隙/DOS 任务需加空带 buffer；默认 \(N/2\)（绝缘体）或 auto（金属）
`nspin` + `starting_magnetization`	化简 2	磁性体系必设
`occupations` + `smearing` + `degauss`	化简 9（金属 BZ）+ 化简 11（SCF）	金属必加 smearing
`input_dft` / 赝势选择	化简 5、6	选 LDA/PBE/PBEsol
`mixing_beta` + `electron_maxstep`	化简 11（SCF）	难收敛体系调整

#	主题	核心
Q1	\(E[\rho]\) 每项含义与量级	\(T\) / \(V_{ext}\) / \(J\) / \(E_{xc}\)；Si 原胞分别 +100 / −300 / +60 / −90 eV，大量抵消成 −230 eV（决定 DFT 精度要求）
Q2	Thomas-Fermi 为什么不行	局域近似太粗；无壳层；Teller 定理证明 TF 不能束缚分子；~10% 动能误差远超化学精度；KS 用轨道算动能避开这个问题
Q3	\(E_{xc}\) 组成和量级	交换 \(E_x\)（~10–30%）+ 关联 \(E_c\)（~1–5%）+ 动能修正；He 原子 \(E_x = -28\) eV, \(E_c = -1.1\) eV；化学性质由 \(E_c\) 决定
Q4	为什么 KS 是特征值等式	变分 + 正交约束的 Lagrange 乘子就是本征值；基组下变成矩阵特征值 \(\mathbf{H}\vec{c} = \varepsilon\vec{c}\)，和 PCA 同构
Q5	\(V_{eff}\) 是什么 / 符号错位	把多体相互作用打包成单粒子势场；\(V_{eff} = V_{ext} + V_H + V_{xc}\)；\(V_{xc} = \delta E_{xc}/\delta\rho\)（能量导数 vs 能量本身）
Q6	KS 轨道的物理意义	严格是辅助对象，实际用于能带/DOS/化学键分析效果好；类比 NN 的 latent representation
Q7	LLM 任务	KS 本身算法细节基本默认；真正任务是 nbnd buffer（带隙/DOS）+ 和化简 2/5/11 交叉参数