HFLD: A Nonempirical London Dispersion-Corrected Hartree-Fock Method for the Quantification and Analysis of Noncovalent Interaction Energies of Large Molecular Systems †

Structural Flexibility is a Decisive Factor in FLP Dihydrogen Cleavage with Tetrahedral Lewis Acids: A Silane Case Study

Dihydrogen activation is the paradigmatic reaction of frustrated Lewis pairs (FLPs). While trigonal-planar Lewis acids have been well established in this transformation, tetrahedral Lewis acids are surprisingly limited. Indeed, several cases were computed as thermodynamically and kinetically feasible but exhibit puzzling discrepancies with experimental results. In the present study, a computational investigation of the factors influencing dihydrogen activation are considered by large ensemble sampling of encounter complexes, deformation energies and the activation strain model for a silicon/nitrogen FLP and compared with a boron/phosphorous FLP. The analysis adds the previously missing dimension of Lewis acids' structural flexibility as a factor that influences preexponential terms beyond pure transition state energies. It sheds light on the origin of "overfrustration" (defined herein), indicates structural constraint in Lewis acids as a linchpin for activation of weak donor substrates, and allows drawing a more refined mechanistic picture of this emblematic reactivity.

Local Energy Decomposition Analysis of London Dispersion Effects: From Simple Model Dimers to Complex Biomolecular Assemblies

ConspectusLondon dispersion (LD) forces are ubiquitous in chemistry, playing a pivotal role in a wide range of chemical processes. For example, they influence the structure of molecular crystals, the selectivity of organocatalytic transformations, and the formation of biomolecular assemblies. Harnessing these forces for chemical applications requires consistent quantification of the LD energy across a broad and diverse spectrum of chemical scenarios. Despite the great progress made in recent years in the development of experimental strategies for LD quantification, quantum chemical methods remain one of the most useful tools in the hand of chemists for the study of these weak interactions. Unfortunately, the accurate quantification of LD effects in complex systems poses many challenges for electronic structure theories. One of the problems stems from the fact that LD forces originate from long-range electronic dynamic correlation, and hence, their rigorous description requires the use of complex, highly correlated wave function-based methods. These methods typically feature a steep scaling with the system size, limiting their applicability to small model systems. Another core challenge lies in disentangling short-range from long-range dynamic correlation, which from a rigorous quantum mechanical perspective is not possible.In this Account, we describe our research endeavors in the development of broadly applicable computational methods for LD quantification in molecular chemistry as well as challenging applications of these schemes in various domains of chemical research. Our strategy lies in the use of local correlation theories to reduce the computational cost associated with complex electronic structure methods while providing at the same time a simple means of decomposition of dynamic correlation into its long-range and short-range components. In particular, the local energy decomposition (LED) scheme at the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) level has emerged as a powerful tool in our research, offering a clear-cut quantitative definition of the LD energy that remains valid across a plethora of different chemical scenarios. Typical applications of this scheme are examined, encompassing protein-ligand interactions and reactivity studies involving many fragments and complex electronic structures. In addition, our research also involves the development of novel cost-effective methodologies, which exploit the LED definition of the LD energy, for LD energy quantification that are, in principle, applicable to systems with thousands of atoms. The Hartree-Fock plus London Dispersion (HFLD) scheme, correcting the HF interaction energy using an approximate CCSD(T)-based LD energy, is a useful, parameter-free electronic structure method for the study of LD effects in systems with hundreds of molecular fragments. However, the usefulness of the LED scheme reaches beyond providing an interpretation of the calculated DLPNO-CCSD(T) or DLPNO-MP2 interaction energies. For example, the dispersion energies obtained from the LED can be fruitfully used in order to parametrize semiempirical dispersion models. We will demonstrate this in the context of an emerging semiempirical method, namely, the Natural Orbital Tied Constructed Hamiltonian (NOTCH) method. NOTCH incorporates LED-derived LD energies and shows promising accuracy at a minimum amount of empiricism. Thus, it holds substantial promise for large and complex systems.

Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

London dispersion (LD) interactions are the main contribution of the attractive part of the van der Waals potential. Even though LD effects are the driving force for molecular aggregation and recognition, the role of these omnipresent interactions in structure and reactivity had been largely underappreciated over decades. However, in the recent years considerable efforts have been made to thoroughly study LD interactions and their potential as a chemical design element for structures and catalysis. This was made possible through a fruitful interplay of theory and experiment. This review highlights recent results and advances in utilizing LD interactions as a structural motif to understand and utilize intra- and intermolecularly LD-stabilized systems. Additionally, we focus on the quantification of LD interactions and their fundamental role in chemical reactions.

