1
|
Oz A, Nitzan A, Hod O, Peralta JE. Electron Dynamics in Open Quantum Systems: The Driven Liouville-von Neumann Methodology within Time-Dependent Density Functional Theory. J Chem Theory Comput 2023; 19:7496-7504. [PMID: 37852250 PMCID: PMC10653109 DOI: 10.1021/acs.jctc.3c00311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Indexed: 10/20/2023]
Abstract
A first-principles approach to describe electron dynamics in open quantum systems driven far from equilibrium via external time-dependent stimuli is introduced. Within this approach, the driven Liouville-von Neumann methodology is used to impose open boundary conditions on finite model systems whose dynamics is described using time-dependent density functional theory. As a proof of concept, the developed methodology is applied to simple spin-compensated model systems, including a hydrogen chain and a graphitic molecular junction. Good agreement between steady-state total currents obtained via direct propagation and those obtained from the self-consistent solution of the corresponding Sylvester equation indicates the validity of the implementation. The capability of the new computational approach to analyze, from first principles, non-equilibrium dynamics of open quantum systems in terms of temporally and spatially resolved current densities is demonstrated. Future extensions of the approach toward the description of dynamical magnetization and decoherence effects are briefly discussed.
Collapse
Affiliation(s)
- Annabelle Oz
- Department
of Physical Chemistry, School of Chemistry, the Raymond and Beverly
Sackler Faculty of Exact Sciences, and the Sackler Center for Computational
Molecular and Materials Science, Tel Aviv
University, Tel Aviv, 6997801, Israel
| | - Abraham Nitzan
- Department
of Physical Chemistry, School of Chemistry, the Raymond and Beverly
Sackler Faculty of Exact Sciences, and the Sackler Center for Computational
Molecular and Materials Science, Tel Aviv
University, Tel Aviv, 6997801, Israel
- Department
of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19103, United States
| | - Oded Hod
- Department
of Physical Chemistry, School of Chemistry, the Raymond and Beverly
Sackler Faculty of Exact Sciences, and the Sackler Center for Computational
Molecular and Materials Science, Tel Aviv
University, Tel Aviv, 6997801, Israel
| | - Juan E. Peralta
- Department
of Physics, Central Michigan University, Mount Pleasant, Michigan 48859, United States
| |
Collapse
|
2
|
Elenewski JE, Wójtowicz G, Rams MM, Zwolak M. Performance of reservoir discretizations in quantum transport simulations. J Chem Phys 2021; 155:124117. [PMID: 34598565 DOI: 10.1063/5.0065799] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Quantum transport simulations often use explicit, yet finite, electronic reservoirs. These should converge to the correct continuum limit, albeit with a trade-off between discretization and computational cost. Here, we study this interplay for extended reservoir simulations, where relaxation maintains a bias or temperature drop across the system. Our analysis begins in the non-interacting limit, where we parameterize different discretizations to compare them on an even footing. For many-body systems, we develop a method to estimate the relaxation that best approximates the continuum by controlling virtual transitions in Kramers turnover for the current. While some discretizations are more efficient for calculating currents, there is little benefit with regard to the overall state of the system. Any gains become marginal for many-body, tensor network simulations, where the relative performance of discretizations varies when sweeping other numerical controls. These results indicate that typical reservoir discretizations have little impact on numerical costs for certain computational tools. The choice of a relaxation parameter is nonetheless crucial, and the method we develop provides a reliable estimate of the optimal relaxation for finite reservoirs.
Collapse
Affiliation(s)
- Justin E Elenewski
- Biophysical and Biomedical Measurement Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Gabriela Wójtowicz
- Jagiellonian University, Institute of Theoretical Physics, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Marek M Rams
- Jagiellonian University, Institute of Theoretical Physics, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Michael Zwolak
- Biophysical and Biomedical Measurement Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
3
|
Zwolak M. Analytic expressions for the steady-state current with finite extended reservoirs. J Chem Phys 2020; 153:224107. [PMID: 33317280 PMCID: PMC8356363 DOI: 10.1063/5.0029223] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Open-system simulations of quantum transport provide a platform for the study of true steady states, Floquet states, and the role of temperature, time dynamics, and fluctuations, among other physical processes. They are rapidly gaining traction, especially techniques that revolve around "extended reservoirs," a collection of a finite number of degrees of freedom with relaxation that maintains a bias or temperature gradient, and have appeared under various guises (e.g., the extended or mesoscopic reservoir, auxiliary master equation, and driven Liouville-von Neumann approaches). Yet, there are still a number of open questions regarding the behavior and convergence of these techniques. Here, we derive general analytical solutions, and associated asymptotic analyses, for the steady-state current driven by finite reservoirs with proportional coupling to the system/junction. In doing so, we present a simplified and unified derivation of the non-interacting and many-body steady-state currents through arbitrary junctions, including outside of proportional coupling. We conjecture that the analytic solution for proportional coupling is the most general of its form for isomodal relaxation (i.e., relaxing proportional coupling will remove the ability to find compact, general analytical expressions for finite reservoirs). These results should be of broad utility in diagnosing the behavior and implementation of extended reservoir and related approaches, including the convergence to the Landauer limit (for non-interacting systems) and the Meir-Wingreen formula (for many-body systems).
Collapse
Affiliation(s)
- Michael Zwolak
- Biophysical and Biomedical Measurement Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| |
Collapse
|
4
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
5
|
Chiang TM, Hsu LY. Quantum transport with electronic relaxation in electrodes: Landauer-type formulas derived from the driven Liouville-von Neumann approach. J Chem Phys 2020; 153:044103. [PMID: 32752664 DOI: 10.1063/5.0007750] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive the exact steady-state solutions for the simplest model systems of resonant tunneling and tunneling with destructive quantum interference from the driven Liouville-von Neumann (DLvN) approach. Under the finite-state lead condition (the two electrodes have finite states), we analyze the asymptotic behavior of the steady-state current in the two limits of electronic relaxation. Under the infinite-state lead condition, the steady-state solutions of the two model systems can be cast as Landauer-type current formulas. According to the formulas, we show that the transmission functions near the resonant peak and the antiresonant dip can be significantly influenced by electronic relaxation in the electrodes. Moreover, under intermediate and strong electronic relaxation conditions, we analytically show that the steady-state current of the DLvN approach dramatically deviates from the Landauer current when destructive quantum interference occurs. In the regime of zero electronic relaxation, our results are reduced to the Landauer formula, indicating that the DLvN approach is equivalent to the Landauer approach when the leads have infinite states without any electronic relaxation.
Collapse
Affiliation(s)
- Tse-Min Chiang
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan
| | - Liang-Yan Hsu
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan
| |
Collapse
|
6
|
Wójtowicz G, Elenewski JE, Rams MM, Zwolak M. Open System Tensor Networks and Kramers' Crossover for Quantum Transport. PHYSICAL REVIEW. A 2020; 101:10.1103/PhysRevA.101.050301. [PMID: 33367191 PMCID: PMC7754794 DOI: 10.1103/physreva.101.050301] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Tensor networks are a powerful tool for many-body ground states with limited entanglement. These methods can nonetheless fail for certain time-dependent processes-such as quantum transport or quenches-where entanglement growth is linear in time. Matrix-product-state decompositions of the resulting out-of-equilibrium states require a bond dimension that grows exponentially, imposing a hard limit on simulation timescales. However, in the case of transport, if the reservoir modes of a closed system are arranged according to their scattering structure, the entanglement growth can be made logarithmic. Here, we apply this ansatz to open systems via extended reservoirs that have explicit relaxation. This enables transport calculations that can access steady states, time dynamics and noise, and periodic driving (e.g., Floquet states). We demonstrate the approach by calculating the transport characteristics of an open, interacting system. These results open a path to scalable and numerically systematic many-body transport calculations with tensor networks.
Collapse
Affiliation(s)
- Gabriela Wójtowicz
- Jagiellonian University, Institute of Theoretical Physics, Lojasiewicza 11, 30-348 Kraków, Poland
| | - Justin E. Elenewski
- Biophysics Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD, USA
- Institute for Research in Electronics and Applied Physics, University of Maryland, College Park, MD, USA
| | - Marek M. Rams
- Jagiellonian University, Institute of Theoretical Physics, Lojasiewicza 11, 30-348 Kraków, Poland
| | - Michael Zwolak
- Biophysics Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD, USA
| |
Collapse
|
7
|
Oz A, Hod O, Nitzan A. Numerical Approach to Nonequilibrium Quantum Thermodynamics: Nonperturbative Treatment of the Driven Resonant Level Model Based on the Driven Liouville von-Neumann Formalism. J Chem Theory Comput 2019; 16:1232-1248. [DOI: 10.1021/acs.jctc.9b00999] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | | | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19103, United States
| |
Collapse
|
8
|
Ridley M, Gull E, Cohen G. Lead geometry and transport statistics in molecular junctions. J Chem Phys 2019; 150:244107. [DOI: 10.1063/1.5096244] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael Ridley
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
- The Raymond and Beverley Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Emanuel Gull
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| | - Guy Cohen
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
- The Raymond and Beverley Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
9
|
Oz I, Hod O, Nitzan A. Evaluation of dynamical properties of open quantum systems using the driven Liouville-von Neumann approach: methodological considerations. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1584338] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Inbal Oz
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, IL, Israel
- The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv, IL, Israel
| | - Oded Hod
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, IL, Israel
- The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv, IL, Israel
| | - Abraham Nitzan
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, IL, Israel
- The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv, IL, Israel
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA, USA
| |
Collapse
|
10
|
Crawford TD, Kumar A, Bazanté AP, Di Remigio R. Reduced‐scaling coupled cluster response theory: Challenges and opportunities. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1406] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- T. Daniel Crawford
- Department of Chemistry Virginia Tech, Blacksburg Virginia
- The Molecular Sciences Software Institute Blacksburg Virginia
| | - Ashutosh Kumar
- Department of Chemistry Virginia Tech, Blacksburg Virginia
| | | | - Roberto Di Remigio
- Department of Chemistry Virginia Tech, Blacksburg Virginia
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry University of Tromsø ‐ The Arctic University of Norway Tromsø Norway
| |
Collapse
|
11
|
You X, Ramakrishna S, Seideman T. Unified theory of plasmon-induced resonance energy transfer and hot electron injection processes for enhanced photocurrent efficiency. J Chem Phys 2018; 149:174304. [PMID: 30408995 DOI: 10.1063/1.5050209] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Plasmons in metal nanoparticles (MNPs) promise to enhance solar energy conversion in semiconductors. Two essential mechanisms of enhancement in the near-field regime are hot electron injection (HEI) and plasmon-induced resonance energy transfer (PIRET). Individual studies of both mechanisms indicate that the PIRET efficiency is limited by the short lifetime of the plasmon, whereas the hot electrons result from the plasmon decay. The development of a unified theory of the coupled HEI and PIRET processes is fundamentally interesting and necessary for making reliable predictions but is complicated by the multiple interactions between various components that participate in the enhancement process. In this paper, we use the model-Hamiltonian approach to develop a combined theoretical framework including both PIRET and HEI. The coupled dynamics as well as the time evolution of hot electron energy distribution are studied. The theory further predicts an interference-induced asymmetry in the spectral dependence of PIRET, which can be used to distinguish it from HEI. As the relative contributions of PIRET and HEI strongly depend on the size of the MNPs, this presents itself as a simple route to control the strength of their contributions. The results presented here can further guide future applications of plasmonic solar energy harvesting.
Collapse
Affiliation(s)
- Xinyuan You
- Graduate Program in Applied Physics, Northwestern University, Evanston, Illinois 60208, USA
| | - S Ramakrishna
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - Tamar Seideman
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
12
|
Elenewski JE, Gruss D, Zwolak M. Communication: Master equations for electron transport: The limits of the Markovian limit. J Chem Phys 2018; 147:151101. [PMID: 29055298 DOI: 10.1063/1.5000747] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Master equations are increasingly popular for the simulation of time-dependent electronic transport in nanoscale devices. Several recent Markovian approaches use "extended reservoirs"-explicit degrees of freedom associated with the electrodes-distinguishing them from many previous classes of master equations. Starting from a Lindblad equation, we develop a common foundation for these approaches. Due to the incorporation of explicit electrode states, these methods do not require a large bias or even "true Markovianity" of the reservoirs. Nonetheless, their predictions are only physically relevant when the Markovian relaxation is weaker than the thermal broadening and when the extended reservoirs are "sufficiently large," in a sense that we quantify. These considerations hold despite complete positivity and respect for Pauli exclusion at any relaxation strength.
Collapse
Affiliation(s)
- Justin E Elenewski
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Daniel Gruss
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Michael Zwolak
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
13
|
Affiliation(s)
- Zhu Ruan
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem , Jerusalem, Israel
| | - Roi Baer
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem , Jerusalem, Israel
| |
Collapse
|
14
|
You X, Ramakrishna S, Seideman T. Origin of Plasmon Lineshape and Enhanced Hot Electron Generation in Metal Nanoparticles. J Phys Chem Lett 2018; 9:141-145. [PMID: 29256610 DOI: 10.1021/acs.jpclett.7b03126] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Plasmon-generated hot carriers are currently being studied intensively for their role in enhancing the efficiency of photovoltaic and photocatalytic processes. Theoretical studies of the hot electrons subsystem have generated insight, but we show that a unified quantum-mechanical treatment of the plasmon and hot electrons reveals new physical phenomena. Instead of a unidirectional energy transfer process in Landau damping, back energy transfer is predicted in small metal nanoparticles (MNPs) within a model-Hamiltonian approach. As a result, the single Lorentzian plasmonic line shape is modulated by a multipeak structure, whose individual line width provides a direct way to probe the electronic dephasing. More importantly, the hot electron generation can be enhanced greatly by matching the incident energy to the peaks of the modulated line shape.
Collapse
Affiliation(s)
- Xinyuan You
- Graduate Program in Applied Physics and ‡Department of Chemistry, Northwestern University , Evanston, Illinois 60208, United States
| | - S Ramakrishna
- Graduate Program in Applied Physics and ‡Department of Chemistry, Northwestern University , Evanston, Illinois 60208, United States
| | - Tamar Seideman
- Graduate Program in Applied Physics and ‡Department of Chemistry, Northwestern University , Evanston, Illinois 60208, United States
| |
Collapse
|
15
|
Gruss D, Smolyanitsky A, Zwolak M. Communication: Relaxation-limited electronic currents in extended reservoir simulations. J Chem Phys 2017; 147:141102. [PMID: 29031250 PMCID: PMC5718372 DOI: 10.1063/1.4997022] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
Open-system approaches are gaining traction in the simulation of charge transport in nanoscale and molecular electronic devices. In particular, "extended reservoir" simulations, where explicit reservoir degrees of freedom are present, allow for the computation of both real-time and steady-state properties but require relaxation of the extended reservoirs. The strength of this relaxation, γ, influences the conductance, giving rise to a "turnover" behavior analogous to Kramers turnover in chemical reaction rates. We derive explicit, general expressions for the weak and strong relaxation limits. For weak relaxation, the conductance increases linearly with γ and every electronic state of the total explicit system contributes to the electronic current according to its "reduced" weight in the two extended reservoir regions. Essentially, this represents two conductors in series-one at each interface with the implicit reservoirs that provide the relaxation. For strong relaxation, a "dual" expression-one with the same functional form-results, except now proportional to 1/γ and dependent on the system of interest's electronic states, reflecting that the strong relaxation is localizing electrons in the extended reservoirs. Higher order behavior (e.g., γ2 or 1/γ2) can occur when there is a gap in the frequency spectrum. Moreover, inhomogeneity in the frequency spacing can give rise to a pseudo-plateau regime. These findings yield a physically motivated approach to diagnosing numerical simulations and understanding the influence of relaxation, and we examine their occurrence in both simple models and a realistic, fluctuating graphene nanoribbon.
Collapse
Affiliation(s)
- Daniel Gruss
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
- Maryland Nanocenter, University of Maryland, College Park, MD 20742, USA
| | - Alex Smolyanitsky
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
| | - Michael Zwolak
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| |
Collapse
|
16
|
Goings JJ, Lestrange PJ, Li X. Real‐time time‐dependent electronic structure theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1341] [Citation(s) in RCA: 84] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
| | | | - Xiaosong Li
- Department of ChemistryUniversity of Washington Seattle WA USA
| |
Collapse
|
17
|
Morzan UN, Ramírez FF, González Lebrero MC, Scherlis DA. Electron transport in real time from first-principles. J Chem Phys 2017; 146:044110. [PMID: 28147541 DOI: 10.1063/1.4974095] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
While the vast majority of calculations reported on molecular conductance have been based on the static non-equilibrium Green's function formalism combined with density functional theory (DFT), in recent years a few time-dependent approaches to transport have started to emerge. Among these, the driven Liouville-von Neumann equation [C. G. Sánchez et al., J. Chem. Phys. 124, 214708 (2006)] is a simple and appealing route relying on a tunable rate parameter, which has been explored in the context of semi-empirical methods. In the present study, we adapt this formulation to a density functional theory framework and analyze its performance. In particular, it is implemented in an efficient all-electron DFT code with Gaussian basis functions, suitable for quantum-dynamics simulations of large molecular systems. At variance with the case of the tight-binding calculations reported in the literature, we find that now the initial perturbation to drive the system out of equilibrium plays a fundamental role in the stability of the electron dynamics. The equation of motion used in previous tight-binding implementations with massive electrodes has to be modified to produce a stable and unidirectional current during time propagation in time-dependent DFT simulations using much smaller leads. Moreover, we propose a procedure to get rid of the dependence of the current-voltage curves on the rate parameter. This method is employed to obtain the current-voltage characteristic of saturated and unsaturated hydrocarbons of different lengths, with very promising prospects.
Collapse
Affiliation(s)
- Uriel N Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| | - Francisco F Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| | - Mariano C González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| | - Damián A Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| |
Collapse
|
18
|
Zelovich T, Hansen T, Liu ZF, Neaton JB, Kronik L, Hod O. Parameter-free driven Liouville-von Neumann approach for time-dependent electronic transport simulations in open quantum systems. J Chem Phys 2017. [DOI: 10.1063/1.4976731] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Affiliation(s)
- Tamar Zelovich
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Thorsten Hansen
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | - Zhen-Fei Liu
- Department of Physics, University of California, Berkeley, Berkeley, California 94720-7300, USA
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, Berkeley, California 94720-7300, USA
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Kavli Energy NanoSciences Institute, Berkeley, California 94720, USA
| | - Leeor Kronik
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Oded Hod
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|