Diversity of protonated mixed pyrene-water clusters investigated by collision induced dissociation.

Protonated mixed pyrene-water clusters, (Py)m(H2O)nH+, where m = [1-3] and n = [1-10], are generated using a cryogenic molecular cluster source. Subsequently, the mass-selected mixed clusters undergo controlled collisions with rare gases, and the resulting fragmentation mass spectra are meticulously analyzed to discern distinct fragmentation channels. Notably, protonated water cluster fragments emerge for n ≥ 3, whereas they are absent for n = 1 and 2. The experimental results are complemented by theoretical calculations of structures and energetics for (Py)(H2O)nH+ with n = [1-4]. These calculations reveal a shift in proton localization, transitioning from the pyrene molecule for n = 1 and 2 to water molecules for n ≥ 3. The results support a formation scenario wherein water molecules attach to protonated pyrene PyH+ seeds, and, by extension, to (Py)2H+ and (Py)3H+ seeds. Various isomers are identified, corresponding to potential protonation sites on the pyrene molecule. Protonated polycyclic aromatic hydrocarbons are likely to be formed in cold, dense interstellar clouds and protoplanetary disks due to the high proton affinity of these species. Our findings show that the presence of protonated PAHs in these environments could lead to the formation of water clusters and mixed carbon-water nanograins, having a potential impact on the water cycle in regions of planet formation.

Fingerprinting fragments of fragile interstellar molecules: dissociation chemistry of pyridine and benzonitrile revealed by infrared spectroscopy and theory.

The cationic fragmentation products in the dissociative ionization of pyridine and benzonitrile have been studied by infrared action spectroscopy in a cryogenic ion trap instrument at the Free-Electron Lasers for Infrared eXperiments (FELIX) Laboratory. A comparison of the experimental vibrational fingerprints of the dominant cationic fragments with those from quantum chemical calculations revealed a diversity of molecular fragment structures. The loss of HCN/HNC is shown to be the major fragmentation channel for both pyridine and benzonitrile. Using the determined structures of the cationic fragments, potential energy surfaces have been calculated to elucidate the nature of the neutral fragment partner. In the fragmentation chemistry of pyridine, multiple non-cyclic structures are formed, whereas the fragmentation of benzonitrile dominantly leads to the formation of cyclic structures. Among the fragments are linear cyano-(di)acetylene˙+, methylene-cyclopropene˙+ and o- and m-benzyne˙+ structures, the latter possible building blocks in interstellar polycyclic aromatic hydrocarbon (PAH) formation chemistry. Molecular dynamics simulations using density functional based tight binding (MD/DFTB) were performed and used to benchmark and elucidate the different fragmentation pathways based on the experimentally determined structures. The implications of the difference in fragments observed for pyridine and benzonitrile are discussed in an astrochemical context.

Formation of the acenaphthylene cation as a common C2H2-loss fragment in dissociative ionization of the PAH isomers anthracene and phenanthrene.

Polycyclic aromatic hydrocarbons (PAHs) are thought to be a major constituent of astrophysical environments, being the carriers of the ubiquitous aromatic infrared bands (AIBs) observed in the spectra of galactic and extra-galactic sources that are irradiated by ultraviolet (UV) photons. Small (2-cycles) PAHs were unambiguously detected in the TMC-1 dark cloud, showing that PAH growth pathways exist even at low temperatures. The processing of PAHs by UV photons also leads to their fragmentation, which has been recognized in recent years as an alternative route to the generally accepted bottom-up chemical pathways for the formation of complex hydrocarbons in UV-rich interstellar regions. Here we consider the C12H8+ ion that is formed in our experiments from the dissociative ionization of the anthracene and phenanthrene (C14H10) molecules. By employing the sensitive action spectroscopic scheme of infrared pre-dissociation (IRPD) in a cryogenic ion trap instrument coupled to the free-electron lasers at the FELIX Laboratory, we have recorded the broadband and narrow line-width gas-phase IR spectra of the fragment ions (C12H8+) and also the reference spectra of three low energy isomers of C12H8+. By comparing the experimental spectra to those obtained from quantum chemical calculations we have identified the dominant structure of the fragment ion formed in the dissociation process to be the acenaphthylene cation for both isomeric precursors. Ab initio molecular dynamics simulations are presented to elucidate the fragmentation process. This result reinforces the dominant role of species containing a pentagonal ring in the photochemistry of small PAHs.

