Potential Energy Function for State X1∑+ of LaH

The present work has derived the potential energy function for the ground state X1∑+of LaH in Murrell-Sorbie function form. In the first place, the electronic state and it′ s rea sonable dissodation limits are correctly determined based on Atomic and Molecular Reaction Statics (AMRS), and then, using the relativistic compact effective potential (RCEP) for La, the equilibrium geometry and dissociation energy for LaH have been calculated by the QCISD method. The calculated results for Re'De' Be' αe' ωe and ωexe are 2.125A, 2.623eV,3.7333cm-1, 0.0723 cm-1, 1461.73cm-1 and 21.383cm-1 respectively, which are in good agreement with experimental or calculated values in references.     

二氧化钚分子的多体展式势能函数

从导出基态PuO2分子的电子状态X5Σ g正确地判断其离解极限出发,采用MP2方法,应用相对论有效原子实模型(RECP)优化出PuO2(X5Σ g)分子稳定构型为线性OPuO(D∞h),其平衡核间距Re=0.18004nm.同时也计算出振动频率,并优化出存在亚稳态的Pu-O-O(C∞v)构型.使用多体项展式理论方法,导出了基态PuO2分子的分析势能函数.该势能表面准确地再现了O-Pu-O(D∞h)平衡结构和亚稳态的Pu-O-O(C∞v)构型.然后根据势能函数等值图讨论了O(3Pg) PuO反应的势能面静态特征.     

PdH2、YH2分子的结构与势能函数

用密度泛函理论的B3LYP方法,对钯和钇原子采用SDD收缩价基函数,氢原子采用6-311＋＋G**全电子基函数,对PdH2和YH2体系的结构进行优化计算,得到PdH2分子最稳态为C2v构型,电子组态为1A1,平衡核间距RPdH=0.1692 nm,键角∠HPdH=29.4°,离解能De=5.5212 eV,基态简正振动频率:ν1(b2)=1470.1 cm－1、ν2(a1)=1007.9 cm－1、ν3(a1)=2907.0 cm－1.YH2分子最稳态也为C2v构型,电子组态2A1,RYH=0.1962 nm,∠HYH=114.3°,De=5.6691 eV,基态简正振动频率:ν1(b2)=1457.9 cm－1、ν2(a1)=476.0 cm－1、ν3(a1)=1506.3 cm－1.由微观过程的可逆性原理分析了分子的可能离解极限.并用多体项展式理论方法分别导出基态PdH2和YH2分子的势能函数,其等值势能面图准确地再现了PdH2和YH2分子的结构特征和离解能,由此讨论了Pd ＋ H2和Y ＋ H2分子反应的势能面静态特征.     

PuC和PuC2的分子结构与势能函数

采用密度泛函B3LYP方法和相对论有效原子实理论模型优化出PuC和PuC2分子稳定构型,其电子状态分别为X5Σ－和X5A2.PuC2分子为C2v构型,其∠CPuC=147.67°,平衡核间距Re=0.22819 nm, 离解能De=5.543 eV, 并计算出谐振动频率:ν1=61.736 cm－1、ν2=229.894 cm－1、ν3=305.582 cm－1.在此基础上,运用多体项展式理论方法,导出了基态PuC2分子的分析势能函数,该势能面准确地再现了PuC2分子的稳定结构,并根据势能面等值图讨论了PuC＋C反应和Pu＋C2反应的势能面静态特征.     

Al(2Pu)+H2(X1Σ+g)的分子反应动力学

基于多体展式方法所导出的ＡｌＨ２（Ｘ＾２Ａ１）分析势能函数,用准经典的Ｍｏｎｔｅ－Ｃａｒｌｏ轨迹法对Ａｌ（＾２Ｐｕ）＋Ｈ２（Ｘ＾１∑＾＋ｇ,ｕ＝ｊ＝０）的分子反应动力学过程进行了计算。结果表明,此反应的主产物为交换反应Ａｌ（＾２Ｐｕ）＋Ｈ２（Ｘ＾１∑＾＋ｇ,ｖ＝ｊ＝０）→ＡｌＨ（Ｘ＾１∑＾＋,Ｖ’,ｊ’）＋Ｈ（＾２Ｓｇ）的ＡｌＨ（Ｘ＾１∑＾＋,ｖ’,ｊ’）没有发现ＡｌＨ２（Ｘ＾２Ａ１）的络合物。而     

