2. 密度泛函理论  Hohenberg-Kohn 定理 1 )基态系统的所有物理性质都由电子密度唯一决定,能量与电子密度为一一映射。 2 )对应于电子密度的变分原理:任意近似电子密度所对应的能量值都大于等于基 态对应的真正密度所决定的能量值。  密度泛函理论 ( Density Functional.

Slides:



Advertisements
Similar presentations
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
Advertisements

第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
1 CASPT2//CASSCF and TDDFT//CASSCF Mapping of the Excited State Isomerization Path of a Minimal Model of the Retinal Chromophore Simona Fantacci, † Annapaola.
信号与系统 第三章 傅里叶变换 东北大学 2017/2/27.
碰撞 两物体互相接触时间极短而互作用力较大
第三章 函数逼近 — 最佳平方逼近.
第五节 微积分基本公式 、变速直线运动中位置函数与速度 函数的联系 二、积分上限函数及其导数 三、牛顿—莱布尼茨公式.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
不确定度的传递与合成 间接测量结果不确定度的评估
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
1.5 场函数的高阶微分运算 1、场函数的三种基本微分运算 标量场的梯度f ,矢量场的散度F 和F 旋度简称 “三度” 运算。
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
5-4 分子性质的自洽场-分子轨道计算(Hartree-Fock)
Geophysical Laboratory
第六章 自旋和角动量 复旦大学 苏汝铿.
柯红卫 兰州大学 Hong-Wei Ke and Xue-Qian Li, arXiv: and  arXiv:
计算机数学基础 主讲老师: 邓辉文.
§2 求导法则 2.1 求导数的四则运算法则 下面分三部分加以证明, 并同时给出相应的推论和例题 .
三、价层电子对互斥理论 基本要点: ABn分子或离子的几何构型取决于与中心A原子的价层电子对数目。 价层电子对=σ键电子对+孤对电子对
第8章 静电场 图为1930年E.O.劳伦斯制成的世界上第一台回旋加速器.
Positive gate bias-Induced Reliability in IGZO TFT
第4章 非线性规划 4.5 约束最优化方法 2019/4/6 山东大学 软件学院.
如何選擇計算方法 量子化學計算方法簡介.
分子建模与模拟导论 主讲: 王延颋,中国科学院理论物理研究所 助教: 石锐,中国科学院理论物理研究所
分子建模与模拟导论 主讲:王延颋, 中国科学院理论物理研究所 助教: 邓礼, 中国科学院理论物理研究所
分子建模与模拟导论 主讲: 王延颋,中国科学院理论物理研究所 助教: 张小虎,中国科学院理论物理研究所
第7讲 自旋与泡利原理.
过程自发变化的判据 能否用下列判据来判断? DU≤0 或 DH≤0 DS≥0.
3.8.1 代数法计算终点误差 终点误差公式和终点误差图及其应用 3.8 酸碱滴定的终点误差
概 率 统 计 主讲教师 叶宏 山东大学数学院.
线 性 代 数 厦门大学线性代数教学组 2019年4月24日6时8分 / 45.
WPT MRC. WPT MRC 由题目引出的几个问题 1.做MRC-WPT的多了,与其他文章的区别是什么? 2.Charging Control的手段是什么? 3.Power Reigon是什么东西?
3. 分子动力学 (Molecular Dynamics,MD) 算法
成绩是怎么算出来的? 16级第一学期半期考试成绩 班级 姓名 语文 数学 英语 政治 历史 地理 物理 化学 生物 总分 1 张三1 115
激光器的速率方程.
第4章 Excel电子表格制作软件 4.4 函数(一).
正切函数的图象和性质 周期函数定义: 一般地,对于函数 (x),如果存在一个非零常数T,使得当x取定义域内的每一个值时,都有
第16讲 相似矩阵与方阵的对角化 主要内容: 1.相似矩阵 2. 方阵的对角化.
第15章 量子力学(quantum mechanics) 初步
第四章 第四节 函数图形的描绘 一、渐近线 二、图形描绘的步骤 三 、作图举例.
An Introduction to Electronic Structure Calculation
相关与回归 非确定关系 在宏观上存在关系,但并未精确到可以用函数关系来表达。青少年身高与年龄,体重与体表面积 非确定关系:
一 测定气体分子速率分布的实验 实验装置 金属蒸汽 显示屏 狭缝 接抽气泵.
Skyrme 能量密度泛函 在重离子熔合反应中的应用
学习任务三 偏导数 结合一元函数的导数学习二元函数的偏导数是非常有用的. 要求了解二元函数的偏导数的定义, 掌握二元函数偏导数的计算.
第18 讲 配合物:晶体场理论.
海报题目 简介: 介绍此项仿真工作的目标和需要解决的问题。 可以添加合适的图片。
I. 第一性计算 (First Principles Calculations)
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
2019/5/20 第三节 高阶导数 1.
热力学第一定律的应用 --理想气体等容过程、定容摩尔热容 --理想气体等压过程 、定压摩尔热容.
§2 方阵的特征值与特征向量.
滤波减速器的体积优化 仵凡 Advanced Design Group.
正弦函数的性质与图像.
第四节 向量的乘积 一、两向量的数量积 二、两向量的向量积.
第三章 从概率分布函数的抽样 (Sampling from Probability Distribution Functions)
第四章 函数的 积分学 第七节 定积分的换元积分法     与分部积分法 一、定积分的换元积分法 二、定积分的分部积分法.
Volterra-Lotka方程 1925年, A. Lotka(美)和V. Volterra(意)给出了第一个两物种间的捕食模型。
热力学与统计物理 金晓峰 复旦大学物理系 /7/27.
本底对汞原子第一激发能测量的影响 钱振宇
第十七讲 密码执行(1).
§2 自由代数 定义19.7:设X是集合,G是一个T-代数,为X到G的函数,若对每个T-代数A和X到A的函数,都存在唯一的G到A的同态映射,使得=,则称G(更严格的说是(G,))是生成集X上的自由T-代数。X中的元素称为生成元。 A变, 变 变, 也变 对给定的 和A,是唯一的.
海报题目 简介: 介绍此项仿真工作的目标和需要解决的问题。 可以添加合适的图片。
Presentation transcript:

2. 密度泛函理论  Hohenberg-Kohn 定理 1 )基态系统的所有物理性质都由电子密度唯一决定,能量与电子密度为一一映射。 2 )对应于电子密度的变分原理:任意近似电子密度所对应的能量值都大于等于基 态对应的真正密度所决定的能量值。  密度泛函理论 ( Density Functional Theory, DFT) 虽然证明了电子密度和基态能量的一一对应关系是存在的,但是两者之间的泛函形 式未知。各种 DFT 的目的就是从不同的简化物理图象出发,给出近似的泛函形式。 E [ρ] = T [ρ] + E ne [ρ] + J [ρ] + K [ρ] E ne [ρ] (核子 - 电子势能)和 J [ρ] (库仑积分)已知, T [ρ] (动能)和 K [ρ] (交换积 分)未知。  Kohn-Sham 理论 DFT 的 HF 理论。给定了未知泛函的形式后,类似 HF 方法,得到准本征态方程 Kohn-Sham 方程。 HF 的计算量,但是自动包括了电子关联的贡献。

 Local Density Methods 假设局域电子密度可以被认为是均匀电子气,或等效地说,电子密度是随空间缓慢 变化的函数。 交换项 Local Density Approximation (LDA) Local Spin Density Approximation (LSDA) 关联项 Vosko , Wilk , and Nusair ( VWN )

