Structural and dynamical properties predicted by reactive force fields simulations for four common pure fluids at liquid and gaseous non-reactive conditions

Four common pure fluids were chosen to elucidate the reliability of reactive force fields in estimating bulk properties of selected molecular systems: CH₄, H₂O, CO₂ and H₂. The pure fluids are not expected to undergo chemical reactions at the conditions chosen for these simulations. The ‘combustion’ ReaxFF was chosen as reactive force field. In the case of water, we also considered the ‘aqueous’ ReaxFF model. The results were compared to data obtained implementing popular classic force fields. In the gas phase, it was found that simulations conducted using the ‘combustion’ ReaxFF formalism yield structural properties in reasonable good agreement with classic simulations for CO₂ and H₂, but not for CH₄ and H₂O. In the liquid phase, ‘combustion’ ReaxFF simulations reproduce reasonably well the structure obtained from classic simulations for CH₄, degrade for CO₂ and H₂, and are rather poor for H₂O. In the gas phase, the simulation results are compared to experimental second virial coefficient data. The ‘combustion’ ReaxFF simulations yield second virial coefficients that are not sufficiently negative for CH₄ and CO₂, and slightly too negative for H₂. The ‘combustion’ ReaxFF parameterisation induces too strong an effective attraction between water molecules, while the ‘aqueous’ ReaxFF yields a second virial coefficient that is in reasonable agreement with experiments. The ‘combustion’ ReaxFF parameterisation yields acceptable self-diffusion coefficients for gas-phase properties of CH₄, CO₂ and H₂. In the liquid phase, the results are good for CO₂, while the self-diffusion coefficient predicted for liquid CH₄ is slower, and that predicted for liquid H₂ is about nine times faster than those expected based on classic simulations. The ‘aqueous’ ReaxFF parameterisation yields good results for both the structure and the diffusion of both liquid and vapour water.