Performance and Self-Consistency of the Generalized Dielectric Dependent Hybrid Functional

Exceptional Spatial Variation of Charge Injection Energies on Plasmonic Surfaces

Charge injection into a molecule on a metallic interface is a key step in many photoactivated reactions. We employ the many-body perturbation theory and compute the hole and electron injection energies for CO2 molecule on an Au nanoparticle with ∼3,000 electrons and compare it to results for idealized infinite surfaces. We demonstrate a surprisingly large variation of the injection energy barrier depending on the precise molecular position on the surface. Multiple "hot spots," characterized by low energy barriers, arise due to the competition between the plasmonic coupling and the degree of hybridization between the molecule and the substrate. The charge injection barrier to the adsorbate on the nanoparticle surface decreases from the facet edge to the facet center. This finding contrasts with the typical picture in which the electric field enhancement on the nanoparticle edges is considered the most critical factor.

Nonempirical Range-Separated Hybrid Functional with Spatially Dependent Screened Exchange

Electronic structure calculations based on density functional theory (DFT) have successfully predicted numerous ground-state properties of a variety of molecules and materials. However, exchange and correlation functionals currently used in the literature, including semilocal and hybrid functionals, are often inaccurate to describe the electronic properties of heterogeneous solids, especially systems composed of building blocks with large dielectric mismatch. Here, we present a dielectric-dependent range-separated hybrid functional, screened-exchange range-separated hybrid (SE-RSH), for the investigation of heterogeneous materials. We define a spatially dependent fraction of exact exchange inspired by the static Coulomb-hole and screened-exchange (COHSEX) approximation used in many-body perturbation theory, and we show that the proposed functional accurately predicts the electronic structure of several nonmetallic interfaces, three- and two-dimensional, pristine, and defective solids and nanoparticles.

Charge Localization in Acene Crystals from Ab Initio Electronic Structure

The performance of Koopmans-compliant hybrid functionals in reproducing the electronic structure of organic crystals is tested for a series of acene crystals. The calculated band gaps are found to be consistent with those achieved with the GW method at a fraction of the computational cost and in excellent accord with the experimental results at room temperature, when including the thermal renormalization. The energetics of excess holes and electrons reveals a struggle between polaronic localization and band-like delocalization. The consequences of these results on the transport properties of acene crystals are discussed.

Is a thin mechanism appropriate for aromatic nitration?

The mechanism of toluene nitration by NO₂BF₄ in dichloromethane solution is investigated by performing advanced ab initio MD simulations of the reaction trajectories, including at full quantum mechanical level the effects of both the solvent and of the counterion. The time evolution of the encounter complex, as well as that of the associated electronic structure, for different trajectories reveals that a single electron transfer step fastly occurs after reactants are accommodated in a common solvation shell, always preceding the formation of the σ-complex. The present results strongly suggest that the regioselectivity of the reaction is spin-density driven and that a thin mechanism, one based on reaction intermediates and transition states, can be appropriate to describe aromatic nitration.

Efficient and improved prediction of the band offsets at semiconductorheterojunctions from meta-GGA density functionals: a benchmark study

Accurate theoretical prediction of the band offsets at interfaces of semiconductor heterostructures can of-ten be quite challenging. Although density functional theory has been reasonably successful to carry outsuch calculations and efficient and accurate semilocal functionals are desirable to reduce the computational cost. In general, the semilocal functionals based on the generalized gradient approximation (GGA) significantly underestimate the bulk band gaps. This, in turn, results in inaccurate estimates of the band offsets at the heterointerfaces. In this paper, we investigate the performance of several advanced meta-GGA functionals in the computational prediction of band offsets at semiconductor heterojunctions. In particular, we investigate the performance of r 2 SCAN (revised strongly-constrained and appropriately-normed functional), rMGGAC (revised semilocal functional based on cuspless hydrogen model and Pauli kinetic energy density functional), mTASK (modified Aschebrock and Kümmel meta-GGA functional), and LMBJ (local modified Becke-Johnson) exchange-correlation functionals. Our results strongly suggest that these meta-GGA functionals for supercell calculations perform quite well, especially, when compared to computationally more demanding GW calculations. We also present band offsets calculated using ionization potentials and electron affinities, as well as band alignment via the branch point energies. Overall, our study shows that the aforementioned meta-GGA functionals can be used within the DFT framework to estimate the band offsets in semiconductor heterostructures with predictive accuracy.

