Source author record

Michele Pavanello

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

15 published item(s)

preprint2026arXiv

Pseudopotentials for Orbital-Free DFT: Capturing Nonlocality and Correcting Functional Approximants

Developing reliable pseudopotentials for orbital-free density functional theory (OF-DFT), especially for transition metals, remains a significant challenge. In this study, we provide a theoretical framework for analyzing pseudization strategies for OF-DFT calculations. From the analysis arises a proposed pseudization method which involves constructing local pseudopotentials by targeting existing Kohn-Sham DFT pseudopotentials through an optimized effective potential procedure. We produce four distinct sets of local pseudopotentials and evaluate their accuracy and transferability on the transition metal elements. Our results indicate a substantial improvement over previously available pseudopotentials. Although current OF-DFT functionals still only reach a qualitative accuracy for transition metals, our newly developed pseudopotentials provide a rigorous framework for further methodological advancements.

preprint2020arXiv

DFTpy: An efficient and object-oriented platform for orbital-free DFT simulations

In silico materials design is hampered by the computational complexity of Kohn-Sham DFT, which scales cubically with the system size. Owing to the development of new-generation kinetic energy density functionals (KEDFs), orbital-free DFT (OFDFT, a linear-scaling method) can now be successfully applied to a large class of semiconductors and such finite systems as quantum dots and metal clusters. In this work, we present DFTpy, an open source software implementing OFDFT written entirely in Python 3 and outsourcing the computationally expensive operations to third-party modules, such as NumPy and SciPy. When fast simulations are in order, DFTpy exploits the fast Fourier transforms (FFTs) from PyFFTW. New-generation, nonlocal and density-dependent-kernel KEDFs are made computationally efficient by employing linear splines and other methods for fast kernel builds. We showcase DFTpy by solving for the electronic structure of a million-atom system of aluminum metal which was computed on a single CPU. The Python 3 implementation is object-oriented, opening the door to easy implementation of new features. As an example, we present a time-dependent OFDFT implementation (hydrodynamic DFT) which we use to compute the spectra of small metal cluster recovering qualitatively the time-dependent Kohn-Sham DFT result. The Python code base allows for easy implementation of APIs. We showcase the combination of DFTpy and ASE for molecular dynamics simulations (NVT) of liquid metals. DFTpy is released under the MIT license.

preprint2019arXiv

Capturing Multireference Excited States by Constrained DFT

The computation of excited electronic states with commonly employed (approximate) methods is challenging, typically yielding states of lower quality than the corresponding ground state for a higher computational cost. In this work, we present a mean field method that extends the previously proposed eXcited Constrained DFT (XCDFT) from single Slater determinants to ensemble 1-RDMs for computing low-lying excited states. The method still retains an associated computational complexity comparable to a semilocal DFT calculation while at the same time is capable of approaching states with multireference character. We benchmark the quality of this method on well-established test sets, finding good descriptions of the electronic structure of multireference states and maintaining an overall accuracy for the predicted excitation energies comparable to semilocal TDDFT.

preprint2019arXiv

Nonlocal Subsystem Density Functional Theory

By invoking a divide-and-conquer strategy, subsystem DFT dramatically reduces the computational cost of large-scale, \textit{ab-initio} electronic structure simulations of molecules and materials. The central ingredient setting subsystem DFT apart from Kohn-Sham DFT is the non-additive kinetic energy functional (NAKE). Currently employed NAKEs are at most semilocal (i.e., they only depend on the electron density and its gradient), and as a result of this approximation, so far only systems composed of weakly interacting subsystems have been successfully tackled. In this work, we advance the state-of-the-art by introducing fully nonlocal NAKEs in subsystem DFT simulations for the first time. A benchmark analysis based on the S22-5 test set shows that nonlocal NAKEs considerably improve the computed interaction energies and electron density compared to commonly employed GGA NAKEs, especially when the inter-subsystem electron density overlap is high. Most importantly, we resolve the long standing problem of too attractive interaction energy curves typically resulting from the use of GGA NAKEs.

preprint2015arXiv

A critical look at methods for calculating charge transfer couplings fast and accurately

We present here a short and subjective review of methods for calculating charge transfer couplings. Although we mostly focus on Density Functional Theory, we discuss a small subset of semiempirical methods as well as the adiabatic-to-diabatic transformation methods typically coupled with wavefunction-based electronic structure calculations. In this work, we will present the reader with a critical assessment of the regimes that can be modeled by the various methods their strengths and weaknesses. In order to give a feeling about the practical aspects of the calculations, we also provide the reader with a practical protocol for running coupling calculations with the recently developed FDE-ET method.

preprint2015arXiv

Performance of Frozen Density Embedding for Modeling Hole Transfer Reactions

