VASP

vasp的分子动力学模拟

VASP2010-01-1502:26:36阅读57评论0字号：大中小

vasp做分子动力学的好处，由于vasp是近些年开发的比较成熟的软件，在做电

子scf速度方面有较好的优势。

缺点：可选系综太少。

尽管如此，对于大多数有关分子动力学的任务还是可以胜任的。

主要使用的系综是NVT和NVE。

下面我将对主要参数进行介绍！

一般做分子动力学的时候都需要较多原子，一般都超过100个。

当原子数多的时候，k点实际就需要较少了。有的时候用一个k点就行，不过这都

需要严格的测试。通常超过200个原子的时候，用一个k点，即Gamma点就可

以了。

INCAR：

EDIFF一般来说，用1E-4或者1E-5都可以，这个参数只是对第一个离子步的

自洽影响大一些，对于长时间的分子动力学的模拟，精度小一点也无所谓，但不

能太小。

IBRION=0分子动力学模拟

IALGO=48一般用48，对于原子数较多，这个优化方式较好。

NSW=1000多少个时间步长。

POTIM=3时间步长,单位fs,通常1到3.

ISIF=2计算外界的压力.

NBLOCK=1多少个时间步长，写一次CONTCAR，CHG和CHGCAR，PCDAT.

KBLOCK=50NBLOCK*KBLOCK个步长写一次XDATCAR.

ISMEAR=-1费米迪拉克分布.

SIGMA=0.05单位:电子伏

NELMIN=8一般用6到8,最小的电子scf数.太少的话,收敛的不好.LREAL=A

APACO=10径向分布函数距离,单位是埃.

NPACO=200径向分布函数插的点数.

LCHARG=F尽量不写电荷密度,否则CHG文件太大.

TEBEG=300初始温度.

TEEND=300终态温度。不设的话，等于TEBEG.

SMASS-3NVEenmble;-1用来做模拟退火。大于0NVT系综。

【转】vasp的分子动力学模拟

★★★★★★★★

小木虫(金币+1):奖励一下，谢谢提供资源

uuv2010(金币+1):您是否可以做一个专题，详细讲讲怎么做？比如第一步需要干什么，第二步需要干什

么，结果怎么分析……如果能做一个这样完整的专题就太好了，不知道您是否有兴趣？2011-07-1318:20:12

uuv2010(金币+1):多谢提供资源2011-07-1617:39:55

uuv2010(金币+5,1ST强帖+1):多谢您的详细讲解！感谢就此专题与大家分享！2011-08-1218:25:12

vasp做分子动力学的好处，由于vasp是近些年开发的比较成熟的软件，在做电子

scf速度方面有较好的优势。

缺点：可选系综太少。

尽管如此，对于大多数有关分子动力学的任务还是可以胜任的。

主要使用的系综是NVT和NVE。

下面我将对主要参数进行介绍！

一般做分子动力学的时候都需要较多原子，一般都超过100个。

当原子数多的时候，k点实际就需要较少了。有的时候用一个k点就行，不过这都

需要严格的测试。通常超过200个原子的时候，用一个k点，即Gamma点就可

以了。

INCAR：

EDIFF一般来说，用1E-4或者1E-5都可以，这个参数只是对第一个离子步的

自洽影响大一些，对于长时间的分子动力学的模拟，精度小一点也无所谓，但不

能太小。

IBRION=0分子动力学模拟

IALGO=48一般用48，对于原子数较多，这个优化方式较好。

NSW=1000多少个时间步长。

POTIM=3时间步长,单位fs,通常1到3.

ISIF=2计算外界的压力.

NBLOCK=1多少个时间步长，写一次CONTCAR，CHG和CHGCAR，PCDAT.

KBLOCK=50NBLOCK*KBLOCK个步长写一次XDATCAR.

ISMEAR=-1费米迪拉克分布.

SIGMA=0.05单位:电子伏

NELMIN=8一般用6到8,最小的电子scf数.太少的话,收敛的不好.LREAL=A

APACO=10径向分布函数距离,单位是埃.

NPACO=200径向分布函数插的点数.

LCHARG=F尽量不写电荷密度,否则CHG文件太大.

