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.

Influence of the internal degrees of freedom of coronene molecules on the nonlinear dynamics of a columnar chain

The nonlinear dynamics of a one-dimensional molecular crystal in the form of a chain of planar coronene molecules is analyzed. Using molecular dynamics, it is shown that a chain of coronene molecules supports acoustic solitons, rotobreathers, and discrete breathers. An increase in the size of planar molecules in a chain leads to an increase in the number of internal degrees of freedom. This results in an increase in the rate of emission of phonons from spatially localized nonlinear excitations and a decrease in their lifetime. Presented results contribute to the understanding of the effect of the rotational and internal vibrational modes of molecules on the nonlinear dynamics of molecular crystals.

Theoretical and computational methodologies for understanding coordination self-assembly complexes

This perspective highlights three theoretical and computational methods to capture the coordination self-assembly processes at the molecular level: quantum chemical modeling, molecular dynamics, and reaction network analysis. These methods cover the different scales from the metal-ligand bond to a more global aspect, and approaches that are best suited to understand the coordination self-assembly from different perspectives are introduced. Theoretical and numerical researches based on these methods are not merely ways of interpreting the experimental studies but complementary to them.

Towards kinetic control of coordination self-assembly: a case study of a Pd3L6 double-walled triangle to predict the outcomes by a reaction network model

Numerical analysis of self-assembly process (NASAP) was performed for a [Pd3L6]6+ double-walled triangle (DWT) complex. With a chemical reaction network and a parameter set of the reaction rate constants obtained from a numerical search in an eighteen-dimensional parameter space to obtain a good fit to the data from the experimental counterpart (quantitative analysis of self-assembly process, QASAP), a refined calculation resulted in a detailed time evolution of each molecular species. Analysis based on those clues revealed dominant self-assembly pathways and a balance between inter- and intramolecular reactions, and enabled prediction of the reaction outcomes depending on the initial stoichiometric ratio under kinetic control.

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.

Solvation of coronene oligomers by para-H2 molecules: the effects of size and shape

The stepwise solvation of various cationic coronene oligomers by para-hydrogen (p-H₂) molecules was computationally investigated using a united-atom model for the p-H₂ molecules and the Silvera-Goldman potential, together with a polarizable description for the interaction with the hydrocarbon molecules. A survey of the energy landscape for oligomers containing between 1 and 4 coronene molecules and possible different conformers was carried out using standard global optimization, the hydrocarbon complex being kept as rigid. The most stable structures provided the starting configuration of systematic path-integral molecular dynamics simulations at 2 K. The variations of the geometric and energetic properties of the solvation shell were determined with increasing number of para-hydrogen molecules. The relative stability of the solvation shell is generally found to be more robustly determined by the energy increment (or dissociation energy) than by geometrical indicators, especially when the oligomers have less ordered structures. In agreement with recent mass spectrometry experiments, the size at which the first solvation shell is complete is found to vary approximately linearly with the oligomer size when the coronene molecules stack together, with a slope that is related to the offset between two successive molecules.

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.

Sphere Encapsulated Monte Carlo: Obtaining Minimum Energy Configurations of Large Aromatic Systems

We introduce a simple global optimization approach that is able to find minimum energy configurations of clusters containing aromatic molecules. The translational and rotational perturbations required in Monte Carlo-based methods often lead to unrealistic configurations within which two or more molecular rings intersect, causing many of the computational steps to be rejected and the optimization process to be inefficient. Here we develop a modification of the basin-hopping global optimization procedure tailored to tackle problems with intersecting molecular rings. Termed the Sphere Encapsulated Monte Carlo (SEMC) method, this method introduces sphere-based rearrangement and minimization steps at each iteration, and its performance is shown through the exploration of potential energy landscapes of polycyclic aromatic hydrocarbon (PAH) clusters, systems of interest in combustion and astrophysics research. The SEMC method provides clusters that are accurate to 5% mean difference of the minimum energy at a 10-fold speed up compared to previous work using advanced molecular dynamics simulations. Importantly, the SEMC method captures key structural characteristics and molecular size partitioning trends as measured by the molecular radial distances and coordination numbers. The advantages of the SEMC method are further highlighted in its application to previously unstudied heterogeneous PAH clusters.

Aggregation Behavior of Model Asphaltenes Revealed from Large-Scale Coarse-Grained Molecular Simulations

