Physisorption and dissociative chemisorption of H2 on sub-nanosized Al$_{13}^{-}$ anion cluster: ab initio study

We have studied the interaction of Al₁₃^-_{13}^{-} anion cluster with H₂. Both the long range interaction and dissociative adsorption have been examined using the established correlated ab initio methods, MP2 and CCSD(T), in conjunction with the augmented correlation consistent basis sets up to aug-cc-pVTZ. The formation of the weakly bound (physisorbed) end-on anion complex Al₁₃^-_{13}^{-}...H₂ is predicted for the interacting Al...H distances of 3.95 ? with the H-H axis pointing towards the ‘hollow’ site of Al₁₃^-_{13}^{-} and binding energy (D_e)D_{e}) of 0.7 kcal/mol at the estimated complete basis set (CBS) limit of CCSD(T). The barrier height for H₂ dissociation on Al₁₃^-_{13}^{-} of 41.6 (42.9) kcal/mol calculated at the ZPVE-corrected CCSD(T)/aug-cc-pVTZ (estimated CCSD(T)/CBS) level is at least twice as large as that evaluated by us for a dissociative adsorption of H₂ on an open-shell Al₁₃ neutral cluster. To our knowledge, this report presents the first “benchmark” quality study of the physisorption and dissociative chemisorption of molecular hydrogen on Al₁₃^-_{13}^{-} anion cluster.     

The keto-enol equilibrium in substituted acetaldehydes: focal-point analysis and ab initio limit

High-level ab initio electronic structure calculations up to the CCSD(T) theory level, including extrapolations to the complete basis set (CBS) limit, resulted in high precision energetics of the tautomeric equilibrium in 2-substituted acetaldehydes (XH₂C-CHO). The CCSD(T)/CBS relative energies of the tautomers were estimated using CCSD(T)/aug-cc-pVTZ, MP3/aug-cc-pVQZ, and MP2/aug-cc-pV5Z calculations with MP2/aug-cc-pVTZ geometries. The relative enol (XHC?=?CHOH) stabilities (ΔE _{e,CCSD(T)/CBS}) were found to be 5.98?±?0.17, ?1.67?±?0.82, 7.64?±?0.21, 8.39?±?0.31, 2.82?±?0.52, 10.27?±?0.39, 9.12?±?0.18, 5.47?±?0.53, 7.50?±?0.43, 10.12?±?0.51, 8.49?±?0.33, and 6.19?±?0.18?kcal?mol^?1 for X?=?BeH, BH₂, CH₃, Cl, CN, F, H, NC, NH₂, OCH₃, OH, and SH, respectively. Inconsistencies between the results of complex/composite energy computations methods Gn/CBS (G2, G3, CBS-4M, and CBS-QB3) and high-level ab initio methods (CCSD(T)/CBS and MP2/CBS) were found. DFT/aug-cc-pVTZ results with B3LYP, PBE0 (PBE1PBE), TPSS, and BMK density functionals were close to the CCSD(T)/CBS levels (MAD?=?1.04?kcal?mol^?1).     

Calculation of accurate imaginary frequencies and tunnelling coefficients for hydrogen abstraction reactions using IRCmax

The effect of level of theory on the imaginary frequency and corresponding tunnelling coefficients has been studied for a test set of hydrogen abstraction reactions: ?CH₂X + CH₃Y → CH₃X + ?CH₂Y for (X,Y) = (H,H), (H,CN), (H,F), (H,Li) and (F,Li). It is found that the imaginary frequency is very sensitive to the level of theory used, with Hartree-Fock (HF) methods severely overestimating the imaginary frequency compared with high-level CCSD(T)/6-311G(d,p) calculations. The errors for the other methods are smaller but nonetheless significant, with MP2 methods overestimating the imaginary frequency and density functional theory (DFT) methods underestimating it. In the case of the HF methods, this leads to errors in the tunnelling coefficient of several orders of magnitude, while for the better DFT and MP2 methods errors of a factor of 2–3 are observed. To address this problem, an IRCmax procedure for estimating the imaginary frequency has been developed and it is found that IRCmax imaginary frequencies calculated with CCSD(T)/6-311G(d,p) single points along a low-level HF/6-31G(d) minimum energy path provide excellent approximations to the high-level values, at a fraction of the computational cost.     

Addition complex and insertion isomers on the potential energy surface of the reaction of indium dimer with water studied with relativistic ECP

