Bimolecular reactions of S 2+ with Ar, H 2 and N 2 : reactivity and dynamics †

The reactivity, energetics and dynamics of bimolecular reactions between S 2+ and three neutral species (Ar, H 2 and N 2 ) have been studied using a position-sensitive coincidence methodology at centre-of-mass collision energies below 6 eV. This is the first study of bimolecular reactions involving S 2+ , a species detected in planetary ionospheres, the interstellar medium, and in anthropogenic manufacturing processes. The reactant dication beam employed consists predominantly of S 2+ in the ground 3 P state, but some excited states are also present. Most of the observed reactions involve the ground state of S 2+ , but the dissociative electron transfer reactions appear to exclusively involve excited states of this atomic dication. We observe exclusively single electron-transfer between S 2+ and Ar, a process which exhibits strong forward scatting typical of the Landau–Zener style dynamics observed for other dicationic electron transfer reactions. Following collisions between S 2+ + H 2 , non-dissociative and dissociative single electron-transfer reactions were detected. The dynamics here show evidence for the formation of a long-lived collision complex, [SH 2 ] 2+ , in the dissociative single electron-transfer channel. The formation of SH + was not observed. In contrast, the collisions of S 2+ + N 2 result in the formation of SN + + N + in addition to the products of single electron-transfer reactions.


Introduction
Doubly charged positive ions (dications) can be found in energised environments such as the ionospheres of planets and their satellites, 1-8 the interstellar medium, [9][10][11] and in anthropogenic plasmas. [12][13][14] Significant bimolecular reactivity is shown following collisions of both atomic and molecular dications with neutral species in the gas-phase. [15][16][17][18][19][20] The lifetimes of atomic dications in planetary ionospheres are expected to be primarily determined by such collisional processes. 21 In addition, despite their usual inherent thermodynamic instability, the metastable electronic states of molecular dications have been shown to possess lifetimes sufficient to allow collisions with other species in ionospheres and elsewhere. 15 Given the above lifetimes, the demonstrated bimolecular reactivity of these dications strongly suggests that dication chemistry can play a role in ionospheric processes; 22 for example, dications could be involved in the chemistry of complex molecule formation through reactions involving carbon chaingrowth. 21,[23][24][25][26] There are numerous examples of the detection of atomic dications in planetary ionospheres. 27 However, it is difficult to unambiguously detect molecular dications in ionospheres using simple mass spectrometry, the usual sampling technique. This difficulty arises because these environments often contain monocations with the same mass to charge ratio as the target dication. 15 The historical absence of a definitive detection of ionospheric molecular dications may account for the relative neglect of these species in models of ionosphere chemistry. 21 Encouragingly, the presence of the CO 2 2+ dication in the ionosphere of Mars has been recently confirmed, the first detection of a molecular dication in such an environment. 28 In order to identify dication reactions of ionospheric interest, laboratory-based experiments to probe dicationic reactivity, along with further in situ detection of dications in ionospheres, are vital. 29 Indeed, recognition of the presence of molecular dications in planetary atmospheres has led to the identification of the role these species can play in atmospheric erosion. [30][31][32][33] To develop further our understanding of dicationic processes, this paper presents an investigation of the reactions of S 2+ with Ar, H 2 and N 2 . This work reveals the reactivity, reaction mechanisms and reaction energetics of S 2+ , extending and enhancing our understanding of the chemistry of dications in environments such as planetary atmospheres. Sulfur-containing species are found in a wide variety of astrophysical environments [34][35][36][37] and on Earth, where they are ubiquitous both in nature and from anthropogenic activities. 38,39 The S 2+ ion has been observed in both terrestrial and non-terrestrial environments. For example, when plasmas involving SF 6 are used for etching, S 2+ can be formed. 40,41 In interstellar clouds, organic compounds containing sulfur have been detected, including species with S-X bonds where X = O, C, N, H, 34 and S 2+ has been observed in nebulae close to sources of high energy photons. [42][43][44] Indeed, the reactions of the sulfur dication are thought to be important in the ionosphere of Io, the satellite of Jupiter, where S 2+ has been detected by the Galileo 27 and Voyager I 45,46 spacecraft, as well as by the Hubble telescope. 47 Io's sulfur-rich atmosphere is derived from this moon's intense volcanic activity which results in the expulsion of SO 2 . 48,49 Hence, S 2+ is thought to be one of the major species in the Io torus, 50 a region of the Jovian system where charged species generated in Io's ionosphere, such as S 2+ , are transferred to Jupiter, dominating both the planets magnetosphere 49,51 and the torus of Europa. 52 Despite this clear relevance of S 2+ to astrophysical and anthropogenic environments, and the longevity of the atomic dication, there are no reports in the literature (to the authors' knowledge) of the bimolecular reactivity of S 2+ dications in the gas phase. Therefore, the investigation of the collisions of S 2+ with neutral molecules presented in this study provides the first experimental information on the reactivity of this ion; that information is particularly pertinent to the chemistry of a number of astrophysical and terrestrial environments.
As a first investigation of the chemistry of S 2+ , we have chosen to study the reactions with three neutral gases that are important in interstellar and atmospheric chemistry: Ar, H 2 and N 2 . Argon is often present in significant concentrations in planetary atmospheres. 53-55 H 2 is the most abundant molecule in the universe and is the dominant species making up the gaseous planets, including Jupiter. [56][57][58] Molecular nitrogen (N 2 ) is important in the atmospheres of terrestrial bodies in the solar system, especially the Earth and Titan where it is the dominant species. 21,22,54,55,[59][60][61][62][63] The reactions resulting from dication collisions with argon have been well studied. Single electron-transfer (SET) was observed in collisions with rare gas dications: Ne 2+ , Ar 2+ , Kr 2+ or Xe 2+ . [64][65][66][67][68][69] Collisions between N 2 2+ and Ar also resulted in the observation of SET reactions. 70 Collisions between Ar and hydrogen-containing dications, such as CHX 2+ (X = F, Cl, Br, I), often result in SET and proton-transfer (PT) reactions; the latter, for example, forming ArH + . [71][72][73] Interactions between Ar and SF x 2+ (x = 2-4) or CF y 2+ (y = 1-3) have been the subject of several studies, with SET and collision-induced dissociation (CID) channels being observed. 74 and CF 3 2+ respectively. 79,80 A variety of dication reactions with H 2 have been studied at low collision energies, often resulting in SET reactions. Reactions between Ar 2+ and H 2 resulted in non-dissociative and dissociative SET (NDSET and DSET respectively). 64,66,81 With some H-containing dications (for example CHCl 2+ ) H + transfer to the neutral H 2 is observed, producing H 3 + via the formation of collision complexes. [82][83][84] Direct hydride (H À ) transfer from the H 2 to the dication is also seen in some systems. 83,85 Many other bond-forming channels have been observed in dication reactions with H 2 , often with evidence that the reactions proceed via complexation. 17,[86][87][88][89] Some of the new bonds  formed in this class of reaction include N-H, 89 H-H, 84  C-H, 82,83,85-88 F-H, 17,87,88 Cl-H, 83 and Br-H. 82 The ubiquity of bond-forming reactions in dication/H 2 collision systems, and the prevalence of S-H bonds in nature, marks out the S 2+ + H 2 collision system as a candidate for the generation of new bonds. The reactions following dication collisions with N 2 have been the subject of several previous studies, including recent work by the current authors. 90 In the Ar 2+ + N 2 collision system, SET, double electron-transfer (DET) and bond-forming pathways were observed. 64,[66][67][68]81,[90][91][92] SET channels were also reported following the interactions of N 2 with CO 2+ and N 2

