Showing 1–50 of 53 results for author: Della Sala, F

Towards adiabatic-connection interpolation model with broader applicability

Authors: Lucian A. Constantin, Szymon Śmiga, Fabio Della Sala

Abstract: The Adiabatic Connection Integrand Interpolation (ACII) method represents a general path for calculating correlation energies in electronic systems within the Den sity Functional Theory. ACII functionals include both exact-exchange and the second-order correlation energy, as well as an interpolating function toward the strictly-correlated electron (SCE) regime. Several interpolating functions have… ▽ More The Adiabatic Connection Integrand Interpolation (ACII) method represents a general path for calculating correlation energies in electronic systems within the Den sity Functional Theory. ACII functionals include both exact-exchange and the second-order correlation energy, as well as an interpolating function toward the strictly-correlated electron (SCE) regime. Several interpolating functions have been proposed in the last years targeting different properties, yet an accurate ACII approach with broad applicability is sti ll missing. Recently, we have proposed an ACII functional that was made accurate for the three-dimensional (3D) uniform electron gas as well as for model metal clusters. In this work we present an ACII functional (named genISI2) which is very accurate for both three-dimensional (3D) and two-dimensional (2D) uniform electron gases and for the q uasi-2D infinite barrier model, where most of the exchange-correlation functionals fail badly, as well as for strongly correlated two-electrons systems. Using the exact-exchange Kohn-Sham orbitals, we have also assessed the genISI2 for various molecular systems, showing a superior performance with respect to the o ther ACII methods for total energies, atomization energies, and ionization potentials. The genISI2 functional can thus find application in a broad range of systems and properties. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas

Authors: Lucian A. Constantin, Subrata Jana, Szymon Śmiga, Fabio Della Sala

Abstract: The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being in principle exact in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct the ge… ▽ More The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being in principle exact in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct the genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z $\leq$ 912 and for the non-relativistic noble atoms with Z $\leq$ 290. For the jellium clusters, the genISI is remarkably accurate, while for the noble atoms, it shows a good performance, similar to other ISI-like methods. Then, the genISI functional can open the path using the ISI-like method in solid-state calculations. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Regularized and Opposite spin-scaled functionals from Møller-Plesset adiabatic connection -- higher accuracy at lower cost

Authors: Kimberly J. Daas, Derk P. Kooi, Nina C. Peters, Eduardo Fabiano, Fabio Della Sala, Paola Gori-Giorgi, Stefan Vuckovic

Abstract: Non-covalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation… ▽ More Non-covalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields $c_{\rm os}κ_{\rm os}$-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies. With the $N^4$ formal scaling, $c_{\rm os}κ_{\rm os}$-SPL2, is competitive or often outperforms more expensive dispersion-corrected double hybrids for NCIs.The accuracy of $c_{\rm os}κ_{\rm os}$-SPL2 particularly shines for anionic halogen bonded complexes, where it surpasses standard dispersion-corrected DFT by a factor of 3 to 5. △ Less

Submitted 7 July, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 12 pages + 5 SI, 8 figures + 6 SI

Gaussian expansion of Yukawa non-local kinetic energy functionals: application to metal clusters

Authors: F. Sarcinella, S. Śmiga, F. Della Sala, E. Fabiano

Abstract: The development of kinetic energy (KE) functionals is one of the current challenges in density functional theory (DFT). The Yukawa non-local KE functionals [Phys. Rev. B 103, 155127 (2021)] have been shown to describe accurately the Lindhard response of the homogeneous electron gas (HEG) directly in the real space, without any step in the reciprocal space. However, the Yukawa kernel employs an exp… ▽ More The development of kinetic energy (KE) functionals is one of the current challenges in density functional theory (DFT). The Yukawa non-local KE functionals [Phys. Rev. B 103, 155127 (2021)] have been shown to describe accurately the Lindhard response of the homogeneous electron gas (HEG) directly in the real space, without any step in the reciprocal space. However, the Yukawa kernel employs an exponential function which cannot be efficiently represented in conventional Gaussian-based quantum chemistry codes. Here, we present an expansion of the Yukawa kernel in Gaussian functions. We show that for the HEG this expansion is independent of the electronic density, and that for general finite systems the accuracy can be easily tuned. Finally, we present results for atomistic sodium clusters of different sizes, showing that simple Yukawa functionals can give superior accuracy as compared to semilocal functionals. △ Less

