arxiv: v2 [astro-ph] 1 Oct 2007

Transcription

1 Mon. Not. R. Astron. So. 000, 1 25 (2007) Printed 20 June 2018 (MN LATEX style file v2.2) Long-term tidal evolution of short-eriod lanets with omanions arxiv: v2 [astro-h] 1 Ot 2007 Rosemary A. Mardling 1 1 Shool of Mathematial Sienes, Monash University, Vitoria, 3800, Australia Aeted... Reeived...; in original form... ABSTRACT Of the fourteen transiting extrasolar lanetary systems for whih radii have been measured, at least three aear to be onsiderably larger than theoretial estimates suggest. It has been roosed by Bodenheimer, Lin & Mardling that undeteted omanions ating to exite the orbital eentriity are resonsible for these oversized lanets, as they find new equilibrium radii in resonse to being tidally heated. In the ase of HD , this hyothesis has been rejeted by some authors beause there is no sign of suh a omanion at the 5 ms 1 level, and beause it is diffiult to say onlusively that the eentriity is non-zero. Transit timing analysis as well as a diret transit searh have further onstrained the existene of very short-eriod omanions, eseially in resonant orbits. Whether or not a omanion is resonsible for the large radius of HD b, almost ertainly some short-eriod systems have omanions whih fore their eentriities to nonzero values. This aer is dediated to quantifying this effet. The eentriity of a short-eriod lanet will only be exited as long as its (nonresonant) omanion s eentriity is non-zero. Here we show that the latter deays on a timesale whih deends on the struture of the interior lanet, a timesale whih is often shorter than the lifetime of the system. This inludes Earth-mass lanets in the habitable zones of some stars. We determine whih onfigurations are aable of sustaining signifiant eentriity for at least the age of the system, and show that these inlude systems with omanion masses as low as a fration of an Earth mass. The orbital arameters of suh omanions are onsistent with reent alulations whih show that the migration roess an indue the formation of low mass lanets external to the orbits of hot Juiters. Systems with inflated lanets are therefore good targets in the searh for terrestrial lanets. Key words: lanetary systems elestial mehanis stellar dynamis methods: analytial lanetary systems: formation 1 INTRODUCTION While most extrasolar lanet disoveries to date have relied on the radial veloity method of detetion with its inherent unertainty in the lanetary mass, a few systems have been observed in transit, thereby breaking this degeneray. Not only does this allow a reasonably aurate determination of the lanetary mass (given a reliable estimate for the stellar mass is available), it allows one to estimate the radius of the lanet, and hene its density. Before the disovery that HD b transited its host star (Charbonneau et al. 2000; Henry et al. 2000), it was by no means lear that the lanets disovered till then, all with masses similar to or more than Juiter s, were gas giants. The extra information gleaned from hotometri transit data has rovided the oortunity to study lanetary struture (Bodenheimer, Lin & Mardling 2001), while setrosoi data have revealed their hemial omosition (Charbonneau et al. 2002). In the former ase, determination of the lanetary radius mardling@si.monash.edu.au

2 2 Rosemary A. Mardling Table 1. Observational data for inflated lanets Planet Period (day) m (M ) m (M J ) a (AU) e obs R (R J ) HD b ± ±0.009 a 1.32±0.03 WASP-1b b ± HAT-P-1b ± ± a Laughlin et al. (2005a). b Stemels et al. (2007). This value is based on very little data (Bakos et al. 2006). allows one to ut onstraints on the ossible existene of a lanetary ore, and this in turn allows one to onstrain the mode of lanet formation, at least for that lanet. In rinile this allows a hoie between the ore aretion model (Safronov 1969; Wetherill 1980; Bodenheimer & Pollak 1986) and the gravitational instability model (Kuier 1951; Cameron 1978; Boss 1997), although it is not lear that the latter should rodue ore-less lanets given aretion of solids lus subsequent gravitational settling almost ertainly follows formation. Either way, suh an exerise requires knowledge of the equilibrium struture of a gas giant lanet, with and without a ore, and with and without external and internal heating. While it an be argued that there are roblems mathing theoretial radii with observation for all transiting lanets (Burrows et al. 2007), three stand out as artiularly exessive. These are HD b (Charbonneau et al. 2000; Henry et al. 2000; Knutson et al. 2007), WASP-1b (Collier Cameron et al. 2006) and HAT-P-1b (Bakos et al. 2006), for whih the lanetary radii and other observationally derived arameters are listed in Table 1. Data for HD and HAT-P-1 are taken from Burrows et al. (2006) and referenes therein, exet for the eentriities, while data for WASP-1 are taken from Stemels et al. (2007). Here m and m are the masses of the star and lanet resetively, a and e are the orbital semimajor axis and eentriity and R is the estimated lanetary radius (R J is Juiter s radius). The eentriity of HD b is from Laughlin et al. (2005a), a aer devoted seifially to disussing this ontentious orbital element, while that of HAT-P- 1b is from Bakos et al. (2006), a value based on very little data. However, as ointed out by the latter authors, HAT-P-1b is further from its arent star than HD b, and is suh that the relatively large estimated value of 0.09 rodues a similar lanetary radius. Various mehanisms have been suggested to aount for the exess radii of the lanets listed in Table 1. Obvious ontenders are stellar irradiation(burrows et al. 2000; Guillot & Showman 2002; Bodenheimer, Laughlin & Lin 2003; Burrows, Sudarsky & Hubbard 2003; Arras & Bildsten 2006; Fortney, Marley & Barnes 2007; Burrows et al. 2007) and tidal foring by a omanion lanet (Bodenheimer, Lin & Mardling 2001; Bodenheimer, Laughlin & Lin 2003), the latter requiring non-zero eentriity of the observed lanetary orbit to be viable. A less obvious mehanism has been roosed by Winn & Holman (2005) in whih obliquity tides are fored if the system an manage to remain in a suitable Cassini state. Both Levrard et al. (2007) and Fabryky, Johnson & Goodman (2007) verify that in rinile it is ossible for this mehanism to oerate as long as the system dissiates energy at a rate whih is too small to inflate the lanet to the observed size. This aer fouses on the mehanism of tidal foring rovided by a omanion lanet, and hene requires that the orbit of the transiting lanet be non-irular. Using a bootstra Monte Carlo tehnique to estimate onfidene intervals for the eentriity of HD b, Laughlin et al. (2005a) find that e < with 99.73% onfidene, e < with 95.4% onfidene, and e < 0.02 with 68.3% onfidene, with longitude of eriastron = 80 o ± 80 o. Their best-fit measurement for the eentriity is e = However, they oint out that beause e has a hard lower bound of zero, any best-fit model to a radial veloity data set will have e 0, with largest estimates orresonding to 90 o and 270 o. Using syntheti data orresonding to a truly irular orbit, suh bootstra trials rodue a distribution in e sae whih is symmetri about these two values, while use of the real radial veloity data rodues a distintly asymmetri distribution whih is biased towards = 90 o. This suggests that in fat e 0, although e = 0 annot be ruled out. Laughlin et al. (2005a) oint out that the fat that the estimate for is entred near 90 o suggests that even if HD b has a non-zero eentriity, timing of the seondary elise will find e os 0, as has reently been found by Deming et al. (2005) in the infrared. Note also that the orbital model of Laughlin et al. (2005a) leads to a stellar jitter estimate of 3.3 ms 1 whih is lose to the redited value of 2.8 ms 1. Thus, at least for HD , we onentrate our attention in this aer on omanions whih rodue stellar reflex veloities not muh larger than this value. The resene of a very short-eriod omanion to HD b has been ruled out by transit-timing analysis for interior and exterior lanets with eriods u to 21 days (orresonding to the 6:1 resonane at 0.16 AU) (Agol & Steffen 2007), as well as by a diret transit searh for lanets with radii down to 2R and eriods u to 14 days (Croll et al. 2007). Here we will onentrate on omanion lanets with semimajor axes greater than these lower bounds. This aer resents a generi analysis for short-eriod systems with external omanion lanets (low-mass or otherwise),

