text:
# Cosmological Perturbations with Multiple Fluids and Fields

## 1 Introduction

Cosmological models with multiple components of underlying energy-momentum content are not only general but also necessarily appear in many cosmological situations including the early and the late evolution stages of our world models. Assuming the gauge degrees of freedom are properly fixed, for the $`n`$ component medium with fluids or fields in Einstein gravity we anticipate a set of $`n`$ coupled second-order differential equations for the scalar-type perturbation, a set of $`n`$ (in general) coupled first-order differential equations for the vector-type perturbation, and one second-order differential equation for the tensor-type perturbation. The situations of the vector- and tensor-type perturbations are rather trivial to handle, and, in fact, the situation is similar even in some classes of generalized versions of gravity theories . However, the scalar-type perturbation in such a multi-component system is naturally more complicated, and often requires numerical methods. Not only the equations are complicated for the scalar-type perturbation, the behavior and the way of handling the equations depend on our choice of the gauge condition which we need to make. In either analytic or numerical handling of the scalar-type perturbations, it is useful and is often necessary to examine the problems from the perspectives of different gauge conditions.

In the following, for the benefit of future studies, we would like to present various useful forms of the scalar-type perturbation equations considering arbitrary numbers of imperfect fluids and scalar fields with general mutual interactions among the components in Einstein gravity. Equations are presented in the gauge-ready forms which allow to take any of our favorite gauge conditions easily, and also allow easy translation of a solution derived in one gauge condition to the ones in other gauge conditions. We also present closed form equations in terms of the curvature and isocurvature perturbation variables in several useful forms, and derive general solutions if available. We recently have presented similar work on systems with multiple minimally coupled scalar fields . In the present work we present a complementary part now including the multiple fluids as well as the interactions among fluids and fields; we also extend  by considering general $`K`$ and $`\mathrm{\Lambda }`$ and the fluid formulation of the field perturbations in §4.

The manuscript is mainly concerned with the formal aspects of perturbations in the multi-component situations. There are some advantages of investigating perturbations in a formal way keeping their most general forms. For example, we emphasize that the basic set of equations for a single component case remains valid even in the multi-component situations with additional equations for the components. Also, by including most general energy-momentum tensor in the fluids system we will show that the same equations remain valid for the mixture of fluids and fields, and sometimes even for the generalized gravity theories. By comparing the fluid system and field system we can clarify the similarity and difference between them. Although, we intentionally have kept the presentation formal, many new results are presented together with the known ones; see §6 for a summary. Previous studies of perturbations in the multi-component fluid system in a similar spirit can be found in . We start by introducing our notations, fundamental equations, and strategy in §2. We set $`c1`$.

## 2 Perturbed world model and basic equations

We consider an arbitrary number of mutually interacting imperfect fluids in Einstein gravity. The gravitational field equation and the energy-momentum conservation equation are

$`G_{ab}8\pi GT_{ab}\mathrm{\Lambda }g_{ab},T_{a;b}^b=0.`$ (1)
The energy-momentum tensor can be decomposed into the sum of the individual one as

$`T_{ab}={\displaystyle \underset{l}{}}T_{(l)ab},`$ (2)
and the energy-momentum conservation gives

$`T_{(i)a;b}^bQ_{(i)a},{\displaystyle \underset{l}{}}Q_{(l)a}=0,`$ (3)
where $`(i)`$ indicates the $`i`$-th component of $`n`$ matters (including fluids and fields) with $`i,j,k,\mathrm{}=1,2,\mathrm{},n`$.

We consider the general perturbations in the FLRW world model. In the spatially homogeneous and isotropic background we are considering, to the linear order the three different types of perturbation decouple from each other and evolve independently. The presence of tensor-type anisotropic stress contributes as the source or sink of the graviational wave of the perturbed spacetime. The rotational perturbation of the individual fluid is described by the angular-momentum conservation where the vector-type anisotropic stress and the mutual interaction terms among fluids can work as the source or sink of the rotation of individual component. The most general situations of vector- and tensor-type perturbations are presented in . In the following we will consider the scalar-type perturbation only.

A metric with the most general scalar-type perturbation in the FLRW background is

$`ds^2=a^2\left(1+2\alpha \right)d\eta ^22a^2\beta _{,\alpha }d\eta dx^\alpha +a^2\left[g_{\alpha \beta }^{(3)}\left(1+2\phi \right)+2\gamma _{,\alpha |\beta }\right]dx^\alpha dx^\beta ,`$ (4)
where $`a(\eta )`$ is the cosmic scale factor. $`\alpha `$, $`\beta `$, $`\gamma `$ and $`\phi `$ are spacetime dependent perturbed order variables. A vertical bar $`|`$ indicates a covariant derivative based on $`g_{\alpha \beta }^{(3)}`$. Considering the FLRW background and the scalar-type perturbation, the energy-momentum tensor can be decomposed into the collective fluid quantities as

$`T_0^0=\left(\overline{\mu }+\delta \mu \right),T_\alpha ^0={\displaystyle \frac{1}{k}}\left(\mu +p\right)v_{,\alpha },T_\beta ^\alpha =\left(\overline{p}+\delta p\right)\delta _\beta ^\alpha +\left({\displaystyle \frac{1}{k^2}}^\alpha _\beta +{\displaystyle \frac{1}{3}}\delta _\beta ^\alpha \right)\pi ^{(s)},`$ (5)
where $`\mu (\overline{\mu }+\delta \mu )`$, $`p(\overline{p}+\delta p)`$, $`v`$, and $`\pi ^{(s)}`$ are the energy density, the isotropic pressure, the velocity or flux related variable, and the anisotropic stress, respectively; $`k`$ is the comoving wavenumber. In terms of the individual matter’s fluid quantities we have

$`\overline{\mu }={\displaystyle \underset{l}{}}\overline{\mu }_{(l)},\delta \mu ={\displaystyle \underset{l}{}}\delta \mu _{(l)},`$ (6)
and similarly for $`\overline{p}`$, $`\delta p`$, $`(\mu +p)v`$, and $`\pi ^{(s)}`$; an overbar indicates the background order quantity, and will be ignored unless necessary. The interaction terms among fluids introduced in eq. (3) are decomposed as

$`Q_{(i)0}a\left[\overline{Q}_{(i)}(1+\alpha )+\delta Q_{(i)}\right],Q_{(i)\alpha }J_{(i),\alpha }.`$ (7)
The equations for the background are:

$`H^2={\displaystyle \frac{8\pi G}{3}}\mu +{\displaystyle \frac{\mathrm{\Lambda }}{3}}{\displaystyle \frac{K}{a^2}},\dot{H}=4\pi G(\mu +p)+{\displaystyle \frac{K}{a^2}},`$ (8)
$`\dot{\mu }_{(i)}+3H\left(\mu _{(i)}+p_{(i)}\right)=Q_{(i)},`$ (9)
where an overdot indicates the time derivative based on $`t`$ defined as $`dtad\eta `$. $`K`$ is the normalized spatial curvature, and $`H\dot{a}/a`$. The two equations in eq. (8) follow from $`G_0^0`$ and $`G_\alpha ^\alpha 3G_0^0`$ components of the field equation in eq. (1), respectively, and eq. (9) follows from eq. (3). By adding eq. (9) over components we have $`\dot{\mu }+3H(\mu +p)=0`$ which also follows from eq. (8).

We present a complete set of equations describing the scalar-type perturbation without fixing the temporal gauge condition, i.e., in a gauge-ready form :

$`\kappa 3\left(\dot{\phi }+H\alpha \right)+{\displaystyle \frac{k^2}{a^2}}\chi ,`$ (10)
$`{\displaystyle \frac{k^23K}{a^2}}\phi +H\kappa =4\pi G\delta \mu ,`$ (11)
$`\kappa {\displaystyle \frac{k^23K}{a^2}}\chi =12\pi G{\displaystyle \frac{a}{k}}(\mu +p)v,`$ (12)
$`\dot{\chi }+H\chi \alpha \phi =8\pi G{\displaystyle \frac{a^2}{k^2}}\pi ^{(s)},`$ (13)
$`\dot{\kappa }+2H\kappa +\left(3\dot{H}{\displaystyle \frac{k^2}{a^2}}\right)\alpha =4\pi G\left(\delta \mu +3\delta p\right),`$ (14)
$`\delta \dot{\mu }_{(i)}+3H\left(\delta \mu _{(i)}+\delta p_{(i)}\right)={\displaystyle \frac{k}{a}}\left(\mu _{(i)}+p_{(i)}\right)v_{(i)}+\dot{\mu }_{(i)}\alpha +\left(\mu _{(i)}+p_{(i)}\right)\kappa +\delta Q_{(i)},`$ (15)
$`{\displaystyle \frac{1}{a^4(\mu _{(i)}+p_{(i)})}}\left[a^4(\mu _{(i)}+p_{(i)})v_{(i)}\right]^{}={\displaystyle \frac{k}{a}}\left[\alpha +{\displaystyle \frac{1}{\mu _{(i)}+p_{(i)}}}\left(\delta p_{(i)}{\displaystyle \frac{2}{3}}{\displaystyle \frac{k^23K}{k^2}}\pi _{(i)}^{(s)}J_{(i)}\right)\right],`$ (16)
where $`\chi a(\beta +a\dot{\gamma })`$. $`\chi `$, $`\phi `$ and $`\kappa `$ correspond to the shear, the three-space curvature and the perturbed expansion of the normal-frame vector field, respectively, see  for details. $`_\alpha `$ and $`\mathrm{\Delta }`$ are the covariant derivative and the Laplacian based on $`g_{\alpha \beta }^{(3)}`$. Equations (10-14) are: the definition of $`\kappa `$, ADM energy constraint ($`G_0^0`$ component of the field equation), momentum constraint ($`G_\alpha ^0`$ component), tracefree part of ADM propagation ($`G_\beta ^\alpha \frac{1}{3}\delta _\beta ^\alpha G_\gamma ^\gamma `$ component), and Raychaudhuri equation ($`G_\gamma ^\gamma G_0^0`$ component), respectively. Equations (15,16) follow from the energy conservation ($`T_{(i)0;b}^b=Q_{(i)0}`$) and the momentum conservation ($`T_{(i)\alpha ;b}^b=Q_{(i)\alpha }`$), respectively. We note that eqs. (15,16) in the same forms are valid even in a class of generalized gravity theories considered in ; see eqs. (46,47) in . Equations (10-14) also remain valid in the generalized gravity where we can reinterprete the fluid quantities as the effective ones; see the Appendix B of .

By adding properly eqs. (15,16) over all components of the fluids, and using properties in eq. (6), we get equations for the collective fluid quantities as:

$`\delta \dot{\mu }+3H\left(\delta \mu +\delta p\right)=(\mu +p)\left(\kappa 3H\alpha {\displaystyle \frac{k}{a}}v\right),`$ (17)
$`{\displaystyle \frac{\left[a^4(\mu +p)v\right]^{}}{a^4(\mu +p)}}={\displaystyle \frac{k}{a}}\left[\alpha +{\displaystyle \frac{1}{\mu +p}}\left(\delta p{\displaystyle \frac{2}{3}}{\displaystyle \frac{k^23K}{k^2}}\pi ^{(s)}\right)\right].`$ (18)
We introduce

$`e_{(i)}\delta p_{(i)}c_{(i)}^2\delta \mu _{(i)},c_{(i)}^2\dot{p}_{(i)}/\dot{\mu }_{(i)},e\delta pc_s^2\delta \mu ,c_s^2\dot{p}/\dot{\mu },`$
$`w_{(i)}p_{(i)}/\mu _{(i)},wp/\mu ,`$ (19)
where $`e_{(i)}`$ and $`e`$ are the entropic perturbations. Equation (15) can be written in another form as

