BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Xinguo Ren

Xinguo Ren appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phcond-mat.str-elcond-mat.otherphysics.atm-clus

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2026arXiv

Influence of Exchange-Correlation Functionals and Neural Network Architectures on Li$^+$-Ion Conductivity in Solid-State Electrolyte from Molecular Dynamics Simulations with Machine-Learning Force Fields

With the rapid advancement of machine learning techniques for materials simulations, machine-learned force fields (MLFFs) have become a powerful tool that complements first-principles calculations by enabling high-accuracy molecular dynamics simulations over extended timescales. Typically, MLFFs are trained on data generated from density functional theory (DFT) using a specific exchange-correlation (XC) functional, with the goal of reproducing DFT-level properties. However, the uncertainties in MLFF-based simulations--arising from variations in both MLFF model architectures and the choice of XC functionals--remain insufficiently understood. In this work, we construct MLFF models of different architectures trained on DFT data from both semilocal and hybrid functionals to describe Li$^+$ diffusion in the solid-state electrolyte Li$_6$PS$_5$Cl. We systematically investigate how different XC functionals influence the Li$^+$ diffusion coefficient. To reduce statistical uncertainty, the mean squared displacements are averaged over 300 independent molecular dynamics (MD) trajectories of 70 ps each, yielding statistical variations below $1\%$. This enables a clear assessment of the respective influences of the functional and the MLFF model. Due to its tendency to underestimate band gaps and migration barriers, the semilocal functional predicts consistently higher Li$^+$ diffusion coefficients, compared to the hybrid functional. Furthermore, comparisons among various neural network methods reveal that the differences in predicted diffusion coefficients arising from different network architectures are of the same order of magnitude as those caused by different functionals, indicating that the choice of the network model itself substantially influences the MLFF predictions. This observation calls from an urgent need for standardized protocols to minimize model-dependent biases in MLFF-based MD.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Density functional theory plus dynamical mean field theory within the framework of linear combination of numerical atomic orbitals: Formulation and benchmarks

The combination of density functional theory with dynamical mean-field theory (DFT+DMFT) has become a powerful first-principles approach to tackle strongly correlated materials in condensed matter physics. The wide use of this approach relies on robust and easy-to-use implementations, and its implementation in various numerical frameworks will increase its applicability on the one hand and help crosscheck the validity of the obtained results on the other. In the work, we develop a formalism within the linear combination of numerical atomic orbital (NAO) basis set framework, which allows for merging NAO-based DFT codes with DMFT quantum impurity solvers. The formalism is implemented by interfacing two NAO-based DFT codes with three DMFT impurity solvers, and its validity is testified by benchmark calculations for a wide range of strongly correlated materials, including 3\textit{d} transition metal compounds, lanthanides, and actinides. Our work not only enables DFT+DMFT calculations using popular and rapidly developing NAO-based DFT code packages, but also facilitates the combination of more advanced beyond-DFT methodologies available in this codes with the DMFT machinery.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

DFT+$U$ within the framework of linear combination of numerical atomic orbitals

We present a formulation and implementation of the DFT+\textit{U} method within the framework of linear combination of numerical atomic orbitals (NAO). Our implementation not only enables single-point total energy and electronic-structure calculations but also provides access to atomic forces and stresses, hence allowing for full structure relaxations of periodic systems. Furthermore, our implementation allows one to deal with non-collinear spin texture, with the spin-orbit coupling (SOC) effect treated self-consistently. The key aspect behind our implementation is a suitable definition of the correlated subspace when multiple atomic orbitals with the same angular momentum are used, and this is addressed via the "Mulliken charge projector" constructed in terms of the first (most localized) atomic orbital within the $d/f$ angular momentum channel. The important Hubbard $U$ and Hund $J$ parameters can be estimated from a screened Coulomb potential of the Yukawa type, with the screening parameter either chosen semi-empirically or determined from the Thomas-Fermi screening model. Benchmark calculations are performed for four late transition metal monoxide bulk systems, i.e., MnO, FeO, CoO, and NiO, and for the 5$d$-electron compounds IrO$_2$. For the former type of systems, we check the performance of our DFT+$U$ implementation for calculating band gaps, magnetic moments, electronic band structures, as well as forces and stresses; for the latter, the efficacy of our DFT+$U$+SOC implementation is assessed. Systematic comparisons with available experimental results, and especially with the results from other implementation schemes are carried out, which demonstrate the validity of our NAO-based DFT+$U$ formalism and implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Importance of exact exchange to the geometric and electronic structures of Cs$_2$$B$$B'$$X_6$ double perovskites