Hydrogenation of C24 Carbon Clusters: Structural Diversity and Energetic Properties.

This work aims at exploring the potential energy surfaces of C24Hn=0,6,12,18,24 using the genetic algorithm in combination with the density functional based tight binding potential. The structural diversity was analyzed using order parameters, in particular the sum of the numbers of 5- and 6-carbon rings R5/6. The most abundant and lowest energy population was designated as the flake population (isomers of variable shapes, large R5/6 values), characterized by an increasing number of spherical isomers when nH/nC increases. Simultaneously, the fraction of the pretzel population (spherical isomers, smaller R5/6 values) increases. The fraction of the cage population (largest R5/6 values) remains extremely minor while the branched population (smallest R5/6 values) remains the highest energy population for all nH/nC ratios. For all C24Hn=0,6,12,18,24 clusters, the evolution of the carbon ring size distribution with energy clearly shows the correlation between the stability and the number of 6-carbon rings. The average values of the ionization potentials of all populations were found to decrease when nH/nC increases, ranging from 7.9 down to 6.4 eV. This trend was correlated to geometric and electronic factors, in particular to carbon hybridization. These results are of astrophysical interest, especially regarding the role of carbon species in the gas ionization.

Infrared Spectroscopy of Chemically Diverse Carbon Clusters: A Data-Driven Approach.

Carbon clusters exhibit a broad diversity of topologies and shapes, encompassing fullerene-like cages, graphene-like flakes, and more disordered pretzel-like and branched structures. Here, we examine computationally their infrared spectra in relation with these structures from a statistical perspective. Individual spectra for broad samples of isomers were determined by means of the self-consistent charge density functional-based tight-binding method, and an interpolation scheme is designed to reproduce the spectral features by regression on a much smaller subset of the sample. This interpolation proceeds by encoding the structures using appropriate descriptors and selecting them through principal component analysis, Gaussian regression or inverse distance weighting providing the nonlinear weighting functions. Metric learning is employed to reduce the global error on a preselected testing set. The interpolated spectra satisfactorily reproduce the specific spectral features and their dependence on the size and shape, enabling quantitative prediction away from the testing set. Finally, the classification of structures within the four proposed families is critically discussed through a statistical analysis of the sample based on iterative label spreading.

Ultrafast electronic relaxations from the S3 state of pyrene.

The ultrafast relaxation occurring in pyrene upon excitation at 4.68 eV was studied in a supersonic gas-jet fs pump-probe experiment. Mass spectrometry and velocity map imaging of photoelectrons produced by probing via multiphoton ionisation at 800 nm reveal that the initially prepared wave packet exhibits a fast relaxation (<80 fs), followed by a slower one of 200 fs. By comparing the propensity rules of photoionisation observed at one color with ab initio calculations, we tentatively assign these two timescales to a first internal conversion to the dark bB3g state followed by a second one to the long lived aB2u first excited state. Vertical excitation energies determined using ab initio Multi-State Complete Active Space 2nd order Perturbation Theory (MS-CASPT2), as well as oscillator strengths between several electronic states, are reported.

Direct Evidence of the Benzylium and Tropylium Cations as the Two Long-Lived Isomers of C7 H7.

Disentangling the isomeric structure of C7 H7 + is a longstanding experimental issue. We report here the full mid-infrared vibrational spectrum of C7 H7 + tagged with Ne obtained with infrared-predissociation spectroscopy at 10 K. Saturation depletion measurements were used to assign the contribution of benzylium and tropylium isomers and demonstrate that no other isomer is involved. Recorded spectral features compare well with density functional theory calculations. This opens perspectives for a better understanding and control of the formation paths leading to either tropylium or benzylium ions.

