A&A 460, 323-329 (2006)
DOI: 10.1051/0004-6361:20065804

Influence of a new potential energy surface on the rotational (de)excitation of H2O by H2 at low temperature

M.-L. Dubernet1 - F. Daniel1 - A. Grosjean2 - A. Faure3 - P. Valiron3 - M. Wernli3 - L. Wiesenfeld3 - C. Rist3 - J. Noga4 - J. Tennyson5


1 - Observatoire de Paris-Meudon, LERMA UMR CNRS 8112, 5 Place Jules Janssen, 92195 Meudon Cedex, France
2 - Laboratoire d'Astrophysique, Observatoire de Besançon, UMR CNRS 6091, Université de Franche-Comté, 41 bis avenue de l'Observatoire, BP 1615, 25010 Besançon Cedex, France
3 - Laboratoire d'Astrophysique, Observatoire de Grenoble, UMR CNRS 5571, Université Joseph Fourier, 38041 Grenoble Cedex 09, France
4 - Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynska dolina CH1, 84215 Bratislava, Slovakia
5 - Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK

Received 12 June 2006 / Accepted 2 August 2006

Abstract
Aims. Using a newly determined 5D potential energy surface for H2-H2O we provide an extended and revised set of rate coefficients for de-excitation of the lowest 10 para- and 10 ortho- rotational levels of H2O by collisions with para-(j=0) and ortho-H2(j=1), for kinetic temperatures from 5 K to 20 K.
Methods. Our close coupling scattering calculations involve a slightly improved set of coupled channels with respect to previous calculations. In addition, we discuss the influence of several features of this new 5D interaction on the rotational excitation cross sections.
Results. The new interaction potential leads to significantly different rate coefficients for collisions with para-H2 (j=0). In particular the de-excitation rate coefficient for the 110 to 101 transition is increased by up to 300% at 5 K. At 20 K this increase is 75%. Rate coefficients for collisions with ortho-H2(j=1) are modified to a lesser extent, by up to 40%. The influence of the new potential on collisions with both para-(j=0) and ortho-H2(j=1) is expected to become less pronounced at higher temperatures.

Key words: molecular data - molecular processes - ISM: molecules

1 Introduction

Water is a key molecule for the chemistry and the energy balance of the gas in cold clouds and star forming regions, thanks to its relatively large abundance and large dipole moment. A wealth of observational data has already been obtained by the Infrared Space Observatory (see for example: Cernicharo et al. 2006; Tsuji 2001; Spinoglio et al. 2001; Cernicharo & Crovisier 2005; Wright et al. 2000), the Submillimeter Wave Astronomy Satellite (Melnick et al. 2000) and the ODIN satellite (Sandqvist et al. 2003; Wilson et al. 2003). In the near future the Heterodyne Instrument for the Far-Infrared (HIFI) will be launched on board the Herschel Space Observatory. It will observe with unprecedented sensitivity spectra of many molecules with an emphasis on water lines in regions from cold molecular clouds to star forming regions. The interpretation of these spectra will rely upon the accuracy of the available collisional excitation rate coefficients that enter into the population balance of the emitting levels of the molecules. In temperature range from 5 K to 1500 K the most abundant collider likely to excite the molecules is the hydrogen molecule, followed by the helium atom. In diffuse clouds and photo-ionized environments, water collisions with atomic H and electrons might also play a significant role. Much work has dealt with the rotational or ro-vibrational excitation of many molecules by He or H2, always trying to improve accuracy (see bibliography in the BASECOL database, Dubernet et al. 2006). Some of the latest calculations (Lique et al. 2005; Dubernet 2005; Daniel et al. 2004,2005b) have been used to model astrophysical observations and have shown that careful calculation of the excitation rate coefficients can influence the interpretation of spectra (Daniel et al. 2006,2005a; Lique & Cernicharo 2006).

Until recently, only calculation of the excitation of water by He was feasible and a series of such calculations were performed in order to determine accurate rate coefficients (Green et al. 1993; Maluendes et al. 1992; Palma et al. 1988b,1989,1988a; Green 1980). These studies involved improving the potential energy surface (PES), improving the quality of the scattering calculations and extending previous work either to a larger ranges of temperature or to more transitions.

A pioneering rigid-body 5D PES was obtained by Phillips et al. (1994) for the excitation of the rotational levels of H2O by H2. Using this PES,  Phillips et al. (1996) computed rate coefficients for temperatures ranging from 20 K to 140 K. Dubernet & Grosjean (2002) and Grosjean et al. (2003) extended this work down to 5 K and pointed out that such low temperature rates are highly sensitive to a proper description of resonances.

