2019/04/30 by Francisco Ramírez, Francisco F. Ramírez, Daniel Dundas +4 · 13 citations
Engineering · Mathematics · Physics and Astronomy · #Classical mechanics #Context (archaeology) #Equations of motion #Formalism (music) #Liouville equation #Master equation #Mathematical physics #Mathematics #Molecular Junctions and Nanostructures #Physics #Pure mathematics #Quantum #Quantum and electron transport phenomena #Quantum mechanics #Surface and Thin Film Phenomena #Von Neumann architecture #cond-mat.mes-hall
paper · pdf · open access · doi:10.1021/acs.jpcc.8b12319
published in The Journal of Physical Chemistry C 123(20), 12542-12555 (American Chemical Society) · 40 pages, 15 figures
openalex publication_date 2019/04/30 · arxiv created 2019/05/10 · arxiv updated 2019/05/14 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/05
The so-called driven Liouville–von Neumann equation is a dynamical formulation to simulate a voltage bias across a molecular system and to model a time-dependent current in a grand-canonical framework. This approach introduces a damping term in the equation of motion that drives the charge to a reference, out of equilibrium density. Originally proposed by Horsfield and co-workers, further work on this scheme has led to different coexisting versions of this equation. On the other hand, the multiple-probe scheme devised by Todorov and collaborators, known as the hairy-probes method, is a formal treatment based on Green’s functions that allows the electrochemical potentials in two regions of an open quantum system to be fixed. In this article, the equations of motion of the hairy-probes formalism are rewritten to show that, under certain conditions, they can assume the same algebraic structure as the driven Liouville–von Neumann equation in the form proposed by Morzan et al. ( J. Chem. Phys. 2017, 146, 044110). In this way, a new formal ground is provided for the latter, identifying the origin of every term. The performances of the different methods are explored using tight-binding time-dependent simulations in three trial structures, designated as ballistic, disordered, and resonant models. In the context of first-principles Hamiltonians, the driven Liouville–von Neumann approach is of special interest, because it does not require the calculation of Green’s functions. Hence, the effects of replacing the reference density based on the Green’s function by one obtained from an applied field are investigated, to gain a deeper understanding of the limitations and the range of applicability of the driven Liouville–von Neumann equation.