TEBEG=300初始温度.

TEEND=300终态温度。不设的话，等于TEBEG.

SMA微信怎么退出 SS-3NVEenmble;-1用来做模拟退火。大于0NVT系综。

/////////////////////////////////////////////////////////////////////////////////////////人事的工作 ///

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

/////////////////////

1)收敛判据的选择

结构弛豫的判据一般有两种选择：能量和力。这两者是相关的，理想情况下，能

量收敛到基态，力也应该是收敛到平衡态的。但是数值计算过程上的差异导致以

二者为判据的收敛速度差异很大，力收敛速度绝大部分情况下都慢于能量收敛速

度。这是因为力的计算是在能量的基础上进行的，能量对坐标的一阶导数得到力。

计算量的增大和误差的传递导致力收敛慢。

到底是以能量为收敛判据，还是以力为收敛判据呢？关心能量的人，觉得以能量为判据就够了；关心力相关量的人，没有选择，只能用力作为收敛标准。对于超

胞体系的结构优化，文献大部分采用Gamma点做单点优化。这个时候即使采用

力为判据（EDIFFG=-0.02），在做静态自洽计算能量的时候，会发现，原本已经

收敛得好好的力在不少敏感位置还是超过了结构优化时设置的标准。这个时候，

是不是该怀疑对超胞仅做Gamma点结构优化的合理性呢？是不是要提高K点密

度再做结构优化呢

在我看来，这取决于所研究的问题的复杂程度。我们的计算从原胞开始，到超胞，

到掺杂结构，到吸附结构，到反应和解离。每一步都在增加复杂程度。结构优化

终点与初始结构是有关的，如果遇到对初始结构敏感的优化，那就头疼了。而且，

还要注意到，催化反应不仅与原子本身及其化学环境有关，还会与几何构型有关。

气固催化反应过程是电子的传递过程，也是分子拆分与重新组合的过程。如果优

化终点的构型不同，可能会导致化学反应的途径上的差异。仅从这一点来看，第

一性原理计算的复杂性，结果上的合理性判断都不是手册上写的那么简单。

对于涉及构型敏感性的结构优化过程，我觉得，以力作为收敛判据更合适。而且

需要在Gamma点优化的基础上再提高K点密度继续优化，直到静态自洽计算时

力也是达到收敛标准的。

(2)结构优化参数设置

结构优化，或者叫弛豫，是后续计算的基础。其收敛性受两个主要因素影响：初

始结构的合理性和弛豫参数的设置。

初始结构

初始结构包括原子堆积方式，和自旋、磁性、电荷、偶极等具有明确物理意义的

模型相关参数。比如掺杂，表面吸附，空位等结构，初始原子的距离，角度等的

设置需要有一定的经验积累。DFT计算短程强相互作用（相对于范德华力），如

果初始距离设置过远（如超过4埃），则明显导致收敛很慢甚至得到不合理的结

果。

比较好的设置方法可以参照键长。比如CO在O顶位的吸附，可以参照CO2中C-O

键长来设置（如增长20%）。也可以参照文献。记住一些常见键长，典型晶体中

原子间距离等参数，有助于提高初始结构设置的合理性。实在不行，可以先在小

体系上测试，然后再放到大体系中算。

弛豫参数

弛豫参数对收敛速度影响很大，这一点在计算工作没有全部铺开时可能不会觉察

到有什么不妥，反正就给NSW设置个“无穷大”的数，最后总会有结果的。但是，

时间是宝贵的，恰当的设置3小时就收敛的结果，不恰当的设置可能要一个白天

加一个黑夜。如果你赶文章或者赶着毕业，你就知道这意味这什么。

结构优化分电子迭代和离子弛豫两个嵌套的过程。电子迭代自洽的速度，有四个

响很大的因素：初始结构的合理性，k点密度，是否考虑自旋和高斯展宽（SIGMA）；

离子弛豫的收敛速度，有三个很大的影响因素：弛豫方法（IBRION）,步长（POTIM）

和收敛判据（EDIFFG）.

一般来说，针对理论催化的计算，初始结构都是不太合理的。因此一开始采用很

