Direct production of molecular oxygen from carbon dioxide and helium ion collisions

The prebiotic mechanism to produce molecular oxygen (O2) in carbon dioxide (CO2)-rich planetary atmospheres is of great importance in understanding astrochemical reactions and is potentially relevant to the origin of life on Earth. Here, we demonstrate that, aside from the direct productions of O2 by photodissociation and dissociative electron attachment, the low-energy ion-molecule reaction between cationic helium in solar winds and molecular CO2 is a noticeable mechanism. Branching ratios of the reaction channels are determined, and their absolute cross-sections are estimated accordingly. The present findings represent a further, indispensable step towards fully understanding the origins of atmospheric O2.

by the measurements of time-of-flight (TOF) mass spectra and ion velocity maps using a crossed-beam apparatus [21][22][23] .
9][20] ), while those technical troubles are settled in the present measurements.On the other hand, the translational energy of He + in the Martian atmosphere spreads over a wide range and depends on its sources, for instance, the velocity of He + is comparable to that of He 2+ in the solar wind if He + is produced by the charge exchange reaction between He 2+ and neutral hydrogen 13 .Besides the thermalenergy collisions, He + + CO 2 reactions at the collision energy of electron volts are particularly interesting, since the products in different electronic states, as shown in Fig. 1, play versatile roles in the subsequent processes.
CO 2 + formed in the resonant charge exchange between CO 2 and He + can be populated in the ro-vibrationally excited states of electronic state C 2Σ g + or the 2 Π u -satellite states, but its spontaneous dissociation is highly preferred due to an electron depletion from the bonding orbital 4σ g or 1π u of CO 2 (refs. 24,25).20]26 ), and as one of their subsequent processes, the photoemission A ͂ 2 Π → X ͂ 2 Σ + of the CO + yield may contribute to the comet tail glow 26 .Such scenarios can be further found in Fig. 1.For instance, the phosphorescence from the O 2 in singlet states (a 1 Δ g , b 1 Σ g + ) was observed in the aurora 27 .The reactions noted in Fig. 1, once authenticated experimentally, should be of fundamental importance toward understanding the oxygen-related processes in planetary atmospheres.
The present experiments were conducted with our home-made crossed-beam apparatus which was recently updated by introducing the three-dimensional velocity map imaging technique of ionic products 22,23 .The TOF mass spectra and velocity distributions of all ionic yields can be obtained simultaneously.The collection efficiency of ionic yields is close to 100%, and each peak in the TOF mass spectrum corresponds to the whole Newton spheres of a certain type of ionic yield.The multiple collisions in the ion flowing or the ion drift of cyclotron resonance [18][19][20] are absent or reduced significantly in the crossed-beam arrangement where the pulsed reactant ion beam is perpendicular to a supersonic target beam.
Due to the energy redistribution in a collision, the endothermic channels (their energetic thresholds can be found in Table S1 of Supplementary Note 1) in Fig. 1 are possibly accessed.For example, around 3.5 eV of the collision energy, the reactions leading to O 2 (X ͂ 3 Σ g ¯) + C + ( 4 P 1/2 ) and 2O( 3 P 2 ) + C + ( 2 P 1/2 ) are allowed in energetics if the collisional or translational energy is completely transformed into the internal energies of the target.However, in the TOF mass spectrometry measurement, the same product, C + as mentioned above, cannot be distinguished from which channel the ion originates.Different ion production channels could be further identified by the ion velocity imaging measurements.In this work, the combinational measurements were accomplished at three collision energies in the reaction center-of-mass (c.m.) coordinate (E c.m. = 1.94, 2.49, and 4.75 eV).

