Anchoring the Absolute Proton Affinity Scale

Koopmans' theorem for acidic protons

The famous Brønsted acidity, which is relevant in many areas of experimental and synthetic chemistry, but also in biochemistry and other areas, is investigated from a new perspective. Nuclear electronic orbital methods, which explicitly account for the quantum character of selected protons, are applied. The resulting orbital energies of the proton wavefunction are interpreted and related to enthalpies of deprotonation and acid strength in analogy to the Koopmans' theorem for electrons. For a set of organic acids, we observe a correlation which indicates the validity of such a NEO-Koopmans' approach and opens up new opportunities for the computational investigation of more complex acidic systems.

Structure, energetics, and spectroscopy of the chromophores of HHe+n, H2He+n, and He+n clusters and their deuterated isotopologues

The linear molecular ions H₂He⁺, HHe+2, and He+3 are the central units (chromophores) of certain He-solvated complexes of the H₂He+n, HHe+n, and He+n families, respectively. These are complexes which do exist, according to mass-spectrometry studies, up to very high n values. Apparently, for some of the H₂He+n and He+n complexes, the linear symmetric tetratomic H₂He+2 and the diatomic He+2 cations, respectively, may also be the central units. In this study, definitive structures, relative energies, zero-point vibrational energies, and (an)harmonic vibrational fundamentals, and, in some cases, overtones and combination bands, are established mostly for the triatomic chromophores. The study is also extended to the deuterated isotopologues D₂He⁺, DHe+2, and D₂He+2. To facilitate and improve the electronic-structure computations performed, new atom-centered, fixed-exponent, Gaussian-type basis sets called MAX, with X = T(3), Q(4), P(5), and H(6), are designed for the H and He atoms. The focal-point-analysis (FPA) technique is employed to determine definitive relative energies with tight uncertainties for reactions involving the molecular ions. The FPA results determined include the 0 K proton and deuteron affinities of the ⁴He atom, 14 875(9) cm^-1 [177.95(11) kJ mol^-1] and 15 229(8) cm^-1 [182.18(10) kJ mol^-1], respectively, the dissociation energies of the He+2 → He⁺ + He, HHe+2 → HHe⁺ + He, and He+3 → He+2 + He reactions, 19 099(13) cm^-1 [228.48(16) kJ mol^-1], 3948(7) cm^-1 [47.23(8) kJ mol^-1], and 1401(12) cm^-1 [16.76(14) kJ mol^-1], respectively, the dissociation energy of the DHe+2 → DHe⁺ + He reaction, 4033(6) cm^-1 [48.25(7) kJ mol^-1], the isomerization energy between the two linear isomers of the [H, He, He]⁺ system, 3828(40) cm^-1 [45.79(48) kJ mol^-1], and the dissociation energies of the H₂He⁺ → H+2 + He and the H₂He+2 → H₂He⁺ + He reactions, 1789(4) cm^-1 [21.40(5) kJ mol^-1] and 435(6) cm^-1 [5.20(7) kJ mol^-1], respectively. The FPA estimates of the first dissociation energy of D₂He⁺ and D₂He+2 are 1986(4) cm^-1 [23.76(5) kJ mol^-1] and 474(5) cm^-1 [5.67(6) kJ mol^-1], respectively. Determining the vibrational fundamentals of the triatomic chromophores with second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) techniques, both built around the Eckart-Watson Hamiltonian, proved unusually challenging. For the species studied, VPT2 has difficulties yielding dependable results, in some cases even for the fundamentals of the H-containing molecular cations, while carefully executed VCI computations yield considerably improved spectroscopic results. In a few cases unusually large anharmonic corrections to the fundamentals, on the order of 15% of the harmonic value, have been observed.

Benchmark ab initio proton affinity of glycine

A systematic conformational search reveals three N- (amino) and eight O- (carbonyl) protonated glycine conformers with benchmark equilibrium(adiabatic) relative energies in the 0.00-7.51(0.00-7.37) and 25.91-31.61(24.45-30.28) kcal mol-1 ranges, respectively. Benchmark ab initio structures of the glycine conformers and its protonated species are obtained at the CCSD(T)-F12b/aug-cc-pVTZ level of theory and the relative energy computations consider basis-set effects up to aug-cc-pVQZ with CCSD(T)-F12b, electron correlation up to CCSDT(Q), core correlation corrections, scalar relativistic effects, and zero-point energy contributions. The best predictions for Boltzmann-averaged 0(298.15) K proton affinities and [298.15 K gas-phase basicities] of glycine are 211.00(212.43)[204.75] and 186.38(187.64)[180.21] kcal mol-1 for N- and O-protonation, respectively, in excellent agreement with experiments.

