Formation and Evolution of Wakes in the Spacetime Generated by a Cosmic String in f(R) Theory of Gravity

Gabriela Coutinho de Carvalho^{1},
Maria Emília Xavier Guimarães^{2},
Pablo Oliveira Mesquita^{3},
Joaquim Lopes Neto^{3}

Show more

1. Introduction

Extended Theories of Gravity have recently received increasingly attention in issues such as dark matter [1] [2] [3] [4] and dark energy [5] [6] [7]. In particular, the $f\left(R\right)$ theories of gravity have been suggested as a possible alternative to explain the late time cosmic speed-up experienced by our universe [8] [9]. Such theories avoid the Ostrogradski’s instability that can otherwise prove to be problematic for general higher derivatives theories [10] [11]. Besides, these theories could be probed by the recent LIGO detections [12] and by gravitational lensing observations [13].

On the other hand, it is well known that several types of topological defects may have been created by the vacuum phase transitions in the early universe [14] [15]. In particular, cosmic strings have been extensively studied in many kinds of alternative gravity theories, notably as scalar-tensor and $f\left(R\right)$ theories where many aspects and applications were developed [16] [17] [18] [19].

In this paper, our main purpose is to study the formation and evolution of wakes in cylindrically symmetric solutions in the framework of $f\left(R\right)$ theories in vacuum. In particular, we aim to explore the propagation of particles and light in a $f\left(R\right)$ cosmic string. We compare our results with the wakes formed by cosmic string solutions obtained in General Relativity and the Scalar-Tensor Theories of Gravity.

This paper is structured as follows. In Section 2, we briefly review the $f\left(R\right)$ theories of gravity (for comprehensive reviews in $f\left(R\right)$ gravity see [20] ). In Section 3, for the sake of completeness, we review the obtention of the cylindrically symmetric solutions in vacuum. We give special attention to the cosmic string solution. Is it important to stress that, in this section, we follow the reference [21] ? In Section 4, the formation and evolution of wakes in the background of a $f\left(R\right)$ are studied and these are our original work. Finally, in the conclusion we summarize our main results and discuss some perspectives.

2. $f\left(R\right)$ Gravity in the Metric Formalism—A Brief Review

In this section, we will mainly follow the references [20].

The action associated with the modified theories of gravity coupled with matter fields is given by:

$S=\frac{1}{2\kappa}{\displaystyle \int {\text{d}}^{4}x\sqrt{-g}f\left(R\right)}+{S}_{m}\mathrm{,}$ (1)

where
$f\left(R\right)$ is an analytical function of the the Ricci scalar, *R*,
$\kappa =8\pi G$ and
${S}_{m}$ corresponds to the action associated with the matter fields. By using the metric formalism, the field equations become:

${G}_{\mu \nu}\equiv {R}_{\mu \nu}-\frac{1}{2}R{g}_{\mu \nu}={T}_{\mu \nu}^{c}+\kappa {\stackrel{\u02dc}{T}}_{\mu \nu}^{m},$ (2)

in which ${T}_{\mu \nu}^{c}$ is the geometric energy-momentum tensor, namely

${T}_{\mu \nu}^{c}=\frac{1}{F\left(R\right)}\left\{\frac{1}{2}{g}_{\mu \nu}\left(f\left(R\right)-F\left(R\right)R\right)+{\nabla}^{\alpha}{\nabla}^{\beta}F\left(R\right)\left({g}_{\alpha \mu}{g}_{\beta \nu}-{g}_{\mu \nu}{g}_{\alpha \beta}\right)\right\}$ (3)

with $F\left(R\right)\equiv \frac{\text{d}f\left(R\right)}{\text{d}R}$.

The standard minimally coupled energy-momentum tensor, ${T}_{\mu \nu}^{m}$, derived from the matter action, is related to ${\stackrel{\u02dc}{T}}_{\mu \nu}^{m}$ by

${\stackrel{\u02dc}{T}}_{\mu \nu}^{m}={T}_{\mu \nu}^{m}/F\left(R\right)\mathrm{.}$ (4)

Thus the field equations can be written as

