Observational constraints on G-corrected holographic dark energy using a Markov chain Monte Carlo method

We constrain holographic dark energy (HDE) with time varying gravitational coupling constant in the framework of the modified Friedmann equations using cosmological data from type Ia supernovae, baryon acoustic oscillations, cosmic microwave background radiation and X-ray gas mass fraction. Applying a Markov Chain Monte Carlo (MCMC) simulation, we obtain the best fit values of the model and cosmological parameters within 1σ confidence level (CL) in a flat universe as: $\varOmega_{b}h^{2}=0.0222^{+0.0018}_{-0.0013}$ , $\varOmega_{c}h^{2}=0.1121^{+0.0110}_{-0.0079}$ , $\alpha_{G}\equiv \dot{G}/(HG) =0.1647^{+0.3547}_{-0.2971}$ and the HDE constant $c=0.9322^{+0.4569}_{-0.5447}$ . Using the best fit values, the equation of state of the dark component at the present time w _d0 at 1σ CL can cross the phantom boundary w=?1.     

Holographic f(T) gravity model

We try to study the corresponding relation between f(T) gravity and holographic dark energy (HDE). A kind of energy density from f(T) is introduced which has the same role as HDE density. A f(T) model according to the HDE model is calculated. We find out a torsion scalar T based on the scalar factor is assumed by Capoziello et al. (Phys. Lett. B 639:135, 2006). The effective torsion equation of state, deceleration parameter of the holographic f(T)-gravity model are calculated.     

The generalized second law of thermodynamics for the interacting in f(T) gravity

We study the validity of the generalized second law (GSL) of gravitational thermodynamics in a non-flat FRW universe containing the interacting in f(T) gravity. We consider that the boundary of the universe to be confined by the dynamical apparent horizon in FRW universe. In general, we discuss the effective equation of state, deceleration parameter and GLS in this framewok. Also, we find that the interacting-term Q modifies these quantities and in particular, the evolution of the total entropy, results in an increases on the GLS of thermodynamic, by a factor $4\pi R_{A}^{3} Q/3$ . By using a viable f(T) gravity with an exponential dependence on the torsion, we develop a model where the interaction term is related to the total energy density of matter. Here, we find that a crossing of phantom divide line is possible for the interacting-f(T) model.     

Reconstruction of f(R), f(T) and f(\mathcal{G}) models inspired by variable deceleration parameter

We study an special law for the deceleration parameter, recently proposed by Akarsu and Dereli, in the context of f(R), f(T) and $f(\mathcal{G})$ theories of modified gravity. This law covers the law of Berman for obtaining exact cosmological models to account for the current acceleration of the universe, and also gives the opportunity to generalize many of the dark energy models having better consistency with the cosmological observations. Our aim is to reconstruct the f(R), f(T) and $f(\mathcal{G})$ models inspired by this law of variable deceleration parameter. Such models may then exhibit better consistency with the cosmological observations.     

Optical spectral decomposition of NGC 4051

Considering the host galaxy contribution, a spectral decomposition method is used to reanalyzed the archive data of optical spectra for a narrow line Seyfert 1 galaxy, NGC 4051. The light curves of the continuum f _λ (5100 Å), and Hβ, He ii, Fe ii emission lines are given. We find strong flux correlations between line emissions of Hβ, He ii, Fe ii and the continuum f _λ (5100 Å). These low-ionization lines (Hβ, Fe ii, He ii) have “inverse” intrinsic Baldwin effects. Using the methods of the cross-correlation function and the Monte Carlo simulation, we find the time delays, with respect to the continuum, are $3.45^{+12.0}_{-0.5}~\mbox{days}$ with the probability of 34 % for the intermediate component of Hβ, $6.45^{+13.0}_{-1.0}~\mbox{days}$ with the probability of 65 % for the intermediate component of He ii. From these intermediate components of Hβ and He ii, the calculated central black hole masses are $0.86^{+4.35}_{-0.33}\times 10^{6}$ and $0.82^{+3.12}_{-0.45}\times 10^{6}~M_{\odot }$ . We also find that the time delays for Fe ii are $9.7^{+3.0}_{-5.0}~\mbox{days}$ with the probability of 36 %, $8.45^{+1.0}_{-2.0}~\mbox{days}$ with the probability of 18 % for the total epochs and “subset 1” data, respectively. It seems that the Fe ii emission region is outside of the Hβ emission region.     

