ROME. III. The Arecibo Search for Star–Planet Interactions at 5 GHz

After nearly three decades of discovery, many exoplanetary systems have been studied and characterized in detail with one important exception: exoplanet magnetism. Although many surveys sought to detect magnetospheric radio emissions from exoplanets to directly measure their magnetic field strengths, they have yet to reveal an unambiguous detection. However, the indirect detection of exoplanet magnetic fields by measuring their influence on their host stars via magnetic star–planet interactions has recently gained prominence as an alternative method of discovery. This third paper of the Radio Observations of Magnetized Exoplanets series presents the results of a targeted radio survey of eight nearby exoplanet-hosting systems that may engage in star–planet interactions. This survey, conducted with the Arecibo radio telescope at ∼5 GHz, has the greatest frequency coverage of any to date while providing millijansky-level sensitivity over <1 s integration times. No exoplanet-induced stellar radio bursts were detected. The orbital phase coverage of candidate systems for magnetic star–planet interactions is described, and the survey results are examined within the context of the plasma flow–obstacle paradigm and searches for star–planet interactions at other wavelengths.


INTRODUCTION
Stars interact with their orbiting exoplanet companions via a variety of mechanisms, but perhaps none is more controversial than star-planet interactions (Lanza 2018).Cuntz et al. (2000) hypothesized that tidal and magnetic interactions between host stars and exoplanets with semimajor axes a 0.1 au might increase stellar activity in several ways.Tidal forcing of the stellar atmosphere and convective zone may increase local or global magnetic activity, while unipolar or dipolar magnetic plasma flow-obstacle interactions may enhance activity in stellar atmospheres (Zarka 2018).Magnetic and tidal effects should be distinguishable by the presence of either one or two enhancements in activity per stellar rotation period, respectively.Both tidal and magnetic activity enhancement should be most pronounced in the Corresponding author: Matthew Route mproute@olemiss.educorona and transition region, and to a lesser extent the chromosphere, due to their comparatively lower densities and proximity to the interacting exoplanet companion (Saar & Cuntz 2001).
Observational support for these interactions originally came from a chromospheric enhancement monitoring campaign of hot jupiter host stars including HD 179949, υ And, and HD 189733 (Shkolnik et al. 2005(Shkolnik et al. , 2008)).Several data sets of Ca II K residuals appeared to show activity enhanced by ∼1%, but since the modulation varied from epoch to epoch and typically did not occur at the orbital phases expected from initial theoretical models, a more complicated energy transport model was required to explain the intermittent activity.Later calculations found that the detected chromospheric activity was several orders of magnitude greater than is available from plasma flow-obstacle interactions (e.g., see Zarka 2007 on HD 179949) calling such detections into question.Similarly, statistical analyses of stars with hot exoplanet companions came to conflicting conclusions, with some failing to find a statistically significant correlation between exoplanets with small semimajor axes and chromospheric activity enhancement (e.g., Canto Martins et al. (2011); Poppenhaegar & Wolk (2014)).Rubenstein & Schaefer (2000) hypothesized that asyet-undetected planetary companions could cause the emission of superflares (E ∼10 33 -10 38 erg) from F and G main-sequence stars.They argued that the properties of the flares emitted from κ Ceti and π 1 UMa were similar to those emitted by RS Canum Venaticorum (RS CVn) stars.They proposed that open field lines from the host star were anchored at the exoplanetary companion, which would wrap the field lines about the host star as the exoplanet orbits it, thereby amplifying the stored magnetic energy.Subsequent interactions of magnetic loops between the star and exoplanet would result in magnetic reconnection events that trigger superflares.They further hypothesized that these flares would be luminous at X-ray to radio wavelengths, in analogy to RS CVn and solar flares.
More recently, attention has focused on the radio emission from M dwarfs that may result from star-planet interactions with potential exoplanet companions, including the M9 dwarf TVLM 513-46546 as observed at 4-5 GHz (Wolszczan & Route 2014;Curiel et al. 2020), the M4.5 dwarf YZ Ceti at 550-900 MHz and 2-4 GHz (Pineda & Villadsen 2023;Trigilio et al. 2023), the M5.5 dwarf Proxima Centauri at 1.1-3.1 GHz (Pérez-Torres et al. 2021), the M4.5 dwarf GJ 1151 at 150 MHz (Vendantham et al. 2020;Blanco-Pozo et al. 2023), and several M1.5-M6.0dwarfs observed at 144 MHz (Callingham et al. 2021).In each case, highly circularly polarized, high brightness temperature emission consistent with electron cyclotron maser (ECM) emission was detected that may result from magnetic starplanet interactions.Alternatively, the observed characteristics at lower frequencies are also consistent with plasma radiation, and several other non-planet-hosting M dwarfs exhibit similar emission (e.g., Stepanov et al. (2001); Lyutikov et al. (2020)) Therefore, radio surveys are a potentially powerful tool to observe coronal flaring and shock processes that may result from magnetic star-planet interactions and diagnose their associated plasma density, local magnetic field, particle acceleration, and wave propagation processes (e.g., Benz & Güdel (2010)).While magnetic star-planet interactions may cause detectable stellar radio emissions, it is important to note that they would do so against a backdrop of a legion of intrinsic stellar radio phenomena.High spatial and temporal resolution observations of the Sun reveal that thermal bremsstrahlung, incoherent gyrosynchrotron, and coherent emissions processes contribute to its myriad radio emissions.At mm to cm wavelengths, gyrosynchrotron dominates, while at dm to m wavelengths, coherent plasma radiation dominates (Bastian et al. 1998).The ECM emission may cause the fine structure within type IV bursts and solar microwave spikes, which range in frequency from ν ∼50 MHz to 10 GHz and ν ∼100 MHz to 5 GHz, respectively (Treumann 2006;Pick & Vilmer 2008).Dulk (1987) summarized how these four radio emission mechanismsthermal bremsstrahlung, gyrosynchrotron, plasma, and ECM -differ in brightness temperature, circular polarization fraction, frequency of operation, and timescale of variability.A search for stellar flaring at 5-10 GHz may be especially promising, since gyrosynchrotron microwave bursts tend to peak within this frequency window, and several types of ECM-generated activity are also present.
Given these considerations, we conducted a survey of nearby exoplanet-hosting systems at Arecibo Observatory (AO) in 2010-2011 to search for enhanced coronal radio emissions that may be linked to magnetic starplanet interactions.A detailed account of the analysis and interpretation of one exoplanet system that we surveyed, HD 189733, was presented in ROME I and II (Route 2019;Route & Looney 2019).Section 2 describes our target selection methodology, instrumentation, and observations.Section 3 presents the detection limits and orbital phase coverage of our survey for potential magnetic star-planet interactions.Section 4 contextualizes these results with respect to searches for magnetic star-planet interactions at other wavelengths and considers the reasons for nondetection.Section 5 concludes by reviewing the significance of this survey and providing suggestions for future work.The results of our survey indicate that the detection and characterization of exoplanet magnetism at radio wavelengths via magnetic star-planet interactions remains an admirable, although difficult, goal.