Stationary points on the lowest singlet and triplet In₂ + H₂O potential energy surfaces (PESs) have been explored using the coupled cluster method, including single and double excitations with perturbative triples (CCSD(T)), and the density functional theory (DFT), employing the effective core potential (ECP) for indium (In), which accounts for scalar relativistic effects, with the triple-zeta quality basis set. The CCSD(T) calculated binding energy and anharmonic ν₂-bending mode frequency for the triplet ground-state addition complex, In₂^… OH₂(³B₁), are consistent with the complex detected in the matrix isolation infrared (IR) spectroscopic study under the thermal conditions. The two minimum energy crossing points between the triplet and the singlet PESs that have been located between the structures of In₂^…OH₂ and the transition state for the O–H bond breakage are not likely to be thermally accessible under the low-temperature matrix conditions. With the CCSD(T)-calculated In₂ + H₂O reaction profile and anharmonic vibrational frequencies for several In₂(H)(OH) insertion product isomers, we support the IR matrix isolation detection (by two experimental groups) of the lowest energy singlet double-bridged In(μ-H)(μ-OH)In isomer. For the proposed two-step mechanism of H₂ elimination from the In₂(H)(OH) species, the estimated energy barriers are also compatible with experiment.     

Reaction mechanism and rate constants of the CH+CH4 reaction: a theoretical study

Ab initio and density functional calculations have been performed to elucidate the mechanism of CH radical insertion into methane. The results show that the reaction can be viewed to occur via two stages. On the first stage, the CH radical approaches methane without large structural changes to acquire proper positioning for the subsequent stage, where H-migration occurs from CH₄ to CH, along with a C–C bond formation. Where the first stage ends and the second begins, a tight transition state was located using the B3LYP/6-311G(d,p) and MP4(SDQ)/6-311++G(d,p) methods. Using a rigid rotor – harmonic oscillator approach within transition state theory, we show that at the MP5/6-311++G(d,p)//MP4(SDQ)/6-311++G(d,p) level the calculated rate constants are in a reasonably good agreement with experiment in a broad temperature range of 145–581 K. Even at low temperatures, the insertion reaction bottleneck is found about the location of the tight transition state, rather than at long separations between the CH and CH₄ reactants. In addition, high level CCSD(T)-F12/CBS calculations of the remainder of the C₂H₅ potential energy surface predict the CH+CH₄ reaction to proceed via the initial insertion step to the ethyl radical which then can emit a hydrogen atom to form highly exothermic C₂H₄+H products.     

A high level computational study of the CH4/CF4 dimer: how does it compare with the CH4/CH4 and CF4/CF4 dimers?

The interaction within the methane–methane (CH₄/CH₄), perfluoromethane–perfluoromethane (CF₄/CF₄) methane–perfluoromethane dimers (CH₄/CF₄) was calculated using the Hartree–Fock (HF) method, multiple orders of Møller–Plesset perturbation theory [MP2, MP3, MP4(DQ), MP4(SDQ), MP4(SDTQ)], and coupled cluster theory [CCSD, CCSD(T)], as well as the PW91, B97D, and M06-2X density functional theory (DFT) functionals. The basis sets of Dunning and coworkers (aug-cc-pVxZ, x?=?D, T, Q), Krishnan and coworkers [6-311++G(d,p), 6-311++G(2d,2p)], and Tsuzuki and coworkers [aug(df, pd)-6-311G(d,p)] were used. Basis set superposition error (BSSE) was corrected via the counterpoise method in all cases. Interaction energies obtained with the MP2 method do not fit with the experimental finding that the methane–perfluoromethane system phase separates at 94.5?K. It was not until the CCSD(T) method was considered that the interaction energy of the methane–perfluoromethane dimer (?0.69?kcal?mol^?1) was found to be intermediate between the methane (?0.51?kcal?mol^?1) and perfluoromethane (?0.78?kcal?mol^?1) dimers. This suggests that a perfluoromethane molecule interacts preferentially with another perfluoromethane (by about 0.09?kcal?mol^?1) than with a methane molecule. At temperatures much lower than the CH₄/CF₄ critical solution temperature of 94.5?K, this energy difference becomes significant and leads perfluoromethane molecules to associate with themselves, forming a phase separation. The DFT functionals yielded erratic results for the three dimers. Further development of DFT is needed in order to model dispersion interactions in hydrocarbon/perfluorocarbon systems.     

