Parameter differentiation and quantum state decomposition for time varying Schrödinger equations

For the unitary operator, solution of the Schrödinger equation corresponding to a time-varying Hamiltonian, the relation between the Magnus and the product of exponentials expansions can be expressed in terms of a system of first-order differential equations in the parameters of the two expansions. A method is proposed to compute such differential equations explicitly and in a closed form.