Very accurate potential energy curve of the LiH molecule

We present very accurate calculations of the ground-state potential energy curve (PEC) of the LiH molecule performed with all-electron explicitly correlated Gaussian functions with shifted centers. The PEC is generated with the variational method involving simultaneous optimization of all Gaussians with an approach employing the analytical first derivatives of the energy with respect to the Gaussian nonlinear parameters (i.e., the exponents and the coordinates of the shifts). The LiH internuclear distance is varied between 1.8 and 40 bohrs. The absolute accuracy of the generated PEC is estimated as not exceeding 0.3 cm(-1). The adiabatic corrections for the four LiH isotopologues, i.e., (7)LiH, (6)LiH, (7)LiD, and (6)LiD, are also calculated and added to the LiH PEC. The aforementioned PECs are then used to calculate the vibrational energies for these systems. The maximum difference between the computed and the experimental vibrational transitions is smaller than 0.9 cm(-1). The contribution of the adiabatic correction to the dissociation energy of (7)LiH molecule is 10.7 cm(-1). The magnitude of this correction shows its importance in calculating the LiH spectroscopic constants. As the estimated contribution of the nonadiabatic and relativistic effects to the ground state dissociation energy is around 0.3 cm(-1), their inclusion in the LiH PEC calculation seems to be the next most important contribution to evaluate in order to improve the accuracy achieved in this work.     

Non-Born-Oppenheimer variational calculation of the ground-state vibrational spectrum of LiH+

Very accurate, rigorous, variational, non-Born-Oppenheimer (non-BO) calculations have been performed for the fully symmetric, bound states of the LiH(+) ion. These states correspond to the ground and excited vibrational states of LiH(+) in the ground (2)Sigma(+) electronic state. The non-BO wave functions of the states have been expanded in terms of spherical N-particle explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance and 5600 Gaussians were used for each state. The calculations that, to our knowledge, are the most accurate ever performed for a diatomic system with three electrons have yielded six bound states. Average interparticle distances and nucleus-nucleus correlation function plots are presented.     

Quasi-classical trajectory study of the adiabatic reactions occurring on the two lowest-lying electronic states of the LiH2+ system

Quasi-classical trajectory calculations have been performed on the adiabatically allowed reactions taking place on the two lowest-lying electronic states of the LiH2+ system, using the ab initio potential energy surfaces of Martinazzo et al. (J. Chem. Phys., 2003, 119, 11 241). These reactions comprise: (i) the exoergic H2 and H2+ formation occurring through LiH+ + H and LiH + H+ collisions in the ground and in the first electronically excited state, respectively; (ii) the endoergic (ground state) LiH+ dissociation induced by collisions with H atoms; and (iii) the endoergic (excited state) Li + H2+ --> LiH + H+ reaction. The topic is of relevance for a better understanding of the lithium chemistry in the early universe. Thermal rate constants for the above reactions have been computed in the temperature range 10-5000 K and found in reasonably good agreement with estimates based on the capture model.     

On the ground state of titanium phosphide, TiP: a theoretical investigation

Using multireference variational and coupled cluster methods in conjunction with very large core-correlation-consistent basis sets, we have confirmed that the ground state of TiP is of (2)Sigma(+) symmetry with the first excited state A (2)Delta no more than 3.5 kcal/mol higher. We also report full potential energy curves, dissociation energies, bond lengths, dipole moments, and the usual spectroscopic constants.     

自优化剩余函数量子Monte Carlo方法

提出剩余函数量子Monte Carlo的一个新算法，这是一个自优化和自改善的过程.与以前的算法相比，本算法中的试探函数的优化是在剩余函数方法中同步进行的，而不是在变分Monte Carlo计算之前.为了优化试探函数，使用一种改进了的速降法，这是一个步长能够自动调节，超线性收敛的优化技术.在这个算法中,还使用了一种新的相关函数，它满足电子与电子以及电子与核奇点条件.此方法已被用于计算H2、LiH、Li2、H2O分子的基态以及CH2的X 3B1态、1 1A1态和2 1A1态的能量值.     

Structural and spectroscopic study of the LiRb molecule beyond the Born-Oppenheimer approximation

