Communication: rigorous calculation of dissociation energies (D0) of the water trimer, (H2O)3 and (D2O)3

Using a recent, full-dimensional, ab initio potential energy surface [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] together with rigorous diffusion Monte Carlo calculations of the zero-point energy of the water trimer, we report dissociation energies, D(0), to form one monomer plus the water dimer and three monomers. The calculations make use of essentially exact zero-point energies for the water trimer, dimer, and monomer, and benchmark values of the electronic dissociation energies, D(e), of the water trimer [J. A. Anderson, K. Crager, L. Fedoroff, and G. S. Tschumper, J. Chem. Phys. 121, 11023 (2004)]. The D(0) results are 3855 and 2726 cm(-1) for the 3H(2)O and H(2)O + (H(2)O)(2) dissociation channels, respectively, and 4206 and 2947 cm(-1) for 3D(2)O and D(2)O + (D(2)O)(2) dissociation channels, respectively. The results have estimated uncertainties of 20 and 30 cm(-1) for the monomer plus dimer and three monomer of dissociation channels, respectively.     

Vibrational spectroscopy of (SO4(2-)).(H2O)n clusters, n=1-5: harmonic and anharmonic calculations and experiment

The vibrational spectroscopy of (SO4(2-)).(H2O)n is studied by theoretical calculations for n=1-5, and the results are compared with experiments for n=3-5. The calculations use both ab initio MP2 and DFT/B3LYP potential energy surfaces. Both harmonic and anharmonic calculations are reported, the latter with the CC-VSCF method. The main findings are the following: (1) With one exception (H2O bending mode), the anharmonicity of the observed transitions, all in the experimental window of 540-1850 cm(-1), is negligible. The computed anharmonic coupling suggests that intramolecular vibrational redistribution does not play any role for the observed linewidths. (2) Comparison with experiment at the harmonic level of computed fundamental frequencies indicates that MP2 is significantly more accurate than DFT/B3LYP for these systems. (3) Strong anharmonic effects are, however, calculated for numerous transitions of these systems, which are outside the present observation window. These include fundamentals as well as combination modes. (4) Combination modes for the n=1 and n=2 clusters are computed. Several relatively strong combination transitions are predicted. These show strong anharmonic effects. (5) An interesting effect of the zero point energy (ZPE) on structure is found for (SO4(2-)).(H2O)(5): The global minimum of the potential energy corresponds to a C(s) structure, but with incorporation of ZPE the lowest energy structure is C2v, in accordance with experiment. (6) No stable structures were found for (OH-).(HSO4-).(H2O)n, for n相似文献   

Communication: quasiclassical trajectory calculations of correlated product-state distributions for the dissociation of (H2O)2 and (D2O)2

Stimulated by recent experiments [B. E. Rocher-Casterline, L. C. Ch'ng, A. K. Mollner, and H. Reisler, J. Chem. Phys. 134, 211101 (2011)], we report quasiclassical trajectory calculations of the dissociation dynamics of the water dimer, (H(2)O)(2) (and also (D(2)O)(2)) using a full-dimensional ab initio potential energy surface. The dissociation is initiated by exciting the H-bonded OH(OD)-stretch, as done experimentally for (H(2)O)(2). Normal mode analysis of the fragment pairs is done and the correlated vibrational populations are obtained by (a) standard histogram binning (HB), (b) harmonic normal-mode energy-based Gaussian binning (GB), and (c) a modified version of (b) using accurate vibrational energies obtained in the Cartesian space. We show that HB allows opening quantum mechanically closed states, whereas GB, especially via (c), gives physically correct results. Dissociation of both (H(2)O)(2) and (D(2)O)(2) mainly produces either fragment in the bending excited (010) state. The H(2)O(J) and D(2)O(J) rotational distributions are similar, peaking at J = 3-5. The computations do not show significant difference between the ro-vibrational distributions of the donor and acceptor fragments. Diffusion Monte Carlo computations are performed for (D(2)O)(2) providing an accurate zero-point energy of 7247 cm(-1), and thus, a benchmark D(0) of 1244 ± 5 cm(-1).     

Ab initio potential energy and dipole moment surfaces of (H2O)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with M?ller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).     

