A kinetic study of the reactions of Ca+ ions with O3, O2, N2, CO2 and H2O

The reactions between Ca(+)(4(2)S(1/2)) and O(3), O(2), N(2), CO(2) and H(2)O were studied using two techniques: the pulsed laser photo-dissociation at 193 nm of an organo-calcium vapour, followed by time-resolved laser-induced fluorescence spectroscopy of Ca(+) at 393.37 nm (Ca(+)(4(2)P(3/2)-4(2)S(1/2))); and the pulsed laser ablation at 532 nm of a calcite target in a fast flow tube, followed by mass spectrometric detection of Ca(+). The rate coefficient for the reaction with O(3) is essentially independent of temperature, k(189-312 K) = (3.9 +/- 1.2) x 10(-10) cm(3) molecule(-1) s(-1), and is about 35% of the Langevin capture frequency. One reason for this is that there is a lack of correlation between the reactant and product potential energy surfaces for near coplanar collisions. The recombination reactions of Ca(+) with O(2), CO(2) and H(2)O were found to be in the fall-off region over the experimental pressure range (1-80 Torr). The data were fitted by RRKM theory combined with quantum calculations on CaO(2)(+), Ca(+).CO(2) and Ca(+).H(2)O, yielding the following results with He as third body when extrapolated from 10(-3)-10(3) Torr and a temperature range of 100-1500 K. For Ca(+) + O(2): log(10)(k(rec,0)/cm(6) molecule(-2) s(-1)) = -26.16 - 1.113log(10)T- 0.056log(10)(2)T, k(rec,infinity) = 1.4 x 10(-10) cm(3) molecule(-1) s(-1), F(c) = 0.56. For Ca(+) + CO(2): log(10)(k(rec,0)/ cm(6) molecule(-2) s(-1)) = -27.94 + 2.204log(10)T- 1.124log(10)(2)T, k(rec,infinity) = 3.5 x 10(-11) cm(3) molecule(-1) s(-1), F(c) = 0.60. For Ca(+) + H(2)O: log(10)(k(rec,0)/ cm(6) molecule(-2) s(-1)) = -23.88 - 1.823log(10)T- 0.063log(10)(2)T, k(rec,infinity) = 7.3 x 10(-11)exp(830 J mol(-1)/RT) cm(3) molecule(-1) s(-1), F(c) = 0.50 (F(c) is the broadening factor). A classical trajectory analysis of the Ca(+) + CO(2) reaction is then used to investigate the small high pressure limiting rate coefficient, which is significantly below the Langevin capture frequency. Finally, the implications of these results for calcium chemistry in the mesosphere are discussed.     

Experimental investigation of the complete inner shell hydration energies of Ca2+: threshold collision-induced dissociation of Ca(2+)(H2O)x Complexes (x = 2-8)

The sequential bond energies of Ca(2+)(H(2)O)(x) complexes, where x = 1-8, are measured by threshold collision-induced dissociation (TCID) in a guided ion beam tandem mass spectrometer. From an electrospray ionization source that produces an initial distribution of Ca(2+)(H(2)O)(x) complexes where x = 6-8, complexes down to x = 2 are formed using an in-source fragmentation technique. Ca(2+)(H(2)O) cannot be formed in this source because charge separation into CaOH(+) and H(3)O(+) is a lower energy pathway than simple water loss from Ca(2+)(H(2)O)(2). The kinetic energy dependent cross sections for dissociation of Ca(2+)(H(2)O)(x) complexes, where x = 2-9, are examined over a wide energy range to monitor all dissociation products and are modeled to obtain 0 and 298 K binding energies. Analysis of both primary and secondary water molecule losses from each sized complex provides thermochemistry for the sequential hydration energies of Ca(2+) for x = 1-8 and the first experimental values for x = 1-4. Additionally, the thermodynamic onsets leading to the charge separation products from Ca(2+)(H(2)O)(2) and Ca(2+)(H(2)O)(3) are determined for the first time. Our experimental results for x = 1-6 agree well with previously calculated binding enthalpies as well as quantum chemical calculations performed here. Agreement for x = 1 is improved when the basis set on calcium includes core correlation.     

Unimolecular reactions of dihydrated alkaline earth metal dications M2+(H2O)2, M = Be, Mg, Ca, Sr, and Ba: salt-bridge mechanism in the proton-transfer reaction M2+(H2O)2 --> MOH+ + H3O

