Numerical Solution of Quantum Tunneling

The quantum hydrodynamics dynamics (QHD) is introduced to study the probability of electron through the tunneling junction. The complex form of the wave function was transformed into polar form, the time evolution equations of C and S were given in the Lagrange framework, and the explicit central differencing scheme was adopted for time integration of the QHD equations. In the Lagrangian framework, the grid points of the system will inevitable become scattered and unstructured during the evolution, the spatially smooth functions C and S are especially suitable for rational function approximation method to enable the implementation of a robust numerical scheme. As an example, the dynamics of a Gaussian wave packet is calculated and results are presented in the figures. In these figures, wave packet penetration of the barrier can be clearly observed, showing the quantum effect obviously.