GGA (见下)中的 PW91 修改了 VWN 的泛函形式:  Gradient Corrected Methods Gradient Corrected or Generalized Gradient Approximation (GGA): 泛函不仅决定 于电子密度,还决定于电子密度的梯度。 交换项 Perdew and Wang (PW86): 修正 LSDA 的泛函形式:加入高阶项。 Becke (B or B88): 正确的能量密度渐进行为。

Becke and Roussel (BR): 加入轨道波函数的导数项。 Perdew and Wang (PW91) 关联项 Lee, Yang, and Parr ( LYP )

Perdew ( P86 ):修正 LSDA 的梯度项。 Perdew and Wang ( PW91 or P91 ):改进 P86 。 其中 在 LSDA 部分已经给出。

Becke ( B95 ):更好地满足一些基本的物理约束。  混合方法 混合 HF 和 DFT 给出的能量项。 Becke 3 parameter functional (B3)

 交换和关联项的组合应用 SVWN = LSDA + VWN BLYP = B88 + LYP BP86 = B88 + P86 BPW91 = B88 + PW91 B3LYP = B3 + LYP B3P86 = B3 + P86 B3PW91 = B3 + PW91  一般而言, GGA 比 LSDA 效果要好得多。 GGA 的计算量与 HF 相仿,但构型和振 动频率的精确度一般要好于 MP2 ,与 CC 可比拟。 DFT 方法总的来说对静电作用描 述得更好一些,而对于范德华作用描述得差一些。  DFT 算法的实现与 HF 相似。基本的 DFT 算法复杂度为 M 4 ,最新计算技术使 DFT 的计算量线性化。  DFT 的最大问题在于没有统一的理论方法系统地提高计算精度,即更复杂的泛函形 式不一定计算精度越高,而是与被研究体系密切相关。  运用 DFT 计算的软件包之一: VASP (Vienna Ab-initio Simulation Package) 应用周期性边界条件以计算较大的体系。

