Towards an atomistic understanding of polymorphism in molecular solids.

Understanding and controlling polymorphism in molecular solids is a major unsolved problem in crystal engineering. While the ability to calculate accurate lattice energies with atomistic modelling provides valuable insight into the associated energy scales, existing methods cannot connect energy differences to the delicate balances of intra- and intermolecular forces that ultimately determine polymorph stability ordering. We report herein a protocol for applying Quantum Chemical Topology (QCT) to study the key intra- and intermolecular interactions in molecular solids, which we use to compare the three known polymorphs of succinic acid including the recently-discovered γ form. QCT provides a rigorous partitioning of the total energy into contributions associated with topological atoms, and a quantitative and chemically intuitive description of the intra- and intermolecular interactions. The newly-proposed Relative Energy Gradient (REG) method ranks atomistic energy terms (steric, electrostatic and exchange) by their importance in constructing the total energy profile for a chemical process. We find that the conformation of the succinic acid molecule is governed by a balance of large and opposing electrostatic interactions, while the H-bond dimerisation is governed by a combination of electrostatics and sterics. In the solids, an atomistic energy balance emerges that governs the contraction, towards the equilibrium geometry, of a molecular cluster representing the bulk crystal. The protocol we put forward is as general as the capabilities of the underlying quantum-mechanical model and it can provide novel perspectives on polymorphism in a wide range of chemical systems.

A FFLUX Water Model: Flexible, Polarizable and with a Multipolar Description of Electrostatics.

Key to progress in molecular simulation is the development of advanced models that go beyond the limitations of traditional force fields that employ a fixed, point charge-based description of electrostatics. Taking water as an example system, the FFLUX framework is shown capable of producing models that are flexible, polarizable and have a multipolar description of the electrostatics. The kriging machine-learning methods used in FFLUX are able to reproduce the intramolecular potential energy surface and multipole moments of a single water molecule with chemical accuracy using as few as 50 training configurations. Molecular dynamics simulations of water clusters (25-216 molecules) using the new FFLUX model reveal that incorporating charge-quadrupole, dipole-dipole, and quadrupole-charge interactions into the description of the electrostatics results in significant changes to the intermolecular structuring of the water molecules. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study of the Halogen Bond with Explicit Analysis of Electron Correlation.

Energy profiles of seven halogen-bonded complexes were analysed with the topological energy partitioning called Interacting Quantum Atoms (IQA) at MP4(SDQ)/6-31+G(2d,2p) level of theory. Explicit interatomic electron correlation energies are included in the analysis. Four complexes combine X2 (X = Cl or F) with HCN or NH3, while the remaining three combine ClF with HCN, NH3 or N2. Each complex was systematically deformed by translating the constituent molecules along its central axis linking X and N, and reoptimising its remaining geometry. The Relative Energy Gradient (REG) method (Theor. Chem. Acc. 2017, 136, 86) then computes which IQA energies most correlate with the total energy during the process of complex formation and further compression beyond the respective equilibrium geometries. It turns out that the covalent energy (i.e., exchange) of the halogen bond, X…N, itself drives the complex formation. When the complexes are compressed from their equilibrium to shorter X…N distance then the intra-atomic energy of N is in charge. When the REG analysis is restricted to electron correlation then the interatomic correlation energy between X and N again drives the complex formation, and the complex compression is best described by the destabilisation of the through-space correlation energy between N and the "outer" halogen.

Atomic Partitioning of the MPn (n = 2, 3, 4) Dynamic Electron Correlation Energy by the Interacting Quantum Atoms Method: A Fast and Accurate Electrostatic Potential Integral Approach.

