Comment on Inverse problem from discrete spectrum in the D = 2 dimensional space (2014 Phys. Scr. 989 085201)

We analyse a method for the construction of the potential-energy function from the moments of the ground-state density. The sum rule on which some expressions are based appear to be wrong, as well as the moments and potential-energy functions derived for some illustrative examples.

Post-processing of the planewave approximation of Schrödinger equations. Part I: linear operators

In this article we prove a priori error estimates for the perturbation-based post-processing of the plane-wave approximation of Schrödinger equations introduced and tested numerically in previous works (Cancès, Dusson, Maday, Stamm and Vohralík, (2014), A perturbation-method-based a posteriori estimator for the planewave discretization of nonlinear Schrödinger equations. C. R. Math., 352, 941--946; Cancès, Dusson, Maday, Stamm and Vohralík, (2016), A perturbation-method-based postprocessing for the planewave discretization of Kohn–Sham models. J. Comput. Phys., 307, 446--459.) We consider here a Schrödinger operator ${{\mathscr{H}} \,}= -\frac{1}{2}\varDelta +{\mathscr{V}}$ on $L^2(\varOmega )$, where $\varOmega $ is a cubic box with periodic boundary conditions and where ${\mathscr{V}}$ is a multiplicative operator by a regular-enough function ${\mathscr{V}}$. The quantities of interest are, on the one hand, the ground-state energy defined as the sum of the lowest $N$ eigenvalues of ${{\mathscr{H}} \,}$, and, on the other hand, the ground-state density matrix that is the spectral projector on the vector space spanned by the associated eigenvectors. Such a problem is central in first-principle molecular simulation, since it corresponds to the so-called linear subproblem in Kohn–Sham density functional theory. Interpreting the exact eigenpairs of ${{\mathscr{H}} \,}$ as perturbations of the numerical eigenpairs obtained by a variational approximation in a plane-wave (i.e., Fourier) basis we compute first-order corrections for the eigenfunctions, which are turned into corrections on the ground-state density matrix. This allows us to increase the accuracy by a factor proportional to the inverse of the kinetic energy cutoff ${E_{\textrm{c}}}^{-1}$ of both the ground-state energy and the ground-state density matrix in Hilbert–Schmidt norm at a low computational extra cost. Indeed, the computation of the corrections only requires the computation of the residual of the solution in a larger plane-wave basis and two fast Fourier transforms per eigenvalue.

Moments of the ground state density for the d-dimensional Fermi gas in an harmonic trap

We consider properties of the ground state density for the [Formula: see text]-dimensional Fermi gas in an harmonic trap. Previous work has shown that the [Formula: see text]-dimensional Fourier transform has a very simple functional form. It is shown that this fact can be used to deduce that the density itself satisfies a third-order linear differential equation, previously known in the literature but from other considerations. It is shown too how this implies a closed form expression for the [Formula: see text]th non-negative integer moments of the density, and a second-order recurrence. Both can be extended to general Re[Formula: see text]. The moments, and the smoothed density, permit expansions in [Formula: see text], where [Formula: see text], with [Formula: see text] denoting the shell label. The moment expansion substituted in the second-order recurrence gives a generalization of the Harer–Zagier recurrence, satisfied by the coefficients of the [Formula: see text] expansion of the moments of the spectral density for the Gaussian unitary ensemble in random matrix theory.

Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules

Prompted by recent reports of large errors in noncovalent interaction (NI) energies obtained from many-body perturbation theory (MBPT), we compare the performance of second-order Møller–Plesset MBPT (MP2), spin-scaled MP2, dispersion-corrected semilocal density functional approximations (DFA), and the post-Kohn–Sham random phase approximation (RPA) for predicting binding energies of supramolecular complexes contained in the S66, L7, and S30L benchmarks. All binding energies are extrapolated to the basis set limit, corrected for basis set superposition errors, and compared to reference results of the domain-based local pair-natural orbital coupled-cluster (DLPNO-CCSD(T)) or better quality. Our results confirm that MP2 severely overestimates binding energies of large complexes, producing relative errors of over 100% for several benchmark compounds. RPA relative errors consistently range between 5-10%, significantly less than reported previously using smaller basis sets, whereas spin-scaled MP2 methods show limitations similar to MP2, albeit less pronounced, and empirically dispersion-corrected DFAs perform almost as well as RPA. Regression analysis reveals a systematic increase of relative MP2 binding energy errors with the system size at a rate of approximately 0.1% per valence electron, whereas the RPA and dispersion-corrected DFA relative errors are virtually independent of the system size. These observations are corroborated by a comparison of computed rotational constants of organic molecules to gas-phase spectroscopy data contained in the ROT34 benchmark. To analyze these results, an asymptotic adiabatic connection symmetry-adapted perturbation theory (AC-SAPT) is developed which uses monomers at full coupling whose ground-state density is constrained to the ground-state density of the complex. Using the fluctuation–dissipation theorem, we obtain a nonperturbative “screened second-order” expression for the dispersion energy in terms of monomer quantities which is exact for non-overlapping subsystems and free of induction terms; a first-order RPA-like approximation to the Hartree, exchange, and correlation kernel recovers the macroscopic Lifshitz limit. The AC-SAPT expansion of the interaction energy is obtained from Taylor expansion of the coupling strength integrand. Explicit expressions for the convergence radius of the AC-SAPT series are derived within RPA and MBPT and numerically evaluated. Whereas the AC-SAPT expansion is always convergent for nondegenerate monomers when RPA is used, it is found to spuriously diverge for second-order MBPT, except for the smallest and least polarizable monomers. The divergence of the AC-SAPT series within MBPT is numerically confirmed within RPA; prior numerical results on the convergence of the SAPT expansion for MBPT methods are revisited and support this conclusion once sufficiently high orders are included. The cause of the failure of MBPT methods for NIs of large systems is missing or incomplete “electrodynamic” screening of the Coulomb interaction due to induced particle–hole pairs between electrons in different monomers, leaving the effective interaction too strong for AC-SAPT to converge. Hence, MBPT cannot be considered reliable for quantitative predictions of NIs, even in moderately polarizable molecules with a few tens of atoms. The failure to accurately account for electrodynamic polarization makes MBPT qualitatively unsuitable for applications such as NIs of nanostructures, macromolecules, and soft materials; more robust non-perturbative approaches such as RPA or coupled cluster methods should be used instead whenever possible.<br>

Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules

Prompted by recent reports of large errors in noncovalent interaction (NI) energies obtained from many-body perturbation theory (MBPT), we compare the performance of second-order Møller–Plesset MBPT (MP2), spin-scaled MP2, dispersion-corrected semilocal density functional approximations (DFA), and the post-Kohn–Sham random phase approximation (RPA) for predicting binding energies of supramolecular complexes contained in the S66, L7, and S30L benchmarks. All binding energies are extrapolated to the basis set limit, corrected for basis set superposition errors, and compared to reference results of the domain-based local pair-natural orbital coupled-cluster (DLPNO-CCSD(T)) or better quality. Our results confirm that MP2 severely overestimates binding energies of large complexes, producing relative errors of over 100% for several benchmark compounds. RPA relative errors consistently range between 5-10%, significantly less than reported previously using smaller basis sets, whereas spin-scaled MP2 methods show limitations similar to MP2, albeit less pronounced, and empirically dispersion-corrected DFAs perform almost as well as RPA. Regression analysis reveals a systematic increase of relative MP2 binding energy errors with the system size at a rate of approximately 0.1% per valence electron, whereas the RPA and dispersion-corrected DFA relative errors are virtually independent of the system size. These observations are corroborated by a comparison of computed rotational constants of organic molecules to gas-phase spectroscopy data contained in the ROT34 benchmark. To analyze these results, an asymptotic adiabatic connection symmetry-adapted perturbation theory (AC-SAPT) is developed which uses monomers at full coupling whose ground-state density is constrained to the ground-state density of the complex. Using the fluctuation–dissipation theorem, we obtain a nonperturbative “screened second-order” expression for the dispersion energy in terms of monomer quantities which is exact for non-overlapping subsystems and free of induction terms; a first-order RPA-like approximation to the Hartree, exchange, and correlation kernel recovers the macroscopic Lifshitz limit. The AC-SAPT expansion of the interaction energy is obtained from Taylor expansion of the coupling strength integrand. Explicit expressions for the convergence radius of the AC-SAPT series are derived within RPA and MBPT and numerically evaluated. Whereas the AC-SAPT expansion is always convergent for nondegenerate monomers when RPA is used, it is found to spuriously diverge for second-order MBPT, except for the smallest and least polarizable monomers. The divergence of the AC-SAPT series within MBPT is numerically confirmed within RPA; prior numerical results on the convergence of the SAPT expansion for MBPT methods are revisited and support this conclusion once sufficiently high orders are included. The cause of the failure of MBPT methods for NIs of large systems is missing or incomplete “electrodynamic” screening of the Coulomb interaction due to induced particle–hole pairs between electrons in different monomers, leaving the effective interaction too strong for AC-SAPT to converge. Hence, MBPT cannot be considered reliable for quantitative predictions of NIs, even in moderately polarizable molecules with a few tens of atoms. The failure to accurately account for electrodynamic polarization makes MBPT qualitatively unsuitable for applications such as NIs of nanostructures, macromolecules, and soft materials; more robust non-perturbative approaches such as RPA or coupled cluster methods should be used instead whenever possible.<br>

Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules

Prompted by recent reports of large errors in noncovalent interaction (NI) energies obtained from many-body perturbation theory (MBPT), we compare the performance of second-order Møller–Plesset MBPT (MP2), spin-scaled MP2, dispersion-corrected semilocal density functional approximations (DFA), and the post-Kohn–Sham random phase approximation (RPA) for predicting binding energies of supramolecular complexes contained in the S66, L7, and S30L benchmarks. All binding energies are extrapolated to the basis set limit, corrected for basis set superposition errors, and compared to reference results of the domain-based local pair-natural orbital coupled-cluster (DLPNO-CCSD(T)) or better quality. Our results confirm that MP2 severely overestimates binding energies of large complexes, producing relative errors of over 100% for several benchmark compounds. RPA relative errors consistently range between 5-10%, significantly less than reported previously using smaller basis sets, whereas spin-scaled MP2 methods show limitations similar to MP2, albeit less pronounced, and empirically dispersion-corrected DFAs perform almost as well as RPA. Regression analysis reveals a systematic increase of relative MP2 binding energy errors with the system size at a rate of approximately 1‰ per valence electron, whereas the RPA and dispersion-corrected DFA relative errors are virtually independent of the system size. These observations are corroborated by a comparison of computed rotational constants of organic molecules to gas-phase spectroscopy data contained in the ROT34 benchmark. To analyze these results, an asymptotic adiabatic connection symmetry-adapted perturbation theory (AC-SAPT) is developed which uses monomers at full coupling whose ground-state density is constrained to the ground-state density of the complex. Using the fluctuation–dissipation theorem, we obtain a nonperturbative “screened second-order” expression for the dispersion energy in terms of monomer quantities which is exact for non-overlapping subsystems and free of induction terms; a first-order RPA-like approximation to the Hartree, exchange, and correlation kernel recovers the macroscopic Lifshitz limit. The AC-SAPT expansion of the interaction energy is obtained from Taylor expansion of the coupling strength integrand. Explicit expressions for the convergence radius of the AC-SAPT series are derived within RPA and MBPT and numerically evaluated. Whereas the AC-SAPT expansion is always convergent for nondegenerate monomers when RPA is used, it is found to spuriously diverge for second-order MBPT, except for the smallest and least polarizable monomers. The divergence of the AC-SAPT series within MBPT is numerically confirmed within RPA; prior numerical results on the convergence of the SAPT expansion for MBPT methods are revisited and support this conclusion once sufficiently high orders are included. The cause of the failure of MBPT methods for NIs of large systems is missing or incomplete “electrodynamic” screening of the Coulomb interaction due to induced particle–hole pairs between electrons in different monomers, leaving the effective interaction too strong for AC-SAPT to converge. Hence, MBPT cannot be considered reliable for quantitative predictions of NIs, even in moderately polarizable molecules with a few tens of atoms. The failure to accurately account for electrodynamic polarization makes MBPT qualitatively unsuitable for applications such as NIs of nanostructures, macromolecules, and soft materials; more robust non-perturbative approaches such as RPA or coupled cluster methods should be used instead whenever possible.<br>

Exploring Baryon Rich Matter with Heavy-Ion Collisions

Collisions of heavy nuclei at (ultra-)relativistic energies provide a fascinating opportunity to re-create various forms of matter in the laboratory. For a short extent of time (10-22 s), matter under extreme conditions of temperature and density can exist. In dedicated experiments, one explores the microscopic structure of strongly interacting matter and its phase diagram. In heavy-ion reactions at SIS18 collision energies, matter is substantially compressed (2–3 times ground-state density), while moderate temperatures are reached (T < 70 MeV). The conditions closely resemble those that prevail, e.g., in neutron star mergers. Matter under such conditions is currently being studied at the High Acceptance DiElecton Spectrometer (HADES). Important topics of the research program are the mechanisms of strangeness production, the emissivity of matter, and the role of baryonic resonances herein. In this contribution, we will focus on the important experimental results obtained by HADES in Au+Au collisions at 2.4 GeV center-of-mass energy. We will also present perspectives for future experiments with HADES and CBM at SIS100, where higher beam energies and intensities will allow for the studies of the first-order deconfinement phase transition and its critical endpoint.

Exact Generalized Kohn-Sham Theory for Hybrid Functionals

p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; line-height: 18.0px; font: 15.8px Helvetica; color: #000000; -webkit-text-stroke: #000000; background-color: #ffffff} span.s1 {font-kerning: none} span.s2 {font-kerning: none; color: #000000} <p>Hybrid functionals have proven to be of immense practical value in density functional theory calculations. While they are often thought to be a heuristic construct, it has been established that this is in fact not the case. Here, we present a rigorous and formally exact generalized Kohn-Sham (GKS) density functional theory of hybrid functionals, in which exact remainder exchange-correlation potentials combine with a fraction of Fock exchange to produce the correct ground state density. Specifically, we generalize the well-known adiabatic con- nection theorem to the case of exact hybrid functional theory and use it to provide a rigorous distinction between multiplicative exchange and correlation components. We examine the exact theory by inverting reference electron densities to obtain exact GKS potentials for hybrid functionals, showing that an equivalent description of the many-electron problem is obtained with any arbitrary global fraction of Fock exchange. We establish the dependence of these exact components on the fraction of Fock exchange and use the observed trends to shed new light on the results of approximate hybrid functional calculations.</p>