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.

Synthesis and characterization of cyclotriphosphazenes containing silicon as single solid-state precursors for the formation of silicon/phosphorus nanostructured materials

The synthesis and characterization of new organosilicon derivatives of N3P3Cl6, N3P3[NH(CH2)3Si(OEt)3]6 (1), N3P3[NH(CH2)3Si(OEt)3]3[NCH3(CH2)3CN]3 (2), and N3P3[NH(CH2)3Si(OEt)3]3[HOC6H4(CH2)CN]3 (3) are reported. Pyrolysis of 1, 2, and 3 in air and at several temperatures results in nanostructured materials whose composition and morphology depend on the temperature of pyrolysis and the substituents of the phosphazenes ring. The products stem from the reaction of SiO2 with P2O5, leading to either crystalline Si5(PO4)6O, SiP2O7 or an amorphous phase as the glass Si5(PO4)6O/3SiO2·2P2O5, depending on the temperature and nature of the trimer precursors. From 1 at 800 °C, core−shell microspheres of SiO2 coated with Si5(PO4)6O are obtained, while in other cases, mesoporous or dense structures are observed. Atomic force microscopy examination after deposition of the materials on monocrystalline silicon wafers evidences morphology strongly dependent on the precursors. Isolated islands of size ∼9 nm are observed from 1, whereas dense nanostructures with a mean height of 13 nm are formed from 3. Brunauer−Emmett−Teller measurements show mesoporous materials with low surface areas. The proposed growth mechanism involves the formation of cross-linking structures and of vacancies by carbonization of the organic matter, where the silicon compounds nucleate. Thus, for the first time, unique silicon nanostructured materials are obtained from cyclic phosphazenes containing silicon.

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.