P^s_2/P^s_0$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.23\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.55\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=3.00\,{\rm Mpc}^{-1}$]{}]{} \[t\]\[b\]\[1\]\[0\][[IR]{}]{} \[t\]\[b\]\[1\]\[0\][[HR]{}]{}

We have simulated the distribution, in both real and redshift space, for $\xb$ values starting from $1.0$ to $0.1$ with an interval of $\Delta \xb = 0.1$. For each value of $\xb$ we have used 12 independent [realizations [^9]]{} of the simulation to estimate the mean and $1\sigma$ error of the power spectrum. The simulations are all carried out at a fixed redshift $z = 8$. Ideally, one should adopt a model for the evolution of $\xb$ with $z$ and simulate each neutral fraction at the appropriate redshift. This, however, makes the results dependent on the model for the redshift evolution of $\xb$ and also makes the computations rather cumbersome. We have sidestepped this issue by keeping $z$ fixed at a value $z = 8$ which is in the redshift range that is considered favourable for observing the EoR 21-cm signal. However due to this approach we are unable account for the evolution for the dark matter density field as well the evolution of the sources of ionizing photons in our simulation. We interpret the simulations for different values of $\xb$ as representing different stages of reionization.

Figure \[fig:h1\_map\] shows the distributions at different stages of reionization in both real and redshift space for a single realization of the simulation. We see that the distributions in the HR and the IR models are nearly indistinguishable in the early stages of reionization ($\xb \geq 0.7$). In the early stages reionization proceeds inside-out in both of these models with the high density regions ionizing first and the low density regions ionizing later. The maps of the two models become quite distinct at the later stages of reionization ($\xb \leq 0.5$). In the IR model we find many small isolated neutral clumps distributed within the regions which are completely ionized in the HR model. These small clumps correspond to the high density regions
P^s_2 / P^s_0 $ ] { } ] { } \[c\]\[r\]\[1\]\[0\][[$k=0.23\,{\rm Mpc}^{-1}$ ] { } ] { } \[c\]\[r\]\[1\]\[0\][[$k=0.55\,{\rm Mpc}^{-1}$ ] { } ] { } \[c\]\[r\]\[1\]\[0\][[$k=3.00\,{\rm Mpc}^{-1}$ ] { } ] { } \[t\]\[b\]\[1\]\[0\][[IR ] { } ] { } \[t\]\[b\]\[1\]\[0\][[HR ] { } ] { } 

 We have simulated the distribution, in both real and redshift space, for $ \xb$ value start from $ 1.0 $ to $ 0.1 $ with an interval of $ \Delta \xb = 0.1$. For each value of $ \xb$ we have used 12 independent [ realization [ ^9 ] ] { } of the simulation to estimate the beggarly and $ 1\sigma$ erroneousness of the power spectrum. The simulation are all carried out at a fixed red shift $ z = 8$. Ideally, one should adopt a model for the evolution of $ \xb$ with $ z$ and imitate each neutral fraction at the appropriate redshift. This, however, makes the results dependent on the exemplar for the redshift evolution of $ \xb$ and also name the computations rather cumbersome. We have sidestepped this topic by keeping $ z$ fixed at a value $ z = 8 $ which is in the redshift range that is study favourable for observing the EoR 21 - cm signal. However due to this approach we are unable score for the development for the dark matter density field as well the evolution of the sources of ionizing photon in our simulation. We translate the simulations for different values of $ \xb$ as represent different stages of reionization. 

 Figure \[fig: h1\_map\ ] shows the distribution at different stages of reionization in both real and redshift quad for a single realization of the simulation. We see that the distributions in the HR and the IR models are closely indistinguishable in the early stages of reionization ($ \xb \geq 0.7 $). In the early stages reionization proceeds inside - out in both of these model with the eminent density regions ionizing foremost and the low density regions ionizing later. The maps of the two models become quite distinct at the later stage of reionization ($ \xb \leq 0.5 $). In the IR model we find oneself many minor isolated impersonal clumps distribute within the regions which are wholly ionized in the HR model. These small clumps match to the high concentration regions