Submitted 3 April, 2023; originally announced April 2023.

Comments: 9 pages, 6 figures

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.

Comments: 40 pages, 6 figures

Noncovalent interactions from models for the Møller-Plesset adiabatic connection

Authors: Timothy J. Daas, Eduardo Fabiano, Fabio Della Sala, Paola Gori-Giorgi, Stefan Vuckovic

Abstract: Given the omnipresence of non-covalent interactions (NCIs), their accurate simulations are of crucial importance across various scientific disciplines. Here we construct accurate models for the description of NCIs by an interpolation along the Møller-Plesset adiabatic connection (MP AC). Our interpolation approximates the correlation energy, by recovering MP2 at small coupling strengths and the co… ▽ More Given the omnipresence of non-covalent interactions (NCIs), their accurate simulations are of crucial importance across various scientific disciplines. Here we construct accurate models for the description of NCIs by an interpolation along the Møller-Plesset adiabatic connection (MP AC). Our interpolation approximates the correlation energy, by recovering MP2 at small coupling strengths and the correct large-coupling strength expansion of the MP AC, recently shown to be a functional of the Hartree-Fock density. Our models are size consistent for fragments with non-degenerate ground states, have the same cost as double hybrids and require no dispersion corrections to capture NCIs accurately. These interpolations greatly reduce large MP2 errors for typical $π$-stacking complexes (e.g., benzene-pyridine dimers) and for the L7 dataset. They are also competitive with state-of-the-art dispersion enhanced functionals and can even significantly outperform them for a variety of datasets, such as CT7 and L7. △ Less

Submitted 10 May, 2021; v1 submitted 10 April, 2021; originally announced April 2021.

Comments: 30 pages, 6 figures + SI (11 pages, 8 figures), added small changes in the main text

Laplacian-Level Quantum Hydrodynamic Theory for Plasmonics

Authors: Henrikh M. Baghramyan, Fabio Della Sala, Cristian Ciracì

Abstract: An accurate description of the optical response of subwavelength metallic particles and nanogap structures is a key problem of plasmonics. Quantum hydrodynamic theory (QHT) has emerged as a powerful method to calculate the optical response of metallic nanoparticles (NPs) since it takes into account nonlocality and spill-out effects. Nevertheless, the absorption spectra of metallic NPs obtained wit… ▽ More An accurate description of the optical response of subwavelength metallic particles and nanogap structures is a key problem of plasmonics. Quantum hydrodynamic theory (QHT) has emerged as a powerful method to calculate the optical response of metallic nanoparticles (NPs) since it takes into account nonlocality and spill-out effects. Nevertheless, the absorption spectra of metallic NPs obtained with conventional QHT, i.e., incorporating Thomas-Fermi (TF) and von Weizsäcker (vW) kinetic energy (KE) contributions, can be affected by several spurious resonances at energies higher than the main localized surface plasmon (LSP). These peaks are not present in reference time-dependent density-functional-theory (TD-DFT) spectra, where, instead, only a broad shoulder exists. Moreover, we show here that these peaks incorrectly reduce the LSP peak intensity and have a strong dependence on the simulation domain size so that a proper calculation of QHT absorption spectra can be problematic. In this article, we introduce a more general QHT method accounting for KE contributions depending on the Laplacian of the electronic density ($q$), thus, beyond the gradient-only dependence of the TFvW functional. We show that employing a KE functional with a term proportional to $q^2$ results in an absorption spectrum free of spurious peaks, with LSP resonance of correct intensity and numerically stable Bennett state. Finally, we present a novel Laplacian-level KE energy functional that is very accurate for the description of the optical properties of NPs with different sizes as well as for dimers. △ Less

Submitted 9 March, 2021; v1 submitted 6 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. X 11, 011049 (2021)

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.

Comments: 16 pages, 10 figures

Nanowire-Intensified MEF in Hybrid Polymer-Plasmonic Electrospun Filaments

Authors: Andrea Camposeo, Radoslaw Jurga, Maria Moffa, Alberto Portone, Francesco Cardarelli, Fabio Della Sala, Cristian Ciracì, Dario Pisignano