Recently, the quantum topological energy partitioning method called interacting quantum atoms (IQA) has been extended to MPn (n = 2, 3, 4) wave functions. This enables the extraction of chemical insight related to dynamic electron correlation. The large computational expense of the IQA-MPn approach is compensated by the advantages that IQA offers compared to older nontopological energy decomposition schemes. This expense is problematic in the construction of a machine learning training set to create kriging models for topological atoms. However, the algorithm presented here markedly accelerates the calculation of atomically partitioned electron correlation energies. Then again, the algorithm cannot calculate pairwise interatomic energies because it applies analytical integrals over whole space (rather than over atomic volumes). However, these pairwise energies are not needed in the quantum topological force field FFLUX, which only uses the energy of an atom interacting with all remaining atoms of the system that it is part of. Thus, it is now feasible to generate accurate and sizeable training sets at MPn level of theory. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Infrared Intensification and Hydrogen Bond Stabilization: Beyond Point Charges.

Infrared band intensification of the A-H bond stretching mode of A-H···B acid-base systems has long been known to be the most spectacular spectral change occurring on hydrogen bonding. A QTAIM/CCTDP model is reported here to quantitatively explain the electronic structure origins of intensification and investigate the correlation between experimental enthalpies of formation and infrared hydrogen bond stretching intensifications amply reported in the literature. Augmented correlation-consistent polarized triple-zeta quantum calculations at the MP2 level were performed on complexes with HF and HCl electron acceptors and HF, HCl, NH3, H2O, HCN, acetonitrile, formic acid, acetaldehyde, and formaldehyde electron donor molecules. The A-H stretching band intensities are calculated to be 3 to 40 times larger than their monomer values. Although the acidic hydrogen atomic charge is important for determining the intensities of HF complexes relative to HCl complexes with the same electron donor, they are not important for infrared intensifications occurring on hydrogen bond formation for a series of bases with a common acid. Charge transfers are found to be the most important factor resulting in the intensifications, but dipolar polarization effects are also significant for each series of complexes. A mechanism involving intra-acid and intermolecular electron transfers as well as atomic polarizations is proposed for understanding the intensifications. The calculated sums of the intermolecular electron transfer and acid dipolar polarization contributions to the dipole moment derivatives for each series of complexes are highly correlated with their enthalpies of formation and H-bond intensifications. This could be related to increasing electron transfer from base to acid that correlates with the calculated hydrogen bonding energies and may be a consequence of the A-H bond elongation on complex formation having amplitudes similar to those expected for the A-H vibration.

Decomposition of Møller-Plesset Energies within the Quantum Theory of Atoms-in-Molecules.

We discuss two main approaches to decompose the Møller-Plesset perturbation theory molecular energies into atomic contributions within the interacting quantum atoms (IQA) formalism, as implemented in the programs Morphy and AIMAll. For this purpose, the so-called intraatomic energies (also known as self-energies) of a representative set of 55 small molecules are compared with each other. The origin of the possible discrepancies between both approaches is analyzed, and linear models linking the two approaches are proposed for each atom type.

The Transferability of Topologically Partitioned Electron Correlation Energies in Water Clusters.

The electronic effects that govern the cohesion of water clusters are complex, demanding the inclusion of N-body, Coulomb, exchange and correlation effects. Here we present a much needed quantitative study of the effect of correlation (and hence dispersion) energy on the stabilization of water clusters. For this purpose we used a topological energy partitioning method called Interacting Quantum Atoms (IQA) to partition water clusters into topological atoms, based on a MP2/6-31G(d,p) wave function, and modified versions of GAUSSIAN09 and the Quantum Chemical Topology (QCT) program MORFI. Most of the cohesion in the water clusters provided by electron correlation comes from intramolecular energy stabilization. Hydrogen bond-related interactions tend to largely cancel each other. Electron correlation energies are transferable in almost all instances within 1 kcal mol-1 . This observed transferability is very important to the further development of the QCT force field FFLUX, especially to the future modelling of liquid water.

Quantum Theory of Atoms in Molecules Charge-Charge Transfer-Dipolar Polarization Classification of Infrared Intensities.

