Minisuperspace quantum cosmology in f(Q) gravity

f(Q) symmetric-teleparallel gravity is considered in view of quantum cosmology. Specifically, we derive cosmological equations for f(Q) models and then investigate the related energy conditions. In the Minisuperspace formalism, the point-like f(Q) Hamiltonian is taken into account. In this framework, we obtain and solve the Wheeler–De Witt equation, thus finding the wave function of the universe in different cases. We show that the Hartle criterion can be applied and classical observable universes occur.


I. INTRODUCTION
Recent observations on supermassive black holes [1,2] and gravitational waves [3] further confirm the validity of General Relativity (GR) as a theory describing gravitational phenomena also in the strong field regime.These results are complementary to the well-known Solar System tests in the weak field limit.Despite these undoubted successes, Einstein's theory cannot be considered as the final theory for the gravitational interaction because it presents shortcomings at UV and IR scales [4][5][6].Furthermore, the impossibility to deal with GR under the same standard as the other interactions is one of the main weakness of the theory.Indeed, at quantum level, UV divergences cannot be canceled out through standard renormalization techniques [7,8].This problem has been firstly addressed by considering the asymptotic safety scenario [9][10][11], according to which there should be nontrivial Gaussian fixed points in the renormalization group flow driving the values of the coupling constants in the UV regime.Another solution was proposed in Ref. [12], where the author extends the action including second-order curvature invariants and makes the resulting theory renormalizable, albeit the price to pay is the loss of unitarity.
Although fixing UV problems in GR is one of the most active research areas nowadays, there is no final answer to the open questions occurring in GR, meaning that no self-consistent theory of Quantum Gravity so far exists.At the astrophysical and cosmological levels, incompatibilities with observations yielded the introduction of dark energy [13] and dark matter [14,15], which, according to the current picture of the universe, should account for the majority of energy-matter content.The main puzzle is that, at the moment, any experimental attempt to address observed large-scale dark components as new fundamental particles has given no final result.
The aforementioned problems are only few examples of the unsolved puzzles exhibited by GR [16].However, such difficulties led to the introduction of several different theories of gravity, modifying e.g. the main assumptions of GR such as the Lorentz Invariance [17][18][19] or the Equivalence Principle [20][21][22], or extending the gravitational action to functions of higher-order curvature invariants [23][24][25][26][27][28].To the latter category belongs the so called f (R) gravity, which extends the Hilbert-Einstein action by including non-linear functions of the scalar curvature.This theory provides several interesting results at any scales.For instance, at galactic scales, a power-law model seems to be capable of fitting the galaxy rotation curves without any dark matter [29][30][31][32].At cosmological scales, also dark energy can be mimicked without introducing any cosmological constant [33][34][35][36][37][38].Moreover, some models are potentially capable of addressing the mass-radius relation of Neutron Stars without introducing exotic equations of state (EoS) [39,40], or the early stage of the inflationary universe without additional scalar fields [6,41].At the quantum level, several models have been constructed with the purpose of addressing UV shortcomings of GR [42][43][44][45][46][47].
Within the large amount of alternatives to Einstein's theory, particular interest has been recently gained by those theories involving torsion or non-metricity.Specifically, it can be shown that gravity can be described using torsion or non-metricity instead of curvature, and the resulting theories are dynamically equivalent to GR [48][49][50][51][52][53].As pointed out in Sec.II, this is due to the fact that theories including either torsion scalar T or non-metricity scalar Q, instead of curvature scalar R, yield exactly the same field equations as GR.As torsion is strictly related to the anti-symmetry of the affine connection with respect to the lowest indexes, non-metricity arise when considering a non-vanishing covariant derivative of metric tensor, namely ∇ µ g αβ = 0 where isometry is not conserved.The advantage of these approaches are, among the others, that Lorentz Invariance and Equivalence Principle are not required a priori like in GR.Furthermore, they can be dealt with within the framework of gauge theories [22].These features could allow, if any, the possibility of a quantum approach to gravity avoiding shortcomings of GR.Furthermore, they seems to solve naturally several problems related to the cosmological dark side [54][55][56][57].
In this paper we focus on the cosmological extension of symmetric teleparallel theory whose action contains a function of the non-metricty scalar Q, i.e. f (Q).More precisely, we consider the function f (Q) = βQ + αQ s , with α, β and s being free parameters.In this way, GR is straightforwardly recovered as soon as α → 0.
In the first part of the paper, we consider a Friedman-Lemaître-Robertson-Walker (FLRW) space-time with non-vanishing spatial curvature and find cosmological solutions of the field equations.Then we study the corresponding energy conditions and finally, recasting the model in the Hamiltonian formalism, we derive the Wheeler-De Witt (WDW) equation to obtain the so called Wave Function of the Universe.Specifically, the latter approach is based on the Arnowitt-Deser-Misner (ADM) formalism [58][59][60][61], where the Hamiltonian of GR is derived by means of a (3 + 1) decomposition of the metric.The final result is a Schrödinger-like equation, the WDW equation, whose solution provides the Wave Function of the Universe.The meaning of the latter is still debated because its probabilistic meaning, as in the standard Quantum Mechanics, cannot be applied.However, several interpretations have been provided [61][62][63][64] and, despite the lack of a direct probabilistic interpretation, it gives important information on the initial conditions and evolution of the universe at early stages.Particularly relevant is the Hartle criterion, according to which an oscillating wave function can be related to classical dynamical systems interpretable as observable universes.
The paper is organized as follows: Sec.II is a summary of Symmetric Teleparallel Equivalent of General Relativity (STEGR) and its modifications.In Sec.III, we consider the cosmology of the f (Q) = βQ+αQ s model, derive and solve the corresponding field equations and study the related energy conditions.The ADM formalism, introduced in Appendix A, is then applied to the modified STEGR model in Sec.IV.The related WDW equation is solved in different cases.In particular, we show that the Hartle selection criterion (see Appendix B) holds and classical observable universes can be recovered.Conclusions are reported in Sec.V.

