Introducing physical relaxation terms in Bloch equations

Positiveness conservation with the splitting scheme

Abstract

The Bloch equation models the evolution of the state of electrons in matter described by a Hamiltonian. To model more physical phenomena we have to introduce phenomenological relaxation terms. The introduction of these terms has to conserve some positiveness properties. The aim of this paper is to review possible relaxation models and to provide insight into how to discretize them properly in view of numerical computations.

Publication
Journal of Computational Physics, 170(2), 603–613 (2001)