Fundamental infrared vibrational transition intensities of gas-phase molecules are sensitive probes of changes in electronic structure accompanying small molecular distortions. Models containing charge, charge transfer, and dipolar polarization effects are necessary for a successful classification of the C-H, C-F, and C-Cl stretching and bending intensities. C-H stretching and in-plane bending vibrations involving sp3 carbon atoms have small equilibrium charge contributions and are accurately modeled by the charge transfer-counterpolarization contribution and its interaction with equilibrium charge movement. Large C-F and C═O stretching intensities have dominant equilibrium charge movement contributions compared to their charge transfer-dipolar polarization ones and are accurately estimated by equilibrium charge and the interaction contribution. The C-F and C-Cl bending modes have charge and charge transfer-dipolar polarization contribution sums that are of similar size but opposite sign to their interaction values resulting in small intensities. Experimental in-plane C-H bends have small average intensities of 12.6 ± 10.4 km mol-1 owing to negligible charge contributions and charge transfer-counterpolarization cancellations, whereas their average out-of-plane experimental intensities are much larger, 65.7 ± 20.0 km mol-1, as charge transfer is zero and only dipolar polarization takes place. The C-F bending intensities have large charge contributions but very small intensities. Their average experimental out-of-plane intensity of 9.9 ± 12.6 km mol-1 arises from the cancellation of large charge contributions by dipolar polarization contributions. The experimental average in-plane C-F bending intensity, 5.8 ± 7.3 km mol-1, is also small owing to charge and charge transfer-counterpolarization sums being canceled by their interaction contributions. Models containing only atomic charges and their fluxes are incapable of describing electronic structure changes for simple molecular distortions that are of interest in classifying infrared intensities. One can expect dipolar polarization effects to also be important for larger distortions of chemical interest.

Atomic polarizations necessary for coherent infrared intensity modeling with theoretical calculations.

The inclusion of atomic polarizations for describing molecular electronic structure changes on vibration is shown to be necessary for coherent infrared intensity modeling. Atomic charges from the ChelpG partition scheme and atomic charges and dipoles from Quantum Theory of Atoms in Molecules (QTAIM) were employed within two different models to describe the stretching and bending vibrational intensities of the C-H, C-F, and C=O groups. The model employing the QTAIM parameters was the Charge-Charge Transfer and Dipolar Polarization model (QTAIM/CCTDP), and the model employing the ChelpG charges was the Equilibrium Charge-Charge Flux (ChelpG/ECCF). The QTAIM/CCTDP models result in characteristic proportions of the charge-charge transfer-dipolar polarization contributions even though their sums giving the total intensities do not discriminate between these vibrations. According to the QTAIM/CCTDP model, the carbon monoxide intensity has electronic structure changes similar to those of the carbonyl stretches whereas they resemble those of the CH stretches for the ChelpG/ECCF model.

Characteristic infrared intensities of carbonyl stretching vibrations.

The experimental infrared fundamental intensities of gas phase carbonyl compounds obtained by the integration of spectral bands in the Pacific Northwest National Laboratory (PNNL) spectral database are in good agreement with the intensities reported by other laboratories having a root mean square error of 27 km mol(-1) or about 13% of the average intensity value. The Quantum Theory of Atoms in Molecules/Charge-Charge Transfer-Counterpolarization (QTAIM/CCTCP) model indicates that the large intensity variation from 61.7 to 415.4 km mol(-1) is largely due to static atomic charge contributions, whereas charge transfer and counterpolarization effects essentially cancel one another leaving only a small net effect. The Characteristic Substituent Shift Model estimates the atomic charge contributions to the carbonyl stretching intensities within 30 km mol(-1) or 10% of the average contribution. However, owing to the size of the 2 × C × CTCP interaction contribution, the total intensities cannot be estimated with this degree of accuracy. The dynamic intensity contributions of the carbon and oxygen atoms account for almost all of the total stretching intensities. These contributions vary over large ranges with the dynamic contributions of carbon being about twice the size of the oxygen ones for a large majority of carbonyls. Although the carbon monoxide molecule has an almost null dipole moment contrary to the very polar bond of the characteristic carbonyl group, its QTAIM/CCTCP model is very similar to those found for the carbonyl compounds.