Complexes of resorcin[4]arene with secondary amines: synthesis, solvent influence on "in-out" structure, and theoretical calculations of non-covalent interactions

Resorcin[4]arenes (R[4]A) are macrocyclic compounds with a cavity structure. Despite a relatively small cavity, these compounds are capable of forming complexes with small organic molecules. The current paper focuses on the synthesis of complexes between R[4]A and secondary aliphatic amines (sec-amines). Through NMR spectroscopy, it was observed that "in-out" complexes are formed depending on the solvent. It was also found that the stoichiometry of the formed complexes depends on the size of the amine molecule. The automated interaction sites screening (aISS) made it possible to generate molecular ensembles of complexes. The geometry of the ensembles was first optimized with the r2scan-3c functional and, finally, the structure with the lowest energy, with the functional PBE0-D4/mTZVPP/CPCM. The Hartree-Fock plus London dispersion (HFLD) method was used for the study of non-covalent interactions (NCI). The calculations lead to the conclusion that a reduction in electrostatic interactions and an increase in exchange and dispersion interactions in CHCl3 in relation to DMSO are the driving forces behind the placement of sec-amine molecules into the R[4]A cavity and the formation of "in" type complexes.

Accurate Evaluation of Dispersion Energies at Coupled Cluster Level to Understand the Substituent Effects in Am(III) and Eu(III) Complexes

The effect of cyclic and aromatic substituents on the complexation behavior of phosphine oxide ligands with Am(III) and Eu(III) was investigated at density functional theory (DFT) and domain-based local pair natural orbital coupled-cluster (DLPNO-CC) levels. Combining DFT with accurate coupled cluster methods, we have evaluated the dispersion energy contributions to the complexation energies for trivalent Am and Eu complexes for the first time. Irrespective of the nature of substituents on the P atom, the electronic structure of the P═O group remains identical in all of the ligands. The study reveals the importance of dispersion interactions during complexation and is estimated to be more significant for Am(III) than for Eu(III) complexes.

Stacked binding of a PET ligand to Alzheimer's tau paired helical filaments

Accumulation of filamentous aggregates of tau protein in the brain is a pathological hallmark of Alzheimer's disease (AD) and many other neurodegenerative tauopathies. The filaments adopt disease-specific cross-β amyloid conformations that self-propagate and are implicated in neuronal loss. Development of molecular diagnostics and therapeutics is of critical importance. However, mechanisms of small molecule binding to the amyloid core is poorly understood. We used cryo-electron microscopy to determine a 2.7 Å structure of AD patient-derived tau paired-helical filaments bound to the PET ligand GTP-1. The compound is bound stoichiometrically at a single site along an exposed cleft of each protofilament in a stacked arrangement matching the fibril symmetry. Multiscale modeling reveals pi-pi aromatic interactions that pair favorably with the small molecule-protein contacts, supporting high specificity and affinity for the AD tau conformation. This binding mode offers critical insight into designing compounds to target different amyloid folds found across neurodegenerative diseases.

Exploring the Limits of Intramolecular London Dispersion Stabilization with Bulky Dispersion Energy Donors in Alkane Solution

We present an experimental study of a cyclooctatetraene-based molecular balance disubstituted with increasingly bulky tert-butyl (tBu), adamantyl (Ad), and diamantyl (Dia) substituents in the 1,4-/1,6-positions for which we determined the valence-bond shift equilibrium in n-hexane (hex), n-octane (oct), and n-dodecane (dod). Computations including implicit and explicit solvation support our temperature-dependent NMR equilibrium measurements indicating that the more sterically crowded 1,6-isomer is always favored, irrespective of solvent, and that the free energy is quite insensitive to substituent size.

Supramolecular interactions between ethylene-bridged oligoureas: nanorings and chains formed by cooperative positive allostery

