Text this: Multi-Domain Modeling of Free-Running Harmonic Frequency Comb Formation in Terahertz Quantum Cascade Lasers