Theoretical study of substitution effects on molecular reorganization energy in organic semiconductors

Periodically Twisted Molecular Wires Based on a Fused Unit for Efficient Intramolecular Hopping Transport

Realizing efficient long-distance intramolecular charge transport based on a hopping mechanism is a key challenge in molecular electronics. In hopping transport, a smaller reorganization energy (λ) and energy difference between hopping sites (ΔEhs) should lead to a smaller activation energy and faster charge transfer. However, the development of π-extended molecules that meet these requirements is challenging. In this study, we successfully synthesized several nanometer-scale π-extended molecules composed of a fused π-conjugated unit as a hopping site for reducing λ. Conformational twists between fused units effectively localize π-conjugation in each unit, contributing to reducing ΔEhs. The expected electronic structures of the oligomers were confirmed using spectroscopic and electrochemical measurements. Single-molecule conductance measurements exhibited higher conductance and lower activation energy than those of nonfused oligothiophenes. First-principles calculations indicated that smaller λ and ΔEhs values explain the high conductance. These results highlight the efficiency of the proposed molecular design for effective intramolecular hopping transport.

Lightening flavin by amination for fluorescent sensing

Monitoring of reactive oxygen species (ROS), such as O2˙-, etc., in organisms is of great significance, not only for their essential role in biological processes, but their excessive production may also result in many diseases. Flavin (FL) is a fluorophore that naturally exists in flavoenzymes, and its fluorescent emission (FE) becomes negligible when reduced. This enables the application of FL derivatives as fluorescent sensors for ROS. We presented a theoretical investigation to address the impact of amino substitution on the photophysical properties of aminoflavins (AmFLs). Resulting from the interplay of electronic and positional effects, amination at C8 enhances the electronic coupling between the ground state and the first singlet excited state by enlarging the adiabatic energy change of the electronic transitions and the emission transition dipole moments, weakens the vibronic coupling by decreasing the contribution of isoalloxazine to the frontier molecular orbitals, redshifts the absorption band, and enhances the fluorescent emission drastically in 8AmFL. The theoretically estimated fluorescent emission intensity of 8AmFL is ∼40 times that of FL, suggesting its potential application as a fluorescent sensor.

Machine-learning-assisted performance improvements for multi-resonance thermally activated delayed fluorescence molecules

With favorable colour purity, multi-resonance thermally activated delayed fluorescence (MR-TADF) molecules exhibit enormous potential in high-definition displays. Due to the relatively small chemical space of MR-TADF molecules, it is challenging to improve molecular performance through domain-specific expertise alone. To address this problem, we focused on optimizing the classic molecule, DABNA-1, using machine learning (ML). Molecular morphing operations were initially employed to generate the adjacent chemical space of DABNA-1. Subsequently, a machine learning model was trained with a limited database and used to predict the properties throughout the generated chemical space. It was confirmed that the top 100 molecules suggested by machine learning present excellent electronic structures, characterized by small reorganization energy and singlet-triplet energy gaps. Our results indicate that the improvement in electronic structures can be elucidated through the view of the molecular orbital (MO). The results also reveal that the top 5 molecules present weaker vibronic peaks of the emission spectrum, demonstrating higher colour purity when compared to DABNA-1. Notably, the M2 molecule presents a high RISC rate, indicating its promising future as a high-efficiency MR-TADF molecule. Our machine-learning-assisted approach facilitates the rapid optimization of classical molecules, addressing a crucial requirement within the organic optoelectronic materials community.

Electric-Field Effects on the Internal Charge Reorganization Energies of Crystalline Organic Semiconductors

The synergistic effects of molecular packing and external electric fields (EEFs, including axial and nonaxial fields) on the internal charge reorganization energies (λ) of typical p-type SMOS have been investigated. Combined quantum and molecular mechanics calculations show that, for all-ring-fused rigid molecules single-molecule approximation and neglect of EEFs are adequate for computing λ, while for nonrigid molecules with inter-ring carbon-carbon (IRCC) linkers, the above simplifications may cause a significant deviation from the actual λ. For nonrigid molecules, solid-state packing can prevent "bad" EEFs (Fz and Fyz) from enhancing λ (adverse to charge transfer), while it allows λ to be greatly reduced (in favor of charge transfer) if "good" EEFs (Fx, Fxy, Fxz and Fxyz) are imposed. Last, a simple strategy that can divide λ into each subring contribution for IRCC-linked molecules has been proposed.