Recently an accurate 9D PES for the deformable H2-H2O system was calculated by Faure et al. (2005a). This PES combined a) conventional 5D and 9D CCSD(T) calculations (coupled cluster with perturbative triples), and b) accurate calibration data using the explicitely correlated CCSD(T)-R12 approach (Noga & Kutzelnigg 1994).

As a first application of the 9D PES, high temperature ( 1500 < T < 4000 K) rate constants for the relaxation of the $\nu_2$ bending mode of H2O were estimated from quasi classical trajectory calculations (Faure et al. 2005a) and the role of rotation in the vibrational relaxation of water was emphasized (Faure et al. 2005b).

Another application of this 9D PES was to construct an accurate 5D PES suitable for inelastic rotational calculations by averaging over H2and H2O ground state vibrational states. As was pointed out in Faure et al. (2005a), this state-averaged PES is very close to a rigid-body PES using state-averaged geometries for H2O and H2.

A goal of this paper is to assess the influence of this new 5D PES on the rotational excitation cross sections. The composite nature of the PES also offers an opportunity to discuss the respective importance of the proper state-averaging (or of the proper choice of a rigid-body geometry) and of the R12 corrections towards the CCSD(T) infinite basis set limit.

A second goal is to use the new 5D PES of Faure et al. (2005a) to provide an extended set of rate coefficients between 5 K and 20 K for de-excitation of the lowest 10 para- and 10 ortho- rotational levels of H2O by collisions with para-(j=0) and ortho-H2(j=1). This set includes more rotational transitions than studied by Dubernet & Grosjean (2002) and Grosjean et al. (2003) with the rigid 5D potential energy surface of Phillips et al. (1994). The present calculations therefore use a larger rotational basis set for water in order to ensure convergence for transitions among the highest rotational levels. Another advance on previous calculations (Phillips et al. 1996; Dubernet & Grosjean 2002; Grosjean et al. 2003) is the use of experimental rotational energies (Dabrowski 1984) for the H2 monomer. Non-equilibrium H2 ortho-to-para ratios are observed in astrophysics so that the rate coefficients for excitation by ortho-H2(j=1) could be of relevance at temperatures as low as 20 K.

Such highly accurate quantum scattering calculations are extremely computer-time consuming above 20 K. Therefore we present results now in order to make them available to the astronomical community in a timely manner. In the 5 K to 20 K range, contributions from the excited rotational states of H2are negligible and it is therefore computationally tractable to use both the best converged basis set and the highest accuracy scattering method. Moreover, the influence of the new PES is expected to be most significant in this low temperature range.

In Sect. 2 we make comparaisons to show the effects of the new PES and in Sect. 3 we describe our methods and present our new results.

2 Influence of the new 5D potential energy surface

Phillips et al. (1994) calculated the interaction potential for 722 points on the surface and for fixed theoretical equilibrium geometries $r_{\rm {HH}} = 1.402 ~a_0$ , $r_{\rm {OH}} = 1.809 ~a_0$ and HOH angle equal to 104.52$^\circ$. The PES was obtained by means of fourth-order pertubation theory (MP4) and was expanded over 48 angular terms. Faure et al. (2005a) followed a three-step procedure: (i) a 5D rigid-rotor PES reference was computed at the CCSD(T) level of theory for fixed theoretical equilibrium geometries $r_{\rm {HH}} = 1.400 ~a_0$, $r_{\rm {OH}} = 1.809 ~a_0$ and HOH angle equal to 104.22$^\circ$; (ii) this reference surface was then calibrated to a few cm-1 accuracy using 812 "high cost'' CCSD(T)-R12 calculations; (iii) the R12-corrected rigid rotor surface was then extended to 9D using a new set of CCSD(T) calculations. The resulting 9D PES contains all relevant information to describe the interaction of all H2O-H2 isotopologues in zero-point or excited vibrational states. This 9D PES is subsequently averaged over the ground vibrational state of H2O and H2 in order to obtain the final average 5D CCSDT(T)-R12 PES. This latter 5D PES involves a conservative angular expansion over 149 angular functions $T_{p_1 q_1 p_2 p} (\Theta, \Phi, \Theta ', \Phi ', \alpha,
\beta, \gamma)$ (according to the conventions in  Phillips et al. (1994)). A subset of 83 functions was selected, see Table 1, using a self-consistent Monte Carlo estimator and was checked to be sufficient for scattering calculations (see details below).

