Space environment and magnetospheric Poynting fluxes of the exoplanet $\tau$ Bo\"otis b

Context: The first tentative detection of a magnetic field on the Hot Jupiter type exoplanet $\tau$ Bo\"otis b was recently reported by Turner et al. (2021). The magnetic field was inferred from observations of circularly-polarized radio emission obtained with the LOFAR telescopes. The radio emission is possibly a consequence of the interaction of the surrounding stellar wind with the planet's magnetic field. Methods: We perform magnetohydrodynamic simulations of the space environment around $\tau$ Bo\"otis b and its interaction with the stellar wind using the PLUTO code. We study the magnetospheric energy fluxes and effects of different magnetic field orientations in order to understand the physical processes which cause energy fluxes leading to the observed radio emission given the proposed magnetic field strength in Turner et al. (2021). Furthermore we study the effect of stellar wind density and pressure on magnetospheric energy fluxes given the uncertainty of extrasolar stellar wind predictions. Results: We find in our simulations that the interaction is most likely super--Alfv\'enic and energy fluxes generated by the stellar wind--planet interaction are consistent with the observed radio powers. Magnetospheric Poynting fluxes are of the order of 1-8 $\times10^{18}$ W for open, semi-open and closed magnetospheres. The Poynting fluxes are energetically consistent with the radio powers in Turner et al. (2021) for a magnetospheric Poynting flux-to-radio efficiency $>10^{-3}$ when the magnetic fields of the planet and star are aligned. In case of lower efficiency factors the magnetospheric emission scenario is according to the parameter space modeled in this study not powerful enough. In case of a magnetic polarity reversal of the host star towards an anti-aligned field configuration, expected radio powers in the magnetospheric emission scenario fall below the observable threshold.


Introduction
Recently, tentative measurements of auroral radio emission from the hot Jupiter exoplanet τ Boötis b were obtained with the Low Frequency Array (LOFAR; Turner et al. 2021).These observations might be considered the strongest evidence so far of an intrinsic magnetic field on a planet outside the Solar System if the emission indeed originates from the planet's vicinity.They then also imply that τ Boötis b possesses a magnetosphere (MS) that interacts with its surrounding stellar wind.The radio observations by Turner et al. (2021), if confirmed, thus help pave the way for the field of extrasolar space physics.In this work we therefore use properties from the observed radio signals to derive new constraints on the space environment around τ Boötis b.
The massive hot Jupiter τ Boötis b (Butler et al. 1997) is a very good candidate for remotely observing a powerful interaction of a stellar wind with an exoplanet's MS for several reasons: At ∼16 pc, the τ Boötis system is relatively close to the Solar System.The planet orbits its host star, τ Boötis A, at a short distance of 0.046 astronomical units (Butler et al. 1997).Additionally, its large mass (>5 M Jupiter ) may cause its exobase to remain close to the planet, leading to a MS not completely filled with dense plasma and thus allowing for radio emission to be produced efficiently and to escape the planet's vicinity (Weber et al. 2018;Daley-Yates & Stevens 2018).
The tentative radio measurements obtained with LOFAR comprise two signals that probably originate from the vicinity of τ Boötis b (Turner et al. 2021).The circularly polarized signals were detected in the 21-30 MHz and 15-21 MHz frequency bands, respectively.The emission possibly originates from gyrating, energetic electrons precipitating toward the planetary polar regions, emitting radio waves generated through the electron cyclotron maser instability (ECMI), which is expected to be the dominant mechanism for exoplanetary radio emission (Zarka 1998;Treumann 2006).From these signals, the planetary magnetic field strength can be inferred directly since the emission frequency corresponds to the local electron gyrofrequency.The existing observations are consistent with expectations for emitted power from the radio-magnetic Bode's law (Zarka et al. 2001(Zarka et al. , 2018;;Zarka 2007), for the polarization (e.g., circular polarization; Zarka 1998;Grießmeier et al. 2005), and for the frequency (i.e., slightly above Earth's ionospheric cutoff; Grießmeier et al. 2007bGrießmeier et al. , 2011;;Griessmeier 2017).The measured radio signal, however, needs multisite follow-up observations, preferably at various radio wavelengths, to confirm and to further constrain the magnetic field environment of τ Boötis b (Turner et al. 2021).
In addition to radio emission (e.g., Grießmeier 2015;Farrell et al. 1999;Zarka et al. 2001;Zarka 2007), other indirect methods have been proposed to detect and constrain the magnetic fields of exoplanets.They include optical signatures in the stellar chromosphere by means of Ca II H&K line excess emission induced by star-planet interactions (SPIs; Cuntz et al. 2000;Cauley et al. 2019;Shkolnik et al. 2003Shkolnik et al. , 2005Shkolnik et al. , 2008) ) and asymmetries in near-UV stellar light curves together with UV absorption signatures caused by transiting planetary bow shocks (Vidotto et al. 2010(Vidotto et al. , 2011;;Llama et al. 2011).The SPIand transit-related observations can lead to false positives (e.g., Turner et al. 2021Turner et al. , 2016;;Kislyakova et al. 2016;Preusse et al. 2006;Lai et al. 2010;Kopp et al. 2011;Miller et al. 2012Miller et al. , 2015;;Bisikalo et al. 2013;Alexander et al. 2016;Gurumath et al. 2018;Route 2019) due to the sets of model assumptions involved in the process.Radio observations, in contrast, can directly constrain the magnetic field amplitude and are therefore less susceptible to false positives (Grießmeier 2015).The success of radio observations has been demonstrated in the past in the Solar System.For example, Jupiter's magnetic field was discovered through radio observations (Franklin & Burke 1958) before spacecraft confirmed it with in situ magnetometer measurements.
Since τ Boötis b may be the first exoplanet with a directly observed magnetic field, it provides a unique opportunity to constrain the space environment around this exoplanet.However, various properties of τ Boötis b are unknowns, such as the radius, size, and extent of its atmosphere above the 1 bar level, as well as stellar wind parameters.τ Boötis A is a solar-like F7 IV-V star (Gray et al. 2001); its coronal temperature and pressure might therefore be comparable to those of the Sun.The coronal base density, and consequently the stellar wind mass loss rate, is the most uncertain free parameter of previous studies of the stellar wind from τ Boötis A (Vidotto et al. 2012;Nicholson et al. 2016).Recently, new constraints on stellar winds of M dwarf stars were determined via the use of astrospherical absorption signatures induced by the interaction of the stellar wind with the interstellar medium (Wood et al. 2021).The question naturally arises of whether stellar wind-planet interactions might also produce observable signatures capable of providing constraints on stellar wind properties, such as density (i.e., the mass loss rate) or pressure (i.e., temperature), which will be addressed in this paper.
The proximity of hot Jupiters to their host stars can potentially cause sub-Alfvénic SPIs, which is expected to produce observable signatures in the stellar atmosphere (e.g., chromospheric emission Ca II H & K line emission; Shkolnik et al. 2003Shkolnik et al. , 2008;;Cauley et al. 2019) or the planetary atmosphere (e.g., auroral radio emission; Cohen et al. 2018;Turnpenney et al. 2018;Bastian et al. 2022;Kavanagh et al. 2021Kavanagh et al. , 2022)).Such magnetic SPIs in exoplanetary systems were excessively studied by means of magnetohydrodynamic (MHD) simulations (e.g., Preusse et al. 2006Preusse et al. , 2007;;Zhilkin & Bisikalo 2020;Varela et al. 2018Varela et al. , 2022)), partly with a focus on the far-field interaction and incorporating self-consistent stellar wind models (e.g., Strugarek et al. 2014Strugarek et al. , 2019a,b;,b;Cohen et al. 2011Cohen et al. , 2014;;Vidotto et al. 2015;Vidotto & Donati 2017).Very little modeling of super-Alfvénic stellar wind-MS interactions has been done to our knowledge despite the fact that a large number of close-in exoplanets might be exposed, at least temporarily, to super-Alfvénic wind conditions (i.e., orbiting outside the Alfvén surface; Zhilkin & Bisikalo 2019).In the case of super-Alfvénic or, more precisely, super-fast magnetosonic stellar wind flows, a bow shock forms upstream of the planetary MS because of the flow being faster than the fastest MHD wave mode.In this case, the planet does not interact with the star since waves are not able to propagate upstream.This might be the case for τ Boötis b, as the planet is likely exposed to a super-fast stellar wind, according to Nicholson et al. (2016).
The generation of radio emission from exoplanets, as well as its properties and dependence on stellar wind and planetary parameters, was studied intensively using numerical simulations (Nichols & Milan 2016;Varela et al. 2016Varela et al. , 2018;;Turnpenney et al. 2020;Daley-Yates & Stevens 2018;Kavanagh et al. 2020) for other or generic exoplanets.However, little to no emphasis was placed on studying the detailed spatial structure and energetics of the magnetospheric Poynting fluxes that ultimately deliver the available electromagnetic energy capable of driving planetary auroral emissions at various wavelengths.
In order to better understand the space environment around τ Boötis b, we performed MHD simulations of the near space environment of τ Boötis b and its magnetic field interacting with the surrounding stellar wind plasma using the PLUTO code.The stellar wind model is based on wind simulations (Vidotto et al. 2012;Nicholson et al. 2016) driven by magnetic surface maps derived from magnetic measurements of τ Boötis A (Marsden et al. 2014;Mengel et al. 2016;Jeffers et al. 2018).The magnetic field estimate of the planet's intrinsic field, based on the tentative magnetic field strengths derived by Turner et al. (2021), is used to model the planetary MS.We specifically aim to better understand the magnetospheric energy fluxes around τ Boötis b and, more generally, hot-Jupiter-type exoplanets that are exposed to similar stellar wind conditions.We also address the question of how stellar wind variability in the time-independent case affects magnetospheric Poynting fluxes and therefore possible radio powers generated by the interaction.
The paper is structured in the following way: An overview of the physical model that describes the plasma interaction of τ Boötis b with the surrounding stellar wind is given in Sect. 2. The numerical setup is summarized in Sect.2.1, and details about the stellar wind model can be found in Sect.2.2.The τ Boötis b model is described in Sect.2.3.In the subsequent section we show our results, starting with a general description of the interaction in Sect.3.1, followed by a study of the spatial structure of Poynting fluxes in Sect.3.2.1.Then we discuss the energetics of the interaction in Sect.3.2.2,where we also compare possible radio emission output with the observations by Turner et al. (2021).The results are followed by a discussion about the role and importance of the stellar wind in powering the energy fluxes in the MS of the exoplanet in Sect.4.1.At last, we discuss possible auroral radio emission and its detectability in the scope of stellar wind variability (Sect.4.2).

