Size-Transferable Prediction of Excited State Properties for Molecular Assemblies with a Machine Learning Exciton Model

Computational modeling of the excited states of molecular aggregates faces significant computational challenges and size heterogeneity. Current machine learning (ML) models, typically trained on specific-sized aggregates, struggle with scalability. We found that the exciton model Hamiltonian of large aggregates can be decomposed into dimer pairs, allowing an ML model trained on dimers to reconstruct Hamiltonians for aggregates of any size. We also proposed a new method to address the phase-correction problem by introducing coupling terms' approximations. Our model accurately predicted the excitation energies of the trimer and tetramer of perylene and tetracene and estimated S1 oscillator strengths of perylene aggregates. Leveraging our ML model, the optical gaps of nanosized perylene aggregates with up to 50 monomers are analyzed, qualitatively revealing the role of different couplings on their size dependency. Future work will explore transferability across different monomers to predict optical properties in heterogeneous assemblies.

Highly Selective On-Surface Dehydrogenative Aromatization of n-Hexyl to Phenyl Substituents

Dehydrogenative aromatization of alkyl substituents represents a powerful approach to aryl-substituted functional molecules. However, the inertness of alkyl groups and the need for harsh reaction conditions accompanied by low product selectivity hamper its widespread applications. Here, we demonstrate the highly selective on-surface thermal aromatization of n-hexyl substituents on polycyclic aromatic hydrocarbons (PAHs) to their phenyl-substituted analogues under mild conditions. After depositing two representative precursor molecules, n-hexyl-substituted hexaphenylbenzene (HPB-Hex) and bianthryl octacarboxylic tetraimide (BATI-Hex), onto a pre-heated Au(111) substrate, dehydroaromatization of the peripheral n-hexyl groups into phenyl rings occurs, following the planarization of hexaphenylbenzene and bianthryl core. This process involves sequential intramolecular C-C bond rotations, dehydrogenation, and cyclodehydrogenation reactions, yielding phenyl-substituted hexabenzocoronene (HBC-Ph) and bisanthene octacarboxylic tetraimides (BSTI-Ph). The reaction sequences were monitored using scanning tunneling microscopy (STM) and bond-resolved non-contact atomic force microscopy (nc-AFM), offering structural proof of both intermediates and final products. These experimental techniques were complemented by density functional theory (DFT) simulations, which facilitated the detection of crucial steps in the conversion of n-hexyl to phenyl groups. Moreover, the effect of alkyl aromatization on the electronic properties of the newly formed aromatic hydrocarbons was elucidated using scanning tunneling spectroscopy (STS).

Roto-Cyclization of 4-Bromopicene in On-Surface Synthesis

Progress toward single-molecule electronics relies on a thorough understanding of local physico-chemical processes and development of synthetic routines for controlled hetero-coupling. We demonstrate a structurally unexpected ring closure process for a homo-coupled 4,4'-bipicenyl, realized in on-surface synthesis. An initial covalent C-C coupling of 4-bromopicene locks at lower temperatures the position and geometrically shields part of 4,4'-bipicenyl. Employing this effect of shielding might offer a path toward controlled stepwise hetero-coupling. At higher temperatures, a thermally activated three-dimensional rotation upon hydrogen dissociation, a dehydrogenative roto-cyclization, lifts the surface-dimensionality restriction, and leads to the formation of a perylene. Thereby, the shielded molecular part becomes accessible again.

Size Isn't Everything: Geometric Tuning in Polycyclic Aromatic Hydrocarbons and Its Implications for Carbon Nanodots

