Absolutely Localized Molecular Orbitals Research Articles

Quantum Mechanics Characterization of Non-Covalent Interaction in Nucleotide Fragments

Accurate calculation of non-covalent interaction energies in nucleotides is crucial for understanding the driving forces governing nucleic acid structure and function, as well as developing advanced molecular mechanics forcefields or machine learning potentials tailored to nucleic acids. Here, we dissect the nucleotides’ structure into three main constituents: nucleobases (A, G, C, T, and U), sugar moieties (ribose and deoxyribose), and phosphate group. The interactions among these fragments and between fragments and water were analyzed. Different quantum mechanical methods were compared for their accuracy in capturing the interaction energy. The non-covalent interaction energy was decomposed into electrostatics, exchange-repulsion, dispersion, and induction using two ab initio methods: Symmetry-Adapted Perturbation Theory (SAPT) and Absolutely Localized Molecular Orbitals (ALMO). These calculations provide a benchmark for different QM methods, in addition to providing a valuable understanding of the roles of various intermolecular forces in hydrogen bonding and aromatic stacking. With SAPT, a higher theory level and/or larger basis set did not necessarily give more accuracy. It is hard to know which combination would be best for a given system. In contrast, ALMO EDA2 did not show dependence on theory level or basis set; additionally, it is faster.

A Free Energy Decomposition Analysis: Insight into Binding Thermodynamics from Absolutely Localized Molecular Orbitals.

Traditional energy decomposition analysis (EDA) methods can provide an interpretive decomposition of non-covalent electronic binding energies. However, by definition, they neglect entropic effects and nuclear contributions to the enthalpy. With the objective of revealing the chemical origins of trends in free energies of binding, we introduce the concept of a Gibbs decomposition analysis (GDA) by coupling the absolutely localized molecular orbital treatment of electrons in non-covalent interactions with the simplest possible quantum rigid rotor-harmonic oscillator treatment of nuclear motion at finite temperature. The resulting pilot GDA is employed to decompose enthalpic and entropic contributions to the free energy of association of the water dimer, fluoride-water dimer, and water binding to an open metal site in the metal-organic framework Cu(I)-MFU-4l. The results show enthalpic trends that generally track the electronic binding energy and entropic trends that reveal the increasing price of loss of translational and rotational degrees freedom with temperature.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Computational investigation of substituent effects on the fluorescence wavelengths of oxyluciferin analogs

Oxyluciferin, which is the light emitter for firefly bioluminescence, has been subjected to extensive chemical modifications to tune its emission wavelength and quantum yield. However, the exact mechanisms for various electron-donating and withdrawing groups to perturb the photophysical properties of oxyluciferin analogs are still not fully understood. To elucidate the substituent effects on the fluorescence wavelength of oxyluciferin analogs, we applied the absolutely localized molecular orbitals (ALMO)-based frontier orbital analysis to assess various types of interactions (i.e. permanent electrostatics/exchange repulsion, polarization, occupied–occupied orbital mixing, virtual–virtual orbital mixing, and charge-transfer) between the oxyluciferin and substituent orbitals. We suggested two distinct mechanisms that can lead to red-shifted oxyluciferin emission wavelength, a design objective that can help increase the tissue penetration of bioluminescence emission. Within the first mechanism, an electron-donating group (such as an amino or dimethylamino group) can contribute its highest occupied molecular orbital (HOMO) to an out-of-phase combination with oxyluciferin’s HOMO, thus raising the HOMO energy of the substituted analog and narrowing its HOMO-LUMO gap. Alternatively, an electron-withdrawing group (such as a nitro or cyano group) can participate in an in-phase virtual–virtual orbital mixing of fragment LUMOs, thus lowering the LUMO energy of the substituted analog. Such an ALMO-based frontier orbital analysis is expected to lead to intuitive principles for designing analogs of not only the oxyluciferin molecule, but also many other functional dyes.

Localized Spin Rotations: A Size-Consistent Approach to Nonorthogonal Configuration Interaction.

