Joseph R. Lane

Are Bond Critical Points Really Critical for Hydrogen Bonding

Atoms in Molecules (AIM) theory is routinely used to assess hydrogen bond formation; however its stringent criteria controversially exclude some systems that otherwise appear to exhibit weak hydrogen bonds. We show that a regional analysis of the reduced density gradient, as provided by the recently introduced Non-Covalent Interactions (NCI) index, transcends AIM theory to deliver a chemically intuitive description of hydrogen bonding for a series of 1,n-alkanediols. This regional definition of interactions overcomes the known caveat of only analyzing electron density critical points. In other words, the NCI approach is a simple and elegant generalization of the bond critical point approach, which raises the title question. Namely, is it the presence of an electron density bond critical point that defines a hydrogen bond or the general topology in the region surrounding it?

A computational study of the oxidation of SO2 to SO3 by gas-phase organic oxidants.

We have studied the oxidation of SO(2) to SO(3) by four peroxyradicals and two carbonyl oxides (Criegee intermediates) using both density functional theory, B3LYP, and explicitly correlated coupled cluster theory, CCSD(T)-F12. All the studied peroxyradicals react very slowly with SO(2) due to energy barriers (activation energies) of around 10 kcal/mol or more. We find that water molecules are not able to catalyze these reactions. The reaction of stabilized Criegee intermediates with SO(2) is predicted to be fast, as the transition states for these oxidation reactions are below the free reactants in energy. The atmospheric relevance of these reactions depends on the lifetimes of the Criegee intermediates, which, at present, is highly uncertain.

CCSDTQ Optimized Geometry of Water Dimer.

The equilibrium geometry of the lowest energy structure of water dimer [(H2O)2] has been investigated using coupled cluster theory. A hierarchy of conventional coupled cluster methods is utilized up to singles doubles triples and quadruples excitations (CCSDTQ). The geometry of (H2O)2 is also optimized using the explicitly correlated coupled cluster singles doubles and perturbative triples [CCSD(T)-F12b] method. Overall, we find that the effect of including excitations beyond CCSD(T) is smaller than inclusion of core-valence correlation and comparable to scalar-relativistic and adiabatic effects.

Identification of the dimethylamine-trimethylamine complex in the gas phase

We have identified the dimethylamine-trimethylamine complex (DMA-TMA) at room temperature in the gas phase. The Fourier transform infrared (FTIR) spectrum of DMA-TMA in the NH-stretching fundamental region was obtained by spectral subtraction of spectra of each monomer. Explicitly correlated coupled cluster calculations were used to determine the minimum energy structure and interaction energy of DMA-TMA. Frequencies and intensities of NH-stretching transitions were also calculated at this level of theory with an anharmonic oscillator local mode model. The fundamental NH-stretching intensity in DMA-TMA is calculated to be approximately 700 times larger than that of the DMA monomer. The measured and calculated intensity is used to determine a room temperature equilibrium constant of DMA-TMA of 1.7 × 10(-3) atm(-1) at 298 K.

Explicit correlation and intermolecular interactions: Investigating carbon dioxide complexes with the CCSD(T)-F12 method

We have optimized the lowest energy structures and calculated interaction energies for the CO(2)-Ar, CO(2)-N(2), CO(2)-CO, CO(2)-H(2)O, and CO(2)-NH(3) dimers with the recently developed explicitly correlated coupled cluster singles doubles and perturbative triples [CCSD(T)]-F12 methods and the associated VXZ-F12 (where X = D,T,Q) basis sets. For a given cardinal number, we find that results obtained with the CCSD(T)-F12 methods are much closer to the CCSD(T) complete basis set limit than the conventional CCSD(T) results. The relatively modest increase in the computational cost between explicit and conventional CCSD(T) is more than compensated for by the impressive accuracy of the CCSD(T)-F12 method. We recommend use of the CCSD(T)-F12 methods in combination with the VXZ-F12 basis sets for the accurate determination of equilibrium geometries and interaction energies of weakly bound electron donor acceptor complexes.

Absolute Intensities of NH-Stretching Transitions in Dimethylamine and Pyrrole

Vibrational spectra of vapor-phase dimethylamine (DMA) and pyrrole have been recorded in the 1000 to 13000 cm(-1) region using long path conventional spectroscopy techniques. We have focused on the absolute intensities of the NH-stretching fundamental and overtone transitions; Δν(NH) = 1-4 regions for DMA and the Δν(NH) = 1-3 regions for pyrrole. In the Δν(NH) = 1-3 regions for DMA, evidence of tunneling splitting associated with the NH-wagging mode is observed. For DMA, the fundamental NH-stretching transition intensity is weaker than the first NH-stretching overtone. Also, the fundamental NH-stretching transition in DMA is much weaker than the fundamental transition in pyrrole. We have used an anharmonic oscillator local mode model with ab initio calculated local mode parameters and dipole moment functions at the CCSD(T)/aug-cc-pVTZ level to calculate the NH-stretching intensities and explain this intensity anomaly in DMA.

Explicit correlation and basis set superposition error: The structure and energy of carbon dioxide dimer

