Time-dependent density functional study of the electronic potential energy curves and excitation spectrum of the oxygen molecule

Orbital energies, ionization potentials, molecular constants, potential energy curves, and the excitation spectrum of O(2) are calculated using time-dependent density functional theory (TDDFT) with Tamm-Dancoff approximation (TDA). The calculated negative highest occupied molecular orbital energy (-epsilon(HOMO)) is compared with the energy difference ionization potential for five exchange correlation functionals consisting of the local density approximation (LDAxc), gradient corrected Becke exchange plus Perdew correlation (B(88X)+P(86C)), gradient regulated asymptotic correction (GRAC), statistical average of orbital potentials (SAOP), and van Leeuwen and Baerends asymptotically correct potential (LB94). The potential energy curves calculated using TDDFT with the TDA at internuclear distances from 1.0 to 1.8 A are divided into three groups according to the electron configurations. The 1pi(u) (4)1pi(g) (2) electron configuration gives rise to the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) states; the 1pi(u) (3)1pi(g) (3) electron configuration gives rise to the c (1)Sigma(u) (-), C (3)Delta(u), and A (3)Sigma(u) (+) states; and the B (3)Sigma(u) (-), A (1)Delta(u), and f (1)Sigma(u) (+) states are determined by the mixing of two or more electron configurations. The excitation spectrum of the oxygen molecule, calculated with the aforementioned exchange correlation functionals, shows that the results are quite sensitive to the choice of functional. The LDAxc and the B(88X)+P(86C) functionals produce similar spectroscopic patterns with a single strongly absorbing band positioned at 19.82 and 19.72 eV, respectively, while the asymptotically corrected exchange correlation functionals of the SAOP and the LB94 varieties yield similar excitation spectra where the computed strongly absorbing band is located at 16.09 and 16.42 eV, respectively. However, all of the exchange correlation functionals yield only one strongly absorbing band (oscillator strength greater than 0.1) in the energy interval of 0-20 eV, which is assigned to a X (3)Sigma(g) (-) to (3)Sigma(u) (-) transition. Furthermore, the oxygen molecule has a rich spectrum in the energy range of 14-20 eV and no spin allowed absorption bands are predicted to be observed in the range of 0-6 eV.     

Ab initio potential-energy surfaces of O2(X3Sigmag -, a1Deltag, b1Sigmag +) +O2 (X3Sigmag -, a1Deltag, b1Sigmag +): mechanism of quenching of O2 (a 1Deltag)

Ab initio computational studies were carried out in order to explore the possible mechanisms of quenching of O(2)(a (1)Delta(g)) by O(2)(X (3)Sigma(g) (-)): the self-quenching of O(2)(a (1)Delta(g)) and other energy-transfer processes involving two O(2) molecules. All eighteen states arising from two O(2) molecules in the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) states are considered. After scans at the state-averaged complete active space self-consistent field method to identify possible regions of crossing between states belonging to different asymptotes, complete active state second-order perturbation theory high-symmetry optimization and low-symmetry scans established that four different minima on the seams of crossing (MSXs), arising between the a+a manifold and the X+b manifold and responsible for self-quenching: O(2)(a (1)Delta(g))+O(2)(a (1)Delta(g))-->O(2)(X (3)Sigma(g) (-))+O(2)(b (1)Sigma(g) (+)), have coplanar C(2h) or C(2v) symmetries and are only 0.45 eV barrier relative to the a+a asymptote and energetically easily accessible. The rate constant for this process was estimated based on the Landau-Zener formalism. The MSXs for quenching of O(2)(a (1)Delta(g)) by the ground state O(2)(X (3)Sigma(g) (-)):O(2)(a (1)Delta(g))+O(2)(X (3)Sigma(g) (-))-->O(2)(X (3)Sigma(g) (-))+O(2)(X (3)Sigma(g) (-)) require higher energies and the process is not likely to be important.     

Ground and lowest-lying electronic states of CoN. A multiconfigurational study