Numerical simulation
In this section, we introduce our physical model and the numerics to describe the interaction of τ Boötis b and its intrinsic magnetic field with its surrounding stellar wind.The MHD model together with the numerical model and coordinate A133, page 2 of 19 F. Elekes and J. Saur: A&A proofs, manuscript no.aa44947-22 system are presented in Sect.2.1.We introduce the stellar wind that is included as boundary condition for the plasma variables in Sect.2.2 followed by the description of parametrizations of physical processes introduced by the planet and its atmosphere in Sect.2.3.

Method
We performed single-fluid ideal, non-resistive and nonviscous MHD simulations using the open-source code PLUTO (v. 4.4) in spherical coordinates (Mignone et al. 2007).The MHD equations to solve are where ρu is the momentum density, u the velocity, ρ the mass density, p t the total pressure (e.g., magnetic and thermal) and p the thermal pressure.B is the magnetic flux density, −u × B in Eq. ( 4) is the electric field in the ideal limit with infinite electrical conductivity.E t is the total energy density, E t = ρe + ρv 2 /2 + B 2 /2µ 0 , and e the specific internal energy.The system is closed by the equation of state in the form p = ρe(γ − 1), where γ is the ratio of specific heats for the adiabatic case.
As for magnetic diffusion we do not include a diffusion term in the induction equation (Eq.( 4)) but point out that numerical diffusion, especially for coarse grids such as in our simulation, introduce numerical diffusion sufficient to allow for reconnection (see Varela et al. 2018).To justify this assumption we performed test simulations incorporating magnetic diffusion and found it to not influence the results of this paper significantly (see Appendix B for a detailed discussion on this topic).
We include plasma production, P, and loss terms, L (Eqs.( 1)-( 3)) to account for photoionization, dissociative recombination together with associated momentum and internal energy transfer between neutral atmospheric and magnetospheric plasma particles as well as ion-neutral collisions.We note that the neutral species is not simulated or altered by the interaction with the ion species.Details on how plasma production and loss are modeled can be found in Sect.2.3.The mass of plasma particles is denoted by m p and m n describes the mass of neutral particles.We assume the plasma to completely consist of ionized hydrogen atoms, m p = m H + .The atmosphere only consists of neutral molecular hydrogen, m n = m H 2 .
The spherical grid consists of 256 non-equidistant radial grid cells as well as 64 and 128 equally spaced angular grid cells in the θ and ϕ dimension, respectively.The radial grid is divided into three regions.From 1 to 1.2 planetary radii (R p ), the grid contains 10 uniform cells.After that from 1.2 to 12 R p the next 150 cells increase in size with a factor of ∼ 1.01 per cell.The last 96 cells from 12 R p toward the outer boundary at 70 R p increase gradually with a factor of ∼1.02.The positive x-axis points parallel to the relative velocity v 0 of the stellar wind in the frame of the planet.The stellar wind magnetic field is assumed to be perpendicular to v 0 and is antiparallel to the z-axis.The y-axis completes the right handed coordinate system.The colatitude, θ, is measured from the positive z-axis, and longitudes, Φ, are measured from the positive y-axis within the xy-plane.The origin is located at the planetary center.We ran all simulations for approximately 3.6 h physical time until a quasi-steady-state is reached in the vicinity of the planet (r < 30).Small fluctuations cannot be avoided, although larger-scale structure and dynamics within the MS remain already almost constant after approximately 2 h physical time.

Stellar wind model
The derived stellar wind parameters from Nicholson et al. (2016) resemble those of the Sun, such as the polytropic index, γ = 1.1 (Van Doorsselaere et al. 2011), and the stellar coronal base temperature, which is not well constrained by observations, is set to 2 × 10 6 K as typical value for the solar coronae (Nicholson et al. 2016;Vidotto et al. 2012;Van Doorsselaere et al. 2011).The magnetic field of τ Boötis A was studied excessively during several epochs and magnetic surface maps as well as several magnetic polarity reversals were observed (Donati et al. 2008;Fares et al. 2009Fares et al. , 2013)).The coronal base density remains an educated guess based on a comparison of emission measure values obtained from X-ray spectra of τ Boötis A (Vidotto et al. 2012;Maggio et al. 2011).Due to the uncertainty of the base density estimate, different stellar wind densities will be investigated separately in the scope of magnetospheric Poynting fluxes and possible radio powers in Sects.3.2.3 and 4.1.
The stellar wind is applied through constant inflow boundary conditions at the upstream hemisphere (Φ = 0 to 180 • ).The magnetic field is assumed to be perpendicular to the relative velocity v 0 of the wind (i.e., parallel to the negative z-axis).The inflow velocity of the plasma, which we call the relative velocity v 0 , is parallel to the x-axis and is composed of the radial velocity of the wind v sw and the orbital velocity of the planet.The adopted plasma parameters of the wind are summarized in Table 1 and were averaged over the several epochs studied by Nicholson et al. (2016).

τ Boötis b model
We assumed a radially symmetric neutral atmosphere with a scale height of H = 4373 km.Thus, the scale height extends over three radial grid cells and, consequently, the neutral atmosphere is sufficiently resolved within the numerical grid.We assume an atmosphere consisting of molecular hydrogen as it is, followed by helium, the most abundant constituent of the Jovian atmosphere (Atreya et al. 2003).The collisional cross section is assumed to be σ in = 2 × 10 −19 m 2 for H + -H 2 collisions with momentum transfer for low-eV relative velocities between the colliding particles (Tabata & Shirai 2000).In our simulations the collision frequency is ν in ≈ 0.5 s −1 , so that ν in = vσ in n(r), where v ≈ v 0 denotes a typical velocity in the system and n n (r) is the atmosphere number density as a function of radial distance from the center, where n n,0 = 8 × 10 12 m −3 is the surface number density.Based on test studies, we found that for n n,0 ≈ 8 × 10 12 m −3 the ionneutral collisions nearly completely bring the incoming plasma flow to a halt in the atmosphere.This results in plasma pile up in form of a shell around the planet.Increasing the density would thus not produce a larger interaction.We used a simplified description of photoionization.We neglected the shadow zone exerted by the planet's body and parameterized plasma production through photoionization using only the radial dependence of the neutral atmosphere density, The radial symmetric ionization partially mimics some night side ionization through electron impact ionization.For the photoionization frequency of hydrogen exposed to a solar-like UV radiation environment at a distance of approximately 0.046 AU from the star we take the value from Kislyakova et al. (2014), ν ion = 6 × 10 −5 s −1 .Plasma loss was introduced through the recombination of hydrogen ions.The loss term therefore depends on the plasma density, Plasma loss is switched off if the plasma density falls below the background density (i.e., n(r, t) ≤ n sw ) as stellar wind ions and electrons recombine significantly slower due to the higher electron temperatures in the stellar wind.Given an electron temperature of roughly T e ≈ 7500 K for a hot Jupiter exoplanet's ionosphere with semimajor axis of 0.046 AU around a Sun-like star derived by Koskinen et al. (2010) and using the formula from Storey & Hummer (1995), we find the hydrogen ion recombination rate, α, to be 5.1 × 10 −19 m 3 s −1 .Further discussion about the underlying assumption about our atmosphere model can be found in Appendix D.
Recent tentative auroral radio measurements from τ Boötis b give a first observational constraint on its magnetic field strength.Turner et al. (2021) found the polar surface magnetic flux density B p to lie between 7.5 and 10.7 G for two right-handed circularly polarized signals.We assume a dipole field and adopt the average value of both Stokes V + signals (Turner et al. 2021), B p = 9.1 G, for our simulations.Furthermore, we study the effect of dipole orientation on the stellar wind-planet interaction through simulating an open (0 • tilt), semi-open (90 • tilt) and closed MS (180 • tilt), where the tilt is measured with respect to the negative z-axis.The various tilts are realized by rotating the stellar background magnetic field accordingly so that the planetary dipole axis is always parallel to the z-axis.Given the strong magnetic variability of τ Boötis A (e.g., several magnetic polarity reversals were observed as well as a chromospheric activity cycle in terms of S-indices of roughly 240 days Donati et al. 2008;Fares et al. 2009Fares et al. , 2013;;Mengel et al. 2016;Mittag et al. 2017;Jeffers et al. 2018) we are also able to study the effect of the host star's magnetic field topology on the stellar wind-planet interaction and associated magnetospheric energy fluxes.
The magnetic field is implemented using the insulatingboundary method by Duling et al. (2014).This method ensures that no radial electric currents exist within the insulating boundary, which we assume to be the planet's neutral atmosphere below its ionosphere.

