Ab initio calculations, potential representation and vibrational dynamics of He2Br2 van der Waals complex

An intermolecular potential energy surface for He(2)Br(2) complex in the ground state is calculated at the levels of fourth-order (MP4) Moller-Plesset and coupled-cluster [CCSD(T)] approximations, using large-core pseudopotential for Br atoms and the aug-cc-pV5Z basis set for He. The surface is characterized by three minima and the minimum energy pathways through them. The global minimum corresponds to a linear He-Br(2)-He configuration, while the two other ones to "police-nightstick" and tetrahedral structures. The corresponding well depths are -90.39/-89.18, -81.23/-80.78 and -74.40/-74.02 cm(-1), respectively, at MP4/CCSD(T) levels of theory. It is found that results obtained by summing three-body parametrized HeBr(2) interactions and the He-He interaction are in very good accord with the corresponding MP4/CSSD(T) configuration energies of the He(2)Br(2). Variational calculations using a sum of three-body interactions are presented to study the bound states of the vdW He(2)Br(2) complex. The binding energy D(0) and the corresponding vibrationally averaged structure are determined for different isomers of the cluster and their comparison with the available experimental data is discussed.     

Theoretical investigation of the He4Br2 conformers

Full dimensional quantum dynamics calculations of the three lowest isomers of the He(4)Br(2) van der Waals molecule in its ground electronic state are reported. The calculations are performed using the multiconfiguration time-dependent Hartree (MCTDH) method and a realistic potential form that includes the sum of three body ab initio coupled-cluster single double triple [CCSD(T)] He-Br(2) interactions plus the He-He and Br-Br interactions. This potential exhibits several multiple minima, with the three lowest ones lying very close in energy, just within 2 cm(-1). Such small differences are also found in the calculated binding energies of the three most stable conformers, indicating the floppiness of the system and, thus, the need of accurate potential forms and quantum full dynamics methods to treat this kind of complexes. The 12 dimensional results reported in this work present benchmark data and, thus, can serve to evaluate approximate methods aiming to describe higher order rare gas-dihalogen (N > 4) complexes. A comparison with previous studies using different potential forms and approaches to the energetics for the He(4)Br(2) cluster is also presented.     

Ab initio vibrational predissociation dynamics of He-I2(B) complex

