High field density-functional-theory based Monte Carlo: 4 H -SiC impact ionization and velocity saturationA full-band ensemble Monte Carlo simulation has been used to study the high-field carrier transport properties of 4H-SiC. The complicated band structure of 4H-SiC requires the consideration of band-to-band tunneling at high electric fields. We have used two models for the band-to-band tunneling; one is based on the overlap test and the other on the solution of the multiband Schrödinger equations. The latter simulations have only been performed for holes in the c-axis direction, since the computer capacity requirement are exceedingly high. Impact-ionization transition rates and phonon scattering rates have been calculated numerically directly from the full band structure. Coupling constants for the phonon interaction have been deduced by fitting of the simulated low-field mobility as a function of lattice temperature to experimental data. Secondary hot electrons generated as a consequence of hole-initiated impact ionization are considered in the study for both models of band-to-band tunneling. When the multiband Schrödinger equation model is used for holes in the c-axis direction, a significant change in the electron energy distribution is found, since the hole impact-ionization rate is very much increased with this model. The secondary electrons increase the average energy of the electron distribution leading to a significant increase in the electron-initiated impact-ionization coefficients. Our simulation results clearly show that both electrons and holes have to be considered in order to understand electron-initiated impact ionization in 4H-SiC.