Coronene Dimer Research Articles

Cooperative Binding and Chirogenesis in an Expanded Perylene Bisimide Cyclophane.

The encapsulation of more than one guest molecule into a synthetic cavity is a highly desirable yet a highly challenging task to achieve for neutral supramolecular hosts in organic media. Herein, we report a neutral perylene bisimide cyclophane, which has a tailored chiral cavity with an interchromophoric distance of 11.2 Å, capable of binding two aromatic guests in a π-stacked fashion. Detailed host-guest binding studies with a series of aromatic guests revealed that the encapsulation of the second guest in this cyclophane is notably more favored than the first one. Accordingly, for the encapsulation of the coronene dimer, a cooperativity factor (α) as high as 485 was observed, which is remarkably high for neutral host-guest systems. Furthermore, a successful chirality transfer, from the chiral host to encapsulated coronenes, resulted in a chiral charge-transfer (CT) complex and the rare observation of circularly polarized emission originating from the CT state for a noncovalent donor-acceptor assembly in solution. The involvement of the CT state also afforded an enhancement in the luminescence dissymmetry factor (glum) value due to its relatively large magnetic transition dipole moment. The 1:2 binding pattern and chirality-transfer were unambiguously verified by single-crystal X-ray diffraction analysis of the host-guest superstructures.

A theoretical exploration of different π-π stacking dimers of coronenes and its substituted analogues

Quantum chemical calculations employing density functional theory have been performed to study the geometry, π-interactions and band gaps of coronene and its mono heteroatom substituted analogues. With the objective of systematic structural modification of coronene for the design of functional materials, the study focuses on the role of substitution of carbon in coronene with heteroatom boron or nitrogen on their interplanar π-interactions, band gaps and stabilization energies. The work intensively explores various π-π stacking interactions in dimers of coronenes, theoretically, paving the possible existence of different polymorphs arising from the dimers. The monomeric planes are stacked together through π-interactions with an interplanar distance of ∼ 3.3 Å to form dimers. A total of 36 dimers of substituted coronenes are formed in addition to their corresponding mirror images. Based on the geometries, the dimers are found to exist in five different types. The calculations show that the band gap of unsubstituted monomer coronene is 3.27 eV and upon substitution of a single carbon atom of coronene with heteroatom nitrogen or boron, the band gaps of monosubstituted coronenes are reduced to 0.81–1.16 eV and are dependent on the position of substitution. The band gaps of the dimers are further reduced compared to those of monomers attributing to π-π stacking interactions. The reduction is more in nitrogen-substituted coronene dimers. Also, the substituted coronene dimers are stabilized by as high as 115 times than those of unsubstituted coronene dimers.

Effects of Iron Addition on the Collision of Polycyclic Aromatic Hydrocarbon Clusters: A Molecular Dynamics Study

In this work, soot particle size distributions in iron-doped premixed ethylene flames are examined using scanning mobility particle sizer measurements. It is found that iron addition promotes the growth in soot particle size, and the enhanced particle coagulation is inferred to be an important reason. To support that, the influence of iron addition on the coagulation of polycyclic aromatic hydrocarbon (PAH) clusters, the analogue of incipient soot particles, is further investigated using molecular dynamics simulations. Based on the results of hundreds of binary head-on collision simulations, the collision between two coronene-Fe-coronene dimers is found to have a significantly higher coagulation efficiency than that between two coronene dimers. However, this enhancement effect weakens as the size of the PAH monomer increases. Although the coagulation efficiency can be increased by iron addition, the collision frequency is almost unaffected, as revealed from the binary off-central collision simulations. Moreover, the simulation results of coronene cluster growth via coagulation show that iron addition promotes coronene cluster growth, leading to larger cluster size, which may explain the larger soot particle size observed in iron-doped flames.

Dimerization of dehydrogenated polycyclic aromatic hydrocarbons on graphene.

Dimerization of polycyclic aromatic hydrocarbons (PAHs) is an important, yet poorly understood, step in the on-surface synthesis of graphene (nanoribbon), soot formation, and growth of carbonaceous dust grains in the interstellar medium (ISM). The on-surface synthesis of graphene and the growth of carbonaceous dust grains in the ISM require the chemical dimerization in which chemical bonds are formed between PAH monomers. An accurate and cheap method of exploring structure rearrangements is needed to reveal the mechanism of chemical dimerization on surfaces. This work has investigated the chemical dimerization of two dehydrogenated PAHs (coronene and pentacene) on graphene via an evolutionary algorithm augmented by machine learning surrogate potentials and a set of customized structure operators. Different dimer structures on surfaces have been successfully located by our structure search methods. Their binding energies are within the experimental errors of temperature programmed desorption measurements. The mechanism of coronene dimer formation on graphene is further studied and discussed.

