First principles simulation of amorphous silicon bulk, interfaces, and nanowires for photovoltaics

Amorphous silicon has become the material of choice for many technologies, with major applications in large area electronics: displays, image sensing and thin film photovoltaic cells. This technology development has occurred because amorphous silicon is a thin film semiconductor that can be deposited on large, low cost substrates using low temperature. In this thesis, classical molecular dynamics and first principles DFT calculations have been performed to generate structural models of amorphous and hydrogenated amorphous silicon and interfaces of amorphous and crystalline silicon, with the ultimate aim of understanding the photovoltaic properties of core-shell crystalline amorphous Si nanowire structures. We have shown, unexpectedly, from the simulations, that our understanding of hydrogenated bulk a-Si needs to be revisited, with our robust finding that when fully saturated with hydrogen, bulk a-Si exhibits a constant optical energy gap, irrespective of the hydrogen concentration in the sample. Unsaturated a-Si:H, with a lower than optimum hydrogen content, shows a smaller optical gap, that increases with hydrogen content until saturation is reached. The mobility gaps obtained from an analysis of the electronic states show similar behavior. We also obtained that the optical and mobility gaps show a volcano curve as the H content is varied from 7% (undersaturation) to 18% (mild oversaturation). In the case of mild over saturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structures shows the extra H atoms in this case form a bridge between neighboring silicon atoms which increases the corresponding Si-Si distance and promotes bond length disorder in the sample. That has the potential to enhance the Staebler-Wronski effect. Planar interface models of amorphous-crystalline silicon have been generated in Si (100), (110) and (111) surfaces. The interface models are characterized by structure, RDF, electronic density of states and optical absorption spectrum. We find that the least stable (100) surface will result in the formation of the thickest amorphous silicon layer, while the most stable (110) surface forms the smallest amorphous region. We calculated for the first time band offsets of a-Si:H/c-Si heterojunctions from first principles and examined the influence of different surface orientations and amorphous layer thickness on the offsets and implications for device performance. The band offsets depend on the amorphous layer thickness and increase with thickness. By controlling the amorphous layer thickness we can potentially optimise the solar cell parameters. Finally, we have successfully generated different amorphous layer thickness of the a-Si/c-Si and a-Si:H/c-Si 5 nm nanowires from heat and quench. We perform structural analysis of the a-Si-/c-Si nanowires. The RDF, Si-Si bond length distributions, and the coordination number distributions of amorphous regions of the nanowires reproduce similar behaviour compared to bulk amorphous silicon. In the final part of this thesis we examine different surface terminating chemical groups, -H, - OH and –NH2 in (001) GeNW. Our work shows that the diameter of Ge nanowires and the nature of surface terminating groups both play a significant role in both the magnitude and the nature of the nanowire band gaps, allowing tuning of the band gap by up to 1.1 eV. We also show for the first time how the nanowire diameter and surface termination shifts the absorption edge in the Ge nanowires to longer wavelengths. Thus, the combination of nanowire diameter and surface chemistry can be effectively utilised to tune the band gaps and thus light absorption properties of small diameter Ge nanowires.

First-principles calculation of carrier-phonon scattering in n-type Si1−xGex alloys

First-principles electronic structure methods are used to find the rates of inelastic intravalley and intervalley n-type carrier scattering in Si1-xGex alloys. Scattering parameters for all relevant Delta and L intra- and intervalley scattering are calculated. The short-wavelength acoustic and the optical phonon modes in the alloy are computed using the random mass approximation, with interatomic forces calculated in the virtual crystal approximation using density functional perturbation theory. Optical phonon and intervalley scattering matrix elements are calculated from these modes of the disordered alloy. It is found that alloy disorder has only a small effect on the overall inelastic intervalley scattering rate at room temperature. Intravalley acoustic scattering rates are calculated within the deformation potential approximation. The acoustic deformation potentials are found directly and the range of validity of the deformation potential approximation verified in long-wavelength frozen phonon calculations. Details of the calculation of elastic alloy scattering rates presented in an earlier paper are also given. Elastic alloy disorder scattering is found to dominate over inelastic scattering, except for almost pure silicon (x approximate to 0) or almost pure germanium (x approximate to 1), where acoustic phonon scattering is predominant. The n-type carrier mobility, calculated from the total (elastic plus inelastic) scattering rate, using the Boltzmann transport equation in the relaxation time approximation, is in excellent agreement with experiments on bulk, unstrained alloys..

