{'fraction_non_alphanumeric': 0.04179760282544009, 'fraction_numerical': 0.03865513501487806, 'mean_word_length': 4.264236568584395, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A Radio-Selected Black Hole X-ray Binary Candidate in the Milky Way Globular Cluster M62

We report the discovery of a candidate stellar-mass black hole in the Milky Way globular cluster M62. We detected the black hole candidate, which we term M62-VLA1, in the core of the cluster using deep radio continuum imaging from the Karl G. Jansky Very Large Array. M62-VLA1 is a faint source, with a flux density of 18.7 +/- 1.9 microJy at 6.2 GHz and a flat radio spectrum (alpha=-0.24 +/- 0.42, for S_nu = nu^alpha). M62 is the second Milky Way cluster with a candidate stellar-mass black hole; unlike the two candidate black holes previously found in the cluster M22, M62-VLA1 is associated with a Chandra X-ray source, supporting its identification as a black hole X-ray binary. Measurements of its radio and X-ray luminosity, while not simultaneous, place M62-VLA1 squarely on the well-established radio--X-ray correlation for stellar-mass black holes. In archival Hubble Space Telescope imaging, M62-VLA1 is coincident with a star near the lower red giant branch. This possible optical counterpart shows a blue excess, H alpha emission, and optical variability. The radio, X-ray, and optical properties of M62-VLA1 are very similar to those for V404 Cyg, one of the best-studied quiescent stellar-mass black holes. We cannot yet rule out alternative scenarios for the radio source, such as a flaring neutron star or background galaxy; future observations are necessary to determine whether M62-VLA1 is indeed an accreting stellar-mass black hole.


INTRODUCTION
For a typical globular cluster (GC) with a present day mass of ∼10 5 -10 6 M ⊙ , hundreds of stellar-mass black holes (BHs) should be born during the first ∼10 Myr after formation (Larson 1984;Kulkarni et al. 1993). These BHs have masses ∼5-20 M ⊙ and mark the endpoints of stars with main-sequence masses 25 M ⊙ (Heger et al. 2003). The specific frequency of X-ray binaries in GCs is ∼100 times larger than the field and provides strong evidence that mass-transferring binaries are dynamically formed with high efficiency in GCs (e.g., Kundu et al. 2002;Pooley et al. 2003). If BHs are present in GCs, a subset should be detectable as accreting binaries.
Yet, soon after the detection of luminous X-ray sources in Milky Way GCs, it became clear that accreting neutron stars, not BHs, dominate luminous cluster X-ray binaries. The principal evidence is the detection of Type I X-ray bursts that are well-explained by thermonuclear runaway events on the surfaces of neutron stars (Lewin & Joss 1981). Of the 18 luminous Galactic GC X-ray sources published so far, all are thought to contain neutron stars, due to the presence of X-ray bursts (in 15; see e.g. Verbunt & Lewin 2006;Altamirano et al. 2012), coherent pulsations (in 4; e.g. Altamirano et al. 2010), optical spectroscopy (van Zyl et al. 2004), or a low radioto-X-ray flux ratio (Bozzo et al. 2011, see below for the rationale).
Theoretical expectations for the presence of BHs in GCs remain uncertain. After birth, assuming the BHs receive small natal kicks, they will rapidly mass segregate to the cluster center in a fraction of a relaxation time and form a subcluster that is dynamically segregated from the rest of the GC. Many BH-BH binaries will then be formed through three-body interactions, and these binaries halt the collapse of the subcluster through interactions with one another. This process will tend to harden the binaries and lead to interactions with large recoil velocities, ejecting many BHs from the GC. Once a sufficient number of BHs have been ejected, the subcluster will no longer be dynamically decoupled from the stars in the core of the cluster, and the efficiency of BH ejection will decrease.
The eventual fate of the remaining BHs is poorly understood. Many papers have suggested that GCs should have at most one BH (or BH-BH binary) remaining at the present day, and that many clusters should have none (Kulkarni et al. 1993;Sigurdsson & Hernquist 1993;Portegies Zwart & McMillan 2000;Kalogera et al. 2004). Some recent papers have argued a few BHs may be retained (Mackey et al. 2008;Moody & Sigurdsson 2009;Aarseth 2012), while two new theory papers suggest that in some circumstances many BHs could survive in GCs (Sippel & Hurley 2013;Morscher et al. 2013). In sum, recent work supports the idea that BH ejection in GCs is less efficient than previously thought.
Over the past ∼6 years, there have been solid BH candidates identified in GCs around other galaxies (e.g., Irwin et al. 2010). These sources are characterized as likely BH binaries because: (i) they have X-ray luminosities far above the Eddington luminosity for a single neutron star, and (ii) they vary significantly on short timescales, making it implausible that the luminosity originates from a superposition of several neutron star X-ray binaries. The most luminous of these objects is associated with a GC in the massive Virgo elliptical NGC 4472, and has a peak L X ∼ 4 × 10 39 erg/s . Optical spectroscopy of the associated GC shows broad (1500 km s −1 ) [O III] emission but no Balmer lines, and the combined data is best explained by a model in which a stellar-mass BH is accreting at a super-Eddington rate from a CO white dwarf (Zepf et al. 2008;Peacock et al. 2012).
The few observed super-Eddington BHs in extragalactic GCs would then simply be those with the most extreme accretion rates. They likely represent only the very tip of the iceberg in terms of BH binaries in GCs. Many BHs with lower accretion rates are almost certain to exist among X-ray sources in GCs, but they are greatly outnumbered by neutron-star binaries and difficult to distinguish from neutron stars using X-ray data alone.
In Strader et al. (2012a), we developed a new strategy for identifying quiescent BH binaries in Milky Way GCs, making use of both radio and X-ray data (see also Maccarone & Knigge 2007). Stellar-mass BHs accreting at low rates have compact jets which emit radio continuum via partially self-absorbed synchrotron emission (Blandford & Königl 1979). Thus they are much more luminous in the radio than neutron stars with comparable X-ray luminosity: L R /L X is ∼ 2 orders of magnitude higher for BHs than neutron stars (Migliari & Fender 2006). Before the recent upgrade to the VLA, the radio emission from a quiescent BH like A0620-00 or V404 Cyg would not have been detectable at high significance at typical GC distances (Gallo et al. 2006). The upgraded VLA can now readily detect the expected flux densities (tens of µJy) in reasonable exposure times.
Using this technique, we discovered two candidate stellar-mass BHs in the cluster core of M22 (Strader et al. 2012a). The sources have flat radio spectra and 6 GHz flux densities of 55-60 µJy. As these sources are not detected in shallow archival Chandra imaging, they cannot yet be placed directly on the L X − L R relation; nevertheless, their overall properties are consistent with those expected from accreting BH binaries.
Here we report the discovery of a BH candidate in a second Galactic GC, M62 (NGC 6266; D = 6.8 kpc; Harris 1996). We term this source M62-VLA1. Unlike the case for the M22 sources, M62-VLA1 has clear X-ray and optical counterparts, and so is the most compelling candidate black hole X-ray binary in a Milky Way GC.
In Section 2, we describe our VLA observations, along with archival Chandra and HST imaging. In Section 3, we present evidence that M62-VLA1 is an accreting stellar-mass BH. Section 4 discusses the host binary system: the binary separation and binary companion. We assess alternative explanations for M62-VLA1 in Section 5. We summarize our findings in Section 6.

VLA Radio Data
We observed M62 with the Karl G. Jansky Very Large Array (VLA) over the time period 2012 Sept 10-16 as part of the program 12B-073 (P.I. Strader). Ten hours were spent observing the cluster, split among seven blocks of 1-1.75 hour duration, yielding a total of 7 hours on source. We observed in C band with 2 GHz total bandwidth and four polarization products. Of the two basebands of 1 GHz width, one was centered at 5.0 GHz and the other at 7.4 GHz. The array was in BnA configuration, giving a resolution of 1.4 ′′ × 1.1 ′′ at 5.0 GHz. The field of view (full-width at half power) of the VLA at the average frequency of 6.2 GHz was ∼7.3 ′ in diameter, significantly larger than the half-light diameter of M62 (1.8 ′ ; Harris 1996).
We observed J1700-2610 as the secondary phase calibrator and J1407+2827 as the polarization leakage calibrator. 3C286 was used as an absolute flux density, bandpass, and polarization angle calibrator. The data were reduced using standard routines in AIPS. Weights from switched power measurements were applied using TYAPL. Each observing block was edited for bad data and interference and then calibrated. For each individual calibrated baseband, the data were concatenated in the uv plane and then imaged with a Briggs robust value of 1. A bright source at the edge of the field called for phase and amplitude self calibration. Figure 1 shows a deep co-added image of both basebands, obtained by smoothing the 7.4 GHz baseband to the resolution of the 5.0 GHz basebands and averaging these together. The rms sensitivity of this co-added image is 2.0 µJy beam −1 .
Flux density measurements were carried out separately on each baseband's image. We fit Gaussians using JMFIT in AIPS assuming the source is point-like. Measured flux densities for M62-VLA 1 are 19.9±3.2 µJy (5.0 GHz) and 18.1 ± 2.3 µJy (7.4 GHz). Assuming the flux densities follow a power law of the form S ν = ν α , then the spectral index of M62-VLA1 is α = −0.24 ± 0.42. In contrast, PSR J1701-3006A was measured to have a 5.0 GHz flux density of 23.6 ± 2.9 µJy, but was not detected at 7.4 GHz (3.7 ± 2.3 µJy), implying α −2.0. The spectrum of M62-VLA1 is inconsistent with the steep spectrum expected for a pulsar.   (Lützgendorf et al. 2013). The other radio source in the central region of M62 is a known pulsar (red circle; Possenti et al. 2003). North is up and east is to the left.
of X-ray sources, with 51 sources detected within the half-mass cluster radius in the Chandra image and only 2-3 expected to belong to the background (Pooley et al. 2003). However, the properties of the individual sources have not been published. M62-VLA1 is the only X-ray source in M62 which has a counterpart in our VLA image ( Figure 1). It is detected in the Chandra image at high significance, with 688 counts or signal-to-noise ratio of 26.
We fit the X-ray spectrum of M62-VLA in the energy range 0.5-7 keV (experiments with other energy ranges yielded very similar results). The foreground absorption was held constant at N H = 2.8 × 10 21 cm −2 , calculated from the reddening estimate for M62 (E(B − V ) = 0.47; Harris 1996) and using the Güver &Özel (2009) relation between reddening and N H . The spectrum was binned so that there were at least 20 counts per channel, and then fit via χ 2 minimization in XSPEC.
The spectrum was not well-fit by an absorbed blackbody (χ 2 /ν = 5.17 with 29 degrees of freedom). How-ever, we did find a reasonable fit with a combined blackbody and power-law spectrum, as observed for accreting neutron stars. The parameters for this model are listed in Table 1; the power law component accounts for about ∼2/3 of the flux and the best-fit blackbody temperature is 1.9 × 10 6 K. A similarly good fit can be achieved by fitting a model neutron star atmosphere plus a power law, using a NSATMOS model in XSPEC (Heinke et al. 2006) and assuming a neutron star mass of 1.4 M ⊙ and radius of 12 km. The best-fit parameters are again listed in Table 1; the luminosity is practically identical to the power-law + blackbody model, but the temperature of the best-fit neutron star atmosphere is lower.
These fits emphasize the well-known conclusion that it is challenging to distinguish BHs from neutron stars on the basis of X-ray spectra alone.

HST Optical Data
There are archival Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) data that cover the core of M62, obtained in 2004 (two years after the Chandra data discussed in §2.2, but more than eight years before the VLA data were observed). These data comprised 880 sec of exposure time in F 435W (B-equivalent; split into three frames), 1170 sec in F 625W (R-equivalent; five frames), and 3610 sec in F 658N (narrow band covering Hα; ten frames).
We photometered the individual frames using DOLPHOT (Dolphin 2000) before producing a matched photometric catalog. We note that M62 has a large and differential foreground reddening with mean E(B − V ) = 0.47 (Alonso-García et al. 2011); all photometry quoted in this paper is uncorrected for reddening.
To match the astrometry of the VLA images, we iteratively corrected the astrometry of the HST images to the ICRS using a large number of 2MASS sources across the images. The overall uncertainty in the astrometric zero point of the image is about 0.15 ′′ . The uncertainty in the position of M62-VLA1 itself is 0.07 ′′ , giving a combined 1σ uncertainty of 0.17 ′′ in matching the HST and VLA data. Unfortunately there are no additional sources that appear in both the HST and VLA or Chandra images, so further checks on the astrometry are not possible. Figure 2 shows the position of M62-VLA1 superposed on the HST F 435W image. Even given the significant astrometric uncertainty, M62-VLA1 is a very close (< 0.05 ′′ ) match to a moderately bright star in the images. This object has mean F 625W = 17.38 mag, and F 435W − F 625W = 1.41. In a F 625W vs. F 435W − F 625W color-magnitude diagram the star falls near the lower red giant branch of M62, but lies ∼0.4 mag blueward of the cluster fiducial ( Figure 3).
The F 625W − F 658N color of the star is ∼0.08 to 0.1 mag redder than the inferred locus of F 435W − F 625W giants, calculated by connecting the principal population of red giant branch stars to that of horizontal branch stars. This photometry is consistent with the presence of Hα emission with an equivalent width of ∼7Å in the F 658N filter. Figures 3 and 4 show color-magnitude and color-color diagrams with the position of this star marked.
This brighter star is flanked by two fainter stars that also fall just within the nominal astrometric error circle (red region in Figure 2). These other stars have photometry consistent with main sequence stars just below the turnoff of M62.
There is some evidence that the brighter candidate optical counterpart star is variable. The ten frames of HST The mean photometry of the M62-VLA1 candidate optical counterpart discussed in §2.3 is plotted as a red square. The horizontal branch would be located at F 625W 16 mag, but most of these stars are saturated in the F 625W photometry. The candidate optical counterpart of M62-VLA1 (red squares; data at two epochs) is ∼0.09 mag more luminous in F 658N than the main locus of stars. One explanation is that this deviation is due to the presence of Hα emission in the F 658N band with an equivalent width of ∼7Å.
F 658N photometry, acquired over 3 hours, show a linear trend to fainter magnitudes, with amplitude ∼0.06 mag over this time period (see Table 2, which lists all photometry of this star, and Figure 5). To assess the significance of this trend, we examined the photometry of stars of comparable brightness (16.5 < F 658N < 17.5 mag). For these comparison stars, on average, there is no evidence for a systematic trend in the photometry among the frames. If we fit linear relations to the F 658N magnitudes as a function of time for all the comparison stars, the inferred amplitude of the trend for the M62-VLA1 candidate counterpart is larger than 98% of the stars. The three F 435W frames cover a comparable time baseline and show a similar trend, although with a somewhat smaller amplitude. Due to saturation, only two of five frames in F 625W have reliable photometry, but these are consistent with the behavior observed in the other filters.
Because of the crowding in the core of M62 and the relatively modest amplitude of the variability, we consider this finding worthy of further investigation but far from conclusive.

