Catalysis on Pristine 2D Materials via Dispersion and Electrostatic Interactions

A semi-automated quantum-mechanical workflow for the generation of molecular monolayers and aggregates

Organic electronics (OE) such as organic light-emitting diodes or organic solar cells represent an important and innovative research area to achieve global goals like environmentally friendly energy production. To accelerate OE material discovery, various computational methods are employed. For the initial generation of structures, a molecular cluster approach is employed. Here, we present a semi-automated workflow for the generation of monolayers and aggregates using the GFNn-xTB methods and composite density functional theory (DFT-3c). Furthermore, we present the novel D11A8MERO dye interaction energy benchmark with high-level coupled cluster reference interaction energies for the assessment of efficient quantum chemical and force-field methods. GFN2-xTB performs similar to low-cost DFT, reaching DFT/mGGA accuracy at two orders of magnitude lower computational cost. As an example application, we investigate the influence of the dye aggregate size on the optical and electrical properties and show that at least four molecules in a cluster model are needed for a qualitatively reasonable description.

Catalyzing epoxy oxygen migration on the basal surface of graphene oxide using strong hydrogen-bond donors

High-level double-hybrid DFT simulations reveal that strong hydrogen-bond-donor catalysts (e.g., ethylene glycol, guanidine, and thiourea) significantly accelerate the migration of epoxy oxygen on the surface of graphene oxide, enhancing the reaction rate by 6-12 orders of magnitude. These results shed light on previously puzzling experimental observations.

Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

London dispersion (LD) interactions are the main contribution of the attractive part of the van der Waals potential. Even though LD effects are the driving force for molecular aggregation and recognition, the role of these omnipresent interactions in structure and reactivity had been largely underappreciated over decades. However, in the recent years considerable efforts have been made to thoroughly study LD interactions and their potential as a chemical design element for structures and catalysis. This was made possible through a fruitful interplay of theory and experiment. This review highlights recent results and advances in utilizing LD interactions as a structural motif to understand and utilize intra- and intermolecularly LD-stabilized systems. Additionally, we focus on the quantification of LD interactions and their fundamental role in chemical reactions.

Mutual induced fit transition structure stabilization of corannulene's bowl-to-bowl inversion in a perylene bisimide cyclophane

Corannulene is known to undergo a fast bowl-to-bowl inversion at r.t. via a planar transition structure (TS). Herein we present the catalysis of this process within a perylene bisimide (PBI) cyclophane composed of chirally twisted, non-planar chromophores, linked by para-xylylene spacers. Variable temperature NMR studies reveal that the bowl-to-bowl inversion is significantly accelerated within the cyclophane template despite the structural non-complementarity between the binding site of the host and the TS of the guest. The observed acceleration corresponds to a decrease in the bowl-to-bowl inversion barrier of 11.6 kJ mol-1 compared to the uncatalyzed process. Comparative binding studies for corannulene (20 π-electrons) and other planar polycyclic aromatic hydrocarbons (PAHs) with 14 to 24 π-electrons were applied to rationalize this barrier reduction. They revealed high binding constants that reach, in tetrachloromethane as a solvent, the picomolar range for the largest guest coronene. Computational models corroborate these experimental results and suggest that both TS stabilization and ground state destabilization contribute to the observed catalytic effect. Hereby, we find a "mutual induced fit" between host and guest in the TS complex, such that mutual geometric adaptation of the energetically favored planar TS and curved π-systems of the host results in an unprecedented non-planar TS of corannulene. Concomitant partial planarization of the PBI units optimizes noncovalent TS stabilization by π-π stacking interactions. This observation of a "mutual induced fit" in the TS of a host-guest complex was further validated experimentally by single crystal X-ray analysis of a host-guest complex with coronene as a qualitative transition state analogue.

Catalytic effect of graphene on the inversion of corannulene using a continuum approach with the Lennard-Jones potential

The catalytic effect of graphene on the corannulene bowl-to-bowl inversion is confirmed in this paper using a pair-wise dispersion interaction model. In particular, a continuum approach together with the Lennard-Jones potential are adopted to determine the interaction energy between corannulene and graphene. These results are consistent with previous quantum chemical studies, which showed that a graphene sheet reduces the barrier height for the bowl-to-bowl inversion in corannulene. However, the results presented here demonstrate, for the first time, that the catalytic activity of graphene can be reproduced using pair-wise dispersion interactions alone. This demonstrates the major role that pair-wise dispersion interactions play in the catalytic activity of graphene.

Noncovalent wedging effect catalyzed the cis to syn transformation of a surface-adsorbed polymer backbone toward an unusual thermodynamically stable supramolecular product

The significant influence of noncovalent interactions on catalytic processes has been recently appreciated but is still in its infancy. In this report, it is found that wedging Me-PTCDI (small-molecule) between the alkyl chains of PffBT4T-2OD (polymer) and a graphite substrate can reduce the energy barrier of flipping over the surface-adsorbed alkylthiophene group from the cis to syn conformation, revealing the catalytic role of Me-PTCDI via a noncovalent wedging effect. The wedging of Me-PTCDI brings the interactions between the alkyl chains and substrate to a very weak level by lifting up the alkyl chains, which eliminates the major hindrance of the flipping process to one main factor: the torsion of the dihedral angles of the thiophene group. The Me-PTCDI/syn PffBT4T-2OD arrangement shows unusual stability compared to the cis one because the syn conformation allows the alkyl chains to construct dense lamella and facilitates interactions between Me-PTCDI and the syn PffBT4T-2OD backbones. The results are helpful for boosting the development of noncovalent catalysis and bottom-up fabrications toward devices functionalized at a molecular level.

