On the geometry of no-boundary instantons in loop quantum cosmology

We study the geometry of Euclidean instantons in loop quantum cosmology (LQC) such as those relevant for the no-boundary proposal. Confining ourselves to the simplest case of a cosmological constant in minisuperspace cosmologies, we analyze solutions of the semiclassical (Euclidean) path integral in LQC. We find that the geometry of LQC instantons have the peculiar feature of an infinite tail which distinguishes them from Einstein gravity. Moreover, due to quantum-geometry corrections, the small-a behaviour of these instantons seem to naturally favor a closing-off of the geometry in a regular fashion as was originally proposed for the no-boundary wavefunction.


Introduction
In this work, we focus on the geometry of these Euclidean instanton solutions in LQC. This is a necessary first step before using such solutions to consider nucleation of universes from nothing and employing the measure provided by the associated wavefunction for predicting the probabilities of physically interesting phenomena for the Lorentzian histories. Our starting point shall be the (Euclidean) path integral for quantum gravity, along with the prescription that the initial conditions are provided by the no-boundary proposal. The main new ingredient, in comparison to the original Hartle-Hawking proposal, shall be the 'effective' action appearing in the path integral derived from LQC, as opposed to the usual Einstein-Hilbert one. (It was established in [4] that by replacing the standard FLRW action by the 'polymerized' version of it, the path integral formulation of LQC, in its phase space realization, retains all the crucial aspects of the quantum geometry which appear in the canonical LQC.) Other than this, our formalism shall be exactly the same as in the original no-boundary proposal: We shall look only at the saddle-point approximation of the Euclidean path integral and consider the wavefunction to be a functional of the value of the scale factor only at the final (spatial) boundary. Moreover, as shall be obvious throughout our paper, we make a minisuperspace approximation for all our calculations and consider the matter content to be only that due to a cosmological constant. The latter approximation ensures that we always have a compact instantonic solution and do not require dealing with subtleties which can give rise to Euclidean wormholes [5]. Since we want to show new features of the LQC instantons with respect to its geometry, as compared to the original Hartle-Hawking ones, these approximations shall help us emphasize our main result without unnecessarily complicating the system. For our purposes, the effective action consists of two main types of quantum corrections specific to LQC -the holonomy and inverse-triad modifications [6]. The first appears due to the fact that there are no quantum operators corresponding to the connection or extrinsic curvature (in other words, the momenta conjugate to the spatial metric) on the kinematic 1 Hilbert space of the theory. On the other hand, there are well-defined operators corresponding to the holonomy (or parallel transport) of the connection [7]. Therefore, one expresses the curvature operator in terms of these holonomies instead of the connection itself. Classically, one can take the limit such that one recovers the expression of the curvature written in terms of the connection from the expression given for the holonomies. However, the geometrical operators in the full loop quantun gravity have discrete spectra for quantities such as area and volume [8] on the kinematical Hilbert space spanned by the spin-network states, rendering taking such a limit unviable. Therefore, one inherits an 'area-gap', in analogy with the minimum energy-gap of the harmonic oscillator, from the full theory in LQC [9]. The main effect of this regularization of the curvature in terms of holonomies, for symmetry-reduced models, lie in replacing the extrinsic curvature by matrix elements of SU(2)-holonomies which are periodic functions of the connection. Specifically, for minisuperspace cosmologies, we have H → sin(δH)/δ, where H is the 1 'Kinematic' here refers to the fact that the Gauss and the (spatial) diffeomorphism constraints have been solved whereas 'physical' would imply the solution of the Hamiltonian constraint as well. For a minisuperspace model, these distinctions are not very important since the only leftover symmetry of the system is time-reparameterization invariance.
Hubble parameter and δ is related to the area-gap.
The discrete spectra of area and volume operators also lead to other type of corrections in LQC. The most significant of them are the inverse-triad corrections which arise from the requirement of having a well-defined operator corresponding to the inverse of some power of the scale factor whose spectra contains the zero eigenvalue. Naively, it is impossible to have a densely-defined operator in such a case. However, using the aforementioned holonomy operators and what is commonly known as the 'Thiemann trick' in the literature, one can express the relation [10] In this definition,ĥ := exp(iδp a ), for the momentum conjugate to the scale factor p a ∝ȧ, is precisely a SU(2)-valued holonomy operator mentioned previously. It is clear from this relation that one can have an operator, whose classical limit is some inverse power of the scale factor on the RHS although we do not require any inverse operator on the LHS [11]. Using this, one gets rid of the singular behaviour of any function which contains some inverse power of a due to the replacement by these aforementioned inverse-triad corrections. Once again, their form for minisuperspace cosmologies is rather simple, as shall be explicitly demonstrated later. Let us briefly summarize our main result. Conceptually, at least in the cosmological constant case, the main effect of the LQC quantum-geometry corrections lies in the small-a behaviour of the Euclidean LQC instantons. As shall be demonstrated, due to the inversetriad corrections, the LQC-modified Friedmann equation is such that the solution tails off to zero at the symmetry point of the theory. The geometry of the LQC instanton will emerge to be quite different from the original Hartle-Hawking proposal with an infinitely stretched tail in Euclidean time; however, their topology remains the same. Moreover, such an infinitely long tail of the instanton (in imaginary time) is not an inherent problem since the only meaningful physical quantity is the probability of nucleation which remains finite for this system. The interesting fact is the quantum-geometry regularization is such that this tail closes the geometry in a regular way without requiring any additional finetuning even though the field equations are heavily modified in LQC. This is suggestive of the fact that the no-boundary proposal is robust and, if anything, such a necessary tail-off of LQC instantons to zero points towards it being more natural in the presence of quantum-geometry corrections.

