Source author record

Andreas Savin

Andreas Savin 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

25works
8topics
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

25 published item(s)

preprint2022arXiv

The effect of uncertainty on building blocks in molecules

Probabilities to find a chosen number of electrons in flexible domains of space are calculated for highly correlated wave functions. Quantum mechanics can produce higher probabilities for chemically relevant arrangements of electrons in these regions. However, the probability to have a given arrangement, e.g., that corresponding to chemical formulas (bonds or atoms), is low although being often maximal. Like in valence bond theory, it is useful to consider alternative distributions of the electrons. Exchanges of electrons should be considered not only between atoms, but also between other types of regions, like those attributed to lone pairs. It is useful to have definitions flexible enough to allow the user to find the reference representations he considers most relevant. We tentatively suggest a tool (the effective number of parties) to help one make the choice.

preprint2022arXiv

Two-particle coalescence conditions revisited

The notion of the n-th order local energy, generated by the n-th power of the Hamiltonian, has been introduced. The n-th order two-particle coalescence conditions have been derived from the requirements that the n-th order local energy at the coalescence point is non-singular and equal to the n-th power of the Hamiltonian eigenvalue. The first condition leads to energy-independent constraints. The second one is state-specific. The analysis has been done using a radial, one-dimensional, model Hamiltonian. The model is valid in the asymptotic region of r ~ 0. The coalescence conditions set the relations between the expansion coefficients of the radial wave function into a power series with respect to r.

preprint2022arXiv

What is the number of electrons in a spatial domain?

We like to attribute a number of electrons to spatial domains (atoms, bonds, ...). However, as a rule, the number of electrons in a spatial domain is not a sharp number. We thus study probabilities for having any number of electrons (between 0 and the total number of electrons in the system) in a given spatial domain. We show that by choosing a domain that maximizes a chosen probability (or is close to it), one obtains higher probabilities for chemically relevant regions. The probability to have a given electronic arrangement, - for example, by attributing a number of electrons to an atomic shell - can be low. It remains so even in the "best" case, i.e, if the spatial domain is chosen to maximize the chosen probability. In other words, the number of electrons in a spatial region significantly fluctuates. The freedom of choosing the number of electrons we are interested in shows that a "chemical" question is not always well-posed. We show it using as an example the KrF2 molecule.

preprint2020arXiv

Acknowledging user requirements for accuracy in computational chemistry benchmarks

Computational chemistry has become an important complement to experimental measurements. In order to choose among the multitude of the existing approximations, it is common to use benchmark data sets, and to issue recommendations based on numbers such as mean absolute errors. We argue, using as an example band gaps calculated with density functional approximations, that a more careful study of the benchmark data is needed, stressing that the user's requirements play a role in the choice of an appropriate method. We also appeal to those who measure data capable of being used as a reference, to publish error estimates. We show how the latter can affect the judgment of approximations used in computational chemistry.

preprint2020arXiv

Probabilistic performance estimators for computational chemistry methods: Systematic Improvement Probability and Ranking Probability Matrix. II. Applications

In the first part of this study (Paper I), we introduced the systematic improvement probability (SIP) as a tool to assess the level of improvement on absolute errors to be expected when switching between two computational chemistry methods. We developed also two indicators based on robust statistics to address the uncertainty of ranking in computational chemistry benchmarks: Pinv , the inversion probability between two values of a statistic, and Pr , the ranking probability matrix. In this second part, these indicators are applied to nine data sets extracted from the recent benchmarking literature. We illustrate also how the correlation between the error sets might contain useful information on the benchmark dataset quality, notably when experimental data are used as reference.

preprint2020arXiv

Using the Gini coefficient to characterize the shape of computational chemistry error distributions

The distribution of errors is a central object in the assesment and benchmarking of computational chemistry methods. The popular and often blind use of the mean unsigned error as a benchmarking statistic leads to ignore distributions features that impact the reliability of the tested methods. We explore how the Gini coefficient offers a global representation of the errors distribution, but, except for extreme values, does not enable an unambiguous diagnostic. We propose to relieve the ambiguity by applying the Gini coefficient to mode-centered error distributions. This version can usefully complement benchmarking statistics and alert on error sets with potentially problematic shapes.

