----------------------------

We use a two-step algorithm for radiative transfer. First, we evolve all simulation volumes solving the radiative transfer equation in a low angular resolution mode simultaneously with the equations of hydrodynamics. To correct for low resolution, we apply a high angular resolution transfer filter iteratively to our solutions at $z=3$ to improve on the equilibrium position of HI, HeI and HeII ionization fronts.

We developed a UVB radiative transfer module for the parallel AMR cosmology N-body Eulerian hydrodynamics code *Enzo*. The fluid flow equations are solved with the PPM scheme [@colella.84] on a comoving grid (see for @oshea......04 for details), and chemical abundances for various ionization and molecular states of hydrogen and helium are solved with a chemical reaction network with 9 species and 28 reactions [@anninos...97; @abel...97]. Each level of grid refinement has twice the spatial resolution of the previous level. We solve the radiative transfer equation with a photon conserving scheme self-consistently at each level of resolution carrying fluxes explicitly from parent grids to all of its subgrids. On each grid, radiative transfer is computed with a timestep which is usually significantly smaller than the hydro timestep and is adjusted adaptively to obtain the best balance between accuracy and the speed of calculations. Due to complexity of the setup we have not included point source radiation in our models yet, and also we limit transport of background radiation to simple sweeps along the xyz-coordinate axis. We then post-process *Enzo* output with much higher angular resolution transfer. Below we describe these two steps in detail.

Stage one: time stepping and parallelization of the base grid
-------------------------------------------------------------

By definition, radiation propagates at the speed of light. In an explicit advection numerical scheme the Courant condition would necessarily require prohibitively small timesteps to guarantee stability. However, our photon conservation technique is inherently stable, and from this standpoint there is no need to take very small timesteps, but accuracy of the solution is an entirely different issue. In many cases where the supply of ionizing photons greatly exceeds the recombination rate I-fronts can propagate at or close to the speed of light, and even in this regime our photon conservation will give an accurate solution as long as the radiation-chemistry timestep is small enough. How do we ensure that we provide timesteps small enough to balance accuracy and speed of calculations?

Let us
---------------------------- 

 We use a two - step algorithm for radiative transfer. First, we develop all model volumes solving the radiative transportation equation in a broken angular resolution manner simultaneously with the equations of hydrodynamics. To correct for abject resolution, we apply a gamey angular resolution transfer filter iteratively to our solution at $ z=3 $ to improve on the equilibrium position of HI, HeI and HeII ionization fronts. 

 We developed a UVB radiative transfer faculty for the parallel AMR cosmology N - body Eulerian hydrodynamics code * Enzo *. The fluid flow equality are solved with the PPM scheme [ @colella.84 ] on a comoving power system (see for @oshea...... 04 for detail), and chemical abundances for various ionization and molecular states of hydrogen and helium are solved with a chemical reaction net with 9 species and 28 reactions [ @anninos... 97; @abel... 97 ]. Each level of grid refinement has twice the spatial settlement of the previous level. We solve the radiative transfer equation with a photon conserving scheme self - consistently at each level of resolution carry fluxes explicitly from parent grids to all of its subgrids. On each power system, radiative transportation is computed with a timestep which is usually significantly smaller than the hydro timestep and is align adaptively to receive the best balance between accuracy and the focal ratio of calculations. Due to complexity of the setup we have not included point source radiation in our model yet, and also we limit transport of background radiation to elementary sweeps along the xyz - coordinate axis. We then post - process * Enzo * output with much high angular resolution transfer. Below we identify these two steps in detail. 

 Stage one: time stepping and parallelization of the nucleotide grid 
 ------------------------------------------------------------- 

 By definition, radiotherapy propagates at the speed of light. In an denotative advection numerical scheme the Courant condition would necessarily require prohibitively small timesteps to undertake stability. However, our photon conservation technique is inherently static, and from this standpoint there is no need to take very small timesteps, but accuracy of the solution is an wholly different issue. In many cases where the provision of ionizing photons greatly exceeds the recombination rate I - fronts can circulate at or near to the speed of light, and even in this regimen our photon conservation will give an accurate solution as long as the radiotherapy - chemistry timestep is small enough. How do we ensure that we provide timesteps humble enough to balance accuracy and speed of calculations? 

 Let us
----------------------------

We kse a two-step algorithm nor radiative trcbsfer. Hirst, ws evolve all simulation volumes solvmng rhe rqdiative transfer equagion in a low angylar eesolution mode simultansluslv xith the equatipns of hydsodynamics. To wofrzct for low resolution, we apply a hidh angukag resolution ttansftr siltsg lteratively to our solutions at $a=3$ to imkrove on the equilobrium position of HI, HeI wnd JeII ionization frlnts.

