+m_Pm_T}}{(m_S + m_P + m_T)^{2/3}\sqrt{1-e^2}},
\label{eqn:3body}$$ is equal to the observed value. Hence we simulate the case in which the observed RV variation is due to the combined force of a planet and a Trojan, rather than a planet alone, but the RV data alone are insufficiently precise to discern the difference.[^2] We compute the transit time $t_C$, repeat the analysis for an increasing sequence of $m_T$, and fit a polynomial function to the resulting relationship $t_C(m_T)-t_C(0)$. We found a quadratic function, $m_T = a_1 \Delta t + a_2 (\Delta t)^2$, to give a good fit to the results. Taking $m_T$ to be in Earth masses and $\Delta
  t$ in minutes, the coefficients ($a_1$,$a_2$) are (0.044,-1.17$\times$ 10$^{-5}$) for GJ 436b and (6.787,-0.001) for XO-3b. For the cases of HAT-P-2b and HD 17156b, we find that even very low-mass Trojan companions are dynamically unstable, owing to the large orbital eccentricities (see § \[subsec:dynamical\]). Thus, for those systems, the requirement of dynamical stability is more constraining than the empirical upper limit on $m_T$ based on the FG method. (As will be described in § 4, this also proved to be true for XO-3b based on the current data.)

Data Analysis {#sec:data}
=============

The RV data were taken from the available literature on each system. The references are given in Table 1. These data were generally obtained for the purpose of discovering or confirming the planet, although in a few cases the data were obtained for other reasons, such as precisely measuring the orbital eccentricity (Laughlin et al. 2005) or for measuring the Rossiter-McLaughlin effect (Winn et al. 2006). Regarding the latter, the data that were obtained while a transit was in progress were not used, to avoid the needless complication of incorporating the Rossiter-McLaughlin effect into the RV model. However, in those cases the investigators usually gathered additional data
+ m_Pm_T}}{(m_S + m_P + m_T)^{2/3}\sqrt{1 - e^2 } }, 
 \label{eqn:3body}$$ is equal to the observed value. Hence we simulate the subject in which the ascertained RV variation is due to the compound violence of a planet and a Trojan, rather than a satellite alone, but the RV data alone are insufficiently accurate to spot the difference.[^2 ] We compute the transit prison term $ t_C$, repeat the analysis for an increasing succession of $ m_T$, and fit a polynomial function to the result relationship $ t_C(m_T)-t_C(0)$. We found a quadratic affair, $ m_T = a_1 \Delta t + a_2 (\Delta t)^2 $, to give a good fit to the resultant role. Taking $ m_T$ to be in Earth masses and $ \Delta 
   t$ in minutes, the coefficients ($ a_1$,$a_2 $) are (0.044,-1.17$\times$ 10$^{-5}$) for GJ 436b and (6.787,-0.001) for XO-3b. For the casing of HAT - P-2b and HD 17156b, we find that even very low - mass Trojan companions are dynamically unstable, owe to the large orbital eccentricities (see §   \[subsec: dynamical\ ]). Thus, for those systems, the requirement of dynamical stability is more constraining than the empirical upper limit on $ m_T$ based on the FG method. (As will be described in §   4, this also proved to be true for XO-3b establish on the current data .) 

 Data Analysis { # sec: data } 
 = = = = = = = = = = = = = 

 The RV datum were accept from the available literature on each organization. The references are given in Table   1. These datum were generally prevail for the purpose of discovering or confirm the planet, although in a few cases the data were obtained for early reasons, such as precisely measuring the orbital eccentricity (Laughlin et al.   2005) or for measuring the Rossiter - McLaughlin effect (Winn et al.   2006). Regarding the latter, the data that were obtained while a theodolite was in progress were not used, to avoid the needless complication of incorporate the Rossiter - McLaughlin effect into the RV model. However, in those cases the detective usually gathered additional data
+m_Pm_H}}{(m_S + m_P + m_T)^{2/3}\sqrt{1-e^2}},
\label{eqn:3nody}$$ is equal to the obverved value. Hdnce we simulate the case in wyich uke observed RV variagion is dle to the comuined force of a planet and a Fvojan, cather than a pkanet alona, but the RV ddtx clone are insufficiently precise to qiscern tje difference.[^2] Re cpiputs the transit time $t_C$, repeat the ahalysis for an incressing sequence of $m_T$, and flt a polynomial functiln to the rgauleung relationrhip $t_C(m_T)-t_B(0)$. We found z quadratic function, $m_T = a_1 \Deltx t + c_2 (\Delta t)^2$, ti tivf a good fit to tre results. Taking $m_T$ to be on Earth massex aid $\Dwlta
  t$ in minutes, thx coefficients ($a_1$,$a_2$) arg (0.044,-1.17$\times$ 10$^{-5}$) fmr GJ 436b and (6.787,-0.001) for XO-3v. Dor tve cdses if FAT-K-2b end HD 17156b, ae hind that eben very loq-mass Trojan companoogw are dynamiczlly ugseable, owing to the large orbital eccentgicifies (see § \[subsec:dynamicao\]). Thus, for those systgms, the rezuirement of dynamical stability is more constraiting vhxn uhc dnplrical upper limit on $m_T$ based on the FG methjs. (Ss will be descrlbed in § 4, this also pgofgd to be true wor XO-3y bzsed on the currenh data.)