QTAIM-Based Characteristic Group Infrared Intensities of Amino Acids and Their Transference to Peptides.

Dynamic atomic contributions (DACs) to the infrared intensities of 14 amino acids have been transferred to three peptide molecules, glycylglycine, trialanine, and the melanocyte-inhibiting factor MIF-1, to estimate the infrared intensities of the most strategic peptide bands. The DACs of the amino acids and infrared intensities of the peptides were determined at the DFT B3LYP/6-311+(d,p) level. The Quantum Theory of Atoms In Molecules (QTAIM) Charge-Charge Transfer-Dipolar Polarization (CCTDP) model at this Density Functional Theory (DFT) level was used to classify the O-H, NH2, N-H, and C═O stretching as well as the NH2 bending characteristic groups for use in the transference procedure. Contrary to the frequencies, the intensities within these groups can have very diverse values, although their electronic structure changes upon vibration have predictable QTAIM behaviors for each group. Compared to the DFT calculated values, the two transferred O-H stretching intensities of the peptides are estimated with a root-mean-square (rms) error of 19.1 km mol-1. Six NH2 symmetric and antisymmetric stretching intensities were determined with a 9.9 km mol-1 error. The eight estimated C═O stretching bands have a rms error of 78.0 km mol-1 or 23.6% of the average DFT peptide C═O intensity of 328.4 km mol-1. The proposed procedure is applicable to experimental infrared intensities if a calibration set of molecules with known atomic polar tensors and normal coordinate transformations is available.

Dynamic atomic contributions to infrared intensities of fundamental bands.

Dynamic atomic intensity contributions to fundamental infrared intensities are defined as the scalar products of dipole moment derivative vectors for atomic displacements and the total dipole derivative vector of the normal mode. Intensities of functional group vibrations of the fluorochloromethanes can be estimated within 6.5 km mol(-1) by displacing only the functional group atoms rather than all the atoms in the molecules. The asymmetric CF2 stretching intensity, calculated to be 126.5 km mol(-1) higher than the symmetric one, is accounted for by an 81.7 km mol(-1) difference owing to the carbon atom displacement and 40.6 km mol(-1) for both fluorine displacements. Within the Quantum Theory of Atoms in Molecules (QTAIM) model differences in atomic polarizations are found to be the most important for explaining the difference in these carbon dynamic intensity contributions. Carbon atom displacements almost completely account for the differences in the symmetric and asymmetric CCl2 stretching intensities of dichloromethane, 103.9 of the total calculated value of 105.2 km mol(-1). Contrary to that found for the CF2 vibrations intramolecular charge transfer provoked by the carbon atom displacement almost exclusively explains this difference. The very similar intensity values of the symmetric and asymmetric CH2 stretching intensities in CH2F2 arise from nearly equal carbon and hydrogen atom contributions for these vibrations. All atomic contributions to the intensities for these vibrations in CH2Cl2 are very small. Sums of dynamic contributions of the individual intensities for all vibrational modes of the molecule are shown to be equal to mass weighted atomic effective charges that can be determined from atomic polar tensors evaluated from experimental infrared intensities and frequencies. Dynamic contributions for individual intensities can also be determined solely from experimental data.

An atom in molecules study of infrared intensity enhancements in fundamental donor stretching bands in hydrogen bond formation.

