Based on all electron relativistic density functional theory calculations, the properties of single oxygen vacancies in TiO 2 nanoparticles (NPs) have been obtained using a suitable representative model consisting of an octahedral (TiO 2 ) 84 nanoparticle of ∼3 nm size terminated with (101) facets. This nanoparticle can be safely considered at the onset of the so-called scalable regime where properties scale linearly with size toward bulklike limit, and hence results can be more directly compared to experiment. A set of reduced Ti 84 O 167 nanoparticles are selected to investigate the geometric, energetic, and electronic properties by using PBE semilocal functional with three different amounts of Fock exchange: 0% (PBE), 12.5% (PBEx), and 25% (PBE0). In particular, using the PBEx hybrid functional, previously validated for bulk anatase and rutile, it is predicted that the highly (three)-coordinated oxygen atom, located in the subsurface, and the least coordinated one at top sites are energetically the most suitable candidate for generating the oxygen vacancy. The subsurface case is in line with conclusions from experiments carried out on (101) single crystal anatase surfaces. The electronic structure of the reduced particles suggests that these would have better photocatalytic activity than their stoichiometric counterparts. Nevertheless, several properties of reduced TiO 2 NPs are strongly affected by the choice of the exchange-correlation functional, implying that, in absence of validation by comparison to experiment, predictions must be taken with caution.