计算解释材料,实验检验解释

从原子到器件信号,计算要选择保留哪些物理,实验要区分哪些解释成立。

从原子到响应,每次简化都有自己的误差

第一性原理从元素、结构和基本相互作用出发计算,不靠拟合目标实验来决定答案。它允许近似。DFT(密度泛函理论)用电子密度描述基态能量。有效理论只保留目标能量与长度尺度需要的自由度,并把其余自由度的影响写进参数和测量算符。

电子激发可以用 GW(电子自能修正)处理,电子与空穴组成的激子可以用 BSE(Bethe–Salpeter 方程)处理。两者常用于修正能谱与光学响应。

原子位置 · 元素 · 电子相互作用 · 外场与环境微观规则可以写出,真实器件的结构仍需确定 DFT:求密度、总能量与辅助轨道猜密度求有效势求轨道更新密度,直到自洽物理近似:交换与相关的泛函 直接多体计算求多电子波函数例如量子 Monte Carlo绕开近似 DFT 泛函 结构与稳定性能量 · 力 · 密度可直接从 DFT 得到 第一性原理导出的有效理论低能轨道 · 屏蔽作用 · 算符物理近似:消去与截断 修正激发与光谱GW 自能 · BSE 激子仍有相关与顶角近似 多体求解平均场 · 涨落 · 数值波函数物理近似:求解器与有限尺寸 可测响应:光谱、磁矩、压缩率、电流…… 所有路径都要检查网格、基组和收敛电阻还要散射与边界 原子、电子、外场与环境微观规则可以写出,结构仍需确定 DFT:求密度、能量与辅助轨道猜密度 → 求有效势 → 求轨道更新密度,反复直到自洽物理近似:交换与相关的泛函 第一性原理导出的有效理论保留低能轨道、作用与测量算符物理近似:消去与截断 多体求解平均场、涨落或数值波函数物理近似:求解器与有限尺寸 可测响应光谱、磁矩、压缩率、电流……旁路:DFT → 结构;GW/BSE → 光谱也可直接求多电子波函数,不必经过 DFT
点击图中步骤,查看它在算什么。DFT 不等于无相互作用,有效理论也不等于经验拟合。 独立图:SVG · PNG

DFT 算的是密度与能量,辅助能带还不是完整的真实光谱

常用的 Kohn–Sham 算法构造一组辅助单电子轨道,使它们给出所需电子密度。电子相互作用通过有效势计入。未知的交换与相关部分通常用近似泛函处理。自洽只说明这个近似内部算一致了。

Kohn–Sham,1965

DFT 实际做了哪些运算?温度如何进入?

给定原子位置,程序先猜电子密度。它由密度算出有效势,再把微分方程离散成矩阵,求辅助轨道与占据数。新轨道给出新密度。程序反复更新,直到密度与能量变化满足阈值。随后可算原子受力,再优化结构。

概念上,总能量由辅助电子的动能、原子对电子的势能、平均静电作用、交换相关能和原子间作用组成。“泛函”就是输入整幅密度分布,输出一个数的规则。精确泛函若已知,基态 DFT 在理论上可以精确。实际泛函、赝势与基组通常有近似。

Mermin 的有限温理论把目标改为热平衡自由能。有限温本身不是禁止计算的障碍。电子热占据、原子热运动与动力学需要分别处理。为积分方便加入的“展宽”,也不总等于真实器件温度。

有效理论为什么能保留第一性原理来源?

低能电子可以虚跃迁到高能态再回来。因此,删掉高能态时必须把它们的影响写回低能描述。对双层石墨烯,这可能留下四个碳轨道、修正后的跃迁与屏蔽相互作用。算光、磁与电流响应时,还必须修正对应的测量算符。

精确消去高能态可定义为能量相关的有效 Hamiltonian(编码能量与时间演化的算符)。实际模型通常只保留少量静态参数、两体作用和有限空间范围。它是可改进的近似。不能仅因为低能能带拟合得好,就认定全部响应准确。

第一性原理导出的模型与直接拟合实验的模型都可以有用,但证据来源不同。用于比较实验准确度时,要说明哪些参数曾用实验校准。

误差如何减少?算力能减少其中哪一部分?
误差来源增加什么仅加算力是否足够
网格、基组、统计与有限尺寸更密采样、更多轨道、更大系统与外推可以改进,需观察收敛
近似交换相关与电子激发更完整的多体方法与交叉核验固定近似的偏差会保留
有效理论截断扩大空间,保留动态作用与测量算符需要扩展模型内容
平均场与波函数限制加入涨落,用更强求解器部分问题可改进,部分遇到指数代价
实际堆垛、应变、散射与接触结构表征、微观散射与器件边界未知输入不能由收敛自动识别
“低密度”和“小能标”具体有多小?