An isotopic-independent highly accurate potential energy surface for CO2 isotopologues and an initial (12)C(16)O2 infrared line list

An isotopic-independent, highly accurate potential energy surface (PES) has been determined for CO(2) by refining a purely ab initio PES with selected, purely experimentally determined rovibrational energy levels. The purely ab initio PES is denoted Ames-0, while the refined PES is denoted Ames-1. Detailed tests are performed to demonstrate the spectroscopic accuracy of the Ames-1 PES. It is shown that Ames-1 yields σ(rms) (root-mean-squares error) = 0.0156 cm(-1) for 6873 J = 0-117 (12)C(16)O(2) experimental energy levels, even though less than 500 (12)C(16)O(2) energy levels were included in the refinement procedure. It is also demonstrated that, without any additional refinement, Ames-1 yields very good agreement for isotopologues. Specifically, for the (12)C(16)O(2) and (13)C(16)O(2) isotopologues, spectroscopic constants G(v) computed from Ames-1 are within ±0.01 and 0.02 cm(-1) of reliable experimentally derived values, while for the (16)O(12)C(18)O, (16)O(12)C(17)O, (16)O(13)C(18)O, (16)O(13)C(17)O, (12)C(18)O(2), (17)O(12)C(18)O, (12)C(17)O(2), (13)C(18)O(2), (13)C(17)O(2), (17)O(13)C(18)O, and (14)C(16)O(2) isotopologues, the differences are between ±0.10 and 0.15 cm(-1). To our knowledge, this is the first time a polyatomic PES has been refined using such high J values, and this has led to new challenges in the refinement procedure. An initial high quality, purely ab initio dipole moment surface (DMS) is constructed and used to generate a 296 K line list. For most bands, experimental IR intensities are well reproduced for (12)C(16)O(2) using Ames-1 and the DMS. For more than 80% of the bands, the experimental intensities are reproduced with σ(rms)(ΔI) < 20% or σ(rms)(ΔI∕δ(obs)) < 5. A few exceptions are analyzed and discussed. Directions for future improvements are discussed, though it is concluded that the current Ames-1 and the DMS should be useful in analyzing and assigning high-resolution laboratory or astronomical spectra.     

Ab initio studies on Al(+)(H(2)O)(n), HAlOH(+)(H(2)O)(n-1), and the size-dependent H(2) elimination reaction

We report computational studies on Al(+)(H(2)O)(n), and HAlOH(+)(H(2)O)(n-1), n = 6-14, by the density functional theory based ab initio molecular dynamics method, employing a planewave basis set with pseudopotentials, and also by conventional methods with Gaussian basis sets. The mechanism for the intracluster H(2) elimination reaction is explored. First, a new size-dependent insertion reaction for the transformation of Al(+)(H(2)O)(n), into HAlOH(+)(H(2)O)(n-1) is discovered for n > or = 8. This is because of the presence of a fairly stable six-water-ring structure in Al(+)(H(2)O)(n) with 12 members, including the Al(+). This structure promotes acidic dissociation and, for n > or = 8, leads to the insertion reaction. Gaussian based BPW91 and MP2 calculations with 6-31G* and 6-31G** basis sets confirmed the existence of such structures and located the transition structures for the insertion reaction. The calculated transition barrier is 10.0 kcal/mol for n = 9 and 7.1 kcal/mol for n = 8 at the MP2/6-31G** level, with zero-point energy corrections. Second, the experimentally observed size-dependent H(2) elimination reaction is related to the conformation of HAlOH(+)(H(2)O)(n-1), instead of Al(+)(H(2)O)(n). As n increases from 6 to 14, the structure of the HAlOH(+)(H(2)O)(n-1) cluster changes into a caged structure, with the Al-H bond buried inside, and protons produced in acidic dissociation could then travel through the H(2)O network to the vicinity of the Al-H bond and react with the hydride H to produce H(2). The structural transformation is completed at n = 13, coincident approximately with the onset of the H(2) elimination reaction. From constrained ab initio MD simulations, we estimated the free energy barrier for the H(2) elimination reaction to be 0.7 eV (16 kcal/mol) at n = 13, 1.5 eV (35 kcal/mol) at n = 12, and 4.5 eV (100 kcal/mol) at n = 8. The existence of transition structures for the H(2) elimination has also been verified by ab initio calculations at the MP2/6-31G** level. Finally, the switch-off of the H(2) elimination for n > 24 is explored and attributed to the diffusion of protons through enlarged hydrogen bonded H(2)O networks, which reduces the probability of finding a proton near the Al-H bond.     

