Extensive theoretical studies of the hydrogen‐bonded complexes (H2O)2, (H2O)2H+, (HF)2, (HF)2H+, F2H−, and (NH3)2

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis

Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.

Intermolecular interaction of diamine-diol binary system: A mini-review

The highly selective chemical reaction of carbon dioxide with organic amines is considered to a mature technology and a feasible initial path for carbon capture. In order to solve the disadvantages of high volatility, equipment corrosion and high energy consumption of traditional organic amines, amine alcohol "mixture based" solution has been developed and showed excellent carbon dioxide absorption capacity, which is due to the positive effect of intermolecular interaction in amine alcohol "mixture based" solution system on thermodynamic properties. However, the influencing factors of the intermolecular force in multicomponent solution system are complex, including the chemical, physical, structural effects. Therefore, it is necessary to comprehensively use a variety of characterization methods to systematically understand the form of intermolecular interaction in multicomponent solution system. This review systematically discusses the determination of intermolecular interactions in diamine-diol multicomponent solutions by three mainstream research methods, theoretical calculation method, spectral method, and thermodynamic method, aiming to provide a theoretical reference for the industrial production, the supplement to experimental data, and construction and understanding of theoretical models of multicomponent solution system.

Nonuniform Proton Transfer and Strong Hydrogen Bonding within Cation, Anion, and Neutral Clusters of Ammonia and Hydrogen Fluoride

Structural studies of ammonium halide nanoparticles can help to reveal fundamental information about the detailed nature of intermolecular forces. This study focuses on small ammonium fluoride clusters, which exhibit complex behavior in comparison to other ammonium halide clusters due to the weak acidity and strong hydrogen-bonding ability of HF. Calculations of optimized structures and binding energies are presented for cation, anion, and neutral clusters using MP2, CCSD(T), FNO-CCSD(T), ωB97M-V, and MN15 methods. The extent to which proton transfer occurs between two given cluster components was quantified using a dimensionless proton-transfer parameter (ξ_PT), leading to a classification of different types of hydrogen bonds within the clusters. Whereas the neutral clusters exhibit a complex transition from ordinary hydrogen bonding to a combination of shared-proton hydrogen bonds and complete proton transfers, the anion and cation systems exhibit a rapid transition toward complete proton transfer from HF to NH₃, with incomplete proton transfer observed only in the smallest anion and cation clusters. Ionic interaction energies of these clusters were also computed and found to exhibit trends which can be interpreted by the size-dependent behavior of ξ_PT. This work extends our understanding of the size-dependent trends in intermolecular forces which govern the formation of anhydrous ammonium halide clusters as well as the relationship between strong hydrogen bonding and proton transfer.

Assessing the Applicability of the Geometric Counterpoise Correction in B2PLYP/Double-ζ Calculations for Thermochemistry, Kinetics, and Noncovalent Interactions

We present a proof-of-concept study of the suitability of Kruse and Grimme’s geometric counterpoise correction (gCP) for basis set superposition errors (BSSEs) in double-hybrid density functional calculations with a double-ζ basis set. The gCP approach only requires geometrical information as an input and no orbital/density information is needed. Therefore, this correction is practically free of any additional cost. gCP is trained against the Boys and Bernardi counterpoise correction across a set of 528 noncovalently bound dimers. We investigate the suitability of the approach for the B2PLYP/def2-SVP level of theory, and reveal error compensation effects—missing London dispersion and the BSSE—associated with B2PLYP/def2-SVP calculations, and present B2PLYP-gCP-D3(BJ)/def2-SVP with the reparametrised DFT-D3(BJ) and gCP corrections as a more balanced alternative. Benchmarking results on the S66x8 benchmark set for noncovalent interactions and the GMTKN55 database for main-group thermochemistry, kinetics, and noncovalent interactions show a statistical improvement of the B2PLYP-gCP-D3(BJ) scheme over plain B2PLYP and B2PLYP-D3(BJ). B2PLYP-D3(BJ) shows significant overestimation of interaction energies, barrier heights with larger deviations from the reference values, and wrong relative stabilities in conformers, all of which can be associated with BSSE. We find that the gCP-corrected method represents a significant improvement over B2PLYP-D3(BJ), particularly for intramolecular noncovalent interactions. These findings encourage future developments of efficient double-hybrid DFT strategies that can be applied when double-hybrid calculations with large basis sets are not feasible due to system size.