Adiabatic and diabatic potential energy curves and the permanent and transition dipole moments of the low-lying electronic states of the LiRb molecule dissociating into Rb(5s, 5p, 4d, 6s, 6p, 5d, 7s, 6d) + Li(2s, 2p) have been investigated. The molecular calculations are performed with an ab initio approach based on nonempirical pseudopotentials for Rb(+) and Li(+) cores, parametrized l-dependent core polarization potentials and full configuration interaction calculations. The derived spectroscopic constants (R(e), D(e), T(e), ω(e), ω(e)x(e), and B(e)) of the ground state and lower excited states are in good agreement with the available theoretical works. However, the 8-10(1)Σ(+), 8-10(3)Σ(+), 6(1,3)Π, and 3(1,3)Δ excited states are studied for the first time. In addition, to the potential energy, accurate permanent and transition dipole moments have been determined for a wide interval of internuclear distances. The permanent dipole moment of LiRb has revealed ionic characters both relating to electron transfer and yielding Li(-)Rb(+) and Li(+)Rb(-) arrangements. The diabatic potential energy for the (1,3)Σ(+), (1,3)Π, and (1,3)Δ symmetries has been performed for this molecule for the first time. The diabatization method is based on variational effective Hamiltonian theory and effective metric, where the adiabatic and diabatic states are connected by an appropriate unitary transformation.     

Optimization and application of lithium parameters for the reactive force field, ReaxFF

To make a practical molecular dynamics (MD) simulation of the large-scale reactive chemical systems of Li-H and Li-C, we have optimized parameters of the reactive force field (ReaxFF) for these systems. The parameters for this force field were obtained from fitting to the results of density functional theory (DFT) calculations on the structures and energy barriers for a number of Li-H and Li-C molecules, including Li(2), LiH, Li(2)H(2), H(3)C-Li, H(3)C-H(2)C-Li, H(2)C=C-LiH, HCCLi, H(6)C(5)-Li, and Li(2)C(2), and to the equations of state and lattice parameters for condensed phases of Li. The accuracy of the developed ReaxFF was also tested by comparison to the dissociation energies of lithium-benzene sandwich compounds and the collision behavior of lithium atoms with a C(60) buckyball.     

Dissociation energy of the water dimer from quantum Monte Carlo calculations