Table 1: The 83 angular coefficients kept in the angular expansion of the potential energy surface of Faure et al. (2005a).

Figure 1 gives a comparison of the leading radial coefficients from Phillips et al. (1994) with the corresponding terms from the PES of Faure et al. (2005a). Faure et al.'s terms V1203 and V2224 are more repulsive by a few cm-1, their isotropic V0000 term is deeper by about 10 cm-1with a slight shift in the minimum, their V2024 term is smoother pointing to an expansion problem in Phillips et al.'s PES.


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[angle=270]{5804fig1.eps}}\end{figure} Figure 1: Comparison of the leading radial coefficients (in cm-1) of the H2O-H2 interaction, i.e. the isotropic term, the dipole-quadrupole term V1023, the quadrupole-quadrupole terms V2024, V2224, from Phillips et al. (1994) (dashed lines) with the corresponding terms from the PES of Faure et al. (2005a) (solid lines) used in the present calculations. Note: the radial coefficients are divided by $4 \pi $, so that the V0000 term represents the isotropic potential.
Open with DEXTER

2.1 Influence of PES on cross sections

In order to probe the influence of the various PES on the scattering results, the rotational state-to-state cross-sections are obtained with exactly the same methods, parameters and molecular basis set sizes as in Dubernet & Grosjean (2002) and Grosjean et al. (2003).

The PES tests are first carried out on the excitation cross sections of ortho-water with para-H2(j=0), at total energies chosen in the energy range necessary to get rate coefficients up to 20 K, i.e. at total energy of 90 cm-1, 200 cm-1 and 410 cm-1. At these total energies the number of open channels of ortho-H2O is respectively 3, 6 and 14. All cross sections but those at 410 cm-1 and those among the 3 lowest transitions at 200 cm-1 are in the resonant regime.

The tests aim to probe independently the influence of 3 features of the average 5D PES (Faure et al. 2005a): the number of terms kept in the angular expansion of the PES, the R12 correction and the vibrational corrections and at comparing these results with cross-sections calculated using the Phillips et al. (1994) PES. These detailed tests were not extended to rate coefficients, due to the cumbersome energy grid required to accurately reproduce the resonant structures in the cross-sections (Dubernet & Grosjean 2002).

For the three total energies we find that the 83 angular terms of Table 1 are sufficient to achieve an accuracy better than 0.1% on cross sections compared to using 149 angular functions; therefore all further tests and final rate calculations use 83 angular terms in the potential expansion.

Figure 2 compares the isotropic term of the H2O-H2PES for the PES of Phillips et al. (1994) for steps 1, 2 and 3 of the PES of Faure et al. (2005a). Our CCSD(T) reference PES (step 1) is qualitatively similar to the MP4 PES of Phillips et al. (1994), and mostly differs by a fully converged radial and angular coverage. The R12 correction (step 2) is weakly anisotropic. It lowers the isotropic term by about 7  cm-1, and changes the V1203 and V2224 terms in a similar fashion to the vibrational correction (step 3) from R=7  a0 to infinity; it has little effect both at short range for the V1203 and V2224 terms and for the V2024 term over the whole distance range. On the other hand the vibrational correction slightly deepens the isotropic term beyond R=7  a0 while it mostly affects the magnitude of all anisotropic terms at shorter range.


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[angle=270]{5804fig2.eps}}\end{figure} Figure 2: Comparison of the isotropic H2O-H2 interaction terms ( $V_{0000}/4\pi $ in cm-1) for various PES: the solid line shows the final average 5D PES, the dotted line gives Phillips et al.'s PES (Phillips et al. 1994), the initial 5D reference PES (Faure et al. 2005a) (step 1) is represented by the long dashed line, and step 2 where the R12 correction has been added is given by the dashed line.
Open with DEXTER

Figure 3 shows cross sections for the excitation of ortho-H2O by para-H2(j=0) obtained for the 4 levels of PES: (1) corresponds to the PES of Phillips et al. (1994); (2) corresponds to the 5D reference PES (Faure et al. 2005a) (step 1); (3) corresponds to step 2 where the R12 correction has been added; (4) gives the cross-section obtained with the full 5D-PES including the vibrational correction.


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[angle=270]{5804fig3.eps...
...}\par\resizebox{8.8cm}{!}{\includegraphics[angle=270]{5804fig5.eps}}\end{figure} Figure 3: Cross sections of the rotational excitation of ortho-H2O by para-H2(j=0) obtained at total energy of (a) 90  cm-1, (b) 200 cm-1 and (c) 410  cm-1 for the 4 levels of PES by increasing accuracy: (1) is Phillips et al.'s PES (Phillips et al. 1994); (2) corresponds to the 5D reference PES (Faure et al. 2005a) (step 1); (3) corresponds to step 2 where the R12 correction has been added; (4) gives the cross-section obtained with the full PES including the vibrational correction.
Open with DEXTER