The Hartle-Hawking proposal revisited
In this section, we first briefly review the geometry of the no-boundary instantons in Einstein gravity. In the process, we also fix our notation for the rest of the paper.

The Wheeler-de Witt equation and boundary conditions
Let us consider Einstein gravity with a scalar field In this paper, we shall exclusively focus on a minisuperspace cosmological model [12] where σ 2 = 2/3π is some normalization constant 2 . By assuming the slow-roll limitφ ≈ 0, the Lagrangian can be simplified as whereV = 16V /9. From this, one can get the conjugate momentum p a = −aȧ/N. The Hamiltonian H is obtained by the usual Legendre transform where On quantization, by replacing p a = −i (d/da), one gets the Wheeler-de Witt equation for the wave function of the scale factor where γ is a constant due to the ambiguity in operator-ordering and In the semi-classical regime, the ambiguity due to operator ordering is not that important and can be ignored in a first approximation [12]. It is straightforward to see that the behavior of the system in the classically allowed region U < 0 (hence, a > 1/ √V ) and in the classically forbidden region U > 0 (hence, a < 1/ √V ) are different. For the classically allowed region, the solution is essentially oscillatory and can be a superposition of in-going and out-going modes; for the classically forbidden region, the solution is a superposition of exponentially growing and decaying modes.
In order to extract a specific solution from these general solutions, one needs to impose boundary conditions. However, quantum cosmology is a closed system in which a set of boundary conditions cannot be determined by the environment external to the setup as is the normal practice. In general, there is no fundamental principle to assign all the boundary conditions necessary to specify the wave function of the universe Ψ. At best, the general consensus is that these boundary conditions need to be supplied as additional fundamental laws of nature. There are two famous, mathematically consistent wavefunctions corresponding to the following boundary conditions: 1. The Hartle-Hawking proposal [2] -If we choose the exponentially growing mode for a < 1/ √V , then the wave function becomes a superposition of in-going and out-going modes for a > 1/ √V .

2.
The tunneling proposal [13] -If we choose the out-going mode for a > 1/ √V , then the wave function becomes a superposition of growing and decaying modes for a < 1/ √V . The probability of the universe nucleating from 'nothing' mainly depends on the contribution from the classically disallowed (hence, quantum) regime. Therefore, the (leading-order contribution to the) probability distribution is approximately where + corresponds to the Hartle-Hawking wavefunction, and − to the tunneling wavefunction, respectively.

