Molecular dynamics simulations on [FePAH]+ π-complexes of astrophysical interest: anharmonic infrared spectroscopy

Theoretical study of electronic structures, magnetic properties, and ultrafast spin manipulation in transition metal adsorbed polycyclic-aromatic-hydrocarbon molecules

We present a first-principles study of the structural, electronic, and magnetic properties of TM(PAH)0/+ (TM = Fe, Co, Ni; PAH = C10H8, C16H10, C24H12, C32H14) complexes and explore the laser-induced spin dynamics as well as their stability with respect to various laser parameters. For each complex, the most stable configuration shows that the TM atom prefers to adsorb at the hollow site of the carbon ring with a slight deviation from the center. The electronic structure and spin localization of the complexes are found to be largely affected by the TM type. Driven by various laser pulses, spin-crossover scenarios are achieved in all structures, while spin-transfer between TM and PAH is achieved in Ni(C10H8), Ni(C16H10), and Ni(C24H12). The influence of the laser energy and chirp on the dynamics is also investigated, providing important information regarding the stability and sensitivity of the dynamical process. All results are believed to reveal the physics nature of the TM-PAH systems, to guide the experimental realization of their ultrafast spin dynamics and thus to promote their applications in future spintronics.

Modeling the infrared cascade spectra of small PAHs: the 11.2 μm band

The profile of the 11.2 μm feature of the infrared (IR) cascade emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules is investigated using a vibrational anharmonic method. Several factors are found to affect the profile including: the energy of the initially absorbed ultraviolet (UV) photon, the density of vibrational states, the anharmonic nature of the vibrational modes, the relative intensities of the vibrational modes, the rotational temperature of the molecule, and blending with nearby features. Each of these factors is explored independently and influence either the red or blue wing of the 11.2 μm feature. The majority impact solely the red wing, with the only factor altering the blue wing being the rotational temperature.

Infrared Spectroscopy of Chemically Diverse Carbon Clusters: A Data-Driven Approach

Carbon clusters exhibit a broad diversity of topologies and shapes, encompassing fullerene-like cages, graphene-like flakes, and more disordered pretzel-like and branched structures. Here, we examine computationally their infrared spectra in relation with these structures from a statistical perspective. Individual spectra for broad samples of isomers were determined by means of the self-consistent charge density functional-based tight-binding method, and an interpolation scheme is designed to reproduce the spectral features by regression on a much smaller subset of the sample. This interpolation proceeds by encoding the structures using appropriate descriptors and selecting them through principal component analysis, Gaussian regression or inverse distance weighting providing the nonlinear weighting functions. Metric learning is employed to reduce the global error on a preselected testing set. The interpolated spectra satisfactorily reproduce the specific spectral features and their dependence on the size and shape, enabling quantitative prediction away from the testing set. Finally, the classification of structures within the four proposed families is critically discussed through a statistical analysis of the sample based on iterative label spreading.

A Spectroscopic View on Cosmic PAH Emission

ConspectusPolycyclic aromatic hydrocarbon molecules (PAHs) are ubiquitously present at high abundances in the Universe. They are detected through their infrared (IR) fluorescence UV pumped by nearby massive stars. Hence, their infrared emission is used to determine the star formation rate in galaxies, one of the key indicators for understanding the evolution of galaxies. Together with fullerenes, PAHs are the largest molecules found in space. They significantly partake in a variety of physical and chemical processes in space, influencing star and planet formation as well as galaxy evolution.Since the IR features from PAHs originate from chemical bonds involving only nearest neighbor atoms, they have only a weak dependence on the size and structure of the molecule, and it is therefore not possible to identify the individual PAH molecules that make up the cosmic PAH family. This strongly hampers the interpretation of their astronomical fingerprints. Despite the lack of identification, constraints can be set on the characteristics of the cosmic PAH family thanks to a joint effort of astronomers, physicists, and chemists.This Account presents the spectroscopic properties of the cosmic PAH emission as well as the intrinsic spectroscopic properties of PAHs and astronomical modeling of the PAH evolution required for the interpretation of the cosmic PAH characteristics. We discuss the observed spectral signatures tracing PAH properties such as charge, size, and structure and highlight the related challenges. We discuss the recent success of anharmonic calculations of PAH infrared absorption and emission spectra and outline the path forward. Finally, we illustrate the importance of models on PAH processing for the interpretation of the astronomical data in terms of the charge balance and PAH destruction.Throughout this Account, we emphasize that huge progress is on the horizon on the astronomical front. Indeed, the world is eagerly awaiting the launch of the James Webb Space Telescope (JWST). With its incredible improvement in spatial resolution, combined with its complete spectral coverage of the PAH infrared emission bands at medium spectral resolution and superb sensitivity, the JWST will revolutionize PAH research. Previous observations could only present spectra averaged over regions with vastly different properties, thus greatly confusing their interpretation. The amazing spatial resolution of JWST will disentangle these different regions. This will allow us to quantify precisely how PAHs are modified by the physical conditions of their host environment and thus trace how PAHs evolve across space. However, this will only be achieved when the necessary (and still missing) fundamental properties of PAHs, outlined in this Account, are known. We strongly encourage you to join this effort.