Adsorption of PAHs on interstellar ice viewed by classical molecular dynamics.

Polycyclic Aromatic Hydrocarbons (PAHs) are a family of molecules which represent the best candidates to explain the observation of one set of features in the Interstellar Medium (ISM): the Aromatic Interstellar Bands (AIBs). They could also contribute to the Diffuse Interstellar Bands (DIBs). In dense molecular clouds, PAHs may condense onto interstellar grains, contributing to the complex chemistry occurring in their icy mantles, composed essentially of water. In this context, the adsorption of various aromatic molecules, from benzene to ovalene, on different ices - both amorphous and crystalline - is investigated by means of classical molecular dynamics simulations. Initially, a systematic parametrization of the electronic charges on the chosen PAHs in these environments is carried out, and benchmarked with reference to free energies of solvation in liquid water. Then we propose a new, rigorous methodology, transferable to any other PAH or molecular species, to evaluate the charges to be applied to the molecule in the gas phase, at interfaces, or in liquid water. Ultimately, the adsorption energies calculated for the various PAHs are used to derive a function predicting the adsorption energy of any PAH on a given ice surface as a function of the number of C and H atoms it contains. For all PAHs studied, the largest adsorption energies are found on the crystalline hexagonal ice surface (Ih). Binding energy maps constructed for each PAH-ice pair give valuable insight into adsorption site densities and the barriers to surface diffusion. One key result is that the amorphous surface offers a smaller number of adsorption sites compared to that of hexagonal ice. A direct correlation between the location of energetically favourable adsorption sites and the presence of dangling H-bonds is also demonstrated using these maps, showing that PAHs adsorb preferentially on sites offering dangling H-bonds. The present work represents a complete description of PAH-ice interaction in the ground electronic state and at low temperature, providing the binding energies and barrier heights necessary to the ongoing improvement of astrochemical models.

Correction: Adsorption of PAHs on interstellar ice viewed by classical molecular dynamics.

Correction for 'Adsorption of PAHs on interstellar ice viewed by classical molecular dynamics' by Eric Michoulier et al., Phys. Chem. Chem. Phys., 2018, DOI: 10.1039/c8cp00593a.

Atomic hydrogen interactions with gas-phase coronene cations: hydrogenation versus fragmentation.

Sequential hydrogenation of polycyclic aromatic hydrocarbon (PAH) cations drives a gradual transition from a planar to a puckered geometry and from an aromatic to an aliphatic electronic structure. The resulting H-induced weakening of the molecular structure together with the exothermic nature of the consecutive H-attachment processes can lead to substantial molecular fragmentation. We have studied H attachment to gas-phase coronene cations in a radiofrequency ion trap using tandem mass spectrometry. With increasing hydrogenation, C2Hi loss and multifragmentation are identified as main de-excitation channels. To understand the dependence of both channels on H-exposure time, we have simulated the molecular stability and fragmentation channels of hydrogenated PAHs using a molecular dynamics approach employing potential energies determined by a density functional based tight binding method. As the coronene fragmentation patterns depend on the balance between energy deposition by H-attachment and the extent of cooling in between subsequent attachment processes, we investigate several scenarios for the energy distribution of hydrogenated PAHs. Good agreement between experiment and simulation is reached, when realistic energy distributions are considered.

Theoretical determination of adsorption and ionisation energies of polycyclic aromatic hydrocarbons on water ice.

In dense interstellar environments, Polycyclic Aromatic Hydrocarbons (PAHs) are likely to condense onto or integrate into water ice mantles covering dust grains. Understanding the role of ice in the photo-induced processes involving adsorbed PAHs is therefore a key issue in astrochemistry. This requires (i) the knowledge of PAH-ice interactions, i.e. PAH-ice adsorption energies and local structures at the PAH-ice interface, as well as (ii) the understanding of the fate of electrons in the PAH-ice system upon excitation. Regarding (i), in this work, we determined the lowest energy structures of PAH-ice systems for a variety of PAHs ranging from naphthalene to ovalene on three types of ice - crystalline (Ih and Ic) and amorphous (low density) - using an explicit description of the electrons and a finite-sized system. The electronic structure was determined using the Self Consistent Charge Density Functional based Tight Binding (SCC-DFTB) scheme with modified Mulliken charges in order to ensure a good description of the studied systems. Regarding (ii), the influence of the interaction with ice on the Vertical Ionisation Potentials (VIPs) of the series of PAHs was determined using the constrained SCC-DFTB scheme benchmarked against correlated wavefunction results for PAH-(H2O)n (n = 1-6, 13) clusters. The results show a deviation equal, at most, to ∼1.4 eV of the VIPs of PAHs adsorbed on ice with respect to the gas phase values. Our results are discussed in the light of experimental data and previous theoretical studies.

