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

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+.

Evolution of the ionisation energy with the stepwise growth of chiral clusters of [4]helicene

Polycyclic aromatic hydrocarbons (PAHs) are widely established as ubiquitous in the interstellar medium (ISM), but considering their prevalence in harsh vacuum environments, the role of ionisation in the formation of PAH clusters is poorly understood, particularly if a chirality-dependent aggregation route is considered. Here we report on photoelectron spectroscopy experiments on [4]helicene clusters performed with a vacuum ultraviolet synchrotron beamline. Aggregates (up to the heptamer) of [4]helicene, the smallest PAH with helical chirality, were produced and investigated with a combined experimental and theoretical approach using several state-of-the-art quantum-chemical methodologies. The ionisation onsets are extracted for each cluster size from the mass-selected photoelectron spectra and compared with calculations of vertical ionisation energies. We explore the complex aggregation topologies emerging from the multitude of isomers formed through clustering of P and M, the two enantiomers of [4]helicene. The very satisfactory benchmarking between experimental ionisation onsets vs. predicted ionisation energies allows the identification of theoretically predicted potential aggregation motifs and corresponding energetic ordering of chiral clusters. Our structural models suggest that a homochiral aggregation route is energetically favoured over heterochiral arrangements with increasing cluster size, hinting at potential symmetry breaking in PAH cluster formation at the scale of small grains.

Addressing electronic and dynamical evolution of molecules and molecular clusters: DFTB simulations of energy relaxation in polycyclic aromatic hydrocarbons

We present a review of the capabilities of the density functional based Tight Binding (DFTB) scheme to address the electronic relaxation and dynamical evolution of molecules and molecular clusters following energy deposition via either collision or photoabsorption. The basics and extensions of DFTB for addressing these systems and in particular their electronic states and their dynamical evolution are reviewed. Applications to PAH molecules and clusters, carbonaceous systems of major interest in astrochemical/astrophysical context, are reported. A variety of processes are examined and discussed such as collisional hydrogenation, fast collisional processes and induced electronic and charge dynamics, collision-induced fragmentation, photo-induced fragmentation, relaxation in high electronic states, electronic-to-vibrational energy conversion and statistical versus non-statistical fragmentation. This review illustrates how simulations may help to unravel different relaxation mechanisms depending on various factors such as the system size, specific electronic structure or excitation conditions, in close connection with experiments.

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

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

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.

Water Attachment onto Size-Selected Cationic Pyrene Clusters

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

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

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

Simulation of Liquids with the Tight-Binding Density-Functional Approach and Improved Atomic Charges

Theoretical description of liquids, especially liquid water, is an ongoing subject with important implications in various domains such as homogeneous catalysis; solvation of molecular, ionic, and biomolecular species; and reactivity. Various formalisms exist to describe liquids, each one displaying its own balance between accuracy and computational cost that defines its range of applications. The present article revisits the ability of the density-functional-based tight-binding (SCC-DFTB) approach to model liquids by focusing on liquid water and liquid benzene under ambient conditions. To do so, we benchmark a recent correction for the SCC-DFTB atomic charges that allows for a drastic improvement of the pair radial distribution functions of liquid water as compared to both experimental data and density-functional theory results performed in the generalized-gradient approximation. We also report the coupling of the deMonNano and i-PI codes to perform path-integral molecular dynamics. This allows us to rationalize the impact of nuclear quantum effects on the SCC-DFTB description of liquid water. This study evidences the rather good ability of SCC-DFTB to describe liquid water and liquid benzene. As the first example of application, we also present results for a benzene molecule solvated in water with the perspectives of further studies devoted to solvent/water interfaces.

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.

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.

Photo-dissociation of naphthalene dimer cations stored in a compact electrostatic ion storage ring

Naphthalene dimer cations [C₁₀H₈]₂ ⁺ have been produced by using an electron cyclotron resonance plasma ion source and stored in a compact electrostatic ion storage ring. We show that the radiative cooling of these cations is much slower than the isolated monomer naphthalene cations. We also report on photo-dissociation studies in the gas phase of naphthalene dimer cations at high internal energy. The dissociation energy is estimated to 0.5 eV in close agreement with previous measurements but a factor of 2 smaller than recent (density functional theory (DFT) and ab initio) theoretical studies. As uncertainties on theory as well as on the experiment cannot be as large as this difference, we conclude that this discrepancy may be due to temperature effects with possible isomerization. As an interpretation of the photo-dissociation spectrum of naphthalene dimer cations, we propose a tentative simple analytical model based on effective Morse potentials. These effective potentials are expected to "average" temperature effects that would apparently result in a smaller energy difference between the fundamental and dissociation states due to the twisting vibration modes of the naphthalene dimer cations.

Density-functional tight-binding approach for metal clusters, nanoparticles, surfaces and bulk: application to silver and gold

Density-functional based tight-binding (DFTB) is an efficient quantum mechanical method that can describe a variety of systems, going from organic and inorganic compounds to metallic and hybrid materials. The present topical review addresses the ability and performance of DFTB to investigate energetic, structural, spectroscopic and dynamical properties of gold and silver materials. After a brief overview of the theoretical basis of DFTB, its parametrization and its transferability, we report its past and recent applications to gold and silver systems, including small clusters, nanoparticles, bulk and surfaces, bare and interacting with various organic and inorganic compounds. The range of applications covered by those studies goes from plasmonics and molecular electronics, to energy conversion and surface chemistry. Finally, perspectives of DFTB in the field of gold and silver surfaces and NPs are outlined.

Size Effect in the Ionization Energy of PAH Clusters

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