cp2k分区控制温度跑QMMM MD，如何选择系综

Ninilol

各位老师好，请教几个关于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

复制代码