Source author record

Adrienn Ruzsinszky

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2021arXiv

First-principles wavevector- and frequency-dependent exchange-correlation kernel for jellium at all densities

We propose a spatially and temporally nonlocal exchange-correlation (xc) kernel for the spin-unpolarized fluid phase of ground-state jellium, for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact xc kernel, to accurately describe the correlation energies of bulk jellium, and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel [A. Ruzsinszky, et al., Phys. Rev. B 101, 245135 (2020)] are maintained in the new revised MCP07 (rMCP07) kernel, while others are added. The revision $f_\mathrm{xc}^\mathrm{rMCP07}(q,ω)$ differs from MCP07 only for non-zero frequencies $ω$. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wavevector $q$ tends to zero, the kernel has a $-4πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.}

preprint2020arXiv

Constraint-based Wavevector- and Frequency-dependent Exchange-Correlation Kernel of the Uniform Electron Gas

According to time-dependent density functional theory, the exact exchange-correlation kernel f$_{xc}$(n, q, $ω$) determines not only the ground-state energy but also the excited-state energies/lifetimes and time-dependent linear density response of an electron gas of uniform density n $=$ 3/(4$π$r$^3_s$). Here we propose a parametrization of this function based upon the satisfaction of exact constraints. For the static ($ω$ = 0) limit, we modify the model of Constantin and Pitarke at small wavevector q to recover the known second-order gradient expansion, plus other changes. For all frequencies $ω$ at q $=$ 0, we use the model of Gross, Kohn, and Iwamoto. A Cauchy integral extends this model to complex $ω$ and implies the standard Kramers-Kronig relations. A scaling relation permits closed forms for not only the imaginary but also the real part of f$_{xc}$ for real $ω$. We then combine these ingredients by damping out the $ω$ dependence at large q in the same way that the q dependence is damped. Away from q $=$ 0 and $ω$ $=$ 0, the correlation contribution to the kernel becomes dominant over exchange, even at r$_s$ $=$ 4, the valence electron density of metallic sodium. The resulting correlation energy from integration over imaginary $ω$ is essentially exact. The plasmon pole of the density response function is found by analytic continuation of f$_{xc}$ to $ω$ just below the real axis, and the resulting plasmon lifetime first decreases from infinity and then increases as q grows from 0 toward the electron-hole continuum. A static charge-density wave is found for r$_s$ $>$ 69, and shown to be associated with softening of the plasmon mode. The exchange-only version of our static kernel confirms Overhauser's 1968 prediction that correlation enhances the charge-density wave.

preprint2020arXiv

Examining the order-of-limits problem and lattice constant performance of the Tao--Mo Functional

In their recent communication [Phys. Rev. Lett., 117, 073001 (2016)] Tao and Mo presented a semi-local density functional derived from the density matrix expansion of the exchange hole localised by a general coordinate transformation. We show that the order-of-limits problem present in the functional, dismissed as harmless in the original publication, causes severe errors in predicted phase transition pressures. We also show that the claim that lattice volume prediction accuracy exceeds that of existing similar functionals was based on comparison to reference data that misses anharmonic zero-point expansion and consequently overestimates accuracy. By highlighting these omissions, we give a more accurate assessment of the Tao-Mo functional and show a simple route to resolving the problems.

preprint2020arXiv

The Fermi-Löwdin self-interaction correction for ionization energies of organic molecules