2+
. 93,94 DSET resulting from the collisions of Ne 2+ + N 2 was shown to occur via two mechanisms, one involving a collision complex. 95 Bond-forming reactions have been observed following the interactions of N 2 with C 4 H 3 2+ and O 2 2+ , resulting in the formation of N-H and N-O bonds. 84,96 Additionally, S-N bondformation has been observed following the gas-phase collisions of S + ions with ammonia. 97 In our experiments, we collide S 2+ ions with the selected neutral species using a Position-Sensitive Coincidence Mass Spectrometer (PSCO-MS), at centre-of-mass (CM) collision energies of less than 6 eV. The PSCO-MS combines coincident detection of product cations with a crossed-beam methodology. The results from our experiments give a detailed insight into the reactivity and dynamics of dication-neutral interactions. For example, SET reactions are observed following the interactions of S 2+ with Ar, N 2 and H 2 , including both non-dissociative and dissociative channels with the molecular targets. The collisions of S 2+ with N 2 also result in the formation of SN + + N + via a collision complex. The dynamics of the DSET channel in the S 2+ + H 2 system also shows evidence for the formation of a collision complex, [SH 2 ] 2+ .

Experimental
Coincidence techniques involve the simultaneous detection of two or more products from a single reactive event. Dication interactions with neutrals often generate pairs of monocations, and these pairs of ions are detected in coincidence in our PSCO-MS experiment, which has been described in detail in the literature. 76,98,99 In brief, a pulsed beam of dications is directed into the field-free source region of a time-of-flight mass spectrometer (TOF-MS). In this region the dications interact with a jet of the neutral reactant. Subsequent application of an extraction voltage to the source region allows the TOF-MS to detect the cation pairs generated from the dication-neutral interactions. The detection of these ions involves recording their arrival time, and position, at a large microchannel-plate detector. From this raw data, a list of flight times and arrival positions of the ions detected in pairs, a two-dimensional mass spectrum, can be generated revealing the different reactive channels. The positional data accompanying the ionic detections yields the relative motion of the products of each reactive event, providing a detailed insight into the mechanisms of each reactive channel. 99 The S 2+ ions used in the experiments described in this paper were generated via electron ionization of H 2 S (CKgas, 99.5%) by 100 eV electrons in a custom-built ion source. Previous investigations show that S 2+ (m/z = 16), HS 2+ (m/z = 16.5), and H 2 S 2+ (m/z = 17) dications are all formed following the bombardment of H 2 S with electrons. 100 The positively charged ions are extracted from the ion source and pass through a hemispherical energy analyser to restrict the translational energy spread of the final S 2+ beam to B0.3 eV. The continuous beam of ions exiting the hemispherical analyser is then pulsed, using a set of electrostatic deflectors, before being accelerated and focussed into a commercial velocity filter. The velocity filter is set to transmit just the 32 S 2+ (m/z = 16) ions. O + ions, generated from the ionization of background gases in the ion source, may also be present in our dication beam because O + has the same m/z ratio as 32 S 2+ . However, reactions involving O + ions will not result in the generation of two cationic products. Therefore, the products of reactions involving these contaminant O + ions will not contribute to our coincidence spectra. The resulting pulsed beam of energy-constrained S 2+ ions is then decelerated to o10 eV in the laboratory frame before entering the source region of the TOF-MS. In the TOF-MS source region the beam of dications is crossed with an effusive jet of the appropriate neutral species: Ar (BOC, 99.998%), H 2 (BOC, 499.995%), or N 2 (BOC, 499.998%). Single-collision conditions 101 are achieved by employing an appropriately low pressure of the neutral collision partner. Under these conditions most dications do not undergo a collision and only a small percentage experience one collision. Such a pressure regime ensures no secondary reactions, due to successive collisions with two neutral species, influence the S 2+ reactivity we observe. An electric field is applied across the TOF-MS source when the dication pulse reaches the centre of this region. This electric field accelerates positively charged species into the second electric field (acceleration region) of the TOF-MS and then on into the flight tube. At the end of the flight tube, the cations are detected by the position-sensitive detector which comprises a chevron-pair of microchannel plates (diameter = 12.7 cm) located in front of a dual delay-line anode. 98 The voltage pulse applied to the source region also starts the ion timing circuitry, to which the signals from the detector provide stop pulses. The experiments in this work employed TOF-MS source fields of 183 V cm À1 . The electrostatic fields employed in the TOF region are homogeneous in order to achieve high mass resolution, therefore our PSCO technique is distinct from 'velocity map imaging'.
Signals from the detector are amplified and discriminated before being passed to a PC-based time-to-digital converter. If two ions are observed in the same TOF cycle, a coincidence event is recorded and each ion's arrival time and impact position on the detector are stored for off-line analysis. The use of single-collision conditions ensures 'false' coincidences are kept to a minimum in our spectra. The ion pairs data can be plotted as a 2D histogram, a 'pairs spectrum', where the time of flights (t 1 , t 2 ) of each ion in the pair are used as the (x, y) coordinates. Peaks in the pairs spectrum readily identify bimolecular reaction channels that form a pair of positively-charged product ions. The subset of events associated with each such peak, events corresponding to an individual reaction channel, can then be selected for further off-line analysis.
As shown in previous work, the positional and time of flight information for each ion of a pair can be used to generate their x, y and z velocity vectors in the laboratory frame; here the z-axis is defined by the principal axis of the TOF-MS. 98 The x and y velocity vectors of an ion are determined from the positional information and flight time; the z vector is determined from the deviation of the observed TOF from the expected TOF of an identical ion with zero initial kinetic energy. The laboratory frame velocities are then converted into the CM frame using the initial dication velocity. 98 Often the pair of monocations resulting from the reaction between a dication and a neutral are accompanied by a neutral species: a three-body reaction. A powerful feature of the PSCO-MS experiment is that the CM velocity of such a neutral product can be determined from the CM velocities of the detected ionic products via conservation of momentum. 98 To reveal the dynamics of a given reaction channel, a CM scattering diagram (Fig. 1) can be generated from the velocities of the product ions. Such CM scattering diagrams are radial histograms that, for each event collected for a given reaction channel, plot the magnitude of the products' CM velocity |w i | as the radial co-ordinate and the scattering angle y between w i and the CM velocity of the incident dication as the angular coordinate. In the kinematics that apply in our experiment, where the dication has significantly more momentum than the neutrals, the velocity of the incident dication is closely oriented with the velocity of the centre of mass. In our CM scattering diagrams, since 01 r y r 1801, the data for one product can be shown in the upper semi-circle of the figure and the data for another product in the lower semi-circle, as the scattering of each ion is azimuthally symmetric. The typical angular resolution of the scattering data achieved by the PSCO apparatus operating with a high source field is 41. 69 It should be noted that in our data treatment, to enable presentation in a two-dimensional figure on the page, we are integrating over the azimuthal angle (in which the scattering is cylindrically symmetric). That is, we are binning events according to their value of y irrespective of the azimuthal scattering angle. This integration over the azimuthal angle results in an isotropic scattering distribution giving a sin(y) intensity distribution in our scattering diagrams, giving our experimental arrangement a lower relative detection efficiency at values of y very close to 01 or 1801. This phenomenon is ubiquitous where integration over the azimuthal angle occurs and is clearly visible in the scattering diagrams, where even strongly forward scattered reactions exhibit a peak in y away from y = 0.
In all the scattering diagrams presented in this work, a logarithmic scale is used for the scattering intensity. Such scales allow the diagrams to reveal subtleties of the low intensity scattering. However, we should be aware that such logarithmic scales can, at first glance, overemphasise the importance of the areas in the scattering diagrams with low intensities. In the ESI † ( Fig. SI 1) we contrast scattering diagrams with linear and logarithmic intensity scales.
For three-body reactions, internal-frame scattering diagrams can be a powerful aid in interpreting the reaction dynamics. In this class of scattering diagram |w i | is again the radial coordinate, but the angular coordinate is now the CM scattering angle with respect to CM velocity of one of the other product species.
From the CM velocities of the product species the total kinetic energy release (KER) T for a given reactive event can also be determined using the individual CM velocities of the products. 98 The exoergicity of the reaction DE can then be determined from T and the CM collision energy, E com : where E products and E reactants are the relative energies of the product and reactant states respectively. If the products lie lower in energy than the reactants, the resulting exoergicity will be positive. Performing this analysis for all the events collected for a given reaction channel provides a histogram of the experimental exoergicities of the detected reactive events (e.g. Fig. 2). Given knowledge of the available electronic states of the reactants and products, the experimental exoergicity spectrum can reveal the electronic states involved in the reaction. Previous work has shown that such experimental exoergicity spectra can clearly resolve states separated by 1 eV. 102 To interpret the experimental exoergicity spectra presented below we need to consider the electronic states of S 2+ populated in our dication beam. Beams of S 2+ ions generated by EI are not well characterised, but electronic state information is, of course, available for the S 2+ dication from atomic spectroscopy. 103,104 Experiments studying the electron ionization of H 2 S have previously reported the formation of S 2+ , but provide little information on the S 2+ states populated. 100,105,106 In our experiments, after the formation of S 2+ from H 2 S via electron ionization in the source region, the S 2+ dications must travel through the apparatus to the interaction region where they can collide with the appropriate neutral species. In the time taken to travel to the interaction region (B100 ms), one would expect S 2+ formed in higher lying excited states to radiatively decay to lower energy, metastable, states, where such a decay is 'allowed'. Hence, it seems reasonable that the S 2+ ions in our beam predominantly exist in one of the [Ne]3s 2 3p 2 states ( 3 P, 1 D or 1 S) from the ground configuration. Such a deduction is in strong accord with experimental determination of the dication states populated when generating the analogous (p 4 ) states of rare gas dications. 107 Here, electronic states from the ground configuration dominate the dication beam, with the associated electronic states having populations approximately in accord with their degeneracy 9[ 3 P]:5 Ground state S 2+ ( 3 P) has an energy of 33.7 eV relative to a neutral S atom in its ground ( 1 S) state. The excited 1 D and 1 S states lie 35.1 and 37.1 eV above the ground state of S respectively. S 2+ also has some higher energy states which are potentially long-lived and could therefore also be present in our dication beam: the 5 S state derived from the [Ne]3s 1 3p 3 configuration, which lies 41.0 eV above the ground state of S, and the 3 F state, derived from the [Ne]3s 2 3p 1 3d 1 configuration, lies 48.8 eV above the ground state of S. Emission lines attributed to decay from the 5 S state of S 2+ have been observed experimentally, including in the atmosphere Io. [108][109][110] These investigations reveal the 5 S state is metastable, with estimated lifetimes on the order of B100 ms, and therefore is possibly a minority species in our S 2+ beam. 111 The 3 F state has a shorter lifetime than 5 S, but could also be present in the dication beam. 103,111 Given the above, and also noting that the cross section for formation of the various S 2+ electronic states will increase with excess energy in this ionizing energy regime favouring the population of the p 2 states, we expect our dication beam to be dominated by ions in the ground 3 P state (60%) with the next most abundant state 1 D (30%), with minor contribution (10%) from the 1 S state. Small quantities of the excited 5 S and 3 F states may be present, but likely with a markedly lower abundance than the 1 S state. We will see below that our experimental exoergicity spectra strongly support this deduced beam composition, with literature exoergicities associated with reactions of the 3 P state accounting for the only observed reaction with Ar, the dominant reaction with H 2 and 98% of the product flux with N 2 .