We have carried out a thorough benchmark of the FDE-ET method for calculating hole transfer couplings. We have considered 10 exchange-correlation functionals, 3 non-additive kinetic energy functionals and 3 basis sets. Overall, we conclude that with a 7% mean relative unsigned error, the PBE functional coupled with the PW91k non-additive Kinetic energy functional and a TZP basis set constitutes the most stable, and accurate level of theory for hole-transfer coupling calculations. The FDE-ET method is found to be an excellent tool for computing diabatic couplings for hole transfer reactions.

preprint2015arXiv

Subsystem Density-Functional Theory as an Effective Tool for Modeling Ground and Excited States, their Dynamics, and Many-Body Interactions

Subsystem Density-Functional Theory (DFT) is an emerging technique for calculating the electronic structure of complex molecular and condensed phase systems. In this topical review, we focus on some recent advances in this field related to the computation of condensed phase systems, their excited states, and the evaluation of many-body interactions between the subsystems. As subsystem DFT is in principle an exact theory, any advance in this field can have a dual role. One is the possible applicability of a resulting method in practical calculations. The other is the possibility of shedding light on some quantum-mechanical phenomenon which is more easily treated by subdividing a supersystem into subsystems. An example of the latter is many-body interactions. In the discussion, we present some recent work from our research group as well as some new results, casting them in the current state-of-the-art in this review as comprehensively as possible.

preprint2015arXiv

Subsystem real-time Time Dependent Density Functional Theory

We present the extension of Frozen Density Embedding (FDE) theory to real-time Time Dependent Density Functional Theory (rt-TDDFT). FDE a is DFT-in-DFT embedding method that allows to partition a larger Kohn-Sham system into a set of smaller, coupled Kohn-Sham systems. Additional to the computational advantage, FDE provides physical insight into the properties of embedded systems and the coupling interactions between them. The extension to rt-TDDFT is done straightforwardly by evolving the Kohn-Sham subsystems in time simultaneously, while updating the embedding potential between the systems at every time step. Two main applications are presented: the explicit excitation energy transfer in real time between subsystems is demonstrated for the case of the Na$_4$ cluster and the effect of the embedding on optical spectra of coupled chromophores. In particular, the importance of including the full dynamic response in the embedding potential is demonstrated.

preprint2015arXiv

Visible line intensities of the triatomic hydrogen ion from experiment and theory

The visible spectrum of H3+ is studied using high-sensitivity action spectroscopy in a cryogenic radiofrequency multipole trap. Advances are made to measure the weak ro-vibrational transitions from the lowest rotational states of H3+ up to high excitation energies providing visible line intensities and, after normalisation to an infrared calibration line, the corresponding Einstein $B$ coefficients. {\it Ab initio} predictions for the Einstein $B$ coefficients are obtained from a highly precise dipole moment surface of H3+ and found to be in excellent agreement, even in the region where states have been classified as chaotic.

preprint2014arXiv

FDE-vdW: A van der Waals Inclusive Subsystem Density-Functional Theory

We present a formally exact van der Waals inclusive electronic structure theory, called FDE-vdW, based on the Frozen Density Embedding formulation of subsystem Density-Functional Theory. In subsystem DFT, the energy functional is composed of subsystem additive and non-additive terms. We show that an appropriate definition of the long-range correlation energy is given by the value of the non-additive correlation functional. This functional is evaluated using the Fluctuation-Dissipation Theorem aided by a formally exact decomposition of the response functions into subsystem contributions. FDE-vdW is derived in detail and several approximate schemes are proposed, which lead to practical implementations of the method. We show that FDE-vdW is Casimir-Polder consistent, i.e. it reduces to the generalized Casimir-Polder formula for asymptotic inter-subsystems separations. Pilot calculations of binding energies of 13 weakly bound complexes singled out from the S22 set show a dramatic improvement upon semilocal subsystem DFT, provided that an appropriate exchange functional is employed. The convergence of FDE-vdW with basis set size is discussed, as well as its dependence on the choice of associated density functional approximant.

preprint2014arXiv

Periodic Subsystem Density-Functional Theory

By partitioning the electron density into subsystem contributions, the Frozen Density Embedding (FDE) formulation of subsystem DFT has recently emerged as a powerful tool for reducing the computational scaling of Kohn--Sham DFT. To date, however, FDE has been employed to molecular systems only. Periodic systems, such as metals, semiconductors, and other crystalline solids have been outside the applicability of FDE, mostly because of the lack of a periodic FDE implementation. To fill this gap, in this work we aim at extening FDE to treat subsystems of molecular and periodic character. This goal is achieved by a dual approach. On one side, the development of a theoretical framework for periodic subsystem DFT. On the other, the realization of the method into a parallel computer code. We find that periodic FDE is capable of reproducing total electron densities and (to a lesser extent) also interaction energies of molecular systems weakly interacting with metallic surfaces. In the pilot calculations considered, we find that FDE fails in those cases where there is appreciable density overlap between the subsystems. Conversely, we find FDE to be in semiquantitative agreement (but still within chemical accuracy) with Kohn--Sham DFT when the inter-subsystem density overlap is low. We also conclude that to make FDE a suitable method for describing molecular adsorption at surfaces, kinetic energy density functionals that go beyond the GGA level must be employed.