(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest-occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. The inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction corrections on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zσ, only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.

preprint2020arXiv

Understanding plasmon dispersion in nearly-free-electron metals: the relevance of exact constraints for novel exchange-correlation kernels within time-dependent density functional theory

Small-wavevector excitations in Coulomb-interacting systems can be decomposed into the high-energy collective longitudinal plasmon and the low-energy single-electron excitations. At the critical wavevector and corresponding frequency where the plasmon branch merges with the single-electron excitation region, the collective energy of the plasmon dissipates into single electron-hole excitations. The jellium model provides a reasonable description of the electron-energy-loss spectrum (EELS) of metals close to the free-electron limit. The random phase approximation (RPA) is exact in the high-density limit but can capture the plasmonic dispersion reasonably even for densities with rs > 1. RPA and all beyond-RPA methods investigated here, result in a wrong infinite plasmon lifetime for a wavevector smaller than the critical one where the plasmon dispersion curve runs into particle-hole excitations. Exchange-correlation kernel corrections to RPA modify the plasmon dispersion curve. There is however a large difference in the construction and form of the kernels investigated earlier. Our current work introduces recent model exchange-only and exchange-correlation kernels and discusses the relevance of some exact constraints in the construction of the kernel. We show that, because the plasmon dispersion samples a range of wavevectors smaller than the range sampled by the correlation energy, different kernels can make a strong difference for the correlation energy and a weak difference for the plasmon dispersion. This work completes our understanding about the plasmon dispersion in realistic metals, such as Cs, where a negative plasmon dispersion has been observed. We find only positive plasmon dispersion in jellium at the density for Cs.

preprint2019arXiv

A simple self-interaction correction to RPA-like correlation energies

The random phase approximation (RPA) is exact for the exchange energy of a many-electron ground state, but RPA makes the correlation energy too negative by about 0.5 eV/electron. That large short-range error, which tends to cancel out of iso-electronic energy differences, is largely corrected by an exchange-correlation kernel, or (as in RPA+) by an additive local or semilocal correction. RPA+ is by construction exact for the homogeneous electron gas, and it is also accurate for the jellium surface. RPA+ often gives realistic total energies for atoms or solids in which spin-polarization corrections are absent or small. RPA and RPA+ also yield realistic singlet binding energy curves for H2 and N2, and thus RPA+ yields correct total energies even for spin-unpolarized atoms with fractional spins and strong correlation, as in stretched H2 or N2. However, RPA and RPA+ can be very wrong for spin-polarized one-electron systems (especially for stretched H2+), and also for the spin-polarization energies of atoms. The spin-polarization energy is often a small part of the total energy of an atom, but important for ionization energies, electron affinities, and the atomization energies of molecules. Here we propose a computationally efficient generalized RPA+ (gRPA+) that changes RPA+ only for spin-polarized systems by making gRPA+ exact for all one-electron densities, in the same simple semilocal way that the correlation energy densities of many meta-generalized gradent approximations are made self-correlation free. By construction, gRPA+ does not degrade the exact RPA+ description of jellium. gRPA+ is found to greatly improve upon RPA and RPA+ for the ionization energies and electron affinities of light atoms. Many versions of RPA with an approximate exchange-correlation kernel fail to be exact for all one-electron densities, and they can also be self-interaction corrected in this way.

preprint2019arXiv

Molecule-surface interaction from van der Waals-corrected semilocal density functionals: the example of thiophene on transition-metal surfaces

Semi-local density functional approximations are widely used. None of them can capture the long-range van der Waals (vdW) attraction between separated subsystems, but they differ remarkably in the extent to which they capture intermediate-range vdW effects responsible for equilibrium bonds between neighboring small closed-shell subsystems. The local density approximation (LDA) often overestimates this effect, while the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) underestimates it. The strongly-constrained and appropriately normed (SCAN) meta-GGA often estimates it well. All of these semi-local functionals require an additive non-local correction such as the revised Vydrov-Van Voorhis 2010 (rVV10) to capture the long-range part. This work reports adsorption energies and the corresponding geometry of the aromatic thiophene (C$_4$H$_4$S) bound to transition metal surfaces. The adsorption process requires a genuine interplay of covalent and weak binding and requires a simultaneously accurate description of surface and adsorption energies with the correct prediction of the adsorption site. All these quantities must come from well balanced short and long-range correlation effects for a universally applicable method for weak interactions with chemical accuracy.

preprint2015arXiv

SCAN: An Efficient Density Functional Yielding Accurate Structures and Energies of Diversely-Bonded Materials

Kohn-Sham density functional theory (DFT) is a widely-used electronic structure theory for materials as well as molecules. DFT is needed especially for large systems, ab initio molecular dynamics, and high-throughput searches for functional materials. DFT's accuracy and computational efficiency are limited by the approximation to its exchange-correlation energy. Currently, the local density approximation (LDA) and generalized gradient approximations (GGAs) dominate materials computation mainly due to their efficiency. We show here that the recently developed non-empirical strongly constrained and appropriately normed (SCAN) meta-GGA improves significantly over LDA and the standard Perdew-Burke-Ernzerhof GGA for geometries and energies of diversely-bonded materials (including covalent, metallic, ionic, hydrogen, and van der Waals bonds) at comparable efficiency. Thus SCAN may be useful even for soft matter. Often SCAN matches or improves upon the accuracy of a computationally expensive hybrid functional, at almost-GGA cost. SCAN is therefore expected to have a broad impact on materials science.

preprint2015arXiv

Strongly Constrained and Appropriately Normed Semilocal Density Functional

The ground-state energy, electron density, and related properties of ordinary matter can be computed efficiently when the exchange-correlation energy as a functional of the density is approximated semilocally. We propose the first meta-GGA (meta-generalized gradient approximation) that is fully constrained, obeying all 17 known exact constraints that a meta-GGA can. It is also exact or nearly exact for a set of appropriate norms, including rare-gas atoms and nonbonded interactions. This SCAN (strongly constrained and appropriately normed) meta-GGA achieves remarkable accuracy for systems where the exact exchange-correlation hole is localized near its electron, and especially for lattice constants and weak interactions.

preprint2014arXiv

Gedanken Densities and Exact Constraints in Density Functional Theory