Photochemistry of Fe:H2O Adducts in Argon Matrixes: A Combined Experimental and Theoretical Study in the Mid-IR and UV-Visible Regions.

The photochemistry of Fe:H2O adducts is of interest in fields as diverse as catalysis and astrochemistry. Industrially, iron can be used as a catalyst to convert H2O to H2, whereas in the interstellar medium it may be an important component of dust grains, influencing the chemistry on their icy surfaces. This study consisted of the deposition and spectral characterization of binary systems of atomic iron with H2O in cryogenic argon matrixes. In this way, we were able to obtain information about the interaction of the two species; we observed the formation of adducts of iron monomers and dimers with water molecules in the mid-IR and UV-visible spectral domains. Upon irradiation with a UV radiation source, the iron species were inserted into the water molecules to form HFeOH and HFe2OH, leading in some cases to the formation of FeO possibly accompanied by the production of H2. DFT and correlated multireference wave function calculations confirmed our attributions. This combination of IR and UV-visible spectroscopy with theoretical calculations allowed us to determine, for the first time, the spectral characteristics of iron adducts and their photoproducts in the UV-visible and in the OH stretching region of the mid-IR domain.

Investigating the importance of edge-structure in the loss of H/H2 of PAH cations: the case of dibenzopyrene isomers.

We present a detailed study of the main dehydrogenation processes of two dibenzopyrene cation (C24H14+) isomers, namely dibenzo(a,e)pyrene (AE+) and dibenzo(a,l)pyrene (AL+). First, action spectroscopy under VUV photons was performed using synchrotron radiation in the 8-20 eV range. We observed lower dissociation thresholds for the non-planar molecule (AL+) than for the planar one (AE+) for the main dissociation pathways: H and 2H/H2 loss. In order to rationalize the experimental results, dissociation paths were investigated by means of density functional theory calculations. In the case of H loss, which is the dominant channel at the lowest energies, the observed difference between the two isomers can be explained by the presence in AL+ of two C-H bonds with considerably lower adiabatic dissociation energies. In both isomers the 2H/H2 loss channels are observed only at about 1 eV higher than H loss. We suggest that this is due to the propensity of bay H atoms to easily form H2. In addition, in the case of AL+, we cannot exclude a competition between 2H and H2 channels. In particular, the formation of a stable dissociation product with a five-membered ring could account for the low energy sequential loss of 2 hydrogens. This work shows the potential role of non-compact PAHs containing bay regions in the production of H2 in space.

Identification of the fragment of the 1-methylpyrene cation by mid-IR spectroscopy.

The fragment of the 1-methylpyrene cation, C17H11+, is expected to exist in two isomeric forms, 1-pyrenemethylium PyrCH2+ and the tropylium containing species PyrC7+. We measured the infrared (IR) action spectrum of cold C17H11+ tagged with Ne using a cryogenic ion trap instrument coupled to the FELIX laser. Comparison of the experimental data with density functional theory calculations allows us to identify the PyrCH2+ isomer in our experiments. The IR Multi-Photon Dissociation spectrum was also recorded following the C2H2 loss channel. Its analysis suggests combined effects of anharmonicity and isomerisation while heating the trapped ions, as shown by molecular dynamics simulations.

Theoretical investigation of the solid-liquid phase transition in protonated water clusters.

