# Water Clusters in Interaction with Corannulene in a Rare Gas Matrix: Structures, Stability and IR Spectra

Héloïse Leboucher, Joëlle Mascetti, Christian Aupetit, Jennifer Noble, Aude Simon

► To cite this version:

Héloïse Leboucher, Joëlle Mascetti, Christian Aupetit, Jennifer Noble, Aude Simon. Water Clusters in Interaction with Corannulene in a Rare Gas Matrix: Structures, Stability and IR Spectra. Photochem. 2022, 2 (2), pp.237-262. 10.3390/photochem2020018. hal-03689280

HAL Id: hal-03689280

https://hal.science/hal-03689280

Submitted on 12 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Water Clusters in Interaction with Corannulene in a Rare Gas Matrix: Structures, Stability and IR Spectra

Heloïse Leboucher<fnref n="1" />,*, Joëlle Mascetti<fnref n="2" />,†, Christian Aupetit<fnref n="2" />, Jennifer A. Noble<fnref n="3" />,† and Aude Simon <fnref n="1" />,* †

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Leboucher, H.; Mascetti, J.; Aupetit, C.; Noble, J.A.; Simon, A. Water Clusters in Interaction with Corannulene in a Rare Gas Matrix: Structures, Stability and IR Spectra.

*Photochem. 2022*, *2*, 237–262.

https://doi.org/10.3390/photochem2020018

Received: 1 March 2022

Accepted: 23 March 2022

Published: 25 March 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** The interaction of polycyclic aromatic hydrocarbons (PAHs) with water is of paramount importance in atmospheric and astrophysical contexts. We report here a combined theoretical and experimental study of corannulene-water interactions in low temperature matrices and of the matrix's influence on the photoreactivity of corannulene with water. The theoretical study was performed using a mixed density functional based tight binding/force field approach to describe the corannulene-water clusters trapped in an argon matrix, together with Born-Oppenheimer molecular dynamics to determine finite-temperature IR spectra. The results are discussed in the light of experimental matrix isolation FTIR spectroscopic data. We show that in the solid phase, $\pi$ isomers of $(C_{20}H_{10})(H_2O)_n$, with $n = 2$ or 3, are energetically favored. These $\pi$ complexes are characterized by small shifts in corannulene vibrational modes and large shifts in water bands. These $\pi$ structures, particularly stable in the case of the water trimer where the water cluster is trapped "inside" the corannulene bowl, may account for the difference in photoreactivity of non-planar-compared to planar-PAHs with water. Indeed, planar PAHs such as pyrene and coronene embedded in H₂O:Ar matrices form $\sigma$ isomers and react with water to form alcohols and quinones under low energy UV irradiation, whereas no photoreactivity was observed for corannulene under the same experimental conditions.

**Keywords:** PAHs; astrochemistry; matrix isolation FTIR spectroscopy; water clusters; corannulene; DFTB; BOMD

## 1. Introduction

The interaction of polycyclic aromatic hydrocarbons (PAHs) with water is of paramount importance in a large number of fields, such as the petroleum industry [1], atmospheric science [2,3], astrophysics and astrochemistry. Due to its allotropy, a large fraction of interstellar carbon is expected to be included in large molecules or in dust grains with various degrees of organization and aromaticity, including fullerenes, diamondoids, PAHs, clusters of PAHs, or amorphous carbon grains, possibly hydrogenated [4]. AstroPAHs in particular are believed to be present in the interstellar medium (ISM), where they account for up to 20% of the total carbon [5,6]. PAHs have been of significant interest since the proposal in the eighties that they were the carriers of the Aromatic Infrared Bands (AIBs), a set of mid-IR emission bands observed in many regions of the ISM [7,8]. This proposal led to many experimental and theoretical spectroscopic studies in order to identify a specific PAH molecule [6]. However, it is only very recently that specific PAH molecules, the two isomers of cyano-naphthalene, have been successfully detected based on their rotational spectra [9].

In interstellar cold dense molecular clouds, where the temperature ranges from 10 to 50 K, PAHs are expected to be frozen on or in icy mantles of dust grains [10]. These ice

<footnote n="1">¹ Laboratoire de Chimie et Physique Quantiques (LCPQ), Fédération FeRMI, UMR5626, Université de Toulouse [UT3] & CNRS, University Paul Sabatier, 118 Route de Narbonne, 31062 Toulouse, France; heloise.leboucher@irsamc.ups-tlse.fr</footnote>

<footnote n="2">² Institut des Sciences Moléculaires (ISM), Université de Bordeaux, UMR5255 CNRS, 351 Cours de la Libération, CEDEX, 33405 Talence, France; joelle.mascetti@u-bordeaux.fr (J.M.); christian.aupetit@u-bordeaux.fr (C.A.)</footnote>

<footnote n="3">³ Laboratoire Physique des Interactions Ioniques et Moléculaires (PIIM), UMR7345, CNRS, Aix-Marseille Université, 13013 Marseille, France; jennifer.noble@univ-amu.fr</footnote>

<footnote>* Correspondence: aude.simon@irsamc.ups-tlse.fr</footnote>
mantles are mainly composed of amorphous solid water (ASW), and allow chemical
reactions to occur after heating or irradiation [11]. It has been shown that PAHs embedded
in ASW, when irradiated by VUV light, undergo erosion processes and lead to oxidation
products such as alcohols and ketones [12] via the formation of cationic intermediates.
Interestingly, the formation of oxygenated PAHs is also observed in ASW when using
low energy UV-visible irradiation [13–16]. It has also been shown that small complexes of
coronene and water isolated in cryogenic argon matrices exhibit photochemical reactions
similar to those observed in ASW [17]. We have previously demonstrated that, when
embedded in a solid matrix, coronene forms σ complexes with water molecules arranged
in interaction with the edge of the planar PAH [18]. Such σ interaction is expected to favor
photochemical reaction through better H/charge transfer to form oxidation products from
planar PAHs [17,19]. These results suggest that a relatively easy formation of oxygenated
PAH photoproducts could account for shifts observed in the AIBs.

