 as the energy-limited approximation [Watson1981]{} and allows a relatively simple analytic description of mass loss rates.

$$\label{masslosseq}
\dot{M}_{\mathrm{e-lim}} \approx \frac{\epsilon \pi F_{\mathrm{XUV}} R_{\mathrm{XUV}}^3}{G M_{\mathrm{p}} K_{\mathrm{tide}}}$$

$$K_{\mathrm{tide}} = (1 - \frac{3}{2 \xi} + \frac{1}{2 \xi^3})$$

$$\xi = \frac{R_{\mathrm{Hill}}}{R_{\mathrm{XUV}}}$$

Equation (\[masslosseq\]) describes our estimate of the mass loss rate based on the formulation from [Erkaev2007]{}. $F_{\mathrm{XUV}}$ is the total flux between $1-1200$ Å, which is given by Ribas(2005) for Sun-like stars. For stars older than 100 Myr, Ribas found that at 1 AU $F_{\mathrm{XUV}} = 29.7 \tau^{-1.23} \: \mathrm{erg \, s^{-1} \, cm^{-2}}$, where $\tau$ is the age of the star in Gyr. Using this power law, we scale the XUV flux to the appropriate age and semi-major axis for each planet in our models. Although Ribas only targeted Sun-like stars, [Sanz-Forcada2010]{} found similar results for a wide range of stellar types from M3 to F7. Hereafter, we will simply refer to the entire $1-1200$ Å$;$ spectrum as XUV.

$R_{\mathrm{XUV}}$ is the planetary radius at which the atmosphere becomes optically thick to XUV photons, which [Murray-Clay2009]{} find occurs at pressures around a nanobar, in the hot Jupiter context. For our work, we assume that the atmosphere is isothermal between the optical and XUV photospheres. This neglects heating from photo-disassociation, which should occur around a $\mathrm{\mu bar}$ [Kempton2012a]{}. However, this effect should be relatively small and if anything will lead to slight underestimate of the mass loss rate. We vary pressure of the XUV photosphere from 0.1 nbar to 10 nbar to include the uncertainty in the structure of the XUV photosphere. For H/He atmosp
as the energy - limited approximation [ Watson1981 ] { } and allows a relatively childlike analytic description of aggregate loss rates. 

 $ $ \label{masslosseq } 
 \dot{M}_{\mathrm{e - lim } } \approx \frac{\epsilon \pi F_{\mathrm{XUV } } R_{\mathrm{XUV}}^3}{G M_{\mathrm{p } } K_{\mathrm{tide}}}$$ 

 $ $ K_{\mathrm{tide } } = (1 - \frac{3}{2 \xi } + \frac{1}{2 \xi^3})$$ 

 $ $ \xi = \frac{R_{\mathrm{Hill}}}{R_{\mathrm{XUV}}}$$ 

 Equation (\[masslosseq\ ]) trace our appraisal of the mass loss rate establish on the formulation from [ Erkaev2007 ] { }. $ F_{\mathrm{XUV}}$ is the total magnetic field between $ 1 - 1200 $ Å, which is given by Ribas(2005) for Sun - like star. For stars older than 100 Myr, Ribas find oneself that at 1 AU $ F_{\mathrm{XUV } } = 29.7 \tau^{-1.23 } \: \mathrm{erg \, s^{-1 } \, cm^{-2}}$, where $ \tau$ is the age of the star in Gyr. Using this power law, we scale the XUV flux to the appropriate historic period and semi - major axis for each planet in our models. Although Ribas entirely targeted Sun - like stars, [ Sanz - Forcada2010 ] { } found alike results for a wide range of stellar types from M3 to F7. hereinafter, we will simply refer to the entire $ 1 - 1200 $ Å$;$ spectrum as XUV. 

 $ R_{\mathrm{XUV}}$ is the planetary radius at which the atmosphere become optically thick to XUV photons, which [ Murray - Clay2009 ] { } find occurs at pressures around a nanobar, in the hot Jupiter context. For our work, we assume that the atmosphere is isothermal between the optical and XUV photosphere. This neglect heating from photograph - dissociation, which should occur around a $ \mathrm{\mu bar}$ [ Kempton2012a ] { }. However, this effect should be relatively small and if anything will lead to slender underestimate of the mass loss rate. We vary atmospheric pressure of the XUV photosphere from 0.1 nbar to 10 nbar to include the uncertainty in the structure of the XUV photosphere. For H / He atmosp
 as the energy-limited approwimation [Watson1981]{} cbd allmws a delativeuy simple analytic descriptiln of mqss loss rates.

