.
ECCC3 Paper 32

INTRODUCTION.

In general, polyatomic radical molecules such as H3., H2F., H3O., NH4., CH5. have been called Rydberg molecules by G.Herzberg [Herzberg 87] . These molecules are viewed by Herzberg to possess, in their ground state, a lone electron which orbits around a positive core, on a low-lying yet diffuse orbital.

The (H3O.) radical has a long and tumultuous career. It has been found, lost, and found again quite a number of times. This species has been first hypothetized to exist in 1907 [Goldsmith 1907] . In several experiments in the thirties, the H3O. was claimed to have been found. In 1964, Magee [Magee 64] suggested that the H3O. radical should be an important transient species in radiation chemistry, formed by the recapture process of a quasifree electron by H3O+ . Sworki [Sworski 64] , from kinetics data, suggested that, in the Gamma-radiolysis of water, H3O. should be the main precursor of H2. After an 85 eV or 100 eV electron radiolysis, in gaz-phase, Melton and coworkers [Melton 66] , [Melton 67] , claimed to find evidence of H3O.. In their experiments, it was H3O+ which was detected by mass spectrometry, after ionization of the neutral or negatively charged parent species. It was supposed that the main source of H3O+ is H3O. ionization. H3O+ might be also a direct product of radiolysis. The measured IP is 10.9 eV. Melton and coworkers hinted that this high value could correspond to the IP of a hydronium radical in an excited state. Nevertheless, Melton and coworkers were able to detect D3O- directly [Melton 67] .

Kongshaun and coworkers [Kongshaun 71] irriadiated acidic solutions with X-rays,Gamma-Rays, H. and Ar atoms generated by an electric discharge. They concluded that H3O. must have been formed.

Martin and Swift [Martin 71] claimed to have isolated the H3O. and D3O. radicals in UV-irradiated aqueous matrices, and presented an ESR quartet spectrum. However their experiment was not reproducible [Wargon 72] , and was latter explained in terms of CH3. impurities [Noda 73] .

At last, in 1984, using a neutralized ion beam spectroscopy technique, Gellene and Porter confirmed a preliminary study [Williams 80] , and brought solid experimental evidence [Gellene 84] of the metastability of D3O.. They found also D3O(D2O). and H3O(H2O)., but found no trace of H3O.. They reported an ionization potential (IP) of 4.3 eV for D3O. , and suggested a 3.4 eV IP lower limit for the hydronium hydrate radical. For comparison, the IP for a sodium atom is 5.1 eV. Latter, Raksit and Porter observed an intriguing oxygen atom isotope effect. [Raksit 87]. With its most naturally abundant isotope 16O, D3O. is metastable. When substituted with the heavier 18O, D3O. is no longer detectable while D518O 2. can still be detected. These surprising isotop effects were recently confirmed in further experiments [Hudgins 94] .

In 1987, Griffith, Harris and Beynon [Griffith 87] , using a neutralization-reionization spectroscopy technique, claimed that the metastability of H3., NH4., H3O. was established. In particular, the lifetime of H3O. was claimed to extend up to 1 microsecond. However, after reinterpretation of their own data, these authors soon retracted their claims.



Picosecond laser spectroscopy experiments were conducted in mixed water/alcohol solvents, with the help of water-specific fluorescent probe molecules, in Robinson's laboratory [Lee 85] [Moore 85] . It was found that a (H2O)4+/-1 cluster was the lowest size water acceptor structure, required both for proton hydration and "hydrated electron" formation. Latter on, Robinson disregarding the traditional hydrated electron cavity model, proposed instead a semi-ionic pair model ( H3O. -- e-)aq of the "hydrated electron", where the hydronium radical plays a major role [Robinson 87] , [Hameka 87] .

Subsequently, Muguet proposed an itinerant radical model [Muguet 92] , [Muguet 94] for the "hydrated electron". The 1.7 eV absorbing species, traditionally known as the "hydrated electron", is no longer viewed as an electron solvated in liquid water cavities, but as a hydrated hydronium radical. The hydronium radical is considered as an itinerant species, constituted with ever-changing oxygen and hydrogen atoms. It is a topological locus, moving across the H-bond network. H3O. is neutral, but because of the high polarizablity of the unpaired electron cloud, it moves like a negatively charged species, in presence of an electric field.

Within femtosecond laser spectroscopy experiments performed in Gauduel's research group [Gauduel 93] , [Gauduel 94a] , [Gauduel 94b] , ( see also our WATOC96 WWW poster ) after a two-photon photoexcitation (2 x 4 eV = 8 eV) of pure liquid water, three spectroscopic species have been discriminated :

A short-lived infrared band (" pre-hydrated electron") species.

A long-lived 1.7 eV band " hydrated electron" (e-)aq species.

A 1.45 eV centered band species. with a large kinetic isotopic H/D substitution effect:

