Near-infrared absorption and radiative cooling of naphthalene dimers (C10H8)2.

The radiative cooling of naphthalene dimer cations, (C10H8)2+ was studied experimentally through action spectroscopy using two different electrostatic ion-beam storage rings, DESIREE in Stockholm and Mini-Ring in Lyon. The spectral characteristics of the charge resonance (CR) band were observed to vary significantly with a storage time of up to 30 seconds in DESIREE. In particular, the position of the CR band shifts to the blue, with specific times (inverse of rates) of 0.64 s and 8.0 s in the 0-5 s and 5-30 s storage time ranges, respectively. These long-time scales indicate that the internal energy distribution of the stored ions evolves by vibrational radiative cooling, which is consistent with the absence of fast radiative cooling via recurrent fluorescence for (C10H8)2+. Density functional based tight binding calculations with local excitations and configuration interactions (DFTB-EXCI) were used to simulate the absorption spectrum for ion temperatures between 10 and 500 K. The evolution of the bandwidth and position with temperature is in qualitative agreement with the experimental findings. Furthermore, these calculations yielded linear temperature dependencies for both the shift and the broadening. Combining the relationship between the CR band position and the ion temperature with the results of the statistical model, we demonstrate that the observed blue shift can be used to determine the radiative cooling rate of (C10H8)2+.

DFTB Simulation of Charged Clusters Using Machine Learning Charge Inference.

We present a modification to self-consistent charge density functional-based tight binding (SCC-DFTB), which allows computation based on approximate atomic charges. We obtain these charges by means of a machine learning (ML) process that combines a Coulomb model with a neural network. This allows us to avoid the SCC cycles in the SCC-DFTB calculation while maintaining its accuracy. The main input of the model is the atomic positions characterized by a set of atom-centered symmetry functions. The charge inference from our ML algorithm is as close as 10-2 units of charge from the exact SCC solution. Our ML-DFTB approach provides a good approximation of the density matrix and of the energy and forces with only a single diagonalization. This is a significant computational saving with respect to the complete SCC algorithm, which allows us to investigate a bigger ensemble of atoms. We show the quality of our approach in the case of charged silicon carbide (SiC) clusters. The ML-DFTB potential energy surface (PES) mimics the SCC-DFTB PES rather well, despite its simplicity. This allows us to obtain the same geometric structure ordering with respect to energy for small clusters. The dissociation barriers for ion emission are well-reproduced, which opens the way to investigating ion field emission and charged cluster stability. The ML-DFTB approach is obviously not limited to charged clusters or SiC materials. It opens a new route to investigate larger clusters than those investigated by standard SCC-DFTB, as well as surface and solid-state chemistry at the atomic level.

Electronic effects in the dissociative ionisation of pyrene clusters.

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

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.

Threshold collision induced dissociation of pyrene cluster cations.

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

Dependence upon charge of the vibrational spectra of small Polycyclic Aromatic Hydrocarbon clusters: the example of pyrene.

Infrared spectra are computed for neutral and cationic clusters of Polycyclic Aromatic Hydrocarbon clusters, namely ( C 16 H 10 ) n = 1 , 4 ( 0 / + ) , using the Density Functional based Tight Binding scheme combined with a Configuration Interaction scheme (DFTB-CI) in the double harmonic approximation. Cross-comparison is carried out with DFT and simple DFTB. Similarly to the monomer cation, the IR spectra of cluster cations are characterized by a depletion of the intensity of the CH stretch modes around 3000 cm-1, with a weak revival for n = 3 and 4. The in-plane CCC modes in the region 1400-2000 cm-1 are enhanced while the CH bending modes in the range 700-1000 cm-1 are significantly weakened with respect to the monomer cation, in particular for n = 2. Finally, soft modes corresponding to diedral fluctuations of the monomers within the central stack of the ion structure, possibly mixed with monomer folding, are also observed in the region 70-120 cm-1.

Electron Diffraction of Pyrene Nanoclusters Embedded in Superfluid Helium Droplets.

We report electron diffraction of pyrene nanoclusters embedded in superfluid helium droplets. Using a least-squares fitting procedure, we have been able to separate the contribution of helium from those of the pyrene nanoclusters and determine the most likely structures for dimers and trimers. We confirm that pyrene dimers form a parallel double-layer structure with an interlayer distance of 3.5 Šand suggest that pyrene trimers form a sandwich structure but that the molecular planes are not completely parallel. The relative contributions of the dimers and trimers are ∼6:1. This work is an extension of our effort of solving structures of biological molecules using serial single-molecule electron diffraction imaging. The success of electron diffraction from an all-light-atom sample embedded in helium droplets offers reassuring evidence of the feasibility of this approach.

Thermal evaporation of pyrene clusters.

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

Structures and Energetics of Neutral and Cationic Pyrene Clusters.

The low energy structures of neutral and cationic pyrene clusters containing up to seven molecules are searched through a global exploration scheme combining parallel tempering Monte Carlo algorithm and local quenches. The potential energies are computed at the density functional based tight binding level for neutrals and configuration interaction density functional based tight binding for cations in order to treat properly the charge resonance. New simplified versions of these schemes are also presented and used during the global exploration. Neutral clusters are shown to be made of compact assemblies of sub-blocs containing up to three units whereas cations present a charged dimer or trimer core surrounded by neutral units. The structural features of the clusters are analyzed and correlated for the cation with the charge distribution. The stability of clusters is also discussed in terms of cohesive and evaporation energies. Adiabatic and vertical ionization potentials are also discussed.

An extended DFTB-CI model for charge-transfer excited states in cationic molecular clusters: model studies versus ab initio calculations in small PAH clusters.

We present an extension of the constrained density functional tight binding scheme combined with configuration interaction (DFTB-CI) to efficiently compute excited states of molecular cluster cations and their oscillator strengths from the ground state. The present extension consists of generalizing the initial model, relying on configurations with holes in the monomer HOMOs only, to configurations involving sub-HOMO holes, allowing for the description of higher excited states. The extended scheme is benchmarked on selected energy pathways with respect to available ab initio and new CASPT2 reference calculations on the benzene, naphthalene and pyrene dimer cations. The ability of the model to describe the potential energy surfaces and the transition dipole moments is discussed. The vertical electronic absorption spectra of the three dimer cations are calculated and compared with the theoretical litterature and available experimental data. Finally, the electronic absorption spectra of low energy isomers of the trimer and tetramer pyrene cluster cations are also predicted.

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

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