adf:temperaturecontrol
差别
这里会显示出您选择的修订版和当前版本之间的差别。
| 两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
| adf:temperaturecontrol [2017/01/15 21:43] – liu.jun | adf:temperaturecontrol [2020/09/01 13:06] (当前版本) – [如何模拟升温-保温-降温-保温等过程] liu.jun | ||
|---|---|---|---|
| 行 3: | 行 3: | ||
| 1，主面板的参数设置如下图所示： | 1，主面板的参数设置如下图所示： | ||
| - | {{ :adf:reaxfftmp01.png |}} | + | {{ :adf:reaxtmp01.png?650 }} |
| 需要注意： | 需要注意： | ||
| - | * 总的步数一定是大于第二步中设置的步数。否则第二步的设置就没有意义了; | + | * Number of Iterations是MD模拟总的步数，一定是大于后面设置温度变化中设置的步数。否则第二步的设置就没有意义了; |
| - | * Main这一页里面设置的温度优先级低于Model > Temperature Regime中温度区间里面设置的温度，因此在温度区间里面设置过的温度，则以温度区间的设置为准，没有设置过的时间里面，以Main里面设置的温度为准； | + | * Main > Temperature里面设置的温度优先级低于Model > Temperature Regime中温度区间里面设置的温度，因此在温度区间里面设置过的温度，则以温度区间的设置为准，没有设置过的时间里面，以Main里面设置的温度为准； |
| - | * 之所以选择NVT系综，是因为这个系综对于微观系统的稳定性比较好。并且，虽然是不断变温的过程，但每个温度台阶，实际上都是一个NVT系综。具体可以参考[[adf: | + | * 之所以选择NVT系综，是因为这个系综对于微观系统的稳定性比较好。并且，虽然是不断变温的过程，但每个温度台阶，实际上都是一个NVT系综。具体可以参考[[adf: |
| 2，设置温度区间： | 2，设置温度区间： | ||
| 行 15: | 行 15: | ||
| 如下图： | 如下图： | ||
| - | {{ : | + | {{ : |
| 第一行设置，表示第0步开始，所有的原子，温度稳定为298K，在大约30fs内调整到该温度，之后温度保持不变，即温度变化梯度（dT）为0。 | 第一行设置，表示第0步开始，所有的原子，温度稳定为298K，在大约30fs内调整到该温度，之后温度保持不变，即温度变化梯度（dT）为0。 | ||
| - | 第二行设置，表示从2500步（也就是625fs）开始，温度从前面的298K升高，升温速率达约0.05K／步，这里的398K实际上不起作用，只有升温梯度0.05K／步起作用，一直保持这个速率升温。 | + | 第二行设置，表示从2500步（也就是625fs）开始，温度从298K升高，升温速率达约0.05K／步，一直保持这个速率升温。那么到10000步时温度为773K。 |
| - | 第三行设置，表示升温到1万步（也就是2500fs）的时候，温度变为398K左右，之后温度变化速率为0K／步。 | + | 第三行设置，表示升温到1万步（也就是2500fs）的时候，温度突变为398K左右，之后温度变化速率为0K/步。 |
| - | 第四行设置，表示从1.3万步（也就是3250fs）开始，温度在前面398K的基础上下降，下降速率为0.5K／步，负号表示下降。 | + | 第四行设置，表示从1.3万步（也就是3250fs）开始，温度从398K下降，下降速率为0.5K/步，负号表示下降。 |
| - | 总结：如果温度升高或者降低，那么温度设置不起作用，速率起作用，起始温度为当前温度；如果保温，那么温度设置起作用，速率设为0。 | + | Damp的设置不太重要，一般设置100fs以内，对结果影响不大。 |
| - | + | =====脚本中的设置===== | |
| - | 如此，得到的温度变化曲线，和上面的设置对照起来看： | + | 图形界面中的这些设置，会自动更新到对应的*.run文件中，最后这个例子里面，*.run文件中，关于控温这部分的字段为： |
| - | + | <code bash> | |
| - | {{ : | + | cat > tregime.in <<eor |
| - | + | # | |
| - | 为了便于理解，我们用另外一个更适合实际用途的例子来演示设置和结果（该例子，大致的估算了升温、降温的截止时间，以便于在何时的时间内设置保温）： | + | #start #nzones at1 at2 T1 Tdamp1 dT1 ... |
| - | + | 0 1 1 1500 298.0 30.0 0.0 | |
| - | {{ : | + | 2500 1 1 1500 398.0 30.0 0.05 |
| - | + | 10000 1 1 1500 398.0 30.0 0.0 | |
| - | {{ : | + | 13000 1 1 1500 298.0 30.0 -0.05 |
| - | <WRAP center round download | + | eor |
| - | + | </code> | |
| - | [[https:// | + | |
| - | </WRAP> | + | |
adf/temperaturecontrol.1484487801.txt.gz · 最后更改: 2017/01/15 21:43 由 liu.jun