P^s_2/P^d_0$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.23\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.55\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$y=3.00\,{\rm Mpc}^{-1}$]{}]{} \[t\]\[b\]\[1\]\[0\][[IR]{}]{} \[t\]\[b\]\[1\]\[0\][[KE]{}]{}

We hate simumated thd distribution, in both real end eedshuft space, for $\xb$ valuer startinh from $1.0$ ro $0.1$ xith an interval of $\Delta \xb = 0.1$. Nor ecci value of $\xb$ wg have used 12 independent [raauivations [^9]]{} of the simulation to estimaee the kewn and $1\sigma$ etror ps ths power spectrum. The simulations ade all barried out at a gixed redshift $z = 8$. Ideally, one should adopt a mofel for the evojytion of $\xb$ dith $z$ and simulate eadh neutral fraction at the apprupriace redshift. Tyis, vowever, makxs the results depcmdent mn the kodel for the vedshmft wvolution of $\xb$ and anso makes the compotations rdtker cumbersome. We havw widesjeppeg thkw irsut bb kseping $z$ hixed at a balue $z = 8$ wyich is in the redsnisn range that js confiqered favourable for observing the EoR 21-bm sjgnal. However due to thus approach we are unwble accotnt for the evolution for the dark matter density fielv xs ccju tje evolution of the sources of ionizing photoga on our simulatiok. We interpret the slmijations for dkffereuf balues of $\xb$ as reoresentyng dufferent ftagrs of reionization.

Figure \[fit:h1\_map\] shows nhe eistributions at dnfferent stayes of reiomization in both real aud redahift space for a sihele realization uf nhe vimulation. We see that the distribuvions in the HR snd thq IR modeld are nearly indistinguidhablg in tve early shages of reionization ($\xb \geq 0.7$). Ii the early sjagas geionizatnon prpceeds inside-jut in both of these iodelr with the high dxnsity regiogs ionizing fhtst and the lmw densiey rwgiobs ionixkng later. The kaps of tkt two modeos become quite dixtivdt at the later stqges of reionizstiun ($\vb \lxq 0.5$). Ig the IR moden we fivc manh small isooxted neutral clumps distsibufed within the regooks which qre compjetely ionizec in the HR model. Uhese vmaul climpf correspond to the high densify regiond
P^s_2/P^s_0$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.23\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.55\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=3.00\,{\rm Mpc}^{-1}$]{}]{} We simulated the in both real values from $1.0$ to with an interval $\Delta \xb = 0.1$. For each of $\xb$ we have used 12 independent [realizations [^9]]{} of the simulation to the mean and $1\sigma$ error of the power spectrum. The simulations are all out a redshift = 8$. Ideally, one should adopt a model for the evolution of $\xb$ with $z$ and each neutral fraction at the appropriate redshift. This, makes the results dependent the model for the redshift of and also the rather We have sidestepped issue by keeping $z$ fixed at a value $z = 8$ which is in the redshift range is considered observing the 21-cm However to this approach unable account for the evolution for density field as well the evolution of the of ionizing in our simulation. We interpret the for different values of $\xb$ as representing different of reionization. Figure \[fig:h1\_map\] shows the distributions at different stages of reionization in both real space for a single of the simulation. see the in HR and IR models are nearly indistinguishable in the early stages of reionization \geq 0.7$). In the early stages reionization proceeds inside-out in these with the high regions ionizing first and low regions ionizing later. The the models at later of reionization ($\xb \leq In the IR model we many small isolated neutral which are completely ionized in the HR model. small clumps correspond to the high density
P^s_2/P^s_0$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.23\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.55\,{\rm Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=3.00\,{\rm mpc}^{-1}$]{}]{} \[t\]\[b\]\[1\]\[0\][[IR]{}]{} \[t\]\[b\]\[1\]\[0\][[Hr]{}]{}

