}{d\rho}, \\
c_{\Gamma\Theta} &= \frac{\alpha}{3 (3c_{123}+2\alpha \gamma_2)}\left[\left( 1-\frac{2c_{13}}{c_{14}}\right)(1+2\gamma_2)  - 6\epsilon_H \gamma_2\left( 1+\frac{2}{3}\gamma_3\right)   \right] + \frac{dP}{d\rho}, \\
c_{\Gamma W} &= \frac{2\gamma_1(1+2\gamma_2)}{3\left( 2\gamma_1-1 \right) }, \\
c_{\Gamma X} & = \frac{4\alpha\gamma_1}{3\left( 2\gamma_1-1 \right)(3c_{123}+2\alpha \gamma_2)}\left[\left(1+ \frac{c_{13}}{c_{14}}\right) (1+2\gamma_2)^2 + \frac{3c_{13}}{\alpha}\left(1+2\gamma_2\left[1-\epsilon_H\left(1+\frac{2}{3}\gamma_3 \right)  \right] \right) \right], \\ \label{c_GammaY}
c_{\Gamma Y} &= \frac{2\gamma_1(1+2\gamma_2)}{9\left( 1-2\gamma_1\right) }.\end{aligned}$$ Note that in and the perturbed fluid variables are those for the dark energy fluid.

In order to ensure these results are truly gauge invariant, we must do the same calculation in the synchronous gauge. However, as mentioned previously, we now have an extra variable, $T$, to deal with. Therefore, let us suppose that in the synchronous gauge we find that $$\begin{aligned}
w_{\mathrm{de}}\Pi^S &= c_{\Pi \Delta} \Delta + c_{\Pi \Theta}\hat{\Theta} + c_{\Pi X}X + c_{\Pi Y}K^2 Y+c_{\Pi T}T, \\
w_{\mathrm{de}}\Gamma &= c_{\Gamma\Delta} \Delta + c_{\Gamma\Theta}\hat{\Theta} + c_{\Gamma W} W + c_{\Gamma X} X + c_{\Gamma Y} K^
} { d\rho }, \\ 
 c_{\Gamma\Theta } & = \frac{\alpha}{3 (3c_{123}+2\alpha \gamma_2)}\left[\left (1-\frac{2c_{13}}{c_{14}}\right)(1 + 2\gamma_2)   - 6\epsilon_H \gamma_2\left (1+\frac{2}{3}\gamma_3\right)    \right ] + \frac{dP}{d\rho }, \\ 
 c_{\Gamma W } & = \frac{2\gamma_1(1 + 2\gamma_2)}{3\left (2\gamma_1 - 1 \right) }, \\ 
 c_{\Gamma X } & = \frac{4\alpha\gamma_1}{3\left (2\gamma_1 - 1 \right)(3c_{123}+2\alpha \gamma_2)}\left[\left(1 + \frac{c_{13}}{c_{14}}\right) (1 + 2\gamma_2)^2 + \frac{3c_{13}}{\alpha}\left(1 + 2\gamma_2\left[1-\epsilon_H\left(1+\frac{2}{3}\gamma_3 \right)   \right ] \right) \right ], \\ \label{c_GammaY } 
 c_{\Gamma Y } & = \frac{2\gamma_1(1 + 2\gamma_2)}{9\left (1 - 2\gamma_1\right) } .\end{aligned}$$ Note that in and the perturbed fluid variables are those for the blue department of energy fluid. 

 In order to ensure these result are sincerely gauge invariant, we must do the like calculation in the synchronous gauge. However, as mentioned previously, we now suffer an extra variable, $ T$, to distribute with. Therefore, let us speculate that in the synchronous bore we find that $ $ \begin{aligned } 
 w_{\mathrm{de}}\Pi^S & = c_{\Pi \Delta } \Delta + c_{\Pi \Theta}\hat{\Theta } + c_{\Pi X}X + c_{\Pi Y}K^2 Y+c_{\Pi T}T, \\ 
 w_{\mathrm{de}}\Gamma & = c_{\Gamma\Delta } \Delta + c_{\Gamma\Theta}\hat{\Theta } + c_{\Gamma W } W + c_{\Gamma X } X + c_{\Gamma Y } K^
}{d\rhl}, \\
c_{\Gamma\Theta} &= \frac{\alpha}{3 (3c_{123}+2\alpha \gamma_2)}\leyr[\left( 1-\hrac{2c_{13}}{c_{14}}\rjght)(1+2\gammx_2)  - 6\epsilon_H \gamma_2\left( 1+\frac{2}{3}\gemma_3\eight)   \right] + \frac{dP}{d\rho}, \\
c_{\Eamma W} &= \vrac{2\gammq_1(1+2\gamnq_2)}{3\left( 2\gamme_1-1 \right) }, \\
c_{\Gammz X} & = \hrac{4\alpha\gamma_1}{3\lgft( 2\gamma_1-1 \richt)(3c_{123}+2\alpha \gammd_2)}\ldfc[\left(1+ \frac{c_{13}}{c_{14}}\right) (1+2\gamma_2)^2 + \frac{3c_{13}}{\alpha}\lqft(1+2\gamms_2\lfft[1-\epsilon_H\lefj(1+\frac{2}{3}\bwmma_3 \ginht)  \right] \right) \right], \\ \label{c_GajmaY}
c_{\Gakma Y} &= \frac{2\gakma_1(1+2\gamma_2)}{9\left( 1-2\gamma_1\right) }.\enf{alihned}$$ Note that in wnd the perjhrbqe fluid varixbles are uhmse for thg dark energy fluid.

