The conventional Ewald expression for the electrostatic energy and forces is recast in a form that can be evaluated to high accuracy in order N log(N) steps using fast Fourier transforms. The fast Fourier Poisson method does not rely on interpolation approaches or Taylor/multipole expansions, and can be easily integrated with conventional molecular dynamics algorithms.