{'fraction_non_alphanumeric': 0.038120416169196654, 'fraction_numerical': 0.03297515492637444, 'mean_word_length': 4.2750037486879595, 'pattern_counts': {'":': 0, '<': 5, '<?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}}
Predicting radio emission from the newborn hot Jupiter V830 Tau and its host star

Magnetised exoplanets are expected to emit at radio frequencies analogously to the radio auroral emission of Earth and Jupiter. We predict the radio emission from V830 Tau b, the youngest (2 Myr) detected exoplanet to date. We model the host star wind using 3DMHD simulations that take into account its surface magnetism. With this, we constrain the local conditions around V830 Tau b that we use to then compute its radio emission. We estimate average radio flux densities of 6 to 24mJy, depending on the assumed radius of the planet (one or two Rjupiter). These radio fluxes are present peaks that are up to twice the average values. We show here that these fluxes are weakly dependent (a factor of 1.8) on the assumed polar planetary magnetic field (10 to 100G), opposed to the maximum frequency of the emission, which ranges from 18 to 240MHz. We also estimate the thermal radio emission from the stellar wind. By comparing our results with VLA and VLBA observations of the system, we constrain the stellar mass-loss rate to be<3e-9 Msun/yr, with likely values between ~1e-12 and 1e-10 Msun/yr. The frequency-dependent extension of the radio-emitting wind is around ~ 3 to 30 Rstar for frequencies in the range of 275 to 50MHz, implying that V830 Tau b, at an orbital distance of 6.1 Rstar, could be embedded in the regions of the host star's wind that are optically thick to radio wavelengths, but not deeply so. Planetary emission can only propagate in the stellar wind plasma if the frequency of the cyclotron emission exceeds the stellar wind plasma frequency. For that, we find that for planetary radio emission to propagate through the host star wind, planetary magnetic field strengths larger than ~1.3 to 13 G are required. The V830 Tau system is a very interesting system for conducting radio observations from both the perspective of radio emission from the planet as well as from the host star's wind.


Introduction
Detecting radio emission from exoplanets would revolutionise the area of exoplanetary studies. First, it would open up a new avenue for the direct detection of exoplanets (Farrell et al. 1999). Second, it would allow us to measure exoplanetary magnetic fields, which so far have only been elusively probed (Shkolnik et al. 2008;Vidotto et al. 2010a;Kislyakova et al. 2014). Exoplanetary magnetic fields can reveal information about a planet's interior structure and dynamics (e.g. whether it has a convecting, electrically conducting interior) and are also believed to play an important role in conditions for habitability, by protecting the planet's surface from energetic cosmic particles, protecting its atmosphere from violent chemical changes, and potentially helping atmospheric retention (Lundin et al. 2007;Grießmeier et al. 2016).
The magnetised planets in the solar system generate radio emission through the electron cyclotron maser instability. The power of this emission is related to the kinetic and/or magnetic powers of the incident solar wind over several orders of magnitude (Farrell et al. 1999;Zarka et al. 2001). Although the physics of this relation are still not well understood, this "radiometric Bode's relation" indicates that planetary radio emission is somehow powered by the interaction between the planetary magnetic field and the solar wind. In analogy to the solar system, it has been suggested that exoplanets may also produce radio emission due to their interaction with the winds of their host stars (Farrell et al. 1999). 1 Calculations based on the solar-system analogy suggest that radio emission from close-in exoplanets should be several orders of magnitude larger than the largest radio emitter of the solar system, Jupiter (e.g., Zarka et al. 2001;Lazio et al. 2004;Grießmeier et al. 2005;Zarka 2007;Vidotto et al. 2010b;Fares et al. 2010;Hess & Zarka 2011). This is because, at the short orbital distances of close-in planets, the local power of the host star wind is significantly higher than, for example, the local power of the solar wind dissipated at the distances of the planets in our solar system. These and many other theoretical predictions have motivated a good number of observational searches of exoplanetary radio emission, yielding mostly negative results (e.g. Bas-tian et al. 2000;Ryabov et al. 2004;Lazio & Farrell 2007;Hallinan et al. 2013, but see also Lecavelier des Etangs et al. 2013;Sirothia et al. 2014, who found ambiguous hints of the existence of exoplanetary radio emissions). There are several reasons for the non-detections, such as a lack of exoplanetary magnetism, a mismatch of the frequency search, and/or a low sensitivity of the present-day instruments (Bastian et al. 2000, Grießmeier et al. 2011.
One might overcome the issue of instrumental sensitivity by turning to exoplanets that, not only orbit closer to their stars, but that orbit stars that have winds that are more powerful than solar-type winds, and whose magnetic fields are also stronger. Stars that are more active than the Sun are indeed expected to host more powerful winds (Wood 2004;See et al. 2014;Johnstone et al. 2015;do Nascimento et al. 2016). In particular, stars in the pre-main sequence with close-in planets would be ideal targets to search for planetary radio emission (Grießmeier et al. 2005;Vidotto et al. 2010b). With the recent detection of the first hot giants orbiting stars younger than 10-20 Myr (V830 Tau b, Donati et al. 2016, and K2-33b, David et al. 2016, TAP26b, Yu et al. 2017, now is an excellent timing for testing theoretical expectations of exoplanetary radio emission. In this paper, we compute the expected radio emission from the newborn hot Jupiter V830 Tau b. With the spectropolarimetric monitoring of V830 Tau, Donati et al. (2016) were able to extract the planetary signature of V830 Tau b, a 2 Myr-old hot Jupiter orbiting at 0.057 au, while, at the same time, reconstructing the large-scale magnetic field topology of the host star. It is precisely the latter measurements that we use in this paper to derive the wind properties of the host star (Section 2). With that, we are then able to derive the expected radio emission of V830 Tau b, starting from the assumption that V830 Tau b is magnetised (Section 3). Given the high densities of the wind of young active stars, these winds can become optically thick to free-free radiation at radio wavelengths. We thus investigate in which circumstances the planet would be embedded in the radio-emitting region of the stellar wind and whether planetary radio emission would be able to propagate through the stellar wind plasma, given that the planetary cyclotron frequency of emission should exceed the plasma frequency of the stellar wind (Section 4). For physicallyreasonable conditions, we show that planetary radio emission from V830 Tau b is expected to peak at 12 mJy, presenting therefore a high potential of detection with LOFAR (the Low-Frequency Array) the upgraded UTR-2 (Ukrainian T-shaped Radio telescope) and GMRT (Giant Metrewave Radio Telescope ,see Section 5). Interestingly, its host star has been detected using VLA (the Karl G. Jansky Very Large Array) at 4.5 and 7.5 GHz and with VLBA (the Very Long Baseline Array) at 8.4 GHz, becoming the first exoplanet host with detected radio emission (Bower et al. 2016).

Stellar wind modelling
We model the wind of V830 Tau by means of three-dimensional (3D) magnetohydrodynamics (MHD) simulations, similarly to Vidotto et al. (2012Vidotto et al. ( , 2015. We refer the reader to these papers for more details of the model, which are summarised next. We use the 3D MHD numerical code BATS-R-US (Powell et al. 1999;Tóth et al. 2012), modified as in Vidotto et al. (2012). BATS-R-US solves the set of ideal MHD equations for the mass density, velocity, magnetic field, and gas pressure. We assume the wind is polytropic with a polytropic index γ = 1.15 and a fully ionised hydrogen wind.
For the physical characteristics of the star, we use a rotation period of 2.741 days, mass of 1M , and radius of R = 2R . The radial part of the stellar magnetic field, anchored at the wind base, is constrained from observations collected in November-December 2015 (Donati et al. 2016). Initially, the field is considered to be in its lowest (potential) state, but as the simulation evolves in time the wind particles self-consistently interact with the magnetic field lines (and vice-versa), removing the field from its initial potential state. At the wind base, we adopt a temperature of 10 6 K and a number density of n 0 = 10 12 cm −3 , which are free parameters of our model. The star rotates as a solid body with a rotation axis along the z-axis. Our grid is the same as described in Vidotto et al. (2014b). The resultant wind solution, obtained self-consistently, is found when the system reaches steady state in the reference frame corotating with the star. For the parameters we adopted, we obtain a wind mass-loss rate ofṀ ∼ 3 × 10 −9 M yr −1 . This mass-loss rate is an upper limit for the case V830 Tau, as we will demonstrate in Section 4. Figure 1 illustrates the output of the wind simulations of V830 Tau. We overplot the (assumed equatorial and circular) orbital radius of V830 Tau b (black circle) and a cut of the Alfvén surface S A at the equatorial plane (xy plane). For the parameters we adopted, the planet's orbit lies almost completely inside S A . The Alfvén surface is expected to expand for lower values of n 0 (Vidotto et al. 2009). This indicates that V830 Tau b is truly orbiting in the sub-Alfvénic regime, if not at all times, at least during most of its orbit. The Alfvén surface represents the boundary between a magnetically-dominated stellar wind (interior to S A ) and a region that is dominated by the wind inertia. The Alfvén surface is relevant for computing angular momentum losses (Weber & Davis 1967) and indicates the types of star-planet interactions. Planets orbiting in the sub-Alfvénic regime can be directly connected to the magnetic field of the star and perturbations can travel to/from the star (e.g. Preusse et al. 2005;Saur et al. 2013;Strugarek et al. 2015). In this case, the planetary magnetosphere is embedded in that of the star and it is believed that there is continuous reconnection of the planet's and stellar magnetic fields, similar to the Ganymede-Jupiter system. When the planet orbits outside S A (in the super-Alfvénic regime), as in the case of the Earth, the magnetospheric cavity surrounding the planet deflects the stellar wind particles. We estimate the size of the magnetosphere of V830 Tau b in the next section.

Exoplanetary radio emission
To calculate the radio emission of V830 Tau b, we use the radiometric Bode's law for the solar system, in which the planetary radio power P radio can be decomposed into a power released from the dissipation of kinetic energy of the stellar wind P k and/or a power released from the dissipation of magnetic energy of the wind P B (e.g. Farrell et al. 1999;Zarka et al. 2001) where η k and η B are efficiency ratios. In the solar system, η k = 1 × 10 −5 and η B = 2 × 10 −3 (Zarka 2007); we assume the same efficiency ratios here. In the past, it has been argued that it was not possible to decide which incident power actually drives the radio power observed from the magnetic planets of the solar system (Zarka 2007). Recently, however, Zarka (2010) argued that the magnetic field is likely to be a determinant for extracting part of the flow power and converting it to energetic particles (see e.g. Jardine & Cameron 2008;Nichols & Milan 2016, for a theoretical basis of this empirical relation). Here, for completeness, we compute the radio powers coming from both kinetic and magnetic energies. As we will see, the latter is much larger than the former and, following the arguments in Zarka (2010), is the preferred method for estimating radio fluxes. The dissipated kinetic and magnetic powers of the impacting wind on the planet are approximated as, respectively, (Appendix A) where B ⊥ is the magnetic field component perpendicular to ∆u.
Here, |∆u| = |u − v K | is the relative velocity between the wind and the Keplerian velocity (v K ) of the planet, and ρ, p , and B are the local density, pressure, and magnetic field intensity of the stellar wind. Neglecting planet thermal pressure, the size of the planet's magnetopause r M can be estimated by pressure balance between the stellar wind total pressure and the planet's magnetic pressure: p tot = B 2 p (r M )/(8π), where p tot = ρ∆u 2 + B 2 8π + p is the sum of the ram, magnetic, and thermal pressures of the wind; B p (r M ) is the intensity of the planet's magnetic field at the nose of the magnetopause. For a dipolar planetary magnetic field: where R p is the planetary radius, R is the radial coordinate centred at the planet, and B p is the polar magnetic field intensity (at the equator, the intensity is 1 2 B p ). Assuming the dipole is aligned with the planetary orbital spin axis, the magnetospheric radius (where R = r M ) is (4) Figure 2a shows the stand-off distance of V830 Tau b's magnetopause calculated using Eq. (4) as the planet orbits around the star. 2 We assume three different values of B p : 10, 50, and 100 G. For comparison, we note that the maximum intensity of Jupiter's magnetic field is 14.3 G (Bagenal 1992). For B p = 10 G, the average magnetospheric radius is ∼ 1.31R p (i.e. very close to the planet surface), while in the second and third cases, r M ∼ 2.2R p and 2.8R p , respectively. If the wind densities (and therefore mass-loss rates) were to decrease by a factor of 100, the magnetospheric sizes would increase by about 20% of the values presented in Figure 2a (in such a case, the planet's orbit would lie completely within the Alfvén surface). This shows that, even in the case of a "high" planetary magnetic field and low wind densities, the magnetosphere of V830 Tau b is expected to be quite small, due to the harsh conditions of the external environment. The radio flux is related to the radio power as where ω is the solid angle of the (hollow) cone of emission, d is the distance to the system (147 pc in the case of V830 Tau), and ∆ f is the emission bandwidth, approximately the cyclotron frequency (Grießmeier et al. 2007). In Appendix A, we show that the radio flux density due to the dissipated wind kinetic and magnetic powers simplify to respectively. Here η k and η B are constants and are related to η k and η B following Eqs. (A.13) and (A.19). The planetary magnetic field dependence is hidden in the function f (α 0 ), which varies between 0 and 3.3 for any planetary magnetic field intensity ( Figure A.1). In particular, we assume that the beaming angle of the radio emission occurs at a ring with colatitude α 0 and thickness δα = 17.5 o Equation (8) shows that the flux released in the dissipation of Poynting flux becomes dominant when the kinetic energy (ρ(∆u) 2 /2) of the impacting wind is smaller than 3200π 2 times the magnetic energy (B 2 ⊥ /8π). This is the case of the V830 Tau system, due to the intense stellar magnetism.
From Eqs. (6) and (7), we see that the radio fluxes are functions of the angular size of the planet (R p /d) and the properties of the stellar wind surrounding the planet. We also note that radio fluxes depend relatively weakly on the magnetic intensity of the planetary dipolar field. This dependence is hidden in f (α 0 ) (cf Appendix A). For example, for polar magnetic fields of 10, 50, and 100 G, respectively, the (northern) polar-cap boundaries are located at colatitudes α 0 of 61 o , 42 o , and 37 o , respectively. For these values of α 0 , f (α 0 ) is 2.2, 1.4, and 1.2, respectively. Therefore, as radio fluxes are proportional to f (α 0 ), a change of field intensities from 10 to 100 G results in a rather small change (a factor of 1.8) in the computed radio fluxes. This is good news, given that exoplanetary magnetism remains observationally elusive (cf Sect. 1). Figure 3 shows the predicted radio flux densities along a planetary year assuming dissipation of kinetic (top) and magnetic (bottom) stellar wind powers. We found that radio emission arising from the dissipation of wind magnetic power is considerably larger (∼ 30 times) than from kinetic power and this difference could become even larger should a lower stellar massloss rate be adopted. This is due to the kinetic-to-magnetic ratio being smaller than 3200π 2 (cf Eq. 8). In Appendix B, we show that our radio predictions for V830 Tau b are robust against a reduction of stellar wind densities and mass-loss rates, given the high magnetism of V830 Tau. We assume here that the planet has a radius of 1 to 2 R jup , but we note that, although the mass of V830 Tau b has been constrained by observations to be about 0.7 M jup (Donati et al. 2016), no constraints on its radius currently exist. Since radio fluxes are proportional to R 2 p , the fluxes calculated for R p = 2R jup are four times larger than for the calculations with R p = R jup . For the case of R p = R jup , the estimated magnetic radio fluxes are on average ∼ 6 mJy, with peaks at 11 mJy. This can be considered as a conservative estimate. Due to its youth and the close distance to the host star, it is likely that the radius of V830 Tau b is larger than 1R jup , since V830 Tau b is likely to be still in the contracting phase (e.g. K2-29 b, which has a similar mass, but orbits an older, 450 Myr-old star, has a radius of 1.19 R jup , Santerne et al. 2016). For the case of R p = 2R jup , the estimated magnetic radio fluxes are on average ∼ 24 mJy, with peaks at 44 mJy. We adopted in Figure 3 a dipolar magnetic field whose polar intensity is 10 G, but again note that an increase in magnetic field intensity by a factor of 10 would result in a change on the predicted radio flux by a factor of 1.8. The planetary magnetic field has, however, a strong influence on the frequency of emission. Since this is cyclotron emission, by measuring its frequency one is able to derive the intensity of the planetary magnetic field (cf Eq. A.3, see Hess & Zarka 2011 for a suggestion on how to conduct this estimation). Figure 2b shows the predicted maximum emission frequency, considering the same polar planetary magnetic fields as in Figure 2a. This frequency is calculated at the border of the polar-cap boundary, which is located at different colatitudes depending on the size of the planet magnetosphere (cf Eq. A.1, Vidotto et al. 2011). The frequencies are 18, 114, and 240 MHz for polar magnetic fields of 10, 50, and 100 G, respectively. The polar-cap boundaries extend down to latitudes of ±29 o , ±48 o , and ±53 o , respectively. The magnetic field intensity at these latitudes is shown in the right axis of Figure 2b. For comparison, at the Earth, the size of the polar cap is about 17 to 20 o (Milan et al. 2009), that is, extending to latitudes of ±70 to ±73 o , showing once more the extreme case of the magnetosphere of V830 Tau b.

Stellar versus exoplanetary radio emissions
If V830 Tau b is indeed a magnetised planet, it can emit at radio wavelengths due to a physical process called electron cyclotron maser instability (Section 3). However, V830 Tau b may not be the only radio emitter in this exoplanetary system. The hot plasma from the stellar wind can also emit at radio wavelengths, although through a different physical process. Here, we estimate radio emission from the host star's wind.

Radio emission from the stellar wind
A thermal ionised plasma emits bremsstrahlung (free-free) radiation across the electromagnetic spectrum. In the case of stellar winds, free-free radio emission is more intense in the innermost regions of the stellar wind, where the densities are higher (Panagia & Felli 1975;Wright & Barlow 1975;Lim & White 1996;Guedel 2002). For large enough densities, these regions can become optically thick to radio wavelengths and, if a planet is embedded in this region, then it is possible that most of the planetary radio emission gets absorbed and does not escape.
A simple way to calculate the radio emission of stellar winds was presented in Panagia & Felli (1975), who assume that the wind is spherically symmetric and isothermal with temperature T . These hypotheses are, however, not true in our simulations: the asymmetric distribution of surface magnetic fields generates an asymmetric stellar wind (Vidotto et al. 2014b). We also adopt a polytropic wind model, such that the wind temperature is not constant. To take into account the asymmetries of the stellar wind, a radiative transfer calculation performed on each cell of the numerical grid and integrated along the line-of-sight would provide more accurate predictions for the radio-emitting wind. This detailed study will be done in a future work. In the present paper, we compute the free-free emission of the inner regions of the wind of V830 Tau using the model developed by Panagia & Felli (1975) coupled to our simulation results. We caution that these computations should be considered as estimates of the radio emission.
The stellar wind radio flux density at a frequency ν is given by (Panagia & Felli 1975) S ν = 10 −29 A(α)R 2 5.624 × 10 −28 I(α)n 2 0 R where the functions I(α) and A(α) are given by and τ c = 3. The wind density is assumed to decay as a power-law with exponent α with n 0 being the wind base density. The distance R ν within which half of the emission S ν is produced is For the wind temperature T in Equations (9) and (13), we use the temperature adopted at the wind base (10 6 K) and we use the base density n 0 = 10 12 cm −3 adopted in the model described in Section 2. To estimate α, we compute a power-law fit for regions within r < 20R , which resulted in α 3.08. 3 With this value of α, the temperature and frequency dependencies of R ν become which is flatter than an emitting wind that has reached asymptotic speed (R ν ∝ ν −0.7 T −0.45 , Guedel 2002). Figure 4a shows the size of the radio-emitting region, where the dotted line indicates the orbital radius of V830 Tau b. We also explore the effects of the wind mass-loss rates on the radio emission of the wind. For this estimate, we scale the wind densities of the model presented in Section 2, such that the mass-loss rates of our parametric study vary between 10 −12 and 3 × 10 −9 M yr −1 . We remind readers that the model presented in Section 2 (our fiducial model) hasṀ ∼ 3 × 10 −9 M yr −1 . We note that R ν gets smaller for higher frequencies and/or lower mass-loss rates. At 10 MHz and forṀ in the range [3 × 10 −12 , 3 × 10 −9 ] M yr −1 , the planet's orbit is always within R ν , that is, the planet is embedded in the radio-emitting wind. For 100 MHz and in the same range ofṀ, R ν ∼ 6 to 90 R and r orb R ν forṀ 3 × 10 −12 M yr −1 . For 1000 MHz, r orb > R ν forṀ 3 × 10 −11 M yr −1 .
Regarding the flux density of the emitting wind, S ν increases for higher frequencies and/or higher mass-loss rates. For all the range ofṀ we investigated and the range of frequencies ([1, 6000] MHz), the flux density of the emitting wind never reaches the mJy level and S ν 0.1 mJy. Given the low values of S ν at low frequencies (S ν 10 −5 mJy for ν 100 MHz), detecting radio emission from V830 Tau's wind at these frequencies with present-day technology is not feasible (Section 5).
By comparing our radio emission estimates with the observations presented in Bower et al. (2016), we can place constraints on the mass-loss rates of V830 Tau. At certain observing epochs, Bower et al. (2016) detected V830 Tau in a flaring state, while at other observing epochs, only upper limits for radio emission of V830 Tau could be derived, namely, < 0.066 and < 0.147 mJy at 6 GHz (VLA observations) and < 0.117 mJy at 8.4 GHz (VLBA). These upper limits are shown as circles in Figure 4b. From that, we infer an upper limit ofṀ 3 × 10 −9 M yr −1 for the mass-loss rate of V830 Tau, which is consistent to the value used in our fiducial model (Section 2).

Can radio emission from V830 Tau b propagate through the stellar wind?
The planetary radio emission can only propagate in the stellar wind plasma if the (maximum) frequency of emission Ω c = eB(α 0 )/(2πm e c) is larger than the stellar wind plasma frequency ω p = [n e e 2 /(πm e )] 1/2 everywhere along the propagation path. In these expressions, n e is the local electron density of the stellar wind (n e = n/2 for a fully ionised hydrogen wind), e and m e are the electron charge and mass, respectively, and B(α 0 ) is the planetary magnetic field at colatitude α 0 (half-aperture of the polarcap boundary). In the scenario in which the planetary emission propagates towards wind lower densities, the condition Ω c > ω p is met when B(α 0 ) > n e 10 5 cm −3 where n e is taken as the electron density at the planetary orbit. Figure 5 shows the minimum planetary magnetic field intensity of V830 Tau b required for the propagation of planetary radio emission through the wind of the host star, as a function of the stellar wind mass-loss rates. Mass-loss rates of weak-lined T Tauri stars have not been observationally constrained. It is expected that these winds are intermediate between those of classical T Tauri stars with dense winds (Ṁ ∼ 10 −10 − 10 −8 M yr −1 , Hartigan et al. 1995;Iguchi & Itoh 2016) and zero-age main sequence stars with less dense winds. Very active, main sequence stars can have mass-loss rates ofṀ ∼ 2 × 10 −12 M yr −1 (Wood et al. 2005). In our work, we estimate an upper limit for the wind of V830 Tau to beṀ 3 × 10 −9 M yr −1 (Figure 4b). Given these values, we speculate that the most probableṀ for the winds of weak-lined T Tauri stars lies in the range ∼ 10 −12 -10 −10 M yr −1 . We use this range of values for the estimates we present next. 4 4 The non-detection of radio emission from stellar winds can be a useful way to observationally constrain mass-loss rates in solar-type stars, which have rarefied winds (e.g. Gaidos et al. 2000;Villadsen et al. 2014;Fichtinger et al. 2017). ForṀ ∼ 10 −12 -10 −10 M yr −1 , the size of the wind-emitting region ranges between R ν ∼ 3 -15 R at 275 MHz and between R ν ∼ 5 -30 R at 50 MHz (cf. Figure 4a). At these frequencies, R ν is comparable to the planetary orbital radius (6.1 R ), implying that V830 Tau b might be embedded in the regions of the stellar wind that are optically thick to radio wavelengths, although not so deeply embedded. Given that the radio flux we computed for the planet (∼ 6 -24 mJy, Figure 3b) is several orders of magnitude larger that the stellar radio flux (∼ 10 −6 -10 −3 mJy, Figure 4b), it is possible that, even after attenuation, a significant fraction of the planetary radio flux can escape. An accurate radiative transfer calculation of the amount of planetary flux that reaches the observer is left for a more detailed future study.
For the same range ofṀ, planetary radio emission can propagate through the stellar wind if the planetary magnetic field strength is 1.3 -13 G ( Figure 5). These minimum values encompass Jupiter's magnetic field intensities and appear to be physically reasonable. The equivalent (minimum) frequency of planetary emission is about 4 -40 MHz. At a frequency of ∼ 50 MHz (see previous paragraph) the wind is (radio) optically thick out to ∼ 5 to 30 R and the planet is not so deeply embedded in the radio-emitting wind. Altogether, these conditions seem very encouraging for planetary radio emission to escape the stellar wind of V830 Tau and be detected at Earth.

Discussion and Conclusions
Based on 3D MHD simulations, we presented here estimates of the radio emission expected for V830 Tau b, the youngest detected exoplanet to date. It orbits an active 2 Myr-old weak-lined T Tauri star at a distance of 0.057 au, which is about 6.1 stellar radii. At this distance, the environment surrounding this exoplanet is quite harsh. Using the observationally reconstructed stellar magnetic field, we simulated the wind of this star, allowing us to infer the conditions of the external ambient medium surrounding V830 Tau b. With these conditions, we then computed the planetary radio emission, using, as analogy, the radio-metric Bode's law derived for the magnetised planets in the solar system. According to this empirical law, there exists a relation between the dissipated power of the stellar wind impacting on the magnetosphere of the planet and the power released in the planetary emission. Our model uses simple approximations and should be considered as an initial attempt at calculating the radio flux of V830 Tau b, which is by far currently the best target for detection of exoplanetary radio emission. Detecting V830 Tau b at radio wavelengths would ultimately help us to refine our models.
We showed that, although the frequency of the radio emission is intimately related to the assumed magnetic field of the planet, the radio fluxes are only weakly dependent on that. The estimated flux densities from dissipated magnetic wind energy are on the order of 6 mJy, with peaks at 11 mJy, for an assumed planetary radius R p = R jup . Given the youth of V830 Tau b, it is likely that this is a lower limit for the planetary radius. Alternatively, for R p = 2R jup , the radio fluxes we estimated peak at 44 mJy, with average emission of 24 mJy. If V830 Tau b were to have a polar magnetic field intensity of 14.3 G (the maximum value of Jupiter's magnetic field), this means that the emission would occur at a frequency of about 28 MHz, originating at latitudes of about 34 o , where the magnetic field intensity is ∼ 10 G. Given that LOFAR sensitivity for a 1-h integration time at 20-40 MHz is 3-30 mJy (Grießmeier et al. 2011), we conclude that V830 Tau b is an excellent target to look for exoplanetary radio emission with LOFAR. Other present-day instruments that have the potential to detect such radio fluxes in the low-frequency range would be the upgraded UTR-2 (Ryabov et al. 2010) and GMRT (Sirothia et al. 2014). In the near future, the Square Kilometre Array (SKA)-low array system is expected to outperform LOFAR in terms of sensitivity limits for low frequency ranges (Zarka et al. 2015) and it would be an ideal tool for detecting exoplanets at radio wavelengths.
In Vidotto et al. (2012Vidotto et al. ( , 2015, the radio flux densities of several hot Jupiters (namely: τ Boo b, HD 46375b, HD 73256b, HD 102195b, HD 130322b, HD 179949b, and HD 189733b) were computed using the same method as here, that is, one in which the data-driven 3D MHD simulations provide the stellar wind characteristics used in the radio flux computations. Compared to the radio flux densities derived in those past works (0.02 to 1 mJy), V830 Tau b presents the best prospect for detecting radio emission (in spite of it not being the closest among the studied systems). This is because V830 Tau is an active star with a large-scale surface magnetic field intensity reaching 700G (Donati et al. , 2016, up to two orders of magnitude larger than the magnetism of the stars studied in Vidotto et al. (2015). The latter are main-sequence solar-type stars, which are significantly older than V830 Tau. In their youth, stars are more active and their magnetism more intense (Vidotto et al. 2014a).
In addition to the variability of radio emission seen along one planetary year (Figs. 2 and 3, see also Fares et al. 2010;Vidotto et al. 2010bVidotto et al. , 2012Vidotto et al. , 2015, variation within a few years timescale is also expected due to the evolution of the host star magnetism. This means that, within a few years, the conditions of the stellar wind surrounding the planet are likely to change, causing also changes in exoplanetary radio emission. This could be challenging for radio detection, which will likely require some monitoring of the system. The star itself may also contribute at radio frequencies (and potentially through flares that could make the planet modulation harder to detect). We computed thermal radio emission generated in the high density regions of the stellar wind. Comparing our estimated values with the non-detections reported in Bower et al. (2016), we were able to place a constraint in the mass-loss rate of V830 Tau:Ṁ 3 × 10 −9 M yr −1 . We argued that the most likely values ofṀ lie between ∼ 10 −12 and 10 −10 M yr −1 . With these values, the radio-emitting wind extends to distances of ∼ 3 to 30 R at frequencies of 50 to 275 MHz. This implies that the orbit of V830 Tau b (at 6.1 R ) might sit inside the radioemitting wind, but does not seem to be deeply embedded. Given that the radio fluxes of the planet are estimated to be ∼ 3 to 7 orders of magnitude larger than the radio flux of the stellar wind, it is possible that even after undergoing absorption by the stellar wind, a significant fraction of the planetary radio emission can escape the optically (radio) thick wind and reach the observer. Altogether, the V830 Tau system is a very encouraging system for conducting radio observations from both the perspective of radio emission from the planet as well as from the host star's wind. Equally important, the lack of detection of emission at radio wavelengths can place important constraints on the characteristics of the yet unobserved winds from weak-lined T Tauri stars.
Finally, exoplanetary radio emission may offer a new avenue for detecting young hot giants. The recent discovery of three of them (V830 Tau b, K2-33b, and TAP26b) may, in fact, indicate that these planets are more frequent than their equivalents around mature Sun-like stars. However, the extreme activity levels of their hosts prevents the use of traditional planet-search techniques. Radio searches can thus boost the findings of young hot Jupiters. It should also give us the chance to probe their magnetic fields from the detected radio frequency, and to simultaneously study stellar winds and star-planet interactions.
Acknowledgements. This work used the BATS-R-US tools developed at the University of Michigan Center for Space Environment Modeling and made available through the NASA Community Coordinated Modeling Center. This work was supported by a grant from the Swiss National Supercomputing Centre (CSCS) under project ID s516. JFD thanks the IDEX initiative of Université Fédérale Toulouse Midi-Pyrénées for awarding a Chaire d'Attractivité in the framework of which the study of V830 Tau was carried out. The authors would like to thank Manuel Guedel and the anonymous referee for their constructive comments, which greatly improved this manuscript. into radio power. On the other hand, for close-in planets orbiting around highly magnetised stars (in the magnetic-pressure dominated region), the radio flux becomes roughly independent of the density, being mainly affected by the local conditions of the stellar magnetic field and the relative velocity of the planet through the stellar wind (Eq. B.4). Since the latter condition is the most likely for V830 Tau b, we conclude that the stellar wind density, an unknown in our numerical simulations, plays a minor role in our estimates of planetary radio fluxes.