第一性原理分子动力学课件

Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,11/7/2009,#,*,第一性原理分子动力学,第一性原理分子动力学,1,第一性,原理,分子动力学的体系及应用,一、第一性原理分子动力学的发展情况,二、考虑电子自旋的第一性原理分子动力学体系,三、考虑电子自旋第一性原理分子动力学的应用,四、两个算例的讨论,第一性原理分子动力学的体系及应用一、第一性原理分子动力学的发,2,一、第一性原理分子动力学的发展情况,用分子动力学方法讨论金属材料的结构相变及力学性质，已经有很多文献报导，并且分子动力学的方法也因势函数的选取有很多种，诸如,Lennard-Jones,势分子动力学、,Morse,势分子动力学方法等。,第一性,原理,分子动力学的体系及应用,针对不同的材料，构建介观条件下的对势，取决于对材料介观结构的深刻理解，这给函势数的构建带来一定的困难，从而给分子动力学的模拟带来困难。,一、第一性原理分子动力学的发展情况 用分子动力学方,3,如果忽略压力改变，系统在热力学过程中可以用恒温恒压系综NPT描述。,第一性原理分子动力学的体系及应用,三、考虑电子自旋第一性原理分子动力学的应用,对应的Lagrangean函数为,波函数 和自旋波函数 满足完整约束条件,在用常规的分子动力学模拟具有周期边界条件的金属晶体时，如果质量中心固定，那么分子速度就是常数，根据Boltzmann 理论热力学过程等效于恒温过程。,第一性原理分子动力学的体系及应用问题,式中第一项是电子运动Hamiltonian量平均值；,第一性原理分子动力学的体系及应用,第一性原理分子动力学的体系及应用,波函数 和自旋波函数 满足完整约束条件,第一性原理分子动力学的体系及应用,考虑电子自旋的温度相关模型，则可以分析材料磁性质随温度变化的关系，同时也可以分析基于磁致伸性质的材料压磁效应。,和 是附加动能和势能项，这就是耦合原子运动和电,Phase transformation for NiTi Alloy,第一性原理分子动力学的体系及应用,第一性,原理,分子动力学的体系及应用问题,当前的一个趋势是从介观研究领域走向微观领域，采用第一性原理从头计算法分子动力学则是一个选项。,1965,年,Kohn,和,Sham,提出,Kohn-Sham,方程,1,，,标志着密度泛函理论的诞生。,1972,年,von Barth,和,Hedin,2,以及,Pant,和,Rajagopal,3,分别提出了自旋密度泛函理论（,Spin Density Functional Theory,SDFT,）。,1985,年,R.Car,和,M.Parrinello,首先提出“分子动力学和密度泛函理论的统一方法”,4,，由此得到三个运动方程：,如果忽略压力改变，系统在热力学过程中可以用恒温恒压系综NPT,4,第一性,原理,分子动力学的体系及应用,与密度泛函理论相似，区别仅在于在势函数中增加了,一项磁相互作用项，因此系统的,Hamiltonian,算符为,与之对应的,Kohn-Sham,方程，增加电子自旋密度,第一性原理分子动力学的体系及应用 与密度泛函理论相似，区别,5,第一性,原理,分子动力学的体系及应用问题,1985,年,William G.Hoover,提出正则动力学：平衡相空间分布的概念,建立了来自非牛顿力学的正则分布。,5,得到,Nos Hamiltonian,运动方程,1984,年,Shuichi Nos,提出了正则系综的分子动力学模拟方法，为第一性原理分子动力学计算提供了可选用的,NVE,，,NHP,和,NPT,系综。,6,第一性原理分子动力学的体系及应用问题 198,6,第一性,原理,分子动力学的体系及应用,当前的第一性原理分子动力学均为不包括电子自旋的分子动力学，也未见报道用将该方法用于材料的微结构分析，尤其是材料磁性质的研究。,二、考虑电子自旋的第一性原理分子动力学的体系,对应于考虑电子自旋的,Hamiltonian,函数，由此可以写出在,Born-Oppenheimer,等能面上的能量泛函,关于电子运动,第一性原理分子动力学的体系及应用 当前的第一,7,第一性,原理,分子动力学的体系及应用问题,关于电子自旋,对应的,Kohn-Sham,方程和,SDFT,方程为,第一性原理分子动力学的体系及应用问题关于电子自旋对应的Koh,8,第一性,原理,分子动力学的体系及应用问题,对应的,Lagrangean,函数为,波函数 和自旋波函数 满足完整约束条件,第一性原理分子动力学的体系及应用问题对应的Lagrangea,9,第一性,原理,分子动力学的体系及应用,对应于,Lagrangean,函数的,Hamiltonian,函数因此写为,式中第一项是电子运动,Hamiltonian,量平均值；第二项是,电子,自旋运动,Hamiltonian,量平均值；第三项是离子间的,Coulomb,作用势。,和 是附加动能和势能项，这就是耦合原子运动和电,子运动以及电子自旋的能量方程。,第一性原理分子动力学的体系及应用对应于Lagrangean函,10,第一性,原理,分子动力学的体系及应用,金属晶体结构，如果忽略热力学过程中外加应力场对晶体的作用，则可以用等温等压,Gibbs,正则系综,NPT,进行描述。,NPT,系综是,NVE,系综的扩展，因此定义相同。,用,S.Nos,正则系综分子动力学方法，可以将,Lagrangean,函数写为,(*),第一性原理分子动力学的体系及应用 金属晶体结构，如果忽,11,Phase transformation for NiTi Alloy,Born-Oppenheimer势能面等能面上的特定点，对恒压情况而言，可以当成分子动力学计算中的运动常量。,第一性原理分子动力学的体系及应用,对应于考虑电子自旋的Hamiltonian函数，由此可以写出在Born-Oppenheimer等能面上的能量泛函,二、考虑电子自旋的第一性原理分子动力学体系,与密度泛函理论相似，区别仅在于在势函数中增加了,中的微缺陷相关，基于第一性原理分子动力学，预计可以探究裂纹产生的微机理，以及相关的磁记忆特性。,第一性原理分子动力学的体系及应用,如果忽略压力改变，系统在热力学过程中可以用恒温恒压系综NPT描述。,第一性原理分子动力学的体系及应用问题,Phase transformation for NiTi,12,第一性,原理,分子动力学的体系及应用,(*),式,定义了一个广义动能,和势能,第一性原理分子动力学的体系及应用(*)式定义了一个广义动能,13,在平衡态 、和 下，,经典动能 可以通过动力学方程产生的轨迹对其求关于时间的平均值，并且用与此系统相适应的归一化条件得到。,变换速度 、和 ，就可以改变温度,。,第一性,原理,分子动力学的体系及应用,当温度,系统达到平衡能量最小状态。,这样建立的模型就是随温度变化的模型。,在平衡态 、,14,第一性,原理,分子动力学的体系及应用,在讨论铁素体和奥氏体相变时，由于各部分存在能量交换，而每部分在热力学平衡态可以用,Gibbs,正则系综理论进行描述。在用常规的分子动力学模拟具有周期边界条件的金属晶体时，如果质量中心固定，那么分子速度就是常数，根据,Boltzmann,理论热力学过程等效于恒温过程。如果忽略压力改变，系统在热力学过程中可以用恒温恒压系综,NPT,描述。这时位于,Born-Oppenheimer,势能面等能面上的特定点，对恒压情况而言，可以当成分子动力学计算中的运动常量。,第一性原理分子动力学的体系及应用 在讨论铁素体和奥氏体,15,第一性原理分子动力学课件,16,第一性,原理,分子动力学的体系及应用,相应的配分函数为,第一性原理分子动力学的体系及应用相应的配分函数为,17,第一性,原理,分子动力学的体系及应用,第一性原理分子动力学的体系及应用,18,第一性,原理,分子动力学的体系及应用,三、考虑电子自旋第一性原理分子动力学的应用,考虑电子自旋第一性原理分子动力学是基于量子力学的算法，因此可以对材料的微观结构进行分析，如相变分析以及磁性计算等。如果某些材料的磁性质是次要的，则可以不考虑电子自旋的影响。考虑电子自旋的温度相关模型，则可以分析材料磁性质随温度变化的关系，同时也可以分析基于磁致伸性质的材料压磁效应。,材料的断裂力学性质，是与材料的缺陷，尤其是材料,中的微缺陷相关，基于第一性原理分子动力学，预计可以探究裂纹产生的微机理，以及相关的磁记忆特性。,第一性原理分子动力学的体系及应用三、考虑电子自旋第一性原理分,19,第一性,原理,分子动力学的体系及应用,四、两个算例的讨论,1.,铁素体与奥氏体相变及磁性的从头计算动力学分析,2.,Ab Initio Molecular Dynamics Study of B2-B19,Phase transformation for NiTi Alloy,第一性原理分子动力学的体系及应用四、两个算例的讨论1.铁素体,20,第一性原理分子动力学课件,21,第一性,原理,分子动力学的体系及应用,谢谢！,第一性原理分子动力学的体系及应用谢谢！,22,