3 Comanions to short-eriod lanets 3 estimating the range of orbital arameters and masses required for suh omanions to be resonsible for inflated radii. An unexeted result is that very low mass lanets are sometimes aable of induing signifiant tidal foring on extremely long timesales. Setion 2 disusses the relationshi between the eentriity and the Q-value of a lanet for a given value of the lanet s luminosity, the Q-value or tidal quality fator (Goldreih & Soter 1966) being a measure of how effiiently the tides dissiate their osillation energy. Setion 3 disusses three-body dynamis in the resene of dissiation (see also Murray & Dermott (2000) and Adams & Laughlin (2006b)). We start by reviewing the seular evolution of oint-mass olanar three-body systems, with and without relativisti effets, deriving exliit formulae for the eriods of variation of the eentriities and their amlitudes. We identify the orbital arameters whih orresond to libration and irulation of the angle η, where is the longitude of eriastron of the omanion, and in artiular find the fixed oints in e η sae. We then show that tidal daming results in the system evolving towards the relevant fixed oint (now a quasi-fixed oint) on three times the tidal daming timesale. One this is ahieved the system ontinues to evolve toward the true fixed oint for whih both eentriities are zero. The timesale for this latter hase may be several orders of magnitude longer than the tidal daming timesale, and an in fat be muh longer than the age of the system. We illustrate the theory for hyothetial omanions to HD , and show that omanion masses as low as a fration of an Earth mass are aable of sustaining non-zero eentriites in the observed lanetary orbit for at least the age of the system, and in some ases at the level required to sustain the observed radii at the urrent eoh. Setion 4 mas out regions in m a arameter sae, where m anda are the mass and semimajor axis of ahyothetial omanion, for whih the eentriities of HD b and HAT-P1b would be at their observed level, as well as for WASP-1b whih for whih the eentriity is urrently unonstrained. Setion 5 alies the seular theory develoed here to determine whih onfigurations have irularized omanion lanets, while Setion 6 resents a summary and a disussion. 2 THE ECCENTRICITY-RADIUS RELATION Various authors have studied the roblem of the resonse of a lanet to internal and external heat soures (eg. Burrows et al. 2000, Bodenheimer, Lin & Mardling 2001, Bodenheimer, Laughlin & Lin 2003, Burrows et al. 2007). In artiular, models inlude lanets with and without ores, with and without stellar irradiation, and with and without an internal heat soure. In all these models the lanet assumes an equilibrium radius, R, and it is the relationshi between this and the rate of energy dissiation in the lanet, L, that is used in an attemt to exlain the range of lanetary radii observed in transiting lanets. In the ase that L is due to tides raised by an external omanion, the rate of internal energy dissiation in the lanet is given by (eg., Yoder & Peale (1981), Bodenheimer, Laughlin & Lin (2003)) ( )( ) Gm m e 2 Ė d =, (1) τ ir a to seond order in the eentriity, where τ ir e ė = 2 21n ( Q k ) (m m ) ( a R ) 5 is the irularization timesale of the orbit, with n the orbital frequeny, k the tidal Love number of the lanet and Q its Q-value or tidal quality fator (Goldreih & Soter 1966). 1 Equating this to the ower, L, needed to maintain an observed lanetary radius gives e in terms of the observable quantities m, R and a (m atually dros out), as well as the theoretial quantities Q, k and L : ( ) 5/4( ) m R 5/2( ) 15/4 ( ) a Q 1/2( ) ( 1/2 ) 1/2 k L e = 0.020, (3) M R J 0.04AU 10 5 k J 10 8 L where k J = 0.34 is the tidal Love number for Juiter. Burrows et al. (2007) alulate L for oreless solar metaliity lanets whih are irradiated by their arent stars, values for whih are listed in Table 2 under L Bu. Bodenheimer, Laughlin & Lin (2003) alulate L seifially for HD b only (in addition to various generi models), and find a value for L ten times higher for a oreless lanet, while for a lanet with a 20M ore, they find L is another fator of ten higher. Note that they take a slightly higher value for R (1.35R J). Values for e orresonding to L = L Bu and Q = 10 5 (see disussion below) are shown in Table 2 as well as values orresonding to L = 10L Bu. The estimated value of e for HD b (listed under e obs ) is more onsistent with L = L Bu than L = 10L Bu while the oosite is true for HAT-P-1b, although the eentriity (2) 1 Note that Goldreih & Soter (1966) use a modified Q-value, Q, whih absorbs the Love number suh that Q = 3Q /2k ; equation (4) in (Bodenheimer, Laughlin & Lin 2003) should have Q instead of Q. Another otential soure of onfusion is the use of the (quadruole) asidal motion onstant k 2 rather than the tidal Love number; these are related by k = 2k 2.

