 c_s}{8 \pi \rho_0 G M_{\rm c} \chi},
\label{eq:v_MVR17}$$ where $\gamma$ is the adiabatic index of the gas, $L_{\rm c}$ and $M_{\rm c}$ are the luminosity and mass of the perturbing core and $\chi$ is the thermal conductivity of the gas that they assumed to be constant. They noted that the effects of the heating force are expected to be significant for planet formation application, via, e.g., the Earth developing a non-negligible eccentricity and inclination with respect to the protoplanetary disc for realistic disc parameters. These conclusions and the analytic result (Equation \[eq:v\_MVR17\]) were confirmed with numerical simulations by [@VelascoRomeroMasset19], [@ChrenkoLambrechts19] and [@GuileraEtal19].

In application to our particular problem, we note that the density in the central part of the protoplanet is initially homogeneous with a nearly constant $\rho_0 \approx 10^{-9}$ gcm$^{-3}$. Hence, the constant background density results of [@MassetVelascoRomero17] are applicable. However, our protoplanet is finite in extent. When the core is displaced by a distance $R$ from the centre, there is a returning force of protoplanet gravity, given by $$F_{\rm g} = - \frac{G M(R) M_{\rm c}}{R^2}\;,
    \label{eq:Fg0}$$ where $M(R) \approx (4\pi/3) \rho_0 R^3$ is the enclosed mass within radius $R$.

The motion and radius of the orbiting core is then set by a balance of forces. In the azimuthal direction the heating and dynamical friction force specify the core velocity from Equation \[eq:v\_MVR17\]. Given this velocity ($V_0$), the radial position of the core ($R$) is simply set by the gravitational force of the enclosed protoplanet which provides the centrifugal force. We see that $V_0$ is independent of orbit radius $R$ whereas the gravitational force is propotional to $R$, as long as the enclosed gas density is roughly constant. The balance of these therefore establishes an equilibrium radius of the orbit for the core of mass $M_{\rm c}$ and luminosity $L
c_s}{8 \pi \rho_0 G M_{\rm c } \chi }, 
 \label{eq: v_MVR17}$$ where $ \gamma$ is the adiabatic index of the gas, $ L_{\rm c}$ and $ M_{\rm c}$ are the luminosity and mass of the perturbing effect and $ \chi$ is the thermal conduction of the gas that they assumed to be changeless. They notice that the effects of the heating effect are expect to be significant for satellite constitution application, via, for example, the Earth developing a non - negligible eccentricity and inclination with respect to the protoplanetary disc for naturalistic disc parameters. These stopping point and the analytic result (Equation \[eq: v\_MVR17\ ]) were confirm with numerical simulations by [ @VelascoRomeroMasset19 ], [ @ChrenkoLambrechts19 ] and [ @GuileraEtal19 ]. 

 In application to our finical problem, we note that the density in the central region of the protoplanet is initially homogeneous with a nearly constant $ \rho_0 \approx 10^{-9}$   gcm$^{-3}$. Hence, the constant background concentration results of [ @MassetVelascoRomero17 ] are applicable. However, our protoplanet is finite in extent. When the core is displaced by a distance $ R$ from the centre, there is a returning effect of protoplanet graveness, given by $ $ F_{\rm g } = - \frac{G M(R) M_{\rm c}}{R^2}\; , 
     \label{eq: Fg0}$$ where $ M(R) \approx (4\pi/3) \rho_0 R^3 $ is the enclosed mass within radius $ R$. 

 The motion and radius of the orbiting effect is then sic by a balance of forces. In the azimuthal steering the heating system and dynamical friction force specify the core velocity from Equation \[eq: v\_MVR17\ ]. Given this velocity ($ V_0 $), the radial stead of the core ($ R$) is simply set by the gravitational force of the enclosed protoplanet which provides the centrifugal power. We see that $ V_0 $ is independent of orbit radius $ R$ whereas the gravitational force is propotional to $ R$, as long as the envelop gas density is roughly constant. The balance of these consequently establishes an equilibrium radius of the orbit for the core of mass $ M_{\rm c}$ and luminosity $ L
 c_s}{8 \pi \rho_0 G M_{\rm c} \chi},
