Solving spin quantum-master equations with matrix continued-fraction   methods: application to superparamagnets
  We implement continued-fraction techniques to solve exactly quantum master equations for a spin with arbitrary S coupled to a (bosonic) thermal bath. The full spin density matrix is obtained, so that along with relaxation and thermoactivation, coherent dynamics is included (precession, tunnel, etc.). The method is applied to study isotropic spins and spins in a bistable anisotropy potential (superparamagnets). We present examples of static response, the dynamical susceptibility including the contribution of the different relaxation modes, and of spin resonance in transverse fields. 
J. L. Garcia-Palacios and D. Zueco