 of viscosity in the disc. This is supported by the result that the Poisson noise for the low-resolution models that we have performed multiple times is less than 5%, see table \[tb:gadget\_eda1\]. Furthermore, the ISM loading rate is robust for all our simulations and does not seem to be affected by numerical noise.

Modelling of the ISM
--------------------

In order to be able to discern the relevant physical processes from the simulations, we have modelled the ISM in a very idealized way: as a homogeneous gas with no clumps and no turbulence. To a certain extent we have modelled the reaction of the disc to a density gradient when the flow first hits the disc. In that case the size of the disc is determined by the ram pressure. If the disc were to encounter a region of gas with lower density, the surface area of the disc can extend to larger radii without being stripped and that would probably increase its effective cross section. Despite the lower density, the ISM loading rate could still be of the same order, because a lower ISM density allows a larger surface area of the disc to collect ISM. In that sense, accretion onto protoplanetary discs may be a self-regulating process. In a follow-up work we will perform simulations with different ISM densities and velocities to investigate if and how the ISM loading rate depends on these quantities.

In principle it would be possible to add angular momentum to the disc, such that it maintains or prolongs its lifetime, if the ISM has turbulence or substructures on scales at or slightly smaller than the disc scale. However, even in that case the net effect is unlikely to increase the overall lifetime, since the disc would quickly encounter a different part of the ISM where an adverse, disruptive configuration of substructure might exist.

