A Radio Counterpart to a Jupiter-mass Binary Object in Orion

Using James Webb Space Telescope near-infrared data of the inner Orion Nebula, Pearson & McCaughrean detected 40 Jupiter-mass binary objects (JuMBOs). These systems are not associated with stars and their components have masses of giant Jupiter-like planets and separations in the plane of the sky of order ∼100 au. The existence of these wide free-floating planetary-mass binaries was unexpected in our current theories of star and planet formation. Here we report the radio continuum (6.1 and 10.0 GHz) Karl G. Jansky Very Large Array detection of a counterpart to JuMBO 24. The radio emission appears to be steady at a level of ∼50 μJy over timescales of days and years. We set an upper limit of ≃15 km s−1 to the velocity of the radio source in the plane of the sky. As in the near-infrared, the radio emission seems to be coming from both components of the binary.


INTRODUCTION
In a new James Webb Space Telescope (JWST) nearinfrared survey of the inner Orion Nebula and Trapezium Cluster, Pearson & McCaughrean (2023) detected 40 Jupiter-Mass Binary Objects (JuMBOS).Additionally, two triple objects were reported.These multiple systems are not associated to stars and their components have masses in the range of 0.6 to 14 Jupiter masses and separations in the plane of the sky between 28 and 384 au (0.07 and 0.99 arcsec at a distance of 388 pc; Kounkel et al. 2017).The existence of these wide planetary-mass binaries is surprising since our current knowledge of star and planet formation cannot account for them.A possible theoretical explanation has been offered by Wang et al. (2023) in terms of an ejection after a close stellar flyby between two solar-like stars.Rogue planets (single planets not gravitationally bound to any star or brown dwarf, also known as free-floating or orphan planets) have been detected in the optical and near-infrared (Lucas & Roche 2000;Zapatero Osorio et al. 2000).This latter class of single sources, however, is compatible with the known physical conditions of forming stellar and planetary systems (Hurley & Shara 2002;Kroupa & Bouvier 2003;Whitworth & Zinnecker 2004;Miret-Roig 2023).We have searched for radio counterparts to the 42 JuMBOS in Orion using sensitive observations from the archive of the Karl G. Jansky Very Large Array (VLA) of the National Radio Astronomy Observatory.
Here we report the radio detection of JuMBO 24 in three epochs distributed across a decade, in 2012, 2018, and 2022.This paper is organized as follows.In Section 2 we present the observations, while in Section 3 we discuss the implications of this detection.Finally, in Section 4 we summarize our conclusions.