$`\dot{\delta }_{(i)}+3H\left(c_{(i)}^2w_{(i)}\right)\delta _{(i)}+3H\left(1+w_{(i)}\right)q_{(i)}\delta _{(i)}`$
$`=\left(1+w_{(i)}\right)\left[\kappa 3H\left(1q_{(i)}\right)\alpha {\displaystyle \frac{k}{a}}v_{(i)}\right]+{\displaystyle \frac{1}{\mu _{(i)}}}\left(3He_{(i)}+\delta Q_{(i)}\right),`$ (20)
where $`Q_{(i)}3H(\mu _{(i)}+p_{(i)})q_{(i)}`$, and we used $`\dot{w}_{(i)}=3H(c_{(i)}^2w_{(i)})(1+w_{(i)})(1q_{(i)})`$; $`\delta _{(i)}\delta \mu _{(i)}/\mu _{(i)}`$.

In the following we briefly summarize the gauge-ready strategy suggested by Bardeen  and elaborated in ; see  for a recent extension. Under the gauge transformation $`\stackrel{~}{x}^a=x^a+\xi ^a`$ with $`\xi ^ta\xi ^0`$ ($`0=\eta `$) the perturbed metric and fluid quantities change as :

$`\stackrel{~}{\alpha }=\alpha \dot{\xi }^t,\stackrel{~}{\phi }=\phi H\xi ^t,\stackrel{~}{\kappa }=\kappa +\left(3\dot{H}+{\displaystyle \frac{\mathrm{\Delta }}{a^2}}\right)\xi ^t,\stackrel{~}{\chi }=\chi \xi ^t,\stackrel{~}{v}=v{\displaystyle \frac{k}{a}}\xi ^t,\stackrel{~}{v}_{(i)}=v_{(i)}{\displaystyle \frac{k}{a}}\xi ^t,`$
$`\delta \stackrel{~}{\mu }=\delta \mu \dot{\mu }\xi ^t,\delta \stackrel{~}{\mu }_{(i)}=\delta \mu _{(i)}\dot{\mu }_{(i)}\xi ^t,\delta \stackrel{~}{p}=\delta p\dot{p}\xi ^t,\delta \stackrel{~}{p}_{(i)}=\delta p_{(i)}\dot{p}_{(i)}\xi ^t,`$
$`\delta \stackrel{~}{Q}_{(i)}=\delta Q_{(i)}\dot{Q}_{(i)}\xi ^t,\stackrel{~}{J}_{(i)}=J_{(i)}+Q_{(i)}\xi ^t.`$ (21)
As the temporal gauge fixing condition we can impose one condition in any of these temporally gauge dependent variables: $`\alpha 0`$ (synchronous gauge), $`\phi 0`$ (uniform-curvature gauge), $`\kappa 0`$ (uniform-expansion gauge), $`\chi 0`$ (zero-shear gauge), $`v/k0`$ (comoving gauge), $`\delta \mu 0`$ (uniform-density gauge), $`\delta p0`$, $`v_{(i)}/k0`$, $`\delta \mu _{(i)}0`$, $`\delta p_{(i)}0`$, etc. By examining these we notice that, except for the synchronous gauge condition (which fixes $`\alpha =0`$), each of the gauge conditions fixes the temporal gauge mode completely. Thus, a variable in such a gauge condition uniquely corresponds to a gauge-invariant combination which combines the variable concerned and the variable used in the gauge condition. As examples, one can recognize the following combinations are gauge-invariant.

$`\delta \mu _v\delta \mu {\displaystyle \frac{a}{k}}\dot{\mu }v,\phi _v\phi {\displaystyle \frac{aH}{k}}v,\phi _\chi \phi H\chi H\chi _\phi ,v_\chi v{\displaystyle \frac{k}{a}}\chi ,`$ (22)
etc. The original study of cosmological perturbation by Lifshitz  was made in the synchronous gauge. The zero-shear gauge and the comoving gauge were introduced by Harrison  and Nariai , respectively. The combination $`\phi _v`$ was first introduced by Lukash  and Bardeen . The gauge-invariant combinations $`\phi _\chi `$ and $`v_\chi `$ are equivalent to $`\phi `$ and $`v`$ in the zero-shear gauge which takes $`\chi 0`$ as the gauge condition, and $`\delta \mu _v`$ is the same as $`\delta \mu `$ in the comoving gauge, etc. In this way, we can systematically construct various gauge-invariant combinations for a given variable. Since there exist many gauge-invariant combinations even for a given variable, say for $`\delta \mu `$ we could have $`\delta \mu _\phi `$, $`\delta \mu _v`$, $`\delta \mu _{v_{(i)}}`$, etc., using our notation we can easily recognize which gauge-invariant combination we are using. Generally, we do not know the suitable gauge condition a priori. The proposal made in  is that we write the set of equation without fixing the (temporal) gauge condition and arrange the equations so that we can implement easily various fundamental gauge conditions: eqs. (10-16) are arranged accordingly. We term it a gauge-ready approach.

## 3 Fluid system

### 3.1 Newtonian analog

The following equations most closely resemble Newtonian hydrodynamic equations in the context of linear perturbation of the Friedmann world model: from eqs. (11,12), eqs. (17,18,12), eqs. (13,18), eqs. (10,12,13), and eq (13), respectively we have

$`{\displaystyle \frac{k^23K}{a^2}}\phi _\chi =4\pi G\delta \mu _v,`$ (23)
$`\delta \dot{\mu }_v+3H\delta \mu _v={\displaystyle \frac{k^23K}{k^2}}\left[{\displaystyle \frac{k}{a}}(\mu +p)v_\chi +2H\pi ^{(s)}\right],`$ (24)
$`\dot{v}_\chi +Hv_\chi ={\displaystyle \frac{k}{a}}\left(\phi _\chi +{\displaystyle \frac{\delta p_v}{\mu +p}}8\pi G{\displaystyle \frac{a^2}{k^2}}\pi ^{(s)}{\displaystyle \frac{2}{3}}{\displaystyle \frac{k^23K}{k^2}}{\displaystyle \frac{\pi ^{(s)}}{\mu +p}}\right),`$ (25)
$`\dot{\phi }_\chi +H\phi _\chi =4\pi G{\displaystyle \frac{a}{k}}\left(\mu +p\right)v_\chi 8\pi GH{\displaystyle \frac{a^2}{k^2}}\pi ^{(s)},`$ (26)
$`\alpha _\chi =\phi _\chi 8\pi G{\displaystyle \frac{a^2}{k^2}}\pi ^{(s)}.`$ (27)
Equations (23-25) can be compared with the Poisson, the continuity, and the Euler (Navier-Stokes) equations in Newtonian theory, respectively. By setting $`p=0=\pi ^{(s)}`$ and ignoring $`K=0`$ we recover the well known equations in the Newtonian perturbation theory with $`\phi _\chi `$, $`\delta \mu _v`$ and $`v_\chi `$ corresponding to the Newtonian potential, density and velocity perturbations, respectively. Notice the appearence of pressure terms in rather unexpected places in eqs. (24,25) and none appearing in eq. (23). Although the zero-shear gauge and the comoving gauge were first studied by Harrison and Nariai, respectively , these equations using mixed gauge-invariant combinations were first presented by Bardeen in eqs. (4.3-4.8) of . Newtonian correspondence of these fully general relativistic equations and variables was investigated in . In Fig. 1 we present the behaviors of $`\phi _\chi `$, $`v_\chi `$ ($`v_{(c)\chi }`$) and $`\delta _v`$ ($`\delta _{(c)v}`$) in the case of three component system with the photon, baryon and the cold dark matter ($`c`$). We also compared the behavior with $`\phi `$, $`v`$ ($`v_{(c)}`$) and $`\delta `$ ($`\delta _{(c)}`$) in the uniform-expansion gauge.

In  we noted that in the sub-horizon scale in matter dominated era $`\phi `$, $`\delta `$, and $`v`$ in the uniform-expansion gauge all behave like $`\phi _\chi `$, $`\delta _v`$, and $`v_\chi `$ which most closely resemble the corresponding Newtonian behavior, see also §84 of . We compared the behaviors in the uniform-expansion gauge with the variables with Newtonian behavior in Fig. 1.

From eqs. (14,17,18) we can derive

$`{\displaystyle \frac{1+w}{a^2H}}\left[{\displaystyle \frac{H^2}{a(\mu +p)}}\left({\displaystyle \frac{a^3\mu }{H}}\delta _v\right)^{}\right]^{}+c_s^2{\displaystyle \frac{k^2}{a^2}}\delta _v={\displaystyle \frac{k^23K}{a^2}}\left[{\displaystyle \frac{e}{\mu }}+2{\displaystyle \frac{1+w}{H}}\left({\displaystyle \frac{a^2H^2}{k^2}}{\displaystyle \frac{\pi ^{(s)}}{\mu +p}}\right)^{}{\displaystyle \frac{2}{3}}{\displaystyle \frac{\pi ^{(s)}}{\mu }}\right].`$ (28)
The multi-component version of this equation will be presented later, see eq. (61). This equation is fully relativistic but most closely resembes the well known Newtonian density perturbation equation; only for pressureless fluid our comoving gauge coincides with the synchronous gauge. $`\delta _v`$ is directly related to the density gradient variable introduced in the covariant approach ; see §IV and the note added in . By using eq. (23), the equation for $`\phi _\chi `$ can be derived simply from eq. (28), see eq. (35).

In the multi-component situation, although eqs. (24,25) are still valid for the total fluid quantities, we can also derive the equations for individual component. From eqs. (15,16,12) and eqs. (16,13), respectively, we have

$`\delta \dot{\mu }_{(i)v_{(i)}}+3H\delta \mu _{(i)v_{(i)}}={\displaystyle \frac{k^23K}{k^2}}\left[{\displaystyle \frac{k}{a}}(\mu _{(i)}+p_{(i)})v_{(i)\chi }+2H\pi _{(i)}^{(s)}\right]`$
$`+12\pi G{\displaystyle \frac{a}{k}}(\mu +p)(\mu _{(i)}+p_{(i)})\left(v_\chi v_{(i)\chi }\right)+\delta Q_{(i)v_{(i)}}+Q_{(i)}\alpha _{v_{(i)}}3HJ_{(i)v_{(i)}},`$ (29)
$`\dot{v}_{(i)\chi }+Hv_{(i)\chi }+{\displaystyle \frac{Q_{(i)}}{\mu _{(i)}+p_{(i)}}}v_{(i)\chi }`$
$`={\displaystyle \frac{k}{a}}\left(\phi _\chi +{\displaystyle \frac{\delta p_{(i)v_{(i)}}}{\mu _{(i)}+p_{(i)}}}8\pi G{\displaystyle \frac{a^2}{k^2}}\pi ^{(s)}{\displaystyle \frac{2}{3}}{\displaystyle \frac{k^23K}{k^2}}{\displaystyle \frac{\pi _{(i)}^{(s)}}{\mu _{(i)}+p_{(i)}}}{\displaystyle \frac{J_{(i)\chi }}{\mu _{(i)}+p_{(i)}}}\right).`$ (30)
By adding eqs. (29,30) over the components properly we have eqs. (24,25). Later we will notice that the above equations are not only valid for multiple fluids but also for the system with additional multiple fields.

### 3.2 Curvature perturbations

We introduce

$`\mathrm{\Phi }\phi _v{\displaystyle \frac{K/a^2}{4\pi G(\mu +p)}}\phi _\chi ,`$ (31)
which becomes $`\phi _v`$, the Lukash variable, for vanishing $`K`$. Using eqs. (22,26) we can show

