Ab initio dynamics and photoionization mass spectrometry reveal ion-molecule pathways from ionized acetylene clusters to benzene cation

Symmetry-breaking dynamics of a photoionized carbon dioxide dimer

Photoionization can initiate structural reorganization of molecular matter and drive formation of new chemical bonds. Here, we used time-resolved extreme ultraviolet (EUV) pump - EUV probe Coulomb explosion imaging of carbon dioxide dimer ionC O 2 2 + dynamics, that combined with ab initio molecular dynamics simulations, revealed unexpected asymmetric structural rearrangement. We show that ionization by the pump pulse induces rearrangement from the slipped-parallel (C2h) geometry of the neutral C O 2 dimer towards a T-shaped (C2v) structure on the ~100 fs timescale, although the most stable slipped-parallel (C2h) structure of the ionic dimer. Moreover, we find that excited states of the ionized C O 2 dimer can exhibit formation of aCO 3 moiety in theC 2 O 4 + complex that can persist even after a suitably time-delayed second photoionization in a metastableC 2 O 4 2 + dication. Our results suggest that charge asymmetry plays an important role in the ionization-induced dynamics in such dimers that are present in C O 2 rich environments.

Magic Numbers and Stabilities of Photoionized Water Clusters: Computational and Experimental Characterization of the Nanosolvated Hydronium Ion

The stability and distributions of small water clusters generated in a supersonic beam expansion are interrogated by tunable vacuum ultraviolet (VUV) radiation generated at a synchrotron. Time-of-flight mass spectrometry reveals enhanced population of various protonated water clusters (H⁺(H₂O)_n) based upon ionization energy and photoionization distance from source, suggesting there are "magic" numbers below the traditional n = 21 that predominates in the literature. These intensity distributions suggest that VUV threshold photoionization (11.0-11.5 eV) of neutral water clusters close to the nozzle exit leads to a different nonequilibrium state compared to a skimmed molecular beam. This results in the appearance of a new magic number at 14. Metadynamics conformer searches coupled with modern density functional calculations are used to identify the global minimum energy structures of protonated water clusters between n = 2 and 21, as well as the manifold of low-lying metastable minima. New lowest energy structures are reported for the cases of n = 5, 6, 11, 12, 16, and 18, and special stability is identified by several measures. These theoretical results are in agreement with the experiments performed in this work in that n = 14 is shown to exhibit additional stability, based on the computed second-order stabilization energy relative to most cluster sizes, though not to the extent of the well-known n = 21 cluster. Other cluster sizes that show some additional energetic stability are n = 7, 9, 12, 17, and 19. To gain insight into the balance between ion-water and water-water interactions as a function of the cluster size, an analysis of the effective two-body interactions (which sum exactly to the total interaction energy) was performed. This analysis reveals a crossover as a function of cluster size between a water-hydronium-dominated regime for small clusters and a water-water-dominated regime for larger clusters around n = 17.

Exploring three-body fragmentation of acetylene trication

The three-body breakup of [C₂H₂]³⁺ formed in collision with Xe⁹⁺ moving at 0.5 atomic units of velocity is studied by using recoil ion momentum spectroscopy. Three-body breakup channels leading to (H⁺, C⁺, CH⁺) and (H⁺, H⁺, C₂ ⁺) fragments are observed in the experiment and their kinetic energy release is measured. The breakup into (H⁺, C⁺, CH⁺) occurs via concerted and sequential modes, whereas the breakup into (H⁺, H⁺, C₂ ⁺) shows only the concerted mode. By collecting events coming exclusively from the sequential breakup leading to (H⁺, C⁺, CH⁺), we have determined the kinetic energy release for the unimolecular fragmentation of the molecular intermediate, [C₂H]²⁺. By using ab initio calculations, the potential energy surface for the lowest electronic state of [C₂H]²⁺ is generated, which shows the existence of a metastable state with two possible dissociation pathways. A discussion on the agreement between our experimental results and these ab initio calculations is presented.

Methane Cluster Fragmentation by Fast Electron Impact

We investigate the fragmentation of the CH4 cluster by fast electron impact at stagnation pressures from 0.5 bar to 16 bar. By measuring the time of flight spectrum (TOF), two types of ions, including (CH4)n−1CH5+ and (CH4)n−2(C2Hm)+, are observed. In the 1D TOF spectrum, it is shown that for the stagnation pressure larger than 4 bar, the former ion is predominant for each n, similar to the previous experimental result. However, as the pressure decreases to 0.5 or 2 bar, the contribution of the C2Hm+ ion is dominant over that of the CH4CH5+ ion. In the 2D coincident TOF spectrum, the above two patterns of ions are also distinguished, and the enhancement of C2Hm+ is observed at 4 bar pressure. The phenomena appearing in 2D and 1D TOF spectra imply that the C2Hm+ ion prefers to survive in a smaller cluster, while the stabilization of the protonated ion needs a more massive cluster environment.