Anharmonic Infrared Spectra of Thermally Excited Pyrene (C16H10): A Combined View of DFT-Based GVPT2 with AnharmonicCaOs, and Approximate DFT Molecular dynamics with DemonNano

The study of the Aromatic Infrared Bands (AIBs) in astronomical environments has opened interesting spectroscopic questions on the effect of anharmonicity on the infrared (IR) spectrum of hot polycyclic aromatic hydrocarbons (PAHs) and related species in isolated conditions. The forthcoming James Webb Space Telescope will unveil unprecedented spatial and spectral details in the AIB spectrum; significant advancement is thus necessary now to model the infrared emission of PAHs, their presumed carriers, with enough detail to exploit the information content of the AIBs. This requires including anharmonicity in such models, and to do so systematically for all species included, requiring a difficult compromise between accuracy and efficiency. We performed a benchmark study to compare the performances of two methods in calculating anharmonic spectra, comparing them to available experimental data. One is a full quantum method, AnharmoniCaOs, relying on an ab initio potential, and the other relies on Molecular Dynamics simulations using a Density Functional based Tight Binding potential. The first one is found to be very accurate and detailed, but it becomes computationally very expensive for increasing temperature; the second is faster and can be used for arbitrarily high temperatures, but is less accurate. Still, its results can be used to model the evolution with temperature of isolated bands. We propose a new recipe to model anharmonic AIB emission using minimal assumptions on the general behaviour of band positions and widths with temperature, which can be defined by a small number of empirical parameters. Modelling accuracy will depend critically on these empirical parameters, allowing for an incremental improvement in model results, as better estimates become gradually available.

SCC-DFTB Parameters for Fe-C Interactions

We present an optimized density-functional tight-binding (DFTB) parameterization for iron-based complexes based on the popular trans3d set of parameters. The transferability of the original and optimized parameterizations is assessed using a set of 50 iron complexes, which include carbonyl, cyanide, polypyridine, and cyclometalated ligands. DFTB-optimized structures predicted using the trans3d parameters show a good agreement with both experimental crystal geometries and density functional theory (DFT)-optimized structures for Fe-N bond lengths. Conversely, Fe-C bond lengths are systematically overestimated. We improve the accuracy of Fe-C interactions by truncating the Fe-O repulsive potential and reparameterizing the Fe-C repulsive potential using a training set of six isolated iron complexes. The new trans3d*-LANLFeC parameter set can produce accurate Fe-C bond lengths in both geometry optimizations and molecular dynamics (MD) simulations, without significantly affecting the accuracy of Fe-N bond lengths. Moreover, the potential energy curves of Fe-C interactions are considerably improved. This improved parameterization may open the door to accurate MD simulations at the DFTB level of theory for large systems containing iron complexes, such as sensitizer-semiconductor assemblies in dye-sensitized solar cells, that are not easily accessible with DFT approaches because of the large number of atoms.

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.

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.

Experimental Approach to the Study of Anharmonicity in the Infrared Spectrum of Pyrene from 14 to 723 K

Quantifying the effect of anharmonicity on the infrared spectrum of large molecules such as polycyclic aromatic hydrocarbons (PAHs) at high temperatures is the focus of a number of theoretical and experimental studies, many of them motivated by astrophysical applications. We recorded the IR spectrum of pyrene C16H10 microcrystals embedded in KBr pellets over a wide range of temperatures (14-723 K) and studied the evolution of band positions, widths, and integrated intensities with temperature. We identified jumps for some of the spectral characteristics of some bands in the 423-473 K range. These were attributed to a change of phase from crystal to molten in condensed pyrene, which appears to affect more strongly bands involving large CH motions. Empirical anharmonicity factors that quantify the linear evolution of band positions and widths with temperature for values larger than ∼150-250 K, depending on the band, were retrieved from both phases and averaged to provide recommended values for these anharmonicity factors. The derived values were found to be consistent with available gas phase data. We conclude about the relevance of the methodology to produce data that can be compared with calculated anharmonic IR spectra and provide input for models that simulate the IR emission of astro-PAHs.

