Atomic self-interaction correction for molecules and solids

We present an atomic orbital based approximate scheme for self-interaction correction (SIC) to the local density approximation of density functional theory. The method, based on the idea of Filippetti and Spaldin [Phys. Rev. B 67, 125109 (2003)], is implemented in a code using localized numerical atomic orbital basis sets and is now suitable for both molecules and extended solids. After deriving the fundamental equations as a non-variational approximation of the self-consistent SIC theory, we present results for a wide range of molecules and insulators. In particular, we investigate the effect of re-scaling the self-interaction correction and we establish a link with the existing atomic-like corrective scheme LDA+U. We find that when no re-scaling is applied, i.e. when we consider the entire atomic correction, the Kohn-Sham HOMO eigenvalue is a rather good approximation to the experimental ionization potential for molecules. Similarly the HOMO eigenvalues of negatively charged molecules reproduce closely the molecular affinities. In contrast a re-scaling of about 50% is necessary to reproduce insulator bandgaps in solids, which otherwise are largely overestimated. The method therefore represents a Kohn-Sham based single-particle theory and offers good prospects for applications where the actual position of the Kohn-Sham eigenvalues is important, such as quantum transport.Comment: 16 pages, 7 figure

How good are recent density functionals for ground and excited states of one-electron systems?

Sun et al. [J. Chem. Phys. 144, 191101 (2016)] suggested that common density functional approximations (DFAs) should exhibit large energy errors for excited states as a necessary consequence of orbital nodality. Motivated by self-interaction corrected density functional calculations on many-electron systems, we continue their study with the exactly solvable 1s, 2p, and 3d states of 36 hydrogenic one-electron ions (H-Kr35+) and demonstrate with self-consistent calculations that state-of-the-art DFAs indeed exhibit large errors for the 2p and 3d excited states. We consider 56 functionals at the local density approximation (LDA), generalized gradient approximation (GGA) as well as meta-GGA levels, also including several hybrid functionals like the recently proposed machine-learned DM21 local hybrid functional. The best non-hybrid functional for the 1s ground state is revTPSS. The 2p and 3d excited states are more difficult for DFAs as Sun et al. predicted, and LDA functionals turn out to yield the most systematic accuracy for these states amongst non-hybrid functionals. The best performance for the three states overall is observed with the BHandH global hybrid GGA functional, which contains 50% Hartree-Fock exchange and 50% LDA exchange. The performance of DM21 is found to be inconsistent, yielding good accuracy for some states and systems and poor accuracy for others. Based on these results, we recommend including a variety of one-electron cations in future training of machine-learned density functionals. (C) 2022 Author(s).Peer reviewe

Self-interaction corrected SCAN functional for molecules and solids in the numeric atom-center orbital framework

Das „Strongly Constrained and Appropriately Normed“ (SCAN) Austausch-Korrelations-Funktional gehört zur Familie der meta-GGA (generalized gradient approximation) Funktionale. Es gibt aber auch Nachteile Zum einen leiden SCAN Rechnungen oft unter numerischen Instabilitäten, wodurch sehr viele Iteration zum Erreichen von Selbst-Konsistenz benötigt werden. Zum anderen leidet SCAN unter dem von GGA Methoden bekannten Selbstwechselwirkung-Fehler. Im ersten Teil der Arbeit habe ich die numerischen Stabilitätsprobleme in SCAN Rechnungen im Rahmen der numerischen Realraum-Integrationsroutinen im Code FHI-aims untersucht. Diese Analyse zeigt, dass die genannte Probleme durch Anwendung von standardisierten Dichte-Mischalgorithmen für die kinetische Energiedichte abgemildert werden können. Dadurch wird auch in SCAN-Rechnungen eine schnelle und stabile Konvergenz zur selbstkonsistenten Lösung ermöglicht. Im zweiten Teil der Arbeit habe ich untersucht, in welchem Rahmen sich der Selbstwechselwirkung-Fehler in SCAN mittels des von Perdew und Zunger vorgeschlagenen Selbstinteraktionskorrekturalgorithmus (PZ-SIC) verringern lässt. Es wurden aber auch Optimierungen für die PZ-SIC Methode entwickelt. Inspiriert von den ursprünglichen Argumenten in der PZ-SIC-Methode und anderen lokalisierten Methoden, wird in dieser Arbeit eine neuartige Randbedingung (orbital density constraint) vorgeschlagen, die sicherstellt, dass die PZ-SIC Orbitale während des Selbstkonsistenzzyklus lokalisiert bleiben. Dies mildert die Anfangswertabhängigkeit deutlich ab und hilft dabei, in die korrekte selbst-konsistente Lösung mit minimaler Energie zu konvergieren, unabhängig davon ob reelle oder komplexe SIC Orbitale verwendet werden. Die in dieser Arbeit getägtigen Entwicklungen und Untersuchungen sind Wegbereiter dafür, in Zukunft mit SIC-SCAN Rechnungen deutlich genauere ab initio Rechnungen mit nur gering höherem Rechenaufwand durchführen zu können.The state-of-the-art “Strongly Constrained and Appropriately Normed” (SCAN) functional pertains to the family of meta-generalized-gradient approximation (meta-GGA) exchange-correlation functionals. Nonetheless, SCAN suffers from some well-documented deficiencies. In the first part of this thesis, I revisited the known numerical instability problems of the SCAN functional in the context of the numerical, real-space integration framework used in the FHI-aims code. This analysis revealed that applying standard density-mixing algorithms to the kinetic energy density attenuates and largely cures these numerical issues. By this means, SCAN calculations converge towards the self-consistent solution as fast and as efficiently as lower-order GGA calculations. In the second part of the thesis, I investigated strategies to alleviate the self-interaction error in SCAN calculations by using the self-interaction correction algorithm proposed by Perdew and Zunger (PZ-SIC). Inspired by the original arguments in PZ-SIC and other localized methods, I introduced a mathematical constraint, i.e., the orbital density constraint, that forces the orbitals to retain their localization throughout the self-consistency cycle. In turn, this alleviates the multiple-solutions problem and facilitates the convergence towards the correct, lowest-energy solution both for complex and real SIC orbitals. The developments and investigations performed in this thesis pave the road towards a more wide-spread use of SIC-SCAN calculations in the future, allowing more accurate predictions within only moderate increases of computational cost