Prediction of pKa s of Late Transition-Metal Hydrides via a QM/QM Approach

Three implicit solvation models, the conductor-like polarizable continuum model (C-PCM), the conductor-like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pK_a values, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM-containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

Riddles of the structure and vibrational dynamics of HO3 resolved near the ab initio limit

The hydridotrioxygen (HO₃) radical has been investigated in many previous theoretical and experimental studies over several decades, originally because of its possible relevance to the tropospheric HO_x cycle but more recently because of its fascinating chemical bonding, geometric structure, and vibrational dynamics. We have executed new, comprehensive research on this vexing molecule via focal point analyses (FPA) to approach the ab initio limit of optimized geometric structures, relative energies, complete quartic force fields, and the entire reaction path for cis-trans isomerization. High-order coupled cluster theory was applied through the CCSDT(Q) and even CCSDTQ(P) levels, and CBS extrapolations were performed using cc-pVXZ (X = 2-6) basis sets. The cis isomer proves to be higher than trans by 0.52 kcal mol^-1, but this energetic ordering is achieved only after the CCSDT(Q) milestone is reached; the barrier for cis → trans isomerization is a minute 0.27 kcal mol^-1. The FPA central r_e(O-O) bond length of trans-HO₃ is astonishingly long (1.670 Å), consistent with the semiexperimental r_e distance we extracted from microwave rotational constants of 10 isotopologues using FPA vibration-rotation interaction constants (α_i). The D₀(HO-O₂) dissociation energy converges to a mere 2.80 ± 0.25 kcal mol^-1. Contrary to expectation for such a weakly bound system, vibrational perturbation theory performs remarkably well with the FPA anharmonic force fields, even for the torsional fundamental near 130 cm^-1. Exact numerical procedures are applied to the potential energy function for the torsional reaction path to obtain energy levels, tunneling rates, and radiative lifetimes. The cis → trans isomerization occurs via tunneling with an inherent half-life of 1.4 × 10^-11 s and 8.6 × 10^-10 s for HO₃ and DO₃, respectively, thus resolving the mystery of why the cis species has not been observed in previous experiments executed in dissipative environments that allow collisional cooling of the trans-HO₃ product. In contrast, the pure ground eigenstate of the cis species in a vacuum is predicted to have a spontaneous radiative lifetime of about 1 h and 5 days for HO₃ and DO₃, respectively.

Ab initio structure and vibration-rotation dynamics of the formyl and isoformyl cations, HCO+/HOC. J Chem Phys 2019;150:154307. [PMID: 31005073 DOI: 10.1063/1.5089718]

Accurate structure and potential energy surface of the formyl and isoformyl cation system, HCO⁺/HOC⁺, in its ground electronic state X̃ ¹Σ⁺ have been determined from ab initio calculations using the coupled-cluster approach in conjunction with the correlation-consistent basis sets up to septuple-zeta quality. Both the isomers are confirmed to be linear at equilibrium, with the total energy minimum of HOC⁺ lying 14 120 cm^-1 above that of HCO⁺ and the HCO⁺ → HOC⁺ isomerization energy barrier being 26 870 cm^-1 (in the Born-Oppenheimer approximation). The equilibrium structural parameters for HCO⁺ are estimated to be r_e(HC) = 1.0919 Å and r_e(CO) = 1.1058 Å, whereas those for HOC⁺ are estimated to be r_e(HO) = 0.9899 Å and r_e(CO) = 1.1544 Å. The vibration-rotation energy levels were predicted for various isotopologues using a variational approach and compared with the experimental data. For the spectroscopically well characterized formyl cation, the observed vibration-rotation energies and the effective rotational constants are reproduced to within about 2.3 cm^-1 and 1.7 MHz, respectively. The role of the core-electron correlation, higher-order valence-electron correlation, scalar relativistic, and adiabatic effects in determining the structure and vibration-rotation dynamics of both the isomers is discussed.

