More Like this

1. In silico prediction of annihilators for triplet–triplet annihilation upconversion via auxiliary-field quantum Monte Carlo

   https://doi.org/10.1039/d0sc03381b  

   Weber, John L. ; Churchill, Emily M. ; Jockusch, Steffen ; Arthur, Evan J. ; Pun, Andrew B. ; Zhang, Shiwei ; Friesner, Richard A. ; Campos, Luis M. ; Reichman, David R. ; Shee, James ( January 2021 , Chemical Science)

   null (Ed.)

   The energy of the lowest-lying triplet state (T1) relative to the ground and first-excited singlet states (S0, S1) plays a critical role in optical multiexcitonic processes of organic chromophores. Focusing on triplet–triplet annihilation (TTA) upconversion, the S0 to T1 energy gap, known as the triplet energy, is difficult to measure experimentally for most molecules of interest. Ab initio predictions can provide a useful alternative, however low-scaling electronic structure methods such as the Kohn–Sham and time-dependent variants of Density Functional Theory (DFT) rely heavily on the fraction of exact exchange chosen for a given functional, and tend to be unreliable when strong electronic correlation is present. Here, we use auxiliary-field quantum Monte Carlo (AFQMC), a scalable electronic structure method capable of accurately describing even strongly correlated molecules, to predict the triplet energies for a series of candidate annihilators for TTA upconversion, including 9,10 substituted anthracenes and substituted benzothiadiazole (BTD) and benzoselenodiazole (BSeD) compounds. We compare our results to predictions from a number of commonly used DFT functionals, as well as DLPNO-CCSD(T 0 ), a localized approximation to coupled cluster with singles, doubles, and perturbative triples. Together with S1 estimates from absorption/emission spectra, which are well-reproduced by TD-DFT calculations employing the range-corrected hybrid functional CAM-B3LYP, we provide predictions regarding the thermodynamic feasibility of upconversion by requiring (a) the measured T1 of the sensitizer exceeds that of the calculated T1 of the candidate annihilator, and (b) twice the T1 of the annihilator exceeds its S1 energetic value. We demonstrate a successful example of in silico discovery of a novel annihilator, phenyl-substituted BTD, and present experimental validation via low temperature phosphorescence and the presence of upconverted blue light emission when coupled to a platinum octaethylporphyrin (PtOEP) sensitizer. The BTD framework thus represents a new class of annihilators for TTA upconversion. Its chemical functionalization, guided by the computational tools utilized herein, provides a promising route towards high energy (violet to near-UV) emission. 

   more » « less
