Jong-Won Song

Long-range corrected density functional calculations of chemical reactions: Redetermination of parameter

Chemical reaction calculations were carried out using the long-range correction (LC) scheme, which improves long-range exchange effects in density functional theory (DFT) [J. Chem. Phys. 115, 3540 (2001); 120, 8425 (2004)]. A new determination of the LC scheme parameter mu was made by a root mean square fit of the percent error in calculated atomization energies. As a result, the parameter mu was optimized as 0.47, which is higher than the previous one (mu=0.33). Using this new parameter mu, LC-DFT was firstly applied to geometry optimizations of the G2 benchmark set molecules. Consequently, this new LC-DFT gave more accurate bond lengths and bond angles than previous LC-DFT and hybrid B3LYP results. Following this result, the authors calculated reaction barrier height energies of benchmark reaction sets, which have been underestimated in conventional DFT calculations. Calculated results showed that LC-DFT provided much more accurate barrier height energies with errors less than half those of previous LC-DFT and B3LYP studies. To test the general validity of the new LC-DFT, the authors finally calculated reaction enthalpies. As a result, they found that the LC scheme using the new mu clearly improved the accuracy of calculated enthalpies. The authors therefore conclude that the insufficient inclusion of long-range exchange effects is responsible for the underestimation of reaction barriers in DFT calculations and that LC-DFT using the new parameter is a powerful tool for theoretically investigating chemical reactions.

On Koopmans' theorem in density functional theory.

This paper clarifies why long-range corrected (LC) density functional theory gives orbital energies quantitatively. First, the highest occupied molecular orbital and the lowest unoccupied molecular orbital energies of typical molecules are compared with the minus vertical ionization potentials (IPs) and electron affinities (EAs), respectively. Consequently, only LC exchange functionals are found to give the orbital energies close to the minus IPs and EAs, while other functionals considerably underestimate them. The reproducibility of orbital energies is hardly affected by the difference in the short-range part of LC functionals. Fractional occupation calculations are then carried out to clarify the reason for the accurate orbital energies of LC functionals. As a result, only LC functionals are found to keep the orbital energies almost constant for fractional occupied orbitals. The direct orbital energy dependence on the fractional occupation is expressed by the exchange self-interaction (SI) energy through the potential derivative of the exchange functional plus the Coulomb SI energy. On the basis of this, the exchange SI energies through the potential derivatives are compared with the minus Coulomb SI energy. Consequently, these are revealed to be cancelled out only by LC functionals except for H, He, and Ne atoms.

An improved long-range corrected hybrid exchange-correlation functional including a short-range Gaussian attenuation (LCgau-BOP)

A new hybrid exchange-correlation functional is presented based on the long-range correction (LC) scheme [H. Iikura et al., J. Chem. Phys. 115, 3540 (2001); Tawada et al., J. Chem. Phys. 120, 8425 (2004)], named LCgau-BOP. The key feature is the use of a two-parameter Gaussian correction to the Coulomb attenuation, which allows a more flexible description of exact exchange at short-range interelectronic separations. The new partitioning preserves 100% exact exchange in the long range, which is known to be important for the success of the LC scheme, with an asymptotic attenuation described by a standard error function with a parameter of 0.42. The LCgau partitioning was optimized for the reproduction of atomization energies over the G2 set and reaction barrier heights over Database/3, and produced results which are superior to B3LYP, CAM-BLYP, and the best LC functionals we are aware of. The results highlight the importance of including a substantial portion of exact exchange in the short range. Using the same parameters, the new functional was tested for the reproduction of geometries, as well as valence, Rydberg and charge-transfer excitations which are known challenges for conventional density functional theory. Our conclusion is that LCgau-BOP can provide a consistently more accurate description of thermochemistries, chemical reactions, and excitation energies than other existing long-range corrected functionals.

Nonlinear optical property calculations of polyynes with long-range corrected hybrid exchange-correlation functionals