4 4 Rosemary A. Mardling Table 2. Theoretial data assuming Q = 10 5 Planet Age (Gyr) τ ir (Gyr) L Bu (L ) e (L Bu ) e(10lbu ) HD b 5.5±1.5 a ±0.009 WASP-1b 2.0±1.0 b HAT-P-1b 3.6±1.0 a ±0.02 eobs a Burrows et al. (2007) b Stemels et al. (2007) value for the latter awaits refinement. Note that in reality L is a funtion of R and m, however, its role in equation (3) is to arameterize various authors models. In both the models of Burrows et al. (2007) and Bodenheimer, Laughlin & Lin (2003), tidal energy is deosited homogeneously throughout the lanet s enveloe and dissiated loally. If, for examle, tidal energy is redominantly dissiated in the outer layers of the lanet, a signifiantly smaller value of L, and hene e, may be required to ahieve the same lanetary radius (as long as the energy is dissiated suffiiently dee below the surfae; Wu (2005) finds that tidal exitation of inertial modes results in energy being dissiated very lose to the surfae resulting in very little strutural adjustment). Alternatively, for a given e one would require larger values of Q for orresondingly smaller values of L. Note that this senario makes it diffiult to draw onlusions about the existene or otherwise of a lanetary ore. There is signifiant unertainty in the theory of tidal daming in stars and lanets, muh of it assoiated with internal struture and (in the ase of gaseous objets) the theory of turbulent dissiation (see Ogilvie & Lin (2004) for an exellent review). In the ase of solid lanets, one an infer daming timesales based on laboratory measurements of Q-values of solid materials, and for the seial ase of the Earth-Moon system, one an atually measure the rate of reession of the Moon. In the ase of Juiter and its satellite system, it is ossible to ut lower and uer bounds on Juiter s Q-value ( : Yoder & Peale 1981). These bounds are inferred from the existene of the Lalae resonane: if Q J were too low Io would have already assed through the 2:1 resonane with Euroa (robably destabilizing the system), while if it were too high it would never have been atured in the first lae. Similarly for Uranus, Q U < 39,000 (Tittemore & Wisdom 1989) while for Netune, 12,000 < Q N < 330,000 (Banfield & Murray 1992). There are various fators whih might be exeted to influene the Q-value of a hot Juiter and how this might differ from that of Juiter. One obvious differene between hot Juiters and our own Juiter are sin rates. While Juiter s sin eriod is 10 hrs, hot Juiters are aroximately synhronous with the orbital motion and hene have eriods around an order of magnitude longer. Goldreih & Soter (1966, 388) oint out that the smaller the mass of a lanet, the lower the Q-value. More reent work by Ogilvie & Lin (2004) and Wu (2005) who onsider tidal foring in slowly rotating lanets (aroriate to synhronously rotating short-eriod lanets) find that Q-values for suh lanets tend to be similar to that of Juiter. Wu (2003) shows that the observed orbital arameters of short-eriod extrasolar lanets indiate that their Q-values lie between and The lower value is inreased if one takes into aount the fat that lanets are likely to arrive at their urrent ositions with signifiantly larger radii that those observed today (see Fig. 1 of Bodenheimer, Lin & Mardling (2001)). Throughout this aer we use mostly Q = 10 5 to illustrate the theory develoed here. 3 LONG-TERM ECCENTRICITY EVOLUTION The estimated ages of HD , HAT-P-1 and WASP-1 are shown in Table 2 together with the irularization timesales of their lanets (assuming Q = 10 5 ). Clearly a mehanism is needed to maintain orbital eentriity if tidal heating is the soure of energy whih maintains the inflated lanetary radii. It seems to be imliitly assumed that the osillatory seular variations in the eentriities of a air of lose lanets (Murray & Dermott 2000) revail indefinitely, even in the resene of tidal dissiation. In this setion we show that the natural endoint of suh evolution is in fat a air of irular orbits, with the strutural roerties of the tidally fored lanet determining (at least in art) the timesale on whih this ours. This evolution ours in three stages (see Figure 3): (1) irulation of the angle between the lines of asides of the two lanets aomanied by the slow osillation of the eentriities at onstant amlitude, together with deline of the mean value of the inner lanet s eentriity until it reahes a (quasi)-fixed value: this ours on the irularization timesale; (2) libration of the angle between the lines of asides aomanied by the slow osillation of the eentriities with reduing amlitude but maintenane of the mean value of the inner eentriity: this ours on twie the irularization timesale; and (3) either (a): a slow non-osillatory deline in both eentriities to zero or (b): maintenane of the inner eentriity for some (often extremely long) eriod of time followed by its inrease to some maximum value followed by an often fast deline of both