Mirror‐Image Cofacial Coronene Dimers Characterized by CD and CPL Spectroscopy: A Twisted Bilayer Nanographene

AbstractThe front cover artwork is provided by Naoki Aratani, Hiroko Yamada, and co‐workers at NAIST (Japan). This picture shows Mt. Fuji, the highest mountain in Japan, reflecting off the water surface in a beautiful symmetrical form. On this image, we superimposed the mirror‐image CPL spectra obtained in this study for a 1,8‐naphthylene‐linked face‐to‐face coronene dimer. Read the full text of the Communication at 10.1002/cptc.202100166.

Front Cover: Mirror‐Image Cofacial Coronene Dimers Characterized by CD and CPL Spectroscopy: A Twisted Bilayer Nanographene (ChemPhotoChem 11/2021)

The Front Cover shows Mt. Fuji, the highest mountain in Japan, reflecting off the water surface in a beautiful symmetrical form. On this image, we superimposed the mirror-image circularly polarized luminescence (CPL) spectra of the 1,8-naphthylene-linked face-to-face coronene dimer reported in this study. More information can be found in the Communication by Hiroko Yamada, Naoki Aratani and co-workers.

Mirror‐Image Cofacial Coronene Dimers Characterized by CD and CPL Spectroscopy: A Twisted Bilayer Nanographene

AbstractBilayer graphene (BLG), which changes its electronic properties significantly depending on the interplanar distance and rotation angle, is a hot research topic. We successfully synthesized discrete BLG models by bridging two coronenes with a 1,8‐naphthalene linkage (C2). Owing to the large rotational energy, the syn‐ and anti‐atropisomers of C2 could be isolated. X‐ray crystallographic analysis revealed that the syn‐form (C2S) is an AB stacking type and the anti‐form (C2A) is a chiral twisted bilayer graphene (TBG) type. Therefore, the anti‐form of C2 was optically resolved and characterized by circular dichroism (CD) and circularly polarized luminescence (CPL) spectroscopy. The Cotton effect was observed in the CD spectra at the longest absorption wavelength, and the corresponding bisignate CPL spectra were observed. The dissymmetry factors |glum| of C2A were found to be 2.0×10−3 at 435 nm and 1.2×10−3 at 450 nm.

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA-ellipticine and buckycatcher-fullerene.

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA-ellipticine intercalation complex, and the buckycatcher-C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller-Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA-ellipticine, and buckycatcher-C60, respectively. Such sizable discrepancy above "chemical accuracy" for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree-Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Continuum Modelling for Interacting Coronene Molecules with a Carbon Nanotube.

The production of single dimensional carbon structures has recently been made easier using carbon nanotubes. We consider here encapsulated coronene molecules, which are flat and circular-shaped polycyclic aromatic hydrocarbons, inside carbon nanotubes. Depending on the radius of the nanotube, certain specific configurations of the coronene molecules can be achieved that give rise to the formation of stacked columns or aid in forming nanoribbons. Due to their symmetrical structure, a coronene molecule may be modelled by three inner circular rings of carbon atoms and one outer circular ring of hydrogen atoms, while the carbon nanotube is modelled as a circular tube. Using the continuous model and the Lennard-Jones potential, we are able to analytically formulate an expression for the potential energy for a coronene dimer and coronene inside a carbon nanotube. Subsequently, stacking of coronene molecules inside a nanotube is investigated. We find that the minimum energy tilt angle of coronenes in a stack differs from that of a single coronene within the same nanotube. More specifically, for both (18, 0) and (19, 0) zigzag carbon nanotube, we find that the minimum energy tilt angles of the single coronene case (≈42 and ≈20 respectively) do not occur in the stack model.

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

Mass spectra of helium nanodroplets doped with H2 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 H2 molecules indicating a commensurate solvation shell that preserves the D6h symmetry of the substrate. No such feature is seen in the abundance of the coronene dimer through tetramer complexed with H2; 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.

Geometric and electronic structures of one-dimensionally polymerized coronene molecules

