Low-lying quartet electronic states of nitrogen dioxide

The environmentally active molecule nitrogen dioxide (NO2) has been systematically studied using high level theoretical methods. The electronic ground state and the low-lying quartet states of NO2 have been investigated. Single reference restricted open-shell self-consistent field (SCF), complete active space SCF (CASSCF), spin-restricted (R) and spin-unrestricted (U) configuration interaction with single and double excitations (CISD), coupled cluster with single and double excitations (CCSD), CCSD with perturbative triple excitations [CCSD(T)], and internally contracted multireference configuration interaction (ICMRCI) methods along with Dunning's correlation consistent polarized valence cc-pVXZ and augmented cc-pVXZ (where X=T,Q,5) basis sets were used in this research. At the aug-cc-pV5Z/UCCSD(T) level the classical adiabatic excitation energies (Te values) of the three lowest-lying quartet excited states were predicted to be 83.3 kcalmol (3.61 eV, 29 200 cm(-1)) for the ? 4A2 state, 93.3 kcalmol (4.05 eV, 32 600 cm(-1)) for the b 4B2 state, and 100.8 kcalmol (4.37 eV, 35 300 cm(-1)) for the c 4A1 state. The quantum mechanical excitation energies (T 0 values) were determined to be 81.6 kcalmol (3.54 eV, 28 500 cm(-1)) for the a 4A2 state and 90.7 kcalmol (3.93 eV, 31 700 cm(-1)) for the b 4B2 state. The lowest quartet linear Renner-Teller 4Pi state gives rise to the a 4A2 state with 112.8 degrees and the b 4B2 state with 124.4 degrees <(ONO) bond angles upon bending. The b state shows some peculiar behavior. Although CASSCF, RCISD, UCISD, RCCSD, UCCSD, and RCCSD(T) methods predicted the presence of a Cs equilibrium geometry (a double minimum 4A' state), SCF, UCCSD(T), and ICMRCI wave functions predicted the C2v structure for the b 4B2 state. The importance of both dynamical and nondynamical correlation treatments for the energy difference between C2v and Cs structures of b state is highlighted in this context. The c 4A1 state is predicted to have a very small bond angle of 85.8 degrees . Potential energy diagrams with respect to the bond angles of the ground state and four quartet states are presented.     

Evidence for different types of water participation in the solvolysis of 1-adamantyl, tert-butyl, and methyl chlorides from density functional theory computations

[structure: see text] The activation energy in the gas phase (deltaE(double dagger)) and the free energy of activation (deltaG(double dagger)) in water solution for the hydrolysis of the monohydrates of methyl chloride (MeCl), tert-butyl chloride (t-BuCl), and 1-adamantyl chloride (AdCl) have been computed with the B3LYP/631-G(d) method and the polarizable continuum (PCM) solvation model. There is a fair agreement between the deltaG(double dagger) values computed by us and the experimental data. The mechanistic implications of our computations are in severe contradiction with conventional representations. Thus, the computed nucleophilic solvent assistance (NSA) for the backside attack of a water molecule in the hydrolysis of MeCl is slightly lower than the corresponding NSA for t-BuCl. Hence, the hydrolysis of both MeCl and t-BuCl takes place mainly according to the classical S(N)2 mechanism. The most relevant difference is that deltaG(double dagger) for the frontside attack of water to t-BuCl is disfavored only by ca. 2 kcal/mol with regard to the backside attack but by ca. 23 kcal/mol in the case of MeCl. The higher solvolysis rate in water of t-BuCl in relation to AdCl is not due to steric factors affecting the specific solvation of the corresponding transition states, but to differential bulk solvent effects, which are accounted for by the PCM model.     

A theoretical study of rate coefficients for the O + NO vibrational relaxation

Quasi-classical trajectories have been integrated to study the vibrational relaxation of the O + NO(v) process as a function of the initial vibrational quantum number for T = 298 K, 1500 K, and 3000 K. Two reliable potential energy surfaces have been employed for the A' and A' doublet states of NO2. The calculated vibrational relaxation rate constants show a nearly v-independent behavior at room temperature and a moderate increase with v for higher temperatures. Although deviating significantly from the recommended values, good agreement with recent experimental results has been obtained. The importance of multi-quantum transitions is also analyzed.     

Nonequilibrium vibrational excitation of OH radicals generated during multibubble cavitation in water