粗糙的优化（EDIFF=0.001，EDIFFG=-0.2），很低的K点密度(Gamma)，不考

虑自旋就可以了，这样NSW<60的设置就比较好。其它参数可以默认。

经过第一轮优化，就可以进入下一步细致的优化了。就我的经验，

EDIFF=1E-4,EDIFFG=-0.05，不考虑自旋，IBRION=2，其它默认，NSW=100;

跑完后可以设置IBRION=1，减小OPTIM（默认为0.5，可以设置0.2）继续优

化。

优化的时候让它自己闷头跑是不对的，经常看看中间过程，根据情况调节优化参

数是可以很好的提高优化速度。这个时候，提交两个以上的任务排队是好的方式，

一个在调整的时候，下一个可以接着运行，不会因为停下当前任务导致机器空闲。

无论结构优化还是静态自洽，电子步的收敛也常常让新手头痛。如果电子步不能

在40步内收敛，要么是参数设置的问题，要么是初始模型太糟糕（糟糕的不是一

点点）。

静态自洽过程电子步不收敛一般是参数设置有问题。这个时候，改变迭代算法

（ALGO），提高高斯展宽（SIGMA增加）,设置自洽延迟（NELMDL）都是不错的

方法。对于大体系比较难收敛的话，可以先调节AMIN,BMIX跑十关于海的古诗 多步，得到电荷

密度和波函数，再重新计算。实在没办法了，可以先放任它跑40步，没有收敛的

迹象的话，停下来，得到电荷密度和波函数后重新计算。一般都能在40步内收敛。

对于离子弛豫过程，不调节关系也不大。开始两个离子再回到从前歌词 步可能要跑满60步（默认

的），后面就会越来越快了。

总的说来，一般入门者，多看手册，多想多理解，多上机实践总结，比较容易提

高到一个熟练操作工的水平。

如果要想做到“精确打击”，做到能在问题始发的时候就立刻采取有效措施来解决，

就需要回归基础理论和计算方法上来了。

(3)优化结果对初始结构和“优化路径”的依赖

原子吸附问题不大，但是小分子吸附，存在初始构型上的差异。slab上水平放置，

还是垂直放置，可能导致收敛结果上的差异。根据H-K理论，理想情况下，优化

得到的应该是全局最小，但在数值计算的时候可能经常碰到不是全局最小的情况。

实际操作中发现，多个不同初始结构优化收敛后在能量和结构上存在一定差异。

为了加快收敛速度，特别是对于表面-分子吸附结构，初始放松约束，比如

EDIFF=1E-3，EDIFFG=-0.3,NSW=30可能是很好的设置。但是下面的情况应当慎

重：

EDIFF=1E-3;

EDIFFG=-0.1;！或者更小

NSW=500；！或者更大

电子步收敛约束较小，而离子步约束偏大，离子步数又很多，这种情况下，可能

导致的结果是结构弛豫到严重未知的区间。

再在这个基础上提高约束来优化，可能就是徒劳的了——结果不可逆转的偏向不正

常的区间。

好的做法，是对初始结构做比较松弛的约束，弛豫离子步NSW应该限制在一个较

小的数值内。EDIFF=1E-3的话，EDIFFG也最好是偏大一些，如-0.3而不是-0.1.这

样可以在较少的步数内达到初步收敛。

对于远离基态的初始结构，一开始在非常松弛的约束下跑若干离子步，时间上带

来的好处是很大的。对于100个原子的体系用vasp做Gamma点优化，如果一开

始就是正常优化（EDIFF=1E-4,EDIFFG=-0.02）设置，开始十个离子步可能都要

花上几个小时。如果这个时候才发现输入文件有错误，那下午的时间就白费了，

顺便带上晚上机器空转。

所以，我习惯的做法，是在初始几步优化后，会用xcrysden检查一下XDATCAR

中的数据，用生成，然后看看弛豫过程是不是按照设想的

那样。后续过程跑完一个收敛过程，就再检查一下。如此这般，才放心

的展开后续计算。

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

/////////////////////////////////////////////////////////////////////////////////////////////////////////////////

结构完全相同的、处于各种运动状态的、各自独立的系统的集合。全称为统计系

