This paper presents the results of a detailed theoretical investigation of the iron line formation NLTE problem in a three-dimensional model of the solar photosphere, which we have obtained from a very recent radiation hydrodynamics simulation of solar surface convection. In this Ðrst paper we have neglected the e †ects of horizontal radiative transfer on the atomic level populations, but we have considered a realistic atomic model for iron that contains hundreds of radiative transitions from the UV to the IR. The self-consistent solutions of the kinetic and transfer equations have been obtained with a new NLTE code, which is based on very efficient iterative methods. We Ðnd that overionization due to the near-UV radiation Ðeld does take place but mainly in the granular atmospheric regions. This wellknown NLTE mechanism tends to produce underpopulation of all the Fe I levels and a very small overexcitation of the Fe II levels. All over the three-dimensional photospheric model Fe II is the dominant ionization stage. We Ðnd signiÐcant LTE versus NLTE discrepancies mainly for the low-excitation Fe I lines. This applies to both the vertically emergent proÐles from the granular regions and also to the spatially averaged proÐles. These discrepancies are due to the line opacity deÐcits that result from the aforementioned underpopulation of the Fe I levels. The emergent proÐles of the low-excitation lines of Fe I are thus weaker in NLTE than in LTE. In particular, the largest errors in the equivalent widths (due to the LTE assumption) are found for the weakest low-excitation lines of Fe I. We also give quantitative estimates of the errors in the temperature structure of semiempirical solar granulation models obtained via the application of LTE inversion techniques to several groups of Fe I lines. For instance, the widely used Fe I 6301 and 6302 lines tend to lead to an overestimation of about 100È200 K in the Ó granular regions but to a similar underestimation in the intergranular plasma.The present paper considers also the case of the Sun observed with low spatial resolution, with particular emphasis on the long-standing iron abundance problem. We show that it is possible to obtain a very good Ðt to the observed spectral line shapes by slightly changing the iron abundance (for both the LTE and NLTE cases). In general, the iron abundance we need for reaching the best NLTE Ðt to observed equivalent widths is 0.074^0.03 dex larger than that needed to obtain the best LTE Ðt. Our most relevant conclusion with regard to the solar iron abundance issue is the following : if NLTE e †ects are fully taken into account in the three-dimensional model of the solar photosphere, we obtain the meteoritic iron abundance value However, if the abundance analysis is done assuming LTE, (A Fe \ 7.50). we Ðnd in close agreement with the recent LTE analysis of Asplund and collaborators. A Fe \ 7.43, Our results do indicate that NLTE e †ects are signiÐcant but not above the 0.1 dex level in the Sun. We consider our NLTE result for the i...