https://github.com/scqubits/scqubits-examples/blob/master/examples/demo_legacy_explorer.ipynb

[ahkatlio]

Define the interaction term

interaction_term = InteractionTerm(
g_strength=0.2,
operator_list=[(0, qbt.n_operator()), (1, osc.creation_operator() + osc.annihilation_operator())]
)

Define the subsystem list

subsystem_list = [qbt, osc]

Get the Hamiltonian

hamiltonian = interaction_term.hamiltonian(subsystem_list)

Print the Hamiltonian

print(hamiltonian)