In this work, we present some experimental and theoretical results concerning corannulene C₂₀H₁₀ in interaction with water clusters embedded in a rare gas matrix. Corannulene is a fully aromatic non-planar PAH [20] of C₅v symmetry, consisting of a five membered carbon cycle surrounded by five six-membered carbon cycles that can be described as a bowl-like molecule with two faces: convex and concave (see Figure 1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Scheme of corannulene C₂₀H₁₀, a hydrogenated cap of buckminsterfullerene of C₅v symmetry, top and side views.</figcaption></figure>

Corannulene can be regarded as a cap of buckminsterfullerene C₆₀—a molecule that was identified in Space more than 10 years ago [21,22]—saturated with H atoms. It can thus reasonably be considered as a possible intermediate in “top-down” or “bottom-up” routes leading to C₆₀ formation. The mechanism for the formation of buckminsterfullerene C₆₀ is still under debate. Following its discovery, observational studies suggested that fullerenes could be produced by photochemical dehydrogenation and isomerisation of amorphous carbon grains [23–26], but the composition of such grains remains elusive. Among current hypotheses, buckminsterfullerene could be formed following decay and rearrangement of dehydrogenated PAHs [27] under the action of UV photons [28] or by collisions with energetic particles [29]. Such a “top-down” mechanism is expected to prevail in low-density environments such as interstellar clouds [30] and the corannulene molecule can be regarded as a possible intermediate in such a sequence. Alternatively, a “bottom-up” mechanism, consisting of the successive growth from smaller carbon chain building blocks, known as the closed network growth mechanism [31], would prevail in the hot, dense envelopes of evolved stars [30]. The need to characterize the structural diversity of carbon clusters, possibly hydrogenated and with amorphous character, that could play a role in the formation of buckminsterfullerene, motivated the development of dedicated experimental and theoretical approaches in order to obtain their spectral features in the UV-visible [32,33] and IR domains [34–37].

Because of its bowl shape, corannulene could be detected in astrophysical environments by radioastronomy, as put forward by Pilleri et al. [38]. The astrophysical interest of corannulene has also motivated experimental and theoretical spectroscopic studies in the IR, Raman and UV-visible domains [39], both in the solid phase and, in particular, in rare
gas matrices [39], in order to propose spectroscopic characteristics that should be searched for in space. Recently, its formation in the gas phase was achieved by Zhao et al. [40] using a "bottom-up" approach involving the reaction of 7-fluoranthenyl (C₁₆H₉) and benzo[ghi]fluoranthen-5-yl (C₁₈H₉) radicals with acetylene (C₂H₂) in conditions mimicking those of carbon-rich circumstellar envelopes. Interestingly, due to its bowl shape that can invert and form host-guest complexes, corannulene also presents a particularly interesting building block in materials chemistry [41].

Studies on the structures, energetics and spectroscopic signatures of corannulene-water complexes are, to the best of our knowledge, scarce. Theoretical studies of water clusters (H₂O)ₙ in interaction with [C₂₀H₁₀]⁺ were performed by Hernández-Rojas et al. [42]. Global optimization using a force field approach was used to determine the geometries of the larger clusters. They showed that, for n = 1 and 2, the water cluster preferentially binds to the side of [C₂₀H₁₀]⁺, whereas larger sized water clusters interact preferentially with the concave face of cationic corannulene. Gas phase neutral corannulene-water complexes have been studied by Pérez et al. [43] using rotational spectroscopy complemented by *ab initio* calculations. They have shown that in (C₂₀H₁₀)(H₂O), the water molecule resides inside the bowl-like structure of corannulene, where it rotates freely around the C₅ axis and that dispersion interactions represent the main contribution to the binding. We have decided to investigate (C₂₀H₁₀)(H₂O)ₙ complexes isolated in argon matrices at 10 K. Our results are compared to those obtained in the gas phase [43] and to our previous study on coronene:water complexes. We detail our theoretical and experimental approaches in Section 2, report our results in Section 3, and discuss the experimental results in light of the calculations in Section 4.

## 2. Materials and Methods

### 2.1. *Computational Methods*

In order to determine the structures, energetics and IR spectra of (C₂₀H₁₀)(H₂O)ₙ (n = 1–3) embedded in an argon matrix, (hereafter quoted as: Ar) we used the same approach as that previously applied to coronene-water clusters (C₂₄H₁₂)(H₂O)ₙ (n = 1, 2) [18]. All calculations were performed with the deMonNano code [44].

The electronic structure of (C₂₀H₁₀)(H₂O)ₙ is described at the self consistent charge density functional based tight binding (SCC-DFTB) level of theory [45] using the Hamiltonian that we modified so as to obtain a good description of PAH-water interactions and used to describe coronene:water aggregates in an Ar matrix [18]. Briefly, an empirical dispersion term was added and the Mulliken charges were replaced by Truhlar Charge Model 3 (CM3) charges [46] in order to improve the description of the polarisation of the bonds. The latter modification revealed important for the description of long range electrostatic interactions [47,48] and for IR spectra [49,50]. This modified SCC-DFTB scheme (denoted hereafter simply as DFTB) was benchmarked against wavefunction and DFT calculations [48,51]. It was shown to provide a good description of the structures and energetics of water clusters and PAHs in interaction with water clusters [48,51,52]. The values of the DFTB shifts induced by the adsorption of PAHs determined for the IR bands of water-ice dangling OH were shown to be in good agreement with experimental data [53]. Focusing more specifically on corannulene:water, when optimising the structure of (C₂₀H₁₀)(H₂O) at the DFTB level, we found that the most stable structure was similar to that found by Pérez et al. [43] in the gas phase, i.e., with the water molecule interacting with the concave face of corannulene. We computed a DFTB intermolecular binding energy of −21 kJ·mol⁻¹ (without zero-point energy) which is in very good agreement with the binding energy of −22.40 kJ·mol⁻¹ determined using symmetry adapted perturbation theory (SAPT) analysis [43]. Besides, additional comparison between gas phase DFTB and DFT structures, energetics and IR spectra for isolated corannulene:water clusters was made (see data base associated to the present paper), revealing similar trends as for benzene:water and coronene:water clusters [48].

The argon-argon interaction is described at the force field (FF) level using the Aziz potential [54]. The DFTB/FF coupling is incorporated via a scheme close to first order degenerate perturbation theory and developed to describe clusters of hydrocarbons and water
trapped inside a rare gas matrix [55,56]. This DFTB/FF scheme was carefully benchmarked for PAH:H₂O:Ar systems, i.e., involving the same interactions as in this work, regarding structures, energetics and IR spectra [18,55,56]. In this work, the Ar matrix is modeled by a large cluster (1139 Ar atoms for the total pure Ar matrix) organised in a face cubic centered (fcc) crystalline structure at the center of which the impurity ((C₂₀H₁₀)(H₂O)ₙ) is inserted after removing a minimal number x of Ar atoms (x is the vacancy size). In the rest of the manuscript, the total system will be referred to as (C₂₀H₁₀)(H₂O)ₙ:Ar.

To obtain the initial geometries, several symmetry planes were considered for inserting the corannulene molecule into the matrix. After determining the most stable situation for the embedded corannulene, all the possible interaction sites of the water molecule with C₂₀H₁₀ were considered. For n = 1, the substitution of one Ar atom by one water molecule was considered for all Ar atoms interacting with the concave and convex faces of corannulene as well as in the [111] insertion layer of corannulene when Ar atoms interact with the H atoms of corannulene. Once the most stable geometries for the three situations are determined, the substitution of all possible Ar atoms interacting with the water monomer in the most stable (C₂₀H₁₀)(H₂O)ₙ:Ar systems by another water molecule are considered. The most stable situations where the dimer interacts with the concave face, convex face, and “in plane” H atoms of corannulene are retained and the interaction of a third water molecule with the existing dimer is considered. Local geometry optimisations are then performed and the relative stabilities of the determined isomers are estimated by comparison of their substitution energy values (E<sub>sub</sub>, defined in ref. [18]) determined as follows: E<sub>sub</sub> is equal to the energy of the optimised total system (E(C₂₀H₁₀)(H₂O)ₙ:Ar) plus the energy of the vacancy (E(Ar<sub>x</sub>)) minus the energy of the optimised total matrix (E(:Ar)) minus the energy of the isolated monomers (E(C₂₀H₁₀) + nE(H₂O)), i.e.,

$$E_{sub} = E((C_{20}H_{10})(H_{2}O)_{n}:Ar) + E(Ar_{x}) - E(:Ar) - E(C_{20}H_{10}) - nE(H_{2}O) \quad (1)$$

The IR spectra of the most stable (lowest value of E<sub>sub</sub>) energy structures for all water cluster sizes in all three configurations in which the water cluster interacts with the concave face, convex face and “in plane” H atoms were first computed in the harmonic approximation. However, finite-temperature effects due to the anharmonicity of the PES may be expected even at low temperature. In this case, an *a priori* more suitable approach is to derive IR spectra from molecular dynamics (MD) simulations computing the Fourier transform of the dipole moment μ autocorrelation function

$$\alpha(\omega) \propto \omega^2 \int_{0}^{+\infty} dt \langle \mu(0) \cdot \mu(t) \rangle e^{i\omega t} \quad (2)$$

Therefore we used the same approach as in the work by Simon et al. [18], consisting in performing Born-Oppenheimer MD (BOMD) simulations in the microcanonical ensemble, computing the electronic structure on-the-fly with the DFTB/FF method. The IR absorption cross sections are determined using Equation (2) with the electrostatic dipole moment computed using the CM3 charges of the corannulene:water clusters, the Ar atoms remaining neutral. The positions of the Ar atoms at the surface of the cluster are frozen to ensure the absence of surface Ar atom motions in the molecular dynamics. The system was first thermalised at the desired temperature (10 K) by means of short simulations (10 ps) in the canonical ensemble using a Nose-Hoover chain of thermostats. These first (NVT) simulations allow the generation of the initial conditions for MD simulations performed in the microcanonical ensemble (NVE). For each corannulene-water cluster, 7 MD simulations of 100 ps in the (NVE) ensemble were run, using a time step of 0.1 fs. The IR spectra were computed for each simulation and averaged. In such conditions (low temperature), convergence is expected on positions but not necessarily on intensities [18], which was checked for the specific systems studied in the present work (see Appendix B, Figure A1). In addition, the quality of the dynamic IR spectra inside the Ar matrix was checked for the water monomer and dimer which, as will be discussed in Section 3.1.2, rotate in the matrix,
increasing the simulation time from 100 ps to 500 ps (see Appendix B, Figure A2). Due to the large number of large systems to be investigated consistently, we used the systematic procedure described hereabove (7 × 100 ps) for all of them.

## 2.2. *Experimental Methods*

Experiments were performed in a high vacuum experimental setup consisting of a stainless steel chamber with a base pressure of $10^{-7}$ mbar, containing a CsBr substrate cooled to 10 K by a closed-cycle He cryostat (Cryophysics Cryodine). The sample temperature was monitored by a Si diode thermometer positioned on the copper substrate holder in proximity to the CsBr window. Deionised ultrapure water was subjected to multiple freeze-pump-thaw cycles under vacuum to remove dissolved gases. Argon and water were premixed in a dosing line with concentrations of water in argon varying from 1:10 to 1:1 and injected into the chamber at a rate of 1 mL min$^{-1}$ in order to form different mixtures containing H$_2$O monomers, dimers, trimers, and larger clusters up to hexamers. Corannulene (TCI, 97%, used without purification) was sublimated by heating the powder to 120 °C in an oven within the high vacuum chamber and co-deposited with the argon:water mixture which was introduced as a diffusion jet flow via a gas nozzle inlet located 2 cm from the CsBr substrate. Depositions lasted for between one and two hours. Infrared spectra of corannulene:H$_2$O:Ar samples were recorded in transmission mode using a Bruker 70 V FTIR spectrometer with a DTGS detector from 4000 to 400 cm$^{-1}$ at a 0.5 cm$^{-1}$ resolution, with each spectrum averaged over 200 scans. Reference spectra of corannulene:Ar and H$_2$O:Ar at different concentrations were also recorded for comparison. All spectra were recorded at 10 K. Depositions were irradiated at 10 K with a mercury lamp ($\lambda$ > 235 nm, average power 150 mW, total fluence 0.72 J m$^{-2}$), and IR spectra measured after a series of irradiations (from 1 to 120 min).

# 3. Results

## 3.1. *Computational Results: Influence of the Matrix*

### 3.1.1. Structures and Energetics

Corannulene was inserted in different symmetry planes of the Ar matrix represented by a finite size fcc cluster initially of 1139 Ar atoms. After removing the minimum number of Ar atoms we found that the most stable situation corresponded to an insertion of the corannulene molecule in the [111] plane with a vacancy of 7 Ar atoms (see Figure 2, $E_{sub} = -2470$ cm$^{-1}$, also see Appendix A, Table A1 for the results of the other symmetry planes). We obtained the same result as for coronene [18] and other PAHs such as naphthalene [57]. The substitution energy is smaller than that computed for coronene using the same model ($E_{sub} = -4810$ cm$^{-1}$ [18]) which can be understood as the non planarity of corannulene leads to more perturbation of the Ar layers close to corannulene.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Optimised structure of C$_{20}$H$_{10}$:Ar obtained at the DFTB/FF level of theory. The total system is reported on the left hand side, the selection of the [111] plane where corannulene is inserted is represented on the right hand side.</figcaption></figure>

In Table 1 are reported the substitution energies obtained using Equation (1) for the most stable isomers of (C$_{20}$H$_{10}$)(H$_{2}$O)$_{n}$:Ar ($n$ = 1–3) of different types i.e., depending
whether the water molecules interact with the concave or convex faces of corannulene or
with its H atoms. The geometries of the most stable isomers are reported in Figures 3–5 for
*n* = 1, 2 and 3 respectively.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Most stable structures for (C<sub>20</sub>H<sub>10</sub>)(H<sub>2</sub>O):Ar obtained at the DFTB/FF level of theory (selection of layers).</figcaption></figure>

In the case of *n* = 1, the three conformers (1)¹side, (2)¹vex and (3)¹cav lie in an energy range of 150 cm⁻¹ (1.8 kJ·mol⁻¹), the most stable one corresponding to the water molecule interacting with an H atom of corannulene in the [111] symmetry plane ((1)¹side). Interestingly, the conformer in which the water molecule interacts with the concave face, known to be the most stable isomer in the gas phase [43], is the highest energy one ((3)¹cav), although the energy difference is small. These subtle differences from the gas phase originate from the constraints of the fcc lattice. For instance, in the (3)¹cav isomer, the water molecule substitutes an Ar atom whose position is not aligned with the C₅V axis of corannulene, and this prevents the water molecule from having the most stabilizing interaction with the corannulene as in the gas phase. When relaxing the geometry of this isomer after removing the matrix environment, the water molecule migrates towards the center of the “cup” and becomes the lowest energy isomer (see last column of Table 1).
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Most stable structures for (C₂₀H₁₀)(H₂O)₂:Ar obtained at the DFTB/FF level of theory (selection of layers).</figcaption></figure>

The perturbation of the [111] Ar layers adjacent to that containing corannulene due to the presence of the impurity can be visualized on the side views of the systems in Figures 3–5, a few Ar atoms being moved out of their original plane. Now comparing with the results obtained for coronene (C₂₄H₁₂, planar PAH) inside an Ar matrix using the same computational approach, the isomer similar to the (1)¹side isomer of the present work, designated as σ isomer in Simon et al. [18], was found more stable than the π isomer by ~800 cm⁻¹ in the case of coronene. This means that the stabilization of the “side” (or “σ”) isomer is enhanced in the case of a planar PAH with respect to the case of corannulene where it appears in competition with the two “π” isomers (2)¹vex and (3)¹cav.

In the case of *n* = 2, the most stable isomer corresponds to the water dimer interacting through two H atoms with the convex face of corannulene ((1)²vex, see Figure 4), the two H atoms pointing towards two central C atoms of corannulene. The corresponding isomer in which the water dimer interacts with the concave face of corannulene ((3)²conv) is found 210 cm⁻¹ (or 2.5 kJ·mol⁻¹) higher in energy. The isomer in which the water dimer interacts with both the concave (resp. convex) face and side H of corannulene is found 190 cm⁻¹ or 2.3 kJ·mol⁻¹ (resp. 240 cm⁻¹ or 2.9 kJ·mol⁻¹) higher in energy. As in the case of *n* = 1, when the geometries of (C₂₀H₁₀)(H₂O)₂ are further relaxed after removal of the Ar matrix, the water dimer tends to move towards the center of the cup when it interacts with the concave face of corannulene and it leads to the stabilisation of such conformers (see last column of Table 1).
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Most stable structures for (C<sub>20</sub>H<sub>10</sub>)(H<sub>2</sub>O)<sub>3</sub>:Ar obtained at the DFTB/FF level of theory (selection of layers).</figcaption></figure>

<table><caption>**Table 1.** Substitution energies of the (C<sub>20</sub>H<sub>10</sub>)(H<sub>2</sub>O)<sub>n</sub> complexes (n ≤ 3) and, for a given stoichiometry, relative energy with respect to the most stable system. In parenthesis are added the relative energy of the corresponding isolated (C<sub>20</sub>H<sub>10</sub>)(H<sub>2</sub>O)<sub>n</sub> aggregate (after removing the Ar matrix) after local DFTB geometry optimisation.</caption><thead><tr><th>n<sub>H<sub>2</sub>O</sub></th><th>Number of Atoms</th><th>Isomer</th><th>E<sub>sub</sub> (cm<sup>-1</sup>)</th><th>&Delta;E (kJ&middot;mol<sup>-1</sup>)</th></tr></thead><tbody><tr><td rowspan="3">1</td><td rowspan="3">1164</td><td>(1)<sup>1side</sup></td><td>-3590</td><td>0.0 (+10.3)</td></tr><tr><td>(2)<sup>1vex</sup></td><td>-3550</td><td>+0.5 (+7.5)</td></tr><tr><td>(3)<sup>1cav</sup></td><td>-3440</td><td>+1.8 (0.0)</td></tr><tr><td rowspan="4">2</td><td rowspan="4">1166</td><td>(1)<sup>2vex</sup></td><td>-5890</td><td>0.0 (+16.5)</td></tr><tr><td>(2)<sup>2cav-side</sup></td><td>-5700</td><td>+2.3 (0.0)</td></tr><tr><td>(3)<sup>2cav</sup></td><td>-5680</td><td>+2.5 (+3.4)</td></tr><tr><td>(4)<sup>2vex-side</sup></td><td>-5650</td><td>+12.9 (+9.5)</td></tr><tr><td rowspan="4">3</td><td rowspan="4">1168</td><td>(1)<sup>3cav</sup></td><td>-9200</td><td>0.0 (0.0)</td></tr><tr><td>(2)<sup>3vex</sup></td><td>-8860</td><td>+4.1 (+10.5)</td></tr><tr><td>(3)<sup>3vex-side</sup></td><td>-8610</td><td>+7.1 (+12.3)</td></tr><tr><td>(4)<sup>3cav-side</sup></td><td>-8420</td><td>+9.3 (+0.4)</td></tr></tbody></table>

In the case of n = 3, the isomer in which the water cyclic trimer, located in the [111] plane adjacent to that of corannulene, with its center approximately on the C<sub>5V</sub> axis of corannulene, interacts with the concave face of corannulene through three H atoms ((1)<sup>3cav</sup>, see Figure 5), is shown to be particularly stable with our model. Contrary to the smaller clusters (n = 1, 2), it is also the case after relaxation of the cluster after removing the Ar matrix as this conformer remains the most stable one in the gas phase. Next isomer, (2)<sup>3vex</sup>,
in which the cyclic trimer interacts with the convex face of corannulene through two H
atoms, is found 340 cm⁻¹ higher in energy. The next two isomers in which the water
trimer adopts a "linear form" interacting with both one face (convex or concave) and an H
atom of corannulene, namely the (3)³vex−side and (4)³cav−side are found respectively 590 and
780 cm⁻¹ (7.1 and 9.3 kJ·mol⁻¹) higher in energy than the stable (1)³cav isomer.

3.1.2. IR Spectra

This section reports the computed IR spectra of all the (C20H10)(H2O)n:Ar isomers for n = 1–3 and the comparison with those of C20H10 and of the corresponding (H2O)n clusters determined at the same level of theory. These simulations aim at getting insights into the influence of the coordination site (concave, convex or side) and of the geometry of the water cluster (linear or cycle)—in the case of the trimer—on the IR spectrum of the water molecule/cluster and corannulene upon coordination. The values of the positions and the shifts of the bands of interest upon coordination are reported in Tables A2–A4 for the harmonic spectra (n = 1 to 3 respectively), whereas Tables 2–4 report the data for finite-temperature (10 K) spectra. For the latter, the maxima were determined using the full width at half maximum (FWHM). Due to the fluctuations of the simulated IR cross sections and the significant broadening and merge of bands in some cases, we can estimate a human error of ±5 cm⁻¹ on the reported wavenumbers. The 10 K dynamic spectra of all the isomers reported in Figures 3–5 are displayed in Figures 6–8.

First of all, we found that, for all clusters, the bands of corannulene were hardly affected
upon coordination, and it is in line with the experimental results. This is clearly the case
of the v<sub>CH</sub> (C-H stretch) bands, and for the δ(CH) + v(CC) bands (combination of in plane
C-H bending and C-C stretches) that have very weak intensities in the computed spectra.
Regarding the γ<sub>CH</sub> (C-H out-of-plane) most modes, their energy can be slightly blue-shifted
for the "side" isomers all the more as the cluster size increases. The most intense γ<sub>CH</sub> band
is shifted by +5 cm<sup>-1</sup> for (1)<sup>1-side</sup>, +9 cm<sup>-1</sup> for (2)<sup>2cav-side</sup>, +8 and +14 cm<sup>-1</sup> for (3)<sup>3vex-side</sup>
and (4)<sup>3cav-side</sup> respectively using the harmonic data. This small blue shift was also observed
for the σ isomer the coronene-water complex [18]. The computed dynamic shifts were found
similar to the harmonic ones. Therefore a side isomer may be distinguished from a convex or
concave isomer through the position of the γ<sub>CH</sub> band as in the case of coronene [18]. It should
be noticed that these shift values are small with respect to the accuracy of the computed
absolute positions. However, the description of intermolecular interactions has been carefully
benchmarked against wavefunction results, insuring their accuracy [48].

In the following, we focus on the influence of the interaction with corannulene on the
IR spectra of (H₂O)ₙ:Ar (n = 1–3).

Regarding the dynamic spectra of the isomers with one water molecule (Figure 6), the
first striking feature is the narrowing of the water bands upon coordination. Indeed, in our
simulations of H₂O:Ar at 10 K, the water molecule vibrates and rotates in the matrix, leading
to the broadening of the resonant symmetric and assymmetric stretching modes (ν₁ and ν₃)
as well as the bending mode (ν₂). Besides, the stretching modes are significantly redshifted
with respect to the harmonic ones, as we reported in our previous work focusing on the IR
spectra of the water molecules and clusters inside a rare gas matrix [56]. As can be seen in
Tables A2 and 2, the shifts are quite variable depending on the isomer. For the (1)¹side isomer,
the ν₁ and ν₃ modes are blueshifted upon coordination, the blue shift being more pronounced in
the anharmonic spectrum (+38 and +35 vs. +4 and +2 cm⁻¹ for the anharmonic and harmonic
ν₁ and ν₃ modes respectively). On the opposite, in the case of the (2)¹vex isomer, the ν₁ and
ν₃ bands are redshifted, the shifts being less pronounced in the dynamic spectra (−36 and
−43 cm⁻¹ vs. −17 and −1 cm⁻¹). The case of (3)¹cav is intermediate as the ν₁ and ν₃ modes are
redshifted at the harmonic level and become slightly blueshifted in the 10 K dynamic spectra
(Figure 6). Regarding the ν₂ mode, it is blueshifted (except in the (1)¹side harmonic spectrum)
increasing from the side to convex and concave isomers. However, the significant narrowing of
the water bands upon coordination may lead to the impossibility to detect the (C₂₀H₁₀)(H₂O):Ar
clusters in the matrix as they can be hindered by the bands of the water molecule.
For larger clusters, only the dynamic spectra will be discussed. Regarding (C₂₀H₁₀)₀,₁ (H₂O)₂:Ar, let us just specify that, regarding the 10 K spectrum of the water dimer, the v₁ and v₃ modes are redshifted with respect to the harmonic spectra. The shifts obtained for the O-H modes for the H atoms interacting with the Ar atoms of the matrix (the v₁ and v₃ most energetic modes, usually designated as v₁a and v₃a, "a" standing for "acceptor") are the most shifted (both by 48 cm⁻¹) whereas the lowest energetic ones (v₁d and v₃d, "d" standing for "donor") are only shifted by 7 and 9 cm⁻¹ respectively. It can be understood as the former interact with the Ar atoms of the matrix and therefore are more sensitive to the loose interaction with Ar expected to lead to anharmonic effects at low temperature. Extremely weak intensity bands also accompany the intense ones, to a lesser extent than for the water monomer because the motion is more hindered by the matrix environment. The water dimer vibrates around its equilibrium position with the possible occurrence of the rotation of the H-acceptor water monomer around the OO axis. Interestingly, the water dimer motion in our simulations appears similar to the experimental evidence of the "acceptor switching and rotation of the water dimer around its O-O axis" by Ceponkus et al. [58].

<table><caption>**Table 2.** Positions (cm⁻¹) of the IR anharmonic bands in matrix of the corannulene, water monomer and the corresponding (C₂₀H₁₀)(H₂O) isomers. For the latter, the shifts of the water modes induced by the adsorption on corannulene are also reported (/). For H₂O, the numbers in parenthesis are the positions determined for the longer simulations (see Appendix B and Figure A2). Experimental data are mentioned for comparison. <sup>a</sup> this work, <sup>b</sup> [59] <sup>c</sup> [60].</caption><thead><tr><th></th><th>Mode</th><th>This Work</th><th>Expt.</th><th>(1)<sup>1side</sup></th><th>(2)<sup>1vex</sup></th><th>(3)<sup>1cav</sup></th></tr></thead><tbody><tr><td rowspan="5">C<sub>20</sub>H<sub>10</sub></td><td>γ(CH)</td><td>811</td><td>837<sup>a</sup></td><td>816</td><td>812</td><td>813</td></tr><tr><td>δ(CH)</td><td>1171</td><td>1139<sup>a</sup></td><td>1180</td><td>1175</td><td>1174</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1270</td><td>1316<sup>a</sup></td><td>1270</td><td>1269</td><td>1269</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1640</td><td>1441<sup>a</sup></td><td>1639</td><td>1641</td><td>1640</td></tr><tr><td>ν(CH)</td><td>1849</td><td>3050<sup>a</sup></td><td>1852</td><td>1849</td><td>1851</td></tr><tr><td rowspan="3">(H<sub>2</sub>O)</td><td>v<sub>1</sub></td><td>3756 (3760)</td><td>3638.5<sup>b</sup>, 3639.2<sup>c</sup></td><td>3794 /+38</td><td>3739 /−17</td><td>3763 /+7</td></tr><tr><td>v<sub>2</sub></td><td>1561 (1563)</td><td>1589<sup>b</sup>, 1589.1<sup>c</sup></td><td>1568 /+7</td><td>1582 /+21</td><td>1593 /+32</td></tr><tr><td>v<sub>3</sub></td><td>4021 (4021)</td><td>3733.5<sup>b</sup>, 3737.2<sup>c</sup></td><td>4056 /+35</td><td>4020 /−1</td><td>4028 /+7</td></tr></tbody></table>

<table><caption>**Table 3.** Positions (cm⁻¹) of the IR anharmonic bands in matrix of the corannulene, water dimer and the corresponding (C₂₀H₁₀)(H₂O)₂. For the latter, the shifts of the water modes induced by the adsorption on corannulene are also reported (/). Experimental data are mentioned for comparison. <sup>a</sup> this work, <sup>b</sup> [58].</caption><thead><tr><th></th><th>Mode</th><th>This Work</th><th>Expt.</th><th>(1)<sup>2vex</sup></th><th>(2)<sup>2cav-side</sup></th><th>(3)<sup>2cav</sup></th><th>(4)<sup>2vex-side</sup></th></tr></thead><tbody><tr><td rowspan="5">C<sub>20</sub>H<sub>10</sub></td><td>γ(CH)</td><td>811</td><td>837<sup>a</sup></td><td>813</td><td>820</td><td>814</td><td>818</td></tr><tr><td>δ(CH)</td><td>1176</td><td>1139<sup>a</sup></td><td>1176</td><td>1174</td><td>1175</td><td>1178</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1270</td><td>1316<sup>a</sup></td><td>1269</td><td>1267</td><td>1268</td><td>1261</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1640</td><td>1441<sup>a</sup></td><td>1641</td><td>1639</td><td>1640</td><td>1639</td></tr><tr><td>ν(CH)</td><td>1849</td><td>3050<sup>a</sup></td><td>1850</td><td>1852</td><td>1849</td><td>1852</td></tr><tr><td rowspan="5">(H<sub>2</sub>O)<sub>2</sub></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>v<sub>1</sub></td><td>3658</td><td>3574<sup>b</sup></td><td>3630 /−28</td><td>3593 /−65</td><td>3624 /−34</td><td>3627 /−31</td></tr><tr><td></td><td>3760</td><td>3633<sup>b</sup></td><td>3748 /−12</td><td>3745 /−15</td><td>3742 /−18</td><td>3754 /−6</td></tr><tr><td>v<sub>2</sub></td><td>1562</td><td>1597<sup>b</sup></td><td>1568 /+6</td><td>1570 /+8</td><td>1585 /+23</td><td>1569 /+7</td></tr><tr><td></td><td>1587</td><td>1610<sup>b</sup></td><td>1597 /+10</td><td>1586 /−1</td><td>1608 /+21</td><td>1580 /−7</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>v<sub>3</sub></td><td>3982</td><td>3708<sup>b</sup></td><td>3926 /−56</td><td>3956 /−26</td><td>3924 /−58</td><td>3982 /0</td></tr><tr><td></td><td></td><td>4024</td><td>3730<sup>b</sup></td><td>4031 /+7</td><td>4013 /−11</td><td>4008 /−16</td><td>4028 /+4</td></tr></tbody></table>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** IR spectra derived from MD/DFTB-FF simulations at 10 K of the three lowest energy isomers of (C₂₀H₁₀)(H₂O):Ar (red), H₂O:Ar (blue), and C₂₀H₁₀:Ar (green).</figcaption></figure>
<table><caption>**Table 4.** Positions (cm⁻¹) of the IR anharmonic bands in matrix of the corannulene, water trimer and the corresponding (C₂₀H₁₀)(H₂O)₃ isomers. In the case of the presence of broad bands, the maxima were determined using the FWHM, which is indicated in parenthesis. The shifts of the water modes induced by the adsorption on corannulene are also reported (/). * The shift was determined using a maximum position of 1579 cm⁻¹ for the v₂ band of (H₂O)₃ estimated considering the structured broad band with a FWHM of 30 cm⁻¹. Experimental data are mentioned for comparison. a this work, b [61].</caption><thead><tr><th>Mode</th><th>This Work</th><th>Expt.</th><th>(1)<sup>3cav</sup></th><th>(2)<sup>3vex</sup></th><th>(3)<sup>3vex-side</sup></th><th>(4)<sup>3cav-side</sup></th></tr></thead><tbody><tr><td rowspan="6">C<sub>20</sub>H<sub>10</sub></td><td>γ(CH)</td><td>811</td><td>837<sup>a</sup></td><td>815</td><td>813</td><td>819</td><td>825</td></tr><tr><td>δ(CH)</td><td>1176</td><td>1139<sup>a</sup></td><td>1175</td><td>1175</td><td>1180</td><td>1176</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1270</td><td>1316<sup>a</sup></td><td>1269</td><td>1273</td><td>1272</td><td>1271</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1640</td><td>1441<sup>a</sup></td><td>1639</td><td>1639</td><td>1639</td><td>1638</td></tr><tr><td>ν(CH)</td><td>1849</td><td>3050<sup>a</sup></td><td>1851</td><td>1849</td><td>1851</td><td>1852</td></tr><tr><td>ν<sub>1</sub></td><td>3596 (50)</td><td>3514<sup>b</sup></td><td>3591 (12) /-5</td><td>3605 (49) /+9</td><td>3554 (12) /-42<br>3595 (10) /-1<br>3721 (10) /+125</td><td>3537 (11) /-59<br>3584 (10) /-12<br>3741 (8) /+145</td></tr><tr><td rowspan="3">(H<sub>2</sub>O)<sub>3</sub></td><td rowspan="3">ν<sub>2</sub></td><td>1566</td><td>1566 (0)</td><td>1548 /-18</td><td rowspan="3">1586 (20) /+7 *</td><td rowspan="3">1577 /+11<br>1592 /+22<br>1594 /+3</td></tr><tr><td>1570</td><td>1571 /+1</td><td>1581 /+11</td></tr><tr><td>1591</td><td>1599 /+8</td><td>1600 /+9</td></tr><tr><td rowspan="3">ν<sub>3</sub></td><td rowspan="3">3948 (26)</td><td rowspan="3">3699<sup>b</sup></td><td rowspan="3">3902 (13) /-46</td><td rowspan="3">3950 (76) /+3</td><td>3944 /-4</td><td>3940 /-8</td></tr><tr><td>3960 /+12</td><td>3954 /+6</td></tr><tr><td>4012 /+64</td><td>4001 /+53</td></tr></tbody></table>

<figure><img src="image_11.png" /><figcaption>(1)²vex</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(2)²cav-side</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>(3)²cav</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>(4)²vex-side</figcaption></figure>

<caption>**Figure 7.** IR spectra derived from MD/DFTB-FF simulations at 10 K of the four lowest energy isomers of (C₂₀H₁₀)(H₂O)₂:Ar (red), (H₂O)₂:Ar (blue), and C₂₀H₁₀:Ar (green).</caption>
<figure><img src="image_15.png" /><figcaption>**Figure 8.** IR spectra derived from MD/DFTB-FF simulations at 10 K of the four lowest energy isomers of (C₂₀H₁₀)(H₂O)₃:Ar (red), (H₂O)₃:Ar (blue), and C₂₀H₁₀:Ar (green).</figcaption></figure>

Regarding the spectra of all C₂₀H₁₀)(H₂O)₂:Ar isomers reported Figure 7, the differences in the shifts of the bands with respect to those of the bare water dimer might help identifying the type of isomer (cav, vex, cav-side or vex-side) formed in the experiments. Upon interaction with corannulene, the (H₂O)₂ v₁ modes are globally redshifted and the shifts values are smaller in the 10 K spectra than in the harmonic ones. The shifts values are dependent on the isomer being smaller for the (1)²vex and (4)²vex-side isomers (see Table 3). The v₃ bands are redshifted for (2)²cav-side and (3)²cav, but is it not the case for the two other isomers for which the higher energy v₃ band is blueshifted. Finally, a difference between the shifts of the v₂ bands between the different isomers is noticeable. For instance, those for (1)²vex and (3)²cav are both blueshifted, the blue shift being larger for the latter.

