Thermodynamic and dynamical properties of filling-control metal-insulator transition (MIT) in the Hubbard model are studied by the operator projection method, especially in two dimensions. This is a non-perturbative analytic approach to many-body systems. The present theory incorporates the Mott-Hubbard, Brinkman-Rice and Slater pictures of the MIT into a unified framework, together with reproducing low-energy narrow band arising from spin-charge fluctuations. At half filling, single-particle spectra A(ω, k) show formation of two Hubbard bands and their antiferromagnetic shadows separated by a Mott gap in the plane of energy ω and momentum k with lowering temperatures. These four bands produce splitting to two low-energy narrow bands and two SDW-like bands in the dispersion. Near half filling, the low-energy narrow band persists at low temperatures. This narrow band has a particularly weak dispersion and large weights around (π, 0) and (0, π) momenta. The velocity of these low-energy excitations is shown to vanish towards the MIT, indicating the mass divergence as in the Brinkman-Rice picture, but most prominently around (π, 0) and (0, π) with strong momentum dependence. This reflects the suppression of the coherence near the MIT. Main structures in A(ω, k) show remarkable similarities to quantum Monte-Carlo results in two dimensions as well as in the one-dimensional Hubbard model. The charge compressibility appears to diverge with decreasing doping concentration in both one and two dimensions in agreement with the exact and quantum Monte-Carlo results. We also discuss implications of the flat dispersion formed near the Fermi level to the observations in high-Tc cuprate superconductors.KEYWORDS: metal-insulator transition, Hubbard model, projection operator §1. IntroductionThe issue of the metal-insulator transitions (MIT) postulated by Mott 1) has opened extensive studies on the subject from both theoretical and experimental aspects.2) Pioneering works by Hubbard 3, 4, 5) and Kanamori 6) played an active part in theoretical studies on strong electron correlations in narrow-band electron systems.7) In spite of a long history of studies on strongly correlated electrons, 7, 2) an appropriate and satisfactory theoretical description of the MIT remains controversial, even for the simple Hubbard model. Particularly, in one and two dimensions, fluctuations become significant near the MIT. This makes it difficult to reliably describe physical properties related to the MIT. This paper aims at developing a formalism of operator projection 8,9,10) to construct a new non-perturbative analytic theory of many-body correlated electrons.11) This method allows us to systematically improve the Hartree-Fock theory, perturbation theories, Hubbard approximations 3, 5) and two-pole approximations. 12, 13) Using this new method, we give static and dynamical properties of the Hubbard model in one and two dimensions at and near half filling.Early theories on transitions from metals to insulators in clean systems are divided into three ...