5 Comanions to short-eriod lanets 5 eentriities to zero. Whether or not the final hase ours via route (a) or (b), it roeeds on a timesale sometimes several orders of magnitude longer than the irularization timesale. In general, the rates of hange of a lanet s orbital elements will have ontributions from omanions if they exist, the relativisti ost-newtonian ontribution to the otential of the star, the tidal and sin bulges of the star and lanet, and tidal dissiation in the atmosheres of the star and the lanet (Mardling & Lin 2002; Adams & Laughlin 2006a,b). For olanar systems, lanetary elements whih evolve are the semimajor axis, the eentriity and the longitude of eriastron. The only effet whih ontributes to a seular hange in the semimajor axis (exet when the system is in a resonane) is tidal dissiation, while for short-eriod systems the asidal motion is dominated by the third body and the relativisti otential of the star. The eentriity in turn is dominated by the third body and tidal daming. We will assume that the orbits are olanar, that the lanet is synhronously rotating with the orbital motion, and that the sin bulge of the star is negligible. For the analysis we will assume that the rate of hange of the lanet s semimajor axis is negligible, while numerial integrations will inlude this effet (see Rasio et al. (1996) and Mardling & Lin (2004) for disussions of this oint). The lan of this setion is as follows: (1) Study the equations governing the seular evolution of three olanar Newtonian oint masses for whih e is small omared to the omanion eentriity and m m, and determine the regions of arameter sae for whih the system librates or irulates. In artiular, determine the fixed oints of this dynamial system; (2) Determine how the results of (1) are modified by the inlusion of ost-newtonian terms in the otential of the star; (3) Introdue tidal daming and show that the system evolves towards doubly irular orbits on timesales desribed above. 3.1 Newtonian three-body seular evolution Consider a olanar oint-mass three-body system with masses m, m and m (for star, lanet and omanion resetively) for whih m m, but with no onstraints on m, and let the star-lanet and star-omanion semimajor axes be a and a resetively. 2 Further, let the orresonding eentriities be e and e, with e 1 but again with no onstraint on e, and let the longitudes of eriastra be and. While no energy is exhanged between the orbits of a (nonresonant) stable system (and hene the semimajor axes remain onstant), to first order in e the time-averaged rates of hange of e, e, and are (Mardling & Lin (2002); also see Mardling (2007a) for the full seular equations to all orders and inluding inlination and reession) ( de dt = 15 m )( a ) 4 16 ne sin( ), (4) m a (1 e 2 5/2 ) de dt = ne ( m m )( a a ) 3 sin( ) (1 e 2 ) 2, d dt and d dt = 3 4 n ( m m )( a a ) 3(1 e) 3/2 [1 5 4 = 3 4 n ( m m )( a a ) 2( 1 e 2 ) 2 [ ( ) ( ) ] a e os( ) a e 1 e 2 ( )( ) ( ) ] a e 1+4e 2 a e (1 e 2 os( ). (7) ) wheren andn aretheorbital frequenies(meanmotions)ofthelanetandomanionresetively.notethatmardling & Lin (2002) and Mardling (2007a) emloy a Legendre exansion in the ratio a /a to derive the rates of hange of the elements and make no assumtions about the mass ratios, eentriities or inlinations, in ontrast to the literal exansion used by Murray & Dermott (2000) whih relies on the smallness of the eentriities and inlinations but holds for all semimajor axes ratios. Equations (4) to (7) are given to otoole order (orresonding to the fourth ower in a /a ); see also Lee & Peale (2003) and Ford & Rasio (2000). (5) (6) Libration versus Cirulation of In order to derive simle onditions under whih the angle irulates or librates, we assume that e is signifiantly greater than e 1 (a situation relevant to the resent study; see Nagasawa, Lin & Ida (2003) for an alternative treatment). Then from (5), ė 0, and we an write ė = W oe sinη (8) 2 In fat, a Jaobian oordinate system is used to derive equations (4) to (7) so that the omanion s orbital elements are atually measured with reset to the entre of mass of the star-lanet system.

6 6 Rosemary A. Mardling and η = W q W o ( e e ) osη, where η and W q and W o are onstants (q for quadruole and o for otoole) defined by W q = 3 ( ) a 3 ( ) [ ( ) ] m 4 n ε 3 a m 1 ε 1 a m a m and W o = n ( a a ) 4 ( m m ) ε 5, where ε = 1 e 2. The fixed oints of (8) and (9) are (e,η) = (e W o/w q,2kπ) when W q > 0, (9) (10) (11) (12) and (e,η) = ( e W o/w q,(2k +1)π) when W q < 0. (13) Here k is any integer. Writing x e e (eq), where e (eq) = e W o/ W q = (5/4)(a /a )e ε 2 and linearizing about η = 0 we have, (14) 1 a /a (m /m )ε 1 ẋ = e W oη, η = or ( W 2 q W o ) x e, (15) η +W 2 q η = 0 (16) whih has solutions η(t) = A os(w qt+β), (17) where A 2 = η(0) 2 + ( 1 e (0)/e (eq) ) 2 and tanβ = ( ) 1 e (0)/e (eq) /η(0). (18) From (16) we onlude that the fixed oint η = 0 is stable. For libration around 2kπ, k 0, add 2kπ to (17). From (15) we have that e (t) = e (eq) [1 A sin(w qt+β)] e (lib) (t) (19) so that the maximum and minimum values of e when a system is librating are aroximately e max,min = (1±A)e (eq). For libration about η = π, also a stable fixed oint, (20) η(t) = π A os(w qt+β), (21) with A and β, e (t) and e max,min again given by (18), (19) and (20) resetively. For libration around (2k +1)π, k 0, add 2kπ to (21). Sine W q > 0 for libration about η = 0 and W q < 0 for libration about η = π, from (10) we have that m m > a a ε 1 for libration about η = 0; m m < a a ε 1 for libration about η = π. (22) Consider libration about η = 0. This an our as long as there are solutions to η = 0, that is, there exist real values for η when osη = (e /e )W q/w o. From (19) and (14) this is (in the linear aroximation) osη = 1 A sin(w qt+β) (23) whose minimum value is 1 A (taking A to be ositive) orresonding to the maximum value η ahieves during a libration yle. But in order for e to be ositive or zero, 0 A 1 from (20). Thus we have the result that libration about η = 0 an our as long as η < π/2. Moreover, from (20) the maximum value of e assoiated with A = 1, and hene the value assoiated with the boundary between libration and irulation, is

