Matching Items (9)
Filtering by
- All Subjects: Monte Carlo method
- Creators: Vasileska, Dragica
Description
The goal of this research work is to develop a particle-based device simulator for modeling strained silicon devices. Two separate modules had to be developed for that purpose: A generic bulk Monte Carlo simulation code which in the long-time limit solves the Boltzmann transport equation for electrons; and an extension to this code that solves for the bulk properties of strained silicon. One scattering table is needed for conventional silicon, whereas, because of the strain breaking the symmetry of the system, three scattering tables are needed for modeling strained silicon material. Simulation results for the average drift velocity and the average electron energy are in close agreement with published data. A Monte Carlo device simulation tool has also been employed to integrate the effects of self-heating into device simulation for Silicon on Insulator devices. The effects of different types of materials for buried oxide layers have been studied. Sapphire, Aluminum Nitride (AlN), Silicon dioxide (SiO2) and Diamond have been used as target materials of interest in the analysis and the effects of varying insulator layer thickness have also been investigated. It was observed that although AlN exhibits the best isothermal behavior, diamond is the best choice when thermal effects are accounted for.
ContributorsQazi, Suleman (Author) / Vasileska, Dragica (Thesis advisor) / Goodnick, Stephen (Committee member) / Tao, Meng (Committee member) / Arizona State University (Publisher)
Created2013
Description
In very small electronic devices the alternate capture and emission of carriers at an individual defect site located at the interface of Si:SiO2 of a MOSFET generates discrete switching in the device conductance referred to as a random telegraph signal (RTS) or random telegraph noise (RTN). In this research work, the integration of random defects positioned across the channel at the Si:SiO2 interface from source end to the drain end in the presence of different random dopant distributions are used to conduct Ensemble Monte-Carlo ( EMC ) based numerical simulation of key device performance metrics for 45 nm gate length MOSFET device. The two main performance parameters that affect RTS based reliability measurements are percentage change in threshold voltage and percentage change in drain current fluctuation in the saturation region. It has been observed as a result of the simulation that changes in both and values moderately decrease as the defect position is gradually moved from source end to the drain end of the channel. Precise analytical device physics based model needs to be developed to explain and assess the EMC simulation based higher VT fluctuations as experienced for trap positions at the source side. A new analytical model has been developed that simultaneously takes account of dopant number variations in the channel and depletion region underneath and carrier mobility fluctuations resulting from fluctuations in surface potential barriers. Comparisons of this new analytical model along with existing analytical models are shown to correlate with 3D EMC simulation based model for assessment of VT fluctuations percentage induced by a single interface trap. With scaling of devices beyond 32 nm node, halo doping at the source and drain are routinely incorporated to combat the threshold voltage roll-off that takes place with effective channel length reduction. As a final study on this regard, 3D EMC simulation method based computations of threshold voltage fluctuations have been performed for varying source and drain halo pocket length to illustrate the threshold voltage fluctuations related reliability problems that have been aggravated by trap positions near the source at the interface compared to conventional 45 nm MOSFET.
ContributorsAshraf, Nabil Shovon (Author) / Vasileska, Dragica (Thesis advisor) / Schroder, Dieter (Committee member) / Goodnick, Stephen (Committee member) / Goryll, Michael (Committee member) / Arizona State University (Publisher)
Created2011
Description
In semiconductor physics, many properties or phenomena of materials can be brought to light through certain changes in the materials. Having a tool to define new material properties so as to highlight certain phenomena greatly increases the ability to understand that phenomena. The generalized Monte Carlo tool allows the user to do that by keeping every parameter used to define a material, within the non-parabolic band approximation, a variable in the control of the user. A material is defined by defining its valleys, energies, valley effective masses and their directions. The types of scattering to be included can also be chosen. The non-parabolic band structure model is used. With the deployment of the generalized Monte Carlo tool onto www.nanoHUB.org the tool will be available to users around the world. This makes it a very useful educational tool that can be incorporated into curriculums. The tool is integrated with Rappture, to allow user-friendly access of the tool. The user can freely define a material in an easy systematic way without having to worry about the coding involved. The output results are automatically graphed and since the code incorporates an analytic band structure model, it is relatively fast. The versatility of the tool has been investigated and has produced results closely matching the experimental values for some common materials. The tool has been uploaded onto www.nanoHUB.org by integrating it with the Rappture interface. By using Rappture as the user interface, one can easily make changes to the current parameter sets to obtain even more accurate results.
ContributorsHathwar, Raghuraj (Author) / Vasileska, Dragica (Thesis advisor) / Goodnick, Stephen M (Committee member) / Saraniti, Marco (Committee member) / Arizona State University (Publisher)
Created2011
Description
ABSTRACT An Ensemble Monte Carlo (EMC) computer code has been developed to simulate, semi-classically, spin-dependent electron transport in quasi two-dimensional (2D) III-V semiconductors. The code accounts for both three-dimensional (3D) and quasi-2D transport, utilizing either 3D or 2D scattering mechanisms, as appropriate. Phonon, alloy, interface roughness, and impurity scattering mechanisms are included, accounting for the Pauli Exclusion Principle via a rejection algorithm. The 2D carrier states are calculated via a self-consistent 1D Schrödinger-3D-Poisson solution in which the charge distribution of the 2D carriers in the quantization direction is taken as the spatial distribution of the squared envelope functions within the Hartree approximation. The wavefunctions, subband energies, and 2D scattering rates are updated periodically by solving a series of 1D Schrödinger wave equations (SWE) over the real-space domain of the device at fixed time intervals. The electrostatic potential is updated by periodically solving the 3D Poisson equation. Spin-polarized transport is modeled via a spin density-matrix formalism that accounts for D'yakanov-Perel (DP) scattering. Also, the code allows for the easy inclusion of additional scattering mechanisms and structural modifications to devices. As an application of the simulator, the current voltage characteristics of an InGaAs/InAlAs HEMT are simulated, corresponding to nanoscale III-V HEMTs currently being fabricated by Intel Corporation. The comparative effects of various scattering parameters, material properties and structural attributes are investigated and compared with experiments where reasonable agreement is obtained. The spatial evolution of spin-polarized carriers in prototypical Spin Field Effect Transistor (SpinFET) devices is then simulated. Studies of the spin coherence times in quasi-2D structures is first investigated and compared to experimental results. It is found that the simulated spin coherence times for GaAs structures are in reasonable agreement with experiment. The SpinFET structure studied is a scaled-down version of the InGaAs/InAlAs HEMT discussed in this work, in which spin-polarized carriers are injected at the source, and the coherence length is studied as a function of gate voltage via the Rashba effect.
ContributorsTierney, Brian David (Author) / Goodnick, Stephen (Thesis advisor) / Ferry, David (Committee member) / Akis, Richard (Committee member) / Saraniti, Marco (Committee member) / Vasileska, Dragica (Committee member) / Arizona State University (Publisher)
Created2011
Description
This dissertation addresses challenges pertaining to multi-junction (MJ) solar cells from material development to device design and characterization. Firstly, among the various methods to improve the energy conversion efficiency of MJ solar cells using, a novel approach proposed recently is to use II-VI (MgZnCd)(SeTe) and III-V (AlGaIn)(AsSb) semiconductors lattice-matched on GaSb or InAs substrates for current-matched subcells with minimal defect densities. CdSe/CdTe superlattices are proposed as a potential candidate for a subcell in the MJ solar cell designs using this material system, and therefore the material properties of the superlattices are studied. The high structural qualities of the superlattices are obtained from high resolution X-ray diffraction measurements and cross-sectional transmission electron microscopy images. The effective bandgap energies of the superlattices obtained from the photoluminescence (PL) measurements vary with the layer thicknesses, and are smaller than the bandgap energies of either the constituent material. Furthermore, The PL peak position measured at the steady state exhibits a blue shift that increases with the excess carrier concentration. These results confirm a strong type-II band edge alignment between CdSe and CdTe. The valence band offset between unstrained CdSe and CdTe is determined as 0.63 eV±0.06 eV by fitting the measured PL peak positions using the Kronig-Penney model. The blue shift in PL peak position is found to be primarily caused by the band bending effect based on self-consistent solutions of the Schrödinger and Poisson equations. Secondly, the design of the contact grid layout is studied to maximize the power output and energy conversion efficiency for concentrator solar cells. Because the conventional minimum power loss method used for the contact design is not accurate in determining the series resistance loss, a method of using a distributed series resistance model to maximize the power output is proposed for the contact design. It is found that the junction recombination loss in addition to the series resistance loss and shadowing loss can significantly affect the contact layout. The optimal finger spacing and maximum efficiency calculated by the two methods are close, and the differences are dependent on the series resistance and saturation currents of solar cells. Lastly, the accurate measurements of external quantum efficiency (EQE) are important for the design and development of MJ solar cells. However, the electrical and optical couplings between the subcells have caused EQE measurement artifacts. In order to interpret the measurement artifacts, DC and small signal models are built for the bias condition and the scan of chopped monochromatic light in the EQE measurements. Characterization methods are developed for the device parameters used in the models. The EQE measurement artifacts are found to be caused by the shunt and luminescence coupling effects, and can be minimized using proper voltage and light biases. Novel measurement methods using a pulse voltage bias or a pulse light bias are invented to eliminate the EQE measurement artifacts. These measurement methods are nondestructive and easy to implement. The pulse voltage bias or pulse light bias is superimposed on the conventional DC voltage and light biases, in order to control the operating points of the subcells and counterbalance the effects of shunt and luminescence coupling. The methods are demonstrated for the first time to effectively eliminate the measurement artifacts.
ContributorsLi, Jingjing (Author) / Zhang, Yong-Hang (Thesis advisor) / Tao, Meng (Committee member) / Schroder, Dieter (Committee member) / Vasileska, Dragica (Committee member) / Arizona State University (Publisher)
Created2012
Description
Thermal effects in nano-scaled devices were reviewed and modeling methodologies to deal with this issue were discussed. The phonon energy balance equations model, being one of the important previous works regarding the modeling of heating effects in nano-scale devices, was derived. Then, detailed description was given on the Monte Carlo (MC) solution of the phonon Boltzmann Transport Equation. The phonon MC solver was developed next as part of this thesis. Simulation results of the thermal conductivity in bulk Si show good agreement with theoretical/experimental values from literature.
ContributorsYoo, Seung Kyung (Author) / Vasileska, Dragica (Thesis advisor) / Ferry, David (Committee member) / Goodnick, Stephen (Committee member) / Arizona State University (Publisher)
Created2015
Description
Self-heating degrades the performance of devices in advanced technology nodes. Understanding of self-heating effects is necessary to improve device performance. Heat generation in these devices occurs at nanometer scales but heat transfer is a microscopic phenomena. Hence a multi-scale modeling approach is required to study the self-heating effects. A state of the art Monte Carlo device simulator and the commercially available Giga 3D tool from Silvaco are used in our study to understand the self heating effects. The Monte Carlo device simulator solves the electrical transport and heat generation for nanometer length scales accurately while the Giga 3D tool solves for thermal transport over micrometer length scales. The approach used is to understand the self-heating effects in a test device structure, composed of a heater and a sensor, fabricated and characterized by IMEC. The heater is the Device Under Test(DUT) and the sensor is used as a probe. Therefore, the heater is biased in the saturation region and the sensor is biased in the sub-threshold regime. Both are planar MOSFETs of gate length equal to 22 nm. The simulated I-V characteristics of the sensor match with the experimental behavior at lower applied drain voltages but differ at higher applied biases.
The self-heating model assumes that the heat transport within the device follows Energy Balance model which may not be accurate. To properly study heat transport within the device, a state of the art Monte Carlo device simulator is necessary. In this regard, the Phonon Monte Carlo(PMC) simulator is developed. Phonons are treated as quasi particles that carry heat energy. Like electrons, phonons obey a corresponding Boltzmann Transport Equation(BTE) which can be used to study their transport. The direct solution of the BTE for phonons is possible, but it is difficult to incorporate all scattering mechanisms. In the Monte Carlo based solution method, it is easier to incorporate different relevant scattering mechanisms. Although the Monte Carlo method is computationally intensive, it provides good insight into the physical nature of the transport problem. Hence Monte Carlo based techniques are used in the present work for studying phonon transport. Monte Carlo simulations require calculating the scattering rates for different scattering processes. In the present work, scattering rates for three phonon interactions are calculated from different approaches presented in the literature. Optical phonons are also included in the transport problem. Finally, the temperature dependence of thermal conductivity for silicon is calculated in the range from 100K to 900K and is compared to available experimental data.
The self-heating model assumes that the heat transport within the device follows Energy Balance model which may not be accurate. To properly study heat transport within the device, a state of the art Monte Carlo device simulator is necessary. In this regard, the Phonon Monte Carlo(PMC) simulator is developed. Phonons are treated as quasi particles that carry heat energy. Like electrons, phonons obey a corresponding Boltzmann Transport Equation(BTE) which can be used to study their transport. The direct solution of the BTE for phonons is possible, but it is difficult to incorporate all scattering mechanisms. In the Monte Carlo based solution method, it is easier to incorporate different relevant scattering mechanisms. Although the Monte Carlo method is computationally intensive, it provides good insight into the physical nature of the transport problem. Hence Monte Carlo based techniques are used in the present work for studying phonon transport. Monte Carlo simulations require calculating the scattering rates for different scattering processes. In the present work, scattering rates for three phonon interactions are calculated from different approaches presented in the literature. Optical phonons are also included in the transport problem. Finally, the temperature dependence of thermal conductivity for silicon is calculated in the range from 100K to 900K and is compared to available experimental data.
ContributorsShaik, Abdul Rawoof (Author) / Vasileska, Dragica (Thesis advisor) / Ferry, David (Committee member) / Goodnick, Stephan (Committee member) / Arizona State University (Publisher)
Created2016
Description
In this work, transport in nanowire materials and nanowire field effect transistors is studied using a full band Monte Carlo simulator within the tight binding basis. Chapter 1 is dedicated to the importance of nanowires and nanoscale devices in present day electronics and the necessity to use a computationally efficient tool to simulate transport in these devices. Chapter 2 discusses the calculation of the full band structure of nanowires based on an atomistic tight binding approach, particularly noting the use of the exact same tight binding parameters for bulk band structures as well as the nanowire band structures. Chapter 3 contains the scattering rate formula for deformation potential, polar optical phonon, ionized impurity and impact ionization scattering in nanowires using Fermi’s golden rule and the tight binding basis to describe the wave functions. A method to calculate the dielectric screening in 1D systems within the tight binding basis is also described. Importantly, the scattering rates of nanowires tends to the bulk scattering rates at high energies, enabling the use of the same parameter set that were fitted to bulk experimental data to be used in the simulation of nanowire transport. A robust and efficient method to model interband tunneling is discussed in chapter 4 and its importance in nanowire transport is highlighted. In chapter 5, energy relaxation of excited electrons is studied for free standing nanowires and cladded nanowires. Finally, in chapter 6, a full band Monte Carlo particle based solver is created which treats confinement in a full quantum way and the current voltage characteristics as well as the subthreshold swing and percentage of ballistic transport is analyzed for an In0.7Ga0.3As junctionless nanowire field effect transistor.
ContributorsHathwar, Raghuraj (Author) / Goodnick, Stephen M (Committee member) / Saraniti, Marco (Committee member) / Vasileska, Dragica (Committee member) / Ferry, David K. (Committee member) / Arizona State University (Publisher)
Created2016
Description
4H-SiC has been widely used in many applications. All of these benefit from its extremely high critical electric field and good electron mobility. For example, 4H-SiC possesses a critical field ten times higher than that of Si, which allows high-voltage blocking layers composed of 4H-SiC to be approximately a tenth the thickness of a comparable Si device. This, in turn, reduces the device on-resistance and power losses while maintaining the same high blocking capability.
Unfortunately, commercial TCAD tools like Sentaurus and Silvaco Atlas are based on the effective mass approximation, while most 4H-SiC devices are not operated under low electric field, so the parabolic-like band approximation does not hold anymore. Hence, to get more accurate and reliable simulation results, full-band analysis is needed. The first step in the development of a full-band device simulator is the calculation of the band structure. In this work, the empirical pseudopotential method (EPM) is adopted. The next task in the sequence is the calculation of the scattering rates. Acoustic, non-polar optical phonon, polar optical phonon and Coulomb scattering are considered. Coulomb scattering is treated in real space using the particle-particle-particle-mesh (P3M) approach. The third task is coupling the bulk full-band solver with a 3D Poisson equation solver to generate a full-band device simulator.
For proof-of-concept of the methodology adopted here, a 3D resistor is simulated first. From the resistor simulations, the low-field electron mobility dependence upon Coulomb scattering in 4H-SiC devices is extracted. The simulated mobility results agree very well with available experimental data. Next, a 3D VDMOS is simulated. The nature of the physical processes occurring in both steady-state and transient conditions are revealed for the two generations of 3D VDMOS devices being considered in the study.
Due to its comprehensive nature, the developed tool serves as a basis for future investigation of 4H-SiC power devices.
Unfortunately, commercial TCAD tools like Sentaurus and Silvaco Atlas are based on the effective mass approximation, while most 4H-SiC devices are not operated under low electric field, so the parabolic-like band approximation does not hold anymore. Hence, to get more accurate and reliable simulation results, full-band analysis is needed. The first step in the development of a full-band device simulator is the calculation of the band structure. In this work, the empirical pseudopotential method (EPM) is adopted. The next task in the sequence is the calculation of the scattering rates. Acoustic, non-polar optical phonon, polar optical phonon and Coulomb scattering are considered. Coulomb scattering is treated in real space using the particle-particle-particle-mesh (P3M) approach. The third task is coupling the bulk full-band solver with a 3D Poisson equation solver to generate a full-band device simulator.
For proof-of-concept of the methodology adopted here, a 3D resistor is simulated first. From the resistor simulations, the low-field electron mobility dependence upon Coulomb scattering in 4H-SiC devices is extracted. The simulated mobility results agree very well with available experimental data. Next, a 3D VDMOS is simulated. The nature of the physical processes occurring in both steady-state and transient conditions are revealed for the two generations of 3D VDMOS devices being considered in the study.
Due to its comprehensive nature, the developed tool serves as a basis for future investigation of 4H-SiC power devices.
ContributorsCheng, Chi-Yin (Author) / Vasileska, Dragica (Thesis advisor) / Goodnick, Stephen M (Thesis advisor) / Ponce, Fernando (Committee member) / Zhao, Yuji (Committee member) / Arizona State University (Publisher)
Created2020