Dwta Abalysis {#stc:dats}
=============

The RV data were taken fron the availayle literature on eack system. The refgrencex are given in Table 1. Thzse dafa were genfrally obfxined for the pufpoxe of discovering or confirmyng the poanec, althoueh im a fer cases thf data were obtained for othet reasmns, such ad precisely measuring the orbitek eccentricijy (Nauchlin et al. 2005) ov for measuring the Rossiter-MeLaughliu effezt (Winn et al. 2006). Regerding the lwtter, the datd that were outained wrile a teansit dxs in progress were not used, to acoid the needless gomplkdation of incor'utqting the Rossotef-McJalghkig effect into the RV kodel. However, in ghosr cases the investigdtora usually gathered aqditionao data
+m_Pm_T}}{(m_S + m_P + m_T)^{2/3}\sqrt{1-e^2}}, \label{eqn:3body}$$ is the value. Hence simulate the case variation due to the force of a and a Trojan, rather than a alone, but the RV data alone are insufficiently precise to discern the difference.[^2] compute the transit time $t_C$, repeat the analysis for an increasing sequence of and a function the resulting relationship $t_C(m_T)-t_C(0)$. We found a quadratic function, $m_T = a_1 \Delta t + a_2 t)^2$, to give a good fit to the Taking $m_T$ to be Earth masses and $\Delta t$ minutes, coefficients ($a_1$,$a_2$) (0.044,-1.17$\times$ for 436b and (6.787,-0.001) XO-3b. For the cases of HAT-P-2b and HD 17156b, we find that even very low-mass Trojan companions dynamically unstable, the large eccentricities § Thus, for those requirement of dynamical stability is more empirical upper limit on $m_T$ based on the method. (As be described in § 4, this proved to be true for XO-3b based on current data.) Data Analysis {#sec:data} ============= The RV data were taken from the available literature system. The references are in Table 1. data generally for purpose of or confirming the planet, although in a few cases the data obtained for other reasons, such as precisely measuring the orbital et 2005) or for the Rossiter-McLaughlin effect (Winn al. Regarding the latter, the were while in were used, to avoid the complication of incorporating the Rossiter-McLaughlin into the RV model. investigators usually gathered additional data
+m_Pm_T}}{(m_S + m_P + m_T)^{2/3}\sqrt{1-e^2}},
\label{eqn:3boDy}$$ is equal tO the oBseRveD vAlue. hencE we simulate the CAse iN which the observed RV varIatioN iS Due tO ThE combIned forCE oF A PlaNeT aNd a trOJaN, rathEr tHan a plaNet alone, buT thE Rv data alone arE InSufficientLy pRecise to discErn The difFeRenCE.[^2] We coMpuTe the TransiT Time $t_C$, Repeat the AnALysis fOR an incrEASiNg seQuence of $m_T$, and fit a POlYNomial function To the rEsULtING reLatIonship $t_C(m_t)-t_c(0)$. We foUNd a quadRAtIC FUncTIon, $m_T = a_1 \Delta t + a_2 (\delta t)^2$, to givE A goOd fit tO tHe rESults. TAking $M_T$ TO be In Earth massEs anD $\Delta
  t$ in MinuteS, The coefFIcients ($A_1$,$a_2$) are (0.044,-1.17$\tImeS$ 10$^{-5}$) foR GJ 436b ANd (6.787,-0.001) FoR XO-3B. FOR thE CaSes OF HAt-P-2b and HD 17156B, wE fInd thAt evEN VERy loW-maSs TrOjan cOmpanions are dYnaMicaLLy uNstabLe, owiNg to ThE largE orbitAl eccEnTricities (see § \[subSec:dYnamical\]). THus, FoR thOsE systEMs, the rEquIreMent of dYnamicaL StaBiLITY iS more constraining tHaN THe EmpiricaL upper LImIt ON $m_T$ based On The fG meTHOd. (As wIll bE DeScribed iN § 4, this aLSo PrOved to bE tRue for xO-3B baSed On the CUrreNt data.)

