cosmic reionization corresponds to the transition from a fully neutral to a highly ionized intergalactic medium ( igm ) , driven by ultraviolet ( uv ) radiation from the first stars and black holes .
the transition is a key milestone in cosmic structure formation , marking the formation of the first luminous objects .
reionization represents the last major epoch of cosmic evolution left to explore .
study of the igm , galaxies , and quasars present during that time is a primary science driver for essentially all future large - area telescopes , at all wavelengths .
recent observations of the gunn - peterson effect , i.e. , ly@xmath4 absorption by the neutral igm , toward the most distant quasars ( @xmath5 ) , and the large scale polarization of the cmb , corresponding to thompson scattering during reionization , have set the first constraints on the reionization process .
these data suggest significant variance in both space and time , starting perhaps as far back as @xmath6 @xcite and extending to @xmath7 @xcite .
current probes of the reionization are limited : present wmap - v data indicates the 5  @xmath8 detection of the e - mode of polarization which rules out any instantaneous reionization at @xmath9 at 3.5  @xmath8 level . for the gunn - peterson effect ,
the igm becomes optically thick to ly@xmath4 absorption for a neutral fraction as small as @xmath10 .
it has been widely recognized that mapping the red - shifted hi 21 cm line has great potential for direct studies of the neutral igm during reionization @xcite .
there are number of upcoming low - frequency arrays whose key science goal is to detect the hi 21 cm signal from the epoch of reionization ( eor ) .
this includes the murchison widefield array [ mwa ] @xcite , precision array to probe epoch of reionization [ paper ] @xcite and low frequency array [ lofar ] @xcite .
one of the major challenges for all of these upcoming arrays will be the removal of the continuum foreground sources in order to detect the signals from the eor . in this paper
we discuss how the radio interferometric imaging techniques are going to affect the foreground source modeling and subsequent removal from the data - set in order to search for the eor signal .
recently , there has been substantial research on foreground source modeling ( @xcite , @xcite , @xcite , etc ) at these low frequencies .
similar effort has also been made in exploring different techniques to remove the foregrounds from the eor data - set by @xcite , @xcite,@xcite , @xcite , @xcite and @xcite .
they primarily focus on the removal of the sources that are fainter than a certain @xmath11 ( @xmath12  jy ) level .
most of these works do not consider the foreground sources brighter than @xmath11 and how accurately they need to be removed from the data - set by some real - time calibration or modeling technique @xcite in the uv - domain . in this paper
we deal with the bright point sources above @xmath11 and the limitations that will be caused due to imperfect removal of such sources .
one of the objectives of this paper is to demonstrate the effect of the frequency dependent side - lobes of the synthesized beam on the foreground source removal strategies .
the effect of other frequency dependent terms like primary beam , ionosphere , etc will be addressed in future papers .    in section-2
, we discuss briefly our choice of model for the eor signal , i.e. the largest expected cosmic stromgren sphere ( css ) which represents the only signature that can be detected in the image domain by the upcoming radio - telescopes .
section-3 presents the detailed array parameters that have been used in the simulations . in section-4
, we discuss the foreground source model that has been used in every simulation performed for this paper .
section-5 outlines the simulation methodology and a possible data reduction procedure that might be followed while processing the raw - data from the upcoming low - frequency radio telescopes in order to extract the eor signal .
this procedure may not be exactly the same as what will be actually implemented for these upcoming telescopes . in section-6
, we discuss the propagation of different forms of error through the radio interferometric data reduction procedure .
we consider two frequency independent errors : i ) position error in the global sky model ( gsm ) that is used to remove the bright sources above the @xmath11 level and ii ) residual calibration error .
we also present the results showing how these errors propagate to the final residual ( foreground - removed ) spectral image - cube .
finally , in the last section we discuss the implications of the results from our simulations and our recommendations for the upcoming low - frequency arrays in order to detect the signal of cosmic reionization .
first generation low frequency arrays like mwa , paper , lofar , etc are meant to detect three potential signatures of eor : i ) 3d power spectrum , ii ) rare and large cosmic stromgren sphere ( css ) and iii ) hi 21 cm forest @xcite .
all of the above signatures require the standard low - frequency radio - interferometric calibration and imaging like self calibration , frequency dependent calibration and source removal .    for simplicity ,
we restrict ourselves to the large and rare cosmic stromgren sphere which are formed around the luminous quasars at the end of reionization .
these css form the only potential eor signature that can be detected in the image domain by the upcoming radio telescopes .
these css are rare and are expected to have a brightness temperature of 20 ( @xmath13 ) mk with a physical size of @xmath14  mpc physical size , where @xmath13 is the neutral fraction of the inter galactic medium ( igm ) . the physical size of these rare css has been derived from the ly-@xmath4 spectra @xcite . according to @xcite , we derive the angular size of the css ( @xmath15 ) from the co - moving radius of css ( @xmath16 ) : @xmath17 this gives the angular scale of the css to be @xmath18  arc - minutes along with a line - width of @xmath19 given by @xcite : @xmath20 the line - width of the css translates to 2.5 mhz .
thus the total flux density of the css is about 0.24 @xmath13 mjy which gives a surface brightness of : @xmath21 where @xmath22 is the size of the synthesized beam for an array with maximum baseline 1.5 km at 158 mhz ( @xmath23 ) .    for a late reionization model , @xcite predict that the @xmath24 field - of - view and 16 mhz bandwidth of mwa , will include atleast one of these large and rare hii regions ( @xmath25 mpc at @xmath26 ) .
moreover , there is also the possibility of finding smaller hii regions ( @xmath27 mpc ) and up to @xmath28 fossil hii regions due to nonactive agn within the same field - of - view , depending on the duty cycle @xcite .
table 1 outlines the basic array parameters that we have adopted .
we note that most of these parameters reflect actual specifications for the upcoming mwa-512 array .
figure 1 shows the array layout for the 512 element array with maximum baseline of 1.5 km
. this might not be the final array design or specification for mwa .
cc no . of tiles & 512 + central frequency & 158 mhz ( @xmath29 ) + field of view & @xmath30 15@xmath31 at 158 mhz .
( @xmath32 ) + synthesized beam & @xmath30 4.5 at 158 mhz .
( @xmath32 ) + effective area per tile & @xmath30 17 @xmath33 + maximum baseline & @xmath30 1.5 km .
+ total bandwidth & 32 mhz + @xmath34 & @xmath30 250 k + channel width & @xmath30 32 khz + thermal noise & @xmath30 15.4 @xmath35  jy / beam +    & ( @xmath36 hours @xmath37 @xmath38 mhz ) +    in order to detect the signal from cosmic reionization , arrays like mwa , paper , lofar , etc will have to overcome the thermal noise limitation @xcite : @xmath39 where @xmath40 denotes the final rms noise in the image from the channel width of @xmath41 after @xmath42 hours of integration .
@xmath34 , @xmath43 and @xmath44 denotes the system temperature , effective collecting area of each element / tile and number of elements / tiles in the array . according to the above equation ,
the thermal noise is 15.4 @xmath35 jy / beam after @xmath45 hours of integration with @xmath46k at these low frequencies is dominated by the sky temperature ( @xmath47 ) @xcite .
we have adopted the above noise value for 158 mhz . in practice , one should quote the noise figure for the lowest frequency edge which corresponds to the highest value for @xmath48 . ] and channel width of 2.5 mhz , which is also the spectral width of a css , as discussed in the previous section
.    most of the upcoming low frequency telescopes will be transit - instruments and will observe a field around its transit .
hence we have used 6 hours of integrations for all the simulations , assuming that the telescopes will observe a field between @xmath49   3 hours in hour angle . here
, we have assumed that the field - of - view of the instrument to be @xmath24 .    in this paper
, we aim to demonstrate the case for the 512-element array , as mentioned in table 1 .
however , the simulations of 512-element are expensive .
hence , we have considered a simpler geometry with 128-element array ( figure 1 ) , with a maximum baseline of 600 meters .
this array design does not represent any of the upcoming telescopes but has been adopted in order to simplify the simulations .
apart from the effective area and synthesized beam value the rest of the specifications for the 128-element array remains the same as in table 1 . in our test simulations , we have used identical residual calibration errors and created separate resultant spectral image cubes for the two different array specifications ( 128-element and 512-element ) .
the rms noise level in identical regions of these two maps differed by a factor of @xmath50 .
the synthesized beam or the psf ( point spread function ) image from these two array configurations shows that the rms noise level in psf side - lobes for the 512-element is a factor of @xmath50 lower than that in the 128-element array .
this confirms that there is a scaling property between the 512-element and the 128-element array .
the same scaling property is also evident in the results from our test simulations with the position errors .
we have exploited this property and have used the 128-element array for all the simulations referred to in the later sections and have subsequently scaled down the rms noise values by the scaling factor , in order to represent the same for the 512-element array .
we also note here that the thermal noise and the css surface brightness values that are quoted throughout the paper corresponds to the 512-element array configuration .
since we are dealing with point sources in our foreground source model the change in the maximum baseline between the two configurations will not affect our conclusions .
the only frequency dependent component in our simulations is the side - lobes of the point spread function ( psf ) . in this section
we discuss the effects of different weighting schemes on the psf shape and its dependence on frequency .
figures 2 and 3 show the effect of three different weighting functions on the psf side - lobes for the 512-element array .
these figures show different plots from three well - known weighting schemes used in the synthesis imaging : natural , uniform and robust .
@xcite showed that the robust weighting scheme smoothly varies between the uniform to natural weighting schemes .
figure 3 shows how the near - by side - lobes are almost coherent for two different channels separated by 32 mhz while the far - away side - lobes become more stochastic in frequency .
the effect of the different weighting scheme shows similar results for the 128-elements as well but with a higher psf side - lobe level .
the rms value computed on the far away side - lobe level of the psf for the 128-element array is about @xmath50 times higher than that in the 512-element array .
we recognize that the choice of optimal weighting scheme depends on a number of parameters based on the array configurations , telescope hardware specifications and also on precise data reduction algorithms used .
hence all the upcoming eor experiments need an optimal weighting scheme for their individual data reduction procedures .
however , exploring the optimal weighting scheme is beyond the scope of this paper . for our work
, we have adopted natural weighting scheme .
psf side - lobe level reduces with increase in hour angle coverage . in our simulations
we have used six hours of observing time to produce each data - set . however , in the actual experiment the observing time for each data - set may be less than that .
hence the psf side - lobe level will also degrade accordingly .
our foreground sky model only includes point sources .
no extended emission from galactic foreground is included as a part of the sky model .
recently , there has been extensive research regarding the foreground source removal dealing with sources below @xmath51  jy level .
most of these analysis implicitly assume that foreground sources above the @xmath11 level are being removed perfectly .
but in reality any imperfect calibration will introduce artifacts in the residual data after source removal .
hence in this paper our main aim is to explore the level of accuracy needed in these calibration procedures in order to ensure the residual errors from the strong foreground source removal do not obscure the detection of the signal from cosmic reionization .
the choice of @xmath11 level is not totally arbitrary and will be discussed in the final section of this paper .
the sky model is derived from the logn - logs distribution of sources and is termed as the global sky model ( gsm ) from now onwards . since our gsm only includes sources above 1 jy , we follow the source count from the 6c survey at mhz @xcite . the source count from the 6c survey
is given as : @xmath52    for a field - of - view of @xmath24 the total number of sources ( @xmath53  jy ) equates to @xmath54 , following the above power - law distribution .
the entire flux range between 1 - @xmath55  jy was divided into several bins .
the source population in each of these bins has been predicted following the above logn - logs distribution . in order to assign fluxes to individual sources inside a bin ,
a gaussian random number generator has been used .
the strongest source in our gsm is @xmath56  jy . the center of the field is chosen such that it coincides with one of the cold spots in the foreground galactic emission seen from the southern hemisphere .
southern sky is chosen since both the upcoming arrays mwa and paper are being constructed in western australia .
the exact field center used for the gsm is 4 hours in right ascension and -26 degree in declination . in order to assign a position to each of these sources within the field - of - view ,
another gaussian random number generator has been used which predicted the offset from the field center for respective sources . in the gsm
all the foreground sources are flat spectrum , i.e. with zero spectral index ( @xmath57 ) .
figure 4 shows the cleaned image of the gsm used for all of the simulations .
the clean algorithm used in this process is a wide - field variant of the well - known clark - clean algorithm , which uses w - projection algorithm for 3-dimensional imaging .    in practice
, the upcoming eor telescopes should be using a similar global sky model ( gsm ) but created from the existing source catalogs available for the part of the sky and frequency ranges of their observations .
these telescopes will use the gsm as the preliminary sky model in order to detect and subsequently remove the bright foreground sources in the observed data - set .
all the simulations have been performed within the casa package . as mentioned before ,
simulations for both 512-elements as well as 128-elements have been conducted and compared in the later sections .
the @xmath24 field - of - view will include @xmath54 bright sources ( @xmath58  jy ) .
the individual flux densities of these foregrounds are @xmath59 time higher than the signal from cosmic reionization that these instruments are aiming to detect .
so the challenge lies in calibration and subsequent removal of such bright sources from the raw data - sets . in order to add to the constraints
, the data rates of most of these telescopes ( e.g. @xmath60  gb / s for mwa ; @xcite ) will not allow them to store the raw visibilities produced by the correlator .
hence real - time calibration and imaging needs to be done in order to reduce the data volume and store the final product in the form of image cubes @xcite .
the critical steps include removal of the bright sources above the @xmath11 level from the data - sets in these iterative rounds of real - time calibration and imaging procedure . as a result
the residual image - cubes will not be dominated by these bright sources and the rest of the foregrounds can be removed in the image domain .
however , the accuracy of the foreground source removal strategies are strongly dependent on the data reduction procedure . the likely data reduction procedure which will be followed by the upcoming telescopes
can be outlined as :    * the raw data - sets from the correlator will go through real - time calibration and subsequent removal of the bright sources based on some global sky model ( gsm ) , down to @xmath11 level , in the uv domain . *
the residual data - sets will be imaged and stored as a cube for the future processing and removal of sources which are below @xmath11 .
however , due to imperfect calibration during the process of real - time calibration and source removal if any systematic error is incorporated in the residual image - cubes they can not be minimized by simply averaging over a large number of such residual image - cubes . hence to understand the level of accuracy required in the calibration procedure to reduce the residual systematic errors in the data - set is critical .    here
, we outline the reduction pathway that have been followed throughout the paper :    * the observed visibilities ( @xmath61 or @xmath62 ) are simulated for 6 hours ( @xmath63  hours in hour - angle ) from the gsm ( with no position or flux error ) and the array configuration from table : 1 . *
model visibilities ( @xmath64 ) are also generated using the same gsm but they are now corrupted with either gsm position errors or with residual calibration errors . *
model visibilities are subtracted from the observed visibilities to residual visibilities ( @xmath65 ) .
this step will be referred to as uvsub @xcite .
* we fit a third order polynomial in frequency to the dirty image - cube ( @xmath66 ) formed from @xmath67 .
this step will be referred to as imlin @xcite .
after the imlin step , we investigate the rms noise in the residual image cube .
we have tried to use different orders of polynomial for the imlin step but our conclusions did not change significantly .
this may be mainly due to the fact that the only frequency dependent term in these simulations is the psf .
however , it should be mentioned in this context that on using higher order polynomials in the imlin step we will be taking out relevant structures at those scales in the frequency space @xcite .
since the reionization signal is a spectral signature with a width @xmath68  mhz ( size of the largest possible css ) , using a very higher order polynomial in the imlin might remove the signature itself .
we have also explored uvlin @xcite , which fits a polynomial in the uv domain to @xmath67 . but uvlin works perfectly only within a small field - of - view depending on the channel width in frequency @xcite .
hence , in our proposed reduction pathway we are using uvsub followed by imlin only .
in this section we describe the effects of different errors arising from radio interferometric data processing and how they propagate and show up as artifacts in the final residual image . here , we consider effects of the errors arising from the imperfect gsm and residual calibration errors on removal of strong foregrounds .
the source positions in the global sky model ( gsm ) are accurate up to a level , e.g. the nvss catalog has an accuracy of @xmath69  arc - seconds in positions of the sources detected @xcite . the real - time calibration procedure for mwa like telescopes might implement some method to iteratively solve for the positions of bright sources taking the initial value from the gsm , e.g. asp - clean @xcite .
an iterative scheme will also converge with some residual position and flux errors in the model .
hence , any error in the final gsm position or flux will contribute to imperfect source removal .
the resultant artifact will limit the dynamic range in the residual data - set and obscure the detection of the faint signal from reionization .
the main aim of the simulations in this section is to determine the accuracy of the gsm ( position and flux ) required in order to reach the required rms noise level .
the gsm position errors considered in the following simulations are actually systematic residual errors after the iterative position determination . here
, we deal with residual gsm errors in position and will also mention briefly the effect of the gsm flux errors in the residual data - set .    in order to introduce the position error in the gsm we have added an error term in the right ascension ( @xmath4 ) of each of the sources in the gsm .
this error term is derived from a gaussian distribution:- @xmath70 where @xmath8 denotes the error level in the position , i.e. for nvss catalog @xmath71  arc - seconds .
the steps in the simulation performed in order to explore the effect of the gsm error in position is as follows :
@xmath61 are simulated .
2 .   predict model visibilities from the gsm model with source - position error ( @xmath72 ) .
3 .   subtract the model visibilities from the observed visibilities to get @xmath73 .
4 .   make @xmath66 and apply imlin as mentioned befor.e    in order to explain the above mentioned steps in detail we refer to figure 5 , which demonstrates the procedure of the reduction followed in the simulations involving gsm errors .
figure 5a shows the dirty map ( peak value : @xmath74  jy / beam ) from the residual visibilities where the model visibilities where corrupted with a gsm position error of 0.01 arc - second .
figure 5b ( peak value : @xmath75  jy / beam ) is generated after applying imlin to the residual image .
the gray - scale color - bars represent the variation of the intensity in two images .
here we discuss the implications of the results from the simulations with the gsm errors .
figure 6a represents the variation in the rms noise level of three regions in : [ 1 ] true image of @xmath61 ( figure 4 ) , [ 2 ] image of @xmath76 after uvsub ( figure 5a ) and [ 3 ] image after imlin ( figure 5b ) .
the gsm position errors used in this case is 0.01 arc - second .    in figure 6b , the rms noise value in three different parts of the field - of - view have been plotted as a function of different residual gsm position errors .
the rms values quoted here are from the images of residual visibilities after uvsub [ 2 ] .
hence the solid line denotes the trend in the decrease of the rms noise in the uvsub image with the decrease in the gsm position error .
the magnitude of the position error ( value of @xmath8 as in equation 6 ) has been varied from 6 arc - seconds ( as in nvss catalog ) down to an extremely low value of @xmath77  arc - seconds .
the original simulations were performed with 128-element array and the resulting rms values are scaled down by factor of 5 in order to represent the same for the 512-element array . this scaling property has already been discussed in details in section-3 of this paper . according to the figure
6b it is evident that an accuracy of 0.1 arc - second in gsm position is needed to achieve the required rms noise in order to detect the reionization signal after removal of strong foreground point sources .
the solid curve in figure 6b represents the mean curve of the rms noise level variation and has a mean slope of @xmath12 in the log - log space .
this means that a decrease in the rms noise level by another order of magnitude is possible on achieving an extra order of magnitude accuracy in the gsm position estimates .
@xcite have shown similar results in demonstrating the effect of the pixelation error on the dynamic range of the image .
recently , @xcite have shown that a image made from a data - set using the vla ( 27-element ) obtains a dynamic range limitation of @xmath78 when the sources are shifted by 0.01 pixel in the image plane .
this result is similar to the result quoted in figure 6b , where it is evident that gsm position error of 1 arc - second ( @xmath79 pixel ) results in an rms noise of @xmath80 jy / beam corresponding to a dynamic range ( dr ) of @xmath81 .
the better psf side - lobe level in 512-element array over vla27-element array causes an extra order of magnitude gain in the dynamic range in the case of the former .    here ,
we have considered only systematic residual gsm position errors which are the same from one day to the another .
this will not reduce after averaging over residual image - cubes from successive days of observations .
the real - time calibration system might implement an iterative scheme ( e.g. asp - clean , @xcite ) to determine the source positions .
however , unless the systematic error determined by any such scheme is smaller than 0.1 arcsecond , averaging over long lengths of time will not mitigate the errors due to imperfect gsm .
the gsm position error limit that we report in this paper is thus also applicable for the specifications of the accuracy that any such scheme delivers .
in this section we discuss the results of introducing gsm flux error in our simulations .
gsm flux errors have been derived from a gaussian random distribution ( similar to equation 6 ) with mean zero and @xmath82 of flux values of respective sources .
similar to the steps of simulations discussed in section 6.1 , the gsm model with flux error is used to corrupt model visibilities .
the residual image cube shows sources with flux values of @xmath83 of the original sources in the same position .
since the gsm contains only sources between 1 jy and @xmath56  jy , the maximum residual flux density in the residual image is @xmath84 jy / beam and can be denoted as sources below the @xmath11 limit @xcite .
the removal of such sources are beyond the scope of this paper .
flux errors might be as high as @xmath85 , as quoted for most of the source catalogs such as nvss .
the residuals will show up as sources with flux density as high as few jy / beam . in this case
the residuals are well above the level of the @xmath11 and perfect removal of such sources in the subsequent imlin step still remains a challenge .
this section deals with the effect of residual calibration error in the foreground source subtraction from the raw data - set coming out of the upcoming eor instruments .
the antenna dependent complex gains can be defined as : @xmath86 where @xmath87  s are the antenna dependent complex gains , @xmath88 are the corrupted visibilities and @xmath89 are the true visibilities or @xmath62 in our case ; @xmath90 denotes time variation .    the model for the complex gains used in the simulations is given by : @xmath91 where @xmath92 , @xmath93 , @xmath94 and @xmath95 are for amplitudes ( @xmath96 ) and phases ( @xmath97 ) .
the gains in both amplitude and phase are derived from individual gaussian random distribution as in equation 6 .
the process of simulation is as follows :    1 .
@xmath61 are simulated .
model visibilities are predicted from a perfect gsm and are corrupted with @xmath98 to compute @xmath99 3 .   in uvsub
we get : @xmath100 .
4 .   make @xmath101 and apply imlin as mentioned before .
the model for the complex antenna - dependent gain errors ( equation 7 ) includes two terms for amplitude ( @xmath92 ) and phase @xmath93 gains , which are constant for each antenna throughout each day of observing ( 6 hours ) . along with them
there are also two small additive random offsets @xmath102 , which vary within a single day of observation .
following @xcite , it has been determined through several simulations that the constant offsets ( @xmath103 ) are the dominant terms compared to the time dependent terms . moreover , since the time variation for @xmath104 are derived from a gaussian distribution similar to equation 6 , the effect of these terms decreases upon averaging over time .
hence , in our simulations we have only included the constant gain errors @xmath105 and redefined equation 8 as : @xmath106    it is certainly evident that some variant of the self calibration algorithm will be implemented in the real - time calibration procedure for the upcoming mwa like telescopes in order to facilitate the bright foreground source removal . in our simulations , we have not implemented any form of calibration algorithm . instead , we have used the measurement equations with residual calibration errors in order to explore the accuracy with which any calibration procedure must work in order to achieve the desired rms noise limit and detect the signal from reionization .
figure 7 demonstrates the three steps followed to simulate the effect of the residual calibration errors in the eor data - set .
figure 7a shows the dirty image of the residual uvsub data and figure 7b shows the image after imlin .
both these figures are from the simulation with @xmath107 amplitude error or @xmath108  degree phase error .
figure 8a shows the variation in the rms noise level in the three regions of the images made using : [ 1 ] @xmath61 ( figure 4 ) , [ 2 ] @xmath76 in uvsub ( figure 7a ) and [ 3 ] after imlin ( figure 7b ) .
it is clear from figure 8a that the rms noise level reached near the bright sources after the uvsub stage is @xmath10 jy / beam , corresponding to a dynamic range of @xmath109 .
@xcite calculated the dynamic range limitation due to a time independent amplitude error in all the baselines that will be introduced to a point source data . if we denote the visibility function ( @xmath110 ) as : @xmath111 where @xmath112 are the amplitude and phase errors respectively , the dynamic range of the image is limited to : @xmath113 where n is the number of antennas . in our case , @xmath114 , @xmath115 and @xmath116  degree .
these values yield a dynamic range of @xmath109 , which agrees with the dynamic range we obtained from our simulations .
figure 8b denotes the rms noise value in different parts of the uvsubed image as a function of different residual calibration errors .
the three different residual calibration errors considered here are @xmath117 , @xmath107 and @xmath118 in amplitude or 0.01 degree , 0.1 degree and 1 degree in phase respectively .
the resultant curve in the log - scale has a slope of @xmath119 .
this indicates that in order to achieve the desired rms noise level to detect the reionization signal , the residual calibration errors should be @xmath120 in amplitude or @xmath121  degree in phase . in practice , it is extremely challenging to achieve such accuracy with a real - time calibration procedure . here ,
we have only considered systematic residual calibration errors which will not vary from one day to the other .    as before ,
these simulations were originally performed for the 128-element array and then the rms noise values were scaled down by a factor @xmath50 following the scaling property that has been discussed in initial sections .
however , we have also repeated the simulations for the 512-element array with residual calibration error of @xmath122 in amplitude or 0.1  degree in phase and found the results consistent with our scaled values .      in the previous section it
is clearly shown that even the low systematic residual calibration errors restrict the dynamic range in the residual image to @xmath78 , whereas the desired dynamic range is @xmath123 in order to detect the signal from reionization .
the obvious next step to follow is to add the uvsubed image cubes from successive days in order to check whether the residual calibration error reduces as @xmath124 , where @xmath125 denotes the number of days over which the uvsub images has been added and each day consists of 6 hours of observations .
we can demonstrate the same from the expression of the residual visibilities ( @xmath100 ) .
if the calibration procedure is able to remove all systematic gain errors from the corrupted data contributing to residual calibration errors that are purely random from day to day , then they will reduce upon averaging .
figure 9 demonstrates the effect of the averaging over 20 days of uvsubed images with the residual calibration errors of @xmath122 in amplitude or @xmath108  degree in phase .
figure 9a denotes the dirty image after the first day and figure 9b denotes the averaged image after 20 days .
the figures demonstrate a significant reduction in the rms noise level after 20 days of averaging .
figure 10a shows the reduction in the rms noise level in different parts of the image on averaging over @xmath126 days .
the solid line denotes the theoretical noise curve that follows @xmath124 .
hence , the reduction in rms noise on averaging over different days is consistent with the theoretical prediction ( @xmath124 ) .    a uvsubed image from each day ( 6 hours )
is also limited by the thermal noise , which reduces down as @xmath124 .
if the residual calibration errors are random between successive data - sets from different days of integration then it should also reduce down as @xmath124 . as a result
, it should then be sufficient to achieve a resultant rms noise level in the images below the daily thermal noise limit . from table 1
we can estimate that the thermal noise limit from 1 day of integration ( 6 hours ) is @xmath127  jy / beam .
figure 9b is the same as figure 8b except that we have plotted the new noise estimates from 6 hours of observation .
this figure now explains that the residual calibration errors should be @xmath128 in amplitude or @xmath129  degree in phase in order to reach below the daily thermal noise limit .
once the rms noise level in the uvsubed image is below the daily thermal noise limit , it can reduce down along with the thermal noise in order to detect the faint reionization signal within the time limit set by the telescope sensitivity ( i.e. @xmath45  hours ) .
the most important caveat in the entire analysis is the fact that we assume that the residual calibration errors are not systematic from one day to the next .
if there are any systematics left over from day to day ( e.g. due to antenna gain errors , antenna primary beam errors , pointing errors , etc ) then those will not go down as @xmath130 .
moreover , even with these assumptions the much needed accuracy in calibration is alarmingly high .
there are few important results and concerns that have come out from the simulations performed in this paper .
the simulations have been performed for both 128-element array and 512-element array .
necessary comparisons between both arrays have been discussed at essential stages of this paper . the comparison between the two configuration
can be well explained by a simple scaling relation which is the ratio between the rms value of the psf side - lobe levels in two configurations .
hence it was sufficient to perform the 512-element simulations only for the optimal value of the gsm position error and residual calibration error .
the choice of @xmath51  jy in our gsm is not totally arbitrary keeping in mind the fact that the real - time calibration procedure needs to model and remove all the sources down to @xmath11 from the raw uv data - sets within a duration constrained by the instrumental specifications . since the errors we have discussed in this paper are dominated by the bright sources inside the field - of - view , our analysis will not change significantly by reducing the @xmath11 level and including more and more fainter sources .        in the section dealing with the gsm position errors , we have only discussed about systematic errors which do not change from one day to the next .
if the calibration procedure involves some iterative position calibration scheme that will evaluate the source position everyday then the errors from those procedure should vary from one day to the other .
hence , there is a possibility of reaching the desired rms noise level even with lesser accuracy in the gsm position estimates .    in practice
, a better gsm can be obtained a priori from a deepsurvey of the specific fields to be observed with mwa , lofar , etc . using the present facilities like gmrt and wsrt @xcite at 150 mhz .
mwa , lofar , etc can also improve the existing gsm after sufficient days of observations .
the simulations with calibration errors assume that the image cubes are formed after each day of observations ( 6 hours ) .
we assume that the calibration errors are incoherent beyond 6 hours and as a result the effect of calibration error will reduce with successive days of integration .
if calibration error have shorter coherence time , then lesser accuracy in the calibration can be tolerated in order to reach below the thermal noise limit for the shorter integrations .
however , the systematics need to be removed at those timescales to make sure that the residual errors are still random from one data - set to the next .
it should be noted that the psf side - lobe level will also be higher due to shorter integration period , which in turn might make it harder to solve and remove the systematics at a shorter timescales .
the data analysis pipeline adapted in this paper resembles mostly the upcoming mwa telescope .
since the raw visibilities can not be stored , the calibration and removal of the brightest sources in the real - time calibration pipeline is one of the major step in foreground subtraction @xcite .
the accuracy in gsm position and residual calibration ( discussed in this paper ) is critical for extraction of the eor signal . in other telescopes like lofar , where it is possible to store the raw visibilities ,
the calibration will be carried out through a combination of real - time processing and off - line reprocessing scheme .
the ability to store the raw visibilities will allow more accurate calibration of the data in the off - line mode .
however , the combination of real - time and off - line calibration pipeline for lofar still needs to achieve similar accuracy in calibration ( discussed in this paper ) in order to reach the desired dynamic range to extract the eor signal .    in the present paper
we have not discussed in great detail the step of polynomial fitting in the image domain .
this is mainly because the lack of frequency dependence of the sky model as well as the instrumental gain model .
our simulations do not include thermal noise . in practice ,
the imlin step should be applied to the residual image - cubes , averaged over total number of days required to beat thermal noise .
hence , it should produce similar decrease in the rms noise level as obtained in this paper on `` noise - less '' data .
the major feedback from all these simulations that we get is that the raw uv - data should be retained , even considering the huge data - rate issue , until the iterative real - time calibration procedures achieve the desired accuracy .
upcoming telescopes like mwa , paper , lofar , etc will also be trying to detect the eor signal in the power spectral domain . hence , a relevant extension of our present work will be to extend these analysis into the power spectral domain .
we are also planning to deal with frequency dependent and direction dependent calibration errors in subsequent publications .
ad and cc are grateful for support from the max - planck society and the alexander von humboldt foundation through the max planck forshungspreise 2005 .
the national radio astronomy observatory is a facility of the national science foundation operated under cooperative agreement by associated universities , inc .
, d.  c. , parsons , a. , bradley , r. , parashare , c. , gugliucci , n. , mastrantonio , e. , herne , d. , lynch , m. , wright , m. , werhimer , d. , carilli , c. , datta , a. , & aguirre , j. 2007 , in bulletin of the american astronomical society , vol .
38 , bulletin of the american astronomical society , 967+    , g. , de bruyn , a.  g. , brentjens , m.  a. , ciardi , b. , harker , g. , jeli , v. , koopmans , l.  v.  e. , labropoulos , p. , offringa , a. , pandey , v.  n. , schaye , j. , thomas , r.  m. , yatawatta , s. , & zaroubi , s. 2009 , , 500 , 965        , d.  s. , schwab , f.  r. , & sramek , r.  a. 1999 , in astronomical society of the pacific conference series , vol .
180 , synthesis imaging in radio astronomy ii , ed .
g.  b. taylor , c.  l. carilli , & r.  a. perley , 127+                        , g. , zaroubi , s. , bernardi , g. , brentjens , m.  a. , de bruyn , a.  g. , ciardi , b. , jeli , v. , koopmans , l.  v.  e. , labropoulos , p. , mellema , g. , offringa , a. , pandey , v.  n. , schaye , j. , thomas , r.  m. , & yatawatta , s. 2009 ,    , v. , zaroubi , s. , labropoulos , p. , thomas , r.  m. , bernardi , g. , brentjens , m.  a. , de bruyn , a.  g. , ciardi , b. , harker , g. , koopmans , l.  v.  e. , pandey , v.  n. , schaye , j. , & yatawatta , s. 2008 , , 389 , 1319    , e. , dunkley , j. , nolta , m.  r. , bennett , c.  l. , gold , b. , hinshaw , g. , jarosik , n. , larson , d. , limon , m. , page , l. , spergel , d.  n. , halpern , m. , hill , r.  s. , kogut , a. , meyer , s.  s. , tucker , g.  s. , weiland , j.  l. , wollack , e. , & wright , e.  l. 2008 , arxiv e - prints    , p. , koopmans , l.  v.  e. , jelic , v. , yatawatta , s. , thomas , r.  m. , bernardi , g. , brentjens , m. , de bruyn , g. , ciardi , b. , harker , g. , offringa , a. , pandey , v.  n. , schaye , j. , & zaroubi , s. 2009 , arxiv e - prints      , c.  j. , cappallo , r.  j. , morales , m.  f. , briggs , f.  h. , benkevitch , l. , bowman , j.  d. , bunton , j.  d. , burns , s. , corey , b.  e. , desouza , l. , doeleman , s.  s. , derome , m. , deshpande , a. , gopalakrishna , m.  r. , greenhill , l.  j. , herne , d. , hewitt , j.  n. , kamini , p.  a. , kasper , j.  c. , kincaid , b.  b. , kocz , j. , kowald , e. , kratzenberg , e. , kumar , d. , lynch , m.  j. , madhavi , s. , matejek , m. , mitchell , d. , morgan , e. , oberoi , d. , ord , s. , pathikulangara , j. , prabu , t. , rogers , a.  e.  e. , roshi , a. , salah , j.  e. , sault , r.  j. , udaya shankar , n. , srivani , k.  s. , stevens , j. , tingay , s. , vaccarella , a. , waterson , m. , wayth , r.  b. , webster , r.  l. , whitney , a.  r. , williams , a. , & williams , c. 2009 , arxiv e - prints              , r.  m. , zaroubi , s. , ciardi , b. , pawlik , a.  h. , labropoulos , p. , jelic , v. , bernardi , g. , brentjens , m.  a. , de bruyn , a.  g. , harker , g.  j.  a. , koopmans , l.  v.  e. , mellema , g. , pandey , v.  n. , schaye , j. , & yatawatta , s. 2008 , arxiv e - prints