$$\label{masrlosseq}
\don{M}_{\mathrm{e-oim}} \epprox \frac{\epsilon \pi F_{\mathrm{XHY}} R_{\machcm{XUV}}^3}{G M_{\mathrm{p}} K_{\mathrm{tige}}}$$

$$K_{\mathrm{tide}} = (1 - \fxac{3}{2 \xi} + \frac{1}{2 \xi^3})$$

$$\xi = \frac{R_{\mathrm{Hill}}}{R_{\mwthrm{XUF}}}$$

Eeuation (\[masslofseq\]) qescdpbts our estimate of the mass loss date baved on the fotmulation from [Erkaev2007]{}. $F_{\matjrm{XKV}}$ is the total flkx between $1-1200$ Å, whyxh is given cy Ribas(2005) for Sun-like sjars. For stars older than 100 Myr, Rkbas yound that qt 1 AK $F_{\mathrm{XUV}} = 29.7 \tal^{-1.23} \: \mathrm{erg \, s^{-1} \, cm^{-2}}$, wvere $\tai$ is the age on the stqr in Gyr. Using this 'ower law, we scale tre XUV flgx to the appropriare age dnd vemi-oqjof asix ror eafh 'lanet in ohr models. Aothough Ribas only uardvyed Sun-like atars, [Fagz-Forcada2010]{} found similar results for a wpde dange of stellar types drom M3 to F7. Hereafter, we will fimply refer to the entire $1-1200$ Å$;$ spectrum as XUV.

$R_{\mathsm{XUV}}$ ks uhc ooajetary radius at which the atmosphere becomes kpuicslly thick to WUV photons, which [Kugrsi-Clay2009]{} find occors at 'deasures around a najobar, ig the hot Jupiuer cpntext. For our work, we assune that the ctmisphere is isotherlal between the optival and XUV photospherer. Thjs neglects heating rfom photo-disassozianion, which should occur around a $\mathrm{\nu bcr}$ [Kemptun2012a]{}. Nowevew, this efffct smmuld be relatively smalp dnd if anyhhing will lead to slight underxxtimate of tne mavs loss xate. Wc vary pressure of the XUV phptosphexe froo 0.1 nbar to 10 nbar vo include tre uncertaintf in the strurture of ehe ZUV photospfdre. For H/He atkosp
as the energy-limited approximation [Watson1981]{} and allows simple description of loss rates. $$\label{masslosseq} R_{\mathrm{XUV}}^3}{G K_{\mathrm{tide}}}$$ $$K_{\mathrm{tide}} = - \frac{3}{2 \xi} \frac{1}{2 \xi^3})$$ $$\xi = \frac{R_{\mathrm{Hill}}}{R_{\mathrm{XUV}}}$$ Equation describes our estimate of the mass loss rate based on the formulation from $F_{\mathrm{XUV}}$ is the total flux between $1-1200$ Å, which is given by Ribas(2005) Sun-like For older 100 Myr, Ribas found that at 1 AU $F_{\mathrm{XUV}} = 29.7 \tau^{-1.23} \: \mathrm{erg \, s^{-1} cm^{-2}}$, where $\tau$ is the age of the in Gyr. Using this law, we scale the XUV to appropriate age semi-major for planet in our Although Ribas only targeted Sun-like stars, [Sanz-Forcada2010]{} found similar results for a wide range of stellar types M3 to we will refer the $1-1200$ Å$;$ spectrum $R_{\mathrm{XUV}}$ is the planetary radius at becomes optically thick to XUV photons, which [Murray-Clay2009]{} occurs at around a nanobar, in the hot context. For our work, we assume that the is isothermal between the optical and XUV photospheres. This neglects heating from photo-disassociation, which should a $\mathrm{\mu bar}$ [Kempton2012a]{}. this effect should relatively and anything lead to underestimate of the mass loss rate. We vary pressure of the photosphere from 0.1 nbar to 10 nbar to include the the of the XUV For H/He atmosp
 as the energy-limited approxiMation [WatsOn1981]{} and AllOws A rElatIvelY simple analytiC DescRiption of mass loss rates.

