Study of the Wigner function boundary conditions at different barrier heights

In this work, we compute the Wigner function from wavefunctions obtained by solving the Schrödinger Equation. Our goal is to investigate issues that we encounter when simulating devices by solving the Wigner transport equation, namely the numerical discrepancies between the boundaries and the neighboring regions. In this paper, we focus on how the boundary conditions are affected by the barrier height in single- and double-barrier devices.