Hydrogen bonding in the mixed HF/HCl dimer: Is it better to give or receive?

The ClH⋯FH and FH⋯ClH configurations of the mixed HF/HCl dimer (where the donor⋯acceptor notation indicates the directionality of the hydrogen bond) as well as the transition state connecting the two configurations have been optimized using MP2 and CCSD(T) with correlation consistent basis sets as large as aug‐cc‐pV(5 + d)Z. Harmonic vibrational frequencies confirmed that both configurations correspond to minima and that the transition state has exactly one imaginary frequency. In addition, anharmonic vibrational frequencies computed with second‐order vibrational perturbation theory (VPT2) are within 6 cm⁻¹ of the available experimental values and deviate by no more than 4 cm⁻¹ for the complexation induced HF frequency shifts. The CCSD(T) electronic energies obtained with the largest basis set indicate that the barrier height is 0.40 kcal mol⁻¹ and the FH⋯ClH configuration lies 0.19 kcal mol⁻¹ below the ClH⋯FH configuration. While only modestly attenuating the barrier height, the inclusion of either the harmonic or anharmonic zero‐point vibrational energy effectively makes both minima isoenergetic, with the ClH⋯FH configuration being lower by only 0.03 kcal mol⁻¹. © 2018 Wiley Periodicals, Inc.     

Theoretical study on structures and vibrational spectra of M+(H2O)Ar (M = Cu,Ag, Au)

A theoretical study on the structures and vibrational spectra of M⁺(H₂O)Ar_0‐1 (M = Cu, Ag, Au) complexes was performed using ab initio method. Geometrical structures, binding energies (BEs), OH stretching vibrational frequencies, and infrared (IR) absorption intensities are investigated in detail for various isomers with Ar atom bound to different binding sites of M⁺(H₂O). CCSD(T) calculations predict that BEs are 14.5, 7.5, and 14.4 kcal/mol for Ar atom bound to the noble metal ion in M⁺(H₂O)Ar (M = Cu, Ag, Au) complexes, respectively, and the corresponding values have been computed to be 1.5, 1.3, and 2.1 kcal/mol when Ar atom attaches to a H atom of water molecule. The former structure is predicted to be more stable than the latter structure. Moreover, when compared with the M⁺(H₂O) species, tagging Ar atom to metal cation yields a minor perturbation on the IR spectra, whereas binding Ar atom to an OH site leads to a large redshift in OH stretching vibrations. The relationships between isomers and vibrational spectra are discussed. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

The HOX⋯SO2 (X=F,Cl, Br,I) Binary Complexes: Implications for Atmospheric Chemistry

Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO₂ binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug-cc-pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug-cc-pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug-cc-pwCVDZ-PP pseudopotential. 27 HOX⋯SO₂ complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO₂ complexes in this study ranged from 1.35 to 3.81 kcal mol⁻¹. The single-interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol⁻¹, while the single-interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol⁻¹, indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO₂ which may guide further research of related systems.     

The Beryllium Pentamer: Trailing an Uneven Sequence of Dissociation Energies

Recent high‐resolution spectroscopic studies by Merritt, Bondybey, and Heaven (Science 2009 , 324, 1548) have heightened the anticipation that small beryllium clusters will soon be observed in the laboratory. Beryllium clusters are important discrete models for the theoretical study of metals. The trigonal bipyramidal Be₅ molecule is studied using high‐level coupled cluster methods. We obtain the optimized geometry, atomization and dissociation energies, and vibrational frequencies. The c～CCSDT(Q) method is employed to compute the atomization and dissociation energies. In this approach, complete basis set (CBS) extrapolations at the CCSD(T) level of theory are combined with an additive correction for the effect of iterative triple and perturbative quadruple excitations. Harmonic vibrational frequencies are obtained using analytic gradients computed at the CCSD(T) level of theory. We report an atomization energy of 129.6 kcal mol^?1 at the trigonal bipyramid global minimum geometry. The Be₅→Be₄+Be dissociation energy is predicted to be 39.5 kcal mol^?1. The analogous dissociation energies for the smaller beryllium clusters are 64.0 kcal mol^?1 (Be₄→Be₃+Be), 24.2 kcal mol^?1 (Be₃→Be₂+Be), and 2.7 kcal mol^?1 (Be₂→Be+Be). The trigonal bipyramidal Be₅ structure has an equatorial–equatorial bond length of 2.000 Å and an axial–equatorial distance of 2.060 Å. Harmonic frequencies of 730, 611, 456, 583, 488, and 338 cm^?1 are obtained at the CCSD(T)/cc‐pCVQZ level of theory. Quadruple excitations are found to make noticeable contributions to the energetics of the pentamer, which exhibits a significant level of static correlation.     

