Improved transfer matrix methods for calculating quantum   transmission-coefficient
  Methods for calculating the transmission coefficient are proposed, all of which arise from improved non-reflecting WKB boundary conditions at the edge of the computational domain in 1-dimensional geometries. In the first, the Schr\"{o}dinger equation is solved numerically while the second is a transfer matrix (TM) algorithm where the potential is approximated by steps, but with the first and last matrix modified to reflect the new boundary condition. Both methods give excellent results with first order WKB boundary conditions. The third uses the transfer matrix method with third order WKB boundary conditions. For the the parabolic potential, the average error for the modified third order TM method reduces by factor of 4100 over the unmodified TM method. 
Debabrata Biswas and Vishal Kumar
quant-ph cond-mat.mes-hall