Novel phenomena of the Hartle-Hawking wave function

We find a novel phenomenon in the solution to the Wheeler-DeWitt equation by solving numerically the equation assuming \(O(4)\)-symmetry and imposing the Hartle-Hawking wave function as a boundary condition. In the slow-roll limit, as expected, the numerical solution gives the most dominant steepest-descent that describes the probability distribution for the initial condition of a universe. The probability is consistent with the Euclidean computations, and the overall shape of the wave function is compatible with analytical approximations, although there exist novel differences in the detailed probability computation. Our approach gives an alternative point of view of the no-boundary wave function from the wave function point of view. Possible interpretations and conceptual issues of this wave function are discussed.