Target Selection
Target systems were chosen based on the potential for their physical properties to generate magnetic starplanet interactions and their location in the sky (Table 1).Physical properties that served as selection criteria include large masses and small semimajor axes, a. Exoplanet mass may alter the potential for starplanet interactions in somewhat different ways depending upon whether a dipolar or unipolar magnetic interaction with the host star occurs.For a strongly magnetized exoplanet, dipolar interactions would cause nearly continuous reconnection between large-and/or small-scale stellar magnetic fields and the exoplanetary magnetic field.In this case, theoretical models indicate that an exoplanet with a larger mass would have a larger convective core and greater thermal flux to drive a magnetic dynamo that results in a stronger magnetic field (Griessmeier et al. 2007;Christensen et al. 2009;Reiners & Christensen 2010).A stronger magnetic field would yield a larger magnetosphere (obstacle) to interact with a strongly magnetized stellar wind (plasma flow).On the other hand, unipolar interactions between a weakly/ unmagnetized exoplanet and a strongly magnetized stellar wind would create a fieldaligned current circuit propagated via magnetohydrodynamic waves.In this case, a more massive exoplanet presents a larger cross-section (obstacle) via its surface or its ionosphere to interact with a strongly magnetized stellar wind (plasma flow; Zarka 2018).Exoplanets close to their host stars may be tidally locked with weakened or nonexistent intrinsic magnetic fields (Zarka 2007).In both cases, keV electrons would be accelerated toward the source of the strong magnetic field, the magnetized star, with a fraction of the electronic kinetic energy radiated away at radio wavelengths.Each targeted system has at least one exoplanet at a distance a 0.1 au from its host star, which is the range considered by Cuntz & Shkolnik (2002) to be especially promising to enhance stellar activity.Our target stars range from spectral types G0 to M2.5 and include 51 Peg and HD 209458, which both appeared in their study, as well as GJ 176, HD 46375, 55 Cnc, GJ 436, HD 102195, and HD 189733.Physical system properties aside, we targeted systems <50 pc from the Sun.Since stellar coronal flare luminosities follow a power-law distribution (e.g., Güdel et al. 2003), increasingly distant sources require stronger and rarer magnetic star-planet interaction induced flares for detection.The single, fixed dish nature of the Arecibo radio telescope also constrained our targets to have declinations of 0 • to +38 • .