We report a study of the electronic dissociation energy of the water dimer using quantum Monte Carlo techniques. We have performed variational quantum Monte Carlo and diffusion quantum Monte Carlo (DMC) calculations of the electronic ground state of the water monomer and dimer using all-electron and pseudopotential approaches. We have used Slater-Jastrow trial wave functions with B3LYP type single-particle orbitals, into which we have incorporated backflow correlations. When backflow correlations are introduced, the total energy of the water monomer decreases by about 4-5 mhartree, yielding a DMC energy of -76.428 30(5) hartree, which is only 10 mhartree above the experimental value. In our pseudopotential DMC calculations, we have compared the total energies of the water monomer and dimer obtained using the locality approximation with those from the variational scheme recently proposed by Casula [Phys. Rev. B 74, 161102(R) (2006)]. The time step errors in the Casula scheme are larger, and the extrapolation of the energy to zero time step always lies above the result obtained with the locality approximation. However, the errors cancel when energy differences are taken, yielding electronic dissociation energies within error bars of each other. The dissociation energies obtained in our various all-electron and pseudopotential calculations range between 5.03(7) and 5.47(9) kcalmol and are in good agreement with experiment. Our calculations give monomer dipole moments which range between 1.897(2) and 1.909(4) D and dimer dipole moments which range between 2.628(6) and 2.672(5) D.     

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H(5)O(2) (+) using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh et al. [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.     

Vibrational energies of LiH2(+) and LiD2(+) in the Ã1Σ+ electronic state

In connection with the recent study of the ground electronic state of the LiH2(+) molecular ion (Kraemer, W. P.; Spirko, V. Chem. Phys. 2006, 330, 190), the adiabatic three-dimensional double-minimum potential enery surface of the first excited electronic state was evaluated, including its two lowest atom-diatom dissociation channels as well as the three-atom complete fragmentation asymptote. Applying the Sutcliffe-Tennyson Hamiltonian for triatomic molecules, the levels of all bound vibrational states and the levels of the states localized in the two energy minimum regions were separately determined. The validity of statistical methods such as the density of states approach and the nearest-neighbor level spacing distribution (NNSD) was tested for the light LiH2(+) ion. Special effort was put into investigating possible effects of a tunnelling motion across the proton-transfer barrier on the vibrational level pattern using the NNSD approach.     

Relativistic calculations of ground and excited states of LiYb molecule for ultracold photoassociation spectroscopy studies

We report a series of quantum-chemical calculations for the ground and some of the low-lying excited states of an isolated LiYb molecule by the spin-orbit multistate complete active space second-order perturbation theory (SO-MS-CASPT2). Potential energy curves, spectroscopic constants, and transition dipole moments (TDMs) at both spin-free and spin-orbit levels are obtained. Large spin-orbit effects especially in the TDMs of the molecular states dissociating to Yb((3)P(0,1,2)) excited states are found. To ensure the reliability of our calculations, we test five types of incremental basis sets and study their effect on the equilibrium distance and dissociation energy of the ground state. We also compare CASPT2 and CCSD(T) results for the ground state spectroscopic constants at the spin-free relativistic level. The discrepancies between the CASPT2 and CCSD(T) results are only 0.01 ? in equilibrium bond distance (R(e)) and 200 cm(-1) in dissociation energy (D(e)). Our CASPT2 calculation in the supermolecular state (R=100 a.u.) with the largest basis set reproduces experimental atomic excitation energies within 3% error. Transition dipole moments of the super molecular state (R=100 a.u.) dissociating to Li((2)P) excited states are quite close to experimental atomic TDMs as compared to the Yb((3)P) and Yb((1)P) excited states. The information obtained from this work would be useful for ultracold photoassociation experiments on LiYb.     

用变分Monte Carlo方法处理分子

对变分量子Ｍｏｎｔｅ Ｃａｒｔｏ方法提出了一种种算法：将传统的Ｈａｒｔｒｅｅ－Ｆｏｅｋ方法与量子Ｍｏｎｔｅ Ｃａｒｌｏ方法有机结合在一起；导出了“局部能”的解析式；使用了一种新的相关函数和随机数发生器。我们用这个新算法计算了Ｈ２、ＬｉＨ、Ｌｉ２、Ｈ２Ｏ、Ｆ２分子的基态和ＣＨ２分子的＾３Ｂ１、＾１Ａ１态的能量。计算结果表明，这个新算法在精度和统计误差两个方面比一般ＶＭＣ过程都要好得多。     

Direct and inverse reactions of LiH+ with He(1S) from quantum calculations: mechanisms and rates

The gas-phase reaction of LiH(+) (X(2)Σ) with He((1)S) atoms, yielding Li(+)He with a small endothermicity for the rotovibrational ground state of the reagents, is analysed using the quantum reactive approach that employs the Negative Imaginary Potential (NIP) scheme discussed earlier in the literature. The dependence of low-T rates on the initial vibrational state of LiH(+) is analysed and the role of low-energy Feshbach resonances is also discussed. The inverse destruction reaction of LiHe(+), a markedly exothermic process, is also investigated and the rates are computed in the same range of temperatures. The possible roles of these reactions in early universe astrophysical networks, in He droplets environments or in cold traps are briefly discussed.     

Accurate Ab Initio Calculations for LiH and its Ions, LiH+ and LiH−

Ab initio calculations were performed for LiH using a pseudopotential approach with CPP corrections and huge basis sets on both atoms. A wide range of ^1,3Σ⁺ electronic adiabatic states have been investigated, from the ground state up to those dissociating into Li(5p)+H. Permanent and transition electric dipole moments are also considered for the first few excited states. Comparison with experiments and recent all-electron calculations, reveals an excellent global accuracy, only the bottom of the ground state being better described by all-electron approaches. Using almost identical basis sets, coupled cluster all-electron calculations are performed for the ground states of LiH⁺, LiH⁻ and LiH. High care has been given to the correct relative position of the asymptotes, allowing for this rather complete set of accurate ab initio data to be useful for further molecular physics studies.     

Hartree-Fock-Heitler-London method. 2. First and second row diatomic hydrides

The Hartree-Fock-Heitler-London, HF-HL, method is a new ab initio approach which variationally combines the Hartree-Fock, HF, and the Heitler-London, HL, approximations, yielding correct dissociation products. Furthermore, the new method accounts for nondynamical correlation and explicitly considers avoided crossing. With the HF-HL model we compute the ground-state potential energy curves for H2 [1Sigma+g], LiH [X 1Sigma+], BeH [2Sigma+], BH [1Sigma+], CH [2Pi], NH [3Sigma-], OH [2Pi], and FH [1Sigma+], obtaining in average 80% of the experimental binding energy with a correct representation of bond breaking. Inclusion of ionic configurations improves the computed binding energy. The computed dipole moment is in agreement with laboratory data. The dynamical and nondynamical correlation energies for atomic and molecular systems with 2-10 electrons are analyzed. For BeH the avoided crossing of the two lowest [2Sigma+] states is considered in detail. The HF-HL function is proposed as the zero-order reference wave function for molecular systems. To account for the dynamical correlation energy a post-HF-HL technique based on multiconfiguration expansions is presented. We have computed the potential energy curves for H2 [1Sigma+g], HeH [2Sigma+], LiH [X1Sigma+], LiH [A1Sigma+], and BeH [2Sigma+]. The corresponding computed binding energies are 109.26 (109.48), 0.01 (0.01), 57.68 (58.00), 24.19 (24.82), and 49.61 (49.83) kcal/mol, with the experimental values given in parentheses. The corresponding total energies are -1.1741, -3.4035, -8.0695, -7.9446, and -15.2452 hartrees, respectively, the best ab initio variational published calculations, H2 excluded.     

An accurate non-Born-Oppenheimer calculation of the first purely vibrational transition in LiH molecule

In this work we study the ground and the first vibrationally excited states of LiH molecule. We performed an extensive nonrelativistic variational calculations of the two states without using the Born-Oppenheimer approximation. The results are analyzed and compared with the data extracted from recent experiments. The 0<--1 transition energy obtained in the calculations converged to a value which is less than a wave number above the transition energy estimated from the available experimental data concerning the LiH rovibrational transitions. We discuss the remaining discrepancy and the procedure used to determine the "experimental" transition frequencies.     

Electron correlation described by extended geminal models: The EXGEM4 and EXGEM5 models

Within the framework of the general extended geminal model, two new approximate models EXGEM 4 and EXGEM 5 are introduced. The models are tested against full CI calculations on the water molecule for three different nuclear configurations and a full CI potential energy curve for the LiH molecule in the ground state. On the basis of these calculations, it is suggested that the models will yield electronic correlation energies with an accuracy of 1–2% of the corresponding full CI result.     

Stability and structure of Li(n)H molecules (n=3-6): experimental and density functional study

The molecules Li(3)H and Li(4)H have been identified in mass-spectrometric measurements over solutions of hydrogen in liquid Li, and the gaseous equilibria of the reactions: Li(3)H+Li=Li(2)H+Li(2), Li(3)H+Li(2)=Li(2)H+Li(3), Li(3)H+Li=LiH+Li(3), Li(3)H+LiH=2Li(2)H, and Li(4)H+Li(2)=Li(3)H+Li(3) have been measured. Density functional calculations of Li(n)H molecules (n=3-6) provide structures, vibrational frequencies, ionization energies, and free energy functions of these molecules, and these are used to estimate the enthalpies of these reactions and the atomization energies of Li(3)H (119.4 kcal/mol) and Li(4)H (151.8 kcal/mol).     

Energies of the first row atoms from quantum Monte Carlo

All-electron variational and diffusion quantum Monte Carlo calculations of the ground state energies of the first row atoms (from Li to Ne) are reported. The authors use trial wave functions of four types: single-determinant Slater-Jastrow wave functions, multideterminant Slater-Jastrow wave functions, single-determinant Slater-Jastrow wave functions with backflow transformations, and multideterminant Slater-Jastrow wave functions with backflow transformations. At the diffusion quantum Monte Carlo level and using their multideterminant Slater-Jastrow wave functions with backflow transformations, they recover 99% or more of the correlation energies for Li, Be, B, C, N, and Ne, 97% for O, and 98% for F.     

The effect of metal cations on the nature of the first electronic transition of liquid water as studied by attenuated total reflection far-ultraviolet spectroscopy

The first electronic transition (?←X?) of liquid water was studied from the perspective of the hydration of cations by analyzing the attenuated total reflection far-ultraviolet (ATR-FUV) spectra of the Group I, II, and XIII metal nitrate electrolyte solutions. The ?←X? transition energies of 1 M electrolyte solutions are higher (Li(+): 8.024 eV and Cs(+): 8.013 eV) than that of pure water (8.010 eV) and linearly correlate with the Gibbs energies of hydration of the cations. The increases in the ?←X? transition energies are mostly attributable to the hydrogen bond formation energies of water molecules in the ground state induced by the presence of the cations. The deviation from the linear relation was observed for the high charge density cations, H(+), Li(+), and Be(2+), which reflects that the electronic energies in the excited states are also perturbed. Quantum chemical calculations show that the ?←X? transition energies of the water-cation complexes depend on the hydration structures of the cations. The calculated ?←X? transition energies of the water molecules hydrating high charge density cations spread more widely than those of the low charge density cations. The calculated transition energy spreads of the water-cation complexes directly correlate with the widths of the ?←X? transition bands measured by ATR-FUV spectroscopy.