preprint2016arXiv

Exchange-correlation functionals via local interpolation along the adiabatic connection

The construction of density-functional approximations is explored by modeling the adiabatic connection em locally, using energy densities defined in terms of the electrostatic potential of the exchange-correlation hole. These local models are more amenable to the construction of size-consistent approximations than their global counterparts. In this work we use accurate input local ingredients to assess the accuracy of range of local interpolation models against accurate exchange-correlation energy densities. The importance of the strictly-correlated electrons (SCE) functional describing the strong coupling limit is emphasized, enabling the corresponding interpolated functionals to treat strong correlation effects. In addition to exploring the performance of such models numerically for the helium and beryllium isoelectronic series and the dissociation of the hydrogen molecule, an approximate analytic model is presented for the initial slope of the local adiabatic connection. Comparisons are made with approaches based on global models and prospects for future approximations based on the local adiabatic connection are discussed.

preprint2016arXiv

Smooth models for the Coulomb potential

Smooth model potentials with parameters selected to reproduce the spectrum of one-electron atoms are used to approximate the singular Coulomb potential. Even when the potentials do not mimic the Coulomb singularity, much of the spectrum is reproduced within the chemical accuracy. For the Hydrogen atom, the smooth approximations to the Coulomb potential are more accurate for higher angular momentum states. The transferability of the model potentials from an attractive interaction (Hydrogen atom) to a repulsive one (Harmonium and the uniform electron gas) is discussed.

preprint2016arXiv

Spooky correlations and unusual van der Waals forces between gapless and near-gapless molecules

We consider the zero-temperature van der Waals interaction between two molecules, each of which has a zero or near-zero electronic gap between a groundstate and the first excited state, using a toy model molecule ( equilateral H3) as an example. We show that the van der Waals energy between two groundstate molecules falls off as D^(-3) instead of the usual D^(-6) dependence, when the molecules are separated by distance D: We show that this is caused by perfect "spooky" correlation between the two fluctuating electric dipoles. The phenomenon is related to, but not the same as, the "resonant" interaction between an electronically excited and a groundstate molecule introduced by Eisenschitz and London in 1930. It is also an example of "type C van der Waals non- additivity" recently introduced by one of us ( Int. J. Quantum Chem. 114, 1157 (2014)). Our toy molecule H3 is not stable, but symmetry considerations suggest that a similar vdW phenomenon may be observable, despite Jahn-Teller effects, in molecules with discrete rotational symmetry and broken inversion symmetry, such as certain metal atom clusters. The motion of the nuclei will need to be included for a definitive analysis of such cases, however.

preprint2015arXiv

Calculating excitation energies by extrapolation along adiabatic connections

In this paper, an alternative method to range-separated linear-response time-dependent density-functional theory and perturbation theory is proposed to improve the estimation of the energies of a physical system from the energies of a partially interacting system. Starting from the analysis of the Taylor expansion of the energies of the partially interacting system around the physical system, we use an extrapolation scheme to improve the estimation of the energies of the physical system at an intermediate point of the range-separated or linear adiabatic connection where either the electron--electron interaction is scaled or only the long-range part of the Coulomb interaction is included. The extrapolation scheme is first applied to the range-separated energies of the helium and beryllium atoms and of the hydrogen molecule at its equilibrium and stretched geometries. It improves significantly the convergence rate of the energies toward their exact limit with respect to the range-separation parameter. The range-separated extrapolation scheme is compared with a similar approach for the linear adiabatic connection, highlighting the relative strengths and weaknesses of each approach.

preprint2015arXiv

Maximum probability domains for the analysis of the microscopic structure of liquids

We introduce the concept of maximum probability domains (MPDs), developed in the context of the analysis of electronic densities, in the study of the microscopic spatial structures of liquids. The idea of locating a particle in a three dimensional region, by determining the domain where the probability of finding that, and only that, particle is maximum, gives an interesting characterization of the local structure of the liquid. The optimization procedure, required for the search of the domain of maximum probability, is carried out by the implementation of the level set method. Results for a couple of case studies are presented, to illustrate the structure of liquid water at ambient conditions and upon increasing pressure from the point of view of MPDs and to compare the information encoded in the solvation shells of sodium in water with, once again, that extracted from the MPDs.