Ab initio potential energy surface and intermolecular vibrational frequency of C3-He complex

The intermolecular potential energy surface for C₃–He complex has been constructed using supermolecular CCSD(T) and MP4 methods. The potential surfaces have been calculated for 27 values of R ranging from 2.8 to 8.0 Å and 19 values of θ equally spaced between 0° and 180°. Both CCSD(T) and MP4 potentials have similar global behaviors. The global minimum in each of the potentials corresponds to the slightly distorted T-shaped geometry. On the basis of these two potentials, the intermolecular vibrational energies and wavefunctions were calculated. The energy level pattern of the vdW vibrational states was predicted for C₃–He complex. The zero point bending motion of this complex has a range of 180°. The calculated fundamental frequency of vdW bending is 3.16 cm^?1 at CCSD(T) level, and 5.38 cm^?1 at the MP4 level. In addition, we have also constructed the intermolecular potential energy surface with C₃ bending coordinate of 160° by using supermolecular CCSD(T) method. Two local minima including arrow-shaped and Y-shaped configurations were determined. The rotational constants of three C₃–He structures including T-shaped, arrow-shaped and Y-shaped configurations at CCSD(T) level were also reported.     

Enthalpy difference between conformations of normal alkanes: effects of basis set and chain length on intramolecular basis set superposition error

The quantum chemistry of conformation equilibrium is a field where great accuracy (better than 100?cal?mol^?1) is needed because the energy difference between molecular conformers rarely exceeds 1000–3000?cal?mol^?1. The conformation equilibrium of straight-chain (normal) alkanes is of particular interest and importance for modern chemistry. In this paper, an extra error source for high-quality ab initio (first principles) and DFT calculations of the conformation equilibrium of normal alkanes, namely the intramolecular basis set superposition error (BSSE), is discussed. In contrast to out-of-plane vibrations in benzene molecules, diffuse functions on carbon and hydrogen atoms were found to greatly reduce the relative BSSE of n-alkanes. The corrections due to the intramolecular BSSE were found to be almost identical for the MP2, MP4, and CCSD(T) levels of theory. Their cancelation is expected when CCSD(T)/CBS (CBS, complete basis set) energies are evaluated by addition schemes. For larger normal alkanes (N?>?12), the magnitude of the BSSE correction was found to be up to three times larger than the relative stability of the conformer; in this case, the basis set superposition error led to a two orders of magnitude difference in conformer abundance. No error cancelation due to the basis set superposition was found. A comparison with amino acid, peptide, and protein data was provided.     

Benchmark calculations with correlated molecular wavefunctions. XIII. Potential energy curves for He2, Ne2 and Ar2 using correlation consistent basis sets through augmented sextuple zeta

The potential energy curves of the rare gas dimers He2, Ne2, and Ar2 have been computed using correlation consistent basis sets ranging from singly augmented aug-cc-pVDZ sets through triply augmented t-aug-cc-pV6Z sets, with the augmented sextuple basis sets being reported herein. Several methods for including electron correlation were investigated, namely Møller—Plesset perturbation theory (MP2, MP3 and MP4) and coupled cluster theory [CCSD and CCSD(T)]. For He2 CCSD(T)/d-aug-cc-pV6Z calculations yield a well depth of 7.35 cm-1 (10.58 K), with an estimated complete basis set (CBS) limit of 7.40 cm-1 (10.65 K). The latter is smaller than the 'exact' well depth (Aziz, R. A., Janzen, A. R., and Moldover, M. R., 1995, Phys. Rev. Lett., 74, 1586) by about 0.2 cm-1 (0.35 K). The Ne2 well depth, computed with the CCSD(T)/d-aug-cc-pV6Z method, is 28.31 cm-1 and the estimated CBS limit is 28.4 cm-1, approximately 1 cm-1 smaller than the empirical potential of Aziz, R. A., and Slaman, M., J., 1989, Chem. Phys., 130, 187. Inclusion of core and core—valence correlation effects has a negligible effect on the Ne2 well depth, decreasing it by only 0.04 cm-1. For Ar2, CCSD(T)/d-aug-cc-pV6Z calculations yield a well depth of 96.2 cm-1. The corresponding HFDID potential of Aziz, R. A., 1993, J. chem. Phys., 99, 4518 predicts of De of 99.7 cm-1. Inclusion of core and core-valence effects in Ar2 increases the well depth and decreases the discrepancy by approximately 1 cm-1.     

