Water in an External Electric Field: Comparing Charge Distribution Methods Using ReaxFF Simulations

Area of Science:

• Computational chemistry
• Materials science
• Physical chemistry

Background:

• Reactive molecular dynamics simulations are crucial for understanding chemical processes under external electric fields.
• The charge equilibration (QEq) method is commonly used for charge distribution in ReaxFF simulations of electric fields.
• QEq's assumption of a contiguous conductor is problematic for insulators and semiconductors.

Purpose of the Study:

• To investigate the limitations of the charge equilibration (QEq) method for modeling external electric fields in reactive systems.
• To evaluate the performance of the atom-condensed Kohn-Sham density functional theory (ACKS2) approach for electric field simulations.
• To compare the physical realism and system-size convergence of QEq and ACKS2 methods.

Main Methods:

• Simulations using the ReaxFF reactive molecular dynamics method.
• Application of the charge equilibration (QEq) method for charge distribution.
• Implementation and application of the atom-condensed Kohn-Sham density functional theory (ACKS2) approach.
• Analysis of charge distributions and system-size convergence for water in an external electric field.

Main Results:

• QEq predicts unphysical charge distributions and properties that do not converge with system size for water in an electric field.
• The ACKS2 approach yields more physically realistic charge distributions.
• Properties calculated using ACKS2 demonstrate convergence with system size, unlike QEq.
• విద్యుత్ క్షేత్రాల ప్రభావాలను అధ్యయనం చేయడానికి ACKS2, QEq పద్ధతుల యొక్క ప్రయోజనాలు మరియు అప్రయోజనాలను ప్రదర్శిస్తుంది.

Conclusions:

• The QEq method has significant limitations when modeling external electric fields, particularly for insulating or semiconducting systems.
• The ACKS2 approach provides a more accurate and physically meaningful description of charge distribution in external electric fields.
• ACKS2 addresses key shortcomings of QEq, offering improved convergence and realism for reactive molecular dynamics simulations.