Density Matrix Modelling of Terahertz Frequency Quantum Cascade Lasers: Steady State Analysis and Maxwell-Bloch Dynamics