A Note on Computation of Relative Magnetic-Helicity Flux Across the Photosphere

A number of investigations of the rate of relative magnetic-helicity transport across the photosphere [ ${{{\rm d} H}/{{\rm d}t}}|_{S}$ ] have reported differences in the estimates computed from two different formulations of the relative-helicity flux-density proxy G _A and G _θ . There have been suggestions that G _θ is a more robust helicity-flux density proxy and that the differences in the estimates of ${{{\rm d} H}/{{\rm d}t}}|_{S}$ are caused by biases in G _A, noise, and/or the boundary conditions. In this note, we prove that the differences are caused by the inconsistent choice of boundary conditions in the explicit or implicit Green’s function [ $\mathcal{G}(\mathbf {x},\mathbf{x}')$ ] used for computing G _A and G _θ when comparing the helicity-flux estimates based on G _A and G _θ . When the boundary conditions in $\mathcal{G}$ are chosen consistently, the two helicity-flux density proxies, [G _A and G _θ ] produce essentially identical results for the rate of helicity transport across the photosphere. They also yield essentially identical results for the rate of helicity transport of the shearing and advection terms separately. Using MHD simulation, HMI observational data, and Monte Carlo simulations of noise we show that this result is robust. Neither the shape of the active region, nor the shape of the boundary, nor data noise causes any difference in the rate of helicity transport computed via G _A and G _θ .     

Ultra-high energy neutrino fluxes from supermassive AGN black holes

We compute the ultra-high energy (UHE) neutrino fluxes from plausible accreting supermassive black holes closely linking to the 377 active galactic nuclei (AGNs). They have well-determined black hole masses collected from the literature. The neutrinos are produced via simple or modified URCA processes, even after the neutrino trapping, in superdense proto-matter medium. The resulting fluxes are ranging from: (1) (quark reactions)— $J^{q}_{\nu\varepsilon}/(\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1})\simeq8.29\times 10^{-16}$ to 3.18×10^?4, with the average $\overline{J}^{q}_{\nu\varepsilon}\simeq5.53\times 10^{-10}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ , where ε _d ～10^?12 is the opening parameter; (2) (pionic reactions)— $J^{\pi}_{\nu\varepsilon} \simeq0.112J^{q}_{\nu\varepsilon}$ , with the average $J^{\pi}_{\nu\varepsilon} \simeq3.66\times 10^{-11}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ ; and (3) (modified URCA processes)— $J^{URCA}_{\nu\varepsilon}\simeq7.39\times10^{-11} J^{q}_{\nu\varepsilon}$ , with the average $\overline{J}^{URCA}_{\nu\varepsilon} \simeq2.41\times10^{-20} \varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ . We conclude that the AGNs are favored as promising pure neutrino sources, because the computed neutrino fluxes are highly beamed along the plane of accretion disk, peaked at high energies and collimated in smaller opening angle θ～ε _d .     

Fourier analysis of the light curves of eclipsing variables,XXIII

The method of evaluating the photometric perturbationsB _2m of eclipsing variables in the frequency domain, developed by Kopal (1959, 1975e, 1978) for an interpretation of mutual eclipses in systems whose components are distorted by axial rotation and mutual tidal action. The aim of the present paper has been to establish explicit expressions for the photometric perturbationB _2m in such systems, regardless of the kind of eclipses and non-integral values ofm. Recently, Kopal (1978) introduced two different kinds of integrals with respect to associated α-functions andI-integrals which have been expressed in terms of certain general types of series that can be easily programmed for automatic computation within seconds of real time on highspeed computers. Following a brief introduction (Section 1) in which the need of this new approach will be expounded, in Section 3 we shall deduce the integral $$\int_0^{\theta \prime } {\tfrac{{\alpha _n^\prime }}{\delta }} d(sin^{2m} \theta )$$ in terms of a certain general type of series and also β-function, which should enable us to evaluate explicit expressions forf _* ^(h) ,f ₁ ^(h) ,f ₂ ^(h) as well asB _2m .     

The seasonal variation of the newtonian constant of gravitation and its relationship with the “classical tests” of general relativity

The ratio between the Earth's perihelion advance (Δθ) _E and the solar gravitational red shift (GRS) (Δø _s e)a ₀/c ² has been rewritten using the assumption that the Newtonian constant of gravitationG varies seasonally and is given by the relationship, first found by Gasanalizade (1992b) for an aphelion-perihelion difference of (ΔG)_a?p . It is concluded that $$\begin{gathered} (\Delta \theta )_E = \frac{{3\pi }}{e}\frac{{(\Delta \phi _{sE} )_{A_0 } }}{{c^2 }}\frac{{(\Delta G)_{a - p} }}{{G_0 }} = 0.038388 \sec {\text{onds}} {\text{of}} {\text{arc}} {\text{per}} {\text{revolution,}} \hfill \\ \frac{{(\Delta G)_{a - p} }}{{G_0 }} = \frac{e}{{3\pi }}\frac{{(\Delta \theta )_E }}{{(\Delta \phi _{sE} )_{A_0 } /c^2 }} = 1.56116 \times 10^{ - 4} . \hfill \\ \end{gathered} $$ The results obtained here can be readily understood by using the Parametrized Post-Newtonian (PPN) formalism, which predicts an anisotropy in the “locally measured” value ofG, and without conflicting with the general relativity.     

Stability of ion-acoustic waves in a pair-ion plasma with a third species of ions: application to cometary plasmas

A popular model of a cometary plasma is hydrogen (H⁺) with positively charged oxygen (O⁺) as a heavier ion component. However, the discovery of negatively charged oxygen (O^?) ions enables one to model a cometary plasma as a pair-ion plasma (of O⁺ and O^?) with hydrogen as a third ion constituent. We have, therefore, studied the stability of the ion-acoustic wave in such a pair-ion plasma with hydrogen and electrons streaming with velocities $V_{d\mathrm{H}^{+}}$ and V _de , respectively, relative to the oxygen ions. We find the calculated frequency of the ion-acoustic wave with this model to be in good agreement with the observed frequencies. The ion-acoustic wave can also be driven unstable by the streaming velocity of the hydrogen ions. The growth rate increases with increasing hydrogen density $n_{\mathrm{H}^{+}}$ , and streaming velocities $V_{d\mathrm{H}^{+}}$ and V _de . It, however, decreases with increasing oxygen ion densities $n_{\mathrm{O}^{+}}$ and $n_{\mathrm{O}^{-}}$ .     

The Chandra X-ray galaxy clusters at z<1.4: constraints on the evolution of L X −T−M g relations

We analyzed the luminosity-temperature-mass of gas (L _X ?T?M _g ) relations for a sample of 21 Chandra galaxy clusters. We used the standard approach (β?model) to evaluate these relations for our sample that differs from other catalogues since it considers galaxy clusters at higher redshifts (0.4<z<1.4). We assumed power-law relations in the form $L_{X} \sim(1 +z)^{A_{L_{X}T}} T^{\beta_{L_{X}T}}$ , $M_{g} \sim(1 + z)^{A_{M_{g}T}} T^{\beta_{M_{g}T}}$ , and $M_{g} \sim(1 + z)^{A_{M_{g}L_{X}}} L^{\beta_{M_{g}L_{X}}}$ . We obtained the following fitting parameters with 68 % confidence level: $A_{L_{X}T} = 1.50 \pm0.23$ , $\beta_{L_{X}T} = 2.55 \pm0.07$ ; $A_{M_{g}T} = -0.58 \pm0.13$ and $\beta_{M_{g}T} = 1.77 \pm0.16$ ; $A_{M_{g}L_{X}} \approx-1.86 \pm0.34$ and $\beta_{M_{g}L_{X}} = 0.73 \pm0.15$ , respectively. We found that the evolution of the M _g ?T relation is small, while the M _g ?L _X relation is strong for the cosmological parameters Ω _m =0.27 and Ω _Λ =0.73. In overall, the clusters at high-z have stronger dependencies between L _X ?T?M _g correlations, than those for clusters at low-z. For most of galaxy clusters (first of all, from MACS and RCS surveys) these results are obtained for the first time.     

A relation between Giorgilli-Galgani's and Deprit's algorithms

If \(T = \sum\nolimits_{i = 1}^\infty {\varepsilon ^i } T_i\) and \(W = \sum\nolimits_{n = 1}^\infty {n\varepsilon ^{n - 1} } W^{\left( n \right)}\) are respectively the generators of Giorgilli-Galgani's and Deprit's transformations, we show that the change of variables generated byT is the inverse of the one generated byW, ifT _i =W ⁽ⁱ⁾ for anyi. The method used is to show that the recurrence which defines the first algorithm can also be obtained with the second one.     

A comparison of theoretical Fe xii emission line strengths with EUV observations of a solar active region

New theoretical electron-density-sensitive Fe xii emission line ratios $$R_1 = I(3s^2 3p^3 {}^4S_{3/2} - 3s3p^4 {}^4P_{5/2} )/I(3s^2 3p^3 {}^2P_{3/2} - 3s3p^4 D_{5/2} )$$ and $$R_2 = I(3s^2 3p^3 {}^2P_{3/2} - 3s3p^4 {}^2D_{5/2} )/I(3s^2 3p^3 {}^4S_{3/2} - 3s3p^2 P_{3/2} )$$ are derived using R-matrix electron impact excitation rate calculations. We have identified the Fexii \(3s^2 3p^3 {}^4S_{3/2} - 3s3p^4 {}^4P_{5/2} ,{\text{ }}3s^2 3p^3 {}^2P_{3/2} - 3s^3 3p^4 {}^2D_{5/2} ,{\text{ }}3s^2 3p^3 S_{3/2} - 3s^2 3p^3 P_{3/2} \) and \(3s^2 3p^3 {}^4S_{3/2} - 3s^2 3p^3 {}^2P_{1/2}\) transitions in an active region spectrum obtained with the Harvard S-055 spectrometer on board Skylab at wavelengths of 364.0, 382.8, 1241.7, and 1349.4 Å, respectively. Electron densities determined from the observed values of R ₁ (log N _e ? 11.0) and R ₂(log N _e ? 11.4) are significantly larger than the typical active region measurements, but are similar to those derived from some active region spectra observed with the Skylab 2082A instrument, which provides observational support for the atomic data adopted in the line ratio calculations, and also for the identification of the Fe xii transitions in the S-055 spectrum. However the observed value of R ₃ = I(1349.4 Å)/I(1241.7 Å) is approximately a factor of two larger than one would expect from theory which, considering that the 1349.4 Å line lies at the edge of the S-055 wavelength coverage, may reflect errors in the instrument efficiency curve. Another possibility is that the 1349.4 Å transition is blended, probably with Si ii 1350.1 Å.     

Sur les solutions Lagrangiennes du problème des trois corps solides avec la loi de Weber

In this work we consider the problem of translational-rotational motion of three solid bodies, for which the elementary particles attract each other according to different Weber's laws for each pair of bodies. This problem represents a special case of the generalized problem of three solids considered in a previous work, (Dubochin, 1974) and it gives an example of the verification of the existence conditions for the Lagrangian solutions. In these solutions, the centers of mass always for m an equilateral triangle. Each body has axial symmetry with the plane of symmetry perpendicular to the axis of symmetry rotates uniformly around this axis, which at any instant stays perpendicular to the plane of the triangle formed by the centers of mass. According to Weber's law (Tisserand, 1896) the elementary particles of two bodiesT _i andT _j (i, j=0, 1, 2) are attracted by forces which are proportional to the function $$F_{ij} (W) = \frac{{f_{ij} }}{{\Delta _{ij^2 } }}\left[ {1 - a_{ij} \dot \Delta _{ij^2 } + 2a_{ij} \Delta _{ij} \ddot \Delta _{ij} } \right]$$ wheref _ij anda _ij (in generalf _ji ≠f _ij anda _ji ≠a _ij ) are functions of the timet, and where the real quantities Δ_ij are the mutual distances between the particles of the bodiesT _i andT _j , and where \(\dot \Delta _{ij} \) and \(\ddot \Delta _{ij} \) are their derivatives with respect to the time. The analysis of the general conditions for the Lagrangian solutions gives the following results for the case of Weber's laws.

1. Only the invariant Lagrangian solutions, (the traingle of the centres of mass does not change in time) are possible in this problem.
2. Besides the conditions (NL) obtained in the case of the Newton-Coulomb law, (all thea _ij are zero), the complementary conditions (WL) must be satisfied.

In particular, if all the bodies are spheres or homogeneous ellipsoids, they must necessarily have the same dimensions, but they can have different masses.     

The planetary spin and rotation period: a modern approach

Using a new approach, we have obtained a formula for calculating the rotation period and radius of planets. In the ordinary gravitomagnetism the gravitational spin (S) orbit (L) coupling, $\vec{L}\cdot\vec{S}\propto L^{2}$ , while our model predicts that $\vec{L}\cdot\vec{S}\propto\frac{m}{M}L^{2}$ , where M and m are the central and orbiting masses, respectively. Hence, planets during their evolution exchange L and S until they reach a final stability at which MS∝mL, or $S\propto\frac{m^{2}}{v}$ , where v is the orbital velocity of the planet. Rotational properties of our planetary system and exoplanets are in agreement with our predictions. The radius (R) and rotational period (D) of tidally locked planet at a distance a from its star, are related by, $D^{2}\propto\sqrt{\frac{M}{m^{3}}}R^{3}$ and that $R\propto\sqrt{\frac {m}{M}}a$ .     

Nonlinear newtonian theory and the effect of time dilatation

The fact that the energy density ρ_g of a static spherically symmetric gravitational field acts as a source of gravity, gives us a harmonic function \(f\left( \varphi \right) = e^{\varphi /c^2 } \) , which is determined by the nonlinear differential equation $$\nabla ^2 \varphi = 4\pi k\rho _g = - \frac{1}{{c^2 }}\left( {\nabla \varphi } \right)^2 $$ Furthermore, we formulate the infinitesimal time-interval between a couple of events measured by two different inertial observers, one in a position with potential φ-i.e., dt _φ and the other in a position with potential φ=0-i.e., dt ₀, as $${\text{d}}t_\varphi = f{\text{d}}t_0 .$$ When the principle of equivalence is satisfied, we obtain the well-known effect of time dilatation.     

Sur un modele explicitant les differents types morphologiques des galaxies

In the now classical Lindblad-Lin density-wave theory, the linearization of the collisionless Boltzmann equation is made by assuming the potential functionU expressed in the formU=U ₀ + \(\tilde U\) +... WhereU ₀ is the background axisymmetric potential and \(\tilde U<< U_0 \) . Then the corresponding density distribution is \(\rho = \rho _0 + \tilde \rho (\tilde \rho<< \rho _0 )\) and the linearized equation connecting \(\tilde U\) and the component \(\tilde f\) of the distribution function is given by $$\frac{{\partial \tilde f}}{{\partial t}} + \upsilon \frac{{\partial \tilde f}}{{\partial x}} - \frac{{\partial U_0 }}{{\partial x}} \cdot \frac{{\partial \tilde f}}{{\partial \upsilon }} = \frac{{\partial \tilde U}}{{\partial x}}\frac{{\partial f_0 }}{{\partial \upsilon }}.$$ One looks for spiral self-consistent solutions which also satisfy Poisson's equation $$\nabla ^2 \tilde U = 4\pi G\tilde \rho = 4\pi G\int {\tilde f d\upsilon .} $$ Lin and Shu (1964) have shown that such solutions exist in special cases. In the present work, we adopt anopposite proceeding. Poisson's equation contains two unknown quantities \(\tilde U\) and \(\tilde \rho \) . It could be completelysolved if a second independent equation connecting \(\tilde U\) and \(\tilde \rho \) was known. Such an equation is hopelesslyobtained by direct observational means; the only way is to postulate it in a mathematical form. In a previouswork, Louise (1981) has shown that Poisson's equation accounted for distances of planets in the solar system(following to the Titius-Bode's law revised by Balsano and Hughes (1979)) if the following relation wasassumed $$\rho ^2 = k\frac{{\tilde U}}{{r^2 }} (k = cte).$$ We now postulate again this relation in order to solve Poisson's equation. Then, $$\nabla ^2 \tilde U - \frac{{\alpha ^2 }}{{r^2 }}\tilde U = 0, (\alpha ^2 = 4\pi Gk).$$ The solution is found in a classical way to be of the form $$\tilde U = cte J_v (pr)e^{ - pz} e^{jn\theta } $$ wheren = integer,p =cte andJ _v (pr) = Bessel function with indexv (v ² =n ² + α²). By use of the Hankel function instead ofJ _v (pr) for large values ofr, the spiral structure is found to be given by $$\tilde U = cte e^{ - pz} e^{j[\Phi _v (r) + n\theta ]} , \Phi _v (r) = pr - \pi /2(v + \tfrac{1}{2}).$$ For small values ofr, \(\tilde U\) = 0: the center of a galaxy is not affected by the density wave which is onlyresponsible of the spiral structure. For various values ofp,n andv, other forms of galaxies can be taken into account: Ring, barred and spiral-barred shapes etc. In order to generalize previous calculations, we further postulateρ ₀ =kU ₀/r ², leading to Poisson'sequation which accounts for the disc population $$\nabla ^2 U_0 - \frac{{\alpha ^2 }}{{r^2 }}U_0 = 0.$$ AsU ₀ is assumed axisymmetrical, the obvious solution is of the form $$U_0 = \frac{{cte}}{{r^v }}e^{ - pz} , \rho _0 = \frac{{cte}}{{r^{2 + v} }}e^{ - pz} .$$ Finally, Poisson's equation is completely solvable under the assumptionρ =k(U/r ². The general solution,valid for both disc and spiral arm populations, becomes $$U = cte e^{ - pz} \left\{ {r^{ - v} + } \right.\left. {cte e^{j[\Phi _v (r) + n\theta ]} } \right\},$$ The density distribution along the O _z axis is supported by Burstein's (1979) observations.     