Results and discussion
The collisions of S 2+ with Ar PSCO-MS spectra were recorded following the collisions of S 2+ with Ar at E cm = 5.3 eV. As with all reactions of S 2+ , there is no prior literature with which to compare our results. The only peak in the coincidence pairs spectrum, S + + Ar + , results from single electron-transfer (SET). Fig. 1 shows the scattering of the S + and Ar + product ions from this SET reaction, revealing strong 'forward' scattering, where the velocity of the S + product is oriented in the same direction as the velocity of the reactant dication, w(S 2+ ). This forward scattering has been commonly observed for other dicationic SET processes and arises from a direct electron-transfer, where the electron is transferred at a significant interspecies separation (3-6 Å, the so-called reaction window). Such processes are typically well described, for a wide range of collision partners by a Landau-Zener (LZ) formalism. 69,99,112,113 This LZ approach, and the general applicability of the concept of the reaction window, has been established by studies involving a wide range of different ions and targets. 112,114 Fig. 2 shows a histogram of the event exoergicities recorded in this SET channel. The bulk of the structure in this spectrum is between 4.8 eV and 7.2 eV. In order to interpret the experimental exoergicity distribution, we must consider the possible reactant and product electronic states involved. We discussed above the fact our beam is primarily composed of dications in the p 2 configuration, and likely dominated by the ground S 2+ ( 3 P) state. It should be noted that in dication-neutral reactions, in this energy regime, the collision energy (E cm ) generally does not couple efficiently to the potential energy surface. 15,112 Therefore, E cm does not usually need to be included when considering the energetics of these reactions. In the ESI, † we present in detail the accessible reactant and product states for this reaction and their associated literature exoergicities, a significant number of which lie in the reaction window (2-7 eV). However, given the arguments above, we see that pathways (a)-(d) of the 3 P and 1 D states of S 2+ populating the ground state of Ar + and the ground and first excited states of S + , correlate well with the experimental exoergicity distribution (Fig. 2). These dication states are expected to dominate our beam and the product states involved are readily accessible. Thus, pathways (a)-(d) contributing the bulk of the signals seems highly probable.
Pathways In summary, the above analysis shows that reactions between S 2+ and Ar are dominated by single electron-transfer with dynamics and energetics that conform well with the established models of dicationic electron transfer. The exoergicities observed support the arguments that our S 2+ beam is primarily composed of ions from the ground (p 2 ) configuration. Of course, small contributions from the higher excited states of S 2+ may be present but would be swamped by the signals from the 3 P and 1 D states.
The collisions of S 2+ with H 2 PSCO-MS spectra were recorded following the collisions of S 2+ and H 2 at E cm = 0.71 eV. Two reaction channels were observed in the coincidence pairs spectrum ( Table 1). The dominant channel, Rxn I, is non-dissociative single electron-transfer (NDSET), producing S + + H 2 + . A dissociative single electrontransfer (DSET) reaction is also observed (Rxn II), resulting in the formation of S + + H + + H. We note, in passing, that the pairs count rate for the S 2+ /H 2 collision system was markedly lower than the experiments involving Ar or N 2 .
There is no peak in the coincidence spectrum resulting from the transfer of two electrons from H 2 to S 2+ , resulting in the formation of H + + H + . Previous work has shown that in dicationneutral systems, facile double electron-transfer (DET) involves a concerted mechanism in which the reactant and product asymptotes have to lie close in energy (o1 eV). 115 Whilst the recombination energy from S 2+ gaining two electrons is enough (431.7 eV) to form H + + H + sequentially from H 2 , the (vertical) double ionization potential of H 2 is B51 eV, significantly higher in energy than is accessible from even the 3 F state of S 2+ . 103,116 Therefore, we would not expect the S 2+ + H 2 collisions in our experiment to result in DET via a concerted mechanism, in agreement with the experimental data. Hence, the absence of a DET channel in this collision system, despite the accessibility of a step-wise mechanism, adds to the increasing evidence that, in dication-neutral systems, DET reactions occur via a concerted two-electron transfer. 115 We note that there is no peak in the coincidence spectrum corresponding to SH + + H + , the products of a bond-forming reaction. SH + is observed in star-forming regions 117 and diffuse interstellar clouds where it is thought to be involved in ionneutral chemistry resulting in the formation of a range of sulfur containing species. 118 The formation of SH + + H + from S 2+ + H 2 is exothermic by B9 eV, therefore we might expect to observe SH + + H + in our experiment, 103,119,120 especially considering the ubiquity of the S-H bond in biological systems, and the propensity of dication + H 2 reactions to result in bondformation. 17,[82][83][84][85][86][87][88][89] The fact we do not observe the formation of SH + + H + could be because in order to transition to the SH + + H + asymptote, the S 2+ + H 2 system has to traverse a region of the potential energy surface with ready access to SET channels. Indeed, the model developed by Herman et al., 86 to explain the competition between SET and bond-forming reactivity in dication-neutral collision systems, shows that the presence of favourable curve crossings that result in SET reduce the probability of bond-forming reactions occurring. The fact we do not observe SH + following the collisions of S 2+ with H 2 suggests that this route is unlikely to contribute to the formation of SH + observed in astrophysical environments.

