The general potential of power-law inflation is as $ V(\phi)\propto \phi^{n} $ with scalar field $\phi$. The behavior of inflation is often known as power-law expansion like $S(\eta)\propto \eta^{1+\beta}$ with $1+\beta<0$. In this paper, the theoretical spectra of relic gravitational waves are compared with the measured strain sensitivity of Advanced LIGO and VIRGO, corresponding to the latest detected events of gravitational waves. The results show tight constraints on $\beta$ and $n$. Also, the obtained constraints indicate that special types of the potential of inflation, prototype, and KKLTI models, which are originated from string theory, are more suitable than other models. Our analysis shows that there exist some more chances for the detection of relic gravitational waves.

The potential $V(\phi)=\lambda \phi^{n}$ is responsible for the inflation of the universe as scalar field $\phi$ oscillates quickly around some point where $V(\phi)$ has a minimum. The end of this stage has an important role on the further evolution stages of the universe. The created particles are responsible for reheating the universe at the end of this stage. The behaviour of the inflation and reheating stages are often known as power law expansion $S(\eta) \varpropto \eta^{1+\beta}$, $S(\eta) \varpropto \eta^{1+\beta_{s}}$ respectively. The reheating temperature ($T_{rh}$) and $\beta_{s}$ give us valuable information about the reheating stage. Recently people have studied about the behaviour of $T_{rh}$ based on slow-roll inflation and initial condition of quantum normalization. It is shown that there is some discrepancy on $T_{rh}$ due to amount of $\beta_{s}$ under the condition of slow-roll inflation and quantum normalization \cite{acq}. Therefore the author is believed in \cite{acq} that the quantum normalization may not be a good initial condition. But it seems that we can remove this discrepancy by determining the appropriate parameter $\beta_{s}$ and hence the obtained temperatures based on the calculated $\beta_{s}$ are in favour of both mentioned conditions. Then from given $\beta_{s}$, we can calculate $T_{rh}$, tensor to scalar ratio $r$ and parameters $\beta, n$ based on the Planck and WMAP-9 data. The observed results of $r, \beta_{s}, \beta$ and $n$ have consistency with their constrains. Also the results of $T_{rh}$ are in agreement with its general range and special range based on the DECIGO and BBO detectors.

The thermal spectrum of relic gravitational waves causes the new amplitude that called `modified amplitude'. Our analysis shows that, there exist some chances for detection of the thermal spectrum in addition to the usual spectrum by Adv.LIGO data of GW150914 and detector based on the maser light (Dml). The behavior of the inflation and reheating stages are often known as power law expansion like $ S(\eta)\propto \eta^{1+\beta}$, $S(\eta)\propto \eta^{1+\beta_s}$ respectively. The $\beta$ and $\beta_s$ have an unique effect on the shape of the spectrum. We find some upper bounds on the $\beta$ and $\beta_s$ by comparison the usual and the thermal spectrum with the Adv.LIGO and Dml. As this result modifies our information about the nature of the evolution of inflation and reheating stages.

The spectrum of thermal gravitational waves is obtained by including the high frequency thermal gravitons created from extra-dimensional effect and is a new feature of the spectrum. The amplitude and spectral energy density of gravitational waves in thermal vacuum state are found enhanced. The amplitude of the waves get modified in the frequency range (10$^{-16}$ -10 $^{8}$ Hz) but the corresponding spectral energy density is less than the upper bound of various estimated results. With the addition of higher frequency thermal waves, the obtained spectral energy density of the wave in thermal vacuum state does not exceed the upper bound put by nucleosynthesis rate. The existence of cosmologically originated thermal gravitational waves due to extra dimension is not ruled out.