Current nonorthogonal configuration interaction (NOCI) methods often use a set of self-consistent field (SCF) states selected based on chemical intuition. However, it may be challenging to track these SCF states across a dissociation profile and the NOCI states recovered may be spin contaminated. In this Article, we propose a method of applying spin rotation on symmetry broken unrestricted Hartree-Fock (sb-UHF) states to generate a basis for NOCI. The dissociation of ethene was examined by localizing spin rotation on each resulting carbene fragment. We show that this gives a size-consistent description of its dissociation and results in spin-pure states at all geometries. The dissociation was also studied with different orbitals, namely, canonical UHF and absolutely localized molecular orbitals (ALMO). Furthermore, we demonstrate that the method can be used to restore spin symmetry of symmetry broken SCF wave functions for molecules of various sizes, marking an improvement over existing NOCI methods.

Consistent inclusion of continuum solvation in energy decomposition analysis: theory and application to molecular CO2 reduction catalysts.

To facilitate computational investigation of intermolecular interactions in the solution phase, we report the development of ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework of energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, all the quantum mechanical states involved in the variational EDA procedure are computed with the presence of solvent environment so that solvation effects are incorporated in the evaluation of all its energy components. After validation on several model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexes that are related to molecular CO2 reduction catalysis. For [FeTPP(CO2-κC)]2− (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO2 binding, –N(CH3)3+ (TMA) and –OH, stabilize the complex via through-structure and through-space mechanisms, respectively. The coulombic interaction between the positively charged TMA group and activated CO2 is found to be largely attenuated by the polar solvent. Furthermore, we also provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO2via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible “pockets” in that electrostatics is unscreened. For the reactant and product complexes associated with the electron transfer from the p-terphenyl radical anion to CO2, we demonstrate that the double terminal substitution of p-terphenyl by electron-withdrawing groups considerably strengthens the binding in the product state while moderately weakens that in the reactant state, which are both dominated by the substituent tuning of the electrostatics component. These applications illustrate that this new extension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantify their impacts on chemical reactivity in solution.

Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer.

Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time-dependent density functional theory, and exorbitant computational cost when combined with high-level wavefunction-based methods. Here, we introduce a density functional theory (DFT)-based diabatization scheme that directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs), which we denote as Δ-ALMO(MSDFT2). We demonstrate that our method, which combines ALMO calculations with the ΔSCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.

Accurate and efficient DFT-based diabatization for hole and electron transfer using absolutely localized molecular orbitals.

Diabatic states and the couplings between them are important for quantifying, elucidating, and predicting the rates and mechanisms of many chemical and biochemical processes. Here, we propose and investigate approaches to accurately compute diabatic couplings from density functional theory (DFT) using absolutely localized molecular orbitals (ALMOs). ALMOs provide an appealing approach to generate variationally optimized diabatic states and obtain their associated forces, which allows for the relaxation of the donor and acceptor orbitals in a way that is internally consistent in how the method treats both the donor and acceptor states. Here, we show that one can obtain more accurate electronic couplings between ALMO-based diabats by employing the symmetrized transition density matrix to evaluate the exchange-correlation contribution. We demonstrate that this approach yields accurate results in comparison to other commonly used DFT-based diabatization methods across a wide array of electron and hole transfer processes occurring in systems ranging from conjugated organic molecules, such as thiophene and pentacene, to DNA base pairs. We also show that this approach yields accurate diabatic couplings even when combined with lower tiers of the DFT hierarchy, opening the door to combining it with quantum dynamics approaches to provide an ab initio treatment of nonadiabatic processes in the condensed phase.

Silver clusters tune up electronic properties of graphene nanoflakes: A comprehensive theoretical study

