Artur F. Izmaylov

Influence of the exchange screening parameter on the performance of screened hybrid functionals

This work reexamines the effect of the exchange screening parameter omega on the performance of the Heyd-Scuseria-Ernzerhof (HSE) screened hybrid functional. We show that variation of the screening parameter influences solid band gaps the most. Other properties such as molecular thermochemistry or lattice constants of solids change little with omega. We recommend a new version of HSE with the screening parameter omega=0.11 bohr(-1) for further use. Compared to the original implementation, the new parametrization yields better thermochemical results and preserves the good accuracy for band gaps and lattice constants in solids.

Can short-range hybrids describe long-range-dependent properties?

Long-range-corrected hybrids, which incorporate all of the long-range exact exchange interaction, improve performance for a host of molecular properties. The long-range portion of exact exchange is both computationally and formally problematic in solids, and screened hybrids therefore eliminate it. While screened hybrids give similar results to their parent global hybrids for many molecular properties, one may worry that they perform poorly for those properties that are improved by the long-range-correction procedure. In this paper, we show that at least for the Heyd-Scuseria-Ernzerhof (HSE) screened hybrid, this is not the case; for many properties improved by long-range-correction, screened hybrids and global hybrids deliver essentially the same results. We suggest that this is because screened hybrids and global hybrids have fundamentally the same many-electron self-interaction error. We also introduce some small revisions to our computational implementation of the HSE screened hybrid, and we recommend these revisions for future applications of HSE.

Efficient evaluation of short-range Hartree-Fock exchange in large molecules and periodic systems

We present an efficient algorithm for the evaluation of short-range Hartree-Fock exchange energies and geometry gradients in Gaussian basis sets. Our method uses a hierarchy of screening levels to eliminate negligible two-electron integrals whose evaluation is the fundamental computational bottleneck of the procedure. By applying our screening technique to the Heyd-Scuseria-Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] short-range Coulomb hybrid density functional, we achieve a computational efficiency comparable with that of standard nonhybrid density functional calculations.

The importance of middle-range Hartree-Fock-type exchange for hybrid density functionals

Hybrid functionals are responsible for much of the utility of modern Kohn-Sham density functional theory. When rigorously applied to solid-state metallic and small band gap systems, however, the slow decay of their nonlocal Hartree-Fock-type exchange makes hybrids computationally challenging and introduces unphysical effects. This can be remedied by using a range-separated hybrid which only keeps short-range nonlocal exchange, as in the functional of Heyd et al. [J. Chem. Phys. 118, 8207 (2003)]. On the other hand, many molecular properties require full long-range nonlocal exchange, which can also be included by means of a range-separated hybrid such as the recently introduced LC-omegaPBE functional [O. A. Vydrov and G. E. Scuseria, J. Chem. Phys. 125, 234109 (2006)]. In this paper, we show that a three-range hybrid which mainly includes middle-range Hartree-Fock-type exchange and neglects long- and short-range Hartree-Fock-type exchange yields excellent accuracy for thermochemistry, barrier heights, and band gaps, emphasizing that the middle-range part of the 1/r potential seems crucial to accurately model these properties.

Assessment of a Middle-Range Hybrid Functional.

While hybrid functionals are largely responsible for the utility of modern Kohn-Sham density functional theory, they are not without their weaknesses. In the solid state, the slow decay of their nonlocal Hartree-Fock-type exchange makes hybrids computationally demanding and can introduce unphysical effects. Both problems can be remedied by a screened hybrid which uses exact exchange only at short-range. Many molecular properties, in contrast, benefit from the inclusion of long-range exact exchange. Recently, the authors reconciled these two seemingly contradictory requirements by introducing the HISS functional [ J. Chem. Phys. 2007 , 127 , 221103 ], which uses exact exchange only in the middle range. In this paper, we expand upon our previous work, benchmarking the performance of the HISS functional for several simple properties and applying it to the dissociation of homonuclear diatomic cations and to the polarizability of linear H2 chains to determine the importance of middle-range exact exchange for these systems, which are expected to be sensitive to the asymptotic exchange potential.

