Infrared-Radio-follow-up Observations for Detection of the Magnetic Radio Emission of Extra-Solar Planets: A New Window to Detect Exoplanets

There are several methods for indirectly detecting exoplanets, such as transit, radial velocity, astrometry, and the conventional gravitational microlensing approach. These methods rely on observing the effects of exoplanets on the emission or motion of observed stars. All these techniques have focused on the optical or infrared domains. However, an alternative method for exoplanet detection via microlensing events involves planets orbiting the source star, creating a binary source system. In this study, we explore a novel approach to detecting and studying exoplanets exclusively through their radio emissions resulting from magnetospheric processes. We propose utilizing the Roman telescope as a survey observer to detect microlensing events. Subsequently, we investigate the potential for detecting planetary radio signals through follow-up observations of these microlensing events in the radio band using the SKA telescope. This method is viable due to the comparable radio emission levels of exoplanets and their parent stars, unlike optical and infrared emissions. We conduct a Monte Carlo simulation to replicate the observations by the Nancy Roman Telescope, followed by a follow-up observation in radio frequencies using the SKA telescope. We determine that approximately 1317 exoplanets exhibit detectable signals by the SKA telescope during the 7-season observations by the Nancy Roman Telescope. This result indicates that such a method cannot only facilitate the direct detection of exoplanets but also enable the measurement of their magnetic field strength through analysis of their radio emissions.


INTRODUCTION
A planetary magnetic field is one of the few observable quantities that contain information about the properties of the planetary interior since the internal dynamo produces it.Further, the extent to which exoplanets retain their atmospheres over geological time scales depends on whether they are exposed directly to stellar winds.Also, planetary magnetic fields could potentially serve as indicators of plate tectonics.Therefore, information about the magnetic fields of extrasolar planets is essential to assess their habitability.So far, our knowledge of planetary magnetic fields is limited to our solar system's planets, which are not representative of the vast array of possible planets, such as hot Jupiters (Lammer et al. 2009;Lingam & Loeb 2021;Cockell et al. 2016;Kane 2021).
But how can we observe and measure the magnetic field of exoplanets?The magnetospheres of Earth, Jupiter, Saturn, Uranus, and Neptune all show evidence of nonthermal continuum radiation emission resulting from the solar wind's interaction with the planet's magnetic field (Zarka 1998).The electron cyclotron maser instability (ECMI) is a phenomenon that occurs in plasmas where electrons emit coherent electromagnetic radiation due to confinement by magnetic field lines along which they travel, gyrate, and accelerate (Melrose & Dulk 1982;Dulk 1985).The ECMI is the most efficient radio emission mechanism, dominating all other emissions (beam-plasma instabilities and others).Based on our knowledge of the solar system planets, the average ECMI is related to the total solar wind power and the magnetic field of the planet by using the Radiometric Bode's Law (RBL) model (Zarka et al. 2001;Zarka 2004Zarka , 2007)).
Therefore, observations of this emission type can indirectly measure and classify the planetary magnetic field.This quantification of the magnetic field gives us constraints on internal structure models and planetary rotation, informing us about the evolutionary history of the planet and host system (Hess et al. 2011).Furthermore, the exoplanets' radio emission may reveal the presence of exomoons around them.For example, it has been shown that part of the Jovian decameter radio emission is "controlled" by the satellite Io (Bigg 1964).Noyola et al. (2014) argues that the Io-Jupiter interaction could also be found in exoplanet-moon pairs, and the exoplanets' radio emissions could be used to detect such systems directly.
A multitude of exoplanet detection techniques rely on indirect optical/infrared wavelength observations (Wright et al. 2012;Knutson et al. 2014;Mandel & Agol 2002;Seager & Hui 2002;Traub & Oppenheimer 2010;Guyon et al. 2005;Gaudi 2012;Bagheri et al. 2019).However, measurements of the planetary magnetic fields require direct observations at the low-frequency radio wavelengths.Despite efforts to detect exoplanets in the radio band via ground-based radio telescopes, no promising detection of planetary radio emissions has been successful (Murphy et al. 2015;Lazio et al. 2018;Turner et al. 2021), primarily due to the low sensitivity of ground-based radio telescopes.Consequently, direct confirmation of the existence of exomoons has so far remained impossible.
In this work, we combine two different exoplanet detection methods to make the radio emission of exoplanets detectable.The novel idea behind this study is that during a microlensing event, light intensity will be magnified at all frequencies.Therefore, a microlensing event can serve as a natural light magnifier.The idea is to use an optical/infrared telescope to detect a microlensing event that would trigger follow-up observations by radio telescopes upon detection.If the source star has at least one exoplanet in microlensing events, we can directly measure the light emitted from the exoplanet (Figure 1).The magnified light of the exoplanet could be observable in the infrared wavelength as we have shown in (Bagheri et al. 2019).However, the radio band is the most desirable wavelength window to detect exoplanets (Vedantham et al. 2020;Viswanath et al. 2020;Foster et al. 2020) since the planetary radio emission, which includes the effects of exomoons, may dominate the signal from the parent star (Cauley et al. 2019;Zarka 2004).This method will enable us to measure the magnetic fields of exoplanets and provide us with information about the existence of exomoons.To show this, we simulate the Nancy Grace Roman Space Telescope (previously named Wide Filed Infrared Space Telescope, WFIRST) as an infrared survey telescope to alert the Square Kilometre Array (SKA) radio telescope to the occurrence of microlensing events for follow-up observations.2. METHODOLOGY

