Search

The field of molecular electronics has emerged from efforts to understand electron propagation through single molecules and to use them in electronic circuits. Serving as a testbed for advanced theoretical methods, it reveals a significant discrepancy between the operational time scales of experiments (static to GHz frequencies) and theoretical models (femtoseconds). Utilizing a recently developed time-linear nonequilibrium Green's functions formalism, we model molecular junctions on experimentally accessible timescales. Our study focuses on the quantum pump effect in a Benzenedithiol molecule connected to two copper electrodes and coupled with cavity photons. By calculating both electric and photonic current responses to an ac bias voltage, we observe pronounced electroluminescence and high harmonic generation in this setup. The mechanism of the latter effect is more analogous to that from solids than from isolated molecules, with even harmonics being suppressed or enhanced depending on the symmetry of the driving field., Comment: 13 pages and 6 figures in the main text, 7 pages and 12 figures in the supporting information

The interaction of electrons with quantized phonons and photons underlies the ultrafast dynamics of systems ranging from molecules to solids, giving rise to a plethora of physical phenomena experimentally accessible using time-resolved techniques. Green's function methods offer an invaluable interpretation tool since scattering mechanisms of growing complexity can be selectively incorporated in the theory. cheers is a general-purpose nonequilibrium Green's function code that implements virtually all known many-body approximations and is designed for first principles studies of ultrafast processes in molecular and model solid state systems. The aims of generality, extensibility, efficiency, and user friendliness of the code are achieved through the underlying theory development and the use of modern software design practices. Here, we motivate the necessity for the creation of such a code and overview its design and capabilities., Comment: 15 pages, 4 figures

We present a time-linear scaling method to simulate open and correlated quantum systems out of equilibrium. The method inherits from many-body perturbation theory the possibility to choose selectively the most relevant scattering processes in the dynamics, thereby paving the way to the real-time characterization of correlated ultrafast phenomena in quantum transport. The open system dynamics is described in terms of an embedding correlator from which the time-dependent current can be calculated using the Meir-Wingreen formula. We show how to efficiently implement our approach through a simple grafting into recently proposed time-linear Green's function methods for closed systems. Electron-electron and electron-phonon interactions can be treated on equal footing while preserving all fundametal conservation laws., Comment: 6 pages and 3 figures in the main text, 4 pages and 2 figures in the supplemental material

We demonstrate the feasibility of the time-linear scaling formulation of the $GW$ method [Phys. Rev. Lett. {\bf 124}, 076601 (2020)] for {\it ab initio} simulations of optically driven two-dimensional materials. The time-dependent $GW$ equations are derived and solved numerically in the basis of Bloch states. We address carrier multiplication and relaxation in photo-excited graphene and find deviations from the typical exponential behavior predicted by the Markovian Boltzmann approach. For resonantly pumped semiconductor we discover a self-sustained screening cascade leading to the Mott transition of coherent excitons. Our results draw attention to the importance of non-Markovian and dynamical screening effects in out-of-equilibrium phenomena., Comment: 6 pages, 4 figures

Diagrammatic perturbation theory is a powerful tool for the investigation of interacting many-body systems, the self-energy operator $\Sigma$ 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.

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.

Ultrashort non-resonant electromagnetic pulses applied to effective one-electron systems may operate on the electronic state as a position or momentum translation operator. As derived here, extension to many-body correlated systems exposes qualitatively new aspects. For instance, to the lowest order in the electric field intensity the action of the pulse is expressible in terms of the two-body reduced density matrix enabling thus to probe various facets of electronic correlations. As an experimental realization we propose a pump-probe scheme in which after a weak, swift "kick" by the non-resonant pulse the survival probability for remaining in the initial state is measured. This probability we correlate to the two-body reduced density matrix. Since the strength of electronic correlation is bond-length sensitive, measuring the survival probability may allow for a direct insight into the bond-dependent two-body correlation in the ground state. As an illustration, full numerical calculations for two molecular systems are provided and different measures of electronic correlations are analyzed., Comment: 4 figures

This work provides a unified theoretical treatment of the single and correlated double-electron emission from a general electronic system. Using Feshbach projection method, the states of interest are selected by the projection operator; the Feshbach-Schur map determines the effective Hamiltonian and the optical potential for the emitted electrons. On the other hand, the nonequilibrium Green's functions method is demonstrated to be a complementary approach and an explicit correspondence between both methods is established. For a self-contained exposition some results on single electron emission are re-derived using both formalisms. New insights and results are obtained for the correlated electron-pair emission: This includes the effective two-electron Hamiltonian, the explicit form of the Feshbach self-energy in terms of the many-body self-energies, and the diagrammatic expansion of the two-particle current. As an illustration of the diagrammatic technique the process of the two-particle emission assisted by the excitation of plasmons is explicitly worked out., Comment: 8 figures

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 figures

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 figures