Polarizabilities (alpha), second-hyperpolarizabilities (gamma), and the gamma scaling factors (c) of polyynes [H-(C[triple bond]C)(n)-H, n = 1-8] were evaluated using the long-range corrected (LC) density functional theory (DFT) and LC-DFT with a short-range Gaussian attenuation (LCgau), as well as high quality wavefunction methods. We show that the c values obtained from LC- and LCgau-DFT are consistent with those from CCSD(T) calculations. Furthermore, the polyyne c values we obtained are seen to be smaller than the c values derived from previously calculated polyene gamma values [Sekino et al., J. Chem. Phys. 126, 014107 (2007)] in all the methods we consider. We compare our results with those obtained experimentally [Shepkov et al., J. Chem. Phys. 120, 6807 (2004).] from end-capped polyynes [i-Pr(3)Si-(C[triple bond]C)(n)-Sii-Pr(3)], which show larger c values for polyynes than polyenes. Our alpha and gamma calculations with i-Pr(3)Si-(C[triple bond]C)(n)-Sii-Pr(3) (n = 4, 5, 6, and 8) show that i-Pr(3)Si- may participate in pi molecular orbital delocalization, which can unexpectedly affect the c value. We also confirm the importance of molecular geometry in these nonlinear optical calculations. We find that while LC- and LCgau-DFT excellently reproduce experimental geometries and bond length alternation (BLA), MP2 optimized geometries have a BLA that is too short to be used for accurate alpha and gamma calculations.

Calculations of Alkane Energies Using Long-Range Corrected DFT Combined with Intramolecular van der Waals Correlation

The isodesmic reaction enthalpies of n-alkanes to ethane, which have so far been known to give systematic errors in density functional theory (DFT) calculations, are successfully reproduced using long-range corrected (LC) DFT combined with a local response dispersion (LRD) method, which reveals that the failure of conventional DFT results from the lack of long-range exchange interactions in exchange functionals as well as intramolecular dispersion and medium-range electron correlations in correlation functionals.

Core-excitation energy calculations with a long-range corrected hybrid exchange-correlation functional including a short-range Gaussian attenuation (LCgau-BOP)

We report the calculations of core-excitation energies of first-row atoms using the time-dependent density functional theory (DFT) and the long-range correction (LC) scheme for exchange-correlation functionals, including LC-BOP, Coulomb-attenuated method BLYP, and our recently developed LCgau-BOP method, which includes a flexible portion of short-range Hartree-Fock (HF) exchange through the inclusion of a Gaussian function in the LC scheme. We show that the LC scheme completely fails to improve the poor accuracy of conventional generalized gradient approximation functionals, while the LCgau scheme gives an accuracy which is an order of magnitude better than BLYP and significantly better than B3LYP. A reoptimization of the two parameters controlling the inclusion of short-range HF exchange in the LCgau method enables the errors to be reduced to the order of 0.1 eV which is competitive with the best DFT methods we are aware of. This reparametrization does not affect the LC scheme and therefore maintains the high accuracy of predicted reaction barrier heights. Moreover, while there is some loss in accuracy in thermochemical predictions compared to the previously optimized LCgau-BOP, rms errors in the atomization energies over the G2 test set are found to be comparable to B3LYP. Finally, we attempt to rationalize the success of the LC and LCgau schemes in terms of the well-known self-interaction error (SIE) of conventional functionals. To estimate the role of the SIE, we examine the total energy calculations for systems with a fractional number of electrons, not only in the highest occupied molecular orbital but also in the 1s-characterized core orbital. Our conclusion is that the inclusion of short-range HF exchange in LC-type functionals can significantly alleviate the problems of the SIE in the core region. In particular, we confirm that the absence of the SIE diagnostics in the core orbital energies correlates with the accurate prediction of core-excitation energies using the newly optimized LCgau approach.

Communication: A new hybrid exchange correlation functional for band-gap calculations using a short-range Gaussian attenuation (Gaussian-Perdue–Burke–Ernzerhof)