Results
In this section, we first present results of our modeling, which provides an overview of the space plasma environment of τ Boötis b (Sect.3.1).Then we study in detail the magnetospheric Poynting fluxes in Sect.3.2.

Structure of the interaction
The simulated plasma velocities and pressures according to the basic model (Table 1) are displayed in Fig. 1  case.The magnetic field tilt θ B is the angle between the external field (parallel to the z-axis) and the planet's magnetic moment.We note that, due to the symmetries chosen in our model studies, the stellar wind and intrinsic magnetic field are not inclined with respect to the z-axis; therefore, we also show projected field lines (black solid lines) in the xz-plane.Color contours denote plasma pressure in µPa (right colorbar).Arrows represent velocity components, their magnitudes are color coded (left color bar).The length of arrows indicate the magnitudes of the shown components.Spatial dimensions are given in units of planetary radii.
The intrinsic magnetic field and its corresponding MS poses an obstacle to the stellar wind flow coming from negative x -direction.The flow outside the MS is super-Alfvénic (M A = 5.36) and super-fast magnetosonic (M f = 1.6), (see Table 1), where , with the sound speed c s = γp sw /ρ sw , polytropic index γ = 1.1 (Nicholson et al. 2016) and Alfvén velocity v A = B sw / √ µ 0 ρ sw .The super-fast interaction enforces a bow shock to be formed roughly 5 R p in front of the planet followed by a fairly thick magnetosheath.Since no wave is able to propagate upstream, the stellar wind plasma is unperturbed until the bow shock.The structure of the MS strongly depends on the internal field orientation as visible in Fig. 1 with an increase of overall MS size toward higher magnetic axis tilts.For the open and semi-open MS (Fig. 1 top and middle plot, respectively) two magnetic lobes form, separated by a thin plasma sheet, where open magnetic field lines connect to the stellar wind field several planetary radii downstream.The day side magnetopause, defined by the location of the last closed field line, lies between 3 and 3.5 R p , while the night side magnetopause is located at roughly 5 R p for the open and semi-open MS, respectively.The downstream side magnetopause is very narrow in the z-direction as expected due to the magnetic field lines convected downstream together with the stellar wind flow and due to the magnetic stresses stretching the magnetic field.The closed MS case (Fig. 1, bottom plot) has a night side magnetopause lying several planetary radii (∼17R p ) downstream (not shown in the plots).While the upstream magnetopause is controlled by the stellar wind thermal and magnetic pressure balanced with those exerted by the planet's surroundings, the downstream MS is influenced by reconnection (i.e., the merging of planetary with stellar wind field lines).Magnetic reconnection is most efficient for a magnetic moment parallel to the ambient field (here the z-axis), and therefore the fraction of open planetary field lines connected to the star decreases significantly with an intrinsic field moment being directed antiparallel to the stellar field.As the stellar wind plasma primarily penetrates the MS along magnetic field lines, the amount of plasma and thermal pressure decreases as well with increasing magnetic axis tilt.
We note that, as seen in Fig. 1 (bottom), the MS is completely closed.This is due to the perfect anti-parallel alignment of the planetary and stellar wind magnetic field.
Within the MS, the flow velocity is strongly reduced and has weak upstream components in the negative x-direction due to magnetic tension exerted on planetary field lines.Magnetic reconnection takes place at the upstream and downstream side where velocities, both within and outside the MS, are strongly enhanced due to acceleration through released magnetic energy.Velocities are slightly larger at the flanks of the MS compared to the upstream side and exceed the initial stellar wind velocity at the downstream side where stellar wind as well as planetary field lines merge together again and accelerate the plasma.
Thermal pressures are strongly enhanced within the magnetosheath, where stellar wind plasma is decelerated abruptly and compressed, so that kinetic energy is converted into heat.Plasma may penetrate the MS along open magnetic field lines in the polar cusps where pressure is enhanced as well.The cusps act as channels for plasma transport into the MS.There is a trend toward lower pressures in the cusps for increasing magnetic axis tilt.This is directly connected to the amount of stellar wind plasma advected toward the planet as the amount of A133, page 5 of 19 injected plasma is related to the ability of magnetic field lines to merge with the ambient field.This becomes increasingly difficult for planetary magnetic moments that have components antiparallel to the ambient field; therefore, the area fraction of open magnetic field lines and thus the size of the plasma injection channel is maximal for a completely open MS.Here pressures up to 160 µPa can be reached while the closed MS case shows pressures up to roughly 90 µPa.

