We developed modeling capabilities, within the Vorpal particle-in-cell code, for three-dimensional simulations of surface effects and electron emission from semiconductor photocathodes. They include calculation of emission probabilities using general, piece-wise continuous, space-time dependent surface potentials, effective mass, and band bending field effects. We applied these models, in combination with previously implemented capabilities for modeling charge generation and transport in diamond, to investigate the emission dependence on applied electric field in the range from approximately 2 MV/m to 17 MV/m along the [100] direction. The simulation results were compared to experimental data. For the considered parameter regime, conservation of transverse electron momentum (in the plane of the emission surface) allows direct emission from only two (parallel to [100]) of the six equivalent lowest conduction band valleys. When the electron affinity χ is the only parameter varied in the simulations, the value χ = 0.31 eV leads to overall qualitative agreement with the probability of emission deduced from experiments. Including band bending in the simulations improves the agreement with the experimental data, particularly at low applied fields, but not significantly. Using surface potentials with different profiles further allows us to investigate the emission as a function of potential barrier height, width, and vacuum level position. However, adding surface patches with different levels of hydrogenation, modeled with position-dependent electron affinity, leads to the closest agreement with the experimental data.