Euclidean path integral
Since there is no fundamental principle to choose a boundary condition for the universe in general, we can look towards the path integral quantization for some guidance. In this paper, we shall take this route of quantizing gravity in the path integral formulation instead of the canonical one outlined above. We believe this is the most elegant presentation of the original Hartle-Hawking proposal [2] and shall, therefore, stick to it for the LQC-modified case [1]. However, as an aside, we note that the Wheeler-de Witt equation mentioned above transforms into a "difference" equation (of finite step-size) in LQC resulting from a modified Hamiltonian constraint on a non-separable Hilbert space arising in the theory. It is completely natural that one might impose boundary conditions analogous to the ones mentioned above for this difference equation in a canonical quantization [14]. Between two hypersurfaces (which take given values of the 3-metric and scalar field specified on it) at the initial boundary (h i , φ i ) and the final one (h f , φ f ), the (Feynmann) propagator is given by where we sum over all geometries, allowing for topology-changes, which have the specified initial and final boundaries. This integration is highly-oscillatory and ill-defined, but its convergence properties can perhaps be improved on introducing a Wick-rotation to Euclidean time dt = −idτ , where now we sum over all Euclidean geometries and corresponding field combinations with the given boundaries (left of Fig. 1). For usual quantum field theories, this corresponds to the ground state wavefunction. Although there is no straightforward way to define the ground state in quantum gravity, but Hartle and Hawking proposed that the above form of the Euclidean path integral may correspond to the ground state wavefunction of the universe.
Let us further assume that the intermediating Euclidean geometries between the initial and final boundaries are disconnected (right of Fig. 1). Especially, if we consider a closed universe (i.e., k = +1 FLRW cosmology), then the wavefunction remains well-defined if we remove the initial boundary This unique boundary condition has been given the so-called nickname no-boundary proposal, because it has no initial boundary.
There is a lot of justifiable controversy whether this Euclidean path integral is a good approximation of the original Lorentzian path integral or not [15,16,17,18]. However, this no-boundary proposal, defined as an Euclidean path integral, is attractive because of several nice properties it possesses.
• The Euclidean path integral can be interpreted as the partition function of a thermal system [19] whereĤ is the quantum Hamiltonian of the system and β is the inverse of the Hawking temperature. The right-hand side can be evaluated in the steepest-descent approximation as where F = E − T S is the Helmholtz free energy, E and S being energy and entropy of the system, respectively. Here, A is the area of the cosmological horizon. Since the ADM energy E is zero for dS space, we can consistently recover the Bekenstein-Hawking entropy formula S = A/4. This reveals that the Euclidean path integral is consistent with the expected thermodynamic properties of gravity.
• Building on this result, one may consider other semi-classical effects in dS space as well. The classical (Lorentzian) equation of motion for a scalar field on a fixed dS background is given byφ where H =ȧ/a. On inserting random 'white noise' in the slow-roll limit as a thermal effect due to Hawking temperature, one obtains the Langevin equation [20,21] where one imposes the white noise conditions: Then, the probability to have the field value φ at t will follow the Fokker-Planck equation [22] ∂P (φ, t) ∂t while the probability to have the field value initially χ at t = 0 will follow the equation In the static limit, a solution that satisfies both of these equations is given by We can interpret this as the tunneling probability of a homogeneous part of a universe that tunnels from the field value χ to φ via stochastic quantum fluctuations when the wavelength is of the order of the Hubble radius.
This wave function is consistent with the Euclidean path integral approximated by the Hawking-Moss instantons [23]. On further normalizing the initial boundary, one can obtain the no-boundary wave function. Therefore, we can conclude that the Euclidean path integral describes the stationary limit, or thermal equilibrium, of quantum fluctuations of the Hubble-scale wavelength modes consistently, whereas in many situations, such a thermal equilibrium is coincident with the ground state of the system [21].
The above physical motivations are reason enough to investigate and apply the Euclidean path integral with the no-boundary condition as the wavefuntion of the universe, not only due to its mathematical simplicity but also due to its self-consistency in the low-energy limit with quantum field theory in curved spacetime.

Semiclassical approximation: instanton solutions
One can calculate the Euclidean path integral in the saddle-point approximation by using Euclidean on-shell solutions, or so-called instantons, as The on-shell solution in pure dS space, on imposing the no-boundary condition a(0) = 0, becomes regular to give where H 2 0 = 8πV /3. This solution reveals that at the a(0) = 0 ("South Pole") point, we need to impose the conditionȧ = 1 from the Hamiltonian constraint.
Inserting this solution, one can evaluate the Euclidean action. In the phase space formulation, the action takes the form On plugging the on-shell condition, H = 0 is automatically satisfied. Hence, This solution can be analytically continued to Lorentzian time for any constant-τ hypersurface, but after the Wick-rotation, the metric is in general complex-valued except at thė a = 0 hypersurface (i.e., τ = π/2H 0 , left of Fig. 2). If we Wick-rotate on this surface, then the metric becomes (right of Fig. 2) We give more mathematical details corresponding to this solution in the following section.