Accurate solid-state band gaps via screened hybrid electronic structure calculations.

The band energy differences of solids calculated with screened hybrid density functionals, such as the functional of Heyd-Scuseria-Ernzerhof (HSE), reproduce experimental band gaps with a high degree of accuracy. This unexpected result is here rationalized by observing that band energy differences obtained from generalized Kohn-Sham calculations with screened (short-range) Hartree-Fock-type exchange approach the excitation energies obtained via time-dependent density functional calculations with the corresponding unscreened functional. The latter are expected to be the accurate predictions of the experimental optical absorption spectra. While the optimum screening parameter (omega) is system dependent, the HSE standard value of omega=0.11 bohr(-1) represents a reasonable compromise across diverse systems.

Nonequilibrium Fermi golden rule for electronic transitions through conical intersections

We consider photoinduced electronic transitions through conical intersections in large molecules. Starting from the linear vibronic model Hamiltonian and treating linear diabatic couplings within the second order cumulant expansion, we have developed a simple analytical expression for the time evolution of electronic populations at finite temperature. The derived expression can be seen as a nonequilibrium generalization of the Fermi golden rule due to a nonequilibrium character of the initial photoinduced nuclear distribution. All parameters in our model are obtained from electronic structure calculations followed by a diabatization procedure. The results of our model are found to agree well with those of quantum dynamics for a test set of systems: fulvene molecule, 2,6-bis(methylene) adamantyl cation, and its dimethyl derivative.

Linear-scaling calculation of static and dynamic polarizabilities in Hartree-Fock and density functional theory for periodic systems

We present a linear-scaling method for analytically calculating static and dynamic polarizabilities with Hartree-Fock and density functional theory, using Gaussian orbitals and periodic boundary conditions. Our approach uses the direct space fast multipole method to evaluate the long-range Coulomb contributions. For exact exchange, we use efficient screening techniques developed for energy calculations. We then demonstrate the capabilities of our approach with benchmark calculations on one-, two-, and three-dimensional systems.

When do we need to account for the geometric phase in excited state dynamics

We investigate the role of the geometric phase (GP) in an internal conversion process when the system changes its electronic state by passing through a conical intersection (CI). Local analysis of a two-dimensional linear vibronic coupling (LVC) model Hamiltonian near the CI shows that the role of the GP is twofold. First, it compensates for a repulsion created by the so-called diagonal Born-Oppenheimer correction. Second, the GP enhances the non-adiabatic transition probability for a wave-packet part that experiences a central collision with the CI. To assess the significance of both GP contributions we propose two indicators that can be computed from parameters of electronic surfaces and initial conditions. To generalize our analysis to N-dimensional systems we introduce a reduction of a general N-dimensional LVC model to an effective 2D LVC model using a mode transformation that preserves short-time dynamics of the original N-dimensional model. Using examples of the bis(methylene) adamantyl and butatriene cations, and the pyrazine molecule we have demonstrated that their effective 2D models reproduce the short-time dynamics of the corresponding full dimensional models, and the introduced indicators are very reliable in assessing GP effects.

The effective local potential method: Implementation for molecules and relation to approximate optimized effective potential techniques

We have recently formulated a new approach, named the effective local potential (ELP) method, for calculating local exchange-correlation potentials for orbital-dependent functionals based on minimizing the variance of the difference between a given nonlocal potential and its desired local counterpart [V. N. Staroverov et al., J. Chem. Phys. 125, 081104 (2006)]. Here we show that under a mildly simplifying assumption of frozen molecular orbitals, the equation defining the ELP has a unique analytic solution which is identical with the expression arising in the localized Hartree-Fock (LHF) and common energy denominator approximations (CEDA) to the optimized effective potential. The ELP procedure differs from the CEDA and LHF in that it yields the target potential as an expansion in auxiliary basis functions. We report extensive calculations of atomic and molecular properties using the frozen-orbital ELP method and its iterative generalization to prove that ELP results agree with the corresponding LHF and CEDA values, as they should. Finally, we make the case for extending the iterative frozen-orbital ELP method to full orbital relaxation.