综。系综是用统计方法描述热力学系统的统计规律性时引入的一个基本概念；系

综是统计理论的一种表述方式；系综并不是实际的物体，构成系综的系统才是实

际物体。

研究气体热运动性质和规律的早期统计理论是气体动理论。统计物理学的研究对

象和研究方法与气体动理论有许多共同之处，为了避免气体动理论研究中的困难，

它不是以分子而是以由大量分子组成的整个热力学系统为统计的个体。系综理论

使统计物理成为普遍的微观统计理论。

系统的一种可能的运动状态，可用相宇中的一个相点表示，随着时间的推移，系

统的运动状态改变了，相应的相点在相宇中运动，描绘出一条轨迹，由大量系统

构成的系综则可表为相宇中大量相点的集合，随着时间的推移，各个相点分别沿

各自的轨迹运动，类似于流体的流动。系综并不是实际的物体，构成系综的系统

才是实际物体。约束条件是由一组外加宏观参量来表示。在平衡统计力学范畴下，

可以用来处理稳定系综。

一、常用系综分类

根据宏观约束条件，系综被分为以下几种：

1.正则系综（canonicalenmble），全称应为“宏观正则系综”，简写为NVT，

即表示具有确定的粒子数（N）、体积（V）、温度（T）。正则系综是蒙特卡罗

方法模拟处理的典型代表。假定N个粒子处在体积为V的盒子内，将其埋入温度

恒为T的热浴中。此时，总能量（E）和系统压强（P）可能在某一平均值附近起

伏变化。励志寓言故事 平衡体系代表封闭系统，与大热源热接触平衡的恒温系统。正则系综的

特征函数是亥姆霍兹自由能F（N,V,T）。

2.微正则系综（micro-canonicalenmble），简写为NVE，即表示具有确定的粒

子数（N）、体积（V）、总能量（E）。微正则系综广泛被应用在分子动力学模拟

中。假定N个粒子处在体积为V的盒子内，并固定总能量（E）。此时，系综的温

度（T）和系统压强（P）可能在某一平均值附近起伏变化。平衡体系为孤立系统，

与外界即无能量交换，也无粒子交换。微正则系综的特征函数是熵S（N,V,E）。

3.等温等压(constant-pressure,constant-temperature)，简写为NPT，即表示具

有确定的粒子数（N）、压强（P）、温度（T）。一般是在蒙特卡罗模拟中实现。

其总能量(E)和系统体积(V)可能存在起伏。体系是可移动系统壁情况下的恒温热

浴。特征函数是吉布斯自由能G(N,P,T)。

4.等压等焓(contant-pressure,constant-enthalpy)，简写为NPH，即表示具有确

定的粒子数（N）、压强（P）、焓（H）。由于由于H=E+PV，故在该系综下进

行模拟时要保持压力与焓值为固定，其调节技术的实现也有一定的难度，这种系

综在实际的分子动力学模拟中已经很少遇到了。

5.巨正则系综（grandcanonicalenmble），简写为VT，即表示具有确定的粒

体积（青铜器鉴定 V）、温度（T）和化学势（）。巨正则系综通常是蒙特卡罗模拟的对象和

手段。此时、系统能量（E）、压强（P）和粒子数（N）会在某一平均值附近有一

个起伏。体系是一个开放系统，与大热源大粒子源热接触平衡而具有恒定的T，。

特征函数是马休(Massieu)函数J（，V，T）。

二、系综调节

系综调节主要是指在进行分子动力学计算过程中，对温度和压力参数的调节，分

为调温技术和调压技术。

1.调温技术

在NVT系综或NPT系综中，即使在NVE系综模拟的平衡态中，也经常调整

温度到期望值。如果希望知道系统的平衡态性质怎样依赖于温度，那么就必须在

不同的温度下进行模拟。

目前实现对温度的调节有4种方式：速度标度、Berendn热浴、Gaussian热浴、

No—Hoover热浴。

2.调压技术

在等压模拟中，可以通过改变模拟原胞的三个方向或一个方向的尺寸来实现体积

的变化．类似于温度控制的方法，也有许多方法用于压力控制，总的来说有以下3

