vector auto
regressive ( var ) models play an important role in contemporary macro  economics , being an example of an approach called the `` dynamic stochastic general equilibrium '' ( dsge ) , which is superseding traditional large  scale macro  econometric forecasting methodologies  @xcite .
the motivation behind them is based on the assertion that more recent values of a variable are more likely to contain useful information about its future movements than the older ones . on the other hand , a standard tool in multivariate time series analysis is vector moving average ( vma ) models , which is really a linear regression of the present value of the time series w.r.t .
the past values of a white noise . a broader class of stochastic processes used in macro
economy comprises both these kinds together in the form of vector auto
regressive moving average ( varma ) models
. these methodologies can capture certain spatial and temporal structures of multidimensional variables which are often neglected in practice ; including them not only results in more accurate estimation , but also leads to models which are more interpretable .
they are widely used by academia and central banks ( _ cf .
_ the european central bank s smets
wouters model for the euro zone  @xcite ) , as they constitute quite a simple version of the dsge equations .
varma models are constructed from a number of univariate arma ( box  jenkins ; see for example  @xcite ) processes , typically coupled with each other . in this paper
, we investigate only a significantly simplified circumstance when there is no coupling between the many arma components .
one may argue that this is too far fetched and will be of no use in describing an economic reality
. however , one may also treat it as a `` zeroth  order hypothesis , '' analogously to the idea of  @xcite in finance , namely that the case with no cross  covariances is considered theoretically , and subsequently compared to some real  world data modeled by a varma process ; any discrepancy between the two will reflect nontrivial cross
covariances present in the system , thus permitting their investigation .
this latter route is taken in this communication .
a challenging and yet increasingly important problem is the estimation of large covariance matrices generated by these stationary processes , since high dimensionality of the data as compared to the sample size is quite common in many statistical problems ( the `` dimensionality curse '' ) .
therefore , an appropriate `` noise cleaning '' procedure has to be implemented , and random matrix theory ( rmt ) provides a natural and efficient outfit for doing that .
in particular , the mean spectral densities ( a.k.a .
`` limiting spectral distributions , '' lsd ) of the pearson estimators of the cross
covariances for the @xmath7 and @xmath8 models , in the relevant high  dimensionality sector and under the full decoupling , have been derived in  @xcite by applying the framework proposed by  @xcite .    in this paper , we suggest that such calculations can be considerably simplified by resorting to a mathematical concept of the free random variables ( frv ) calculus  @xcite , succinctly introduced in sec .  [ s : doublycorrelatedwishartensemblesandfreerandomvariables ] .
our general frv formula  @xcite allows not only to rediscover , which much less strain , the two fourth  order polynomial equations obtained in  @xcite in the @xmath7 and @xmath8 cases , but also to derive a sixth  order equation ( [ eq : varmaoneonemainequation ] ) which produces the mean spectral density for a more involved @xmath6 model .
the results are verified by numerical simulations , which show a perfect agreement .
this is all done in sec .
[ s : varmafromfreerandomvariables ] .
@xmath2 stochastic processes , as we will see below , fall within quite a general set  up encountered in many areas of science where a probabilistic nature of multiple degrees of freedom evolving in time is relevant , for example , multivariate time series analysis in finance , applied macro  econometrics and engineering . to describe this framework ,
consider a situation of @xmath3 time
dependent random variables which are measured at @xmath4 consecutive time moments ( separated by some time interval @xmath9 ) ; let be the value od the @xmath10th ( @xmath11 ) random number at the @xmath12th time moment ( @xmath13 ) ; together , they make up a rectangular @xmath14 matrix @xmath15 . in what usually would be the first approximation ,
each is supposed to be drawn from a gaussian probability distribution .
we will also assume that they have mean values zero , .
these degrees of freedom may in principle display mutual correlations .
a set of correlated zero
mean gaussian numbers is fully characterized by the two  point covariance function , if the underlying stochastic process generating these numbers is stationary .
linear stochastic processes , including @xmath2 , belong to this category .
we will restrict our attention to an even narrower class where the cross  correlations between different variables and the auto  correlations between different time moments are factorized , _
i.e. _ , @xmath16 in this setting , the inter
variable covariances do not change in time ( and are described by an @xmath17 cross  covariance matrix @xmath18 ) , and also the temporal covariances are identical for all the numbers ( and are included in a @xmath19 auto  covariance matrix @xmath20 ; both these matrices are symmetric and positive  definite ) . the gaussian probability measure with this structure of covariances
is known from textbooks , @xmath21_{i j } y_{j b } \left [ \mathbf{a}^{- 1 } \right]_{b a } \right ) { \mathrm{d}}\mathbf{y } = \ ] ] @xmath22 where the normalization constant , and the integration measure , while the letters `` c.g . ''
stand for `` correlated gaussian . ''
now , a standard way to approach correlated gaussian random numbers is to recall that they can always be decomposed as linear combinations of uncorrelated gaussian degrees of freedom ; indeed , this is achieved through the transformation @xmath23 where the square roots of the covariance matrices , necessary to facilitate the transition , exist due to the positive  definiteness of @xmath18 and @xmath20 ; the new normalization reads .      an essential problem in multivariate analysis is to determine ( estimate ) the covariance matrices @xmath18 and @xmath20 from given @xmath3 time series of length @xmath4 of the realizations of our random variables . for simplicity
, we do not distinguish in notation between random numbers , _
i.e. _ , the population , and their realizations in actual experiments , _ i.e. _ , the sample .
since the realized cross  covariance between degrees @xmath10 and @xmath24 at the same time @xmath12 is , the simplest method to estimate the today s cross  covariance is to compute the time average , @xmath25 this is usually named the `` pearson estimator '' , up to the prefactor which depending on the context is @xmath26 or @xmath27 .
other estimators might be introduced , such as between distinct degrees of freedom at separate time moments ( `` time  delayed estimators '' ) , or with certain decreasing weights given to older measurements to reflect their growing obsolescence ( `` weighted estimators '' ) , but we will not investigate them in this article .
furthermore , in the last equality in ( [ eq : estimatorcdefinition ] ) , we cast @xmath28 through the uncorrelated gaussian numbers contained in , the price to pay for this being that the covariance matrices now enter into the expression for @xmath28 , making it more complicated ; this will be the form used hereafter .
the random matrix @xmath28 is called a `` doubly correlated wishart ensemble ''  @xcite .
let us also mention that the auto  covariance matrix @xmath20 can be estimated through .
however , it is verified that this object carries identical information to the one contained in @xmath28 ( it is `` dual '' to @xmath28 ) , and therefore may safely be discarded . indeed , these two estimators have same non  zero eigenvalues ( modulo an overall rescaling by @xmath29 ) , and the larger one has @xmath30 additional zero modes .
any historical estimator is inevitably marred by the measurement noise ; it will reflect the true covariances only to a certain degree , with a superimposed broadening due to the finiteness of the time series .
more precisely , there are @xmath31 independent elements in @xmath18 , to be estimated from @xmath32 measured quantities @xmath15 , hence the estimation accuracy will depend on the `` rectangularity ratio , '' @xmath33 the closer @xmath29 to zero , the more truthful the estimate .
this is a cornerstone of classical multivariate analysis .
unfortunately , a practical situation will typically feature a large number of variables sampled over a comparably big number of time snapshots , so that we may approximately talk about the `` thermodynamical limit , '' @xmath34 on the other hand , it is exactly this limit in which the frv calculus ( see the subsection below for its brief elucidation ) can be applied ; hence , the challenge of de  noising is somewhat counterbalanced by the computationally powerful frv techniques .
any study of a ( real symmetric @xmath36 ) random matrix @xmath37 will most surely include a fundamental question about the average values of its ( real ) eigenvalues .
they are concisely encoded in the `` mean spectral density , '' @xmath38 here the expectation map @xmath39 is understood to be taken w.r.t .
the probability measure @xmath40 of the random matrix .
we will always have this distribution rotationally ( _ i.e. _ , , with @xmath41 orthogonal ) invariant , and hence the full information about @xmath37 resides in its eigenvalues , distributed on average according to ( [ eq : meanspectraldensitydefinition ] ) .    on the practical side , it is more convenient to work with either of the two equivalent objects , @xmath42 referred to as the `` green s function '' ( or the `` resolvent '' ) and the `` @xmath35transform '' of @xmath37 .
the latter is also called the `` moments generating function , '' since if the `` moments '' of @xmath37 exist , it can be expanded into a power series around @xmath43 as .
it should however be underlined that even for probability measures disallowing such an expansion ( heavy  tailed distributions , preeminent in finance , being an example ) , the quantities ( [ eq : greenfunctionandmtransformdefinition ] ) still manage to entirely capture the spectral properties of @xmath37 ; hence the name `` @xmath35transform '' more appropriate , in addition to being more compact .
we will show that for our purposes ( multiplication of random matrices ; see par .
[ sss : thentransformandfreerandomvariables ] ) the @xmath35transform serves better than the green s function .
however , it is customary to write the relationship between ( [ eq : meanspectraldensitydefinition ] ) and ( [ eq : greenfunctionandmtransformdefinition ] ) in terms of this latter , @xmath44 resulting from a well  known formula for generalized functions , .
the doubly correlated wishart ensemble @xmath28 ( [ eq : estimatorcdefinition ] ) may be viewed as a product of several random and non  random matrices .
the general problem of multiplying random matrices seems formidable .
in classical probability theory , it can be effectively handled in the special situation when the random terms are independent : then , the exponential map reduces it to the addition problem of independent random numbers , solved by considering the logarithm of the characteristic functions of the respective pdfs , which proves to be additive . in matrix probability theory ,
a crucial insight came from d.  voiculescu and coworkers and r.  speicher  @xcite , who showed how to parallel the commutative construction in the noncommutative world .
it starts with the notion of `` freeness , '' which basically comprises probabilistic independence together with a lack of any directional correlation between two random matrices .
this nontrivial new property happens to be the right extension of classical independence , as it allows for an efficient algorithm of multiplying free random variables ( frv ) , which we state below :    step 1 : : :    suppose we have two random matrices , and , mutually free .
their    spectral properties are best wrought into the    @xmath35transforms    ( [ eq : greenfunctionandmtransformdefinition ] ) , and .
step 2 : : :    the critical maneuver is to turn attention to the functional inverses    of these @xmath35transforms , the so  called    `` @xmath3transforms , ''    @xmath45 step 3 : : :    the @xmath3transforms submit to a very straightforward rule    upon multiplying free random matrices ( the `` frv multiplication    law '' ) , @xmath46 step 4 : : :    finally , it remains to functionally invert the resulting    @xmath3transform to gain the @xmath35transform of the    product , , and consequently , all the spectral properties via    formula  ( [ eq : sokhotskyformula ] ) .
it is stunning that such a simple prescription ( relying on the choice of the @xmath35transform as the carrier of the mean spectral information , and the construction of its functional inverse , the @xmath3transform , which essentially multiplies under taking the free product ) resolves the multiplication problem for free random noncommutative objects .
let us just mention that the addition problem may be tackled along similar lines : in this case , the green s function should be exploited , its functional inverse considered (; it is sometimes called the `` blue s function ''  @xcite ) , which obeys the `` frv addition law , '' , for two free random matrices . in this paper , we do not resort to using this addition formula , even though our problem could be approached through it as well .
let us also remark that in the original mathematical formulations  @xcite of these frames , a slightly different language is employed : instead of the @xmath3transform , the `` @xmath47transform '' is found convenient , , while in place of the blue s function , one engages the `` @xmath48transform , '' .
they fulfil simpler laws , and , respectively .      the innate potential of the frv multiplication algorithm ( [ eq : frvmultiplicationlaw ] )
is surely revealed when inspecting the doubly correlated wishart random matrix ( [ eq : estimatorcdefinition ] ) .
this has been done in detain in  @xcite , so we will only accentuate the main results here , referring the reader to the original paper for a thorough explanation .
the idea is that one uses twice the cyclic property of the trace ( which permits cyclic shifts in the order of the terms ) , and twice the frv multiplication law ( [ eq : frvmultiplicationlaw ] ) ( to break the @xmath3transforms of products of matrices down to their constituents ) , in order to reduce the problem to solving the uncorrelated wishart ensemble .
this last model is further simplified , again by the cyclic property and the frv multiplication rule applied once , to the standard @xmath49 random matrix squared ( and the projector , designed to chip the rectangle off the square @xmath49 ) , whose properties are firmly established .
let us sketch the derivation , @xmath50 @xmath51 @xmath52    this is the basic formula .
since the spectral properties of @xmath28 are given by its @xmath35transform , , it is more pedagogical to recast ( [ eq : doublycorrelatedwishartensemblethemainequation1 ] ) as an equation for the unknown @xmath35 , @xmath53 it provides a means for computing the mean spectral density of a doubly correlated wishart random matrix once the `` true '' covariance matrices @xmath18 and @xmath20 are given .    in this communication ,
only a particular instance of this fundamental formula is applied , namely with an arbitrary auto  covariance matrix @xmath20 , but with trivial cross
covariances , . using that , equation ( [ eq : doublycorrelatedwishartensemblethemainequation2 ] ) thins out to @xmath54 which will be strongly exploited below .
let us mention that these equalities ( [ eq : doublycorrelatedwishartensemblethemainequation2 ] ) , ( [ eq : doublycorrelatedwishartensemblethemainequation3 ] ) have been derived through other , more tedious , techniques ( the planar feynman
diagrammatic expansion , the replica trick ) in  @xcite .
in what follows , we will assume that the @xmath0 , @xmath1 , or @xmath2 stochastic processes are covariance ( weak ) stationary ; for details , we refer to  @xcite .
it implies certain restrictions on their parameters , but we will not bother with this issue in the current work .
another consequence is that the processes display some interesting features , such as invertibility .
for all this , we must in particular take both @xmath3 and @xmath4 large from the start , with their ratio @xmath55 fixed ( [ eq : thermodynamicallimit ] ) .
more precisely , we stretch the range of the @xmath12index from minus to plus infinity .
this means that all the finite  size effects ( appearing at the ends of the time series ) are readily disregarded .
in particular , there is no need to care about initial conditions for the processes , and all the recurrence relations are assumed to continue to the infinite past .
we consider a situation when @xmath3 stochastic variables evolve according to identical independent @xmath0 ( vector moving average ) processes , which we sample over a time span of @xmath4 moments .
this is a simple generalization of the standard univariate weak  stationary moving average @xmath56 . in such a setting , the value of the @xmath10th ( @xmath11 ) random variable at time moment @xmath12 ( @xmath13 ) can be expressed as @xmath57 here all the s are iid standard ( mean zero , variance one ) gaussian random numbers ( white noise ) , .
the s are some @xmath58 real constants ; importantly , they do not depend on the index @xmath10 , which reflects the fact that the processes are identical and independent ( no `` spatial '' covariances among the variables ) .
the rank @xmath59 of the process is a positive integer .      in order to handle such a process ( [ eq : vmaqdefinition ] ) , notice that the s , being linear combinations of uncorrelated gaussian numbers , must also be gaussian random variables , albeit displaying some correlations
therefore , to fully characterize these variables , it is sufficient to calculate their two  point covariance function ; this is straightforwardly done ( see appendix  [ ss : theautocovariancematrixforvmaq ] for details ) , @xmath60 where @xmath61 in other words , the cross  covariance matrix is trivial , ( no correlations between different variables ) , while the auto  covariance matrix , responsible for temporal correlations , can be called `` @xmath62diagonal . ''
in the course of this article , we will use several different auto  covariance matrices , and for brevity , we decide to label them with superscripts ; their definitions are all collected in appendix  [ ss : alistofthevariousautocovariancematricesused ] .
for example , in the simplest case of @xmath7 , it is tri  diagonal , @xmath63      such an infinite matrix ( [ eq : vmaqautocovariancematrix ] ) is translationally invariant ( as announced , it is one of the implications of the weak stationarity ) , _ i.e. _ , the value of any of its entries depends only on the distance between its indices , ; specifically , , for @xmath64 , and .
hence , it is convenient to rewrite this matrix in the fourier space , @xmath65    in this representation , the @xmath35transform of is readily obtained  @xcite , @xmath66 this integral can be evaluated by the method of residues for any value of @xmath59 , which we do in appendix  [ ss : themtransformoftheautocovariancematrixforvmaq ] , where also we print the general result ( [ eq : vmaqgreenfunctionoftheautocovariancematrix ] ) . in particular , for @xmath67 , @xmath68 where the square roots are principal .      we will be interested in investigating the spectral properties of the pearson estimator ( [ eq : estimatorcdefinition ] ) . the @xmath35transform of this correlated wishart random matrix , , can be retrieved from equation ( [ eq : doublycorrelatedwishartensemblethemainequation3 ] ) .
we could write it for any @xmath59 using ( [ eq : vmaqgreenfunctionoftheautocovariancematrix ] ) , but we will restrict ourselves only to @xmath67 , in which case the substitution of ( [ eq : vmaonemomentsgeneratingfunctionfora ] ) leads to a fourth  order polynomial ( ferrari ) equation for the unknown @xmath35 ,
@xmath69 @xmath70 @xmath71 the frv technique allowed us therefore to find this equation in a matter of a few lines of a simple algebraic computation .
it has already been derived in  @xcite , and ( [ eq : vmaonemainequation ] ) may be verified to coincide with the version given in that paper . in  @xcite
, the pertinent equation is printed before ( a.6 ) , and to compare the two , one needs to change their variables into ours according to @xmath72 , @xmath73 , and .
the last equality means that and @xmath74 of  @xcite correspond in our language to the green s functions and , respectively , where is the pearson estimator dual to @xmath28 .
as mentioned , a quick extension to the case of arbitrary @xmath59 is possible , however the resulting equations for @xmath35 will be significantly more complicated ; for instance , for @xmath75 , a lengthy ninth  order polynomial equation is discovered .        a set  up of @xmath3 identical and independent @xmath1 ( vector auto  regressive ) processes
is somewhat akin to ( [ eq : vmaqdefinition ] ) , _ i.e. _ , we consider @xmath3 decoupled copies of a standard univariate @xmath76 process , @xmath77 it is again described by the demeaned and standardized gaussian white noise ( which triggers the stochastic evolution ) , as well as @xmath58 real constants , , with @xmath78 . as announced before , the time stretches to the past infinity , so no initial condition is necessary .
although at first sight ( [ eq : varqdefinition ] ) may appear to be a more involved recurrence relation for the s , it is actually easily reduced to the @xmath0 case : it remains to remark that if one exchanges the s with the s , one precisely arrives at the @xmath0 process with the constants , , @xmath78 .
in other words , the auto  covariance matrix of the @xmath1 process ( [ eq : varqdefinition ] ) is simply the inverse of the auto  covariance matrix of the corresponding @xmath0 process with the described modification of the parameters , @xmath79 this inverse exists thanks to the weak stationarity supposition .
the fourier transform of the auto  covariance matrix of @xmath1 is therefore a ( number ) inverse of its counterpart for @xmath0 with its parameters appropriately changed , @xmath80 where @xmath81 and where we define .    in order to find the @xmath35transform of the inverse matrix , , one employs a general result , true for any ( real symmetric ) random matrix @xmath37 , and obtainable through an easy algebra , @xmath82 since the quantity is known for any @xmath59 ( [ eq : vmaqgreenfunctionoftheautocovariancematrix ] ) , hence is via ( [ eq : mhinversefrommh ] ) , but we will not print it explicitly .
let us just give it for @xmath67 , in which case ( [ eq : mhinversefrommh ] ) and ( [ eq : vmaonemomentsgeneratingfunctionfora ] ) yield @xmath83      despite being somewhat outside of the main line of thought of this article , an interesting question would be to search for an explicit expression for the auto  covariance matrix from its fourier transform ( [ eq : varqautocovariancematrixfourier ] ) , @xmath84 where we exploited the fact that must be translationally invariant , .
this computation would shed light on the structure of temporal correlations present in a var setting .
this integral is evaluated by the method of residues in a very similar manner to the one shown in appendix  [ ss : themtransformoftheautocovariancematrixforvmaq ] , and we do this in appendix  [ ss : theautocovariancematrixforvarq ] .
we discover that the auto  covariance matrix is a sum of @xmath59 exponential decays , @xmath85 where are constants , and are the characteristic times ( [ eq : varqautocovariancematrixcharacteristictimes ] ) , @xmath86 ; these constituents are given explicitly in ( [ eq : varqautocovariancematrix2 ] ) .
this is a well  known fact , nevertheless we wanted to establish it again within our approach .
for example , for @xmath67 , the auto  covariance matrix of @xmath8 is one exponential decay , @xmath87 where we assumed for simplicity ( the formula can be easily extended to all values of ) .
having found an expression for the @xmath35transform of the auto  covariance matrix of a @xmath1 ( [ eq : mhinversefrommh ] ) , ( [ eq : vmaqgreenfunctionoftheautocovariancematrix ] ) , we may proceed to investigate the equation ( [ eq : doublycorrelatedwishartensemblethemainequation3 ] ) for the @xmath35transform of the correlated wishart random matrix ( [ eq : estimatorcdefinition ] ) .
we will do this explicitly only for @xmath67 , when ( [ eq : varonemomentsgeneratingfunctionfora ] ) leads to a fourth  order ( ferrari ) polynomial equation for the unknown @xmath35 , @xmath88 @xmath89 this equation has been derived by another method in  @xcite , and our result confirms their equation ( a.8 ) , with the change in notation , @xmath72 , @xmath73 , @xmath90 .
the two types of processes which we elaborated on above , and , can be combined into one stochastic process called , @xmath91 now it is a straightforward and well  known observation ( which can be verified by a direct calculation ) that the auto  covariance matrix of this process is simply the product ( in any order ) of the auto  covariance matrices of the var and vma pieces ; more precisely , @xmath92 where corresponds to the generic model ( [ eq : vmaqautocovariancematrix ] ) , while denotes the auto  covariance matrix of with a slightly different modification of the parameters compared to the previously used , namely , , for .
we have thus already made use here of the fact that the auto  covariance matrix of a var process is the inverse of the auto  covariance matrix of a certain corresponding vma process ( [ eq : varqfromvmaq ] ) , but the new change in parameters necessary in moving from var to vma has effectively w.r.t .
what we had before ( [ eq : varqfromvmaq ] ) ; it is understandable : this `` missing '' is now included in the matrix of the other process .     of the cross  covariances in the @xmath6 process computed numerically from the sixth  order polynomial equation ( [ eq : varmaoneonemainequation ] ) , for various values of the process parameters .
the scale of these parameters is determined by choosing everywhere .
recall that the theoretical formula ( [ eq : varmaoneonemainequation ] ) is valid in the thermodynamical limit ( [ eq : thermodynamicallimit ] ) of @xmath93 , with @xmath94 kept finite . + up left : we set the remaining varma parameters to , , while the rectangularity ratio takes the values @xmath95 ( the purple line ) , @xmath96 ( red ) , @xmath97 ( magenta ) , @xmath98 ( pink ) ; each one is @xmath99 times smaller than the preceding one .
we observe how the graphs become increasingly peaked ( narrower and taller ) around @xmath100 as @xmath29 decreases , which reflects the movement of the estimator @xmath28 toward its underlying value .
+ up right : we fix @xmath101 and consider the two varma parameters equal to each other , with the values ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath104 ( pink ) .
+ down left : we hold @xmath101 and , and modify ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; for this last value , the @xmath6 model reduces to @xmath8 . + down right : similarly , but this time we assign @xmath101 and , while changing ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; this last value corresponds to @xmath7.,title="fig:",width=321 ]   of the cross  covariances in the @xmath6 process computed numerically from the sixth  order polynomial equation ( [ eq : varmaoneonemainequation ] ) , for various values of the process parameters .
the scale of these parameters is determined by choosing everywhere . recall that the theoretical formula ( [ eq : varmaoneonemainequation ] ) is valid in the thermodynamical limit ( [ eq : thermodynamicallimit ] ) of @xmath93 , with @xmath94 kept finite . + up left : we set the remaining varma parameters to , , while the rectangularity ratio takes the values @xmath95 ( the purple line ) , @xmath96 ( red ) , @xmath97 ( magenta ) , @xmath98 ( pink ) ; each one is @xmath99 times smaller than the preceding one .
we observe how the graphs become increasingly peaked ( narrower and taller ) around @xmath100 as @xmath29 decreases , which reflects the movement of the estimator @xmath28 toward its underlying value .
+ up right : we fix @xmath101 and consider the two varma parameters equal to each other , with the values ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath104 ( pink ) .
+ down left : we hold @xmath101 and , and modify ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; for this last value , the @xmath6 model reduces to @xmath8 .
+ down right : similarly , but this time we assign @xmath101 and , while changing ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; this last value corresponds to @xmath7.,title="fig:",width=321 ]   of the cross  covariances in the @xmath6 process computed numerically from the sixth  order polynomial equation ( [ eq : varmaoneonemainequation ] ) , for various values of the process parameters .
the scale of these parameters is determined by choosing everywhere .
recall that the theoretical formula ( [ eq : varmaoneonemainequation ] ) is valid in the thermodynamical limit ( [ eq : thermodynamicallimit ] ) of @xmath93 , with @xmath94 kept finite . + up left : we set the remaining varma parameters to , , while the rectangularity ratio takes the values @xmath95 ( the purple line ) , @xmath96 ( red ) , @xmath97 ( magenta ) , @xmath98 ( pink ) ; each one is @xmath99 times smaller than the preceding one .
we observe how the graphs become increasingly peaked ( narrower and taller ) around @xmath100 as @xmath29 decreases , which reflects the movement of the estimator @xmath28 toward its underlying value .
+ up right : we fix @xmath101 and consider the two varma parameters equal to each other , with the values ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath104 ( pink ) .
+ down left : we hold @xmath101 and , and modify ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; for this last value , the @xmath6 model reduces to @xmath8 . + down right : similarly , but this time we assign @xmath101 and , while changing ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; this last value corresponds to @xmath7.,title="fig:",width=321 ]   of the cross  covariances in the @xmath6 process computed numerically from the sixth  order polynomial equation ( [ eq : varmaoneonemainequation ] ) , for various values of the process parameters .
the scale of these parameters is determined by choosing everywhere .
recall that the theoretical formula ( [ eq : varmaoneonemainequation ] ) is valid in the thermodynamical limit ( [ eq : thermodynamicallimit ] ) of @xmath93 , with @xmath94 kept finite . + up left : we set the remaining varma parameters to , , while the rectangularity ratio takes the values @xmath95 ( the purple line ) , @xmath96 ( red ) , @xmath97 ( magenta ) , @xmath98 ( pink ) ; each one is @xmath99 times smaller than the preceding one .
we observe how the graphs become increasingly peaked ( narrower and taller ) around @xmath100 as @xmath29 decreases , which reflects the movement of the estimator @xmath28 toward its underlying value .
+ up right : we fix @xmath101 and consider the two varma parameters equal to each other , with the values ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath104 ( pink ) .
+ down left : we hold @xmath101 and , and modify ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; for this last value , the @xmath6 model reduces to @xmath8 . + down right : similarly , but this time we assign @xmath101 and , while changing ( purple ) , @xmath102 ( red ) , @xmath103 ( magenta ) , @xmath105 ( pink ) ; this last value corresponds to @xmath7.,title="fig:",width=321 ]    ) ( the dashed red lines ) .
the conformity is nearly perfect .
we generate the matrices @xmath15 of sizes @xmath106 , @xmath107 ( _ i.e. _ , @xmath101 ) from the @xmath6 process with the parameters , ,
. the monte carlo simulation is repeated @xmath108 ( left ) or @xmath109 ( right ) times ; in this latter case , a significant improvement in the quality of the agreement is seen .
one notices finite  size effects at the edges of the spectrum ( `` leaking out '' of eigenvalues ) : in the numerical simulations , @xmath3 and @xmath4 are obviously finite , while equation ( [ eq : varmaoneonemainequation ] ) is legitimate in the thermodynamical limit ( [ eq : thermodynamicallimit ] ) only , hence the small discrepancies ; by enlarging the chosen dimensions @xmath110 one would diminish this fallout.,title="fig:",width=321 ] ) ( the dashed red lines ) .
the conformity is nearly perfect .
we generate the matrices @xmath15 of sizes @xmath106 , @xmath107 ( _ i.e. _ , @xmath101 ) from the @xmath6 process with the parameters , ,
. the monte carlo simulation is repeated @xmath108 ( left ) or @xmath109 ( right ) times ; in this latter case , a significant improvement in the quality of the agreement is seen .
one notices finite  size effects at the edges of the spectrum ( `` leaking out '' of eigenvalues ) : in the numerical simulations , @xmath3 and @xmath4 are obviously finite , while equation ( [ eq : varmaoneonemainequation ] ) is legitimate in the thermodynamical limit ( [ eq : thermodynamicallimit ] ) only , hence the small discrepancies ; by enlarging the chosen dimensions @xmath110 one would diminish this fallout.,title="fig:",width=321 ]      the fourier transform of the auto  covariance matrix of ( [ eq : varmaq1q2fromvarq1vmaq2 ] ) is simply the product of the respective fourier transforms ( [ eq : vmaqautocovariancematrixfourier ] ) and ( [ eq : varqautocovariancematrixfourier ] ) , @xmath111 where @xmath112 where we recall . for instance , for @xmath6
( it is described by three constants , , , ) , one explicitly has @xmath113    the @xmath35transform of can consequently be derived from the general formula ( [ eq : momentsgeneratingfunctionforatranslationallyinvariant ] ) .
we will evaluate here the pertinent integral only for the simplest @xmath6 process , even though an arbitrary case may be handled by the technique of residues , @xmath114      one might again attempt to track the structure of temporal covariances in a varma process .
this can be done either by the inverse fourier transform of ( [ eq : varmaq1q2autocovariancematrixfourier ] ) , or through a direct computation based on the recurrence relation ( [ eq : varmaq1q2definition ] ) ( importantly , adhering to the assumption that it stretches to the past infinity ) .
let us print the result just for @xmath6 , @xmath115 where for simplicity .
this is an exponential decay , with the characteristic time of the var piece , with an additional term on the diagonal .
expression ( [ eq : varmaoneonemomentsgeneratingfunctionforafive ] ) , along with the fundamental frv formula ( [ eq : doublycorrelatedwishartensemblethemainequation3 ] ) , allow us to write the equation satisfied by the @xmath35transform of the pearson estimator ( [ eq : estimatorcdefinition ] ) of the cross  covariances in the @xmath6 process ; it happens to be polynomial of order six , and we print it ( [ eq : varmaoneonemainequation ] ) in appendix  [ ss : theequationforthemtransformofthepearsonestimatorofthecovariancesforvarmaoneone ] .
it may be solved numerically , a proper solution chosen ( the one which leads to a sensible density : real , positive  definite , normalized to unity ) , and finally , the mean spectral density derived from ( [ eq : sokhotskyformula ] ) .
we show the shapes of this density for a variety of the values of the parameters @xmath29 , , , in fig .
[ fig : varmaoneonetheory ] .
moreover , in order to test the result ( [ eq : varmaoneonemainequation ] ) , and more broadly , to further establish our frv framework in the guise of formula ( [ eq : doublycorrelatedwishartensemblethemainequation3 ] ) , the theoretical form of the density is compared to monte carlo simulations in fig .
[ fig : varmaoneonetheoryplusmc ] ; they remain in excellent concord .
these are the main findings of this article .
in this paper we attempted to advertise the power and flexibility of the free random variables calculus for multivariate stochastic processes of the varma type .
the frv calculus is ideally suited for multidimensional time series problems , provided the dimensions of the underlying matrices are large .
the operational procedures are simple , algebraic and transparent .
the structure of the final formula which relates the moments generating function of the population covariance and the sample covariance allows one to easily derive eigenvalue density of the sample covariance .
we in detail illustrated how this procedure works for @xmath6 , confronted the theoretical prediction with numerical data obtained by monte carlo simulations of the varma process and observed a perfect agreement .
the frv calculus is not restricted to gaussian variables .
it also works for non
gaussian processes , including those with heavy  tailed increments belonging to the lvy basin of attraction , where the moments do not exist . since the majority of data collected nowadays is naturally stored in the form of huge matrices , we believe that the frv technique is the most natural candidate for the matrix  valued `` probability calculus '' that can provide efficient algorithms for cleaning ( de  noising ) large sets of data and unraveling essential but hidden correlations .    this work has been supported by the polish ministry of science grant no .
n  n202  229137 ( 20092012 ) .
aj acknowledges the support of clico ltd .
in this appendix , we sketch a proof of the formula ( [ eq : vmaqautocovariancematrix ] ) for the auto  covariance matrix of the @xmath0 process .
as mentioned , since the random variables are centered gaussian , this matrix alone suffices to completely capture all their properties .
we set @xmath116 ; the dependence on this index may be dropped as there are no correlations here .
we use the definition ( [ eq : vmaqdefinition ] ) of @xmath0 , as well as the auto  covariance structure of the white noise , .
this leads to @xmath117 the double sum is symmetrized , the index @xmath118 replaced by @xmath119 , @xmath120 and the order of the sums interchanged ( an elegant method for this is explained in  @xcite ) , @xmath121 which , upon splitting the sum over @xmath122 into three pieces ( from @xmath123 to @xmath124 , @xmath125 , and from @xmath126 to @xmath59 ) , is quickly seen to coincide with ( [ eq : vmaqautocovariancematrix ] ) .        * by we denote the auto  covariance matrix of the @xmath0 process with the generic constants , with @xmath127 , as defined in ( [ eq : vmaqdefinition ] ) . * by we denote the auto  covariance matrix of the @xmath0 process with the constants , , where @xmath78 . * by we denote the auto  covariance matrix of the @xmath1 process with the generic constants , , with @xmath78 , as defined in ( [ eq : varqdefinition ] )
. there holds ( [ eq : varqfromvmaq ] ) . * by we denote the auto  covariance matrix of the @xmath128 process with the constants , , where . * by we denote the auto  covariance matrix of the @xmath2 process with the generic constants , , with and , according to the definition ( [ eq : varmaq1q2definition ] ) .
there is ( [ eq : varmaq1q2fromvarq1vmaq2 ] ) , where in the latter piece .
[ [ ss : themtransformoftheautocovariancematrixforvmaq ] ] the @xmath35transform of the auto  covariance matrix for @xmath0 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    we will derive here the @xmath35transform ( [ eq : momentsgeneratingfunctionforatranslationallyinvariant ] ) of the auto  covariance matrix of an arbitrary @xmath0 process , using the expression for its fourier transform ( [ eq : vmaqautocovariancematrixfourier ] ) .
it is a little simpler to consider the green s function , @xmath129 where the integration range has been halved due to the evenness of the integrand .
this integral is performed with help of the change of variables @xmath130 .
the measure , when @xmath131 $ ] , reads .
a basic observation is that the denominator of the integrand is a linear combination of @xmath132 , for @xmath133 , and each such a cosine can be cast as a polynomial of order @xmath122 in @xmath134 through the de moivre formula .
hence , the denominator is a polynomial of order @xmath59 in @xmath134 , @xmath135 where the s are the @xmath59 roots ( which we assume to be single ) , and @xmath136 is the coefficient at . using the method of residues ,
one readily finds @xmath137 where the two square roots on the r.h.s .
are principal .
this is an explicit formula for the green s function of , provided one has factorized the order@xmath59 polynomial ( [ eq : vmaqgreenfunctionoftheautocovariancematrixdenominator ] ) .
let us argue now that the fourier transform ( [ eq : varqautocovariancematrixfourier ] ) leads to the auto  covariance matrix of @xmath1 ( [ eq : varqautocovariancematrixfromitsfouriertransform ] ) of the form of a sum of exponential decays ( [ eq : varqautocovariancematrix1 ] ) , and let us give precise expressions for the constants and the characteristic times , @xmath86 .
we proceed by the technique of residues , analogously to appendix  [ ss : themtransformoftheautocovariancematrixforvmaq ] , however this time with aid of another variable , , related to the previously used through @xmath138 .
the integration measure is @xmath139 , and the integration path is counterclockwise around the centered unit circle .
the denominator of the integrand is a polynomial of order @xmath59 in the variable @xmath134 , having thus some @xmath59 roots , @xmath78 .
therefore , there are @xmath140 corresponding solutions for the variable @xmath141 , with a half of them inside the integration path and a half outside ; let be the solutions to with the absolute values less than @xmath126 .
only them contribute to the integral , and their residues straightforwardly give @xmath142 this is indeed @xmath59 exponents , @xmath86 .
remark that the solutions may be complex , hence this is really @xmath59 different exponential decays , with the characteristic times @xmath143 ( these times are positive as the roots have the absolute values less than @xmath126 ) , possibly modulated by sinusoidal oscillations when a root has an imaginary part .    for example , for @xmath67
there is one exponential decay ( [ eq : varoneautocovariancematrix ] ) , while for @xmath75 , one obtains either two exponential decays ( the two roots are real and different ) , or one exponential decay modulated by oscillations ( the two roots are complex and mutually conjugate ) , _ etc . _
[ [ ss : theequationforthemtransformofthepearsonestimatorofthecovariancesforvarmaoneone ] ] the equation for the @xmath35transform of the pearson estimator of the covariances for @xmath6 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    the sixth  order polynomial equation obeyed by in the case of @xmath6 reads , @xmath144 @xmath145 @xmath146 @xmath147 @xmath148 @xmath149 @xmath150 @xmath151 @xmath152 @xmath153 @xmath154 @xmath155 @xmath156 @xmath157 @xmath158 @xmath159 @xmath160 @xmath161 @xmath162 @xmath163 this equation in a ` mathematica ` file can be obtained from the authors upon request .
99 sims c. a. , _ macroeconomics and reality _ , econometrica * 48 * ( 1980 ) 1 .
smets f. , wouters r. , _ an estimated stochastic dynamic general equilibrium model of the euro area _ , european central bank working paper series , working paper * 171 * , august 2002 [ ` http://www.ecb.int ` ] .
box g. , jenkins g. m. , reinsel g. , _ time series analysis : forecasting and control _
, prentice  hall , englewood cliffs , nj , 1994 .
laloux l. , cizeau p. , bouchaud j.p . , potters m. , _ noise dressing of financial correlation matrices _ , phys .
* 83 * ( 1999 ) 1467 [ ` arxiv : cond - mat/9810255 ` ] .
plerou v. , gopikrishnan p. , rosenow b. , amaral l. a. n. , stanley h. e. , _ universal and non  universal properties of cross  correlations in financial time series _ , phys .
* 83 * ( 1999 ) 1471 [ ` arxiv : cond - mat/9902283 ` ] .
jin b. , wang c. , miao b. , huang m.n .
l. , _ limiting spectral distribution of large  dimensional sample covariance matrices generated by varma _
, journal of multivariate analysis * 100 * ( 2009 ) 2112 .
bai z. d. , silverstein j. w. , _ spectral analysis of large dimensional random matrices _ , science press , beijing , 2006 .
voiculescu d. v. , dykema k. j. , nica a. , _ free random variables _ , crm monograph series , vol .
1 , am . math .
providence , 1992 .
speicher r. , _ multiplicative functions on the lattice of non  crossing partitions and free convolution _ , math .
* 298 * ( 1994 ) 611 .
burda z. , jarosz a. , jurkiewicz j. , nowak m. a. , papp g. , zahed i. , _ applying free random variables to random matrix analysis of financial data _ , submitted to quantitative finance .
wishart j. , _ the generalized product moment distribution in samples from a normal multivariate population _
, biometrika * a 20 * ( 1928 ) 32 .
zee a. , _ law of addition in random matrix theory _ ,
nucl . phys .
* b 474 * ( 1996 ) 726 [ ` arxiv : cond - mat/9602146 ` ] .
janik r. a. , nowak m. a. , papp g. , zahed i. , _ various shades of blue s functions _ , acta phys
* b 28 * ( 1997 ) 2949 [ ` arxiv : hep - th/9710103 ` ] .
burda z. , grlich a. , jarosz a. , jurkiewicz j. , _ signal and noise in correlation matrix _
, physica * a 343 * ( 2004 ) 295 [ ` arxiv : cond - mat/0305627 ` ] .
burda z. , jurkiewicz j. , _ signal and noise in financial correlation matrices _ , physica * a 344 * ( 2004 ) 67 [ ` arxiv : cond - mat/0312496 ` ] .
burda z. , jurkiewicz j. , wacaw b. , _ spectral moments of correlated wishart matrices _ , phys .
e 71 * ( 2005 ) 026111 [ ` arxiv : cond - mat/0405263 ` ] .
burda z. , grlich a. , jurkiewicz j. , wacaw b. , _ correlated wishart matrices and critical horizons _ , eur .
j. * b 49 * ( 2006 ) 319 [ ` arxiv : cond - mat/0508341 ` ] .
burda z. , jurkiewicz j. , wacaw b. , _ eigenvalue density of empirical covariance matrix for correlated samples _
, acta phys . pol .
* b 36 * ( 2005 ) 2641 [ ` arxiv : cond - mat/0508451 ` ] .
ltkepohl h. , _ new introduction to multiple time series analysis _ , springer verlag , berlin , 2005 .
graham r. , knuth d. , patashnik o. , _ concrete mathematics : a foundation for computer science _
, addison
wesley , 1994 .


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 . ) .