In pure light water, the 1.45 eV band begins to appear within the pump pulse, and reaches its maximum 500 fs after and has disappeared in 2 ps. According to pseudo-first order kinetics after pump/probe signal deconvolution, the appearance time (T=1/k ) is 130 fs and the mean lifetime is 340 fs, with an error bar of +/- 20 fs.

In pure deuterated water, the 1.45eV band begins also to appear within the pump pulse, and but reaches its maximum 750 fs after and has disappeared in 4 ps. The appearance time (T=1/k) is 320 fs while the mean lifetime is 750 fs

In concentrated acidic solutions [Gauduel 90] , a near IR band in the 1.45 eV region has also been observed, which has been written conveniently as an ( H3O+ -- e-)hyd encounter pair while no hypothesis concerning this latter structure were proposed.

The hydronium radical has also been the object of numerous theoretical investigations. Conway,Bockris and Linton [Conway 55] theorized that the radical H3O. as a transient species must exist in order to explain proton conductance. Bernstein [Bernstein 63] tried to show that HA radicals might be thermodynamically stable, where A is a saturated proton acceptor. He suggested reactions with H. atoms, in order to form such radicals. With a one-center STO basis set, Bishop [Bishop 66] computed H3O. to be thermodynamically unstable. Melton and Joy [Melton 67] , also with a one-center STO expansion, found an equilateral triangle geometry for H3O. which was unstable against dissociation. At the ROHF level, Gangi and Bader [Gangi 71] found a pyramidal C3v geometry with a small inversion barrier of 2 kcal/mole, and an IP of 4.6eV. Their basis set was complemented by a single diffuse s-type gaussian centered on the oxygen atom. They predict that the H3O. radical might be marginally metastable at low temperature with a dissociation barrier of only 6.6 kcal/Mole. At the INDO level [Efskind 72] , a planar geometry was obtained, and the IP was found to be 7.45 eV. Furthermore the H3O. radical was predicted to be stable against dissociation ( -40 kcal/mole). ESR parameters were computed at the INDO and CNDO levels [Chuvylkin 72] . At the UHF level, M.C.R. Symons and coworkers [Claxton 73] found a pyramidal geometry and also found large positive hyperfine coupling constants on the protons. Magee [Magee 64] proposed a pathway for the formation of the hydrated electron involving the formation of H3O.. The hydronium radical is first formed through the capture of a quasifree electron by H3O+. This step would be followed by an autoionization of H3O. which would leads to the solvation of the ionized electron into a liquid water cavity. Webster [Webster 75] investigated the ability of the H3O+ cation to capture an extra electron. The resulting radical was computed to be slightly nonplanar. Magee's suggestion [Magee 64] also stimulated an ab initio study by Maurice Schwartz [Schwartz 76] who found that the transition to the lowest excited state corresponded to an energy of 1.87 eV. The IP is 4.75 eV. His basis set is supplemented by two s and two p diffuse gaussians on the oxygen atom, and one extra diffuse s-type and p-type gaussians on each hydrogen atom. Strangely enough, this interesting paper is never quoted.

Niblaeus ,Roos and Siegbahn [Niblaeus 77] performed several computations both at the SCF/UHF and CI levels. Their largest basis set comprises a valence gaussian set with polarization functions, supplemented by only two s diffuse gaussians centered on the oxygen atom. H3O. was determined to be thermodynamically unstable, the dissociation products being H2O and H.. The H3O. radical corresponds to a local minimum, one pathway to the formation of H 3O , being electron attachment to H3O+. It is interesting to mention absolute energies since it gives an idea of the quality of the basis set and the method. Their best energies for H3O. were, at the UHF level -76.499,800 Hartrees and with a direct CI (10,074 CSFs) method -76.715,324 Hartrees. These authors computed a dissociation barrier of 4.6 kcal/Mole at the UHF level and a dissociation barrier of only 3.4 kcal/Mole at the CI level. They concluded that the existence of a metastable state with a measurable lifetime was unlikeky.

Raynor and Herschbach [Raynor 82] computed the transition energies as well as the transition dipole moments of electronic transitions between the ground state and excited states as well as between excited states for various radicals. Each electronic state was calculated with respect to a frozen cationic core. Their cation basis set is supplemented by five s, five p,five d diffuse STOs centered on the oxygen atom. For H3O., they predicted two transitions at 550 nm (2.25 eV) and 661 nm (1.8 eV) with large Einstein emission coefficients 6.48 107 s-1 and 5.64 107 s-1 respectively, to the ground state.