$`\mathrm{\Phi }={\displaystyle \frac{H^2}{4\pi G(\mu +p)a}}\left({\displaystyle \frac{a}{H}}\phi _\chi \right)^{}+{\displaystyle \frac{2H^2}{\mu +p}}{\displaystyle \frac{a^2}{k^2}}\pi ^{(s)}.`$ (32)
On the other hand, from eqs. (28,23,32) or eqs. (23,26,25) we have

$`\dot{\mathrm{\Phi }}={\displaystyle \frac{Hc_s^2}{4\pi G(\mu +p)}}{\displaystyle \frac{k^2}{a^2}}\phi _\chi {\displaystyle \frac{H}{\mu +p}}\left(e{\displaystyle \frac{2}{3}}\pi ^{(s)}\right).`$ (33)
Thus, in a pressureless case, if we could ignore $`e`$ and $`\pi ^{(s)}`$, $`\mathrm{\Phi }(𝐱,t)=C(𝐱)`$ is an exact solution valid in general scale . Combining eqs. (32,33) we have closed form equations for $`\mathrm{\Phi }`$ and $`\phi _\chi `$:

$`{\displaystyle \frac{H^2c_s^2}{(\mu +p)a^3}}\left[{\displaystyle \frac{(\mu +p)a^3}{H^2c_s^2}}\dot{\mathrm{\Phi }}\right]^{}+c_s^2{\displaystyle \frac{k^2}{a^2}}\mathrm{\Phi }={\displaystyle \frac{H^2c_s^2}{(\mu +p)a^3}}\left[{\displaystyle \frac{a^3}{c_s^2H}}\left(e{\displaystyle \frac{2}{3}}\pi ^{(s)}\right)\right]^{}+2H^2{\displaystyle \frac{c_s^2}{\mu +p}}\pi ^{(s)},`$ (34)
$`{\displaystyle \frac{\mu +p}{H}}\left[{\displaystyle \frac{H^2}{a(\mu +p)}}\left({\displaystyle \frac{a}{H}}\phi _\chi \right)^{}\right]^{}+c_s^2{\displaystyle \frac{k^2}{a^2}}\phi _\chi =4\pi G\left(e{\displaystyle \frac{2}{3}}\pi ^{(s)}\right){\displaystyle \frac{4\pi G(\mu +p)}{H}}\left({\displaystyle \frac{2H^2}{\mu +p}}{\displaystyle \frac{a^2}{k^2}}\pi ^{(s)}\right)^{}.`$ (35)
Equation (34) is valid for $`c_s^20`$; for $`c_s^2=0`$ we have eq. (33) instead.

In the case of multiple ideal fluids without direct interactions among them, later we will see that the differences in the sound velocities cause nonvanishing $`e`$ term which mixes the curvature perturbation $`\mathrm{\Phi }`$ with the isocurvature modes, see eq. (46). In §4 we will see that a minimally coupled scalar field has a nonvanishing entropic perturbation, see eq. (76). Thus, we can derive a similar equation for $`\mathrm{\Phi }`$ as in eq. (84) which is the form valid even in the multi-component case.

Although introduced in a different way, a variable proportional to $`\mathrm{\Phi }`$ was first introduced by Field and Shepley in 1968 , see also . In the single component ideal fluid case, using $`vz\mathrm{\Phi }`$ with $`za\sqrt{\mu +p}/(c_sH)`$ and $`\frac{}{\eta }`$, eq. (34) can be written as

$`v^{\prime \prime }+\left(c_s^2k^2z^{\prime \prime }/z\right)v=0,`$ (36)
which is a well known equation in the cosmological perturbation, . This equation can be found in eq. (45) of  by Field and Shepley. The variable $`H`$ in eq. (43) of  and $`\varphi `$ in eq. (3.3) of  are proportional to $`v`$, and $`\varphi `$ in eq. (122) of  is the same as $`\mathrm{\Phi }`$<sup>1</sup><sup>1</sup>1 Ignorant of this rich history, and even forgetting our own contribution unfortunately, we have rediscovered this important variable $`\mathrm{\Phi }`$ in . .

The $`\zeta `$ variable introduced in  is the same as $`\phi _\delta \phi +\frac{\delta \mu }{3(\mu +p)}`$. From eqs. (10,17) we can derive

$`\dot{\phi }_\delta ={\displaystyle \frac{k}{3a}}v_\chi {\displaystyle \frac{He}{\mu +p}}.`$ (37)
Using eqs. (13,18) we have

$`\ddot{\phi }_\delta +(23c_s^2)H\dot{\phi }_\delta +c_s^2{\displaystyle \frac{k^2}{a^2}}\phi _\delta =\left({\displaystyle \frac{1}{3}}+c_s^2\right){\displaystyle \frac{k^2}{a^2}}\phi _\chi `$
$`\left({\displaystyle \frac{He}{\mu +p}}\right)^{}\left[(23c_s^2)H^2+{\displaystyle \frac{k^2}{3a^2}}\right]{\displaystyle \frac{e}{\mu +p}}+{\displaystyle \frac{8\pi G}{3}}\pi ^{(s)}+{\displaystyle \frac{2}{9}}{\displaystyle \frac{k^23K}{a^2}}{\displaystyle \frac{\pi ^{(s)}}{\mu +p}}.`$ (38)
We can also derive the closed form equation for $`\phi _\delta `$. From the definition of $`\phi _\delta `$, evaluating it in the comoving gauge, and using eqs. (23,31) we have

$`\phi _\delta (k,t)=\mathrm{\Phi }+{\displaystyle \frac{k^2}{a^2}}{\displaystyle \frac{1}{12\pi G(\mu +p)}}\phi _\chi ,`$ (39)
where $`\mathrm{\Phi }`$ and $`\phi _\chi `$ are given above. For $`\phi _\kappa `$ from the definition of $`\phi _\kappa `$ and using eq. (11) and eq. (12) we have

$`\phi _\kappa =\left(1+{\displaystyle \frac{k^23K}{12\pi G(\mu +p)a^2}}\right)^1\phi _\delta =\phi _\chi \left(1+{\displaystyle \frac{k^23K}{12\pi G(\mu +p)a^2}}\right)^1{\displaystyle \frac{aH}{k}}v_\chi .`$ (40)
Using the first expression in eq. (38) we can guess the equation for $`\phi _\kappa `$. All the equations above are generally valid in the multi-component situation including fluids and fields.

In the single component ideal fluid (later we will see that it also applies for multiple ideal fluids under an adiabatic condition, thus vanishing $`e`$) with ignorable $`c_s^2k^2`$ term compared with $`z^{\prime \prime }/z`$ term in eq. (36), thus valid in the super-sound-horizon scale, we have general solutions

$`\mathrm{\Phi }(k,t)=C(k)d(k){\displaystyle \frac{k^2}{4\pi G}}{\displaystyle ^t}{\displaystyle \frac{c_s^2H^2}{a^3(\mu +p)}}𝑑t,`$ (41)
$`\phi _\chi (k,t)=4\pi GC(k){\displaystyle \frac{H}{a}}{\displaystyle ^t}{\displaystyle \frac{a(\mu +p)}{H^2}}𝑑t+d(k){\displaystyle \frac{H}{a}},`$ (42)
where $`C(k)`$ and $`d(k)`$ are the coefficients of the growing and decaying solutions, respectively. The coefficients are matched using eqs. (32,33). In order to derive $`\mathrm{\Phi }`$ from $`\phi _\chi `$ and vice versa to proper order, we need to consider the next leading order terms in the growing solution of $`\mathrm{\Phi }`$ and the decaying solution of $`\phi _\chi `$; these follow trivially from eqs. (34,35) as

$`\mathrm{\Phi }=C(k)\left\{1+k^2\left[{\displaystyle ^t}{\displaystyle \frac{a(\mu +p)}{H^2}}\left({\displaystyle ^t}{\displaystyle \frac{c_s^2H^2}{a^3(\mu +p)}}𝑑t\right)𝑑t{\displaystyle ^t}{\displaystyle \frac{a(\mu +p)}{H^2}}𝑑t{\displaystyle ^t}{\displaystyle \frac{c_s^2H^2}{a^3(\mu +p)}}𝑑t\right]\right\}`$
$`d(k){\displaystyle \frac{k^2}{4\pi G}}{\displaystyle ^t}{\displaystyle \frac{c_s^2H^2}{a^3(\mu +p)}}𝑑t,`$
$`\phi _\chi =4\pi GC(k){\displaystyle \frac{H}{a}}{\displaystyle ^t}{\displaystyle \frac{a(\mu +p)}{H^2}}𝑑t`$
$`+d(k){\displaystyle \frac{H}{a}}\left\{1+k^2\left[{\displaystyle ^t}\left({\displaystyle ^t}{\displaystyle \frac{a(\mu +p)}{H^2}}𝑑t\right){\displaystyle \frac{c_s^2H^2}{a^3(\mu +p)}}𝑑t{\displaystyle ^t}{\displaystyle \frac{a(\mu +p)}{H^2}}𝑑t{\displaystyle ^t}{\displaystyle \frac{c_s^2H^2}{a^3(\mu +p)}}𝑑t\right]\right\}.`$ (43)
Equation (41) includes $`c_s^2=0`$ limit; thus for a pressureless fluid we simply have $`\mathrm{\Phi }=C`$, . Solution for $`\phi _\delta `$ follows from eq. (39). Solution for $`\phi _\kappa `$ follows using the first expression in eq. (40). We emphasize that these solutions are valid in the super-sound-horizon (thus virtually in all scales in the matter dominated era).

Notice that the decaying solutions of $`\mathrm{\Phi }`$, $`\phi _\delta `$, and $`\phi _\kappa `$ (ignoring $`K`$ in this case) are higher order in the large-scale expansion compared with the one of $`\phi _\chi `$. In the super-horizon scale we have

$`\phi _v=\phi _\delta =\phi _\kappa =C,`$ (44)
thus time independent, whereas $`\phi _\chi `$ evolves in time according to eq. (42). As mentioned in §3.1, despite its complex behavior compared with the curvature variables in other gauges, $`\phi _\chi `$ most closely resembles the behavior of the perturbed Newtonian gravitational potential . In the $`K=0=\mathrm{\Lambda }`$ background dominated by an ideal fluid with constant $`w`$, in the super-sound horizon, the growing mode behaves as

$`\phi _\chi ={\displaystyle \frac{3+3w}{5+3w}}C,`$ (45)
thus $`\phi _\chi =\frac{2}{3}\phi _v`$ and $`\frac{3}{5}\phi _v`$ in the radiation and the matter dominated eras, respectively. We show the evolutions in Figs. 2 and 3. The amplitudes of $`\phi _v`$, $`\phi _\kappa `$ and $`\phi _\chi `$ drop as the scale comes inside the sound-horizon; compare Figs. 2 and 3.

In the sub-horizon scale in matter dominated era we have the following behaviors:

(i) $`\phi _v`$ remains conserved.

(ii) $`\phi _\delta `$ approaches $`\frac{1}{3}\delta _v`$. This is because the second term in eq. (39) dominates the first, thus $`\phi _\delta (\phi _v+\frac{1}{3(1+w)}\delta _v)\frac{1}{3}\delta _v`$ (in the sub-horizon $`K`$ term is not important); it also follows from eq. (38).

(iii) $`\phi _\kappa `$ approaches $`\phi _\chi `$. This follows from the second expression in eq. (40) and using eq. (26).

We compared these behaviors in Fig. 2. Oscillatory behaviors of $`\phi _v`$ and others occur as the scale comes inside horizon during the radiation dominated era where the sound-horizon scale is comparable to the (visual) horizon scale. Since the above analytic solutions are valid for the super-sound-horizon we do not anticipate the above behaviors in (ii) and (iii) to be valid inside horizon in the radiation dominated era. The numerical results confirm (ii) and (iii) even in such a situation: Case (iii) is more understandable by examining the second expression in eq. (40) and eq. (26) without using the large-scale solution.

