Krylov-subspace approach for the efficient control of quantum many-body dynamics

Abstract

The gradient ascent pulse engineering (GRAPE) algorithm is a celebrated control algorithm with excellent converging rates, owing to a piecewise-constant ansatz for the control function that allows for cheap objective gradients. However, the computational effort involved in the exact simulation of quantum dynamics quickly becomes a bottleneck limiting the control of large systems. In this paper, we experiment with a modified version of GRAPE that uses Krylov approximations (K-GRAPE) to deal efficiently with high-dimensional state spaces. Even though the number of parameters required by an arbitrary control task scale linearly with the dimension of the system, we find a constant elementary computational effort (the effort per parameter). Since the elementary effort of GRAPE is superquadratic, this speed up allows us to reach dimensions far beyond. The performance of the K-GRAPE algorithm is benchmarked in the paradigmatic XXZ spin-chain model.