“…The evolution over a short time interval (of extension δ t ) of the nuclear wave packet Ψ H2 ( R , t ) representing the H 2 molecule is then performed numerically in two successive steps. In the first step, the population loss due to the ionization of the neutral molecule is taken into account by introducing the instantaneous ionization rate W H2 ( R , t ), according to Ψ̅ normalH 2 ( R , t ) = false[ 1 − W normalH 2 ( R , t ) δ t false] 1 / 2 Ψ normalH 2 ( R , t ) To evaluate the ionization rate of H 2 in an intense laser field, we use the molecular ionization model , derived from the PPT approach, − as described in detail in Vigneau et al for instance. Finally, in a second step, the resulting wave packet Ψ̅ normalH 2 ( R , t ) is propagated numerically according to Ψ normalH 2 ( R , t + δ t ) = Û normalH 2 ( t + δ t ← t ) Ψ̅ normalH 2 ( R , t ) with the evolution operator Û normalH 2 ( t + δ ...…”