In order to simulate the time-dynamics during the onset of the drive pulse, the system was also modeled numerically using QuTiP [1]}, [2]}. The model, accounting for multiple qubits, is a modification of the model described in [3]} to add the direct 'dipole-dipole' coupling term between the qubits \(\mathcal {G}_{ij}\) :
\(H =H_{qubits}+H_{qubit-drive}+H_{qubit-qubit}\) 
\(H_{qubits} = \sum _{i=1}^2\omega _{i} \left|e_i\right> \left<e_i\right|\) 
\(H_{qubit-drive} = \sum _{i=1}^2\left(\epsilon _i \sigma _i^{+}+h.c.\right)\) 
\(H_{\textnormal {qubit-qubit}} = \sum _{j> i=1}^2\left( J_{ij} + \mathcal {G}_{ij}\right) \left( \sigma _i^- \sigma _j^+ + \sigma _i^+\sigma _j^-\right)\) 
\(\dot{\rho }= -\frac{i}{\hbar }\left[H,\rho \right] + \sum _{i,j=1}^2 \gamma _{ij}\left[ \sigma _i^- \rho \sigma _j^+ -\frac{1}{2}\lbrace \sigma _i^+\sigma _j^-, \rho \rbrace \right].\) 
