I typically work in Matlab, but recently I've been considering switching to QuTiP and I'm curious if it's easy to simulate a certain type of problem in QuTiP.
I know QuTiP can handle time-dependent Hamiltonians, but can it handle pulses of light entering a ensemble of atoms? Such systems that are coupled to an external field that follows equations of motion unrelated to the system's Hamiltonian. (In my specific case, the atoms evolve with the Optical-Bloch equations, but the field that drives the atoms evolves with Maxwell-Schrodinger equations which are coupled to the OBE.)
For a pulse of light enter an atomic medium, as the pulse updates density matrices along its propagation path, the buildup of coherence influences the pulse itself. So not only is your Hamiltonian time-dependent due to the rabi-frequencies evolving in time, but you have to update your E-field during propagation as it travels through the medium.
This is a pretty common type of problem solved in quantum optics, so I can imagine someone else has tried to do such a thing in QuTip, but I haven't had any success finding anything on it. Can anyone point me to an example of a free-space simulation of a pulse of light entering an ensemble using QuTiP? And if such a thing hasn't been solved before, is it worth my time looking into using the toolbox for help or just coding it from scratch?