                                            2                                  ↪
H = ω₁⋅ad₁⋅a₁ + ω₂⋅ad₂⋅a₂ + EC₃₃⋅(-ng₃ + n₃)  + EL₁⋅(-φ₁ - φ₂)⋅-φₑ₁/2 + EL₂⋅(- ↪

↪                                                       
↪ φ₁ + φ₂)⋅-φₑ₁/2 - EJ₁⋅cos(-φ₃ + φ₁) - EJ₂⋅cos(φ₁ + φ₃)
------------------------------------------------------------
mode 1:       harmonic       φ₁ = zp₁⋅(a₁ + ad₁)       ω₁/2⋅π = 3.22489       zp₁ = 2.49
mode 2:       harmonic       φ₂ = zp₂⋅(a₂ + ad₂)       ω₂/2⋅π = 0.39497       zp₂ = 0.872
mode 3:       charge       ng₃ = 0
------------------------------------------------------------
parameters:       EC₃₃ = 0.296       EL₁ = 0.13       EL₂ = 0.13       EJ₁ = 5.0       EJ₂ = 5.0
loops:       φₑ₁/2⋅π = 0.5