Microlensing events
Microlensing, a powerful tool in the study of astrophysics, involves the gravitational bending of light by massive objects, typically stars, within the lensing galaxy.This phenomenon, predicted by Einstein's theory of general relativity, manifests as a temporary increase in the brightness of a background source, such as a distant star, as the foreground lensing object passes in front of it.The gravitational field of the foreground object acts as a lens, magnifying and distorting the light from the background source.In lensing, the magnification, denoted by A, represents the factor by which the flux of the source is amplified during the microlensing event.Therefore it is given by the ratio of the image size to the source size: For point source magnification in microlensing, where the source size is negligible compared to the Einstein radius (R E ), the magnification can be simplified as: where u is the impact parameter, defined as the closest approach of the source to the lens normalized by the Einstein radius (Gaudi 2012).The lensing system may comprise either a single object (single-lens) or two objects (binary-lens), with the possibility of additional components (Daněk & Heyrovskỳ 2019).The presence of two lenses introduces additional complexity to the gravitational lensing phenomenon, leading to more intricate and often asymmetric light curves.Source systems can also consist of one or more objects.The conventional technique utilized in gravitational microlensing for exoplanet detection relies on the premise that the lens star hosts an exoplanet, constituting a binary lens system.Alternatively, another microlensing approach involves a scenario where the planet orbits the source star, resulting in a binary source configuration (Bagheri et al. 2019) as illustrated in Figure 1.In this study, we simulate binary-source microlensing events by assuming both single-lens and binary-lens configurations.We consider the finite source effect for the source star while neglecting limb darkening.

Source and Lens Stars Parameters
To characterize the star mass density within the Milky Way galaxy, we rely on the Bésançon model, extensively discussed in previous works (Spergel et al. 2015;Robin et al. 2003;Gardner et al. 2014).This model presents the distribution of matter in our galaxy as a superposition of eight thin disk structures with varying ages, a thick disk component, and a central (old) bar structure consisting of two components.Our choice is the updated model proposed by (Gardner et al. 2014).We select the model parameters fitted to a two ellipsoid bar.This refined model, utilizing star distribution data from the Hipparcos catalog, has proven effective in interpreting microlensing data from the EROS collaboration in the direction of the spiral arms (Moniez et al. 2017).When generating binary lenses, the mass ratio between the two lenses is determined from the distribution function proposed by (Duquennoy & Mayor 1991).We choose the semi-major axis s of the binary orbit of microlenses from the Öpik's law where the distribution function for the primary-secondary distance is proportional to ρ(s) = dN/ds ∝ s −1 in the range of [0.6, 30] au.