Recent developments in light-emitting carbon nanodots and molecular organic semiconductors have seen renewed interest in the properties of polycyclic aromatic hydrocarbons (PAHs) as a family. The networks of delocalized π electrons in sp2-hybridized carbon grant PAHs light-emissive properties right across the visible spectrum. However, the mechanistic understanding of their emission energy has been limited due to the ground state-focused methods of determination. This computational chemistry work, therefore, seeks to validate existing rules and elucidate new features and characteristics of PAHs that influence their emissions. Predictions based on (time-dependent) density functional theory account for the full 3-dimensional electronic structure of ground and excited states and reveal that twisting and near-degeneracies strongly influence emission spectra and may therefore be used to tune the color of PAHs and, hence, carbon nanodots. We particularly note that the influence of twisting goes beyond torsional destabilization of the ground-state and geometric relaxation of the excited state, with a third contribution associated with the electric transition dipole. Symmetries and peri-condensation may also have an effect, but this could not be statistically confirmed. In pursuing this goal, we demonstrate that with minimal changes to molecular size, the entire visible spectrum may be spanned by geometric modification alone; we have also provided a first estimate of emission energy for 35 molecules currently lacking published emission spectra as well as clear guidelines for when more sophisticated computational techniques are required to predict the properties of PAHs accurately.

Substitution Effect of a Single Nitrogen Atom on π-Electronic Systems of Linear Polycyclic Aromatic Hydrocarbons (PAHs): Theoretically Visualized Coexistence of Mono- and Polycyclic π-Electron Delocalization

We theoretically investigated the nitrogen substitution effect on the molecular structure and π-electron delocalization in linear nitrogen-substituted polycyclic aromatic hydrocarbons (N-PAHs). Based on the optimized molecular structures and magnetic field-induced parameters of fused bi- and tricyclic linear N-PAHs, we found that the local π-electron delocalization of subcycles (e.g., mono- and bicyclic constituent moieties) in linear N-PAHs is preserved, despite deviation from ideal structures of parent monocycles. The introduction of a fused five-membered ring with a pyrrolic N atom (N-5MR) in linear N-PAHs significantly perturbs the π-electronic condition of the neighboring fused six-membered ring (6MR). Monocyclic pyrrole exhibits substantial bond length alternations, strongly influencing the π-electronic systems of both the fused N-5MR and 6MR in linear N-PAHs, depending on the location of shared covalent bonds. A fused six-membered ring with a graphitic N atom in an indolizine moiety cannot generate monocyclic π-electron delocalization but instead contributes to the formation of polycyclic π-electron delocalization. This is evidenced by bifurcated diatropic ring currents induced by an external magnetic field. In conclusion, the satisfaction of Hückel's 4n + 2 rule for both mono- and polycycles is crucial for understanding the overall π-electron delocalization. It is crucial to consider the unique characteristics of the three types of substituted N atoms and the spatial arrangement of 5MR and 6MR in N-PAHs.

Wood-based superblack

Light is a powerful and sustainable resource, but it can be detrimental to the performance and longevity of optical devices. Materials with near-zero light reflectance, i.e. superblack materials, are sought to improve the performance of several light-centered technologies. Here we report a simple top-down strategy, guided by computational methods, to develop robust superblack materials following metal-free wood delignification and carbonization (1500 °C). Subwavelength severed cells evolve under shrinkage stresses, yielding vertically aligned carbon microfiber arrays with a thickness of ~100 µm and light reflectance as low as 0.36% and independent of the incidence angle. The formation of such structures is rationalized based on delignification method, lignin content, carbonization temperature and wood density. Moreover, our measurements indicate a laser beam reflectivity lower than commercial light stoppers in current use. Overall, the wood-based superblack material is introduced as a mechanically robust surrogate for microfabricated carbon nanotube arrays.

Challenges in the Direct Detection of Chirality-induced Spin Selectivity: Investigation of Foldamer-based Donor-acceptor Dyads

Over the past two decades, the chirality-induced spin selectivity (CISS) effect was reported in several experiments disclosing a unique connection between chirality and electron spin. Recent theoretical works highlighted time-resolved Electron Paramagnetic Resonance (trEPR) as a powerful tool to directly detect the spin polarization resulting from CISS. Here, we report a first attempt to detect CISS at the molecular level by linking the pyrene electron donor to the fullerene acceptor with chiral peptide bridges of different length and electric dipole moment. The dyads are investigated by an array of techniques, including cyclic voltammetry, steady-state and transient optical spectroscopies, and trEPR. Despite the promising energy alignment of the electronic levels, our multi-technique analysis reveals no evidence of electron transfer (ET), highlighting the challenges of spectroscopic detection of CISS. However, the analysis allows the formulation of guidelines for the design of chiral organic model systems suitable to directly probe CISS-polarized ET.