Results and discussion
The TOF mass spectrum at E c.m. = 4.75 eV, as a typical example, is shown in Fig. 2a, in which the C + and O 2 + ions are observed clearly but their relative intensities are much lower than those of other ionic yields.Such low signals of C + and O 2 + certainly brought difficulties in the previous measurements [18][19][20]26 . Her, we further obtained the branching ratios of the reactions and plotted their collision energy dependences in Fig. 2b.A comparison between the present branching ratios and the previous ones [18][19][20] is presented in Table S2 of Supplementary Note 2.
In the reactions at the higher collision energies (10-10 3 eV), the production ratio γ (C + /CO + ) was found to gradually increase with the enhancement of He + kinetic energy 24 .Here one can find that the branching ratios of the reactions denoted with Eqs. ( 4) and ( 5) decrease with the E c.m. enhancement.Such different collision-energy dependences should be attributed to reaction dynamics, e.g., diversities of collisional trajectories and dissociation pathways.Usually, the time scale of a charge exchange at the higher collision energy could be shorter than that of molecular vibrational motions.O 2 can only be formed via the slow atomic roaming 6 or the bending motion with a large amplitude of CO 2 (ref. 7).Thereby, similar scenarios should exist for the O 2 or O 2 + product of the He + + CO 2 reactions at the lower collision energy.
The velocity maps of C + and O 2 + products in Fig. 3 correspond to the ionic three-dimensional Newton spheres projected on the detector or reaction plane, rather than the central slices of the spheres.Nevertheless, we still can find the collision-energy dependences of the velocity distributions.Figure 3a indicates that the C + yields are primarily located around the target CO 2 at the lower collision energies, but most of them are distributed in the forward direction (namely the CO 2 flying direction) at E c.m. = 4.75 eV.In contrast to the C + velocity maps, as shown in Fig. 3b, most of the O 2 + ions are always distributed around CO 2 and become crowded at E c.m. = 4.75 eV.The forward scattering preference and the product distributions around the reactant CO 2 indicate that the reactions would like to proceed in the largeparameter collisions, namely, the collisional or translational energy is hardly transferred to the target in such a contactless reaction.Therefore, the higher-threshold channel leading to the fragments 2 O( 3 P 2 ) + C + ( 2 P 1/2 ) is unfavorable at E c.m. = 4.75 eV, due to the inefficient translational-to-internal energy transformation.The isotropic angular distributions (red circles in Fig. 3) represent the C + yield of the CO 2 + dissociation where CO 2 + is assumed to be rebounded toward the c.m. position and receives 4.75 eV in a head-on collision.What we observed here obviously deviates from this hypothetical model.
The O 2 coproduct is likely in b 1 Σ g + at the lower collision energies, primarily due to near energy-resonance of the channel to O 2 (b 1 Σ g + ) + C + ( 2 P 1/2 ) (see Fig. 1) and partially according to the velocity u(C + ) distributions in Fig. S1a.The O 2 coproduct could be in X 3Σ g ¯or A ̃3Σ u + state at the higher collision energy 4.75 eV.The u(O 2 + ) profiles in Fig. S1b, by contrast, exhibit a preference for the ground-state product O 2 + (X 2Π g ).O 2 or O 2 + in other states could be produced as the minor yields.More importantly, the absolute cross sections or production efficiencies of O 2 and O 2 + are the fundamental data of atmospheric reaction networks, but very few are available so far [18][19][20]24,26,28,29 . Extrapolatig the total cross-sections of the He + + CO 2 reaction (all channels) from the higher collision energies 28,29 to the low-energy range and according to the branching ratios obtained in this work, we further derive the absolute cross-sections of about 10 −22 m 2 of the reactions denoted with Eqs. ( 4) and ( 5) in a collision energy range of 1-5 eV (more details can be found in Supplementary Note 3).
1 Asymptotic energies of three charge-exchange reactions between He + and CO 2 .The product energy is scaled, concerning the zero value for the resonant charge exchange of the motionless reagents He + ( 2 S) + CO 2 (X 1Σ g ).Besides the exothermic pathways, the endothermic processes could proceed if the collision energy is availably transformed.These cross-sections are a few hundred times larger than that of the direct production of O 2 via DEA 7 .
As discussed previously, oxygen-atom roaming and molecular bending are two feasible pathways to produce O 2 from CO 2 , where the former usually results in the nearly isotropic angular distributions of the products while a distinct anisotropy can be observed for the latter 6,7 .As proposed above, in the largeparameter collisional reaction, the He atom formed by the prompt charge exchange should be a spectator in the subsequent dissociation.To further reveal this unimolecular-like process, we make a coordinate transformation from the reaction c.m. one (used in Fig. 3) to the CO 2 reaction coordinate (more details can be found in Supplementary Note 4).As illustrated in Fig. S2, the CO 2 reaction coordinate is also defined with respect to the collision axis, but its scattering angle is different from that in the reaction c.m. coordinate.Figure 4 shows the angular and velocity distributions of C + and O 2 + in the CO 2 reaction coordinate.
Figure 4a exhibits the nearly isotropic distributions of C + at the lower E c.m. values and a clearly forward preference at 4.75 eV, where O 2 and C + yields are assumed to be in specific states.As indicated in Fig. S3, this collision-energy dependence scarcely relies on the products' states.A cartoon in the bottom panel of Fig. S3b illustrates a possible molecular orientation, namely, the dissociation undergoes the bended conformation which is almost perpendicular to the collision axis.It is reasonable that such a stereodynamics feature becomes indistinct in the slower collisional reaction.All O 2 + angular distributions in Fig. 4c are almost isotropic, which could be attributed to the roaming dynamics, namely, the roaming O + abstracts the oxygen atom from the CO moiety.The velocities of C + and O 2 + in Fig. 4b and d indicate the lower values at the higher collision energy, due to an enhancement of the products in the higher states.
To further understand the dissociation dynamics of CO 2 + , we plot the potential energy surfaces of C ̃2Σ g + state in Fig. S4 of Supplementary Note 5. Different shapes or patterns of above potential energy surfaces imply that the O 2 + or O 2 formation hardly undergoes a combinational motion of molecular bending and symmetric bond stretching, but possibly benefits from an asymmetric stretching.The asymmetric stretching should be also responsible for the predominance of the CO + and O + products observed in this work and the previous studies [18][19][20]24,26,28,29 . However, itwould be a tough task to simulate the dynamic trajectories leading to O 2 or O 2 + because of their much lower yields.Furthermore, the nonadiabatic process, for instance, coupling between C ̃2Σ g + state and 2 Π u -satellite states of CO 2 +25 , may be involved in the reactions.