7 Comanions to short-eriod lanets 7 Figure 1. Evolution of e and η for an HD tye system (arameters given in text). The full seular equations were integrated to rodue these urves. (a) shows libratory and irulatory behaviour onsistent with e (eq) = and e boundary = 0.126, while (b) demonstrates that η is almost onstant for e (0) not too lose to e boundary. The dashed urve in both lots orresonds to e (0) = e boundary. The libration and irulation eriods are onsistent with the analytial value of yr. e boundary = 2e (eq) whih also holds for libration about η = π. In the latter ase libration an our as long as η π < π/2. These results are onsistent with numerially integrating the full equations (equations (4) to (7) but with the full deendene on e ; see Mardling (2007a)), an examle of whih is resented in Figure 1 for a system with m = 0.64M J, m = 1M J, m = 1.14M, a = AU, a = 0.4 AU and e = 0.3. The values for the inner lanet are aroriate for HD b. Figure 1(a) lots e against η for η(0) = 0 and 12 values of e (0) ranging from 0.06 to 0.2. Equation (14) gives e (eq) Figure 1(b) shows that for e > e boundary is not too lose to e boundary. Putting η = W q gives η(t) = W qt+η(0), so that from (8), ė = W osin(w qt+η(0)) and (24) = and e boundary = , η irulates with a frequeny that is aroximately onstant ( W q) as long as e e (t) = e (0)+(W o/w q)[os(w qt+η(0)) osη(0)] e (ir) (t). (26) Therefore if W q > 0, the maximum and minimum values of e orresond to η = 2kπ and η = (2k +1)π resetively, where k is an integer, and are given by e max,min = e (0)±e (eq) [1 osη(0)], while for W q < 0, these are e max,min = e (0)±e (eq) [1±osη(0)] (28) orresonding to η = (2k +1)π and η = 2kπ resetively. From (19) and (26) we see that both the libration and irulation eriods are given by 2π/W q. For the examle shown in Figure 1, 2π/W q = yr. The atual time variation of e an be obtained by forming de /de from (4) and (5), whih uon integration gives [ e (t) = 1 ( βe (t) 2 C ) ] 1/2 2, (29) where β = 1 2 a a ( m m ) and C = βe (0) 2 1 e (0) 2, (30) and e (t) is given by (19) or (26) aording to whether η librates or irulates. Summarizing these results, we have that libration will our about the fixed oint (e,η) = (e (eq),0) at frequeny W q > 0 when m a > ε 1 and 0 e < 2e (eq) and π/2 < η < π/2, (31) m a while libration will our about the fixed oint (e,η) = (e (eq) m m <,π) at frequeny W q < 0 when a a ε 1 and 0 e < 2e (eq) and π/2 < η < 3π/2. (32) (25) (27)

8 8 Rosemary A. Mardling Figure 2. Same as Figure (1) but with W GR inluded. Libration and irulation of η is more raid and the maximum amlitude for librating systems is redued in this ase. The dashed urves orresond to e = e boundary. The libration and irulation eriods are onsistent with the analytial value of 5700 yr. Similar results hold for k 0. If neither of the onditions (31) or (32) are met, η will irulate at frequeny W q with the maximumandminimumvalues ofe givenby(27)when m /m > ε 1 a/a (W q > 0),and (28)when m /m < ε 1 a/a (W q < 0). 3.2 Relativisti Effets Contributions to the otential of the star from ost-newtonian terms affet the asidal advane only. The orbit-averaged ontribution to η from the innermost orbit is (eg. Mardling & Lin, 2002) ( W GR = 3n na ) 2, (33) 1 e 2 where is the seed of light, so that η = W q +W GR W o ( e e ) osη W q W o ( e e ) osη, (34) where W q = 3 ( ) 3 ( ) [ a m 4 n ε 3 a m 1 ( ) ] a m ε 1 +γε 3, (35) a m and γ = 4(n a /) 2 (m /m )(a /a ) 3 is, to first-order in e, the ratio of W GR to the quadruole ontribution to. One effet of W GR on the system is to shorten the eentriity modulation and asidal eriod (omare Figure 2 with Figure 1). In this examle, W GR W q in ontrast to Merury s orbit around the Sun for whih W GR is only 7% of the ontributions to W q from the other lanets in the Solar System. The equilibrium eentriity is now given by 3 e (eq) = e W o/ W q = (5/4)(a /a )e ε 2 1 a /a (m /m )ε 1 +γε 3. (36) Thus when a /a (m /m )ε 1 < 1, the effet of W GR is to derease the equilibrium eentriity, an examle of whih is shown in Figure 2, while for a /a (m /m )ε 1 > 1 the oosite is true. Mardling & Lin (2004) show that the former an be vital for the survival of short eriod terrestrial lanets with omanions. The analogues of (31) and (32) are that libration will our about the fixed oint (e,η) = (e (eq) where e (eq) m m >,0) at frequeny W q > 0, is now given by (36), when a ε 1 /(1+γε 3 a ) and 0 e < 2e (eq) and π/2 < η < π/2, (37) while libration will our about the fixed oint (e,η) = (e (eq) m m <,π) at frequeny W q < 0 when a a ε 1 /(1+γε 3 ) and 0 e < 2e (eq) and π/2 < η < 3π/2. (38) 3 Note that equation (4) in Mardling & Lin (2004) is in error by a fator of 2.

9 Comanions to short-eriod lanets 9 If neither of the onditions (37) or (38) are met, then when m /m > ε 1 a/a /(1+γε 3 ), η will irulate at frequeny W q in the ositive diretion with the maximum and minimum values of e given by (27), while when m /m < ε 1 a/a /(1+γε 3 ), η will irulate in the negative diretion at the same frequeny with maximum and minimum values of e given by (28). The time variations of η, e and e are given by the equations set out in the revious setion, with e (eq) given by (36). The form of (36) suggests that for a range of arameters, e (eq) an be large. Singular values of e (eq) exist for some ombinations of arameters, although when the otoole term in (7) is inluded in (9) these are relaed by finite, although still relatively large values (see Nagasawa, Lin & Ida(2003)). As we will show in Setion 3.3.4, this has imortant onsequenes for eentriity foring by low-mass omanions. 3.3 Tidal daming and quasi-fixed-oint behaviour While tidal dissiation in the innermost lanet diretly affets the rate of hange of its eentriity, it also indiretly affets the rate of hange of the omanion lanet s eentriity (reall that the revious analysis assumed e was onstant), as well as the angle between the asidal lines, η. 4 The relevant system of equations now beomes ė = W oe sinη W T e, ė = W e sinη, η = W q W o ( e e ) osη, (39) (40) (41) where from (2), W T = 21 2 n ( k Q )( m m )( R a ) 5 = τ 1 ir (42) and from (7), W = n ( m m )( a a ) 3ε 4, where again we have ignored the otoole term (roortional to e /e ) in (7). The damed autonomous system (39)-(41) evolves to the formal fixed oint (e,e,η) = (0,0,η ), where η is given below and is suh that η (mod2π) 0 or η π (mod2π) 0 for W q > 0 and W q < 0 resetively. However, there are three timesales assoiated with this evolution: (1) The system irulates until e, e and η are suh that libration onditions are met (see (31) and (32): the system is atured by the assoiated fixed oint). This ours on the irularization timesale of the inner orbit, 5 τ ir; (2) The system librates about and evolves towards the (now quasi-)fixed oint (e,e,η) = (e,e,η ), where e e (0) and e e (eq) (equation (36)), on the timesale 2τ ir (e and e will be roerly defined below); (3) The system evolves smoothly towards the formal fixed oint (e,e,η) = (0,0,η ) on some timesale τ τ ir. The details for these hases are as follows, where for the first two we assume that ė ė so that e e (0). (43) Cirulation Following Setion 3.1.1, the irulation hase an be aroximated with η = W q. Then (39) beomes ė +W T e = W oe sin(w qt+η(0)) (44) whih has solution e (t) = Ce W Tt +σe (eq) os(w qt+η(0)) (45) where σ = sgn [ ] W q, C = e(0) σe (eq) osη(0) and we have assumed that W T W q. This hase is shown in blue in Figure 3 and evolves on the timesale 1/W T = τ ir. Notie that the amlitude of variation is not affeted by dissiation and is given by e (eq). Also note that Q is set artifiially low in order to learly demonstrate how this hase roeeds. 4 The ontributions to η from the tidal and sin bulges of the star and lanets are negligible omared to the terms inluded here. 5 This hase is skied if the initial onditions are suh that the system starts librating immediately.