Instrumentation and Observations
The targeted exoplanet systems were surveyed during AO program A2471, conducted from 2010 January 6 to 2011 September 7.Each exoplanetary system was observed for a total of 1-2 hr, sometimes across multiple epochs, although no single epoch was longer than 2 hr (Table 2 1 ).The ∼2 hr observing session limitation represents the maximum time for a target to transit the field of view of the fixed dish of the Arecibo radio telescope.During these observing sessions, 20 s calibration on-off scans using a local oscillator were interlaced between 10 minute target scans.
Although the 305 m dish is fixed, the William E. Gordon radio telescope at AO offered several advantages relative to other facilities, including its exquisite sensitivity, accurate polarization measurements over short time scales, and large, ∼1 GHz bandpass, which was exemplary in 2010-2011.The observing system consisted of the C-band receiver 2 and the recently commissioned Mock spectrometer 3 .The system temperature ranged from 25-32 K, with antenna gains of 5.5-9.0K Jy −1 .The half-power beamwidth varied from 0.97' × 1.09' (azimuth × zenith angle) at 4.500 GHz central frequency to 0.79' × 0.92' at 5.400 GHz (Chris Salter, personal communication).Full Stokes parameters were computed by the Mock spectrometer array based on the dual-linear polarization input signals from the receiver.This array consists of seven field-programmable gate array (FPGA) Fast Fourier Transform (FFT) Mock spectrometers, which were individually tuned to central frequencies of 4. 325, 4.466, 4.608, 4.750, 4.892, 5.034, and 5.176 GHz, resulting in an ∼1 GHz simultaneous bandpass stretching from 4.239 to 5.262 GHz.Each 8192channel spectrometer yields a 172 MHz bandpass that overlaps the frequency range of each neighboring spectrometer by ∼30 MHz.The spectrometers were configured for 0.1 s sampling.
Prior to analysis, data were rebinned to a frequency (time) resolution of 83.6 kHz (0.9 s) to improve the signal to noise ratio characteristics.Radio frequency interference (RFI) was reduced through an iterative statistical process detailed in Route (2013).Data analysis consisted of the creation of dynamic spectra (timefrequency spectrograms) and bandpass-averaged time series graphs of flux density for each Stokes parameter in each spectrometer.Since one of the hallmarks of the two most likely types of radio emission, ECM and gyrosynchrotron, is significant circular polarization, we searched Stokes V dynamic spectra for 10% bursts of emission, which are readily apparent when present (e.g., Figure 1).AO is insensitive to quiescent radio emission due to its local calibration procedure and confusion limitations; thus, only rapid (∼minutes), circularly polarized radio bursts could be detected.Candidate bursts were identified by surpassing a 3σ noise threshold in their Stokes V time series graph.The 1σ instrumental sensitivity was ∼0.1-1.2 mJy, with per-scan sensitivity distribution properties described in Route (2017).Once identified, candidate burst Stokes Q and U components were examined for behavior indicative of RFI, and their morphology was assessed against a library of RFI artifacts acquired during our previous surveys at AO to eliminate false positives (Figure 1).