Poynting fluxes and aurorae
We are interested in understanding the electromagnetic coupling of the stellar wind with the MS of the exoplanet, its atmosphere and ionosphere.The energy fluxes associated with these electromagnetic coupling processes provide the energetics for the auroral emission from the exoplanet's MS at radio and other wavelengths.Therefore, we study the Poynting flux to calculate the maximum available electromagnetic energy fluxes.We describe the spatial structure of magnetospheric Poynting fluxes in Sect.3.2.1.Then we study the energetics of the interaction and effects of stellar wind variability on magnetospheric energetics in the subsequent Sects.3.2.2 and 3.2.3.The Poynting vector, S, can be rewritten in the ideal MHD case using the convective electrical field (e.g., Saur et al. 2013), A133, page 6 of 19 F. Elekes and J. Saur: A&A proofs, manuscript no.aa44947-22 which is bodily carried by the plasma flow perpendicular to magnetic field lines, denoted by the perpendicular velocity u ⊥ .The Poynting flux describes the transport of magnetic enthalpy, which is a factor of two larger than the magnetic energy density, B 2 /2µ 0 (e.g., Saur et al. 2013).In the remainder of this work, we mostly present Poynting fluxes, but need to consider the factor of two when we compare magnetic energy densities with thermal (i.e., internal energy) or kinetic energy densities based on their flows.For both, the open and closed MS case, flow velocities are strongly reduced at the upstream (ϕ = 0-180 degrees) and downstream (ϕ = 180-360 degrees) side down to speeds below 10 km s −1 .This occurs due to interaction with the intrinsic magnetic field and momentum transfer with the neutral atmosphere.The OCFB is located at roughly θ ≈ 30 • and 150 • in the open MS case.Due to the perfectly anti-parallel configuration of the stellar wind and planetary magnetic field no open planetary field lines exist in the closed MS model.This has also been observed in sub-Alfvénic simulations using this field topology (Ip et al. 2004;Strugarek et al. 2015).
Open magnetosphere model (Fig. 2  Closed magnetosphere model (Fig. 2 right).In the closed MS model, highest velocities can be found near the planetary poles confined to an area below 25 degrees colatitude and similar in the south.These high velocities are caused by tension on high latitude closed field lines that are strongly stretched toward the downstream side by the stellar wind and reach up to 17 planetary radii.Poynting fluxes oriented away from the planet are confined to narrow bands encircling the high latitude polar regions between 40 and 80 degrees colatitude and similar in the south.At the upstream side Poynting fluxes vanish near the equatorial regions due to plasma flow being mainly aligned with planetary field lines.Inward oriented Poynting fluxes occur near the polar axis slightly shifted toward the downstream side.
We now study the Poynting fluxes parallel to the unperturbed background magnetic field because in the Solar System MSs they are considered the root energy fluxes from which a small fraction can be converted into auroral radio emission.Poynting fluxes provide the energy from which wave-particle interaction can draw energy to accelerate electrons (e.g., for Jupiter Hill 2001;Saur et al. 2021).The resulting energetic electrons then can be subject to the electron maser instability (Treumann 2006;Zarka 2007).The interaction of the stellar wind with τ Boötis b's magnetic field perturbs the magnetic and electric field, which causes the Poynting fluxes.To quantitatively assess the associated Poynting flux, we therefore use the unperturbed magnetic background field B 0 = B p,0 + B sw (i.e., the initial dipole and stellar wind field) to calculate the Poynting flux on this field, S • e B 0 .The unit vector e B 0 points in the direction of unperturbed magnetic field lines.These projections give insight on where electromagnetic energy is transported through either propagating magnetic disturbances (i.e., Alfvén waves) or convection.The bottom row of Fig. 2

Energetics of the interaction
To estimate the total available Poynting flux, which serves as the root energy flux, we assume for simplicity that the radio emission is generated in a shell 1 R p above the surface of the exoplanet.This particular choice is inspired by the fact that radio emission around Jupiter and other Solar System planets arises from altitudes about 1R p (or larger) above the planet's surface (e.g., Zarka 1998;Hess & Zarka 2011) where strong electron acceleration takes place (e.g., for Jupiter Mauk et al. 2020).Poynting fluxes within the MS of τ Boötis b only vary little as a function of distance from the planet (see Appendix A for a discussion on the choice of r).
Available electromagnetic power for possible conversion into electron acceleration and radio emission is given by the divergence of the Poynting flux in this shell with volume V, where A shell is the surface area of the shell and n the surface normal vector.To investigate the maximal Poynting flux that can be dissipated in the shell we assume that the Poynting flux entering Notes.Integrated Poynting fluxes over a sphere with radius r = 2R p . (a) Magnetospheric Poynting flux (Eq.( 11)). (b) Magnetospheric Poynting flux parallel to B 0 (Eq.( 12)).
the shell from above or below is dissipated within the shell.For mathematical simplicity, we further let the thickness of the shell grow infinitesimally small such that with A sphere the area of the sphere located at 2 R P from the center.In physical terms it means that the possible dissipation in the shell can be supplied with energy fluxes from below the shell (i.e., coming from the planet's ionosphere) or from above the shell (i.e., coming from the MS or stellar wind).
Ultimately, the energy flux is coming from the stellar wind, but the energy flux can be reflected or converted in the ionosphere and can be redirected away from the planet again.This integrated Poynting flux serves as a proxy for maximum available electromagnetic energy dissipated within an auroral acceleration region.Alternatively, we integrate the components of the Poynting flux parallel to the unperturbed magnetic field, B 0 = B − δB, where δB denotes the magnetic field perturbation generated by the interaction.These Poynting fluxes take into account the energy flux of which a fraction can directly contribute to particle acceleration and powering the ECMI-driven emission, where e B 0 is the unit vector pointing in direction of B 0 .We refer to this Poynting flux component as the auroral Poynting flux.As opposed to in Eq. ( 11), P a,|| (Eq.( 12)) serves as a more realistic estimator for calculating auroral energy dissipation since Eq. ( 11) includes a significant contribution of convected energy, which is likely not converted into particle acceleration.Table 2 summarizes integrated Poynting fluxes according to Eq. ( 11) (third column) for all three intrinsic magnetic field orientations.The 4th column shows integrated Poynting fluxes along the unperturbed field (Eq.( 12)).Integrated Poynting fluxes range from 3.5 × 10 18 down to 6.9 × 10 17 W for the open toward the closed MS model.Poynting fluxes along the unperturbed field (Eq.( 12)) amount to ∼9 × 10 17 and 10 17 W for the open and closed MS, respectively.The effect of magnetic topology on convected energy within the MS is therefore significant as the powers differ by almost one order of magnitude.Magnetic stress due to the stellar wind interaction can work on the MSs less strongly if the MS is closed, thus giving rise to weaker flows and therefore weaker convected Poynting fluxes.The trend is similar for Poynting fluxes along B 0 , P a,|| , but here the powers are reduced by almost an order of magnitude below the integrated total Poynting fluxes P a .

Influence of stellar wind variability on magnetospheric energetics
For modeling the space environment of τ Boötis b, the properties of its surrounding stellar wind carry very large uncertainties, in particular the stellar wind density.In Nicholson et al. (2016) and Vidotto et al. (2012), the coronal base density was estimated by choosing the electron density so that it can reproduce electron measure observations of τ Boötis A. The energy fluxes within the MS are powered by and limited by the maximum incident power of the stellar wind flow transferring onto the magnetospheric obstacle.Zarka (2007) found that the observed radio power of Solar System planets is nearly a constant fraction of the incident kinetic and magnetic energy convected through the obstacle's cross section, πR 2 mp , where R mp is the magnetospheric stand-off distance of the magnetized planet.This energy is utilized in perturbing the topology of the planets magnetic field, which in turn results in currents induced by changes in magnetic flux.Therefore, the incident power controls the energetics within the MS.The magnetic Poynting flux, P B , and the kinetic energy flux, P kin , convected through the obstacle's cross section can be calculated as follows: Additionally, the thermal energy flux should be considered as well as it cannot be neglected for close-in orbits where stellar wind temperature, T , pressure and density are high, with n sw being the stellar wind particle density and v 0 denoting the incident stellar wind velocity.
The magnetopause distance R mp can be obtained from an equilibrium between stellar wind and planetary ram (p ram = ρ sw v 2 0 ), magnetic (p B = B 2 /2µ 0 ) and thermal pressure.Both, the magnetospheric thermal and ram pressures are considered negligible, thus p ram,sw + p B,sw + p therm,sw = p B,pl , where the subscript sw stands for stellar wind and pl for planet.The magnetopause distance (or magnetospheric stand-off distance) can then be calculated from All parameters can be found in Table 1.The parameter v 0 refers to the relative velocity between the stellar wind and planet and B p to the planetary surface magnetic field at the equator.A certain fraction of the total incident power, intersecting with the magnetopause can eventually be converted for the generation of radio emission within the MS.The fraction of total incident energy, ϵ, that may result in radio emission is expected to range from 10 −5 to ∼3 × 10 −3 (i.e., P radio = ϵP i ) in the radiometric Bode law (see Fig. 3; Zarka 2007).We point out that various efficiencies for converting incident energy flux into electromagnetic radiation are discussed in the literature.For example, the efficiency of conversion from magnetospheric, auroral Poynting fluxes to radio emission, which accounts for A133, page 8 of 19 the efficiency of electron acceleration through wave-particle interaction and the efficiency of the electron-cyclotron maser, should be separated from the generic efficiency factor obtained from the radiometric Bode law (Zarka 2007; see Fig. 3).For Jupiter's radio emission the efficiency for conversion from magnetospheric, auroral Ponyting fluxes to radio emission is roughly 0.3-3 × 10 −4 (Saur et al. 2021).We denote this efficiency by ϵ a (Fig. 3).
As the stellar wind density is the most uncertain parameter we performed simulations with densities ranging from 0.05 • ρ 0 to 100 • ρ 0 (see Table 1 for the basic model).To get an understanding on how stellar wind variability affects the structure of the MS we show xz-plane slices similar to Fig. 1 for the two extreme cases (0.05 • ρ 0 and 100 • ρ 0 ) in Fig. C.1.We do not solve a self-consistent stellar wind model but instead follow the solar wind solution of Parker (1958) where the solution of the solar wind velocity v(r) is independent of the coronal base density n c,0 .In this solution, stellar mass and base temperature control v(r) and T (r), where r is the distance from the Sun.For simplicity of the parameter study of this subsection, we chose an isothermal approach and changed the density together with the pressure, p 0 (and therefore T ), according to p 0 ∝ ρ 0 (see Eq. ( 15)).We therefore kept the temperature constant, and consequently, according to Parker (1958), the velocity does not change.Given the average stellar mass loss rate of τ Boötis A of Ṁ ≈ 2.3 × 10 −12 M ⊙ yr −1 estimated by Nicholson et al. (2016) (see also our basic model, Table 1) the parameter range of stellar wind densities considered in this parameter study translates to mass loss rates between 1.15 × 10 −13 M ⊙ yr −1 and 2.3 × 10 −10 M ⊙ yr −1 since Ṁ ∝ ρ 0 .
We integrated the Poynting flux along the unperturbed field over a spherical shell with radius 2R p (e.g., Eq. ( 12)) in order to obtain an understanding of how much incident energy flux is eventually converted to auroral Poynting fluxes.Resulting powers are shown in Fig. 4 as a function of ρ 0 (and p 0 ).The simulated convected energy fluxes follow the trend of incident energy flux estimates (red solid line in Fig. 4) but are reduced to fractions of the total incident energy flux, P i , between 15 and 20% for the open MS and between 1 and 5% for the closed MS.Changes to  1 for the basic stellar wind model).The powers are calculated using Eqs.( 13)-( 15) and ( 17).The solid red line represents the sum of the three energy components.The MS stand-off distance, R mp (solid blue line) is calculated with the same set of ρ 0 and p 0 using Eq. ( 16).The stand-off distance can be obtained from the right y-axis, given in terms of planetary radii, R p .The magenta arrow indicates the convected power in the limiting case where the plasma density in the astrosphere approaches zero and the planet is only exposed to the stellar dipole field (see Sect. 4.1.2for discussion of this case).Red crosses show simulated Poynting fluxes along the unperturbed magnetic field over a spherical shell with radius 2R p for the closed MS (Eq.( 12)).Red stars indicate simulated Poynting fluxes for the closed MS.Bottom: Transfer function T a = P a /P i as a function of ρ 0 and p 0 (i.e., the conversion from incident energy to auroral Poynting fluxes).The scaling behavior according to power laws of T a is indicated by spectral indices above the lines.
the stellar wind density ρ 0 (and in the same manner p 0 ) affect the incident power inflicted on πR 2 mp but also influence the magnetospheric cross section in an opposite manner, as it can be seen in Fig. 4 (blue solid line).The magnetospheric stand-off distance scales according to R mp ∝ (ρ sw v 2 0 + p sw ) −1/6 and the incident energy flux with P i ∝ ρ sw + p sw ; therefore, the incident energy flux increase dominates over the effect of a shrinking MS due to increasing thermal and kinetic pressure.This is also validated by our simulation results (Fig. 4), implying an approximately linear scaling of auroral Poynting fluxes with ρ 0 and p 0 at least in the regime between 3 × 10 11 and 3 × 10 13 H + m −3 .Below the point where stellar wind magnetic energy dominates over A133, page 9 of 19 A&A 671, A133 (2023) thermal and kinetic energy near 10 11 H + m −3 , auroral Poynting fluxes seem to saturate near 2-3 × 10 17 W (open MS) and near 1-2 × 10 16 W (closed MS).Above 3 × 10 13 H + m −3 the increase of auroral Poynting fluxes with ρ 0 (and p 0 ) deviates further from the course of incident flux, implying a saturation toward 10 19 W (open MS).This, however, has to be validated further through future simulations.

Discussion
In this section, we discuss the importance of the stellar wind on magnetospheric energetics (Sect.4.1) and on possible auroral radio emission (Sect.4.2).

Impact of the stellar wind on magnetospheric energetics
In the following sections, we study the conversion of incident to dissipated power within the MS (Sect.4.1.1)as a function of stellar wind density and pressure.We also discuss the limiting case of an absent stellar wind (Sect.4.1.2).

Stellar wind variability, its effect on magnetospheric energetics, and the scaling behavior of auroral Poynting fluxes
We separate the considered stellar wind density and pressure range introduced in Sect.3.2.3,Fig. 4 in two regimes: Regime 1 ranges from a vanishing stellar wind up to a density at roughly 10 11 H + m −3 , where kinetic and thermal energy fluxes fall below the persistent magnetic energy flux, which dominates the flow (compare red curves in Fig. 4).Above roughly 10 11 H + m −3 the flow is super-Alfvénic (M A ≈ 2) and super-fast (M f ≈ 1).The interaction is super-Alfvénic for the whole parameter space used in our simulations and sub-fast only for the lowest simulated density (ρ sw = 7 × 10 10 H + m −3 , M A ≈ 1.2).The incident energy nearly stagnates below ρ 0 = 10 10 H + m −3 (red dotted line).Below this point the incident energy flux asymptotically approaches its minimum at 8 × 10 17 W, as we assume that only the plasma density decrease but the incident magnetic field is kept constant.In this regime, it can be expected that the stellar wind magnetic field solution transitions from the Parker solution (e.g., B ∝ r −2 ) to a pure stellar multipole (here dipole) solution (e.g., B sw = B star ∝ r −3 ) with decreasing stellar wind density.Eventually, when the stellar wind density hypothetically approaches zero, only the dipolar stellar magnetic field interacts with the planetary magnetic field.This limiting case will be separately discussed in Sect.4.1.2.Regime 2 ranges from roughly 10 11 H + m −3 up to arbitrarily high stellar wind densities.Here kinetic and thermal energy fluxes dominate the flow.We now focus on this regime.Considering the total energy flux convected through the magnetospheric cross section πR 2 mp , P total (red solid line in Fig. 4), we observe a nearly constant efficiency of conversion from incident stellar wind energy toward magnetospheric Poynting fluxes at auroral altitudes (we assumed r ≈ 2R p ) with increasing density and pressure.We calculate the transfer function T a as the conversion ratio from total incident energy flux P i to the simulated auroral Poynting fluxes parallel to the unperturbed field (Eq.( 12)), P a (red crosses and stars in Fig. 4), within the MS, such that T a = P a /P i (see Fig. 3 for a schematic illustrating the role of T a ).The transfer function also contains information on the magnetic topology and thus the efficiency of reconnection.The transfer function is displayed in the lower panel of Fig. 4. For the open MS, auroral Poynting fluxes decrease as a function of ρ 0 and p 0 according to an approximately constant ratio up to ∼3 × 10 13 H + m −3 .For higher densities and pressures, the transfer function scales with an exponent of ∼−0.4,indicating a decrease of efficiency for conversion from incident to auroral energy fluxes.As the MS is increasingly compressed due to stronger ram and thermal pressures, the magnetopause eventually crosses the spherical shell with r = 2R p after a critical density of ∼3 × 10 13 H + m −3 and pressure of 3 × 10 −4 Pa.
For the closed MS, T a behaves differently.The transfer function instead increases slightly from 7 × 10 10 to ∼3 × 10 13 H + m −3 following a power law with exponent ∼0.6.The opposite behavior compared to the open MS transfer function might be a consequence of a geometry where less reconnection occurs.The stronger ram and thermal pressure exert a stronger tension on planetary field lines, which in turn release stronger energy fluxes during reconnection.

Beyond the MHD limit
The more the stellar wind density decreases, the emptier the heliosphere of τ Boötis A becomes.In analogy with the solar wind, the stellar wind density upstream of the MS of τ Boötis b may change by orders of magnitude.The solar wind density is observed to vary by more than two orders of magnitude (see, e.g., Chané et al. 2012) such that Earth's bow shock can disappear and the Earth develops Alfvén wings.In the hypothetical limit when the density approaches zero, τ Boötis b will still be exposed to the stellar magnetic field B τ (r) (which then decreases according to B τ (r) ∝ r −3 instead of B sw ∝ r −2 ) and will propagate through it.Therefore, the incoming Poynting flux of the star does not, in contrast to kinetic and thermal energy flux, vanish and is maintained by the relative motion between the stellar and planetary magnetic field.In the case of an empty heliosphere, the interaction around τ Boötis b is no longer magnetohydrodynamic; it turns electromagnetic.Then the movement of τ Boötis b within the external magnetic field of the star is a unipolar or homopolar interaction (i.e., a moving conductive object with external field similar to a current generator in classical electromagnetism).In the case of a stellar magnetic field rotating at the same speed as the planet orbits around its host star (i.e., in the case of total orbital and spin synchronization), nonexistent changes in magnetic flux lead to a system where no work can be done by the magnetic fields.Therefore, no currents are induced and the MS remains energetically silent.Although it is expected that close-in exoplanets are in nearly tidally locked rotation, Murray & Dermott (2000) suggest in their textbook that completely synchronous rotation might not be possible if the planet has no permanent magnetic quadrupole moment and its eccentricity is not zero, which is the case for τ Boötis b (Wang & Ford 2011).Taking the rotation period of τ Boötis A, P τ = 3.1 ± 0.1 days (Brown et al. 2021;Mengel et al. 2016), and the sidereal rotation period of τ Boötis b, P orb = 3.31 days (Butler et al. 1997;Wang & Ford 2011), we can calculate the relative velocity between the stellar and planetary magnetic field v 0 = 2πa × (1/P τ − 1/P orb ) ≈ 10.4 km s −1 , where a = 0.046 AU is the semimajor axis of τ Boötis b.We use an average surface magnetic flux density of τ Boötis A obtained by Marsden et al. (2014); Mengel et al. (2016); Jeffers et al. (2018), B τ,0 = 2.4 G, and calculate the flux density at 0.046 AU using the dipole formula, B τ (r = a) = B τ,0 • (a/R τ ) −3 ≈ 0.73 G, where R τ = 1.43 R ⊙ is the stellar radius (Bonfanti et al. 2016).The magnetospheric stand-off distance in this case is defined purely by the balance of stellar and planetary magnetic pressure: A133, page 10 of 19 F. Elekes and J. Saur: A&A proofs, manuscript no.aa44947-22 Now we can calculate the stellar magnetic power convected on the MS using This value is included in Fig. 4 as magenta arrow.We can conclude the following: (a) The maximum possible magnetospheric stand-off distance is reached in complete absence of a stellar wind and if the stellar and planetary magnetic fields are antiparallel.If the stellar and planetary fields are parallel (i.e.open MS) the planetary field lines are all connected to the stellar field and no magnetopause can be determined.In the closed MS case the magnetopause lies at roughly 8.6 planetary radii in the upstream direction.
(b) Even if there is no stellar wind, the magnetic interaction between the stellar and planetary magnetic field still has the potential to drive an interaction with an available power limit of roughly 8 × 10 15 W due to the relative motion of τ Boötis b in the stellar magnetic field.Radio emission would still be possible, although it would be very weak; corresponding radio fluxes at Earth's position would be far below today's telescope sensitivity limit.Considering Poynting flux-to-radio power efficiencies between 10 −4 and 10 −2 , radio powers can reach values between 10 11 and 10 14 W.These emitted powers exceed the strongest radio sources within the Solar System by several orders of magnitude, with Jupiter's aurora being the strongest radio emitter (P radio ≈ 10 10 -10 11 W, Zarka 2007), although Jupiter's emission is, in contrast to close-in exoplanets, powered by internal, rotationally driven mechanisms.The transition from a MS interaction with a stellar wind field (according to the Parker solution, Parker 1958) to an interaction with a pure dipolar stellar magnetic field goes with an energetic transition followed by a decrease of maximum emitted radio power.This might pose a possible opportunity for constraining stellar wind densities in the future.Solving a self-consistent stellar wind model and comparing auroral Poynting fluxes for different stellar wind base densities could reveal the critical density range where the transition from a stellar magnetic-field-dominated electrodynamic interaction to a stellar-wind-dominated MHD interaction takes place.Comparing the magnetospheric Poynting fluxes and corresponding radio powers with possible future observations could reveal if the stellar wind density lies below or above the critical density.

Magnetospheric Poynting fluxes and auroral radio emission
In this section, we discuss possible radio emission scenarios (Sect.4.2.1),study how magnetic topology of the interaction as well as stellar wind variability affects auroral radio emission output (Sect.4.2.2).We also discuss the possibility of a sub-Alfvénic emission scenario (Sect.4.2.3) as well as a rotation-driven MS of τ Boötis b in Sect.4.2.4,followed by a discussion on possible source regions and radio frequencies of auroral emission in Sect.4.2.5.

On the different radio emission scenarios
There are several scenarios capable of generating observable radio emission that must be distinguished from each other: (1) If the stellar wind is sub-Alfvénic, Alfvén waves are able to propagate upstream toward the stellar atmosphere along Alfvén wings and possibly drive electron acceleration and radio emission in the stellar vicinity.Local radio emission within the MS can also be generated in this scenario.We refer to this scenario as sub-Alfvénic emission scenario.This scenario is discussed briefly in Sect.4.2.3.
(2) If the stellar wind is super-Alfvénic, no MHD wave is able to propagate upstream.The stellar wind-MS interaction, however, drives Poynting fluxes within the MS, which may, to some extent, generate auroral radio emission (see Fig. 3).We refer to this scenario as the magnetospheric emission scenario.Due to the stellar wind being super-Alfvénic for all simulations, we focus on the magnetospheric emission (Sect.4.2.2).
(3) In a rotation-dominated MS scenario the rotating planet and its magnetic field causes corotation of magnetospheric plasma that, at some point, breaks down due to conservation of angular momentum if radial mass transfer takes place.This corotation breakdown exerts magnetic stresses on the field lines that are the root cause of auroral Poynting fluxes, which in turn drive auroral radio emission.This scenario is discussed in Sect.4.2.4.

Effects of magnetic field tilt and stellar wind variability on auroral radio emission: Magnetospheric emission scenario
Figure 5 shows radio powers as a function of magnetic axis tilt.Radio powers are obtained by multiplying integrated auroral Poynting fluxes (i.e., Eq. ( 12)), which serve as a proxy for the maximum available electromagnetic energy that is transported along magnetic field lines, by efficiency factors for converting magnetospheric Poynting fluxes to radio power, ϵ, ranging from 10 −4 to 10 −2 .This range covers proposed (Zarka 2007) and observed efficiency factors (e.g., ϵ ≈ 10 −4 for Jupiter, Saur et al. 2021).The modeled magnetic field tilt can also be interpreted as stellar magnetic field orientation within this work, allowing us to study the effect of varying stellar magnetic field polarity on magnetospheric Poynting fluxes and limits for associated radio emission.Radio powers within the limits inferred from observations by Turner et al. (2021) 2019) is plotted as a yellow line.The vertical gray line marks the basic model (Table 1).
than in the Jovian MS (Saur et al. 2021).Electric fields generated by reconnection between stellar wind and planetary magnetic field lines are expected to contribute significantly to powering electron acceleration and therefore the ECMI (Jardine & Collier Cameron 2008).In our studies we find reconnection to indirectly play an important role (Fig. 5) because auroral Poynting fluxes and consequently radio powers drop by nearly an order of magnitude from an open to a closed MS.This is due to magnetic stress exerted by the stellar wind interaction being less strong for closed MSs.The polarity of τ Boötis A's magnetic field switches every approximate 360 days (Fares et al. 2013).Shorter cycles in magnetic activity levels (by means of S-indices) were also observed (Mengel et al. 2016).A difference of half an order of magnitude to almost an order of magnitude can therefore be caused by a polarity reversal of τ Boötis A's magnetic field.This results in radio emission whose observability is expected to fluctuate periodically in a nearly 1-year cycle.We note that the stellar wind magnetic field strength was kept constant in our parameter study, although in reality the field strength may vary strongly and influence produced radio emission significantly (See et al. 2015).The emitted radio flux observed at Earth's position can be calculated with (Grießmeier et al. 2005(Grießmeier et al. , 2007b) where Ω is the solid angle of the beam and δν the emission bandwidth that is approximately equal to the maximum gyrofrequency (Grießmeier et al. 2007b), ν g,e ≈ 24 MHz.We assume a solid angle of Ω = 1.6 sr, similar to Jupiter's decametric radio emission (Zarka et al. 2004).The distance to the τ Boötis system is 15.6 pc.We calculate the radio flux for both, the open and closed MS model, as a function of ρ 0 and p 0 with radio efficiencies ϵ a = 10 −4 -10 −2 .The results are displayed in Fig.  (2021).Horizontal gray lines indicate theoretical sensitivity limits of the LOFAR telescopes for 20MHz.As stated by Turner et al. (2019), the realistic sensitivity might be slightly lower for circularly polarized (Stokes V) signals.We therefore include the expected sensitivity calculated by Turner et al. (2019) as yellow line.The results in Fig. 6 indicate that radio efficiencies between ∼3 × 10 −3 and ∼10 −2 are most consistent with the tentative observations (Turner et al. 2021).The efficiency accounts for several steps from conversion of auroral Poynting fluxes to radio emission (e.g., wave-particle interaction, electron acceleration and ECMI), and therefore an efficiency on the order of 1-10% might be unrealistic.The efficiency for Jupiter's auroral emission is roughly ϵ a = 0.3-3 × 10 −4 (Saur et al. 2021), and therefore ϵ a = 10 −2 might be too high.Moreover, high plasma densities within the MS injected by the dense stellar wind and due to strong irradiation, which results in high ionization rates and inflated atmospheres (e.g., for ν And b, see Erkaev et al. 2022), may further decrease the ECMI efficiency or even prevent it (Weber et al. 2017(Weber et al. , 2018;;Daley-Yates & Stevens 2018).
Assuming the radio efficiency to lie near 10 −3 , the radio flux from a closed MS falls below the detection threshold (yellow line).Therefore, in the case of a polarity reversal of τ Boöts A's magnetic field (i.e., from aligned with the planetary field to antialigned), the radio signal would not be observable anymore in the case of radio efficiency equal or below ∼3 × 10 −3 .In the case of ϵ a ≈ 10 −4 all radio fluxes for the basic model fall below the sensitivity limit.The observability increases, however, if stellar wind density and pressure rises, rendering ϵ a = 10 −3 -10 −2 to possible efficiencies to observe emission from open and closed MSs.Additionally, the ECMI efficiency (Treumann 2006;Weber et al. 2017) as well as efficiency of electron acceleration through wave-particle interaction decreases dramatically with increasing plasma density (Saur et al. 2018), making the higher density and pressure regime a less likely scenario to explain the tentative observations.As the pressure rises, the magnetopause is getting closer to the planet, reducing the space of magnetospheric diluted plasma regions between the magnetopause and atmosphere where radio emission might occur.We therefore conclude that the basic model (vertical gray line) and slightly different configurations represent the most likely scenarios if the emission is indeed generated by stellar wind-driven auroral Poynting fluxes.In this case radio emission is only observable, if the stellar wind and planetary magnetic fields are aligned (i.e., the MS is open).Given the high efficiencies (ϵ > 10 −3 ) needed by our model in order to generate radio emission that is consistent with the tentative observations, the magnetospheric emission scenario might not be energetic enough to explain the observations.

