这个是降温的 退火还需要 平衡一段时间 把SMASS设置为大于0 然后把始末温度都设置成 你要的温度
可否把退火的步骤写下来？我只前做过退火，可后来发现这个算法不适合我做的题目。曾经在一片PRB上看到有人是这样做的： 1.在T=2000K平衡一段时间（比如2ps) 2.在T=2000 ---->T=1500K降温
3.在T=1500K平衡一段时间（比如2ps) 4.然后在降温 ....
这在模拟固液相变的时候很成功。但是对于固体-固体相变是否合适不清楚，但是也有人做过。
思路是这样的：通过不断升温-降温-升温-降温....直到找到一种新的结构为止。
理论上来说，通过模拟退火是可以完全找到一种物质的平衡构型的，之所以不太用，可能和计算量有关。毕竟需要构建超原胞，而且收敛速度不快。
uuv2010(金币+2): 多谢分享 2011-08-12 18:22:02
MD计算中，比如模拟时长1ps，步长为1fs，每一飞秒后写离子坐标信息，结束后，会写入1000步坐标信息，outcar2poscal就是从OUTCAR中提取这些信息，产生1000个POSCAR。而CONTCAR中只保留了MD结束时的离子坐标信息。
挺有意思。
一般来说，想要得到每一步的离子位置信息都是通过XDATCAR得到的。 可以在INCAR中设置，NBLOCK=1，KBLOCK=1。
uuv2010(金币+2): 多谢分享 2011-08-12 18: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
我知道输出信息有原子坐标及其三维结构，且由原子坐标、速度等运动信息可以得到系统的热力学性质及其动力学性质，想请教大家具体怎么样对OUTCAR\\XDATCAR文件进行操作以得到这些结构参数及径向分布函数信息呢？再次表示感谢...
一般OUTCAR里头会有坐标信息，选择倒数N个构型，利用linux指令提取坐标信息POSITION，通过一定数学处理，一般是编个小程序吧，就能算出径向分布了
相关推荐：
loading