Non-dissociative single electron-transfer (NDSET)
The most intense product channel we observe following collisions between S 2+ and H 2 is NDSET, forming S + and H 2 + . Fig. 3 shows the CM scattering of the S + and H 2 + products in this NDSET channel. It should be noted that because the H 2 + is significantly lighter than the S + fragment, it carries away a large part of the energy release, and the low S + velocity leads to increased uncertainty in the scattering angles and velocity of the S + product. Fig. 3 shows that the velocity of the S + product ion is broadly oriented with the velocity of the incident dication, w(S 2+ ). Conversely, the scattering of the H 2 + product ion is broadly oriented with w(H 2 ), anti-parallel to w(S 2+ ). This form of scattering was also noted in the S 2+ + Ar SET reaction, and, as mentioned in the discussion of that channel, such kinematics are consistent with direct electron-transfer occurring at a significant interspecies separation (3-6 Å). Such scattering patterns are commonly observed for other NDSET processes and are generally well represented by a LZ formalism. 69,99,112,113 The experimental exoergicity distribution of this NDSET channel, shown in Fig. 4, is centred at B6.7 eV, with a full width at half maximum (FWHM) from 4.8-8.6 eV. To rationalise the experimental exoergicity spectrum in this NDSET reaction we must again consider the possible electronic states of the reactant and product species. For this collision system, energetic data are readily available. Given the evidence from the S 2+ + Ar reaction, and the arguments presented in the experimental section, we would expect the reactions to involve S 2+ in the 3s 2 3p 2 states. 103 The H 2 , emitted as an effusive beam, will be in its ground vibronic state,   are seen to be significantly less favoured. 131 Therefore, the involvement of excited electronic H 2 + states in this channel is unlikely. From energetic considerations, we see that there are many possible reaction pathways that result in literature exoergicities that fall within the range of the observed experimental exoergicities. The reactant dication, S 2+ , could be in any of the S 2+ 3s 2 3p 2 electronic states ( 3 P, 1 D and 1 S). The H 2 + product ion can be formed in ground state, X 2 S g + , with a wide range of vibrational excitation: with an energy of anywhere from 15.4 eV (X 2 S g + v = 0) to the dissociation limit at B18.1 eV relative to H 2 (X 1 S g + ).
The S + product ion can be formed in its three lowest energy states: 4 S, 2 D and 2 P. Previous studies have shown that the vibrational state distribution of product molecular monocations, formed in dication-neutral SET reactions, are dominated by vertical (Frank-Condon style) transitions from the relevant reactant (dication or neutral molecule). 93 Photoelectron spectra of H 2 show that H 2 + is primarily populated from H 2 with an energy of oB17.5 eV with respect to the neutral ground state, and it is likely that the H 2 + formed from H 2 in this NDSET channel will also be formed in this energy range. 121,122,126,127 Given this H 2 + energy range, and the likely dominance of S 2+ ( 3 P) states in the beam, we propose that pathways (e) and (f) involving the 3 P state are primarily involved in this NDSET channel. The listed literature exoergicities for these pathways assume the population of H 2 + ( 2 S g + ) v = 0. However, as noted above, it is likely that a spread of H 2 + vibrational states will be formed. The population of excited vibrational levels of H 2 + will, of course, spread the observed experimental exoergicity to lower values. The literature exoergicities resulting from pathways (e) and (f) are marked on Fig. 4, with arrows showing the energy spread due to the potential vibrational excitation of H 2 + , and provide a good match to the observed experimental exoergicity distribution (Fig. 4). Of course, analogous channels associated with the other states ( 1 D, 1 S) from the 3p 2 configuration of S 2+ , which we expect to be minor components in the S 2+ beam are also likely to contribute here. Indeed, analogous channels to (e) and (f) for these S 2+ excited states nicely account for the lower intensity signals up to 11.3 eV in the experimental exoergicity spectrum. However, as for the reaction with Ar, we can easily account for the bulk of the observed exoergicities in terms of the S 2+ states we argue dominate the dication beam populating readily accessible states of the product ions. The dynamics we extract from each pair of ions making up the S + + H + peak in the coincidence spectrum reveal that the reaction proceeds via DSET. The established mechanism for dication + neutral DSET reactions is that an initial LZ style electron-transfer generates one (or both) of the product monocations in a dissociative state. 15 The nascent dissociative states subsequently fragment to yield the observed products. Fig. 5a shows the CM scattering diagram for the S + and H + products from this DSET channel. As with the NDSET channel discussed above, it should be noted that the relatively low S + velocity leads to increased uncertainty in the scattering angles of the S + . The S + is broadly forward scattered but with a tail towards higher angles when compared with the NDSET reaction (Fig. 3), scattering mirrored in the H + signals. This scattering pattern is indicative of a longer-lived association between the reactant species than in the standard NDSET process: a [SH 2 ] 2+ collision complex. Therefore, this DSET reaction appears not to follow the standard model of dication + neutral electron-transfer where an initial LZ style (long-distance) electron-transfer populates one of the products in a dissociative electronic state. 95,132,133 As discussed, such an LZ style electron-transfer reaction should lead to a more prominent forward scattering pattern than is observed here. Indeed, although not common, SET involving complexation has been observed before in dication/neutral reactions. 90,95,132,133 Electron ionization experiments studying H 2 S show evidence for H 2 S 2+ states that are energetically accessible to the S 2+ + H 2 reactants in this system, and that dissociate to form S + + H 2 + . 100 These H 2 S 2+ states could be involved as a short-lived intermediates in this DSET channel. Fig. 5b shows the scattering of the H + and H fragments relative to the velocity of the S + product ion. We refer the reader back to Section 3 for information on the construction of these internal frame scattering diagrams. The H + and H fragments are both clearly backward scattered, away from the S + ion and their scattering is distributed about a clear H þÃ 2 precursor velocity which is indicated in Fig. 5b. This scattering motif, which has been observed before for many DSET reactions, 112,114,115,134 clearly shows the [SH 2 ] 2+ complex dissociates into S + + H þÃ 2 before the fragmentation of H þÃ 2 . The H + fragment ion distribution is slightly more energetically backward scattered than that of the H fragment, likely because the H þÃ 2 ion fragments within the electric field of the S + ion, resulting in the H + experiencing an additional acceleration.
The experimental exoergicity observed from this DSET channel is shown in Fig. 6. The experimental exoergicity distribution has a broad maximum, centred at B10 eV, with a FWHM from 6.5-15.5 eV. In order to give the exoergicities observed in this channel, the S 2+ dication must be in an electronically excited state, thus the dominant species we expect in the dication beam S 2+ ( 3 P) cannot contribute and the experimental exoergicity spectrum reveals the excited S 2+ states in the beam. Even the coupling of the collision energy (0.71 eV) into the energetics does not allow reactions of S 2+ ( 3 P) to generate the observed exoergicities. Representative reactions of the relevant S 2+ excited states (g)-(j) which fit the observed experimental exoergicity distribution are given below. Variations on these processes are also compatible with the exoergicity data ( Fig. 6) involving, for example, the formation of excited states of S + . The literature exoergicities resulting from pathways (g)-(j) are marked on the experimental exoergicity spectrum in Fig. 6. There is an exoergic channel for the 1 D state (Fig. 6, (g)), but the maximum exoergicity possible for the highest energy state of the p 2 configuration ( 1 S) is 8.6 eV (marked by (h) in Fig. 6) and there are clearly events with experimental exoergicities markedly above this value; such exoergicities can only be generated by electronic states of S 2+ from excited configurations. That is, the demanding energetics allow the contribution of these minor beam components to be distinguished from that of the p 2 configuration. Above we considered a minor beam component as comprising the metastable 5 S and 3 F states. Reaction of S 2+ ( 5 S) can give exoergicities up to B12.5 eV, and the involvement of 3 F can give exoergicities of up to B20.4 eV. These exoergicities nicely encompass the higher end of the observed experimental exoergicity distribution (Fig. 6, (i)-(j)) and indeed the reaction of the 3 F state correlates nicely with a weak high exoergicity feature in the spectrum (Fig. 6, (j)).
Despite the assignments of the product states not being definitive, it is clear that the DSET reaction with H 2 involves the excited states of S 2+ we predicted to be present in the dication beam. The contribution of these excited states dominates the product flux in this reaction channel as the demanding energetics mean the dominant dication state in the ion beam, the 3 P ground state of S 2+ , cannot contribute. Given our deductions concerning the dominance of the 3 P state in the beam, the involvement of excited S 2+ in this channel explains the low signal rate noted above. Applying a powerful approach we have used in previous work to study the energetics of DSET reactions we can, on an event-by-event basis, determine the velocity of the precursor H þÃ 2 ion from the measured velocity of the S + ion using conservation of momentum. The S + and H þÃ 2 velocities then allow us to construct an exoergicity spectrum for the initial electron transfer step to reveal the electronic states involved. Subsequently, again on an event by event basis we can examine the velocities of the H and H + products in the frame of the H þÃ 2 precursor velocity to derive a spectrum representing the energies involved in the dissociation of the H þÃ 2 ion. This analysis assumes the formation of S + + H þÃ 2 can be treated as a two-body reaction, dynamically distinct from the subsequent dissociation of the H þÃ 2 . Given the observed scattering dynamics (Fig. 5b) this is clearly a good, but not perfect approximation. Applying this approach yields the experimental exoergicity distribution for dissociation of H þÃ 2 to H + + H shown in Fig. 7.
In 1986, Cordaro et al. 125 studied the dissociation of H 2 + generated by electron ionization, observing experimental exoergicities from 2-8 eV. Strathdee and Browning 135 observed similar kinetic energy releases, ranging from 4-10 eV, using 26.9 eV photoionization. In their study of H 2 electron ionization, Köllmann 136 reported H + kinetic energy distributions that correspond to kinetic energy releases ranging from 0-24 eV resulting from H 2 + dissociation. The electronic states attributed to these kinetic energy releases were the high vibrational states of H 2 + (X 2 S g + ) and H 2 + (A 2 S u + ). Other experiments studying photoand electron-ionization of H 2 also resulted in the observation of kinetic energy releases from H 2 + dissociation of up to B20 eV. [137][138][139] The bulk of the experimental exoergicity distribution resulting from H 2 + dissociation observed in our experiment (Fig. 7) matches well with Cordaro et al. 125 and Strathdee and Browning, 135 and the shoulder towards higher experimental exoergicities is in agreement with the distributions observed by Köllmann and others. [136][137][138][139] The agreement between the exoergicity release resulting from H 2 + dissociation observed in our experiment and those previously reported confirms that the mechanism responsible for the formation of S + + H + in this channel is DSET via H þÃ 2 . As mentioned, if H 2 + is formed at energies over B18.1 eV above the ground state of the neutral H 2 molecule, it will dissociate. 124,125 Our analysis of the experimental exoergicity ( Fig. 6) above suggests that only the lowest energy H state ( 2 S) is formed in this dissociation. To yield experimental exoergicities in the range observed in Figure 7, H 2 + therefore needs to be primarily formed with energies between 19-24 eV relative to H 2 (X 1 S g + ). Thus, the H 2 + ion could be populated in the A 2 S u + state. As noted above, the geometry of the A 2 S u + state in the 19-24 eV energy range is significantly different to the neutral molecule and would therefore not be expected to be efficiently formed in a vertical ionizing transition. Moreover, the formation of H 2 + (A 2 S u + ) from H 2 (X 1 S g + ) is a two-electron transition, which, as discussed above, is thought to be significantly less likely than a one-electron process. 131 However, both forming the H 2 + ion with a markedly different geometry to the neutral species, and the movement of two electrons are more feasible if the reaction proceeds via a collision complex, as suggested by the scattering we observe. Alternatively, H 2 + could be formed in the continuum of the X 2 S g + state, above the 18.1 eV dissociation asymptote. The involement of either the H 2 + (A 2 S u + ) or H 2 + (X 2 S g + ) states is in good agreement with the H 2 + dissociation studies discussed above. 125,136 As discussed above, the exoergicity of the initial electrontransfer step of this DSET reaction can also be evaluated using the precursor velocity of the H 2 + ion and w(S + ). Using this method, an estimate of the experimental exoergicity distribution for the initial electron-transfer step (S 2+ + H 2 -S + + H þÃ 2 ) was determined and is shown in Fig. 8. The experimental exoergicity distribution has a broad peak centred at B7 eV with a FWHM from 1-12 eV. The observed experimental exoergicity distribution exceeds the typical range for LZ style processes (2-7 eV), 112,114 providing further indirect evidence that this DSET reaction proceeds via a collision complex, as suggested by the scattering.
The experimental exoergicity distribution determined for the initial electron-transfer step (Fig. 8) aligns satisfyingly with the literature exoergicities expected from our above assessment of the likely electronic states involved in this channel. These include the excited states of S 2+ ( 1 D, 1 S, 5 S and 3 F), the lowest three energy states of S + ( 4 S, 2 D and 2 P), and H 2 + (A 2 S u + ) or the continuum of the H 2 + (X 2 S g + ) state. The possible involvement of this large number of states, and the range of vibrational excitation available for H 2 + , results in too many accessible pathways to determine which specific states are involved. In summary, DSET occurs following the collisions of S 2+ + H 2 , resulting in the formation of S + + H + (+ H). This channel involves S 2+ in its electronically excited states derived from the 3s 2 3p 2 electron configuration ( 1 D, 1 S) and also shows evidence for the presence of higher energy S 2+ states in the beam: 5 S and 3 F, resulting from the 3s 1 3p 3 and 3s 2 3p 1 3d 1 configurations respectively. This DSET channel involves the formation of S + primarily in its lowest electronic state ( 4 S), with contributions from the first two excited states, 2 D and 2 P. The nascent H 2 + ion is generated at initial energies of between 19 eV and 24 eV relative to H 2 (X 1 S g + ) before dissociating to the lowest energy asymptote: H + + H( 2 S).  The collisions of S 2+ + N 2 PSCO-MS spectra were recorded following the collisions of S 2+ and N 2 at E cm = 4.7 eV. In the coincidence pairs spectrum three reaction channels which produce pairs of product ions were observed, and are shown in Table 2. The dominant channel, Rxn III, is NDSET, producing S + + N 2 + . A DSET reaction (Rxn IV), results in the formation of S + + N + + N. Finally, a bond-forming channel (Rxn V), producing SN + + N + , is detected. Whilst the PSCO-MS apparatus is optimised for the detection of pairs of ions, we also record a simple mass spectrum of single ion arrivals. We do not detect any structures present in this mass spectrum ions in the m/z = 23 region that could be indicative of SN 2+ or the m/z = 30 region that could be indicative of SN 2 2+ . However, we note that the sensitivity of our experiment is optimised for detecting ion pairs involving energetic monocations, not product dications. It is notable that there is no peak in the coincidence spectrum corresponding to N + + N + , the products of a dissociative DET reaction. However, the absence of a DET reaction is not surprising as the polarisability of S (2.9 Â 10 À24 cm 3 ) is greater than N 2 (1.7403 Â 10 À24 cm 3 ). 140 In this situation, under the concerted model of DET, 115,133 we would not expect the S 2+ + N 2 and S + N 2 2+ curves to cross even if the S 2+ + N 2 -S + N 2 2+ reaction was exothermic; this situation has been shown before to strongly disfavour DET. Added to this, DET between S 2+ and N 2 would be endothermic for S 2+ states with energies below S 2+ ( 3 F), the states that we have deduced make up the bulk of our S 2+ beam.
Non-dissociative single electron-transfer (NDSET) Fig. 9 shows the CM scattering diagram for the S + + N 2 + product ions of the NDSET reaction, S 2+ + N 2 -S + + N 2 + . A strong forward scattering pattern is observed, where the S + is oriented in the same direction as w(S 2+ ), whilst the N 2 + product is oriented in the same direction as w(N 2 ). Note the peaking of the signals in Fig. 9 slightly away from y = 01 and 1801 as discussed in the experimental section. This scattering pattern is typical of a direct electron-transfer process, occurring at a relatively large interspecies separation (3-6 Å) as discussed above. These electron-transfer dynamics are well represented by a Landau-Zener formalism. 69,99,112,113 Fig. 10 shows a histogram for the experimental exoergicities recorded in this NDSET reaction channel. The spectrum shows a structure with a maximum at 5.6 eV, and a FWHM from 4.3-6.9 eV. The observed experimental exoergicities are typical for such a LZ style SET reaction. 112,114 As discussed above, we have good evidence that the dication beam is primarily comprised of S 2+ in the 3s 2 3p 2 configuration, dominated by the 3 P state. The N 2 reactant will be in its ground vibronic state, X 1 S g + v = 0. The ground state of N 2 + , X 2 S g + , has an energy of 15.58 eV relative to N 2 (X 1 S g + ). 141,142 If N 2 + is formed with energies above B24.3 eV relative to N 2 (X 1 S g + ) (corresponding to N 2 + (C 2 S u + v = 3)), it will be unstable to dissociation to the N + ( 3 P) + N( 4 S) asymptote. 103,143 Photoionization studies have shown that N 2 + states generated with more than 24.3 eV, relative to N 2 (X 1 S g + ), have dissociation lifetimes less than the timescale of our experiment. Therefore, we do not expect any contribution to the N 2 + counts observed in this channel from N 2 + states formed with more than 24.3 eV of energy with respect to N 2 . [143][144][145][146] Electron-transfer occurring between N 2 and any of the 3 P, 1 D and 1 S electronic states of S 2+ can result in the formation of S + in any of the states from the ground state ( 4 S) up to the 4 F state, derived from the [Ne] 3s 2 3p 2 3d 1 configuration. a The signal-to-noise ratio in Rxn IV is low and therefore it is difficult to extract meaningful energetic information from this channel.  Considering the above constraints, we find that there are many possible electronic pathways that could be involved in the formation of S + + N 2 + . However, as discussed above we expect the major component of our S 2+ beam to be ions in the 3 P state. Possible pathways involving the 3 P state are shown by (k)-(n), and their resulting literature exoergicities are shown on Fig. 10. These channels show N 2 + is likely formed in one of the stable X 2 S g + , A 2 P u or B 2 S u + states, which reassuringly are those which dominate photoelectron spectra, 141 and the S + products are primarily the lowest two electronic states: 4 S and 2 D. Although these contributions from the dication 3 P state (pathways (k)-(n)) can account for the signals we observe, especially when varying degrees of vibrational excitation of the product N 2 + ion is considered which will broaden the observed exoergicities, excited S 2+ states could also be contributing to the product flux in this NDSET channel. However, given energetically favourable channels for the reaction of the 3 P state, which lie in the reaction window, and the expected dominance of this state in the ion beam, we feel confident that pathways (k)-(n) account for the bulk of the events we detect in this reactive channel. The peak in the coincidence pairs spectrum corresponding to the detection of ions with m/z = 32 and m/z = 14 is due to the formation of S + and N + . The peak is weak, so it is only possible to probe the dynamics and exoergicities in an indicative manner. In order to discriminate against some of the stray counts arising from background noise, which become problematic in this channel due to the low number of counts, any events that gave unrealistic experimental exoergicities were removed before the analysis. Fig. 11a shows the CM scattering diagram for the S + and N + products resulting from this channel. The S + ions are broadly forward scattered, however, there is scattering to higher angles. Conversely, the N + ions are broadly backscattered, with some ions scattered towards lower angles. One might expect a higher propensity of forward scattering if this channel resulted from a direct, LZ style electron-transfer. The presence of scattering to higher angles could be evidence of the involvement of a collision complex, however, due to the low intensity signal in this channel, it is difficult to be conclusive. Fig. 11b shows the internal frame scattering of the N + and N products relative to the S + product ion. The N + and N ions are scattered around the expected precursor velocity of an N þÃ 2 intermediate, confirming the mechanism of this channel as DSET. SET occurs, producing N þÃ 2 in a dissociative state (and S + ), which then dissociates to form N + + N. In Fig. 11b the N + is scattered away from the S + with a higher velocity than the N fragment. This asymmetric scattering of the N + and N fragments is because the N þÃ 2 species dissociates within the electric field of the S + ion, and Coulomb repulsion then accelerates only the N + species. Such behaviour has been noted before for DSET channels in other dicationic collision systems. 90,115,133 The experimental exoergicity distribution resulting from this DSET channel reveals a broad structure between 0.5 eV and 11.0 eV. Note that due to the low intensity, the energy resolution is low, and hence any conclusions we draw from this data must be tentative. It is not common to see energy releases of less than 2 eV for reactions involving LZ style ET, however, as discussed above, the scattering hints that a more complicated reaction mechanism is involved. The formation of N + + S + + N resulting from the ground state of S 2+ ( 3 P) is endothermic, therefore the occurrence of this channel again confirms the existence of electronically excited states of S 2+ in the dication beam. Indeed, the exclusive involvement of higher energy S 2+ states undoubtedly accounts for the low intensity of this DSET channel. The involvement of S 2+ ( 1 S) can result in exoergicities in this DSET channel of up to 2.4 eV, significantly lower than the bulk of the observed experimental exoergicity distribution. However, the involvement of the 5 S and 3 F states of S 2+ allows events with exoergicities of up to 6.3 eV and 14.1 eV respectively. Given the high experimental exoergicities observed in this DSET channel, the S + , N + and N fragments are likely all formed in their ground or lower energy excited electronic states.
The experimental exoergicity of the N þÃ 2 dissociation, calculated from the N 2 + precursor velocity, has structure from 0.3-4.0 eV. This range of experimental exoergicities is comparable to those observed from N 2 + dissociation in the DSET channel in our previous investigation of Ar 2+ + N 2 collisions. 90 As mentioned above, N 2 + states generated with more than 24.3 eV above N 2 (X 1 S g + ) are unstable to dissociation and will dissociate within the lifetime of our experiment. Considering the experimental exoergicity of the complete reaction and the N 2 + dissociation exoergicity, N 2 + is likely formed with energies of 25-29 eV above N 2 (X 1 S g + ) and primarily dissociates to the lowest energy N + + N asymptote, with both N + and N in their ground electronic states.