无隙、忽略三角翘曲的抛物线双层石墨烯,载流子密度 10¹⁰ 与 10¹¹ cm⁻² 对应费米能量约 0.35 与 3.5 meV(毫电子伏)。这只是量级换算。有位移场、关联作用或三角翘曲时,要重新算实际能谱。

三角翘曲会把低能费米线拆成多个口袋,其参考能标约 1 meV。1 K 的热能是 0.086 meV。因此 0.1 meV 相对 eV(电子伏)级原子能量很小,却相当于上述 0.35 meV 的约 29%,足以影响低温响应。相排序是否翻转,要与候选相自由能差比较。

小绝对能量本身可以计算。第一性原理文献算出过约 24 μeV(微电子伏)的本征自旋轨道劈裂。这个数值不等于对现实器件已认证同样小的总误差。

McCann–Koshino 综述,式 95 · Konschuh 等,2012

凝聚态计算研究材料的状态、激发和响应

共同对象是原子或格点、相互作用、外场与温度。研究者求量子态或热平衡分布,再计算能量、关联和测量信号。研究方向按问题划分,方法跨方向使用。DFT、量子 Monte Carlo 和张量网络不是互斥的研究领域。

材料结构能量 · 力 · 声子稳定性、缺陷与材料设计 电子与光能带 · 寿命 · 激子光谱、磁响应与拓扑 强关联与新相关联 · 序参量 · 自由能磁性、超导与量子 Hall 共同计算对象原子/格点 + 相互作用 + 外场 + 温度量子态或热分布 → 能量、关联、激发与响应 输运与器件电流 · 散射 · 空间势电阻、噪声与非线性响应 温度与动力学热分布 · 原子运动 · 演化相变、热输运与泵浦响应 计算方法泛函 · 算法 · 数据模型扩大可算范围、控制误差 同一材料可以同时属于多个方向 常见资源跨度 · 规划量级 工作站小型有效模型 一个或数个节点常规电子结构 小集群至超算精细多体、大材料与扫描 具体实例见下表。单点计算与完整研究的资源不同。 共同计算对象原子/格点、相互作用、外场、温度求量子态或热分布再求能量、关联、激发与响应 材料结构能量、力、声子稳定性与材料设计 电子与光能带、寿命、激子光谱与拓扑响应 强关联与新相关联、序参量磁性、超导、Hall 输运与器件电流、散射、空间势电阻与噪声 温度与动力学热分布、运动、演化相变与泵浦响应 计算方法泛函、算法、数据规模与误差控制 方法跨方向使用同一材料可以同时属于多个方向 常见资源跨度 · 规划量级工作站:小型有效模型一个或数个节点:常规电子结构小集群至超算:精细多体、大材料与扫描
点击方向,查看典型问题与方法。声子是晶格振动的量子,序参量是区分有序状态的量。 独立图:SVG · PNG

材料只是输入,真正要计算的是它如何组织与响应

常见输出包括总能量、电子密度、相关函数、激发谱与响应函数。相关函数描述两个位置或时刻之间的联系。响应函数描述外加电场、磁场或光改变后,电流、磁化或吸收如何变化。

通常用多少算力?先分清一次任务与完整研究

简单能带与有效模型可在工作站上运行。复杂电子结构常用一个或多个计算节点。高精度多体波函数和大型材料计算会使用数十到数千张 GPU(图形处理器)。CPU(中央处理器)也用于这些任务。领域没有统一的平均成本。下列是公开任务的实测点,不能按电子数直接横向比较。V100 与 A100 是两代 GPU 型号。

公开算例机器与耗时任务边界
钨原胞,密集动量网格32 CPU 核 × 134 分钟单个振动波矢响应与前置电子态
431 原子、6034 电子24 V100 GPU × 335 秒有限元 DFT 的一次基态计算
8191 原子、114674 电子3600 V100 GPU × 3250 秒大型有限元 DFT 基态计算
石墨烯,8 个碳原子、48 全电子32 A100 GPU × 3 天直接多体神经波函数的一次训练

CPU 核时 = CPU 核数 × 小时。GPU 时 = GPU 数 × 小时。不同硬件、方法与精度不能据此直接比较。结构优化、尺寸外推、外场扫描和重复验证需另计。

有限元 DFT 大规模计算 · DeepSolid 直接波函数计算 · 算力原始来源与限定