Circumpentacene with Open-Shell Singlet Diradical Character

Circumacenes (CAs) are a distinctive type of benzenoid polycyclic aromatic hydrocarbons where an acene unit is completely enclosed by a layer of outer fused benzene rings. Despite their unique structures, the synthesis of CAs is challenging, and until recently, the largest CA molecule synthesized was circumanthracene. In this study, we report the successful synthesis of an extended circumpentacene derivative 1, which represents the largest CA molecule synthesized to date. Its structure was confirmed by X-ray crystallographic analysis and its electronic properties were systematically investigated by both experiments and theoretical calculations. It shows a unique open-shell diradical character due to the existence of extended zigzag edges, with a moderate diradical character index (y0 =39.7 %) and a small singlet-triplet energy gap (ΔES-T =-4.47 kcal/mol). It exhibits a dominant local aromatic character with π-electrons delocalized in the individual aromatic sextet rings. It has a small HOMO-LUMO energy gap and displays amphoteric redox behavior. The electronic structures of its dication and dianion can be considered as doubly charged structures in which two coronene units are fused with a central aromatic benzene ring. This study provides a new route toward stable multizigzag-edged graphene-like molecules with open-shell di/polyradical character.

Optimized attenuated interaction: Enabling stochastic Bethe-Salpeter spectra for large systems

