Extending and simplifying the electronegativity equalization method1Dedicated to Professor Herman van Bekkum on the occasion of his 65th birthday.1

Block Chemistry for Accurate Modeling of Epoxy Resins

Accurate molecular modeling of the physical and chemical behavior of highly cross-linked epoxy resins at the atomistic scale is important for the design of new property-optimized materials. However, a systematic approach to parametrizing and characterizing these systems in molecular dynamics is missing. We therefore present a unified scheme to derive atomic charges for amine-based epoxy resins, in agreement with the AMBER force field, based on defining reactive fragments─blocks─building the network. The approach is applicable to all stages of curing from pure liquid to gelation to fully cured glass. We utilize this approach to study DGEBA/DDS epoxy systems, incorporating dynamic topology changes into atomistic molecular dynamics simulations of the curing reaction with 127,000 atoms. We study size effects in our simulations and predict the gel point utilizing a rigorous percolation theory to recover accurately the experimental data. Furthermore, we observe excellent agreement between the estimated and the experimentally determined glass transition temperatures as a function of curing rate. Finally, we demonstrate the quality of our model by the prediction of the elastic modulus based on uniaxial tensile tests. The presented scheme paves the way for a broadly consistent approach for modeling and characterizing all amine-based epoxy resins.

Unifying Charge-Flow Polarization Models

We show that several models where electric polarization in molecular systems is modeled by charge-flow between atoms can all be considered as different manifestations of a general underlying mathematical structure. The models can be classified according to whether they employ atomic or bond parameters and whether they employ atom/bond hardness or softness. We show that an ab initio calculated charge response kernel can be considered as the inverse screened Coulombic matrix projected onto the zero-charge subspace, and this may provide a method for deriving charge screening functions to be used in force fields. The analysis suggests that some models contain redundancies, and we argue that a parameterization of charge-flow models in terms of bond softness is preferable as it depends on local quantities and decay to zero upon bond dissociation, while bond hardness depends on global quantities and increases toward infinity upon bond dissociation.

Redox-Transition from Irreversible to Reversible Vitamin C by Pore Confinement in Microporous Carbon Network

Enhancement of redox-reversibility in electroactive species has been studied because of fundamental interest and their importance for energy storage systems. Various electroactive molecules suffer from redox-irreversible behavior, and this is a critical reason for their exclusion as redox electrolytes in energy storage systems. In this article, we fully demonstrated that ascorbic acid (ASC), which is an abundant but redox-irreversible molecule, can become redox-reversible when it is confined in microporous carbon regimes. From a theoretical perspective, redox-reversibility in an electrochemical reaction coupled with an irreversible chemical process can be greatly enhanced due to kinetic acceleration toward the inverse direction of the chemical reaction by accumulation of products in the nanoconfined regime. However, the kinetic acceleration in a nanoconfined domain shows limitations for enhancing the redox-reversibility, which indicates that stabilization of the species undergoing an irreversible chemical process is another important factor for redox-reversibility enhancement. The origin of nanoporous confinement of ASC and its enhanced redox-reversibility was rationalized by molecular dynamics simulations. We found that ASC-clusters of a fully protonated ASC and its conjugated base formed inside carbon pores, which would be a main driving force for its confinement in microporous carbon networks. Lastly, we demonstrated a prototype energy storage device using redox-reversible ASC in microporous carbon as the half electrode, which shows the feasibility of ASC as a possible redox electrolyte in an aqueous energy storage system.

Generalized Method for Charge-Transfer Equilibration in Reactive Molecular Dynamics

Variable charge models (e.g., electronegativity equalization method (EEM), charge equilibration (QEq), electrostatic plus (ES+)) used in reactive molecular dynamics simulations often inherently impose a global charge transfer between atoms (approximating each system as an ideal metal). Consequently, most surface processes (e.g., adsorption, desorption, deposition, sputtering) are affected, potentially causing dubious dynamics. This issue has been addressed by certain split charge variants (i.e., split charge equilibration (SQE), redoxSQE) through a distance-dependent bond hardness, by the atomic charge ACKS2 and QTPIE models, which are based on the Kohn-Sham density functional theory, as well as by an electronegativity screening extension to the QEq model (approximating each system as an ideal insulator). In a brief review of the QEq and the QTPIE model, their applicability for studying surface interactions is assessed in this work. Following this evaluation, a revised generalization of the QEq and QTPIE models is proposed and formulated, called the charge-transfer equilibration model or in short the QTE model. This method is based on the equilibration of charge-transfer variables, which locally constrain the split charge transfer per unit time (i.e., due to overlapping orbitals) without any kind of bond hardness specification. Furthermore, a formalism relying solely on atomic charges is obtained by a respective transformation, employing an extended Lagrangian method. We moreover propose a mirror boundary condition and its implementation to accelerate surface investigations. The models proposed in this work facilitate reactive molecular dynamics simulations, which describe various materials and surface phenomena appropriately.