The unimolecular reactivity of M(2+)(H(2)O)(2), M = Be, Mg, Ca, Sr, and Ba, is investigated by density functional theory. Dissociation of the complex occurs either by proton transfer to form singly charged metal hydroxide, MOH(+), and protonated water, H(3)O(+), or by loss of water to form M(2+)(H(2)O) and H(2)O. Charge transfer from water to the metal forming H(2)O(+) and M(+)(H(2)O) is not favorable for any of the metal complexes. The relative energetics of these processes are dominated by the metal dication size. Formation of MOH(+) proceeds first by one water ligand moving to the second solvation shell followed by proton transfer to this second-shell water molecule and subsequent Coulomb explosion. These hydroxide formation reactions are exothermic with activation energies that are comparable to the water binding energy for the larger metals. This results in a competition between proton transfer and loss of a water molecule. The arrangement with one water ligand in the second solvation shell is a local minimum on the potential energy surface for all metals except Be. The two transition states separating this intermediate from the reactant and the products are identified. The second transition state determines the height of the activation barrier and corresponds to a M(2+)-OH(-)-H(3)O(+) "salt-bridge" structure. The computed B3LYP energy of this structure can be quantitatively reproduced by a simple ionic model in which Lewis charges are localized on individual atoms. This salt-bridge arrangement lowers the activation energy of the proton-transfer reaction by providing a loophole on the potential energy surface for the escape of H(3)O(+). Similar salt-bridge mechanisms may be involved in a number of proton-transfer reactions in small solvated metal ion complexes, as well as in other ionic reactions.     

A kinetic study of Ca-containing ions reacting with O, O2, CO2 and H2O: implications for calcium ion chemistry in the upper atmosphere

A series of gas-phase reactions involving molecular Ca-containing ions was studied by the pulsed laser ablation of a calcite target to produce Ca(+) in a fast flow of He, followed by the addition of reagents downstream and detection of ions by quadrupole mass spectrometry. Most of the reactions that were studied are important for describing the chemistry of meteor-ablated calcium in the earth's upper atmosphere. The following rate coefficients were measured: k(CaO(+) + O --> Ca(+) + O(2)) = (4.2 +/- 2.8) x 10(-11) at 197 K and (6.3 +/- 3.0) x 10(-11) at 294 K; k(CaO(+) + CO --> Ca(+) + CO(2), 294 K) = (2.8 +/- 1.5) x 10(-10); k(Ca(+).CO(2) + O(2) --> CaO(2)(+) + CO(2), 294 K) = (1.2 +/- 0.5) x10(-10); k(Ca(+).CO(2) + H(2)O --> Ca(+).H(2)O + CO(2)) = (13.0 +/- 4.0) x 10(-10); and k(Ca(+).H(2)O + O(2) --> CaO(2)(+) + H(2)O, 294 K) = (4.0 +/- 2.5) x 10(-10) cm(3) molecule(-1) s(-1). The quoted uncertainties are a combination of the 1sigma standard errors in the kinetic data and the systematic errors in the models used to extract the rate coefficients. Rate coefficients were also obtained for the following recombination (also termed association) reactions in He bath gas: k(Ca(+).CO(2) + CO(2) --> Ca(+).(CO(2))(2), 294 K) = (2.6 +/- 1.0) x 10(-29); k(Ca(+).H(2)O + H(2)O --> Ca(+).(H(2)O)(2)) = (1.6 +/- 1.1) x 10(-27); and k(CaO(2)(+) + O(2) --> CaO(2)(+).O(2)) < 1 x 10(-31) cm(6) molecule(-2) s(-1). These recombination rate coefficients, as well as those for the ligand-switching reactions listed above, were then interpreted using a combination of high level quantum chemistry calculations and RRKM theory using an inverse Laplace transform solution of the master equation. The surprisingly slow reaction between CaO(+) and O was explained using quantum chemistry calculations on the lowest (2)A', (2)A' and (4)A' potential energy surfaces. These calculations indicate that reaction mostly occurs on the (2)A' surface, leading to production of Ca(+)((2)S) + O(2)((1)Delta(g)). The importance of this reaction for controlling the lifetime of Ca(+) in the upper mesosphere and lower thermosphere is then discussed.     

Theoretical study of the gas-phase chemiionization reactions La + O and La + O2