Ethylene-bridged oligoureas are dynamic foldamers in which the polarity of a coherent chain of intramolecular hydrogen bonds may be controlled by intra- or intermolecular interactions with hydrogen-bond donors or acceptors. In this paper, we describe the way that supramolecular interactions between ethylene-bridged oligoureas bearing a 3,5-bis(trifluoromethyl)phenylurea (BTMP) terminus leads to higher-order structures both in the crystalline state and in solution. The oligoureas self-assemble by head-to-tail hydrogen bonding interactions to form either supramolecular 'nanorings' with cyclic hydrogen bond chain directionality, or supramolecular helical chains of hydrogen bonds. The self-assembly process features a cascade of cooperative positive allostery, in which each intermolecular hydrogen bond formation at the BTMP terminus switches the native hydrogen bond chain directionality of monomers, favouring further assembly. Monomers with a benzyl urea terminus self-assemble into nanorings, whereas monomers with a N-ethyl urea terminus form helical chains. In the crystal state, parallel helices have identical handedness and polarity, whereas antiparallel helices have opposite handedness. The overall dipole moment of crystals is zero due to the antiparallel arrangements of local dipoles in the crystal packing. Supramolecular interactions in solution were also examined by VT and DOSY NMR spectroscopy, up to the point of crystal formation. The size of higher aggregates in dichloromethane was estimated by their hydrodynamic radius. The relative orientation of the monomers within the aggregates, determined by 2D ROESY NMR, was the same as in the crystals, where syn-orientations lead to the formation of rings and anti-orientations result in chains. Overall, the switch of hydrogen bond polarity propagates intermolecularly in crystal and solution states, constituting an example of intermolecular communication within supramolecular polymers.

Janus faced fluorocyclohexanes for supramolecular assembly: synthesis and solid state structures of equatorial mono-, di- and tri alkylated cyclohexanes and with tri-axial C-F bonds to impart polarity

Concise and general synthesis protocols are reported to generate all-syn mono-, di- and tri-alkylated cyclohexanes where a single fluorine is located on the remaining carbons of the ring. The alkyl groups are positioned to lie equatorially and to have triaxial C-F bonds imparting polarity to these ring systems. Intermolecular electrostatic interactions in the solid-state structure of the trialkylated systems are explored and the resultant supramolecular order opens up prospects for design in soft materials.

Ab initio development of generalized Lennard-Jones (Mie) force fields for predictions of thermodynamic properties in advanced molecular-based SAFT equations of state

A methodology for obtaining molecular parameters of a modified statistical associating fluid theory for variable-range interactions of Mie form (SAFT-VR Mie) equation of state (EoS) from ab initio calculations is proposed for non-associative species that can be modeled as single spherical segments. The methodology provides a strategy to map interatomic or intermolecular potentials obtained from ab initio quantum-chemistry calculations to the corresponding Mie potentials that can be used within the SAFT-VR Mie EoS. The inclusion of corrections for quantum and many-body effects allows for an excellent, fully predictive description of the vapor-liquid envelope and other bulk thermodynamic properties of noble gases; this description is of similar or superior quality to that obtained using SAFT-VR Mie with parameters regressed in the traditional way using experimental thermodynamic-property data. The methodology is extended to an anisotropic species, methane, where similar levels of accuracy are obtained. The efficacy of using less-accurate quantum-chemistry methods in this methodology is explored, showing that these methods do not provide satisfactory results, although we note that the description is nevertheless substantially better than those obtained using the conductor-like screening model for describing real solvents (COSMO-RS), the only other fully predictive ab initio method currently available. Overall, the reliance on thermophysical data is completely dispensed with, providing the first extensible, wholly predictive SAFT-type EoSs.

Open-Shell Variant of the London Dispersion-Corrected Hartree-Fock Method (HFLD) for the Quantification and Analysis of Noncovalent Interaction Energies

The London dispersion (LD)-corrected Hartree–Fock (HF) method (HFLD) is an ab initio approach for the quantification and analysis of noncovalent interactions (NCIs) in large systems that is based on the domain-based local pair natural orbital coupled-cluster (DLPNO-CC) theory. In the original HFLD paper, we discussed the implementation, accuracy, and efficiency of its closed-shell variant. Herein, an extension of this method to open-shell molecular systems is presented. Its accuracy is tested on challenging benchmark sets for NCIs, using CCSD(T) energies at the estimated complete basis set limit as reference. The HFLD scheme was found to be as accurate as the best-performing dispersion-corrected exchange-correlation functionals, while being nonempirical and equally efficient. In addition, it can be combined with the well-established local energy decomposition (LED) for the analysis of NCIs, thus yielding additional physical insights.

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Unexpected triaxial preferences in some all-syn 1,3,5-trifluorocyclohexanes

Theory and solution NMR indicate that all-syn 1,3,5-trifluorocyclohexane adopts the expected tri-equatorial conformation, however in the solid state the more polar triaxial conformation is observed.