We investigate the lead-free halide double perovskites (HDPs) Cs$ _2BB'X_6$ ($B$=Ag, Na; $B'$=In, Bi; $X$=Cl, Br) via first-principles calculations. We find that both the geometric and electric structures of the HDPs obtained by the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional are much better than those of the Perdew-Burke-Ernzerhof (PBE) functional. Importantly, we find that the electronic structures of DHPs are very sensitive to their geometries, especially the $B$-$X$ bond lengths. As a consequence, the electronic structures calculated by the HSE functional using the PBE optimized geometries may still significantly underestimate the band gaps, whereas the calculations on the HSE optimized geometries provide much more satisfactory results. The sensitivity of the band gaps of the DHPs to their geometries opens a promising path for the band structure engineering via doping and alloying. This work therefore provides an useful guideline for further improvement of HDPs materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ultrafast shift current dynamics in WS$_{2}$ monolayer

The shift current effect, in materials lacking inversion symmetry, may potentially allow the performance of photovoltaics to surpass the Shockley-Queisser limit for traditional p-n junction-based photovoltaics. Although the shift-current effect has been studied from first-principles via second-order perturbation theory, an understanding of the dynamics of hot carriers is still lacking. We investigate the dynamics of the shift current in monolayer WS$_{2}$ via real-time propagation time-dependent density functional theory (rt-TDDFT). We find that the shift current can be generated within 10 - 20 fs after turning on the lights and dissipates within approximately a few tens of femtoseconds after turning off the lights. This property can be used for ultrafast photon detection. This work provides an important step toward understanding the dynamics of shift-current effects, which is crucial for device applications.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

All-electron periodic $G_0W_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks

We present an all-electron, periodic {\GnWn} implementation within the numerical atomic orbital (NAO) basis framework. A localized variant of the resolution-of-the-identity (RI) approximation is employed to significantly reduce the computational cost of evaluating and storing the two-electron Coulomb repulsion integrals. We demonstrate that the error arising from localized RI approximation can be reduced to an insignificant level by enhancing the set of auxiliary basis functions, used to expand the products of two single-particle NAOs. An efficient algorithm is introduced to deal with the Coulomb singularity in the Brillouin zone sampling that is suitable for the NAO framework. We perform systematic convergence tests and identify a set of computational parameters, which can serve as the default choice for most practical purposes. Benchmark calculations are carried out for a set of prototypical semiconductors and insulators, and compared to independent reference values obtained from an independent $G_0W_0$ implementation based on linearized augmented plane waves (LAPW) plus high-energy localized orbitals (HLOs) basis set, as well as experimental results. With a moderate (FHI-aims \textit{tier} 2) NAO basis set, our $G_0W_0$ calculations produce band gaps that typically lie in between the standard LAPW and the LAPW+HLO results. Complementing \textit{tier} 2 with highly localized Slater-type orbitals (STOs), we find that the obtained band gaps show an overall convergence towards the LAPW+HLO results. The algorithms and techniques developed in this work pave the way for efficient implementations of correlated methods within the NAO framework.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Phase Stability of the Argon Crystal: A First-Principles Study Based on Random Phase Approximation plus Renormalized Single Excitation Corrections

The energy differences between the face-centered cubic (FCC) and hexagonal closed packed (HCP) structures of the argon (Ar) crystal are studied using the first-principles electronic-structure approach at the level of random phase approximation (RPA) plus renormalized single exitation (rSE) correction. By treating both structures at equal footing (i.e., employing the same computational supercell and k grid sampling), our RPA+rSE calculations show that, at zero temperature, the FCC structure is lower in energy than the HCP structure over a wide pressure range. The influence of zero-point energy (ZPE) is also studied and and it is found that ZPE only plays a secondary role in determining the relative stability of the two structures, whereas the electron correlation effect dominates. We further examine the equation of states in the high pressure regime, and our RPA+rSE results, complemented with phonon contributions, show excellent agreement with available experimental data. Finally, by computing the Gibbs free energies for both the FCC and HCP at different temperatures, we are able to generate a $T$-$P$ phase diagram for the Ar crystal, disclosing the pressure-temperature range for each phase. Our calculations show that the FCC phase has a slightly larger entropy and volume than HCP phase at the temperature and pressure condition.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ab Initio Bethe-Salpeter Equation Approach to Neutral Excitations in Molecules with Numeric Atom-Centered Orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as "Thiel's set". Literature reference data based on Gaussian-type orbitals are reproduced to about one meV precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted "tier2+aug2" is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large-scale ab initio simulations based on systematically improvable atomic basis