The La + O and La + O 2 chemiionization reactions have been investigated with quantum chemical methods. For La + O 2(X (3)Sigma g) and La + O 2(a (1)Delta g), the chemiionization reaction La + O 2 --> LaO 2 (+) + e (-) has been shown to be endothermic and does not contribute to the experimental chemielectron spectra. For the La + O 2(X (3)Sigma g) reaction conditions, chemielectrons are produced by La + O 2 --> LaO + O, followed by La + O --> LaO (+) + e (-). This is supported by the same chemielectron band, arising from La + O --> LaO (+) + e (-), being observed from both the La + O( (3)P) and La + O 2(X (3)Sigma g) reaction conditions. For La + O 2(a (1)Delta g), a chemielectron band with higher electron kinetic energy than that obtained from La + O 2(X (3)Sigma g) is observed. This is attributed to production of O( (1)D) from the reaction La + O 2(a (1)Delta g) --> LaO + O( (1)D), followed by chemiionization via the reaction La + O( (1)D) --> LaO (+) + e (-). Potential energy curves are computed for a number of states of LaO, LaO* and LaO (+) to establish mechanisms for the observed La + O --> LaO (+) + e (-) chemiionization reactions.     

The Zeeman tuning of the A 6Sigma+-X 6Sigma+ transition of chromium monohydride

The magnetic tuning of the low-rotational levels of the A (6)Sigma(+) (v = 1 and 0) states of chromium monohydride, (52)CrH, have been experimentally investigated using optical spectroscopy of the (0, 0) and (1, 0) bands of the A (6)Sigma(+)-X (6)Sigma(+) transition. The tuning of the numerous low-rotational lines in the A (6)Sigma(+)-X (6)Sigma(+) (0, 0) band can be accurately modeled using a single set of g-factors (g(S) and g(l)) which are close to the expected values. In contrast, the g-factors for the A (6)Sigma(+) (v = 1) state required to model the magnetic tuning of low-rotational lines in the A (6)Sigma(+)-X (6)Sigma(+) (1, 0) band are strongly dependent upon rotational and fine structure component and the determined effective values for g(S) deviate significantly from 2.002. Interpretation of the quantum level variation of g(S) is presented. The magnetic hyperfine structure of the (0, 0) and (1, 0) bands of the A (6)Sigma(+)-X (6)Sigma(+) transition is analyzed to produce proton Fermi contact, b(F) and dipolar, c, magnetic hyperfine parameters of 19(1) MHz and 34(5) MHz for the A (6)Sigma(+) (v = 0) state and 21(2) MHz and 30(7) MHz for the A (6)Sigma(+) (v = 1) state.     

A kinetic study of the reactions FeO+ + O, Fe+.N2 + O, Fe+.O2 + O and FeO+ + CO: implications for sporadic E layers in the upper atmosphere

These gas-phase reactions were studied by pulsed laser ablation of an iron target to produce Fe(+) in a fast flow tube, with detection of the ions by quadrupole mass spectrometry. Fe(+).N(2) and Fe(+).O(2) were produced by injecting N(2) and O(2), respectively, into the flow tube. FeO(+) was produced from Fe(+) by addition of N(2)O, or by ligand-switching from Fe(+).N(2) following the addition of atomic O. The following rate coefficients were measured: k(FeO(+) + O --> Fe(+) + O(2), 186-294 K) = (3.2 +/- 1.5) x 10(-11); k(Fe(+).N(2) + O --> FeO(+)+ N(2), 294 K) = (4.6 +/- 2.5) x 10(-10); k(Fe(+).O(2) + O --> FeO(+) + O(2), 294 K) = (6.3 +/- 2.7) x 10(-11); and k(FeO(+) + CO --> Fe(+) + CO(2), 294 K) = (1.59 +/- 0.34) x 10(-10) cm(3) molecule(-1) s(-1), where the quoted uncertainties are a combination of the 1sigma standard errors in the kinetic data and the systematic experimental errors. The surprisingly slow reaction between FeO(+) and O is examined using ab initio quantum calculations of the relevant potential energy surfaces. The importance of this reaction for controlling the lifetime of sporadic E layers is then demonstrated using a model of the upper mesosphere and lower thermosphere.     

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.     

A kinetic study of the reactions of Fe+ with N2O, N2, O2, CO2 and H2O, and the ligand-switching reactions Fe+.X + Y --> Fe+.Y + X (X = N2, O2, CO2; Y = O2, H2O)