We develop an improved stochastic formalism for the Bethe-Salpeter equation (BSE), based on an exact separation of the effective-interaction W into two parts, W = (W - vW) + vW, where the latter is formally any translationally invariant interaction, vW(r - r'). When optimizing the fit of the exchange kernel vW to W, using a stochastic sampling W, the difference W - vW becomes quite small. Then, in the main BSE routine, this small difference is stochastically sampled. The number of stochastic samples needed for an accurate spectrum is then largely independent of system size. While the method is formally cubic in scaling, the scaling prefactor is small due to the constant number of stochastic orbitals needed for sampling W.

Aggregation-Induced Emission of Contorted Polycondensed Aromatic Hydrocarbons Made by Edge Extension Using a Palladium-Catalyzed Cyclopentannulation Reaction

Contorted polycyclic aromatic hydrocarbons (PAHs), CPA1-2 and CPB1-2, bearing peripheral five-membered rings were synthesized employing a palladium-catalyzed cyclopentannulation reaction using specially designed diaryl acetylene synthons TPE and TPEN with commercially available dibromo- anthracene DBA and bianthracene DBBA derivatives. The resulting target compounds CPA1-2 and CPB1-2 were isolated in excellent yield and found to be highly soluble in common organic solvents, which allowed for their structural characterization and investigation of the photophysical properties, disclosing their aggregation-induced emission (AIE) properties in THF at selective concentration ranges of water fractions in the solvent mixture. Examination of the contorted PAH structures by means of density functional theory (DFT) revealed higher electronic conjugation in the more rigid and planar anthracene-containing CPA1-2 derivatives when compared to the twisted bianthracene-bearing moieties CBPA1-2 with HOMO-LUMO bandgaps (ΔE) of ∼2.32 eV for the former PAHs and ∼2.78 eV for the latter ones.

Optical Characterization of Bitumen: Inspecting the Quantum Size Effect in the Nanostructured Phase

Bitumen (BIT) is an oil sub-product with many applications. A variety of literature reports the analysis of its optical properties, as the fluorescence and the Fourier transform infrared (FT-IR) optical transmittance (OT). We have performed photoluminescence (PL) on the visible region of the spectrum and OT on the visible-NIR in solutions of BIT in toluene at different concentrations. The results revealed a nonlinear, concentration-dependent effect attributed to intermolecular solvent-solute interactions. PL analysis of pure BIT as a function of the temperature and the laser pumping power pointed to the existence of nanosized crystalline inclusions on the amorphous BIT matrix with noticeable optical properties. The inference was confronted with the results of X-ray diffraction studies and literature reports. The possibility of the occurrence of a quantum size effect governing the luminescence is considered.

Optical Properties of Polycyclic Aromatic Hydrocarbon Clusters with Oxygenated Substitutions: A Theoretical Analysis

The introduction of functional groups at high coverage levels can have significant impacts on the band structures of polycyclic aromatic hydrocarbon (PAH) clusters. The HOMO-LUMO gaps are highly sensitive to the type and distribution of functional groups. An in-house method is proposed to build PAH (naphthalene, pyrene, coronene, and ovalene) clusters with surface functionalization of -OH, -COOH and -CHO groups using the DFT method. The -CHO groups are found to reduce the gap value the most, but exceptions exist due to the spatial distribution of functional groups. Considering the impact of -CHO groups only, we can approximate that the impact of functional groups lies in the range of 0.14-0.89 eV. Applying further analysis on the possible energy number of energy transitions of substituted PAH clusters, it is shown that PAH clusters with oxygenated functions still behave like an indirect band gap material. The coupling effect of PAH stacking and PAH size is also addressed. A simple expression is proposed to estimate the bandgap of a mixed system using the HOMO and LUMO energy of the two components. Further attempts are made to interpret recent experiments from the impact of PAH stacking, PAH size, and functional groups.

Low-Dimensional Porous Carbon Networks Using Single-/Triple-Coupling Polycyclic Hydrocarbon Precursors

Polycyclic hydrocarbons (PHs) share the same hexagonal structure of sp² carbons as graphene but possess an energy gap due to quantum confinement effect. PHs can be synthesized by a bottom-up strategy starting from small building blocks covalently bonded into large 2D organic sheets. Further investigation of the role of the covalent bonding/coupling ways on their electronic properties is needed. Here, we demonstrate a surface-mediated synthesis of hexa-peri-hexabenzocoronene (HBC) and its extended HBC oligomers (dimers, trimers, and tetramers) via single- and triple-coupling ways and reveal the implication of different covalent bonding on their electronic properties. High-resolution low-temperature scanning tunneling microscopy and noncontact atomic force microscopy are employed to in situ determine the atomic structures of as-synthesized HBC oligomers. Scanning tunneling spectroscopy measurements show that the length extension of HBC oligomers narrows the energy gap between highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). Furthermore, the energy gaps of triple-coupling HBC oligomers are smaller and decrease more significantly than that of the single-coupling ones. We hypothesize that the triple coupling promotes a more effective delocalization of π-electrons than the single coupling, according to density functional theory calculations. We also demonstrate that the HBC oligomers can further extend across the substrate steps to achieve conjugated polymers and large-area porous carbon networks.

DFT calculation in design of near-infrared absorbing nitrogen-doped graphene quantum dots

The near-infrared light (NIR) absorption of nitrogen-doped graphene quantum dots (NGQDs) containing different N-doping sites is systematically investigated with density functional theory (DFT) and time-dependent density functional theory (TD-DFT) calculations with Perdew-Burke-Ernzerhof (PBE) functionals. The results show that the ultra-small HOMO-LUMO gaps (0.3-1.0 eV) of various N-doping structures (graphitic, amino, and pyridinic at center, and graphitic at edge) are attributed to the spin-polarization of the energy states, which effectively enhances the NIR absorption for NGQDs. Overall, the graphitic N-doping structure exhibits the best NIR absorption. Moreover, the electron attraction effect of the different N-sites is found to be crucial for the LUMO level, where stronger electron attraction lowers the LUMO energy. This work provides critical insight in further design of NGQDs for NIR absorption.

Modifications of EHPDB Physical Properties through Doping with Fe2O3 Nanoparticles (Part II)