为什么小有效模型有时比数千电子 DFT 更难?

DFT 主要处理辅助单电子轨道。直接多体求解要描述电子配置之间的量子关联。后者的状态空间可随轨道数指数增长。一个小模型的所有多电子配置,仍可能比大材料的近似单电子计算更昂贵。

例如固定半填充、固定自旋分量的 Hubbard 模型,16、20、24 个格点,仅保存一个双精度复数状态向量就约需 2.65 GB、546 GB、117 TB。实际求解还需要多个向量与矩阵运算。这是精确对角化的规模示例,不是所有算法都必须支付的代价。

有符号问题的 Monte Carlo,采样代价还会受到正负贡献抵消影响。张量网络与神经波函数能避开部分难点,但精度仍要用尺寸、表示能力和方法交叉检查。

测到新信号以后,还要判别内部发生了什么

实验学家制备、表征和控制样品,测电流、电压、吸收、磁化或局部信号。器件把内部量子状态转成可读取的响应。理论与计算把候选机制转成可检验预测。一条曲线通常不能唯一确定一个相。实验学家也会建模,理论学家也会计算,角色彼此重叠。

制备与表征样品、栅极、接触、结构知道实际测的是什么 改变实验条件电场、磁场、温度、光控制系统受到的扰动 读出信号电流、电压、吸收、电容测量有误差与器件影响 找到可重复的现象峰、平台、高阻区、临界温度、对称性变化“出现异常”还不是“知道是哪一种相” 理论与计算提出候选机制模型 → 内部状态 → 多个观测预测同一信号可以有不同来源 实验设计判别测量第二种探针、新条件、新样品检验没有参与拟合的预测 检验结果更新模型,也决定下一次实验 制备、表征并控制器件结构、栅极、接触、电场、磁场、温度 读出信号电流、电压、吸收、电容等 找到可重复的现象峰、平台、高阻区、临界温度异常信号不唯一确定内部状态 候选机制与可检验预测模型连接内部状态与多个观测理论、计算、实验共同建模 判别测量新条件、新样品、第二种探针检验没有参与拟合的预测检验结果更新解释与下一次实验
模型解释看不见的自由度,判别测量减少解释的歧义。 独立图:SVG · PNG

判别实验问的是:哪一个解释还能预测新数据?

把同一组参数用于不同探针和新条件,解释才受到更强约束。例如双层石墨烯中的高阻区,可以进一步测温度依赖、Hall 响应、压缩率、光谱或倾斜磁场。各探针对能带、关联与无序的敏感性不同。

实验学家还做什么?为什么需要理论与计算?

实验不仅测器件。还包括材料生长、结构与局部探针、光谱、散射、磁测量、超快动力学和低温高压。器件是其中一种组织样品、控制条件并读取信号的方法。

模型能说明哪些自由度重要、哪些对称性受到破坏,以及不同机制会给出什么预测。计算能把模型变成定量曲线,并检查在实际参数下是否足够大。实验能确认器件输入、筛掉机制,并发现现有模型没有解释的事实。

例如高阻区可以来自能带开隙、关联有序、无序局域化或接触。单独拟合一条电阻曲线不能唯一选择机制。额外探针和跨条件预测能增加证据,但仍需说明排除了哪些替代解释。

双层石墨烯的文献支持改进,尚未认证全面超过实验

普通 AB(Bernal 堆垛)双层石墨烯的已有工作显示:电子结构修正、有效模型和响应公式都可能带来大变化。能隙是指定激发的最低能量门槛,定义随观测量不同。谷 g 因子是磁场引起两谷能量差的响应系数。更密的数值网格不能代替更完整的物理描述。

加了什么计算怎样变化证据能支持什么
GW 修正电子自能准粒子能隙比 DFT 增大 56–81%电子结构近似可影响几十 meV
有效模型加入层内极化修复低掺杂能隙约两倍的差别降阶可以改进,旧简单模型漏了物理
补齐激子磁响应p 激子谷 g 因子 13.03 → 1.79,实验 1.4±0.8响应公式漏项可造成数量级偏差
高精度有效模型求解量子 Hall 激活隙 5.6 K,实验 5.1±0.2 K定量成功,仍不是完整第一性原理认证
调研中,哪些判断被具体证据修正?

2026 激子结果是预印本。其两种假定 hBN(六方氮化硼封装层)堆垛给 s 激子谷 g 因子 13.48 和 20.32。更接近实验的结构不是对真实器件结构的独立确认。BSE 激子能量收敛至 5 meV,也不能当成 g 因子的总误差。

原始文献与完整调研证据