On combining the conductor-like screening model and optimally tuned range-separated hybrid density functionals

Screened optimally tuned range separated hybrid functional for solvated low bandgap molecular systems

The description of electronic properties of low bandgap molecular system is often performed by using density functional theory (DFT) and time dependent (TD) DFT calculations with the optimally tuned range-separated hybrid (OT-RSH) functional, as it contains the necessary ingredients to reliably predict charge transfer excitations. However, the range separating parameter (ω) is system-dependent and its optimization, including the chemical environment, is intricate. Refaely-Abramson et al. demonstrated that the gap renormalization in molecular crystals, a ground state property, can be represented by an OT-RSH functional screened by ɛstatic [Phys. Rev. B 88, 081204(R) (2013)], the zero frequency scalar dielectric constant. In this study, we propose the use of an OT-RSH functional screened by the scalar dielectric constant in the high frequency limit (OT-sRSH), ɛ∞, an appropriate constraint for vertical ionization energies or excitations in a dielectric environment. We have performed calculations for S,N-heteroacene derivatives in tetrahydrofuran and dichloromethane. The "unscreened" OT-RSH functional tends to underestimate experimental ionization potentials (IPs) and optical gaps (Egs) by up to 1.5 and 0.5 eV, respectively. In contrast, OT-sRSH functional calculations underestimate IPs and Egs by only 0.4 and 0.2 eV. We also compared the OT-sRSH results to explicitly solvated OT-RSH functional calculations for oligothiophenes in dioxane, benzene in ammonia, and methylene blue in water. We observe that both the approaches perform similarly for weakly interacting intermolecular systems and deviate for solvent-solute interacting systems, as expected. In conclusion, the OT-sRSH functional can describe molecular systems with environmental polarization effects accurately, a step toward describing realistic molecular systems.

Multistate, Polarizable QM/MM Embedding Scheme Based on the Direct Reaction Field Method: Solvatochromic Shifts, Analytical Gradients and Optimizations of Conical Intersections in Solution

We recently introduced a polarizable embedding scheme based on an integral-exact reformulation of the direct reaction field method (IEDRF) that accounts for the differential solvation of ground and excited states in QM/MM simulations. The polarization and dispersion interactions between the quantum-mechanical (QM) and molecular-mechanical (MM) regions are described by the DRF Hamiltonian, while the Pauli repulsion between explicitly treated QM electrons and the implicit electron density around MM atoms is modeled with effective core potentials. A single Hamiltonian is used for all electronic states so that Born-Oppenheimer states belonging to the same geometry are orthogonal and state crossings are well-defined. In this work, we describe the implementation of the method using graphical processing unit acceleration in TeraChem, where it is combined with multiple electronic structure methods, including Hartree-Fock, time-dependent density functional theory, and complete active space self-consistent field. In contrast with older implementations of the DRF method, integrals of the polarization operators are evaluated exactly. Expressions for ingredients needed to construct analytical gradients and nonadiabatic coupling vectors are derived and tested by optimizing a conical intersection between two excited states in the presence of a polarizable solvent shell. The method is applied to estimate the solvent shifts of absorption energies of a series of donor-acceptor dyes having low-lying charge-transfer states. Even for a nonpolar solvent such as n-hexane, the inclusion of its static polarizability leads to non-negligible shifts that improve the agreement to essentially quantitative levels (0.03 eV) with full-system calculations. Good agreement with the positions of the experimental absorption maxima measured in solution is also observed.

Tailoring the Weight of Surface and Intralayer Edge States to Control LUMO Energies

The energies of the frontier molecular orbitals determine the optoelectronic properties in organic films, which are crucial for their application, and strongly depend on the morphology and supramolecular structure. The impact of the latter two properties on the electronic energy levels relies primarily on nearest-neighbor interactions, which are difficult to study due to their nanoscale nature and heterogeneity. Here, an automated method is presented for fabricating thin films with a tailored ratio of surface to bulk sites and a controlled extension of domain edges, both of which are used to control nearest-neighbor interactions. This method uses a Langmuir-Schaefer-type rolling transfer of Langmuir layers (rtLL) to minimize flow during the deposition of rigid Langmuir layers composed of π-conjugated molecules. Using UV-vis absorption spectroscopy, atomic force microscopy, and transmission electron microscopy, it is shown that the rtLL method advances the deposition of multi-Langmuir layers and enables the production of films with defined morphology. The variation in nearest-neighbor interactions is thus achieved and the resulting systematically tuned lowest unoccupied molecular orbital (LUMO) energies (determined via square-wave voltammetry) enable the establishment of a model that functionally relates the LUMO energies to a morphological descriptor, allowing for the prediction of the range of accessible LUMO energies.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Bifacial Dye Membranes: Ultrathin and Free-Standing although not Being Covalently Bound

