Thick brane structures in generalized hybrid metric-Palatini gravity

In this work, we study 5-dimensional braneworld scenarios in the scalar-tensor representation of the generalized hybrid metric-Palatini gravitational theory. We start by considering a model for a brane supported purely by the gravitational scalar fields of the theory and then consider other distinct cases where the models are also supported by an additional matter scalar field. We investigate the stability of the gravity sector and show that the models are all robust against small fluctuations of the metric. In particular, in the presence of the additional scalar field, we find that the profile of the gravitational zero mode may be controlled by the parameters of the model, being also capable of developing internal structure.


Introduction
Since its development in 1915, General Relativity (GR) has been a very successful theory, at least for local tests. However, the discovery of the late-time cosmic acceleration [1,2] has spurred research in modified dynamics at large scales. In this context, the possibility that the Hilbert-Einstein term is supplemented with more general combinations of curvature invariants has been extensively explored [3][4][5][6][7][8][9][10][11]. Indeed, it was soon found that the usual metric formulation, which considers that the metric is the fundamental field, differs generically from its Palatini (or metric-affine) counterpart [12], where here the metric and the connection are assumed to be the two fundamental fields of the theory. The metric a e-mail: joaoluis92@gmail.com (corresponding author) b e-mail: ferreira.douglasdaf@gmail.com c e-mail: dbazeia@gmail.com d e-mail: fslobo@fc.ul.pt approach leads to higher-order derivative equations, contrary to the field equations in the Palatini formulation, which are second-order. We do emphasize that the purely f (R) metric and metric-affine formalisms coincide in GR, however, the two formalisms lead to different results considering more generic functions f (R) [12]. For instance, the scalar-tensor representation of f (R) gravity is useful to illustrate the differences between the metric and Palatini approaches. In the metric formalism, the scalar field behaves as a dynamical field, which satisfies a modified Klein-Gordon equation with self-interactions that essentially depend on the form of f (R). However, the scalar field φ should have a very low mass, which implies a long interaction range, in order to have an impact at large scales; nevertheless, light scalars do indeed have an impact at smaller scales, and are strongly constrained by local observations, at the laboratory and Solar System scales, unless screening mechanisms are invoked [13][14][15]. In the scalar-tensor representation of the Palatini formalism, the scalar field satisfies an algebraic function of the trace of the matter stress-energy tensor, which lead to undesired gradient instabilities [16][17][18][19].
However, these difficulties may be avoided within a hybrid variation of these theories, in which the purely metric Einstein-Hilbert action is supplemented with a metricaffine correction term [20][21][22][23][24]. More specifically, an interesting aspect of these theories is the possibility of generating long-range forces without conflicts with the local tests and without invoking screening mechanisms. The possibility of expressing these hybrid f (R) metric-Palatini theories using a scalar-tensor representation simplifies the analysis of the field equations and the construction of solutions. It is interesting to note that this theory, which in the linear approach takes the form R + f (R), one retains the positive results through the Einstein-Hilbert term R and the addi-tional gravitational corrections are given by the metric-affine f (R) component, where the Palatini curvature scalar R is constructed in terms of an independent connection. A wide variety of applications of this hybrid metric-Palatini theory has been explored, namely, in considering that dark matter is a geometric effect of modified gravity [25,26], in exploring the Cauchy problem [27] and the Noether symmetries [28], black hole, wormhole and stellar solutions [29][30][31][32][33][34], the Einstein static Universe [35], string-like configurations [36,37], and thick branes [38], among others [39].
The linear hybrid metric-Palatini theory can be further generalized, where the gravitational action depends on a general function of both the metric and Palatini curvature scalars [40]. This extension has also received attention with a plethora of applications, namely, cosmological solutions [41,42], weak-field phenomenology [43], wormhole [44] and black hole [45] solutions. In this work, we consider the possibility to study braneworld structures in the generalized hybrid metric-Palatini gravity. In this context, as is wellknown, the Randall-Sundrum braneworld model originally proposed [46] was soon generalized [47][48][49] to describe thick braneworld scenarios in the presence of scalar fields, which are included to source the warped 5-dimensional AdS 5 geometry with a single extra spatial dimension of infinite extent. Soon after, several authors investigated the braneworld scenario with a diversity of motivations; see, e.g., Refs.  and references therein for many distinct possibilities to implement investigations on thick braneworlds, including the absence of scalar fields, the presence of two scalar fields with standard dynamics, the case of tachyonic and other generalized dynamics, asymmetric thick brane, and also the f (R), Gauss-Bonnet and several other possibilities of modifications of GR [5][6][7]9].
In order to deal with the 5-dimensional generalized hybrid metric-Palatini braneworld scenario, the study developed in the present work is outlined in the following manner: in Sect. 2 we introduce the model and write the equations of motion. Furthermore, in Sect. 3 we consider the warped geometry with a single extra dimension of infinite extent, rewrite the equations of motion on general grounds and show how to get to the thick braneworld scenarios in the case of standard GR. We also describe distinct solutions of current interest, in the absence and in the presence of an extra field χ . In Sect. 4, we discuss the robustness of the geometric sector in the generalized hybrid metric-Palatini braneworld theory. Finally, in Sect. 5, we summarize our results and conclude. Throughout this work, we implement all the calculations considering a system of units for which c = 1, using capital latin indexes {M, N , ...} running from 0 to 4 and greek indexes {μ, ν, ...} running from 0 to 3.

