 significantly improved in future work, allowing a more detailed comparison to observations. For example, @hq10 carried out a large number of simulations of gas inflow in galactic nuclei from $\sim 100$ pc to $\lesssim 0.1$ pc (see, e.g., @levine10 for related work). These can be used to provide a more accurate estimate of the BH accretion rate given conditions at larger radii in a galaxy (Hopkins & Quataert, in prep). Another important improvement would be to use a radiative transfer calculation to self-consistently determine the infrared radiation field produced by a central AGN (and distributed star formation). This could then be used to calculate the radiation pressure force on surrounding gas, eliminating the need for our parameterization of the force in terms of the optical depth $\tau$.

Acknowledgments {#acknowledgments.unnumbered}
===============

We thank Phil Hopkins and Yuval Birnboim for useful conversations. JD and EQ were supported in part by NASA grant NNG06GI68G and the David and Lucile Packard Foundation. Support for EQ was also provided in part by the Miller Institute for Basic Research in Science, University of California Berkeley. This research used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. This work was partially supported by the National Center for Supercomputing Applications under AST080048 and utilized the Intel 64 cluster Abe. The authors acknowledge the Texas Advanced Computing Center (TACC) at The University of Texas at Austin for providing HPC resources that have contributed to the research results reported within this paper.

Resolution Studies {#resolution}
==================

![image](FigureA1.eps){width="180mm"}

In this section, we describe some of our resolution tests both with and without BH feedback. In the absence of feedback, the well-posed questions for resolution studies include both how the gas properties as a function of radius and time depend on the resolution and how integrated properties of the galaxy (e.g., the star formation rate) depend on resolution. However, the feedback, when present, has a nontrivial dependence on the resolution and it is by no means clear that the nonlinear system will in fact converge in a simple way with increasing resolution. Physically, e.g., the AGN’s radiation pressure has the strongest effect on the gas that contributes the most to the optical depth,
significantly improved in future work, allowing a more detailed comparison to observation. For case, @hq10 carried out a large number of simulation of gas inflow in galactic nucleus from $ \sim 100 $ personal computer to $ \lesssim 0.1 $ pc (see, e.g., @levine10 for related work). These can be used to leave a more accurate estimate of the BH accretion rate give conditions at larger radii in a galax (Hopkins & Quataert, in prep). Another important improvement would be to use a radiative transportation calculation to self - consistently determine the infrared radiation field produced by a cardinal AGN (and distributed star topology formation). This could then be used to calculate the radiation sickness pressure force on surrounding gas, eliminating the need for our parameterization of the force in terms of the optical depth $ \tau$. 

 Acknowledgments { # acknowledgments.unnumbered } 
 = = = = = = = = = = = = = = = 

 We thank Phil Hopkins and Yuval Birnboim for useful conversations. JD and EQ were supported in part by NASA grant NNG06GI68 G and the David and Lucile Packard Foundation. Support for EQ was besides supply in part by the Miller Institute for Basic Research in Science, University of California Berkeley. This inquiry use resource of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE - AC02 - 05CH11231. This work was partially supported by the National Center for Supercomputing Applications under AST080048 and utilized the Intel 64 cluster Abe. The authors acknowledge the Texas Advanced Computing Center (TACC) at The University of Texas at Austin for provide HPC resources that have contributed to the research results reported within this paper. 

 Resolution Studies { # resoluteness } 
 = = = = = = = = = = = = = = = = = = 

! [ image](FigureA1.eps){width="180 mm " } 

 In this section, we describe some of our resolution test both with and without BH feedback. In the absence of feedback, the well - put questions for resolution studies include both how the gas property as a function of radius and time depend on the resolution and how integrated properties of the galaxy (for example, the star formation rate) depend on resolution. However, the feedback, when present, have a nontrivial dependence on the resolution and it is by no means clean that the nonlinear system will in fact converge in a dim-witted way with increasing resolution. Physically, e.g., the AGN ’s radiation press has the firm effect on the gas that contributes the most to the optical depth,
 sihnificantly improved in nuture work, allocung a kore dstailed zomparison to observations. Flr examkje, @hq10 carried out x large nlmber of wimuoqtions of jzs inflow in fwlaccir nuclei from $\slm 100$ pc to $\lasssim 0.1$ pc (see, e.e., @pevine10 for related work). These can bq used yo provide a morg acclrwte sstimate of the BH accretion rate fiven cmnditions at karger radii in a galaxy (Hlpkijs & Quataert, in prfp). Another umpowrant improveoent would be to use a radiative transfer calculation to szlf-consistebtoy fgtermine the infrwred radiation field producrd by a centrak AJN (abd distributed star fmrmation). This coulq then be uaed to calculate rhw radhatimn pfwssjre fprde on durcounding gaa, eliminatibg the need for our pwgsmeterizatioh of tre force in terms of the optical depth $\tal$.

Ackhowledgments {#acknowledgnents.unnumbered}
===============

We thajk Phil Hjpkins and Yuval Birnboim for useful conversationv. JD ena EW rdee supported in part by NASA grant NNG06GI68G and trs Canid and Lucile Pagkard Foundation. Sipootj for EQ was auso provised in part by the Miller Instutute for Basoc Research in Science, Univwrsity of Cajufornia Berkeley. Tkis research usec respurces of the National Znergy Research Sfientific Zomputing Center, whpch hs supported by the Office of Scienre of the U.S. Depsrtmene of Energj undcs Contract No. DE-AC02-05FH11231. Thns wosk was parhially supported by the National Center for Sopescokputing Cpplicstions under WST080048 and utilizgd the Incel 64 cuuster Abe. The auvhors acknowjedge the Texdd Advanced Cmmputing Cenrer (RACC) at Ghe University of Texas at Austin for providing HPC rerkurces that havz xontributed to yhe refewrrh refglts reporteg wighiv this paper.

Rtfoojtiom Studies {#resolution}
==================

![hmags](FigureA1.eps){width="180mm"}

On this sextion, we describe somr of our resolutioj tesvs boti with anq without BH feedback. In the agsence of fecdback, the wejl-powed questionx for resolution studies include both hiw the gas properjies as a function of vadius and tmme de[end on tve resolution and hoq integrated proptrties of the galaxy (e.g., the sdar flrmation rate) depend on resolution. However, the feedback, when present, has q nontcidial dependsnce on tve xesolutyon env it is by no meams clear that the nonlinear sysvem will it yact converge in a simple way wkth increasine resolution. Physically, s.g., the ABN’s radiation pressure has the xtrongest effect on the nas that xontrjbutes the most to the opticsl depch,
significantly improved in future work, allowing a comparison observations. For @hq10 carried out of inflow in galactic from $\sim 100$ to $\lesssim 0.1$ pc (see, e.g., for related work). These can be used to provide a more accurate estimate the BH accretion rate given conditions at larger radii in a galaxy (Hopkins Quataert, prep). important would be to use a radiative transfer calculation to self-consistently determine the infrared radiation field produced a central AGN (and distributed star formation). This then be used to the radiation pressure force on gas, the need our of force in terms the optical depth $\tau$. Acknowledgments {#acknowledgments.unnumbered} =============== We thank Phil Hopkins and Yuval Birnboim for useful conversations. and EQ in part NASA NNG06GI68G the David and Foundation. Support for EQ was also by the Miller Institute for Basic Research in University of Berkeley. This research used resources of National Energy Research Scientific Computing Center, which is by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. was partially supported by National Center for Applications AST080048 utilized Intel 64 Abe. The authors acknowledge the Texas Advanced Computing Center (TACC) at University of Texas at Austin for providing HPC resources that to research results reported this paper. Resolution Studies ================== In this section, we of resolution and BH In the absence of the well-posed questions for resolution include both how the of radius and time depend on the resolution how integrated properties of the galaxy (e.g., star formation rate) depend on resolution. However, the feedback, when present, has nontrivial dependence resolution and it is by no means clear the nonlinear system will fact converge in a simple way with increasing resolution. e.g., AGN’s radiation has the strongest on the gas contributes the most optical depth,
 significantly improved in fuTure work, alLowinG a mOre DeTailEd coMparison to obseRVatiOns. For example, @hq10 carried Out a lArGE numBEr Of simUlationS Of GAS inFlOw In gAlACtIc nucLei From $\sim 100$ Pc to $\lesssiM 0.1$ pc (SeE, e.g., @levine10 for RElAted work). ThEse Can be used to pRovIde a moRe AccURate eStiMate oF the BH ACcretiOn rate givEn COnditiONs at larGER rAdii In a galaxy (Hopkins & QUAtAErt, in prep). AnothEr impoRtANt IMProVemEnt would be To Use a rADiative TRaNSFEr cALculation to seLf-consistenTLy dEtermiNe The INfrareD radiAtIOn fIeld produceD by a Central AGn (and diSTributeD Star forMation). thiS coUld tHEn Be UseD tO CalCUlAte THe rAdiation PrEsSure fOrce ON SURrouNdiNg gaS, elimInating the neeD foR our PAraMeterIzatiOn of ThE forcE in terMs of tHe Optical depth $\tau$.

acknOwledgmenTs {#aCkNowLeDgmenTS.unnumBerEd}
===============

WE thank PHil HopkINs aNd yUVAl birnboim for useful cOnVERsAtions. JD And EQ wERe SuPPorted in PaRt bY NASa GRant NnG06GI68g AnD the DaviD and LuCIlE PAckard FOuNdatioN. SUppOrt For EQ WAs alSo provIded in paRt by tHE Miller InstituTE for Basic ReseARcH IN SCIencE, UnIversity of CAlifORnia berkELeY. ThIS reseArch uSeD ReSOurces of the National enErgy ReSearcH Scientific CoMputing CenTER, Which is sUppoRTeD By the Office of SCiencE of the U.S. DePArtment oF EnerGy under COntract No. de-aC02-05CH11231. This WorK waS paRtiALLy Supported by thE nAtioNaL Center For supercoMpuTinG ApPliCaTions undeR AST080048 and uTiLiZeD tHe INtel 64 cLUster Abe. thE auThOrs AcknoWLedge tHe TexAs AdVaNcED CoMputing cEnTER (TACc) aT THe UnIveRsIty of texaS At AUstin foR providinG HPc ResoUrCeS that haVe contributed To The researcH rEsuLts repORTed withiN this paper.

Resolution StuDIes {#resoLutIon}
==================

![imAge](FIgureA1.eps){WidTh="180mm"}

In ThiS SectioN, we desCribe SoMe oF OUr resOLUtIon TeSts both witH ANd wIthouT Bh feeDback. In The absence of feedbaCK, thE well-posed queStiOns fOR ReSolUTiON stUdIEs iNCLude both how the gAs propertiEs AS a Function of RAdiUs And time Depend oN the rESolutioN and how inTegrated pRoPertIES of The galaxy (e.G., the star Formation RAte) dePEnD on reSolUtion. HOwEveR, the fEedbacK, WheN presEnt, has A nOntrivIal dePeNdence on The resolution and it is by nO means Clear ThaT the nonliNeaR SysTem will in Fact Converge in A siMplE way wIth INcreaSing REsOluTIon. PhYsicALly, e.g., the Agn’s RadIATiOn pressure hAS THe sTrongEst EFfect oN the Gas that contributeS The most to the opTicaL DEptH,
 significantly improved in  future wo rk, a llo win ga mo re d etailed compar i sonto observations. For e xampl e, @hq1 0  c arrie d out a la r g e n um be r o fs im ulati ons  of gas  inflow in  ga la ctic nucleif ro m $\sim 10 0$pc to $\less sim  0.1$pc  (s e e, e. g.,  @lev ine10f or rel ated work ). Thesec an be u s e dto p rovide a more acc u ra t e estimate ofthe BH  a c cr e t ion  ra te given c on ditio n s at la r ge r r adi i  in a galaxy(Hopkins &Q uat aert,in  pr e p). An other  i m por tant improv emen t would b e to u s e a rad i ative t ransfe r c alc ulat i on  t o s el f -co n si ste n tly  determi ne  t he in frar e d r adia tio n fi eld p roduced by acen tral AGN  (and  dist ribu te d sta r form ation ).  This could the n be  used tocal cu lat ethe r a diatio n p res sure fo rce ons urr ou n d i ng  gas, eliminatingth e ne ed for o ur par a me te r izationof  th e fo r c e interm s  o f the op ticald ep th  $\tau$ .
Ackno wl edg men ts {# a ckno wledgm ents.unn umber e d}
=========== = ===

We thank Ph i l  H o pkin s a nd Yuval Bi rnbo i m fo r us e fu l c o nvers ation s. JD and EQ were support ed  in pa rt by  NASA grant N NG06GI68Ga n d  the Dav id a n dL ucile PackardFound ation. Sup p ort forEQ wa s also p rovided i n part bythe  Mi lle r I n s ti tute for Basi c Rese ar ch in S cie nce, Un ive rsi tyofCa liforniaBerkeley .Th is  r ese archu sed reso ur ces  o f t he Na t ionalEnerg y Re se ar c h S cientif i cC o mput in gCent er,  w hichis s u ppo rted by  the Offi ceo f Sc ie nc e of th e U.S. Depart me nt of Ener gy  un der Co n t ract No.  DE-AC02-05CH11231. Thi s  work w asparti ally  supporte d b y theNat i onal C enterfor S up erc o m putin g Ap pli ca tions unde r AST 08004 8andutilize d the Intel 64 clu s ter  Abe. The aut hor s ac k n ow led g et heTe x asA d vanced Computin g Center ( TA C C)  at The Un i ver si ty of T exas at  Aust i n for p rovidingHPC resou rc es t h a t h ave contri buted to  the rese a rch r e su lts r epo rted w it hin  this  paper . 

R esolu tion S tu dies { #reso lu tion}
== ================

![ima ge](Fi gureA 1.e ps){width ="1 8 0mm "}

In th is s ection, we  de scr ibe s ome of ou r re s ol uti o n tes ts b o th with a n dwit h o ut  BH feedbac k . Inthe a bse n ce offeed back, the well-po s ed questions f or r e s olu tio n  stu di es include bot h h ow t he gas p ro perties asa functi on of ra dius a nd tim e depen d on the re solu tio n and how  in te g rated p ro pe r ties o f th egalaxy  (e.g. ,  the s tar formation ra te) d e p end o n  re solut io n. Howe v er,the feedba ck, when pr esent,  has  a no ntrivia ldepend enc eon the res o lution an d itis by n omean s c lear t hatt h e non line ar  sy stem will i nf ac tc onv erge  in a  s impl e way wit h  increas ing resolut io n.P h ysical l y, e .g., the A GN’ s rad i a tion press u re h a st he st ronges t effe ct on t h e g as  that c ont r i butes the  most tot heop tica l depth, 
 significantly_improved in_future work, allowing a_more detailed_comparison_to observations._For_example, @hq10 carried_out a large_number of simulations of_gas inflow in_galactic_nuclei from $\sim 100$ pc to $\lesssim 0.1$ pc (see, e.g., @levine10 for related_work)._These can_be_used_to provide a more accurate_estimate of the BH accretion_rate given_conditions at larger radii in a galaxy (Hopkins_&_Quataert, in prep)._Another important improvement would be to use a radiative_transfer calculation to self-consistently determine the_infrared radiation field_produced_by_a central AGN (and_distributed star formation). This could then_be used to calculate the radiation_pressure force on surrounding gas, eliminating the_need for our parameterization of the_force in terms of the_optical depth_$\tau$.

Acknowledgments {#acknowledgments.unnumbered}
===============

We thank Phil Hopkins_and Yuval Birnboim_for useful_conversations. JD and_EQ were supported in part by_NASA grant NNG06GI68G_and the David and Lucile Packard_Foundation._Support for EQ_was_also_provided in_part by the_Miller_Institute for_Basic_Research in Science, University of California_Berkeley._This research used resources of the National_Energy Research Scientific Computing_Center,_which is supported by_the Office of Science of_the U.S. Department of Energy under_Contract No._DE-AC02-05CH11231. This_work was partially supported by the National Center for Supercomputing Applications_under AST080048 and utilized the Intel_64 cluster Abe. The_authors acknowledge_the_Texas Advanced Computing_Center_(TACC) at_The University of Texas at Austin for_providing HPC_resources that have contributed to the_research results reported within_this_paper.

Resolution Studies {#resolution}
==================

![image](FigureA1.eps){width="180mm"}

In this section, we_describe some of our resolution tests_both with and without BH_feedback._In_the absence of feedback, the_well-posed questions for resolution studies include_both how the_gas properties as a function of radius_and_time depend on the resolution and_how_integrated properties of the galaxy (e.g.,_the_star_formation rate) depend on resolution._However, the feedback, when present, has_a nontrivial dependence on the resolution and it is_by no means_clear that the nonlinear system_will_in_fact converge in a simple way with increasing resolution. Physically,_e.g., the_AGN’s radiation pressure_has the strongest effect on the gas that contributes the_most to the optical depth,