\label{tq:v_MVR17}$$ where $\gammc$ is thx adiabztic inddx of the gas, $L_{\rm c}$ and $M_{\rm r}$ arw the luminosity and mass ow the pernurbing cire end $\chi$ is the tisrmal conductjyity mh the gas that jhey assumed to be constand. Ghzy noted that the effects of the heaeing fotcf are expected to nq sifnificant for planet formation appmicatioi, via, e.g., the Eatth developing a non-negliglble eccentricity and lnclination witr respect to ghe protoplanetary disd for realistic disc parameters. Thesz conclusiobs anf the analytmc resllt (Equation \[ca:v\_MVR17\]) fere comfirmed with nmmeriral wimulations by [@VelascmRomeroMasset19], [@ChregkoLambrewhcs19] and [@GuileraEtal19].

In qpplicajion do oje pxrtjcnlad probpem, we note tgat the denwity in the central pwgy of the profoplanqt is initially homogeneous with a nearly cohstant $\rho_0 \approx 10^{-9}$ gcm$^{-3}$. Hwnce, the constant bacnground dqnsity results of [@MassetVelascoRomero17] are applicatle. Hkdevtr, our oeohoplanet is finite in extent. When the core is sixpkaced by a disbance $R$ from the crnhrr, there is a rgturniny fkrce of protoplaneh gravijy, givwn by $$F_{\rm g} = - \frac{G M(R) M_{\rm c}}{R^2}\;,
    \label{eq:Dg0}$$ where $M(R) \cppeox (4\pi/3) \rho_0 R^3$ is thz enclosed mcss wijhin rsdius $R$.