Quantum state-to-state dynamics for the quenching process of Br(2P1/2) + H2(v(i) = 0, 1, j(i) = 0)

Quantum state-to-state dynamics for the quenching process Br((2)P(1/2)) + H(2)(v(i) = 0, 1, j(i) = 0) → Br((2)P(3/2)) + H(2)(v(f), j(f)) has been studied based on two-state model on the recent coupled potential energy surfaces. It was found that the quenching probabilities have some oscillatory structures due to the interference of reflected flux in the Br((2)P(1/2)) + H(2) and Br((2)P(3/2)) + H(2) channels by repulsive potential in the near-resonant electronic-to-vibrational energy transfer process. The final vibrational state resolved integral cross sections were found to be dominated by the quenching process Br((2)P(1/2)) + H(2)(v) → Br((2)P(3/2)) + H(2)(v+1) and the nonadiabatic reaction probabilities for Br((2)P(1/2)) + H(2)(v = 0, 1, j(i) = 0) are quite small, which are consistent with previous theoretical and experimental results. Our calculated total quenching rate constant for Br((2)P(1/2)) + H(2)(v(i) = 0, j(i) = 0) at room temperature is in good agreement with the available experimental data.     

基态UC2分子的结构和势能函数

采用密度泛函理论(DFT)的B3LYP方法和相对论有效原子实势理论模型(RECP),对UC2分子可能的结构进行优化计算,得到UC2分子稳定构型为角形C-U-C(C2v);由微观可逆性原理,判断了UC2分子的离解极限;并且导出了基态UC2分子(X 5B1)的多体项展式势能函数,其势能面等值图展现了C-U-C(C2v)稳定结构;根据势能面等值图,讨论了C+UC(X 3П)反应和U+C2(X 1∑+g)反应的势能面静态特征.     

Proton formation in 2+1 resonance enhanced multiphoton excitation of HCl and HBr via (Omega=0) Rydberg and ion-pair states

Molecular beam cooled HCl was state selected by two-photon excitation of the V (1) summation operator(0(+)) [v=9,11-13,15], E (1) summation operator(0(+)) [v=0], and g (3) summation operator(-)(0(+)) [v=0] states through either the Q(0) or Q(1) lines of the respective (1,3) summation operator(0(+))<--<--X (1) summation operator(0(+)) transition. Similarly, HBr was excited to the V (1) summation operator(0(+)) [v=m+3, m+5-m+8], E (1) summation operator(0(+)) [v=0], and H (1) summation operator(0(+)) [v=0] states through the Q(0) or Q(1) lines. Following absorption of a third photon, protons were formed by three different mechanisms and detected using velocity map imaging. (1) H(*)(n=2) was formed in coincidence with (2)P(i) halogen atoms and subsequently ionized. For HCl, photodissociation into H(*)(n=2)+Cl((2)P(12)) was dominant over the formation of Cl((2)P(32)) and was attributed to parallel excitation of the repulsive [(2) (2)Pi4llambda] superexcited (Omega=0) states. For HBr, the Br((2)P(32))Br((2)P(12)) ratio decreases with increasing excitation energy. This indicates that both the [(3) (2)Pi(12)5llambda] and the [B (2) summation operator5llambda] superexcited (Omega=0) states contribute to the formation of H(*)(n=2). (2) For selected intermediate states HCl was found to dissociate into the H(+)+Cl(-) ion pair with over 20% relative yield. A mechanism is proposed by which a bound [A (2) summation operatornlsigma] (1) summation operator(0(+)) superexcited state acts as a gateway state to dissociation into the ion pair. (3) For all intermediate states, protons were formed by dissociation of HX(+)[v(+)] following a parallel, DeltaOmega=0, excitation. The quantum yield for the dissociation process was obtained using previously reported photoionization efficiency data and was found to peak at v(+)=6-7 for HCl and v(+)=12 for HBr. This is consistent with excitation of the repulsive A(2) summation operator(12) and (2) (2)Pi states of HCl(+), and the (3) (2)Pi state of HBr(+). Rotational alignment of the Omega=0(+) intermediate states is evident from the angular distribution of the excited H(*)(n=2) photofragments. This effect has been observed previously and was used here to verify the reliability of the measured spatial anisotropy parameters.     