$F\left(R\right){R}_{\mu \nu}-\frac{1}{2}f\left(R\right){g}_{\mu \nu}-{\nabla}_{\mu}{\nabla}_{\nu}F\left(R\right)+{g}_{\mu \nu}\square F\left(R\right)=\kappa {T}_{\mu \nu}^{m}\mathrm{.}$ (5)

Taking the trace of the above equation we get

$F\left(R\right)R-2f\left(R\right)+3\square F\left(R\right)=\kappa {T}^{m}\mathrm{,}$ (6)

which express a further scalar degree of freedom that arises in the modified theory. Through this equation it is possible to express $f\left(R\right)$ in terms of its derivatives and the trace of the matter energy-momentum tensor, as follows

$f\left(R\right)=\frac{1}{2}\left(F\left(R\right)R+3\square F\left(R\right)-\kappa {T}^{m}\right)\mathrm{.}$ (7)

Substituting the above expression into (5) we obtain

$F\left(R\right){R}_{\mu \nu}-{\nabla}_{\mu}{\nabla}_{\nu}F\left(R\right)-\kappa {\stackrel{\u02dc}{T}}_{\mu \nu}^{m}=\frac{{g}_{\mu \nu}}{4}\left[F\left(R\right)R-\text{\hspace{0.17em}}\square F\left(R\right)-\kappa {T}^{m}\right]\mathrm{.}$ (8)

From this expression we can see that the combination below

${C}_{\mu}=\frac{F\left(R\right){R}_{\mu \mu}-{\nabla}_{\mu}{\nabla}_{\mu}F\left(R\right)-\kappa {T}_{\mu \mu}^{m}}{{g}_{\mu \mu}}\mathrm{,}$ (9)

with fixed indices, is independent of the corresponding index. So, the following relation

${C}_{\mu}-{C}_{\nu}=\mathrm{0,}$ (10)

holds for all $\mu $ and $\nu $.

3. Vacuum Cylindrically Symmetric Solutions in $f\left(R\right)$ Gravity: A Brief Review

In this section we review how to derive the field equations associated with cosmic string system in the context of modified theories of gravity following Azadi *et al*. [21].

The $f\left(R\right)$ field equations in vacuum are given by

$F\left(R\right){R}_{\mu \nu}-\frac{1}{2}f\left(R\right){g}_{\mu \nu}-\left({\nabla}_{\mu}{\nabla}_{\nu}-{g}_{\mu \nu}\square \right)F\left(R\right)=0$ (11)

Taking the trace of Equation (11), we have

$F\left(R\right)R-2f\left(R\right)+3\square F\left(R\right)=0.$ (12)

Now, since we are interested in obtaining static solutions with cylindrical symmetry in vacuum, we will work with a general metric in Weyl coordinates $\left(t\mathrm{,}r\mathrm{,}\varphi \mathrm{,}z\right)$ given by

${g}_{\mu \nu}=diag\left(-{\text{e}}^{2k-2u}\mathrm{,}{\text{e}}^{2k-2u}\mathrm{,}{\omega}^{2}{\text{e}}^{-2u}\mathrm{,}{\text{e}}^{2u}\right)\mathrm{,}$ (13)

where
$k\mathrm{,}u$ and
$\omega $ are functions of *r* only.

The non-zero components of the Ricci tensor are:

$\begin{array}{l}{R}_{00}={k}^{\u2033}-{u}^{\u2033}+\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}-\frac{{u}^{\prime}{\omega}^{\prime}}{\omega}\\ {R}_{11}=-{k}^{\u2033}+{u}^{\u2033}-\frac{{\omega}^{\u2033}}{\omega}+\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}-\frac{{u}^{\prime}{\omega}^{\prime}}{\omega}-2{{u}^{\prime}}^{2}\\ {R}_{22}={\text{e}}^{-2k}\left(\omega {\omega}^{\prime}{u}^{\prime}-\omega {\omega}^{\u2033}+{\omega}^{2}{u}^{\u2033}\right)\\ {R}_{33}={\text{e}}^{4u-2k}\left(-{u}^{\u2033}-\frac{{u}^{\prime}{\omega}^{\prime}}{\omega}\right),\end{array}$ (14)