The aim of our study was to analyze the influence of various concentrations of γ-Fe₂O₃ nanoparticles on the physical properties of the liquid crystalline ferroelectric SmC* phase, as well as to check the effect of introducing nanoparticles in the LC matrix on their properties in the prepared five nanocomposites. UV-vis spectroscopy showed that the admixture reduced the absorption of nanocomposites in the UV range, additional absorption bands appeared, and all nanocomposites were transparent in the range of 500–850 nm. The molecular dynamics in particular phases of the nanocomposites were investigated by the dielectric spectroscopy method, and it was found that nanoparticles caused a significant increase in the dielectric constant at low frequencies, a strong modification of the dielectric processes in the SmC* phase, and the emergence of new relaxation processes for the highest dopant concentrations. SQUID magnetometry allowed us to determine the magnetic nature of the nanoparticles used, and to show that the blocked state of nanoparticles was preserved in nanocomposites (hysteresis loops were also registered in the ferroelectric SmC* phase). The dependence of the coercive field on the admixture concentration and the widening of the hysteresis loop in nanocomposites in relation to pure nanoparticles were also found. In turn, the FT-MIR spectroscopy method was used to check the influence of the impurity concentration on the formation/disappearance or modification of the absorption bands, and the modification of both the FWHM and the maximum positions for the four selected vibrations in the MIR range, as well as the discontinuous behavior of these parameters at the phase transitions, were found.

Solid-phase excitation-emission matrix spectroscopy for chemical analysis of combustion aerosols

Exposure to ultrafine combustion aerosols such as particulate matter (PM) from residential woodburning, forest fires, cigarette smoke, and traffic emission have been linked to adverse health outcomes. Excitation-emission matrix (EEM) spectroscopy presents a sensitive and cost-effective alternative for analysis of PM organic fraction. However, as with other analytical chemistry methods, the miniaturization is hindered by a solvent extraction step and a need for benchtop instrumentation. We present a methodology for collecting and in-situ analysis of airborne nanoparticles that eliminates labor-intensive sample preparation and miniaturizes the detection platform. Nanoparticles are electrostatically collected onto a transparent substrate coated with solid-phase (SP) solvent-polydimethylsiloxane (PDMS). The PM organic fraction is extracted into PDMS and analyzed in-situ, thus avoiding liquid-phase extraction. In the SP-EEM analysis, we evaluated external and internal excitation schemes. Internal excitation shows the lowest scattering interference but leads to signal masking from PDMS fluorescence for λ<250nm. The external excitation EEM spectra are dependent on the excitation light incident angle; ranges of 30-40° and 55-65° show the best results. SP-EEM spectra of woodsmoke and cigarette smoke samples are in good agreement with the EEM spectra of liquid-phase extracts. The SP-EEM technique can be used to develop wearable sensors for exposure assessments and environmental monitoring.

On the Aromatic Stabilization of Fused Polycyclic Aromatic Hydrocarbons

The thermodynamic properties and band gap energies were evaluated for six ortho- and peri-fused polycyclic aromatic hydrocarbons (PAHs): triphenylene; benzo[a]pyrene; benzo[e]pyrene; perylene; benzo[ghi]perylene; coronene. The standard molar enthalpies of formation in the crystalline state and the standard molar enthalpies of sublimation were measured by high precision combustion calorimetry and Knudsen effusion methodology, respectively. The combination of the molar enthalpies of formation in the crystalline state with the respective enthalpies of sublimation was used to evaluate the energetics of the progressive peri-fusion of the aromatic moieties from triphenylene to coronene aiming to investigate the hypothetical superaromaticity character of coronene. The linear trend of the enthalpy of formation in crystalline and gaseous phases in the series (from benzo[e]pyrene to coronene) is an irrefutable indication of a non-superaromaticity character of coronene. High accurate thermodynamic properties of sublimation (volatility, enthalpy, and entropy of sublimation) were derived by the measurement of vapor pressures as a function of temperature, using a Knudsen/quartz crystal effusion methodology. Furthermore, the π-electronic conjugation of these compounds was explored by evaluation of the optical band gaps along with this series of compounds. The morphology of perylene, benzo[ghi]perylene, and coronene thin films, deposited by physical vapor deposition onto transparent conductive oxide substrates (ITO and FTO), was used to analyze the nucleation and growth mechanisms. The morphologies observed were found to be related to the cohesive energy and entropy of the bulk.