The sonoluminescence (SL) spectra of OH(A(2)Σ(+)) excited state produced during the sonolysis of water sparged with argon were measured and analyzed at various ultrasonic frequencies (20, 204, 362, 609, and 1057 kHz) in order to determine the intrabubble conditions created by multibubble cavitation. The relative populations of the OH(A(2)Σ(+)) v' = 1-4 vibrational states as well as the vibronic temperatures (T(v), T(e)) have been calculated after deconvolution of the SL spectra. The results of this study provide evidence for nonequilibrium plasma formation during sonolysis of water in the presence of argon. At low ultrasonic frequency (20 kHz), a weakly excited plasma with Brau vibrational distribution is formed (T(e) ~ 0.7 eV and T(v) ~ 5000 K). By contrast, at high-frequency ultrasound, the plasma inside the collapsing bubbles exhibits Treanor behavior typical for strong vibrational excitation. The T(e) and T(v) values increase with ultrasonic frequency, reaching T(e) ~ 1 eV and T(v) ~ 9800 K at 1057 kHz.     

High-resolution millimeter wave spectroscopy and multichannel quantum defect theory of the hyperfine structure in high Rydberg states of molecular hydrogen H2

Experimental and theoretical methodologies have been developed to determine the hyperfine structure of molecular ions from detailed studies of the Rydberg spectrum and have been tested on molecular hydrogen. The hyperfine structure in l=0-3 Rydberg states of H2 located below the X 2Sigmag+(v+=0,N+=1) ground state of ortho H2+ has been measured in the range of principal quantum number n=50-65 at sub-MHz resolution by millimeter wave spectroscopy following laser excitation to np and nd Rydberg states using a variety of single-photon and multiphoton excitation sequences. The np1(1), nd1(1), and the nf1(0-3) Rydberg states were found to be metastable and to have lifetimes of more than 5 micros beyond n=50. Members of other series, such as the nd1(2), nd1(3), and the np1(0) series, were found to have lifetimes of more than 1 mus. Local perturbations induced by low-n Rydberg states belonging to series converging on rovibrationally excited levels of H2+ reduce the lifetimes in narrow ranges of n values. The hyperfine structure is strongly dependent on the value of the orbital angular momentum l. In the penetrating s and p states at n approximately 50 the exchange interaction dominates over the hyperfine interaction and the levels can be labeled by the total electron spin angular momentum quantum number S (S=0 or 1). In the less penetrating d and f Rydberg states, the hyperfine interaction between the core nuclear and electron spins is larger than the exchange interaction and the Rydberg states are of mixed singlet and triplet character. A procedure based on the Stark effect and on the systematic analysis of selection rules and combination differences was developed to determine the orbital and the total angular momentum quantum numbers l and F and to construct an energy map of p and f Rydberg levels between n=54 and 64 with relative positions of an accuracy of better than 1 MHz. Multichannel quantum defect theory (MQDT) was extended to treat the hyperfine structure in molecular Rydberg states and was used to analyze the observed hyperfine structure of the p and f Rydberg states of H2. The frame transformation between the Born-Oppenheimer channels described by the angular momentum coupling scheme (abetaJ) and the asymptotic channels described by the (e[bbetaS+]) coupling scheme was derived and enables an elegant treatment of all intermediate coupling cases. Purely ab initio quantum defect theory reproduced the experimentally determined positions to within 40 MHz for the p levels and 13 MHz for the f levels. By slight adjustments of the quantum defect functions and their energy dependences and by consideration of the p-f interaction, of the singlet-triplet splittings of the f levels, and of the departure of the ionic levels from pure coupling case (bbetaS+), the agreement between theory and experiment could be improved to 600 kHz. By comparing the results of MQDT calculations of the hyperfine structure of f Rydberg levels with those of coupled equations calculations, the frame transformation approximation of MQDT was shown to be accurate to within 300 kHz. The extrapolated ionic hyperfine structure of the X 2Sigmag+(v+=0,N+=1) ionic level corresponds to the ab initio prediciton of Babb and Dalgarno [Phys. Rev. A 46, R5317 (1992)] within the experimental error.     

Electronic spectrum of the AlC(2) radical

