计算化学公社

标题: cp2k分区控制温度跑QMMM MD，如何选择系综 [打印本页]
作者
Author: Ninilol    时间: 2025-12-9 15:52
标题: cp2k分区控制温度跑QMMM MD，如何选择系综
各位老师好，请教几个关于cp2k分区控制温度跑QMMM MD的问题：
我跑的体系是一个蛋白和小分子的QMMM MD的体系，希望提高QM区的温度，但是MM区的温度依旧保持300K。我根据cp2k官方的manual设置langevin控制温度，inp文件如下代码块的部分，参考：Langevin Dynamics — CP2K documentation，但是langevin跑的是NVT系综。而我注意到在帖子：求助：分子动力学中系综如何选择 - 分子模拟 (Molecular Modeling) - 计算化学公社中，sob老师提到分子动力学若没有特殊理由，就不该用NVT，应该NPT到底。所以我这样的体系是不是不能使用cp2k的langevin控温的做法？如果不合适的话，有没有其他的办法可以分区控温跑QMMM MD？
  1. &MOTION
  2. &PRINT
  3.   &RESTART
  4.    &EACH
  5.     MD 1
  6.    &END EACH
  7.   &END RESTART
  8. &RESTART_HISTORY
  9.   &EACH
  10.   MD 100
  11.   &END EACH
  12. &END RESTART_HISTORY
  13. &TRAJECTORY
  14.    &EACH
  15.     MD 20
  16.    &END EACH
  17. &END TRAJECTORY
  18. #&VELOCITIES
  19.   #FILENAME ./MWW-TiO2-15-b-grimme-MD
  20.   #&END VELOCITIES
  21. &END PRINT
  22.   &MD
  23.     ENSEMBLE Langevin
  24.     STEPS 40000
  25.     TIMESTEP 0.2
  26.     TEMPERATURE 310.0
  27.         &LANGEVIN
  28.                 Gamma        0.1
  29.         &END LANGEVIN
  30.         &THERMAL_REGION
  31.                 DO_LANGEVIN_DEFAULT        .TRUE.
  32.                 &DEFINE_REGION
  33.                         LIST        3121..3144
  34.                         DO_LANGEVIN        .TRUE.
  35.                         TEMPERATURE       310
  36.                 &END DEFINE_REGION
  37.                 &PRINT
  38.                         &LANGEVIN_REGIONS
  39.                         &End LANGEVIN_REGIONS
  40.                         &TEMPERATURE
  41.                         &End TEMPERATURE
  42.                 &END PRINT
  43.         &END THERMAL_REGION
  44.     &THERMOSTAT
  45.       TYPE CSVR
  46.       &CSVR
  47.         TIMECON 10.
  48.       &END CSVR
  49.       REGION DEFINED
  50.       &DEFINE_REGION
  51.         MM_SUBSYS ATOMIC
  52.       &END DEFINE_REGION
  53.       &DEFINE_REGION
  54.         QM_SUBSYS ATOMIC
  55.       &END DEFINE_REGION
  56.     &END THERMOSTAT
  57.     &PRINT
  58.       FORCE_LAST
  59.     &END PRINT
  60.   &END MD
  61. &END MOTION
复制代码







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