Researcher profile

L. A. Constantin

L. A. Constantin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2016arXiv

Kinetic and Exchange Energy Densities near the Nucleus

We investigate the behavior of the kinetic and the exchange energy densities near the nuclear cusp of atomic systems. Considering hydrogenic orbitals, we derive analytical expressions near the nucleus, for single shells, as well as in the semiclassical limit of large non-relativistic neutral atoms. We show that a model based on the helium iso-electronic series is very accurate, as also confirmed by numerical calculations on real atoms up to two thousands electrons. Based on this model, we propose non-local density-dependent ingredients that are suitable for the description of the kinetic and exchange energy densities in the region close to the nucleus. These non-local ingredients are invariant under the uniform scaling of the density, and they can be used in the construction of non-local exchange-correlation and kinetic functionals.

preprint2016arXiv

Semiclassical atom theory applied to solid-state physics

Using the semiclassical neutral atom theory, we extend to fourth order the modified gradient expansion of the exchange energy of density functional theory. This expansion can be applied both to large atoms and solid-state problems. Moreover, we show that it can be employed to construct a simple and non-empirical generalized gradient approximation (GGA) exchange-correlation functional competitive with state-of-the-art GGAs for solids, but also reasonably accurate for large atoms and ordinary chemistry.

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.

preprint2015arXiv

Assessment of the TCA functional in computational chemistry and solid-state physics

We assess the Tognetti-Cortona-Adamo (TCA) generalized gradient approximation correlation functional [J. Chem. Phys. 128:034101 (2008)] for a variety of electronic systems. We find that, even if the TCA functional is not exact for the uniform electron gas, it is very accurate for the jellium surface correlation energies and it gives a realistic description of the quantum oscillations and surface effects of various jellium clusters, that are important model systems in computational chemistry and solid-state physics. When the TCA correlation is combined with the non-empirical PBEint, Wu-Cohen, and PBEsol$_b$ exchange functionals, the resulting exchange-correlation approximations provide good performances for a broad palette of systems and properties, being reasonably accurate for thermochemistry and geometry of molecules, transition metal complexes, non-covalent interactions,equilibrium lattice constants, bulk moduli, and cohesive energies of solids.

preprint2015arXiv

Global hybrids from the semiclassical atom theory satisfying the local density linear response

We propose global hybrid approximations of the exchange-correlation (XC) energy functional which reproduce well the modified fourth-order gradient expansion of the exchange energy in the semiclassical limit of many-electron neutral atoms and recover the full local density approximation (LDA) linear response. These XC functionals represent the hybrid versions of the APBE functional [Phys. Rev. Lett. 106, 186406, (2011)] yet employing an additional correlation functional which uses the localization concept of the correlation energy density to improve the compatibility with the Hartree-Fock exchange as well as the coupling-constant-resolved XC potential energy. Broad energetical and structural testings, including thermochemistry and geometry, transition metal complexes, non-covalent interactions, gold clusters and small gold-molecule interfaces, as well as an analysis of the hybrid parameters, show that our construction is quite robust. In particular, our testing shows that the resulting hybrid, including 20\% of Hartree-Fock exchange and named hAPBE, performs remarkably well for a broad palette of systems and properties, being generally better than popular hybrids (PBE0 and B3LYP). Semi-empirical dispersion corrections are also provided.

preprint2015arXiv

Gradient-dependent upper bound for the exchange-correlation energy and application to density functional theory

We propose a simple gradient-dependent bound for the exchange-correlation energy (sLL), based on the recent non-local bound derived by Lewin and Lieb. We show that sLL is equivalent to the original Lieb-Oxford bound in rapidly-varying density cases but it is tighter for slowly-varying density systems. To show the utility of the sLL bound we apply it to the construction of simple semilocal and non-local exchange and correlation functionals. In both cases improved results, with respect to the use of Lieb-Oxford bound, are obtained showing the power of the sLL bound.

preprint2015arXiv

Subsystem density functional theory with meta generalized gradient approximation exchange-correlation functionals

