The wide-band limit is a commonly used approximation to analyze transport through nanoscale devices. In this work we investigate its applicability to the study of charge and heat transport through molecular break junctions exposed to voltage biases and temperature gradients. We find by comparative simulations that while the wide-band-limit approximation faithfully describes the long-time charge and heat transport, it fails to characterize the short-time behavior of the junction. In particular, we show that the charge current flowing through the device shows a discontinuity when a temperature gradient is applied, while the energy flow is discontinuous when a voltage bias is switched on and even diverges when the junction is exposed to both a temperature gradient and a voltage bias. We provide an explanation for this pathological behavior and propose two possible solutions to this problem.

Sudden ionisation of a relatively large molecule can initiate a correlation-driven process dubbed charge migration, where the electron density distribution is expected to rapidly move along the molecular backbone. Capturing this few-femtosecond or attosecond charge redistribution would represent the real-time observation of electron correlation in a molecule with the enticing prospect of following the energy flow from a single excited electron to the other coupled electrons in the system. Here, we report a time-resolved study of the correlation-driven charge migration process occurring in the nucleic-acid base adenine after ionisation with a 15–35 eV attosecond pulse. We find that the production of intact doubly charged adenine – via a shortly-delayed laser-induced second ionisation event – represents the signature of a charge inflation mechanism resulting from many-body excitation. This conclusion is supported by first-principles time-dependent simulations. These findings may contribute to the control of molecular reactivity at the electronic, few-femtosecond time scale.

The Auger decay is a relevant recombination channel during the first few femtoseconds of molecular targets impinged by attosecond XUV or soft X-ray pulses. Including this mechanism in timedependent simulations of charge-migration processes is a difficult task, and Auger scatterings are often ignored altogether. In this work we present an advance of the current state-of-the-art by putting forward a real-time approach based on nonequilibrium Green's functions suitable for firstprinciples calculations of molecules with tens of active electrons. To demonstrate the accuracy of the method we report comparisons against accurate grid simulations of one-dimensional systems. We also predict a highly asymmetric profile of the Auger wavepacket, with a long tail exhibiting ripples temporally spaced by the inverse of the Auger energy. arXiv:1806.03043v1 [cond-mat.mes-hall]

We have recently proposed a Nonequilibrium Green's Function (NEGF) approach to include Auger decay processes in the ultrafast charge dynamics of photoionized molecules. Within the so called Generalized Kadanoff-Baym Ansatz the fundamental unknowns of the NEGF equations are the reduced one-particle density matrix of bound electrons and the occupations of the continuum states. Both unknowns are one-time functions like the density in Time-Dependent Functional Theory (TDDFT). In this work we assess the accuracy of the approach against Configuration Interaction (CI) calculations in one-dimensional model systems. Our results show that NEGF correctly captures qualitative and quantitative features of the relaxation dynamics provided that the energy of the Auger electron is much larger than the Coulomb repulsion between two holes in the valence shells. For the accuracy of the results dynamical electron-electron correlations or, equivalently, memory effects play a pivotal role. The combination of our NEGF approach with the Sham-Schlüter equation may provide useful insights for the development of TDDFT exchange-correlation potentials with a history dependence.

In the nonequilibrium Green's function approach, the approximation of the correlation self-energy at the second-Born level is of particular interest, since it allows for a maximal speed-up in computational scaling when used together with the Generalized Kadanoff-Baym Ansatz for the Green's function. The present day numerical time-propagation algorithms for the Green's function are able to tackle first principles simulations of atoms and molecules, but they are limited to relatively small systems due to unfavourable scaling of self-energy diagrams with respect to the basis size. We propose an efficient computation of the self-energy diagrams by using tensorcontraction operations to transform the internal summations into functions of external low-level linear algebra libraries. We discuss the achieved computational speed-up in transient electron dynamics in selected molecular systems.

The dynamics of biologically relevant molecules exposed to ionizing radiation contains many facets and spans several orders of magnitude in time and energy. In the extreme ultraviolet (XUV) spectral range, multi-electronic phenomena and bands of correlated states with inner-valence holes must be accounted for in addition to a plethora of vibrational modes and available dissociation channels. To be able to track changes in charge density and bond length during ultrafast reactions is an important endeavour towards more general abilities to simulate and control photochemical processes, possibly inspired by those that have evolved biologically. By using attosecond XUV pulses extending up to 35 eV and few-femtosecond near-infrared pulses, we have previously time-resolved correlated electronic dynamics and charge migration occurring in the biologically relevant molecule adenine after XUV-induced sudden ionization. Here, by using additional experimental data, we comprehensively report on both electronic and vibrational dynamics of this nucleobase in an energy range little explored to date with high temporal resolution. The time-dependent yields of parent and fragment ions in mass spectra are analysed to extract exponential time constants and oscillation periods. Together with time-dependent density functional theory and ab-initio Green’s function methods, we identify different vibrational and electronic processes. Beyond providing further insights into the XUV-induced dynamics of an important nucleobase, our work demonstrates that yields of specific dissociation outcomes can be influenced by sufficiently well-timed ultrashort pulses, therefore providing a new route for the control of the multi-electronic and dissociative dynamics of a DNA building block.

The search for exchange-correlation functionals going beyond the adiabatic approximation has always been a challenging task for time-dependent density-functional theory. Starting from known results and using symmetry properties, we put forward a nonadiabatic exchange-correlation functional for lattice models describing a generic transport setup. We show that this functional reduces to known results for a single quantum dot connected to one or two reservoirs and furthermore yields the adiabatic local-density approximation in the static limit. Finally, we analyze the features of the exchange-correlation potential and the physics it describes in a linear chain connected to two reservoirs where the transport is induced by a bias voltage applied to the reservoirs. We find that the Coulomb blockade is correctly described for a half-filled chain, while additional effects arise as the doping of the chain changes.

No abstract