第一性原理计算(first-principles calculation)是基于密度泛函理论(Density Functional Theory,简称DFT)的一种量子力学计算方法,它是建立在波恩-奥本海默绝热近似、非相对论近似以及单电子近似的基础上的计算理论。第一性原理计算的基本思想是将多个原子构成的体系看成是由多个电子和原子核组成的系统,只需要最基本物理常数:m0、e、h、c、kB,不采用任何实验得到的数据,以及经验、半经验参数来求解薛定愕方程。第一性原理计算能够从微观原子层面上揭示材料性质的本质规律,具有较高的预测精度,广泛应用于材料科学、化学、物理、催化、生命科学等领域。
第一性原理常见的体系包括但不限于金属单质、金属氧化物、合金、半导体材料、二维材料、钙钛矿、MOF、COF、表面结构、异质结构等。
常用的软件有:VASP, CASTEP, CP2K, QE等。
可以的计算的内容包括但不限于:
结构性质:晶格常数、原子坐标、键长、键角、XRD、表面能、界面能、形成能
电子性质:电荷密度、态密度、能带,投影态密度、投影能带、COHP、静电势、ELF、功函数
力学性质:弹性常数、弹性模量、体积模量、德拜温度、原子振动频率、压电张量
吸附体系:吸附能、结合能、差分电荷密度、bader电荷、mulliken电荷
催化性质:OER、ORR、HER、CO2RR、NRR、NO3RR、NORR、SRR、迁移能垒、过渡态搜索
光学性质:介电函数、折射率、吸收光谱、反射光谱、能量损失函数、激子结合能
声子性质:声子谱、声子态密度
磁 性 质:磁矩、磁化率
自由能台阶图
迁移能垒
能带
态密度
吸附能
过渡态
差分电荷密度
功函数
声子谱
光学性质
弹性性质
磁性质
电荷密度
空位形成能
火山图
COOP/COHP
声子态密度
电荷布居
电子局域密度函数
内建电场
布拜图
相图
d带中心
静电势
自旋密度
界面能
介电常数
剪切模量
能带电位匹配
泊松比
表面能
第一性原理分子动力学(AIMD)
自由能台阶图用于描述催化反应各个基元步骤在反应过程中的能量变化,是研究反应机理和判断速控步骤的重要工具。通过DFT计算不同中间体的吸附自由能,并结合零点能、熵和温度修正,可以构建完整的反应自由能路径。
该分析方法常用于HER、OER、ORR、CO2还原、NRR等催化体系研究,可帮助判断哪一步最困难、理论过电位有多大、不同催化剂活性差异来自哪里。对于催化机理分析和高水平SCI论文写作,自由能台阶图是非常经典且有说服力的结果之一。
迁移能垒用于表征原子、离子或分子从一个位置移动到另一个位置时所需克服的能量障碍,是研究扩散能力和动力学性能的重要参数。通常结合NEB等方法计算得到,可以反映迁移路径是否通畅、传输是否容易发生。
该性质广泛应用于锂离子电池、钠离子电池、固态电解质、气体扩散、表面迁移和缺陷迁移研究。通过迁移能垒分析,可以判断离子传导快慢、扩散限制步骤以及不同掺杂或缺陷对迁移行为的影响,是研究储能材料和动力学过程的重要计算内容。
体现结构的电子性质的常用方法之一。
能带结构是指材料中所有电子能级的分布情况,包括价带和导带,价带顶和导带底之间的禁带宽度即为带隙。其中,价带顶的能量位置决定了材料的化学反应性质,导带底的能量位置决定了材料的电子传输性质。
态密度反映材料在不同能量范围内电子态的分布情况,是分析电子结构的核心结果之一。总态密度可展示整体电子态分布,分波态密度则能进一步区分不同原子和轨道对电子结构的贡献。
态密度常用于研究导电性、掺杂效应、缺陷调控、催化活性来源以及成键特征等问题。通过观察费米能级附近的态密度变化,可以判断材料电子活性强弱;通过PDOS还可分析反应活性位点和成键特征。它是能带结构的重要补充,也是电子结构论文中的基础图件。
吸附能用于衡量分子、原子或反应中间体在材料表面吸附的强弱,是研究表面反应、催化活性和界面相互作用的基础指标。通过比较吸附前后体系总能量的变化,可比较不同吸附构型或位点的相对稳定性,并分析不同位点、不同构型的吸附优劣。
吸附能广泛用于催化剂筛选、气体传感、储能材料、污染物去除等研究中。它能够帮助回答“活性物种是否容易吸附”“哪个位点更稳定”等关键问题,是表面科学和催化研究中最常见的计算项目之一。
过渡态是基元反应路径中能量最高、最不稳定的关键结构,代表体系从反应物转化为产物时必须跨越的临界状态。通过DFT过渡态搜索,可以得到反应过程中的关键构型,并用于构建完整的反应机理。
过渡态计算常用于表面催化、分子反应、解离过程、成键断键过程和扩散过程研究。它能够帮助判断反应是否容易发生、哪条反应路径更优、不同催化剂对反应活化的促进效果如何。对于机理类SCI论文而言,过渡态结果能显著增强结论的深度和说服力。
差分电荷密度用于展示体系相互作用前后电子密度的重新分布情况,可直观反映电子积累和耗尽区域,是分析界面作用、吸附成键和电荷转移趋势的重要方法。通常通过吸附体系电荷密度减去各独立组分电荷密度得到。
该分析常用于分子吸附、异质结界面、电催化中间体结合和掺杂调控研究。通过差分电荷密度图,可以定性观察电子重新分布方向;若需进一步定量分析电荷转移,通常可结合Bader、Hirshfeld等电荷分析方法。
其大小标志着束缚电子的强弱,功函数越小,电子越容易逸出材料表面。
声子谱用于描述材料中晶格振动模式及其频率分布,是判断结构动力学稳定性和研究热学性质的重要工具。通过计算声子色散关系和声子态密度,可以分析材料在不同波矢下的振动行为,并判断是否存在虚频。
声子谱常用于晶体材料、二维材料和新型功能材料研究。若体系不存在虚频,通常说明结构在动力学上稳定;若出现明显虚频,则可能意味着结构不稳定或存在相变趋势。此外,声子谱还能为热容、熵、自由能和热导率等性质分析提供基础支持。
光学性质用于表征材料与光相互作用时的响应行为,包括吸收、反射、折射、介电函数、吸收系数和能量损失函数等。通过DFT及相关方法,可以预测材料在不同波长范围内的光学响应特征。
该性质常用于光催化、光电器件、太阳能材料、发光材料和透明导电材料研究。通过光学性质分析,可以判断材料对可见光或紫外光的吸收能力、是否适合用于光响应应用,以及掺杂或缺陷是否能够改善光吸收范围。对于研究材料光电性能的论文,光学计算是十分常见的重要内容。
弹性性质反映材料在外力作用下抵抗形变并恢复原状的能力,是评估材料力学性能的重要指标。通过DFT计算弹性常数、体模量、剪切模量、杨氏模量和泊松比等参数,可以系统分析材料的刚性、抗压/抗剪能力、延展性趋势和各向异性。
该性质广泛用于结构材料、陶瓷材料、二维材料、能源材料和功能薄膜研究。通过弹性性质分析,可以判断材料是否容易变形、是否具有机械稳定性、适不适合作为器件基底或工程材料。对于兼顾性能与稳定性的SCI研究,弹性参数是非常重要的理论支撑结果。
体现结构的磁各向异性、磁导率等
磁性质是指物质在磁场中表现出的一系列性质,包括磁化、磁导率、磁滞、磁各向异性等。物质的磁性质与其电子结构密切相关,尤其是未成对的电子自旋。通过第一性原理计算,可以对结构的磁序、自旋密度、磁矩等进行理论计算分析。
空位形成能用于衡量材料中形成原子空位缺陷所需的能量,是研究缺陷稳定性和缺陷浓度的重要指标。形成能越低,说明该类缺陷越容易形成,在实际材料中出现的可能性也越高。
该性质广泛用于氧空位、金属空位、阴离子空位等缺陷工程研究,尤其在催化、储能、半导体和氧化物材料领域应用非常广泛。通过空位形成能计算,可以分析哪种缺陷更稳定、缺陷对电子结构和催化活性的影响,以及材料在不同环境下是否容易产生缺陷,是研究缺陷调控机制的重要工具。
理解催化活性与结构性质之间的关系。
火山图是一种用来描述催化反应速率与反应物浓度之间关系的图表。它基于阿伦尼乌斯方程(Arrhenius equation),该方程描述了化学反应速率与温度和反应物活化能之间的关系。火山图的名称来源于其图表形状,类似于火山的轮廓,
HOMO-LUMO分析是研究分子前线轨道性质的经典方法,其中HOMO代表最高占据分子轨道,LUMO代表最低未占据分子轨道。它们反映了分子最容易失去电子和接受电子的区域,其能级差可作为分子稳定性、电子跃迁难易和反应活性的参考指标。
该性质常用于分析分子稳定性、电荷转移能力、光电响应特征及反应活性。通过HOMO和LUMO轨道分布,可以判断反应可能发生的位置;通过HOMO-LUMO能隙,则可比较不同分子的活性强弱和激发难易程度。对于有机分子、配合物和光功能分子研究,这是一项非常基础且常用的计算内容。
电子局域函数用于表征电子在空间中的局域化程度,是分析化学键类型和电子分布特征的重要工具。通过ELF图,可以辅助分析共价性、电子局域化程度、孤对电子和金属键特征。
该性质常用于研究成键机制、界面相互作用、缺陷调控和催化活性中心特征。相比仅看键长或电荷转移,ELF能够更清楚地展示电子是否集中在某一化学键附近,或者是否存在明显的局域电子区域,因此在解释复杂材料中成键本质方面非常有帮助,也是高水平理论研究中常见的图像结果之一。
影响电子和空穴在异质结中的输运。
异质结内建电场是指在两种不同半导体材料接触时,由于能带差异而在界面附近形成的电场。这种电场是由于两种材料中的电子和空穴在热力学平衡状态下的化学势不同而产生的。
判断在特定条件下哪种结构是稳定的。
Pourbaix图,也称为电位-pH图或布拜图,是一种用于表示水溶液中不同氧化态的化学物质稳定性的图表。它是由比利时化学家Maurice Pourbaix于1940年代提出的,因此以他的名字命名。Pourbaix图在电化学、腐蚀科学、矿物学和环境化学等领域有着广泛的应用。
了解材料在不同条件下的相变行为。
结构的相图是材料科学中用来描述材料在不同条件(如温度、压力)下相变(固相、液相、气相等)的图形表示。它通常包括一个或多个坐标轴,表示温度和/或压力,以及相变区域、相边界和相平衡点。
d带中心是描述过渡金属材料d电子态相对费米能级位置的重要参数,常用于分析金属表面对吸附物的结合能力和催化活性趋势。在许多过渡金属催化体系中,d带中心越接近费米能级,表面与反应中间体的相互作用往往越强。
该分析常用于电催化、热催化、合金催化和单原子催化研究。通过d带中心可以解释掺杂、应变、缺陷或组分调控对吸附强度和催化活性的影响。但该规律属于经验描述符,实际分析中通常需要结合吸附能、PDOS和反应自由能共同讨论。
静电势用于描述分子或孤立体系周围空间中电荷分布所产生的电势变化,是分析反应活性位点和分子表面电性特征的重要工具。通过计算静电势分布,可以直观识别体系中电子富集区和电子贫化区,从而判断哪些位置更容易受到亲核或亲电试剂进攻。
该性质常用于催化活性位分析、分子识别、主客体相互作用、药物分子反应位点预测等研究。通过ESP图或表面映射结果,可以帮助解释为什么某些原子更容易参与反应、不同取代基如何改变分子反应性,是分子体系机理分析中非常常见的一类电子结构结果。
界面能用于描述两种材料在接触形成界面时的能量代价,是研究异质结、复合材料和多相界面稳定性的重要指标。界面能越低,通常说明界面结合越稳定,更容易形成高质量的界面结构。
该性质常用于异质结催化、二维材料堆叠、金属/半导体接触、电池界面和复合材料研究。通过界面能计算,可以判断不同界面构型的稳定性、比较不同取向或堆垛方式的优劣,并解释界面结合强弱与性能之间的关系。对于研究界面调控、载流子分离和复合材料稳定性的SCI论文,界面能具有重要参考价值。
材料在受到压缩或拉伸时的变形特性。
泊松比(Poisson’s Ratio)是一个描述材料在受到压缩或拉伸时横向(垂直于受力方向)形变程度的无量纲物理量。它是材料力学性质的一个重要参数,用于表征材料的弹性行为。泊松比定义为横向形变与轴向形变之比,即横向形变与轴向形变之比。
表面能用于表征形成某一晶体表面所需的能量,是判断不同晶面稳定性的重要参数。一般来说,表面能越低,说明该表面越稳定、越容易在实际晶体中暴露出来。
表面能广泛应用于晶体生长、纳米材料形貌控制、催化剂暴露晶面筛选和界面工程研究。通过表面能计算,可以分析不同晶面的稳定顺序,解释实验中观察到的颗粒形貌,并帮助选择更有可能具有高活性的暴露面。对于表面催化和晶面效应相关SCI论文,表面能是非常基础而关键的理论指标。