SOL-2 THEORY: iterative dQ and Hessian calculations added