Stability of the triangular Lagrange points beyond Gascheau��s value

We examine the stability of the triangular Lagrange points L ₄ and L ₅ for secondary masses larger than the Gascheau??s value ${\mu_{\rm G}= (1-\sqrt{23/27}/2)= 0.0385208\ldots}$ (also known as the Routh value) in the restricted, planar circular three-body problem. Above that limit the triangular Lagrange points are linearly unstable. Here we show that between??? _G and ${\mu \approx 0.039}$ , the L ₄ and L ₅ points are globally stable in the sense that a particle released at those points at zero velocity (in the corotating frame) remains in the vicinity of those points for an indefinite time. We also show that there exists a family of stable periodic orbits surrounding L ₄ or L ₅ for ${\mu \ge \mu_G}$ . We show that??? _G is actually the first value of a series ${\mu_0 (=\mu_G), \mu_1,\ldots, \mu_i,\ldots}$ corresponding to successive period doublings of the orbits, which exhibit ${1, 2, \ldots, 2^i,\ldots}$ cycles around L ₄ or L ₅. Those orbits follow a Feigenbaum cascade leading to disappearance into chaos at a value ${\mu_\infty = 0.0463004\ldots}$ which generalizes Gascheau??s work.     

Some identities for spheroidal harmonics

The spheroidal harmonics expressions $$\left[ {P_{2k}^{2s} \left( {i\xi } \right)P_{2k - 2r}^{2s} \left( \eta \right) - P_{2k - 2r}^{2s} \left( {i\xi } \right)P_{2k}^{2s} \left( \eta \right)} \right]e^{i2s\theta } $$ and $$\left[ {\eta ^2 P_{2k}^{2s} \left( {i\xi } \right)P_{2k - 2r}^{2s} \left( \eta \right) + \xi ^2 P_{2k - 2r}^{2s} \left( {i\xi } \right)P_{2k}^{2s} \left( \eta \right)} \right]e^{i2s\theta } $$ , have ξ²+η² as a factor. A method is presented for obtaining for these two expressions the coefficient of ξ²+η² in the form of a linear combination of terms of the formP _2m ^2s (iξ)P _2n ^2s (η)e ^i2sθ. Explicit formulae are exhibited for the casesr=1, 2, 3 and any positive or zero integersk ands. Such identities are useful in gravitational potential theory for ellipsoidal distributions when matching Legendre function expansions are employed.     

The unification of the inflation with late-time acceleration in Born-Infeld-f(R) gravity

We study accelerating dynamics from Born-Infeld-f(R) gravity in a simplified conformal approach without matter. In Makarenko et al. (arXiv:1404.2850 [gr-qc], 2011b) it was derived eventually any Dark Energy cosmology from above theory. In this Letter we apply the technique of Makarenko et al. (arXiv:1404.2850 [gr-qc], 2011b) to show that Born-Infeld-f(R) gravity may describe very realistic universe admitting the unification of early-time inflation with late-time acceleration. Specifically, the evolution with periodic as well as non-periodic behavior is considered with possibility to cross the phantom-divide at early or late-times.