An ab initio CI investigation of structure and bonding in LiC2H2 complexes

The structures and energies of various LiC₂H₂ complexes have been investigated by means of ab initio molecular orbital calculations. Analytic SCF gradients were employed with a double-ζ basis set to locate and characterize stationary points on the energy surface. Single-point CI calculations using a double-ζ + diffuse and polarization basis set have been carried out at the DZ + P SCF stationary points. With the highest-level theory, the Li—vinylidene complex and the cis bridged adduct are found to be the most favorable arrangements, the former complex being slightly more stable by about 2 kcal mol⁻¹. These molecules are bound respectively by about 5 and 3 kcal mole⁻¹ relative to infinitely separated lithium plus acetylene. Harmonic vibrational frequencies are also reported and confirm the existence of the cis LiC₂H₂ species recently observed in a solid argon matrix.     

Formic acid catalyzed gas-phase reaction of H2O with SO3 and the reverse reaction: a theoretical study

The formic acid catalyzed gas‐phase reaction between H₂O and SO₃ and its reverse reaction are respectively investigated by means of quantum chemical calculations at the CCSD(T)//B3LYP/cc‐pv(T+d)z and CCSD(T)//MP2/aug‐cc‐pv(T+d)z levels of theory. Remarkably, the activation energy relative to the reactants for the reaction of H₂O with SO₃ is lowered through formic acid catalysis from 15.97 kcal mol^?1 to ?15.12 and ?14.83 kcal mol^?1 for the formed H₂O ??? SO₃ complex plus HCOOH and the formed H₂O ??? HCOOH complex plus SO₃, respectively, at the CCSD(T)//MP2/aug‐cc‐pv(T+d)z level. For the reverse reaction, the energy barrier for decomposition of sulfuric acid is reduced to ?3.07 kcal mol^?1 from 35.82 kcal mol^?1 with the aid of formic acid. The results show that formic acid plays a strong catalytic role in facilitating the formation and decomposition of sulfuric acid. The rate constant of the SO₃+H₂O reaction with formic acid is 10⁵ times greater than that of the corresponding reaction with water dimer. The calculated rate constant for the HCOOH+H₂SO₄ reaction is about 10^?13 cm³ molecule^?1 s^?1 in the temperature range 200–280 K. The results of the present investigation show that formic acid plays a crucial role in the cycle between SO₃ and H₂SO₄ in atmospheric chemistry.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

Exploring new species on the [H,S, Se,Cl] potential energy surface

