We propose an iterative variational quantum algorithm to simulate the time evolution of arbitrary initial states within a given subspace. The algorithm compresses the Trotter circuit into a shorter-depth parameterized circuit, which is optimized simultaneously over multiple initial states in a single training process using fidelity-based cost functions. After the whole training procedure, we provide an efficiently computable lower bound on the fidelities for arbitrary states within the subspace, which guarantees the performance of the algorithm in the worst-case training scenario. We also show our cost function exhibits a barren-plateau-free region near the initial parameters at each iteration in the training landscape. The experimental demonstration of the algorithm is presented through the simulation of a 2-qubit Ising model on an IBMQ processor. As a demonstration for a larger system, a simulation of a 10-qubit Ising model is also provided.
Park et al. (Mon,) studied this question.