The keto-enol equilibrium in substituted acetaldehydes: focal-point analysis and ab initio limit

High-level ab initio electronic structure calculations up to the CCSD(T) theory level, including extrapolations to the complete basis set (CBS) limit, resulted in high precision energetics of the tautomeric equilibrium in 2-substituted acetaldehydes (XH₂C-CHO). The CCSD(T)/CBS relative energies of the tautomers were estimated using CCSD(T)/aug-cc-pVTZ, MP3/aug-cc-pVQZ, and MP2/aug-cc-pV5Z calculations with MP2/aug-cc-pVTZ geometries. The relative enol (XHC?=?CHOH) stabilities (ΔE _{e,CCSD(T)/CBS}) were found to be 5.98?±?0.17, ?1.67?±?0.82, 7.64?±?0.21, 8.39?±?0.31, 2.82?±?0.52, 10.27?±?0.39, 9.12?±?0.18, 5.47?±?0.53, 7.50?±?0.43, 10.12?±?0.51, 8.49?±?0.33, and 6.19?±?0.18?kcal?mol^?1 for X?=?BeH, BH₂, CH₃, Cl, CN, F, H, NC, NH₂, OCH₃, OH, and SH, respectively. Inconsistencies between the results of complex/composite energy computations methods Gn/CBS (G2, G3, CBS-4M, and CBS-QB3) and high-level ab initio methods (CCSD(T)/CBS and MP2/CBS) were found. DFT/aug-cc-pVTZ results with B3LYP, PBE0 (PBE1PBE), TPSS, and BMK density functionals were close to the CCSD(T)/CBS levels (MAD?=?1.04?kcal?mol^?1).     

The gas phase oxidation of HCOOH by Cl and NH2 radicals. Proton coupled electron transfer versus hydrogen atom transfer*

ABSTRACT

The reaction of formic acid (HCOOH) with chlorine atom and amidogen radical (NH₂) have been investigated using high level theoretical methods such BH&HLYP, MP2, QCISD, and CCSD(T) with the 6–311?+?G(2df,2p), aug-cc-pVTZ, aug-cc-pVQZ and extrapolation to CBS basis sets. The abstraction of the acidic and formyl hydrogen atoms of the acid by the two radicals has been considered, and the different reactions proceed either by a proton coupled electron transfer (pcet) and hydrogen atom transfer (hat) mechanisms. Our calculated rate constant at 298?K for the reaction with Cl is 1.14?×?10^?13?cm³?molecule^?1?s^?1 in good agreement with the experimental value 1.8?±?0.12/2.0?×?10^?13?cm³?molecule^?1?s^?1 and the reaction proceeds exclusively by abstraction of the formyl hydrogen atom, via hat mechanism, producing HOCO+ClH. The calculated rate constant, at 298?K, for the reaction with NH₂ is 1.71?×?10^?15?cm³?molecule^?1?s^?1, and the reaction goes through the abstraction of the acidic hydrogen atom, via a pcet mechanism, leading to the formation of HCOO+NH₃.     

Atmospheric oxidation chemistry of 1–methoxy 2–propyl acetate initiated by OH radicals: kinetics and mechanisms

Kinetics and mechanism of the gas-phase reaction of CH₃C(O)OCH(CH₃)CH₂OCH₃ (MPA) with OH radicals in the presence of O₂ and NO have been investigated theoretically by performing a high and reliable level of theory, viz., CCSD(T)/6-311?+?G(d,p)//BH&HLYP/6-311++G(d,p)?+?0.9335×ZPE. The calculations predict that the H-abstraction from the ?CH₂?O? position of MPA is the most facile channel, which leads to the formation of the corresponding alkoxy radicals CH₃C(O)OCH(CH₃)C(O ?)HOCH₃ under atmospheric conditions. This activated radicals CH₃C(O)OCH(CH₃)C(O ?)HOCH₃ will undergo further rearrangement, fragmentation and oxidative reactions and predominantly leads to the formation of various products (methyl formate HC(O)OCH₃ and acetic anhydride CH₃C(O)OC(O)CH₃). In the presence of water, acetic anhydride can convert into acetic acid CH₃C(O)OH via the hydrolysis reaction. The calculated total rate constants over the temperature range 263–372?K are used to derive a negative activation energy (E_a= ?5.88 kJ/mol) and an pre-exponential factor (A?=?1.78×10^?12 cm³ molecule^?1 s^?1). The obtained Arrhenius parameters presented here are in strong agreement with the experimental values. Moreover, the temperature dependence of the total rate constant over a temperature range of 263?1000?K can be described by k?=?5.60 × 10^?14×(T/298?K)^3.4×exp(1725.7?K/T) cm³ molecule^?1 s^?1.     

