-
Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit
Authors:
S. Śmiga,
F. Della Sala,
P. Gori-Giorgi,
E. Fabiano
Abstract:
Adiabatic connection models (ACMs), which interpolate between the limits of weak and strong interaction, are powerful tools to build accurate exchange-correlation functionals. If the exact weak-interaction expansion from second-order perturbation theory is included, a self-consistent implementation of these functionals is challenging and still absent in the literature. In this work we fill this ga…
▽ More
Adiabatic connection models (ACMs), which interpolate between the limits of weak and strong interaction, are powerful tools to build accurate exchange-correlation functionals. If the exact weak-interaction expansion from second-order perturbation theory is included, a self-consistent implementation of these functionals is challenging and still absent in the literature. In this work we fill this gap by presenting a fully self-consistent-field (SCF) implementation of some popular ACM functionals. While using second-order perturbation theory at weak interactions, we have also introduced new generalised gradient approximations (GGA's), beyond the usual point-charge-plus-continuum model, for the first two leading terms at strong interactions, which are crucial to ensure robustness and reliability. We then assess the SCF-ACM functionals for molecular systems and for prototypical strong-correlation problems. We find that they perform well for both the total energy and the electronic density and that the impact of SCF orbitals is directly connected to the accuracy of the ACM functional form. For the H$_2$ dissociation the SCF-ACM functionals yield significant improvements with respect to standard functionals, also thanks to the use of the new GGA's for the strong-coupling functionals.
△ Less
Submitted 7 September, 2022; v1 submitted 23 February, 2022;
originally announced February 2022.
-
Methods to generate the reference total and Pauli kinetic potentials
Authors:
Szymon Śmiga,
Sylwia Siecińska,
Eduardo Fabiano
Abstract:
We have derived a new method which allows to compute the full and the Pauli reference kinetic potentials for atoms and molecules in a real space representation. This is done by applying the optimized effective potential (OEP) method to {the} Kohn-Sham non-interacting kinetic energy expression. Additionally, we have also derived a simplified OEP variant based on the common energy denominator approx…
▽ More
We have derived a new method which allows to compute the full and the Pauli reference kinetic potentials for atoms and molecules in a real space representation. This is done by applying the optimized effective potential (OEP) method to {the} Kohn-Sham non-interacting kinetic energy expression. Additionally, we have also derived a simplified OEP variant based on the common energy denominator approximation which has proven to give much more stable and robust results than the original OEP one. Moreover, we have also proved that at the solution point our approach is formally equivalent to the commonly used Bartolotti-Acharya formula.
△ Less
Submitted 7 May, 2020;
originally announced May 2020.
-
Controlling the Functional Properties of Oligothiophene Crystalline Nano/Micro-Fibers via Tailoring of the Self-Assembling Molecular Precursors
Authors:
Francesca Di Maria,
Mattia Zangoli,
Massimo Gazzano,
Eduardo Fabiano,
Denis Gentili,
Alberto Zanelli,
Andrea Fermi,
Giacomo Bergamini,
Davide Bonifazi,
Andrea Perinot,
Mario Caironi,
Raffaello Mazzaro,
Vittorio Morandi,
Giuseppe Gigli,
Andrea Liscio,
Giovanna Barbarella
Abstract:
Oligothiophenes are pi-conjugated semiconducting and fluorescent molecules whose self-assembly properties are widely investigated for application in organic electronics, optoelectronics, biophotonics and sensing. We report here an approach to the preparation of crystalline oli-gothiophene nano/micro-fibers based on the use of a sulfur overrich quaterthiophene building block, -T4S4-, containing in…
▽ More
Oligothiophenes are pi-conjugated semiconducting and fluorescent molecules whose self-assembly properties are widely investigated for application in organic electronics, optoelectronics, biophotonics and sensing. We report here an approach to the preparation of crystalline oli-gothiophene nano/micro-fibers based on the use of a sulfur overrich quaterthiophene building block, -T4S4-, containing in its covalent network all the information needed to promote the di-rectional, pi-pi stacking driven, self-assembly of Ar-T4S4-Ar oligomers into fibers with hierar-chical supramolecular arrangement from nano- to microscale. We show that when Ar varies from unsubstituted thiophene to thiophene substituted with electron withdrawing groups, a wide redistribution of the molecular electronic charge takes place without substantially affecting the aggregation modalities of the oligomer. In this way a structurally comparable series of fibers is obtained having progressively varying optical properties, redox potentials, photoconductivity and type of prevailing charge carriers (from p- to n-type). A thorough characterization of the fi-bers based on SEM, CD, CV, X-ray diffraction, UV-vis and PL spectroscopies, photoconductivi-ty and KPFM measurements is reported. With the aid of DFT calculations, combined with X-ray data, a model accounting for the growth of the fibers from molecular to nano- and microscale is proposed. We believe that the simple strategy outlined in this study allows to establish a straightforward correlation between the molecular structure of the components and the function-al properties of the corresponding self-assembled nano/micro-fibers.
△ Less
Submitted 5 December, 2019;
originally announced December 2019.
-
Strong-interaction limit of an adiabatic connection in Hartree-Fock theory
Authors:
Michael Seidl,
Sara Giarrusso,
Stefan Vuckovic,
Eduardo Fabiano,
Paola Gori-Giorgi
Abstract:
We show that the leading term in the strong-interaction limit of the adiabatic connection that has as weak-interaction expansion the Moeller-Plesset perturbation theory can be fully determined from a functional of the Hartree-Fock density. We analyze this functional and highlight similarities and differences with the strong-interaction limit of the density-fixed adiabatic connection case of Kohn-S…
▽ More
We show that the leading term in the strong-interaction limit of the adiabatic connection that has as weak-interaction expansion the Moeller-Plesset perturbation theory can be fully determined from a functional of the Hartree-Fock density. We analyze this functional and highlight similarities and differences with the strong-interaction limit of the density-fixed adiabatic connection case of Kohn-Sham density functional theory.
△ Less
Submitted 3 December, 2018; v1 submitted 26 October, 2018;
originally announced October 2018.
-
Investigation of the exchange-correlation potential of functionals based on the adiabatic connection interpolation
Authors:
E. Fabiano,
S. Śmiga,
S. Giarrusso,
T. J. Daas,
F. Della Sala,
I. Grabowski,
P. Gori-Giorgi
Abstract:
We have studied the correlation potentials produced by various adiabatic connection models (ACM) for several atoms and molecules. The results have been compared to accurate reference potentials (coupled cluster and quantum Monte Carlo results) as well as to state-of-the-art ab initio DFT approaches. We have found that all the ACMs yield correlation potentials that exhibit a correct behavior, quite…
▽ More
We have studied the correlation potentials produced by various adiabatic connection models (ACM) for several atoms and molecules. The results have been compared to accurate reference potentials (coupled cluster and quantum Monte Carlo results) as well as to state-of-the-art ab initio DFT approaches. We have found that all the ACMs yield correlation potentials that exhibit a correct behavior, quite resembling scaled second-order Görling-Levy (GL2) potentials, and including most of the physically meaningful features of the accurate reference data. The behavior and contribution of the strong-interaction limit potentials has also been investigated and discussed.
△ Less
Submitted 19 October, 2018;
originally announced October 2018.
-
Semilocal Pauli-Gaussian Kinetic Functionals for Orbital-Free Density Functional Theory Calculations of Solids
Authors:
L. A. Constantin,
E. Fabiano,
F. Della Sala
Abstract:
Kinetic energy (KE) approximations are key elements in orbital-free density functional theory. To date, the use of non-local functionals, possibly employing system dependent parameters, has been considered mandatory in order to obtain satisfactory accuracy for different solid-state systems, whereas semilocal approximations are generally regarded as unfit to this aim. Here, we show that instead pro…
▽ More
Kinetic energy (KE) approximations are key elements in orbital-free density functional theory. To date, the use of non-local functionals, possibly employing system dependent parameters, has been considered mandatory in order to obtain satisfactory accuracy for different solid-state systems, whereas semilocal approximations are generally regarded as unfit to this aim. Here, we show that instead properly constructed semilocal approximations, the Pauli-Gaussian (PG) KE functionals, especially at the Laplacian-level of theory, can indeed achieve similar accuracy as non-local functionals and can be accurate for both metals and semiconductors, without the need of system-dependent parameters.
△ Less
Submitted 24 July, 2018; v1 submitted 5 February, 2018;
originally announced February 2018.
-
Assessment of interaction-strength interpolation formulas for gold and silver clusters
Authors:
S. Giarrusso,
P. Gori-Giorgi,
F. Della Sala,
E. Fabiano
Abstract:
The performance of functionals based on the idea of interpolating between the weak and the strong-interaction limits the global adiabatic-connection integrand is carefully studied for the challenging case of noble-metal clusters. Different interpolation formulas are considered and various features of this approach are analyzed. It is found that these functionals, when used as a correlation correct…
▽ More
The performance of functionals based on the idea of interpolating between the weak and the strong-interaction limits the global adiabatic-connection integrand is carefully studied for the challenging case of noble-metal clusters. Different interpolation formulas are considered and various features of this approach are analyzed. It is found that these functionals, when used as a correlation correction to Hartree-Fock, are quite robust for the description of atomization energies, while performing less well for ionization potentials. Future directions that can be envisaged from this study and a previous one on main group chemistry are discussed.
△ Less
Submitted 30 January, 2018;
originally announced January 2018.
-
Electrostatic mechanophores in tuneable light-emitting piezo-polymer nanowires
Authors:
Luana Persano,
Andrea Camposeo,
Aleksandr V. Terentjevs,
Fabio Della Sala,
Eduardo Fabiano,
Martina Montinaro,
Dario Pisignano
Abstract:
Electromechanical coupling through piezoelectric polymer chains allows the emission of organic molecules in active nanowires to be tuned. This effect is evidenced by highly bendable arrays of counter-ion dye-doped nanowires made of a poly(vinylidenefluoride) copolymer. A reversible redshift of the dye emission is found upon the application of dynamic stress during highly accurate bending experimen…
▽ More
Electromechanical coupling through piezoelectric polymer chains allows the emission of organic molecules in active nanowires to be tuned. This effect is evidenced by highly bendable arrays of counter-ion dye-doped nanowires made of a poly(vinylidenefluoride) copolymer. A reversible redshift of the dye emission is found upon the application of dynamic stress during highly accurate bending experiments. By density functional theory calculations it is found that these photophysical properties are associated with mechanical stresses applied to electrostatically interacting molecular systems, namely to counterion-mediated states that involve light-emitting molecules as well as charged regions of piezoelectric polymer chains. These systems are an electrostatic class of supramolecular functional stress-sensitive units, which might impart new functionalities in hybrid molecular nanosystems and anisotropic nanostructures for sensing devices and soft robotics.
△ Less
Submitted 2 June, 2017;
originally announced June 2017.
-
Jellium-with-gap model applied to semilocal kinetic functionals
Authors:
L. A. Constantin,
E. Fabiano,
S. Śmiga,
F. Della Sala
Abstract:
We investigate a highly-nonlocal generalization of the Lindhard function, given by the jellium-with-gap model. We find a band-gap-dependent gradient expansion of the kinetic energy, which performs noticeably well for large atoms. Using the static linear response theory and the simplest semilocal model for the local band gap, we derive a non-empirical generalized gradient approximation (GGA) of the…
▽ More
We investigate a highly-nonlocal generalization of the Lindhard function, given by the jellium-with-gap model. We find a band-gap-dependent gradient expansion of the kinetic energy, which performs noticeably well for large atoms. Using the static linear response theory and the simplest semilocal model for the local band gap, we derive a non-empirical generalized gradient approximation (GGA) of the kinetic energy. This GGA kinetic energy functional is remarkably accurate for the description of weakly interacting molecular systems within the subsystem formulation of Density Functional Theory.
△ Less
Submitted 17 May, 2017;
originally announced May 2017.
-
Laplacian-dependent models of the kinetic energy density: Applications in subsystem density functional theory with meta-generalized gradient approximation functionals
Authors:
S. Śmiga,
E. Fabiano,
L. A. Constantin,
F. Della Sala
Abstract:
The development of semilocal models for the kinetic energy density (KED) is an important topic in density functional theory (DFT). This is especially true for subsystem DFT, where these models are necessary to construct the required non-additive embedding contributions. In particular, these models can also be efficiently employed to replace the exact KED in meta-Generalized Gradient Approximation…
▽ More
The development of semilocal models for the kinetic energy density (KED) is an important topic in density functional theory (DFT). This is especially true for subsystem DFT, where these models are necessary to construct the required non-additive embedding contributions. In particular, these models can also be efficiently employed to replace the exact KED in meta-Generalized Gradient Approximation (meta-GGA) exchange-correlation functionals allowing to extend the subsystem DFT applicability to the meta-GGA level of theory. Here, we present a two-dimensional scan of semilocal KED models as linear functionals of the reduced gradient and of the reduced Laplacian, for atoms and weakly-bound molecular systems. We find that several models can perform well but in any case the Laplacian contribution is extremely important to model the local features of the KED. Indeed a simple model constructed as the sum of Thomas-Fermi KED and 1/6 of the Laplacian of the density yields the best accuracy for atoms and weakly-bound molecular systems. These KED models are tested within subsystem DFT with various meta-GGA exchange-correlation functionals for non-bonded systems, showing a good accuracy of the method.
△ Less
Submitted 14 February, 2017;
originally announced February 2017.
-
The interaction-strength interpolation method for main-group chemistry: benchmarking, limitations, and perspectives
Authors:
E. Fabiano,
P. Gori-Giorgi,
M. Seidl,
F. Della Sala
Abstract:
We have tested the original interaction-strength-interpolation (ISI) exchange-correlation functional for main group chemistry. The ISI functional is based on an interpolation between the weak and strong coupling limits and includes exact-exchange as well as the Görling-Levy second-order energy. We have analyzed in detail the basis-set dependence of the ISI functional, its dependence on the ground-…
▽ More
We have tested the original interaction-strength-interpolation (ISI) exchange-correlation functional for main group chemistry. The ISI functional is based on an interpolation between the weak and strong coupling limits and includes exact-exchange as well as the Görling-Levy second-order energy. We have analyzed in detail the basis-set dependence of the ISI functional, its dependence on the ground-state orbitals, and the influence of the size-consistency problem. We show and explain some of the expected limitations of the ISI functional (i.e. for atomization energies), but also unexpected results, such as the good performance for the interaction energy of dispersion-bonded complexes when the ISI correlation is used as a correction to Hartree-Fock.
△ Less
Submitted 20 September, 2016;
originally announced September 2016.
-
Hartree potential dependent exchange functional
Authors:
L. A. Constantin,
E. Fabiano,
F. Della Sala
Abstract:
We introduce a novel non-local ingredient for the construction of exchange density functionals: the reduced Hartree parameter, which is invariant under the uniform scaling of the density and represents the exact exchange enhancement factor for one- and two-electron systems. The reduced Hartree parameter is used together with the conventional meta-generalized gradient approximation (meta-GGA) semil…
▽ More
We introduce a novel non-local ingredient for the construction of exchange density functionals: the reduced Hartree parameter, which is invariant under the uniform scaling of the density and represents the exact exchange enhancement factor for one- and two-electron systems. The reduced Hartree parameter is used together with the conventional meta-generalized gradient approximation (meta-GGA) semilocal ingredients (i.e. the electron density, its gradient and the kinetic energy density) to construct a new generation exchange functional, termed u-meta-GGA. This u-meta-GGA functional is exact for {the exchange of} any one- and two-electron systems, is size-consistent and non-empirical, satisfies the uniform density scaling relation, and recovers the modified gradient expansion derived from the semiclassical atom theory. For atoms, ions, jellium spheres, and molecules, it shows a good accuracy, being often better than meta-GGA exchange functionals. Our construction validates the use of the reduced Hartree ingredient in exchange-correlation functional development, opening the way to an additional rung in the Jacob's ladder classification of non-empirical density functionals.
△ Less
Submitted 31 August, 2016;
originally announced August 2016.
-
Kinetic and Exchange Energy Densities near the Nucleus
Authors:
L. A. Constantin,
E. Fabiano,
F. Della Sala
Abstract:
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 b…
▽ More
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.
△ Less
Submitted 31 March, 2016;
originally announced March 2016.
-
Semiclassical atom theory applied to solid-state physics
Authors:
L. A. Constantin,
A. Terentjevs,
F. Della Sala,
P. Cortona,
E. Fabiano
Abstract:
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 w…
▽ More
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.
△ Less
Submitted 25 January, 2016;
originally announced January 2016.
-
Semilocal density functional theory with correct surface asymptotics
Authors:
L. A. Constantin,
E. Fabiano,
J. M. Pitarke,
F. Della Sala
Abstract:
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 desc…
▽ More
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.
△ Less
Submitted 30 March, 2016; v1 submitted 21 January, 2016;
originally announced January 2016.
-
Assessment of the TCA functional in computational chemistry and solid-state physics
Authors:
E. Fabiano,
L. A. Constantin,
A. Terentjevs,
F. Della Sala,
P. Cortona
Abstract:
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 e…
▽ More
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.
△ Less
Submitted 6 October, 2015;
originally announced October 2015.
-
Accurate non-covalent interaction energies via an efficient MP2 scaling procedure
Authors:
E. Fabiano,
F. Della Sala,
I. Grabowski
Abstract:
Using the observed proportionality of CCSD(T) and MP2 correlation interaction energies [I. Grabowski, E. Fabiano, F. Della Sala, Phys. Chem. Chem. Phys. 15, 15485 (2013)] we propose a simple scaling procedure to compute accurate interaction energies of non-covalent complexes. Our method makes use of MP2 and CCSD(T) correlation energies, computed in relatively small basis sets, and fitted scaling c…
▽ More
Using the observed proportionality of CCSD(T) and MP2 correlation interaction energies [I. Grabowski, E. Fabiano, F. Della Sala, Phys. Chem. Chem. Phys. 15, 15485 (2013)] we propose a simple scaling procedure to compute accurate interaction energies of non-covalent complexes. Our method makes use of MP2 and CCSD(T) correlation energies, computed in relatively small basis sets, and fitted scaling coefficients to yield interaction energies of almost complete basis set limit CCSD(T) quality. Thanks to the good transferability of the scaling coefficients involved in the calculations, good results can be easily obtained for different intermolecular distances.
△ Less
Submitted 13 July, 2015;
originally announced July 2015.
-
Subsystem density functional theory with meta generalized gradient approximation exchange-correlation functionals
Authors:
S. Śmiga,
E. Fabiano,
S. Laricchia,
L. A. Constantin,
F. Della Sala
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 May, 2015;
originally announced May 2015.
-
Global hybrids from the semiclassical atom theory satisfying the local density linear response
Authors:
E. Fabiano,
L. A. Constantin,
P. Cortona,
F. Della Sala
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 January, 2015;
originally announced January 2015.
-
Kohn-Sham Kinetic Energy Density in the Nuclear and Asymptotic Regions: Deviations from the Von Weizsäcker Behavior and Applications to Density Functionals
Authors:
F. Della Sala,
E. Fabiano,
L. A. Constantin
Abstract:
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, a…
▽ More
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.
△ Less
Submitted 14 November, 2014;
originally announced November 2014.
-
Gradient-dependent upper bound for the exchange-correlation energy and application to density functional theory
Authors:
L. A. Constantin,
E. Fabiano,
A. Terentjevs,
F. Della Sala
Abstract:
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-…
▽ More
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.
△ Less
Submitted 21 January, 2015; v1 submitted 6 November, 2014;
originally announced November 2014.
-
Orbital-dependent second-order scaled-opposite-spin correlation functionals in the optimized effective potential method
Authors:
I. Grabowski,
E. Fabiano,
A. Teale,
S. Śmiga,
A. Buksztel,
F. Della Sala
Abstract:
The performance of correlated optimized effective potential (OEP) functionals based on the spin-resolved second-order correlation energy is analyzed. The relative importance of singly- and doubly- excited contributions as well as the effect of scaling the same- and opposite- spin components is investigated in detail comparing OEP results with Kohn-Sham (KS) quantities determined via an inversion p…
▽ More
The performance of correlated optimized effective potential (OEP) functionals based on the spin-resolved second-order correlation energy is analyzed. The relative importance of singly- and doubly- excited contributions as well as the effect of scaling the same- and opposite- spin components is investigated in detail comparing OEP results with Kohn-Sham (KS) quantities determined via an inversion procedure using accurate ab initio electronic densities. Special attention is dedicated in particular to the recently proposed scaled-opposite-spin OEP functional [I. Grabowski, E. Fabiano and F. Della Sala, Phys. Rev. B, 87, 075103, (2013)] which is the most advantageous from a computational point of view. We find that for high accuracy a careful, system dependent, selection of the scaling coefficient is required. We analyze several size-extensive approaches for this selection. Finally, we find that a composite approach, named OEP2-SOSh, based on a post-SCF rescaling of the correlation energy can yield high accuracy for many properties, being comparable with the most accurate OEP procedures previously reported in the literature but at substantially reduced computational effort.
△ Less
Submitted 29 July, 2014;
originally announced July 2014.
-
Wave-function and density functional theory studies of dihydrogen complexes
Authors:
E. Fabiano,
L. A. Constantin,
F. Della Sala
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 July, 2014;
originally announced July 2014.
-
Generalized gradient approximation correlation energy functionals based on the uniform electron gas with gap model
Authors:
E. Fabiano,
P. E. Trevisanutto,
A. Terentjevs,
L. A. Constantin
Abstract:
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 densit…
▽ More
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.
△ Less
Submitted 14 April, 2014;
originally announced April 2014.
-
Frozen density embedding with non-integer subsystems' particle numbers
Authors:
E. Fabiano,
S. Laricchia,
F. Della Sala
Abstract:
We extend the frozen density embedding theory to non-integer subsystems' particles numbers. Different features of this formulation are discussed, with special concern for approximate embedding calculations. In particular, we highlight the relation between the non-integer particle-number partition scheme and the resulting embedding errors. Finally, we provide a discussion of the implications of the…
▽ More
We extend the frozen density embedding theory to non-integer subsystems' particles numbers. Different features of this formulation are discussed, with special concern for approximate embedding calculations. In particular, we highlight the relation between the non-integer particle-number partition scheme and the resulting embedding errors. Finally, we provide a discussion of the implications of the present theory for the derivative discontinuity issue and the calculation of chemical reactivity descriptors.
△ Less
Submitted 18 March, 2014;
originally announced March 2014.
-
Laplacian-level kinetic energy approximations based on the fourth-order gradient expansion: Global assessment and application to the subsystem formulation of density functional theory
Authors:
S. Laricchia,
L. A. Constantin,
E. Fabiano,
F. Della Sala
Abstract:
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-…
▽ More
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.
△ Less
Submitted 18 March, 2014;
originally announced March 2014.
-
A Density Difference Based Analysis of Orbital--Dependent Exchange--Correlation Functionals
Authors:
I. Grabowski,
A. M. Teale,
E. Fabiano,
S. Smiga,
A. Buksztel,
F. Della Sala
Abstract:
We present a density difference based analysis for a range of orbital--dependent Kohn--Sham functionals. Results for atoms, some members of the neon isoelectronic series and small molecules are reported and compared with ab initio wave-function calculations. Particular attention is paid to the quality of approximations to the exchange--only optimized effective potential (OEP) approach: we consider…
▽ More
We present a density difference based analysis for a range of orbital--dependent Kohn--Sham functionals. Results for atoms, some members of the neon isoelectronic series and small molecules are reported and compared with ab initio wave-function calculations. Particular attention is paid to the quality of approximations to the exchange--only optimized effective potential (OEP) approach: we consider both the Localized Hartree Fock as well as the Krieger-Li-Iafrate methods. Analysis of density differences at the exchange--only level reveals the impact the approximations have on the resulting electronic densities. These differences are further quantified in terms of the ground state energies, frontier orbital energy differences and highest occupied orbital energies obtained. At the correlated level an OEP approach based on a perturbative second--order correlation energy expression is shown to deliver results comparable with those from traditional wave function approaches, making it suitable for use as a benchmark against which to compare standard density--functional approximations.
△ Less
Submitted 6 November, 2013;
originally announced November 2013.
-
Construction of a general semilocal exchange-correlation hole model: Application to nonempirical meta-GGA functionals
Authors:
L. A. Constantin,
E. Fabiano,
F. Della Sala
Abstract:
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.…
▽ More
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.
△ Less
Submitted 16 September, 2013;
originally announced September 2013.
-
Simple non-empirical procedure for spin-component-scaled MP2 methods applied to the calculation of dissociation energy curve of noncovalently-interacting systems
Authors:
I. Grabowski,
E. Fabiano,
F. Della Sala
Abstract:
We present a simple and non-empirical method to determine optimal scaling coefficients, within the (spin-component)-scaled MP2 approach, for calculating intermolecular potential energies of noncovalently-interacting systems. The method is based on an observed proportionality between (spin-component) MP2 and CCSD(T) energies for a wide range of intermolecular distances and allows to compute with hi…
▽ More
We present a simple and non-empirical method to determine optimal scaling coefficients, within the (spin-component)-scaled MP2 approach, for calculating intermolecular potential energies of noncovalently-interacting systems. The method is based on an observed proportionality between (spin-component) MP2 and CCSD(T) energies for a wide range of intermolecular distances and allows to compute with high accuracy a large portion of the dissociation curve at the cost of a single CCSD(T) calculation. The accuracy of the present procedure is assessed for a series of noncovalently-interacting test systems: the obtained results reproduce CCSD(T) quality in all cases and definitely outperform conventional MP2, CCSD and SCS-MP2 results. The difficult case of the Beryllium dimer is also considered.
△ Less
Submitted 28 August, 2013;
originally announced August 2013.
-
Meta-GGA exchange-correlation functional with a balanced treatment of nonlocality
Authors:
L. A. Constantin,
E. Fabiano,
F. Della Sala
Abstract:
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 problem…
▽ More
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.
△ Less
Submitted 16 May, 2013;
originally announced May 2013.
-
Semilocal and Hybrid Density Embedding Calculations of Ground-State Charge-Transfer Complexes
Authors:
S. Laricchia,
E. Fabiano,
F. Della Sala
Abstract:
We apply the frozen density embedding method, using a full relaxation of embedded densities through a freeze-and-thaw procedure, to study the electronic structure of several benchmark ground-state charge-transfer complexes, in order to assess the merits and limitations of the approach for this class of systems. The calculations are performed using both semilocal and hybrid exchange-correlation (XC…
▽ More
We apply the frozen density embedding method, using a full relaxation of embedded densities through a freeze-and-thaw procedure, to study the electronic structure of several benchmark ground-state charge-transfer complexes, in order to assess the merits and limitations of the approach for this class of systems. The calculations are performed using both semilocal and hybrid exchange-correlation (XC) functionals. The results show that embedding calculations using semilocal XC functionals yield rather large deviations with respect to the corresponding supermolecular calculations. Due to a large error cancellation effect, however, they can often provide a relatively good description of the electronic structure of charge-transfer complexes, in contrast to supermolecular calculations performed at the same level of theory. On the contrary, when hybrid XC functionals are employed, both embedding and supermolecular calculations agree very well with each other and with the reference benchmark results.
In conclusion, for the study of ground-state charge-transfer complexes via embedding calculations hybrid XC functionals are the method of choice due to their higher reliability and superior performance.
△ Less
Submitted 16 April, 2013;
originally announced April 2013.
-
Accurate ionization potential of gold anionic clusters from density functional theory and many-body perturbation theory
Authors:
A. Tanwar,
E. Fabiano,
P. E. Trevisanutto,
L. Chiodo,
F. Della Sala
Abstract:
We present a theoretical study of the ionization potential in small anionic gold clusters, using density functional theory, with and without exact-exchange, and many body perturbation theory, namely the G0W0 approach. We find that G0W0 is the best approach and correctly describes the first ionization potential with an accuracy of about 0.1 eV.
We present a theoretical study of the ionization potential in small anionic gold clusters, using density functional theory, with and without exact-exchange, and many body perturbation theory, namely the G0W0 approach. We find that G0W0 is the best approach and correctly describes the first ionization potential with an accuracy of about 0.1 eV.
△ Less
Submitted 16 April, 2013;
originally announced April 2013.
-
Accuracy of basis-set extrapolation schemes for DFT-RPA correlation energies in molecular calculations
Authors:
E. Fabiano,
F. Della Sala
Abstract:
We construct a reference benchmark set for atomic and molecular random-phase-approximation (RPA) correlation energies in a density functional theory (DFT) framework at the complete basis set limit. This set is used to evaluate the accuracy of some popular extrapolation schemes for RPA all-electron molecular calculations. The results indicate that for absolute energies accurate results, clearly out…
▽ More
We construct a reference benchmark set for atomic and molecular random-phase-approximation (RPA) correlation energies in a density functional theory (DFT) framework at the complete basis set limit. This set is used to evaluate the accuracy of some popular extrapolation schemes for RPA all-electron molecular calculations. The results indicate that for absolute energies accurate results, clearly outperforming raw data, are achievable with two-point extrapolation schemes based on quintuple- and sextuple-zeta basis sets. Moreover, we show that results in good agreement with the benchmark can also be also obtained by using a semiempirical extrapolation procedure based on quadruple- and quintuple-zeta basis sets. Finally, we analyze the performance of different extrapolation schemes for atomization energies.
△ Less
Submitted 21 February, 2013;
originally announced February 2013.
-
Spin-dependent gradient correction for more accurate atomization energies of molecules
Authors:
L. A. Constantin,
E. Fabiano,
F. Della Sala
Abstract:
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…
▽ More
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.
△ Less
Submitted 31 January, 2013;
originally announced January 2013.
-
Testing the broad applicability of the PBEint GGA functional and its one-parameter hybrid form
Authors:
E. Fabiano,
L. A. Constantin,
F. Della Sala
Abstract:
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 f…
▽ More
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.
△ Less
Submitted 24 January, 2013;
originally announced January 2013.
-
Relevance of coordinate and particle-number scaling in density functional theory
Authors:
Eduardo Fabiano,
Lucian A. Constantin
Abstract:
We discuss a $β$-dependent family of electronic density scalings of the form $n_λ(\R)=λ^{3β+1}\; n(λ^β\R)$ in the context of density functional theory. In particular, we consider the following special cases: the Thomas-Fermi scaling ($β=1/3$ and $λ\gg 1$), which is crucial for the semiclassical theory of neutral atoms; the uniform-electron-gas scaling ($β=-1/3$ and $λ\gg 1$), that is important in…
▽ More
We discuss a $β$-dependent family of electronic density scalings of the form $n_λ(\R)=λ^{3β+1}\; n(λ^β\R)$ in the context of density functional theory. In particular, we consider the following special cases: the Thomas-Fermi scaling ($β=1/3$ and $λ\gg 1$), which is crucial for the semiclassical theory of neutral atoms; the uniform-electron-gas scaling ($β=-1/3$ and $λ\gg 1$), that is important in the semiclassical theory of metallic clusters; the homogeneous density scaling ($β=0$) which can be related to the self-interaction problem in density functional theory when $λ\leq 1$; the fractional scaling ($β=1$ and $λ\leq 1$), that is important for atom and molecule fragmentation; and the strong-correlation scaling ($β=-1$ and $λ\gg 1$) that is important to describe the strong correlation limit.
The results of our work provide evidence for the importance of this family of scalings in semiclassical and quantum theory of electronic systems, and indicate that these scaling properties must be considered as important constraints in the construction of new approximate density functionals. We also show, using the uniform-electron-gas scaling, that the curvature energy of metallic clusters is related to the second-order gradient expansion of kinetic and exchange-correlation energies.
△ Less
Submitted 17 January, 2013;
originally announced January 2013.
-
Phonons and gluons in the crystalline color superconducting phase of QCD
Authors:
R. Casalbuoni,
E. Fabiano,
R. Gatto,
M. Mannarelli,
G. Nardulli
Abstract:
The High Density Effective Theory formalism is used to calculate the low energy properties of the phonons and gluons in the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase of two-flavor QCD. The effective phonon Lagrangian for the cubic crystal structure, which is favored according to a recent study, depends, at the second order in the derivatives, on three parameters which we calculate in this pa…
▽ More
The High Density Effective Theory formalism is used to calculate the low energy properties of the phonons and gluons in the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase of two-flavor QCD. The effective phonon Lagrangian for the cubic crystal structure, which is favored according to a recent study, depends, at the second order in the derivatives, on three parameters which we calculate in this paper. We also compute for soft momenta the effective lagrangian for the gluons of the unbroken SU(2)_c group, both for a single plane wave structure and for the cube. In both cases the Meissner gluon mass vanishes as in the case of complete isotropy; on the other hand there is a partial Debye screening due to the existence of blocking regions on the Fermi spheres. For the single plane wave structure the gluon dielectric tensor is non isotropic, while it is isotropic for the cubic crystal, in spite of the intrinsic anisotropy of the structure.
△ Less
Submitted 13 August, 2002;
originally announced August 2002.