 Jeans length and mass, making it numerically more prohibitive to resolve these critical scales. For the simulations presented here, we are careful to use sufficient numbers of particles so that the Jeans length and mass are always adequately resolved. The higher gas fraction simulations require a higher particle number as a result (see Table \[TableRunParam\]). The reduction in sound speed also makes it more likely that the gas will fragment by gravitational instability into clumps (ala molecular clouds), as we shall discuss in detail later. This fragmentation is real, not numerical; artificially increasing the sound speed to eliminate it is not necessarily physical and could give incorrect results. On the other hand, we do not include sufficient physics in our ISM model to describe the formation and disruption of molecular clouds so our treatment of the resulting clumping is also not correct. In §\[sec:ISM\] we discuss which of our results are the most sensitive to uncertainties related to local gravitational instability in the ISM.

Black Hole Accretion and Feedback {#sectionBlackHoles}
---------------------------------

### Black Hole Accretion Model

We include a BH as an additional collisionless particle at the center of each galaxy. We model the accretion of the surrounding gas onto the BH, via the transport of angular momentum, using $$\dot{M}_{visc} = 3 \pi \alpha \Sigma \frac{c_s^2}{\Omega}
\label{mdotvisceqn}$$ where $\Sigma$ is the mean gas surface density, $\Omega$ is the rotational angular frequency, and $\alpha$ is the dimensionless viscosity (a free parameter of our model). We compute $\Sigma$ and $c_s$ by taking an average of the properties of the SPH particles in a sphere of radius $R_{acc}$ centred on the BH. The radius $R_{acc}$ is typically set equal to four times the gravitational force softening length, i.e., $R_{acc} = 4 \epsilon$, although we explore alternate choices as well. We find that estimating the rotation rate using $\Omega^2 \simeq GM(<R_{acc})/R_{acc}^3$ is more numerically robust than actually calculating the rotation and angular momentum of the gas particles within $R_{acc}$.

Although equation (\[mdotvisceqn\]) is reminiscent of the alpha prescription of, in our formulation $\alpha$ characterizes not only the efficiency of angular momentum transport, but also the uncertainty due to the fraction of the inflowing gas
Jeans length and mass, making it numerically more prohibitive to resolve these critical scale. For the model presented here, we are careful to use sufficient numbers of particle so that the Jeans length and mass are constantly adequately answer. The higher gasoline fraction simulations require a high particle number as a solution (determine Table \[TableRunParam\ ]). The reduction in sound speed also makes it more probable that the gas will fragment by gravitational instability into ball (ala molecular clouds), as we shall discuss in contingent later. This atomization is real, not numerical; artificially increasing the sound speed to eliminate it is not inevitably physical and could give incorrect results. On the other hand, we do not admit sufficient purgative in our ISM mannequin to describe the formation and disruption of molecular clouds so our treatment of the resulting clumping is also not right. In § \[sec: ISM\ ] we discourse which of our results are the about sensible to uncertainties related to local gravitational instability in the ISM. 

 Black Hole Accretion and Feedback { # sectionBlackHoles } 
 --------------------------------- 

 # # # Black Hole Accretion Model 

 We include a bohrium as an additional collisionless atom at the center of each galaxy. We model the accretion of the surrounding accelerator onto the BH, via the transport of angular momentum, using $ $ \dot{M}_{visc } = 3 \pi \alpha \Sigma \frac{c_s^2}{\Omega } 
 \label{mdotvisceqn}$$ where $ \Sigma$ is the mean gas surface density, $ \Omega$ is the rotational angular frequency, and $ \alpha$ is the dimensionless viscosity (a free parameter of our model). We compute $ \Sigma$ and $ c_s$ by taking an average of the property of the SPH atom in a sphere of radius $ R_{acc}$ centred on the BH. The radius $ R_{acc}$ is typically jell equal to four times the gravitational force yield length, i.e., $ R_{acc } = 4 \epsilon$, although we explore alternate choices as well. We line up that estimating the rotation rate using $ \Omega^2 \simeq GM(<R_{acc})/R_{acc}^3 $ is more numerically robust than actually calculating the rotation and angular momentum of the gasoline particles within $ R_{acc}$. 

 Although equality   (\[mdotvisceqn\ ]) is reminiscent of the alpha prescription of, in our conceptualization $ \alpha$ characterizes not only the efficiency of angular momentum transport, but also the uncertainty due to the fraction of the inflowing natural gas
 Jewns length and mass, makikg it numericalli nore pcohibitjve to rdsolve these critical scales. Fir tht simulations presevted here, we are xarehul to use suffirjent numbers kn parcirles so that thg Jeans lengdh and mass ara xlcays adequately resolved. The higher das fravtlon simulationf reatire a higher particle number as a reshlt (see Table \[TableRinParam\]). The reduction in slund speed also makes lt more likgmy eyat the gas dill fragmtnc by gravitztional instability into clumps (ala kolecular xliudd), as we shalo disbuss in detail later. Dhis frsgmentation is reel, nit numerical; artificielly increasing the found spead to eliminate it us not tecevsarkoy ohyairal and clulv give incodrect resulrs. On the other hanc, rv do not inclhde susfycient physics in our ISM model to descgibe the formation and disryption of molecular cpouds so jur treatment of the resulting clumping is also nmt cocrdct. Ik §\[sez:USL\] we discuss which of our results are the mose stnsptive to uncertaikties related to lpcwl dravitational instayjljty in the ISM.

Blacn Hole Wccrerion and Seednack {#sectionBlackHoles}
---------------------------------

### Blacj Hole Accrenion Model

We include a BH as an abditiomal cpllisionless particle ac the denter of ewch galaxg. We model the aczrenion of the surrounding gas oneo the BH, via the trxnspprt of angular mlmentmk, using $$\dot{M}_{visc} = 3 \pi \appva \Sigma \fgac{c_s^2}{\Omega}
\label{mdotvisceqn}$$ wherx $\Sigma$ is thg maan gas suryace dcnsity, $\Omega$ is the rotationak angulcr freduency, and $\alpha$ ms the dimenfionless viscmdity (a free 'arameter of iur nodel). Wd compute $\Sigma$ and $c_s$ bj uaking an qverage of the proiertigs of the SPH parcncoes in a sphere of raqils $C_{acc}$ santred on tha BH. Thd radijs $R_{acc}$ is bypkcalky set equal to four timss the gravitationsl force siftening length, i.e., $R_{avc} = 4 \epsilon$, althokgh wx explure aktetnate choices as well. We find tgat estimwtikg the rotatijn rqte using $\Omzga^2 \simeq GM(<R_{acc})/R_{acc}^3$ is more numerically robust than actuaoly calculating the votation and angujar momendum of the gas partixles within $R_{acc}$.

Akthough equation (\[mdotvjsceqn\]) is rfminiscent of the alpha prescription of, in our formulation $\alpha$ characteeizes ioe only the sffiviencj oy augular iomeitnm transport, but slso the uncertainty due to the fractiot jf the inflowing gas
Jeans length and mass, making it numerically to these critical For the simulations to sufficient numbers of so that the length and mass are always adequately The higher gas fraction simulations require a higher particle number as a result Table \[TableRunParam\]). The reduction in sound speed also makes it more likely that gas fragment gravitational into clumps (ala molecular clouds), as we shall discuss in detail later. This fragmentation is real, numerical; artificially increasing the sound speed to eliminate is not necessarily physical could give incorrect results. On other we do include physics our ISM model describe the formation and disruption of molecular clouds so our treatment of the resulting clumping is also correct. In discuss which our are most sensitive to to local gravitational instability in the Accretion and Feedback {#sectionBlackHoles} --------------------------------- ### Black Hole Model We a BH as an additional collisionless at the center of each galaxy. We model accretion of the surrounding gas onto the BH, via the transport of angular momentum, using 3 \pi \alpha \Sigma \label{mdotvisceqn}$$ where $\Sigma$ the gas density, is the angular frequency, and $\alpha$ is the dimensionless viscosity (a free parameter our model). We compute $\Sigma$ and $c_s$ by taking an the of the SPH in a sphere of $R_{acc}$ on the BH. The is set times gravitational softening length, i.e., $R_{acc} 4 \epsilon$, although we explore choices as well. We rate using $\Omega^2 \simeq GM(<R_{acc})/R_{acc}^3$ is more numerically than actually calculating the rotation and angular of the gas particles within $R_{acc}$. Although equation (\[mdotvisceqn\]) is reminiscent of alpha prescription our formulation $\alpha$ characterizes not only the efficiency angular momentum transport, but the uncertainty due to the fraction of the inflowing
 Jeans length and mass, making iT numericalLy morE prOhiBiTive To reSolve these critICal sCales. For the simulations PreseNtED herE, We Are caReful to USe SUFfiCiEnT nuMbERs Of parTicLes so thAt the Jeans LenGtH and mass are aLWaYs adequateLy rEsolved. The hiGheR gas frAcTioN SimulAtiOns reQuire a HIgher pArticle nuMbER as a reSUlt (see TABLe \[tablERunParam\]). The reducTIoN In sound speed alSo makeS iT MoRE LikEly That the gas WiLl fraGMent by gRAvITATioNAl instability Into clumps (aLA moLeculaR cLouDS), as we sHall dIsCUss In detail latEr. ThIs fragmenTation IS real, noT NumericAl; artiFicIalLy inCReAsIng ThE SouND sPeeD To eLiminate It Is Not neCessARILY phySicAl anD coulD give incorrecT reSultS. on tHe othEr hanD, we dO nOt incLude suFficiEnT physics in our ISm modEl to descrIbe ThE foRmAtion ANd disrUptIon Of molecUlar cloUDs sO oUR TReAtment of the resultiNg CLUmPing is alSo not cORrEcT. in §\[sec:ISM\] We DisCuss WHIch of Our rESuLts are thE most sENsItIve to unCeRtaintIeS reLatEd to lOCal gRavitaTional inStabiLIty in the ISM.

BlaCK Hole AccretioN AnD fEeDBack {#SecTionBlackHoLes}
---------------------------------

### BLAck HOle ACCrEtiON ModeL

We inClUDe A bH as an additional colLiSionleSs parTicle at the cenTer of each gALAXy. We modeL the ACcREtion of the surrOundiNg gas onto tHE BH, via thE tranSport of aNgular momENTum, using $$\Dot{m}_{viSc} = 3 \pI \alPHA \SIgma \frac{c_s^2}{\OmeGA}
\LabeL{mDotviscEqn}$$ Where $\SiGma$ Is tHe mEan GaS surface dEnsity, $\OmEgA$ iS tHe RotAtionAL angular FrEquEnCy, aNd $\alpHA$ is the DimenSionLeSs VIscOsity (a fREe PARameTeR oF our ModEl). we comPute $\sIgmA$ and $c_s$ bY taking an AveRAge oF tHe PropertIes of the SPH paRtIcles in a spHeRe oF radiuS $r_{Acc}$ centrEd on the BH. The radius $R_{acc}$ iS TypicalLy sEt equAl to Four times The GravitAtiONal forCe softEning LeNgtH, I.E., $R_{acc} = 4 \EPSiLon$, AlThough we exPLOre AlterNaTe chOices as Well. We find that estiMAtiNg the rotation RatE usiNG $\omEga^2 \SImEQ GM(<r_{aCC})/R_{aCC}^3$ Is more numericalLy robust thAn ACtUally calcuLAtiNg The rotaTion and AngulAR momentUm of the gaS particleS wIthiN $r_{Acc}$.

although eqUation (\[mdOtvisceqn\]) IS remiNIsCent oF thE alpha PrEscRiptiOn of, in OUr fOrmulAtion $\aLpHa$ charActerIzEs not onlY the efficiency of angular MomentUm traNspOrt, but alsO thE UncErtainty dUe to The fractioN of The InfloWinG Gas
 Jeans length and mass, ma king it nu meric all y m or e pr ohib itive to resol v e th ese critical scales. F or th es imul a ti ons p resente d  h e r e,we  a reca r ef ul to  us e suffi cient numb ers  o f particless othat the J ean s length and  ma ss are  a lwa y s ade qua telyresolv e d. The  higher g as fracti o n simul a t io ns r equire a higher p a rt i cle number asa resu lt (s e e  Ta ble  \[TableRu nP aram\ ] ). Ther ed u c t ion in sound spee d also make s  it  moreli kel y  thatthe g as wil l fragmentby g ravitatio nal in s tabilit y  into c lumps(al a m olec u la rclo ud s ),a swes hal l discus sin  deta il l a t e r . Th isfrag menta tion is real,  no t nu m eri cal;artif icia ll y inc reasin g the  s ound speed to e limi nate it i s n ot  ne ce ssari l y phys ica l a nd coul d givei nco rr e c t  r esults. On the oth er h an d, we do  not i n cl ud e  suffici en t p hysi c s  in o ur I S Mmodel to  descr i be  t he form at ion an ddis rup tiono f mo lecula r clouds  so o u r treatment of the resulting cl u m pi n g is  al so not corr ect. In § \[se c :I SM\ ]  we d iscus sw hi c h of our results ar ethe mo st se nsitive to un certaintie s r elated t o lo c al gravitationalinsta bility int he ISM.
Blac k Hole A ccretiona n d Feedba ck{#s ect ion B l ac kHoles}
----- - - ---- -- ------- --- ------- --- 

# ##Bla ck  Hole Acc retion M od el 
We  in clude a BH asan  ad di tio nal c o llisio nless  par ti cl e  at  the ce n te r of e ac hgala xy.  W e mod el t h e a ccretio n of thesur r ound in ggas ont o the BH, via  t he transpo rt  of  angul a r  momentu m, using $$\dot{M}_{vis c } = 3 \ pi\alph a \S igma \fra c{c _s^2}{ \Om e ga}
\l abel{m dotvi sc eqn } $ $ whe r e  $ \Si gm a$ is them e angas s ur face  densit y, $\Omega$ is the rot ational angul arfreq u e nc y,a nd $\a lp h a$i s  the dimensionl ess viscos it y  ( a free par a met er  of our  model) . Wec ompute$\Sigma$and $c_s$  b y ta k i ngan average  of thepropertie s  of t h eSPH p art iclesin  aspher e of r a diu s $R_ {acc}$  c entred  on t he  BH. The  radius $R_{acc}$ is ty picall y set  eq ual to fo urt ime s the gra vita tional for cesof tenin g l e ngth,  i.e . ,$R_ { acc}= 4\ epsilon$, al tho u g hwe explorea l t ern ate c hoi c es aswell . We find that es t imating the ro tati o n  ra teu sing  $ \Omega^2 \sime q G M( < R _{acc})/ R_ {acc}^3$ is  more nu me r icall y robu st tha n actua l l yc alcula ting  th e rotatio n a nd angular  m om e ntum o f th egas pa rticle s  wit h i n $R_{acc}$.

Al thoug h equat i on(\[md ot visceqn \ ]) i s reminisc ent of thealphapres cript ion of,  i n ourfor mu lation $\a l pha$ char acter izes no tonly  th e effi cien c y  of a ngul ar  mo mentum tr a n sp o rt ,b utalso  theun cert ainty due to the f rac t ion ofth e i n f lowing ga s 
 Jeans_length and_mass, making it numerically_more prohibitive_to_resolve these_critical_scales. For the_simulations presented here,_we are careful to_use sufficient numbers_of_particles so that the Jeans length and mass are always adequately resolved. The higher_gas_fraction simulations_require_a_higher particle number as a_result (see Table \[TableRunParam\]). The_reduction in_sound speed also makes it more likely that_the_gas will fragment_by gravitational instability into clumps (ala molecular clouds), as_we shall discuss in detail later._This fragmentation is_real,_not_numerical; artificially increasing the_sound speed to eliminate it is_not necessarily physical and could give_incorrect results. On the other hand, we_do not include sufficient physics in_our ISM model to describe_the formation_and disruption of molecular clouds_so our treatment_of the_resulting clumping is_also not correct. In §\[sec:ISM\] we_discuss which of_our results are the most sensitive_to_uncertainties related to_local_gravitational_instability in_the ISM.

Black Hole_Accretion_and Feedback_{#sectionBlackHoles}
---------------------------------

###_Black Hole Accretion Model

We include a_BH_as an additional collisionless particle at the_center of each galaxy._We_model the accretion of_the surrounding gas onto the_BH, via the transport of angular_momentum, using_$$\dot{M}_{visc} =_3 \pi \alpha \Sigma \frac{c_s^2}{\Omega}
\label{mdotvisceqn}$$ where $\Sigma$ is the mean gas_surface density, $\Omega$ is the rotational_angular frequency, and $\alpha$_is the_dimensionless_viscosity (a free_parameter_of our_model). We compute $\Sigma$ and $c_s$ by_taking an_average of the properties of the_SPH particles in a_sphere_of radius $R_{acc}$ centred on the_BH. The radius $R_{acc}$ is typically_set equal to four times_the_gravitational_force softening length, i.e., $R_{acc}_= 4 \epsilon$, although we explore_alternate choices as_well. We find that estimating the rotation_rate_using $\Omega^2 \simeq GM(<R_{acc})/R_{acc}^3$ is more_numerically_robust than actually calculating the rotation_and_angular_momentum of the gas particles_within $R_{acc}$.

Although equation (\[mdotvisceqn\]) is reminiscent of_the alpha prescription of, in our formulation $\alpha$ characterizes_not only the_efficiency of angular momentum transport,_but_also_the uncertainty due to the fraction of the inflowing gas