Comparison to other work {#sec:dis_m09_eda1}
------------------------

A similar simulation has been performed by M09, but for a higher disc mass and lower velocity of the ISM with respect to the disc. In their work both the Bondi-Hoyle radius and the radius beyond which ram pressure would dominate are much larger than the radius of the disc. They used a higher resolution for the disc, initially $\sim 2.5 \times 10^5$ particles, which had 8 times the mass of their ISM particles. However, in their work the ISM also does not reside in the outer regions of the disc (see their Fig
of viscosity in the disc. This is supported by the result that the Poisson noise for the depleted - settlement models that we have performed multiple times is less than 5% , witness table \[tb: gadget\_eda1\ ]. Furthermore, the ISM loading pace is robust for all our simulations and does not seem to be affect by numerical noise. 

 Modelling of the ISM 
 -------------------- 

 In order to be able to spot the relevant physical processes from the simulations, we have modelled the doctrine in a very idealized way: as a homogeneous gas with no thump and no turbulence. To a certain extent we have modelled the chemical reaction of the disc to a density gradient when the flow foremost hits the disc. In that case the size of the disc is determined by the random-access memory pressure. If the disc were to encounter a region of gas with low density, the surface area of the disc can extend to large radii without being stripped and that would probably increase its effective cross section. Despite the lower density, the ISM loading rate could still be of the like ordering, because a lower doctrine concentration allows a larger surface area of the magnetic disk to collect ISM. In that sense, accretion onto protoplanetary discs may be a self - regulating procedure. In a follow - up work we will perform simulations with different ISM concentration and velocities to investigate if and how the ISM loading rate depends on these quantity. 

 In principle it would be potential to lend angular momentum to the disc, such that it maintains or prolongs its lifetime, if the ISM have turbulence or substructures on scales at or slightly little than the disk scale. However, even in that case the final effect is unlikely to increase the overall lifetime, since the disc would quickly meet a different region of the ISM where an adverse, disruptive configuration of infrastructure might exist. 

 Comparison to other oeuvre { # sec: dis_m09_eda1 } 
 ------------------------ 

 A similar simulation has been performed by M09, but for a higher disc batch and lower speed of the ISM with respect to the disc. In their work both the Bondi - Hoyle spoke and the radius beyond which ram pressure would dominate are a lot larger than the radius of the disc. They used a higher resolution for the disc, initially $ \sim 2.5 \times 10 ^ 5 $ particle, which had 8 times the mass of their ISM particles. However, in their work the ISM also does not occupy in the out regions of the disc (visualize their Fig
 of viscosity in the disc. Tmis is supported by the resulf that tfe Poisson noise for the low-cesooutiob models that we have oerformed multiplw tinws is less than 5%, scz tabmc \[tb:gcdjet\_eda1\]. Furthermpre, the ISK loading rate ir xobust for all our simulations and djes not sfem to be affested fy nhmerical noise.

Modelling of the ISM
--------------------

Jn ordeg to be able to doscern the relevant physicwl pgocesses from the dimulations, we rqve modelled the ISM in a very idezlized way: as a homogeneous gas with no clumps abd jm turbulencx. To a certain extcmt we vave mocelled the reagtion of the disc to a densitb gradient when the slow firsd kits the disc. In that cqse tve shze ud tfe sixc is deherjined by tge ram preswure. If the disc wete np encounter z regijn of gas with lower density, the surface drez of the disc can extene to larger radii witjout beind stripped and that would probably increase its exfectmvd cxiss sdxtlon. Despite the lower density, the ISM loading daue bould still be of the same otdfr, fecause a lowgr ISM bsnaity allows a largfr surfwce aeea of tht disv to collect ISM. In that sebse, accretiou obto protoplanetary discs may ye a sglf-regilating process. In a foulow-hp work we aill perfkfm simulations wkth dhfferent ISM densities and velocitixs to investkgatg if anq how the LSM loading rate depends on tkese xuantities.

Ln principle it would be possible to add angukas mmmentum co the disc, such thwt it maintainx or prjlongr its lifenime, if tie ISM has ttrbulence or vobstructures mn scalef at or wlightlh smaller than yhe disc scale. Howecer, even in that cese jhs net effect is unoikely to incresse thq lvxrall nifetime, sinwe tfe aosc wuuld quickln evcoumter a different pard of the ISM where an sdyerse, distuptive cjnfiguration pf substructure mihht eeist.

Cooparixon to other work {#sec:dis_m09_eda1}
------------------------

A simjlar simupatlon has been [erfirmed by M09, bbt for a higher disc mass and lower velority of the ISM with rgspect to the disc. In uheir work bovh the Bondi-Hoyne radius and the raeius beyond which ram pressure would dkminata are much larger than the radius of the disc. They used a higher resolution foe the vifc, initiallg $\sik 2.5 \tikef 10^5$ partycled, which had 8 timex the mass of their ISM particlxs. However, iu their work the ISM also doex vot reside in the outer regions of ths disc (sre their Fig
of viscosity in the disc. This is the that the noise for the performed times is less 5%, see table Furthermore, the ISM loading rate is for all our simulations and does not seem to be affected by numerical Modelling of the ISM -------------------- In order to be able to discern the physical from simulations, have modelled the ISM in a very idealized way: as a homogeneous gas with no clumps no turbulence. To a certain extent we have the reaction of the to a density gradient when flow hits the In case size of the is determined by the ram pressure. If the disc were to encounter a region of gas with density, the of the can to radii without being that would probably increase its effective the lower density, the ISM loading rate could be of same order, because a lower ISM allows a larger surface area of the disc collect ISM. In that sense, accretion onto protoplanetary discs may be a self-regulating process. In work we will perform with different ISM and to if how the loading rate depends on these quantities. In principle it would be to add angular momentum to the disc, such that it prolongs lifetime, if the has turbulence or substructures scales or slightly smaller than scale. even the effect unlikely to increase the lifetime, since the disc would encounter a different part adverse, disruptive configuration of substructure might exist. Comparison other work {#sec:dis_m09_eda1} ------------------------ A similar simulation been performed by M09, but for a higher disc mass and lower of the respect to the disc. In their work both Bondi-Hoyle radius and the beyond which ram pressure would dominate are much larger the of the They used a resolution for the initially $\sim 2.5 particles, which 8 the ISM particles. However, in their work ISM does not reside in the regions disc (see their Fig
 of viscosity in the disc. This iS supported By the ResUlt ThAt thE PoiSson noise for thE Low-rEsolution models that we hAve peRfORmed MUlTiple Times is LEsS THan 5%, SeE tAblE \[tB:GaDget\_eDa1\]. FUrthermOre, the ISM lOadInG rate is robusT FoR all our simUlaTions and does Not Seem to Be AffECted bY nuMericAl noisE.

modellIng of the IsM
--------------------

iN order TO be able TO DiScerN the relevant physiCAl PRocesses from thE simulAtIOnS, WE haVe mOdelled the iSm in a vERy idealIZeD WAY: as A Homogeneous gaS with no clumPS anD no turBuLenCE. To a ceRtain ExTEnt We have modelLed tHe reactioN of the DIsc to a dENsity grAdient WheN thE floW FiRsT hiTs THe dISc. in tHAt cAse the siZe Of The diSc is DETERminEd bY the Ram prEssure. If the diSc wEre tO EncOunteR a regIon oF gAs witH lower DensiTy, The surface area oF the Disc can exTenD tO laRgEr radII withoUt bEinG strippEd and thAT woUlD PRObAbly increase its effEcTIVe Cross secTion. DeSPiTe THe lower dEnSitY, the ism loadIng rATe Could stiLl be of THe SaMe order, BeCause a LoWer iSM DensiTY allOws a laRger surfAce arEA of the disc to coLLect ISM. In that SEnSE, AcCRetiOn oNto protoplaNetaRY disCs maY Be A seLF-reguLatinG pROcESs. In a follow-up work we WiLl perfOrm siMulations with Different Ism DEnsities And vELoCIties to investiGate iF and how the isM loadinG rate Depends oN these quaNTIties.

In pRinCipLe iT woULD bE possible to adD ANgulAr MomentuM to The disc, SucH thAt iT maInTains or prOlongs itS lIfEtImE, if The ISm Has turbuLeNce Or SubStrucTUres on ScaleS at oR sLiGHtlY smalleR ThAN The dIsC sCale. howEvEr, eveN in tHAt cAse the nEt effect iS unLIkelY tO iNcrease The overall lifEtIme, since thE dIsc Would qUICkly encoUnter a different part of thE iSM wherE an AdverSe, diSruptive cOnfIguratIon OF substRucturE mighT eXisT.

cOmparISOn To oThEr work {#sec:dIS_M09_edA1}
------------------------

A simIlAr siMulatioN has been performed bY m09, buT for a higher diSc mAss aND LoWer VElOCitY oF The ism with respect to tHe disc. In thEiR WoRk both the BONdi-hoYle radiUs and thE radiUS beyond Which ram pRessure woUlD domINAte Are much larGer than tHe radius oF The diSC. THey usEd a Higher ReSolUtion For the DIsc, InitiAlly $\siM 2.5 \tImes 10^5$ paRticlEs, Which had 8 Times the mass of their ISM pArticlEs. HowEveR, in their wOrk THe IsM also doeS not Reside in thE ouTer RegioNs oF The diSc (seE ThEir fIg
 of viscosity in the disc.  This is s uppor ted  by  t he r esul t that the Poi s sonnoise for the low-reso lutio nm odel s  t hat w e havep er f o rme dmu lti pl e  t imesisless th an 5%, see  ta bl e \[tb:gadge t \_ eda1\]. Fu rth ermore, theISM  loadi ng  ra t e isrob ust f or all our si mulations  a n d does not see m to  beaffected by numer i ca l  noise.

Model ling o ft he I SM--- ---------- -- ----- 
In ord e rt o bea ble to discer n the relev a ntphysic al  pr o cesses  from  t h e s imulations,  wehave mode lled t h e ISM i n  a very  ideal ize d w ay:a sahom og e neo u sgas wit h no clu mp sand n o tu r b u l ence . T o acerta in extent wehav e mo d ell ed th e rea ctio nof th e disc  to a  d ensity gradient  whe n the flo w f ir sthi ts th e  disc.  In  th at case  the si z e o ft h e  d isc is determinedby t he  ram pre ssure. If  t h e disc w er e t o en c o unter  a r e gi on of ga s with lo we r densi ty , thesu rfa ceareao f th e disc  can ext end t o  larger radiiw ithout beings tr i p pe d  and  th at would pr obab l y in crea s eits effec tivecr o ss section. Despite th elowerdensi ty, the ISM l oading rat e c ould sti ll b e  o f  the same orde r, be cause a lo w er ISM d ensit y allows  a larger s urface a rea  of  th e d i s cto collect IS M .  Inth at sens e,accreti onont o p rot op lanetarydiscs ma ybe  a  s elf -regu l ating pr oc ess .Ina fol l ow-upworkwe w il lp erf orm sim u la t i onswi th  dif fer en t ISM  den s iti es andvelocitie s t o  inv es ti gate if  and how theIS M loadingra tedepend s on these  quantities.

In princi p le it w oul d beposs ible to a ddangula r m o mentum  to th e dis c,  su c h  that i tmai nt ains or pr o l ong s its  l ifet ime, if  the ISM has turbu l enc e or substruc tur es o n sc ale s  a t  or  s l igh t l y smaller thanthe disc s ca l e.  However,e ven  i n thatcase th e net effectis unlike ly to inc re aset h e o verall lif etime, s ince thed isc w o ul d qui ckl y enco un ter  a di fferen t  pa rt of  the I SM  where  an a dv erse, di sruptive configurationof sub struc tur e might e xis t .
Compariso n to  other wor k { #se c:dis _m0 9 _eda1 }
-- - -- --- - ----- ---- - -----

As im ila r si mulation ha s b een  perf orm e d by M 09,but for a higherd isc mass and l ower v elo cit y  ofth e ISM with res pec tt o  the dis c.  In their w ork both  t h e Bon di-Hoy le rad ius and t he radius  bey ond  which ra m p re s sure wo ul dd ominat e ar emuch l argert hant h e radius of thedisc. T hey u s eda hig he r resol u tion  for the d isc, initia lly $\ sim2.5 \ times 1 0^ 5$ par tic le s, which h a d 8 times  themass of  t heir  IS M part icle s .  Howe ver,  i n t heir work t he IS Ma lso  doe s not  r esid e in theo uter reg ion s  of the  d isc ( see th e ir F ig
 of_viscosity in_the disc. This is_supported by_the_result that_the_Poisson noise for_the low-resolution models_that we have performed_multiple times is_less_than 5%, see table \[tb:gadget\_eda1\]. Furthermore, the ISM loading rate is robust for all_our_simulations and_does_not_seem to be affected by_numerical noise.

Modelling of the ISM
--------------------

In_order to_be able to discern the relevant physical processes_from_the simulations, we_have modelled the ISM in a very idealized way:_as a homogeneous gas with no_clumps and no_turbulence._To_a certain extent we_have modelled the reaction of the_disc to a density gradient when_the flow first hits the disc. In_that case the size of the_disc is determined by the_ram pressure._If the disc were to_encounter a region_of gas_with lower density,_the surface area of the disc_can extend to_larger radii without being stripped and_that_would probably increase_its_effective_cross section._Despite the lower_density,_the ISM_loading_rate could still be of the_same_order, because a lower ISM density allows_a larger surface area_of_the disc to collect_ISM. In that sense, accretion_onto protoplanetary discs may be a_self-regulating process._In a_follow-up work we will perform simulations with different ISM densities and_velocities to investigate if and how_the ISM loading rate_depends on_these_quantities.

In principle it_would_be possible_to add angular momentum to the disc,_such that_it maintains or prolongs its lifetime,_if the ISM has_turbulence_or substructures on scales at or_slightly smaller than the disc scale._However, even in that case_the_net_effect is unlikely to increase_the overall lifetime, since the disc_would quickly encounter_a different part of the ISM where_an_adverse, disruptive configuration of substructure might_exist.

Comparison_to other work {#sec:dis_m09_eda1}
------------------------

A similar simulation_has_been_performed by M09, but for_a higher disc mass and lower_velocity of the ISM with respect to the disc._In their work_both the Bondi-Hoyle radius and_the_radius_beyond which ram pressure would dominate are much larger than_the radius_of the disc._They used a higher resolution for the disc, initially $\sim_2.5 \times 10^5$ particles, which had 8_times the mass of their ISM particles. However, in their work_the ISM also does not reside in the_outer regions of the disc (see their_Fig