Optimal control of interacting quantum systems based on the first-order Magnus approximation: Application to multiple dipole-dipole coupled molecular rotors

We develop a methodology for performing approximate optimal control simulations for quantum systems with multiple interacting degrees of freedom. The quantum dynamics are modeled using the first-order Magnus approximation in the interaction picture, where the interactions between different degrees of freedom are treated as the perturbation. We present a numerical procedure for implementing this approximation, which leverages the separability of the zeroth-order time evolution operator and the pairwise nature of common interactions for a reduced computational cost. This formulation of the first-order Magnus approximation is suitable to be combined with gradient-free methods for control field optimization; to this end, we adopt a Stochastic Hill Climbing algorithm. The associated computational costs are analyzed and compared with those of the exact simulation in the large N limit. For numerical illustrations, we perform approximate optimal control simulations for systems of two and three dipole-dipole coupled molecular rotors under the influence of a global control field. For the two rotor system, we optimize fields for both orientation and entanglement control objectives. For the three rotor system, we optimize fields for orienting rotors in the same direction as well as in opposite directions.