$$\Label{MaSSlosSEq}
\Dot{M}_{\mAthrm{e-lIM}} \aPPRox \FrAc{\EpsIlON \pI F_{\matHrm{xUV}} R_{\matHrm{XUV}}^3}{G M_{\maThrM{p}} k_{\mathrm{tide}}}$$

$$K_{\MAtHrm{tide}} = (1 - \fraC{3}{2 \xi} + \Frac{1}{2 \xi^3})$$

$$\xi = \frac{r_{\maThrm{HiLl}}}{r_{\maTHrm{XUv}}}$$

EqUatioN (\[masslOSseq\]) deScribes ouR eSTimate OF the masS LOsS ratE based on the formulATiON from [Erkaev2007]{}. $F_{\maThrm{XUv}}$ iS ThE TOtaL flUx between $1-1200$ Å, WhIch is GIven by RIBaS(2005) FOR SuN-Like stars. For sTars older thAN 100 MyR, Ribas FoUnd THat at 1 Au $F_{\matHrM{xUV}} = 29.7 \Tau^{-1.23} \: \mathrm{erG \, s^{-1} \, cm^{-2}}$, Where $\tau$ iS the agE Of the stAR in Gyr. USing thIs pOweR law, WE sCaLe tHe xuV fLUx To tHE apPropriatE aGe And seMi-maJOR AXis fOr eAch pLanet In our models. AlThoUgh RIBas Only tArgetEd SuN-lIke stArs, [SanZ-ForcAdA2010]{} found similar reSultS for a wide RanGe Of sTeLlar tYPes froM M3 tO F7. HEreafteR, we will SImpLy REFEr To the entire $1-1200$ Å$;$ spectrUm AS xUv.

$R_{\mathrm{xUV}}$ is tHE pLaNEtary radIuS at WhicH THe atmOsphERe Becomes oPticalLY tHiCk to XUV PhOtons, wHiCh [MUrrAy-ClaY2009]{} Find Occurs At pressuRes arOUnd a nanobar, in tHE hot Jupiter coNTeXT. foR Our wOrk, We assume thaT the ATmosPherE Is IsoTHermaL betwEeN ThE Optical and XUV photosPhEres. ThIs negLects heating fRom photo-diSASSociatioN, whiCH sHOuld occur arounD a $\matHrm{\mu bar}$ [KeMPton2012a]{}. HowEver, tHis effecT should be RELatively SmaLl aNd iF anYTHiNg will lead to sLIGht uNdErestimAte Of the maSs lOss RatE. We VaRy pressurE of the XUv pHoToSpHerE from 0.1 NBar to 10 nbaR tO inClUde The unCErtainTy in tHe stRuCtURe oF the XUV PHoTOSpheRe. foR H/He AtmOsP
 as the energy-limited app roximation  [Wat son 198 1] {} a nd a llows a relati v elysimple analytic descri ption  o f  mas s  l oss r ates.
$ $\ l a bel {m as slo ss e q} 
\dot {M} _{\math rm{e-lim}}  \a pp rox \frac{\e p si lon \pi F_ {\m athrm{XUV}}R_{ \mathr m{ XUV } }^3}{ G M _{\ma thrm{p } } K_{\ mathrm{ti de } }}$$
$ $K_{\ma t h rm {tid e}} = (1 - \frac{ 3 }{ 2  \xi} + \frac{ 1}{2 \ xi ^ 3} ) $ $
$$\ xi = \frac {R _{\ma t hrm{Hil l }} } { R _{\ m athrm{XUV}}}$ $

