Charge equilibration for molecular dynamics simulations
01 Apr 1991-The Journal of Physical Chemistry (American Chemical Society)-Vol. 95, Iss: 8, pp 3358-3363
TL;DR: In this paper, an approach for predicting charge distributions in molecules for use in molecular dynamics simulations is presented, where an atomic chemical potential is constructed by using these quantities plus shielded electrostatic interactions between all charges.
Abstract: We report here an approach for predicting charge distributions in molecules for use in molecular dynamics simulations. The input data are experimental atomic ionization potentials, electron affinities, and atomic radii. An atomic chemical potential is constructed by using these quantities plus shielded electrostatic interactions between all charges. Requiring equal chemical potentials leads to equilibrium charges that depend upon geometry. This charge equilibration (QEq) approach leads to charges in excellent agreement with experimental dipole moments and with the atomic charges obtained from the electrostatic potentials of accurate ab initio calculations. QEq can be used to predict charges for any polymer, ceramic, semiconductor, or biological system, allowing extension of molecular dynamics studies to broad classes of new systems. The charges depend upon environment and change during molecular dynamics calculations. We indicate how this approach can also be used to predict infrared intensities, dielectric constants, and other charge-related properties.
Citations
More filters
TL;DR: In this paper, a force field for large-scale reactive chemical systems (1000s of atoms) is proposed. But the force field does not have Coulomb and Morse potentials to describe nonbond interactions between all atoms.
Abstract: To make practical the molecular dynamics simulation of large scale reactive chemical systems (1000s of atoms), we developed ReaxFF, a force field for reactive systems. ReaxFF uses a general relationship between bond distance and bond order on one hand and between bond order and bond energy on the other hand that leads to proper dissociation of bonds to separated atoms. Other valence terms present in the force field (angle and torsion) are defined in terms of the same bond orders so that all these terms go to zero smoothly as bonds break. In addition, ReaxFF has Coulomb and Morse (van der Waals) potentials to describe nonbond interactions between all atoms (no exclusions). These nonbond interactions are shielded at short range so that the Coulomb and van der Waals interactions become constant as Rij → 0. We report here the ReaxFF for hydrocarbons. The parameters were derived from quantum chemical calculations on bond dissociation and reactions of small molecules plus heat of formation and geometry data for...
4,455 citations
TL;DR: This chapter discusses the development of DFT as a tool for Calculating Atomic andMolecular Properties and its applications, as well as some of the fundamental and Computational aspects.
Abstract: I. Introduction: Conceptual vs Fundamental andComputational Aspects of DFT1793II. Fundamental and Computational Aspects of DFT 1795A. The Basics of DFT: The Hohenberg−KohnTheorems1795B. DFT as a Tool for Calculating Atomic andMolecular Properties: The Kohn−ShamEquations1796C. Electronic Chemical Potential andElectronegativity: Bridging Computational andConceptual DFT1797III. DFT-Based Concepts and Principles 1798A. General Scheme: Nalewajski’s ChargeSensitivity Analysis1798B. Concepts and Their Calculation 18001. Electronegativity and the ElectronicChemical Potential18002. Global Hardness and Softness 18023. The Electronic Fukui Function, LocalSoftness, and Softness Kernel18074. Local Hardness and Hardness Kernel 18135. The Molecular Shape FunctionsSimilarity 18146. The Nuclear Fukui Function and ItsDerivatives18167. Spin-Polarized Generalizations 18198. Solvent Effects 18209. Time Evolution of Reactivity Indices 1821C. Principles 18221. Sanderson’s Electronegativity EqualizationPrinciple18222. Pearson’s Hard and Soft Acids andBases Principle18253. The Maximum Hardness Principle 1829IV. Applications 1833A. Atoms and Functional Groups 1833B. Molecular Properties 18381. Dipole Moment, Hardness, Softness, andRelated Properties18382. Conformation 18403. Aromaticity 1840C. Reactivity 18421. Introduction 18422. Comparison of Intramolecular ReactivitySequences18443. Comparison of Intermolecular ReactivitySequences18494. Excited States 1857D. Clusters and Catalysis 1858V. Conclusions 1860VI. Glossary of Most Important Symbols andAcronyms1860VII. Acknowledgments 1861VIII. Note Added in Proof 1862IX. References 1865
3,890 citations
TL;DR: A continuous surface charge (CSC) approach is introduced that leads to a smooth and robust formalism for the PCM models and achieves a clear separation between "model" and "cavity" which, together with simple generalizations of modern integral codes, is all that is required for an extensible and efficient implementation of thePCM models.
Abstract: Continuum solvation models are appealing because of the simplified yet accurate description they provide of the solvent effect on a solute, described either by quantum mechanical or classical methods. The polarizable continuum model (PCM) family of solvation models is among the most widely used, although their application has been hampered by discontinuities and singularities arising from the discretization of the integral equations at the solute-solvent interface. In this contribution we introduce a continuous surface charge (CSC) approach that leads to a smooth and robust formalism for the PCM models. We start from the scheme proposed over ten years ago by York and Karplus and we generalize it in various ways, including the extension to analytic second derivatives with respect to atomic positions. We propose an optimal discrete representation of the integral operators required for the determination of the apparent surface charge. We achieve a clear separation between “model” and “cavity” which, together with simple generalizations of modern integral codes, is all that is required for an extensible and efficient implementation of the PCM models. Following this approach we are now able to introduce solvent effects on energies, structures, and vibrational frequencies (analytical first and second derivatives with respect to atomic coordinates), magnetic properties (derivatives with respect of magnetic field using GIAOs), and in the calculation more complex properties like frequency-dependent Raman activities, vibrational circular dichroism, and Raman optical activity.
2,033 citations
TL;DR: The General Utility Lattice Program (GULP) as discussed by the authors has been extended to include the ability to simulate polymers and surfaces, as well as adding many other new features, and the current status of the program is fully documented.
Abstract: The General Utility Lattice Program (GULP) has been extended to include the ability to simulate polymers and surfaces, as well as adding many other new features, and the current status of the program is fully documented. Both the background theory is described, as well as providing a concise review of some of the previous applications in order to demonstrate the range of its use. Examples are presented of work performed using the new compatibilities of the software, including the calculation of Born effective charges, mechanical properties as a function of applied pressure, calculation of frequency-dependent dielectric data, surface reconstructions of calcite and the performance of a linear-scaling algorithm for bond-order potentials.
1,987 citations
Cites methods from "Charge equilibration for molecular ..."
...The alternative, as used by Rappe and Goddard in their QEq method [41], is to use a damped Coulomb potential that allows for the fact that at short distances the interaction arises from the overlap of electron density, rather than from just simple point ions....
[...]
TL;DR: Several of the fundamental algorithms used in LAMMPS are described along with the design strategies which have made it flexible for both users and developers, and some capabilities recently added to the code which were enabled by this flexibility are highlighted.
1,956 citations