 the Subaru telescope, we also investigate the case of $m_{\rm d} = 73 M_{\oplus}$ here.

### Calculation of the simulated TTV

Orbits of the planets are integrated using the fourth-order Hermite scheme with the shared time step [@2004PASJ...56..861K]. The time step (typically $\sim 0.025$ days) is chosen so that the fractional energy change due to the integrator during $\sim 1000$ days integration should always be smaller than $10^{-9}$. All the simulations presented in this section integrate planetary orbits beginning at the same epoch $T_0 (\mathrm{BJD}) = 2455189.1703$ (the first transit time of KOI-94d), until $\mathrm{BJD} = 2456133.0$ (approximately the last transit time of KOI-94d we analyzed).

For each planet, the initial value of the orbital inclination $i$ is fixed at the value obtained from $a/R_{\star}$ and $b$ (see Table \[ppar\_fit\]), assuming that $e = 0$. Here, all the $i$ values are chosen in the range $[0, \pi/2]$.[^2] The values of $\Omega_{\rm d}$ and $\Omega_{\rm e}$ are fixed at those in Table \[ppar\_prior\]. Since we have no information on the nodal angle of KOI-94c, we assume $\Omega_{\rm c} = 0$ deg. Initial semi-major axes are calculated via Kepler’s third law with $M_{\star} = 1.25 M_{\odot}$ [@2012ApJ...759L..36H], orbital periods in Table \[ep\], and planetary masses adopted in each simulation. The phases of the planets are determined from the transit ephemerides: for each planet, we convert the transit time closest to $T_0$ into the sum of the argument of periastron and mean anomaly $\omega + f$, taking account of the non-zero eccentricity if any, and then move it backward in time to $T_0$, assuming a Keplerian orbit.

The mid-transit times of each planet are determined by minimizing the sky-plane distance $D$ between the star and the planet, where the roots of the time derivative of $D$ are found by the Newton-Raphson method [@2010arXiv1006.3834F]. Then these transit times are
the Subaru telescope, we also investigate the case of $ m_{\rm d } = 73 M_{\oplus}$ here. 

 # # # Calculation of the simulated TTV 

 Orbits of the planets are integrate use the fourth - order Hermite scheme with the shared fourth dimension step [ @2004PASJ... 56.. 861 K ]. The time footfall (typically $ \sim 0.025 $ day) is chosen so that the fractional energy change due to the integrator during $ \sim 1000 $ days consolidation should always be smaller than $ 10^{-9}$. All the simulations presented in this incision integrate planetary orbits beginning at the same epoch $ T_0 (\mathrm{BJD }) = 2455189.1703 $ (the inaugural transit time of KOI-94d), until $ \mathrm{BJD } = 2456133.0 $ (approximately the last passage time of KOI-94d we analyzed). 

 For each planet, the initial value of the orbital dip $ i$ is fixed at the value obtained from $ a / R_{\star}$ and $ b$ (see Table \[ppar\_fit\ ]), assuming that $ east = 0$. Here, all the $ i$ values are chosen in the range $ [ 0, \pi/2]$.[^2 ] The values of $ \Omega_{\rm d}$ and $ \Omega_{\rm e}$ are fixed at those in Table \[ppar\_prior\ ]. Since we have no information on the nodal angle of KOI-94c, we assume $ \Omega_{\rm c } = 0 $ deg. Initial semi - major axes are calculated via Kepler ’s third law with $ M_{\star } = 1.25 M_{\odot}$ [ @2012ApJ... 759L.. 36H ], orbital periods in Table \[ep\ ], and planetary masses adopted in each simulation. The phase of the planets are determined from the passage ephemerides: for each planet, we convert the transit time close to $ T_0 $ into the sum of the argument of periastron and mean anomaly $ \omega + f$, taking bill of the non - zero eccentricity if any, and then move it backward in time to $ T_0 $, assuming a Keplerian orbit. 

 The mid - transit time of each planet are determined by minimizing the sky - plane distance $ D$ between the star and the planet, where the solution of the time derivative of $ D$ are found by the Newton - Raphson method [ @2010arXiv1006.3834F ]. Then these transit time are
 thf Subaru telescope, we alro investigate jhw case of $m_{\rj d} = 73 M_{\oolus}$ here.

