Source author record

J. M. Pitarke

J. M. Pitarke 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

14works
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

14 published item(s)

preprint2022arXiv

Semilocal approximations to the Kohn-Sham exchange potential as applied to a metal surface

Several semilocal exchange potentials usually employed in the framework of density-functional theory (DFT) are tested and compared with their exact counterpart, the exchange Optimized Effective Potential (OEP), as applied to the jellium-slab model of a metal-vacuum interface. Driven by their explicit dependence on the ground-state density, its gradient, and its kinetic-energy density, the three analyzed semilocal exchange potentials approach their respective asymptotic limits faster than in the case of the OEP, all of them having an asymptotic scaling of the form $-α\,e^2/z + V_{\infty}$, with $α< 1$. Here we provide the leading analytic asymptotics of the three model potentials under study, and we find that none of them exhibits the exact OEP slab asymptotics $-\;e^2/z$. While the so-called Becke-Roussel potential's leading asymptote is close to its exact OEP counterpart, the other two model potentials under study approach a material-dependent positive constant value far into the vacuum, resulting in considerably overestimated ionization potentials.

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.

preprint2016arXiv

Semilocal density functional theory with correct surface asymptotics

Semilocal density functional theory is the most used computational method for electronic structure calculations in theoretical solid-state physics and quantum chemistry of large systems, providing good accuracy with a very attractive computational cost. Nevertheless, because of the non-locality of the exchange-correlation hole outside a metal surface, it was always considered inappropriate to describe the correct surface asymptotics. Here, we derive, within the semilocal density functional theory formalism, an exact condition for the image-like surface asymptotics of both the exchange-correlation energy per particle and potential. We show that this condition can be easily incorporated into a practical computational tool, at the simple meta-generalized-gradient approximation level of theory. Using this tool, we also show that the Airy-gas model exhibits asymptotic properties that are closely related to the ones at metal surfaces. This result highlights the relevance of the linear effective potential model to the metal surface asymptotics.

preprint2013arXiv

Acoustic plasmons in extrinsic free-standing graphene

An acoustic plasmon is predicted to occur, in addition to the conventional two-dimensional (2D) plasmon, as the collective motion of a system of two types of electronic carriers coexisting in the very same 2D band of extrinsic (doped or gated) graphene. The origin of this novel mode resides in the strong anisotropy that is present in the graphene band structure near the Dirac point. This anisotropy allows for the coexistence of carriers moving with two distinct Fermi velocities along the Gamma-K direction, which leads to two modes of collective oscillation: one mode in which the two types of electrons oscillate in phase with one another [this is the conventional 2D graphene plasmon, which at long wavelengths (q->0) has the same dispersion, q^1/2, as the conventional 2D plasmon of a 2D free electron gas], and the other mode found here corresponding to a low-frequency acoustic oscillation [whose energy exhibits at long wavelengths a linear dependence on the 2D wavenumber q] in which the two types of electrons oscillate out of phase. If this prediction is confirmed experimentally, it will represent the first realization of acoustic plasmons originated in the collective motion of a system of two types of carriers coexisting within the very same band.

preprint2012arXiv

Momentum-space finite-size corrections for Quantum-Monte-Carlo calculations

Extended solids are frequently simulated as finite systems with periodic boundary conditions, which due to the long-range nature of the Coulomb interaction may lead to slowly decaying finite- size errors. In the case of Quantum-Monte-Carlo simulations, which are based on real space, both real-space and momentum-space solutions to this problem exist. Here, we describe a hybrid method which using real-space data models the spherically averaged structure factor in momentum space. We show that (i) by integration our hybrid method exactly maps onto the real-space model periodic Coulomb-interaction (MPC) method and (ii) therefore our method combines the best of both worlds (real-space and momentum-space). One can use known momentum-resolved behavior to improve convergence where MPC fails (e.g., at surface-like systems). In contrast to pure momentum-space methods, our method only deals with a simple single-valued function and, hence, better lends itself to interpolation with exact small-momentum data as no directional information is needed. By virtue of integration, the resulting finite-size corrections can be written as an addition to MPC.

preprint2011arXiv

Adiabatic-Connection-Fluctuation-Dissipation approach to the long-range behavior of the exchange-correlation energy at metal surfaces: A numerical study for jellium slabs

A still open issue in many-body theory is the asymptotic behavior of the exchange-correlation energy and potential in the vacuum region of a metal surface. Here we report a numerical study of the position-dependent exchange-correlation energy for jellium slabs, as obtained by combining the formally exact adiabatic-connection-fluctuation-dissipation theorem with either time-dependent density-functional theory or an inhomogeneous Singwi-Tosi-Land-Sjölander approach. We find that the inclusion of correlation allows to obtain well-converged semi-infinite-jellium results (independent of the slab thickness) that exhibit an image-like asymptotic behavior close to the classical image potential $V_{im}(z)=-e^2/4z$.

preprint2010arXiv

An efficient method for the Quantum Monte Carlo evaluation of the static density-response function of a many-electron system

In a recent Letter we introduced Hellmann-Feynman operator sampling in diffusion Monte Carlo calculations. Here we derive, by evaluating the second derivative of the total energy, an efficient method for the calculation of the static density-response function of a many-electron system. Our analysis of the effect of the nodes suggests that correlation is described correctly and we find that the effect of the nodes can be dealt with.

preprint2010arXiv

Benchmark Quantum Monte Carlo calculations of the ground-state kinetic, interaction, and total energy of the three-dimensional electron gas

