# Non-Oscillatory Power Spectrum From States of Low Energy in Kinetically Dominated Early Universes

^{1}Departamento de Física Teórica and IPARCOS, Universidad Complutense de Madrid, Madrid, Spain^{2}Departamento de Física Teórica y Del Cosmos, Universidad de Granada, Granada, Spain

Recently, States of Low Energy (SLEs) have been proposed as viable vacuum states of primordial perturbations within Loop Quantum Cosmology (LQC). In this work we investigate the effect of the high curvature region of LQC on the definition of SLEs. Shifting the support of the test function that defines them away from this regime results in primordial power spectra of perturbations closer to those of the so-called Non-oscillatory (NO) vacuum, which is another viable choice of initial conditions previously introduced in the LQC context. Furthermore, through a comparison with the Hadamard-like SLEs, we prove that the NO vacuum is of Hadamard type as well.

## 1 Introduction

In a previous work (Martín-Benito et al., 2021), we have proposed the States of Low Energy (SLEs) introduced in (Olbermann, 2007) as viable candidates for the vacuum state of cosmological perturbations in Loop Quantum Cosmology (LQC). We were motivated by the fact that they were proven to be Hadamard states that minimized the regularized energy density when smeared along the time-like curve of an isotropic observer via a test function. Furthermore, they had been shown to provide a qualitative behavior in the ultraviolet (UV) and infrared regimes of the primordial power spectra of scalar and tensor perturbations that agrees with observations in models where a period of kinetic dominance precedes inflation (Banerjee and Niedermaier, 2020), which is the case in LQC. However, in (Martín-Benito et al., 2021) we have only considered test functions that could be seen as natural choices within LQC, namely, ones with support on the high curvature regime. As long as this is the case, we have shown that the ambiguity in the introduction of the test function is surpassed in this context, in the sense that the resulting SLE and power spectra seem to be very insensitive to its shape and support, provided it is wide enough.

In this work, we investigate the effect of shifting the test function away from the high curvature regime. Firstly, this provides a more complete analysis of the SLEs and the ambiguity of the test function. Secondly, this allows us to distinguish in the primordial power spectra the consequences coming directly from LQC corrections and those related to having a period of kinetic dominance prior to inflation, which can also be obtained in a classical scenario. We will show that if the test function ignores the Planckian region, the effect in the resulting SLE is appreciable. Furthermore, in the power spectra the oscillations that were previously found for lower wave numbers are now dampened.

This motivates us to compare our results with those found in the LQC literature that adopts as initial conditions for the perturbations the so-called non-oscillatory (NO) vacuum state (de Blas and Olmedo, 2016; Castelló Gomar et al., 2017). As the name suggests, this state is precisely defined to minimize mode by mode the amplitude of the oscillations of the primordial power spectra in a given time interval. It turns out that this minimization in time is reflected in a minimization of oscillations in the *k* domain of the power spectra. This NO prescription has been motivated as well as a good candidate for the vacuum of the perturbations (Elizaga Navascués et al., 2020). One question that so far remained unanswered is whether this NO vacuum is or not of Hadamard type. In this work, by comparing with SLEs, we show that indeed this is the case. To do so, we resort to their UV expansions, obtained for SLEs in (Banerjee and Niedermaier, 2020) and for the NO vacuum in (Elizaga Navascués et al., 2020).

This manuscript is organized as follows. In Section 3 we review the application of SLEs in LQC as presented in (Martín-Benito et al., 2021). In Section 4 we explore the consequence of excluding the high curvature regime from the test function, computing the corresponding SLEs and power spectra at the end of inflation. Section 5 is devoted to a proof that the NO vacuum is Hadamard, based on the comparison with SLEs in the UV limit. Finally, we conclude in Section 6 with a discussion and closing remarks.

Throughout we adopt Planck units *c* = *ℏ* = *G* = 1 for numerical computations, keeping factors of *G* in expressions.

## 2 Cosmological Perturbations and States of Low Energy in LQC