For the largest cross sections, differences from 1% to 15% occur between Phillips et al.'s PES (Phillips et al. 1994) and the 5D reference PES (Faure et al. 2005a) (step 1), the R12 correction (step 2) brings differences between 3% and 30% compared to the 5D reference PES, the effect decreasing with increasing energy. The vibrational correction (step 3) induces the largest difference, i.e. up to a factor of 2 with respect to the 5D reference PES at low energy and around 20-30% for the strongest transitions at high energy. Overall the final average 5D PES (Faure et al. 2005a) induces differences in the largest cross sections of up to 40% at 410 cm-1, up to 70% at 200 cm-1, up to a factor 2 at 90 cm-1. Away from the resonance region these large relative variations correspond to small absolute differences (about 1 Å2) on small values of cross sections for collisions with para-H2(j=0) as can be seen in Fig. 3.

The absolute differences found for collisions with ortho-H2(j=1) are about the same order of magnitude, but excitation cross sections by ortho-H2(j=1) are generally an order of magnitude larger than cross sections with para-H2(j=0). Therefore the overall effect of the new PES on the excitation cross section with ortho-H2(j=1) is not large. But it is found that at a total energy of 265 cm-1 and 321.706 cm-1 (equivalent respectively to 143.294 cm-1 and 200 cm-1 for para-H2(j=0)) both the R12 and the vibrational correction have a similar effect on excitation of cross sections of water by ortho-H2(j=1). This is consistent with the fact that the two corrections change the anisotropic terms V1203, V2224 in a similar way for the long range part of the potential. For collisions with para-H2(j=0), there is no straightforward interpretation of the variation of the cross sections with the R12 and vibrational corrections because of the role of j(H2) = 2 closed channels which account for about 20% to 30% of the cross sections. Nevertheless the generally stronger influence on cross section of the vibrational correction to the PES over the R12 correction seems to be related to the larger vibrationally corrected V1203, V2224 and V2024 terms for R < 7  a0.

2.2 Influence of PES on rate coefficients

In order to probe the influence of the new 5D PES (Faure et al. 2005a) we calculate rate coefficients for the excitation of ortho-H2O with para-H2(j=0) and ortho-H2(j=1) from 5 K to 20 K, using the same scattering methods as in Dubernet & Grosjean (2002) and  Grosjean et al. (2003).

Table 2 gives these rate coefficients at 20 K for some transitions, their ratios at 5 K, 10 K, 20 K with the Dubernet and co-workers rate coefficients (Dubernet & Grosjean 2002; Grosjean et al. 2003) (ratios (3)/(2)) and with Phillips et al. (1996) results (ratios (3)/(1)) and the ratios of the rate coefficients of Dubernet & Grosjean (2002); Grosjean et al. (2003) to those of Phillips et al. (1996) (ratios (2)/(1)).

Table 2: Ratios of de-excitation rates for ortho-H2O colliding with para-(j=0)and ortho-H2(j=1). Part a illustrates the effects at 20 K of an improved treatment of both scattering and PES (Col. (3)/(1)) and of the sole improvement of the scattering treatment (Col. (2)/(1)). The levels are labelled with jK-1K1. The new rates are given for reference (Col. (3)) and slightly differ from the final rates in Table 3 by the collisional treatment (see text). Here (1) stands for calculations of Phillips et al. (1996); (2) for calculations of Dubernet at al. (2002) and Grosjean et al. (2003) using the same PES as in (1) and an improved treatment of scattering; and (3) for the present calculations using the new PES by Faure et al. (2005a) and the same treatment of scattering as in (2). Part b illustrates the influence of the sole PES at 5, 12 and 20 K. Conventions are similar to those of upper part.

