1. Introduction
We present a formal mathematical argument demonstrating that a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmology with a finite-dimensional parameter space cannot, in general, satisfy the combination of one global early-time integral constraint, one functional late-time constraint, and a set of discrete normalisation constraints. This result is structural and does not rely on any specific observational data.
The overdetermination of FLRW can be diagnosed empirically: several independent probes already hint at this implicitly through their mutually inconsistent preferences for the form of
[
1,
2,
3,
4]. However, none of these observations make the structural conflict explicit. A formal demonstration requires a statement about the dimensionality of the model’s image relative to the space of constraints. For this reason, our argument is framed in terms of a finite-dimensional parameter manifold embedded in an infinite-dimensional Banach space of observables. An empirical analysis cannot establish the impossibility of an early–late FLRW mapping; the geometric argument can.
2. Observational Incompatibility
Empirical results from the Alcock-Paczyński (AP) sector [
5,
6] - in particular the extended tomography AP test applied to SDSS data [
7,
8] - show that using low-redshift geometric information alone, the preferred equation of state satisfies
specifically
in a tension of
with the flat
CDM value
required by the Planck CMB measurement which fixes the early-time integral constraint
together with the acoustic scale
. In a FLRW cosmology, the flat
CDM value forces a specific combination of
and pins the late–time curvature sector to behave as an exact cosmological constant. In particular, the CMB constraints sharply enforce
However, prominent shape-based AP measurements [
8] require a derivative history
whose functional form is incompatible with the
requirement. Matching the observed AP distortion at
demands a deviation
that corresponds precisely to the shift
. If one enforces
, the predicted AP ratio acquires the wrong redshift dependence; if one enforces
, the early–time integral
fails, breaking consistency with the CMB acoustic scale. Because FLRW ties the local derivative constraint
and the global integral constraint
through the same density parameters, no single choice of
reconciles both simultaneously.
Results from the TDCOSMO Collaboration [
4], provide an independent, JWST-calibrated confirmation that this early–late discrepancy persists even under maximally conservative systematics. The authors note that their results cross the
threshold, a consensus-level admission that the contradiction is not an artefact of imprecise data. They further remark that
“there is no obvious leading contender to reconcile the measurements in tension,” an acknowledgement that parameter-level adjustments within
CDM have been exhausted. Consequently, despite refined kinematics, extended samples, and hierarchical modelling, the inferred late-time value of
remains incompatible with the Planck-calibrated early-time geometry.
The recent ACT DR6 results strengthens this situation [
3]. ACT+BAO likelihoods tightly confine the late-time curvature to be indistinguishable from a cosmological constant, sharply restricting any fluid-based freedom in
. Quantitatively, the ACT results enforce
with negligible allowed evolution, because any deviation in the curvature sector alters the CMB–to–BAO distance ladder and is immediately excluded. Thus the entire class of exotic-fluid parametrisations -
models, CPL forms, phantom branches, early or dynamical dark energy, unified dark fluids - lose the degrees of freedom required to adjust the AP shape without simultaneously violating the observed curvature constraint.
In combination, the observational situation is therefore overdetermined under the FLRW mapping.
3. Geometric Origin of the Incompatibility
The observational conflict described above is not a statistical fluctuation but a manifestation of a structural obstruction. Its origin lies in a simple but far-reaching geometric fact: the FLRW model maps a
finite-dimensional parameter space into an
infinite-dimensional constraint space of observables. Once the early–time and late–time quantities are fixed, the model possesses no remaining functional freedom to alter the AP shape. The result is therefore not merely a tension but a failure of the mapping architecture itself. We now make this proposal more precise [
9,
10,
11].
Let
P be the finite-dimensional parameter space of flat FLRW cosmology (e.g.
). The model defines a smooth map
where the constraint space
contains all observational quantities entering CMB, chronometer, and AP analyses. Crucially,
so that
contains at least one function-space component, making it infinite-dimensional.
By contrast, the image is a smooth manifold of dimension at most . Therefore, a tiny finite-dimensional surface is embedded inside a vast infinite-dimensional Banach space of possible observational outcomes.
This is the structural origin of the incompatibility. The CMB integral constraints and the AP derivative constraints do not demand points in ; they demand functions with specific shapes across a redshift interval. The FLRW map, with its finite degrees of freedom, cannot generically reproduce these shapes once early–late sectors are simultaneously imposed. The obstruction is therefore functional, not statistical.
To make this explicit, let
denote the observational regions compatible with Planck and AP respectively. When ACT fixes the late–time curvature sector to a cosmological-constant form, these two regions select
different functional behaviours of the expansion history. Because
has dimension
while these constraints live in a space with infinitely many directions of variation, the image manifold cannot generically intersect both sets simultaneously.
The obstruction is therefore not that the manifold passes “far” from one region while approaching the other. The deeper problem is that
and the functional requirements imposed by and point in different directions in the ambient space. FLRW cannot move in those directions because those directions do not exist within .
Thus no point of
can satisfy both constraint sets simultaneously. Joint-likelihood fits [
8] compute, in effect, the point of the finite-dimensional image manifold closest to two non-intersecting infinite-dimensional constraint regions. The resulting parameter vector is therefore not a consistent cosmology but the projection of incompatible functional requirements onto an insufficiently flexible manifold.
4. Statement of the Incompatibility Theorem
The structural origin of the early–late tension can be expressed in a single geometric statement:
This equality encapsulates the full content of the incompatibility theorem. The symbols appearing in it refer to the following objects.
- (1)
The Constraint Space C.
All possible combinations of cosmological observations are treated as points in the infinite-dimensional Banach space
whose components correspond respectively to (a) an early-time integral constraint such as the CMB acoustic scale
, (b) a continuous Alcock–Paczyński function
defined on a redshift interval
, and (c) a finite collection of late-time distance measurements
. The functional component
ensures that
C is infinite-dimensional.
- (2)
The FLRW Map .
A cosmological parameter vector () is sent by the FLRW model to its predicted observational triple: .
Since P is finite-dimensional (with dimension ), its image is an n-dimensional submanifold of the infinite-dimensional space C. This submanifold represents all observational outcomes compatible with any FLRW cosmology.
- (3)
The Measure .
The symbol
denotes any reasonable measure on
C (Gaussian, quasi-invariant, Wiener, etc.). A key geometric fact is that finite-dimensional submanifolds have measure zero in infinite-dimensional Banach spaces and follows from the fact that Gaussian (or, more generally, quasi-invariant) measures on such spaces assign zero measure to finite-dimensional affine subspaces [
12,
13]. This is the functional-analytic analogue of a line having zero area in a plane.
- (4)
Consequence.
Since is finite-dimensional while C is infinite-dimensional, the FLRW prediction set occupies zero measure inside the space of all possible observational constraints: . Thus the observational outcomes permitted by the FLRW metric form a geometrically negligible subset of constraint space. Equivalently, the chance that a randomly chosen triple lies on the FLRW image manifold is identically zero.
5. Framework, Definitions, and Postulates
We now show formally, using the manifold structure of the FLRW parameter map , that the observational inconsistency reflects a genuine geometric incompatibility: no point of the image manifold can satisfy both constraint regions simultaneously, and the inconsistency cannot therefore be attributable to statistical noise or insufficient data alone.
The functional-analytic setting adopted here traces back to Banach’s foundational formulation of complete normed spaces [
14]. Measure-theoretic properties of infinite-dimensional Banach spaces, including the fact that finite-dimensional submanifolds have Gaussian measure zero, follow the standard modern treatment in Bogachev [
12].
Definition 1 (FLRW Parameter Space). We consider flat FLRW cosmologies specified by a parameter vector . For each , the expansion history is a smooth function for . The space of admissible FLRW histories is thus an n-dimensional manifold of functions.
Definition 2 (Comoving Angular Diameter Distance)
. For a given expansion history , define
This functional depends on the entire inverse expansion history on .
Definition 3 (Alcock–Paczynski Functional)
. Define the AP distortion functional
This observable mixes local information about with integrated information about its past behaviour.
We now state three abstract postulates representing the structure of early- and late-time geometric constraints.
Postulate 1 (Early-Time Integral Constraint)
. There exists a constant such that
where (e.g., the last-scattering epoch). This is a single, global constraint on the definite integral of .
Postulate 2 (Late-Time Functional Constraint)
. There exists a target function defined on a nontrivial interval with such that
In the idealised limit of fine redshift resolution and negligible uncertainty, this imposes an uncountable family of pointwise constraints on .
Postulate 3 (Late-Time Discrete Normalization Constraints)
. There exist redshifts with and constants such that
These impose N additional discrete constraints on the distance scale at low redshift.
6. Theorem: FLRW Structural Incompatibility
Theorem 1 (FLRW Structural Incompatibility)
. Let be a flat FLRW expansion history with . Under Postulates 1–3, the system of constraints
is generically overdetermined.
In particular, unless the triple lies on an n-dimensional submanifold of the infinite-dimensional constraint space , the preimage of this triple under the FLRW map is empty.
Proof. Consider the map
with domain
and codomain
where
denotes the Banach space of continuous functions on
.
The image is (at most) an n-dimensional submanifold of the infinite-dimensional space .
Postulate 1 fixes one component in the factor.
Postulate 2 selects one element in the factor.
Postulate 3 fixes N points in the factor.
Together, these postulates specify a target point
For
to be realized within the FLRW framework, it must lie on the finite-dimensional submanifold
. However, in an infinite-dimensional space, any finite-dimensional submanifold has
measure zero [
12]
Therefore, the FLRW family is generically incapable of satisfying Postulates 1–3 simultaneously. □
6.1. Error-Tolerant Version of the Incompatibility
In practice, observations do not specify a single constraint triple , but an allowed neighbourhood determined by measurement uncertainties. We now show that the measure-zero obstruction persists even when we allow a finite tolerance.
Let
be the separable Banach constraint space as above, endowed with a non-degenerate Gaussian measure
(for example, the product of a Gaussian measure on the function-space component
and Euclidean Gaussian measures on the finite-dimensional factors) [
12]. Let
be the FLRW parameter map, where
P is a finite-dimensional smooth manifold.
For
define the
–thickening of the image
Theorem 2 (Error-tolerant incompatibility)
. Assume is with . Then for every ,
In particular, for –almost every constraint triple and every fixed , there is no parameter vector for which the FLRW predictions lie within the error tolerance:
Sketch of proof. Since
P is finite-dimensional and
is
, the image
is contained in a countable union of finite-dimensional
embedded submanifolds of
(by taking a countable atlas of
P and using the separability of
). For a Gaussian measure on an infinite-dimensional separable Banach space, any finite-dimensional
submanifold has measure zero [
12].
Furthermore, the –tubular neighbourhood of such a submanifold can be covered by a countable union of Lipschitz images of for some finite k (using a partition of unity and local parametrizations of the manifold and its normal bundle). Since Lipschitz images of are also –null (as they are still finite-dimensional in measure-theoretic sense under the Gaussian measure), the thickening remains a –null set.
Now let
By definition,
, so
Hence, for
–almost every
, no choice of parameters
p achieves
□
Operationally, Theorem 2 states that inflating the observational target from a single triple to any fixed-width error region (a ball of radius in ) does not rescue the FLRW map: the set of observational triples that can be matched even within the uncertainties still has measure zero in the full constraint space. Measurement errors enlarge the target, but a finite-dimensional cosmology remains too small to hit it except on a –null set.
7. Remarks
The incompatibility already follows from Postulates 1 and 2. Postulate 3 further restricts the compatible set, but is not required for the structural overdetermination.
Finite Data Sampling The conclusion of overdetermination persists if Postulate 2 is weakened to require agreement only at a finite set of M redshifts. In this case the constraint space reduces to , while the FLRW parameter space remains n-dimensional. Whenever - the situation for all modern surveys - the system is overdetermined even without invoking the idealised functional limit.
Fine-Tuning Interpretation The theorem implies that a simultaneous solution to Postulates 1-3 requires the target triple to lie on the n-dimensional submanifold within an infinite-dimensional constraint space. Such alignment constitutes a fine-tuning condition of co-dimension . In the absence of this special tuning, the system has no solution within the FLRW family. This clarifies that the structural incompatibility is not a statistical anomaly but a generic feature of the FLRW mapping.
On the status of and Dark Energy Although represents an unobserved physical component with negative pressure - dark energy, it contributes only a single scalar parameter to the cosmological map. Its inclusion therefore, leaves the image of the FLRW model finite-dimensional, and thus of measure zero. In particular, may alter the dynamics of the expansion but does not introduce a functional degree of freedom capable of relaxing the rigid early–late coupling responsible for the incompatibility.
Interpretation The theorem does not assert that FLRW fails for particular data, nor that current tensions arise from statistical fluctuations. It states that the rigidity of the current FLRW map - a finite-dimensional parameterisation attempting to satisfy both an early-time integral constraint and a late-time functional derivative constraint - makes compatibility a non-generic event of measure zero. The incompatibility is geometric, not statistical. (
See Appendix A & Appendix B for further clarifications)
8. Outlook
Formulating the early–late tension within a functional-analytic framework clarifies both the origin of the incompatibility and suggests a route to its resolution. The obstruction proved here follows from a single geometric fact: the FLRW model defines a smooth map from a finite-dimensional parameter manifold P into the infinite-dimensional constraint space . The image is therefore a finite-dimensional submanifold of and, for generic observational data , has no preimage under . This places the tension on structural rather than statistical grounds.
A crucial implication is that adding further FLRW parameters cannot resolve the incompatibility: such parameters merely reposition points within the same finite-dimensional image , because they alter the numerical inputs to the map without changing its functional form. By contrast, resolving the obstruction requires structural degrees of freedom - that modify the map itself rather than its coordinates.
Accordingly, restoring compatibility requires enlarging the image of by extending the function space, and defining a correspondingly deformed map whose additional elements genuinely alter the functional form of the mapping while still arising from a single coherent expansion history.
A structural function or functions may be introduced phenomenologically at first as effective geometric degrees of freedom that relax the classical rigidity of FLRW while preserving its early-universe limit. Their physical interpretation - whether geometric, dynamical, or topological - may emerge only once a consistent structural extension has been identified.
Functional analysis not only diagnoses the incompatibility but also supplies the geometric constraint that any viable extended FLRW theory must satisfy: it must introduce, at minimum, additional structural degrees of freedom required to overcome the mapping obstruction established here.
9. Conclusion
We have established a purely geometric obstruction within the FLRW framework: the model imposes a rigid coupling between the global integral of and its derivative manifestation at low redshift through the functional . We have demonstrated that the FLRW metric, the geometric foundation of the standard cosmological model, imposes a mathematical rigidity that is incompatible with the existence of three independent classes of geometric observations: the CMB acoustic scale (fixing the high-z integral), tomographic AP measurements (fixing the low-z derivative), and late-time distance indicators (fixing the amplitude). No function H(z) within this framework can satisfy all three constraints simultaneously. The incompatibility follows directly from the finite dimensionality of the FLRW parameter space relative to the infinite-dimensional constraint space of observables. A finite-dimensional image cannot, except on a set of measure zero, satisfy independent functional constraints drawn from distinct observational sectors. It follows that any consistent reformulation of the cosmological mapping must introduce at a minimum a functional degree of freedom to enlarge the image. This requirement should be regarded as a structural condition on any future extension of the standard model.
Our result formalises the current emerging anomalies in early-late FLRW mapping. It explains why no amount of greater precision in future surveys within an unextended FLRW framework will resolve the tension.
Appendix A. Clarifications on the Structural Incompatibility Theorem
The following notes address common questions or misconceptions that may arise in response to Theorem 1. Each item provides a concise clarification of the mathematical structure of the argument.
-
Q1.
-
Is the theorem trivial? Is it merely the statement that a model with few parameters becomes overdetermined by enough data?
A. No. The theorem concerns the geometry of the map , which includes both a global integral functional and a continuous functional. The structural incompatibility arises from the fact that the FLRW image is an n-dimensional submanifold of an infinite-dimensional space. This is not a counting of data points but a geometric limitation of the FLRW mapping itself.
-
Q2.
-
What if the AP functional is effectively low-dimensional? Could it reduce to a few principal components?
A. Even in the finite-sampling case, the constraint space becomes . Overdetermination arises whenever , which is the case for all modern datasets. Thus the structural mismatch persists without taking the idealised limit of a full functional constraint.
-
Q3.
-
Are the constraints really independent? Could the early-time integral and late-time AP shape be related?
A. They are independent by construction. The early-time constraint fixes a definite integral of up to , while depends only on H and its integral up to . No functional identity links them for general FLRW histories; any such identity would have to be explicitly demonstrated.
-
Q4.
-
Could FLRW be enlarged by adding curvature or more fluids? Do more parameters remove the tension?
A. Adding curvature increases n by at most one. Adding smooth dark-energy fluids enlarges but does not alter the fact that FLRW backgrounds constitute a finite-dimensional family of functions. The theorem applies to the entire class of finite-parameter FLRW expansions. Avoiding it requires modifying the metric ansatz itself.
-
Q5.
-
What about modified gravity? Do or Horndeski theories increase the freedom of the background?
A. No. In any metric theory with FLRW symmetry, the background is still determined by a finite number of parameters. Modified gravity does not generate an infinite-dimensional family of background expansions, so the theorem applies unchanged.
-
Q6.
-
Does the theorem assume infinite precision? Is the functional constraint valid only in the zero-uncertainty limit?
A. No. The finite-sampling form of the AP constraint is treated explicitly in Remark 2. The structural incompatibility arises whenever the number of independent low-redshift constraints exceeds the dimension n of the FLRW parameter space. The idealised functional limit is used only to describe the geometry of the constraint space.
-
Q7.
-
Could the universe lie exactly on the FLRW-compatible manifold? Is fine-tuning a valid escape?
A. This requires the target triple to lie on an n-dimensional submanifold of the infinite-dimensional space . Such points form a set of measure zero and correspond to extreme fine-tuning. The theorem therefore, treats compatibility as non-generic.
A natural concern is whether a sufficiently tuned choice of cosmological parameters might restore consistency within the FLRW framework. This intuition arises from the usual statistical practice of adjusting the small number of FLRW parameters to accommodate diverse datasets. However, the incompatibility demonstrated in the preceding formal proof is not parametric but geometric. Once the CMB calibration fixes the global integral constraint, the FLRW ansatz determines the full late-time functional form of . The AP sector imposes an effectively infinite-dimensional constraint, and the image of a finite-dimensional parameter manifold cannot, except on a set of measure zero, intersect it. Consequently, no amount of parameter tuning can reconcile the two constraint classes; the obstruction arises from the rigidity of the FLRW mapping itself.
-
Q8.
-
FLRW is successful empirically. Does that contradict the theorem?
A. No. FLRW can successfully describe individual observational sectors (e.g. early-time geometry or specific late-time probes). The theorem concerns the joint imposition of early-time integral constraints and late-time derivative constraints. It states that the FLRW map cannot generically satisfy both simultaneously.
-
Q9.
-
What if the FLRW map has hidden degeneracies or unusual Jacobians? Could this enlarge its image?
A. No. A map from an n-dimensional manifold has image of dimension at most n, even if the map is non-injective. Non-injectivity cannot increase the dimensionality of the image. This is a standard result in differential geometry and ensures the theorem’s robustness.
-
Q10.
-
Isn’t this just saying FLRW is rigid? Is the result philosophical rather than technical?
A. The rigidity is a mathematical property of the FLRW line element. The theorem formally demonstrates that this rigidity has a concrete consequence: the combined early-time, functional late-time, and discrete normalisation constraints generically have no common preimage in the finite-dimensional FLRW family. This is a technical geometric result, not a philosophical claim.
Appendix B. Conceptual Clarification of the Geometric Setting
Our argument employs geometric language that lies slightly outside the standard cosmology toolkit, which typically focuses on parametric fits, Fisher matrices, and likelihood surfaces. In that familiar context, a manifold is often visualized as a flexible two or three-dimensional surface embedded in . This appendix clarifies a key conceptual point that can arise when the argument is first encountered: why the FLRW image manifold and an observed constraint point cannot, in general, be made to intersect.
Dimensionality and Measure
The image of the FLRW map has dimension equal to the number of cosmological parameters (at most –7 in typical models). The constraint space is infinite-dimensional because of its function-space component .
In an infinite-dimensional Banach space, any finite-dimensional submanifold has
measure zero. This is not a visual metaphor - it is a rigorous statement in geometry [
12,
13]. Therefore a randomly chosen point in such a space lies on a given finite-dimensional submanifold with probability zero.
Why the “Flexible Sheet” Intuition Fails
The intuitive picture of two finite-dimensional surfaces in that can be bent or deformed to intersect does not apply here. The actual setting is fundamentally different:
A finite-dimensional manifold (the image of FLRW parameters),
embedded in an infinite-dimensional Banach space (the space of constraint triples ).
The embedding map is fixed by the physics of the FLRW metric and distance integrals; it is not a freely deformable geometric object. The impossibility of intersection for a generic constraint point is not assumed but forced by the structure of the problem:
The finite dimensionality of ,
The functional (infinite-dimensional) nature of the AP constraint,
The fixed acoustic scale calibration (the early-time integral),
The lack of additional functional degrees of freedom in the FLRW ansatz.
This clarifies why our theorem is structural rather than empirical. The mismatch arises not from noisy data or insufficient model complexity, but from the intrinsic geometric limitations of the FLRW framework when confronted with independent early and late-universe constraints. Consequently, when modern surveys provide high-dimensional (or functional) late-time constraints, the probability that these align exactly with a Planck-calibrated FLRW model is zero - not small, but mathematically zero. This is the essence of the Redshift Catastrophe.
References
- Collaboration, D. BAO constraints from baryon acoustic oscillations in galaxy, quasar and Ly-α forest tracers: DESI Data Release 1. Journal of Cosmology and Astroparticle Physics 2024, 2024, P000.
- Scolnic, D.; et al. Pantheon+ Type Ia Supernova Distance Measurements. Zenodo, Dataset,urlhttps://github.com/PantheonPlusSH0ES/DataRelease, 2025.
- Calabrese, E.; Hill, J.C.; Jense, H.T.; La Posta, A.; Abril-Cabezas, I.; Addison, G.E.; Ade, P.A.; Aiola, S.; Alford, T.; Alonso, D.; et al. The Atacama Cosmology Telescope: DR6 constraints on extended cosmological models. Journal of Cosmology and Astroparticle Physics 2025, 2025, 063. [CrossRef]
- Birrer, S.; Buckley-Geer, E.J.; Cappellari, M.; Courbin, F.; Dux, F.; Fassnacht, C.D.; Frieman, J.A.; Galan, A.; Gilman, D.; Huang, X.Y.; et al. TDCOSMO 2025: Cosmological constraints from strong lensing time delays. Astron. Astrophys. 2025, 704, A63. [CrossRef]
- Alcock, C.; Paczynski, B. An evolution free test for non-zero cosmological constant. Nature 1979, 281, 358–359. [CrossRef]
- Ballinger, W.E.; Peacock, J.A.; Heavens, A.F. Measuring the cosmological constant with redshift surveys. Monthly Notices of the Royal Astronomical Society 1996, 282, 877–888. [CrossRef]
- Dong, F.; Park, C.; Hong, S.E.; Kim, J.; Hwang, H.S.; Park, H.; Appleby, S. Tomographic Alcock–Paczyński Test with Redshift-space Correlation Function: Evidence for the Dark Energy Equation-of-state Parameter w > −1. The Astrophysical Journal 2023, 953, 98. [CrossRef]
- Dong, X.; et al. Constraints on dynamical dark energy from the extended Alcock–Paczynski test with SDSS data. Astrophysical Journal 2025. in press.
- Lang, S. Fundamentals of Differential Geometry; Springer New York, 1999. [CrossRef]
- Hamilton, R.S. The inverse function theorem of Nash and Moser. Bulletin of the American Mathematical Society 1982, 7, 65–222. [CrossRef]
- Abraham, R.; Marsden, J.E.; Ratiu, T. Manifolds, Tensor Analysis, and Applications, 2nd ed.; Vol. 75, Applied Mathematical Sciences, Springer, 1988. [CrossRef]
- Bogachev, V.I. Measure Theory; Springer Berlin Heidelberg, 2007. [CrossRef]
- Kuo, H.H., Gaussian measures in Banach spaces. In Gaussian Measures in Banach Spaces; Springer Berlin Heidelberg, 1975; p. 1–109. [CrossRef]
- Banach, S. Théorie des Opérations Linéaires; Monografje Matematyczne: Warsaw, 1932. Reprinted by Chelsea Publishing Company, New York.
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).