Smooth periodic gauge satisfying crystal symmetry and periodicity to study high-harmonic generation in solids

Intense lasers can easily drive nonadiabatic transitions of excited electron wave packets across the Brillouin zones, thus transition dipole moments (TDM) between energy bands of solids should be continuous, satisfying crystal symmetry, and periodic at zone boundaries. While current ab initio algorithms are powerful in calculating band structures of solids, they all introduced random phases into the eigenfunctions at each crystal momentum k. Here we show how to choose a "smooth-periodic"gauge where TDMs can be smooth versus k, preserving crystal symmetry, as well as maintaining periodic at boundaries. The symmetry properties of TDMs with respect to k ensure the absence of even-order harmonics from MgO with inversion symmetry, while the TDM in the "smooth-periodic"gauge for broken-symmetry ZnO is responsible for even harmonics that were underestimated in previous simulations. These results reveal the importance of correctly treating the complex TDMs that satisfy crystal symmetry and continuous across zone boundaries in nonlinear laser-solid interactions, which has been elusive in most theories so far.