Deep PSF photometric catalog of the VVV survey data

The Vista Variables in the Via Lactea (VVV) survey has performed a multi-epoch near-infrared imaging of the inner Galactic plane. High-fidelity photometric catalogs are needed to utilize the data. We aim at producing a deep, point-spread-function (PSF) photometric catalog for the VVV survey J, H, and Ks band data. Specifically, we aim at taking advantage of all the epochs of the survey to reach high limiting magnitudes. We develop an automatic PSF-fitting pipeline based on the DaoPHOT algorithm and perform photometry on the stacked VVV images in J, H, and Ks bands. We present a PSF photometric catalog in the Vega system that contains about 926 million sources in the J, H, and Ks filters. About 10% of the sources are flagged as possible spurious detections. The 5 sigma limiting magnitudes of the sources with high reliability are about 20.8, 19.5, and 18.7 mag in the J, H, and Ks band, respectively, depending on the local crowding condition. Our photometric catalog reaches on average about one magnitude deeper than the previously released PSF DoPHOT photometric catalog. It also includes less spurious detections. There are significant differences in the brightnesses of faint sources between our catalog and the previously released one. The likely origin of these differences is in the different photometric algorithms that are utilized; it is not straightforward to assess which catalog is more accurate in which situations. Our new catalog is beneficial especially for science goals that require high limiting magnitudes; our catalog reaches such in fields that have a relatively uniform source number density. Overall, the limiting magnitudes and completeness are different in the fields with different crowding conditions.

