High-order harmonic generation from solids using Houston States

Abstract

The Houston states (Krieger and Iafrate in Phys Rev B 33: 5494 1986), which are instantaneous eigenstates of the time-dependent Hamiltonian, are known as an appropriate basis function to observe the role of each band in high-order harmonic generation (HHG) in solids. However, due to the sharp blow-up of the transition matrix elements between the Houston states (Wu et al. in Phys Rev A 91: 043839 2015), the Houston state treatment is known to be highly numerically unstable if more than two conduction bands are included. We develop a smooth variable discretization (SVD) method with Houston states by exploiting the fact that the wave function varies smoothly in time to solve the time-dependent Schrodinger equation (TDSE) for a solid exposed to a strong laser field. We demonstrate that the SVD method with the Houston states is numerically stable and can avoid the blow-up. Using the SVD method with the Houston states, we investigate the role of each band for the multiple plateau structure in HHG, and we find that strongly coupled bands should be treated as a whole and that the traditional classification of the currents into interband and intraband currents is not a suitable way to understand the plateaus in the HHG spectrum.