Effect of CO and H adsorption on the compositional structure of binary nanoalloys via DFT modeling

Theoretical methods (DFT/B3LYP and MP2) have been used to optimize the geometries of the Al _x Sb _y (x + y = 3, 5) clusters and their anions. Single point energy computations at CCSD(T) level have also been performed using the optimized B3LYP and MP2 structures. The basis sets used for Al and Sb atoms are 6-311+G(2d) and LANL2DZdp ECP, respectively. Harmonic vibrational frequency computations were carried out to confirm the nature of the stationary points. We report the structural and spectroscopic parameters of the named clusters. We also report the relative energy of the clusters, the vertical electron detachment energy, the adiabatic electron detachment energy and the adiabatic electron affinity. The most stable structures at the CCSD(T)//MP2 level are, the D _∞h linear structure (AlSb₂) and the C _2v V-bent structure (AlSb ₂ ^? ), the C _2v V-bent structure (Al₂Sb and its anion), the C _2v edge-capped tetrahedron (Al₂Sb₃ and its anion), the C _2v trigonal bipyramidal structure (Al₃Sb₂ and its anion), the C _4v square pyramidal (AlSb₄) and a C _2v ground structure for its anion, the C _2v planar trapezoidal structure (Al₄Sb) and the C _2v edge-capped tetrahedron (Al₄Sb^?). The adiabatic electron affinities calculated at the CCSD(T)//MP2 level are 2.17 eV (AlSb₂), 2.17 eV (Al₂Sb), 2.38 eV (Al₂Sb₃), 2.76 eV (Al₃Sb₂), 2.21 eV (AlSb₄) and 2.03 eV (Al₄Sb). The findings of this research are analysed, discussed and compared with the analogous picnogenides clusters.     

Ab initio potential energy surface and intermolecular vibrational frequencies of C3–Ar complex

The intermolecular potential energy surfaces for C₃–Ar have been calculated by supermolecular CCSD(T) and MP4 methods. The MP4 and CCSD(T) potentials have similar global behaviours. Their global minima all correspond to the slightly distorted T-shaped geometries. From these two potentials, the intermolecular vibrational energies and wavefunctions were calculated. The energy level patterns of the vdW vibrational states were predicted for the C₃–Ar complex. The zero point bending motion of this complex has a range of approximately 60°. The calculated transition frequencies of vdW bending agree well with available experimental data.     

Theoretical study of B3O radical isomers and their interconversion pathways

Despite the fact that B₃O is the second simplest B _n O radical after BO, a controversy recently emerged concerning the molecular structure of its global minimum. Two recent theoretical groups predicted the linear quartet BBBO to be the ground isomer. By contrast, another recent theoretical group reported that B₃O has a doublet B₃-ring ground structure. Moreover, larger B _n O clusters usually have low-lying B₃-ring isomers. In order to determine the accurate energetic competition between linear and cyclic structures in both the doublet and quartet, and to understand the detailed isomerism between various isomers, which is vital for understanding the formation mechanism of B₃O, we report the first potential energy surface (PES) study of B₃O at various computational levels, including CCSD(T)/6-311+G(2df), CCSD(T)/aug-cc-pVTZ, CCSD(T)/aug-cc-pVQZ and G3B3 for the single-point energy, as well as B3LYP/6-311+G(d) and QCISD/6-311+G(d) for geometrical optimisation. It is shown that the isomers in the quartet state are all thermodynamically more stable than the corresponding doublet ones, and on both the quartet and doublet PESs, the linear form has the lowest energy. Therefore, our study on both linear and cyclic isomers shows that the linear quartet BBBO ⁴ 01 is definitively the ground isomer. Although being much less stable than the quartet linear BBBO global minimum by >20 kcal mol^?1, five cyclic isomers exist as local minima, with the bi-cyclic structure ⁴ 02 possessing the smallest barrier of around 15 kcal mol^?1. The dissociation energies for direct combination processes B₃ + O, B₂ + BO and B + B₂O are discussed. The present work may be helpful in obtaining a deep understanding of the doping and oxidation process of pure B _n clusters.     

硼、氮取代二聚苯中碳原子导致层间相互作用本质变化研究