We present a first-principles computer code package (ABACUS) that is based on density functional theory and numerical atomic basis sets. Theoretical foundations and numerical techniques used in the code are described, with focus on the accuracy and transferability of the hierarchical atomic basis sets as generated using a scheme proposed by Chen, Guo and He [J. Phys.:Condens. Matter \textbf{22}, 445501 (2010)]. Benchmark results are presented for a variety of systems include molecules, solids, surfaces, and defects. All results show that the ABACUS package with its associated atomic basis sets is an efficient and reliable tool for simulating both small and large-scale materials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-Consistent Green Function Embedding for Advanced Electronic Structure Methods Based on a Dynamical Mean-Field Concept

We present an embedding scheme for periodic systems that facilitates the treatment of the physically important part (here the unit cell) with advanced electronic-structure methods, that are computationally too expensive for periodic systems. The rest of the periodic system is treated with computationally less demanding approaches, e.g., Kohn-Sham density-functional theory, in a self- consistent manner. Our scheme is based on the concept of dynamical mean-field theory (DMFT) formulated in terms of Green functions. In contrast to the original DMFT formulation for correlated model Hamiltonians, we here consider the unit cell as local embedded cluster in a first-principles way, that includes all electronic degrees of freedom. Our real-space dynamical mean-field embedding (RDMFE) scheme features two nested Dyson equations, one for the embedded cluster and another for the periodic surrounding. The total energy is computed from the resulting Green functions. The performance of our scheme is demonstrated by treating the embedded region with hybrid functionals and many-body perturbation theory in the GW approach for simple bulk systems. The total energy and the density of states converge rapidly with respect to the computational parameters and approach their bulk limit with increasing cluster (i.e., unit cell) size.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Static correlation and electron localization in molecular dimers from the self-consistent RPA and GW approximation

We investigate static correlation and delocalization errors in the self-consistent GW and random-phase approximation (RPA) by studying molecular dissociation of the H_2 and LiH molecules. Although both approximations contain topologically identical diagrams, the non-locality and frequency dependence of the GW self-energy crucially influence the different energy contributions to the total energy as compared to the use of a static local potential in the RPA. The latter leads to significantly larger correlation energies which allow for a better description of static correlation at intermediate bond distances. The substantial error found in GW is further analyzed by comparing spin-restricted and spin-unrestricted calculations. At large but finite nuclear separation their difference gives an estimate of the so-called fractional spin error normally determined only in the dissociation limit. Furthermore, a calculation of the dipole moment of the LiH molecule at dissociation reveals a large delocalization error in GW making the fractional charge error comparable to the RPA. The analyses are supplemented by explicit formulae for the GW Green's function and total energy of a simplified two-level model providing additional insights into the dissociation limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strong charge and spin fluctuations in La$_2$O$_3$Fe$_2$Se$_2$

The electronic structure and magnetic properties of the strongly correlated material La$_2$O$_3$Fe$_2$Se$_2$ are studied by using both the density function theory plus $U$ (DFT+$U$) method and the DFT plus Gutzwiller (DFT+G) variational method. The ground-state magnetic structure of this material obtained with DFT+$U$ is consistent with recent experiments, but its band gap is significantly overestimated by DFT+$U$, even with a small Hubbard $U$ value. In contrast, the DFT+G method yields a band gap of 0.1 - 0.2 eV, in excellent agreement with experiment. Detailed analysis shows that the electronic and magnetic properties of of La$_2$O$_3$Fe$_2$Se$_2$ are strongly affected by charge and spin fluctuations which are missing in the DFT+$U$ method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The origin of hyperferroelectricity in Li$B$O$_3$ ($B$=V, Nb, Ta, Os)

The electronic and structural properties of Li$B$O$_3$ ($B$=V, Nb, Ta, Os) are investigated via first-principles methods. We show that Li$B$O$_3$ belong to the recently proposed hyperferroelectrics, i.e., they all have unstable longitudinal optic phonon modes. Especially, the ferroelectric-like instability in the metal LiOsO$_3$, whose optical dielectric constant goes to infinity, is a limiting case of hyperferroelectrics. Via an effective Hamiltonian, we further show that, in contrast to normal proper ferroelectricity, in which the ferroelectric instability usually comes from long-range coulomb interactions, the hyperferroelectric instability is due to the structure instability driven by short-range interactions. This could happen in systems with large ion size mismatches, which therefore provides a useful guidance in searching for novel hyperferroelectrics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First-Principles Description of Charge Transfer in Donor-Acceptor Compounds from Self-Consistent Many-Body Perturbation Theory

