| 查看: 6384 | 回复: 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开始上升 多谢~ |
回复此楼» 猜你喜欢
2027广东省杰青
已经有5人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子动力学的缺点
已经有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人回复
2楼2015-04-08 08:53:00
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
3楼2015-04-08 08:56:40
4楼2015-04-10 13:27:11
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
5楼2015-04-12 18:59:57
6楼2015-04-12 21:20:15
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
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326.5小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2015-04-13 12:09:09
happyzyf
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 962.5
- 帖子: 95
- 在线: 95.5小时
- 虫号: 1259643
- 注册: 2011-04-09
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
9楼2015-04-13 15:09:18
|
http://muchong.com/bbs/viewthread.php?tid=7106381&fpage=1&target=blank http://muchong.com/bbs/viewthread.php?tid=7017928&fpage=1&target=blank；楼主这个其实我也是新手一般大家都是先NVT跑一下，然后NVE跑，但是NVT中SMASS选哪个，一般取>0,详细的话估计的看它的原代码. |
» 本帖已获得的红花（最新10朵）
happyzyf10楼2015-04-14 08:53:06
