recently there has been a considerable interest in trying to obtain atomic decompositions of the space @xmath6 .
these decompositions are usually obtained in terms of frames generated by a family of functions translated on a regular grid , and dilated by powers of a dilation matrix .
the uniformity of the grid and the structure of the dilations can be exploited to obtain very sharp results . for irregular grids and unstructured dilations or
if dilations are replaced by other transformations the situation is more complex and requires different techniques .
one method is to use the regular case and try to obtain perturbations of the grid that preserve the frame structure .
another possibility is to obtain irregular samples of the continuous transform , that have the required properties .    in this article
we study frame decompositions of the space @xmath4 using translations of a family of functions on irregular grids , and arbitrary dilations , and we even replace dilations by other transformations .
our approach is different and very general , allowing quite general constructions .
we prove the existence of smooth time - frequency frame atoms in several variables .
the setting includes as particular cases , wavelet frames on irregular lattices and with a set of dilations or transformations that do not have a group structure .
another particular case are non - harmonic gabor frames with non - uniform covering of the euclidean space .
it also leads to new constructions of wavelet and gabor frames with regular lattice translates .
one of the nice features of the proposed method is that it unifies different atomic decompositions .
for the case of regular lattices guido weiss and his group @xcite developed a very fundamental program to characterize a large class of decompositions of @xmath4 through certain equations that the generators must satisfy .
this is an important attempt to unify gabor and wavelets decompositions .
other fundamental construction of mra wavelet frames on regular lattices can also be found in @xcite , @xcite , @xcite .
our methods can be used to produce a substantial part of these systems .
a set @xmath7 is a wavelet set if the inverse fourier transform of the characteristic function of the set is a wavelet .
wavelet sets , frame wavelet sets and methods for constructing such sets have been studied recently @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite @xcite , @xcite .
our methods give constructions of wavelet sets with translations on irregular grids .
let @xmath8 and @xmath9 be countable index sets .
we consider families of functions @xmath10 and discrete sets @xmath11 such that the collection @xmath12 form a frame for @xmath4 .
the wavelet case is obtained when @xmath13 with a an expansive matrix and @xmath14 a fixed atom .
we want to stress here that our constructions are much more general , allowing for example a different invertible ( not necessarily expansive ) matrix @xmath15 for each @xmath16 for the case of orthogonal wavelets , yang wang @xcite has recently considered wavelet sets associated with arbitrary families of invertible matrices and irregular sets of translates .
he gave conditions for the existence of such wavelet sets and related them to spectral pairs .
irregular wavelet and gabor frames also have been studied as perturbations of uniform ( lattice translate ) frames and also as sampling of the continuous wavelet/ gabor transform .
see @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite,@xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , @xcite .
the approach in this article can be considered in the spirit of the classic construction in 1 dimension of smooth regular tight frames done by daubechies , grossmann and meyer in @xcite .
they found , for the case of uniform lattices , general conditions on a compactly supported smooth function @xmath17 , in order that it generates a tight gabor frame of @xmath18 . in the affine case they found necessary and sufficient conditions for a band limited function in order that it generates a smooth wavelet frame .
see also @xcite .
there were other related attempts to obtain atomic decompositions of functional spaces using very general systems .
see for example @xcite , @xcite in the context of locally compact groups .
this paper is organized as follows : section [ notation ] introduces the notation and some preliminaries .
section [ wc ] presents a theorem on wavelet construction on arbitrary , sufficiently dense , but otherwise irregular grids and with arbitrary dilation or even invertible transformation matrices .
specific constructions of such wavelets are obtained in section [ ew ] , first in the 1-d case and then in the multidimensional case .
a general theory of frame atomic decomposition of @xmath4 is obtained in section [ gen - res ] . using the concept of outer frame , reconstruction formulas for these atomic decompositions
are obtained in section [ rf ] .
throughout the paper @xmath8 and @xmath9 will denote countable index sets , and @xmath19 will stand for the function @xmath20 we will use @xmath21 to denote the lebesgue measure of a measurable set @xmath22 .    a set @xmath23 of measurable functions on @xmath24 is called a _ riesz partition of unity _ ( * rpu * ) , if there exist constants @xmath25 such that @xmath26    let @xmath27 be a family of measurable subsets of @xmath5 .
a _ riesz partition of unity * associated to @xmath28 * _ , is a set @xmath23 of measurable functions , such that    1 .
@xmath29 2 .
there exist constants @xmath25 such that @xmath30    * if @xmath31 is a rpu , then @xmath32 is also a rpu . *
if @xmath33 , we will say that @xmath31 is a _ regular partition of unity_. * if the sets in @xmath28 are essentially disjoint ( i.e. @xmath34 ) , the family @xmath35 will yield a regular partition of unity associated to @xmath28 .
* every rpu @xmath31 can be normalized to obtain a regular partition of unity by considering @xmath36 * given a family @xmath37 of measurable sets on @xmath5 , define @xmath38 where @xmath39 is the cardinal of the set @xmath40 .
the value @xmath41 is called the _ covering index _ of @xmath28 .
we now recall the definition of frame for a given close subspace @xmath42 of @xmath4 .
a set of functions @xmath43 is a _ frame _ for @xmath42 if @xmath44 and there exist constants @xmath45 , such that @xmath46    for a measurable set @xmath47 we will denote by @xmath48 the functions that have support in @xmath49 , and by @xmath50 the functions whose fourier transform has support in @xmath49 , i.e. @xmath51    we will need also the following definition .
let @xmath52 be a hilbert space .
a collection of vectors @xmath53 is an _
outer frame _ for a closed subspace @xmath42 of @xmath52 , if @xmath54 is a frame for @xmath42 , where @xmath55 is the orthogonal projection onto @xmath42 , or equivalently , there exist constants @xmath45 , such that @xmath56    related definitions to the concept of outer frames appear in @xcite , @xcite .    throughout the paper ,
we will use the following immediate and very useful fact about frames : + if @xmath57 is a frame for @xmath48 , and @xmath58 , then @xmath57 is an outer frame for @xmath59 .
our first results concerns the construction of wavelets @xmath60 with translates on the arbitrary irregular grid @xmath61 and with an arbitrary countable family of invertible @xmath2 matrices @xmath3 .    [ wavelets ] [ wavelet ] let @xmath62 be a set of finite measure , @xmath17 a function in @xmath4 and @xmath63 a family of invertible matrices .    for each @xmath64 set @xmath65 and let @xmath66 .
assume that @xmath28 is a covering of @xmath5 , @xmath67 is a rpu with bounds @xmath68 and @xmath69 and that @xmath70 .
consider @xmath71 such that for each @xmath64 , the set @xmath72 forms a frame for @xmath73 with lower and upper frame bounds @xmath74 and @xmath75 respectively . if @xmath76 and @xmath77 , then the collection @xmath78 is a wavelet frame of @xmath4 with bounds @xmath79 and @xmath80 , generated by a single function @xmath81 , where @xmath81 is the inverse fourier transform of @xmath17 .
since for each @xmath64 we have that @xmath82 forms a frame for @xmath73 with lower and upper frame bounds @xmath74 and @xmath75 respectively , an application of part  [ coro-2 ] of corollary [ transl ] for the matrix @xmath83 shows that @xmath84 forms a frame of @xmath85 with the same bounds . from the definition of @xmath86 , @xmath87 is then a frame for @xmath85 with frame bounds @xmath88 and @xmath89 .    on the other side , if @xmath90 then @xmath67 is associated to @xmath91 .
so , we can apply proposition [ genexpdil ] to @xmath92 and @xmath93 , to conclude that , @xmath94 forms a frame of @xmath4 with lower frame bound @xmath95 and upper frame bound @xmath96 .
this gives that , @xmath97 forms a frame of @xmath4 with frame bounds @xmath79 ( or @xmath98 ) and @xmath80 .
the theorem now follows from an application of the inverse fourier transform .    1 .
the set of matrices @xmath99 can be arbitrary and need not have a group structure .
the set @xmath100 can also be chosen to have a simple structure .
for example , @xmath101 , @xmath102 where @xmath103 is a rotation and @xmath104 a dilation matrix , will be used to construct directional wavelets .
an even simpler example is @xmath105 , @xmath106 , where @xmath107 is an invertible matrix which gives a construction of wavelet frames on @xmath5 .
3 .   note that @xmath17 does not need to be compactly supported .
we will use the theorem above to construct specific examples of wavelets , e.g. , directional wavelets , isotropic wavelets , etc .
* interesting particular cases of theorem  [ wavelet ] .
*    1 .   @xmath108 .
+ let @xmath109 be such that @xmath110 is a frame for @xmath73 with frame bounds @xmath111 and @xmath112
. then @xmath113 forms a wavelet frame of @xmath4 with bounds @xmath79 and @xmath80 .
2 .   @xmath114 , with @xmath115 .
+ each of the following sets are wavelet frames of @xmath4 with bounds @xmath79 and @xmath80 : @xmath116    1 .
if the set @xmath28 is a tiling of @xmath5 then the wavelets constructed above are shannon - like wavelets , thus not well localized in space . to obtain well localized space - frequency wavelets ,
@xmath117 must be constructed to be a smooth partition of unity , e.g. , at least @xmath118 as demonstrated in the examples in section [ ew ] , below
reconstruction formulas for such wavelet frames are developed in section [ rf ] .
to be able to use theorem [ wavelet ] to construct concrete examples of wavelet frames on irregular grids , we first need to construct exponential frames @xmath119 ( also called fourier frames ) for @xmath120 .
exponential frames play a central role in sampling theory for paley - wiener spaces ( also known as spaces of band - limited functions ) .
the density of a set @xmath121 and separateness of the points in @xmath122 play a fundamental role for finding exponential frames @xmath119 for @xmath120 .
a sequence @xmath123 is separated if @xmath124    there are many notions for the density of a set @xmath122 .
we start with three definitions that are due to beurling .    1 .
a lower uniform density @xmath125 of a separated sequence @xmath126 is defined as @xmath127 where @xmath128^d}\right)$ ] , where @xmath129 denotes the cardinal of the set @xmath130 .
an upper uniform density @xmath131 of a separated sequence @xmath122 is defined as @xmath132 where @xmath133^d}\right)$ ] .
3 .   if @xmath134 , then @xmath122 is said to have uniform beurling density @xmath135 .
the limits in the definitions of @xmath125 and @xmath131 exist ( see @xcite ) .    as an example , let @xmath136 be separated and assume that there exists @xmath137 such that @xmath138 , for all @xmath139
. then @xmath140 .
for the one dimensional case , beurling proved the following theorem .
( beurling ) [ b1 ] let @xmath141 be separated , @xmath142 and @xmath143 $ ] . if @xmath144 then @xmath145 is a frame for @xmath146 .
this previous result however is only valid in one dimension . for higher dimensions ,
beurling introduced the following notion :    the gap @xmath147 of the set @xmath148 is defined as @xmath149    equivalently , the gap @xmath147 can be defined as @xmath150    it is not difficult to show that if @xmath122 has gap @xmath147 , then @xmath151 . for a separated set @xmath122 , and for the case where @xmath152 is the ball @xmath153 of radius @xmath154 centered at the origin , beurling @xcite proved the following result :    [ beurlingthm ] let @xmath126 be separated , and @xmath155 .
if @xmath156 , then @xmath145 is a frame for @xmath146 .    for a very clear exposition of some of the beurling density results
see @xcite .      the construction in the following theorem ( which is a particular case of theorem  [ wavelet ] ) , generalizes a similar result of @xcite to the irregular case .
see also @xcite .
[ 1dwvlt ] let @xmath157\cup [ 1/2,1]$ ] , @xmath158 , and let @xmath159 be a real valued function such that @xmath160 $ ] , @xmath161 , and @xmath162 on @xmath163 $ ] .
assume that for each @xmath164 , the sequence @xmath165 is separated and that @xmath166 .
then for each @xmath167 , the set @xmath168 is a frame of @xmath169 , where @xmath170 .
if furthermore the sets @xmath171 are chosen such that the frame bounds @xmath74 and @xmath75 satisfy @xmath172 and @xmath173 , then the set @xmath174 where @xmath175 is a wavelet frame for @xmath18 .    the wavelet frame constructed in the theorem above is of the form + @xmath174 which is slightly different form than the one constructed in theorem  [ wavelet ] .
this discrepancy is due to a convenient choice of the irregular set @xmath176 that we have adopted in the statement of the theorem above .
note that the wavelets constructed in the theorem above are real and symmetric .
actually , if one wants @xmath81 with good decay , @xmath177 can be easily constructed to be @xmath178 , @xmath179 , even @xmath180 .    as a corollary , if we choose the sampling sets @xmath181 to be nested , i.e. , @xmath182 , we get    [ 1dwvltc1 ] let @xmath183 , @xmath184 , and @xmath177 be as in theorem [ 1dwvlt ] .
assume that the sequences @xmath185 are separated and such that @xmath186 .
if @xmath187 , then for each @xmath167 , the set @xmath188 is a frame for @xmath189 . if furthermore the frame bounds @xmath74 and @xmath75 satisfy @xmath172 and @xmath190 , then the set @xmath191 where @xmath192 is a wavelet frame for @xmath18 .
since @xmath193 we have that @xmath194)\ge 2 \#(x_{j}\cap [ -r , r]-1 ) $ ] .
thus @xmath195 . but @xmath187 , therefore @xmath196 .
the corollary then follows directly from theorem [ 1dwvlt ] .    from theorems [ 1dwvlt ] and [ wavelet ]
, we immediately get the following corollary .
[ 1dwvltc2 ] let @xmath183 , @xmath184 , and @xmath197 be as in theorem [ 1dwvlt ] .
assume that the set @xmath198 is separated and that @xmath199 .
then , the set @xmath200 is a frame of @xmath169 , and the set of functions @xmath201 where @xmath192 is a wavelet frame for @xmath18 .      1 .
_ shannon - type wavelet frames _ : we use corollary [ 1dwvltc1 ] , with @xmath202 and @xmath203}$ ] to get wavelet frames of the form @xmath204 these wavelets are not well localized since the decay at @xmath205 is @xmath206 .
shannon - type wavelet bases _ : if we choose @xmath181 such that @xmath207 , then by kadec s @xmath208-theorem , we immediately get that @xmath209 constructed above form a wavelet riesz basis for @xmath18 .
3 .   _ well localized wavelet frames _ : for faster decay of the wavelet frames , we choose @xmath177 to be a smoother function . let @xmath210}\ast\cdots\ast\chi_{[0,1]}$ ] be the b - spline of degree @xmath211 ( note that @xmath212 $ ] ) .
let @xmath213 , and @xmath214 .
then we get a wavelet frame of the form @xmath215 for this case the wavelets decay as @xmath216 .      1 .   for @xmath217 ,
let @xmath218 , and let @xmath219 .
if @xmath220 and @xmath221 , then using proposition  [ rdframes ] below for product frames , the set @xmath222 form a frame for @xmath223 ^ 2}$ ] .
let @xmath224 , and @xmath225 , then @xmath226 .
we can then use theorem [ wavelet ] to construct wavelet frames for @xmath227 : * _ shannon - type radial wavelets _ : let @xmath228 , then @xmath17 is radial .
thus the function @xmath81 defined as @xmath229 satisfies @xmath230 , where @xmath231 .
we then construct the shannon - like wavelet frame for @xmath227 , as in theorem [ wavelet ] .
a related construction of non - separable radial shannon - type frame wavelets and multiwavelets can be found in @xcite , and @xcite . * _ well localized radial wavelets _ : to construct wavelet frames with polynomial decay in space , we let @xmath232 , and construct the wavelet frames using theorem [ wavelet ] ( see figure [ fig2 ] ) .
+   radial wavelet frames that are well localized in space .
, width=188 ] 2 .
the points @xmath233 lie on an irregular grid of the form @xmath234 .
however , we may be interested in points @xmath235 that do not lie on irregular grids of the form @xmath234 . for this case , the same constructions above can be used to form wavelets frame for @xmath227 , as long as the gap @xmath236 .
3 .   as in the 1-d examples
above , we can also use corollary [ 1dwvltc1 ] to construct wavelets on irregular grids satisfying @xmath237 .
directional wavelet frames _ : [ examples-4 ] we can easily construct directional wavelet frames as follows : let @xmath238 be a region defined by @xmath239 , and define @xmath240 .
let @xmath241 , and @xmath103 be the matrix of a rotation by an angle @xmath242 .
let @xmath81 be such that @xmath243 , then we obtain the wavelet frame for @xmath227 of the form @xmath244 .
the index @xmath245 codes for the resolution of the wavelet , while the index @xmath246 codes for four possible directions .
thus the wavelet frame coefficients encode time scale as well as directional information .
clearly one can choose any number of directions and adapt the previous construction .
an obvious modification as shown in figure [ fig1 ] , yields wavelet frames with polynomial decay
. +   well - localized directional wavelet : the regions @xmath247 and @xmath248 that can be used to construct well localized directional wavelets .
, width=188 ] + very nice constructions of smooth directional wavelet frames on regular grids were obtained before in @xcite .
_ spiral _ in this example we will define a dilation covering by spiral annulus sectors . + let @xmath249 , and @xmath250 the spiral curve defined by @xmath251 for @xmath252 define @xmath253 to be the rotation of angle @xmath254 : @xmath255 $ ] .
the curve @xmath250 satisfies : @xmath256 note that for positive @xmath254 the matrix @xmath257 is expansive .
+ now we are ready to define the covering elements . set @xmath258 and @xmath259 , for some integer @xmath260 so that @xmath261 .
define the spiral annulus sector @xmath262 ( see figure [ fig3 ] ) .
so @xmath263 is compact and @xmath264 is a disjoint covering of @xmath265 .
+ choose @xmath266 sufficiently small and @xmath17 a smooth function that does not vanish in @xmath263 and with support in @xmath267 .
define @xmath268 .
select a separated set @xmath269 such that @xmath270 .
+ the set @xmath271 form a wavelet frame of @xmath227 generated by a single wavelet @xmath81 that is band - limited , with good decay and directional in frequency .
+   spiral wavelet frames .
, width=188 ] 6 .
obviously , all the constructions above can be generalized to @xmath5 for any dimension @xmath272 .
some of the wavelet frames may be associated with mras .
for example , the so called shannon wavelet frame constructed above is associated with the shannon mra @xmath273\ } $ ] , @xmath274 . in general
however , the precise relation needs further investigation .
in this section we will develop a general method for constructing time - frequency frame atoms in several variables .
this construction allows us to construct the previously introduced wavelet frames on irregular grids and with arbitrary dilation matrices or other types of transformations .
it also allows us to construct non - harmonic gabor frames on non - uniform coverings of @xmath5 as described in section [ ctfa ] below .
let @xmath275 and @xmath276 be non - empty open subsets of @xmath5 , @xmath277 is an invertible @xmath278 map , with @xmath278 inverse @xmath279 , i.e. , @xmath280 is a @xmath278 homeomorphism with @xmath278 inverse @xmath279 .
define @xmath281 and @xmath282 , where @xmath283 denotes the derivative of @xmath280 at @xmath284 , @xmath285 $ ] .
we have the following proposition :    [ chvar ] with the above notation , assume that @xmath254 is positive and @xmath286 is finite , then if @xmath57 is a frame for @xmath59 with frame bounds @xmath111 and @xmath112 , then @xmath287 is a frame for @xmath288 with frame bounds @xmath289 and @xmath290 .
assume that @xmath291 .
then for @xmath64 , by an application of the change of variables formula @xmath292 since @xmath293 is finite and bounded away from zero , @xmath294 is in @xmath59 .    using that @xmath295 is a frame of @xmath296 with bounds @xmath297 we have @xmath298 now ,
applying again the change of variables theorem , we obtain @xmath299 we have that @xmath300 = \left(\det[t^{\prime}(y)]\right)^{-1}$ ] for all @xmath301 .
thus @xmath302 that is , from , and we have @xmath303 this proves the upper inequality for @xmath304 .
the lower inequality is obtained in the same way , with the obvious modifications .
proposition [ chvar ] establishes that @xmath280 defines an isomorphism @xmath305 defined by @xmath306 , thereby transforming frames into frames . by taking @xmath280 to be a translation or a dilation
, we get the well known result :    [ transl ] let @xmath307 be an open subset of @xmath5 .
let @xmath308 be any point , and @xmath115 an invertible matrix .
we have ,    1 .
* translation * [ coro-1 ] @xmath57 is a frame for @xmath48 with bounds @xmath111 and @xmath112 , if and only if @xmath309 is a frame for @xmath310 with the same bounds . 2 .
* dilation * [ coro-2 ] @xmath57 is a frame for @xmath48 with bounds @xmath111 and @xmath112 , if and only if @xmath311 is a frame for @xmath312 with the same bounds .
in fact the corollary remains true if we only assume that @xmath183 is measurable .
part  [ coro-1 ] is a direct application of the proposition , for the case that @xmath313 .    for part  [ coro-2 ] , and
the transformation @xmath314 , the proposition tells us that @xmath315 is a frame with frame bounds @xmath316 and @xmath317 , therefore dividing each function by @xmath318 we obtain the result .
for the next theorem , we need to introduce some definitions .
let @xmath319 be a rpu with bounds @xmath68 and @xmath69 .
for each @xmath320 set @xmath321    let @xmath322 , and define @xmath323 for a given @xmath324 , we discard all those @xmath167 such that @xmath325 has measure zero
. note that if @xmath326 , then we can only claim that @xmath327 is a rpu associated to @xmath328 with constants @xmath324 and @xmath69 .
[ unionframe ] let @xmath329 .    1 .   [ unionframe-1 ]
assume that @xmath330 is a frame for @xmath331 with lower and upper frame bounds @xmath74 and @xmath75 respectively . if @xmath76 and @xmath332 , then @xmath333 is a frame for @xmath4 , with frame bounds @xmath334 and @xmath335 .
[ unionframe-2 ] assume that for each @xmath336 , @xmath330 is a frame for @xmath337 with lower and upper frame bounds @xmath74 and @xmath75 respectively . if @xmath338 and @xmath339 , then @xmath340 is a frame for @xmath341 , with frame bounds @xmath342 and @xmath335 .    1 .   given @xmath343 and denoting by @xmath344
we will first show that @xmath345 for @xmath346 where we used dominated convergence for the interchange of the integral with the sum .
the other inequality is analogous .
+ for each @xmath64 and @xmath343 , we use the fact that @xmath347 is a frame for @xmath331 , and that @xmath348 to obtain @xmath349 so summing over @xmath167 @xmath350 2 .   for the second case
we observe that the set @xmath351 forms a rpu associated to @xmath352 with bounds @xmath324 and @xmath69 , for if @xmath353 , @xmath354 furthermore , by hypothesis , @xmath355 these subspaces @xmath337 correspond to the subspaces @xmath331 defined in for the rpu @xmath356 , and therefore we can use the previous result ( applied to @xmath357 instead of @xmath4 ) to conclude that @xmath358 is a frame for @xmath359 .
the proof is complete by noting that @xmath360 .
* note that in the previous theorem , instead of choosing a frame for the subspaces @xmath331 , we could have chosen any collection of functions of @xmath4 that form an outer frame for @xmath331 . * if @xmath17 is a bounded function and @xmath361 , then it is easy to see that @xmath362 if and only if @xmath363 .
so the spaces @xmath331 defined in will coincide in most of the cases with @xmath364 .    as a very important particular case of the previous theorem
, we have the following corollaries .
[ unionframe - cor ] let @xmath365 be a family of subsets of @xmath5 , not necessarily disjoint , and let @xmath366 be a rpu with constants @xmath367 and @xmath69 not necessarily associated to @xmath28 .
assume that @xmath330 is a frame for @xmath368 with lower and upper frame bounds @xmath74 and @xmath75 respectively . if @xmath76 and @xmath332 , then    1 .
[ unionframe - cor-1 ] if @xmath67 is associated to @xmath28 ( i.e. @xmath369 ) , then @xmath333 is a frame for @xmath370 , with frame bounds @xmath334 and @xmath335 .
[ unionframe - cor-2 ] if instead @xmath371 , then also @xmath333 is a frame for @xmath372 , with frame bounds @xmath342 and @xmath335 .
denote by @xmath373 the orthogonal projection on the subspace @xmath331 defined in .
since @xmath374 forms a frame for @xmath368 with bounds @xmath74 and @xmath75 , then @xmath375 forms a frame for @xmath331 with the same bounds .
so , by theorem  [ unionframe ] , @xmath376 forms a frame of @xmath377 with frame bounds @xmath79 and @xmath80 .
we obtain part @xmath378 of the corollary observing that : @xmath379    the second part is a consequence of the fact that @xmath380 is a rpu associated to @xmath28 with bounds @xmath324 and @xmath69 and that @xmath381 .
therefore we can apply the second part of theorem [ unionframe ] .
* for the first case , the set @xmath382 is not necessarily a frame for @xmath383 or even for @xmath384 , even though @xmath385 is a frame for @xmath386 .
* in the second case , the subspace @xmath331 is contained in @xmath387
. * note that if @xmath67 is a regular partition of unity , then the frame bounds for the frame constructed in the corollary , are @xmath111 and @xmath112 .
[ new - cor ] let @xmath388 be a rpu with constants @xmath367 and @xmath69 and @xmath389 and @xmath390be coverings of @xmath5 such that @xmath391 and @xmath392 for all @xmath64 and some constant @xmath393 .
assume that @xmath330 is a frame for @xmath368 with lower and upper frame bounds @xmath74 and @xmath75 respectively . if @xmath76 and @xmath332 , then @xmath333 is a frame for @xmath6 , with frame bounds @xmath342 and @xmath335 .    as in the proof of the theorem , we have the inequalities @xmath394 on the other side , @xmath395 and then @xmath396 and @xmath397    the following proposition is a direct application of fubini s theorem , and allows us to construct frames in a product space . by induction the theorem
can be generalized to hold for any finite number of factors .
[ rdframes ] let @xmath398 and @xmath399 be measurable subsets of @xmath400 and @xmath401 respectively , and let @xmath402 and @xmath403 be frames for @xmath404 and @xmath405 with frame bounds @xmath406 and @xmath407 . then @xmath408 is a frame of @xmath409 with frame bounds @xmath410 .    for any function @xmath411 in @xmath409 ,
@xmath412 which yields the lower frame bound @xmath413 .
the upper frame bound @xmath414 can be obtained in a similar fashion .
we now particularize our previous results to frames of the form @xmath415 , where @xmath17 is a fixed function . using the fourier transform ,
these types of frames allow us to construct wavelets @xmath416 with translates on the arbitrary irregular grid @xmath122 and with an arbitrary countable family of invertible @xmath2 matrices @xmath3 ( cf .
theorem  [ wavelet ] ) .
first as a particular case of corollary [ unionframe - cor ] , we obtain the following proposition .
[ genexpdil ] assume that @xmath417 forms a covering of @xmath5 , and let @xmath418 be a rpu with bounds @xmath68 and @xmath69 associated to @xmath28 .
assume also that @xmath419 is a frame for @xmath368 with lower and upper frame bounds @xmath74 and @xmath75 respectively .
if @xmath76 and @xmath77 , then @xmath420 is a frame for @xmath4 with frame bounds @xmath79 and @xmath80 .
if @xmath421 , @xmath86 , and @xmath422 are chosen such that @xmath423 , @xmath252 , @xmath424 and @xmath425 , where @xmath250 is a lattice in @xmath5 , then we obtain the standard gabor or weyl - heisenberg frames .
thus in general , the construction above can be viewed as non - harmonic gabor frames with variable windows @xmath422 .
the following wavelet frame construction is a direct application of the previous proposition .
taking @xmath426 , the set @xmath427 is a wavelet frame for @xmath4 with frame bounds @xmath79 and @xmath80 .
therefore , there exists a dual frame @xmath428 such that @xmath429    we can also use proposition [ genexpdil ] and beurling theorem [ b1 ] , to obtain a non - harmonic gabor frame as the following example shows :    [ gf ] [ non - harmonic gabor frames ] let @xmath430 $ ] and let @xmath431}\ast\chi_{[0,1]}\ast\chi_{[0,1]}\ast\chi_{[0,1]}$ ] be the b - spline of degree @xmath432 .
clearly , @xmath433 .
let @xmath185 be sets in @xmath434 chosen such that for each @xmath167 , @xmath435 and such that the frame bounds @xmath74 and @xmath75 satisfy @xmath172 and @xmath173 . then using proposition [ genexpdil ] and beurling theorem [ b1 ] , we obtain a non - harmonic gabor frame of the form @xmath436 .
obviously , we can also use a non - uniform partition of @xmath434 and get generalized non - harmonic gabor frames as discussed earlier .
if @xmath62 be a measurable subset of @xmath5 , and if the family of sets @xmath28 is defined by means of expanding or contracting @xmath183 , then we obtain the theorems of section  [ wc ] as particular cases of proposition  [ genexpdil ] .      by corollaries  [ transl ]  and  [ unionframe - cor ]
, we can construct a frame for @xmath437 starting from a frame for @xmath438 where @xmath104 is _ any _ subset of @xmath5 with nonempty interior . hence to build a frame for @xmath439 it is enough to start with a frame for @xmath440 , where @xmath275 is an open disk . specifically , since any bounded measurable set can be covered by a finite number of translates of @xmath275 , we can use corollary  [ transl](1 ) to find a frame for @xmath439 .
we can also expand @xmath275 until it covers @xmath183 and use corollary  [ transl](2 ) .
this shows that there are many ways to construct frames for @xmath439 starting from a frame for @xmath440 .
obviously , the particular construction will depend on the application .
we will describe now two particular constructions .    1 .
given a frame for @xmath441 where @xmath104 is a measurable subset of @xmath5 , we will construct a frame for @xmath442 .
let @xmath250 be a regular lattice in @xmath5 ( i.e. @xmath443 , where @xmath103 is an invertible @xmath2 matrix with real entries ) , and let @xmath104 be a measurable subset of @xmath5 such that @xmath444 with a finite covering index .
let @xmath183 be a measurable subset of @xmath5 and define @xmath445 for @xmath446 .
let @xmath447 . by corollary [ transl ] ,
if @xmath448 is a frame for @xmath449 , then @xmath450 is a frame for @xmath451 .
hence , @xmath452 is also a frame for @xmath453 .
therefore , by corollary [ unionframe - cor ] , @xmath454 is a frame for @xmath73 .
as an example , when @xmath183 is a measurable subset of @xmath5 , @xmath455^d$ ] , and @xmath456 , we have that if @xmath457^d } : k \in k\}$ ] is a frame for @xmath458^d)$ ] , then @xmath459 is a frame for @xmath73 ( recall that @xmath460^d + \gamma)$ ] ) .
+ it is easy to see that the covering requirement @xmath461 in the previous construction is not restrictive .
specifically , we only need @xmath462 , where @xmath254 is any positive real number . furthermore
, the construction remains valid if for each @xmath463 we choose a _ different _ set @xmath464 such that @xmath465 is a frame for @xmath449 , and @xmath466 and @xmath467 , where @xmath468 and @xmath469 are the lower and upper frame bounds of @xmath470 .
finally , note that if @xmath471 we obtain a frame for @xmath4 .
2 .   if @xmath183 is bounded we can construct frames for @xmath442 using theorem  [ beurlingthm ] .
let @xmath472 , and @xmath473 such that @xmath474 .
let @xmath475 be such that @xmath476 .
then using beurlings theorem ( [ beurlingthm ] ) we obtain that @xmath477 is a frame of @xmath478 .
so , using corollary  [ transl](1 ) , @xmath479 is an outer frame of @xmath442 , and @xmath480 is a frame for @xmath442 .      in view of the previous results
, we will be interested in constructing particular kinds of riesz partitions of unity associated to special coverings of the space .
the next results provide the necessary tools to accomplish this task .
if @xmath107 is a @xmath481 matrix , we will say that @xmath107 is _ expansive _ , if @xmath482 for every eigenvalue @xmath356 of @xmath107 .
we will use the following known result ( see for example @xcite , pg .
let @xmath40 be in @xmath483 and @xmath266 .
there exists a matrix norm @xmath484 such that @xmath485 where @xmath486 is the spectral radius of the matrix @xmath40 , and there exists a norm @xmath487 in @xmath488 such that @xmath489    as a consequence of this lemma , if @xmath107 is an expansive matrix , then there exists a norm @xmath490 in @xmath488 such that @xmath491 and therefore @xmath492 in particular for every @xmath493 , @xmath494    [ prop * * ] let @xmath107 be a @xmath2 expansive matrix and @xmath495 a bounded set such that    1 .
there exists @xmath266 such that @xmath496 .
[ prop**1 ] 2 .
@xmath497.[prop**2 ]    then @xmath498 , the covering index of the family @xmath499 , is finite , i.e. there exists an integer @xmath500 such that @xmath501 .    for @xmath502
define @xmath503 since @xmath504 , it is enough to prove that @xmath505 and @xmath506 are uniformly bounded in @xmath5 .
we will see that @xmath505 is uniformly bounded .
a similar argument proves the claim for @xmath506 .
it is easy to see that @xmath507 thus , if @xmath505 is bounded on @xmath276 , then , by ( [ prop**2 ] ) @xmath505 is bounded in @xmath265 with the same bound .
now , using ( [ prop**1 ] ) , we see that there exist @xmath508 such that @xmath509 , for all @xmath510 .
fix @xmath510 .
if @xmath511 , then there exists @xmath512 such that @xmath513 and @xmath514 and therefore @xmath515 .
but since @xmath516 , necessarily there exists @xmath517 such that for every @xmath518 . hence for every @xmath510 , @xmath519 .
the next proposition shows the construction of a large class of riesz partitions of unity , for families of sets obtained by dilation of a compact set .
[ prop2 ] let @xmath62 be a compact set and @xmath107 a @xmath2 expansive matrix such that    1 .
@xmath520 [ prop2 - 1 ] 2 .
[ prop2 - 3]@xmath521 .
let @xmath17 be any measurable function , and @xmath522 some constants such that    1 .
@xmath523 2 .
[ prop2 - 5 ] @xmath524 on @xmath183 3 .
@xmath525 on @xmath526 , where @xmath527 and @xmath528    then the family of functions @xmath529 is a rpu associated to @xmath530 .
if @xmath531 , by ( [ prop2 - 3 ] ) there exists @xmath164 and @xmath532 such that @xmath533 , so by ( [ prop2 - 5 ] ) , @xmath534 now by proposition  [ prop * * ] , the covering index @xmath535 of the family @xmath536 is finite . using ( [ prop2 - 1 ] ) and ( [ prop2 - 3 ] )
, we see that @xmath537 and since @xmath538 which proves the proposition .
this proposition generalizes easily to the case where we replace for each @xmath167 , @xmath539 by an invertible matrix @xmath15 in such a way that @xmath540 is a covering of @xmath265 with finite index .
the next lemma shows , that the assumption of having a compact set that covers @xmath5 by dilations is actually necessary , if one wants to construct a rpu associated to a covering of @xmath5 of the form @xmath541 , where @xmath276 is a bounded open set .    if @xmath107 is invertible , and @xmath276 is a bounded open set such that    1 .
[ lem - i ] there exists @xmath266 with @xmath542 2 .
[ lem - ii ] @xmath497    then there exists a compact set @xmath543 such that @xmath544 .
it is clearly enough to find @xmath183 such that @xmath545 .
let @xmath546 we will prove ( by the contradiction ) that for some @xmath547 , @xmath548 covers @xmath276 by dilations by @xmath107 .
assume that for each @xmath547 , there exists @xmath549 such that @xmath550 .
since @xmath551 is compact , and @xmath552 , there exists a subsequence @xmath553 and @xmath554 such that @xmath555 . by ( [ lem - i ] ) and ( [ lem - ii ] ) there exists @xmath517 such that @xmath556 , and since @xmath557 is open , @xmath558 for @xmath559 .
set @xmath560 and @xmath561 , then @xmath562 for @xmath563 .
choose @xmath564 small enough such that @xmath565 , and so @xmath566 .
let @xmath567 be such that if @xmath568 , then @xmath569 .
then for @xmath570 , and @xmath571 we have @xmath572 and thus @xmath573 so @xmath574 for all @xmath211 such that @xmath575 .
this contradicts our assumption that @xmath576 , since , if @xmath568 and @xmath577 , we have @xmath578 , and so @xmath579 .
we show next that it is not difficult to construct bounded sets that cover @xmath5 by dilations .
let @xmath495 be a bounded set such that @xmath580 , and @xmath107 an expansive @xmath481 matrix .
set @xmath581 , then @xmath582 is a covering of @xmath265 with finite covering index . furthermore , if @xmath583 then the sets @xmath584 are disjoint .
choose @xmath266 such that @xmath585 .
let @xmath586 be an arbitrary point . by equation   @xmath587 , for some norm @xmath487 in @xmath488 .
so there exists a positive integer @xmath211 such that @xmath588 , @xmath589 , i.e. , @xmath590 .    since @xmath276 is bounded , and @xmath591 , there exists @xmath592 such that @xmath593 and @xmath594 , i.e. , @xmath595 .
the finiteness of the covering index follows from proposition  [ prop * * ] , and if @xmath583 the disjointness property follows immediately .    *
if we want to obtain coverings of @xmath596 with compact sets , we can choose the set @xmath49 where @xmath183 is the set constructed in the previous lemma .
clearly the sets @xmath597 cover @xmath596 , and if @xmath598 , then @xmath597 are almost disjoint , i.e.  @xmath599 . *
even though for _ any _ set @xmath276 it is true that @xmath600 , in general it is not true that the family @xmath601 is disjoint , as the following example shows : let @xmath602 $ ] and @xmath276 the rectangle @xmath603\times[-1,1]$ ] .
we can combine the results of this section with theorem  [ wavelet ] , to obtain the following recipe to construct smooth , wavelet frames of @xmath5 associated to a single dilation matrix @xmath107 and an irregular grid .
the wavelets obtained by this method can be constructed to have polynomial decay of any degree , as exemplified in section [ wf1d ] .
* recipe : *    * select a bounded set @xmath495 such that @xmath580 and @xmath604 .
* select a function @xmath17 of class @xmath605 such that @xmath606 on @xmath607 and @xmath608 for some small @xmath266 . *
consider a set @xmath609 , such that @xmath122 is separated and @xmath610 , where @xmath611 .
then the following collection is a frame of @xmath4 : @xmath612 where @xmath81 is the inverse fourier transform of @xmath17 .
we first note that although the set @xmath613 in theorem  [ wavelet ] is a wavelet frame for @xmath4 , it is not in general true that for a fixed @xmath167 the set @xmath614 is a frame , unless the wavelet frame is of shannon - type , hence , not well - localized .
thus , for well - localized wavelets , it appears that the reconstruction of a function @xmath343 from the wavelet coefficients @xmath615 can not be obtained in a stable way by first reconstructing at each level @xmath167 and then obtaining @xmath411 by summing over all levels @xmath167 .
however what follows shows that this is still possible :      @xmath618 this , together with the assumptions that @xmath619 is a frame for @xmath620 for each @xmath167 , and that @xmath621 , permits to reconstruct @xmath622 using a dual frame @xmath623 of @xmath72 .
note that , since @xmath624 is unitary then @xmath625 is a dual frame of @xmath626 for @xmath368 .
thus , it is always possible to reconstruct each @xmath627 in a stable way and then obtain @xmath411 by summing up over all levels @xmath167 .
one drawback is that the dual frame may not be well - localized , since each @xmath628 may be discontinuous at the boundary of @xmath183 . to treat this problem
we simply note that @xmath629 multiplying both sides by @xmath422 we obtain    @xmath630 .
if we choose @xmath17 to be in @xmath631 , and if @xmath628 is in @xmath632 , then @xmath633 will be in @xmath631 and therefore decays polynomially in space .
hence the partial sums of the series ( [ recfor ] ) will have good convergence properties .
we can then sum equation   over @xmath64 to obtain @xmath634 , and then divide by @xmath635 to obtain @xmath636 .
c.  k.  chui , w.  he , j.  stckler and q.  sun , _ compactly supported tight affine frames with integer dilations and maximum vanishing momemnts _ , advances computational math .
, special issue on frames , * 18 * ( 2003 ) , no .  2 , 159187 .
m.  papadakis , g.  gogoshin , i.  a.  kakadiaris , d.  j.  kouri , and d.  k.  hoffman , _ non - separable radial frame multiresolution analysis in multidimensions and isotropic fast wavelet algorithms _ , spie - wavelet x , to appear .


