Sodium laser guide star (LGS) becomes an essential part in modern astronomical adaptive optics system, especially for the next generation extremely large ground based telescope. The LGS technology has experienced the developmental stages as natural guide star, Rayleigh LGS, sodium LGS and constellation of LGS. The sky coverage is still limited in that the LGS cannot not be used to detect the tip/tilt aberrations. While the polychromatic laser guide star (PLGS) is one of the most effective ways to enlarge the sky coverage to 100%. Previous simulation models are insufficient for the accurate calculation of the return flux, especially for the simulation model of PLGS which is generated by one-photon excitation of mesospheric sodium atoms. The simulation model based on Bloch equations proposed in this paper can be used to compute the return flux of one-photon excited PLGS precisely. Doppler broadening, beam atom exchanging, collisions and recoil are taken into account in the model. The return flux is validated by the return efficiency. The simulation results indicate that with one-photon excitation of sodium atoms, a return efficiency of 330 nm is minimum compared with those of other wavelengths; the saturation power density will decrease with recoil increasing and increase with collision rate increasing; an optimal line-width exists up to maximum the photon return efficiency. In the best case, when the power density is 10 W/m2 at the sodium layer, the maximum return efficiency at 330 nm is 0.907 photons/s/sr/atom/(W/m2) with an optimal laser line-width of 18 MHz.