RESULTS
No radio bursts of any type were detected from any exoplanet-hosting system in our survey.We emphasize that this survey was sensitive only to 10% circularly polarized bursts of ECM or gyrosynchrotron radio emission of several minutes duration and was entirely insensitive to quiescent (unpolarized and/or slowly varying) emission.However, we can compute 3σ detection limits on the radio emission from target stars using the standard deviation (σ) of the frequency-integrated time series of the cleanest Mock spectrometer, centered at 4.466 GHz (Table 3).Comprehensive, frequency-dependent sensitivity limits on a system-by-system basis are available in Route (2013).A comparison of these detection limits to several types of stellar flares is depicted in Figure 2.
Incidentally, since each surveyed planetary system fits entirely within the half-power beamwidth, not only do these detection limits constrain the luminosity of stellar flares induced by magnetic star-planet interactions, but they also constrain the detectability of magnetospheric radio emissions from exoplanets within those systems.Planetary systems with massive exoplanets close to their parent stars are thought to be especially promising targets for the detection of magnetospheric, or auroral, planetary radio emission due to the larger incident Poynting flux and denser stellar wind impinging on their magnetospheres (Zarka 2018).Although this survey failed to detect any such radio emissions, the detection limits in terms of maximum Jovian ECMinduced emissions are also reported in Table 3.

Orbital Phase Coverage of Companions for Star-Planet Interactions
One important consideration in our search for magnetic star-planet interactions is the orbital phase coverage of our observations.Borrowing from the nomenclature for hot Jupiters (e.g., Dawson & Johnson (2018)), we compute the orbital phase coverage for the hot (P <10 days) exoplanet in each system, where T 0 represents the beginning or end of the observation, T t represents the measured (or estimated) midpoint of transit (or periapse), P is the exoplanet orbital period, and E orb is the orbital cycle.Orbital phase 0.0 denotes when transit or periapse would occur.The results for hot exoplanets are shown in Table 4 and Figure  3.
Even among the shortest-period exoplanets, the orbital phase coverage is quite modest since many targets were observed for ∼1 hr, and only four targeted systems were observed for >1.5 hr.The peak orbital phase coverage occurs for 55 Cnc e (∆φ orb =0.08), although the coverage and timing of the observations for HD 189733 are also significant (Section 4.1.2).We can also calculate the orbital phase coverage for the warm (10 days≤ P <200 days) exoplanet 55 Cnc b, which is a tiny ∆φ orb =0.004.Thus, despite the excellent sensitivity of AO, modest orbital phase coverage only loosely constrains the incidence of star-planet interactions within the surveyed systems.No change was detected to a level of 0.001.Similarly, we did not detect any sign of magnetic star-planet interactions in this system, although our observations only occur at φ orb ∼0.68.This orbital phase is somewhat removed from φ orb ∼0.5, which theory indicates may experience augmented activity due to tidal star-planet interactions.
55 Cnc.-The 55 Cnc binary system consists of a K0 primary, an M4.5 secondary, and five exoplanets that circle the primary.The closest orbiting exoplanet, the sub-Neptune 55 Cnc e, may be close enough to its host star to cause star-planet interactions, although its magnetosphere would likely present a very small cross section to the stellar wind for plasma flow-obstacle interactions.While Folsom et al. (2020) determined that 55 Cnc e likely orbits within its parent star's Alfvén radius, Morris et al. (2021) determined that a planetary magnetic field of B ∼1 G would generate ∼10 5 less energy than required to explain the photometric phase variations observed in the system by CHEOPS.Similarly, we found no evidence for star-planet induced coronal emissions in this system.51 Peg.-Saar & Cuntz (2001) found no sign of starplanet interactions in the Ca II infrared triplet (IRT) down to ∼3.2%.Shkolnik et al. (2003) reported significant nightly activity for 51 Peg which they suggested might be cyclical and phase-shifted based on ∼10 nights of Ca II H and K observations.Subsequent observations found no variability down to 0.001 (Shkolnik et al. 2005).We found no evidence for magnetic star-planet induced coronal activity either at φ orb ∼0.22 or at φ orb ∼0.43, near the anticipated tidal activity enhancement zone.
HD 209458.-Aswith 51 Peg b, preliminary monitoring of chromospheric Ca II H and K activity suggested nightly activity which could be periodic in nature, but was later found to lack variability down to the level of 0.001 (Shkolnik et al. 2003(Shkolnik et al. , 2005)).Alternatively, Jensen et al. (2012) found an anomalous shifting of the broad component relative to the narrow component of the Hα line in transmission spectra during transit and secondary eclipse.They suggested that stellar chromospheric activity enhanced by star-planet interactions may alter the Hα line shape, although this phenomenon may also point to a data reduction artifact.Our radio observations do not find evidence of magnetic star-planet induced coronal activity within this system at φ orb ∼0.01 or φ orb ∼0.87.Note that enhanced activity near φ orb =0.0 would be anticipated in both tidal and magnetic interaction cases.