A series of reactions involving Fe(+) ions were studied by the pulsed laser ablation of an iron target, with detection of ions by quadrupole mass spectrometry at the downstream end of a fast flow tube. The reactions of Fe(+) with N(2)O, N(2) and O(2) were studied in order to benchmark this new technique. Extending measurements of the rate coefficient for Fe(+) + N(2)O from 773 K to 185 K shows that the reaction exhibits marked non-Arrhenius behaviour, which appears to be explained by excitation of the N(2)O bending vibrational modes. The recombination of Fe(+) with CO(2) and H(2)O in He was then studied over a range of pressure and temperature. The data were fitted by RRKM theory combined with ab initio quantum calculations on Fe(+).CO(2) and Fe(+).H(2)O, yielding the following results (120-400 K and 0-10(3) Torr). For Fe(+) + CO(2): k(rec,0) = 1.0 x 10(-29) (T/300 K)(-2.31) cm(6) molecule(-2) s(-1); k(rec,infinity) = 8.1 x 10(-10) cm(3) molecule(-1) s(-1). For Fe(+) + H(2)O: k(rec,0) = 5.3 x 10(-29) (T/300 K)(-2.02) cm(6) molecule(-2) s(-1); k(rec,infinity) = 2.1 x 10(-9) (T/300 K)(-0.41) cm(3) molecule(-1) s(-1). The uncertainty in these rate coefficients is determined using a Monte Carlo procedure. A series of exothermic ligand-switching reactions were also studied at 294 K: k(Fe(+).N(2) + O(2)) = (3.17 +/- 0.41) x 10(-10), k(Fe(+).CO(2) + O(2)) = (2.16 +/- 0.35) x 10(-10), k(Fe(+).N(2) + H(2)O) = (1.25 +/- 0.14) x 10(-9) and k(Fe(+).O(2) + H(2)O) = (8.79 +/- 1.30) x 10(-10) cm(3) molecule(-1) s(-1), which are all between 36 and 52% of their theoretical upper limits calculated from long-range capture theory. Finally, the role of these reactions in the chemistry of meteor-ablated iron in the upper atmosphere is discussed. The removal rates of Fe(+) by N(2), O(2), CO(2) and H(2)O at 90 km altitude are approximately 0.1, 0.07, 3 x 10(-4) and 1 x 10(-6) s(-1), respectively. The initially formed Fe(+).N(2) and Fe(+).O(2) are converted into the H(2)O complex at approximately 0.05 s(-1). Fe(+).H(2)O should therefore be the most abundant single-ligand Fe(+) complex in the mesosphere below 90 km.     

Theoretical study on the group 2 atoms + N2O reactions

The electronic structure aspects of the M (1S,3P) + N2O(X 1sigma+) (M = Be, Mg, Ca) reactions are investigated using the CASSCF/MRMP2 (complete active space SCF and the multireference M?ller-Plesset perturbation theory of the second order) computational methodology. The lowest adiabatic 1 1A' and 1 3A' potential energy surfaces (PESs) favor the bending dissociation mechanism of N2O in all studied cases. The rate-limiting channels are determined by the classical barriers that decrease in the series Be (8.9) > Mg (7.0) > Ca (1.2) kcal/mol, whereas the spin-forbidden reaction channels are found to be less important. A comparison with elaborated kinetic results (Plane et al. J. Phys. Chem. 1990, 94, 5255; Gas-Phase Metal Reactions; Elsevier: Amsterdam, 1992; Vinckier et al. J. Phys. Chem. A 1999, 103, 5328) on the Ca (1S) + N2O(X 1sigma+) reaction is presented, and the differences in the kinetic behavior of the title reactions are discussed. Our results also indicate that the techniques based on the multiconfigurational wave functions are unavoidable if a correct topology of the PESs governing these reactions is needed.     

Kinetic study of the reaction Ca+ + N2O from 188 to 1207 K