Intrinsic Photostability in Dithiolonaphthalimide Achieved by Disulfide Bond-Induced Excited-State Quenching

Disulfide bridges common in proteins show excellent photostability achieved by ultrafast internal conversion and maintain the stability of the tertiary structure. When disulfide bonds exist in aromatic compounds, the rigid chemical structure may affect the cleavage and reforming dynamics of disulfide bonds. In this work, a model compound with a disulfide five-membered-ring structure, 4,5-dithiolo-N-(2,6-dimethylphenyl)-1,8-naphthalimide (DTDPNI), is selected to elaborate the effect of disulfide modification on the excited-state deactivation mechanism. Quantum chemical calculations show that the S-S stretching leads to a dramatic decrease in the energy gap between the S1 and S0 states, similar to the situation in 1,2-dithiane. Due to the efficient nonradiative process, the excited-state lifetime of DTDPNI resolved by ultrafast spectroscopy is determined to be ∼20 ps. It is found that the excellent photostability is achieved by ultrafast excited-state quenching induced by the S-S stretching, rather than the cleavage of the disulfide bond; even the disulfide bridge is in a very rigid aromatic molecular system.

Computational studies on nitrogen (N)-substituted 2,6-diphenylanthracene: a novel precursor of organic field effect transistor materials

The investigation shows that nitrogen (N)-substituted π-conjugated semiconductor materials have improved optical and electronic performance and work efficiently in organic field-effect transistors (OFETs).

Rational design of organic semiconductors with low internal reorganization energies for hole and electron transport: position effect of aza-substitution in phenalenyl derivatives

Amphoteric-redox phenalenyl radical (PLY) is a suitable candidate used to design ambipolar organic materials. Because the singly occupied nonbonding molecular orbital (NBMO) of PLY has a perfect local nonbonding character, its internal reorganization energy (λ) for transporting holes (λ⁺) or electrons (λ^-) is known to be small. Herein, PLY is employed to study the position effect of the aza group on the λ. By adding or extracting an electron from the NBMO, the bond length alterations can be minute. Therefore, the PLY derivatives are also an excellent candidate to study the contributions from the bond angle alterations to the λ. Substituting the aza groups at the β- or α-positions of PLY shows two different trends. When consecutively substituting the aza group at the three β-positions of PLY, the λs are consistently decreased. Contrarily, a series of double functionalization of aza groups at the four α-positions of PLY, the λs are increased. It is because the local bonding or antibonding character in frontier orbitals (FMO) is observed in α2N-PLY and α4N-PLY. As the FMOs of the three β-substituted PLYs and α6N-PLY have perfect local nonbonding character, we found the bond angle alterations are the main contributors of λ. The λs for most aza-PLYs were smaller than 100 meV. Thus, we propose a design rule for substituting aza groups on the parent molecules with strong local nonbonding character in their FMOs. Based on the adiabatic ionization potential and electron affinity, two π-extended PLY derivatives with small λ were recommended for fabricating air-stable ambipolar OFET.

Active discovery of organic semiconductors

The versatility of organic molecules generates a rich design space for organic semiconductors (OSCs) considered for electronics applications. Offering unparalleled promise for materials discovery, the vastness of this design space also dictates efficient search strategies. Here, we present an active machine learning (AML) approach that explores an unlimited search space through consecutive application of molecular morphing operations. Evaluating the suitability of OSC candidates on the basis of charge injection and mobility descriptors, the approach successively queries predictive-quality first-principles calculations to build a refining surrogate model. The AML approach is optimized in a truncated test space, providing deep methodological insight by visualizing it as a chemical space network. Significantly outperforming a conventional computational funnel, the optimized AML approach rapidly identifies well-known and hitherto unknown molecular OSC candidates with superior charge conduction properties. Most importantly, it constantly finds further candidates with highest efficiency while continuing its exploration of the endless design space.

FB-REDA: fragment-based decomposition analysis of the reorganization energy for organic semiconductors