where prime (*'*) indicates derivative with respect to *r*. Therefore, the scalar curvature is

$R=-2{\text{e}}^{2u}\left(\frac{-\omega {u}^{\u2033}+\omega {k}^{\u2033}-{u}^{\prime}{\omega}^{\prime}+{\omega}^{\u2033}+\omega {{u}^{\prime}}^{2}}{\omega {\text{e}}^{2k}}\right),$ (15)

Replacing $f\left(R\right)$ of Equation (12) in Equation (11), finally we get:

$\frac{F{R}_{\mu \nu}-{\nabla}_{\mu}{\nabla}_{\nu}F}{{g}_{\mu \nu}}=\frac{1}{4}\left(FR-\text{\hspace{0.17em}}\square F\left(R\right)\right)\mathrm{.}$ (16)

Defining ${A}_{\mu}$ (which is the equivalent of ${C}_{\mu}$ in the absence of matter) as

${A}_{\mu}=\frac{F{R}_{\mu \mu}-{\nabla}_{\mu}{\nabla}_{\mu}F}{{g}_{\mu \mu}}\mathrm{,}$ (17)

we can easily see that

${A}_{\mu}=\frac{1}{4}\left(FR-\text{\hspace{0.17em}}\square F\left(R\right)\right)\mathrm{,}$ (18)

which is a scalar quantity. This means that ${A}_{\mu}={A}_{\nu}$ for any $\mu \mathrm{,}\nu $, which implies that we can replace Equation (11) by ${A}_{t}={A}_{r}$, ${A}_{t}={A}_{\varphi}$ and ${A}_{t}={A}_{z}$. After some straightforward calculation, we finally get, respectively

$-{F}^{\u2033}+2{F}^{\prime}\left({k}^{\prime}-{u}^{\prime}\right)+F\left(-\frac{2{k}^{\prime}{\omega}^{\prime}}{\omega}+\frac{{\omega}^{\u2033}}{\omega}+2{{u}^{\prime}}^{2}\right)=0$ (19)

$F{\omega}^{2}\left(-{k}^{\u2033}-\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}+\frac{{\omega}^{\u2033}}{\omega}\right)+{F}^{\prime}\left(\omega {\omega}^{\prime}-{\omega}^{2}{k}^{\prime}\right)=0$ (20)

$F\left(-{k}^{\u2033}+2{u}^{\u2033}-\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}+\frac{2{\omega}^{\prime}{u}^{\prime}}{\omega}\right)+{F}^{\prime}\left({k}^{\prime}-2{u}^{\prime}\right)=0.$ (21)

Therefore, any group of functions $k\mathrm{,}u$ and $\omega $ which satisfies the equations above is a solution of the modified gravity equations in vacuum.

Since these equations are highly nonlinear, we will consider the particular case where $R=0$. This is justified by the fact that we are interested in cosmic string solution and the external metric of a cosmic string is locally flat.

Field Equations Solutions for the Special Case *R*=0

Let us start deriving Equation (11) with respect to *r*,

$R{F}^{\prime}-{R}^{\prime}F+3{\left(\square F\right)}^{\prime}=0.$ (22)

In the particular case where *R*= constant, Equation (22) implies that
${F}^{\prime}=0$. As a consequence, Equations (17-19) become

$2{{u}^{\prime}}^{2}+\frac{{\omega}^{\u2033}}{\omega}-2\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}=0$ (23)

${k}^{\u2033}+\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}-\frac{{\omega}^{\u2033}}{\omega}=0$ (24)

$2{u}^{\u2033}+2\frac{{\omega}^{\prime}{u}^{\prime}}{\omega}-\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}-{{k}^{\prime}}^{\prime}=0$ (25)

In order to solve this non-linear equations system, let us first sum up Equations (22) and (23). In doing this, we get

${{u}^{\prime}}^{\prime}+\frac{{\omega}^{\prime}{u}^{\prime}}{\omega}-\frac{1}{2}\frac{{\omega}^{\u2033}}{\omega}=0.$ (26)

Defining ${u}^{\prime}=g\left(r\right)$ in Equation (26), we get

$\omega {g}^{\prime}+{\omega}^{\prime}g=\frac{1}{2}{\omega}^{\u2033}$ (27)

which can be rewritten as

$\frac{\text{d}\left(\omega g\right)}{\text{d}r}=\frac{1}{2}\frac{\text{d}}{\text{d}r}\left(\frac{\text{d}\omega}{\text{d}r}\right).$ (28)