preprint2014arXiv

Excitation energies along a range-separated adiabatic connection

We present a study of the variation of total energies and excitationenergies along a range-separated adiabatic connection. This connectionlinks the non-interacting Kohn-Sham electronic system to the physicalinteracting system by progressively switching on theelectron-electron interactions whilst simultaneously adjusting aone-electron effective potential so as to keep the ground-statedensity constant. The interactions are introduced in arange-dependent manner, first introducing predominantly long-range,and then all-range, interactions as the physical system is approached,as opposed to the conventional adiabatic connection where theinteractions are introduced by globally scaling the standard Coulomb interaction.Reference data are reported for the He and Be atoms and the H2molecule, obtained by calculating the short-range effective potentialat the full configuration-interaction level using Lieb'sLegendre-transform approach. As the strength of the electron-electroninteractions increases, the excitation energies, calculated for thepartially interacting systems along the adiabatic connection, offerincreasingly accurate approximations to the exact excitation energies.Importantly, the excitation energies calculated at an intermediatepoint of the adiabatic connection are much better approximations tothe exact excitation energies than are the corresponding Kohn-Shamexcitation energies. This is particularly evident in situationsinvolving strong static correlation effects and states with multipleexcitation character, such as the dissociating H2 molecule. Theseresults highlight the utility of long-range interacting referencesystems as a starting point for the calculation of excitation energiesand are of interest for developing and analyzing practical approximaterange-separated density-functional methodologies.

preprint2014arXiv

Excited states from range-separated density-functional perturbation theory

We explore the possibility of calculating electronic excited states by using perturbation theory along a range-separated adiabatic connection. Starting from the energies of a partially interacting Hamiltonian, a first-order correction is defined with two variants of perturbation theory: a straight-forward perturbation theory, and an extension of the G{ö}rling--Levy one that has the advantage of keeping the ground-state density constant at each order in the perturbation. Only the first, simpler, variant is tested here on the helium and beryllium atoms and on the dihydrogene molecule. The first-order correction within this perturbation theory improves significantly the total ground-and excited-state energies of the different systems. However, the excitation energies are mostly deterio-rated with respect to the zeroth-order ones, which may be explained by the fact that the ionization energy is no longer correct for all interaction strengths. The second variant of the perturbation theory should improve these results but has not been tested yet along the range-separated adiabatic connection.

preprint2014arXiv

Range-separated density-functional theory with random phase approximation applied to noncovalent intermolecular interactions

Range-separated methods combining a short-range density functional with long-range random phase approximations RPAs with or without exchange response kernel are tested on rare-gas dimers and the S22 benchmark set of weakly interacting complexes of Jurecka et al. [Phys. Chem. Chem. Phys. 8, 1985 (2006 )] . The methods are also compared to full-range RPA approaches. Both range separation and inclusion of the Hartree-Fock exchange kernel largely improve the accuracy of intermolecular interaction energies. The best results are obtained with the method called RSH+RPAx, which yields interaction energies for the S22 set with an estimated mean absolute error of about 0.5-0.6 kcal/mol, corresponding to a mean absolute percentage error of about 7%-9% depending on the reference interaction energies used. In particular, the RSH+RPAx method is found to be overall more accurate than the range-separated method based on long-range second-order Møller-Plesset MP2 perturbation theory RSH+MP2 .

preprint2013arXiv

Electronic excitations from a linear-response range-separated hybrid scheme