Ion-molecule reactions involving metallic species play a central role in the chemistry of planetary ionospheres and in many combustion processes. The kinetics of the Ca(+) + N(2)O --> CaO(+) + N(2) reaction was studied by the pulsed multiphoton dissociation at 193 nm of organo-calcium vapor in the presence of N(2)O, followed by time-resolved laser-induced fluorescence spectroscopy of Ca(+) at 393.37 nm (4(2)P(3/2) <-- 4(2)S(1/2)). This yielded k(188-1207 K) = 5.45 x 10(-11) (T/300 K)(0.53) exp(282 K/T) cm(3) molecule(-1) s(-1), with an estimated accuracy of +/-13% (188-600 K) and +/-27% (600-1207 K). The temperature dependence of this barrierless reaction, with a minimum in the rate coefficient between 400 and 600 K, appears to be explained by the role of N(2)O vibrational excitation. This is examined using a classical trajectory treatment on a potential energy surface calculated at the B3LYP/6-311+g(2d,p) level of theory.     

Theoretical analysis of structural, energetic, electronic, and defect properties of Li2O

The structural, energetic, and electronic properties of stoichiometric and defective Li(2)O were studied theoretically. The reliability of the Perdew-Wang method in the framework of density functional theory (DFT), and of two DFT/Hartree-Fock hybrid methods (PW1PW and B3LYP), was examined by comparison of calculated and available experimental data. Atom-centered orbitals and plane waves were used as basis functions for the crystalline orbitals. For both cases, the basis set dependence of calculated properties was investigated. With most of the methods, good agreement with the experimental Li(2)O lattice parameter and cohesive energy was obtained. In accordance with experiment, the analysis of electronic properties shows that Li(2)O is a wide gap insulator. Among the considered methods, the hybrid methods PW1PW and B3LYP give the best agreement with experiment for the band gap. The formation of an isolated cation vacancy defect and an F center in Li(2)O were studied. The effect of local relaxation on the calculated defect formation energies and the defect-induced changes of electronic properties were investigated and compared to available experimental results. The migration of a Li(+) ion in Li(2)O bulk was investigated. The activation energy for the migration of a Li(+) ion from its regular tetrahedral site to an adjacent cation vacancy was calculated, including the effect of local relaxation. The calculated activation barriers, 0.27-0.33 eV, are in excellent agreement with experiment.     

H-bond network in amino acid cocrystals with H2O or H2O2. The DFT study of serine-H2O and serine-H2O2

The structure, IR spectrum, and H-bond network in the serine-H(2)O and serine-H(2)O(2) crystals were studied using DFT computations with periodic boundary conditions. Two different basis sets were used: the all-electron Gaussian-type orbital basis set and the plane wave basis set. Computed frequencies of the IR-active vibrations of the titled crystals are quite different in the range of 10-100 cm(-1). Harmonic approximation fails to reproduce IR active bands in the 2500-2800 frequency region of serine-H(2)O and serine-H(2)O(2). The bands around 2500 and 2700 cm(-1) do exist in the anharmonic IR spectra and are caused by the first overtone of the OH bending vibrations of H(2)O and a combination vibration of the symmetric and asymmetric bendings of H(2)O(2). The quantum-topological analysis of the crystalline electron density enables us to describe quantitatively the H-bond network. It is much more complex in the title crystals than in a serine crystal. Appearance of water leads to an increase of the energy of the amino acid-amino acid interactions, up to ~50 kJ/mol. The energy of the amino acid-water H-bonds is ~30 kJ/mol. The H(2)O/H(2)O(2) substitution does not change the H-bond network; however, the energy of the amino acid-H(2)O(2) contacts increases up to 60 kJ/mol. This is caused by the fact that H(2)O(2) is a much better proton donor than H(2)O in the title crystals.     

A kinetic study of Mg+ and Mg-containing ions reacting with O3, O2, N2, CO2, N2O and H2O: implications for magnesium ion chemistry in the upper atmosphere

