A crossed molecular beam investigation of the N(2D) + pyridine reaction and implications for prebiotic chemistry

Vacuum-ultraviolet irradiation of pyridine:acetylene ices relevant to Titan astrochemistry

Nitrogen-containing polycyclic aromatic hydrocarbons (NPAHs) are important molecules for astrochemistry and prebiotic chemistry, as their occurrence spans from interstellar molecular clouds to planetary systems. Their formation has been previously explored in gas phase experiments, but the role of solid-state chemical reactions in their formation under cryogenic conditions remains elusive. Here, we explore the formation of NPAHs through vacuum ultraviolet (VUV) irradiation of pyridine:acetylene ices in amorphous and co-crystalline phases, with the aim to simulate conditions relevant to the interstellar medium and Titan's atmosphere. Our results show that the synthesis of ethynylpyridines from VUV-irradiated pyridine:acetylene amorphous ices is achievable at 18 K. In the co-crystal phase, photolysis at 110 K leads to the formation of NPAHs such as quinolizinium+ and precursors, reflecting a dynamical system under our conditions. In contrast, irradiation at 90 K under stable conditions did not produce volatile photoproducts. These results suggest that such chemical processes can occur in Titan's atmosphere and potentially in its stratosphere, where the co-condensation of these molecules can form composite ices. Concurrently, the formation of stable co-crystals can influence the depletion rates of pyridine, which suggests that these structures can be preserved and potentially delivered to Titan's surface. Our findings provide insights into the molecular diversity and chemical evolution of organic matter on Titan, crucial for future space exploration missions, such as the Dragonfly mission, which may uncover higher-order organics derived from pyridine precursors on Titan's surface.

Reactant Discovery with an Ab Initio Nanoreactor: Exploration of Astrophysical N-Heterocycle Precursors and Formation Pathways

The incorporation of nitrogen atoms into cyclic compounds is essential for terrestrial life; nitrogen-containing (N-)heterocycles make up DNA and RNA nucleobases, several amino acids, B vitamins, porphyrins, and other components of biomolecules. The discovery of these molecules on meteorites with non-terrestrial isotopic abundances supports the hypothesis of exogenous delivery of prebiotic material to early Earth; however, there has been no detection of these species in interstellar environments, indicating that there is a need for greater knowledge of their astrochemical formation and destruction pathways. Here, we present results of simulations of gas-phase pyrrole and pyridine formation from an ab initio nanoreactor, a first-principles molecular dynamics simulation method that accelerates reaction discovery by applying non-equilibrium forces that are agnostic to individual reaction coordinates. Using the nanoreactor in a retrosynthetic mode, starting with the N-heterocycle of interest and a radical leaving group, then considering the discovered reaction pathways in reverse, a rich landscape of N-heterocycle-forming reactivity can be found. Several of these reaction pathways, when mapped to their corresponding minimum energy paths, correspond to novel barrierless formation pathways for pyridine and pyrrole, starting from both detected and hypothesized astrochemical precursors. This study demonstrates how first-principles reaction discovery can build mechanistic knowledge in astrochemical environments as well as in early Earth models such as Titan's atmosphere where N-heterocycles have been tentatively detected.

Unveiling the Reaction Mechanism of the N(2D) + Pyridine Reaction: Ring-Contraction versus 7-Membered-Ring Formation Channels

