Molecularly induced order promotes charge separation through delocalized charge-transfer states at donor-acceptor heterojunctions.

The energetic landscape at the interface between electron donating and accepting molecular materials favors efficient conversion of intermolecular charge-transfer (CT) states into free charge carriers (FCC) in high-performance organic solar cells. Here, we elucidate how interfacial energetics, charge generation and radiative recombination are affected by molecular arrangement. We experimentally determine the CT dissociation properties of a series of model, small molecule donor-acceptor blends, where the used acceptors (B2PYMPM, B3PYMPM and B4PYMPM) differ only in the nitrogen position of their lateral pyridine rings. We find that the formation of an ordered, face-on molecular packing in B4PYMPM is beneficial to efficient, field-independent charge separation, leading to fill factors above 70% in photovoltaic devices. This is rationalized by a comprehensive computational protocol showing that, compared to the more amorphous and isotropically oriented B2PYMPM, the higher structural order of B4PYMPM molecules leads to more delocalized CT states. Furthermore, we find no correlation between the quantum efficiency of FCC radiative recombination and the bound or unbound nature of the CT states. This work highlights the importance of structural ordering at donor-acceptor interfaces for efficient FCC generation and shows that less bound CT states do not preclude efficient radiative recombination.

Characterizing Counterion-Dependent Aggregation of Rhodamine B by Classical Molecular Dynamics Simulations.

The aggregation in a solution of charged dyes such as Rhodamine B (RB) is significantly affected by the type of counterion, which can determine the self-assembled structure that in turn modulates the optical properties. RB aggregation can be boosted by hydrophobic and bulky fluorinated tetraphenylborate counterions, such as F5TPB, with the formation of nanoparticles whose fluorescence quantum yield (FQY) is affected by the degree of fluorination. Here, we developed a classical force field (FF) based on the standard generalized Amber parameters that allows modeling the self-assembling process of RB/F5TPB systems in water, consistent with experimental evidence. Namely, the classical MD simulations employing the re-parametrized FF reproduce the formation of nanoparticles in the RB/F5TPB system, while in the presence of iodide counterions, only RB dimeric species can be formed. Within the large, self-assembled RB/F5TPB aggregates, the occurrence of an H-type RB-RB dimer can be observed, a species that is expected to quench RB fluorescence, in agreement with the experimental data of FQY. The outcome provides atomistic details on the role of the bulky F5TPB counterion as a spacer, with the developed classical FF representing a step towards reliable modeling of dye aggregation in RB-based materials.

Predicting the Second-Order Nonlinear Optical Responses of Organic Materials: The Role of Dynamics.

The last 30 years have witnessed an ever-growing application of computational chemistry for rationalizing the nonlinear optical (NLO) responses of organic chromophores. More specifically, quantum chemical calculations proved highly helpful in gaining fundamental insights into the factors governing the magnitude and character of molecular first hyperpolarizabilities (ß), be they either intrinsic to the chromophore molecular structure and arising from symmetry, chemical substitution, or π-electron delocalization, or induced by external contributions such as the laser probe or solvation and polarization effects. Most theoretical reports assumed a rigid picture of the investigated systems, the NLO responses being computed solely at the most stable geometry of the chromophores. Yet, recent developments combining classical molecular dynamics (MD) simulations and DFT calculations have evidenced the significant role of structural fluctuations, which may induce broad distributions of NLO responses, and even generate them in some instances.This Account presents recent case studies in which theoretical simulations have highlighted these effects. The discussion specifically focuses on the simulation of the second-order NLO properties that can be measured experimentally either from Hyper-Rayleigh Scattering (HRS) or Electric-Field Induced Second Harmonic Generation (EFISHG). More general but technical topics concerning several aspects of the calculations of hyperpolarizabilities are instead discussed in the Supporting Information.Selected examples include organic chromophores, photochromic systems, and ionic complexes in the liquid phase, for which the effects of explicit solvation, concentration, and chromophore aggregation are emphasized, as well as large flexible systems such as peptide chains and pyrimidine-based helical polymers, in which the relative variations of the responses were shown to be several times larger than their average values. The impact of geometrical fluctuations is also illustrated for supramolecular architectures with the examples of nanoparticles formed by organic dipolar dyes in water solution, whose soft nature allows for large shape variations translating into huge fluctuations in time of their NLO response, and of self-assembled monolayers (SAMs) based on indolino-oxazolidine or azobenzene switches, in which the geometrical distortions of the photochromic molecules, as well as their orientational and positional disorder within the SAMs, highly impact their NLO response and contrast upon switching. Finally, the effects of the rigidity and fluidity of the surrounding are evidenced for NLO dyes inserted in phospholipid bilayers.

