Properly consider deuteron in various initialisations and methods

Fixes a few bugs where deuteron was not being properly considered. The only significant one is in RichSellmeirFunc as this is the only one actively used in production.