The bulge represents the best compromise between old and massive Galactic component, and as such its study is a valuable opportunity to understand how the bulk of the Milky Way formed and evolved. In addition, being the only bulge in which we can individually resolve stars in all evolutionary sequences, the properties of its stellar content provide crucial insights on the formation of bulges at large. We aim at providing a detailed and comprehensive census of the Milky Way bulge stellar populations by producing deep and accurate photometric catalogs of the inner ~300 sqdeg of the Galaxy. We perform DAOPHOT/ALLFRAME PSF-fitting photometry of multi-epochs J and Ks images provided by the VVV survey to obtain deep photometric catalogs. Artificial star experiments have been conducted on all images to properly assess the completeness and the accuracy of the photometric measurements. We present a photometric database containing nearly 600 million stars across the bulge area surveyed by the VVV. Through the comparison of derived color-magnitude diagrams of selected fields representative of different levels of extinction and crowding, we show the quality, completeness and depth of the new catalogs. With the exception of the fields located along the plane, this new photometry samples stars down to ~1-2 mag below the MS-TO with unprecedented accuracy. To demonstrate the tremendous potential inherent to this new dataset, we give few examples of possible applications such as: i) star counts studies through the dataset completeness map; ii) surface brightness map; and iii) cross-correlation with Gaia DR2. The database presented here represents an invaluable collection for the whole community, and we encourage its exploitation. The photometric catalogs including completeness information are publicly available through the ESO Science Archive.

Recent observational programmes are providing a global view of the Milky Way bulge that serves as template for detailed comparison with models and extragalactic bulges. A number of surveys (i.e. VVV, GIBS, GES, ARGOS, BRAVA, APOGEE) are producing comprehensive and detailed extinction, metallicity, kinematics and stellar density maps of the Galactic bulge with unprecedented accuracy. However, the still missing key ingredient is the distribution of stellar ages across the bulge. To overcome this limitation, we aim to age-date the stellar population in several bulge fields with the ultimate goal of deriving an age map of the Bulge. This paper presents the methodology and the first results obtained for a field along the Bulge minor axis, at $b=-6^\circ$. We use a new PSF-fitting photometry of the VISTA Variables in the V\'{i}a L\'{a}ctea (VVV) survey data to construct deep color-magnitude diagrams of the bulge stellar population down to $\sim$ 2 mag below the Main Sequence turnoff. We find the bulk of the bulge stellar population in the observed field along the minor axis to be at least older than $\sim$ 7.5 Gyr. In particular, when the metallicity distribution function spectroscopically derived by GIBS is used, the best fit to the data is obtained with a combination of synthetic populations with ages in between $\sim$ 7.5 Gyr and 11 Gyr. However, the fraction of stars younger than $\sim$ 10 Gyr strongly depends upon the number of Blue Straggler Stars present in the bulge. Simulations show that the observed color-magnitude diagram of the bulge in the field along the minor axis is incompatible with the presence of a conspicuous population of intermediate-age/young (i.e. $\lesssim 5$ Gyr) stars.