Simulation of semiconductor devices with a local numerical approach

Abstract A numerical solution of the Drift-Diffusion Model for simulation of semiconductor devices based on the local meshless numerical method is presented. Numerical difficulties inherited from convection-dominated processes and high gradients near junctions typically results in oscillations within the solution. The difficulties can be alleviated by artificial dissipation schemes or by other stabilization approaches that often require a complex computation to improve the solution convergence. We applied a simple numerical approach with a local coupling and without special treatments of nonlinearities. The proposed approach is straightforward to implement and is suitable for parallel execution. We demonstrate the efficiency of the proposed methodology on a simulation of PN junction. The results are compared against previously published data with a good agreement achieved. The applicability of the proposed methodology is confirmed with the simulation of extended tests with more complicated geometries and more intense dynamics. The computational efficiency is demonstrated through the measurement of execution time and speedup on shared memory computer architecture.

[1]  R. Courant,et al.  On the solution of nonlinear hyperbolic differential equations by finite differences , 1952 .

[2]  Nir Tessler,et al.  Two-dimensional simulation of polymer field-effect transistor , 2001 .

[3]  Paul Lin,et al.  Improving multigrid performance for unstructured mesh drift–diffusion simulations on 147,000 cores , 2012 .

[4]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[5]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[6]  Sajeev John,et al.  Finite difference discretization of semiconductor drift-diffusion equations for nanowire solar cells , 2012, Comput. Phys. Commun..

[7]  Mateo Valero,et al.  Moving from petaflops to petadata , 2013, CACM.

[8]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[9]  Hamou Sadat,et al.  A new meshless approach for three dimensional fluid flow and related heat transfer problems , 2012, 1809.02502.

[10]  S. M. Sze.  Physics of semiconductor devices /2nd edition/ , 1981 .

[11]  B. Nayroles,et al.  Generalizing the finite element method: Diffuse approximation and diffuse elements , 1992 .

[12]  Satya N. Atluri,et al.  MESHLESS LOCAL PETROV-GALERKIN (MLPG) METHOD FOR CONVECTION DIFFUSION PROBLEMS , 2000 .

[13]  Evan F. Bollig,et al.  Solution to PDEs using radial basis function finite-differences (RBF-FD) on multiple GPUs , 2012, J. Comput. Phys..

[14]  Guirong Liu,et al.  A point interpolation meshless method based on radial basis functions , 2002 .

[15]  B. Šarler,et al.  H-Adaptive Local Radial Basis Function Collocation Meshless Method , 2011 .

[16]  A. Arefmanesh,et al.  Buoyancy-driven fluid flow and heat transfer in a square cavity with a wavy baffle—Meshless numerical analysis , 2013 .

[17]  Roman Trobec,et al.  Computational complexity and parallelization of the meshless local Petrov-Galerkin method , 2009 .

[18]  S. Sivoththaman,et al.  Modeling and simulation of spherical solar cells , 2005, Conference Record of the Thirty-first IEEE Photovoltaic Specialists Conference, 2005..

[19]  G. Kosec,et al.  Local strong form meshless method on multiple GraphicsProcessing Units , 2013 .

[20]  Randolph E. Bank,et al.  A new discretization scheme for the semiconductor current continuity equations , 1989, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[21]  F. Brezzi,et al.  A relationship between stabilized finite element methods and the Galerkin method with bubble functions , 1992 .

[22]  S. Selberherr,et al.  A review of hydrodynamic and energy-transport models for semiconductor device simulation , 2003, Proc. IEEE.

[23]  H. Power,et al.  A Scalable Meshless Formulation Based on RBF Hermitian Interpolation for 3D Nonlinear Heat Conduction Problems , 2010 .

[24]  Y.-M. Guo.  An Overrange Collocation Method , 2011 .

[25]  Numerical simulation study of current–voltage characteristics of a Schottky diode with inverse doped surface layer , 2013 .

[26]  M. Shur.  Physics of Semiconductor Devices , 1969 .

[27]  R. W. Kelsall,et al.  The Monte Carlo method for semiconductor device simulation , 1995 .

[28]  S. C. Fan,et al.  Local multiquadric approximation for solving boundary value problems , 2003 .

[29]  Graham F. Carey,et al.  Advanced Numerical Methods and Software Approaches for Semiconductor Device Simulation , 2000, VLSI Design.

[30]  Božidar Šarler,et al.  From Global to Local Radial Basis Function Collocation Method for Transport Phenomena , 2007 .

[31]  R. Trobec,et al.  Comparison of local weak and strong form meshless methods for 2-D diffusion equation , 2012 .

[32]  Thong-See Lee.  A false transient approach to steady state solution of fluid flow through vascular constrictions , 1991 .

[33]  Roman Trobec,et al.  A Parallel Meshless Numerical Approach for the Solution of Transport Phenomena , 2012 .

[34]  S. Sannasiraj,et al.  SPH Simulation of shallow water wavepropagation , 2013 .

[35]  Arnaud G. Malan,et al.  An artificial compressibility CBS method for modelling heat transfer and fluid flow in heterogeneous porous materials , 2011 .

[36]  Yong Duan,et al.  A numerical scheme for nonlinear Schrödinger equation by MQ quasi-interpolation , 2013 .

[37]  Asen Asenov,et al.  Enriched residual free bubbles for semiconductor device simulation , 2012 .

[38]  S. M. Sze,et al.  Physics of semiconductor devices , 1969 .