and the lack of numerical simulations, the details of relaxation processes and virialization in MOND are much less known. Thus, before discussing the specific properties of the collapse end-products we present a general overview of the time evolution of the virial quantities in our simulations, postponing to Section \[secphsp\] a more detailed description of the phase-space evolution. In particular, in Fig. \[figtime\] we show the time evolution of $2 K/|W|$, $K$, $W$, and $K+W$ for simulations D, M1, and N. In the diagrams time is normalized to ${t_{\rm dyn}}$, so plots referring to different simulations are directly comparable (the values of ${t_{\rm dyn}}$ in time units for the five simulations are given in Section \[secic\]). In simulation N (right column) we find the well known behavior of Newtonian dissipationless collapses: $2 K/|W|$ has a peak, then oscillates, and eventually converges to the equilibrium value $2 K/|W|=1$; the total energy $K+W$ is nicely conserved during the collapse, though it presents a secular drift, a well known feature of time integration in N-body codes. The time evolution of the same quantities is significantly different in a dMOND simulation (left column). In particular, the virial ratio $2 K/|W|$ quickly becomes close to one, 