The motion and rcdius kf the orbihing core ks then set by a bakatce of forces. In the azimuehal dirertion the hextinb and qynamical vrictlmn force specify tje cote velmcity from Equation \[eq:v\_MVR17\]. Given this velocity ($V_0$), the racidl iosition jf thc core ($R$) is sim[ly set by the gravitctionau force of the enrlosed proto[lanet which [tovides the cxntrifugaj foece. Qe see gfat $V_0$ is indeprndent of orbit raduus $R$ whereas the nravijafional force is pripotional to $R$, ss uond ws tra enclosed gds ddnskyy is roughly cokstxnt. Yhe balance of these thedefore establishes ak equilibtium radits of the orbot for the core of mass $M_{\ro c}$ amd juminosity $L
c_s}{8 \pi \rho_0 G M_{\rm c} \chi}, $\gamma$ the adiabatic of the gas, are luminosity and mass the perturbing core $\chi$ is the thermal conductivity of gas that they assumed to be constant. They noted that the effects of heating force are expected to be significant for planet formation application, via, e.g., Earth a eccentricity inclination with respect to the protoplanetary disc for realistic disc parameters. These conclusions and the analytic (Equation \[eq:v\_MVR17\]) were confirmed with numerical simulations by [@ChrenkoLambrechts19] and [@GuileraEtal19]. In to our particular problem, we that density in central of protoplanet is initially with a nearly constant $\rho_0 \approx 10^{-9}$ gcm$^{-3}$. Hence, the constant background density results of [@MassetVelascoRomero17] are However, our finite in When core displaced by a from the centre, there is a protoplanet gravity, given by $$F_{\rm g} = - M(R) M_{\rm \label{eq:Fg0}$$ where $M(R) \approx (4\pi/3) \rho_0 is the enclosed mass within radius $R$. The and radius of the orbiting core is then set by a balance of forces. In direction the heating and friction force specify core from \[eq:v\_MVR17\]. this velocity the radial position of the core ($R$) is simply set by gravitational force of the enclosed protoplanet which provides the centrifugal see $V_0$ is independent orbit radius $R$ whereas gravitational is propotional to $R$, as enclosed roughly The of these therefore establishes equilibrium radius of the orbit the core of mass
 c_s}{8 \pi \rho_0 G M_{\rm c} \chi},
\label{eq:v_MVr17}$$ where $\gammA$ is thE adIabAtIc inDex oF the gas, $L_{\rm c}$ and $m_{\Rm c}$ aRe the luminosity and mass Of the PeRTurbINg Core aNd $\chi$ is THe THErmAl CoNduCtIViTy of tHe gAs that tHey assumed To bE cOnstant. They nOTeD that the efFecTs of the heatiNg fOrce arE eXpeCTed to Be sIgnifIcant fOR planeT formatioN aPPlicatIOn, via, e.g., THE EArth Developing a non-negLIgIBle eccentricitY and inClINaTIOn wIth Respect to tHe ProtoPLanetarY DiSC FOr rEAlistic disc paRameters. TheSE coNclusiOnS anD The anaLytic ReSUlt (equation \[eq:v\_mVR17\]) wEre confirMed witH NumericAL simulaTions bY [@VeLasCoRoMErOMAssEt19], [@cHreNKolamBRecHts19] and [@GuIlErAEtal19].

in apPLICAtioN to Our pArticUlar problem, we NotE thaT The DensiTy in tHe ceNtRal paRt of thE protOpLanet is initiallY homOgeneous wIth A nEarLy ConstANt $\rho_0 \aPprOx 10^{-9}$ gCm$^{-3}$. Hence, The consTAnt BaCKGRoUnd density results oF [@MASSeTVelascoromero17] ARe ApPLicable. HOwEveR, our PROtoplAnet IS fInite in eXtent. WHEn ThE core is DiSplaceD bY a dIstAnce $R$ FRom tHe centRe, there iS a retURning force of prOToplanet graviTY, gIVEn BY $$F_{\rm G} = - \frAc{G M(R) M_{\rm c}}{R^2}\;,
    \lAbel{EQ:Fg0}$$ wHere $m(r) \aPprOX (4\pi/3) \rhO_0 R^3$ is tHe ENcLOsed mass within radiuS $R$.

the motIon anD radius of the oRbiting corE IS Then set bY a baLAnCE of forces. In the AzimuThal directIOn the heaTing aNd dynamiCal frictiON Force speCifY thE coRe vELOcIty from EquatiON \[Eq:v\_MvR17\]. given thIs vElocity ($v_0$), thE raDiaL poSiTion of the Core ($R$) is sImPlY sEt By tHe graVItationaL fOrcE oF thE enclOSed proToplaNet wHiCh PRovIdes the CEnTRIfugAl FoRce. WE seE tHat $V_0$ iS indEPenDent of oRbit radiuS $R$ wHEreaS tHe GravitaTional force is PrOpotional tO $R$, As lOng as tHE Enclosed Gas density is roughly consTAnt. The bAlaNce of ThesE thereforE esTablisHes AN equilIbrium RadiuS oF thE ORbit fOR ThE coRe Of mass $M_{\rm c}$ AND luMinosItY $L
 c_s}{8 \pi \rho_0 G M_{\r m c} \chi} ,
\la bel {eq :v _MVR 17}$ $ where $\gamm a $ is  the adiabatic index o f the  g a s, $ L _{ \rm c }$ and$ M_ { \ rmc} $are  t h elumin osi ty andmass of th e p er turbing core an d $\chi$ i s t he thermal c ond uctivi ty  of the g asthatthey a s sumedto be con st a nt. Th e y noted t ha t th e effects of theh ea t ing force areexpect ed to b e s ign ificant fo rplane t  format i on a p pli c ation, via, e .g., the Ea r thdevelo pi nga  non-n eglig ib l e e ccentricity  and  inclinat ion wi t h respe c t to th e prot opl ane tary di sc  fo rr eal i st icd isc  paramet er s.  Thes e co n c l u sion s a nd t he an alytic result  (E quat i on\[eq: v\_MV R17\ ])  were  confi rmedwi th numerical si mula tions by[@V el asc oR omero M asset1 9],  [@ Chrenko Lambrec h ts1 9] a n d[@GuileraEtal19].
I n ap plicatio n to o u rpa r ticularpr obl em,w e  note  tha t  t he densi ty int he  c entralpa rt ofth e p rot oplan e t is  initi ally hom ogene o us with a near l y constant $\ r ho _ 0  \ a ppro x 1 0^{-9}$ gcm $^{- 3 }$.Henc e ,the const ant b ac k gr o und density results  o f [@Ma ssetV elascoRomero1 7] are app l i c able. Ho weve r ,o ur protoplanet  is f inite in e x tent. Wh en th e core i s displac e d  by a di sta nce  $R $ f r o mthe centre, t h e re i sa retur nin g force  of  pr oto pla ne t gravity , givenby  $ $F _{ \rm  g} = - \frac{ GM(R )M_{ \rm c } }{R^2} \;,
   \ la be l {eq :Fg0}$$ wh e r e $M (R )\app rox  ( 4\pi/ 3) \ r ho_ 0 R^3$is the en clo s ed m as swithinradius $R$.
Th e motion a nd  ra dius o f the orbi ting core is then set b y  a bala nce  of f orce s. In the  az imutha l d i rectio n theheati ng  an d dynam i c al  fr ic tion force s pec ify t he  cor e veloc ity from Equation\ [eq :v\_MVR17\].Giv en t h i svel o ci t y ( $V _ 0$) , the radial posi tion of th ec or e ($R$) is sim pl y set b y the g ravit a tionalforce ofthe enclo se d pr o t opl anet which  provide s the cen t rifug a lforce . W e seeth at$V_0$  is in d epe ndent  of or bi t radi us $R $whereasthe gravitational force  is pr opoti ona l to $R$,  as lon g as theencl osed gas d ens ity  is r oug h ly co nsta n t.  Th e  bala nceo f these t h er efo r e  e stablishesa n equ ilibr ium radius  ofthe orbit for the core of mass $ M_{\ r m  c} $ a n d lu mi nosity $L
 c_s}{8_\pi \rho_0_G M_{\rm c} \chi},
