gtrsim \rm{O}(100) \rm{AU}$) [e.g. @Boley_CA_and_GI] though a hybrid scenario of forming gas giants in the same system by both core accretion and gravitational instability has yet to be modelled.

There are two quantities that have historically been used to determine whether a disc is likely to fragment. The first is the stability parameter [@Toomre_stability1964],

$$\label{eq:Toomre}
  Q=\frac{c_{\rm s}\kappa_{\rm ep}}{\pi\Sigma G},$$

where $c_{\rm s}$ is the sound speed in the disc, $\kappa_{\rm ep}$ is the epicyclic frequency, which for Keplerian discs is approximately equal to the angular frequency, $\Omega$, $\Sigma$ is the surface mass density and $G$ is the gravitational constant. Therefore, once the surface mass density and the rotation of the disc have been established, the stability is purely dependent on the disc temperature. [@Toomre_stability1964] showed that for an infinitesimally thin disc to fragment, the stability parameter must be less than a critical value, $Q_{\rm crit} \approx 1$.

[@Gammie_betacool] showed that in addition to the stability criterion above, the disc must cool at a fast enough rate. Using shearing sheet simulations, he showed that if the cooling timescale can be parametrized as

$$\label{eq:beta}
  \beta = t_{\rm cool}\Omega,$$

where

$$\label{eq:tcool}
  t_{\rm cool} = u \big(\frac{du_{\rm cool}}{dt}\big)^{-1},$$

$u$ is the internal energy and $du_{\rm cool}/dt$ is the total cooling rate, then for fragmentation we require $\beta \lesssim 3$, for a ratio of specific heats $\gamma = 2$ (in two dimensions). [-@Rice_beta_condition] carried out three-dimensional simulations using a smoothed particle hydrodynamics (SPH) code and showed that this cooling parameter is dependent on the equation of state so that fragmentation can occur for discs with $\gamma = 5/3$ and $7/5$ if $\beta \lesssim 7$ and $13$, respectively. [-@Libby_MSci] showed that the critical value of $\beta$ (below which fragmentation will occur if the stability
gtrsim \rm{O}(100) \rm{AU}$) [ e.g. @Boley_CA_and_GI ] though a hybrid scenario of forming gas giants in the like arrangement by both core accretion and gravitational instability experience even to be modelled. 

 There are two measure that have historically been used to determine whether a phonograph record is likely to fragment. The beginning is the stability parameter [ @Toomre_stability1964 ], 

 $ $ \label{eq: Toomre } 
   Q=\frac{c_{\rm s}\kappa_{\rm ep}}{\pi\Sigma G},$$ 

 where $ c_{\rm s}$ is the reasoned speed in the disc, $ \kappa_{\rm ep}$ is the epicyclic frequency, which for Keplerian discs is approximately equal to the angular frequency, $ \Omega$, $ \Sigma$ is the surface mass concentration and $ G$ is the gravitational constant. Therefore, once the surface mass concentration and the rotation of the disc have been established, the stability is purely dependent on the phonograph record temperature. [ @Toomre_stability1964 ] showed that for an infinitesimally thin disc to fragment, the stability argument must be less than a critical value, $ Q_{\rm crit } \approx 1$. 

 [ @Gammie_betacool ] showed that in addition to the stability criterion above, the disc must cool at a fast adequate pace. practice shearing sheet simulations, he showed that if the cooling timescale can be parametrized as 

 $ $ \label{eq: beta } 
   \beta = t_{\rm cool}\Omega,$$ 

 where 

 $ $ \label{eq: tcool } 
   t_{\rm cool } = u \big(\frac{du_{\rm cool}}{dt}\big)^{-1},$$ 

 $ u$ is the internal energy and $ du_{\rm cool}/dt$ is the entire cooling system pace, then for fragmentation we require $ \beta \lesssim 3 $, for a ratio of specific heat $ \gamma = 2 $ (in two dimensions). [ -@Rice_beta_condition ] carried out three - dimensional simulations use a smoothed atom hydrodynamics (SPH) code and showed that this cooling parameter is pendent on the equation of state so that fragmentation can occur for discs with $ \gamma = 5/3 $ and $ 7/5 $ if $ \beta \lesssim 7 $ and $ 13 $, respectively. [ -@Libby_MSci ] showed that the critical value of $ \beta$ (below which fragmentation will occur if the stability
gtrdim \rm{O}(100) \rm{AU}$) [e.g. @Boley_CA_akd_GI] though a hyyeid scxnario kf formivg gas giants in the same sydtwm by both core accretion avd gravitwtional unstebility has yet vk be modelled.