We investigate charge transfer in prototypical molecular donor-acceptor compounds using hybrid density functional theory (DFT) and the GW approximation at the perturbative level (G0W0) and at full self-consistency (sc-GW). For the systems considered here, no charge transfer should be expected at large intermolecular separation according to photoemission experiment and accurate quantum-chemistry calculations. The capability of hybrid exchange-correlation functionals of reproducing this feature depends critically on the fraction of exact exchange $α$, as for small values of $α$ spurious fractional charge transfer is observed between the donor and the acceptor. G0W0 based on hybrid DFT yields the correct alignment of the frontier orbitals for all values of $α$. However, G0W0 has no capacity to alter the ground-state properties of the system, because of its perturbative nature. The electron density in donor-acceptor compounds thus remains incorrect for small $α$ values. In sc-GW, where the Green's function is obtained from the iterative solution of the Dyson equation, the electron density is updated and reflects the correct description of the level alignment at the GW level, demonstrating the importance of self-consistent many-body approaches for the description of ground- and excited-state properties in donor-acceptor systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bond Breaking and Bond Formation: How Electron Correlation is Captured in Many-Body Perturbation Theory and Density-Functional Theory

For the paradigmatic case of H2-dissociation we compare state-of-the-art many-body perturbation theory (MBPT) in the GW approximation and density-functional theory (DFT) in the exact-exchange plus random-phase approximation for the correlation energy (EX+cRPA). For an unbiased comparison and to prevent spurious starting point effects both approaches are iterated to full self-consistency (i.e. sc-RPA and sc-GW). The exchange-correlation diagrams in both approaches are topologically identical, but in sc-RPA they are evaluated with non-interacting and in sc-GW with interacting Green functions. This has a profound consequence for the dissociation region, where sc-RPA is superior to sc-GW. We argue that for a given diagrammatic expansion, the DFT framework outperforms the many-body framework when it comes to bond-breaking. We attribute this to the difference in the correlation energy rather than the treatment of the kinetic energy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schrödinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self-consistent GW: All-electron implementation with localized basis functions

This paper describes an all-electron implementation of the self-consistent GW (sc-GW) approach -- i.e. based on the solution of the Dyson equation -- in an all-electron numeric atom-centered orbital (NAO) basis set. We cast Hedin's equations into a matrix form that is suitable for numerical calculations by means of i) the resolution of identity technique to handle 4-center integrals; and ii) a basis representation for the imaginary-frequency dependence of dynamical operators. In contrast to perturbative G0W0, sc-GW provides a consistent framework for ground- and excited-state properties and facilitates an unbiased assessment of the GW approximation. For excited-states, we benchmark sc-GW for five molecules relevant for organic photovoltaic applications: thiophene, benzothiazole, 1,2,5-thiadiazole, naphthalene, and tetrathiafulvalene. At self-consistency, the quasi-particle energies are found to be in good agreement with experiment and, on average, more accurate than G0W0 based on Hartree-Fock (HF) or density-functional theory with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. Based on the Galitskii-Migdal total energy, structural properties are investigated for a set of diatomic molecules. For binding energies, bond lengths, and vibrational frequencies sc-GW and G0W0 achieve a comparable performance, which is, however, not as good as that of exact-exchange plus correlation in the random-phase approximation (EX+cRPA) and its advancement to renormalized second-order perturbation theory (rPT2). Finally, the improved description of dipole moments for a small set of diatomic molecules demonstrates the quality of the sc-GW ground state density.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Benchmark of GW Methods for Azabenzenes: Is the GW Approximation Good Enough?

Many-body perturbation theory in the GW approximation is a useful method for describing electronic properties associated with charged excitations. A hierarchy of GW methods exists, starting from non-self-consistent G0W0, through partial self-consistency in the eigenvalues (ev-scGW) and in the Green function (scGW0), to fully self-consistent GW (scGW). Here, we assess the performance of these methods for benzene, pyridine, and the diazines. The quasiparticle spectra are compared to photoemission spectroscopy (PES) experiments with respect to all measured particle removal energies and the ordering of the frontier orbitals. We find that the accuracy of the calculated spectra does not match the expectations based on their level of self-consistency. In particular, for certain starting points G0W0 and scGW0 provide spectra in better agreement with the PES than scGW.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Assessment of correlation energies based on the random-phase approximation