Despite the relevance of the reactions of the prototypical nitrogen-containing six-membered aromatic molecule (N-heterocyclic) of pyridine (C6H5N) in environmental science, astrochemistry, planetary science, prebiotic chemistry, and materials science, few experimental/theoretical studies exist on the bimolecular reactions involving pyridine and neutral atomic/molecular radicals. We report a combined experimental and theoretical study on the elementary reaction of pyridine with excited nitrogen atoms, N(2D), aimed at providing information about the primary reaction products and their branching fractions (BFs). From previous crossed molecular beam (CMB) experiments with mass-spectrometric detection and present synergistic calculations of the reactive potential energy surface (PES) and product BFs we have unveiled the reaction mechanism. It is found that the reaction proceeds via N(2D) barrierless addition to pyridine that, via bridged intermediates followed by N atom "sliding" into the ring, leads to 7-membered-ring structures. They further evolve, mainly via ring-contraction mechanisms toward 5-membered-ring radical products and, to a smaller extent, via H-displacement mechanisms toward 7-membered-ring isomeric products and their isomers. Using the theoretical statistical estimates, an improved fit of the experimental data previously reported has been obtained, leading to the following results for the dominant product channels: C4H4N (pyrrolyl) + HCN (BF = 0.61 ± 0.20), C3H3N2 (1H-imidazolyl/1H-pyrazolyl) + C2H2 (BF = 0.11 ± 0.06), and C5H4N2 (7-membered-ring molecules or pyrrole carbonitriles) + H (BF = 0.28 ± 0.10). The ring-contraction product channels C4H4N (pyrrolyl) + HCN, C3H3N2 (1H-imidazolyl) + C2H2, C3H3N2 (1H-pyrazolyl) + C2H2, and C5H5 (cyclopentadienyl) + N2 have statistical BFs of 0.54, 0.09, 0.11, and 0.07, respectively. Among the H-displacement channels, the cyclic-CHCHCHCHNCN + H channel and cyclic-CHCHCHCHCN2 + H are theoretically predicted to have a comparable BF (0.07 and 0.06, respectively), while the other isomeric 7-membered-ring molecule + H channel has a BF of 0.03. Pyrrole-carbonitriles and 1H-ethynyl-1H-imidazole (+ H) isomeric channels have an overall BF of 0.03. Implications for the chemistry of Saturn's moon Titan and prebiotic chemistry, as well as for understanding the N-doping of graphene or carbon nanotubes, are noted.

Synthesis and characterization of azaborepin radicals in solid neon through boron-mediated C-N bond cleavage of pyridine

The reactivity of pyridine is a complex topic due to its unique electronic structure. The reactions of atomic boron with pyridine molecules in solid neon have been investigated using matrix isolation infrared absorption spectroscopy. Three products (marked as A, B, and C) were observed and characterized through 10B, D and 15N isotopic substitution pyridine regents as well as quantum chemical calculations. In the reaction, the ground-state boron atom can attack the lone pair electrons of the nitrogen atom in the pyridine molecule, resulting in the formation of a 1-boropyridinyl radical (A). Alternatively, addition to the aromatic π-system of pyridine can occur in a [1,4] type, leading to the formation of a B[η2(1,4)-C5H5N] complex (B). Under UV-visible light (280 < λ < 580 nm) irradiation, these two compounds can further undergo photo-isomerization to form BN-embedded seven-membered azaborepin compounds (C). The observation of species A, B, and the subsequent photo-isomerization to species C is consistent with theoretical predictions, indicating that these reactions are kinetically favorable. This research provides valuable insights into the future design and synthesis of corresponding BN heterocyclic derivatives.

Nonadiabatic quantum dynamics explores non-monotonic photodissociation branching of N2 into the N(4S) + N(2D) and N(4S) + N(2P) product channels

Vacuum ultraviolet (VUV) photodissociation of N2 molecules is a source of reactive N atoms in the interstellar medium. In the energy range of VUV optical excitation of N2, the N-N triple bond cleavage leads to three types of atoms: ground-state N(4S) and excited-state N(2P) and N(2D). The latter is the highest reactive and it is believed to be the primary participant in reactions with hydrocarbons in Titan's atmosphere. Experimental studies have observed a non-monotonic energy dependence and non-statistical character of the photodissociation of N2. This implies different dissociation pathways and final atomic products for different wavelength regions in the sunlight spectrum. We here apply ab initio quantum chemical and nonadiabatic quantum dynamical techniques to follow the path of an electronic state from the excitation of a particular singlet 1Σ+u and 1Πu vibronic level of N2 to its dissociation into different atomic products. We simulate dynamics for two isotopomers of the nitrogen molecule, 14N2 and 14N15N for which experimental data on the branching are available. Our computations capture the non-monotonic energy dependence of the photodissociation branching ratios in the energy range 108 000-116 000 cm-1. Tracing the quantum dynamics in a bunch of electronic states enables us to identify the key components that determine the efficacy of singlet to triplet population transfer and therefore predissociation lifetimes and branching ratios for different energy regions.

An experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction with implications for the photochemical models of Titan

We report on a combined experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction, which is of relevance in the aromatic chemistry of the atmosphere of Titan. Experimentally, the reaction was studied (i) under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (Ec) of 31.8 kJ mol-1 to determine the primary products, their branching fractions (BFs), and the reaction micromechanism, and (ii) in a continuous supersonic flow reactor to determine the rate constant as a function of temperature from 50 K to 296 K. Theoretically, electronic structure calculations of the doublet C6H6N potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction mechanism. The reaction is found to proceed via barrierless addition of N(2D) to the aromatic ring of C6H6, followed by formation of several cyclic (five-, six-, and seven-membered ring) and linear isomeric C6H6N intermediates that can undergo unimolecular decomposition to bimolecular products. Statistical estimates of product BFs on the theoretical PES were carried out under the conditions of the CMB experiments and at the temperatures relevant for Titan's atmosphere. In all conditions the ring-contraction channel leading to C5H5 (cyclopentadienyl) + HCN is dominant, while minor contributions come from the channels leading to o-C6H5N (o-N-cycloheptatriene radical) + H, C4H4N (pyrrolyl) + C2H2 (acetylene), C5H5CN (cyano-cyclopentadiene) + H, and p-C6H5N + H. Rate constants (which are close to the gas kinetic limit at all temperatures, with the recommended value of 2.19 ± 0.30 × 10-10 cm3 s-1 over the 50-296 K range) and BFs have been used in a photochemical model of Titan's atmosphere to simulate the effect of the title reaction on the species abundances as a function of the altitude.

Reaction N(2D) + CH2CCH2 (Allene): An Experimental and Theoretical Investigation and Implications for the Photochemical Models of Titan

We report on a combined experimental and theoretical investigation of the N(2D) + CH2CCH2 (allene) reaction of relevance in the atmospheric chemistry of Titan. Experimentally, the reaction was investigated (i) under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (E c) of 33 kJ/mol to determine the primary products and the reaction micromechanism and (ii) in a continuous supersonic flow reactor to determine the rate constant as a function of temperature from 50 to 296 K. Theoretically, electronic structure calculations of the doublet C3H4N potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction mechanism. The reaction is found to proceed via barrierless addition of N(2D) to one of the two equivalent carbon-carbon double bonds of CH2CCH2, followed by the formation of several cyclic and linear isomeric C3H4N intermediates that can undergo unimolecular decomposition to bimolecular products with elimination of H, CH3, HCN, HNC, and CN. The kinetic experiments confirm the barrierless nature of the reaction through the measurement of rate constants close to the gas-kinetic rate at all temperatures. Statistical estimates of product branching fractions (BFs) on the theoretical PES were carried out under the conditions of the CMB experiments at room temperature and at temperatures (94 and 175 K) relevant for Titan. Up to 14 competing product channels were statistically predicted with the main ones at E c = 33 kJ/mol being formation of cyclic-CH2C(N)CH + H (BF = 87.0%) followed by CHCCHNH + H (BF = 10.5%) and CH2CCNH + H (BF = 1.4%) the other 11 possible channels being negligible (BFs ranging from 0 to 0.5%). BFs under the other conditions are essentially unchanged. Experimental dynamical information could only be obtained on the overall H-displacement channel, while other possible channels could not be confirmed within the sensitivity of the method. This is also in line with theoretical predictions as the other possible channels are predicted to be negligible, including the HCN/HNC + C2H3 (vinyl) channels (overall BF < 1%). The dynamics and product distributions are dramatically different with respect to those observed in the isomeric reaction N(2D) + CH3CCH (propyne), where at a similar E c the main product channels are CH2NH (methanimine) + C2H (BF = 41%), c-C(N)CH + CH3 (BF = 32%), and CH2CHCN (vinyl cyanide) + H (BF = 12%). Rate coefficients (the recommended value is 1.7 (±0.2) × 10-10 cm3 s-1 over the 50-300 K range) and BFs have been used in a photochemical model of Titan's atmosphere to simulate the effect of the title reaction on the species abundance (including any new products formed) as a function of the altitude.

The N(2D) + CH2CHCN (Vinyl Cyanide) Reaction: A Combined Crossed Molecular Beam and Theoretical Study and Implications for the Atmosphere of Titan

The reaction of electronically excited nitrogen atoms, N(²D), with vinyl cyanide, CH₂CHCN, has been investigated under single-collision conditions by the crossed molecular beam (CMB) scattering method with mass spectrometric detection and time-of-flight (TOF) analysis at the collision energy, E_c, of 31.4 kJ/mol. Synergistic electronic structure calculations of the doublet potential energy surface (PES) have been performed to assist in the interpretation of the experimental results and characterize the overall reaction micromechanism. Statistical (Rice–Ramsperger–Kassel–Marcus, RRKM) calculations of product branching fractions (BFs) on the theoretical PES have been carried out at different values of temperature, including the one corresponding to the temperature (175 K) of Titan’s stratosphere and at a total energy corresponding to the E_c of the CMB experiment. According to our theoretical calculations, the reaction is found to proceed via barrierless addition of N(²D) to the carbon–carbon double bond of CH₂=CH–CN, followed by the formation of cyclic and linear intermediates that can undergo H, CN, and HCN elimination. In competition, the N(²D) addition to the CN group is also possible via a submerged barrier, leading ultimately to N₂ + C₃H₃ formation, the most exothermic of all possible channels. Product angular and TOF distributions have been recorded for the H-displacement channels leading to the formation of a variety of possible C₃H₂N₂ isomeric products. Experimentally, no evidence of CN, HCN, and N₂ forming channels was observed. These findings were corroborated by the theory, which predicts a variety of competing product channels, following N(²D) addition to the double bond, with the main ones, at E_c = 31.4 kJ/mol, being six isomeric H forming channels: c-CH(N)CHCN + H (BF = 35.0%), c-CHNCHCN + H (BF = 28.1%), CH₂NCCN + H (BF = 26.3%), c-CH₂(N)CCN(cyano-azirine) + H (BF = 7.4%), trans-HNCCHCN + H (BF = 1.6%), and cis-HNCCHCN + H (BF = 1.3%), while C–C bond breaking channels leading to c-CH₂(N)CH(2H-azirine) + CN and c-CH₂(N)C + HCN are predicted to be negligible (0.02% and 0.2%, respectively). The highly exothermic N₂ + CH₂CCH (propargyl) channel is also predicted to be negligible because of the very high isomerization barrier from the initial addition intermediate to the precursor intermediate able to lead to products. The predicted product BFs are found to have, in general, a very weak energy dependence. The above cyclic and linear products containing an additional C–N bond could be potential precursors of more complex, N-rich organic molecules that contribute to the formation of the aerosols on Titan’s upper atmosphere. Overall, the results are expected to have a significant impact on the gas-phase chemistry of Titan’s atmosphere and should be properly included in the photochemical models.

Rotational and Vibrational Spectra of the Pyridyl Radicals: A Coupled-Cluster Study

Pyridyl is a prototypical nitrogen-containing aromatic radical that may be a key intermediate in the formation of nitrogen-containing aromatic molecules under astrophysical conditions. On meteorites, a variety of complex molecules with nitrogen-containing rings have been detected with nonterrestrial isotopic abundances, and larger nitrogen-containing polycyclic aromatic hydrocarbons (PANHs) have been proposed to be responsible for certain unidentified infrared emission bands in the interstellar medium. In this work, the three isomers of pyridyl (2-, 3-, and 4-pyridyl) have been investigated with coupled cluster methods. For each species, structures were optimized at the CCSD(T)/cc-pwCVTZ level of theory and force fields were calculated at the CCSD(T)/ANO0 level of theory. Second-order vibrational perturbation theory (VPT2) was used to derive anharmonic vibrational frequencies and vibrationally corrected rotational constants, and resonances among vibrational states below 3500 cm^-1 were treated variationally with the VPT2+K method. The results yield a complete set of spectroscopic parameters needed to simulate the pure rotational spectrum of each isomer, including electron-spin, spin-spin, and nuclear hyperfine interactions, and the calculated hyperfine parameters agree well with the limited available data from electron paramagnetic resonance spectroscopy. For the handful of experimentally measured vibrational frequencies determined from photoelectron spectroscopy and matrix isolation spectroscopy, the typical agreement is comparable to experimental uncertainty. The predicted parameters for rotational spectroscopy reported here can guide new experimental investigations into the yet-unobserved rotational spectra of these radicals.

