H2, HD, and D2 inside C60: coupled translation-rotation eigenstates of the endohedral molecules from quantum five-dimensional calculations

We have performed rigorous quantum five-dimensional (5D) calculations of the translation-rotation (T-R) energy levels and wave functions of H(2), HD, and D(2) inside C(60). This work is an extension of our earlier investigation of the quantum T-R dynamics of H(2)@C(60) [M. Xu et al., J. Chem. Phys. 128, 011101 (2008)] and uses the same computational methodology. Two 5D intermolecular potential energy surfaces (PESs) were employed, differing considerably in their well depths and the degree of confinement of the hydrogen molecule. Our calculations revealed pronounced sensitivity of the endohedral T-R dynamics to the differences in the interaction potentials, and to the large variations in the masses and the rotational constants of H(2), HD, and D(2). The T-R levels vary significantly in their energies and ordering on the two PESs, as well as from one isotopomer to another. Nevertheless, they all display the same distinctive patterns of degeneracies, which can be qualitatively understood and assigned in terms the model which combines the isotropic three-dimensional harmonic oscillator, the rigid rotor, and the coupling between the orbital and the rotational angular momenta of H(2)/HD/D(2). The quantum number j associated with the rotation of H(2), HD, and D(2) was found to be a good quantum number for H(2) and D(2) on both PESs, while most of the T-R levels of HD exhibit strong mixing of two or more rotational basis functions with different j values.     

Liquid-phase force field and dipole moment derivatives with respect to internal coordinates of benzene

This paper presents an analysis of the infrared vibrational intensities found for C(6)H(6), C(6)D(6) and C(6)H(5)D in the liquid phase, motivated in part by the quite marked intensity differences between the fundamentals of C(6)H(6) and C(6)D(6) in the liquid, and between corresponding vibrations in the liquid and gas phases. The analysis is carried out under the harmonic approximation and results from a determination of the force field for liquid C(6)H(6), C(6)D(6) and C(6)H(5)D. The force constants for the liquid-phase are presented and compared to those in the literature for the gas-phase. Previously reported experimental intensities are used along with the eigenvectors of the force field analysis to determine the dipole moment derivatives with respect to symmetry and internal coordinates. The dipole moment derivatives with respect to internal coordinates obtained are partial differentialmicro/ partial differentials=0.38+/-0.02DebyeA(-1), partial differentialmicro/ partial differentialt=0.24+/-0.01, partial differentialmicro/ partial differentialbeta=0.26+/-0.01, and partial differentialmicro/ partial differentialgamma=0.64+/-0.03DebyeA(-1). There is very little difference between the dipole moment derivatives with respect to internal coordinates obtained from non-linear least squares fitting of the two D(6h) isotopomers and those obtained from non-linear least squares fitting of the three isotopomers. The results show that there is significant intensity sharing in the CH stretch region of C(6)H(5)D between the fundamental and combination bands.     

Quantum statistical and wave packet studies of insertion reactions of S(1D) with H2, HD, and D2

A thorough theoretical investigation of the reactions between S(1D) and various hydrogen isotopomers (H2, D2, and HD) has been carried out using a recent ab initio potential energy surface. State-resolved integral and differential cross sections, thermal rate constants, and their dependence on energy or temperature were obtained from quantum mechanical capture probabilities within a statistical model. For comparison, the J=0 reaction probabilities were also computed using an exact wave packet method. The statistical results are in excellent agreement with available exact differential and integral cross sections. The comparison with experimental results shows that the agreement is reasonably good in general, but some significant differences exist, particularly for the SD/SH branching ratio in the S(1D)+HD reaction.     

Rotational spectrum of the dimethyl ether-acetylene complex: evidence for an effective C2v geometry