Theoretical investigation of electronic and optical properties of nitrogen doped triangular shaped graphene quantum dots

The electronic and optical properties of graphene quantum dots can be significantly tailored by doping it with heteroatoms, thus extending its potential applications. In this work, we have employed time-dependent density functional theory to systematically explore the effect of introduction of nitrogen atoms in varying concentration at pyridinic and graphitic configuration in armchair and zigzag-edged triangular shaped graphene quantum dots (TQDs) of different sizes. Our results indicate that the electronic band-gap in these N-doped systems can be effectively tuned by varying the configuration as well as concentration of dopants and nature of edge-termination. The variation of electronic band-gap is critically determined by the localized/delocalized nature of molecular orbitals and presence of additional energy levels due to dopant nitrogen atoms. However, the significance of these extra energy levels in modulating the optical properties (appearance of characteristic N-dopant absorption peaks) becomes conspicuous only for specific configuration and concentration of nitrogen atoms. In addition, our studies have attributed the strong dependence of blue/red-shift of absorption spectra and variation in the peak profile to position as well as concentration of dopant atoms and edge-termination pattern. Further, it is observed that the effect of increasing size of TQDs on the strength of most intense absorption peak of pyridinic N-doped TQDs is remarkably different from graphitic N-doped systems. This selective manipulation of optical properties in TQDs due to different N-doping pattern can open up new frontiers for rational design of novel optoelectronic devices.

Excitation-Emission Matrix Spectroscopy for Analysis of Chemical Composition of Combustion Generated Particulate Matter

Analysis of particulate matter (PM) is important for the assessment of human exposures to potentially harmful agents, notably combustion-generated PM. Specifically, polycyclic aromatic hydrocarbons (PAHs) found in ultrafine PM have been linked to cardiovascular diseases and carcinogenic and mutagenic effects. In this study, we quantify the presence and concentrations of PAHs with lower molecular weight (LMW, 126 < MW < 202) and higher molecular weight (HMW, 226 < MW < 302), i.e., smaller and larger than Pyrene, in combustion-generated PM using excitation-emission matrix (EEM) fluorescence spectroscopy. Laboratory combustion PM samples were generated in a laminar diffusion inverted gravity flame reactor (IGFR) operated on ethylene and ethane. Fuel dilution by Ar in 0% to 90% range controlled the flame temperature. The colder flames result in lower PM yields however, the PM PAH content increases significantly. Temperature thresholds for PM transition from low to high organic carbon content were characterized based on the maximum flame temperature (T_max,c ∼ 1791 to 1857 K) and the highest soot luminosity region temperature (T*_c ∼ 1600 to 1650K). Principal component regression (PCR) analysis of the EEM spectra of IGFR samples correlates to GCMS data with R² = 0.988 for LMW and 0.998 for HMW PAHs. PCR-EEM analysis trained on the IGFR samples was applied to PM samples from woodsmoke and diesel exhaust, the model accurately predicts HMW PAH concentrations with R² = 0.976 and overestimates LMW PAHs.

Synthesis of Circumpyrene by Alkyne Benzannulation of Brominated Dibenzo[hi,st]ovalene

A transition-metal catalyzed alkyne benzannulation allowed an unprecedented synthesis of circumpyrene, starting from 3,11-dibromo-6,14-dimesityldibenzo[hi,st]ovalene (DBOV). The circumpyrene was characterized by a combination of NMR, mass spectrometry, and single-crystal X-ray diffraction analysis, revealing its multizigzag-edged structure. Two newly introduced C═C bonds in circumpyrene strongly perturbed the electronic structures of DBOV, as evidenced by increased optical and electrochemical energy gaps. This is in good agreement with an increased number of Clar's sextets as well as a decreased number of π-electrons in the conjugation pathway of circumpyrene, according to anisotropy of the induced current density (ACID) calculations. The present approach opens a new avenue to multizigzag-edged nanographenes and offers insights into their (opto)electronic properties.