A far-ultraviolet-driven photoevaporation flow observed in a protoplanetary disk

Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photodissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, which affects planet formation within the disks. We report James Webb Space Telescope and Atacama Large Millimeter Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modeling their kinematics and excitation allowed us to constrain the physical conditions within the gas. We quantified the mass-loss rate induced by the FUV irradiation and found that it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk.

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.

Threshold collision induced dissociation of protonated water clusters

We report threshold collision induced dissociation experiments on protonated water clusters thermalized at low temperature for sizes n = 19-23. Fragmentation cross sections are recorded as a function of the collision energy and analyzed with a statistical model. This model allows us to account for dissociation cascades and provides values for the dissociation energies of each cluster. These values, averaging around 0.47 eV, are in good agreement with theoretical predictions at various levels of theory. Furthermore, the dissociation energies show a trend for the n = 21 magic and n = 22 anti-magic numbers relative to their neighbours, which is also in agreement with theory. These results provide further evidence to resolve the disagreement between previously published experimental values. A careful quantitative treatment of cascade dissociation in this model introduces interdependence between the dissociation energies of neighboring sizes, which reduces the number of free fitting parameters and improves both reliability and uncertainties on absolute dissociation energies deduced from experiments.

Formation of the methyl cation by photochemistry in a protoplanetary disk

Forty years ago, it was proposed that gas-phase organic chemistry in the interstellar medium can be initiated by the methyl cation CH3+ (refs. 1-3), but so far it has not been observed outside the Solar System4,5. Alternative routes involving processes on grain surfaces have been invoked6,7. Here we report James Webb Space Telescope observations of CH3+ in a protoplanetary disk in the Orion star-forming region. We find that gas-phase organic chemistry is activated by ultraviolet irradiation.

Atmospheric molecular blobs shape up circumstellar envelopes of AGB stars

During their thermally pulsing phase, asymptotic giant branch (AGB) stars eject material that forms extended dusty envelopes1. Visible polarimetric imaging found clumpy dust clouds within two stellar radii of several oxygen-rich stars2-6. Inhomogeneous molecular gas has also been observed in multiple emission lines within several stellar radii of different oxygen-rich stars, including W Hya and Mira7-10. At the stellar surface level, infrared images have shown intricate structures around the carbon semiregular variable R Scl and in the S-type star π1 Gru11,12. Infrared images have also shown clumpy dust structures within a few stellar radii of the prototypical carbon AGB star IRC+10°216 (refs. 13,14), and studies of molecular gas distribution beyond the dust formation zone have also shown complex circumstellar structures15. Because of the lack of sufficient spatial resolution, however, the distribution of molecular gas in the stellar atmosphere and the dust formation zone of AGB carbon stars is not known, nor is how it is subsequently expelled. Here we report observations with a resolution of one stellar radius of the recently formed dust and molecular gas in the atmosphere of IRC+10°216. Lines of HCN, SiS and SiC2 appear at different radii and in different clumps, which we interpret as large convective cells in the photosphere, as seen in Betelgeuse16. The convective cells coalesce with pulsation, causing anisotropies that, together with companions17,18, shape its circumstellar envelope.

Electronic effects in the dissociative ionisation of pyrene clusters

We present a combined experimental and theoretical study on the dissociative ionisation of clusters of pyrene. We measured the experimental appearance energies in the photon energy range 7.2-12.0 eV of the fragments formed from neutral monomer loss for clusters up to the hexamer. The results obtained show a deviation from statistical dissociation. From electronic structure calculations, we suggest that the role of excited states must be considered in the interpretation of experimental results, even in these relatively large systems. Non-statistical effects in the dissociative ionization process of polycyclic aromatic hydrocarbon (PAH) clusters may have an impact on the assessment of mechanisms determining the stability of these clusters in astrophysical environments.

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.

Isomer Differentiation of Trapped C16H10+ Using Low-Energy Collisions and Visible/VUV Photons

Polycyclic aromatic hydrocarbons are major species in astrophysical environments, and this motivates their study in samples of astrochemical interest such as meteorites and laboratory analogues of stardust. Molecular analyses of carbonaceous matter in these samples show a dominant peak at m/z = 202.078 corresponding to C₁₆H₁₀. Obtaining information on the associated isomeric structures is a challenge for the molecular analysis of samples available in very small quantities (mg or less). Here we show that coupling laser desorption ionization mass spectrometry with ion trapping opens up the possibility of unraveling isomers by activating ion fragmentation via collisions or photon absorption. We report the best criteria for differentiating isomers with comparable dissociation energies, namely pyrene, fluoranthene, and 9-ethynylphenanthrene, on the basis of the parent dissociation curve and the ratio of dehydrogenation channels. Photoabsorption schemes (multiple photon absorption in the visible range and single photon absorption at 10.5 eV) are more effective in differentiating these isomers than activation by low energy collisions. The impact of the activation scheme on the fragmentation kinetics and dehydrogenation pathways is discussed. By analyzing the 10.5 eV photodissociation measurements with a simple kinetic model, we were able to derive a branching ratio for the H and 2H/H₂ loss channels of the parent ions. The results suggest a role in the formation of H₂ for bay hydrogens that are present in both fluoranthene and 9-ethynylphenanthrene. In addition, we suggest for the latter the presence of a highly competitive 2H loss channel, possibly associated with the formation of a pentagonal ring.

Water Attachment onto Size-Selected Cationic Pyrene Clusters

We report measurements of the attachment rates of water molecules onto mass-selected cationic pyrene clusters for size from n = 4 to 13 pyrene units and for different collision energies. Comparison of the attachment rates with the collision rates measured in collision-induced dissociation experiments provides access to the values of the sticking coefficient. The strong dependence of the attachment rates on size and collision energy is rationalized through a model in which we use a Langevin-type collision rate and adjust on experimental data the statistical dissociation rate of the water molecule from the cluster after attachment. This allows us to extrapolate our results to the conditions of isolation and long time scales encountered in astrophysical environments.

Anharmonic Infrared Spectra of Thermally Excited Pyrene (C16H10): A Combined View of DFT-Based GVPT2 with AnharmonicCaOs, and Approximate DFT Molecular dynamics with DemonNano