The rotational spectra for five isotopomers of the 1:1 weakly bound complex formed between dimethyl ether (DME) and acetylene (HCCH) have been measured by Fourier transform microwave spectroscopy. The experimental rotational constants, planar moments, and dipole moment components are consistent with a floppy complex possessing an effective C2v structure in which the hydrogen atom of acetylene is hydrogen bonded to the oxygen atom of dimethyl ether with an intermolecular H...O separation of 2.08(3) A. Experimental rotational constants for the normal isotopic species are A = 10382.5(17) MHz, B = 1535.7187(18) MHz, and C = 1328.3990(17) MHz and the dipole moment components are mua= mutotal = 1.91(10) D. Ab initio calculations at the MP2/6-311++G(2d,2p) level indicate that the energy barrier for motion of the HCCH subunit between the lone pairs of the DME, via a C2v intermediate structure, is very low (approximately 0.29 kJ mol(-1)). Inclusion of basis set superposition error and zero point energy corrections to the energies of four stationary points located on the potential energy surface shows that the relative stabilities are particularly sensitive to these corrections. The ab initio optimizations give rotational constants for the C2v structure of A = 10066 MHz, B = 1496 MHz, and C = 1324 MHz, and a dipole moment of mua= mu(total) = 2.12 D, in reasonable agreement with the experimentally determined values. The structural parameters and energetics of the DME-HCCH complex will be discussed and compared to similar complexes such as H2O-HCCH.     

Reactions of C(1D) with H2 and its deuterated isotopomers, a wave packet study

Using a Chebyshev wave packet method, total and state-resolved reaction probabilities (J=0) were calculated for the reactions of C(1D) with various hydrogen isotopomers (H2, D2, and HD, nu i=0, j i=0) on a recent ab initio potential energy surface. For all the isotopic variants, it was found that the initial state specified reaction probabilities have no energy threshold and are strongly oscillatory, indicative of the involvement of long-lived resonances in this barrierless reaction. The J=0 product vibrational and rotational distributions for all three isotopic reactions, and the CH/CD branching ratio for the C+HD reaction, show strong dependence on the collision energy, further underscoring the important role played by the resonances. The generally decaying vibrational distributions and highly excited rotational distributions, which corroborate an insertion mechanism, and the dominance of the CD+H channel in the C+HD reaction are consistent with existing experimental observations. Initial state specified integral cross sections and rate constants were estimated using a capture model. The estimated rate constants were found to be close and in the order kHD>kH2>kD2. Finally, a method to calculate branching ratio in the C+HD reaction is proposed.     

Rotational spectrum and structure of the carbonyl sulfide-trifluoromethane weakly bound dimer

Pure rotational spectra of five isotopomers of the 1:1 weakly bound complex formed between carbonyl sulfide and trifluoromethane (TFM) have been measured using Fourier transform microwave spectroscopy. The experimental rotational constants and dipole moment components are consistent with a structure of C(s) symmetry in which the dipole moment vectors of OCS and HCF(3) are aligned antiparallel and at an angle of about 40 degrees and with a center of mass separation of 3.965(26) A. The derived H...O distance is 2.90(5) A, which is up to 0.6 A longer than is seen in other similar TFM complexes exhibiting C-H...O interactions. Ab initio calculations at the MP2/6-311++G(2d,2p) level give a structure with rotational constants that are in reasonable agreement with those of the normal isotopic species.     

Radiative transition probabilities, lifetimes and dipole moments for the vibrational levels of the X1Sigma+ ground state of 39K85Rb

Using a potential energy curve (based primarily on the RKR potential of Amiot and Verges [J. Chem. Phys. 112, 7068 (2000)]) and a dipole moment function (based primarily on ab initio calculations of Park et al. [Chem. Phys. 257, 135 (2000)]), we have calculated radiative transition probabilities (Einstein A coefficients), radiative lifetimes, and dipole moment expectation values involving all vibrational levels (for several rotational quantum numbers) of the X1Sigma+ ground state of 39K85Rb. We observe that the radiative lifetimes of vibrationally excited levels, in particular, are approximately 10(3)-10(6) seconds, far too long to be significant in most ultracold experiments involving 39K85Rb or its isotopomers. Comparison with other molecules (LiH and HF) suggests that simple scaling (A approximately mu2nu3 approximately tau(-1)) will predict similarly long lifetimes for many other heteronuclear molecules, e.g., RbCs.     

Quantum vibrational analysis and infrared spectra of microhydrated sodium ions using an ab initio potential

We present a full-dimensional potential energy surface and a dipole moment surface (DMS) for hydrated sodium ion. These surfaces are based on an n-body expansion for both the potential energy and the dipole moment, truncated at the two-body level for the H(2)O-Na(+) interaction and also for the DMS. The water-water interaction is truncated at the three-body level. The new full-dimensional two-body H(2)O-Na(+) potential is a fit to roughly 20,000 coupled-cluster single double (triple)/aug-cc-pVTZ energies. Properties of this two-body potential and the potential describing (H(2)O)(n)Na(+) clusters, with n up to 4 are given. We then report anharmonic, coupled vibrational calculations with the "local-monomer model" to obtain infrared spectra and also 0 K radial distribution functions for these clusters. Some comparisons are made with the recent infrared predissociation spectroscopy experiments of Miller and Lisy [J. Am. Chem. Soc. 130, 15381 (2008).].     