Geometry of the Hartle-Hawking instantons
We give a more detailed derivation of the schematics described for the minisuperspace model in the previous section. The Friedmann equation (or the Hamiltonian constraint) for the k = 1 FLRW universe, in Euclidean time, is given by where we use the metric ds 2 = N 2 (η)dη 2 + a 2 (η) dΩ 2 3 and set G = 1 throughout 3 . η denotes the Euclidean time parameter while reserving t for Lorentzian time, as before. We rewrite the other relevant equation which is the scalar field equation (also, in Euclidean time) 3 Note the new normalization chosen here to facilitate comparison with LQC later on.
First of all, let us make a gauge choice and fix the lapse function N = 1. As pointed out in [24], this can be rigorously achieved by introducing the complex variable τ (η) = η 0 dη ′ N(η ′ ). Given any lapse function, the variable τ defines a complex contour on the τplane. Once we rewrite the above equations in terms of the variable τ , the task of finding the no-boundary instantons is to solve these equations for the pair of complex analytic functions a(τ ) and φ(τ ), given the appropriate boundary conditions.
The set of equations, in terms of this new variable, can be expressed aṡ where a dot refers to a derivative with respect to τ and the Hubble parameter H :=ȧ/a. The fact that we can write the RHS of (29) as a function of the scale factor alone is only possible for the simplest case of a massless scalar or a cosmological constant. There is always the Raychaudhuri equation involving the second derivative of a but only two of these three equations are linearly-independent. For our purposes of examining the onshell Euclidean instantons, required for estimating the path integral by its saddle-points, considering these two equations is sufficient. The usual procedure is to solve the above equations for the 'no-boundary' boundary conditions [25]: a(0) = 0 andφ(0) = 0. The first condition is a requirement that the geometry must close in a regular fashion while the second is a necessary condition for keeping the solution for φ regular as a → 0 [5]. It is often customary to quote another conditionȧ(0) = 1; however, this is the consequence of the Friedmann equation. In general, requiring that the scale factor and the scalar field take some fixed value on the final surface, (a(τ f ) = b, φ(τ f ) = χ), and some fixed value initially, (a(0) = 0,φ(0) = 0), exhausts all the conditions necessary to give a unique solution. The value of the derivative of the scalar field must be fixed from the scalar field equation (30) to be zero while the value of the scalar field at the 'South Pole' -φ(0) = φ 0 -gives the oneparameter family of instantonic solutions which satisfy the no-boundary proposal [24,25]. Additional tunings are necessary to ensure the classicality of our universe at late-times, the details of which are unimportant for our purposes (see [24,26]). If we take the pure gravity model, in the absence of any scalar field, one can analytically solve for the Euclidean instanton to find that there is a O(5)−symmetric solution given by a(τ ) = 3/Λ sin Λ/3 τ (this is Eqn. (23) above written in the new normalization). In the presence of a scalar field, one requires that the potential is sufficiently flat, i.e. of the inflationary type, for the solution to be regular. In the case of a slowly varying potential, the solution for a(τ ) is a deformed version of the sine function. However, for the massless scalar field (i.e. in the absence of any potential term at all), there are no compact instantons which would give rise to a nontrivial universe. This is obvious from the fact that in this case, the scalar field is in the "no roll" condition and the energy density of the universe is trivial. However, this scenario, which is beyond the scope of this paper, leads to new solutions for the no-boundary proposal in the presence of LQC corrections [1,27].
Before going on to the LQC instantons, let us revisit the geometry of these Hartle-Hawking instantons in Einstein gravity. (This is schematically shown in the right panel of Fig. 1.) Restricting to the case of pure gravity is already sufficient to illustrate its salient features. For Lorentzian signatures, one has the usual dS solution of Einstein's equations with a positive cosmological constant as ds 2 L = −dt 2 + a 2 L (t)dΩ 2 , with a L (t) = 3/Λ cosh Λ/3 t (equivalent to (26) above in the new normalization). To get the O(5) invariant Euclidean instanton from this result, one analytically continues t such that τ = τ f + it, where τ f is the point where a(τ ) reaches its maximum value. In other words, one gets a Lorentzian dS spacetime from a Euclidean hemisphere (right of Fig. 2), by matching the two hypersurfaces across the zero-extrinsic curvature (ȧ = 0) 'bounce' surface. The effective potential, V(a), goes to zero on this surface. Such a sharp transition between a real Euclidean half-sphere and a real Lorentzian part is only possible for the simplest example of a pure cosmological constant considered in this paper. In general, in the presence of an inflationary-type potential, the transition would be in terms of 'fuzzy' Euclidean instantons [24,25], whereby the solutions would be complex in some parts [28]. These details, however, are not important for us while focusing near the 'South Pole' to exhibit the general 'shuttle-cock' type shape of the no-boundary instanton in Euclidean gravity.
For the pure dS model, the Friedmann equation, in Euclidean time, takes the forṁ a 2 = 1 − Λa 2 /3 which clearly shows that as a → 0, one getsȧ → 1. The geometric interpretation of this result goes as follows. The Euclidean 4-metric, possessing the O(5) symmetry, ds 2 = dτ 2 + a 2 dΩ 2 3 has to smoothly close-off in a regular manner into flat space (written in spherical coordinates) ds 2 = dr 2 + r 2 dΩ 2 3 . For this to happen, one has to identify a(τ ) ∼ τ as a → 0. This suggests thatȧ → 1 in this limit, as is required from the Hamiltonian constraint. However, as mentioned before, this requirement for the derivative of the scale factor automatically follows from the constraint and is not part of the no-boundary condition.
Let us make one last comment before presenting our new results for the LQC instantons. The no-boundary initial condition is simply that the geometry closes off smoothly, as denoted by a(0) = 0. This shall be important later on for the LQC instantons. If we try to solve the Friedmann equation, we get the (famous) unique solution (23) only on imposing the no-boundary condition. Of course, choosing the 'initial' point at τ = 0 is only for convenience. A priori, there is no need for such a condition to be satisfied by the modified field equations in LQC. However, as we shall show in the LQC case, the initial condition that compact (Euclidean) instantons in LQC go off to a → 0 is favored naturally even in the presence of quantum-geometry corrections, at least in the pure gravity case. This is the main result of our work which we shall elaborate on in the following sections.

