Displaced Harmonic Oscillator Model in QuTipr

[Eric]

Hello QuTip Users,

I want to simulate molecular electronic absorption with QuTip. Say we want to construct a shifted 1D harmonic oscillator excited state (3 eV above ground state, 65 fs period, dimensionless displacement of 2). I am interested exciting this system with 10 femtosecond gaussian pulses.

From the documentation, I suspect that the excited state potential can be represented as follows: displace(N,2).dag() * destroy(N).dag() * destroy(N) * displace(N,2) + excitation_energy. How does one build the full hamiltonian with the ground state, excited state, and optical couplings using the tensor command? Any guidance is much appreciated.

Best,

Eric