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

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.

Accurate SCC-DFTB Parametrization of Liquid Water with Improved Atomic Charges and Iterative Boltzmann Inversion

This work presents improvements of the description of liquid water within the self-consistent-charge density-functional based tight-binding scheme combining the use of Weighted Mulliken (WMull) charges and optimized O-H repulsive potential through the iterative Boltzmann inversion (IBI) process. The quality of the newly developed models is validated considering pair radial distribution functions (RDFs), as well as other structural, energetic, thermodynamic, and dynamic properties. The use of WMull charges certainly improves the agreement with experimental data, however leading to over-structured RDFs at short distance, that can be further improved by considering an optimized O-H repulsive potential obtained by the IBI process. Three different schemes were used to optimize this potential: (i) optimization including short O-H distances. This led to accurate RDFs as well as improved self-diffusion coefficient and heat of vaporization, while the proton transfer energy barrier is severely deteriorated; (ii) optimization starting at long distance. The proton transfer energy barrier is recovered while the heat of vaporization is deteriorated and the O-H RDF is less accurate at short distance; (iii) optimization within the path-integral molecular dynamics scheme which allows us to exclude nuclear quantum effects from the repulsive potential. The latter potential, in conjunction with the WMull improved atomic charges, provides similar results as (i) for structural, dynamic, and thermodynamic properties while recovering a large part of the proton transfer energy barrier. It therefore offers a good compromise to study both dynamic properties and chemistry within liquid water at a quantum chemical level.

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.

How does tuning the van der Waals bonding strength affect adsorbate structure?

Organic molecular thin-films are employed for manufacturing a wide variety of electronic devices, including memory devices and transistors. A precise description of the atomic-scale interactions in aromatic carbon systems is of paramount importance for the design of organic thin-films and carbon-based nanomaterials. Here we investigate the binding and structure of pyrazine on graphite with neutron diffraction and spin-echo measurements. Diffraction data of the ordered phase of deuterated pyrazine, (C₄D₄N₂), adsorbed on the graphite (0001) basal plane surface are compared to scattering simulations and complemented by van der Waals corrected density functional theory calculations. The lattice constant of pyrazine on graphite is found to be (6.06 ± 0.02) Å. Compared to benzene (C₆D₆) adsorption on graphite, the pyrazine overlayer appears to be much more thermodynamically stable, up to 320 K, and continues in layer-by-layer growth. Both findings suggest a direct correlation between the intensity of van der Waals bonding and the stability of the self-assembled overlayer because the nitrogen atoms in the six-membered ring of pyrazine increase the van der Waals bonding in comparison to benzene, which only contains carbon atoms.

A combined theoretical and experimental study of small anthracene-water clusters

Water-cluster interactions with polycyclic aromatic hydrocarbons (PAHs) are of paramount interest in many chemical and biological processes. We report a study of anthracene monomers and dimers with water (up to four)-cluster systems utilizing molecular beam vacuum-UV photoionization mass spectrometry and density functional calculations. Structural loss in photoionization efficiency curves when adding water indicates that various isomers are generated, while theory indicates only a slight shift in energy in photoionization states of different isomers. Calculations reveal that the energetic tendency of water is to remain clustered and not to disperse around the PAH. Theoretically, we observe water confinement exclusively in the case of four water clusters and only when the anthracenes are in a cross configuration due to optimal OH⋯π interactions, indicating dependence on the size and structure of the PAH. Furthermore theory sheds light on the structural changes that occur in water upon ionization of anthracene, due to the optimal interactions of the resulting hole and water hydrogen atoms.

Periodic DFTB for Supported Clusters: Implementation and Application on Benzene Dimers Deposited on Graphene

The interest for properties of clusters deposited on surfaces has grown in recent years. In this framework, the Density Functional based Tight Binding (DFTB) method appears as a promising tool due to its ability to treat extended systems at the quantum level with a low computational cost. We report the implementation of periodic boundary conditions for DFTB within the deMonNano code with k-points formalism and corrections for intermolecular interactions. The quality of DFTB results is evaluated by comparison with dispersion-corrected DFT calculations. Optimized lattice properties for a graphene sheet and graphite bulk are in agreement with reference data. The deposition of both benzene monomer and dimers on graphene are investigated and the observed trends are similar at the DFT and DFTB levels. Moreover, interaction energies are of similar orders of magnitude for these two levels of calculation. This study has evidenced the high stability of a structure made of two benzene molecules deposited close to each other on the graphene sheet. This work demonstrates the ability of the new implementation to investigate surface-deposited molecular clusters properties.

Using Excimeric Fluorescence to Study How the Cooling Rate Determines the Behavior of Naphthalenes in Freeze-Concentrated Solutions: Vitrification and Crystallization

We utilized fluorescence spectroscopy to learn about the molecular arrangement of naphthalene (Np) and 1-methylnaphthalene (MeNp) in frozen aqueous solutions. The freezing induces pronounced compound aggregation in the freeze-concentrated solution (FCS) in between the ice grains. The fluorescence spectroscopy revealed prevalent formation of a vitrified solution and minor crystallization of aromatic compounds. The FCS is shown as a specific environment, differing significantly from not only the pure compounds but also the ice surfaces. The results indicate marked disparity between the behavior of the Np and the MeNp; the cooling rate has a major impact on the former but not on the latter. The spectrum of the Np solution frozen at a faster cooling rate (ca 20 K/min) exhibited a temperature-dependent spectral behavior, whereas the spectrum of the solution frozen at a slower rate (ca 2 K/min) did not alter before melting. We interpret the observation through considering the varied composition of the FCS: Fast freezing leads to a higher water content expressed by the plasticizing effect, allowing molecular rearrangement, while slow cooling produces a more concentrated and drier environment. The experiments were conceived as generalizable for environmentally relevant pollutants and human-made freezing.

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.

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.

PAH/PAH(CF3 )n Donor/Acceptor Charge-Transfer Complexes in Solution and in Solid-State Co-Crystals

A solution, solid-state, and computational study is reported of polycyclic aromatic hydrocarbon PAH/PAH(CF3 )n donor/acceptor (D/A) charge-transfer complexes that involve six PAH(CF3 )n acceptors with known gas-phase electron affinities that range from 2.11(2) to 2.805(15) eV and four PAH donors, including seven CT co-crystal X-ray structures that exhibit hexagonal arrays of mixed π-stacks with 1/1, 1/2, or 2/1 D/A stoichiometries (PAH=anthracene, azulene, coronene, perylene, pyrene, triphenylene; n=5, 6). These are the first D/A CT complexes with PAH(CF3 )n acceptors to be studied in detail. The nine D/A combinations were chosen to allow several structural and electronic comparisons to be made, providing new insights about controlling D/A interactions and the structures of CT co-crystals. The comparisons include, among others, CT complexes of the same PAH(CF3 )n acceptor with four PAH donors and CT complexes of the same donor with four PAH(CF3 )n acceptors. All nine CT complexes exhibit charge-transfer bands in solution with λmax between 467 and 600 nm. A plot of E(λmax ) versus [IE(donor)-EA(acceptor)] for the nine CT complexes studied is linear with a slope of 0.72±0.03 eV eV-1 . This plot is the first of its kind for CT complexes with structurally related donors and acceptors for which precise experimental gas-phase IEs and EAs are known. It demonstrates that conclusions based on the common assumption that the slope of a CT E(λmax ) versus [IE-EA] plot is unity may be incorrect in at least some cases and should be reconsidered.