In this section we will briefly review the dynamics of cosmological perturbations in LQC through its hybrid approach, as well as the definition of SLEs in this context, as exposed in (Martín-Benito et al., 2021). Let us start by considering the spatially flat FLRW model with scale factor *a*, minimally coupled to the scalar field *ϕ* subject to the potential *V*(*ϕ*), which will drive inflation. Cosmological perturbations are usually described by scalar and tensor gauge invariant perturbations *I* denotes the two possible polarizations of tensor perturbations. Expanding in Fourier modes

where the dot denotes derivative with respect to cosmological time *t*, and *ω*_{k} depends on the quantization. It is common to work with the rescaled fields *η*, such that *dη* = *dt*/*a*. Then, we find the equations of motion of the Fourier modes of these fields, *u*_{k} and *μ*_{k} respectively, to be.

where the prime denotes derivative with respect to conformal time *η* and *s*^{(s)} (*η*) and *s*^{(t)} (*η*) are the time-dependent mass terms of scalar and tensor modes respectively. From the hybrid approach to LQC, one can write these as functions of the background variables *a*, *ρ* (inflaton energy density), *P* (inflaton pressure) and the inflaton potential *V*(*ϕ*) as (Elizaga Navascués et al., 2018a):

where

To simplify notation, in the following we will use *s*(*η*) to refer generically to both of them, as our comments apply equally to both scalar and tensor modes. When doing so, for simplicity, we will refer only to *u* as everything is analogous for *μ*^{I}. It is easy to find that *s*(*η*) can be related to

Generally, there are no analytical solutions to such equations of motion, and results have to be obtained numerically, given initial conditions *u*_{k} (0),

where *D*_{k} is a positive function and *C*_{k} any real function. Once defined, the perturbations can be evolved until a time *η*_{end} during inflation when all the scales of interest have crossed the horizon. The primordial power spectra of the comoving curvature perturbation

are evaluated at *η* = *η*_{end}. The choice of initial conditions amounts to a choice of vacuum state for the perturbations. In this context, there is no notion of a unique natural vacuum. Indeed, several proposals have been made of initial vacua within the LQC framework (Agullo et al., 2015; de Blas and Olmedo, 2016; Ashtekar and Gupt, 2017; Elizaga Navascués et al., 2020) that result in primordial power spectra compatible with observations. In these analyses, initial conditions are set at the LQC bounce where the scale factor of the geometry reaches a minimum, and then it starts expanding. At this bounce, the spacetime curvature reaches a maximum value of the order of the Planck scale. The work of (Martín-Benito et al., 2021) applied the SLE construction defined in (Olbermann, 2007) to this context. These are defined as the states that minimize the energy density smeared along a time-like curve, specified by a test function *f*. In the following we summarize this procedure, adapted to our notation (namely working with *u* and *μ* and in conformal time). For further details we refer the reader to (Olbermann, 2007; Martín-Benito et al., 2021). Given a fiducial solution *v* to the equation of motion (3), the SLE associated to the test function *f*(*η*) is found through the Bogoliubov transformation

where the Bogoliubov coefficients *α*(*k*) and *β*(*k*) are found uniquely (up to a phase) to be.

with.

Note that these quantities carry a dependence on the test function *f*. Indeed, as remarked, Eq. 10 defines the SLE associated to this *f*. This introduces an ambiguity in the procedure, which has been explored within the LQC approach in (Martín-Benito et al., 2021). In that work, only natural choices for *f* within this framework were considered, whose support thus included the bounce of LQC. In this current investigation, we will consider test functions that exclude it.

## 3 Effect of the Bounce in Sles

In (Martín-Benito et al., 2021) we have shown that there are two families of test functions that can be seen as natural choices for the smearing function within LQC, and that provide SLEs that are very insensitive to their particular form. Namely, we have found that for a test function supported around the bounce of LQC the resulting SLE does not qualitatively depend on its shape or support, as long as it is wide enough. In the case of a test function supported on the expanding branch only, from the bounce onward, in (Martín-Benito et al., 2021) only the case of a steep (but smooth) step function was investigated, in order to fully retain the contributions coming from the bounce. In this case, the SLE remains insensitive to the size of the support as long as it is wide enough. The resulting power spectrum inherits this independence on the choice of test function, and coincidentally shows good agreement with the one of a second order adiabatic vacuum state.

In this section we explore the consequences of not including the bounce in the support of the test function. This way we will be able to study also the effect of the shape of the test function when supported only on the expanding branch away from the high-curvature regime. This will allow us to provide a comparison with an analogous classical scenario of an FLRW model with a period of kinetic dominance prior to inflation.

Let us start by considering the smooth step function *f*^{2} plotted in Figure 1, supported in the interval

where *δ* determines the ramping up, with a smaller *δ* resulting in a steeper step, and *S* is the auxiliary function:

**FIGURE 1**. The smooth step function defined in (15), represented in terms of its parameters: initial and final points, *η*_{i} and *η*_{f} respectively, and *δ*, which controls the ramping up.

Figure 2 shows the initial conditions, parametrized through (8), corresponding to the SLE obtained for scalar perturbations when considering the test function (15), with *η*_{0} = 0, 1, 10 and 100 Planck seconds after the bounce, with *η*_{f} fixed at the onset of inflation, and for a sharp step of *δ* ∼ 0.06. The case of *η*_{0} = 0 corresponds to the one analysed in (Martín-Benito et al., 2021). The effect of excluding the bounce is immediately noticed as soon as the support of the test function is moved one Planck second into the expanding branch. If we push the initial time further into the future, the change is gradually decreased, and for *η*_{0} = 100 we see some convergence. The corresponding figure for tensor modes is omitted since the initial conditions are essentially the same, as discussed in (Martín-Benito et al., 2021). Within this family of test functions that exclude the bounce, we have also investigated the consequences of changing their shape. In all these cases, we find that, as the starting point moves further away from the bounce, the SLE becomes more insensitive to the shape of the test function. For this reason, below, we will focus our comments on the four step functions defined above, as they already show the different qualitative behaviors one may obtain from different test functions in this scenario.

**FIGURE 2**. Initial conditions in terms of *D*_{k} and *C*_{k}, as constructed in (8), for scalar modes at the bounce corresponding to the SLEs obtained with the window function (Eq. 15) covering the expanding branch until the onset of inflation with starting points: *η*_{0} = 0 (solid gray line), *η*_{0} = 1 (dashed red line), *η*_{0} = 10 (dotted green line) and *η*_{0} = 100 (dotted-dashed blue line). The scale of *k* is in Planck units. All computations were performed for a quadratic potential *V*(*ϕ*) = *m*^{2}*ϕ*^{2}/2, with *m* = 1.2 × 10^{–6} and with the value of the scalar field at the bounce fixed to *ϕ*_{B} = 1.225 (toy value). For tensor modes, the resulting SLE at the bounce shows no significant qualitative differences.

Figure 3 shows the corresponding primordial power spectra for scalar and tensor perturbations, computed through (9). Here, the effect of removing the bounce is evident. As the support of the test function is pushed further away from the high curvature regime, the oscillations in the power spectra are gradually dampened. It is interesting to note that, in fact, as Figure 4 shows, the power spectra are pushed towards those obtained from the non-oscillatory (NO) vacuum state defined in (de Blas and Olmedo, 2016), which is constructed by minimizing the oscillations in time of the power spectrum of perturbations for the whole expanding branch, including the bounce. We further note that the case where the support of the test function starts at *η*_{0} = 100 will essentially correspond to that obtained by using the SLE as the vacuum state of primordial perturbations in a classical FLRW model with a period of kinetic dominance prior to inflation. However, for smaller *η*_{0}, SLEs show oscillations in *k* at and below scales comparable to those of the curvature at that initial time. Then we can conclude that the oscillations that appear in the power spectra when including the high curvature region (for instance the bounce of LQC) in the support of the test function open an interesting observational window.

**FIGURE 3**. Power spectra of the comoving curvature perturbation *η*_{0} = 0 (solid gray line), *η*_{0} = 1 (dashed red line), *η*_{0} = 10 (dotted green line) and *η*_{0} = 100 (dotted-dashed blue line). The scale of *k* is in Planck units. All computations were performed for a quadratic potential *V*(*ϕ*) = *m*^{2}*ϕ*^{2}/2, with *m* = 1.2 × 10^{–6} and with the value of the scalar field at the bounce fixed to *ϕ*_{B} = 1.225 (toy value).

**FIGURE 4**. Comparison between the power spectra of the comoving curvature perturbation *η*_{0} = 100 (dotted-dashed blue line) and to the NO vacuum (solid black line). The scale of *k* is in Planck units. All computations were performed for a quadratic potential *V*(*ϕ*) = *m*^{2}*ϕ*^{2}/2, with *m* = 1.2 × 10^{−6} and with the value of the scalar field at the bounce fixed to *ϕ*_{B} = 1.225 (toy value).

For completion, we added an appendix where we apply the SLE and NO vacuum prescriptions in a classical Universe dominated by the kinetic energy of the scalar field. We discuss the situations in which they agree with the natural choice for vacuum state considered in (Contaldi et al., 2003).

## 4 Comparison Between Sle and No Vacuum

One remarkable property of the SLEs that is usually not explicitly proven for other vacua proposals is that they are of Hadamard type (Olbermann, 2007; Banerjee and Niedermaier, 2020). This guarantees that computations such as that of the expectation value of the renormalized stress-energy tensor will be well defined. On the other hand, the NO vacuum has only been proven to behave in the ultraviolet (UV) asymptotic regime as a high order adiabatic state, at least of fourth order (Elizaga Navascués et al., 2018b). Indeed, considering two adiabatic states of orders *n* and *m*, one can compute the *β* coefficient of the Bogoliubov transformation between the two:

and find that in the UV |*β*| decays with *k*^{−l−2}, where *l* = min (*n*, *m*). In the case of the comparison of the NO vacuum with an nth-order adiabatic one, it was found that |*β*| ∼ *k*^{−2−n} at least up to *n* = 4, which shows that the NO vacuum is the highest order one of the two. As a Hadamard type vacuum is an infinite order adiabatic state, this is an indication that the NO vacuum might be as well, though a stronger proof would be desirable. In this section, we will provide one, through a comparison with the (Hadamard-like) SLEs.

To simplify the comparison, let us write the UV expansions of both the SLE and the NO state as

where ∼ means the behavior in the large *k* regime. The NO vacuum state has recently been analysed analytically in (Elizaga Navascués et al., 2020). In particular, that work has found that the state admits the UV asymptotic expansion (18) with:

where

and the *γ*_{n} coefficients are given by the iterative relation

with *γ*_{0} = *s*(*η*) and *γ*_{−n} = 0 for all *n* > 0. With this expansion, we are able to compute the NO state up to any order in 1/*k* easily. Actually, one can check by direct inspection that

where the *G*_{n} are determined recursively by

with *G*_{0} = 1. Remarkably, in (Banerjee and Niedermaier, 2020), the SLEs are found to have the same asymptotic expansion (22), regardless of the choice of the test function.

Therefore, the *β* coefficients of the Bogoliubov transformation between the SLE and NO vacuum are identically zero in the UV. Thus, we conclude that the NO vacuum is of Hadamard type since it displays exactly the same short-distance structure as the SLEs.

## 5 Conclusion and Discussion

SLEs have recently been proposed as a suitable choice for the vacuum state of perturbations in LQC (Martín-Benito et al., 2021), where the dependence of the state on the test function was explored. For that investigation, only test functions that included the bounce of LQC were analysed, as they are natural choices within this framework. In this work, we investigate the effect of pushing the support of the test function away from the bounce and indeed from the high curvature regime. In addition to offering a more complete analysis of SLEs within LQC, this allows us to disentangle the effects coming from quantum corrections to the dynamics, which are important in the high curvature regime, from those that arise from having only a period of kinetic dominance prior to inflation, which can be found also in classical inflationary models, and is not a direct consequence of the quantum nature of geometry.

We have found that whether the support of the test function includes the high curvature regime or not has a greater influence on the resulting SLE than any other parameter of the test function that has been studied previously. Indeed, in (Martín-Benito et al., 2021), we had already shown that as long as the support of the test function includes the high curvature regime and it is wide enough, the SLE is very insensitive to its shape and support. In this work, we have shown that as soon as the test function is pushed away from the bounce, the SLE suffers a big shift, which then converges as the test function is pushed further away from the high curvature regime. We have also found that in this case, when convergence with respect to the support is reached, the SLE is again insensitive to the shape of the test function. Furthermore, through the computation of the power spectra of perturbations at the end of inflation, we see that as the test function is shifted away from the high curvature region the oscillations found for lower modes (scales comparable to those of the curvature at initial time) are gradually dampened, and the spectra are pushed to those of the NO vacuum state introduced in (de Blas and Olmedo, 2016). Then it is safe to conclude that these oscillations are in fact a consequence of the corrections coming from LQC, which opens an interesting observational window into signatures from LQC in observations of the CMB. For instance, the enhancement of power at super-Hubble scales in the primordial power spectrum of scalar perturbations has been proposed, together with large scale non-Gaussianities, as a mechanism to explain several anomalies in the CMB (Agullo et al., 2021a; Agullo et al., 2021b). From this perspective, the power spectrum provided by SLEs prescription when including the bounce is physically relevant. On the other hand, the NO-like power spectra show a lack of power at large scales of primordial origin that can, on the one hand, alleviate some tensions in the CMB (Ashtekar et al., 2020; Ashtekar et al., 2021), and on the other hand, ease the *trans*-Planckian issues on these scenarios (Brandenberger and Martin, 2013; Ashtekar and Gupt, 2017). See (Li et al., 2021) for a recent comparison on further details about different proposals in LQC. However, a detailed analysis of all this requires a rigorous investigation that we leave for future work.

Finally, the fact that SLEs are proven to be Hadamard is a great advantage that most proposals don’t enjoy. Typically, this property is difficult to prove explicitly. One strategy, that may be enough for practical purposes, is to compare a state with an adiabatic one of increasing (finite) order, and show, through the *β* coefficients of the Bogoliubov transformation between the two states, that the state in question is always of higher order than the adiabatic one considered. This shows that it is at least a very high order adiabatic state, and since a Hadamard state is an adiabatic state of infinite order, then most likely so is the proposed state. However, we now have a family of states, namely SLEs, that are explicitly Hadamard. Therefore, the *β* coefficients of the transformation between any Hadamard state and any SLE should decrease faster than any power of the wave number. We have applied this reasoning to the NO vacuum state, that had previously been shown to be at least of fourth order (Elizaga Navascués et al., 2018b). We find that, in the ultraviolet limit of large wave numbers, the asymptotic expansion that the NO vacuum satisfies (found in (Elizaga Navascués et al., 2020)) agrees exactly with that of the SLE (Banerjee and Niedermaier, 2020) (no matter the test function chosen to define it). As a consequence, the *β* coefficients of the transformation between the two will be identically zero in the ultraviolet. In other words, the NO vacuum has the same short-distance structure than the SLEs, which proves that the NO vacuum is of Hadamard type as well.

## Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

## Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

## Funding

This work is supported by the Spanish Government through the projects FIS2017-86497-C2-2-P and PID2019-105943GB-I00 (with FEDER contribution). RN acknowledges financial support from Fundação para a Ciência e a Tecnologia (FCT) through the research grant SFRH/BD/143525/2019. JO acknowledges the Operative Program FEDER2014-2020 Junta de Andalucía-Consejería de Economía y Conocimiento under project E-FQM-262-UGR18 by Universidad de Granada.

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Acknowledgments

The authors would like to thank Beatriz Elizaga Navascués and Guillermo Mena Marugán for useful discussions.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2021.702543/full#supplementary-material

## References

Agullo, I., Nelson, W., and Ashtekar, A. (2015). Preferred Instantaneous Vacuum for Linear Scalar fields in Cosmological Space-Times. *Phys. Rev. D* 91, 064051. doi:10.1103/PhysRevD.91.064051

Agullo, I., Kranas, D., and Sreenath, V. (2021). Anomalies in the CMB from a Cosmic Bounce. *Gen. Relativ Gravit.* 53, 17. doi:10.1007/s10714-020-02778-9

Agullo, I., Kranas, D., and Sreenath, V. (2021). Large Scale Anomalies in the CMB and Non-gaussianity in Bouncing Cosmologies. *Class. Quan. Grav.* 38, 065010. doi:10.1088/1361-6382/abc521

Ashtekar, A., Gupt, B., and Sreenath, V. (2021). *Cosmic Tango between the Very Small and the Very Large: Addressing CMB Anomalies through Loop Quantum Cosmology*.

Ashtekar, A., and Gupt, B. (2017). Initial Conditions for Cosmological Perturbations. *Class. Quan. Grav.* 34, 035004. doi:10.1088/1361-6382/aa52d4

Ashtekar, A., Gupt, B., Jeong, D., and Sreenath, V. (2020). Alleviating the Tension in the Cosmic Microwave Background Using Planck-Scale Physics. *Phys. Rev. Lett.* 125, 051302. doi:10.1103/PhysRevLett.125.051302

Banerjee, R., and Niedermaier, M. (2020). Bonus Properties of States of Low Energy. *J. Math. Phys.* 61, 103511. doi:10.1063/5.0019311

Brandenberger, R. H., and Martin, J. (2013). Trans-Planckian Issues for Inflationary Cosmology. *Class. Quan. Grav.* 30, 113001. doi:10.1088/0264-9381/30/11/113001

Castelló Gomar, L., Mena Marugán, G. A., Martín de Blas, D., and Olmedo, J. (2017). Hybrid Loop Quantum Cosmology and Predictions for the Cosmic Microwave Background. *Phys. Rev. D* 96, 103528. doi:10.1103/PhysRevD.96.103528

Contaldi, C. R., Peloso, M., Kofman, L., and Linde, A. (2003). Suppressing the Lower Multipoles in the CMB Anisotropies. *J. Cosmol. Astropart. Phys.* 2003, 002. doi:10.1088/1475-7516/2003/07/002

de Blas, D. M., and Olmedo, J. (2016). Primordial Power Spectra for Scalar Perturbations in Loop Quantum Cosmology. *J. Cosmol. Astropart. Phys.* 2016, 029. doi:10.1088/1475-7516/2016/06/029

Elizaga Navascués, B., Martin de Blas, D., and Mena Marugán, G. A. (2018). Time-dependent Mass of Cosmological Perturbations in the Hybrid and Dressed Metric Approaches to Loop Quantum Cosmology. *Phys. Rev. D* 97, 043523.

Elizaga Navascués, B., Marugán, G. A. M., and Prado, S. (2020). Non-oscillating Power Spectra in Loop Quantum Cosmology. *Class. Quant. Grav.* 38, 035001. doi:10.1088/1361-6382/abc6bb

Elizaga Navascués, B., Martín de Blas, D., and Mena Marugán, G. (2018). The Vacuum State of Primordial Fluctuations in Hybrid Loop Quantum Cosmology. *Universe* 4, 98. doi:10.3390/universe4100098

Li, B.-F., Singh, P., and Wang, A. (2021). Phenomenological Implications of Modified Loop Cosmologies: an Overview. *Front. Astron. Space Sci.* 8, 701417. doi:10.3389/fspas.2021.701417

Martín-Benito, M., Neves, R. B., and Olmedo, J. (2021). *States of Low Energy in Bouncing Inflationary Scenarios in Loop Quantum Cosmology*.

Keywords: loop quantum cosmology, inflation, primordial perturbations, power spectra, quantum field theory in curved spacetimes, states of low energy, non-oscillatory vacuum state

Citation: Martín-Benito M, Neves RB and Olmedo J (2021) Non-Oscillatory Power Spectrum From States of Low Energy in Kinetically Dominated Early Universes. *Front. Astron. Space Sci.* 8:702543. doi: 10.3389/fspas.2021.702543

Received: 29 April 2021; Accepted: 28 July 2021;

Published: 17 August 2021.

Edited by:

Francesca Vidotto, Western University, CanadaReviewed by:

Jorge Pullin, Louisiana State University, United StatesSayantan Choudhury, National Institute of Science Education and Research (NISER), India

Copyright © 2021 Martín-Benito , Neves and Olmedo . This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mercedes Martín-Benito , m.martin.benito@ucm.es