Temperature-Dependent Structural Phase Transition in Rubrene Single Crystals: The Missing Piece from the Charge Mobility Puzzle?

Accurate structural models for rubrene, the benchmark organic semiconductor, derived from synchrotron X-ray data in the temperature range of 100-300 K, show that its cofacially stacked tetracene backbone units remain blocked with respect to each other upon cooling to 200 K and start to slip below that temperature. The release of the blocked slippage occurs at approximately the same temperature as the hole mobility crossover. The blocking between 200 and 300 K is caused by a negative correlation between the relatively small thermal expansion along the crystallographic b-axis and the relatively large widening of the angle between herringbone-stacked tetracene units. DFT calculations reveal that this blocked slippage is accompanied by a discontinuity in the variation with temperature of the electronic couplings associated with hole transport between cofacially stacked tetracene backbones.

Direct Triple Annulations: A Way to Design Large Triazastarphenes with Intertwined Hexagonal Packing.

A new straightforward synthetic strategy has been elaborated to achieve star-shaped triazatrinaphthylene and, for the first time, triazatrianthrylene derivatives. Their solution- and solid-state properties were thoroughly characterized by cyclic voltammetry, UV-vis absorption spectroscopy, X-ray diffraction, and density functional theory calculations. Original hexagonal molecular arrangements were found in the crystal phase, which opens a new pathway for designing materials with improved three-dimensional charge-transport properties.

In Silico Ultrafast Nonlinear Spectroscopy Meets Experiments: The Case of Perylene Bisimide Dye.

Spectroscopy simulations are of paramount importance for the interpretation of experimental electronic spectra, the disentangling of overlapping spectral features, and the tracing of the microscopic origin of the observed signals. Linear and nonlinear simulations are based on the results drawn from electronic structure calculations that provide the necessary parameterization of the molecular systems probed by light. Here, we investigate the applicability of excited-state properties obtained from linear-response time-dependent density functional theory (TDDFT) in the description of nonlinear spectra by employing the pseudowavefunction approach and compare them with benchmarks from highly accurate RASSCF/RASPT2 calculations and with high temporal resolution experimental results. As a test case, we consider the prediction of femtosecond transient absorption and two-dimensional electronic spectroscopy of a perylene bisimide dye in solution. We find that experimental signals are well reproduced by both theoretical approaches, showing that the computationally cheaper TDDFT can be a suitable option for the simulation of nonlinear spectroscopy of molecular systems that are too large to be treated with higher-level RASSCF/RASPT2 methods.

Self-assembling, structure and nonlinear optical properties of fluorescent organic nanoparticles in water.

Owing to their intense emission, low toxicity and solubility in aqueous medium, fluorescent organic nanoparticles (FONs) have emerged as promising alternatives to inorganic ones for the realization of exogenous probes for bioimaging applications. However, the intimate structure of FONs in solution, as well as the role played by intermolecular interactions on their optical properties, remains challenging to study. Following a recent Second-Harmonic Scattering (SHS) investigation led by two of us [Daniel et al., ACS Photonics, 2015, 2, 1209], we report herein a computational study of the structural organization and second-order nonlinear optical (NLO) properties of FONs based on dipolar chromophores incorporating a hydrophobic triphenylamine electron-donating unit and a slightly hydrophilic aldehyde electron-withdrawing unit at their extremities. Molecular dynamics simulations of the FON formation in water are associated with quantum chemical calculations, to provide insight into the molecular aggregation process, the molecular orientation of the dipolar dyes within the nanoparticles, and the dynamical behavior of their NLO properties. Moreover, the impact of intermolecular interactions on the NLO responses of the FONs is investigated by employing the tight-binding version of the recently developed simplified time-dependent density functional theory (sTD-DFT) approach, allowing the all-atom quantum mechanics treatment of nanoparticles.

Peptoids as a Chiral Stationary Phase for Liquid Chromatography: Insights from Molecular Dynamics Simulations.