Layers of aligned dyes are key to photo-driven charge separation in dye sensitized solar cells, but cannot be exploited as rectifying membranes in photocatalysis to separate half-cells because they are not sufficiently stable. While impressive work on the fabrication of stable noncovalent membranes has been recently demonstrated, these membranes are inherently suffering from non-uniform orientation of the constituting dyes. To stabilize layers made from uniformly assembled and aligned dyes, they can be covalently cross-linked via functional groups or via chromophores at the expense of their optical properties. Here stable membranes from established dyes are reported that do not need to be elaborately functionalized nor do their chromophores need to be destroyed. These membranes are free-standing, although being only non-covalently linked. To enable uniform dye-alignment, Langmuir layers made from linear, water-insoluble dyes are used. That water-soluble charge transfer dyes adsorb onto and intercalate into the Langmuir layer from the aqueous subphase, thus yielding free-standing, molecularly thin membranes are demonstrated. The developed bifacial layers consist almost entirely of π-conjugated units and thus can conduct charges and can be further engineered for optoelectronic and photocatalytic applications.

Langmuir-Blodgett Films of Diketopyrrolopyrroles with Tunable Amphiphilicity

In this work, we present the formation of H- and J-aggregates of amphiphilic centrosymmetric diketopyrrolopyrroles containing aliphatic or aromatic amino groups. The inherent amphiphilicity of these dyes predestines their assembly at interfaces to form ordered supramolecular structures. In this work, we employed the Langmuir-Blodgett (LB) technique to generate, manipulate, and deposit such supramolecular structures. The aforementioned amines provide an additional means to control the formation of the supramolecular assemblies. In the resulting LB films, both H- and J-aggregates of the dyes can be realized, leading to very broad absorption spectra. In contrast to many reports on H- and J-aggregates, the interactions between the symmetric diketopyrrolopyrroles are controlled via interface assembly and π-stacking and not by dipolar interactions. We show that in the case of the aliphatic, but not for the aromatic amine functionalization, the usage of an acidic subphase enables the transition from H- to J-aggregate-dominated LB films via an increase in the surface pressure during deposition.

Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory

The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.

Charge-transfer electronic states in organic solar cells: a TDDFT study

The prediction of new organic photovoltaic materials in organic solar cells (OSCs) must include a precise description of charge-transfer states because they are involved in electron-transfer processes such as charge separation and charge recombination which govern the device efficiency. Also, as the experimental performance of an optoelectronic device is measured for nonequilibrium nanostructures, computational approaches need models that can incorporate morphology effects. Usually, this aspect is treated by molecular dynamics simulation (MDS) methodologies; however, methodologies and formalisms to calculate the electron-transfer processes are still controversial and sometimes do not connect their information with the phase morphologies. In this work we propose a simple and fast characterization of electron-transfer processes to find the rate constants by analysing the distribution of vertical excitation energies of both local excitation (LE) and charge-transfer (CT) states using TD-DFT calculations in the donor-acceptor pair structures which were extracted from MDS. This proposal assumes that conformational changes are prevented and equilibria are not achieved while the electron-transfer events take effect, and thus the only pathway that connects the LE and CT states is their surface crossing point where an ideal distribution might exist. Different density functionals and dialectric models were tested. The results indicate a close relationship between the proposal and experimental data for electron-transfer events, suggesting the application of this method in the rational design of new photovoltaic materials.

Macroscopic quantum electrodynamics and density functional theory approaches to dispersion interactions between fullerenes

The processing and material properties of commercial organic semiconductors, for e.g. fullerenes is largely controlled by their precise arrangements, specially intermolecular symmetries, distances and orientations, more specifically, molecular polarisabilities. These supramolecular parameters heavily influence their electronic structure, thereby determining molecular photophysics and therefore dictating their usability as n-type semiconductors. In this article we evaluate van der Waals potentials of a fullerene dimer model system using two approaches: (a) Density Functional Theory and, (b) Macroscopic Quantum Electrodynamics, which is particularly suited for describing long-range van der Waals interactions. Essentially, we determine and explain the model symmetry, distance and rotational dependencies on binding energies and spectral changes. The resultant spectral tuning is compared using both methods showing correspondence within the constraints placed by the different model assumptions. We envision that the application of macroscopic methods and structure/property relationships laid forward in this article will find use in fundamental supramolecular electronics.