Secreted CLIC3 drives cancer progression through its glutathione-dependent oxidoreductase activity

The secretome of cancer and stromal cells generates a microenvironment that contributes to tumour cell invasion and angiogenesis. Here we compare the secretome of human mammary normal and cancer-associated fibroblasts (CAFs). We discover that the chloride intracellular channel protein 3 (CLIC3) is an abundant component of the CAF secretome. Secreted CLIC3 promotes invasive behaviour of endothelial cells to drive angiogenesis and increases invasiveness of cancer cells both in vivo and in 3D cell culture models, and this requires active transglutaminase-2 (TGM2). CLIC3 acts as a glutathione-dependent oxidoreductase that reduces TGM2 and regulates TGM2 binding to its cofactors. Finally, CLIC3 is also secreted by cancer cells, is abundant in the stromal and tumour compartments of aggressive ovarian cancers and its levels correlate with poor clinical outcome. This work reveals a previously undescribed invasive mechanism whereby the secretion of a glutathione-dependent oxidoreductase drives angiogenesis and cancer progression by promoting TGM2-dependent invasion

Origin of the Pseudogap in High-Temperature Cuprate Superconductors

Cuprate high-temperature superconductors exhibit a pseudogap in the normal state that decreases monotonically with increasing hole doping and closes at x \approx 0.19 holes per planar CuO2 while the superconducting doping range is 0.05 < x < 0.27 with optimal Tc at x \approx 0.16. Using ab initio quantum calculations at the level that leads to accurate band gaps, we found that four-Cu-site plaquettes are created in the vicinity of dopants. At x \approx 0.05 the plaquettes percolate, so that the Cu dx2y2/O p{\sigma} orbitals inside the plaquettes now form a band of states along the percolating swath. This leads to metallic conductivity and below Tc to superconductivity. Plaquettes disconnected from the percolating swath are found to have degenerate states at the Fermi level that split and lead to the pseudogap. The pseudogap can be calculated by simply counting the spatial distribution of isolated plaquettes, leading to an excellent fit to experiment. This provides strong evidence in favor of inhomogeneous plaquettes in cuprates.Comment: 24 pages (4 pages main text plus 20 pages supplement

Tuning the graphene band gap by thermodynamic control of molecular self-assembly on graphene

Recent interest in functionalised graphene has been motivated by the prospect of creating a two-dimensional semiconductor with a tuneable band gap. Various approaches to band gap engineering have been made over the last decade, one of which is chemical functionalisation. However, the patterning of molecular adsorption onto graphene has proved to be difficult, as grown structures tend to be stochastic in nature. In the first part of this work, a predictive physical model of the self-assembly of halogenated carbene layers on graphene is suggested. Self-assembly of the adsorbed layer is found to be governed by a combination of the curvature of the graphene sheet, local distortions, as introduced by molecular adsorption, and short-range intermolecular repulsion. The thermodynamics of bidental covalent molecular adsorption and the resultant electronic structure are computed using Density Functional Theory. It is predicted that a direct band gap is opened that is tuneable by varying coverages and is dependent on the ripple ampli- tude. This provides a mechanism for the controlled engineering of graphene’s electronic structure and thus its use in semiconductor technologies. In the second part of this work, the formation of intrinsic ripples in graphene sheets under isotropic compression is examined. An isotropic compression of graphene is shown to induce a structural deformation on the basis of Density Functional Perturbation Theory. Static instabilities, indicated by imaginary fre- quency phonon modes, are induced in the high symmetry Γ – K (zigzag) and Γ – M (armchair) directions by an isotropic compressive strain of the graphene sheet. The wavelength of the unstable modes (ripples) is directly related to the magnitude of the strain and remarkably insensitive to the direction of propagation in the 2D lattice. These calculations further suggest that the formation energy of the ripple is isotropic for lower strains and becomes anisotropic for larger strains. This is a result of graphene’s elastic property, which is depen- dent on direction and strain. Within the quasi harmonic approximation this is combined with the observation that molecular adsorption energies depend strongly on curvature to suggest a strategy for generating ordered overlayers in order to tune the functional properties of graphene. Based on the results of this work, we can conclude that (pre-)rippled graphene sheets can be used to direct molecular adsorption in order to form specific patterns by tuning the thermodynamic equilibrium of the addition reaction of small (organic) molecules.Open Acces