Reactions between Mg(+) and O(3), O(2), N(2), CO(2) and N(2)O were studied using the pulsed laser photo-dissociation at 193 nm of Mg(C(5)H(7)O(2))(2) vapour, followed by time-resolved laser-induced fluorescence of Mg(+) at 279.6 nm (Mg(+)(3(2)P(3/2)-3(2)S(1/2))). The rate coefficient for the reaction Mg(+) + O(3) is at the Langevin capture rate coefficient and independent of temperature, k(190-340 K) = (1.17 ± 0.19) × 10(-9) cm(3) molecule(-1) s(-1) (1σ error). The reaction MgO(+) + O(3) is also fast, k(295 K) = (8.5 ± 1.5) × 10(-10) cm(3) molecule(-1) s(-1), and produces Mg(+) + 2O(2) with a branching ratio of (0.35 ± 0.21), the major channel forming MgO(2)(+) + O(2). Rate data for Mg(+) recombination reactions yielded the following low-pressure limiting rate coefficients: k(Mg(+) + N(2)) = 2.7 × 10(-31) (T/300 K)(-1.88); k(Mg(+) + O(2)) = 4.1 × 10(-31) (T/300 K)(-1.65); k(Mg(+) + CO(2)) = 7.3 × 10(-30) (T/300 K)(-1.59); k(Mg(+) + N(2)O) = 1.9 × 10(-30) (T/300 K)(-2.51) cm(6) molecule(-2) s(-1), with 1σ errors of ±15%. Reactions involving molecular Mg-containing ions were then studied at 295 K by the pulsed laser ablation of a magnesite target in a fast flow tube, with mass spectrometric detection. Rate coefficients for the following ligand-switching reactions were measured: k(Mg(+)·CO(2) + H(2)O → Mg(+)·H(2)O + CO(2)) = (5.1 ± 0.9) × 10(-11); k(MgO(2)(+) + H(2)O → Mg(+)·H(2)O + O(2)) = (1.9 ± 0.6) × 10(-11); k(Mg(+)·N(2) + O(2)→ Mg(+)·O(2) + N(2)) = (3.5 ± 1.5) × 10(-12) cm(3) molecule(-1) s(-1). Low-pressure limiting rate coefficients were obtained for the following recombination reactions in He: k(MgO(2)(+) + O(2)) = 9.0 × 10(-30) (T/300 K)(-3.80); k(Mg(+)·CO(2) + CO(2)) = 2.3 × 10(-29) (T/300 K)(-5.08); k(Mg(+)·H(2)O + H(2)O) = 3.0 × 10(-28) (T/300 K)(-3.96); k(MgO(2)(+) + N(2)) = 4.7 × 10(-30) (T/300 K)(-3.75); k(MgO(2)(+) + CO(2)) = 6.6 × 10(-29) (T/300 K)(-4.18); k(Mg(+)·H(2)O + O(2)) = 1.2 × 10(-27) (T/300 K)(-4.13) cm(6) molecule(-2) s(-1). The implications of these results for magnesium ion chemistry in the atmosphere are discussed.     

Simulation of Ca2+ and Mg2+ solvation using polarizable atomic multipole potential

The alkaline earth metals calcium and magnesium are critically involved in many biomolecular processes. To understand the hydration thermodynamics of these ions, we have performed molecular dynamics simulations using a polarizable potential. Particle-mesh Ewald for point multipoles has been applied to the calculation of electrostatic interactions. The parameters in this model have been determined from an ab initio quantum mechanical calculation of dimer interactions between ions and water. Two methods for ion solvation free energy calculation, free energy perturbation, and the Bennett acceptance ratio have been compared. Both predict results consistent with other theoretical estimations while the Bennett approach leads to a much smaller statistical error. Based on the Born theory and the ion-oxygen radial distribution functions, we estimate the effective size of the ions in solution, concluding that K(+) > Na(+) congruent with Ca(2+) > Mg(2+). There appears to be much stronger perturbation in water structure, dynamics, and dipole moment around the divalent cations than the monovalent K(+) and Na(+). The average water coordination numbers for Ca(2+) and Mg(2+) are 7.3 and 6, respectively. The lifetime of water molecules in the first solvation shell of Mg(2+) is on the order of hundreds of picoseconds, in contrast to only few picoseconds for Ca(2+), K(+), or Na(+).     

Resonant States of H3+ and D2H+

Vibrational resonances for H(3) (+) and D(2)H(+), as well as H(3) (+) at J=3, are calculated using a complex absorbing potential (CAP) method with an automated procedure to find stability points in the complex plane. Two different CAP functional forms and different CAP extents are used to analyze the consistency of the results. Calculations are performed using discrete variable representation continuum basis elements calculated to high levels of accuracy by diagonalizing large, dense, Hamiltonian matrices. For D(2)H(+), two energy regions are analyzed: the one where D(2)+H(+) is the only dissociation product and the one where HD+D(+) can also be formed. Branching ratios are obtained in the latter case by using different CAPs. It is shown that H(3) (+) and D(2)H(+) support some narrow Feshbach-type resonances but that higher angular momentum states must be studied to model the pre-dissociation spectrum recorded by Carrington and co-workers [J. Chem. Phys. 98, 1073 (1993)].     