H/D isotope effect on porphine and porphycene molecules with multicomponent hybrid density functional theory

To analyze the H/D isotope effect on porphine and porphycene molecules including the protonic/deuteronic quantum nature and electron correlation efficiently, the authors have developed the new scheme of the multicomponent hybrid density functional theory [MC_(HF+DFT)]. The optimized geometries of porphine, porphycene, and these deuterated isotopomers by our MC_(HF+DFT) method are in good agreement with the experimental "high-symmetric" structures, contrary to the "low-symmetric" geometries optimized by pure multicomponent Hartree-Fock method. The optimized geometries for HD-porphine and HD-porphycene molecules, in which an inner hydrogen is replaced to a deuterium, are found to be low symmetric. Such drastic geometrical change induces the electronic polarization, and gives rise to the slight dipole moment values in these HD species. Their results clearly indicate that the difference of the nuclear quantum nature between inner proton and inner deuteron directly influences the molecular geometry and electronic structure.     

Microwave and ab initio studies of rare gas-methane van der Waals complexes

Rotational spectra of the weakly bound Kr-methane van der Waals complex were recorded using a pulsed molecular beam Fourier transform microwave spectrometer in the range from 3.5 to 18 GHz. Spectra of 25 isotopomers of Kr-methane were assigned and analyzed. For isotopomers containing CH4, 13CH4, and CD4, two sets of transitions with K = 0 and one with K = 1 were recorded, correlating to the j = 0, 1, and 2 rotational levels of free methane, respectively (j is the rotational angular momentum quantum number of the methane monomer). For isotopomers containing CH3D and CHD3, two K = 0 components were recorded, correlating to the j(k) = 0(0) and 1(1) rotational levels of free methane (k corresponds to the projection of j onto the C3 axis of CH3D and CHD3). The obtained spectroscopic results were used to derive van der Waals bond distance R, van der Waals stretching frequency nu(s), and the corresponding stretching force constant k(s). Nuclear spin statistical weights of individual states were obtained from molecular symmetry group analyses and were compared with the observed relative transition intensities. The tentatively assigned j = 2 transitions were more intense than predicted from symmetry considerations. This is attributed to a relatively large effective dipole moment of this state, supported by ab initio dipole moment calculations. Ab initio potential energy calculations of Kr-CH4 and Ar-CH4 were done at the coupled cluster level of theory, with single and double excitations and perturbative inclusion of triple excitations, using the aug-cc-pVTZ basis set supplemented with bond functions. The theoretical results show that the angular dynamics of the dimer does not change significantly when the binding partner of methane changes from Ar to Kr. The dipole moment of Ar-CH4 was calculated at various configurations, providing a qualitative explanation for the unsuccessful spectral searches for rotational transitions of Ar-CH4.     

An ab initio study of the structure, torsional potential energy function, and electric properties of disilane, ethane, and their deuterated isotopomers

Highly accurate ab initio computations of the molecular structure and properties, torsional potential energy function, and harmonic force field of disilane and ethane have been carried out. Equilibrium parameters as well as vibrational corrections have been evaluated. In addition, for these systems a vibrational averaging procedure has been employed for calculating the dipole moment of molecules which have no permanent dipole moment, i.e., SiH(3)SiD(3) and CH(3)CD(3). The molecular and spectroscopic properties calculated for ethane and its isotopomers provide a calibration against known experimental data, allowing us to estimate the reliability of our computed results for disilane for which there is much less experimental data. The goal of the present study is to predict the molecular parameters, with estimated uncertainties, that determine the microwave spectrum of SiH(3)SiD(3).     

Ab initio potential energy and dipole moment surfaces of (H2O)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with M?ller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).     

Supersonically cooled hydronium ions in a slit-jet discharge: high-resolution infrared spectroscopy and tunneling dynamics of HD2O+