10 10 Rosemary A. Mardling Libration The quasi-fixed oint (e,e,η ) is obtained by utting ė = 0 and η = 0 in (39) and (41) resetively. While there is only one true equilibrium value of e (that is, e = 0), we take e to be the value of e when the libration amlitude of e dros below some seified small value (see next setion). Eliminating e /e between (39) and (41) then gives sinη = WT W q osη WT W q < 0, (46) where the seond ste omes from the fat that sgn(osη ) = sgn(w q) and W T/ W q 1. Thus for W q > 0, η (mod2π) W T/W q while for W q < 0, (η π)(mod2π) W T/W q. For HD b with a 0.1M J omanion at 0.4 AU, η = radians, deending on the value of e (taking Q = 10 5 ). The value of e is given by (36) with e = e. Linearizing about (e,e,η ) and utting x = e e and y = η η rodues the system ẋ = e W oy W T x, ẏ = (W 2 q/e W o)x, (47) (48) whih yields the damed librating solution e (t) = e +Ae W Tt/2 os(w qt+δ) (49) and η(t) = η Ae W Tt/2 sin(w qt+δ). (50) For the uroses of this analysis, the onstants A and δ an be taken to be suh that the librating hase is entered at the first ourrene of e = 0. Writing the time at whih this ours as t = t lib, (49) and (50) beome e (t) = e { 1 e W T (t t lib )/2 os[w q(t t lib )] } (51) and η(t) = η e W T(t t lib )/2 sin[w q(t t lib )] (52) In ontrast to the irulating hase, the amlitude of osillation of e (t) diminishes on a timesale 2/W T = 2τ ir as the system aroahes the quasi-fixed oint. Figure 3 shows this behaviour in red. For the analysis above we have assumed e is onstant. However, as Figure 3() illustrates, e slowly delines and has no analogue of e (eq) (exet e = 0). The following setion studies the effet of this on the long-term evolution of the system Long-term evolution One the quasi-fixed oint is reahed, the omanion s eentriity evolves aording to (40) with e = e W o/ W q and sinη = W T/ W q : ė = WWoWT ( 25 = 16 = e W 2 q )( m )( ) 5/2 a WT m a ( ) λ e τ ir F(e ) where ( )( )( ) 25 5/2 m a λ = 16 m a and e F(e ) F(e ) = ε 3 (1 αε 1 +γε 3 ) 2 (54) with α = a /a (m /m ) (reall that ε = 1 e 2 ). Note that α and γ are effetively onstant throughout the evolution sine the semimajor axes evolve on timesales muh longer than the age of the system (see Figure 5(d)). One an obtain an aroximate timesale on whih e (and hene e ) evolves by utting ε ε(e ) ε, where e is yet to be formally defined. Equation (53) then gives (53) e (t) = e (t )e (t t)/τ, (55) where τ = (F(e )/λ)τ ir (56)

11 Comanions to short-eriod lanets 11 Q = π ϖ η, ϖ, ϖ ϖ η= ϖ ϖ 14 π 0 Figure 3. Cature hase for HD b with a 0.1M J omanion at 0.4 AU with e = 0.4. The full seular equations inluding evolution of a and a are integrated (Mardling & Lin 2004). The Q-value of the observed lanet is set artifiially low at Q = 100 in order to illustrate the roess (timesales are linearly roortional to Q and values for e are not very sensitive to the value of Q; for this examle, the relative hange in e is less than 2% over four orders of magnitude in Q ). Blue urves orresond to the irulation hase, red to the libration hase and green to the slow evolution or quasi-fixed oint hase (see text for more disussion). In anel () notie that the relative hange in e is small, justifying (at least for this examle) taking e onst during this hase in the analysis. The diamond symbols orresond to the oints (t,e ), (t,e ) and (t,η ) in anels (b), () and (d) resetively. Panel (d) shows the asidal advane of the inner and outer orbits, with loking onto. This ours irresetive of the mass of the outer body beause e /e 1 (see equation (7)). and t orresonds to the beginning of the slow evolution hase whih we hoose to be the first instane of η η < ǫ for all t > t. Here ǫ 1 is an arbitrarily hosen arameter and is taken to be ǫ = 10 4 for the studies resented in this aer. This then serves as definitions of e and e, that is, e = e (t ) and e = e (t ). The transition to the slow hase at time t is indiated in Figures 3(b)-(d) by a diamond symbol. Note that τ τ ir for most stable systems, exet when F(e ) is small (see Setion 3.3.4). The aroximate solution (55) is the blue dot-dashed urve in Figure 4(a). Keeing the full deendene on e in (53) and integrating gives f(e ;α,γ) f(e ;α,γ) = λ(t t )/τ ir, where f(x;α,γ) = 1 3 γx γx4 +(α 3γ)x x 2[ α 2 42 γ( 3x 4 11x ) ( α x 2 4 ) +γ 2 g(x) ] +(γ α+1) 2 lnx [ (γ α) 2 +1 ] ( ) ln 1+ 1 x 2 (57) (58) with g(x) = ( 35x 8 185x x 4 506x ) /315. (59) This solution is shown as a dashed red urve in Figure 4(a) and is learly suerior to the estimate (55) whih underestimates the timesale on whih the system evolves one asidal loking has taken lae. We define the true timesale on whih irularization of the omanion orbit takes lae (the e-folding time) to be the time orresonding to e /e, that is, from (57), τ true = [ 3+λ 1 (f(e ;α,γ) f(e /e;α,γ)) ] τ ir, (60)