Although many of the results in this section were presented in our own previous work in , which were mainly concerned with the single component fluid case, we emphasize that by keeping the most general form of fluid quantities the equations are valid even in the context of the multiple fluids and fields as well. Figures 2 and 3 show how the conserved quantities in the single component case are affected by the multi-component nature of the more realistic situations.

### 3.3 Isocurvature perturbations

We emphasize that eqs. (23-40) are generally valid for the multi-component situation. In the multi-component case $`e`$ has additional contributions besides the entropic perturbation from the individual component. Following a pioneering study by Kodama and Sasaki, using eq. (19) we decompose

$`e=e_{\mathrm{rel}}+e_{\mathrm{int}},e_{\mathrm{int}}{\displaystyle \underset{k}{}}e_{(k)},`$
$`e_{\mathrm{rel}}{\displaystyle \underset{k}{}}(c_{(k)}^2c_s^2)\delta \mu _{(k)}={\displaystyle \frac{1}{2}}{\displaystyle \underset{k,l}{}}{\displaystyle \frac{(\mu _{(k)}+p_{(k)})(\mu _{(l)}+p_{(l)})}{\mu +p}}(c_{(k)}^2c_{(l)}^2)S_{(kl)}+{\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}q_{(k)}c_{(k)}^2\delta \mu ,`$ (46)
where we introduced the following combinations

$`S_{(ij)}{\displaystyle \frac{\delta \mu _{(i)}}{\mu _{(i)}+p_{(i)}}}{\displaystyle \frac{\delta \mu _{(j)}}{\mu _{(j)}+p_{(j)}}}.`$ (47)
$`S_{(ij)}`$ is gauge-invariant for $`q_{(i)}=0`$, see eq. (21). Our definition differs from a gauge-invariant definition based on the comoving gauge ($`v/k=0`$) originally introduced in , coinciding only for $`q_{(i)}=0`$; for another gauge-invariant definition see .

Thus, we notice that the RHSs of eqs. (34,35) vanish for the following situations:

(i) single ideal fluid ($`e=0=\pi ^{(s)}`$),

(ii) multiple pressureless ideal fluids, thus $`c_{(i)}^2=0`$,

(iii) multiple ideal fluids ($`e_{(i)}=0=\pi _{(i)}`$) with the $`c_{(i)}^2=c_{(j)}^2`$, implying $`c_{(i)}^2=c_s^2`$, thus $`e_{\mathrm{rel}}=0`$.

(iv) multiple ideal fluids with adiabatic perturbation, i.e., $`S_{(ij)}=0`$ for $`q_{(i)}=0`$; for $`q_{(i)}0`$ it is not meaningful to set $`S_{(ij)}=0`$ without imposing the gauge condition.

In these situations the large-scale asymptotic solutions in eqs. (41,42) remain valid, thus $`\mathrm{\Phi }`$ is generally conserved in the super-sound-horizon scale considering general $`K`$, $`\mathrm{\Lambda }`$ and time-varying $`p(\mu )`$. The case (iv) shows that the strictly adiabatic mode ($`S_{(ij)}=0`$) is conserved in the large-scale limit for general number of fluid components.

In terms of $`S_{(ij)}`$ and a gauge-invariant combination $`v_{(ij)}(v_{(i)}v_{(j)})`$, from eqs. (20,16) we can derive

$`\dot{S}_{(ij)}+{\displaystyle \frac{3}{2}}H\{[q_{(i)}(1+c_{(i)}^2)+q_{(j)}(1+c_{(j)}^2)]S_{(ij)}`$
$`+[q_{(i)}(1+c_{(i)}^2)q_{(j)}(1+c_{(j)}^2)]{\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}(S_{(ik)}+S_{(jk)})\}`$
$`={\displaystyle \frac{k}{a}}v_{(ij)}3H\left[q_{(i)}(1+c_{(i)}^2)q_{(j)}(1+c_{(j)}^2)\right]{\displaystyle \frac{\delta \mu }{\mu +p}}+3H(q_{(i)}q_{(j)})\alpha 3He_{(ij)}+\delta Q_{(ij)},`$ (48)
$`\dot{v}_{(ij)}+Hv_{(ij)}{\displaystyle \frac{3}{2}}H\{[c_{(i)}^2+c_{(j)}^2q_{(i)}(1+c_{(i)}^2)q_{(j)}(1+c_{(j)}^2)]v_{(ij)}`$
$`+[c_{(i)}^2c_{(j)}^2q_{(i)}(1+c_{(i)}^2)+q_{(j)}(1+c_{(j)}^2)]{\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}(v_{(ik)}+v_{(jk)})\}`$
$`=3H[c_{(i)}^2c_{(j)}^2q_{(i)}(1+c_{(i)}^2)+q_{(j)}(1+c_{(j)}^2)]v+{\displaystyle \frac{k}{a}}[(c_{(i)}^2c_{(j)}^2){\displaystyle \frac{\delta \mu }{\mu +p}}+{\displaystyle \frac{1}{2}}(c_{(i)}^2+c_{(j)}^2)S_{(ij)}`$
$`+{\displaystyle \frac{1}{2}}(c_{(i)}^2c_{(j)}^2){\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}(S_{(ik)}+S_{(jk)})+e_{(ij)}{\displaystyle \frac{2}{3}}{\displaystyle \frac{k^23K}{k^2}}\pi _{(ij)}^{(s)}J_{(ij)}],`$ (49)
where $`e_{(ij)}`$, $`\pi _{(ij)}^{(s)}`$, $`\delta Q_{(ij)}`$, and $`J_{(ij)}`$ are defined similarly as $`S_{(ij)}`$:

$`e_{(ij)}{\displaystyle \frac{e_{(i)}}{\mu _{(i)}+p_{(i)}}}{\displaystyle \frac{e_{(j)}}{\mu _{(j)}+p_{(j)}}},\pi _{(ij)}^{(s)}{\displaystyle \frac{\pi _{(i)}^{(s)}}{\mu _{(i)}+p_{(i)}}}{\displaystyle \frac{\pi _{(j)}^{(s)}}{\mu _{(j)}+p_{(j)}}},`$
$`\delta Q_{(ij)}{\displaystyle \frac{\delta Q_{(i)}}{\mu _{(i)}+p_{(i)}}}{\displaystyle \frac{\delta Q_{(j)}}{\mu _{(j)}+p_{(j)}}},J_{(ij)}{\displaystyle \frac{J_{(i)}}{\mu _{(i)}+p_{(i)}}}{\displaystyle \frac{J_{(j)}}{\mu _{(j)}+p_{(j)}}}.`$ (50)
Although $`S_{(ij)}`$, $`\delta Q_{(ij)}`$ and $`J_{(ij)}`$ are not gauge-invariant for nonvanishing $`q_{(i)}`$, eqs. (48,49) are written in a gauge-ready form; these equations were derived in eqs. (35,36) of . Under the comoving gauge condition based on $`v/k=0`$, removing $`\alpha _v`$ using eq. (18), we can show that eqs. (48,49) reduce to (the corrected) eqs. (5.53,5.57) in ; typographical errors in  were corrected in eqs. (A.37,A.38) of . We emphasize that, since eqs. (15,16) are valid in the context of generalized gravity theories considered in , eqs. (48,49) in the same forms remain valid in the context of the generalized gravity theories as well.

In the case of multiple field system, we have nonvanishing entropic perturbation of the individual component. Later we will show that the relative entropic perturbation $`e_{(ij)}`$ causes a curious situation of making eq. (49) simple, and in fact making it identically satisfied, and eq. (48) alone gives the second order system; see eqs. (78-80).

For $`q_{(i)}=0`$, $`S_{(ij)}`$ is gauge-invariant. In such a case eqs. (48,49,28,46) form a closed set of equations for $`S_{(ij)}`$ and $`\delta _v`$; $`\delta \mu `$ and $`v`$ term in eq. (49) combine to give $`\delta \mu _v`$. In this way, $`\delta \mu _v`$ and $`S_{(ij)}`$ may provide the curvature and isocurvature perturbations, respectively. Using eqs. (23,33) we can express $`\delta \mu _v`$ in terms $`\mathrm{\Phi }`$ and $`\phi _\chi `$ as

$`\delta \mu _v={\displaystyle \frac{k^23K}{4\pi Ga^2}}\phi _\chi ={\displaystyle \frac{13K/k^2}{c_s^2}}\left({\displaystyle \frac{\mu +p}{H}}\dot{\mathrm{\Phi }}+e{\displaystyle \frac{2}{3}}\pi ^{(s)}\right),`$ (51)
where $`e`$ is given in eq. (46), and eqs. (34,35) provide the equation for $`\mathrm{\Phi }`$ and $`\phi _\chi `$. In this way, $`\mathrm{\Phi }`$ (or $`\phi _\chi `$) and $`S_{(ij)}`$ may provide the curvature and isocurvature perturbations. Our use of these nomenclature is somewhat misleading in the sense that, actually, the conditions $`S_{(ij)}=0`$ make the perturbation adiabatic, and the condition $`\mathrm{\Phi }=0`$ (or $`\phi _\chi =0`$) makes the perturbations isocurvature. \[In the literature, often $`\phi _\chi `$ is chosen to indicate the curvature mode. $`\phi _\chi `$ is simply related to $`\delta \mu _v`$ as in eqs. (23), and it is related to $`\mathrm{\Phi }`$ as in eq. (32).\] We consider $`\mathrm{\Phi }`$ (or $`\phi _\chi `$) and $`S_{(ij)}`$ characterize the curvature (adiabatic) and isocurvature perturbations, respectively. For nonvanishing $`q_{(i)}`$ we still have a freedom to impose a temporal gauge-condition depending on the situation, and we can use eqs. (10-14) to replace the metric and matter variables in terms of the relevant curvature and isocurvature perturbation variables. For $`n`$-component fluids we can derive $`n`$-set of coupled second-order differential equations.

In the case we could ignore the mutual interaction among the fluids for the background (thus, set $`q_{(i)}=0`$) eqs. (48,49) give

$`\ddot{S}_{(ij)}+\left[2{\displaystyle \frac{3}{2}}(c_{(i)}^2+c_{(j)}^2)\right]H\dot{S}_{(ij)}{\displaystyle \frac{3}{2}}H(c_{(i)}^2c_{(j)}^2){\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}\left(\dot{S}_{(ik)}+\dot{S}_{(jk)}\right)`$
$`+{\displaystyle \frac{1}{2}}{\displaystyle \frac{k^2}{a^2}}\left[(c_{(i)}^2+c_{(j)}^2)S_{(ij)}+(c_{(i)}^2c_{(j)}^2){\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}\left(S_{(ik)}+S_{(jk)}\right)\right]`$
$`={\displaystyle \frac{k^2}{a^2}}(c_{(i)}^2c_{(j)}^2){\displaystyle \frac{\delta \mu _v}{\mu +p}}\left(3He_{(ij)}\delta Q_{(ij)}\right)^{}\left[2{\displaystyle \frac{3}{2}}(c_{(i)}^2+c_{(j)}^2)\right]H\left(3He_{(ij)}\delta Q_{(ij)}\right)`$
$`+{\displaystyle \frac{3}{2}}H(c_{(i)}^2c_{(j)}^2){\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}\left[3H(e_{(ik)}+e_{(jk)})\delta Q_{(ik)}\delta Q_{(jk)}\right]`$
$`{\displaystyle \frac{k^2}{a^2}}\left(e_{(ij)}{\displaystyle \frac{2}{3}}{\displaystyle \frac{k^23K}{k^2}}\pi _{(ij)}^{(s)}J_{(ij)}\right),`$ (52)
where $`\delta \mu _v`$ is given in eq. (51); equations for $`\delta \mu _v`$, $`\mathrm{\Phi }`$ and $`\phi _\chi `$ are given in eqs. (28,34,35). Assuming no mutual interaction among fluids (thus, $`\delta Q_{(i)}=0=J_{(i)}`$) and ideal fluids (thus, $`e_{(i)}=0=\pi _{(i)}^{(s)}`$), the term in the RHS of eq. (52) becomes

