Mathematica code:

R[n_] := (SeedRandom[n]; RandomReal[])
RE[n_] := (SeedRandom[n]; RandomReal[ExponentialDistribution[4]])

G[Q_, t_] :=
 Graphics[
  Table[
   {PointSize[ .06* (Sqrt[RE[Q*2 i^2] + RE[Q*2 i + 1]^2])],
    Opacity[.2 + .4*Sum[Exp[-(t - (R[Q*i] + j))^2/.1], {j, {-1, 0, 1}}]],
    Point[
       Table[
         RotationTransform[(k + (-1)^(i)*t)*2 Pi/8]
           [{RE[Q*2 i],RE[Q*2 i +1]}],
       {k, 0, 7, 1}]]},
 {i, 1, 90, 1}],
 PlotRange -> 1.5, ImageSize -> 500]

Table[
 ListAnimate[
  Table[
   G[Q,t],
  {t, 0, .95, .05}],
 AnimationRunning -> False],
{Q, {4, 8, 11, 13, 15, 18 }]