Research Article | Open Access

Sumanta Chakraborty, "Field Equations for Lovelock Gravity: An Alternative Route", *Advances in High Energy Physics*, vol. 2018, Article ID 6509045, 6 pages, 2018. https://doi.org/10.1155/2018/6509045

# Field Equations for Lovelock Gravity: An Alternative Route

**Academic Editor:**Izzet Sakalli

#### Abstract

We present an alternative derivation of the gravitational field equations for Lovelock gravity starting from Newton’s law, which is closer in spirit to the thermodynamic description of gravity. As a warm up exercise, we have explicitly demonstrated that, projecting the Riemann curvature tensor appropriately and taking a cue from Poisson’s equation, Einstein’s equations immediately follow. The above derivation naturally generalizes to Lovelock gravity theories where an appropriate curvature tensor satisfying the symmetries as well as the Bianchi derivative properties of the Riemann tensor has to be used. Interestingly, in the above derivation, the thermodynamic route to gravitational field equations, suited for null hypersurfaces, emerges quiet naturally.

#### 1. Introduction

Equivalence principle acts as the guiding lighthouse to understand how matter fields behave in a curved spacetime background. Unfortunately, there exists no such principle which helps to answer the opposite, namely, how matter fields curve the spacetime [1–7]. Lack of such principle resulted in the vast landscape of various alternative gravity theories among which general relativity remains the most useful one. Despite this large variety of alternative gravitational theories one can do better by imposing some physical requirements on the systems of interest. In particular, the restriction to the class of theories having utmost second-order derivatives of the metric is a judicious one, since this helps to overcome the well known Ostrogradski instability by evading the existence of any ghost modes in the theory [8–10]. Surprisingly enough the above criteria turn out to be very interesting ones, as they single out a very specific class of unique gravitational theories, known in the literature as the Lanczos-Lovelock models of gravity [11–14]. In particular, all these models satisfy the criterion that Bianchi derivative of the field tensor identically vanishes [15, 16]. One can add to this list a large number of interesting properties; a few among them are as follows: (a) the action functional for Lanczos-Lovelock gravity can be broken into a bulk part and a surface contribution, which is closely related to the Wald entropy associated with the black hole horizon [17]; (b) the field equations in Lanczos-Lovelock gravity can be expressed as a thermodynamic identity on an arbitrary null surface [18–26]; (c) the difference between a suitably defined surface and bulk degrees of freedom in the context of Lanczos-Lovelock gravity can be interpreted as the evolution of the spacetime [27]; (d) the surface and bulk terms in the Lanczos-Lovelock gravity are related by a holographic relation between them [28, 29].

In the standard textbook treatment, one first introduces the gravitational Lagrangian , which is for general relativity, while the Lagrangian is a polynomial in the Riemann tensor for general Lanczos-Lovelock gravity. Variation of the action functional due to arbitrary variation of the metric, with appropriate boundary conditions [30, 31], leads to the corresponding field equations for gravity, namely, the rule by which matter tells the spacetime how to curve (see also [32]). However, as emphasized earlier there is no physical principle which guides us to the precise mathematical structure of the gravitational action; what is more the field equations obtained equate matter which is intrinsically* quantum* with spacetime geometry describing gravity* classically*. Despite this conceptual discomfort, there exist two additional representations of the gravitational dynamics, which are intrinsically observer dependent. The first one uses timelike observers with four-velocity and equates two quantities where this observer measures in the geometric as well as in the matter sector, leading to a scalar equationHere stands for the gravitational constant in four dimensions and enforcing this equation for all observers with four-velocity leading to Einstein’s equations . The other approach uses null vectors (i.e., ) and leads to the following scalar equation:In this situation as well validity of the above expression for all null vectors , along with two times contracted Bianchi identity and covariant conservation of matter energy momentum tensor, amounts to furnishing the ten components of Einstein’s equations. Following [33], in this article, we would like to derive (1) and (2) from a geometrical point of view for Lanczos-Lovelock theories of gravity. We would also argue why the route taken in this work is a natural one compared to the standard derivation of Einstein’s equations from the gravitational action. We have organized the paper as follows: In Section 2 we briefly review the derivation of Einstein’s equations, which will be generalized in Section 3 to Lanczos-Lovelock theories of gravity. Finally we provide the concluding remarks before presenting a derivation regarding the curvature tensor associated with the Lanczos-Lovelock theory of gravity in Appendix.