12 12 Rosemary A. Mardling Q = 100 Q = 100 Figure 4. Continuation of numerial solution in Figure 3 showing the long-term evolution of (a) e and (b) e (blak urves, sitting underneath red urves), again with Q = 100 for illustration (note the y-axis sales are different to those in Figure 3). This orresonds to 125 Gyr of evolution when Q = 10 5 (omare this to τ ir = Gyr). Also lotted in (a) is aroximate solution (55) (blue dot-dash urve) and the analyti solution (58) (dashed red urve). The aroximate solution underestimates the timesale on whih e evolves by a fator of 1.8. Also lotted in (b) is the analyti exression for e (t) (red urve) using (36) with e = e (t) from (58). The analyti solutions are essentially indistinguishable from the numerial solutions. where we have ut t = 3τ ir (1 for the irulation hase + 2 for the libration hase). As e slowly dereases, e maintains a value given by (36). Any dearture from this quasi-equilibrium value is quikly damed out on the irularization timesale τ ir. Thus we an use (57) and (36) to alulate analytially e (t) and e (t) during the slow hase, and hene to determine the urrent value of e using the estimated age of the system. The analyti solution for e (t) is shown in Figure 4(b) (red urve) and is indistinguishable from the numerial solution The inreasing e regime: very low-mass omanions The solution for the examle given in Figure 4 is essentially a modified deaying exonential. However, there is another regime of the funtion f(x; α, γ) whih behaves very differently and is assoiated with very low-mass omanions whih, for a large art of the life of suh systems, fore e to inrease to some maximum value before deaying to zero at the end of the tidally evolving hase. An examle of this is shown in Figure 5 for the HD system with a omanion of mass m = M J = 4.2M at a = 0.3 AU and an initial omanion eentriity of 0.7, a high value hosen to learly illustrate the theory. The system enters the slow hase with e = 0.705, e = and η 43π. The average reflex veloity of the star due to the omanion is 0.64ms 1. See the figure ation for a detailed desrition of the evolution of the system. The red dashed urves again reresent the analyti solutions (58) and (36), learly demonstrating their reditive ower. In fat, the analyti solution for e asses through a singularity as the denominator in (36) asses through zero. The full numerial solution avoids this, but still manages to inrease to a signifiant level on a reasonably long timesale. For Q = 10 5, assage through the singularity ours at 2.2 Gyr, well before the system age of 5.5 Gyr, while for Q = , assage through this oint ours now. Note that this system is stable when the full three-body (non-seular) equations are integrated inluding ost-newtonian relativisti terms, tidal and sin bulges, as well as tidal daming (Mardling & Lin 2002). Although systems whih ass through the singularity suffer a raid deline in e one they have done so, they tend to maintain the value of e that they had at the beginning of the slow hase before this haens (see also Figure 6). This is in ontrast to systems whih don t ass through the singularity; these simly deay aroximately exonentially. To determine the arts of arameter sae for whih e neessarily asses through the singularity and therefore maintains e before it does so, and in artiular whih systems ass through the singularity after the urrent age of the system, we start by studying the behaviour of the denominator of (36). Defining (e ;α,γ) = 1 αε 1 +γε 3, (61) where (again), α = a /a (m /m ), Figure 7 lots this funtion against m for seifi values of a and for e = 0.8 and e = 0. There are three regimes of interest. As the system evolves from e down to e = 0, either (a) (and hene W q) remains negative so that osη 1 always (blue arrow); (b) asses through zero so that η goes through a short eriod of irulation as osη asses from -1 to 1 (yellow bakground, red arrows) or () is ositive always and osη 1 (green arrows). In the ase shown here (e = 0.8), systems ross = 0 for a finite range of values of m when a < 0.29; when 0.29 < a < 0.52, all systems with m less than some ritial value ross = 0; when a > 0.52 no systems ross = 0. Figure 7(b), for examle, shows that all systems with m < 0.25M J and a = 0.3 (and e = 0.8) neessarily evolve through

13 HD : a = 0.3, e = 0.7, m = 4.2 M, = 100 Comanions to short-eriod lanets 13 Q e e η 40 a π π Figure 5. A terrestrial-mass lanet an sustain the observed eentriity of HD but not for long enough in this examle where m = 4.2M, at least for Q < Reall that the atual timesale is roortional to Q so that 6 Myr with Q = 100 orresonds to 6 Gyr with Q = The blak urves show the numerial solutions while the red dashed urves are given by the analytial solutions (36) for e and (58) for e. (a): The evolution of e following ature. After effetively maintaining the value e for some 2 Gyr (for Q = 10 5 ), e inreases to a maximum of around 0.04 before raidly delining to around and then slowly delining to zero. The inrease in e orresonds to a singularity in the analyti solution whih fairly aurately follows the full solution. The fat that the analytial solution is reliable onfirms this deendene. Thus, for Q = 10 5, the maximum value of e ours around 2.7 Gyr. Reall also that the value of e at the beginning of the slow hase is quite insensitive to e (0), η(0) and Q. (b): The evolution of e following ature. Again the analyti solution is a good aroximation to the full solution. The blue dot-dashed urve is the aroximate analyti solution given by (55). Thus when a system must ass through the singularity, the timesale given by (56) overestimates the irularization timesale in ontrast to the revious examle. (): The evolution of η =. During the re-ature hase, η irulates bakwards before being atured at η = 43π (reall that when W q < 0, osη 1). As the system asses through the singularity, that is, as W q asses through zero to beome ositive (see equation (35)), η irulates bak u to 40π so that osη 1. The brief irulatory eriod is aomanied by osillations in e and e whih an be seen in (a) and (b). (d): The evolution of a. The full (seular) numerial solution inludes the evolution of both semimajor axes; here we see that over the entire evolution of this system, a hanges by less than 4%. = 0. The ranges above vary with different values of e ; urves orresonding to e < 0.8 lie between the two shown. What remains to be determined, however, is whih of these systems maintain a non-zero value of e for at least as long as the urrent age of the system before they do so. One this is done, the arameter ranges should be further restrited by determining whih systems are stable (Mardling 2007b). Figure 8 lots e (now), e (t ) and e (now) against omanion mass for the HD system for various values of a and e. Here now refers to t = 5.5 Gyr, and the quantities are alulated using (36) and (58). The green-filled irle orresonds to the system shown in Figure 5 at t = 5.5 Gyr for Q = 10 5 (t = 5.5 Myr for Q = 10 2 ), while the red-filled irles orresond to Figure 6 at the same time. Note the extremely low omanion masses in Figures 6(a) and (b): m = 1.9M and m = 0.4M resetively. The maximum values of e ahieved are (a) at 5.8 Gyr and (b) at 6.0 Gyr. The orresonding values of e are 0.27 and At 5.5 Gyr, the values for e are (a) and (b) 0.004, with orresonding values for e of 0.41 and Thus senario (a) reresents a ossible solution for HD for a value of Q around 10 5,

