Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: random tessellations * st*able with respect to * it*erations ( stit ) have recently been introduced in @xcite and @xcite as a new model for random tessellation in @xmath3 and have quickly attracted considerable interest in modern stochastic geometry as well fitting the growing demand for non - trivial and flexible but mathematically tractable tessellation models . the stit tessellations may be interpreted as outcome of a random cell division process , which makes them very attractive for applications , see for example @xcite . other potential applications include mathematical modelling of systems of cracks , joints or fissures in rock or the so - called craquele of thin layers . a general approach to random cell division processes has recently appeared in @xcite and the construction of stit tessellations can roughly be described as follows . at first , we fix a compact and convex window @xmath1 in which the construction is carried out . for simplicity we assume @xmath4 to be a polytope , as this implies that the resulting tessellation has polytopal cells with probability one . next , an exponentially random life time is assigned to @xmath4 , whereby the parameter of the distribution is given as a constant multiple of the integral - geometric mean width of @xmath4 , see section [ secstit ] for details . upon expiry of this life time a @xmath4 hitting uniform random hyperplane is chosen , is introduced in @xmath4 and divides the window into two polytopal sub - cells @xmath5 and @xmath6 . the construction continues now recursively and independently in both of these sub - cells , where the newly introduced hyperplanes are always chopped - off by the boundaries of their mother - cells . the whole construction is continued until some deterministic time threshold @xmath0 is reached . regarded in time , the construction can by interpreted as a pure - jump markov process on the space of tessellations of the window @xmath4 . + our construction shares some common features with random fragmentation processes or branching markov chains in the sense of @xcite . the cells of the resulting tessellation within @xmath4 can be regarded as particles in a suitable polish space and the dynamics of the particles is non - interacting in a sense that different particles ( cells of the tessellation ) have independent evolutions , which is indeed the case in our construction . moreover , whenever a particle dies it is replaced by exactly two new particles , namely the two newly generated sub - cells . in addition , the life times of the particles are exponentially distributed , as assumed in @xcite , but in general not independent in contrast to the fragmentation theory . + in the recent paper @xcite the authors have introduced a new technique for studying the geometric properties of stit tessellations based on martingales and the general theory of martingale problems for pure - jump type markov processes . in particular , with these new developments , the variance of the total surface area of a stationary and isotropic iteration stable random tessellation @xmath7 in @xmath1 ( a random stit tessellation ) has been determined by integral - geometric means and the corresponding central limit theory has been established . strikingly , as already signalled by results in the special two dimensional case @xcite , and as confirmed by the results of the present paper , it turns out that in a certain rather strong sense the asymptotic behaviour of the surface area process dominates and fully determines the asymptotic geometry of the stit tessellation . in particular the surface area variance is _ the _ basic second - order parameter of the tessellation and all second - order characteristics of functionals considered in our work can be reduced to it , likewise non - trivial functional limits in law exist for the stit surface area process whereas the limits in law for other natural related processes arise as deterministic functionals of the corresponding surface process . another crucial phenomenon arising in this context , as first noted in @xcite , is that the asymptotic theories for dimension @xmath8 and @xmath9 differ strongly in many important aspects . we have studied the planar case @xmath8 in the recent separate paper @xcite and thus only consider @xmath9 in the present one . + the purpose of this paper is to establish a second - order and limit theory for integral - geometric characteristics of stationary and isotropic stit tessellations in dimensions higher than two . more precisely , the characteristics studied in this work are the cumulative intrinsic volumes of all orders for the collections of so - called maximal polytopes of @xmath7 . these are in codimension @xmath10 the cell - separating facets introduced during the random cell division process described above , constituting the basic building blocks of a stit tessellation as discussed in detail below . we shall provide explicit as well as asymptotic variance expressions for these parameters of the random tessellation @xmath11 as @xmath12 tends to infinity for a sequence @xmath13 of expanding convex windows . further , we will find the covariance measures for random lower - dimensional face measures generated by @xmath14 finally , we shall also give the corresponding convergence in law statements , obtaining non - gaussian limits for the studied case @xmath15 as opposed to the classical gaussian limits arising for @xmath8 , see @xcite . +   + the paper is organized as follows : in section [ secstit ] we define stit tessellations , specializing to the stationary and isotropic set - up in the focus of this paper , and we discuss their basic properties that are needed in our arguments . next , in section [ sec2ndord ] we calculate the variances of the cumulative intrinsic volumes of all orders for stit tessellations . this includes both , exact formulae and asymptotic analysis upon letting the window size grow to infinity . further , in section [ sec2ndmeas ] we extend the second - order analysis to the level of lower - dimensional face measures induced by stit tessellations , thus taking into account not only the numeric characteristics but also the spatial profile of the stit face processes . finally , in section [ secclt ] we develop the corresponding functional limit theory with non - gaussian limit processes for @xmath16 in order to keep the paper self - contained , we will recall important faces from @xcite and sketch some of their proofs for the readers convenience . the purpose of this section is to provide a short self - contained discussion of stit tessellations in @xmath3 as studied in this paper and to summarize their basic properties for easy reference . we will restrict to the stationary and isotropic case in subsection [ subsecmnwc ] below , specializing to the scope of the paper . the general reference for these  by now classical  properties throughout this section is @xcite . next , in subsection [ subsecmart ] we will discuss the martingale tools developed in @xcite and underlying our present theory . finally , in subsection [ subsec1stord ] we provide certain useful mean value relationships for intrinsic volumes in context of stit tessellations . we start with a compact and convex polytope @xmath1 in which our construction is carried out and denote by @xmath17 the standard isometry - invariant measure on the space @xmath18 of ( affine ) hyperplanes in @xmath3 normalized so as to induce unit surface intensity on @xmath19 ( note that in our earlier papers @xcite and @xcite this measure has been denoted by @xmath20 ) . we call @xmath17 the _ driving measure _ of the construction . assign now to @xmath4 an exponentially random lifetime with parameter @xmath21)$ ] , where @xmath22:=\{h\in{\cal h}:\ h\cap w\neq\emptyset\}\ ] ] is the set of hyperplanes hitting @xmath4 . upon expiry of this random life time , a random hyperplane is chosen according to the distribution @xmath21)^{-1}\lambda(\cdot\cap[w])$ ] , is introduced in @xmath4 and is chopped off by its boundary . this is , the window @xmath4 splits into the two polyhedral sub - cells @xmath5 and @xmath6 that are separated by the introduced hyperplane piece . the construction continues now recursively and independently in @xmath5 and @xmath6 and is stopped if some previously fixed deterministic time threshold @xmath23 is reached . our assumptions ensure that the cells of the tessellation constructed until time @xmath23 are convex polyhedra in @xmath4 with probability one . they are denoted by @xmath24 and we denote by @xmath7 the random closed set in @xmath4 that consists of the union of cell - boundaries of cells constructed until time @xmath0 , see figure [ figure ] . the construction of @xmath7 is referred to as the _ mnw - construction _ after the names of its inventors  mecke , nagel and weiss  and the tessellation @xmath7 itself is called a random _ stit tessellation_. the abbreviation stit comes from the crucial property enjoyed by @xmath25 namely that of being able under eration , for which we refer to @xcite or @xcite . we call the cell separating @xmath26-dimensional faces the _ @xmath26-dimensional maximal polytopes _ and denote the collection of all such polytopes of @xmath7 by @xmath27 . moreover , we introduce the set @xmath28 of @xmath29-dimensional faces of @xmath26-dimensional maximal polytopes of @xmath7 by @xmath30 where by @xmath31 we mean the set of all @xmath29-dimensional faces of the @xmath26-dimensional polytope @xmath32 , @xmath33 . the maximal polytopes are often referred to as _ i - polytopes _  this terminology originates from the historically first considered particular case @xmath8 where the maximal polytopes ( which are just line segments in this case ) assumed shapes similar to the literal @xmath34 +         it is an important observation that the spatio - temporal construction of @xmath7 satisfies the markov property in time parameter @xmath0 , which means that @xmath35 where @xmath36 denotes the operation of iteration of tessellations , see @xcite , @xcite or @xcite . to make the paper more self - contained , we recall now the definition of the meaning of the operation @xmath36 . to this end , we regard the tessellation @xmath7 as a frame or primary tessellation and associate with each cell @xmath37 an i.i.d . copy @xmath38 of @xmath39 . another tessellation @xmath40 of @xmath4 is now defined as @xmath41 i.e. we consider the union of @xmath7 and the cut - outs of @xmath38 s within the cells @xmath42 of the primary tessellation ( note that in the definition we have used the interpretation of a tessellation as a random closed subset of @xmath4 ) . we say that @xmath40 is the _ iteration _ of @xmath7 with @xmath39 . the remarkable property of our tessellations constructed by the mnw - process is that the outcome @xmath43 coincides in law with @xmath44 , i.e. with the continuation of the mnw - construction until time @xmath45 . thus , it is the same either to continue the mnw - construction from @xmath0 until time @xmath45 or to perform at time @xmath0 an iteration of @xmath7 with @xmath39 . + the local properties established so far , can be extended to the whole space , since the random tessellations @xmath7 are _ consistent _ in @xmath4 , by which we mean that for any fixed @xmath23 and @xmath46 it holds @xmath47 where @xmath48 stands for equality in distribution . this implies  in view of the consistency theorem ( * ? ? ? 2.3.1 )  the existence of the whole - space tessellation @xmath49 such that @xmath50 for each compact convex @xmath1 . + it directly follows from the markov property that the random tessellations @xmath49 are stable with respect to the operation @xmath36 , i.e. @xmath51 where @xmath52 means the dilation with a factor @xmath53 . this property explains also the abbreviation stit , because the last equation is a classical probabilistic stability relation . the random tessellations @xmath49 share another important property , namely that the _ intersection _ of @xmath49 with a @xmath29-dimensional plane @xmath54 , @xmath55 , is again an iteration stable random tessellation . more precisely , we have @xmath56 and the integral - geometric constant @xmath57 is given by @xmath58 see eq . ( 3.29 t ) in @xcite . moreover , it is easy to see from the properties of the capacity functional of @xmath49 , compare with ( * ? ? ? 5(ii ) ) , that stit tessellations have the following _ scaling property _ : @xmath59 i.e. the tessellation @xmath49 has , upon rescaling with the factor @xmath60 the same distribution as @xmath61 , that is the stit tessellation with surface intensity @xmath10 . + we close this section by mentioning that the random tessellations @xmath49 have another important property , namely that @xmath49 has _ poisson typical cells _ or poisson cells for short , see lemma 3 in @xcite or part ( 1 ) of proposition [ propxx ] below . this is to say , the typical cell of @xmath49 has the same distribution as the typical cell of a stationary and isotropic poisson hyperplane tessellation with surface density @xmath0 ( see chap . 10.3 in @xcite for references about the classical poisson hyperplane model ) . in particular , this fact combined with the intersection property from above shows that the typical cell of the lower dimensional stit tessellation @xmath62 , with @xmath54 as above , has the same distribution as the typical cell of a poisson hyperplane tessellation in @xmath63 with surface density @xmath64 . more formally , we define the cell intensity measure @xmath65 of the stit tessellation @xmath7 by @xmath66 and its @xmath26-dimensional maximal polytope intensity measure @xmath67 by @xmath68 where @xmath69 stands for the uni mass dirac measure concentrated at @xmath70 . moreover , we let @xmath71 be the cell intensity measure and @xmath72 be the @xmath26-face intensity measure of a poisson hyperplane tessellation @xmath73 within @xmath1 having intensity measure @xmath74 , which are defined similarly to @xmath65 and @xmath67 above ( clearly , the @xmath26-dimensional maximal polytopes have to be replaced by the set of @xmath26-faces of @xmath73 ) . these definitions bring us in the position to reformulate special cases of theorems 1 and 2 from @xcite adapted to our later purposes :    [ propxx ] we have @xmath75      as already noted above , the mnw - construction of iteration stable random tessellations @xmath7 in finite volumes @xmath1 enjoys a markov property in the continuous time parameter @xmath76 in our previous work @xcite we have used this fact combined with the classical theory of martingale problems for pure jump markov processes to construct a class of natural martingales associated to the mnw - process . in this paper we only need a part of that theory . to formulate it we let @xmath77 be a bounded and measurable functional on the space of @xmath26-polytopes in @xmath3 and we denote by @xmath78 the collection of @xmath26-dimensional maximal polytopes of a tessellation @xmath79 standing for a generic realization of @xmath7 for some @xmath80 write @xmath81 any hyperplane @xmath82 $ ] hitting the window @xmath4 is tessellated by the intersection with @xmath79 and we denote by @xmath83 the set of all @xmath26-dimensional cells of this tessellation and introduce @xmath84}\sum_{f\in{\operatorname{cells}}(y\cap h)}\phi(f)\lambda(dh).\ ] ] it is also convenient to introduce the bar notation for centered versions of these quantities with @xmath85 that is to say @xmath86 and likewise @xmath87 with this notation , in view of the results developed in @xcite , we have    [ martprop ] for bounded and measurable functionals @xmath77 and @xmath88 on the space of @xmath26-polytopes in @xmath3 , the stochastic processes @xmath89 and @xmath90ds\label{mart2}\end{aligned}\ ] ] are martingales with respect to the filtration @xmath91 induced by @xmath92 .    [ [ sketch - of - a - proof ] ] sketch of a proof : + + + + + + + + + + + + + + + + + +    in order to make the paper self - contained , we give here the main idea of the proof of proposition [ martprop ] , although there is some overlap with @xcite . + at first , it is a direct consequence of the mnw - construction that the generator @xmath93 of the pure - jump markov process @xmath7 , @xmath23 , is given by @xmath94}\sum_{f\in{\operatorname{cells}}(y\cap h)}[f(y\cup f)-f(y)]\lambda(dh),\ ] ] where @xmath79 stands for some instant of @xmath7 and @xmath95 is a bounded and measurable function on the space of tessellations of @xmath4 . applying now the classical dynkin formula ( see for example lemma 19.21 in @xcite ) with @xmath96 gives the martingale property of the random process ( [ mart1 ] ) . ( it should be mentioned that so defined @xmath95 need not to be bounded in general . however , this technical difficulty can be overcome with a suitable localization argument . ) + for the second statement we consider the time - augmented random process @xmath97 , @xmath23 , which has generator @xmath98 given by @xmath99(y)+\left[{\partial\over\partial t}g(y,\cdot)\right](t)\ ] ] for appropriate functions @xmath100 . now , using dynkin s formula now for the product process @xmath97 and with @xmath101 , @xmath102 , gives ( again after a suitable localization argument ) the martingale property of the random process @xmath103 for another functional @xmath88 we can apply the latter property for @xmath104 and @xmath105 and subtract the two results , which shows the martingale property of the random process ( [ mart2]).@xmath106      in this subsection we discuss certain basic first - order consequences of proposition [ martprop ] to be of use for our future reference . we will write @xmath107 for the intrinsic volume of order @xmath29 with @xmath108 ( for the standard definition of these functionals we refer to @xcite and the references cited therein ) . further , write @xmath109 note now that whenever a new facet @xmath32 splits a cell @xmath42 into @xmath110 and @xmath111 of @xmath7 giving rise to a new tessellation @xmath112 , we have @xmath113 since @xmath107 has the valuation property . consequently , constructing the tessellation @xmath102 by successive cell splits we easily get @xmath114 it is our aim to relate @xmath115 and @xmath116 . for this purpose we use a special case of _ crofton s formula _ from classical integral geometry , which reads @xmath117}v_j(k\cap h)\lambda(dh)=\gamma_{j+1}v_{j+1}(k),\label{eqcrofton}\ ] ] where @xmath118 is a convex body and where the constant @xmath57 is given by ( [ gammaj ] ) , see ( * ? ? ? applying this formula to @xmath115 yields @xmath119}\sum_{f\in{\operatorname{cells}}(y\cap h)}v_j(f)\lambda(dh)=\int_{[w]}\sum_{c\in{\operatorname{cells}}(y)}v_j(c\cap h)\lambda(dh)\\   & = & \sum_{c\in{\operatorname{cells}}(y)}\gamma_{j+1}v_{j+1}(c)=\gamma_{j+1 } f_{j+1}(y).\label{avjtechnicaltool}\end{aligned}\ ] ] to streamline the notation , we shall write @xmath120 and @xmath121 below . note that , upon combining ( [ fjsj ] ) and ( [ avjtechnicaltool ] ) , it follows from ( [ mart1 ] ) with @xmath122 that @xmath123 are all @xmath91-martingales for @xmath124 in particular , bearing in mind that @xmath125 we conclude that @xmath126 is a @xmath91-martingale , in particular @xmath127 , where we have used @xmath128 . the latter equation is extended by    [ propofirstorder ] for @xmath129 we have @xmath130 with the convention that @xmath131 and @xmath132 .    [ [ proof ] ] proof : + + + + + +    taking expectations in ( [ martexpl1ord ] ) we get @xmath133 continuing recursively by applying ( [ fjsj ] ) and ( [ avjtechnicaltool ] ) we end up with @xmath134 however , the definitions of @xmath39 and @xmath135 imply that @xmath136 for any @xmath137 . thus , @xmath138 and the result follows immediately.@xmath106 +   + it is interesting to compare the mean value formula from the last proposition with that from @xcite , section 3.2 . we denote by @xmath139 the density of the @xmath29-th intrinsic volume of the collection of @xmath26-dimensional maximal polytopes of @xmath49 , i.e. @xmath140 for arbitrary compact windows @xmath1 with @xmath141 , see ( * ? ? ? using now proposition [ propofirstorder ] we get @xmath142 because of the homogeneity of intrinsic volumes . on the other hand , we have shown in @xcite that @xmath143 where here and in the sequel we denote by @xmath144 , @xmath145 , the volume of the @xmath146-dimensional unit ball . indeed , these two values are identical , since @xmath147    [ rema1 ] our first - order formula from the last proposition contains the extra boundary correction term @xmath148 in contrast to our mean value formula ( [ xxxeqmean ] ) from @xcite . these additional terms come from the fact that we consider maximal polytopes possibly chopped off by the boundary of the domain @xmath4 rather than the points of an associated center function for full facets in @xmath149 thus , it may happen for instance that in two neighbouring regions one observes two facets which can coalesce when putting these regions together into one volume , whence the lower - order finite volume corrections arise . in this section we develop a full second - order theory for intrinsic volumes of stationary and isotropic random stit tessellations . this is first done exactly in subsection [ subsecexact ] for @xmath150 inside a bounded convex observation window @xmath1 . then , in subsection [ subsecasympt ] we derive the corresponding asymptotic expressions for @xmath151      to proceed with second - order calculations we shall use the notation already introduced in subsection [ subsec1stord ] above . observe first that in view of ( [ fjsj ] ) the relation ( [ avjtechnicaltool ] ) simplifies and we obtain @xmath152 recalling that @xmath120 and @xmath153 putting likewise @xmath154 and using ( [ mart2 ] ) in proposition [ martprop ] with @xmath155 and @xmath156 , we see that @xmath157 ds\ ] ] is a martingale with respect to the filtration @xmath91 induced by @xmath158 . substituting ( [ baraj ] ) , this implies that @xmath159 ds\ ] ] is an @xmath91-martingale as well . this is a crucial formula because , upon taking expectations , it allows us to express @xmath160 in terms of corresponding covariances with indices @xmath161 and @xmath162 in other words , we get @xmath163 ds.\ ] ] it is important to observe that this recursion terminates because @xmath164 , which allows us to provide an explicit expression for the covariances @xmath160 . to obtain the desired formula , denote @xmath165 for any @xmath166 \to { \bbb r}$ ] and @xmath167 for which the iterated integral exists . then we claim that for @xmath168 it holds @xmath169 @xmath170 with the convention that @xmath171 and @xmath172 whereas this can be readily verified by a straightforward induction in view of ( [ principalcov ] ) , there is a more natural way to see it . in fact , in the course of recursive applications of ( [ principalcov ] ) the covariance @xmath173 is represented in terms of @xmath174 and @xmath175 which can be interpreted in terms of lattice walks on pairs of indices from @xmath176 to @xmath177 where only steps @xmath178 and @xmath179 are allowed , each step @xmath178 involving multiplication by @xmath180 plus integral iteration and each step @xmath179 resulting in multiplication by @xmath181 plus integral iteration . there are @xmath182 such walks , whence the first line in ( [ auxeq ] ) follows . the second line of ( [ auxeq ] ) is due to additional @xmath183 terms in ( [ principalcov ] ) , which are born at all times @xmath184 of the lattice walk discussed above , with @xmath185 standing for the number of @xmath178 steps and @xmath186 for the number of @xmath179 steps . note that no additional term is born at @xmath187 though , as it corresponds to having reached the ( co)variance @xmath188 , which does not get expanded any further . in ( [ auxeq ] ) we have substituted @xmath189 and @xmath190 for notational convenience . + applying now ( [ principalcov ] ) for @xmath191 and using that @xmath164 we see that @xmath192 whence we finally get from ( [ auxeq ] ) the following exact variance expression :    [ thmconvexact ] the covariance between the intrinsic volume processes @xmath193 and @xmath194 for @xmath168 of a stationary and isotropic random stit tessellation @xmath7 is given by @xmath195 @xmath196 with the usual convention that @xmath171 and @xmath172    unfortunately , we are not able to make the covariance in the above theorem any more explicit . this is due to the presence of mixed moments @xmath197 whose evaluation is technically related to the problem of providing a general formula for joint moments of lower - order intrinsic volumes of @xmath26-dimensional poisson cells ( arising as hyper - planar sections of @xmath7 , which are in addition possibly chopped off by the boundary of the window @xmath4 ) which is not currently available in required generality up to the best of our knowledge . fortunately though , the offending mixed moments turn out to be of negligible order in large window size asymptotics and we are able to provide fully explicit asymptotic formulae in subsection [ subsecasympt ] below . + we would like to point out that in the special case @xmath8 our theorem [ thmconvexact ] reduces to theorem 1 and corollary 1 in @xcite where we have studied variances and central limit theory for maximal edge ( i - segment ) lengths and vertex counts and where we could give fully explicit exact variance formulae exploiting the particular features of the planar setting . in corollary 2 ibidem we also have provided asymptotic expressions for these variances for sequences of growing observation windows . for this reason , in our asymptotic considerations below we will restrict to the case @xmath9 where the general asymptotic expressions are essentially different from those arising in the exceptional planar case . this is due to the variance dichotomy established in @xcite and mentioned in section [ secintro ] above . let @xmath1 be a compact convex window and consider the sequence @xmath198 . we will write from now on @xmath199 instead of @xmath200 and likewise @xmath201 for @xmath202 in order to emphasize the dependence on @xmath12 . our main interest in this section is to derive from theorem [ thmconvexact ] an asymptotic expression for the covariances @xmath203 as @xmath204 . we start with the following    [ propaux1 ] for @xmath205 with @xmath0 fixed we have @xmath206    [ [ proof-1 ] ] proof : + + + + + +    observe first that , in view of ( [ itint ] ) , the integral @xmath207 only involves integration with respect to the variable @xmath208 and does not affect the order in @xmath12 . thus , it is sufficient to prove the proposition for @xmath209 and @xmath210 which we shall henceforth assume without loss of generality . we have @xmath211 where @xmath212 and @xmath213 to provide a bound for @xmath214 we need some additional notation . write @xmath215 for the common value of @xmath216 with @xmath217 ranging through hyperplanes in @xmath3 and where @xmath218 the typical cell of the sectional stit tessellation @xmath219 . using the scaling property ( [ stitscaling ] ) and the homogeneity of the intrinsic volumes we readily get @xmath220 indeed , @xmath221 to proceed , write , recalling ( [ adef ] ) , @xmath222}\sum_{f\in{\operatorname{cells}}(y(s , w_r)\cap h)}v_{d-1-m}(f)v_{d-1-n}(f)\lambda(dh)\\ \nonumber & \leq &   { \bbb e}\int_{[w_r]}\sum_{f\in{\operatorname{cells}}(y(s)\cap h),\ ; f \cap w_r \neq \emptyset }    v_{d-1-m}(f)v_{d-1-n}(f)\lambda(dh)\\ \nonumber & = & \int_{[w_r ] } \frac{\varsigma_{m , n;s}}{\varsigma_{0,d-1;s } }    { \bbb e}{\operatorname{vol}}_{d-1}({\operatorname{typicalcell}}(y(s ) \cap h ) \oplus ( -(w_r \cap h ) ) ) \lambda(dh),\end{aligned}\ ] ] where @xmath223 stands for the usual minkowski addition . recalling that @xmath224 is poisson , see section [ subsecmnwc ] , we readily get for @xmath225 , @xmath226 and thus we conclude in view of ( [ bound1 ] ) that , for @xmath227 @xmath228 consequently , @xmath229 next , we find a bound for @xmath230 to this end , we note that during the time interval @xmath231 $ ] there are @xmath232 cell splits within @xmath233 in the course of the mnw - construction and hence the expectations of sums @xmath234 are of order @xmath235 uniformly in hyperplanes @xmath236 $ ] and in @xmath237.$ ] thus , using that @xmath238 ) = o(r)$ ] and recalling the definition ( [ adef ] ) we obtain @xmath239 putting ( [ bounds1 ] ) and ( [ bounds2 ] ) together completes the proof of proposition . @xmath106 +   + this implies that asymptotically , as @xmath204 , the terms appearing in theorem [ thmconvexact ] with @xmath240 are of order at most @xmath241 and thus negligible compared with the leading @xmath242-term . for this reason , the asymptotic behaviour of @xmath243 is dominated by the term with @xmath244 , which is of order @xmath245 where by @xmath246 we mean quantities bounded both from below and above by multiplicities of the argument in the usual landau notation . for this dominating case proposition [ propaux1 ] is refined by    [ propasympt ] we have for @xmath168 @xmath247    [ [ proof-2 ] ] proof : + + + + + +    first , recall from thm . 4 in @xcite that the variance of the total surface area of @xmath7 equals @xmath248 the main argument for ( [ totalvar ] ) reads as follows : at first , we use the fact that stit tessellations have poisson typical cells and find @xmath249}\sum_{f\in{\operatorname{cells}}(y(s , w)\cap h)}v_{d-1}^2(f)\lambda(dh)\\ \nonumber & = & \int_{[w]}\int_{w\cap h}\int_{w\cap h}{\bf 1}[x , y\ \text{are in the same cell of}\ y(s , w)\cap h]dxdy\lambda(dh)\\ \nonumber & = & \int_{[w]}\int_{w\cap h}\int_{w\cap h}e^{-s\lambda([\overline{xy}])}dxdy\lambda(dh),\end{aligned}\ ] ] where @xmath250 $ ] stands for the usual indicator function and @xmath251 for the line segment connecting @xmath252 and @xmath253 . using now ( [ mart2 ] ) with @xmath254 , taking expectations and noting that the mixed term vanishes leads to @xmath255}\int_{w\cap h}\int_{w\cap h}e^{-s\lambda([\overline{xy}])}dxdy\lambda(dh)\\ \nonumber & = & \int_{[w]}\int_{w\cap h}\int_{w\cap h}{1-e^{-t\lambda([\overline{xy}])}\over\lambda([\overline{xy}])}dxdy\lambda(dh).\end{aligned}\ ] ] the latter integral can be transformed into ( [ totalvar ] ) by using an integral - geometric formula of blaschke - petkantschin - type and the crofton formula ( [ eqcrofton ] )  see @xcite for details . + recall now ( [ vard1expr ] ) and write @xmath256 where the last equality follows by the scaling property ( [ stitscaling ] ) of stit tessellations . thus , @xmath257 now , the binomial theorem implies that there exists some constant @xmath258 such that @xmath259 with @xmath260 . in view of ( [ eeee1 ] ) this proves the desired result.@xmath106 +   + consequently , by combining propositions [ propaux1 ] and [ propasympt ] with theorem [ thmconvexact ] we get    [ corasycov ] asymptotically , as @xmath204 we have for @xmath168 @xmath261 @xmath262 where @xmath263 denotes the @xmath264-energy of @xmath4 , i.e. @xmath265    the affine blaschke - petkantschin formula ( * ? ? ? 7.2.7 ) can be used to provide an integral - geometric expression for the measure - geometric energy functional @xmath266 . in fact , we have @xmath267 where @xmath268 denotes the @xmath26-st chord power integral of the convex body @xmath4 in the sense of @xcite , page 363 , i.e. @xmath269 where @xmath270 denotes the space of lines in @xmath3 with invariant measure @xmath271 and @xmath272 is the volume of the @xmath273-dimensional unit ball . this means that in the asymptotic covariance formula from corollary [ corasycov ] the dependence on the geometry of @xmath4 is encoded by the non - additive @xmath263 or equivalently by @xmath268 . we will from now on use the representation in terms of chord power integrals , as it allows an easier comparison with other tessellation models . summarizing , this yields the following    asymptotically , as @xmath204 , we have for @xmath168 @xmath274 and for @xmath275 @xmath276    especially for the practically relevant case @xmath277 , we have for @xmath278 and @xmath279 @xmath280 in general , @xmath268 or equivalently @xmath263 is rather difficult to evaluate explicitly . however , for the unit ball @xmath281 in @xmath3 we have by applying @xcite , theorem 8.6.6 ( with a corrected constant ) @xmath282 for example for @xmath277 this gives us the value @xmath283 . for the interpretation of computer simulations it is of particular interest to evaluate the chord power integral @xmath284 for a @xmath285-dimensional cube @xmath286 with side length @xmath287 . it can be shown that the numerical value of @xmath284 is given by @xmath288 . in this section we focus our attention on the spatial pair - correlation structure for the processes of maximal polytopes of arbitrary dimensionalities induced by @xmath289 arising for @xmath108 as random @xmath29-volume measures concentrated on the union of all @xmath29-faces of @xmath26-dimensional maximal polytopes of @xmath49 . the nature of the so defined face measures is somewhat different than that of cumulative intrinsic volume processes considered above , in particular the face measures keep track not just of the cumulative numeric characteristics of stit tessellations but also of their spatial profile , moreover it should be emphasized that in general the total mass of a face measure may be quite unrelated to the cumulative intrinsic volume of the corresponding order , as for example @xmath290 , the number of @xmath26-dimensional maximal polytopes of @xmath49 in @xmath1 , is not deterministically related to the number of vertices of @xmath7 as soon as @xmath15 etc . however , we decided to consider the face measures in this paper as they are of interest in their own right and supplement our results from the other sections , showing the power and versatility of the methods developed in this paper . + for @xmath291 we consider the ( random ) @xmath29-th order face measure @xmath292 where @xmath293 is the @xmath294-truncated @xmath29-volume measure @xmath295(a ) = { \operatorname{vol}}_j(a \cap e),\ ; a \subseteq { \bbb r}^d , $ ] whereas @xmath296 where by @xmath31 we mean the collection of all @xmath29-faces of the @xmath26-dimensional polytope @xmath32 . we also abuse the notation by putting @xmath297 we shall be interested in the covariance measures @xmath298 , on @xmath299 given by @xmath300 for all bounded measurable @xmath301 with bounded support , where the standard duality notation @xmath302 is used for integration @xmath303 denote by @xmath304 the law of the typical @xmath26-dimensional maximal polytope , of the stit tessellation @xmath49 and let @xmath305 be the corresponding facet density . the following proposition is crucial for this section :    [ covrecursionlemma ] for @xmath306 we have @xmath307 @xmath308 ds\ ] ] with @xmath309 and @xmath310 given as in ( [ gammaj ] ) . [ [ proof-3 ] ] proof : + + + + + +    for bounded measurable and boundedly supported @xmath311 and @xmath312 and for @xmath313 consider the facet functionals @xmath314 and @xmath315 defined analogously . then , choosing some compact convex @xmath4 containing the supports of @xmath311 and @xmath316 in its interior we see that , recalling ( [ sigmadef ] ) and ( [ adef ] ) , @xmath317 } \sum_{f \in { \operatorname{cells}}(y \cap h ) } j^g_i(f ) \lambda(dh ) =        \sum_{e \in { \operatorname{maxpolytopes}}_{i+1}(y ) } \left\langle g ,   \int_{[w ] } { \operatorname{vol}}_i(\cdot \cap ( e \cap h ) ) \lambda(dh )        \right\rangle\ ] ] @xmath318 where we have used the crofton formula ( [ eqcrofton ] ) . applying now ( [ mart2 ] ) with @xmath319 and @xmath320 we get upon taking expectations @xmath321 @xmath322 ds\ ] ] and thus , in view of ( [ aformeas ] ) , @xmath323 @xmath324 ds.\ ] ] using that @xmath325 and similar relationships , we end up with @xmath326 @xmath327   ds \right\rangle.\ ] ] thus , to establish ( [ covmeasrecursion ] ) it is enough to show that @xmath328 to prove ( [ covaux2 ] ) use ( [ stitsection ] ) and write , recalling ( [ adef ] ) , @xmath329 } { \bbb e } \sum_{f \in { \operatorname{cells}}(y(s , w ) \cap h ) } j^g_i j^h_j(f)\lambda(dh)\\   & = & \int_{[w ] } \int j^g_i j^h_j(f ) { \bbb m}^{y(s , w ) \cap h}(df ) \lambda(dh),\label{eqxxx }   \end{aligned}\ ] ] where @xmath330 is the cell intensity measure for the sectional stit tessellation @xmath331 , compare with ( [ cellintensitymeasure ] ) . using proposition [ propxx ] ( a ) we are led to @xmath332 } \int j^g_i j^h_j(f ) { \bbb m}^{y(s , w ) \cap h}(df ) \lambda(dh)=\int_{[w ] } \int j^g_i j^h_j(f ) { \bbb m}^{{\operatorname{pht}}(s , w ) \cap h}(df ) \lambda(dh).\ ] ] now , using slivnyak s theorem ( cf . 3.2.5 ) ) we obtain for any bounded measurable function @xmath77 on the space of @xmath26-dimensional polytopes @xmath333}\int\phi(f){\bbb m}^{{\operatorname{pht}}(1,w)\cap h}(df)\lambda(dh ) & = & \int\int_{[c]}\phi(c\cap h)\lambda(dh){\bbb m}^{{\operatorname{pht}}(1,w)}(dc)\\ \nonumber & = & \int\phi(f){\bbb f}_{d-1}^{{\operatorname{pht}}(1,w)},\end{aligned}\ ] ] where @xmath334 is the cell intensity measure of @xmath335 , whereas @xmath336 is the cell intensity measure and @xmath337 is the @xmath26-face intensity measure of a poisson hyperplane tessellation @xmath338 within @xmath4 having intensity measure @xmath17 . thus , replacing @xmath17 by @xmath339 , @xmath340 , we arrive at @xmath332}\int\phi(f){\bbb m}^{{\operatorname{pht}}(s , w)\cap h}(df)\lambda(dh)={1\over s}\int\phi(f){\bbb f}_{d-1}^{{\operatorname{pht}}(s , w)}(df),\ ] ] whence continuing ( [ eqxxx ] ) we find @xmath341 hence , by applying proposition [ propxx ] ( b ) we get from ( [ aaux1 ] ) , @xmath342 where , recall , @xmath343 is the @xmath26-dimensional maximal polytope intensity measure for @xmath7 , see ( [ faceintensitymeasure ] ) . in view of ( [ aaux2 ] ) the required relation ( [ covaux2 ] ) follows now directly by the definition of typical @xmath26-dimensional maximal polytope ( * ? ? ? * chap . 4.1 , chap . 10 ) and the definition of @xmath344 upon taking into account that the supports of both @xmath311 and @xmath316 are contained in the interior of @xmath345 putting ( [ covaux1 ] ) and ( [ covaux2 ] ) together yields ( [ covmeasrecursion ] ) and thus completes the proof . @xmath106 +   + to proceed note that @xmath346 is a constant measure and therefore any covariances involving it vanish . thus , using ( [ covmeasrecursion ] ) in proposition [ covrecursionlemma ] and arguing as in the derivation of theorem [ thmconvexact ] from the crucial relation ( [ principalcov ] ) we obtain    [ thmcovmeas ] for @xmath347 we have @xmath348 @xmath349          { \cal i}^{k+l - m - n}\left(\lambda^{(d-1)}_{y(\cdot ) } { \bbb q}^{y(\cdot)}_{d-1};t \right)(df ) dx\ ] ] with the usual convention that @xmath171 and @xmath172    it is interesting to note that for the particular case @xmath350 we get    [ covcor1 ] we have @xmath351{\bbb q}^{y(t)}_{d-1}(df ) dx.\ ] ] this means that for @xmath350 the covariance measure of the surface area process @xmath352 coincides with that of the surface area process induced by the homogeneous and isotropic boolean model with grain distribution @xmath304 and with grain density @xmath353 a result first established by weiss , ohser and nagel @xcite in the special planar case by completely different methods . recall further from theorem 3 in @xcite that @xmath354 where @xmath355 stands for the distribution of the typical facet of the poisson hyperplane tessellation @xmath356 with surface intensity @xmath208 ( this is to say , the hyperplane process has intensity measure @xmath339 ) . but this random polytope has the same distribution @xmath357 as the typical cell of a poisson hyperplane tessellation in @xmath358 with surface intensity @xmath359 with @xmath360 given by ( [ gammaj ] ) . thus , @xmath361{\bbb q}_{{\operatorname{typicalcell}}}^{{\operatorname{pht}}(\gamma_{d-1}s,{\bbb r}^{d-1})}(dc)dxds\ ] ] and the covariance measure is reduced to known quantities . + for @xmath362 the situation becomes more complicated and mixtures of typical cell distributions corresponding to different time moments arise in the right hand side of ( [ covmeasfullexpr ] ) . as in section [ sec2ndord ] above , also here explicit calculations are precluded for @xmath362 due to the lack of known formulae for mixed moments of general order intrinsic volumes of poisson cells . in this section we present a central limit theory for the suitably rescaled intrinsic volume processes @xmath363 in strong contrast to the two - dimensional case considered in @xcite , where a classical gaussian limit behaviour is observed , for the case @xmath364 in focus of this paper the situation is very different and non - gaussian limits arise . in this context , to proceed with a full discussion below , we recapitulate first some facts already known from @xcite in a way specialized for our present purposes . define the rescaled intrinsic volume processes @xmath365}$ ] for @xmath366 by @xmath367.\ ] ] note that the shift by @xmath368 in time argument here is of technical importance as placing the time origin just after the very initial _ big bang _ phase of the mnw - construction , where the dominating fluctuations of @xmath369 arise , but where nothing of asymptotic significance happens for @xmath370 since the big bang phase evolution for @xmath369 has been considered in full detail in ( * ? ? ? 5 ) , in this paper we only focus on the later phase @xmath371.$ ] putting together the present section with ( * ? ? ? 5 ) reveals remarkable richness of the complete asymptotic picture for stit tessellations in large windows . + we begin with the observation that follows by the theory developed in subsection 5.3 of @xcite :    [ vd1 ] the process @xmath372}$ ] converges in law , as @xmath373 in the space @xmath374 $ ] of right - continuous functions with left - hand limits ( cdlg ) on @xmath375 $ ] endowed with the usual skorokhod @xmath376-topology ( * ? ? ? vi.1 ) to the constant process @xmath377 where @xmath378 is a certain non - gaussian square - integrable random variable with variance @xmath379      the relation ( 85 ) in @xcite implies that @xmath380 converges in law to @xmath381 as @xmath382 on the other hand , recall that @xmath383 is a martingale in view of ( [ martexpl1ord ] ) . moreover , by ( [ totalvar ] ) and the definition of @xmath383 combined with the scaling property of stit tessellations we have @xmath384 ||x - y||}}{||x - y||^2 }         dx dy    \end{aligned}\ ] ] and hence @xmath385 tends to @xmath386 as @xmath382 consequently , the asymptotic constancy of the limit process follows now by doob s @xmath387-maximal inequality ( * ? ? ? * thm 3.8(iv ) ) which completes our argument . @xmath106 +   +    [ nongaussexpl ] the convergence in law and non - gaussianity results referred to above come from the paper @xcite . in order to keep the present work formally self - contained we give here a very brief sketch of the main arguments . + the convergence is readily guaranteed by the martingale convergence theory as ibidem . the crucial point is the non - gaussianity . here , the idea , as described in detail in section 5.3 of @xcite , relies on showing that @xmath381 has its tails much heavier than normal random variables . we proceed , roughly speaking , by constructing a suitable class of _ initial frame tessellations _ of the body @xmath4 with the properties that    * the number of frame hyperplanes hitting @xmath4 equals @xmath388 , * the frame hyperplanes intersect within @xmath4 only very seldom , which implies that @xmath4 gets subdivided into @xmath389 cells ( recall the landau notation ) . this can be achieved by choosing @xmath273 principal directions and keeping all hyperplanes in the collection under construction approximately aligned to these directions . importantly , it turns out that this way we can get the lower bound @xmath390 for the probability that the real initial frame arising in the very initial phase for @xmath391  usually called _ big bang _ , see below  upon suitable spatial re - scaling does fall into the desired class . now , denoting the last event by @xmath392 , we see that on @xmath392 the total deviation @xmath393 gets decomposed into    * a sum of independent contributions @xmath394 coming from all respective cells @xmath395 of the frame tessellation , * a correction term of order @xmath389 due to the centred joint contribution of all hyperplanes of the frame tessellation ( observe that each such hyperplane necessarily contributes @xmath396 , whereas @xmath397 ) .    finally we can use the knowledge of @xmath398 and @xmath399 and properly exploit the independence of the random variables @xmath394 , to conclude that on @xmath392 we can find a further sub - event @xmath400 with @xmath401 such that @xmath402 and hence @xmath403 , which could not hold if @xmath381 were gaussian . we refer the reader to @xcite for further details . [ ivlimtheo ] the vector @xmath404}$ ] converges in law , as @xmath373 in the space @xmath405;{\bbb r}^d)$ ] of @xmath3-valued cdlg functions on @xmath375 $ ] endowed with the usual skorokhod @xmath376-topology ( * ? ? ? * chap.vi.1 ) , to the stochastic process @xmath406.\ ] ]    before proceeding with the proof we discuss certain striking features of the phenomenon described in theorem [ ivlimtheo ] . namely , although all intrinsic volume processes @xmath407 exhibit fluctuations of order @xmath408 the nature of these fluctuations differs very much between @xmath409 and @xmath410    * as shown in subsection 5.3 in @xcite and in proposition [ vd1 ] above , the leading - order deviations of @xmath411 in large @xmath12 asymptotics arise very early in the course of the mnw - construction , in its initial stages usually referred to as the _ big bang phase_. here , this is the time period @xmath412.$ ] during the later stages of the construction , i.e. say the time interval @xmath413,$ ] the variance increase is of lower order and any newly arising fluctuations are negligible compared to those originating from the _ big bang_. in the asymptotic picture this means that the initial fluctuation remains _ frozen _ throughout the rest of the dynamics , whence the constant limit for @xmath372}$ ] as @xmath382 * in contrast , the leading - order deviations of @xmath414 for @xmath415 arise and cumulate constantly in @xmath0 with deterministic polynomial rates depending on @xmath416 times the initial _ big bang _ fluctuation of the process @xmath417 which , in this sense , _ stores the entire randomness _ of the intrinsic volume vector . the mechanism determining the dependence of fluctuations of intrinsic volume processes of orders @xmath415 given those for @xmath418 and the resulting form of the limit process ( [ limproc ] ) will be discussed in the sequel . the crucial step of the proof relies on considering for each @xmath108 the auxiliary process @xmath419 which is a centred @xmath91-martingale by ( [ mart1 ] ) and which is the same as @xmath420 by ( [ baraj ] ) . the idea is to show that for @xmath421 the processes @xmath422 are of negligible order in large @xmath12 asymptotics . indeed , upon squaring and taking expectations we get @xmath423 using that @xmath424 as follows by the martingale property of @xmath425 and inserting this to ( [ varsigmaaux ] ) we are led to @xmath426 @xmath427 with the last two terms cancelling . thus , @xmath428 whence , by ( [ principalcov ] ) , @xmath429 in view of proposition [ propaux1 ] for @xmath421 we have @xmath430\ ] ] and thus the relation ( [ varsigmaaux2 ] ) implies @xmath431\ ] ] for @xmath432 hence finally @xmath433 }   ( \hat\sigma_{v_j;t}^{w_r})^2 = o(r^{2(d-1)})\ ] ] for @xmath421 by doob s @xmath387-maximal inequality ( * ? ? ? * thm 3.8(iv ) ) applied to the martingale @xmath434 defined as in ( [ sigmaaux ] , [ sigmaaux2 ] ) with @xmath4 replaced by @xmath233 there according to our usual notational convention . + with ( [ finalnegl ] ) it is now easy to complete the proof . indeed , since the normalization in the definition of intrinsic volume processes @xmath435 involves a prefactor @xmath436 the relation ( [ finalnegl ] ) allows us to recursively substitute @xmath437 for @xmath438 , as soon as @xmath421 , without affecting the large @xmath439 asymptotics in law ( note that the technically motivated shift by @xmath440 in time argument of the rescaled process @xmath435 is asymptotically negligible as inducing only a negligible @xmath387-difference precisely calculated in theorem [ thmconvexact ] ) . the application of such recursive substitutions combined with proposition [ vd1 ] completes the proof of theorem [ ivlimtheo ] . @xmath106 +   + it is crucial to emphasize at this point that , in intuitive terms , under the normalization of @xmath441 @xmath442 the mechanism governing the rise of fluctuations of intrinsic volume processes of orders @xmath443 given those for @xmath135 reduces effectively to the simple approximation @xmath444 and its recursive application . of course this _ simple approximation _ follows itself by rather non - trivial arguments above . clearly , this description only characterizes the leading order fluctuations as considered in theorem [ ivlimtheo ] although our tools should allow a more delicate characterization of non - leading lower order fluctuations as well , see e.g. subsection 5.1 in @xcite for the particular case of @xmath445      we are grateful to claudia redenbach and joachim ohser who provided the two pictures of stit tessellations . we would also like to thank the anonymous referees for their comments and suggestions as well as the editor - in - chief for his help . + the first author was supported by the polish minister of science and higher education grant n n201 385234 ( 2008 - 2010 ) . the second author was supported by the swiss national science foundation , grant pp002 - 114715/1 . t. schreiber , c. thle , second - order properties and central limit theory for the vertex process of iteration infinitely divisible and iteration stable random tessellations in the plane , adv * 42 * , 913935 ( 2010 ) .
Abstract:

