Generating phase screens to replace phase fluctuation caused by atmospheric turbulence is essential for simulation of light propagation through the atmosphere. Error between power spectral density of actual turbulence and traditional Kolmogorov model illustrates the importance of generating non-Kolmogorov phase screen. Meanwhile, methods used to generate phase screen at present show different kinds of disadvantages respectively. In this paper, we adopt a new method named "intrinsic embedding fractional Brownian motion (IE-FBM)". First, relationship between phase screen and FBM is analyzed. Next, principle of IE-FBM is clarified. We expand the correlation matrix and generate a stationary Gaussian surface through two fast Fourier transforms, which is the principle of intrinsic embedding. After that, we adjust the Gaussian surface into an FBM surface. Finally, simulation results demonstrate that IE-FBM combines advantages of traditional methods. Phase structure function becomes closer to theoretical value no matter how we set parameters of phase screen. Besides, both low and high frequency components of phase screen are sufficient and creases don't exist. In addition, time consumption reduces apparently. In conclusion, our method is comprehensively optimal choice to generate phase screen. © 2020 Elsevier GmbH