In order to evsure these resolfs dre truly geuge igvariant, we must do dhe samr calculation ln thx sybchronous gauge. Howevxr, as mentioned prevyously, we nkw have an extra caeiablg, $T$, tm dexo wkth. Tiersfore, pet us supposs that in tye synchronous gaugt wq find that $$\befin{alidnqd}
w_{\mathrm{de}}\Pi^S &= c_{\Pi \Delta} \Delta + c_{\Pi \Thtta}\haf{\Theta} + c_{\Pi X}X + c_{\Pi Y}K^2 T+c_{\Pi T}T, \\
w_{\mathrm{de}}\Gammw &= c_{\Gamma\Qelta} \Delta + c_{\Gamma\Theta}\hat{\Theta} + c_{\Gamma W} W + c_{\Gakma X} B + e_{\Nwona Y} K^
}{d\rho}, \\ c_{\Gamma\Theta} &= \frac{\alpha}{3 (3c_{123}+2\alpha \gamma_2)}\left[\left( 6\epsilon_H 1+\frac{2}{3}\gamma_3\right) \right] \frac{dP}{d\rho}, \\ c_{\Gamma }, c_{\Gamma X} & \frac{4\alpha\gamma_1}{3\left( 2\gamma_1-1 \right)(3c_{123}+2\alpha \frac{c_{13}}{c_{14}}\right) (1+2\gamma_2)^2 + \frac{3c_{13}}{\alpha}\left(1+2\gamma_2\left[1-\epsilon_H\left(1+\frac{2}{3}\gamma_3 \right) \right] \right], \\ \label{c_GammaY} c_{\Gamma Y} &= \frac{2\gamma_1(1+2\gamma_2)}{9\left( 1-2\gamma_1\right) }.\end{aligned}$$ Note that in and perturbed fluid variables are those for the dark energy fluid. In order to these are gauge we must do the same calculation in the synchronous gauge. However, as mentioned previously, we now an extra variable, $T$, to deal with. Therefore, us suppose that in synchronous gauge we find that w_{\mathrm{de}}\Pi^S c_{\Pi \Delta} + \Theta}\hat{\Theta} c_{\Pi X}X + Y}K^2 Y+c_{\Pi T}T, \\ w_{\mathrm{de}}\Gamma &= c_{\Gamma\Delta} \Delta + c_{\Gamma\Theta}\hat{\Theta} + c_{\Gamma W} W + c_{\Gamma X} + c_{\Gamma
}{d\rho}, \\
c_{\Gamma\Theta} &= \frac{\alpha}{3 (3c_{123}+2\Alpha \gamma_2)}\Left[\lEft( 1-\FraC{2c_{13}}{C_{14}}\rigHt)(1+2\gaMma_2)  - 6\epsilon_H \gamMA_2\lefT( 1+\frac{2}{3}\gamma_3\right)   \right] + \frAc{dP}{d\RhO}, \\
C_{\GamMA W} &= \Frac{2\gAmma_1(1+2\gamMA_2)}{3\lEFT( 2\gaMmA_1-1 \rIghT) }, \\
c_{\gAmMa X} & = \frAc{4\aLpha\gamMa_1}{3\left( 2\gammA_1-1 \riGhT)(3c_{123}+2\alpha \gamma_2)}\LEfT[\left(1+ \frac{c_{13}}{C_{14}}\riGht) (1+2\gamma_2)^2 + \frac{3C_{13}}{\alPha}\lefT(1+2\gAmmA_2\Left[1-\ePsiLon_H\lEft(1+\fraC{2}{3}\Gamma_3 \rIght)  \right] \RiGHt) \righT], \\ \Label{c_GAMMay}
c_{\GaMma Y} &= \frac{2\gamma_1(1+2\gammA_2)}{9\LeFT( 1-2\gamma_1\right) }.\end{AligneD}$$ NOTe THAt iN anD the perturBeD fluiD VariablES aRE THosE For the dark eneRgy fluid.

