Viscoelastic relaxation time is an important concept to characterize the viscoelastic response of materials, which is directly related to the interactions among the microscopic atoms of materials. Few studies have focused on the methods of characterizing viscoelastic relaxation time. To investigate how to represent viscoelastic relaxation time effectively, the viscoelastic relaxation times of the monoatomic Lennard-Jones system on 22 conditions in a range of T^ * = 0.85–5, ρ^* = 0.85–1, ε = 0.97–1, and σ = 0.8–1.3 are discussed from a microscopic perspective by the equilibrium molecular dynamics methods. Static viscoelasticity (viscosity η^*, high-frequency shear modulus G_\infty^* ) is calculated by the Green-Kubo formula, and the Fourier transform is applied to the calculation of dynamic viscoelasticity (storage modulus G'^* and loss modulus G''^* ). On this basis, the viscoelastic characteristic relaxation time ( \tau _\rmMD^*), Maxwell relaxation time ( \tau _\rmMaxwell^*) and the lifetime of the state of local atomic connectivity ( \tau _\rmLC^*) are calculated. The viscoelastic characteristic relaxation time \tau _\rmMD^*, defined when the two responses crossover, is the key measure of the period of such a stimulus when the storage modulus (elasticity) equals the loss modulus (viscosity). Maxwell relaxation time \tau _\rmMaxwell^* = \eta ^*/G_\infty ^*, where η^* is the static viscosity under infinitely low stimulus frequency (i.e., zero shear rate), G_\infty^* is the instantaneous shear modulus under infinitely high stimulus frequency, and \tau _\rmLC^* is the time it takes for an atom to lose or gain one nearest neighbor. The result is observed that \tau _\rmLC^* is closer to \tau _\rmMD^* than \tau _\rmMaxwell^*. But the calculation of \tau _\rmLC^* needs to take into count the trajectories of all atoms in a certain time range, which takes a lot of time and computing resources. Finally, in order to characterize viscoelastic relaxation time more easily, Kramers’ rate theory is used to describe the dissociation and association of atoms, according to the radial distribution functions. And a method of predicting the viscoelasticity of the monoatomic Lennard-Jones system is proposed and established. The comparison of all the viscoelastic relaxation times obtained above shows that \tau _\rmMaxwell^* is quite different from \tau _\rmMD^* at low temperature in the monoatomic Lennard-Jones system. Compared with \tau _\rmMaxwell^*, \tau _\rmLC^* is close to \tau _\rmMD^*. But the calculation of \tau _\rmLC^* requires a lot of time and computing resources. Most importantly, the relaxation time calculated by our proposed method is closer to \tau _\rmMD^*. The method of predicting the viscoelastic relaxation time of the monoatomic Lennard-Jones system is accurate and reliable, which provides a new idea for studying the viscoelastic relaxation time of materials.