Abstract: Hybrid polymer-plasmonic nanostructures might combine high enhancement of localized fields from metal nanoparticles with light confinement and long-range transport in subwavelength dielectric structures. Here we report on the complex behavior of fluorophores coupling to Au nanoparticles within polymer nanowires, which features localized metal-enhanced fluorescence (MEF) with unique characteristics… ▽ More Hybrid polymer-plasmonic nanostructures might combine high enhancement of localized fields from metal nanoparticles with light confinement and long-range transport in subwavelength dielectric structures. Here we report on the complex behavior of fluorophores coupling to Au nanoparticles within polymer nanowires, which features localized metal-enhanced fluorescence (MEF) with unique characteristics compared to conventional structures. The intensification effect when the particle is placed in the organic filaments is remarkably higher with respect to thin films of comparable thickness, thus highlighting a specific, nanowire-related enhancement of MEF effects. A dependence on the confinement volume in the dielectric nanowire is also evidenced, with MEF significantly increasing upon reducing the wire diameter. These findings are rationalized by finite element simulations, predicting a position-dependent enhancement of the quantum yield of fluorophores embedded in the fibers. Calculation of the ensemble-averaged fluorescence enhancement unveils the possibility of strongly enhancing the overall emission intensity for structures with size twice the diameter of the embedded metal particles. These new, hybrid fluorescent systems with localized enhanced emission, as well as the general Nanowire-Intensified MEF effect associated to them, are highly relevant for developing nanoscale light-emitting devices with high efficiency and inter-coupled through nanofiber networks, highly sensitive optical sensors, and novel laser architectures. △ Less

Submitted 15 April, 2018; originally announced April 2018.