We study linear-response time-dependent density-functional theory (DFT) based on the single-determinant range-separated hybrid (RSH) scheme, i.e. combining a long-range Hartree-Fock exchange kernel with a short-range DFT exchange-correlation kernel, for calculating electronic excitation energies of molecular systems. It is an alternative to the long-range correction (LC) scheme which has a standard full-range DFT correlation kernel instead of only a short-range one. We discuss the local-density approximation (LDA) to the short-range exchange-correlation kernel, and assess the performance of the linear-response RSH scheme for singlet-singlet and singlet-triplet valence and Rydberg excitations in the N2, CO, H2CO, C2H4, and C6H6 molecules, and for the first charge-transfer excitation in the C2H4-C2F4 dimer. The introduction of long-range HF exchange corrects the underestimation of charge-transfer and high-lying Rydberg excitation energies obtained with standard (semi)local density-functional approximations, but also leads to underestimated excitation energies to low-lying spin-triplet valence states which can be cured by the Tamm-Dancoff approximation. This work thus suggests that the present linear-response RSH scheme is a reasonable starting approximation for describing electronic excitation energies, even before adding an explicit treatment of long-range correlation.

preprint2012arXiv

A multiconfigurational hybrid density-functional theory

We propose a multiconfigurational hybrid density-functional theory which rigorously combines a multiconfiguration self-consistent-field calculation with a density-functional approximation based on a linear decomposition of the electron-electron interaction. This gives a straightforward extension of the usual hybrid approximations by essentially adding a fraction λof exact static correlation in addition to the fraction λof exact exchange. Test calculations on the cycloaddition reactions of ozone with ethylene or acetylene and the dissociation of diatomic molecules with the Perdew-Burke-Ernzerhof (PBE) and Becke-Lee-Yang-Parr (BLYP) density functionals show that a good value of λis 0.25, as in the usual hybrid approximations. The results suggest that the proposed multiconfigurational hybrid approximations can improve over usual density-functional calculations for situations with strong static correlation effects.

preprint2011arXiv

Closed-shell ring coupled cluster doubles theory with range separation applied on weak intermolecular interactions