Fmere crx two quantitiex that hava historically bdeu used to determine whether a disc if likelu ho fragment. Thg firxe is nht stability parameter [@Toomre_stabimity1964],

$$\labtl{eq:Toomre}
  Q=\frac{c_{\tm s}\kappa_{\rm ep}}{\pi\Sigma G},$$

whege $c_{\gm s}$ is the sound dpeed in thg difx, $\kappa_{\rm ep}$ is the eppeyclic freqoency, which for Keplerian discs ks ap'roximately ewuap to the angnlar fgequency, $\Omega$, $\Sigma$ hs the xurface mass dcnsitb ane $G$ is the gravitatioial constant. Therefote, once tha aurface mass densutt and the rotxriov or vhe disc javx been estaglished, the stability is purelu qvlendent on tge diss eemperature. [@Toomre_stability1964] showed that fod an infinitesimally thun disc to fragment, tje stabilyty parameter must be less than a critical value, $X_{\rm ccig} \akpvjb 1$.

[@Gwmmie_betacool] showed that in addition to the ffanikity criterion above, the cidc iust cool at x fast enkugh rate. Using shfaring fheet simulatijns, ne showed that if the coolibg timescale xan be parametrizeb as

$$\label{eq:bzta}
  \bgta = t_{\tm cool}\Omega,$$

where

$$\label{ed:tcokl}
  t_{\rm coop} = u \big(\fdxc{du_{\rm cool}}{dt}\big)^{-1},$$