A Computational and Experimental View of Hydrogen Bonding in Glycerol Water Clusters

Polyol-water clusters provide a template to probe ionization and solvation processes of paramount interest in atmospheric and interstellar chemistry. We generate glycerol water clusters in a continuous supersonic jet expansion and interrogate the neutral species with synchrotron-based tunable vacuum ultraviolet photoionization mass spectrometry. A series of glycerol fragments (m/z 44, 61, 62, and 74) clustered with water are observed. A judicious combination of backing pressure, nozzle temperature, and water vapor pressure allows for tuning the mol % of glycerol. The recorded appearance energies of the water cluster series m/z 62 and 74 are similar to that observed in pure glycerol, while the m/z 61 series shows a dependence on cluster composition. Furthermore, this series also tracks the water concentration of the beam. Theoretical calculations on neutral and ionized clusters visualize the hydrogen bond network in these water clusters and provide an assessment of the number of glycerol-glycerol, glycerol-water, and water-water hydrogen bonds in the cluster, as well as their interaction energies. This method of bond counting and interaction energy assessment explains the changes in the mass spectrum as a function of mol % and offers a glimpse of the disruption of the hydrogen bond network in glycerol-water clusters. The calculations also reveal interesting barrierless chemical processes in photoionized glycerol water clusters that are either activated or do not occur without the presence of water. Examples include spontaneous intramolecular proton transfer within glycerol to form a distonic ion, nonactivated breaking of a C-C bond, and spontaneous proton transfer from glycerol to water. These results appear relevant to radiation-induced chemical processing of alcohol-water ices in the interstellar medium.

Excited state photochemically driven surface formation of benzene from acetylene ices on Pluto and in the outer solar system

NASA's New Horizons mission unveiled a diverse landscape of Pluto's surface with massive regions being neutral in color, while others like Cthulhu Macula range from golden-yellow to reddish comprising up to half of Pluto's carbon budget. Here, we demonstrate in laboratory experiments merged with electronic structure calculations that the photolysis of solid acetylene - the most abundant precipitate on Pluto's surface - by low energy ultraviolet photons efficiently synthesizes benzene and polycyclic aromatic hydrocarbons via excited state photochemistry thus providing critical molecular building blocks for the colored surface material. Since low energy photons deliver doses to Pluto's surface exceeding those from cosmic rays by six orders of magnitude, these processes may significantly contribute to the coloration of Pluto's surface and of hydrocarbon-covered surfaces of Solar System bodies such as Triton in general. This discovery critically enhances our perception of the distribution of aromatic molecules and carbon throughout our Solar System.

Cyclotrimerization of Acetylene under Thermal Conditions: Gas-Phase Kinetics of V+ and Fe+ + C2H2

The kinetics of successive reactions of acetylene (C₂H₂) initiated on either vanadium or iron atomic cations have been investigated under thermal conditions using the variable-ion source and temperature-adjustable selected-ion flow tube apparatus. Consistent with the literature results, the reaction of Fe⁺ + C₂H₂ primarily yields Fe⁺(m/z = (C₂H₂)₃); however, analysis via quantum chemical calculations and statistical modeling shows that the mechanism does not form benzene upon the third acetylene addition. The kinetics are more consistent with successive addition of three acetylene molecules, yielding Fe⁺(C₂H₂)₃, followed by an addition of a fourth acetylene molecule, initiating cyclotrimerization, yielding either Fe⁺(C₂H₂) + neutral benzene or Fe⁺(Bz) + acetylene, where Bz is a benzene ligand. In contrast, the reaction of V⁺ + C₂H₂ yields products via successive associations V⁺(m/z = (C₂H₂)_n) either with or without a bimolecular step involving loss of one H₂ and V⁺C₂(m/z = (C₂H₂)_m), where n and m extend at least up to 11 under conditions of 0.32 Torr at 300 K. Stabilized V⁺(Bz) is not a significant intermediate in the association mechanism. We propose a plausible mechanism for the generation of neutral benzene in this reaction and compare with the Fe⁺ results. The reaction steps that produce benzene result in turnover of the single-atom catalyst, and the large hydrocarbons produced that remain associated to the catalyst are proposed to be polycyclic aromatic hydrocarbons.

The Effect of Cluster Size on the Intra-Cluster Ionic Polymerization Process