Thermochemical properties of small open-shell systems: experimental and high-level ab initio results for NH2 and

The first adiabatic ionization energy and the first singlet–triplet splitting of the amidogen radical (NH₂) have been determined by high-level ab initio quantum chemistry based on the coupled-cluster approach (90?041 and 10?319?cm^?1, respectively) and by high-resolution pulsed-field-ionization zero-kinetic-energy (PFI-ZEKE) photoelectron spectroscopy (90 083.8 ± 1.0 and 10 222.0 ± 1.3?cm^?1, respectively). A comparison between the theoretical and experimental values demonstrates the predictive powers of high-level ab initio theory in the derivation of the thermochemical properties of small molecular systems. The absolute accuracy of better than 100?cm^?1 alleviates the experimental search for the relevant spectral features.     

Theoretical study of rearrangements in water dimer and trimer

The global minimum and transition states for the acceptor-tunnelling, donor-acceptor interchange and bifurcation tunnelling rearrangements of the water dimer, and the single-flip, bifurcation and concerted proton transfer processes in the water trimer have been reinvestigated. Our analysis of the tunnelling splittings and spectroscopy is based on ab initio calculations at the computational level of second-order M?ller-Plesset (MP2) theory with basis sets of aug-cc-pVXZ quality (X = D, T, Q for the dimer; X = D, T for the trimer). In both water dimer and trimer, the binding energy, barrier heights, intermonomer distances, and harmonic frequencies converge smoothly as the size of the basis set increases. In the water dimer, the binding energy was evaluated as 5.09kcal mol^?1, while the activation energies are 0.52 (acceptor-tunnelling) 0.79 (donor-acceptor interchange), and 1.94 kcal mol^?1 (bifurcation tunnelling) at the MP2/aug-cc-pVQZ level. In the water trimer, the binding energy was evaluated as 16.29 kcal mol^?1, while the activation energies are 0.28 (single-flip), 2.34 (bifurcation), and 26.36 (proton transfer) kcal mol^?1 at the MP2/aug-cc-pVTZ level.     

Mass spectrometric and ab initio study of the interaction between 9-methylguanine and amino acid amide group

The temperature dependent field ionization mass spectrometry method combined with ab initio calculations was used to determine the interaction energies and the structures of 9-methylguanine-acrylamide dimers. Acrylamide mimics the side chain amide group of the natural amino acids asparagine and glutamine. The experimental enthalpy of the dimer formation derived from the van't Hoff plot is ?59.5 ± 3.8 kJ mol^?1. The value is higher than interaction energies between acrylamide and other nucleic acid bases which were determined to be ?57.0 for 1-methylcytosine, ?52.0 for 9-methyladenine, and ?40.6 kJ mol^?1 for 1-methyl-uracil. In total, eight hydrogen bonded dimers formed by the three lowest energy 9-methylguanine tautomers and acrylamide were found in the quantum chemical calculations performed at the DFT/B3LYP/6-31++G?? and MP2/6-31++G?? levels of theory. The relative stability and the interaction energies of the dimers were calculated accounting for the basis set superposition error and the zero-point vibrational energy correction. The lowest energy dimer found in the calculations is formed by acrylamide (Ac) with the keto tautomer of 9-methylguanine (Gk). It is stabilized by two intermolecular H bonds, C6=O(Gk) · · · H—N(Ac) and Nl—H(Gk) · · ·O(Ac), and it is more stable than the second lowest energy dimer by ≈ 25 kJ mol^?1. The calculated interaction energies of the lowest energy 9-methylguanine-acrylamide dimer are ?65.0 kJ mol^?1 and ?67.7 kJ mol^?1 at the MP2 and DFT levels of theory, respectively. The experimental enthalpy of the dimer formation is in good agreement with both the calculated interaction energies of the GkAc dimer and much higher than the interaction energies calculated for all other 9-methylguanine-acrylamide dimers. This proved that only one dimer was present in the experimental samples. To verify whether acrylamide is a good model of the amino acid-amide group, we performed direct calculations of the 9-methylguanine-glutamine dimers at the same levels of theory as used for the complexes involving acrylamide. The interaction energies found for the lowest energy 9-methylguanine-glutamine dimer are ?65.1 kJ mon^?1 (MP2/6-31++G??) and ?66.2 kJ mol^?1 (DFT/B3LYP/6-31++G??) and these values are very close (within 0.5 kJ mol^?1) to the interaction energies obtained for the 9-methylguanine-acrylamide dimers.     