On the basis of density functional theory, we studied the energetics, geometries, and electronic structures of dicoronylene molecules (coronene dimers) and one-dimensionally polymerized coronene molecules. Our calculations showed that the relative stability of the coronene dimers and polymers strongly depends on the mutual arrangement of molecules with respect to the polymer or dimer directions: the coronene dimer and the polymer with six-membered rings in the intermolecular bonds are the most stable structures among the three possible intermolecular arrangements. The dimer possesses a large energy gap between the highest occupied and lowest unoccupied states, which is narrower than that of the isolated coronene. The polymers are semiconductors with a moderate direct band gap and a large indirect gap for the polymers with six- and four-membered rings and that with five-membered rings, respectively, indicating the possibility of the electronic structure control of polymers via the intermolecular arrangement.

Beyond Energies: Geometries of Nonbonded Molecular Complexes as Metrics for Assessing Electronic Structure Approaches.

Electronic structure approaches for calculating intermolecular interactions have traditionally been benchmarked almost exclusively on the basis of energy-centric metrics. Herein, we explore the idea of utilizing a metric related to geometry. On a diverse series of noncovalently interacting systems of different sizes, from the water dimer to the coronene dimer, we evaluate a variety of electronic structure approximations with respect to their abilities to reproduce coupled-cluster-level geometries. Specifically, we examine Hartree-Fock, second-order Møller-Plesset perturbation theory (MP2), attenuated MP2, scaled MP2, and a number of density functionals, many of which include empirical or nonempirical van der Waals dispersion corrections. We find a number of trends that transcend system size and interaction type. For instance, functionals incorporating VV10 nonlocal correlation tend to yield highly accurate geometries; ωB97X-V and B97M-V, in particular, stand out. We establish that intermolecular distance, as measured by, e.g., the center-of-mass separation of two molecules, is the geometric parameter that deviates most profoundly among the various methods. This property of the equilibrium intermolecular separation, coupled with its accessibility via a small series of well-defined single-point calculations, makes it an ideal metric for the development and evaluation of electronic structure methods.

Mapping the genome of meta-generalized gradient approximation density functionals: the search for B97M-V.

A meta-generalized gradient approximation density functional paired with the VV10 nonlocal correlation functional is presented. The functional form is selected from more than 10(10) choices carved out of a functional space of almost 10(40) possibilities. Raw data come from training a vast number of candidate functional forms on a comprehensive training set of 1095 data points and testing the resulting fits on a comprehensive primary test set of 1153 data points. Functional forms are ranked based on their ability to reproduce the data in both the training and primary test sets with minimum empiricism, and filtered based on a set of physical constraints and an often-overlooked condition of satisfactory numerical precision with medium-sized integration grids. The resulting optimal functional form has 4 linear exchange parameters, 4 linear same-spin correlation parameters, and 4 linear opposite-spin correlation parameters, for a total of 12 fitted parameters. The final density functional, B97M-V, is further assessed on a secondary test set of 212 data points, applied to several large systems including the coronene dimer and water clusters, tested for the accurate prediction of intramolecular and intermolecular geometries, verified to have a readily attainable basis set limit, and checked for grid sensitivity. Compared to existing density functionals, B97M-V is remarkably accurate for non-bonded interactions and very satisfactory for thermochemical quantities such as atomization energies, but inherits the demonstrable limitations of existing local density functionals for barrier heights.

Third-order nonlinear optical properties of one-dimensional open-shell molecular aggregates composed of phenalenyl radicals.

The impact of intermolecular interactions on the third-order nonlinear optical (NLO) properties of open-shell molecular aggregates has been elucidated by considering one-dimensional aggregates of π-π stacked phenalenyl radicals with different intermolecular distances and the long-range corrected spin-unrestricted density functional theory method. In the phenalenyl dimer, which can be considered as a diradicaloid system, the diradical character strongly depends on the intermolecular distance, and the larger the intermolecular distance is, the larger the diradical character becomes. Then, around the equilibrium stacking distance that corresponds to an intermediate diradical character, its second hyperpolarizability (γ) is maximized and its value per monomer exhibits about a 30-fold enhancement with respect to the isolated phenalenyl monomer. This suggests that equilibrium is an optimal compromise between localization and delocalization of the radical electron pairs in such pancake bonding. No such effect was observed for the closed-shell coronene dimer. Moreover, when going from the dimer (diradical) to the tetramer (tetraradical), the γ-enhancement ratio increases nonlinearly with the aggregate size, whereas switching from the singlet to the highest spin (quintet) state causes a significant reduction of γ. Finally, for the tetramer, another one-order enhancement of γ is achieved for the dicationic singlet relative to its singlet neutral state. These results demonstrate the key role of intermolecular π-π stacking interactions and charge in open-shell (supra)molecular systems to achieve enhanced third-order NLO properties.