The random-phase approximation to the ground state correlation energy (RPA) in combination with exact exchange (EX) has brought Kohn-Sham (KS) density functional theory one step closer towards a universal, "general purpose first principles method". In an effort to systematically assess the influence of several correlation energy contributions beyond RPA, this work presents dissociation energies of small molecules and solids, activation energies for hydrogen transfer and non-hydrogen transfer reactions, as well as reaction energies for a number of common test sets. We benchmark EX+RPA and several flavors of energy functionals going beyond it: second-order screened exchange (SOSEX), single excitation (SE) corrections, renormalized single excitation (rSE) corrections, as well as their combinations. Both the single excitation correction as well as the SOSEX contribution to the correlation energy significantly improve upon the notorious tendency of EX+RPA to underbind. Surprisingly, activation energies obtained using EX+RPA based on a KS reference alone are remarkably accurate. RPA+SOSEX+rSE provides an equal level of accuracy for reaction as well as activation energies and overall gives the most balanced performance, which makes it applicable to a wide range of systems and chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Density-functional Theory for f electron Systems: the α-γ Phase Transition in Cerium

The isostructural α-γ phase transition in cerium is analyzed using density-functional theory with different exchange-correlation functionals, in particular the PBE0 hybrid functional and the exact- exchange plus correlation in the random-phase approximation [(EX+cRPA)@PBE0] approach. We show that the Hartree-Fock exchange part of the hybrid functional actuates two distinct solutions at zero temperature that can be associated with the α and γ phases of cerium. However, despite the relatively good structural and magnetic properties, PBE0 predicts the γ phase to be the stable phase at ambient pressure and zero temperature, in contradiction with low temperature experiments. EX+cRPA reverses the energetic ordering, which emphasizes the importance of correlation for rare- earth systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions

Efficient implementations of electronic structure methods are essential for first-principles modeling of molecules and solids. We here present a particularly efficient common framework for methods beyond semilocal density-functional theory, including Hartree-Fock (HF), hybrid density functionals, random-phase approximation (RPA), second-order Møller-Plesset perturbation theory (MP2), and the $GW$ method. This computational framework allows us to use compact and accurate numeric atom-centered orbitals (popular in many implementations of semilocal density-functional theory) as basis functions. The essence of our framework is to employ the "resolution of identity (RI)" technique to facilitate the treatment of both the two-electron Coulomb repulsion integrals (required in all these approaches) as well as the linear density-response function (required for RPA and $GW$). This is possible because these quantities can be expressed in terms of products of single-particle basis functions, which can in turn be expanded in a set of auxiliary basis functions (ABFs). The construction of ABFs lies at the heart of the RI technique, and here we propose a simple prescription for constructing the ABFs which can be applied regardless of whether the underlying radial functions have a specific analytical shape (e.g., Gaussian) or are numerically tabulated. We demonstrate the accuracy of our RI implementation for Gaussian and NAO basis functions, as well as the convergence behavior of our NAO basis sets for the above-mentioned methods. Benchmark results are presented for the ionization energies of 50 selected atoms and molecules from the G2 ion test set as obtained with $GW$ and MP2 self-energy methods, and the G2-I atomization energies as well as the S22 molecular interaction energies as obtained with the RPA method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unified description of ground and excited states of finite systems: the self-consistent GW approach

GW calculations with fully self-consistent G and W -- based on the iterative solution of the Dyson equation -- provide an approach for consistently describing ground and excited states on the same quantum mechanical level. We show that for the systems considered here self-consistent GW reaches the same final Green function regardless of the initial reference state. Self-consistency systematically improves ionization energies and total energies of closed shell systems compared to G_0W_0 based on Hartree-Fock and (semi)local density-functional theory. These improvements also translate to the electron density as exemplified by an improved description of dipole moments and permit us to assess the quality of ground state properties such as bond lengths and vibrational frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Beyond the Random Phase Approximation for the Electron Correlation Energy: The Importance of Single Excitations

The random phase approximation (RPA) for the electron correlation energy, combined with the exact-exchange energy, represents the state-of-the-art exchange-correlation functional within density-functional theory (DFT). However, the standard RPA practice -- evaluating both the exact-exchange and the RPA correlation energy using local or semilocal Kohn-Sham (KS) orbitals -- leads to a systematic underbinding of molecules and solids. Here we demonstrate that this behavior is largely corrected by adding a "single excitation" (SE) contribution, so far not included in the standard RPA scheme. A similar improvement can also be achieved by replacing the non-self-consistent exact-exchange total energy by the corresponding self-consistent Hartree-Fock total energy, while retaining the RPA correlation energy evaluated using Kohn-Sham orbitals. Both schemes achieve chemical accuracy for a standard benchmark set of non-covalent intermolecular interactions.

0 citations0 reviewsNo code linkedOpen access