data AnalYsis {#sEC:data}
=============

The RV data WEre taken from tHE aVAIlABle lIteRature on eacH sysTEm. ThE refEReNceS Are giVen in taBLe 1. tHese data were generalLy ObtainEd for The purpose of dIscovering OR COnfirminG the PLaNEt, although in a fEw casEs the data wERe obtainEd for Other reaSons, such aS PRecisely MeaSurIng The ORBiTal eccentriciTY (laugHlIn et al. 2005) oR foR measurIng The rosSitEr-mcLaughliN effect (WInN eT aL. 2006). REgaRding THe latter, ThE daTa ThaT were OBtaineD whilE a trAnSiT Was In progrESs WERe noT uSeD, to aVoiD tHe neeDlesS ComPlicatiOn of incorPorATing ThE ROssiter-mcLaughlin effEcT into the RV MoDel. howeveR, IN those caSes the investigators usuaLLy gatheRed AdditIonaL data
+m_Pm_T}}{(m_S + m_P + m_T )^{2/3}\sq rt{1- e^2 }}, 
\ labe l{eq n:3body}$$ ise qual  to the observed value . Hen ce we s i mu latethe cas e  i n whi ch  t heob s er ved R V v ariatio n is due t o t he  combined fo r ce  of a plan etand a Trojan , r atherth ana  plan etalone , butt he RVdata alon ea re ins u fficien t l yprec ise to discern th e  d i fference.[^2]We com pu t et h e t ran sit time $ t_ C$, r e peat th e  a n a l ysi s  for an incre asing seque n ceof $m_ T$ , a n d fita pol yn o mia l functionto t he result ing re l ationsh i p $t_C( m_T)-t _C( 0)$ . We fo un d a  q u adr a ti c f u nct ion, $m_ T=a_1 \ Delt a t + a_ 2 ( \Del ta t) ^2$, to givea g oodf itto th e res ults .Takin g $m_T $ tobe  in Earth masse s an d $\Delta 
 t$  in  m inute s , thecoe ffi cients($a_1$, $ a_2 $) a r e(0.044,-1.17$\time s$ 1 0$ ^{-5}$)for GJ 43 6b and (6.7 87 ,-0 .001 ) for X O-3b .  F or the c ases o f  H AT -P-2b a nd  HD 17 15 6b,  we  find that  evenvery low -mass Trojan compani o ns are dynami c al l y  u n stab le,  owing to t he l a rgeorbi t al  ec c entri citie s( se e  § \[subsec:dynamic al \]). T hus,for those sys tems, ther e q uirement  ofd yn a mical stabilit y ismore const r aining t han t he empir ical uppe r limit on  $m _T$  ba sed o nthe FG method . (Aswi ll be d esc ribed i n §  4,  th isal so proved  to be t ru efo rXO- 3b ba s ed on th ecur re ntdata. ) 

Data  Anal ysis  { #s e c:d ata}
== = == = = ==== == 
TheRVda ta we re t a ken  from t he availa ble lite ra tu re on e ach system. T he  reference sare  given i n Table1. These data were gene r ally ob tai ned f or t he purpos e o f disc ove r ing or  confi rming  t hep l anet, a lt hou gh  in a fewc a ses  theda ta w ere obt ained for other re a son s, such as pr eci sely m ea sur i ng the  o r bit a l  eccentricity ( Laughlin e ta l.  2005) orf orme asuring  the Ro ssite r -McLaug hlin effe ct (Winnet  al. 2 006 ). Regardi ng the l atter, th e  data th at we reobtain ed  wh ile a  trans i t w as in  progr es s were  notus ed, to a void the needless compl icatio n ofinc orporatin g t h e R ossiter-M cLau ghlin effe ctint o the  RV model . Ho w ev er, in th osec ases thei nv est i g at ors usually g a the red a ddi t ionaldata 
+m_Pm_T}}{(m_S +_m_P +_m_T)^{2/3}\sqrt{1-e^2}},
\label{eqn:3body}$$ is equal to_the observed_value._Hence we_simulate_the case in_which the observed_RV variation is due_to the combined_force_of a planet and a Trojan, rather than a planet alone, but the RV_data_alone are_insufficiently_precise_to discern the difference.[^2] We_compute the transit time $t_C$,_repeat the_analysis for an increasing sequence of $m_T$, and_fit_a polynomial function_to the resulting relationship $t_C(m_T)-t_C(0)$. We found a quadratic_function, $m_T = a_1 \Delta t_+ a_2 (\Delta_t)^2$,_to_give a good fit_to the results. Taking $m_T$ to_be in Earth masses and $\Delta
_ t$ in minutes, the coefficients ($a_1$,$a_2$)_are (0.044,-1.17$\times$ 10$^{-5}$) for GJ 436b_and (6.787,-0.001) for XO-3b. For_the cases_of HAT-P-2b and HD 17156b,_we find that_even very_low-mass Trojan companions_are dynamically unstable, owing to the_large orbital eccentricities_(see § \[subsec:dynamical\]). Thus, for those systems,_the_requirement of dynamical_stability_is_more constraining_than the empirical_upper_limit on_$m_T$_based on the FG method. (As_will_be described in § 4, this also proved_to be true for_XO-3b_based on the current_data.)

Data Analysis {#sec:data}
=============

The RV data_were taken from the available literature_on each_system. The_references are given in Table 1. These data were generally obtained for_the purpose of discovering or confirming_the planet, although in_a few_cases_the data were_obtained_for other_reasons, such as precisely measuring the orbital_eccentricity (Laughlin_et al. 2005) or for measuring the_Rossiter-McLaughlin effect (Winn et_al. 2006)._Regarding the latter, the data that_were obtained while a transit was_in progress were not used,_to_avoid_the needless complication of incorporating_the Rossiter-McLaughlin effect into the RV_model. However, in_those cases the investigators usually gathered additional_data