Abstract
Molecular quasiparticle and excitation energies determine essentially the spectral characteristics measured in various spectroscopic experiments. Accurate prediction of these energies has been rather challenging for ground-state density functional methods, because the commonly adopted density function approximations suffer from delocalization error. In this work, by presuming a quantitative correspondence between the quasiparticle energies and the generalized Kohn–Sham orbital energies, and employing a previously developed global scaling correction approach, we achieve substantially improved prediction of molecular quasiparticle and excitation energies. In addition, we also extend our previous study on temporary anions in resonant states, which are associated with negative molecular electron affinities. The proposed approach does not require any explicit self-consistent field calculation on the excited-state species, and is thus highly efficient and convenient for practical purposes.
1. Introduction
Density function theory (DFT) (Hohenberg and Kohn, ) has made great success in practical calculations for ground-state electronic properties because of its outstanding balance between accuracy and computational cost. In the Kohn–Sham (KS) scheme of DFT (Hohenberg and Kohn, ; Kohn and Sham, ), the effective single-particle equations can be written as (by omitting the spin indices and adopting the atomic units)
Here, vext(r) is the external potential, vH(r) is the Hartree potential, vxc(r) is the local exchange-correlation (XC) potential, and {ϕm(r)} and {εm} are the KS/generalized KS (GKS) orbitals and their eigenvalues, respectively. In the GKS scheme, vxc(r) is replaced by a non-local potential, . The KS equations can be solved self-consistently to produce the ground-state energy and charge density. However, it is challenging to apply conventional ground-state density functional methods to calculate excited-state-related properties, such as the quasiparticle (QP) energies and the electronic excitation energies, which will be introduced as follows.
When an electronic system is perturbed by incoming photons or electrons, in order to preserve a single-particle picture, the concept of QP is often adopted. In a direct photoemission experiment, an electron on a molecule absorbs the energy of a photon and gets excited away from the molecule. Such a process leaves a quasihole in the molecule whose energy level is renormalized by the presence of the other electrons. Similarly, in an inverse photoemission experiment, an electron attaches to a molecule by emitting a photon, which leads to the formation of a quasielectron whose energy level is influenced by the existing electrons in the molecule (Onida et al., ).
The actual QP energies and wavefunctions can be obtained by solving the QP equations as follows (Hedin, ; Aulbur et al., ),
Here, {ψm(r)} and {ωm} are the QP wavefunctions and energies, respectively, and Σ is a non-local and energy-dependent self-energy operator, with the imaginary part of its eigenvalues giving the lifetime of the QPs. To enable practical calculations, an approximate scheme for Σ is to be employed. The most widely used scheme is the GW approximation (Hedin, ; Hybertsen and Louie, ; Aulbur et al., ; Dvorak et al., ). Therefore, regarding the calculation of QP energies, many-body perturbation theory within the GW approximation has become a popular method at present (Hedin, ; Hybertsen and Louie, ; Louie and Hybertsen, ; Aulbur et al., ; Onida et al., ; Dvorak et al., ). However, the somewhat large computational cost makes it difficult to apply the GW method to complex systems. Thus, a highly efficient and accurate method for the prediction of QP energies is sought for.
It is tempting to relate the KS/GKS orbital energies to QP energies, because the KS and GKS schemes are in conformity with an effective single-electron description. However, with conventional density functional approximations (DFAs), such as the local density approximation (LDA) (Slater, ; Vosko et al., ), generalized gradient approximations (GGAs), and hybrid functionals, the calculated KS/GKS orbital energies usually deviate severely from the QP energies. Such deviations have also led to significant underestimation of band gaps, which is largely due to the delocalization error associated with the DFAs (Cohen et al., ). In the exact DFT, the ground-state energy of a system with a fractional number of electrons, E0(N+n) (subscript 0 denotes the ground state corresponding to the fixed vext), should satisfy the Perdew–Parr–Levy–Balduz (PPLB) condition (Perdew et al., , ; Yang et al., ): E0(N + n) = (1 − n)E0(N) + nE0(N + 1), where 0 < n < 1 is a fractional number.
The PPLB condition infers that and , where I ≡ E0(N − 1) − E0(N) and A ≡ E0(N) − E0(N + 1) are the vertical ionization potential (IP) and electron affinity (EA) of the N-electron system, respectively. It has been proved (Cohen et al., ; Yang et al., ) that if the XC energy is an explicit and differentiable functional of the electron density or the KS reduced density matrix, we have and , where εHOMO and εLUMO are the energies of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), respectively. Therefore, if the PPLB condition can be satisfied, we should have I = −εHOMO and A = −εLUMO.
Within the framework of ground-state DFT, a rigorous mapping between the other remaining KS/GKS orbital energies apart from the HOMO and LUMO and the QP energies has not been established. However, in practice the Koopmans-like relations have been proposed and adopted by many authors (Hill et al., ; Coropceanu et al., ; Vargas et al., ; Bartlett, ; Gritsenko and Baerends, ; Tsuneda et al., ; Dauth et al., ; Körzdörfer et al., ; Baerends et al., ; Bartlett and Ranasinghe, ; Puschnig et al., ; Ranasinghe et al., ; Thierbach et al., ). These relations have the form of and . Here, the index i (a) denotes the occupied (virtual) KS/GKS orbital of the N-electron system from (to) which an electron is deprived (added), with () being the corresponding vertical IP (EA). It is easily recognized that these vertical IPs and EAs coincide with the energies of quasiholes and quasielectrons, i.e., and , respectively. Computationally, approximating QP energies by KS/GKS orbital energies has the advantage of requiring only a single self-consistent field (SCF) calculation for the ground state of the interested molecule.
The excited-state properties of molecular systems can be probed by photon absorption experiments (Onida et al., ). However, theoretical characterization of the excited states is rather challenging because the excited electron and the resulting hole cannot be treated separately. Numerous methods have been developed for the calculation of excitation energies. Coupled cluster (CC) (Schreiber et al., ; Silva-Junior et al., ; Winter et al., ; Wang et al., ; Dreuw and Wormit, ; Jacquemin et al., ) and multi-reference methods (Andersson et al., ; Potts et al., ; Slavicek and Martinez, ; Hoyer et al., ) are able to describe electronic excited states with a high accuracy. However, the expensive computational cost makes the application of these methods to large systems rather difficult. As a straightforward extension of the GW approach (Hedin, ; Hybertsen and Louie, ; Onida et al., ), the Bethe–Salpeter equation (BSE) (Rohlfing and Louie, ; Onida et al., ; Jacquemin et al., ) provides another method for the calculation of excited states, which is however also quite expensive. The time-dependent DFT (TDDFT) (Runge and Gross, ; Casida, ) is in principle an exact extension of the ground-state DFT, and it has been widely employed to study neutral excitations. Despite its success, TDDFT faces several challenges, such as double excitation character, multi-reference problems, and high-spin excited states (Ipatov et al., ; Laurent and Jacquemin, ; Santoro and Jacquemin, ).
Recently, a simple method has been proposed, which attempts to acquire excitation energies by using only the KS/GKS orbital energies of the molecular cations calculated by ground-state DFT (Haiduke and Bartlett, ; Mei et al., ). Such a method is referred to as the QE-DFT (QP energies from DFT), which has been employed to describe excited-state potential energy surfaces and conical intersections (Mei and Yang, ). Following the idea of QE-DFT, molecular excitation energies have also been expressed by KS/GKS orbital energies obtained with long-range corrected functionals (Hirao et al., ). Details about QE-DFT are to be presented in section 2.1.
In addition to neutral molecules, in this work we also consider the resonance states of temporary anions. A temporary anion has an energy higher than that of the neutral species, and thus its EA has a negative value. This means the anion is unstable and lasts only a short time. Although temporary anions cannot be studied by traditional spectroscopic techniques, they can be observed via resonances (sharp variations) in the cross-sections of various electron scattering processes (Jordan and Burrow, ). In the gas phase, the resonances can be identified by the electron transmission spectroscopy (Sanche and Schulz, ; Schulz, ; Jordan and Burrow, ). Since these resonances belong to the continuous part of the spectrum, they cannot be addressed by conventional electronic structure methods for bound states. A number of theoretical methods have been proposed to tackle the problem of temporary anions. For instance, it has been proposed that the attractive components of electron-molecule interaction are combined with a long-range repulsive potential to produce a barrier, behind which the excess electron can be temporarily trapped (Jordan et al., ). Moreover, the negative EAs have been studied by the GW method (Hedin, ; Hybertsen and Louie, ; Govoni and Galli, ), the electron-propagator methods (Longo et al., ; Ortiz, ; Dolgounitcheva et al., ), and the equation-of-motion coupled cluster (EOM-CC) approach (Stanton and Bartlett, ; Nooijen and Bartlett, ; Dutta et al., ; Jagau et al., ; Skomorowski et al., ; Ma et al., ), again with considerable computational cost.
In order to describe the unbound resonance states within the DFT approach, Tozer and De Proft () have proposed an approximate approach to evaluate the EA based on the KS frontier orbital energies (Kohn and Sham, ) and the accurate IP. Zhang et al. () have used directly the negative of GKS eigenvalue of the neutral ground-state molecule as an approximation of EA corresponding to the resonance state of the anion. The good accuracy was made possible because of the use of the global scaling correction (GSC) (Zheng et al., ), which will be introduced later. At the same time, another method has been developed to evaluate the negative EA from the GKS eigenvalue of the neutral ground states (Carmona-Espíndola et al., ). Different from GSC, this method is designed to impose the derivative discontinuity of the exact XC potential. Our work (Zhang et al., ) proceeded the work of Mei et al. () and that of Haiduke and Bartlett (Haiduke and Bartlett, ) in the direct use of GKS eigenvalues of the N-electron ground state to approximate the excited state energy of the corresponding (N + 1)-electron system, with the excited state of the (N + 1)-electron system being a unbound resonance state.
For achieving an accurate prediction of QP energies with ground-state density functional methods, it is crucial to reduce the delocalization error associated with the adopted DFA. Enormous efforts have been made, which have led to the development of the GSC (Zheng et al., ) and local scaling correction (LSC) (Li et al., ) approaches, which alleviate substantially the delocalization error of various DFAs for systems involving global and local fractional electron distributions, respectively. This is done by imposing explicitly the PPLB condition on the form of DFA. Recently, a localized orbital scaling correction (LOSC) (Li et al., ; Su et al., ) has been constructed to join the merits of GSC and LSC. The LOSC approach is capable of correcting the energy, energy derivative, and electron density of any finite system in a self-consistent and size-consistent manner. In particular, the LOSC approach has been applied in conjunction with the QE-DFT to predict QP and excitation energies of molecules (Mei et al., ).
In this work, we revisit the non-empirical GSC approach (Zheng et al., , , ; Zhang et al., ) and explore the possibility of using it to achieve an accurate prediction of QP and excitation energies. With a perturbative treatment of the orbital relaxation induced by the addition (removal) of an infinitesimal amount of electron to (from) a molecule, the GSC approach has been demonstrated to improve systematically the prediction of KS frontier orbital energies and band gaps of molecules (Zhang et al., ). Based on the idea of QE-DFT, we will extend the scope of GSC from the frontier orbitals to the other KS/GKS orbitals.
The remainder of this paper is organized as follows. In section 2, we present the QE-DFT method to calculate the QP and excitation energies within the framework of ground-state DFT, as well as the GSC approach to achieve the accurate KS/GKS orbital energies. In section 3, numerical results of the QP energies, electronic excitation energies, and resonance energies are presented and discussed. Finally, we summarize this work in section 4.
2. Methodology
2.1. QE-DFT Method for the Calculation of QP, Excitation, and Resonance Energies
In the QE-DFT method, the following Koopmans-like relations are adopted, which use the energies of occupied and virtual KS/GKS orbitals to approximate the quasihole and quasielectron energies, respectively.
Here, {εi(N)} and {εa(N)} are the occupied and virtual orbital energies of the N-electron system, respectively. Ea(N + 1) denotes the energy of the (N + 1)-electron system formed by adding an excess electron to the ath virtual orbital of the N-electron system at its ground state. Note that the subscript a refers to the N-electron system, and the value of such an orbital index may vary in the (N + 1)-electron system because of the relaxation and re-ordering of the orbitals upon the perturbation induced by electron addition. A similar argument applies to the ith occupied orbital of the N-electron system.
From Equation (3), it is obvious that the excitation energy of an electron from the HOMO to a virtual orbital of the N-electron system, which corresponds to the ath orbital of the (N − 1)-electron system, can be calculated as follows (Haiduke and Bartlett, ; Mei et al., ):
Here, in the second equality of Equation (4), we have chosen to use the (N − 1)-electron system as a reference system. This means the electronic excitation from the HOMO to a virtual orbital can be regarded as consisting of two steps: first an electron is removed from the HOMO of the N-electron system, giving rise to an (N − 1)-electron system in its ground state, and then an excess electron is put to the ath virtual orbital of the (N − 1)-electron system, which is energetically higher than the frontier orbitals, resulting in an excited N-electron system. Accordingly, E0(N − 1) is the ground-state energy of the (N − 1)-electron system, and Ea(N) denotes the energy of the N-electron system that is finally obtained. Thus, Equation (4) can describe the excitation of an electron from the HOMO to any virtual KS/GKS orbital (LUMO and above), as long as the orbital finds its counterpart in the (N − 1)-electron system.
Specifically, if we presume the (N − 1)-electron reference system contains one more spin-α electrons than spin-β electrons, the first triplet-state excitation energy of the N-electron system is calculated by
Higher triplet-state excitation energies can be calculated similarly.
We now consider the first singlet excited state formed by adding a spin-β electron to the ath virtual orbital of the ground-state (N − 1)-electron system. It is well-known that the accurate calculation of open-shell singlet states is quite challenging for the density functional methods within the KS/GKS scheme. This is because the electronic wavefunction naturally involves more than one Slater determinant, and such a multireference character is hardly captured by the presently used DFAs due to their intrinsic static correlation error (Cohen et al., ). Moreover, an open-shell singlet wavefunction in the form of a single Slater determinant of KS/GKS orbitals is not an eigenstate of the total spin operator. In practice, people have attempted to circumvent the problem of static correlation error by explicitly using more than one Slater determinant. For instance, the singlet-state energy of an N-electron system has been written as (Ess et al., )
Here, EM(N) represents the energy of a single-Slater-determinant wavefunction with the excited spin-β electron occupying the virtual orbital. The second term on the right-hand side is a correction to the singlet-state energy, which accounts for the spin contamination of the single-Slater-determinant wavefunction, with χ being a parameter. The singlet-state excitation energy of the N-electron is thus obtained as
To improve the accuracy of ΔES, a spin purification procedure has been proposed (Ziegler et al., ), which amounts to χ = 1 in Equation (6). Specifically, the first singlet-state excitation energy is calculated by
Obviously, with the QE-DFT method, the calculation of excitation energies requires the SCF calculations to be performed explicitly only for the ground-state (N − 1)-electron system.
Regarding temporary anions, we only consider the scenario that the LUMO of the neutral molecule is already an unbound orbital, which corresponds to a negative EA. Consequently, addition of an excess electron to the LUMO gives rise to a resonant state. Traditionally, the molecular EA is obtained by performing SCF calculations separately for the neutral and anionic systems and taking the energy difference between them. This is referred to the ΔSCF method. However, in practice it is extremely difficult to carry out an SCF calculation for the anionic species if it is in a resonant state.
Since the LUMO is a frontier orbital, the PPLB condition holds exactly, and thus the negative EA can be obtained directly from the positive LUMO energy via the following equality:
By using Equation (9), the SCF calculation on the temporary anion that is potentially problematic is no longer needed.
2.2. GSC Approach for the Accurate Prediction of KS/GKS Orbital Energies
From section 2.1, the prediction of QP and excitation energies is transformed to the accurate calculation of KS/GKS orbital energies. To this end, we employ a non-empirical GSC approach developed in our previous works Zheng et al. (); Zhang et al. () to reduce the delocalization error of some frequently used DFAs. It has been demonstrated that the GSC approach greatly improves the accuracy of the frontier KS/GKS orbitals. In the following, we shall go beyond the frontier orbitals and extend the application of GSC to all the KS/GKS orbitals.
In the KS or GKS scheme, the total electronic energy in the ground state is E0(N) = Ts + Vext + J + Exc. With the KS/GKS orbitals fixed as the electron number is varied, the KS kinetic energy Ts and external potential energy Vext are linear in ρ(r), while the electron Coulomb energy J is quadratic, and Exc is usually non-linear in ρ(r). The GSC approach establishes a linear energy function that satisfies the PPLB condition,
by linearizing both J and Exc with respect to the fractional electron number n. The difference between Ẽ0 and E0 is just the GSC for the energy:
Here, can express explicitly by the electron density and some other quantities, where ϕm(r) and nm are the mth KS/GKS orbital and electron occupation number, respectively. For simplicity, the spin indices are omitted.
The addition of the n fractional electron to the LUMO presents a perturbation to the N-electron system, and the change in electron density in response to such a perturbation is
where and are the first- and second-order Fukui functions (Parr and Yang, ; Yang et al., ; Yang and Parr, ), respectively. Accordingly, the relaxation of KS/GKS orbitals upon the addition of n fractional electron can be expanded in a perturbative series as , with being the kth-order orbital relaxation. Thus, the Fukui functions can be expressed explicitly in terms of orbital relaxation. For instance, the first-order Fukui function is
Here, the subscript f denotes the frontier orbital, with f = LUMO (f = HOMO) in the case of electron addition (removal). The explicit forms of orbital relaxation up to the third order have been derived and provided in Zhang et al. (), with all the perturbation Hamiltonian matrices determined by a self-consistent process. Ultimately, all orders of orbital relaxation and Fukui quantities are expressed in terms of {ϕm(r)} and {εm} of the N-electron system. The scaling correction to the frontier orbital energy is then evaluated by the Janak's theorem (Janak, ) in a post-SCF manner,
where is the kth-order correction to the frontier orbital energy.
An accurate prediction of molecular IP and EA has been achieved by employing the GSC approach (Zheng et al., , ; Zhang et al., , ) via
In practical calculations, the perturbative series needs to be truncated at a certain order. It is worth pointing out that the accuracy of the GSC does not necessarily increase with further inclusion of higher order orbital relaxation. This is because the present form of GSC only treats the exchange energy Ex, while the correlation energy Ec is presumed to be much smaller and hence its correction is omitted. However, the correction to Ec may have a comparable magnitude to the high-order corrections to Ex. For instance, regarding the prediction of EA, while the inclusion of first-order orbital relaxation is found optimal for the LDA and GGA (such as BLYP, Becke, ; Lee et al., ), the inclusion of orbital relaxation up to second-order is most favorable for the hybrid functional B3LYP (Lee et al., ; Becke, ).
To extend the GSC approach beyond the frontier KS/GKS orbitals, we presume that the PPLB condition could be generalized to the following energy linearity relation:
The GSC to the energy of the (N + n)-electron system is
where Ea(N + n) is the energy of the (N + n)-electron system in an excited state, since the n fractional electron is now added to the ath virtual orbital of the N-electron system. Similarly, the changes of electron density and KS/GKS orbitals in response to the perturbation caused by the electron addition process, as well as their contributions to , are calculated by using the self-consistent perturbation theory presented in Zhang et al. (). This finally gives rise to the GSC to the KS/GKS orbital energies:
Likewise, for the scenario that n fractional electron is deprived from the ith occupied orbital of the N-electron system, we have
With the QE-DFT method, we can now use the scaling corrected KS/GKS orbital energies to approximate the QP energies and the related vertical IPs and EAs as follows:
3. Results and Discussions
3.1. QP Energies
3.1.1. Quasihole Energies of Molecules
Because of the lack of highly accurate experimental or theoretical data for the molecular quasielectron energies (except for those associated with the LUMOs), in this work we only compare the calculated quasihole energies that are associated with the occupied KS/GKS orbitals to the reference data available in the literature.
We first look into 56 quasihole energies of 12 molecules by calculating the scaling corrected orbital energies, and make comparison with experimentally measured vertical IPs. The examined molecules are cyanogen, CO, acetylene, water, ethylene, ammonia, acetonitrile, fluoromethane, benzene, naphthalene, furan, O2, and formic acid, which exhibit diversified geometric and electronic features. Among these molecules, the geometries of benzene and naphthalene are extracted from Mei et al. (), while the structures of the other molecules are optimized with the B3LYP/6-311+g** method by using the Gaussian09 package (Frisch et al., ).
The GSC approach presented in section 2.2 is employed to correct the orbital energies calculated by various mainstream DFAs, including the LDA, the GGAs (BLYP and PBE, Perdew et al., ), and the hybrid functional B3LYP. For these DFAs, the orbital relaxation up to the second order is considered for calculating the scaling corrections of the occupied orbital energies. The GSC approach is implemented in an in-house built quantum chemistry software package QM4D (Hu et al., ).
Figure 1 compares the averaged deviations of the calculated and from the quasihole energies {ωi} extracted from the experimentally measured vertical IPs. It is shown clearly that the mean absolute errors (MAEs) associated with the original DFAs are as large as several eVs, while by applying the GSC approach, the MAEs are substantially reduced to less than 0.5 eV. Take the B3LYP functional as an example. It yields an MAE of 3.05 eV, which is the smallest among all the uncorrected DFAs, and the MAE is greatly reduced to 0.28 eV by using the GSC-B3LYP. If instead the orbital relaxation is treated up to the first and third order, the MAE becomes 0.74 eV and 0.43 eV, respectively. The dependence on the order of orbital relaxation is consistent with the trend observed in our previous work (Zhang et al., ).
Figure 1
In a previous study by Chong et al. (), 10 out of 12 molecules examined in Figure 1 (without benzene and naphthalene) have been investigated by calculating their KS/GKS orbital energies by using an approximate XC potential obtained with the statistical averaging of (model) orbital potentials (SAOP). For these 10 molecules, the MAE reported in Chong et al. () is 0.38 eV, while the GSC-B3LYP yields a somewhat smaller MAE of 0.28 eV, albeit the different molecular geometries and basis sets adopted.
The comparison between the individual orbital energies {−εi} calculated by B3LYP and GSC-B3LYP and the experimentally measured vertical IPs is depicted in Figure 2. It is apparent that the uncorrected orbital energies deviate systematically and significantly from the experimental QP energies, while such deviations are largely alleviated by applying the GSC approach.
Figure 2
3.1.2. Photoemission Spectra
The QP energies can also be extracted from the peak positions of experimentally measured photoemission spectra (PES). We employ the QE-DFT to study the PES of 14 molecules. The same molecular geometries and basis set (cc-pVTZ, Dunning, ; Woon and Dunning Jr, ) as those adopted in Mei et al. () are used here. The PES are simulated by setting the energy of each KS/GKS orbital as the center of a QP peak, and assuming all QP peaks have the same amplitude and are broadened by the same Gaussian function with λ = 0.2eV.
Figure 3 depicts the experimentally measured and theoretically simulated PES of a maleic anhydride and a benzoquinone, while those of the other 12 molecules are presented in Supplementary Material. Clearly, both the PBE and B3LYP yield considerable errors in the peak positions of the simulated PES. More specifically, they tend to predict much too high quasihole energies and too low quasielectron energies. This is because the uncorrected DFAs (PBE and B3LYP) suffer from delocalization error, as they violate the rigorous PPLB condition and the extended energy linearity relation.
Figure 3
The use of GSC improves significantly the simulated PES. For GSC-PBE, the orbital relaxation is considered up to the first and second order for the virtual and occupied KS/GKS orbitals, respectively; while for GSC-B3LYP, the orbital relaxation is included up to second order for all the KS/GKS orbitals. From the comparison shown in Figure 3, it is evident that the GSC-DFAs achieve at least the same level of accuracy as the results of GW method (Knight et al., ). Moreover, the computational cost of the QE-DFT method by using a GSC-DFA is supposedly much cheaper than that of the GW method, because the former requires only a single SCF calculation at the DFT level.
3.2. Energies of Low-Lying Excited States
We now turn to the energies of low-lying excited states of molecules. By employing the QE-DFT method, we carry out calculations on 48 low excitation energies of the 16 molecules investigated previously in Mei et al. (). The cationic species of all these molecules indeed contain one more spin-α electrons than spin-β electrons, and hence their triplet and singlet excitation energies are computed by using equations 5 and 7, respectively. Since the calculations involve only the virtual KS/GKS orbitals of the cations, the orbital relaxation is considered up to the second order for GSC-B3LYP and up to the first order for other GSC-DFAs, respectively.
Figure 4 compares the MAEs of different types of excitation energies calculated by various DFAs, and the detailed results can be found in the Supplementary Material.
Figure 4
Intriguingly, for the lowest (HOMO-to-LUMO) triplet excitations, the uncorrected DFAs yield reasonably accurate excitation energies, and the application of the GSC approach does not lead to any improvement. In particular, B3LYP yields an MAE as small as 0.17 eV for the T1 excitations. Such an appealing accuracy is likely due to the cancellation of delocalization error. Equation (5) involves the difference between a pair of virtual orbital energies. Thus, when two virtual orbitals are close in energy, their associated delocalization errors are expected to cancel out (Mei et al., ). Consequently, the GSC approach does not help. Such an error cancellation mechanism becomes less favorable for higher excitations. For instance, as displayed in Figure 4, the uncorrected DFAs tend to yield a larger MAE for the T2 excitations, and applying the GSC indeed leads to improved accuracy. The latter is because the scaling correction to each individual QP energy starts to take effect.
For the S1 excitations, the GSC-DFAs yield MAEs that are somewhat larger than the original DFAs. This is because a second type of systematic error of the DFAs, the static correlation error, becomes prominent and significantly affects the calculated ΔES1. Taking the LDA functional as an example. If the first singlet excited state is described by a single Slater determinant, i.e., by setting χ = 0 in Equation (6) and (7), the MAEs of ΔES1 are 1.44eV and 1.33eV for the original LDA and GSC-LDA, respectively. In contrast, after adopting the spin purification procedure (by setting χ = 1), the MAEs reduce to 0.73eV and 0.92eV for the original LDA and GSC-LDA, respectively. Therefore, the spin purification procedure indeed diminishes the MAE of ΔES1 by circumventing the problem of static correlation error. However, the somewhat larger MAE of the GSC-LDA seems to indicate that the spin purification formula is not entirely compatible with the present GSC scheme. For higher singlet excitations, the static correlation error becomes much less significant, as signified by the much smaller second term on the right-hand side of Equation (6). Consequently, the MAE of ΔES2 experiences a rather minor change by invoking the spin purification. For instance, the MAE increases slightly from 0.68eV to 0.70eV with the original LDA, while it reduces slightly from 0.57eV to 0.53eV with the GSC-LDA.
Among all the DFAs examined, the GSC-B3LYP functional achieves an optimal performance for all the low-lying excitations studied. The overall accuracy of GSC-B3LYP is comparable to the TDDFT-B3LYP. This affirms that it is entirely possible and practical to access excited-state properties of molecules within the framework of ground-state DFT.
We further extend our test to cover three transition metal atoms M (M = Cu, Ag, Au) and nine transition metal compounds MX (X = F, Cl, Br). The calculated results are presented in Table 1 along with the CCSD(T) results and the available experimental data. The small MAEs between the calculated results and experimental data further affirm the applicability of the GSC approach.
Table 1
| Final state | Exp. | CCSD(T) | B3LYP | GSC-B3LYP | ||
|---|---|---|---|---|---|---|
| IP | Cu | d10 | 7.73 | 7.69 | 5.35 | 7.90 |
| Ag | d10 | 7.58 | 7.52 | 5.30 | 7.80 | |
| Au | d10 | 9.23 | 9.13 | 6.61 | 9.28 | |
| MAEa | 0.07 | 2.43 | 0.14 | |||
| T1 | Cu | d10p1 | 3.81 | 3.87 | 4.92 | 3.94 |
| Ag | d10p1 | 3.84 | 3.74 | 4.80 | 3.87 | |
| Au | d10p1 | 4.95 | 5.01 | 6.06 | 5.00 | |
| CuF | 3Σ+ | 1.81 | 1.81 | 2.73 | 2.05 | |
| CuCl | 3Σ+ | 2.35 | 2.43 | 3.24 | 2.35 | |
| CuBr | 3Σ+ | 2.54 | 2.56 | 3.19 | 2.35 | |
| AgF | 13Σ+ | 3.09 | 3.27 | 3.26 | 2.54 | |
| AgCl | 13Σ+ | N/A | 3.50 | 3.66 | 2.97 | |
| AgBr | 13Σ+ | N/A | 3.37 | 3.43 | 2.83 | |
| AuF | 13Σ+ | N/A | 2.08 | 2.17 | 1.80 | |
| AuCl | 13Σ+ | N/A | 2.53 | 2.61 | 2.14 | |
| AuBr | 13Σ+ | N/A | 2.60 | 2.52 | 2.09 | |
| MAEa | 0.07 | 0.83 | 0.17 | |||
| MAEb | 0.50 | 0.31 |
Calculated and experimental ionization potentials (IPs) and T1 excitation energies of various transition metal atoms and compounds.
The experimental excitation energies and the CCSD(T) calculation results are extracted from Guichemerre et al. (). The experimental geometries (Guichemerre et al., ) and the def2-TZVPD basis set (Schäfer et al., ; Rappoport and Furche, ) are used for all the calculations.
The MAE is between the calculated values and the experimental data.
The MAE is between the DFT and the CCSD(T) results.
3.3. Resonance Energies of Temporary Anions
For a temporary anion in a resonant state, the corresponding neutral molecule has a negative EA, for which the conventional ΔSCF method often yields problematic results. This is because the choice of an appropriate basis set is difficult for the SCF calculation of a temporary anion. On the one hand, the energy of a temporary anion is rather sensitive to the inclusion of diffuse basis functions (Guerra, ). On the other hand, the diffuse basis functions may artificially delocalize the excess electron (Cohen et al., , ), and thus result in incorrect electron density distribution.
Alternatively, using the scaling corrected LUMO energy to determine the energy of the temporary anion has made impressive progress. It has been demonstrated that the GSC-PBE functional predicts highly accurate negative EAs by using Equation (16) (Zhang et al., ). For a set of 38 molecules proposed in Tozer and De Proft (), the resulting MAE is as small as 0.18 eV with the aug-cc-pVTZ basis set. Recently, a similar accuracy has been reached by the explicit inclusion of derivative discontinuity in the GGA exchange potential (Carmona-Espíndola et al., ). In this section, we extend our calculation to 26 new molecules that are beyond the above mentioned works, and hence expand the test set to a total of 64 molecules. The molecular geometries are optimized at the B3LYP/6-311+G** level with the Gaussian09 suite of programs (Frisch et al., ). For the GSC approach, the relaxation of KS/GKS orbitals is considered up to second-order for GSC-B3LYP, and to first-order for other GSC-DFAs, respectively.
Table 2 lists the experimental and calculated EAs of the newly added 26 molecules. The experimental data are extracted from Jordan and Burrow (), Chiu et al. (), and Ng et al. (), while the theoretical data take either the values of −εLUMO (or −εa if it is the ath virtual orbital that is related to the resonant state) or the energy difference between the neutral and anionic species (the ΔSCF method). More details are given in the Supplementary Material.
Table 2
| Molecule | Exp. | LDA | BLYP | B3LYP | PBE | GSC- LDA | GSC- BLYP | GSC- B3LYP | GSC- PBE | ΔSCF- PBE | ΔSCF- B3LYP |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Monofluoroethylene | −1.91 | 1.25 | 0.95 | 0.31 | 1.03 | −1.97 | −2.18 | −2.06 | −2.15 | −0.47 | −0.52 |
| trans-1,2-difluoroethylene | −1.84 | 1.32 | 1.01 | 0.39 | 1.07 | −1.95 | −2.18 | −2.09 | −2.16 | −0.50 | −0.58 |
| cis-1,2-difluoroethylenea | −2.18 | 1.14 | 0.82 | 0.19 | 0.89 | −2.18 | −2.42 | −2.34 | −2.40 | −0.36 | −0.40 |
| 1,1-Difluoroethylenea | −2.39 | 1.09 | 0.79 | 0.17 | 0.84 | −2.15 | −2.32 | −2.09 | −2.33 | −0.42 | −0.47 |
| Trifluoroethylenea | −2.45 | 1.06 | 0.74 | 0.11 | 0.77 | −2.30 | −2.53 | −2.38 | −2.53 | −0.41 | −0.45 |
| Tetrafluoroethyleneb | −3.00 | 0.82 | 0.47 | 0.21 | 0.49 | −2.71 | −3.04 | −3.07 | −3.03 | −0.89 | −0.91 |
| Nitrogen | −2.20 | 2.18 | 1.88 | 0.98 | 1.92 | −2.20 | −2.41 | −2.40 | −2.41 | N/A | −1.83 |
| Formaldehyde | −0.86 | 2.91 | 2.57 | 1.75 | 2.66 | −0.92 | −1.18 | −1.14 | −1.14 | N/A | −0.46 |
| Butadiene | −0.62 | 2.12 | 1.76 | 1.18 | 1.90 | −0.61 | −0.93 | −0.85 | −0.82 | N/A | −0.73 |
| Biphenyl | −0.30 | 2.04 | 1.63 | 1.13 | 1.80 | 0.01 | −0.40 | −0.39 | −0.24 | N/A | −0.37 |
| Trichloromethane | −0.35 | 2.47 | 2.22 | 1.55 | 2.31 | −0.24 | −0.44 | −0.46 | −0.38 | −0.15 | −0.26 |
| Dichlorofluoromethane | −0.96 | 1.99 | 1.74 | 1.08 | 1.81 | −0.90 | −1.06 | −1.00 | −1.04 | −0.36 | −0.44 |
| Dichlorodifluoromethane | −0.98 | 2.37 | 2.11 | 1.43 | 2.17 | −0.61 | −0.80 | −0.79 | −0.77 | −0.42 | −0.48 |
| Dichloromethane | −1.23 | 1.74 | 1.52 | 0.90 | 1.59 | −1.02 | −1.08 | −0.89 | −1.09 | −0.31 | −0.38 |
| Benzene | −1.15 | 1.44 | 1.06 | 0.50 | 1.21 | −1.18 | −1.52 | −1.48 | −1.40 | −0.36 | −0.42 |
| CO | −1.80 | 2.24 | 1.94 | 1.12 | 2.00 | −1.89 | −2.07 | −1.96 | −2.08 | −1.05 | −1.11 |
| Cyanogen | −0.58 | 3.87 | 3.48 | 2.84 | 3.61 | −0.48 | 0.12 | 0.23 | 0.23 | 0.21 | 0.29 |
| Propyned | −2.95 | 0.13 | 0.02 | −0.40 | 0.01 | −2.29 | −1.83 | −1.23 | −2.13 | −0.40 | −0.47 |
| Butadiyne | −1.00 | 2.09 | 1.73 | 1.16 | 1.87 | −0.77 | −1.07 | −0.96 | −0.97 | −0.25 | −0.36 |
| Tetramethylethylenee | −2.27 | 0.42 | 0.20 | −0.31 | 0.28 | −1.81 | −1.79 | −1.48 | −1.86 | −0.34 | −0.41 |
| Acetylenec | −2.60 | 0.57 | 0.32 | −0.28 | 0.39 | −2.51 | −2.58 | −2.36 | −2.61 | −0.46 | −0.53 |
| Acrylonitrile | −0.21 | 3.00 | 2.62 | 2.01 | 2.76 | −0.01 | −0.35 | −0.27 | −0.24 | 0.02 | −0.16 |
| 1,4-Cyclohexadienea | −1.75 | 1.05 | 0.67 | 0.34 | 0.81 | −1.45 | −1.80 | −1.89 | −1.68 | −0.34 | −0.56 |
| Toluene | −1.11 | 1.39 | 1.01 | 0.46 | 1.16 | −1.11 | −1.45 | −1.43 | −1.33 | −0.34 | −0.42 |
| Ethylbenzene | −1.17 | 1.37 | 0.99 | 0.47 | 1.14 | −1.07 | −1.43 | −0.90 | −1.31 | −0.28 | −0.37 |
| Isopropylbenzene | −1.08 | 1.39 | 0.99 | 0.44 | 1.16 | −1.01 | −1.38 | −1.41 | −1.25 | −0.26 | −0.34 |
| MAE | 3.17 | 2.85 | 2.23 | 2.95 | 0.17 | 0.25 | 0.30 | 0.21 | 1.22 | 1.01 |
Experimental and calculated electron affinities (EAs) of 26 new molecules that are not included in Tozer and De Proft () and Zhang et al. ().
The calculated EAs are obtained by using the ΔSCF method, or take the values of the uncorrected or scaling corrected . All energies are in units of eV. The aug-cc-pVTZ basis set is adopted for all the calculated data listed in this table.
The −εLUMO+1 calculated with the GSC-B3LYP is taken as the EA of this molecule.
The −εLUMO+1 calculated with the GSC-B3LYP, GSC-BLYP, and GSC-PBE are taken as the EA of this molecule.
The −εLUMO+2 calculated with the GSC-B3LYP is taken as the EA of this molecule.
The −εLUMO+2 calculated with the GSC-B3LYP and −εLUMO+1 calculated with other DFAs are taken as the EA of this molecule.
The −εLUMO+3 calculated with the GSC-B3LYP and −εLUMO+1 calculated with other DFAs are taken as the EA of this molecule.
Figure 5 visualizes the MAEs of the calculated EAs of the extended set of 64 molecules. Obviously, the application of the GSC approach greatly improves the accuracy of the virtual orbital energies (particularly the εLUMO). The MAE is reduced from several eVs with the original DFAs to less than 0.5 eV with the GSC-DFAs. Moreover, the MAE is further reduced by adopting a more diffuse basis set. This is because a more complete basis set is more favorable for a perturbative treatment of scaling correction and orbital relaxation. The lowest MAE reached for the whole extended set is 0.14 eV with the GSC-LDA.
Figure 5
As already been pointed out in Zhang et al. (), the use of a very diffuse basis set (such as aug-cc-pVTZ) may give rise to highly delocalized virtual KS/GKS orbitals with energies close to the molecular chemical potential. These orbitals are actually not relevant to the resonant state of the temporary anion of our interest, and should be left out of theoretical analysis. Therefore, we need to choose carefully the virtual orbital, which is genuinely pertinent to the formation of the temporary anion. For instance, in the case of a cis-butene molecule, the few lowest virtual orbitals calculated at the B3LYP/cc-pVTZ and B3LYP/aug-cc-pVTZ levels are depicted in Figure 6. Apparently, with the B3LYP/aug-cc-pVTZ method, the three lowest virtual orbitals (from LUMO to LUMO+2) are rather diffuse. Occupation on any of these orbitals by an excess electron will lead to an unbound state. Therefore, these orbitals are not relevant to the formation of the temporary anion. By scrutinizing the spatial distribution of the virtual KS/GKS orbitals, it is recognized that ϕLUMO+3(r) would give rise to the resonant state of the temporary anion, as it exhibits a same shape as ϕLUMO(r) obtained with the cc-pVTZ basis set. In such a case, instead of using Equation (16), the EA is predicted by . Similarly, with the B3LYP/aug-cc-pVTZ method there are some other molecules for which a virtual orbital other than the LUMO should be chosen. The virtual orbital pertinent to the temporary anion is ϕLUMO+1(r) for 9 molecules (aniline, propene, CO2, guanine, 1,4-cyclohexadiene, cis-1,2-difluoroethylene, 1,1-difluoroethylene, trifluoroethylene and tetrafluoroethylene), ϕLUMO+2(r) for 3 molecules (trimethylethylene, propyne and acetylene), ϕLUMO+3(r) for 3 molecules (pyrrole, trans-butene, and tetramethylethylene), and ϕLUMO+4(r) for one molecule (cyclohexene).
Figure 6
As shown in Figure 5, unlike the QE-DFT method, increasing the size of basis set does not improve the accuracy of the ΔSCF method. This is because through the SCF calculation of the molecular anion by using a diffuse basis set, the excess electron is more inclined to reside on the delocalized orbital, which has a lower energy. Consequently, it is difficult to have the excess electron correctly occupying the virtual orbital that is pertinent to the resonant state of temporary anion. In contrast, the QE-DFT method in conjunction with the GSC approach does not require an SCF calculation for the anionic species, and is clearly more favorable for the prediction of resonance energies of temporary anions.
4. Conclusion
To summarize, we have calculated the QP, excitation, and resonance energies of molecules by employing the QE-DFT method. A non-empirical GSC approach is used to reduce the delocalization error associated with the DFAs by imposing an energy linearity condition for systems with a fractional number of electrons. The accuracy of the results obtained in this work with the GSC-DFAs is overall similar to that achieved in Mei et al. () by the LOSC method (Li et al., ). For instance, the GSC-B3LYP yields an MAE of 0.36 eV and a mean sign error (MSE) of −0.16 eV for the 48 excitation energies of 16 molecules (see section 3.2). These errors are slightly smaller than the MAE of 0.49 eV and the MSE of −0.19 eV resulted by the LOSC-B3LYP method (Mei et al., ). The marginal superiority in the performance of the GSC is because of the explicit treatment of the relaxation of KS/GKS orbitals upon electron addition or removal. Relaxation of KS/GKS orbitals and electron density could also be included in the LOSC calculations (Su et al., ), which will further improve the accuracy of predicted QP energies. Moreover, as the size of the molecule increases to a certain extent, the energies at integer electron numbers may become less accurate. For such systems, the correction offered by the GSC approach may be inadequate, and the LOSC method with size-consistent corrections (Su et al., ; Yang et al., ) to DFA should be used.
For the various DFAs considered in this paper, the GSC-B3LYP yields the overall best performance. Our calculation results achieve at least the same level of accuracy as some more expensive methods, such as the GW method for QP energies, the TDDFT method for excitation energies, and the EOM-CC method for resonance energies. This thus affirms that it is entirely possible and practical to study excited-state properties within the framework of ground-state DFT. Despite the promising results, the prediction of singlet excitation energies still has plenty of room for improvement. This is because another source of error associated with the DFAs, the static correlation error, comes into play, which may be corrected by imposing a constancy condition on systems with fractional spins (Cohen et al., ). Further work is needed along this direction.
Statements
Data availability statement
All datasets generated for this study are included in the article/Supplementary Material.
Author contributions
XZ and WY conceived the project. XY conducted the numerical calculations. XY, XZ, and WY analyzed the data and WY wrote the paper. All authors contributed to the article and approved the submitted version.
Funding
XY and XZ acknowledge the support from the Ministry of Science and Technology of China (Grant Nos. 2016YFA0400900 and 2016YFA0200600), the National Natural Science Foundation of China (Grant No. 21973086), and the Ministry of Education of China (111 Project Grant No. B18051). WY has been supported by the National Science Foundation (CHE-1900338) and the National Institute of General Medical Sciences of National Institute of Health (R01-GM061870).
Acknowledgments
The support from the Supercomputing Center of University of Science and Technology of China is appreciated.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.588808/full#supplementary-material
References
1
AnderssonK.MalmqvistP. A.RoosB. O.SadlejA. J.WolinskiK. (1990). Second-order perturbation theory with a CASSCF reference function. J. Phys. Chem. 94, 5483–5488. 10.1021/j100377a012
2
AulburW. G.JönssonL.WilkinsJ. W. (2000). Quasiparticle calculations in solids. Solid State Phys. 54, 1–218. 10.1016/S0081-1947(08)60248-9
3
BaerendsE.GritsenkoO.Van MeerR. (2013). The Kohn-Sham gap, the fundamental gap and the optical gap: the physical meaning of occupied and virtual Kohn-Sham orbital energies. Phys. Chem. Chem. Phys. 15, 16408–16425. 10.1039/c3cp52547c
4
BartlettR. J. (2009). Towards an exact correlated orbital theory for electrons. Chem. Phys. Lett. 484, 1–9. 10.1016/j.cplett.2009.10.053
5
BartlettR. J.RanasingheD. S. (2017). The power of exact conditions in electronic structure theory. Chem. Phys. Lett. 669, 54–70. 10.1016/j.cplett.2016.12.017
6
BeckeA. D. (1988). Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A38:3098. 10.1103/PhysRevA.38.3098
7
BeckeA. D. (1993). Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. 10.1063/1.464913
8
BrundleC.RobinM.KueblerN. (1972). Perfluoro effect in photoelectron spectroscopy. II. Aromatic molecules. J. Am. Chem. Soc. 94, 1466–1475. 10.1021/ja00760a008
9
Carmona-EspíndolaJ.GázquezJ. L.VelaA.TrickeyS. B. (2020). Negative electron affinities and derivative discontinuity contribution from a generalized gradient approximation exchange functional. J. Phys. Chem. A124, 1334–1342. 10.1021/acs.jpca.9b10956
10
CasidaM. E. (1995). Time-dependent density functional response theory for molecules, in Recent Advances in Density Functional Methods: (Part I), ed ChongD. P. (Singapore: World Scientific), 155–192.
11
ChiuN. S.BurrowP. D.JordanK. D. (1979). Temporary anions of the fluoroethylenes. Chem. Phys. Lett. 68, 121–126. 10.1016/0009-2614(79)80082-2
12
ChongD. P.GritsenkoO. V.BaerendsE. J. (2002). Interpretation of the Kohn-Sham orbital energies as approximate vertical ionization potentials. J. Chem. Phys. 116, 1760–1772. 10.1063/1.1430255
13
CohenA. J.Mori-SánchezP.YangW. (2008a). Fractional charge perspective on the band gap in density-functional theory. Phys. Rev. B77:115123. 10.1103/PhysRevB.77.115123
14
CohenA. J.Mori-SánchezP.YangW. (2008b). Fractional spins and static correlation error in density functional theory. J. Chem. Phys. 129:121104. 10.1063/1.2987202
15
CohenA. J.Mori-SánchezP.YangW. (2008c). Insights into current limitations of density functional theory. Science321, 792–794. 10.1126/science.1158722
16
CohenA. J.Mori-SánchezP.YangW. (2012). Challenges for density functional theory. Chem. Rev. 112, 289–320. 10.1021/cr200107z
17
CoropceanuV.MalagoliM.da Silva FilhoD.GruhnN.BillT.BrédasJ. (2002). Hole-and electron-vibrational couplings in oligoacene crystals: intramolecular contributions. Phys. Rev. Lett. 89:275503. 10.1103/PhysRevLett.89.275503
18
DauthM.KörzdörferT.KümmelS.ZiroffJ.WiessnerM.SchöllA.et al. (2011). Orbital density reconstruction for molecules. Phys. Rev. Lett. 107:193002. 10.1103/PhysRevLett.107.193002
19
DolgounitchevaO.Díaz-TinocoM.ZakrzewskiV. G.RichardR. M.MaromN.SherrillC. D.et al. (2016). Accurate ionization potentials and electron affinities of acceptor molecules IV: electron-propagator methods. J. Chem. Theory Comput. 12, 627–637. 10.1021/acs.jctc.5b00872
20
DreuwA.WormitM. (2015). The algebraic diagrammatic construction scheme for the polarization propagator for the calculation of excited states. Wiley Interdiscipl. Rev. Comput. Mol. Sci. 5, 82–95. 10.1002/wcms.1206
21
DunningT. H.Jr. (1989). Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 90, 1007–1023. 10.1063/1.456153
22
DuttaA. K.GuptaJ.PathakH.VavalN.PalS. (2014). Partitioned EOMEA-MBPT (2): an efficient N5 scaling method for calculation of electron affinities. J. Chem. Theory Comput. 10, 1923–1933. 10.1021/ct4009409
23
DvorakM.ChenX.-J.WuZ. (2014). Quasiparticle energies and excitonic effects in dense solid hydrogen near metallization. Phys. Rev. B90:035103. 10.1103/PhysRevB.90.035103
24
EssD. H.JohnsonE. R.HuX.YangW. (2011). Singlet-triplet energy gaps for diradicals from fractional-spin density-functional theory. J. Phys. Chem. A115, 76–83. 10.1021/jp109280y
25
FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.et al. (2009). Gaussian 09, RevisionA.01. Wallingford, CT: Gaussian, Inc.
26
GovoniM.GalliG. (2018). GW100: Comparison of methods and accuracy of results obtained with the WEST code. J. Chem. Theory Comput. 14, 1895–1909. 10.1021/acs.jctc.7b00952
27
GritsenkoO.BaerendsE. J. (2009). The analog of Koopmans' theorem for virtual Kohn-Sham orbital energies. Can. J. Chem. 87, 1383–1391. 10.1139/V09-088
28
GuerraM. (1990). On the use of diffuse functions for estimating negative electron affinities with LCAO methods. Chem. Phys. Lett. 167, 315–319. 10.1016/0009-2614(90)87174-P
29
GuichemerreM.ChambaudG.StollH. (2002). Electronic structure and spectroscopy of monohalides of metals of group I-B. Chem. Phys. 280, 71–102. 10.1016/S0301-0104(02)00510-4
30
HaidukeR. L. A.BartlettR. J. (2018). Communication: can excitation energies be obtained from orbital energies in a correlated orbital theory?J. Chem. Phys. 149:131101. 10.1063/1.5052442
31
HedinL. (1965). New method for calculating the one-particle Green's function with application to the electron-gas problem. Phys. Rev. 139:A796. 10.1103/PhysRev.139.A796
32
HillI.KahnA.CornilJ.Dos SantosD.BrédasJ. (2000). Occupied and unoccupied electronic levels in organic π-conjugated molecules: comparison between experiment and theory. Chem. Phys. Lett. 317, 444–450. 10.1016/S0009-2614(99)01384-6
33
HiraoK.ChanB.SongJ.-W.BhattaraiK.TewaryS. (2020). Excitation energies expressed as orbital energies of Kohn-Sham density functional theory with long-range corrected functionals. J. Comput. Chem. 41, 1368–1383. 10.1002/jcc.26181
34
HohenbergP.KohnW. (1964). Inhomogeneous electron gas. Phys. Rev. 136:B864. 10.1103/PhysRev.136.B864
35
HoyerC. E.GhoshS.TruhlarD. G.GagliardiL. (2016). Multiconfiguration pair-density functional theory is as accurate as CASPT2 for electronic excitation. J. Phys. Chem. Lett. 7, 586–591. 10.1021/acs.jpclett.5b02773
36
HuX.HuH.ZhengX.ZengX.WuP.PengD.et al. (2020). An In-house Program for QM/MM Simulations. Available online at: http://www.qm4d.info
37
HybertsenM. S.LouieS. G. (1986). Electron correlation in semiconductors and insulators: band gaps and quasiparticle energies. Phys. Rev. B34:5390. 10.1103/PhysRevB.34.5390
38
IpatovA.CordovaF.DoriolL. J.CasidaM. E. (2009). Excited-state spin-contamination in time-dependent density-functional theory for molecules with open-shell ground states. J. Mol. Struct. Theochem. 914, 60–73. 10.1016/j.theochem.2009.07.036
39
JacqueminD.DucheminI.BlaseX. (2015). 0-0 energies using hybrid schemes: Benchmarks of TD-DFT, CIS(D), ADC(2), CC2, and BSE/GW formalisms for 80 real-life compounds. J. Chem. Theory Comput. 11, 5340–5359. 10.1021/acs.jctc.5b00619
40
JacqueminD.DucheminI.BlaseX. (2017). Is the Bethe-Salpeter formalism accurate for excitation energies? Comparisons with TD-DFT, CASPT2, and EOM-CCSD. J. Phys. Chem. Lett. 8, 1524–1529. 10.1021/acs.jpclett.7b00381
41
JagauT.BravayaK. B.KrylovA. I. (2017). Extending quantum chemistry of bound states to electronic resonances. Annu. Rev. Phys. Chem. 68, 525–553. 10.1146/annurev-physchem-052516-050622
42
JanakJ. (1978). Proof that in density-functional theory. Phys. Rev. B18:7165. 10.1103/PhysRevB.18.7165
43
JordanK. D.BurrowP. D. (1978). Studies of the temporary anion states of unsaturated hydrocarbons by electron transmission spectroscopy. ChemInform10, 341–348. 10.1021/ar50129a004
44
JordanK. D.BurrowP. D. (1987). Temporary anion states of polyatomic hydrocarbons. Chem. Rev. 87, 557–588. 10.1021/cr00079a005
45
JordanK. D.VooraV. K.SimonsJ. (2014). Negative electron affinities from conventional electronic structure methods. Theor. Chem. Acc. 133:1445. 10.1007/s00214-014-1445-1
46
KendallR. A.DunningT. H.Jr.HarrisonR. J. (1992). Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 96, 6796–6806. 10.1063/1.462569
47
KnightJ. W.WangX.GallandiL.DolgounitchevaO.RenX.OrtizJ. V.et al. (2016). Accurate ionization potentials and electron affinities of acceptor molecules III: a benchmark of GW methods. J. Chem. Theory Comput. 12, 615–626. 10.1021/acs.jctc.5b00871
48
KohnW.ShamL. J. (1965). Self-consistent equations including exchange and correlation effects. Phys. Rev. 140:A1133. 10.1103/PhysRev.140.A1133
49
KörzdörferT.ParrishR. M.MaromN.SearsJ. S.SherrillC. D.BrédasJ.-L. (2012). Assessment of the performance of tuned range-separated hybrid density functionals in predicting accurate quasiparticle spectra. Phys. Rev. B86:205110. 10.1103/PhysRevB.86.205110
50
LaurentA. D.JacqueminD. (2013). TD-DFT benchmarks: a review. Int. J. Quant. Chem. 113, 2019–2039. 10.1002/qua.24438
51
LeeC.YangW.ParrR. G. (1988). Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B37:785. 10.1103/PhysRevB.37.785
52
LiC.ZhengX.CohenA. J.Mori-SánchezP.YangW. (2015). Local scaling correction for reducing delocalization error in density functional approximations. Phys. Rev. Lett. 114:053001. 10.1103/PhysRevLett.114.053001
53
LiC.ZhengX.SuN. Q.YangW. (2017). Localized orbital scaling correction for systematic elimination of delocalization error in density functional approximations. Nat. Sci. Rev. 5, 203–215. 10.1093/nsr/nwx111
54
LongoR.ChampagneB.ÖhrnY. (1995). Electron propagator theory and application. Theoret. Chim. Acta90, 397–419. 10.1007/BF01113544
55
LouieS. G.HybertsenM. S. (1987). Theory of quasiparticle energies: band gaps and excitation spectra in solids. Int. J. Quant. Chem. 32, 31–44. 10.1002/qua.560320706
56
MaF.WangZ.GuoM.WangF. (2020). Approximate equation-of-motion coupled-cluster methods for electron affinities of closed-shell molecules. J. Chem. Phys. 152:124111. 10.1063/1.5142736
57
MeiY.LiC.SuN. Q.YangW. (2019). Approximating quasiparticle and excitation energies from ground state generalized Kohn-Sham calculations. J. Phys. Chem. A123, 666–673. 10.1021/acs.jpca.8b10380
58
MeiY.YangW. (2019). Excited-state potential energy surfaces, conical intersections, and analytical gradients from ground-state density functional theory. J. Phys. Chem. Lett. 10, 2538–2545. 10.1021/acs.jpclett.9b00712
59
NgL.BalajiV.JordanK. D. (1983). Measurement of the vertical electron affinities of cyanogen and 2,4-hexadiyne. Chem. Phys. Lett. 101, 171–176. 10.1016/0009-2614(83)87365-5
60
NooijenM.BartlettR. J. (1995). Equation of motion coupled cluster method for electron attachment. J. Chem. Phys. 102, 3629–3647. 10.1063/1.468592
61
OnidaG.ReiningL.RubioA. (2002). Electronic excitations: density-functional versus many-body Green-function approaches. Rev. Modern Phys. 74:601. 10.1103/RevModPhys.74.601
62
OrtizJ. V. (2013). Electron propagator theory: an approach to prediction and interpretation in quantum chemistry. Wiley Interdiscipl. Rev. Comput. Mol. Sci. 3, 123–142. 10.1002/wcms.1116
63
ParrR. G.YangW. (1984). Density functional approach to the frontier-electron theory of chemical reactivity. J. Am. Chem. Soc. 106, 4049–4050. 10.1021/ja00326a036
64
PerdewJ. P.BurkeK.ErnzerhofM. (1996). Generalized gradient approximation made simple. Phys. Rev. Lett. 77:3865. 10.1103/PhysRevLett.77.3865
65
PerdewJ. P.ParrR. G.LevyM.BalduzJ. L.Jr. (1982). Density-functional theory for fractional particle number: derivative discontinuities of the energy. Phys. Rev. Lett. 49:1691. 10.1103/PhysRevLett.49.1691
66
PerdewJ. P.RuzsinszkyA.CsonkaG. I.VydrovO. A.ScuseriaG. E.StaroverovV. N.et al. (2007). Exchange and correlation in open systems of fluctuating electron number. Phys. Rev. A76:040501. 10.1103/PhysRevA.76.040501
67
PottsD. M.TaylorC. M.ChaudhuriR. K.FreedK. F. (2001). The improved virtual orbital-complete active space configuration interaction method, a “packageable” efficient ab initio many-body method for describing electronically excited states. J. Chem. Phys. 114, 2592–2600. 10.1063/1.1337053
68
PuschnigP.BoeseA.WillenbockelM.MeyerM.LüftnerD.ReinischE.et al. (2017). Energy ordering of molecular orbitals. J. Phys. Chem. Lett. 8, 208–213. 10.1021/acs.jpclett.6b02517
69
RanasingheD. S.MargrafJ. T.JinY.BartlettR. J. (2017). Does the ionization potential condition employed in QTP functionals mitigate the self-interaction error?J. Chem. Phys. 146, 034102–034102. 10.1063/1.4973727
70
RappoportD.FurcheF. (2010). Property-optimized Gaussian basis sets for molecular response calculations. J. Chem. Phys. 133:134105. 10.1063/1.3484283
71
RohlfingM.LouieS. G. (2000). Electron-hole excitations and optical spectra from first principles. Phys. Rev. B62:4927. 10.1103/PhysRevB.62.4927
72
RungeE.GrossE. K. (1984). Density-functional theory for time-dependent systems. Phys. Rev. Lett. 52:997. 10.1103/PhysRevLett.52.997
73
SancheL.SchulzG. (1972). Electron transmission spectroscopy: rare gases. Phys. Rev. A5:1672. 10.1103/PhysRevA.5.1672
74
SantoroF.JacqueminD. (2016). Going beyond the vertical approximation with time-dependent density functional theory. Wiley Interdiscipl. Rev. Comput. Mol. Sci. 6, 460–486. 10.1002/wcms.1260
75
SchäferA.HuberC.AhlrichsR. (1994). Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr. J. Chem. Phys. 100, 5829–5835. 10.1063/1.467146
76
SchmidtW. (1977). Photoelectron spectra of polynuclear aromatics. V. Correlations with ultraviolet absorption spectra in the catacondensed series. J. Chem. Phys. 66, 828–845. 10.1063/1.433961
77
SchreiberM.Silva-JuniorM. R.SauerS. P.ThielW. (2008). Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3. J. Chem. Phys. 128:134110. 10.1063/1.2889385
78
SchulzG. J. (1973). Resonances in electron impact on atoms. Rev. Modern Phys. 45:378. 10.1103/RevModPhys.45.378
79
Silva-JuniorM. R.SchreiberM.SauerS. P.ThielW. (2008). Benchmarks for electronically excited states: time-dependent density functional theory and density functional theory based multireference configuration interaction. J. Chem. Phys.129:104103. 10.1063/1.2973541
80
SkomorowskiW.GulaniaS.KrylovA. I. (2018). Bound and continuum-embedded states of cyanopolyyne anions. Phys. Chem. Chem. Phys. 20, 4805–4817. 10.1039/C7CP08227D
81
SlaterJ. C. (1951). A simplification of the Hartree-Fock method. Phys. Rev. 81:385. 10.1103/PhysRev.81.385
82
SlavicekP.MartinezT. J. (2010). Ab initio floating occupation molecular orbital-complete active space configuration interaction: an efficient approximation to CASSCF. J. Chem. Phys. 132, 234102–234102. 10.1063/1.3436501
83
StantonJ. F.BartlettR. J. (1993). The equation of motion coupled-cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties. J. Chem. Phys. 98, 7029–7039. 10.1063/1.464746
84
SuN. Q.MahlerA.YangW. (2020). Preserving symmetry and degeneracy in the localized orbital scaling correction approach. J. Chem. Phys. Lett. 11, 1528–1535. 10.1021/acs.jpclett.9b03888
85
ThierbachA.NeissC.GallandiL.MaromN.KörzdörferT.GörlingA. (2017). Accurate valence ionization energies from Kohn-Sham eigenvalues with the help of potential adjustors. J. Chem. Theory Comput. 13, 4726–4740. 10.1021/acs.jctc.7b00490
86
TozerD. J.De ProftF. (2005). Computation of the hardness and the problem of negative electron affinities in density functional theory. J. Phys. Chem. A109, 8923–8929. 10.1021/jp053504y
87
TsunedaT.SongJ.-W.SuzukiS.HiraoK. (2010). On Koopmans' theorem in density functional theory. J. Chem. Phys. 133:174101. 10.1063/1.3491272
88
VargasR.GarzaJ.CedilloA. (2005). Koopmans-like approximation in the Kohn-Sham method and the impact of the frozen core approximation on the computation of the reactivity parameters of the density functional theory. J. Phys. Chem. A109, 8880–8892. 10.1021/jp052111w
89
VoskoS. H.WilkL.NusairM. (1980). Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis. UC Irvine58, 1200–1211. 10.1139/p80-159
90
WangZ.TuZ.WangF. (2014). Equation-of-motion coupled-cluster theory for excitation energies of closed-shell systems with spin-orbit coupling. J. Chem. Theory Comput. 10, 5567–5576. 10.1021/ct500854m
91
WinterN. O.GrafN. K.LeutwylerS.HättigC. (2013). Benchmarks for 0-0 transitions of aromatic organic molecules: DFT/B3LYP, ADC (2), CC2, SOS-CC2 and SCS-CC2 compared to high-resolution gas-phase data. Phys. Chem. Chem. Phys. 15, 6623–6630. 10.1039/C2CP42694C
92
WoonD. E.DunningT. H.Jr (1993). Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 98, 1358–1371. 10.1063/1.464303
93
YangW.CohenA. J.Mori-SánchezP. (2012). Derivative discontinuity, bandgap and lowest unoccupied molecular orbital in density functional theory. J. Chem. Phys. 136:204111. 10.1063/1.3702391
94
YangW.ParrR. G. (1985). Hardness, softness, and the Fukui function in the electronic theory of metals and catalysis. Proc. Natl. Acad. Sci. U.S.A. 82, 6723–6726. 10.1073/pnas.82.20.6723
95
YangW.ParrR. G.PucciR. (1984). Electron density, Kohn-Sham frontier orbitals, and Fukui functions. J. Chem. Phys. 81, 2862–2863. 10.1063/1.447964
96
YangW.ZhangY.AyersP. W. (2000). Degenerate ground states and a fractional number of electrons in density and reduced density matrix functional theory. Phys. Rev. Lett. 84:5172. 10.1103/PhysRevLett.84.5172
97
YangX.HeZ.ZhengX. (2020). Unit cell consistency of maximally localized Wannier functions. Electron. Struct. 2:014001. 10.1088/2516-1075/ab5e5a
98
ZhangD.YangX.ZhengX.YangW. (2018). Accurate density functional prediction of molecular electron affinity with the scaling corrected Kohn-Sham frontier orbital energies. Mol. Phys. 116, 927–934. 10.1080/00268976.2017.1382738
99
ZhangD.ZhengX.LiC.YangW. (2015). Orbital relaxation effects on Kohn-Sham frontier orbital energies in density functional theory. J. Chem. Phys. 142:154113. 10.1063/1.4918347
100
ZhengX.CohenA. J.Mori-SánchezP.HuX.YangW. (2011). Improving band gap prediction in density functional theory from molecules to solids. Phys. Rev. Lett. 107:026403. 10.1103/PhysRevLett.107.026403
101
ZhengX.LiC.ZhangD.YangW. (2015). Scaling correction approaches for reducing delocalization error in density functional approximations. Sci. China Chem. 58, 1825–1844. 10.1007/s11426-015-5501-z
102
ZhengX.ZhouT.YangW. (2013). A nonempirical scaling correction approach for density functional methods involving substantial amount of Hartree-Fock exchange. J. Chem. Phys. 138:174105. 10.1063/1.4801922
103
ZieglerT.RaukA.BaerendsE. J. (1977). On the calculation of multiplet energies by the Hartree-Fock-Slater method. Theor. Chem. Acc. 43, 261–271. 10.1007/BF00551551
Summary
Keywords
density functional theory, delocalization error, scaling correction approach, quasiparticle energies, electronic excitation energies, electron affinity
Citation
Yang X, Zheng X and Yang W (2020) Density Functional Prediction of Quasiparticle, Excitation, and Resonance Energies of Molecules With a Global Scaling Correction Approach. Front. Chem. 8:588808. doi: 10.3389/fchem.2020.588808
Received
29 July 2020
Accepted
23 September 2020
Published
08 December 2020
Volume
8 - 2020
Edited by
Mohan Chen, Peking University, China
Reviewed by
Bingbing Suo, Northwest University, China; Jun Cheng, Xiamen University, China
Updates
Copyright
© 2020 Yang, Zheng and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiao Zheng xz58@ustc.edu.cnWeitao Yang weitao.yang@duke.edu
This article was submitted to Theoretical and Computational Chemistry, a section of the journal Frontiers in Chemistry
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.