计算化学公社

标题: Al/PTFE复合材料在空气中燃烧机理的ReaxFF力场参数求助 [打印本页]

作者
Author: LaplaceL 时间: 2020-8-28 11:06
标题: Al/PTFE复合材料在空气中燃烧机理的ReaxFF力场参数求助
各位大佬好
我目前在做Al/PTFE复合材料在空气中燃烧机理的相关研究，想借助分子动力学模拟对某些性能参数做一些解释，在学习LAMMPS过程中，一直找不到合适的力场文件（Al/O/C/H/F)，在看文献过程中看到了“We used the ReaxFF40,41 reactive force field to investigate the chemical reactions of GO, GF, and Al with GO and GF. Because the entire system contains three different types of materials, multiple ReaxFF
parameters were merged into a single force field (i.e., Al/C/H/O, C/ H/F, and F-O interactions employed by Hong and van Duin,42 Rahnamoun and van Duin,43 and Kim et al.,44 respectively).这么一句话，作者将几个力场中的参数合并到一个力场中，只是怎么实现的呢？附上文章链接：https://dx.doi.org/10.1021/acsami.9b20397 。

同时请各位大佬帮我看看，这么in文件有没有问题，哪里需要修改呢？

############################################################################
# LAMMPS Input file for calculating for Al2O3 and -(CF2-CF2)- system    #
# wuth a ReaxFF potential                                              #
############################################################################
####Initialization##########################################################
units          real
atom_style       charge                                                    #模拟条件初始化#
boundary       p p p
neighbor       2.5 bin
neigh_modify    every 10 delay 0 check no
####define box and create atoms#############################################
read_data       data.Al2O3                                              #构建或读取初始模型#
####define interatomic potential############################################
####             reax args: hbcut hbnewflag tripflag precision
pair_style       reax/c NULL checkqeq yes                                  #反应性立场#
pair_coeff       * * ffield.reax C C F H C C O Al                         #反应力场参数输入，参数文件名为ffield.reax#
####define the calculation of the reaction energy###########################
compute reax all pair reax/c                                              #计算与反应有关的全部能量#
####define variable#########################################################
variable eb    equal c_reax[1]                                           #定义变量#
variable ea    equal c_reax[2]
variable elp    equal c_reax[3]
variable emol equal c_reax[4]
variable ev    equal c_reax[5]
variable epen equal c_reax[6]
variable ecoa equal c_reax[7]
variable ehb    equal c_reax[8]
variable et    equal c_reax[9]
variable eco    equal c_reax[10]
variable ew    equal c_reax[11]
variable ep    equal c_reax[12]
variable efi    equal c_reax[13]
variable eqeq equal c_reax[14]
####define initial velocity#################################################
velocity       all create 298.15 4928459                                  #给初速度#
####room temperature relaxation#############################################
fix             1 all reax/c/bonds 1000 bonds.reax                      #输出每个轨迹步系统内原子连接其他原子信息#
fix             2 all nvt temp 298.15 298.15 10.0
fix             3 all qeq/reax 1 0.0 10.0 1e-6 reax/c                   #对系统电荷平衡校正计算#
dump             1 all custom 1000 trj.cool id type q x y z                #输出系统内原子轨迹坐标#
dump_modify    1 sort 1
restart          100 restart1 restart2
thermo       1000
thermo_style    custom step cpu temp evdwl ecoul epair etotal press &
               v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb &
            v_et v_eco v_ew v_ep v_efi v_eqeq                         #输出系统各及与化学反应相关的宏观性质数据#
