Molecular dynamics in low-spin excited states

Benchmarking Charge-Transfer Excited States in TADF Emitters: ΔDFT Outperforms TD-DFT for Emission Energies

Charge-transfer (CT) excited states are crucial to organic light-emitting diodes (OLEDs), particularly to those based on thermally activated delayed fluorescence (TADF). However, accurately modeling CT states remains challenging, even with modern implementations of (time-dependent) density functional theory [(TD-)DFT], especially in a dielectric environment. To identify shortcomings and improve the methodology, we previously established the STGABS27 benchmark set with highly accurate experimental references for the adiabatic energy gap between the lowest singlet and triplet excited states (ΔEST). Here, we diversify this set to the STGABS27-EMS benchmark by including experimental emission energies (Eem) and use this new set to (re)-evaluate various DFT-based approaches. Surprisingly, these tests demonstrate that a state-specific (un)restricted open-shell Kohn-Sham (U/ROKS) DFT coupled with a polarizable continuum model for perturbative state-specific nonequilibrium solvation (ptSS-PCM) provides exceptional accuracy for predicting Eem over a wide range of density functionals. In contrast, the main workhorse of the field, Tamm-Dancoff-approximated TD-DFT (TDA-DFT) paired with the same ptSS-PCM, is distinctly less accurate and strongly functional-dependent. More importantly, while TDA-DFT requires the choice of two very different density functionals for good performance on either ΔEST or Eem, the time-independent U/ROKS/PCM approaches deliver excellent accuracy for both quantities with a wide variety of functionals.

ΔDFT Predicts Inverted Singlet-Triplet Gaps with Chemical Accuracy at a Fraction of the Cost of Wave Function-Based Approaches

Efficient OLEDs need to quickly convert singlet and triplet excitons into photons. Molecules with an inverted singlet-triplet energy gap (INVEST) are promising candidates for this task. However, typical INVEST molecules have drawbacks like too low oscillator strengths and excitation energies. High-throughput screening could identify suitable INVEST molecules, but existing methods are problematic: The workhorse method TD-DFT cannot reproduce gap inversion, while wave function-based methods are too slow. This study proposes a state-specific method based on unrestricted Kohn-Sham DFT with common hybrid functionals. Tuned on the new INVEST15 benchmark set, this method achieves an error of less than 1 kcal/mol, which is traced back to error cancellation between spin contamination and dynamic correlation. Applied to the larger and structurally diverse NAH159 set in a black-box fashion, the method maintains a small error (1.2 kcal/mol) and accurately predicts gap signs in 83% of cases, confirming its robustness and suitability for screening workflows.

MiMiC: A high-performance framework for multiscale molecular dynamics simulations

MiMiC is a framework for performing multiscale simulations in which loosely coupled external programs describe individual subsystems at different resolutions and levels of theory. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data (MPMD) paradigm, serving as an intermediary responsible for fast data exchange and interactions between the subsystems. The main goal of MiMiC is to avoid interfering with the underlying parallelization of the external programs, including the operability on hybrid architectures (e.g., CPU/GPU), and keep their setup and execution as close as possible to the original. At the moment, MiMiC offers an efficient implementation of electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) that has demonstrated unprecedented parallel scaling in simulations of large biomolecules using CPMD and GROMACS as QM and MM engines, respectively. However, as it is designed for high flexibility with general multiscale models in mind, it can be straightforwardly extended beyond QM/MM. In this article, we illustrate the software design and the features of the framework, which make it a compelling choice for multiscale simulations in the upcoming era of exascale high-performance computing.

Orbital-Optimized Versus Time-Dependent Density Functional Calculations of Intramolecular Charge Transfer Excited States

The performance of time-independent, orbital-optimized calculations of excited states is assessed with respect to charge transfer excitations in organic molecules in comparison to the linear-response time-dependent density functional theory (TD-DFT) approach. A direct optimization method to converge on saddle points of the electronic energy surface is used to carry out calculations with the local density approximation (LDA) and the generalized gradient approximation (GGA) functionals PBE and BLYP for a set of 27 excitations in 15 molecules. The time-independent approach is fully variational and provides a relaxed excited state electron density from which the extent of charge transfer is quantified. The TD-DFT calculations are generally found to provide larger charge transfer distances compared to the orbital-optimized calculations, even when including orbital relaxation effects with the Z-vector method. While the error on the excitation energy relative to theoretical best estimates is found to increase with the extent of charge transfer up to ca. -2 eV for TD-DFT, no correlation is observed for the orbital-optimized approach. The orbital-optimized calculations with the LDA and the GGA functionals provide a mean absolute error of ∼0.7 eV, outperforming TD-DFT with both local and global hybrid functionals for excitations with a long-range charge transfer character. Orbital-optimized calculations with the global hybrid functional B3LYP and the range-separated hybrid functional CAM-B3LYP on a selection of states with short- and long-range charge transfer indicate that inclusion of exact exchange has a small effect on the charge transfer distance, while it significantly improves the excitation energy, with the best-performing functional CAM-B3LYP providing an absolute error typically around 0.15 eV.

