I'm studying what happens to a wave packet if it propagates in a medium under linear optics conditions. In my equations I have the wave vector $k=\frac{\omega}{c}n(\omega)$ (often called $\beta$) which is developed in Taylor series as:
$\beta(\omega)\simeq\beta(\omega_0)+\left.\frac{d\beta}{d\omega}\right|_{\omega_0}(\omega-\omega_0)+\left.\frac{d^2\beta}{d\omega^2}\right|_{\omega_0}(\omega-\omega_0)^2=\beta(\omega_0)+\beta_1(\omega-\omega_0)+\beta_2(\omega-\omega_0)^2$
If we consider only $\beta_1$, we get:
$\beta_1=\left.\left[\frac{d}{d\omega}\frac{\omega}{c}n(\omega)\right]\right|_{\omega_0}=\frac{n(\omega_0)}{c}+\frac{\omega_0}{c}\left.\frac{dn}{d\omega}\right|_{\omega_0}$
where the first element is clearly the inverse of the phase velocity $v_p$ of a wave having frequency $\omega_0$, the second one is a "corrective" term and all together gives the inverse of the group velocity $\beta_1=\frac{1}{v_g}$. From RP Photonics "The first-order term contains the inverse group velocity (i.e., the group delay per unit length) and describes an overall time delay without an effect on the pulse shape". So, if I suppose $\beta_2=0$ I get a travelling pulse that won't be subjected to time broadening. It is also clear that, if $n$ does not depend on frequency (like in vacuum) we get $v_p=v_g$; vice versa, if $n=n(\omega) \Rightarrow v_g\neq v_p$, so the shape will move at a different velocity with respect its phase.
Now, my question is the following:
Suppose we have the second case, that is a Gaussian pulse whose envelope moves at a different group velocity with respect to the phase velocity as shown at this link: (https://en.wikipedia.org/wiki/File:Wave_packet_propagation_(phase_faster_than_group,_nondispersive).gif). For this to happen I have to suppose that $n$ depends on frequency so, from the previous equation, the two velocities are different. But if $n=n(\omega)$, this means that each wave composing the Gaussian packet will move at a different speed as shown in Figure 1 at this link (https://www.rp-photonics.com/group_velocity.html) where the phase fronts of different frequency components propagate with different velocities faster than group velocity and the shape of the pulse isn't broadening. From my book I read "The group velocity dispersion (GVD) is given by the dependence of $v_g$ from $\omega$: inside the envelope, different frequencies moves at different velocity and this is the cause of the time broadening of the pulse".
If velocities are different, why the pulse isn't broadening? And, if different velocities are possible without having a time broadening of the pulse, what is the physical cause of that broadening?
EDIT:
After a lot of research I have cleared up some of the doubts concerning the problem I mentioned above. I was helped a lot by this PDF (link: http://sharif.edu/~kavehvash/Group_Phase_Velocity.pdf) which explains what happens to the phase and group velocities according to the dispersion relation. In particular, two cases are distinguished:
The dispersion relation is linear: this causes group velocity and phase velocity to be equal. This happens because in this condition it should be $n=n_0$ and from the formula of $\beta_1$ we get $v_g=v_p$.
The dispersion relation is non-linear: in this case, the two velocities are always different from each other (unless you have special restrictive conditions). This is due to the fact that $n=n(\omega)$ and each component moves at a different speed. From the PDF: "A very important consequence of this is that our initial wave package broadens out with time because the partial waves forming it gradually move out of phase with each other" and this because, if $v_p=v_p(\omega)$, then $v_g=v_g(\omega)$. The only alternative for which $v_g$ does not depend on $\omega$ and does not broaden is that it is $0$, i.e. there is a standing wave.
If what I have written is correct I am left with the last doubt: is it possible to have $v_g\neq v_p$ without obtaining a time broadening as shown in the following two links?
Link2 (Figure 1): https://www.rp-photonics.com/group_velocity.html
If yes, how?