The VVV images contain hundreds of millions of stars above the detection limits of the survey.Obtaining photometric measurements for all these sources poses a fundamental challenge for data processing and analysis because performing accurate photometry in the crowded conditions of the Galactic plane is very difficult.Various approaches to confronting this challenge exist; the choice between them is not trivial and depends on the qualities desired for the data (see, e.g., discussion in Sects.3.1 and 3.2.3).This basic problem is crucially enhanced by the sheer amount of data in the VVV survey; handling them requires fully automated, parallelized procedures for all relevant tasks.Indeed, various photometric techniques have been applied to the VVV survey data in the past (Saito et al. 2012a;Mauro et al. 2013;Alonso-García et al. 2018;Surot et al. 2019a,b).Recently, the VVV science team has published the first point spread function (PSF) photometric catalog for the whole VVV survey (Alonso-García et al. 2018).The catalog contains approximately 846 million sources with detections in at least three filter bands of Z, Y, J, H, and K s ; detections in the J, H, and K s bands have been made for about 534 million of these sources (see Sect. 2 for further details).
In this paper, we develop a new PSF photometric pipeline and apply it to the VVV survey J-, H-, and K s -band data.
Our work is motivated by the interest to reach deeper limiting magnitudes than those of the catalog by Alonso-García et al. (2018).This interest is driven by our goal to perform dust extinction mapping in the Galactic plane in the future (e.g., Kainulainen et al. 2011Kainulainen et al. , 2013;;Kainulainen & Tan 2013;Butler et al. 2014;Mattern et al. 2018; for that purpose, detecting as many sources as possible is desired.We show that this can indeed be achieved by exploiting different data reduction and photometry schemes: our final photometric catalog reaches on average roughly one magnitude deeper than the catalog of Alonso-García et al. (2018).Here, we describe the automated pipeline that performs the photometry, present its results, and release the catalog to the community.

VVV survey data
Here, we use the near-infrared imaging data of the VVV survey.The VVV survey covers the Galactic bulge (−10 ≤ l ≤ 10, −10 ≤ b ≤ 5) and part of the adjacent Galactic plane (−65 ≤ l ≤ −10, −2 ≤ b ≤ 2).Below, we describe the data to the degree required to understand the implementation of our data reduction and photometry procedures (described in Sect.3).For the detailed description of the survey and the data, we refer to the comprehensive survey papers cited in the text below.
The survey instrument, VIRCAM, is equipped with 16 detectors 2048 × 2048 pixels2 in size, with a pixel scale of 0.339 .The detectors are arranged in a 4 × 4 array with gaps along the X and Y axes.A single pointing with the detector array is called a pawprint and consists of 16 single detector images (SDIs), covering ∼0.595 deg 2 on the sky.As with all VISTA observations, the basic observational unit is the observation block (OB1 ).There are two types of OB in the VVV survey: multi-filter, single-epoch OBs (JHK s and ZY OBs), and K s -only variabilitymonitoring OBs.The VVV OBs use the "Tile6n+Jitter2u" observing pattern, which means there are six pawprint positions and two small offsets at each pawprint position.Figure 1 shows the "Tile6n" dithering pattern for one VISTA detector.Two exposures at each pawprint position can be combined to construct a stacked pawprint.Six stacked pawprints can be combined to construct a tile that covers a contiguous field of ∼1.687 deg 2 .This value is slightly larger than the average size of a VVV tile (∼1.64 deg 2 ) because the adjacent VVV tiles share some overlaps.More information about the instruments and observing strategy can be found in Sutherland et al. (2015), Minniti et al. (2010), andSaito et al. (2012a).Altogether, the coverage of the VVV survey consists of 348 tiles, including 196 tiles in the bulge (tile names start with "b"), and 152 tiles in the disk (tile names start with "d").All tile centers and numbers can be found in Saito et al. (2012a).
Table 1 shows the exposure time of VVV OBs.Each pawprint is a co-addition of n exposures and therefore the exposure time per pawprint (Exp p ) is equal to the DIT (detector integration time) multiplied by the NDIT (number of detector integration time).In the bulge area, the Exp p is 2 × 6 s, 4 s, and 4 s for J, H, and K s , respectively.In the disk area, the Exp p is 2 × 10 s for J, H, and K s for the multi-filter, single-epoch OBs, and 4 s for the K s -only variability-monitoring OBs.The stacked pawprints are constructed by combining two small-offset jitter pawprints and thus the exposure time per stacked pawprint (Exp sp ) is 2 × Exp p .In one epoch, six stacked pawprints can be combined to obtain a tile image.However, considering the large The VVV data have been processed using VISTA data flow system (VDFS) pipeline at the Cambridge Astronomical Survey Unit (CASU 2 ; Irwin et al. 2004;Cross et al. 2009;Lewis et al. 2010;Saito et al. 2012a).The process of VVV data reduction and one approach to calibration and photometry can be found in Saito et al. (2012a) and González-Fernández et al. (2018).Here we only summarize the main data-reduction process.First, the dark current is subtracted from each raw image and then the linearity is corrected for every detector.The flat-field correction is also applied by dividing by the twilight flats.Second, the sky background model is constructed with the exposures in a given filter within several concatenated OBs.Considering the variability of the near-infrared sky, at least 24 input pawprints closest in time (within ∼30 min) are combined to produce the background sky model that is then subtracted from each input pawprint.Finally the jitter stacking is conducted to construct the stacked pawprints after correcting the stripe pattern introduced by the detector readout electronics.
We use the JHK s images and source catalogs from Data Release 4 (DR4) provided by the VISTA Science Archive (VSA3 ; Hambly et al. 2008;Cross et al. 2012).We note that VSA DR4 includes one photometric catalog based on aperture photometry.The 5σ limiting magnitude of the catalog Notes. (a) Detector integration time. (b) Number of DITs. (c) Exposure time per pawprint. (d) Exposure time per stacked pawprint.
is K s ∼ 17−18 mag for most tiles (Saito et al. 2012a), and it drops to K s ∼ 15−16 mag in the crowded fields close to the Galactic center (Saito et al. 2012a).Recently, Alonso-García et al. ( 2018) performed PSF photometry on the multi-epoch VVV images using DoPHOT (Mateo & Schechter 1989;Schechter et al. 1993) and used the average fluxes over all epochs to derive the brightnesses of the detected sources.The 5σ limiting magnitude of the catalog is K s ∼ 18 mag in the VVV disk, K s ∼ 17.5 mag in the VVV innermost bulge, and K s ∼ 18 mag in the VVV outermost bulge.This DoPHOT PSF photometric catalog is also publicly available through the VSA DR4.

Methods: PSF photometry
In this section, we first justify our choices for the photometric techniques (Sect.3.1) and then describe the details of our PSF photometry procedure (Sect.3.2).The resulting catalog is presented in Sect. 4. To avoid ambiguity, we introduce the nomenclature used.The multi-epoch stacked pawprints are combined to construct the deep images.We refer to them as deep pawprints.Meanwhile, we refer to a SDI from a stacked pawprint as a stacked SDI.We combine the stacked SDIs of the same detector from all multi-epoch stacked pawprints together to obtain the deep SDIs (see Sect. 3.2.1).

DaoPHOT vs. DoPHOT
Aperture photometry is known to be relatively inaccurate for crowded fields; this is definitely so for the sensitive VVV data, especially at the latitudes close to the Galactic plane and in the Galactic bulge.To achieve accurate photometry in the areas of high stellar density, and to reach deeper limiting magnitude than with the aperture photometry, we perform PSF photometry of the VVV data.
We considered two commonly used PSF-fitting programs, namely DoPHOT and DaoPHOT (Stetson 1987).DoPHOT uses an analytic function to model the stellar PSF.It is unable to fit more than two stars simultaneously although DoPHOT considers double stars.The main advantage of DoPHOT is that it runs very quickly and in an entirely automated manner.DaoPHOT combines both the analytic and empirical methods to model the stellar PSF.It fits the central part of a bright star with an analytic function and then uses look-up tables to do corrections from the analytic function to match the observed stellar profile.Once a PSF model is derived, DaoPHOT tries to fit many stars simultaneously.To do that, DaoPHOT first separates all stars into groups and then performs PSF-fitting for all stars in a group simultaneously.The disadvantage of DaoPHOT is that it requires careful configuration of input parameters and is relatively slow.
Fig. 2. Flow chart of our automatic PSF-fitting pipeline that is used to do the PSF photometry on the VVV survey imaging data.Mauro et al. (2013) compared the DaoPHOT photometric results with DoPHOT and aperture photometry using the VVV survey data of the globular cluster M 22.These latter authors found that DaoPHOT can return a higher detection rate of faint sources and therefore goes deeper than DoPHOT and aperture photometry.Because we wish to include as many sources as possible in the extinction mapping, we decided to use DaoPHOT to construct a new, deep PSF photometric catalog of the VVV data.

Automatic PSF-fitting pipeline
We developed an automatic pipeline to perform PSF photometry on the VVV survey imaging data.The pipeline uses the DaoPHOT algorithm and is mainly written in IDL and Python.The pipeline is adapted to run in multi-core mode, which significantly reduces the required CPU time.Figure 2 shows the flow chart of the pipeline.We give the detailed description of the pipeline in the following sections (Sects.3.2.1-3.2.4).

Image stacking
To reach as deep a photometry as possible, we stack together multi-epoch images for each filter, that is, J, H, and K s .In principle, there are two options for stacking, one starting with tile images and the other starting with stacked pawprints.Considering that each detector of the detector array potentially has a different PSF and that seeing conditions can vary between images, the PSF is expected to vary across a tile.This makes it difficult to model the PSF with a universal function.To avoid the problems arising from the PSF variation, we decided to work with the stacked pawprints (as also suggested by Mauro et al. 2013).
In practice, we download "good" (deprecated=0)4 stacked pawprints of multiple epochs for each tile region from VSA.For the n epoch observations in one tile region for one filter, there are usually 6 × n stacked pawprints.We split each stacked pawprint into 16 stacked SDIs.In each filter band, about 6×n stacked SDIs corresponding to the same detector are checked and then combined using the SWARP software (Bertin et al. 2002) after excluding the stacked SDIs obtained in bad weather conditions.The detailed process is as follows.First, we calculate a common center based on the centers of the stacked SDIs that correspond to the same detector.We then reproject and resample each stacked SDI onto a common reference frame defined by the common center; this is done using the bilinear interpolation based on the astrometric WCS information of the stacked pawprints provided by CASU.We also subtract a constant back-ground from each stacked SDI during the above resampling.Finally, the original resampled stacked SDIs are co-added to the deep SDIs with SWARP.We note that the median WCS rms of the CASU astrometric solution is ∼70 mas (Saito et al. 2012b).The position uncertainty introduced by the SWARP resampling process is ∼0.03 pixels, corresponding to ∼10 mas.Thus, the position accuracy of the deep SDIs is dominated by the CASU astrometric solution.We finally obtain 16 deep SDIs in each filter band of each tile region.As an example, Fig. 3 (top panels) shows the deep SDIs for one detector of the tile "b214" in J, H, and K s bands.To avoid the non-negligible field distortion across a tile image (Sutherland et al. 2015;Irwin et al. 2004), we did not mosaic these 16 deep SDIs to a stacked tile.
We note that the CASU data-reduction pipeline used fast bilinear interpolation to obtain the stacked pawprints.We also use bilinear interpolation to stack the multi-epoch stacked SDIs.However, the bilinear interpolation can introduce spatially correlated noise in the final images (as pointed out by Meingast et al. 2016).Such a phenomenon is illustrated in Fig. 3 (bottom panels), which shows the background rms maps of deep SDIs for one detector of tile "b214".We can see obvious Moiré patterns in J and H background rms maps.Unfortunately, the spatially correlated noise introduced by the bilinear interpolation can also affect photometry.It can produce Moiré patterns in the color space, which would further transfer into systematic patterns in extinction maps (the extinction mapping is based on colors of stars).Our solution is to lessen the effect of these patterns during the photometric calibration process (explained in Sect.3.2.3).While this produces satisfactory results for our goal, complete A85, page 4 of 19 elimination of the spatially correlated noise would require a re-processing of all VVV survey raw data using more suitable interpolation methods (e.g., a higher-order resampling kernel as suggested by Meingast et al. 2016); this is clearly beyond the scope of our paper.

Source detection and photometry
The source detection and photometry is performed on the deep SDIs.To better handle the large amount of data and parallelise the photometry, we split each deep SDI into pieces of ∼1000 × 1000 pixels.The PSF photometry is performed on each piece with PyRAF 5 that is a Python front for running IRAF 6 tasks.The source detection is performed using the DAOFIND task with the signal-to-noise ratio threshold of S /N > 3. The PSF function is modeled with the PSF task of DaoPHOT.Finally, we use ALLSTAR task to obtain the photometric results.
To illustrate the quality of the process, Fig. 4 shows two partial regions of the tiles "d003" and "b305" before and after PSF subtraction in K s band.The tile "d003" is located in the VVV disk region and represents a relatively uncrowded field.The tile "b305" represents a highly crowded field close to the Galactic center.The example shows that most stars have been well fitted and subtracted from the stacked images.The residuals around bright stars in the PSF-subtracted images are mainly due to the saturation and nonlinearity effects.
For each piece of the deep SDIs, an instrumental PSF-fitting photometric catalog is obtained from the output of the DaoPHOT ALLSTAR task.We apply the astrometric WCS information recorded in each corresponding deep SDI (see Sect. 3.2.1).We create a single-band catalog by cross-matching the WCS coordinates of the sources in all pieces of the 16 detectors of a VVV tile.This is done using the STILTS package (Taylor 2006) with a tolerance of 0 .5. We then merge the single-band catalogs using STILTS with a tolerance of 1 (Saito et al. 2012a).For bandmerging about 93-98% (with a mean value of ∼96%) of the sources are matched within 0 .5.

Photometric calibration
The absolute photometric calibration was obtained by a comparison of the instrumental magnitudes of relatively isolated, unsaturated bright sources with their counterparts in the released DoPHOT PSF photometric catalog (Alonso-García et al. 2018).As mentioned in Sect.3.2.1,our DaoPHOT photometry is affected by the noise patterns introduced by bilinear interpolation.Figure 5 shows J − H color maps made from our DaoPHOT catalog before the photometric calibration and from the released DoPHOT catalog.We can see obvious Moiré patterns in the J − H color map of our DaoPHOT photometry.There are no Moiré patterns in the DoPHOT-based color map.This can be understood by considering the differences in our data treatment.Alonso-García et al. (2018) performed DoPHOT PSF photometry on every stacked SDI, which are also affected by the noise patterns.However, using the photometric uncertainties as weights, these latter authors averaged the brightnesses of the sources detected in more than one stacked SDI.This averaging process dealt with the noise patterns.Since the noise patterns are not present in the DoPHOT catalog, we decided to use it to calibrate our DaoPHOT photometry and remove, to a certain extent, the Moiré patterns.The following general transformation equations were used in this process: In these equations, z, c, and b are the zeropoints, coefficients of color terms, and coefficients of functions that are related to the positions, respectively.P(x, y) describes the magnitude or color difference between the DoPHOT and DaoPHOT photometries as a function of the image position and therefore also traces the Moiré patterns.P(x, y) is constructed from a smoothed differential between the source magnitudes, or colors, in the DoPHOT and DaoPHOT photometries.For example, P J (x, y) is obtained by smoothing the [J DoP − J DaoP ] of sources using a Gaussian kernel with the FWHM of 10 -30 .The value of the FWHM depends on the source number density of the field; we require at least five sources within each beam to facilitate the smoothing.In practice, the calibration proceeds as follows for each tile.(1) We crossmatch our DaoPHOT photometry with DoPHOT catalog using a tolerance of 0.34 (∼1 pixel).Considering the standard deviations of J DoP − J DaoP , H DoP − H DaoP , and sources with (2) We construct P(x, y) with the matched sources.Figure 6 shows the P(x, y) constructed in tile "b214".(3) We select the matched unsaturated sources with uncertainties of < 0.05 at J, H, and K s .Using these high-quality matched sources, a linear regression is performed to solve the transformation Eqs. ( 1)-( 3).(4) We apply the obtained zeropoints and coefficients to our DaoPHOT catalog to obtain the calibrated photometric results.We repeat the above four steps 5-10 times as long as the Moiré patterns keep decreasing and then obtain the final photometry.To illustrate the result, Fig. 5 shows the J − H color distribution of sources in our DaoPHOT catalog after the iterative photometric calibration.The Moiré patterns in the J − H color space have been almost completely removed.
Figure 7 shows a comparison between DoPHOT photometry and our calibrated DaoPHOT photometry for sources in the tiles "d003" and "b305".They match at the bright end, but there are systematic differences (∼0.1-0.3 mag) at the faint end.DoPHOT tends to yield brighter photometry than DaoPHOT for the faint sources.Considering that the average source number density in the bulge tile "b305" is higher than that in the disk tile "d003", it seems that this systematic difference is larger for more crowded fields.This systematic difference between DoPHOT and DaoPHOT for faint stars has been previously seen in studies of stellar clusters (Friel & Geisler 1991;Janes & Heasley 1993;Hill et al. 1998).In particular, Friel & Geisler (1991) compared DoPHOT and DaoPHOT in a cluster field and an uncrowded field.They found that the systematic difference between DoPHOT and DaoPHOT is small in the uncrowded field.However, in the cluster field the difference is a function of position relative to the center of cluster, increasing towards the cluster center.This surroundings of stars without considering the PSF-fitting process.This usually results in higher sky estimations by DaoPHOT than by DoPHOT (Friel & Geisler 1991;Schechter et al. 1993;Hill et al. 1998).Overall, the two photometry packages adopt different philosophies and use different algorithms.As a result, it is difficult to judge which one is "better" without a direct calibration of the faint stars.Adding to this problem, experiments on simulated frames with artificial stars have found that neither DoPHOT nor DaoPHOT can recover the true flux of faint stars in crowded fields (Schechter et al. 1993).Therefore, we here note that there are systematic differences between the released DoPHOT catalog and our DaoPHOT photometry for the faint sources, even though we use the bright sources of the DoPHOT catalog to calibrate our DaoPHOT catalog.In Sect.4.2.2, we revisit this issue and quantify the systematic difference between the DaoPHOT and DoPHOT photometric results through a detailed comparison.
After the calibration, our DaoPHOT catalog is in the VISTA photometric system.To date, the 2MASS survey (Skrutskie et al. 2006) is still the main all-sky near-infrared photometric survey.Many models and data in the literature are based on the 2MASS photometric system.Therefore, to enable convenient comparison with previous works in literature, we decided to provide our catalog in both VISTA and 2MASS systems.To make the transfer to the 2MASS system, we re-calculated the transformations ourselves for each tile using the method described in Soto et al. (2013).Figure 8 shows the comparison between 2MASS and our DaoPHOT photometric catalog in tiles "d003" and "b305".We also adopt a magnitude limit in the bright end of our DaoPHOT catalog in order to remove the sources that are saturated or affected by residual nonlinearity.This limit is (J, H, K s ) = (13.8,12.8, 12.8) mag in VVV disk tiles and (13.0, 11.8, 11.8) mag in bulge tiles.The sources brighter than the limit in any band are replaced by 2MASS photometry.We note that the method suggested by Soto et al. (2013) does not consider the effect of interstellar extinction.The transformations calculated without the interstellar extinction correction could introduce additional bias (González-Fernández et al. 2018).Therefore, in our catalog, we provide the photometry in both the VISTA and 2MASS systems.The users can use their own transformations to calculate the magnitudes of sources in other photometric systems.

Flagging of spurious detections
Saturated sources usually cause a number of spurious detections around them.To flag the spurious detections, the following two steps are performed.We firstly estimate the mean sky level (skymode) for each piece image and define that an image pixel is affected by saturation if its pixel value is higher than a threshold of skymode + 0.8 × (sat-skymode), A85, page 7 of 19 where sat is the saturation level of the detector as given in the VIRCAM user manual7 .We also estimate the mean FWHM of all stars in each piece image and then search pixels around each star within the mean FWHM radius.If there are any surrounding pixels affected by saturation, the star is marked as a source affected by saturation.The second step is the source flagging based on the photometric uncertainties.Mauro et al. (2013) suggested that the spurious detections can be rejected with the sigma-clipping method based on the magnitude-error relation.Figure 9 illustrates the procedure.For the photometric uncertainties in m 1 band and magnitudes in m 2 band, we use an exponential function to fit the σ m 1 − m 2 relation: The corresponding threshold can be obtained as Any source that has a magnitude of m2 and an error greater than σ thres m1 (m2) is flagged as a spurious detection.For the three bands, that is, J, H, and K s , there are nine possible combinations.We note that we only fit the relations of σ J − J, σ H − H, and σ K s − K s and obtain σ thres J (J), σ thres H (H), and σ thres K s (K s ).For the other combinations we simply assume Figure 10 shows the spatial distribution of the spurious detections in a partial region of the tile "d003".We also mark the 2MASS sources that have J < 13.8 mag or H < 12.8 mag or K s < 12.8 mag.We found that the flagged spurious detections are mainly the false detections around the bright sources that are affected by saturation or nonlinearity.There are also some very faint sources with large photometric uncertainties.We note that the flagging procedure inevitably flags some real detections as spurious sources.We found in Fig. 10 that some real detections that are erroneously flagged as spurious detections are sources with significant differences in the photometric uncertainties among the different bands.Overall, the process flags ∼5−27% of the sources in the different tiles as spurious detections with a median value of ∼10%.In total, about 73 million sources are flagged as spurious detections.
A85, page 8 of 19 The process above mainly flags sources with relatively large uncertainties.Because we combine the multi-epoch data and perform the photometry on the combined images, we have lost the flux variability of the sources and the variable stars could have large photometric uncertainties.Therefore, the variable stars could be flagged as spurious detections in our DaoPHOT catalog.To look into the magnitude of this effect, we crossmatched the flagged spurious detections with the International Variable Star Index (VSX8 ) catalog (Watson et al. 2006) that is an up-todate database collecting the known and suspected variable stars from the literature.To account for the large astrometric uncertainties of sources in some of the older publications, we also use a large matching tolerance of 10 (Drake et al. 2013).Finally, 127 555 (<0.2%) variable stars are identified in our flagged spurious detections.Of course, the VSX catalog is not a uniform catalog.The future systematic and uniform searches of variables with the VVV data will be helpful to quantify the fraction of variable stars among our flagged spurious detections (Medina et al. 2018).

Results and discussion: PSF photometric catalog
With our pipeline, we detect about 926 million sources in the VVV survey area and obtain their photometry.About 342 million sources are in the disk area and about 584 million in the bulge area.We note that we only keep the sources that have detections in all three of the J, H, and K s bands.If we exclude the sources flagged as spurious detections (see Sect.   .Spatial distribution of spurious detections in a partial region (∼50 × 50 ) of tile "d003".The background is part of a deep SDI at K s and the spurious detections are marked with red circles.We also mark the 2MASS sources that have J < 13.8 mag or H < 12.8 mag or K s < 12.8 mag with cyan pluses.
Our final catalog contains 29 parameters, including the WCS coordinates, photometric magnitudes and uncertainties, and the quality control parameters.An overview of the available columns of the catalog is given in Table 2.The provided coordinates ([RAJ2000, DEJ2000]) are in the equatorial system.The photometric magnitudes ([J_mag_2MASS, H_mag_2MASS, K s _mag_2MASS]) are in the 2MASS system while the magnitudes ([J_mag_VISTA, H_mag_VISTA, K s _mag_VISTA]) are A85, page 9 of 19 in the catalog can be used to filter out the possible spurious detections (see Sect. 3.2.4).In this paper, a spurious detection is defined as a source for which any of the nine flags has a value of 1. Obviously, this is a conservative definition; the potential users can also adopt their own definitions, for example only using [JJ_flag, HH_flag, K s K s _flag] equal 1 to define the spurious detections.
In the following, we first describe the limiting magnitudes and completeness of our catalog and then compare it with the previously released PSF photometric catalog by Alonso-García et al. (2018).We emphasize that the limiting magnitudes and completeness of our catalog are calculated based on the photometric magnitudes in the 2MASS system while the comparison with the DoPHOT catalog is based on the photometric magnitudes in the VISTA system.

Limiting magnitudes and completeness
The 5σ limiting magnitudes are estimated using the sources with photometric uncertainties of σ < 0.2 mag after excluding spurious detections.We split each tile into many ∼10 × 10 subregions and adopt the magnitudes of the faintest source with σ < 0.2 mag as the 5σ limiting magnitudes of each subregion.The 5σ limiting magnitudes of the tile are then adopted as the median values of the limiting magnitudes of all subregions.Figures 11-13 (top panels) show the 5σ magnitude limits of all VVV tiles.Due to the different exposure time, image quality, and crowding of the fields, the 5σ magnitude limits in J, H, and K s bands are in the ranges of 19.6-21.5, 18.5-20.1, and 18.1-19.3 mag with median values of 20.8, 19.5, and 18.7 mag, respectively.To provide a first-order measurement of completeness at each VVV tile, we constructed single-band brightness distributions in J, H, and K s bands for each tile using sources after excluding spurious detections.In practice, we also split each tile into many ∼10 × 10 subregions and make histograms of magnitudes with a bin size of 0.2 mag in each filter band for each subregion.The completeness of each subregion is estimated roughly as the magnitude of the histogram peak.The completeness of each tile is obtained by averaging the completeness values of all subregions.Figures 11-13 (bottom panels) show the completeness of each VVV tile in J, H, and K s .This results in the completeness limits of our DaoPHOT photometric catalog of J ∼ 17.5−19.9,H ∼ 15.8-18.7,and K s ∼ 14.4-18.2mag.For the VVV disk area, the mean value of completeness is K s ∼ 17.6 mag with a standard deviation of 0.2 mag.However, for the VVV bulge area, the completeness is highly variable.The outer bulge region has the completeness of K s ∼ 17.7 mag while the completeness of inner bulge region drops from K s ∼ 17 mag to only 14.4 mag.The mean value of completeness for the whole VVV bulge area is K s ∼ 17.4 mag with a standard deviation of 0.6 mag.

Comparison with the released VVV PSF photometry
In Sect.3.2.3we present a comparison of the brightnesses between our DaoPHOT PSF photometry and the DoPHOT photometry released in Alonso-García et al. (2018).Here, we additionally compare the photometries by examining the resulting luminosity functions, photometric difference, color-magnitude diagrams (CMDs), and the color excess maps.We perform the comparison in six ∼10 × 10 regions that are located in the disk and bulge areas, cover different Galactic latitudes, and thus represent regions with different stellar fields and source number densities.

Luminosity distribution
Because we performed DaoPHOT photometry on the deep SDIs that have longer total exposure times, we expect to detect more faint sources than is included in the released DoPHOT photometric catalogs.Figure 14 shows the K s band luminosity distributions for the six test fields.We show the K s magnitude distributions of all detected DaoPHOT sources and of the highquality (hq) DaoPHOT sources after excluding the spurious detections (see Sect. 3.1).Our DaoPHOT photometric catalog includes more faint sources and reaches ∼1 mag deeper than the released DoPHOT catalog.As an example, Fig. 15 shows the sources detected both in our DaoPHOT catalog and the released DoPHOT catalog in a small region of a deep SDI in one test field (region 5).The figure also shows the DaoPHOTonly sources, that is, those detected only in our DaoPHOT catalog.
We note that in some regions in the magnitude range of K s ∼ 16-17 mag the DoPHOT catalog includes more sources than our DaoPHOT catalog (e.g., region 1, 3, 4, and 5).To look further into this, we cross-matched the released DoPHOT catalog with our DaoPHOT catalog in each test field and identified "DoPHOT-only" sources.We found that most of the DoPHOTonly sources are rejected DaoPHOT sources which failed to be fitted by the DaoPHOT PSF models.For example, in region 5 there are 76 999 DoPHOT and 75 155 DaoPHOT sources, and 21 460 (28%) DoPHOT sources have no DaoPHOT counterparts.Of these 21 460 DoPHOT-only sources, 15 241 (71%) sources are in the DaoPHOT rejected-source catalog and are therefore removed as spurious detections.Figure 16 (right panel) shows the distribution of some of these 15241 sources on the K s deep SDI.However, there are 6219 (8%) DoPHOT sources that are neither in the DaoPHOT source catalog nor in the DaoPHOT rejected-source catalog.We examined these sources visually on the images.Figure 16 (left panel) shows some of A85, page 10 of 19 the 6219 sources on the K s deep SDI.These sources appear to be very faint and/or highly blended with other sources.Therefore, we suspect that most of them are also spurious detections.These discrepancies could be due to the different algorithms used in DoPHOT and DaoPHOT, as pointed out by Schechter et al. (1993).DaoPHOT detects objects using a threshold of signalto-noise ratio (S/N) and then applies PSF-fitting to the detected objects.During the PSF-fitting process, some objects could fail the fitting and be rejected.In contrast, DoPHOT detects objects using a strict S/N threshold and fits all detected objects without rejection.Therefore, if a low S/N threshold is adopted, DaoPHOT is likely to detect objects without a well-defined S/N limit, but DoPHOT is inclined to include more spurious detections.We conclude that our DaoPHOT photometric catalog can reach about 1 mag deeper and includes less spurious detections than the released DoPHOT catalog.

Photometric differences
In Sect.3.2.3,we compared the photometric magnitudes between our DaoPHOT catalog and the released DoPHOT catalog and found a systematic difference at the faint end.In this section, we investigate these differences further in six test fields that represent different crowding conditions.We emphasize that in this and the subsequent sections (Sects.4.2.3 and 4.2.4), we only consider DoPHOT sources that have detections at J, H, and K s simultaneously and DaoPHOT sources that are not spurious detections (see Sect. 3.2.4).We also require that the uncertainties of the DoPHOT and DaoPHOT sources be less than 0.35 mag in the J, H, and K s bands.
We checked the photometric difference between DaoPHOT and DoPHOT source in six test fields and found that there are obvious systematic differences at the faint ends in all six test fields.To quantify the systematic differences, Fig. 17 shows the moving medians of ∆J, ∆H, and ∆K s between DaoPHOT and DoPHOT sources as functions of J, H, and K s magnitudes, respectively.In the Galactic disk area (regions 1, 2, and 3), the systematic differences appear at J ∼ 17, H ∼ 16, and K s ∼ 15.5 mag and can reach up to ∆J ∼ 0.2, ∆H ∼ 0.2, and ∆K s ∼ 0.3 mag.In region 4, which is located in the outer bulge, the systematic differences appear at J ∼ 16, H ∼ 15.5, and K s ∼ 15 mag and reach up to ∆J ∼ 0.4, ∆H ∼ 0.3, and ∆K s ∼ 0.3 mag at J ∼ 18.5, H ∼ 18, and K s ∼ 17.5 mag, respectively.In region 5, located in the inner bulge, the systematic differences appear at J ∼ 15, H ∼ 14, and K s ∼ 14 mag and reach up to ∆J ∼ 0.5, ∆H ∼ 0.4, and ∆K s ∼ 0.4 mag at J ∼ 18, H ∼ 17, and K s ∼ 17 mag, respectively.In region 6, located in the Galactic center, the systematic differences appear at J ∼ 17, H ∼ 15, and K s ∼ 13.5 mag and reach up to ∆J ∼ 0.15, ∆H ∼ 0.25, and ∆K s ∼ 0.3 mag at J ∼ 19, H ∼ 17.5, and K s ∼ 16.5 mag, respectively.Therefore, the differences depend on the location in the Galactic disk and bulge, caused mainly by the local crowding conditions.We also note that the source number density is different in different filter bands in one location.For example, the average source number density of region 6 is similar to that of region 5 at K s band, but significantly smaller than that of region 5 at J band due to the heavy interstellar extinction towards region 6.
We also investigate the photometric colors between DaoPHOT and DoPHOT catalogs.There are systematic color A85, page 11 of 19 differences at the faint end, especially in the crowded fields such as regions 5 and 6 which are located close to the Galactic center.Figure 18 shows the moving medians of ∆JK s , ∆JH, and ∆HK s between DoPHOT and DaoPHOT catalogs as functions of K s magnitudes in six test fields.In regions 1, 2, and 3, there are no significant differences in the [J − H] color.However, there are differences in [H − K s ] and [J − K s ] colors that appear at K s ∼ 16 mag and can reach up to ∼0.1 mag at K s ∼ 17.5 mag.In region 4, there is no systematic difference for any of the three colors.In region 5, the differences in all [J − K s ], [J − H], and [H − K s ] appear at K s ∼ 15 mag and can reach up to ∼0.1, ∼0.05, and ∼0.05 mag at K s ∼17 mag, respectively.In region 6, differences in all colors appear at K s ∼ 14.5 mag and can reach up to ∼0.2, ∼0.1, and ∼0.1 mag at K s ∼ 16-17 mag, respectively.
Overall, the systematic photometric magnitude difference between our DaoPHOT catalog and the released DoPHOT catalog is different in the different locations of the Galactic disk and bulge, mainly depending on the local source number density.The systematic color difference is less than 0.1 mag and appears at K s ∼ 16 mag in the Galactic disk area; it reaches up to 0.2 mag and appears at K s ∼ 14-15 mag in the bulge area.In general, the systematic difference of colors between DaoPHOT and DoPHOT catalogs is smaller than that of brightnesses.

Detectable features in the CMDs
The CMD is a powerful tool to identify the stellar populations in the Galactic disk and bulge (Saito et al. 2012b;Soto et al. 2013;Alonso-García et al. 2018).Figures 19 and 20 compare the CMDs constructed with our DaoPHOT catalog and the released DoPHOT catalog in six test fields.We can see obvious features in the CMDs.To understand the stellar populations corresponding to these detectable features, we obtained synthetic CMDs for the fields using the Besançon Galactic model (Robin et al. 2003).We only show the most populated sequences with different color symbols, that is, the main sequence (MS) disk stars and the giants and subgiants in the disk and bulge.
Figure 19 shows the CMDs in regions 1, 2, and 3.In region 1, there are two main features in the CMD.The feature that extends from K s ∼ 13 and [J − K s ] ∼ 0.3 down to K s ∼ 19.5 and [J − K s ] ∼ 2 corresponds to the disk MS and disk subgiant population.The other feature that occupies the color space of [J − K s ] ∼ 1.4 and K s ∼ 14 corresponds to the disk giants.In region 2 and 3, there are three main features.The features with relative blue ([J −K s ] ∼ 1-2) and red ([J −K s ] ∼ 3-5) colors correspond to disk MS and disk giants, respectively.There is also another feature between them, extending from [J − K s ] ∼ 1.5 and K s ∼ 14 down to [J − K s ] ∼ 4 and K s ∼ 18.We note that most of disk subgiants are located in this region in the synthetic CMDs.
Figure 20 shows the CMDs for regions 4, 5, and 6.The CMDs of regions 4 and 5 have similar patterns.The different stellar populations are mixed up in the faint end, but there are two distinct branches in the bright end.The blue branch that extends from [J − K s ] ∼ 0.3 and K s ∼ 13 down to [J − K s ] ∼ 0.7-1 and K s ∼ 16 mainly corresponds to the disk MS stars.The red branch that extends from [J − K s ] ∼ 1-1.5 and K s ∼ 12 down to [J − K s ] ∼ 0.8-1.2 and K s ∼ 16 consists of disk giants ( 20%) and bulge giants ( 70%).When investigating the CMD of region 6, the features become significantly blurred, which is A85, page 12 of 19 due to the heavy interstellar extinction.Compared with the synthetic CMD, it seems that the different stellar populations such as disk MS stars, disk subgiants, and bulge giants can be roughly separated from one another.There is also an overdensity in the CMD of region 6 at [J − K s ] ∼ 4.5, K s ∼ 14.5.Based on the synthetic CMD it mainly consists of disk giants (∼20%) and bulge giants (∼80%).
The CMDs constructed with DaoPHOT and DoPHOT are relatively similar.To quantify the differences between them, we show the differential CMDs in Figs.19 and 20.These were obtained by subtracting the number density CMD map built with DoPHOT from that built with DaoPHOT.Indeed, because our DaoPHOT photometry can reach higher limiting magnitudes than the released DoPHOT photometry, there are positive overdensities at the faint ends of the differential CMDs.There are also negative overdensities at K s ∼ 16-17 mag in some differential CMDs, which means that the DoPHOT catalog includes more sources than the DaoPHOT catalog in this range.We already investigated the excess in Sect.4.2.1 and suggested that it results from the DoPHOT catalog including more spurious detections than our DaoPHOT catalog.Except for the above two features, there are no other significant patterns in the differential CMDs.Therefore, there is no significant difference between the CMDs constructed with DaoPHOT and DoPHOT catalogs in the range K s 16 mag.

Color excess maps
We plan to use our DaoPHOT catalogs in the future to derive extinction maps for molecular clouds in the Galactic plane.For this, it is important to investigate the difference of extinction maps constructed with DaoPHOT and DoPHOT catalogs.However, performing extinction mapping for molecular clouds is not a straightforward process, especially when clouds are located at distances larger than approximately 1 kpc.This is mainly due to the contamination by foreground sources that are located between the clouds and the observer (e.g., Kainulainen et al. 2011).Therefore, in this section we only compare the color excess maps derived with DaoPHOT and DoPHOT catalogs in six ∼1 • × 1 • test fields.To keep the test simple, we do not try to isolate molecular clouds or exclude any foreground contamination.The six ∼1 • × 1 • test fields have the same centers as the six previously used ∼10 × 10 test fields.
The color excess of each star can be obtained with where [J − K s ] intrinsic is the intrinsic color that is usually estimated as the average color of stars in a nearby control field that is free from extinction.In practice, we first adopt an arbitrary value as the intrinsic color and then calculate the color excess of each star.The color excess maps are obtained by smoothing the color excesses with a Gaussian kernel of FWHM = 30 in width.Based on these maps, we select one control field that has relatively low color excess; Fig. 21 shows the locations of the control fields.Second, we use the average [J − K s ] color of the sources in each control field as the intrinsic color and then obtain the final color excess maps.We note that it is very difficult to select the control fields that are free from extinction in the Galactic plane.However, our purpose here is to investigate the difference of the color excess maps built with DaoPHOT and DoPHOT catalogs, not to study extinction itself.Contamination of the control field does not affect the relative values of color excesses derived from the DaoPHOT and DoPHOT catalogs.
We produce color excess maps using both DaoPHOT and DoPHOT catalogs (see Fig. 21 for the DaoPHOT maps).To compare the maps, we obtain differential color excess maps by subtracting the DaoPHOT-based map from the DoPHOT-based map. Figure 22 shows the differential maps.The differential maps show patterns close to the background rms of DaoPHOT and DoPHOT color excess maps.The median values and the standard deviations of the six differential maps are −0.02± 0.04, −0.01 ± 0.07, −0.02 ± 0.09, −0.01 ± 0.03, 0.03 ± 0.03, and 0.06 ± 0.18.Thus, there is no significant systematic difference between the DaoPHOT and DoPHOT color excess maps for the test fields except region 6 which is one of the most crowded fields in the Galactic bulge.In region 6, the systematic difference between the DaoPHOT and DoPHOT color excess maps is about 0.06, corresponding to the visual extinction of ∼0.35 mag if using the extinction law suggested by Rieke & Lebofsky (1985).The above analysis is based on color excess maps made at a spatial resolution of 30 .We repeated the analysis with maps made at the resolutions of 60 and 15 .The change of spatial resolution does not affect the aforementioned conclusions.

Conclusions of the comparisons
Overall, there are systematic differences in the photometric magnitudes between DaoPHOT and DoPHOT catalogs beginning at K s ∼ 14-16 mag.These differences can reach up to ∆K s ∼ 0.1-0.4mag at K s ∼ 17 mag, depending on the crowding of the fields.However, the systematic photometric color differences are less than 0.1 mag except for the most crowded fields close to the Galactic center.Moreover, there is no significant difference between the color excess maps built with DaoPHOT and DoPHOT sources in the uncrowded fields.In very crowded fields, the difference in the colors of the two catalogs can be up to 0.2 mag and result in the systematic offset of ∼0.06 (∆A V ∼ 0.35 mag) between the color excess maps based on them.
In principle, the DaoPHOT algorithm is more suitable for crowded fields than the DoPHOT algorithm, crowded fields (Schechter et al. 1993;Becker et al. 2007).Thus, our DaoPHOT catalog should perform better in crowded fields than the released DoPHOT catalog.However, we must emphasize that our DaoPHOT catalog is obtained based on common settings for general photometric parameters; it is difficult to set common parameters for fields with different crowding conditions.Potential users can perform their own PSF photometry using dedicated parameter configurations with our PSF-fitting pipeline should they find that our catalog does not satisfy their scientific requirements.We will release our pipeline to the community in the near future.

Summary and conclusions
In this paper, we present a new PSF photometric catalog based on the DaoPHOT algorithm for the VVV survey data.We also compare the catalog to the PSF photometric catalog recently released by Alonso-García et al. (2018).The main results are summarized as follows: -We developed an automatic PSF-fitting pipeline based on the DaoPHOT algorithm.To reach the highest possible limiting magnitudes, we stack the multi-epoch VVV images and perform photometry on the stacked images.The pipeline mainly uses the PyRAF package and is adapted to run in a multi-core mode.We expect to use the pipeline for other survey data in the future.-Our catalog contains about 926 million sources with detections at J, H, and K s bands.About 10% of the sources are flagged as possible spurious detections.Based on the sources that are not flagged, we estimate the 5σ limiting magnitudes to be about 20.8, 19.5, and 18.7 at J, H, and K s , respectively.Using the peak of the brightness distributions as a first-order measure of completeness (after excluding spurious detections), we estimate the completeness limits of our catalog to be about 19.0, 18.0, and 17.5 mag in J, H, and K s bands, respectively.-The comparison of our DaoPHOT catalog with the previously released DoPHOT catalog suggests that on average A85, page 17 of 19 Developing accurate and efficient photometric pipelines is necessary to make use of the approximately 1 billion sources detected by the VVV survey.This is not a trivial task and the properties of the outcome depend on the chosen approach.Our deep PSF photometric catalog can be exploited for a multitude of scientific purposes.In the near future, we plan to use it to derive a new dust extinction map for molecular clouds in the Galactic plane and to investigate the internal structure of the clouds.

Fig. 1 .
Fig. 1.Sketch of "Tile6n" observing pattern for one VISTA detector.The six dithering pointings are labeled with the red filled circles.The numerals marked on red circles and the arrows show the dithering sequence.The filled green box represents the field of view of a detector and the solid blue square marks the total area covered by six dithering images.

Fig. 5 .Fig. 6 .
Fig. 5. J − H color images of sources in tile "b214" of our DaoPHOT catalog before the photometric calibration (left), of the released DoPHOT catalog (middle, Alonso-García et al. 2018), and of our DaoPHOT catalog after the photometric calibration (right).The images are produced by smoothing the source J − H colors with a Gaussian kernel (FWHM = 30 ).

Fig. 7 .
Fig. 7. Released DoPHOT photometry vs. our DaoPHOT photometry in tiles "d003" (left panels) and "b305" (right panels).The backgrounds are source number density maps.The green curves show moving medians with a box width of 0.2 mag while the red lines mark the zero difference value.

Fig. 8 .
Fig.8.2MASS vs. our DaoPHOT photometry in tiles "d003" (left panels) and "b305" (right panels).The backgrounds show source number density.The green curves show moving medians with a box width of 0.2 mag while the red horizontal lines mark the zero difference value between the 2MASS and our DaoPHOT catalogs.The blue vertical lines show the magnitude limits of (J, H, K s ) = (13.8,12.8, 12.8) mag and (13.0, 11.8, 11.8) mag for VVV disk tiles and bulge tiles, individually.The source brighter than above limits have been replaced by 2MASS photometric results.

Fig. 9 .
Fig. 9. Relations between photometric uncertainties and magnitudes for sources in a partial region of "d003" (∼3 × 3 central region of tile "d003").The gray dots represent all detected sources.The flagged spurious detections are marked with green dots.The blue curves show the exponential fittings for all the sources while the red curves mark the corresponding thresholds of photometric uncertainties that are used to flag the possible spurious detections.
3.1), there are about 853 million sources remaining, from which about 319 and 534 million are in the disk and bulge area, respectively.As a comparison, the released DoPHOT catalog includes about 180 and 354 million sources with J, H, K detections in the disk and bulge area.

Fig. 10
Fig.10.Spatial distribution of spurious detections in a partial region (∼50 × 50 ) of tile "d003".The background is part of a deep SDI at K s and the spurious detections are marked with red circles.We also mark the 2MASS sources that have J < 13.8 mag or H < 12.8 mag or K s < 12.8 mag with cyan pluses.

Fig. 11 .
Fig. 11.5σ magnitude limits (top panel) and completeness limits (bottom panel) for each VVV tile in the J band.

Fig. 12 .
Fig. 12. 5σ magnitude limits (top panel) and completeness limits (bottom panel) for each VVV tile in the H band.

Fig. 13 .
Fig. 13.5σ magnitude limits (top panel) and completeness limits (bottom panel) for each VVV tile in the K s band.

Fig. 14 .Fig. 15 .Fig. 17 .
Fig. 14.Luminosity distributions in the K s band in six ∼10 × 10 regions.All detected DaoPHOT sources are shown with black solid lines while the hq DaoPHOT sources after removing possible spurious detections are labeled with red dotted lines.The DoPHOT sources retrieved from the released PSF photometric catalog (Alonso-García et al. 2018) are marked with blue dash-dotted lines.

Fig. 18 .Fig
Fig. 18.Moving medians of the photometric color difference of ∆JK s (left panel), ∆JH (middle panel), and ∆HK s (right panel) between DaoPHOT and DoPHOT catalogs as a function of K s band magnitude for the six ∼10 × 10 test fields.The different regions are shown with different color curves.

Table 1 .
Exposure times of VVV OBs.

Table 2 .
Available data columns of our DaoPHOT photometric catalog for the VVV survey Notes.The full catalog is available at the CDS. in the VISTA system.All the photometric magnitudes are in the Vega system.The parameters N iter ([J_iter, H_iter, K s _iter]), sharpness ([J_sharpness, H_sharpness, K s _sharpness]), and χ 2 ([J_chi, H_chi, K s _chi]) describe the PSF-fitting quality.The N iter values give the number of iterations used to fit the stars.
The χ 2 refers to the PSF-fitting goodness: χ 2 = 1 means a perfect fit and a poor fit has a value far from 1.The sharpness is a shape measurement of source detections on the image: sharpness=0 means that the detected source is round.The nine spurious detection flags ([JJ_flag, HH_flag, K s K s _flag, JH_flag, JK s _flag, HJ_flag, HK s _flag, K s J_flag, K s H_flag])