We deveoopeq a UVB radiagive transfer module fkr the parallel AMR cosmology N-cody Zulerian hyeridyjdmics code *Xnzo*. Tre fluid flow equatimns are solved with tme PPK sxheme [@colella.84] on a cokoving grid (see fot @oshea......04 fos betails), and chemical qbyndanwes xor xqriuus ipnjzatioj aid moleculad states of hydrogen and heliuk wgr solved witg a chqmycal reaction network with 9 species and 28 rsactions [@anninos...97; @abel...97]. Eqch level of grid reflnement hws twice the spatial resolution of the previous lavel. Xe soovc thd rwdiative transfer equation with a photon consqdvonn scheme self-confistently ay favr level of rerolution darrying fluxes exolicitli from parent gwids to all of its subgrids. On wach grid, rabiarive transfer is clmputed witk a tikestel which is usually signnficanfly smaller than the fydro timestep avd ps agjusted adaptively to obtayn the bewt bclance bdtwegn accuwacy and tje spcad of calculations. Due jo com[lexity of the setup we have not included point source tadhatpon in oux modeks yet, and alfo we limit trcnsport jf bazkground rzdiatioi to simple fweeps along dje xyz-coordiiate axis. We rhen post-pruzess *Enzo* outpit with mlck higher qngular resolution trxhsfer. Below we btsctube these two xteos yn dxtail.

Fdage one: tima stdppkmg ana parallelidatkon pf the base grid
-------------------------------------------------------------

By dafinjtion, radiation prppwgates aj the speqd of light. Im an explicit adveftion nuoericsl fcheme the Courant condition wkuld necedsavily require [rohlbityvely smalk timesteps to guarantee stability. Howeter, our photon conservarion technique is iuhtrently stablx, and srom this standpoint there is no need to take yery small timesteps, bht accgracy of the solution is an entirely different issue. In many cases where the sypply of ionizing lhotpns ggeacly exceeqs tie recombination rste I-fronts can propagate at or close tm che speed of light, and even im ghis regime oor photon conservation wilm give am accurate solution as long as yhe radiation-chemistry tlmestep iw smaml enolfh. How do we ensure yhat wz provude timwsteis small enougf tp nalance acfurary cnd speed of calculations?

Let us
---------------------------- We use a two-step algorithm for First, evolve all volumes solving the low resolution mode simultaneously the equations of To correct for low resolution, we a high angular resolution transfer filter iteratively to our solutions at $z=3$ to on the equilibrium position of HI, HeI and HeII ionization fronts. We developed UVB transfer for parallel AMR cosmology N-body Eulerian hydrodynamics code *Enzo*. The fluid flow equations are solved with the scheme [@colella.84] on a comoving grid (see for for details), and chemical for various ionization and molecular of and helium solved a reaction network with species and 28 reactions [@anninos...97; @abel...97]. Each level of grid refinement has twice the spatial resolution of previous level. the radiative equation a conserving scheme self-consistently level of resolution carrying fluxes explicitly to all of its subgrids. On each grid, transfer is with a timestep which is usually smaller than the hydro timestep and is adjusted to obtain the best balance between accuracy and the speed of calculations. Due to complexity setup we have not point source radiation our yet, also limit transport background radiation to simple sweeps along the xyz-coordinate axis. We then *Enzo* output with much higher angular resolution transfer. Below we two in detail. Stage time stepping and parallelization the grid ------------------------------------------------------------- By definition, at speed an advection scheme the Courant condition necessarily require prohibitively small timesteps guarantee stability. However, our stable, and from this standpoint there is no to take very small timesteps, but accuracy the solution is an entirely different issue. In many cases where the of ionizing exceeds the recombination rate I-fronts can propagate at close to the speed light, and even in this regime our photon conservation give accurate solution long as the timestep is small How do we we provide small to speed of calculations? Let us
----------------------------

We use a two-step algorithm for Radiative tRansfEr. FIrsT, wE evoLve aLl simulation voLUmes Solving the radiative traNsfer EqUAtioN In A low aNgular rESoLUTioN mOdE siMuLTaNeousLy wIth the eQuations of HydRoDynamics. To coRReCt for low reSolUtion, we apply A hiGh anguLaR reSOlutiOn tRansfEr filtER iteraTively to oUr SOlutioNS at $z=3$ to iMPRoVe on The equilibrium posITiON of HI, HeI and HeIi ionizAtIOn FROntS.

We Developed a uVb radiATive traNSfER MOduLE for the parallEl AMR cosmolOGy N-Body EuLeRiaN HydrodYnamiCs COde *enzo*. The fluiD floW equationS are soLVed with THe PPM scHeme [@coLelLa.84] oN a coMOvInG grId (SEe fOR @oSheA......04 For Details), aNd ChEmicaL abuNDANCes fOr vArioUs ionIzation and molEcuLar sTAteS of hyDrogeN and HeLium aRe solvEd witH a Chemical reactioN netWork with 9 sPecIeS anD 28 rEactiONs [@anniNos...97; @AbeL...97]. Each leVel of grID reFiNEMEnT has twice the spatiaL rESOlUtion of tHe prevIOuS lEVel. We solVe The RadiATIve trAnsfER eQuation wIth a phOToN cOnserviNg Scheme SeLf-cOnsIstenTLy at Each leVel of resOlutiON carrying fluxeS Explicitly froM PaRENt GRids To aLl of its subgRids. oN eacH griD, RaDiaTIve trAnsfeR iS CoMPuted with a timestep wHiCh is usUally Significantly Smaller thaN THE hydro tiMestEP aND is adjusted adaPtiveLy to obtain THe best baLance Between aCcuracy anD THe speed oF caLcuLatIonS. dUe To complexity oF THe seTuP we have Not IncludeD poInt SouRce RaDiation in Our modelS yEt, AnD aLso We limIT transpoRt Of bAcKgrOund rADiatioN to siMple SwEePS alOng the xYZ-cOORdinAtE aXis. WE thEn Post-pRoceSS *EnZo* outpuT with much HigHEr anGuLaR resoluTion transfer. BElOw we descriBe TheSe two sTEPs in detaIl.

Stage one: time stepping aND parallEliZatioN of tHe base griD
-------------------------------------------------------------

By DefiniTioN, RadiatIon proPagatEs At tHE Speed OF LiGht. in An explicit ADVecTion nUmEricAl schemE the Courant conditiON woUld necessarilY reQuirE PRoHibITiVEly SmALl tIMEsteps to guarantEe stabilitY. HOWeVer, our photON coNsErvatioN techniQue is INherentLy stable, aNd from thiS sTandPOInt There is no nEed to takE very smalL TimesTEpS, but aCcuRacy of ThE soLutioN is an eNTirEly diFferenT iSsue. In Many cAsEs where tHe supply of ionizing photoNs greaTly exCeeDs the recoMbiNAtiOn rate I-frOnts Can propagaTe aT or Close To tHE speeD of lIGhT, anD Even iN thiS Regime our PHoTon CONsErvation wilL GIVe aN accuRatE SolutiOn as Long as the radiatioN-Chemistry timesTep iS SMalL enOUgh. HOw Do we ensure that We pRoVIDe timestEpS small enougH to balanCe ACcuraCy and sPeed of CalculaTIOnS?

let us
-------------------------- --

We use  a tw o-s tep  a lgor ithm  for radiative tran sfer. First, we evolve  allsi m ulat i on  volu mes sol v in g the  r ad iat iv e  t ransf erequatio n in a low  an gu lar resoluti o nmode simul tan eously withthe  equat io nso f hyd rod ynami cs. To correc t for low  r e soluti o n, we a p p ly  a h igh angular resol u ti o n transfer fil ter it er a ti v e lytoour soluti on s at$ z=3$ to im p r o veo n the equilib rium positi o n o f HI,He I a n d HeII  ioni za t ion  fronts.

W e de veloped a  UVB r a diative transfe r modu lefor  the pa ra lle lA MRc os mol o gyN-body E ul er ian h ydro d y n a mics  co de * Enzo* . The fluid f low  equ a tio ns ar e sol vedwi th th e PPMschem e[@colella.84] o n acomovinggri d(se efor @ o shea.. ... .04  for de tails), and  c h e m ic al abundances forva r i ou s ioniza tion a n dmo l ecular s ta tes  ofh y droge n an d  h elium ar e solv e dwi th a ch em ical r ea cti onnetwo r k wi th 9 s pecies a nd 28 reactions [@an n inos...97; @a b el . . .9 7 ]. E ach  level of g ridr efin emen t  h ast wicethe s pa t ia l  resolution of thepr evious  leve l. We solve t he radiati v e transfer  equ a ti o n with a photo n con serving sc h eme self -cons istently  at eachl e vel of r eso lut ion  ca r r yi ng fluxes exp l i citl yfrom pa ren t grids  to  al l o f i ts  subgrids . On eac hgr id ,rad iativ e  transfe risco mpu ted w i th a t imest ep w hi ch isusually si g n ific an tl y sm all er  than  the hyd ro time step andisa djus te dadaptiv ely to obtain  t he best ba la nce  betwe e n  accurac y and the speed of calc u lations . D ue to  com plexity o f t he set upw e have  not i nclud ed  po i n t sou r c erad ia tion in ou r mod els y et , an d alsowe limit transport ofbackground ra dia tion t osim p le swe ep s  al o n g the xyz-coord inate axis .W ethen post- p roc es s *Enzo * outpu t wit h  much h igher ang ular reso lu tion t ran sfer. Belo w we des cribe the s e two st eps i n d etail. 
Sta ge on e: tim e  st eppin g andpa rallel izati on  of thebase grid
------------- ------ ----- --- --------- --- - --- --------- ---- -----

Bydef ini tion,  ra d iatio n pr o pa gat e s atthes peed of l i gh t.I n  a n explicita d v ect ion n ume r ical s chem e the Courant con d ition would ne cess a r ily  re q uire  p rohibitively s mal lt i mestepsto  guaranteestabilit y. Howev er, ou r phot on cons e r va t ion te chni que  is inher ent ly stable,  a nd from t hisst andpoi nt the r e is n o need to take v ery s m a ll ti m est eps,bu t accur a cy o f the solu tion is anentire ly d iffer ent iss ue . In m any  c ases where the suppl y ofionizin gphot ons  great ly e x c eedsthere com binationr a te I- fr o nts  can  prop ag ateat or clo s e to the  sp e ed of l ig ht, a nd eve n  i n this regim e o ur ph o t on conserv a tion wi l l giv e an a ccurat e solut i onas  long a s t h e  radiatio n-chemist r y t im este p is sma ll  eno ug h.Ho w do we ensurethat weprovidet i mestep s s m alle noughto  balance acc uracy  a n d speed o f c alculat i ons ?
Let  us
----------------------------

We use_a two-step_algorithm for radiative transfer._First, we_evolve_all simulation_volumes_solving the radiative_transfer equation in_a low angular resolution_mode simultaneously with_the_equations of hydrodynamics. To correct for low resolution, we apply a high angular resolution_transfer_filter iteratively_to_our_solutions at $z=3$ to improve_on the equilibrium position of_HI, HeI_and HeII ionization fronts.

We developed a UVB radiative_transfer_module for the_parallel AMR cosmology N-body Eulerian hydrodynamics code *Enzo*. The_fluid flow equations are solved with_the PPM scheme_[@colella.84]_on_a comoving grid (see_for @oshea......04 for details), and chemical_abundances for various ionization and molecular_states of hydrogen and helium are solved_with a chemical reaction network with_9 species and 28 reactions_[@anninos...97; @abel...97]._Each level of grid refinement_has twice the_spatial resolution_of the previous_level. We solve the radiative transfer_equation with a_photon conserving scheme self-consistently at each_level_of resolution carrying_fluxes_explicitly_from parent_grids to all_of_its subgrids._On_each grid, radiative transfer is computed_with_a timestep which is usually significantly smaller_than the hydro timestep_and_is adjusted adaptively to_obtain the best balance between_accuracy and the speed of calculations._Due to_complexity of_the setup we have not included point source radiation in our_models yet, and also we limit_transport of background radiation_to simple_sweeps_along the xyz-coordinate_axis._We then_post-process *Enzo* output with much higher angular_resolution transfer._Below we describe these two steps_in detail.

Stage one: time_stepping_and parallelization of the base grid
-------------------------------------------------------------

By_definition, radiation propagates at the speed_of light. In an explicit_advection_numerical_scheme the Courant condition would_necessarily require prohibitively small timesteps to_guarantee stability. However,_our photon conservation technique is inherently stable,_and_from this standpoint there is no_need_to take very small timesteps, but_accuracy_of_the solution is an entirely_different issue. In many cases where_the supply of ionizing photons greatly exceeds the recombination_rate I-fronts can_propagate at or close to_the_speed_of light, and even in this regime our photon conservation_will give_an accurate solution_as long as the radiation-chemistry timestep is small enough. How_do we ensure that we provide timesteps_small enough to balance accuracy and speed of calculations?

Let us