This investigation is concerned with the characterization of seleno‐sulfide‐halogen model systems, the isomerization processes, and the dissociation into diatomic fragment channels on the [H, S, Se, Cl] potential energy surface. Structural, energetic, and vibrational data were obtained at the CCSD(T) and MP2 levels of theory with the series of correlation consistent basis sets and extrapolated to the complete basis set (CBS) limit. For the frequencies, additional computations were performed to include the contribution of anharmonic effects, and for the determination of the heats of formation, important corrections incorporating core‐valence correlation effects and relativistic effects (scalar and spin‐orbit) were taken into account. CCSD(T)/CBS relative stability (kcal mol^?1) follows the order: HSSeCl (0.0), HSeSCl (8.80), SSeHCl (23.52), and SeSHCl (25.87). The cis‐rotational barrier for the two lowest isomers is practically identical (10.14 and 10.09 kcal mol^?1), whereas for the trans barrier, we obtained 9.25 (HSSeCl) and 8.45 (HSeSCl) kcal mol^?1. Dissociation of HSSeCl (HSeSCl) into HS (HSe) + SeCl (SCl) requires 59.70 (56.30) kcal mol^?1. For the most stable isomer, we predict a value of the heat of formation at 298.15 K of 2.53 kcal mol^?1. One of the outcomes of this research is that the MP2 results are consistent with those of CCSD(T). The MP2 method turns out to be a reliable alternative for a first exploration of larger catenated species, although it accounts for a lesser fraction of correlation effects. © 2012 Wiley Periodicals, Inc.     

Silaacetylene: A possible target for experimental studies

The equilibrium geometries and transition states for interconversion of the CSiH₂ isomers in the singlet electronic ground state are optimized at the MP2 and CCSD(T) levels of theory using a TZ2P basis set. The heats of formation, vibrational frequencies, infrared intensities, and rotational constants are also predicted. There are three energy minima on the CSiH₂ potential energy surface. Energy calculations at CCSD(T)/TZ2P(fd) + ZPE predict that the global energy minimum is silavinylidene (1), which is 34.1 kcal mol⁻¹ lower in energy than trans-bent silaacetylene (2) and 84.1 kcal mol⁻¹ more stable than the vinylidene isomer (3). The barrier for rearrangement 2→1 is calculated at the same level of theory to be 5.1 kcal mol⁻¹, while for the rearrangement 3→2 a barrier of 2.7 kcal mol⁻¹ is predicted. The natural bond orbital (NBO) population scheme indicates a clear polarization of the C(SINGLE BOND)Si bonds toward the carbon end. A significant ionic contribution to the C(SINGLE BOND)Si bonds of 1 and 2 is suggested by the NBO analysis. The C(SINGLE BOND)Si bond length of trans-bent silaacetylene (2) is longer than previously calculated [1.665 Å at CCSD(T)/TZ2P)]. The calculated carbon-silicon bond length of 2 is in the middle between the C(SINGLE BOND)Si double bond length of 1 (1.721 Å) and the C(SINGLE BOND)Si triple bond of the linear form HCSiH (4), which is 1.604 Å. Structure 4 is a higher-order saddle point on the potential energy surface. © 1996 by John Wiley & Sons, Inc.     

Competing Excited-State Hydrogen and Proton-Transfer Processes in 6-Azaindole-S3,4 and 2,6-Diazaindole-S3,4 Clusters (S=H2O,NH3)

Excited state hydrogen (ESHT) and proton (ESPT) transfer reaction pathways in the three and four solvent clusters of 6-azaindole (6AI-S_3,4) and 2,6-diazaindole (26DAI-S_3,4)(S=H₂O, NH₃) were computationally investigated to understand the fate of photo-excited biomolecules. The ESHT energy barriers in (H₂O)₃ complexes (39.6–41.3 kJmol⁻¹) were decreased in (H₂O)₄ complexes (23.1–20.2 kJmol⁻¹). Lengthening the solvent chain lowered the barrier because of the relaxed transition states geometries with reduced angular strains. Replacing the water molecule with ammonia drastically decreased the energy barriers to 21.4–21.3 kJmol⁻¹ in (NH₃)₃ complexes and 8.1–9.5 kJ mol⁻¹ in (NH₃)₄ complexes. The transition states were identified as H_a atom attached to the first solvent molecule. The formation of stronger hydrogen bonds in (NH₃)_3,4 complexes resulted in facile ESHT reaction than that in the (H₂O)_3,4 complexes. The ESPT energy barriers in 6AI-S_3,4 and 26DAI-S_3,4 were found to range between 40–73 kJmol⁻¹. The above values were significantly higher than that of the ESHT processes and hence are considered as a minor channel in the process. The effect of N(2) insertion was explored for the very first time in the isolated solvent clusters using local vibrational mode analysis. In DAI-S₄, the higher K_a(H_a⋯S_a) values depicted the increased photoacidity of the N(1)-H_a group which may facilitate the hydrogen transfer reaction. However, the increased N(6)⋯H_b bond length elevated the reaction barriers. Therefore, in the ESHT reaction channel, the co-existence of two competing factors led to a marginal/no change in the overall energy barriers due to the N(2) insertion. In the ESPT reaction pathway, the energy barriers showed notable increase upon N(2) insertion because of the increased N(6)⋯H_b bond length.     