#### 2. Einstein’s Equations following Newton’s Path

To set the stage for Lanczos-Lovelock theories of gravity we would like to briefly review the corresponding situation in general relativity, which will mainly follow from [33]. In Newtonian theory one describes gravity using the gravitational potential . Given the matter density , the dynamics of the gravitational field is being determined through Poisson’s equation . When one invokes principle of equivalence, the potential is readily identified with components of the metric, in particular, . One then interprets Poisson’s equation in metric language as , where matter energy density is interpreted as the time-time component of the divergence free, symmetric, matter energy momentum tensor . In order to have a relativistic generalization one naively makes the following replacement . This immediately suggests looking for second derivatives of the metric which is a second rank tensor as well as divergence free, paving the way to the Einstein tensor. On the other hand, for the right hand side the natural choice being the matter energy momentum tensor , one ends up equating Einstein’s tensor with the matter energy momentum tensor resulting in Einstein’s equations (for a more detailed discussion, see [33]).

However, we could have just followed Newton’s path and look for relativistic generalization of Poisson’s equation itself, possibly leading to a scalar equation describing gravity in general relativity. As a first step one must realize that the right hand side of Poisson’s equation is intrinsically observer dependent, as it is the energy density of some matter field measured by an observer with some four-velocity . Given the symmetric and conserved matter energy momentum tensor , the matter energy density as measured by an observer with four-velocity is . Hence the same observer dependence must continue to exist in general relativity and the right hand side of the desired general relativistic equation for gravity should be [33].

To derive the left hand side, that is, analogue of , we note that this requires deriving a scalar object which involves two spatial derivatives of the metric (since this is what is responsible for the term in nonrelativistic limit) and necessarily depends on the four-velocity . The only tensor depending on two derivatives of the metric that can be constructed corresponds to the curvature tensor . In order to obtain double spatial derivatives acting on the metric one has to project all the components of the curvature tensor on the plane orthogonal to using and hence construct a scalar thereof. To keep generality, instead of looking for a timelike vector , we will concentrate on a particular vector , with norm . Then one can immediately introduce the following tensor:for which the final property ensures that it is a projector. Then one can define a projected Riemann curvature with respect to the vector by projecting all the indices of the Riemann curvature tensor, leading to (by the very definition, the Riemann tensor has a generic form in local inertial frame; thus the above projection ensures that has only double spatial derivatives when is a timelike vector; all time derivatives appearing in are single in nature and hence vanish in the local inertial frame all together; this is the prime motivation of introduction of this projected Riemann tensor; the same can be ascertained from the Gauss-Codazzi equation as well; see [1, 5])The scalar constructed out of this projected Riemann tensor becomesNote that we have thrown away another term involving , since due to antisymmetry properties of the Riemann tensor the above quantity identically vanishes. Furthermore, the above scalar* by* construction has only two derivatives of the metric and is completely spatial in the instantaneous rest frame of the timelike observer. Hence is the relativistic generalization of and thus must be equal to the corresponding matter energy density, which in this case corresponds to , where is the matter energy momentum tensor. Thus finally one obtains the following equation:Surprisingly, through the above exercise we have achieved two results in one go. Firstly, for , that is, for unit normalized timelike vectors we get back the equation . Then requiring these equations to hold for all timelike observers will lead to . At this stage it will be worthwhile to mention that a similar approach as the above one was taken in [34] to derive Einstein’s equations. However the key difference between the approach in [34] and in this work is the use of the projection tensor judiciously. For our approach the use of projection tensor is of utmost importance in contrast to [34].

Finally the limit would lead to the null version of the above equation, which reads . In this case besides demanding the validity of the above equation for all null vectors, it is important to use contracted Bianchi identity as well as covariant conservation of energy momentum tensor, leading to . Note that in the second situation a cosmological constant has automatically come into existence.

#### 3. Newton Leads the Way to Lovelock