The ratios (3)/(2) give the effect of using the new PES of Faure et al. (2005a) only, the ratios (2)/(1) show the difference induced by different dynamical calculations, while the ratios (3)/(1) include the effect of both the new PES (Faure et al. 2005a) and of the scattering calculations described elsewhere (Dubernet & Grosjean 2002; Grosjean et al. 2003). The rate coefficients for H2O-para-H2(j=0)presented in Dubernet & Grosjean (2002) were surprisingly different from those presented earlier by Phillips et al. (1996) at 20 K (in spite of using the same PES). In particular, the excitation rate coefficient for the 101-110 transition of water was more than 50% larger than the Phillips et al. (1996) result. This was thought to be a consequence of an inadequately fine energy grid used by  Phillips et al. (1996) in integrating cross-sections over a Maxwellian distribution of kinetic energies. The results of Grosjean et al. (2003) calculated between 5 K and 20 K for ortho-H2(j=1) did not present any major difference at 20 K to the Phillips et al. (1996) results. These results are summarized by ratios (2)/(1) in Table 2.

We find that variations due to the PES can be as large as a factor 4 at 5 K for the de-excitation of the 110 to 101 transition by para-H2(j=0), then decreasing to a factor 1.76 at 20 K. The (3)/(2) ratios of the transition 303 to 221 show similar behaviors with a factor of 2.74 at 5 K. For other transitions with para-H2(j=0) the (3)/(2) ratios are not sensitive to temperature and are either increased by 80-100% or decreased by about 20%, noting that changes due to the PES do not occur uniformally among transitions.

Another feature is that collisions with ortho-H2(j=1) are found to be much less sensitive to the PES than those with para-H2(j=0). Due to the absence of the dipole-quadrupole coupling term, the cross sections involving para-H2(j=0) are about an order of magnitude smaller, and are deeply affected by the strong resonance regime in the 5-20 K temperature range (Dubernet & Grosjean 2002). While the resonance regime is very sensitive to the details of the PES, strongly affecting the para-H2(j=0) cross-sections and rates, it plays a smaller role in relative value for ortho-H2(j=1), limiting the dependance of the ortho-H2(j=1) rates on the PES. Consequently, the variations due to the new PES varies between a few and 40% for collisions with ortho-H2(j=1) between 5 K and 20 K.

3 New rotational rate coefficients between 5 K and 20 K

For astrophysical applications a new set of effective de-excitation rate coefficients is now calculated using the average 5D PES of Faure et al. (2005a) and slightly improved scattering calculations compared to the procedure used in the previous sections and described in our previous publications (Dubernet & Grosjean 2002; Grosjean et al. 2003). We use close-coupling calculations for the whole energy range, the H2 energy levels are now experimental energies (Dabrowski 1984) and the basis set for H2O has been increased to $j_{\rm H_2O} =8$ with a cut-off procedure in energy in order to achieve better accuracy for the last 2 channels. This larger rotational basis set has no influence on the transitions previously considered by  Dubernet & Grosjean (2002) and Grosjean et al. (2003). This can be verified by comparing the rate coefficients set (3) of Table 2 with the final rate coefficients of Table 3. The basis sets for para/ortho-H2 still contain respectively $j(\rm H_2)=0,~2$ and $j(\rm H_2)=1,~3$. Extreme care was taken to reproduce all resonances correctly and to span the whole energy range necessary to obtain reliable rate coefficients up to 20 K.

