in - network function computation is one of the fundamental paradigms that increases the efficiency of sensor networks vs - a - vs conventional data networks .
sensor nodes , in addition to sensing and communication capabilities , are often equipped with basic computational capabilities . depending on the task for which they are deployed
, a sensor network can be viewed as a distributed platform for collecting and computing a specific function of the sensor data .
for example , a sensor network for environment monitoring may only be concerned with keeping track of the average temperature and humidity in a region .
similarly ` alarm ' networks , such as those for detecting forest fires , require only the maximum temperature .
the baseline approach for performing such tasks is to aggregate all the data at a central node and then perform offline computations ; the premise of in - network computation is that distributed computation schemes can provide sizable improvement in the performance of the network .
however , from the perspective of designing network algorithms , in - network function computation poses a greater challenge than data networks as the freedom to combine and compress packets , as long as the desired information is preserved , destroys the flow conservation laws central to data networks .
the network has a lot more flexibility , so much so as to make quantifying its performance much more challenging @xcite .
our focus in this paper is to develop a queue - based framework for such systems , and use it to design and analyze network algorithms . by network algorithms ,
we refer to cross layer algorithms that jointly perform the following tasks :    1 .   _ aggregating _ the data at nodes via in - network computation , 2 .
_ routing _ packets between nodes , and 3 .
_ scheduling _ links between nodes for packet transmission .
cross - layer algorithms for data networks , although very successful in both theory and increasingly in real - system implementation , are concerned only with the scheduling and routing aspects .
hence , there is a need for a new framework and new algorithms for in - network function computation in sensor networks . keeping in mind the lessons learnt from the success of data networks ,
our aim is to design network algorithms that are _ dynamic _
( i.e. , the algorithm should not be designed assuming static network parameters , but rather , use the network state to adaptively learn the network parameters ) , _ robust _
( i.e. , the algorithm adapts to temporal changes in traffic and network topology ) , _ capable of dealing with a large class of functions _
( i.e. , if the function being computed by the network changes , then one should only need to make minor changes to the scheduling and routing algorithms ) , and _
generalizable to all network topologies_.    due to the wide range of potential applications , there are many existing models for such networks , and many different perspectives from which they are analyzed .
some representative works in this regard are as follows :    * the pioneering work of giridhar and kumar @xcite considers the function computation problem from the point of view of the capacity scaling framework of gupta and kumar @xcite .
in particular , they quantify scaling bounds for certain classes of functions ( symmetric , divisible , type - sensitive , and type - threshold ) under the protocol model of wireless communications and for collocated graphs and random geometric graphs . *
other papers consider the function computation problem from the point of view of information theory @xcite and communication complexity @xcite , characterizing various metrics ( refresh rate , number of messages , etc . ) for different functions in terms of certain properties of the graph , the function to be computed , and the underlying data correlation .
all the above works take an idealized ` bottom - up ' approach to determine the fundamental limits of the problem , and hence are not directly suitable for designing practical network algorithms . *
similar in spirit to the above papers , another approach is to study function computation from the perspective of source coding@xcite .
these works characterize bounds and show the existence of coding schemes for noiseless , wireline networks . as with the previous algorithms , these policies tend to be idealized , using more complex coding - based schemes instead of simple routing and aggregation ( we will later show that such simple strategies suffice for optimal in - network computation of a number of functions of interest ) ; further , these papers do not have explicitly defined policies but rather existence results for such policies .
* in contrast to the ` bottom - up ' approach of all the above works , krishnamachari et al .
@xcite adopt a more ` top - down ' approach whereby they formulate network models that abstract out some of the complexity while allowing quantification of performance gains ( in their case , energy and delay ) .
their models do not , however , achieve the optimal throughput and also do not allow for the design of dynamic network algorithms .
* an alternate model of sensor networks is to assume that nodes are capable of in - network _ compression _ , wherein only the compression ( and not merging ) of flows is permitted .
for example , baek et al .
@xcite consider routing algorithms for power savings in hierarchical sensor networks .
similarly , sharma et al .
@xcite design energy - efficient queue - based algorithms under the assumption that the only operation allowed in addition to routing and scheduling is compression of packets at the source node .
the queue - based model for data networks has proved to be an essential tool in designing provably - efficient algorithms for such systems .
this model has provided a common framework for understanding various aspects of data network performance such as throughput @xcite , delay @xcite , flow utility maximization @xcite , network utility maximization @xcite , distributed algorithms @xcite , among others ( for an overview , refer to @xcite ) . in addition , these algorithms have been implemented in real systems @xcite , including in sensor networks @xcite , with good results . however ,
these algorithms are designed for data networks , and can not exploit any potential benefit from in - network computation .
more recently , this framework has been extended to fork - and - join networks with fixed routing@xcite , and resource allocation in processing networks@xcite .
using fixed routing in a network usually leads to suboptimal operations as the routes may not be designed to optimize the network performance ; in general , even choosing the single best fixed route can perform arbitrarily worse than with dynamic routing ( see example in section  [ sec : refrate ] ) .
further , static routing is not robust to temporal changes in the network .
however , introducing dynamic routing with in - network computation destroys the flow conservation equations that exist in data networks and networks with fixed flows , as the flow out of a node depends both on inflow as well as ( dynamic ) packet aggregation at that node .
thus , there is a need to come up with a new queue - based framework and algorithms for efficient function computation in sensor networks , and our paper is a step in this direction .
our main contributions in this paper are as follows :    * we identify a class of functions , the fully - multiplexible or fmux functions , for which we provide a tight characterization of the maximum refresh rate with in - network computation , i.e. , the maximum rate at which the sensors can generate data such that the computation can be performed in a stable manner .
more formally , we show that for these functions , if the refresh rate exceeds a certain graph parameter ( the stochastic min - min - cut , which we define formally in section [ sec : refrate ] ) , then the system is transient under _ any _ algorithm , whereas for any rate lower than this parameter , we construct a policy that can stabilize the system . * leveraging the results of massouli et al .
@xcite on _ broadcasting _ , we obtain a wireline routing algorithm for _ aggregation _ via in - network computation of fmux functions in directed acyclic graphs .
our approach is based on the observation that broadcasting and aggregation are in some sense , duals , of each other .
more technically , the duality occurs between isolation of packets in aggregation ( i.e. , a packet does not have neighboring packets to aggregate with ) and that of multiple receptions of the same packet ( from different neighbors ) in broadcasting . by suitably modifying the approach in @xcite
, we are able to develop an in - network aggregation algorithm for which routing is completely decentralized , and simple random packet forwarding and aggregation suffices for throughput - optimality . * for general wireless networks
we develop dynamic algorithms based on a centralized allocation of routes ( dynamic flow splitting ) and maxweight - type scheduling .
in particular , we show that loading rounds on trees in a greedy manner ( whereby an incoming round is loaded on the least weighted aggregation tree ) , coupled with an appropriate scheduling rule , is throughput - optimal for computing fmux functions .
the analysis of this algorithm is unique in that in addition to an appropriate lyapunov function , it requires the construction of appropriate tree - packings of the network graph in order to show the throughput - optimality of this routing scheme .
notation : throughout the paper , we use calligraphic fonts ( @xmath1 , etc . ) to denote sets and the corresponding capital letter ( @xmath2 , etc . ) to denote their cardinality .
we interchangeably use @xmath3 or @xmath4 for adding elements to sets , and @xmath5 for deleting elements from sets , and sometimes for brevity of exposition , use the element @xmath6 to denote the singleton set @xmath7 when the meaning is clear from the context ( in particular , for a set @xmath8 and element @xmath6 , @xmath9 ) .
we also use the shorthand notation @xmath10\triangleq\{1,2,\ldots , n\}$ ] .
in this section we describe the system model we study in the rest of the paper . at a high level ,
the system consists of a network of @xmath11 nodes , one of which is the data aggregator and the rest are sensors .
sensor nodes are capable of three tasks : sensing the environment , transmitting to and receiving data from other nodes , and performing computations on the data .
the sensors are assumed to _ sense the environment in a synchronous manner _ , and the overall purpose of the system is to compute a specific function of the synchronously generated sensor data and forward it to the aggregator .
further , the function computation is assumed to be done in a repeated manner , and the metric used to quantify the efficacy of an algorithm is the maximum synchronous rate at which the sensors can generate data such that the required function of the data can be forwarded to the aggregator in a stable manner .
this rate is henceforth referred to as the _ maximum refresh rate _ of the network .
before we describe the queueing framework for function computation , we first outline the general communication model that we consider in this work .
this model is the same as that considered for studying data networks @xcite . in the next section , we will outline the modifications we make in order to capture the in - network computation aspect of a sensor network .    *
communication graph : * we model the topology of the sensor network as a directed graph @xmath12 , where @xmath13 is a set of @xmath11 nodes , and @xmath14 is a set of @xmath15 directed links which determine the connectivity between nodes .
there is a special node , @xmath16 , referred to as the _ aggregator _ , and the rest of the nodes in @xmath13 are sensor nodes .
directed link @xmath17 represents that there exists a communication channel from node @xmath6 to node @xmath18 ( in wireline this corresponds to a physical channel , while in wireless it represents the fact that the nodes are in radio range ) .
* transmission model : * following the convention in literature @xcite , we consider a continuous time system in case of wireline systems , whereas in the case of wireless networks , we assume that time is slotted , and state all rates in bits per slot . in wireline networks ,
we define a vector of link rates @xmath19 ; one bit is assumed to traverse a link @xmath20 with a random transit time with distribution _ _ exponential__@xmath21 .
the transit times are independent across links and across packets crossing the same link .    for wireless networks ,
we make the following assumptions / definitions @xcite :    * we assume that the channels between nodes are constant ( but can extend to time varying channels with added notation , see @xcite ) .
the wireless nature of the network is reflected in the interference constraints . * for transmission schedule @xmath22
, @xmath23 denotes the link - rate vector of transmission rates over the links under the chosen schedule .
* @xmath24 is the set of valid schedules that obey the interference constraints ( henceforth referred to as independent sets ) .
@xmath25 is said to be _
admissible _ if the link - rates can be achieved simultaneously in a time slot .
@xmath26 is the set of all admissible @xmath25 and is assumed to be time invariant as stated above .
further , we assume that @xmath27 .
* @xmath28 is said to be _ obtainable _ if @xmath29 , the convex hull of @xmath30 .
an obtainable link - rate vector can be achieved by time sharing over admissible link - rate vectors .
* from the definition of the convex hull , we have that for every obtainable rate vector @xmath29 , there exists a probability measure @xmath31 over @xmath32 such that @xmath33 .
the vectors @xmath34 are called static service split ( or sss ) rules@xcite , and represent time sharing fractions between different independent sets in order to achieve the rate @xmath28 .    up to this point ,
the system is identical to one considered for data networks . to highlight the unique features of a physical sensor network performing function computation ( and how they affect the modeling of such a system ) , we consider the following example . in the process
, we also indicate the gains achievable via in - network computation versus data download and processing at the aggregator .
[ cols= " < " , ]      finally we present the complete dynamic algorithm for fmux computation .
the algorithm separates the routing and scheduling components as follows : when a round of packets arrives in the network , we first ` load ' all packets of the round on an aggregation tree ( thereby fixing the routing ) ; next , in each time slot , scheduling is done according to a modified maxweight policy .
the routing is performed using a _ greedy tree - loading policy _ , wherein all incoming rounds in a time slot are loaded on the tree with smallest sum useful - queue , i.e. , least number of useful packets .
formally , we have :    [ algo : treechoosealgo ] calculate @xmath35\right)$ ] for all @xmath36 .
+ find the minimum loaded tree @xmath37 $ ] as : @xmath38=\arg\min_{\tau\in\mathcal{t } } w_{\tau}[t].\ ] ] assign all incoming rounds to aggregation tree @xmath37 $ ] .
the scheduling algorithm is similar to the maxweight policy@xcite , in that it picks a maximum independent set with weights given by the product of the rate and the maximum queue across an edge .
formally we have the following algorithm :    [ algo : compbpalgo1 ] place packets arriving on tree @xmath39 at node @xmath40 in @xmath41 $ ] for non - leaf nodes , and @xmath42 $ ] for leaf nodes .
+ calculate @xmath43=\max_{\tau\in\mathcal{t}:(i , j)\in\tau}q_i^{\tau , u}[t]$ ] .
also define @xmath44 $ ] as the tree which maximizes @xmath43 $ ] .
+ compute schedule @xmath45 $ ] as : @xmath46=\arg\max_{\mathbf{c}\in\gamma}\sum_{(i , j)\in\mathcal{l}}p_{ij}[t]c_{ij}[t].\ ] ] + consider link @xmath47 .
if @xmath48>0 $ ] , then transmit the first @xmath49 , q_i^{\tau^*(i , j)[t],u}[t])$ ] packets from @xmath50,nu}[t]$ ] .
for the sake of completeness , we note that in all the above algorithms , tie - breaking rules as well as the service discipline ( i.e. , among a set of multiple packets suitable for transmission , which one gets priority ) are assumed to be random ; this is done for the sake of convenience , and we note that there are many possible tie - breaking rules and service disciplines which would suffice .
now we can state and prove the throughput optimality of this algorithm .
[ thm : treesss ] the dynamic queue based policy consisting of greedy tree loading ( algorithm [ algo : treechoosealgo ] ) and maxweight scheduling ( algorithm [ algo : compbpalgo1 ] ) stabilizes the system for any refresh rate @xmath51 that is less than the maximum refresh rate @xmath52 .    before proceeding further , we point out a particular novel aspect of the proof of this theorem .
similar to previous papers @xcite , we use a quadratic lyapunov function for showing stability ; however our technique for bounding the lyapunov drift is quite different from those used for point - to - point data .
the difficulty arises from the fact that although edmonds theorem guarantees the existence of an optimal tree - packing for the network , the trees in this optimal packing are unknown to the algorithm ; consequently it is unclear whether routing over all trees could lead to instability via packet accumulation on trees not involved in the optimal packing .
we circumvent this by showing the existence of some intermediate tree packings between the optimal and the desired refresh rates , which allow uniform bounding of the lyapunov drift .
we now present the complete proof .
we define a candidate lyapunov function @xmath53 $ ] as @xmath54=\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]^2,\end{aligned}\ ] ] with corresponding lyapunov drift given by @xmath55=\mathbb{e}[v[t+1]-v[t]|\mathbf{q}[t]].\ ] ] similar to before , we have that @xmath53\geq 0 $ ] for all states of the system , and that @xmath56<\infty$ ] . we now need to show that given @xmath57 , there exists @xmath58 such that if @xmath59>q_{\max}$ ] for some @xmath40 , then @xmath56<-\delta$ ] .
now we have @xmath55=\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}\mathbb{e}\left[(\delta q_i^{\tau , u}[t])^2 + 2q_i^{\tau , u}[t]\delta q_i^{\tau , u}[t]|\mathbf{q}[t]\right],\ ] ] and defining @xmath60 $ ] to be arrival of useful packets on tree @xmath39 to node @xmath40 , we have @xmath61=a_i^{\tau , u}[t]-d_{(i , p^{\tau}(i))}^{\tau}[t],\ ] ] and thus @xmath62)^2\leq m_a + ( lc_{\max})^2 $ ] ( due to external arrivals plus inter - node transmissions ) .
let @xmath63 .
then we have @xmath55\leq m_2 + 2\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]\mathbb{e}\left[a_i^{\tau , u}[t]-d_{(i , p^{\tau}(i))}^{\tau}[t]|\mathbf{q}[t]\right].\ ] ]    from the definition of @xmath52 , we know that there exists an optimal rate point @xmath64 and the corresponding optimal sss rule @xmath65 that maximizes the min - mincut .
consider now a refresh rate @xmath51 less than the @xmath52 , such that @xmath66 .
note that the algorithm can potentially split the incoming flow @xmath51 over _ every spanning tree _ of the network , in order to dynamically arrive at the optimal packing . to uniformly bound the lyapunov drift
, we first need to construct two tree packings : an ` achievable ' packing @xmath67 such that @xmath68 which serves as a proxy for the flow - splitting , and a ` near - optimal ' packing @xmath69 such that @xmath70 and further which has the property that @xmath71 _ uniformly over all spanning trees _ ( for some @xmath72 which we define below ) .
we do so as follows .
assume that there exists @xmath73 such that if any edge @xmath20 is scheduled alone ( i.e. @xmath74 ) , then @xmath75 ( this is simply a formal definition of existence of a link ) .
we can now perturb the optimal sss rule to get a new rate point @xmath76 with the following two properties :    1 .
every edge @xmath20 has capacity @xmath77 .
2 .   the min - mincut of the network at the rate point @xmath78 is @xmath79 .
this helps ensure that the ` near - optimal ' tree packing @xmath80 can have some mass on each edge of the graph .    to construct the perturbed sss rule @xmath81 ,
consider the optimal sss rule @xmath82 .
we define @xmath83 ( i.e. , the set of independent sets that have some mass under @xmath65 ) and @xmath84 ( which is @xmath85 as the cardinality of @xmath86 is finite ) .
now we reduce each @xmath87 by @xmath88 .
this reduces the min - mincut by at most @xmath89 . to see this , note that the capacity of any edge @xmath90 reduces from @xmath91 to @xmath92 where : @xmath93 further , the maximum number of edges across the min - mincut is bounded by @xmath15 .
thus the min - mincut of the network at the rate point @xmath78 is @xmath79 .
next , suppose @xmath94 is the set of edges with zero flow under @xmath65 .
we now complete the definition of the perturbed sss rule @xmath81 ( using the fact that singleton edges are valid independent sets ) as follows : @xmath95 to see that this is a valid sss rule , note that @xmath96 , which is the weight we have distributed equally over all links in @xmath94 .
the rate point under this sss rule is henceforth denoted as @xmath97 .
then for edges in @xmath94 we have @xmath98 .
now since there are only @xmath15 edges , each with positive capacity @xmath92 , therefore there exists some @xmath99 such that every edge @xmath20 has capacity @xmath100 under sss rule @xmath81 .
finally , applying edmonds theorem ( theorem  [ thm : edmonds ] ) on the network under @xmath81 , we get a packing @xmath101 such that we have @xmath102 before proceeding further , we need the following definitions :    * @xmath103 .
* @xmath104 .
* @xmath105 . * @xmath106 .
note that @xmath107 as @xmath108 and the packing is not tight on the finite set @xmath109 .
similarly , @xmath110 .
finally we can construct the tree packings ( on the network under sss rule @xmath81 ) that we need to bound the lyapunov drift :    1 .
the ` achievable ' tree packing , @xmath111 is defined as : @xmath112 then clearly @xmath113 is a packing ( as we are only removing mass from a valid packing ) and further : @xmath114 2 .
the ` near - optimal ' tree packing , @xmath115 is defined as : @xmath116 first we need to show that this is a valid tree packing .
to see this , note that the maximum load added on any edge is bounded by @xmath117 ( since in the worst case , all the trees in @xmath118 can contain some edge ) . for any edge in @xmath109 , this is less than the slack ( @xmath119 by definition ) that was already present .
for an edge in @xmath120 , we know at least one tree in @xmath121 contained it ( as every edge in the graph has positive capacity under the sss rule @xmath81 ) , and hence we subtract a load of at least @xmath122 , which is again greater than the amount of load we add .
thus @xmath115 is a valid packing .
+ further we have that @xmath123 .
+    in addition , defining @xmath124 we get that @xmath125 .
thus we have constructed the two tree packings we need .
we now return to bounding the lyapunov drift . from above
, we have @xmath55\leq m_2 + 2\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]\mathbb{e}\left[a_i^{\tau , u}[t]-d_{(i , p^{\tau}(i))}^{\tau}[t]|\mathbf{q}[t]\right].\ ] ] now , let @xmath126 $ ] be the rate for packets on aggregation tree @xmath39 on link @xmath127 allocated by the policy in time slot @xmath128 ( thus @xmath129=c_{(i , p^{\tau}(i))}[t]$ ] ) .
then we have @xmath130\mathbb{e}&\left[d_{(i , p^{\tau}(i))}^{\tau}[t]|\mathbf{q}[t]\right]\\ & = \sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]\mathbb{e}\left[c_{ip^{\tau}(i)}^{\tau}[t]-\max\{c_{ip^{\tau}(i)}^{\tau}[t]-q^{\tau , u}_i[t],0\}|\mathbf{q}[t]\right]\\ & \geq\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]\mathbb{e}\left[c_{ip^{\tau}(i)}^{\tau}[t]-\max\{c_{\max}-q^{\tau , u}_i[t],0\}|\mathbf{q}[t]\right]\\ & \geq\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]\mathbb{e}\left[c_{(i , p^{\tau}(i))}^{\tau}[t]\right]-nl|\mathcal{t}|c_{\max}^2.\end{aligned}\ ] ] further , from the definition of the policy , we know that @xmath130\mathbb{e}&\left[c_{(i , p^{\tau}(i))}^{\tau}[t]|\mathbf{q}[t]\right]\\ & = \mathbb{e}\left[\max_{\mathbf{c}\in\gamma}\sum_{(i , j)\in\mathcal{l}}\max_{\tau\in\mathcal{t}:(i , j)\in\tau}\left\{q_i^{\tau , u}[t]\right\}c_{ij}[t]|\mathbf{q}[t]\right]\\ & \geq\max_{\mathbf{c}\in\gamma}\sum_{(i , j)\in\mathcal{l}}\max_{\tau\in\mathcal{t}:(i , j)\in\tau}\left\{q_i^{\tau , u}[t]\right\}\mathbb{e}\left[c_{ij}[t]|\mathbf{q}[t]\right]\\ & \geq\max_{\mathbf{c}\in\mathcal{ch}(\gamma)}\sum_{(i , j)\in\mathcal{l}}\max_{\tau\in\mathcal{t}:(i , j)\in\tau}\left\{q_i^{\tau , u}[t]\right\}\mathbb{e}\left[c_{ij}[t]|\mathbf{q}[t]\right]\\ & \geq\max_{\mathbf{c}\in\mathcal{ch}(\gamma)}\sum_{(i , j)\in\mathcal{l}}\sum_{\tau\in\mathcal{t}:(i , j)\in\tau}q_i^{\tau , u}[t]\mathbb{e}\left[c_{ij}^{\tau}[t]|\mathbf{q}[t]\right]\\ & \mbox{(where $ c_{ij}^{\tau}[t]$ is any tree - packing of a given $ \mathbf{c}\in\mathcal{ch}(\gamma)$)}\\ & \geq\sum_{(i , j)\in\mathcal{l}}\sum_{\tau\in\mathcal{t}:(i , j)\in\tau}q_i^{\tau , u}[t]\hat{c}_{ij}^{\tau},\end{aligned}\ ] ] where for any edge @xmath47 , @xmath131 represents any valid split of @xmath132 between trees lying on that edge , i.e. , @xmath133 .
in particular , therefore , we can use the tree packing @xmath115 to get @xmath134\mathbb{e}\left[c_{(i , p^{\tau}(i))}^{\tau}[t]|\mathbf{q}[t]\right ] \geq\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]\hat{\lambda}^{\tau}.\end{aligned}\ ] ] combining inequalities , and defining @xmath135 we get @xmath136\leq m_3 + 2\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t]\left(\mathbb{e}\left[a_i^{\tau}[t]|\mathbf{q}[t]\right]-\hat{\lambda}^{\tau}\right).\end{aligned}\ ] ]    finally , define @xmath137 $ ] to be the rate of rounds arriving on tree @xmath39 . then from the greedy round - tree assignment algorithm , and using the fact that each round results in exactly one useful packet at each node , we get    _ i_q_i^,u[t ] & + & = _ i_q_i^,u[t ] + & = + & _ \{a^[t]}_i_q_i^,u[t ] + & _
i_q_i^,u[t]_.    thus we get @xmath136&\leq m_3 - 2\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t](\hat{\lambda}_{\tau}-\lambda_{\tau}')\\ & \leq m_3 - 2\epsilon_4\sum_{i\in\mathcal{n}}\sum_{\tau\in\mathcal{t}}q_i^{\tau , u}[t].\end{aligned}\ ] ] in order to have @xmath56<-\delta$ ] if @xmath138\geq q_i^{\tau , u}[t]>q_{\max}$ ] for some @xmath139 , we can choose @xmath140 .
thus @xmath53 $ ] is a valid lyapunov function and by foster s theorem , our policy is stabilizing for any @xmath141 .
we have presented a queueing - based framework for in - network function computation .
we have used this framework to gain insights into designing dynamic and distributed algorithms for in - network function computation and to quantify the performance gains over data - download .
we have focused on a class of functions , the fmux functions , which exhibit maximum compression on aggregation , and for which we have used the parity and max functions as representative examples . for such functions we have developed scheduling and routing algorithms under different settings . for wireline networks ,
we have extended the random routing scheme of massouli et al .
@xcite for aggregation . for wireless networks , we have provided a fixed - routing via dynamic flow splitting along with maxweight - like scheduling , which is shown to be throughput - optimal .
the wireless algorithm , as presented , requires routing on all aggregation trees in order to achieve throughput optimality ; this may not be practical in many networks due to the potentially exponential number of trees .
however , as we showed in the example with the complete graph , one can obtain optimal tree packings with a much smaller number of trees ( of the order of @xmath15 ) and one direction of future work is to show how such trees can be selected using simple rules in different networks .    generalizing these algorithms to deal with a broader class of functions , as well as studying the performance of the algorithms with respect to other metrics ( delay , energy consumption , among others ) are other topics for future work .
this work was supported in part by afosr under grant fa9550 - 09 - 1 - 0317 , nsf grants cns-0519535/0519401 ( collaborative grant ) and cns-0964391 .
10    a.  giridhar and p.  r. kumar , `` computing and communicating functions over sensor networks , '' _ ieee journal on selected areas in communications _ ,
23 , no .  4 , pp .
755764 , 2005 .
p.  gupta and p.  r. kumar , `` the capacity of wireless networks , '' _ ieee transactions on information theory _ , vol .
46 , no .  2 , pp .
388404 , 2000 .
a.  orlitsky and j.  r. roche , `` coding for computing , '' _ ieee transactions on information theory _ , vol .
47 , no .  3 , pp .
903917 , 2001 .
n.  ma , p.  ishwar , and p.  gupta , `` information - theoretic bounds for multiround function computation in collocated networks , '' _ corr _ , vol .
abs/0901.2356 , 2009 .
h.  kowshik and p.  r. kumar , `` optimal computation of symmetric boolean functions in tree networks , '' in _ ieee international symposium on information theory  isit 2010 _ , ieee , july 2010 .
r.  appuswamy , m.  franceschetti , n.  karamchandani , and k.  zeger , `` network coding for computing part i : cut - set bounds , '' _ corr _ , vol .
abs/0912.2820 , 2009 .
t.  moscibroda , `` the worst - case capacity of wireless sensor networks , '' in _ internat .
workshop on inform .
processing in sensor networks ( ipsn ) _ , pp .  110 , 2007 .
b.  krishnamachari , d.  estrin , and s.  b. wicker , `` the impact of data aggregation in wireless sensor networks , '' in _ proc .
the 26th ieee internat .
conf . distributed computing systems ( icdcs ) _ , pp .  575578 , 2002 .
n.  karamchandani , l.  keller , c.  fragouli , and m.  franceschetti , `` function computation via subspace coding , '' in _ proc .
information theory ( isit ) _ , ieee , july 2010 .
s.  j. baek , g.  de  veciana , and x.  su , `` minimizing energy consumption in large - scale sensor networks through distributed data compression and hierarchical aggregation , '' _ ieee journal on selected areas in communications _ ,
22 , no .  6 , pp .
11301140 , 2004 .
a.  sharma , l.  golubchik , r.  govindan , and m.  j. neely , `` dynamic data compression in multi - hop wireless networks , '' in _ proc .
acm sigmetrics / performance conf .
_ , pp .  145156 , 2009 .
l.  tassiulas and a.  ephremides , `` stability properties of constrained queueing systems and scheduling policies for maximum throughput in multihop radio networks , '' _ ieee transactions on automatic control _ , vol .  4 , pp .
19361948 , december 1992 .
m.  andrews , k.  kumaran , k.  ramanan , a.  l. stolyar , r.  vijayakumar , and p.  whiting , `` scheduling in a queueing system with asynchronously varying service rates , '' _ probability in engineering and informational sciences _ , vol .  14 , pp .  191217 , 2004 .
l.  bui , r.  srikant , and a.  l. stolyar , `` novel architectures and algorithms for delay reduction in back - pressure scheduling and routing , '' in _ proc .
ieee infocom .
_ , pp .  29362940 , 2009 .
m.  j. neely , `` delay analysis for maximal scheduling with flow control in wireless networks with bursty traffic , '' _ ieee / acm trans .
_ , vol .
17 , no .  4 , pp .
11461159 , 2009 .
x.  lin , n.  shroff , and r.  srikant , `` a tutorial on cross - layer optimization in wireless networks , '' _ ieee journal on selected areas in communications _ , pp .  14521463 , august 2006 .
l.  bui , r.  srikant , and a.  l. stolyar , `` optimal resource allocation for multicast flows in multihop wireless networks , '' _ proc .
acm sigmetrics / performance evaluation review conf .
_ , vol .
35 , no .  3 , p.  43
, 2007 .
l.  lin , n.  b. shroff , and r.  srikant , `` energy - aware routing in sensor networks : a large system approach , '' _ ad hoc netw .
_ , vol .  5 , no .  6 , pp .
818831 , 2007 .
l.  ying , s.  shakkottai , and a.  reddy , `` on combining shortest - path and back - pressure routing over multihop wireless networks , '' in _ proc .
ieee infocom .
_ , pp .  16741682 , 2009 .
a.  l. stolyar , `` maximizing queueing network utility subject to stability : greedy primal - dual algorithm , '' _ queueing syst .
_ , vol .
50 , no .  4 , pp .
401457 , 2005 .
l.  jiang , d.  shah , j.  shin , and j.  c. walrand , `` distributed random access algorithm : scheduling and congesion control , '' _ corr _ , vol .
abs/0907.1266 , 2009 .
l.  georgiadis , m.  j. neely , and l.  tassiulas , `` resource allocation and cross - layer control in wireless networks , '' _ foundations and trends in networking _ , vol .  1 , pp .  1144 , 2006 .
u.  akyol , m.  andrews , p.  gupta , j.  d. hobby , i.  saniee , and a.  l. stolyar , `` joint scheduling and congestion control in mobile ad - hoc networks , '' in _ proc .
ieee infocom .
_ , pp .  619627 , 2008 .
j.  ryu , v.  bhargava , n.  paine , and s.  shakkottai , `` design , implementation and evaluation of back - pressure routing / rate control for intermittently connected networks , '' in _ acm int .
conf . on mobile computing and networking ( mobicom ) _ ,
s.  moeller , a.  sridharan , b.  krishnamachari , and o.  gnawali , `` routing without routes : the backpressure collection protocol , '' in _ internat .
workshop on inform .
processing in sensor networks ( ipsn ) _ , pp .  279290 , 2010 .
h.  zhao , c.  h. xia , z.  liu , and d.  towsley , `` a unified modeling framework for distributed resource allocation of general fork and join processing networks , '' in _ proc .
acm sigmetrics conf .
_ , 2010 .
l.  jiang , and j.  c. walrand , `` stable and utility - maximizing scheduling for stochastic processing networks , '' _ proc .
allerton conf .
communication , control , and computing _
, 2009 .
l.  massouli , a.  twigg , c.  gkantsidis , and p.  rodriguez , `` randomized decentralized broadcasting algorithms , '' in _ proc .
ieee infocom .
_ , pp .  10731081 , 2007 .
h.  kowshik and p.  r. kumar , `` zero - error function computation in sensor networks , '' in _ proc .
conf . on decision and control
_ , pp .  37873792 , 2009 .
s.  madden , m.  j.  franklin , j.  m.  hellerstein , and w.  hong , `` tag : a tiny aggregation service for ad - hoc sensor networks , '' in _ sigops oper .
_ , vol .
36 , pp .  131146 , 2002 .
j.  edmonds , `` edge - disjoint branchings . '' in _ combinatorial algorithms _ , 1972 .
we now present the complete proof for the throughput - optimality of algorithm [ algo : wlinerpfalgo ] in directed acyclic graphs .
since the proof closely follows the proof of massouli et al.@xcite , we do not go into complete details , but try mainly to highlight the modifications we make in order to perform aggregation rather than broadcast .
first we need a lemma that ensures that under the useful packet transmission rule , each round of packets follows a spanning tree .
recall that the footprint of a round of packets is defined as the set of nodes in which the packets of that round is present .
further , recall that a set @xmath8 is said to be a valid footprint set if each node in @xmath8 has a path to @xmath142 in the subgraph induced by @xmath8 ; the collection of such sets is denoted by @xmath143 .
we assume throughout that @xmath144 , for otherwise the min - mincut is @xmath145 .
note that since we operate in continuous time , only one packet transmission ocurs at a given time with probability @xmath146 ; further , we require that the local state information is available at the time of making routing decision .
now we have the following lemma :      since the underlying graph is directed acyclic , we re - label the nodes as @xmath149 according to their topological ordering , where node @xmath145 is the aggregator @xmath142 , and all edges are from a higher numbered node to a lower numbered node .
further , given a round of packets on a valid footprint set @xmath8 , we have that each node @xmath150 has at least one route to @xmath142 using only nodes in @xmath8 ; for short , we refer to such a route as a _ path from @xmath151 to @xmath142 in @xmath8_.    since we are operating in continuous time , with probability @xmath146 only one packet transmission occurs at a given time .
now suppose a useful packet is transmitted on edge @xmath152 , where @xmath153 , resulting in a new footprint set @xmath154 . for @xmath148 to be a valid footprint ,
we need that even after the transmission , each node @xmath155 has a path to @xmath142 in @xmath148 . to do this
, we need to consider a partition of the nodes in @xmath148 into @xmath156 classes :    * node @xmath155 such that @xmath157 in the topological order : due to the topological ordering property , a path from @xmath151 to @xmath142 in @xmath8 is clearly unaffected by the packet transmission from @xmath158 to @xmath40 . *
node @xmath159 such that there exists a path from @xmath151 to @xmath142 in @xmath8 which does not include @xmath158 : such a path is also unaffected by the packet transmission from @xmath158 to @xmath40 and hence is still present in @xmath148 .
* node @xmath159 such that _ all paths _ from @xmath151 to @xmath142 in @xmath8 pass through @xmath158 : we show by contradiction that this case is impossible under the rules of useful packet transmission .
for any path from @xmath151 to @xmath142 in @xmath8 , let @xmath160 be the node immediately before @xmath158 ( i.e. , the path is @xmath161 ) .
then @xmath162 has no path to @xmath142 in @xmath8 that does not pass through @xmath158 , for otherwise we have a path from @xmath151 to @xmath162 , and then to @xmath142 , which does not pass through @xmath158 .
this means that @xmath162 becomes isolated upon transmission of packet from @xmath158 to @xmath40 , which violates the non - isolation condition of useful packet forwarding .      the main idea behind the proof in @xcite was to define the ` footprint counter ' variables to represent the state of the system , and considering an appropriate function of these that allowed translating the local decisions of the nodes in terms of global graph parameters .
in order to modify the proof for broadcast , we defined a similar collection of counter variables in section [ sec : algo1 ] , and now define their associated dynamics as follows .
* arrival of new round : @xmath163 ( this corresponds to adding a packet to the queue with footprint @xmath13 , as a packet of the new round is simultaneously generated at all the nodes ) .
* completion of packet transfer : this is only for active packets , i.e. , those currently under transmission . for active packet @xmath164 with corresponding @xmath165 and @xmath166 , we have : @xmath167 ( the first equation corresponds to removing the edge over which packet transmission was completed , and also updating the footprint of the packet to include the new node .
the second updates the list of idle packets in case there is no other instance of this packet being transmitted . ) .
* initiation of a new transfer at an idle link .
the new packet is selected uniformly at random among the set of useful packets at the node .
if @xmath168 , then a new packet transfer is formally described as follows : * * select a useful packet of an idle round with footprint @xmath169 , with probability @xmath170 select a useful packet of an active round @xmath164 with @xmath171 with probability @xmath172 * * if idle packet with footprint @xmath8 is selected : @xmath173 , with @xmath174 .
if packet of active round @xmath175 is selected , then @xmath176    we note here that the node itself does not need to know these global counters to perform packet selection ; rather , this emerges from the use of the random useful packet forwarding rule .
the idea of relating the local packet selection rule to the global counters is crucial in proving the optimality of the algorithm . the local rules for checking whether a packet is useful or not corresponds to selecting packets whose global footprint obeys certain properties ; picking a useful packet uniformly at random therefore corresponds to picking a packet from such a useful global footprint with a probability proportional to the corresponding counter variable .
observe that in order to determine the flow into a footprint set @xmath8 , we need to consider the collection of sets which include @xmath8 and have one extra node .
we now define the fluid limits of the system .
this is similar in spirit to the fluid limit of the system in @xcite , so we try to use similar notation .
the existence of the limit also follows immediately from their convergence results , so we omit it due to lack of space and refer interested readers to @xcite for technical details .
* @xmath178 @xmath179 * work conservation : at almost every @xmath180 is differentiable and if @xmath181 ( where @xmath182 is the fluid trajectory associated with @xmath183 ) , then we have @xmath184 * @xmath185 are non - decreasing , lipschitz continuous , with lipschitz constant @xmath186 , and @xmath187 is @xmath186-lipschitz .
for any @xmath188 set of all fluid trajectories with initial condition @xmath189 , and further , we define @xmath190 as the state of the mc with initial conditions @xmath191 , @xmath192 .
now , as in @xcite , for a sequence of initial conditions @xmath193 s.t . for a sequence of positive numbers @xmath194 and
the limit @xmath195 exists in @xmath196 , we have that @xmath197 : @xmath198}||y^n(t)-f(t)||\geq\epsilon]=0.\ ] ]      next we define the candidate lyapunov function that we use to analyze the stability of the system . in @xcite ,
the function was defined in terms of queues ( or counters ) that counted all the packets whose footprint was contained inside a set @xmath8 .
the advantage of these queues for studying broadcast was that their rate of increase was controlled by external arrivals to the system , while they were drained due to transfers across the cut defined by the set @xmath8 .    for the purpose of studying aggregation ,
we need to identify an equivalent set of queues to reflect the unique dynamics of the system .
in particular , we consider for each set @xmath8 a queue of all rounds whose footprints _ are not entirely contained within @xmath8_. these queues ( counters ) exhibit similar properties to the ones considered for broadcast in that every incoming round is counted by all these queues ( as every node in the network generates a packet ) , while the drain of these queues is controlled by flow across the cut defined by the set @xmath8 .
formally , we have the following theorem :      as in @xcite , before proving this theorem we first we need a combinatorial lemma .
this lemma and its proof parallels a corresponding lemma in @xcite , with modifications to deal with aggregation and the @xmath206 counter variables we have defined above .        note that lemma [ lem : counting ] does not depend on the algorithm , or the fluid model in any way .
it is a pure combinatorial property of the way that the quantities are defined .
in other words , any function mapping the sets @xmath147 to @xmath217 obeys the lemma _ for any @xmath218_. later we use the ability to control @xmath219 to obtain uniform bounds on the lyapunov drift .    for the first condition ,
consider @xmath169 such that @xmath220 then we have @xmath221 and thus @xmath222 however , from the definition of the @xmath223 , we have that @xmath224 for all @xmath225 .
hence we have that @xmath226    for the second condition , consider @xmath147 such that@xmath227 .
further , consider set @xmath148 such that @xmath228 and satisfying @xmath229 then we have @xmath230 thus for our condition , we need @xmath231 and noting the fact that @xmath223 are increasing with @xmath40 , it is sufficient to ensure @xmath232 this in fact holds with equality because of our choice of @xmath223 . thus , given any @xmath218 , we can construct @xmath223 such that the two conditions hold .
( proof of theorem [ thm : lyap ] ) given @xmath218 , we define @xmath223 as in lemma [ lem : counting ] . then , or any @xmath233 , if @xmath234 is a set which belongs to arg - max of @xmath235 , then @xmath236 ( unless all the fluid sample paths are identically @xmath145 ) .
next we use the optimality of @xmath234 to obtain some relations between @xmath237 and the weight across its cut - edges .
@xmath238 such that @xmath17 , we have from the contrapositive of the first condition of lemma [ lem : counting ] ( as @xmath234 is in the arg - max ) that @xmath239 similarly from condition @xmath240 , @xmath241 such that @xmath242 , we have that @xmath243 now we have @xmath244\\ = & \lambda-\sum_{v\in s^ * , u\notin s^*}c_{uv}\left[1-\sum_{s'\nsubseteq s^ * , v\in s ' , u\notin s'}\frac{x_{s'+u}}{x_{+u - v}}\right],\\ & \mbox{(from defn of fluid trajectories)}\\ \leq & \lambda-\sum_{u\in s^ * , v\notin s^*}c_{uv}+\sum_{u\in s^ * , v\notin s^*}c_{uv}\sum_{s'\nsubseteq s^ * , v\in s ' , u\notin s'}\alpha\\ & \mbox{(from previous observation)}\\ \leq & \lambda-\sum_{v\in s^ * , u\notin s^*}c_{uv}+\max_{(u , v)\in \mathcal{l}}c_{uv}|\mathcal{l}|2^k\alpha.\end{aligned}\ ] ] if we choose @xmath219 and @xmath245 as follows : @xmath246 then we get that , for all @xmath247 , @xmath248 to argue that this implies negative drift of the lyapunov function , i.e. @xmath249 , we observe that by definition @xmath250 . finally , using the lipschitz continuity of the trajectories ,
it is sufficient to show this property holds for the sets @xmath247 .


infrared luminous galaxies ( irlgs ) , most of whose huge luminosities are radiated in the far - infrared ( 40500 @xmath0 m ) as dust emission ( @xmath2 @xmath4 10@xmath5 ergs s@xmath6 ) , @xmath7 75 km s@xmath6 mpc@xmath6 and @xmath8 @xmath7 0.5 are used throughout this paper . ]
are a significant population at the bright end of galaxy luminosity function ( soifer et al .
two types of activity , active galactic nucleus ( agn ) and starburst , are thought to contribute to the far - infrared emission of irlgs .
the estimate of the relative energetic contributions of these kinds of activity is an important issue for understanding the nature of irlgs .
furthermore , since the bulk of the cosmic submillimeter background emission can be explained by irlgs at high - redshift ( blain et al .
1999 ) , understanding the energy sources of nearby irlgs is closely related to the origin of the background emission .
both agn and starburst activity produces strong emission lines . since many emission lines , such as narrow emission lines of hydrogen or h@xmath9 molecules ,
can be caused by both kinds of activity , it is not easy to distinguish their origin . the polycyclic aromatic hydrocarbon ( pah ) emission features can be a powerful tool for separating these kinds of activity and estimating their relative importance because the features are associated only with starburst activity and not with agn activity ( moorwood 1986 ) . using the pah emission feature at 7.7 @xmath0 m , systematic studies of the energy sources of irlgs have been reported based on the _ infrared space observatory _ ( _ iso _ ) spectra at 5.811.6 @xmath0 m ( lutz et al .
1998 ; genzel et al .
1998 ; rigopoulou et al .
the observed ratios between the 7.7 @xmath0 m pah peak flux and 7.7 @xmath0 m continuum flux of irlgs were estimated to be as high as those of starburst galaxies with lower far - infrared luminosities , leading these authors to suggest that starburst activity is a dominant energy source for most of the irlgs they observed . however , the 7.7 @xmath0 m pah to far - infrared luminosity ratios of irlgs are roughly half of those of starburst galaxies ( genzel & cesarsky 2000 ) .
this indicates that observed starburst activity can account for only half of the far - infrared luminosities for irlgs , although higher dust extinction toward starburst activity in irlgs could make the intrinsic contribution of starburst activity higher .    in the 5.811.6 @xmath0 m spectra of irlgs ,
particularly obscured ones , not only the pah emission , but also broad silicate dust absorption , is likely to exist at 813 @xmath0 m ( mathis 1990 ) .
limited wavelength coverage of most _ iso _ spectra ( 5.811.6 @xmath0 m ) does not allow the determination of a secure continuum level for many irlgs and thus the estimate of the pah flux could be quantitatively uncertain , as mentioned in rigopoulou et al . ( 1999 ) .
in fact , the 7.7 @xmath0 m pah emission flux of arp 220 , an irlg for which an _ iso _ spectrum with wider wavelength coverage ( 516 @xmath0 m ) is available , could be considerably smaller than that estimated by rigopoulou et al .
, if the continuum level is determined in a different way ( elbaz et al .
1998 ; dudley 1999 ) . furthermore ,
in the case of another irlg , iras 05189@xmath102524 , the reported equivalent widths of the 7.7 @xmath0 m pah emission feature based on the _ iso _ spectrum differ by a factor of @xmath113 according to different authors ( clavel et al .
2000 ; laureijs et al .
2000 ) .
spectroscopy in the range 34 @xmath0 m can be a powerful method for investigating the energy source of irlgs , by avoiding the above uncertainty .
first , since dust extinction at 34 @xmath0 m is estimated to be as small as that at 78 @xmath0 m
( @xmath12/@xmath13 @xmath3 0.8 ; lutz et al .
1996 ) , we can detect signs of obscured energy sources .
next , we can diagnose the relative importance of agn and starburst activity using spectral features in the range 34 @xmath0 m .
if an irlg is powered by starburst activity , pah emission should be detected at 3.3 @xmath0 m and sometimes at 3.4 @xmath0 m ( tokunaga et al .
if a source is powered by obscured agn activity , a 3.4 @xmath0 m absorption feature of carbonaceous dust grains should be detected ( pendleton et al . 1994 ) , as found in some sources ( e.g. , ngc 1068 ; bridger , wright & geballe 1994 ; imanishi et al .
if a source is powered by unobscured agn activity , a 34 @xmath0 m spectrum should be nearly featureless .
figure 1 summarizes the expected spectral shapes for the above three cases and for the cases of composites of agn and starburst activity . finally , since a continuum level both longward and shortward of these features at 3.33.4 @xmath0 m is observable in the @xmath14-band ( 2.84.2 @xmath0 m ) atmospheric window of earth for nearby ( _ z _
@xmath15 0.18 ) sources , there is no serious uncertainty in determining a continuum level .
ground - based spectroscopy usually makes use of small slit widths , typically less than a few arcsec .
one arcsec corresponds to greater than a few 100 pc in irlgs at z @xmath11 0.02 .
there is currently an increasing amount of evidence based on high spatial resolution imaging studies of infrared 225 @xmath0 m continuum emission that the bulk of the huge dust emission luminosities of irlgs come from compact ( @xmath15 a few 100 pc ) regions powered either by agn and/or compact starburst activity , with little contribution from extended regions ( surace & sanders 1999 ; surace , sanders , & evans 2000 ; scoville et al .
2000 ; soifer et al .
thus , ground - based 34 @xmath0 m spectroscopy can fully diagnose the properties of dust emission from the energetically important compact regions in most irlgs .
the present study reports ground - based spectroscopy in the wavelength range 34 @xmath0 m of nine irlgs to investigate their energy sources .
to achieve signal - to - noise ratios higher than 1020 with a few hours of on - source integration time using the 3.8-m united kingdom infrared telescope ( ukirt ) , irlgs brighter than 11 mag at @xmath14 or @xmath16 were selected .
the targets are summarized in table 1 .
a total of nine irlgs were observed .
far - infrared luminosities are higher than @xmath310@xmath5 ergs s@xmath6 .
in addition to the irlgs , ngc 253 , a nearby ( @xmath33 mpc ; tully 1988 ) galaxy , was also observed .
the emission properties of ngc 253 can be explained by starburst activity only ( e.g. , rieke et al .
1980 ; engelbracht et al . 1998 ) , and no strong agn activity has been found . although its far - infrared luminosity ( table 1 ) is more than an order of magnitude smaller than those of the irlgs , this source is powered by compact ( @xmath3100 pc ) starburst activity ( kalas & wynn - williams 1994 ) .
ngc 253 is thus an important nearby source for studying the nature of compact starburst activity found in irlgs , and its 34 @xmath0 m spectrum can be used as a template of compact starburst activity .
we used the cooled grating spectrometer ( cgs4 ; mountain et al .
1990 ) to obtain 34 @xmath0 m spectra of the irlgs and ngc 253 with ukirt on mauna kea , hawaii .
an observing log is summarized in table 2 .
sky conditions were photometric throughout the observations .
the detector was a 256 @xmath17 256 insb array .
the 40 l mm@xmath6 grating with 2 pixels wide slit (= 1@xmath18 ) was used .
the resulting spectral resolution was @xmath3750 at 3.5 @xmath0 m .
spectra were obtained toward the flux peak at 34 @xmath0 m .
a standard telescope nodding technique was employed along the slit to subtract the signal from the sky .
the nodding amplitude was @xmath312@xmath19 except for ngc 253 , for which it was @xmath324@xmath19 .
the position angles of the slit used are summarized in table 2 . for irlgs that have clear double nuclei with a small separation ( @xmath15 several arcsec ; arp 220 , ngc 6240 , mrk 273 , and iras 08572 + 3915 ; scoville et al .
2000 ) , the position angles were set so that signals from both the nuclei were observed simultaneously .
a- to g - type standard stars ( table 2 ) were observed with almost the same airmass as the target sources to correct for the transmission of earth s atmosphere . for target sources
whose standard stars are a - type , data points at 3.74 @xmath0 m and 4.05 @xmath0 m were removed from final spectra because these data points are affected by , respectively , strong pf@xmath20 and br@xmath21 absorption lines in the spectra of a - type stars .
the @xmath14- or @xmath16-band magnitudes of standard stars were estimated from their @xmath22-band magnitudes , by adopting the @xmath23 or @xmath24 colors of the stellar types of individual standard stars ( tokunaga 2000 ) .
standard data analysis procedures were employed using iraf .
first , data of bad pixels were replaced with the interpolated signals of the surrounding pixels .
next , bias was subtracted from the obtained frames and the frames were divided by a flat image . then , the spectra of the targets and the standard stars were extracted . since a secondary nucleus was not clearly recognizable as a distinct nucleus in all cases , spectra were extracted by integrating signals over 8@xmath1920@xmath19 along the slit ( table 3 ) . wavelength calibration was performed using an argon or a krypton lamp , and is believed to be accurate within 0.005 @xmath0 m .
the galaxy spectra were divided by those of the standard stars , and were multiplied by the spectra of blackbodies with temperatures corresponding to individual standard stars ( table 2 ) .
after flux calibration based on the adopted standard star fluxes , the final spectra were produced .
the flux - calibrated spectra are shown in figure 2 .
spectra for some irlgs were available in the literature , and the new spectra generally agree with previously published spectra .
detailed comparisons of our spectra with previous ones are found in appendix a.    the comparison of the magnitudes of the cgs4 data with those in the literature is summarized in table 3 .
the magnitudes of irlgs based on the cgs4 data roughly agree with photometric magnitudes measured with larger apertures , suggesting that nearly all of previously measured 34 @xmath0 m emission is also detected in the cgs4 data .
the spectra of ic 694 , arp 220 , ngc 6240 , and mrk 273 are dominated by the 3.3 @xmath0 m pah emission and no detectable 3.4 @xmath0 m carbonaceous dust absorption , as in the case of ngc 253 .
the spectrum of mrk 231 shows the 3.3 @xmath0 m pah emission and possible 3.4 @xmath0 m carbonaceous dust absorption , but the reality of the latter feature is currently doubtful .
the spectrum of iras 05189@xmath102524 shows both 3.3 @xmath0 m pah emission and detectable 3.4 @xmath0 m carbonaceous absorption with an optical depth of @xmath25 @xmath3 0.04 .
the spectrum of iras 08572 + 3915 is dominated by the 3.4 @xmath0 m carbonaceous dust absorption feature ( @xmath25 @xmath3 0.9 ) . neither detectable 3.3 @xmath0 m pah emission nor the 3.4 @xmath0 m carbonaceous dust absorption is found in the spectra of iras 20460 + 1925 and iras 23060 + 0505 . the fluxes and luminosities of the 3.3 @xmath0 m pah emission feature are summarized in table 4 .
in addition to the above features , pf@xmath20 emission lines ( rest wavelength = 3.74 @xmath0 m ) are clearly detected in the spectra of ic 694 and ngc 253 .
the observed fluxes are 2 @xmath17 10@xmath26 ergs s@xmath6 @xmath27 and 6 @xmath17 10@xmath26 ergs s@xmath6 @xmath27 for ic 694 and ngc 253 , respectively . in the spectrum of ngc 6240
, we attribute the flux excess at 3.71 @xmath0 m as h@xmath9 00 s(15 ) emission line ( rest wavelength = 3.626 @xmath0 m ) .
its flux is 1 @xmath17 10@xmath26 ergs s@xmath6 @xmath27 .
ngc 6240 is known to show strong h@xmath9 emission lines whose strength ratios can be approximated by thermal excitation with the temperature of @xmath32000 k ( sugai et al .
the comparison with the h@xmath9 10 s(0 ) flux ( 1.8 @xmath17 10@xmath28 ergs s@xmath6 @xmath27 ) measured with a 3@xmath19 @xmath17 6@xmath19 aperture ( goldader et al .
1995 ) provides the 00 s(15 ) to 10 s(1 ) flux ratio of 0.06 .
this ratio is similar to that found in shock gas in the orion molecular outflow where the h@xmath9 line ratios are roughly fitted by 2000 k thermal emission ( brand et al .
in irlgs , agn activity and compact starburst activity are likely to coexist spatially close to each other .
if energetic ionizing photons from agn activity penetrate the compact starburst regions and destroy pahs ( voit 1992 ) that otherwise should produce the strong 3.3 @xmath0 m emission feature , the 3.3 @xmath0 m pah luminosity could be small in spite of the presence of strong compact starburst activity .
however , in irlgs , the destruction of pahs in the starburst regions by energetic photons from agn activity probably does not occur because of the large column density and high pressures in the interstellar medium ( maloney 1999 ) . in the following discussion
, it is assumed that compact starburst activity in the irlgs produces the 3.3 @xmath0 m pah emission intrinsically as strongly as the star - formation activity in starburst galaxies with lower far - infrared luminosities , namely , in proportion to the far - infrared luminnosities .      for starburst - dominated galaxies , the rest - frame equivalent widths of the 3.3 @xmath0 m pah emission feature ( ew@xmath29 ) do not change significantly as a result of dust extinction because flux attenuation of the 3.3 @xmath0 m pah emission and 34 @xmath0 m continuum is similar .
thus very small ew@xmath30 values are most likely caused by the contribution of agn activity to the 34 @xmath0 m continuum emission .
the rest - frame ew@xmath29 values of the nine irlgs are summarized in table 4 .
the value for ngc 253 is 120 nm , which is nearly a median value for starburst - dominated galaxies ( moorwood 1986 ) .
if the scatter of the values is taken into account and irlgs with ew@xmath29 = 60240 nm are regarded as candidate starburst - dominated galaxies , then only ic 694 , arp 220 , and ngc 6240 have ew@xmath29 values as high as those of starburst galaxies , and are thus candidate starburst - dominated galaxies .
mrk 273 shows a factor of @xmath312 smaller value than those of starburst galaxies and thus agn activity is estimated to contribute @xmath3150% of the 34 @xmath0 m continuum flux .
the remaining five irlgs ( mrk 231 , iras 05189@xmath102524 , iras 08572 + 3915 , iras 23060 + 0505 , and iras 20460 + 1925 ) have ew@xmath29 ratios more than an order of magnitude smaller , indicating that @xmath1190% of the 34 @xmath0 m continuum flux originates in agn activity . since the far - infrared to 34 @xmath0 m flux ratio of dust emission heated by agn activity through a radiative transfer mechanism depends strongly on dust radial density distribution ( ivezic , nenkova & elitzur 1999 ) , conversion from the agn contribution at 34 @xmath0 m to that at far - infrared wavelengths
is not certain .
however , the ew@xmath29 values for these six irlgs ( mrk 273 and the five irlgs ) give no evidence that they are powered predominantly by compact starburst activity .      in @xmath325.1 ,
three irlgs have ew@xmath29 values as high as those of starburst galaxies .
however , large ew@xmath29 values of the irlgs do not necessarily mean that their far - infrared luminosities are powered predominantly by starburst activity .
although highly obscured agn activity can contribute to the observed far - infrared luminosity because of negligible extinction in this wavelength range , the contribution of the highly obscured agn to the 34 @xmath0 m flux can be highly attenuated . as shown in figure 1@xmath33 , the apparent spectral shapes of galaxies powered both by starburst and highly obscured agn activity could be nearly the same as those powered only by starburst activity ( fig .
1@xmath34 ) .    a straightforward means of estimating the energy contribution of compact starburst activity is to estimate the absolute 3.3 @xmath0 m pah luminosity and compare it with the far - infrared luminosity .
the 3.3 @xmath0 m pah emission to far - infrared luminosity ratios ( @xmath1/@xmath2 ) of the irlgs are summarized in table 4 .
the ratios of nearby starburst galaxies with lower far - infrared luminosity are estimated to be @xmath31 @xmath17 10@xmath35 ( mouri et al .
1990 ) , although the scatter in this ratio is difficult to estimate because , in the case of nearby starburst galaxies , the ratios are affected by different aperture sizes between 34 @xmath0 m and far - infrared data ( mouri et al .
the observed @xmath1/@xmath2 ratios for all nine irlgs are less than 3 @xmath17 10@xmath36 , a factor of @xmath113 smaller than those of starburst galaxies . the smaller ratios for the six irlgs with small ew@xmath29 values can be explained naturally by a dominant contribution of agn activity , while a few mechanisms can explain the smaller ratios in the three irlgs with large ew@xmath29 values ( ic 694 , arp 220 , and ngc 6240 ) .
if highly obscured agn activity contributes significantly to the far - infrared luminosity but little to the 34 @xmath0 m flux due to high flux attenuation , then the large ew@xmath29 values and the small @xmath1/@xmath2 ratios can be explained .
this is the first possible mechanism .
a second mechanism is the effect of possible signal loss in the cgs4 data
. the possible signal loss for ngc 6240 could be as high as 0.5 mag or a factor of 1.6 , while that for ic 694 and arp 220 is a factor of @xmath371.2 .
provided that spectral shape outside the cgs4 slit is the same as that inside it , the possible slit loss could decrease the @xmath1/@xmath2 ratios only slightly .
a third possible mechanism is the effect of dust extinction .
the @xmath1/@xmath2 ratios for starburst - dominated irlgs could be smaller than those of starburst galaxies with lower far - infrared luminosities if dust extinction toward compact starburst activity in the irlgs is much higher .
for ic 694 , arp 220 , and ngc 6240 , dust extinction toward starburst regions in the case of a foreground dust screen model could be as high as @xmath38 @xmath3 20 mag ( sugai et al . 1999 ; alonso - herrero et al .
2000 ) , @xmath38 @xmath3 45 mag ( rigopoulou et al . 1999 , fig .
7 ) , and @xmath38 @xmath3 7 mag ( @xmath39 = 0.8 mag ; sugai et al .
1997 ) , respectively . by adopting a standard extinction curve ( @xmath40 @xmath3 0.060.07 @xmath17 @xmath38 ; rieke & lebofsky 1985 ; lutz et al . 1996 ) and taking possible signal loss into account , the _ intrinsic _ @xmath1/@xmath2 ratios after correction for dust extinction could be as high as 47 @xmath17 10@xmath36 for all three .
provided that the @xmath1/@xmath2 ratios of control starburst galaxy samples suffer the effect of dust extinction with @xmath38 @xmath3 10 mag ( rigopoulou et al .
1999 , fig .
7 ) , the intrinsic @xmath1/@xmath2 ratios of starburst - dominated galaxies are estimated to be 1.71.9 @xmath17 10@xmath35 . the intrinsic @xmath1/@xmath2 ratios for the three irlgs are a factor of @xmath33 smaller . for reference , the 3.3 @xmath0 m pah flux of ngc 253 within our slit is measured to be 2.8 @xmath17 10@xmath41 ergs s@xmath6 @xmath27 .
the @xmath14-band magnitude inside our slit is 2.7 mag fainter than the total magnitude of this galaxy ( table 3 ) . if the 34 @xmath0 m spectrum shape outside the slit is the same as that inside it , and if @xmath38 = 14 mag is adopted as dust extinction toward starburst regions in this galaxy ( rigopoulou et al .
1999 , fig .
7 ) , then the intrinsic @xmath1/@xmath2 ratio is 1.31.5 @xmath17 10@xmath35 , consistent with the value of 1.71.9 @xmath17 10@xmath35 within @xmath325% .
although dust extinction for arp 220 was derived based on mid - infrared data , the dust extinction for ic 694 and ngc 6240 was estimated based on near - infrared data .
the dust extinction estimate based on near - infrared data could be smaller than that based on mid - infrared data ; the latter estimate is thought to be more reliable , particularly in the case of high dust extinction ( genzel et al .
1998 ) . from this viewpoint
, it is possible that the intrinsic @xmath1/@xmath2 ratios for ic 694 and ngc 6240 could be higher than our estimate .
further , when dust is mixed with emission sources , dust extinction correction factor could be higher than that for a foreground screen dust model ( e.g. , genzel et al .
since the scatter of the intrinsic @xmath1/@xmath2 ratios for starburst galaxies is also poorly known , the energy source of these three irlgs can not be determined confidently from our diagnoses .
the energy sources of the nine irlgs derived from the 34 @xmath0 m spectra is summarized in table 5 . for ic 694 , arp 220 , and ngc 6240 , which show strong 3.3 @xmath0 m pah emission in equivalent width , it is not clear whether compact starburst activity is the only important energy source or a significant fraction of the luminosity is powered by obscured agn activity ( figs .
1@xmath34 and 1@xmath33 ) .
mrk 231 and mrk 273 , which show weak 3.3 @xmath0 m pah emission and no detectable 3.4 @xmath0 m carbonaceous dust absorption , are suggested to be powered by either unobscured or obscured agn activity ( figs .
1@xmath42 and 1@xmath33 ) .
iras 05189@xmath102524 , which shows weak 3.3 @xmath0 m pah emission and detectable ( but weak ) 3.4 @xmath0 m carbonaceous dust absorption , is thought to be powered by moderately obscured agn activity ( fig .
1@xmath33 ) .
the 34 @xmath0 m spectrum of iras 08572 + 3915 is typical of highly dust obscured agns , and thus this source is most likely powered by highly embedded agn activity ( fig .
1@xmath43 ) , as further supported by more detail study based on 320 @xmath0 m spectroscopy ( dudley & wynn - williams 1997 ; imanishi & ueno 2000 ) .
the featureless 34 @xmath0 m spectra of iras 23060 + 0505 and iras 20460 + 1925 suggest that these sources are powered by less obscured agn activity ( fig .
1@xmath44 ) .
for arp 220 , ngc 6240 , mrk 273 , mrk 231 , iras 05189@xmath102524 , and iras 23060 + 0505 , the 7.7 @xmath0 m pah peak flux to 7.7 @xmath0 m continuum flux ratios were used to diagnose their energy source ( rigopoulou et al . 1999 ; laureijs et al . 2000 ) . these ratios and
the equivalent widths are a similar ( but not identical ) concept , and both provide the pah strength relative to the continuum emission .
the results of the 34 @xmath0 m spectra and the _ iso _ spectra are compared in table 5 . as mentioned in @xmath32 1 , a continuum determination based on the _ iso _
5.811.6 @xmath0 m spectra is not easy . when data at 11 @xmath0 m were used to determine a continuum level ( rigopoulou et al .
1999 ) , the 7.7 @xmath0 m pah flux for sources with a strong silicate dust absorption feature could be overestimated . since a different amount of effect of dust extinction at different wavelengths does not change the relative pah strengths significantly ( @xmath32 5.1 ) , the comparison of the strength of the 3.3 @xmath0 m and 7.7 @xmath0 m pah emission relative to the continuum could provide useful information on the effect of the silicate dust absorption in the _ iso _ spectra .    the relative pah strengths ( ew@xmath29 and 7.7 @xmath0 m pah peak to continuum flux ratio ) for mrk 231 , iras 05189@xmath102524 , and iras 23060 + 0505 are much smaller than those of starburst galaxies in both wavelengths .
the relative pah strength for ngc 6240 is the same within 30% as the typical value of starburst - dominated galaxies at both wavelengths .
in fact , for ngc 6240 , no indication of strong silicate dust absorption was found in ground - based 813 @xmath0 m spectra ( roche et al .
1991 ; dudley 1999 ) .
however , for arp 220 , the 34 @xmath0 m spectrum provides the ew@xmath29 value of @xmath380 nm , slightly smaller than those of starburst - dominated galaxies ( @xmath3120 nm ) , while the _ iso _ spectrum provides the 7.7 @xmath0 m pah peak - to - continuum flux ratio of @xmath34.2 , significantly higher than the average value of starburst galaxies ( @xmath33.0 ) .
also for mrk 273 , the 34 @xmath0 m spectrum gives ew@xmath29 @xmath335 nm , a factor of @xmath33 smaller than starburst galaxies , while the _ iso _ spectrum gives the 7.7 @xmath0 m pah peak - to - continuum flux ratio @xmath31.9 , only @xmath340% smaller . for arp 220 and
mrk 273 , dudley ( 1999 ) argued that silicate dust absorption plays an important role based on ground - based 813 @xmath0 m spectra .
the discrepancy may be caused by the effect of the silicate dust absorption .
the above comparison suggests that , for at least a some fraction of irlgs , the 7.7 @xmath0 m pah flux estimated by rigopoulou et al .
( 1999 ) could be too large because of the use of data affected by silicate dust absorption as a continuum level . in this case
, the presence of the silicate dust absorption is underestimated .
highly obscured agns show spectra characterized by weak pah emission and strong silicate dust absorption .
virtually no irlgs powered by highly obscured agns were found based on the _ iso _ data ( genzel et al .
1998 ; rigopoulou et al .
1999 ) , which may be the result of the determination of the continuum level .    finally , for iras 05189@xmath102524 , the flux and equivalent width of the 3.3 @xmath0 m pah emission estimated by an _ iso _
spectrum obtained using a larger ( 24@xmath19 @xmath17 24@xmath19 ) aperture ( clavel et al .
2000 ) are a factor of @xmath36 larger than the estimate in this paper .
if the larger values are the case , the bulk of the 3.3 @xmath0 m pah emission detected in the _ iso _ spectrum must originate in extended star - forming regions outside the cgs4 slit .
the extended star - formation regions could explain nearly all of the far - infrared luminosity of this source , if the value of @xmath1/@xmath2 = 1.71.9 @xmath17 10@xmath35 is adopted .
however , this source shows no extended continuum emission at 225 @xmath0 m ( scoville et al .
2000 ; soifer et al .
if the extended regions are fully responsible for the far - infrared emission at 60 @xmath0 m but contribute to , say , less than 5% of the 25 @xmath0 m flux ( soifer et al . 2000 ) , then the 60 @xmath0 m to 25 @xmath0 m flux ratio in jy is @xmath1180 , much larger than any nearby star - forming galaxies ( rice et al .
1988 ) . since the physical conditions in the extended regions of iras 05189@xmath102524
( not nuclear compact regions that could have an extreme environment in density and so on ) are expected not to be quite different from those in nearby star - forming galaxies , such a sudden flux increase from the mid - infrared ( @xmath1525 @xmath0 m ) to far - infrared ( @xmath1140 @xmath0 m ) only in the case of iras 05189@xmath102524 is very unlikely .
while the 3.3 @xmath0 m pah flux in this paper is estimated by tracing the spectral profile ( see the caption of table 4 for more details ) , clavel et al . simply sum up all the flux above the adopted continuum level in the pre - defined wavelength range ( j. clavel , private communication ) .
the _ iso _
spectrum provided by dr .
clavel ( private communication ) shows a flux excess above an adopted continuum level at @xmath33.3 @xmath0 m , but it has lower spectral resolution ( r @xmath3 100 ) and larger scatter in data points at 34 @xmath0 m than the cgs4 spectrum .
furthermore , although the wavelength of the 3.3 @xmath0 m pah emission feature usually extends from 3.24 @xmath0 m to 3.34 @xmath0 m in the rest - frame ( tokunaga et al . 1991 ; this work ) , the flux excess in the _ iso _ spectrum is found from 3.25 @xmath0 m to 3.40 @xmath0 m ( not 3.34 @xmath0 m ) . if this excess is actually due to the 3.3 @xmath0 m pah emission feature , its spectral shape in the extended regions of iras 05189@xmath102524 is quite different from typical ones .
more sensitive and higher spectral resolution spectroscopy with a wide aperture is necessary to address this issue .
we tentatively consider that the extended regions contribute little to the total far - infrared luminosity of iras 05189@xmath102524 .
soifer et al . ( 2000 ) suggested based on their 825 @xmath0 m imaging data that agn activity is the dominant energy source of mrk 231 , iras 05189@xmath102524 , and iras 08572 + 3915 .
our suggestions are the same for these three sources .
soifer et al . were unable to distinguish the energy source of mrk 273 and arp 220 .
we argue that agn activity is the main energy source of mrk 273 , while the energy source of arp 220 is not clear . as a whole ,
the suggestions in this paper are consistent with those by soifer et al .      in the @xmath152 @xmath0 m spectra of mrk 231 , iras 05189@xmath102524 , iras 23060 + 0505 , and iras 20460 + 1925 ,
a broad ( @xmath112000 km s@xmath6 in full width at half - maximum ) component of neutral hydrogen recombination line emission was detected , and their intrinsic luminosities are high enough to account for the far - infrared luminosities resulting from agn activity alone ( veilleux , sanders , & kim 1999b ) .
our suggestions about the energy source for these four irlgs are consistent with theirs .      among the nine irlgs ,
five ( mrk 231 , iras 05189@xmath102524 , iras 08572 + 3915 , iras 23060 + 0505 , and iras 20460 + 1925 ) have warm ( _ iras _ 25 @xmath0 m to 60 @xmath0 m flux ratio @xmath46/@xmath47 @xmath11 0.2 ) far - infrared color ( table 1 ) .
the above three infrared diagnoses have provided a growing amount of evidence that these warm irlgs are powered predominantly by agn activity ( table 5 ) . however , warm irlgs constitute only @xmath320% of total irlgs , and most irlgs have cool ( @xmath46/@xmath47 @xmath15 0.2 ) far - infrared color ( kim & sanders 1998 ) . it was argued based on the estimate of the 7.7 @xmath0 m pah strengths in the _ iso _ spectra that these cool irlgs are powered predominantly by starburst activity , with little agn contribution ( genzel & cesarsky 2000 ) . among the nine irlgs ,
three ( arp 220 , ngc 6240 , and mrk 273 ) are cool irlgs and are classified as starburst - dominated irlgs based on the 7.7 @xmath0 m pah strength ( table 5 ) . however , for arp 220 and mrk 273 , our 34 @xmath0 m spectra suggest that the 7.7 @xmath0 m pah strengths , and thereby the magnitudes of starburst activity , may be overestimated based on the _ iso _ spectra . in fact , for mrk 273 , although the relative 7.7 @xmath0 m pah emission strength suggests this source is powered by starburst activity ( rigopoulou et al .
1999 ) , a different diagnosis based on high excitation lines in the _ iso _ spectra indicates that agn activity is a dominant energy source ( genzel et al .
1998 ) . for ngc 6240
, x - ray data suggest that highly embedded and highly luminous agn activity is an important energy source ( vignati et al .
the contribution of agn activity in the cool irlgs may be higher than the estimate based on the 7.7 @xmath0 m pah emission strength in the _ iso _ spectra .
high spatial resolution radio images have been used to investigate the relative importance between agn and starburst activity in irlgs . a compact radio core , which was thought to be strong evidence for agn activity , is no longer definitive evidence for agn activity because both agn and compact starburst activity could explain observed properties of the compact radio core ( smith , lonsdale , & lonsdale 1998 ) .
however , irlgs whose radio emission morphology is extended without a compact core are thought to be powered by starburst activity .
compact radio cores were detected in ic 694 , arp 220 , ngc 6240 , mrk 273 , mrk 231 , and iras 20460 + 1925 ( lonsdale , lonsdale , & smith 1992 ; lonsdale , smith , & lonsdale 1993 ; heisler et al .
1998 ; kewley et al .
2000 ) , but not in iras 05189 + 2524 , iras 08572 + 3915 , and iras 23060 + 0505 ( lonsdale et al .
1993 ; heisler et al .
it is puzzling that no radio core emission was detected in the three warm irlgs , for which more than one infrared diagnoses consistently suggest that agn activity is a dominant energy source ( table 5 ) .      a high fraction ( 6/9 ) of irlgs in this sample
are suggested to be powered by agn activity .
this is thought to be caused by the selection effect of picking irlgs that are bright at 34 @xmath0 m . at a given far - infrared flux , the 34 @xmath0 m flux of less obscured
agns is generally higher than that of starburst - dominated galaxies ( schmitt et al .
1997 ) . among the nine irlgs , 5/9 ( 55% )
are optically classified as seyfert types ( table 1 ) and 4/9 ( 44% ; mrk 231 , iras 05189@xmath102524 , iras 23060 + 1925 , and iras 23060 + 0505 ) show a detectable broad emission line at @xmath152 @xmath0 m ( that is , agn activity is less obscured ) .
these fractions are significantly higher than those in a larger sample of irlgs ( @xmath320% ; veilleux , kim , & sanders 1999a ; veilleux et al .
furthermore , the fraction of warm irlgs in this sample ( 5/9 ; 55% ) is much higher than that in a larger sample of irlgs ( @xmath320% ; kim & sanders 1998 ) .
observations of a larger number of fainter irlgs at @xmath14 are clearly necessary to evaluate the relative role between agn and starburst activity in irlgs as a whole .
the energy sources of nine irlgs were diagnosed using the ew@xmath29 values and the @xmath1/@xmath2 ratios .
the following main results were found :    1 .   mrk 231 , mrk 273 , and iras 05189@xmath102524 show detectable 3.3 @xmath0 m pah emission in their 34 @xmath0 m spectra . however , the pah emission strengths are too small to account for the bulk of their far - infrared luminosities by starburst activity , indicating that these sources are powered by agn activity .
no pah emission was detected in the spectra of iras 08572 + 3915 , iras 23060 + 0505 , and iras 20460 + 1925 .
it was suggested that iras 08572 + 3915 is powered by highly obscured agn activity , while iras 23060 + 0505 and iras 20460 + 1925 are powered by less - obscured agn activity .
3 .   for ic 694 ,
arp 220 , and ngc 6240 , while the ew@xmath29 values are as large as those of starburst galaxies , the @xmath1/@xmath2 ratios after correction for screen dust extinction are a factor of @xmath33 smaller . whether these three irlgs are powered by starburst activity only or obscured
agn activity contributes significantly is not clear because actual dust extinction correction factor could be higher than our estimate and the scatter of the intrinsic @xmath1/@xmath2 ratios for starburst activity is poorly known .
our energy diagnoses reinforced the previous suggestions that agn activity is the dominant energy source in warm irlgs .
we thank dr .
t. kerr , dr .
j. davies , t. carroll , and t. wold for their support during the ukirt observing runs .
we are grateful to dr .
d. b. sanders and the anonymous referee for their useful comments .
a. t. tokunaga and h. ando give mi the opportunity to work at the university of hawaii . l.
good kindly proofread this paper , and k. teramura modified fig .
the united kingdom infrared telescope is operated by the joint astronomy centre on behalf of the u.k . particle physics and astronomy research council .
mi is financially supported by the japan society for the promotion of science during his stay at the university of hawaii .
ccd gratefully acknowledges the national research council ( usa ) , the naval research laboratory , and the office of naval research for their generous support , and the sponsorship of dr .
j. fischer .
our new spectra ( mostly higher signal - to - noise ratios and spectral resolution ) are compared with old spectra in the literature .      a flux - calibrated 34 @xmath0 m spectrum taken with a 2@xmath487 aperture
was presented by dennefeld & desert ( 1990 ) . the overall spectral shape and flux level in the new spectrum
are consistent with the previous spectrum . using narrow band filters , satyapal et al .
( 1999 ) estimated the flux of the 3.3 @xmath0 m pah emission to be 8.5 @xmath17 10@xmath28 ergs s@xmath6 @xmath27 within a 7@xmath19 aperture .
this is @xmath340% higher than our estimate , but this difference does not affect our discussion significantly .
low - resolution spectra taken using a circularly variable filter with a 8@xmath487 aperture were presented by rieke et al .
they showed strong 3.3 @xmath0 m pah emission .
the new spectra show slightly smaller continuum flux levels but similar spectral shapes .      a flux - calibrated 34 @xmath0 m spectrum taken with a 3@xmath488 aperture was presented by imanishi et al .
the new spectrum provides much higher signal - to - noise ratios .
the spectral shape and flux level are consistent with the previous spectrum .
a flux - uncalibrated 34 @xmath0 m spectrum was presented by wright et al .
it showed weak 3.3 @xmath0 m pah emission and weak 3.4 @xmath0 m carbonaceous dust absorption with an optical depth of @xmath30.05 .
the new spectrum is consistent with the previous spectrum in the spectral shape and the optical depth of the carbonaceous dust absorption .
a flux - uncalibrated 34 @xmath0 m spectrum was presented by wright et al .
it was dominated by the 3.4 @xmath0 m carbonaceous dust absorption feature , and its optical depth was estimated to be @xmath30.9 ( pendleton 1996 ) .
both the spectral shape and the optical depth in the new spectrum are consistent with the previous spectrum .
a low - resolution spectrum taken using a circularly variable filter with a 7@xmath485 aperture was presented by moorwood ( 1986 ) .
both the new and old spectra show the pah emission at 3.3 @xmath0 m and 3.4 @xmath0 m .
the flux level in our new spectrum is smaller than the old spectrum because of smaller aperture size .
alonso - herrero , a. , rieke , g. h. , rieke , m. j. , & scoville , n. z. 2000 , apj , 532 , 845 armus , l. , heckman , t. m. , & miley , g. k. 1989 , apj , 347 , 727 blain , a. w. , kneib , j. -p . ,
ivison , r. j. , & smail , i. 1999 , apj , 512 , l87 brand , p. w. j. l. , moorhouse , a. , burton , m. g. , geballe , t. r. , bird , m. , & wade , r. 1988 , apj , 334 , l103 bridger , a. , wright , c. s. , & geballe , t. r. 1994 , in mclean i. , ed . , infrared astronomy with arrays : the next generation , kluwer , dordrecht , p.537 clavel , j. , et al .
2000 , a&a , 357 , 839 dennefeld , m. , & desert , f. x. 1990 , a&a , 227 , 379 dudley , c. c. 1999 , mnras , 307 , 553 dudley , c. c. , wynn - williams c. g. 1997 , apj , 488 , 720 elbaz , d. , et al . 1998 ,
proceedings 34th liege international astrophysics colloquium ` the next generation space telescope : science drivers and technological challenges ' liege , belgium , p.47 , ( astro - ph/9807209 ) engelbracht , c. w. , rieke , m. j. , rieke , g. h. , kelly , d. m. , & achtermann , j. m. 1998 , apj , 505 , 639 frogel , j. a. , gillett , f. c. , terndrup , d. m. , & vader , j. p. 1989 , apj , 343 , 672 genzel , r. , & cesarsky , c. j. 2000 , ara&a , in press ( astro - ph/0002184 ) genzel , r. , weitzel , l. , tacconi - garman , l. e. , blietz , m. , cameron , m. , krabbe , a. , & lutz , d. 1995 , apj , 444 , 129 genzel , r. , et al .
1998 , apj , 498 , 579 glass , i. s. , & moorwood , a. f. m. 1985 , mnras , 214 , 429 goldader , j. d. , joseph , r. d. , doyon , r. , & sanders , d. b. 1995 , apj , 444 , 97 heisler , c. a. , norris , r. p. , jauncey .
d. l. , reynolds , j. e. , & king , e. a. 1998 , mnras , 300 , 1111 hill , g. j. , wynn - williams , c. g. , & becklin , e. e. 1987 , apj , 316 , l11 imanishi , m. , & ueno , s. 2000 , apj , 535 , 626 imanishi , m. , terada , h. , goto , m. , & maihara , t. 1998 , pasj , 50 , 399 imanishi , m. , terada , h. , sugiyama , k. , motohara , k. , goto , m. , & maihara , t. 1997 , pasj , 49 , 69 ivezic , z. , nenkova , m. , & elitzur , m. 1999 , user manual for dusty , university of kentucky internal report , accessible at http://www.pa.uky.edu/@xmath49moche/dusty/ ( astro - ph/9910475 ) joy , m. , telesco , c. m. , decher , r. , lester , d. f. , harvey , p. m. , rickard , l. j. , & bushouse , h. 1989 , apj , 339 , 100 kalas , p. , & wynn - williams , c. g. 1994 , apj , 434 , 546 kewley , l. j. , heisler , c. a. , dopita , m. a. , sutherland , r. , norris , r. p. , reynolds , j. , & lumsden , s. 2000 , apj , 530 , 704 kim , d. -c .
, & sanders , d. b. 1998 , apjs , 119 , 41 laureijs , r. j. , et al .
2000 , a&a , in press ( astro - ph/0005076 ) lonsdale , c. j. , lonsdale , c. j. , & smith , h. e. 1992 , apj , 391 , 629 lonsdale , c. j. , smith , h. e. , & lonsdale , c. j. 1993 , apj , 405 , l9 lutz , d. , et al .
1996 , a&a , 315 , l269 lutz , d. , spoon , h. w. w. , rigopoulou , d. , moorwood , a. f. m. , genzel , r. 1998 , apj , 505 , l103 maloney , p. r. 1999 , ap&ss , in press , the ringberg workshop , `` ultraluminous galaxies : monsters or babies '' ( astro - ph/9903275 ) mathis , j. s. 1990 , ara&a , 28 , 37 moorwood , a. f. m. 1986 , a&a , 166 , 4 mountain , c. m. , robertson , d. j. , lee , t. j. , & wade , r. 1990 , proc .
spie , 1235 , 25 mouri , h. , kawara , k. , taniguchi , y. , & nishida , m. 1990 , apj , 356 , l39 pendleton , y. 1996 , in the cosmic dust connection , greenberg , j. m. , ed . , kluwer academic publishers , p.71 pendleton , y. j. , sandford , s. a. , allamandola , l. j. , tielens , a. g. g. m. , & sellgren , k. 1994 , apj , 437 , 683 rice , w. , lonsdale , c. j. , soifer , b. t. , neugebauer , g. , kopan , e. l. , lloyd , l. a. , de jong , t. , & habing , h. j. 1988 , apjs , 68 , 91 rieke , g. h. , & lebofsky , m. j. 1985 , apj , 288 , 618 rieke , g. h. , curti , r. m. , black , j. h. , kailey , w. f. , mcalary , c. w. , lebofsky , m. j. , & elston , r. 1985 , apj , 290 , 116 rieke , g. h. , lebofsky , m. j. , thompson , r. i. , low , f. j. , & tokunaga , a. t. 1980 , apj , 238 , 24 rigopoulou , d. , spoon , h. w. w. , genzel , r. , lutz , d. , moorwood , a. f. m. , & tran , q. d. 1999 , aj , 118 , 2625 roche , p. f. , aitken , d. k. , smith , c. h. , & ward , m. j. 1991 , mnras , 248 , 606 rudy , r. j. , levan , p. d. , & rodriguez - espinosa , j. m. 1982 , aj , 87 , 598 sanders , d. b. , & mirabel , i. f. 1996 , ara&a , 34 , 749 satyapal , s. , watson , d. m. , forrest , w. j. , pipher , j. l. , fischer , j. , greenhouse , m. a. , smith , h. a. , & woodward , c. e. 1999 , apj , 516 , 704 schmitt , h. r. , kinney , a. l. , calzetti , d. , & bergmann , t. s. 1997 , aj , 114 , 592 scoville , n. z. , et al .
2000 , aj , 119 , 991 smith , h. e. , lonsdale , c. j. , & lonsdale , c. j. 1998 , apj , 492 , 137 soifer , b. t. , sanders , d. b. , madore , b. f. , neugebauer , g. , danielson , g. e. , elias , j. h. , lonsdale , c. j. , & rice , w. l. 1987 , apj , 320 , 238 soifer , b. t. , et al .
2000 , aj , 119 , 509 sugai , h. , malkan , m. a. , ward , m. j. , davies , r. i. , & mclean , i. s. 1997 , apj , 481 , 186 sugai , h. , davies , r. i. , malkan , m. a. , mclean , i. s. , usuda , t. , & ward , m. j. 1999 , apj , 527 , 778 surace , j. a. , & sanders , d. b. 1999 , apj , 512 , 162 surace , j. a. , sanders , d. b. , & evans , a. s. 2000 , apj , 529 , 170 tokunaga , a. t. , sellgren , k. , smith , r. g. , nagata , t. , sakata , a. , & nakada , y. 1991 , apj , 380 , 452 tokunaga , a. t. 2000 , in allen s astrophysical quantities , ed .
a. n. cox ( 4th ed : aip press : springer ) , chapter 7 , p.143 tully , r. b. 1988 , nearby galaxies catalog , cambridge university press , cambridge , p.10 vader , j. p. , frogel , j. a. , terndrup , d. m. , & heisler , c. a. 1993 , aj , 106 , 1743 veilleux , s. , kim , d. -c .
, & sanders , d. b. 1999a , apj , 522 , 113 veilleux , s. , kim , d. -c . , sanders , d. b. , mazzarella , j. m. , & soifer , b. t. 1995 , apjs , 98 , 171 veilleux , s. , sanders , d. b. , & kim , d. -c
1999b , apj , 522 , 139 vignati , p. , et al .
1999 , a&a , 349 , l57 voit , g. m. 1992 , mnras , 258 , 841 wright , g. s. , bridger , a. , geballe , t. r. , & pendleton , y. 1996 , new extragalactic perspectives in the new south africa , kluwer academic publishers , block , d. l. , & greenberg , j. m. , ( eds ) , p.143 zhou , s. , wynn - williams , c. g. , & sanders , d. b. 1993 , apj , 409 , 149       note .
( 1 ) : object name . col . ( 2 ) : energy source derived by this work
. `` agn '' and `` sb '' mean agn - dominated and starburst - dominated , respectively .
`` sb+agn '' means both kinds of activity contribute roughly equally .
( 3 ) : energy source derived based on the ratio between the 7.7 @xmath0 m pah peak flux and 7.7 @xmath0 m continuum flux .
the ratios are shown in parentheses .
adopted from rigopoulou et al .
( 1999 ) except for iras 05189@xmath102524 for which the ratio is adopted from laureijs et al .
the average ratios for starburst galaxies and agns are @xmath33.0 and @xmath30.2 , respectively ( rigopoulou et al . 1999 ) . those with the ratios of @xmath111.0 and @xmath151.0 are classified as `` agn '' and `` sb '' , respectively ( rigopoulou et al . 1999 ) .
( 4 ) : energy source derived from high spatial resolution 825 @xmath0 m imaging ( soifer et al .
( 5 ) : energy source derived from 2 @xmath0 m spectroscopy ( veilleux et al .
1999b ) . those with a detected broad emission line whose luminosity is high enough to account for the far - infrared luminosity by agn activity are classified as `` agn ''