An electronic transition of the AlC2 radical (C2v structure) has been observed using laser-induced fluorescence spectroscopy. The molecule was prepared in a supersonic expansion by ablation of an aluminum rod in the presence of acetylene gas. A spectrum was recorded in the 451-453 nm region and assigned to the C 2B2-X 2A1 system (T0 = 22,102.7 cm(-1)) based on a rotational analysis and agreement with calculated molecular parameters and excitation energies. Ab initio results obtained using couple cluster methods are in accord with previous theoretical work which concludes that ground-state AlC2 possesses a T-shaped C2v 2A1 geometry, with the linear 2Sigma+ AlCC isomer 0.70 eV higher in energy. A fit of the experimental spectrum yields rotational constants in the ground and electronically excited states that are in reasonable agreement with the calculated values: A' = 1.7093(107), B' = 0.4052(50), C' = 0.3228(49) cm(-1) for the X 2A1 state, and A' = 1.5621(137), B' = 0.4028(46), C' = 0.3201(54) cm(-1) for C 2B2. Variation in individual fluorescence lifetimes suggests that the emitting C 2B2 state undergoes rovibronic mixing with lower lying electronic states.     

Observation of the nd 1Delta(g) (n = 6, 7, and 8) Rydberg states of Na2 by optical-optical double resonance spectroscopy: L uncoupling and perturbations

The nd (1)Delta(g) (n = 6, 7, and 8) Rydberg states of Na(2) correlating with the asymptotic limits of Na(3s) + Na(nd) have been observed using high-resolution cw optical-optical double resonance spectroscopy corresponding to the rovibrational transitions X (1)Sigma(g)(+)(v("),J(")) + h nu(pump) --> B (1)Pi(u)(v('),J(')) + h nu(probe) --> nd (1)Delta(g)(v,J). Totals of 104, 83, and 45 identified rovibrational e/f-parity levels in the ranges of v = 0-11, 11 < or = J < or = 83; v = 0-10, 11 < or = J < or = 83; and v = 0-10, 11 < or = J < or = 65, have been assigned to the 6d (1)Delta(g), 7d (1)Delta(g), and 8d (1)Delta(g) states, respectively. Using the observed quantum levels, molecular constants were determined from the Dunham fits of the e-parity levels and the Rydberg-Klein-Rees potential curves were constructed for the nd (1)Delta(g)(n = 6-8) states. The characteristics of the estimated Lambda-doubling splitting constants (q(0), q(v), and mu) with n(= 5-8) of the nd (1)Delta(g) series have been explored. Detailed investigations reveal that the nd (1)Delta(g)(n = 6-8) states involve L uncoupling from the internuclear axis and each of these states is affected by an asymmetric perturbation caused by the up and down adjacent states. The rotational-branch intensity and position anomalies in the observed spectra of the nd (1)Delta(g) series (n = 5-8) of Na(2) lead to the conclusion that due to the effects of the L-uncoupling perturbations, the same l complexes approaching the same ion-core limits result in the same l-mixing processes which lead to the formation of the supercomplexes due to the anisotropy of the molecular-ion [Na(2)(+)(3s)] field. This would open up opportunities to study the effects of L uncoupling and perturbations in the nd series and high Rydberg states of other alkali dimers.     

d-orbital effects on stereochemical non-rigidity: twisted Ti(IV) intramolecular dynamics

The isomerization dynamics of tris-catecholate complexes have been investigated by variable-temperature NMR methods, demonstrating that the intramolecular racemization of Delta and Lambda enantiomers of d0 Ti(IV) is facile and faster than that of d10 Ga(III) and Ge(IV) analogues. Activation parameters for the racemization of K2[Ti2(3)] (H(2)2 = 2,3-dihydroxy-N,N'-diisopropylterephthalamide) were determined from line shape analysis of 1H NMR spectra [methanol-d4: deltaH++ = 47(1) kJ/mol; deltaS++ = -34(4) J/mol K; deltaG++(298) = 57(3) kJ/mol; DMF-d7: deltaH++ = 55(1) kJ/mol; deltaS++ = -16(4) J/mol K; deltaG++(298) = 59(3) kJ/mol; D2O (pD* = 8.6, 20% MeOD): deltaH++ = 48(3) kJ/mol; deltaS++ = -28(10) J/mol K; deltaG++(298) = 56(3) kJ/mol]. The study of K2[Ti4(3)] (H(2)4 = 2,3-dihydroxy-N-tert-butyl-N'-benzylterephthalamide) reveals two distinct isomerization processes: faster racemization of mer-[Ti4(3)]2- by way of a Bailar twist mechanism (D3h transition state) [T(c) approximately 242 K, methanol-d4], and a slower merright harpoon over left harpoonfac [Ti4(3)]2- isomerization by way of a Ray-Dutt mechanism (C2v transition state) [T(c) approximately 281 K, methanol-d4]. The solution behavior of the Ti(IV) complexes mirrors that reported previously for analogous Ga(III) complexes, while that of analogous Ge(IV) complexes was too inert to be detected by 1H NMR up to 400 K. These experimental findings are augmented by DFT calculations of the ML3 ground states and Bailar and Ray-Dutt transition states, which correctly predict the relative kinetic barriers of complexes of the three metal ions, in addition to faithfully reproducing the ground-state structures. Orbital calculations support the conclusion that participation of the Ti(IV) d orbitals in ligand bonding contributes to the greater stabilization of the prismatic Ti(IV) transition states.     

