Dans un contexte où les sciences lasers évoluent rapidement, il est essentiel de pouvoir résoudre l’équation de Schödinger pour des systèmes multiélectroniques dans le régime non perturbatif. Nous proposons un algorithme de dynamique multi-électronique corrélée qui repose sur une partition de type Feshbach. La forme Adams de cette partition combinée avec un propagateur Cayley-Crank-Nicholson permet à cet algorithme de conserver la norme de la fonction d’onde et d’avoir un bon comportement de convergence. Nous détaillons la mise en oeuvre de cette méthodologie pour la dynamique d’ionisation en champ intense, d’abord au niveau mono-électronique, puis au niveau de N-électrons. Afin d’illustrer cette méthodologie, nous avons étudié l’ionisation de la molécule H₂ à noyaux fixes lorsqu’elle est soumise à un champ laser XUV et attoseconde.