Re: [lammps-users] Transition from EEM to QEq when implementing ReaxFF in LAMMPS
Georg Heinze <georg.heinze@...6617...>
Mon, 03 Jul 2017 14:18:15 +0200
I am still wondering about the different charge equilibration methods
in LAMMPS and the Reax Code. The main difference between my results
and the results of Khalilov et al. is long range Coulomb repulsion
between oxygen molecules approaching a silicon nanowire and oxygen
atoms adsorbed on the wire. This repulsion is happening in my
simulation because of the global, instantaneous charge redistribution
done by QEq. The results of Khalilov et al. do not seem to show the
repulsion even though EEM also works globally and instantaneously.
Does anybody of you know how one can explain this difference?
Zitat von Georg Heinze <georg.heinze@...6617...>:
I am currently simulating the oxidation of silicon nanowires using
ReaxFF with Parameters by Buehler et al. (Phys. Rev. Lett. 96, 2006)
in LAMMPS. My first approach is to replicate calculations by
Khalilov et al. (Nanoscale 5, 2013) carried out using the Reax Code.
When comparing my data to the data obtained by Khalilov et al. I
noticed a difference that is most likely related to charges. Charge
equalization in LAMMPS is realized via QEq while in the Reax Code it
is realized via EEM. In the original QEq paper by Rappe and Goddard
(J. Phys. Chem. 95, 1991) it is mentioned that the charges obtained
via EEM are a factor 3-6 small in comparison to charges obtained via
QEq. In the original ReaxFF paper by van Duin et al. (J. Phys. Chem.
A 105, 2001) it is argued that QEq and EEM are very similar and that
the main difference between the two can be compensated by adjusting
a shielding parameter in the Coulomb Term of ReaxFF. The question I
have is whether and how the difference between QEq and EEM was
considered when ReaxFF was implemented in LAMMPS.