A NEW BLACK HOLE CANDIDATE
The most convincing evidence supporting a BH identification for M62-VLA1 is its ratio of X-ray luminosity (L X ) to radio luminosity (L R ). Figure 6 plots M62-VLA1 on the standard X-ray-radio correlation (Gallo et al. 2006). Its mean flux density of 19 µJy corresponds to an equivalent 8.4 GHz radio luminosity of L R = 9 × 10 27 erg s −1 at the distance of M62, assuming a flat spectrum (α = 0; if the measured α = −0.24 is used instead, then L R = 8 × 10 27 erg s −1 ). A spectral fit to the X-ray data gives a luminosity of L X = 7 × 10 31 erg s −1 over 3-9 keV. This particular radio frequency and Xray spectral range are chosen because they are the most commonly used in L X − L R relations in the literature.
With the important caveat that variability could be present between the X-ray and radio epochs, M62-VLA1 sits squarely on the BH L X − L R relation. The source appears to be a doppelganger for the quiescent BH V404 Cyg in its X-ray and radio luminosities (Miller-Jones et al. 2009). The radio luminosity for M62-VLA1 is 2 orders of magnitude higher than expected for accreting neutron stars or white dwarfs.
Both the radio and X-ray spectra of M62-VLA1 are consistent with an accreting stellar-mass BH. The radio spectral index is consistent with being flat (α = −0.2 ± 0.4), similar to the radio indices of known low-luminosity accreting stellar-mass BHs (α = 0.0-0.2; Fender 2001; Gallo et al. 2005). However, the radio spectral index could also be consistent with a flaring neutron star, which can have a radio spectrum that varies from α ≈ −0.6 as expected for optically-thin synchrotron emission to slightly inverted (α > 0; e.g., Marti et al. 1992;Fender 1997;Moore et al. 2000;Miller-Jones et al. 2010).
As described in Section 2.2, we also cannot rule out a neutron star origin for M62-VLA1 from the X-ray spectra alone. Either a single power-law or a power-law + thermal model is a good match to the X-ray data. The power-law spectral fit for M62-VLA1 (Γ = 2.50 ± 0.14) is comparable to that of the quiescent BH binaries V404 Cyg (Γ = 2.17 ± 0.13; Corbel et al. 2008), A0620-00 (Γ = 2.1 +0.5 −0.4 ; Gallo et al. 2006), and MAXI J1659-152 (Γ = 2.5 ± 0.3; Jonker et al. 2012), although it is slightly softer than commonly observed. The X-ray spectrum of M62-VLA1 is consistent with the claim that the lowluminosity BHs display a softening of their X-ray spectra in quiescence (Corbel et al. 2006;Wu & Gu 2008).
At a projected distance of 3.7 ′′ from the photometric center (0.12 pc at the distance of M62), the separation of M62-VLA1 is too large to be consistent with an intermediate-mass BH (∼10 2 − 10 4 M ⊙ ): such sources would be expected to be 1 ′′ from the center (Chatterjee et al. 2002; Gallo et al. (2006) (L R ∝ L 0.58 X , normalized by a least-squares fit to the simultaneous detections with L X < 2 × 10 34 erg s −1 ). The blue lines are two possible correlations for accreting neutron stars (Migliari & Fender 2006). The solid green line is the maximum radio continuum luminosity detected for accreting white dwarfs (Körding et al. 2008(Körding et al. , 2011 tometric center could be offset from the dynamical center, but other evidence disfavors an intermediatemass BH interpretation. For example, the fundamental plane of accreting BHs (Merloni et al. 2003) implies that an intermediate-mass BH should have a much higher L R /L X than a stellar-mass BH.

THE HOST BINARY SYSTEM
Assuming that M62-VLA1 is a stellar-mass BH in M62, the intracluster gas density is very unlikely to be high enough to produce the observed radio and X-ray luminosity via Bondi accretion. Therefore, M62-VLA1 must be accreting from a binary companion. Assuming a nominal BH mass of 10 M ⊙ , M62-VLA1 would be accreting at ∼ (2 − 3) × 10 −7 L edd .
The relatively high quiescent X-ray luminosity of M62-VLA1, coupled with the period-luminosity correlation for BH binaries (Garcia et al. 2001;Reynolds & Miller 2011), suggests a longer orbital period than expected for a dwarf secondary. A giant donor is more likely. This would be yet another similarity between M62-VLA1 and V404 Cyg (which has a long orbital period for a BH binary at 6.5 days; Casares et al. 1992).
Consistent with this interpretation, the HST data discussed in §2.3 reveals a candidate optical counterpart consistent with a star on the lower red giant branch of M62. However, the F 435W − F 625W color of the star is unusually blue for an M62 giant, which could indicate the contribution of thermal emission from an accretion disk around the black hole. The photometry also suggests the presence of Hα emission, which could arise from an accretion disk.
In fact, the properties of the possible optical counterpart to M62-VLA1-including its unusual color-present a reasonable resemblance to V404 Cyg. Casares et al. (1993) estimate that for V404 Cyg, ∼26% and ∼9% of the B and R-band light, respectively, is attributable to accretion. Given these values, which are uncertain by perhaps a factor of two, the star is more than 0.2 mag bluer in B − R than it would be without the presence of an accretion disk (compare this to the ∼0.4 mag offset inferred for M62-VLA1). The equivalent width of Hα emission in V404 Cyg is ∼20Å, although somewhat time variable (Hynes et al. 2002(Hynes et al. , 2009. From our broadband and F 658N photometry of M62-VLA1, we estimate a rough Hα equivalent width of ∼7Å. Optical variability often accompanies accreting binaries (see Edmonds et al. 2003 for binaries in 47 Tucanae), and the HST variability observed in all three filters (including the F 658N filter covering Hα; see §2.3 and Fig-ure 5) supports the identification of the candidate giant counterpart with the radio source M62-VLA1.
For comparison, consider V404 Cyg, which exhibits variability on a wide range of time scales, from minutes to days (Casares et al. 1993;Pavlenko et al. 1996;Shahbaz et al. 2003;Zurita et al. 2004). In addition to its 6.5-day orbital period that produces ellipsoidal variations, it shows a ∼6-hr quasi-periodic oscillation with variable amplitude, 0.05-0.25 mag (Pavlenko et al. 1996); a similar mode might explain the variability plotted in Figure 5. In addition, when V404 Cyg varies the Hα and optical continuum usually scale together (Hynes et al. 2002); the behavior is also seen in the proposed optical counterpart to M62-VLA1 ( Figure 5).
The association between M62-VLA1 and this optical source could be bolstered by further high-resolution photometry, follow-up spectroscopy, and a proper-motion measurement. It should be kept in mind that M62-VLA1 could instead be associated with a star not detected in the HST images, such as an M dwarf-although this scenario might be disfavored, as the shorter orbital period would lead to a lower mass accretion rate and thus a lower quiescent X-ray luminosity than observed. The system could also, in principle, be a white dwarf-BH binary at an orbital period such that its quiescent luminosity would match the observations; the first extragalactic BH X-ray binary candidate in a GC is likely to be such a system Zepf et al. 2008).

ALTERNATIVE EXPLANATIONS
While the properties of M62-VLA1 are well-described by an accreting stellar-mass BH in M62, here we consider other possible scenarios which might produce this source. Additional discussion for some of these possibilities can be found in the Supplementary Information of Strader et al. (2012a).

Accreting Neutron Star or White Dwarf
As shown in Figure 6, neutron stars or white dwarfs accreting from binary companions are expected to have much lower L R /L X and L R than M62-VLA1 (Fuerst et al. 1986;Migliari & Fender 2006;Körding et al. 2008Körding et al. , 2011Byckling et al. 2010). An exception is the class of symbiotic stars, which are white dwarfs accreting from luminous red giants. Radio emission in symbiotic stars is usually thermal, originating from the ionized red giant wind. Therefore, high L R /L X values are observed in these systems (log L R /L X ≈ −2 to −1, defined for a radio frequency of 6.2 GHz and X-ray energy range 0.5-8 keV; Maccarone et al. 2012). This ratio is significantly higher than observed for M62-VLA1 (log L R /L X = −4.7). Symbiotic stars also show inverted radio spectra (α ≈ 0.6; Seaquist & Taylor 1990), inconsistent with M62-VLA1. Finally, the giant star which has been tentatively associated with M62-VLA1 is also less optically luminous than typical symbiotic stars. For these reasons, combined with the fact that there are no known symbiotic stars in GCs, we consider it unlikely that M62-VLA1 is a symbiotic star.
We can offer a contrived scenario that could make the observations of M62-VLA1 consistent with that for an accreting neutron star, if we assume the source is strongly variable. The radio and X-ray data on M62 were taken 10 years apart. The "high" radio luminosity could represent a strong flare on a neutron star, while the accompanying rise in L X may have remained below the detection threshold of X-ray all-sky monitors. The Chandra data would then need to represent a quiescent phase. We note that the duty cycle of neutron star flares is probably low but is poorly constrained.
Strong neutron star outbursts are usually accompanied by a rapid rise and decline in radio luminosity on short timescales (factor of 3 in ∼4 days; Moore et al. 2000;Tudose et al. 2009;Miller-Jones et al. 2010). To test this predicted evolution, we imaged each of the VLA observing blocks separately and measured the flux density of M62-VLA1 in six distinct epochs spanning a week (the second and third observing blocks were scheduled directly adjacent to one another and combined into one epoch); the light curve is plotted in Figure 7. M62-VLA1 varies by a factor of 3, although its variations are not in the form of a single flare. Instead, the flux density drops over the first two days to a non-detection during the third epoch, and then rises again to reach a maximum in the last epoch. The light curve of M62-VLA1 therefore appears to be inconsistent with a single strong outburst of a neutron star, and is more consistent with the "burbling" observed for quiescent stellar-mass black holes (Miller-Jones et al. 2008). Still, for a secure interpretation of M62-VLA1, simultaneous X-ray and radio data are highly desirable, and are being pursued.

Pulsar or Supernova Remnant
Millisecond pulsars are often found in the cores of GCs, but they display steep radio spectra (α −1; Kramer et al. 1999), inconsistent with the measured spectral index of M62-VLA1.
Pulsars in the field occasionally drive a wind that interacts with ambient material, termed a pulsar wind nebula. These objects have flat radio spectra, but again are more radio loud than M62-VLA1 (log L R /L X ≈ 0; Maccarone et al. 2012). In addition, they generally have L X > 10 34 erg s −1 , large diameters, short lifetimes, and dense surroundings, all inconsistent with observations of M62-VLA1 (Gaensler & Slane 2006). A similar analysis would apply to "redback" systems in which the material originates in an ablated companion (Roberts 2013).
Like pulsar wind nebulae, supernova remnants should be short lived, angularly extended, and located in dense environments. Supernova remnants in our Milky Way are more radio luminous than M62-VLA1 (>0.1 Jy) and have steeper radio spectra (α ≈ −0.7; Green 2009). For these reasons, we exclude pulsars, pulsar wind nebulae, and supernova remnants as likely identities of M62-VLA1.

Planetary Nebula
Planetary nebulae emit optically-thin thermal emission at radio wavelengths, and would therefore show a flat spectral index consistent with M62-VLA1. However, planetary nebulae have low X-ray luminosities (∼10 30 erg s −1 ; Montez et al. 2010), much fainter than M62-VLA1. In addition, many planetary nebulae are bright in Hα, but no nebulosity is observed around M62-VLA in archival HST /WFPC2 images obtained in the F656N filter.

Foreground Active Star
Foreground (D 100 pc) coronally-active late-type stars can be bright in both the X-ray and radio regimes, show flat radio spectra, and have log L R /L X ≈ −6 to −2 (Berger et al. 2010;Maccarone et al. 2012). The candidate optical counterpart of M62-VLA1 is a K giant, not a very late-type dwarf. On the other hand, this match could be spurious and an active dwarf star might escape undetected in the optical imaging if sufficiently distant ( 50 pc). To produce detectable radio and Xray emission at these distances, the star would need to be flaring (Berger et al. 2010). However, a flaring star should also be strongly circularly polarized (Berger 2002;Hallinan et al. 2008), and there is no evidence that M62-VLA1 is circularly polarized (3σ upper limit 22%).

Background Galaxy
It is possible, but very unlikely, that M62-VLA1 is a background radio galaxy. Galaxies exhibit a wide range of log L R /L X from 0 to −5 (Maccarone et al. 2012). In addition, they show a diversity of radio spectral slopes, although steeper indices than measured for M62-VLA1 are typical (α ≈ −0.7; Fomalont et al. 2002;Randall et al. 2012). Background galaxies with flat radio spectra are either active galactic nuclei with partially self-absorbed jets, or star-forming galaxies in which the combination of synchrotron and thermal emission produces a flat composite spectrum. Background objects with flux densities as faint as M62-VLA1 (∼20 µJy) have not yet been characterized in detail.
We use calculations similar to those in Strader et al. (2012a) to determine the expected number of background sources with properties comparable to those of M62-VLA1. Fomalont et al. (2002) estimate that there are 0.3 sources per sq. arcmin with 8.4 GHz flux densities in the range 10-30 µJy. We would therefore expect ∼0.01 sources in the 6.6 ′′ radius core of M62 with similar flux densities as M62-VLA1. This number decreases to ∼0.003 when the spectral index is considered; Fomalont et al. (2002) find that only 30% of the 47 faint radio sources (< 35 µJy) in their sample have α ≥ −0.4, measured from VLA 1.4 and 8.4 GHz data. The lack of a visible background galaxy further reduces the probability that M62-VLA1 is a background object: if we conservatively assume that we could detect optical counterparts to I ≈ 21.5, about a third of the faint flat-spectrum radio sources in the Fomalont et al. sample would have detectable optical matches. On the basis of the radio and optical data, we therefore estimate that the probability that M62-VLA1 is a background galaxy is < 0.003.
To date we have only investigated four GCs with deep VLA imaging (Strader et al. 2012b, plus M62); two of these GCs have candidate stellar-mass BHs. Therefore it is unlikely that the detection of a BH look-alike in M62 is due to a "multiple trials" effect. Nonetheless, a direct elimination of the background-galaxy possibility is still needed. As in the accreting neutron star scenario, simultaneous radio and X-ray observations will help distinguish between a stellar-mass BH and a background galaxy. The fundamental plane of BH activity suggests that L R /L X ∝ M 0.38 BH (Merloni et al. 2003). If simultaneous data show a similar ratio as currently observed, then M62-VLA1 is unlikely to be a supermassive BH in a background active galactic nucleus. In addition, deeper HST photometry would allow us to search for faint background galaxies coincident with M62-VLA1. Finally, deep highresolution radio imaging with the High Sensitivity Array would allow us to test if M62-VLA1 shares the proper motion of M62 itself (Dinescu et al. 2003).

CONCLUSIONS AND IMPLICATIONS
M62-VLA1 is the only radio source in the core of M62 down to flux density levels of ∼6 µJy at 5-10 GHz. We have argued that it is likely to be a stellar-mass BH accreting from a binary companion. This scenario is supported by its X-ray and radio luminosities and spectra. The secondary may be a giant star near the base of the red giant branch, as suggested by HST imaging and the L X -period relation for known BHs.
While all of our data are consistent with a stellar-mass BH interpretation, other explanations are possible (although unlikely). The most plausible alternative scenarios are (i) an accreting neutron star X-ray binary that is in quiescence during the Chandra observations but flaring during the VLA observations; or (ii) a background galaxy. Future observations-principally simultaneous X-ray and radio imaging-will conclusively distinguish between these scenarios.
M62 is now the second Milky Way GC with significant evidence for stellar-mass BHs (see Strader et al. 2012a for discussion of the BH candidates in M22). We have identified BH candidates in two out of four GCs with deep VLA radio imaging. While the sample is certainly small, these findings suggest that stellar-mass BHs may be present in many GCs. However, the current sample is not representative of the mass or core density distribution of the Galactic GC system, so our findings are not yet easily generalized. We note that both M62 and M22 are massive GCs, but their structural parameters are quite different. While M22 has a large core radius, M62 has a relatively dense, compact core. This contrast may imply that there are no clear structural diagnostics of the presence BHs in GCs (in contrast with theoretical findings by e.g., Mackey et al. 2008;Sippel & Hurley 2013;Morscher et al. 2013).
It is worth emphasizing that the BHs which can be detected in accreting binary systems are likely to represent only a small fraction of the population of BHs in a star cluster, since single BHs and those in binaries not undergoing mass transfer may dominate the BH population in a GC (e.g., Ivanova et al. 2010).
If indeed many GCs host substantial populations of BHs, this represents an important shift in our view of stellar remnants in GCs. The potential implications are large: GCs would become good hunting grounds for new stellar-mass BHs, of which few are still confirmed; BHs in GCs are likely to be more massive than those in the field, opening up a new area of BH parameter space; it would be possible to determine accurate physical parameters for a significant sample of BHs, since the distances to GCs are accurately known; additional quiescent BH systems would offer new opportunities to investigate the physics of accretion at the lowest luminosities; and if GCs have multiple BHs, the likelihood of the dynamical formation of BH-BH and BH-pulsar systems will increase, improving prospects for the detection of gravitational waves and stringent tests of general relativity. Deep radio imaging of additional GCs, combined with detailed follow-up observations of the candidate BHs in M62 and M22, will determine whether GCs-long thought to be among the least likely hosts for BHs-turn out to be among the best.
We thank S. Sigurdsson and N. Ivanova for useful comments. L. Chomiuk is a Jansky Fellow of the National Radio Astronomy Observatory. J. C. A. Miller-Jones acknowledges support from an Australian Research Council Discovery Grant (DP120102393). C. Heinke is supported by NSERC and an Ingenuity New Faculty Award. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The search for an optical counterpart was based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). The scientific results reported in this article are based to a significant degree on data obtained from the Chandra Data Archive.