In view of the above result, one immediately asks for the corresponding situation in Lovelock gravity: can field equations of Lovelock gravity be derived following the above procedure? This is what we will explore in this section and shall show that one can indeed derive the field equations for Lovelock gravity following an equivalent procedure. Before jumping into the details let us mention some structural aspects of Lovelock gravity. Briefly speaking, Lovelock gravity corresponds to a class of gravitational Lagrangians which are polynomial in the Riemann curvature tensor yielding field equations which are of second order in the metric. The analytical form of such a polynomial (also called a pure Lovelock term) of order involves Riemann curvature tensors contracted appropriately, such thatThe above relation defines the tensor associated with the th order Lanczos-Lovelock gravity, having all the symmetries of the Riemann tensor with the following algebraic structure:The tensor satisfies an additional criterion , which ensures that the field equations derived from this Lagrangian are of second order. In what follows we will exclusively concentrate on the th order Lanczos-Lovelock gravity and hence shall remove the symbol “” from the superscripts of various geometrical expressions.

One can now start from Poisson’s equation and fix the right hand side to be , which is the matter energy density associated with the vector field . In order to get the left hand side we must construct an appropriate curvature tensor suited for the Lovelock gravity and project it using the projector introduced in (3) before constructing a scalar out of it. There are two possible choices for such a curvature tensor among which we will discuss the simpler one in the following, while the complicated one is deferred to Appendix 4. The curvature tensor described here was first introduced in [15, 16] and is defined as follows:where is the th order Lovelock polynomial and . Since the th order Lanczos-Lovelock Lagrangian depends on powers of Riemann, the above definition for ensures that it depends on powers of Riemann tensor and hence exactly coincides with (8). Furthermore, the tensor can be easily generalized to the full Lovelock polynomial by just adding over different values, but we will concentrate on a single term in the full Lanczos-Lovelock Lagrangian. Also note that the above defined “Lovelock” Riemann tensor has all the symmetries of the original “Einstein” Riemann tensor . Further it satisfies the contracted Bianchi identity, which will be sufficient for our purpose. Given the above th order “Lovelock” Riemann, one can project it in the plane orthogonal to and obtain the following scalar:Explicit evaluation of the above scalar can be performed keeping in mind that the tensor is antisymmetric under exchange of the indices and , respectively, leading toThus equating (11) to the matter energy density through we finally obtainFor unit normalized timelike vectors, and , leading to , where is the analogue of Einstein tensor in Lovelock gravity. On the other hand, for the null vectors we arrive at , where the tensor is the analogue of Ricci tensor in Lovelock gravity. Hence even in the case of Lovelock gravity, if one assumes that holds for all timelike observers, the field equations for Lovelock gravity follow. While in the null case, besides demanding the validity of for all null vectors, one has to use the Bianchi identity associated with Lovelock theories as well as covariant conservation of matter energy momentum tensor to arrive at the Lovelock field equations , which inherit the cosmological constant as well. Note that the above result has been derived in the context of th order Lanczos-Lovelock gravity, which can be generalized to the general Lanczos-Lovelock Lagrangian in a straightforward manner. This results in . Interestingly, the cosmological constant besides being generated as an integration constant of the field equations also gets contribution from the term in the Lanczos-Lovelock gravity. Therefore one may choose this term in the Lagrangian appropriately to arrive at the present small value of the cosmological constant. Therefore we conclude that the derivation of field equations for gravity can always be achieved starting from Poisson’s equation and subsequently projecting a suitable curvature tensor, whether it is Einstein gravity or Lovelock.

#### 4. Concluding Remarks

By equivalence principle gravity manifests itself by curving the spacetime, which the material particles follow. In particular one can invoke special relativity in locally freely falling frame and hence write down the laws of motion in curvilinear coordinates, thus describing motion in curved spacetime. The notion of locally freely falling observer brings in intrinsic observer dependence in the theory and introduces observers for whom a local spacetime region is causally inaccessible, known as a local Rindler observer. Remarkably the local vacuum state of a test quantum field (as fit for local inertial observers) will appear as thermal to the local Rindler observer [35, 36]. If any matter field (characterized by matter energy momentum tensor ) crosses the local Rindler horizon, it will appear to be thermalized by the Rindler observer (since the matter will take infinite time to reach the horizon) and the corresponding heat density is being given by (for a perfect fluid the above quantity is given by , which by Gibbs-Duhem relation is the matter heat density), where is the null normal to the horizon. Note that the above heat density for matter is invariant under the transformation .