种技术：Berendn方法、Anderson方法、Parrinello-Rahman方法。

三、系综选择

原则上巨正则系综应用最广，但却不一定是最方便的。因为可以看到三种系综的

演化过程既是约束解除的过程，却是以增加变量为代价的，这也就增加了数学上

的复杂性。因此一般情况下如果不需求解。，则不必使用巨正则系综。

系综选择的基本原则为：

1.微正则系综能够简单的求得近独立，全同，定域粒子系统，并且每个粒子只

能有两个不同的可能状态，例如简单的铁磁，顺磁模型

2.微正则系综难求的系统，可用正则系综求解

3.当微正则和正则系综均难求时，可用巨正则系综求解

///////////////////////////////////////////////cad怎么移动 /////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

////////////////////////////////日本军部 ////////////////////////////////////////////////////////////

/////////////////////

系综（enmble）：在一定的宏观条件下，大量性质和结构完全相同的、处于各

种运动状态的、各自独立的系统的集合。全称为统计系综。系综是用统计方法描

述热力学系统的统计规律性时引入的一个基本概念；系综是统计理论的一种表述

方式。

目录

系综的性质

研究对象

常用的三个系综

原理编辑本段系综的性质

系综是假想的概念，并不是真实的客观实体。真正的实体是组成系综的一个

个系统，这些系统具有完全相同的力学性质。每个系统的微观状态可能相

同，也可能不同，但是处于平衡状态时，系综的平均值应该是确定的。编辑本段

研究对象

研究气体热运动性质和规律的早期统计理论是气体动理论。统计物理学的研

究对象和研究方法与气体动理论有许多共同之处，为了避免气体动理论研究中的

困难，它不是以分子而是以由大量分子组成的整个热力学系统为统计的个体。系

综理论使统计物理成为普遍的微观统计理论。编辑本段常用的三个系综

J.W.吉布斯把整个系统作为统计的个体，提出研究大量系统构成的系综在

相宇中的分布，克服了气体动理论的困难，建立了统计物理。在平衡态统计理论

中，对于能量和粒子数固定的孤立系统，采用微正则系综（NVE）；对于可以和大

热源交换能量但粒子数固定的系统，采用正则系综（NVT）；对于可以和大热源

交换能量和粒子的系统，采用巨正则系综(mVT)。这是三种常用的系综，各系综在

相宇中的分布密度函数均已得出。量子统计与经典统计的研究对象和研究方法相

同，在量子统计中系综概念仍然适用。区别在于量子统计认为微观粒子的运动遵

循量子力学规律而不是经典力学规律，微观运动状态具有不连续性，需用量子态

而不是相宇来描述。编辑本段原理

系统的一种可能的运动状态，可用相宇中的一个相点表示，随着时间的推移，

系统的运动状态改变了，相应的相点在相宇中运动，描绘出一条轨迹，由大量系

统构成的系综则可表为相宇中大量相点的集合，随着时间的推移，各个相点分别

沿各自的轨迹运动，类似于流体的流动。若系统具有s个自由度，则相宇是

以s个广义坐标p（详写为p、p2……ps）和s个广义动量q(详写为q1、q2……qs)

为直角坐标构成的2s维空间。在相宇内任一点（p，q）附近单位相体积元内的相

