Numerical Solution of the Constrained Wigner Equation

Kosik, Robert; Cervenka, J.; Kosina, H.

doi:10.23919/sispad49475.2020.9241624

Cited by 1 publication

(2 citation statements)

References 10 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In [8], a spectral method for the constrained sigma equation with an absorbing potential layer has been presented. This amounts to a discrete solution of the constrained Wigner equation in k-space.…”

Section: Simulation Resultsmentioning

confidence: 99%

See 1 more Smart Citation

Numerical constraints and non-spatial open boundary conditions for the Wigner equation

2021

Self Cite

View full text Add to dashboard Cite

We discuss boundary value problems for the characteristic stationary von Neumann equation (stationary sigma equation) and the stationary Wigner equation in a single spatial dimension. The two equations are related by a Fourier transform in the non-spatial coordinate. In general, a solution to the characteristic equation does not produce a corresponding Wigner solution as the Fourier transform will not exist. Solution of the stationary Wigner equation on a shifted k-grid gives unphysical results. Results showing a negative differential resistance in IV-curves of resonant tunneling diodes using Frensley’s method are a numerical artefact from using upwinding on a coarse grid. We introduce the integro-differential sigma equation which avoids distributional parts at $$k=0$$ k = 0 in the Wigner transform. The Wigner equation for $$k=0$$ k = 0 represents an algebraic constraint needed to avoid poles in the solution at $$k=0$$ k = 0 . We impose the inverse Fourier transform of the integrability constraint in the integro-differential sigma equation. After a cutoff, we find that this gives fully homogeneous boundary conditions in the non-spatial coordinate which is overdetermined. Employing an absorbing potential layer double homogeneous boundary conditions are naturally fulfilled. Simulation results for resonant tunneling diodes from solving the constrained sigma equation in the least squares sense with an absorbing potential reproduce results from the quantum transmitting boundary with high accuracy. We discuss the zero bias case where also good agreement is found. In conclusion, we argue that properly formulated open boundary conditions have to be imposed on non-spatial boundaries in the sigma equation both in the stationary and the transient case. When solving the Wigner equation, an absorbing potential layer has to be employed.

show abstract

Section: Simulation Resultsmentioning

confidence: 99%

“…use the method described in [8] which avoids the explicit building of normal equations. The numerical method is evaluated by simulating a GaAs/AlGaAs RTD and comparison with the quantum transmitting boundary method which we take as the reference model.…”

Section: Simulation Resultsmentioning

confidence: 99%