We consider quantization of the gravity-scalar field system in the minisuperspace approximation. It turns out that in the gauge fixed deparametrized theory where the scale factor plays the role of time, the Hamiltonian can be uniquely defined without any ordering ambiguity as the square root of a self-adjoint operator. Moreover, the Hamiltonian degenerates to zero and the Schrödinger equation becomes well behaved as the scale factor vanishes. Therefore, there is no technical or physical obstruction for the initial wave-function of the universe to be an arbitrary vector in the Hilbert space, which demonstrates the severeness of the initial condition problem in quantum cosmology.Although quantum mechanics has a lot of intriguing features, the evolution dictated by the Schrödinger equation is fully deterministic, which is no more different than classical wave propagation. The wave-function is uniquely determined in time once the initial state is somehow specified at the beginning. In principle, the initial state can be arbitrarily chosen by the external agent preparing the system inasmuch as the freedom of choosing initial conditions in classical mechanics.The situation seems to be different in an approach to quantum cosmology based on the Wheeler-DeWitt (WDW) equation due to its timeless characteristic. Some would argue that timelessness is a fundamental property of quantum gravity. One then hopes to impose certain (mathematically and physically motivated) boundary conditions for the WDW equation in quantum cosmology that are supposed to fix the wave-function of the universe uniquely. This, however, is not easily achievable and there are different viable suggestions like the no boundary [1] or the tunneling [2] proposals (see also [3,4] for other alternatives and [5] for an attempt to determine the initial state by referring to the special dynamics of the loop quantum cosmology).More recently, a Lorentzian path integral method based on the Picard-Lefschetz theory has been put forward as the basis of quantum cosmology and it has been claimed that both the no boundary and the tunneling proposals yield unsuppressed perturbations, thus they are problematic in the Lorentzian framework [6][7][8][9]. This has raised a debate in the literature, see [10][11][12][13][14][15][16][17]. The wave-function obtained by the Lorentzian path integral solves the WDW equation, at least in the minisuperspace model when the proper-time gauge is employed [18], and the recent discussion is again centered around the WDW approach.General relativity has many peculiarities that clearly distinguish it from a standard gauge theory of internal symmetries. To begin with the spacetime itself becomes dynamical rather than being a fixed background. Somehow relatedly, the "structure constants" of the gauge algebra of coordinate transformations turn out to be field dependent. Moreover, the Hamiltonian of general relativity vanishes and it becomes a constraint of the theory.Nevertheless, all these properties do not imply a perfect timelessness and indee...