AtomicChargeCalculator: interactive web-based calculation of atomic charges in large biomolecular complexes and drug-like molecules

Background

Partial atomic charges are a well-established concept, useful in understanding and modeling the chemical behavior of molecules, from simple compounds, to large biomolecular complexes with many reactive sites.

Results

This paper introduces AtomicChargeCalculator (ACC), a web-based application for the calculation and analysis of atomic charges which respond to changes in molecular conformation and chemical environment. ACC relies on an empirical method to rapidly compute atomic charges with accuracy comparable to quantum mechanical approaches. Due to its efficient implementation, ACC can handle any type of molecular system, regardless of size and chemical complexity, from drug-like molecules to biomacromolecular complexes with hundreds of thousands of atoms. ACC writes out atomic charges into common molecular structure files, and offers interactive facilities for statistical analysis and comparison of the results, in both tabular and graphical form.

Conclusions

Due to high customizability and speed, easy streamlining and the unified platform for calculation and analysis, ACC caters to all fields of life sciences, from drug design to nanocarriers. ACC is freely available via the Internet at http://ncbr.muni.cz/ACC.

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-015-0099-x) contains supplementary material, which is available to authorized users.

Direct computation of parameters for accurate polarizable force fields

We present an improved electronic linear response model to incorporate polarization and charge-transfer effects in polarizable force fields. This model is a generalization of the Atom-Condensed Kohn-Sham Density Functional Theory (DFT), approximated to second order (ACKS2): it can now be defined with any underlying variational theory (next to KS-DFT) and it can include atomic multipoles and off-center basis functions. Parameters in this model are computed efficiently as expectation values of an electronic wavefunction, obviating the need for their calibration, regularization, and manual tuning. In the limit of a complete density and potential basis set in the ACKS2 model, the linear response properties of the underlying theory for a given molecular geometry are reproduced exactly. A numerical validation with a test set of 110 molecules shows that very accurate models can already be obtained with fluctuating charges and dipoles. These features greatly facilitate the development of polarizable force fields.

Can the electronegativity equalization method predict spectroscopic properties?

The electronegativity equalization method is classically used as a method allowing the fast generation of atomic charges using a set of calibrated parameters and provided knowledge of the molecular structure. Recently, it has started being used for the calculation of other reactivity descriptors and for the development of polarizable and reactive force fields. For such applications, it is of interest to know whether the method, through the inclusion of the molecular geometry in the Taylor expansion of the energy, would also allow sufficiently accurate predictions of spectroscopic data. In this work, relevant quantities for IR spectroscopy are considered, namely the dipole derivatives and the Cartesian Hessian. Despite careful calibration of parameters for this specific task, it is shown that the current models yield insufficiently accurate results.

An atomistic-based chemophysical environment for evaluating asphalt oxidation and antioxidants

Asphalt binders in service conditions are subject to oxidative aging that involves the reactions between oxygen molecules and the component species of bulk asphalt. As a result, significant alterations can occur to the desired physical and/or mechanical properties of asphalt. A common practice to alleviate asphalt aging has been to employ different chemical additives or modifiers as antioxidants. The current state of knowledge in asphalt oxidation and antioxidant evaluation is centered on determining the degradation of asphalt physical properties, mainly the viscosity and ductility. Such practices, although meeting direct engineering needs, do not contribute to the fundamental understanding of the aging and anti-oxidation mechanisms, and thereby developing anti-aging strategies. From this standpoint, this study was initiated to study the chemical and physical bases of asphalt oxidation, as well as the anti-oxidation mechanisms of bio-based antioxidants using the coniferyl-alcohol lignin as an example. A quantum chemistry (QC) based chemophysical environment is developed, in which the various chemical reactions between asphalt component species and oxygen, as well as the incurred physical changes are studied. X-ray photoelectron spectroscopy (XPS) was used to validate the modified and unmodified asphalt models.

Accurate conformation-dependent molecular electrostatic potentials for high-throughput in silico drug discovery

