Calculating the Hamiltonian for frequency-dependent squeezing

using SystemLevelHamiltonian

filter_cav = cavity(:filter)
sqz_cav = squeezing_cavity(:sqz)
ifo = radiation_pressure_cavity(:ifo)

sys = concatenate(:feqdepsqz,[sqz_cav,filter_cav,ifo])
sys = feedbackreduce(sys,:Out_sqz,:In_filter)
sys = feedbackreduce(sys,:Out_filter,:In_ifo)

SLH(:feqdepsqz, [:In_sqz], [:Out_ifo], [1.0;;], Vector{QuantumCumulants.QAdd}[[(sqrt(κ_filter)*(a_filter)+sqrt(κ_sqz)*(a_sqz)+sqrt(κ_ifo)*(a_ifo))]], ((0 + 1im)*ϵ_sqz*(a_sqz′*a_sqz′)+(0 - 1im)*ϵ_sqz*(a_sqz*a_sqz)+Δ_filter*(a_filter′*a_filter)+Δ_ifo*(a_ifo′*a_ifo)+Ω_ifo*(b_ifo′*b_ifo)+-g_ifo*(a_ifo′*a_ifo*b_ifo′)+-g_ifo*(a_ifo′*a_ifo*b_ifo)+(0.0 - 0.5im)*sqrt(κ_filter)*sqrt(κ_sqz)*(a_sqz*a_filter′)+(-0.0 + 0.5im)*sqrt(κ_filter)*sqrt(κ_sqz)*(a_sqz′*a_filter)+(0.0 - 0.5im)*sqrt(κ_filter)*sqrt(κ_ifo)*(a_filter*a_ifo′)+(0.0 - 0.5im)*sqrt(κ_ifo)*sqrt(κ_sqz)*(a_sqz*a_ifo′)+(-0.0 + 0.5im)*sqrt(κ_filter)*sqrt(κ_ifo)*(a_filter′*a_ifo)+(-0.0 + 0.5im)*sqrt(κ_ifo)*sqrt(κ_sqz)*(a_sqz′*a_ifo)))

sys.H

\[1 i \epsilon_{sqz} a_sqz^\dagger a_sqz^\dagger -1 i \epsilon_{sqz} a_{sqz} a_{sqz} + \Delta_{filter} a_filter^\dagger a_{filter} + \Delta_{ifo} a_ifo^\dagger a_{ifo} + \Omega_{ifo} b_ifo^\dagger b_{ifo} -1 g_{ifo} a_ifo^\dagger a_{ifo} b_ifo^\dagger -1 g_{ifo} a_ifo^\dagger a_{ifo} b_{ifo} -0.5 i \sqrt{\kappa_{filter}} \sqrt{\kappa_{sqz}} a_{sqz} a_filter^\dagger + 0.5 i \sqrt{\kappa_{filter}} \sqrt{\kappa_{sqz}} a_sqz^\dagger a_{filter} -0.5 i \sqrt{\kappa_{filter}} \sqrt{\kappa_{ifo}} a_{filter} a_ifo^\dagger -0.5 i \sqrt{\kappa_{ifo}} \sqrt{\kappa_{sqz}} a_{sqz} a_ifo^\dagger + 0.5 i \sqrt{\kappa_{filter}} \sqrt{\kappa_{ifo}} a_filter^\dagger a_{ifo} + 0.5 i \sqrt{\kappa_{ifo}} \sqrt{\kappa_{sqz}} a_sqz^\dagger a_{ifo}\]

operators(sys)

Set{Any} with 8 elements:
  b_ifo′
  a_ifo′
  a_ifo
  a_filter′
  a_sqz′
  a_sqz
  a_filter
  b_ifo

parameters(sys)

Set{Any} with 8 elements:
  κ_sqz
  ϵ_sqz
  g_ifo
  κ_ifo
  Ω_ifo
  Δ_filter
  κ_filter
  Δ_ifo