2012
DOI: 10.1063/1.4767460
|View full text |Cite
|
Sign up to set email alerts
|

An efficient solution of Liouville-von Neumann equation that is applicable to zero and finite temperatures

Abstract: Application of quantum dissipation theory to electronic dynamics has been limited to model systems with few energy levels, and its numerical solutions are mostly restricted to high temperatures. A highly accurate and efficient numerical algorithm, which is based on the Chebyshev spectral method, is developed to integrate a single-particle Liouville-von Neumann equation, and the two long-standing limitations of quantum dissipation theory are resolved in the context of quantum transport. Its computational time s… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
3
1

Citation Types

1
31
0

Year Published

2013
2013
2021
2021

Publication Types

Select...
8
1

Relationship

1
8

Authors

Journals

citations
Cited by 32 publications
(32 citation statements)
references
References 21 publications
1
31
0
Order By: Relevance
“…30 and similar to the method of Ref. 51, truncating at a certain K ch limits the approach to a simulation time dependent on the truncation order K ch . The uniform Chebyshev decomposition scheme shows an exponential convergence, generating extremely accurate results given the criteria for convergence are fulfilled 36 .…”
Section: Functionsmentioning
confidence: 99%
“…30 and similar to the method of Ref. 51, truncating at a certain K ch limits the approach to a simulation time dependent on the truncation order K ch . The uniform Chebyshev decomposition scheme shows an exponential convergence, generating extremely accurate results given the criteria for convergence are fulfilled 36 .…”
Section: Functionsmentioning
confidence: 99%
“…It should be noted that the density operators for a closed system should satisfy the following requirements throughout the whole evolution [3]: (a) Hermitian, so that all the probabilities are real, (b) trace-preserving, because the sum of the probabilities over any complete set must be one, (c) positive, otherwise some probabilities might be negative. Different methods such as Liouville-von Neumann equation have been developed for the solution of the spin density operators [2][3][4][5][6][7]. In the special case of two-level systems, the time evolution of the density matrix elements can be described by the classical Bloch equation [8].…”
Section: Introductionmentioning
confidence: 99%
“…10,22 Beyond the WBL, the RSDM based HEOM have been applied to simulate time-dependent quantum transport in tightbinding model system. 23,24 However, in the derivation of RSDM based HEOM, a localized orthonormal basis set is assumed. In practical first principles calculation, atomic orbital basis, which are the eigenfunctions of the single-electron Hamiltonian of individual atoms, are commonly used due to their localized nature and clear chemical meaning.…”
Section: Introductionmentioning
confidence: 99%