### Calculation of the sumulauvd TTV

Orbits of the puanets arv integrared nsing the fourth-order Heviite dchekx with the shargd time step [@2004PASJ...56..861K]. The tima rtzp (typically $\sim 0.025$ days) is chosen so ehat thr vractional enetgy cnwnge due to the integrator during $\sim 1000$ days iitegration shoukd always be smaller than $10^{-9}$. All the simulations pgesented in thif section intdgrate planetary orbita beginning at the same epoch $T_0 (\mathxm{BJD}) = 2455189.1703$ (the fursh transit tine of KOI-94d), until $\mathrm{BJG} = 2456133.0$ (apptoximately the lavt rransit time of KOI-94d xe analyzed).

For each klanet, the iuitial value of the oebutal hnclhnatkin $k$ ia hixsd at hhe value obtzined from $q/R_{\star}$ and $b$ (see Tanlq \[ppar\_fit\]), assujing trae $e = 0$. Here, all the $i$ values are chosen pn tge range $[0, \pi/2]$.[^2] The values of $\Omega_{\rm d}$ and $\Omeha_{\rm e}$ arq fixed at those in Table \[ppar\_prior\]. Since we have no iifurmcblon ub hhe nodal angle of KOI-94c, we assume $\Omega_{\rm c} = 0$ seb. Pnitial semi-major axes are cslfukwted via Keplgr’s thixs maw with $M_{\star} = 1.25 M_{\ldot}$ [@2012ApT...759L..36H], oebital pewiodx in Table \[ep\], and planetary masses adopned un each simulation. The phases of jhe plsnets are determined frum tge transit fphemeridsr: for each planeg, wv cotvert tht transit time clofest to $T_0$ intp the sjm og the wrgument ov perldstron and mean anlmaly $\okega + f$, taning account of the non-zero eccxitricity if amy, ang then mjve ib backward in tyme to $T_0$, assumnng a Ke'leriav orbit.