We explore different variants of the random phase approximation (RPA) to the correlation energy derived from closed-shell ring-diagram approximations to coupled cluster doubles theory. We implement these variants in range-separated density-functional theory, i.e. by combining the long-range random phase approximations with short-range density-functional approximations. We perform tests on the rare-gas dimers He_2, Ne_2, and Ar_2, and on the weakly interacting molecular complexes of the S22 set of Jurecka et al. [Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The two best variants correspond to the ones originally proposed by Szabo and Ostlund [J. Chem. Phys. 67, 4351 (1977)]. With range separation, they reach mean absolute errors on the equilibrium interaction energies of the S22 set of about 0.4 kcal/mol, corresponding to mean absolute percentage errors of about 4 %, with the aug-cc-pVDZ basis set.

preprint2011arXiv

Double-hybrid density-functional theory made rigorous

We provide a rigorous derivation of a class of double-hybrid approximations, combining Hartree-Fock exchange and second-order Moller-Plesset correlation with a semilocal exchange-correlation density functional. These double-hybrid approximations contain only one empirical parameter and use a density-scaled correlation energy functional. Neglecting density scaling leads to a one-parameter version of the standard double-hybrid approximations. We assess the performance of these double-hybrid schemes on representative test sets of atomization energies and reaction barrier heights, and we compare to other hybrid approximations, including range-separated hybrids. Our best one-parameter double-hybrid approximation, called 1DH-BLYP, roughly reproduces the two parameters of the standard B2-PLYP or B2GP-PLYP double-hybrid approximations, which shows that these methods are not only empirically close to an optimum for general chemical applications but are also theoretically supported.

preprint2010arXiv

Range-separated density-functional theory with random phase approximation: detailed formalism and illustrative applications

Using Green-function many-body theory, we present the details of a formally exact adiabatic-connection fluctuation-dissipation density-functional theory based on range separation, which was sketched in Toulouse, Gerber, Jansen, Savin and Angyan, Phys. Rev. Lett. 102, 096404 (2009). Range-separated density-functional theory approaches combining short-range density functional approximations with long-range random phase approximations (RPA) are then obtained as well-identified approximations on the long-range Green-function self-energy. Range-separated RPA-type schemes with or without long-range Hartree-Fock exchange response kernel are assessed on rare-gas and alkaline-earth dimers, and compared to range-separated second-order perturbation theory and range-separated coupled-cluster theory.

preprint2010arXiv

Symmetry Breaking of Relativistic Multiconfiguration Methods in the Nonrelativistic Limit

The multiconfiguration Dirac-Fock method allows to calculate the state of relativistic electrons in atoms or molecules. This method has been known for a long time to provide certain wrong predictions in the nonrelativistic limit. We study in full mathematical details the nonlinear model obtained in the nonrelativistic limit for Be-like atoms. We show that the method with sp+pd configurations in the J=1 sector leads to a symmetry breaking phenomenon in the sense that the ground state is never an eigenvector of L^2 or S^2. We thereby complement and clarify some previous studies.

preprint2009arXiv

Full Configuration Interaction wave function as a formal solution to the Optimized Effective Potential and Kohn-Sham models in finite basis sets

Using finite basis sets, it is shown how to construct a local Hamiltonian, such that one of its infinitely many degenerate eigenfunctions is the ground state full configuration interaction (FCI) wave function in that basis set. Formally, the local potential of this Hamiltonian is the optimized effective potential and the exact Kohn-Sham potential at the same time, because the FCI wave function yields the exact ground-state density and energy. It is not the aim of this paper to provide a new algorithm for obtaining FCI wave functions. A new potential is the goal.

preprint2008arXiv

Range separation combined with the Overhauser model: Application to the H$_2$ molecule along the dissociation curve

The combination of density-functional theory with other approaches to the many-electron problem through the separation of the electron-electron interaction into a short-range and a long-range contribution (range separation) is a successful strategy, which is raising more and more interest in recent years. We focus here on a range-separated method in which only the short-range correlation energy needs to be approximated, and we model it within the "extended Overhauser approach". We consider the paradigmatic case of the H$_2$ molecule along the dissociation curve, finding encouraging results. By means of very accurate variational wavefunctions, we also study how the effective electron-electron interaction appearing in the Overhauser model should be in order to yield the exact correlation energy for standard Kohn-Sham density functional theory.

preprint2008arXiv

Study of the discontinuity of the exchange-correlation potential in an exactly soluble case

It was found by Perdew, Parr, Levy, and Balduz [Phys. Rev. Lett. {\bf 49}, 1691 (1982)] and by Sham and Schlüter [Phys. Rev. Lett. {\bf 51}, 1884 (1983)] that the exact Kohn-Sham exchange-correlation potential of an open system may jump discontinuosly as the particle number crosses an integer, with important physical consequences. Recently, Sagvolden and Perdew [Phys. Rev. A {\bf 77}, 012517 (2008)] have analyzed the discontinuity of the exchange-correlation potential as the particle number crosses one, with an illustration that uses a model density for the H$^-$ ion. In this work, we extend their analysis to the case in which the external potential is the simple harmonic confinement, choosing spring-constant values for which the two-electron hamiltonian has an analytic solution. This way, we can obtain the exact, analytic, exchange and correlation potentials for particle number fluctuating between zero and two, illustrating the discontinuity as the particle number crosses one without introducing any model or approximation. We also discuss exchange and correlation separately.

preprint2006arXiv

Scaling relations, virial theorem and energy densities for long-range and short-range density functionals

We analyze a decomposition of the Coulomb electron-electron interaction into a long-range and a short-range part in the framework of density functional theory, deriving some scaling relations and the corresponding virial theorem. We study the behavior of the local density approximation in the high-density limit for the long-range and the short-range functionals by carrying out a detailed analysis of the correlation energy of a uniform electron gas interacting via a long-range only electron-electron repulsion. Possible definitions of exchange and correlation energy densities are discussed and clarified with some examples.

preprint2006arXiv

Short-range exchange-correlation energy of a uniform electron gas with modified electron-electron interaction

We calculate the short-range exchange-correlation energy of the uniform electron gas with two modified electron-electron interactions. While the short-range exchange functionals are calculated analytically, Coupled-Cluster and Fermi-hypernetted chain calculations are carried out for the correlation energy and the results are fitted to an analytical parametrization. These data enable to construct the local density approximation corresponding to these modified interactions.