Finally, the dynamic spectra of (C₂₀H₁₀)₀,₁(H₂O)₀,₃:Ar are reported in Figure 8. Regarding the water modes of (H₂O)₃, each type of mode leads to three broad bands more or less structured centered at 3596, 1569 and 3948 cm⁻¹ for v₁ to v₃. During the simulation, the water trimer vibrates around its equilibrium position with occasional flipping of the O-H bonds- not involved in intermolecular O-H bonds- from one side of the plane of the OOO cycle to the other side. As in the case of one and two water molecules, the water IR bands become narrower due to a hindrance of their motion upon interaction with corannulene, although it is not as relevant in the case of (2)³vex where the water trimer globally maintains its structure. It is all the more verified for the “cav” isomers (1)³cav and (4)³cav-side, reflecting
stronger corannulene-water interactions than for the “vex” isomers. A significant difference is noticeable between the “side” isomers and the others. In the former, a band can be found at ~3720 and 3740 cm⁻¹, that was assigned to the O-H stretching mode for the bond interacting with the PAH [51]. In contrast with the other stoichiometries, the stretching modes do not appear systematically reshifted. The $v_1$ mode is hardly redshifted in the case of the most stable (1)³cav isomers whereas it is slightly blueshifted in the (2)³vex spectrum. The $v_3$ mode is redshifted for (1)³cav but leads to mostly blueshifted bands for the three other isomers. The $v_2$ modes are also globally blueshifted for all isomers, the amplitude of the shifts being dependent on the isomer. As in the case of the other stoichiometries, the differences in the spectra may allow us to distinguish specific isomers formed in the experiment.

