Raise an Error for missing density values
I've just had a case where I didn't define density for a GHB package. This results in the GHB being ignored (since one of its values is missing).
So I think the solution is never giving density a default value of None. I'm guessing density is also required even when concentration isn't (for multiple systems).