On the Use of Popular Basis Sets: Impact of the Intramolecular Basis Set Superposition Error

The magnitude of intramolecular basis set superposition error (BSSE) is revealed via computing systematic trends in molecular properties. This type of error is largely neglected in the study of the chemical properties of small molecules and it has historically been analyzed just in the study of large molecules and processes dominated by non-covalent interactions (typically dimerization or molecular complexation and recognition events). In this work we try to provide proof of the broader prevalence of this error, which permeates all types of electronic structure calculations, particularly when employing insufficiently large basis sets.

Performance of polarization-consistent vs. correlation-consistent basis sets for CCSD(T) prediction of water dimer interaction energy

Detailed study of Jensen’s polarization-consistent vs. Dunning’s correlation-consistent basis set families performance on the extrapolation of raw and counterpoise-corrected interaction energies of water dimer using coupled cluster with single, double, and perturbative correction for connected triple excitations (CCSD(T)) in the complete basis set (CBS) limit are reported. Both 3-parameter exponential and 2-parameter inverse-power fits vs. the cardinal number of basis set, as well as the number of basis functions were analyzed and compared with one of the most extensive CCSD(T) results reported recently. The obtained results for both Jensen- and Dunning-type basis sets underestimate raw interaction energy by less than 0.136 kcal/mol with respect to the reference value of − 4.98065 kcal/mol. The use of counterpoise correction further improves (closer to the reference value) interaction energy. Asymptotic convergence of 3-parameter fitted interaction energy with respect to both cardinal number of basis set and the number of basis functions are closer to the reference value at the CBS limit than other fitting approaches considered here. Separate fits of Hartree-Fock and correlation interaction energy with 3-parameter formula additionally improved the results, and the smallest CBS deviation from the reference value is about 0.001 kcal/mol (underestimated) for CCSD(T)/aug-cc-pVXZ calculations. However, Jensen’s basis set underestimates such value to 0.012 kcal/mol. No improvement was observed for using the number of basis functions instead of cardinal number for fitting.

Intramolecular hydrogen bonding in malonaldehyde and its radical analogues

High level Brueckner doubles with triples correction method-based ab initio calculations have been used to investigate the nature of intramolecular hydrogen bonding and intramolecular hydrogen atom transfer in cis-malonaldehyde (MA) and its radical analogues. The radicals considered here are the ones that correspond to the homolytic cleavage of C-H bonds in cis-MA. The results suggest that cis-MA and its radical analogues, cis-MA_RS, and cis-MA_RA, both exist in planar geometry. The calculated intramolecular O-H⋯O=C bond in cis-MA is shorter than that in the radical analogues. The intramolecular hydrogen bond in cis-MA is stronger than in its radicals by at least 3.0 kcal/mol. The stability of a cis-malonaldehyde radical correlates with the extent of electron spin delocalization; cis-MA_RA, in which the radical spin is more delocalized, is the most stable MA radical, whereas cis-MA_RS, in which the radical spin is strongly localized, is the least stable radical. The natural bond orbital analysis indicates that the intramolecular hydrogen bonding (O⋯H⋯O) in cis-malonaldehyde radicals is stabilized by the interaction between the lone pair orbitals of donor oxygen and the σ^* orbital of acceptor O-H bond (n → σ^*_OH). The calculated barriers indicate that the intramolecular proton transfer in cis-MA involves 2.2 kcal/mol lower barrier than that in cis-MA_RS.

Push it to the limit: Characterizing the convergence of common sequences of basis sets for intermolecular interactions as described by density functional theory