HD 189733
The proximity of the HD 189733 system to the Earth, the brightness of the primary star, and the fortuitous alignment of its hot Jupiter companion that results in a transit every 2.219 days have made the system a particularly inviting target for study.With a chromospheric activity index, S-value=0.525,HD 189733 A is among the 10% most active K dwarfs in the the California and Carnegie Planet Search Project.Transmission spectroscopy of the frequent transits has resulted in the discovery of the exoplanet atmosphere's composition, thermal properties, and evaporation.Multiwavelength observations of the system have revealed a wealth of stellar activity from optical to X-ray wavelengths, including spots and flares at all rotational and exoplanetary orbital phases.Compellingly, several research groups announced the discovery of enhanced chromospheric, transition region, and coronal activity that may be phased with exoplanet orbital phases φ orb ∼0.5-1.0 and caused by magnetic and/or tidal star-planet interactions.In one case, it was furthermore suggested that the primary star may actively accrete evaporated material from HD 189733 b (see ROME I and references therein).
In ROME I, we analyzed AO radio observations to search for enhanced coronal activity associated with such star-planet interactions.We also analyzed Microvariability and Oscillations of Stars (MOST), Automated Photoelectric Telescope (APT), and Wise photometry to search for periodicities suggestive of stellar hemispheres that contain more darkened starspots than average, which may point to enhanced photospheric activity associated with these interactions.We also compiled all observations and reported instances of stellar spotting and flaring at radio through X-ray wavelengths to attempt to discover enhancements phased to the orbital period of the exoplanet.Our results could find no statistically significant enhancement, including at the orbital phases where other research groups thought starplanet interactions may occur.In ROME II, we again leveraged MOST, APT, and Wise photometry to search for photometric bright features that may be related to accretion structures, using the T Tauri paradigm as a template.After computing plasma parameters such as density, temperature, accretion rate, and flare structure sizes, we determined an accretion rate of Ṁ ∼ 10 9 to 10 11 g s −1 , which would result in undetectably low activity at all wavelengths.Since the publication of ROME II, new flare observations and several reanalyses of existing optical and UV data sets have been presented, which we will address in a forthcoming publication (e.g., Bourrier et al. (2020) and references therein).
However, for the purposes of this paper, we focus on research that may directly influence the ability of HD 189733 A or its exoplanet companion to produce radio emission.Based on their convective dynamo model, Reiners & Christensen (2010) estimated that HD 189733 b would have a detectable radio flux density of Φ radio ∼57 mJy up to an ECM cutoff frequency of ν ∼39 MHz.Zaghoo & Collins (2018) used a revised interior conductivity profile to estimate peak ECM emission at ν ∼15 MHz with a radio flux density of Φ radio ∼20 mJy.Kavanagh et al. (2019) leveraged the 3D magnetohydrodynamic modeling code BATS-R-US to investigate the prospects for the detection of radio emissions from HD 189733 b.They used the magnetic field topology maps of HD 189733 A derived from several epochs of Zeeman Doppler Imaging (ZDI) observations from Fares et al. (2017) as inputs.They found that an exoplanetary polar magnetic field of B ∼10 G would be required to create ECM emission at ν ∼ 25 MHz that is higher than the plasma frequency of the stellar wind environment and thus could escape to an observer with Φ radio ∼50-130 mJy.Our radio observations do not shed any light on these predictions, as no exoplanet magnetospheric activity was observed.We also find no evidence for magnetic star-planet induced coronal activity within this system at φ orb ∼0.57-0.61(ROME I).

