Exponentiation of Parametric Hamiltonians via Unitary interpolation

The effort to generate matrix exponentials and associated differentials, required to determine the time evolution of quantum systems, frequently constrains the evaluation of problems in quantum control theory, variational circuit compilation, or Monte Carlo sampling. We introduce two ideas for the time-efficient approximation of matrix exponentials of linear multiparametric Hamiltonians. We modify the Suzuki-Trotter product formula from an approximation to an interpolation scheme to improve both accuracy and walltime. This allows us to achieve high fidelities within a single interpolation step, which can be computed directly from cached matrices. Furthermore, we define the interpolation on a grid of system parameters, and show that the interpolation infidelity converges with fourth-order accuracy in the number of interpolation bins.