Those effective rotational inelastic rate coefficients are given by the sum of the inelastic rate coefficients (Eq. (1)) over the final j2' states for a given initial j2:


 \begin{displaymath}\hat{R}_{j_2}(j\alpha \rightarrow j'\alpha')(T)\!=\!\sum_{j'_2}
R(j\alpha,j_2 \rightarrow j'\alpha',j'_2)(T),
\end{displaymath} (1)

and the state-to-state rotational inelastic rate coefficients are the Boltzmann thermal averages of the inelastic cross sections:

 \begin{displaymath}R(\beta \rightarrow \beta')(T) = \left({\frac{8}{\pi\mu}}\rig...
...igma _{\beta\rightarrow \beta'}(E)E
{\rm e}^{-E/kT}{\rm d}E,
\end{displaymath} (2)

where $\beta\equiv j\alpha,j_2$ and $\beta'\equiv j'\alpha',j'_2$, unprimed and primed quantum numbers label initial and final states of the molecules, j and j2 are the rotational angular momenta of H2O and H2 and $\alpha$ specifies the other H2O quantum numbers (e.g. K-1, K1), E is the kinetic energy and k is the Boltzmann constant.

In the temperature range between 5 and 20 K, our effective rate coefficients follow the principle of detailed balance within the calculation error, because there are no transitions to the j2 = 3 rotational level of ortho-H2, and transitions to the j2 =2rotational level of para-H2 are negligeable. Therefore the excitation effective rate coefficients can be obtained by detailed balance using the energy levels given by Green et al. (1993) and listed in the BASECOL database (Dubernet et al. 2006). But in general only the $R(j\alpha, j_2 \rightarrow j'\alpha', j'_2)$ and the $R(j\alpha\rightarrow j'\alpha')$ (rate coefficients summed over all final j'2 states and averaged over initial j2 states) satisfy the usual detailed balance relations between forward and reverse rate coefficients, and the $\hat{R}_{j_2}(j\alpha \rightarrow j'\alpha')$do not since they do not involve a complete thermal average (Phillips et al. 1996).

Table 3: Effective de-excitation rate coefficients for Eq. (2) (in cm3s-1) for para-H2O with para-H2 (j2=0) and ortho-H2 (j2=1). The levels are labelled with jK-1K1. The levels are listed by increasing energy.

Tables 3 presents a sub-set of these new effective rate coefficients calculated among the lowest 10 levels of ortho-H2O and para-H2O owing to collisions with para-H2(j=0) and ortho-H2(j=1) at kinetic temperatures ranging from 5 to 20 K. The full set of effective rate coefficients, state-to-state rate coefficients and energy levels used in the calculations can be retrieved from the BASECOL database (Dubernet et al. 2006) or obtained on request from one of the authors (MLD). The BASECOL database (Dubernet et al. 2006) also provides Einstein coefficients taken from either the JPL (Pickett et al. 1998) or the CDMS (Müller et al. 2005) catalogs. The labelling of transitions for the Einstein coefficients is consistant with the labelling of energy levels and rate coefficients.

4 Concluding remarks

A new average 5D potential energy surface has been used to calculate new collisional de-excitation rate coefficients for the lowest 10 levels of ortho and para-water in collision with ortho (j=1) and para-H2(j=0) for kinetic temperatures from 5 to 20 K. The new PES of Faure et al. (2005a) leads to a significant re-evaluation of the rate coefficients for the excitation of H2O by para-H2(j=0). Indeed for the 110 to 101 transition observed by the SWAS satellite, the new PES increases the de-excitation rate coefficient by a factor of 4 at 5 K and by 75% at 20 K compared to the previous results of Dubernet et al. (Dubernet & Grosjean 2002; Grosjean et al. 2003), leading to a total increase at 20 K of 185% with respect to the  Phillips et al. (1996) results. The de-excitation rate coefficients of other transitions induced by para-H2 can increase by as much as 100% or decrease by 30%. For collision with ortho-H2(j=1) the new PES has a smaller effect on the de-excitation rate coefficients: a maximum change of 40%.

Although fairly accurate, the present 5D PES may still deviate from the CCSD(T) limit by about 1-2 cm-1 and has no provision for high excitations beyond the CCSD(T) limit. Recent investigation of CO-H2 interactions by Noga et al. (2006) indicate that the missing triple and quadruple excitations lower the absolute minimum by about 3 cm-1 and change the anisotropic behaviour. However the CO molecule, with its triple bond and small dipole, might be an extreme case in this respect. The implications of such minor PES inaccuracies for the low temperature H2-H2O rates is hard to estimate.

Our next goal is to obtain the best possible accuracy for the whole temperature range up to 500 K and for all levels of interest, using close coupling and coupled states methods. Quasi classical calculations are also feasible and provide a fair agreement with quantal calculations down to 100 K (Faure et al. 2006). However they are unreliable for estimating the smallest rates, i.e. below 10-11 cm3 s-1, and quantal calculations therefore are preferred.

Higher temperature, up to 1000 K, rotational rate coefficients will also be obtained using IOS calculations checked against coupled states calculations, and will be compared with vibrational close coupling-IOS (VCC-IOS) calculations in order to test the influence of the first bending state of water.

Acknowledgements
Most scattering calculations were performed at the IDRIS-CNRS and CINES under project 2005 04 1472. This research was supported by the CNRS national program "Physique et Chimie du Milieu Interstellaire'' and by the FP6 Research Training Network "Molecular Universe'', contract Number: MRTN-CT-200'-512302. F.D. and M.W. were supported by the Ministère de l'Enseignement Supérieur et de la Recherche. F.D. wishes to thank Spanish MAE-AECI 2004/2005 for a fellowship grant from Program IIa.

References

 

Copyright ESO 2006