Search

The current flux density is a vector field that can be used to describe theoretically how electrons flow in a system out-of-equilibrium. In this work, we unequivocally demonstrate that the signal obtained from time-resolved X-ray scattering does not only map the time-evolution of the electronic charge distribution, but also encodes information about the associated electronic current flux density. We show how the electronic current flux density qualitatively maps the distribution of electronic momenta and reveals the underlying mechanism of ultrafast charge migration processes, while also providing quantitative information about the timescales of electronic coherences., Comment: 6 pages, 4 figures

To assist the design of novel, highly efficient molecular junctions, a deep understanding of the precise charge transport mechanisms through these devices is of prime importance. In the present contribution, we describe a procedure to investigate spatially-resolved electron transport through a nanojunction from first principles, at the example of a nitro-substituted oligo-(phenylene ethynylene) covalently bound to graphene nanoribbon leads. Recently, we demonstrated that the conductivity of this single-molecule-graphene-nanoribbon junction can be switched quantitatively and reversibly upon application of a static electric field in a top gate position, in the spirit of a traditional field effect transistor [J. Phys. Chem. C, 2016, 120, 28808-28819]. The propensity of the central oligomer unit to align with the external field was found to induce a damped rotational motion and to cause an interruption of the conjugated $\pi$-system, thereby drastically reducing the conductance through the nanojunction. In the current work, we use the driven Liouville-von-Neumann (DLvN) approach for time-dependent electronic transport calculations to simulate the electronic current dynamics under time-dependent potential biases for the two logical states of the nanojunction. Our quantum dynamical simulations rely on a novel localization procedure using an orthonormal set of molecular orbitals obtained from a standard density functional theory calculation to generate a localized representation for the different parts of the molecular junction. The transparent DLvN formalism allows us to directly access the density matrix and to reconstruct the time-dependent electronic current density, unraveling unique mechanistic details of the electron transport.

In this contribution, we extend our framework for analyzing and visualizing correlated many-electron dynamics to non-variational, highly scalable electronic structure method. Specifically, an explicitly time-dependent electronic wave packet is written as a linear combination of $N$-electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time-dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open-source Python program detCI@ORBKIT, which extends the capabilities of our recently published post-processing toolbox [J. Comput. Chem. 37 (2016) 1511]. From the output of standard quantum chemistry packages using atom-centered Gaussian-type basis functions, the framework exploits the multi-determinental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one-electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser-driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher-level method provided a judicious choice of functional is made. Broadband excitation of a medium-sized organic chromophore further demonstrates the scalability of the method. In addition, the time-dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance.

Graphene pseudohalides are natural candidates for use in molecular sensing due to their greater chemical activity as compared to both graphene halides and pristine graphene. Though their study is still in its infancy, being hindered until recently by the unavailability of both selective and efficient procedures for their synthesis, they promise to considerably widen the application potential of chemically modified graphenes. Herein, we employ vdW-DFT to study the structural and electronic properties of two selected graphene pseudohalides namely cyanographone and isocyanographone and investigate the potential use of the latter as a chemical sensor via electron transport calculations.

The aim of the present contribution is to provide a framework for analyzing and visualizing the correlated many-electron dynamics of molecular systems, where an explicitly time-dependent electronic wave packet is represented as a linear combination of $N$-electron wave functions. The central quantity of interest is the electronic flux density, which contains all information about the transient electronic density, the associated phase, and their temporal evolution. It is computed from the associated one-electron operator by reducing the multi-determinantal, many-electron wave packet using the Slater-Condon rules. Here, we introduce a general tool for post-processing multi-determinant configuration-interaction wave functions obtained at various levels of theory. It is tailored to extract directly the data from the output of standard quantum chemistry packages using atom-centered Gaussian-type basis functions. The procedure is implemented in the open-source Python program detCI@ORBKIT, which shares and builds upon the modular design of our recently published post-processing toolbox [J. Comput. Chem. 37 (2016) 1511]. The new procedure is applied to ultrafast charge migration processes in different molecular systems, demonstrating its broad applicability. Convergence of the $N$-electron dynamics with respect to the electronic structure theory level and basis set size is investigated. This provides an assessment of the robustness of qualitative and quantitative statements that can be made concerning dynamical features observed in charge migration simulations.

ORBKIT is a toolbox for post-processing electronic structure calculations based on a highly modular and portable Python architecture. The program allows computing a multitude of electronic properties of molecular systems on arbitrary spatial grids from the basis set representation of its electronic wavefunction, as well as several grid-independent properties. The required data can be extracted directly from the standard output of a large number of quantum chemical programs. ORBKIT can be used as a standalone program to determine standard quantities, for example, the electron density, molecular orbitals, and derivatives thereof. The cornerstone of ORBKIT is its modular structure. The existing basic functions can be arranged in an individual way and can be easily extended by user-written modules to determine any other desired quantities. ORBKIT offers multiple output formats that can be processed by common visualization tools (VMD, Molden, etc.). Additionally, ORBKIT possesses routines to order molecular orbitals computed at different nuclear configurations according to their electronic character and to interpolate the wavefunction between these configurations. The program is open-source under GNU-LGPLv3 license and freely available at http://sourceforge.net/projects/orbkit/. This article provides an overview of ORBKIT with particular focus on its capabilities and applicability, and includes several example calculations.

The Born-Oppenheimer approximation leads to the counterintuitive result of a vanishing electronic flux density upon vibrational dynamics in the electronic ground state. To circumvent this long known issue, we propose using pairwise anti-symmetrically translated vibronic densities to generate a symmetric electronic density that can be forced to satisfy the continuity equation approximately. The so-called Born-Oppenheimer broken symmetry ansatz yields all components of the flux density simultaneously while requiring only knowledge about the nuclear quantum dynamics on the electronic adiabatic ground state potential energy surface. The underlying minimization procedure is transparent and computationally inexpensive, and the solution can be computed from the standard output of any quantum chemistry program. Taylor series expansion reveals that the implicit electron dynamics originates from non-adiabatic coupling to the explicit Born-Oppenheimer nuclear dynamics. The new approach is applied to the ${\rm H}_2^+$ molecular ion vibrating in its ${}^2\Sigma^+_g$ ground state. The electronic flux density is found to have the correct nodal structure and symmetry properties at all times.