Conclusions
Based on the previous and present findings, we make a conclusive remark on the origins of atmospheric O 2 .Since the production efficiency of the photodissociation of CO 2 was not high enough to reproduce the measured data 30 , the other production sources of those O 2 molecules must be considered.The He + + CO 2 reaction, particularly in the low collision energy range, is an important mechanism but was ignored before this study.Up to date, we have known different mechanisms for the origin of atmospheric O 2 .The O atoms produced in the ion-molecule reaction, DEA, and photodissociation can also participate in the formation of O 2 .Some contributions could be from the O + and O 2 + ions by the electron capture, however, these ions are likely responsible for the oxygen escaping 16,17 .All-inclusive networks of the atmospheric and interstellar oxygen-related reactions are expected to be established by more and more laboratory studies and astronomical observations.Methods Experiments.Our three-dimensional velocity map imaging apparatus has been described elsewhere 22 , in which the product ion detector consists of a set of multichannel plates and a Delay-Line-Detector (DLD ϕ = 80 mm, commercially available from Roentdek).The ion source that previously produced Ar + ions 21 was updated recently to produce lighter ions (such as He + ).In the present experiment, the energy spreads (ΔE/E ~7%) of the He + bunches are comparable to that of the Ar + bunches.The x−y positions (the plane of the detector) and the time of flight (along z direction) of each ionic signal are recorded simultaneously with the DLD, and then a Newton sphere of a certain type of ionic yield is rebuilt with these data 22,23 .In Fig. 3, the twodimensional images are plotted with the raw data of whole Newton spheres of the ionic products projected on the detector, where one pixel corresponds to ca. 200 m/s of C + or O 2 + .On the other hand, the TOF mass spectrum can be directly obtained from the measurement of whole Newton spheres, thus the ion collection efficiency is nearly 100%.High-purity (99.99%) samples (He and CO 2 ) are available commercially and used without further purification.The He + reactant ions are produced in the electron impacts at 30 eV.The kinetic energies and the energy spreads of the reactants He + and CO 2 (its ΔE/E ~6%) were calibrated or determined prior to the experiments.During the experiments, the reaction chamber is evacuated under a steady vacuum condition of 2.0 × 10 −7 mbar and the reactions happen in the field-free region.The working frequency is 10 kHz, and we need dozens of hours to record the velocity images at a given collision energy.
Calculations.Equilibrium geometry and vibrational frequencies of the ground-state CO 2 were calculated at the density functional theory level wb97xd/def2TZVP.This molecular geometry was used as the initial point in scanning the potential energy surfaces along the vibrational motions, i.e., symmetric stretching Q1, bending Q2, and asymmetric stretching Q3.Table S4 shows the corresponding values of real bond length and bond angle for each Q1-Q3.The equation of motion method based on the coupled cluster method limited to singles and doubles excitations and the basis set cc-pVTZ [31][32][33][34] (encoded in Guassion16 suit of programs 35 ) were used here to obtain the two-dimensional potential energy surfaces of CO 2 + (C 2 Σ g + ).

Fig. 2
Fig. 2 Ionic yields of the collisional reaction between He + and CO 2 at 4.75 eV (a) and branching ratios (b). a Time-of-flight (TOF) mass spectrum is plotted with the real ion counts.b Branching ratios (with the statistic uncertainties of the measurements) of O 2 + and C + yields are amplified to be more visible.

Fig. 3 Fig. 4
Fig. 3 Velocity map images of C + (a) and O 2 + (b).Each velocity image is the projection of three-dimensional Newton spheres of C + or O 2 + to the ion detector.Red circles represent the isotropic distributions for O 2 (X 3Σ g ¯) + C +(2 P 1/2 ) or O 2 +( X 2Π g ) + C( 3 P 0 ) (outer circle) and O 2 (X 3Σ g ¯) + C + ( 4 P 1/2 ) or O 2 + (a 4 Π u ) + C( 3 P 0 ) (inner circle) yields that are assumed to be produced in photodissociation of CO 2 + at the center-of-mass (c.m.) of the collision.White lines with arrows show the coordinate transformation from the laboratory (the reactants' velocities are represented with v, broken lines) and the collisional c.m. coordinate (the reactants' velocities are represented with u, solid lines).The ionic signals in the left and right sides respective to c.m. are defined as the forward and backward scatterings.