Threshold-photoelectron-spectroscopy-coincidence study of the double photoionization of HBr

A threshold-photoelectron-coincidence spectrum of HBr has been recorded in the 32.2-35.8 eV photon energy range, with a resolution of approximately 0.01 eV, using a synchrotron radiation source. The X (3)Sigma(-) and a (1)Delta(2) states of the HBr(2+) dication are clearly observed in the spectrum, while there is no clear evidence for the formation of the b (1)Sigma(+) electronic state. For the first two states, the vibrational states v=0-3 have been resolved, while for the ground X (3)Sigma(-) state also spin-orbit splitting has been detected. The results appear in good agreement with previous experimental observations. A comparison with theoretical predictions indicates the role of "noncovalent" contributions to the interaction between the two atomic partners for the formation of metastable states.     

A theoretical and experimental study of the SO2(2+) dication

The double photoionization spectrum of SO2 has been measured using the TOF-PEPECO technique and contains one resolved band. Detailed electronic structure calculations and experimental comparisons allow the resolved band to be identified as the A 1A2 state of the SO2(2+) dication, with its adiabatic ionization energy at 35.284+/-0.02 eV. According to the most accurate calculations, the ground state level of SO2(2+) must be located near 33.48 eV, well below the range accessed by vertical transitions from neutral SO2. Transient SO2 (2+) molecules detected by mass spectrometry may be identified either as the sharp levels of the A 1A2 state or as ground state levels populated by nonvertical ionization pathways.     

Ab initio study of the low lying electronic states of ZnF and ZnF-

Highly correlated ab initio calculations have been performed for an accurate determination of the electronic structure and of the spectroscopy of the low lying electronic states of the ZnF system. Using effective core pseudopotentials and aug-cc-pVQZ basis sets for both atoms, the potential curves, the dipole moment functions, and the transition dipole moments between relevant electronic states have been calculated at the multireference-configuration-interaction level. The spectroscopic constants calculated for the X(2)Sigma(+) ground state are in good agreement with the most recent theoretical and experimental values. It is shown that, besides the X(2)Sigma(+) ground state, the B(2)Sigma(+), the C(2)Pi, and the D(2)Sigma(+) states are bound. The A(2)Pi state, which has been mentioned in previous works, is not bound but its potential presents a shoulder in the Franck-Condon region of the X(2)Sigma(+) ground state. All of the low lying quartet states are found to be repulsive. The absorption transitions from the v=0 level of the X(2)Sigma(+) ground state toward the three bound states have been evaluated and the spectra are presented. The potential energy of the ZnF(-) molecular anion has been determined in the vicinity of its equilibrium geometry and the electronic affinity of ZnF (EA=1.843 eV with the zero energy point correction) has been calculated in agreement with the photoelectron spectroscopy experiments.     

Theoretical study of electronic states of N22+in an intense radiation field

The Floquet states of N(2) (2+) created by the interactions of the six lowest singlet (1 (1)Sigma(g) (+), 1 (1)Delta(g), 2 (1)Sigma(g) (+), 1 (1)Pi(u), 1 (1)Pi(g), and 1 (1)Sigma(u) (-)) states of the dication with intense (0.4 x 10(13) Wcm(2)) radiation have been studied using the recently developed multireference configuration interaction method with single and double excitations (MRCISD)-based approach. The adiabatic Floquet state coinciding near its minimum with the initial X (1)Sigma(g) (+) ground state and asymptotically correlating with A (1)Pi(u) (m = -1), i.e., with one less photon in the dressed state, is expected to be metastable, as is the ground state in the absence of a field, but to support up to the v(max) = 12 quasibound vibrational level in comparison with v(max) = 11 in the parent field-free X (1)Sigma(g) (+) ground state. The tunneling lifetimes of the highest vibrational levels in this adiabatic Floquet state are predicted to be several orders longer than those in the parent field-free state. Analysis of the complete basis set limit extrapolated MRCISD potential energy curve of the field-free X (1)Sigma(g) (+) state of N(2) (2+) calculated in the present work (R(e) = 1.130 A, omega(e) = 2011 cm(-1), omega(e)x(e) = 26.1 cm(-1)) is in good agreement with spectroscopic experimental data. Calculations on the field-free A (1)Pi(u) state (T(e) = 12 106 cm(-1), R(e) = 1.252 A, omega(e) = 1438 cm(-1), omega(e)x(e) = 23.5 cm(-1)) generally support earlier theoretical work and do not support reported experimental values.     