The gas-phase pyrolysis of methyl azidoformate in the absence and presence of water: a theoretical study

The potential energy surface profiles for the gas-phase pyrolysis of methyl azidoformate (MA, CH₃OC(O)N₃) in the absence and presence of one water molecule have been investigated by ab initio methods at CCSD(T)/6-311++G(2df,2pd)//MP2(full)/6-311++G(d,p)+0.95×ZPE levels of theory. Three types of mechanisms are discussed for the gas-phase decomposition of CH₃OC(O)N₃. Ab initio calculations show that a four-membered-ring intermediate can be formed by the stepwise routes. The resulting intermediate can undergo two competitive decomposition channels to generate the major products CO₂?+?CH₂?=?NH and HNCO?+?HC(O)H. The calculated results are in qualitative agreement with the observed experimental data. However, CH₃ONCO can be produced from the Curtius-type rearrangement route. This is an intriguing finding in this study. Moreover, the effect of one water molecule on the gas-phase pyrolysis of MA has been also explored. We find that the relative energy of the hydrated transition states is effectively lowered when water is added to the reaction. However, the estimated rate constant at 625?K for the naked reaction is about 30 times faster than the reaction with water. Thus, a single water molecule cannot play an important role in the thermal decomposition of MA.     

Ab initio determination of the torsional spectra of acetic acid

The torsional potential energy surface and the favourite geometries of acetic acid are determined with MP4/cc-p VTZ ab initio calculations. The molecule shows two planar trans and cis conformers whose energy difference is 1882.7 cm^?1. Both minimum energy geometries are separated by a barrier of 4432.1 cm^?1. The most stable trans-conformer shows a quite low methyl torsion barrier of 169.8 cm^?1. The roto-torsional energy levels have been calculated up to J = 10. The two torsional fundamental frequencies of the trans-conformer, the methyl and the OH torsion are 82.857 (A²) and 77.050cm^?1 (E) and 568.532 (A²) and 568.418cm^?1 (E). The V₃ barrier causes a splitting of 0.315cm^?1 in the ground vibrational state where the quartic centrifugal distortion constants have been predicted to be D_J = 90.4kHz, D_JK = ?301.5kHz and D_K = 165.4kHz. Finally, the far-infrared spectra of two isotopomers have been simulated from ab initio calculations.     

A computational study of N2···HHeF: comparison with N2···HArF and N2···HKrF

A weakly bound linear complex of N₂ and HHeF was found to be stable with respect to the constituent monomers by ab initio calculations at various levels of theory (MP2, MP3, MP4(SDQ) and QCISD) using a 6-311++G(2d,2p) basis set. The complex N₂···HHeF was found to have a zero-point vibrational energy corrected binding energy of 14.5?kJ mol^?1 (QCISD) and exhibits a large harmonic vibrational frequency blue shift of 375?cm^?1 for the He–H stretching vibration mode, with a diminished infrared intensity for this mode on formation of the complex. The frequency shift for this mode was also found to be very sensitive to the level of theory employed for the calculation, and is rationalized by considering intermolecular electrostatic and charge-transfer effects. The results for N₂···HHeF are compared with corresponding results for the related complexes N₂···HArF and N₂···HKrF, both of which contain the same proton acceptor molecule.     

Quantum chemical studies of the solvation of the hydroxide ion

Abstract