## 3.2. Experiments

In this section, we present the FTIR spectra of water, corannulene, and water:corannulene mixtures recorded in argon matrices at 10 K, where water is present in the form of monomers, dimers, trimers and larger aggregates, depending on the H₂O:Ar ratio. The data will be compared to the computational results presented in the above section and discussed in the next section.

The IR spectra of water in solid argon have been previously studied in detail by several authors (see, for example, refs. [59–67]). Our experimental spectra agree with literature spectra and our assignments are based on them. The assignments of water monomers and water aggregates formed in our deposits will allow us to identify perturbed bands (shifts) and/or additional bands induced by the presence of corannulene. Figure A3 (a reproduction of Figure 3 in ref. [18]) presents the FTIR spectra of different H₂O:Ar deposits at 10 K with increasing concentrations of water: (a) H₂O:Ar = 1:500, where only monomers and dimers are present; (b) H₂O:Ar = 1:50, where larger aggregates (from trimers to hexamers) appear; and (c) H₂O:Ar = 1:25, where a mixture of water aggregates and water ice (ASW) is observed. These spectra have been discussed in detail in our previous paper [18] and Table A5 lists the main band positions of water aggregates (H₂O)ₙ, with n = 1–6, as well as those of some water complexes H₂O:CO₂, H₂O:N₂ and (H₂O)₂:N₂.

Although corannulene and coronene are both aromatic, they present marked differences in structure and reactivity [20]. Due to its C5v symmetry, and its bowl-like non-planar structure, the corannulene molecule exhibits a richer IR spectrum than coronene, with more active modes and larger intensities. Table A6 lists the vibrational modes observed and calculated in this work, compared to experimental data from the NASA Ames PAH IR spectral database [68–70] and previous study by Rouillé et al. [39]. A summary of the main experimental bands is reported in Tables 2–4 where the frequencies are compared to the IR anharmonic band positions calculated in this work.

Figure 9 presents the FTIR spectra of corannulene in Ar and in H₂O:Ar (1:50) in two spectral regions of interest: 2000–400 cm⁻¹, where we can observe the main corannulene bands plus the deformation modes $v_2$ of water; and 4000–2800 cm⁻¹, where we can observe all OH stretching modes of water monomers and aggregates plus the CH stretching bands of corannulene. Despite several attempts, it was impossible to obtain a pure corannulene:Ar spectrum without any traces of water; a small amount of water monomer (lines at 3756 and 3778 cm⁻¹) and trace of (H₂O)₂ (line at 3708 cm⁻¹) are always present in the experimental corannulene:Ar spectra. No change is observed in water monomer lines, nor in corannulene lines in this spectrum, so we consider that the trace amount of water in the matrix is not interacting with the corannulene.

If we now consider depositions with intentionally added water H₂O:Ar 1:50 (Figure 9, red traces labelled c), we can see that, apart from a minor broadening, the vibrational bands of corannulene remain almost unchanged in frequency compared to corannulene in argon. Contrary to what was observed in the case of coronene [17], where the most intense out-of-plane $\gamma$(CH) mode was blue-shifted by up to 12 cm⁻¹ upon interaction with water aggregates, no or very small blue shifts (a few cm⁻¹ i.e., within the spectral resolution) are observed for the $\gamma$(CH) mode of corannulene at 837 cm⁻¹. The main changes in the FTIR
spectrum appear in the OH stretching modes of water, especially in the bands assigned to water dimers and trimers. Small blue shifts and broadenings are observed for the bands located at 3700 (trimer), 3706 (dimer) and 3712 (monomer) cm⁻¹. No change is detected below 3500 cm⁻¹, in the region of the stretching modes of (H₂O)ₙ with n = 4–6.

<figure><img src="image_16.png" /><figcaption>**Figure 9.** Experimental spectra measured in the wavelength ranges (a) 4000–2800 cm⁻¹ and (b) 2000–400 cm⁻¹. Spectra in both panels are as follows: a—Corannulene in argon (green trace), b—water in argon (blue trace), and c—corannulene and water in argon (red trace).</figcaption></figure>

Figure 10 presents the difference spectrum between the corannulene:H₂O:Ar deposition minus the H₂O:Ar deposition at 10 K in the 3800–3200 cm⁻¹ and 1800–1500 cm⁻¹ ranges, where stretching and bending water modes are, respectively, located. A new band appears at 3727 cm⁻¹, as well as two broad, weak bands at 3603 and 3613 cm⁻¹. These features are in proximity to a weak band in the H₂O:Ar spectrum assigned to a symmetric stretching mode of the water trimer at 3611 cm⁻¹. A weak band is also observed at 3536 cm⁻¹, in the region of the symmetric stretching modes of water dimers (3574 cm⁻¹) and trimers (3518 cm⁻¹). Additionally, a new band clearly appears at 1598 cm⁻¹, in the HOH bending region of water dimers at 1593 cm⁻¹ and water trimers at 1602 cm⁻¹. We can thus conclude that molecular interaction between water dimers and/or trimers with corannulene are observed in the FTIR spectra.