We analyze the methodology and the performance of subsystem density functional theory (DFT) with meta-generalized gradient approximation (meta-GGA) exchange-correlation functionals for non-bonded systems. Meta-GGA functionals depend on the Kohn-Sham kinetic energy density (KED), which is not known as an explicit functional of the density. Therefore, they cannot be directly applied in subsystem DFT calculations. We propose a Laplacian-level approximation to the KED which overcomes the problem and provides a simple and accurate way to apply meta-GGA exchange-correlation functionals in subsystem DFT calculations. The so obtained density and energy errors, with respect to the corresponding supermolecular calculations, are comparable with conventional approaches, depending almost exclusively on the approximations in the non-additive kinetic embedding term. An embedding energy error decomposition explains the accuracy of our method.

preprint2014arXiv

Generalized gradient approximation correlation energy functionals based on the uniform electron gas with gap model

We study the uniform electron gas with a gap model in the context of density functional theory. Based on this analysis, we construct two local gap models that realize generalized gradient approximation (GGA) correlation functionals satisfying numerous exact constraints for the correlation energy. The first one, named GAPc, fulfills the full second-order correlation gradient expansion at any density regime, is very accurate for jellium surfaces, comparable to state-of-the-art GGAs for atomic systems and molecular systems, and well compatible with known semilocal exchange. The second functional, named GAPloc, is satisfying the same exact conditions, except that the second-order gradient expansion is sacrificed for a better behavior under the Thomas-Fermi scaling and a more realistic correlation energy density of the hellium atom. The GAPloc functional displays a high accuracy for atomic correlation energies, still preserving a reasonable behavior for jellium surfaces. Moreover, it shows a higher compatibility with the Hartree-Fock exchange than other semilocal correlation functionals. This feature is explained in terms of the real-space analysis of the GAPloc correlation energy.

preprint2014arXiv

Kohn-Sham Kinetic Energy Density in the Nuclear and Asymptotic Regions: Deviations from the Von Weizsäcker Behavior and Applications to Density Functionals

We show that the Kohn-Sham positive-definite kinetic energy (KE) density significantly differs from the von Weizsäcker (VW) one at the nuclear cusp as well as in the asymptotic region. At the nuclear cusp, the VW functional is shown to be linear and the contribution of p-type orbitals to the KE density is theoretically derived and numerically demonstrated in the limit of infinite nuclear charge, as well in the semiclassical limit of neutral large atoms. In the latter case, it reaches 12 of the KE density. In the asymptotic region we find new exact constraints for meta Generalized Gradient Approximation (meta-GGA) exchange functionals: with an exchange enhancement factor proportional to $\sqrtα$, where $α$ is the common meta-GGA ingredient, both the exchange energy density and the potential are proportional to the exact ones. In addition, this describes exactly the large-gradient limit of quasi-two dimensional systems.

preprint2014arXiv

Laplacian-level kinetic energy approximations based on the fourth-order gradient expansion: Global assessment and application to the subsystem formulation of density functional theory