The study of the Aromatic Infrared Bands (AIBs) in astronomical environments has opened interesting spectroscopic questions on the effect of anharmonicity on the infrared (IR) spectrum of hot polycyclic aromatic hydrocarbons (PAHs) and related species in isolated conditions. The forthcoming James Webb Space Telescope will unveil unprecedented spatial and spectral details in the AIB spectrum; significant advancement is thus necessary now to model the infrared emission of PAHs, their presumed carriers, with enough detail to exploit the information content of the AIBs. This requires including anharmonicity in such models, and to do so systematically for all species included, requiring a difficult compromise between accuracy and efficiency. We performed a benchmark study to compare the performances of two methods in calculating anharmonic spectra, comparing them to available experimental data. One is a full quantum method, AnharmoniCaOs, relying on an ab initio potential, and the other relies on Molecular Dynamics simulations using a Density Functional based Tight Binding potential. The first one is found to be very accurate and detailed, but it becomes computationally very expensive for increasing temperature; the second is faster and can be used for arbitrarily high temperatures, but is less accurate. Still, its results can be used to model the evolution with temperature of isolated bands. We propose a new recipe to model anharmonic AIB emission using minimal assumptions on the general behaviour of band positions and widths with temperature, which can be defined by a small number of empirical parameters. Modelling accuracy will depend critically on these empirical parameters, allowing for an incremental improvement in model results, as better estimates become gradually available.

The impact and recovery of asteroid 2018 LA

The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to triangulate the location of the fireball's main flare over the Central Kalahari Game Reserve. 23 meteorites were recovered. A consortium study of eight of these classifies Motopi Pan as a HED polymict breccia derived from howardite, cumulate and basaltic eucrite, and diogenite lithologies. Before impact, 2018 LA was a solid rock of ~156 cm diameter with high bulk density ~2.85 g/cm³, a relatively low albedo pv ~ 0.25, no significant opposition effect on the asteroid brightness, and an impact kinetic energy of ~0.2 kt. The orbit of 2018 LA is consistent with an origin at Vesta (or its Vestoids) and delivery into an Earth-impacting orbit via the v₆ resonance. The impact that ejected 2018 LA in an orbit towards Earth occurred 22.8 ± 3.8 Ma ago. Zircons record a concordant U-Pb age of 4563 ± 11 Ma and a consistent ²⁰⁷Pb/²⁰⁶Pb age of 4563 ± 6 Ma. A much younger Pb-Pb phosphate resetting age of 4234 ± 41 Ma was found. From this impact chronology, we discuss what is the possible source crater of Motopi Pan and the age of Vesta's Veneneia impact basin.

Impact of metals on (star)dust chemistry: a laboratory astrophysics approach

Laboratory experiments are essential in exploring the mechanisms involved in stardust formation. One key question is how a metal is incorporated into dust for an environment rich in elements involved in stardust formation (C, H, O, Si). To address experimentally this question we have used a radiofrequency cold plasma reactor in which cyclic organosilicon dust formation is observed. Metallic (silver) atoms were injected in the plasma during the dust nucleation phase to study their incorporation in the dust. The experiments show formation of silver nanoparticles (~15 nm) under conditions in which organosilicon dust of size 200 nm or less is grown. The presence of AgSiO bonds, revealed by infrared spectroscopy, suggests the presence of junctions between the metallic nanoparticles and the organosilicon dust. Even after annealing we could not conclude on the formation of silver silicates, emphasizing that most of silver is included in the metallic nanoparticles. The molecular analysis performed by laser mass spectrometry exhibits a complex chemistry leading to a variety of molecules including large hydrocarbons and organometallic species. In order to gain insights into the involved chemical molecular pathways, the reactivity of silver atoms/ions with acetylene was studied in a laser vaporization source. Key organometallic species, Ag n C2H m (n=1-3; m=0-2), were identified and their structures and energetic data computed using density functional theory. This allows us to propose that molecular Ag-C seeds promote the formation of Ag clusters but also catalyze hydrocarbon growth. Throughout the article, we show how the developed methodology can be used to characterize the incorporation of metal atoms both in the molecular and dust phases. The presence of silver species in the plasma was motivated by objectives finding their application in other research fields than astrochemistry. Still, the reported methodology is a demonstration laying down the ground for future studies on metals of astrophysical interest such as iron.

Experimental and theoretical study of photo-dissociation spectroscopy of pyrene dimer radical cations stored in a compact electrostatic ion storage ring

In this paper, we present an experimental and theoretical study of the photo-dissociation of free-flying dimer radical cations of pyrene (C16H10)2+. Experimentally, the dimers were produced in the plasma of an electron cyclotron resonance ion source and stored in an electrostatic ion storage ring, the Mini-Ring for times up to 10 ms and the photo-dissociation spectrum was recorded in the 400 to 2000 nm range. Two broad absorption bands were observed at 550 (2.25 eV) and 1560 nm (0.79 eV), respectively. Theoretical simulations of the absorption spectrum as a function of the temperature were performed using the Density Functional based Tight Binding approach within the Extended Configuration Interaction scheme (DFTB-EXCI) to determine the electronic structure. The simulation involved all excited electronic states correlated asymptotically with the five lowest excited states D1-D5 of the monomer cation and a Monte Carlo exploration of the electronic ground state potential energy surface. The simulations exhibit three major bands at 1.0, 2.1 and 2.8 eV respectively. They allow assigning the experimental band at 1560 nm to absorption by the charge resonance (CR) excited state correlated with the ground state of the monomer D0. The band at 550 nm is tentatively attributed to dimer states correlated with excited states D2-D4, in the monomer cation. Simulations also show that the CR band broadens and shifts towards longer wavelength with increasing temperature. It results from the dependence on the geometry of the energy gap between the ground state and the lowest excited state. The comparison of the experimental spectrum with theoretical spectra at various temperatures allows us to estimate the temperature of the stored (C16H10)2+ in the 300-400 K range, which is also in line with the expected temperatures of the ions deduced from the analysis of the natural decay curve.

Molecular content of nascent soot: Family characterization using two-step laser desorption laser ionization mass spectrometry

Molecules constituting nascent soot particles have been analyzed by two-step laser desorption laser ionization mass spectrometry. Three samples have been collected from a slightly sooting ethylene/air premixed flame with the aim to investigate soot composition in the transition from nucleated to just-grown soot particles. Sampling locations have been selected based on the evolution of the particle size distribution along the flame axis. The mass spectrometric results point to a strong evolution of the molecular composition. Just-nucleated soot is rich in polycyclic aromatic hydrocarbons (PAHs) dominated by medium sizes from 18 to 40 carbon atoms but containing sizes as large as 90 carbon atoms. Most abundant PAHs are in the form of peri-condensed structures. The presence of a large fraction of odd numbered carbon species shows that pentagonal cycles are a common feature of the detected population. Increasing the distance from the burner outlet, i.e., the particle residence time in flame, leads to an evolution of the chemical composition of this population with a major contribution of carbon clusters including also fullerenes up to about 160 carbon atoms. Our data support a scenario in which large PAHs containing pentagonal rings evolve very efficiently upon thermal processing by a series of dehydrogenation and isomerization processes to form fullerenes. This chemistry happens in the early steps of soot growth showing that carbonization is already active at this stage. © 2020 The Authors. Published by Elsevier Inc. on behalf of The Combustion Institute. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Astrochemical relevance of VUV ionization of large PAH cations. Astron Astrophys 2020;641:A98. [PMID: 33154599 PMCID: PMC7116310 DOI: 10.1051/0004-6361/202038139]