Exoplanet Parameters and their Occurrence
In our simulations, we categorize the characteristic parameters of exoplanets into two primary groups: rocky planets and Jovian planets.According to the conventional planet formation theory, rocky planets are situated within the snow line of their parent star (Kennedy et al. 2006;Kennedy & Kenyon 2008).The snowline of a star, also known as the frost line or ice line, refers to the distance from the star where volatile compounds, such as water, can condense into solid form due to the decrease in temperature.Beyond the snowline, the temperature drops low enough for these volatile compounds to freeze, forming icy grains or particles (Kennedy & Kenyon 2008).This snowline can be expressed with the snowline of the Sun as where M ⋆ is the mass of the parent star (Gould et al. 2010).For Jovian planets, as they can migrate towards their parent star from distant orbits (as discussed by (Murray et al. 1998;Papaloizou & Terquem 2005), we set the maximum distance from the parent star at 10 au.This range is consistent with observations from MOA microlensing studies, indicating at least one bound planet per star within this range, represented as a = 0.01 − 10 au (Sumi et al. 2011).We assume that the distribution of exoplanet semi-major axes follows the Öpik's.
The planetary mass range in our simulation is 0.002 − 10M J (Schlaufman 2018).For planets with masses > 100 M ⊕ orbiting host stars with masses greater than 0.6 solar masses, we utilized the correlation between giant planet occurrence and host star mass and metallicity as outlined by (Johnson et al. 2010;Fulton et al. 2021).Following the occurrence dependencies discussed in (Fulton et al. 2021;Santerne et al. 2016,?;Wittenmyer et al. 2020;Winn & Petigura 2024), we adopted their result of planet occurrence.We summarized the frequencies of various exoplanet types from our simulation in Table 1.
Generally, exoplanets exhibit a wide range of mass-radius relationships, reflecting their diverse compositions and evolutionary histories.However, a commonly used empirical relationship for estimating the radius of exoplanets based on their mass and incident flux is the following power-law scaling relation where F is the bolometric incident flux from the parent star and can be calculated by where R ⋆ is the stellar radius, T eff is the effective stellar temperature, a is the semi-major axis, and e is the orbital eccentricity (Weiss et al. 2013).1. Summary of the planet occurrence frequencies in our simulation, as a function of planet mass and period.

Estimation of Exoplanets Radio Emission
To simulate the detection process involving the Roman telescope for surveying microlensing events and subsequent follow-up observations with the SKA telescope, a model is necessary to estimate the emission of exoplanets based on their characteristics and the attributes of their parent stars.In general, the flux received from a planet contains thermal radiation due to its intrinsic temperature, as well as the reflection radiation from the parent star and radio emission if it has a magnetic field due to the electron cyclotron maser instability (ECMI) (Melrose & Dulk 1982;Dulk 1985).The thermal and reflection radiation can be readily computed from the planet's characteristics and its parent star (Bagheri et al. 2019).The ECMI radiation is characterized by strong circular polarization and high anisotropy.ECMI stands out as the most efficient mechanism for radio emission, superseding other mechanisms like beam-plasma instabilities.Consequently, it plays a dominant role in the generation of radio emissions.This dominance is exemplified in the brightness of auroral radio emissions from planets in the solar system, particularly Jupiter, which rival the intensity of solar radio bursts.
Based on observations of the magnetized Solar system planets, we know that all auroral radio emissions power is related to incident energy flux of the stellar wind (Zarka et al. 2001;Zarka 2004Zarka , 2007;;Zarka et al. 2018).In all cases in our solar system, about 2 × 10 −3 of the electrons' energy goes into radio waves (Zarka et al. 2018).Therefore, the average ECMI radiation can be linked to both the total power of the solar wind and the magnetic field strength of the planet, utilizing the RBL model (Zarka et al. 2001;Zarka 2004Zarka , 2007)).The gyrofrequency can be expressed mathematically as: where ν is the emitted pick frequency, e is the charge of the emitting particles (in this case, electrons), B is the magnetic field strength, and m is the mass of the emitting particles (electron mass in this case).This equation shows that the emitted frequency is directly proportional to the magnetic field strength of the planet.Therefore, observations of this emission type can provide an indirect means of measuring and classifying the planetary magnetic field.This quantification of the magnetic field gives us constraints on internal structure models and planetary rotation, informing us about the evolutionary history of the planet and host system (Hess et al. 2011).
The RBL model is frequently used to estimate the radio emission of extrasolar planets (Farrell et al. 1999;Zarka et al. 2001;Joseph et al. 2004;Zarka 2007;Christensen 2010).The simplicity of this model has some advantages, but it does not provide a complete picture of all processes involved, in particular for close-in exoplanets as discussed in (Bagheri et al. 2024b) (for Examples of studies that extend beyond the RBL model for calculating radio emission (Vidotto & Donati 2017;Li et al. 2023;Dong et al. 2018;Vidotto 2017)).However, the RBL model facilitates comparative analyses between the radio emissions of exoplanets and those observed within the solar system, thereby enhancing our contextual understanding of planetary magnetospheric dynamics.This capacity for comparative analysis underscores the utility of the RBL model in advancing our comprehension of exoplanetary environments and their potential habitability while also serving as a foundation for future investigations into the broader astrophysical implications of planetary radio emissions.Thus, in this paper, we follow the approach outlined by (Joseph et al. 2004) to estimate the radio flux density of exoplanets.As described by (Farrell et al. 1999), the radio power from an exoplanet is related to the incident power of the stellar winds.Following (Farrell et al. 1999), the Jovian decametric component is considered at least partly related to the solar wind kinetic energy input and is used as the base power level to give where a is the semi-major axis of the planet's orbit in au, Ṁ ion is the stellar ionized mass-loss rate, v ∞ is the terminal velocity of the stellar wind, and ω is the corotation speed.For the terminal velocity of stellar wind, we use where R ⋆ and M ⋆ are the host star radius and mass in solar radius and mass units (O'Gorman et al. 2018).Equation ( 8) serves as a mere approximation for low-mass stars, however, recent study (Mesquita & Vidotto 2020) indicates that the terminal velocity of M dwarf stellar winds either aligns well with predictions from Parker wind models, in which stellar winds are driven by radiation pressure mechanism, or exceed these model estimates.As a result, employing this equation may yield conservative estimations of stellar winds terminal velocity and the overall count of planets with detectable signals.Then, the radio flux density (S ν ) emitted by an exoplanet can be estimated using the following equation where P rad is the total radio emission power of the exoplanet (7), d is the distance between the exoplanet and the observer, and ν is the frequency of the radio emission, where R p is the planetary radius.Substituting (9) into (7) gives where ω J , M J , and R J , are Jupiter's corotation speed, mass, and radius and Ω is the beaming solid angle of the emission.In deriving this expression, we have followed (Farrell et al. 1999) and assumed that the planet will emit ECMI emission between the frequencies 0.3ν c and ν c , where ν c is the maximum radiation frequency.
The planetary corotational speed (ω) is a crucial parameter in estimating radio flux from exoplanets.Exoplanets orbiting closely around their host stars experience significant tidal dissipation, resulting in considerably slowed corotation or complete tidal locking.This phenomenon has a profound impact on the production of auroral radiation due to the permanent orientation of one side of the planet towards the star, causing asymmetrical ionospheric conductance (Zarka et al. 2001;Seager & Hui 2002).To consider the effect of tidal locking, we assume all the exoplanets with orbits < 0.1 au are tidally locked, so the orbital period is equal to their corotation period.For all other exoplanets, we use the Darwin-Radau relation to relate corotation and oblateness (Murray & Dermott 2000) where C is 0.4 for rocky planets and for gas giant planets C = 0.25 (Hubbard 1984).In this equation, f is the planet's oblateness; the oblateness of planets refers to the degree to which a planet deviates from a perfect sphere, taking on a slightly flattened shape due to the centrifugal force generated by its corotation.This deviation from sphericity is more pronounced for rapidly rotating planets.The oblateness of a planet is quantified by where R eq and R pole are the equatorial and polar radii of the planet.In our simulation, for the Jovian planet with orbital distance a > 0.1 au, we use the oblateness of Jupiter (f = 0.064), and for rocky planets, we use the Earth's value which is f = 0.00335 (Barnes & Fortney 2003).In Our simulation we also assume no planets rotate faster than Jupiter.
In this paper, the arrangement considered for the lens, source star, and planet entails that the observed light curve of microlensing events results from a combination of two distinct light curves: one originating from the source star and the other from the planet.Consequently, the final light curve results from merging these two light curves.To calculate host stars' radio flux at the desired frequency, we use the relationship between the luminosity of the host star in the X band and the radio band.Stars emit X-rays primarily through high-energy processes such as magnetic reconnection events, coronal heating, and interactions between fast-moving charged particles and the stellar atmosphere.The X-ray luminosity of a star is thus influenced by its magnetic activity level, mass, age, and evolutionary stage.
On the other hand, radio emission from stars is more associated with non-thermal processes, such as ECMI.The strength of radio emission depends on factors including the density and strength of the stellar magnetic field, the presence of energetic particles, and the efficiency of acceleration mechanisms.While there is no universal scaling relationship between X-ray and radio luminosity for all types of stars, correlations have been observed within specific stellar populations, such as active stars with strong magnetic fields or those undergoing rapid evolution.One can use empirical relationships based on observations for the radio luminosity of main sequence stars.One such relationship is the Gudel-Benz relation, which relates the X-ray and radio luminosities: where L X is the X-ray luminosity and L R is the radio luminosity (Guedel & Benz 1993).Notably, ECMI radiation exhibits strong circular polarization, whereas solar or stellar plasma radiation lacks polarization and occurs sporadically.Additionally, ECMI radiation is beamed anisotropically, resulting in significant modulation by the planetary rotation (Zarka et al. 2014).Therefore, polarization and temporal variations should enable the differentiation between stellar and exoplanetary radio emissions, as gravitational lensing does not change the polarization of the light rays.