Koopmans Spectral Functionals in Periodic Boundary Conditions

Koopmans spectral functionals aim to describe simultaneously ground-state properties and charged excitations of atoms, molecules, nanostructures, and periodic crystals. This is achieved by augmenting standard density functionals with simple but physically motivated orbital-density-dependent corrections. These corrections act on a set of localized orbitals that, in periodic systems, resemble maximally localized Wannier functions. At variance with the original, direct supercell implementation (Phys. Rev. X 2018, 8, 021051), we discuss here (i) the complex but efficient formalism required for a periodic boundary code using explicit Brillouin zone sampling and (ii) the calculation of the screened Koopmans corrections with density functional perturbation theory. In addition to delivering improved scaling with system size, the present development makes the calculation of band structures with Koopmans functionals straightforward. The implementation in the open-source Quantum ESPRESSO distribution and the application to prototypical insulating and semiconducting systems are presented and discussed.

Solid-state performance of a meta-GGA screened hybrid density functional constructed from Pauli kinetic enhancement factor dependent semilocal exchange hole

The semilocal form of the exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, that depends only on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is used subsequently to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional.

Green's Function Formulation of Quantum Defect Embedding Theory

We present a Green's function formulation of the quantum defect embedding theory (QDET) where a double counting scheme is rigorously derived within the G₀W₀ approximation. We then show the robustness of our methodology by applying the theory with the newly derived scheme to several defects in diamond. Additionally, we discuss a strategy to obtain converged results as a function of the size and composition of the active space. Our results show that QDET is a promising approach to investigate strongly correlated states of defects in solids.

vdW-DF-ahcx: a range-separated van der Waals density functional hybrid

Hybrid density functionals replace a fraction of an underlying generalized-gradient approximation (GGA) exchange description with a Fock-exchange component. Range-separated hybrids (RSHs) also effectively screen the Fock-exchange component and thus open the door for characterizations of metals and adsorption at metal surfaces. The RSHs are traditionally based on a robust GGA, such as PBE (Perdew J Pet al1996Phys. Rev. Lett.773865), for example, as implemented in the HSE design (Heyd Jet al2003J. Chem. Phys.1188207). Here we define an analytical-hole (Henderson T Met al2008J. Chem. Phys.128194105) consistent-exchange RSH extension to the van der Waals density functional (vdW-DF) method (Berland Ket al2015Rep. Prog. Phys.78066501), launching vdW-DF-ahcx. We characterize the GGA-type exchange in the vdW-DF-cx version (Berland K and Hyldgaard P 2014Phys. Rev. B89035412), isolate the short-ranged exchange component, and define the new vdW-DF hybrid. We find that the performance vdW-DF-ahcx compares favorably to (dispersion-corrected) HSE for descriptions of bulk (broad molecular) properties. We also find that it provides accurate descriptions of noble-metal surface properties, including CO adsorption.

Computational design of quantum defects in two-dimensional materials

Missing atoms or atom substitutions (point defects) in crystal lattices in two-dimensional (2D) materials are potential hosts for emerging quantum technologies, such as single-photon emitters and spin quantum bits (qubits). First-principles-guided design of quantum defects in 2D materials is paving the way for rational spin qubit discovery. Here we discuss the frontier of first-principles theory development and the challenges in predicting the critical physical properties of point defects in 2D materials for quantum information technology, in particular for optoelectronic and spin-optotronic properties. Strong many-body interactions at reduced dimensionality require advanced electronic structure methods beyond mean-field theory. The great challenges for developing theoretical methods that are appropriate for strongly correlated defect states, as well as general approaches for predicting spin relaxation and the decoherence time of spin defects, are yet to be addressed.

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe-Salpeter equation (BSE) for electron-hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

Koopmans-Compliant Functionals and Potentials and Their Application to the GW100 Test Set