Ab initio potential energy and dipole moment surfaces for H5O2 +

Full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for H(5)O(2) (+). Tens of thousands of coupled-cluster [CCSD(T)] and second-order Moller-Plesset (MP2) calculations of electronic energies, using aug-cc-pVTZ basis, were done. The energies were fit very precisely in terms of all the internuclear distances, using standard least-square procedures, however, with a fitting basis that satisfies permutational symmetry with respect to like atoms. The H(5)O(2) (+) PES is a fit to 48 189 CCSD(T) energies, containing 7962 polynomial coefficients. The PES has a rms fitting error of 34.9 cm(-1) for the entire data set up to 110 000 cm(-1). This surface can describe various internal floppy motions, including the H atom exchanges, monomer inversions, and monomer torsions. First- and higher-order saddle points have been located on the surface and compared with available previous theoretical work. In addition, the PES dissociates correctly (and symmetrically) to H(2)O+H(3)O(+), with D(e)=11 923.8 cm(-1). Geometrical and vibrational properties of the monomer fragments are presented. The corresponding global DMS fit (MP2 based) involves 3844 polynomial coefficients and also dissociates correctly.     

Kinetic studies of atmospherically relevant silicon chemistry. Part III: Reactions of Si+ and SiO+ with O3, and Si+ with O2

Silicon ions are generated in the Earth's upper atmosphere by hyperthermal collisions of material ablated from incoming meteoroids with atmospheric molecules, and from charge transfer of silicon-bearing neutral species with major atmospheric ions. Reported Si(+) number density vs. height profiles show a sharp decrease below 95 km, which has been commonly attributed to the fast reaction with H(2)O. Here we report rate coefficients and branching ratios of the reactions of Si(+) and SiO(+) with O(3), measured using a flow tube with a laser ablation source and detection of ions by quadrupole mass spectrometry. The results obtained are (2σ uncertainty): k(Si(+) + O(3), 298 K) = (6.5 ± 2.1) × 10(-10) cm(3) molecule(-1) s(-1), with three product channels (branching ratios): SiO(+) + O(2) (0.52 ± 0.24), SiO + O(2)(+) (0.48 ± 0.24), and SiO(2)(+) + O (<0.1); k(SiO(+) + O(3), 298 K) = (6 ± 4) × 10(-10) cm(3) molecule(-1) s(-1), where the major products (branching ratio ≥ 0.95) are SiO(2) + O(2)(+). Reactions (1) and (2) therefore have the unusual ability to neutralise silicon directly, as well as forming molecular ions which can undergo dissociative recombination with electrons. These reactions, along with the recently reported reaction between Si(+) and O(2)((1)Δ(g)), largely explain the disappearance of Si(+) below 95 km in the atmosphere, relative to other major meteoric ions such as Fe(+) and Mg(+). The rate coefficient of the Si(+) + O(2) + He reaction was measured to be k(298 K) = (9.0±1.3) × 10(-30) cm(6) molecule(-2) s(-1), in agreement with previous measurements. The SiO(2)(+) species produced from this reaction, which could be vibrationally excited, is observed to charge transfer at a relatively slow rate with O(2), with a rate constant of k(298 K) = (1.5 ± 1.0) × 10(-13) cm(3) molecule(-1) s(-1).     

MRCI investigation of different isomers of Ni2O2H2(+)

The geometrical and electronic structures of different isomers of Ni(2)O(2)H(2)(+) are investigated by multireference configuration interaction (MRCI) calculations using natural atomic orbital basis sets. The lowest-lying isomer, Ni(2)(OH)(2)(+), has a rhombic shape with two OH groups bridging the Ni atoms. The next isomer in energetic order with a relative energy of 0.29 eV consists of a linear NiONi(OH(2))(+) chain. Other structures with a rhombic shape, (NiH)(2)O(2)(+), with H bound to the Ni atoms have considerably higher energies, above 4 eV. Especially the low-lying isomers are characterised by a large number of low-lying electronic terms. The product Ni(2)O(2)H(2)(+) of the reaction of Ni(2)O(2)(+) with small alkanes is likely to have the rhombic Ni(2)(OH)(2)(+) structure. The reaction energy of the reaction Ni(2)O(2)(+) + H(2)→ Ni(2)(OH)(2)(+) is estimated to be about -3.5 eV.