The Rayleigh-Benard theory by Grossmann and Lohse [J. Fluid Mech. 407, 27 (2000)] is extended towards very large Prandtl numbers P r. The Nusselt number N u is found here to be independent of P r. However, for fixed Rayleigh numbers Ra > 10 10 a maximum around P r ≈ 2 in the N u(P r)-dependence is predicted which is absent for lower Ra. We moreover offer the full functional dependences of N u(Ra, P r) and Re(Ra, P r) within this extended theory, rather than only giving the limiting power laws as done in ref. [1]. This enables us to more realistically describe the transitions between the various scaling regimes, including their widths.In thermal convection, the control parameters are the Rayleigh number Ra and the Prandtl number P r. The system responds with the Nusselt number Nu (the dimensionless heat flux) and the Reynolds number Re (the dimensionless large scale velocity). The key question is to understand the dependences Nu(Ra, P r) and Re(Ra, P r). In experiments, traditionally the Prandtl number was more or less kept fixed [2][3][4]. However, the recent experiments in the vicinity of the critical point of helium gas [5,6] and of SF 6 [7] or with various alcohols [8] allow to vary both Ra and P r and thus to explore a larger domain of the Ra−P r parameter space of Rayleigh-Benard (RB) convection, in particular that for P r ≫ 1. While the experiments of Steinberg's group [7] suggest a decreasing Nusselt number with increasing P r, namely Nu = 0.22Ra 0.3±0.03 P r −0.2±0.04 in 10 9 ≤ Ra ≤ 10 14 and 1 ≤ P r ≤ 93, the experiments of the Ahlers group suggest a saturation of Nu with increasing P r for fixed Ra, at least up to Ra = 10 10 [9]. The same saturation (at fixed Ra = 6 · 10 5 ) is found in the numerical simulations by Verzicco and Camussi [10] and Herring and Kerr [11].The large P r regime of the latest experiments has not been covered by the recent theory on thermal convection by Grossmann and Lohse (GL, [1]), which otherwise does pretty well in accounting for various measurements. In particular, it explains the low P r measurements of Cioni et al.[4] (P r = 0.025), the low P r numerics which reveal Nu ∼ P r 0.14 for fixed Ra [10, 11], and the above mentioned experiments by Niemela et al. [6] and Xu et al. [8].In the present paper we extend the GL theory in a natural way to the regime of very large P r, on which no statement has been made in the original paper [1]. We find Nu to be independent of P r in that regime. We in addition present the complete functional dependences Nu(Ra, P r) and Re(Ra, P r) within the GL theory, rather than only giving the limiting power laws and superpositions of those as was done in [1]. This enables us to more realistically describe the transitions between the various scaling regimes found already in [1].Approach: To make this paper selfcontained we very briefly recapitulate the key idea of the GL theory, which is to decompose in the volume averages of the energy dissipation rate ǫ u and the thermal dissipation rate ǫ θ into their boundary layer (BL) and bulk contributions...