$`{\displaystyle \frac{k^2}{a^2}}(c_{(i)}^2c_{(j)}^2){\displaystyle \frac{\delta \mu _v}{\mu +p}},`$ (53)
which vanishes for

(i) pressureless fluids,

(ii) fluids with $`c_{(i)}^2=c_{(j)}^2`$,

(iii) $`\delta \mu _v`$ = 0,

(iv) most importantly, by comparing eqs. (52,51,28,46) or eqs. (52,51,34,46), we notice that, for $`K=0`$, the RHS becomes negligible in the large-scale limit. Therefore, in the large-scale limit the curvature mode in the RHS of eq. (52) does not contribute to the isocurvature modes, thus, isocurvature modes decouple from the curvature one; this result was known by Kodama and Sasaki, see eq. (2.16) in ; for another argument in a two-component situation, see . This should be compared with the scalar field system where we expect less strong decoupling in general; see below eq. (83) and .

In the cases of (i) and (ii), assuming ideal fluids ($`e_{(ij)}=0=\pi _{(ij)}^{(s)}`$) and without interactions among fluids ($`\delta Q_{(ij)}=0=J_{(ij)}`$), eq. (52) becomes

$`\ddot{S}_{(ij)}+\left(23c_s^2\right)H\dot{S}_{(ij)}+c_s^2{\displaystyle \frac{k^2}{a^2}}S_{(ij)}=0,`$ (54)
where $`c_s^2c_{(i)}^2`$. In such a case the RHS of eq. (28) vanishes, thus we have a general large-scale solution for $`\delta _v`$

$`\delta _v={\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_{(k)}\delta _{(k)v}=(k^23K){\displaystyle \frac{H}{a^3\mu }}\left[C(k){\displaystyle ^t}{\displaystyle \frac{a(\mu +p)}{H^2}}𝑑t+{\displaystyle \frac{1}{4\pi G}}d(k)\right],`$ (55)
where the coefficients are matched using eqs. (23,42). If $`w=w_{(i)}=\mathrm{constant}`$ we have $`c_s^2=w`$, in the large-scale limit the general solution of eq. (54) becomes

$`S_{(ij)}=C_{ij}(k)+D_{ij}(k){\displaystyle ^t}{\displaystyle \frac{dt}{a^{23w}}}.`$ (56)
We note that although equations in terms of $`\delta _v`$ and $`S_{(ij)}`$ are decoupled in this case, the equations directly in terms of $`\delta _{(i)}`$ are generally coupled even in the same situation, see eq. (61) as an example; similar situations in the scalar field system can be found below eq. (90).

In the two-component ideal fluid system without mutual interactions, using $`SS_{(12)}`$, eq. (52) becomes

$`\ddot{S}+H\left(23c_z^2\right)\dot{S}+c_z^2{\displaystyle \frac{k^2}{a^2}}S={\displaystyle \frac{k^2}{a^2}}(c_{(1)}^2c_{(2)}^2){\displaystyle \frac{\delta \mu _v}{\mu +p}},`$
$`c_z^2{\displaystyle \frac{c_{(2)}^2(\mu _{(1)}+p_{(1)})+c_{(1)}^2(\mu _{(2)}+p_{(2)})}{\mu +p}}.`$ (57)
Various asymptotic solutions in the case of radiation (ideal fluid with $`p_{(r)}=\frac{1}{3}\mu _{(r)}`$) plus dust (pressureless ideal fluid with $`p_{(c)}=0`$) system were studied thoroughly by Kodama and Sasaki in .

The adiabatic and the various isocurvature initial conditions in a realistic four component system with the cold dark matter, the massless neutrinos and the tighly coupled baryon and photon in the radiation dominated limit are presented in .

### 3.4 Component equations

In this subsection we assume no interaction among fluids, thus $`q_{(i)}=0`$ and $`\delta Q_{(i)}=0=J_{(i)}`$. We have

$`\phi _{\delta _{(i)}}\phi H{\displaystyle \frac{\delta \mu _{(i)}}{\dot{\mu }_{(i)}}}=\phi +{\displaystyle \frac{\delta \mu _{(i)}}{3(\mu _{(i)}+p_{(i)})}},`$ (58)
where in the second step we used $`q_{(i)}=0`$. From eqs. (10,15) we can derive

$`\dot{\phi }_{\delta _{(i)}}={\displaystyle \frac{k}{3a}}v_{(i)\chi }{\displaystyle \frac{He_{(i)}}{\mu _{(i)}+p_{(i)}}}.`$ (59)
Using eqs. (13,16) we can derive

$`\ddot{\phi }_{\delta _{(i)}}+H(23c_{(i)}^2)\dot{\phi }_{\delta _{(i)}}+c_{(i)}^2{\displaystyle \frac{k^2}{a^2}}\phi _{\delta _{(i)}}=\left({\displaystyle \frac{1}{3}}+c_{(i)}^2\right){\displaystyle \frac{k^2}{a^2}}\phi _\chi `$
$`\left({\displaystyle \frac{He_{(i)}}{\mu _{(i)}+p_{(i)}}}\right)^{}\left[(23c_{(i)}^2)H^2+{\displaystyle \frac{k^2}{3a^2}}\right]{\displaystyle \frac{e_{(i)}}{\mu _{(i)}+p_{(i)}}}+{\displaystyle \frac{8\pi G}{3}}\pi ^{(s)}+{\displaystyle \frac{2}{9}}{\displaystyle \frac{k^23K}{a^2}}{\displaystyle \frac{\pi _{(i)}^{(s)}}{\mu _{(i)}+p_{(i)}}}.`$ (60)
Using $`\phi _{\delta _{(i)}}=\frac{\delta _{(i)\phi }}{3(1+w_{(i)})}`$ eq. (60) can be directly converted to a set of equations for $`\delta _{(i)\phi }`$, i.e., equations for $`\delta _{(i)}`$ in the uniform-curvature gauge. We plot the evolution of $`\phi _{\delta _{(c)}}`$ for the cold dark matter $`(c)`$ in Figs. 2 and 3. During the super-horizon scale $`\phi _{\delta _{(c)}}`$ and $`\phi _\delta `$ are conserved like $`\phi _v`$ and others. As they come inside the horizon they are dominated by the density parts; in the matter dominated era using eqs. (60,38) we can show $`\phi _{\delta _{(c)}}`$ and $`\phi _\delta `$ behave as $`\frac{1}{3}\delta _v\frac{1}{3}\delta _{(c)v}\frac{1}{3}\delta _{(c)v_{(c)}}`$.

It is also convenient to have equations for the individual component of density perturbation in a gauge-invariant combination with the velocity of the same component, i.e., in terms of $`\delta _{(i)v_{(i)}}\delta _{(i)}+3H(a/k)(1+w_{(i)})v_{(i)}`$. By evaluating eqs. (14,15,16) in the $`v_{(i)}=0`$ gauge condition we can derive

$`\ddot{\delta }_{(i)v_{(i)}}+\left[2+3(c_{(i)}^22w_{(i)})\right]H\dot{\delta }_{(i)v_{(i)}}+\left[c_{(i)}^2{\displaystyle \frac{k^2}{a^2}}3\dot{H}(w_{(i)}+c_{(i)}^2)+3H^2(3c_{(i)}^25w_{(i)})\right]\delta _{(i)v_{(i)}}`$
$`4\pi G(1+w_{(i)}){\displaystyle \underset{k}{}}\mu _{(k)}(1+3c_{(k)}^2)\delta _{(k)v_{(i)}}`$
$`={\displaystyle \frac{12\pi G}{\mu _{(i)}}}{\displaystyle \underset{k}{}}\left[(\mu _{(i)}+p_{(i)})e_{(k)}\left(\mu _{(k)}+p_{(k)}\right)e_{(i)}\right]{\displaystyle \frac{k^23K}{a^2}}{\displaystyle \frac{e_{(i)}}{\mu _{(i)}}}`$
$`{\displaystyle \frac{k^23K}{k^2}}\left[2{\displaystyle \frac{1+w_{(i)}}{a^2H}}\left({\displaystyle \frac{a^2H^2\pi _{(i)}^{(s)}}{\mu _{(i)}+p_{(i)}}}\right)^{}{\displaystyle \frac{2}{3}}{\displaystyle \frac{k^2}{a^2}}{\displaystyle \frac{\pi _{(i)}^{(s)}}{\mu _{(i)}}}\right],`$ (61)
with

$`4\pi G(1+w_{(i)}){\displaystyle \underset{k}{}}\mu _{(k)}(1+3c_{(k)}^2)\delta _{(k)v_{(i)}}=4\pi G(1+w_{(i)}){\displaystyle \underset{k}{}}\mu _{(k)}(1+3c_{(k)}^2)\delta _{(k)v_{(k)}}+{\displaystyle \frac{12\pi GH}{3\dot{H}k^2/a^2}}{\displaystyle \underset{k}{}}`$
$`\times \mu _{(k)}(1+3c_{(k)}^2)\left[(1+w_{(k)})\left(\dot{\delta }_{(i)v_{(i)}}3Hw_{(i)}\delta _{(i)v_{(i)}}\right)(1+w_{(i)})\left(\dot{\delta }_{(k)v_{(k)}}3Hw_{(k)}\delta _{(k)v_{(k)}}\right)\right],`$ (62)
where we used

$`\left(3\dot{H}{\displaystyle \frac{k^2}{a^2}}\right){\displaystyle \frac{a}{k}}v_{(ij)}={\displaystyle \frac{1}{1+w_{(i)}}}\left(\dot{\delta }_{(i)v_{(i)}}3Hw_{(i)}\delta _{(i)v_{(i)}}\right){\displaystyle \frac{1}{1+w_{(j)}}}\left(\dot{\delta }_{(j)v_{(j)}}3Hw_{(j)}\delta _{(j)v_{(j)}}\right),`$ (63)
which follows from eq. (15) using $`\kappa _{v_{(i)}}=\kappa +(3\dot{H}k^2/a^2)(a/k)v_{(i)}`$. We have

$`\delta _v={\displaystyle \frac{1}{\mu }}{\displaystyle \underset{k}{}}\mu _{(k)}\delta _{(k)v_{(k)}},S_{(ij)}={\displaystyle \frac{\delta \mu _{(i)v_{(i)}}}{\mu _{(i)}+p_{(i)}}}{\displaystyle \frac{\delta \mu _{(j)v_{(j)}}}{\mu _{(j)}+p_{(j)}}}3H{\displaystyle \frac{a}{k}}v_{(ij)}.`$ (64)
Equations (61,62) provide a closed set of second-order differential equations in terms of $`\delta _{(i)v_{(i)}}`$’s. This includes the case with one additional scalar field, see eq. (76) where $`\delta \mu _{(i)v_{(i)}}=\delta \mu _{(i)\delta \varphi _{(i)}}=\dot{\varphi }_{(i)}^2\alpha _{\delta \varphi _{(i)}}`$; for multiple scalar fields we also have nonvanishing $`\delta Q_{(i)}`$ in eq. (77). In the single component situation $`\delta _{(i)v_{(i)}}`$ becomes $`\delta _v`$ which behaves similarly as in the Newtonian situation . In the single component limit eq. (61) reduces to eq. (28). Using eqs. (10-18) one can also derive a closed form equations in terms of $`\delta _{(i)v}`$.