We investigate theoretically the relative time delay of photoelectrons originating from the different subshells (2s and 2p) of neon. This quantity was measured via attosecond streaking and studied theoretically by Schultze et al. [Science 328, 1658 (2010)]. A substantial discrepancy was found between the measured and the calculated values of the relative time delay. Several theoretical studies has been put forward to resolve this issue, e.g. by including correlation effects. In the present paper we explore the directional dependence of the photoelectron emission and the consequences for the inferred time delay. Our quantum mechanical calculations for an electron subject to laser fields and an effective single particle potential show that the time delay is indeed strongly angular dependent. Compared to strict forward emission we find that accounting for emission within a cone of 45 deg aperture, leads to a substantially increase of the relative time delay., Comment: This paper has been withdrawn by the author due some crucial errors regarding the summation over all m states

For tracing the spatiotemporal evolution of electronic systems, we suggest and analyze theoretically a setup that exploits the excellent spatial resolution based on scanning tunneling microscopy techniques combined with the temporal resolution of femtosecond pump-probe photoelectron spectroscopy. As an example, we consider the laser-induced, local vibrational dynamics of a surface-adsorbed molecule. The photoelectrons released by a laser pulse can be collected by the scanning tip and utilized to access the spatiotemporal dynamics. Our proof-of-principle calculations are based on the solution of the time-dependent Schrooedinger equation supported by the ab initio computation of the matrix elements determining the dynamics., Comment: Supporting Materials and Web-Enhanced Features under http://pubs.acs.org

We inspect the initial and the long time evolution of excitations a Fermi liquids by analyzing the time behavior of the electron spectral function. Focusing on the short-time limit we study the electron-boson model for the homogenous electron gas and apply the first order (in boson propagator) cumulant expansion of the electron Green's function. In addition to a quadratic decay in time upon triggering the excitation, we identify non-analytic terms in the time expansion similar to those found in the Fermi edge singularity phenomenon. We also demonstrate that the exponential decay in time in the long-time limit is inconsistent with the GW approximation for the self-energy. The background for this is the Paley-Wiener theorem of complex analysis. To reconcile with the Fermi liquid behavior an inclusion of higher order diagrams (in the screened Coulomb interaction) is required.

Second harmonic generation by spherical nanoparticles is a non-local optical process that can also be viewed as the result of the non-linear response of the thin interface layer. The classical electrodynamic description, based e.g. on the non-linear Mie theory, entails the knowledge of the dielectric function and the surface non-linear optical susceptibility, both quantities are usually assumed to be predetermined, for instance from experiment. We propose here an approach based on the semi-classical approximation for the quantum sum-over-states expression that allows to capture the second-order optical process from first principles. A key input is the electronic density, which can be obtained from effective single particle approaches such as the density-functional theory in the local density implementation. We show that the resulting integral equations can be solved very efficiently rendering thus the treatment of macroscopic systems. As an illustration we present numerical results for the magic Na2869 cluster., Comment: 11 pages, 5 figures

The collective response of matter is ubiquitous and widely exploited, e.g. in plasmonic, optical and electronic devices. Here we trace on an attosecond time scale the birth of collective excitations in a finite system and find distinct new features in this regime. Combining quantum chemical computation with quantum kinetic methods we calculate the time-dependent light absorption and refraction in fullerene that serve as indicators for the emergence of collective modes. We explain the numerically calculated novel transient features by an analytical model and point out the relevance for ultra-fast photonic and electronic applications. A scheme is proposed to measure the predicted effects via the emergent attosecond metrology., Comment: 11 pages, 3 figures, accepted in Phys. Rev. A

We present ab initio calculations of the quasiparticle decay times in a Buckminsterfullerene based on the many-body perturbation theory. A particularly lucid representation arises when the broadening of the quasiparticle states is plotted in the angular momentum and energy coordinates. In this representation the main spectroscopic features of the fullerene consist of two occupied nearly parabolic bands, and delocalized plane-wave-like unoccupied states with a few long-lived electronic states (the superatom molecular orbitals, SAMOs) embedded in the continuum of Fermi-liquid states. SAMOs have been recently uncovered experimentally by M. Feng, J. Zhao, and H. Petek [Science 320, 359 (2008)] using scanning tunneling spectroscopy. The present calculations offer an explanation of their unusual stability and unveil their long-lived nature making them good candidates for applications in the molecular electronics. From the fundamental point of view these states illustrate a concept of the Fock-space localization [B. L. Altshuler, Y. Gefen, A. Kamenev, and L. S. Levitov, Phys. Rev. Lett. 78, 2803 (1997)] with properties drastically different from the Fermi-liquid excitations.

The initial stages of the quasiparticle decay in a Fermi liquid are governed by a time-scale distinct from the scattering rates as derived from the Fermi golden rule approach. We show that the initial decay is nonexponential and that it is determined by the zeroth spectral moment of the electron self-energy. We analyzed numerically a number of approximations for the self-energy by comparing with exact configuration interaction calculations for small finite system with fragmented states. A numerically simple approach for computing the spreading of the quasiparticle states for large systems is devised., Comment: manuscript and supplement files