求助：CP2K中使用界面模型（加真空层），如果需要考虑热膨胀，是选择NPT还是NVT系综？

XWH · 发表于 Post on 2021-11-24 23:25:03

CP2K模拟界面或表面模型时，三维周期或者Z方向非周期的话，考虑热膨胀，是使用NVT还是NPT系综比较好？如果选的是三维周期+真空层模拟NPT可行吗？

sobereva · 发表于 Post on 2021-11-26 11:36:35

热膨胀当然得NPT，要不然盒子尺寸都不变怎么体现晶胞的膨胀
平行于界面方向考虑周期性，垂直的方向不考虑周期性

XWH · 发表于 Post on 2021-11-26 21:15:11

本帖最后由 XWH 于 2021-11-26 21:17 编辑

sobereva 发表于 2021-11-26 11:36
热膨胀当然得NPT，要不然盒子尺寸都不变怎么体现晶胞的膨胀
平行于界面方向考虑周期性，垂直的方向不考虑 ...

您好，请问一下，如果在使用NPT-I且Z方向非周期时，出现Stress Tensor not tested for 2D systems，这种错误的话，我应该在STRESS_TENSOR 使用那种方式比较好点？或者应该是其他什么方法消除这种问题？
出现错误时使用的代码如下：
&FORCE_EVAL
STRESS_TENSOR ANALYTICAL
  METHOD Quickstep
  &SUBSYS
&CELL
   A    5.82640000    0.00000000    0.00000000
   B -2.92121033    5.04117716    0.00000000
   C    0.00000000    0.00000000 19.72530000
      ALPHA_BETA_GAMMA 90 90 120.0910
   PERIODIC XY #Direction of applied PBC (geometry aspect)
……
&POISSON
   PERIODIC XY #Direction(s) of PBC for calculating electrostatics
   PSOLVER ANALYTIC #The way to solve Poisson equation
&END POISSON

wolfli369 · 发表于 Post on 2021-11-27 09:47:36

XWH 发表于 2021-11-26 21:15
您好，请问一下，如果在使用NPT-I且Z方向非周期时，出现Stress Tensor not tested for 2D systems，这种 ...

https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL.html考虑一下PSOLVER参数

XWH · 发表于 Post on 2021-11-28 00:03:40

wolfli369 发表于 2021-11-27 09:47
https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL.html考虑一下PSOLVER参数

感谢！

sobereva · 发表于 Post on 2022-3-17 07:24:46

XWH 发表于 2021-11-26 21:15
您好，请问一下，如果在使用NPT-I且Z方向非周期时，出现Stress Tensor not tested for 2D systems，这种 ...

照常用三维周期性优化，但Z方向通过&CELL_OPT里的CONSTRAINT Z固定尺寸保证真空区的存在。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助：CP2K中使用界面模型（加真空层），如果需要考虑热膨胀，是选择NPT还是NVT系综？

浏览过的版块