DFT study on the hydrogen bonds of phenol–cyclohexanone and phenol–H2O2 in the Baeyer–Villiger oxidation

Hydrogen bonds of phenol–cyclohexanone and phenol–H₂O₂ in the studied Baeyer–Villiger (B–V) oxidation have been investigated by HF, B3LYP, and MP2 methods with various basis sets. The accurate single‐point energies were performed using CCSD(T)/6‐31+G(d,p) and CCSD(T)/aug‐cc‐pVDZ on the optimized geometries of MP2/6‐31+G(d,p). It has been confirmed that B3LYP/6‐31+G(d,p) could be used to study such hydrogen bonds. Energetic analysis of complexes was carried out using the Xantheas method with BSSE corrected by CP method. Orbital energy order (ε) illuminated that phenol with good hydrogen donor‐acceptor property can interact with cyclohexanone or H₂O₂ to form hydrogen bound complexes, and the binding energies (BE) range from ?4.38 to ?14.06 kcal mol^?1. NBO analysis indicated that the redistribution of atomic charges in the complexes facilitated nucleophilic attack of H₂O₂ on cyclohexanone. The calculated results match remarkably well with the experimental phenomena. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

The proton‐coupled proton transfer mechanism,H2O catalysis,and hydrogen tunneling effects in the reaction of HNCH2 with HCOOH in the interstellar medium

The reaction HNCH₂ + HCOOH → H₂NCH₂COOH is supposed to be an important reaction related to the possible origin of amino acids on the early Earth. We find that it has an energy barrier of 87.37 kcal mol⁻¹ obtained with MP2/6‐311+G** in the gas phase, but it is likely enhanced to occur in the interstellar medium (ISM) through a proton‐coupled proton transfer reaction, initiated by HNCH₂ coupled with H₂⁺, H₃⁺, or H₃⁺O. H₂⁺, H₃⁺, and H₃⁺O serve as a donor of energy in the coupled reactions. H⁺, which is a key species to the coupled reactions, further, plays a catalytic role in reducing a barrier up to 14.14 kcal mol⁻¹. In the coupled reaction with H₃⁺O, H₂O, which can seize, transport, and deliver a proton from HCOOH to H₂NCH₂⁺, reduces a barrier up to 14.96 kcal mol⁻¹. A significant hydrogen‐tunneling pathway is predicted by the temperature dependences of k_H^CVT/SCT, calculated using the small curvature tunneling (SCT) approximation and canonical variational transition state theory (CVT). Hydrogen tunneling is another important mechanism to make the reaction happen in the ISM. The achieved results can be applied to discuss the origin of amino acids from the materials of the Earth itself. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

London Dispersion Effects in a Distannene/Tristannane Equilibrium: Energies of their Interconversion and the Suppression of the Monomeric Stannylene Intermediate

