Yeah, I think there's an error. When calculating the vector of Rabi frequencies omega_vec , we multiply by fact, but that is set to zero for any Fock state for which n - sideband_index < 0 . But I think that implicitly assumes that we're starting in the internal ground state. If we started in the internal excited state, then the Rabi frequency for the Fock state |n = 0> should not be zero when driving a red sideband. We can probably fix this bug in different ways. One would be to correctly calculate omega_vec depending on which state we start in.
Error is here https://github.com/OxIonics/ionics_fits/blob/c5be47da9542e5d64423ef6b7f29cbe49ed27f3e/ionics_fits/models/laser_rabi.py#L193
From @AUTProgram