standard propagating gravitational waves with nonzero friction term. The widely discussed inequality $D_L^{({\rm gw})}\neq D_L^{({\rm em})}$ induced by the friction term inspire our work. After quantizing the fields under slow expansion approximation (equivalently short wavelength approximation), we point out the quantum nature of $D_L^{({\rm gw})}\neq D_L^{({\rm em})}$ is the particle number non-conservation of graviton. For the inflation with de Sitter background, we obtain the full solution of the perturbation equation, and quantitatively calculate the influence of the friction term on the spectrum index.

Acknowledgements {#acknowledgements.unnumbered}
================

This work was supported by the National Natural Science Foundation of China under Grants No. 11633001 and the Strategic Priority Research Program of the Chinese Academy of Sciences, Grant No. XDB23000000.

Appendix: Whittaker function {#appendix-whittaker-function.unnumbered}
============================

Definition of the Whittaker equation and functions {#sec:AppA}
--------------------------------------------------

As in [@Gradshteyn2007], the Whittaker equation reads $$\frac{{{\rm d}}^2W}{{{\rm d}}z^2}+\left(-\frac{1}{4}+\frac{\lambda}{z}+\frac{1/4-\mu^2}{z^2}\right)W=0,$$ and two linear independent solutions are $$\begin{aligned}
  M_{\lambda,\mu}(z)&=z^{\mu+\frac{1}{2}}e^{-z/2}{}_1F_1(\mu-\lambda+\frac{1}{2};2\mu+1;z),\label{eq:39}\\
  W_{\lambda,\mu}(z)&=\frac{\Gamma(-2\mu)}{\Gamma(\frac{1}{2}-\mu-\lambda)}M_{\lambda,\mu}(z)\nonumber\\
  &\qquad\quad+\frac{\Gamma(2\mu)}{\Gamma(\frac{1}{2}+\mu-\lambda)}M_{\lambda,-\mu}(z),\label{eq:40}\end{aligned}$$ where the confluent hypergeometric function $$\label{eq:41}
  {}_1F_1(\alpha;\gamma;z)=1+\frac{\alpha}{\gamma}\frac{z}{1!}+\frac{\alpha(\alpha+1)}{\gamma(\gamma+1)}\frac{z^
standard propagating gravitational waves with nonzero friction term. The wide discourse inequality $ D_L^{({\rm gw})}\neq D_L^{({\rm em})}$ induced by the friction condition cheer our work. After quantizing the field under slow expansion approximation (equivalently short wavelength approximation), we point out the quantum nature of $ D_L^{({\rm gw})}\neq D_L^{({\rm em})}$ is the atom number non - conservation of graviton. For the inflation with de Sitter background, we obtain the wide solution of the perturbation equation, and quantitatively calculate the influence of the friction term on the spectrum exponent. 

 Acknowledgements { # acknowledgements.unnumbered } 
 = = = = = = = = = = = = = = = = 

 This work was supported by the National Natural Science Foundation of China under Grants No. 11633001 and the Strategic Priority Research Program of the Chinese Academy of Sciences, Grant No. XDB23000000. 

 Appendix: Whittaker routine { # appendix-whittaker-function.unnumbered } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 Definition of the Whittaker equation and functions { # sec: AppA } 
 -------------------------------------------------- 

 As in [ @Gradshteyn2007 ], the Whittaker equality reads $ $ \frac{{{\rm d}}^2W}{{{\rm d}}z^2}+\left(-\frac{1}{4}+\frac{\lambda}{z}+\frac{1/4-\mu^2}{z^2}\right)W=0,$$ and two analogue independent solutions are $ $ \begin{aligned } 
   M_{\lambda,\mu}(z)&=z^{\mu+\frac{1}{2}}e^{-z/2}{}_1F_1(\mu-\lambda+\frac{1}{2};2\mu+1;z),\label{eq:39}\\ 
   W_{\lambda,\mu}(z)&=\frac{\Gamma(-2\mu)}{\Gamma(\frac{1}{2}-\mu-\lambda)}M_{\lambda,\mu}(z)\nonumber\\ 
   & \qquad\quad+\frac{\Gamma(2\mu)}{\Gamma(\frac{1}{2}+\mu-\lambda)}M_{\lambda,-\mu}(z),\label{eq:40}\end{aligned}$$ where the confluent hypergeometric function $ $ \label{eq:41 } 
   { } _ 1F_1(\alpha;\gamma;z)=1+\frac{\alpha}{\gamma}\frac{z}{1!}+\frac{\alpha(\alpha+1)}{\gamma(\gamma+1)}\frac{z^
stajdard propagating gravitxtional waves wnrh nonvero frjction tdrm. The widely discussed ineeuqlity $D_L^{({\rm gw})}\neq D_L^{({\rm em})}$ inauced by nhe frictuon uerm inspire our xkrk. Aftcx quahbiziny vhe fields undet slow expanvion approximadiun (equivalently short wavelength apprjximatipn), we point out jhe qlagtum nature of $D_L^{({\rm gw})}\neq D_L^{({\rm em})}$ is tge partpcle number non-comservation of graviton. For the inflation with de Sitter bacjgrotbd, we obtain the full solution of jhe perturbation equation, and quxntitctively calxuoatf the influeice of the frictiok term mn the xpectrum index.

Sckiowlwdgements {#acknowledgekents.unnumbered}
================

Thif work wav aupported by the Barionan Nadurau Sckende Fkundatlon of China hnder Grantw No. 11633001 and the Strattgis Priority Ressarch [rjgram of the Chinese Academy of Sciencev, Gdant No. XDB23000000.

Appendix: Whirtaker function {#appenfix-whittaher-function.unnumbered}
============================

Definition of the Whittaker aquatmov aue funzrilns {#sec:AppA}
--------------------------------------------------

As in [@Gradshteyn2007], the Whittaker equwfipn reads $$\frac{{{\rm b}}^2W}{{{\rm d}}z^2}+\left(-\frac{1}{4}+\ftaf{\lsibda}{z}+\frac{1/4-\mu^2}{z^2}\rkght)W=0,$$ chd two linear indepejdent sjlutiins are $$\btgin{akigned}
  M_{\lambda,\mu}(z)&=z^{\mu+\frac{1}{2}}e^{-z/2}{}_1D_1(\mu-\lambda+\frab{1}{2};2\mu+1;z),\oabel{eq:39}\\
  W_{\lambda,\mu}(v)&=\frac{\Gamma(-2\mu)}{\Yamma(\ftac{1}{2}-\mu-\lsmbda)}M_{\lambda,\mu}(z)\nonumber\\
  &\qqhad\quad+\frac{\Hamma(2\mu)}{\Gajoa(\frac{1}{2}+\mu-\lambda)}M_{\lxmbca,-\ku}(z),\label{td:40}\end{aligned}$$ where the confouenc hypergdomejric fugction $$\labfl{eq:41}
  {}_1F_1(\alpha;\gamma;z)=1+\frac{\aloha}{\galmd}\frac{z}{1!}+\frac{\wlpha(\alpha+1)}{\gamma(\gamma+1)}\frac{z^
standard propagating gravitational waves with nonzero friction widely inequality $D_L^{({\rm D_L^{({\rm em})}$ induced our After quantizing the under slow expansion (equivalently short wavelength approximation), we point the quantum nature of $D_L^{({\rm gw})}\neq D_L^{({\rm em})}$ is the particle number non-conservation graviton. For the inflation with de Sitter background, we obtain the full solution the equation, quantitatively the influence of the friction term on the spectrum index. Acknowledgements {#acknowledgements.unnumbered} ================ This work was by the National Natural Science Foundation of China Grants No. 11633001 and Strategic Priority Research Program of Chinese of Sciences, No. Appendix: function {#appendix-whittaker-function.unnumbered} ============================ of the Whittaker equation and functions {#sec:AppA} -------------------------------------------------- As in [@Gradshteyn2007], the Whittaker equation reads $$\frac{{{\rm d}}^2W}{{{\rm and two solutions are M_{\lambda,\mu}(z)&=z^{\mu+\frac{1}{2}}e^{-z/2}{}_1F_1(\mu-\lambda+\frac{1}{2};2\mu+1;z),\label{eq:39}\\ &\qquad\quad+\frac{\Gamma(2\mu)}{\Gamma(\frac{1}{2}+\mu-\lambda)}M_{\lambda,-\mu}(z),\label{eq:40}\end{aligned}$$ the confluent hypergeometric {}_1F_1(\alpha;\gamma;z)=1+\frac{\alpha}{\gamma}\frac{z}{1!}+\frac{\alpha(\alpha+1)}{\gamma(\gamma+1)}\frac{z^
standard propagating gravitAtional wavEs witH noNzeRo FricTion Term. The widely dIScusSed inequality $D_L^{({\rm gw})}\neq d_L^{({\rm eM})}$ iNDuceD By The frIction tERm INSpiRe OuR woRk. aFtEr quaNtiZing the Fields undeR slOw Expansion appROxImation (equIvaLently short wAveLength ApProXImatiOn), wE poinT out thE QuantuM nature of $d_L^{({\RM gw})}\neq d_l^{({\rm em})}$ is THE pArtiCle number non-conseRVaTIon of graviton. FOr the iNfLAtION wiTh dE Sitter bacKgRound, WE obtain THe FULL soLUtion of the perTurbation eqUAtiOn, and qUaNtiTAtivelY calcUlATe tHe influence Of thE friction Term on THe spectRUm index.

acknowLedGemEnts {#ACkNoWleDgEMenTS.uNnuMBerEd}
================

This woRk WaS suppOrteD BY THe NaTioNal NAturaL Science FoundAtiOn of cHinA undeR GranTs No. 11633001 AnD the STrategIc PriOrIty Research ProgRam oF the ChineSe ACaDemY oF ScieNCes, GraNt NO. XDb23000000.

AppendIx: WhittAKer FuNCTIoN {#appendix-whittaker-FuNCTiOn.unnumbEred}
============================

DeFInItIOn of the WHiTtaKer eQUAtion And fUNcTions {#sec:appA}
--------------------------------------------------

As IN [@GRaDshteyn2007], ThE WhittAkEr eQuaTion rEAds $$\fRac{{{\rm d}}^2w}{{{\rm d}}z^2}+\lefT(-\frac{1}{4}+\FRac{\lambda}{z}+\frac{1/4-\MU^2}{z^2}\right)W=0,$$ and twO LiNEAr INdepEndEnt solutionS are $$\BEgin{AligNEd}
  m_{\laMBda,\mu}(Z)&=z^{\mu+\fRaC{1}{2}}E^{-z/2}{}_1f_1(\Mu-\lambda+\frac{1}{2};2\mu+1;z),\labeL{eQ:39}\\
  W_{\lambDa,\mu}(z)&=\Frac{\Gamma(-2\mu)}{\GaMma(\frac{1}{2}-\mu-\lAMBDa)}M_{\lambdA,\mu}(z)\NOnUMber\\
  &\qquad\quad+\fRac{\GaMma(2\mu)}{\Gamma(\FRac{1}{2}+\mu-\lamBda)}M_{\lAmbda,-\mu}(z),\Label{eq:40}\enD{ALigned}$$ whEre The ConFluENT hYpergeometric FUNctiOn $$\Label{eq:41}
  {}_1f_1(\alPha;\gammA;z)=1+\fRac{\AlpHa}{\gAmMa}\frac{z}{1!}+\frAc{\alpha(\aLpHa+1)}{\GaMmA(\gaMma+1)}\frAC{z^
standard propagating gravi tational w aveswit h n on zero  fri ction term. Th e  wid ely discussed inequali ty $D _L ^ {({\ r mgw})} \neq D_ L ^{ ( { \rm  e m} )}$  i n du ced b y t he fric tion termins pi re our work. Af ter quanti zin g the fields  un der sl ow  ex p ansio n a pprox imatio n  (equi valentlysh o rt wav e lengtha p pr oxim ation), we pointo ut the quantum na ture o f$ D_ L ^ {({ \rm  gw})}\neq  D _L^{( { \rm em} ) }$ i s  th e  particle num ber non-con s erv ationof  gr a viton.  Forth e  in flation wit h de  Sitter b ackgro u nd, weo btain t he ful l s olu tion of  t hepe r tur b at ion equ ation, a nd  q uanti tati v e l y  cal cul atethe i nfluence of t hefric t ion  term  on t he s pe ctrum  index .

Ac kn owledgements {# ackn owledgeme nts .u nnu mb ered} ====== === === ====

T his wor k  wa ss u p po rted by the Nation al N at ural Sci ence F o un da t ion of C hi naunde r Grant s No .  1 1633001and th e  S tr ategicPr iority  R ese arc h Pro g ramof the  Chinese  Acad e my of Sciences ,  Grant No. XD B 23 0 0 00 0 0.
App endix: Whit take r  fun ctio n  { #ap p endix -whit ta k er - function.unnumbered }====== ===== ============= ====

Defi n i t ion of t he W h it t aker equationand f unctions { # sec:AppA }
--- -------- --------- - - -------- --- --- --- --- - - -- ----

As in [ @ G rads ht eyn2007 ],the Whi tta ker  eq uat io n reads $ $\frac{{ {\ rm  d }} ^2W }{{{\ r m d}}z^2 }+ \le ft (-\ frac{ 1 }{4}+\ frac{ \lam bd a} { z}+ \frac{1 / 4- \ m u^2} {z ^2 }\ri ght )W =0,$$  and two  linear  independ ent solu ti on s are $ $\begin{align ed }
  M_{\la mb da, \mu}(z ) & =z^{\mu+ \frac{1}{2}}e^{-z/2}{}_ 1 F_1(\mu -\l ambda +\fr ac{1}{2}; 2\m u+1;z) ,\l a bel{eq :39}\\ 
  W_ {\ lam b d a,\mu } ( z) &=\ fr ac{\Gamma( - 2 \mu )}{\G am ma(\ frac{1} {2}-\mu-\lambda)}M _ {\l ambda,\mu}(z) \no numb e r \\ 
 & \q q uad \q u ad+ \ f rac{\Gamma(2\mu )}{\Gamma( \f r ac {1}{2}+\mu - \la mb da)}M_{ \lambda ,-\mu } (z),\la bel{eq:40 }\end{ali gn ed}$ $ whe re the con fluent h ypergeome t ric f u nc tion$$\ label{ eq :41 }
  { }_1F_1 ( \al pha;\ gamma; z) =1+\fr ac{\a lp ha}{\gam ma}\frac{z}{1!}+\frac{\ alpha( \alph a+1 )}{\gamma (\g a mma +1)}\frac {z^
standard propagating_gravitational waves_with nonzero friction term._The widely_discussed_inequality $D_L^{({\rm_gw})}\neq_D_L^{({\rm em})}$ induced_by the friction_term inspire our work._After quantizing the_fields_under slow expansion approximation (equivalently short wavelength approximation), we point out the quantum nature_of_$D_L^{({\rm gw})}\neq_D_L^{({\rm_em})}$_is the particle number non-conservation_of graviton. For the inflation_with de_Sitter background, we obtain the full solution of_the_perturbation equation, and_quantitatively calculate the influence of the friction term on_the spectrum index.

Acknowledgements {#acknowledgements.unnumbered}
================

This work was_supported by the_National_Natural_Science Foundation of China_under Grants No. 11633001 and the_Strategic Priority Research Program of the_Chinese Academy of Sciences, Grant No. XDB23000000.

Appendix:_Whittaker function {#appendix-whittaker-function.unnumbered}
============================

Definition of the Whittaker_equation and functions {#sec:AppA}
--------------------------------------------------

As in_[@Gradshteyn2007], the_Whittaker equation reads $$\frac{{{\rm d}}^2W}{{{\rm_d}}z^2}+\left(-\frac{1}{4}+\frac{\lambda}{z}+\frac{1/4-\mu^2}{z^2}\right)W=0,$$ and two_linear independent_solutions are $$\begin{aligned}
_ M_{\lambda,\mu}(z)&=z^{\mu+\frac{1}{2}}e^{-z/2}{}_1F_1(\mu-\lambda+\frac{1}{2};2\mu+1;z),\label{eq:39}\\
  W_{\lambda,\mu}(z)&=\frac{\Gamma(-2\mu)}{\Gamma(\frac{1}{2}-\mu-\lambda)}M_{\lambda,\mu}(z)\nonumber\\
  &\qquad\quad+\frac{\Gamma(2\mu)}{\Gamma(\frac{1}{2}+\mu-\lambda)}M_{\lambda,-\mu}(z),\label{eq:40}\end{aligned}$$_where the confluent_hypergeometric function $$\label{eq:41}
  {}_1F_1(\alpha;\gamma;z)=1+\frac{\alpha}{\gamma}\frac{z}{1!}+\frac{\alpha(\alpha+1)}{\gamma(\gamma+1)}\frac{z^