Conserving Approximations in Time-Dependent Density Functional Theory

In the present work we propose a theory for obtaining successively better approximations to the linear response functions of time-dependent density or current-density functional theory. The new technique is based on the variational approach to many-body perturbation theory (MBPT) as developed during the sixties and later expanded by us in the mid nineties. Due to this feature the resulting response functions obey a large number of conservation laws such as particle and momentum conservation and sum rules. The quality of the obtained results is governed by the physical processes built in through MBPT but also by the choice of variational expressions. We here present several conserving response functions of different sophistication to be used in the calculation of the optical response of solids and nano-scale systems.Comment: 11 pages, 4 figures, revised versio

Invariance of the Kohn (sloshing) mode in a conserving theory

It is proven that the center of mass (COM or Kohn) oscillation of a many-body system in a harmonic trap coincides with the motion of a single particle as long as conserving approximations are applied to treat the interactions. The two conditions formulated by Kadanoff and Baym \cite{kb-book} are shown to be sufficient to preserve the COM mode. The result equally applies to zero and finite temperature, as well as to nonequilibrium situations, and to the linear and nonlinear response regimes

Many-body Green's function theory for electron-phonon interactions: ground state properties of the Holstein dimer

We study ground-state properties of a two-site, two-electron Holstein model describing two molecules coupled indirectly via electron-phonon interaction by using both exact diagonalization and self-consistent diagrammatic many-body perturbation theory. The Hartree and self-consistent Born approximations used in the present work are studied at different levels of self-consistency. The governing equations are shown to exhibit multiple solutions when the electron-phonon interaction is sufficiently strong whereas at smaller interactions only a single solution is found. The additional solutions at larger electron-phonon couplings correspond to symmetry-broken states with inhomogeneous electron densities. A comparison to exact results indicates that this symmetry breaking is strongly correlated with the formation of a bipolaron state in which the two electrons prefer to reside on the same molecule. The results further show that the Hartree and partially self-consistent Born solutions obtained by enforcing symmetry do not compare well with exact energetics, while the fully self-consistent Born approximation improves the qualitative and quantitative agreement with exact results in the same symmetric case. This together with a presented natural occupation number analysis supports the conclusion that the fully self-consistent approximation describes partially the bipolaron crossover. These results contribute to better understanding how these approximations cope with the strong localizing effect of the electron-phonon interaction.Comment: 9 figures, corrected typo

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.Comment: 9 pages, 2 figure

On the Executability of Interactive Computation

The model of interactive Turing machines (ITMs) has been proposed to characterise which stream translations are interactively computable; the model of reactive Turing machines (RTMs) has been proposed to characterise which behaviours are reactively executable. In this article we provide a comparison of the two models. We show, on the one hand, that the behaviour exhibited by ITMs is reactively executable, and, on the other hand, that the stream translations naturally associated with RTMs are interactively computable. We conclude from these results that the theory of reactive executability subsumes the theory of interactive computability. Inspired by the existing model of ITMs with advice, which provides a model of evolving computation, we also consider RTMs with advice and we establish that a facility of advice considerably upgrades the behavioural expressiveness of RTMs: every countable transition system can be simulated by some RTM with advice up to a fine notion of behavioural equivalence.Comment: 15 pages, 0 figure

Diagrammatic expansion for positive spectral functions beyond GW: Application to vertex corrections in the electron gas

We present a diagrammatic approach to construct self-energy approximations within many-body perturbation theory with positive spectral properties. The method cures the problem of negative spectral functions which arises from a straightforward inclusion of vertex diagrams beyond the GW approximation. Our approach consists of a two-steps procedure: we first express the approximate many-body self-energy as a product of half-diagrams and then identify the minimal number of half-diagrams to add in order to form a perfect square. The resulting self-energy is an unconventional sum of self-energy diagrams in which the internal lines of half a diagram are time-ordered Green's functions whereas those of the other half are anti-time-ordered Green's functions, and the lines joining the two halves are either lesser or greater Green's functions. The theory is developed using noninteracting Green's functions and subsequently extended to self-consistent Green's functions. Issues related to the conserving properties of diagrammatic approximations with positive spectral functions are also addressed. As a major application of the formalism we derive the minimal set of additional diagrams to make positive the spectral function of the GW approximation with lowest-order vertex corrections and screened interactions. The method is then applied to vertex corrections in the three-dimensional homogeneous electron gas by using a combination of analytical frequency integrations and numerical Monte-Carlo momentum integrations to evaluate the diagrams.Comment: 19 pages, 19 figure

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.Comment: 16 pages, 13 figure

Diagrammatic expansion for positive density-response spectra: Application to the electron gas

In a recent paper [Phys. Rev. B 90, 115134 (2014)] we put forward a diagrammatic expansion for the self-energy which guarantees the positivity of the spectral function. In this work we extend the theory to the density response function. We write the generic diagram for the density-response spectrum as the sum of partitions. In a partition the original diagram is evaluated using time-ordered Green's functions (GF) on the left-half of the diagram, antitime-ordered GF on the right-half of the diagram and lesser or greater GF gluing the two halves. As there exist more than one way to cut a diagram in two halves, to every diagram corresponds more than one partition. We recognize that the most convenient diagrammatic objects for constructing a theory of positive spectra are the half-diagrams. Diagrammatic approximations obtained by summing the squares of half-diagrams do indeed correspond to a combination of partitions which, by construction, yield a positive spectrum. We develop the theory using bare GF and subsequently extend it to dressed GF. We further prove a connection between the positivity of the spectral function and the analytic properties of the polarizability. The general theory is illustrated with several examples and then applied to solve the long-standing problem of including vertex corrections without altering the positivity of the spectrum. In fact already the first-order vertex diagram, relevant to the study of gradient expansion, Friedel oscillations, etc., leads to spectra which are negative in certain frequency domain. We find that the simplest approximation to cure this deficiency is given by the sum of the zero-th order bubble diagram, the first-order vertex diagram and a partition of the second-order ladder diagram. We evaluate this approximation in the 3D homogeneous electron gas and show the positivity of the spectrum for all frequencies and densities.Comment: 19 pages, 19 figure

Vertex corrections for positive-definite spectral functions of simple metals

We present a systematic study of vertex corrections in the homogeneous electron gas at metallic densities. The vertex diagrams are built using a recently proposed positive-definite diagrammatic expansion for the spectral function. The vertex function not only provides corrections to the well known plasmon and particle-hole scatterings, but also gives rise to new physical processes such as generation of two plasmon excitations or the decay of the one-particle state into a two-particles-one-hole state. By an efficient Monte Carlo momentum integration we are able to show that the additional scattering channels are responsible for the bandwidth reduction observed in photoemission experiments on bulk sodium, appearance of the secondary plasmon satellite below the Fermi level, and a substantial redistribution of spectral weights. The feasibility of the approach for first-principles band-structure calculations is also discussed