adf:viscosity

差别

这里会显示出您选择的修订版和当前版本之间的差别。

到此差别页面的链接

--- adf:viscosity [2022/02/17 09:36] – [模型] liu.jun
+++ adf:viscosity [2023/08/30 10:16] (当前版本) – [模型] liu.jun
@@ 行 8: / 行 8: @@
   * //T//为温度
   * //P<sub>αβ</sub>//为压力张量的非对角项
-  * //<P<sub>αβ</sub>(0)P<sub>αβ</sub>(t)>//即自相关函数对所有时间原点、所有非对角分量的平均
+  * //<P<sub>αβ</sub>(0)P<sub>αβ</sub>(t)>//即压力张量自相关函数对所有时间原点、所有非对角分量的平均
-因为粘度的计算不涉及反应，因此除了使用ReaxFF外，也可以使用Force Field模块中其他力场（例如GAFF力场）进行分子动力学模拟。分子动力学的参数设置大同小异，只是Main面板的参数设置有所差别。
+因为粘度的计算不涉及反应，因此除了使用ReaxFF外，也可以使用Force Field模块中其他力场（例如GAFF力场），甚至也可以使用第一性原理进行分子动力学模拟，从而类似地生成粘度信息。分子动力学的参数设置大同小异，只是Main面板的参数设置有所差别。
 这里以苯为例，室温下苯的粘度实验值为0.5 mPa s。
@@ 行 18: / 行 18: @@
 {{ :adf:viscosity01.png?650 }}
-注意体系的大小，会严重影响粘度的计算数值，因此**体系规模应该大一些，这样粘度数值可靠性强一些**（统计振荡减弱）。
+注意体系的大小，会严重影响粘度的计算数值，因此**体系规模应该大一些，这样粘度数值可靠性强一些**（统涨洛误差减弱）。
 =====ReaxFF分子动力学参数设置=====
 基本的设置（周期性一定要确保为Bulk）：
 {{ :adf:viscosity02.png?650 }}
-由于不关心化学反应，因此步长不需要太小，1fs就可以了，但是**运动轨迹对粘度非常重要，因此每个5fs就需要保存一次轨迹（甚至可以更小）**，关于动量、角动量守恒的限制也可以去掉（起始温度也可以不设置，因为后面Thermostat会设置温度）：
+由于不关心化学反应，因此步长不需要太小，1fs就可以了，但是**运动轨迹对粘度非常重要，因此每隔5fs就需要保存一次轨迹（甚至可以更小）**，关于动量、角动量守恒的限制也可以去掉（起始温度也可以不设置，因为后面Thermostat会设置温度）：
 {{ :adf:viscosity03.png?650 }}
 温度的设置：
@@ 行 41: / 行 41: @@
   * Steps设置为300-12001
   * Property设置为Pressure Tensor
-  * Atoms:在Movie窗口Ctrl A选中所有原子，然后点击此处的➕
+  * Atoms:在Movie窗口Ctrl A选中所有原子，然后点击此处的➕。这里也可以留白不处理，也表示将所有原子纳入考虑
-  * Vector elements设置为4 5 6，表示P<sub>yz</sub>、P<sub>xz</sub>、P<sub>xy</sub> (当然1、2、3表示P<sub>xx</sub>、P<sub>yy</sub>、P<sub>zz</sub>)
+  * Vector elements设置为4 5 6，表示P<sub>yz</sub>、P<sub>xz</sub>、P<sub>xy</sub> (当然1、2、3表示P<sub>xx</sub>、P<sub>yy</sub>、P<sub>zz</sub>)，不过粘度只跟非对角项有关，所以这里必须填入4 5 6
   * Max ACF Step建议设置为3000帧，即大致为15 ps左右
@@ 行 59: / 行 59: @@
 与实验符合的极好。
+=====补充：QSPR预测粘度=====
+AMS的COSMO-RS模块，也能预测液体粘度：[[adf:viscositycosmo-rs]]

adf/viscosity.1645061788.txt.gz · 最后更改: 2022/02/17 09:36 由 liu.jun

© 2014-2022 费米科技（京ICP备14023855号）