Potential Reasons for Nondetection of Star-Planet Interactions
There are several potential reasons for the nondetection of coronal radio emissions induced by magnetic star-planet interactions, including sensitivity, observing frequency, observation duration, and viewing geometry.Considering sensitivity first, ECM-induced stellar flaring would likely be detectable in all cases (Figure 2), as was found in our study comparing the activity of HD 189733 A to active M dwarf ECM flares (ROME I).Similar to the findings of that work, only luminous gyrosynchrotron flares might be visible from some targets within our survey, as calculated from the Güdel-Benz relationship and their measured X-ray emissions.It is also possible that the radio counterparts of extremely luminous optical superflares may be visible from our survey (Schaefer et al. 2000).As mentioned in Section 1, it has been hypothesized that such flares may be powered by the magnetic interaction of close-in giant exoplanets.However, it is nontrivial to estimate the radio luminosity of flares based on their optical emissions alone (e.g., Priest 2014).If such flares have radio luminosities comparable to their optical emission, L ∼5×10 31 erg s −1 , they would be readily detectable by our survey.
Figure 2 indicates that a single order-of-magnitude increase in sensitivity should make the stellar flares of the nearest active stars, such as HD 189733 A, accessible.An improvement in sensitivity of ∼10 1 -10 2 would also enable the characterization of coronal activity from older, less active stars such as the Sun.However, given the difficulty of untangling star-planet induced flaring from intrinsic coronal activity (i.e.ROME I), searches for the former phenomenon will require high-cadence, multiepoch observations.
Next we consider the frequency range of our survey with respect to coronal magnetic fields.While the mean solar photospheric magnetic field is B ∼3-10 G, within active regions, this field can reach B ∼100 G, and within sunspot umbrae, field strengths as high as B ∼6 kG have been measured (Priest 2014).Within the solar corona, magnetic fields as high as B ∼2.5 kG have been measured above photospheric active regions.Gyroresonance, gyrosynchrotron, and ECM emissions all contribute to solar coronal radio emissions that range in frequency from ν =1 to 20 GHz (Aschwanden 2005).Among these, only activity caused by the latter two processes would be detectable at great distances.Thus, the frequency range of our survey was well suited for the de-tection of stellar coronal radio emissions associated with both intrinsic and star-planet induced activity.
Companion exoplanets may induce enhanced photospheric, chromospheric, and coronal emissions, including radio emissions, in their host stars either once per orbit due to magnetic effects or twice per orbit due to tidal effects.Thus, the host star rotation, hot exoplanet orbital, and system beat periods are all of interest to untangle the effects of intrinsic stellar activity from activity induced by star-planet interactions.However, hot Jupiter orbital periods are shorter than their host star rotation periods, making the orbital phase coverage the timescale of interest.As described in Section 3.1, the maximum orbital phase coverage for our hot Jupiter targets is ∆φ orb ≤0.08 (Table 4).Thus, to untangle the effects of magnetic star-planet induced activity from intrinsic stellar activity would require nearly complete orbital phase coverage over multiple orbits.As ROME I and ROME II indicated, even for as well-studied a system as HD 189733, it remains unclear whether starplanet interactions occur there.
Finally, we consider the system viewing geometry.Gyrosynchrotron radiation, ubiquitous in solar flares, is generated by mildly relativistic electrons with Lorentz factors of 1 γ 5. Emission is increasingly beamed at higher energies, thereby reducing the opening angles of emission from 90 • to 10 • with respect to the local magnetic field, which can be quite complicated in active regions (Bastian et al. 1998).In situ measurements from the magnetized solar system planets and maser numerical simulations indicate that ECM emission is beamed into a hollow cone with half-apex angles of 30 • to 90 • (Zarka 1998;Treumann 2006), while numerical simulations of ultracool dwarf ECM emission suggest that half-apex angles of ∼1 • may be possible (Lynch et al. 2015).This beaming pattern could critically constrain the detection of such emission from active regions participating in star-planet interactions.However, given that stellar active regions likely span the entire stellar surface or may be restricted to large activity bands (e.g., McIntosh et al. (2014)) visible from many vantage points, it is unlikely that viewing geometry plays a major role in reducing star-planet interaction detectability.On the other hand, the beaming of ECM emission may contribute to the nondetection of magnetospheric radio emissions from the exoplanet companions in the surveyed systems.