在MP2/aug-cc-pVTZ水平下优化了二聚苯及其碳被硼、氮原子取代的几何构型;计算了这些二聚体的CCSD(T)/CBS相互作用能;并用SAPT2+/aug-cc-pVDZ分析了相互作用能成分.探索了硼、氮取代二聚苯中碳原子所导致层间相互作用本质变化.结果表明:稳定构型而言,从苯(C_6H_6)、1-氮-2-硼杂苯(BNC_4H_6)、1,3-二氮-2,4-二硼杂苯(B_2N_2C_2H_6)二聚体的平行移位(PD)构型向无机苯(B_3N_3H_6)二聚体的夹心(S)构型转变,其中C_6H_6和B_3N_3H_6二聚体稳定构型的堆叠型式分别与之相对应的石墨烯(GE)和六方氮化硼(h-BN)2D层间材料的堆叠型式相一致.硼、氮原子取代二聚苯中的碳原子后使其相互作用能增大,其中BNC_4H_6和B_2N_2C_2H_6二聚体的相互作用能增大较为明显.所研究体系二聚体稳定构型均以色散能为主导、静电能次之、诱导能相对较小.硼、氮取代二聚苯中碳原子后其静电能对总吸引能的贡献明显增大.     

Density functional study of structural and electronic properties of Al<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript>P (2 ≤ <Emphasis Type="Italic">n</Emphasis> ≤ 12) clusters

Low-lying equilibrium geometric structures of Phosphorus-doped aluminum cluster Al _n P (n = 2–12) clusters obtained by an all-electron linear combination of atomic orbital approach, within spin-polarized density functional theory, are reported. The binding energy, dissociation energy, and stability of these clusters are studied within the local spin density approximation (LSDA) and the three-parameter hybrid generalized gradient approximation (GGA) due to Becke-Lee-Yang-Parr (B3LYP). Ionization potentials, electron affinities, hardness, and static polarizabilities are calculated for the ground-state structures within the GGA. It is observed that symmetric structures with the P atom occupying a peripheral position are lowest-energy geometries of Al _n P (n = 2, 4–11), while the P impurities of Al₃P and Al₁₂P prefer to occupy internal sites in the aluminum clusters. Generalized gradient approximation extends bond lengths as compared to the LSDA lengths. The odd-even oscillations in the dissociation energy, the second differences in energy, the HOMO–LUMO gaps, the ionization potential, the electron affinity, and the hardness are more pronounced within both GGA and LSDA. The stability analysis based on the energies clearly shows the clusters with an even number of valence electrons are more stable than clusters with odd number of valence electrons.     

Molecular structure and IR spectra of bromomethanes by DFT and post-Hartree-Fock MP2 and CCSD(T) calculations

The molecular parameters (geometries, rotational constants, dipole moments) and vibrational IR spectra (harmonic wavenumbers, absolute intensities) of bromomethanes (CH₃Br, CH₂Br₂, CHBr₃, CBr₄) are predicted by a density functional theory with the hybrid Becke3-LYP functional (DFT) and post-Hartree-Fock methods (MP2, CCSD(T)) using a 6-311G(2d,2p)-type basis set. The MP2 calculations are carried out with different numbers of frozen core orbitals to find how the number of bromine orbitals used for electron correlation influences the predicted molecular parameters and IR spectra of the species in question. Three options were used: (a) all electrons (full), with both the core and valence orbitals considered; (b) partial frozen core option (pfc), when the orbitals up to 3p of bromine were frozen; and (c) full frozen core option (ffc), when all core orbitals up to 3d were frozen. The CCSD(T) calculations for geometric parameters were carried out with both the pfc and ffc options, while for the prediction of the IR spectra only the ffc option was used. In addition, the calculations at the DFT and MP2(pfc) levels with inclusion of f functions on carbon and bromine atoms in bromomethanes (and also the CCSD(T)(pfc) calculations for CH₃Br) were carried out to predict the changes in the geometric parameters and/or vibrational IR spectra of the molecules upon inclusion of f functions The geometries of bromomethanes (particularly the CBr bond lengths) are predicted better by the DFT and CCSD(T) calculations when the f functions (in particular on bromine atom) are included, while the MP2 calculations without f functions are good enough for correct predictions of the molecular geometries. The molecular parameters and vibrational IR spectra of bromomethanes in question and their deuterated species predicted by the DFT, MP2(ffc) and CCSD(T)(ffc) with the 6-311G(2d,2p) basis set agree well with the available experimental data.     

Theoretical investigation of the oxidation pathways of the Cl-initiated reaction of 2-methyl-3-buten-2-ol

