Unphysical features in the application of the Boltzmann collision operator in the time-dependent modeling of quantum transport

Abstract

In this work, the use of the Boltzmann collision operator for dissipative quantum transport is analyzed. Its mathematical role on the description of the time-evolution of the density matrix during a collision can be understood as processes of adding and subtracting states. We show that unphysical results can be present in quantum simulations when the old states (that built the density matrix associated to an open system before the collision) are different from the additional states generated by the Boltzmann collision operator. As a consequence of the Fermi Golden rule, the new generated sates are usually eigenstates of the momentum or kinetic energy. Then, the different time-evolutions of old and new states involved in a collision process can originate negative values of the charge density, even longer after the collision. This unphysical feature disappears when the Boltzmann collision operator generates states that were already present in the density matrix of the quantum system before the collision. Following these ideas, in this paper, we introduce an algorithm that models phonon-electron interactions through the Boltzmann collision operator without negative values of the charge density. The model only requires the exact knowledge, at all times, of the states that build the density matrix of the open system.

Time-dependent quantum Monte Carlo simulation of electron devices with two-dimensional Dirac materials: A genuine terahertz signature for graphene

“…The limits of the reciprocal space {k x , k z } are selected indirectly by the occupation function f sum (E) in Eq. (12) in the main text. That is, the maximum value of the wave vector components, k x,max and k z,max , must be selected large enough to be sure that f sum (E(k x,max )) = f sum (E(k z,max )) ≈ 0.…”

“…In addition, because of the inherent quantum contextuality 14,15 , the predictions of these continuously measured system, in principle, depend on the type of measuring apparatus implemented in each model 16 . For the particular THz predictions of electron device invoked here, in addition, the relevant electrical current is the total current which is the sum of the conduction (flux of particles) plus the displacement (time-derivative of the electric field) components 12 . The displacement current, which is usually negligible for DC predictions, can no longer be ignored for THz predictions.…”

Identifying weak values with intrinsic dynamical properties in modal theories

“…This is not always true when computing probability with other more orthodox tools. For example, the Wigner function distribution, which provides simultaneous information on the position x = g and momentum p = s, has negative values [46]. It is in this sense that the Wigner distribution function is a quasiprobability distribution to account for its negatives values.…”

Simulation of Electronic Quantum Devices: Failure of Semiclassical Models

“…The two terms in the Wigner transport equation (Equation (1)) are then intrinsically conflicting: the deterministic one accounts for the QW subband structure via its confinement potential V(r), while the scattering one does not. This contradiction, recently pointed out also within the density-matrix formalism [71], is responsible for the dissipation slowdown as well as for the wrong thermalization dynamics. Indeed, the steady-state solution of the bulk-like Boltzmann collision term in Equation 10is just the Fermi-Dirac distribution.…”