Hi Lammps users,
I’m trying to simulate the reaction between CaO and CO2 using reaxff. I used the primitive cell of CaO to expand the cell and then created a cylindrical channel in the middle, filling it with CO2 using Packmol. However, no matter how I modify the model, the pressure always shows as nan during the minimization process. Even if the minimization process succeeds, the pressure still shows as nan after several thousand steps in the subsequent NVT simulation.
My data file and the input file are as follows,
Input
units real
boundary p p p
atom_style charge
timestep 0.5
neighbor 2.0 bin
neigh_modify every 10 delay 0 check yes
read_data fin.lmp
set type 1 charge +2.0
set type 2 charge -2.0
set type 3 charge +4.0
pair_style reax/c NULL
pair_coeff * * ffield.reax.CHO C O Ca
fix 1 all qeq/reax 1 0.0 10 1e-6 reax/c
fix 2 all reax/c/bonds 1000 bonds.reax
min_style cg
minimize 1.0e-4 1.0e-6 100 1000
velocity all create 300.0 825757 mom yes rot yes dist gaussian
thermo 100
thermo_style custom step press pe temp ke
fix 3 all nvt temp 300 500 100
dump 2 all atom 100 dump.xyz
run 4000000
Data file
HEADER
8940 atoms
3 atom types
0.000000000000 59.561050000000 xlo xhi
0.000000000000 57.493000000000 ylo yhi
0.000000000000 61.046000000000 zlo zhi
Masses
1 40.07800000 # Ca
2 15.99900000 # O
3 12.01100000 # C
Atoms # charge
1 1 0.000000 1.230000000000 1.168000000000 1.185000000000
2 1 0.000000 1.182000000000 3.687000000000 3.625000000000
3 1 0.000000 3.565000000000 1.209000000000 3.573000000000
4 1 0.000000 3.590000000000 3.605000000000 1.312000000000
5 2 0.000000 3.713000000000 3.600000000000 3.638000000000
6 2 0.000000 3.633000000000 1.124000000000 1.152000000000
7 2 0.000000 1.125000000000 3.666000000000 1.240000000000
8 2 0.000000 1.165000000000 1.263000000000 3.577000000000
9 1 0.000000 5.943000000000 1.183000000000 1.313000000000
10 1 0.000000 6.004000000000 3.684000000000 3.593000000000
11 1 0.000000 8.365000000000 1.273000000000 3.602000000000
12 1 0.000000 8.371000000000 3.611000000000 1.228000000000
13 2 0.000000 8.396000000000 3.593000000000 3.600000000000
14 2 0.000000 8.348000000000 1.107000000000 1.197000000000
15 2 0.000000 5.978000000000 3.694000000000 1.237000000000
16 2 0.000000 6.071000000000 1.249000000000 3.686000000000
……
8641 3 0.000000 26.806000000000 35.156000000000 34.617000000000
8642 2 0.000000 25.918000000000 35.522000000000 33.885000000000
8643 2 0.000000 27.695000000000 34.791000000000 35.348000000000
8644 3 0.000000 52.560000000000 37.948000000000 16.923000000000
8645 2 0.000000 52.284000000000 38.234000000000 15.782000000000
8646 2 0.000000 52.835000000000 37.662000000000 18.064000000000
8647 3 0.000000 22.005000000000 42.718000000000 31.317000000000
8648 2 0.000000 22.132000000000 42.579000000000 30.124000000000
8649 2 0.000000 21.877000000000 42.858000000000 32.510000000000
8650 3 0.000000 1.898000000000 41.753000000000 20.749000000000
8651 2 0.000000 2.258000000000 40.826000000000 20.063000000000
8652 2 0.000000 1.539000000000 42.679000000000 21.436000000000