Source author record

Klaas J. H. Giesbertz

Klaas J. H. Giesbertz 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

9works
11topics
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

9 published item(s)

preprint2022arXiv

Secondary kinetic peak in the Kohn-Sham potential and its connection to the response step

We consider a prototypical 1D model Hamiltonian for a stretched heteronuclear molecule and construct individual components of the corresponding KS potential, namely: the kinetic, the N - 1, and the conditional potentials. These components show very special features, such as peaks and steps, in regions where the density is drastically low. Some of these features are quite well known, whereas others, such as a secondary peak in the kinetic potential or a second bump in the conditional potential, are less or not known at all. We discuss these features building on the analytical model treated in J. Chem. Theory Comput. 14, 4151 (2018). In particular, we provide an explanation for the underlying mechanism which determines the appearance of both peaks in the kinetic potential and elucidate why these peaks delineate the region over which the plateau structure, due to the N - 1 potential, stretches. We assess the validity of the Heitler-London Ansatz at large but finite internuclear distance, showing that, if optimal orbitals are used, this model is an excellent approximation to the exact wavefunction. Notably, we find that the second natural orbital presents an extra node very far out on the side of the more electronegative atom.

preprint2020arXiv

Large coupling-strength expansion of the Møller-Plesset adiabatic connection: From paradigmatic cases to variational expressions for the leading terms

We study in detail the first three leading terms of the large coupling-strength limit of the adiabatic connection that has as weak-interaction expansion the Møller-Plesset perturbation theory. We first focus on the H atom, both in the spin-polarized and the spin-unpolarized case, reporting numerical and analytical results. In particular, we derive an asymptotic equation that turns out to have simple analytical solutions for certain channels. The asymptotic H atom solution for the spin-unpolarized case is then shown to be variationally optimal for the many-electron spin-restricted closed-shell case, providing expressions for the large coupling-strength density functionals up to the third leading order. We also analyze the H2 molecule and the uniform electron gas.

preprint2020arXiv

Light-matter hybrid-orbital-based first-principles methods: the influence of the polariton statistics

A detailed understanding of strong matter-photon interactions requires first-principle methods that can solve the fundamental Pauli-Fierz Hamiltonian of non-relativistic quantum electrodynamics efficiently. A possible way to extend well-established electronic-structure methods to this situation is to embed the Pauli-Fierz Hamiltonian in a higher-dimensional light-matter hybrid auxiliary configuration space. In this work we show the importance of the resulting hybrid Fermi-Bose statistics of the polaritons, which are the new fundamental particles of the ``photon-dressed'' Pauli-Fierz Hamiltonian for systems in cavities. We show that violations of these statistics can lead to unphysical results. We present an efficient way to ensure the proper symmetry of the underlying wave functions by enforcing representability conditions on the dressed one-body reduced density matrix. We further present a general prescription how to extend a given first-principles approach to polaritons and as an example introduce polaritonic Hartree-Fock theory. While being a single-reference method in polariton space, polaritonic Hartree-Fock is a multi-reference method in the electronic space, i.e. it describes electronic correlations. We also discuss possible applications to polaritonic QEDFT. We apply this theory to a lattice model and find that the more delocalized the bound-state wave function of the particles is, the stronger it reacts to photons. The main reason is that within a small energy range many states with different electronic configurations are available as opposed to a strongly bound (and hence energetically separated) ground-state wave function. This indicates that under certain conditions coupling to the quantum vacuum of a cavity can indeed modify ground state properties.

preprint2016arXiv

Avoiding the 4-index transformation in one-body reduced density matrix functional calculations for separable functionals

One of the major computational bottlenecks in one-body reduced density matrix (1RDM) functional theory is the evaluation of approximate 1RDM functionals and their derivatives. The reason is that more advanced approximate functionals are almost exclusively defined in the natural orbital basis, so a 4-index transformation of the two-electron integrals appears to be unavoidable. I will show that this is not the case and that so-called separable functionals can be evaluated much more efficiently, i.e. only at cubic cost in the basis size. Since most approximate functionals are actually separable, this new algorithm is an important development to make 1RDM functional theory calculations feasible for large electronic systems.

preprint2016arXiv

Invertibility of the retarded response functions for initial mixed states: application to one-body reduced density matrix functional theory

