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

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.

A Higher-Order Accurate Operator Splitting Spectral Method for The Wigner-Poisson System

An accurate description of 2-D quantum transport in a double-gate metal oxide semiconductor filed effect transistor (dgMOSFET) requires a high-resolution solver to a coupled system of the 4-D Wigner equation and 2-D Poisson equation. In this paper, we propose an operator splitting spectral method to evolve such Wigner-Poisson (WP) system in 4-D phase space with high accuracy. After an operator splitting of the Wigner equation, the resulting two sub-equations can be solved analytically with spectral approximation in phase space. Meanwhile, we adopt a Chebyshev spectral method to solve the Poisson equation. Spectral convergence in phase space and a fourth-order accuracy in time are both numerically verified. Finally, we apply the proposed solver into simulating dgMOSFET, develop the steady states from long-time simulations and obtain numerically converged current-voltage (I-V) curves.

Solving the Wigner Equation with Signed Particles Monte Carlo for Chemically Relevant Potentials

This paper presents the Signed Particles Monte Carlo algorithm for the solution of the transient Wigner equation for potentials relevant in chemical physics. Benchmarks include the harmonic and the double well potentials.

Solving the Wigner Equation with Signed Particles Monte Carlo for Chemically Relevant Potentials

This paper presents the Signed Particles Monte Carlo algorithm for the solution of the transient Wigner equation for potentials relevant in chemical physics. Benchmarks include the harmonic and the double well potentials.