Graphene-silver composites show a synergic effect of both components with enhanced electronic, adsorption and catalytic properties. We present here the detailed theoretical study on geometric, thermodynamic and electronic properties of silver clusters Agn (n = 2–10) adsorbed on graphene nanoflakes as a model for silver-graphene composites. First, a benchmark study is performed on Agn…coronene complexes (n = 2, 4) by evaluating various functional against reported high level ab-initio method for the selection of best functional for further calculations. Then, density functional theory study at the best chosen level is performed to inquire how clustering of adsorbate attached with carbon surface can impact geometric, electronic and thermodynamic properties of silver-graphene composites. Furthermore, the interaction energies are corrected for dispersion and basis set superposition error in order to get more realistic estimate of interaction energies. In progression, we investigated the dependence of interaction energy upon cluster size and dimensionality. For Agn…coronene, the binding energies increase monotonically as the number of atoms in the cluster increases. However, for Agn…dodecabenzocoronene, the interaction energies show a sudden drop at Ag7…dodecabenzocoronene. The difference in behavior between Agn…dodecabenzocoronene and Agn…coronene for interaction energy is attributed to edge effect present in coronene. Furthermore, coplanar orientations of the silver clusters on polyaromatic hydrocarbons have higher interaction energies. The 2D-3D transition of the silver cluster is observed for Ag7…coronene. CDA analysis reveals that backdonation from dodecabenzocoronene to silver clusters is more dominant in Ag4–Ag7 whereas other clusters have dominant charge transfer to dodecabenzocoronene. Charge transfer analysis via constrained density functional theory (CDFT) reveals significant energy lowering from charge donor-acceptor interactions. Further energy decomposition analysis (EDA) using absolutely-localized molecular orbitals (ALMO) suggests that, while charge-transfer is important in these systems, the sum of electrostatics, Pauli repulsion, and London dispersion are far more significant as silver cluster size increases. NCI results also support the presence of van der Waals and electrostatic interactions. The hollow top is most favorite interaction sites compared to bridge top and head top binding sites over coronene lattice. The ionization potential, electron affinity, frontier molecular orbital analysis, chemical hardness, softness, chemical potential, and Fermi levels are deliberated to verify the stability and reactivity of most stable silver-graphene composites. The remarkable outcome of the current findings makes our silver-graphene composite as a potential applicant for high performance electronic and catalytic devices.

A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach.

A novel method for the accurate and efficient calculation of interaction energies in weakly bound complexes composed of a large number of molecules is presented. The new ALMO+RPAd method circumvents the prohibitive scaling of coupled cluster singles and doubles while still providing similar accuracy across a diverse range of weakly bound chemical systems. Linear-scaling procedures for the Fock build are given utilizing absolutely localized molecular orbitals (ALMOs), resulting in the a priori exclusion of basis set superposition errors. A bespoke data structure and algorithm using density fitting are described, leading to linear scaling for the storage and computation of the two-electron integrals. Electron correlation is included through a new, linear-scaling pairwise local random phase approximation approach, including exchange interactions, and decomposed into purely dispersive excitations (RPAxd). Collectively, these allow meaningful decomposition of the interaction energy into physically distinct contributions: electrostatic, polarization, charge transfer, and dispersion. Comparison with symmetry-adapted perturbation theory shows good qualitative agreement. Tests on various dimers and the S66 benchmark set demonstrate results within 0.5 kcal mol-1 of coupled cluster singles and doubles results. On a large cluster of water molecules, we achieve calculations involving over 3500 orbital and 12,000 auxiliary basis functions in under 10 min on a single CPU core.

Energy Decomposition Analysis for Excimers Using Absolutely Localized Molecular Orbitals within Time-Dependent Density Functional Theory and Configuration Interaction with Single Excitations.