Modeling the Adsorption of Polycyclic Aromatic Hydrocarbons on Graphynes: An Improved Lennard-Jones Formulation

Research on the development of theoretical methodologies for modeling noncovalent interactions governing the adsorption of polycyclic aromatic hydrocarbons (PAHs) on graphene and other two-dimensional materials is being intensely pursued in recent times. Highly accurate empirical potentials have emerged as a viable alternative to first-principles calculations for performing large-scale simulations. Herein, we report exploration of the potential energy surfaces for the adsorption of cata-condensed and peri-condensed PAHs on graphynes (GYs) using the improved Lennard-Jones (ILJ) potential. Initially, the ILJ potential is parametrized against benchmark electronic structure calculations performed on a selected set of PAH-GY complexes using dispersion-corrected density functional theory. The accuracy of the parametrization scheme is then assessed by a comparison of the adsorption features predicted from the ILJ potential with those computed using electronic structure calculations. The potential energy profiles as well as the single point energy calculations and geometry reoptimizations performed on the minimum-energy configurations predicted by the ILJ potential for a broader range of PAH-GY complexes provided a validation of the parametrization scheme. Finally, by an extrapolation of the PAH adsorption energies on various GYs, we estimated the interlayer cohesion energies for the van der Waals bilayer heterostructures of GYs with graphene to be in the range of 25-50 meV/atom.

An anisotropic dressed pairwise potential model for the adsorption of noble gases on boron nitride sheets

Development of empirical potentials with accurate parameterization is indispensable while modeling large-scale systems. Herein, we report accurate parameterization of an anisotropic dressed pairwise potential model (PPM) for probing the adsorption of noble gases, He, Ne, Ar and Kr on boron nitride sheets. For the noble gas binding on B₄₈N₄₈H₂₄, we carried out a least-squares fit analysis of the dispersion and dispersionless contributions of the interaction potential separately. The transferability of the parameters for a range of molecular model systems of boron nitride is further established. The dressed PPM is then used in conjunction with a global optimization technique, namely particle swarm optimization (PSO) to assess the possibility of performing large-scale simulations with the PPM-PSO methodology. The results obtained for the adsorption of 2-5 noble gases on BN sheets establish the proof-of-concept, encouraging the pursuit of large-scale simulations using the PPM-PSO approach.

Graphene-induced planarization of cyclooctatetraene derivatives

Stable equilibrium compounds containing a planar antiaromatic cyclooctatetraene (COT) ring are promising candidates for organic electronic devices such as organic semiconductor transistors. The planarization of COT by incorporation into rigid planar π-systems, as well as by oxidation or reduction has attracted considerable attention in recent years. Using dispersion-corrected density functional theory calculations, we explore an alternative approach of planarizing COT derivatives by adsorption onto graphene. We show that strong π-π stacking interactions between graphene and COT derivatives induce a planar structure with an antiaromatic central COT ring. In addition to being reversible, this strategy provides a novel approach for planarizing COT without the need for incorporation into a rigid structure, atomic substitution, oxidation, or reduction.

π-π Catalysis in Carbon Flatland-Flipping [8]Annulene on Graphene

Noncovalent interactions are an integral part of the modern catalysis toolbox. Although stronger noncovalent interactions such as hydrogen bonding are commonly the main driving force of catalysis, π-π interactions typically provide smaller additional stabilizations, for example, to afford selectivity enhancements. Here, it is shown computationally that pristine graphene flakes may efficiently catalyze the skeletal inversions of various benzannulated cyclooctatetraene derivatives, providing an example of a catalytic process driven solely by π-π stacking interactions. Hereby, the catalytic effect results from disproportionate shape complementarity between catalyst and transition structure compared with catalyst and reactant. An energy decomposition analysis reveals electrostatic and, especially with increasing system size, to a larger extent, dispersion interactions as the origin of stabilization.

Flat corannulene: when a transition state becomes a stable molecule

Flat corannulene has been considered so far only as a transition state of the bowl-to-bowl inversion process. This study was driven by the prediction that substituents with strong steric repulsion could destabilize the bowl-shaped conformation of this molecule to such an extent that the highly unstable planar geometry would become an isolable molecule. To examine the substituents' effect on the corannulene bowl depth, optimized structures for the highly-congested decakis(t-butylsulfido)corannulene were calculated. The computations, performed with both the M06-2X/def2-TZVP and the B3LYP/def2-TZVP methods (the latter with and without Grimme's D3 dispersion correction), predict that this molecule can achieve two minimum structures: a flat carbon framework and a bowl-shaped structure, which are very close in energy. This rather unusual compound was easily synthesized from decachlorocorannulene under mild reaction conditions, and X-ray crystallographic studies gave similar results to the theoretical predictions. This compound crystallized in two different polymorphs, one exhibiting a completely flat corannulene core and the other having a bowl-shaped conformation.