Coronene encapsulation in single-walled carbon nanotubes: stacked columns, peapods, and nanoribbons.

Encapsulation of coronene inside single-walled carbon nanotubes (SWNTs) was studied under various conditions. Under high vacuum, two main types of molecular encapsulation were observed by using transmission electron microscopy: coronene dimers and molecular stacking columns perpendicular or tilted (45-60°) with regard to the axis of the SWNTs. A relatively small number of short nanoribbons or polymerized coronene molecular chains were observed. However, experiments performed under an argon atmosphere (0.17 MPa) revealed reactions between the coronene molecules and the formation of hydrogen-terminated graphene nanoribbons. It was also observed that the morphology of the encapsulated products depend on the diameter of the SWNTs. The experimental results are explained by using density functional theory calculations through the energies of the coronene molecules inside the SWNTs, which depend on the orientation of the molecules and the diameter of the tubes.

Density functional theory study of the interaction of vinyl radical, ethyne, and ethene with benzene, aimed to define an affordable computational level to investigate stability trends in large van der Waals complexes

Our purpose is to identify a computational level sufficiently dependable and affordable to assess trends in the interaction of a variety of radical or closed shell unsaturated hydro-carbons A adsorbed on soot platelet models B. These systems, of environmental interest, would unavoidably have rather large sizes, thus prompting to explore in this paper the performances of relatively low-level computational methods and compare them with higher-level reference results. To this end, the interaction of three complexes between non-polar species, vinyl radical, ethyne, or ethene (A) with benzene (B) is studied, since these species, involved themselves in growth processes of polycyclic aromatic hydrocarbons (PAHs) and soot particles, are small enough to allow high-level reference calculations of the interaction energy ΔEAB. Counterpoise-corrected interaction energies ΔEAB are used at all stages. (1) Density Functional Theory (DFT) unconstrained optimizations of the A-B complexes are carried out, using the B3LYP-D, ωB97X-D, and M06-2X functionals, with six basis sets: 6-31G(d), 6-311 (2d,p), and 6-311++G(3df,3pd); aug-cc-pVDZ and aug-cc-pVTZ; N07T. (2) Then, unconstrained optimizations by Møller-Plesset second order Perturbation Theory (MP2), with each basis set, allow subsequent single point Coupled Cluster Singles Doubles and perturbative estimate of the Triples energy computations with the same basis sets [CCSD(T)//MP2]. (3) Based on an additivity assumption of (i) the estimated MP2 energy at the complete basis set limit [EMP2/CBS] and (ii) the higher-order correlation energy effects in passing from MP2 to CCSD(T) at the aug-cc-pVTZ basis set, ΔECC-MP, a CCSD(T)/CBS estimate is obtained and taken as a computational energy reference. At DFT, variations in ΔEAB with basis set are not large for the title molecules, and the three functionals perform rather satisfactorily even with rather small basis sets [6-31G(d) and N07T], exhibiting deviation from the computational reference of less than 1 kcal mol(-1). The zero-point vibrational energy corrected estimates Δ(EAB+ZPE), obtained with the three functionals and the 6-31G(d) and N07T basis sets, are compared with experimental D0 measures, when available. In particular, this comparison is finally extended to the naphthalene and coronene dimers and to three π-π associations of different PAHs (R, made by 10, 16, or 24 C atoms) and P (80 C atoms).

ωB97X-V: A 10-parameter, range-separated hybrid, generalized gradient approximation density functional with nonlocal correlation, designed by a survival-of-the-fittest strategy

A 10-parameter, range-separated hybrid (RSH), generalized gradient approximation (GGA) density functional with nonlocal correlation (VV10) is presented. Instead of truncating the B97-type power series inhomogeneity correction factors (ICF) for the exchange, same-spin correlation, and opposite-spin correlation functionals uniformly, all 16,383 combinations of the linear parameters up to fourth order (m = 4) are considered. These functionals are individually fit to a training set and the resulting parameters are validated on a primary test set in order to identify the 3 optimal ICF expansions. Through this procedure, it is discovered that the functional that performs best on the training and primary test sets has 7 linear parameters, with 3 additional nonlinear parameters from range-separation and nonlocal correlation. The resulting density functional, ωB97X-V, is further assessed on a secondary test set, the parallel-displaced coronene dimer, as well as several geometry datasets. Furthermore, the basis set dependence and integration grid sensitivity of ωB97X-V are analyzed and documented in order to facilitate the use of the functional.

The accuracy of quantum chemical methods for large noncovalent complexes.

We evaluate the performance of the most widely used wavefunction, density functional theory, and semiempirical methods for the description of noncovalent interactions in a set of larger, mostly dispersion-stabilized noncovalent complexes (the L7 data set). The methods tested include MP2, MP3, SCS-MP2, SCS(MI)-MP2, MP2.5, MP2.X, MP2C, DFT-D, DFT-D3 (B3-LYP-D3, B-LYP-D3, TPSS-D3, PW6B95-D3, M06-2X-D3) and M06-2X, and semiempirical methods augmented with dispersion and hydrogen bonding corrections: SCC-DFTB-D, PM6-D, PM6-DH2 and PM6-D3H4. The test complexes are the octadecane dimer, the guanine trimer, the circumcoronene…adenine dimer, the coronene dimer, the guanine-cytosine dimer, the circumcoronene…guanine-cytosine dimer, and an amyloid fragment trimer containing phenylalanine residues. The best performing method is MP2.5 with relative root mean square deviation (rRMSD) of 4 %. It can thus be recommended as an alternative to the CCSD(T)/CBS (alternatively QCISD(T)/CBS) benchmark for molecular systems which exceed current computational capacity. The second best non-DFT method is MP2C with rRMSD of 8 %. A method with the most favorable "accuracy/cost" ratio belongs to the DFT family: BLYP-D3, with an rRMSD of 8 %. Semiempirical methods deliver less accurate results (the rRMSD exceeds 25 %). Nevertheless, their absolute errors are close to some much more expensive methods such as M06-2X, MP2 or SCS(MI)-MP2, and thus their price/performance ratio is excellent.

Attenuated second-order Møller–Plesset perturbation theory: performance within the aug-cc-pVTZ basis

Attenuated second-order Møller-Plesset perturbation theory (MP2) within the finite aug-cc-pVTZ (aTZ) basis set is developed for inter- and intra-molecular non-bonded interactions. A single attenuation parameter is optimized on the S66 database of 66 intermolecular interactions, leading to a very large RMS error reduction by a factor of greater than 5 relative to standard MP2/aTZ. Attenuation introduces an error of opposite sign to basis set superposition error (BSSE) and overestimation of dispersion interactions in finite basis MP2. A variety of tests including the S22 set, conformer energies of peptides, alkanes, sugars, sulfate-water clusters, and the coronene dimer establish the transferability of the MP2(terfc, aTZ) model to other inter and intra-molecular interactions. Direct comparisons against attenuation in the smaller aug-cc-pVDZ basis shows that MP2(terfc, aTZ) often significantly outperforms MP2(terfc, aDZ), although at higher computational cost. MP2(terfc, aDZ) and MP2(terfc, aTZ) often out-perform MP2 at the complete basis set limit. Comparison of the two attenuated MP2 models against each other and against attenuation using non-augmented basis sets gives insight into the error cancellation responsible for their remarkable success.

Buckyplates and Buckybowls: Examining the Effects of Curvature on π–π Interactions

π-π interactions are integral to many areas of chemistry, biochemistry, and materials science. Here we use electronic structure theory to analyze how π-π interactions change as the π-systems are curved in model complexes based on coronene and corannulene dimers. Curvature redistributes electronic charge in the π-cloud and creates a dipole moment in these systems, leading to enhanced intermolecular electrostatic interactions in the concave-convex (nested) geometries that are the focus of this work. Curvature of both monomers also has a geometric effect on the interaction by decreasing the average C-C distance between monomers and by increasing the magnitude of both favorable London dispersion interactions and unfavorable exchange-repulsion interactions. Overall, increasing curvature in nested π-π interactions leads to more favorable interaction energies regardless of the native state of the monomers, except at short distances where the most highly curved systems are less favorable as exchange repulsion terms begin to dominate the interaction.