Force related atomic multipoles in planar molecules. Derivation of atomic quadrupole and octupole moments

Atomic multipoles as defined by current methods generally do not account for forces in molecules that arise from external electrostatic fields. It is pointed out that such forces and the electrostatic potential that the molecule itself generates are both determined by the molecular multipolar tensors. The latter constitute therefore the fundamental molecular constants that determine the molecular electrostatics apart from polarization. In general the multipolar tensors include contributions from the atomic multipoles and their fluxes. In planar molecules, however, the perpendicular charge flux is zero by symmetry. This gives rise to a (previously introduced) formalism that extracts analytical, force-related, atomic multipoles from the molecular multipolar tensors. This formalism is extended in this work to include force-related (FR) atomic quadrupoles and octupoles in planar molecules. The properties of the FR atomic multipoles, including their perpendicular fluxes, are discussed and some formal theoretical and computational advantages that characterize them are indicated. As an example, the electrostatics of OCS, including the molecular electrostatic potential and the forces on the nuclei due to an external point charge, is discussed.     

Atomic multipoles and perpendicular electrostatic forces in diatomic and planar molecules

Current methods for assigning atomic multipoles focus on reproduction of the molecular electrostatic potential. Another aspect of electrostatic interaction, which is usually overlooked, is the forces that an external electric field exerts on the nuclei of a molecule. In a self-consistent theory, both the electrostatic potential and force should be accounted for. However, in general it is not easy to meet this requirement for the force. For planar molecules, though, a formal solution is available in terms of atomic multipoles that are extracted from the molecular multipolar tensors. These Force-Related (FR) atomic multipoles are discussed in detail for some typical diatomics and planar polyatomics, and are shown to provide a solid uniform framework for treating both aspects of the electrostatics. In contrast, the commonly used potential-derived charges (i.e., the atomic charges obtained by fitting the electrostatic potential) can yield large deviations with respect to electrostatic forces on the nuclei, even when the electrostatic potential is very well reproduced.     

Atomic forces for geometry‐dependent point multipole and Gaussian multipole models

In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry‐dependent multipole models. In this study, atomic force expressions for geometry‐dependent multipoles are presented for use in simulations of flexible molecules. The atomic forces are derived by first proposing a new general expression for Wigner function derivatives . The force equations can be applied to electrostatic models based on atomic point multipoles or Gaussian multipole charge density. Hydrogen‐bonded dimers are used to test the intermolecular electrostatic energies and atomic forces calculated by geometry‐dependent multipoles fit to the ab initio electrostatic potential. The electrostatic energies and forces are compared with their reference ab initio values. It is shown that both static and geometry‐dependent multipole models are able to reproduce total molecular forces and torques with respect to ab initio, whereas geometry‐dependent multipoles are needed to reproduce ab initio atomic forces. The expressions for atomic force can be used in simulations of flexible molecules with atomic multipoles. In addition, the results presented in this work should lead to further development of next generation force fields composed of geometry‐dependent multipole models. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Atomic multipoles: Electrostatic potential fit,local reference axis systems,and conformational dependence

Currently, all standard force fields for biomolecular simulations use point charges to model intermolecular electrostatic interactions. This is a fast and simple approach but has deficiencies when the electrostatic potential (ESP) is compared to that from ab initio methods. Here, we show how atomic multipoles can be rigorously implemented into common biomolecular force fields. For this, a comprehensive set of local reference axis systems is introduced, which represents a universal solution for treating atom‐centered multipoles for all small organic molecules and proteins. Furthermore, we introduce a new method for fitting atomic multipole moments to the quantum mechanically derived ESP. This methods yields a 50–90% error reduction compared to both point charges fit to the ESP and multipoles directly calculated from the ab initio electron density. It is shown that it is necessary to directly fit the multipole moments of conformational ensembles to the ESP. Ignoring the conformational dependence or averaging over parameters from different conformations dramatically deteriorates the results obtained with atomic multipole moments, rendering multipoles worse than partial charges. © 2012 Wiley Periodicals, Inc.     

Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides

The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.     

Torque and atomic forces for Cartesian tensor atomic multipoles with an application to crystal unit cell optimization

New equations for torque and atomic force are derived for use in flexible molecule force fields with atomic multipoles. The expressions are based on Cartesian tensors with arbitrary multipole rank. The standard method for rotating Cartesian tensor multipoles and calculating torque is to first represent the tensor with n indexes and 3ⁿ redundant components. In this work, new expressions for directly rotating the unique (n + 1)(n + 2)/2 Cartesian tensor multipole components Θ_pqr are given by introducing Cartesian tensor rotation matrix elements X( R ). A polynomial expression and a recursion relation for X( R ) are derived. For comparison, the analogous rotation matrix for spherical tensor multipoles are the Wigner functions D( R ). The expressions for X( R ) are used to derive simple equations for torque and atomic force. The torque and atomic force equations are applied to the geometry optimization of small molecule crystal unit cells. In addition, a discussion of computational efficiency as a function of increasing multipole rank is given for Cartesian tensors. © 2016 Wiley Periodicals, Inc.     

Dynamic multipole polarizabilities of atomic systems

Two previous schemes for the calculation of dynamic multipole polarizabilities of atomic systems, namely, a hydrodynamic analogy to quantum mechanics and a perturbation expansion of the energy in multipoles, follow also from the conventional time-dependent perturbation scheme.     

Effective polarization energy of the naphthalene molecular crystal: a study on the polarizable force field