Polyaromatic hydrocarbons (PAHs) are widespread in the interstellar medium (ISM). The abundance and relevance of PAHs call for a clear understanding of their formation mechanisms, which, to date, have not been completely deciphered. Of particular interest is the formation of benzene, the basic building block of PAHs. It has been shown that the ionization of neutral clusters can lead to an intra-cluster ionic polymerization process that results in molecular growth. Ab-initio molecular dynamics (AIMD) studies in clusters consisting of 3-6 units of acetylene modeling ionization events under ISM conditions have shown maximum aggregation of three acetylene molecules forming bonded C₆H₆⁺ species; the larger the number of acetylene molecules, the higher the production of C₆H₆⁺. These results lead to the question of whether clusters larger than those studied thus far promote aggregation beyond three acetylene units and whether larger clusters can result in higher C₆H₆⁺ production. In this study, we report results from AIMD simulations modeling the ionization of 10 and 20 acetylene clusters. The simulations show aggregation of up to four acetylene units producing bonded C₈H₈⁺. Interestingly, C₈H₈⁺ bicyclic species were identified, setting a precedent for their astrochemical identification. Comparable reactivity rates were shown with 10 and 20 acetylene clusters.

Molecular dynamics reveals formation path of benzonitrile and other molecules in conditions relevant to the interstellar medium

Polycyclic aromatic hydrocarbons and polycyclic aromatic nitrogen heterocycles are believed to be widespread in different areas of the interstellar medium. However, the astronomical detection of specific aromatic molecules is extremely challenging. As a result, only a few aromatic molecules have been identified, and very little is known about how they are formed in different areas of the interstellar medium. Recently, McGuire et al. [Science 359, 202-205 (2018)] detected the simple aromatic molecule benzonitrile in Taurus Molecular Cloud-1. Although benzonitrile has been observed, the molecular mechanism for its formation is still unknown. In this study, we use quantum chemistry and ab initio molecular dynamics to model ionization processes of van der Waals clusters containing cyanoacetylene and acetylene molecules. We demonstrate computationally that the clusters' ionization leads to molecular formation. For pure cyanoacetylene clusters, we observe bond formation among two and three monomer units, whereas in mixed clusters, bond formation can also occur in up to four units. We show that the large amount of energy available to the system after ionization can lead to barrier crossing and the formation of complex molecules. Our study reveals the rich chemistry that is observed upon ionization of the clusters, with a wide variety of molecules being formed. Benzonitrile is among the observed molecules, and we study the potential energy path for its formation. These results also offer insights that can guide astronomers in their search for aromatic molecules in the interstellar medium.

Formation of covalently bound C4H4 + upon electron-impact ionization of acetylene dimer

We investigate the formation mechanisms of covalently bound C₄H₄ ⁺ cations from direct ionization of hydrogen bonded dimers of acetylene molecules through fragment ion and electron coincident momentum spectroscopy and quantum chemistry calculations. The measurements of momenta and energies of two outgoing electrons and one ion in triple-coincidence allow us to assign the ionization channels associated with different ionic fragments. The measured binding energy spectra show that the formation of C₄H₄ ⁺ can be attributed to the ionization of the outermost 1π_u orbital of acetylene. The kinetic energy distributions of the ionic fragments indicate that the C₄H₄ ⁺ ions originate from direct ionization of acetylene dimers while ions resulting from the fragmentation of larger clusters would obtain significantly larger momenta. The formation of C₄H₄ ⁺ through the evaporation mechanism in larger clusters is not identified in the present experiments. The calculated potential energy curves show a potential well for the electronic ground state of (C₂H₂)2⁺, supporting that the ionization of (C₂H₂)₂ dimers can form stable C₂H₂⋅C₂H₂ ⁺(1π_u ^-1) cations. Further transition state analysis and ab initio molecular dynamics simulations reveal a detailed picture of the formation dynamics. After ionization of (C₂H₂)₂, the system undergoes a significant rearrangement of the structure involving, in particular, C-C bond formation and hydrogen migrations, leading to different C₄4⁺ isomers.

Molecular growth upon ionization of van der Waals clusters containing HCCH and HCN is a pathway to prebiotic molecules

The growth mechanisms of organic molecules in an ionizing environment such as the interstellar medium are not completely understood. Here we examine by means of ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) computations the possibility of bond formation and molecular growth upon ionization of van der Waals clusters of pure HCN clusters, and mixed clusters of HCN and HCCH, both of which are widespread in the interstellar medium. Ionization of van der Waals clusters can potentially lead to growth in low temperature and low-density environments. Our results show, that upon ionization of the pure HCN clusters, strongly bound stable structures are formed that contain NH bonds, and growth beyond pairwise HCN molecules is seen only in a small percentage of cases. In contrast, mixed clusters, where HCCH is preferentially ionized over HCN, can grow up to 3 or 4 units long with new carbon-carbon and carbon-nitrogen covalent bonds. Moreover, cyclic molecules formed, such as the radical cation of pyridine, which is a prebiotic molecule. The results presented here are significant as they provide a feasible pathway for molecular growth of small organic molecules containing both carbon and nitrogen in cold and relatively denser environments such as in dense molecular clouds but closer to the photo-dissociation regions, and protoplanetary disks. In the mechanism we propose, first, a neutral van der Waals cluster is formed. Once the cluster is formed it can undergo photoionization which leads to chemical reactivity without any reaction barrier.