2. Semiclassics: The hidden theory behind the success of DFT

   Okun, Pavel ; Burke, Kieron ( May 2021 , ArXivorg)

   We argue that the success of DFT can be understood in terms of a semiclassical expansion around a very specific limit. This limit was identified long ago by Lieb and Simon for the total electronic energy of a system. This is a universal limit of all electronic structure: atoms, molecules, and solids. For the total energy, Thomas-Fermi theory becomes relatively exact in the limit. The limit can also be studied for much simpler model systems, including non-interacting fermions in a one-dimensional well, where the WKB approximation applies for individual eigenvalues and eigenfunctions. Summation techniques lead to energies and densities that are functionals of the potential. We consider several examples in one dimension (fermions in a box, in a harmonic well, in a linear half-well, and in the Pöschl-Teller well. The effects of higher dimension are also illustrated with the three-dimensional harmonic well and the Bohr atom, non-interacting fermions in a Coulomb well. Modern density functional calculations use the Kohn-Sham scheme almost exclusively. The same semiclassical limit can be studied for the Kohn-Sham kinetic energy, for the exchange energy, and for the correlation energy. For all three, the local density approximation appears to become relatively exact in this limit. Recent work, both analytic and numerical, explores how this limit is approached, in an effort to deduce the leading corrections to the local approximation. A simple scheme, using the Euler-Maclaurin summation formula, is the result of many different attempts at this problem. In very simple cases, the correction formulas are much more accurate than standard density functionals. Several functionals are already in widespread use in both chemistry and materials that incorporate these limits, and prospects for the future are discussed. 

   more » « less
3. Static polarizabilities within the generalized Kohn–Sham semicanonical projected random phase approximation (GKS-spRPA)

   https://doi.org/10.1063/5.0103664  

   Balasubramani, Sree Ganesh ; Voora, Vamsee K. ; Furche, Filipp ( October 2022 , The Journal of Chemical Physics)

   An analytical implementation of static dipole polarizabilities within the generalized Kohn–Sham semicanonical projected random phase approximation (GKS-spRPA) method for spin-restricted closed-shell and spin-unrestricted open-shell references is presented. General second-order analytical derivatives of the GKS-spRPA energy functional are derived using a Lagrangian approach. By resolution-of-the-identity and complex frequency integration methods, an asymptotic [Formula: see text] scaling of operation count and [Formula: see text] scaling of storage is realized, i.e., the computational requirements are comparable to those for GKS-spRPA ground state energies. GKS-spRPA polarizabilities are assessed for small molecules, conjugated long-chain hydrocarbons, metallocenes, and metal clusters, by comparison against Hartree–Fock (HF), semilocal density functional approximations (DFAs), second-order Møller–Plesset perturbation theory, range-separated hybrids, and experimental data. For conjugated polydiacetylene and polybutatriene oligomers, GKS-spRPA effectively addresses the “overpolarization” problem of semilocal DFAs and the somewhat erratic behavior of post-PBE RPA polarizabilities without empirical adjustments. The ensemble averaged GKS-spRPA polarizabilities of sodium clusters (Na n for n = 2, 3, …, 10) exhibit a mean absolute deviation comparable to PBE with significantly fewer outliers than HF. In conclusion, analytical second-order derivatives of GKS-spRPA energies provide a computationally viable and consistent approach to molecular polarizabilities, including systems prohibitive for other methods due to their size and/or electronic structure. 

   more » « less
4. Assessment of electronic structure methods for the determination of the ground spin states of Fe( ii ), Fe( iii ) and Fe( iv ) complexes

   https://doi.org/10.1039/C7CP01263B  

   Verma, Pragya ; Varga, Zoltan ; Klein, Johannes E. ; Cramer, Christopher J. ; Que, Lawrence ; Truhlar, Donald G. ( January 2017 , Physical Chemistry Chemical Physics)

   Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe 2+ ion, gaseous FeO and 14 iron complexes using Kohn–Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe( ii ), five are Fe( iii ) and two are Fe( iv ) complexes. These are calculated using 20 exchange–correlation functionals. In particular, we use a local spin density approximation (LSDA) – GVWN5, four generalized gradient approximations (GGAs) – BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) – GAM and N12, two meta-GGAs – M06-L and M11-L, a meta-NGA – MN15-L, five hybrid GGAs – B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs – M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA – MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe 2+ ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. The local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state. 

   more » « less
5. Accelerate stochastic calculation of random-phase approximation correlation energy difference with an atom-based correlated sampling

   https://doi.org/10.1088/2516-1075/abde94  

   Chi, Y-C. ; Huang, C. ( March 2021 , Electronic structure)

   de Jong, Bert ; Nieminen, Risto (Ed.)

   A kernel polynomial method is developed to calculate the random phase approximation (RPA) correlation energy. In the method, the RPA correlation energy is formulated in terms of the matrix that is the product of the Coulomb potential and the density linear response functions. The integration over the matrix's eigenvalues is calculated by expanding the density of states of the matrix in terms of the Chebyshev polynomials. The coefficients in the expansion are obtained through stochastic sampling. Since it is often the energy difference between two systems that is of much interest in practice, another focus of this work is to develop a correlated sampling scheme to accelerate the convergence of the stochastic calculations of the RPA correlation energy difference between two similar systems. The scheme is termed the atom-based correlated sampling (ACS). The performance of ACS is examined by calculating the isomerization energy of acetone to 2-propenol and the energy of the water–gas shift reaction. Using ACS, the convergences of these two examples are accelerated by 3.6 and 4.5 times, respectively. The methods developed in this work are expected to be useful for calculating RPA-level reaction energies for the reactions that take place in local regions, such as calculating the adsorption energies of molecules on transition metal surfaces for modeling surface catalysis. 

   more » « less