Structures and stabilities of PAH clusters solvated by water aggregates: The case of the pyrene dimer

Although clusters made of polycyclic aromatic hydrocarbon and water monomers are relevant objects in both atmospheric and astrophysical science, little is known about their energetic and structural properties. In this work, we perform global explorations of the potential energy landscapes of neutral clusters made of two pyrene units and one to ten water molecules using a density-functional-based tight-binding (DFTB) potential followed by local optimizations at the density-functional theory level. We discuss the binding energies with respect to various dissociation channels. It shows that cohesion energies of the water clusters interacting with a pyrene dimer are larger than those of the pure water clusters, reaching for the largest clusters an asymptotic limit similar to that of pure water clusters and that, although the hexamer and octamer can be considered magic numbers for isolated water clusters, it is not the case anymore when they are interacting with a pyrene dimer. Ionization potentials are also computed by making use of the configuration interaction extension of DFTB, and we show that in cations, the charge is mostly carried by the pyrene molecules.

Infrared Spectroscopy and Photochemistry of Anthracoronene in Cosmic Water Ice

We present a laboratory study of the polycyclic aromatic hydrocarbon (PAH) anthracoronene (AntCor, C₃₆H₁₈) in simulated interstellar ices in order to determine its possible contribution to the broad infrared absorption bands in the 5-8 μm wavelength interval. The Fourier transform infrared (FTIR) spectrum of AntCor, codeposited with water ice, was collected. The FTIR spectrum of the sample irradiated with ultraviolet photons was also collected. Unirradiated and UV-irradiated AntCor embedded in water ice have not been studied before; therefore, the molecule's band positions and intensities were compared to published data on AntCor in an argon matrix and theoretical calculations (DFT), as well as the published results of its parent molecules, coronene and anthracene, in water ice. The experimental band strengths for unirradiated AntCor exhibit variability as a function of PAH:H₂O concentration, with two distinct groupings of band intensities. AntCor clustering occurs for all concentrations and has a significant effect on PAH degradation rates and photoproduct variability. Near-IR spectra of irradiated AntCor samples show that AntCor⁺ production increases as the concentration of AntCor in water ice decreases. Photoproduct bands are assigned to AntCor⁺, cationic alcohols, protonated AntCor, and ketones. We report the rate constants of the photoproduct production for the 1:1280 AntCor:H₂O concentration. CO₂ production from AntCor is much less than what was previously reported for Ant and Cor and exhibits two distinct regimes as a function of AntCor:H₂O concentration. The contribution of AntCor photoproducts to astronomical spectra can be estimated by comparison with the observed intensities in the 7.4-8.0 μm range.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

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-(H₂O)_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.

Microhydration of PAH+ cations: evolution of hydration network in naphthalene+-(H2O) n clusters (n ≤ 5)

The interaction of polycyclic aromatic hydrocarbon molecules with water (H2O = W) is of fundamental importance in chemistry and biology. Herein, size-selected microhydrated naphthalene cation nanoclusters, Np+-W n (n ≤ 5), are characterized by infrared photodissociation (IRPD) spectroscopy in the C-H and O-H stretch range to follow the stepwise evolution of the hydration network around this prototypical PAH+ cation. The IRPD spectra are highly sensitive to the hydration structure and are analyzed by dispersion-corrected density functional theory calculations (B3LYP-D3/aug-cc-pVTZ) to determine the predominant structural isomers. For n = 1, W forms a bifurcated CH···O ionic hydrogen bond (H-bond) to two acidic CH protons of the bicyclic ring. For n ≥ 2, the formation of H-bonded solvent networks dominates over interior ion solvation, because of strong cooperativity in the former case. For n ≥ 3, cyclic W n solvent structures are attached to the CH protons of Np+. However, while for n = 3 the W3 ring binds in the CH···O plane to Np+, for n ≥ 4 the cyclic W n clusters are additionally stabilized by stacking interactions, leading to sandwich-type configurations. No intracluster proton transfer from Np+ to the W n solvent is observed in the studied size range (n ≤ 5), because of the high proton affinity of the naphthyl radical compared to W n . This is different from microhydrated benzene+ clusters, (Bz-W n )+, for which proton transfer is energetically favorable for n ≥ 4 due to the much lower proton affinity of the phenyl radical. Hence, because of the presence of polycyclic rings, the interaction of PAH+ cations with W is qualitatively different from that of monocyclic Bz+ with respect to interaction strength, structure of the hydration shell, and chemical reactivity. These differences are rationalized and quantified by quantum chemical analysis using the natural bond orbital (NBO) and noncovalent interaction (NCI) approaches.

Matrix Isolation Spectroscopy: Aqueous p-Toluenesulfonic Acid Solvation

Interaction between p-toluenesulfonic acid (pTSA) and water is studied at -20 °C in a CCl₄ matrix. In CCl₄ water exists as monomers with restricted rotational motion about its symmetry axis. Additionally, CCl₄ is transparent in the hydrogen-bonded region; CCl₄ thus constitutes an excellent ambient thermal energy matrix isolation medium for diagnosing interactions with water. Introducing pTSA-nH₂O gives rise to two narrow resonances at 3642 cm^-1 and at 2835 cm^-1 plus a broad 3000-3550 cm^-1 absorption. In addition, negative monomer symmetric and asymmetric stretch features relative to nominally dry CCl₄ indicate that fewer water monomers exist in the cooled (-20 °C) acid solution than in room-temperature anhydrous CCl₄. The negative peaks along with the broad absorption band indicate that water monomers are incorporated into clusters. The 3642 cm^-1 resonance is assigned to the OH-π interaction with a cluster containing many water molecules per acid molecule. The 2835 cm^-1 resonance is assigned to the (S-)O-H stretch of pTSA-dihydrate. The coexistence of these two species provides insights into interactions in this acid-water CCl₄ system.

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:H₂O adducts is of interest in fields as diverse as catalysis and astrochemistry. Industrially, iron can be used as a catalyst to convert H₂O to H₂, 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 H₂O 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 HFe₂OH, leading in some cases to the formation of FeO possibly accompanied by the production of H₂. 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.