The lowest-lying X1Sigma+, a3Phi, b3II, c5Delta, A1Phi, and B1II electronic states of CoN have been investigated at the ab initio MRCI and MS-CASPT2 levels, with extended atomic basis sets and inclusion of scalar relativistic effects. Among the singlet states, the A1Phi and B1II states have been described for the first time. Potential energy curves, excitation energies, spectroscopic constants, and bonding character for all states are reported. Comparison with other early transition-metal nitrides (ScN, TiN, VN, and CrN), isoelectronic (NiC) and isovalent (RhN and IrN) species has been made, besides analyzing the B1II <=> X1+ electronic transition in terms of Franck-Condon factors, Einstein coefficients, and radiative lifetimes. At both levels of theory, the following energetic order has been obtained: X1Sigma+, a3Phi, b3II, c5Delta, A1Phi, and B1II, with good agreement with experimental results. In contrast, previous DFT and MRCI calculations predicted the ground state to be the 5Delta state.     

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.     

Observation and calculation of the Cs2 2 3Delta(1g) and b 3Pi(0u) states

The Cs(2) 2 (3)Delta(1g) and b (3)Pi(0u) states have been observed by infrared-infrared double resonance spectroscopy for the first time. 221 2 (3)Delta(1g)<--A (1)Sigma(u) (+)<--X (1)Sigma(g) (+) double resonance lines have been assigned to transitions into the 2 (3)Delta(1g) v=6-13 vibrational levels. Resolved fluorescence into the b (3)Pi(0u) v(')=0-48 levels has been recorded. Molecular constants and potential energy curves are determined by the global fit of the entire set of the experimental data. Theoretical potential energy curves of the 2 (3)Delta(g) and b (3)Pi(u) states have been determined in the framework of the pseudopotential method and are compared with the experimental results.     

Ab initio study of the electronic structure and bonding of aluminum nitride

For the diatomic aluminum nitride (AlN), we have constructed potential energy curves for 45 states employing multi-reference variational methods and quantitative basis sets. Thirty-six states are relatively strongly bound, five present local minima, and four are of repulsive nature. Almost all states are of intense multi-reference character rendering their calculation and interpretation quite problematic. Our tentative assignment of the ground state is 3Pi, while a 3Sigma- state is above by less than 1 kcal/mol. Our best estimate for the binding energy of the X3Pi state is D0 = 56.0 +/- 0.5 kcal/mol at re = 1.783 A, in good agreement with the experimental values of D = 66 +/- 9 kcal/mol and re = 1.7864 A. The binding energy of the A3Sigma- state is very similar to the X state because they both correlate to the ground-state atoms, but the bond distance of the former is 0.13 A longer. The first seven states can be tagged as follows: X3Pi, A3Sigma-, a1Sigma+, b1Pi, c1Delta, B3Sigma+, and d1Sigma+, a rather definitive order with the exception of X and A states.     

Low-lying singlet and triplet electronic states of RhB

The low-lying XSigma+, a3Delta, A1Delta, b3Sigma+, B1Pi, c3Pi, C1Phi, D1Sigma+, E1Pi, d3Phi, and e3Pi electronic states of RhB have been investigated at the ab initio level, using the multistate multiconfigurational second-order perturbation (MS-CASPT2) theory, with extended atomic basis sets and inclusion of scalar relativistic effects. Among the eleven electronic states included in this work, only three (the X1Sigma+, D1Sigma+, and E1Pi states) have been investigated experimentally. Potential energy curves, spectroscopic constants, dipole moments, binding energies, and chemical bonding aspects are presented for all electronic states.     

Theoretical study of the gas-phase chemiionization reactions La + O and La + O2

The La + O and La + O 2 chemiionization reactions have been investigated with quantum chemical methods. For La + O 2(X (3)Sigma g) and La + O 2(a (1)Delta g), the chemiionization reaction La + O 2 --> LaO 2 (+) + e (-) has been shown to be endothermic and does not contribute to the experimental chemielectron spectra. For the La + O 2(X (3)Sigma g) reaction conditions, chemielectrons are produced by La + O 2 --> LaO + O, followed by La + O --> LaO (+) + e (-). This is supported by the same chemielectron band, arising from La + O --> LaO (+) + e (-), being observed from both the La + O( (3)P) and La + O 2(X (3)Sigma g) reaction conditions. For La + O 2(a (1)Delta g), a chemielectron band with higher electron kinetic energy than that obtained from La + O 2(X (3)Sigma g) is observed. This is attributed to production of O( (1)D) from the reaction La + O 2(a (1)Delta g) --> LaO + O( (1)D), followed by chemiionization via the reaction La + O( (1)D) --> LaO (+) + e (-). Potential energy curves are computed for a number of states of LaO, LaO* and LaO (+) to establish mechanisms for the observed La + O --> LaO (+) + e (-) chemiionization reactions.     

Spectroscopic properties and potential energy curves of low-lying electronic states of RuC

The RuC molecule has been a challenging species due to the open-shell nature of Ru resulting in a large number of low-lying electronic states. We have carried out state-of-the-art calculations using the complete active space multiconfiguration self-consistent field followed by multireference configuration interaction methods that included up to 18 million configurations, in conjunction with relativistic effects. We have computed 29 low-lying electronic states of RuC with different spin multiplicities and spatial symmetries with energy separations less than 38,000 cm(-1). We find two very closely low-lying electronic states for RuC, viz., 1Sigma+ and 3Delta with the 1Sigma+ being stabilized at higher levels of theory. Our computed spectroscopic constants and dipole moments are in good agreement with experiment although we have reported more electronic states than those that have been observed experimentally. Our computations reveal a strongly bound 1Sigma+ state with a large dipole moment which is most likely the experimentally observed ground state and an energetically close 3Delta state with a smaller dipole moment. Overall our computed spectroscopic constants of the excited states with energy separations less than 18,000 cm(-1) agree quite well with those of the corresponding observed states.     

Spin-orbit coupling in O2(v)+O2 collisions: a new energy transfer mechanism

A reduced dimensionality model is used to study the relaxation of highly vibrationally excited O(2)(X (3)Sigma(g) (-),v>/=20) in collisions with O(2)(X (3)Sigma(g) (-),v=0). Spin-orbit coupled potential energy surfaces are employed to incorporate the vibrational-to-electronic energy transfer mechanism involving the O(2)(a (1)Delta(g)) and O(2)(b (1)Sigma(g) (+)) excited states. The transition probabilities obtained show a sharp increase for v>/=26 providing the first direct evidence of the important role played by the electronic energy transfer processes in the depletion of O(2)(X (3)Sigma(g) (-),v>/=26).     

Coupled cluster investigation on the low-lying electronic states of CuCN and CuNC and the ground state barrier to isomerization

The observation of several metal cyanides and isocyanides in interstellar space has raised much interest these molecules. Optimum molecular structures, harmonic vibrational frequencies, and dipole moments of the ground electronic states (X1Sigma+), triplet excited states, and open shell singlet excited states of CuCN and CuNC were determined using different levels of nonrelativistic and scalar relativistic (Douglas-Kroll) [Ann. Phys. 82, 89 (1979)] coupled cluster theory in conjunction with atomic natural orbital basis sets and correlation consistent basis sets. For the relativistic computations the specially contracted correlation consistent Douglas-Kroll (DK) basis sets were used. Moreover, barriers to isomerization from CuCN to CuNC were computed. The predicted structures of the X1Sigma+ state for CuCN are re(Cu-C)=1.826 A and re(C-N)=1.167 A, at the most sophisticated level of theory, the scalar relativistic DK-CCSD(T)/cc-pVQZ(DK) method. These results are in excellent agreement with the experimentally determined Cu-C bond length of 1.829 A and C-N bond distance of 1.162 A. At the same level of theory, the zero-point corrected barrier to isomerization from CuCN to CuNC is estimated to be 14.7 kcal mol(-1), and the cyanide is more stable than the isocyanide by 11.5 kcal mol(-1). For both CuCN and CuNC the 3Sigma+ state is the lowest lying excited electronic state. At the DK-CCSD/cc-pVQZ(DK) level of theory, the energetic ordering of excited states of CuCN and CuNC is X1Sigma+相似文献   

Renner-Teller and spin-orbit vibronic coupling effects in linear triatomic molecules with a half-filled pi shell

The vibronic and spin-orbit-induced interactions among the (3)Sigma(-), (1)Delta, and (1)Sigma(+) electronic states arising from a half-filled pi orbital of a linear triatomic molecule are considered, employing the microscopic (Breit-Pauli) spin-orbit coupling operator. The 6 x 6 Hamiltonian matrix is derived in a diabatic spin-orbital electronic basis set, including terms up to fourth order in the expansion of the molecular Hamiltonian in the bending normal coordinate about the linear geometry. The symmetry properties of the Hamiltonian are analyzed. Aside from the nonrelativistic fourth-order Renner-Teller vibronic coupling within the (1)Delta state and the second-order nonrelativistic vibronic coupling between the (1)Sigma(+) and (1)Delta states, there exist zeroth-order, first-order, as well as third-order vibronic coupling terms of spin-orbit origin. The latter are absent when the phenomenological expression for the spin-orbit coupling operator is used instead of the microscopic form. The effects of the nonrelativistic and spin-orbit-induced vibronic coupling mechanisms on the (3)Sigma(-), (1)Delta, and (1)Sigma(+) adiabatic potential energy surfaces as well as on the spin-vibronic energy levels are discussed for selected parameter values.     

Dissociation of the OCS+ ion in low-lying electronic states studied using multiconfiguration second-order perturbation theory

Complete active space self-consistent-field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with atomic natural orbital basis sets were performed to investigate the S-loss direct dissociation of the 1 2Pi(X 2Pi), 2 2Pi(A 2Pi), 1 2Sigma+(B 2Sigma+), 1 4Sigma-, 1 2Sigma-, and 1 2Delta states of the OCS+ ion and the predissociations of the 1 2Pi, 2 2Pi, and 1 2Sigma+ states. Our calculations indicate that the S-loss dissociation products of the OCS(+) ion in the six states are the ground-state CO molecule plus the S+ ion in different electronic states. The CASPT2//CASSCF potential energy curves were calculated for the S-loss dissociation from the six states. The calculations indicate that the dissociation of the 1 4Sigma- state leads to the CO + S+ (4Su) products representing the first dissociation limit; the dissociations of the 1 2Pi, 1 2Sigma-, and 1 2Delta states lead to the CO + S+(2Du) products representing the second dissociation limit; and the dissociations of the 2 2Pi and 1 2Sigma+ states lead to the CO + S+(2Pu) products representing the third dissociation limit. Seams of the 1 2Pi-1 4Sigma-, 2 2Pi-1 4Sigma-, 2 2Pi-1 2Sigma-, 2 2Pi-1 2Delta, and 1 2Sigma(+)-1 4Sigma- potential energy surface intersections were calculated at the CASPT2 level, and the minima along the seams were located. The calculations indicate that within the experimental energy range (15.07-16.0 eV) the 2 2Pi(A 2Pi) state can be predissociated by 1 4Sigma- forming the S+(4Su) ion and can undergo internal conversion to 1 2Pi followed by the direct dissociation of 1 2Pi forming S+(2Du) and that within the experimental energy range (16.04-16.54 eV) the 1 2Sigma+(B 2Sigma+) state can be predissociated by 1 4Sigma- forming the S+(4Su) ion and can undergo internal conversion to 2 2Pi followed by the predissociation of 2 2Pi by 1 2Sigma- and 1 2Delta forming the S+(2Du) ion. These indications are in line with the experimental fact that both the 4Su and 2Du states of the S+ ion can be formed from the 2 2Pi and 1 2Sigma+ states of the OCS+ ion.     

Theoretical study on the electronic states of NaLi

Configuration interaction calculations have been carried out on electronic states of the NaLi molecule and the cation NaLi(+). Potential energy curves are presented for the lowest nine (1)Sigma(+), seven (1)Pi, four (1)Delta, eight (3)Sigma(+), seven (3)Pi, and four (3)Delta states of NaLi as well as for the lowest ten (2)Sigma(+), six (2)Pi, and two (2)Delta states of NaLi(+). The results of the present many-electron configuration interaction calculations on the cation are in support of previous core-polarization effective potential calculations. The present calculations on the NaLi molecule are complementary to previous theoretical work on this system, including recently observed electronic states that had not been calculated previously as well as an investigation of nonadiabatic effects leading to spectral perturbations. Furthermore, ab initio potential energy curves of the neutral and the ground state of the cation are employed to determine quantum defect that may be employed to generate potential energy curves for nd and (n+1)p (for n>3) Rydberg states of NaLi. The present results on the 3 (1)Pi and 4 (1)Pi states are in good agreement with recent experimental work, whereas on the basis of theoretical data, the recently observed state 5 (1)Pi is better described as 6 (1)Pi.     

Electronic spectra of linear isoelectronic clusters C2n+1S and C2n+1Cl+ (n = 0-4): an ab initio study

Structures and stabilities of linear carbon chains C2n+1S and C2n+1Cl+ (n=0-4) in their ground states have been investigated by the CCSD and B3LYP approaches. The CASSCF calculations have been used to determine geometries of selected excited states of both isoelectronic series. Linear C2n+1S cluster has a cumulenic carbon framework, whereas its isoelectronic C2n+1Cl+ has a dominant character of acetylenic structure in the vicinity of terminal Cl. The vertical excitation energies of low-lying excited states have been calculated by the CASPT2 method. Calculations show that the excitation energies have nonlinear size dependence. The 2(1)Sigma+<--X1Sigma+ transition energy in C2n+1S has a limit of 1.78 eV, as the chain size is long enough. The predicted vertical excitation energies for relatively strong 1(1)Pi<--X1Sigma+ and 2(1)Sigma+<--X1Sigma+ transitions are in reasonable agreement with available experimental values. The spin-orbit effect on the spin-forbidden transition in both series is generally small, and the enhancement of the spin-forbidden transition by spin-orbit coupling exhibits geometrical and electronic structural dependence.     

Charge transfer in the electron donor-acceptor complex BH3NH3

As a simple yet strongly binding electron donor-acceptor (EDA) complex, BH(3)NH(3) serves as a good example to study the electron pair donor-acceptor complexes. We employed both the ab initio valence bond (VB) and block-localized wave function (BLW) methods to explore the electron transfer from NH(3) to BH(3). Conventionally, EDA complexes have been described by two diabatic states: one neutral state and one ionic charge-transferred state. Ab initio VB self-consistent field (VBSCF) computations generate the energy profiles of the two diabatic states together with the adiabatic (ground) state. Our calculations evidently demonstrated that the electron transfer between NH(3) and BH(3) falls in the abnormal regime where the reorganization energy is less than the exoergicity of the reaction. The nature of the NH(3)-BH(3) interaction is probed by an energy decomposition scheme based on the BLW method. We found that the variation of the charge-transfer energy with the donor-acceptor distance is insensitive to the computation levels and basis sets, but the estimation of the amount of electron transferred heavily depends on the population analysis procedures. The recent resurgence of interest in the nature of the rotation barrier in ethane prompted us to analyze the conformational change of BH(3)NH(3), which is an isoelectronic system with ethane. We found that the preference of the staggered structure over the eclipsed structure of BH(3)NH(3) is dominated by the Pauli exchange repulsion.     

The highly anharmonic BH5 potential energy surface characterized in the ab initio limit

The strong sensitivity to level of theory of the salient features of the ground state potential energy surface of BH(5) has been overcome by rigorously converged ab initio computations employing correlation-consistent basis sets cc-p(C)VXZ (X=2-6), explicitly correlated R12 corrections, and coupled-cluster theory complete through quadruple excitations (CCSDTQ). Extrapolations via our focal point method yield a C(s)-symmetry global minimum of BH(3)-H(2) type featuring interfragment B-H distances of (1.401, 1.414) A, an H(2) bond length elongated to 0.803 A, and a BH(3)+H(2) dissociation energy D(e)(D(0))=6.6 (1.2) kcal mol(-1). The classical barriers for H(2) internal rotation and hydrogen scrambling are 0.07 and 5.81 kcal mol(-1), respectively, above the BH(5) minimum. Our thermochemical computations yield Delta(f)H(0) ( degrees )[BH(5)(g)]=-111.3+/-0.2 kcal mol(-1)+Delta(f)H(0) ( degrees )[B(g)], which is limited in accuracy only by persistent uncertainties in the enthalpy of formation of gaseous boron. As a first step in investigating the extremely anharmonic 12-dimensional vibrational dynamics of BH(5), a complete quartic force field has been computed at the all-electron cc-pCVQZ CCSD(T) level of theory. Previous matrix isolation infrared assignments of the nu(2) and nu(9) stretching modes of BH(5) compare favorably with our computed vibrational fundamentals, but the experimental assignment of the nu(6) bending mode of the BH(3) moiety is not supported by computed isotopic shifts.     

Electronic states and potential energy curves of molybdenum carbide and its ions

The potential energy curves and spectroscopic constants of the ground and 29 low-lying excited states of MoC with different spin and spatial symmetries within 48 000 cm(-1) have been investigated. We have used the complete active space multiconfiguration self-consistent field methodology, followed by multireference configuration interaction (MRCI) methods. Relativistic effects were considered with the aid of relativistic effective core potentials in conjunction with these methods. The results are in agreement with previous studies that determined the ground state as X (3)Sigma(-). At the MRCISD+Q level, the transition energies to the 1 (3)Delta and 4 (1)Delta states are 3430 and 8048 cm(-1), respectively, in fair agreement with the results obtained by DaBell et al. [J. Chem. Phy. 114, 2938 (2001)], namely, 4003 and 7834 cm(-1), respectively. The three band systems located at 18 611, 20 700, and 22 520 cm(-1) observed by Brugh et al. [J. Chem. Phy. 109, 7851 (1998)] were attributed to the excited 11 (3)Sigma(-), 14 (3)Pi, and 15 (1)Pi states respectively. At the MRCISD level, these states are 17 560, 20 836, and 20 952 cm(-1) above the ground state respectively. We have also identified a (3)Pi state lying 14 309 cm(-1) above the ground state. The ground states of the molecular ions are predicted to be (4)Sigma(-) and (2)Delta for MoC(-) and MoC(+), respectively.     

Equilibrium structure and torsional barrier of BH3NH3

Born-Oppenheimer equilibrium structures, r(e)(BO), of the electronic ground state of the borazane (BH3NH3) molecule of C3v point-group symmetry are computed ab initio using the CCSD(T) method with basis sets up to quintuple-zeta quality. Inclusion of the counterpoise correction and extrapolation of the structural parameters to the complete basis set limit yield a best estimate of r(e)(BO) of BH3NH3. The anharmonic force field of BH3NH3, computed at the CCSD(T) level of theory with a basis set of triple-zeta quality, allows the determination of semi-experimental equilibrium rotational constants, which in turn result in a semi-experimental equilibrium structure, r(e)(SE). The r(e)(BO) and r(e)(SE) structures are in excellent agreement, indicating the validity of the methods used for their determination. The empirical mass-dependent structure, r(m)(1), of BH3NH3 is also determined. Although it is inferior in quality to the previous two structures, it is much more accurate than the standard empirical r0 and r(s) structures reported earlier for BH3NH3. The semi-experimental r(e)(SE) as well as the empirical r(m)(1) structures determined are based on experimental ground-state rotational constants available from the literature for nine isotopologues of borazane. The effective barrier to the internal rotation of BH3NH3, a molecule isoelectronic with CH3CH3, has been computed ab initio, employing the focal-point analysis (FPA) approach, to be 699 +/- 11 cm(-1). This compares favorably with an empirical redetermination of the effective barrier based on the above r(e)(SE) structure, V3 = 718(17) cm(-1).