Research Article | Open Access
M. Farasat Shamir, "Anisotropic Universe in Gravity", Advances in High Energy Physics, vol. 2017, Article ID 6378904, 10 pages, 2017. https://doi.org/10.1155/2017/6378904
Anisotropic Universe in Gravity
This paper is devoted to investigating the recently introduced theory of gravity, where is the Gauss-Bonnet term and is the trace of the energy-momentum tensor. For this purpose, anisotropic background is chosen and a power law gravity model is used to find the exact solutions of field equations. In particular, a general solution is obtained which is further used to reconstruct some important solutions in cosmological contexts. The physical quantities like energy density, pressure, and equation of state parameter are calculated. A Starobinsky-like model is proposed which is used to analyze the behavior of universe for different values of equation of state parameter. It is concluded that presence of term in the bivariate function may give many cosmologically important solutions of the field equations.
Some important modifications of general relativity (GR) have been proposed in the last two decades. The mostly discussed theories are and theories of gravity ( is the Ricci scalar and is the trace of energy-momentum tensor) which have been treated most seriously [1–14]. However, recently cosmology has been severely challenged . On the other hand, modified Gauss-Bonnet (GB) gravity is another theory which has gained popularity in the last few years [16–18]. It is also known as theory of gravity, where is a generic function of GB invariant . GB term plays an important role as it may avoid ghost contributions and helps in regularizing the gravitational action . Thus, to save theories, one may include terms. A theory with a similar idea has been recently proposed named as gravity . Some interesting work has been done in the recent past using modified GB theories.
Anisotropic compact stars in modified gravity have been discussed by Abbas et al. . Houndjo et al.  found the exact solutions of field equations using cylindrical symmetry and it was concluded that there existed seven families of exact solutions for three different forms of models. The exact cylindrically symmetric solutions of modified field equations recovered cosmic string space-time . A more generalized version of GB gravity known as gravity has also been discussed widely. Wu and Ma  investigated the spherically symmetric solutions at low energy where the weak-field and slow-motion limit of gravity was developed. Laurentis et al.  discussed cosmological inflation in theory. Sharif and Ikram  examined warm inflation in theory of gravity using scalar fields for the Friedmann-Robertson-Walker (FRW) universe model. Conserved quantities have been recently explored in FRW background using Noether symmetry . Sharif and Fatima  discussed the role of GB term for the early and late-time accelerating phases of the universe by considering two viable models. García et al.  explored energy conditions to prove the viability of some gravity models. gravity energy conditions have been recently explored where the WEC was used along with the recent estimated values of cosmological parameters to determine the viability of some specific choices of gravity models .
In the context of gravity, Sharif and Ikram proved that the massive test particles followed nongeodesic lines of geometry due to the presence of extra force and examined the energy conditions for flat FRW universe . The same authors  used reconstruction techniques to reproduce the cosmic evolution corresponding to de Sitter universe, power law solutions, and phantom/nonphantom eras in this theory. In a recent paper , theory of gravity was discussed using Noether symmetry approach. Two specific models were studied to determine the conserved quantities and it was concluded that the well-known de Sitter solution could be reconstructed for some specific choice of gravity model. In a recent paper, Sharif and Ikram  studied the wormhole solutions using power law gravity models and it was shown that traversable wormhole solutions were physically acceptable in theory of gravity. In another work , Noether symmetry methodology has been used to study some cosmologically important gravity models with anisotropic background. It is concluded that the specific models of modified GB gravity may be used to reconstruct CDM cosmology without involving any cosmological constant. For some particular choices of gravity models, it is anticipated that this theory may explain the late-time cosmic acceleration. Thus it seems interesting to explore further the modified gravity. Moreover, in comparison with gravity, the presence of the matter term in the bivariate function may give many solutions of the field equations and the theory may support the accelerated expansion of universe under certain conditions for the model under consideration.
In this paper, we are focussed to investigate the dynamics of gravity with anisotropic background. It is well-known that the isotropic models are among the best choices to study large scale structure of the universe. Moreover, according to the cosmological observations including the cosmic microwave background (CMB) radiation, the current universe is isotropic. However, it is believed that the early universe may not have been exactly uniform. Also the local anisotropies that we observe today in galaxies and super clusters also motivate us to model the universe with anisotropic background. Bianchi type models are among the simplest models with anisotropic background. In particular, the investigation of Bianchi type universe in context of modified theories is interesting. In this work, we are interested to explore gravity using locally rotationally symmetric (LRS) Bianchi type space-time. We find the exact solutions of the LRS Bianchi type field equations in theory of gravity. In particular, a general solution with power law gravity model is reported. The plan of the paper is as follows: some basics of gravity and field equations are discussed in Section 2. Section 3 is devoted to exploring the exact solutions of modified field equations. Section 4 is used to reconstruct some important cosmological solutions. Final remarks are given in the last section.
2. Some Basics of Gravity
The modified GB gravity is given by the action 
Here and denote the GB term and the trace of the energy-momentum tensor, respectively, whereas is the standard matter Lagrangian, is the Ricci scalar, is the determinant of metric tensor, and is a coupling constant. The field equations can be obtained by varying the action equation (1) with respect to the metric tensor where the symbols involved have their usual meanings and and the subscript or in the functions denotes the partial derivatives. It would be interesting to notice that if we substitute in (2), then the field equations of gravity are recovered. Moreover, the case reduces the modified field equations to the usual GR equations. For the sake of simplicity, from now onwards we consider , , and so forth. The trace of (2) gives
It may be noticed that this relates , , and differentially and not algebraically as in GR, where . This indicates that the modified field equations may admit many solutions in addition to other modified theories and GR. The covariant divergence of (2) is given bywhich is not zero. It is due to the presence of higher order derivatives of the energy-momentum tensor that are naturally present in the field equations. Thus the theory might be plagued by divergences at astrophysical scales. This seems to be an issue with some other higher order derivatives theories as well that include higher order terms of energy-momentum tensor. However, to deal with the issue, one can put some constraints to (4) to obtain standard conservation equation . Here we take the spatially homogeneous, anisotropic, LRS Bianchi type space-timewhere and are cosmic scale factors. Moreover, we consider that the universe is composed of perfect fluidwhere and denote the energy density and pressure of the fluid, respectively. The Ricci scalar and GB invariant for (5) are given aswhere the overdot denotes the derivative with respect to the time coordinate. Now we define some textbook physical quantities. The average scale factor and average Hubble parameter for the model under consideration take the form
These are three highly nonlinear and difficult differential equations with five unknowns. Thus we need an additional constraint to investigate any exact solution. Here we may consider a physical condition that shear scalar is proportional to expansion scalar which provideswhere is an arbitrary real number. In literature [35–38], many authors explored the exact solutions of field equations using this condition. Thus, using (12), field equations (11) take the form
Now we investigate the exact solutions of these field equations.
3. Exact Solutions of Modified Field Equations
We consider the model aswhere and are arbitrary constants. Further, we choose in power law form; that is,
This model has already been proposed by Cognola et al.  and it is interesting because the chances of the appearing of Big-Rip singularity vanish using this model. Subtraction of (14) and (15) yields
For simplicity and without loss of any generality, we choose so that (18) takes the form
After inserting the value of GB invariant (20) reduces to differential equations with three unknowns , , and . It would be worthwhile to mention here that many solutions can be found using (20). Here we consider the power law form; that is,where and are arbitrary constants. Using this in (20), we obtain a constraint equation
This equation is satisfied for such that
Thus, corresponding to two roots of this equation, we obtain two choices of models:where and are integration constants. The first model recovers the usual GB gravity for and . The second model with square root term is important as it leads to a viable inflation in the presence of massive scalar field . It is mentioned that different forms of can be assumed to reconstruct the solutions. However, we propose only two models for the present analysis.
The behaviour of energy density plus pressure of universe can be seen from Figure 1(a) for the model with . It is clear that as time grows which further suggests EoS parameter . This is interesting as the phantom-like dark energy is found to be in the region where . The universe with phantom dark energy ends up with a finite time future singularity known as cosmic doomsday or Big Rip [40, 41]. Moreover, accelerated expansion of universe is described when [42–44]. Figure 1(b) depicts the behavior of for radiation universe.
Using EoS parameter , it follows that
Using (28), we may analyze the behavior of universe by choosing different values of EoS parameter.(i) (dust fluid). In this case, the energy density of the universe turns out to be where The graphical behavior of energy density for both the roots is shown in Figures 2(a) and 2(b). The first root gives positive energy density but it approaches zero as grows. However, the second root provides unrealistic behavior as energy density approaches at later times.(ii) (subrelativistic fluid). Here the energy density of the universe turns out to be In this case, the first root, Figure 3(a), gives positive energy density having similar limiting value as in the dust case. However, the second root, Figure 3(b), gives negative energy density which corresponds to negative pressure universe. It is interesting because it is believed that dark energy has negative pressure causing accelerated expansion of universe.(iii) (radiation fluid). In this case, (28) becomes linear and the energy density of the universe turns out to be The graphical behavior of energy density in this case is shown in Figure 1(b).(iv) (ultrarelativistic fluid). Here the energy density of the universe turns out to be Here both the roots give positive energy density. However, the first root shows a decreasing behavior, Figure 4(a), while the second root gives increasing trend as the time passes, Figure 4(b).(v) (stiff fluid). Here the energy density of the universe turns out to be The graphical behavior in this case is similar to that of ultrarelativistic fluid universe, Figures 5(a) and 5(b).
Thus, the solution metric takes the form
Thus many physical solutions with different values of EoS parameter are possible. This is due to the fact that the presence of term in the bivariate function may give many constraints which ultimately yield different solutions of the field equations and the theory may support the accelerated expansion of universe under certain conditions for the model under consideration.
3.1. Discussion of Some Important Cosmological Parameters
The Ricci scalar and GB invariant for the solution metric (35) take the form
It is obvious that the singularities exist at and . The average Hubble parameter, average scale factor, and volume scale factor of universe take the form
These parameters attain either zero value or tend to infinity at the points of singularity. In particular, the average scale factor is zero at the initial epoch and hence the model has a point type singularity . The redshift for a distant source is directly related to the scale factor of the universe at the time when the photons were emitted from the source. The scale factor and redshift are related through the equationwhere is the present values of the scale factor. Using (37), we getwhere represents the present values of Hubble’s parameter. We can write the value of Hubble’s parameter in terms of redshift parameter as
The deceleration, jerk, and snap parameters for the solution are given by
The expansion scalar and shear scalar become
It is mentioned here that the isotropy condition as is also satisfied in this case. It is also evident that for small values of even when is not very large. This indicates that transition to isotropy is also possible for some suitable values of and (other than ). Thus, anisotropic universe isotropizes in general and FRW limit does exist with power law model under consideration. It can be seen from (37) and (42) that the volume of universe is zero at while the expansion scalar is infinite, which suggests that the universe started its evolution with zero volume at , that is, big bang scenario.
4. Reconstruction of Some Important Cosmological Solutions
Here we discuss some special cases to reconstruct some important cosmological solutions.
4.1. Flat FRW Solution
For a special case when and , space-time (35) takes the formwhich is the solution of well-known flat FRW metric. In this case attains an imaginary value which shows that the model is not viable for flat FRW solution. However, when , using (25), EoS parameter turns out to be
One can further analyze using different possibilities and hence in this case the reconstruction flat FRW solution is justified.
4.2. Kasner Type Solution
We can recover interesting Kasner type solution by putting in (35)
It is similar to the well-known Kasner’s metric  and one can obtain the exact space-time after redefining the parameters. Here EoS parameter takes the form
Hence our both values of incorporate Kasner type solution.
4.3. Exponential Law Solutions
It is to be noticed that (20) has an exponential solution of the formwith the constraint equationwhere and are arbitrary constants. Here EoS parameter becomeswhich is independent of parameter and hence any viable power law gravity model can be used with an appropriate model. The solution metric in this case turns out to be
The average Hubble parameter is zero for this solution. All other dynamical parameters expansion scalar , shear scalar , and volume scale factor of universe are constant here.
5. Final Remarks
The main purpose of this work is to investigate newly introduced modified GB theory, namely, gravity. For this purpose, anisotropic background is chosen. LRS Bianchi type space-time is the simplest model and has been used frequently in different contexts. To our knowledge, this is the first attempt to investigate the exact solutions of gravity for LRS Bianchi type space-time. Moreover, we have not used any conventional assumption like constant deceleration parameter to investigate the exact solutions. It is mentioned here that the field equations are highly nonlinear and it is due to the inclusion of bivariate function in the standard action. So we assume that the shear scalar is proportional to the expansion scalar which gives , where , are the metric coefficients and is an arbitrary constant. We consider the model , where and are arbitrary constants. Moreover, we have used a power law form gravity model already available in literature . The interesting feature of this model is that the chances of the appearing of Big-Rip singularity are minimized. Further, the viability of this model has already been discussed in different cosmological contexts [47–49]. We have also proposed a Starobinsky-like model; that is, , where is an arbitrary constant. The interesting aspect of the Starobinsky-like model is that it provides many exact solutions and one can choose the best fit solutions according to the requirements.
Using the modified field equations, we have formulated a general differential equation (20) which can further be used to investigate exact solutions. Mainly we used power law and exponential forms of metric coefficients to explore the exact solutions of modified field equations. It is shown that two gravity models are associated with the power law solution. The first model recovers the usual GB gravity for and . However, the second model with square root term looks important as it leads to a viable inflation in the presence of massive scalar field . It is mentioned here that different forms of can be assumed to reconstruct the solutions. However, we have discussed only two models for the present analysis. Firstly, we choose linear gravity model. The expressions for energy density and pressure of universe are defined for anisotropy parameters and with . The graphical behavior (see Figure 1(a)) shows that as time grows which implies that EoS parameter . This is interesting as the phantom-like dark energy is found to be in the region where . The universe with phantom dark energy ends up with a finite time future singularity known as cosmic doomsday or Big Rip [40, 41]. Moreover, accelerated expansion of universe with de Sitter type evolution is described when [42–44]. However, the second model proposed as Starobinsky-like model involves squared term of . This gives a quadratic equation in and we have analyzed the behavior of universe by choosing different values of EoS parameter. In particular, analysis is given for (dust fluid), (subrelativistic fluid), (radiation fluid), (ultrarelativistic fluid), and (stiff fluid). The graphical behavior of energy density for cases is discussed. In particular, in case of subrelativistic fluid, one root corresponds to negative pressure universe. It is interesting because it is believed that dark energy has negative pressure causing accelerated expansion of universe. Thus, it is concluded that presence of term in the bivariate function may give many solutions of the field equations and the theory supports the accelerated expansion of universe under certain conditions for the model under consideration.
Lastly, we have reconstructed some important cosmological solutions. The first solution corresponds to flat FRW space-time. This solution is valid for describing the stiff fluid universe for the first model. However, the solution is not physical for the second model as is imaginary. The second solution provides the well-known Kasner universe and it gives the value of anisotropy parameter . The third solution is obtained by exponential law assumption. It gives the average Hubble parameter zero and all other dynamical parameters like expansion scalar , shear scalar , and volume scale factor of universe constant. EoS parameter is independent of the parameter for this solution and we get an explicit expression for involving . Thus any power law gravity model may be used with an appropriate model to reconstruct some important cosmological solutions.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
- A. A. Starobinsky, “A new type of isotropic cosmological models without singularity,” Physics Letters B, vol. 91, no. 1, pp. 99–102, 1980.
- A. A. Starobinsky, “Disappearing cosmological constant in f(R) gravity,” Jounal of Experimental and Theoretical Physics, vol. 86, pp. 157–163, 2007.
- S. A. Appleby and R. A. Battye, “Do consistent models mimic general relativity plus Λ?” Physics Letters B, vol. 654, no. 1-2, pp. 7–12, 2007.
- S. A. Appleby, R. A. Battye, and A. A. Starobinsky, “Curing singularities in cosmological evolution of gravity,” Journal of Cosmology and Astroparticle Physics, vol. 1006, 2010.
- S. Nojiri and S. D. Odintsov, “Final state and thermodynamics of a dark energy universe,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 70, no. 10, Article ID 103522, 2004.
- T. Harko, F. S. N. Lobo, S. Nojiri, and S. D. Odintsov, “ gravity,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 84, no. 2, Article ID 024020, 2011.
- S. Nojiri and S. D. Odintsov, “Introduction to modified gravity and gravitational alternative for dark energy,” International Journal of Geometric Methods in Modern Physics, vol. 4, no. 1, pp. 115–145, 2007.
- K. Bamba, S. Capozziello, S. Nojiri, and S. D. Odintsov, “Dark energy cosmology: the equivalent description via different theoretical models and cosmography tests,” Astrophysics and Space Science, vol. 342, no. 1, pp. 155–228, 2012.
- M. Sharif and M. F. Shamir, “Exact solutions of Bianchi-type i and v spacetimes in the theory of gravity,” Classical and Quantum Gravity, vol. 26, no. 23, Article ID 235020, 2009.
- M. Sharif and M. F. Shamir, “Non-vacuum Bianchi types I and V in gravity,” General Relativity and Gravitation, vol. 42, no. 11, pp. 2643–2655, 2010.
- W. Hu and I. Sawicki, “Models of cosmic acceleration that evade solar system tests,” Physical Review D, vol. 76, Article ID 064004, 2007.
- A. De Felice and S. J. Tsujikawa, “ theories,” Living Reviews in Relativity, vol. 13, p. 3, 2010.
- T. P. Sotiriou and V. Faraoni, “ theories of gravity,” Reviews of Modern Physics, vol. 82, no. 1, pp. 451–497, 2010.
- T. Clifton, P. G. Ferreira, A. Padilla, and C. Skordis, “Modified gravity and cosmology,” Physics Reports, vol. 513, no. 1-3, pp. 1–189, 2012.
- H. E. Velten and T. R. Caramês, “Cosmological inviability of gravity,” Physical Review D, vol. 95, Article ID 123536, 2017.
- S. Nojiri and S. D. Odintsov, “Modified Gauss–Bonnet theory as gravitational alternative for dark energy,” Physics Letters B, vol. 631, no. 1-2, pp. 1–6, 2005.
- S. Nojiri, S. D. Odintsov, and O. . Gorbunova, “Dark energy problem: from phantom theory to modified Gauss-Bonnet gravity,” Journal of Physics A: Mathematical and General, vol. 39, no. 21, pp. 6627–6633, 2006.
- G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, and S. Zerbini, “Dark energy in modified Gauss-Bonnet gravity: late-time acceleration and the hierarchy problem,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 73, Article ID 084007, 2006.
- T. Chiba, “Generalized gravity and a ghost,” Journal of Cosmology and Astroparticle Physics, no. 3, pp. 107–113, 2005.
- M. Sharif and A. Ikram, “Energy conditions in gravity,” The European Physical Journal Special Topics, vol. 76, p. 640, 2016.
- G. Abbas, D. Momeni, M. Aamir Ali, R. Myrzakulov, and S. Qaisar, “Anisotropic compact stars in gravity,” Astrophysics and Space Science, vol. 357, no. 2, article 158, 2015.
- M. J. S. Houndjo, M. E. Rodrigues, D. Momeni, and R. Myrzakulov, “Exploring cylindrical solutions in modified gravity,” Canadian Journal of Physics, vol. 92, no. 12, pp. 1528–1540, 2014.
- M. E. Rodrigues, M. J. S. Houndjo, D. Momeni, and R. Myrzakulov, “A type of Levi–Civita solution in modified Gauss–Bonnet gravity,” Canadian Journal of Physics, vol. 92, no. 12, pp. 173–176, 2014.
- B. Wu and B.-Q. Ma, “Spherically symmetric solution of f(R,G) gravity at low energy,” Physical Review D, vol. 92, no. 4, Article ID 044012, 2015.
- M. Laurentis, M. Paolella, and S. Capozziello, “Cosmological inflation in F(R, G) gravity,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 91, no. 8, Article ID 083531, 2015.
- M. Sharif and A. Ikram, “Warm inflation in f(G) theory of gravity,” Journal of Experimental and Theoretical Physics, vol. 123, no. 1, pp. 40–50, 2016.
- M. Sharif and H. I. Fatima, “Noether symmetries in f(G) gravity,” Journal of Experimental and Theoretical Physics, vol. 122, no. 1, pp. 104–112, 2016.
- M. Sharif and H. I. Fatima, “Built-in inflation in f(G) gravity,” International Journal of Modern Physics D: Gravitation, Astrophysics, Cosmology, vol. 25, no. 1, Article ID 1650011, 12 pages, 2016.
- N. M. García, F. S. N. Lobo, J. P. Mimoso, and T. Harko, “F(G) modified gravity and the energy conditions,” Journal of Physics: Conference Series, vol. 314, no. 1, Article ID 012056, 2011.
- K. Atazadeh and F. Darabi, “Energy conditions in f(R,G) gravity,” General Relativity and Gravitation, vol. 46, no. 2, pp. 1–14, 2014.
- M. Sharif and A. Ikram, “Stability analysis of some reconstructed cosmological models in f(G,T) gravity,” Physics of the Dark Universe, vol. 17, pp. 1–9, 2017.
- M. Shamir and M. Ahmad, “Noether symmetry approach in gravity,” European Physical Journal, vol. 77, 55 pages, 2017.
- M. Sharif and A. Ikram, “Existence of Static Wormholes in Gravity,” https://arxiv.org/abs/1707.05162.
- M. Shamir and M. Ahmad, “Some exact solutions in gravity via Noether symmetries,” Modern Physics Letters, vol. A32, Article ID 1750086, 16 pages, 2017.
- X.-X. Wang, “Bianchi type-III string cosmological model with bulk viscosity in general relativity,” Chinese Physics Letters, vol. 22, no. 1, pp. 29–32, 2005.
- K. S. Thorne, “Primordial element formation, primordial magnetic fields, and the isotropy of the universe,” Astrophysical Journal, vol. 148, article 51, 1967.
- C. B. Collins and S. W. Hawking, “Why is the universe isotropic?” The Astrophysical Journal, vol. 180, pp. 317–334, 1973.
- M. Sharif and M. Zubair, “Dynamics of Bianchi I universe with magnetized anisotropic Dark Energy,” Astrophysics and Space Science, vol. 330, no. 1-2, pp. 399–405, 2010.
- S. Myrzakul, R. Myrzakulov, and L. Sebastiani, “Chaotic inflation in higher derivative gravity theories,” The European Physical Journal C, vol. 75, article 111, 2015.
- A. A. Starobinsky, “Future and origin of our Universe: modern view,” Gravitation and Cosmology, vol. 6, no. 3, pp. 157–163, 2000.
- R. R. Caldwell, “A phantom menace? Cosmological consequences of a dark energy component with super-negative equation of state,” Physics Letters, vol. 545, no. 1-2, pp. 23–29, 2002.
- J. Hogan, “Unseen Universe: welcome to the dark side,” Nature, vol. 448, no. 7151, pp. 240–245, 2007.
- P. S. Corasaniti, M. Kunz, D. Parkinson, E. J. Copeland, and B. A. Bassett, “The Foundations of observing dark energy dynamics with the Wilkinson Microwave Anisotropy Probe,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 70, Article ID 083006, 14 pages, 2004.
- J. Weller and A. M. Lewis, “Large-scale cosmic microwave background anisotropies and dark energy,” Monthly Notices of the Royal Astronomical Society, vol. 346, no. 3, pp. 987–993, 2003.
- M. A. MacCallum, “A class of homogeneous cosmological models III. Asymptotic behaviour,” Communications in Mathematical Physics, vol. 20, pp. 57–84, 1971.
- M. Cataldo and S. D. Campo, “Viscous cosmology in the Kasner metric,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 61, Article ID 128301, 2000.
- G. Cognola, M. Gastaldi, and S. Zerbini, “On the stability of a class of modified gravitational models,” International Journal of Theoretical Physics, vol. 47, no. 4, pp. 898–910, 2008.
- A. De Felice and S. Tsujikawa, “Construction of cosmologically viable gravity models,” Physics Letters B, vol. 675, no. 1, pp. 1–8, 2009.
- K. Bamba, S. D. Odintsov, L. Sebastiani, and S. Zerbini, “Finite-time future singularities in modified Gauss-Bonnet and gravity and singularity avoidance,” The European Physical Journal C, vol. 67, no. 1, pp. 295–310, 2010.
Copyright © 2017 M. Farasat Shamir. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP3.