用户工具

站点工具


adf:setchageformolecule

差别

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

到此差别页面的链接

两侧同时换到之前的修订记录前一修订版
后一修订版
前一修订版
adf:setchageformolecule [2020/10/11 17:58] – [一、为分子指定带电量] liu.junadf:setchageformolecule [2024/05/15 21:17] (当前版本) liu.jun
行 1: 行 1:
 ====== 如何为某部分原子(例如某个分子)设置固定电荷====== ====== 如何为某部分原子(例如某个分子)设置固定电荷======
- +注意该功能设置的电在整个分子动力学模拟过程中,都是保持不变的。 
-=====一、为分子指定带电量===== +  * [[adf:setchargeforreaxff2020|]] 
-有的时候我们希望在ReaxFF的模拟过程中,希望某些分子上带上一定的电量。这可以通过在run文件的geo部分设置相应关键字进行指定。我们以一个水分子和一个SO<sub>3</sub><sup>2-</sup>的体系为例说明。 +  * [[adf:setchargeforreaxff2019|2019及以前版本]]
- +
-我们希望SO<sub>3</sub>这四个原子总共带-2那么***.run文件中cat > geo <<eor到eor之间**,需要增加几行,**修改后的*.run文件**如下所示,是这一段增加几行以后的情况: +
- +
-<code bash> +
-......  +
-......  +
-*.run文件其他内容省略...  +
-cat > geo <<eor +
-BIOGRF 200 +
-MOLCHARGE 4 7 -2 +
-MOLCHARGE 1 3 0 +
-DESCRP chargeformolecules +
-REMARK Created by ADFinput +
-CRYSTX    10.00000   10.00000   10.00000   90.00000   90.00000   90.00000 +
-HETATM     1 O                  -1.13087   0.90935  -0.15491 O      1 1  0.0  +
-HETATM     2 H                  -1.06778   0.12386  -0.71599 H      1 1  0.0  +
-HETATM     3 H                  -0.46757   0.93562   0.54899 H      1 1  0.0  +
-HETATM     4 S                   1.19586  -1.51753   0.46753 S      1 1  0.0  +
-HETATM     5 O                   1.05268  -0.43047   1.46261 O      1 1  0.0  +
-HETATM     6 O                   0.32616  -1.51702  -0.72999 O      1 1  0.0  +
-HETATM     7 O                   2.23479  -2.55822   0.63107 O      1 1  0.0  +
-END +
- +
-eor +
-*.run文件其他内容省略...  +
-......  +
-</code> +
- +
-其中: +
-  * 1~3原子是水分子的,总带电量为0,因此设置为:MOLCHARGE 1 3 0 +
-  * 4~7原子是SO3的,总带电量为-2,因此设置为:MOLCHARGE 4 7 -2 +
-  * 必须为所有原子指定,不然会报错。例如本例中如果只设置了MOLCHARGE 4 7 -2,没有设置MOLCHARGE 1 3 0,就会报错。 +
-  * 该情况仅适用于动力学过程中电荷保持不变情况;如果分解了、和其他原子成键了,这个设定就没有太大意义了。 +
- +
-=====二、可能的用途===== +
-为孤立H原子或者OH、酸根设置带电量,可一定程度上模拟酸碱度环境影响。+
adf/setchageformolecule.1602410303.txt.gz · 最后更改: 2020/10/11 17:58 由 liu.jun

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