Diptarka Hait

Prediction of Excited-State Energies and Singlet–Triplet Gaps of Charge-Transfer States Using a Restricted Open-Shell Kohn–Sham Approach

Organic molecules with charge-transfer (CT) excited states are widely used in industry and are especially attractive as candidates for fabrication of energy efficient OLEDs, as they can harvest energy from nonradiative triplets by means of thermally activated delayed fluorescence (TADF). It is therefore useful to have computational protocols for accurate estimation of their electronic spectra in order to screen candidate molecules for OLED applications. However, it is difficult to predict the photophysical properties of TADF molecules with LR-TDDFT, as semilocal LR-TDDFT is incapable of accurately modeling CT states. Herein, we study absorption energies, emission energies, zero-zero transition energies, and singlet-triplet gaps of TADF molecules using a restricted open-shell Kohn-Sham (ROKS) approach instead and discover that ROKS calculations with semilocal hybrid functionals are in good agreement with experiments-unlike TDDFT, which significantly underestimates energy gaps. We also propose a cheap computational protocol for studying excited states with large CT character that is found to give good agreement with experimental results without having to perform any excited-state geometry optimizations.

Condensed phase electron transfer beyond the Condon approximation

Condensed phase electron transfer problems are often simplified by making the Condon approximation: the approximation that the coupling connecting two charge-transfer diabatic states is a constant. Unfortunately, the Condon approximation does not predict the existence of conical intersections, which are ubiquitous in both gas-phase and condensed-phase photochemical dynamics. In this paper, we develop a formalism to treat condensed-phase dynamics beyond the Condon approximation. We show that even for an extremely simple test system, hexaaquairon(ii)/hexaaquairon(iii) self-exchange in water, the electronic coupling is expected to fluctuate rapidly and non-Condon effects must be considered to obtain quantitatively accurate ultrafast nonequilibrium dynamics. As diabatic couplings are expected to fluctuate substantially in many condensed-phase electron transfer systems, non-Condon effects may be essential to quantitatively capture accurate short-time dynamics.

How Accurate Is Density Functional Theory at Predicting Dipole Moments? An Assessment Using a New Database of 200 Benchmark Values

Dipole moments are a simple, global measure of the accuracy of the electron density of a polar molecule. Dipole moments also affect the interactions of a molecule with other molecules as well as electric fields. To directly assess the accuracy of modern density functionals for calculating dipole moments, we have developed a database of 200 benchmark dipole moments, using coupled cluster theory through triple excitations, extrapolated to the complete basis set limit. This new database is used to assess the performance of 88 popular or recently developed density functionals. The results suggest that double hybrid functionals perform the best, yielding dipole moments within about 3.6-4.5% regularized RMS error versus the reference values-which is not very different from the 4% regularized RMS error produced by coupled cluster singles and doubles. Many hybrid functionals also perform quite well, generating regularized RMS errors in the 5-6% range. Some functionals, however, exhibit large outliers, and local functionals in general perform less well than hybrids or double hybrids.

A hybrid memory kernel approach for condensed phase non-adiabatic dynamics

The spin-boson model is a simplified Hamiltonian often used to study non-adiabatic dynamics in large condensed phase systems, even though it has not been solved in a fully analytic fashion. Herein, we present an exact analytic expression for the dynamics of the spin-boson model in the infinitely slow-bath limit and generalize it to approximate dynamics for faster baths. We achieve the latter by developing a hybrid approach that combines the exact slow-bath result with the popular non-interacting blip approximation (NIBA) method to generate a memory kernel that is formally exact to second-order in the diabatic coupling but also contains higher-order contributions approximated from the second-order term alone. This kernel has the same computational complexity as the NIBA, but is found to yield dramatically superior dynamics in regimes where the NIBA breaks down-such as systems with large diabatic coupling or energy bias. This indicates that this hybrid approach could be used to cheaply incorporate higher-order effects into second-order methods and could potentially be generalized to develop alternate kernel resummation schemes.

Communication: xDH double hybrid functionals can be qualitatively incorrect for non-equilibrium geometries: Dipole moment inversion and barriers to radical-radical association using XYG3 and XYGJ-OS

Double hybrid (DH) density functionals are amongst the most accurate density functional approximations developed so far, largely due to the incorporation of correlation effects from unoccupied orbitals via second order perturbation theory (PT2). The xDH family of DH functionals calculate energy directly from orbitals optimized by a lower level approach like B3LYP, without self-consistent optimization. XYG3 and XYGJ-OS are two widely used xDH functionals that are known to be quite accurate at equilibrium geometries. Here, we show that the XYG3 and XYGJ-OS functionals can be ill behaved for stretched bonds well beyond the Coulson-Fischer point, predicting unphysical dipole moments and humps in potential energy curves for some simple systems like the hydrogen fluoride molecule. Numerical experiments and analysis show that these failures are not due to PT2. Instead, a large mismatch at stretched bond-lengths between the reference B3LYP orbitals and the optimized orbitals associated with the non-PT2 part of XYG3 leads to an unphysically large non-Hellman-Feynman contribution to first order properties like forces and electron densities.

Bimolecular Reaction Dynamics in the Phenyl–Silane System: Exploring the Prototype of a Radical Substitution Mechanism

We present a combined experimental and theoretical investigation of the bimolecular gas-phase reaction of the phenyl radical (C6H5) with silane (SiH4) under single collision conditions to investigate the chemical dynamics of forming phenylsilane (C6H5SiH3) via a bimolecular radical substitution mechanism at a tetracoordinated silicon atom. Verified by electronic structure and quasiclassical trajectory calculations, the replacement of a single carbon atom in methane by silicon lowers the barrier to substitution, thus defying conventional wisdom that tetracoordinated hydrides undergo preferentially hydrogen abstraction. This reaction mechanism provides fundamental insights into the hitherto unexplored gas-phase chemical dynamics of radical substitution reactions of mononuclear main group hydrides under single collision conditions and highlights the distinct reactivity of silicon compared to its isovalent carbon. This mechanism might be also involved in the synthesis of cyanosilane (SiH3CN) and methylsilane (CH3SiH3) probed in the circumstellar envelope of the carbon star IRC+10216.

Delocalization Errors in Density Functional Theory are Essentially Quadratic in Fractional Occupation Number

Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps, and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C12H14 polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge-based data point could therefore permit fitting/tuning of functionals with lower delocalization error.