High Harmonic Generation in a Kronig-Penney Model Solid

Abstract

In 2010 high harmonic generation (HHG) in solids was first observed where high order harmonics of a strong laser field's frequency were observed. HHG in solids is now a rapidly developing field that allows for exciting applications like fully solid state attosecond XUV sources and new ultrafast resolution imaging techniques for quantum dynamics in solids.

HHG in solids has been explained by two mechanisms: an interband mechanism, due to polarization associated with separate energy bands, and an intraband mechanism that results from nonlinearities and population changes associated with each individual band. While interband HHG has been seen in wide bandwidth semiconductors, intraband HHG has been observed in narrow bandwidth dielectrics. There has not yet been an explanation of the alternation of mechanisms with material differences. The main goal of this thesis is to attempt to provide a better understanding of the most important mechanisms and where they prevail. Although numerical modelling of HHG requires consideration of multiple energy bands, a two-band model consisting only of a valence band and a single conduction band can explain the most important mechanisms. This model requires a given material's band gap between its valence and conduction bands as well as dipole matrix elements between the bands. In this thesis we follow the Kronig-Penney model to develop a 1D delta-function potential model of solids to obtain these properties required of the two-band model. We implement this in a Wannier quasi-classical (WQC) model of interband HHG in semiconductors that explains the dominant dynamics leading to such through quasi-classical real space electron-hole pair trajectories. Although HHG in solids can be explained to be the result of a resonant process in which an electron-hole pair is generated in the first step, there are also virtual transition processes that lack consideration. These processes do not conserve energy and correspond to transitions to conduction bands resulting from field induced distortions of the ground state. We use methodology introduced by Keldysh for optical field ionization of atoms and solids along with the 1D delta-function potential model to quantify how both resonant and virtual transitions lead to HHG in solids for wide and low bandwidth solids.