Terahertz gain on shallow donor transitions in silicon