Implemented Jasche et al.'s ionization.   Still some bugs...