For pressureless ideal fluids ($`w_{(i)}=0`$ and $`e_{(i)}=0=\pi _{(i)}^{(s)}`$) eq. (61) gives

$`\ddot{\delta }_{(i)}+2H\dot{\delta }_{(i)}+c_{(i)}^2{\displaystyle \frac{k^2}{a^2}}\delta _{(i)}4\pi G{\displaystyle \underset{k}{}}\mu _{(k)}\delta _{(k)}=0,`$ (65)
under the $`v_{(i)}=0`$ gauge condition; thus the complete set of equations is made under mixed gauge conditions. Using eqs. (10-18) we can show that eq. (65) is valid under the $`v_{(j)}=0`$ or $`v=0`$ gauge conditions as well. Thus we can regard eq. (65) is valid under any single comoving gauge condition of $`v_{(j)}=0`$ or $`v=0`$; the $`v_{(c)}=0`$ gauge is the same as the synchronous gauge without remaining gauge mode. In fact, one can derive eq. (65) in the Newtonian context.

In the case of ideal fluids with one component a pressureless fluid ($`w_{(c)}=0`$), eq. (61) for $`i=c`$ becomes

$`\ddot{\delta }_{(c)v_{(c)}}+2H\dot{\delta }_{(c)v_{(c)}}4\pi G{\displaystyle \underset{k}{}}\mu _{(k)}(1+3c_{(k)}^2)\delta _{(k)v_{(k)}}=0.`$ (66)
Solutions in the case of radiation and dust were studied in .

## 4 Scalar field system

Considering our previous study of the multiple scalar field system in  the following can be considered as a supplement to . Besides considering general $`K`$ and $`\mathrm{\Lambda }`$ the equations are presented in the context of fluid formulation using fluid quantities; in fact, part of the fluid formulation was presented in Sec. 4.4.2 of . By presenting the scalar field system together with the real fluid system we can compare the similarity and difference of the two systems more easily.

A system of minimally coupled scalar fields can be re-interpreted as a system of fluids with special fluid quantities. For an arbitrary number of minimally coupled scalar fields we have

$`T_{ab}={\displaystyle \underset{k}{}}\left(\varphi _{(k),a}\varphi _{(k),b}{\displaystyle \frac{1}{2}}g_{ab}\varphi _{(k)}^{;c}\varphi _{(k),c}\right)Vg_{ab},`$ (67)
where $`\varphi _{(i)}`$ indicates $`i`$-th scalar field with $`i,j,k,\mathrm{}=1,2,\mathrm{}n`$; $`V=V(\varphi _{(k)})=V(\varphi _{(1)},\varphi _{(2)},\mathrm{},\varphi _{(n)})`$. Additionally we have equations of motion for the scalar fields as

$`\varphi _{(i)a}^{;a}V_{,(i)}=0,`$ (68)
where $`V_{,(i)}V/(\varphi _{(i)})`$.

We decompose the scalar fields as $`\varphi _{(i)}=\overline{\varphi }_{(i)}+\delta \varphi _{(i)}`$. To the background order, from eqs. (5,67,68) we have

$`\mu ={\displaystyle \frac{1}{2}}{\displaystyle \underset{(k)}{}}\dot{\varphi }_{(k)}^2+V,p={\displaystyle \frac{1}{2}}{\displaystyle \underset{(k)}{}}\dot{\varphi }_{(k)}^2V,`$ (69)
$`\ddot{\varphi }_{(i)}+3H\dot{\varphi }_{(i)}+V_{,(i)}=0.`$ (70)
Equations (8,9,69,70) provide a complete set for the background evolution. In order to match with eq. (9) it is convenient to introduce fluid quantities of the individual field as

$`\mu _{(i)}+p_{(i)}\dot{\varphi }_{(i)}^2,\dot{\mu }_{(i)}3H\left(\mu _{(i)}+p_{(i)}\right).`$ (71)
In this way we have $`Q_{(i)}=0`$ by definition. From eq. (71) we have $`\dot{\mu }_{(i)}=(\ddot{\varphi }_{(i)}+V_{,(i)})\dot{\varphi }_{(i)}`$ and $`\dot{p}_{(i)}=(\ddot{\varphi }_{(i)}V_{,(i)})\dot{\varphi }_{(i)}`$.

To the perturbed order, from eqs. (5,67) we have

$`\delta \mu ={\displaystyle \underset{k}{}}\left(\dot{\varphi }_{(k)}\delta \dot{\varphi }_{(k)}\dot{\varphi }_{(k)}^2\alpha +V_{,(k)}\delta \varphi _{(k)}\right),\delta p={\displaystyle \underset{k}{}}\left(\dot{\varphi }_{(k)}\delta \dot{\varphi }_{(k)}\dot{\varphi }_{(k)}^2\alpha V_{,(k)}\delta \varphi _{(k)}\right),`$
$`(\mu +p)v={\displaystyle \frac{k}{a}}{\displaystyle \underset{k}{}}\dot{\varphi }_{(k)}\delta \varphi _{(k)},\pi ^{(s)}=0.`$ (72)
In the mixture of the fluids and fields, the effective fluid quantities of the scalar fields in eqs. (69,72) should be added to the total fluid quantities in eq. (6). Equation (68) gives

$`\delta \ddot{\varphi }_{(i)}+3H\delta \dot{\varphi }_{(i)}+{\displaystyle \frac{k^2}{a^2}}\delta \varphi _{(i)}+{\displaystyle \underset{j}{}}V_{,(i)(j)}\delta \varphi _{(j)}=\dot{\varphi }_{(i)}\left(\kappa +\dot{\alpha }\right)+\left(2\ddot{\varphi }_{(i)}+3H\dot{\varphi }_{(i)}\right)\alpha .`$ (73)
Equations (10-14,72,73) provide a complete set. Under the gauge transformation we have

$`\delta \stackrel{~}{\varphi }_{(i)}=\delta \varphi _{(i)}\dot{\varphi }_{(i)}\xi ^t.`$ (74)
Thus, in a single component case, the uniform-field gauge $`\delta \varphi 0`$ gives $`v=0`$ which is the comoving gauge condition .

In order to match with eqs. (15,16) it is convenient to introduce fluid quantities of the individual field as

$`\delta \mu _{(i)}\dot{\varphi }_{(i)}\delta \dot{\varphi }_{(i)}\dot{\varphi }_{(i)}^2\alpha +V_{,(i)}\delta \varphi _{(i)},\delta p_{(i)}\dot{\varphi }_{(i)}\delta \dot{\varphi }_{(i)}\dot{\varphi }_{(i)}^2\alpha V_{,(i)}\delta \varphi _{(i)},v_{(i)}{\displaystyle \frac{k}{a}}{\displaystyle \frac{\delta \varphi _{(i)}}{\dot{\varphi }_{(i)}}},\pi _{(i)}^{(s)}0.`$ (75)
In this way, we can show $`e_{(i)}`$ in eq. (19) becomes

$`e_{(i)}=(1c_{(i)}^2)\delta \mu _{(i)v_{(i)}}.`$ (76)
With these definitions of fluid quantities eqs. (15,16) are satisfied with the following results

$`\delta Q_{(i)}={\displaystyle \underset{k}{}}V_{,(i)(k)}\left(\dot{\varphi }_{(k)}\delta \varphi _{(i)}\dot{\varphi }_{(i)}\delta \varphi _{(k)}\right),J_{(i)}=0,`$ (77)
where we used eq. (73) to derive eq. (15). From eq. (76) we can derive

$`e_{(ij)}={\displaystyle \frac{1}{2}}(2c_{(i)}^2c_{(j)}^2)\left(S_{(ij)}+3H{\displaystyle \frac{a}{k}}v_{(ij)}\right)`$
$`(c_{(i)}^2c_{(j)}^2)\left\{{\displaystyle \frac{\delta \mu }{\mu +p}}+3H{\displaystyle \frac{a}{k}}v+{\displaystyle \frac{1}{2}}{\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}\left[S_{(ik)}+S_{(jk)}+3H{\displaystyle \frac{a}{k}}\left(v_{(ik)}+v_{(jk)}\right)\right]\right\}.`$ (78)
Thus, eq. (49) simplifies to become $`\dot{v}_{(ij)}2Hv_{(ij)}=(k/a)S_{(ij)}`$ which is satisfied identically; these were found in . Introducing a gauge-invariant combination

$`\delta \varphi _{(ij)}{\displaystyle \frac{\delta \varphi _{(i)}}{\dot{\varphi }_{(i)}}}{\displaystyle \frac{\delta \varphi _{(j)}}{\dot{\varphi }_{(j)}}},`$ (79)
we have $`v_{(ij)}=(k/a)\delta \varphi _{(ij)}`$, and eq. (48) leads to

$`\delta \ddot{\varphi }_{(ij)}{\displaystyle \frac{3}{2}}H\left[(c_{(i)}^2+c_{(j)}^2)\delta \dot{\varphi }_{(ij)}+(c_{(i)}^2c_{(j)}^2){\displaystyle \underset{k}{}}{\displaystyle \frac{\mu _{(k)}+p_{(k)}}{\mu +p}}\left(\delta \dot{\varphi }_{(ik)}+\delta \dot{\varphi }_{(jk)}\right)\right]+\left(3\dot{H}+{\displaystyle \frac{k^2}{a^2}}\right)\delta \varphi _{(ij)}\delta Q_{(ij)}`$
$`=3H(c_{(i)}^2c_{(j)}^2){\displaystyle \frac{\delta \mu _v}{\mu +p}},`$ (80)
where, from eq. (77) we have

$`\delta Q_{(ij)}={\displaystyle \underset{k}{}}\left(V_{,(i)(k)}{\displaystyle \frac{\dot{\varphi }_{(k)}}{\dot{\varphi }_{(i)}}}\delta \varphi _{(ik)}V_{,(j)(k)}{\displaystyle \frac{\dot{\varphi }_{(k)}}{\dot{\varphi }_{(j)}}}\delta \varphi _{(jk)}\right).`$ (81)
Using $`c_{(i)}^2=1+2V_{,(i)}/(3H\dot{\varphi }_{(i)})`$ we have

$`{\displaystyle \frac{1}{a^3\dot{\varphi }_{(i)}\dot{\varphi }_{(j)}}}\left(a^3\dot{\varphi }_{(i)}\dot{\varphi }_{(j)}\delta \dot{\varphi }_{(ij)}\right)^{}\left({\displaystyle \frac{V_{,(i)}}{\dot{\varphi }_{(i)}}}{\displaystyle \frac{V_{,(j)}}{\dot{\varphi }_{(j)}}}\right){\displaystyle \underset{k}{}}{\displaystyle \frac{\dot{\varphi }_{(k)}^2}{\mu +p}}\left(\delta \dot{\varphi }_{(ik)}+\delta \dot{\varphi }_{(jk)}\right)`$
$`+\left(3\dot{H}+{\displaystyle \frac{k^2}{a^2}}\right)\delta \varphi _{(ij)}{\displaystyle \underset{k}{}}\left(V_{,(i)(k)}{\displaystyle \frac{\dot{\varphi }_{(k)}}{\dot{\varphi }_{(i)}}}\delta \varphi _{(ik)}V_{,(j)(k)}{\displaystyle \frac{\dot{\varphi }_{(k)}}{\dot{\varphi }_{(j)}}}\delta \varphi _{(jk)}\right)=2\left({\displaystyle \frac{V_{,(i)}}{\dot{\varphi }_{(i)}}}{\displaystyle \frac{V_{,(j)}}{\dot{\varphi }_{(j)}}}\right){\displaystyle \frac{\delta \mu _v}{\mu +p}}.`$ (82)
Using eq. (23), assuming $`K=0`$, eq. (82) becomes eq. (26) in ; in the present case it is valid for general $`K`$. From eqs. (19,72) we have $`e=(1c_s^2)\delta \mu _v2_kV_{,(k)}\delta \varphi _{(k)v}`$. Using eqs. (21,74) we have

