Modeling Self-Starting Harmonic Mode Locking in Terahertz Quantum Cascade Lasers