Gaussian process models for preliminary low-thrust trajectory optimization