It turns out the NSE solver works just fine at lower temperatures, and the assumed mass fractions screw it up. Just commented out the assumed values below T=0.5 to let the NSE solver act everywhere.