点数目D（p，q，t）称为密度函数。D(p，q，t）在整个相宇的积分等于全部相

点数，即等于系综所包含的全部系统数N，与时间t无关。定义(p,q,t)=D(p,q,t)/N，

称为系综的概率密度函数。（p，q，t）dpdq表示在t时刻出现在（p，q）点附

近相体积元dpdq内的相点数在全部相点数中所占的比值，即表示任一系统在t

时刻其运动状态处于（p，q）附近的相体积元dpdq内的概率。显然，概率密度

函数（p，q，t）满足归一化条件∫（p,q,t）dpdq=1。统计物理学的认为

系统的任意宏观量I（t）是相应微观量L（p，q）在一定宏观条件下对系统一切可

能的微观运动状态的统计平均值，即I（t）＝∫L（p，q）（p，q，t）dpdq。由

此可见，经典统计物理的基本课题是确定各种条件下系综的概率密度函数（p，

q，t），确定后，即可对相应的热力学系统的宏观性质作出统计描述。这就是统

计系综的方法。（p，q，t）的具体形式与系统所处的宏观状态有关。如果

系统处于平衡态，则＝（p，q）不显含时间t，在平衡态的系综理论中，由能

量、体积和粒子数都固定的系统构成的统计系综称为微正则系综；由与温度恒定

的大热源接触，具有确定粒子数和体积的系统构成的统计系综称为正则系综；由

与温度恒定的大热源和化学势恒定的大粒子源接触，具有确定体积的系统构成的

统计系综称为巨正则系综；由与温度恒定的大热源接触并通过无摩擦的活塞与恒

压强源接触，具有确定粒子数的系统构成的统计系综称为等温等压系综。上述各

种统计系综都有各自的概率密度函数。在微正则系综中，系统处于所有可能的微

观状态上的概率都相等，即概率密度是不随时间改变的常数，这就是等概率原理。

等概率原理是平衡态统计物理的基本假设，它的正确性由它的推论与实际相符而

得到肯定。由微正则系统可以推导出其它系综的概率分布函数的形式。

微正则系综是由许多具有相同能量，粒子数，体积的体系的集合。它是统计力学

系综的一种。其配分函数是在能量E_0上的能态密度。微正则系综是个简并度

下的正则系综，正则系综可以被分开进入子系综，每个子系综被对应到可能的能

量值且自身为另一些微正则系综。

[

LasteditedbyGina88on2011-8-3at15:46]

uuv2010(金币+2):多谢分享2011-08-1218:22:13

这是一个例子，也许对你有提示作用

#------INCAR------

IBRION=0#分子动力学ISTART=0#fromscratch

IALGO=48#如果发现自恰不收敛的情况可换成38。

NSW=5000#时间步长，一般在500步左右就能达到平衡，但要想充分的话，

最好3000以上。

NBLOCK=1;KBLOCK=50#每个NBLOCK步数后，对关联函数和DOS会被

计算，然后离子位置信息写

到相应的XDATCAR。在KBLOCK*NBLOCK步数后，对关联函数和DOS会被写

到PCDAT和DOSCAR。

ISYM=0#取消对称性。

SMASS=-3#-3NVEenmble,-1模拟退火.>0NVTenmble

POTIM=2.0#时间步长。一般1到3.这个最好跑一下NVE系综，看一下总

能是否守恒。

ISMEAR=-1;SIGMA=0.1#-1是电子费米达拉克占据，主要用来处理电子含

温。LREAL=A

ISIF=2#就算压力。

NPACO=200,APACO=10#计算对关联分布的。具体可以看manual，很简单。

NELMIN=6#电子最少的自恰步数，一般是4到8。

NBANDS=700#能带数，这个最好自己设，否则容易出问题。

LCHARG=F#尽量不写CHG，太占硬盘资源。

TEBEG=1000；TEEND=1000#初始温度和末态温度，默认是两者相当，所以

一般设TEBEG就可以。

##------END---------

PREC一般也不用设，默认即可。ENCUT会等于ENMAX。我能想起来的参数就

这些，有问题可以直接问

我。因为分子动力学是一个统计概念，当然原子越多越好的。所以我认为最少也

得算100个原子。一般100个

原子的话，k点用一个就可以。注意的是在分子动力学里，PSTRESS是没有的用

的。因为没有NPH或NPT

系综。再有分子动力学常用的几个概念就是，对关联函数，这个实验可测。写在

PCDAT里。速度自关联函数，

这个可以自己看一下公式，很简单，可以自己写程序做。

做分子动力学最好用vasp5，因为很多参数都不用十二月英文 测了，4.6好像还得测，

否则慢。还有实际vasp在输出温度是需要修正的，这个可以看看manual，很简

单，如果原子数多的话可以不

修正。还有压力也需要修正，但通常影响不大，可以直接用vasp的。

youzhizhe(金币+1):谢谢分享。如果是转载，请注明转载地址。2011-08-1309:53:45

cenwanglai:这个是goldfish专家的帖子2011-11-1722:13:23

很长一段时间以来，一直有人在问BOMD与CPMD究竟有什么不同？这里我就简

单说一些二者的区别。

实际上，我想大部分人理解和接触的第一原理分子动力学方法以CPMD居多。

CPMD，就是Car和Parrinello两个人作出的基于密度泛函的分子动力学方法，其

特点是在引入电子虚拟质量，将电子运动耦合到了运动方程中，每一步分子动力

学计算后，对电子结构的计算就不再需要自洽场迭代这一过程，因此可以大大节

省计算资源，在计算机还不是那么好的80，90年代是弥足珍贵的。CPMD的建立

开创了第一原理分子动力学方法的新时代，使其真正开始了实用化进程。

而BOMD，顾名思义，就是Born-Oppenheimer分子动力学。Born-Oppenheimer

近似，也就是绝热近似，指的是将电子和离子的求解分离开来，只处理离子的动

力学部分，而认为电子可以快速跟上电子的运动。其特点是每一步分子动力学计

算之后都需要对电子结构进行自洽场迭代，使电子达到基态。正式由于这个自洽

场迭代过程需要的计算量巨大，致使其一直没有达到广泛运用，直到90年代中后

期，计算机技术的发展，才使BOMD开始逐步被人们重视。

CPMD和BOMD各有优劣。CPMD虽然计算速度更快，不需要进行自洽场迭代计

算，计算量小，但是由于计算中并没有使体系真正达到基态，而只是尽量靠近基

态，因此其准确性对电子的虚拟质量这个参数的选取依赖程度很大，一旦计算参

数不对，得到的体系很可能远远偏离真实的势能面，得不到正确的动力学轨迹。

同时，为了保证其尽量靠近基态，分子动力学时间步长一般选得较小。而BOMD

虽然计算量大，但是由于每一步都保证系统达到基态，因此其分子动力学步长可

以取得较大，一般1fs到5fs都有可能。综合来说，如果能结合二者的优势，第一

原理分子动力学计算效率将大大提高，这也成为近年来其发展的重要方向。

偶曾经用CASTEP的MD算了一个晶格的结构相变问题，目前打算把手上那个表面的东西搞完后

就做点QMD的东西。现在偶每天除了看点文献外就是马丁的《ElectronicStructure》和《AbInitio

MolecularDynamics:BasicTheoryandAdvancedMethods》两本书对比看，感触颇多啊！

SMASS=1,2,3虽说都是NVT系统，但是，还是有区别的，SMASS是控制温度震荡的频率的。

要选择合适的值，以防止No-热浴和离子的运动脱耦，导致非正则的系宗。

这是我以前模拟退火的一个例子

SYSTEM=TiNiENCUT=550

ISTART=0;ICHARG=2

ISMEAR=-5VOSKOWN=1

NSW=500;IBRION=0SMASS=-1

ISIF=0;POTIM=0.2

NWRITE=1ISYM=0

TEBEG=500;TEEND=0IWAVPR=2

PREC=Accurate

一般做分子动力学的时候都需要较多原子，一般都超过100个。

当原子数多的时候，k点实际就需要较少了。有的时候用一个k点就行，不过这都需要严格的测

试。通常超过200个原子的时候，用一个k点，即Gamma点就可以了。

INCAR：

EDIFF一般来说，用1E-4或者1E-5都可以，这个参数只是对第一个离子步的自洽影响大一些，

对于长时间的分子动力学的模拟，精度小一点也无所谓，但不能太小。

IBRION=0分子动力学模拟

IALGO=48一般用48，对于原子数较多，这个优化方式较好。

NSW=1000多少个时间步长。

POTIM=3时间步长,单位fs,通常1到3.

ISIF=2计算外界的压力.

NBLOCK=1多少个时间步长，写一次CONTCAR，CHG和CHGCAR，PCDAT.

KBLOCK=50NBLOCK*KBLOCK个步长写一次XDATCAR.

ISMEAR=-1费米迪拉克分布.

SIGMA=0.05单位:电子伏

NELMIN=8一般用6到8,最小的电子scf数.太少的话,收敛的不好.LREAL=A

APACO=10径向分布函数距离,单位是埃.

NPACO=200径向分布函数插的点数.

LCHARG=F尽量不写电荷密度,否则CHG文件太大.

TEBEG=300初始温度.

TEEND=300终态温度。不设的话，等于TEBEG.

SMASS-3NVEenmble;-1用来做模拟退火。大于0NVT系综。

SMASS=1,2,3是没有区别的。都是NVTenmble。

SMASS只要是大于0就是NVT系综。

学习了，但是要指出你的一点小错误。

NBLOCK=1指的是多少个离子步（时间步长只是对于MD而言，对于其它的计算，NBLOCK也

是起控制作用的）写一次XDATCAR

KBLOCK=50NBLOCK*KBLOCK个离子步写一次PCDAT.

CONTCAR是每个离子步之后都会写出来的，但是会用新的把老的覆盖

CHG是在每10个离子步写一次，不会覆盖

CHGCAR是在任务正常结束之后才写的。

在vasp的说明书中关于模拟退火有这么一段话：

Usuallyasimulatedannealingrunismoreefficientifallmassareequal,since

thentheenergydissipatesmorequicklybetweendifferentvibrationalmodes.

tition

functionsremainsunaffectedbyachangeoftheionicmass.

我想知道这个POMASS是不是真的可以随便改，比如我想算H在Fe里的扩散，

我把H的POMASS设置成和Fe的一样行不行？

谢谢！

看了各位的回复让我壮了胆子试了一试，果然对能量没什么影响，计算效率稍有

提高。

得到的启示是：以后还是不要去钻研这种小东西，多花些时间在设计模型上才是

正道！！！

谢谢各位！

这个是降温的退火还需要平衡一段时间把SMASS设置为大于0然后把始末温度都设置成你

要的温度

可否把退火的步骤写下来？我只前做过退火，可后来发现这个算法不适合我做的题目。曾经在一

片PRB上看到有人是这样做的：

1.在T=2000K平衡一段时间（比如2ps)