Fully atomistic simulations of models of asphaltenes in simple solvents have allowed the study of trends in aggregation phenomena to understand the underlying role played by molecular structure. The detail included at this scale of molecular modeling is, however, at odds with the required spatial and temporal resolution needed to fully understand asphaltene aggregation. The computational cost required to explore the relevant scales can be reduced by employing coarse-grained (CG) models, which consist of lumping a few atoms into a single segment that is characterized by effective interactions. In this work, CG force fields developed via the statistical associating fluid theory (SAFT-γ) [ Müller , E. A. ; Jackson , G. Annu. Rev. Chem. Biomol. Eng. 5 , 2014 , 405 - 427 ] equation of state (EoS) provide a reliable pathway to link the molecular description with macroscopic thermophysical data. A recent modification of the SAFT-VR EoS [ Müller , E. A. ; Mejía , A. Langmuir 33 , 2017 , 11518 - 11529 ], which allows for the parameterization of homonuclear rings, is selected as the starting point to develop CG models for polycyclic aromatic hydrocarbons. The new aromatic-core models, along with others published for simpler organic molecules, are adopted for the construction of asphaltene models by combining different chemical moieties in a group-contribution fashion. We apply the procedure to two previously reported asphaltene models and perform molecular dynamics simulations to validate the coarse-grained representation against benchmark systems of 27 asphaltenes in a pure solvent (toluene or heptane) described in a fully atomistic fashion. An excellent match between both levels of description is observed for the cluster size, radii of gyration, and relative-shape-anisotropy-factor distributions. We exploit the advantages of the CG representation by simulating systems containing up to 2000 asphaltene molecules in an explicit solvent investigating the effect of asphaltene concentration, solvent composition, and temperature on aggregation. By studying large systems facilitated by the use of CG models, we observe stable continuous distributions of molecular aggregates at conditions away from the two-phase precipitation point. As a further example application, a widely accepted interpretation of cluster-size distributions in asphaltenic systems is challenged by performing system-size tests, reversibility checks, and a time-dependence analysis. The proposed coarse-graining procedure is seen to be general and predictive and, hence, can be applied to other asphaltenic molecular structures.

Coarse-grained modeling of the nucleation of polycyclic aromatic hydrocarbons into soot precursors

The aggregation and physical growth of polycyclic aromatic hydrocarbon (PAH) molecules was simulated using a coarse-grained (CG) approach based on the Paramonov-Yaliraki (PY) potential and a stochastic Monte Carlo framework, following earlier efforts in which the structure [Phys. Chem. Chem. Phys., 2016, 18, 13736] and equilibrium thermodynamics [Phys. Chem. Chem. Phys., 2017, 19, 1884] were investigated and critically compared to the predictions of all-atom models. Homomolecular and heteromolecular assemblies of pyrene, coronene, and circumcoronene were considered at various temperatures and compositions, and the distributions of aggregation products were characterized. Under the simulated conditions, and in agreement with earlier studies, the clusters are rather small and, in the case of pyrene-rich systems, only formed below 1000 K. The clusters obtained by spontaneous aggregation of isolated molecules are statistically analysed. For the selected sizes of tetramers and octamers, broad distributions of isomers are obtained with a clear entropic stabilization. In heteronuclear assemblies, our results suggest a minor spontaneous segregation towards pure and equi concentrations at variance with purely statistical expectations.

Structural Correlations and Percolation in Twisted Perylene Diimides Using a Simple Anisotropic Coarse-Grained Model

Large, twisted, and fused conjugated molecular architectures have begun to appear more prominently in the organic semiconductor literature. From a modeling perspective, such structures present a challenge to conventional simulation techniques; atomistic resolutions are computationally inefficient, while traditional isotropic coarse-grained models do not capture the inherent anisotropies of the molecules. In this work, we develop a simple coarse-grained model that explicitly incorporates the anisotropy of these molecular architectures, thereby providing a route toward analyzing π-stacking, and thus qualitative electronic structure, at a computationally efficient coarse-grained resolution. Our simple coarse-grained model maintains relative orientations of conjugated rings, as well as inter-ring dihedrals, that are critical for understanding electronic and excitonic transport in bulk systems. We apply this model to understand structural correlations in several recently synthesized perylene diimide (PDI)-based organic semiconductors. Twisted and nonplanar molecular architectures are found to promote amorphous morphologies while maintaining local π-stacking. A graph theoretical network analysis demonstrates that these twisted molecules are more likely to form percolating three-dimensional pathways for charge motion than strictly planar molecules, which show connectivity in only one dimension.

The structure of coronene cluster ions inferred from H2 uptake in the gas phase

Mass spectra of helium nanodroplets doped with H₂ and coronene feature anomalies in the ion abundance that reveal anomalies in the energetics of adsorption sites. The coronene monomer ion strongly adsorbs up to n = 38 H₂ molecules indicating a commensurate solvation shell that preserves the D_6h symmetry of the substrate. No such feature is seen in the abundance of the coronene dimer through tetramer complexed with H₂; this observation rules out a vertical columnar structure. Instead we see evidence for a columnar structure in which adjacent coronenes are displaced in parallel, forming terraces that offer additional strong adsorption sites. The experimental value for the number of adsorption sites per terrace, approximately six, barely depends on the number of coronene molecules. The displacement estimated from this number exceeds the value reported in several theoretical studies of the bare, neutral coronene dimer.