Stark coefficients for highly excited rovibrational states of H2O

Quantum beat spectroscopy is combined with triple-resonance vibrational overtone excitation to measure the Stark coefficients (SCs) of the water molecule for 28 rovibrational levels lying from 27,600 to 41,000 cm(-1). These data provide a stringent test for assessing the accuracy of the available potential energy surfaces (PESs) and dipole moment surfaces (DMSs) of this benchmark molecule in this energy region, which is inaccessible by direct absorption. SCs, calculated using the combination of a high accuracy, spectroscopically determined PES and a recent ab initio DMS, are within the 1% accuracy of available experimental data for levels below 25,000 cm(-1), and within 4.5% for coefficients associated with levels up to 35,000 cm(-1). However, the error in the computed coefficients is over 60% for the very high rovibrational states lying just below the lowest dissociation threshold, due, it seems, to lack of a high accuracy PES in this region. The comparative analysis suggests further steps, which may bring the theoretical predictions closer to the experimental accuracy.     

Microwave spectrum of the H2DO+ ion: inversion-rotation transitions and inversion splitting

Inversion-rotation spectral lines of the monodeuterated hydronium ion, H(2)DO(+), have been observed by a source-modulation spectrometer in the millimeter- to submillimeter-wave region. The ion was generated by a hollow-cathode discharge in a gas mixture of H(2)O and D(2)O. Nine inversion-rotation lines were measured precisely for the lowest pair of inversion doublets in the frequency region from 210 to 720 GHz. The measured lines were analyzed to derive rotational constants in the inversion-doublet states and inversion splitting. The inversion splitting in the ground state was determined to be 1,215,866(410) MHz, that is, 40.5569(137) cm(-1), where the numbers in parentheses give probable uncertainties estimated from the Jacobian matrix of the assumed centrifugal distortion constants of the inversion-doublet states. The determined inversion splitting is off by -0.58 cm(-1) from the predicted value of 41.14 cm(-1) by Rayamaki et al. using high-order coupled cluster ab initio calculations [J. Chem. Phys. 118, 10929 (2003)], and by 0.039 cm(-1) from the observed value of 40.518(10) cm(-1) by Dong and Nesbitt using high-resolution jet-cooled infrared spectroscopy [J. Chem. Phys. 125, 144311 (2006)] beyond the quoted uncertainty. The most astronomically important transition 0(00)(-)-1(0)(+) for the ortho species was measured at 673,257.024(31) MHz, which could be used as a radioastronomical probe investigating interstellar chemistry of deuterium fractionation in space.     

Ab initio Quantum-Mechanical Calculations of Molecular Structure and Conformations of 2,2-Dichloroethanal in the Ground and Excited Lowest Triplet States

Ab initio calculations were carried out to investigate the molecular structure of 2,2-dichloroethanal (DCE, CHCl₂CHO) in the ground (S ₀) and excited lowest triplet (₁) states. It is found that electronic excitation of DCE from the S ₀ to T ₁ state occurs with top rotations and a loss of planarity of the carbonyl fragments. Six minima corresponding to three pairs of enantiomers were found on the potential energy surface (PES) of the DCE molecule in the ₁ state. Based on the PES calculated (by the UHF and CASSCF methods in a 6-31G** basis) for DCE in the ₁ state, the one-dimensional torsional and inversion problems and the two-dimensional torsional-inversion problems are solved. A comparison of the results has revealed a relationship between the torsional and inversion motions.     

Microwave spectrum of the HD2O+ ion: inversion-rotation transitions and inversion splitting

Inversion-rotation spectral lines of the dideuterated hydronium ion, HD2O+, have been observed by a source-modulation millimeter- to submillimeter-wave spectrometer. The ion was generated by a hollow-cathode discharge in a gas mixture of D2O and H2O in a free-space cell. Ten inversion-rotation lines were measured precisely for the lowest pair of inversion doublets in the frequency region from 380 to 730 GHz. The observed lines include the most astronomically important transitions, 0(00) (-)-1(10)+ for the para species at 380 538.031(32) MHz and 1(01) (-)-1(11)+ for the ortho species at 728 420.189(34) MHz, which could be used as a radio astronomical probe investigating interstellar chemistry of deuterium fractionation. An analysis of the measured lines has yielded the rotational constants in the ground doublet states and the inversion splitting. The inversion splitting in the ground state was determined to be 808 866(34) MHz, that is, 26.980 87(113) cm(-1), where the numbers in parentheses give uncertainties estimated from the Jacobian matrix of the assumed centrifugal distortion constants. The determined inversion splitting is off by -0.51 cm(-1) from the predicted value of 27.49 cm(-1) by Rajamaki et al. using high-order coupled cluster ab initio calculation [J. Chem. Phys. 118, 10929 (2003)], and by -0.0510 cm(-1) from the observed value of 27.0318(72) cm(-1) by Dong et al. using high-resolution jet-cooled infrared spectroscopy [J. Chem. Phys. 122, 224301 (2005)] beyond the quoted uncertainty.     