polarization energy of the localized charge in organic solids consists of electronic polarization energy, permanent electrostatic interactions, and inter/intra molecular relaxation energies. The effective electronic polarization energies for an electron/hole carrier were successfully estimated by AMOEBA polarizable force field in naphthalene molecular crystals. Both electronic polarization energy and permanent electrostatic interaction were in agreement with the preview experimental values. In addition, the influence of the multipoles from different distributed mutipole analysis (DMA) fitting options on the electrostatic interactions are discussed in this paper. We found that the multipoles obtained from Gauss-Hermite quadrature without diffuse function or grid-based quadrature with 0.325 Å H atomic radius will give reasonable electronic polarization energies and permanent interactions for electron and hole carriers.     

Consistent treatment of inter- and intramolecular polarization in molecular mechanics calculations

A protocol is described for the treatment of molecular polarization in force field calculations. The resulting model is consistent in that both inter- and intramolecular polarization are handled within a single scheme. An analytical formula for removing intramolecular polarization from a set of atomic multipoles for an arbitrary static structure or conformation is given. With the help of the intramolecular polarization, these permanent atomic multipoles can then be applied in modeling alternative conformations of a molecule. Equipped with this simple technique, one can derive transferable electrostatic parameters for peptides and proteins using flexible model compounds such as dipeptides. The proposed procedure is tested for its ability to describe the electrostatic potential around various configurations of the N-methylacetamide dimer. The effect of different intramolecular polarization schemes on the accuracy of a force field model of the electrostatic potential of alanine dipeptide is investigated. A group-based scheme for including direct intramolecular polarization is shown to be most successful in accounting for the conformational dependence of electrostatic potentials.     

First-principles calculation of local atomic polarizabilities

Common methods of determining atomic polarizabilities suffer from the inclusion of nonlocal effects such as charge polarization. A new method is described for determining fully ab initio atomic polarizabilities based on calculating the response of atomic multipoles to the local electrostatic potential. The localized atomic polarizabilities are then used to calculate induction energies that are compared to ab initio induction energies to test their usefulness in practical applications. These polarizabilities are shown to be an improvement over the corresponding molecular polarizabilities, in terms of both absolute accuracy and the convergence of the multipolar induction series. The transferability of localized polarizabilities for the alkane series is also discussed.     

Reconciling alternate methods for the determination of charge distributions: a probabilistic approach to high-dimensional least-squares approximations

We propose extensions and improvements of the statistical analysis of distributed multipoles (SADM) algorithm put forth by Chipot et al. (Mol Phys 94: 881–895, 1998) for the derivation of distributed atomic multipoles from the quantum-mechanical electrostatic potential. The method is mathematically extended to general least-squares problems and provides an alternative approximation method in cases where the original least-squares problem is computationally not tractable, either because of its ill-posedness or its high-dimensionality. The solution is approximated employing a Monte Carlo method that takes the average of a random variable defined as the solutions of random small least-squares problems drawn as subsystems of the original problem. The conditions that ensure convergence and consistency of the method are discussed, along with an analysis of the computational cost in specific instances.     

Atomic partitioning of two-center potentials for slater basis

In a previous paper we have shown that the long-range potential generated by the two-center distributions can be written as a sum of two multipolar expansions centered at the terminal atoms and defined so that the two series of multipoles be minimal. In this paper we show that both charge distributions and short-range potentials can also be separated in atomic contributions in a way consistent with the above mentioned partition of the long-range potentials.     

Microwave assisted micellar extraction coupled with solid phase microextraction for the determination of organochlorine pesticides in soil samples

Microwave assisted micellar extraction (MAME) coupled with solid phase microextraction (SPME) and HPLC-UV determination have been used for the determination of five organochlorine pesticides from agricultural soil samples. A non-ionic surfactant, Polyoxyethlylene 10 Lauryl Ether was used, and the different variables for the optimization of MAME and SPME procedures were studied. This method was applied successfully to the determination of these pesticides in several kinds of agricultural soil samples with different characteristics. Most of the compounds studied can be recovered in good yields with R.S.D. lower than 9% and detection limit ranged between 56-96 ng g⁻¹ for the pesticides studied.     

Cumulative atomic multipole representation of the molecular charge distribution and its basis set dependence

A simple procedure to decompose the theoretical molecular charge distribution into cumulative atomic multipoles supplementing any population analysis scheme has been described and tested for a number of molecules in extended basis sets. This approach may be applied to describe local charge distributions in neutral as well as charged systems and also leads to a simplified point-charge model conserving the local anisotropy of the atomic charge distribution in molecules. Such an approach may be useful in estimating intermolecular interactions, representing the molecular environment in solvent effect or enzyme catalytic activity studies, evaluation of molecular electrostatic potentials or tracing the quality of basis set functions.     

Electrostatic representations of molecular groups

The treatment of electrostatic interactions between molecular groups is based on quantum mechanical electron density distributions. A general method is suggested using tesseral harmonic multipole expansions and representing each set of point multipoles by a cluster of discrete charges small with respect to atomic dimensions. The procedure is carried out for the amino group with an electron density derived from the wave function of ammonia. We evaluate the optical rotation of two aminopyrrolidones as a test case. Expansions are compared for several approximate wave functions for ammonia.     

An alternative parametrization for the theory of complex spectra

The aim of this paper is to investigate a new parametrization for the matrix, set up on the states of an open (nuclear or atomic) shell, of any two-body interaction expandable in multipoles. Such a parametrization is particularly convenient for expressing the term energies in a simple form and for producing high degeneracies among the terms. Further, the parametrization considered allows the two-body interaction to be decomposed into a Dirac interaction and a remaining part.