London Dispersion Helps Refine Steric A-Values: Dispersion Energy Donor Scales

We suggest a scale of dispersion energy donors (DEDs) that allows for direct comparisons with steric effects. This scale is based on the classic A-values and allows groups to reorient to minimize strain, thereby providing an advantage over raw group polarizabilities. The A-value can no longer be considered purely a steric factor. Even for groups that do not participate in charge transfer or electrostatic interactions, the A-value includes Pauli repulsion (steric hindrance) and attractive London dispersion (LD) interactions. Although the common assumption is that, at the distances found in monosubstituted cyclohexanes, steric demands are the key factors influencing conformer preferences, we show in this computational study that there is a non-negligible LD part. We use this system to build a DED scale and a complementary steric scale. These scales are quantitatively comparable, as they are based on the same system, and allow for comparison of the two competing interactions in experimentally relevant settings. In addition, we show that LD interactions can be used to explain puzzling data regarding relative group sizes.

Local energy decomposition analysis and molecular properties of encapsulated methane in fullerene (CH4@C60)

Methane has been successfully encapsulated within cages of C₆₀ fullerene, which is an appropriate model system to study confinement effects. Its chemistry and physics are also relevant for theoretical model descriptions. Here we provide insights into intermolecular interactions and predicted spectroscopic responses of the CH₄@C₆₀ complex and compared them with results from other methods and with data from the literature. Local energy decomposition analysis (LED) within the domain-based local pair natural orbital coupled cluster singles, doubles, and perturbative triples (DLPNO-CCSD(T)) framework was used, and an efficient protocol for studies of endohedral complexes of fullerenes is proposed. This approach allowed us to assess energies in relation to electronic and geometric preparation, electrostatics, exchange, and London dispersion for the CH₄@C₆₀ endohedral complex. The calculated stabilization energy of CH₄ inside the C₆₀ fullerene was -13.5 kcal mol^-1 and its magnitude was significantly larger than the latent heat of evaporation of CH₄. Evaluation of vibrational frequencies and polarizabilities of the CH₄@C₆₀ complex revealed that the infrared (IR) and Raman bands of the endohedral CH₄ were essentially "silent" due to the dielectric screening effect of C₆₀, which acted as a molecular Faraday cage. Absorption spectra in the UV-vis domain and ionization potentials of C₆₀ and CH₄@C₆₀ were predicted. They were almost identical. The calculated ¹H/¹³C NMR shifts and spin-spin coupling constants were in very good agreement with experimental data. In addition, reference DLPNO-CCSD(T) interaction energies for complexes with noble gases (Ng@C₆₀; Ng = He, Ne, Ar, Kr) were calculated. The values were compared with those derived from supramolecular MP2/SCS-MP2 calculations and estimates with London-type formulas by Pyykkö and coworkers [Phys. Chem. Chem. Phys., 2010, 12, 6187-6203], and with values derived from DFT-based symmetry-adapted perturbation theory (DFT-SAPT) by Hesselmann & Korona [Phys. Chem. Chem. Phys., 2011, 13, 732-743]. Selected points at the potential energy surface of the endohedral He₂@C₆₀ trimer were considered. In contrast to previous theoretical attempts with the DFT/MP2/SCS-MP2/DFT-SAPT methods, our calculations at the DLPNO-CCSD(T) level of theory predicted the He₂@C₆₀ trimer to be thermodynamically stable, which is in agreement with experimental observations.

Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

Herein, we provide new insights into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of a large portion of human B-DNA by using advanced quantum mechanical methods. Our results indicate that (i) the effect of non-neighboring bases on the inter-strand interaction is negligibly small, (ii) London dispersion effects are essential for the stability of the duplex structure, (iii) the largest contribution to the stability of the duplex structure is the Watson-Crick base pairing - consistent with previous computational investigations, (iv) the effect of stacking between adjacent bases is relatively small but still essential for the duplex structure stability and (v) there are no cooperativity effects between intra-strand stacking and inter-strand base pairing interactions. These results are consistent with atomic force microscope measurements and provide the first theoretical validation of nearest neighbor approaches for predicting thermodynamic data of arbitrary DNA sequences.

Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals

This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N⁶) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.

Fragment-Based Local Coupled Cluster Embedding Approach for the Quantification and Analysis of Noncovalent Interactions: Exploring the Many-Body Expansion of the Local Coupled Cluster Energy