At this stage, one can ask a natural question, “what about heat density of gravity?”. Surprisingly, one can answer the same in the above setting. Considering a general null surface it turns out that one can interpret as the heat density of the spacetime. This originates from the fact that one can have a one to one correspondence between and the viscous dissipation term , where is the shear of the null congruence and is its expansion, with and being shear and bulk viscous coefficients. Thus the term is related to heating of the spacetime [37]. Given this thermodynamic backdrop, it is clear that the Einstein’s equations when written as (2) not only yield the geometrical input of the gravitational theory but are also physically well-motivated since the equality of (2) can be thought of as an equilibrium situation, where the heat produced by gravity is being compensated by that of matter.

Thus the equation arises more naturally from the relativistic generalization of Newton’s law and the usefulness of writing Einstein’s equations in this manner stems from the fact that one might interpret both sides of these equations independently and the equations themselves follow due to a balancing act performed by spacetime itself [33, 38, 39]. We would like to reiterate that the field equations derived in this context are purely geometrical and follow Newton’s path. One first realizes that energy density associated with any material body is intrinsically observer dependent and surprisingly one can construct a tensor (again dependent on observer) which contains spatial derivatives of the metric alone. Keeping this as a curved spacetime generalization of Newton’s law one uniquely arrives at Einstein’s equations when one makes use of all observers (or all the null surfaces). This shows that the most natural generalization of Newton’s law to curved spacetime is (of course, leading to Einstein’s equations, but at a secondary level) bolstering the claim that gravity is intrinsically a thermodynamic phenomenon.

#### Appendix

#### An Alternative Riemann Tensor for Lovelock Gravity

In Lovelock gravity it is possible to define two tensors having the symmetry properties of Riemann and satisfying Bianchi identity. The first one corresponds to defining a tensor for th order Lovelock polynomial by multiplying such curvature tensors and then an alternating tensor of rank as [40]Note that the above tensor is completely antisymmetric in both upper and lower indices.

Let us now project all the indices on a lower dimensional spacelike hypersurface using the projection tensor , such that one obtains another tensor, but whose inner product with the normal identically vanishes. Hence one arrives atOne can immediately construct a scalar out of the projected tensor, leading toThere would be two terms contributing to the above expression, one when all the projectors are related by Kronecker deltas and when one of the projector is replaced by while all the others are replaced by Kronecker deltas. In any other case, for example if two projectors are replaced by the normal, then it would identically vanish, thanks to the completely antisymmetric nature of . Thus finally we obtainOne can now use the following identities:as well assuch that one arrives atwhere stands for the analogue of Ricci tensor associated with the Lovelock gravitational action. Hence this particular Lovelock Riemann tensor reproduces the Lovelock field equations for gravity as well if the procedure outlined above is being followed. The above exercise explicitly demonstrates the robustness of the idea presented here.

#### Conflicts of Interest

The author states that there are no conflicts of interest.

#### Acknowledgments

Sumanta Chakraborty gratefully acknowledges the help of T. Padmanabhan for suggesting this project and also for helpful discussions throughout it. He also thanks Naresh Dadhich, Kinjalk Lochan, and Krishnamohan Parattu for various fruitful discussions. A part of this work was completed while the author was visiting Albert Einstein Institute in Golm, Germany, and the author gratefully acknowledges the warm hospitality there. Research of Sumanta Chakraborty is supported by the SERB-NPDF Grant (PDF/2016/001589) from DST, Government of India.

#### References