In oRDer To ensuRe TheSE resulTs are TrULy gAuge invariaNt, we Must do the Same caLCulatioN In the syNchronOus GauGe. HoWEvEr, As mEnTIonED pRevIOusLy, we now hAvE aN extrA varIABLE, $T$, to DeaL witH. TherEfore, let us supPosE thaT In tHe synChronOus gAuGe we fInd thaT $$\begiN{aLigned}
w_{\mathrm{de}}\pi^S &= c_{\pi \Delta} \DeLta + C_{\PI \ThEtA}\hat{\THEta} + c_{\Pi x}X + c_{\pi Y}k^2 Y+c_{\Pi T}T, \\
W_{\mathrm{DE}}\GaMmA &= C_{\gAmMa\Delta} \Delta + c_{\Gamma\thETA}\hAt{\Theta} + c_{\gamma W} w + C_{\GAmMA X} X + c_{\GammA Y} k^
}{d\rho}, \\
c_{\Gamma\The ta} &= \fr ac{\a lph a}{ 3(3c_ {123 }+2\alpha \gam m a_2) }\left[\left( 1-\frac{ 2c_{1 3} } {c_{ 1 4} }\rig ht)(1+2 \ ga m m a_2 ) -  6\ ep s il on_H\ga mma_2\l eft( 1+\fr ac{ 2} {3}\gamma_3\ r ig ht)   \rig ht]  + \frac{dP} {d\ rho},\\ 
c_ { \Gamm a W } &=\frac{ 2 \gamma _1(1+2\ga mm a _2)}{3 \ left( 2 \ g am ma_1 -1 \right) }, \\c _{ \ Gamma X} & = \ frac{4 \a l ph a \ gam ma_ 1}{3\left(  2 \gamm a _1-1 \r i gh t ) ( 3c_ { 123}+2\alpha\gamma_2)}\ l eft [\left (1 + \ f rac{c_ {13}} {c _ {14 }}\right) ( 1+2\ gamma_2)^ 2 + \f r ac{3c_{ 1 3}}{\al pha}\l eft (1+ 2\ga m ma _2 \le ft [ 1-\ e ps ilo n _H\ left(1+\ fr ac {2}{3 }\ga m m a _ 3 \r igh t) \righ t] \right) \r igh t],\ \ \ label {c_Ga mmaY }c_{\G amma Y } &=\f rac{2\gamma_1(1 +2\g amma_2)}{ 9\l ef t(1- 2\gam m a_1\ri ght ) } .\end{a ligned} $ $ N ot e t ha t in and the pertu rb e d  f luid var iables ar et hose for  t hedark e nergy  flu i d. 

In ord er toe ns ur e these  r esults  a retru ly ga u ge i nvaria nt, we m ust d o  the same calc u lation in the sy n c hr o nous  ga uge. Howeve r, a s  men tion e dpre v iousl y, we  n o wh ave an extra variab le , $T$,  to d eal with. The refore, le t u s suppos e th a ti n the synchron ous g auge we fi n d that $ $\beg in{align ed}
w_{\m a t hrm{de}} \Pi ^S&=c_{ \ P i\Delta} \Delt a + c_ {\ Pi \The ta} \hat{\T het a}+ c _{\ Pi  X}X + c_ {\Pi Y}K ^2  Y +c _{ \Pi  T}T, \\
w_{\m at hrm {d e}} \Gamm a  &= c_ {\Gam ma\D el ta }  \D elta +c _{ \ G amma \T he ta}\ hat {\ Theta } +c _{\ Gamma W } W + c_{ \Ga m ma X }X+ c_{\G amma Y} K^
}{d\rho}, \\
c_{\Gamma\Theta}_&= \frac{\alpha}{3_(3c_{123}+2\alpha \gamma_2)}\left[\left( 1-\frac{2c_{13}}{c_{14}}\right)(1+2\gamma_2) _- 6\epsilon_H_\gamma_2\left(_1+\frac{2}{3}\gamma_3\right) __\right] + \frac{dP}{d\rho},_\\
c_{\Gamma W} &=_\frac{2\gamma_1(1+2\gamma_2)}{3\left( 2\gamma_1-1 \right) },_\\
c_{\Gamma X} &_=_\frac{4\alpha\gamma_1}{3\left( 2\gamma_1-1 \right)(3c_{123}+2\alpha \gamma_2)}\left[\left(1+ \frac{c_{13}}{c_{14}}\right) (1+2\gamma_2)^2 + \frac{3c_{13}}{\alpha}\left(1+2\gamma_2\left[1-\epsilon_H\left(1+\frac{2}{3}\gamma_3 \right)  \right] \right) \right], \\_\label{c_GammaY}
c_{\Gamma_Y} &=_\frac{2\gamma_1(1+2\gamma_2)}{9\left(_1-2\gamma_1\right)_}.\end{aligned}$$ Note that in and_the perturbed fluid variables are_those for_the dark energy fluid.

In order to ensure these_results_are truly gauge_invariant, we must do the same calculation in the_synchronous gauge. However, as mentioned previously,_we now have_an_extra_variable, $T$, to deal_with. Therefore, let us suppose that_in the synchronous gauge we find_that $$\begin{aligned}
w_{\mathrm{de}}\Pi^S &= c_{\Pi \Delta} \Delta +_c_{\Pi \Theta}\hat{\Theta} + c_{\Pi X}X +_c_{\Pi Y}K^2 Y+c_{\Pi T}T, \\
w_{\mathrm{de}}\Gamma_&= c_{\Gamma\Delta}_\Delta + c_{\Gamma\Theta}\hat{\Theta} + c_{\Gamma_W} W +_c_{\Gamma X}_X + c_{\Gamma_Y} K^