We have computed a new grid of evolutionary subdwarf B star (sdB) models from the start of central He burning, taking into account atomic diffusion due to radiative levitation, gravitational settling, concentration diffusion, and thermal diffusion. We have computed the non-adiabatic pulsation properties of the models and present the predicted p-mode and g-mode instability strips. In previous studies of the sdB instability strips, artificial abundance enhancements of Fe and Ni were introduced in the pulsation driving layers. In our models, the abundance enhancements of Fe and Ni occur naturally, eradicating the need to use artificial enhancements. We find that the abundance increases of Fe and Ni were previously underestimated and show that the instability strip predicted by our simulations solves the so-called blue edge problem of the subdwarf B star g-mode instability strip. The hottest known g-mode pulsator, KIC 10139564, now resides well within the instability strip even when only modes with low spherical degrees (l ≤ 2) are considered.