Three-dimensional Potential Energy Surface and Bound States of the Ar2-Ne Complex

The first three-dimensional interaction potential energy surface (PES) of the Ar₂-Ne complex is developed using the single and double excitation coupled cluster theory with noniterative treatment of triple excitations CCSD(T). The aug-cc-pVQZ basis sets are employed for all atoms, including an additional (3s3p2d2f1g) set of midpoint bond functions. The calculated single point energies are fitted to an analytic two-dimensional potential model at each of seven fixed r_Ar2 values. The seven model potentials are then used to construct the three-dimensional PES by interpolating along (r—r_e) using a sixth-order polynomial. The PES is used in the following rovibrational energy levels calculations. The comparisons of theoretical transition frequencies and spectroscopic constants with the experimental results are given.     

Large amplitude out-of-plane vibrations of 1,3-benzodioxole in the S0 and S1 states: an analysis of fluorescence and excitation spectra by ab initio calculations

We report the analytical expressions of the two-dimensional potential energy surfaces (PES) spanned by the puckering and flapping vibrations in the S0 and S1 states of 1,3-benzodioxole (BDO). Both PES are obtained from S0 and S1 energies computed on a grid of 2500 molecular geometries at the CASPT2 level. Both the S0 and S1 PES are anharmonic, and the planar geometry corresponds to a barrier that separates two minima at nonplanar geometries along the puckering/flapping deformations. Eigenvalues and eigenvectors of the mixed puckering/flapping modes are calculated by the Meyer flexible model. Improved vibronic levels, in better agreement with the observed spectra, are obtained by suitably optimized CASPT2 surfaces. To assign the lower-energy (0-500 cm(-1)) portion of emission and absorption spectra, we evaluate the band intensities by estimating the Franck-Condon factors between the puckering/flapping eigenvectors of the S0 and S1 states. From these calculations, we obtain a satisfactory assignment of the ground state IR spectra and of the fluorescence excitation spectrum. Both assignments are supported by the analysis of the vibrational structures of several single vibronic level (SVL) fluorescence spectra. The successful interpretation of these spectra shows that the S0 and S1 PES that we derive for BDO are substantially correct. The barrier heights in the two states are similar: 125.7 and 190.4 cm(-1) in S0 and in S1, respectively. In S0, the barrier is associated essentially with the puckering motion. In S1, it involves to a considerable extent also the flapping coordinate, whose vibrational frequency is much lower in S1 than in S0. This fact introduces a substantial Duschinsky effect in the S0-S1 transitions of BDO.     

Intramolecular proton transfer in malonaldehyde: accurate multilayer multi-configurational time-dependent Hartree calculations

Full-dimensional (multilayer) multi-configurational time-dependent Hartree calculations studying the intramolecular proton transfer in malonaldehyde based on a recent potential energy surface (PES) [Wang et al., J. Chem. Phys. 128, 224314 (2008)] are presented. The most accurate calculations yield a ground state tunneling splitting of 23.8 cm(-1) and a zero point energy of 14,678 cm(-1). Extensive convergence tests indicate an error margin of the quantum dynamics calculations for the tunneling splitting of about 0.2 cm(-1). These results are to be compared with the experimental value of the tunneling splitting of 21.58 cm(-1) and results of Monte Carlo calculations of Wang et al. on the same PES which yielded a zero point energy of 14,677.9 cm(-1) with statistical errors of 2-3 cm(-1) and a tunneling splitting of 21.6 cm(-1). The present data includes contributions resulting from the vibrational angular momenta to the tunneling splitting and the zero point energy of 0.2 cm(-1) and 2.4 cm(-1), respectively, which have been computed using a perturbative approach.     

Highly accurate potential-energy and dipole moment surfaces for vibrational state calculations of methane