NBN- and BNB-Phenalenyls: the Yin and Yang of Heteroatom-doped π Systems

NBN- and BNB-doped phenalenyls are isoelectronic to phenalenyl anions and cations, respectively. They represent a pair of complementary molecules that have essentially identical structures but opposite properties as electron donors and acceptors. The NBN-phenalenyls 1-4 considered here were prepared from N,N'-dimethyl-1,8-diaminonaphthalene and readily available boron-containing building blocks (i. e., BH3⋅SMe2 (1), p-CF3-C6H4B(OH)2 (2), C6H5B(OH)2 (3), or MesBCl2/iPr2NEt (4)). Treatment of 1 with 4-Me2N-2,6-Me2-C6H2Li gave the corresponding NBN derivative 5. The BNB-phenalenyl 6 was synthesized from 1,8-naphthalenediyl-bridged diborane(6), PhNH2, and MesMgBr. A computational study reveals that the photoemission of 1, 4, and 5 originates from locally excited (LE) states at the NBN-phenalenyl fragments, while that of 2 is dominated by charge transfer (CT) from the NBN-phenalenyl to the p-CF3-C6H4 fragment. Depending on the dihedral angle θ between its Ph and NBN planes, compound 3 emits mainly from a less polar LE (θ >55°) or more polar CT state (θ <55°). In turn, the energetic preference for either state is governed by the polarity of the solvent used. An equimolar aggregate of the NBN- and BNB-phenalenyls 3 and 6 (in THF/H2O) shows a distinct red-shifted emission compared to that of the individual components, which originates from an intermolecular CT state.

Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions

X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn-Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. We also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O3. This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

Probing C-I bond fission in the UV photochemistry of 2-iodothiophene with core-to-valence transient absorption spectroscopy

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C4H3IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Gradual absorption changes up to ∼0.5 to 1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30 fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C-I bond, enabling a nonadiabatic transition to a dissociative π→σC-I* state. For the subsequent fragmentation, a relatively narrow bond-length region along the C-I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments.

Nonadiabatic Derivative Couplings Calculated Using Information of Potential Energy Surfaces without Wavefunctions: Ab Initio and Machine Learning Implementations

In this work, we implemented an approximate algorithm for calculating nonadiabatic coupling matrix elements (NACMEs) of a polyatomic system with ab initio methods and machine learning (ML) models. Utilizing this algorithm, one can calculate NACMEs using only the information of potential energy surfaces (PESs), i.e., energies, and gradients as well as Hessian matrix elements. We used a realistic system, namely CH₂NH, to compare NACMEs calculated by this approximate PES-based algorithm and the accurate wavefunction-based algorithm. Our results show that this approximate PES-based algorithm can give very accurate results comparable to the wavefunction-based algorithm except at energetically degenerate points, i.e., conical intersections. We also tested a machine learning (ML)-trained model with this approximate PES-based algorithm, which also supplied similarly accurate NACMEs but more efficiently. The advantage of this PES-based algorithm is its significant potential to combine with electronic structure methods that do not implement wavefunction-based algorithms, low-scaling energy-based fragment methods, etc., and in particular efficient ML models, to compute NACMEs. The present work could encourage further research on nonadiabatic processes of large systems simulated by ab initio nonadiabatic dynamics simulation methods in which NACMEs are always required.

Eco-friendly inorganic molecular novel antiperovskites for light-emitting application

The development of perovskite light-emitting diodes (PeLEDs) has progressed rapidly over the past several years, with high external quantum efficiencies exceeding 20%. However, the deployment of PeLEDs in commercial devices still faces severe challenges, such as environmental pollution, instability and low photoluminescence quantum yields (PLQYs). In this work, we perform high-throughput calculations to exhaustively search the unexplored and eco-friendly novel antiperovskite space (formula: X₃B[MN₄], with octahedron [BX₆] and tetrahedron [MN₄]). The novel antiperovskites have a unique structure whereby a tetrahedron can be embedded into an octahedral skeleton as a light-emitting center causing a space confinement effect, leading to the characteristics of a low-dimensional electronic structure, which then makes these materials potential light-emitting material candidates with a high PLQY and light-emitting stability. Under the guidance of newly derived tolerance, octahedral, and tetrahedral factors, 266 stable candidates are successfully screened out from 6320 compounds. Moreover, the antiperovskite materials Ba₃I_0.5F_0.5(SbS₄), Ca₃O(SnO₄), Ba₃F_0.5I_0.5(InSe₄), Ba₃O_0.5S_0.5(ZrS₄), Ca₃O(TiO₄), and Rb₃Cl_0.5I_0.5(ZnI₄) possess an appropriate bandgap, thermodynamic and kinetic stability, and excellent electronic and optical properties, making them promising light-emitting materials.

