Analytic first-order properties from explicitly correlated many-body perturbation theory and Gaussian geminal basis