We report variational and diffusion Quantum Monte Carlo ground-state energies of the three-dimensional electron gas using a model periodic Coulomb interaction and backflow corrections for N=54, 102, 178, and 226 electrons. We remove finite-size effects by extrapolation and we find lower energies than previously reported. Using the Hellman-Feynman operator sampling method introduced in Phys. Rev. Lett. 99, 126406 (2007), we compute accurately, within the fixed-node pproximation, the separate kinetic and interaction contributions to the total ground-state energy. The difference between the interaction energies obtained from the original Slater-determinant nodes and the backflow-displaced nodes is found to be considerably larger than the difference between the corresponding kinetic energies.

preprint2010arXiv

Localized versus extended systems in density-functional theory: some lessons from the Kohn-Sham exact exchange potential

A long-standing puzzle in density-functional theory is the issue of the long-range behavior of the Kohn-Sham exchange-correlation potential at metal surfaces. As an important step towards its solution, it is proved here, through a rigurouos asymptotic analysis and accurate numerical solution of the Optimized-Effective-Potential integral equation, that the Kohn-Sham exact exchange potential decays as $\ln(z)/z$ far into the vacuum side of an {\it extended} semi-infinite jellium. In contrast to the situation in {\it localized} systems, like atoms, molecules, and slabs, this dominant contribution does not arise from the so-called Slater potential. This exact-exchange result provides a strong constraint on the suitability of approximate correlation-energy functionals.

preprint2009arXiv

Position-dependent exact-exchange energy for slabs and semi-infinite jellium

The position-dependent exact-exchange energy per particle $\varepsilon_x(z)$ (defined as the interaction between a given electron at $z$ and its exact-exchange hole) at metal surfaces is investigated, by using either jellium slabs or the semi-infinite (SI) jellium model. For jellium slabs, we prove analytically and numerically that in the vacuum region far away from the surface $\varepsilon_{x}^{\text{Slab}}(z \to \infty) \to - e^{2}/2z$, {\it independent} of the bulk electron density, which is exactly half the corresponding exact-exchange potential $V_{x}(z \to \infty) \to - e^2/z$ [Phys. Rev. Lett. {\bf 97}, 026802 (2006)] of density-functional theory, as occurs in the case of finite systems. The fitting of $\varepsilon_{x}^{\text{Slab}}(z)$ to a physically motivated image-like expression is feasible, but the resulting location of the image plane shows strong finite-size oscillations every time a slab discrete energy level becomes occupied. For a semi-infinite jellium, the asymptotic behavior of $\varepsilon_{x}^{\text{SI}}(z)$ is somehow different. As in the case of jellium slabs $\varepsilon_{x}^{\text{SI}}(z \to \infty)$ has an image-like behavior of the form $\propto - e^2/z$, but now with a density-dependent coefficient that in general differs from the slab universal coefficient 1/2. Our numerical estimates for this coefficient agree with two previous analytical estimates for the same. For an arbitrary finite thickness of a jellium slab, we find that the asymptotic limits of $\varepsilon_{x}^{\text{Slab}}(z)$ and $\varepsilon_{x}^{\text{SI}}(z)$ only coincide in the low-density limit ($r_s \to \infty$), where the density-dependent coefficient of the semi-infinite jellium approaches the slab {\it universal} coefficient 1/2.

preprint2008arXiv

Band structure effects on the Be(0001) acoustic-surface-plasmon energy dispersion

We report first-principles calculations of acoustic surface plasmons on the (0001) surface of Be, as obtained in the random-phase approximation of many-body theory. The energy dispersion of these collective excitations has been obtained along two symmetry directions. Our results show a considerable anisotropy of acoustic surface plasmons, and underline the capability of experimental measurements of these plasmons to {\it map} the electron-hole excitation spectrum of the quasi two-dimensional Shockley surface state band that is present on the Be(0001) surface.

preprint2006arXiv

Quantum Monte Carlo modelling of the spherically averaged structure factor of a many-electron system

The interaction and exchange-correlation contributions to the ground-state energy of an arbitrary many-electron system can be obtained from a spherical average of the wavevector-dependent diagonal structure factor (SF). We model the continuous-k spherically averaged SF using quantum Monte Carlo calculations in finite simulation cells. We thus derive a method that allows to substantially reduce the troublesome Coulomb finite-size errors that are usually present in ground-state energy calculations. To demonstrate this, we perform variational Monte Carlo calculations of the interaction energy of the homogeneous electron gas. The method is, however, equally applicable to arbitrary inhomogeneous systems.

preprint2004arXiv

Nonlinear, band-structure, and surface effects in the interaction of charged particles with solids

A survey is presented of various aspects of the interaction of charged particles with solids. In the framework of many-body perturbation theory, we study the nonlinear interaction of charged particles with a free gas of interacting electrons; in particular, nonlinear corrections to the stopping power are analyzed, and special emphasis is made on the separate contributions that are originated in the excitation of either electron-hole pairs, single plasmons, or double plasmons. Ab initio calculations of the electronic energy loss of ions moving in real solids are also presented, and the energy loss of charged particles interacting with simple metal surfaces is addressed.

preprint1997arXiv

Effective medium theory of the optical properties of aligned carbon nanotubes

We present an effective medium theory in order to analyze the reported optical properties of aligned carbon nanotube films (W.A. de Heer et al., Science 268, 845 (1995)). This methodology is based on photonic band structure calculations and allows treatment of complex media consisting of particles that interact strongly. We also develop a simple Maxwell-Garnett type approach for studying this system. In comparing the results of both mean field theories, we demonstrate that the inclusion of the full electromagnetic coupling between the nanotubes, as our numerical scheme does, is necessary for a complete explanation of the experimental data.