We study numerically the Bloch electron wave-packet dynamics in periodic potentials to simulate laser-solid interactions. We introduce a quasi-classical model in the k space combined with the energy band structure to understand the high-order harmonic generation (HHG) process occurring in a subcycle timescale. This model interprets the multiple plateau structure in HHG spectra well and the linear dependence of cutoff energies on the amplitude of vector potential of the laser fields. It also predicts the emission time of HHG, which agrees well with the results by solving the timedependent Schrödinger equation (TDSE). It provides a scheme to reconstruct the energy dispersion relations in Brillouin zone and to control the trajectories of HHG by varying the shape of laser pulses. This model is instructive for experimental measurements.PACS numbers: 42.65. Ky, 42.65.Re, 72.20.Ht High-order harmonic generation (HHG) in atomic and molecular systems in the gas phase has been well studied theoretically and experimentally [1,2]. It can be understood by a semi-classical three-step model [3]. The bound electron may be ionized by tunneling, then driven by the external laser field. When it recombines with the parent ion, harmonics are emitted. The cutoff energy is around I p +3.17U p (U p is the ponderomotive energy, which is proportional to A 2 0 , where A 0 is the amplitude of the vector potential of the laser fields). The HHG has resulted in the birth of attosecond (1 as = 10 −18 s) pulses [4,5] and new imaging tools, such as molecular tomography [6] and spectroscopy [7]. Recent experiments have demonstrated that light-solid interactions offer a wide range of other phenomena and applications to be explored [8][9][10], including HHG from solid-state materials [11][12][13][14][15].Experimental results present a recollision feature in two-color laser fields [16], which is similar to HHG from the gas phase. However, it is also found that the cutoff of HHG from the solid phase depends on the strength E 0 of the laser fields linearly [11], rather than quadratically in the gas phase. Theoretical studies [17,18] show that the cutoff energy of HHG from the solid phase also depends on the laser wavelength λ linearly, rather than λ 2 in the gas phase. Multi-plateau structure in solid HHG has also been found theoretically [17,18] This work introduces a quasi-classical model to investigate the electron dynamic processes under the laser fields in the wave vector k space, which is similar to the three-step model for HHG generated from the atomic and molecular systems in the coordinate space [3].It is also described into three steps: Zener tunneling, electron wave packet oscillation in conduction bands, and an instantaneous electron-hole pair recombination. Since the electron and hole are delocalized, it is unnecessary for them to return to their original positions to emit HHG. The interband nonresonant and resonant Zener tunneling had been studied previously [21,[24][25][26]. The electronhole recollisions have also been observed in e...