We present a fragment-based decomposition analysis tool (FB-REDA) for the reorganisation energy (λ). This tool delivers insights on how to rationally design low-λ organic semiconductors. The contribution of the fragment vibrational modes to the reorganization energy is exploited to identity the individual contributions of the molecular building blocks. The usefulness of the approach is demonstrated by offering three strategies to reduce the reorganization energy of a promising dopant-free hole transport material (TPA1PM, λ = 213 meV). A reduction of nearly 50% (TPD3PM, λ = 108 meV) is achieved. The proposed design principles are likely transferable to other organic semiconductors exploiting common molecular building blocks.

Molecular vibrations reduce the maximum achievable photovoltage in organic solar cells

The low-energy edge of optical absorption spectra is critical for the performance of solar cells, but is not well understood in the case of organic solar cells (OSCs). We study the microscopic origin of exciton bands in molecular blends and investigate their role in OSCs. We simulate the temperature dependence of the excitonic density of states and low-energy absorption features, including low-frequency molecular vibrations and multi-exciton hybridisation. For model donor-acceptor blends featuring charge-transfer excitons, our simulations agree very well with temperature-dependent experimental absorption spectra. We unveil that the quantum effect of zero-point vibrations, mediated by electron-phonon interaction, causes a substantial exciton bandwidth and reduces the open-circuit voltage, which is predicted from electronic and vibronic molecular parameters. This effect is surprisingly strong at room temperature and can substantially limit the OSC's efficiency. Strategies to reduce these vibration-induced voltage losses are discussed for a larger set of systems and different heterojunction geometries.

A computational study of anisotropic charge transport in air-stable fluorinated benzobisbenzothiophene (FBBBT) derivatives

A computational study of anisotropical charge transport properties of fluorinated benzobisbenzothiohphene derivatives (FBBBT) is presented. The values of IP_adia of all FBBBTs are found in the range of 6.00-6.20 eV inferring the fact that the investigated compounds have ambient air-stability. In addition, the energy levels of FBBBT s are found to be lower than those of benzobisbenzothiophene (BBBT) compound indicating higher charge carrier stability in the former. Hirshfield surface analyses showed that, in all the studied compounds, the principal identifiable interaction were mostly due to F⋯H and H⋯H intermolecular couplings with no contribution from S⋯S bondings. The calculated maximum μ_hole(μ_elec) value of the compounds FBBBT-a and FBBBT-b was found to be 0.483 (0.794) cm²V^- 1s^- 1 and 0.688 (0.542) cm²V^- 1s^- 1 respectively in the direction of transistor channel (Φ = 93.39 ^∘(273.30^∘) for FBBBT-a and Φ = 92.24 ^∘/272.72 ^∘ for FBBBT-b). For FBBBT-c, the maximum μ_elec(μ_hole) value of 0.933 (0.233) cm²V^- 1s^- 1 appeared for Φ = 0 ^∘/179.90 ^∘. In addition, the compounds FBBBT-a and FBBBT-b possess two additional fluorine atoms attached at the X positions in the backbone, which result in an increment in μ_elec values (1.4 times and 0.78 times higher than μ_hole) in these two compounds at a particular crystal direction.

Design of singlet fission chromophores with cyclic (alkyl)(amino) carbene building blocks

We use MRSF-TDDFT and NEVPT2 methods to design singlet fission chromophores with the building blocks of cyclic (alkyl)(amino)carbenes (CAACs). CAAC dimers with C₂, C₄, and p-phenylene spacers are considered. The substitutions with trifluoromethyls and fluorine atoms at the α C position are investigated. The electronegative substituents enhance the π accepting capability of the α C while maintaining it as a quaternary C atom. The phenylene-connected dimers with the two substitutions are identified as promising candidates for singlet fission chromophores. The cylindrically symmetric C₂ and C₄ spacers allow for substantial structural reorganizations in the S₀-to-S₁ and S₀-to-T₁ excitations. Although the two substituted dimers with the C₄ spacer satisfy (or very close to satisfy) the primary thermodynamics criterion for singlet fission, the significant structural reorganizations result in high barriers so that the fission is kinetically unfavorable.

Deciphering the Multifarious Charge-Transport Behaviour of Crystalline Propeller-Shaped Triphenylamine Analogues