II. THE SYMMETRIC TELEPARALLEL EQUIVALENTS OF GENERAL RELATIVITY
As mentioned in the introduction, several possible extensions/modifications of GR can be proposed; a large part of them relaxes some assumptions of GR.This is the case of theories violating the Lorentz Invariance, or assuming general affine connections which can be non-metric compatible.See [53] for a discussion.In the latter case, one of the possible modifications consists in relaxing the assumption of symmetric connection with respect to the low indexes.As a consequence, the Equivalence Principle is not required at the foundation of the theory and torsion occurs in the description of the space-time dynamics [50].In this way, the affine connection Γ α µν can be written in terms of the standard Levi-Civita connection Γα µν plus an additional contribution related to the torsion; clearly, GR is recovered as soon as the connection is assumed to be of the Levi-Civita one.In a space-time with non-vanishing curvature and torsion, the general connection can be written as where K α µν is the contorsion tensor defined as If one chooses the Teleparallel Gauge, where curvature vanishes identically, the gravitational action turns out to be equivalent to the Hilbert-Einstein one and the resulting theory is called Teleparallel Equivalent of General Relativity (TEGR).In particular, setting Γ α µν = 0, Eq. ( 1) becomes Thus, by defining the torsion superpotential and the torsion scalar, respecively as it turns out that the action is dynamically equivalent to the Hilbert-Einstein one and, therefore, it leads to the same field equations.It is worth noticing that, using tetrad fields, TEGR can be recast as a gauge theory of translation group in the flat tangent space-time [65] Similar considerations apply relaxing another assumption of GR, that is the metricity principle.Imposing the covariant derivative of metric to be non-vanishing, i.e. ∇ α g µν ≡ Q αµν , the affine connection turns out to contain an additional contribution, called disformation tensor and denoted by L α µν .In a general space-time with curvature and non-metricity, the connection reads as From the above definitions, it is possible to introduce the so called non-metricity scalar Q as or, in a more compact form, as with Now, to obtain a theory dynamically equivalent to GR, one can choose the so called coincident gauge, where the total affine connection vanishes, i.e.Γ α µν = 0. Therefore, the disformation tensor turns out to be equivalent to the Levi-Civita connection, up to a sign: As a consequence, the STEGR action, differing from the Hilbert-Einstein one by a boundary term, is This is due to the fact that the scalar curvature (written in terms of the Levi-Civita connection) and the non-metricity scalar are equivalent up to a total divergence.Specifically, considering Eqs. ( 9)-( 11), one gets with ∇α being the covariant derivative written in terms of the Levi-Civita connection.Also notice that, in the coincident gauge, the non-metricity scalar can be written in terms of the Levi-Civita connection and the metric as [48]: According to the above considerations, the gravitational interaction can be equivalently described either by curvature, or by torsion, or by non-metricity.In summary, the general connection containing all possible contributions reads Einstein's gravity is recovered for K α µν + L α µν = 0, so that connection is Levi-Civita.In the teleparallel gauge, the Levi-Civita contribution turns out to be equivalent to the contorsion tensor, so that the space-time can be described by torsion only.In the coincident gauge, where K α µν = Γ α µν = 0, non-metricity is the only non-vanishing component.The actions corresponding to these three theories are totally equivalent up to a boundary term.This means that the related field equations are exactly the same.

