Nonequilibrium Green’s function calculation for four-level scheme terahertz quantum cascade lasers | Synapse