Researcher profile

R. van Leeuwen

R. van Leeuwen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2020arXiv

Dynamically screened vertex correction to $GW$

Diagrammatic perturbation theory is a powerful tool for the investigation of interacting many-body systems, the self-energy operator $Σ$ encoding all the variety of scattering processes. In the simplest scenario of correlated electrons described by the $GW$ approximation for the electron self-energy, a particle transfers a part of its energy to neutral excitations. Higher-order (in screened Coulomb interaction $W$) self-energy diagrams lead to improved electron spectral functions (SF) by taking more complicated scattering channels into account and by adding corrections to lower order self-energy terms. However, they also may lead to unphysical negative spectral function. The resolution of this difficulty has been demonstrated in our previous works. The main idea is to represent the self-energy operator in a Fermi Golden rule form which leads to the manifestly positive definite SF and allows for a very efficient numerical algorithm. So far, the method has only been applied to 3D electron gas, which is a paradigmatic system, but a rather simple one. Here, we systematically extend the method to 2D including realistic systems such as mono and bilayer graphene. We focus on one of the most important vertex function effects involving the exchange of two particles in the final state. We demonstrate that it should be evaluated with the proper screening and discuss its influence on the quasiparticle properties.

preprint2014arXiv

Compact two-electron wave function for bond dissociation and Van der Waals interactions: A natural amplitude assessment

Electron correlations in molecules can be divided in short range dynamical correlations, long range Van der Waals type interactions and near degeneracy static correlations. In this work we analyze for a one-dimensional model of a two-electron system how these three types of correlations can be incorporated in a simple wave function of restricted functional form consisting of an orbital product multiplied by a single correlation function $f(r_{12})$ depending on the interelectronic distance $r_{12}$. Since the three types of correlations mentioned lead to different signatures in terms of the natural orbital (NO) amplitudes in two-electron systems we make an analysis of the wave function in terms of the NO amplitudes for a model system of a diatomic molecule. In our numerical implementation we fully optimize the orbitals and the correlation function on a spatial grid without restrictions on their functional form. Due to this particular form of the wave function, we can prove that none of the amplitudes vanishes and moreover that it displays a distinct sign pattern and a series of avoided crossings as a function of the bond distance in agreement with the exact solution. This shows that the wave function Ansatz correctly incorporates the long range Van der Waals interactions. We further show that the approximate wave function gives an excellent binding curve and is able to describe static correlations. We show that in order to do this the correlation function $f(r_{12})$ needs to diverge for large $r_{12}$ at large internuclear distances while for shorter bond distances it increases as a function of $r_{12}$ to a maximum value after which it decays exponentially. We further give a physical interpretation of this behavior.

preprint2014arXiv

Time-domain response of atomically thin $\mathrm{MoS_2}$ nanomechanical resonators

We measure the energy relaxation rate of single- and few-layer molybdenum disulphide ($\mathrm{MoS_2}$) nanomechanical resonators by detecting the resonator ring-down. Recent experiments on these devices show a remarkably low quality (Q)-factor when taking spectrum measurements at room temperature. The origin of the low spectral Q-factor is an open question, and it has been proposed that besides dissipative processes, frequency fluctuations contribute significantly to the resonance line-width. The spectral measurements performed thus far however, do not allow one to distinguish these two processes. Here, we use time-domain measurements to quantify the dissipation. We compare the Q-factor obtained from the ring-down measurements to those obtained from the thermal noise spectrum and from the frequency response of the driven device. In few-layer and single-layer $\mathrm{MoS_2}$ resonators the two are in close agreement, which demonstrates that the spectral line-width in $\mathrm{MoS_2}$ membranes at room temperature is limited by dissipation, and that excess spectral broadening plays a negligible role.

preprint2013arXiv

Exact Density-Functionals with Initial-State Dependence and Memory

