Under classical mechanics, the general equation of particle motion in the periodic field is derived. In the dampless case, the existence possibility of the higher-order harmonic radiation is explored by using Bessel function expansion of a generalized trigonometrical function and the multi-scale method. In the damping case, the critical properties and a chaotic behavior are discussed by the Melnikov method. The results show that the use of a higher-order harmonic radiation of non-relativistic particles as a short-wavelength laser source is perfectly possible, and the system's critical condition is related to its parameters. Only by adjusting parameters suitablely, the stable higher-order harmonic radiation with bigger intensity can be obtained.In the 70 s of the 20 th century, Kumakhov [1] discovered the channeling radiation which was produced by the interaction between the relativistic particles with strong lattice field (1000 T), and the radiation energy can enter the X-region or -region due to the Doppler effect. Since the 1980s, it has been predicted that the spontaneous channeling radiation can be converted into coherent radiation, and various possible solutions [2,3] have been proposed.It is known that when the relativistic particles interact with the material (crystal), the high energy radiation will be generated, but is needs a high-energy accelerator to provide the particle beams. However, due to the high cost of highenergy accelerator, the scheme is limited. Later, people discovered that the higher-order harmonic radiation of nonrelativistic particles was expected to reach this goal [4,5] . The researches showed [6][7][8] that the non-relativistic particles in the periodic field might generate a stronger higher-order harmonic radiation, and be able to enter X-region or soft X-region.One of the key problems is to find or design a periodic (potential) field for controlling the particle motion. M Z Shao et al [9] pointed out that the doped superlattice is a typical periodic material. motion. To ensure the stable higher-order harmonic radiation, in this paper we first derive the general equation of particle motion in the system consisting of doped superlattice and an external electric field along the growth direction of the superlattice. Based on the equation, the possibility of the higher-order harmonic radiation for dampless case is discussed by using the Bessel function expansion of the geralized trigonometric function and multi-scales method. Then, in damping case the stability, critical properties and chaotic behavior of the system are discussed by the Melnikov method.A typical periodic medium is the superlattice, and we are concerned about the doped superlattice. In the growth direction, the conduction band bottom of the substrate material is modulated periodically by a sine-like quantum well due to alternating doping [9] . The potential well depth can be controlled by doping concentration, while its width can be adjusted by the intrinsic thickness. If the beam of charged particles (e.g. electrons) al...