Full-dimensional ab initio potential-energy surface (PES) and dipole moment surface are constructed for a methane molecule at the CCSD(T)/cc-pVTZ and MP2/cc-pVTZ levels of theory, respectively, by the modified Shepard interpolation method based on the fourth-order Taylor expansion [MSI(4th)]. The reference points for the interpolation have been set in the coupling region of CH symmetric and antisymmetric stretching modes so as to reproduce the vibrational energy levels related to CH stretching vibrations. The vibrational configuration-interaction calculations have been performed to obtain the energy levels and the absorption intensities up to 9000 cm(-1) with the use of MSI(4th)-PES. The calculated fundamental frequencies and low-lying vibrational energy levels show that MSI(4th) is superior to the widely employed quartic force field, giving a better agreement with the experimental values. The absorption bands of overtones as well as combination bands, which are caused by purely anharmonic effects, have been obtained up to 9000 cm(-1). Strongly coupled states with visible intensity have been found in the 6500-9000 cm(-1) region where the experimental data are still lacking.     

Measurements of line strengths in the 2nu1 band of the HO2 radical using laser photolysis/continuous wave cavity ring-down spectroscopy (cw-CRDS)

Absolute absorption cross sections of the absorption spectrum of the 2nu1 band of the HO2 radical in the near-IR region were measured by continuous wave cavity ring-down spectroscopy (cw-CRDS) coupled to laser photolysis in the wavelength range 6604-6696 cm(-1) with a resolution better than 0.003 cm(-1). Absolute absorption cross sections were obtained by measuring the decay of the HO2 self-reaction, and they are given for the 100 most intense lines. The most important absorption feature in this wavelength range was found at 6638.20 cm(-1), exhibiting an absorption cross section of sigma = 2.72 x 10(-19) cm2 at 50 Torr He. Using this absorption line, we obtain a detection limit for the HO2 radical at 50 Torr of 6.5 x 10(10) cm(-3).     

Variational calculation of specific, highly excited vibrational states in DFCO: comparison with experimental data

A stimulated emission pumping spectra of jet-cooled DFCO performed by Crane et al. (J. Mol. Spectrosc. 1997, 183, 273) has provided a great number of ro-vibrational lines up to 9000 cm(-1) of excitation energy. By combining a Jacobi-Wilson (JW) approach with a Davidson scheme, we calculate the lines provided by the experiment up to 9000 cm(-1) using an ab initio global potential energy surface (PES) developed by Kato et al. (J. Chem. Phys. 1997, 107, 6114). Comparisons between experimental and calculated data provide a critical test of the quality of the PES used. We show that the variational calculated energies can be efficiently corrected by taking into account the error observed for the A' fundamental transitions nu(i) (i = 1, ..., 5) and the first overtone 2nu(6). A detailed analysis of the eigenstates obtained by the calculation allows one to quantify the coupling between the different modes. Such an information is essential to understand and predict the energy flow through a DFCO molecule that is initially excited.     

The Role of High Excitations in Constructing Sub-spectroscopic Accuracy Intermolecular Potential of He-HCN: Critically Examined by the High-Resolution Spectra with Resonance States

Interpreting high-resolution rovibrational spectra of weakly bound complexes commonly requires spectroscopic accuracy (<1 cm^-1) potential energy surfaces (PES). Constructing high-accuracy ab initio PES relies on the high-level electronic structure approaches and the accurate physical models to represent the potentials. The coupled cluster approaches including single and double excitations with a perturbational estimate of triple excitations (CCSD(T)) have been termed the "gold standard" of electronic structure theory, and widely used in generating intermolecular interaction energies for most van der Waals complexes. However, for HCN-He complex, the observed millimeter-wave spectroscopy with high-excited resonance states has not been assigned and interpreted even on the ab initio PES computed at CCSD(T) level of theory with the complete basis set (CBS) limit. In this work, an effective three-dimensional ab initio PES for HCN-He, which explicitly incorporates dependence on the Q₁ (C-H) normal-mode coordinate of the HCN monomer has been calculated at the CCSD(T)/CBS level. The post-CCSD(T) interaction energy has been examined and included in our PES. Analytic two-dimensional PESs are obtained by least-squares fitting vibrationally averaged interaction energies for v₁(C-H)=0, and 1 to the Morse/Long-Range potential function form with root-mean-square deviations (RMSD) smaller than 0.011 cm^-1. The role and significance of the post-CCSD(T) interaction energy contribution are clearly illustrated by comparison with the predicted rovibrational energy levels. With or without post-CCSD(T) corrections, the value of dissociation limit (D₀) is 8.919 or 9.403 cm^-1, respectively. The predicted millimeter-wave transitions and intensities from the PES with post-CCSD(T) excitation corrections are in good agreement with the available experimental data with RMS discrepancy of 0.072 cm^-1. Moreover, the infrared spectrum for HCN-He complex is predicted for the first time. These results will serve as a good starting point and provide reliable guidance for future infrared studies of HCN doped in (He)_n clusters.     