The mechanism and products of the reaction of 2-methyl-3-buten-2-ol (MBO232) with Cl atoms in the presence of O₂ have been elucidated by performing high-level quantum chemistry calculations. The geometries of the reactants, intermediates, transition states, and products are optimized at the MP2(full)/6-311G(d,?p) level, and their single-point energies are refined at the CCSD(T)/6-311?+?G(d,?p) level. The potential energy surface profiles have been constructed at the CCSD(T)/6-311?+?G(d,?p)//MP2(full)/6-311G(d,?p)?+?0.95?×?ZPE level of theory, and the possible channels involved in the reaction are also discussed. The calculations indicate that the reaction predominantly proceeds via the addition of Cl atoms to the double bond rather than the direct abstraction of the H atoms in MBO232. The nascent adducts (CH₃)₂C(OH)CHCH₂Cl (IM1) and (CH₃)₂C(OH)CHClCH₂ (IM2) do not undergo subsequent isomerization and dissociation reactions, but rather react with O₂. The theoretical results show that the major products are CH₂ClCHO and CH₃C(O)CH₃ for the reaction of MBO232?+?Cl in the presence of O₂, which is in good agreement with the experimental finding.     

Structures and stabilities of CBHz (z ≤ 8) and CxB3-xHz (x = 1, 2, z ≤ 14): prediction of novel organo-boron radicals

Hydrogenated boron-carbon clusters, i.e. organo-borons, have received considerable attention both theoretically and experimentally. Herein, using a topology searching strategy, we systematically explore the structures and stabilities of small organo-borons with CBH_z (z ≤ 8) and C_xB_3-xH_z (x = 1, 2,z≤ 14) stoichiometry, with particular interests in the intrinsic stabilities of the organo-boron radicals. At the CCSD(T)/aug-cc-pVQZ//B3LYP/aug-cc-pVTZ and CCSD(T)/aug-cc-pVQZ//MP2/aug-cc-pVTZ levels, the stabilities of these global minimum organo-boron species were evaluated by considering dissociation pathways and the binding energy per atom. Aside from the five already studied radicals (CBH₂, CBH₄, C₂BH₂, C₂BH₄ and CB₂H₃), we predict six novel radicals, i.e. CBH₆, C₂BH₆, C₂BH₈, CB₂H, CB₂H₅ and CB₂H₇, which could be detected under suitable circumstances. However, observation of the highly hydrogenated CB₂H₉ radical is much less likely due to its minute stability towards H-extrusion. The computationally determined stable/meta-stable maximum hydrogenation numbers for CB, C₂B and CB₂ (6, 8 and 8, respectively) are in excellent agreement with a simple electron-counting model for C_xB_y chains. The newly predicted organo-boron radicals await future laboratory verification.     

The hydrogenated aluminium trimer: a theoretical examination of the formation and interconversion pathways

Five doublet isomers of the Al₃H₂ cluster lying within a narrow range of 5 kcal/mol, along with the isomerization transition states connecting them, have been located with the coupled-cluster CCSD(T) and DFT methods. The two most stable doublet structures, the C_2v planar including the two Hs bound terminally and C₁ non-planar showing one H in terminal site and the other in threefold site are found to be essentially degenerate. Although the reaction of Al₃ with H₂ to yield Al₃H₂ is found to be significantly exothermic, by 23.5 kcal/mol, this hydrogenation is impeded by a considerable kinetic barrier of 16 kcal/mol. Our result is consistent with the observed lack of reactivity of Al_n towards H₂(D₂) for n=3 under thermal conditions [3]. The quartet Al₃H₂ isomers are predicted to lie 16–21 kcal/mol higher in energy than the doublet analogues. Further dimerization of Al₃H₂ to form Al₆H₄ has also been examined. Electronic supplementary material Supplementary Online Material     

BeH，BeD，BeT分子基态(X2Σ+)的结构与势能函数

采用量子力学从头算方法，运用二次组态相互作用QCISD(T)/aug-cc-pVTZ和电子相关单双耦合簇CCSD(T)/6-311++G(3df，2pd)研究了BeH，BeD，BeT分子基态的结构与势能函数，计算出了这些分子的光谱数据(ω_e，ω_eχ_e，B_e，α_e，D_e)，结果与实验光谱数据吻合较好.这表明上述分子基态的势能函数可用经修正的Murrell-Sorbie+c₆函数来表示. 关键词： BeH BeD BeT分子基态 分子结构 势能函数