Low-lying isomers and finite temperature behavior of (H2O)6 -

(H2O)(6) (-) appears as a "magic" number water cluster in (H2O)(n) (-) mass spectra. The structure of the (H2O)(6) (-) isomer dominating the experimental population has been established only recently [N. I. Hammer et al., J. Phys. Chem. A 109, 7896 (2005)], and the most noteworthy characteristic of this isomer is the localization of the excess electron in the vicinity of a double-acceptor monomer. In the present work, we use a quantum Drude model to characterize the low-energy isomers and the finite temperature properties of (H2O)(6) (-). Comparison with ab initio calculations shows that the use of a water model employing distributed polarizabilities and distributed repulsive sites is necessary to correctly reproduce the energy ordering of the low-lying isomers. Both the simulations and the ab initio calculations predict that there are several isomers of (H2O)(6) (-) significantly lower in energy than the experimentally observed species, suggesting that the experimental distribution is far from equilibrium.     

Photodissociation spectroscopy of Zn+(H2O) and Zn+(D2O)

We report on a study of the photodissociation spectroscopy of weakly bound Zn+(H2O) and Zn+(D2O) complexes. The work is supported by ab initio electronic structure calculations of the ground and low-lying excited energy surfaces. We assign two molecular absorption bands in the near UV correlating to Zn+ (4s-4p)-based transitions, and identify vibrational progressions associated with both intermolecular and intramolecular vibrational modes of the cluster. Partially resolved rotational structure is consistent with a C(2V) equilibrium complex geometry. Experimental spectroscopic constants are in very good agreement with ab initio theoretical predictions. Results are compared with previous work on main group and transition metal ion-H2O clusters.     

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.     

H2O broadening of a CO2 line and its nearest neighbors near 6360 cm(-1)

Remote sensing of CO(2) requires high-fidelity reference data of spectral line parameters to be successful. The 6360 cm(-1) region is commonly used by satellites, field campaigns, and point-source gas sensors because it contains well-characterized and relatively isolated transitions of appropriate line strengths for atmospheric applications. However, the presence of gases other than CO(2), N(2), and O(2) can be a source of uncertainty for atmospheric measurements. Near 6360 cm(-1), there are numerous H(2)O and HDO transitions. Water makes up approximately 1-4% of Earth's lower atmosphere and can interfere with remote sensing measurements by (1) appearing as a direct spectral interference or (2) acting as a foreign broadener for CO(2) lines. The primary goal of this work was to quantify H(2)O broadening of CO(2) through precision spectroscopy measurements on the R16e transition at 6359.967 cm(-1) and its two nearest neighbors. A secondary goal was to assess the accuracy of H(2)O reference line parameters in the HITRAN 2008 database for spectrally removing typical levels of moisture from air samples containing approximately 400 ppm of CO(2).     

The mechanism of proton exchange: guided ion beam studies of the reactions, H(H2O)n+ (n=1-4)+D2O and D(D2O)n+ (n=1-4)+H2O

Reactions of protonated water clusters, H(H(2)O)(n) (+) (n=1-4) with D(2)O and their "mirror" reactions, D(D(2)O)(n) (+) (n=1-4) with H(2)O, are studied using guided-ion beam mass spectrometry. Absolute reaction cross sections are determined as a function of collision energy from thermal energy to over 10 eV. At low collision energies, we observe reactions in which H(2)O and D(2)O molecules are interchanged and reactions where H-D exchange has occurred. As the collision energy is increased, the H-D exchange products decrease and the water exchange products become dominant. At high collision energies, processes in which one or more water molecules are lost from the reactant ions become important, with simple collision-induced dissociation processes, i.e., those without H-D exchange, being dominant. Threshold energies of endothermic channels are measured and used to determine binding energies of the proton bound complexes, which are consistent with those determined by thermal equilibrium measurements and previous collision-induced dissociation studies. A kinetic scheme that relies only on the ratio of isomerization and dissociation rate constants successfully accounts for the kinetic energy dependence observed in the branching ratios for H-D and water exchange products in all systems. Rice-Ramsperger-Kassel-Marcus theory and ab initio calculations confirm the feasibility and establish the details of this kinetic model.     

