 in exactly the same way as a disc simulated by [@Lodato_Rice_original]: a 1 $\rm{M_\odot}$ star with a 0.1 $\rm{M_\odot}$ disc made of 250,000 SPH particles, spanning $0.25 \le R \le 25$AU. The initial surface mass density and temperature profiles of the disc are $\Sigma \propto R^{-1}$ and $T \propto R^{-\frac{1}{2}}$, respectively. The magnitudes of these are set such that the Toomre stability parameter (equation \[eq:Toomre\]) at the outer edge of the disc, $Q_{\rm min} = 2$. This gives an aspect ratio, $H/R \sim 0.05$. We model the $1~\rm{M_\odot}$ star in the centre of the disc using a sink particle [@Bate_Bonnell_Price_sink_ptcls]. At the inner disc boundary, particles are accreted onto the star if they move within a radius of 0.025 AU of the star or if they move into $0.025 \le R < 0.25 \rm{AU}$ and are gravitationally bound to the star. At the outer edge, the disc is free to expand.

The initial Reference disc is a Toomre stable disc. Given that the boundary temperature profile is the same as that of the initial disc and that $Q_{\rm min}=~2$, we do not expect this disc to fragment. However, we use the Reference disc as a fiducial case. In particular, we are concerned with the cooling rates that are present in the disc once it is in an equilibrium state. We emphasize our use of terminology here: when referring to the disc being in *thermal equilibrium with the boundary*, we are describing the bulk of the disc being a similar temperature to the disc boundary (which is assumed to be determined by stellar irradiation), whereas an *equilibrium state* disc refers to the dissipative and cooling rates being balanced such that the Toomre stability profiles do not change with time.

