the MTJ STO and avoid convergence issues, enabling a reliable MTJ STO model. Moreover, efficient simulations of both the stand-alone MTJ STO model and the MTJ STO model combined with CMOS circuits, are demonstrated.

Comprehensive and Compact MTJ STO Model in Verilog-A
====================================================

As detailed in Part I [@Tingsu_2014TEDI], the characteristics of an MTJ STO include the DC operating point, operating frequency $\omega_\text{g}$, output peak power $P(\omega)$, and linewidth $2 \Delta \omega$ (the full width at half-maximum). These characteristics vary greatly as the biasing condition changes. The biasing condition for the MTJ STO is typically composed of the amplitude ${H}_\text{ext}$ and the in-plane angle $\phi_\text{ext}$ of the external magnetic field, as well as the applied DC current $I_\text{DC}$.

The complete Verilog-A code of the comprehensive and compact MTJ STO model is available from [@model_link].

Computational efficiency of the comprehensive model
---------------------------------------------------

To fully implement the MTJ STO model in Verilog-A, solving Eq.(4a, 4b) in [@Tingsu_2014TEDI] for the angle and magnitude of the effective field is realized solely by using Verilog-A. The effective field angle (see Eq.(4a) in [@Tingsu_2014TEDI]) is solved numerically using the fixed-point iteration method, based on which the effective field magnitude can be simply obtained using Eq.(4b) in 