our current hypothesis is that the galaxies and black holes observed today originated over 13 gyr ago , growing from seeds of primordial density perturbations .
one can test this hypothesis by studying the star formation rate ( sfr ) history , from the epoch of reionization ( eor ) at redshifts @xmath21 , through the peak era at @xmath22 ( hopkins & beacom 2006 ) down to the present epoch .
the feedback of ionizing radiation , kinetic energy , and heavy elements leaves imprints on early stars , supernovae , and galaxies , providing a  fossil record " that can be detected through abundances in galactic halo stars and the intergalactic medium ( igm ) and in the distributions of mass , metallicity , and luminosity of galaxies .
determining when and how the universe was reionized by these early sources have been important questions for decades ( gunn & peterson 1965 ; sunyaev 1977 ; robertson et  al .
it has been suggested that igm reionization was complete by @xmath23 ( fan et  al .
2001 ; gnedin & fan 2006 ; fan et  al .
2006 ; hu & cowie 2006 ) , based on strong ly@xmath2  absorption from neutral hydrogen along lines of sight to qsos at @xmath24 .
becker et  al .
( 2007 ) and songaila ( 2004 ) used transmission of the ly@xmath2(and ly@xmath25 ) forest out to @xmath26 and @xmath27 , respectively , to suggest a smoothly decreasing ionization rate toward higher redshifts .
recent surveys of high - redshift galaxies and ly@xmath2  emitters ( bouwens et  al .   2011a ; ouchi et  al .
2010 ; kashikawa et  al .  2011 ; ono et  al .  2012
; schenker et  al .
2012 ) infer an increasing igm neutral fraction from the declining populations between @xmath3 .
further evidence for an increasing neutral fraction comes from the decreasing sizes of ionized  near zones " associated with quasars between @xmath28 and @xmath29 ( carilli et  al .
2010 ) and from the ly@xmath2  damping wing in the transmission profile toward the newly discovered quasar at @xmath30 ( mortlock et  al .
these studies all suggest that the igm is becoming increasingly neutral between @xmath31 , marking the end of cosmic reionization when ionized regions overlap and percolate .
whether the epoch of full reionization occurs at @xmath0 is still not ascertained .
a contrasting estimate of the eor comes from the measured optical depth , @xmath32 , to electron scattering of the cosmic microwave background ( cmb ) in _ wmap-7
_ observations ( larson et  al .  2011 ;
komatsu et  al .
the error bars come from the central 68% in the marginalized cumulative distribution . using additional cosmological parameter constraints , they infer single - epoch reionization at @xmath33 .
although such a high redshift could be explained with @xmath34cdm simulations and modeled sfr histories ( choudhury & ferrara 2005 ; trac & cen 2007 ; shull & venkatesan 2008 ) , the cmb observations are at variance with optical surveys that suggest late reionization , unless reionization is a process that extends to higher redshifts .    for a fully ionized igm , including both h and he , the optical depth @xmath14 for @xmath35 ( see section 2.1 ) . in the cmb analysis ( larson et  al .
2011 ) , marginalization of @xmath36 with other cosmological parameters allows the possibility of lower optical depth , with a reionization epoch as low as @xmath37 at 95% c.l .
one can also invoke a partially ionized igm at @xmath38 , as discussed by many groups ( cen 2003 ; venkatesan et  al .
2003 ; ricotti & ostriker 2004 ; benson et  al .
2006 ; shull & venkatesan 2008 ) .
however , even with the recent progress in finding high-@xmath12 galaxies , we still do not know whether galaxies are the sole agents of reionization
. current observations of high-@xmath12 galaxies leave open several ionization scenarios , some involving simple hydrogen reionization at @xmath39 and  reionization at @xmath40 , and others with more complex ionization histories ( bolton & haehnelt 2007 ; venkatesan , tumlinson & shull 2003 ; cen 2003 ) that depend on sfrs at @xmath41 .
shull & venkatesan ( 2008 ) demonstrated how the cmb optical depth constrains the sfr and igm metallicity history at @xmath17 , and trenti & shull ( 2010 ) quantified the metallicity - driven transition from population  iii ( metal - free stars ) to population  ii ( stars formed from metal - enriched gas ) .
several recent observations provide valuable constraints on the luminosity function of high-@xmath12 galaxies .
the number density of galaxies appears to drop rapidly at @xmath17 ( bouwens et  al .
2009 , 2010a , b , 2011a , b , c ) . with a comoving sfr density @xmath42 at @xmath43 ( gonzlez et  al .  2010 ; bouwens et  al .
2011a ) , the observed galaxies do not produce enough ionizing photons in the lyman continuum ( lyc ) to maintain a photoionized igm against recombinations .
however , the luminosity function is steep , and the total lyc budget requires extrapolation to low - luminosity galaxies ( trenti et  al .  2010 ; bouwens et  al .
moreover , the conversion from sfr to lyc production rate relies on insecure calibrations from theoretical models and comparison with high - mass , low - metallicity stars .
we revisit the calculation of lyc photon production and assess the high-@xmath12 galaxy contribution to reionization .
we also analyze several factors , such as the photon escape fraction ( @xmath11 ) , igm clumping factor ( @xmath44 ) , and electron temperature ( @xmath45 ) , which enter the calculation of the  critical star formation rate " ( @xmath46 ) necessary to maintain a photoionized igm .    in section 2 , we calculate @xmath46 ( @xmath47 mpc@xmath48 ) in a filamentary igm , equating the production rate of lyman continuum ( lyc ) photons with the hydrogen recombination rate .
the photoionization rate depends on the mass function of stellar populations , their evolutionary tracks and stellar atmospheres , and the escape fraction , @xmath11 of lyc photons away from their galactic sources .
the recombination rate depends on the density and temperature of the igm , properties we explore with cosmological simulations . in section 3
, we give our results for the critical sfr at @xmath49 and present our new sfr simulator , a user - friendly interface for calculating @xmath50 and @xmath51 . in section 4 , we discuss the implications for the hydrogen eor .
consistency between high - redshift galaxies and cmb optical depth appears to require @xmath0 and a partially ionized igm at @xmath17 .
the peak signal from redshifted 21-cm emission would likely occur during the heating period between @xmath52 ( 145163 mhz ) when the hydrogen neutral fraction @xmath53 ( pritchard et  al .  2010 ; lidz et  al .
we denote by @xmath1 ( @xmath54 ) the global star formation rate per co - moving volume . using a simple argument ( madau et  al .
1999 ) , balancing photoionization with radiative recombination , we estimate the critical sfr to maintain igm photoionization at @xmath17 , assuming that the lyc photons are produced by populations of massive ( ob - type ) stars .
because the mass in collapsed objects ( clusters , groups , galaxies ) is still small at high redshift , the igm contains most of the cosmological baryons , at mean density @xmath55 for a hubble constant denoted @xmath56 , we adopt the wmap-7 ( plus bao + @xmath57 ) parameters , @xmath58 and @xmath59 ( komatsu et  al .
2011 ) relative to a critical density @xmath60 g  cm@xmath48 . from the corresponding helium mass fraction @xmath62 ( peimbert et  al .
2007 ) , we adopt a mean hydrogen number density , @xmath63    in a fully ionized igm , the cmb optical depth back to @xmath64 can be written as the integral of @xmath65 , the electron density times the thomson cross section along proper length , @xmath66 \ ; dz     \ ; , \ ] ] for a standard @xmath34cdm cosmology ( @xmath67 ) with @xmath68^{1/2}$ ] .
this integral can be done analytically ( shull & venkatesan 2008 ) , @xmath69       \left [ \ { \omega_m ( 1+z_{\rm rei})^3 + \omega_{\lambda } \}^{1/2 } - 1 \right ] \ ; .\ ] ] in the high - redshift limit , when @xmath70 , this expression simplifies to @xmath71 ( 1+z_{\rm rei})^{3/2 } \approx ( 0.0521 ) \left [ \frac { ( 1+z_{\rm rei})}{8 } \right ] ^{3/2}\ ] ] independent of @xmath72 to lowest order . the helium and electron densities are written @xmath73 and @xmath74 for singly ionized helium , where @xmath75 by number . to these formulae , we add @xmath76 , from electrons donated by  reionized at @xmath77 ( shull et  al .
helium therefore contributes @xmath788% to @xmath36 , and a fully ionized igm produces @xmath14 , 0.060 , and 0.070 back to redshifts @xmath79 , 8 , and 9 , respectively .
a comoving volume of 1 mpc@xmath80 contains @xmath81 hydrogen atoms .
our simple ionization criterion requires a sfr density that produces a number of lyc photons equal to @xmath82 over a hydrogen recombination time , @xmath83^{-1}$ ] .
the hydrogen case - b recombination rate coefficient ( osterbrock & ferland 2006 ) is @xmath84 , scaled to an igm temperature @xmath85 . for typical igm ionization histories and photoelectric heating rates ,
numerical simulations predict that diffuse photoionized filaments of hydrogen have temperatures ranging from 5000  k to 20,000  k ( dav et  al .  2001 ; smith et  al .
these are consistent with temperatures inferred from observations of the ly@xmath2  forest at @xmath86 ( becker et  al .
2011 ) .
owing to gravitational instability , a realistic igm is inhomogenous and filamentary .
semi - analytical models of the reionization of the universe often adopt a  clumping factor " , @xmath87 , to account for inhomogeneity in estimates of the enhanced recombination rate in denser igm filaments .
the clumping factor therefore plays an important role in computing the critical sfr density needed to maintain the reionization of the universe .
the clumping factor is also used in numerical simulations to implement `` sub - grid physics '' , in which changes in the density field occur on scales below the resolution of the simulation and are also approximated by the factor @xmath44 ( gnedin & ostriker 1997 ; madau et  al .  1999 ; miralda - escud et  al .
2000 ; miralda - escud 2003 ; kohler et  al .
2007 ) .
we correct the recombination time for density variations scaled to a fiducial @xmath88 , found in the simulations described below . at @xmath43 ,
the igm filaments have electron density @xmath89 ^ 3 c_h$ ] , and the characteristic times for hydrogen recombination and hubble expansion are , @xmath90^{-3 }   \ ; , \\
t_{\rm h }    & \approx & [ h_0 \omega_m^{1/2 } ( 1+z)^{3/2}]^{-1 } \approx ( 1.18~{\rm gyr } )                     \left [ \frac { ( 1+z)}{8 } \right]^{-3/2 }   \ ;   .
\end{aligned}\ ] ] in our calculations , we express the reionization criterion as @xmath91 , where @xmath46 is the critical sfr density ( @xmath54 ) and @xmath92 is the conversion factor from @xmath1 to the lyc production rate ( see section 2.2 ) .
we define @xmath11 as the fraction of lyc photons that escape from their galactic sources into the igm ( dove & shull 1994 ) .
recent statistical estimates ( nestor et  al .
2011 ) suggest that @xmath93 for an ensemble of 26 lyman - break galaxies and 130 ly@xmath2 emitters at @xmath94 , and it could be higher for the lower - mass galaxies that likely dominate the escaping lyc at @xmath24 ( fernandez & shull 2011 ) . the lyc production efficiency , @xmath92 , is expressed in units @xmath95 photons per @xmath96 of star formation , since typical massive stars emit @xmath97 lyc photons over their lifetime . to evaluate @xmath92
, we convert the sfr ( by mass ) into numbers of ob - stars and compute the total number of ionizing photons produced by a star of mass @xmath98 over its lifetime .
we then integrate over an imf , @xmath99 , with a range @xmath100 .
the standard mass range is 0.1 @xmath96 to 100 @xmath96 , but changes in the mass range and imf slope will affect the lyc production substantially .
for example , differences in the imf have been associated with higher jeans masses in low - metallicity gas in the high - redshift igm ( abel et  al .
2002 ; bromm & loeb 2003 ) , and tumlinson ( 2007 ) and smith et  al .
( 2009 ) noted the potential influence of cmb temperature on modes of high - redshift star formation .
consequently , most cosmological simulations or calculations include a metallicity - induced imf transition ( trenti & shull 2010 ) between low - metallicity population  iii star formation and metal - enhanced population  ii .      the number of lyc photons produced per total mass in star formation
depends on the imf of the stellar population and is given by the conversion factor @xmath92 .
we calculate this conversion by integrating the total number of lyc photons produced over the entire mass in star formation .
@xmath101 here , @xmath102 is the lifetime - integrated number of lyc photons as a function of mass calculated from stellar atmosphere models and evolutionary tracks .
the imf , @xmath103 , is integrated over the mass range @xmath100 , where @xmath98 is expressed in solar units . in the sfr simulator discussed in section 3.3 , the user can choose between a normal and broken imf power law .
the lower integration limit in the numerator , @xmath104 , is the mass at which stars no longer produce significant amounts of lyc photons .    in our calculator , we use and compare two models that calculate @xmath102 . first , using stellar atmospheres and evolutionary tracks ( r.  s.  sutherland & j.  m.  shull , unpublished ) , we find that , over its main - sequence and post - main - sequence lifetime , an ob star of mass @xmath98 produces a total number @xmath105 of ionizing photons , where @xmath106 over the mass range @xmath107 and for metallicities @xmath108 .
we have fitted our results to the form @xmath109 , for @xmath110 ( the mass @xmath104 defines the effective lower limit for stars that produce significant numbers of lyc photons ) . for metallicities in the range @xmath111 ( where @xmath112 is the solar metallicity ) and for masses @xmath113
, the fitted coefficients are @xmath114 and @xmath115 .
thus , for @xmath116 , we have @xmath117 and @xmath118 , while for @xmath119 we have @xmath120 and @xmath121 . inserting these values for @xmath102 into equation ( 8)
, we integrate these lyc photon yields over the imf , to derive the conversion coefficient , @xmath92 , from _ total mass _ in star formation to _ total number _ of lyc photons produced ( in units of @xmath95 photons per @xmath96 ) .
the integrals in eq .
( 8) can be done analytically as functions of the imf parameters and metallicity . for a salpeter imf ( @xmath122 , @xmath123 )
we find : @xmath124 ( for @xmath125 ) , @xmath126 ( for @xmath127 ) , and @xmath128 ( for @xmath129 ) .
adopting a typical ( low-@xmath130 ) value @xmath131 ( @xmath132 photons per @xmath96 ) , we can solve for the critical sfr for reionization , scaled to clumping factor @xmath6 , escape fraction @xmath7 , and gas temperature @xmath133 : @xmath134^{3 }      \left [ \frac { c_h/3 } { f_{\rm esc}/0.2 } \right ]     \left [ \frac { 0.004 } { q_{\rm lyc } } \right ] t_4^{-0.845 } \ ; .\ ] ] our chosen value of @xmath88 is consistent with recent downward revisions ( pawlik et  al .
2009 ) and with our numerical simulations discussed in section 2.3 .
escape fractions @xmath135 have been inferred from observations at @xmath40 ( shapley et  al .
2006 ; nestor et  al .
2011 ) and theoretical expectations ( fernandez & shull 2011 ) . equation ( 9 ) agrees with prior estimates ( madau et  al .
1999 ) when adjusted for our new scaling factors , in particular the ratio @xmath136 ) .
their earlier paper assumed @xmath137 , @xmath138 , @xmath139 , @xmath140 , and @xmath141 .
our fiducial redshift has increased from @xmath141 to @xmath142 , appropriate for the new discoveries of high - redshift galaxy candidates for reionization .
expressed with the same coefficients in eq .
( 9 ) , their coefficient would be nearly the same , @xmath143 at @xmath142 .
one of the improvements in our formulation is to better identify the dependences on the physical parameters ( @xmath44 , @xmath11 , @xmath45 ) and the sfr - to - lyc conversion factor ( @xmath144 , which can change with different imfs and atmospheres .
a related calculation is the production rate of ionizing ( lyc ) photons per unit volume , needed to balance hydrogen recombinations . with the same assumptions as above ,
this is @xmath145 ^ 3 t_{4}^{-0.845 }     \left ( \frac { c_h}{3 } \right )    \ ;   .\ ] ] for standard high - mass stars ( o7  v , solar metallicity ) each with lyc photon luminosity @xmath146  s@xmath147 , this rate corresponds to an o - star space density @xmath148 at @xmath142 .
schaerer ( 2002 , 2003 ) also computed models of population iii and low - metallicity stars based on non - lte model atmospheres and new stellar evolution tracks and evolutionary synthesis models .
the lifetime total number of lyc photons produced per star of mass @xmath98 is @xmath149 . for stars of mass parameter @xmath150 ,
the number of ionizing photons , @xmath151 , emitted per second per star , is given by : @xmath152     = \left\ { \begin{array}{lll }      43.61 + 4.90x - 0.83x^2 &    z = 0 ,   &    9 - 500 m_\odot , \\      39.29 + 8.55x    &    z=0 ,     &    5 - 9 m_\odot,\\      27.80 + 30.68x-14.80x^2 + 2.50x^3    & z = 0.02 z_\odot , &    7 - 150 m_\odot ,
\\      27.89 + 27.75x-11.87x^2 + 1.73x^3    &    z= z_\odot , &    7 - 120 m_\odot \\ \end{array }       \right .
\label{eq : qvalues}\ ] ]    and the star s lifetime @xmath153 is given by :    @xmath154   = \left\ { \begin{array}{ll }      9.785 - 3.759x+1.413x^2 - 0.186x^3   &    z = 0 , \\
9.59 - 2.79x+0.63x^2   &    z = 0.02 z_\odot , \\
9.986 - 3.497x+0.894x^2    &    z = z_\odot \end{array }       \right .
\label{eq : lifetimes}\ ] ]    sample values of @xmath92 and @xmath46 for different imfs are shown in table  1 , for both model atmospheres .
one can obtain different values of @xmath92 for power - law imfs by varying their high - mass slope @xmath25 and the minimum and maximum masses .
increases in @xmath92 translate into _ decreases _ in critical sfr .
for example , at @xmath155 , if the minimum mass is raised to @xmath156 with a salpeter slope ( @xmath157 ) , the lyc production factor rises to @xmath158 , some 2.5 times higher than the equivalent model , @xmath128 for @xmath159 .
if the imf is flatter , @xmath160 instead of 2.35 , with minimum mass fixed at @xmath159 , @xmath161 ( 4.6 times higher ) . at the upper end of the imf ,
if one increases @xmath162 from 100 to 200 @xmath96 , keeping @xmath157 and @xmath163 , one finds @xmath164 ( 1.4 times higher ) .
all of these variations can be explored with our reionization / sfr simulator , described in section 3.3 .
the critical sfr of @xmath165 can be understood from simple arithmetic . over the 386 myr recombination time
( eq .  [ 6 ] ) with @xmath6 , @xmath7 , and @xmath142 , approximately 20 million stars are formed per comoving mpc@xmath80 in a salpeter imf ( @xmath166 ) with mean stellar mass @xmath167 . of these stars , a small fraction , @xmath168 ,
are massive ob stars ( @xmath169 . at an average of @xmath170 lyc photons and escape fraction @xmath7 , these @xmath7814,000 ob stars will produce a net ( escaping ) @xmath171 lyc photons over their lifetime .
these photons are sufficient to ionize @xmath81 hydrogen atoms mpc@xmath48 .
the simulations used in this study were performed with ` enzo ` , an eulerian adaptive mesh - refinement ( amr ) , hydrodynamical + n - body code ( bryan & norman 1997 ;
oshea et  al .  2004 , 2005 ) .
smith et  al .
( 2011 ) enhanced this code by adding new modules for star formation , primordial chemistry , and cooling rates consistent with ionizing radiation , metal transport , and feedback .
the ionizing background is spatially constant and optically thin , but variable in redshift . at this stage , we have not implemented radiative transfer or spectral filtering by the igm .    for the clumping factor calculation , we ran a simulation on a 50@xmath172  mpc static grid (  unigrid " ) cube with @xmath173 cells , denoted as run 50_1024_2 in table  1 of smith et  al .
the radiative heating from the ionizing background plays an important role in determining the properties of the filamentary structure . as a filament
is ionized and heated , its density drops and its temperature rises ; both effects reduce the recombination rate . to study these effects
, we ran four moderate - resolution ( 50@xmath172  mpc unigrid cube with @xmath174 cells ) simulations with varying ionizing backgrounds , summarized in table  2 .
after initial submission of this paper , we ran a @xmath175 simulation , to check convergence and assess the  cosmic variance " among eight sub - volumes of the @xmath173 and @xmath175 simulations .
the standard uv background was taken from haardt & madau ( 2001 ) , although we also explore new computations of high-@xmath12 sfrs by trenti et  al .
( 2010 ) and haardt & madau ( 2012 ) .
these four simulations were : ( 1 ) no photoionizing background ; ( 2 ) uv background ramped up from @xmath176 to @xmath177 ( run 50_512_2 from smith et  al .
2011 ) ; ( 3 ) uv background ramped up from @xmath178 to @xmath179 ; and ( 4 ) uv background ramped up from @xmath178 to @xmath179 , but twice as strong as in ( 3 ) .
post - processing of the simulations was performed using the data analysis and visualization package , ` yt ` , documented by turk et  al .
( 2011 ) .    for regions of ionized hydrogen of density @xmath180
, we calculate the clumping factor , @xmath44 . using two different methods .
the first calculation , which has been used in some earlier studies , uses density weighting . in this
 density field " ( df ) method we define @xmath181 and average the density fields ( quantities @xmath182 , denoting either @xmath180 or @xmath183 ) where parameter averages are computed by summing over grid cells ( @xmath184 ) , subject to various `` cuts '' on the igm overdensity and gas temperature and weighted by factors , @xmath185 , @xmath186 in the second calculation , we compare the local recombination rate to the global average recombination rate , averaged over density and temperature in cells , @xmath187 where again @xmath188 denotes a weighted average and @xmath189 is the case - b radiative recombination rate coefficient for hydrogen , as tabulated by osterbrock & ferland ( 2006 ) .
we refer to this as the  recombination rate " ( rr ) method .    because the purpose of the clumping factor is to correct for an enhanced recombination rate , we believe @xmath190 to be a more appropriate representation .
the recombinations that are important to removing ionizing photons occur in the filamentary structure of the igm , and the clumping factor should only be calculated in these structures . to assess the critical sfr necessary to maintain an ionized medium , we focus on grid cells that are significantly ionized .
because a negligible amount of recombination occurs in cells containing mostly neutral gas , these cells should not contribute to the clumping factor .
if we do not exclude neutral gas , the clumping factor is extremely high ( @xmath191 ) before the ionizing background is turned on .
we find large density gradients between the neutral gas ( uniformly distributed over the simulation ) and the ionized gas . because the clumping factor is a measure of the inhomogeneity of the medium , a large density gradient yields a large value of @xmath44 .
if low - density voids are included in the calculation , the larger density gradient leads to an overestimate of the clumping . by setting both upper and lower density thresholds
, we can exclude collapsed halos and low - density voids from our calculations .
previous studies ( miralda - escud 2000 ; miralda - escud et  al .
2003 ; pawlik et  al .
2009 ) addressed these issues by setting a density threshold that excludes collapsed halos , but they did not set a lower limit to exclude voids from their calculations . to explore these effects in a filamentary igm
, we make various cuts of our data in baryon overdensity ( @xmath192 ) and in temperature , metallicity , and hydrogen ionization fraction ( @xmath193 ) . in our standard formulation ,
we include only those cells that meet the following criteria : @xmath194 , @xmath195  k , @xmath196 , and @xmath197 .
we believe this  data cut " adequately represents unenriched igm lying on an adiabat ( see figure  19 of smith et  al .
we also explore the clumping factor with no lower density threshold ( total range @xmath198 ) .
our results , presented in section 3.1 , show a small increase in @xmath199 from the wider range in densities by including low - density cells with @xmath200 . however , these low - density voids do not contribute substantially to the recombination rate .    in summary ,
our prescriptions for calculating the clumping factor yield a more physical representation of the enhanced recombination rate that is an important component to many reionization models . by not assuming a fully ionized medium and by specifically following @xmath180 , we are able to exclude denser neutral gas that does not contribute appreciably to the recombination rate .
we make simple assumptions on the reionization process and reionization history ( section 3.2 ) , turning on the ionizing radiation field at redshifts @xmath142 or @xmath201 and following the thermal history arising from photoelectric heating , radiative cooling , and pressure smoothing ( pawlik et  al .
the metal - line and molecular cooling , metal transport , and feedback included in our simulations also allow us to accurately represent the thermodynamics of the gas , which has been shown to have a significant effect on the evolution of the clumping factor .
future models will include discrete sources and radiative transfer , accounting for temperature increases arising from photo - heating with spectral hardening ( abel & haehnelt 1999 ) . because our current simulations employ a spatially constant ionizing radiation field , we anticipate carrying out these more realistic situations .
early studies of igm clumping adopted high values , @xmath202 at @xmath86 ( gnedin & ostriker 1997 ) . as discussed earlier , we believe these values are too high for the ionized igm filaments , which have expanded as a result of the heat deposited by lyc photons . the differences between observations and inferred critical sfr densities can largely be attributed to this high clumping factor ( sawicki & thompson 2006 ; bouwens et  al .
more recent studies have trended towards less clumping .
raievi & theuns ( 2011 ) argue that using a global clumping factor overestimates the recombination rate , and that local values should be used instead . in this study
, we calculate the clumping factor for a series of high - resolution cosmological simulations for ionized hydrogen and helium that explore how the photoheating of an ionizing background affects the igm thermodynamics ( density and temperature ) .
we impose criteria in overdensity , temperature , metallicity , and ionization fraction to constrain our calculations to igm filaments where recombinations and ionization fronts are most important .
the clumping in simulation @xmath203 is calculated as a function of redshift from @xmath204 to @xmath205 for both @xmath206 and @xmath190 , weighted by overdensity and by volume .
note that the unigrid simulations automatically yield volume weighting , when averaged over cells .
for the remainder of our simulations , we find a power - law overdensity distribution , with an average form @xmath207 .
weighting these cells by overdensity gives undue emphasis on small numbers of high - density cells .
this method is also biased , since it does not calculate the clumping factor where most of the recombinations are occurring .
figure 1 compares the df and rr methods and shows the difference between the different weights used in our calculations . in both methods ,
weighting by overdensity results in higher clumping factors .
the recombination - rate method ( eq .  [ 15 ] ) results in a lower clumping factor because it accounts for both density and temperature effects from the ionizing background on the clumping averages .
photoelectric heating during photoionization causes the filaments to expand , lowering the density .
the increased electron temperature also reduces the radiative recombination rate coefficient . the df method ( eq .  [ 13 ] ) includes no dependence on the recombination rate coefficient , and therefore only counts the effect of photoheating on the filament density .
see section  3.2 for further discussion .    in our current simulations
, we turn on a spatially constant uv background at @xmath142 or @xmath201 .
the effects of different  density cuts " in the summation over cells ( see section 2.3 ) are shown in figure 2 .
we find little difference at @xmath208 ( turn - on of ionizing radiation ) , and see a small increase in @xmath44 at @xmath209 when we include the lower density cells with @xmath200 . in future , more refined simulations with radiative transfer , we expect photoionization to be initiated in high - density regions , where the stars are located .
once the uv background is turned on , lower - density regions are easily photoionized .
the photoelectric heating of the uv background causes the filaments to expand and become diffuse , resulting in a lower clumping factor .
the elevated temperature also reduces the recombination rate coefficient , which in turn lowers @xmath44 when calculated by the rr method .
when calculating a mean recombination rate , only the recombining ionized gas is relevant .
by @xmath210 , the igm is nearly completely ionized and only small regions of  remain .
as noted earlier , in more refined simulations , the remaining  would be largely self - shielded from the uv background and likely to reside in clumps of high density yet unreached by ionization fronts .
figure 3 explores the convergence and cosmic variance among different simulations . in two panels ,
we compare results from our @xmath174 simulations with larger simulations with @xmath173 and @xmath175 cells . in each panel , we show the average clumping factors , together with those in eight sub - volumes .
the average values of @xmath44 agree in the @xmath173 and @xmath175 simulations , and the small ( @xmath211% ) variance among these sub - volumes suggests that the @xmath173 and @xmath175 simulations are converged . over the redshift range @xmath212 in the @xmath175 simulation , the clumping factor is well fitted by a power - law , @xmath213 ^{-1.1 }    \ ;   , \ ] ] representing a slow rise in clumping to lower redshift , after the turn - on of ionizing radiation .      to study the effect of radiative heating from the ionizing background on filamentary structure , we ran a suite of moderate - resolution simulations with different ionizing backgrounds ( see table  2 for details ) .
the redshift at which the uv background is turned on affects the clumping factor , causing it to drop and then recover at lower redshift .
pawlik et  al .
( 2009 ) attribute this effect to _ jeans filtering _ , where the photo - heating raises the cosmological jeans mass , preventing further accretion onto low - mass halos and smoothing out small - scale density fluctuations .
the photoheating also heats the filaments we are focusing on , which lowers the hydrogen recombination rate coefficient , @xmath214 , and causes the filaments to expand ; both of these effects reduce the clumping factor . without a uv background , the clumping continues to rise as the filaments gravitationally collapse .
pawlik et  al .
( 2009 ) also claim that the clumping factor at @xmath177 is insensitive to when the background is turned on , as long as it is turned on at @xmath215 . when an ionizing background is introduced , photoheating acts as a positive feedback to reionization by lowering the clumping factor and making it easier to stay ionized .
this same photoheating mechanism also suppresses star formation in low - mass halos , which in turn lowers the ionizing photon production rate by star - forming galaxies and acts as a negative feedback to reionization ( pawlik et  al .  2009 ) .
these thermodynamic processes affect clumping and structure and emphasize the importance of carefully modeling the strength of feedback processes and their effects on jeans mass .
we have explored what happens when the ionizing radiation turned on earlier , especially during the interval @xmath3 marking the transition from a neutral to fully ionized igm . in a series of @xmath174 simulations ,
we explore the influence of turn - on of photoionizing radiation between redshifts @xmath142 and @xmath201 .
figure  4 shows @xmath44 , computed for the moderate - resolution simulations via the density - field method and weighted by volume .
we do not plot simulation @xmath216 , since it is identical to simulation @xmath217 .
as in the results of pawlik et  al .
( 2009 ) , we find that photoheating from the ionizing background results in a decrease in the clumping factor .
the clumping factors fall along two tracks : a higher track at redshifts above the turn - on of the uv background , and a lower one after turn - on .
by @xmath141 , we find nearly identical clumping factors for both backgrounds .
after a substantial recovery time , the redshift when the background is turned on is not important . before this recovery ,
the clumping factor of the earlier background is lower by a factor of @xmath782 , resulting in earlier reionization .
one can compare the strengths of feedback to reionization , where  positive feedback " lowers the clumping factor and  negative feedback " suppresses star formation .
at @xmath218 the total stellar mass ( sfr density ) of simulation @xmath219 is 1.19 ( 1.13 ) times lower than that of simulation @xmath220 , while the clumping factor is 1.64 times lower .
for simulation @xmath217 , the total stellar mass ( sfr density ) and clumping factor are 1.52 ( 1.49 ) and 1.66 , respectively , times lower than those of simulation @xmath220 .
this suggests that the positive feedback introduced by the background is greater than the negative feedback .
however , the stellar mass ( sfr density ) does not recover in the same manner as the clumping factor .
once photoheating suppresses the formation of small - mass halos , the hubble flow takes over and prevents them from collapsing and forming stars .
therefore , the redshift at which the background is turned on determines whether the positive or negative feedback dominates and whether the ionizing background will cause reionization to be accelerated or delayed .      in connection with this project ,
we have developed a user interface for calculating the critical sfr density @xmath221 needed to maintain the igm ionization at a given redshift .
the software computes the effects of variations in the stellar imf ( slope , mass - range ) and model atmospheres , and the redshift evolution of metallicity and gas thermodynamics ( density , temperature , coupling to the cmb ) .
the clumping factor and lyc escape fraction are free parameters in the calculator .
our simulations find ranges of @xmath222 depending on redshift , overdensity , and thermal phase of the ( photoionized or shock - heated ) igm .
for the new @xmath175 simuation , the global mean clumping factor is @xmath223 , with a power - law fit @xmath224^{-1.1}$ ] for redshifts between @xmath212 .
this calculator is a useful tool for determining the population of galaxies responsible for reionization .
the critical sfr per co - moving volume ( eq . [ 9 ] )
is obtained by balancing the production rate of lyc photons with the number of hydrogen recombinations . here
, @xmath225 is the clumping factor of ionized hydrogen , @xmath11 is the escape fraction of lyc photons from their host galaxies , and @xmath226 is the temperature scaled to @xmath227  k. the conversion factor , @xmath92 , from stellar mass to total number of lyc photons produced , is regulated by the imf and model atmospheres .
the user has the option of controlling these parameters to determine the resulting critical sfr density , subject to several observational constraints .
this simulator can be accessed at http://casa.colorado.edu/~harnessa/sfrcalculator with an ` html ` interface for easy use .
two standard tests of the sfr use the simulator to calculate the ionization histories of  and  and compare them to the ionized volume filling factor , @xmath228 , and the cmb optical depth , @xmath51 .
the calculator derives @xmath51 by integrating the differential form of eq .
( 3 ) for various sfr histories and parameters .
the average evolution of @xmath229 is found by numerical integration of the rate equation ( madau , haardt , & rees 1999 ) expressing the sources and sinks of ionized zones , @xmath230 the source term , @xmath231 , represents the net ionizing photon production rate , computed from selected models of the sfr density , @xmath1 . here , @xmath232 is the mean hydrogen number density , and @xmath233 is the hydrogen recombination timescale , which depends on @xmath44 as shown in eq .
the igm is assumed to be fully ionized when @xmath234 .
we integrate a similar equation for @xmath235 to follow  photoionization in the qso 4-ryd continuum .
we use the qso emissivities at 1 ryd ( haardt & madau 2012 ) , extrapolated to 4 ryd assuming a spectrum with specific flux @xmath236 .
figure  5 shows ionization histories , quantified by @xmath228 and @xmath51 , for various values of clumping factor , escape fraction , igm temperature , and lyc - production efficiencies ( @xmath92 ) .
we adopt the sfr history from trenti et  al .
( 2010 ) with an evolving luminosity function .
our standard model adopts @xmath6 , @xmath237 , @xmath238 , and @xmath131 ( @xmath239 photons / m@xmath240 ) .
different curves show the effects of changing @xmath44 and @xmath11 , including two models in which these parameters evolve with redshift . with these parameters , we are typically able to complete reionization by @xmath43 , consistent with observations of gunn - peterson troughs , redshift evolution in ly@xmath2  emitters , and igm neutral fraction
. the model with a constant @xmath241 does not complete the ionization until @xmath242 , which is far too late .
thus , we be believe that lyc escape fractions must be considerably larger ( @xmath243% ) perhaps evolving to higher values at @xmath24 shown by green and magenta curves .
figure 6 compares the dependence of @xmath228 and @xmath51 on sfr histories , computed with an evolving luminosity function ( trenti et  al .
2010 ) and new calculations ( haardt & madau 2012 ) of the star formation rate density , @xmath244 .
figure 7 compares the effects of two choices of model atmospheres , with a fixed sfr history from trenti et  al .
the on - line calculator provides ionization fractions , @xmath228 and @xmath245 , together with cmb optical depth , @xmath51 . in the simulator
, users can select other igm parameters and sfr histories .
the main difference between the two sfr models lies in the assumptions at @xmath246 .
haardt & madau ( 2012 ) rely on an empirical extrapolation of the star formation rate as a function of redshift , while trenti et  al .
( 2010 ) adopt a physically motivated model based on the evolution of the dark - matter halo mass function .
the two approaches are similar at @xmath247 , but differ significantly at higher redshift , where an empirical extrapolation does not capture the sharp drop in the number density of galaxies observed at @xmath248 ( see figure 8 in oesch et  al .
2011 ) . as a consequence ,
the reionization history from the haardt & madau ( 2012 ) model is more extended at high @xmath12 , especially when the efficiency of reionization is increased because of evolving clumping factor and escape fraction ( our preferred models , green and magenta curves in fig .
the two models yield quite different predictions for the duration of reionization , defined as the redshift interval , @xmath249 , over which @xmath229 evolves from 20% to 80% ionized .
our preferred sfr models ( green and magenta lines in fig .  5 ) have @xmath19 ( from @xmath250 to @xmath251 ) , whereas the haardt - madau sfr histories exhibit a more extended interval , @xmath252 ( from @xmath253 to @xmath43 ) .
this difference in @xmath249 could be tested by upcoming 21-cm experiments ; see bowman & rodgers ( 2010 ) for an initial constraint , @xmath254 .
figure 8 illustrates a third constraint on the reionization epoch ( pritchard et  al .
2010 ; lidz et  al .
2011 ) comparing sfr histories with estimates of the ionizing background at @xmath255 ( fan et  al .
2006 ; bolton & haehnelt 2007 ; haardt & madau 2012 ) . the lyc co - moving emissivity ( in photons s@xmath147
mpc@xmath48 ) , defined as @xmath231 , can be related to the ionizing background at @xmath256 , using recent estimates of the hydrogen photoionization rate , @xmath257 , from haardt & madau ( 2012 ) and the lyc mean free path , @xmath258 , from songaila & cowie ( 2010 ) .
the lyc emissivity is proportional to the star formation rate density , computed from our halo mass function model ( trenti et  al .
2010 ) , integrated down to absolute magnitudes @xmath259 ( bouwens et  al .
2011a ) or to @xmath260 , the faint limit suggested by trenti et  al .
we adopt a fiducial lyc production parameter @xmath131 , corresponding to @xmath239 lyc photons produced per @xmath96 of star formation , and we use two different models for lyc escape fraction , @xmath11 ( constant at 20% and varying with redshift ) . for quantitative values
, we assume an ionizing background with specific intensity , @xmath261 ( in units erg  @xmath262  s@xmath147  sr@xmath147  hz@xmath147 ) with power - law index @xmath263 at energies above @xmath264 ryd . the hydrogen photoionization rate is @xmath265 $ ] for a hydrogen photoionization cross section @xmath266 with @xmath267  @xmath268 . for this spectrum ,
the frequency - integrated ionizing intensity is @xmath269 , and the lyc photon flux ( photons  @xmath262  s@xmath147 ) integrated over all solid angles is @xmath270 .
because we analyze the _ photon _ mean - free path , we normalize @xmath271 to @xmath272 and @xmath273 , @xmath274 by approximating @xmath272 as the product of the lyc emissivity and mean - free path , we arrive at the calibrations : @xmath275      \\
\phi_{\rm lyc } & = &   \frac { ( \alpha+3)}{\alpha }   \frac { \gamma_{\rm hi } } { \sigma_0 }            = ( 1.98 \times 10 ^ 5~{\rm photons~cm}^{-2}~{\rm s}^{-1 } )              \left [ \frac { \gamma_{\rm hi } } { 5 \times 10^{-13}~{\rm s}^{-1 } } \right ]    \ ;   .
\end{aligned}\ ] ] finally , we relate @xmath272 to the star - formation rate @xmath1 and ionization rate @xmath273 , @xmath276 \left [ \frac { 0.1 }   { f_{\rm esc } } \right ]                \left [ \frac   { 0.004}{q_{\rm lyc } } \right ] \left [ \frac { 9.8~{\rm pmpc } } { \lambda_{\rm hi } } \right ]                \left [ \frac { 6.5}{1+z } \right]^{3 }        \ ; .\ ] ] here , we have scaled the parameters to the same values assumed by lidz et  al .
( 2011 ) , namely @xmath277 , @xmath278 , an ionizing spectral index @xmath279 , and redshift @xmath280 , at which songaila & cowie ( 2010 ) fit @xmath281 proper mpc .
our parameter @xmath131 corresponds to the lidz et  al .
( 2011 ) lyc photon production calibration , @xmath282 per @xmath47 of star formation .
however , our coefficient , @xmath283 , is slightly larger than their value , @xmath284 , an effect that may arise from our different method of relating @xmath273 to the ionizing radiation field and sfr .    more careful examination of these parameters suggests that , at @xmath285 , the hydrogen ionization rate @xmath286 , given in table 3 of haardt & madau ( 2012 ) , and the observed mean free path , @xmath258 , at @xmath280 may be closer to 6 proper mpc ( see figure 10 of songaila & cowie 2010 ) .
rescaling to those two parameters , we find a similar coefficient of @xmath287 .
this sfr density at @xmath255 ( see also figure 8) is comparable to the critical value needed to maintain reionization at @xmath142 , but the observed rates and mean free paths are declining rapidly with redshift .
all three constraints on sfr suggest that full reionization is more likely to occur at redshift @xmath0 than at @xmath288 .
( uv background turned on at @xmath176 ) with different weights .
red lines : clumping calculated via density field ( df ) method .
blue lines : clumping via recombination rate ( rr ) method ( see section 2.2 ) .
solid lines correspond to weighting by volume , and dashed lines correspond to weighting by baryon overdensity .
we believe the rr method , with volume weighting , is a more accurate measure of clumping . ]
simulation , using the rr method in equation ( 15 ) .
solid line : clumping factor summed over cells with overdensity between @xmath289 , excluding low - density voids with @xmath200 . dashed line : including all cells with @xmath198 .
]     agree for @xmath173 and @xmath175 simulations , with variance computed over eight sub - volumes in our simulations .
top panel : eight @xmath290 sub - volumes in @xmath175 run .
bottom panel : eight @xmath174 sub - volumes in @xmath173 run . in each panel , our original @xmath174 run is shown as dotted black line , and values for @xmath175 and @xmath173 runs as heavy black lines .
large variations in @xmath44 at high redshifts arise from small numbers of cells prior to turn - on of the ionizing background at @xmath201 ( top ) and @xmath142 ( bottom ) . for the @xmath175 simulation ( top panel )
, the clumping factor is well - fitted by @xmath224^{-1.1}$ ] between @xmath212 .
, title="fig : " ]   agree for @xmath173 and @xmath175 simulations , with variance computed over eight sub - volumes in our simulations .
top panel : eight @xmath290 sub - volumes in @xmath175 run .
bottom panel : eight @xmath174 sub - volumes in @xmath173 run . in each panel
, our original @xmath174 run is shown as dotted black line , and values for @xmath175 and @xmath173 runs as heavy black lines .
large variations in @xmath44 at high redshifts arise from small numbers of cells prior to turn - on of the ionizing background at @xmath201 ( top ) and @xmath142 ( bottom ) . for the @xmath175 simulation ( top panel ) ,
the clumping factor is well - fitted by @xmath224^{-1.1}$ ] between @xmath212 .
, title="fig : " ]     igm simulations , with sums weighted by volume ( see eq . [ 14 ] ) .
red solid line has uv background ( haardt & madau 2001 ) turned on at @xmath178 ( simulation @xmath217 ) , blue dashed line at @xmath176 ( simulation @xmath219 ) , and green dot - dashed line has no background ( simulation @xmath220 ) .
we observe two tracks for @xmath199 : a high - track for redshifts above turn - on of the uv background , and a lower track after turn - on . during most of the reionization epoch , from @xmath21
, @xmath44 lies between 1.5 and 3 . ]     and cmb optical depth @xmath51 versus redshift ( eqs .   [ 3 ] and [ 17 ] ) computed with @xmath291 k , @xmath131 ( @xmath132 photons per @xmath96 ) and using sfr history and evolving luminosity function from trenti et  al .
( 2010 ) integrated down to @xmath260 .
the igm is assumed to be fully ionized when @xmath234 .
blue line : constant @xmath6 and @xmath241 ( this model reionizes too late to be viable ) .
red line : constant @xmath6 and @xmath7 .
green line : constant @xmath6 but redshift - dependent @xmath292 ( haardt & madau 2012 ) .
magenta line : redshift - dependent @xmath293 ( pawlik et  al .
2009 ) and @xmath294 . our two preferred models with variable @xmath44 or @xmath11 ( green , magenta ) naturally produce @xmath0 .
solid black circles indicate redshifts of 20% , 50% , and 80% ionization ; the duration @xmath19 is defined between 20% and 80% points . ]    , with same color codes as fig .  5 , comparing two models for sfr history .
solid lines show sfr history from trenti et  al .
( 2010 ) , and dashed lines show sfr from equation ( 53 ) in haardt & madau ( 2012 ) , which is larger and more extended at higher redshifts , particularly for @xmath295 . bottom : corresponding cmb optical depths @xmath51 .
additional optical depth may arise from sources at @xmath17 .
, title="fig : " ] , with same color codes as fig .  5 , comparing two models for sfr history .
solid lines show sfr history from trenti et  al .
( 2010 ) , and dashed lines show sfr from equation ( 53 ) in haardt & madau ( 2012 ) , which is larger and more extended at higher redshifts , particularly for @xmath295 . bottom : corresponding cmb optical depths @xmath51 .
additional optical depth may arise from sources at @xmath17 .
, title="fig : " ]    , @xmath245 , and @xmath51 from our on - line sfr / reionization simulator .
solid lines show model atmospheres labeled as sutherland & shull ( s&s ) , while dashed lines are for schaerer ( 2002 , 2003 ) , assuming @xmath296 .
models assume sfr history of trenti et  al .
( 2010 ) and parameters with constant @xmath6 and redshift - dependent @xmath292 ( haardt & madau 2012 ) .
dot - dashed line shows  ionization history , @xmath245 , computed for qso ( 4-ryd continuum ) ionization as described in section 3.3 . in these models , is fully ionized by @xmath43 and  by @xmath297 .
additional optical depth may arise from sources at @xmath17 . ]    ) from halo mass function model ( trenti et  al .
2010 ) integrated for galaxies down to @xmath259 ( solid blue line ) and to @xmath298 ( dotted green line ) .
red data points are from bouwens et  al .
( 2011a ) for galaxies down to @xmath259 .
solid data point at @xmath285 is from ionizing background constraint ( eq .
21 ) , using haaradt & madau ( 2012 ) model of @xmath299 and lyc mean free paths from ( songaila & cowie 2010 ) .
bottom : co - moving emissivity , @xmath300 , of lyc photons , corresponding to above sfrs , integrated to @xmath259 and @xmath298 .
we use a lyc production parameter @xmath131 ( @xmath301 lyc photons  s@xmath147 per @xmath96  yr@xmath147 ) and two different escape - fraction models : @xmath7 ( blue curves and data points ) and @xmath302 ( red curves and data points ) . , title="fig : " ] ) from halo mass function model ( trenti et  al .
2010 ) integrated for galaxies down to @xmath259 ( solid blue line ) and to @xmath298 ( dotted green line ) .
red data points are from bouwens et  al .
( 2011a ) for galaxies down to @xmath259 .
solid data point at @xmath285 is from ionizing background constraint ( eq .  21 ) , using haaradt & madau ( 2012 ) model of @xmath299 and lyc mean free paths from ( songaila & cowie 2010 ) .
bottom : co - moving emissivity , @xmath300 , of lyc photons , corresponding to above sfrs , integrated to @xmath259 and @xmath298 .
we use a lyc production parameter @xmath131 ( @xmath301 lyc photons  s@xmath147 per @xmath96  yr@xmath147 ) and two different escape - fraction models : @xmath7 ( blue curves and data points ) and @xmath302 ( red curves and data points ) . , title="fig : " ]
the major results of our study can be summarized with the following points :    1 .
we calculated the critical star formation rate required to maintain a photoionized igm , incorporating four free parameters ( @xmath44 , @xmath11 , @xmath45 , @xmath303 ) that control the rates of lyc photon production and radiative recombination .
our best estimate at @xmath142 is @xmath304 ^ 3   ( c_h/3)(0.2/f_{\rm esc } ) t_4^{-0.845}$ ] for fiducial values of igm clumping factor @xmath305 , lyc escape fraction @xmath7 , temperature @xmath8  k , and standard imfs and low - metallicity stellar atmosphere ( @xmath131 ) .
an epoch of full reionization at @xmath0 is consistent with recent optical / ir measurements of sfr history and a rising igm neutral fraction at @xmath3 , marking the tail end of reionization .
these observations include the decrease in numbers of high-@xmath12 galaxies and ly@xmath2  emitters and igm damping - wing intrusion into the ly@xmath2  transmission profiles in high-@xmath12 qsos .
our newly developed sfr and reionization calculator , now available on - line at + http://casa.colorado.edu/~harnessa/sfrcalculator , allows users to calculate the ionization history , critical sfr , and cmb optical depth , and to assess whether observed sfrs , imfs , and other parameters are consistent with igm reionization .
4 .   reconciling late reionization at @xmath0 with @xmath13 of the cmb
likely requires an epoch of partial ionization .
a fully ionized igm back to @xmath79 produces @xmath306 , and an additional optical depth , @xmath16 , may arise from early sources of uv / x - ray photons at @xmath17 .
alternatively , one can appeal to the likelihood contours from wmap , which allow optical depths as low as @xmath307 ( 95% c.l . )
the _ planck _ experiment may clarify the situation in several years .
if the eor is more complex , as we suggest , then redshifted 21-cm experiments should focus on the interval @xmath308 , corresponding to frequencies 124167 mhz , for the maximum signal of igm heating ( @xmath309  mk ) produced when the hydrogen neutral fraction @xmath53 ( pritchard et  al .
2010 ) .
we conclude by speculating about which future observations can best constrain the eor .
additional data from from the _ planck _ mission should provide confirmation of the cmb optical depth with smaller error bars .
this information will constrain the additional amount of ionization at @xmath17 and narrow the range @xmath310 produced by high - redshift sources in the partially ionized igm . ongoing surveys for high-@xmath12 galaxies , ly@xmath2  emitters , and qso near - zone sizes will better quantify the rise of neutral fraction , @xmath311 at @xmath312 . on the theoretical front
, we are running larger igm simulations on @xmath175 and @xmath313 unigrids and will add discrete sources of ionizing radiation and radiative transfer in order to capture the heating and clumping more accurately .
we will also include source turn - on at @xmath208 to test the decrease and recovery of @xmath44 , as seen in figures 3 and 4 .
finally , we plan to carry out more detailed modeling of the  ( 21-cm ) signal , coupled to the kinetic and spin temperatures driven by heating at @xmath17 .
as described earlier , the duration of the reionization transition and the 21-cm emission during the heating phase could provide discriminants of various sfr histories at @xmath24 .
this work was supported by grants to the astrophysical theory program ( nnx07-ag77 g from nasa and ast07 - 07474 from nsf ) at the university of colorado boulder .
we thank joanna dunkley , chris carilli , richard ellis , and piero madau for useful discussions on reionization and cmb optical depth statistics .
we are grateful to the referee for suggesting additional model comparisons with the ionizing background at @xmath256 .
cccccccc 0.1 & 100 & 2.35 & @xmath314 & 0.00236 & 0.00286 & 0.0306 & 0.0253 + 0.1 & 100 & 2.35 & @xmath315 & 0.00397 & 0.00558 & 0.0181 & 0.0129 + 0.1 & 100 & 2.35 & @xmath316 & 0.00401 & 0.00752 & 0.0180 & 0.0096 + 0.1 & 100 & 2.35 & @xmath317 & 0.00365 & @xmath318 & 0.0197 & @xmath318
+ 0.1 & 100 & 2.35 & @xmath319 & 0.00383 & @xmath318 & 0.0188 & @xmath318 + 1.0 & 100 & 2.35 & @xmath319 & 0.00976 & @xmath318 & 0.0074 & @xmath318 + 0.1 & 100 & 2.00 & @xmath319 & 0.01267 & @xmath318 & 0.0057 & @xmath318 + 0.1 & 200 & 2.35 & @xmath319 & 0.00543 & @xmath318 & 0.0133 & @xmath318 +    lcccc [ tbl : sims ] @xmath203 & 50 & 1024 & 7 & 1 + @xmath219 & 50 & 512 & 7 & 1 + @xmath217 & 50 & 512 & 9 & 1 + @xmath216 & 50 & 512 & 9 & 2 + @xmath220 & 50 & 512 & n / a & 0 + @xmath320 & 50 & 1536 & 9 & 1 +