CONCLUSION
In this third installment in the Radio Observations of Magnetized Exoplanets (ROME) series, we step back from analyzing the details of a particularly interesting exoplanetary system that may engage in star-planet interactions, HD 189733 A/B/b (ROME I and II), and present the complete results of our 2010-2011 survey that sought to detect and characterize stellar radio emissions that may be caused by magnetic star-planet interactions with close-in exoplanets.These interactions may induce enhanced coronal, transition region, chromospheric, and photospheric activity in their host stars.Such activity may generate powerful ECM or gyrosynchrotron flares phased to the exoplanet orbit.
The survey targeted eight systems that have exoplanet companions within a 0.1 au of the host star.This survey leveraged the exquisitely sensitive 305 m Arecibo radio telescope tuned to 4.2-5.2GHz to search for strong stellar magnetic fields and probe the optimal frequency range for the detection of ECM and gyrosynchrotron solar radio flares.Dynamic spectra in all four Stokes parameters from the Mock spectrometers were used to search for 10% circularly polarized ECM and gyrosynchrotron emission and to mitigate RFI.All targets were observed for 2 hrs, the maximum time for them to transit the fixed dish.
Our survey did not detect stellar radio flares induced by magnetic star-planet interactions.The 3σ sensitivity for targeted systems ranged from 0.98 to 2.1 mJy, resulting in detection limits of ν L ν = 5.1×10 23 erg s −1 to 1.4×10 25 erg s −1 .These luminosity detection thresholds correspond to modest ECM flares, such as those that appear on the M3.5 AD Leo, and powerful gyrosynchrotron flares, such as from the young active star ǫ Eri.It has been hypothesized that extremely luminous optical superflares from FGK stars might be caused by starplanet interactions, which have luminosities L ∼5×10 31 erg s −1 .Although our survey was sensitive to radio emissions of energy output equal to these optical superflares, it is unclear whether the radio counterparts to such flares would be detectable with current instrumentation.Incidentally, these luminosity detection thresholds also constrain the high-frequency magnetospheric (auroral) emissions that may come from the exoplanets themselves.
Several targets, including HD 46375, 55 Cnc, 51 Peg, HD 209458, and HD 189733, have previously been investigated for signs of chromospheric and photospheric activity enhancement due to potential tidal or magnetic star-planet interactions.Yet in each case, subsequent work challenged the initial putative discoveries of increased magnetic activity phased with the orbital or beat periods of their hot Jupiter companions.Similarly, our search for enhanced coronal activity did not detect signs of magnetic star-planet interactions.The inability to detect intrinsic stellar flaring, much less enhanced coronal activity from star-planet interactions at radio wavelengths, indicates that future efforts will require an ∼10 1 -10 2 increase in sensitivity.
These survey results provide guidance as to how to enhance the scientific return of future surveys for starplanet interactions at radio wavelengths.First, increases in sensitivity of 1-2 orders of magnitude are required to detect magnetic star-planet interactions.Second, future surveys should simultaneously observe as wide a bandpass as possible given the great uncertainty in flare magnetic field strengths that may be generated by starplanet interactions.Third, future surveys should include polarimetry to distinguish among and characterize various types of astrophysical emission, as well as an RFI mitigation strategy.Fourth, multiepoch, highcadence observations are required to untangle the effects of intrinsic stellar activity on differentially rotating stars from enhanced activity induced by tidal and magnetic star-planet interactions.Finally, high-fidelity theoretical simulations are required to model the influence of star-planet interactions on stellar atmospheres and assess the prospects for observing enhanced photospheric, chromospheric, transition region, and coronal activity within hot exoplanet systems.16) Martins et al. (2015).
a Also known as HD 285968 b.
Figure 1.Stokes V dynamic spectra and bandpass-integrated time series from our exoplanet (left) and ultracool dwarf (right) radio surveys at AO. Bandpass corrections, which consist of applying a polynomial fit from the central ∼80% of each spectral window to the ∼10% on each edge and renormalizing the flux density, mitigate subband discontinuities (scalloping) in sensitivity.Horizontal gray bars near 4.8 and 5.15 GHz represent the excision of strong RFI.Even though a patch of left circularly polarized emission occurs at t ∼425 s in a dynamic spectrum from GJ 176 at left, follow-up analysis revealed that this feature is also strongly linearly polarized and appears in multiple data sets of different objects.The emission is therefore an RFI artifact.At right is a real, ∼70% circularly polarized radio flare from the T6.5 ultracool dwarf 2MASS J10475385+2124234 (Route & Wolszczan 2012, 2013).2019).While these emissions are generated by incoherent gyrosynchrotron flaring, flares observed at AD Leo are created by the ECM process and are more readily detectable due to their coherence and beaming (Stepanov et al. 2001).We establish a maximum luminosity for the "superflares" category using the 15 day B-and V -band flare observations of the F8 IV star 5 Serpentis, although it is nontrivial to relate the optical properties of the flares to those at radio wavelengths (Schaefer et al. 2000;Priest 2014).The dashed and dotted lines mark the detection thresholds of our most and least sensitive observations, those at GJ 176 and HD 209458, respectively, relative to these flare energies.
Exoplanet Systems of Special Interest 4.1.1.Targets with Previous Magnetic Star-Planet Interaction Observations HD 46375.-Shkolniket al. (2005) searched for periodic star-planet induced chromospheric variability in the Ca II H and K lines (3968 and 3933 Å) among HD 46375 and several other program stars.