We haVe sImuLaTed tHe diStribution, in boTH reaL and redshift space, for $\xb$ ValueS sTArtiNG fRom $1.0$ to $0.1$ With an iNTeRVAl oF $\DElTa \xB = 0.1$. FOR eAch vaLue Of $\xb$ we hAve used 12 indEpeNdEnt [realizatiONs [^9]]{} Of the simulAtiOn to estimate The Mean anD $1\sIgmA$ Error Of tHe powEr specTRum. The SimulatioNs ARe all cARried ouT AT a FixeD redshift $z = 8$. Ideally, ONe SHould adopt a modEl for tHe EVoLUTioN of $\Xb$ with $z$ and SiMulatE Each neuTRaL FRActIOn at the appropRiate redshiFT. ThIs, howeVeR, maKEs the rEsultS dEPenDent on the moDel fOr the redsHift evOLution oF $\Xb$ and alSo makeS thE coMputATiOnS raThER cuMBeRsoME. We Have sideStEpPed thIs isSUE BY keePinG $z$ fiXed at A value $z = 8$ which iS in The rEDshIft raNge thAt is CoNsideRed favOurabLe For observing the eoR 21-cM signal. HoWevEr Due To This aPProach We aRe uNable acCount foR The EvOLUTiOn for the dark matter DeNSItY field as Well thE EvOlUTion of thE sOurCes oF IOniziNg phOToNs in our sImulatIOn. we InterprEt The simUlAtiOns For diFFereNt valuEs of $\xb$ as RepreSEnting differenT Stages of reionIZaTIOn.

fIgurE \[fiG:h1\_map\] shows tHe diSTribUtioNS aT diFFerenT stagEs OF rEIonization in both reaL aNd redsHift sPace for a singlE realizatiON OF the simuLatiON. WE See that the distRibutIons in the Hr And the IR ModelS are nearLy indistiNGUishable In tHe eArlY stAGEs Of reionizatioN ($\XB \geq 0.7$). in The earlY stAges reiOniZatIon ProCeEds inside-Out in botH oF tHeSe ModEls wiTH the high DeNsiTy RegIons iONizing First And tHe LoW DenSity regIOnS IOnizInG lAter. the MaPs of tHe twO ModEls becoMe quite diStiNCt at ThE lAter staGes of reionizaTiOn ($\xb \leq 0.5$). In tHe iR mOdel we FINd many smAll isolated neutral clumpS DistribUteD withIn thE regions wHicH are coMplETely ioNized iN the Hr mOdeL. tHese sMALl CluMpS corresponD TO thE high DeNsitY regionS
P^s_2/P^s_0$]{}]{} \[c\]\[ r\]\[1\]\[ 0\][[ $k= 0.2 3\ ,{\r m Mp c}^{-1}$]{}]{} \[c\ ]\[r\]\[1\]\[0\][[$k=0 .55\, {\ r m Mp c }^ {-1}$ ]{}]{}\ [c \ ] \[r \] \[ 1\] \[ 0 \] [[$k= 3.0 0\,{\rm  Mpc}^{-1} $]{ }] {} \[t\]\[b\ ] \[ 1\]\[0\][[ IR] {}]{} \[t\]\ [b\ ]\[1\] \[ 0\] [ [HR]{ }]{ }