Sub-Alfvénic emission scenario
Although there is no sub-Alfvénic interaction within the parameter space we considered, the possibility of such an interaction and its consequences on possible radio emission should not be neglected.By choosing a stellar wind density of ρ sw = 0.03 ρ 0 we find an Alfvénic Mach number of M A ≈ 0.9.In this case, Alfvén waves may propagate back to the star through Alfvén wings connecting the planetary magnetic field with the star.The electromagnetic energy channeled through this flux tube can be calculated using the model from Saur et al. (2013), where θ = 0 • is the angle that describes the deviation of the flow from being perpendicular to the stellar wind magnetic A133, page 12 of 19 F. Elekes and J. Saur: A&A proofs, manuscript no.aa44947-22 field, R mp ≈ 5 R p the magnetospheric stand-of distance and ᾱ the interaction strength.Due to the planet presumably possessing an ionosphere, which favors a strong plasma interaction, we chose ᾱ ≈ 1.Using v 0 and B sw from our basic model (Table 1), the energy flux channeled through Alfvén wings parallel to magnetic field lines amounts to P S = 1.2 × 10 18 W.This energy flux may contribute to electron acceleration and consequently to ECMI-driven radio emission in the space environment near the star.Taking the same range of Poynting flux to radio efficiencies (10 −4 -10 −2 ), the expected radio power ranges from roughly 10 14 to 10 16 W, which partially overlaps with the range of radio power deduced by Turner et al. (2021).The emission frequency would, however, be much lower compared to the signals observed by Turner et al. (2021) due to the low magnetic field strength of τ Boötis A (Turner et al. 2021).Additionally, Turner et al. (2021) state that stellar emission would be less strongly circularly polarized compared to planetary emission.Nevertheless the root energy flux able to possibly drive radio emission near the stellar atmosphere from such a sub-Alfvénic interaction exceeds the auroral Poynting fluxes in the magnetospheric emission scenario by roughly half an order of magnitude.Consequently this scenario might provide energetically a possible explanation for the observed radio emission if the interaction is indeed sub-Alfvénic.This scenario requires the stellar wind density to be reduced by at least a factor of ∼30 compared to our basic model.This, however, could well be within the realm of possible densities given the observed density variations of the solar wind even though the intervals of variability are short (Chané et al. 2012).