We analytically construct the wave function that, for a given initial state, produces a prescribed density for a quantum ring with two non-interacting particles in a singlet state. In this case the initial state is completely determined by the initial density, the initial time-derivative of the density and a single integer that characterizes the (angular) momentum of the system. We then give an exact analytic expression for the exchange-correlation potential that relates two non-interacting systems with different initial states. This is used to demonstrate how the Kohn-Sham procedure predicts the density of a reference system without the need of solving the reference system's Schrödinger equation. We further numerically construct the exchange-correlation potential for an analytically solvable system of two electrons on a quantum ring with a squared cosine two-body interaction. For the same case we derive an explicit analytic expression for the exchange-correlation kernel and analyze its frequency-dependence (memory) in detail. We compare the result to simple adiabatic approximations and investigate the single-pole approximation. These approximations fail to describe the doubly-excited states, but perform well in describing the singly-excited states.

preprint2013arXiv

Natural occupation numbers: When do they vanish?

The non-vanishing of the natural orbital occupation numbers of the one-particle density matrix of many-body systems has important consequences for the existence of a density matrix-potential mapping for nonlocal potentials in reduced density matrix functional theory and for the validity of the extended Koopmans' Theorem. On the basis of Weyl's theorem we give a connection between the differentiability properties of the ground state wave function and the rate at which the natural occupations approach zero when ordered as a descending series. We show, in particular, that the presence of a Coulomb cusp in the wave function leads, in general, to a power law decay of the natural occupations, whereas infinitely differentiable wave-functions typically have natural occupations that decay exponentially. We analyze for a number of explicit examples of two-particle systems that in case the wave function is non-analytic at its spatial diagonal (for instance, due to the presence of a Coulomb cusp) the natural orbital occupations are non-vanishing. We further derive a more general criterium for the non-vanishing of NO occupations for two-particle wave functions with a certain separability structure. On the basis of this criterium we show that for a two-particle system of harmonically confined electrons with a Coulombic interaction (the so-called Hookium) the natural orbital occupations never vanish.

preprint2013arXiv

Strongly coupled modes in a weakly driven micromechanical resonator

We demonstrate strong coupling between the flexural vibration modes of a clamped-clamped micromechanical resonator vibrating at low amplitudes. This coupling enables the direct measurement of the frequency response via amplitude- and phase modulation schemes using the fundamental mode as a mechanical detector. In the linear regime, a frequency shift of $\mathrm{0.8\,Hz}$ is observed for a mode with a line width of $\mathrm{5.8\,Hz}$ in vacuum. The measured response is well-described by the analytical model based on the Euler-Bernoulli beam including tension. Calculations predict an upper limit for the room-temperature Q-factor of $\mathrm{4.5\times10^5}$ for our top-down fabricated micromechanical beam resonators.

preprint2012arXiv

Density-potential mappings in quantum dynamics

In a recent letter [Europhys. Lett. 95, 13001 (2011)] the question of whether the density of a time-dependent quantum system determines its external potential was reformulated as a fixed point problem. This idea was used to generalize the existence and uniqueness theorems underlying time-dependent density functional theory. In this work we extend this proof to allow for more general norms and provide a numerical implementation of the fixed-point iteration scheme. We focus on the one-dimensional case as it allows for a more in-depth analysis using singular Sturm-Liouville theory and at the same time provides an easy visualization of the numerical applications in space and time. We give an explicit relation between the boundary conditions on the density and the convergence properties of the fixed-point procedure via the spectral properties of the associated Sturm-Liouville operator. We show precisely under which conditions discrete and continuous spectra arise and give explicit examples. These conditions are then used to show that in the most physically relevant cases the fixed point procedure converges. This is further demonstrated with an example.

preprint2012arXiv

Local Control and v-Representability of Correlated Quantum Dynamics

We present a local control scheme to construct the external potential v that, for a given initial state, produces a prescribed time-dependent density in an interacting quantum many-body system. The numerical method is efficient and stable even for large and rapid density variations irrespective of the initial state and the interactions. The method can at the same time be used to answer fundamental v-representability questions in density-functional theory. In particular, in the absence of interactions, it allows us to construct the exact time-dependent Kohn-Sham potential for arbitrary initial states. We illustrate the method in a correlated one-dimensional two-electron system with different interactions, initial states and densities. For a Kohn-Sham system with a correlated initial state we demonstrate the interplay between memory and initial-state dependence as well as the failure of any adiabatic approximation.

preprint2011arXiv

Comparative study of many-body perturbation theory and time-dependent density functional theory in the out-of-equilibrium Anderson model