Action and field equations
The generalized hybrid metric-Palatini theory [40,41] in 4+1 dimensional gravity is described by an action functional S of the form where κ 2 ≡ 8π G 5 , G 5 is the 5-dimensional Newtonian constant, is a 5-dimensional spacetime manifold on which we define the coordinate set x M , g is the determinant of the metric g M N , f is an arbitrary function of the Ricci scalar, where∇ denotes the covariant derivative written in terms of the independent connectionˆ P M N . Recalling that √ −g represents a scalar density of weight 1, we have that∇ P √ −g = 0, and thus Eq. (3) implies the existence of a new metric tensor h M N = (∂ f /∂R) g M N , conformally related to the metric g M N with a conformal factor ∂ f /∂R, for which the connec-tionˆ P M N is the Levi-Civita connection, i.e., one can writê where ∂ M denotes partial derivatives. The conformal relation between the metrics h M N and g M N implies that the two Ricci scalars R M N and R M N , assumed a priori as independent, are in fact related by the expression where ∇ M denotes covariant derivatives written in terms of the connection P M N , = ∇ P ∇ P is the D'Alembert oper-ator, and the subscripted f R denote derivatives of f with respect to R. An equivalent scalar-tensor representation of Eq. (1) can be obtained via the definition of two scalar fields ϕ and ψ and a scalar potential V (ϕ, ψ) as Inserting both Eqs. (6) and (7) in Eq. (1), using the trace of Eq. (5) to cancel the factor R, writing f R = −ψ and neglecting the boundary term proportional to ψ, we arrive to the scalar-tensor representation of the generalized hybrid metric-Palatini gravity as Without further consideration of the matter fields χ , there are three independent variables in Eq. (8), these are the metric g M N and the two scalar fields ϕ and ψ. Varying Eq. (8) with respect to the metric g M N yields the modified field equations where G M N is the Einstein's tensor, and T M N is the stressenergy tensor defined in the usual manner as Two equations of motion describing the dynamics of the scalar fields ϕ and ψ can also be obtained as follows. We start with a variation of Eq. (8) with respect to ϕ and ψ. The resultant equations of motion will depend on the Ricci scalar R. Taking the trace of Eq. (9), we are able to cancel the terms depending on R but we are left with a system of two coupled differential equations both depending on ϕ and ψ. These two equations can be algebraically manipulated to isolate the terms ϕ and ψ, so to isolate the dynamics of each scalar field in its own equation. The resultant equations of motion for ϕ and ψ are respectively, where the subscripted potentials V ϕ and V ψ denote derivatives of the potential V (ϕ, ψ) with respect to the scalar fields ϕ and ψ, respectively, and T = g M N T M N is the trace of the stress-energy tensor.