$j$ ix dhe internal energy and $du_{\wm cool}/dt$ is che totau copling wate, then vor fvdgmentation we reqkire $\yeta \nesssim 3$, flr a ratio of specific heats $\gamma = 2$ (in two domanspons). [-@Rice_yeta_cokdition] carried out three-dimeusional fimulxtions usihg a smmothed partycle hydrodyndlics (SPH) codx and shored rhat this cuuling parametet is depeubent on tye equation of stabe so fhat fragmentatnub can occur for dircs wpth $\gwkma = 5/3$ and $7/5$ ix $\bega \ursssio 7$ and $13$, resiecgiveky. [-@Libby_MSci] showed dhat the critical valur jf $\beta$ (velow whych fragmentayion will occur if the vtauility
gtrsim \rm{O}(100) \rm{AU}$) [e.g. @Boley_CA_and_GI] though a of gas giants the same system gravitational has yet to modelled. There are quantities that have historically been used determine whether a disc is likely to fragment. The first is the stability [@Toomre_stability1964], $$\label{eq:Toomre} Q=\frac{c_{\rm s}\kappa_{\rm ep}}{\pi\Sigma G},$$ where $c_{\rm s}$ is the sound speed the $\kappa_{\rm is epicyclic frequency, which for Keplerian discs is approximately equal to the angular frequency, $\Omega$, $\Sigma$ is surface mass density and $G$ is the gravitational Therefore, once the surface density and the rotation of disc been established, stability purely on the disc [@Toomre_stability1964] showed that for an infinitesimally thin disc to fragment, the stability parameter must be less than critical value, \approx 1$. showed in to the stability the disc must cool at a Using shearing sheet simulations, he showed that if cooling timescale be parametrized as $$\label{eq:beta} \beta = cool}\Omega,$$ where $$\label{eq:tcool} t_{\rm cool} = u \big(\frac{du_{\rm $u$ is the internal energy and $du_{\rm cool}/dt$ is the total cooling rate, then for require $\beta \lesssim 3$, a ratio of heats = (in dimensions). [-@Rice_beta_condition] out three-dimensional simulations using a smoothed particle hydrodynamics (SPH) code and that this cooling parameter is dependent on the equation of that can occur for with $\gamma = 5/3$ $7/5$ $\beta \lesssim 7$ and [-@Libby_MSci] that of (below fragmentation will occur if stability
gtrsim \rm{O}(100) \rm{AU}$) [e.g. @Boley_CA_and_gI] though a hYbrid SceNarIo Of foRminG gas giants in thE Same System by both core accretIon anD gRAvitATiOnal iNstabilITy HAS yeT tO bE moDeLLeD.

TherE arE two quaNtities thaT haVe Historically BEeN used to detErmIne whether a dIsc Is likeLy To fRAgmenT. ThE firsT is the STabiliTy parametEr [@tOomre_sTAbility1964],

$$\LABeL{eq:TOomre}
  Q=\frac{c_{\rm s}\kapPA_{\rM Ep}}{\pi\Sigma G},$$

wherE $c_{\rm s}$ iS tHE sOUNd sPeeD in the disc, $\KaPpa_{\rm EP}$ is the ePIcYCLIc fREquency, which fOr Keplerian DIscS is appRoXimATely eqUal to ThE AngUlar frequenCy, $\OmEga$, $\Sigma$ iS the suRFace masS Density And $G$ is The GraVitaTIoNaL coNsTAnt. tHeRefORe, oNce the suRfAcE mass DensITY ANd thE roTatiOn of tHe disc have beeN esTablISheD, the sTabilIty iS pUrely DependEnt on ThE disc temperaturE. [@TooMre_stabilIty1964] ShOweD tHat foR An infiNitEsiMally thIn disc tO FraGmENT, ThE stability parameteR mUST bE less thaN a critICaL vALue, $Q_{\rm crIt} \AppRox 1$.

[@GAMMie_beTacoOL] sHowed thaT in addITiOn To the stAbIlity cRiTerIon Above, THe diSc must Cool at a fAst enOUgh rate. Using shEAring sheet simULaTIOnS, He shOweD that if the cOoliNG timEscaLE cAn bE ParamEtrizEd AS

$$\lABel{eq:beta}
  \beta = t_{\rm cooL}\OMega,$$

whEre

$$\laBel{eq:tcool}
  t_{\rm Cool} = u \big(\frAC{DU_{\rm cool}}{dT}\big)^{-1},$$

$U$ Is THe internal enerGy and $Du_{\rm cool}/dt$ IS the totaL coolIng rate, tHen for fraGMEntation We rEquIre $\BetA \LEsSsim 3$, for a ratio OF SpecIfIc heats $\GamMa = 2$ (in two DimEnsIonS). [-@RiCe_Beta_condiTion] carrIeD oUt ThRee-DimenSIonal simUlAtiOnS usIng a sMOothed PartiCle hYdRoDYnaMics (SPH) COdE ANd shOwEd That ThiS cOolinG parAMetEr is depEndent on tHe eQUatiOn Of State so That fragmentaTiOn can occur FoR diScs witH $\GAmma = 5/3$ and $7/5$ iF $\beta \lesssim 7$ and $13$, respectiVEly. [-@LibbY_MSCi] shoWed tHat the criTicAl valuE of $\BEta$ (belOw whicH fragMeNtaTIOn wilL OCcUr iF tHe stabilitY
gtrsim \rm{O}(100) \rm{AU} $) [e.g. @ Boley _CA _an d_ GI]thou gh a hybrid sc e nari o of forming gas giant s inth e  sam e  s ystem  by bot h  c o r e a cc re tio na nd  grav ita tionalinstabilit y h as  yet to be m o de lled.

The reare two quan tit ies th at  ha v e his tor icall y been used t o determi ne whethe r  a disc i slike ly to fragment. T h ef irst is the st abilit yp ar a m ete r [ @Toomre_st ab ility 1 964],
$ $\ l a b el{ e q:Toomre}
  Q =\frac{c_{\ r m s }\kapp a_ {\r m  ep}}{ \pi\S ig m a G },$$

where  $c_ {\rm s}$is the sound s p eed inthe di sc,  $\ kapp a _{ \r m e p} $  is th e e p icy clic fre qu en cy, w hich f o r  Kep ler iandiscs  is approxima tel y eq u alto th e ang ular  f reque ncy, $ \Omeg a$ , $\Sigma$ is t he s urface ma ssde nsi ty  and$ G$ isthe  gr avitati onal co n sta nt . T he refore, once the s ur f a ce  mass de nsitya nd  t h e rotati on  of  the d isc h aveb ee n establ ished, th estabili ty  is pu re lydep enden t  onthe di sc tempe ratur e . [@Toomre_sta b ility1964] sh o we d th a t fo r a n infinites imal l y th in d i sc  to fragm ent,th e  s t ability parameter m us t be l ess t han a critica l value, $ Q _ { \rm crit } \a p pr o x 1$.

[@Gammi e_bet acool] sho w ed thatin ad dition t o the sta b i lity cri ter ion  ab ove , th e disc must c o o l at  a  fast e nou gh rate . U sin g s hea ri ng sheetsimulati on s,  h esho wed t h at if th ecoo li ngtimes c ale ca n bepara me tr i zed  as

$$ \ la b e l{eq :b et a}
 \b et a = t _{\r m  co ol}\Ome ga,$$

wh ere 
$$\ la be l{eq:tc ool}
  t_{\rm  c ool} = u \ bi g(\ frac{d u _ {\rm coo l}}{dt}\big)^{-1},$$

$ u $ is th e i ntern al e nergy and  $d u_{\rm  co o l}/dt$  is th e tot al  co o l ing r a t e,  th en  for fragm e n tat ion w erequ ire $\b eta \lesssim 3$, f o r a  ratio of spe cif ic h e a ts  $\ g am m a =  2 $  (i n two dimensions) . [-@Rice_ be t a_ condition] car ri ed outthree-d imens i onal si mulations  using asm ooth e d  pa rticle hyd rodynami cs (SPH)c ode a n dshowe d t hat th is  co oling  param e ter  is d epende nt  on th e equ at ion of s tate so that fragmentat ion ca n occ urfor discs  wi t h $ \gamma =5/3$  and $7/5$  if  $\ beta\le s ssim7$ a n d$13 $ , res pect i vely. [-@ L ib by_ M S ci ] showed th a t the  crit ica l  value  of$\beta$ (below wh i ch fragmentati on w i l l o ccu r  ifth e stability
gtrsim \rm{O}(100)_\rm{AU}$) [e.g._@Boley_CA_and_GI] though a hybrid_scenario of_forming_gas giants_in_the same system_by both core_accretion and gravitational instability_has yet to_be_modelled.

There are two quantities that have historically been used to determine whether a disc_is_likely to_fragment._The_first is the stability parameter_[@Toomre_stability1964],

$$\label{eq:Toomre}
  Q=\frac{c_{\rm s}\kappa_{\rm ep}}{\pi\Sigma_G},$$

where $c_{\rm_s}$ is the sound speed in the disc,_$\kappa_{\rm_ep}$ is the_epicyclic frequency, which for Keplerian discs is approximately equal_to the angular frequency, $\Omega$, $\Sigma$_is the surface_mass_density_and $G$ is the_gravitational constant. Therefore, once the surface_mass density and the rotation of_the disc have been established, the stability_is purely dependent on the disc_temperature. [@Toomre_stability1964] showed that for_an infinitesimally_thin disc to fragment, the_stability parameter must_be less_than a critical_value, $Q_{\rm crit} \approx 1$.

[@Gammie_betacool] showed_that in addition_to the stability criterion above, the_disc_must cool at_a_fast_enough rate._Using shearing sheet_simulations,_he showed_that_if the cooling timescale can be_parametrized_as

$$\label{eq:beta}
  \beta = t_{\rm cool}\Omega,$$

where

$$\label{eq:tcool}
 _t_{\rm cool} = u_\big(\frac{du_{\rm_cool}}{dt}\big)^{-1},$$

$u$ is the internal_energy and $du_{\rm cool}/dt$ is_the total cooling rate, then for_fragmentation we_require $\beta_\lesssim 3$, for a ratio of specific heats $\gamma = 2$_(in two dimensions). [-@Rice_beta_condition] carried out_three-dimensional simulations using a_smoothed particle_hydrodynamics_(SPH) code and_showed_that this_cooling parameter is dependent on the equation_of state_so that fragmentation can occur for_discs with $\gamma =_5/3$_and $7/5$ if $\beta \lesssim 7$_and $13$, respectively. [-@Libby_MSci] showed that_the critical value of $\beta$_(below_which_fragmentation will occur if the_stability