We study time-dependent electron transport through an Anderson model. The electronic interactions on the impurity site are included via the self-energy approximations at Hartree-Fock (HF), second Born (2B), GW, and T-Matrix level as well as within a time-dependent density functional (TDDFT) scheme based on the adiabatic Bethe-Ansatz local density approximation (ABALDA) for the exchange correlation potential. The Anderson model is driven out of equilibrium by applying a bias to the leads and its nonequilibrium dynamics is determined by real-time propagation. The time-dependent currents and densities are compared to benchmark results obtained with the time-dependent density matrix renormalization group (tDMRG) method. Many-body perturbation theory beyond HF gives results in close agreement with tDMRG especially within the 2B approximation. We find that the TDDFT approach with the ABALDA approximation produces accurate results for the densities on the impurity site but overestimates the currents. This problem is found to have its origin in an overestimation of the lead densities which indicates that the exchange correlation potential must attain nonzero values in the leads.

preprint2011arXiv

Kadanoff-Baym approach to double-excitations in finite systems

We benchmark many-body perturbation theory by studying neutral, as well as non-neutral, excitations of finite lattice systems. The neutral excitation spectra are obtained by time-propagating the Kadanoff-Baym equations in the Hartree-Fock and second Born approximations. Our method is equivalent to solving the Bethe-Salpeter equation with a high-level kernel while respecting self-consistently, which guarantees the fulfillment of a frequency sum rule. As a result, we find that a time-local method, such as Hartree-Fock, can give incomplete spectra, while already the second Born, which is the simplest time-nonlocal approximation, reproduces well most of the additional excitations, which are characterized as double-excitations.

preprint2011arXiv

On the Kirzhnits gradient expansion in two dimensions

We derive the semiclassical Kirzhnits expansion of the D-dimensional one-particle density matrix up to the second order in $\hbar$. We focus on the two-dimensional (2D) case and show that all the gradient corrections both to the 2D one-particle density and to the kinetic energy density vanish. However, the 2D Kirzhnits expansion satisfies the consistency criterion of Gross and Proetto [J. Chem. Theory Comput. 5, 844 (2009)] for the functional derivatives of the density and the noninteracting kinetic energy with respect to the Kohn-Sham potential. Finally we show that the gradient correction to the exchange energy diverges in agreement with the previous linear-response study.

preprint2011arXiv

Real-time switching between multiple steady-states in quantum transport

We study transport through an interacting model system consisting of a central correlated site coupled to finite bandwidth tight-binding leads, which are considered as effectively noninteracting. Its nonequilibrium properties are determined by real-time propagation of the Kadanoff-Baym equations after applying a bias voltage to the system. The electronic interactions on the central site are incorporated by means of self-energy approximations at Hartree-Fock, second Born and GW level. We investigate the conditions under which multiple steady-state solutions occur within different self-energy approximations, and analyze in detail the nature of these states from an analysis of their spectral functions. At the Hartree-Fock level at least two stable steady-state solutions with different densities and currents can be found. By applying a gate voltage-pulse at a given time we are able to switch between these solutions. With the same parameters we find only one steady-state solution when the self-consistent second Born and GW approximations are considered. We therefore conclude that treatment of many-body interactions beyond mean-field can destroy bistability and lead to qualitatively different results as compared those at mean-field level.

preprint2008arXiv

Nonequilibrium Green's functions approach to strongly correlated few-electron quantum dots

The effect of electron-electron scattering on the equilibrium properties of few-electron quantum dots is investigated by means of nonequilibrium Green's functions theory. The ground and equilibrium state is self-consistently computed from the Matsubara (imaginary time) Green's function for the spatially inhomogeneous quantum dot system whose constituent charge carriers are treated as spin-polarized. To include correlations, the Dyson equation is solved, starting from a Hartree-Fock reference state, within a conserving (second order) self-energy approximation where direct and exchange contributions to the electron-electron interaction are included on the same footing. We present results for the zero and finite temperature charge carrier density, the orbital-resolved distribution functions and the self-consistent total energies and spectral functions for isotropic, two-dimensional parabolic confinement as well as for the limit of large anisotropy--quasi-one-dimensional entrapment. For the considered quantum dots with N=2, 3 and 6 electrons, the analysis comprises the crossover from Fermi gas/liquid (at large carrier density) to Wigner molecule or crystal behavior (in the low-density limit).