Exploring Molecular Energy Landscapes by Coupling the DFTB Potential with a Tree-Based Stochastic Algorithm: Investigation of the Conformational Diversity of Phthalates.

Exploring the global energy landscape of relatively large molecules at the quantum level is a challenging problem. In this work, we report the coupling of a nonredundant conformational space exploration method, namely, the robotics-inspired iterative global exploration and local optimization (IGLOO) algorithm, with the quantum-chemical density functional tight binding (DFTB) potential. The application of this fast and efficient computational approach to three close-sized molecules of the phthalate family (DBP, BBP, and DEHP) showed that they present different conformational landscapes. These differences have been rationalized by making use of descriptors based on distances and dihedral angles. Coulomb interactions, steric hindrance, and dispersive interactions have been found to drive the geometric properties. A strong correlation has been evidenced between the two dihedral angles describing the side-chain orientation of the phthalate molecules. Our approach identifies low-energy minima without prior knowledge of the potential energy surface, paving the way for future investigations into transition paths and states.

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.

Temperature driven transformations of glycine molecules embedded in interstellar ice.

The formation of glycine amino acid on ice grains in space raises fundamental questions about glycine chemistry in interstellar media. In this work, we studied glycine conformational space and the related tautomerization mechanisms in water media by means of QM/MM molecular dynamics simulations of four glycine conformational isomers (cc, ct, tc, and tt). Interstellar low density amorphous (LDA) ice and T = 20 K were considered as representative for a cold interstellar ice environment, while temperatures of 250 and 450 K were included to model rapid local heating in the ice. In addition to the LDA environment, water clusters with 4, 17, and 27 H2O molecules were subjected to QM/MM dynamics simulations that allowed glycine tautomerization behaviour to be evaluated in water surface-like environments. The tautomerization processes were found to be strongly dependent on the number of water molecules and specific isomer structure. All the glycine isomers mostly preserve their canonical "neutral" conformations under interstellar conditions.

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

Probing the electronic structure and ground state symmetry of gas phase C60+ via VUV photoionization and comparison with theory.

Recently, some of us reviewed and studied the photoionization dynamics of C60 that are of great interest to the astrochemical community as four of the diffuse interstellar bands (DIBs) have been assigned to electronic transitions in the C60+ cation. Our previous analysis of the threshold photoelectron spectrum (TPES) of C60 [Hrodmarsson et al., Phys. Chem. Chem. Phys. 22, 13880-13892 (2020)] appeared to give indication of D3d ground state symmetry, in contrast to theoretical predictions of D5d symmetry. Here, we revisit our original measurements taking account of a previous theoretical spectrum presented in the work of Manini et al., Phys. Rev. Lett. 91(19), 196402 (2003), obtained within a vibronic model parametrized on density functional theory/local-density approximation electronic structure involving all hg Jahn-Teller active modes, which couple to the 2Hu components of the ground state of the C60+ cation. By reanalyzing our measured TPES of the ground state of the C60 Buckminsterfullerene, we find a striking resemblance to the theoretical spectrum calculated in the work of Manini et al., Phys. Rev. Lett. 91(19), 196402 (2003), and we provide assignments for many of the hg modes. In order to obtain deeper insights into the temperature effects and possible anharmonicity effects, we provide complementary modeling of the photoelectron spectrum via classical molecular dynamics (MD) involving density functional based tight binding (DFTB) computations of the electronic structure for both C60 and C60+. The validity of the DFTB modeling is first checked vs the IR spectra of both species which are well established from IR spectroscopic studies. To aid the interpretation of our measured TPES and the comparisons to the ab initio spectrum we showcase the complementarity of utilizing MD calculations to predict the PES evolution at high temperatures expected in our experiment. The comparison with the theoretical spectrum presented in the work of Manini et al., Phys. Rev. Lett. 91(19), 196402 (2003), furthermore, provides further evidence for a D5d symmetric ground state of the C60+ cation in the gas phase, in complement to IR spectroscopy in frozen noble gas matrices. This not only allows us to assign the first adiabatic ionization transition and thus determine the ionization energy of C60 with greater accuracy than has been achieved at 7.598 ± 0.005 eV, but we also assign the two lowest excited states (2E1u and 2E2u) which are visible in our TPES. Finally, we discuss the energetics of additional DIBs that could be assigned to C60+ in the future.

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.