In [J. Chem. Phys. 143, 054102 (2015)] I have derived conditions to characterize the kernel of the retarded response function, under the assumption that the initial state is a ground state. In this article I demonstrate its generalization to mixed states (ensembles). To make the proof work, the weights in the ensemble need to be decreasing for increasing energies of the pure states from which the mixed state is constructed. The resulting conditions are not easy to verify, but under the additional assumptions that the ensemble weights are directly related to the energies and that the full spectrum of the Hamiltonian participates in the ensemble, it is shown that potentials only belong to the kernel of the retarded response function if they commute with the initial Hamiltonian. These additional assumptions are valid for thermodynamic ensembles, which makes this result also physically relevant. The conditions on the potentials for the thermodynamic ensembles are much stronger than in the pure state (zero temperature) case, leading to a much less involved kernel when the conditions are applied to the retarded one-body reduced density matrix response function.

preprint2015arXiv

Invertibility of retarded response functions for Laplace transformable potentials: application to one-body reduced density matrix functional theory

A theorem for the invertibility of arbitrary response functions is presented under the following conditions: the time-dependence of the potentials should be Laplace transformable and the initial state should be a ground state, though it might be degenerate. This theorem provides a rigorous foundation for all density-functional-like theories in the time-dependent linear response regime. Especially for time-dependent one-body reduced density matrix (1RDM) functional theory this is an important step forward, since a solid foundation has currently been lacking. The theorem is equally valid for static response functions in the non-degenerate case, so can be used to characterize the uniqueness of the potential in the ground state version of the corresponding density-functional-like theory. Such a classification of the uniqueness of the non-local potential in ground state 1RDM functional theory has been lacking for decades. With the aid of presented invertibility theorem presented here, a complete classification of the non-uniqueness of the non-local potential in 1RDM functional theory can be given for the first time.

preprint2014arXiv

Response calculations based on an independent particle system with the exact one-particle density matrix: polarizabilities

Recently, we have demonstrated that the problems finding a suitable adiabatic approximation in time-dependent one-body reduced density matrix functional theory can be remedied by introducing an additional degree of freedom to describe the system: the phase of the natural orbitals [Phys. Rev. Lett. 105, 013002 (2010), J. Chem. Phys. 133, 174119 (2010)]. In this article we will show in detail how the frequency-dependent response equations give the proper static limit ($ω\to0$), including the perturbation in the chemical potential, which is required in static response theory to ensure the correct number of particles. Additionally we show results for the polarizability for H$_2$ and compare the performance of two different two-electron functionals: the phase-including Löwdin-Shull functional and the density matrix form of the Löwdin-Shull functional.

preprint2013arXiv

Long-range interactions and the sign of natural amplitudes in two-electron systems

In singlet two-electron systems the natural occupation numbers of the one-particle reduced density matrix are given as squares of the natural amplitudes which are defined as the expansion coefficients of the two-electron wave function in a natural orbital basis. In this work we relate the sign of the natural amplitudes to the nature of the two-body interaction. We show that long-range Coulomb-type interactions are responsible for the appearance of positive amplitudes and give both analytical and numerical examples that illustrate how the long-distance structure of the wave function affects these amplitudes. We further demonstrate that the amplitudes show an avoided crossing behavior as function of a parameter in the Hamiltonian and use this feature to show that these amplitudes never become zero, except for special interactions in which infinitely many of them can become zero simultaneously when changing the interaction strength. This mechanism of avoided crossings provides an alternative argument for the non-vanishing of the natural occupation numbers in Coulomb systems.

preprint2013arXiv

Towards non-local density functionals by explicit modeling of the exchange-correlation hole in inhomogeneous systems

We put forward new approach for the development of a non-local density functional by a direct modeling of the shape of exchange-correlation (xc) hole in inhomogeneous systems. The functional is aimed at giving an accurate xc-energy and an accurate corresponding xc-potential even in difficult near-degeneracy situations such as molecular bond breaking. In particular we demand that: (1) the xc hole properly contains -1 electron, (2) the xc-potential has the asymptotic -1/r behavior outside finite systems and (3) the xc-potential has the correct step structure related to the derivative discontinuities of the xc-energy functional. None of the currently existing functionals satisfies all these requirements. These demands are achieved by screening the exchange hole in such a way that the pair-correlation function is symmetric and satisfies the sum-rule. These two features immediately imply (1) and (2) while the explicit dependence of the exchange hole on the Kohn-Sham orbitals implies (3). Preliminary calculations show an improved physical description of the dissociating hydrogen molecule. Though the total energy is still far from perfect, the binding curve from our non-local density functional provides a significant improvement over the local density approximation.