Chemical bond-formation
The peak in the pairs spectrum corresponding to the ions m/z = 46 and m/z = 14 is due to the detection of SN + and N + : a bondforming reaction (Rxn V). The SN + molecular ion is relevant to a range of astrophysical environments, including in comets and nebulae. 147-151 SN + is also important terrestrially, in man-made plasmas. 152 Fig. 12 shows the CM scattering for the SN + and N + product ions in this channel, revealing an isotropic scattering pattern. Such scattering is strong evidence for the involvement of an intermediate, long-lived, collision complex, [SN 2 ] 2+ . The collision complex exists for at least several rotations, before fragmenting into SN + and N + , as shown in (o).
The experimental exoergicity of Rxn V, determined from the product ion velocities, is shown in Fig. 13. The experimental exoergicity distribution has a maximum at 2.0 eV with a FWHM from 0.5-4.5 eV. In order to rationalise the observed experimental exoergicities, we must again consider the possible electronic states of the reactant and product species. Due to its astrophysical relevance, the SN + species has been the subject of experimental study, 150,153,154 and several theoretical investigations have probed the electronic states of SN + . 151,[155][156][157] In this work we use the electronic state structure calculated by Ben Yaghlane & Hochlaf (2009), 151 as it provides the most comprehensive data. The SN + ground state (X 1 S + ) lies B6.3 eV below the lowest energy dissociation asymptote of S + ( 4 S) + N( 4 S), and B8.3 eV below the second lowest energy dissociation asymptote of S + ( 2 D) + N( 4 S). There are also many other bound SN + electronic states lying below these dissociation asymptotes, several with well depths sufficient to sustain vibrational excitation. Again, we assume that the dication beam is primarily comprised of the ground configuration, dominated by the 3 P state, and N 2 is in its ground vibronic state ( 1 S g + ). These energetics mean that N + can therefore be formed in any of its four lowest energy electronic states. Considering the above, there are a large number of pathways that match the observed experimental exoergicity distribution of this bond-forming channel. Two likely pathways which match well with the exoergicity spectrum and involve the ground state of S 2+ , the dominant ion in our beam, are shown by  The 1 D and 1 S states of S 2+ could also be involved, as well as the first three excited states of N + ( 1 D, 1 S and 5 S). Moreover, it is possible that a number of the electronic states of SN + that are stable to dissociation could be populated in this channel.
To summarise, this channel involves a chemical bondforming reaction between S 2+ + N 2 , producing SN + + N + . An isotropic scattering pattern is observed, clearly showing the involvement of a collision complex, [SN 2 ] 2+ , typical for a dicationic bond-forming reaction. The observed experimental exoergicity distribution is in good agreement with the population of known electronic states of SN + , 151 and the involvement of the S 2+ states that we know dominate our S 2+ beam.