Jet-cooled high-resolution infrared spectra of partially deuterated hydronium ion (HD2O+) in the O-H stretch region (nu3 band) are obtained for the first time, exploiting the high ion densities, long absorption path lengths, and concentration modulation capabilities of the slit-jet discharge spectrometer. Least-squares analysis with a Watson asymmetric top Hamiltonian yields rovibrational constants and provides high level tests of ab initio molecular structure predictions. Transitions out of both the lower (nu3(+)<--0(+)) and the upper (nu3(-)<--0(-)) tunneling levels, as well as transitions across the tunneling gap (nu3(-)<--0(+)) are observed. The nu3(-)<--0(+) transitions in HD2O+ acquire oscillator strength by loss of D(3h) symmetry, and permit both ground-state-[27.0318(72) cm(-1)] and excited-state-[17.7612(54) cm(-1)]-tunneling splittings to be determined to spectroscopic precision from a single rovibrational band. The splittings and band origins calculated with recent high level ab initio six-dimensional potential surface predictions for H3O+ and isotopomers [X. C. Huang, S. Carter, and J. M. Bowman, J. Chem. Phys. 118, 5431 (2003); T. Rajamaki, A. Miani, and L. Halonen, J. Chem. Phys. 118, 10929 (2003)] are in very good agreement with the current experimental results.     

Vibrational ground state properties of H(5)(+) and its isotopomers from diffusion Monte Carlo calculations

Diffusion Monte Carlo computations, with and without importance sampling, of the zero-point properties of H(5)(+) and its isotopomers using a recent high accuracy global potential energy surface are presented. The global minimum of the potential possesses C(2v) symmetry, but the calculations predict a D(2d) geometry for zero-point averaged structure of H(5)(+) with one H atom "in the middle" between two HH diatoms. The predicted zero-point geometries of the deuterated forms have H in the middle preferred over D in the middle and for a nonsymmetric arrangement of D atoms the preferred arrangement is one which maximizes the number of D as the triatomic ion. We speculate on the consequences of these preferences in scattering of H(2)+H(3)(+) and isotopomers at low energies, such as those in the interstellar medium.     

An interpretation of the enhancement of the water dipole moment due to the presence of other water molecules

The dipole moment of the gas phase water monomer is 1.85 D. When solvated in bulk water, the dipole moment of an individual water molecule is observed to be enhanced to the much larger value of 2.9 +/- 0.6 D. To understand the origin of this dipole moment enhancement, the effective fragment potential (EFP) method is used to solvate an ab initio water molecule to predict the dipole moments for various cluster sizes. The dipole moment as a function of cluster size, nH 2O, is investigated [for n = 6-20 (even n), 26, 32, 41, and 50]. Localized charge distributions are used in conjunction with localized molecular orbitals to interpret the dipole moment enhancement. These calculations suggest that the enhancement of the dipole moment originates from the decrease of the angle between the dipole vectors of the lone pairs on oxygen as the number of hydrogen bonds to that oxygen increases. Thus, the decreased angle, and the consequent increase in water dipole moment, is most likely to occur in environments with a larger number of hydrogen bonds, such as the center of a cluster of water molecules.     

Quantum dynamics of H2, D2, and HD in the small dodecahedral cage of clathrate hydrate: evaluating H2-water nanocage interaction potentials by comparison of theory with inelastic neutron scattering experiments