Low-Energy Transformation Pathways between Naphthalene Isomers.

The transformation pathways between low-energy naphthalene isomers are studied by investigating the topology of the energy landscape of this astrophysically relevant molecule. The threshold algorithm is used to identify the minima basins of the isomers on the potential energy surface of the system and to evaluate the probability flows between them. The transition pathways between the different basins and the associated probabilities were investigated for several lid energies up to 11 eV, this value being close to the highest photon energy in the interstellar medium (13.6 eV). More than a hundred isomers were identified and a set of 23 minima was selected among them, on the basis of their energy and probability of occurrence. The return probabilities of these 23 minima and the transition probabilities between them were computed for several lid energies up to 11 eV. The first connection appeared at 3.5 eV while all minima were found to be connected at 9.5 eV. The local density of state was also sampled inside their respective basins. This work gives insight into both energy and entropic barriers separating the different basins, which also provides information about the transition regions of the energy landscape.

Collision-induced dissociation of protonated uracil water clusters probed by molecular dynamics simulations.

Collision-induced dissociation experiments of hydrated molecular species can provide a wealth of important information. However, they often need a theoretical support to extract chemical information. In the present article, in order to provide a detailed description of recent experimental measurements [Braud et al., J. Chem. Phys., 2019, 150, 014303], collision simulations between low-energy protonated uracil water clusters (H2O)1-7,11,12UH+ and an Ar atom were performed using a quantum mechanics/molecular mechanics formalism based on the self-consistent-charge density-functional based tight-binding method. The theoretical proportion of formed neutral vs. protonated uracil containing clusters, total fragmentation cross sections as well as the mass spectra of charged fragments are consistent with the experimental data which highlights the accuracy of the present simulations. They allow to probe which fragments are formed on the short time scale and rationalize the location of the excess proton on these fragments. We demonstrate that this latter property is highly influenced by the nature of the aggregate undergoing the collision. Analyses of the time evolution of the fragments populations and of their relative abundances demonstrate that, up to 7 water molecules, a direct dissociation mechanism occurs after collision whereas for 11 and 12 water molecules a statistical mechanism is more likely to participate. Although scarce in the literature, the present simulations appear as a useful tool to complement collision-induced dissociation experiments of hydrated molecular species.

Hydrogenation of C24 Carbon Clusters: Structural Diversity and Energetic Properties.

This work aims at exploring the potential energy surfaces of C24Hn=0,6,12,18,24 using the genetic algorithm in combination with the density functional based tight binding potential. The structural diversity was analyzed using order parameters, in particular the sum of the numbers of 5- and 6-carbon rings R5/6. The most abundant and lowest energy population was designated as the flake population (isomers of variable shapes, large R5/6 values), characterized by an increasing number of spherical isomers when nH/nC increases. Simultaneously, the fraction of the pretzel population (spherical isomers, smaller R5/6 values) increases. The fraction of the cage population (largest R5/6 values) remains extremely minor while the branched population (smallest R5/6 values) remains the highest energy population for all nH/nC ratios. For all C24Hn=0,6,12,18,24 clusters, the evolution of the carbon ring size distribution with energy clearly shows the correlation between the stability and the number of 6-carbon rings. The average values of the ionization potentials of all populations were found to decrease when nH/nC increases, ranging from 7.9 down to 6.4 eV. This trend was correlated to geometric and electronic factors, in particular to carbon hybridization. These results are of astrophysical interest, especially regarding the role of carbon species in the gas ionization.

Theoretical investigation of the electronic relaxation in highly excited chrysene and tetracene: The effect of armchair vs zigzag edge.