Application of multireference linearized coupled-cluster theory to atomic and molecular systems

Applications of the multireference linearized coupled-cluster single-doubles (MRLCCSD) to atomic and molecular systems have been carried out. MRLCCSD is exploited to calculate the ground-state energies of HF, H2O, NH3, CH4, N2, BF, and C2with basis sets, cc-pVDZ, cc-pVTZ and cc-pVQZ. The equilibrium bond lengths and vibration frequencies of HF, HCl, Li2, LiH, LiF, LiBr, BH, and AlF are computed with MRLCCSD and compared with the experimental data. The electron affinities of F and CH as well as the proton affinities of H2O and NH3are also calculated with MRLCCSD. These results are compared with the results produced with second-order perturbation theory, linearized coupled-cluster doubles (LCCD), coupled-cluster doubles (CCD), coupled-cluster singles and doubles (CCSD), CCSD with perturbative triples correction (CCSD(T)). It is shown that all results obtained with MRLCCSD are reliable and accurate.

Definitive thermochemistry and kinetics of the interconversions among conformers of n-butane and n-pentane

The focal-point analysis (FPA) technique is used for the definitive characterization of conformational interconversion parameters, including activation energy barriers, activation free energies, and kinetic rate coefficients at 298 K, of two n-alkanes, n-butane, and n-pentane, yielding the first complete analysis of their interconversion kinetics. The FPA implementation developed in this study is based on geometry optimizations and harmonic frequency computations carried out with density functional theory methods and single-point energy computations up to the CCSD(T) level of electronic structure theory using atom-centered Gaussian basis sets as large as cc-pV5Z. The anharmonic vibrational computations are carried out, at the MP2/6-31G* level of theory. Reflecting the convergence behavior of the Gibbs free-energy terms and the interconversion parameters, well-defined uncertainties, mostly neglected in previous theoretical studies, are provided. Finally, the effect of these uncertainties on the concentrations of the conformers of n-butane and n-pentane is examined via a global Monte-Carlo uncertainty analysis. © 2017 Wiley Periodicals, Inc.

Theoretical spectroscopic parameters for isotopic variants of HCO+ and HOC. J Chem Phys 2017;147:114111. [PMID: 28938835 DOI: 10.1063/1.4998467]

Theoretical spectroscopic parameters are derived for all isotopologues of HCO⁺ and HOC⁺ involving H, D, ¹⁶O, ¹⁷O, ¹⁸O, ¹²C, and ¹³C by means of a two-step procedure. Full-dimensional rovibrational calculations are first carried out to obtain numerically exact rovibrational energies for J = 0-15 in both parities. Effective spectroscopic constants for the vibrational ground state, ν₁, ν₂, and ν₃ are determined by fitting the calculated rovibrational energies to appropriate spectroscopic Hamiltonians. Combining our vibration-rotation corrections with the available experimental ground-state rotational constants, we also derive the new estimate for the equilibrium structure of HCO⁺, r_e(CH) = 1.091 98 Å and r_e(CO) = 1.105 62 Å, and for the equilibrium structure of HOC⁺, r_e(HO) = 0.990 48 Å and r_e(CO) = 1.154 47 Å. Regarding the spectroscopic parameters, our estimates are in excellent agreement with available experimental results for the isotopic variants of both HCO⁺ and HOC⁺: the agreement for the rotational constants B_v is within 3 MHz, for the quartic centrifugal distortion constants D_v within 1 kHz, and for the effective ℓ-doubling constants q_v within 2 MHz. We thus expect that our results can provide useful assistance in analyzing expected observations of the rare isotopologues of HCO⁺ and HOC⁺ that are not yet experimentally known.

Mechanisms of the Formation of Adenine, Guanine, and Their Analogues in UV-Irradiated Mixed NH3:H2O Molecular Ices Containing Purine