Full-dimensional time-dependent wave packet dynamics of H2 + D2 reaction

Collision induced dissociation (CID), four center reaction (4C), and single exchange reaction (SE) in H(2) (v(1) = high) + D(2) (v(2) = low) were studied by means of time-dependent wave packet approach within a full-dimensional model. Initial state-selected total reaction probabilities for the three competitive processes have been computed on two realistic global potential energy surfaces of Aguado-Suárez-Paniagua and Boothroyd-Martin-Keogh-Peterson (BMKP) with the total angular momentum J = 0. The role of both vibrationally excited and rotationally excited reagents was examined by varying the initial vibrational and rotational states. The vibrational excitation of the hot diatom gives an enhancement effect on the CID process, while the vibrational excitation of the cold diatom gives an inhibition effect. The rotational excitation of both reagents has a significant effect on the reaction process. The 4C and SE probabilities are at least one order of magnitude smaller than the CID probabilities over the energy range considered. Isotope substitution effects were also studied by substituting the collider D(2) by H(2) and HD on the BMKP potential energy surfaces. The CID process is most efficient for the H(2) + D(2) combination and least efficient for the H(2) + H(2) combination and is different for the 4C and SE processes.     

A full dimensional time-dependent wave packet study for the H4 four-center, collision induced dissociation, and single exchange reactions: reaction probabilities for J=0

A time-dependent initial state selected wave packet method has been developed to study the H2(v(1)=10-11,j1=0)+H2'(v2=0,j2=0)-->HH'+HH' four-center (4C) reaction, and two other competing reactions: the H2+H2'-->H+H+H2' collision induced dissociation (CID) and the H2+H2'-->H+HH'+H' single exchange (SE) reaction, in full six dimensions. Initial state-specific total reaction probabilities for these three competing reactions are presented for total angular momentum J=0 and the effects of reagent vibration on reactions are examined. It is found that (a) the CID process is the dominant process over the whole energy range considered in this study, but the 4C and SE processes also have non-negligible probabilities; (b) the SE process has a lower threshold energy than the 4C process, but the SE probability increases slower than the 4C probability as collision energy increases; (c) the vibrational excitation of H2(v1) is much more efficient than translational motion for promoting these processes, in particular to the CID process.     

Ground state potential energy curve and dissociation energy of MgH

