Non-linear two-frequency analysis of SIS mixers through harmonic balance

We present a numerical procedure for simulating the non-linear behaviour of Superconductor-Insulator-Superconductor (SIS) tunnel-junction mixers. Mixers of this kind are used extensively in high-performance submillimetre-wave astronomical receivers. The procedure is complete in that both the signal and local-oscillator voltages can be at high levels, with significant amounts of harmonic content; also, the full quantum mechanical behaviour of the tunnel junction is taken into account. The initial simulations show that multi-tone analysis can be carried out in a numerically efficient manner, and that saturation, and the attendant generation of IF harmonics, can be modelled rigorously. This work constitutes the first full non-linear model of an SIS mixer, and the theoretical and numerical procedures described, will be of great importance when designing the next generation of heterodyne astronomical receivers.