To understand and model the solvation of the hydroxide ion, OH(H₂O)^? _n clusters, n = 1?5, are studied using ab initio quantum chemical techniques, largely at the MP2 level of theory using a double zeta plus polarization functions basis extended by diffuse functions. Energies and vibrational frequencies, together with thermodynamic quantities such as enthalpies, entropies and Gibbs free energies, are computed. This permits comparison with experimental estimates of the successive thermodynamic changes associated with the reaction OH(H₂O)^? _n + H₂O → OH(H₂O)^? _n+1. The theoretical values are in good agreement with experiment. The free energy of hydration of OH^? is modelled by a composite discrete-continuum method where the effects of the first hydration shell (n = 3) are obtained from the gas phase cluster calculation, while the long-range effects are modelled using self consistent reaction field theory, namely by calculating the solvation energy of OH(H₂O)^? _n in a dielectric continuum. The best estimate of the solvation (free) energy at 298 K is ?84·5 kcal mol^?1, compared to the experimental value of ?102·8 kcal mol^?1.     

Kinetic and mechanistic study of the atmospheric reaction of MBO331 with Cl atoms

The present work deals with the reaction of 3-methyl-3-buten-1-ol (MBO331) with Cl atoms, which has been investigated by gas chromatography with flame ionization detection (GC-FID) at atmospheric pressure in N₂ or air, using the relative rate technique. The rate constant reaction at 298?±?1?K was found to be (5.01?±?0.70)?×?10^?10?cm³ molecule^?1?s^?1, using cyclohexane, octane and 1-butene as a reference compounds. The temperature dependence for the reaction was studied within the 298?333?K range. Additionally, a product identification under atmospheric conditions has been performed for the first time by GC-MS, with 3-methyl-3-butenal, methacrolein and chloroacetone being observed as degradation products. A theoretical study on the reaction at the QCISD(T)/6-311G**//MP2/6-311G** level was also carried out to obtain more information on the mechanism. From the theoretical study it can be predicted that Cl addition to the double bond proceeds through lower energy barriers than H-abstraction pathways and therefore is energetically favoured. Finally, atmospheric implications of the results obtained are discussed.     

Explicitly correlated benchmark calculations on C8H8 isomer energy separations: how accurate are DFT,double-hybrid,and composite ab initio procedures?

Accurate isomerization energies are obtained for a set of 45 C₈H₈ isomers by means of the high-level, ab initio W1-F12 thermochemical protocol. The 45 isomers involve a range of hydrocarbon functional groups, including (linear and cyclic) polyacetylene, polyyne, and cumulene moieties, as well as aromatic, anti-aromatic, and highly-strained rings. Performance of a variety of DFT functionals for the isomerization energies is evaluated. This proves to be a challenging test: only six of the 56 tested functionals attain root mean square deviations (RMSDs) below 3?kcal?mol^?1 (the performance of MP2), namely: 2.9 (B972-D), 2.8 (PW6B95), 2.7 (B3PW91-D), 2.2 (PWPB95-D3), 2.1 (ωB97X-D), and 1.2 (DSD-PBEP86) kcal?mol^?1. Isomers involving highly-strained fused rings or long cumulenic chains provide a ‘torture test’ for most functionals. Finally, we evaluate the performance of composite procedures (e.g. G4, G4(MP2), CBS-QB3, and CBS-APNO), as well as that of standard ab initio procedures (e.g. MP2, SCS-MP2, MP4, CCSD, and SCS-CCSD). Both connected triples and post-MP4 singles and doubles are important for accurate results. SCS-MP2 actually outperforms MP4(SDQ) for this problem, while SCS-MP3 yields similar performance as CCSD and slightly bests MP4. All the tested empirical composite procedures show excellent performance with RMSDs below 1?kcal?mol^?1.     

A refined estimate of the bond length of methane

The atmospheric reaction of H₂S with Cl was investigated using high level ab initio calculations and Canonical Variational Transition State Theory (CVTST). The adduct formation step is the dynamical bottleneck, and the rate constant was calculated to be 1.2 × 10^?9 cm³ molecule^?1 s^?1, which is around ten times greater than the upper experimental value. Additional ab initio classical trajectory calculations show that the adduct formed in the initial collision can easily dissociate, recrossing the variational transition state. The stabilization of this species depends on the vibrational excitation of H₂S molecule, which requires an almost collinear SH-Cl collision. These dynamical effects provide an explanation for the substantial error in the rate constant obtained using CVTST.     

Direct ab initio dynamics calculations of the reaction rate for the hydrogen abstraction reaction of NCO with CH4 and C2H6