We present an improved energy decomposition analysis (EDA) scheme for understanding intermolecular interactions in delocalized excited states, especially in excimers. In the EDA procedure, excited states are treated with linear response theory such as configuration interaction singles (CIS) or time-dependent density functional theory (TDDFT), and absolutely localized molecular orbitals (ALMOs) are used to define the intermediate (frozen, excitonic coupling, and polarized) states. The intermolecular interaction energy is thereby separated into frozen, excitonic splitting, polarization, and charge transfer contributions. The excitonic splitting term describes the delocalization effect as two or more degenerate local excitations coupled with each other, which is often an important binding force in excimers. A maximum overlap state-tracking procedure is introduced to connect the initial fragment excitations to the constrained intermediate states and finally to the unconstrained delocalized states of the complex. The EDA scheme is applied to several excimer systems, including the He2* and Ne2* noble gas excimers, the doubly hydrogen-bonded 2-pyridone dimer, and the aromatic benzene and perylene excimers. We are able to gain some useful insights into the role each term is playing in the formation of these excimers, and the resulting method may also be useful for understanding a range of other complexes in excited states.

Efficient Implementation of NOCI-MP2 Using the Resolution of the Identity Approximation with Application to Charged Dimers and Long C-C Bonds in Ethane Derivatives.

An efficient implementation of the perturb-then-diagonalize nonorthogonal configuration interaction method with second-order Møller-Plesset perturbation theory (NOCI-MP2) is presented. Relative to other low scaling multireference perturbation theories, NOCI-MP2 often requires a much smaller active space because of the use of nonorthogonal reference configurations. Reworking the NOCI-MP2 equations with the resolution of the identity (RI) approximation enables the method to have the same memory requirements and computational scaling as single reference RI-MP2. The working equations are extended to include single substitutions as required when the reference determinants do not satisfy the Hartree-Fock equations. A detailed computational algorithm is presented along with timings to establish the performance of the implementation. NOCI-MP2 is applied to the binding energy and charge resonance energy in dication and monocation π dimers, as well as didiamantane ethane, and hexaphenylethane. A well-defined set of nonorthogonal determinants are obtained using absolutely localized molecular orbitals (ALMOs), as solutions to the self-consistent field for molecular interactions (SCF-MI) equations corresponding to covalent and ionic determinants. Agreement with experimental information where available, and other multireference methods, is satisfactory, with the use of an 0.3 au level shift to guard against large MP2 amplitudes. For didiamantane ethane and hexaphenylethane, large dispersion forces help stabilize the molecules despite the steric repulsion. By contrast, in the case of hexaphenylethane, the energy penalty from the geometric distortion of the fragments significantly weakens the bond.

Understanding Non-Covalent Interactions: Correlated Energy Decomposition Analysis and Applications to Halogen Bonding

Non-covalent interactions play a primordial role in chemistry. Beyond their quantification, the detailed understanding of their physical processes is necessary to rationalize chemical trends and improve designs of chemical systems. Energy decomposition analyses allow detailed insight into non-covalent interactions by extracting electrostatics, Pauli repulsion, polarization, dispersion and charge transfer components from interaction energies. Recent work has demonstrated that electronic correlation influenced significantly all of these energy components, whereas previous decompositions only partitioned correlation between dispersion and charge transfer. The MP2 energy decomposition analysis with Absolutely Localized Molecular Orbitals (MP2 ALMO-EDA) takes these results fully into account and offers a correlation correction for each extracted component. A recent detailed investigation of the CCSD dispersion energy showed that a small number of virtual orbitals is sufficient to describe dispersion interactions accurately in the long-range, which potentially offers a basis-set independent definition of dispersion. Finally, we present an application of MP2 ALMO-EDA to a series of unusual halogen bonding complexes where charge transfer dominates over the electrostatic σ-hole interaction.

On the Computational Characterization of Charge-Transfer Effects in Noncovalently Bound Molecular Complexes.

Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a "CT-free" state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the "CT-free" state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal-carbonyl complexes, and complexes formed by water and metal cations. For most of these systems, CDFT yields a much smaller equilibrium CT energy compared to that given by the ALMO-based definition. This is mainly because the CDFT population constraint does not fully inhibit CT, which means that the CDFT "CT-free" state is in fact CT-contaminated. Examples of this contamination include (i) matching forward and backward donation (e.g., formic acid dimer) and (ii) unidirectional CT without changing fragment populations. The magnitude of the latter effect is quantified in systems such as the water dimer by employing a 3-space density constraint in addition to the orbital constraint. Furthermore, by means of the adiabatic EDA, it is shown that several observable effects of CT, such as the "pyramidalization" of the planar BH3 molecule upon the complexation with Lewis bases, already appear on the "CT-free" CDFT surface. These results reveal the essential distinctions between the ALMO and CDFT definitions of CT and suggest that the former is more consistent with accepted understanding of the role of CT in intermolecular binding.