No-boundary instantons in LQC
In [1], it was shown how the effective LQC action modifies the instantons in the theory even for a simple cosmological constant. Moreover, this leads to slight enhancement of the probability of nucleation of the dS universe from nothing due to the LQC corrections. The general shape of the instanton is reproduced in Fig 2. However, we shall only be interested in the small-a behaviour of this LQC instanton. In particular, what emerges to be intriguing is the infinite tail of the instanton. Although this infinite tail is in Euclidean time, and therefore not physically relevant directly, it does have certain distinguishing features which we shall demonstrate below.
We begin with the modified Friedmann in LQC [29] due to the quantum geometry corrections mentioned in the Introduction.
whereρ is the contribution from a positive cosmological constant.
where v(a) := K 3 with V 0 = 16π 2 . Although we have adhered to the conventions of [29], we have generalized their results by adding in non-perturbative expressions for the inverse-triad corrections. On first look, the above set of equations look rather complicated due to the different terms involved. Here, f (a) represents the inverse-triad corrections whereas holonomy modifications show up in ρ 1 and ρ 2 . Importantly, note that due to the presence of the inverse-triad corrections, it is always possible to impose the no-boundary condition a → 0. However, to gain some intuition into the modified equation, let us begin by setting the holonomy corrections to zero for simplicity. This would be like taking the area gap ∆ to zero. In this limit, ∆ → 0, we geṫ It is easy to check that in the large a ≫ 1 limit, one gets f (a) ≈ v(a), and therefore we get back the usual Friedmann equation for a closed universe. However, in the a ≪ 1 limit, one cannot make such an approximation. Instead, in this limit, we get f (a) ≈ v(a) 2 . Our aim in this work is not to solve for those instantons which extremizes the Euclidean path integral but rather to examine its small-a behaviour. Therefore, considering v(a) ∝ a 3 and reinstating the holonomy modifications, we get the leading order term for a ≈ 0 aṡ for some constant C > 0. To obtain this result, we notice that the leading order term comes from theρ 2 term in (31) whereas the remaining terms are subdominant. This is a term which arises only in the quantum-corrected Friedmann equation (there is no term quadratic in the energy density in the classical Friedmann equation). Theρ 2 term comes with an additional minus sign which leads to a C > 0. Moreover, note that the dominant contribution in the classical case comes from the curvature term (1/a 2 ) whereas that term, contained in ρ 1 , is now sub-dominant. As already argued in the previous section, the essential requirement of the no-boundary condition is the geometry should be closed off in a regular manner and the condition onȧ should follow from the Hamiltonian constraint. In the LQC case, the modified Hamiltonian constraint impliesȧ = 0 instead of 1. Nevertheless, V(a) remains regular even in this case. The above findings for no-boundary instantons in LQC in quite remarkable. In order to appreciate this properly, let us make a few comments. Firstly, note that there was no reason that the modified Friedmann equation have to allow for the a → 0 limit to be imposed consistently. It could easily have been that this limit is singular in LQC. To illustrate this, let us consider only the holonomy modifications while ignoring the inversetriad ones. Typically, for the Lorentzian effective trajectories such an approximation is completely justified and valid even near the 'bounce' surface. In this case, the RHS of (31) has a singular term coming from 'ρ 1 ρ 2 ' (proportional to ∆ 2 ), which is absent in the classical  Figure 4: The black curve is an example ofȧ 2 for Λ = 1, G = 1, and l P l = 0.1, where the red dashed curve is the limit of the Einstein gravity with the same a max that satisfieṡ a max = 0. Right is the behavior near the a = 0 limit.
case. However, luckily for us, when one considers Euclidean histories as is required for our case, one cannot ignore the inverse-triad corrections any longer. Secondly, the structure of the modified Friedmann equation is such that the resulting instantons remain regular for all values of τ even on imposing the no-boundary condition. For a counterexample, imagine if the form of the equation was such thatȧ 2 ∝ a n with n > 2, in that case, the limit τ → 0 would have been singular and there would not have been consistent no-boundary instantons in the theory. Moreover, these inverse-a modifications not only play a crucial role in ensuring that the no-boundary condition can be imposed but also modify the geometry of these instantons to distinguish them from the Einstein gravity. It could also have been the case that the inverse-a modifications are such that one still gets the same condition forȧ at the South Pole. In that case, although the explicit solutions of the instantons would have been different, there would have been no difference in the geometry of LQC and Einstein gravity instantons. The quantum-geometry corrections in LQC conspire to ensure that we have no-boundary instantons in the theory with such a geometry that is tapers off to the symmetry point in a novel fashion.
Concretely, the small-a solution for a is given by a 0 e cτ . Obviously, this implies that the point τ = 0 is not a good point to impose the no-boundary condition. Rather both (a(τ ),ȧ(τ )) goes to zero as τ → −∞. However, this does not represent any difficulty since this infinite stretching is in Euclidean time and is, thus, not physically relevant directly. This novel feature of LQC instantons can be seen from Fig. 3, where the tail of the compact instanton is stretched infinitely, asymptotically tapering off to zero. The tail does contribute to the probability of nucleation of the Lorentzian dS universe, although the path integral and consequently the probability remains finite and well-defined in spite of this infinitely stretched geometry.