The kinetics of the hydrogen abstraction reactions NCO + CH₄ (R1) and NCO + C₂H₆ (R2) have been studied over a wide temperature range. The minimum energy paths (MEPs) were calculated at the MP2/cc-pVDZ level and single-point calculations were refined at the G3MP2 level. The rate constants for the title reactions were calculated using canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions. The fitted three-parameter formulae are k ₁ = 2.52 × 10^?22 T ^3.46 exp(2466/T) and k ₂ = 9.8 × 10^?22 T ^3.2 exp(411.8/T) cm³ molecule^?1 s^?1 for (R1) and (R2), respectively. The calculated rate constants were found to be in good agreement with the available experimental data. Deuterium kinetic isotope effects were also investigated. Both reactions show a significant kinetic isotope effect in the low-temperature range.     

DFT study on abstraction reaction mechanism of oh radical with 2‐methoxyphenol

Reaction mechanism of 2‐methoxyphenol (2MP) (guaiacol) with OH radical has been performed using density functional theory methods BH&HLYP and MPW1K method with 6‐311++G(d,p) basis set. Single‐point energy calculations were done using CCSD(T)/6‐311++G(d,p). The theoretical results reveal that the hydrogen abstraction from methoxy group is found to be the dominant reaction channel with an energy barrier of 9.31 kcal/mol. Also, time‐dependent density functional theory calculations have been performed using BH&HLYP/6‐311++G(d,p) level of theory, and the results reveal that the reactions occur in ground state than the excited state. The results of reaction force profile indicate that structural rearrangements are most influential with high percentage than the relaxation process. The calculated theoretical rate constants (12.19 × 10^?11 cm³ molecule^?1 s^?1) are in good agreement with the experimental rate constant. The atmospheric lifetime of 2‐methoxyphenol with respect to OH radicals is 2.27 hours, which implies that OH radical plays an important role in the degradation of 2MP. The Wiberg bond index of the abstraction reaction reveals that the bond order is concerted, partially synchronic. The reactant‐like transition state satisfies Hammond postulate, which eventually results in an exothermic reaction, and the product‐like transition state reveals in endothermic nature.     

Enthalpy difference between conformations of normal alkanes: effects of basis set and chain length on intramolecular basis set superposition error

The quantum chemistry of conformation equilibrium is a field where great accuracy (better than 100?cal?mol^?1) is needed because the energy difference between molecular conformers rarely exceeds 1000–3000?cal?mol^?1. The conformation equilibrium of straight-chain (normal) alkanes is of particular interest and importance for modern chemistry. In this paper, an extra error source for high-quality ab initio (first principles) and DFT calculations of the conformation equilibrium of normal alkanes, namely the intramolecular basis set superposition error (BSSE), is discussed. In contrast to out-of-plane vibrations in benzene molecules, diffuse functions on carbon and hydrogen atoms were found to greatly reduce the relative BSSE of n-alkanes. The corrections due to the intramolecular BSSE were found to be almost identical for the MP2, MP4, and CCSD(T) levels of theory. Their cancelation is expected when CCSD(T)/CBS (CBS, complete basis set) energies are evaluated by addition schemes. For larger normal alkanes (N?>?12), the magnitude of the BSSE correction was found to be up to three times larger than the relative stability of the conformer; in this case, the basis set superposition error led to a two orders of magnitude difference in conformer abundance. No error cancelation due to the basis set superposition was found. A comparison with amino acid, peptide, and protein data was provided.     

Theoretical investigations of the gas phase reaction of limonene (C10H16) with OH radical

The rate coefficients of hydroxyl radical (OH) reaction with limonene were computed using canonical variational transition state theory with small-curvature tunnelling between 275 and 400 K. The geometries and frequencies of all the stationary points are calculated using hybrid density functional theory methods M06-2X and MPWB1K with 6-31+G(d,p), 6-311++G(d,p), and 6-311+G(2df,2p) basis sets. Both addition and abstraction channels of the title reaction were explored. The rate coefficients obtained over the temperature range of 275–400 K were used to derive the Arrhenius expressions: k(T) = 4.06×10^?34 T^7.07 exp[4515/T] and k(T) = 7.37×10^?25 T^3.9 exp[3169/T] cm³ molecule^?1 s^?1 at M06-2X/6-311+G(2df,2p) and MPWB1K/6-311+G(2df,2p) levels of theory, respectively. Kinetic study indicated that addition reactions are major contributors to the total reaction in the studied temperature range. The atmospheric lifetime (τ) of limonene due to its reactions with various tropospheric oxidants was calculated and concluded that limonene is lost in the atmosphere within a few hours after it is released. The ozone production potential of limonene was computed to be (14–18) ppm, which indicated that degradation of limonene would lead to a significant amount of ozone production in the troposphere.     

