An improved Hamiltonian constraint operator is introduced in loop quantum cosmology. Quantum dynamics of the spatially flat, isotropic model with a massless scalar field is then studied in detail using analytical and numerical methods. The scalar field continues to serve as 'emergent time', the big bang is again replaced by a quantum bounce, and quantum evolution remains deterministic across the deep Planck regime. However, while with the Hamiltonian constraint used so far in loop quantum cosmology the quantum bounce can occur even at low matter densities, with the new Hamiltonian constraint it occurs only at a Planck-scale density. Thus, the new quantum dynamics retains the attractive features of current evolutions in loop quantum cosmology but, at the same time, cures their main weakness.

Some long-standing issues concerning the quantum nature of the big bang are resolved in the context of homogeneous isotropic models with a scalar field. Specifically, the known results on the resolution of the big-bang singularity in loop quantum cosmology are significantly extended as follows: (i) the scalar field is shown to serve as an internal clock, thereby providing a detailed realization of the "emergent time" idea; (ii) the physical Hilbert space, Dirac observables, and semiclassical states are constructed rigorously; (iii) the Hamiltonian constraint is solved numerically to show that the big bang is replaced by a big bounce. Thanks to the nonperturbative, background independent methods, unlike in other approaches the quantum evolution is deterministic across the deep Planck regime.

Analytical and numerical methods are developed to analyze the quantum nature of the big bang in the setting of loop quantum cosmology. They enable one to explore the effects of quantum geometry both on the gravitational and matter sectors and significantly extend the known results on the resolution of the big bang singularity. Specifically, the following results are established for the homogeneous isotropic model with a massless scalar field: i) the scalar field is shown to serve as an internal clock, thereby providing a detailed realization of the 'emergent time' idea; ii) the physical Hilbert space, Dirac observables and semi-classical states are constructed rigorously; iii) the Hamiltonian constraint is solved numerically to show that the big bang is replaced by a big bounce. Thanks to the non-perturbative, background independent methods, unlike in other approaches the quantum evolution is deterministic across the deep Planck regime. Our constructions also provide a conceptual framework and technical tools which can be used in more general models. In this sense, they provide foundations for analyzing physical issues associated with the Planck regime of loop quantum cosmology as a whole.

The closed, k=1, FRW model coupled to a massless scalar field is investigated in the framework of loop quantum cosmology using analytical and numerical methods. As in the k = 0 case, the scalar field can be again used as emergent time to construct the physical Hilbert space and introduce Dirac observables. The resulting framework is then used to address a major challenge of quantum cosmology: resolving the bigbang singularity while retaining agreement with general relativity at large scales. It is shown that the framework fulfills this task. In particular, for states which are semi-classical at some late time, the big-bang is replaced by a quantum bounce and a recollapse occurs at the value of the scale factor predicted by classical general relativity. Thus, the 'difficulties' pointed out by Green and Unruh in the k=1 case do not arise in a more systematic treatment. As in k = 0 models, quantum dynamics is deterministic across the deep Planck regime. However, because it also retains the classical recollapse, in contrast to the k = 0 case one is now led to a cyclic model. Finally, we clarify some issues raised by Laguna's recent work addressed to computational physicists.

To every axi-symmetric isolated horizon we associate two sets of numbers, M n and J n with n = 0, 1, 2, . . ., representing its mass and angular momentum multipoles. They provide a diffeomorphism invariant characterization of the horizon geometry. Physically, they can be thought of as the 'source multipoles' of black holes in equilibrium. These structures have a variety of potential applications ranging from equations of motion of black holes and numerical relativity to quantum gravity.