With the aim of systematically characterizing the convergence of common families of basis sets such that general recommendations for basis sets can be made, we have tested a wide variety of basis sets against complete-basis binding energies across the S22 set of intermolecular interactions-noncovalent interactions of small and medium-sized molecules consisting of first- and second-row atoms-with three distinct density functional approximations: SPW92, a form of local-density approximation; B3LYP, a global hybrid generalized gradient approximation; and B97M-V, a meta-generalized gradient approximation with nonlocal correlation. We have found that it is remarkably difficult to reach the basis set limit; for the methods and systems examined, the most complete basis is Jensen's pc-4. The Dunning correlation-consistent sequence of basis sets converges slowly relative to the Jensen sequence. The Karlsruhe basis sets are quite cost effective, particularly when a correction for basis set superposition error is applied: counterpoise-corrected def2-SVPD binding energies are better than corresponding energies computed in comparably sized Dunning and Jensen bases, and on par with uncorrected results in basis sets 3-4 times larger. These trends are exhibited regardless of the level of density functional approximation employed. A sense of the magnitude of the intrinsic incompleteness error of each basis set not only provides a foundation for guiding basis set choice in future studies but also facilitates quantitative comparison of existing studies on similar types of systems.

Testing the CP-correction procedure with different DFT methods on H-bonding complexes of κ-carrabiose with water molecules

Interaction of water molecules with κ-carrabiose disaccharide, within three H-bonding complexes, was investigated. Particular interest was focused on the way with which the BSSE correction has to be performed. Two strategies were used, either performing BSSE correction during or after optimization. For this aim, several DFT-functionals (hybrid GGA and hybrid meta-GGA) and 6-31 + G* basis set were considered. The results demonstrated the uselessness of including of BSSE-CP correction during optimization for all complexes. From a structural point of view, a proper H-bonding description was obtained using the PBE0 functional for all complexes. The basis set effect on the BSSE using B3LYP functional was also investigated. The reliability of B3LYP/6-31 + G** and B3LYP/6-31++G** models for the complexes involving one or two water molecules was reported while the use of B3LYP/6-311 + G** or B3LYP/6-311++G** levels was shown to be more appropriate for larger complexes equivalent to that involving three water molecules. CP-corrected interaction energies were demonstrated to be closer to CBS-4 M interaction energies than the uncorrected ones. Graphical abstract Functional and basis set effects on BSSE.

Molecularly Tuning the Radicaloid N-H···O═C Hydrogen Bond

Substituent effects on the open shell N-H···O═C hydrogen-bond has never been reported. This study examines how 12 functional groups composed of electron donating groups (EDG), halogen atoms and electron withdrawing groups (EWG) affect the N-H···O═C hydrogen-bond properties in a six-membered cyclic model system of O═C(Y)-CH═C(X)N-H. It is found that group effects on this open shell H-bonding system are significant and have predictive trends when X = H and Y is varied. When Y is an EDG, the N-H···O═C hydrogen-bond is strengthened; and when Y is an EWG, the bond is weakened; whereas the variation in electronic properties of X group do not exhibit a significant impact upon the hydrogen bond strength. The structural impact of the stronger N-H···O═C hydrogen-bond are (1) shorter H and O distance, r(H···O) and (2) a longer N-H bond length, r(NH). The stronger N-H···O═C hydrogen-bond also acts to pull the H and O in toward one another which has an effect on the bond angles. Our findings show that there is a linear relationship between hydrogen-bond angle and N-H···O═C hydrogen-bond energy in this unusual H-bonding system. In addition, there is a linear correlation of the r(H···O) and the hydrogen bond energy. A short r(H···O) distance corresponds to a large hydrogen bond energy when Y is varied. The observed trends and findings have been validated using three different methods (UB3LYP, M06-2X, and UMP2) with two different basis sets.

Quantum chemical MP2 results on some hydrates of cytosine: binding sites, energies and the first hydration shell