We investigated the formation mechanisms of the nucleobases adenine and guanine and the nucleobase analogues hypoxanthine, xanthine, isoguanine, and 2,6-diaminopurine in a UV-irradiated mixed 10:1 H₂O:NH₃ ice seeded with precursor purine by using ab initio and density functional theory computations. Our quantum chemical investigations suggest that a multistep reaction mechanism involving purine cation, hydroxyl and amino radicals, together with water and ammonia, explains the experimentally obtained products in an independent study. The relative abundances of these products appear to largely follow from relative thermodynamic stabilities. The key role of the purine cation is likely to be the reason why purine is not functionalized in pure ammonia ice, where cations are promptly neutralized by free electrons from NH₃ ionization. Amine group addition to purine is slightly favored over hydroxyl group attachment based on energetics, but hydroxyl is much more abundant due to higher abundance of H₂O. The amino group is preferentially attached to the 6 position, giving 6-aminopurine, that is, adenine, while the hydroxyl group is preferentially attached to the 2 position, leading to 2-hydroxypurine. A second substitution by hydroxyl or amino group occurs at either the 6 or the 2 position depending on the first substitution. Given that H₂O is far more abundant than NH₃ in the experimentally studied ices (as well as based on interstellar abundances), xanthine and isoguanine are expected to be the most abundant bi-substituted photoproducts. Key Words: Astrophysical ice-Abiotic organic synthesis-Nucleic acids-Origin of life-RNA world. Astrobiology 17, 771-785.

Enhanced Basicity of Push-Pull Nitrogen Bases in the Gas Phase

Nitrogen bases containing one or more pushing amino-group(s) directly linked to a pulling cyano, imino, or phosphoimino group, as well as those in which the pushing and pulling moieties are separated by a conjugated spacer (C═X)_n, where X is CH or N, display an exceptionally strong basicity. The n-π conjugation between the pushing and pulling groups in such systems lowers the basicity of the pushing amino-group(s) and increases the basicity of the pulling cyano, imino, or phosphoimino group. In the gas phase, most of the so-called push-pull nitrogen bases exhibit a very high basicity. This paper presents an analysis of the exceptional gas-phase basicity, mostly in terms of experimental data, in relation with structure and conjugation of various subfamilies of push-pull nitrogen bases: nitriles, azoles, azines, amidines, guanidines, vinamidines, biguanides, and phosphazenes. The strong basicity of biomolecules containing a push-pull nitrogen substructure, such as bioamines, amino acids, and peptides containing push-pull side chains, nucleobases, and their nucleosides and nucleotides, is also analyzed. Progress and perspectives of experimental determinations of GBs and PAs of highly basic compounds, termed as "superbases", are presented and benchmarked on the basis of theoretical calculations on existing or hypothetical molecules.

Sub-Doppler infrared spectroscopy and formation dynamics of triacetylene in a slit supersonic expansion

Infrared spectroscopy and formation dynamics of triacetylene are investigated in a slit jet supersonic discharge and probed with sub-Doppler resolution (≈60 MHz) on the fundamental antisymmetric CH stretch mode (ν5). The triacetylene is generated in the throat of the discharge by sequential attack of ethynyl radical with acetyelene and diacetylene: (i) HCCH → HCC + H, (ii) HCC + HCCH → HCCCCH + H, (iii) HCC + HCCCCH → HCCCCCCH + H, cooled rapidly in the slit expansion to 15 K, and probed by near shot-noise-limited absorption sensitivity with a tunable difference-frequency infrared laser. The combination of jet cooled temperatures (Trot = 15 K) and low spectral congestion permits (i) analysis of rotationally avoided crossings in the ν5 band ascribed to Coriolis interactions, as well as (ii) first detection of ν5 Π-Π hot band progressions built on the ν12 sym CC bend and definitively assigned via state-of-the-art ab initio vibration-rotation interaction parameters (αi), which make for interesting comparison with recent spectroscopic studies of Doney et al. [J. Mol. Spectrosc. 316, 54 (2015)]. The combined data provide direct evidence for significantly non-equilibrium populations in the CC bending manifold, dynamically consistent with a strongly bent radical intermediate and transition states for forming triacetylene product. The presence of intense triacetylene signals under cold, low density slit jet conditions provides support for (i) barrierless addition of HCC with HCCCCH and (ii) a high quantum yield for HCCCCCCH formation. Complete basis set calculations for energetics [CCSD(T)-f12/VnZ-f12, n = 2,3] and frequencies [CCSD(T)-f12/VdZ-f12] are presented for both radical intermediate and transition state species, predicting collision stabilization in the slit jet expansion to be competitive with unimolecular decomposition with increasing polyyne chain length.