Integrating Equation (28), we obtain

$g\left(r\right)={u}^{\prime}=\frac{1}{2}\frac{{\omega}^{\prime}+{c}_{2}}{\omega},$ (29)

where ${c}_{2}$ is a constant to be determined later.

Now, subtracting Equation (23) from Equation (22), we get

${k}^{\u2033}+\frac{{k}^{\prime}{\omega}^{\prime}}{\omega}-\frac{{\omega}^{\u2033}}{\omega}=0,$ (30)

and we find

${k}^{\prime}=\frac{{\omega}^{\prime}+{c}_{1}}{\omega}$ (31)

where ${c}_{1}$ is a constant to be determined later.

Replacing Equations (29) and (31) in Equation (21), we get a differential equation for the function $\omega $, which is

$\frac{1}{2}{\left(\frac{{\omega}^{\prime}+{c}_{2}}{\omega}\right)}^{2}+\frac{{\omega}^{\u2033}}{\omega}=2{\omega}^{\prime}\frac{{\omega}^{\prime}+{c}_{1}}{{\omega}^{2}}$ (32)

In order to solve our equations, we will make the hypothesis that
$\omega \left(r\right)$ is a linear function of *r*. This is justified by cosmic string solutions in either General Relativity or Scalar-tensor gravities. Therefore,
${\omega}^{\u2033}=0$. Hence, we have

$u={c}_{3}\pm \sqrt{\frac{{c}_{5}}{{c}_{6}}}\mathrm{ln}\omega $ (33)

$k={c}_{4}+\frac{{c}_{5}}{{c}_{6}}\mathrm{ln}\frac{\omega}{{c}_{6}}$ (34)

$\omega ={c}_{6}r+{c}_{7}$ (35)

In order to satisfy Equation (32), the constants ${c}_{5}$ and ${c}_{6}$ must obey the following relations

${c}_{2}=\left(2\sqrt{\frac{{c}_{5}}{{c}_{6}}}-1\right){c}_{6}$

${c}_{1}={c}_{5}-{c}_{6}.$ (36)

We can easily see that the metric functions (33-35) satisfy Equation (32) in the particular case where $R=0$.

Redefining the quantities ${\stackrel{\u02dc}{c}}_{4}={c}_{4}-\left({c}_{5}/{c}_{6}\right)\mathrm{ln}{c}_{6}$ and $\rho =\omega ={c}_{6}r$ and making ${c}_{7}=0$ without any loss of generality, we can write down the metric in Weyl coordinates as

$\text{d}{s}^{2}={\text{e}}^{-2\left({c}_{3}\pm \sqrt{\frac{{c}_{5}}{{c}_{6}}}\mathrm{ln}\rho \right)}\left[{\text{e}}^{2\left({\stackrel{\u02dc}{c}}_{4}+\frac{{c}_{5}}{{c}_{6}}\mathrm{ln}\rho \right)}\left(\frac{\text{d}{\rho}^{2}}{{c}_{6}^{2}}-\text{d}{t}^{2}\right)+{\rho}^{2}\text{d}{\varphi}^{2}\right]+{\text{e}}^{2\left({c}_{3}\pm \sqrt{\frac{{c}_{5}}{{c}_{6}}}\mathrm{ln}\rho \right)}\text{d}{z}^{2}.$ (37)

Defining the quantities $m=\sqrt{\frac{{c}_{5}}{{c}_{6}}}$ and $A=\frac{{\text{e}}^{{\stackrel{\u02dc}{c}}_{4}-{c}_{3}}}{{c}_{6}}$ and defining new coordinates such as

$\stackrel{\u02dc}{t}={c}_{6}{A}^{\frac{1}{m\left(m\mp 1\right)+1}}t$

$\stackrel{\u02dc}{\rho}={A}^{\frac{1}{m\left(m\mp 1\right)+1}}\rho $

$\stackrel{\u02dc}{\varphi}={\text{e}}^{-{c}_{3}}{A}^{-\frac{1\mp m}{m\left(m\mp 1\right)+1}}\varphi $ (38)

$\stackrel{\u02dc}{z}={\text{e}}^{{c}_{3}}{A}^{\frac{\mp m}{m\left(m\mp 1\right)+1}}z\text{\hspace{0.05em}}\text{\hspace{0.05em}}.$ (39)