Figure 2 .
Figure 2. Survey νLν radio luminosity detection limits compared to known stellar flare luminosities.For each object, the spectral type(s) are given in parentheses.The flare radio luminosity ranges for HD 189733 A, ǫ Eri, and the interacting binary BY Dra AB are calculated from their X-ray flare luminosity ranges via the Güdel-Benz relationship (Benz & Güdel 2010; Route 2019).While these emissions are generated by incoherent gyrosynchrotron flaring, flares observed at AD Leo are created by the ECM process and are more readily detectable due to their coherence and beaming(Stepanov et al. 2001).We establish a maximum luminosity for the "superflares" category using the 15 day B-and V -band flare observations of the F8 IV star 5 Serpentis, although it is nontrivial to relate the optical properties of the flares to those at radio wavelengths(Schaefer et al. 2000;Priest 2014).The dashed and dotted lines mark the detection thresholds of our most and least sensitive observations, those at GJ 176 and HD 209458, respectively, relative to these flare energies.

Figure 3 .
Figure 3. Orbital phase coverage of hot exoplanets (P ≤10 d) to search for magnetic star-planet interactions.Maximum orbital phase coverage occurs for 55 Cnc e, (∆φ orb =0.08).GJ 176 b, HD 209458 b, and 51 Peg b each have two observing epochs depicted.The orbital phase coverage of HD 189733 b occurs during phases where previous investigators suggested that they found increased stellar activity induced by potential magnetic star-planet interactions (ROME I) and is therefore significant.

Table 1 .
Survey Target PropertiesNote-Exoplanet-hosting system properties.In the property references column (rightmost column), bold, underlined, and italicized numerals denote discovery, semimajor axis, and orbital period references, respectively.The final number in the column in normal font provides the companion object's mass reference.All distances are from the Gaia Collaboration (2021).

Table 2 .
Observations ListThe scan range column denotes the continuous sequence of on-source, calibration-on, and calibration-off scans that focused on the listed target on a given day.For example, file a2471.20101218.b0s1g0.12500.fitsprovides the on-source data to create the dynamic spectra from the first, lowest-frequency Mock spectrometer (b0) during the GJ 176 observation (scan 12500) that occurred on 2010 December 19 (Figure1, left).Asterisks in the scan range column denote an observing session that immediately continued into the following day.

Table 3 .
(Zarka et al. 2004esultsNote-Detection limits and radio luminosities of target systems.aLuminosity in terms of the average power output of Jupiter's ECM-induced decametric (DAM) during maximum solar activity (as opposed to average power or peak power), or L J,rad = 8.2×10 17 erg s −1(Zarka et al. 2004).
(Route & Wolszczan 2013ensity detection limit presented inRoute (2019), with luminosity updated due to revised distance to HD 189733 (Gaia Collaboration 2021).cRadioemissionflux density detection limits for all HR 8799 exoplanets as reported in our first ultracool dwarf radio survey(Route & Wolszczan 2013).