Neutral Dissociation of Pyridine Evoked by Irradiation of Ionized Atomic and Molecular Hydrogen Beams

The interactions of ions with molecules and the determination of their dissociation patterns are challenging endeavors of fundamental importance for theoretical and experimental science. In particular, the investigations on bond-breaking and new bond-forming processes triggered by the ionic impact may shed light on the stellar wind interaction with interstellar media, ionic beam irradiations of the living cells, ion-track nanotechnology, radiation hardness analysis of materials, and focused ion beam etching, deposition, and lithography. Due to its vital role in the natural environment, the pyridine molecule has become the subject of both basic and applied research in recent years. Therefore, dissociation of the gas phase pyridine (C5H5N) into neutral excited atomic and molecular fragments following protons (H+) and dihydrogen cations (H2+) impact has been investigated experimentally in the 5-1000 eV energy range. The collision-induced emission spectroscopy has been exploited to detect luminescence in the wavelength range from 190 to 520 nm at the different kinetic energies of both cations. High-resolution optical fragmentation spectra reveal emission bands due to the CH(A2Δ→X2Πr; B2Σ+→X2Πr; C2Σ+→X2Πr) and CN(B2Σ+→X2Σ+) transitions as well as atomic H and C lines. Their spectral line shapes and qualitative band intensities are examined in detail. The analysis shows that the H2+ irradiation enhances pyridine ring fragmentation and creates various fragments more pronounced than H+ cations. The plausible collisional processes and fragmentation pathways leading to the identified products are discussed and compared with the latest results obtained in cation-induced fragmentation of pyridine.

The Reaction N(2D) + CH3CCH (Methylacetylene): A Combined Crossed Molecular Beams and Theoretical Investigation and Implications for the Atmosphere of Titan

The reaction of excited nitrogen atoms N(²D) with CH₃CCH (methylacetylene) was investigated under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (E_c) of 31.0 kJ/mol. Synergistic electronic structure calculations of the doublet potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction micromechanism. Theoretically, the reaction is found to proceed via a barrierless addition of N(²D) to the carbon–carbon triple bond of CH₃CCH and an insertion of N(²D) into the CH bond of the methyl group, followed by the formation of cyclic and linear intermediates that can undergo H, CH₃, and C₂H elimination or isomerize to other intermediates before unimolecularly decaying to a variety of products. Kinetic calculations for addition and insertion mechanisms and statistical (Rice-Ramsperger-Kassel-Marcus) computations of product branching fractions (BFs) on the theoretical PES were performed at different values of total energy, including the one corresponding to the temperature (175 K) of Titan’s stratosphere and that of the CMB experiment. Up to 14 competing product channels were statistically predicted, with the main ones, at E_c = 31.0 kJ/mol, being the formation of CH₂NH (methanimine) + C₂H (ethylidyne) (BF = 0.41), c-C(N)CH + CH₃ (BF = 0.32), CH₂CHCN (acrylonitrile) + H (BF = 0.12), and c-CH₂C(N)CH + H (BF = 0.04). Of the 14 possible channels, seven correspond to H displacement channels of different exothermicity, for a total H channel BF of ∼0.25 at E_c = 31.0 kJ/mol. Experimentally, dynamical information could only be obtained about the overall H channels. In particular, the experiment corroborates the formation of acrylonitrile + H, which is the most exothermic of all 14 reaction channels and is theoretically calculated to be the dominant H-forming channel (BF = 0.12). The products containing a novel C–N bond could be potential precursors to form other nitriles (C₂N₂, C₃N) or more complex organic species containing N atoms in planetary atmospheres, such as those of Titan and Pluto. Overall, the results are expected to have a potentially significant impact on the understanding of the gas-phase chemistry of Titan’s atmosphere and the modeling of that atmosphere.