Vibrational modes ascribed to the stretching of X-H bonds from donor monomers (HXdonor) in complexes presenting hydrogen bonds (HF···HF, HCl···HCl, HCN···HCN, HNC···HNC, HCN···HF, HF···HCl and H2O···HF) exhibit large (4 to 7 times) infrared intensity increments during complexation according to CCSD/cc-pVQZ-mod calculations. These intensity increases are explained by the charge-charge flux-dipole flux (CCFDF) model based on multipoles from the Quantum Theory of Atoms in Molecules (QTAIM) as resulting from a reinforcing interaction between two contributions to the dipole moment derivatives with respect to the vibrational displacements: charge and charge flux. As such, variations that occur in their intensity cross terms in hydrogen bond formation correlate nicely with the intensity enhancements. These stretching modes of HXdonor bonds can be approximately modeled by sole displacement of the positively charged hydrogens towards the acceptor terminal atom with concomitant electronic charge transfers in the opposite direction that are larger than those occurring for the H atom displacements of their isolated donor molecules. This analysis indicates that the charge-charge flux interaction reinforcement on H-bond complexation is associated with variations of atomic charge fluxes in both parent molecules and small electronic charge transfers between them. The QTAIM/CCFDF model also indicates that atomic dipole flux contributions do not play a significant role in these intensity enhancements.

Atomic charge transfer-counter polarization effects determine infrared CH intensities of hydrocarbons: a quantum theory of atoms in molecules model.

Atomic charge transfer-counter polarization effects determine most of the infrared fundamental CH intensities of simple hydrocarbons, methane, ethylene, ethane, propyne, cyclopropane and allene. The quantum theory of atoms in molecules/charge-charge flux-dipole flux model predicted the values of 30 CH intensities ranging from 0 to 123 km mol(-1) with a root mean square (rms) error of only 4.2 km mol(-1) without including a specific equilibrium atomic charge term. Sums of the contributions from terms involving charge flux and/or dipole flux averaged 20.3 km mol(-1), about ten times larger than the average charge contribution of 2.0 km mol(-1). The only notable exceptions are the CH stretching and bending intensities of acetylene and two of the propyne vibrations for hydrogens bound to sp hybridized carbon atoms. Calculations were carried out at four quantum levels, MP2/6-311++G(3d,3p), MP2/cc-pVTZ, QCISD/6-311++G(3d,3p) and QCISD/cc-pVTZ. The results calculated at the QCISD level are the most accurate among the four with root mean square errors of 4.7 and 5.0 km mol(-1) for the 6-311++G(3d,3p) and cc-pVTZ basis sets. These values are close to the estimated aggregate experimental error of the hydrocarbon intensities, 4.0 km mol(-1). The atomic charge transfer-counter polarization effect is much larger than the charge effect for the results of all four quantum levels. Charge transfer-counter polarization effects are expected to also be important in vibrations of more polar molecules for which equilibrium charge contributions can be large.

Quantum theory of atoms in molecules/charge-charge flux-dipole flux models for fundamental vibrational intensity changes on H-bond formation of water and hydrogen fluoride.

The Quantum Theory of Atoms In Molecules/Charge-Charge Flux-Dipole Flux (QTAIM/CCFDF) model has been used to investigate the electronic structure variations associated with intensity changes on dimerization for the vibrations of the water and hydrogen fluoride dimers as well as in the water-hydrogen fluoride complex. QCISD/cc-pVTZ wave functions applied in the QTAIM/CCFDF model accurately provide the fundamental band intensities of water and its dimer predicting symmetric and antisymmetric stretching intensity increases for the donor unit of 159 and 47 km mol(-1) on H-bond formation compared with the experimental values of 141 and 53 km mol(-1). The symmetric stretching of the proton donor water in the dimer has intensity contributions parallel and perpendicular to its C2v axis. The largest calculated increase of 107 km mol(-1) is perpendicular to this axis and owes to equilibrium atomic charge displacements on vibration. Charge flux decreases occurring parallel and perpendicular to this axis result in 42 and 40 km mol(-1) total intensity increases for the symmetric and antisymmetric stretches, respectively. These decreases in charge flux result in intensity enhancements because of the interaction contributions to the intensities between charge flux and the other quantities. Even though dipole flux contributions are much smaller than the charge and charge flux ones in both monomer and dimer water they are important for calculating the total intensity values for their stretching vibrations since the charge-charge flux interaction term cancels the charge and charge flux contributions. The QTAIM/CCFDF hydrogen-bonded stretching intensity strengthening of 321 km mol(-1) on HF dimerization and 592 km mol(-1) on HF:H2O complexation can essentially be explained by charge, charge flux and their interaction cross term. Atomic contributions to the intensities are also calculated. The bridge hydrogen atomic contributions alone explain 145, 237, and 574 km mol(-1) of the H-bond stretching intensity enhancements for the water and HF dimers and their heterodimer compared with total increments of 149, 321, and 592 km mol(-1), respectively.