INDUCED RADIO EMISSION BY EXOMOONS
There are more than 200 moons in our solar system, many of which orbit gas giants like Jupiter and Saturn.Despite significant efforts, the successful detection of exomoons remains elusive (Narang et al. 2023).If an exomoon resides within the magnetosphere of its host planet, it can influence the planet's radio emissions by injecting plasma that accelerates electrons within the magnetosphere.A notable example is Io, a moon characterized by intense volcanic activity within Jupiter's magnetosphere.This activity generates a tenuous atmosphere around Io (Lopes et al. 2007), contributing ions to Jupiter's magnetosphere and generating a plasma (Reich et al. 2009).These ions in motion accelerate electrons, ultimately leading to the formation of an ECMI along the field lines, resulting in radio emissions known as Io-DAM (see Figure 1-a in (Cecconi 2014)) (Crary 1997;Mauk et al. 2001).Similar effects are observed with other Galilean moons and moons orbiting Saturn, such as Titan and Enceladus.Consequently, it is reasonable to anticipate moon-planet interactions in other stellar systems (Green et al. 2021).ECMI radio emissions from exomoon magnetosphere have not had a confirmed detection, but there are three candidate systems under study (Turner et al. 2021).
To simulate the induced radio emission by exomoons, one should know the probability of having a moon such as Io.The results of N-body simulations in (Heller et al. 2014;Sasaki et al. 2010) show that 80% of Jovian planets have moons with masses M m /M p ≈ 10 −5 .However, stellar tidal forces serve to slow down a planet's rotation and in conjunction with tidal migration, contribute to removing moons; At short a, migration of moons is significant and occurs quickly, which can lead to orbital instabilities (Sucerquia et al. 2020).This process limits the lifetimes of larger moons in extrasolar planetary systems (Barnes & O'brien 2002).Therefore, in our simulation, we assume 80% of Jovian planets with masses M p ≥ 0.1 M J with a ≥ 0.1au have one moon.So exoplanets with the semi-major axis of less than 0.1 au do not have any moons, which is consistent with assuming they are tidally-locked.Noyola et al. (2014) proposed that interactions akin to those between Io and Jupiter's magnetosphere could also occur in exoplanet-moon pairs.The frequency of the moon-induced radio emission is slightly different, where ν ci is the frequency of radio emission induced by a moon, such as Io, a m is the same-major axis of the moon's orbit, and ν c can be obtained by (10).Equation ( 16) is asymptotic, approaching a maximum value of ν c as a m becomes very large (Noyola et al. 2016).Assuming the plasma density of the exomoon is equal to the Io's plasma density, the optimal distance for the exomoon to have the maximum of the induced radio emission can be obtained by (Noyola et al. 2014) Therefore, the total radio flux induced by a moon is where d is the distance of star-planet-moon system to Earth, ρ m is the moon's plasma density, R m is the moon's radius, and β is the efficiency coefficient of converting Joule dissipation to radio waves (we assume β ≈ 1% as is in the Jupiter-Io system (Zarka et al. 2001)).B m is the magnetic field affecting the exomoon, and so , where B p is the magnetic field of the exoplanet at its poles.The plasma speed, V 0 , is computed assuming it corotates with the planet's magnetic field.Therefore, where G is the gravitational constant, M p is the planet's mass, and ω p is the planet's corotation speed.Finally, ρ c is the critical density which is defined as Noyola et al. 2014).In our simulation, we assume IO's plasma density and IO's radius for ρ m and R m .We also assume Ω = 0.2 (Noyola et al. 2014).In our simulation, the induced radio emission by an Io-like moon accounts for a maximum of 0.2% of the total radiation emitted by the exoplanets.