Unraveling substituent effects on frontier orbitals of conjugated molecules using an absolutely localized molecular orbital based analysis.

It is common to introduce electron-donating or electron-withdrawing substituent groups into functional conjugated molecules (such as dyes) to tune their electronic structure properties (such as frontier orbital energy levels) and photophysical properties (such as absorption and emission wavelengths). However, there lacks a generally applicable tool that can unravel the underlying interactions between orbitals from a substrate molecule and those from its substituents in modern electronic structure calculations, despite the long history of qualitative molecular orbital theory. In this work, the absolutely localized molecular orbitals (ALMO) based analysis is extended to analyze the effects of substituent groups on the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of a given system. This provides a bottom-up avenue towards quantification of effects from distinct physical origins (e.g. permanent electrostatics/Pauli repulsion, mutual polarization, inter-fragment orbital mixing). For the example case of prodan (a typical dye molecule), it is found that inter-fragment orbital mixing plays a key role in narrowing the HOMO-LUMO gap of the naphthalene core. Specifically, an out-of-phase mixing of high-lying occupied orbitals on the naphthalene core and the dimethylamino group leads to an elevated HOMO, whereas an in-phase combination of LUMOs on the naphthalene core and the propionyl group lowers the LUMO energy of the entire molecule. We expect this ALMO-based analysis to bridge the gap between concepts from qualitative orbital interaction analysis and quantitative electronic structure calculations.

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis.

The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH3, NH4+, CN-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.

Simulating the absorption spectra of helium clusters (N = 70, 150, 231, 300) using a charge transfer correction to superposition of fragment single excitations.

Simulations of the n = 2 absorption spectra of HeN (N = 70, 150, 231, 300) clusters are reported, with nuclear configurations sampled by path integral molecular dynamics. The electronic structure is treated by a new approach, ALMO-CIS+CT, which is a formulation of configuration interaction singles (CIS) based on absolutely localized molecular orbitals (ALMOs). The method generalizes the previously reported ALMO-CIS model [K. D. Closser et al. J. Chem. Theory Comput. 11, 5791 (2015)] to include spatially localized charge transfer (CT) effects. It is designed to recover large numbers of excited states in atomic and molecular clusters, such as the entire n = 2 Rydberg band in helium clusters. ALMO-CIS+CT is shown to recover most of the error caused by neglecting charge transfer in ALMO-CIS and has comparable accuracy to standard CIS for helium clusters. For the n = 2 band, CT stabilizes states towards the blue edge by up to 0.5 eV. ALMO-CIS+CT retains the formal cubic scaling of ALMO-CIS with respect to system size. With improvements to the implementation over that originally reported for ALMO-CIS, ALMO-CIS+CT is able to treat helium clusters with hundreds of atoms using modest computing resources. A detailed simulation of the absorption spectra associated with the 2s and 2p bands of helium clusters up to 300 atoms is reported, using path integral molecular dynamics with a spherical boundary condition to generate atomic configurations at 3 K. The main features of experimentally reported fluorescence excitation spectra for helium clusters are reproduced.

Efficient Implementation of Energy Decomposition Analysis for Second-Order Møller-Plesset Perturbation Theory and Application to Anion-π Interactions.