First-principles calculation of P-type alloy scattering in Si1-xGex

The p-type carrier scattering rate due to alloy disorder in Si1-xGex alloys is obtained from first principles. The required alloy scattering matrix elements are calculated from the energy splitting of the valence bands, which arise when one average host atom is replaced by a Ge or Si atom in supercells containing up to 128 atoms. Alloy scattering within the valence bands is found to be characterized by a single scattering parameter. The hole mobility is calculated from the scattering rate using the Boltzmann transport equation in the relaxation time approximation. The results are in good agreement with experiments on bulk, unstrained alloys..

First-principles calculation of alloy scattering in GexSi1-x

First-principles electronic structure methods are used to find the rates of intravalley and intervalley n-type carrier scattering due to alloy disorder in Si1-xGex alloys. The required alloy scattering matrix elements are calculated from the energy splitting of nearly degenerate Bloch states which arises when one average host atom is replaced by a Ge or Si atom in supercells containing up to 128 atoms. Scattering parameters for all relevant Delta and L intravalley and intervalley alloy scattering are calculated. Atomic relaxation is found to have a substantial effect on the scattering parameters. f-type intervalley scattering between Delta valleys is found to be comparable to other scattering channels. The n-type carrier mobility, calculated from the scattering rate using the Boltzmann transport equation in the relaxation time approximation, is in excellent agreement with experiments on bulk, unstrained alloys.

A tight-binding analysis of the electronic properties of III-nitride semiconductors

This thesis divides into two distinct parts, both of which are underpinned by the tight-binding model. The first part covers our implementation of the tight-binding model in conjunction with the Berry phase theory of electronic polarisation to probe the atomistic origins of spontaneous polarisation and piezoelectricity as well as attempting to accurately calculate the values and coefficients associated with these phenomena. We first develop an analytic model for the polarisation of a one-dimensional linear chain of atoms. We compare the zincblende and ideal wurtzite structures in terms of effective charges, spontaneous polarisation and piezoelectric coefficients, within a first nearest neighbour tight-binding model. We further compare these to real wurtzite structures and conclude that accurate quantitative results are beyond the scope of this model but qualitative trends can still be described. The second part of this thesis deals with implementing the tight-binding model to investigate the effect of local alloy fluctuations in bulk AlGaN alloys and InGaN quantum wells. We calculate the band gap evolution of Al1_xGaxN across the full composition range and compare it to experiment as well as fitting bowing parameters to the band gap as well as to the conduction band and valence band edges. We also investigate the wavefunction character of the valence band edge to determine the composition at which the optical polarisation switches in Al1_xGaxN alloys. Finally, we examine electron and hole localisation in InGaN quantum wells. We show how the built-in field localises the carriers along the c-axis and how local alloy fluctuations strongly localise the highest hole states in the c-plane, while the electrons remain delocalised in the c-plane. We show how this localisation affects the charge density overlap and also investigate the effect of well width fluctuations on the localisation of the electrons.

Giant piezoresistance in silicon-germanium alloys

We use first-principles electronic structure methods to show that the piezoresistive strain gauge factor of single-crystalline bulk n-type silicon-germanium alloys at carefully controlled composition can reach values of G = 500, three times larger than that of silicon, the most sensitive such material used in industry today. At cryogenic temperatures of 4 K we find gauge factors of G = 135 000, 13 times larger than that observed in Si whiskers. The improved piezoresistance is achieved by tuning the scattering of carriers between different (Delta and L) conduction band valleys by controlling the alloy composition and strain configuration.