Protonated water clusters have received a lot of attention as they offer tools to bridge the gap between molecular and bulk scales of water. However, their properties are still not fully understood and deserve further theoretical and experimental investigations. In this work, we simulate the caloric curves of protonated water clusters (H2O)nH+ (n = 20-23). These curves, which have recently been measured experimentally, are characteristic of the phase changes occurring in the aggregates with respect to temperature. The present simulations are achieved by combining parallel-tempering molecular dynamics and the self-consistent-charge density-functional based tight-binding approach and are focused on a restricted size range around (H2O)21H+ which presents singular properties. The shape of the experimental caloric curves and their size dependence are satisfactorily reproduced by the simulations which allows us to further provide a description of the phase transition in terms of structural modifications, dynamics of water molecules and proton mobility. Similar to the experiments, we observe that (H2O)21H+ exhibits a sharper phase transition than the neighbouring size clusters, which can be traced back to both structural and dynamic peculiarities.

Structural Characterization of Sulfur-Containing Water Clusters Using a Density-Functional Based Tight-Binding Approach.

A global optimization search of low-energy isomers is carried out to investigate the structural and stability properties of sulfur-containing water clusters, including both (H2O)nSO42- and (H2O)nH2SO4 aggregates. The systematic optimization algorithm involves a combination of parallel-tempering molecular dynamics and periodic gradient-driven quenches with energy and energy-gradient calculations performed using the Self-Consistent-Charge Density-Functional based Tight-Binding (SCC-DFTB) scheme. Comparisons with new MP2 and DFT calculations on the smallest systems and previous ab initio investigations of the literature show that the SCC-DFTB approach provides a fairly accurate description of both neutral and ionic species, comparable to that of DFT. Structural and stability features of larger sulfur-containing clusters, with up to 20 water molecules, are also determined using the SCC-DFTB scheme. The interest of this work is 2-fold: (i) the benchmark on small species demonstrates the ability of SCC-DFTB to describe complex potential energy landscapes involving hydrogen-bonds and proton transfers; (ii) it opens the way to the study of large clusters that can hardly be performed within ab initio approaches.

Phase changes of the water hexamer and octamer in the gas phase and adsorbed on polycyclic aromatic hydrocarbons.

We investigate thermodynamic properties of small water clusters adsorbed on polycyclic aromatic hydrocarbons (PAHs), which are relevant systems in the context of astrophysical and atmospheric chemistry. We present heat capacity curves computed from parallel-tempering molecular dynamics and Monte Carlo simulations that were performed using the self-consistent-charge density-functional based tight-binding method. These curves are characteristic of the phase changes occurring in the aggregates and provide useful information on the evolution of the interaction between the water molecules and the PAHs as a function of temperature. After benchmarking our approach on the water hexamer and octamer in the gas phase, we present some results for these same clusters adsorbed on coronene and circumcoronene. When compared to the curves obtained for the isolated water clusters, the phase change temperature significantly decreases for the (H2O)8-PAH clusters whereas it depends on the nature of the PAH in the case of the hexamer. We analyse these differences as connected to the relative energies of the optimized characteristic isomers and to their dynamical behavior. We also evidence the population changes of the various cluster isomers as a function of temperature.

Electronic Spectroscopy of [FePAH](+) Complexes in the Region of the Diffuse Interstellar Bands: Multireference Wave Function Studies on [FeC6H6](+).

The low-energy states and electronic spectrum in the near-infrared-visible region of [FeC6H6](+) are studied by theoretical approaches. An exhaustive exploration of the potential energy surface of [FeC6H6](+) is performed using the density functional theory method. The ground state is found to be a (4)A1 state. The structures of the lowest energy states ((4)A2 and (4)A1) are used to perform multireference wave function calculations by means of the multistate complete active space with perturbation at the second order method. Contrary to the density functional theory results ((4)A1 ground state), multireference perturbative calculations show that the (4)A2 state is the ground state. The vertical electronic spectrum is computed and compared with the astronomical diffuse interstellar bands, a set of near-infrared-visible bands detected on the extinction curve in our and other galaxies. Many transitions are found in this domain, corresponding to d → d, d → 4s, or d → π* excitations, but few are allowed and, if they are, their oscillation strengths are small. Even though some band positions could match some of the observed bands, the relative intensities do not fit, making the contribution of the [Fe-C6H6](+) complexes to the diffuse interstellar bands questionable. This work, however, lays the foundation for the studies of polycyclic aromatic hydrocarbons (PAHs) complexed to Fe cations that are more likely to possess d → π* and π → π* transitions in the diffuse interstellar bands domain. PAH ligands indeed possess a larger number of π and π* orbitals, respectively, higher and lower in energy than those of C6H6, which are expected to lead to lower energy d → π* and π → π* transitions in [FePAH](+) than in [FeC6H6](+) complexes.