-r +r  CeO 2 reduces SO 3 immediately  CuO has comparable energy costs for SO 3 reduction and SO 2 oxidation  SO 3 may contaminate CuO and TiO 2 surfaces by forming sulfate Energy Differences for Reduction of SO 3 * * Y. Wang, S. Rashkeev et al. to be submitted. Catalysis to accelerate: SO 3 -> SO 2 + ½ O 2

3. 第一性计算的应用举例 3.1. 数量分析( Population Analysis )  确定每个原子上有效的电子数目(一般不是整数)。一个重要的应用是给定每个 原子上的部分电荷( partial charge ),作为全原子模拟时经验力场的一部分。  基于原子轨道(基矢)的分析: Mulliken Population Analysis: 正交不归一的基矢。 Löwdin Population Analysis: 正交归一的基矢。  基于静电势( Electrostatic Potential, ESP )的分析 : 把分子周围由范德华半径至两到三倍的距离的三维空间范围离散化成格点,由第 一性计算得出格点上的静电势,用最小平方拟合法决定每个原子上的部分电荷。 Partial charges of Kapton unit

3.2. 化学反应过渡态( Transition State , TS )计算  反应势能面( Potential Energy Surface, PES ):除反应坐标之外的其它自由 度上系统都处在最低能量态。 R P 鞍点  过渡态理论( Tranistion State Theory, TST ):假设沿反应坐标的所有点都处 在热力学平衡态,因而系统处在某一状态的几率服从玻尔兹曼分布。  鞍点( Saddle Point ):即 TS ,沿反应坐标的极大值点。  Arrhenius Law :宏观化学反应速率,是 TS 和反应物之间的吉布斯自由能差值。 h 是 Planck 常数, k B 是 Boltzmann 常数。 R 是气体常数, T 是温度。  化学反应平衡常数: 是反应物和生成物之间的吉布斯自由能差值。

II. 分子建模  目的:把整个原子作为一个质点进行模拟(全原子模拟, Atomistic Simulation 或 All-atom Simulation ,也叫做 Force Field Method 或 Molecular Mechanics) , 以减少计算自由度,加大可计算的体系的空间和时间尺度,简化数据处理和 分析。  方法:一般的做法是根据原子间相互作用的物理特性,预先设定一个有待定参数的 二体或多体的相互作用的经验势的函数形式,然后根据第一性计算的数据或 实验结果拟合经验势的参数。  误差:因为描述体系的自由度被大大减少,全原子模型不可能重建系统的所有性 质。拟合参数时,往往选择一组最关心的物理性质进行拟合,以求误差尽 量小,而放松对其它性质的要求。所以要根据待研究的物理问题适当选取 全原子模型。