We  haves imulat ed the di st r ibutio n , in bo t h  r ealand redshift spac e ,f or $\xb$ value s star ti n gf r om$1. 0$ to $0.1 $witha n inter v al o f  $\ D elta \xb = 0. 1$. For eac h  va lue of  $ \xb $  we ha ve us ed 12independent  [re alization s [^9] ] {} of t h e simul ationtoest imat e  t he  me an and $1 \si g ma$  error o fth e pow er s p e c t rum.  Th e si mulat ions are allcar ried out  at a  fixe d re ds hift$z = 8 $. Id ea lly, one should  ado pt a mode l f or  th eevolu t ion of  $\ xb$  with $ z$ ands imu la t e ea ch neutral fractio na t  t he appro priate re ds h ift. Thi s,  ho weve r ,  make s th e  r esults d epende n ton  the mo de l forth e r eds hifte volu tion o f $\xb$and a l so makes the c o mputations ra t he r cu m bers ome . We have s ides t eppe d th i siss u e bykeepi ng $z $  fixed at a value $ z= 8$ w hichis in the red shift rang e t hat is c onsi d er e d favourable f or ob serving th e  EoR 21- cm si gnal. Ho wever due t o this a ppr oac h w e a r e  u nable account f or t he  evolut ion  for th e d ark  ma tte rdensity f ield aswe ll  t he  ev oluti o n of the  s our ce s o f ion i zing p hoton s in  o ur sim ulation .  W e inte rp re t th e s im ulati onsf ordiffere nt values  of $\xb $as  repres enting differ en t stages o frei onizat i o n.

Figu re \[fig:h1\_map\] show s  the di str ibuti onsat differ ent  stage s o f  reion izatio n inbo thr e al an d re dsh if t space fo r a s ingle  r eali zationof the simulation. Wesee that thedis trib u t io nsi nt heHR and t he IR models ar e nearly i nd i st inguishabl e  in  t he earl y stage s ofr eioniza tion ($\x b \geq 0. 7$ ). I n the  early sta ges reio nizationp rocee d sinsid e-o ut inbo thof th ese mo d els  with  the h ig h dens ity r eg ions ion izing first and the low  densi ty re gio ns ionizi ngl ate r. The ma ps o f the twomod els  beco meq uitedist i nc t a t  thelate r  stages o f  r eio n i za tion ($\xb\ l e q 0 .5$).  In the IR  mod el we find many s m all isolated n eutr a l  cl ump s  dis tr ibuted withinthe  r e g ions whi ch  are comple tely ion iz e d inthe HR  model . These s ma l l clum ps c orr espond to  th eh igh den si ty region s
P^s_2/P^s_0$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.23\,{\rm_Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=0.55\,{\rm_Mpc}^{-1}$]{}]{} \[c\]\[r\]\[1\]\[0\][[$k=3.00\,{\rm Mpc}^{-1}$]{}]{} \[t\]\[b\]\[1\]\[0\][[IR]{}]{}_\[t\]\[b\]\[1\]\[0\][[HR]{}]{}

We have_simulated_the distribution,_in_both real and_redshift space, for_$\xb$ values starting from_$1.0$ to $0.1$_with_an interval of $\Delta \xb = 0.1$. For each value of $\xb$ we have_used_12 independent_[realizations_[^9]]{}_of the simulation to estimate_the mean and $1\sigma$ error_of the_power spectrum. The simulations are all carried out_at_a fixed redshift_$z = 8$. Ideally, one should adopt a model_for the evolution of $\xb$ with_$z$ and simulate_each_neutral_fraction at the appropriate_redshift. This, however, makes the results_dependent on the model for the_redshift evolution of $\xb$ and also makes_the computations rather cumbersome. We have_sidestepped this issue by keeping_$z$ fixed_at a value $z =_8$ which is_in the_redshift range that_is considered favourable for observing the_EoR 21-cm signal._However due to this approach we_are_unable account for_the_evolution_for the_dark matter density_field_as well_the_evolution of the sources of ionizing_photons_in our simulation. We interpret the simulations_for different values of_$\xb$_as representing different stages_of reionization.

Figure \[fig:h1\_map\] shows the_distributions at different stages of reionization_in both_real and_redshift space for a single realization of the simulation. We see_that the distributions in the HR_and the IR models_are nearly_indistinguishable_in the early_stages_of reionization_($\xb \geq 0.7$). In the early stages_reionization proceeds_inside-out in both of these models_with the high density_regions_ionizing first and the low density_regions ionizing later. The maps of_the two models become quite_distinct_at_the later stages of reionization_($\xb \leq 0.5$). In the IR_model we find_many small isolated neutral clumps distributed within_the_regions which are completely ionized in_the_HR model. These small clumps correspond_to_the_high density regions