Design of oxygen reduction bimetallic catalysts: ab-initio-derived thermodynamic guidelines

The Gibbs free energies of key elementary steps for the electrocatalytic oxygen reduction reaction (ORR) are calculated with B3LYP type of density functional theory: O2 + M + H+ + e- (0 eV) --> HOO-M (deltaG1), HOO-M + M --> HO-M + O-M (deltaG2), O2 + 2M + H+ + e- (0 eV) --> O-M + HO-M (deltaG3), and HO-M + O-M + 3H+ + 3e- (0 eV) --> 2H2O + 2M (deltaG4), where H+ is modeled as H3(+)O(H2O)3 and M stands for the adsorption site of a metal catalyst modeled by a single metal atom as well as by an M3 cluster. Taking Pt as a reference, deltaG4 is plotted against deltaG1 for 17 metals from groups V to XII. It is found that no single metal has both deltaG1 and deltaG4 more negative than Pt, although some of them have either more negative deltaG1 or more negative deltaG4. This enables us to explain thermodynamically why no other single metal catalyzes the ORR as effectively as Pt does. Moreover, a thermodynamic analysis reveals that the signs of delta deltaG (the difference between deltaG of other metals and deltaG of Pt) strongly correlate with the valence electronic structure of metals, i.e., delta deltaG1 < 0 and delta deltaG4 > 0 for metals M with vacant valence d orbitals, whereas delta deltaG1 > 0 and delta deltaG4 < 0 for metals M' with fully occupied valence d orbitals. Thus, a simple thermodynamic rule for the design of bimetallic catalysts for the ORR is proposed: couple a metal M (delta deltaG1 < 0) with a second metal M' (delta deltaG4 < 0) to form an alloy catalyst MM'3. The rationale behind this selection is based on M being more efficient for the rate-determining step, i.e., for the formation of the adsorbed species M-OOH, while M' can enhance the reductions of O and OH in the last three electron-transfer steps.     

Approximation function for the direct calculation of rate constants and Gibbs activation energies of enantiomerization of racemic mixtures from chromatographic parameters in dynamic chromatography

