Photoionization and transient Wannier-Stark ladder in silicon: First-principles simulations versus Keldysh theory

Nonlinear photoionization of dielectrics and semiconductors is widely treated in the framework of the Keldysh theory whose validity is limited to photon energies that are small compared to the band gap and relatively low laser intensities. The time-dependent density functional theory (TDDFT) simulations, which are free of these limitations, enable one to gain insight into nonequilibrium dynamics of the electronic structure. Here we apply TDDFT to investigate the photoionization of silicon crystal by ultrashort laser pulses in a wide range of laser wavelengths and intensities and compare the results with predictions of the Keldysh theory. Photoionization rates derived from the simulations considerably exceed the data obtained with the Keldysh theory within the validity range of the latter. Possible reasons for the discrepancy are discussed and we provide fundamental data on the photoionization rates beyond the limits of the Keldysh theory. By investigating the features of the Stark shift as a function of photon energy and laser field strength, a manifestation of the transient Wannier-Stark ladder states is revealed, which become blurred with increasing laser field strength. Finally, it is shown that the TDDFT simulations can potentially provide reliable data on the electron damping time that is of high importance for large-scale modeling.