We have performed rigorous quantum five-dimensional (5D) calculations and analysis of the translation-rotation (T-R) energy levels of one H(2), D(2), and HD molecule inside the small dodecahedral (H(2)O)(20) cage of the structure II clathrate hydrate, which was treated as rigid. The H(2)- cage intermolecular potential energy surface (PES) used previously in the molecular dynamics simulations of the hydrogen hydrates [Alavi et al., J. Chem. Phys. 123, 024507 (2005)] was employed. This PES, denoted here as SPC/E, combines an effective, empirical water-water pair potential [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and electrostatic interactions between the partial charges placed on H(2)O and H(2). The 5D T-R eigenstates of HD were calculated also on another 5D H(2)-cage PES denoted PA-D, used by us earlier to investigate the quantum T-R dynamics of H(2) and D(2) in the small cage [Xu et al., J. Phys. Chem. B 110, 24806 (2006)]. In the PA-D PES, the hydrogen-water pair potential is described by the ab initio 5D PES of the isolated H(2)-H(2)O dimer. The quality of the SPC/E and the PA-D H(2)-cage PESs was tested by direct comparison of the T-R excitation energies calculated on them to the results of two recent inelastic neutron scattering (INS) studies of H(2) and HD inside the small clathrate cage. The translational fundamental and overtone excitations, as well as the triplet splittings of the j=0-->j=1 rotational transitions, of H(2) and HD in the small cage calculated on the SPC/E PES agree very well with the INS results and represent a significant improvement over the results computed on the PA-D PES. Our calculations on the SPC/E PES also make predictions about several spectroscopic observables for the encapsulated H(2), D(2), and HD, which have not been measured yet.     

Direct visual evidence for quinoidal charge delocalization in poly-p-phenylene cation radical

Recently, X-ray crystallographic evidence of quinoidal charge delocalization in poly-p-phenylene cation radicals was reported [Banerjee, M. et al., J. Am. Chem. Soc. 2007, 129, 8070]. In this paper, direct visual evidence for quinoidal charge delocalization in quaterphylene (QP) is shown with three-dimensional (3D) charge difference densities. It was revealed that the extra positive charge mainly localized on the two center units at the ground state, while the extra positive charge will delocalize to the two outer units upon electronic state transitions by photoexcitation. The 2D plots together with the corresponding charge difference densities were interpreted as large distance-charge oscillations, implying that in the positive species upon excitation a nearly free oscillating motion of a hole occurs. For the QP cation radical, the transition dipole moment of S1 represents mesoscopic dipole antennae.     

Parallel tempering simulations of the 13-center Lennard-Jones dipole-dipole cluster (muD=0-->0.5 a.u.)

We investigate the thermodynamic behavior of the thirteen center uniform Lennard-Jones dipole-dipole cluster [(LJDD)(13)] for a wide range of dipole moment strengths. We find a relatively wide range of potential parameters where solid-solid coexistence manifests itself. Using structural characterization methods we determine the shape of the few isomers that contribute to the solid-solid coexistence region. The thermal distributions of the size of the net dipole moment are broad even at the coldest temperatures of the simulation where the (LJDD)(13) cluster is solid.     

Ab initio potential and dipole moment surfaces for water. II. Local-monomer calculations of the infrared spectra of water clusters

We employ recent flexible ab initio potential energy and dipole surfaces [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] to the calculation of IR spectra of the intramolecular modes of water clusters. We use a quantum approach that begins with a partitioned normal-mode analysis of perturbed monomers, and then obtains solutions of the corresponding Schro?dinger equations for the fully coupled intramolecular modes of each perturbed monomer. For water clusters, these modes are the two stretches and the bend. This approach is tested against benchmark calculations for the water dimer and trimer and then applied to the water clusters (H(2)O)(n) for n = 6-10 and n = 20. Comparisons of the spectra are made with previous ab initio harmonic and empirical potential calculations and available experiments.     

The ground-state rotational spectrum and molecular geometry of ethynylstannane

The ground-state rotational spectra of 24 isotopomers of ethynylstannane have been observed by pulsed-jet, Fourier-transform microwave spectroscopy. The spectroscopic constants, B(0,)D(J) and D(JK) are reported for symmetric-top isotopomers H(3)(n)Sn(12)C(12)CH, where n = 116, 117, 118, 119, 120, 122 and 124, D(3)(n)Sn(12)C(12)CH, where n = 116, 118, 120, 122 and 124, H(3)(n)Sn(13)C(12) CH and H(3)(n)Sn(12)C(13)CH , where n = 116,118 and 120, and H(3)(n)Sn(12)C(12)CD, where n = 116, 118 and 120. In addition, the values of A(0), B(0), C(0), Delta(J) and Delta(JK) were obtained for the three asymmetric-top isotopomers DH(2)(n)Sn(12)C(12)CH, where n = 116, 118 and 120. Hyperfine structure was resolved and assigned in the transitions of the isotopomers H(3)(n)SnCCD, where n = 116, 118 and 120, and in the isotopomers H(3)(117)SnCCH and H(3)(119)SnCCH. In the former group, the hyperfine structure arises from D nuclear quadrupole coupling while in the latter group its origin lies in the spin-rotation coupling of the I = 1/2 Sn nuclear spin to the rotational motion. For these isotopomers, D nuclear quadrupole and spin-rotation coupling constants are determined where appropriate. The rotational constants obtained for the 24 isotopomers of H(3)SnCCH were used to obtain the following types of molecular geometry for ethynylstannane: r(0), r(s), and r(m).