Formation of coronene:water complexes: FTIR study in argon matrices and theoretical characterisation

Coronene preferentially forms σ-type rather than π-type complexes with water in low temperature argon matrices.

Cationic Methylene-Pyrene Isomers and Isomerization Pathways: Finite Temperature Theoretical Studies

This paper provides spectral characterizations of the two isomers of the 1-methylenepyrene cation, namely, the 1-pyrenemethylium and a pyrene-like isomer owing a tropylium cycle. Both are possible photodissociation products of the 1-methylpyrene cation and were proposed as potential contributors to the diffuse interstellar bands. In that respect, vibrational and electronic spectra are computed for the optimized structures at the density functional theory (DFT) and time-dependent (TD-)DFT levels. Finite temperature effects on these spectra are estimated from molecular dynamics simulations within the density functional-based tight-binding (DFTB) and TD-DFTB frameworks, these methods being first benchmarked against DFT and TD-DFT calculations. The computed spectra allow discrimination of the two isomers. When the temperature increases, bands are observed to redshift and merge. The isomerization mechanism is investigated with the metadynamics technique, a biased dynamics scheme allowing to probe reaction mechanisms with high energy barriers by investigating the free energy surface at various temperatures. Four pathways with similar barrier heights (3.5-4 eV) are found, showing that the interconversion process would only occur in interstellar clouds under photoactivation. The present study opens the way to simulations on larger methyl- and methylenePAHs of astrophysical interest and their experimental investigation.

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.

Photoelectron and electron momentum spectroscopy of tetrahydrofuran from a molecular dynamical perspective

The results of experimental studies of the valence electronic structure of tetrahydrofuran employing He I photoelectron spectroscopy as well as Electron Momentum Spectroscopy (EMS) have been reinterpreted on the basis of Molecular Dynamical simulations employing the classical MM3 force field and large-scale quantum mechanical simulations employing Born-Oppenheimer Molecular Dynamics in conjunction with the dispersion corrected ωB97XD exchange-correlation functional. Analysis of the produced atomic trajectories demonstrates the importance of thermal deviations from the lowest energy path for pseudorotation, in the form of considerable variations of the ring-puckering amplitude. These deviations are found to have a significant influence on several outer-valence electron momentum distributions, as well as on the He I photoelectron spectrum.

Vibrational spectroscopy and molecular dynamics of water monomers and dimers adsorbed on polycyclic aromatic hydrocarbons

This paper reports structures, energetics, dynamics and spectroscopy of H2O and (H2O)2 systems adsorbed on coronene (C24H12), a compact polycyclic aromatic hydrocarbon (PAH). On-the-fly Born-Oppenheimer molecular dynamics simulations are performed for temperatures T varying from 10 to 300 K, on a potential energy surface obtained within the self-consistent-charge density-functional based tight-binding (SCC-DFTB) approach. Anharmonic infrared (IR) spectra are extracted from these simulations. We first benchmark the SCC-DFTB semi-empirical hamiltonian vs. DFT (Density Functional Theory) calculations that include dispersion, on (C6H6)(H2O)1,2 small complexes. We find that charge corrections and inclusion of dispersion contributions in DFTB are necessary to obtain consistent structures, energetics and IR spectra. Using this Hamiltonian, the structures, energetics and IR features of the low-energy isomers of (C24H12)(H2O)1,2 are found to be similar to the DFT ones, with evidence for a stabilizing edge-coordination. The temperature dependence of the motions of H2O and (H2O)2 on the surface of C24H12 is analysed, revealing ultra-fast periodic motion. The water dimer starts diffusing at a higher temperature than the water monomer (150 K vs. 10 K respectively), which appears to be consistent with the binding energies. Qualitative and quantitative analyses of the effects of T on the IR spectra are performed. Anharmonic factors in particular are derived and it is shown that they can be used as signatures for the presence of PAH-water complexes. Finally, this paper lays the foundations for the studies of larger (PAH)m(H2O)n clusters, that can be treated with the efficient computational approach benchmarked in this paper.