CONTEXT

As a part of interstellar dust, polycyclic aromatic hydrocarbons (PAHs) are processed by the interaction with vacuum ultra-violet (VUV) photons that are emitted by hot young stars. This interaction leads to the emission of the well-known aromatic infrared bands but also of electrons, which can significantly contribute to the heating of the interstellar gas.

AIMS

Our aim is to investigate the impact of molecular size on the photoionization properties of cationic PAHs.

METHODS

Trapped PAH cations of sizes between 30 and 48 carbon atoms were submitted to VUV photons in the range of 9 to 20 eV from the DESIRS beamline at the synchrotron SOLEIL. All resulting photoproducts including dications and fragment cations were mass-analyzed and recorded as a function of photon energy.

RESULTS

Photoionization is found to be predominant over dissociation at all energies, which differs from an earlier study on smaller PAHs. The photoionization branching ratio reaches 0.98 at 20 eV for the largest studied PAH. The photoionization threshold is observed to be between 9.1 and 10.2 eV, in agreement with the evolution of the ionization potential with size. Ionization cross sections were indirectly obtained and photoionization yields extracted from their ratio with theoretical photoabsorption cross sections, which were calculated using time-dependent density functional theory. An analytical function was derived to calculate this yield for a given molecular size.

CONCLUSIONS

Large PAH cations could be efficiently ionized in H i regions and provide a contribution to the heating of the gas by photoelectric effect. Also, at the border of or in H ii regions, PAHs could be exposed to photons of energy higher than 13.6 eV. Our work provides recipes to be used in astronomical models to quantify these points.

Threshold collision induced dissociation of pyrene cluster cations

We report threshold collision induced dissociation experiments on cationic pyrene clusters, for sizes n = 2-6. Fragmentation cross sections are recorded as a function of the collision energy and analyzed with a statistical model. This model can account for the dissociation cascades and provides values for the dissociation energies. These values, of the order of 0.7 eV-1 eV, are in excellent agreement with those previously derived from thermal evaporation. They confirm the charge resonance stability enhancement predicted by theoretical calculations. In addition, remarkable agreement is obtained with theoretical predictions for the two smaller sizes n = 2 and 3. For the larger sizes, the agreement remains good, although the theoretical values obtained for the most stable structures are systematically higher by 0.2 eV. This offset could be attributed to approximations in the calculations. Still, there is an indication in the results of an incomplete description of the role of isomerization and/or direct dissociation upon collisions. Finally, by-product clusters containing dehydrogenated species are found to dissociate at energies comparable to the non-dehydrogenated ones, which shows no evidence for covalent bonds within the clusters.

The Chemistry of Cosmic Dust Analogues from C, C2, and C2H2 in C-Rich Circumstellar Envelopes

Interstellar carbonaceous dust is mainly formed in the innermost regions of circumstellar envelopes around carbon-rich asymptotic giant branch (AGB) stars. In these highly chemically stratified regions, atomic and diatomic carbon, along with acetylene are the most abundant species after H2 and CO. In a previous study, we addressed the chemistry of carbon (C and C2) with H2 showing that acetylene and aliphatic species form efficiently in the dust formation region of carbon-rich AGBs whereas aromatics do not. Still, acetylene is known to be a key ingredient in the formation of linear polyacetylenic chains, benzene and polycyclic aromatic hydrocarbons (PAHs), as shown by previous experiments. However, these experiments have not considered the chemistry of carbon (C and C2) with C2H2. In this work, by employing a sufficient amount of acetylene, we investigate its gas-phase interaction with atomic and diatomic carbon. We show that the chemistry involved produces linear polyacetylenic chains, benzene and other PAHs, which are observed with high abundances in the early evolutionary phase of planetary nebulae. More importantly, we have found a non-negligible amount of pure and hydrogenated carbon clusters as well as aromatics with aliphatic substitutions, both being a direct consequence of the addition of atomic carbon. The incorporation of alkyl substituents into aromatics can be rationalized by a mechanism involving hydrogen abstraction followed by methyl addition. All the species detected in gas phase are incorporated into the nanometric sized dust analogues, which consist of a complex mixture of sp, sp2 and sp3 hydrocarbons with amorphous morphology.

Prevalence of non-aromatic carbonaceous molecules in the inner regions of circumstellar envelopes

Evolved stars are a foundry of chemical complexity, gas and dust that provides the building blocks of planets and life, and dust nucleation first occurs in their photosphere. Despite their importance, the circumstellar regions enveloping these stars remain hidden to many observations, thus dust formation processes are still poorly understood. Laboratory astrophysics provides complementary routes to unveil these chemical processes, but most experiments rely on combustion or plasma decomposition of molecular precursors under physical conditions far removed from those in space. We have built an ultra-high vacuum machine combining atomic gas aggregation with advanced in-situ characterization techniques to reproduce and characterize the bottom-up dust formation process. We show that carbonaceous dust analogues formed from low-pressure gas-phase condensation of C atoms in a hydrogen atmosphere, in a C/H2 ratio similar to that reported for evolved stars, leads to the formation of amorphous C nanograins and aliphatic C-clusters. Aromatic species or fullerenes do not form effectively under these conditions, raising implications for the revision of the chemical mechanisms taking place in circumstellar envelopes.

Learning mid-IR emission spectra of polycyclic aromatic hydrocarbon populations from observations

CONTEXT

The James Webb Space Telescope (JWST) will deliver an unprecedented quantity of high-quality spectral data over the 0.6-28 μm range. It will combine sensitivity, spectral resolution, and spatial resolution. Specific tools are required to provide efficient scientific analysis of such large data sets.

AIMS

Our aim is to illustrate the potential of unsupervised learning methods to get insights into chemical variations in the populations that carry the aromatic infrared bands (AIBs), more specifically polycyclic aromatic hydrocarbon (PAH) species and carbonaceous very small grains (VSGs).

METHODS

We present a method based on linear fitting and blind signal separation for extracting representative spectra for a spectral data set. The method is fast and robust, which ensures its applicability to JWST spectral cubes. We tested this method on a sample of ISO-SWS data, which resemble most closely the JWST spectra in terms of spectral resolution and coverage.

RESULTS

