| 查看: 5904 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
happyzyf铜虫 (小有名气)
|
[求助]
vasp计算分子动力学温度变化不正常 已有1人参与
|
|
|
各位大侠好！我再用vasp的分子动力学做一个模拟退火的过程，感觉结果有问题，OSZICAR里面的温度是从0开始的。INCAR设置如下 NPAR = 10 LPLANE = .TRUE. LREAL = AUTO ICHARG = 1 ENCUT = 300 PREC = Normal ALGO = VeryFast NELMIN = 4 NELM = 150 EDIFF = 1E-4 ISMEAR = 0 SIGMA = 0.1 # ISPIN = 2 # LWAVE = .FALSE. # LCHARG = .FALSE. IBRION = 0 NSW = 5000 POTIM = 2 SMASS = 3 TEBEG = 774 TEEND = 300 NBLOCK = 1 OSZICAR中温度从0开始上升，最高到1400多，然后开始下降，这中间还有震荡，不知道是什么原因。后来我改了INCAR里面的参数，固定体系始末温度都在300K，想让它平衡一段时间看看，设置如下 NPAR = 10 LPLANE = .TRUE. LREAL = AUTO ICHARG = 1 ENCUT = 300 PREC = Normal ALGO = VeryFast NELMIN = 4 NELM = 150 EDIFF = 1E-4 ISMEAR = 0 # SIGMA = 0.1 # ISPIN = 2 LWAVE = .FALSE. LCHARG = .FALSE. IBRION = 0 NSW = 3000 POTIM = 1 SMASS = 2 TEBEG = 300 TEEND = 300 NBLOCK = 1 KBLOCK = 50 OSZICAR里面温度仍然是从0开始上升 多谢~ |
回复此楼» 猜你喜欢
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有7人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子动力学的缺点
已经有14人回复
- 感觉VASP计算缺陷存在很大问题
已经有17人回复
- Vasp相比经典分子动力学的优缺点
已经有11人回复
- 分子动力学半径模拟
已经有3人回复
- vasp分子动力学报错
已经有16人回复
- 请教一个vasp分子动力学计算的问题
已经有9人回复
- VASP计算表面吸附 但是能量总是不收敛/。。。
已经有9人回复
- 提交VASP任务，总是排队，不让计算！！
已经有30人回复
- 相同的vasp版本，不同机器上计算同一个体系，能量相差近10ev
已经有17人回复
- vasp计算单个氧原子和分子的能量
已经有20人回复
- VASP计算缺陷形成能是否还要对帯隙进行修正？
已经有5人回复
- VASP分子动力学优化的INCAR如下，怎么CONTCAR和POSCAR没有什么变化啊？
已经有13人回复
- 分子动力学入门
已经有19人回复
- 【任务征集】VASP跑分子动力学数据处理方法
已经有20人回复
- 【活动】vasp计算分子动力学讨论合集及专贴（筹划中）
已经有33人回复
- 【求助】关于vasp进行分子动力学模拟的求助
已经有12人回复
- 【求助】分子动力学中体系温度不断增加?
已经有16人回复
- 【求助】请教vasp计算恒温第一原理分子动力学
已经有9人回复
- 【求助】vasp分子动力学模拟求助
已经有15人回复
- 【原创】vasp的分子动力学模拟
已经有44人回复
- 【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能
已经有19人回复
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
3楼2015-04-08 08:56:40
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
5楼2015-04-12 18:59:57
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
|
手册中说“For SMASS>=0 a canonical ensemble is simulated using the algorithm of Nosé. The Nosé mass controls the frequency of the temperature oscillations during the simulation (see [1, 2, 3]. For SMASS=0 Nosé-mass corresponding to period of 40 time steps will be chosen. The Nosé-mass should be set so that the induced temperature fluctuation show approximately the same frequencies as the typical 'phonon'-frequencies for the specific system. For liquids something like 'phonon'-frequencies might be obtained from the spectrum of the velocity auto-correlation function. If the ionic-frequencies differ by an order of magnitude from the frequencies of the induce temperature fluctuations Nosé thermostat and ionic movement might decouple leading to a non canonical ensemble. The frequency of the approximate temperature fluctuations induced by the Nosé-thermostat can be found on the OUTCAR file.”也就是说“SMASS大于等于0时，表示是用用Nose算法来模拟NVT系综。Nose mass控制了模拟过程中的温度针纺的频率。SMASS=0时，默认步长40 time step，需设置Nose mass使得温度震荡频率与体系声子频率一致。对于液体，可以从速度自修正函数中得到类似声子的频率。离子的频率与声子频率差一个数量级会导致系综不再是正则系综。温度震荡的近似频率可以从OUTCAR中得到。”但是里面仍然没有区分SMASS=0,1,2,3的情况。还有一个问题就是OUTCAR里面没找到手册里说的温度震荡的频率，麻烦告知一下好吗？谢谢~ |
7楼2015-04-13 09:17:56
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
9楼2015-04-13 15:09:18
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
送红花一朵 11楼2015-04-14 12:13:49