Gas-Phase Affinity Scales for Typical Ionic Liquid Moieties Determined by using Cooks' Kinetic Method

Gas-phase affinity studies based on cations and anions commonly present in ionic liquid structures, give quantitative information about the magnitude of the interactions holding the two species together when ILs are formed. They also provide clues on how these interactions depend on the nature of the cationic and anionic moieties. In the present work, mass spectrometric experiments, performed using electrospray ionization quadrupole ion-trap and Fourier transform ion cyclotron resonance mass spectrometry, were used to obtain two affinity scales by Cooks' kinetic method: one scale for the various cations for the bis(trifluoromethylsulfonyl)imide anion, [NTf2 ](-) , and another for the different anions for the 1-butyl-3-methylimidazolium cation, [C4 mim](+) . The obtained results are compared with previously reported data and discussed in terms of the structural characteristics of the different cationic and anionic species.

Active Thermochemical Tables: Sequential Bond Dissociation Enthalpies of Methane, Ethane, and Methanol and the Related Thermochemistry

Active Thermochemical Tables (ATcT) thermochemistry for the sequential bond dissociations of methane, ethane, and methanol systems were obtained by analyzing and solving a very large thermochemical network (TN). Values for all possible C-H, C-C, C-O, and O-H bond dissociation enthalpies at 298.15 K (BDE298) and bond dissociation energies at 0 K (D0) are presented. The corresponding ATcT standard gas-phase enthalpies of formation of the resulting CHn, n = 4-0 species (methane, methyl, methylene, methylidyne, and carbon atom), C2Hn, n = 6-0 species (ethane, ethyl, ethylene, ethylidene, vinyl, ethylidyne, acetylene, vinylidene, ethynyl, and ethynylene), and COHn, n = 4-0 species (methanol, hydroxymethyl, methoxy, formaldehyde, hydroxymethylene, formyl, isoformyl, and carbon monoxide) are also presented. The ATcT thermochemistry of carbon dioxide, water, hydroxyl, and carbon, oxygen, and hydrogen atoms is also included, together with the sequential BDEs of CO2 and H2O. The provenances of the ATcT enthalpies of formation, which are quite distributed and involve a large number of relevant determinations, are analyzed by variance decomposition and discussed in terms of principal contributions. The underlying reasons for periodic appearances of remarkably low and/or unusually high BDEs, alternating along the dissociation sequences, are analyzed and quantitatively rationalized. The present ATcT results are the most accurate thermochemical values currently available for these species.

On the protonation of water

Dissociative photoionization onsets of water and water dimer, measured by Imaging Photoelectron Photoion Coincidence (iPEPICO) Spectroscopy, are used in a floating thermochemical cycle to determine the proton affinity of water with unprecedented accuracy, as confirmed by state-of-the-artab initioquantum-chemical calculations.

Effects of collisional and vibrational velocity on proton and deuteron transfer in the reaction of HOD+ with CO

Reaction of HOD(+) with CO was studied over the collision energy (E(col)) range between 0.18 and 2.87 eV, for HOD(+) in its ground state and with one quantum in each of its vibrational modes: (001)--predominantly OH stretch; (010)--bend, and (100)--predominately OD stretch. In addition to integral cross sections, product recoil velocity distributions were also measured for each initial condition. The dominant reactions are near-thermoneutral proton and deuteron transfer, generating HCO(+) and DCO(+) product ions by a predominantly direct mechanism. The HCO(+) and DCO(+) channels occur with a combined efficiency of 76% for ground state HOD(+) at our lowest E(col), increasing to 94% for E(col) around 0.33 eV, then falling at high E(col) to ~40%. The HCO(+) and DCO(+) channels have a complicated dependence on the HOD(+) vibrational state. Excitation of the OH or OD stretch modes enhances H(+) or D(+) transfer, respectively, and inhibits D(+) or H(+) transfer. Bend excitation preferentially enhances H(+) transfer, with no effect on D(+) transfer. There is no coupling of energy initially in any HOD(+) vibrational mode to recoil velocity of either of the product ions, even at low E(col) where vibrational excitation doubles or triples the energy available to products. The results suggest that transfer of H or D atoms is enhanced if the atom in question has a high vibrational velocity, and that this effect offsets what is otherwise a general inhibition of reactivity by added energy. HOCO(+) + D and DOCO(+) + H products are also observed, but as minor channels despite being barrierless and exoergic. These channels appear to be complex mediated at low E(col), essentially vanish at intermediate E(col), then reappear with a direct reaction mechanism at high E(col).

