Prediction of conformationally dependent atomic multipole moments in carbohydrates

The conformational flexibility of carbohydrates is challenging within the field of computational chemistry. This flexibility causes the electron density to change, which leads to fluctuating atomic multipole moments. Quantum Chemical Topology (QCT) allows for the partitioning of an “atom in a molecule,” thus localizing electron density to finite atomic domains, which permits the unambiguous evaluation of atomic multipole moments. By selecting an ensemble of physically realistic conformers of a chemical system, one evaluates the various multipole moments at defined points in configuration space. The subsequent implementation of the machine learning method kriging delivers the evaluation of an analytical function, which smoothly interpolates between these points. This allows for the prediction of atomic multipole moments at new points in conformational space, not trained for but within prediction range. In this work, we demonstrate that the carbohydrates erythrose and threose are amenable to the above methodology. We investigate how kriging models respond when the training ensemble incorporating multiple energy minima and their environment in conformational space. Additionally, we evaluate the gains in predictive capacity of our models as the size of the training ensemble increases. We believe this approach to be entirely novel within the field of carbohydrates. For a modest training set size of 600, more than 90% of the external test configurations have an error in the total (predicted) electrostatic energy (relative to ab initio) of maximum 1 kJ mol^?1 for open chains and just over 90% an error of maximum 4 kJ mol^?1 for rings. © 2015 Wiley Periodicals, Inc.     

Calculation of electrostatic interaction energies in molecular dimers from atomic multipole moments obtained by different methods of electron density partitioning

Accurate and fast evaluation of electrostatic interactions in molecular systems is still one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry, including molecular recognition, protein modeling and drug design. One of the most convenient and accurate approaches is based on a Buckingham-type approximation that uses the multipole moment expansion of molecular/atomic charge distributions. In the mid-1980s it was shown that the pseudoatom model commonly used in experimental X-ray charge density studies can be easily combined with the Buckingham-type approach for calculation of electrostatic interactions, plus atom-atom potentials for evaluation of the total interaction energies in molecular systems. While many such studies have been reported, little attention has been paid to the accuracy of evaluation of the purely electrostatic interactions as errors may be absorbed in the semiempirical atom-atom potentials that have to be used to account for exchange repulsion and dispersion forces. This study is aimed at the evaluation of the accuracy of the calculation of electrostatic interaction energies with the Buckingham approach. To eliminate experimental uncertainties, the atomic moments are based on theoretical single-molecule electron densities calculated at various levels of theory. The electrostatic interaction energies for a total of 11 dimers of alpha-glycine, N-acetylglycine and L-(+)-lactic acid structures calculated according to Buckingham with pseudoatom, stockholder and atoms-in-molecules moments are compared with those evaluated with the Morokuma-Ziegler energy decomposition scheme. For alpha-glycine a comparison with direct "pixel-by-pixel" integration method, recently developed Gavezzotti, is also made. It is found that the theoretical pseudoatom moments combined with the Buckingham model do predict the correct relative electrostatic interactions energies, although the absolute interaction energies are underestimated in some cases. The good agreement between electrostatic interaction energies computed with Morokuma-Ziegler partitioning, Gavezzotti's method, and the Buckingham approach with atoms-in-molecules moments demonstrates that reliable and accurate evaluation of electrostatic interactions in molecular systems of considerable complexity is now feasible.     

Low cost prediction of relative stabilities of hydrogen bonded complexes from atomic multipole moments for overly short intermolecular distances

The relative stability of biologically relevant, hydrogen bonded complexes with shortened distances can be assessed at low cost by the electrostatic multipole term alone more successfully than by ab initio methods. These results imply that atomic multipole moments may help improve ligand–receptor ranking predictions, particularly in cases where accurate structural data are not available. © 2013 Wiley Periodicals, Inc.     

Point atomic multipole moments for simulation of electrostatic potential and field in all-siliceous zeolites

Calibration method of atomic multipole moments (AMMs) is presented with respect to geometries of all-siliceous zeolite models obtained with X-ray diffraction (XRD) methods. Mulliken atomic charges and AMMs are calculated for all-siliceous types possessing small size elementary unit cells at the hybrid density functional theory (DFT) (B3LYP) and general gradient approximation (GGA) Perdew-Burke-Ernzerhof (PBE) levels and then used to fit the dependences versus geometry variables for the Mulliken charges and versus special coordinate for the AMMs. Fitted and exact charges and AMMs are used to compute electrostatic potential (EP) and electric field (EF) for all-siliceous zeolites with CRYSTAL. A possibility of application of the point AMMs to quantum mechanical/molecular mechanics computations or classic simulation of physical adsorption is evaluated. The considered models expand over wide range of structural parameters and could be applied even to amorphous all-siliceous systems.     

Analysis of nuclear quantum effects on hydrogen bonding

The impact of nuclear quantum effects on hydrogen bonding is investigated for a series of hydrogen fluoride (HF)n clusters and a partially solvated fluoride anion, F-(H2O). The nuclear quantum effects are included using the path integral formalism in conjunction with the Car-Parrinello molecular dynamics (PICPMD) method and using the second-order vibrational perturbation theory (VPT2) approach. For the HF clusters, a directional change in the impact of nuclear quantum effects on the hydrogen-bonding strength is observed as the clusters evolve toward the condensed phase. Specifically, the inclusion of nuclear quantum effects increases the F-F distances for the (HF)n=2-4 clusters and decreases the F-F distances for the (HF)n>4 clusters. This directional change occurs because the enhanced electrostatic interactions between the HF monomers become more dominant than the zero point energy effects of librational modes as the size of the HF clusters increases. For the F-(H2O) system, the inclusion of nuclear quantum effects decreases the F-O distance and strengthens the hydrogen bonding interaction between the fluoride anion and the water molecule because of enhanced electrostatic interactions. The vibrationally averaged 19F shielding constant for F-(H2O) is significantly lower than the value for the equilibrium geometry, indicating that the electronic density on the fluorine decreases as a result of the quantum delocalization of the shared hydrogen. Deuteration of this system leads to an increase in the vibrationally averaged F-O distance and nuclear magnetic shielding constant because of the smaller degree of quantum delocalization for deuterium.     