Numerical results
We give some sample numerical solutions for the LQC instantons, going beyond the smalla limit, to illustrate our claims regarding their infinite tail. Figs. 4 and 5 show a typical shape of the effective potential V LQC (a) and its solution a(τ ) for the HR, respectively. This solution demonstrates that the instanton is indeed infinitely stretched (Fig. 3). The approximate behavior of the nucleation probability is of the form (Fig. 6) with a model dependent positive constant c, where A = 4πa 2 max . As mentioned before, a max denotes the max value the instanton takes in the Euclidean regime, on which surface we analytically continue to the Lorentzian regime. The above constant c can easily be absorbed away in the normalization of the probability measure and the only relevant correction due to LQC, over the Einstein-Hilbert value, comes from the parameter d. This parameter d can be approximately expanded by d ≃ 8.7 × l 2 P l /γ (Fig. 7) in terms of the fundamental parameters -Planck length and Immirzi parameter -of the theory, as shown via numerical reconstruction.

Robustness of the no-boundary condition
In Einstein gravity, when considering a positive cosmological constant as the only matter source, the only solution which is regular at the South Pole is given by (23). The crucial point here is that one gets this solution for the instanton in Einstein gravity from the (Euclidean) Friedmann equation on imposing the no-boundary condition. Interestingly, as was shown in previous section, the small-a behaviour for the LQC instantons following the modified Friedmann equation is always of the form  Figure 7: Left: By varying l P l , one can calculate d ∝ l 2 P l numerically. Right: By varying γ (with l P l = 0.1), one can see a linear dependence d ∝ 1/γ. We can numerically conclude that d ≃ 8.7 × l 2 P l /γ.
for some constant C. This shows that as τ → −∞, we get a → 0 implying a natural implementation of the no-boundary condition in LQC instantons, as a result of quantumgeometry regularizations, at least for the simplest case of a cosmological constant. In this sense, this shows that the no-boundary condition is robust and more natural when the Einstein-Hilbert action is augmented by LQG corrections. The Euclidean instantons, in this pure gravity scenario, are always going to be compact and there is no risk of an Euclidean wormhole forming for some values of the parameter space. This conclusion is true both for Einstein gravity as well as the LQC case. Therefore, one always gets compact Euclidean instantons for the pure dS case in both cases. The overall conceptual picture of the nucleation of the universe from nothing is also the same in both cases. One quantitative difference is due to the modified equations in LQC: The 'bounce' surface (ȧ = 0), as predicted by the effective trajectories of sharply-peaked semiclassical states, in LQC is (slightly) different from the hypersurface joining the Euclidean and Lorentzian parts in Einstein gravity. However, the qualitative behaviour remains the same and this difference is reflected in enhancing the probability of nucleation of the universe manifesting as a next-to-leading-order correction in LQC. As mentioned in [1], these type of terms can also appear in Einstein gravity from going to higher order corrections and there should be a competition between the two terms in LQC -one appearing from radiative corrections and the other from the inherent quantum geometry. It is thus difficult to unambiguously state that LQC enhances the tunneling amplitude for such no-boundary universes.

