计算化学公社

标题: cp2k中关于NPT_F和NPT_I的选用和相应参数设置 [打印本页]

作者
Author: nusiew 时间: 2021-10-26 09:33
标题: cp2k中关于NPT_F和NPT_I的选用和相应参数设置
本帖最后由 nusiew 于 2021-10-26 09:38 编辑

大家好，我想用NPT来跑熔体体系的AIMD，看到cp2k的NPT关键词有NPT_F和NPT_I两个，如下，官网说NPT_I适用于各向同性的体系，气体、液体以及非晶物质等各向同性的体系是不是都是用NPT_I这个关键词？各向异性的要用NPT_F这个关键词？另外，如果我想要模拟1个标准大气压下的熔体，压强是不是设为1bar就可以了。最后，如下，inp文件除了要在&MOTION部分对ENSEMBLE和PRESSSURE进行声明外，inp文件的其他部分还需要写出那些关键词来匹配NPT的设置呀？谢谢大家！

NPT_F
constant temperature and pressure using a flexible cell
NPT_I
constant temperature and pressure using an isotropic cell

&MOTION
  &MD
ENSEMBLE NPT_I
STEPS 1000
TIMESTEP 1
TEMPERATURE 1800
&BAROSTAT
   PRESSURE 1
   TIMECON 1000
&END BAROSTAT
&THERMOSTAT
   &NOSE
      LENGTH 1
      YOSHIDA 3
      TIMECON 1000
      MTS 2
   &END NOSE
&END
  &END MD
&END MOTION

作者
Author: sobereva 时间: 2021-10-26 14:26
是

是

作者
Author: nusiew 时间: 2021-10-26 17:17

sobereva 发表于 2021-10-26 14:26
是

是

请问熔体体系跑NPT之前需要先跑一个NVT吗？这样做有没有什么道理呢？

作者
Author: ydeyun 时间: 2021-10-27 10:17
本帖最后由 ydeyun 于 2021-10-27 10:18 编辑

nusiew 发表于 2021-10-26 17:17
请问熔体体系跑NPT之前需要先跑一个NVT吗？这样做有没有什么道理呢？

社长之前的帖子讲过（http://bbs.keinsci.com/thread-22765-1-1.html）

一般平衡体系的动力学只需要三个过程：
能量极小化
预平衡，或者叫平衡相（equilibrium）。用NPT系综
产生相（production）。用NPT系综
预平衡和产生相放到一起作为一个动力学模拟也可以。但由于通常平衡相用berendsen压浴而产生相用PR压浴，所以通常分开做。

作者
Author: nusiew 时间: 2021-10-27 10:56

sobereva 发表于 2021-10-26 14:26
是

是

我发现用NPT_I跑出来的结果中，盒子在x, y, z三个方向上变化的数值是同步一样的？请问你们跑出来的是这样的吗？
(, 下载次数 Times of downloads: 21)

作者
Author: sobereva 时间: 2021-10-27 13:00

nusiew 发表于 2021-10-26 17:17
请问熔体体系跑NPT之前需要先跑一个NVT吗？这样做有没有什么道理呢？

无数人做NPT之前做NVT都是完全多余的，我在论坛回复得有十几次这种问题了，没有特殊理由不要做NVT，纯粹多此一举（除非初始结构真是垃圾得不得了、有特别严重不合理接触）

作者
Author: sobereva 时间: 2021-10-27 13:01

nusiew 发表于 2021-10-27 10:56
我发现用NPT_I跑出来的结果中，盒子在x, y, z三个方向上变化的数值是同步一样的？请问你们跑出来的是这样 ...

isotropic本来就是这样，不这样就不叫isotropic了

作者
Author: JCenter 时间: 2025-5-25 21:29

sobereva 发表于 2021-10-27 13:01
isotropic本来就是这样，不这样就不叫isotropic了

卢老师，有几个问题想请教您：
1.CP2K跑NPT时正常结束，未存在单步不收敛情况，是不是CELL_REF就没必要用了。我看我的盒子体积变化在250 立方埃左右。
2.CELL_REF使用后，格点质量发生变化，计算的能量等性质估计都受影响。那这种情况下，哪些计算量存在不准，不能使用和讨论呢。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)