We have developed a new hybrid functional [Gaussian-Perdue-Burke-Ernzerhof (Gau-PBE)] that is suitable for the calculation of solid state bandgaps using a periodic boundary condition. The characteristic of this functional is the use of a Gaussian attenuation scheme (Gau) to include a short-range Hartree-Fock (HF) exchange. This new functional can perform barrier height calculations with an accuracy comparable to the middle-range hybrid functional and bandgap calculations with an accuracy comparable to the Heyd-Scuseria-Ernzerhof (HSE) functional. However, the point is that the performance can be achieved using a Gaussian HF exchange, while the HISS functional calculates twice HF-exchange integrations using an error function to obtain both performances. In addition, Gau-PBE functional can decrease the time cost of bandgap calculations by an average of 40% compared to the HSE functional.

Long-range corrected density functional theory with accelerated Hartree-Fock exchange integration using a two-Gaussian operator [LC-ωPBE(2Gau)]

Since the advent of hybrid functional in 1993, it has become a main quantum chemical tool for the calculation of energies and properties of molecular systems. Following the introduction of long-range corrected hybrid scheme for density functional theory a decade later, the applicability of the hybrid functional has been further amplified due to the resulting increased performance on orbital energy, excitation energy, non-linear optical property, barrier height, and so on. Nevertheless, the high cost associated with the evaluation of Hartree-Fock (HF) exchange integrals remains a bottleneck for the broader and more active applications of hybrid functionals to large molecular and periodic systems. Here, we propose a very simple yet efficient method for the computation of long-range corrected hybrid scheme. It uses a modified two-Gaussian attenuating operator instead of the error function for the long-range HF exchange integral. As a result, the two-Gaussian HF operator, which mimics the shape of the error function operator, reduces computational time dramatically (e.g., about 14 times acceleration in C diamond calculation using periodic boundary condition) and enables lower scaling with system size, while maintaining the improved features of the long-range corrected density functional theory.

A semiempirical long-range corrected exchange correlation functional including a short-range Gaussian attenuation (LCgau-B97)

We applied an improved long‐range correction scheme including a short‐range Gaussian attenuation (LCgau) to the Becke97 (B97) exchange correlation functional. In the optimization of LCgau‐B97 functional, the linear parameters are determined by least squares fitting. Optimizing μ parameter (0.2) that controls long‐range portion of Hartree‐Fock (HF) exchange to excitation energies of large molecules (Chai and Head‐Gordon, J Chem Phys 2008, 128, 084106) and additional short‐range Gaussian parameters (a = 0.15 and k = 0.9) that controls HF exchange inclusion ranging from short‐range to mid‐range (0.5–3 Å) to ground state properties achieved high performances of LCgau‐B97 simultaneously on both ground state and excited state properties, which is better than other tested semiempirical density functional theory (DFT) functionals, such as ωB97, ωB97X, BMK, and M0x‐family. We also found that while a small μ value (∼0.2) in LC‐DFT is appropriate to the local excitation and intramolecular charge‐transfer excitation energies, a larger μ value (0.42) is desirable in the Rydberg excitation‐energy calculations.

Gaussian attenuation hybrid scheme applied to the Ernzerhof-Perdew exchange hole model (Gau-PBEh)

Recently, we developed a Gaussian attenuation (Gau) scheme for solid-state bandgap calculation that uses a two-electron Gaussian function operator to include short-range Hartree-Fock exchange and combined it with the long-range Perdew-Burke-Ernzerhof (PBE) exchange correlation functional (Gau-PBE). Here, we apply the Ernzerhof-Perdew exchange hole (EP) model of PBE (PBEh) as a long-range density functional theory (DFT) exchange part to the Gau scheme (Gau-PBEh). We found that applying the EP model to the Gau scheme improves atomization energies and solid-state lattice constants and that the exact exchange included using the Gau scheme plays a critical role in simultaneously reproducing solid-state bandgaps and barrier heights. In addition, Gau-PBEh takes nearly the same computation time for bandgap calculations as Gau-PBE, implying less than 60% of the time taken in Heyd-Scuseria-Ernzerhof hybrid DFT functional calculations.