loos/srDFT_GW
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

merge bib files

Browse Source
This commit is contained in:
Pierre-Francois Loos 2019-10-24 20:33:01 +02:00
parent aefa690a1e
commit e49b178c75
4 changed files with 329 additions and 12077 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

326
Manuscript/GW-srDFT.bib
View File

@ -1,13 +1,314 @@
%% This BibTeX bibliography file was created using BibDesk. %% This BibTeX bibliography file was created using BibDesk.
%% http://bibdesk.sourceforge.net/ %% http://bibdesk.sourceforge.net/
%% Created for Pierre-Francois Loos at 2019-10-23 10:10:27 +0200 %% Created for Pierre-Francois Loos at 2019-10-24 20:30:56 +0200
%% Saved with string encoding Unicode (UTF-8) %% Saved with string encoding Unicode (UTF-8)
@article{Dunning_1989,
Author = {Thom H. Dunning},
Date-Added = {2019-10-24 20:22:40 +0200},
Date-Modified = {2019-10-24 20:22:47 +0200},
Journal = {J. Chem. Phys.},
Pages = {1007},
Title = {Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen},
Volume = {90},
Year = {1989}}
@article{Tew_2007,
Author = {D. P. Tew and W. Klopper and C. Neiss and C. Hattig},
Date-Added = {2019-10-24 20:19:01 +0200},
Date-Modified = {2019-10-24 20:19:11 +0200},
Doi = {10.1039/b617230j},
Journal = {Phys. Chem. Chem. Phys.},
Pages = {1921},
Title = {Quintuple-{{$\zeta$}} Quality Coupled-Cluster Correlation Energies With Triple-{{$\zeta$}} Basis Sets},
Volume = {9},
Year = {2007},
Bdsk-Url-1 = {https://doi.org/10.1039/b617230j}}
@article{GinTewGarAla-JCTC-18,
Author = {Giner, Emmanuel and Tew, David P. and Garniron, Yann and Alavi, Ali},
Date-Added = {2019-10-24 20:19:01 +0200},
Date-Modified = {2019-10-24 20:19:01 +0200},
Doi = {10.1021/acs.jctc.8b00591},
Eprint = {https://doi.org/10.1021/acs.jctc.8b00591},
Journal = {Journal of Chemical Theory and Computation},
Note = {PMID: 30347156},
Number = {12},
Pages = {6240-6252},
Title = {Interplay between Electronic Correlation and Metal--Ligand Delocalization in the Spectroscopy of Transition Metal Compounds: Case Study on a Series of Planar Cu2+ Complexes},
Url = {https://doi.org/10.1021/acs.jctc.8b00591},
Volume = {14},
Year = {2018},
Bdsk-Url-1 = {https://doi.org/10.1021/acs.jctc.8b00591}}
@article{Hattig_2012,
Author = {C. Hattig and W. Klopper and A. Kohn and D. P. Tew},
Date-Added = {2019-10-24 20:19:01 +0200},
Date-Modified = {2019-10-24 20:19:29 +0200},
Doi = {10.1021/cr200168z},
Journal = {Chem. Rev.},
Pages = {4},
Title = {Explicitly Correlated Electrons in Molecules},
Volume = {112},
Year = {2012},
Bdsk-Url-1 = {https://doi.org/10.1021/cr200168z}}
@article{Kesharwani_2018,
Author = {M. K. Kesharwani and N. Sylvetsky and A. Kohn and D. P. Tew and Jan M. L. Martin},
Date-Added = {2019-10-24 20:19:01 +0200},
Date-Modified = {2019-10-24 20:19:22 +0200},
Doi = {10.1063/1.5048665},
Journal = {J. Chem. Phys.},
Pages = {154109},
Title = {Do CCSD And Approximate CCSD-F12 Variants Converge To The Same Basis Set Limits? The Case Of Atomization Energies And Approximate Ccsd-F12 Variants Converge To The Same Basis Set Limits? The Case Of Atomization Energies},
Volume = {149},
Year = {2018},
Bdsk-Url-1 = {https://doi.org/10.1063/1.5048665}}
@article{BooCleAlaTew-JCP-2012,
Author = {G. H. Booth and D. Cleland and A. Alavi and D. P. Tew},
Date-Added = {2019-10-24 20:19:01 +0200},
Date-Modified = {2019-10-24 20:19:01 +0200},
Doi = {10.1063/1.4762445},
Journal = {J. Chem. Phys.},
Pages = {164112},
Title = {An Explicitly Correlated Approach To Basis Set Incompleteness In Full Configuration Interaction Quantum Monte Carlo},
Volume = {137},
Year = {2012},
Bdsk-Url-1 = {https://doi.org/10.1063/1.4762445}}
@article{Tew08,
Author = {D. P. Tew},
Date-Added = {2019-10-24 20:18:51 +0200},
Date-Modified = {2019-10-24 20:18:51 +0200},
Journal = {J. Chem. Phys.},
Pages = {014104},
Title = {Second-order coalescence conditions of molecular wave functions},
Volume = {129},
Year = {2008}}
@article{Gruneis_2017,
Author = {A. Gr\"uneis and S. Hirata and Y.-Y. Ohnishi and S. Ten-no},
Date-Added = {2019-10-24 20:18:21 +0200},
Date-Modified = {2019-10-24 20:18:44 +0200},
Doi = {10.1063/1.4976974},
Journal = {J. Chem. Phys.},
Pages = {080901},
Title = {Perspective: Explicitly Correlated Electronic Structure Theory For Complex Systems},
Volume = {146},
Year = {2017},
Bdsk-Url-1 = {https://doi.org/10.1063/1.4976974}}
@article{Tenno_2012b,
Author = {S. Ten-no},
Date-Added = {2019-10-24 20:16:31 +0200},
Date-Modified = {2019-10-24 20:29:50 +0200},
Journal = {Theor. Chem. Acc.},
Pages = {1070},
Title = {Explicitly Correlated Wave Functions: Summary And Perspective},
Volume = {131},
Year = {2012}}
@article{Tenno_2012a,
Author = {S. Ten-no and J. Noga},
Date-Added = {2019-10-24 20:16:31 +0200},
Date-Modified = {2019-10-24 20:29:37 +0200},
Journal = {WIREs Comput. Mol. Sci.},
Pages = {114},
Title = {Explicitly Correlated Electronic Structure Theory From R12/F12 Ansatze},
Volume = {2},
Year = {2012}}
@article{Kong_2012,
Author = {L. Kong and F. A. Bischo and E. F. Valeev},
Date-Added = {2019-10-24 20:16:04 +0200},
Date-Modified = {2019-10-24 20:25:42 +0200},
Journal = {Chem. Rev.},
Pages = {75},
Title = {Explicitly Correlated R12/F12 Methods for Electronic Structure},
Volume = {112},
Year = {2012}}
@article{Noga_1994,
Author = {J. Noga and W. Kutzelnigg},
Date-Added = {2019-10-24 20:15:35 +0200},
Date-Modified = {2019-10-24 20:26:32 +0200},
Journal = {J. Chem. Phys.},
Pages = {7738},
Title = {Coupled Cluster Theory That Takes Care Of The Correlation Cusp By Inclusion Of Linear Terms In The Interelectronic Coordinates},
Volume = {101},
Year = {1994}}
@article{Hylleraas_1929,
Author = {E. A. Hylleraas},
Date-Added = {2019-10-24 20:14:55 +0200},
Date-Modified = {2019-10-24 20:18:09 +0200},
Journal = {Z. Phys.},
Pages = {347},
Title = {Neue Berechnung der Energie des Heliums im Grundzustande, sowie des tiefsten Terms von Ortho-Helium},
Volume = {54},
Year = {1929}}
@article{Kutzelnigg_1985,
Author = {W. Kutzelnigg},
Date-Added = {2019-10-24 20:14:19 +0200},
Date-Modified = {2019-10-24 20:25:59 +0200},
Journal = {Theor. Chim. Acta},
Pages = {445},
Title = {R12-Dependent Terms In The Wave Function As Closed Sums Of Partial Wave Amplitudes For Large L},
Volume = {68},
Year = {1985}}
@article{Kutzelnigg_1991,
Author = {W. Kutzelnigg and W. Klopper},
Date-Added = {2019-10-24 20:14:19 +0200},
Date-Modified = {2019-10-24 20:26:11 +0200},
Journal = {J. Chem. Phys.},
Pages = {1985},
Title = {Wave Functions With Terms Linear In The Interelectronic Coordinates To Take Care Of The Correlation Cusp. I. General Theory},
Volume = {94},
Year = {1991}}
@article{Kato_1957,
Author = {T. Kato},
Date-Added = {2019-10-24 20:13:25 +0200},
Date-Modified = {2019-10-24 20:25:17 +0200},
Journal = {Commun. Pure Appl. Math.},
Pages = {151},
Title = {On The Eigenfunctions Of Many-Particle Systems In Quantum Mechanics},
Volume = {10},
Year = {1957}}
@article{Gori-Giorgi_2006,
Author = {P. Gori-Giorgi and A. Savin},
Date-Added = {2019-10-24 20:11:03 +0200},
Date-Modified = {2019-10-24 20:24:22 +0200},
Journal = {Phys. Rev. A},
Pages = {032506},
Title = {Properties Of Short-Range And Long-Range Correlation Energy Density Functionals From Electron-Electron Coalescence},
Volume = {73},
Year = {2006}}
@article{Toulouse_2004,
Author = {J. Toulouse and F. Colonna and A. Savin},
Date-Added = {2019-10-24 20:10:53 +0200},
Date-Modified = {2019-10-24 20:19:46 +0200},
Journal = {Phys. Rev. A},
Keywords = {density functional theory; wave functions; exchange interactions (electron); electron correlations},
Number = {6},
Pages = {062505},
Publisher = {APS},
Title = {Long-range--short-range separation of the electron-electron interaction in density-functional theory},
Volume = {70},
Year = {2004}}
@article{Paziani_2006,
Author = {Simone Paziani and Saverio Moroni and Paola Gori-Giorgi and Giovanni B. Bachelet},
Date-Added = {2019-10-24 20:10:05 +0200},
Date-Modified = {2019-10-24 20:27:32 +0200},
Journal = {Phys. Rev. B},
Pages = {155111},
Title = {Local-Spin-Density Functional For Multideterminant Density Functional Theory},
Volume = {73},
Year = {2006}}
@article{Toulouse_2005,
Author = {J. Toulouse and P. Gori-Giorgi and A. Savin},
Date-Added = {2019-10-24 20:09:32 +0200},
Date-Modified = {2019-10-24 20:30:35 +0200},
Journal = {Theor. Chem. Acc.},
Pages = {305},
Title = {A Short-Range Correlation Energy Density Functional With Multi-Determinantal Reference},
Volume = {114},
Year = {2005}}
@book{Martin_2016,
Author = {Richard M. Martin and Lucia Reining and David M. Ceperley},
Date-Added = {2019-10-24 20:06:32 +0200},
Date-Modified = {2019-10-24 20:07:05 +0200},
Publisher = {Cambridge University Press},
Title = {Interacting Electrons: Theory and Computational Approaches},
Year = {2016}}
@article{Ferte_2019,
Author = {Anthony Fert\'e and Emmanuel Giner and Julien Toulouse},
Date-Added = {2019-10-24 20:05:01 +0200},
Date-Modified = {2019-10-24 20:12:16 +0200},
Doi = {10.1063/1.5082638},
Journal = {J. Chem. Phys.},
Pages = {084103},
Title = {Range-separated multideterminant density-functional theory with a short-range correlation functional of the on-top pair density},
Volume = {150},
Year = {2019},
Bdsk-Url-1 = {https://doi.org/10.1063/1.5082638}}
@article{Loos_2019,
Author = {P. F. Loos and B. Pradines and A. Scemama and J. Toulouse and E. Giner},
Date-Added = {2019-10-24 20:02:39 +0200},
Date-Modified = {2019-10-24 20:02:43 +0200},
Doi = {10.1021/acs.jpclett.9b01176},
Journal = {J. Phys. Chem. Lett.},
Pages = {2931--2937},
Title = {A Density-Based Basis-Set Correction for Wave Function Theory},
Volume = {10},
Year = {2019},
Bdsk-Url-1 = {https://doi.org/10.1021/acs.jpclett.9b01176}}
@article{Giner_2019,
Author = {E. Giner and A. Scemama and J. Toulouse and P. F. Loos},
Date-Added = {2019-10-24 20:02:05 +0200},
Date-Modified = {2019-10-24 20:02:10 +0200},
Journal = {J. Chem. Phys.},
Pages = {144118},
Title = {Chemically Accurate Excitation Energies With Small Basis Sets},
Volume = {151},
Year = {2019}}
@article{Lieb_1983,
Author = {E. H. Lieb},
Date-Added = {2019-10-24 20:00:53 +0200},
Date-Modified = {2019-10-24 20:01:02 +0200},
Journal = {Int. J. Quantum Chem.},
Pages = {243},
Title = {Density Functionals For Coulomb Systems},
Volume = {{24}},
Year = {1983}}
@article{Levy_1982,
Author = {M. Levy},
Date-Added = {2019-10-24 20:00:30 +0200},
Date-Modified = {2019-10-24 20:00:41 +0200},
Journal = {Phys. Rev. A},
Pages = {1200},
Title = {Electron Densities In Search Of Hamiltonians},
Volume = {26},
Year = {1982}}
@article{Levy_1979,
Author = {M. Levy},
Date-Added = {2019-10-24 20:00:01 +0200},
Date-Modified = {2019-10-24 20:00:07 +0200},
Journal = {Proc. Natl. Acad. Sci. U.S.A.},
Pages = {6062},
Title = {Universal Variational Functionals Of Electron Densities, First-Order Density Matrices, And Natural Spin-Orbitals And Solution Of The V-Representability Problem},
Volume = {76},
Year = {1979}}
@article{Giner_2018,
Author = {Emmanuel Giner and Barth\'elemy Pradines and Anthony Fert\'e and Roland Assaraf and Andreas Savin and Julien Toulouse},
Date-Added = {2019-10-24 19:59:35 +0200},
Date-Modified = {2019-10-24 19:59:43 +0200},
Journal = {J. Chem. Phys.},
Pages = {194301},
Title = {Curing basis-set convergence of wave-function theory using density-functional theory: A systematically improvable approach},
Volume = {149},
Year = {2018}}
@article{Duchemin_2017, @article{Duchemin_2017,
Author = {I. Duchemin and J. Li and X. Blase}, Author = {I. Duchemin and J. Li and X. Blase},
Date-Added = {2019-10-23 10:08:33 +0200}, Date-Added = {2019-10-23 10:08:33 +0200},
@ -17,7 +318,8 @@
Pages = {1199}, Pages = {1199},
Title = {Hybrid and Constrained Resolution-of-Identity Techniques for Coulomb Integrals}, Title = {Hybrid and Constrained Resolution-of-Identity Techniques for Coulomb Integrals},
Volume = {13}, Volume = {13},
Year = {2017}} Year = {2017},
Bdsk-Url-1 = {https://doi.org/10.1021/acs.jctc.6b01215}}
@article{Rojas_1995, @article{Rojas_1995,
Author = {H. N. Rojas and R. W. Godby and R. J. Needs}, Author = {H. N. Rojas and R. W. Godby and R. J. Needs},
@ -28,7 +330,8 @@
Pages = {1827}, Pages = {1827},
Title = {Space-Time Method for Ab Initio Calculations of Self-Energies and Dielectric Response Functions of Solids}, Title = {Space-Time Method for Ab Initio Calculations of Self-Energies and Dielectric Response Functions of Solids},
Volume = {74}, Volume = {74},
Year = {1995}} Year = {1995},
Bdsk-Url-1 = {https://doi.org/10.1103/PhysRevLett.74.1827}}
@article{Duchemin_2019, @article{Duchemin_2019,
Author = {Ivan Duchemin and Xavier Blase}, Author = {Ivan Duchemin and Xavier Blase},
@ -39,7 +342,8 @@
Pages = {174120}, Pages = {174120},
Title = {Separable Resolution-of-the-Identity with All-Electron Gaussian Bases: Application to Cubic-scaling RPA}, Title = {Separable Resolution-of-the-Identity with All-Electron Gaussian Bases: Application to Cubic-scaling RPA},
Volume = {150}, Volume = {150},
Year = {2019}} Year = {2019},
Bdsk-Url-1 = {https://doi.org/10.1063/1.5090605}}
@article{Weigend_2003a, @article{Weigend_2003a,
Author = {Weigend, Florian and Furche, Filipp and Ahlrichs, Reinhart}, Author = {Weigend, Florian and Furche, Filipp and Ahlrichs, Reinhart},
@ -758,20 +1062,6 @@
Year = {2012}, Year = {2012},
Bdsk-Url-1 = {https://dx.doi.org/10.1103/PhysRevB.85.035106}} Bdsk-Url-1 = {https://dx.doi.org/10.1103/PhysRevB.85.035106}}
@book{Martin_2016,
Abstract = {"Recent progress in the theory and computation of electronic structure is bringing an unprecedented level of capability for research. Many-body methods are becoming essential tools vital for quantitative calculations and understanding materials phenomena in physics, chemistry, materials science and other fields. This book provides a unified exposition of the most-used tools: many-body perturbation theory, dynamical mean field theory and quantum Monte Carlo simulations. Each topic is introduced with a less technical overview for a broad readership, followed by in-depth descriptions and mathematical formulation"--},
Address = {New York, NY},
Author = {Martin, Richard M. and Reining, Lucia and Ceperley, David},
Date-Added = {2018-03-06 19:42:15 +0000},
Date-Modified = {2018-03-06 19:42:15 +0000},
Isbn = {978-0-521-87150-1},
Keywords = {Electronic structure,Electrons,Many-body problem,Monte Carlo method,Perturbation (Quantum dynamics),Quantum theory},
Lccn = {QC176.8.E4 M368 2016},
Publisher = {{Cambridge University Press}},
Shorttitle = {Interacting Electrons},
Title = {Interacting Electrons: Theory and Computational Approaches},
Year = {2016}}
@article{Phillips_2014, @article{Phillips_2014,
Author = {Phillips, Jordan J. and Zgid, Dominika}, Author = {Phillips, Jordan J. and Zgid, Dominika},
Date-Added = {2018-03-06 19:42:15 +0000}, Date-Added = {2018-03-06 19:42:15 +0000},

48
Manuscript/GW-srDFT.tex
View File

@ -32,14 +32,6 @@
\newcommand{\PFL}[1]{\titou{(\underline{\bf PFL}: #1)}} \newcommand{\PFL}[1]{\titou{(\underline{\bf PFL}: #1)}}
\newcommand{\AS}[1]{\toto{(\underline{\bf AS}: #1)}} \newcommand{\AS}[1]{\toto{(\underline{\bf AS}: #1)}}
\usepackage{hyperref}
\hypersetup{
colorlinks=true,
linkcolor=blue,
filecolor=blue,
urlcolor=blue,
citecolor=blue
}
\newcommand{\mc}{\multicolumn} \newcommand{\mc}{\multicolumn}
\newcommand{\fnm}{\footnotemark} \newcommand{\fnm}{\footnotemark}
\newcommand{\fnt}{\footnotetext} \newcommand{\fnt}{\footnotetext}
@ -184,7 +176,7 @@ We also compute the ionization potentials of the five canonical nucleobases (ade
\section{Introduction} \section{Introduction}
\label{sec:intro} \label{sec:intro}
%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%
The purpose of many-body perturbation theory (MBPT) based on Green functions is to solve the formidable many-body problem by adding the electron-electron Coulomb interaction perturbatively starting from an independent-particle model. \cite{MarReiCep-BOOK-16} In this approach, the \textit{screening} of the Coulomb interaction is an essential quantity that is responsible for a rich variety of phenomena that would be otherwise absent (such as quasiparticle satellites and lifetimes). \cite{Aryasetiawan_1998, Onida_2002, Reining_2017} The purpose of many-body perturbation theory (MBPT) based on Green functions is to solve the formidable many-body problem by adding the electron-electron Coulomb interaction perturbatively starting from an independent-particle model. \cite{Martin_2016} In this approach, the \textit{screening} of the Coulomb interaction is an essential quantity that is responsible for a rich variety of phenomena that would be otherwise absent (such as quasiparticle satellites and lifetimes). \cite{Aryasetiawan_1998, Onida_2002, Reining_2017}
%\jt{Is it the screened Coulomb interaction which is responsible for these or more generally the Coulomb interaction? I would say the second.} %\jt{Is it the screened Coulomb interaction which is responsible for these or more generally the Coulomb interaction? I would say the second.}
The so-called {\GW} approximation is the workhorse of MBPT and has a long and successful history in the calculation of the electronic structure of solids. \cite{Aryasetiawan_1998, Onida_2002, Reining_2017} The so-called {\GW} approximation is the workhorse of MBPT and has a long and successful history in the calculation of the electronic structure of solids. \cite{Aryasetiawan_1998, Onida_2002, Reining_2017}
@ -224,17 +216,17 @@ Although obviously starting-point dependent, \cite{Bruneval_2013, Jacquemin_2016
For finite systems such as atoms and molecules, partially \cite{Hybertsen_1986, Shishkin_2007, Blase_2011, Faber_2011} or fully self-consistent \cite{Caruso_2012, Caruso_2013, Caruso_2013a, Caruso_2013b} {\GW} methods have shown great promise. \cite{Ke_2011, Blase_2011, Faber_2011, Koval_2014, Hung_2016, Blase_2018, Jacquemin_2017} For finite systems such as atoms and molecules, partially \cite{Hybertsen_1986, Shishkin_2007, Blase_2011, Faber_2011} or fully self-consistent \cite{Caruso_2012, Caruso_2013, Caruso_2013a, Caruso_2013b} {\GW} methods have shown great promise. \cite{Ke_2011, Blase_2011, Faber_2011, Koval_2014, Hung_2016, Blase_2018, Jacquemin_2017}
Similar to other electron correlation methods, MBPT methods suffer from the usual slow convergence of energetic properties with respect to the size of the one-electron basis set. Similar to other electron correlation methods, MBPT methods suffer from the usual slow convergence of energetic properties with respect to the size of the one-electron basis set.
This can be tracked down to the lack of explicit electron-electron terms modeling the infamous electron-electron coalescence point (also known as Kato cusp \cite{Kat-CPAM-57}) and, more specifically, the Coulomb correlation hole around it. This can be tracked down to the lack of explicit electron-electron terms modeling the infamous electron-electron coalescence point (also known as Kato cusp \cite{Kato_1957}) and, more specifically, the Coulomb correlation hole around it.
Pioneered by Hylleraas \cite{Hyl-ZP-29} in the 1930's and popularized in the 1990's by Kutzelnigg and coworkers \cite{Kut-TCA-85, NogKut-JCP-94, KutKlo-JCP-91} (and subsequently others \cite{KonBisVal-CR-12, HatKloKohTew-CR-12, TenNog-WIREs-12, GruHirOhnTen-JCP-17}), the so-called F12 methods overcome this slow convergence by employing geminal basis functions that closely resemble the correlation holes in electronic wave functions. Pioneered by Hylleraas \cite{Hylleraas_1929} in the 1930's and popularized in the 1990's by Kutzelnigg and coworkers \cite{Kutzelnigg_1985, Noga_1994, Kutzelnigg_1991} (and subsequently others \cite{Kong_2012, Hattig_2012, Tenno_2012a, Tenno_2012b, Gruneis_2017}), the so-called F12 methods overcome this slow convergence by employing geminal basis functions that closely resemble the correlation holes in electronic wave functions.
F12 methods are now routinely employed in computational chemistry and provide robust tools for electronic structure calculations where small basis sets may be used to obtain near complete basis set (CBS) limit accuracy. \cite{TewKloNeiHat-PCCP-07} F12 methods are now routinely employed in computational chemistry and provide robust tools for electronic structure calculations where small basis sets may be used to obtain near complete basis set (CBS) limit accuracy. \cite{Tew_2007}
The basis-set correction presented here follow a different route, and relies on the range-separated density-functional theory (RS-DFT) formalism to capture, thanks to a short-range correlation functional, the missing part of the short-range correlation effects. \cite{GinPraFerAssSavTou-JCP-18} The basis-set correction presented here follow a different route, and relies on the range-separated density-functional theory (RS-DFT) formalism to capture, thanks to a short-range correlation functional, the missing part of the short-range correlation effects. \cite{Giner_2018}
As shown in recent studies on both ground- and excited-state properties, \cite{LooPraSceTouGin-JPCL-19, GinSceTouLoo-JCP-19} similar to F12 methods, it significantly speeds up the convergence of energetics towards the CBS limit while avoiding the usage of the large auxiliary basis sets that are used in F12 methods to avoid the numerous three- and four-electron integrals. \cite{KonBisVal-CR-12, HatKloKohTew-CR-12, TenNog-WIREs-12, GruHirOhnTen-JCP-17, Barca_2018} As shown in recent studies on both ground- and excited-state properties, \cite{Loos_2019, Giner_2019} similar to F12 methods, it significantly speeds up the convergence of energetics towards the CBS limit while avoiding the usage of the large auxiliary basis sets that are used in F12 methods to avoid the numerous three- and four-electron integrals. \cite{Kong_2012, Hattig_2012, Tenno_2012a, Tenno_2012b, Gruneis_2017, Barca_2018}
Explicitly correlated F12 correction schemes have been derived for second-order Green function methods (GF2) \cite{SzaboBook, Casida_1989, Casida_1991, Stefanucci_2013, Ortiz_2013, Phillips_2014, Phillips_2015, Rusakov_2014, Rusakov_2016, Hirata_2015, Hirata_2017, Loos_2018} by Ten-no and coworkers \cite{Ohnishi_2016, Johnson_2018} and Valeev and coworkers. \cite{Pavosevic_2017, Teke_2019} Explicitly correlated F12 correction schemes have been derived for second-order Green function methods (GF2) \cite{SzaboBook, Casida_1989, Casida_1991, Stefanucci_2013, Ortiz_2013, Phillips_2014, Phillips_2015, Rusakov_2014, Rusakov_2016, Hirata_2015, Hirata_2017, Loos_2018} by Ten-no and coworkers \cite{Ohnishi_2016, Johnson_2018} and Valeev and coworkers. \cite{Pavosevic_2017, Teke_2019}
However, to the best of our knowledge, a F12-based correction for {\GW} has not been designed yet. However, to the best of our knowledge, a F12-based correction for {\GW} has not been designed yet.
In the present manuscript, we illustrate the performance of the density-based basis set correction developed in Refs.~\onlinecite{GinPraFerAssSavTou-JCP-18, LooPraSceTouGin-JPCL-19, GinSceTouLoo-JCP-19} on ionization potentials obtained within {\GOWO}. In the present manuscript, we illustrate the performance of the density-based basis set correction developed in Refs.~\onlinecite{Giner_2018, Loos_2019, Giner_2019} on ionization potentials obtained within {\GOWO}.
Note that the present basis set correction can be straightforwardly applied to other properties (\textit{e.g.}, electron affinities and fundamental gap), as well as other flavors of (self-consistent) {\GW} or Green function-based methods, such as GF2 (and its higher-order variants). Note that the present basis set correction can be straightforwardly applied to other properties (\textit{e.g.}, electron affinities and fundamental gap), as well as other flavors of (self-consistent) {\GW} or Green function-based methods, such as GF2 (and its higher-order variants).
The paper is organized as follows. The paper is organized as follows.
@ -252,7 +244,7 @@ Unless otherwise stated, atomic units are used throughout.
\subsection{MBPT with DFT basis set correction} \subsection{MBPT with DFT basis set correction}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Following Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}, we start by defining, for a $\Ne$-electron system with nuclei-electron potential $\vne(\br{})$, the approximate ground-state energy for one-electron densities $\n{}{}$ which are ``representable'' in a finite basis set $\Bas$ Following Ref.~\onlinecite{Giner_2018}, we start by defining, for a $\Ne$-electron system with nuclei-electron potential $\vne(\br{})$, the approximate ground-state energy for one-electron densities $\n{}{}$ which are ``representable'' in a finite basis set $\Bas$
\begin{equation} \begin{equation}
\E{0}{\Bas} = \min_{\n{}{} \in \cD^\Bas} \qty{ \F{}{}[n] + \int \vne(\br{}) \n{}{}(\br{}) \dbr{} }, \E{0}{\Bas} = \min_{\n{}{} \in \cD^\Bas} \qty{ \F{}{}[n] + \int \vne(\br{}) \n{}{}(\br{}) \dbr{} },
\label{eq:E0B} \label{eq:E0B}
@ -262,7 +254,7 @@ In this expression,
\begin{equation} \begin{equation}
\F{}{}[n] = \min_{\wf{}{} \rightsquigarrow \n{}{}} \mel*{\wf{}{}}{\hT + \hWee{}}{\wf{}{}} \F{}{}[n] = \min_{\wf{}{} \rightsquigarrow \n{}{}} \mel*{\wf{}{}}{\hT + \hWee{}}{\wf{}{}}
\end{equation} \end{equation}
is the exact Levy-Lieb universal density functional~\cite{Lev-PNAS-79, Lev-PRA-82, Lie-IJQC-83}, where the notation $\wf{}{} \rightsquigarrow \n{}{}$ in Eq.~\eqref{eq:E0B} states that $\wf{}{}$ yields the one-electron density $\n{}{}$. is the exact Levy-Lieb universal density functional, \cite{Levy_1979, Levy_1982, Lieb_1983} where the notation $\wf{}{} \rightsquigarrow \n{}{}$ in Eq.~\eqref{eq:E0B} states that $\wf{}{}$ yields the one-electron density $\n{}{}$.
$\hT$ and $\hWee{}$ are the kinetic and electron-electron interaction operators. $\hT$ and $\hWee{}$ are the kinetic and electron-electron interaction operators.
The exact Levy-Lieb universal density functional is then decomposed as The exact Levy-Lieb universal density functional is then decomposed as
\begin{equation} \begin{equation}
@ -273,13 +265,13 @@ where $\F{}{\Bas}[\n{}{}]$ is the Levy-Lieb density functional with wave functio
\begin{equation} \begin{equation}
\F{}{\Bas}[\n{}{}] = \min_{\wf{}{\Bas} \rightsquigarrow \n{}{}} \mel*{\wf{}{\Bas}}{ \hT + \hWee{}}{\wf{}{\Bas}}, \F{}{\Bas}[\n{}{}] = \min_{\wf{}{\Bas} \rightsquigarrow \n{}{}} \mel*{\wf{}{\Bas}}{ \hT + \hWee{}}{\wf{}{\Bas}},
\end{equation} \end{equation}
and $\bE{}{\Bas}[\n{}{}]$ is the complementary basis-correction density functional. \cite{GinPraFerAssSavTou-JCP-18} and $\bE{}{\Bas}[\n{}{}]$ is the complementary basis-correction density functional. \cite{Giner_2018}
In the present work, instead of using wave-function methods for calculating $\F{}{\Bas}[\n{}{}]$, we use Green-function methods. We assume that there exists a functional $\Omega^\Bas[\G{}{\Bas}]$ of $\Ne$-representable one-electron Green functions $\G{}{\Bas}(\br{},\br{}',\omega)$ representable in the basis set $\Bas$ and yielding the density $n$ which gives $\F{}{\Bas}[\n{}{}]$ at a stationary point In the present work, instead of using wave-function methods for calculating $\F{}{\Bas}[\n{}{}]$, we use Green-function methods. We assume that there exists a functional $\Omega^\Bas[\G{}{\Bas}]$ of $\Ne$-representable one-electron Green functions $\G{}{\Bas}(\br{},\br{}',\omega)$ representable in the basis set $\Bas$ and yielding the density $n$ which gives $\F{}{\Bas}[\n{}{}]$ at a stationary point
\begin{equation} \begin{equation}
\F{}{\Bas}[\n{}{}] = \stat{\G{}{\Bas} \rightsquigarrow \n{}{}} \Omega^\Bas[\G{}{\Bas}]. \F{}{\Bas}[\n{}{}] = \stat{\G{}{\Bas} \rightsquigarrow \n{}{}} \Omega^\Bas[\G{}{\Bas}].
\label{eq:FBn} \label{eq:FBn}
\end{equation} \end{equation}
The reason why we use a stationary condition rather than a minimization condition is that only a stationary property is generally known for functionals of the Green function. For example, we can choose for $\Omega^\Bas[G]$ a Klein-like energy functional (see, \textit{e.g.}, Refs.~\onlinecite{SteLee-BOOK-13, MarReiCep-BOOK-16, DahLee-JCP-05, DahLeeBar-IJQC-05, DahLeeBar-PRA-06}) The reason why we use a stationary condition rather than a minimization condition is that only a stationary property is generally known for functionals of the Green function. For example, we can choose for $\Omega^\Bas[G]$ a Klein-like energy functional (see, \textit{e.g.}, Refs.~\onlinecite{Stefanucci_2013, Martin_2016, Dahlen_2005, Dahlen_2005a, Dahlen_2006})
\begin{equation} \begin{equation}
\Omega^\Bas[\G{}{}] = \Tr[\ln( - \G{}{} ) ] - \Tr[ (\Gs{\Bas})^{-1} \G{}{} - 1 ] + \Phi_\Hxc^\Bas[\G{}{}], \Omega^\Bas[\G{}{}] = \Tr[\ln( - \G{}{} ) ] - \Tr[ (\Gs{\Bas})^{-1} \G{}{} - 1 ] + \Phi_\Hxc^\Bas[\G{}{}],
\label{eq:OmegaB} \label{eq:OmegaB}
@ -431,9 +423,9 @@ where $\eGOWO{\HOMO}$ and $\eGOWO{\LUMO}$ are the HOMO and LUMO orbital energies
%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%
The frequency-independent local self-energy $\bSig{}{\Bas}[\n{}{}](\br{},\br{}') = \bpot{}{\Bas}[\n{}{}](\br{}) \delta(\br{}-\br{}')$ originates from the functional derivative of complementary basis-correction density functionals $\bpot{}{\Bas}[\n{}{}](\br{}) = \delta \bE{}{\Bas}[\n{}{}] / \delta \n{}{}(\br{})$. The frequency-independent local self-energy $\bSig{}{\Bas}[\n{}{}](\br{},\br{}') = \bpot{}{\Bas}[\n{}{}](\br{}) \delta(\br{}-\br{}')$ originates from the functional derivative of complementary basis-correction density functionals $\bpot{}{\Bas}[\n{}{}](\br{}) = \delta \bE{}{\Bas}[\n{}{}] / \delta \n{}{}(\br{})$.
Here, we employ two types of complementary, short-range correlation functionals $\bE{}{\Bas}[\n{}{}]$: a short-range local-density approximation ($\srLDA$) functional with multideterminant reference \cite{TouGorSav-TCA-05,PazMorGorBac-PRB-06} and a short-range Perdew-Burke-Ernzerhof ($\srPBE$) correlation functional \cite{FerGinTou-JCP-19,LooPraSceTouGin-JPCL-19} which interpolates between the usual PBE functional \cite{PerBurErn-PRL-96} at $\mu = 0$ and the exact large-$\mu$ behavior~\cite{TouColSav-PRA-04, GorSav-PRA-06, PazMorGorBac-PRB-06} using the on-top pair density from the uniform-electron gas. \cite{LooPraSceTouGin-JPCL-19} Here, we employ two types of complementary, short-range correlation functionals $\bE{}{\Bas}[\n{}{}]$: a short-range local-density approximation ($\srLDA$) functional with multideterminant reference \cite{Toulouse_2005, Paziani_2006} and a short-range Perdew-Burke-Ernzerhof ($\srPBE$) correlation functional \cite{Ferte_2019, Loos_2019} which interpolates between the usual PBE functional \cite{Perdew_1996} at $\mu = 0$ and the exact large-$\mu$ behavior~\cite{Toulouse_2004, Gori-Giorgi_2006, Paziani_2006} using the on-top pair density from the uniform-electron gas. \cite{Loos_2019}
Additionally to the one-electron density, these RS-DFT functionals requires a range-separation function $\rsmu{}{\Bas}(\br{})$ which automatically adapts to the spatial non-homogeneity of the basis-set incompleteness error and is computed using the opposite-spin on-top pair density. Additionally to the one-electron density, these RS-DFT functionals requires a range-separation function $\rsmu{}{\Bas}(\br{})$ which automatically adapts to the spatial non-homogeneity of the basis-set incompleteness error and is computed using the one-electron density and opposite-spin on-top pair density originating from the HF or KS orbitals obtained in the basis set $\Bas$.
We refer the interested reader to Refs.~\onlinecite{GinPraFerAssSavTou-JCP-18, LooPraSceTouGin-JPCL-19, GinSceTouLoo-JCP-19} where our procedure is thoroughly detailed and the explicit expressions of these two short-range correlation functionals are provided. We refer the interested reader to Refs.~\onlinecite{Giner_2018, Loos_2019, Giner_2019} where our procedure is thoroughly detailed and the explicit expressions of these two short-range correlation functionals are provided.
The basis set corrected {\GOWO} quasiparticle energies are thus given by The basis set corrected {\GOWO} quasiparticle energies are thus given by
\begin{equation} \begin{equation}
@ -546,11 +538,11 @@ All the geometries have been extracted from the GW100 set. \cite{vanSetten_2015}
Unless otherwise stated, all the {\GOWO} calculations have been performed with the MOLGW software developed by Bruneval and coworkers. \cite{Bruneval_2016a} Unless otherwise stated, all the {\GOWO} calculations have been performed with the MOLGW software developed by Bruneval and coworkers. \cite{Bruneval_2016a}
The HF, PBE, and PBE0 calculations as well as the srLDA and srPBE basis set corrections have been computed with Quantum Package, \cite{QP2} which by default uses the SG-2 quadrature grid for the numerical integrations. The HF, PBE, and PBE0 calculations as well as the srLDA and srPBE basis set corrections have been computed with Quantum Package, \cite{QP2} which by default uses the SG-2 quadrature grid for the numerical integrations.
Frozen-core (FC) calculations are systematically performed. Frozen-core (FC) calculations are systematically performed.
The FC density-based basis set correction~\cite{LooPraSceTouGin-JPCL-19} is used consistently with the FC approximation in the {\GOWO} calculations. The FC density-based basis set correction~\cite{Loos_2019} is used consistently with the FC approximation in the {\GOWO} calculations.
The {\GOWO} quasiparticle energies have been obtained ``graphically'', \textit{i.e.}, by solving the non-linear, frequency-dependent quasiparticle equation \eqref{eq:QP-G0W0} (without linearization). The {\GOWO} quasiparticle energies have been obtained ``graphically'', \textit{i.e.}, by solving the non-linear, frequency-dependent quasiparticle equation \eqref{eq:QP-G0W0} (without linearization).
Moreover, the infinitesimal $\eta$ in Eq.~\eqref{eq:SigC} has been set to zero. Moreover, the infinitesimal $\eta$ in Eq.~\eqref{eq:SigC} has been set to zero.
Compared to the conventional $\order*{\Nbas^6}$ computational cost of {\GW} (where $\Nbas$ is the number of basis functions), the present basis set correction represents a marginal additional cost as further discussed in Refs.~\onlinecite{LooPraSceTouGin-JPCL-19, GinSceTouLoo-JCP-19}. Compared to the conventional $\order*{\Nbas^6}$ computational cost of {\GW} (where $\Nbas$ is the number of basis functions), the present basis set correction represents a marginal additional cost as further discussed in Refs.~\onlinecite{Loos_2019, Giner_2019}.
Note, however, that the formal $\order*{\Nbas^6}$ cost of {\GW} can be significantly reduced thanks to resolution-of-the-identity techniques \cite{vanSetten_2013, Bruneval_2016, Duchemin_2017} and other tricks. \cite{Rojas_1995, Duchemin_2019} Note, however, that the formal $\order*{\Nbas^6}$ cost of {\GW} can be significantly reduced thanks to resolution-of-the-identity techniques \cite{vanSetten_2013, Bruneval_2016, Duchemin_2017} and other tricks. \cite{Rojas_1995, Duchemin_2019}
%%% FIG 1 %%% %%% FIG 1 %%%
@ -559,7 +551,7 @@ Note, however, that the formal $\order*{\Nbas^6}$ cost of {\GW} can be significa
\hspace{1cm} \hspace{1cm}
\includegraphics[width=0.45\linewidth]{IP_G0W0PBE0_H2O} \includegraphics[width=0.45\linewidth]{IP_G0W0PBE0_H2O}
\caption{ \caption{
IP (in eV) of the water molecule computed at the {\GOWO} (black circles), {\GOWO}+srLDA (red squares), and {\GOWO}+srPBE (blue diamonds) levels of theory with increasingly large Dunning's basis sets \cite{Dun-JCP-89} (cc-pVDZ, cc-pVTZ, cc-pVQZ, and cc-pV5Z) with two different starting points: HF (left) and PBE0 (right). IP (in eV) of the water molecule computed at the {\GOWO} (black circles), {\GOWO}+srLDA (red squares), and {\GOWO}+srPBE (blue diamonds) levels of theory with increasingly large Dunning's basis sets \cite{Dunning_1989} (cc-pVDZ, cc-pVTZ, cc-pVQZ, and cc-pV5Z) with two different starting points: HF (left) and PBE0 (right).
The thick black line represents the CBS value obtained by extrapolation (see text for more details). The thick black line represents the CBS value obtained by extrapolation (see text for more details).
The green area corresponds to chemical accuracy (\textit{i.e.}, error below $1$ {\kcal} or $0.043$ eV). The green area corresponds to chemical accuracy (\textit{i.e.}, error below $1$ {\kcal} or $0.043$ eV).
\label{fig:IP_G0W0_H2O} \label{fig:IP_G0W0_H2O}
@ -604,7 +596,7 @@ Very similar conclusions are drawn at the {\GOWO}@{\PBEO} level (see Table \ref{
For example, at the {\GOWO}@PBE0+srLDA/cc-pVQZ level, the MAD is only $0.02$ eV with a maximum error as small as $0.09$ eV. For example, at the {\GOWO}@PBE0+srLDA/cc-pVQZ level, the MAD is only $0.02$ eV with a maximum error as small as $0.09$ eV.
It is worth pointing out that, for ground-state properties such as atomization and correlation energies, the density-based correction brought a more significant basis set reduction. It is worth pointing out that, for ground-state properties such as atomization and correlation energies, the density-based correction brought a more significant basis set reduction.
For example, we evidenced in Ref.~\onlinecite{LooPraSceTouGin-JPCL-19} that quintuple-$\zeta$ quality atomization and correlation energies are recovered with triple-$\zeta$ basis sets. For example, we evidenced in Ref.~\onlinecite{Loos_2019} that quintuple-$\zeta$ quality atomization and correlation energies are recovered with triple-$\zeta$ basis sets.
Here, the overall gain seems to be less important. Here, the overall gain seems to be less important.
The potential reasons for this could be: i) potential-based DFT corrections are usually less accurate than the ones based directly on energies, and ii) because the present scheme only corrects the basis set incompleteness error originating from the electron-electron cusp, some incompleteness remains at the HF or KS level. The potential reasons for this could be: i) potential-based DFT corrections are usually less accurate than the ones based directly on energies, and ii) because the present scheme only corrects the basis set incompleteness error originating from the electron-electron cusp, some incompleteness remains at the HF or KS level.
@ -703,7 +695,7 @@ This is quite remarkable as the number of basis functions jumps from $371$ to $7
\section{Conclusion} \section{Conclusion}
\label{sec:conclusion} \label{sec:conclusion}
%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%
In the present manuscript, we have shown that the density-based basis set correction developed by some of the authors in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18} and applied recently to ground- and excited-state properties \cite{LooPraSceTouGin-JPCL-19, GinSceTouLoo-JCP-19} can also be successfully applied to Green function methods such as {\GW}. In the present manuscript, we have shown that the density-based basis set correction developed by some of the authors in Ref.~\onlinecite{Giner_2018} and applied recently to ground- and excited-state properties \cite{Loos_2019, Giner_2019} can also be successfully applied to Green function methods such as {\GW}.
In particular, we have evidenced that the present basis set correction (which relies on LDA- or PBE-based short-range correlation functionals) significantly speeds up the convergence of IPs for small and larger molecules towards the CBS limit. In particular, we have evidenced that the present basis set correction (which relies on LDA- or PBE-based short-range correlation functionals) significantly speeds up the convergence of IPs for small and larger molecules towards the CBS limit.
These findings have been observed for different {\GW} starting points (HF, PBE, and PBE0). These findings have been observed for different {\GW} starting points (HF, PBE, and PBE0).
@ -727,6 +719,6 @@ This work has been supported through the EUR grant NanoX ANR-17-EURE-0009 in the
\end{acknowledgements} \end{acknowledgements}
%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%
\bibliography{GW-srDFT,GW-srDFT-control,biblio} \bibliography{GW-srDFT,GW-srDFT-control}
\end{document} \end{document}

2
Manuscript/SI/GW-srDFT-SI.tex
View File

@ -200,6 +200,6 @@
\end{figure*} \end{figure*}
\bibliography{../GW-srDFT,../GW-srDFT-control,../biblio} \bibliography{../GW-srDFT,../GW-srDFT-control}
\end{document} \end{document}

12030
Manuscript/biblio.bib
View File

File diff suppressed because it is too large Load Diff