How accessible is atomic charge information from infrared intensities? A QTAIM/CCFDF interpretation.

Infrared fundamental intensities calculated by the quantum theory of atoms in molecules/charge-charge flux-dipole flux (QTAIM/CCFDF) method have been partitioned into charge, charge flux, and dipole flux contributions as well as their charge-charge flux, charge-dipole flux, and charge flux-dipole flux interaction contributions. The interaction contributions can be positive or negative and do not depend on molecular orientations in coordinate systems or normal coordinate phase definitions, as do CCFDF dipole moment derivative contributions. If interactions are positive, their corresponding dipole moment derivative contributions have the same polarity reinforcing the total intensity estimates whereas negative contributions indicate opposite polarities and lower CCFDF intensities. Intensity partitioning is carried out for the normal coordinates of acetylene, ethylene, ethane, all the chlorofluoromethanes, the X(2)CY (X = F, Cl; Y = O, S) molecules, the difluoro- and dichloroethylenes and BF(3). QTAIM/CCFDF calculated intensities with optimized quantum levels agree within 11.3 km mol(-1) of the experimental values. The CH stretching and in-plane bending vibrations are characterized by significant charge flux, dipole flux, and charge flux-dipole flux interaction contributions with the negative interaction tending to cancel the individual contributions resulting in vary small intensity values. CF stretching and bending vibrations have large charge, charge-charge flux, and charge-dipole flux contributions for which the two interaction contributions tend to cancel one another. The experimental CF stretching intensities can be estimated to within 31.7 km mol(-1) or 16.3% by a sum of these three contributions. However, the charge contribution alone is not successful at quantitatively estimating these CF intensities. Although the CCl stretching vibrations have significant charge-charge flux and charge-dipole flux contributions, like those of the CF stretches, both of these interaction contributions have opposite signs for these two types of vibrations.

QTAIM charge-charge flux-dipole flux interpretation of electronegativity and potential models of the fluorochloromethane mean dipole moment derivatives.

Infrared fundamental vibrational intensities and quantum theory atoms in molecules (QTAIM) charge-charge flux-dipole flux (CCFDF) contributions to the polar tensors of the fluorochloromethanes have been calculated at the QCISD/cc-pVTZ level. A root-mean-square error of 20.0 km mol(-1) has been found compared to an experimental error estimate of 14.4 and 21.1 km mol(-1) for MP2/6-311++G(3d,3p) results. The errors in the QCISD polar tensor elements and mean dipole moment derivatives are 0.059 e when compared with the experimental values. Both theoretical levels provide results showing that the dynamical charge and dipole fluxes provide significant contributions to the mean dipole moment derivatives and tend to be of opposite signs canceling one another. Although the experimental mean dipole moment derivative values suggest that all the fluorochloromethane molecules have electronic structures consistent with a simple electronegativity model with transferable atomic charges for their terminal atoms, the QTAIM/CCFDF models confirm this only for the fluoromethanes. Whereas the fluorine atom does not suffer a saturation effect in its capacity to drain electronic charge from carbon atoms that are attached to other fluorine and chlorine atoms, the zero flux electronic charge of the chlorine atom depends on the number and kind of the other substituent atoms. Both the QTAIM carbon charges (r = 0.990) and mean dipole moment derivatives (r = 0.996) are found to obey Siegbahn's potential model for carbon 1s electron ionization energies at the QCISD/cc-pVTZ level. The latter is a consequence of the carbon mean derivatives obeying the electronegativity model and not necessarily to their similarities with atomic charges. Atomic dipole contributions to the neighboring atom electrostatic potentials of the fluorochloromethanes are found to be of comparable size to the atomic charge contributions and increase the accuracy of Siegbahn's model for the QTAIM charge model results. Substitution effects of the hydrogen, fluorine, and chlorine atoms on the charge and dipole flux QTAIM contributions are found to be additive for the mean dipole derivatives of the fluorochloromethanes.

