Gaussian processes can model complex experimental data and interatomic interactions while enabling rigorous uncertainty quantification and propagation within a Bayesian framework. We apply state-of-the-art GP methods, including local [1], spectral, and non-stationary kernel design [2], to develop next-generation molecular simulation tools.