VLA Observations
We analyzed the VLA data of the three projects listed in Table 1.The data were calibrated in the standard manner using the CASA (Common Astronomy Software Applications; McMullin et al. 2007) package of NRAO and the pipeline provided for VLA 1 observations.A robust weighting of 0 (Briggs 1995) was used in the imaging to optimize the compromise between sensitivity and angular resolution.The images were made with a (u, v) range set to > 100 kλ, to suppress structures larger than 2 ′′ that limit the signal-to-noise ratio in the Orion field.
The project SD0630 consists of five sessions made on 2012 September 30 and October 2, 3, 4, and 5.The projects 17A-069 and 22A-285 were made on single sessions on 2018 May 11 and 2022 July 2, respectively.In Table 1 we also give the mean epoch of the observations, the configurations used, the central frequency and the bandwidth of the observations, the parameters of the synthesized beam, the radio flux density and position of JuMBO 24, and the 4-σ upper limit to the absolute circular polarization.The gain calibrator, J0541−0541, was the same for all observations so the astrometry can be directly compared more easily.The absolute flux calibrator was J0137+331 (3C 48) for projects SD0630 and 22A-285, and J0521+1638 (3C 138) for project 17A-069.
We find that the near-infrared source JuMBO 24 coincides within 0.05 arcsec with a radio continuum source present in all three VLA images.The position of JuMBO 24 given by Pearson & McCaughrean (2023) is RA(J2000) = 05 h 35 m 19.s 50616, DEC(J2000) = −05 • 23 ′ 39. ′′ 7303.The positions of the radio counterpart source for the three epochs are given in Table 1.This radio source first appeared in the list of Forbrich et al. (2016, their source 441).JuMBO 24 was originally detected in near-infrared Hubble Space Telescope observations in the F160W (1.4-1.8 µm) and F110W (0.8-1.4 µm) filters, with magnitudes of 17.11 and 18.33, respectively (Luhman et al. 2000, their source 274).Slesnick et al. (2004) reported a K-band magnitude of 19.51 (their source 728).These authors used near-infrared photometry and spectroscopy as well as evolutionary stellar models (D'Antona & Mazzitelli 1994;Baraffe et al. 1998) to determine an infrared spectral type of M5.5±1.5 and a mass of ∼ 52 M J , under the assumption of a single star.The JWST data and the interpretation of Pearson & McCaughrean (2023) indicate that JuMBO 24 is a binary with both components having the same giant-planet mass of 11.5 M J and a separation of 28 au in the plane of the sky.JuMBO 24 stands out in the list of these objects because it has the largest total mass (23 M J ) and the closest separation in the plane of the sky (28 au).Also, the primary of JuMBO 24 has the second smallest extinction, A V = 3.6, behind only JuMBO 27, that has A V = 2.4.

DISCUSSION
The first centimeter detection of a single planetary mass object, SIMP J01365662+0933473, was achieved by Kao et al. (2018).The first identification of a centimeter continuum source with a planetary-mass binary object presented here is also of great relevance by the information the radio data can provide.We exemplify these possibilities in what follows.Further work, particularly in relation to the origin of the radio emission is deferred to forthcoming papers.

The Probability of a Random Association
In this section we estimate how likely it is that the positional agreement of JuMBO 24 with one of the radio sources in the region is fortuitous and not real.In their study of Orion, Forbrich et al. (2016) detected a total of 556 radio sources in a region of 13. ′ 6 × 13. ′ 6.Most of the sources, however, are enclosed in a region with radius of ∼ 3 ′ (see Figure 5 of Forbrich et al. (2016)).Of the 556 sources, 358 have flux densities equal or larger than 50 µJy, the characteristic value for the radio source proposed to be associated with JuMBO 24.These results give a solid angle density of 3.5 × 10 −3 radio sources per arcsec 2 .
The most accurate position for the radio source is that obtained from the multiepoch project SD0630.Convolving the statistical errors of the position of the radio source (Table 1) with the systematic positional error of 0. ′′ 01, as estimated by Dzib et al. (2017), we obtain a ±1 − σ error box of 0. ′′ 027 × 0. ′′ 031.The a priori probability of a source randomly falling in this error box is then 2.9 × 10 −6 .However, the coincidence of the radio source with any of the 42 JuMBOs in Orion would have been considered significant, so assuming that these possible coincidences are independent events we have to multiply the probability of a given JuMBO source being spatially associated times 42.We finally get an a priori probability of one JuMBO randomly coinciding in position with one radio source to be only 1.2 × 10 −4 , supporting a real association.Note-Columns 4 and 5 are the central frequency and bandwidth of the observations.The position given in columns 8 and 9 show only the seconds of the RA and the arcseconds of the DEC, measured from RA(J2000) = 05 h 35 m and DEC(J2000) = −05 • 23 ′ , respectively.Column 10 gives the 4-σ upper limit to the absolute circular polarization.

Upper Limits to the Proper Motion
Using the radio positions determined at the three epochs we can search for the proper motions of the source.As noted before, the fact that the same nearby (∼ 1. • 5) gain calibrator (J0541−0541) was used in all observations implies a sound astrometry.A least-squares fit to the positions given in Table 1, after adding in quadrature a systematic error of 0. ′′ 01, yields the following values: This result implies that there are no significant proper motions at a 3-σ level of ∼8 mas yr −1 , that corresponds to a transverse velocity of ∼ 15 km s −1 at a distance of 388 pc (Kounkel et al. 2017).These upper limits rule out large velocities in the plane of the sky, although it should be pointed out that significant velocities could be present in the radial direction.Our measurements can be registered to an approximate rest frame for the Orion Cluster by subtracting the average proper motion of the radio stars in the core of Orion (Dzib et al. 2017), µ α cos δ = +1.07± 0.09 mas yr −1 ; µ δ = −0.84± 0.16 mas yr −1 .We obtain: We conclude that JuMBO 24 is not moving at large velocities compared to the radio stars in Orion.van Elteren et al. (2019) have made numerical simulations for the evolution of young and dense stellar clusters.They find that single free-floating planets ejected from their original system after a strong encounter with another star will have velocities typically three times those of the stars in the cluster.Dzib et al. (2017) find a typical velocity dispersion of ∼ 2 − 3 km s −1 along a plane-of-thesky coordinate for the radio-emitting stars in the core of the Orion Nebula Cluster.Assuming that the results of van Elteren et al. ( 2019) can be extended to free-floating binary planets, we expect velocity dispersions of ∼ 6 − 9 km s −1 for the JuMBOs as well as for the ∼500 single planetary-mass objects reported in Orion by Pearson & McCaughrean (2023).Future highly accurate determinations of the kinematics of the Orion JuMBOs will help understand the origin of these unexpected sources.

Time Monitoring
The two measurements at 6.1 GHz given in Table 1 are separated by a decade and, within the noise, they are consistent with being equal.Another evidence that weighs in favor of a steady flux is the analysis of the five different epochs observed in project SD0630.In Table 2, we give the epochs and flux densities for each individual observing session, where we can see that all flux densities are consistent within noise with a value of ≃ 45 µJy.While more intensive observations are needed to establish firmly the time behavior of the radio emission from JuMBO 24, all existing data are consistent with a steady flux density over timescales of days to years.

Nature of the Radio Emission
The objects with known radio emission that are more similar to the components of the JuMBOs are the ultracool dwarfs (Kao & Sebastian Pineda 2022).The radio emission from the latter objects has two main components: a periodically bursting emission arising in aurorae and a slowly varying, quiescent emission with a low degree of circular polarization arising from a structure similar in morphology to the Jovian radiation belts (Climent et al. 2023;Kao et al. 2023).The bursting component is circularly polarized and it is believed to be produced by the electron cyclotron maser instability (Hallinan et al. 2006), the same coherent emission process that creates the auroral radio emission from the planets of the Solar System (Pineda et al. 2017).On the other hand, the quiescent component is usually interpreted as being caused by radio emission from relativistic electrons giving rise to gyrosynchrotron or synchrotron mechanisms, according to how relativistic the electrons are.Now it is known that this quiescent emission is produced in radiation belts (Climent et al. 2023;Kao et al. 2023).Given that the radio emission from JuMBO 24 is steady in time, one would expect it to be produced by relativistic electrons spiraling in the magnetic fields of a radiation belt.
The spectral index of the radio emission is an indicator of the nature of the radiation.The observations of SD0630 were made using two bandwidths of 1 GHz each, centered at 4.8 and 7.4 GHz, respectively.We obtain that the flux densities at these two frequencies are 53.8±6.0 and 43.0±5.8µJy, respectively.Assuming that the continuum spectrum is described by a power law, S ν ∝ ν α , we derive a spectral index of α = −0.52 ± 0.40.We can also derive a spectral index assuming that the source has a steady flux density and combining the values given in Table 1.From these data we derive α = +0.50 ± 0.43.The weighted mean and rms of these two determinations is α = −0.05±0.51.Given the large uncertainty, it is difficult to discriminate between possible emission mechanisms, but we can rule out an optically-thick blackbody in the Rayleigh-Jeans regime (α = 2) and the negative steep indices characteristic of pulsars (α ≤ -2).
We did not detect circular polarization at any epoch.However, the upper limits are not very stringent.The presence of circular polarization is a strong indicator of either electron cyclotron maser instability, cyclotron or gyrosynchrotron radiation but, of course, its absence does not rule out these mechanisms.
As noted above, the ultracool dwarfs are the closest thing to giant planets that have been detected in the radio.The radio luminosities of nine quiescent binary ultracool dwarf systems discussed by Kao & Sebastian Pineda (2022) are in the range of 10 11.9−14.6erg s −1 Hz −1 .JuMBO 24 has a radio luminosity of 10 16.0 erg s −1 Hz −1 .We then conclude that the radio luminosity of JuMBO 24 is orders of magnitude larger than naively expected.

Position and Morphology of the Radio and Infrared Sources
In Figure 1 we show an image of the radio continuum emission obtained from concatenating all the data of project SD0630.We can see that the radio source coincides closely with the near-infrared position given by Pearson & McCaughrean (2023).A Gaussian ellipsoid fit to the radio image using the task IMFIT of CASA gives a deconvolved angular size of 0. ′′ 179 ± 0. ′′ 045 × 0. ′′ 143 ± 0. ′′ 074; 60 • ± 89 • .Although the radio source is only marginally resolved, these angular dimensions and the orientation of the major axis are consistent with the value obtained from Gaussian ellipsoid fits to the JWST images.For instance a fit to the image at 2.77 µm results in 0. ′′ 170 ± 0. ′′ 002 × 0. ′′ 132 ± 0. ′′ 002; 97 • ± 2 • .The shorter wavelength JWST observations show that the cause of the extension along the east-west direction is the binarity of the source, which is comprised of two objects separated by about 100 mas precisely along the east-west direction.This result suggests that, as in the near-infrared, there is significant radio emission from both components of the binary.This result is similar to that obtained for the compact ultra-cool dwarf binary VHS 1256−1257AB, where the optical and radio emissions show compatible dimensions and orientation (Rodriguez et al. 2023).

CONCLUSIONS
We report the detection of radio continuum emission at 6 to 10 GHz associated with the recently identified Jupiter-Mass Binary Object JuMBO 24 in the Orion Nebula Cluster.Within the errors, the radio emission appears steady on timescales from days to years and exhibits no circular polarization.No significant proper motion is detected between the three observing epochs, ruling out large velocities (≥ 15 km s −1 ) in the plane of the sky.The radio emission is marginally resolved in the same direction as the infrared source detected by the JWST, suggesting that the radio emission comes from a combination of the two planetary mass objects.Additional radio observations are necessary to pin down the nature of the radio emission mechanism.

Figure 1 .
Figure 1.VLA image of the 6.1 GHz emission from JuMBO 24, made concatenating all the data of project SD0630.The contours start at, and increase by steps of, 9 µJy beam −1 .The color scale indicates intensity as shown by the color bar to the right.The white rectangle indicates the position of JuMBO 24 as reported by Pearson & McCaughrean (2023): RA(J2000) = 05 h 35 m 19.s 50616, DEC(J2000) = −05 • 23 ′ 39. ′′ 7303.The red ellipse shows the angular size and orientation of the target in the 2.77 µm JWST images (0. ′′ 170 × 0. ′′ 132; 97 • ).The synthesized beam of the radio map is shown at the bottom left of the image.

Table 2 .
Flux Densities of JuMBO 24 at 6.1