In doing this, the metric reduces to a very simple form^{1}

$\text{d}{s}^{2}={\stackrel{\u02dc}{\rho}}^{2m\left(m\mp 1\right)}\left(-\text{d}{\stackrel{\u02dc}{t}}^{2}+\text{d}{\stackrel{\u02dc}{\rho}}^{2}\right)+{\stackrel{\u02dc}{\rho}}^{2\left(1\mp m\right)}\text{d}{\stackrel{\u02dc}{\varphi}}^{2}+{\stackrel{\u02dc}{\rho}}^{\pm 2m}\text{d}{\stackrel{\u02dc}{z}}^{2}$ (40)

Applying the complex transformation $\stackrel{\u02dc}{t}\to i\stackrel{\u02dc}{z}$ and $\stackrel{\u02dc}{z}\to i\stackrel{\u02dc}{t}$, we get a well known metric [22]

$\text{d}{s}^{2}={\stackrel{\u02dc}{\rho}}^{2m\left(m\mp 1\right)}\left(\text{d}{\stackrel{\u02dc}{z}}^{2}+\text{d}{\stackrel{\u02dc}{\rho}}^{2}\right)+{\stackrel{\u02dc}{\rho}}^{2\left(1\mp m\right)}\text{d}{\stackrel{\u02dc}{\varphi}}^{2}-{\stackrel{\u02dc}{\rho}}^{\pm 2m}\text{d}{\stackrel{\u02dc}{t}}^{2},$ (41)

which, apart from the sigh $\mp $ is pretty much the same as the Levi-Civita static cylindrically symmetric solution in General Relativity [23].

When
$m=0$, the spacetime (41) becomes^{2}

$\text{d}{s}^{2}=-\text{d}{t}^{2}+\text{d}{\rho}^{2}+{\rho}^{2}{\text{e}}^{-2{c}_{4}}{c}_{6}^{2}\text{d}{\varphi}^{2}+\text{d}{z}^{2}.$ (42)

It is very easy to see that this spacetime is locally flat but not globally Euclidean. This spacetime is conical, as we can see in Figure 1, with a deficit angle equal to

$\delta \varphi =2\pi \left(1-{\text{e}}^{-2{c}_{4}}{c}_{6}\right),$ (43)

as long as ${\text{e}}^{-2{c}_{4}}{c}_{6}<1$ which imposes some constraints on the constants ${c}_{4}$ and ${c}_{6}$.

4. Formation and Evolution of Wakes in the Spacetime of a $f\left(R\right)$ Cosmic String

In this section we will study the formation and evolution of some structures when a cosmic string moves through a region containing baryonic matter.

4.1. The Formation of Wakes

Let us suppose that the cosmic string moves with a constant velocity
${v}_{c}$ in the *x*-axis. Since
${g}_{tt}$ is constant in (42), the string does not exert any gravitational force on test particles. However, test particles do suffer a perturbation when passing through a cosmic string. To see that, let us consider that one is in a comoving frame in which the string is at rest and, by a Lorentz transformation, the baryonic matter is approaching the string with the velocity
${v}_{c}$. From Figure 2, we can see how the velocity of the test particles are perturbed and we can calculate this perturbation. Of course, it depends on the deficit angle and, hence, on the parameters of the
$f\left(R\right)$ theory

${u}_{i}=\delta {v}_{c}\approx \frac{\delta \varphi}{2}\gamma {v}_{c}={\text{e}}^{-2{c}_{4}}{c}_{6}\pi \gamma {v}_{c},$ (44)

where $\gamma ={\left(1-{v}_{c}^{2}\right)}^{-1/2}$, considering $c=1$.

In this way, particles which move in regions where $y>0$ may collide with particles which move in regions where $y<0$ after passing through the string and form stable structures called “wakes”, see Figure 3.

4.2. The Evolution of Wakes: The Zel’dovich Approximation

Let us now make a quantitative description of the accretion problem using the Zel’dovich approximation [24], which consists in considering the Newtonian

Figure 1. The Conical Spacetime.

Figure 2. (a) A cosmic string moving with constant velocity *v _{c}*. (b) Particles moving with constant velocity

Figure 3. Particles collision with impact parameters equal to *R* and *R'*.

