The auxiliary field diffusion Monte Carlo method with a path constraint is described and results shown for the neutron matter equation of state and spin susceptibility. Results with novel constraining wave functions which include spin-orbit induced backflow or BCS pairing correlations are discussed.