Rotationally dominated magnetospheric emission scenario
One of the possible generators for auroral emissions is radial mass transport within the MS as it is the case for Jupiter (e.g., Hill 2001;Zarka et al. 2018Zarka et al. , 2021)).Plasma is moving radially outward due to centrifugal forces within an exoplanet's MS.
For small distances from the planet (or stellar host) the plasma corotates with the host.After a certain distance L in units of planetary radii, L × R p , the corotation breaks down due to conservation of angular momentum (i.e., the plasma orbits its host with a smaller angular velocity than those from the rotation of the host).The relative velocity exerts magnetic stresses on the frozen-in field lines, the tension accelerates plasma along the field lines back to the host.The distance of corotation breakdown can be estimated using the colatitude of the region, where auroral emission occurs, θ a , following Hill (2001), sin(θ a ) = L−1/2 .Saur et al. ( 2021) derived a so-called auroral power potential for magnetized rotation-dominated hosts, where the subscript host refers in our case to the planet τ Boötis b.The hosts magnetic flux density is denoted by B host , Ω host ≈ 2π/P orb = 2.2 × 10 −5 rad s −1 is its angular velocity and R host = R p the planet's radius.We note that the discussion in this subsection is only based on theoretical considerations since we neglected planetary rotation in our MHD model.Following Saur et al. (2021) we can calculate the energy flux due to mass transport, P mag , by using where Σ p denotes the Pedersen conductance.We assume a value of 1 S and 10 S that are of the same order of magnitude as observed conductances near Jupiter (e.g., Hill 2001;Strobel & Atreya 1983;Hinson et al. 1998).For a hot Jupiter exoplanet the ionospheric Pedersen conductivity might be larger due to the larger ionization (Koskinen et al. 2010).However, the height of the resulting ionospheric layer and thus the mobility of the electrons and ions within the atmosphere also play an important role for the values of its conductivity.To our knowledge no dedicated study for the conductance of τ Boötis b is available.We also note that in Eq. ( 22) the distance of corotation breakdown L is proportional to Σ 1/4 (Hill 2001).Thus, the P mag in Eq. ( 23) depends effectively weaker on Σ p (i.e., P mag is proportional to Σ 1/2 p ).In our simulations we integrate the ion component of the Pedersen conductivity along the z-axis at the poles from the surface (r = R p ) to r = 1.1 R p where the plasma density peaks and find the height integrated Pedersen conductance to be on the order of 10 S. We find, with aurora occurring at field lines with L-parameter smaller than ∼4 that have to lie within the MS (e.g., L ≈ 3.5), the auroral power potential to be S pot ≈ 0.5 W m −2 S −1 and the magnetic power to lie between P mag ≈ 10 14 W and 10 15 W. The auroral power potential of Jupiter, for comparison, is roughly 600 W m −2 S −1 .The resulting energy flux between the MS and ionosphere lies near the Jovian value (3.1 × 10 14 W, Saur et al. 2021).This result is several orders of magnitude below the power of magnetospheric Poynting fluxes powered by the stellar wind-planet interaction in our model (e.g., 10 17 -10 18 W, Table .2).Therefore, we can safely assume that the MS of τ Boötis b is not rotationally dominated, which is in agreement with the expected slow rotation nearly synchronized with its orbit around the star.