Observation Strategy
The NASA Roman Space Telescope is a next-generation space observatory designed to investigate various astrophysical phenomena using infrared observations.Named in honor of the astronomer Nancy Grace Roman, often regarded as the "Mother of Hubble" for her pivotal role in the development of the Hubble Space Telescope (Roman 2019), this advanced instrument is slated to launch in the mid-2020s.With its wide field of view, the telescope is designed to conduct large-scale surveys to explore dark energy, dark matter, and the formation and evolution of galaxies with unprecedented detail (Rose et al. 2021).The Roman Space Telescope will also be crucial in advancing our knowledge of exoplanets using microlensing events.Roman W149 filter (0.9 − 2 µm) will have a significantly low detection threshold, with a zero-point magnitude of 27.61.The microlensing survey of Roman will monitor 1.97deg 2 of the Galactic bulge (in the direction of b = −1.5 and l = 0.5 in the Galactic coordinate) with 15-minute cadence, over six 72-day per season (Spiegel et al. 2005), potentially detecting thousands of exoplanets via the perturbations that they produce on the microlensing light curves (Bagheri et al. 2019).Since microlensing events usually have a duration of a few days up to a few weeks, we propose to use the Roman telescope as a survey observer to detect a microlensing event and then alert that event to a ground-based radio telescope for a follow-up observation of that microlensing event in the radio band.The most suitable ground-based telescope to detect planetary radio emission is the Square Kilometre Array (SKA) telescope.As the world's largest radio telescope, the SKA comprises two arrays: SKA1-Low, optimized for low-frequency observations between 50 MHz and 350 MHz, and SKA1-Mid, designed for mid-frequency observations between 350 MHz and 14 GHz.The vast collecting area of the SKA, totaling one square kilometer when fully operational, enables it to detect faint radio signals from cosmic phenomena spanning a wide range of scales, from nearby planets to distant galaxies and beyond.Furthermore, the SKA's innovative design incorporates advanced signal processing techniques and data analysis algorithms, maximizing its scientific output while minimizing data volumes and processing requirements.