Answer: stationary and isotropic iteration stable random tessellations are considered , which can be constructed by a random process of cell division . the collection of maximal polytopes at a fixed time @xmath0 within a convex window @xmath1 is regarded and formulas for mean values , variances , as well as a characterization of certain covariance measures are proved . the focus is on the case @xmath2 , which is different from the planar one , treated separately in @xcite . moreover , a multivariate limit theorem for the vector of suitably rescaled intrinsic volumes is established , leading in each component  in sharp contrast to the situation in the plane  to a non - gaussian limit . # 1    * key words : * central limit theory ; integral geometry ; intrinsic volumes ; iteration / nesting ; markov process ; martingale ; random tessellation ; stochastic stability ; stochastic geometry + * msc ( 2000 ) : * primary : 60d05 ; secondary : 52a22 ; 60f05


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: dust emission gives us crucial information to understand both the history of cosmic star formation in galaxies and that of supermassive black hole ( smbh ) growth in galactic centers . intense star formation produces a great amount of dust , which makes its activity invisible in the ultraviolet band but visible in the infrared band . similarly , the central engines of active galactic nuclei ( agns ) are surrounded by dusty `` tori '' @xcite . since optical and ultraviolet lights are very easily obscured by the torus , the most complete search for agns can be made by detecting hard x - rays from the central engine or infrared emission from the heated dust . utilizing x - ray data of the chandra deep field south , @xcite reported that the fraction of `` compton - thick '' agns ( those with absorption column densities of @xmath6@xmath7 ) increases with redshift from @xmath8 to @xmath9 . this fact implies that the level of dust - obscuration of agn may be correlated to the star formation activity , which has a peak at @xmath10 @xcite . indeed , @xcite gathered all available _ spitzer _ spectra of local compton - thick agns and found strong starburst ( sb ) features in their average infrared spectrum , suggesting that sb activity and highly obscured tori are somehow coupled each other . @xcite also reported the same trend for buried agns with small opening angles of tori (  new type  agns ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) detected in the _ swift_/bat survey . although hard x - rays are useful to search for obscured agns , extremely heavily buried agns with @xmath11@xmath7 are difficult to be detected due to flux attenuation by repeated comptonization even at energies above 10 kev ( e.g. , * ? ? ? * ; * ? ? ? candidates of galaxies that host such buried agns are infrared galaxies , defined as those having large infrared luminosities ( ) with @xmath12 . among infrared galaxies , those in the range @xmath13 are called luminous infrared galaxies ( lirgs ; * ? ? ? * ) and more luminous ones in the range @xmath14 are called ultraluminous infrared galaxies ( ulirgs ; * ? ? ? their bolometric luminosities are dominated by the infrared emission , suggesting that very luminous heating sources are surrounded by dust and then the heated dust re - emits in the infrared band . the hidden energy sources are believed to be either sb , agn , or both . disentangling the energy sources of infrared galaxies is crucial to unveil the history of dust obscured star formation and smbh growth in the universe @xcite . @xcite studied agn activity in u / lirgs on the basis of optical observations . they showed that the fraction of agns increases with infrared luminosities in the range of @xmath15 . one caveat of their optical diagnostics is that one can not identify buried agns embedded in dust - tori with very small opening angles where no significant narrow line regions can be formed . unlike the agns which are optically - classified type-2 seyfert , such buried agns are elusive to conventional optical spectroscopy @xcite . after the launch of _ spitzer _ , the agn activity in local u / lirgs has been also addressed in the infrared band . @xcite observed ulirgs with _ spitzer_/irs 58  @xmath0 m spectroscopy , and investigated both the sb and agn contributions by fitting sb / agn templates to the spectra . they found the increasing agn significance with infrared luminosity , confirming the same trend reported from the optical studies . @xcite observed the 538  @xmath0 m spectra of lirgs with _ spitzer_. they decomposed the spectra into sb and agn components using a sb galaxy template and clumpy torus models @xcite , which are supported by observations ( e.g. , * ? ? ? * ; * ? ? ? they confirmed that the agn energetic contribution to the total power increases with infrared luminosity from @xmath16 . infrared 2.55 @xmath0 m spectroscopy is another powerful tool to study optically - elusive buried agns in u / lirgs . one advantage of this band is that dust extinction is much lower than in the optical band @xcite and is similar to that in the 513 @xmath0 m band @xcite . another great advantage is that sb and agn activity can be distinguished based on the spectral features . first , strong polycyclic aromatic hydrocarbon ( pah ) emission , a pure star formation tracer , is located at 3.3  @xmath0 m in sb galaxies , while pure agns exhibit pah - feature free spectra due to the effects that 1 ) the x - ray emission destroying the pahs , and 2 ) the strong continuum emission originating from agn hot - dust diminishes this feature @xcite . sb galaxies generally show large 3.3  @xmath0 m pah equivalent widths of @xmath17  nm , which never go down below 40  nm . thus , objects with @xmath18  nm may be classified as galaxies that harbor buried agns . in sb / agn composite galaxies , strong pah emission can be observed because the agn can not destroy the pah molecules located in the outer ( @xmath19  pc ) region from the central engine due to shielding by dust and gas . however , strong agn continuum makes ew@xmath20 smaller , well below @xmath21  nm . using this method , many authors reported the strong sign of buried agns in ulirgs @xcite and lirgs @xcite . another key feature in the 2.55  @xmath0 m band is the continuum slope @xmath22 ( @xmath23 ) . sb galaxies have blue continuum slopes ( @xmath24 ) in this band due to the contribution of the stellar photospheric continuum . conversely , galaxies with agns have hot - dust emission heated by agns , which produces a red continuum ( @xmath25 ) . this method was first reported by @xcite and successfully applied to ulirgs @xcite and lirgs @xcite . @xcite combined the two methods ( @xmath26  nm and/or @xmath27 ) to find buried agns from a sample of u / lirgs showing no agn signatures in the optical spectra . they found that the fraction of buried agns increases with infrared luminosity in the range of @xmath28 . in this paper , we expand the studies of u / lirgs by @xcite into a lower infrared luminosity range , by including normal infrared galaxies ( irgs ) with @xmath29 , where comprehensive studies are still missing due to the limited sample . this is because l / irgs have spatially more extended infrared emission ( @xmath30  kpc , corresponding to @xmath31several arcseconds at @xmath32 ) compared with ulirgs @xcite , and therefore slit spectroscopy with ground based telescopes could miss their extended emission . slit - less spectroscopy of infrared camera ( irc ; * ? ? ? * ) on board _ akari _ @xcite with 1@xmath33 aperture opened a new window to probe such extended emission of l / irgs . throughout this paper , we adopt @xmath34  km s@xmath35mpc@xmath35 , @xmath36 , and @xmath37 for consistency with the previous publications @xcite . our motivation is to search for buried agns from optically non - seyfert infrared galaxies and investigate their properties over a wider range of infrared luminosity than in previous studies . an ideal sample can be selected from unbiased catalogs of infrared galaxies with various infrared luminosities at redshifts of @xmath38 , where 3.3  @xmath0 m pah is detectable in the _ akari_/irc 2.55.0  @xmath0 m spectra . to this end , we first gather u / lirgs with high far - infrared fluxes from _ iras _ catalogs . our ulirgs are mainly selected from the bright ulirg catalog by @xcite , consisting of 41  sources 25 out of which are optically non - seyfert galaxies . these ulirgs have _ iras _ 60  @xmath0 m fluxes @xmath39  jy and far - infrared ( 40120  @xmath0 m ) luminosities @xmath40  km s@xmath35mpc@xmath35 . therefore , in our study with @xmath41  km s@xmath35mpc@xmath35 , some galaxies reach @xmath42 . ] . in addition to these sources , we also gather 54 fainter , local ( @xmath38 ) optically non - seyfert ulirgs from the literature to increase the sample size . the lirgs are mainly taken from the lirg catalog of @xcite , which is based on the _ iras _ bright galaxy catalog . the catalog consists of 61 sources with the criteria of _ iras _ 60  @xmath0 m fluxes @xmath43  jy and far - infrared ( 40400  @xmath0 m ) luminosities @xmath44 . we select 41 optically non - seyfert galaxies out of this sample . to further increase the number of l / irg targets , we also utilize the catalog by @xcite . they compiled 76 _ iso _ detected sources in the _ iras _ 12  @xmath0 m galaxy catalog , which gathered all galaxies at galactic latitudes of @xmath45 down to the _ iras _ 12  @xmath0 m flux limit of @xmath46  jy @xcite . since agn tori have peak emission at the mir band , 12  @xmath0 m galaxy catalog could tend to include more agns than in far - infrared selected samples . indeed , a majority of the @xcite catalog ( 64 out of 76 sources ) are optically seyfert galaxies . therefore , we gather the 12 optically non - seyfert starburst galaxies from it . in summary , we select a sufficiently large number of non - seyfert infrared galaxies ( 132 sources ) as our `` parent '' sample to be followed - up with _ akari_. the redshifts of all the sources in the above catalogs are distributed within @xmath38 .    , title="fig : " ] + observation was conducted with the irc infrared spectrograph @xcite on board _ akari _ @xcite . the spectra in the 2.55  @xmath0 m band were taken with the ng grism mode . this mode achieves an effective spectral resolution of @xmath47 at 3.6  @xmath0 m , which is sufficient for detecting and tracing the profile of @xmath48  @xmath0 m pah emission . the irc has a @xmath49 arcmin@xmath50 window with a pixel scale of @xmath51  arcsec@xmath50 . all the data were taken as a part of the _ akari _ mission program called `` agnul '' ( pi : t.  nakagawa ) . the observation settings were same as those described in chapter  3 of @xcite . one to seven pointings were assigned for each source according to the brightness . the total on - source exposure time was @xmath52  minutes per each pointing . we used ircz4 ( phase 3 ; post liquid - he mission ) observing mode , where one pointing is composed of eight or nine independent frames @xcite . this mode successfully removes cosmic - ray contaminations even for sources observed by a single pointing . for sources observed with multiple pointings , we combined all the data to achieve the best signal - to - noise ( s / n ) ratio after excluding the data sets with particularly bad quality . we often had to discard data observed in later epochs of _ akari _ s phase 3 operation due to the increasing background signal . spectral analysis was conducted in a standard manner by using the idl data reduction package for _ akari _ irc spectra . the dark - subtraction , linearity correction , and flat - field correction were performed using this idl package . some infrared galaxies with low luminosities ( @xmath53 ) have spatially extended emission . hence , we varied the aperture size for spectral extraction according to the actual signal profile of each source . the background was estimated from both sides of the spectral dispersion direction of the target and was subtracted from the source . the idl package for _ akari _ also performed wavelength and absolute flux calibration . the accuracy of the wavelength calibration is @xmath54  pixel ( @xmath55  @xmath0 m ) , and that of the flux calibration is @xmath56% around the central wavelength of the spectra and up to @xmath57% at the edges close to 2.5  @xmath0 m and 5.0  @xmath0 m .    out of the 132 sources in our parent sample , we observed 37 objects in total from 2008 june until 2010 february with _ akari_. after the epoch , _ akari _ fell into the stage where the data highly suffered from the significant background emission . as the result , we could not observe the remaining 95 targets . also , due to _ akari _ s sun - synchronous polar orbit , higher visibilities were achieved for objects at higher ecliptic latitudes . these _ akari _ s observational constraints make it difficult for us to construct a uniformly flux - limited sample over the entire sky with _ akari_. however , the observable targets are essentially randomly selected from the parent sample with no selection biases regarding the physical nature of the galaxies . hence , we regard that our sample has no obvious bias . next , we also set the criterion that the averaged s / n ratio for each spectral channel must be higher than 5 in order to apply our spectral deconvolution analysis . among the 37 targets observed with _ , five sources were discarded due to their low s / n ratio and the other 10 sources were too faint to be detected . 14 out of 15 non - detected sources are fainter ulirgs gathered from the literature . these sources have 60  @xmath0 m flux  @xmath58  jy , which is roughly 10 times fainter than that of the detected sources ( @xmath59  jy ) . supposing that their infrared spectral energy distribution does not change drastically from that of the detected sources , we can regard that the most of the non - detected sources are intrinsically too faint to be observed with _ figure  1 shows the redshift distribution of the finally selected 22 infrared galaxies . lr*7c eso  286-ig19 & 1122176 - 001 & 0.043 & 0.28 & 1.90 & 12.71 & 10.58 & 12.00 & 2comp + ic  5135 & 1120236 - 001 & 0.016 & 0.63 & 2.14 & 16.67 & 26.27 & 11.33 & 2comp + iras  03068 - 5346 & 1122105 - 001 & 0.074 & 0.25 & 0.20 & 3.43 & 3.92 & 11.98 & 2comp + iras 03538 - 6432 & 1122106 - 001 & 0.30 & 0.25 & 0.25 & 0.96 & 1.54 & 13.01 & 2comp + iras 10494 + 4424 & 1122096 - 001 & 0.092 & 0.25 & 0.25 & 3.36 & 5.60 & 12.22 & 2comp + iras  17028 + 5817 & 1122101 - 001 & 0.11 & 0.25 & 0.08 & 2.49 & 4.05 & 12.22 & 2comp + irasf  07353 + 2903 & 1120200 - 001 & 0.33 & 0.08 & 0.10 & 0.22 & 0.64 & 12.63 & 2comp + mcg  + 02 - 04 - 025 & 1122123 - 001 & 0.031 & 0.37 & 1.46 & 11.13 & 10.29 & 11.67 & 2comp + mcg  + 08 - 23 - 097 & 1122130 - 001 & 0.029 & 0.25 & 0.53 & 5.08 & 8.10 & 11.34 & 2comp + mcg  + 10 - 19 - 057 & 1122132 - 001 & 0.031 & 0.40 & 1.92 & 11.35 & 10.81 & 11.70 & 2comp + mrk  1490 & 1122134 - 001 & 0.026 & 0.25 & 0.86 & 6.21 & 8.38 & 11.30 & 3comp + mrk  551 & 1122171 - 001 & 0.050 & 0.26 & 0.80 & 4.63 & 6.13 & 11.80 & 2comp + mrk  848 & 1122136 - 001 & 0.040 & 0.32 & 1.52 & 9.38 & 10.26 & 11.86 & 2comp + ngc  2339 & 1120230 - 001 & 0.0074 & 0.53 & 2.11 & 18.96 & 32.24 & 10.69 & 3comp + ngc  2388 & 1120231 - 001 & 0.014 & 0.49 & 1.99 & 16.21 & 23.09 & 11.15 & 3comp + ngc  4102 & 1120232 - 001 & 0.0028 & 1.45 & 6.83 & 47.37 & 68.50 & 10.25 & 2comp + ngc  4194 & 1122091 - 001 & 0.0083 & 0.86 & 4.36 & 22.79 & 25.94 & 10.88 & 3comp + ngc  4818 & 1120234 - 001 & 0.0036 & 0.85 & 3.88 & 20.26 & 26.55 & 10.10 & 3comp + ngc  520 & 1120229 - 001 & 0.0076 & 0.78 & 2.83 & 31.52 & 48.40 & 10.91 & 3comp + ngc  6285 & 1122138 - 001 & 0.019 & 0.30 & 0.44 & 7.43 & 23.58 & 11.21 & 3comp + ngc  838 & 1122173 - 001 & 0.013 & 0.59 & 1.83 & 0.40 & 17.94 & 10.74 & 3comp + zw  453.062 & 1122140 - 001 & 0.025 & 0.29 & 0.57 & 7.59 & 11.04 & 11.34 & 2comp +   + for understanding the properties of buried agns in infrared galaxies , it is crucial to analyze the infrared spectra by decomposing sb and agn components . the 2.55.0  @xmath0 m band covered by _ akari _ irc has unique advantages for detecting agn signatures . as introduced in section  1 , one promising method for finding buried agns is the strength ( equivalent width ) of the 3.3  @xmath0 m pah emission as proposed by @xcite . pure sb galaxies show generally strong 3.3  @xmath0 m pah emission with ew@xmath60  nm because plenty of sb - produced uv photons excite pah molecules @xcite . by contrast , if an agn resides in the galaxy , the pah feature becomes weaker , ew@xmath61  nm @xcite . @xcite also used the red continuum slope ( @xmath27 for @xmath62 ) for finding buried agns . in this paper , we develop the previous works not only to find buried agns but also to quantitatively estimate their energy contribution by performing detailed spectral fitting with a continuum model consisting of a stellar component and dust components from agn and/or sb . emission and absorption line features are also included in the model . the agn emission in the near- to mid - infrared originates from the radiation reprocessed by hot dust in the torus . recent studies of dusty torus models reproduce the observed nuclear mid - infrared emission well by assuming that the torus media have a clumpy structure ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) . however , some authors mentioned that there are difficulties in fitting their near - infrared spectral energy distribution ( sed ) @xcite . this is due to the complicated degeneracies including a possible existence of an extra hot - dust ( @xmath63  k ) component originating from the vicinity of the agn @xcite , and/or contamination from the host galaxies , and/or the extinction of the torus emission by interstellar matter in the host galaxies . in this work , we take the simplest approach to approximate the torus dust emission by a single blackbody component , considering our limited wavelength range . presumably , the dusty torus has a continuous dust temperature distribution and produces emission peaked around @xmath510  @xmath0 m . hence , our single blackbody fit should trace its peak temperature at @xmath5300  k. accordingly , we set a conservative upper limit on its temperature as @xmath64  k @xcite . the stellar photospheric emission in the host galaxy contributes to the total radiation in the 2.55.0  @xmath0 m spectral range . this component produces a decreasing continuum flux at @xmath65  @xmath0 m @xcite . we first determine the temperature @xmath66 by fitting the 2mass @xmath67 and @xmath68 photometric data with a blackbody model , and fix the temperature at its best - fit value when analyzing the _ akari _ 2.55.0  @xmath0 m spectra . for mrk  551 and mrk  848 , we can not find any 2mass fluxes . therefore , we use the @xmath67 and @xmath69 photometric data from the literature @xcite instead of the 2mass data . sb produces a large amount of dust , which generally has lower temperatures than those of the agn - heated dust , typically @xmath70  k @xcite . as already suggested by @xcite and many other studies , such cold dust with @xmath71  k does not significantly contribute to the 2.55.0  @xmath0 m spectra , and therefore we ignore the contribution from the sb dust emission . however , some sb galaxies show excess emission from very hot dust with a temparature of @xmath72  k , requring at least two blackbody components . this very hot dust component of @xmath72  k would not originate from the dusty torus , from which blackbody radiation of @xmath5300  k is generally expected ( see section  4.1 . ) . one possible origin is the emission mainly heated by massive stars in hii regions @xcite . if such emission is required from the data , we model it with a single blackbody with a temperature of @xmath73  k in a range of @xmath74  k as suggested by @xcite . lc*8c eso  286-ig19 & @xmath75 & @xmath76 & @xmath77 & @xmath78 & n & y + ic  5135 & @xmath79 & @xmath80 & @xmath77 & @xmath81 & n & y + iras  03068 - 5346 & @xmath82 & @xmath83 & @xmath77 & @xmath84 & n & y + iras  03538 - 6432 & @xmath85 & @xmath86 & @xmath77 & @xmath87 & y & y + iras  10494 + 4424 & @xmath88 & @xmath89 & @xmath77 & @xmath90 & n & y + iras  17028 + 5817 & @xmath91 & @xmath92 & @xmath77 & @xmath93 & n & y + irasf  07353 + 2903 & @xmath94 & @xmath95 & @xmath77 & @xmath96 & y & y + mcg  + 02 - 04 - 025 & @xmath97 & @xmath98 & @xmath77 & @xmath99 & n & y + mcg  + 08 - 23 - 097 & @xmath100 & @xmath101 & @xmath77 & @xmath102 & n & y + mcg  + 10 - 19 - 057 & @xmath103 & @xmath104 & @xmath77 & @xmath105 & n & y + mrk  1490 & @xmath106 & @xmath107 & @xmath108 & @xmath109 & n & n + mrk  551 & @xmath110 & @xmath111 & @xmath77 & @xmath112 & n & y + mrk  848 & @xmath113 & @xmath114 & @xmath77 & @xmath115 & n & y + ngc  2339 & @xmath116 & @xmath117 & @xmath118 & @xmath109 & n & n + ngc  2388 & @xmath119 & @xmath120 & @xmath121 & @xmath109 & n & n + ngc  4102 & @xmath122 & @xmath123 & @xmath77 & @xmath124 & y & y + ngc  4194 & @xmath125 & @xmath126 & @xmath127 & @xmath109 & n & n + ngc  4818 & @xmath128 & @xmath129 & @xmath130 & @xmath109 & n & n + ngc  520 & @xmath131 & @xmath132 & @xmath133 & @xmath109 & n & n + ngc  6285 & @xmath134 & @xmath135 & @xmath136 & @xmath109 & n & n + ngc  838 & @xmath137 & @xmath138 & @xmath139 & @xmath109 & n & n + zw  453.062 & @xmath140 & @xmath141 & @xmath77 & @xmath142 & n & y +   +      the 2.55.0  @xmath0 m spectra of infrared galaxies contain various emission / absorption line features . the previous studies of the _ akari _ `` agnul '' program @xcite reported three pah lines at 3.3 , 3.4 , and 3.5  @xmath0 m , hydrogen brakett series at 4.05  @xmath0 m br@xmath143 and 2.63  @xmath0 m br@xmath144 , and pfund series at 4.65  @xmath0 m pf@xmath144 and 3.74  @xmath0 m pf@xmath145 as emission lines . they also reported absorption features at 3.1  @xmath0 m by h@xmath146o ice , at 4.29  @xmath0 m by co@xmath146 , and at 4.67  @xmath0 m by co. because pah emission is believed to originate from solid state molecules , the line profile becomes lorentzian @xcite . for simplicity , however , we model all the lines including those of pah by gaussian profiles , following @xcite . the central wavelength and line width are fixed within a spectral resolution at appropriate values reported in the literature . lc*10c eso  286-ig19 & @xmath147 & @xmath148 & @xmath149 & @xmath150 & @xmath151 & @xmath152 & @xmath153 & @xmath154 & @xmath155 & @xmath156 + ic  5135 & @xmath157 & @xmath148 & @xmath158 & @xmath159 & @xmath160 & @xmath161 & @xmath153 & @xmath162 & @xmath163 & @xmath164 + iras  03068 - 5346 & @xmath165 & @xmath148 & @xmath166 & @xmath167 & @xmath168 & @xmath169 & @xmath153 & @xmath170 & @xmath171 & @xmath172 + iras  03538 - 6432 & @xmath173 & @xmath148 & @xmath174 & @xmath175 & @xmath176 & @xmath177 & @xmath153 & @xmath178 & @xmath179 & @xmath180 + iras  10494 + 4424 & @xmath181 & @xmath148 & @xmath182 & @xmath183 & @xmath184 & @xmath185 & @xmath153 & @xmath186 & @xmath187 & @xmath188 + iras  17028 + 5817 & @xmath189 & @xmath148 & @xmath190 & @xmath191 & @xmath192 & @xmath193 & @xmath153 & @xmath194 & @xmath195 & @xmath196 + irasf  07353 + 2903 & @xmath197 & @xmath148 & @xmath198 & @xmath199 & @xmath200 & @xmath201 & @xmath153 & @xmath202 & @xmath203 & @xmath204 + mcg  + 02 - 04 - 025 & @xmath205 & @xmath148 & @xmath206 & @xmath207 & @xmath208 & @xmath209 & @xmath153 & @xmath210 & @xmath211 & @xmath212 + mcg  + 08 - 23 - 097 & @xmath213 & @xmath148 & @xmath214 & @xmath215 & @xmath216 & @xmath217 & @xmath153 & @xmath218 & @xmath219 & @xmath220 + mcg  + 10 - 19 - 057 & @xmath221 & @xmath148 & @xmath222 & @xmath223 & @xmath224 & @xmath225 & @xmath153 & @xmath226 & @xmath227 & @xmath228 + mrk  1490 & @xmath229 & @xmath230 & @xmath148 & @xmath231 & @xmath232 & @xmath233 & @xmath234 & @xmath153 & @xmath235 & @xmath77 + mrk  551 & @xmath236 & @xmath148 & @xmath237 & @xmath238 & @xmath239 & @xmath240 & @xmath153 & @xmath241 & @xmath242 & @xmath243 + mrk  848 & @xmath244 & @xmath148 & @xmath245 & @xmath246 & @xmath247 & @xmath248 & @xmath153 & @xmath249 & @xmath250 & @xmath251 + ngc  2339 & @xmath252 & @xmath253 & @xmath148 & @xmath254 & @xmath232 & @xmath255 & @xmath256 & @xmath153 & @xmath257 & @xmath77 + ngc  2388 & @xmath258 & @xmath259 & @xmath148 & @xmath260 & @xmath232 & @xmath261 & @xmath262 & @xmath153 & @xmath263 & @xmath77 + ngc  4102 & @xmath264 & @xmath148 & @xmath265 & @xmath266 & @xmath267 & @xmath268 & @xmath153 & @xmath269 & @xmath270 & @xmath271 + ngc  4194 & @xmath272 & @xmath273 & @xmath148 & @xmath274 & @xmath232 & @xmath275 & @xmath276 & @xmath153 & @xmath277 & @xmath77 + ngc  4818 & @xmath278 & @xmath279 & @xmath148 & @xmath280 & @xmath232 & @xmath281 & @xmath282 & @xmath153 & @xmath283 & @xmath77 + ngc  520 & @xmath284 & @xmath285 & @xmath148 & @xmath286 & @xmath232 & @xmath287 & @xmath288 & @xmath153 & @xmath289 & @xmath77 + ngc  6285 & @xmath259 & @xmath290 & @xmath148 & @xmath291 & @xmath232 & @xmath292 & @xmath293 & @xmath153 & @xmath294 & @xmath77 + ngc  838 & @xmath295 & @xmath296 & @xmath148 & @xmath297 & @xmath232 & @xmath298 & @xmath299 & @xmath153 & @xmath300 & @xmath77 + zw  453.062 & @xmath301 & @xmath148 & @xmath302 & @xmath303 & @xmath304 & @xmath305 & @xmath153 & @xmath306 & @xmath307 & @xmath308 +   +      we thus construct a spectral model consisting of the multiple continuum components over which the emission / absorption lines are superposed . for the continuum we consider three blackbody components , ( 1 ) a direct stellar component @xmath309 , ( 2 ) a dust component from agn with a temperature below 800  k @xmath310 , and ( 3 ) a very hot ( 800  k1200  k ) dust component from hii regions @xmath311 . each blackbody component has two parameters , normalization ( @xmath312 ) and temperature ( @xmath313 ) , which , except for @xmath314 , are left as free parameters in the spectral fit . previous studies suggest that emission from warm dust may be partly absorbed by outer cold dust ( e.g. , * ? ? ? * ; * ? ? ? * ) . to take this into account , we apply an extinction correction by incorporating @xmath315 into the blackbody components . we adopt the dust extinction index of @xmath316 @xcite as @xmath317 the normalization value @xmath318 can be estimated from the optical thickness of 3.1  @xmath0 m ice ( @xmath319 ) , which is calculated by using the deconvolved continuum from the dust ( @xmath320 ) . @xcite estimated the relationship between the dust extinction @xmath321 and @xmath319 with @xmath322 where @xmath323 is the fraction of dust that is covered with an ice mantle . the maximum value of @xmath323 is derived to be 0.3 in the core of u / lirgs in @xcite . from the relation between the optical thickness and extinction with @xmath324 and eq . ( 2 ) with @xmath325 , @xmath326 is derived . by combining eq . ( 1 ) and ( 3 ) , at the @xmath327 band ( 0.55  @xmath0 m ) , @xmath328 accordingly , we model ths observed component as @xmath329 , where @xmath330 . first , we adopt only two continuum components of direct stellar emission and agn / sb dust emission to fit the observed spectra , because the contribution of dust emission from hii regions is not always required from the data . this model ( 2 component model ) is written as : @xmath331 where @xmath332 represents each line s profile as described in section  4.3 . in some cases , mainly for sources showing little contribution from @xmath333 , the 2-component model given by equation ( 6 ) does not fit the data well . then , we add the @xmath334 component to the above model : @xmath335 we adopt this model ( 3 component model ) only when the fitting result is significantly improved from the 2-component model with a @xmath336% confidence level on the basis of an f - test . the fitting model we adopt for each source is summarized in table  1 . as shown in figure  2 , all sources that require the three components are galaxies showing blue contina in their spectra ( e.g. , ngc  4818 ) . as many authors ( e.g. , * ? ? ? * ; * ? ? ? * ) suggested , these galaxies do not show buried agn signs but are more likely normal sb galaxies . this is consistent with our assumption that the extra hot ( @xmath337  k ) component should originate from    hii regions , not from the agn torus . +    [ tbp ] figure  2 displays the observed _ akari _ spectra of all the 22 sources ( red dots with flux error bars ) overplotted with their best - fit models ( blue line ) . each line / continuum component is also plotted in the figure . as noticed , a majority of sources show 3.3  @xmath0 m pah emission features . table  2 and  3 summarize the fitting results .    for finding buried agns , we apply two agn diagnostics : ( 1 ) pah emission line and ( 2 ) torus - dust continuum . the pah diagnostic is based on ew@xmath20 as discussed in section  4 ; we identify a buried agn if it shows a small 3.3  @xmath0 m pah equivalent width with ew@xmath61  nm . the torus - dust diagnostic is based on whether or not there is a contribution from a hot dust component from the torus ; we set the criterion of identifying buried agn as @xmath338  k.    in the left panel of figure  3 , we present the two agn diagnostic diagram applied for our sample , where ulirgs / lirgs / irgs are plotted as pink circles / red triangles / brown squares , respectively . there is a clear boundary of distribution of the dust temperature around @xmath339  k. this supports our criterion that agns should have high dust temperatures ( @xmath338  k ) and are well distinguishable from sbs . interestingly , many infrared galaxies with large pah equivalent widths ( ew@xmath340  nm ) show high dust temperatures ( @xmath341  k ) . one reason could be that because of the large aperture of _ akari_/irc and less luminous nucleus emission in these galaxies , the spectra contain relatively large contribution of 3.3  @xmath0 m pah emission in the host galaxies where pah molecules are not destroyed by agn x - ray photons . this effect works to increase the pah equivalent width , and hence the pah diagnostic could miss buried agns . by contrast , the hot - dust diagnostic is not subject to the aperture effect because the continuum emission can be properly decomposed by the spectral fit . in fact , the torus - dust agn diagnostic can recover all the buried agns identified by the pah diagnostic ; in other words , there is no source with a low dust temperature ( @xmath342  k ) _ and _ a small pah equivalent width ( ew@xmath61  nm ) in our sample . this supports the superiority of using the hot - dust diagnostic for finding agns completely . +    the right panel of figure  3 shows the plot of ew@xmath20 versus dust luminosity ( @xmath343 ) tabulated in table  2 and table  3 . in the figure , we plot only the sources with @xmath344  k. this is because the 2.55.0  @xmath0 m spectral range is not sensitive to the dust emission with low temperature @xmath345  k. in this case , only an upper limit of the dust luminosity can be derived . this is the reason why the dust luminosity with @xmath346  k is very small , @xmath347  erg s@xmath35 . the figure shows that all the sources with @xmath348  k have significant dust emission with @xmath349  erg s@xmath35 . thus , the alternative criterion for detecting agns can be @xmath350  erg s@xmath35 instead of @xmath348  k . in summary , we conclude that the hot - dust agn diagnostic is a better method for finding buried agns in the _ akari _ irc band . +      one motivation of our study is to reveal how frequently buried agns exist in infrared galaxies . figure  4 shows the fraction of buried agn as a function of infrared luminosity . the left and right panels show the results based on the pah agn diagnostic and on the torus - dust agn diagnostic , respectively . as discussed above , the latter diagnostic is more complete than the former for finding buried agns from infrared galaxies . therefore , we discuss the result of the right panel hereafter .    in total , we detect 14 buried agns out of the 22 infrared galaxies without any agn signs in the optical band . we also find a clear trend that the fraction of buried agn increases with infrared luminosity . while only 17% ( 1 out of 6 sources ) of the irgs contain buried agns , the u / lirgs contain them almost ubiquitously ( @xmath351% for lirgs and @xmath352% for ulirgs ) . our results for u / lirgs are consistent with those obtained by @xcite within the statistical uncertainties . although most of the irgs are selected from the 12  @xmath0 m galaxy catalog , which is sensitive to torus - dust emission , an agn is detected only from one source out of the five targets . the high agn fraction in u / lirgs may be related with their high merger rates . there is strong observational evidence that most of u / lirgs have experienced merging or strong galaxy - galaxy interactions @xcite . the mergers can induce strong sb within the galaxies , which also make galactic gas fall into the galactic center and trigger agn activity obscured by dust @xcite . @xcite show that infrared galaxies with signs of recent merger drastically increase at @xmath353 toward a higher luminosity range . this result is in accordance with our finding that the buried agn fraction drastically increases from @xmath354 .     as a function of infrared luminosity . open circles / triangles / squares represent ulirgs / lirgs / irgs with agn signs , respectively . filled triangle / square represents lirgs / irgs without agn signs , respectively . note that all ulirgs in our sample are diagnosed to have buried agns . [ fig-5],title="fig : " ] +    we confirm the dependence of the agn fraction on luminosity by investigating the ratio between the 3.3  @xmath0 m pah luminosity ( @xmath355 ) and total infrared luminosity @xmath356 . pure sb galaxies generally have a constant ratio @xmath357 . if an agn exists inside the galaxy , the ratio becomes smaller because x - rays from the agn destroy pah molecules and the infrared continuum luminosity is increased by the torus - dust emission . figure  5 plots @xmath358 as a function of infrared luminosity for our sample . as noted , it decreases with infrared luminosity , especially at @xmath359 . we also find that the sources with small @xmath360 ratios match to those with buried agn signs , supporting the above idea . note that there is a possibility that this trend may be caused by the destruction of small pah molecules by major mergers @xcite rather than by agns . nevertheless , since all the sources with @xmath361 show buried agn signs , the conservative criterion @xmath362 could be phenomenologically used as another diagnostic for finding buried agns . our second motivation is to determine the energetic importance of agns in these infrared galaxies . as discussed in section  5.1 , the luminosity contribution from dust blackbody component ( @xmath363 ) predominantly originates from agn - heated dust , and the luminosity contribution of sb - heated dust to @xmath363 is negligible . thus , using the best - fit blackbody parameters obtained by the spectral fit , we can estimate the total luminosity from agn - heated dust and hence its ratio to the total infrared luminosity . figure  6 shows the averaged fraction of agn power to the total infrared luminosity for irgs , lirgs , and ulirgs . as noticed , the agn contribution increases with infrared luminosity , while the values are quite small , @xmath364% in irgs , @xmath365% in lirgs , and @xmath366% in ulirgs . this suggests that the bulk of the infrared emission originates from sb , not from agns in these galaxies . one caveat for this result is that we can not completely exclude the possibility that the estimated luminosity of @xmath367 might be partially originated from the torus - dust if it has an extreme temperature distribution . we thus calculate the averaged @xmath368 ratio to check if it could have a significant contribution . we find @xmath369% , suggesting that the main results above would not be affected , even if the observed very hot ( @xmath370  k ) dust component were totally due to the agn emission . another caveat for interpreting this result is the selection bias for our sample , which consists of optically non - seyfert galaxies . previous studies based on samples including optically seyfert agns found larger average agn contribution to the total infrared luminosity , 1015% in lirgs @xcite and 1540% in ulirgs @xcite . our results are in good agreement with those by @xcite within the errors , who obtained the agn to total luminosity ratio of 68% in lirgs and 1119% in ulirgs , using a similar sample to ours that preferentially includes optically non - seyfert infrared galaxies . they calculated the agn contribution by performing sed fit with the decompir package @xcite . this method requires infrared sed data covering the far - infrared band , while our method only uses the _ akari _ 2.55  @xmath0 m spectra . this supports that our simple assumption in the spectral model ( ie , single blackbody for dust emission ) works well in estimating the agn luminosity . this simple method will have a great advantage in the era of _ jwst _ , because we can apply our method to high-@xmath371 galaxies by observing the rest 2.55.0  @xmath0 m band ( 5.010.0  @xmath0 m at @xmath372 and 7.515.0  @xmath0 m at @xmath373 ) with _ jwst_/miri ( 528  @xmath0 m ) . other methods based on the _ spitzer _ bandpass ( 535  @xmath0 m ) would have a difficulty in studying high-@xmath371 objects ; for instance , the observed spectral range corresponds to 15105  @xmath0 m at @xmath373 , which only the _ spica _ mission can cover . thus , our method is achievable for studying infrared galaxies in the distant universe even before the launch of _ spica _ ( from 2025 ) . , title="fig : " ] +    combining our result with the total infrared luminosity density produced by u / lirgs in the local universe ( @xmath374  @xmath375/mpc@xmath376 and @xmath377  @xmath375/mpc@xmath376 ) @xcite , we estimate the local agn luminosity density in the infrared band ( @xmath378 ) . we find @xmath379/mpc@xmath380 for lirgs and @xmath381/mpc@xmath380 for ulirgs . these values are about an order of magnitude lower than the estimate by @xcite . this is because @xcite simply assumed that the infrared luminosity of galaxies classified as agns entirely originates from agns , while we quantitatively take into account the fraction of agn contribution to the total infrared luminosity of infrared galaxies . figure  7 shows the infrared luminosity density as a function of redshift derived by @xcite together with our results of agn the infrared luminosity density . in the era of _ jwst _ , we can fill in the agn infrared luminosity density at high redshifts thanks to their improved sensitivity in the mid - infrared band .    ) of all galaxies , pink filled circles that of ulirgs ( @xmath382 ) , and red filled triangles that of lirgs ( @xmath383 ) , taken from @xcite . pale - pink filled circles and brown filled triangles represent the agn infrared luminosity density of ulirgs and that of lirgs estimated by our work . [ fig-7],title="fig : " ] +      comparison between the starburst and agn activity gives crucial information to understand the process of galaxy - smbh co - evolution . for optically selected seyfert galaxies and qsos , many authors found that the luminosity of star formation ( @xmath384 ) and that of agn ( @xmath385 ) are well correlated with each other @xcite . we can derive a sf luminosity as @xmath386 . the agn - heated dust luminosity ( @xmath363 ) summarized in table 3 can be converted to an agn bolometric luminosity ( @xmath387 ) by using two relations given in @xcite and @xcite , respectively , @xmath388 here @xmath389 represents the intrinsic luminosity of the accretion disk integrated in the optical to x - ray band . assuming @xmath390 , we obtain @xmath391 noted that the estimated bolometric agn luminosity ( @xmath392 ) derived in equation  ( 10 ) can be overestimated . the relations from equation  ( 8) to ( 10 ) are derived from the sample of non - hidden seyfert agns , while our sources are infrared galaxies with highly obscured ( = hidden ) agns . therefore , our sources in this study could have higher covering fraction of dusty - torus than those of seyfert agns . this means the ratio of mid - infrared over bolometric luminosity in our sources could be larger than those of seyfert agns . figure  8 shows the @xmath384 versus @xmath392 relations obtained from our sample . the average luminosities of u / lirg are also plotted in the same figure . the dotted line corresponds to the linear correlation obtained from optically selected seyfert galaxies obtained by @xcite . interestingly , all the sources except two ( iras f07353 + 2903 and eso  286-ig19 ) are located above the seyfert line . this result may reflect the difference in the evolutional stage of galaxies , from a pure sb phase ( u / lirgs ) to an unburied agn phase ( seyferts ) . as discussed in section 5.2 , u / lirgs have recently experienced or are facing mergers , which result in strong sb activity from the epoch @xmath393 to @xmath394 in figure  8 . the sb also triggers agn activity , which continues from @xmath394 to @xmath395 . finally , the sb activity is gradually weakened due to the shortage of gas and feedback from the agn , and the supply of matter from the host galaxy onto the central engine also decreases , from @xmath395 to @xmath396 . in this framework of galaxy evolution , the evolutional stage of ulirgs should be a later phase of mergers than that of lirgs . the picture is in good agreement with observations that the morphology of ulirgs are more compact than that of lirgs , suggesting that ulirgs have just finished the merger process while that of lirgs is still on - going . this scenario suggests that infrared galaxies containing buried agns could be an earlier evolutional stage of agn , which will be evolved to normal ( unburied ) seyferts / qsos . this scenario becomes much more reasonable if our bolometric agn luminosity ( @xmath392 ) is overestimated as discussed in the previous paragraph because all the points are shifted to the left in the figure . while our sample size is still small , the new agn diagnostics developed in this paper can be applied to other sources . further studies using a well - defined complete sample of infrared galaxies should be urged to complete our view of agn evolution . we obtained the _ akari _ 2.55.0  @xmath0 m spectra of 22 infrared galaxies at @xmath397 selected mainly from three infrared catalogs : bright ulirg catalog @xcite , _ iras _ 12  @xmath0 m galaxy catalog @xcite , and _ iras _ bright galaxy catalog @xcite . this band includes unique spectral features for finding `` buried '' agns , the 3.3  @xmath0 m pah emission line and red continuum spectra originated from agn torus emission . we perform detailed spectral fitting to these data by properly modeling the continuum and emission / absorption features . this enables us to decompose the continuum spectra into sb and agn torus - dust components . the results are summarized as follows .    1 .   using the agn diagnostics based on the torus - dust continuum ( @xmath398  k or @xmath399  erg s@xmath35 ) and 3.3  @xmath0 m pah emission strength ( ew@xmath61  nm ) , we find that 14 out of the 22 infrared galaxies have buried agns . we confirm the trend that the buried agn fraction increases with infrared luminosity : @xmath400% for normal irgs , @xmath401% for lirgs , and @xmath402% for ulirgs in our sample . this suggests that the presence buried agns is ubiquitous in u / lirgs . the ratio @xmath358 decreases with infrared luminosity . this is opposed to the previous reports that @xmath403 have a constant value @xcite . this decrease is possibly due to three effects : ( 1 ) destruction of pah molecules by x - rays emitted from the agn , ( 2 ) pah destruction by merger processes , and ( 3 ) increase of agn contribution to the total infrared luminosity . the energy contribution from the agn torus - dust emission to the total infrared luminosity also increases with infrared luminosity , but it only reaches @xmath57% in lirgs and @xmath520% even in ulirgs . this suggests that majority of the total infrared luminosity in u / lirgs originates from sb , _ not _ from agn . 4 .   combining the above results with the luminosity function of infrared galaxies derived by @xcite , we estimate the agn luminosity density of @xmath404/mpc@xmath380 for lirgs and @xmath405/mpc@xmath380 for ulirgs in the local universe . @xmath384 versus @xmath392 luminosity plot for our sample shows that infrared galaxies with buried agns are located above the luminosity correlation line derived from optically - selected seyfert galaxies . we interpret that infrared galaxies could be an early phase of the evolutional track of agns . we are grateful to kenichi yano and rika yamada for their help with the data reduction of _ we also thank claudio ricci and kenta matsuoka for discussions . the research is based on observations with _ akari _ , a jaxa project with the participation of esa . this work was partly supported by the grant - in - aid for jsps fellows for young researchers ( k.i . ) and for scientific research 23540273 ( m.i . ) and 26400228 ( y.u . ) .
Abstract:

Answer: we present the result of a systematic infrared 2.55  @xmath0 m spectroscopic study of 22 nearby infrared galaxies over a wide infrared luminosity range ( @xmath1 ) obtained from _ akari _ infrared camera ( irc ) . the unique band of the _ akari _ irc spectroscopy enables us to access both the 3.3  @xmath0 m polycyclic aromatic hydrocarbon ( pah ) emission feature from star forming activity and the continuum of torus - dust emission heated by an active galactic nucleus ( agn ) . applying our agn diagnostics to the _ akari _ spectra , we discover 14 buried agns . the large fraction of buried agns suggests that agn activity behind the dust is almost ubiquitous in ultra-/luminous infrared galaxies ( u / lirgs ) . we also find that both the fraction and energy contribution of buried agns increase with infrared luminosity from @xmath2 to @xmath3 , including normal infrared galaxies with @xmath4 . the energy contribution from agns in the total infrared luminosity is only @xmath5 7% in lirgs and @xmath5  20% in ulirgs , suggesting that the majority of the infrared luminosity originates from starburst activity . using the pah emission , we investigate the luminosity relation between star formation and agn . we find that these infrared galaxies exhibit higher star formation rates than optically selected seyfert galaxies with the same agn luminosities , implying that infrared galaxies could be an early evolutionary phase of agn .