Exciton Nature of Plasma Phase Transition in Warm Dense Fluid Hydrogen: ROKS Simulation

The transition of warm dense fluid hydrogen from an insulator to a conducting state at pressures of about 20-400 GPa and temperatures of 500-5000 K has been the subject of active scientific research over the past few decades. However, various experimental and theoretical methods do not provide consistent results. In this work, we have applied the restricted open-shell Kohn-Sham (ROKS) method for first principles molecular dynamics of dense hydrogen after thermal excitation to the first singlet excited state. The Wannier localization method has allowed us to analyze the exciton dynamics in this system. The model shows that a key mechanism of the transition is associated with the dissociation of electron-hole pairs, which allows explaining several stages of the transition of fluid H₂ from molecular state to plasma. This mechanism is able to give a quantitative description of several experimental results as well as to resolve the discrepancies between experimental studies.

Symmetry of Identical Particles, Modern Achievements in the Pauli Exclusion Principle, in Superconductivity and in Some Other Phenomena

In this review, the modern achievements in studies of the Pauli exclusion principle (PEP) and the properties of the identical particle systems when PEP is not fulfilled are discussed. The validity of conception of the spin in the framework of density functional theory (DFT) approaches is analyzed. The modern state of the recently discovered Fe-based superconductors is discussed in detail. These materials belong to the paramagnetic semimetal family and become superconductors upon doping. Recently, in 2020, room-temperature superconductivity was realized. However, from the following discussion in the SC community, it was not evident that the results of room-temperature superconductivity have been repeated by other laboratories. Thus, the question “is room temperature really achieved?” is still open. In the concluding remarks, we present the explanation of why the PEP limitations on the symmetry of identical particles system exist in nature, and following from it, some important consequences.

Birth of the Hydrated Electron via Charge-Transfer-to-Solvent Excitation of Aqueous Iodide

