计算化学公社

标题: 求助:CP2K中使用界面模型(加真空层),如果需要考虑热膨胀,是选择NPT还是NVT系综? [打印本页]

作者
Author:
XWH    时间: 2021-11-24 23:25
标题: 求助:CP2K中使用界面模型(加真空层),如果需要考虑热膨胀,是选择NPT还是NVT系综?
CP2K模拟界面或表面模型时,三维周期或者Z方向非周期的话,考虑热膨胀,是使用NVT还是NPT系综比较好?如果选的是三维周期+真空层模拟NPT可行吗?

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

作者
Author:
XWH    时间: 2021-11-26 21:15
本帖最后由 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
作者
Author:
wolfli369    时间: 2021-11-27 09:47
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参数
作者
Author:
XWH    时间: 2021-11-28 00:03
wolfli369 发表于 2021-11-27 09:47
https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL.html考虑一下PSOLVER参数

感谢!
作者
Author:
sobereva    时间: 2022-3-17 07:24
XWH 发表于 2021-11-26 21:15
您好,请问一下,如果在使用NPT-I且Z方向非周期时,出现Stress Tensor not tested for 2D systems,这种 ...

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3