Ab initio rotation-vibration spectra of HCN and HNC

We have calculated an ab initio HCN/HNC linelist for all transitions up to J= 25 and 18000 cm(-1) above the zero point energy. This linelist contains more than 200 million lines each with frequencies and transition dipoles. The linelist has been calculated using our semi-global HCN/HNC VQZANO + PES and dipole moment surface, which were reported in van Mourik et al. (J. Chem. Phys. 115 (2001) 3706). With this linelist we synthesise absorption spectra of HCN and HNC at 298 K and we present the band centre and band transition dipoles for the bands which are major features in these spectra. Several of the HCN bands and many of the HNC bands have not been previously studied. Our line intensities reproduce via fully ab initio methods the unusual intensity structure of the HCN CN stretch fundamental (00(0)1) for the first time and also the forbidden (02(2)0) HCN bending overtone. We also compare the J = 1-->0 pure rotational transition dipole in the HCN/HNC ground and vibrationally excited states with experimental and existing ab initio results.     

Continuous-wave cavity ringdown spectroscopy of the 8nu polyad of water in the 25,195-25,340 cm(-1) range

State-of-the-art experiments and calculations are used to record and assign the data obtained in the weakly absorbing blue energy region of the H2O spectrum. Continuous-wave cavity ringdown absorption spectroscopy with Doppler resolution is used to probe the range from 25,195 to 25,470 cm(-1) with an absorption sensitivity of approximately 1 parts per 10(9) (ppb)/cm. 62 lines of the polyad nu(OH)=8 are reported, of which 43 are assigned using variational nuclear calculations. The study includes absorption line intensities (in the range of 10(-28)-10(-26) cmmolecule) for all lines and self-broadening pressure coefficient for a few lines. The newly obtained energy levels are also reported.     

Jet cooled spectroscopy of H2DO+: Barrier heights and isotope-dependent tunneling dynamics from H3O+ to D3O+

The first high resolution spectroscopic data for jet cooled H2DO+ are reported, specifically via infrared laser direct absorption in the OH stretching region with a slit supersonic jet discharge source. Transitions sampling upper (0-) and lower (0+) tunneling states for both symmetric (nu1+ <-- 0+, nu1- <-- 0-, and nu1- <-- 0+) and antisymmetric (nu3+ <-- 0+ and nu3- <-- 0-) OH stretching bands are observed, where +/- refers to wave function reflection symmetry with respect to the planar umbrella mode transition state. The spectra can be well fitted to a Watson asymmetric top Hamiltonian, revealing band origins and rotational constants for benchmark comparison with high-level ab initio theory. Of particular importance are detection and assignment of the relatively weak band (nu1- <-- 0+) that crosses the inversion tunneling gap, which is optically forbidden in H3O+ or D3O+, but weakly allowed in H2DO+ by lowering of the tunneling transition state symmetry from D(3h) to C(2v). In conjunction with other H2DO+ bands, this permits determination of the tunneling splittings to within spectroscopic precision for each of the ground [40.518(10) cm(-1)], nu1 = 1 [32.666(6) cm(-1)], and nu3 = 1 [25.399(11) cm(-1)] states. A one-dimensional zero-point energy corrected potential along the tunneling coordinate is constructed from high-level ab initio CCSD(T) calculations (AVnZ, n = 3,4,5) and extrapolated to the complete basis set limit to extract tunneling splittings via a vibrationally adiabatic treatment. Perturbative scaling of the potential to match splittings for all four isotopomers permits an experimental estimate of DeltaV0 = 652.9(6) cm(-1) for the tunneling barrier, in good agreement with full six-dimensional ab initio results of Rajamaki, Miani, and Halonen (RMH) [J. Chem. Phys. 118, 10929 (2003)]. (DeltaV0 (RMH) = 650 cm(-1)). The 30%-50% decrease in tunneling splitting observed upon nu1 and nu3 vibrational excitations arises from an increase in OH stretch frequencies at the planar transition state, highlighting the transition between sp2 and sp3 hybridizations of the OHD bonds as a function of inversion bending angle.     