A. The f (Q) extension
As pointed out in [66], non-metricity, torsion and curvature can be thought as a geometric Trinity of Gravity, since all of them leads to totally equivalent theories by different formalism.Nonetheless, the equivalence among the three models also implies that TEGR and STEGR suffer the same shortcomings as GR at cosmological and astrophysical levels.For this reason, in analogy to GR extensions like f (R) gravity, f (T ) and f (Q) extensions are considered in the literature [53,[67][68][69][70][71][72][73][74].However, it is worth pointing out that actions involving f (R), f (T ), and f (Q) are not equivalent, though they can be made equivalent by introducing the corresponding boundary terms [75][76][77].As a matter of fact, f (R) gravity is different with respect to the f (T ) and f (Q) theories, since the corresponding boundary terms B T and B Q play a non-trivial role in the dynamics.
For instance, f (T ) and f (Q) gravity lead to second-order field equations [22], while f (R) gravity gives fourth-order field equations.Moreover, under proper assumptions, f (T ) and f (Q) models can avoid bad ghosts and lead to unitary theories [78].Interestingly, equivalences among the three theories can occur when considering the so called Noether Symmetry Approach [79][80][81][82].Specifically, when the actions are selected by Noether symmetries, it turns out that the related symmetry generators can be made equivalent under proper definitions of the free parameters [83].
The features of f (T ) and f (Q) gravity are today investigated in the literature through the applications to cosmology and astrophysics.
For instance, in [55], authors show that a f (T ) power-law model can fit the galaxy rotation curve.In [84], the H 0 tension is addressed.In [85], an EoS, derived from f (T ) can address the dark energy issue.In [86] a model-independent way to solve the modified Friedman equations, in the framework of f (T ) cosmology, is proposed.
Regarding f (Q) gravity, in Ref. [87] the Big Bang Nucleosynthesis is studied, while, in [88], slow-roll inflation is considered.In [71], bouncing cosmological models are considered within the context of f (Q) gravity.In [89], authors investigate the polarization of gravitational waves.In [90], f (Q) gravity is showed to be potentially capable of challenging the ΛCDM model and, in [91], a model-independent reconstruction of f (Q) cosmology is reported.In [92], wormhole solutions are provided by static and spherically symmetric backgrounds.To conclude a wide class of phenomena and models can be investigated under the standard of f (Q) formalism.
Considering dynamics generated only by non-metricity, the variation of f (Q) action yields the field equations [93] with T µ ν being the matter stress-energy tensor.Here we want to derive exact solutions of f (Q) field equations, in a cosmological background with non-vanishing spatial curvature.After, we will develop the Quantum Cosmology formalism to find out the Wave Function of the Universe.Such a wave function, according to the Hartle criterion, allows to recover the classical evolution of observable universes.
Let us consider the modified symmetric teleparallel action in a cosmological background.We are going to adopt the Lagrange Multipliers method to find a point-like Lagrangian and the related equations of motion.The Lagrange multipliers method is a general mathematical technique used to deal with constraints into optimization or variational problems.In the context of modified gravity, the Lagrange multipliers emerge as constraints highlighting possible further degrees of freedom of theories different with respect to GR.Specifically, by introducing them, one can construct modified action functionals.Variations of these actions with respect to the Lagrange multipliers and the canonical variables yield the modified equations of motion.Lagrangians can be determined through this procedure, as they are derived from the modified action functionals [94].It is worth to point out that the specific form of the Lagrangian, obtained by the Lagrange multipliers, strictly depends on the particular theory of gravity.In the specific case of this paper, we are going to consider the cosmological expression of the non-metricity scalar, which infers the Lagrange multiplier λ, appearing as a constraint in the point-like Lagrangian.To this purpose, let us start from the action containing a function of the non-metricity scalar Q.We can get a point-like Lagrangian considering a FLRW metric of the form with k being the spatial curvature and a the scale factor.Using the Lagrange multipliers method and the cosmological expression of Q with respect to the space-time (19), the action, up to a total factor, can be recast as The Lagrange multiplier λ, as standard, can be found by varying the action with respect to the non-metricity scalar, providing then the cosmological point-like Lagrangian turns out to be From the above Lagrangian it is possible to obtain three partial differential equations, namely the two Euler-Lagrange equations with respect to a and Q: The system must be implemented by the so called energy condition, i.e.
It is important to stress that dynamics gives the evolution of the scale factor.This is because the equation with respect to the non-metricity scalar is a constraint imposed by the Lagrange multiplier.Furthermore, the zero energy condition, directly inferred from the Einstein equations, is another constraint completing the dynamical system.In particular, the vanishing energy condition implies that the related Hamiltonian is zero on the constraint surface due to time reparameterization invariance [95].As a result, the variational approach, naturally providing two different Friedman equations by the variation with respect to the metric, turns out to be completely equivalent to the Lagrangian approach.
On the other hand, the equation with respect to the scale factor corresponds to the (1,1) Einstein equation, while the Euler-Lagrange equation with respect to Q provides the cosmological expression of the non-metricity scalar.Specifically, being we have and then which is the definition of the non-metricity scalar in FLRW space-times.The situation is completely equivalent to that in f (T ) teleparallel equivalent gravity, where the point-like Lagrangian results reduced and independent of the time derivative Ṫ , with T being the torsion scalar [96].This reduction mechanism, related to the Lagrange multipliers, is present in several modified gravity models, see e.g.[79].Albeit, in principle, the non-metricity scalar can be written as a function of a and ȧ, the Lagrangian formalism implies that Q and a have to be treated as separated fields and their relation is thus recovered by the Euler-Lagrange equation with respect to Q, namely Eq. ( 24).This approach is standard in the literature (see e.g.[97][98][99][100]) and will be used to reduce dynamics.
Eqs. ( 23) and ( 25) read, respectively, as For k = 0, the above equations reduce to Let us now consider the function f (Q) = βQ + αQ s , with s, α, β, real constants.It is the straightforward generalization of STEGR.The Lagrangian becomes Starting from Eq. ( 32), we can find cosmological solutions with generic spatial curvature k and, therefore, study the energy conditions in terms of the free parameters.

A. Cosmological Solutions
For α and β different from zero, the general solution of the Euler-Lagrange equations can be found analytically only for k = 0. We have with Considering s = 2, we have a Starobinsky equivalent model in non-metric gravity with a scale factor of the form The related inflationary model has been studied in [88].In order to obtain an accelerating behavior, we must set β/α > 0.
Notice, moreover, that the solution ( 33) is obtained by imposing s = 1, 1/2.As a matter of fact, for the latter value, field equations are identically satisfied independently of the scale factor.The ΛCDM model is clearly recovered for s = 0 and β = −1/2, where Eq. ( 33) can be recast as 1.The Case α = 0 Analytic solutions with non-vanishing spatial curvature are present only if either α = 0 or β = 0.In the former case, we have with c 1 being an integration constant.As in GR, when the spatial curvature vanishes, the cosmological scale factor becomes trivial in vacuum.

The Case β = 0
Let us now solve system (29) by setting β = 0.The only analytic scale factor solving the field equations is If k = 0, only trivial scale factors occur as solutions of the cosmological equations.Moreover, in the special case s = 1/2, the Euler-Lagrange equations turn out to be trivially satisfied for any scale factor.Also notice that this solution only depends on the spatial curvature k, unlike f (R) = αR s gravity, where the cosmological dynamics is determined by the value of the parameter s.In both cases, GR is recovered for s = 1 and little deviations can be studied by setting s = 1 + , with 1.In this case, the function f (Q) = αQ 1+ can be expanded around = 0, providing and the additional logarithmic contribution can behave like an effective cosmological constant.

B. Energy Conditions
The energy condition for this model can be studied for different values of the free parameters s and α.In order to recover the gravitational coupling, we set β = −1/2, so that the RHS of the field equations ( 29) can be though as effective energy density ρ and pressure p provided by geometry, whose expressions are with H = ȧ/a being the Hubble parameter.Therefore, the EoS parameter w reads Notice that by plugging the solution (33) into Eq.( 40), the energy density and the pressure turn out to be constant and opposite.Specifically, they are Hence, the EoS parameter is constantly equal to −1, independently of the time in which it is evaluated.This means that the additional contribution αQ s behaves like an effective cosmological constant.Clearly, in this case, the Null Energy Condition and the Dominant Energy Condition vanish identically, while the Weak Energy Condition and the Strong Energy condition can be violated depending on the values assumed by s and α.For instance, choosing a coupling constant ranging from 0 up to 20 time the GR coupling constant, the WEC is always violated, meaning that f (Q) gravity can act as an exotic fluid.In other words, the further gravitational degrees of freedom resulting from f (Q) affect dynamics giving a further fluid to source the field equations.See also [101,102] for a discussion on generalized energy conditions Quantum Cosmology is an important tool to investigate features of the early universe.In particular, it is an approach needed to fix the initial conditions from where dynamical systems, describing observable universes, come out after some prescriptions are satisfied [62][63][64].In view of reproducing the standard of Quantum Mechanics results, the Hamiltonian formalism and a quantization procedure have to be pursued.
The Quantum Cosmology approach can be applied to different theories of gravity [103][104][105] and results are particularly useful in the framework of the so called Minisuperspace formalism [106].Minisuperspaces are suitable reductions of the 3-metrics configuration space (the Superspace) where cosmological dynamics can be constructed.
The application of this procedure allows to obtain the WDW equation from the related Hamiltonian and then the so called Wave Function of the Universe, a quantity connected to the probability to achieve configurations of cosmological variables which can be interpreted (or not) as initial conditions of the universe (or of the "universes" in a Many Worlds Interpretation of Quantum Mechanics).
Despite the meaning of the Wave Function of the Universe remains uncertain, numerous explanations have been proposed throughout the years.For instance, according to the Many World Interpretation, the wave function originates from quantum measurements that occur simultaneously in different universes, without any collapse of the wave function [61].Another perspective was presented by Hawking, who suggested that the wave function is connected to the probability that early universe evolves into our classical universe [62,63].In [64,107], the emergence of the universe from "nothing" is explored, viewing the universe from a quantum standpoint.This viewpoint helps in solving the challenge of determining the boundary condition, where it is asserted that the fundamental laws of physics must be intrinsic and not imposed externally [108].In this interpretative framework, J. B. Hartle proposed a criterion for extracting information from the wave function based on its behavior in the late-time epoch.Specifically, according to the Hartle criterion, the wave function should exhibit oscillatory behavior within the classically allowed region.Oscillations of the wave function mean that cosmological parameters are correlated.As stated by Hartle, this feature corresponds to describe classical universes [109].This interpretation of the wave function is analogue to that provided by non-relativistic Quantum Mechanics, where the solution of the Schrödinger equation for a particle inside a potential barrier, i.e. in the classically forbidden region, yields a decaying exponential behavior while the wave function outside the barrier gives rise to an oscillating behavior.Adopting the Hartle criterion, in the semiclassical limit, allows us to express the wave function in terms of the classical action S as ψ ∼ e iS and then deriving classical cosmological solutions (See the Appendix B for details).
In the specific case of f (Q) cosmology, after a Legendre transformation of Lagrangian (32), the related Hamiltonian is where is the conjugate momentum related to the scale factor a(t).Notice that there is no momentum related to the non-metricity scalar, due to the reason discussed in Sec.III.In fact, though the Minisuperspace S = {a, Q} is two-dimensional, the resulting tangent space T S = {a, ȧ, Q} is three-dimensional.The Hamiltonian in Eq. ( 43) encompasses the two variables a and Q, treated as separate quantities.Their correlation is reestablished through the Hamilton-Jacobi equations, equivalent to the Euler-Lagrange equations by construction.In this scenario, adopting the Lagrangian formalism becomes crucial, as the Hamiltonian can be simply derived through a straightforward Legendre transformation and the problem of infinite-dimensional Superspace can be overcome by a related Minisuperspace (see Appendix A for details).
In the general case, that is for α, β = 0, the WDW equation can be exactly solved only by setting k = 0.In such a case, the Hamiltonian turns out to be and by solving the differential equation Hψ = 0, one gets where we have defined Here c 1 and c 2 are dimensionless complex coefficients and J ± 1 6 a 3 p(Q) are Bessell functions of first kind.Setting, for instance, s = 2, we have f (Q) = βQ + αQ 2 , and then ) This is the symmetric teleparallel analogue of the Starobinsky model whose inflationary counterpart has been studied in Ref. [88].It is worth noticing that, for large values of a, the wave function keeps an oscillating behaviour, while in the early stages its trend appears exponential-like.The plot of the Wave Function of the Universe in terms of a is reported in Fig. 1.This is in agreement with the Hartle Criterion [110] by which it is possible to recover the classical behavior of the universe.Specifically, the first step is to adopt the semiclassical limit [111] and recast the wave function as ψ ∼ e iS0 , with S 0 being the Hamilton-Jacobi action.Since we are interested in finding classical solutions, we can take the limit of Eq. ( 46) in the late times, where the combination of the Bessel functions is oscillating.Specifically, up to total factors, the classical action can be identified as the argument of the Bessel functions, namely where c 0 is a constant parameter.Formally, from the above action, we get the conjugate momenta: where, from the discussion in Sec.III, it results π Q = 0.
Comparing the above momenta with those provided by Lagrangian (32), we obtain the following equations: where Eq. ( 52) is an algebraic constraint resulting from the degenerate momentum π Q = 0.The cosmological solution is in agreement with those obtained in Sec.III where the scale factor evolves exponentially and the non-metricity scalar is a constant.Moreover, from the relation occurring between the scale factor and the non-metricity scalar, it is possible to constrain the parameter c 0 .
It is worth noticing that, using the Hartle Criterion, we can recover classical solutions by asking for the existence of Noether symmetries, see e.g.[24,103,104,112].In this case, such a criterion demonstrates to be a selection rule.See Appendix B.
Let us now focus on the cases where either α = 0 or β = 0. Setting α = 0 in Lagrangian (32), the Hamiltonian becomes and the related WDW equation can be solved even without imposing k = 0.It provides with i being the imaginary unit and D are parabolic cylinder functions.The wave function becomes trivial in the case of a spatially-flat universe with k = 0. From Eq. ( 55), by applying the Hartle Criterion and pursuing the same approach as in the previous subsection, we get classical solutions.In particular, as showed in Fig. 2, we notice that the parabolic cylinder functions asymptotically behave like the exponential function; therefore, taking the real part and considering that the wave function consists of a linear combinations of functions with subscript −1/2, the Hamilton-Jacobi action, up to a total factor, can be written as: with a conjugate momentum given by By equating Eq. ( 57) with we get the differential equations whose solution is in agreement with results reported in Eq. (37).Once again, it is worth stressing that the Hartle Criterion can be used to find classical solutions without directly solving the field equations.Another interesting case, which deserves to be studied separately, is the case β = 0. Here, the only analytic solution occurs for k = 0, where the Hamiltonian becomes and, therefore, the wave function turns out to be This case is interesting because it is possible to suitably change the variables to introduce a cyclic coordinate into the Minisuperspace derived from the Noether symmetries [79].Specifically, we can operate the coordinate transformation {a, Q} −→ {z, w} , which explicitly is the Lagrangian becomes which is cyclic with respect to z.After a Legendre transformation, we obtain where π z is the conjugate momentum of the new cyclic variable z(t).Now, since z is cyclic, we have meaning that π z is a constant: Therefore, the Wave Function of the Universe is a solution of the following differential equations: where the first is the WDW equation and the second is related to the momentum conservation.The general solution is which is an oscillating wave function satisfying the Hartle Criterion as follows.Specifically, adopting the semiclassical limit [111] and recasting ψ in terms of the action S 0 as the Hamilton-Jacobi equations can provide back the equations of motion for the new variables.More precisely, the system of Hamilton-Jacobi equations It can be easily shown that the system (72) is equivalent to the system of Euler-Lagrange equations with respect to a and Q, when coming back to the old variables.As a matter of facts, by considering the inverse transformations of Eq. ( 64), that is the second equation of (74) becomes from which namely the Euler-Lagrange equation with respect to Q. Equivalently, after few manipulations, the first equation of ( 74) takes the form: which is the equation of motion with respect to the scale factor.Notice finally that, also here, the Wave Function of the Universe is oscillating and therefore the Hartle Criterion is satisfied.

