Strang splitting schemes for N-level Bloch models
Numerical Analysis
2019-09-25 v1 Numerical Analysis
Computational Physics
Abstract
We extend to the N-level Bloch model the splitting scheme which use exact numerical solutions of sub-equations. These exact solutions involve matrix exponentials which we want to avoid to calculate at each time step. We use Newton interpolation to reduce the computational cost. The resulting scheme is nonstandard and preserves all qualitative properties of the Bloch equations. We show numerical simulations to compare this approach with a few other schemes.
Cite
@article{arxiv.1909.10932,
title = {Strang splitting schemes for N-level Bloch models},
author = {Marc Songolo and Brigitte Bidégaray-Fesquet},
journal= {arXiv preprint arXiv:1909.10932},
year = {2019}
}