1. Lennard-Jones ( LJ )势  最常用的描述原子间范德华力的经验势。最广泛使用的是 12-6 LJ :  惰性气体的原子间相互作用仅用 LJ 就基本可以完全描述。 氩原子之间的相互作用, wikipedia  截断距离( cutoff distance ):对于 短程作用,大于 cutoff 的贡献是常数。 三维空间中,以上积分收敛的为短程作 用,发散的为长程作用。

 约化单位的换算  约化单位( reduced unit ) 数值模拟时使用的内部单位,需要乘上常数才能对应于实际体系的真实物理单位 (国际单位制, SI-units )。 1 )给定四个基本物理量的单位 : 长度 L ,质量 M ,时间 t ,电荷电量 Q 2 )计算其它物理量: 能量 温度 压力 质量密度 数量密度 介电常数 其中 k B 是 Boltzmann 常数, N A 是 Avogadro 常数。

因为金属中的价电子可以自由运动, 所以一般要用多体作用描述金属体系的力场。 2. 金属体系的力场  Glue Model  EAM (Embedded Atom Model) 其中 r ij 是两个原子间的距离,   是类型为  和  的原子之间的二体势,   是类型 为  的原子 j 产生的电子电量密度在 i 处的值, F  是一个嵌入函数,代表把类型为  的原子 i 嵌入电子云中需要的能量。 可以用于合金体系。以被广泛应用于多种金属及其合金。 只适用于单一金属。较好地平衡了表面和内部的结构和能量。

 成键作用( Bonded Interactions ): Bonds , Valence Angles , Dihedral Angles (Torsional Angles), Improper Dihedral Angles 3. 化学和生物体系的力场  非成键作用( Nonbonded Interactions ):范德华力和静电力 Pande Group in Stanford

Bond Dihedral Angle Valence Angle Improper Dihedral Angle Kapton unit

4. 粗粒化方法  类似于从第一性原理层面到全原子层面,粗粒化方法力图从全原子层面进一 步简化到粗粒化( coarse-graining , CG )层面,以期大大提高计算的时间和 空间尺度。  困难在于全原子层面上,原子间相互作用并不集中在局部。而在第一性层面 上,电子及其相互作用基本局限在相应的原子核周围。  不同的粗粒化方法着重于重建不同的物性,如结构或扩散特性等。  一些粗粒化方法假定作用势的函数形式,然后用全原子模拟的结果定参数。 另一类从结构函数( RDF )出发,反推出作用势。我们的方法从全原子作用 势出发,通过数学变换较严格地得到粗粒化力场。

粗 粒 化 方 法 I : Multiscale Coarse-Graining (MS-CG) *  事先做全原子模拟,得到全原子力场。  假设中心二体力的粗粒化力场形式。  最小化全原子力场和粗粒化力场之间的差值。 * W. Noid, P. Liu, Y. Wang et al. J. Chem. Phys. 128, (2008). 从全原子力场出发严格建立粗粒化力场,不预先设定力场的函数形式。  大大减少计算自由度。  系统扩散过程加快。  可以在粗粒化层面去除一些原子自由度 (如水分子)  保证较好地重建结构性质  理论中没有考虑可移植性!

差值 : 每个粗粒化点 : 有效力场形式 : 中心二体及线性假设 多维多项式,只有一个极值点! 从全原子模拟中获得 线 性 近 似 表 示 粗 粒 化 力 场线 性 近 似 表 示 粗 粒 化 力 场线 性 近 似 表 示 粗 粒 化 力 场线 性 近 似 表 示 粗 粒 化 力 场

差值 : 变分原理 :  用共轭梯度法由 Ψ 和 g d 求解力场。  先减去长程静电力  同时求解成键和非成键作用力 共 轭 梯 度 法 求 解 力 场 极 值共 轭 梯 度 法 求 解 力 场 极 值共 轭 梯 度 法 求 解 力 场 极 值共 轭 梯 度 法 求 解 力 场 极 值  只有一个极小值解 !  Ψ 可以用于判断最优粗粒化策略

粗 粒 化 方 法 II : Effective Force Coarse-Graining (EF-CG) * * Y. Wang, W. Noid, P. Liu, G. A. Voth Phys. Chem. Chem. Phys. 11, 2002 (2009).  显式计算原子间作用力  具有比 MS-CG 更好的可移植性