Comments: 29 pages, 12 figures, Small (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.

Comments: 19 pages, 3 figures

Journal ref: J. Phys. Chem. Lett. 9, 4385-4390 (2018)

Restoring size consistency of approximate functionals constructed from the adiabatic connection

Authors: S. Vuckovic, P. Gori-Giorgi, F. Della Sala, E. Fabiano

Abstract: Approximate exchange-correlation functionals built by modeling in a non-linear way the adiabatic connection (AC) integrand of density functional theory have many attractive features, being virtually parameters-free and satisfying different exact properties, but they also have a fundamental flaw: they violate the size-consistency condition, crucial to evaluate interaction energies of molecular syst… ▽ More Approximate exchange-correlation functionals built by modeling in a non-linear way the adiabatic connection (AC) integrand of density functional theory have many attractive features, being virtually parameters-free and satisfying different exact properties, but they also have a fundamental flaw: they violate the size-consistency condition, crucial to evaluate interaction energies of molecular systems. We show that size consistency in the AC-based functionals can be restored in a very simple way at no extra computational cost. Results on a large set of benchmark molecular interaction energies show that functionals based on the interaction strength interpolation approximations are significantly more accurate than the second-order perturbation theory. △ Less

Submitted 24 May, 2018; v1 submitted 6 February, 2018; originally announced February 2018.

Comments: 6 pages, 3 figures

Non-local kinetic energy functional from the Jellium-with-gap model: applications to Orbital-Free Density Functional Theory

Authors: L. A. Constantin, E. Fabiano, F. Della Sala

Abstract: Orbital-Free Density Functional Theory (OF-DFT) promises to describe the electronic structure of very large quantum systems, being its computational cost linear with the system size. However, the OF-DFT accuracy strongly depends on the approximation made for the kinetic energy (KE) functional. To date, the most accurate KE functionals are non-local functionals based on the linear-response kernel o… ▽ More Orbital-Free Density Functional Theory (OF-DFT) promises to describe the electronic structure of very large quantum systems, being its computational cost linear with the system size. However, the OF-DFT accuracy strongly depends on the approximation made for the kinetic energy (KE) functional. To date, the most accurate KE functionals are non-local functionals based on the linear-response kernel of the homogeneous electron gas, i.e. the jellium model. Here, we use the linear-response kernel of the jellium-with-gap model, to construct a simple non-local KE functional (named KGAP) which depends on the band gap energy. In the limit of vanishing energy-gap (i.e. in the case of metals), the KGAP is equivalent to the Smargiassi-Madden (SM) functional, which is accurate for metals. For a series of semiconductors (with different energy-gaps), the KGAP performs much better than SM, and results are close to the state-of-the-art functionals with complicated density-dependent kernels. △ Less

Submitted 23 May, 2018; v1 submitted 5 February, 2018; originally announced February 2018.

Comments: 11 pages, 4 figure

Journal ref: Phys. Rev. B 97, 205137 (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.

Comments: 16 pages, 6 figures

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.

Comments: 32 pages, 18 figures

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.

Comments: 9 pages, 2 figures

Journal ref: Phys. Rev. B 95, 115153 (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.

Comments: 11 pages, 5 figures

Journal ref: J. Chem. Phys. 146, 064105 (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.

Comments: 20 pages, 20 figures

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.

Comments: 14 pages, 9 figures

Journal ref: J. Chem. Phys. 145, 084110 (2016)

Emerging giant resonant exciton induced by Ta-substitution in anatase TiO$_{2}$: a tunable correlation effect

Authors: Z. Yong, P. E. Trevisanutto, L. Chiodo, I. Santoso, A. R. Barman, T. C. Asmara, S. Dhar, A. Kotlov, A. Terentjevs, F. Della Sala, V. Olevano, M. Rübhausen, T. Venkatesan, A. Rusydi

Abstract: Titanium dioxide (TiO$_2$) has rich physical properties with potential implications in both fundamental physics and new applications. Up-to-date, the main focus of applied research is to tune its optical properties, which is usually done via doping and/or nano-engineering. However, understanding the role of $d$-electrons in materials and possible functionalization of $d$-electron properties are st… ▽ More Titanium dioxide (TiO$_2$) has rich physical properties with potential implications in both fundamental physics and new applications. Up-to-date, the main focus of applied research is to tune its optical properties, which is usually done via doping and/or nano-engineering. However, understanding the role of $d$-electrons in materials and possible functionalization of $d$-electron properties are still major challenges. Herewith, within a combination of an innovative experimental technique, high energy optical conductivity, and of the state-of-the-art {\it ab initio} electronic structure calculations, we report an emerging, novel resonant exciton in the deep ultraviolet region of the optical response. The resonant exciton evolves upon low concentration Ta-substitution in anatase TiO$_{2}$ films. It is surprisingly robust and related to strong electron-electron and electron-hole interactions. The $d$- and $f$- orbitals localization, due to Ta-substitution, plays an unexpected role, activating strong electronic correlations and dominating the optical response under photoexcitation. Our results shed light on a new optical phenomenon in anatase TiO$_{2}$ films and on the possibility of tuning electronic properties by Ta substitution. △ Less

Submitted 22 May, 2016; originally announced May 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.

Comments: 11 pages, 7 figures

Enhancement of radiative processes in nanofibers with embedded plasmonic nanoparticles

Authors: Radoslaw Jurga, Fabio Della Sala, Dario Pisignano, Cristian Ciracì

Abstract: Efficient manipulation and long distance transport of single-photons is a key component in nanoscale quantum optics. In this letter, we study the emission properties of an individual light emitter placed into a nanofiber and coupled to a metallic nanoparticle. We find that plasmonic field enhancement together with the nanofiber optical confinement uniquely and synergistically contribute to an over… ▽ More Efficient manipulation and long distance transport of single-photons is a key component in nanoscale quantum optics. In this letter, we study the emission properties of an individual light emitter placed into a nanofiber and coupled to a metallic nanoparticle. We find that plasmonic field enhancement together with the nanofiber optical confinement uniquely and synergistically contribute to an overall increase of emission rates as well as quantum yields. △ Less

Submitted 3 February, 2016; originally announced February 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.

Comments: 10 pages, 7 figures

Journal ref: Phys. Rev. B 93, 045126 (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.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 93, 115127 (2016)

Quantum Hydrodynamic Theory for Plasmonics: Impact of the Electron Density Tail

Authors: Cristian Ciracì, Fabio Della Sala

Abstract: Multiscale plasmonic systems e.g. extended metallic nanostructures with sub-nanometer inter-distances) play a key role in the development of next-generation nano-photonic devices. An accurate modeling of the optical interactions in these systems requires an accurate description of both quantum effects and far-field properties. Classical electromagnetism can only describe the latter, while Time… ▽ More Multiscale plasmonic systems e.g. extended metallic nanostructures with sub-nanometer inter-distances) play a key role in the development of next-generation nano-photonic devices. An accurate modeling of the optical interactions in these systems requires an accurate description of both quantum effects and far-field properties. Classical electromagnetism can only describe the latter, while Time-Dependent Density Functional Theory (TD-DFT) can provide a full first-principles quantum treatment. However, TD-DFT becomes computationally prohibitive for sizes that exceed few nanometers, which are irrelevant for practical applications. In this article, we introduce a method based on the quantum hydrodynamic theory (QHT), which includes non-local contributions of the kinetic energy and the correct asymptotic description of the electron density. We show that our QHT method can predict both plasmon energy and spill-out effects in metal nanoparticles in excellent agreement with TD-DFT predictions, thus allowing a reliable and efficient calculations of both quantum and far-field properties in multiscale plasmonic systems. △ Less

Submitted 9 May, 2016; v1 submitted 7 January, 2016; originally announced January 2016.

Journal ref: Phys. Rev. B 93, 205405 (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.

Comments: 14 pages, 6 figures

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.

Comments: 8 pages, 4 figures

Journal ref: Chem. Phys. Lett. 635, 262 (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.

Comments: 14 pages, 3 figures

Journal ref: J. Chem. Phys. 142, 154121 (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.

Comments: 12 pages, 4 figures

Journal ref: J. Chem. Theory Comput. 11, 122 (2015)

Optical conductivity renormalization of graphene on SrTiO$_{3}$ due to resonant excitonic effects mediated by Ti 3\textit{d} orbitals

Authors: Pranjal Kumar Gogoi, Paolo E. Trevisanutto, Ming Yang, Iman Santoso, Teguh Citra Asmara, Aleksandrs Terentjevs, Fabio Della Sala, Mark B. H. Breese, T. Venkatesan, Yuan Ping Feng, Kian Ping Loh, Antonio H. Castro Neto, Andrivo Rusydi

Abstract: We present evidence of a drastic renormalization of the optical conductivity of graphene on SrTiO$_3$ resulting in almost full transparency in the ultraviolet region. These findings are attributed to resonant excitonic effects further supported by \emph{ab initio} Bethe-Salpeter equation and density functional theory calculations. The ($π$,$π$*)-orbitals of graphene and Ti-3\textit{d} $t_{2g}$ orb… ▽ More We present evidence of a drastic renormalization of the optical conductivity of graphene on SrTiO$_3$ resulting in almost full transparency in the ultraviolet region. These findings are attributed to resonant excitonic effects further supported by \emph{ab initio} Bethe-Salpeter equation and density functional theory calculations. The ($π$,$π$*)-orbitals of graphene and Ti-3\textit{d} $t_{2g}$ orbitals of SrTiO$_3$ are strongly hybridized and the interactions of electron-hole states residing in those orbitals play dominant role in the graphene optical conductivity. These interactions are present much below the optical band gap of bulk SrTiO$_3$. These results open a possibility of manipulating interaction strengths in graphene via \textit {d}-orbitals which could be crucial for optical applications. △ Less

Submitted 3 December, 2014; originally announced December 2014.

Comments: 9 pages, 11 figures, Phys. Rev. B (accepted)

Journal ref: Phys. Rev. B 91, 035424 (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.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. B 91, 035126 (2015)

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.

Comments: 5 pages, 2 figures

Journal ref: Phys. Rev. B 91, 041120(R) (2015)

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.

Comments: 14 pages, 4 figures

Journal ref: J. Chem. Phys. 141, 024113 (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.

Comments: 15 pages, 7 figures

Journal ref: J. Chem. Theory Comput. 10, 3151 (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.

Comments: 14 pages, 4 figures

Journal ref: J. Chem. Phys. 140, 114101 (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.

Comments: 19 pages, 6 figures

Journal ref: J. Chem. Theory. Comput. 10, 164 (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.

Comments: 12 pages, 3 figures

Journal ref: Mol. Phys. 112, 700 (2014)

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.

Comments: 11 pages, 10 figures

Journal ref: Phys. Rev. B 88, 125112 (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.

Comments: 9 pages, 7 figures

Journal ref: Phys. Chem. Chem. Phys. 15, 15485 (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.

Comments: 10 pages, 4 figures

Journal ref: J. Chem. Theory Comput. 9, 2256 (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.

Comments: 14 pages, 7 figures

Journal ref: J. Chem. Phys. 138,124112 (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.

Comments: 6 pages, 1 figure

Journal ref: Eur. Phys. J. B 86, 161 (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.

Comments: 10 pages, 2 figures

Journal ref: Theor. Chem. Acc. 131, 1278 (2012)

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.

Comments: 10 pages, 6 figures

Journal ref: J. Chem. Phys. 137, 194105 (2012)

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.

Comments: 13 pages, 3 figures. Accepted for publication in Int. J. Quant. Chem. (2013), DOI: 10.1002/qua.24042

Journal ref: Int. J. Quantum Chem. 113, 673 (2013)

Optical spectra of solids obtained by time-dependent density-functional theory with the jellium-with-gap model exchange-correlation kernel

Authors: Paolo E. Trevisanutto, Aleksandrs Terentjevs, Lucian A. Constantin, Valerio Olevano, Fabio Della Sala

Abstract: Within the framework of ab initio time-dependent density-functional theory (TD-DFT), we propose a static approximation to the exchange-correlation kernel based on the jellium-with-gap model. This kernel accounts for electron-hole interactions and it is able to address both strongly bound excitons and weak excitonic effects. TD-DFT absorption spectra of several bulk materials (both semiconductor an… ▽ More Within the framework of ab initio time-dependent density-functional theory (TD-DFT), we propose a static approximation to the exchange-correlation kernel based on the jellium-with-gap model. This kernel accounts for electron-hole interactions and it is able to address both strongly bound excitons and weak excitonic effects. TD-DFT absorption spectra of several bulk materials (both semiconductor and insulators) are reproduced in very good agreement with the experiments and with a low computational cost. △ Less

Submitted 26 October, 2012; originally announced October 2012.

Comments: 5 pages, 3 figures, 1 table

Journal ref: Phys. Rev. B 87, 205143 (2013)

Range separated functionals in the density functional based tight binding method: Formalism

Authors: Thomas A. Niehaus, Fabio Della Sala

Abstract: A generalization of the density-functional based tight-binding method (DFTB) for the use with range-separated exchange-correlation functionals is presented. It is based on the Generalized Kohn-Sham (GKS) formalism and employs the density matrix as basic variable in the expansion of the energy functional, in contrast to the traditional DFTB scheme. The GKS-TB equations are derived and appropriate i… ▽ More A generalization of the density-functional based tight-binding method (DFTB) for the use with range-separated exchange-correlation functionals is presented. It is based on the Generalized Kohn-Sham (GKS) formalism and employs the density matrix as basic variable in the expansion of the energy functional, in contrast to the traditional DFTB scheme. The GKS-TB equations are derived and appropriate integral approximations are discussed in detail. Implementation issues and numerical aspects of the new scheme are also covered. △ Less

Submitted 8 November, 2011; originally announced November 2011.

Electrostatic Field Driven Alignment of Organic Oligomers on ZnO Surfaces

Authors: F. Della Sala, S. Blumstengel, F. Henneberger

Abstract: We study the physisorption of organic oligomers on the ZnO($10\bar{1}0$) surface using first-principles density-functional theory and non-empirical embedding methods. We find that both in-plane location and orientation of the molecules are completely determined by the coupling of their quadrupole moments to the periodic dipolar electric field present at the semiconductor surface. The adsorption is… ▽ More We study the physisorption of organic oligomers on the ZnO($10\bar{1}0$) surface using first-principles density-functional theory and non-empirical embedding methods. We find that both in-plane location and orientation of the molecules are completely determined by the coupling of their quadrupole moments to the periodic dipolar electric field present at the semiconductor surface. The adsorption is associated with the formation of a molecular dipole moment perpendicular to the surface, which bears an unexpected linear relation to the molecule-substrate interaction energy. Long oligomers such as sexiphenyl become well-aligned with stabilization energies of several 100 meV along rows of positive electric field, in full agreement with recent experiments. These findings define a new route towards the realization of highly-ordered self-assembled arrays of oligomers/polymers on ZnO($10\bar{1}0$) and similar surfaces. △ Less

Submitted 21 January, 2011; originally announced January 2011.

Comments: 4 pages, 4 figures

Journal ref: Phys. Rev. Lett. 107, 146401 (2011)

Absorption and Fluorescence Properties of Oligothiophene Biomarkers from Long-Range-Corrected Time-Dependent Density Functional Theory

Authors: Bryan M. Wong, Manuel Piacenza, Fabio Della Sala

Abstract: The absorption and fluorescence properties in a class of oligothiophene push-pull biomarkers are investigated with a long-range-corrected (LC) density functional method. Using linear response time-dependent density functional theory (TDDFT), we calculate excitation energies, fluorescence energies, oscillator strengths, and excited-state dipole moments. To benchmark and assess the quality of the… ▽ More The absorption and fluorescence properties in a class of oligothiophene push-pull biomarkers are investigated with a long-range-corrected (LC) density functional method. Using linear response time-dependent density functional theory (TDDFT), we calculate excitation energies, fluorescence energies, oscillator strengths, and excited-state dipole moments. To benchmark and assess the quality of the LC-TDDFT formalism, an extensive comparison is made between LC-BLYP excitation energies and approximate coupled cluster singles and doubles (CC2) calculations. When using a properly-optimized value of the range parameter, "mu", we find that the LC technique provides an accurate description of charge-transfer excitations as a function of biomarker size and chemical functionalization. In contrast, we find that re-optimizing the fraction of Hartree Fock exchange in conventional hybrid functionals still yields an inconsistent description of excitation energies and oscillator strengths for the two lowest excited states in our series of biomarkers. The results of the present study emphasize the importance of a distance-dependent contribution of exchange in TDDFT for investigating excited-state properties. △ Less

Submitted 24 April, 2009; originally announced April 2009.

Comments: Accepted by Physical Chemistry Chemical Physics

Journal ref: Physical Chemistry Chemical Physics, 11, 4498 (2009)

Entanglement of electrons in interacting molecules

Authors: Tina A. C. Maiolo, Fabio Della Sala, Luigi Martina, Giulio Soliani

Abstract: Quantum entanglement is a concept commonly used with reference to the existence of certain correlations in quantum systems that have no classical interpretation. It is a useful resource to enhance the mutual information of memory channels or to accelerate some quantum processes as, for example, the factorization in Shor's Algorithm. Moreover, entanglement is a physical observable directly measur… ▽ More Quantum entanglement is a concept commonly used with reference to the existence of certain correlations in quantum systems that have no classical interpretation. It is a useful resource to enhance the mutual information of memory channels or to accelerate some quantum processes as, for example, the factorization in Shor's Algorithm. Moreover, entanglement is a physical observable directly measured by the von Neumann entropy of the system. We have used this concept in order to give a physical meaning to the electron correlation energy in systems of interacting electrons. The electronic correlation is not directly observable, since it is defined as the difference between the exact ground state energy of the many--electrons Schroedinger equation and the Hartree--Fock energy. We have calculated the correlation energy and compared with the entanglement, as functions of the nucleus--nucleus separation using, for the hydrogen molecule, the Configuration Interaction method. Then, in the same spirit, we have analyzed a dimer of ethylene, which represents the simplest organic conjugate system, changing the relative orientation and distance of the molecules, in order to obtain the configuration corresponding to maximum entanglement. △ Less

Submitted 27 October, 2006; originally announced October 2006.

Comments: 15 pages, 7 figures, standard latex

Quasiparticle energies for large molecules: a tight-binding GW approach

Authors: T. A. Niehaus, M. Rohlfing, F. Della Sala, A. Di Carlo, Th. Frauenheim

Abstract: We present a tight-binding based GW approach for the calculation of quasiparticle energy levels in confined systems such as molecules. Key quantities in the GW formalism like the microscopic dielectric function or the screened Coulomb interaction are expressed in a minimal basis of spherically averaged atomic orbitals. All necessary integrals are either precalculated or approximated without reso… ▽ More We present a tight-binding based GW approach for the calculation of quasiparticle energy levels in confined systems such as molecules. Key quantities in the GW formalism like the microscopic dielectric function or the screened Coulomb interaction are expressed in a minimal basis of spherically averaged atomic orbitals. All necessary integrals are either precalculated or approximated without resorting to empirical data. The method is validated against first principles results for benzene and anthracene, where good agreement is found for levels close to the frontier orbitals. Further, the size dependence of the quasiparticle gap is studied for conformers of the polyacenes ($C_{4n+2}H_{2n+4}$) up to n = 30. △ Less

Submitted 1 November, 2004; originally announced November 2004.

Comments: 10 pages, 5 eps figures submitted to Phys. Rev. A

Journal ref: Phys. Rev. A, vol. 71, art.no. 022508, 2005.