Non-adiabatic molecular dynamics of neutral chrysene and tetracene molecules is investigated using Tully's fewest switches surface hopping algorithm coupled to the time-dependent density functional based tight-binding (TD-DFTB) method for electronic structure calculations. We first assess the performance of two DFTB parameter sets based on the computed TD-DFTB absorption spectra. The main focus is given to the analysis of the electronic relaxation from the brightest excited state following absorption of a UV photon. We determine the dynamical relaxation times and discuss the underlying mechanisms. Our results show that the electronic population of the brightest excited singlet state in armchair-edge chrysene decays an order-of-magnitude faster than the one in zigzag-edge tetracene. This is correlated with a qualitatively similar difference of energy gaps between the brightest state and the state lying just below in energy, which is also consistent with our previous study on polyacenes.

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.

Evidencing the relationship between isomer spectra and melting: the 20- and 55-atom silver and gold cluster cases.

The present work highlights the links between melting properties and structural excitation spectra of small gold and silver clusters. The heat capacity curves are computed for Ag20, Au20, Ag55, Au55 and their ions, using a parallel-tempering molecular dynamics scheme to explore the density functional based tight binding (DFTB) potential energy surfaces and the multiple histogram method. It is found that clusters having very symmetric lowest energy structures (Au20, Ag55 and their ions) present sharp or relatively sharp solid-to-liquid transitions and large melting temperatures, important structural excitation energies and a discrete isomer spectrum. Opposite trends are observed for less ordered clusters (Ag20, Au55 and their ions). Regarding the structural evolution with temperature, very symmetric clusters exhibit minor evolution up to the starting melting temperature. The present study also highlights that, in contrast with the case of Au20, a single electron excess or deficiency is not determinant regarding the melting characteristics, even quantitatively, for clusters containing 55 atoms, for gold as for silver.

Non-adiabatic molecular dynamics investigation of the size dependence of the electronic relaxation in polyacenes.

The Tully's fewest switches surface hopping algorithm is implemented within the framework of the time-dependent density functional based tight binding method (TD-DFTB) to simulate the energy relaxation following absorption of a UV photon by polycyclic aromatic hydrocarbons (PAHs). This approach is used to study the size effect on the ultrafast dynamics in excited states for a special class of PAH species called polyacenes. We determine the dynamical relaxation times and discuss the underlying mechanisms. Our results show that there is a striking alternation in decay times of the brightest singlet state for neutral polyacenes with 3 to 6 aromatic cycles. The alternation corresponds to an order-of-magnitude variation between roughly 10 and 100 fs and is correlated with a qualitatively similar alternation of energy gaps between the brightest state and the state lying just below in energy.

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.

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.

Atomic hydrogen interactions with gas-phase coronene cations: hydrogenation versus fragmentation.

Sequential hydrogenation of polycyclic aromatic hydrocarbon (PAH) cations drives a gradual transition from a planar to a puckered geometry and from an aromatic to an aliphatic electronic structure. The resulting H-induced weakening of the molecular structure together with the exothermic nature of the consecutive H-attachment processes can lead to substantial molecular fragmentation. We have studied H attachment to gas-phase coronene cations in a radiofrequency ion trap using tandem mass spectrometry. With increasing hydrogenation, C2Hi loss and multifragmentation are identified as main de-excitation channels. To understand the dependence of both channels on H-exposure time, we have simulated the molecular stability and fragmentation channels of hydrogenated PAHs using a molecular dynamics approach employing potential energies determined by a density functional based tight binding method. As the coronene fragmentation patterns depend on the balance between energy deposition by H-attachment and the extent of cooling in between subsequent attachment processes, we investigate several scenarios for the energy distribution of hydrogenated PAHs. Good agreement between experiment and simulation is reached, when realistic energy distributions are considered.

Theoretical determination of adsorption and ionisation energies of polycyclic aromatic hydrocarbons on water ice.