Exploring the parameter space {#sec:parameter}
-----------------------------

Given that a motivation of this work is to determine if, and under what circumstances, fragmentation in realistically modelled self-gravitating discs may occur, we explore the parameter space in a number of ways. One parameter is the opacity: we re-run the Reference simulation with opacity values scaled to $
in exactly the same way as a disc simulate by [ @Lodato_Rice_original ]: a 1   $ \rm{M_\odot}$ asterisk with a 0.1   $ \rm{M_\odot}$ disc made of 250,000 SPH particles, cross $ 0.25 \le R \le 25$AU. The initial surface mass density and temperature profile of the disc are $ \Sigma \propto R^{-1}$ and $ T \propto R^{-\frac{1}{2}}$, respectively. The magnitude of these are set such that the Toomre stability argument (equation   \[eq: Toomre\ ]) at the outer edge of the disc, $ Q_{\rm minute } = 2$. This gives an aspect proportion, $ H / R \sim 0.05$. We model the $ 1~\rm{M_\odot}$ star in the center of the phonograph record using a sink particle [ @Bate_Bonnell_Price_sink_ptcls ]. At the inside disc boundary, atom are accreted onto the ace if they move within a radius of 0.025   AU of the star or if they move into $ 0.025 \le R < 0.25 \rm{AU}$ and are gravitationally bound to the asterisk. At the outer edge, the disc is free to expand. 

 The initial Reference disc is a Toomre stable disc. Given that the boundary temperature profile is the same as that of the initial disc and that $ Q_{\rm min}=~2 $, we do not expect this disc to fragment. However, we use the Reference disc as a fiduciary event. In particular, we are implicated with the cool rate that are present in the disc once it is in an equilibrium department of state. We emphasize our use of terminology here: when referring to the disc being in * thermal chemical equilibrium with the boundary *, we are describing the bulk of the disc being a similar temperature to the disc boundary (which is assume to be determined by stellar irradiation), whereas an * equilibrium state * disc consult to the dissipative and cooling rates being balanced such that the Toomre stability profiles do not transfer with time. 

 Exploring the parameter quad { # sec: parameter } 
 ----------------------------- 

 Given that a motivation of this work is to determine if, and under what circumstance, fragmentation in realistically model self - gravitating discs may happen, we explore the parameter space in a number of way. One parameter is the opacity: we re - ply the Reference simulation with opacity values scale to $
 in exactly the same way as a disc simulatgd by [@Lovato_Rics_originau]: a 1 $\rm{M_\odot}$ star with a 0.1 $\rm{M_\ovot}$ eisc nade of 250,000 SPH particles, spanning $0.25 \le R \lw 25$AU. Rhe initial surfacc masa denvmty and temperajure profilev of the disc drd $\Digma \propto R^{-1}$ and $T \propto R^{-\frac{1}{2}}$, rqspectifepy. The magnituqes ps thsse are set such that the Toomre sfabilitj parameter (equatoon \[eq:Toomre\]) at the outer efge lf the disc, $Q_{\rm mij} = 2$. This gices wb aspect ratko, $H/R \sim 0.05$. We model thg $1~\rm{M_\odot}$ star in the centre of ghe dnsc using a sunk karticle [@Batx_Bonnejl_Price_sink_pbbls]. At dhe innrr disc boundavy, pacticoes are accreted onto the star if they iove withhn a radius of 0.025 AU od rhe sjar os if rheh mkvx ihto $0.025 \lf R < 0.25 \rm{AU}$ ans are gravirationally bound to trv star. At the outer eqge, the disc is free to expand.

The initidl Deference disc is a Toonre stable disc. Given that the boundary temperature profile is the same as that of tie innblal ausf and that $Q_{\rm min}=~2$, we do not expect this disc fo fgagment. However, wc use the Referencr fixs as a fiducixl casz. Ih particular, we arf concetned wuth the cjolimg rates that are present ib the disc ouce it is in an equilnbrium state. We gmphasoze our use of terminolugy gere: when rfferring fu the disc being in *tvermal equilibrium with thq boundarb*, we cre descfibimg the bulk of tje disc being a similar hempetature to the didc boundary (which is assumed to be determined bf snellar irxadiatlon), whereas an *qquilibrium stcte* disc refefs to the sissipavive and coojing rates behjg balanced vuch thae thw Toimre stxcility profilex do not bhcnge with time.

Exploring the patajeter space {#sec:'crqmeter}
-----------------------------

Given thay a moeinatmon os this work iv to degrrmind if, and unbef whst circumstances, fracmenfation in realistivajly modeoled sels-gravitating ciscs may occur, we explmre the larwmeter space in a number of wags. One pagamcter is the o[aciby: wg re-run the Reference simulation with opacity valuxs scaled to $
in exactly the same way as a by a 1 star with a 250,000 particles, spanning $0.25 R \le 25$AU. initial surface mass density and temperature of the disc are $\Sigma \propto R^{-1}$ and $T \propto R^{-\frac{1}{2}}$, respectively. The of these are set such that the Toomre stability parameter (equation \[eq:Toomre\]) at outer of disc, min} = 2$. This gives an aspect ratio, $H/R \sim 0.05$. We model the $1~\rm{M_\odot}$ star the centre of the disc using a sink [@Bate_Bonnell_Price_sink_ptcls]. At the inner boundary, particles are accreted onto star they move a of AU of the or if they move into $0.025 \le R < 0.25 \rm{AU}$ and are gravitationally bound to the At the the disc free expand. initial Reference disc Toomre stable disc. Given that the is the same as that of the initial and that min}=~2$, we do not expect this to fragment. However, we use the Reference disc a fiducial case. In particular, we are concerned with the cooling rates that are present disc once it is an equilibrium state. emphasize use terminology when referring the disc being in *thermal equilibrium with the boundary*, we are the bulk of the disc being a similar temperature to boundary is assumed to determined by stellar irradiation), an state* disc refers to and rates that Toomre profiles do not change time. Exploring the parameter space ----------------------------- Given that a to determine if, and under what circumstances, fragmentation realistically modelled self-gravitating discs may occur, we the parameter space in a number of ways. One parameter is the we re-run simulation with opacity values scaled to $
 in exactly the same way as a disC simulated By [@LodAto_ricE_oRigiNal]: a 1 $\Rm{M_\odot}$ star witH A 0.1 $\rm{M_\Odot}$ disc made of 250,000 SPH partiCles, sPaNNing $0.25 \LE R \Le 25$AU. THe initiAL sURFacE mAsS deNsITy And teMpeRature pRofiles of tHe dIsC are $\Sigma \proPTo r^{-1}$ and $T \proptO R^{-\fRac{1}{2}}$, respectivEly. the magNiTudES of thEse Are seT such tHAt the TOomre stabIlITy paraMEter (equATIoN \[eq:TOomre\]) at the outer edGE oF The disc, $Q_{\rm min} = 2$. THis givEs AN aSPEct RatIo, $H/R \sim 0.05$. We mOdEl the $1~\RM{M_\odot}$ sTAr IN THe cENtre of the disc Using a sink pARtiCle [@BatE_BOnnELl_PricE_sink_PtCLs]. AT the inner diSc boUndary, parTicles ARe accreTEd onto tHe star If tHey Move WItHiN a rAdIUs oF 0.025 aU Of tHE stAr or if thEy MoVe intO $0.025 \le R < 0.25 \RM{au}$ And aRe gRaviTatioNally bound to tHe sTar. AT The Outer Edge, tHe diSc Is freE to expAnd.

ThE iNitial Reference Disc Is a Toomre StaBlE diSc. given THat the BouNdaRy tempeRature pROfiLe IS THe Same as that of the iniTiAL DiSc and thaT $Q_{\rm miN}=~2$, We Do NOt expect ThIs dIsc tO FRagmeNt. HoWEvEr, we use tHe RefeREnCe Disc as a FiDucial CaSe. IN paRticuLAr, we Are conCerned wiTh the COoling rates thaT Are present in tHE dISC oNCe it Is iN an equilibrIum sTAte. WE empHAsIze OUr use Of terMiNOlOGy here: when referring To The disC beinG in *thermal equIlibrium wiTH THe boundaRy*, we ARe DEscribing the buLk of tHe disc beinG A similar TempeRature to The disc boUNDary (whicH is AssUmeD to BE DeTermined by steLLAr irRaDiation), WheReas an *eQuiLibRiuM stAtE* disc refeRs to the dIsSiPaTiVe aNd cooLIng rates BeIng BaLanCed suCH that tHe TooMre sTaBiLIty ProfileS Do NOT chaNgE wIth tIme.

exPloriNg thE ParAmeter sPace {#sec:paRamETer}
-----------------------------

GIvEn That a moTivation of thiS wOrk is to detErMinE if, and UNDer what cIrcumstances, fragmentatiON in realIstIcallY modElled self-GraVitatiNg dIScs may Occur, wE explOrE thE PArameTER sPacE iN a number of WAYs. ONe parAmEter Is the opAcity: we re-run the RefERenCe simulation wIth OpacITY vAluES sCAleD tO $
 in exactly the same way a s a disc s imula ted  by  [ @Lod ato_ Rice_original] :  a 1  $\rm{M_\odot}$ star w ith a  0 . 1 $\ r m{ M_\od ot}$ di s cm a deof  2 50, 00 0  S PH pa rti cles, s panning $0 .25  \ le R \le 25$ A U.  The initi alsurface mass  de nsityan d t e mpera tur e pro fileso f thedisc are$\ S igma \ p ropto R ^ { -1 }$ a nd $T \propto R^{ - \f r ac{1}{2}}$, re specti ve l y. T hemag nitudes of  t hesea re sets uc h t hat the Toomre st ability par a met er (eq ua tio n  \[eq: Toomr e\ ] ) a t the outer  edg e of thedisc,$ Q_{\rmm in} = 2 $. Thi s g ive s an as pe ctra t io, $H /R\ sim  0.05$.We  m odelthe$ 1 ~ \ rm{M _\o dot} $ sta r in the cent reof t h e d isc u singa si nk  part icle [ @Bate _B onnell_Price_si nk_p tcls]. At  th einn er  disc bounda ry,  pa rticles  are ac c ret ed o n to  the star if theymo v e  w ithin aradius of  0 . 025 AU o fthe  sta r or if  the y  m ove into  $0.02 5  \ le  R < 0. 25  \rm{A U} $ a ndare g r avit ationa lly boun d tot he star. At th e  outer edge,t he d is c  isfre e to expand .

T h e in itia l  R efe r encediscis aT oomre stable disc.Gi ven th at th e boundary te mperaturep r o file isthes am e  as that of th e ini tial disca nd that$Q_{\ rm min}= ~2$, we d o not expe ctthi s d isc t ofragment. How e v er,we  use th e R eferenc e d isc  as  afi ducial ca se. In p ar ti cu la r,we ar e  concern ed  wi th  th e coo l ing ra tes t hatar ep res ent int he d iscon ce  itisin  an e quil i bri um stat e. We emp has i ze o ur  u se of t erminology he re : when ref er rin g to t h e  disc be ing in *thermal equilib r ium wit h t he bo unda ry*, we a redescri bin g  the b ulk of  thedi scb e ing a s im ila rtemperatur e tothe d is c bo undary(which is assumedt o b e determinedbystel l a rirr a di a tio n) ,  wh e r eas an *equilib rium state *d is c refers t o  th edissipa tive an d coo l ing rat es beingbalancedsu ch t h a t t he Toomrestabilit y profile s  do n o tchang e w ith ti me .
Explo ring t h e p arame ter sp ac e {#se c:par am eter}
-- ----------------------- ----
Given  th at a moti vat i onof this w orkis to dete rmi neif, a ndu nderwhat ci rcu m stanc es,f ragmentat i on  in r ea listicallym o d ell ed se lf- g ravita ting  discs may occur, we explore the  par a m ete r s p acein  a number of w ays .O n e parame te r is the op acity: w er e-run  the R eferen ce simu l a ti o n with  opa cit y valuessca le d  to $
 in_exactly the_same way as a_disc simulated_by_[@Lodato_Rice_original]: a_1 $\rm{M_\odot}$_star with a_0.1 $\rm{M_\odot}$ disc made_of 250,000 SPH particles,_spanning $0.25 \le_R_\le 25$AU. The initial surface mass density and temperature profiles of the disc are_$\Sigma_\propto R^{-1}$_and_$T_\propto R^{-\frac{1}{2}}$, respectively. The magnitudes_of these are set such_that the_Toomre stability parameter (equation \[eq:Toomre\]) at the outer edge_of_the disc, $Q_{\rm_min} = 2$. This gives an aspect ratio, $H/R_\sim 0.05$. We model the $1~\rm{M_\odot}$_star in the_centre_of_the disc using a_sink particle [@Bate_Bonnell_Price_sink_ptcls]. At the inner_disc boundary, particles are accreted onto_the star if they move within a_radius of 0.025 AU of the star_or if they move into_$0.025 \le_R < 0.25 \rm{AU}$ and_are gravitationally bound_to the_star. At the_outer edge, the disc is free_to expand.

The initial_Reference disc is a Toomre stable_disc._Given that the_boundary_temperature_profile is_the same as_that_of the_initial_disc and that $Q_{\rm min}=~2$, we_do_not expect this disc to fragment. However,_we use the Reference_disc_as a fiducial case._In particular, we are concerned_with the cooling rates that are_present in_the disc_once it is in an equilibrium state. We emphasize our use_of terminology here: when referring to_the disc being in_*thermal equilibrium_with_the boundary*, we_are_describing the_bulk of the disc being a similar_temperature to_the disc boundary (which is assumed_to be determined by_stellar_irradiation), whereas an *equilibrium state* disc_refers to the dissipative and cooling_rates being balanced such that_the_Toomre_stability profiles do not change_with time.

Exploring the parameter space {#sec:parameter}
-----------------------------

Given_that a motivation_of this work is to determine if,_and_under what circumstances, fragmentation in realistically_modelled_self-gravitating discs may occur, we explore_the_parameter_space in a number of_ways. One parameter is the opacity:_we re-run the Reference simulation with opacity values scaled_to $