An approximation function for enantioselective dynamic chromatography of racemic mixtures of interconverting enantiomers has been derived that allows the direct calculation of enantiomerization rate constants (k1 and k(-1)) and Gibbs activation energies of enantiomerization, deltaG++ , from chromatographic parameters, i.e., retention times of the enantiomers A and B ((t(A)R and t(B)R), peak widths at half height (WA and wB) and the relative plateau height (hplateau), without computer simulation. The reaction rate constants of enantiomerization, k(-1), obtained with this approximation function, have been validated by comparison with a simulated dataset of 15,625 chromatograms. The mean, standard deviation and confidence interval show a high correlation between the approximated and simulated rate constants. The average deviation from the Gibbs activation enthalpy of enantiomerization, deltaG++, has been estimated to be as small as about +/- 0.11 RT.     

Laser induced amplified spontaneous emission from the B2Pi, L2Pi, and I2Sigma+ valence states of NO

Amplified spontaneous emission (ASE) from single rovibrational levels of valence (non-Rydberg) states of NO molecules has been investigated. The B2Pi (v=24 and 25), L2Pi (v=5 and 6), and I2Sigma+ (v=6) levels have been populated through laser optical-optical double resonance excitation via the Rydberg A2Sigma+ state. Term values for the 2Pi states have been determined with an accuracy of +/-0.03 cm(-1). Analyses of rotationally resolved dispersed ASE spectra in the near infrared region have shown that all the lower states belonged to the Rydberg states. The valence approximately Rydberg coupling in the upper manifolds has driven ASE systems from the valence to the Rydberg levels where they benefit from the strong intensities of inter-Rydberg transitions with Deltav=0. The experimentally predicted valence approximately Rydberg interactions have been compared with theoretical treatments.     

Computations of 36 tautomer/isomer equilibria of different lactams

Gas-phase energies of 36 tautomer/isomer pairs of 18 six-membered N-heterocyclic compounds were computed quantum chemically. Among the considered B3LYP, BH&HLYP, BH&HLYP(G), and PW6B95 DFT functionals, the latter two provide accurate tautomer/isomer pair energies with root-mean-square deviations (rmsd) relative to experiments of 0.2 and 0.3 kcal/mol, respectively. Since only few (namely five) experimental data are available, 15 tautomer/isomer pair energies were computed with the very precise QCISD(T)(quadruple-ζ) method serving as reference. Relative to this reference the PW6B95 DFT functional is slightly superior to the BH&HLYP(G) functional, yielding an rmsd of 0.7 and 0.8 kcal/mol, respectively. In contrast to BH&HLYP(G), the PW6B95 DFT functional yields also accurate tautomer/isomer pair energies if zwitterionic structures are involved. The tautomer/isomer pair states possess different amounts of aromaticity. This is characterized by nucleus-independent chemical shift (NICS) values. The tautomer/isomer pair reference energies, from which the energies computed with PW6B95 are subtracted, correlate linearly with the corresponding differences in the NICS values. This correlation is used to construct a correction term for the pair energies computed with PW6B95, yielding tautomer/isomer pair energies with rmsd of 0.3 kcal/mol with respect to the more CPU time demanding QCISD(T)(quadruple-ζ) method.     

Computational studies of (HIO3) isomers and the HO2 + IO reaction pathways

The geometries, harmonic vibrational frequencies, relative energetics, and enthalpies of formation of (HIO3) isomers have been examined using quantum mechanical methods. At all levels of theory employed, MP2, B3LYP, and CCSD(T), the lowest energy structure is found to be the HOIO2 form, which shows particular stability. The two isomers HOOOI and HOOIO are closely located at the CCSD(T) level of theory. The higher energy structure is HIO3. The interisomerization transition states have been determined, along with the transition states involved in the various pathways of the reaction HO2 + IO.     

A new ab initio potential energy surface for studying vibrational relaxation in NO(v) + NO collisions

A new ab initio potential energy surface for the ground state of the NO-NO system has been calculated within a reduced dimensionality model. We find an unusually large vibrational dependence of the interaction potential which explains previous spectroscopic observations. The potential can be used to model vibrational energy transfer, and here we perform quantum scattering calculations of the vibrational relaxation of NO(v). We show that the vibrational relaxation for v = 1 is 4 orders of magnitude larger than that for the related O(2)(v) + O(2) system without having to invoke nonadiabatic mechanisms as had been suggested in the past. For highly vibrationally excited states, we predict a strong dependence of the rates on the vibrational quantum number as has been observed experimentally, although there remain important quantitative differences. The importance of a chemically bound isomer on the relaxation mechanism is analyzed, and we conclude it does not play a role for the values of v considered in the experiment. Finally, the intriguing negative temperature dependence of the vibrational relaxation rate constants observed in experiments was studied using an statistical model to include the presence of many asymptotically degenerate spin-orbit states.     

Ab initio calculations of the electronic states of AsH2 including dissociation characteristics

Multireference configuration interaction calculations have been carried out for low-lying electronic states of AsH(2). Bending potentials for the ten lowest states of AsH(2) are obtained in C(2v) symmetry for As-H distances fixed at the the ground state equilibrium value of 2.845 a(0), as well as for the minimum energy path constrained to R(1) = R(2). The calculated equilibrium geometries for the X?(2)B(1) ground state and the A?(2)A(1) excited state agree very well with the previous experimental and theoretical results, whereas the data for the higher-lying states are obtained for the first time. Asymmetric potential energy surface (PES) cuts (at R(1) = 2.845 a(0), θ = 90.7°) and two-dimensional (2D) PESs for the lowest three states are also new. The calculated ab initio data are used for analysis of possible AsH(2) photodissociation channels and predissociation effects. It is shown that the A?(2)A(1)-X?(2)B(1) transition dipole moment decreases with increasing bending angle, which influences the intensity distribution in the A?(0,0,0)→X? emission spectrum (v(2)' bending series), shifting its maximum to smaller v(2)' quantum numbers.     

Quantum-admixture model of high-spin <--> low-spin transition for ferrous complex molecules

A quantum-admixture model for the d(6) configuration ferrous complex molecules with the high-spin <--> low-spin transition has been established by using the unified crystal-field-coupling (UCFC) scheme. A general study has been made on the spin transition of octahedrally coordinated d(6) complexes, and a special application has been given to an Fe(II) compound Fe(II)(TRIM)(2)(PhCO(2))(ClO(4)). The results show the following: (i) The quantum picture of the spin transition of a d(6) system, such as Fe(II), is much more complex than a simple transition between the pure (5)T(2g) and (1)A(1g) states as usually understood. In practice, owing to spin-orbit coupling, spin is no longer a good quantum number and there is no longer a pure (5)T(2g) or (1)A(1g) state. Each of them splits into substates and each substate is a linear combination of various multiplets. The high-spin --> low-spin transition of an octahedrally coordinated d(6) ion is practically the crossover of the two lowest substates of (5)T(2g) at the critical point. (ii) At the spin-transition critical point the magnetic moment mu(eff) approximately 5.22 mu(B), which is obviously different from the simple average of the mu(eff) values of high-spin and low-spin states but near the saturation value. (iii) The calculation of the effective molecular magnetic moment mu(eff) for an octahedrally coordinated Fe(II) ion shows that the mu(eff)-T curve is in good agreement with Lemercier et al.'s experiment and both the low-spin value mu(eff) = 0.51 mu(B) and the high-spin value mu(eff) = 5.4 mu(B) are comparable with the experimental values 0.76 mu(B) and 5.4 mu(B), respectively. (iv) The T dependence of the crystal field parameter Dq in the spin-transition region is approximately linear.     

Excited states of phosphorescent platinum(II) complexes containing N wedge C wedge N-coordinating tridentate ligands: spectroscopic investigations and time-dependent density functional theory calculations

The absorption and emission spectra of the Pt(II) complexes containing N wedge C wedge N-coordinating tridentate ligands, platinum(II) 1,3-di(2-pyridyl)benzene chloride [Pt(dpb)Cl] and platinum(II) 3,5-di(2-pyridyl)toluene chloride [Pt(dpt)Cl], together with their corresponding free ligands, 1,3-di(2-pyridyl)benzene (dpbH) and 3,5-di(2-pyridyl)toluene (dptH), have been analyzed by density functional theory (DFT) for the ground state and time-dependent DFT (TDDFT) for the excited states. T(1)(A(1)) and S(1)(B(2)) of the complexes (in C(2)(v) symmetry) were assigned on the basis of the calculated excitation energies as well as comparison of the experimental spectroscopic properties and the calculated states' characteristics. The calculated excitation energies for T(1) and S(1) of the complexes as well as those for T(1) of the free ligands were in good agreement with their observed values within 600 cm(-1). The d-pi* characters of the excited states were evaluated from the change in electron densities between the ground and excited states by Mulliken population analysis; values of 25% for T(1) and 32% for S(1) were obtained for both complexes. The calculated values of d-pi* character were found to be consistent with the reported emission lifetimes as well as the observed emission energy shifts from the corresponding free ligands. Most spectroscopic properties of the complexes and the free ligands, which include solvatochromic shift, Stokes shifts, methyl substitution shifts, and emission spectra profiles, were well explained from the calculation results.     

Water dimers in the atmosphere III: equilibrium constant from a flexible potential

We present new results for the water dimer equilibrium constant K(p)(T) in the range 190-390 K, using a flexible potential energy surface fitted to spectroscopical data. The increased numerical complexity due to explicit consideration of the monomer vibrations is handled via an adiabatic (6 + 6)d decoupling between intra- and intermolecular modes. The convergence of the canonical partition function of the dimer is ensured by computing all energy levels up to dissociation for total angular momentum values J = 0-5 and using an extrapolation scheme to higher values. The newly calculated values for K(p)(T) are in very good agreement with available experimental data at room temperature. At higher temperatures, an analysis of the convergence of the partition function reveals that quasi-bound states are likely to contribute to the equilibrium constant. Additional thermodynamical quantities (deltaG, deltaH, deltaS, and C(p)) have also been determined and fit to quadratic expressions a + bT + cT2.