Multi-frequency averaging and uniform accuracy towards numerical approximations for a Bloch model

We are interested in numerically solving a transitional model derived from the Bloch model. The Bloch equation describes the time evolution of the density matrix of a quantum system forced by an electromagnetic wave. In a high frequency and low amplitude regime, it asymptotically reduces to a non-stiff rate equation. As a middle ground, the transitional model governs the diagonal part of the density matrix. It fits in a general setting of linear problems with a high-frequency quasi-periodic forcing and an exponentially decaying forcing. The numerical resolution of such problems is challenging. Adapting high-order averaging techniques to this setting, we separate the slow (rate) dynamics from the fast (oscillatory and decay) dynamics to derive a new micro-macro problem. We derive estimates for the size of the micro part of the decomposition, and of its time derivatives, showing that this new problem is non-stiff. As such, we may solve this micro-macro problem with uniform accuracy using standard numerical schemes. To validate this approach, we present numerical results first on a toy problem and then on the transitional Bloch model.


page 1

page 2

page 3

page 4


Uniformly accurate numerical schemes for a class of dissipative systems

We consider a class of relaxation problems mixing slow and fast variatio...

Uniformly accurate schemes for oscillatory stochastic differential equations

In this work, we adapt the micro-macro methodology to stochastic differe...

High-order implicit time integration scheme with controllable numerical dissipation based on mixed-order Padé expansions

A single-step high-order implicit time integration scheme with controlla...

A time multiscale decomposition in cyclic elasto-plasticity

For the numerical simulation of time-dependent problems, recent works su...

Higher-order generalized-α methods for hyperbolic problems

The generalized-α time-marching method provides second-order accuracy in...

Numerical geometric acoustics: an eikonal-based approach for modeling sound propagation in 3D environments

We present algorithms for solving high-frequency acoustic scattering pro...

High-order quasi-Helmholtz Projectors: Definition, Analyses, Algorithms

The accuracy of the electric field integral equation (EFIE) can be subst...

Please sign up or login with your details

Forgot password? Click here to reset