accretion problem in an expanding Universe by means of the method of linear approximation.

Suppose that the wake is formed at ${t}_{i}>{t}_{eq}$, where ${t}_{eq}$ is the time in which matter starts to dominate over radiation. The physical trajectory of a particle can be written as

$r\left(q\mathrm{,}t\right)=a\left(t\right)\left(q-\psi \left(q\mathrm{,}t\right)\right)\mathrm{,}$ (45)

where $q$ is the unperturbed comoving position and $\psi $ is the comoving displacement developed as a consequence of the gravitational attraction induced by the wake on the particle.

The equation of motion in the Newtonian limit is given by

$\stackrel{\xa8}{r}=-\frac{\partial \Phi \left(r,t\right)}{\partial \stackrel{\to}{r}},$ (46)

where
$\Phi \left(r,t\right)$ is the Newtonian potential which obeys the Poisson equation^{3}

${\nabla}^{2}\Phi =4\pi G\rho \left(r\mathrm{,}t\right)\mathrm{.}$ (47)

The matter density $\rho \left(r\mathrm{,}t\right)$ is determined in terms of the background density ${\rho}_{0}\left(t\right)$ as [24]

$\rho \left(r\mathrm{,}t\right){\text{d}}^{3}r={a}^{3}\left(t\right){\rho}_{0}\left(t\right){\text{d}}^{3}q$ (48)

which gives

$\rho \left(r\mathrm{,}t\right)\simeq {\rho}_{0}\left(1+\frac{\partial {\psi}_{r}\left(q\mathrm{,}t\right)}{\partial q}\right)\mathrm{,}$ (49)

where $r\mathrm{,}q$ and ${\psi}_{r}$ are the radial components of these quantities.

Replacing Equation (49) in the Poisson Equation (47), we obtain

$\frac{\partial \Phi}{\partial \stackrel{\to}{r}}=4\pi G\left[\frac{{\rho}_{0}\left(t\right)}{3}r+{\rho}_{0}\left(t\right)a\left(t\right)\psi \left(q\mathrm{,}t\right)\right]\mathrm{.}$ (50)

If we replace the relation $\stackrel{\xa8}{a}/a=-4\pi G{\rho}_{0}/3$ in the equation above, we get the linearized equation for $\psi $

$\stackrel{\xa8}{\psi}+2\frac{\stackrel{\dot{}}{a}}{a}\stackrel{\dot{}}{\psi}+3\frac{\stackrel{\xa8}{a}}{a}\psi =0\text{\hspace{0.05em}}\text{\hspace{0.05em}}.$ (51)

Since we are working in the matter era, $a\left(t\right)\propto {t}^{2/3}$. Hence, Equation (51) becomes

$\stackrel{\xa8}{\psi}+\frac{4}{3t}\stackrel{\dot{}}{\psi}-\frac{2}{3{t}^{2}}\psi =0\text{\hspace{0.05em}}\text{\hspace{0.05em}}.$ (52)

The equation above is the well-known Euler equation. Applying appropriate conditions such as $\psi \left({t}_{i}\right)=0$ and $\stackrel{\dot{}}{\psi}\left({t}_{i}\right)=-{u}_{i}$ the solution is

$\psi \left(x\mathrm{,}t\right)=\frac{3}{2}\left[\frac{{u}_{i}{t}_{i}^{2}}{t}-{u}_{i}{t}_{i}{\left(t/{t}_{i}\right)}^{2/3}\right]\mathrm{.}$ (53)

The comoving coordinate $q\left(t\right)$ can be calculated using the fact that $\stackrel{\dot{}}{r}=0$, which means that, eventually, the particle stops expanding with the Hubble flow and starts to collapse onto the wake. Therefore, we get

$q\left(t\right)=-\frac{6}{5}\left[\frac{{u}_{i}{t}_{i}^{2}}{t}-{u}_{i}{t}_{i}{\left(t/{t}_{i}\right)}^{2/3}\right]\mathrm{.}$ (54)

Hence, we are now able to compute the wake’s thickness $d\left(t\right)$ and surface density $\sigma \left(t\right)$ [25]

$d\left(t\right)\approx 2q\left(t\right){\left(\frac{t}{{t}_{i}}\right)}^{2/3}$ (55)