Detectability of the Planetary Radio Signal
The detectability of the radio signal of exoplanets depends on, first, the detection of the microlensing event by the Roman telescope, and, second, the frequency and the intensity of the radio emission of the star-planet system.To address the first condition, we calculate the host star and planet magnitude in Roman's W149 filter.For the flux received from the host stars, we consider both the distance-dependent effects of the source stars and the reddening of their apparent magnitudes due to interstellar dust.Using the host stars' position and stellar type, we estimate the extinction caused by interstellar dust along the line of sight, using the comprehensive 3D extinction map provided by (Marshall et al. 2006).We use the relationships outlined in (Nishiyama et al. 2008(Nishiyama et al. , 2009) ) to convert the extinction values from V bands to the band of W149 filter (Bagheri et al. 2019;Penny et al. 2019).The flux received from a planet comprises both thermal radiation originating from the planet's intrinsic temperature and reflected radiation from its parent star.The planet's temperature can be calculated by assuming that the planet's thermal emission adheres to blackbody radiation principles (López-Morales & Seager 2007).This calculation incorporates the absorption of radiation from the parent star by the planet, which is then re-emitted following Boltzmann's law, where T p indicates the temperature of the planet, A B denotes the albedo (0.15 and 0.52 for rocky and Jovian planets), while f characterizes the proportion of re-radiated energy absorbed by the planet (Harrington et al. 2006;Knutson et al. 2007).Since hot Jupiters are the most frequently detected, we adopt a value of 2/3 for f across all planets.Assuming blackbody radiation, the thermal flux (F th ) at the specified frequency ν can be expressed using Planck's law as follows: where, I(ν, T p ) denotes the emitted power per unit area of the emitting surface, per unit solid angle, and per unit frequency.Integrating it across the energy range of the W149 filter and over a hemisphere yields the thermal flux of the planet as observed by Roman.The flux of the planet resulting from the reflection (F ref ) of the parent star's light is determined by where g(Φ) denotes the fraction of the illuminated portion of the planet visible to the observer, and A g represents the geometric albedo.We utilize a geometric albedo formula: A g = 2/3 A B (Bagheri et al. 2019).While this relationship may not hold true for all phase angles, its alteration does not significantly impact results.This is because the thermal radiation emitted by planets is much weaker than the flux of host stars in the Roman W149 filter.Thus, any changes in this relationship would only affect the total flux received by the Roman telescope at the order of 10 − 4 at most.Consequently, the overall flux received from the planet (F p ) can be calculated as the combined sum of the thermal flux and the reflection flux.In the final stage, we add realistic noise.The total noise consists of the intrinsic Poisson fluctuation in the flux received from the microlensing event in a single exposure, the intrinsic Poisson fluctuation within the point spread function (PSF) representing the background sky flux in a single exposure, the read-out and dark noise (Bagheri et al. 2019).Finally, we determine the range within which the noisy signal exhibits magnitudes below the zero-point magnitude of Roman and above the saturation magnitude of 14.8 (Penny et al. 2019).
For the follow-up observation by the SKA, we consider the SKA1-low and SKA1-mid frequency ranges up to 890 MHz.The image noise is given by, where S D = 2.5 is a degradation factor relative to the natural array sensitivity for the specific target Gaussian FWHM resolution of the image, η S = 0.9 is a system efficiency that takes account of the finite correlator efficiency, and η pol = 2 is the number of contributing polarizations and for SKA1 with orthogonal linear polarizations (Braun et al. 2019).
The imaging sensitivity for SKA1-Low and SKA1-mid up to 890 MHz is listed in  2. Image sensitivity of SKA1-Low and SKA1-mid within the indicated frequency bands for continuum observations (Braun et al. 2019) We use the SKA image sensitivity for different ranges of frequencies as summarized in Table 2. Therefore, any event with flux density greater than the SKA sensitivity for at least 4 hours of imaging is considered a detected event.