Conclusion
It is an old expectation that mathematical consistency alone shall be sufficient to derive the boundary conditions in quantum cosmology [30]. The form of the LQC instantons suggest that it might indeed be possible to identify typical smooth initial conditions due to quantum-geometry corrections. This demonstrates that one of the most fundamental proposal for the initial condition in quantum cosmology can appear naturally in another -LQC -due to putative corrections coming from quantum geometry. The introduction of the no-boundary proposal in LQC has also opened new physical possibilities for the latter. Instead of replacing the big bang singularity with a deterministic bounce, as is predicted by semiclassical states in some restricted models of LQC, this opens up the opportunity to allow for Euclidean trajectories leading to a bubble nucleation of our universe. The semiclassical saddle-point approximation of the no-boundary proposal is distinct from the semiclassical sharply-peaked states in LQC and therefore the former acts as an example of how a new state can unveil novel features in a well-established theory. At this point, it is difficult to compare the probability of a bounce versus that of tunneling of the universe from nothing. However, the no-boundary proposal also provides a (set of) natural initial conditions for considering inhomogeneous perturbations in LQC leading to effects observable from early-universe cosmology without having to resort to ad hoc choices for the initial state.
Regarding the geometry of no-boundary instantons in LQC, we have demonstrated that the feature of having an infinite tail distinguishes these instantons from the original Hartle-Hawking ones in Einstein gravity. We end our discussion with a few caveats. Firstly, it has been pointed out recently that the Euclidean path integral in gravity is not a good approximation for the original Lorentzian path integral due to several conceptual reasons [15]. However, even if one works with the Lorentzian path integral and applies a different mathematical trick (Pecard-Lefshetz theory) to improve its convergence, the resulting theory typically has runaway perturbations due to the old conformal factor problem in gravity [16]. Interestingly, LQC can come to the rescue of the no-boundary proposal [31], written as a Lorentzian path integral, even in this case. However, the main physical effect from LQC responsible for this is 'dynamical signature-change' [32], something we have ignored in this work as a first pass. The second caveat is regarding the fact that our discussions were limited to the case of a pure cosmological constant in this paper. Indeed, the more interesting physical scenario is that of having a scalar field in some potential. However, preliminary investigations have already revealed that the solution space of the no-boundary wavefunction is greatly enhanced for such a system in the presence of LQG corrections. Finally, one can ask how physical is the fact that the tail of these instantons are stretched to infinity? As already mentioned, this is only true in Euclidean time and therefore not directly meaningful. However, it might even be possible that for some different gauge choice (i.e. N = 1), one can even avoid such an infinite stretching altogether. Nevertheless, all the interesting effects of having such a geometry as explained in this paper would still be valid in this case. Most importantly, no matter what the gauge choice, the remarkable conclusion that the modified Friedmann equation in LQC not only allows for the no-boundary condition to be imposed but also somehow makes it more natural seems to be robust and points towards a new paradigm in quantum cosmology merging these two mainstream approaches.