Properties and source regions of possible auroral radio emission
Auroral radio emission is mainly characterized by signal strength and frequency.The peak frequency is defined by the local electron gyrofrequency, where e is the electron charge, m e its mass and B the local magnetic flux density (Farrell et al. 1999;Grießmeier et al. 2007a,b).We already discussed possible radio powers from τ Boötis b in Sects.3.2.2 and 4.2.2, now we focus on possible signal frequencies.Turner et al. (2021) observed signals within the 21-30 MHz range.The magnetic flux density inferred from the maximal signal frequency by the authors is almost certainly underestimated since ECMI-driven auroral radio emission typically arises from altitudes above the polar regions of the planet's atmosphere (Zarka 1998;Hess & Zarka 2011) due to three reasons: (a) High ionospheric electron densities cause the electron plasma frequency, ν p,e , to exceed the local gyrofrequency, ν g,e , thus prohibiting the transmission of emitted waves; (b) radio wave amplification through the ECMI mechanism works most efficiently in regions where ν g,e >> ν p,e (Treumann 2006;Weber et al. 2017) and (c) the electrons are usually accelerated in regions where the plasma density along a field line is the smallest (or B/n is maximum; e.g., Saur et al. 2018).There wave particle interaction to accelerate electrons is maximum efficient.The electron plasma frequency is given by where n e = n i is the particle density in a quasi-neutral plasma and ϵ 0 the vacuum permittivity.The source regions of Jupiter's A133, page 13 of 19 A&A 671, A133 (2023) radio emission lie several planetary radii above the surface (e.g., Zarka 1998) where strongest electron acceleration is expected (e.g., Mauk et al. 2017).The maximum gyrofrequency in our model (directly above the magnetic poles) is ∼24 MHz.Auroral emission, however, typically arises from near the OCFB.In our simulations we indeed see that strongest Poynting fluxes are confined to narrow bands encircling the OCFB (or magnetopause) near L = 3 and L = 3.5 shell field lines (Fig. 2).Since source regions of ECMI induced emission typically lie at higher altitudes the magnetic flux density and thus emission frequency is even further reduced.The conditions (a-c) for efficient ECMIdriven radio emission might particularly be fulfilled for regions with low plasma density at the night side of the planet where the shadow exerted by the planet prevents photoionization.Since τ Boötis b is likely in tidally locked rotation around its host star, the night side might exhibit relatively low temperatures, favoring a recombination of electrons and hydrogen ions as well as further electron capture processes.In the shadow zone (i.e., permanent night side of the planet), ECMI induced emission might be more likely and might as well occur in lower altitudes, where the magnetic field is stronger (and thus emission frequency is higher).

Conclusions
In this study, we have modeled the MS of the hot Jupiter exoplanet τ Boötis b and its interaction with the stellar wind.In order to study magnetospheric energetics by means of Poynting fluxes, we performed MHD simulations of the near space environment of τ Boötis b with the stellar wind being injected into the simulation domain.The magnetic field of τ Boötis b was constrained by tentative observations of auroral radio emission (Turner et al. 2021).The stellar wind model is based on simulations where observed magnetic surface maps of τ Boötis A were utilized (Nicholson et al. 2016;Vidotto et al. 2012).We investigated magnetospheric energetics, available electromagnetic power, and limits of possible radio emission originating from above the planet's polar regions.We also explored the possibility of a rotation-dominated magnetospheric emission scenario as well as a sub-Alfvénic (stellar) emission scenario.The effect of stellar wind variability in terms of density, pressure, and magnetic field orientation on magnetospheric energetics was additionally investigated.We find that the stellar wind-MS interaction is super-fast and super-Alfvénic for the majority of stellar wind densities and pressures considered in this study (from 1.4 × 10 11 H + m −3 and 2.5 × 10 −6 Pa to arbitrarily high stellar wind densities and pressures).For the minimum density and pressure considered in this work (ρ sw = 7 × 10 10 H + m −3 and p sw = 1.3 × 10 −5 Pa), the interaction is super-Alfvénic and sub-fast.Maximum available electromagnetic energy convected within the MS amounts to 3.5 × 10 18 W, 1.8 × 10 18 W, and 7 × 10 We present Poynting flux maps of the planet's near space environment.The strongest Poynting fluxes are confined to narrow bands encircling the OCFBs (i.e., the magnetopause) that are displaced toward the planet's downstream side due to magnetic stresses exerted by the stellar wind and its high plasma β.Electromagnetic energy transported parallel to unperturbed field lines is also confined to bands near the magnetopause and the polar open-field-line regions.
Considering the stellar wind density as a free parameter, we chose values between 0.05 and 100 times the basic density, ρ sw = 1.4 × 10 12 H + m −3 (Table 1), and changed the stellar wind pressure in the same manner in order to keep the wind isothermal.Comparing the incident total energy flux convected through the magnetospheric cross section, we find that the transfer function, T a (i.e., the amount of stellar wind incident energy flux converted to magnetospheric Poynting fluxes parallel to unperturbed field lines) amounts to a fraction of roughly 20 % for the open and 2 % for the closed MS.We find that T a ∝ ρ −0.06 0 and T a ∝ ρ 0.57 0 for the open and closed MS, respectively, in the density range 7 × 10 11 -2×10 13 H + m −3 .Increasing the stellar wind ram and thermal pressure and thus the incident energy flux causes the MS to shrink.The reduced capability of the MS for receiving stellar wind energy is outperformed by the increasing pressure and density.The scaling of T a and thus the energy output of the MS nearly stagnates from that point on: the magnetopause becomes so small that the planetary atmosphere starts to interact with the pre-bow-shock and magnetosheath flow.
We investigated different radio emission scenarios.A rotation-dominated magnetospheric emission scenario can be ruled out due to the expected long rotation period of τ Boötis b.The estimated magnetic power generated by the planet's rotation is on the order of P mag ≈ 10 14−15 W. As this would require a very high Poynting flux-to-radio efficiency of 10 −2 , we expect that this scenario is not powerful enough to drive the tentatively observed radio emission.
We considered several efficiencies for the conversion of auroral Poynting fluxes to radio emission and compared the results with observations together with theoretical limits based on the tentative observations by Turner et al. (2021).We find that the generation of radio emission near τ Boötis b requires relatively high efficiencies (ϵ a ≈ 1-5 × 10 −3 ) compared to estimates from Jupiter's rotationally driven aurora (up to ∼3 × 10 −4 ; Saur et al. 2021) in order to explain the observed emission if we assume the emission to originate from the MS of τ Boötis b.Furthermore, we find a strong dependence of radio emission observability on MS topology.In the case of a polarity reversal of τ Boötis A's magnetic field from an aligned to anti-aligned configuration (with respect to the planetary field), the estimated radio power falls below the observational limit.We therefore expect an on-off nature of detectable radio emission in phase with τ Boötis A's magnetic cycle.
We also considered the case of a sub-Alfvénic interaction, which requires the stellar wind density to be at least a factor of 30 lower compared to our basic model.The Poynting flux channeled toward the host star due to this interaction amounts to ∼10 18 W, which slightly exceeds the power converted by the magnetospheric emission scenario.
A133, page 1 of 19 Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0),which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.This article is published in open access under the Subscribe to Open model.Subscribe to A&A to support open access publication.
Fig. 1.Velocity fields (colored arrows, left colorbars) and plasma pressure (color contours, right colorbars) in the xz-plane for the open MS (θ B = 0 • , top), semi-open MS (θ B = 90 • , middle), and closed (θ B = 180 • , bottom) MS cases.Projected stellar wind magnetic field lines are indicated as solid black lines within the xz-plane parallel to the ambient magnetic field.Closed and open magnetospheric field lines are colored in magenta.