In dense interstellar environments, Polycyclic Aromatic Hydrocarbons (PAHs) are likely to condense onto or integrate into water ice mantles covering dust grains. Understanding the role of ice in the photo-induced processes involving adsorbed PAHs is therefore a key issue in astrochemistry. This requires (i) the knowledge of PAH-ice interactions, i.e. PAH-ice adsorption energies and local structures at the PAH-ice interface, as well as (ii) the understanding of the fate of electrons in the PAH-ice system upon excitation. Regarding (i), in this work, we determined the lowest energy structures of PAH-ice systems for a variety of PAHs ranging from naphthalene to ovalene on three types of ice - crystalline (Ih and Ic) and amorphous (low density) - using an explicit description of the electrons and a finite-sized system. The electronic structure was determined using the Self Consistent Charge Density Functional based Tight Binding (SCC-DFTB) scheme with modified Mulliken charges in order to ensure a good description of the studied systems. Regarding (ii), the influence of the interaction with ice on the Vertical Ionisation Potentials (VIPs) of the series of PAHs was determined using the constrained SCC-DFTB scheme benchmarked against correlated wavefunction results for PAH-(H2O)n (n = 1-6, 13) clusters. The results show a deviation equal, at most, to ∼1.4 eV of the VIPs of PAHs adsorbed on ice with respect to the gas phase values. Our results are discussed in the light of experimental data and previous theoretical studies.

Melting of the Au20 Gold Cluster: Does Charge Matter?

We investigate the dependence upon the charge of the heat capacities of magic gold cluster Au20 obtained from density functional-based tight binding theory within parallel tempering molecular dynamics and the multiple histogram method. The melting temperatures, determined from the heat capacity curves, are found to be 1102 K for neutral Au20 and only 866 and 826 K for Au20+ and Au20-. Both the canonical and the microcanonical aspects of the transition are discussed. A convex intruder, associated with a negative heat capacity, is evidenced in the microcanonical entropy in the case of the neutral cluster but is absent in the melting processes of the ions. The present work shows that a single charge quantitatively affects the thermal properties of the gold 20mer.

Au147 nanoparticles: Ordered or amorphous?

Structural aspects of the Au147 cluster have been investigated through a density functional based tight binding global optimization involving a parallel tempering molecular dynamics scheme with quenching followed by geometries relaxation at the Density Functional Theory (DFT) level. The focus is put on the competition between relaxed ordered regular geometries and disordered (or amorphous) structures. The present work shows that Au147 amorphous geometries are relevant low energy candidates and are likely to contribute in finite temperature dynamics and thermodynamics. The structure of the amorphous-like isomers is discussed from the anisotropy parameters, the atomic coordinations, the radial and pair distribution functions, the IR spectra, and the vibrational DOS. With respect to the regular structures, the amorphous geometries are shown to be characterized by a larger number of surface atoms, a less dense volume with reduced coordination number per atom, a propensity to increase the dimension of flat facets at the surface, and a stronger anisotropy. Moreover, all amorphous clusters have similar IR spectra, almost continuous with active frequencies over the whole spectral range, while symmetric clusters are characterized by a few lines with large intensities.

Theoretical investigation of the solid-liquid phase transition in protonated water clusters.

Protonated water clusters have received a lot of attention as they offer tools to bridge the gap between molecular and bulk scales of water. However, their properties are still not fully understood and deserve further theoretical and experimental investigations. In this work, we simulate the caloric curves of protonated water clusters (H2O)nH+ (n = 20-23). These curves, which have recently been measured experimentally, are characteristic of the phase changes occurring in the aggregates with respect to temperature. The present simulations are achieved by combining parallel-tempering molecular dynamics and the self-consistent-charge density-functional based tight-binding approach and are focused on a restricted size range around (H2O)21H+ which presents singular properties. The shape of the experimental caloric curves and their size dependence are satisfactorily reproduced by the simulations which allows us to further provide a description of the phase transition in terms of structural modifications, dynamics of water molecules and proton mobility. Similar to the experiments, we observe that (H2O)21H+ exhibits a sharper phase transition than the neighbouring size clusters, which can be traced back to both structural and dynamic peculiarities.