Koopmans-compliant (KC) functionals have been shown to provide accurate spectral properties through a generalized condition of piecewise linearity of the total energy as a function of the fractional addition/removal of an electron to/from any orbital. We analyze the performance of different KC functionals on a large and standardized set of 100 molecules, the GW100 test set, comparing vertical ionization potentials (taken as opposite of the orbital energies) to those obtained from accurate quantum chemistry methods, and to experimental results. We find excellent agreement, with a mean absolute error of 0.20 eV for the KIPZ functional on the first ionization potential, which is state-of-the-art for both density functional theory (DFT)-based calculations and many-body perturbation theory. We highlight similarities and differences between KC functionals and other electronic-structure approaches, such as dielectric-dependent hybrid functionals and Green's function methods, both from a theoretical and from a practical point of view, arguing that KC potentials can be considered as local and orbital-dependent approximations to the electronic self-energy, already including approximate vertex corrections.

Screened hybrid meta-GGA exchange–correlation functionals for extended systems

Screened Hartree–Fock exchange integrated with semilocal exchange–correlation functionals often proficiently predict several solid-state properties.

Dielectric Screening Meets Optimally Tuned Density Functionals

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

Optical Absorbance Enhancement in PbS QD/Cinnamate Ligand Complexes

We studied the optical absorption enhancement in colloidal suspensions of PbS quantum dots (QD) upon ligand exchange from oleate to a series of cinnamate ligands. By combining experiments and ab initio simulations, we elucidate physical parameters that govern the optical absorption enhancement. We find that, within the cinnamate/PbS QD system, the optical absorption enhancement scales linearly with the electronic gap of the ligand, indicating that the ligand/QD coupling occurs equally efficient between the QD and ligand HOMO and their respective LUMO levels. Disruption of the conjugation that connects the aromatic ring and its substituents to the QD core causes a reduction of the electronic coupling. Our results further support the notion that the ligand/QD complex should be considered as a distinct chemical system with emergent behavior rather than a QD core with ligands whose sole purpose is to passivate surface dangling bonds and prevent agglomeration.

GW100: Comparison of Methods and Accuracy of Results Obtained with the WEST Code

The reproducibility of calculations carried out within many-body perturbation theory at the G₀ W₀ level is assessed for 100 closed shell molecules and compared to that of density functional theory. We consider vertical ionization potentials (VIP) and electron affinities (VEA) obtained with five different codes: BerkeleyGW, FHI-aims, TURBOMOLE, VASP, and WEST. We review the approximations and parameters that control the accuracy of G₀ W₀ results in each code, and we discuss in detail the effect of extrapolation techniques for the parameters entering the WEST code. Differences between the VIP and VEA computed with the various codes are within ∼60 and ∼120 meV, respectively, which is up to four times larger than in the case of the best results obtained with DFT codes. Vertical ionization potentials are validated against experiment and CCSD(T) quantum chemistry results showing a mean absolute relative error of ∼4% for data obtained with WEST. Our analysis of the differences between localized orbitals and plane-wave implementations points out molecules containing Cu, I, Ga, and Xe as major sources of discrepancies, which call for a re-evaluation of the pseudopotentials used for these systems in G₀ W₀ calculations.

Accuracy of dielectric-dependent hybrid functionals in the prediction of optoelectronic properties of metal oxide semiconductors: a comprehensive comparison with many-body GW and experiments

Understanding the electronic structure of metal oxide semiconductors is crucial to their numerous technological applications, such as photoelectrochemical water splitting and solar cells. The needed experimental and theoretical knowledge goes beyond that of pristine bulk crystals, and must include the effects of surfaces and interfaces, as well as those due to the presence of intrinsic defects (e.g. oxygen vacancies), or dopants for band engineering. In this review, we present an account of the recent efforts in predicting and understanding the optoelectronic properties of oxides using ab initio theoretical methods. In particular, we discuss the performance of recently developed dielectric-dependent hybrid functionals, providing a comparison against the results of many-body GW calculations, including G ₀ W ₀ as well as more refined approaches, such as quasiparticle self-consistent GW. We summarize results in the recent literature for the band gap, the band level alignment at surfaces, and optical transition energies in defective oxides, including wide gap oxide semiconductors and transition metal oxides. Correlated transition metal oxides are also discussed. For each method, we describe successes and drawbacks, emphasizing the challenges faced by the development of improved theoretical approaches. The theoretical section is preceded by a critical overview of the main experimental techniques needed to characterize the optoelectronic properties of semiconductors, including absorption and reflection spectroscopy, photoemission, and scanning tunneling spectroscopy (STS).