The electrostatic potential generated by topological atoms. II. Inverse multipole moments

Quantum chemical topology defines finite atoms, whose bounded electron density generates a well-defined electrostatic potential. A multipole expansion based on spherical tensors provides a potential that is formally convergent outside the divergence sphere. Part I of this series [P. L. A. Popelier and M. Rafat, Chem. Phys. Lett.376, 148 (2003)] showed that a continuous multipole expansion expands the convergence region, thereby allowing the electrostatic potential to be evaluated at short range. Here, we propose a different method, based on "inverse" multipole moments, enabling an expansion that converges everywhere. These moments are defined by inverse (i.e., negative) powers of the magnitude of the position vector describing the electron density inside the atom. We illustrate this technique on nitrogen in N(2), oxygen in H(2)O, and oxygen in the phenolic group of the amino acid tyrosine. The proposed method constitutes a considerable advance over the method presented in Part I.     

Fundamental relationships in the theory of electric multipole moments and multipole polarizabilities in static fields

New derivations are given of equations relating molecular electric multipole moments and polarizabilities of general order to the electrostatic energy. The unabridged moment convention is shown to yield relatively simple relations between derivatives of the energy with respect to field gradients and the multipole moments and polarizabilities. Care is taken to distinguish various forms of these derivatives, and one form leads to a proof of a general symmetry of polarizability tensors with respect to permutations of rank indices. The condition of internal equilibrium is shown to be fundamental to the existence of this symmetry. The transformation of multipole moment and polarizability tensors under translation of the coordinate origin is expressed in relatively simple general form. The traceless multipole and polarizability tensors of Buckingham and McLean and Yoshimine are obtained as linear combinations of the unabridged tensors and their traces.     

The effects of dissolved halide anions on hydrogen bonding in liquid water

It is widely believed that the addition of salts to water engenders structural changes in the hydrogen-bond network well beyond the adjacent shell of solvating molecules. Classification of many ions as "structure makers" and "structure breakers" has been based in part on corresponding changes in the vibrational spectra (Raman and IR). Here we show that changes in O-H vibrational spectra induced by the alkali halides in liquid water result instead from the actions of ions' electric fields on adjacent water molecules. Computer simulations that accurately reproduce our experimental measurements suggest that the statistics of hydrogen-bond strengths are only weakly modified beyond this first solvation shell.     

Cumulative coordinates for approximations of high‐order atomic multipole moments in aluminosilicate and aluminophosphate sieves*

A scheme to obtain approximate analytical functions for the atomic distributed multipole moments of the crystallographically different atoms within aluminosilicate and aluminophosphate sieves is discussed. Respective atomic multipole moments are derived within the CRYSTAL95 ab initio periodic Hartree–Fock code with different basis sets, from minimal STO‐3G to 6‐21G*. In order to illustrate the possible applications, distributed analyses are carried out for various structural models from all‐siliceous zeolites and aluminophosphates with ratio Al/P=1 to hydrogen forms of aluminosilicates. Simple approximate forms based on charge and geometry coordinates are proposed for the high‐order moments of each atom, which are further required for the calculation of the electrostatic field within the structures. The possibility to use this analytical approach to evaluate the electrostatic field within embedded cluster models is also shortly discussed. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 83: 70–85, 2001     

Thermoplastic elastomers by hydrogen bonding 4. Influence of hydrogen bonding on the temperature dependence of the viscoelastic properties

The temperature dependence of the viscoelastic properties of thermoreversible polybutadiene networks based on hydrogen bond linkages is analyzed from the logarithmic shift factors loga _T . For binary hydrogen bond complexes thermorheologically simple behavior is observed. The temperature dependence of loga _T is described by the Williams-Landel-Ferry (WLF) equation. The thermoreversible linkages cause an increase in the apparent activation enthalpy of flow which is related to the number of complexing sites in the polymer. Thermorheologically complex behavior is observed in a system with more complex association.     

Electrical effects on the vibrational transitions of hydrogen fluoride due to hydrogen bonding and applied fields

Derivative Hartree-Fock (DHF) theory has been used to study the transition dipole moment of hydrogen fluoride in several hydrogen-bonded complexes and in the presence of applied fields. DHF is an open-ended, analytical means for finding energy derivatives with respect to any number of parameters, and large basis set calculations going through the seventh dipole hyperpolarizability (ninth derivative) are reported. Using multipole moments, multipole polarizabilities and hyperpolarizabilities, the intermolecular electrical influence on vibrational transitions is analyzed.     

Miscible blends through hydrogen bonding: effects on polymer properties

Miscible blends through hydrogen bonding have been intensively studied. The effects of a variety of miscible hydrogen bonded polymer blends on properties such as thermal and thermal oxidative stability, moisture sensitivity, modulus and glass transition temperature are discussed. In addition, the preparation of semi-interpenetrating polymer networks (IPNs) and studies of the effect of crosslinking on the miscibility in hydrogen bonded polymer blends are reviewed.