The atom-centered partial charges-approximation is commonly used in current molecular modeling tools as a computationally inexpensive alternative to quantum mechanics for modeling electrostatics. Even today, the use of partial charges remains useful despite significant advances in improving the efficiency of ab initio methods. Here, we report on new parameters for the EEM and SFKEEM electronegativity equalization-based methods for rapidly determining partial charges that will accurately model the electrostatic potential of flexible molecules. The developed parameters cover most pharmaceutically relevant chemistries, and charges obtained using these parameters reproduce the B3LYP/cc-pVTZ reference electrostatic potential of a set of FDA-approved drug molecules at best to an average accuracy of 13 +/- 4 kJ mol(-1); thus, equipped with these parameters electronegativity equalization-based methods rival the current best non-quantum mechanical methods, such as AM1-BCC, in accuracy, yet incur a lower computational cost. Software implementations of EEM and SFKEEM, including the developed parameters, are included in the conformer-generation tool BALLOON, available free of charge at http://web.abo.fi/fak/mnf/bkf/research/johnson/software.php.

Self-consistent charge equilibration method and its application to Au13Na(n) (n = 1,10) clusters

We modified the charge equilibration method (QEq) of Rappe et al. by including the third and fourth order terms in the series expansion for the energy of charged atoms. This leads to a self-consistent scheme for obtaining charges and total energy in a molecule. We combined the modified QEq with a scaled morse potential (SMP) to get a new potential, SMP/QEq. The SMP/QEq method allows to do realistic geometry optimizations and calculate reliable charge distributions and ionization energies. The performance of SMP/QEq was fully tested by comparing to high level theories and experiments. We established a connection between first-principles DFT and SMP/QEq which gives insight into SMP/QEq and suggests ways to improve it.

A combined electronegativity equalization and electrostatic potential fit method for the determination of atomic point charges

We report an approach for the determination of atomic monopoles of macromolecular systems using connectivity and geometry parameters alone. The method is appropriate also for the calculation of charge distributions based on the quantum mechanically determined wave function and does not suffer from the mathematical instability of other electrostatic potential fit methods.

Optimized and parallelized implementation of the electronegativity equalization method and the atom-bond electronegativity equalization method

The most common way to calculate charge distribution in a molecule is ab initio quantum mechanics (QM). Some faster alternatives to QM have also been developed, the so-called "equalization methods" EEM and ABEEM, which are based on DFT. We have implemented and optimized the EEM and ABEEM methods and created the EEM SOLVER and ABEEM SOLVER programs. It has been found that the most time-consuming part of equalization methods is the reduction of the matrix belonging to the equation system generated by the method. Therefore, for both methods this part was replaced by the parallel algorithm WIRS and implemented within the PVM environment. The parallelized versions of the programs EEM SOLVER and ABEEM SOLVER showed promising results, especially on a single computer with several processors (compact PVM). The implemented programs are available through the Web page http://ncbr.chemi.muni.cz/~n19n/eem_abeem.

Dynamic Charge Equilibration-Morse stretch force field: application to energetics of pure silica zeolites

We present the Dynamic Charge Equilibration (DQEq) method for a self-consistent treatment of charge transfer in force field modeling, where atomic charges are designed to reproduce electrostatic potentials calculated quantum mechanically. Force fields coupled with DQEq allow charges to readjust as geometry changes in classical simulations, using appropriate algorithms for periodic boundary conditions. The full electrostatic energy functional is used to derive the corresponding forces and the second derivatives (hessian) for vibrational calculations. Using DQEq electrostatics, we develop a simple nonbond force field for simulation of silica molecular sieves, where nonelectrostatic interactions are described by two-body Morse stretch terms. Energy minimization calculations with the new force field yield accurate unit cell geometries for siliceous zeolites. Relative enthalpies with respect to quartz and third-law entropies calculated from harmonic vibrational analysis agree very well with available calorimetric data: calculated SiO(2) enthalpies relative to alpha-quartz are within 2 kJ/mol and entropies at 298 K are within 3 J/mol K of the respective experimental values. Contributions from the zero point energy and vibrational degrees of freedom were found to be only about 1 kJ/mol for the free energy of mutual transformations between microporous silica polymorphs. The approach presented here can be applied to interfaces and other oxides as well and it is suitable for development of force fields for accurate modeling of geometry and energetics of microporous and mesoporous materials, while providing a realistic description of electrostatic fields near surfaces and inside pores of adsorbents and catalysts.