V. CONCLUSIONS
In this work we considered the modified STEGR model f (Q) = βQ + αQ s and investigated the related Minisuperspace Quantum Cosmology after discussing the FLRW cosmology and the related energy conditions.The WEC is identically violated and the further degrees of freedom of f (Q) gravity can be dealt as a geometric perfect fluid behaving as dark energy.
The starting point for Quantum Cosmology is the Hamiltonian formulation.From this, it is possible to achieve the WDW equation and then the Wave Function of the Universe.It can provide information on the early stages of the universe after the application of the Hartle Criterion.Interestingly, we found that the wave function trend is oscillating only for large scale factors, that is in the classical zone related to observable universes, while it keeps an exponential trend at early times.This confirms the validity of the Hartle Criterion for this model.In this regard, we showed that a suitable change of variables allows to recast the Wave Function of the Universe as ψ → e iS0 , with S 0 being a classical action.Therefore, by means of the Hamilton-Jacobi equations, the cosmological solutions can be straightforwardly recovered even without solving the field equations directly.Clearly, due to the absence of a self-consistent theory of Quantum Gravity, this approach represents a first step and a sort of effective model, which can be useful for cosmological purposes at early times.Moreover, modifications of STEGR can be extremely helpful to fix shortcomings suffered by GR, since inflationary behaviors and dark energy effects arise without introducing any cosmological constant, but provided by additional geometric contributions in the field equations.This prescription can be also applied to the f (R) extension of GR but, unlike f (R) gravity, modified STEGR leads to second-order field equations, which can be handled more easily.From this point of view, f (Q) gravity can be useful to fix several problems of GR at different scales of energy.
If the system shows symmetries, then it is possible to introduce cyclic variables and conserved momenta.Therefore, if π i q i are the momenta related to the Minisuperspace variables q i and j i the conserved quantities, we have with q i being cyclic.Once the conjugate momenta are promoted to quantum operator, the above relation becomes If the system carries m symmetries (and corresponding conserved quantities j 1 , j 2 . . .j m ), Eq. (B2) and the WDW equation provide a system of m + 1 differential equations: Since j i are conserved quantities, the wave function will be oscillating in the Minisuperspace of variables q i ; in particular, the latter m equations of the above system provide the following solution [105,112]: where m is the number of symmetries, are the directions where symmetries do not exist, n is the total dimension of Minisuperspace.Notice that also the component |χ > can generally depend on the conserved quantities j i , though its functional behavior depends on the specific case considered.In summary, the Hartle Criterion can be related to symmetries, which, in turn, can select classical solutions which can be interpreted as observable universes.

Figure 1 :
Figure 1: Wave Function of the Universe as a function of the scale factor