Quantum chemical intensity determinations of overlapped gas phase infrared bands.

The largest source of experimental error in determining gas phase fundamental infrared intensities arises from the separation of overlapped bands. Quantum chemical calculations at the QCISD/cc-pVTZ and QCISD/aug-cc-pVTZ levels were carried out on four simple hydrocarbons and the fluoro- and chloromethanes with the aim of accurate overlapped band separation. Fundamental vibrational intensity results were compared with individual empirical intensity estimates reported for overlapped band systems. Root mean square differences of 3.7 km mol-1 are found between the experimental and QCISD/cc-pVTZ values for nine overlapped bands of the hydrocarbons and 11.8 km mol-1 for the QCISD/aug-cc-pVTZ values for 12 overlapped bands of the fluoro- and chloromethanes. These values correspond to 14% and 18% of the average hydrocarbon and halomethane intensity values. Previous experimental separation errors were estimated to be quite larger, between 20% and 50%. As quantum calculations are continuously being refined one can expect more accurate band separation results in the future.

MP2-IQA: upscaling the analysis of topologically partitioned electron correlation.

When electronic correlation energy is partitioned topologically, a detailed picture of its distribution emerges, both within atoms and between any two atoms. This methodology allows one to study dispersion beyond its more narrow definition in long-range Rayleigh-Schrödinger perturbation theory. The interacting quantum atoms (IQA) method was applied to MP2/6-31G(d,p) (uncontracted) wave functions of a wide variety of systems: glycine…water (hydration), the ethene dimer (π-π interactions), benzene (aromaticity), cyclobutadiene (antiaromaticity), and NH3BH3 (dative bond). Through the study of molecular complexes it turns out that dispersion energy is either important to a system's stabilization (for the C2H4 dimer) or not important (for Gly…H2O). We have also discovered that the delocalization in benzene lowers the strength of Coulomb repulsion in the bonds, which has been quantified for the first time through IQA. Finally, we showed that the nature of the dative bond is much different from that of a regular covalent bond as it is not destabilized by electronic correlation. Finally, the conclusions obtained for these archetypical systems have implications for the future of the quantum topological force field FFLUX in the simulation of larger systems. Graphical abstract Atomic and bond dynamic correlation energies are now available thanks to IQA. Larger molecules can now be accessed to include resonance and solvation of FFLUX force field.

Probing the robustness of the charge-charge transfer-dipolar polarization model and infrared intensities.

The robustness of the QTAIM charge-charge transfer-dipolar polarization parameters for the CH, CF, and CCl stretching and bending distortions of the fluoro- and chloromethanes was determined comparing results calculated at three quantum levels, MP2/6-311G++(3d,3p), QCISD/cc-pVTZ, and QCISD/aug-cc-pVTZ. The correlation coefficients between the MP2/6-311G++G(d,p) and QCISD/cc-pVTZ results with those of QCISD/aug-cc-pVTZ intensities are excellent, 0.934 and 0.988, respectively, showing that the parameters converge with increasing quality of the quantum levels. In spite of numerical differences, the interpretation of the electronic structure changes occurring for these vibrations is the same for all three quantum levels. Accurate determination of charge transfer-counterpolarization effects is important for properly describing electron density changes for small molecular distortions.