Reaction of {LiC₆H₂−2,4,6-Cyp₃⋅Et₂O}₂ (Cyp=cyclopentyl) ( 1 ) of the new dispersion energy donor (DED) ligand, 2,4,6-triscyclopentylphenyl with SnCl₂ afforded a mixture of the distannene {Sn(C₆H₂−2,4,6-Cyp₃)₂}₂ ( 2 ), and the cyclotristannane {Sn(C₆H₂−2,4,6-Cyp₃)₂}₃ ( 3 ). 2 is favored in solution at higher temperature (345 K or above) whereas 3 is preferred near 298 K. Van't Hoff analysis revealed the 3 to 2 conversion has a ΔH=33.36 kcal mol⁻¹ and ΔS=0.102 kcal mol⁻¹ K⁻¹, which gives a ΔG^{300 K}=+2.86 kcal mol⁻¹, showing that the conversion of 3 to 2 is an endergonic process. Computational studies show that DED stabilization in 3 is −28.5 kcal mol⁻¹ per {Sn(C₆H₂−2,4,6-Cyp₃)₂ unit, which exceeds the DED energy in 2 of −16.3 kcal mol⁻¹ per unit. The data clearly show that dispersion interactions are the main arbiter of the 3 to 2 equilibrium. Both 2 and 3 possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results).     

Effect of aluminum vacancies on the H2O2 or H2O interaction with a gamma-AlOOH surface. A solid-state DFT study

The adsorption of a single H₂O₂ or H₂O molecule on a family of periodic slab models of γ-AlOOH is studied by solid-state DFT. The single H₂O₂ or Н₂О molecule interacts with the perfect (010) slab by intermolecular hydrogen bonds (H-bonds). In the models of γ-AlOOH with oxygen and aluminum vacancies, H₂O₂ or Н₂О also forms covalent O∙∙∙Al bonds. The energies of covalent O∙∙∙Al and H-bonds are estimated by a combined approach based on simultaneous consideration of the total binding energies with BSSE correction and empirical schemes of the Н-bond energy evaluation. The O∙∙∙Al bond energy ranges from ~75 to ~156 kJ mol⁻¹. The total energy of H-bond interactions in the case of H₂O₂ exceeds that of Н₂О by ~30 kJ mol⁻¹ for all considered slab models. In contrast to Н₂О, a H₂O₂ molecule always forms two H-bonds as the proton donor. The energy of these bonds noticeably increase on defect γ-AlOOH surfaces in comparison with the perfect slab due to formation of short (strong) H-bonds by adsorbed H₂O₂.     

B−Hb←:X (X= N,O, P,S, F,Cl, Br) interactions: A density functional study

Density functional theory calculation on B H_b←:X interaction (X= N, O, P, S, F, Cl, Br) is performed. HOMO energy predicts the feasibility of such complexation. Steric and electronic effects play significant role on geometry of the complexes. Interaction energy suggests that the interaction is moderately strong (< 5.00 kcal mol⁻¹) in nature and the complexes are stable in both gas and solvent phase. Electron donating group on: X facilitates the interaction whereas electron withdrawing group impedes the same. MP2 and CCSD(T) calculations further confirm the suitability of ωB97X‐D and M06‐2X functional for studying such interactions. Dispersive interaction is the primary mode of interaction in stabilizing the complexes. ¹H NMR and IR study are also performed. Thermochemical analysis advocates exothermic nature of complexation.     

A structural and theoretical study of the intermolecular interactions in 8‐hydroxyquinolinium‐7‐carboxylate monohydrate

In the title compound, C₁₀H₇NO₃·H₂O, the zwitterionic organic molecules and the water molecules are connected by N—H...O and O—H...O hydrogen bonds to form ribbons, and π–π stacking interactions expand these ribbons into a three‐dimensional net. The energies of these hydrogen bonds adopt values typical for mildly weak interactions (3.33–7.75 kcal mol⁻¹; 1 kcal mol⁻¹ = 4.184 kJ mol⁻¹). The total π–π stacking interactions between aromatic molecules can be classified as mildly strong (energies of 15.3 and 33.9 kcal mol⁻¹), and they are made up of multiple constituent π–π interactions between six‐membered rings. The short intermolecular C—H...O contact between two zwitterionic molecules is nonbonding in character.     