Talbi and Saxon [Talbi 89] performed a thorough study of the PES of the ground state and the first four doublet excited states. Their basis set is supplemented by two s and two p diffuse gaussians centered only on the oxygen atom. In all their MCSCF and CI computations, the reference configuration carries the most important weight (0.94-0.95). Their ground state geometry is pyramidal, with an energy of -76.556,186 Hartrees at the MCSCF level, comprising up to 17000 configurations state functions (CSFs), and of -76.727,647 Hartrees at the MRSDCI level (up to 118,000 CSFs). An Einstein coefficient of 5.46 107 s-1 was calculated for the 3s --> 3p transition. The radiative lifetime of the 3p excited state was estimated around 18 ns. The ionization potential (IP) of the 3s ground state is 5.27 eV, while the IP of the 3p excited state is 2.92 eV. The dissociation barrier for H3O. ---> H2O + H. was found to be 3.58 kcal/mole. Taking into account the zero-point vibrational energies at the local minimum and at the transition point, the barrier is reduced to 0.40 kcal/mole ( D3O 1.32 kcal/mole). Talbi and Saxon could not explain why D3O. might be metastable, contrary to H3O.. Another question which was left unanswered is why D318O. is not detected in experiments, although it includes a heavier oxygen isotope 18O nucleus which implies a slightly lower zero-point vibrational energy. After establishing a correlation diagram, Talbi and Saxon concluded that the dissociation barrier does not arise from configuration mixing or avoided crossing, but originates from a de-Rydbergization process [Mulliken 77] .

McLoughin and Gellene [McLoughin 92] took into account vibrational energy and, rather interestingly, tried to explain hydronium radical metastability with the help of the concept of "adiabatic molecules" as set forward by Pollak [Pollak 83] . A dynamical barrier contribution might appear, if the vibrational energy increases along the dissociation pathway. Towards this goal, these authors computed the energy at more than 300 different geometries at the MP3 level, in order to establish a pointwise harmonic fit of the potential energy surface (PES) in the vicinity of the dissociation pathway. Their valence basis set is supplemented by three s and two p diffuse gaussians on the oxygen atom, and three s-type and one p-type diffuse gaussians on each hydrogen atom. In the vibrational ground state, the estimated tunneling lifetimes in the order of 10-13,10-12 s are inconsistent with experimental lifetime. A dynamical barrier, high enough to permit a microsecond lifetime, is predicted to appear only when the asymmetric strech normal mode ( 1444 cm-1 ) reaches its fourth excited state. In this unlikely high overtone, H3O. is predicted to feature a lifetime at least four times larger than D3O., in contradiction with experimental data. In regards to oxygen isotop effects, attention was brought on the sensitivity of one cubic and one quartic potential constants to oxygen isotopic substitution.

It is quite important to mention that MRCI calculations predict that the hydronium radical should be formed after excitation of a linear water dimer. The interplay between the various excited states appears rather complicated [Sosa 93] .

Recently, at a fixed geometry, within a dielectric cavity, we have performed a series of computations of MCSCF-RF (Reaction Field) ground state H3O. wavefunctions. Then, starting from these wavefunctions, using an electron propagator method ( RPA approximation ), we have been able to compute the excitation spectrum of (H3O.)hyd. [Muguet 96a] , [Muguet 96b] , ( see also our WATOC96 WWW poster ). We predicted that the hydrated (H3O.)hyd radical should feature an absorption band, which is, in regard to both shape and energy range, quite similar to the experimentally observed "hydrated electron" absorption band. This computational result brings support to Muguet's itinerant radical model. Our results suggest also that the excited state radical (H3O.)hyd* could be assigned to the 1.45 eV absorbing species.

Very recently, at a fixed geometry, within a dielectric cavity, we have performed a series of computations of MCSCF-RF (Reaction Field) ground state H3O. wavefunctions. Then, starting from these wavefunctions, using an electron propagator method ( RPA approximation ), we have been able to compute the excitation spectrum of (H3O.)hyd. [Muguet 96a] , [Muguet 96b] , ( see also our WATOC96 WWW poster ). We found that the predicted hydrated (H3O.)hyd radical and experimental (e-)aq absorption bands feature similar shapes, within the same spectral 1.7 eV (720nm) range. Our results suggest also that the excited state radical (H3O.)hyd* could be assigned to the 1.45 eV absorbing species.

Therefore, it appears that, as Talbi and Saxon concluded, there are still "substantial gaps in our understanding of the H3O. system". both in the gas phase and in the liquid phase.



We feel that the Rydberg radical denomination brings in itself the seed of possible misconceptions.

1/ Rydberg excited states are hydrogen-like states [Robin 75] which feature an atomic-like sharp spectra.

2/Highly excited Rydberg states have a quasi-classical character, and the extra electron may be conceived as circling around the positive core like in a planetary orbit.

Low-lying Rydberg ground states are quite different from highly excited Rydberg states.

In this paper, several questions are of a more specific concern :
What is the exact nature of the so called low lying Rydberg orbital ?.
Can we intuitively understand and visualize the deRydbergisation process involved in the dissociation process ?.
Are the basis sets employed so far, large enough ?.


Top page Last page Next page