Peptoids are peptide regioisomers with attractive structural tunability in terms of sequence and three-dimensional arrangement. Peptoids are foreseen to have a great potential for many diverse applications, including their utilization as a chiral stationary phase in chromatography. To achieve chiral recognition, a chiral side chain is required to allow specific interactions with a given enantiomer from a racemic mixture. One of the most studied chiral stationary phases, built with (S)-N-1-phenylethyl (Nspe) units, was shown to be successful in resolving racemic mixtures of binaphthyl derivatives. However, there is currently no description at the atomic scale of the factors favoring its enantioselectivity. Here, we take advantage of steered molecular dynamics simulations to mimic the elution process at the atomic scale and present evidence that the predominantly right-handed helical conformation of Nspe peptoids and their ability to form stronger hydrogen bonds with the (S) enantiomer are responsible for the chiral recognition of the popular chiral probe 2,2'-bihydroxy-1,1'-binaphthyl.

Thin-Film Engineering of Solution-Processable n-Type Silicon Phthalocyanines for Organic Thin-Film Transistors.

Metal and metalloid phthalocyanines are an abundant and established class of materials widely used in the dye and pigment industry as well as in commercial photoreceptors. Silicon phthalocyanines (SiPcs) are among the highest-performing n-type semiconductor materials in this family when used in organic thin-film transistors (OTFTs) as their performance and solid-state arrangement are often increased through axial substitution. Herein, we study eight axially substituted SiPcs and their integration into solution-processed n-type OTFTs. Electrical characterization of the OTFTs, combined with atomic force microscopy (AFM), determined that the length of the alkyl chain affects device performance and thin-film morphology. The effects of high-temperature annealing and spin coating time on film formation, two key processing steps for fabrication of OTFTs, were investigated by grazing-incidence wide-angle X-ray scattering (GIWAXS) and X-ray diffraction (XRD) to elucidate the relationship between thin-film microstructure and device performance. Thermal annealing was shown to change both film crystallinity and SiPc molecular orientation relative to the substrate surface. Spin time affected film crystallinity, morphology, and interplanar d-spacing, thus ultimately modifying device performance. Of the eight materials studied, bis(tri-n-butylsilyl oxide) SiPc exhibited the greatest electron field-effect mobility (0.028 cm2 V-1 s-1, a threshold voltage of 17.6 V) of all reported solution-processed SiPc derivatives.

Fate of Low-Lying Charge-Transfer Excited States in a Donor:Acceptor Blend with a Large Energy Offset.

In an effort to gain a comprehensive picture of the interfacial states in bulk heterojunction solar cells, we provide a combined experimental-theoretical analysis of the energetics and dynamics of low-lying electronic charge-transfer (CT) states in donor:acceptor blends with a large frontier orbital energy offset. By varying the blend composition and temperature, we unravel the static and dynamic contributions to the disordered density of states (DOS) of the CT-state manifold and assess their recombination to the ground state. Namely, we find that static disorder (conformational and electrostatic) shapes the CT DOS and that fast nonradiative recombination crops the low-energy tail of the distribution probed by external quantum efficiency (EQE) measurements (thereby largely contributing to voltage losses). Our results then question the standard practice of extracting microscopic parameters such as exciton energy and energetic disorder from EQE.

Orientation dependent molecular electrostatics drives efficient charge generation in homojunction organic solar cells.

Organic solar cells usually utilise a heterojunction between electron-donating (D) and electron-accepting (A) materials to split excitons into charges. However, the use of D-A blends intrinsically limits the photovoltage and introduces morphological instability. Here, we demonstrate that polycrystalline films of chemically identical molecules offer a promising alternative and show that photoexcitation of α-sexithiophene (α-6T) films results in efficient charge generation. This leads to α-6T based homojunction organic solar cells with an external quantum efficiency reaching up to 44% and an open-circuit voltage of 1.61 V. Morphological, photoemission, and modelling studies show that boundaries between α-6T crystalline domains with different orientations generate an electrostatic landscape with an interfacial energy offset of 0.4 eV, which promotes the formation of hybridised exciton/charge-transfer states at the interface, dissociating efficiently into free charges. Our findings open new avenues for organic solar cell design where material energetics are tuned through molecular electrostatic engineering and mesoscale structural control.

Surface-Confined Macrocyclization via Dynamic Covalent Chemistry.

