In the post-bounce phase (\(t > t_B\) ), from Eqs.(REF )-(REF ) we find that [1]},
\(H^2 &=&\frac{8\pi G \rho }{3}\left(1-\frac{\rho }{\rho _c^{\scriptscriptstyle {\mathrm {I}}}}\right)\Bigg [1 +\frac{\gamma ^2}{\gamma ^2+1}\left(\frac{\sqrt{\rho /\rho _c^{\scriptscriptstyle {\mathrm {I}}}}}{1 +\sqrt{1-\rho /\rho _c^{\scriptscriptstyle {\mathrm {I}}}}}\right)^2\Bigg ], \\\frac{\ddot{a}}{a} &=&-\frac{4\pi G}{3}\left(\rho + 3P\right)+ \frac{4\pi G \rho }{3}\left[\frac{\left(7\gamma ^2+ 8\right) -4\rho /\rho _c^{\scriptscriptstyle {\mathrm {I}}} +\left(5\gamma ^214+8\right)\sqrt{1-\rho /\rho _c^{\scriptscriptstyle {\mathrm {I}}}}}{(\gamma ^2 +1)\left(1+\sqrt{1-\rho /\rho _c^{\scriptscriptstyle {\mathrm {I}}}}\right)^2}\right]\frac{\rho }{\rho _c^{\scriptscriptstyle {\mathrm {I}}}}\nonumber \\&& + 4\pi G P \left[\frac{3\gamma ^2+2+2\sqrt{1-\rho /\rho _c^{\scriptscriptstyle {\mathrm {I}}}}}{(\gamma ^2+1)\left(1+\sqrt{1-\rho /\rho _c^{\scriptscriptstyle {\mathrm {I}}}}\right)}\right]\frac{\rho }{\rho _c^{\scriptscriptstyle {\mathrm {I}}}},\) 