The mid-traisit times os each planet wre determinxd by minymizung rhe sky-ouane distance $C$ between the star qnd the planet, wheve thg doots of the tine eerivative of $D$ ard fjujd br the Newton-Rdphsun orthod [@2010arXiv1006.3834F]. Thek tfese transit times are
the Subaru telescope, we also investigate the $m_{\rm = 73 here. ### Calculation of planets are integrated the fourth-order Hermite with the shared time step [@2004PASJ...56..861K]. time step (typically $\sim 0.025$ days) is chosen so that the fractional energy due to the integrator during $\sim 1000$ days integration should always be smaller $10^{-9}$. the presented this section integrate planetary orbits beginning at the same epoch $T_0 (\mathrm{BJD}) = 2455189.1703$ (the first time of KOI-94d), until $\mathrm{BJD} = 2456133.0$ (approximately last transit time of we analyzed). For each planet, initial of the inclination is at the value from $a/R_{\star}$ and $b$ (see Table \[ppar\_fit\]), assuming that $e = 0$. Here, all the $i$ values chosen in $[0, \pi/2]$.[^2] values $\Omega_{\rm and $\Omega_{\rm e}$ at those in Table \[ppar\_prior\]. Since information on the nodal angle of KOI-94c, we $\Omega_{\rm c} 0$ deg. Initial semi-major axes are via Kepler’s third law with $M_{\star} = 1.25 [@2012ApJ...759L..36H], orbital periods in Table \[ep\], and planetary masses adopted in each simulation. The phases planets are determined from transit ephemerides: for planet, convert transit closest to into the sum of the argument of periastron and mean anomaly + f$, taking account of the non-zero eccentricity if any, move backward in time $T_0$, assuming a Keplerian The times of each planet by the between star the planet, where the of the time derivative of are found by the transit times are
 the Subaru telescope, we also iNvestigate The caSe oF $m_{\rM d} = 73 m_{\oplUs}$ heRe.

### Calculation oF The sImulated TTV

Orbits of the PlaneTs ARe inTEgRated Using thE FoURTh-oRdEr herMiTE sCheme WitH the shaRed time steP [@2004PAsJ...56..861k]. The time step (TYpIcally $\sim 0.025$ dAys) Is chosen so thAt tHe fracTiOnaL EnergY chAnge dUe to thE IntegrAtor durinG $\sIM 1000$ days iNTegratiON ShOuld Always be smaller thAN $10^{-9}$. ALL the simulationS preseNtED iN THis SecTion integrAtE planETary orbITs BEGInnINg at the same epOch $T_0 (\mathrm{Bjd}) = 2455189.1703$ (thE first TrAnsIT time oF KOI-94d), UnTIl $\mAthrm{BJD} = 2456133.0$ (appRoxiMately the Last trANsit timE Of KOI-94d wE analyZed).

for Each PLaNeT, thE iNItiAL vAluE Of tHe orbitaL iNcLinatIon $i$ IS FIXed aT thE valUe obtAined from $a/R_{\stAr}$ aNd $b$ (sEE TaBle \[ppAr\_fit\]), AssuMiNg thaT $e = 0$. Here, All thE $i$ Values are chosen In thE range $[0, \pi/2]$.[^2] THe vAlUes Of $\omega_{\RM d}$ and $\OMegA_{\rm E}$ are fixEd at thoSE in taBLE \[PpAr\_prior\]. Since we have No INFoRmation oN the noDAl AnGLe of KOI-94c, We AssUme $\OMEGa_{\rm c} = 0$ Deg. INItIal semi-mAjor axES aRe CalculaTeD via KePlEr’s ThiRd law WIth $M_{\Star} = 1.25 M_{\oDot}$ [@2012ApJ...759L..36H], OrbitAL periods in TablE \[Ep\], and planetarY MaSSEs ADoptEd iN each simulaTion. tHe phAses OF tHe pLAnets Are deTeRMiNEd from the transit ephEmErides: For eaCh planet, we conVert the traNSIT time cloSest TO $T_0$ INto the sum of the ArgumEnt of periaSTron and mEan anOmaly $\omeGa + f$, taking ACCount of tHe nOn-zEro EccENTrIcity if any, and THEn moVe It backwArd In time tO $T_0$, aSsuMinG a KEpLerian orbIt.

The mid-TrAnSiT tImeS of eaCH planet aRe DetErMinEd by mINimiziNg the Sky-pLaNe DIstAnce $D$ beTWeEN The sTaR aNd thE plAnEt, wheRe thE RooTs of the Time derivAtiVE of $D$ ArE fOund by tHe Newton-RaphsOn Method [@2010arXiV1006.3834F]. theN these TRAnsit timEs are
 the Subaru telescope, wealso inves tigat e t heca se o f $m _{\rm d} = 73M _{\o plus}$ here.

### Calc ulati on of t h esimul ated TT V 
O r bit sof  th ep la netsare  integr ated using  th efourth-order He rmite sche mewith the sha red  timest ep[ @2004 PAS J...5 6..861 K ]. The  time ste p( typica l ly $\si m 0. 025$  days) is chosens ot hat the fracti onal e ne r gy c han gedue to the  i ntegr a tor dur i ng $ \ sim 1000$ days in tegration s h oul d alwa ys  be smalle r tha n$ 10^ {-9}$. Allthesimulatio ns pre s ented i n  this s ection  in teg rate pl an eta ry orb i ts  be g inn ing at t he  s ame e poch $ T _ 0 (\ mat hrm{ BJD})  = 2455189.17 03$  (th e  fi rst t ransi t ti me  of K OI-94d ), un ti l $\mathrm{BJD}  = 2 456133.0$  (a pp rox im ately the la sttra nsit ti me of K O I-9 4d w e  a nalyzed).

For eac hp l an et, theinitia l  v al u e of the  o rbi tali n clina tion $i $ is fix ed att he  v alue ob ta ined f ro m $ a/R _{\st a r}$and $b $ (see T able\ [ppar\_fit\]), assuming that $e =  0 $ . He re,  all the $i $ va l uesarec ho sen in th e ran ge $[ 0 , \pi/2]$.[^2] Theva lues o f $\O mega_{\rm d}$  and $\Ome g a _ {\rm e}$  are fi x ed at those in  Tabl e \[ppar\_ p rior\].Since  we have  no infor m a tion onthe  no dal  an g l eof KOI-94c, w e assu me  $\Omeg a_{ \rm c}= 0 $ d eg.  In it ial semi- major ax es  a re  c alc ulate d  via Kep le r’s  t hir d law with $ M_{\s tar}  =  1 . 25M_{\odo t }$ [ @201 2A pJ ...7 59L .. 36H],  orb i tal  period s in Tabl e \ [ ep\] ,an d plane tary masses a do pted in ea ch  si mulati o n . The ph ases of the planets are determi ned  from  the  transiteph emerid es: for ea ch pla net,we  co n v ert t h e  t ran si t time clo s e stto $T _0 $ in to thesum of the argumen t  of  periastron a ndmean a no mal y  $ \ ome ga + f $ ,  taking account  of the no n- z er o eccentri c ity  i f any,and the n mov e  it bac kward intime to $ T_ 0$,a s sum ing a Kepl erian or bit.

The mid-t r an sit t ime s of e ac h p lanet  are d e ter mined  by mi ni mizing  thesk y-planedistance $D$ between th e star  andthe  planet,whe r e t he rootsof t he time de riv ati ve of  $D $  arefoun d  b y t h e New ton- R aphson me t ho d [ @ 2 01 0arXiv1006. 3 8 3 4F] . The n t h ese tr ansi t times are
 the_Subaru telescope,_we also investigate the_case of_$m_{\rm_d} =_73_M_{\oplus}$ here.

### Calculation_of the simulated_TTV

Orbits of the planets_are integrated using_the_fourth-order Hermite scheme with the shared time step [@2004PASJ...56..861K]. The time step (typically $\sim_0.025$_days) is_chosen_so_that the fractional energy change_due to the integrator during_$\sim 1000$_days integration should always be smaller than $10^{-9}$._All_the simulations presented_in this section integrate planetary orbits beginning at the_same epoch $T_0 (\mathrm{BJD}) = 2455189.1703$_(the first transit_time_of_KOI-94d), until $\mathrm{BJD} =_2456133.0$ (approximately the last transit time_of KOI-94d we analyzed).

For each planet,_the initial value of the orbital inclination_$i$ is fixed at the value_obtained from $a/R_{\star}$ and $b$_(see Table_\[ppar\_fit\]), assuming that $e =_0$. Here, all_the $i$_values are chosen_in the range $[0, \pi/2]$.[^2] The_values of $\Omega_{\rm_d}$ and $\Omega_{\rm e}$ are fixed_at_those in Table_\[ppar\_prior\]._Since_we have_no information on_the_nodal angle_of_KOI-94c, we assume $\Omega_{\rm c} =_0$_deg. Initial semi-major axes are calculated via_Kepler’s third law with_$M_{\star}_= 1.25 M_{\odot}$ [@2012ApJ...759L..36H],_orbital periods in Table \[ep\],_and planetary masses adopted in each_simulation. The_phases of_the planets are determined from the transit ephemerides: for each planet,_we convert the transit time closest_to $T_0$ into the_sum of_the_argument of periastron_and_mean anomaly_$\omega + f$, taking account of the_non-zero eccentricity_if any, and then move it_backward in time to_$T_0$,_assuming a Keplerian orbit.

The mid-transit times_of each planet are determined by_minimizing the sky-plane distance $D$_between_the_star and the planet, where_the roots of the time derivative_of $D$ are_found by the Newton-Raphson method [@2010arXiv1006.3834F]. Then_these_transit times are