Issue 
A&A
Volume 557, September 2013



Article Number  A31  
Number of page(s)  13  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201321790  
Published online  20 August 2013 
Starplanet magnetic interaction and evaporation of planetary atmospheres
INAFOsservatorio Astrofisico di Catania, via S. Sofia 78, 95123 Catania, Italy
email: nuccio.lanza@oact.inaf.it
Received: 29 April 2013
Accepted: 2 July 2013
Context. Stars interact with their closein planets through radiation, gravitation and magnetic fields.
Aims. We investigate the energy input to a planetary atmosphere by reconnection between stellar and planetary magnetic fields and compare it to the energy input of the extreme ultraviolet (EUV) radiation field of the star.
Methods. We quantify the power released by magnetic reconnection at the boundary of the planetary magnetosphere that is conveyed to the atmosphere by accelerated electrons. We introduce simple models to evaluate the energy spectrum of the accelerated electrons and the energy dissipated in the atmospheric layers in the polar region of the planet upon which they impinge. A simple transonic isothermal wind flow along field lines is considered to estimate the increase in mass loss rate in comparison with a planet irradiated only by the EUV flux of its host star.
Results. We find that energetic electrons can reach levels down to column densities of 10^{23}−10^{25} m^{2}, comparable with or deeper than EUV photons, and increase the mass loss rate up to a factor of 30−50 in closein (< 0.10 AU), massive (≥1.5 Jupiter masses) planets. Mass loss rates up to (0.5−1.0) × 10^{9} kg s^{1} are found for atmospheres heated by electrons accelerated by magnetic reconnection at the boundary of planetary magnetospheres. On the other hand, average mass loss rates up to (0.3−1.0) × 10^{10} kg s^{1} are found in the case of magnetic loops interconnecting the planet with the star.
Conclusions. The starplanet magnetic interaction provides a remarkable source of energy for planetary atmospheres, generally comparable with or exceeding that of stellar EUV radiation for closein planets. Therefore, it must be included in models of chemical evolution or evaporation of planetary atmospheres as well as in modelling of light curves of transiting planets at UV wavelengths.
Key words: planetstar interactions / planets and satellites: atmospheres / planets and satellites: magnetic fields / stars: magnetic field / magnetic reconnection / acceleration of particles
© ESO, 2013
1. Introduction
Several hundreds of extrasolar planets are presently known, about 20 percent of which are Jupitermass objects orbiting their host stars within 0.1 AU (hot Jupiters, hereafter HJs). Since latetype stars have outer convection zones that host hydromagnetic dynamos, their magnetic fields are expected to interact with their HJs, especially if the magnetic field of the planets is of the order of a few Gauss or greater. Shkolnik et al. (2005, 2008) found evidence of a chromospheric hot spot rotating in phase with the orbital period of the HJ in HD 179949 and υ And during some observing seasons, while in other seasons the chromospheric flux showed only the usual rotational modulation. For instance, observing υ And in 2009, Poppenhaeger et al. (2011) found only a rotational modulation of the chromospheric flux. Absence of a planetary signature in the case of HD 179949 is reported by Scandariato et al. (2013) in the fall of 2009, while Gurdemir et al. (2012) claim the presence of a planetaryinduced hot spot in 2006. Therefore, the chromospheric flux enhancement associated with the planet is a nonstationary phenomenon, observed only in ≈30−40 percent of the seasons. The power irradiated by the hot spot falls in the range 10^{20}−10^{21} W (Shkolnik et al. 2005), i.e., a few percent of the total chromospheric flux.
Evidence of a corresponding enhancement in the corona has been provided by Saar et al. (2008) who observed a modulation of the Xray emission with the orbital phase in HD 179949. Kashyap et al. (2008) claimed that stars with massive planets closer than 0.15 AU have an Xray coronal luminosity ≈2 times greater than stars with distant planets (>1.5 AU). However, Poppenhaeger et al. (2010) and Poppenhaeger & Schmitt (2011) questioned such a result and explained the apparent correlation between HJs and enhanced Xray emission as a consequence of selection effects. They investigated the Xray spectral energy distribution of six planethosting stars finding that most of the emitting plasma was at temperatures of ≈1 MK, although two stars showed also a significant contribution from a hotter component at ≈4−5 MK. Shkolnik (2013) has investigated samples of stars with closein (a < 0.1 AU) and farout (0.5 < a < 2 AU) planets finding marginal evidence (i.e., at a 2σ level) of an enhanced far ultraviolet emission in the former.
Evidence of photospheric signatures associated with a possible magnetic starplanet interaction are reviewed by, e.g., Lanza (2011) considering also their possible connection with chromospheric and coronal signatures. Recently, Herrero et al. (2013) have discussed further observations supporting the phenomenon.
Theoretical models to account for chromospheric hot spots have been proposed by, e.g., McIvor et al. (2006), Preusse et al. (2006), and Lanza (2008). Those models do not address the origin of the energy dissipated in the hot spots limiting themselves to explaining the phase lag between the spots and the planets. The numerical simulations by Cohen et al. (2011a,b) indicate that a sufficient power can be released in the case of active stars, such as HD 189733, while Lanza (2009, 2012) suggests that energy can be released not only at the site where coronal and planetary fields reconnect, but also in the whole interconnecting loops that can be formed after reconnection (see Sect. 3.2 below).
The outer atmospheres of HJs can be probed through observations of transiting planets. Specifically, the depth of transit varies as a function of wavelength when an extended atmosphere absorbs the stellar flux. For HD 209458 and HD 189733, the absorption in the Lymanα line indicates an atmosphere extending beyond the Roche lobe and evaporating from the planet (e.g., VidalMadjar et al. 2003; Lecavelier Des Etangs et al. 2010). Recent observations by Lecavelier des Etangs et al. (2012) reveal a strong temporal variability of the evaporation rate in the case of the active star HD 189733, probably produced by a variation of the extreme ultraviolet (EUV) flux of the star. Atoms and ions of several chemical species, e.g., O, C, Si, Na, have also been observed (e.g., Sing et al. 2008; Linsky et al. 2010; VidalMadjar et al. 2011). The estimated mass loss rate ranges between 10^{6} and 10^{8} kg s^{1} and the velocity of the evaporation flow is of the order of ~10 km s^{1} with a temperature of ≈(1−2) × 10^{4} K close to the planet, although velocities up to ~100−150 km s^{1} and temperature of ≈10^{5} K have been inferred at several planetary radii (e.g., Lecavelier Des Etangs et al. 2008; Linsky et al. 2010; Lecavelier des Etangs et al. 2012). With these evaporation rates, a planet having a mass comparable with that of Jupiter can survive till the end of the mainsequence evolution of its host star. However, the situation can be different for less massive planets that experience an higher heating rate during the early mainsequence evolution of their host stars, especially in the first 0.5 − 1 Gyr when they are remarkably more active. Lecavelier Des Etangs (2007) has estimated the lifetime of planets of different mass finding that the observed planets are massive enough to survive evaporation, while no planet is observed in the range of mass that would imply a complete evaporation within a time interval of ~5 Gyr.
The situation can change if the planet has a magnetic field comparable with that of Jupiter or larger because the field can trap the ionized gas that is collisionally coupled to the neutral atoms thus preventing the evaporation of the atmosphere. Trammell et al. (2011) and Adams (2011) have investigated the effect of the planetary field on the evaporation considering a simple isothermal wind model. Assuming that the planet has a dipolelike field with the axis of symmetry coincident with the spin axis and normal to the orbital plane, they find that the field prevents evaporation in the lowlatitude region where the hot plasma is frozen to closed field lines. At high latitudes, there are open field lines, but the effect of the tidal force can still reduce or even stop the evaporation in extreme cases (Trammell et al. 2011). Neglecting the effects of tidal forces, the plasma can freely evaporate at high latitudes, i.e., for a colatitude θ ≤ θ_{0}, where θ_{0} ≃ 20° − 30° in the model of Adams (2011). In those models, the energy to accelerate the planetary wind comes from the EUV radiation emitted by the host star. However, the flux reaching the atmospheric layers from which the wind is launched is significantly reduced at high latitudes owing to the oblique trajectory of the light rays. The dilution is proportional to sinθ in comparison with normal incidence. Therefore, it is questionable whether a wind coming from the polar caps of the planet and powered by the stellar EUV flux can account for the observed mass loss rates (see Sect. 5 in Adams 2011).
In the present work, we investigate an additional source of energy to power the evaporation flow, that is nonthermal electrons accelerated at the reconnection site between the stellar and planetary fields. By analogy with solar flares, nonthermal electrons convey a significant fraction of the energy released in the reconnection process to the footpoints of the field lines in the polar regions of the planet. We consider the power released by magnetic reconnection and compare it with the EUV stellar flux in Sect. 2. In Sect. 3, we introduce a simple model for the electron acceleration, consider the dissipation of energy in the outer planetary atmosphere, and estimate the evaporation rate. In Sect. 4, we present an application to the sample of transiting planets and discuss the results in Sect. 5.
2. Importance of magneticpowered evaporation
The minimum power needed to support an evaporation rate of Ṁ_{p} can be estimated from the work done against the planet gravity as GM_{p}Ṁ_{p}/R_{p}, where G is the gravitation constant, M_{p} the mass of the planet, and R_{p} its radius. If an EUV flux F_{EUV} induces the evaporation, the available power is: (1)that allows us to estimate the maximum evaporation rate as (cf., e.g., Watson et al. 1981; MurrayClay et al. 2009): (2)by assuming that all the EUV energy flux is spent to accelerate the flow. If evaporation is powered by energy released by magnetic reconnection, the maximum is given by an analogous formula where the magnetic power P_{mag} appears instead of .
In principle, two different mechanisms can be invoked to produce a release of magnetic energy (e.g., Lanza 2009, 2012): a) reconnection between stellar and planetary magnetic fields at the magnetospheric boundary; b) relaxation of a stressed magnetic loop interconnecting stellar and planetary fields. Figure 1 provides a sketch of the magnetic field configurations associated with these mechanisms that will be further discussed in Sect. 3.1. The power made available by mechanism a) is remarkably smaller than that of b) (cf. Sect. 3.2). Therefore, to adopt a conservative approach in the present exploratory investigation, we consider only the power released by magnetic reconnection at the boundary of the planetary magnetosphere. The power released in the relaxation of an interconnecting loop will be considered in Sect. 3.7.
We assume that the stellar magnetic field B_{s} varies with distance r_{s} from the centre of the star as: (3)where B_{0} is the field at the pole of the star, R_{s} the star radius, and the exponent 2 ≤ s ≤ 3 depends on the geometry of the field lines. Specifically, for a radially directed field, as in the case of a strong stellar wind that blows the field lines nearly radially away from the star, s = 2, while for a potential dipole field s = 3. Values in between these two extrema are also possible, as in the case of the nonlinear forcefree coronal field discussed in Sect. 3.3 of Lanza (2012). The magnetic field of the planet is assumed to be dipolar with the axis perpendicular to the plane of the orbit.
The power P_{mag} released by magnetic reconnection is estimated according to the model in Sect. 2.2 of Lanza (2012): (4)where μ is the magnetic permeability of the plasma, B_{p} the planetary magnetic field at its pole, M_{s} ≫ M_{p} the mass of the star, and a the semimajor axis of the orbit assumed to be circular.
Fig. 1
Upper panel: a meridional section of the stellar and planetary magnetic field lines in a magnetohydrostatic regime in which the stellar field is forcefree and the planetary field is potential. In this case the flux systems of the two bodies are topologically separated (see Sect. 3.1). The field lines of the star and the planet can release energy by reconnecting at the boundary between their magnetospheres that is highlighted with a red dashed contour. Lower panel: the magnetic field of the starplanet system in a magnetohydrostatic regime in which the stellar field is potential. In this case, field lines interconnecting the two bodies exist. Interconnecting magnetic loops can be stressed by moving their footpoints and accumulate energy that can be successively released by a relaxation to the minimumenergy potential field configuration (see Sect. 3.1 for details). 
2.1. Stellar EUV fluxes
Unfortunately, photons with a wavelength shorter than 91.2 nm are strongly absorbed by neutral hydrogen in the interstellar medium, thus only very limited information is available on the fluxes of stars in the wavelength range 10 ≲ λ ≲ 90 nm that provides most of the energy for photoevaporation of planetary atmospheres. Lecavelier Des Etangs (2007) and other authors have extrapolated the EUV fluxes from the observed Xray fluxes, assuming the same EUVtoXray flux ratio of the Sun. An alternative is to build models of stellar coronal emission, calibrated with solar observations and the few stellar measurements presently available, to compute fluxes in the EUV domain (SanzForcada et al. 2011). Although the method of combining models and observations is likely to provide more accurate results, here we prefer to use the simpler approach of Lecavelier Des Etangs (2007). Specifically, we adopt the EUV fluxes given in his Sect. 3 as representative of the irradiation of stars of different spectral types in the middle of their mainsequence evolution.
2.2. Stellar and planetary magnetic fields
The surface magnetic fields of some planethosting stars have been measured using spectropolarimetric techniques. For Ftype stars, poloidal fields up to ~10 G have been found in the case of HD 179949 (Fares et al. 2012) and τ Bootis (Donati et al. 2008; Fares et al. 2009). For Gtype stars, we have a global poloidal field of 1 − 1.5 G in the Sun, while fields up to 3 − 5 G are observed in stars with rotation periods in the range 20−23 days. Stars with a rotation period of 9−12 days show fields of ≈40−50 G (Petit et al. 2008). In the active Ktype star HD 189733, that has a rotation period of ~12 days, Moutou et al. (2007) measured fields up to ~40 G. Comparable or stronger fields are present at the surface of Mtype dwarfs (see Donati & Landstreet 2009,for more information). In view of these observations and to be conservative, we adopt a stellar magnetic field intensity B_{0} = 10,5,2,10 G for the ranges of effective temperatures T_{eff} ≥ 6300 K, 6000 ≤ T_{eff} < 6300 K, 5300 ≤ T_{eff} < 6000 K, and T_{eff} < 5300 K, respectively. They correspond to the spectral type ranges F6F7, F8F9, G, and K plus M of Sect. 3 of Lecavelier Des Etangs (2007), respectively.
Planetary magnetic fields have not been measured yet, although a single radio flux measurement at 150 MHz suggests a field intensity of ≈50 G in HATP11b to be confirmed by further observations (Lecavelier des Etangs et al. 2013). Note that in the case of Jupiter, the field strength derived from its radio emission is ~14 G. An alternative way to estimate planetary fields is provided by the observations of magnetospheric bow shocks in transiting planets. Interpreting the observations of the transit of WASP12 in the UV domain with this model, Vidotto et al. (2010) estimate an upper limit to the planetary field of ~24 G. Theoretical models computed by Christensen et al. (2009) and Reiners & Christensen (2010) indicate that magnetic fields between ~10 and ~100 G should be typical of planets with a mass between ~1 and ~10 Jupiter masses, with lower values characteristic of older planets because the internal heat flux that powers their magnetohydrodynamic dynamos decreases with age. Batygin & Stevenson (2010) find similar fields by scaling according to the planet Elsasser number computed by assuming tidal locking between its rotation and orbital revolution (see also Rüdiger & Hollerbach 2004,for a similar approach in the case of solar system planets). Therefore, it seems reasonable to adopt B_{p} = 10 G for all the massive planets considered in the present investigation.
2.3. Comparison of the input powers
The ratio P_{mag}/P_{EUV}, as obtained from Eqs. (4) and (1), is plotted vs. the orbital semimajor axis in Fig. 2 for the transiting planets known on 25 March 2013 whose parameters have been extracted from the database in http://www.exoplanets.org. We conservatively assume that only 30 percent of the power given by Eq. (4) is converted into accelerated electrons and that only half of them reach the planetary atmosphere, i.e., only 1/6 of the reconnection power is available for evaporation (cf. Sect. 3.3). Note that the ratio P_{mag}/P_{EUV} does not depend on the radius of the planet. We consider three stellar magnetic field geometries as parametrized by s = 2,2.5, and 3. Considering the 200 planets within 0.1 AU, the magnetic reconnection power exceeds the EUV power in 81.5, 22.5, and 4 percent of the cases for s = 2,2.5, and 3, respectively. Considering the 115 planets with a ≤ 0.05 AU, the percentages rise to 92.1, 26.1, and 6.9, respectively. The four planets in the lower left corners of the plots are the three planets of Kepler42 and GJ 1214b that orbit lowmass stars for which the ratio a/R_{s} is significantly greater than in the case of the other closein planets thus making their P_{mag} remarkably small.
The mean values of the ratio P_{mag}/P_{EUV} are ~2.5, ~0.7, and ~0.2 for s = 2,2.5, and 3, respectively. Therefore, the magneticinduced evaporation should not significantly modify the predictions by Lecavelier Des Etangs (2007) on the timescale required for a complete evaporation of a given planet, except possibly for s = 2. However, in a planet endowed with a significant magnetic field, evaporation by EUV irradiation is prevented by the closed field lines in the lowlatitude region and only at high latitudes (say, ≳60°) the flow is allowed to leave the planet along open field lines (cf. Adams 2011; Trammell et al. 2011). Therefore, the EUVinduced evaporation is significantly less efficient than assumed by Lecavelier Des Etangs (2007), if planetary magnetic fields are not negligible. The magneticinduced evaporation may compensate for this reduction if s ≲ 2.5.
We conclude that the power released by magnetic starplanet interaction can play a relevant role in the evaporation of several closein planets and is therefore worth of a detailed investigation. A first step in this direction is undertaken in the next sections of this work.
Fig. 2
Upper panel: the ratio P_{mag}/P_{EUV} vs. the orbital semimajor axis a for the known transiting planets computed for a stellar radial magnetic field (s = 2). The value of P_{mag} given by Eq. (4) has been multiplied by 1/6 (see the text). The dotted horizontal line corresponds to P_{mag} = P_{EUV}. Middle panel: the same as the upper panel, but for s = 2.5 corresponding to a nonlinear forcefree stellar field. Lower panel: the same as the upper panel, but for s = 3, corresponding to a potential dipole stellar field. 
3. Model
3.1. The stellar corona
When the kinetic energy, gravity and thermal pressure p of the plasma are much smaller than the magnetic pressure B^{2}/2μ, the field rules the dynamics and the energy balance of the plasma. In this case, a magnetohydrostatic configuration is characterized by the vanishing of the Lorentz force, i.e., J × B = 0, where J = μ^{1}∇ × B is the current density. In other words, ∇ × B = αB, where the force freeparameter α is constant along magnetic field lines, as immediately follows by taking the divergence of both sides of the equation (see, e.g., Lanza 2012,for more details).
A potential magnetic field has α = 0 along all its field lines and represents the minimum energy state compatible with the boundary conditions. The stellar field is not in general in a potential state because the magnetic stresses at the stellar photosphere store energy in the field itself making it deviate from its minimumenergy state. The relaxation of the field to its potential state is generally prevented by the conservation of magnetic helicity (see, Lanza 2009, 2012,for details). On the other hand, the magnetic field of the planet is in a potential state (α = 0) because currents are generally not allowed to circulate in the nearly neutral planetary atmosphere. Therefore, in the considered magnetohydrostatic regime, a magnetic field line cannot connect the planet with the stellar surface, except when the stellar field is potential, because α is constant along field lines.
If the stellar field is nonpotential, the flux systems of the star and the planet are topologically separated and their magnetic field lines come into contact on the boundary of the planetary magnetosphere(cf. Fig. 1, upper panel). In general, the field lines belonging to the two systems are not parallel on the magnetospheric boundary and the orbital motion of the planet pushes them into contact on one side of the surface forcing continuous magnetic reconnection. Considering that the Alfvén velocity in the stellar corona is generally one order of magnitude larger than the orbital velocity of the planet and that the ratio β ≡ 2μp/B^{2} ≪ 1 (Lanza 2008, 2009, 2012), the magnetic pressures of the stellar and planetary fields are in equilibrium at each point of the boundary, i.e., , where B_{s}(r_{s}) is the stellar field at position r_{s} with respect to the barycentre of the star, and B_{p}(r_{p}) the planetary field at position r_{p} with respect to the barycentre of the planet.
When the magnetic helicity of the stellar field is negligible, the coronal field can relax to its minimumenergy potential state and in this case there are magnetic field lines interconnecting the surface of the star with the planet(cf. Fig. 1, lower panel). Close to the planet, we approximate the potential field as in Sect. 3 of Adams (2011) because the radii of the star and the planet are much smaller than their separation (cf. also Adams et al. 2011).
Considering a spherical polar reference frame (O,r_{s},θ,ϕ) with the origin O at the barycentre of the star and the polar axis ẑ along the stellar spin axis, a stationary flow of the coronal plasma is ruled by the equation: (5)where ρ is the density, v the velocity of the plasma, and the total gravitational plus centrifugal potential, with Ω the angular velocity of rotation of the star, and θ the colatitude measured from its North pole. In the regime β ≪ 1, the plasma flows along magnetic field lines, that is v = vŝ, where ŝ is a unit vector in the direction of the magnetic field, i.e., B = Bŝ. Since thermal conduction is very efficient along field lines at temperatures of the order of ~ 10^{6} K (Priest 1984), we assume that the coronal temperature T_{c} is constant along a given field line, i.e., ∂T_{c}/∂s = 0, and that the plasma follows the ideal gas law: , where is the gas constant and the mean molecular weight. Under these hypotheses, Eq. (5) can be integrated along a given field line by considering that (v·∇)v = ∇(v^{2}/2) + (∇ × v) × v, and gives: (6)with the potential Φ given by: (7)where ρ_{0} is the coronal density at the footpoint of the field line, the ratio of the centrifugal to the gravitational potential on the equator of the star, and Φ_{0} the potential at the footpoint of the field line on the stellar surface. Equation (6) assumes that the velocity vanishes at the footpoint of the field line.
3.2. Interaction between stellar and planetary fields
When the stellar field is in a forcefree nonpotential state, it interacts with the planetary field at the boundary of the planetary magnetosphere that is the surface of discontinuity between the two fields as discussed in Sect. 3.1. Lanza (2009, 2012) has computed the power released by magnetic reconnection at the magnetospheric boundary and its dependence on stellar and planetary parameters considering different models for the stellar coronal field. The maximum dissipated power is found to range from 10^{17} to 10^{19} W, even considering rather extreme values of the parameters. Such values are insufficient by at least one order of magnitude to account for the power radiated by the chromospheric hot spots observed by Shkolnik et al. (2005, 2008) and Gurdemir et al. (2012).
In the present approach, Alfvén waves have been filtered out by assuming that the Alfvén velocity is much larger than the orbital velocity of the planet. However, even without this assumption, it is not possible to account for powers of the order of 10^{20} − 10^{21} W, as found by Saur et al. (2013) by considering the energy fluxes in the socalled Alfvén wing model.
To account for greater powers, we consider a magnetic loop interconnecting the surface of the star with the planet that is steadily stressed by the orbital motion of the planet. The accumulation of energy drives the field out of its initial potential state as in the case of a solar magnetic arcade that is stressed by shear motions at its footpoints. At a given point, the arcade losses equilibrium and the field erupts producing a flare while the magnetic helicity is taken away by a coronal mass ejection allowing the field to relax again to a nearly potential state (e.g., Flyer et al. 2004; Zhang et al. 2006; Zhang & Flyer 2008). The same process is expected to occur in the case of an interconnecting starplanet loop. The effectiveness of the eruption process in taking away magnetic helicity is crucial in this case because, if helicity accumulates in the field, it will lead to a separation of the stellar and planetary flux systems, i.e., the disappearance of the interconnecting loop (cf. Sect. 3.1).
Assuming a sequence of cycles consisting of stress accumulations and eruptions, the average energy dissipation rate in the interconnecting loop must be equal to the average flux of the Poynting vector across its base. In other words, the mean available power is: (8)where f_{AP} is the fraction of the planetary hemisphere crossed by the field lines of the interconnecting loop, E = −v_{0} × B_{p} the electric field, v_{0} the relative orbital velocity, and B_{p} the magnetic field at the base of the loop, located on the surface of the planet. In terms of the surface field of the planet, we have: (9)The fraction of the planetary surface magnetically connected to the stellar field is given by Eq. (27) of Adams (2011): (10)where ζ ≡ B_{s}(a)/B_{p} is the ratio of the magnetic field of the star to that of the planet at the distance a of the planet (cf. Eq. (3)). Note that our parameter ζ is denoted as β in Adams (2011), the symbol being changed here to avoid confusion with the plasma β. Equation (9) is valid when the timescale for the field relaxation to the potential state τ_{rel} is comparable with or shorter than the timescale for magnetic energy buildup by the stress produced by the orbital motion of the planet τ_{acc}. The typical value of τ_{rel} is comparable with the Alfvén transit time along the interconnecting loop and is of the order of ~10^{3} s for closein planets (cf., e.g., the final part of Sects. 2.3 and 3.1 of Lanza 2012), while τ_{acc} ~ 2R_{p}/v_{0} is the timescale for the orbiting planet to cross the base of the interconnecting loop and is of the order of 10^{3} − 10^{4} s for typical R_{p} ~ 10^{8} m and v_{0} ~ 10^{4} − 10^{5} m s^{1}, thus satisfying the requirement τ_{rel} ≤ τ_{acc}.
For a field strength B_{0} ~ B_{p} ~ 10 G, a relative velocity v_{0} ~ 10^{4} − 10^{5} m s^{1}, a planetary radius R_{p} ~ 10^{8} m, and f_{AP} ~ 0.1 − 0.2, we find P ~ 10^{20} − 10^{21} W, that is sufficient to explain the power radiated by chromospheric hot spots associated with HJs. That given by Eq. (9) is indeed the maximum available power that one can expect from magnetic starplanet interaction given the hypotheses adopted here.The energy released in the interconnecting loop is ultimately provided by the orbital motion of the planet, but its dynamical effect on the orbit is negligible in comparison with, e.g., tidal interactions because such a kind of loop is present only for a limited fraction of the time (≈ 10 − 20 percent; cf., e.g., Sect. 4.2 of Lanza 2011).
The constant phase lag between the planet and the hot spot (Shkolnik et al. 2005) can be explained by assuming that the coronal field is initially in an axisymmetric twisted forcefree regime (Lanza 2008, 2009) before reconnecting with the planetary field. During the loss of equilibrium occurring in the reconnection phase, the interconnecting field gets rid of its magnetic helicity by erupting twisted flux tubes as in solar coronal mass ejections, thus relaxing to a nearly potential configuration. The footpoints of the interconnecting loop on the stellar surface cannot move because the field is frozen to the photospheric plasma. In this way, the subsequent stress is applied to a nearly potential loop with the initial phase lag between its stellar and planetary footpoints almost unchanged. If the planet crosses a sequence of loops of an unperturbed forcefree axisymmetric field with an interaction timescale short in comparison with its orbital period, an approximately constant phase lag between the planet and the chromospheric hot spot is expected. If the stellar field has a significant nonaxisymmetric component, the time dependence of the energy release becomes remarkably complicated and neither a constant phase lag nor a simple dependence of the irradiated power on the orbital phase is expected (cf., Lanza 2012,for details).
The energy dissipation in a stressed interconnecting loop is not localized at the boundary of the planetary magnetosphere, as in the case of the interaction of the separated starplanet flux systems. Applying a largescale shearing motion to the footpoints of a solar coronal loop whose magnetic field lines are steadily intermixed by turbulent photospheric motions, Janse & Low (2009, 2010) find that a spatially dense distribution of current sheets is produced extending over a large fraction of the loop volume. Given the similarity with our case, we conclude that the dissipation of magnetic energy occurs over most of the volume of the interconnecting loop. The loop sections having the largest field intensities, i.e., those close to the stellar and the planetary surfaces, are expected to have the highest energy release rates per unit volume.
3.3. Electron acceleration
Since the energy released by magnetic reconnection is mainly transported into the planetary atmosphere by accelerated electrons, we now consider their acceleration process in some detail. Electron acceleration is characteristic of solar flares and reveals itself through the emission of pulses of hard Xrays from the footpoints of flaring loops where the particles impinge upon the dense chromospheric plasma and emit hard Xray bremsstrahlung radiation. The acceleration occurs at the reconnection sites inside flaring loops or magnetic arcades generally located up in the solar corona (see, e.g., Aschwanden 2002, for a review). In a large flare, the total number of electrons accelerated to energies of 20 keV or higher is of the order of 10^{38} − 10^{39} and the timescale of acceleration is as short as ≈300 ms for an individual pulse (Miller et al. 1996). The total energy converted into accelerated electrons is a significant fraction of the total flare energy generally ranging from ≈30 to 60 percent (e. g., Holman 2005; Mann & Warmuth 2011). On the other hand, the energy of other particles (protons, ions) is generally lower or comparable with that of the electrons (see Aschwanden 2002, Sect. 7.4), therefore we shall consider only electrons for simplicity.
Several different acceleration mechanisms have been proposed to occur in the reconnection sites of solar flares (cf., e.g., Miller et al. 1996; Aschwanden 2002; Petrosian 2012). Here we adopt a simple heuristic approach to relate the spectrum of the accelerated electrons to the average magnetic field and particle density in the reconnection region, inspired by the discussion in Sect. 4.1.1 of Aschwanden (2002). The density spectrum of the accelerated electrons is denoted by , where is the number of particles with kinetic energy between E_{0} and E_{0} + dE_{0} per unit volume. We assume that the spectrum has the form: (11)where E_{c} is the socalled cutoff energy, δ the index of the powerlaw distribution for energies beyond the cutoff, and A a normalization constant (see below).A positive and finite value of the total electron energy per unit volume requires δ > 2 (cf. Eq. (13)). The powerlaw dependence of the energy spectrum follows from the spectrum of hard Xrays observed during the impulsive phase of solar flares (see, e.g., Holman et al. 2011). The energy cutoff E_{c} is generally estimated to fall between 15 and 30 keV, but a unbiased determination of its range is difficult because the thermal Xray component emitted by the very hot plasma produced during a flare hides the change in slope of the Xray bremsstrahlung spectrum when the cutoff is at low energies. The powerlaw index δ usually ranges between 4 and 6. The energy distribution below the cutoff is assumed to be flat because a sharp or peaked cutoff or a positive slope would lead to instabilities that rapidly reestablish a flat distribution, as discussed in Sect. 6 of Holman et al. (2011).
If we denote by n_{T} the number density of the accelerated electrons, we have: (12)yielding A = (n_{T}/E_{c})(δ − 1)/δ. The total energy per unit volume is: (13)The source of energy for electron acceleration is the coronal magnetic field. On average, some fraction f of the magnetic energy per unit volume B^{2}/2μ is converted into n_{T} accelerated electrons with a total energy E_{T}. If g is the fraction of accelerated electrons, i.e., n_{T} = gn_{e}, where n_{e} is the electron number density, we have: (14)Assuming in a typical solar flare B ~ 100 G, n_{e} ~ 10^{16} m^{3} (Miller et al. 1996), E_{c} ~ 25 keV, and f ~ 0.3, we obtain from Eq. (14), g ~ 0.1 − 0.3, that is a significant fraction of the coronal plasma electrons is accelerated to suprathermal energies in a solar flare. Conversely, when the fraction of accelerated electrons g, the energy conversion fraction f, the powerlaw index δ, the electron density n_{e}, and the field strength B are given, we can use Eq. (14) to estimate the cutoff energy E_{c}. This approach will be applied in Sects. 3.4 and 4 to estimate E_{c} for the electrons accelerated at the boundary of the planetary magnetosphere. The magnetic field intensity B is given by Eq. (3), while the number density n_{e} follows from Eq. (6) assuming that the plasma is fully ionized.
3.4. Energy dissipation in the outer planetary atmosphere
After their acceleration, highenergy electrons travel along magnetic field lines reaching the stellar chromosphere at one footpoint and the planetary atmosphere at the other footpoint. The energization of the stellar atmosphere has been considered by, e.g., Gu & Suzuki (2009), thus we focus on the energy dissipated into the planetary atmosphere. We follow Holman et al. (2011) to determine the energy lost by the electrons during their travel through the atmospheric layers. The atmosphere is considered to be isothermal with a temperature T (see Sect. 3.5 for justification) and the tidal gravitational potential of the star is neglected. For simplicity, we assume a constant ionization fraction x ≡ n_{p}/n, where n is the total number density of protons, including both free protons of density n_{p} and those bound into hydrogen atoms of density n_{H}, i.e., n = n_{p} + n_{H}. The layers from which the planetary evaporation starts have a typical density n ≈ 10^{15} − 10^{17} m^{3} and a temperature of ≈10^{4} K (MurrayClay et al. 2009; Trammell et al. 2011; Adams 2011), giving a pressure of 10^{9} − 10^{7} bar. In those layers, the chemical composition is dominated by hydrogen (cf., e.g., Venot et al. 2012). Therefore, it is justified to assume an atmosphere of pure hydrogen. Our model applies to the highlatitude zones of the planet where the footpoints of the interacting magnetic field lines are located and the magnetic field can be considered to be vertically directed. The speed of the evaporation flow is generally smaller than the sound speed in the layers where most of the energy is deposited, therefore the stratification is not remarkably different from hydrostatic. In other words, we assume that the number density of protons in the atmosphere varies with the height z from its base layer as: n(z) = n_{p0}exp( − z/H), where n_{p0} is the base density and is the atmosphere scale height, with the gravitational acceleration at the surface of the planet and the mean molecular weight for a pure hydrogen partiallyionized atmosphere .
The loss of energy dE experienced by a single electron of energy E crossing a layer between the heights z and z + dz is (see Brown 1973; Holman et al. 2011,Sect. 4.2): (15)where K′ = 2πc^{4}e^{4}Λ, with c the speed of light, e the electron charge, Λ ≡ Λ_{ee} − Λ_{eH}, λ ≡ Λ_{eH}/Λ, with Λ_{ee} and Λ_{eH} the Coulomb logarithms for electronelectron and electronhydrogen atom collisions, respectively. Since Λ_{ee} ≃ 20 and Λ_{eH} ≃ 7.1 are almost constant in the temperature and density regime of interest, Λ ≃ 12.9 and λ ≃ 0.55. Therefore, the constant K′ ≃ 1.87 × 10^{22}(Λ/12.9) keV^{2} m^{2}. Assuming that the ionization fraction x is constant, Eq. (15) can be integrated to give the energy of an electron at height z: (16)where E_{0} is the energy of the particle when it leaves the acceleration site, K′′ ≡ (λ + x)K′, and the column density of the atmosphere down to the level z. If ℱ(E_{0})dE_{0}, is the number of electrons with initial energy between E_{0} and E_{0} + dE_{0} that goes through the unit area on top of the atmosphere per unit time (electron flux density), the total power dissipated in the layer between z and z + dz per unit area is: (17)where is the minimum energy required for an electron to reach the level z without being stopped by the energy losses in the upper layers of the atmosphere, and Q_{c} the energy loss per proton. Assuming that there is no significant energy loss for the electrons during the travel from their acceleration site to the limit of the atmosphere (an assumption supported by solar flare observations; cf., e.g., Aschwanden 2002,Sect. 6), the flux density spectrum at the limit of the atmosphere is , where is the velocity of the electrons of energy E_{0} and mass m_{e}. From Eq. (11), we obtain: (18)On the other hand, the energy flux density of accelerated electrons at the limit of the atmosphere is: E_{0}ℱ_{0}(E_{0}). Denoting by F_{0} the total energy flux carried by accelerated electrons (energy reaching the top of the atmosphere per unit surface and time) and integrating over the energy E_{0}, we find: (19)that gives the number density of accelerated electrons n_{T} from F_{0}, the cutoff energy E_{c}, and the slope of the spectrum δ.
Substituting Eq. (18) into Eq. (17), we can compute the energy dissipated by accelerated electrons into a given layer of the atmosphere. We express the integrand in terms of the nondimensional quantities w ≡ E_{0}/E_{c} and α_{e} ≡ E_{s}(z)/E_{c} and integrate by parts to eliminate the numerical divergence of the integrand for E = E_{s}, obtaining: (20)where . These expressions can be integrated numerically by truncating the upper limit of the improper integrals to some value w_{max} ≫ max(1,α_{e}) since the error introduced by this approximation does not exceed (w_{max})^{− δ}/δ.
When s = 2, in the layer where the planetary evaporation flow is launched α_{e} ≪ 1 in several cases. In this limit, we develop the square root factors in the above integrals by Newton’s binomial expansion and integrate the corresponding power series. Limiting ourselves to the terms up to the second order in α_{e}, we find: (21)For δ = 4 − 6 and α_{e} ≪ 1, the dissipation rate is dominated by the first two terms, thus Q_{c} is approximately constant because has a weak dependence on the column density in this limit.
The above treatment is adequate as far as the energy of most of the electrons is nonrelativistic, i.e., E_{c} ≲ 0.2m_{e}c^{2} ~ 100 keV, where m_{e}c^{2} is the electron rest energy (e.g., Emslie 1978; Leach & Petrosian 1981). However, when s = 2 in ~7 percent of the cases E_{c} exceeds the above limit requiring a relativistic treatment of the collisional energy loss and ionization rates (cf. Sect. 4). This fraction increases if we consider parameters for the stellar corona that lead to a lower particle density or a stronger magnetic field in the reconnection region than assumed in Sect. 4. Therefore, we treat the case of relativistic electrons in Appendix A.
3.5. Ionization and thermal balance of the outer atmosphere
We now consider the response of the outer planetary atmosphere to the accelerated electrons and the stellar ultraviolet radiation. The ionization fraction x = n_{p}/n is given in equilibrium by a balance between ionization (both radiative and collisional) and recombination processes. The plasma is electrically neutral, so the number densities of the electrons and protons are equal, i.e., n_{e} = n_{p}. Adopting the approximations introduced in Sect. 8 of Trammell et al. (2011), the ionization balance can be written as: (22)where α_{R} ≃ 2.6 × 10^{19}(T/10^{4}K)^{0.8} m^{3} s^{1} is the recombination rate and the photoionization rate per hydrogen atom, with N_{cH} = (1 − x)N_{c} the atomic hydrogen column density from the given point to the star, ν the frequency of the ionizing photon, f(ν) the photon flux per unit frequency interval, simply scaled from that of the Sun (Claire et al. 2012), σ_{pi}(ν) = σ_{0}(ν/ν_{0})^{3} the radiative ionization cross section of the hydrogen, hν_{0} = E_{i} = 13.6 eV the ionization potential, and σ_{0} = 6.3 × 10^{22} m^{2} the cross section at the ionization threshold; , is the ionization rate per unit hydrogen atom by electron collisions, where a_{0} = 5.3 × 10^{11} m is the Bohr radius, the ratio of the collisional ionization crosssection to the area corresponding to the Bohr radius, the number density of the incident electrons of energy E, v(E) their velocity, and their flux density. It can be computed from the flux density at the limit of the atmosphere as given by Eq. (18) taking into account the energy losses as specified by Eq. (16), i.e.: (23)Equation (23) follows from the continuity of the electron flux through the atmosphere, i.e., , and Eq. (16).
The dependence of the nondimensional crosssection on the energy E of the colliding electrons is given by Rudge (1968) and Dere (2007). In our case, the energy of the electrons is generally much greater than E_{i} and the collisional cross section scales as for E ≳ 10 E_{i}; for E_{i} ≤ E ≲ 10 E_{i}, we assume C equal to its mean value, i.e., ~0.6. Note that for E/E_{i} ~ 25, , therefore the lowenergy part of the electron spectrum produces most of the collisional ionizations.
Once the ionization fraction has been calculated in a given layer as a function of the incident EUV flux f(ν) and the accelerated electron flux , we can consider the energy balance. Since the atmosphere is approximately isothermal (see below) and the gas can move only along magnetic field lines, owing to the collisional coupling between ions, electrons, and neutral atoms, thermal conduction and convection can be neglected. The thermal balance in equilibrium is established between the heating sources (EUV flux and accelerated electrons) and the Lymanα radiation that is the most important contributor to the radiative losses in the considered range of temperature, i.e., T ~ (1 − 2) × 10^{4} K, for an optically thin plasma (cf. Dere et al. 2009, and references therein). Lymanα emission is excited by collisions of thermal electrons with hydrogen atoms. Assuming that all the energy of the photoelectrons produced by photoionization is converted into heat, the heating rate by the EUV flux per hydrogen atom is . The heating per unit volume by accelerated electrons is given by Eqs. (17) and (20). Therefore, the thermal balance equation is: (24)where L_{α}(T) = 2.9 × 10^{32}(T/10^{4}K)^{0.5}exp( − 1.184 × 10^{5}K/T) J m^{3} s^{1} is the radiative loss function of the Lymanα line (Dalgarno & McCray 1972), n_{H} = (1 − x)n, and n_{e} = n_{p} = xn. This expression for L_{α} is valid for T ≲ 1.2 × 10^{4} K. Note the exponential dependence on temperature that makes L_{α} increase rapidly with increasing temperature. In other words, a remarkable increase of the heating can be balanced by a modest increase of the temperature of the emitting plasma. On the other hand, an increase of the temperature increases the ionization, thus reducing the Lymanα emissivity for T > (1 − 2) × 10^{4} K. The interplay between these effects makes the Lymanα radiation a sort of thermostat that keeps the temperature of the upper atmospheric layers approximately constant at ≈10^{4} K for a wide range of heating rates (cf., e.g., Fig. 3 in Cox & Tucker 1969).
3.6. A simple isothermal wind model
An illustrative application of the theory introduced above is provided by a model of a transonic isothermal wind, i.e., with a uniform temperature T and ionization fraction x from the launching layer of the flow to the sonic point and beyond. Adams (2011) has introduced such a model that we use in this Section.
The layer from which the wind is accelerated corresponds approximately to the layer where the optical depth for the absorption of EUV photons by hydrogen atoms is unity. As a matter of fact, a fully consistent model should include the energization of deeper layers by nonthermal electrons. This would imply the solution of the full energy equation in an optically thick plasma that makes the problem much more complicated. Therefore, we limit ourselves to the assumption that the base of the planetary wind flow is at the level z_{0} where the optical depth (25)where ⟨ σ_{EUV} ⟩ = 2.0 × 10^{22} m^{2} is the mean absorption crosssection for the EUV radiation. In such a way, we can apply Eq. (24), valid for an optically thin plasma at Lymanα wavelength. The column density down to the level z in an isothermal, partially ionized atmosphere is: N_{c}(z) = n_{p0}Hexp( − z/H), where n_{p0} is the density at the base of the atmosphere (level z = 0). We assume n_{p0} = 2 × 10^{16} m^{3} and verify that our results are not significantly dependent on that value.
The flux of EUV radiation producing the heating and ionization of the planetary atmosphere has been specified in Sect. 2.1 according to Sect. 3 of Lecavelier Des Etangs (2007). The energy flux F_{0} of nonthermal electrons at the boundary of the atmosphere, required to compute the electron spectrum with the theory of Sects. 3.3 and 3.4, is given by: (26)where the power P_{mag} released in the starplanet interaction is given by Eq. (4) and f_{AP} is the fraction of the planetary surface magnetically connected to the star (cf. Eq. 10). In Eq. (26) a factor 1/6 is included to account for an efficiency of ~30 percent in electron acceleration and to consider that only half of the electrons stream toward the planet (cf. Sect. 2.3). To estimate the cutoff energy E_{c}, we apply Eq. (14), with f = 0.3, g = 0.1, and δ = 6, while Eq. (19) is applied to find n_{T} from F_{0}, thus completely specifying the electron energy spectrum.
For given EUV and nonthermal electron fluxes, we solve iteratively for the density, the temperature, and the ionization fraction at the base level z_{0}. Specifically, we start by guessing T ~ 10^{4} K and x ~ 0.5 at z_{0} and compute a first guess for z_{0} from Eq. (25); then we compute the density at z = z_{0} from n(z) = n_{p0}exp( − z/H), solve Eq. (22) to find a refined value of the ionization fraction, and Eq. (24) for a refined value of the temperature (from the Lymanα radiative loss function) that are used to estimate a better value of z_{0} by Eq. (25). By iterating this procedure, we converge to consistent values for z_{0}, n(z_{0}), T(z_{0}), and x(z_{0}) giving τ(z_{0}) = 1. Considering the parameters of the transiting exoplanets as listed in the website www.exoplanets.org on 25 March 2013, and excluding those with a mass below 0.05 Jupiter masses and a semimajor axis a > 0.15 AU, we find that the above procedure fails to converge only in ~1 percent of the cases.
Once the physical parameters at the base of the planetary wind are known, the mass loss rate is obtained by multiplying the nondimensional mass loss rate ṁ given by Eq. (64) of Adams (2011) by m_{proton}n(z_{0})R^{2}(z_{0})a_{s}(z_{0}), where m_{proton} = 1.67 × 10^{27} kg is the mass of the proton, R(z_{0}) the distance of the base of the wind from the centre of the planet, and a_{s}(z_{0}) the sound speed at level z_{0}. For convenience, Eq. (64) of Adams (2011) is reproduced here as: (27)where is a nondimensional measure of the depth of the planet gravitational well and the equation is valid for ζ ≲ 0.01. Equation (27) assumes that the gravitational field of the planet effectively limits the evaporation rate as it is the case for a Jovian mass body. However, when less massive planets are considered, the parameter b becomes lower because the planet radius changes only slightly, thus reducing the energy required to accelerate the flow and moving the sonic point closer to the planet. When b = 3 the sonic point along the polar streamline reaches the surface of the planet (cf. Eq. (50) of Adams 2011) and Eq. (27) is no longer valid. In this regime, the wind can escape freely along field lines with an initial speed comparable with the sound speed (cf. the similar case of the free photoevaporation of circumstellar discs in Störzer & Hollenbach 1999; Adams et al. 2004). In conclusion, we assume that Eq. (27) is valid for b > 3, while we assume that ṁ = 6.45 ζ^{1/3} (i.e., is constant at its maximum value reached at b = 3) when b ≤ 3 to avoid a discontinuity in the variation of the mass loss rate as a function of the parameter b.
Our simple wind model is based on an approximate treatment of the energy balance through the flow. Therefore, it comes of no surprise that it satisfies energy conservation only in order of magnitude. Specifically, we find that in ≈20 percent of the cases, the mass loss rate exceeds its energylimited value: (28)where F_{tot} = F_{EUV} + F_{0} is the total energy input to the planetary atmosphere per unit surface and time as given by the sum of the EUV and nonthermal electron fluxes (cf., e.g., Eq. (14) of Lecavelier Des Etangs 2007). The discrepancy is always smaller than a factor of 3 − 4, so we can trust our model in order of magnitude and apply it to compare cases with and without the contribution of nonthermal electron energization.
3.7. Evaporation rate in the case of an interconnecting loop
In Sect. 3.2, the power released in a stressed magnetic loop interconnecting the planet with the star was estimated and values about 2 − 3 orders of magnitude larger than in the case of reconnection at the boundary of the planetary magnetosphere were obtained. With such a strong energization, the response of the planetary atmosphere cannot be predicted by our simple isothermal wind model in view of the limitations discussed in Sect. 3.6. Moreover, since the release of magnetic energy is probably distributed all along the loop with the highest powers where the field is stronger (cf. Sect. 3.2), a prediction of the electron energy spectrum is rather complicated. On the other hand, a simple energylimited model seems to be more appropriate to estimate the evaporation rate for a preliminary investigation of the effects on planetary atmospheres. Thus we apply Eq. (2) to estimate the mass loss rate in the case of a stressed interconnecting loop. In this case, the energy flux is the sum of the EUV flux and of the (much greater) energy flux coming from magnetic dissipation that falls on the fraction f_{AP} of the planetary surface as given by Eq. (10).
4. Results
In Fig. 3 we plot the cutoff energy E_{c} of the spectrum of nonthermal electrons accelerated at the boundary of the planetary magnetosphere. The cutoff energy is computed with the model in Sect. 3.3 assuming a stellar coronal base density of 5 × 10^{14} m^{3}, a coronal temperature T_{c} = 10^{6} K (cf., e.g., Poppenhaeger et al. 2010), and a spectral index δ = 6; the stellar field decays with distance as that of a dipole (s = 3 in Eq. (3)). The accumulation of the plotted points along discrete sequences is due to the discrete values adopted for the stellar magnetic field as indicated by different symbols (cf. Sect. 2.2). The dependence of E_{c} on the mass of the planet is quite weak and is not plotted.
Closein planets are characterized by values of E_{c} that reach up to ~10 keV for s = 3. The maximum cutoff energy reached in the case of a radial field dependence with s = 2 in Eq. (3) is ~120 keV. This allows the penetration of electrons into the planetary atmosphere down to column densities ranging from 10^{23} to 10^{25} m^{2} in the case of the most energetic spectra for s = 3 and s = 2, respectively. For s = 3, the penetration depth is comparable with the level where most of the EUV photons are absorbed corresponding to a column density of ~10^{22} − 10^{23} m^{2}. Only for the closest planets in the case with s = 2, E_{c} ~ 100 keV and the layer where most of the electron energy is deposited is remarkably deeper than that where EUV photons release most of their energy. In this regime, the parameter α_{e} ≡ E_{s}/E_{c} ≪ 1 in the layer where the EUV photons are absorbed and the planetary wind is launched, thus allowing us to apply Eq. (21) that gives an input energy rate approximately proportional to the number density n(z).
Fig. 3
Cutoff energy of nonthermal electron spectrum vs. orbital semimajor axis for our sample of transiting planets for a spectrum with an index δ = 6 and a stellar field decaying with distance with s = 3 (cf. Eq. (3)). Solid dots refer to a stellar field B_{0} ≥ 10 G, while open diamonds indicate B_{0} < 10 G (cf. Sect. 2.2). 
The ratio of the mass loss rate produced by nonthermal electrons and EUV flux to that induced by EUV flux only is plotted in Fig. 4 vs. the orbital semimajor axis for different values of the parameter s determining the radial dependence of the stellar field in Eq. (3). We plot only transiting planets with semimajor axis a < 0.15 AU, mass M_{p} > 0.05 M_{J}, with M_{J} being the mass of Jupiter, and ζ ≤ 0.01. In the cases in which Ṁ_{p} exceeds the energylimited value given by Eq. (28), the value of the mass loss rate has been fixed at that value (cf. Sect. 3.6).
Figure 4 can be compared with Fig. 2 that showed the ratio of the corresponding input powers. Taking into account the limitations of our simplified model for planetary evaporation, we see that in general the ratio of the mass loss rates is comparable with or is a few times larger than that of the corresponding input powers. Magneticallyinduced effects are larger for closein planets and for smaller values of s, as expected on the basis of Fig. 2, and become small for a ≳ 0.1 AU. The mean increase in the mass loss rate is by a factor of ~2.2 for s = 2, by ~ 2.1, for s = 2.5, and by ~1.3 for s = 3, with a maximum relative increase by a factor of 30 − 50 in the most extreme cases corresponding to closerin planets with mass ≳1.5 M_{J} (see below).
Fig. 4
Upper panel: ratio of the mass loss rate induced by nonthermal electrons and EUV photons to that induced by EUV photons only vs. the semimajor axis of the orbit for our sample of transiting planets (see the text). The parameter s = 2.0 specifies the radial dependence of the stellar magnetic field according to Eq. (3). Middle panel: the same as the upper panel, but for s = 2.5. Lower panel: the same as the upper panel, but for s = 3.0. 
Fig. 5
Upper panel: mass loss rate for a sample of transiting planets vs. their mass (see the text). Open diamonds indicate the loss rate induced by EUV flux only, while filled dots indicate the case when the contribution by nonthermal electrons accelerated at the boundary of the planetary magnetosphere is added. The parameter s = 2.0 specifies the radial dependence of the stellar magnetic field according to Eq. (3). Middle panel: the same as the upper panel, but for s = 2.5. Lower panel: the same as the upper panel, but for s = 3.0. 
We plot in Fig. 5 the mass loss rate Ṁ_{p} vs. the planet mass M_{p} with EUV energization only (open diamonds) and with both EUV and nonthermal electron energization (solid dots) for different values of s, respectively. The overall dependence of Ṁ_{p} on M_{p} is a consequence of the proportionality of ṁ to b^{3}exp( − b) in Eq. (27) and the small variation in planet radius and sound speed at the base of the wind that makes b approximately proportional to M_{p}. A similar dependence is apparent in Fig. 9 of Adams (2011), although his mass loss rates are more than one order of magnitude smaller than ours because of the systematically smaller EUV fluxes adopted by him and the assumption that only ≈0.5 of the EUV flux goes into wind acceleration. Our values seem to be in better agreement with the observations of HD 209458b and HD 189733b for which Ṁ_{p} ≈ 10^{7} and ≈10^{8} kg s^{1} have been reported, respectively (see Sect. 1).
In none of the considered cases the magnetic increase of the mass loss rate has a dramatic effect on the lifetime of the planets as estimated by τ_{ev} = M_{p}/Ṁ_{p}. It is plotted in Fig. 6 vs. the planet mass. Assuming a constant mass loss rate, the cumulative effect after ~5 Gyr amounts to a total mass loss of at most ≈10 percent.
The nonthermal electron energy input can significantly lift the level corresponding to unity optical depth at UV wavelengths. This affects the depth of the transit as observed in UV lines. Although a detailed comparison with the observations is outside the scope of the present investigation, we plot the relative variation of the radius corresponding to τ = 1 at the Lymanα wavelength vs. the planet mass in Fig. 7. For the systems with small magnetic effects, the variation is negligible, but for planets subject to a significant additional heating and having mass M_{p} ≲ 0.5 M_{J}, the increase can reach up to 30 − 40 percent, implying a significantly deeper transit. Therefore, the additional effects induced by starplanet magnetic interaction can be important for modelling transit observations in the UV.
All the above estimates on the relative importance of nonthermal electron vs. EUV energy fluxes are conservative because we have considered an efficiency of ~30 percent in the conversion of magnetic energy into accelerated electrons with only half of the electron flux impinging upon the planet (cf. Eq. (26)). On the other hand, we do not apply any dilution factor to the EUV flux, although the acceleration of the wind happens only in the highlatitude regions of the planet. Depending on the parameter ζ in Eq. (10), the effective EUV flux can be reduced by a factor of 2−5.
Fig. 6
Lifetime for the considered sample of transiting planets vs. their mass. The lifetime is estimated considering the maximum mass loss rate, i.e., by including the contribution of nonthermal electrons accelerated by reconnection at the boundary of the planetary magnetosphere. 
Fig. 7
Relative variation of the radius corresponding to optical depth τ = 1 in the centre of the Lymanα line for normal incidence across the planetary atmosphere vs. the planet mass. The plotted case corresponds to the slowest decay of the stellar field with distance (s = 2.0) that provides an upper limit to the predicted variation. 
Finally, we consider the possibility of a strong atmospheric heating in those systems in which there is a magnetic loop interconnecting the star with the planet. With the energylimited approach of Sect. 3.7, the average mass loss rate and the lifetime of the planets vs. mass are shown in Fig. 8 that has been computed assuming a constant planetary field B_{p} = 10 G and a constant stellar field B_{0} = 2 G with a radial dependence with s = 2.5. The timeaveraged mass loss rate is about one order of magnitude higher than in the case of magnetic reconnection localized at the boundary of the planetary magnetosphere. The instantaneous mass loss rate can be one order of magnitude greater than the average plotted in Fig. 8, computed by assuming that an interconnecting loop is present only ~10 percent of the time. With the adopted parameters, the lifetime of the transiting planets in our sample is always longer than 5 Gyr, indicating that a more frequent or stronger interaction is unlikely.
Fig. 8
Upper panel: mass loss rate averaged over time (in kg s^{1}) vs. planet mass in the case of the dissipation of magnetic energy in a loop interconnecting the star with the planet. Lower panel: the estimated lifetime of evaporating planets vs. their mass. The horizontal dashed line marks an evaporation lifetime of 5 Gyr corresponding to the mean estimated age of solarlike mainsequence stars. 
5. Conclusions
We have investigated starplanet magnetic interaction in closein planets and estimated its impact on the evaporation of their atmospheres. We have considered the energy provided by reconnection between planetary and stellar coronal fields, finding that it is comparable with or larger than the energy input due to stellar EUV radiation for closein (a ≲ 0.10 AU) planets. A simple model to derive the energy spectrum of the electrons accelerated at the boundary of the planetary magnetosphere is introduced and the dissipation of their energy into the planet atmosphere is treated in detail. In most closein planets, accelerated electrons can reach atmospheric layers down to column densities of 10^{23}−10^{25} m^{2}, i.e., comparable with or significantly deeper than the layers where EUV photons are absorbed, corresponding to column densities of 10^{22}−10^{23} m^{2}. The heating by accelerated electrons is not uniformly distributed over the planetary surface, but is localized close to the magnetic poles at the footpoints of magnetic fields lines, i.e., over an area of ≈10−20 percent of the whole surface (Adams 2011).
The theory introduced above can be applied to model evaporation or chemical reactions in the outer layers of a planetary atmosphere. Here, we limit ourselves to a simple application to the isothermal transonic wind model of Adams (2011) and Trammell et al. (2011). The magnetic contribution is found to be important in closeby (a ≲ 0.10 AU), massive (M_{p} ≳ 1.5 M_{J}) planets for which the evaporation rate can be increased up to a factor of 30 − 50 in extreme cases. This is due to the additional heating of the atmosphere by the accelerated electrons that increases the temperature and the sound speed thus reducing the value of the parameter b appearing in Eq. (27). In turn, the exponential dependence on b produces a remarkable increase of the mass loss rate with the additional heating. Moreover, this additional atmospheric heating produces an increase of the apparent radius of the planet at EUV wavelengths that is relevant, e.g., to modelling the transit profile in the Lymanα line.
In addition to the energy released by magnetic reconnection at the boundary of the planetary magnetosphere, we consider energy dissipation in a loop interconnecting the stellar and planetary fields. Such loops can be invoked to account for the hot chromospheric spots observed in some seasons in a few systems hosting hot Jupiters (Shkolnik et al. 2005, 2008). The average mass loss rates predicted in that case reach up to ≈(0.3 − 1.0) × 10^{10} kg s^{1} and can have a remarkable impact on the lifetimes of closein planets, possibly explaining the observed increase of the mean planet mass with the inverse of the semimajor axis in systems with closein planets.
In the present investigation, we have considered stars in the middle of their mainsequence evolution. Both surface magnetic fields and EUV fluxes are larger in young stars, but the ratio of their contributions to the heating of a planetary atmosphere is only weakly dependent on stellar age. Specifically, SanzForcada et al. (2011) propose that F_{EUV} scales as τ^{1.24}, where τ is the stellar age, while the magnetic field strength can be assumed to scale as B_{0} ∝ τ^{0.5} in order to account for Skumanich’s law of stellar rotation evolution (cf., e.g., Kawaler 1988). Since , if B_{p}/B_{0} is assumed to be constant, the ratio P_{mag}/P_{EUV} ∝ τ^{0.24}, only weakly dependent on stellar age. However, in premainsequence stars the EUV flux tend to saturate, while surface fields B_{0} ~ 10^{3} G are observed (e.g., Yang & JohnsKrull 2011). In this phase, magnetic heating dominates and leads to a remarkable evaporation of closein planets, although limited to a few tens of Myr. Therefore, a detailed investigation is needed to predict the impact on the present mass distribution of closein planets.
Another characteristic of stellar magnetic fields is their variability on widely different timescales that can account for the remarkable variation in mass loss rate observed in, e.g., HD 189733b (Lecavelier des Etangs et al. 2012). On the other hand, on the basis of the observed Xray flux variability, the relative amplitude of the EUV flux variability is predicted to be smaller in stars more active than the Sun (cf., e.g., Sect. 6.1.3 of Claire et al. 2012), thus making it difficult to account for a large variability in the mass loss rate if evaporation is powered solely by the EUV flux.
Acknowledgments
The author wishes to thank an anonymous referee for several valuable comments that allowed him to improve this work. He is grateful to Drs. A. S. Bonomo, A. Collier Cameron, S. Desidera, R. Gratton, I. Pagano, E. Shkolnik, and A. Vidotto for interesting discussions on starplanet interaction and planetary evaporation in systems hosting hot Jupiters. The use of solar spectra obtained from the Virtual Planetary Laboratory web site (http://depts.washington.edu/naivpl/content/models/solarflux/) is also gratefully acknowledged.
References
 Adams, F. C. 2011, ApJ, 730, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Adams, F. C., Hollenbach, D., Laughlin, G., & Gorti, U. 2004, ApJ, 611, 360 [NASA ADS] [CrossRef] [Google Scholar]
 Adams, F. C., Cai, M. J., Galli, D., Lizano, S., & Shu, F. H. 2011, ApJ, 743, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J. 2002, Space Sci. Rev., 101, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Batygin, K., & Stevenson, D. J. 2010, ApJ, 714, L238 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, J. C. 1973, Sol. Phys., 28, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Christensen, U. R., Holzwarth, V., & Reiners, A. 2009, Nature, 457, 167 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Claire, M. W., Sheets, J., Cohen, M., et al. 2012, ApJ, 757, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Cohen, O., Kashyap, V. L., Drake, J. J., Sokolov, I. V., & Gombosi, T. I. 2011a, ApJ, 738, 166 [NASA ADS] [CrossRef] [Google Scholar]
 Cohen, O., Kashyap, V. L., Drake, J. J., et al. 2011b, ApJ, 733, 67 [NASA ADS] [CrossRef] [Google Scholar]
 Cox, D. P., & Tucker, W. H. 1969, ApJ, 157, 1157 [NASA ADS] [CrossRef] [Google Scholar]
 Dalgarno, A., & McCray, R. A. 1972, ARA&A, 10, 375 [NASA ADS] [CrossRef] [Google Scholar]
 Dere, K. P. 2007, A&A, 466, 771 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dere, K. P., Landi, E., Young, P. R., et al. 2009, A&A, 498, 915 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Donati, J.F., & Landstreet, J. D. 2009, ARA&A, 47, 333 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Donati, J.F., Moutou, C., Farès, R., et al. 2008, MNRAS, 385, 1179 [NASA ADS] [CrossRef] [Google Scholar]
 Emslie, A. G. 1978, ApJ, 224, 241 [NASA ADS] [CrossRef] [Google Scholar]
 Fares, R., Donati, J.F., Moutou, C., et al. 2009, MNRAS, 398, 1383 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Fares, R., Donati, J.F., Moutou, C., et al. 2012, MNRAS, 423, 1006 [NASA ADS] [CrossRef] [Google Scholar]
 Flyer, N., Fornberg, B., Thomas, S., & Low, B. C. 2004, ApJ, 606, 1210 [NASA ADS] [CrossRef] [Google Scholar]
 Gu, P.G., & Suzuki, T. K. 2009, ApJ, 705, 1189 [NASA ADS] [CrossRef] [Google Scholar]
 Gurdemir, L., Redfield, S., & Cuntz, M. 2012, PASA, 29, 141 [Google Scholar]
 Haug, E. 2004, A&A, 423, 793 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Herrero, E., Lanza, A. F., Ribas, I., Jordi, C., Morales, J. C. 2013, A&A, 553, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Holman, G. D. 2005, Adv. Space Res., 35, 1669 [NASA ADS] [CrossRef] [Google Scholar]
 Holman, G. D., Aschwanden, M. J., Aurass, H., et al. 2011, Space Sci. Rev., 159, 107 [NASA ADS] [CrossRef] [Google Scholar]
 Janse, Å. M., & Low, B. C. 2009, ApJ, 690, 1089 [NASA ADS] [CrossRef] [Google Scholar]
 Janse, Å. M., & Low, B. C. 2010, ApJ, 722, 1844 [NASA ADS] [CrossRef] [Google Scholar]
 Kashyap, V. L., Drake, J. J., & Saar, S. H. 2008, ApJ, 687, 1339 [NASA ADS] [CrossRef] [Google Scholar]
 Kawaler, S. D. 1988, ApJ, 333, 236 [NASA ADS] [CrossRef] [Google Scholar]
 Lanza, A. F. 2008, A&A, 487, 1163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lanza, A. F. 2009, A&A, 505, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lanza, A. F. 2011, Ap&SS, 336, 303 [NASA ADS] [CrossRef] [Google Scholar]
 Lanza, A. F. 2012, A&A, 544, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leach, J., & Petrosian, V. 1981, ApJ, 251, 781 [NASA ADS] [CrossRef] [Google Scholar]
 Lecavelier Des Etangs, A. 2007, A&A, 461, 1185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lecavelier Des Etangs, A., VidalMadjar, A., & Desert, J.M. 2008, Nature, 456, 1 [Google Scholar]
 Lecavelier Des Etangs, A., Ehrenreich, D., VidalMadjar, A., et al. 2010, A&A, 514, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lecavelier des Etangs, A., Bourrier, V., Wheatley, P. J., et al. 2012, A&A, 543, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lecavelier des Etangs, A., Sirothia, S. K., GopalKrishna, & Zarka, P. 2013, A&A, 552, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Linsky, J. L., Yang, H., France, K., et al. 2010, ApJ, 717, 1291 [NASA ADS] [CrossRef] [Google Scholar]
 Mann, G., & Warmuth, A. 2011, A&A, 528, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 McIvor, T., Jardine, M., & Holzwarth, V. 2006, MNRAS, 367, L1 [NASA ADS] [CrossRef] [Google Scholar]
 McTiernan, J. M., & Petrosian, V. 1990, ApJ, 359, 524 [NASA ADS] [CrossRef] [Google Scholar]
 Miller, J. A., Larosa, T. N., & Moore, R. L. 1996, ApJ, 461, 445 [NASA ADS] [CrossRef] [Google Scholar]
 Moutou, C., Donati, J.F., Savalle, R., et al. 2007, A&A, 473, 651 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 MurrayClay, R. A., Chiang, E. I., & Murray, N. 2009, ApJ, 693, 23 [Google Scholar]
 Petit, P., Dintrans, B., Solanki, S. K., et al. 2008, MNRAS, 388, 80 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Petrosian, V. 2012, Space Sci. Rev., 173, 535 [NASA ADS] [CrossRef] [Google Scholar]
 Poppenhaeger, K., & Schmitt, J. H. M. M. 2011, ApJ, 735, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Poppenhaeger, K., Robrade, J., & Schmitt, J. H. M. M. 2010, A&A, 515, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Poppenhaeger, K., Lenz, L. F., Reiners, A., Schmitt, J. H. M. M., & Shkolnik, E. 2011, A&A, 528, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Preusse, S., Kopp, A., Büchner, J., & Motschmann, U. 2006, A&A, 460, 317 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Priest, E. R. 1984, Solar Magnetohydrodynamics, Geophysics and Astrophysics Monographs (Dordrecht: Reidel) [Google Scholar]
 Reiners, A., & Christensen, U. R. 2010, A&A, 522, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rudge, M. R. H. 1968, Rev. Mod. Phys., 40, 564 [NASA ADS] [CrossRef] [Google Scholar]
 Rüdiger, G., & Hollerbach, R. 2004, The Magnetic Universe: Geophysical and Astrophysical Dynamo Theory (WileyVCH), Chap. 2 [Google Scholar]
 Saar, S. H., Cuntz, M., Kashyap, V. L., & Hall, J. C. 2008, IAU Symp., 249, 79 [NASA ADS] [Google Scholar]
 SanzForcada, J., Micela, G., Ribas, I., et al. 2011, A&A, 532, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saur, J., Grambusch, T., Duling, S., Neubauer, F. M., & Simon, S. 2013, A&A, 552, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Scandariato, G., Maggio, A., Lanza, A. F., et al. 2013, A&A, 552, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schlickeiser, R., Pohl, M., & Vainio, R. 2003, ApJ, 596, 840 [NASA ADS] [CrossRef] [Google Scholar]
 Shkolnik, E. L. 2013, ApJ, 766, 9 [NASA ADS] [CrossRef] [Google Scholar]
 Shkolnik, E., Walker, G. A. H., Bohlender, D. A., Gu, P.G., & Kürster, M. 2005, ApJ, 622, 1075 [NASA ADS] [CrossRef] [Google Scholar]
 Shkolnik, E., Bohlender, D. A., Walker, G. A. H., & Collier Cameron, A. 2008, ApJ, 676, 628 [NASA ADS] [CrossRef] [Google Scholar]
 Sing, D. K., VidalMadjar, A., Lecavelier des Etangs, A., et al. 2008, ApJ, 686, 667 [Google Scholar]
 Störzer, H., & Hollenbach, D. 1999, ApJ, 515, 669 [Google Scholar]
 Trammell, G. B., Arras, P., & Li, Z.Y. 2011, ApJ, 728, 152 [NASA ADS] [CrossRef] [Google Scholar]
 Venot, O., Hébrard, E., Agúndez, M., et al. 2012, A&A, 546, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 VidalMadjar, A., Lecavelier des Etangs, A., Désert, J.M., et al. 2003, Nature, 422, 143 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 VidalMadjar, A., Sing, D. K., Lecavelier Des Etangs, A., et al. 2011, A&A, 527, A110 [Google Scholar]
 Vidotto, A. A., Jardine, M., & Helling, C. 2010, ApJ, 722, L168 [NASA ADS] [CrossRef] [Google Scholar]
 Watson, A. J., Donahue, T. M., & Walker, J. C. G. 1981, Icarus, 48, 150 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, H., & JohnsKrull, C. M. 2011, ApJ, 729, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, M., & Flyer, N. 2008, ApJ, 683, 1160 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, M., Flyer, N., & Low, B. C. 2006, ApJ, 644, 575 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Treatment of relativistic electrons
When the cutoff energy E_{c} ≳ 0.2 m_{e}c^{2}, the energy loss of the highenergy electrons during their motion through the plasma must be treated in a relativistic way. For energies up to 10 MeV, losses by bremsstrahlung, synchrotron radiation, and Compton scattering are at least two orders of magnitude smaller than the collisional loss considered here (cf., e.g., Fig. 1 in Haug 2004; McTiernan & Petrosian 1990). In the case of a fully ionized hydrogen plasma, the relativistic energy loss rate per electron is (e.g. Eq. (10) in Haug 2004): (A.1)where E is the kinetic energy of the particle, r_{0} = e^{2}/m_{e}c^{2} the classic electron radius, with m_{e} the electron mass and c the speed of light, Λ_{ee} the Coulomb logarithm for electronelectron collisions, n the total number density of protons, and β ≡ v/c with v the electron speed, not to be confused with the plasma β used before. The kinetic energy is expressed in units of the rest energy as: E ≡ ϵ m_{e}c^{2}, with: (A.2)Considering that K′ ≡ 2πe^{4}Λ_{ee} and dz = cβdt, Eq. (A.1) can be recast in the form: (A.3)that represents the relativistic generalization of Eq. (15) and reduces to it when ϵ ≪ 1. On the other hand, for highly relativistic particles ϵ ≫ 1 and dE/dz = −2(K′/m_{e}c^{2})n becomes independent of the energy. Note that our simple treatment of partial ionization by including a factor (λ + x) in the r.h.s. of Eq. (15) is not rigorously valid in the relativistic regime (cf., e.g., McTiernan & Petrosian 1990), but we shall extend its validity given that the introduced error in the energy loss rate does not exceed a factor of 4 − 5 that can be absorbed into the other uncertainties of our model.
To compute the energy loss per proton Q_{c}, we need to specify the electron flux that is in the relativistic regime. For electron energies ϵ ≤ 0.2 we shall continue to use Eq. (17), while for greater energies we shall neglect the modification of the spectrum produced by the loss of energy in the layers above the level z because it is negligible at high energies. Specifically, to stop an electron with an energy E ≃ 0.2 m_{e}c^{2} ~ 100 keV, we need a column density of ~2.7 × 10^{25} m^{2}, that is much larger than the column density of 10^{22} − 10^{23} m^{2} corresponding to the layer from which the evaporation flow is launched. Therefore, the adopted expression for Q_{c} in the relativistic regime taking into account Eq. (A.3) is: (A.4)where K′′, , α_{e} ≡ E_{s}(z)/E_{c}, and w = E_{0}/E_{c} have already been introduced in the case of Eq. (17), while α_{r} ≡ m_{e}c^{2}/E_{c}. In the relativistic regime, the total electron energy flux used to compute n_{T} is given by the expression: (A.5)where ϵ_{c} ≡ E_{c}/m_{e}c^{2} and (A.6)All the integrals in Eqs. (A.4) and (A.6) are evaluated numerically. Note that the expression for the normalization constant A of the electron spectrum derived from Eq. (12) is still valid in the relativistic regime, thus we can determine the cutoff energy E_{c} as in the classic regime.
The asymptotic expression for the collisional ionization cross section is not modified in the relativistic regime (cf., e.g., Eq. (12) in Schlickeiser et al. 2003). Therefore, we continue to use our formulae by adopting the relativistic expression for the electron flux, i.e., .
All Figures
Fig. 1
Upper panel: a meridional section of the stellar and planetary magnetic field lines in a magnetohydrostatic regime in which the stellar field is forcefree and the planetary field is potential. In this case the flux systems of the two bodies are topologically separated (see Sect. 3.1). The field lines of the star and the planet can release energy by reconnecting at the boundary between their magnetospheres that is highlighted with a red dashed contour. Lower panel: the magnetic field of the starplanet system in a magnetohydrostatic regime in which the stellar field is potential. In this case, field lines interconnecting the two bodies exist. Interconnecting magnetic loops can be stressed by moving their footpoints and accumulate energy that can be successively released by a relaxation to the minimumenergy potential field configuration (see Sect. 3.1 for details). 

In the text 
Fig. 2
Upper panel: the ratio P_{mag}/P_{EUV} vs. the orbital semimajor axis a for the known transiting planets computed for a stellar radial magnetic field (s = 2). The value of P_{mag} given by Eq. (4) has been multiplied by 1/6 (see the text). The dotted horizontal line corresponds to P_{mag} = P_{EUV}. Middle panel: the same as the upper panel, but for s = 2.5 corresponding to a nonlinear forcefree stellar field. Lower panel: the same as the upper panel, but for s = 3, corresponding to a potential dipole stellar field. 

In the text 
Fig. 3
Cutoff energy of nonthermal electron spectrum vs. orbital semimajor axis for our sample of transiting planets for a spectrum with an index δ = 6 and a stellar field decaying with distance with s = 3 (cf. Eq. (3)). Solid dots refer to a stellar field B_{0} ≥ 10 G, while open diamonds indicate B_{0} < 10 G (cf. Sect. 2.2). 

In the text 
Fig. 4
Upper panel: ratio of the mass loss rate induced by nonthermal electrons and EUV photons to that induced by EUV photons only vs. the semimajor axis of the orbit for our sample of transiting planets (see the text). The parameter s = 2.0 specifies the radial dependence of the stellar magnetic field according to Eq. (3). Middle panel: the same as the upper panel, but for s = 2.5. Lower panel: the same as the upper panel, but for s = 3.0. 

In the text 
Fig. 5
Upper panel: mass loss rate for a sample of transiting planets vs. their mass (see the text). Open diamonds indicate the loss rate induced by EUV flux only, while filled dots indicate the case when the contribution by nonthermal electrons accelerated at the boundary of the planetary magnetosphere is added. The parameter s = 2.0 specifies the radial dependence of the stellar magnetic field according to Eq. (3). Middle panel: the same as the upper panel, but for s = 2.5. Lower panel: the same as the upper panel, but for s = 3.0. 

In the text 
Fig. 6
Lifetime for the considered sample of transiting planets vs. their mass. The lifetime is estimated considering the maximum mass loss rate, i.e., by including the contribution of nonthermal electrons accelerated by reconnection at the boundary of the planetary magnetosphere. 

In the text 
Fig. 7
Relative variation of the radius corresponding to optical depth τ = 1 in the centre of the Lymanα line for normal incidence across the planetary atmosphere vs. the planet mass. The plotted case corresponds to the slowest decay of the stellar field with distance (s = 2.0) that provides an upper limit to the predicted variation. 

In the text 
Fig. 8
Upper panel: mass loss rate averaged over time (in kg s^{1}) vs. planet mass in the case of the dissipation of magnetic energy in a loop interconnecting the star with the planet. Lower panel: the estimated lifetime of evaporating planets vs. their mass. The horizontal dashed line marks an evaporation lifetime of 5 Gyr corresponding to the mean estimated age of solarlike mainsequence stars. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.