$\sigma \left(t\right)\approx {\rho}_{0}\left(t\right)d\left(t\right).$ (56)

Finally, we obtain

$d\left(t\right)\approx \frac{12}{5}{a}_{0}\pi \gamma {v}_{c}{\left[\frac{{t}^{4/3}}{{t}_{i}^{1/3}}-\frac{{t}_{i}^{4/3}}{{t}^{1/3}}\right]}^{2/3}$ (57)

$\sigma \approx \frac{12}{5}{\rho}_{0}{a}_{0}\pi \gamma {v}_{c}{\left[\frac{{t}^{4/3}}{{t}_{i}^{1/3}}-\frac{{t}_{i}^{4/3}}{{t}^{1/3}}\right]}^{2/3}$ (58)

where ${a}_{0}={\text{e}}^{-{c}_{4}}{c}_{6}$, which means that the wake’s thickness and density depend on the parameters of the $f\left(R\right)$ theory.

5. Final Remarks

In this work we considered cylindrically symmetric solutions in the framework of $f\left(R\right)$ theory of gravity. These solutions were obtained in vacuum regime and in the particular case where $R=0$. A cosmic string solution was of special interest and we studied the formation and evolution of wakes in this spacetime. Comparing our results with those obtained previously in the literature [25] [26], both in General Relativity and in Scalar-Tensor theories, respectively, we can see that they resemble the GR wakes instead of the scalar-tensor ones, as we would expect since ${g}_{tt}$ is constant and there is no gravitational force exerted by the $f\left(R\right)$ cosmic string in the same way as the GR cosmic string. However, for a precise comparison and further numerical evaluation, we must consider the internal cosmic string matter configuration because we need to determine the metric constants ${c}_{4}\mathrm{,}{c}_{6}$. In particular, they must obey the GUT cosmic string order of magnitude for all parameters.

As we expected, all wake’s physical quantities depended on the parameters of the particular theory of gravity under consideration. But, again, in order to make a quantitative evaluation, we need to consider not a vacuum solution but a full energy-momentum tensor for the internal cosmic string configuration in the same way as [16] [27]. This is under consideration now.

As a future perspective of this work, we plan to study the $f\left(R\right)$ cosmic string as a generator of the rotational curves in galaxies [1]. This work will come as a forthcoming paper.

Acknowledgments

P. O. Mesquita and M. E. X. Guimarães would like to thank PIBIC/CNPq (Conselho Nacional de Desenvolvimento Cientfico e Tecnológico) for a support during the preparation of this work.

NOTES

^{1}The calculations are long but straightforward.

^{2}We will suppress the^{ }symbol from now on because it will not cause any confusion.

^{3}In order to avoid confusion, we have changed the radial coordinate for *r* instead of
$\rho $ which represents now the matter density.

References

[1] Rubin, V. and Ford Jr., W.K. (1970) Rotation of the Andromeda Nebula from a Spectroscopic Survey of Emission Regions. The Astrophysical Journal, 159, 379.

https://doi.org/10.1086/150317

[2] Page, L., et al. (2007) Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Polarization Analysis. The Astrophysical Journal, 170, 335.

[3] Capozziello, S. and de Laurentis, M. (2012) The Dark Matter Problem from f(R) Gravity Viewpoint. Annalen der Physik, 524, 545-578.

https://doi.org/10.1002/andp.201200109

[4] Di Bari, P., King, S.F. and Merle, A. (2013) Dark Radiation or Warm Dark Matter from Long Lived Particle Decays in the Light of Planck. Physics Letters B, 724, 77-83.

https://doi.org/10.1016/j.physletb.2013.06.003

[5] Riess, A.G., et al. (2004) Type Ia Supernova Discoveries at z > 1 from the Hubble Space Telescope: Evidence for Past Deceleration and Constraints on Dark Energy Evolution. The Astrophysical Journal, 607, 665.

[6] Santos, J., et al. (2007) Energy Conditions in f(R) Gravity. Physical Review D, 76, Article ID: 083513.

[7] Costa, S.S., et al. (2018) Dynamical Analysis on f(R,G) Cosmology. Classical and Quantum Gravity, 35, 7.

[8] Carrol, S.M., Duvvuri, V., Trodden, M. and Turner, M.S. (2004) Is Cosmic Speed-Up Due to New Gravitational Physics? Physical Review D, 70, Article ID: 043528.