Theoretical study of anharmonic resonances in HBS+

A theoretical study of spin-rovibronic levels of HBS⁺ in the × ²II ground state is reported. Vibronic energy levels are calculated variationally using highly accurate full dimensional ab initio potential energy surfaces. All experimental levels are reproduced within 15 cm^?1, and predictions of rovibronic levels up to 4000 cm^?1 for four isotopomers are provided. Anharmonic resonances induced by the Renner-Teller effect and strong spin-orbit coupling are analysed in detail.     

Conformational stability from Raman spectra,r0 structural parameters,and vibrational assignment of methylcyclohexane

The Raman spectra (3500–50 cm⁻¹) of the liquid and solid methylcyclohexane and the infrared spectra of the gas and solid methylcyclohexane have been recorded. The Raman band at 754 cm⁻¹ in the liquid has been confidently assigned to the less stable axial conformer and its intensity was recorded as a function of temperature from 25 to −95 °C. By the utilization of 15 different temperatures, the enthalpy difference between the more stable chair‐equatorial conformer and the chair‐axial form was determined to be 712 ± 71 cm⁻¹ (8.50 ± 0.84 kJ/mol). The ab initio predicted value of 710 cm⁻¹ (8.50 kJ/mol) from the MP2(full)/6‐311G(2d,2p) calculations with and without diffuse functions is in excellent agreement. The harmonic force fields, infrared intensities, Raman activities, depolarization ratios, and vibrational wavenumbers have been obtained for both conformers from MP2(full)/6‐31G(d) ab initio calculations. With two scaling factors of 0.88 for the C‐H stretches and 0.9 for the remaining ones, the fundamental wavenumbers have been predicted and along with the depolarization values and infrared band contours (B‐type for A″ modes) a complete vibrational assignment has been made for the chair‐equatorial conformer. Predicted r₀ structural parameters have been provided from adjusted parameters from ab initio MP2(full)/6‐311+G(d,p) calculations. The results are discussed and compared with the corresponding properties of some similar molecules. Copyright © 2009 John Wiley & Sons, Ltd.     

A high level computational study of the CH4/CF4 dimer: how does it compare with the CH4/CH4 and CF4/CF4 dimers?

The interaction within the methane–methane (CH₄/CH₄), perfluoromethane–perfluoromethane (CF₄/CF₄) methane–perfluoromethane dimers (CH₄/CF₄) was calculated using the Hartree–Fock (HF) method, multiple orders of Møller–Plesset perturbation theory [MP2, MP3, MP4(DQ), MP4(SDQ), MP4(SDTQ)], and coupled cluster theory [CCSD, CCSD(T)], as well as the PW91, B97D, and M06-2X density functional theory (DFT) functionals. The basis sets of Dunning and coworkers (aug-cc-pVxZ, x?=?D, T, Q), Krishnan and coworkers [6-311++G(d,p), 6-311++G(2d,2p)], and Tsuzuki and coworkers [aug(df, pd)-6-311G(d,p)] were used. Basis set superposition error (BSSE) was corrected via the counterpoise method in all cases. Interaction energies obtained with the MP2 method do not fit with the experimental finding that the methane–perfluoromethane system phase separates at 94.5?K. It was not until the CCSD(T) method was considered that the interaction energy of the methane–perfluoromethane dimer (?0.69?kcal?mol^?1) was found to be intermediate between the methane (?0.51?kcal?mol^?1) and perfluoromethane (?0.78?kcal?mol^?1) dimers. This suggests that a perfluoromethane molecule interacts preferentially with another perfluoromethane (by about 0.09?kcal?mol^?1) than with a methane molecule. At temperatures much lower than the CH₄/CF₄ critical solution temperature of 94.5?K, this energy difference becomes significant and leads perfluoromethane molecules to associate with themselves, forming a phase separation. The DFT functionals yielded erratic results for the three dimers. Further development of DFT is needed in order to model dispersion interactions in hydrocarbon/perfluorocarbon systems.