Static and flat brane models with scalar field matter
In this section, we will consider matter to be described by a single dynamical scalar field χ with an associated interaction potential U (χ ). The matter action describing this distribution of matter is thus The stress-energy tensor T M N associated with this matter distribution can be computed via a variation of Eq. (13) with respect to the scalar field χ and the definition of T M N provided in Eq. (10). We thus obtain Furthermore, and similarly to the scalar fields ϕ and ψ, an equation of motion for the field χ can be obtained varying Eq. (13) with respect to χ . As the field χ is minimally coupled to the metric g M N , a dynamical equation for χ is immediately obtained with no need for further manipulations. The resultant equation is where the subscripted potential U χ denotes a derivative of the potential U with respect to the scalar field χ . Our complete system of equations thus consists of Eqs. (9), (11), (12) and (15), with the stress-energy tensor T M N given by Eq. (14). For the purpose of this paper, let us consider the static 5-dimensional line element where A (y) is called the warp function, η μν is the 4dimensional Minkowski metric given by η μν = diag (−1, 1, 1, 1), and y represents the extra 5th dimension of infinite extent. We shall assume that both the gravitational scalar fields ϕ = ϕ (y) and ψ = ψ (y) and the matter scalar field χ = χ (y) are constant throughout the 4-dimensional spacetime and vary solely across the extra dimension y. Furthermore, given the isotropy of the 4-dimensional part of the metric, only two independent field equations arise. Inserting these assumptions and the metric from Eq. (16) into the system of Eqs. (9), (11), (12) and (15) yields for the two independent field equations, for the gravitational scalar fields ϕ and ψ, respectively, and for the matter scalar field χ . One can prove that Eq. (21) is not an independent equation for the system in the following manner: take a derivative with respect to y of Eq. (18), use Eq. (17) to cancel the terms depending on ζ , then use Eqs. (19) and (20) to cancel the terms depending on ϕ and ψ respectively, and finally use Eq. (18) itself to cancel the term depending on V (ϕ, ψ). The result of the stated algebraic manipulations is Eq. (21), thus proving its dependence on the remaining equations. Consequently, the system of Eqs. (17) to (20) fully describes the system in study. The system consists of four independent equations for the six independent variables A, ϕ, ψ, χ , V and U .

The standard GR case
The well-known GR standard case can be obtained from the system of Eqs. (17) to (21). This particular case corresponds to a form of the function f (R, R) as Using the definition of the scalar field provided in Eq. (6), one verifies that ϕ = 1 and ψ = 0. Also, from Eq. (7), one verifies that in this case V = 0. Thus, to obtain the standard case, one takes the limit ϕ → 1, ψ → 0, and V → 0 of the modified field equations in Eqs. (17) and (18), which become respectively where we have used Eq. (24) to cancel the dependency of Eq. (23) in A . On the other hand, the particular form of the function f for the standard case given by Eq. (22) does not allow for an equivalent scalar-tensor representation of the theory due to the fact that the determinant of its Jacobian matrix vanishes. This implies that the relationship between the scalar fields ϕ and ψ with R and R is degenerate. The equations of motion for the scalar fields, i.e., Eqs. (19) and (20) are thus effectively removed from the system. Furthermore, one can prove that the equation of motion for the scalar field χ given by Eq. (21) is not independent of the two equations above. To do so, one takes the derivative of Eq. (24) with respect to y and uses Eq. (23) to cancel the term depending on A , obtaining Eq. (21) as a result.
It is common to encounter Eqs. (23) and (24) in the literature with κ 2 = 2. To ease the comparison between our results in the upcoming sections and the available literature, we shall consider κ 2 = 2 from this point onward.