https://doi.org/10.1103/PhysRevD.70.043528

[9] Fay, S., Tavakol, R. and Tsujikawa, S. (2007) f(R) Gravity Theories in Palatini Formalism: Cosmological Dynamics and Observational Constraints. Physical Review D, 74, Article ID: 063509.

[10] Ostrogradski, M. (1850) Mémoires sur les équations différentielles, relatives au problème des isopérimètres. Mémoires de l’Académie Impériale des Sciences de St-Pétersbourg VI, 4, 385.

[11] Woodard, R.P. (2007) Ostrogradsky’s Theorem on Hamiltonian Instability. Lecture Notes in Physics, 720, 403.

[12] Corda, C. (2009) Interferometric Detection of Gravitational Waves: The Definitive Test for General Relativity. International Journal of Modern Physics D, 18, 2275-2282.

https://doi.org/10.1142/S0218271809015904

[13] Liu, H., Wang, X., Li, H. and Ma, Y. (2017) Distinguishing f(R) Theories from General Relativity by Gravitational Lensing Effect. The European Physical Journal C, 77, Article No. 723.

https://doi.org/10.1140/epjc/s10052-017-5308-5

[14] Kibble, T.W. (1976) Topology of Cosmic Domains and Strings. Journal of Physics A: Mathematical and General, 9, 1378.

https://doi.org/10.1088/0305-4470/9/8/029

[15] Vilenkin, A. and Shellard, E.P. (1994) Cosmic String and Other Topological Defects. Cambridge University Press, Cambridge.

[16] Guimarães, M.E.X. (1997) Cosmic String in Scalar-Tensor Gravities. Classical and Quantum Gravity, 14, 435.

https://doi.org/10.1088/0264-9381/14/2/016

[17] Peter, P., Guimarães, M.E.X. and de Andrade, V.C. (2003) Dilatonic Current-Carrying Cosmic Strings. Physical Review D, 67, Article ID: 123509.

https://doi.org/10.1103/PhysRevD.67.123509

[18] Oliveira, A.L.N. and Guimarães, M.E.X. (2003) Wakes in Dilatonic Current-Carrying Cosmic Strings. Physical Review D, 67, 123514.

https://doi.org/10.1103/PhysRevD.67.123514

[19] Rosa, T.O., Guimarães, M.E.X. and Neto, J.L. (2019) On an Exact Cylindrically Symmetric Solution in a Born-Infeld Type Theory of Gravity. Journal of High Energy Physics, Gravitation and Cosmology, 5, 711-718.

https://doi.org/10.4236/jhepgc.2019.53038

[20] Sotiriou, T.P. and Faraoni, V. (2010) f(R) Theories of Gravity. Reviews of Modern Physics, 82, 451.

https://doi.org/10.1103/RevModPhys.82.451

[21] Azadi, A., Momeni, D. and Nouri-Zonoz, M. (2008) Cylindrical Solutions in Metric f(R) Gravity. Physics Letters B, 670, 210.

[22] Hiscock, W.A. (1985) Exact Gravitational Field of a String. Physical Review D, 31, 3288.

https://doi.org/10.1103/PhysRevD.31.3288

[23] Levi-Civita, T. (1917) ds2 Einsteiniani in Campi Newtoniani. Rendiconti. Accademia Nazionale dei Lincei, 26, 307.

[24] Aguirre, A.N. and Brandenberger, R.H. (1995) Accretion of Hot Dark Matter onto Slowly Moving Cosmic Strings. International Journal of Modern Physics D, 4, 711-721.

https://doi.org/10.1142/S0218271895000478

[25] Vachaspati, T. (1986) Cosmic Strings and the Large-Scale Structure of the Universe. Physical Review Letters, 57, 1655.

https://doi.org/10.1103/PhysRevLett.57.1655

[26] Masalskiene, S.R.M. and Guimarães, M.E.X. (2000) Cosmic String Wakes in Scalar-Tensor Gravities. Classical and Quantum Gravity, 17, 3055.

https://doi.org/10.1088/0264-9381/17/15/313

[27] Garfinkle, D. (1985) General Relativistic Strings. Physical Review D, 32, 1323.