Abstract. This article presents the derivation of a semi-classical model of electromagnetic-wave propagation in a non centro-symmetric crystal. It consists of Maxwell's equations for the wave field coupled with a version of Bloch's equations which takes fully into account the discrete symmetry group of the crystal. The model is specialized in the case of a KDP crystal for which information about the dipolar moments at the Bloch level can be recovered from the macroscopic dispersion properties of the material.Mathematics Subject Classification. 78A60, 81V80.

We investigate the dynamics of a chain of oscillators coupled by fullynonlinear interaction potentials. This class of models includes Newton's cradle with Hertzian contact interactions between neighbors. By means of multiple-scale analysis, we give a rigorous asymptotic description of small amplitude solutions over large times. The envelope equation leading to approximate solutions is a discrete p-Schrödinger equation. Our results include the existence of long-lived breather solutions to the original model. For a large class of localized initial conditions, we also estimate the maximal decay of small amplitude solutions over long times.

In this paper, we present splitting schemes for the two-level Bloch model. After proposing two ways to split the Bloch equation, we show that it is possible in each case to generate exact numerical solutions of the obtained sub-equations. These exact solutions involve matrix exponentials which can be expensive to compute. Here, for [Formula: see text] matrices we develop equivalent formulations which reduce the computational cost. These splitting schemes are nonstandard ones and conserve all the physical properties (Hermicity, positiveness and trace) of Bloch equations. In addition, they are explicit, making effective their implementation when coupled with the Maxwell’s equations.