New high-resolution visible emission spectra of the MgH molecule have been recorded with high signal-to-noise ratios using a Fourier transform spectrometer. Many bands of the A 2Pi-->X 2Sigma+ and B' 2Sigma+-->X 2Sigma+ electronic transitions of 24MgH were analyzed; the new data span the v' = 0-3 levels of the A 2Pi and B'2Sigma+ excited states and the v'=0-11 levels of the X 2Sigma+ ground electronic state. The vibration-rotation energy levels of the perturbed A 2Pi and B' 2Sigma+ states were fitted as individual term values, while those of the X 2Sigma+ ground state were fitted using the direct-potential-fit approach. A new analytic potential energy function that imposes the theoretically correct attractive potential at long-range, and a radial Hamiltonian that includes the spin-rotation interaction were employed, and a significantly improved value for the ground state dissociation energy of MgH was obtained. The v'=11 level of the X 2Sigma+ ground electronic state was found to be the highest bound vibrational level of 24MgH, lying only about 13 cm(-1) below the dissociation asymptote. The equilibrium dissociation energy for the X 2Sigma+ ground state of 24MgH has been determined to be De=11104.7+/-0.5 cm(-1) (1.37681+/-0.00006 eV), whereas the zero-point energy (v'=0) is 739.11+/-0.01 cm(-1). The zero-point dissociation energy is therefore D0=10365.6+/-0.5 cm(-1) (1.28517+/-0.00006 eV). The uncertainty in the new experimental dissociation energy of MgH is more than 2 orders of magnitude smaller than that for the best value available in the literature. MgH is now the only hydride molecule other than H2 itself for which all bound vibrational levels of the ground electronic state are observed experimentally and for which the dissociation energy is determined with subwavenumber accuracy.     

A new analytical potential energy surface for the singlet state of He2H+

The analytic potential energy surface (APES) for the exchange reaction of HeH(+) (X(1)Σ(+)) + He at the lowest singlet state 1(1)A(∕) has been built. The APES is expressed as Aguado-Paniagua function based on the many-body expansion. Using the adaptive non-linear least-squares algorithm, the APES is fitted from 15 682 ab initio energy points calculated with the multireference configuration interaction calculation with a large d-aug-cc-pV5Z basis set. To testify the new APES, we calculate the integral cross sections for He + H(+)He (v = 0, 1, 2, j = 0) → HeH(+) + He by means of quasi-classical trajectory and compare them with the previous result in literature.     

Ab initio dynamic study of the reaction of Cl2LaR (R=H, CH3) with H2

In this paper, a comparison between "static" and "dynamic" determination of the thermodynamic (DeltarF degrees) and kinetic data (DeltarF#) for the reaction of Cl2LaR (R=H, CH3) and H2 is given. A difference is obtained in the case of the reaction between Cl2LaH and H2 and can be attributed to a failure of the "static" approach based on the harmonic approximation. The influence of the zero point energy correction is also analyzed but does not explain the 30% difference between the two calculated activation energies. The influence of the flatness of the potential energy surface around the transition state is proved as no such an effect is observed for the reaction of Cl2LaCH3 and H2.     

基态H2O+的全域多体势能面

基态H₂O⁺分子离子是重要的星际反应中间体.用MRCISD/cc-pVTZ先计算H₂O⁺势能面,然后拟合成多体势能函数.全域势能面能很好地描述基态H₂O⁺的行为.与两条反应通道的实验能变和其构型翻转反应能垒比较,表明计算结果与实验值相符.对计算点完成了力场的多项式拟合,对其振动模式进行了分析     

Ultracold collisions and reactions of vibrationally excited OH radicals with oxygen atoms

We report a quantum dynamics study of O + OH (v = 1, j = 0) collisions on its ground electronic state, employing two different potential energy surfaces: the DIMKP surface by Kendrick and Pack, and the XXZLG surface by Xu et al. A time-independent quantum mechanical method based on hyperspherical coordinates has been adopted for the dynamics calculations. Energy-dependent probabilities and rate coefficients are computed for the elastic, inelastic, and reactive channels over the collision energy range E(coll) = 10(-10)-0.35 eV, for J = 0 total angular momentum. Initial state-selected reaction rate coefficients are also calculated from the J = 0 reaction probabilities by applying a J-shifting approximation, for temperatures in the range T = 10(-6)-700 K. Our results show that the dynamics of the collisional process and its outcome are strongly influenced by long-range forces, and chemical reactivity is found to be sensitive to the choice of the potential energy surface. For O + OH (v = 1, j = 0) collisions at low temperatures, vibrational relaxation of OH competes with reactive scattering. Since long-range interactions can facilitate vibrational relaxation processes, we find that the DIMKP potential (which explicitly includes van der Waals dispersion terms) favours vibrational relaxation over chemical reaction at low temperatures. On the DIMKP potential in the ultracold regime, the reaction rate coefficient for O + OH (v = 1, j = 0) is found to be a factor of thirteen lower than that for O + OH (v = 0, j = 0). This significantly high reactivity of OH (v = 0, j = 0), compared to that of OH (v = 1, j = 0), is attributed to enhancement caused by the presence of a HO(2) quasibound state (scattering resonance) with energy near the O + OH (v = 0, j = 0) dissociation threshold. In contrast, the XXZLG potential does not contain explicit van der Waals terms, being just an extrapolation by a nearly constant function at large O-OH distances. Therefore, long-range potential couplings are absent in calculations using the XXZLG surface, which does not induce vibrational relaxation as efficiently as the DIMKP potential. The XXZLG potential leads to a slightly higher reactivity (a factor of 1.4 higher) for O + OH (v = 1, j = 0) compared to that for O + OH (v = 0, j = 0) at ultracold temperatures. Overall, both potential surfaces yield comparable values of reaction rate coefficients at low temperatures for the O + OH (v = 1, j = 0) reaction.     

Theoretical investigation on the electronic and geometric structure of GaN2+ and GaN4+

The electronic and geometric structures of gallium dinitride cation, GaN2+ and gallium tetranitride cation, GaN4+ were systematically studied by employing density functional theory (DFT-B3LYP) and perturbation theory (MP2, MP4) in conjunction with large basis sets, (aug-)cc-pVxZ, x = T, Q. A total of 7 structures for GaN2+ and 24 for GaN4+ were identified, corresponding to minima, transition states, and saddle points. We report geometries and dissociation energies for all the above structures as well as potential energy profiles, potential energy surfaces, and bonding mechanisms for some low-lying electronic states. The calculated dissociation energy (De) of the ground state of GaN2+, X1Sigma+, is 5.6 kcal/mol with respect to Ga+(1S) + N2(X1Sigmag+) and that of the excited state, ?3Pi, is 24.8 kcal/mol with respect to Ga+(3P) + N2(X1Sigmag+). The ground state and the first excited minimum of GaN4+ are of 1A1(C2v) and 3B1(C2v) symmetry with corresponding De of 11.0 and 43.7 kcal/mol with respect to Ga+(1S) + 2N2(X1Sigmag+) for X1A1 and Ga+(3P) + 2N2(X1Sigmag+) for 3B1.     

Homoleptic tetrahydrometalate anions MH(4)(-) (M = Sc,Y, La). Matrix infrared spectra and DFT calculations

Laser-ablated Sc, Y, and La atoms react with molecular hydrogen upon condensation in excess argon, neon, and deuterium to produce the metal dihydride molecules and dihydrogen complexes MH(2) and (H(2))MH(2). The homoleptic tetrahydrometalate anions ScH(4)(-), YH(4)(-), and LaH(4)(-) are formed by electron capture and identified by isotopic substitution (D(2), HD, and H(2) + D(2) mixtures). Doping with CCl(4) to serve as an electron trap virtually eliminates the anion bands, and further supports the anion identifications. The observed vibrational frequencies are in agreement with the results of density functional theory calculations, which predict electron affinities in the 2.8-2.4 eV range for the (H(2))ScH(2), (H(2))YH(2), and (H(2))LaH(2) complexes, and indicate high stability for the MH(4)(-) (M = Sc, La, Y) anions and suggest the promise of synthesis on a larger scale for use as reducing agents.     

State-to-state scattering of D2 from Cu(100) and Pd(111)

Results from state resolved experiments are presented for the interaction of D2(v=1,J=2) with Cu(100) and Pd(111). The reflected molecules were probed using quantum state specific spectroscopy. For D2 scattered from Cu(100) the vibrational survival probability and some transition inelastic probabilities were measured for incident energies from 70-200 meV. The survival probability was found to be larger then that found previously for H2(v=1) scattered from the same surface; these differences are discussed in terms of the lower zero point energy and smaller vibrational energy spacings of D2. D2 translational energy exchange was studied for several different scattering channels and interpreted using simple classical calculations. The survival probability was also measured for D2(v=1) scattered from Pd(111) at one incident energy. Pd is reactive for D2 dissociation and this survival probability was measured to be small and also to be much smaller than that for H2(v=1) under similar conditions. Vibrational relaxation channels were studied for D2 scattering from both Cu(100) and Pd(111). The vibrational relaxation probability on both surfaces was also found to be smaller than that measured for comparable channels for H2. The smaller survival probability and vibrational relaxation probability for D2 on Pd(111) cannot be easily accounted for by the difference in zero point energy and vibrational energy spacings.     

Quantum and classical studies of the O(3P) + H2(v = 0-3,j = 0) --> OH + H reaction using benchmark potential surfaces

We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.     

Exact quantum scattering study of the H + HS reaction on a new ab initio potential energy surface H2S (3A")

We present an exact quantum dynamical study and quasi-classical trajectory (QCT) calculations for the exchange and abstraction processes for the H + HS reaction. These calculations were based on a newly constructed high-quality potential energy surface for the lowest triplet state of H(2)S ((3)A"). The ab initio single-point energies were computed using complete active space self-consistent field and multi-reference configuration interaction method with a basis set of aug-cc-pV5Z. The time-dependent wave packet (TDWP) method was used to calculate the total reaction probabilities and integral cross sections over the collision energy (E(col)) range of 0.0-2.0 eV for the reactant HS initially at the ground state and the first vibrationally excited state. It was found that the initial vibrational excitation of HS enhances both abstraction and exchange processes. In addition, a good agreement is found between QCT and TDWP reaction probabilities at the total momentum J = 0 as a function of collision energy for the H + HS (v = 0, j = 0) reaction.