A theoretical study of the fine and hyperfine interactions in the NCO and CNO radicals

The geometries, the harmonic vibrational frequencies, and the Renner-Teller parameter have been reported for the NCO(+)(X (3)Sigma(-)), NCO(X (2)Pi,A (2)Sigma(+),B (2)Pi,2 (2)Sigma(+)), NCO(-)(X (1)Sigma(+)), CNO(+)(X), CNO(X (2)Pi,A (2)Sigma(+),B (2)Pi,2 (2)Sigma(+)), and CNO(-)(X (1)Sigma(+)) systems at the full valence-complete active space self-consistent-field (fv-CASSCF) level of theory. The (2)Pi electronic states of the NCO and CNO radicals have two distinct real vibrational frequencies for the bending modes and these states are subject to the type A Renner-Teller effect. The total energy of CNO(+) without zero point energy correction of the linear geometry is approximately 31 cm(-1) higher than the bent geometry at the fv-CASSCF level and the inversion barrier vanishes after the zero point energy correction; therefore, the ground state of the CNO(+) may possess a quasilinear geometry. The spin-orbit coupling constants estimated using atomic mean field Hamiltonian at the fv-CASSCF level of theory are in better agreement with the experimental values. The excitation energies, the electron affinity, and the ionization potential have been computed at the complete active space second order perturbation theory (CASPT2) and the multireference singles and doubles configuration (MRSD-CI) levels of theory. The computed values of the electric hyperfine coupling constants for the (14)N atom in the ground state of the NCO radical agree well with the experimental data. The magnetic hyperfine coupling constants (HFCC's) have been estimated employing the configuration selected MRSD-CI and the multireference singles configuration interaction (MRS-CI) methods using iterative natural orbitals (ino) as one particle basis. Sufficiently accurate value of the isotropic contribution to the HFCC's can be obtained using an MRS-CI-ino procedure.     

Ground state potential energy curve and dissociation energy of MgH

New high-resolution visible emission spectra of the MgH molecule have been recorded with high signal-to-noise ratios using a Fourier transform spectrometer. Many bands of the A 2Pi-->X 2Sigma+ and B' 2Sigma+-->X 2Sigma+ electronic transitions of 24MgH were analyzed; the new data span the v' = 0-3 levels of the A 2Pi and B'2Sigma+ excited states and the v'=0-11 levels of the X 2Sigma+ ground electronic state. The vibration-rotation energy levels of the perturbed A 2Pi and B' 2Sigma+ states were fitted as individual term values, while those of the X 2Sigma+ ground state were fitted using the direct-potential-fit approach. A new analytic potential energy function that imposes the theoretically correct attractive potential at long-range, and a radial Hamiltonian that includes the spin-rotation interaction were employed, and a significantly improved value for the ground state dissociation energy of MgH was obtained. The v'=11 level of the X 2Sigma+ ground electronic state was found to be the highest bound vibrational level of 24MgH, lying only about 13 cm(-1) below the dissociation asymptote. The equilibrium dissociation energy for the X 2Sigma+ ground state of 24MgH has been determined to be De=11104.7+/-0.5 cm(-1) (1.37681+/-0.00006 eV), whereas the zero-point energy (v'=0) is 739.11+/-0.01 cm(-1). The zero-point dissociation energy is therefore D0=10365.6+/-0.5 cm(-1) (1.28517+/-0.00006 eV). The uncertainty in the new experimental dissociation energy of MgH is more than 2 orders of magnitude smaller than that for the best value available in the literature. MgH is now the only hydride molecule other than H2 itself for which all bound vibrational levels of the ground electronic state are observed experimentally and for which the dissociation energy is determined with subwavenumber accuracy.     

Electronic states and spectroscopic properties of SiTe and SiTe+

Ab initio based configuration interaction calculations have been carried out to study the low-lying electronic states and spectroscopic properties of the heaviest nonradioactive silicon chalcogenide molecule and its monopositive ion. Spectroscopic constants and potential energy curves of states of both SiTe and SiTe+ within 5 eV are reported. The calculated dissociation energies of SiTe and SiTe+ are 4.41 and 3.52 eV, respectively. Effects of the spin-orbit coupling on the electronic spectrum of both the species are studied in detail. The spin-orbit splitting between the two components of the ground state of SiTe+ is estimated to be 1880 cm(-1). Transitions such as 0+ (II)-X1Sigma(+)0+, 0+ (III)-X1Sigma(+)0+, E1Sigma(+)0+ -X1Sigma(+)0+, and A1Pi1-X1Sigma(+)0+ are predicted to be strong in SiTe. The radiative lifetime of the A1Pi state is less than a microsecond. The X(2)2Pi(1/2)-X(1)2Pi(3/2) transition in SiTe+ is allowed due to spin-orbit mixing. However, it is weak in intensity with a partial lifetime for the X2 state of about 108 ms. The electric dipole moments of both SiTe and SiTe+ in their low-lying states are calculated. The vertical ionization energies for the ionization of the ground-state SiTe to different ionic states are also reported.     

Photodissociation spectroscopy of stored CH+ and CD+ ions: analysis of the b 3Sigma(-)-a 3Pi system

We have measured the photodissociation spectrum of CH(+) and CD(+) molecular ions, stored as fast (MeV) ion beams in the heavy-ion storage ring TSR. Several b (3)Sigma(-)-a (3)Pi bands were observed as strong resonances because a large fraction of the ions in the metastable a (3)Pi(v=0) state were pumped to b (3)Sigma(-) levels and predissociated via the c (3)Sigma(+) state into C(+) and H(D) fragments. From a rotational analysis of the 2-0, 3-0, and 4-0 bands in CH(+) and the 3-0 and 4-0 bands in CD(+), we derive spectroscopic constants for these levels and also revise a previous analysis of the 0-0 and 1-0 bands in CH(+). Combining all data delivers new, significantly adjusted equilibrium constants for the b (3)Sigma(-) and a (3)Pi electronic states. Apart from the spectroscopic analysis, we estimate the predissociation rates of the upper b (3)Sigma(-) vibrational levels in CH(+) and compare them to a model. For the initial rovibrational distribution of the stored metastable CH(+) molecules, the data indicate a faster vibrational cooling than derived before, and rotational cooling at a rate similar to the X (1)Sigma(+) ground state. New aspects of the spin-forbidden a (3)Pi-X (1)Sigma(+) radiative decay are discussed. Finally, we predict b (3)Sigma(-)-a (3)Pi absorption and a (3)Pi-X (1)Sigma(+) emission lines through which CH(+) in the metastable a (3)Pi(v=0) state might be detectable in astrophysical environments.     

Experimental and theoretical study of the electronic spectrum of BeAl

The electronic structure of BeAl was investigated by laser induced fluorescence and resonance enhanced multiphoton ionization spectroscopy. BeAl was formed by pulsed laser ablation of a Be/Al alloy in the presence of helium carrier gas, followed by a free jet expansion into vacuum. In agreement with recent ab initio studies, the molecule was found to have a (2)Pi(1/2) ground state. Transitions to two low lying electronic states, (2)(2)Pi(1/2)(v') <-- X (2)Pi(1/2) (v' = 0) and (1)(2)Delta(v') <-- X (2)Pi(1/2) (v' = 0,1), were observed and rotationally analyzed. An additional band system, identified as (4)(2)Sigma(+)(v') <-- X (2)Pi(1/2), was found in the 28 000-30 100 cm(-1) energy range. This transition exhibited an unusual pattern of vibrational levels resulting from an avoided crossing with the (5)(2)Sigma(+) electronic state. New multi-reference configuration interaction calculations were carried out to facilitate the interpretation of the UV bands.An ionization energy of 48 124(80) cm(-1) was determined for BeAl from photoionization efficiency (PIE) measurements. Fine structure in the PIE curve was attributed to resonances with Rydberg series correlating with vibrationally excited states of the BeAl(+) ion. Analysis of this structure yielded a vibrational frequency of 240(20) cm(-1) for the cation.     

Spectroscopy and metastability of CO2 2+ molecular ions

The spectroscopy and metastability of the carbon dioxide doubly charged ion, the CO(2) (2+) dication, have been studied with photoionization experiments: time-of-flight photoelectron photoelectron coincidence (TOF-PEPECO), threshold photoelectrons coincidence (TPEsCO), and threshold photoelectrons and ion coincidence (TPEsCO ion coincidence) spectroscopies. Vibrational structure is observed in TOF-PEPECO and TPEsCO spectra of the ground and first two excited states. The vibrational structure is dominated by the symmetric stretch except in the TPEsCO spectrum of the ground state where an antisymmetric stretch progression is observed. All three vibrational frequencies are deduced for the ground state and symmetric stretch and bending frequencies are deduced for the first two excited states. Some vibrational structure of higher electronic states is also observed. The threshold for double ionization of carbon dioxide is reported as 37.340+/-0.010 eV. The fragmentation of energy selected CO(2) (2+) ions has been investigated with TPEsCO ion coincidence spectroscopy. A band of metastable states from approximately 38.7 to approximately 41 eV above the ground state of neutral CO(2) has been observed in the experimental time window of approximately 0.1-2.3 mus with a tendency towards shorter lifetimes at higher energies. It is proposed that the metastability is due to slow spin forbidden conversion from bound excited singlet states to unbound continuum states of the triplet ground state. Another result of this investigation is the observation of CO(+)+O(+) formation in indirect dissociative double photoionization below the threshold for formation of CO(2) (2+). The threshold for CO(+)+O(+) formation is found to be 35.56+/-0.10 eV or lower, which is more than 2 eV lower than previous measurements.     

A new mechanism for the production of highly vibrationally excited OH in the mesosphere: an ab initio study of the reactions of O2(A 3Sigmau+ and A' 3Deltau)+H

In an attempt to explain the observed nightglow emission from OH(v=10) in the mesosphere that has the energy greater than the exothermicity of the H+O(3) reaction, potential energy surfaces were calculated for reactions of high lying electronic states of O(2)(A (3)Sigma(u) (+) and A' (3)Delta(u)) with atomic hydrogen H((2)S) to produce the ground state products OH((2)Pi)+O((3)P). From collinear two-dimensional scans, several adiabatic and nonadiabatic pathways have been identified. Multiconfigurational single and double excitation configuration interaction calculations show that the adiabatic pathways on a (4)Delta potential surface from O(2)(A' (3)Delta)+H and a (4)Sigma(+) potential surface from O(2)(A (3)Sigma(u) (+))+H are the most favorable, with the zero-point corrected barrier heights of as low as 0.191 and 0.182 eV, respectively, and the reactions are fast. The transition states for these pathways are collinear and early, and the reaction coordinate suggests that the potential energy release of ca. 3.8 eV (larger than the energy required to excite OH to v=10) is likely to favor high vibrational excitation.     

Spectroscopy of the ground and low-lying excited states of ThO+

The ThO(+) cation is of interest as it is a useful prototype for experimental and theoretical studies of bonding in a simple actinide compound. Formally the ground state of ThO(+) has the configuration Th(3+)(7s)O(2-), where there is a single unpaired electron associated with a closed-shell Th(4+)-ion core. The first tier of excited states above the X (2)Sigma(+) ground state is expected to be 1 (2)Delta, 1 (2)Pi, and 2 (2)Sigma(+) derived from the Th(3+)(6d)O(2-) configuration. Spectroscopic observations of ThO(+) using the pulsed field ionization-zero kinetic-energy photoelectron technique are reported here. Rotationally resolved spectra were recorded for the X (2)Sigma(+), 1 (2)Delta, and 1 (2)Pi states. Extensive vibrational progressions were observed. Surprisingly, it was found that ionization of ThO decreases the dissociation energy, while increasing the vibrational frequency and decreasing the bond length. Accurate values for the ionization energies of ThO [53 253.8(2) cm(-1)] and Th [50 868.71(8) cm(-1)] were determined as part of this investigation.     

Vacuum ultraviolet photoionization of C3

Photoionization efficiency (PIE) curves for C(3) molecules produced by laser ablation are measured from 11.0 to 13.5 eV with tunable vacuum ultraviolet undulator radiation. A step in the PIE curve versus photon energy, obtained with N(2) as the carrier gas, supports the conclusion of very effective cooling of C(3) to its linear (1)Sigma(g)(+) ground state. The second step observed in the PIE curve versus photon energy could be the first experimental evidence of the C(3)(+)((2)Sigma(g)(+)) excited state. The experimental results, complemented by ab initio calculations, suggest a state-to-state vertical ionization energy of 11.70 +/- 0.05 eV between the C(3)(X(1)Sigma(g)(+)) and the C(3)(+)(X(2)Sigma(u)(+)) states. An ionization energy of 11.61 +/- 0.07 eV between the neutral and ionic ground states of C(3) is deduced using the data together with our calculations. Accurate ab initio calculations are performed for both linear and bent geometries on the lowest doublet electronic states of C(3)(+) using Configuration Interaction (CI) approaches and large basis sets. These calculations confirm that C(3)(+) is bent in its electronic ground state, which is separated by a small potential barrier from the (2)Sigma(u)(+) minimum. The gradual increase at the onset of the PIE curve suggests a geometry change between the ground neutral and cationic states. The energies between several doublet states of the ion are theoretically determined to be 0.81, 1.49, and 1.98 eV between the (2)Sigma(u)(+) and the (2)Sigma(g)(+),( 2)Pi(u), (2)Pi(g) excited states of C(3)(+), respectively.     

Intruder states in multireference perturbation theory: the ground state of manganese dimer

A detailed analysis of a severe intruder state problem in the multistate multireference perturbation theory (MS-MRPT) calculations on the ground state of manganese dimer is presented. An enormous number of detected intruder states (> 5000) do not permit finding even an approximate shape of the X(1)Sigma(g) (+) potential energy curve. The intruder states are explicitly demonstrated to originate from quasidegeneracies in the zeroth-order Hamiltonian spectrum. The electronic configurations responsible for appearance of the quasidegeneracies are identified as single and double excitations from the active orbitals to the external orbitals. It is shown that the quasidegeneracy problem can be completely eliminated using shift techniques despite of its severity. The resultant curves are smooth and continuous. Unfortunately, strong dependence of the spectroscopic parameters of the X(1)Sigma(g) (+) state on the shift parameter is observed. This finding rises serious controversies regarding validity of employing shift techniques for solving the intruder state problem in MS-MRPT. Various alternative approaches of removing intruder states (e.g., modification of the basis set or changing the active space) are tested. None of these conventional techniques is able to fully avoid the quasidegeneracies. We believe that the MS-MRPT calculations on the three lowest A(g) states of manganese dimer constitute a perfect benchmark case for studying the behavior of MRPT in extreme situations.     

Reaction dynamics of Si(3PJ)+O2-->SiO(X 1 Sigma+)+O studied by a crossed-beam laser-induced fluorescence technique

Oxidation reaction of the ground state Si atom was studied by using a crossed molecular beam technique at 13.0 kJ/mol of collision energy. The Si atomic beam was generated by laser vaporization and crossed with the oxygen molecular beam at right angle. Products at the crossing region were detected by the laser-induced fluorescence (LIF). The LIF of SiO(A 1 Pi-X 1 Sigma+) was used to determine the vibrational state distribution of the electronic ground state, SiO(X 1 Sigma+). The determined distribution was inverted with the maximum population at v"=4, and in good agreement with the recent quasiclassical trajectory calculation on the singlet potential energy surface. The agreement suggested that an abstraction mechanism is dominant at the collision energy studied here. One of the counterproducts, O(3PJ), was also observed by the vacuum ultraviolet LIF and the distribution of the spin-orbit levels were determined. The formation of O(3PJ) was consistent with the significant population of v"=7 and 8 states of SiO, which could be explained by the presence of the triplet product channel with higher exothermicity.     

Ab initio study on the ground and low-lying excited states of cesium iodide (CsI)

Potential energy curves (PECs) for the ground and low-lying excited states of the cesium iodide (CsI) molecule have been calculated using the internally contracted multireference configuration interaction calculation with single and double excitation method with the relativistic pseudopotentials. PECs for seven Lambda-S states, X 1Sigma+, 2 1Sigma+, 3Sigma+, 1Pi, and 3Pi are first calculated and then those for 13 Omega states are obtained by diagonalizing the matrix of the electronic Hamiltonian H(el) plus the effective one-electron spin-orbit (SO) Hamiltonian H(SO). Spectroscopic constants for the calculated ground X 0+-state PEC with the Davidson correction are found to agree well with the experiment. Transition dipole moments (TDMs) between X 0 and the other Omega states are also obtained and the TDM between X 0+ and A 0+ is predicted to be the largest and that between X 0+ and B 0+ is the second largest around the equilibrium internuclear distance. The TDMs between X 0+ and the Omega=1 states are estimated to be nonzero, but they are notably small as compared with those between the 0+ states. Finally, vibrational levels of the X 0+ PEC for the two isotopic analogs, (133)CsI and (135)CsI, are numerically obtained to investigate the isotope effect on the vibrational-level shift. It has been found that the maximized available isotope shift is approximately 30 cm(-1) around nu=136.     

Theoretical study of isomerization and dissociation of acetylene dication in the ground and excited electronic states

Ab initio calculations employing the configuration interaction method including Davidson's corrections for quadruple excitations have been carried out to unravel the dissociation mechanism of acetylene dication in various electronic states and to elucidate ultrafast acetylene-vinylidene isomerization recently observed experimentally. Both in the ground triplet and the lowest singlet electronic states of C2H2(2+) the proton migration barrier is shown to remain high, in the range of 50 kcal/mol. On the other hand, the barrier in the excited 2 3A" and 1 3A' states decreases to about 15 and 34 kcal/mol, respectively, indicating that the ultrafast proton migration is possible in these states, especially, in 2 3A", even at relatively low available vibrational energies. Rice-Ramsperger-Kassel-Marcus calculations of individual reaction-rate constants and product branching ratios indicate that if C2H(2)2+ dissociates from the ground triplet state, the major reaction products should be CCH+(3Sigma-)+H+ followed by CH+(3Pi)+CH+(1Sigma+) and with a minor contribution (approximately 1%) of C2H+(2A1)+C+(2P). In the lowest singlet state, C2H+(2A1)+C+(2P) are the major dissociation products at low available energies when the other channels are closed, whereas at Eint>5 eV, the CCH+(1A')+H+ products have the largest branching ratio, up to 70% and higher, that of CH+(1Sigma+)+CH+(1Sigma+) is in the range of 25%-27%, and the yield of C2H++C+ is only 2%-3%. The calculated product branching ratios at Eint approximately 17 eV are in qualitative agreement with the available experimental data. The appearance thresholds calculated for the CCH++H+, CH++CH+, and C2H++C+ products are 34.25, 35.12, and 34.55 eV. The results of calculations in the presence of strong electric field show that the field can make the vinylidene isomer unstable and the proton elimination spontaneous, but is unlikely to significantly reduce the barrier for the acetylene-vinylidene isomerization and to render the acetylene configuration unstable or metastable with respect to proton migration.     

New view of the ICN A continuum using photoelectron spectroscopy of ICN-

Negative-ion photoelectron spectroscopy of ICN(-) (X??(2)Σ(+)) reveals transitions to the ground electronic state (X??(1)Σ(+)) of ICN as well as the first five excited states ((3)Π(2), (3)Π(1), Π(0(-) ) (3), Π(0(+) ) (3), and (1)Π(1)) that make up the ICN A continuum. By starting from the equilibrium geometry of the anion, photoelectron spectroscopy characterizes the electronic structure of ICN at an elongated I-C bond length of 2.65 A?. Because of this bond elongation, the lowest three excited states of ICN ((3)Π(2), (3)Π(1), and Π(0(-) ) (3)) are resolved for the first time in the photoelectron spectrum. In addition, the spectrum has a structured peak that arises from the frequently studied conical intersection between the Π(0(+) ) (3) and (1)Π(1) states. The assignment of the spectrum is aided by MR-SO-CISD calculations of the potential energy surfaces for the anion and neutral ICN electronic states, along with calculations of the vibrational levels supported by these states. Through thermochemical cycles involving spectrally narrow transitions to the excited states of ICN, we determine the electron affinity, EA(ICN), to be 1.34(5) (+0.04∕-0.02) eV and the anion dissociation energy, D(0)(X??(2)Σ(+) I-CN(-)), to be 0.83 (+0.04/-0.02) eV.     

Non-Born-Oppenheimer variational calculation of the ground-state vibrational spectrum of LiH+

Very accurate, rigorous, variational, non-Born-Oppenheimer (non-BO) calculations have been performed for the fully symmetric, bound states of the LiH(+) ion. These states correspond to the ground and excited vibrational states of LiH(+) in the ground (2)Sigma(+) electronic state. The non-BO wave functions of the states have been expanded in terms of spherical N-particle explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance and 5600 Gaussians were used for each state. The calculations that, to our knowledge, are the most accurate ever performed for a diatomic system with three electrons have yielded six bound states. Average interparticle distances and nucleus-nucleus correlation function plots are presented.