preprint2014arXiv

Quantifying Environmental Effects on the Decay of Hole Transfer Couplings in Biosystems

In the past two decades, many research groups worldwide have tried to understand and categorize simple regimes in the charge transfer of such biological systems as DNA. Theoretically speaking, the lack of exact theories for electron-nuclear dynamics on one side, and poor quality of the parameters needed by model Hamiltonians and nonadiabatic dynamics alike (such as couplings and site energies) on the other, are the two main difficulties for an appropriate description of the charge transfer phenomena. In this work, we present an application of a previously benchmarked and linear-scaling subsystem DFT method for the calculation of couplings, site energies and superexchange decay factors ($β$) of several biological donor-acceptor dyads, as well as double stranded DNA oligomers comprised of up to 5 base pairs. The calculations are all-electron, and provide a clear view of the role of the environment on superexchange couplings in DNA - they follow experimental trends and confirm previous semiempirical calculations. The subsystem DFT method is proven to be an excellent tool for long-range, bridge-mediated coupling and site energy calculations of embedded molecular systems.

preprint2013arXiv

An Accurate and Linear Scaling Method to Calculate Charge-Transfer Excitation Energies and Diabatic Couplings

Quantum--Mechanical methods that are both computationally fast and accurate are not yet available for electronic excitations having charge transfer character. In this work, we present a significant step forward towards this goal for those charge transfer excitations that take place between non-covalently bound molecules. In particular, we present a method that scales linearly with the number of non-covalently bound molecules in the system and is based on a two-pronged approach: The molecular electronic structure of broken-symmetry charge-localized states is obtained with the Frozen Density Embedding formulation of subsystem Density-Functional Theory; subsequently, in a post-SCF calculation, the full-electron Hamiltonian and overlap matrix elements among the charge-localized states are evaluated with an algorithm which takes full advantage of the subsystem DFT density partitioning technique. The method is benchmarked against Coupled-Cluster calculations and achieves chemical accuracy for the systems considered for intermolecular separations ranging from hydrogen-bond distances to tents of Ångstroms. Numerical examples are provided for molecular clusters comprised of up to 56 non-covalently bound molecules.

preprint2013arXiv

Calculating Hyperfine Couplings in Large Ionic Crystals Containing Hundreds of QM Atoms: Subsystem DFT is the Key

We present an application of the linear scaling Frozen Density Embedding (FDE) formulation of subsystem DFT to the calculation of isotropic hyperfine coupling constants (hfccs) of atoms belonging to a guanine radical cation embedded in a guanine hydrochloride monohydrate crystal. The model systems considered range from an isolated guanine to a 15,000 atom QM/MM cluster where the QM region is comprised of 36 protonated guanine cations, 36 chlorine anions and 42 water molecules. Our calculations show that the embedding effects of the surrounding crystal cannot be reproduced neither by small model systems nor by a pure QM/MM procedure. Instead, a large QM region is needed to fully capture the complicated nature of the embedding effects in this system. The unprecedented system size for a relativistic all-electron isotropic hfccs calculation can be approached in this work because the local nature of the electronic structure of the organic crystals considered is fully captured by the FDE approach.

preprint2013arXiv

On the Subsystem Formulation of Linear-Response Time-Dependent DFT

A new and thorough derivation of linear-response subsystem time-dependent density functional theory (TD-DFT) is presented and analyzed in detail. Two equivalent derivations are presented and naturally yield self-consistent subsystem TD-DFT equations. One reduces to the subsystem TD-DFT formalism of Neugebauer [J. Chem. Phys. 126, 134116 (2007)10.1063/1.2713754]. The other yields Dyson type equations involving three types of subsystem response functions: coupled, uncoupled, and Kohn-Sham. The Dyson type equations for subsystem TD-DFT are derived here for the first time. The response function formalism reveals previously hidden qualities and complications of subsystem TD-DFT compared with the regular TD-DFT of the supersystem. For example, analysis of the pole structure of the subsystem response functions shows that each function contains information about the electronic spectrum of the entire supersystem. In addition, comparison of the subsystem and supersystem response functions shows that, while the correlated response is subsystem additive, the Kohn-Sham response is not. Comparison with the non-subjective partition DFT theory shows that this non-additivity is largely an artifact introduced by the subjective nature of the density partitioning in subsystem DFT.