Three-dimensional quantum mechanical calculations on the vibrational predissociation dynamics of HeI2 B state complex are performed using a potential energy surface accurately fitted to unrestricted open-shell coupled cluster ab initio data, further enabling extrapolation for large I2 bond lengths. A Lanczos iterative method with an optimized complex absorbing potential is used to determine energies and lifetimes of the vibrationally predissociating He,I2(B,v') complex for v'相似文献   

Vibrational and rotational cooling of NO+ in collisions with He

A quantum mechanical investigation of the vibrational and rotational deactivation of NO(+) in collisions with He atoms in the cold and ultracold regime is presented. Ab initio potential energy calculations are carried out at BCCD(T) level and a new global 3D potential energy surface (PES) is obtained by fitting ab initio points within the reproducing kernel Hilbert space method. As a first test of this PES the bound state energies of the (3)He-NO(+) and (4)He-NO(+) complexes are calculated and compared to previous rigid rotor calculations. The efficiency of the vibrational and the rotational cooling of this molecular ion using a buffer gas of helium is then investigated by performing close coupling scattering calculations for collision energy ranging from 10(-6) to 2000 cm(-1). The calculations are performed for the two isotopes (3)He and (4)He and the results are compared to the available experimental data.     

Ab initio potential energy surface, infrared spectra, and dynamics of the ion-molecule complexes between Br- and H2, D2, and HD

A three-dimensional potential energy surface (PES) for the Br(-)-H(2) complex is computed using the ab initio CCSD(T) method and an extended basis set. The PES has two equivalent minima at the linear geometries (equilibrium interfragment distance R(e)=3.34 A and interaction energy D(e)=670 cm(-1)) separated by the barrier at the T-shaped configuration (interfragment distance R(e)=4.03 A and barrier height of 610 cm(-1)). Ab initio points are fitted to a flexible analytical form and used in the variational rovibrational energy level calculations. Simulated infrared spectra of the Br(-)-H(2) and Br(-)-D(2) complexes in the monomer stretching excitation region are in good agreement with the measured ones. Nonstatistical intensity ratios of the complexes of para- and ortho-monomers are qualitatively explained by monomer ligand exchange reactions. Predissociation of the complexes containing vibrationally excited monomers is analyzed and shown to proceed through the near-resonant vibration-to-rotation energy transfer. For complexes involving Br(-) and the HD monomer, two energetically low-lying states are predicted, corresponding predominately to the Br(-)-DH and Br(-)-HD isomeric forms. The results demonstrate the close similarity of the bromide containing complexes to their analogs containing the chloride anion.     

Cartesian coupled coherent states simulations: Ne(n)Br2 dissociation as a test case

In this article, we describe coupled coherent states (CCS) simulations of vibrational predissociation of weakly bounded complexes. The CCS method is implemented in the Cartesian frame in a manner that is similar to classical molecular dynamics. The calculated lifetimes of the vibrationally excited Ne-Br(2)(ν) complexes agree with experiment and previous calculations. Although the CCS method is, in principle, a fully quantum approach, in practice it typically becomes a semiclassical technique at long times. This is especially true following dissociation events. Consequently, it is very difficult to converge the quantum calculations of the final Br(2) vibrational distributions after predissociation and of the autocorrelation functions. However, the main advantage of the method is that it can be applied with relative ease to determine the lifetimes of larger complexes and, in order to demonstrate this, preliminary results for tetra- and penta-atomic clusters are reported.     

Ab initio characterization of the Ne-I2 van der Waals complex: intermolecular potentials and vibrational bound states

A theoretical study of the potential energy surface and bound states is performed for the ground state of the NeI(2) van der Waals (vdW) complex. The three-dimensional interaction energies are obtained from ab initio coupled-cluster, coupled-cluster single double (triple)/complete basis set, calculations using large basis sets, of quadruple- through quintuple-zeta quality, in conjunction with relativistic effective core potentials for the heavy iodine atoms. For the analytical representation of the surface two different schemes, based on fitting and interpolation surface generation techniques, are employed. The surface shows a double-minimum topology for linear and T-shaped configurations. Full variational quantum mechanical calculations are carried out using the model surfaces, and the vibrationally averaged structures and energetics for the NeI(2) isomers are determined. The accuracy of the potential energy surfaces is validated by a comparison between the present results and the corresponding experimental data available. In lieu of more experimental measurements, we also report our results/predictions on higher bound vibrational vdW levels, and the influence of the employed surface on them is discussed.     

Efficiency considerations in the construction of interpolated potential energy surfaces for the calculation of quantum observables by diffusion Monte Carlo

A modified Shepard interpolation scheme is used to construct global potential energy surfaces (PES) in order to calculate quantum observables--vibrationally averaged internal coordinates, fully anharmonic zero-point energies and nuclear radial distribution functions--for a prototypical loosely bound molecular system, the water dimer. The efficiency of PES construction is examined with respect to (a) the method used to sample configurational space, (b) the method used to choose which points to add to the PES data set, and (c) the use of either a one- or two-part weight function. The most efficient method for constructing the PES is found to require a quantum sampling regime, a combination of both h-weight and rms methods for choosing data points and use of the two-part weight function in the interpolation. Using this regime, the quantum diffusion Monte Carlo zero-point energy converges to the exact result within addition of 50 data points. The vibrationally averaged O-O distance and O-O radial distribution function, however, converge more slowly and require addition of over 500 data points. The methods presented here are expected to be applicable to both other loosely bound complexes as well as tightly bound molecular species. When combined with high quality ab initio calculations, these methods should be able to accurately characterize the PES of such species.     

Angle-energy distributions of Penning ions in crossed molecular beams. IV. He*(21 S,2 3S)+H2-->He+H2+ +e-

Relative doubly differential cross sections for the Penning ionization of H(2) by spin-state-selected metastable He (1s2s) are reported at center-of-mass collision energies E of 3.1 and 4.2 kcal/mol in a crossed supersonic beam experiment employing a rotatable mass spectrometer detector. The measurements are sufficiently dense in velocity space as to avoid having to functionalize the differential cross sections in order to transform the intensities into the c.m. The H(2) (+) product is scattered sharply forward, c.m. Deltatheta<10 degrees half-width at half-maximum, with respect to the incident direction of H(2) at both energies for both spin states. On the average the products have lost energy upon recoil, mean recoil energy E(')相似文献   

The COHe system revisited: V—T collisional relaxation from AB initio potentials

A new potential energy surface (PES), recently computed within an ab initio scheme for the CO + He system, is here used to calculate vibrationally inelastic cross sections via a coupled quantum treatment of the dynamics. Results are obtained with different decoupling approximations (BS and IOS) and compared with those yielded by a previously available PES [3] that showed too strong an attractive behaviour. The effect of anisotropic components in the potential is also analysed and discussed for both PES and appears to play a noticeable role only at low collision energies. Relaxation rates are also computed and used to obtain relaxation times within the Landau—Teller (LT) scheme. Comparison with experiments is very satisfactory at high temperatures and medium temperatures (3000 ? T ? 200 K) while low T behaviour is overestimated by the present calculations.     

A theoretical study of He2ICl van der Waals cluster

The structure, energetics, and dynamics of He2ICl complex in its ground state are studied by means of ab initio electronic structure and quantum-mechanical calculations. Interaction energies for selected He2ICl configurations are calculated at the coupled-cluster [CCSD(T)] level of theory using a large-core pseudopotential for the I atom and the aug-cc-pVTZ and aug-cc-pV5Z basis sets for the Cl and He atoms, respectively. The surface is characterized around its lower five minima and the minimum energy pathways through them. The global minimum of the potential corresponds to a "police-nightstick (1)" configuration, the second one to a linear, the next one to tetrahedral configuration, and the following two to "bifork" and "police-nightstick (2)" structures, with well depths of -99.12, -97.42, -88.32, -85.84, and -78.54 cm(-1), respectively. An analytical form based on the sum of the three-body parametrized HeICl interactions plus the He-He interaction is found to represent very well the tetra-atomic CSSD(T) results. The present potential expression is employed to perform variational five-dimensional quantum-mechanical calculations to study the vibrational bound states of the van der Waals He2ICl complex. Results for total angular momentum J = 0 provide the binding energy D0 and the corresponding vibrationally averaged structure for different isomers of the cluster. Comparison of these results with recent experimental observations further justifies the potential used in this work.     

Collision-energy-resolved Penning ionization electron spectroscopy of bromomethanes (CH3Br, CH2Br2, and CHBr3) by collision with He*(2(3)S) metastable atoms

Ionization of bromomethanes (CH3Br, CH2Br2, and CHBr3) upon collision with metastable He*(2(3)S) atoms has been studied by means of collision-energy-resolved Penning ionization electron spectroscopy. Lone-pair (nBr) orbitals of Br4p characters have larger ionization cross sections than sigma(C-Br) orbitals. The collision-energy dependence of the partial ionization cross sections shows that the interaction potential between the molecule and the He*(2(3)S) atom is highly anisotropic around CH3Br or CH2Br2, while isotropic attractive interactions are found for CHBr3. Bands observed at electron energies of approximately 2 eV in the He*(2(3)S) Penning ionization electron spectra (PIES) of CH2Br2 and CHBr3 have no counterpart in ultraviolet (He I) photoionization spectra and theoretical (third-order algebraic diagrammatic construction) one-electron and shake-up ionization spectra. Energy analysis of the processes involved demonstrates that these bands and further bands overlapping with sigma(C-Br) or piCH2 levels are related to autoionization of dissociating (He+ - Br-) pairs. Similarly, a band at an electron energy of approximately 1 eV in the He*(2(3)S) PIES spectra of CH3Br has been ascribed to autoionizing Br** atoms released by dissociation of (unidentified) excited states of the target molecule. A further autoionization (S) band can be discerned at approximately 1 eV below the lone-pair nBr bands in the He*(2(3)S) PIES spectrum of CHBr3. This band has been ascribed to the decay of autoionizing Rydberg states of the target molecule (M**) into vibrationally excited states of the molecular ion. It was found that for this transition, the interaction potential that prevails in the entrance channel is merely attractive.     

On the extent of intramolecular hydrogen bonding in gas-phase and hydrated 1,2-ethanediol

We investigate the quantum dynamical nature of hydrogen bonding in 1,2-ethanediol and monohydrated 1,2-ethanediol using different levels of ab initio theory. Global full-dimensional potential energy surfaces were constructed from PW91/cc-pVDZ, B3LYP/cc-pVDZ, and MP2/cc-pVDZ ab initio data for gas-phase and monohydrated 1,2-ethanediol, using a modified Shepard interpolation scheme. Zero-point energies and nuclear vibrational wave functions were calculated on these surfaces using the quantum diffusion Monte Carlo algorithm. The nature of intra- and intermolecular hydrogen bonding in these molecules was investigated by considering a ground-state nuclear vibrational wavefunction with reduced complete nuclear permutation and inversion (CNPI) symmetry. Separate wavefunction histograms were determined from the ground-state nuclear vibrational wavefunction by projection into bondlength coordinates. The O-H and O-O wavefunction histograms and vibrationally averaged distances were then used to probe the extent of intra- and intermolecular hydrogen bonding. From these data, we conclude that gas-phase ethanediol may possess a weak hydrogen bond, with a relatively short O-O distance but no detectable proton delocalization. Monohydrated ethanediol was found to exhibit no intramolecular hydrogen bonding but instead possessed two intermolecular hydrogen bonds, indicated by both shortening of the O-O distance and significant proton delocalization. The degree of proton delocalization and shortening of the vibrationally averaged O-O distance was found to be dependent on the ab initio method used to generate the potential energy surface (PES) data set.     

Reduced dimensionality quantum scattering calculations on the F + CH(4) --> FH + CH(3) reaction

The existence of recently observed scattering resonances in the hydrogen abstraction reaction F + CH4 --> FH + CH3 was investigated using the reduced dimensionality rotating line umbrella (RLU) quantum scattering model and employing an analytical potential energy surface, PES-2006, recently developed by our group. The calculations were performed in hyperspherical coordinates. The wells found in the hyperspherical adiabats, the oscillatory pattern in the cumulative and state-to-state reaction probabilities, the forward/backward predominance in the differential cross section at a collision energy of 1.8 kcal mol(-1), and the dramatic change of the scattering angle with energy are related to scattering resonances, and they are assigned to a quasi-bound complex on the vibrationally adiabatic ground-state potential.     

Stabilization and rovibronic spectra of the T-shaped and linear ground-state conformers of a weakly bound rare-gas-homonuclear dihalogen complex: He...Br2

Laser-induced fluorescence spectra of Br(2) entrained in a He supersonic expansion have been recorded in the Br(2) B-X, 8-0, 12-0, and 21-0 spectral regions at varying downstream distances, and thus different temperature regimes. Features associated with transitions of the T-shaped and linear He...Br(2)(X,nu(") = 0) complexes are identified. The changes in the relative intensities of the T-shaped and linear features with cooling in the expansion indicate that the linear conformer is energetically more stable than the T-shaped conformer. A He + Br(2)(X,nu(") = 0) ab initio potential-energy surface, computed at the coupled cluster level of theory with a large, flexible basis set, is used to calculate the binding energies of the two conformers, 15.8 and 16.5 cm(-1) for the T-shaped and linear complexes, respectively. This potential and an excited-state potential [M. P. de Lara-Castells, A. A. Buchachenko, G. Delgado-Barrio, and P. Villareal, J. Chem. Phys. 120, 2182 (2004)] are used to calculate the excitation spectra of He...(79)Br(2)(X,nu(") = 0) in the Br(2) B-X, 12-0 region. The calculated spectra are used to make spectral assignments and to determine the energies of the excited-state intermolecular vibrational levels accessed in the observed transitions. Temperature-dependent laser-induced fluorescence spectra and a simple thermodynamic model [D. S. Boucher, J. P. Darr, M. D. Bradke, R. A. Loomis, and A. B. McCoy, Phys. Chem. Chem. Phys. 6, 5275 (2004)] are used to estimate that the linear conformer is 0.4(2) cm(-1) more strongly bound than the T-shaped conformer. Two-laser action spectroscopy experiments reveal that the binding energy of the linear He...(79)Br(2)(X,nu(") = 0) conformer is 17.0(8) cm(-1), and that of the T-shaped He...(79)Br(2)(X,nu(") = 0) conformer is then 16.6(8) cm(-1), in good agreement with the calculated values.     

CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule

The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed.     

Theoretical investigation of the Na+ + H2 system. III. Trajectory calculations of vibrationally inelastic collision processes

On the basis of an analytical potential energy surface for the electronic ground state of the Na⁺ + H₂ system reported recently, extensive trajectory calculations have been performed to study the collision dynamics of vibrationally inelastic processes at total energies up to ～3 eV. Special attention is given to the relative efficiacy of translational and rotational energy, respectively, in promoting vibrational energy transfer. Vibrational transitions are found to be substantially enhanced by initial molecular rotation. Furthermore, the applicability of simple models is discussed.     

Bound state potential energy surface construction: ab initio zero-point energies and vibrationally averaged rotational constants

Collins' method of interpolating a potential energy surface (PES) from quantum chemical calculations for reactive systems (Jordan, M. J. T.; Thompson, K. C.; Collins, M. A. J. Chem. Phys. 1995, 102, 5647. Thompson, K. C.; Jordan, M. J. T.; Collins, M. A. J. Chem. Phys. 1998, 108, 8302. Bettens, R. P. A.; Collins, M. A. J. Chem. Phys. 1999, 111, 816) has been applied to a bound state problem. The interpolation method has been combined for the first time with quantum diffusion Monte Carlo calculations to obtain an accurate ground state zero-point energy, the vibrationally average rotational constants, and the vibrationally averaged internal coordinates. In particular, the system studied was fluoromethane using a composite method approximating the QCISD(T)/6-311++G(2df,2p) level of theory. The approach adopted in this work (a) is fully automated, (b) is fully ab initio, (c) includes all nine nuclear degrees of freedom, (d) requires no assumption of the functional form of the PES, (e) possesses the full symmetry of the system, (f) does not involve fitting any parameters of any kind, and (g) is generally applicable to any system amenable to quantum chemical calculations and Collins' interpolation method. The calculated zero-point energy agrees to within 0.2% of its current best estimate. A0 and B0 are within 0.9 and 0.3%, respectively, of experiment.