Four representative spectra were extracted. Their main characteristics appear consistent with previous studies with populations dominated by cationic PAHs, neutral PAHs, evaporating VSGs, and large ionized PAHs, known as the PAH ^x population. In addition, the 3 μm range, which is considered here for the first time in a blind signal separation (BSS) method, reveals the presence of aliphatics connected to neutral PAHs. Each representative spectrum is found to carry second-order spectral signatures (e.g., small bands), which are connected with the underlying chemical diversity of populations. However, the precise attribution of theses signatures remains limited by the combined small size and heterogeneity of the sample of astronomical spectra available in this study.

CONCLUSIONS

The upcoming JWST data will allow us to overcome this limitation. The large data sets of hyperspectral images provided by JWST analysed with the proposed method, which is fast and robust, will open promising perspectives for our understanding of the chemical evolution of the AIB carriers.

Thermal evaporation of pyrene clusters

This work presents a study of the thermal evaporation and stability of pyrene (C16H10)n clusters. Thermal evaporation rates of positively charged mass-selected clusters are measured for sizes in the range n = 3-40 pyrene units. The experimental setup consists of a gas aggregation source, a thermalization chamber, and a time of flight mass spectrometer. A microcanonical Phase Space Theory (PST) simulation is used to determine the dissociation energies of pyrene clusters by fitting the experimental breakdown curves. Calculations using the Density Functional based Tight Binding combined with a Configuration Interaction (CI-DFTB) model and a hierarchical optimization scheme are also performed in the range n = 2-7 to determine the harmonic frequencies and a theoretical estimation of the dissociation energies. The frequencies are used in the calculations of the density of states needed in the PST simulations, assuming an extrapolation scheme for clusters larger than 7 units. Using the PST model with a minimal set of adjustable parameters, we obtain good fits of the experimental breakdown curves over the full studied size range. The approximations inherent to the PST simulation and the influence of the used parameters are carefully estimated. The derived dissociation energies show significant variations over the studied size range. Compared with neutral clusters, significantly higher values of the dissociation energies are obtained for the smaller sizes and attributed to charge resonance in line with CI-DFTB calculations.

Multi-scale investigation in the frequency domain of Ar/HMDSO dusty plasma with pulsed injection of HMDSO

