用户工具

站点工具


adf:bondboost

差别

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

到此差别页面的链接

两侧同时换到之前的修订记录前一修订版
后一修订版
前一修订版
adf:bondboost [2022/03/04 20:25] – [基于4种Region定义Bond Boost参数] liu.junadf:bondboost [2022/05/06 19:50] (当前版本) – [基于4种Region定义Bond Boost参数] liu.jun
行 64: 行 64:
 {{ :adf:2020bond_boost_10.png?650 }} {{ :adf:2020bond_boost_10.png?650 }}
  
-注意: +**Bond Boost参数** 
-  * Bond Boost lifetime,指每次作用多长时间?这里设置成了3万步,一般可以设置为2000或3000。 +  * Bond Boost lifetime,指每次Restraints作用持续步数上限。这里设置成了3万步,一般可以设置为2000或3000。该时间结束后,会再次检测是否有符合条件的环、链(环、链含义见下文),以施加Restraints。因此该值设置的较大,则Bond boost工作频率就比较低。 
-  * Number of instance,指同时加速反应数。 +  * Number of instance,指同时Restraints环、链(参看下文后有助于理解此处),该数值设置越大,计算速度越慢 
-  * 本例中,这些“原子对”的选择是循环、首尾相连的:1~2、2~3、3~4、4~1。但实际上也可以是链,例如1~2、2~3、3~4。系统会监测这些环或者链的实际情况,如果发现有符合置条件(在R<sub>min</sub>、R<sub>max</sub>范围内)的链或者环存在,下方所定义的Restraints(施加外力)就会生。如果有多满足条件,会选择其中距离之和最短的环或链施加外力+ 
 +**Atom types with region参数:** 
 +  * 本例中,这些“原子对”是循环、首尾相连的“环”:1~2、2~3、3~4、4~1。但也可以是,例如1~2、2~3、3~4。系统会监测环或者链相关“原子对”的实际距离,如果发现设定的几组R<sub>min</sub>、R<sub>max</sub>条件全部同时满足,下方所定义的Restraints(施加外力)就会作用。如果有多组环或者链满足条件,会选择其中距离之和最短的n(Number of instance)施加Restraints
   * 本例中,定义的是环,因此5即1。   * 本例中,定义的是环,因此5即1。
  
-然后设置具体的力的大小,以及方向(推/拉)+**Restraints参数**
  
 {{ :adf:2020bond_boost_11.png?650 }} {{ :adf:2020bond_boost_11.png?650 }}
  
-  * R<sub>0</sub>指目标键长,例如O-H这里设置目标键长是1.5埃,当环或链条件满足时,检查O-H距离,如果大于目标键长,将产生拉近效果,如果小于键长,将产生推开的效果因为这里定义的条件中,区间是3.0~5.5埃时,才会施加外力,因此实际上只会产生“拉力”,没有机会产生推力(因为距离小于1.5埃时,条件已经不具备了)。这种情况下,目标键长并不需要精确,估算一个值即可,因为这种情况下它主要影响外力的大小,对实际键长没有影响。设置适当的话,对反应机理也没有明显影响,只是提高碰撞几率+  * R<sub>0</sub>指目标键长,例如O-H这里设置目标键长是1.5埃,当环或链条件满足时,检查O-H具体距离,如果大于目标键长,将产生拉,如果小于键长,将产生推本例中,因为定义的条件中,当O-H距离在3.0~5.5埃时,才会施加外力,因此实际上只会产生“拉力”,没有机会产生推力(因为距离小于1.5埃时,条件已经不具备了)。这种情况下,目标键长并不需要精确,估算一个值即可,因为这种情况下它主要影响外力的大小,对实际键长没有影响。
   * 类似地,环、链条件满足时,2~3之间会产生外加推力,因为底部2~3的目标键长是2.5,比R<sub>max</sub>=1.5还大。   * 类似地,环、链条件满足时,2~3之间会产生外加推力,因为底部2~3的目标键长是2.5,比R<sub>max</sub>=1.5还大。
-  * 后面两个数值也是力大小公式中的数值,其中一个是正比项,另一个是指数项。公式参考:https://www.scm.com/doc/OldReaxFF/BondBoost.html +  * 后面两个数值也是力大小公式中的数值,其中F∞是正比项,k是指数项,调整F∞可以线性的调整力的大小。公式参考:https://www.scm.com/doc/OldReaxFF/BondBoost.html 
-  * **<color red>注意:Restraints定义的原子对,不需要与上面定义的环、链顺序一致,甚至可以没有任何关系。例如环定义为1~2,2~3,3~4,4~1,Restraints也可以直接施加到2~4上。</color>**这种情况的话,就要非常注意R<sub>0</sub>的精确设置了,因为环、链条件没有关于2~4的,因此任何键长都可能存在的,那推力、拉力都有可能产生。+  * **<color red>Restraints定义的原子对,不需要与上面定义的环、链顺序一致,甚至可以没有任何关系。例如环定义为1~2,2~3,3~4,4~1,Restraints也可以直接施加到2~4上。</color>**这种情况的话,就要非常注意R<sub>0</sub>的精确设置了,因为环、链条件没有关于2~4的,因此任何键长都可能存在的,那推力、拉力都有可能产生。
  
 +用户可以设置多组Bond boost,例如,一个键链搭配一个Restraints,一个键环搭配另一个Restraints。数量没有限制。
 =====一般性分子动力学参数设置===== =====一般性分子动力学参数设置=====
 {{ :adf:2020bond_boost_12.png?400 }} {{ :adf:2020bond_boost_12.png?400 }}
adf/bondboost.1646396717.txt.gz · 最后更改: 2022/03/04 20:25 由 liu.jun

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