thermo_modify    flush yes
timestep       0.1                                                       #模拟步长#
run             1000000                                                 #模拟步数#
####以上过程为常温下驰豫####################################################
undump          1
unfix          2
unfix          3
####heating process to 3000k################################################
fix             4 all nvt temp 298.15 3000.15 10.0                      #采用NVT系综模拟#
fix             5 all qeq/reax 1 0.0 10.0 1e-6 reax/c                   #对系统电荷平衡校正计算#
dump             2 all custom 1000 trj.heating id type q x y z             #输出系统内原子轨迹坐标#
dump_modify    2 sort 1
restart          100 restart3 restart4
thermo       1000
thermo_style    custom step cpu temp evdwl ecoul epair etotal press &
               v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb &
            v_et v_eco v_ew v_ep v_efi v_eqeq                         #输出系统各及与化学反应相关的宏观性质数据#
thermo_modify    flush yes
timestep       1                                                       #模拟步长，可根据需要更改#
run             1000000                                                 #模拟步数，可根据自己情况更改#
####以上过程为线性升温至3000k###############################################
undump          2
unfix          4
unfix          5
####3000k simulation########################################################
fix             6 all nvt temp 3000.15 3000.15 10.0
fix             7 all temp/rescale 50 3000.15 3000.15 10.0 1.0
fix             8 all qeq/reax 1 0.0 10.0 1e-6 reax/c
dump             3 all custom 1000 trj.Al2O3_CF2CF2 id type q x y z
dump_modify    3 sort 1
restart          100 restart5 restart6
thermo       1000
thermo_style    custom step cpu temp evdwl ecoul epair etotal press &
               v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb &
            v_et v_eco v_ew v_ep v_efi v_eqeq
thermo_modify    flush yes
timestep       1
run             10000000
####以上过程为3000k条件下模拟###############################################

作者
Author: dummy1 时间: 2020-8-29 12:54
合并力场是需要各自的力场训练集。当然如果其中一个力场的所有参数不变的话它的训练集就不需要了（比如Al/C/H/O力场的参数不变，你只需要建立F与Al/C/H/O所有原子对的相互作用的训练集就可以了）。

作者
Author: SAI 时间: 2020-8-29 13:38
这是做产气材料吗

作者
Author: LaplaceL 时间: 2020-8-30 16:55

dummy1 发表于 2020-8-29 12:54
合并力场是需要各自的力场训练集。当然如果其中一个力场的所有参数不变的话它的训练集就不需要了（比如Al/C ...

这样子呀，感谢感谢

作者
Author: LaplaceL 时间: 2020-8-30 16:57

SAI 发表于 2020-8-29 13:38
这是做产气材料吗

作燃料的哈哈

作者
Author: zxs1127 时间: 2020-9-8 15:32
请问楼主最后的产物的数据统计是怎么实现的呢

作者
Author: ljpss43 时间: 2021-6-28 09:50
同学你好，我也在做Al-F-C-H-O系统的ReaxFF反应分子模拟研究，看到你发的帖子我备受启发，能否加微信进一步交流呢？我的微信号：18845105293。期待回复！

作者
Author: fuzsong 时间: 2022-1-7 16:19
您好，请咨询一下Al/O/C/H/F reaxff 力场问题。我是做界面摩擦化学的，自己做了些拟合但是感觉不理想。不知道您这Al/O/C/H/F reaxff 力场函数方不方便提供一下，我借鉴一下。非常感谢。fuzsong@163.com

作者
Author: funok 时间: 2022-1-14 22:06

zxs1127 发表于 2020-9-8 15:32
请问楼主最后的产物的数据统计是怎么实现的呢

购买 VARxMD就可以了

作者
Author: 黑色桃花 时间: 2022-1-28 11:51

dummy1 发表于 2020-8-29 12:54
合并力场是需要各自的力场训练集。当然如果其中一个力场的所有参数不变的话它的训练集就不需要了（比如Al/C ...

老师，那反应力场到底怎么开发呢，现在想学学具体操作，方便的话，您赐教

作者
Author: wtsadl000 时间: 2022-2-28 21:56
想问一下最新成果怎么样了，有没有发表相关论文

欢迎光临计算化学公社 (http://bbs.keinsci.com/)