A combination of time-resolved optical emission spectroscopy measurements and collisional-radiative modeling is used to investigate the phenomena occurring over multiple time scales in the frequency domain of a low-pressure, axially-asymmetric capacitively-coupled RF argon plasma with pulsed injection of hexamethyldisiloxane (HMDSO, Si₂O(CH₃)₆). The collisional-radiative model developed here considers the population of argon 1s and all ten 2p levels (in Paschen's notation). The presence of HMDSO in the plasma is accounted for in the model by quenching of the argon 1s states by species generated by plasma processing of HMDSO, including HMDSO-15 (Si₂O(CH₃)₅), acetylene (C₂H₂) and methane (CH₄). Detailed analysis of the relative populations of Ar 2p states reveals cyclic evolutions of the electron temperature, electron density and quenching frequency that are shown to be linked to the kinetics of dust formation in Ar/HMDSO plasmas. Penning ionization of HMDSO and its fragments is found to be an important source of electrons for the plasma maintenance. It is at the origin of the cyclic formation/disappearance of the dust cloud, without attenuation of the phenomenon, as long as the pulsed injection of HMDSO is sustained. The multi-scale approach used in this study further reveals the straightforward relation of the frequency of HMDSO pulsed injection, in particular the HMDSO duty cycle, with the frequency of dust formation/disappearance cycle.

Experimental Approach to the Study of Anharmonicity in the Infrared Spectrum of Pyrene from 14 to 723 K

Quantifying the effect of anharmonicity on the infrared spectrum of large molecules such as polycyclic aromatic hydrocarbons (PAHs) at high temperatures is the focus of a number of theoretical and experimental studies, many of them motivated by astrophysical applications. We recorded the IR spectrum of pyrene C16H10 microcrystals embedded in KBr pellets over a wide range of temperatures (14-723 K) and studied the evolution of band positions, widths, and integrated intensities with temperature. We identified jumps for some of the spectral characteristics of some bands in the 423-473 K range. These were attributed to a change of phase from crystal to molten in condensed pyrene, which appears to affect more strongly bands involving large CH motions. Empirical anharmonicity factors that quantify the linear evolution of band positions and widths with temperature for values larger than ∼150-250 K, depending on the band, were retrieved from both phases and averaged to provide recommended values for these anharmonicity factors. The derived values were found to be consistent with available gas phase data. We conclude about the relevance of the methodology to produce data that can be compared with calculated anharmonic IR spectra and provide input for models that simulate the IR emission of astro-PAHs.

Anharmonic vibrational spectroscopy of polycyclic aromatic hydrocarbons (PAHs)

While powerful techniques exist to accurately account for anharmonicity in vibrational molecular spectroscopy, they are computationally very expensive and cannot be routinely employed for large species and/or at non-zero vibrational temperatures. Motivated by the study of Polycyclic Aromatic Hydrocarbon (PAH) emission in space, we developed a new code, which takes into account all modes and can describe all infrared transitions including bands becoming active due to resonances as well as overtone, combination, and difference bands. In this article, we describe the methodology that was implemented and discuss how the main difficulties were overcome, so as to keep the problem tractable. Benchmarking with high-level calculations was performed on a small molecule. We carried out specific convergence tests on two prototypical PAHs, pyrene (C₁₆H₁₀) and coronene (C₂₄H₁₂), aiming at optimising tunable parameters to achieve both acceptable accuracy and computational costs for this class of molecules. We then report the results obtained at 0 K for pyrene and coronene, comparing the calculated spectra with available experimental data. The theoretical band positions were found to be significantly improved compared to harmonic density functional theory calculations. The band intensities are in reasonable agreement with experiments, the main limitation being the accuracy of the underlying calculations of the quartic force field. This is a first step toward calculating moderately high-temperature spectra of PAHs and other similarly rigid molecules using Monte Carlo sampling.

Direct Evidence of the Benzylium and Tropylium Cations as the Two Long-Lived Isomers of C7 H7. Chemphyschem 2018;19:3182-3185. [PMID: 30238585 PMCID: PMC6420061 DOI: 10.1002/cphc.201800744]

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.

Structure of photodissociation fronts in star-forming regions revealed by observations of high-J CO emission lines with Herschel

CONTEXT

In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines.

AIMS

We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces.

METHODS

We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code.

RESULTS

A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to J up =23 in the Orion Bar (J up =19 in NGC 7023), can only originate from small structures of typical thickness of a few 10-3 pc and at high thermal pressures (Pth ~ 108 K cm-3).

CONCLUSIONS

Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.

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 (C₂₄H₁₄⁺) 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/H₂ 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/H₂ 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 H₂. In addition, in the case of AL⁺, we cannot exclude a competition between 2H and H₂ 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 H₂ in space.

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

The fragment of the 1-methylpyrene cation, C 17 H 11 + , is expected to exist in two isomeric forms, 1-pyrenemethylium PyrCH 2 + and the tropylium containing species PyrC 7 + . We measured the infrared (IR) action spectrum of cold C 17 H 11 + 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 PyrCH 2 + 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.

Unimolecular reaction energies for polycyclic aromatic hydrocarbon ions

Imaging photoelectron photoion coincidence spectroscopy was employed to explore the unimolecular dissociation of the ionized polycyclic aromatic hydrocarbons (PAHs) acenaphthylene, fluorene, cyclopenta[d,e,f]phenanthrene, pyrene, perylene, fluoranthene, dibenzo[a,e]pyrene, dibenzo[a,l]pyrene, coronene and corannulene. The primary reaction is always hydrogen atom loss, with the smaller species also exhibiting loss of C2H2 to varying extents. Combined with previous work on smaller PAH ions, trends in the reaction energies (E0) for loss of H from sp2-C and sp3-C centres, along with hydrocarbon molecule loss were found as a function of the number of carbon atoms in the ionized PAHs ranging in size from naphthalene to coronene. In the case of molecules which possessed at least one sp3-C centre, the activation energy for the loss of an H atom from this site was 2.34 eV, with the exception of cyclopenta[d,e,f]phenanthrene (CPP) ions, for which the E0 was 3.44 ± 0.86 eV due to steric constraints. The hydrogen loss from PAH cations and from their H-loss fragments exhibits two trends, depending on the number of unpaired electrons. For the loss of the first hydrogen atom, the energy is consistently ca. 4.40 eV, while the threshold to lose the second hydrogen atom is much lower at ca. 3.16 eV. The only exception was for the dibenzo[a,l]pyrene cation, which has a unique structure due to steric constraints, resulting in a low H loss reaction energy of 2.85 eV. If C2H2 is lost directly from the precursor cation, the energy required for this dissociation is 4.16 eV. No other fragmentation channels were observed over a large enough sample set for trends to be extrapolated, though data on CH3 and C4H2 loss obtained in previous studies is included for completeness. The dissociation reactions were also studied by collision induced dissociation after ionization by atmospheric pressure chemical ionization. When modeled with a simple temperature-based theory for the post-collision internal energy distribution, there was reasonable agreement between the two sets of data.

Using radio astronomical receivers for molecular spectroscopic characterization in astrochemical laboratory simulations: A proof of concept

We present a proof of concept on the coupling of radio astronomical receivers and spectrometers with chemical reactors and the performances of the resulting setup for spectroscopy and chemical simulations in laboratory astrophysics. Several experiments including cold plasma generation and UV photochemistry were performed in a 40 cm long gas cell placed in the beam path of the Aries 40 m radio telescope receivers operating in the 41-49 GHz frequency range interfaced with fast Fourier transform spectrometers providing 2 GHz bandwidth and 38 kHz resolution. The impedance matching of the cell windows has been studied using different materials. The choice of the material and its thickness was critical to obtain a sensitivity identical to that of standard radio astronomical observations. Spectroscopic signals arising from very low partial pressures of CH₃OH, CH₃CH₂OH, HCOOH, OCS, CS, SO₂ (<10^-3 mbar) were detected in a few seconds. Fast data acquisition was achieved allowing for kinetic measurements in fragmentation experiments using electron impact or UV irradiation. Time evolution of chemical reactions involving OCS, O₂ and CS₂ was also observed demonstrating that reactive species, such as CS, can be maintained with high abundance in the gas phase during these experiments.

The ESO Diffuse Interstellar Bands Large Exploration Survey: EDIBLES I. Project description, survey sample and quality assessment

The carriers of the diffuse interstellar bands (DIBs) are largely unidentified molecules ubiquitously present in the interstellar medium (ISM). After decades of study, two strong and possibly three weak near-infrared DIBs have recently been attributed to the [Formula: see text] fullerene based on observational and laboratory measurements. There is great promise for the identification of the over 400 other known DIBs, as this result could provide chemical hints towards other possible carriers. In an effort to systematically study the properties of the DIB carriers, we have initiated a new large-scale observational survey: the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES). The main objective is to build on and extend existing DIB surveys to make a major step forward in characterising the physical and chemical conditions for a statistically significant sample of interstellar lines-of-sight, with the goal to reverse-engineer key molecular properties of the DIB carriers. EDIBLES is a filler Large Programme using the Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope at Paranal, Chile. It is designed to provide an observationally unbiased view of the presence and behaviour of the DIBs towards early-spectral type stars whose lines-of-sight probe the diffuse-to-translucent ISM. Such a complete dataset will provide a deep census of the atomic and molecular content, physical conditions, chemical abundances and elemental depletion levels for each sightline. Achieving these goals requires a homogeneous set of high-quality data in terms of resolution (R ~ 70 000 - 100 000), sensitivity (S/N up to 1000 per resolution element), and spectral coverage (305-1042 nm), as well as a large sample size (100+ sightlines). In this first paper the goals, objectives and methodology of the EDIBLES programme are described and an initial assessment of the data is provided.

Detection of Buckminsterfullerene emission in the diffuse interstellar medium

Emission of fullerenes in their infrared vibrational bands has been detected in space near hot stars. The proposed attribution of the diffuse interstellar bands at 9577 and 9632 Å to electronic transitions of the buckminsterfullerene cation (i.e. [Formula: see text]) was recently supported by new laboratory data, confirming the presence of this species in the diffuse interstellar medium (ISM). In this letter, we present the detection, also in the diffuse ISM, of the 17.4 and 18.9 μm emission bands commonly attributed to vibrational bands of neutral C60. According to classical models that compute the charge state of large molecules in space, C60 is expected to be mostly neutral in the diffuse ISM. This is in agreement with the abundances of diffuse C60 we derive here from observations. We also find that C60 is less abundant in the diffuse ISM than in star-forming regions, supporting the theory that C60 can be formed in these regions.

Size Effect in the Ionization Energy of PAH Clusters

We report the first experimental measurement of the near-threshold photoionization spectra of polycyclic aromatic hydrocarbon clusters made of pyrene C16H10 and coronene C24H12, obtained using imaging photoelectron-photoion coincidence spectrometry with a VUV synchrotron beamline. The experimental results of the ionization energy are compared to calculated ones obtained from simulations using dedicated electronic structure treatment for large ionized molecular clusters. Experiment and theory consistently find a decrease of the ionization energy with cluster size. The inclusion of temperature effects in the simulations leads to a lowering of this energy and to quantitative agreement with the experiment. In the case of pyrene, both theory and experiment show a discontinuity in the IE trend for the hexamer. This work demonstrates the ability of the models to describe the electronic structure of PAH clusters and suggests that these species are ionized in astronomical environments where they are thought to be present.