14 14 Rosemary A. Mardling HD : e (0) = 0.7, Q = 100 a = 0.3 AU, m = 1.9 M a = 0.35 AU, m = 0.4 M Figure 6. Similar to Figure 5 but with lower values for m. Passage through = 0 ours later in the life of the system (omare anel (a) to Figure 5(a)), while e is smaller. For (a), the maximum value of e = is reahed at 5.8 Gyr for Q = 10 5 (ie., 5.8 Myr for Q = 100), while at 5.5 Gyr, e = For (b), e = at 5.96 Gyr, while e = at 5.5 Gyr. e *= 0 e*= 0.8 Figure 7. Behaviour of the denominator of (36),, used to determine whih onfigurations have e (eq) neessarily assing through a singularity during their evolution. See text for disussion. Note that some of the systems reresented here are not stable when the full three-body equations are integrated with aroriate daming (Mardling & Lin 2002). while senario (b) requires a Q-value around 10% smaller in order that the maximum value of e ours about now (reall that this also deends on the age estimate for the star). Figure 8 demonstrates that high initial values of e are required in order that extremely low-mass omanions be resonsible for maintaining the eentriity of inflated lanets (see Setion 6 for a disussion on this oint). Figure 8, together with Figures 5 and 6 and equation (60), suggest that for given values of a and e, the time a system takes to ahieve its maximum value of e (given it asses through = 0) inreases as the omanion mass dereases. Figure 9 lots this time, T evol, as a funtion of m for various values of a and e (solid urves). Also lotted is τ true (equation (60)) whih is the e-folding time for systems whih don t go through = 0 but simly deay aroximately exonentially (dashed urves). The small values of T evol orresonding to, for examle, systems for whih e = 0.5 and a = 0.34 with m < 0.002M J (the dark blue urve whih aears oinident with the x-axis), omared with the large values of T evol for m > 0.003M J, simly indiate that assage through = 0 ours early in the life of the system, while deay following assage is very slow. Figure 10 illustrates this with the above examle. The blue diamond indiates the nominal beginning of the slow hase (as defined in the disussion following equation (56)). There is little differene between Figures 10(a) and (b) exet that 10(b) does not meet the slow hase riterion until after the system asses through = 0. It is of interest to ask whether or not it is ossible for a system to ski the librating hase and evolve diretly to the true fixed oint (e,e ) = (0,0). Referring to equations (37) and (38) we see that for e < 2e (eq) a system will librate as long as π/2 < η < π/2 or π/2 < η < 3π/2. Sine a irulating system must ass through all values of η, one e is ermanently below 2e (eq) it has no hoie but to librate until a loal quasi-fixed oint is reahed.

15 Comanions to short-eriod lanets 15 HD : age = 5.5 Gyr, Q 5 = 10 e (now) e (now) e (now) e (now) Figure 8. The values of e and e (blak solid and red-dash urves resetively) at t = 5.5 Gyr, and the value of e at the beginning of the slow hase (e (t ): blue dot-dash urves), all with Q = From Figures 5 and 6 it is evident that systems whih ass through = 0 in their lifetime aroximately maintain their value of e at t = t before they do so. Thus, for examle, in (a) all systems with m < 0.004M J will have the same finite (but small) eentriities they had when they entered the slow hase. For a small range of masses (0.004 < m /M J < 0.006) e (now) > e (the blak urve sits above the blue urve), while for greater masses e(now) e. Note that formally these urves have no uer bound at the oint orresonding to = 0, however the finite resolution of the lots suggest erroneously that they do. In reality there does exist an uer bound when the seular equations are diretly integrated as is evident in Figures 5 and 6. It is lear from this figure that high initial values of e are required to have very low-mass omanions resonsible for maintaining the eentriity of inflated lanets, although lower values of e are required for lower values of a (see Figure 13 and Setion 6 for a disussion on this oint). 4 CONSTRAINTS ON ORBITAL PARAMETERS OF COMPANIONS We are now in the osition to onstrain the orbital arameters and masses of ossible omanion lanets to HD b, HAT-P1b and WASP-1b. Given estimated ages, our aim is to find initial onditions and omanion masses for whih these lanets have eentriities in the ranges observed (WASP-1b is urrently unonstrained), and in the ase of HD , rodues a stellar reflex veloity less than 10ms 1 (an uer limit whih gives omanion masses less than the uer bound of 0.3M J ited in Laughlin et al. (2005a)). Aknowledging again that these systems are unlikely to be exatly olanar (but see Figure 15 whih shows that moderate inlination has little effet), we assume zero mutual inlination and vary the arameters a, m and e (0). Our results orresond to the ase e (0) = 0.1 and η = 170 o, but are not very sensitive to using other values of these arameters (ie, they evolve to the quasi-fixed oint (e,e,η ) initially with a variation in e of at most around 10%). Note also that while our aim is to determine whether or not urrent-day solutions exist (eseially for HD whih is the most onstrained), we an also say something about the orbital elements at earlier times and hene (autiously) onstrain the initial onditions.