We test Laplacian-level meta-generalized gradient approximation (meta-GGA) non-interacting kinetic energy functionals based on the fourth-order gradient expansion (GE4). We consider several well known Laplacian-level meta-GGAs from literature (bare GE4, modified GE4, and the MGGA functional of Perdew and Constantin [Phys. Rev. B \textbf{75},155109 (2007)]), as well as two newly designed Laplacian-level kinetic energy functionals (named L0.4 and L0.6). First, a general assessment of the different functionals is performed, testing them for model systems (one-electron densities, Hooke's atom and different jellium systems), atomic and molecular kinetic energies as well as for their behavior with respect to density-scaling transformations. Finally, we assess, for the first time, the performance of the different functionals for Subsystem Density Functional Theory (DFT) calculations on non-covalently interacting systems. We find that the different Laplacian-level meta-GGA kinetic functionals may improve the description of different properties of electronic systems but no clear overall advantage is found over the best GGA functionals. Concerning Subsystem DFT calculations, the here proposed L0.4 kinetic energy functional is competitive with state-of-the-art GGAs, whereas all other Laplacian-level functionals fail badly. The performance of the Laplacian-level functionals is rationalized thanks to a two-dimensional reduced-gradient and reduced-Laplacian decomposition of the non-additive kinetic energy density.

preprint2014arXiv

Wave-function and density functional theory studies of dihydrogen complexes

We performed a benchmark study on a series of dihydrogen bond complexes and constructed a set of reference bond distances and interaction energies. The test set was employed to assess the performance of several wave-function correlated and density functional theory methods. We found that second-order correlation methods describe relatively well the dihydrogen complexes. However, for high accuracy inclusion of triple contributions is important. On the other hand, none of the considered density functional methods can simultaneously yield accurate bond lengths and interaction energies. However, we found that improved results can be obtained by the inclusion of non-local exchange contributions.

preprint2013arXiv

Construction of a general semilocal exchange-correlation hole model: Application to nonempirical meta-GGA functionals

Using a reverse-engineering method we construct a meta-generalized gradient approximation (meta-GGA) angle-averaged exchange-correlation hole model which has a general applicability. It satisfies known exact hole constraints and can exactly recover the exchange-correlation energy density of any reasonable meta-GGA exchange-correlation energy functional satisfying a minimal set of exact properties. The hole model is applied to several non-empirical meta-GGA functionals: the Tao-Perdew-Staroverov-Scuseria (TPSS), the revised TPSS (revTPSS) and the recently Balanced LOCalization (BLOC) meta-GGA of L.A. Constantin, E. Fabiano, and F. Della Sala, (J. Chem. Theory Comput. 9, 2256 (2013)). The empirical M06-L meta-GGA functional is also considered. Real-space analyses of atoms and ions as well as wave-vector analyses of jellium surface energies, show that the meta-GGA hole models, in particular the BLOC one, are very realistic and can reproduce many features of benchmark XC holes. In addition, the BLOC hole model can be used to estimate with good accuracy the Coulomb hole radius of small atoms and ions. Thus, the proposed meta-GGA hole models provide a valuable tool to validate in detail existing meta-GGA functionals, and can be further used in the development of DFT methods beyond the semilocal level of theory.

preprint2013arXiv

Meta-GGA exchange-correlation functional with a balanced treatment of nonlocality

We construct a meta-generalized-gradient approximation which properly balances the nonlocality contributions to the exchange and correlation at the semilocal level. This non-empirical functional shows good accuracy for a broad palette of properties (thermochemistry, structural properties) and systems (molecules, metal clusters, surfaces and bulk solids). The accuracy for several well known problems in electronic structure calculations, such as the bending potential of the silver trimer and the dimensional crossover of anionic gold clusters, is also demonstrated. The inclusion of empirical dispersion corrections is finally discussed and analyzed.

preprint2013arXiv

Spin-dependent gradient correction for more accurate atomization energies of molecules

We discuss, simplify, and improve the spin-dependent correction of L.A. Constantin et al., Phys. Rev. B 84, 233103, for atomization energies, and develop a density parameter of the form $v\propto |\nabla n|/n^{10/9}$, found from the statistical ensemble of one-electron densities. The here constructed exchange-correlation generalized gradient approximations (GGAs), named zvPBEsol and zvPBEint, show a broad applicability, and a good accuracy for many applications, because these corrected functionals significantly improve the atomization and binding energies of molecular systems, without worsening the behavior of the original functionals (PBEsol and PBEint) for other properties. This spin-dependent correction is also applied to meta-GGA dynamical correlation functionals combined with exact-exchange; in this case a significant (about 30%) improvement in atomization energies of small molecules is found.

preprint2013arXiv

Testing the broad applicability of the PBEint GGA functional and its one-parameter hybrid form

We review the performance of the PBEint GGA functional (Phys. Rev. B 2010, 82, 113104) recently proposed to improve the description of hybrid interfaces, and we introduce its one-parameter hybrid form (hPBEint). We consider different well established benchmarks for energetic and structural properties of molecular and solid-state systems as well as model systems and newly developed benchmark sets for dipole moments and metal-molecule interactions. We find that PBEint and hPBEint (with 16.67% Hartree-Fock exchange) yield the overall best performance, working well for most of the considered properties and systems and showing a balanced behavior for different problems. In particular, due to their well-balanced accuracy, they perform well for the description of hybrid metal-molecule interfaces.

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.