Herschel survey and modelling of externally-illuminated photoevaporating protoplanetary disks

CONTEXT

Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far.

AIMS

We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons (i.e. 6 < E < 13.6 eV).

METHODS

We present a unique Herschel survey and new ALMA observations of four externally-illuminated photoevaporating disks (a.k.a. proplyds). For the analysis of these data, we developed a 1D model of the photodissociation region (PDR) of a proplyd, based on the Meudon PDR code and we computed the far infrared line emission.

RESULTS

With this model, we successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about 106 cm-3 and local gas temperatures of about 1000 K. Our modelling suggests that all studied disks are found in a transitional regime resulting from the interplay between several heating and cooling processes that we identify. These differ from those dominating in classical PDRs i.e. grain photo-electric effect and cooling by [OI] and [CII] FIR lines. This specific energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to keep the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From the physical parameters derived from our best-fit models, we estimate mass-loss rates - of the order of 10-7 M⊙/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume photoevaporation in the supercritical regime where the evaporation flow is launched from the disk surface at sound speed.

CONCLUSIONS

We have identified the energetic regime regulating FUV-photoevaporation in proplyds. This regime could be implemented into models of the dynamical evolution of protoplanetary disks.

Identification of PAH Isomeric Structure in Cosmic Dust Analogues: the AROMA setup

We developed a new analytical experimental setup called AROMA (Astrochemistry Research of Organics with Molecular Analyzer) that combines laser desorption/ionization techniques with ion trap mass spectrometry. We report here on the ability of the apparatus to detect aromatic species in complex materials of astrophysical interests and characterize their structures. A limit of detection of 100 femto-grams has been achieved using pure polycyclic aromatic hydrocarbon (PAH) samples, which corresponds to 2x108 molecules in the case of coronene (C24H12). We detected the PAH distribution in the Murchison meteorite, which is made of a complex mixture of extraterrestrial organic compounds. In addition, collision induced dissociation experiments were performed on selected species detected in Murchison, which led to the first firm identification of pyrene and its methylated derivatives in this sample.

The growth of carbon chains in IRC +10216 mapped with ALMA

Linear carbon chains are common in various types of astronomical molecular sources. Possible formation mechanisms involve both bottom-up and top-down routes. We have carried out a combined observational and modeling study of the formation of carbon chains in the C-star envelope IRC +10216, where the polymerization of acetylene and hydrogen cyanide induced by ultraviolet photons can drive the formation of linear carbon chains of increasing length. We have used ALMA to map the emission of λ 3 mm rotational lines of the hydrocarbon radicals C2H, C4H, and C6H, and the CN-containing species CN, C3N, HC3N, and HC5N with an angular resolution of ~1″. The spatial distribution of all these species is a hollow, 5-10″ wide, spherical shell located at a radius of 10-20″ from the star, with no appreciable emission close to the star. Our observations resolve the broad shell of carbon chains into thinner sub-shells which are 1-2″ wide and not fully concentric, indicating that the mass loss process has been discontinuous and not fully isotropic. The radial distributions of the species mapped reveal subtle differences: while the hydrocarbon radicals have very similar radial distributions, the CN-containing species show more diverse distributions, with HC3N appearing earlier in the expansion and the radical CN extending later than the rest of the species. The observed morphology can be rationalized by a chemical model in which the growth of polyynes is mainly produced by rapid gas-phase chemical reactions of C2H and C4H radicals with unsaturated hydrocarbons, while cyanopolyynes are mainly formed from polyynes in gas-phase reactions with CN and C3N radicals.

Spatial distribution of FIR rotationally excited CH+ and OH emission lines in the Orion Bar PDR

CONTEXT

The methylidyne cation (CH⁺) and hydroxyl (OH) are key molecules in the warm interstellar chemistry, but their formation and excitation mechanisms are not well understood. Their abundance and excitation are predicted to be enhanced by the presence of vibrationally excited H₂ or hot gas (~500-1000 K) in photodissociation regions with high incident FUV radiation field. The excitation may also originate in dense gas (> 10⁵ cm^-3) followed by nonreactive collisions with H₂, H, and electrons. Previous observations of the Orion Bar suggest that the rotationally excited CH⁺ and OH correlate with the excited CO, a tracer of dense and warm gas, and formation pumping contributes to CH⁺ excitation.

AIMS

Our goal is to examine the spatial distribution of the rotationally excited CH⁺ and OH emission lines in the Orion Bar in order to establish their physical origin and main formation and excitation mechanisms.

METHODS

We present spatially sampled maps of the CH⁺ J=3-2 transition at 119.8 µm and the OH Λ-doublet at 84 µm in the Orion Bar over an area of 110″×110″ with Herschel (PACS). We compare the spatial distribution of these molecules with those of their chemical precursors, C⁺, O and H₂, and tracers of warm and dense gas (high-J CO). We assess the spatial variation of CH⁺ J=2-1 velocity-resolved line profile at 1669 GHz with Herschel HIFI spectrometer observations.

RESULTS

The OH and especially CH⁺ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region at the edge of the Bar. While notably similar, the differences in the CH⁺ and OH morphologies indicate that CH⁺ formation and excitation are strongly related to the observed vibrationally excited H₂. This, together with the observed broad CH⁺ line widths, indicates that formation pumping contributes to the excitation of this reactive molecular ion. Interestingly, the peak of the rotationally excited OH 84 µm emission coincides with a bright young object, proplyd 244-440, which shows that OH can be an excellent tracer of UV-irradiated dense gas.

CONCLUSIONS

The spatial distribution of CH⁺ and OH revealed in our maps is consistent with previous modeling studies. Both formation pumping and nonreactive collisions in a UV-irradiated dense gas are important CH⁺ J=3-2 excitation processes. The excitation of the OH Λ-doublet at 84 µm is mainly sensitive to the temperature and density.

Cooling of isolated anthracene cations probed with photons of different wavelengths in the Mini-Ring

We report on a direct measurement of the Internal Energy Distribution (IED) shift rate of an initially hot polycyclic aromatic hydrocarbon (PAH) molecular ensemble, anthracene cations (C₁₄H₁₀⁺). The ions were produced in an electron cyclotron resonance (ECR) ion source and stored in an electrostatic ion storage ring, the Mini-Ring. Laser pulses of two wavelengths were sent successively to merge the stored ion bunch at different storage times to enhance the neutral fragment yield due to fast laser induced dissociation. Using this technique, we have been able to determine directly the energy shift rate of the IED, without involving any theoretical simulation or any assumption on dissociation rates, cooling rates, or the initial IED. Theoretical energy shift rates have been estimated from the evolution of simulated IEDs by taking into account the effects of the unimolecular dissociation and two radiative decay mechanisms: the Poincaré fluorescence and the infrared vibrational emission. The comparison between the experimental results and the model provides new evidence of the important role of the Poincaré fluorescence in the overall cooling process of anthracene cations. Although in the short time range the commonly accepted intuition says that the cooling would result mostly from the dissociation of the hottest ions (depletion cooling), we demonstrate that the Poincaré fluorescence is the dominant contribution (about 85%) to the net cooling effect.

