Acknowledgements

I would like to thank Paul Sherwood for his support during this work and Alexey Sokol and Gargi Dutta (UCL) for their help with the QM/MM cluster example. Huub van Dam (PNNL, USA) provided the silicate structure used for the Hessian benchmark calculations.