Solution without matter (χ = 0)
For simplicity, let us start by considering a brane model supported solely by the scalar fields ϕ and ψ, i.e., without the matter field χ , which results in Eq. (21) being automatically satisfied. Furthermore, we shall assume that the dependency of the potential V in the scalar fields ϕ and ψ is of the form V (ϕ − ψ), as commonly considered in literature [41,44]. This choice is made so that the potential is a function of the coupling between R and the scalar fields ϕ and ψ in Eq. (8). In this case, the partial derivatives of V become related by V ϕ = −V ψ ≡V . Subtracting Eq. (18) from Eq. (17) and imposing the assumptions mentioned above leads to On the other hand, the equations for the scalar fields given in Eqs. (19) and (20) become Finally, an equation relating V toV arises from a derivative of V (ϕ (y) , ψ (y)) with respect to y using the chain rule, which results in The system of Eqs. (25) to (28) consists of a system of four independent equations to the five unknowns ϕ, ψ, V , V , and A. Thus, the system is under-determined and an extra constraint must be imposed to close the system. We choose to specify the form of the warp function A to be the usual where A 0 and k are constants and A 0 in particular must be defined positive. Let us focus on even solutions for the scalar fields ϕ and ψ (and consequently V ). These solutions must satisfy the boundary conditions ϕ (0) = ψ (0) = 0. Furthermore, let us denote ϕ (0) = ϕ 0 and ψ (0) = ψ 0 their values at the origin. In the same way, we will have V (0) = 0 and V (0) = V 0 . Inserting these boundary conditions into Eq. (25) yields In order to preserve the positivity of the factor ϕ − ψ in Eq.
(8), we should impose ϕ 0 > ψ 0 . This ensures that at y = 0 we have ϕ (0) > ψ (0). We shall look for solutions for which the positiveness of this factor is maintained throughout the entire range of y.
As can be seen from Eq. (27), imposing a boundary condition ψ (0) = 0 will also set ψ (0) = 0, and thus the only possible even solution for ψ is the constant solution ψ = ψ 0 . We shall thus consider the constant ψ 0 as a free parameter and do not provide plots for the function ψ (y). On the other hand, the numerical solutions for ϕ (y) are plotted in Fig. 1. These solutions satisfy the boundary condition ϕ (0) = 0, grow outwards from y = 0 as expected from the positiveness of ϕ arising from Eq. (30) with ψ = 0, and approach a constant value for |y| 1. The shape of these solutions is independent of the boundary conditions imposed for V and ψ, but they should only be considered as long as ϕ 0 > ψ 0 .
Numerical solutions for the potential V (y) are provided in Fig. 2. We see in the left panel that for V 0 1 the potential has a positive concavity at the origin but eventually turns over and approaches asymptotically a negative constant. On the other hand, for V 0 5 the potential is always increasing and approaches asymptotically a positive constant. The transition between these two behaviors occurs at V 0 ∼ 3.5. A more detailed view of this transition is provided in the right panel.
The solutions provided in this section are consistent with the required localization of gravity; see Sect. 4 for further details.
3.3 Solution with matter 1: Ansatz for the matter field χ and the matter potential U In this section, we shall obtain a solution for a thick brain in the presence of the matter field χ . In order to close the system of Eqs. (17) to (20) one has to impose two constraints into the system. We choose to set the forms of the potential U (χ ) and the scalar field χ as respectively, where W (χ ) is called the super-potential of χ and takes the form The motivation behind the forms of the mater field χ , the potential U and the super-potential W selected in Eqs. (31) to (33) is their close connection to the standard GR case, which have motivated these forms to be widely used in the literature. Inserting Eqs. (31) to (33) into Eq. (21) provides a differential equation for the warp function A which can be solved immediately to yield a solution of the form where A 0 is a dimensionless constant which is usually set in such a way that the warp function has a value at the origin A (0) = 0. For this purpose, we set A 0 = −1/9. Similarly to the previous case, we assume that the dependency of the potential V in ϕ and ψ is of the form V (ϕ − ψ) in such a way that its derivatives become V ϕ = −V ψ =V . Subtracting Eq. (18) from Eq. (17) and imposing the assumptions mentioned above leads to The equations for the scalar fields given in Eqs. (19) and (  respectively. Note that the previous relationship between V andV given in Eq. (28) still hold, as in this section we have imposed the same constraint on V and definedV in the same way.
The set of Eqs. (35) to (37), along with the relationship provided in Eq. (28) constitute a system of four equations to the four unknowns ϕ, ψ, V andV , and thus the system is determined. Again, we focus on even solutions for ϕ and ψ (and consequently V ). To do so, we impose the same boundary conditions at the origin, i.e., ϕ (0) = ψ (0) = 0, ϕ (0) = ϕ 0 , and ψ (0) = ψ 0 . As a consequence, we have also V (0) = 0. Inserting these boundary conditions into Eq. (35) leads to As we have already mentioned, in order to preserve the positivity of the factor ϕ − ψ we should impose ϕ 0 > ψ 0 . However, to ensure that ϕ (0) > ψ (0), we must also guarantee that ϕ − ψ > 1, thus preserving the positiveness of the factor ϕ − ψ throughout the whole range of y.
Similarly to the case without matter, imposing a boundary condition ψ (0) = 0 in Eq. (37) will consequently set ψ (0) = 0, and thus the only possible even solution for ψ is the constant solution ψ = ψ 0 . Therefore, ψ 0 takes again the role of a free parameter in the problem. The numerical solutions for the scalar field ϕ are plotted in Fig. 3. These solutions satisfy the boundary condition ϕ (0) = 0 and grow outwards from the origin, approaching a constant asymptotic value for |y| 1. The shape of these solutions is independent of the boundary conditions for V and ψ but should only be considered for ϕ 0 − ψ 0 > 1. Finally, numerical solutions for the potential V (y) are given in Fig. 4. Similarly to the case without matter, the potential always presents a positive concavity at the origin but, depending on the boundary conditions V 0 and ϕ 0 , its behavior might eventually turn over and approach negative values (for V 0 5 with ϕ 0 = 10) or grow throughout the whole range of y (for V 0 15 with ϕ 0 = 10). The turning point between the two behaviors occurs at roughly V 0 ∼ 10.
The solutions provided in this section are consistent with the required localization of gravity, as we further describe in Sect. 4.

Solution with matter 2: Ansatz for the warp function A and the scalar field ϕ
To complete the analysis, in this section we will derive another solution for a thick brane in the presence of a matter field χ , thus still having the freedom to impose two constraints to close the system of Eqs. (17) to (20). In this case, we shall leave both the matter field χ and its potential U as where A 0 and ϕ 0 are constants defined positive and k is an arbitrary constant. Following the same reasoning as in previous cases, we assume that the dependency of the potential V in ϕ and ψ is of the form V (ϕ − ψ) in such a way that its derivatives become V ϕ = −V ψ =V . Subtracting Eq. (18) from Eq. (17) and inserting Eqs. (39) and (40) yields Inserting the same Eqs. (39) and (40) into the Eqs. (19) to (21) yields the equations for the scalar fields respectively. Note that the previous relationship between V andV given in Eq. (28) still hold, as in this section we have imposed the same constraint on V and definedV in the same way. The system of Eqs. (41) to (44), along with Eq. (28) is a system of five equations to the five unknowns U , V ,V , ψ, and χ , and thus it is determined. A complete analytic solution for this system can be found through the following reasoning. We will not write explicitly the forms of the solutions due to their lengthy character. Instead we will provide the respective plots.
Focusing on even solutions for the scalar field ψ, we have to impose the boundary condition ψ (0) = 0 into Eq. (43). As a consequence, we obtain as before that ψ (0) = 0 and thus the only possible even solution for ψ is the constant solution ψ = ψ 0 , that we take to be a free parameter of the system. Now, inserting this result into Eq. (41) yields a decoupled differential equation for χ which can be solved immediatly. The solutions of this equation will be real functions if the parameters satisfy the following relation: where ψ c is the critical value. The solution obtained for χ is analytic and can be written in terms of elliptic functions. We provide a plot of this solution for different combinations of parameters in Fig. 5. For values of ψ 0 close to the critical value ψ c given in Eq. (45), the matter field χ has a triple step shape, which becomes less evident with a decrease in ψ 0 . Inserting the previous solution for χ into Eq. (44) allows us to solve for U (χ ), which can then be transformed into U (y) via the chain rule and integrated to obtain a solutions for U (y) as We plot this solution for different combinations of parameters in Fig. 6. The shape of the potential does not change dramatically, but it gets steeper with a decrease in ψ 0 . Finally, Eq. (41) can be solved with respect toV and the corresponding solution can be inserted in Eq. (28) and directly integrated to obtain solutions for the potential V , which are These solutions are plotted in Fig. 7. Unlike the previous cases studied where a change in the parameters of the problem could change the shape of the potential drastically, in this case the general shape of the potential remains the same because we are restricted by the inequality in Eq. (45).
The solutions provided in this section are consistent with the required localization of gravity; see Sect. 4 for further details.

Metric perturbations
In this section, we will consider small perturbations of the metric in the form In addition, we will assume small perturbations on the scalar field χ such that χ → χ +χ(x, y). In this case, one can shown that the (μν) components of the linearized field equations are where (4) = η μν ∂ μ ∂ ν . By using the transverse and traceless gauge (∂ μ H μν = 0 and η μν H μν = H = 0), the metric perturbation decouples from the scalar and the above equation reduces to In order to better understand how this equation relates to stability, we will transform it into a Schrödinger-like equation. First, we make a coordinate transformation dy = e A dz which makes the metric (16) conformally flat. In the new variable z, the Eq. (50) can be written as Now, we remove the term with first order derivative in Eq. (51) by the following redefinition of the tensor perturbation: Thus, Eq. (51) takes the form where we have defined Finally, we perform the decompositionĤ μν (x, z) = ξ μν (x)H (z) to obtain the Klein-Gordon equation (4) ξ μν (x) = m 2 ξ μν (x) and the Schrödinger-like equation where with operators Q and Q † given by This factorization shows that the Schrödinger-like equation cannot support states with negative eigenvalues, i.e., m 2 ≥ 0. Thus, the system is stable against small perturbations of the metric. The zero mode is obtained by performing Q H 0 (z) = 0. The result is where N is a normalization constant. To ensure that 4dimensional gravity can be recovered on the brane, the zero mode must be normalizable. The normalization constant can thus be obtained from For the model without matter from Sect. 3.2, it can be show that the integral in Eq. (60) is finite and the normalization constant N can be computed for different combinations of parameters. It can also be shown that the shape of the potential u (y) given in Eq. (56) is independent of the boundary conditions considered for ϕ, ψ and V , being only affected by changes in the warp function A. Thus, in Fig. 8 we display both the potential u (y) and the zero modes H 0 (y) for different values of k, keeping the remaining free parameters constant. Similarly to what happens in the non-generalized version of the hybrid metric-Palatini gravity, the shape of the potential u is always a double well. Given the finiteness of the integral of the zero mode H 0 , the graviton zero mode can thus be localized on the brane.
For the models featuring a matter field χ in Sect. 3.3 and Sect. 3.4, the integral in Eq. (60) is again finite and the normalization constant N can be computed for different combi- , under an appropriate choice of parameters, this change of shape in the potential implies that the graviton zero mode H 0 is localized in a wider region around the origin, suggesting that the braneworld seems to support internal structure. Given the finiteness of the integral of H 0 over y, one concludes that for both the cases with matter the generalized hybrid metric-Palatini braneworld is stable against fluctuations in the metric.

Summary and discussion
As is well-known, in the standard braneworld scenario, the presence of braneworld configurations with internal structure appeared, in the presence of two scalar fields, in the form of a Bloch brane [56]. This suggests that one could investigate braneworld scenarios within the scalar-tensor representation of the generalized hybrid-Palatini gravity, which consists of two distinct scalar fields, and was in fact explored in the present work. Indeed, an interesting possibility would be to consider two extra fields with standard dynamics, as considered in [56]. Another possibility would be to consider the two extra fields with modified dynamics, as recently considered in [75]. The third model investigated in Ref. [75], in particular, was shown to engender an interesting mechanism to induce internal structure into the brane, so it can also be used in the present context to show its efficiency concerning the induction of internal structure in the brane.
More specifically, in this work, we studied several 5dimensional braneworld scenarios within the context of the generalized hybrid metric-Palatini gravity. We investigated two distinct cases, namely, the presence and the absence of the extra scalar field, χ . In the absence of the field χ , we found an interesting solution, and in the presence of χ , we studied two distinct possibilities. One of the models in the presence of a matter field χ showed that variations of parameters do not add new qualitative effects in the profile of the braneworld solutions. However, the third case studied in Sect. 3.4 showed an interesting effect which was also found in [38], concerning the spreading of the zero mode. Although this new effect is different, it suggests the possibility to study braneworld configurations that engender internal structure.
The presence of internal structure has also appeared before, in the context of modified theories of gravity with non-constant curvature [60,65,67], so this can also be considered in the present context of generalized metric-Palatini gravity. Moreover, is would be of interest to understand how the mechanism used in [62,75] to make the brane asymmetric can be extended to the novel scenario described in the present work. Another issue of current interest concerns the entrapment of fermions and gauge fields inside the brane within the generalized metric-Palatini gravity context explored in the present work. The above issues deserve further investigations, and we are now considering some possibilities, hoping to report on them in the near future. Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and no experimental data has been listed.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .