Hi lammps users,
I’m trying to simulate a reactive hydrocarbon system at high temperature, I use the command:
fix spec all reax/c/species 1 1 1000 spec_C3flame.out
fix bond all reax/c/bonds 1000 bond_C3flame.reaxc
to obtain the bond order and the species information. I think we can derive all species from bond order, is that right? However, I find the species derived from the bond order file are different from the species file. For example, the species file of timestep 777000 are:
Timestep No_Moles No_Specs C15H9 C4H3 C18H10 CH3 CO C5H3 H2 C16H9 C2H C47H24 C2H2O C45H15 C24H14 C8H3 C2 H2O C6H4 C29H14 C20H10 C72H37O C28H12 CH4 C18H11 C73H35 C67H31O C18H9 C24H9 C17H10 CH2O C40H18 C55H26 C20H12 C27H9 C2H2 C4H2 C2H4 C22H11 C5H4 C22H12 C25H13 H C24H12 C4H4 C2H3 C4H5
777000 119 45 1 2 6 2 7 1 44 1 3 1 1 1 1 1 3 9 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 3 2 1 1 1 3 1 1 1
As seen in attachment (“part_of_species.out”), the number of molecules is 119. However, from the bond order information, the number of molecules is 111, as seen in attachment (“part_of bond”). By the way, the number of H radicals is 1 in the “part_of_species.out”, but we cannot find any atom with “nb” is 0 in the “part_of_bond”. So, I believe species information from the “reax/c/species” and “fix reax/c/bonds”, but I cannot understand the reason, can anyone help me? Thanks.
Any help is appreciated, and thank you in advance.
Jiang Xiao
part_of_species.out (465 Bytes)
part_of_bond (157 KB)