RESULTS
We simulate 150,000 single-lens and binary-lens microlensing events.Our Monte Carlo simulations indicate the likelihood of detecting radio emission of rocky planets using either single-lens or binary-lens configurations is exceedingly low with the SKA telescope.Consequently, Jovian planets exhibit the highest detection probabilities in both scenarios.This outcome primarily stems from their stronger magnetic field and higher frequency of radio emission relative to rocky planets, resulting in a substantial portion of their radiation falling within the SKA sensitivity.Simulations of Roman-SKA data predict, under the assumption of single-lens configuration, that the number of detecting a Jovian planet with orbital periods of 2 days to 10 years is about 3.54 %.For the binary-lens simulation, this number increases to 8.25 % due to the high magnification of caustic crossings.
Figure 2 shows our simulation's distributions of planetary parameters.Exoplanets with detected radio fluxes are in orange.As depicted in Figure 2-a, it's notable that the majority of exoplanets exhibiting detectable signals are situated at orbital distances lower than 2 au.However, for planets with higher masses, the orbital distance can extend up to 8 au.The planetary radio flux in terms of magnification of the microlensing events is shown in Figure 2-b.Exoplanets within the shaded region exhibit emission frequencies below 50 MHz, rendering their signals undetectable despite their radio emissions being of high intensity and possessing a significant magnification factor during microlensing events.Figures 2-d and 2-e represent the distance of the star-planet system relative to the exoplanet mass and the magnification of the microlensing events.Notably, the range extends up to 10 kpc for exoplanets with detectable signals, underscoring the importance of this observation method.These exoplanets are beyond the reach of current technologies for detection in any other wavelength bands using conventional methods.Identifying such star-planet systems can enhance our understanding of planet formation rates within the galactic bulge and facilitate the detection of extragalactic exoplanets.
The detection efficiency in terms of mass and semi-major axis of the exoplanets is illustrated in Figure 3. Detected planets exhibit masses ranging from M p = 1 to 10M J and semi-major axes spanning from a = 0.01 to 10 au.
We conducted simulations targeting coordinates l = 0.5 • and b = −1.5 • .In this direction, as reported by (Penny et al. 2019), the estimated number of source stars detectable by the Roman telescope with magnitudes below 25 is 240 × 10 6 .Hence, the number of microlensing events that could be detected via Roman observation can be given by: where T obs is the observational time, N ⋆ is the number of background stars detectable by Roman, ϵ(t E ) is the Roman efficiency for detecting microlensing events with the duration of t E , and τ is the optical depth (Bagheri et al. 2019).
The optical depth, averaged over the distances of the source stars, is τ = 2.4 × 10 −6 (Mróz et al. 2019), and our simulation yields an average value of < ϵ/t E >= 0.035.Therefore, around 27000 microlensing events can be detected by Roman observation (Penny et al. 2019).Consequently, approximately 1317 exoplanets could be detected with Roman-SKA observations toward the Galactic bulge.However, it may not be feasible to follow up on all Roman's detected microlensing events; nonetheless, a 10-20% follow-up of these events could significantly impact exoplanet sciences and the measurement of exoplanetary magnetic fields.These findings are summarized in Table 3.Our comprehension of planetary magnetic fields remains largely confined to the solar system.Expanding our knowledge and observation of exoplanetary magnetic fields could provide valuable insights into their internal compositions and dynamics, which are otherwise challenging to ascertain.The detection of radio emissions from confirmed exoplanets has been investigated in multiple studies (e.g., Bagheri et al. (2024a); Li et al. (2023).In this study, we investigate an alternative approach to detecting exoplanets through their radio emissions magnified by microlensing events.Our findings reveal that the planetary radio emissions of approximately 5% of microlensing events could be detected through follow-up observations with Roman-SKA (equivalent to ∼ 8% for binary-lens and ∼ 3.5% for single-lens events).This finding gains significance when considering that, to date, only 210 exoplanets have been detected via microlensing events1 .These detections primarily occurred due to caustic crossings in binary-lens events, making follow-up observations of these planets nearly impracticable.In contrast, we consider exoplanets as companions within the source system of microlensing events, enabling us to conduct follow-up observations on the star-planet system.
Our results underscore the potential of infrared-radio follow-up observations to enhance the rates of direct exoplanet detections while yielding valuable insights into the magnetic fields and internal structure of exoplanets.The observation strategy outlined in this study demonstrates heightened sensitivity to detecting the radio emissions of giant and close-in exoplanets.Observing the radio emissions of close-in exoplanets broadens our comprehension of planetary responses to their space environment and evolutionary trajectories.Hot Jupiters' magnetospheres serve as compelling subjects for study, highlighting the intricate complex interplay between magnetic fields, stellar wind, and atmospheric escape processes.Observing radio emissions of close-in exoplanets and further investigations into these magnetospheric properties will deepen our understanding of exoplanetary systems and their diverse environments.Such insights will significantly contribute to our comprehension of planet formation, migration, and the potential habitability of exoplanets.
The outcomes of this project are important for the prospects of Exoplanet radio emission detection by the contemporary and planned surveys as well as the next generation of Radio telescopes.In a wider view, the results of this work pave the way toward actual infrared-radio follow-up observations of exoplanets/exomoons using the method introduced in this study.The knowledge of the magnetic field and the magnetospheric emissions of exoplanets are considered a science frontier for the next decade, as highlighted in the Origins, Worlds, and Life Planetary Science & Astrobiology Decadal Survey report.Two Priority Science Question Topics identified in this survey report incorporate aspects of planetary magnetic fields and their interplay with the solar wind.These are identified as Q6: "Solid body atmospheres, exospheres, magnetospheres, and climate evolution" and Q12.7: "Exoplanets, Giant planet structure and evolution."(NationalAcademies of Sciences et al. 2022).

Figure 1 .
Figure 1.A schematic view of a binary-source microlensing event.

Figure 2 .
Figure 2. Characteristics of the exoplanets within the simulated event.Exoplanets exhibiting detectable radio signals are highlighted in orange; a) distribution of semi-major axis as a function of planetary mass, b) planetary radio flux vs. magnification of the microlensing events, c) planetary radio flux as a function of the frequency of the planetary radio emissions, d) distance of the star-planet system vs. planetary mass, and e) distance of the star-planet system vs. magnification of the microlensing events.

Figure 3 .
Figure 3.The detection efficiency of the planetary radio emission in terms of the mass and semi-major axis of the exoplanet in our Monte Carlo simulation.
Table 2 by assuming a continuum observation with fractional bandwidth of ∆ν/ν c ≈ 0.3, together with an integration time ∆τ = 1 hour.