A detailed quantum chemical investigation was undertaken to obtain the structure and energetics of cytosine hydrates Cyt·nH2O, with n = 1 to 7. The MP2(fc)/aug-cc-pVDZ level was used as the standard, with some DFT (B3LYP) and coupled cluster calculations, as well as calculations with the aug-cc-pVTZ basis set added for comparison. In a systematic search for microhydrated forms of cytosine, we have found that several structures have not yet been reported in the literature. The energies of different isomers, as well as binding energies are compared. When predicting the stability of a complex, we suggest using a scheme where the water molecules are extracted from a finite model of bulk water. Finally, based on energetic data, we suggest a rational definition of the first hydration shell; with this definition, it contains just six water molecules.

An overlooked yet ubiquitous fluoride congenitor: binding bifluoride in triazolophanes using computer-aided design

Despite its ubiquity during the binding and sensing of fluoride, the role of bifluoride (HF2(-)) and its binding properties are almost always overlooked. Here, we give one of the first examinations of bifluoride recognition in which we use computer-aided design to modify the cavity shape of triazolophanes to better match with HF2(-). Computational investigation indicates that HF2(-) and Cl(-) should have similar binding affinities to the parent triazolophane in the gas phase. Evaluation of the binding geometries revealed a preference for binding of the linear HF2(-) along the north-south axis with a smaller Boltzmann weighted population aligned east-west and all states being accessed rapidly through in-plane precessional rotations of the anion. While the (1)H NMR spectroscopy studies are consistent with the calculated structural aspects, binding affinities in solution were determined to be significantly smaller for the bifluoride than the chloride. Computed geometries suggested that a 20° tilting of the bifluoride (stemming from the cavity size) could account for the 25-fold difference between the two binding affinities, HF2(-) < Cl(-). Structural variations to the triazolophane's geometry and electronic modifications to the network of hydrogen bond donors were subsequently screened in a stepwise manner using density functional theory calculations to yield a final design that eliminates the tilting. Correspondingly, the bifluoride's binding affinity (K ∼ 10(6) M(-1)) increased and was also found to remain equal to chloride in the gas and solution phases. The new oblate cavity appeared to hold the HF2(-) in a single east-west arrangement. Our findings demonstrate the promising ability of computer-aided design to fine-tune the structural and electronic match in anion receptors as a means to control the arrangement and binding strength of a desired guest.

Strength of hydrogen bonds of water depends on local environment

In-depth knowledge of water-water potential is important for devising and evaluating simple water models if they are to accurately describe water properties and reflect various solvation phenomena. Water-water potential depends upon inter-molecular distance, relative orientation of water molecules, and also local environment. When placed at a favorable distance in a favorable orientation, water molecules exhibit a particularly strong attractive interaction called hydrogen bond. Although hydrogen bond is very important for its effects on the elements of life, industrial applications, and bulk water properties, there is no scientific consensus on its true nature and origin. Using quantum-mechanical methods, hydrogen bond strength was calculated in different local environments. A simple empirical linear relationship was discovered between maximum hydrogen bond strength and the number of water molecules in the local environment. The local environment effect was shown to be considerable even on the second coordination shell. Additionally, a negative linear correlation was found between maximum hydrogen bond strength and the distance, at which it was observed. These results provide novel insights into the nature of hydrogen bonding.

Complexation of formaldoxime with water. Infrared matrix isolation and theoretical studies

The 1:1, 1:2 and 2:1 formaldoxime-water complexes isolated in the argon matrices have been studied by help of FTIR spectroscopy and MP2/6-311++G(2d,2p) method. The calculations predicted the stability of the three CH(2)NOH···H(2)O isomeric complexes, three CH(2)NOH···(H(2)O)(2) ones and one (CH(2)NOH)(2)···H(2)O complex. The analysis of the experimental spectra and their comparison with theoretical ones indicated that both the 1:1 and 1:2 complexes trapped in solid argon have the most stable cyclic structures stabilized by the O-H···O and O-H···N bonds between the formaldoxime and water molecules. In the 1:2 complex formaldoxime interacts with the water dimer, one H(2)O molecule acts as a proton acceptor for the OH group of formaldoxime whereas the second H(2)O molecule acts as a proton donor toward the nitrogen atom of the formaldoxime molecule. In the (CH(2)NOH)(2)···H(2)O complex the OH group of the water molecule acts as a proton donor toward one of the oxygen atoms of the formaldoxime cyclic dimer.