A collection of para-substituted propeller-shaped triphenylamine (TPA) derivatives have been computationally investigated for charge-transport characteristics exhibited by the derivatives by using the Marcus-Hush formalism. The various substituents chosen herein, with features that range from electron withdrawing to electron donating in nature, play a key role in defining the reorganisation energy and electronic coupling properties of the TPA derivatives. The TPA moiety is expected to possess weak electronic coupling on the basis of poor orbital overlap upon aggregation, owing to the restriction imposed by the propeller shape of the TPA core. However, the substituent groups attached to the TPA core can significantly dictate the crystal-packing motif of the TPA derivatives, wherein the variety of noncovalent intermolecular interactions subsequently generated drive the packing arrangement and influence electronic coupling between the neighbouring orbitals. Intermolecular interactions in the crystalline architecture of TPA derivatives were probed by using Hirshfeld and quantum theory of atoms-in-molecules techniques. Furthermore, symmetry-adapted perturbation theory analysis of the TPA analogues has revealed that a periodic arrangement of energetically stable dimers with significant electronic coupling is essential to contribute high charge-carrier mobility to the overall crystal.

A quantitative structure-property study of reorganization energy for known p-type organic semiconductors

Intramolecular reorganization energy (RE), which quantifies the electron-phonon coupling strength, is an important charge transport parameter for the theoretical characterization of molecular organic semiconductors (OSCs). On a small scale, the accurate calculation of the RE is trivial; however, for large-scale screening, faster approaches are desirable. We investigate the structure-property relations and present a quantitative structure-property relationship study to facilitate the computation of RE from molecular structure. To this end, we generated a compound set of 171, which was derived from known p-type OSCs built from moieties such as acenes, thiophenes, and pentalenes. We show that simple structural descriptors such as the number of atoms, rings or rotatable bonds only weakly correlate with the RE. On the other hand, we show that regression models based on a more comprehensive representation of the molecules such as SMILES-based molecular signatures and geometry-based molecular transforms can predict the RE with a coefficient of determination of 0.7 and a mean absolute error of 40 meV in the library, in which the RE ranges from 76 to 480 meV. Our analysis indicates that a more extensive compound set for training is necessary for more predictive models.

Theoretical investigations into the charge transfer properties of thiophene α-substituted naphthodithiophene diimides: excellent n-channel and ambipolar organic semiconductors

A theoretical study was carried out to investigate the electronic structures and the charge transport properties of a series of naphthodithiophene diimide (NDTI) thiophene α-substituted derivatives NDTI-X using density functional theory and classical Marcus charge transfer theory. This study deeply revealed the structure-property relationships by analyzing the intermolecular interactions in crystal structures of C8-NDTI and C8-NDTI-Cl thoroughly by using the Hirshfeld surface, QTAIM theories and symmetry-adapted perturbation theory (SAPT). Our results suggested that a 2-D brick-like π-stacking structure makes C8-NDTI-Cl a more excellent n-type semiconducting material with μ_max-e of 2.554 cm² V^-1 s^-1 than C8-NDTI with a herringbone-like slipped π-stacking motif. In addition, the calculated results showed that by modifying the thiophene α-positions of NDTI with electron-withdrawing substituents, -F, -Cl and -CN, low-lying LUMO energy levels and a high adiabatic electron affinity EA(a) can be obtained; while introducing electron-donating groups, benzene (-B), thiophene (-T), benzo[b]thiophene (-BT) and naphtha[2,3-b]thiophene (-NT), expanded the molecular π-conjugated backbone, and narrow band gaps, high EA(a) and small reorganization energies can be obtained. Theoretical simulations predict that NDTI-CN is an excellent air-stable n-type organic semiconducting material with an average electron mobility μ_e of up to 1.743 cm² V^-1 s^-1. Owing to their high EA(a), moderate adiabatic ionization potential IP(a) as well as small hole and electron reorganization energies, NDTI-BT and NDTI-NT are two well-balanced air-stable ambipolar semiconducting materials. The theoretical average hole/electron mobilities are as high as 2.708/3.739 cm² V^-1 s^-1 for C8-NDTI-NT and 1.597/2.350 cm² V^-1 s^-1 for C8-NDTI-BT, respectively.

Theoretical study on electron structure and charge transport properties of tetraazapentacene derivatives