Surface-confined synthesis is a promising approach to build complex molecular nanostructures including macrocycles. However, despite the recent advances in on-surface macrocyclization under ultrahigh vacuum, selective synthesis of monodisperse and multicomponent macrocycles remains a challenge. Here, we report on an on-surface formation of [6 + 6] Schiff-base macrocycles via dynamic covalent chemistry. The macrocycles form two-dimensional crystalline domains on the micrometer scale, enabled by dynamic conversion of open-chain oligomers into well-defined ∼3.0 nm hexagonal macrocycles. We further show that by tailoring the length of the alkyl substituents, it is possible to control which of three possible products-oligomers, macrocycles, or polymers-will form at the surface. In situ scanning tunneling microscopy imaging combined with density functional theory calculations and molecular dynamics simulations unravel the synergistic effect of surface confinement and solvent in leading to preferential on-surface macrocyclization.

Highly emissive excitons with reduced exchange energy in thermally activated delayed fluorescent molecules.

Unlike conventional thermally activated delayed fluorescence chromophores, boron-centered azatriangulene-like molecules combine a small excited-state singlet-triplet energy gap with high oscillator strengths and minor reorganization energies. Here, using highly correlated quantum-chemical calculations, we report this is driven by short-range reorganization of the electron density taking place upon electronic excitation of these multi-resonant structures. Based on this finding, we design a series of π-extended boron- and nitrogen-doped nanographenes as promising candidates for efficient thermally activated delayed fluorescence emitters with concomitantly decreased singlet-triplet energy gaps, improved oscillator strengths and core rigidity compared to previously reported structures, permitting both emission color purity and tunability across the visible spectrum.

Application of Rubrene Air-Gap Transistors as Sensitive MEMS Physical Sensors.

Micro-electromechanical systems (MEMS) made of organic materials have attracted efforts for the development a new generation of physical, chemical, and biological sensors, for which the electromechanical sensitivity is the current major concern. Here, we present an organic MEMS made of a rubrene single-crystal air-gap transistor. Applying mechanical pressure on the semiconductor results in high variations in drain current: an unparalleled gauge factor above 4000 has been measured experimentally. Such a high sensitivity is induced by the modulation of charge injection at the interface between the gold electrode and the rubrene semiconductor as an unusual transducing effect. Applying these devices to the detection of acoustic pressure shows that force down to 230 nN can be measured with a resolution of 40 nN. This study demonstrates that MEMS based on rubrene air-gap transistors constitute a step forward in the development of high-performance flexible sensors.

Pentacene Crystal Growth on Silica and Layer-Dependent Step-Edge Barrier from Atomistic Simulations.

Understanding and controlling the growth of organic crystals deposited from the vapor phase is important for fundamental materials science and necessary for applications in pharmaceutical and organic electronics industries. Here, this process is studied for the paradigmatic case of pentacene on silica by means of a specifically tailored computational approach inspired by the experimental vapor deposition process. This scheme is able to reproduce the early stages of the thin-film formation, characterized by a quasi layer-by-layer growth, thus showcasing its potential as a tool complementary to experimental techniques for investigating organic crystals. Crystalline islands of standing molecules are formed at a critical coverage, as a result of a collective reorientation of disordered aggregates of flat-lying molecules. The growth then proceeds by sequential attachment of molecules at the cluster and then terrace edges. Free-energy calculations allowed us to characterize the step-edge barrier for descending the terraces, a fundamental parameter for growth models for which only indirect experimental measurements are available. The barrier is found to be layer-dependent (approximately 1 kcal/mol for the first monolayer on silica, 2 kcal/mol for the second monolayer) and to extend over a distance comparable with the molecular length.

Second-order nonlinear optical properties of Stenhouse photoswitches: insights from density functional theory.

We report the first investigation of the second-order nonlinear optical (NLO) properties of donor-acceptor Stenhouse adducts (DASAs), an emerging class of colored photochromes that undergo photoswitching with visible light to a colorless form. By using time-dependent density functional theory, we provide insights into the relationships linking the nature of the chemical substituents to the amplitude and contrasts of the NLO response. Solvent and frequency dispersion effects are also analyzed. The calculations predict that DASAs behave as high contrast NLO switches, a finding that extends their potential applications to photo-responsive NLO materials and devices.

Deep-Blue Oxadiazole-Containing Thermally Activated Delayed Fluorescence Emitters for Organic Light-Emitting Diodes.

