An algorithm for Monte Carlo simulation of bremsstrahlung emission by electrons based on the framework of Super MC is presented in this paper with efficient and accurate methods to sample the angular distribution and energy of bremsstrahlung photons. The photon energy is sampled according to scaled energy-loss differential cross sections tabulated by Seltzer and Berger. A novel hybrid model for photon angular distribution by low-and high-energy incident electrons is developed. The model uses Tsai’s full form of angular distribution function with atomic form factors for high-energy incident electrons. For electrons o500 ke V, a simple efficient and accurate analytical distribution function is developed, using adjustable parameters determined from the fitting of numerical values of the shape functions tabulated by Kissel et al. The efficiency of sampling photon energy is 80%. Our angular sampling algorithm for high-energy electron bremsstrahlung based on Tsai distribution function is very efficient(sampling efficiency*70%) in the useful photon energy range. An algorithm for Monte Carlo simulation of bremsstrahlung emission by electrons based on the framework of Super MC is presented in this paper with efficient and accurate methods to sample the angular distribution and energy of bremsstrahlung photons. The photon energy is sampled according to scaled energy-loss differential cross sections tabulated by Seltzer and Berger. A novel hybrid model for photon angular distribution by low-and high-energy incident electrons is developed. The model uses Tsai’s full form of angular distribution function with atomic form factors for high-energy incident electrons. For electrons o 500 ke V, a simple efficient and accurate analytical distribution function is developed, using adjustable parameters determined from the fitting of numerical values ​​of the shape functions tabulated by Kissel et al. The efficiency of sampling photon energy is 80%. Our angular sampling algorithm for high-energy electron bremsstrahlung based on Tsai distribution function is v ery efficient (sampling efficiency * 70%) in the useful photon energy range.