Refactor creation of Gaussian charges
Use gaussian_charges()
Add a test for NaCl