$`e=(1c_s^2)\delta \mu _v{\displaystyle \frac{2}{\mu +p}}{\displaystyle \underset{k,l}{}}V_{,(k)}\dot{\varphi }_{(k)}\dot{\varphi }_{(l)}^2\delta \varphi _{(kl)}.`$
Thus, eq. (51) gives

$`\delta \mu _v={\displaystyle \frac{k^23K}{4\pi Ga^2}}\phi _\chi ={\displaystyle \frac{13K/k^2}{13(1c_s^2)K/k^2}}\left({\displaystyle \frac{\mu +p}{H}}\dot{\mathrm{\Phi }}{\displaystyle \frac{2}{\mu +p}}{\displaystyle \underset{k,l}{}}V_{,(k)}\dot{\varphi }_{(k)}\dot{\varphi }_{(l)}^2\delta \varphi _{(kl)}\right).`$ (83)
The RHS of eq. (82) can be compared with the one of eq. (52), i.e., eq. (53). As mentioned below eq. (53), in the large-scale limit isocurvature modes decouple from the curvature one for the multiple fluid system; in terms of $`\phi _\chi `$ eq. (53) shows $`k^4/a^4`$ factor. However, it is different for the multiple field system; in terms of $`\phi _\chi `$ eq. (82) shows $`k^4/a^4`$ factor. Equation (34) remains valid with vanishing $`\pi ^{(s)}`$, thus

$`{\displaystyle \frac{H^2c_A^2}{(\mu +p)a^3}}\left[{\displaystyle \frac{(\mu +p)a^3}{H^2c_A^2}}\dot{\mathrm{\Phi }}\right]^{}+c_A^2{\displaystyle \frac{k^2}{a^2}}\mathrm{\Phi }={\displaystyle \frac{2H^2c_A^2}{(\mu +p)a^3}}\left[{\displaystyle \frac{a^3}{(\mu +p)Hc_A^2}}{\displaystyle \underset{k,l}{}}V_{,(k)}\dot{\varphi }_{(k)}\dot{\varphi }_{(l)}^2\delta \varphi _{(kl)}\right]^{},`$ (84)
where $`c_A^213(1c_s^2)K/k^2`$. This equation can be compared with eq. (34). Since $`k^20`$ does not imply $`c_A^2k^20`$, contrary to the ideal fluid situation, we have to be careful in examining the large-scale limit for $`K0`$ case. From eq. (28) or from eq. (61), using eq. (76), we can derive

$`{\displaystyle \frac{1+w}{a^2H}}\left[{\displaystyle \frac{H^2}{a(\mu +p)}}\left({\displaystyle \frac{a^3\mu }{H}}\delta _v\right)^{}\right]^{}+c_A^2{\displaystyle \frac{k^2}{a^2}}\delta _v={\displaystyle \frac{k^23K}{a^2}}{\displaystyle \frac{1}{\mu }}{\displaystyle \frac{2}{\mu +p}}{\displaystyle \underset{k,l}{}}V_{,(k)}\dot{\varphi }_{(k)}\dot{\varphi }_{(l)}^2\delta \varphi _{(kl)}.`$ (85)
Equations (82-85) provide a complete set in terms of the curvature mode $`\mathrm{\Phi }`$ (or $`\delta _v`$, or $`\phi _\chi `$ using eq. (23)) and isocurvature modes $`\delta \varphi _{(ij)}`$. Various analytic solutions for the curvature and the isocurvature modes in a system of scalar fields were studied in  assuming $`K=0=\mathrm{\Lambda }`$. Here we have considered general $`K`$ and $`\mathrm{\Lambda }`$ in the background FLRW world model and presented the fluid formulation of the system. Literatures on the subject can be found in ; some selected ones are in , and for recent additions, see .

### 4.1 Slow-roll: linear-order solutions

We introduce the following slow-roll parameters (for different definitions, see )

$`ϵ{\displaystyle \frac{\dot{H}}{H^2}},ϵ_i{\displaystyle \frac{\ddot{\varphi }_{(i)}}{H\dot{\varphi }_{(i)}}},ϵ_{ij}{\displaystyle \frac{V_{,(i)(j)}}{3H^2}}.`$ (86)
In §4 of  we have presented various solutions of eqs. (82,84) to the zero-th order slow-roll limit of $`ϵ_i`$, i.e., ignoring $`ϵ_i`$, thus $`V_{,(i)}=3H\dot{\varphi }_{(i)}`$. In this case eqs. (82,84) decouple from each other: see eqs. (29,30) in . We consider $`K=0=\mathrm{\Lambda }`$. The assisted inflation based on multiple fields each with an exponential potential  has $`ϵ_i=ϵ=\mathrm{constant}`$. Perturbations in the assisted inflation can be managed in simple way and were studied in .

In addition to zero-th order in $`ϵ_i`$, by further assuming zero-th order in $`ϵ`$, we have

$`{\displaystyle \frac{1}{a^3}}\left(a^3\dot{\phi }_v\right)^{}+{\displaystyle \frac{k^2}{a^2}}\phi _v=0,{\displaystyle \frac{1}{a^3}}\left(a^3\delta \dot{\varphi }_{(ij)}\right)^{}+{\displaystyle \frac{k^2}{a^2}}\delta \varphi _{(ij)}=0,`$ (87)
where we used that $`V_{,(i)(k)}`$ terms in eq. (82) are of the order $`ϵ`$ and $`ϵ_i`$. Thus, the curvature and the isocurvature modes are decoupled, and in the large-scale limit we have the general solutions

$`\phi _v(k,t)=C(k)\stackrel{~}{D}(k){\displaystyle ^t}{\displaystyle \frac{1}{a^3}}𝑑t,\delta \varphi _{(ij)}(k,t)=C_{ij}(k)D_{ij}(k){\displaystyle ^t}{\displaystyle \frac{1}{a^3}}𝑑t.`$ (88)
Notice that the non-transient solutions remain constant in time. From eqs. (22,72) we have

$`\phi _v={\displaystyle \frac{H}{\mu +p}}{\displaystyle \underset{k}{}}\dot{\varphi }_{(k)}\delta \varphi _{(k)\phi }.`$ (89)
In the single component scalar field, equation for $`\delta \varphi _\phi `$, which is $`\delta \varphi `$ in the uniform-curvature gauge, resembles most closely the scalar field equation in the given (i.e., without accompanying metric fluctustions) curved background . Equation (17) in  can be written using the slow-roll parameters as

$`\delta \ddot{\varphi }_{(i)\phi }+3H\delta \dot{\varphi }_{(i)\phi }+{\displaystyle \frac{k^2}{a^2}}\delta \varphi _{(i)\phi }={\displaystyle \underset{k}{}}\left[8\pi G\dot{\varphi }_{(i)}\dot{\varphi }_{(k)}\left(3ϵ+ϵ_i+ϵ_k\right)3H^2ϵ_{ik}\right]\delta \varphi _{(k)\phi },`$ (90)
which is still exact. Thus, notice that although equations for $`\phi _v`$ and $`\delta \varphi _{(ij)}`$ are decoupled to the zero-th order in the slow-roll parameters, equations for $`\delta \varphi _{(i)\phi }`$ are generally coupled to the same order.

Now, to the linear order in the slow-roll parameters, and considering the non-transient solutions of eq. (88), eqs. (84,82) give

$`{\displaystyle \frac{1}{a^3}}\left(a^3\dot{\phi }_v^{(1)}\right)^{}+{\displaystyle \frac{k^2}{a^2}}\phi _v^{(1)}={\displaystyle \frac{2H^2}{(\mu +p)^2}}{\displaystyle \underset{j,k}{}}{\displaystyle \frac{1}{a^3}}\left(a^3ϵ_j\right)^{}\dot{\varphi }_{(j)}^2\dot{\varphi }_{(k)}^2\delta \varphi _{(jk)}^{(0)}𝒮^{(1)},`$ (91)
$`{\displaystyle \frac{1}{a^3}}\left(a^3\delta \dot{\varphi }_{(ij)}^{(1)}\right)^{}+{\displaystyle \frac{k^2}{a^2}}\delta \varphi _{(ij)}^{(1)}=3H^2\left[ϵ\delta \varphi _{(ij)}^{(0)}+{\displaystyle \underset{k}{}}\left(ϵ_{ik}{\displaystyle \frac{\dot{\varphi }_{(k)}}{\dot{\varphi }_{(i)}}}\delta \varphi _{(ik)}^{(0)}ϵ_{jk}{\displaystyle \frac{\dot{\varphi }_{(k)}}{\dot{\varphi }_{(j)}}}\delta \varphi _{(jk)}^{(0)}\right)\right]𝒮_{ij}^{(1)},`$ (92)
where superscripts $`(0)`$ and $`(1)`$ indicate the variables to the zero-th order and the linear order in the slow-roll parameters, respectively. Apparently, we can extend this perturbative procedure to the higher order in the slow-roll parameters. Therefore, to the linear order in the slow-roll parameters the general solutions in the large-scale limit become

$`\phi _v^{(1)}={\displaystyle ^t}a^3\left({\displaystyle ^t}{\displaystyle \frac{1}{a^3}}𝑑t\right)𝒮^{(1)}𝑑t+\left({\displaystyle ^t}{\displaystyle \frac{1}{a^3}}𝑑t\right){\displaystyle ^t}a^3𝒮^{(1)}𝑑t,`$ (93)
$`\delta \varphi _{(ij)}^{(1)}={\displaystyle ^t}a^3\left({\displaystyle ^t}{\displaystyle \frac{1}{a^3}}𝑑t\right)𝒮_{ij}^{(1)}𝑑t+\left({\displaystyle ^t}{\displaystyle \frac{1}{a^3}}𝑑t\right){\displaystyle ^t}a^3𝒮_{ij}^{(1)}𝑑t.`$ (94)
We note that to the linear order in the slow-roll parameter $`ϵ_i`$ and considering the non-transient solution in the large-scale limit, the RHS of eq. (82), using eq. (83), vanishes. Thus, in such a case the isocurvature modes decouple from the curvature one, whereas the curvature mode is coupled with the isocurvature ones, see eqs. (91,92). Studies in more general situation with nonlinear sigma type coupling are presented in  keeping covariant forms in field space.

## 5 Fluid and field system

The formulations presented in §2 and 3 are valid for a mixture of arbitrary numbers of fluids and scalar fields. Thus, we will not rewrite the complete set of equations again.

As an example, in this section we consider a two-component system of a fluid and a field with a general interaction term between them. As the Lagrangian we consider

$`=\sqrt{g}\left[{\displaystyle \frac{1}{16\pi G}}R{\displaystyle \frac{1}{2}}\varphi ^{;c}\varphi _{,c}V(\varphi )+L_{(f)}\right].`$ (95)
$`L_{(f)}`$ is the fluid part Lagrangian which also depends on the scalar field as $`L_{(f)}=L_{(f)}(\mathrm{fluid},g_{ab},\varphi )`$. The variation with respect to $`g_{ab}`$ leads to eq. (1) with an identification $`T_{ab}=T_{(f)ab}+T_{(\varphi )ab}`$. The variation with respect to $`\varphi `$ leads to the equation of motion for $`\varphi `$