We have investigated the slipped parallel and t-shaped structures of carbon dioxide dimer [(CO(2))(2)] using both conventional and explicitly correlated coupled cluster methods, inclusive and exclusive of counterpoise (CP) correction. We have determined the geometry of both structures with conventional coupled cluster singles doubles and perturbative triples theory [CCSD(T)] and explicitly correlated cluster singles doubles and perturbative triples theory [CCSD(T)-F12b] at the complete basis set (CBS) limits using custom optimization routines. Consistent with previous investigations, we find that the slipped parallel structure corresponds to the global minimum and is 1.09 kJ mol(-1) lower in energy. For a given cardinal number, the optimized geometries and interaction energies of (CO(2))(2) obtained with the explicitly correlated CCSD(T)-F12b method are closer to the CBS limit than the corresponding conventional CCSD(T) results. Furthermore, the magnitude of basis set superposition error (BSSE) in the CCSD(T)-F12b optimized geometries and interaction energies is appreciably smaller than the magnitude of BSSE in the conventional CCSD(T) results. We decompose the CCSD(T) and CCSD(T)-F12b interaction energies into the constituent HF or HF CABS, CCSD or CCSD-F12b, and (T) contributions. We find that the complementary auxiliary basis set (CABS) singles correction and the F12b approximation significantly reduce the magnitude of BSSE at the HF and CCSD levels of theory, respectively. For a given cardinal number, we find that non-CP corrected, unscaled triples CCSD(T)-F12b/VXZ-F12 interaction energies are in overall best agreement with the CBS limit.

Calculated Electronic Transitions in Sulfuric Acid and Implications for Its Photodissociation in the Atmosphere

We have calculated electronic transitions for sulfuric acid in the ultraviolet region using a hierarchy of coupled cluster response functions and correlation consistent basis sets. Our calculations indicate that the lowest energy singlet transition occurs at 8.42 eV with an oscillator strength of 0.01. The lowest energy triplet state occurs at 8.24 eV. Thus, the cross section of sulfuric acid in the actinic region is likely to be very small and smaller than the upper limit put on this cross section by previous experimental investigations. We estimate the cross section of sulfuric acid in the atmospherically relevant Lyman-alpha region ( approximately 10.2 eV) to be approximately 6 x 10 (-17) cm (2) molecule (-1), a value approximately 30 times larger than the speculative value used in previous atmospheric simulations. We have calculated the J values for photodissociation of sulfuric acid with absorption of visible, UV, and Lyman-alpha radiation, at altitudes between 30 and 100 km. We find that the dominant photodissociation mechanism of sulfuric acid below 70 km is absorption in the visible region by OH stretching overtone transitions, whereas above 70 km, absorption of Lyman-alpha radiation by high energy Rydberg excited states is the favored mechanism. The low lying electronic transitions of sulfuric acid in the UV region do not contribute significantly to its dissociation at any altitude.

Quantifying cooperative intermolecular interactions for improved carbon dioxide capture materials.

We have optimized the geometry and calculated interaction energies for over 100 different complexes of CO(2) with various combinations of electron accepting (Lewis acid) and electron donating (Lewis base) molecules. We have used the recently developed explicitly correlated coupled cluster singles doubles and perturbative triples [CCSD(T)-F12] methods and the associated VXZ-F12 (where X = D,T,Q) basis sets. We observe only modest changes in the geometric parameters of CO(2) upon complexation, which suggests that the geometry of CO(2) adsorbed in a nanoporous material should be similar to that of CO(2) in gas phase. When CO(2) forms a complex with two Lewis acids via the two electron rich terminal oxygen atoms, the interaction energy is less than twice what would be expected for the same complex involving a single Lewis acid. We consider a series of complexes that exhibit simultaneous CO(2)-Lewis acid and CO(2)-Lewis base intermolecular interactions, with total interaction energies spanning 14.1-105.9 kJ mol(-1). For these cooperative complexes, we find that the total interaction energy is greater than the sum of the interaction energies of the constituent complexes. Furthermore, the intermolecular distances of the cooperative complexes are contracted as compared to the constituent complexes. We suggest that metal-organic-framework or similar nanoporous materials could be designed with adsorption sites specifically tailored for CO(2) to allow cooperative intermolecular interactions, facilitating enhanced CO(2) adsorption.

Intramolecular interactions in 2-aminoethanol and 3-aminopropanol.

Gas-phase vibrational spectra of 2-aminoethanol and 3-aminopropanol were recorded up to the third OH-stretching overtone using Fourier transform infrared spectroscopy, cavity ringdown spectroscopy, and intracavity laser photoacoustic spectroscopy. The experimental investigation was supplemented by local mode calculations, and the intramolecular interactions were investigated using atoms in molecules (AIM) and noncovalent interactions (NCI) theories. All calculations were performed at the CCSD(T)-F12a/VDZ-F12 level of theory. For both compounds the most abundant conformer has a structure that allows for hydrogen bond interaction from the OH group to the N atom of the amino group (OH-N). The spectra show signals from both hydrogen bonded and free OH stretches, implying the presence of several conformers. We observe hydrogen-bond-like interactions in both compounds. The red shift of the bonded OH-stretching frequency and intensity enhancement of the fundamental transition suggest that the hydrogen bond interaction is more pronounced in 3-aminopropanol. AIM analysis supports the presence of a hydrogen bond in 3-aminopropanol but not in 2-aminoethanol, whereas NCI analysis shows hydrogen bonding in both compounds with the stronger interaction found in 3-aminopropanol.