An ab initio study of the structures and enthalpies of the hydrogen-bonded complexes of the acids H2O,H2S,HCN, and HCl with the anions OH−, SH−, CN−, and Cl−

Ab initio calculations at second-order Møller-Plesset perturbation theory with the 6-31 + G(d,p) basis set have been performed to determine the equilibrium structures and energies of a series of negative-ion hydrogen-bonded complexes with H₂O, H₂S, HCN, and HCl as proton donors and OH^–, SH^–, CN^–, and Cl^– as proton acceptors. The computed stabilization enthalpies of these complexes are in agreement to within the experimental error of 1 kcal mol^–1 with the gas-phase hydrogen bond enthalpies, except for HOHOH^–, in which case the difference is 1.8 kcal mol^–1. The structures of these complexes exhibit linear hydrogen bonds and directed lone pairs of electrons except for complexes with H₂O as the proton donor, in which cases the hydrogen bonds deviate slightly from linearity. All of the complexes have equilibrium structures in which the hydrogen-bonded proton is nonsymmetrically bound, although the symmetric structures of HOHOH^– and ClHCl^– are only slightly less bound than the equilibrium structures. MP2/6-31 + G(d,p) hydrogen bond energies calculated at optimized MP2/B-31 + G(d,p) and at optimized HF/6-31G(d) geometries are similar. Using HF/6-31G(d) frequencies to evaluate zero-point and thermal vibrational energies does not introduce significant error into the computed hydrogen bond enthalpies of these complexes provided that the hydrogen-bonded proton is definitely nonsymmetrically bound at both Hartree-Fock and MP2.     

Studies on the oxidation mechanism of H2S based on direct examination of the key reactions

By conducting an excimer laser photolysis (193 and 248 nm) behind shock waves, three elementary reactions important in the oxidation of H₂S have been examined, where, H, O, and S atoms have been monitored by the atomic resonance absorption spectrometry. For HS + O₂ → products (1), the rate constants evaluated by numerical simulations are summarized as: k₁ = 3.1 × 10⁻¹¹exp|-75 kJ mol⁻¹/RT| cm³molecule⁻¹s⁻¹ (T = 1400-1850 K) with an uncertainty factor of about 2. Direct measurements of the rate constants for S + O2 → SO + O (2), and SO + O₂ → SO₂ + O (3) yield k₂ = (2.5 ± 0.6) × 10⁻¹¹ exp|-(15.3 ± 2.5) kJ mol⁻¹/RT| cm³molecule⁻¹s⁻¹ (T = 980-1610 K) and, k₃ = (1.7 ± 0.9) × 10⁻¹² exp|-(34 ± 11) kJ mol⁻¹/RT| cm³molecule⁻¹s⁻¹ (T = 1130-1640 K), respectively. By summarizing these data together with the recent experimental results on the H(SINGLE BOND)S(SINGLE BOND)O reaction systems, a new kinetic model for the H₂S oxidation process is constructed. It is found that this simple reaction scheme is consistent with the experimental result on the induction time of SO₂ formation obtained by Bradley and Dobson. © 1997 John Wiley & Sons, Inc. Int J Chem Kinet 29: 57–66, 1997.     

Quantitative decomposition of resonance-assisted hydrogen bond energy in β-diketones into resonance and hydrogen bonding (π- and σ-) components using molecular tailoring and function-based approaches

Using the molecular tailoring and function-based approaches allows one to divide the energy of the O─H⋯O═C resonance-assisted hydrogen bond in a series of the β-diketones into resonance and hydrogen bonding components. The magnitude of the resonance component is assessed as about 6 kcal mol⁻¹. This value increases by ca. 1 kcal mol⁻¹ on going from the weak to strong resonance-assisted hydrogen bonding. The magnitude of the hydrogen bonding component varies in the wide range from 2 to 20 kcal mol⁻¹ depending on the structure of the β-diketone in question.