Energy was introduced into the corannulene:H₂O:Ar samples via photons. Low energy UV irradiations of our deposition did not induce any change in the spectrum, contrary to what was previously observed for coronene:H₂O:Ar matrices at 10 K, where the interaction of coronene with water aggregates (H₂O)ₙ with n up to 3 led to the formation of oxygenated products [17].
<figure><img src="image_17.png" /><figcaption>**Figure 10.** Experimental difference spectra in the wavelength ranges (a) 3800–3200 cm⁻¹ and (b) 1800–1500 cm⁻¹. Spectra b (water, blue) and c (corannulene plus water, red) of Figure 9 are reproduced in dot-dashed lines, while the solid black line represents their difference spectrum (i.e., spectrum c minus spectrum b). New bands discussed in the text are indicated with red arrows at the following frequencies: 3727, 3613, 3603, 3536, and 1598 cm⁻¹.</figcaption></figure>

**4. Discussion**

In this section, we discuss the assignment of the experimental spectra in the light of the calculated finite-temperature spectra (Tables 2–4 for corannulene in interaction with, respectively, water monomers, dimers and trimers). As the absolute positions of the calculated frequencies are not strictly accurate, we will use the calculated shifts to identify the (C₁₀H₂₀)(H₂O)ₙ complexes formed in our experiments. Indeed, regarding the computed absolute positions, we are aware that they are not accurate. It is due to several reasons: (i) the DFTB method is a DFT-based approximate method (truncated hamiltonian, minimal valence basis set...). However, the DFTB/FF method used in the present work, quite unique to our knowledge, is mandatory to be able to perform on-the-fly BOMD simulations taking explicitely into account all atoms (including the atoms of the rare gas matrix). The possible presence of defaults is not taken into account. (ii) 0 K vibrational anharmonicities and nuclear quantum effects, expected to be significant in the case of water clusters, are not taken into account, and it is a major drawback of our computational studies. However, given the size of the systems and the presence of highly symmetric corannulene, an approach such as vibrational perturbational theory using a wavefunction method to describe the electronic structure in conjunction with a large basis set, which would allow to obtain accurate 0 K IR spectra for such systems, is prohibited. In our approach however, as specified in Section 3.1.2, intermolecular interactions are expected to be better described than intramolecular bonds as they have been carefully benchmarked against wavefunction results [48].

The new band observed at 3727 cm⁻¹ may be due to a blue shift of 15 cm⁻¹ of the ν₃ mode of the water monomer (3712 cm⁻¹) in (C₁₀H₂₀)(H₂O) "side" or "cav" geometries
(see Table 2), or to a blue shift of 21 cm⁻¹ of the ν₃ mode of the water dimer (3706 cm⁻¹) in (C₁₀H₂₀)(H₂O)₂vex (see Table 3). These two values being located between those calculated for these two complexes (7 to 35 cm⁻¹), it is impossible to conclude at this stage. However, it seems unlikely that (C₁₀H₂₀)(H₂O) with "side" structure has been formed, as no large red shift has been observed on the ν(CH) of corannulene, whereas (C₁₀H₂₀)(H₂O) in the "cav" form, calculated to induce only a small blue shift in the ν(CH) modes, could be a possible candidate (see Table 2).

If we now consider the new bands observed at 3603 and 3613 cm⁻¹, they may be redshifted from the ν₃ modes of the water dimer (3706 cm⁻¹, shift of 93 cm⁻¹) or trimer (3700 cm⁻¹, shift of 87 cm⁻¹), or blueshifted from the ν₁ mode of the water dimer (3574 cm⁻¹, shifts of 29 and 39 cm⁻¹). From Tables 3 and 4, it appears that no blue shifts are calculated for the ν₁ modes of any (C₁₀H₂₀)(H₂O)₂ structures nor for (C₁₀H₂₀)(H₂O)₃cav, whereas large red shifts are expected for the ν₃ modes of these complexes (from 46 to 58 cm⁻¹). So, the more likely assignment for these two new bands is to (C₁₀H₂₀)(H₂O)₂ "vex" or "cav" forms and to (C₁₀H₂₀)(H₂O)₃cav. For the same reasons, the weak band observed at 3536 cm⁻¹ must be red shifted from the ν₁ mode of the water dimer (3574 cm⁻¹, shift 38 cm⁻¹) in the (C₁₀H₂₀)(H₂O)₂ "vex" or "cav" forms, as only a small blue shift (5 cm⁻¹, see Table 4) is expected from the ν₁ mode of the water trimer observed at 3518 cm⁻¹.

Finally, the new band at 1598 cm⁻¹ could be due to a blue shift of 6 or 5 cm⁻¹ of the ν₂ mode of, respectively, the water monomer observed at 1592 cm⁻¹ or of the dimer at 1593 cm⁻¹, or could be due to a red shift of 4 cm⁻¹ of the ν₂ mode of the water trimer at 1602 cm⁻¹. Comparison with calculated shifts in Tables 3 and 4 shows that only blue shifts are expected for all ν₂ modes in all complexes, whatever the stoichiometry and structures. However, this blue shift is calculated to be 7 cm⁻¹ in (C₁₀H₂₀)(H₂O)side, but more than 30 cm⁻¹ in (C₁₀H₂₀)(H₂O)cav. As we have already eliminated (C₁₀H₂₀)(H₂O)side due to the absence of shifts in the ν(CH) bands of corannulene, we thus assign this new band to (C₁₀H₂₀)(H₂O)₂vex.

From the analysis of our experimental spectra—which is complicated by the presence of excess corannulene and water aggregates (n = 1 to 6)—and thanks to our calculations, we can conclude that: (1) no (C₁₀H₂₀)(H₂O) complex has been detected in our experimental conditions. This can be due to the effect of argon environment, or to a better reactivity of corannulene with small aggregates (dimer, trimer) compared to the water monomer, (2) the (C₁₀H₂₀)(H₂O)₃ complex with the water trimer interacting with the concave face of corannulene has been identified, which is in line with the fact that it is the most stable calculated structure, and (3) the (C₁₀H₂₀)(H₂O)₂ complex is also formed, but it is impossible to firmly characterize its structure, which may be "vex" or "cav". We have not seen any experimental evidence for the formation of complexes with larger aggregates (n = 4 to 6) on the difference spectra. Even if their formation cannot be completely excluded, it is hardly probable that they could be formed, because the water tetramer and larger cluster spectra are quite different from those of water dimer and trimer and located lower in energy.

## 5. Conclusions

Our study has shown that the weak PAH-water complexes formed in rare gas matrices may be different from those obtained in the gas phase. Our prior computations have shown that embedding in the matrix stabilises σ isomers in the case of planar aromatics [18], while for the clusters involving the non-planar corannulene studied in this work, π isomers remain favored in the cases of the water dimer and trimer. Such geometrical differences may account for the differences in photo-reactivity of these systems: indeed, pyrene and coronene react with water molecules to form alcohols and quinones under low energy irradiation [13,14] while corannulene does not, even when embedded in ASW (private comm. from J. Mascetti). The photoreactivity in planar PAHs could occur through the population of low energy charge transfer PAH⁺-H₂O⁻ electronic excited states that are present when the water molecules interact through their oxygen atoms with the H atoms of the PAH [17,19]. Therefore, from the modeling reported in the present work, the most stable form of (C₂₀H₁₀)(H₂O)₃ in an argon matrix does not appear favorable to the presence
of such low energy excited states. The investigation of the electronic excited states of (C20H10)(H2O)n:Ar systems (n = 1–3) will be the object of future work. Our studies illustrate the dependence of the mutual orientation of the reactants on the nature of bimolecular photo-redox reactions as evoked in a recent review [71]. The differentiation between flat and curved PAHs may be of importance in UV absorption and electron transfer reactions in PAH:water molecular systems. From an astrophysical point of view, even if it is still premature to project these results directly to PAHs of astro-relevant sizes embedded in water ice, there is probably no difference in the photoreactivity between flat and curved PAHs when using high energy UV photons that lead to PAH ionization, but the influence of geometry should be taken into account when using low energy UV photons that induce charge transfer reactions.

**Author Contributions:** All authors contributed to the work. H.L. and A.S. performed calculations; C.A. and J.M. ran experiments. All authors have discussed the results. Draft preparation and writing were done by H.L., A.S., J.A.N. and J.M. All authors have read and agreed to the submitted version of the manuscript.

