New 2MASS Near-infrared Photometry for Globular Clusters in M31
  We present 2MASS $JHK_{\rm s}$ photometry for 913 star clusters and
candidates in the field of M31, which are selected from the latest Revised
Bologna Catalog of M31 globular clusters (GCs) and candidates. The photometric
measurements in this paper supplement this catalog, and provide a most
comprehensive and homogeneous photometric catalog for M31 GCs in the $JHK_{\rm
s}$ bandpasses. In general, our photometry is consistent with previous
measurements. The globular cluster luminosity function (GCLF) peaks for the
confirmed GCs derived by fitting a $t_5$ distribution using maximum likelihood
method are: $J_0 = 15.348_{-0.208}^{+0.206}$, $H_0 = 14.703_{-0.180}^{+0.176}$,
and ${K_{\rm s}}_0 = 14.534_{-0.146}^{+0.142}$, all of which agree well with
previous studies. The GCLFs are different between metal-rich (MR) and
metal-poor (MP), inner and outer subpopulations, as that MP clusters are
fainter than their MR counterparts, and the inner clusters are brighter than
the outer ones, which confirm previous results. The NIR colors of the GC
candidates are on average redder than those of the confirmed GCs, which lead to
an obscure bimodal distribution of the color indices. The relation of
$(V-K_{\rm s})_0$ and metallicity shows a notable departure from linearity,
with a shallower slope towards the redder end. The color-magnitude diagram
(CMD) and color-color diagram show that many GC candidates are located out of
the evolutionary tracks, suggesting that some of them may be false M31 GC
candidates. The CMD also shows that the initial mass function of M31 GCs covers
a large range, and the majority of the clusters have initial masses between
$10^3$ and $10^6$ $M_{\odot}$.