$`\varphi _c^{;c}V_{,\varphi }=L_{(f),\varphi }\mathrm{\Gamma }.`$ (96)
Thus, from eq. (3) we have $`Q_{(f)a}=Q_{(\varphi )a}`$ and

$`T_{(f)a;b}^bQ_{(f)a}=L_{(f),\varphi }\varphi _{,a}`$ (97)
where we used $`i,j,\mathrm{}=1,2f,\varphi `$. As we decompose $`\mathrm{\Gamma }=\overline{\mathrm{\Gamma }}+\delta \mathrm{\Gamma }`$, from eq. (7) we have

$`\overline{Q}_{(f)}=\overline{\mathrm{\Gamma }}\dot{\overline{\varphi }},\delta Q_{(f)}=\delta \mathrm{\Gamma }\dot{\varphi }+\mathrm{\Gamma }(\delta \dot{\varphi }\dot{\varphi }\alpha ),J_{(f)}=\mathrm{\Gamma }\delta \varphi .`$ (98)
We have $`\overline{Q}_{(\varphi )}=\overline{Q}_{(f)}`$ and similarly for $`\delta Q_{(\varphi )}`$ and $`J_{(\varphi )}`$.

By interpreting the fluid quantities properly, equations in §2 and 3 remain valid. For the background we set $`\mu =\mu _{(f)}+\mu _{(\varphi )}`$ and similarly for $`p`$ with

$`\mu _{(\varphi )}={\displaystyle \frac{1}{2}}\dot{\varphi }^2+V,p_{(\varphi )}={\displaystyle \frac{1}{2}}\dot{\varphi }^2V.`$ (99)
Equation (8) remains valid. Equation (9) is valid for $`(i)=(f)`$, and eq. (96) leads to

$`\ddot{\varphi }+3H\dot{\varphi }+V_{,\varphi }=\mathrm{\Gamma }.`$ (100)
For the perturbations we set $`\delta \mu =\delta \mu _{(f)}+\delta \mu _{(\varphi )}`$ and similarly for $`\delta p`$, $`(\mu +p)v`$, and $`\pi ^{(s)}`$, with

$`\delta \mu _{(\varphi )}=\dot{\varphi }\delta \dot{\varphi }\dot{\varphi }^2\alpha +V_{,\varphi }\delta \varphi ,\delta p_{(\varphi )}=\dot{\varphi }\delta \dot{\varphi }\dot{\varphi }^2\alpha V_{,\varphi }\delta \varphi ,v_{(\varphi )}={\displaystyle \frac{k}{a}}{\displaystyle \frac{\delta \varphi }{\dot{\varphi }}},\pi _{(\varphi )}^{(s)}=0.`$ (101)
Equations (10-14) are valid for collective fluid quantities. Equations (15,16) are valid for the fluid with $`(i)=(f)`$, and eq. (96) leads to

$`\delta \ddot{\varphi }+3H\delta \dot{\varphi }+{\displaystyle \frac{k^2}{a^2}}\delta \varphi +V_{,\varphi \varphi }\delta \varphi =\dot{\varphi }\left(\kappa +\dot{\alpha }\right)+\left(2\ddot{\varphi }+3H\dot{\varphi }\right)\alpha \delta \mathrm{\Gamma }.`$ (102)
After choosing the temporal gauge condition we can derive a fourth-order differential equation, or a set of coupled two second-order differential equations. As an example, we derive one form of such a set of equations based on a gauge condition $`v_{(f)}/k0`$ which is the comoving gauge based on the fluid velocity; in the following we sim ply write $`v_{(f)}=0`$ as the gauge condition. Under our gauge condition $`v_{(f)}=0`$ the perturbation variables are gauge invariant, i.e., $`\delta _{(f)}(\delta \mu _{(f)}/\mu _{(f)})=\delta _{(f)v_{(f)}}`$ and $`\delta \varphi =\delta \varphi _{v_{(f)}}`$, where, from eqs. (22,74), we have:

$`\delta _{(f)v_{(f)}}\delta _{(f)}{\displaystyle \frac{\dot{\mu }_{(f)}}{\mu _{(f)}}}{\displaystyle \frac{a}{k}}v_{(f)},\delta \varphi _{v_{(f)}}\delta \varphi \dot{\varphi }{\displaystyle \frac{a}{k}}v_{(f)}.`$ (103)
We assume an ideal fluid with $`wp_{(f)}/\mu _{(f)}=\mathrm{constant}`$. Using eq. (16) for the fluid component we can express $`\alpha `$ in terms of $`\delta _{(f)}`$ and $`\delta \varphi `$. Using eq. (15) we can express $`\kappa `$ in terms of $`\delta _{(f)}`$ and $`\delta \varphi `$. Thus, eqs. (14,102) lead to

$`\ddot{\delta }_{(f)}+\left(23w\right)H\dot{\delta }_{(f)}+\left[w{\displaystyle \frac{k^2}{a^2}}6w\left(\dot{H}+H^2\right)4\pi G\mu _{(f)}(1+w)(1+3w)\right]\delta _{(f)}`$
$`=8\pi G\left[2w\dot{\varphi }^2\delta _{(f)}+(1+w)\left(2\dot{\varphi }\delta \dot{\varphi }V_{,\varphi }\delta \varphi \right)\right]`$
$`{\displaystyle \frac{1}{a^2}}\left\{a^2\left[{\displaystyle \frac{\mathrm{\Gamma }\dot{\varphi }\delta _{(f)}}{\mu _{(f)}}}{\displaystyle \frac{\dot{\varphi }\delta \mathrm{\Gamma }}{\mu _{(f)}}}{\displaystyle \frac{\mathrm{\Gamma }}{\mu _{(f)}}}\left(\delta \dot{\varphi }+3H\delta \varphi \right)\right]\right\}^{}+\left(3\dot{H}+16\pi G\dot{\varphi }^2{\displaystyle \frac{k^2}{a^2}}\right){\displaystyle \frac{\mathrm{\Gamma }}{\mu _{(f)}}}\delta \varphi ,`$ (104)
$`\delta \ddot{\varphi }+3H\delta \dot{\varphi }+\left({\displaystyle \frac{k^2}{a^2}}+V_{,\varphi \varphi }\right)\delta \varphi ={\displaystyle \frac{1}{1+w}}\left[(1w)\dot{\varphi }\dot{\delta }_{(f)}2w\left(\ddot{\varphi }+3H\dot{\varphi }\right)\delta _{(f)}\right]+{\displaystyle \frac{\mathrm{\Gamma }\dot{\varphi }^2}{(1+w)\mu _{(f)}}}\delta _{(f)}`$
$`\left[1+{\displaystyle \frac{\dot{\varphi }^2}{(1+w)\mu _{(f)}}}\right]\delta \mathrm{\Gamma }+{\displaystyle \frac{1}{1+w}}\left[{\displaystyle \frac{\mathrm{\Gamma }\dot{\varphi }}{\mu _{(f)}}}\delta \dot{\varphi }\dot{\varphi }\left({\displaystyle \frac{\mathrm{\Gamma }\delta \varphi }{\mu _{(f)}}}\right)^{}+2\left(\ddot{\varphi }+3H\dot{\varphi }\right){\displaystyle \frac{\mathrm{\Gamma }\delta \varphi }{\mu _{(f)}}}\right].`$ (105)
Although the RHS of eq. (104) contains $`\delta \ddot{\varphi }`$ term, it can be replaced by using eq. (105). Notice that eqs. (104,105) are valid for general $`K`$ and $`\mathrm{\Lambda }`$. The LHS of eq. (104) is the familiar density perturbation in the comoving gauge derived in , and the LHS of of eq. (105) is also the familiar perturbed scalar field equation without the metric perturbation; compare with eq. (102). Thus, our $`v_{(f)}=0`$ gauge choice allows simple equations for the uncoupled parts of the system.

This set of equations or eqs. (61,76) will be useful to handle the following situations:

(i) Warm inflation scenario with nonvanishing interaction term, $`\mathrm{\Gamma }`$ and $`\delta \mathrm{\Gamma }`$, between the field and the radiation ($`w=\frac{1}{3}`$) .

(ii) Time-varying cosmological constant simulated using the scalar field, often called a quintessence ; in this case we may ignore the direct interaction between the fluid and the field $`\mathrm{\Gamma }=0=\delta \mathrm{\Gamma }`$. An application of eqs. (104,105) to a system with an exponential type field potential is made in .

## 6 Summary

In this paper we have investigated aspects of scalar-type cosmological perturbation in the context of multiple numbers of mutually interacting imperfect fluids and minimally coupled scalar fields in Einstein gravity. Equations are presented using the curvature ($`\mathrm{\Phi }`$ or $`\phi _\chi `$) and isocurvature ($`S_{(ij)}`$ or $`\delta \varphi _{(ij)}`$) perturbation variables. It looks there exists no clear consensus about the curvature/isocurvature decomposition of perturbation in the literature. Since the equations are all coupled, such a decomposition may have clearer meaning in setting up the initial conditions. In §3 and 4 we have shown that the equations can be meaningfully classified into the our definitions of the curvature ($`\mathrm{\Phi }`$ or $`\phi _\chi `$) and isocurvature ($`S_{(ij)}`$ or $`\delta \varphi _{(ij)}`$) modes; see below eq. (51). We have shown that either $`\mathrm{\Phi }`$ or $`\phi _\chi `$ can characterize the curvature mode depending on situations.

We have presented the equations in a gauge-ready form where we have a freedom to choose one temporal gauge condition from the variables in eqs. (21,74) or the linear combinations of them. In the general situation, eqs. (10-18) provide a complete set. In the case of scalar fields we also have eqs. (73,72). In terms of the curvature and isocurvature perturbations, eqs. (34,48,49,46,51) with eqs. (10-14) provide a complete set, and in the special case with multiple fields eqs. (80-84) or eqs. (61,76) provide complete sets. In §5 we have shown how to consider the case with mutual interaction among fluids and fields.

Some new results found in this paper are the following:

(i) Gauge-ready formulation of the fluid-field system in Einstein gravity, §2, 4 and 5.

(ii) Several useful forms of the individual fluid equations, eqs. (29,30) and §3.4.

(iii) Analyses of various curvature perturbations in realistic situations of multi-component system, §3.2.

(iv) Equations (48,49) are valid even in a class of generalized gravity theories, see below eq. (16) and eq. (50).

(v) Extension of fluid formulation of the multiple field system, §4.

(vi) In the scalar field system the isocurvature modes are less decoupled from the curvature mode in the large-scale limit compared with the fluids system, see below eq. (83). The couplig term vanishes to the linear order in slow-roll expansion, see below eq. (94).

(vii) Solutions for $`c_{(i)}^2=c_{(j)}^2`$ fluids, eqs. (54-56).

(viii) Solutions valid to the linear order in the slow-roll parameters, §4.1.

The sets of equations in §5 may be useful to analyze the evolution of structures in the world models with quintessence and in the warm inflation scenario. However, these sets of equations are the ones derived in certain gauge conditions, and when we encounter new problems we believe it is best to go back to the original set of equations in the gauge-ready form and see whether we could gather any new perspective from other gauge conditions as well: the set is in eqs. (10-18,73). Corresponding set of equations in a gauge-ready form which is applicable in a wide class of generalized gravity theories including contributions from the kinetic components based on Boltzmann equations is presented in . Specific applications will be made in future occasions.

### Acknowledgments

We thank Winfried Zimdahl for suggesting the study and for his useful discussions and comments during the work. We also wish to thank Stefan Groot Nibbelink and David Wands for useful comments. HN was supported by grant No. 2000-0-113-001-3 from the Basic Research Program of the Korea Science and Engineering Foundation. JH was supported by the Korea Research Foundation Grants (KRF-2000-013-DA004 and 2000-015-DP0080).