**Funding:** This research was funded by the French National Agency for Research (ANR PARCS project ANR-13-BS008-0005) and supported by the Programme National Physique et Chimie du Milieu Insterstellaire (PCMI) of CNRS/INSU with INC/INP co-funded by CEA and CNES, and by the French research network EMIE (Edifices Moléculaires Isolés et Environnés, GDR 3533).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data published in the present paper are made public in the archived datasets on the Zenodo website (https://doi.org/10.5281/zenodo.6368179) published on 18 March 2022.

**Acknowledgments:** H.L. and A.S. thank the computing mesocenter CALMIP (UMS CNRS 3667) for generous allocation of computer resources (p17002). J.M. and C.A. acknowledge the participation of internship students Maxime Gerbeaud-Lassau and Francisco Salvador in the experimental part of this study.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Abbreviations

<table><tr><td>AIBs</td><td>aromatic infrared bands</td></tr><tr><td>ASW</td><td>amorphous solid water</td></tr><tr><td>BOMD</td><td>Born-Oppenheimer molecular dynamics</td></tr><tr><td>DFTB</td><td>density functional based tight binding</td></tr><tr><td>FF</td><td>force field</td></tr><tr><td>IR</td><td>infrared</td></tr><tr><td>FTIR</td><td>Fourier transform infrared</td></tr><tr><td>ISM</td><td>interstellar medium</td></tr><tr><td>PAH</td><td>polycyclic aromatic hydrocarbon</td></tr><tr><td>UV</td><td>ultraviolet</td></tr></table>

## Appendix A. Energetics: Insertion of Corannulene in the Matrix

Different insertion planes for corannulene inside the Ar matrix were tried. The substitution energies are reported in Table A1.
<table><caption>**Table A1.** Substitution energies for the insertion of the corannulene in different cristallographic planes.</caption><thead><tr><th>Plane</th><th>Number of Ar Removed</th><th>E<sub>sub</sub> (cm<sup>-1</sup>)</th></tr></thead><tbody><tr><td>(100)</td><td>9</td><td>1400</td></tr><tr><td>(110)</td><td>5</td><td>3910</td></tr><tr><td>(111)</td><td>7</td><td>-2470</td></tr><tr><td>Cavity</td><td>13</td><td>3080</td></tr></tbody></table>

## Appendix B. Convergence of IR Dynamic Spectra

<figure><img src="image_18.png" /><figcaption>(a) $v_1$ mode for $(3)^{1cav}$</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>(b) $v_3$ mode for $(H_2O)_3$</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>(c) $v_3$ mode for $(1)^{3cav}$</figcaption></figure>

**Figure A1.** IR dynamic spectra at 10 K obtained following the procedure described in Section 2.1 for $(3)^{1cav}$ (a), $(H_2O)_3$ (b) and $(1)^{3cav}$ (c). Zooms on the $v_1$ (a) and $v_3$ (b,c) modes are presented as an illustration of the convergence. In the upper panel of each plot we present the spectra averaged over 7 × 100 ps (red), 5 × 100 ps (blue) and 3 × 100 ps (green). In the lower panel of each plot we show the difference spectra between 7 × 100 ps and 5 × 100 ps dynamics (red) and between 7 × 100 ps and 3 × 100 ps MDs (blue). The resolution of all spectra is 0.4 cm⁻¹.
<figure><img src="image_21.png" /><figcaption>**Figure A2.** IR dynamic spectra at 10 K for H₂O:Ar (a) and (H₂O)₂:Ar (b). The spectra of Figures 6 and 7, resulting from 7 × 100 ps MDs, are reported in blue. Those resulting from longer simulations, 7 × 500 ps for H₂O:Ar and 3 × 500 ps for (H₂O)₂:Ar, are reported in red.</figcaption></figure>

## Appendix C. Harmonic Data

<table><caption>**Table A2.** Positions (cm⁻¹) of the IR harmonic bands in matrix of the corannulene, water monomer and the corresponding (C₂₀H₁₀)(H₂O)₁ isomers. Harmonic IR intensities are reported in parenthesis, as well as shift values (//). Experimental data are also reported for comparison. a this work, b [59], c [60].</caption><thead><tr><th></th><th>Mode</th><th>This Work</th><th>Expt.</th><th>(1)<sup>1side</sup></th><th>(2)<sup>1vex</sup></th><th>(3)<sup>1cav</sup></th></tr></thead><tbody><tr><td rowspan="6">C₂₀H₁₀</td><td>γ(CH)</td><td>812 (147)</td><td>837<sup>a</sup></td><td>817 (135)</td><td>814 (149)</td><td>813 (153)</td></tr><tr><td>δ(CH)</td><td>1175 (2)</td><td>1139<sup>a</sup></td><td>1178 (3)</td><td>1177 (3)</td><td>1176 (3)</td></tr><tr><td></td><td>1271 (&lt;1)</td><td></td><td>1271(&lt;1)</td><td>1271(&lt;1)</td><td>1271 (&lt;1)</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1646 (6)</td><td>1316<sup>a</sup></td><td>1644 (6)</td><td>1646 (3)</td><td>1645 (4)</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1857 (3)</td><td>1441<sup>a</sup></td><td>1856 (3)</td><td>1855 (3)</td><td>1856 (2)</td></tr><tr><td>ν(CH)</td><td>2999* (60)</td><td>3050<sup>a</sup></td><td>2999 (62)</td><td>3002 (55)</td><td>3007 (45)</td></tr><tr><td rowspan="3">(H₂O)</td><td>ν₁</td><td>3807 (39)</td><td>3638.5<sup>b</sup>, 3639.2<sup>c</sup></td><td>3811 (138) / +4</td><td>3771 (52) / −36</td><td>3765 (57) / −42</td></tr><tr><td>ν₂</td><td>1571 (140)</td><td>1589<sup>b</sup>, 1589.1<sup>c</sup></td><td>1568 (289) / −3</td><td>1588 (241) / +17</td><td>1594 (178) / +23</td></tr><tr><td>ν₃</td><td>4081 (33)</td><td>3733.5<sup>b</sup>, 3737.2<sup>c</sup></td><td>4083 (83) / +2</td><td>4038 (21) / −43</td><td>4030 (10) / −51</td></tr></tbody></table>
<table><caption>**Table A3.** Positions (cm⁻¹) of the IR harmonic bands in matrix of the corannulene, water dimer and the corresponding (C₂₀H₁₀)(H₂O)₂ isomers. The shifts of the water modes upon interaction with corannulene are also reported (/). Experimental data are mentioned for comparison. <sup>a</sup> this work, <sup>b</sup> [58].</caption><thead><tr><th></th><th>Mode</th><th>This work</th><th>Expt.</th><th>(1)<sup>2vex</sup></th><th>(2)<sup>2cav-side</sup></th><th>(3)<sup>2cav</sup></th><th>(4)<sup>2vex-side</sup></th></tr></thead><tbody><tr><td rowspan="5">C<sub>20</sub>H<sub>10</sub></td><td>γ(CH)</td><td>812</td><td>837<sup>a</sup></td><td>815</td><td>821</td><td>816</td><td>819</td></tr><tr><td>δ(CH)</td><td>1175</td><td>1139<sup>a</sup></td><td>1174</td><td>1173</td><td>1175</td><td>1179</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1271</td><td>1316<sup>a</sup></td><td>1271</td><td>1271</td><td>1270</td><td>1271</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1646</td><td>1441<sup>a</sup></td><td>1646</td><td>1644</td><td>1645</td><td>1644</td></tr><tr><td>ν(CH)</td><td>1857</td><td>3050<sup>a</sup></td><td>1853</td><td>1855</td><td>1854</td><td>1856</td></tr><tr><td rowspan="5">(H<sub>2</sub>O)<sub>2</sub></td><td>ν<sub>1</sub></td><td>3665</td><td>3574<sup>b</sup></td><td>3628 / -37</td><td>3585 / -80</td><td>3639 / -26</td><td>3613 / -52</td></tr><tr><td>ν<sub>2</sub></td><td>3802</td><td>3633<sup>b</sup></td><td>3745 / -57</td><td>3757 / -45</td><td>3757 / -45</td><td>3752 / -50</td></tr><tr><td>ν<sub>3</sub></td><td>1568</td><td>1597<sup>b</sup></td><td>1566 / -2</td><td>1573 / +5</td><td>1586 / +18</td><td>1570 / +2</td></tr><tr><td>ν<sub>1</sub></td><td>1593</td><td>1610<sup>b</sup></td><td>1597 / +4</td><td>1593 (0)</td><td>1609 / +16</td><td>1581 / -12</td></tr><tr><td>ν<sub>3</sub></td><td>3991</td><td>3708<sup>b</sup></td><td>3925 / -66</td><td>3962 / -29</td><td>3941 / -50</td><td>3983 / -8</td></tr><tr><td></td><td></td><td>4072</td><td>3730<sup>b</sup></td><td>4029 / -43</td><td>4018 / -54</td><td>4021 / -51</td><td>4028 / -44</td></tr></tbody></table>

<table><caption>**Table A4.** Positions (cm⁻¹) of the IR harmonic bands in matrix of the corannulene, water trimer and the corresponding (C₂₀H₁₀)(H₂O)₃. The shifts of the water modes upon interaction with corannulene are also reported (/). Experimental data are mentioned for comparison. <sup>a</sup> this work, <sup>b</sup> [61].</caption><thead><tr><th></th><th>Mode</th><th>This Work</th><th>Expt.</th><th>(1)<sup>3cav</sup></th><th>(2)<sup>3vex</sup></th><th>(3)<sup>3vex-side</sup></th><th>(4)<sup>3cav-side</sup></th></tr></thead><tbody><tr><td rowspan="5">C<sub>20</sub>H<sub>10</sub></td><td>γ(CH)</td><td>812</td><td>837<sup>a</sup></td><td>816</td><td>814</td><td>820</td><td>826</td></tr><tr><td>δ(CH)</td><td>1175</td><td>1139<sup>a</sup></td><td>1174</td><td>1174</td><td>1178</td><td>1177</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1271</td><td>1316<sup>a</sup></td><td>1269</td><td>1271</td><td>1271</td><td>1273</td></tr><tr><td>δ(CH) + ν(CC)</td><td>1646</td><td>1441<sup>a</sup></td><td>1645</td><td>1645</td><td>1644</td><td>1642</td></tr><tr><td>ν(CH)</td><td>1857</td><td>3050<sup>a</sup></td><td>1854</td><td>1853</td><td>1855</td><td>1856</td></tr><tr><td rowspan="5">(H<sub>2</sub>O)<sub>3</sub></td><td>ν<sub>1</sub></td><td>3540</td><td>3514<sup>b</sup></td><td>3514 / -26</td><td>3518 / -22</td><td>3563 / +23</td><td>3536 / -4</td></tr><tr><td>ν<sub>2</sub></td><td>3597</td><td>3514<sup>b</sup></td><td>3580 / -17</td><td>3576 / -21</td><td>3602 / +5</td><td>3595 / -2</td></tr><tr><td>ν<sub>3</sub></td><td>3621</td><td>3514<sup>b</sup></td><td>3581 / -40</td><td>3599 / -22</td><td>3737 / +116</td><td>3748 / +127</td></tr><tr><td>ν<sub>1</sub></td><td>1568</td><td>1564 / -4</td><td>1546 / -22</td><td>1579 / +11</td><td>1577 / +9</td><td>1578 / +6</td></tr><tr><td>ν<sub>2</sub></td><td>1575</td><td>1601<sup>b</sup></td><td>1569 / -6</td><td>1580 / +5</td><td>1587 / +12</td><td>1581 / +6</td></tr><tr><td></td><td>ν<sub>3</sub></td><td>1595</td><td>1601<sup>b</sup></td><td>1600 / +5</td><td>1602 / +7</td><td>1592 / -3</td><td>1595 (0)</td></tr><tr><td></td><td>ν<sub>1</sub></td><td>3960</td><td>3900 / -60</td><td>3894 / -66</td><td>3959 / -1</td><td>3951 / -9</td><td>3951 / +9</td></tr><tr><td></td><td>ν<sub>2</sub></td><td>3964</td><td>3699<sup>b</sup></td><td>3904 / -60</td><td>3925 / -39</td><td>3972 / +8</td><td>3965 / +1</td></tr><tr><td></td><td>ν<sub>3</sub></td><td>3971</td><td>3699<sup>b</sup></td><td>3913 / -58</td><td>3959 / -12</td><td>4017 / +46</td><td>4008 / +37</td></tr></tbody></table>

## Appendix D. Ftir Bands of Water Monomers, Dimers, Trimers and Larger Aggregates in Argon Matrices at 10 K

<table><caption>**Table A5.** From ref. [18] Band positions, in cm⁻¹, of (H₂O)ₙ water species with n = 1 to 6. Attributions are based on previous works [59–62,64–67].</caption><thead><tr><th>n</th><th>Rovibrational Mode</th><th>Our Results (cm<sup>-1</sup>)</th><th>Exp Values (cm<sup>-1</sup>) from Litt.</th><th>Refs.</th></tr></thead><tbody><tr><td rowspan="5">1</td><td>ν<sub>2</sub> + ν<sub>3</sub></td><td>5346.1, 5325.5<br>5293.2, 5280.5</td><td>5345.9, 5325.4<br>5295.4, 5280.6</td><td></td></tr><tr><td>ν<sub>3</sub></td><td>3776.9, 3756.5<br>3724.8, 3711.4</td><td>3776.4, 3756.6<br>3724.9, 3711.3</td><td></td></tr><tr><td>ν<sub>1</sub></td><td>3669.9, 3653.6</td><td>3669.7, 3653.5</td><td rowspan="3">[59,60,62]</td></tr><tr><td>ν<sub>2</sub> RTC</td><td>1662.5, 1657.7</td><td>1661.4, 1657.2</td></tr><tr><td>ν<sub>2</sub></td><td>1636.3, 1624.2<br>1608.1, 1573.3<br>1556.9</td><td>1636.5, 1623.8<br>1607.9, 1573.1<br>1556.7</td></tr><tr><td>ν<sub>2</sub> NRM</td><td>1590.7</td><td>1589.2</td><td></td></tr></tbody></table>
<table><caption>Table A5. *Cont.*</caption><thead><tr><th>n</th><th>Rovibrational Mode</th><th>Our Results (cm<sup>-1</sup>)</th><th>Exp Values (cm<sup>-1</sup>) from Litt.</th><th>Refs.</th></tr></thead><tbody><tr><td rowspan="6">2</td><td>v<sub>3</sub> PD</td><td>3707.9</td><td>3708</td><td></td></tr><tr><td>v<sub>3</sub> PA</td><td>3736.1</td><td>3737.8</td><td></td></tr><tr><td>v<sub>1</sub> PA</td><td>3633.1</td><td>3633.1</td><td></td></tr><tr><td>v<sub>1</sub> PD</td><td>3573.8</td><td>3574</td><td></td></tr><tr><td>v<sub>2</sub> PD</td><td>1611.9</td><td>1610.6</td><td></td></tr><tr><td>v<sub>2</sub> PA</td><td>1593.2</td><td>1593.1</td><td></td></tr><tr><td rowspan="5">3</td><td>v<sub>3</sub> PD</td><td>3699.8</td><td>3700</td><td></td></tr><tr><td>v<sub>1</sub> PA</td><td>3611.1</td><td>3612</td><td></td></tr><tr><td>v<sub>1</sub> PD</td><td>3527.9</td><td>3528</td><td></td></tr><tr><td>v<sub>1</sub> PD</td><td>3517.4</td><td>3516</td><td>[59,61,64,65]</td></tr><tr><td>v<sub>2</sub> PA</td><td>1601.9</td><td>1602</td><td></td></tr><tr><td>4</td><td></td><td>3372.9</td><td>3374</td><td></td></tr><tr><td>5</td><td></td><td>3326.7</td><td>3327</td><td></td></tr><tr><td rowspan="6">6</td><td></td><td>3541.3</td><td>3540</td><td></td></tr><tr><td></td><td>3444.4</td><td>3445</td><td></td></tr><tr><td></td><td>3392.2</td><td>3391</td><td>[64,67]</td></tr><tr><td></td><td>3210.3</td><td>3212</td><td></td></tr><tr><td></td><td>3186.2</td><td>3195.2</td><td></td></tr><tr><td></td><td>3151.5</td><td>3150</td><td>[62]</td></tr><tr><td rowspan="2">H<sub>2</sub>O-CO<sub>2</sub></td><td>v<sub>1</sub></td><td>3638.4</td><td>3638</td><td rowspan="2">[60]</td></tr><tr><td>v<sub>2</sub></td><td>1589.5</td><td>1589.5</td></tr><tr><td rowspan="2">H<sub>2</sub>O-N<sub>2</sub></td><td>v<sub>3</sub></td><td>3730.3</td><td>3729.6</td><td rowspan="4">[66]</td></tr><tr><td>v<sub>1</sub></td><td>3640.2</td><td>3640.2</td></tr><tr><td rowspan="2">(H<sub>2</sub>O)<sub>2</sub>-N<sub>2</sub></td><td>v<sub>1</sub> PD</td><td>3566.8</td><td>3566.6</td></tr><tr><td>v<sub>1</sub> PD</td><td>3563.1</td><td>3563.5</td></tr></tbody></table>

<figure><img src="image_22.png" /><figcaption>**Figure A3.** From ref. [18] Experimental FTIR spectra of H₂O:Ar depositions at 10 K (a) 1:500 (black, only water monomers and dimers are present); (b) 1:50 (red, larger water aggregates appear (trimers, tetramers, pentamers and hexamers); and (c) 1:25 (blue, a mixture of water aggregates and ASW).</figcaption></figure>
## Appendix E. FTIR Bands of Corannulene in Argon Matrices at 10 K

<table><caption>**Table A6.** Observed bands positions (cm⁻¹) and relative intensities (vs: very strong, s: strong, ms: medium strong, m: medium, mw: medium weak, w: weak) of corannulene isolated in argon matrices at 10 K and comparison with results by Rouillé et al. [39] in an argon matrix at 12 K Note that spectra of Rouillé et al. [39] also contained small amounts of water.</caption><thead><tr><th>Vib. Modes</th><th>Wavenumbers (cm<sup>-1</sup>, This Work)</th><th>from Ref. [39]</th></tr></thead><tbody><tr><td rowspan="5">ν(C-H)</td><td>3064.5 m</td><td></td></tr><tr><td>3056.1 m</td><td></td></tr><tr><td>3050.5 m</td><td></td></tr><tr><td>3043.0 m</td><td></td></tr><tr><td>3037.4 m</td><td></td></tr><tr><td rowspan="2">Combination modes</td><td>1906 w</td><td></td></tr><tr><td>1893 w</td><td></td></tr><tr><td rowspan="6">ν(C-C) and δ(C-H) in-plane bending</td><td>1441.0 m</td><td>1442.7</td></tr><tr><td></td><td>1428.8</td></tr><tr><td>1316.3 m</td><td>1316.9</td></tr><tr><td></td><td>1239.3</td></tr><tr><td>1138.8 mw</td><td>1138.6</td></tr><tr><td></td><td>1025.2</td></tr><tr><td rowspan="6">Out-of-plane bending γ(C-H))</td><td>837.1 vs</td><td>837.2</td></tr><tr><td>818.3 w</td><td>817.5</td></tr><tr><td>749.8 w</td><td>748.7</td></tr><tr><td>661.8 s</td><td>661.5</td></tr><tr><td>547.4 m</td><td></td></tr><tr><td>446.8 w</td><td></td></tr><tr><td></td><td>405.9 w</td><td></td></tr></tbody></table>

## References

1. Ghosh, A.K.; Chaudhuri, P.; Kumar, B.; Panja, S.S. Review on aggregation of asphaltene vis-a-vis spectroscopic studies *Fuel* **2016**, *185*, 541. [CrossRef]

2. Callén, M.; de la Cruz, M.; López, J.; Murillo, R.; Navarro, M.; Mastral, A. Some inferences on the mechanism of atmospheric gas/particle partitioning of polycyclic aromatic hydrocarbons (PAH) at Zaragoza (Spain). *Chemosphere* **2008**, *73*, 1357–1365. [CrossRef] [PubMed]

3. Vione, D.; Maurino, V.; Minero, C.; Pelizzetti, E.; Harrison, M.A.J.; Olariu, R.I.; Arsene, C. Photochemical reactions in the tropospheric aqueous phase and on particulate matter. *Chem. Soc. Rev.* **2006**, *35*, 441–453. [CrossRef] [PubMed]

4. Dartois, E. Interstellar Carbon Dust. C **2019**, *5*, 80. [CrossRef]

5. Allamandola, L.J.; Tielens, A.G.; Barker, J.R. Interstellar polycyclic aromatic hydrocarbons: The infrared emission bands, the excitation/emission mechanism, and the astrophysical implications. *Astrophys. J. Suppl. Ser.* **1989**, *71*, 733–775. [CrossRef] [PubMed]

6. Tielens, A. Interstellar Polycyclic Aromatic Hydrocarbon Molecules. *Ann. Rev. Astron. Astrophys.* **2008**, *46*, 289–337. [CrossRef]

7. Allamandola, L.J.; Tielens, A.G.G.M.; Barker, J.R. Polycyclic Aromatic-Hydrocarbons and the Unidentified Infrared-Emission Bands-Auto Exhaust Along the Milky-Way. *Astrophys. J.* **1985**, *290*, L25–L28. [CrossRef]

8. Léger, A.; Puget, J.L. Identification of the 'unidentified' IR emission features of interstellar dust? *Astron. Astrophys.* **1984**, *137*, L5–L8.

9. McGuire, B.A.; Loomis, R.A.; Burkhardt, A.M.; Lee, K.L.K.; Shingledecker, C.N.; Charnley, S.B.; Cooke, I.R.; Cordiner, M.A.; Herbst, E.; Kalenskii, S.; et al. Detection of two interstellar polycyclic aromatic hydrocarbons via spectral matched filtering. *Science* **2021**, *371*, 1265–1269. [CrossRef]

10. Hardegree-Ullman, E.E.; Gudipati, M.S.; Boogert, A.C.A.; Lignell, H.; Allamandola, L.J.; Stapelfeldt, K.R.; Werner, M. Laboratory determination of the infrared band strengths of pyrene frozen in water ice: Implications for the composition of interstellar ices. *Astrophys. J.* **2014**, *784*, 172. [CrossRef]

11. Tielens, A.G.G.M. *The Physics and Chemistry of the Interstellar Medium*; Cambridge University Press: Cambridge, UK, 2005.

12. Cook, A.M.; Ricca, A.; Mattioda, A.L.; Bouwman, J.; Roser, J.; Linnartz, H.; Bregman, J.; Allamandola, L.J. Photochemistry of polycyclic aromatic hydrocarbons in cosmic water ice: The role of ionization and concentration. *Astrophys. J.* **2015**, *799*, 14. [CrossRef]
13. Guennoun, Z.; Aupetit, C.; Mascetti, J. Photochemistry of coronene with water at 10 K: First tentative identification by infrared spectroscopy of oxygen containing coronene products. *Phys. Chem. Chem. Phys.* **2011**, *13*, 7340–7347. [CrossRef] [PubMed]

14. Guennoun, Z.; Aupetit, C.; Mascetti, J. Photochemistry of Pyrene with Water at Low Temperature: Study of Atmospherical and Astrochemical Interest. *J. Phys. Chem. A* **2011**, *115*, 1844–1852. [CrossRef] [PubMed]

15. Noble, J.A.; Michoulier, E.; Aupetit, C.; Mascetti, J. Influence of ice structure on the soft UV photochemistry of PAHs embedded in solid water. *Astron. Astrophys.* **2020**, *644*, A22. [CrossRef]

16. Lignell, A.; Tenelanda-Osorio, L.I.; Gudipati, M.S. Visible-light photoionization of aromatic molecules in water-ice: Organic chemistry across the universe with less energy. *Chem. Phys. Lett.* **2021**, *778*, 138814. [CrossRef]

17. Noble, J.A.; Jouvet, C.; Aupetit, C.; Moudens, A.; Mascetti, J. Efficient photochemistry of coronene:water complexes. *Astron. Astrophys.* **2017**, *599*, A124. [CrossRef]

18. Simon, A.; Noble, J.A.; Rouaut, G.; Moudens, A.; Aupetit, C.; Iftner, C.; Mascetti, J. Formation of coronene: Water complexes: FTIR study in argon matrices and theoretical characterisation. *Phys. Chem. Chem. Phys.* **2017**, *19*, 8516–8529. [CrossRef]

19. Amor, N.B.; Michoulier, E.; Simon, A. Electronic excited states of benzene in interaction with water clusters: Influence of structure and size. *Theor. Chem. Accounts* **2021**, *140*, 70. [CrossRef]

20. Deng, Y.; Yu, D.; Cao, X.; Liu, L.; Rong, C.; Lu, T.; Liu, S. Structure, aromaticity and reactivity of corannulene and its analogues: A conceptual density functional theory and density functional reactivity theory study. *Mol. Phys.* **2018**, *116*, 956–968. [CrossRef]

21. Cami, J.; Bernard-Salas, J.; Peeters, E.; Malek, S.E. Detection of C60 and C70 in a Young Planetary Nebula. *Science* **2010**, *329*, 1180–1182. [CrossRef]

22. Sellgren, K.; Werner, M.W.; Ingalls, J.G.; Smith, J.D.T.; Carleton, T.M.; Joblin, C. C60 in Reflection Nebulae. *Astrophys. J. Lett.* **2010**, *722*, L54. [CrossRef]

23. García-Hernández, D.A.; Manchado, A.; García-Lario, P.; Stanghellini, L.; Villaver, E.; Shaw, R.A.; Szczerba, R.; Perea-Calderon, J.V. Formation of fullerenes in H-containing planetary nebulae. *Astrophys. J. Lett.* **2010**, *724*, L39. [CrossRef]

24. García-Hernández, D.A.; Kameswara Rao, N.; Lambert, D.L. Are C-60 molecules detectable in circumstellar shells of R coronae berealis stars? *Astrophys. J.* **2011**, *729*, 126. [CrossRef]

25. García-Hernández, D.A.; Iglesias-Groth, S.; Acosta-Pulido, J.A.; Manchado, A.; García-Lario, P.; Stanghellini, L.; Villaver, E.; Shaw, R.A.; Cataldo, F. The formation of fullerenes: Clues from new C-60, C-70, and (possible) planar C-24 detections in magellanic cloud planetary nebulae. *Astrophys. J. Lett.* **2011**, *737*, L30. [CrossRef]

26. García-Hernández, D.A.; Villaver, E.; García-Lario, P.; Acosta-Pulido, J.A.; Manchado, A.; Stanghellini, L.; Shaw, R.A.; Cataldo, F. Infrared study of fullerene planetary nebulae. *Astrophys. J.* **2012**, *760*, 107. [CrossRef]

27. Berné, O.; Montillaud, J.; Joblin, C. Top-down formation of fullerenes in the interstellar medium. *Astron. Astrophys.* **2015**, *577*, A133. [CrossRef]

28. Zhen, J.; Castellanos, P.; Paardekooper, D.M.; Linnartz, H.; Tielens, A.G.G.M. Laboratory formation of fullerenes from PAHs: Top-down interstellar chemistry. *Astrophys. J.* **2014**, *797*, L30. [CrossRef]

29. Omont, A. Interstellar fullerene compounds and diffuse interstellar bands. *Astron. Astrophys.* **2016**, *590*, A52. [CrossRef]

30. Berné, O.; Tielens, A.G.G.M. Formation of buckminsterfullerene (C-60) in interstellar space. *Proc. Natl. Acad. Sci. USA* **2012**, *109*, 401. [CrossRef]

31. Dunk, P.W.; Kaiser, N.K.; Hendrickson, C.L.; Quinn, J.P.; Ewels, C.P.; Nakanishi, Y.; Sasaki, Y.; Shinohara, H.; Marshall, A.G.; Kroto, H.W. Closed network growth of fullerenes. *Nat. Comm.* **2012**, *3*, 855. [CrossRef]

32. Gavilan, L.; Alata, I.; Le, K.; Pino, T.; Giuliani, A.; Dartois, E. VUV spectroscopy of carbon dust analogs: Contribution to interstellar extinction. *Astron. Astrophys.* **2016**, *586*, A106. [CrossRef]

33. Dubosq, C.; Calvo, F.; Rapacioli, M.; Dartois, E.; Pino, T.; Falvo, C.; Simon, A. Quantum modeling of the optical spectra of carbon clusters structural families and relation to the interstellar extinction UV-bump. *Astron. Astrophys.* **2020**, *634*, A62. [CrossRef]

34. Carpentier, Y.; Feraud, G.; Dartois, E.; Brunetto, R.; Charon, E.; Cao, A.T.; d'Hendecourt, L.; Brechignac, P.; Rouzaud, J.N.; Pino, T. Nanostructuration of carbonaceous dust as seen through the positions of the 6.2 and 7.7 µm AIBs. *Astron. Astrophys.* **2012**, *548*, A40. [CrossRef]

35. Gavilan, L.; Le, K.C.; Pino, T.; Alata, I.; Giuliani, A.; Dartois, E. Polyaromatic disordered carbon grains as carriers of the UV bump: Far-UV to mid-IR spectroscopy of laboratory analogs. *Astron. Astrophys.* **2017**, *607*, A73. [CrossRef]

36. Dubosq, C.; Falvo, C.; Calvo, F.; Rapacioli, M.; Parneix, P.; Pino, T.; Simon, A. Mapping the structural diversity of C60 carbon clusters and their infrared response. *Astron. Astrophys.* **2019**, *625*, L11. [CrossRef]

37. Calvo, F.; Simon, A.; Parneix, P.; Falvo, C.; Dubosq, C. Infrared Spectroscopy of Chemically Diverse Carbon Clusters: A Data-Driven Approach. *J. Phys. Chem. A* **2021**, *125*, 5509–5518. [CrossRef]

38. Pilleri, P.; Herberth, D.; Giesen, T.F.; Gerin, M.; Joblin, C.; Mulas, G.; Malloci, G.; Grabow, J.U.; Brünken, S.; Surin, L.; et al. Search for corannulene (C20H10) in the Red Rectangle. *Mon. Not. R. Astron. Soc.* **2009**, *397*, 1053–1060. [CrossRef]

39. Rouillé, G.; Jäger, C.; Steglich, M.; Huisken, F.; Henning, T.; Theumer, G.; Bauer, I.; Knölker, H.J. IR, Raman, and UV/Vis Spectra of Corannulene for Use in Possible Interstellar Identification. *ChemPhysChem* **2008**, *9*, 2085–2091. [CrossRef]

40. Zhao, L.; Doddipatla, S.; Kaiser, R.I.; Lu, W.; Kostko, O.; Ahmed, M.; Tuli, L.B.; Morozov, A.N.; Howlader, A.H.; Wnuk, S.F.; et al. Gas-phase synthesis of corannulene—A molecular building block of fullerenes. *Phys. Chem. Chem. Phys.* **2021**, *23*, 5740–5749. [CrossRef]
41. Stuparu, M.C. Corannulene: A Curved Polyarene Building Block for the Construction of Functional Materials. *Acc. Chem. Res.* **2021**, *54*, 2858–2870. [CrossRef]

42. Hernández-Rojas, J.; Calvo, F.; Rabilloud, F.; Bretón, J.; Gomez Llorente, J.M. Modeling Water Clusters on Cationic Carbonaceous Seeds. *J. Phys. Chem. A* **2010**, *114*, 7267–7274. [CrossRef] [PubMed]

43. Pérez, C.; Steber, A.L.; Rijs, A.M.; Temelso, B.; Shields, G.C.; Lopez, J.C.; Kisiel, Z.; Schnell, M. Corannulene and its complex with water: A tiny cup of water. *Phys. Chem. Chem. Phys.* **2017**, *19*, 14214–14223. [CrossRef] [PubMed]

44. Heine, T.; Rapacioli, M.; Patchkovskii, S.; Frenzel, J.; Koster, A.; Calaminici, P.; Duarte, H.A.; Escalante, S.; Flores-Moreno, R.; Goursot, A.; et al. deMonNano. 2009. Available online: http://demon-nano.ups-tlse.fr/ (accessed on 27 February 2022).

45. Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-consistent-charge density-functional tight-binding method for simulations of complex material properties. *Phys. Rev. B* **1998**, *58*, 7260–7268. [CrossRef]

46. Li, J.; Zhu, T.; Cramer, C.; Truhlar, D. New Class IV Charge Model for Extracting Accurate Partial Charges from Wave Functions. *J. Phys. Chem. A* **1998**, *102*, 1820–1831. [CrossRef]

47. Rapacioli, M.; Spiegelman, F.; Talbi, D.; Mineva, T.; Goursot, A.; Heine, T.; Seifert, G. Correction for dispersion and Coulombic interactions in molecular clusters with density functional derived methods: Application to polycyclic aromatic hydrocarbon clusters. *J. Chem. Phys.* **2009**, *130*, 244304. [CrossRef]

48. Simon, A.; Rapacioli, M.; Mascetti, J.; Spiegelman, F. Vibrational spectroscopy and molecular dynamics of water monomers and dimers adsorbed on polycyclic aromatic hydrocarbons. *Phys. Chem. Chem. Phys.* **2012**, *14*, 6771–6786. [CrossRef]

49. Kalinowski, J.; Lesyng, B.; Thompson, J.; Cramer, C.; Truhlar, D. Class IV Charge Model for the Self-Consistent Charge Density-Functional Tight-Binding Method. *J. Phys. Chem. A* **2004**, *108*, 2545–2549. [CrossRef]

50. Joalland, B.; Rapacioli, M.; Simon, A.; Joblin, C.; Marsden, C.J.; Spiegelman, F. Molecular Dynamics Simulations of Anharmonic Infrared Spectra of [SiPAH](+) pi-Complexes. *J. Phys. Chem. A* **2010**, *114*, 5846–5854. [CrossRef]

51. Simon, A.; Spiegelman, F. Water clusters adsorbed on polycyclic aromatic hydrocarbons: Energetics and conformational dynamics. *J. Chem. Phys.* **2013**, *138*, 194309. [CrossRef]

52. Simon, A.; Spiegelman, F. Conformational dynamics and finite-temperature infrared spectra of the water octamer adsorbed on coronene. *Comput. Theor. Chem.* **2013**, *102*, 54–61. [CrossRef]

53. Michoulier, E.; Toubin, C.; Simon, A.; Mascetti, J.; Aupetit, C.; Noble, J.A. Perturbation of the Surface of Amorphous Solid Water by the Adsorption of Polycyclic Aromatic Hydrocarbons. *J. Phys. Chem. C* **2020**, *124*, 2994–3001. [CrossRef]

54. Aziz, R.A. A highly accurate interatomic potential for argon. *J. Chem. Phys.* **1993**, *99*, 4518–4525. [CrossRef]

55. Iftner, C.; Simon, A.; Korchagina, K.; Rapacioli, M.; Spiegelman, F. A density functional tight binding/force field approach to the interaction of molecules with rare gas clusters: Application to $(C_6H_6)^{+/0}Ar_n$ clusters. *J. Chem. Phys.* **2014**, *140*, 034301. [CrossRef] [PubMed]

56. Simon, A.; Iftner, C.; Mascetti, J.; Spiegelman, F. 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. *J. Phys. Chem. A* **2015**, *119*, 2449. [CrossRef]

57. Crépin, C.; Pujo, P.d.; Bouvier, B.; Brenner, V.; Millié, P. A simulation of naphthalene matrix isolation: Comparison with experiments. *Chem. Phys.* **2001**, *272*, 243–258. [CrossRef]

58. Ceponkus, J.; Uvdal, P.; Nelander, B. Acceptor Switching and Axial Rotation of the Water Dimer in Matrices, observed by Infrared Spectroscopy. *J. Chem. Phys.* **2010**, *133*, 074301. [CrossRef]

59. Perchard, J. Anharmonicity and Hydrogen Bonding. III. Analysis of the Near Infrared Spectrum of Water Trapped in Argon Matrix. *Chem. Phys.* **2001**, *273*, 217–233. [CrossRef]

60. Michaut, X.; Vasserot, A.M.; Abouaf-Marguin, L. Temperature and Time Effects on the Rovibrational Structure of Fundamentals of H₂O Trapped in Solid Argon: Hindered Rotation and RTC Satellite. *Vib. Spec.* **2004**, *34*, 83–93. [CrossRef]

61. Ceponkus, J.; Karlstro, G.; Nelander, B. Intermolecular Vibrations of the Water Trimer, a Matrix Isolation Study. *J. Chem. Phys. A* **2005**, *109*, 7859–7864. [CrossRef]

62. Bentwood, R.; Barnes, A.; Orville-Thomas, W. Studies of intermolecular interactions by matrix isolation vibrational spectroscopy: Self-association of water. *J. Mol. Spectrosc.* **1980**, *84*, 391–404. [CrossRef]

63. Buck, U.; Huisken, F. Infrared Spectroscopy of Size-Selected Water and Methanol Clusters. *Chem. Rev.* **2000**, *100*, 3863–3890. [CrossRef]

64. Ceponkus, J.; Uvdal, P.; Nelander, B. Water Tetramer, Pentamer, and Hexamer in Inert Matrices. *J. Phys. Chem. A* **2012**, *116*, 4842–4850. [CrossRef]

65. Ceponkus, J.; Engdahl, A.; Uvdal, P.; Nelander, B. Structure and Dynamics of Small Water Clusters, Trapped in Inert Matrices. *Chem. Phys. Lett.* **2013**, *581*, 1–9. [CrossRef]

66. Coussan, S.; Loutellier, A.; Perchard, J.; Racine, S.; Bouteiller, Y. Matrix isolation infrared spectroscopy and DFT calculations of complexes between water and nitrogen. *J. Mol. Struct.* **1998**, *471*, 37–47. [CrossRef]

67. Hirabayashi, S.; Yamada, K.M. Infrared Spectra and Structure of Water Clusters Trapped in Argon and Krypton Matrices. *J. Mol. Struct.* **2006**, *795*, 78–83. [CrossRef]
68. Bauschlicher, C.W.; Ricca, A.; Boersma, C.; Allamandola, L.J. The NASA Ames PAH IR Spectroscopic Database: Computational Version 3.00 with Updated Content and the Introduction of Multiple Scaling Factors. *Astrophys. J. Suppl. Ser.* **2018**, *234*, 32. [CrossRef]

69. Boersma, C.; Bauschlicher, C.W.; Ricca, A.; Mattioda, A.L.; Cami, J.; Peeters, E.; de Armas, F.S.; Saborido, G.P.; Hudgins, D.M.; Allamandola, L.J. The NASA AMES PAH spectroscopic database version 2.00: Updated content, web site and on off line tools. *Astrophys. J. Suppl. Ser.* **2014**, *211*, 8. [CrossRef]

70. Mattioda, A.L.; Hudgins, D.M.; Boersma, C.; Bauschlicher, C.W.; Ricca, A.; Cami, J.; Peeters, E.; de Armas, F.S.; Saborido, G.P.; Allamandola, L.J. The NASA Ames PAH IR Spectroscopic Database: The Laboratory Spectra. *Astrophys. J. Suppl. Ser.* **2020**, *251*, 22. [CrossRef]

71. Vauthey, E. Elucidating the Mechanism of Bimolecular Photoinduced Electron Transfer Reactions. *J. Phys. Chem. B* **2022**, *126*, 778–788. [CrossRef]