Three dimensional quantum dynamics of (H-, H2) and its isotopic variants

We present the results of a time-dependent quantum mechanical investigation using centrifugal sudden approximation in the form of reaction probability as a function of collision energy (E(trans)) in the range 0.3-3.0 eV for a range of total angular momentum (J) values and the excitation function sigma(E(trans)) for the exchange reaction H(-) + H(2) (v = 0, j = 0) --> H(2) + H(-) and its isotopic variants in three dimensions on an accurate ab initio potential energy surface published recently (J. Chem. Phys. 2004, 121, 9343). The excitation function results are shown to be in excellent agreement with those obtained from crossed beam measurements by Zimmer and Linder for H(-) + D(2) collisions for energies below the threshold for electron detachment channel and somewhat larger than the most recent results of Haufler et al. for (H(-), D(2)) and (D(-), H(2)) collisions.     

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.     

The relaxation of OH (v=1) and OD (v=1) by H2O and D2O at temperatures from 251 to 390 K

We report rate coefficients for the relaxation of OH(v=1) and OD(v=1) by H2O and D2O as a function of temperature between 251 and 390 K. All four rate coefficients exhibit a negative dependence on temperature. In Arrhenius form, the rate coefficients for relaxation (in units of 10(-12) cm3 molecule-1 s-1) can be expressed as: for OH(v=1)+H2O between 263 and 390 K: k=(2.4+/-0.9) exp((460+/-115)/T); for OH(v=1)+D2O between 256 and 371 K: k=(0.49+/-0.16) exp((610+/-90)/T); for OD(v=1)+H2O between 251 and 371 K: k=(0.92+/-0.16) exp((485+/-48)/T); for OD(v=1)+D2O between 253 and 366 K: k=(2.57+/-0.09) exp((342+/-10)/T). Rate coefficients at (297+/-1 K) are also reported for the relaxation of OH(v=2) by D2O and the relaxation of OD(v=2) by H2O and D2O. The results are discussed in terms of a mechanism involving the formation of hydrogen-bonded complexes in which intramolecular vibrational energy redistribution can occur at rates competitive with re-dissociation to the initial collision partners in their original vibrational states. New ab initio calculations on the H2O-HO system have been performed which, inter alia, yield vibrational frequencies for all four complexes: H2O-HO, D2O-HO, H2O-DO and D2O-DO. These data are then employed, adapting a formalism due to Troe (J. Troe, J. Chem. Phys., 1977, 66, 4758), in order to estimate the rates of intramolecular energy transfer from the OH (OD) vibration to other modes in the complexes in order to explain the measured relaxation rates-assuming that relaxation proceeds via the hydrogen-bonded complexes.     

Automated exploration of stable isomers of H(+)(H(2)O)(n) (n = 5-7) via ab initio calculations: An application of the anharmonic downward distortion following algorithm

Isomers of protonated water clusters H(+)(H(2)O)(n) (n = 5-7) have been explored on ab initio potential energy surfaces by means of the anharmonic downward distortion following algorithm. Totally, 9, 24, and 131 isomers for n = 5, 6, and 7, respectively, were located by the automatic exploration, and all of known important isomers previously reported by conventional geometry optimization approaches have been included in the present results. Moreover, structure transitions depending on n and temperature, which were observed by experimental studies, could be reproduced via thermodynamic simulation on the basis of the superposition approach and the present isomer sets.     

Adsorption of O, H, OH, and H2O on Ag(100)

The adsorption of H(2)O and its dissociation products, O, H, and OH, on Ag(100) has been studied using an ab initio embedding method. Results at different sites (atop, bridge, and hollow) are presented. The four-fold hollow site is found to be the most stable adsorption site for O, H, and OH, and the calculated adsorption energies are 87.1, 42.7, and 76.2 kcal mol(-1), respectively. The adsorption energy of water at the atop and bridge sites is almost identical with values of 11.1 and 12.0 kcal mol(-1), respectively. The formation of adsorbed OH species by adsorption of water on oxygen-precovered Ag(100) is predicted to be exothermic by 36 kcal mol(-1).     