Kinetics of the CN + CS2 and CN + SO2 reactions

Diode infrared laser absorption spectroscopy was used to measure the rate constant (k(1)) of the CN + CS(2) reaction for the first time. k(1) was determined to be substantially pressure dependent with a value k(1) = (7.1 ± 0.2 to 41.9 ± 2.9) × 10(-12) cm(3) molecule(-1) s(-1) over 2-40 Torr at 298 K. The potential energy surface (PES) of the reaction was calculated using an ab initio method at B3LYP/6-311++G(d, p)//CCSD(T)/6-311++G(d, p) level of theory. Both experimental and computational results suggest that collision stabilization of the adduct NCSCS may dominate the reaction. The rate constant of the CN + SO(2) reaction was measured to be very slow with an upper limit of k(2) ≤ 3.1 × 10(-14) cm(3) molecule(-1) s(-1), in disagreement with an earlier reported measurement. The PES of this reaction reveals an entrance barrier against formation of the low energy adduct NCOSO, in agreement with the experimental result.     

Rovibrational spectra of ammonia. I. Unprecedented accuracy of a potential energy surface used with nonadiabatic corrections

In this work, we build upon our previous work on the theoretical spectroscopy of ammonia, NH(3). Compared to our 2008 study, we include more physics in our rovibrational calculations and more experimental data in the refinement procedure, and these enable us to produce a potential energy surface (PES) of unprecedented accuracy. We call this the HSL-2 PES. The additional physics we include is a second-order correction for the breakdown of the Born-Oppenheimer approximation, and we find it to be critical for improved results. By including experimental data for higher rotational levels in the refinement procedure, we were able to greatly reduce our systematic errors for the rotational dependence of our predictions. These additions together lead to a significantly improved total angular momentum (J) dependence in our computed rovibrational energies. The root-mean-square error between our predictions using the HSL-2 PES and the reliable energy levels from the HITRAN database for J = 0-6 and J = 7∕8 for (14)NH(3) is only 0.015 cm(-1) and 0.020∕0.023 cm(-1), respectively. The root-mean-square errors for the characteristic inversion splittings are approximately 1∕3 smaller than those for energy levels. The root-mean-square error for the 6002 J = 0-8 transition energies is 0.020 cm(-1). Overall, for J = 0-8, the spectroscopic data computed with HSL-2 is roughly an order of magnitude more accurate relative to our previous best ammonia PES (denoted HSL-1). These impressive numbers are eclipsed only by the root-mean-square error between our predictions for purely rotational transition energies of (15)NH(3) and the highly accurate Cologne database (CDMS): 0.00034 cm(-1) (10 MHz), in other words, 2 orders of magnitude smaller. In addition, we identify a deficiency in the (15)NH(3) energy levels determined from a model of the experimental data.     

The absorption spectrum of water near 750 nm by CW-CRDS: contribution to the search of water dimer absorption

The absorption spectrum of natural water vapour around 750 nm has been recorded with a typical sensitivity of 3 x 10(-10) cm(-1) using a cw cavity ring down spectroscopy set up based on a Ti:sapphire laser. The 13 312.4-13 377.7 cm(-1) spectral interval was chosen as it corresponds to the region where water dimer absorption was recently measured (K. Pfeisticker et al., Science, 2003, 300, 2078-2080). The line parameters (wavenumber and intensity) of a total of 286 lines of water vapor were measured by a one by one fit of the lines to a Voigt profile. For the main water isotopologue, 276 lines were measured with line intensities as weak as 5 x 10(-29) cm molecule(-1)i.e. about 50 times smaller than the weakest H(2)16O line intensities included in the 2004 edition of the HITRAN database. On the basis of the predictions of Schwenke and Partridge, all but 16 lines could be assigned to different isotopologues of water (H(2)16O, H(2)18O, and HD16O) present in natural abundance in the sample. A total of 272 energy levels of H(2)16O were determined and rovibrationally assigned to 18 upper vibrational states. Half of them had not been reported previously. The importance of the additional absorbance resulting from the observation of many new weak lines is discussed in relation to the detection of water dimer absorption and compared to the absorbance predicted by Schwenke and Partridge. The quality of the line parameters of water monomer is shown to be of crucial importance to identify the absorbance of the water dimer in the considered region.     

Vibrations in the B4 rhombic structure