Conclusions
In this work we have studied the gas-phase reactivity of S 2+ with Ar, N 2 and H 2 . This, to our knowledge, is the first study of the bimolecular reactivity of S 2+ . The experiments indicate that the dication beam is composed predominantly of S 2+ ions resulting from the 3s 2 3p 2 ground electronic configuration. Most of the observed reaction channels can be accounted for by the involvement of S 2+ in its ground state, 3 P. However, the dissociative electron transfer reactivity reveals the involvement of higher energy S 2+ states.
The SET reaction between S 2+ and Ar displays dynamics and energetics typical of a Landau-Zener style electron transfer reaction. In the reactions between S 2+ + H 2 , NDSET and DSET were observed. The dynamics here show evidence for the formation of a long-lived association between the reactants, a collision complex [SH 2 ] 2+ . Perhaps surprisingly, the formation of SH + , the potential product of a chemical bond-forming channel, was not observed. The reactions of S 2+ + H 2 are especially relevant as S 2+ has been observed in the outer atmosphere of Jupiter, where H 2 is the dominant chemical species. The collisions of S 2+ + N 2 result in both dissociative and nondissociative SET reactions as well as a bond-forming channel, the latter resulting in the formation of SN + + N + . This bondforming reaction proceeds via a collision complex.

Conflicts of interest
There are no conflicts to declare.