Correlated sampling within the shifted contour auxiliary field Monte Carlo method, implemented using plane waves and pseudopotentials, allows computation of electronic forces on nuclei, potential energy differences, geometric and vibrotational spectroscopic constants. This is exemplified on the N-2 molecule, where it is demonstrated that it is possible to accurately compute forces, dissociation energies, bond length parameters, and harmonic frequencies.