2.在T=2000---->T=1500K降温

3.在T=1500K平衡一段时间（比如2ps)

4.然后在降温....

这在模拟固液相变的时候很成功。但是对于固体-固体相变是否合适不清楚，但是也有人做过。

思路是这样的：通过不断升温-降温-升温-降温....直到找到一种新的结构为止。

理论上来说，通过模拟退火是可以完全找到一种物质的平衡构型的，之所以不太用，可能和计算

量有关。毕竟需要构建超原胞，而且收敛速度不快。

uuv2010(金币+2):多谢分享2011-08-1218:22:02

MD计算中，比如模拟时长1ps，步长为1fs，每一飞秒后写离子坐标信息，结束

后，会写入1000步坐标信息，outcar2poscal就是从OUTCAR中提取这些信息，

产生1000个POSCAR。而CONTCAR中只保留了MD结束时的离子坐标信息。

挺有意思。

一般来说，想要得到每一步的离子位置信息都是通过XDATCAR得到的。

可以在INCAR中设置，NBLOCK=1，KBLOCK=1。

uuv2010(金币+2):多谢分享2011-08-1218:22:25

这个东西是看你需要的，还是看手册比较务实

随便给你整一个吧（仅参考）

ISMEAR=0

SIGMA=0.2

GGA=91

ISIF=2

ENMAX=400

IBRION=0

VOSKOWN=1

LMAXPAW=0

NPLANE=.TRUE.

LCHARG=.FALSE.

LREAL=A

NELMIN=4

BMIX=2.0

NSW=350

NBLOCK=1;KBLOCK=100

SMASS=-1

POTIM=2.00

TEBEG=1000;TEEND=300

APACO=10.0

我知道输出信息有原子坐标及其三维结构，且由原子坐标、速度等运动信息可以得到系统的热力

学性质及其动力学性质，想请教大家具体怎么样对OUTCARXDATCAR文件进行操作以得到这些

结构参数及径向分布函数信息呢？再次表示感谢...

一般OUTCAR里头会有坐标信息，选择倒数N个构型，利用linux指令提取坐标信息POSITION，

通过一定数学处理，一般是编个小程序吧，就能算出径向分布了