- T. Padmanabhan,
*Gravitation: Foundations and Frontiers*, Cambridge University Press, Cambridge, UK, 2010. - L. D. Landau and E. M. Lifshitz,
*The Classical Theory of Fields*, vol. 2 of*Course of Theoretical Physics Series*, Butterworth-Heinemann, 4th edition, 1980. View at: MathSciNet - A. S. Eddington,
*The Mathematicl Theory of Relativity*, Cambridge University Press, Cambridge, UK, 2nd edition, 1960. View at: MathSciNet - W. K. Misner, S. Thorne, and J. A . Wheeler,
*Gravitation*, W. H. Freeman and Company, 3rd edition, 1973. View at: MathSciNet - E. Poisson,
*A Relativist Toolkit*, Cambridge University Press, Cambridge, UK, 1st edition, 2004. View at: Publisher Site | MathSciNet - R. M. Wald,
*General Relativity*, University of Chicago Press, 1st edition, 1984. View at: Publisher Site | MathSciNet - S. W. Hawking and G. F. R. Ellis,
*The Large Scale Structure of Space-Time*, Cambridge Monographs on Mathematical Physics, Cambridge University Press, London, UK, 1973. View at: MathSciNet - T.-j. Chen, M. Fasiello, E. A. Lim, and A. J. Tolley, “Higher derivative theories with constraints: exorcising Ostrogradski's ghost,”
*Journal of Cosmology and Astroparticle Physics*, vol. 2013, no. 2, 2013. View at: Google Scholar | MathSciNet - H. Motohashi, K. Noui, T. Suyama, M. Yamaguchi, and D. Langlois, “Healthy degenerate theories with higher derivatives,”
*Journal of Cosmology and Astroparticle Physics*, vol. 2016, 2016. View at: Publisher Site | Google Scholar - R. Klein and D. Roest, “Exorcising the Ostrogradsky ghost in coupled systems,”
*Journal of High Energy Physics*, no. 7, 2016. View at: Google Scholar | MathSciNet - C. Lanczos, “A remarkable property of the Riemann-Christoffel tensor in four dimensions,”
*Annals of Mathematics*, vol. 39, no. 4, pp. 842–850, 1938. View at: Publisher Site | Google Scholar - C. Lanczos, “Electricity as a natural property of Riemannian geometry,”
*Physical Review A: Atomic, Molecular and Optical Physics*, vol. 39, no. 4, pp. 716–736, 1932. View at: Publisher Site | Google Scholar - D. Lovelock, “The Einstein tensor and its generalizations,”
*Journal of Mathematical Physics*, vol. 12, pp. 498–501, 1971. View at: Publisher Site | Google Scholar | MathSciNet - T. Padmanabhan and D. Kothawala, “Lanczos-Lovelock models of gravity,”
*Physics Reports*, vol. 531, no. 3, pp. 115–171, 2013. View at: Publisher Site | Google Scholar | MathSciNet - N. Dadhich, “Characterization of the Lovelock gravity by Bianchi derivative,”
*Pramana—Journal of Physics*, vol. 74, no. 6, pp. 875–882, 2010. View at: Publisher Site | Google Scholar - X. O. Camanho and N. Dadhich, “On Lovelock analogs of the Riemann tensor,”
*The European Physical Journal C*, vol. 76, no. 3, 2016. View at: Publisher Site | Google Scholar - D. Kothawala, T. Padmanabhan, and S. Sarkar, “Is gravitational entropy quantized?”
*Physical Review D: Particles, Fields, Gravitation and Cosmology*, vol. 78, no. 10, Article ID 104018, 5 pages, 2008. View at: Publisher Site | Google Scholar | MathSciNet - A. Paranjape, S. Sarkar, and T. Padmanabhan, “Thermodynamic route to field equations in Lanczos-LOVelock gravity,”
*Physical Review D: Particles, Fields, Gravitation and Cosmology*, vol. 74, no. 10, Article ID 104015, 2006. View at: Publisher Site | Google Scholar | MathSciNet - M. Akbar and R. G. Cai, “Friedmann equations of FRW universe in scalar–tensor gravity, f(R) gravity and first law of thermodynamics,”
*Physics Letters B*, vol. 635, no. 1, pp. 7–10, 2006. View at: Publisher Site | Google Scholar - R. G. Cai and S. P. Kim, “First law of thermodynamics and Friedmann equations of Friedmann-Robertson-Walker universe,”
*Journal of High Energy Physics*, vol. 2005, article 050, 2005. View at: Publisher Site | Google Scholar - X. H. Ge, “First law of thermodynamics and Friedmann-like equations in braneworld cosmology,”
*Physics Letters B*, vol. 651, no. 1, pp. 49–53, 2007. View at: Publisher Site | Google Scholar | MathSciNet - Y. Gong and A. Wang, “Friedmann Equations and Thermodynamics of Apparent Horizons,”
*Physical Review Letters*, vol. 99, no. 21, 2007. View at: Publisher Site | Google Scholar - S. Wu, B. Wang, and G. Yang, “Thermodynamics on the apparent horizon in generalized gravity theories,”
*Nuclear Physics B*, vol. 799, no. 3, pp. 330–344, 2008. View at: Publisher Site | Google Scholar - S. Chakraborty, “Lanczos-Lovelock gravity from a thermodynamic perspective,”
*Journal of High Energy Physics*, vol. 29, 2015. View at: Google Scholar | MathSciNet - D. Hansen, D. Kubiznak, and R. Mann, “Horizon Thermodynamics from Einsteins Equation of State,”
*Physics Letters B*, vol. 771, pp. 277-–280, 2017. View at: Google Scholar - R. Dey, S. Liberati, and A. Mohd, “Higher derivative gravity: Field equation as the equation of state,”
*Physical Review D: Particles, Fields, Gravitation and Cosmology*, vol. 94, no. 4, 2016. View at: Publisher Site | Google Scholar - T. Padmanabhan, “Equipartition of energy in the horizon degrees of freedom and the emergence of gravity,”
*Modern Physics Letters A*, vol. 25, no. 14, pp. 1129–1136, 2010. View at: Publisher Site | Google Scholar | MathSciNet - G. A. Mukhopadhyay and T. Padmanabhan, “Holography of gravitational action functionals,”
*Physical Review D: Particles, Fields, Gravitation and Cosmology*, vol. 74, no. 12, Article ID 124023, 15 pages, 2006. View at: Publisher Site | Google Scholar - D. Dieks, J. van Dongen, and S. de Haro, “Emergence in holographic scenarios for gravity,”
*Studies in History and Philosophy of Science Part B: Studies in History and Philosophy of Modern Physics*, vol. 52, pp. 203–216, 2015. View at: Publisher Site | Google Scholar - K. Parattu, S. Chakraborty, B. R. Majhi, and T. Padmanabhan, “A boundary term for the gravitational action with null boundaries,”
*General Relativity and Gravitation*, vol. 48, no. 7, Art. 94, 28 pages, 2016. View at: Publisher Site | Google Scholar | MathSciNet - S. Chakraborty, “Boundary Terms of the Einstein–Hilbert Action,” in
*Gravity and the Quantum*, vol. 187 of*Fundamental Theories of Physics*, pp. 43–59, Springer International Publishing, 2017. View at: Publisher Site | Google Scholar - S. Deser, “Gravity from self-interaction redux,”
*General Relativity and Gravitation*, vol. 42, no. 3, pp. 641–646, 2010. View at: Publisher Site | Google Scholar | MathSciNet - T. Padmanabhan, “The atoms of space, gravity and the cosmological constant,”
*International Journal of Modern Physics D: Gravitation, Astrophysics, Cosmology*, vol. 25, no. 7, Article ID 1630020, 33 pages, 2016. View at: Publisher Site | Google Scholar | MathSciNet - N. Straumann,
*General relativity*, Springer, Zurich, Switzerland, 2013. - P. C. Davies, “Scalar production in Schwarzschild and Rindler metrics,”
*Journal of Physics A: Mathematical and General*, vol. 8, no. 4, pp. 609–616, 1975. View at: Publisher Site | Google Scholar - W. G. Unruh, “Notes on black-hole evaporation,”
*Physical Review D: Particles, Fields, Gravitation and Cosmology*, vol. 14, no. 4, pp. 870–892, 1976. View at: Publisher Site | Google Scholar - S. Chakraborty and T. Padmanabhan, “Thermodynamical interpretation of the geometrical variables associated with null surfaces,”
*Physical Review D: Particles, Fields, Gravitation and Cosmology*, vol. 92, no. 10, Article ID 104011, 30 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet - Z. Wang, W. Ai, H. Chen, and J. Deng, “Cosmological model from emergence of space,”
*Physical Review D: Particles, Fields, Gravitation and Cosmology*, vol. 92, no. 2, 2015. View at: Publisher Site | Google Scholar - C. E. Navia, “On the Davies-Unruh effect in a wide range of temperatures,” https://arxiv.org/abs/1701.03442. View at: Google Scholar
- D. Kastor, “The Riemann-Lovelock curvature tensor,”
*Classical and Quantum Gravity*, vol. 29, no. 15, 155007, 9 pages, 2012. View at: Publisher Site | Google Scholar | MathSciNet

#### Copyright

Copyright © 2018 Sumanta Chakraborty. 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 SCOAP^{3}.