Water clusters in an argon matrix: infrared spectra from molecular dynamics simulations with a self-consistent charge density functional-based tight binding/force-field potential.

The present theoretical study aims at investigating the effects of an argon matrix on the structures, energetics, dynamics, and infrared (IR) spectra of small water clusters (H2O)n (n = 1-6). The potential energy surface is obtained from a hybrid self-consistent charge density functional-based tight binding/force-field approach (SCC-DFTB/FF) in which the water clusters are treated at the SCC-DFTB level and the matrix is modeled at the FF level by a cluster consisting of ∼340 Ar atoms with a face centered cubic (fcc) structure, namely (H2O)n/Ar. With respect to a pure FF scheme, this allows a quantum description of the molecular system embedded in the matrix, along with all-atom geometry optimization and molecular dynamics (MD) simulations of the (H2O)n/Ar system. Finite-temperature IR spectra are derived from the MD simulations. The SCC-DFTB/FF scheme is first benchmarked on (H2O)Arn clusters against correlated wave function results and DFT calculations performed in the present work, and against FF data available in the literature. Regarding (H2O)n/Ar systems, the geometries of the water clusters are found to adapt to the fcc environment, possibly leading to intermolecular distortion and matrix perturbation. Several energetical quantities are estimated to characterize the water clusters in the matrix. In the particular case of the water hexamer, substitution and insertion energies for the prism, bag, and cage are found to be lower than that for the 6-member ring isomer. Finite-temperature MD simulations show that the water monomer has a quasifree rotation motion at 13 K, in agreement with experimental data. In the case of the water dimer, the only large-amplitude motion is a distortion-rotation intermolecular motion, whereas only vibration motions around the nuclei equilibrium positions are observed for clusters with larger sizes. Regarding the IR spectra, we find that the matrix environment leads to redshifts of the stretching modes and almost no shift of the bending modes. This is in agreement with experimental data. Furthermore, in the case of the water monomer and dimer, the magnitudes of the computed shifts are in fair agreement with the experimental values. The complex case of the water hexamer, which presents several low-energy isomers, is discussed.

Cationic Methylene-Pyrene Isomers and Isomerization Pathways: Finite Temperature Theoretical Studies.

This paper provides spectral characterizations of the two isomers of the 1-methylenepyrene cation, namely, the 1-pyrenemethylium and a pyrene-like isomer owing a tropylium cycle. Both are possible photodissociation products of the 1-methylpyrene cation and were proposed as potential contributors to the diffuse interstellar bands. In that respect, vibrational and electronic spectra are computed for the optimized structures at the density functional theory (DFT) and time-dependent (TD-)DFT levels. Finite temperature effects on these spectra are estimated from molecular dynamics simulations within the density functional-based tight-binding (DFTB) and TD-DFTB frameworks, these methods being first benchmarked against DFT and TD-DFT calculations. The computed spectra allow discrimination of the two isomers. When the temperature increases, bands are observed to redshift and merge. The isomerization mechanism is investigated with the metadynamics technique, a biased dynamics scheme allowing to probe reaction mechanisms with high energy barriers by investigating the free energy surface at various temperatures. Four pathways with similar barrier heights (3.5-4 eV) are found, showing that the interconversion process would only occur in interstellar clouds under photoactivation. The present study opens the way to simulations on larger methyl- and methylenePAHs of astrophysical interest and their experimental investigation.