By Means of Marcus electron transfer theory, the charge transport properties of tetraazapentacene (4N-PEN) derivatives were systematically explored. The reorganization energies were studied by both adiabatic potential-energy surfaces and normal mode analysis. The charge diffusion constants were evaluated from the random walk simulation. From the perspective of homology modeling, a selected 4N-PEN derivative without experimental crystal structure was built into three kinds of possible packing modes with reference to its relative analogues and then fully optimized. The calculated results show that the charge transport property for the same kind of systems strongly depends on the packing mode, and the π···stacking is more beneficial for electron transport of 4N-PEN derivatives. Meanwhile, the 4N-PEN derivatives have larger electron transfer integrals and lower energy levels of the lowest unoccupied molecular orbitals as well as smaller electron reorganization energies, which provides a three-in-one advantage for electron transport. Fascinatingly, the data obtained from the hopping and band models both suggest that the 4N-PEN derivatives have the intrinsic property of electron transport. Thus, the 4N-PEN derivatives have the potential for preparing n-type organic semiconductors.

Effect of five-membered ring and heteroatom substitution on charge transport properties of perylene discotic derivatives: A theoretical approach

Density functional theory calculations were carried out to investigate the evolvement of charge transport properties of a set of new discotic systems as a function of ring and heteroatom (B, Si, S, and Se) substitution on the basic structure of perylene. The replacement of six-membered rings by five-membered rings in the reference compound has shown a prominent effect on the electron reorganization energy that decreases ∼0.2 eV from perylene to the new carbon five-membered ring derivative. Heteroatom substitution with boron also revealed to lower the LUMO energy level and increase the electron affinity, therefore lowering the electron injection barrier compared to perylene. Since the rate of the charge transfer between two molecules in columnar discotic systems is strongly dependent on the orientation of the stacked cores, the total energy and transfer integral of a dimer as a disc is rotated with respect to the other along the stacking axis have been predicted. Aimed at obtaining a more realistic approach to the bulk structure, the molecular geometry of clusters made up of five discs was fully optimized, and charge transfer rate and mobilities were estimated for charge transport along a one dimensional pathway. Heteroatom substitution with selenium yields electron transfer integral values ∼0.3 eV with a relative disc orientation of 25°, which is the preferred angle according to the dimer energy profile. All the results indicate that the tetraselenium-substituted derivative, not synthetized so far, could be a promising candidate among those studied in this work for the fabrication of n-type semiconductors based on columnar discotic liquid crystals materials.

From charge transport parameters to charge mobility in organic semiconductors through multiscale simulation

This review introduces the development and application of a multiscale approach to assess the charge mobility for organic semiconductors, which combines quantum chemistry, Kinetic Monte Carlo (KMC), and molecular dynamics (MD) simulations. This approach is especially applicable in describing a large class of organic semiconductors with intermolecular electronic coupling (V) much less than intramolecular charge reorganization energy (λ), a situation where the band description fails obviously. The charge transport is modeled as successive charge hopping from one molecule to another. We highlight the quantum nuclear tunneling effect in the charge transfer, beyond the semiclassical Marcus theory. Such an effect is essential for interpreting the "paradoxical" experimental finding that optical measurement indicated "local charge" while electrical measurement indicated "bandlike". Coupled MD and KMC simulations demonstrated that the dynamic disorder caused by intermolecular vibration has negligible effect on the carrier mobility. We further apply the approach for molecular design of n-type materials and for rationalization of experimental results. The charge reorganization energy is analyzed through decomposition into internal coordinates relaxation, so that chemical structure contributions to the intramolecular electron-phonon interaction are revealed and give helpful indication to reduce the charge reorganization energy.

DFT study of vibronic properties of d8 (Ni-, Pd-, and Pt-) phthalocyanines

By means of density functional theory, we have studied the electronic structure and vibronic properties of single neutral NiPc, PdPc, and PtPc molecules and their singly and doubly ionized cations and anions. In particular, the vibronic couplings and reorganization energies of all systems are compared. Partitioning of the reorganization energy, corresponding to the photoelectron spectra of the first and second ionizations of studied molecules, into normal mode contributions shows that the major contributions are due to several vibrational modes with a(1g) symmetry and energies lower than 1600 cm(-1). The results reveal that the reorganization energy due to the singly positive ionization in the studied molecules is up to about one order of magnitude less than other reorganization energies. This makes these metal phthalocyanines, from the perspective of intramolecular reorganization energies, attractive as electron donor for intramolecular electron transfer in electron acceptor-donor systems.