Pulse profiles of millisecond pulsars and their Fourier amplitudes
  Approximate analytical formulae are derived for the pulse profile produced by small hot spots on a rapidly rotating neutron star. Its Fourier amplitudes and phases are calculated. The proposed formalism takes into account gravitational bending of light, Doppler effect, anisotropy of emission, and time delays. Its accuracy is checked with exact numerical calculations. 
Juri Poutanen, Andrei M. Beloborodov