Equation (\[ masslo ss eq\ ] ) desc ribes  o u r e stimate ofthemass loss  rateb ased on the for mulati onfro m [E r ka ev 200 7] { }.$ F_ {\m a thr m{XUV}}$  i sthe t otal f l u x be twe en $ 1-120 0$ Å, which i s g iven byRibas (2005 ) fo rSun-l ike st ars.Fo r stars older t han100 Myr,Rib as  fo un d tha t  at 1AU$F_ {\mathr m{XUV}} = 2 9. 7 \ ta u^{-1.23} \: \math rm { e rg  \, s^{- 1} \,c m^ {- 2 }}$, whe re  $\ tau$ i s the  age of  the sta r in G y r.  U sing th is  power  l aw,  we  scal e  the  XUV f lux to t he ap p ropriate age a n d semi-majora xi s fo r  eac h p lanet in ou r mo d els.  Alt h ou ghR ibasonlyta r ge t ed Sun-like stars,[S anz-Fo rcada 2010]{} found  similar r e s u lts fora wi d er ange of stella r typ es from M3 to F7. H ereaf ter, wewill simp l y  refer t o t heent ire $ 1- 1200$ Å$;$ sp e c trum  a s XUV.
$R _{\math rm{ XUV }}$  is  t he planet ary radi us  a twh ich  thea tmospher ebec om esoptic a lly th ick t o XU Vph o ton s, whic h  [ M u rray -C la y200 9]{ }findoccu r s a t press ures arou nda  nan ob ar , in th e hot Jupiter  c ontext. Fo rour  work, w e assume  that the atmosphere is isother mal  betw eenthe optic aland XU V p h otosph eres.Thisne gle c t s hea t i ng  fr om  photo-dis a s soc iatio n,  whi ch shou ld occur around a$ \ma thrm{\mu bar} $ [ Kemp t o n2 012 a ]{ } . H ow e ver , this effect sho uld be rel at i ve ly small a n d i fanythin g willleadt o sligh t underes timate of  t he m a s s l oss rate.We varypressureo f the XU V pho tos pherefr om0.1 n bar to 10nbarto inc lu de the  unce rt ainty in  the structure of the X UV pho tosph ere . For H/H e a t mos p
 as_the energy-limited_approximation [Watson1981]{} and allows_a relatively_simple_analytic description_of_mass loss rates.

$$\label{masslosseq}
\dot{M}_{\mathrm{e-lim}}_\approx \frac{\epsilon \pi_F_{\mathrm{XUV}} R_{\mathrm{XUV}}^3}{G M_{\mathrm{p}} K_{\mathrm{tide}}}$$

$$K_{\mathrm{tide}}_= (1 -_\frac{3}{2_\xi} + \frac{1}{2 \xi^3})$$

$$\xi = \frac{R_{\mathrm{Hill}}}{R_{\mathrm{XUV}}}$$

Equation (\[masslosseq\]) describes our estimate of the mass loss_rate_based on_the_formulation_from [Erkaev2007]{}. $F_{\mathrm{XUV}}$ is the_total flux between $1-1200$ Å,_which is_given by Ribas(2005) for Sun-like stars. For stars_older_than 100 Myr,_Ribas found that at 1 AU $F_{\mathrm{XUV}} = 29.7_\tau^{-1.23} \: \mathrm{erg \, s^{-1} \,_cm^{-2}}$, where $\tau$_is_the_age of the star_in Gyr. Using this power law,_we scale the XUV flux to_the appropriate age and semi-major axis for_each planet in our models. Although_Ribas only targeted Sun-like stars,_[Sanz-Forcada2010]{} found_similar results for a wide_range of stellar_types from_M3 to F7._Hereafter, we will simply refer to_the entire $1-1200$_Å$;$ spectrum as XUV.

$R_{\mathrm{XUV}}$ is the_planetary_radius at which_the_atmosphere_becomes optically_thick to XUV_photons,_which [Murray-Clay2009]{}_find_occurs at pressures around a nanobar,_in_the hot Jupiter context. For our work,_we assume that the_atmosphere_is isothermal between the_optical and XUV photospheres. This_neglects heating from photo-disassociation, which should_occur around_a $\mathrm{\mu_bar}$ [Kempton2012a]{}. However, this effect should be relatively small and if_anything will lead to slight underestimate_of the mass loss_rate. We_vary_pressure of the_XUV_photosphere from_0.1 nbar to 10 nbar to include_the uncertainty_in the structure of the XUV_photosphere. For H/He atmosp