Fig. 2 .
Fig. 2. Mercator projections of the Poynting flux (upper row), plasma velocity (middle row), and absolute values of Poynting flux components parallel to the unperturbed planetary field plus the small background field from stellar wind (bottom row).The results are shown at an altitude of one planetary radius above the surface.The left column displays maps for the open MS case (θ B = 0 • ), the right column for the closed MS case (θ B = 180 • ).Arrows indicate normalized angular components, and color contours denote radial components.Dashed red lines indicate the location of the OCFB.
3.2.1.Spatial structureWe first describe the spatial structure of the Poynting fluxes and plasma velocities within the MS as displayed in the top and middle row of Fig. 2, respectively.The plots show Mercator projections of the angular vector components over a spherical shell with radius 2 R p .The figure shows the colatitude and longitude as well as the orientation of vectorial quantities and their magnitudes.Positive values indicate radial components pointing away from the planet.Also shown are the open-closed planetary field line boundaries (OCFBs).Magnetic field lines with both foot points on the planetary surface are closed field lines.Each field line that has only one foot point on the planet is an open field line.The OCFB separates areas with open from areas with closed field lines.Thus, the OCFB also represents the magnetopause at that specific radial location.Both open and closed MS cases are shown left and right, respectively.
left).The very narrow vertical extent of the downstream closed field line region of the open MS is caused by magnetic tension due to the magnetized stellar wind.Highest velocities are found within the open field line region mainly at the downstream side where plasma is accelerated downstream through magnetic tension on open lines.Strong Poynting fluxes occur where plasma velocities have strong components perpendicular to the magnetic field.They are found within the open field line region mainly at the downstream side with outward directed Poynting fluxes.Comparatively strong Poynting fluxes, but directed toward the planet, are located on the upstream side near the magnetopause.Within the closed field line region and especially near the equator Poynting fluxes mostly vanish.
shows |S • e B 0 | for the open MS (left) and closed MS (right).We note that only absolute values are shown in the plots in order to clearly identify zero or near-zero power densities.Strongest energy transport along unperturbed field lines occurs over narrow bands encircling the polar open field line regions at the flanks of the planet where velocities are nearly perpendicular to the magnetic field as seen in Fig. 2.Moreover, the spatial structure of Poynting fluxes along B 0 is strictly symmetric with respect to the equator (at θ = 90 • ).A significant amount of energy is transported parallel to the unperturbed field within the polar open field line regions in the open MS case.Parallel energy fluxes reach values up to 10 W m −2 at the flanks of the planet just outside the closed field line regions.Poynting fluxes up to 9 W m −2 are found at the downstream side, above the OCFB.For both, open and closed MS model, strongest convected energy can be found extensively in high latitude regions due to high velocities perpendicular to the magnetic field.Here the planetary field lines are most mobile in a sense that they are bent over toward the downstream side by the stellar wind.For the closed MS parallel Poynting fluxes up to roughly 9 W m −2 can be found directly at the planetary poles slightly shifted toward the upstream side.At lower latitudes parallel Poynting fluxes up to 6 W m −2 are confined to narrow bands at the flanks of the planet.Auroral emission is expected to be strong where Poynting fluxes are large, hence near the OCFB (e.g., mostly confined to the L = 3-3.5 shell at the upstream side) and in the polar regions for both MS models.They vanish completely along the equator.Generally said Poynting fluxes are significantly weaker and confined to the small polar regions for the closed MS model compared to the open MS case.In the open MS model strong parallel Poynting fluxes cover the whole open field line area with their maximum at the flanks of the planet in contrast to the closed MS where the regions of strongest parallel Poynting fluxes are partitioned into smaller areas around the planetary poles.

Fig. 3 .
Fig. 3. Schematic illustrating the several steps from incident stellar wind energy flux toward auroral radio emission.The transfer function, T a (see Sect. 4.1.1)describes the conversion from incident stellar energy to auroral Poynting fluxes.The conversion efficiency from auroral Poynting fluxes (Eq.(12)) to radio emission, ϵ a , implicitly includes the efficiency of electron acceleration and the ECMI mechanism.The steps within the pink shaded area are not included in our model.Brown arrows indicate physical processes, blue arrows denote model parameters that quantify energy conversion, and the magenta arrow the radiometric scaling law.

Fig. 4 .
Fig. 4. Incident energy fluxes and auroral Poynting fluxes as a function of stellar wind density and pressure.Top: analytically calculated incident kinetic (dashed red line), Poynting (dotted red line), and thermal (dashed-dotted red line) energy fluxes convected through the magnetospheric cross section, πR 2 mp , as a function of stellar wind density at the position of τ Boötis b.The pressure varies accordingly as p 0 ∝ ρ 0 (see Table1for the basic stellar wind model).The powers are calculated using Eqs.(13)-(15) and (17).The solid red line represents the sum of the three energy components.The MS stand-off distance, R mp (solid blue line) is calculated with the same set of ρ 0 and p 0 using Eq.(16).The stand-off distance can be obtained from the right y-axis, given in terms of planetary radii, R p .The magenta arrow indicates the convected power in the limiting case where the plasma density in the astrosphere approaches zero and the planet is only exposed to the stellar dipole field (see Sect. 4.1.2for discussion of this case).Red crosses show simulated Poynting fluxes along the unperturbed magnetic field over a spherical shell with radius 2R p for the closed MS (Eq.(12)).Red stars indicate simulated Poynting fluxes for the closed MS.Bottom: Transfer function T a = P a /P i as a function of ρ 0 and p 0 (i.e., the conversion from incident energy to auroral Poynting fluxes).The scaling behavior according to power laws of T a is indicated by spectral indices above the lines.

Fig. 5 .
Fig. 5. Expected radio powers as a function of planetary magnetic axis tilt using auroral-to-radio power conversion efficiencies between 10 −4 and 10 −2 .The auroral Poynting fluxes, S • B 0 , are integrated over a spherical shell with radius 2R p .The gray shaded area represents the observational limits given by Turner et al. (2021).

Fig. 6 .
Fig. 6.Radio flux (Eq.(20)) as a function of stellar wind density and pressure for different efficiency factors.Colored solid and dashed lines represent fluxes for the open and closed MS models.Observational limits (Turner et al. 2021) are indicated by the gray shaded area.Horizontal gray lines display theoretical sensitivity limits of the LOFAR telescope.The real sensitivity for Stokes V signals obtained from Turner et al. (2019) is plotted as a yellow line.The vertical gray line marks the basic model (Table1).
6. Solid and dashed colored lines represent radio fluxes originating from open and closed MSs, respectively.The gray shaded area again denotes the range of observed radio fluxes form Turner et al.
17 W for an open, semiopen, and closed MS configuration.Auroral Poynting fluxes at altitudes of 1R p above the planet reach powers of 8.7 × 10 17 W down to 1.1 × 10 17 W for the open and closed MS, indicating a strong decrease in stellar wind energy converted to auroral Poynting fluxes as a function of magnetic field tilt (either stellar or planetary).

Table 2 .
Integrated magnetospheric Poynting fluxes for different magnetic field topologies.ModelDipole tilt P a (W) (a) P a ||B 0 (W)(b) lie within the gray shaded area.It is visible that efficiency factors in the range of ϵ ≈ (0.3 − 1) × 10 −2 deliver radio powers most consistent with observations if the MS is open or at least semi open given the basic model (Table1).This indicates that the efficiency of auroral Poynting fluxes driving electron acceleration and the electron cyclotron maser emission may be higher in the MS of τ Boötis b