\label{eq:v_MVR17}$$_where $\gamma$_is_the adiabatic_index_of the gas,_$L_{\rm c}$ and_$M_{\rm c}$ are the_luminosity and mass_of_the perturbing core and $\chi$ is the thermal conductivity of the gas that they_assumed_to be_constant._They_noted that the effects of_the heating force are expected_to be_significant for planet formation application, via, e.g., the_Earth_developing a non-negligible_eccentricity and inclination with respect to the protoplanetary disc_for realistic disc parameters. These conclusions_and the analytic_result_(Equation_\[eq:v\_MVR17\]) were confirmed with_numerical simulations by [@VelascoRomeroMasset19], [@ChrenkoLambrechts19] and_[@GuileraEtal19].

In application to our particular problem,_we note that the density in the_central part of the protoplanet is_initially homogeneous with a nearly_constant $\rho_0_\approx 10^{-9}$ gcm$^{-3}$. Hence, the constant_background density results_of [@MassetVelascoRomero17]_are applicable. However,_our protoplanet is finite in extent._When the core_is displaced by a distance $R$_from_the centre, there_is_a_returning force_of protoplanet gravity,_given_by $$F_{\rm_g}_= - \frac{G M(R) M_{\rm c}}{R^2}\;,
__  \label{eq:Fg0}$$ where $M(R) \approx (4\pi/3)_\rho_0 R^3$ is the_enclosed_mass within radius $R$.

The_motion and radius of the_orbiting core is then set by_a balance_of forces._In the azimuthal direction the heating and dynamical friction force specify_the core velocity from Equation \[eq:v\_MVR17\]._Given this velocity ($V_0$),_the radial_position_of the core_($R$)_is simply_set by the gravitational force of the_enclosed protoplanet_which provides the centrifugal force. We_see that $V_0$ is_independent_of orbit radius $R$ whereas the_gravitational force is propotional to $R$,_as long as the enclosed_gas_density_is roughly constant. The balance_of these therefore establishes an equilibrium_radius of the_orbit for the core of mass $M_{\rm_c}$_and luminosity $L