A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I^-(aq), but this initial step in the genesis of e^-(aq) has never been simulated directly using ab initio molecular dynamics. We report the first such simulations, combining ground- and excited-state simulations of I^-(aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e^-(aq) starting from the CTTS state: I^-(aq) + hν → I^-*(aq) → I^•(aq) + e^-(aq). Notably, the equilibrated ground state of e^-(aq) evolves from I^-*(aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron-Withdrawing Groups

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the ¹ππ* (S₂) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron-withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the ¹nπ* (S₁) and ³ππ* (T₁) states on a 1.6 ± 0.4 ps time scale, which is similar to the 1.5 ps time scale earlier observed for acetylacetone [ J. Am. Chem. Soc. 2017, 139, 16576-16583, DOI: 10.1021/jacs.7b07532]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

Unraveling the Nature of Hydrogen Bonds of "Proton Sponges" Based on Car-Parrinello and Metadynamics Approaches

The nature of intra- and intermolecular non-covalent interactions was studied in four naphthalene derivatives commonly referred to as "proton sponges". Special attention was paid to an intramolecular hydrogen bond present in the protonated form of the compounds. The unsubstituted "proton sponge" served as a reference structure to study the substituent influence on the hydrogen bond (HB) properties. We selected three compounds substituted by methoxy, amino, and nitro groups. The presence of the substituents either retained the parent symmetry or rendered the compounds asymmetric. In order to reveal the non-covalent interaction properties, the Hirshfeld surface (HS) was computed for the crystal structures of the studied compounds. Next, quantum-chemical simulations were performed in vacuo and in the crystalline phase. Car-Parrinello molecular dynamics (CPMD), Path Integral Molecular Dynamics (PIMD), and metadynamics were employed to investigate the time-evolution changes of metric parameters and free energy profile in both phases. Additionally, for selected snapshots obtained from the CPMD trajectories, non-covalent interactions and electronic structure were studied. Quantum theory of atoms in molecules (QTAIM) and the Density Overlap Regions Indicator (DORI) were applied for this purpose. It was found based on Hirshfeld surfaces that, besides intramolecular hydrogen bonds, other non-covalent interactions are present and have a strong impact on the crystal structure organization. The CPMD results obtained in both phases showed frequent proton transfer phenomena. The proton was strongly delocalized in the applied time-scale and temperature, especially in the PIMD framework. The use of metadynamics allowed for tracing the free energy profiles and confirming that the hydrogen bonds present in "proton sponges" are Low-Barrier Hydrogen Bonds (LBHBs). The electronic and topological analysis quantitatively described the temperature dependence and time-evolution changes of the electronic structure. The covalency of the hydrogen bonds was estimated based on QTAIM analysis. It was found that strong hydrogen bonds show greater covalency, which is additionally determined by the proton position in the hydrogen bridge.

Classical Nuclear Motion: Comparison to Approaches with Quantum Mechanical Nuclear Motion

Ab initio molecular dynamics combines a classical description of nuclear motion with a density-functional description of the electronic cloud. This approach nicely describes chemical reactions. A possible conclusion is that a quantum mechanical description of nuclear motion is not needed. Using Occam’s razor, this means that, being the simpler approach, classical nuclear motion is preferable. In this paper, it is claimed that nuclear motion is classical, and this hypothesis will be tested in comparison to methods with quantum mechanical nuclear motion. In particular, we apply ab initio molecular dynamics to two photoreactions involving hydrogen. Hydrogen, as the lightest element, is often assumed to show quantum mechanical tunneling. We will see that the classical picture is fully sufficient. The quantum mechanical view leads to phenomena that are difficult to understand, such as the entanglement of nuclear motion. In contrast, it is easy to understand the simple classical picture which assumes that nuclear motion is steady and uniform unless a force is acting. Of course, such a hypothesis must be verified for many systems and phenomena, and this paper is one more step in this direction.

ReaxFF-based nonadiabatic dynamics method for azobenzene derivatives

ReaxFF reactive force fields have been parameterized for the ground and first excited states of azobenzene and its derivatives. In addition, an extended set of ab initio reference data ensures wide applicability, including to azosystems in complex environments. Based on the optimized force fields, nonadiabatic surface hopping simulations produce photoisomerization quantum yields and decay times of azobenzene, both in the gas phase and in n-hexane solution, in reasonable agreement with higher level theory and experiment. The transferability to other azo-compounds is illustrated for different arylazopyrazoles as well as ethylene-bridged azobenzene. Moreover, it has been shown that the model can be easily extended to adsorbates on metal surfaces. The simulation of the ring-opening of cyclobutene triggered by the photoisomerization of azobenzene in a macrocycle highlights the advantages of a reactive force field model.

Either Accurate Singlet-Triplet Gaps or Excited-State Structures: Testing and Understanding the Performance of TD-DFT for TADF Emitters

The energy gap between the lowest singlet and triplet excited states (ΔE_ST) is a key property of thermally activated delayed fluorescence (TADF) emitters, where these states are dominated by charge-transfer (CT) character. Despite its well-known shortcomings concerning CT states, time-dependent density functional theory (TD-DFT) is widely used to predict this gap and study TADF. Moreover, polar CT states exhibit a strong interaction with their molecular environment, which further complicates their computational description. Addressing these two major challenges, this work studies the performance of Tamm-Dancoff-approximated TD-DFT (TDA-DFT) on the recent STGABS27 benchmark set,¹ exploring different strategies to include orbital and structural relaxation, as well as dielectric embedding. The results show that the best-performing strategy is to calculate ΔE_ST at the ground-state structure using functionals with a surprisingly small amount of Fock exchange of ≈10% and without a (complete) solvent model. However, as this approach heavily relies on error cancellation to mimic dielectric relaxation, it is not robust and exhibits large systematic deviations in excited state energies, state characters, and structures. More rigorous approaches, including state-specific solvation, do not share these systematic deviations, but their predicted ΔE_ST values exhibit larger statistical errors. We thus conclude that for the description of CT states in dielectric environments, none of the tested TDA-DFT methods is competitive with the recently presented ROKS/PCM approach regarding robustness, accuracy, and computational efficiency.

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core-valence separated (FC-CVS) EOM-CCSD approach - a successful variant of EOM-CCSD tailored towards core excitations.

Non-iterative Method for Constructing Valence Antibonding Molecular Orbitals and a Molecule-adapted Minimum Basis

While bonding molecular orbitals exhibit constructive interference relative to atomic orbitals, antibonding orbitals show destructive interference. When full localization of occupied orbitals into bonds is possible, bonding and antibonding orbitals exist in 1:1 correspondence with each other. Antibonding orbitals play an important role in chemistry because they are frontier orbitals that determine orbital interactions, as well as much of the response of the bonding orbital to perturbations. In this work, we present an efficient method to construct antibonding orbitals by finding the orbital that yields the maximum opposite spin pair correlation amplitude in second order perturbation theory (AB2) and compare it with other techniques with increasing the size of the basis set. We conclude the AB2 antibonding orbitals are a more robust alternative to the Sano orbitals as initial guesses for valence bond calculations, due to having a useful basis set limit. The AB2 orbitals are also useful for efficiently constructing an active space, and work as good initial guesses for valence excited states. In addition, when combined with the localized occupied orbitals, and relocalized, the result is a set of molecule-adapted minimal basis functions that is built without any reference to atomic orbitals of the free atom. As examples, they are applied to population analysis of halogenated methane derivatives, H-Be-Cl, and \ce{SF6} where they show some advantages relative to good alternative methods.

Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level

A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.

Computing x-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn-Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ∼0.6 eV rms error for the K-edges of C-F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

Revisiting the Performance of Time-Dependent Density Functional Theory for Electronic Excitations: Assessment of 43 Popular and Recently Developed Functionals from Rungs One to Four

In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between computational efficiency and accuracy. The functionals ωB97X-D and BMK are found to offer the best performance overall with a root-mean square error (RMSE) of around 0.27 eV, better than the computationally more demanding CIS(D) wave function method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for time-dependent density functional theory excitation energies, though hybrid meta generalized-gradient approximations (meta-GGAs) are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.

Relativistic Orbital-Optimized Density Functional Theory for Accurate Core-Level Spectroscopy

Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital-optimized density functional theory (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding a ∼0.5 eV root-mean-square error versus experiment with the modern SCAN (and related) functionals. This marks a significant improvement over the >50 eV deviations that are typical for the popular time-dependent DFT (TDDFT) approach. Consequently, experimental spectra are quite well reproduced by OO-DFT/X2C, sans empirical shifts for alignment. OO-DFT/X2C combines high accuracy with ground state DFT cost and is thus a promising route for computing core-level spectra of third period elements. We also explored K and L edges of 3d transition metals to identify limitations of the OO-DFT/X2C approach in modeling the spectra of heavier atoms.

Electric Field Effects on Photoelectrochemical Water Splitting: Perspectives and Outlook

The grand challenges in renewable energy lie in our ability to comprehend efficient energy conversion systems, together with dealing with the problem of intermittency via scalable energy storage systems. Relatively little progress has been made on this at grid scale and two overriding challenges still need to be addressed: (i) limiting damage to the environment and (ii) the question of environmentally friendly energy conversion. The present review focuses on a novel route for producing hydrogen, the ultimate clean fuel, from the Sun, and renewable energy source. Hydrogen can be produced by light-driven photoelectrochemical (PEC) water splitting, but it is very inefficient; rather, we focus here on how electric fields can be applied to metal oxide/water systems in tailoring the interplay with their intrinsic electric fields, and in how this can alter and boost PEC activity, drawing both on experiment and non-equilibrium molecular simulation.

A theoretical perspective of the ultrafast transient absorption dynamics of CsPbBr3

Transient absorption spectra (TAS) of lead halide perovskites can provide important insights into the nature of the photoexcited state dynamics of this prototypical class of materials. Here, we perform ground and excited state molecular dynamics (MD) simulations within a restricted open shell Kohn-Sham (ROKS) approach in order to interpret the characteristic features of the TAS of CsPbBr₃ . Our results reveal that properties such as the finite temperature band gap, the Stokes shift, and therefore, also the TAS are strongly size-dependent. Our TAS simulations show an early positive red-shifted feature on the fs scale that can be explained by geometric relaxation in the excited state. As excited-state processes can crucially affect the electronic properties of this class of photoactive materials, our observations are an important ingredient for further optimization of lead halide based optoelectronic devices.

PCM-ROKS for the Description of Charge-Transfer States in Solution: Singlet-Triplet Gaps with Chemical Accuracy from Open-Shell Kohn-Sham Reaction-Field Calculations

The adiabatic energy gap between the lowest singlet and triplet excited states ΔE_ST is a central property of thermally activated delayed fluorescence (TADF) emitters. Since these states are dominated by a charge-transfer character, causing strong orbital-relaxation and environmental effects, an accurate prediction of ΔE_ST is very challenging, even with modern quantum-chemical excited-state methods. Addressing this major challenge, we present an approach that combines spin-unrestricted (UKS) and restricted open-shell Kohn-Sham (ROKS) self-consistent field calculations with a polarizable-continuum model and range-separated hybrid functionals. Tests on a new representative benchmark set of 27 TADF emitters with accurately known ΔE_ST values termed STGABS27 reveal a robust and unprecedented performance with a mean absolute deviation of only 0.025 eV (∼0.5 kcal/mol) and few deviations greater than 0.05 eV (∼1 kcal/mol), even in electronically challenging cases. Requiring only two geometry optimizations per molecule at the ROKS/UKS level in a compact double-ζ basis, the approach is computationally efficient and can routinely be applied to molecules with more than 100 atoms.

Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules

Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the M_S = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Initial Maximum Overlap Method for Large Systems by the Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Technique

Quantum chemistry offers a large variety of methods to treat excited states. Many of them are based on a multireference wave function ansatz and are therefore characterized by an intrinsic complexity and high computational costs. To overcome these drawbacks and also some limitations of simpler single-reference approaches (e.g., configuration interaction singles and time-dependent density functional theory), the single-determinant Δself-consistent field-initial maximum overlap method (ΔSCF-IMOM) has been proposed. This strategy substitutes the aufbau principle with a criterion that occupies molecular orbitals at successive SCF iterations on the basis of their maximum overlap with a proper set of guess orbitals for the target excited state. In this way, it prevents the SCF process to collapse to the ground state wave function and provides excited state single Slater determinant solutions to the SCF equations. Here, we propose to extend the applicability of the IMOM to the treatment of localized excited states of large systems. To accomplish this task, we coupled it with the QM/ELMO (quantum mechanics/extremely localized molecular orbitals) strategy, a quantum mechanical embedding method in which the most chemically relevant part of the system is treated with traditional quantum chemical approaches, while the rest is described by extremely localized molecular orbitals transferred from recently constructed libraries or proper model molecules. After presenting the theoretical foundations of the new IMOM/ELMO technique, in this paper, we will show and discuss the results of preliminary test calculations carried out on both model systems (i.e., decanoic acid, decene, decapentaene, and solvated acrolein) and a system of biological interest (flavin mononucleotide in the flavodoxin protein). We observed that, for localized excited states, the new IMOM/ELMO method provides reliable results, and it reproduces the outcomes of fully IMOM calculations within the chemical accuracy threshold (i.e., 0.043 eV) by including only a limited number of atoms in the QM region. Furthermore, the first application of our embedding technique to a larger biological system gave completely plausible results in line with those obtained through more traditional quantum mechanical methods, thus opening the possibility of using the new approach in future investigations of photobiology problems.

Substituent Effect in the Cation Radicals of Monosubstituted Benzenes

In 30 monosubstituted benzene cation radicals, studied at the ωB97XD/aug-cc-pVTZ level, the phenyl rings usually adopt a compressed form, but a differently compressed form-equivalent to an elongated one-may coexist. The computational and literature ionization potentials are well correlated. The geometrical and magnetic aromaticity, estimated using HOMA and NICS indices, show the systems to be structurally aromatic but magnetically antiaromatic or only weakly aromatic. The partial charge is split between the substituent and ring and varies the most at C(ipso). In the ring, the spin is 70%, concentrated equally at the C(ipso) and C(p) atoms. The sEDA(D) and pEDA(D) descriptors of the substituent effect in cation radicals, respectively, were determined. In cation radicals, the substituent effect on the σ-electron system is like that in the ground state. The effect on the π-electron systems is long-range, and its propagation in the radical quinone-like ring is unlike that in the neutral molecules. The pEDA(D) descriptor correlates well with the partial spin at C(ipso) and C(p) and weakly with the HOMA(D) index. The correlation of the spin at the ring π-electron system and the pEDA(D) descriptor shows that the electron charge supplied to the ring π-electron system and the spin flow oppositely.

Orbital Optimized Density Functional Theory for Electronic Excited States

Density functional theory (DFT) based modeling of electronic excited states is of importance for investigation of the photophysical/photochemical properties and spectroscopic characterization of large systems. The widely used linear response time-dependent DFT (TDDFT) approach is, however, not effective at modeling many types of excited states, including (but not limited to) charge-transfer states, doubly excited states, and core-level excitations. In this perspective, we discuss state-specific orbital optimized (OO) DFT approaches as an alterative to TDDFT for electronic excited states. We motivate the use of OO-DFT methods and discuss reasons behind their relatively restricted historical usage (vs TDDFT). We subsequently highlight modern developments that address these factors and allow efficient and reliable OO-DFT computations. Several successful applications of OO-DFT for challenging electronic excitations are also presented, indicating their practical efficacy. OO-DFT approaches are thus increasingly becoming a useful route for computing excited states of large chemical systems. We conclude by discussing the limitations and challenges still facing OO-DFT methods, as well as some potential avenues for addressing them.

Quantum-mechanical hydration plays critical role in the stability of firefly oxyluciferin isomers: State-of-the-art calculations of the excited states

Stabilizing mechanisms of three possible isomers (phenolate-keto, phenolate-enol, and phenol-enolate) of the oxyluciferin anion hydrated with quantum explicit water molecules in the first singlet excited state were investigated using first-principles Born-Oppenheimer molecular dynamics simulations for up to 1.8 ns (or 3.7 × 10⁶ MD steps), revealing that the surrounding water molecules were distributed to form clear single-layered structures for phenolate-keto and multi-layered structures for phenolate-enol and phenol-enolate isomers. The isomers employed different stabilizing mechanisms compared to the ground state. Only the phenolate-keto isomer became attracted to the water molecules in its excited state and was stabilized by increasing the number of hydrogen bonds with nearby water molecules. The most stable isomer in the excited state was the phenolate-keto, and the phenolate-enol and phenol-enolate isomers were higher in energy by ∼0.38 eV and 0.57 eV, respectively, than the phenolate-keto. This was in contrast to the case of ground state in which the phenolate-enol was the most stable isomer.

Third-Order Møller-Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error [Bertels, L. W.; J. Phys. Chem. Lett. 2019, 10, 4170 4176]. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N⁶) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

Neutral excitation density-functional theory: an efficient and variational first-principles method for simulating neutral excitations in molecules

We introduce neutral excitation density-functional theory (XDFT), a computationally light, generally applicable, first-principles technique for calculating neutral electronic excitations. The concept is to generalise constrained density functional theory to free it from any assumptions about the spatial confinement of electrons and holes, but to maintain all the advantages of a variational method. The task of calculating the lowest excited state of a given symmetry is thereby simplified to one of performing a simple, low-cost sequence of coupled DFT calculations. We demonstrate the efficacy of the method by calculating the lowest single-particle singlet and triplet excitation energies in the well-known Thiel molecular test set, with results which are in good agreement with linear-response time-dependent density functional theory (LR-TDDFT). Furthermore, we show that XDFT can successfully capture two-electron excitations, in principle, offering a flexible approach to target specific effects beyond state-of-the-art adiabatic-kernel LR-TDDFT. Overall the method makes optical gaps and electron-hole binding energies readily accessible at a computational cost and scaling comparable to that of standard density functional theory. Owing to its multiple qualities beneficial to high-throughput studies where the optical gap is of particular interest; namely broad applicability, low computational demand, and ease of implementation and automation, XDFT presents as a viable candidate for research within materials discovery and informatics frameworks.

Exciton transfer free energy from Car-Parrinello molecular dynamics

A computational approach is presented which allows the calculation of free energies profiles for exciton transfer processes within the framework of ab initio molecular dynamics (AIMD) simulations, sampling both the electronic and the nuclear degrees of freedom. To achieve this, restraining potentials are imposed on the centres of maximally localized Wannier orbitals. The resulting quantum-mechanical orbital forces are derived analytically and implemented in an AIMD program. In analogy to classical umbrella sampling techniques, these restraints are used to control an exciton transfer by incrementally moving the Wannier centres corresponding to the electron-hole pair along a suitable reaction coordinate. The new method is applied to study exciton transfer between two stacked penta(3-methylthiophene) molecules as a function of intermolecular distance. From the resulting free energy profiles, exciton transfer rates and diffusion constants are estimated, which prove to be in line with experimental results.

Substituent Effect in the First Excited Triplet State of Monosubstituted Benzenes

The structure of 30 monosubstituted benzenes in the first excited triplet T₁ state was optimized with both unrestricted (U) and restricted open shell (RO) approximations combined with the ωB97XD/aug-cc-pVTZ basis method. The substituents exhibited diverse σ- and π-electron-donating and/or -withdrawing groups. Two different positions of the substituents are observed in the studied compounds in the T₁ state: one distorted from the plane and the other coplanar with a quinoidal ring. The majority of the substituents are π-electron donating in the first group while π-electron withdrawing in the second one. Basically, U- and RO-ωB97XD approximations yield concordant results except for the B-substituents and a few of the planar groups. In the T₁ state, the studied molecules are not aromatic, yet aromaticity estimated using the HOMA (harmonic oscillator model of aromaticity) index increases from ca. -0.2 to ca. 0.4 with substituent distortion, while in the S₁ state, they are only slightly less aromatic than in the ground state (HOMA ≈0.8 vs ≈1.0, respectively). Unexpectedly, the sEDA(T₁) and pEDA(T₁) substituent effect descriptors do not correlate with analogous parameters for the ground and first excited singlet states. This is because in the T₁ state, the geometry of the ring changes dramatically and the sEDA(T₁) and pEDA(T₁) descriptors do not characterize only the functional group but the entire molecule. Thus, they cannot provide useful scales for the substituents in the T₁ states. We found that the spin density in the T₁ states is accumulated at the C_ipso and C_p atoms, and with the substituent deformation angle, it nonlinearly increases at the former while decreases at the latter. It appeared that the gap between singly unoccupied molecular orbital and singly occupied molecular orbital (SUMO-SOMO) is determined by the change of the SOMO energy because the former is essentially constant. For the nonplanar structures, SOMO correlates with the torsion angle of the substituent and the ground-state pEDA(S₀) descriptor of the π-electron-donating substituents ranging from 0.02 to 0.2 e. Finally, shapes of the SOMO-1 instead of SOMO frontier orbitals in the T₁ state somehow resemble the highest occupied molecular orbital ones of the S₀ and S₁ states. For several planar systems, the shape of the U- and RO-density functional theory-calculated SOMO-1 orbitals differs substantially.

Tuning Electronic Structure in Layered Hybrid Perovskites with Organic Spacer Substitution

Two-dimensional layered halide organic perovskites (LHOPs) are promising candidates for many optoelectronic applications due to their interesting tunable properties. They provide a unique opportunity to control energy and charge dynamics via the independent tunability of the energy levels within the perovskite and the organic spacer for various optoelectronic applications. In the perovskite layer alone, one can replace the Pb (Sn), the halide (X = I, Br, Cl), the organic component, and the number of layers between the organic spacer layers. In addition, there are many possibilities for organic spacer layers between the perovskite layers, making it difficult for experimental methods to comprehensively explore such an extensive combinatorial space. Of particular technological interest is alignment of electronic levels between the perovskite layer and the organic spacer layer, leading to desired transfer of energy or charge carriers between perovskite and organic components. For example, as band edge absorption is almost entirely attributed to the perovskite layer, one way to demonstrate energy transfer is to observe triplet emission from organic spacers. State-of-the-art computational chemistry tools can be used to predict the properties of many stoichiometries in search for LHOPs that have the most promising electronic-structure features. In this first-principles study, we survey a group of π-conjugated organic spacer candidates for use in triplet light-emitting LHOPs. Utilizing density functional theory (DFT) and time-dependent DFT, we calculate the first singlet (S₁) and triplet (T₁) excitation energy in the ground-state geometry and the first triplet excitation energy in the excited-triplet-state relaxed geometry (T₁^*). By comparing these energies to the known lowest exciton energy level of Pb_nX_3n+1 perovskite layers (X = I, Br, Cl), we can identify organic spacer and perovskite layer pairings for possible transfer of Wannier excitons from the inorganic perovskite lattice to spin-triplet Frenkel excitons located on the organic cation. We successfully identify ten organic spacer candidates for possible pairing with perovskite layers of specific halide composition to achieve triplet light emission across the visible energy range. Molecular dynamics simulations predict that finite temperatures and perovskite environment have little influence on the average excitation energies of the two common organic spacers naphthylethylammonium (NEA) and phenelethylammonium (PEA). We find significant thermal broadening up to 0.5 eV of the optical excitation energies appearing due to finite temperature effects. The findings herein provide insights into alignment of electronic levels of the conjugated organic spacer with the layer.

Inorganic antimony halide hybrids with broad yellow emissions

Lead halide perovskites exhibit unexceptionable photoelectric properties. However, these materials are unsatisfactory in terms of stability and toxicity. Herein, we report Rb₇Sb₃Cl₁₆ as a new kind of lead free perovskite variants. This material can be easily obtained through hydrothermal reactions. The composition is determined through structure refinement, elemental analysis and X-ray photoelectron spectra. Rb₇Sb₃Cl₁₆ exhibits a broad yellow emission at 560 nm, with a Stokes shift of 175 nm and a photoluminescence quantum yield (PLQY) around 26%. Rb₇Sb₃Cl₁₆ also shows good thermal and water stability due to its inorganic composition. White light-emitting diodes (LEDs) are constructed by combining Rb₇Sb₃Cl₁₆ as yellow phosphors, our previously reported Cs₂SnCl₆:2.75%Bi as blue phosphors, and commercial UV LED chips as the excitation source, producing a white light with the Commission Internationale de'Eclairage (CIE) color coordinates at (0.39, 0.38).

Deep Learning for Nonadiabatic Excited-State Dynamics

In this work we show that deep learning (DL) can be used for exploring complex and highly nonlinear multistate potential energy surfaces of polyatomic molecules and related nonadiabatic dynamics. Our DL is based on deep neural networks (DNNs), which are used as accurate representations of the CASSCF ground- and excited-state potential energy surfaces (PESs) of CH₂NH. After geometries near conical intersection are included in the training set, the DNN models accurately reproduce excited-state topological structures; photoisomerization paths; and, importantly, conical intersections. We have also demonstrated that the results from nonadiabatic dynamics run with the DNN models are very close to those from the dynamics run with the pure ab initio method. The present work should encourage further studies of using machine learning methods to explore excited-state potential energy surfaces and nonadiabatic dynamics of polyatomic molecules.

A Versatile Computational Strategy To Characterize the Free-Energy Landscape of Excited States in Oligofluorenes

The π-rich rings of conjugated polymers and molecular rotors shape their typical properties, allowing a variety of chemical and photoresponsive phenomena. Herein, we present a computational method in the framework of classical simulations to estimate the free-energy gap between ground and excited states of oligofluorenes, which is a class of molecular rotors widely used in optoelectronic devices, because of the inner torsional rotation triggered by light irradiation. We devised multiple sets of free-energy simulations in combination with free-energy perturbation theory to predict the free-energy gap between the ground state and the first excited state. The computed excitation energies show good agreement with experiments. The approach presented herein allows one to achieve comprehensive sampling of the conformational landscape and accurate estimates of the excited state free-energy landscapes at the same time.