A series of four novel deep-blue to sky-blue thermally activated delayed fluorescence (TADF) emitters (2CzdOXDMe, 2CzdOXD4MeOPh, 2CzdOXDPh, and 2CzdOXD4CF3Ph) have been synthesized and characterized. These oxadiazole-based emitters demonstrated bluer emission compared with the reference emitter 2CzPN thanks to the weaker acceptor strength of the oxadiazole moieties. The oxadiazole compounds doped in hosts (mCP and PPT) emitted from 435 to 474 nm with photoluminescence quantum yields ranging from 14-55%. The emitters possess singlet-triplet excited-state energy gaps (Δ EST) between 0.25 and 0.46 eV resulting in delayed components ranging from 4.8 to 25.8 ms. The OLED device with 2CzdOXD4CF3Ph shows a maximum external quantum efficiency of 11.2% with a sky-blue emission at CIE of (0.17, 0.25), while the device with 2CzdOXD4MeOPh shows a maximum external quantum efficiency of 6.6% with a deep-blue emission at CIE of (0.15, 0.11).

Nonlinear optical responses of self-assembled monolayers functionalized with indolino-oxazolidine photoswitches.

A computational approach combining molecular dynamic simulations and density functional theory (DFT) calculations is implemented to evaluate the second-order nonlinear optical (NLO) responses of photoresponsive self-assembled monolayers (SAMs) based on indolino-oxazolidine molecular switches. These numerical simulations provide a complete atomistic picture of the morphology of the SAMs, revealing a high degree of positional disorder and an almost isotropic orientation of the chromophores. Subsequent DFT calculations, carried out to evaluate the average first hyperpolarizability of indolino-oxazolidine switches within the SAM, predict that the structural disorder does not significantly reduce the NLO contrast compared to that of the isolated molecules. Chromophores in the SAM can assume a limited number of specific conformations, due to the high rotational barrier that characterize the conjugated bonds along the indolino/oxazolidine-dyene-thiophene sequence. A notable exception is the rotation about the thiophene-thioalkyl bond, which is not only almost free, but also strongly correlated with the magnitude of the first hyperpolarizability. Controlling this rotation by chemical design could thus be a viable strategy to optimize the SAMs NLO response and the performance of photoresponsive devices based on indolino/oxazolidine switches.

Multi-dimensional charge transport in supramolecular helical foldamer assemblies.

Aromatic foldamers are bioinspired architectures whose potential use in materials remains largely unexplored. Here we report our investigation of vertical and horizontal charge transport over long distances in helical oligo-quinolinecarboxamide foldamers organized as single monolayers on Au or SiO2. Conductive atomic force microscopy showed that vertical conductivity is efficient and that it displays a low attenuation with foldamer length (0.06 Å-1). In contrast, horizontal charge transport is found to be negligible, demonstrating the strong anisotropy of foldamer monolayers. Kinetic Monte Carlo calculations were used to probe the mechanism of charge transport in these helical molecules and revealed the presence of intramolecular through-space charge transfer integrals approaching those found in pentacene and rubrene crystals, in line with experimental results. Kinetic Monte Carlo simulations of charge hopping along the foldamer chain evidence the strong contribution of multiple 1D and 3D pathways in these architectures and their dependence on conformational order. These findings show that helical foldamer architectures may provide a route for achieving charge transport over long distance by combining multiple charge transport pathways.

Energetic fluctuations in amorphous semiconducting polymers: Impact on charge-carrier mobility.

We present a computational approach to model hole transport in an amorphous semiconducting fluorene-triphenylamine copolymer (TFB), which is based on the combination of molecular dynamics to predict the morphology of the oligomeric system and Kinetic Monte Carlo (KMC), parameterized with quantum chemistry calculations, to simulate hole transport. Carrying out a systematic comparison with available experimental results, we discuss the role that different transport parameters play in the KMC simulation and in particular the dynamic nature of positional and energetic disorder on the temperature and electric field dependence of charge mobility. It emerges that a semi-quantitative agreement with experiments is found only when the dynamic nature of the disorder is taken into account. This study establishes a clear link between microscopic quantities and macroscopic hole mobility for TFB and provides substantial evidence of the importance of incorporating fluctuations, at the molecular level, to obtain results that are in good agreement with temperature and electric field-dependent experimental mobilities. Our work makes a step forward towards the application of nanoscale theoretical schemes as a tool for predictive material screening.