Compression and ablation of the photo-irradiated molecular cloud the Orion Bar

The Orion Bar is the archetypal edge-on molecular cloud surface illuminated by strong ultraviolet radiation from nearby massive stars. Owing to the close distance to Orion (about 1,350 light-year), the effects of stellar feedback on the parental cloud can be studied in detail. Visible-light observations of the Bar1 show that the transition between the hot ionised gas and the warm neutral atomic gas (the ionisation front) is spatially well separated from the transition from atomic to molecular gas (the dissociation front): about 15 arcseconds or 6,200 astronomical units. (One astronomical unit is the Earth-Sun distance.) Static equilibrium models2,3 used to interpret previous far-infrared and radio observations of the neutral gas in the Bar4,5,6 (typically at 10-20 arcsecond resolution) predict an inhomogeneous cloud structure consisting of dense clumps embedded in a lower density extended gas component. Here we report 1 arcsecond resolution millimetre-wave images that allow us to resolve the molecular cloud surface and constrain the gas density and temperature structures at small spatial scales. In contrast to stationary model predictions7,8,9, there is no appreciable offset between the peak of the H₂ vibrational emission (delineating the H/H₂ transition) and the edge of the observed CO and HCO⁺ emission. This implies that the H/H₂ and C⁺/C/CO transition zones are very close. These observations reveal a fragmented ridge of high-density substructures, photo-ablative gas flows and instabilities at the molecular cloud surface. They suggest that the cloud edge has been compressed by a high-pressure wave that currently moves into the molecular cloud. The images demonstrate that dynamical and nonequilibrium effects are important. Thus, they should be included in any realistic description of irradiated interstellar matter.

VUV photo-processing of PAH cations: quantitative study on the ionization versus fragmentation processes

Interstellar polycyclic aromatic hydrocarbons (PAHs) are strongly affected by the absorption of vacuum ultraviolet (VUV) photons in the interstellar medium (ISM), yet the branching ratio between ionization and fragmentation is poorly studied. This is crucial for the stability and charge state of PAHs in the ISM in different environments, affecting in turn the chemistry, the energy balance, and the contribution of PAHs to the extinction and emission curves. We studied the interaction of PAH cations with VUV photons in the 7 - 20 eV range from the synchrotron SOLEIL beamline, DESIRS. We recorded by action spectroscopy the relative intensities of photo-fragmentation and photo-ionization for a set of eight PAH cations ranging in size from 14 to 24 carbon atoms, with different structures. At photon energies below ~13.6 eV fragmentation dominates for the smaller species, while for larger species ionization is immediately competitive after the second ionization potential (IP). At higher photon energies, all species behave similarly, the ionization yield gradually increases, leveling off between 0.8 and 0.9 at ~18 eV. Among isomers, PAH structure appears to mainly affect the fragmentation cross section, but not the ionization cross section. We also measured the second IP for all species and the third IP for two of them, all are in good agreement with theoretical ones confirming that PAH cations can be further ionized in the diffuse ISM. Determining actual PAH dication abundances in the ISM will require detailed modeling. Our measured photo-ionization yields for several PAH cations provide a necessary ingredient for such models.

An optical spectrum of a large isolated gas-phase PAH cation: C78H26

A gas-phase optical spectrum of a large polycyclic aromatic hydrocarbon (PAH) cation - C78H26+- in the 410-610 nm range is presented. This large all-benzenoid PAH should be large enough to be stable with respect to photodissociation in the harsh conditions prevailing in the interstellar medium (ISM). The spectrum is obtained via multi-photon dissociation (MPD) spectroscopy of cationic C78H26 stored in the Fourier Transform Ion Cyclotron Resonance (FT-ICR) cell using the radiation from a mid-band optical parametric oscillator (OPO) laser. The experimental spectrum shows two main absorption peaks at 431 nm and 516 nm, in good agreement with a theoretical spectrum computed via time-dependent density functional theory (TD-DFT). DFT calculations indicate that the equilibrium geometry, with the absolute minimum energy, is of lowered, nonplanar C2 symmetry instead of the more symmetric planar D2h symmetry that is usually the minimum for similar PAHs of smaller size. This kind of slightly broken symmetry could produce some of the fine structure observed in some diffuse interstellar bands (DIBs). It can also favor the folding of C78H26+ fragments and ultimately the formation of fullerenes. This study opens up the possibility to identify the most promising candidates for DIBs amongst large cationic PAHs.

Polycyclic aromatic hydrocarbons and molecular hydrogen in oxygen-rich planetary nebulae: the case of NGC 6720

Evolved stars are primary sources for the formation of polycyclic aromatic hydrocarbons (PAHs) and dust grains. Their circumstellar chemistry is usually designated as either oxygen-rich or carbon-rich, although dual-dust chemistry objects, whose infrared spectra reveal both silicate- and carbon-dust features, are also known. The exact origin and nature of this dual-dust chemistry is not yet understood. Spitzer-IRS mid-infrared spectroscopic imaging of the nearby, oxygen-rich planetary nebula NGC 6720 reveals the presence of the 11.3 μm aromatic (PAH) emission band. It is attributed to emission from neutral PAHs, since no band is observed in the 7-8 μm range. The spatial distribution of PAHs is found to closely follow that of the warm clumpy molecular hydrogen emission. Emission from both neutral PAHs and warm H2 is likely to arise from photo-dissociation regions associated with dense knots that are located within the main ring. The presence of PAHs together with the previously derived high abundance of free carbon (relative to CO) suggest that the local conditions in an oxygen-rich environment can also become conducive to in-situ formation of large carbonaceous molecules, such as PAHs, via a bottom-up chemical pathway. In this scenario, the same stellar source can enrich the interstellar medium with both oxygen-rich dust and large carbonaceous molecules.

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.

Mixed aliphatic and aromatic composition of evaporating very small grains in NGC 7023 revealed by the 3.4/3.3 μm ratio

CONTEXT

A chemical scenario was proposed for photon-dominated regions (PDRs) according to which UV photons from nearby stars lead to the evaporation of very small grains (VSGs) and the production of gas-phase polycyclic aromatic hydrocarbons (PAHs).

AIMS

