P. Berggren,J. Fransson
Abstract:We consider the motion of and Josephson current through a mechanically oscillating superconducting island asymmetrically embedded in a Josephson junction. The electromechanical coupling is provided by distance dependent tunneling rates between the electrodes and the island. The system asymmetry, resulting from the geometrical configuration, leads, for weak coupling, to an equation of the mechanical motion that reduces to the well-known Duffing equation. At zero bias voltage the island motion is determined by the homogenous Duffing equation that opens up two separate regions of solutions depending on the superconducting phases. The island either moves under influence of an anharmonic single well potential, or is governed by a double well potential that allows for off-center oscillations. Under applied bias voltage the island equation of motion turns into a modified Duffing equation, with time dependent coefficients, that demonstrate both quasi periodic and chaotic behavior.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: to study the coupling effect between mechanical oscillators and electron tunneling in a double - Josephson - junction system. Specifically, the author explores the changes in the equation of motion and the Josephson current when the superconducting island undergoes mechanical oscillation in an asymmetric geometric configuration.
### Main problem decomposition:
1. **Non - linear dynamic behavior**:
- The author focuses on how to describe the mechanical motion of the superconducting island and its influence on the Josephson current. In particular, under zero - bias and finite - bias conditions, the equation of motion of the superconducting island will exhibit different characteristics.
- In the case of zero - bias, the motion of the superconducting island is controlled by the Duffing equation, which depends on the superconducting phase difference. The Duffing equation is a non - linear differential equation widely used to describe the non - linear restoring force in small - scale resonators.
2. **Chaotic behavior**:
- When a bias is applied, the equation of motion of the system becomes a modified Duffing equation with time - dependent coefficients, which may lead to quasi - periodic or chaotic behavior. The author analyzed the system behavior under different bias conditions through numerical simulation and Poincaré mapping, and revealed the existence of chaotic solutions.
3. **Application prospects**:
- The research results are of great significance for the applications of nano - electro - mechanical systems (NEMS), such as high - sensitivity detection devices, quantum information processing, etc. In particular, the study of non - linear effects and chaotic behavior is helpful for developing new devices such as low - noise amplifiers and mechanical qubits.
### Mathematical models and formulas:
- **Duffing equation**:
\[
m_I \ddot{u} + (k_I + k_D (\cos(\omega_{J,R} t + \phi_R) - 1)) u - \frac{k_D}{2m_I R^2} (\cos(\omega_{J,R} t + \phi_R) - 1) u^3 = F_L(t) / m_I
\]
where \( m_I \) is the mass of the superconducting island, \( k_I \) is the spring constant, \( k_D=\frac{J_R \alpha}{e R} \), \( F_L(t)=-\frac{J_L \alpha}{e} (\cos(\omega_{J,L} t + \phi_L) - 1) \).
- **Simplified form under zero - bias**:
\[
\ddot{u} + A u - B u^3 = 0
\]
where \( A = \frac{k_I + k_D (\cos(\omega_{J,R} t + \phi_R) - 1)}{m_I} \), \( B=\frac{k_D}{2 m_I R^2} (\cos(\omega_{J,R} t + \phi_R) - 1) \).
- **Chaotic behavior under finite - bias**:
\[
m_I \ddot{u} + \gamma(t) \dot{u} + (k_I + k_D (\cos(\omega_{J,R} t + \phi_R) - 1)) u - \frac{k_D}{2m_I R^2} (\cos(\omega_{J,R} t + \phi_R) - 1) u^3 = F_L(t) / m_I
\]
where \( \gamma(t) \) is the time - dependent damping term.
Through these mathematical models, the author analyzed in detail the dynamic behavior of the system under different conditions and revealed the importance of non - linear and chaotic phenomena in nano - electro - mechanical systems.