Herein, we introduce a fragment-based local coupled cluster embedding approach for the accurate quantification and analysis of noncovalent interactions in molecular aggregates. Our scheme combines two different expansions of the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) energy: the many-body expansion (MBE) and the local energy decomposition (LED). The low-order terms in the MBE are initially computed in the presence of an environment that is treated at a low level of theory. Then, LED is used to decompose the energy of each term in the embedded MBE into additive fragment and fragment-pairwise contributions. This information is used to quantify the total energy of the system while providing at the same time in-depth insights into the nature and cooperativity of noncovalent interactions. Two different approaches are introduced and tested, in which the environment is treated at different levels of theory: the local coupled cluster in the Hartree-Fock (LCC-in-HF) method, in which the environment is treated at the HF level; and the electrostatically embedded local coupled cluster method (LCC-in-EE), in which the environment is replaced by point charges. Both schemes are designed to preserve as much as possible the accuracy of the parent local coupled cluster method for total energies, while being embarrassingly parallel and less memory intensive. These schemes appear to be particularly promising for the study of large and complex molecular aggregates at the coupled cluster level, such as condensed phase systems and protein-ligand interactions.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Unraveling individual host-guest interactions in molecular recognition from first principles quantum mechanics: Insights into the nature of nicotinic acetylcholine receptor agonist binding

Drug binding to a protein target is governed by a complex pattern of noncovalent interactions between the ligand and the residues in the protein's binding pocket. Here we introduce a generally applicable, parameter-free, computational method that allows for the identification, quantification, and analysis of the key ligand-residue interactions responsible for molecular recognition. Our strategy relies on Local Energy Decomposition analysis at the "gold-standard" coupled cluster DLPNO-CCSD(T) level. In the study case shown in this paper, nicotine and imidacloprid binding to the nicotinic acetylcholine receptor, our approach provides new insights into how individual amino acids in the active site determine sensitivity and selectivity of the ligands, extending and refining classical pharmacophore hypotheses. By inference, the method is applicable to any kind of host/guest interactions with potential applications in industrial biocatalysis and protein engineering.

Extrapolation to the Limit of a Complete Pair Natural Orbital Space in Local Coupled-Cluster Calculations

The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method allows one to perform single point energy calculations for systems with hundreds of atoms while retaining essentially the accuracy of its canonical counterpart, with errors that are typically smaller than 1 kcal/mol for relative energies. Crucial to the accuracy and efficiency of the method is a proper definition of the virtual space in which the coupled-cluster equations are solved, which is spanned by a highly compact set of pair natural orbitals (PNOs) that are specific for each electron pair. The dimension of the PNO space is controlled by the T_CutPNO threshold: only PNOs with an occupation number greater than T_CutPNO are included in the correlation space of a given electron pair, whilst the remaining PNOs are discarded. To keep the error of the method small, a conservative T_CutPNO value is used in standard DLPNO-CCSD(T) calculations. This often leads to unnecessarily large PNO spaces, which limits the efficiency of the method. Herein, we introduce a new computational strategy to approach the complete PNO space limit (for a given basis set) that consists in extrapolating the results obtained with different T_CutPNO values. The method is validated on the GMTKN55 set using canonical CCSD(T) data as the reference. Our results demonstrate that a simple two-point extrapolation scheme can be used to significantly increase the efficiency and accuracy of DLPNO-CCSD(T) calculations, thus extending the range of applicability of the technique.

Protobranching as repulsion-induced attraction: a prototype for geminal stabilization

Noncovalent interactions are traditionally defined within the context of their attractive components, such as electrostatics and dispersion. Sources of molecular strain are derived through the destabilization of Coulombic and exchange repulsion. Due to this binary designation, the underlying origin of geminal stability with respect to alkanes (referred to as protobranching) has been an active subject for debate between these competing perspectives. We recast this stabilization as a complementary (Gestalt) interaction between dispersion and exchange repulsion, each impacting the other. We use triplet hydrogen and argon dimer as foundational van der Waals adducts to develop a procedure for the visualization and quantification of both exchange repulsion, Δρ_SCF, and medium-range correlation, ΔΔρ, as perturbations in electron density. We use the framework of the DFT-D3 correction to reproduce the shape of the dispersion potential at medium range and successfully model the trend in stability for the eighteen isomers of octane with a diverse series of functionals: BLYP, B3LYP, BP86, PBE, and PBE0. Collectively, our findings show that protobranching is a manifestation of steric repulsion-reduction in vibrational enthalpy and medium-range electron correlation.