Our goal is to achieve better insight into the composition and evolution of evaporating very small grains (eVSGs) and PAHs through analyzing the infrared (IR) aliphatic and aromatic emission bands.

METHODS

We combined spectro-imagery in the near- and mid-IR to study the spatial evolution of the emission bands in the prototypical PDR NGC 7023. We used near-IR spectra obtained with the IRC instrument onboard AKARI to trace the evolution of the 3.3 μm and 3.4 μm bands, which are associated with aromatic and aliphatic C-H bonds on PAHs. The spectral fitting involved an additional broad feature centered at 3.45 μm that is often referred to as the plateau. Mid-IR observations obtained with the IRS instrument onboard the Spitzer Space Telescope were used to distinguish the signatures of eVSGs and neutral and cationic PAHs. We correlated the spatial evolution of all these bands with the intensity of the UV field given in units of the Habing field G0 to explore how their carriers are processed.

RESULTS

The intensity of the 3.45 μm plateau shows an excellent correlation with that of the 3.3 μm aromatic band (correlation coefficient R = 0.95) and a relatively poor correlation with the aliphatic 3.4 μm band (R=0.77). This indicates that the 3.45 μm feature is dominated by the emission from aromatic bonds. We show that the ratio of the 3.4 μm and 3.3 μm band intensity (I3.4/I3.3) decreases by a factor of 4 at the PDR interface from the more UV-shielded layers (G0 ~ 150, I3.4/I3.3 = 0.13) to the more exposed layers (G0 > 1 × 104, I3.4/I3.3 = 0.03). The intensity of the 3.3 μm band relative to the total neutral PAH intensity shows an overall increase with G0, associated with an increase of both the hardness of the UV field and the H abundance. In contrast, the intensity of the 3.4 μm band relative to the total neutral PAH intensity decreases with G0, showing that their carriers are actively destroyed by UV irradiation and are not efficiently regenerated. The transition region between the aliphatic and aromatic material is found to correspond spatially with the transition zone between neutral PAHs and eVSGs.

CONCLUSIONS

We conclude that the photo-processing of eVSGs leads to the production of PAHs with attached aliphatic sidegroups that are revealed by the 3.4 μm emission band. Our analysis provides evidence for the presence of very small grains of mixed aromatic and aliphatic composition in PDRs.

LABORATORY PHOTO-CHEMISTRY OF PAHS: IONIZATION VERSUS FRAGMENTATION

Interstellar Polycyclic Aromatic Hydrocarbons (PAH) are expected to be strongly processed by Vacuum Ultra-Violet (VUV) photons. Here, we report experimental studies on the ionization and fragmentation of coronene (C24H12), ovalene (C32H14) and hexa-peri-hexabenzocoronene (HBC; C42H18) cations by exposure to synchrotron radiation in the range of 8-40 eV. The results show that for small PAH cations such as coronene, fragmentation (H-loss) is more important than ionization. However, as the size increases, ionization becomes more and more important and for the HBC cation, ionization dominates. These results are discussed and it is concluded that, for large PAHs, fragmentation only becomes important when the photon energy has reached the highest ionization potential accessible. This implies that PAHs are even more photo-stable than previously thought. The implications of this experimental study for the photo-chemical evolution of PAHs in the interstellar medium (ISM) are briefly discussed.

Top-down formation of fullerenes in the interstellar medium

Fullerenes have been recently detected in various circumstellar and interstellar environments, raising the question of their formation pathway. It has been proposed that they can form at the low densities found in the interstellar medium by the photo-chemical processing of large polycyclic aromatic hydrocarbons (PAHs). Following our previous work on the evolution of PAHs in the NGC 7023 reflection nebula, we evaluate, using photochemical modeling, the possibility that the PAH C66H20 (i.e. circumovalene) can lead to the formation of C60 upon irradiation by ultraviolet photons. The chemical pathway involves full dehydrogenation of C66H20, folding into a floppy closed cage and shrinking of the cage by loss of C2 units until it reaches the symmetric C60 molecule. At 10" from the illuminating star and with realistic molecular parameters, the model predicts that 100% of C66H20 is converted into C60 in ~ 105 years, a timescale comparable to the age of the nebula. Shrinking appears to be the kinetically limiting step of the whole process. Hence, PAHs larger than C66H20 are unlikely to contribute significantly to the formation of C60, while PAHs containing between 60 and 66 C atoms should contribute to the formation of C60 with shorter timescales, and PAHs containing less than 60 C atoms will be destroyed. Assuming a classical size distribution for the PAH precursors, our model predicts absolute abundances of C60 are up to several 10-4 of the elemental carbon, i.e. less than a percent of the typical interstellar PAH abundance, which is consistent with observational studies. According to our model, once formed, C60 can survive much longer (> 107 years for radiation fields below G0 = 104) than other fullerenes because of the remarkable stability of the C60 molecule at high internal energies. Hence, a natural consequence is that C60 is more abundant than other fullerenes in highly irradiated environments.

Dissociation of the anthracene radical cation: a comparative look at iPEPICO and collision-induced dissociation mass spectrometry results

The dissociation of the anthracene radical cation has been studied using two different methods: imaging photoelectron photoion coincidence spectrometry (iPEPCO) and atmospheric pressure chemical ionization-collision induced dissociation mass spectrometry (APCI-CID). Four reactions were investigated: (R1) C14H10(+•) → C14H9(+) + H, (R2) C14H9(+) → C14H8(+•) + H, (R3) C14H10(+•) → C12H8(+•) + C2H2 and (R4) C14H10(+•) → C10H8(+•) + C4H2. An attempt was made to assign structures to each fragment ion, and although there is still room for debate whether for the C12H8(+•) fragment ion is a cyclobuta[b]naphthalene or a biphenylene cation, our modeling results and calculations appear to suggest the more likely structure is cyclobuta[b]naphthalene. The results from the iPEPICO fitting of the dissociation of ionized anthracene are E0 = 4.28 ± 0.30 eV (R1), 2.71 ± 0.20 eV (R2), and 4.20 ± 0.30 eV (average of reaction R3) whereas the Δ(‡)S values (in J K(-1) mol(-1)) are 12 ± 15 (R1), 0 ± 15 (R2), and either 7 ± 10 (using cyclobuta[b]naphthalene ion fragment in reaction R3) or 22 ± 10 (using the biphenylene ion fragment in reaction R3). Modeling of the APCI-CID breakdown diagrams required an estimate of the postcollision internal energy distribution, which was arbitrarily assumed to correspond to a Boltzmann distribution in this study. One goal of this work was to determine if this assumption yields satisfactory energetics in agreement with the more constrained and theoretically vetted iPEPICO results. In the end, it did, with the APCI-CID results being similar.