Approximations to the exact density functional for the exchange-correlation energy of a many-electron ground state can be constructed by satisfying constraints that are universal, i.e., valid for all electron densities. Gedanken densities are designed for the purpose of this construction, but need not be realistic. The uniform electron gas is an old gedanken density. Here, we propose a spherical two-electron gedanken density in which the dimensionless density gradient can be an arbitrary positive constant wherever the density is non-zero. The Lieb-Oxford lower bound on the exchange energy can be satisfied within a generalized gradient approximation (GGA) by bounding its enhancement factor or simplest GGA exchange-energy density. This enhancement-factor bound is well known to be sufficient, but our gedanken density shows that it is also necessary. The conventional exact exchange-energy density satisfies no such local bound, but energy densities are not unique, and the simplest GGA exchange-energy density is not an approximation to it. We further derive a strongly and optimally tightened bound on the exchange enhancement factor of a two-electron density, which is satisfied by the local density approximation but is violated by all published GGA's or meta-GGA's. Finally, some consequences of the non-uniform density-scaling behavior for the asymptotics of the exchange enhancement factor of a GGA or meta-GGA are given.

preprint2014arXiv

Self Interaction Corrected Density Functional Theory with Unitary Invariance: Applications to Molecules

Standard spin-density functionals for the exchange-correlation energy of a many-electron ground state make serious self-interaction errors which can be corrected by the Perdew-Zunger self-interaction correction (SIC). We propose a size-extensive construction of SIC orbitals which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional. The SIC orbitals are constructed from a unitary transformation that is explicitly dependent on the non-interacting one-particle density matrix. When this SIC is applied to the local spin-density approximation, improvements are found for the atomization energies of molecules.

preprint2013arXiv

Density Functionals that Recognize Covalent, Metallic, and Weak Bonds

Computationally-efficient semilocal approximations of density functional theory at the level of the local spin density approximation (LSDA) or generalized gradient approximation (GGA) poorly describe weak interactions. We show improved descriptions for weak bonds (without loss of accuracy for strong ones) from a newly-developed semilocal meta-GGA (MGGA), by applying it to molecules, surfaces, and solids. We argue that this improvement comes from using the right MGGA dimensionless ingredient to recognize all types of orbital overlap.

preprint2012arXiv

Effect of orbital-overlap dependence in density functionals

The semilocal meta generalized gradient approximation (MGGA) for the exchange-correlation functional of Kohn-Sham (KS) density functional theory can yield accurate ground-state energies simultaneously for atoms, molecules, surfaces, and solids, due to the inclusion of kinetic energy density as an input. We study for the first time the effect and importance of the dependence of MGGA on the kinetic energy density through the dimensionless inhomogeneity parameter, $α$, that characterizes the extent of orbital overlap. This leads to a simple and wholly new MGGA exchange functional, which interpolates between the single-orbital regime, where $α=0$, and the slowly varying density regime, where $α\approx 1$, and then extrapolates to $α\to \infty$. When combined with a variant of the Perdew-Burke-Erzerhof (PBE) GGA correlation, the resulting MGGA performs equally well for atoms, molecules, surfaces, and solids.

preprint2009arXiv

Assessing the Performance of Recent Density Functionals for Bulk Solids

We assess the performance of recent density functionals for the exchange-correlation energy of a nonmolecular solid, by applying accurate calculations with the GAUSSIAN, BAND, and VASP codes to a test set of 24 solid metals and non-metals. The functionals tested are the modified Perdew-Burke-Ernzerhof generalized gradient approximation (PBEsol GGA), the second-order GGA (SOGGA), and the Armiento-Mattsson 2005 (AM05) GGA. For completeness, we also test more-standard functionals: the local density approximation, the original PBE GGA, and the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA. We find that the recent density functionals for solids reach a high accuracy for bulk properties (lattice constant and bulk modulus). For the cohesive energy, PBE is better than PBEsol overall, as expected, but PBEsol is actually better for the alkali metals and alkali halides. For fair comparison of calculated and experimental results, we consider the zero-point phonon and finite-temperature effects ignored by many workers. We show how Gaussian basis sets and inaccurate experimental reference data may affect the rating of the quality of the functionals. The results show that PBEsol and AM05 perform somewhat differently from each other for alkali metal, alkaline earth metal and alkali halide crystals (where the maximum value of the reduced density gradient is about 2), but perform very similarly for most of the other solids (where it is often about 1). Our explanation for this is consistent with the importance of exchange-correlation nonlocality in regions of core-valence overlap.

preprint2009arXiv

Kinetic energy density functionals from the Airy gas, with an application to the atomization kinetic energies of molecules

We construct and study several semilocal density functional approximations for the positive Kohn-Sham kinetic energy density. These functionals fit the kinetic energy density of the Airy gas and they can be accurate for integrated kinetic energies of atoms, molecules, jellium clusters and jellium surfaces. We find that these functionals are the most accurate ones for atomization kinetic energies of molecules and for fragmentation of jellium clusters. We also report that local and semilocal kinetic energy functionals can show "binding" when the density of a spin unrestricted Kohn-Sham calculation is used.