Microsolvation of hydrogen sulfide: exploration of H2S.(H2O)n and SH-.H3O+.(H2O)n-1 (n = 5-7) cluster structures on ab initio potential energy surfaces by the scaled hypersphere search method

The scaled hypersphere search method was applied to ab initio potential energy surfaces of the H2S.(H2O)n/SH-.H3O+.(H2O)n-1 system with n = 5-7. Local minima databases including 121, 326, and 553 structures for n = 5-7, respectively, were obtained based on calculations at the MP2/6-311++G(3df,2p)//B3LYP/6-31+G** level. In these small cluster sizes, the SH-.H3O+.(H2O)n-1 type is still unstable relative to the H2S.(H2O)n type, and the global minima for H2S.(H2O)n are very similar to those of pure water clusters of (H2O)n+1. Thermodynamic simulations based on the present databases showed a structure transition from the well-mixed (H2O)n+1-like global minimum at low temperatures to unmixed complexes between H2S and (H2O)n at high temperatures.     

Communication: determination of the bond dissociation energy (D0) of the water dimer, (H2O)2, by velocity map imaging

The bond dissociation energy (D(0)) of the water dimer is determined by using state-to-state vibrational predissociation measurements following excitation of the bound OH stretch fundamental of the donor unit of the dimer. Velocity map imaging and resonance-enhanced multiphoton ionization (REMPI) are used to determine pair-correlated product velocity and translational energy distributions. H(2)O fragments are detected in the ground vibrational (000) and the first excited bending (010) states by 2 + 1 REMPI via the C? (1)B(1) (000) ← X? (1)A(1) (000 and 010) transitions. The fragments' velocity and center-of-mass translational energy distributions are determined from images of selected rovibrational levels of H(2)O. An accurate value for D(0) is obtained by fitting both the structure in the images and the maximum velocity of the fragments. This value, D(0) = 1105 ± 10 cm(-1) (13.2 ± 0.12 kJ/mol), is in excellent agreement with the recent theoretical value of D(0) = 1103 ± 4 cm(-1) (13.2 ± 0.05 kJ∕mol) suggested as a benchmark by Shank et al. [J. Chem. Phys. 130, 144314 (2009)].     

Quantum dynamics of proton migration in H2O dications: H2+ formation on ultrafast timescales

Irradiation of isolated water molecules by few-cycle pulses of intense infrared laser light can give rise to ultrafast rearrangement resulting in formation of the H(2) (+) ion. Such unimolecular reactions occur on the potential energy surface of the H(2)O(2+) dication that is accessed when peak laser intensities in the 10(15) W cm(-2) range and pulse durations as short as 9-10 fs are used; ion yields of ~1.5% relative to the H(2)O(+) ion are measured. We also study such reactions by means of time-dependent wavepacket dynamics on an ab initio potential energy surface of the dication and show that a proton, generated from O-H bond rupture, migrates towards the H-atom, and forms vibrationally excited H(2)(+) in a well-defined spatial zone.     

A six-dimensional H(2)-H(2) potential energy surface for bound state spectroscopy

We present a six-dimensional potential energy surface for the (H(2))(2) dimer based on coupled-cluster electronic structure calculations employing large atom-centered Gaussian basis sets and a small set of midbond functions at the dimer's center of mass. The surface is intended to describe accurately the bound and quasibound states of the dimers (H(2))(2), (D(2))(2), and H(2)-D(2) that correlate with H(2) or D(2) monomers in the rovibrational levels (v,j)=(0,0), (0,2), (1,0), and (1,2). We employ a close-coupled approach to compute the energies of these bound and quasibound dimer states using our potential energy surface, and compare the computed energies for infrared and Raman transitions involving these states with experimentally measured transition energies. We use four of the experimentally measured dimer transition energies to make two empirical adjustments to the ab initio potential energy surface; the adjusted surface gives computed transition energies for 56 experimentally observed transitions that agree with experiment to within 0.036 cm(-1). For 26 of the 56 transitions, the agreement between the computed and measured transition energies is within the quoted experimental uncertainty. Finally, we use our potential energy surface to predict the energies of another 34 not-yet-observed infrared and Raman transitions for the three dimers.