INTERACTION OF EMODIN WITH DNA BASES: A DENSITY FUNCTIONAL THEORY

In this study, we present work on the physicochemical interaction between the anticancer drug molecule Emodin (ED) and DNA. Comprehending the physicochemical properties of this drug besides the mechanism by which it interacts with DNA should eventually permit the rational design of novel anticancer or antiviral drugs. The final purpose is the clarification of this novel class of drugs as potential pharmaceutical agents. The properties of the isolated intercalator ED and its stacking interactions with adenine⋯thymine (AT) and guanine⋯cytosine (GC) (nucleic acid base pairs) in face-to-face and face-to-back models were studied by means of the density functional tightbinding (DFTB) method. This method was an approximate version of the density functional theory (DFT) method and it includes London dispersion energy. The molecular modeling of the complex formed between ED and DNA indicated that this complex was capable of contributing to the formation of a constant intercalation site. The results exhibit that ED changes affect DNA structure with reference to bond lengths, bond angles, torsion angles, and charges.

Cation-anion hydrogen bonds: a new class of hydrogen bonds that extends their strength beyond the covalent limit. A theoretical characterization

The existence of O-H···O hydrogen bonds having a strength within the -80 to -210 kcal/mol range, that is, in the range of strength of covalent bonds and well beyond the so-called covalent limit (-50 kcal/mol), is reported on complexes where the O-H proton donor and O acceptor groups are located in ions of opposite sign. A complete analysis of short distance O-H···O hydrogen bonds between charged fragments was performed for cases where the OH and O groups are both located on charged molecules. It shows that these interactions (a) are nonsymmetrical for the O-H and H···O distances, (b) have a noncovalent H···O bond critical point, and (c) have a strong and energetically stable electrostatic component when the OH and O groups are located in oppositely charged molecules. These cation-anion O-H···O interactions are energetically stable, satisfy the usual topology for hydrogen bonds, HBs, and also have the same directionality found in other HBs. Therefore, they should be considered as a new class of HBs, the cation-anion hydrogen bonds.

Counterpoise correction is not useful for short and Van der Waals distances but may be useful at long range

This article investigates the errors in supermolecule calculations for the helium dimer. In a full CI calculation, there are two errors. One is the basis set superposition error (BSSE), the other is the basis set convergence error (BSCE). Both of the errors arise from the incompleteness of the basis set. These two errors make opposite contributions to the interaction energies. The BSCE is by far the largest error in the short range and larger than (but much closer to) BSSE around the Van der Waals minimum. Only at the long range, the BSSE becomes the larger error. The BSCE and BSSE largely cancel each other over the Van der Waals well. Accordingly, it may be recommended to not include the BSSE for the calculation of the potential energy curve from short distance till well beyond the Van der Waals minimum, but it may be recommended to include the BSSE correction if an accurate tail behavior is required. Only if the calculation has used a very large basis set, one can refrain from including the counterpoise correction in the full potential range. These results are based on full CI calculations with the aug-cc-pVXZ (X = D, T, Q, 5) basis sets.

A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

We evaluate the performance of ten functionals (B3LYP, M05, M05-2X, M06, M06-2X, B2PLYP, B2PLYPD, X3LYP, B97D, and MPWB1K) in combination with 16 basis sets ranging in complexity from 6-31G(d) to aug-cc-pV5Z for the calculation of the H-bonded water dimer with the goal of defining which combinations of functionals and basis sets provide a combination of economy and accuracy for H-bonded systems. We have compared the results to the best non-density functional theory (non-DFT) molecular orbital (MO) calculations and to experimental results. Several of the smaller basis sets lead to qualitatively incorrect geometries when optimized on a normal potential energy surface (PES). This problem disappears when the optimization is performed on a counterpoise (CP) corrected PES. The calculated interaction energies (ΔEs) with the largest basis sets vary from -4.42 (B97D) to -5.19 (B2PLYPD) kcal/mol for the different functionals. Small basis sets generally predict stronger interactions than the large ones. We found that, because of error compensation, the smaller basis sets gave the best results (in comparison to experimental and high-level non-DFT MO calculations) when combined with a functional that predicts a weak interaction with the largest basis set. As many applications are complex systems and require economical calculations, we suggest the following functional/basis set combinations in order of increasing complexity and cost: (1) D95(d,p) with B3LYP, B97D, M06, or MPWB1k; (2) 6-311G(d,p) with B3LYP; (3) D95++(d,p) with B3LYP, B97D, or MPWB1K; (4) 6-311++G(d,p) with B3LYP or B97D; and (5) aug-cc-pVDZ with M05-2X, M06-2X, or X3LYP.

