All electron relativistic density functional theory (DFT) based calculations using numerical atom-centered orbitals have been carried out to explore the relative stability, atomic, and electronic structure of a series of stoichiometric TiO2 anatase nanoparticles explicitly containing up to 1365 atoms as a function of size and morphology. The nanoparticles under scrutiny exhibit octahedral or truncated octahedral structures and span the 1-6 nm diameter size range. Initial structures were obtained using the Wulff construction, thus exhibiting the most stable (101) and (001) anatase surfaces. Final structures were obtained from geometry optimization with full relaxation of all structural parameters using both generalized gradient approximation (GGA) and hybrid density functionals. Results show that, for nanoparticles of a similar size, octahedral and truncated octahedral morphologies have comparable energetic stabilities. The electronic structure properties exhibit a clear trend converging:to the bulk values as the size of the nanoparticles increases but with a marked influence of the density functional employed. Our results suggest that electronic structure properties, and hence reactivity, for the largest anatase nanoparticles considered in this study will be similar to those exhibited by even larger mesoscale particles or by bulk systems. Finally, we present compelling evidence that anatase nanoparticles become effectively bulklike when reaching a size of similar to 20 nm diameter.