Energy decomposition analysis (EDA) is a widely used tool for extracting physical and chemical insights from electronic structure calculations of intermolecular interactions, as well as for the development of advanced force fields for describing those interactions. Recently, the absolutely localized molecular orbital (ALMO) EDA has been extended from the self-consistent field level to the second-order Møller-Plesset (MP2) theory level. This paper reports an efficient implementation of the MP2 ALMO-EDA that scales optimally, employs the resolution of the identity (RI) approximation for post-SCF matrix elements, and is shared-memory parallel. The algorithms necessary to achieve this implementation are described in detail. Performance tests using the aug-cc-pVTZ basis set for water clusters of up to 10 molecules are reported. The timings suggest that the MP2 ALMO-EDA is computationally feasible whenever MP2 energy calculations themselves are feasible, and the cost is dominated by the SCF itself in this size regime. The MP2 ALMO-EDA is applied to study the origin of substituent effects in anion-π interactions between chloride and benzene and mono- through hexafluorobenzene. The effect of fluoro substituents was primarily to change the frozen interaction. Detailed analysis supports the interpretation that anion-π interactions are favorable because of electrostatic interaction with the substituents.

Energy decomposition analysis in an adiabatic picture.

Energy decomposition analysis (EDA) of electronic structure calculations has facilitated quantitative understanding of diverse intermolecular interactions. Nevertheless, such analyses are usually performed at a single geometry and thus decompose a "single-point" interaction energy. As a result, the influence of the physically meaningful EDA components on the molecular structure and other properties are not directly obtained. To address this gap, the absolutely localized molecular orbital (ALMO)-EDA is reformulated in an adiabatic picture, where the frozen, polarization, and charge transfer energy contributions are defined as energy differences between the stationary points on different potential energy surfaces (PESs), which are accessed by geometry optimizations at the frozen, polarized and fully relaxed levels of density functional theory (DFT). Other molecular properties such as vibrational frequencies can thus be obtained at the stationary points on each PES. We apply the adiabatic ALMO-EDA to different configurations of the water dimer, the water-Cl- and water-Mg2+/Ca2+ complexes, metallocenes (Fe2+, Ni2+, Cu2+, Zn2+), and the ammonia-borane complex. This method appears to be very useful for unraveling how physical effects such as polarization and charge transfer modulate changes in molecular properties induced by intermolecular interactions. As an example of the insight obtained, we find that a linear hydrogen bond geometry for the water dimer is preferred even without the presence of polarization and charge transfer, while the red shift in the OH stretch frequency is primarily a charge transfer effect; by contrast, a near-linear geometry for the water-chloride hydrogen bond is achieved only when charge transfer is allowed.

Assessing Ion-Water Interactions in the AMOEBA Force Field Using Energy Decomposition Analysis of Electronic Structure Calculations.

AMOEBA is a molecular mechanics force field that addresses some of the shortcomings of a fixed partial charge model, by including permanent atomic point multipoles through quadrupoles, as well as many-body polarization through the use of point inducible dipoles. In this work, we investigate how well AMOEBA formulates its non-bonded interactions, and how it implicitly incorporates quantum mechanical effects such as charge penetration (CP) and charge transfer (CT), for water-water and water-ion interactions. We find that AMOEBA's total interaction energies, as a function of distance and over angular scans for the water dimer and for a range of water-monovalent cations, agree well with an advanced density functional theory (DFT) model, whereas the water-halides and water-divalent cations show significant disagreement with the DFT result, especially in the compressed region when the two fragments overlap. We use a second-generation energy decomposition analysis (EDA) scheme based on absolutely localized molecular orbitals (ALMOs) to show that in the best cases AMOEBA relies on cancellation of errors by softening of the van der Waals (vdW) wall to balance permanent electrostatics that are too unfavorable, thereby compensating for the missing CP effect. CT, as another important stabilizing effect not explicitly taken into account in AMOEBA, is also found to be incorporated by the softened vdW interaction. For the water-halides and water-divalent cations, this compensatory approach is not as well executed by AMOEBA over all distances and angles, wherein permanent electrostatics remains too unfavorable and polarization is overdamped in the former while overestimated in the latter. We conclude that the DFT-based EDA approach can help refine a next-generation AMOEBA model that either realizes a better cancellation of errors for problematic cases like those illustrated here, or serves to guide the parametrization of explicit functional forms for short-range contributions from CP and/or CT.