A double minimum six-dimensional potential energy surface (PES) is determined in symmetry coordinates for the most stable rhombic (D2h) B4 isomer in its 1Ag electronic ground state by fitting to energies calculated ab initio. The PES exhibits a barrier to the D4h square structure of 255 cm(-1). The vibrational levels (J=0) are calculated variationally using an approach which involves the Watson kinetic energy operator expressed in normal coordinates. The pattern of about 65 vibrational levels up to 1600 cm(-1) for all stable isotopomers is analyzed. Analogous to the inversion in ammonia-like molecules, the rhombus rearrangements lead to splittings of the vibrational levels. In B4 it is the B1g (D4h) mode which distorts the square molecule to its planar rhombic form. The anharmonic fundamental vibrational transitions of 11B4 are calculated to be (splittings in parentheses): G(0)=2352(22) cm(-1), nu1(A1g)=1136(24) cm(-1), nu2(B1g)=209(144) cm(-1), nu3(B2g)=1198(19) cm(-1), nu4(B2u)=271(24) cm(-1), and nu5(Eu)=1030(166) cm(-1) (D4h notation). Their variations in all stable isotopomers were investigated. Due to the presence of strong anharmonic resonances between the B1g in-plane distortion and the B2u out-of-plane bending modes, the higher overtones and combination levels are difficult to assign unequivocally.     

Ultraviolet absorption and vibrational spectra of 2-fluoro-5-bromopyridine

The ultraviolet absorption spectrum in the range 340-185 nm in the vapour and solution phase has been measured for 2-fluoro-5-bromopyridine. Three fairly intense band systems identified as the pi* <-- pi transitions II, III and IV have been observed. A detailed vibronic analysis of the vapor and solution spectra is presented. The first system of bands is resolved into about sixty-two distinct vibronic bands in the vapour-phase spectrum. The 0,0 band is located at 35944 cm(-1). Two well-developed progressions, in which the excited state frequencies nu'25 (283 cm(-1)) and nu'19 (550 cm(-1)) are excited by several quanta, have been observed. The corresponding excited state vibrational and anharmonicity constants are found to be omega'i = 292 cm(-1), x'ii = 4.5 cm(-1) (i = 25) and omega'i = 563.8 cm(-1), x'ii = 6.9 cm(-1) (i = 19). The other two band systems show no vibronic structure, the band maxima being located at 48346 and 52701 cm(-1), respectively. The oscillator strength of the band systems in different solutions and the excited state dipole moments associated with the first two transitions have been determined by the solvent-shift method. The infrared spectrum in the region 4000-130 cm(-1) and the laser Raman spectrum of the molecule in the liquid state have been measured and a complete vibrational assignment of the observed frequencies is given. A correlation of the ground and excited state fundamental frequencies observed in the UV absorption spectrum with the Raman or infrared frequencies is presented.     

Calibration-quality adiabatic potential energy surfaces for H3(+) and its isotopologues

Calibration-quality ab initio adiabatic potential energy surfaces (PES) have been determined for all isotopologues of the molecular ion H(3)(+). The underlying Born-Oppenheimer electronic structure computations used optimized explicitly correlated shifted Gaussian functions. The surfaces include diagonal Born-Oppenheimer corrections computed from the accurate electronic wave functions. A fit to the 41,655 ab initio points is presented which gives a standard deviation better than 0.1 cm(-1) when restricted to the points up to 6000 cm(-1) above the first dissociation asymptote. Nuclear motion calculations utilizing this PES, called GLH3P, and an exact kinetic energy operator given in orthogonal internal coordinates are presented. The ro-vibrational transition frequencies for H(3)(+), H(2)D(+), and HD(2)(+) are compared with high resolution measurements. The most sophisticated and complete procedure employed to compute ro-vibrational energy levels, which makes explicit allowance for the inclusion of non-adiabatic effects, reproduces all the known ro-vibrational levels of the H(3)(+) isotopologues considered to better than 0.2 cm(-1). This represents a significant (order-of-magnitude) improvement compared to previous studies of transitions in the visible. Careful treatment of linear geometries is important for high frequency transitions and leads to new assignments for some of the previously observed lines. Prospects for further investigations of non-adiabatic effects in the H(3)(+) isotopologues are discussed. In short, the paper presents (a) an extremely accurate global potential energy surface of H(3)(+) resulting from high accuracy ab initio computations and global fit, (b) very accurate nuclear motion calculations of all available experimental line data up to 16,000 cm(-1), and (c) results suggest that we can predict accurately the lines of H(3)(+) towards dissociation and thus facilitate their experimental observation.