A Numerical Treatment of the Rf SQUID: I. General Properties and Noise Energy

Reinhold Kleiner,Dieter Koelle,John Clarke
DOI: https://doi.org/10.1007/s10909-007-9511-x
2007-01-30
Abstract:We investigate the characteristics and noise performance of rf Superconducting Quantum Interference Devices (SQUIDs) by solving the corresponding Langevin equations numerically and optimizing the model parameters with respect to noise energy. After introducing the basic concepts of the numerical simulations, we give a detailed discussion of the performance of the SQUID as a function of all relevant parameters. The best performance is obtained in the crossover region between the dispersive and dissipative regimes, characterized by an inductance parameter \beta_L'\equiv 2\pi L I_0/\Phi_0\approx 1; L is the loop inductance, I_0 the critical current of the Josephson junction, and \Phi_0 the flux quantum. In this regime, which is not well explored by previous analytical approaches, the lowest (intrinsic) values of noise energy are a factor of about 2 above previous estimates based on analytical approaches. However, several other analytical predictions, such as the inverse proportionality of the noise energy on the tank circuit quality factor and the square of the coupling coefficient between the tank circuit and the SQUID loop, could not be well reproduced. The optimized intrinsic noise energy of the rf SQUID is superior to that of the dc SQUID at all temperatures. Although for technologically achievable parameters this advantage shrinks, particularly at low thermal fluctuation levels, we give an example for realistic parameters that leads to a noise energy comparable to that of the dc SQUID even in this regime.
Superconductivity
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the numerical simulation and optimization of the characteristics and noise performance of radio - frequency superconducting quantum interference devices (rf SQUID). Specifically, the author carried out a detailed numerical simulation of the dynamic behavior of rf SQUID by solving the corresponding Langevin equation, and optimized the model parameters for noise energy. The focus of the research is on the performance of rf SQUID in different working modes, especially in the transition region between the dispersive regime and the dissipative regime, which has not been fully explored by previous analytical methods. ### Main research questions: 1. **Comparison between numerical simulation and analytical methods**: Verify and supplement the existing analytical theories through numerical simulation, especially in the transition region (i.e., the region where \( L\beta' \approx 1 \)), where analytical methods may be inaccurate or inapplicable. 2. **Optimization of noise performance**: Find the conditions that can achieve the lowest noise energy by optimizing the model parameters. This includes optimizing parameters such as the coupling coefficient \(\alpha\), the driving frequency \(f_d\), and the quality factor \(Q_0\). 3. **Comparison with dc SQUID**: Evaluate the noise performance of rf SQUID at different temperatures and parameters, and compare it with dc SQUID to explore whether rf SQUID has advantages under certain conditions. ### Key parameters and formulas: - **Normalized parameters**: - Normalized current \( i = I / I_0 \) - Normalized magnetic flux \( \phi = \Phi / \Phi_0 \) - Normalized voltage \( u = U / (I_0 R) \) - Normalized time \( \tau = t / (R I_0 / 2 \pi \Phi_0) \) - **Key parameters**: - \( \gamma_L = L_T / L \): Inductance ratio - \( L\beta' = 2 \pi L I_0 / \Phi_0 \): Normalized inductance parameter - \( \beta_c = 2 \pi C R I_0 / \Phi_0^2 \): Stewart - McCumber parameter - \( \alpha = M / \sqrt{L L_T} \): Coupling coefficient - \( Q_0 = \omega_0 L_T / R_T \): Quality factor - \( f_0 = \omega_0 / (2 \pi) \): Resonance frequency - \( \Gamma = k_B T / (I_0 \Phi_0) \): Noise parameter - **Main equations**: - **Josephson junction equation**: \[ I_J + R \dot{\delta} + \frac{\Phi_0}{2 \pi} \ddot{\delta} = I_N + I_0 \sin(\delta) \] - **Total magnetic flux relationship**: \[ \Phi = \Phi_{ext} + M I + L J \] - **Normalized equation**: \[ \ddot{\delta} + \beta_c \dot{\delta} + \sin(\delta) = \alpha \left[ \gamma_L \left( \frac{\phi_T}{\alpha} - \frac{\phi_{ext}}{L\beta'} \right) \right] \] ### Conclusion: Through numerical simulation, the author found that rf SQUID can achieve very low noise energy in the transition region where \( L\beta' \approx 1 \), and this result is better than previous analytical predictions. In addition, the optimized rf SQUID has low - frequency noise.