Benchmarking Experimental and Computational Thermochemical Data: A Case Study of the Butane Conformers

Due to its crucial importance, numerous studies have been conducted to determine the enthalpy difference between the conformers of butane. However, it is shown here that the most reliable experimental values are biased due to the statistical model utilized during the evaluation of the raw experimental data. In this study, using the appropriate statistical model, both the experimental expectation values and the associated uncertainties are revised. For the 133-196 and 223-297 K temperature ranges, 668 ± 20 and 653 ± 125 cal mol(-1), respectively, are recommended as reference values. Furthermore, to show that present-day quantum chemistry is a favorable alternative to experimental techniques in the determination of enthalpy differences of conformers, a focal-point analysis, based on coupled-cluster electronic structure computations, has been performed that included contributions of up to perturbative quadruple excitations as well as small correction terms beyond the Born-Oppenheimer and nonrelativistic approximations. For the 133-196 and 223-297 K temperature ranges, in exceptional agreement with the corresponding revised experimental data, our computations yielded 668 ± 3 and 650 ± 6 cal mol(-1), respectively. The most reliable enthalpy difference values for 0 and 298.15 K are also provided by the computational approach, 680.9 ± 2.5 and 647.4 ± 7.0 cal mol(-1), respectively.

Differentiation of isomeric hydroxypyridine N-oxides using metal complexation and electrospray ionization mass spectrometry

Differentiation between two isomers of hydroxypyridine N-oxide according to the metal cation adducts generated by electrospray ionization (ESI) was investigated for different metal cations, namely Mg(II), Al(III), Ca(II), Sc(III), Fe(III), Co(II), Ni(II), Cu(II), Zn(II), Ga(III), besides the diatomic cation VO(IV). Protonated molecules of the isomeric hydroxypyridine N-oxides as well as the singly/doubly charged adducts formed from neutral or deprotonated ligands and a doubly/triply charged cation were produced in the gas phase using ESI, recording mass spectra with different metal ions for each isomer. While complex formation was successful for 2-hydroxypyridine N-oxide with trivalent ions, in the case of 3-hydroxypyridine N-oxide, only peaks related to the protonated molecule were present. On the other hand, divalent cations formed specific species for each isomer, giving characteristic spectra in every case. Hence, differentiation was possible irrespective of the metal cation utilized. In addition, quantum chemical calculations at the B3LYP/6-31+G(d,p) level of theory were performed in order to gain insight into the different complexation of calcium(II) with the isomers of hydroxypyridine N-oxide. The relative stability in the gas phase of the neutral complexes of calcium made up of two ligands, as well as the singly charged and doubly charged complexes, was investigated. The results of these calculations improved the understanding of the differences observed in the mass spectra obtained for each isomer.

Accurate benchmark calculations on the gas-phase basicities of small molecules

Accurate benchmark calculations of gas-phase basicities of small molecules are presented and compared with available experimental results. The optimized geometries and thermochemical analyses were obtained from MP2/aug-cc-pVTZ calculations. Two different ab initio electron-correlated methods MP2 and CCSD(T) were employed for subsequent gas-phase basicity calculations, and the single-point energies were extrapolated to the complete basis set (CBS) limit. The overall accuracy for different ab initio methods is compared, and the accuracy in descending order is CCSD(T)_CBS > CCSD(T)/aug-cc-pVDZ > (MP2/aug-cc-pVQZ approximately MP2_CBS) > HF/aug-cc-pVQZ. The best root-mean-squared-error obtained was 1.0 kcal mol(-1) at the CCSD(T)_CBS//MP2/aug-cc-pVTZ level for a test set of 41 molecules. Clearly, accurate calculations for the electron correlation energy are important for the theoretical prediction of molecular gas-phase basicities. However, conformational effects were also found to be relevant in several instances when more complicated molecules were examined.