Exciton energy transfer reveals spectral signatures of excited states in clusters

Electronic excitation and concomitant energy transfer leading to Penning ionization in argon-acetylene clusters generated in a supersonic expansion are investigated with synchrotron-based photoionization mass spectrometry and electronic structure calculations. Spectral features in the photoionization efficiency of the mixed argon-acetylene clusters reveal a blue shift from the 2P1/2 and 2P3/2 excited states of atomic argon. Analysis of this feature suggests that excited states of argon clusters transfer energy to acetylene, resulting in its ionization and successive evaporation of argon. Theoretically calculated Arn (n = 2-6) cluster spectra are in excellent agreement with experimental observations, and provide insight into the structure and ionization dynamics of the clusters. A comparison between argon-acetylene and argon-water clusters reveals that argon solvates water better, allowing for higher-order excitons and Rydberg states to be populated. These results are explained by theoretical calculations of respective binding energies and structures.

Penning spectroscopy and structure of acetylene oligomers in He nanodroplets

Embedded atoms or molecules in a photoexcited He nanodroplet are well-known to be ionized through inter-atomic relaxation in a Penning process. In this work, we investigate the Penning ionization of acetylene oligomers occurring from the photoexcitation bands of He nanodroplets. In close analogy to conventional Penning electron spectroscopy by thermal atomic collisions, the n = 2 photoexcitation band plays the role of the metastable atomic 1s2s 3,1S He*. This facilitates electron spectroscopy of acetylene aggregates in the sub-Kelvin He environment, providing the following insight into their structure: the molecules in the dopant cluster are loosely bound van der Waals complexes rather than forming covalent compounds. In addition, this work reveals a Penning process stemming from the n = 4 band where charge-transfer from autoionized He in the droplets is known to be the dominant relaxation channel. This allows for excited states of the remnant dopant oligomer Penning-ions to be studied. Hence, we demonstrate Penning ionization electron spectroscopy of doped droplets as an effective technique for investigating dopant oligomers which are easily formed by attachment to the host cluster.

Quantum Chemical Exploration of Intermolecular Reactions of Acetylene

Quantum chemical explorations of potential energy surfaces showed that acetylene produces various products starting from molecular arrays in short distances of 1.3-2.5 Å. Arrays of (C₂ H₂ )₂ gave cyclobutadiene, tetrahedrane, and acetylene dimers. Arrays of (C₂ H₂ )₃ gave benzene, prismane, benzvalene, Dewar benzene, and acetylene trimers. Arrays of (C₂ H₂ )₄ gave cubane, cyclooctatetranene, and acetylene tetramers. Different forms of initial arrays yielded different sets of products; a parallel array of two monomers gave cyclobutadiene, whereas a cross array gave tetrahedrane. Initial molecular arrays with unusually close contacts were estimated to require local forces of 1-9 × 10^-8  N. © 2019 Wiley Periodicals, Inc.

From atoms to aerosols: probing clusters and nanoparticles with synchrotron based mass spectrometry and X-ray spectroscopy

Synchrotron radiation provides insight into spectroscopy and dynamics in clusters and nanoparticles.

Development of an experimental set-up for low-temperature spectroscopic studies of matrix-isolated molecules and molecular ices using synchrotron radiation

An experimental set-up for studying photophysics and photochemistry of molecules in an inert gas medium (matrix-isolated) and in the ice phase at low temperatures has been developed and commissioned at the Photophysics beamline, Indus-1 synchrotron radiation source. This end-station uses an in-house-developed closed-cycle cryostat for achieving cryo-temperatures (∼10 K). Synchrotron radiation from the Photophysics beamline is used as the source of UV-VUV photons and the system is equipped with a Fourier transform infrared spectrometer for characterization of the molecular species formed at low temperature. Various individual components of the end-station like closed-cycle cryostat, experimental chamber, gas mixing and deposition systems are tested to ascertain that the desired performance criteria are satisfied. The performance of the composite system after integration with the Photophysics beamline is evaluated by recording IR and UV-VUV photoabsorption spectra of sulfur dioxide at low temperatures (10 K), both in the ice phase as well as isolated in argon matrices. Results obtained are in good agreement with earlier literature, thus validating the satisfactory performance of the system. As an off-shoot of the study, the VUV absorption spectrum of matrix-isolated SO₂ in argon matrix up to 10.2 eV is reported here for the first time. This experimental end-station will provide new opportunities to study photon-induced reactions in molecules of environmental, astrochemical and industrial importance. Details of the design, development and initial experimental results obtained are presented.