Assessment of the Performance of DFT and DFT-D Methods for Describing Distance Dependence of Hydrogen-Bonded Interactions

Noncovalent interactions such as hydrogen bonds, van der Waals forces, and π-π interactions play important roles influencing the structure, stability, and dynamic properties of biomolecules including DNA and RNA base pairs. In an effort to better understand the fundamental physics of hydrogen bonding (H-bonding), we investigate the distance dependence of interaction energies in the prototype bimolecular complexes of formic acid, formamide, and formamidine. Potential energy curves along the H-bonding dissociation coordinate are examined both by establishing reference CCSD(T) interaction energies extrapolated to the complete basis set limit and by assessing the performance of the density functional methods B3LYP, PBE, PBE0, B970, PB86, M05-2X, and M06-2X and empirical dispersion corrected methods B3LYP-D3, PBE-D3, PBE0-D3, B970-D2, BP86-D3, and ωB97X-D, with basis sets 6-311++G(3df,3pd), aug-cc-pVDZ, and aug-cc-pVTZ. Although H-bonding interactions are dominated by electrostatics, it is necessary to properly account for dispersion interactions to obtain accurate energetics. In order to quantitatively probe the nature of hydrogen bonding interactions as a function of distance, we decompose the interaction energy curves into physically meaningful components with symmetry-adapted perturbation theory (SAPT). The SAPT results confirm that the contribution of dispersion and induction are significant at and near equilibrium, although electrostatics dominate. Among the DFT/DFT-D techniques, the best overall results are obtained utilizing counterpoise-corrected ωB97X-D with the aug-cc-pVDZ basis set.

Complexes in the photocatalytic reaction of CO(2) and H(2)O: theoretical studies

Complexes (H₂O/CO₂, e–(H₂O/CO₂) and h⁺–(H₂O/CO₂)) in the reaction system of CO₂ photoreduction with H₂O were researched by B3LYP and MP2 methods along with natural bond orbital (NBO) analysis. Geometries of these complexes were optimized and frequencies analysis performed. H₂O/CO₂ captured photo-induced electron and hole produced e–(H₂O/CO₂) and h⁺–(H₂O/CO₂), respectively. The results revealed that CO₂ and H₂O molecules could be activated by the photo-induced electrons and holes, and each of these complexes possessed two isomers. Due to the effect of photo-induced electrons, the bond length of C=O and H-O were lengthened, while H-O bonds were shortened, influenced by holes. The infrared (IR) adsorption frequencies of these complexes were different from that of CO₂ and H₂O, which might be attributed to the synergistic effect and which could not be captured experimentally.

Does ammonia hydrogen bond?

Spectroscopic characterizations of the stereochemistry of complexes of ammonia (NH(3)) have strongly confirmed some long-held ideas about the weak interactions of NH(3) while casting doubt on others. As expected, NH(3) is observed to be a nearly universal proton acceptor, accepting hydrogen bonds from even some of the weakest proton donors. Surprisingly, no evidence has been found to support the view that NH(3) acts as a proton donor through hydrogen bonding. A critical evaluation of the work that has been done to gather such evidence, as well as of earlier work involving condensed-phase observations, suggests that NH(3) might well be best described as a powerful hydrogen-bond acceptor with little propensity to donate hydrogen bonds.