Tracing the transition from covalent to non-covalent functionalization of pyrene through C-, N-, and O-based ionic and radical substrates using quantum mechanical calculations

Pyrene is one of the widely investigated aromatic hydrocarbons given its unique optical and electronic properties. Modulating inherent characteristics of pyrene via covalent or non-covalent functionalization has been attractive for a wide variety of advanced biomedical and other device applications. In this study, we have reported the functionalization of pyrene via C, N, and O based ionic and radical substrates, and emphasized the transition of covalent to non-covalent functionalization through making the modulation in the substrate. As expected, strong interactions were observed for cationic substrates, however, anionic substrates also exhibited a competitive binding strength. For instance, methyl and phenyl substituted CH₃ complexes exhibited IEs in the range of -17 kcal mol^-1 to -127 kcal mol^-1 and -14 kcal mol^-1 to -95 kcal mol^-1 and for cationic and anionic substrates, respectively. The analysis of topological parameters showed that un-substituted cationic, anionic, and radical substrates interact with pyrene via covalent interactions, and further become non-covalent upon methylation and phenylation of the substrates. In cationic complexes, the polarisation component is observed to be dominating the interactions, whereas highly competitive contributions from polarization and exchange components were observed in anionic and radical complexes. The contribution of the dispersion component increases with an increase in the degree of methylation and phenylation of the substrate, and starts dominating once the interactions become non-covalent in nature.

Interaction of Aromatic Molecules with Forsterite: Accuracy of the Periodic DFT-D4 Method

Density functional theory (DFT) has provided deep atomic-level insights into the adsorption behavior of aromatic molecules on solid surfaces. However, modeling the surface phenomena of large molecules on mineral surfaces with accurate plane wave methods (PW) can be orders of magnitude more computationally expensive than localized atomic orbitals (LCAO) methods. In the present work, we propose a less costly approach based on the DFT-D4 method (PBE-D4), using LCAO, to study the interactions of aromatic molecules with the {010} forsterite (Mg2SiO4) surface for their relevance in astrochemistry. We studied the interaction of benzene with the pristine {010} forsterite surface and with transition-metal cations (Fe2+ and Ni2+) using PBE-D4 and a vdW-inclusive density functional (Dion, Rydberg, Schröder, Langreth, and Lundqvist (DRSLL)) with LCAO methods. PBE-D4 shows good agreement with coupled-cluster methods (CCSD(T)) for the binding energy trend of cation complexes and with PW methods for the binding energy of benzene on the forsterite surface with a difference of about 0.03 eV. The basis set superposition error (BSSE) correction is shown to be essential to ensure a correct estimation of the binding energies even when large basis sets are employed for single-point calculations of the optimized structures with smaller basis sets. We also studied the interaction of naphthalene and benzocoronene on pristine and transition-metal-doped {010} forsterite surfaces as a test case for PBE-D4. Yielding results that are in good agreement with the plane wave methods with a difference of about 0.02-0.17 eV, the PBE-D4 method is demonstrated to be effective in unraveling the binding structures and the energetic trends of aromatic molecules on pristine and transition-metal-doped forsterite mineral surfaces. Furthermore, PBE-D4 results are in good agreement with its predecessor PBE-D3(BJM) and with the vdW-inclusive density functionals, as long as transition metals are not involved. Hence, PBE-D4/CP-DZP has been proven to be a robust theory level to study the interaction of aromatic molecules on mineral surfaces.

On-surface synthesis of planar acenes via regioselective aryl-aryl coupling

The reaction of 2,2'-dibromo-biphenyl on a Ag(111) surface leads to the formation of planar acenes with a high-regioselectivity rather than nonplanar saddle-shaped tetraphenylene as the typical product in solution chemistry. The regioselective aryl-aryl coupling reaction is attributed to the hydrogen repulsion between the reactants on the confined two-dimensional surface.

Superhydrogenation of pentacene: the reactivity of zigzag-edges

Investigating the hydrogenation of carbonaceous materials is of interest in a wide range of research areas including electronic device development, hydrogen storage, and, in particular, astrocatalytic formation of molecular hydrogen in the universe. Polycyclic Aromatic Hydrocarbons (PAHs) are ubiquitous in space, locking up close to 15% of the elementary carbon. We have used thermal desorption measurements to study the hydrogenation sequence of pentacene from adding one additional H to the fully hydrogenated pentacene species. The experiments reveal that hydrogenated species with an even number of excess H atoms are highly preferred over hydrogenated species with an odd number of H atoms. In addition, the experiments show that specific hydrogenation states of pentacene with 2, 4, 6, 10, 16 and 22 extra H atoms are preferred over other even numbers. We have investigated the structural stability and activation energy barriers for the superhydrogenation of pentacene using Density Functional Theory. The results reveal a preferential hydrogenation pattern set by the activation energy barriers of the hydrogenation steps. Based on these studies, we formulate simple concepts governing the hydrogenation that apply equally well for different PAHs.

Hydrogenated polycyclic aromatic hydrocarbons: isomerism and aromaticity

A simple model based on adjacency matrices is introduced to study the stability of hydrogenated polycyclic aromatic hydrocarbons. Aromaticity governs their relative stability having the most stable isomers the higher number of non-hydrogenated rings.