from being a scientific curiosity only a few years ago , energy harvesting ( eh ) is well on its way to becoming a game - changing technology in the field of self - sustainable , autonomous wireless networked systems .
a major factor that has contributed to the growth of the eh market is the evolution of ultra - low power electronics , which can run on the minuscule amounts of energy supplied by typical solar , vibration or thermal energy harvesters @xcite .
a number of companies are already offering system solutions consisting exclusively of eh sensor nodes @xcite . stimulated by these advances , eh - powered wireless communications have attracted growing interest in recent years @xcite .
different from traditional communication systems , eh from environmental sources shifts the paradigm on resource allocation from reducing energy consumption to the most efficient utilization of opportunistic energy .
energy availability constraints are imposed such that the energy accumulatively consumed up to any time can not exceed what has been accumulatively harvested so far .
considering these new type of constraints , optimal transmission policies were characterized for point - to - point channels in @xcite , for broadcast channels in @xcite , for multi - access channels in @xcite , for interference channels in @xcite , for two - hop relay channels in @xcite , for transmitters with non - ideal circuit - power consumption @xcite , and for systems with battery imperfections in @xcite .
existing works @xcite on eh communications mostly addressed offline optimizations , where the eh profiles were assumed to be known _ a priori_. in practical scenarios , complete predictability of eh profiles is clearly an over - simplified assumption .
relying on past realizations of eh processes and certain statistics of their future evolutions , @xcite developed some heuristic online algorithms , which , however , lack strong analytical performance guarantees . by modeling the eh and/or data processes as markov processes
, online optimizations were cast as markov decision problems ( mdp ) and numerically solved with dynamic programming tools in @xcite .
however , the well - known `` curse - of - dimensionality '' with such solutions precludes their application for all but the simplest practical networks .
leveraging stochastic optimization tools , a few low - complexity online schemes were developed in @xcite .
considering the energy - aware routing with energy replenishment at nodes , @xcite proposed an algorithm that achieves a logarithmic competitive ratio and is asymptotically optimal as the network size grows . a simple asymptotically optimal joint energy allocation and routing scheme was developed for rechargeable sensor networks with static and non - interfering ( i.e. , orthogonal ) links in @xcite . relying on the lyapunov optimization techniques , @xcite developed and analyzed the utility - optimal resource scheduling schemes for general eh - powered multi - hop wireless networks , while @xcite addressed the effect of finite energy and data storage capacities on such resource allocation tasks .
the schemes in @xcite assumed ideal energy storage devices ( i.e. , batteries ) in use . under this assumption
, the energy - queue sizes at the batteries can play the role of `` stochastic '' lagrange multipliers to develop a dual - subgradient based solver to the intended problems .
however , the imperfections with practical batteries could disable this approach . in this paper
, we consider a practical battery model accounting for finite battery capacity , energy ( dis-)charging loss , and energy dissipation over time . by integrating and generalizing the lyapunov optimization techniques in @xcite , we re -
establish a systematic framework to develop and analyze the stochastic online control schemes for eh wireless networks with such imperfect batteries .
specifically , we propose a data - backpressure based scheduling and degenerated energy - queue based power allocation scheme that can maximize the utility of data rates for eh multi - hop wireless networks , without requiring any statistical knowledge of the stochastic channel , data - traffic , and eh processes .
different from @xcite where an eh admission mechanism is performed to ensure finite energy queues , we apply the sample path analysis in @xcite to derive the conditions that the proposed scheme is feasible for any given finite battery capacities without eh admission , which can help fully exploit the available harvested energy .
in addition , we rigorously establish the performance guarantees of the proposed scheme in form of sub - optimality bounds in the presence of practical battery imperfections .
numerical results demonstrate that the proposed scheme can significantly outperform the existing alternatives .
the rest of the paper is organized as follows .
the system models are described in section ii .
the proposed dynamic resource management scheme is developed and analyzed in section iii .
numerical results are provided in section iv , followed by concluding remarks in section v.
as with @xcite , we consider a general eh multi - hop wireless network that operates in slotted time . for convenience ,
the slot duration is normalized to unity ; thus , the terms `` energy '' and `` power '' can be sometimes used interchangeably .
the network is represented by a graph @xmath0 , where @xmath1 denotes the set of network nodes , and @xmath2 , \ ; n , m \in { \cal n}\}$ ] collects the directed links between nodes . for each node @xmath3 , define two sets of neighbor nodes @xmath4 \in { \cal l}\}$ ] , and @xmath5 \in { \cal l}\}$ ] . further define @xmath6 as the maximum in- and out - degree for nodes in the network .
the network delivers packets for data flows ( called commodity in @xcite ) indexed by their destination nodes @xmath7 . per time
slot @xmath8 , a data admission is implemented by the network to decide the number @xmath9 of packets for flow @xmath7 that can be newly admitted at node @xmath10 .
we assume that @xmath11 with some finite @xmath12 at all time .
every node in the network is assumed to be capable of harvesting energy from the environmental ( e.g. , renewable ) sources to power its transmissions .
the amount of harvested energy is clearly random and intermittent over time .
let @xmath13 denote the amount of harvested energy by node @xmath10 at time slot @xmath8 , and let @xmath14 $ ] be called the _ energy state _ at @xmath8 .
we assume that @xmath15 takes values in some finite set and there exists @xmath16 such that @xmath17      per slot @xmath8 , let @xmath18 denote time - varying , random _ channel state _ , which in general can be a @xmath19-by-@xmath19 matrix where the @xmath20 component denotes the channel condition between nodes @xmath10 and @xmath21 .
we assume that @xmath18 take values in some finite set for all time .
given @xmath18 , the network allocates a power vector @xmath22}(t ) , \ ; \forall [ n , m ] \in { \cal l}]$ ] for data transmissions over links , where @xmath23}(t)$ ] denotes the power allocated to node @xmath10 for link @xmath24 $ ] at time @xmath8 .
we assume that each node has a peak power constraint such that @xmath25}(t ) \leq p_{\max } , \quad \forall n\ ] ] with a constant @xmath26 at all time .
given the channel state @xmath18 and the power allocation vector @xmath27 , the transmission rate over the link @xmath24 $ ] is dictated by a rate - power function @xmath28}(t ) = \mu_{[n , m ] } ( \boldsymbol{s}(t ) , \boldsymbol{p}(t)).\ ] ] for any bounded power allocation @xmath27 , we also assume that there exists a finite constant @xmath29 such that @xmath30}(t ) \leq \mu_{\max } , \quad \forall [ n , m ] \in { \cal l}\ ] ] for all time under any channel state @xmath18 .
now let @xmath31}^c(t)$ ] denote the rate allocated to the data flow @xmath7 over link @xmath24 $ ] at time @xmath8 .
it is clear that we have : @xmath32}^c(t ) \leq \mu_{[n , m]}(t ) , \quad \forall [ n , m].\ ] ]        every node in the network has an energy storage device , i.e. , battery , to save the harvested energy .
consider a practical battery with : i ) a finite capacity , ii ) ( dis-)charging loss , and iii ) energy degeneration .
let @xmath41 denote the battery capacity , @xmath42 $ ] the ( dis-)charging efficiency ( e.g. , @xmath43 means that only 90% of the charged or discharged energy is useful while other is conversion loss ) , and @xmath44 $ ] the storage efficiency ( e.g. , @xmath45 means that 10% of the stored energy will be `` leaked '' over a slot , even in the absence of discharging ) .    we can model the battery using an energy queue .
let @xmath46 denote the energy queue size , which indicates the amount of the energy left in the battery of node @xmath10 at time @xmath8 ; and let @xmath47 $ ] .
as the data transmissions are powered by the harvested energy stored in the batteries , the power allocation vector @xmath27 must satisfy the following `` energy availability '' constraint : @xmath48}(t ) \leq \xi \eta e_n(t ) , \quad \forall n\ ] ] where the product @xmath49 captures the discharging loss and energy degeneration .    due to the energy availability constraint ( [ eq.ea ] )
, the energy queue @xmath46 evolves according to the dynamic equation : @xmath50}(t)}{\xi } + \xi e_n(t ) , \label{eq.e}\\      & 0 \leq e_n(t ) \leq e_{\max } \label{eq.emax}\end{aligned}\ ] ] with @xmath51 , @xmath52 .
define the time - average rate for data flow @xmath7 that is admitted into node @xmath10 , as @xmath53 where the expectation are taken over all sources of randomness .
each flow @xmath7 is associated with a utility function @xmath54 , which is assumed to be increasing , continuously differentiable , and strictly concave .
let @xmath55 denote the maximum first derivative of @xmath56 , and define @xmath57 which is assumed to be finite .    note that the energy state @xmath15 and the channel state @xmath18 are random processes in our model .
the eh wireless network to be controlled is thus a stochastic system .
the goal is to design an online resource management scheme that chooses the data admission amounts @xmath58 $ ] , the power allocations @xmath59}(t ) , \forall [ n , m]]$ ] , as well as the routing and scheduling decisions @xmath60}^c(t),\forall [ n , m],\forall c]$ ] per slot @xmath8 , so as to maximize the aggregate utility of time - average data rates subject to ( s. t. ) network operation constraints . upon defining @xmath61
, we wish to solve @xmath62 note that here the constraints ( [ eq.rn])([eq.emax ] ) are in fact implicitly required to hold for any realization of the underlying random state @xmath63 per slot @xmath8 .
the problem ( [ eq.prob ] ) is a stochastic optimization problem .
the problem is challenging as the optimization variables are coupled over time due to the queue dynamics and energy availability constraints in ( [ eq.q])([eq.e ] ) . such problems usually have to be solved by dynamic programming ( dp ) , which suffers from a curse of dimensionality and requires significant knowledge of network statistics .
differently , we next integrate and generalize the lyapunov optimization techniques in @xcite to develop a low - complexity online control algorithm , which can be proven to yield a feasible near - optimal solution for ( [ eq.prob ] ) under conditions , without requiring any statistical knowledge of stochastic eh and channel processes .      to start , we assume that the rate - power function in ( [ eq.mu ] ) satisfies the following two properties for any given channel state @xmath64 :    for two power allocation vectors @xmath65 and @xmath66 , where @xmath66 is obtained by changing any single component @xmath23}$ ] to zero , we have :    * @xmath67}(\boldsymbol{s } , \boldsymbol{p } ) \leq \mu_{[n',m']}(\boldsymbol{s } , \boldsymbol{p}')$ ] , @xmath68 \neq [ n , m]$ ] ; * @xmath69}(\boldsymbol{s } , \boldsymbol{p } ) - \mu_{[n , m]}(\boldsymbol{s } , \boldsymbol{p } ' ) \leq \delta_1 p_{[n , m]}$ ] , for a finite constant @xmath70 .    for two power allocation vectors
@xmath65 and @xmath66 , where @xmath23}'=p_{[n , m]}+\frac{\delta p}{|{\cal n}_n^o|}$ ] , @xmath71 , and @xmath72}'=p_{[n',m]}$ ] , @xmath73 , we have :    * @xmath31}(\boldsymbol{s } , \boldsymbol{p } ) \leq \mu_{[n , m]}(\boldsymbol{s } , \boldsymbol{p}')$ ] , @xmath74 ; * @xmath75}(\boldsymbol{s } , \boldsymbol{p } ) - \mu_{[n',m]}(\boldsymbol{s } , \boldsymbol{p } ' ) ] \leq \delta_2 \delta p$ ] , for a finite constant @xmath76 .
property 1 was one of the keys for asymptotic optimality analysis of the esa scheme in @xcite . here
, properties 1 and 2 will be the keys for our feasibility and optimality gap analysis in the sequel .
these properties are actually satisfied by most rate - power functions , e.g. , when the rate function has finite directional derivatives with respect to power , and the rates do not improve with increased interference . to see it ,
two intriguing examples are provided below .
_ example 1 ( interference channels ) _ : consider the general interference channel case .
let @xmath77}$ ] denote the channel coefficient from node @xmath10 to node @xmath21 ; i.e. , @xmath78 } , \forall n , m \in { \cal n}]$ ] .
the rate function is then @xmath30}= \log \left(1 + \frac{|h_{[n , m]}|^2p_{[n , m]}}{\sum_{[n',m ' ] \neq [ n , m]}(|h_{[n',m]}|^2p_{[n',m ' ] } ) + \sigma^2}\right)\ ] ] where @xmath79 is the noise variance at nodes .
it can be readily derived @xmath80}|^2}{\sigma^2 } , \forall [ n , m ] , \;\forall \boldsymbol{s}\ } \\
\delta_2 = \max\{\frac{\sum_{n ' \neq n}\sum_{m ' \in { \cal n}_{n'}^o } |h_{[n , m']}|^2   } { \sigma^2 } , \;\forall \boldsymbol{s}\ } \end{cases}\ ] ] clearly , @xmath81 and @xmath82 are finite as long as all channel gains @xmath83}|^2 $ ] are bounded at all time .
this is in fact also the necessary condition for @xmath29 .
_ example 2 ( orthogonal channels ) _ : if the wireless links do not interfere with each other ( e.g. , when the adjacent nodes operate on orthogonal frequency bands ) , then the rate function is simply @xmath30}= \log \left(1 + \frac { |h_{[n , m]}|^2p_{[n , m ] } } { \sigma^2}\right).\ ] ] we readily have @xmath80}|^2}{\sigma^2 } , \forall [ n , m ] , \;\forall \boldsymbol{s}\ } \\      \delta_2 = 0 \end{cases}\ ] ] again , @xmath81 is finite if all channel gains @xmath83}|^2 $ ] are bounded at all time .
we assume the following two conditions for the system parameters in development of the proposed algorithm : @xmath84    condition ( [ eq.a1 ] ) is a necessary condition to maintain the stability of the energy queues @xmath46 for every sample path .
if @xmath85 , i.e. , the maximum energy arrival is deterministically greater than the largest energy departure possible , then there exists a sample path of energy queue @xmath46 that grows unbounded .
the condition ( [ eq.a1 ] ) is thus required for establishing the sample path result in the ensuing feasibility analysis . on the other hand , condition ( [ eq.a2 ] ) dictates that the battery capacity is large enough to accommodate the largest possible charging / discharging range .
this then makes the system `` controllable '' by the proposed online scheme . ) and ( [ eq.a2 ] ) , we require @xmath86 .
this could be ensured by `` shaping '' or `` clipping '' the energy arrivals by some external mechanisms . ]
our algorithm depends on two algorithmic parameters , namely a `` queue perturbation '' parameter @xmath87 and a weight parameter @xmath88 .
derived from the feasibility requirement of the proposed algorithm ( see the proof of proposition 1 in the sequel ) , any pair @xmath89 that satisfies the following conditions can be used : @xmath90 where    & v^ : = ; [ eq.gv3 ] + & ^ : = + _ 1 g _ v ; [ eq.gv1 ] + & ^ : = - _ 2 g _
v. [ eq.gv2 ]    note that the interval for @xmath88 in ( [ eq.gv ] ) is well - defined under the condition ( [ eq.a2 ] ) , and the interval for @xmath87 is valid when @xmath91 .
we now present the proposed algorithm :    * initialization * : select a pair of @xmath89 satisfying ( [ eq.gv ] ) , and a constant @xmath92 .    at every time slot @xmath8 , observe states @xmath63 , and queues @xmath93 , then determine @xmath94 $ ] , @xmath95}^*(t ) , \forall [ n , m]]$ ] , and @xmath96}^{c*}(t),\forall [ n , m],\forall c]$ ] as follows .
some comments are in order .    * different from the esa algorithm in @xcite
, there is no eh admission mechanism in the proposed algorithm ; the available harvested energy could be then fully capitalized on for data transmission .
* the perturbed energy queue - size @xmath107 is weighted by @xmath108 in the problem ( [ eq.pa ] ) to determine the optimal power allocation .
these weighted are used to account for the battery degeneration and discharging loss . *
solving the power allocation problem ( [ eq.pa ] ) requires centralized control and can be np - hard . for the general ad - hoc networks with interference channels ,
a maximal - weighted - matching based polynomial - time solver @xcite could be used to approximate the optimal scheduling within a constant factor in a distributed manner .
the proposed algorithm is an online scheme , which dynamically makes _ instantaneous _ greedy control decisions for the stochastic system under consideration , without a - priori knowledge of any statistics of the underlying random processes .
we next show that the proposed scheme can also yield a feasible and asymptotically near - optimal solution for the problem ( [ eq.prob ] ) of interest under conditions .
note that in the proposed algorithm , energy availability constraints ( [ eq.ea ] ) and the bounded energy queue constraints ( [ eq.emax ] ) are ignored .
it is then not clear whether the algorithm is feasible for the problem ( [ eq.prob ] ) .
yet , we will show that by using any pair @xmath89 in ( [ eq.gv ] ) and @xmath109 , it is guaranteed that the online control policy produced by the proposed algorithm is a feasible one for ( [ eq.prob ] ) under the conditions ( [ eq.a1])([eq.a2 ] ) .    to this end , we first show the following lemma .    the power allocation specified by the proposed algorithm obeys : i ) @xmath110}^*(t)=0 $ ] , if @xmath111 ; and ii ) @xmath110}^*(t)=p_{\max}$ ] , if @xmath112 , @xmath52 .
see appendix a.    lemma 1 reveals partial characteristics of the proposed dynamic policy .
specifically , when the energy queue at node @xmath10 is large enough , peak power can be afforded for its data transmissions ; i.e. , @xmath110}^*(t)=p_{\max}$ ] .
on the other hand , when the energy queue at node @xmath10 is small enough , no power should be allocated ; i.e. , @xmath110}^*(t)=0 $ ] . based on the structure in lemma 1
, we can then establish the following result .
under the conditions ( [ eq.a1])([eq.a2 ] ) , the proposed algorithm guarantees : i ) @xmath110}^*(t)=0 $ ] , if @xmath113 , and ii ) @xmath114 , @xmath52 , @xmath39 .
see appendix b.    proposition 1 establishes that node @xmath10 allocates nonzero power over any of its outgoing links only when its energy queue - size @xmath115 in the proposed scheme .
hence , all the power allocation decisions are feasible , i.e. , the energy availability constraint ( [ eq.ea ] ) is indeed _
redundant _ and can be ignored in the problem ( [ eq.pa ] ) .
furthermore , the energy queue - sizes @xmath46 resulted from the proposed scheme are guaranteed to be bounded within @xmath116 $ ] .
these two results together imply that the proposed algorithm with proper selection of @xmath89 and @xmath117 can always yield a feasible control policy for ( [ eq.prob ] ) under the conditions ( [ eq.a1])([eq.a2 ] ) .
note that proposition 1 is a _ sample path _ result ; i.e. , it holds for every time slots under _ arbitrary _ , even non - stationary , @xmath63 processes .      to facilitate the optimality analysis , we assume for now that the random process for @xmath63 is i.i.d . over time slots .
define @xmath118 , and @xmath119}(t)\ } \notag\end{aligned}\ ] ] since @xmath46 , @xmath110}(t)$ ] , and @xmath13 are all bounded , it follows from ( [ eq.e ] ) that @xmath120 as @xmath121 $ ] , @xmath39 , ( [ relax1 ] ) then implies @xmath122    for the queue dynamics in ( [ eq.q ] ) , we can derive similar time - average constraints @xcite : @xmath123}^c \leq \sum_{m \in { \cal n}_n^o } \bar{\mu}_{[n , m]}^c   \quad \forall n , c\ ] ] where @xmath124}^c:=\lim_{t\rightarrow \infty } \frac{1}{t } \sum_{t=0}^{t-1 } \mathbb{e}\{\mu_{[n , m]}^c(t)\}.\ ] ]    consider the following problem : @xmath125 note that the queue dynamic constraints ( [ eq.e ] ) and ( [ eq.q ] ) , which need to be performed per realization per slot @xmath8 , are replaced by the corresponding time - average constraints ( [ relax2 ] ) and ( [ relax3 ] ) , and the constraints ( [ eq.ea ] ) and ( [ eq.emax ] ) are ignored
. it can be shown that the problem ( [ eq.relax ] ) is a relaxed version of ( [ eq.prob ] ) . specifically , any feasible solution of ( [ eq.prob ] ) , satisfying ( [ eq.q])([eq.emax ] ) , @xmath39 , also satisfies ( [ relax2 ] ) and ( [ relax3 ] ) , due to the relevant optimization variables .
it then follows that @xmath126 .    in ( [ eq.relax ] )
, the optimization variables are `` decoupled '' across time slots due to the removal of time - coupling constraints ( [ eq.q])([eq.emax ] ) .
this problem has an easy - to - characterize stationary optimal control policy as formally stated in the next lemma .    if @xmath63 is i.i.d .
, there exists a stationary control policy @xmath127 that is a pure ( possibly randomized ) function of the current @xmath63 , while satisfying ( [ eq.rn])([eq.muc ] ) , and providing the following guarantees per @xmath8 : @xmath128}^{stat}(t)\}}{\xi } \leq ( 1-\eta)e_{\max } , \ ; \forall n \\      & \bar{r}_n^{c , stat } + \sum_{m \in { \cal n}_n^i } \mathbb{e}\{\mu_{[m , n]}^{c , stat}(t)\ } \leq \sum_{m \in { \cal n}_n^o}\mathbb{e}\{\mu_{[n , m]}^{c , stat}(t)\ } ,    \ ; \forall n , c \end{split}\ ] ] where @xmath129 , @xmath23}^{stat}(t)$ ] , and @xmath31}^{c , stat}(t)$ ] denote the data admission , power and rate allocation by policy @xmath127 , and expectations are taken over the randomization of @xmath63 and ( possibly ) @xmath127 .
the proof argument is similar to the proof of ( * ? ? ?
* theorem 4.5 ) ; hence , it is omitted for brevity .
lemma 2 in fact holds for many non - i.i.d .
scenarios as well .
generalizations to other stationary processes , or even to non - stationary processes , can be found in @xcite .
it is worth noting that ( [ eq.opt ] ) not only assures that the stationary control policy @xmath127 achieves the optimal utility for ( [ eq.relax ] ) , but also guarantees that the admitted data rate @xmath129 per slot @xmath8 is equal to the optimal time - average rate @xmath130 ( due to the stationarity of @xmath63 and @xmath127 ) .
this will be important to establish the following result :    suppose that conditions ( [ eq.a1])([eq.gv ] ) hold , and @xmath131 is i.i.d . over slots .
let @xmath132 be the time - average admitted rate vector achieved by the proposed algorithm up to time @xmath133 , i.e. , @xmath134 .
then @xmath135 where the constant @xmath136 with @xmath137 ^ 2 , [ -\xi e_{\max}+(1-\eta)\gamma]^2\ } , \notag \\
b_3 & = \eta(1-\eta)\max\{(e_{\max}-\gamma)^2 , \gamma^2\ } ,   \notag\end{aligned}\ ] ] and @xmath138 is the optimal value of ( [ eq.prob ] ) under any feasible control algorithm , even if that relies on knowing future random realizations .    see appendix c.    proposition 2 asserts that the proposed algorithm asymptotically yields a time - average utility with an optimality gap smaller than @xmath139 .
the proposed scheme is in fact a modified version of the queue - length based stochastic optimization scheme ( e.g. , the esa in @xcite ) , where the `` perturbed '' queue lengths play the role of `` stochastic '' lagrange multipliers with a dual - subgradient solver to the problem of interest .
the gap @xmath140 is inherited from the underlying stochastic subgradient method . on the other hand ,
the gap @xmath141 is due to the combined effect of energy - queue perturbation and battery imperfections , while the gap @xmath142 is incurred by the battery degeneration .
based on propositions 1 and 2 , it is now possible to arrive at our main result .    [ them : main ] suppose that conditions ( [ eq.a1])([eq.gv ] ) hold and @xmath131 is i.i.d . over slots
. then the proposed algorithm yields a feasible dynamic control scheme for  , which has an optimality gap @xmath139 ; i.e. , @xmath143 where @xmath134 and @xmath144 is given by ( [ eq.b ] ) .
interesting comments on the minimum optimality gap with the proposed algorithm are in order .
the results in theorem 1 can be generalized to the more general non - i.i.d .
case where the the energy state @xmath15 and channel state @xmath18 both evolve according to some finite - state irreducible , aperiodic markov chains .
note that the feasibility of the proposed algorithm is a sample path result which does not depend on the assumption on the random processes per remark 2 ; i.e. , it holds also for this non - i.i.d . case . on the other hand ,
the performance guarantee ( i.e. , optimality gap result ) can be obtained by applying the so - called delayed lyapunov drift technique , i.e. , the method of analyzing regenerative cycles of the markov random processes ; see ( * ? ? ?
* theorem 3 ) and ( * ? ? ? * theorem 2 ) .
in this section , numerical results are provided to evaluate the proposed algorithm .
a simple network is considered in fig .
1 , where the nodes 1 - 4 collect data and send data to the sink node 7 through relay nodes 5 and 6 @xcite .    in simulations
, we assume imperfect batteries at nodes , with storage efficiency @xmath164 , and ( dis-)charging efficiency @xmath165 .
table [ tab : param ] lists the values for @xmath166 ( the maximum in- and out - degree for nodes in the network ) , @xmath167 ( the maximum packets that can be newly admitted ) , @xmath168 ( the peak power ) , @xmath169 ( the maximum rate over all the links ) and @xmath170 ( the battery capacity ) .
the utility function is selected as : @xmath171 .
suppose that all the links are independent with each other , implying @xmath172 .
the link state can be either good or bad with equal probability .
one unit of power can deliver two packets when the link state is good , while it can be only used to transmit one packet upon bad link state .
we also assume that the harvested energy @xmath13 is i.i.d . for each node ; @xmath13 is either @xmath173 or 0 with equal probability . as a result
, we have : @xmath174 , @xmath175 and @xmath176 . to satisfy the constraint ( [ eq.a1 ] ) , we set @xmath177 , unless otherwise specified .
each simulation is run for 1200 slots and each result is obtained by averaging over 10 independent runs .
2 depicts the achieved utility for the proposed algorithm with different @xmath87 , when the weight parameter @xmath178 , @xmath179 and @xmath180 .
it is shown that the maximum utility is achieved when the  queue perturbation " parameter @xmath87 is set to @xmath181 in this case .
this is because power allocation policy ( 17 ) tends to save energy if a large @xmath87 is chosen , leading to a large data queue size @xmath182 . on the other hand ,
data admission policy ( [ eq.da ] ) dictates that @xmath183 , if @xmath184 , @xmath185 .
this implies that large @xmath182 reduces rate utility . as a result ,
the utility decreases as @xmath87 grows .
for this reason , we set @xmath186 unless otherwise specified in all the subsequent simulations .
note that when @xmath187 and @xmath188 , we often have @xmath23}(t)=0 , \forall n , \forall t$ ] ; i.e. , no data are allowed to be transmitted over links . in this case , to maintain the queue stability , network utility becomes almost 0 .
[ fig : weight ] shows the network utility of the proposed algorithm with different weight parameter @xmath88 , when @xmath179 and @xmath180 .
it is observed that the network utility decreases after reaching the maximum value for a given @xmath189 .
clearly , the proposed algorithm does not perform well with a small @xmath88 . on the other hand ,
a large @xmath88 is not a good choice in the presence of battery imperfections . to see it ,
we use @xmath181 in ( 14 ) and set @xmath31}(t ) = s_{[n , m]}(t ) p_{[n , m]}(t)$ ] in ( 17 ) .
( @xmath190}(t)$ ] is 2 with good link state , while it is 1 when the link state is bad . )
then power allocation problem ( 17 ) is a linear program , where the weight for @xmath23}(t)$ ] is given by @xmath101}(t)s_{[n , m]}(t)+ \eta e_n(t)-p_{\max}-\delta_1 g_{\max } v$ ] .
clearly , we have @xmath23}(t)=0 $ ] , if its weight is less than zero .
therefore , a large @xmath88 implies that high energy is required to render @xmath23}(t)>0 $ ] for transmission .
however , it is hard to maintain high energy in the presence of energy degeneration . as a result , when @xmath187 and @xmath191 , utility is almost 0 .
the results are consist with our discussion in remark 4 .
the optimality gap with the proposed scheme does not monotonically decrease with @xmath88 any more in the presence of battery imperfections .
determination of the best @xmath88 for the algorithm may refer to the solution of the problem formulated in remark 4-ii ) .    .parameters configuration . [ cols="^,^,^,^,^",options="header " , ]         .,width=336,height=192 ]    .,width=336,height=192 ]        the data and energy queue processes with the proposed algorithm are plotted in fig . [
fig : samplepro ] , when @xmath179 , @xmath192 and @xmath178 .
clearly , the three data queue lines evolve within the region [ 0 , @xmath193 ; i.e. , the stability of the data network is maintained .
in addition , all energy queues satisfy @xmath194 .
this clearly corroborates the feasibility of the proposed scheme , as stated in proposition 1 .
two baseline schemes are chosen to gauge the performance of the proposed algorithm .
the first baseline scheme is esa in @xcite , where perfect batteries are assumed and an eh admission is performed to ensure finite energy queues .
the second scheme is a heuristic greedy algorithm , where a node with the largest data queue backlog is first allowed to transmit data over its best outgoing links ( with best link quality ) .
the node with second largest data queue then selects a best link to transmit , among all its outgoing links that do not conflict with the existing links .
( a link conflicts with the existing links if it shares the same transmit or receive node with any existing links . )
this continues until no nodes or no links can be selected .
if a link @xmath24 $ ] is selected , the power @xmath23}(t)$ ] is set as the maximum value that satisfies both the peak power constraint ( 2 ) and the energy availability constraint ( 6 ) ; @xmath31}(t)$ ] is then determined correspondingly .
finally , the greedy algorithm sets @xmath195 if @xmath196 , while @xmath197}(t)$ ] if @xmath198 .
it is clear that the nodes have finite data queue backlogs ; i.e. , the system is stable .
the feasibility of the greedy algorithm is also guaranteed since the energy availability constraint ( 6 ) is taken into account in power allocation policy .
[ fig : results ] compares the performance of the three algorithms with different @xmath173 , when @xmath199 , @xmath192 and @xmath178 .
it is shown that the proposed algorithm evidently outperforms the esa and the greedy algorithm for any given @xmath173 .
since the greedy algorithm in fact schedules the links in a time division multi - access ( tdma ) manner , the resultant utility is low in general .
in addition , the required energy for implementation of the greedy policy is actually also low ; hence , its performance changes only slightly for different @xmath173 .
the proposed algorithm achieves higher utility than esa due to two reasons .
the first reason is that esa can not make use of all available energy because of its eh admission mechanism , while the proposed algorithm harvested all available energy .
the utilization of available energy ( the ratio of actually harvested energy to available energy ) for esa is depicted in fig .
[ fig : uti ] . for esa ,
node @xmath10 stops harvesting energy when @xmath46 energy reaches @xmath200 @xcite . clearly , given a larger @xmath173 ( i.e.
, renewable energy is abundant ) , esa might waste more available energy , leading to a significant performance loss . on the other hand
, esa also has a performance loss for small @xmath173 case since it does not take into account the battery imperfections .
for instance , when @xmath201 , the utilization of available energy for esa is 100% ; yet , the utility of the proposed algorithm is still 17.2% larger than that of esa in this case .
in this paper , a dynamic resource allocation task was considered for general eh wireless networks . taking into account imperfect finite - capacity energy storage devices ,
a stochastic optimization was formulated to maximize the long - term utility subject to the energy availability constraints .
capitalizing on lyapunov optimization technique , an online control algorithm was proposed to make data admission , power allocation and routing decisions , without requiring any statistical knowledge of channel , data - traffic , and eh processes .
it was shown that the proposed algorithm can efficiently utilize the harvested energy to provide a feasible and asymptotically near - optimal control solution for general eh networks .
using the perturbation weight @xmath92 in queue - backpressure , it was shown in ( * ? ? ? * theorem 2 ) that the data queue sizes satisfy : @xmath202 , @xmath98 , @xmath39 .
note that this is a sample path result which holds for arbitrary @xmath15 and @xmath18 processes .      *
if @xmath205 , suppose @xmath110}^*(t ) > 0 $ ] ; i.e. , there exists a certain @xmath206}^*(t)>0 $ ] with @xmath207 .
we can create a new power vector @xmath66 by setting only @xmath206}'=0 $ ] and all other components remain the same as in @xmath102 .
let @xmath208 denote the objective function of ( [ eq.pa ] ) .
it readily follows that @xmath209}(t)[\mu_{[n , m]}(\boldsymbol{p}^*(t ) ) - \mu_{[n , m]}(\boldsymbol{p } ' ) ] \notag\\          &   ~~~~~ + \frac{\eta}{\xi } ( e_n(t)-\gamma ) p_{[n,\check{m}]}^*(t ) \notag \\          &   \leq w_{[n,\check{m}]}(t)[\mu_{[n,\check{m}]}(\boldsymbol{p}^*(t ) ) - \mu_{[n,\check{m}]}(\boldsymbol{p } ' ) ] \notag \\          &   ~~~~~ + \frac{\eta}{\xi } ( e_n(t)-\gamma ) p_{[n,\check{m}]}^*(t ) \notag \\          &   < [ g_{\max}v -d_{\max}\mu_{\max}]^+ \delta_1
p_{[n,\check{m}]}^*(t ) - \delta_1 g_{\max } v p_{[n,\check{m}]}^*(t ) \notag \\          &   < 0 \notag      \end{aligned}\ ] ] where the first inequality is due to property 1-i ) , and the second inequality is due to property 1-ii ) , as well as the bounds @xmath101}(t ) \leq [ g_{\max}v -d_{\max}\mu_{\max}]^+$ ] and @xmath205 .
this shows that @xmath102 can not have been the optimal power allocation for ( [ eq.pa ] ) , leading to a contradiction .
hence , part i ) in the lemma must hold . * if @xmath210 , suppose @xmath110}^*(t)<p_{\max}$ ] .
let @xmath211}^*(t)$ ] .
we can create a new power vector @xmath66 by setting @xmath23}'=p_{[n , m]}^*(t)+\frac{\delta p}{|{\cal n}_n^o|}$ ] , @xmath71 , and @xmath72}'=p_{[n',m]}^*(t)$ ] , @xmath73 .
it readily follows that @xmath209}(t)[\mu_{[n , m]}(\boldsymbol{p}^*(t ) ) - \mu_{[n , m]}(\boldsymbol{p } ' ) ] \notag\\          &   ~~~~~ - \frac{\eta}{\xi } ( e_n(t)-\gamma ) \delta p \notag \\          &   \leq \sum_{n ' \neq n } \sum_{m \in { \cal n}_{n'}^o } w_{[n',m]}(t)[\mu_{[n',m]}(\boldsymbol{p}^*(t ) ) - \mu_{[n',m]}(\boldsymbol{p } ' ) ] \notag \\          &   ~~~~~ - \frac{\eta}{\xi } ( e_n(t)-\gamma ) \delta p \notag \\          &   < [ g_{\max}v -d_{\max}\mu_{\max}]^+ \delta_2 \delta p - \delta_2 g_{\max } v \delta p \notag \\          &   < 0 \notag      \end{aligned}\ ] ] where the first inequality is due to property 2-i ) , and the second inequality is due to property 2-ii ) , as well as the bounds @xmath101}(t ) \leq [ g_{\max}v -d_{\max}\mu_{\max}]^+$ ] and @xmath210 .
this leads to a contradiction .
hence , part ii ) in the lemma must also hold .
the proof for part ii ) proceeds by induction .
first , recall that @xmath215 , @xmath52 , and suppose that it holds @xmath121 $ ] , @xmath52 , at slot @xmath8 .
we will show the bounds hold for @xmath216 , @xmath52 , as well in subsequent instances .      * if @xmath219 , then it follows from lemma 1 that @xmath110}^*(t)= 0 $ ] . from ( [ eq.e ] )
, we have * * @xmath220 ; * * @xmath221 , due to @xmath222 , @xmath223 , and the definition of @xmath217 in ( [ eq.gv2 ] ) . * if @xmath224 $ ] , then @xmath225}^*(t ) \leq p_{\max}$ ] .
we therefore have * * @xmath226 , due to @xmath227 and the definition of @xmath181 in ( [ eq.gv1 ] ) ; * * @xmath228 , due to @xmath223 as with c1-ii ) . *
if @xmath229 $ ] , it follows from lemma 1 that @xmath110}^*(t ) = p_{\max}$ ] .
we have * * @xmath230 , due to @xmath227 as with c2-i ) ; * * @xmath231 , due to condition ( [ eq.a1 ] ) .
consider the queues under the proposed algorithm . by squaring both sides of ( [ eq.q ] ) , and using the fact that @xmath232^+)^2 \leq x^2 $ ] , we have @xmath233 ^ 2-[q_n^c(t)]^2 \notag \\      & \leq \left[\sum_{m \in { \cal n}_n^o } \mu_{[n , m]}^{c*}(t ) \right]^2 + \left[\sum_{m \in { \cal n}_n^i } \mu_{[m , n]}^{c*}(t ) + r_n^{c*}(t ) \right]^2 \notag \\      & ~~~~ -2q_n^c(t)\left
[ \sum_{m \in { \cal n}_n^o } \mu_{[n , m]}^{c*}(t ) - \sum_{m \in { \cal n}_n^i } \mu_{[m , n]}^{c*}(t ) - r_n^{c*}(t)\right ] . \notag\end{aligned}\ ] ] recall that @xmath234}^{c*}(t ) \leq d_{\max } \mu_{\max}$ ] and @xmath235}^{c*}(t ) + r_n^{c*}(t ) \leq d_{\max } \mu_{\max}+r_{\max}$ ] . upon defining @xmath236
, we readily have : @xmath237 ^ 2-[q_n^c(t)]^2 ) \notag \\      & \leq \tilde{b}_1 - q_n^c(t)\left [ \sum_{m \in { \cal n}_n^o } \mu_{[n , m]}^{c*}(t ) - \sum_{m \in { \cal n}_n^i } \mu_{[m , n]}^{c*}(t ) - r_n^{c*}(t)\right ] .
\label{eq.dq}\end{aligned}\ ] ]    similarly , we can also derive @xmath238 ^ 2-[e_n(t)-\gamma]^2 ) \notag \\      & \leq -\frac{1}{2}(1-\eta^2 ) [ e_n(t)-\gamma]^2 \notag \\      & ~~~~ + \frac{1}{2 } \left[\frac{\sum_{m \in { \cal n}_n^o } p_{[n , m]}^*(t)}{\xi } - \xi e_n(t ) + ( 1-\eta)\gamma\right]^2 \notag\\      & ~~~~ -\eta(e_n(t)-\gamma)\left[\frac{\sum_{m \in { \cal n}_n^o } p_{[n , m]}^*(t)}{\xi } - \xi e_n(t ) + ( 1-\eta)\gamma\right ] \notag \\      & \leq b_2 - \frac{\eta}{\xi}(e_n(t)-\gamma ) \sum_{m \in { \cal n}_n^o } p_{[n , m]}^*(t ) \notag \\      & ~~~~ + \eta ( e_n(t)-\gamma)[\xi e_n(t ) - ( 1-\eta)\gamma ] \label{eq.de}\end{aligned}\ ] ] where @xmath239 ^ 2 , [ -\xi e_{\max}+(1-\eta)\gamma]^2\}$ ] .    consider the lyapunov function @xmath240 ^ 2 + \frac{1}{2 } \sum_n [ e_n(t)-\gamma]^2 $ ] . by summing ( [ eq.dq ] ) over all @xmath241 and ( [ eq.de ] ) over all @xmath10 , and by defining @xmath242
, we have @xmath243}^{c*}(t ) - \sum_{m \in { \cal n}_n^i } \mu_{[m , n]}^{c*}(t ) - r_n^{c*}(t)\right ] \notag \\      & ~~~~ - \sum_n \frac{\eta}{\xi}(e_n(t)-\gamma ) \sum_{m \in { \cal n}_n^o } p_{[n , m]}^*(t )   \notag \\      & ~~~~ + \sum_n \eta ( e_n(t)-\gamma)[\xi e_n(t ) - ( 1-\eta)\gamma]\end{aligned}\ ] ] furthermore , we have @xmath244}^{c*}(t ) \notag \\      & \leq \tilde{b } - \sum_{n , c } [ v u_n^c(r_n^{c*}(t ) ) - q_n^c(t )   r_n^{c*}(t ) ] \notag \\      & ~~~~ - \sum_n \bigl[\sum_c \sum_{m \in { \cal n}_n^o } \mu_{[n , m]}^{c*}(t ) [ q_n^c(t ) -   q_m^c(t ) - \theta]\bigr .
\notag \\      & ~~~~~~~~~~~~~~ \bigl .
+ \frac{\eta}{\xi}(e_n(t)-\gamma ) \sum_{m \in { \cal n}_n^o } p_{[n , m]}^*(t ) \bigr ]   \notag \\      & ~~~~ + \sum_n \eta ( e_n(t)-\gamma)[\xi e_n(t ) - ( 1-\eta)\gamma ] \notag \\      & \leq   \tilde{b } - \sum_{n , c } [ v u_n^c(r_n^{c , stat}(t ) ) - q_n^c(t )   r_n^{c , stat}(t ) ] \notag \\      & ~~~~ - \sum_n \bigl[\sum_c \sum_{m \in { \cal n}_n^o } \mu_{[n , m]}^{c , stat}(t ) [ q_n^c(t ) -   q_m^c(t ) - \theta]\bigr .
\notag \\      & ~~~~~~~~~~~~~~ \bigl .
+ \frac{\eta}{\xi}(e_n(t)-\gamma ) \sum_{m \in { \cal n}_n^o } p_{[n , m]}^{stat}(t ) \bigr ]   \notag \\      & ~~~~ + \sum_n \eta ( e_n(t)-\gamma)[\xi e_n(t ) - ( 1-\eta)\gamma ] \label{eq.rhs}\end{aligned}\ ] ] where the second inequality is due to the fact that the proposed algorithm in fact minimizes the right - hand - side of ( [ eq.rhs ] ) .    taking expectations for ( [ eq.rhs ] )
, we arrive at : @xmath245}^{c*}(t)\bigr\ } \notag \\      & \leq   \tilde{b } - v \sum_{n , c } u_n^c(\bar{r}_n^{c , stat}(t ) ) \notag \\      & ~~ + q_n^c(t ) \bigl [ \bar{r}_n^{c , stat } + \sum_{m \in { \cal n}_n^i } \mathbb{e}\{\mu_{[m , n]}^{c , stat}(t)\ } - \sum_{m \in { \cal n}_n^o}\mathbb{e}\{\mu_{[n , m]}^{c , stat}(t)\ } \big ] \notag \\      & ~~ + \sum_n \eta ( e_n(t)-\gamma ) \bigl[\xi \bar{e}_n - \frac{\sum_{m \in { \cal n}_n^o } p_{[n , m]}^{stat}(t)}{\xi } - ( 1-\eta)\gamma \bigr ]   \notag \\      & ~~ + \theta \sum_n \sum_c \sum_{m \in { \cal n}_n^o } \mathbb{e}\{\mu_{[n , m]}^{c , stat}(t)\ } \notag \\      & \leq \tilde{b }   - v \tilde{u}^{opt } + n \eta(1-\eta)\max\{(e_{\max}-\gamma)^2 , \gamma^2\ } \notag \\      & ~~ + n^2 \theta d_{\max}\mu_{\max}\end{aligned}\ ] ] where the second inequality is due to ( [ eq.opt ] ) in lemma 2 , @xmath121 $ ] by proposition 1 , and @xmath246}^{c , stat}(t)\ } \leq d_{\max}\mu_{\max}$ ] .    with @xmath247 and noticing @xmath248}^{c*}(t)\ } \geq 0 $ ] , we then have @xmath249 summing over all @xmath250 , we have @xmath251 which leads to @xmath252 using jensen s inequality , we see that @xmath253 the proposition follows by taking the limit @xmath254 and using the definition of @xmath132 .
o. ozel , k. tutuncuoglu , y. jing , s. ulukus , and a. yener , `` transmission with energy harvesting nodes in fading wireless channels : optimal policies , '' _ ieee j. sel .
areas commun .
8 , pp . 17321743 , sep . 2011 .
o. ozel , y. jing , s. ulukus , `` optimal broadcast scheduling for an energy harvesting rechargeable transmitter with a finite capacity battery , '' _ ieee trans .
wireless commun .
11 , no . 6 , pp
. 21932203 , jun . 2012 .
o. ozel , k. shahzad , and s. ulukus , `` optimal energy allocation for energy harvesting transmitters with hybrid energy storage and processing cost , '' _ ieee trans .
signal process .
12 , pp . 32323245 , 2014 .
x. wang , z. nan , and t. chen , `` optimal mimo broadcasting for energy harvesting transmitter with non - ideal circuit power consumption , '' _ ieee trans .
wireless commun .
5 , pp . 25002512 , may 2015 .
b. devillers and d. gunduz , `` a general framework for the optimization of energy harvesting communication systems with battery imperfections , '' _ j. commun
2 , pp . 130139 , apr .
2012 .
r. srivastava and c. e. koksal , `` basic performance limits and tradeoffs in energy - harvesting sensor nodes with finite data and energy storage , '' _ ieee / acm trans .
4 , pp . 10491062 , 2013 .
l. lin , n. shroff , and r. srikant , `` asymptotically optimal energy - aware routing for multihop wireless networks with renewable energy sources , '' _ ieee / acm trans .
5 , pp . 10211034 , 2007 .
s. chen , p. sinha , n. shroff , and c. joo , `` a simple asymptocially optimal joint energy allocation and routing schemes in rechargable sensor networks , '' _ ieee / acm trans .
netw . _ , vol .
4 , pp . 13251336 , aug . 2014 .
z. mao , c. koksal , and n. shroff , `` near optimal power and rate control of multi - hop sensor networks with energy replenishment : basic limitations with finite energy and data storage , '' _ ieee trans .
4 , pp . 815829 , 2012 .
j. qin , y. chow , j. yang , and r. rajagopal , `` online modified greedy algorithm for storage control under uncertainty , '' _ ieee trans .
power syst . _ ,
2015 ( accepted ) , [ online ] .
available : http://arxiv.org/pdf/1405.7789v4.pdf .


it is generally agreed that low - frequency ( radio ) emission from relativistic jets in active galactic nuclei ( agn ) is synchrotron emission from nonthermal electrons . however , it is still debated whether the high - frequency ( x - ray & @xmath0-ray ) emission is due solely to compton scattering by the electrons responsible for the low - frequency synchrotron emission or due to hadronic processes resulting from protons being co - accelerated with the electrons ( e.g. , * ? ? ?
if hadronic models are valid , agn would be an expected source of neutrinos , as the proton interactions responsible for the production of neutral pions which result in @xmath0-rays also produce charged pions , which result in secondary electrons , positrons , and neutrinos ( e.g. , * ? ? ?
can we expect agn to be detectable neutrino sources at the level that starburst systems of similar @xmath0-ray luminosity are expected to be ?
observations by the space - based @xmath0-ray telescope _ fermi _ and future observations by the hard x - ray space telescope nustar may help to further resolve this question .    building on the approach of previous models ( e.g. *
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , we developed a semi - analytic model for cosmic ray interactions for the starburst galaxy m82 ( * ? ? ?
* hereafter yegz ) . in developing the model
, we emphasized the important role of galactic winds while limiting the number of free parameters and keeping our models observationally informed .
now , we aim to determine whether the model is applicable to other starburst systems and explore benefits and limitations of such models due to the presence of agn .
the starburst galaxy ngc 253 has been detected in @xmath0-rays by both space - based ( _ fermi _ ) and ground - based ( hess ) telescopes and possibly contains an active central supermassive black hole ( smbh ) @xcite .
seyfert 2 galaxies ngc 1068 and ngc 4945 have also been detected with _
fermi_. though ngc 1068 ( @xmath1 mpc ;
nasa / ipac extragalactic database ) is more distant than the nearby ngc 253 ( @xmath2 mpc ; * ? ? ?
* ) and ngc 4945 ( @xmath3 mpc ; * ? ? ?
* ) , all three have @xmath0-ray fluxes of the same order of magnitude @xcite .    in this paper
, we determine whether our model is applicable to a combined starburst / agn environment .
we select the starburst galaxy ngc 253 with a possible low - activity agn and the high - activity seyfert galaxy ngc 1068 . applying our model to ngc 253 and ngc 1068
, we seek the best - fit cosmic ray interaction model which accurately reproduces both radio and @xmath0-ray spectra within the confines of the observed properties of the interstellar medium ( ism ) and surrounding radiation field .
additionally , we make a prediction of the neutrino fluxes from both galaxies .
llclc distance & 3.9 mpc & 1 & 14.4 mpc & 7 + central molecular zone ( cmz ) radius & 150 pc & 2 & 200 pc & 8 + molecular gas mass & @xmath4 @xmath5 & 3 , 4 & @xmath6 @xmath5 & 9 + average ism density & 500 - 1400 @xmath7 & & 250 @xmath7 & + ir luminosity & @xmath8 @xmath9 & 5 & @xmath10 @xmath9 & 10 + radiation field energy density & 500 - 2000 ev  @xmath7 & & @xmath11 ev  @xmath7 & + sn explosion rate & 0.1 yr@xmath12 & 6 & 0.07 yr@xmath12 & 11 + sn explosion energy & 10@xmath13 ergs & & 10@xmath13 ergs & + sn energy transferred to cr & 4 - 20% & & 10% & + ratio of primary protons & 50 & & 50 & +  to electrons ( n@xmath14/n@xmath15 ) & & & & + slope of primary cr & 2.2/2.3 & & 2.0 & +  source function & & & & +    the next section describes the observed properties of ngc 253 and ngc 1068 . section 3 details how the population of energetic particles was computed .
section 4 contains the results of our model for ngc 253 and ngc 1068 . in section 5 ,
we compare and contrast the results for ngc 253 and ngc 1068 and discuss predictions for the neutrino flux and the resulting implications . in section 6 , we present concluding remarks .
ngc 253 is a giant , barred spiral galaxy with a central starburst region .
tidal interactions are the cause of many starbursts in massive galaxies , and while there is no obvious companion or trigger for the starburst in ngc 253 , kinematic studies suggest ngc 253 may have been involved in a past merger @xcite .
models indicate that the current high level of nuclear star formation has been ongoing for at least 20 - 30 myr @xcite .
the center of ngc 253 is host to several compact objects which have been observed over a variety of wavelengths @xcite .
the strongest non - thermal , compact radio source ( th2 ) has long been considered a possible agn candidate @xcite .
additionally , x - ray emission from the nucleus may be powered by an obscured agn or a combination of supernova remnants and fluorescent line emission from molecular clouds @xcite . however , new observations with nustar and _ chandra _ are unable to rule out the possibility that these compact x - ray sources are ultraluminous x - ray sources ( ulxs ) and not a low - luminosity agn @xcite
.    one of the other major features of ngc 253 is the galactic wind emanating from the nuclear region ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
this galactic wind has been detected in both h@xmath16 and x - rays .
observations of the position and size of the outflow at both wavelengths match up almost exactly @xcite .
optical observations give an outflow velocity on the order of a few 100 km  s@xmath12 ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
additionally , @xcite suggest that , like star formation , the wind is quasi - steady .    due to the obscured nature of the starburst nucleus of ngc 253 ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , determining the properties of the nuclear region has been difficult .
radio observations yield supernova rates ranging from 0.03 - 0.3 yr@xmath12 @xcite . in this work ,
we adopt a rate of 0.1 yr@xmath12 based on @xcite ( see table 1 for additional parameters ) .
additionally , estimates of the molecular gas mass range from @xmath17 @xmath5 ( * ? ? ?
* and references therein ) .
we adopt a value of @xmath18 @xmath5 for the molecular gas mass for this paper , based on work by @xcite .
we adopted an ionized gas mass for ngc 253 by scaling from the ratio of ionized to molecular gas in m82 , which gives a mass of @xmath19 @xmath5 .
ngc 253 is well observed in the radio for both extended emission and the central starburst core .
observations by @xcite and @xcite at 333 mhz and from 1 to 7 ghz show that the radio spectrum turns down only slightly in the core at low frequencies , unlike the radio spectrum for the core of m82 , which is completely absorbed at low frequencies .
high frequency observations from @xcite at 18.5 and 22 ghz and from @xcite at 94 ghz require a free - free emission contribution to the flux in addition to non - thermal synchrotron emission . observed radio and @xmath0-ray spectra are presented with model spectra in section 4 .
ngc 253 has been observed in @xmath0-rays at both gev and tev energies by _ fermi _ and hess .
@xcite present four data points and two upper limits over 0.2 to 200 gev from analysis of 30 months of _ fermi _ data .
@xcite present a combined analysis of five data points and one upper limit over 0.2 tev to 30 tev from hess data with the _ fermi _ data and fit a single power - law to both data sets , resulting in a spectral index of @xmath20 .
@xcite present a more current reduction of the _ fermi _ data which we use for comparison with our models instead of the _ fermi _ data shown in @xcite .
ngc 1068 is the closest , brightest seyfert galaxy .
the total bolometric luminosity of the galaxy is @xmath21 @xmath9 ; half of which is accounted for by the circumnuclear disk ( cnd ) and agn @xcite .
the galaxy supports a kiloparsec - scale , non - thermal radio jet @xcite and a compact radio nucleus @xcite .
the complex chemical and kinematic properties of the nuclear region have been well - studied at multiple wavelengths @xcite .
the circumnuclear starburst ring surrounding the 2.3 kpc stellar bar @xcite is illuminated by star formation in the past 10 - 40 myr @xcite .
the stellar population inside the cnd , the inner @xmath22200 pc , developed in two bursts at 30 and 300 myr ago @xcite .
@xcite give the star formation rate ( sfr ) for the inner galaxy as @xmath22 10 @xmath5 yr@xmath12 .
this gives an upper bound on the supernova rate of @xmath220.07 yr@xmath12 for the nuclear region .
we convert from sfr to supernova rate with by analysis similar to @xcite .
both starburst regions are rich in molecular gas with masses of @xmath23 @xmath5 for the circumnuclear ring and @xmath24 @xmath5 for the cnd @xcite .
co , hcn , and hco@xmath25 observations show that the gas is moderately dense ( @xmath26 @xmath7 ) with a kinetic temperature of @xmath27 k @xcite .      as the nearest and brightest seyfert galaxy , ngc 1068 has been extensively studied in the radio . on the kiloparsec - scale , radio observations have mapped the star - forming spiral arms and the radio jets .
there are extensive observations specifically targeting the agn on the parsec - scale @xcite .
however , as even observations of the inner hundred parsecs have a dominant contribution from the agn @xcite , specifically the radio jet and counter - jet , there are essentially no observations of the central molecular zone ( cmz ) alone .
radio maps show little evidence of ongoing star - formation in the central @xmath22150 pc .
@xmath0-ray data for ngc 1068 are more limited than for ngc 253 . while there are four data points and two upper limits at gev energies , there is only a single upper limit for tev energies from hess @xcite .
the _ fermi _ data range from 0.2 gev to 300 gev and the hess upper limit is at 0.2 tev .
diffusion is thought to be of minor importance in star - forming environments with high gas densities and galactic winds . as such ,
following the approach in yegz , a cosmic ray spectrum can be calculated from a single quantity including the source function ( dependent on supernova rate ) and the lifetime . because we only consider radiative and advective losses , our model is not applicable to spatially extended regions , such as the inner spiral arms in ngc 1068 , where cosmic ray diffusion is important .
thus , we apply our model specifically to the cmzs of galaxies .
supernovae are the assumed drivers of cosmic ray acceleration , and so we adopt a power - law source function ( @xmath28 ) such that @xmath29 where @xmath30 is the supernova rate , @xmath31 is the volume of the starburst region , @xmath32 is the fraction of the supernova energy transferred to cosmic rays , and @xmath33 is 10@xmath13 ergs , the typical energy from a supernova explosion .
thus , the spectrum for primary protons is given by @xmath34 where @xmath35 is the minimum cosmic ray energy , here taken to be @xmath36 gev .
@xmath37 is the cosmic ray lifetime composed of an energy loss timescale ( due to radiative and collisional losses ) and an energy independent timescale : @xmath38 where @xmath39 is the scale height of the starburst region and @xmath40 is the speed of the particles in the wind of the starburst region ( see also * ? ? ?
energy losses include ionization , pion production , bremsstrahlung , synchrotron emission , and the inverse compton effect .
pion production from proton - proton collisions is the dominant energy loss for cosmic ray protons above the 1.3 gev threshold .
as the lifetimes of both charged and neutral pions are very short , pions quickly decay into secondary electrons and positrons , neutrinos , and @xmath0-rays .
the primary proton source function can be used to calculate the source function for charged and neutral pions @xcite , and the pion source function can then be used to calculate the spectrum of secondary electrons and positrons @xcite .
one of the main assumptions of the model is that cosmic rays sample the mean density of the interstellar medium .
the density of the interstellar medium sampled by the cosmic rays is vital to the model as it affects pion production and thus both the @xmath0-ray and radio spectra .
the majority of the ism is hot , low density gas in which pockets of cold , very dense molecular gas and warm ionized gas can be found . in @xcite , we assumed a simple two phase model of cold , molecular gas and hot , ionized gas of negligible density .
we found that for a wide range of cosmic ray injection conditions , magnetic field properties , and wind speeds , the cosmic rays sample the mean density of the ism to within a factor of two .
the effect of a third phase of dense , warm , ionized gas on cosmic ray sampling will be considered in future work .    from cosmic ray proton and primary and secondary electron / positron populations
, we can calculate nonthermal radio synchrotron emission and @xmath0-ray emission by neutral pion decay , bremsstrahlung , and inverse compton ( see yegz for further details ) for comparison with observed data .
primary variables include spectral index ( @xmath41 ) , cosmic ray acceleration efficiency ( @xmath32 ) , molecular gas mass ( @xmath42 ) , magnetic field strength ( @xmath43 ) , wind ( advection ) speed ( @xmath40 ) , ionized gas density ( @xmath44 ) , and interstellar radiation field energy density ( @xmath45 ) .      in yegz
, we found that observations of the starburst core in m82 showed a complete turn down of the radio emission at low frequencies .
we attributed this to severe free - free absorption in the core and adopted a model for the starburst core with a cylindrical wall of warm , ionized gas surrounding a lower density hot medium . because some observations at low frequencies did not have enough resolution to distinguish the core from the halo in m82
, we also included a component for radio emission from the halo .
observations for ngc 253 do not show the same severe free - free absorption at low frequencies as in m82 .
however , there is significant free - free emission at high frequencies , as is also the case in m82 .
thus , we expect that at least some fraction of the radio synchrotron emission is still absorbed by ionized gas . to account for this , we keep the thin shell geometry ( with shell width @xmath46 and radius @xmath47 ) for the nucleus , but assume only a small fraction ( @xmath48 ) of the synchrotron emission is absorbed by the ionized gas . then , the radiative intensity of the absorbed synchrotron emission is given by @xmath49 where @xmath50 is the radius of the hot , diffuse gas ( @xmath51 ) , @xmath52 is the emission coefficient for the hot gas ( @xmath53 ) , and @xmath54 is the absorption coefficient for the warm , ionized gas ( @xmath55 ) .
the majority of the radio spectrum is unabsorbed synchrotron emission from the hot , diffuse gas : @xmath56 additionally , we have some portion of the radio spectrum composed of free - free emission from the ionized gas clouds : @xmath57 where @xmath58 is the emission coefficient for the warm , ionized gas ( @xmath59 ) .
thus , the total emergent intensity from the starburst region is @xmath60 for an absorption fraction of @xmath61 , this reduces to @xmath62 which is the expression for the radio emission that we used for modeling m82 ( yegz ) .
note that unlike in yegz , we do not include a disk / halo component for radio emission in this paper as radio observations of the disk are separable from the halo ( e.g. , * ? ? ?
* ; * ? ? ?
we anticipate that inverse compton scattering is a vital process for producing @xmath0-rays in ngc 253 and ngc 1068 .
for inverse compton scattering , the @xmath0-ray source function is given by @xcite @xmath63 with @xcite @xmath64,\ ] ] where @xmath65 is the energy of the resulting @xmath0-ray , @xmath66 is the energy of the incident photon , @xmath0 is the energy of the electron .
here , the combined cosmic ray electron / positron spectrum , @xmath67 , is in units of @xmath7 .
the function @xmath68 is part of the klein - nishina cross section and is given by @xcite @xmath69 where @xmath70 for the blackbody spectrum , @xmath71 , we use an isotropic , diluted , modified blackbody spectrum ( * ? ? ?
* and references therein ) @xmath72 where @xmath73 is a spatial dilution factor ( given by the normalization @xmath74 ) and @xmath75 corresponds to @xmath76 hz .
as with our previous model for m82 ( yegz ) , we test a range of parameters to find the best fits to both the radio and @xmath0-ray spectra for the starburst core of ngc 253 ( see table 2 ) .
we choose to vary parameters such as magnetic field strength ( @xmath77 @xmath78 g ) and radiation field energy density ( @xmath79 ev  @xmath7 ) as they are not well - constrained by observations .
we vary acceleration efficiency ( @xmath80 ) as it is equivalent to varying the supernova rate and test two spectral indices ( @xmath81 ) . as the observed @xmath0-ray spectrum for ngc 253 is steeper than that of m82 ,
we choose larger spectral indices and do not test for @xmath82 .
we also test a variety of wind ( advection ) speeds ( @xmath83 km  s@xmath12 ) in order to study the no - wind case and wind speeds similar to those derived from optical and x - ray observations .
ll magnetic field strength ( @xmath43 ) & 50 - 350 @xmath78 g + wind ( advection ) speed ( @xmath40 ) & 0 - 1000 km  s@xmath12 + ionized gas density ( @xmath44 ) & 50 - 500 @xmath7 + absorption fraction ( @xmath48 ) & 0.1 - 1.0 +    in regards to the molecular gas mass , there is a wide a range of observed values and there is not yet a clear consensus on which is most accurate .
varying the gas mass would have a similar ( but not identical ) effect to varying the acceleration efficiency or supernova rate . in the case of the cosmic ray electrons
, varying molecular gas mass affects the losses due to ionization and bremsstrahlung but not inverse comtpon or synchrotron emission .
though we adopt the value of @xmath18 @xmath5 , there is likely at least a factor of two uncertainty in this value . as such
, we also test models at @xmath84 @xmath5 . while we vary the molecular gas mass , we keep the ionized gas mass constant . however , we do vary the ionized gas density ( @xmath85 @xmath7 ) , and thus the filling fraction of ionized gas , as it has a significant impact on free - free absorption and emission in the radio spectrum .    to distinguish between models , we use @xmath86 tests to compare to the observed data .
we can see from the contour plots of the @xmath86 values that the @xmath0-ray and radio data constrain the value for the magnetic field strength and wind speed in different ways ( see figure 1 ) .
there is a degeneracy for the magnetic field strength and wind speed for the radio data such that as the magnetic field strength increases , an increase in wind speed allows for a similarly good fit .
the reverse is true for the @xmath0-rays . in the radio spectrum ,
an increase in magnetic field strength results in an increase in synchrotron emission such that the number of electrons must be reduced ( see yegz for further explanation ) .
however , in the @xmath0-ray spectrum , an increase in the magnetic field strength leads to a reduction of the inverse compton flux as fewer electrons are available as a result of the shorter energy loss lifetimes .
thus , the wind speed must be reduced to allow for more electrons to be available to produce inverse compton emission .    because the results are essentially non - complementary , we were not able to find a joint solution ( a best - fit model for both the @xmath0-ray and radio spectra based on a common set of parameters ) for the radio and @xmath0-ray data assuming @xmath18 @xmath5 for the molecular mass .
this is highlighted in figure 2 , where the combined reduced @xmath86 result is dominated by the radio data . for the radio spectrum ,
optimal wind speeds are essentially the no wind case ( @xmath40 ) which is not supported by observations of ngc 253 . for the @xmath0-ray spectrum , optimal wind speeds range from @xmath87 km  s@xmath12 , which is slightly higher than optical determinations of the wind speeds .
additionally , while @xmath0-ray models allow for a wide range of magnetic field strengths , radio models constrain the magnetic field strength to @xmath88 @xmath78 g .
these field strengths correspond to @xmath222000 to @xmath223000 ev  @xmath7 .
cosmic ray energy densities are nearly an order of magnitude below this , @xmath22200 to @xmath22300 ev  @xmath7 for @xmath89 . while the magnetic field energy density lies at the high end of the range of the assumed radiation field energy densities , the cosmic ray energy densities are within a factor of two of the lower end of the range .
taken at face value , these model results suggest that for the system the magnetic field and radiation field are in approximate equipartition while the cosmic rays are sub - equipartition .
in addition to limiting the possible magnetic field strength and wind speed values , the @xmath0-ray data severely constrain the possible combinations of spectral index and acceleration efficiency . while the radio data is fit equally well by all combinations of acceleration efficiency and spectral index , the @xmath0-ray spectrum is easily overestimated because of the large radiation fields and high densities in starburst regions . as such
, the majority of the combinations of acceleration efficiency and spectral index do not result in a minimized ( best - fit ) @xmath86 value for the @xmath0-ray spectrum .
so , to keep the combination of pion decay and inverse compton emission in check , our model rules out acceleration efficiencies of @xmath90 and @xmath91 for both spectral indices .
we also rule out the spectral index of @xmath92 for all acceleration efficiencies , as the best - fit models underestimate the @xmath0-ray spectrum at tev energies ( see figure 3 ) .
thus , the only combination of acceleration efficiency and spectral index which both minimizes and fits the tev energy data is @xmath89 , @xmath93 .
equivalently , a larger acceleration efficiency pared with a lower supernova rate would have the same result , as we constrain their product .
the goodness of our @xmath86 fits also change slightly depending on the radiation field energy density .
the best - fit @xmath0-ray spectrum occurs for an energy density of 2000 ev  @xmath7 for a spectral index of @xmath89 .
the resulting @xmath0-ray spectrum is dominated by pion decay emission but still has a significant contribution from inverse compton emission at tev energies ( see figure 3 ) .    aside from the parameters
discussed thus far , we also vary the fraction of synchrotron emission absorbed in ionized gas clouds ( see table 2 ) .
observations show that at low frequencies , the radio spectrum for ngc 253 levels off , unlike the radio spectrum for m82 which turns over at low frequencies . in the case of m82 , to fit the radio spectrum properly , we required that all of the synchrotron emission from the hot , diffuse medium was absorbed in ionized gas clouds such that @xmath94 ( yegz ) .
however , as the radio spectrum for ngc 253 only flattens , we varied the absorption fraction to find what absorption fraction was required to accurately fit the spectrum .
we found that models with an absorption fraction of up to @xmath95 agreed with all compared radio data points ( see figure 4 ) , including models with no absorption fraction ( @xmath96 ) .
the best - fit models were for an absorption fraction of @xmath97 , which was selected for use for the remainder of our models .
in yegz , we found that while m82 was an electron calorimeter , it was , at best , only a partial ( @xmath2250% ) proton calorimeter .
for ngc 253 , we find that it is an electron calorimeter and that it can be a proton calorimeter , depending upon the selected parameters .
we based our conclusions on a comparison between energy - loss and advection timescales . in ngc 253 , we find that the average wind speeds for the best - fit @xmath0-ray models are on the order of @xmath87 km  s@xmath12 . the best - fit wind speed decreases depending on the assumed magnetic field strength . comparing timescales , we find that ngc 253 is a @xmath2250% proton calorimeter .
however , the best - fits for the radio spectrum were for the no wind case . in this case , ngc 253 would be a complete proton calorimeter ( cf .
while the @xmath0-ray data were readily fit with our model , the radio spectrum presented much more of a challenge .
as previously stated , we were not able to find a joint solution for the @xmath0-ray and radio models .
while our best - fit for the @xmath0-ray spectrum had a reduced @xmath86 value of 2.0 , the best - fit for the radio spectrum had a reduced @xmath86 value of 50 .
in addition to the preferred wind speed being unphysical , the best - fit model does not match up well with the @xcite data .
although we did not explicitly vary molecular gas mass in our models , we did originally test a lower gas mass .
models tested with a lower gas mass of @xmath98 @xmath5 resulted in significantly better radio fits ( see figures 2 and 4 ) .
our best - fit model with the lower mass has a reduced @xmath86 value of 12 .
this reduction in @xmath86 with the lower mass is due to higher bremsstrahlung losses for higher masses .
thus , for the lower mass , more electrons are available to produce synchrotron emission at the low frequencies instead of bremsstrahlung .
the discrepancy between fits for the different gas masses is , in large part , due to the limitations of the one - zone model .
while there are many computational advantages to assuming a simple single zone , there are also disadvantages . in yegz
, we noted that by assuming a single ionized gas density in our model , we were unable to fit both the high - frequency radio data and the low - frequency turn down . in ngc 253
, our one zone model has reached the calorimeter limit for the higher gas mass . as such , while we are producing the same amount of secondary electrons and positrons at both masses , there are fewer electrons / positrons available to produce inverse compton and synchrotron emission due to an increase in bremsstrahlung .
the cmz dust temperature is the key to modeling the @xmath0-ray observations for ngc 1068 as it determines the radiation field that inverse compton emission depends on .
observations by @xcite show a blackbody spectrum with temperatures in the range of @xmath99 for the inner cnd .
when assuming a radiation field from dust with @xmath100 k , the photon number is significantly decreased such that our models produce negligible inverse compton @xmath0-ray emission . however , this dust temperature is attributed to the dusty torus of the agn nucleus and likely does not dominate the larger , surrounding cmz . as such , we assume that the dust temperature in the cmz is on par with the temperatures of molecular gas in the region , @xmath22100 k , and we use this to determine the radiation field spectrum .    as with ngc 253
, we intended to test a variety of different sets of parameters with which to model ngc 1068 .
however , we found ngc 1068 significantly harder to model than ngc 253 .
the upper bound on the supernova rate produces a @xmath0-ray spectrum that is lower by a factor of only a few ( see figure 5 ) .
however , a lower bound on the supernova rate results in a @xmath0-ray spectrum that is nearly two order of magnitudes lower than the observed data .
because we were underestimating the @xmath0-ray emission , we selected parameters to maximize the inverse compton emission ( a magnetic field strength of @xmath101 @xmath78 g and a radiation field energy density of @xmath102 ev  @xmath7 ) and the pion decay emission and bremsstrahlung ( a wind speed of @xmath103 km  s@xmath12 ) . even selecting parameters to augment the @xmath0-ray emission , without invoking an extra source of cosmic rays , we were not able to produce a model which agrees with the _ fermi _ observations to better than a factor of a few .    further complicating matters is the radio spectrum for ngc 1068 .
while the galaxy has been extensively observed in the radio spectrum , the presence of a radio jet greatly overshadows any emission not originating from the agn or its jets .
ultimately , we chose to compare our radio models against a few different radio observations .
first , we plot our radio models against radio observations of the agn core ( s1 ) found in @xcite . though the observations of the agn are on the parsec - scale , they do highlight the fact that they are of a fundamentally different nature than the emission from star - forming regions .
we also plot radio observations of the region labeled s2 @xcite , which could be the counter - jet as seen through the star - forming disk .
additionally , we use the 18  cm radio map in @xcite to compare with our models .
while the outermost contour includes the spiral arms , the second contour should give us an upper limit on the radio continuum in the cmz excluding the component from the agn and the jets .
taking the brightness level ( of 74 mjy / beam ) and multiplying by the area of the cmz , we obtain an upper limit of @xmath22 13 mjy ( see red squared in figure 5 ) .
while nearly all of the energy of the cosmic ray electrons goes into inverse compton emission , some small fraction of the energy is lost due to bremsstrahlung and synchrotron emission . assuming a moderate starburst magnetic field strength of 200  @xmath78 g results in a significant amount of synchrotron emission .
as the upper limit ( red square ) that we plot is nearly an order of magnitude below the model , we must be significantly overestimating the synchrotron emission in the cmz .
thus , while assuming the upper limit on the supernova rate results in a @xmath0-ray spectrum which comes within a factor of a few of the observed data , the radio spectrum disagrees significantly with observations ( see figure 5 ) .
previously , we developed and tested a model for cosmic ray interactions in starburst nuclei ( yegz ) . in testing our model against other nuclei and combined starburst / agn environments ,
we have applied our updated model to two giant , barred spiral galaxies ; specifically , we have modeled both ngc 253 , with an unambiguous starburst nucleus , and ngc 1068 , the archetypal seyfert 2 galaxy .
the @xmath0-ray spectrum for ngc 253 is well fit with a combination of emission from neutral pion decay , bremsstrahlung , and inverse compton , requiring both hadronic and leptonic emission mechanisms .    while we find a joint solution for the @xmath0-ray and radio spectra for the lower molecular gas mass for the ngc 253 cmz , we were unable to do the same at a higher gas mass ( @xmath18 @xmath5 )
. @xcite find best - fit models for both the @xmath0-ray and radio spectra for similarly large densities ( @xmath104 @xmath7 ) .
however , they do not include wind speed as a free parameter and fix their advection / diffusion timescale at 10 myr ( which is equivalent to a wind speed of @xmath22 5 km  s@xmath12 for a half - thickness of 50 pc ) .
this would allow for a good fit to the radio spectrum at higher masses and their lower radiation field would also allow for a similarly good fit for the @xmath0-ray spectrum .
other models include those of @xcite who find joint solutions when assuming parameters closer to our lower gas mass and wind speeds between 300 and 600 km  s@xmath12 .
we believe that the data support our choice of gas masses in the @xmath4 @xmath5 range as well as wind speeds of @xmath105 km  s@xmath12 ( see * ? ? ?
* ; * ? ? ?
* ) and that our models therefore contain astrophysically reasonable parameters .    though located nearly four times further away , ngc 1068 has a @xmath0-ray flux comparable to that of ngc 253 . in comparing our model with observations , we find that the inner starburst region in ngc 1068 does not supply a sufficient cosmic ray population to produce the observed @xmath0-ray emission by itself .
this agrees with analysis of the supernova rate , total gas mass , and @xmath0-ray luminosity by @xcite for ngc 1068 in comparison to other starburst galaxies . while a scenario which boosts the cosmic ray population in the cmz , such as a diffusion of cosmic rays from the circumnuclear ring into the cmz , is possible , it is most likely that the high - energy emission results from interactions in the jet .
an agn jet origin of the @xmath0-ray emission in ngc 1068 agrees with calculations of @xcite who propose a leptonic external inverse compton ( eic ) model for the high - energy emission .
while proton - proton interactions dominate as a cosmic ray energy loss mechanism in the ism , proton-@xmath0 interactions will be dominant in the jet due to the lack of molecular gas necessary for proton - proton interactions .
proton - proton interactions result in both charged and neutral pions . while neutral pions decay quickly into @xmath0-rays , charged pions decay into muons and neutrinos followed by electrons and positrons and more neutrinos
. thus , galaxies whose @xmath0-ray spectra are dominated by or are significantly contributed to by pion decay will also be sources of neutrinos .
as proton-@xmath0 interactions still result in pion production , agn could also be a source for neutrinos as starbursts are , if @xmath0-ray emission in the jet is hadronic .    while a detection of neutrinos from agn would help to determine the origin of @xmath0-ray emission
, such a detection will also provide vital clues for starburst galaxies .
observations of neutrinos provide direct information about the nature of the @xmath0-ray spectrum and the cosmic ray proton population .
better constraining the contribution of neutral pion decay to the @xmath0-ray spectrum , also better constrains the contribution of bremsstrahlung and inverse compton emission . in the case of inverse compton , this would be another diagnostic for determining the radiation field , which can be uncertain in such environments .    to make a neutrino flux prediction , we start with the source function for muon neutrinos which results from charged pion decay ( @xmath106 and @xmath107 ) depends on the source function for pions : @xmath108 where @xmath109 is the source functions for pions ( see yegz for details ) and @xmath110 @xcite .
we found that the @xmath0-ray spectra for ngc 253 had a significant contribution from pion decay .
while this results in a significant neutrino flux ( see figure 6 ) , the steep spectrum of ngc 253 means that m82 is a more likely candidate for eventual detection of the energetic neutrinos for which icecube has good sensitivity @xcite .
though all of our models for ngc 1068 underestimated the @xmath0-ray emission , we find that the corresponding neutrino flux from the starburst region is comparable to that of ngc 253 at tev energies ( see figure 6 ) . however , if the high - energy emission in the jet is hadronic , ngc 1068 may still be a detectable source of neutrinos , particularly if the flat @xmath0-ray spectrum seen in figure 5 continues to higher energies . at energies below 1 gev , models predict a significant difference in the spectra for inverse compton and pion decay @xmath0-ray emission @xcite .
thus , in the future , hard x - ray observations , such as those by nustar , should be able to distinguish between leptonic and hadronic emission mechanisms for @xmath0-ray production in such galaxies and determine whether these galaxies are likely to be detectable in neutrinos .
in this paper , we explore the differences between starburst galactic nuclei and agn and their impact on cosmic ray interactions . as a typical nuclear starburst environment , ngc 253
can be fit with our model that includes inverse compton emission .
we find that for an assumed gas mass of @xmath111 @xmath5 , the @xmath0-ray and radio spectra can not uniquely determine the magnetic field strength or wind speed .
additionally , the @xmath0-ray spectrum also severely limits the possible combinations of spectral index ( to @xmath89 ) and acceleration efficiency ( to @xmath93 ) due to the powerful combination of effects due to inverse compton and pion decay emission .
past models for ngc 253 have assumed a fixed wind ( advection ) speed which presupposes calorimetry ( e.g. * ? ? ?
based on our results , a wind is a critical component for determining the degree to which a starburst system is a calorimeter . while the best - fits to the @xmath0-ray spectrum result in a @xmath2250% proton calorimeter , fits to the radio spectrum result in a complete proton calorimeter
. additionally , our fits result in systems below equipartition ( in regards to the cosmic ray energy densities ) with strong magnetic fields ( see also * ? ? ?
the lack of a solution which fits the @xmath0-ray and radio spectra together with a single set of parameters and the problems in fitting the radio spectrum highlight the limitations of this one - zone model .
we , however , can obtain a good fit to @xmath0-rays and the radio by assuming a lower gas mass .
understanding of the ism conditions is vital to applying these models to starbursts . as such , while the model works well for a late phase starburst such as m82 , it does not as easily fit an early phase starburst such as ngc 253 .    as an agn environment containing a possible starburst region , ngc 1068
is not fit with our self - consistent starburst model .
while we find a few models that come close to predicting the observed @xmath0-ray spectrum for ngc 1068 , all of our models underestimate the high - energy emission and we are unable to find any models of the inner starburst region that agree with the observed data .
thus , we find that the cosmic ray population from the starburst region alone is not enough to accurately fit the observed @xmath0-ray emission and even those models in rough agreement with high - energy observations result in an overestimation of the radio spectrum .
we conclude that the observed @xmath0-ray emission likely occurs in the agn or its jets .    additionally , we find that starburst galaxies such as ngc 253 are potential sources for high energy neutrino detection
. however , agn with high @xmath112 such as ngc 1068 will not be luminous neutrino sources unless their high - energy emission is of hadronic origin .
this work was supported in part by nsf ast-0907837 , nsf phy-0821899 ( to the center for magnetic self - organization in laboratory and astrophysical plasmas ) , and nsf phy-0969061 ( to the icecube collaboration ) .
we thank alberto bolatto for conversations concerning ngc 253 .
we also thank francis halzen for his help and support .
abbasi , r. , abdou , y. , abu - zayyad , t. , et al .
2011 , , 732 , 18 abramowski , a. , acero , f. , aharonian , f. , et al .
2012 , , 757 , 158 ackermann , m. , ajello , m. , allafort , a. , et al .
2012 , , 755 , 164 anchordoqui , l.  a. , hooper , d. , sarkar , s. , & taylor , a.  m.  2008 , aph , 29 , 1 blumenthal , g.  r. , & gould , r.  j.  1970 , rev .
, 42 , 273 bock , j.  j. , neugebauer , g. , matthews , k. , et al .
2000 , , 120 , 2904 boettcher , e. , zweibel , e.  g. , yoast - hull , t.  m. , & gallagher , j.  s. , iii  2013 , , 779 , 12 carilli , c.  l.  1996 , , 305 , 402 condon , j.  j.  1992 , , 30 , 575 davidge , t.  j.  2010 , , 725 , 1342 davies , r.  i. , sugai , h. , & ward , m.  j.  1998 , , 300 , 388 domingo - santamara , e. , & torres , d.  f.  2005 , , 444 , 403 engelbracht , c.  w. , rieke , m.  j. , rieke , g.  h. , kelly , d.  m. , & achtermann , j.  m.  1998 , , 505 , 639 exposito , j. , gratadour , d. , clnet , y. , & rouan , d.  2011 , , 533 , 63 fabbiano , g. , & trinchieri , g.  1984 , , 286 , 491 forbes , d.  a. , polehampton , e. , stevens , i.  r. , brodie , j.  p. , & ward , m.  j.  2000 , , 312 , 689 gallimore , j.  f. , baum , s.  a. , & odea , c.  p.  1996a , , 464 , 198 gallimore , j.  f. , baum , s.  a. , odea , c.  p. , & pedlar , a.  1996b , , 458 , 136 gallimore , j.  f. , baum , s.  a. , & odea , c.  p.  2004
, , 613 , 794 hailey - dunsheath , s. , nikola , t. , stacey , g.  j. , et al .
2008 , , 689 , l109 hailey - dunsheath , s. , sturm , e. , fischer , j. , et al .  2012 , , 755 , 57 halzen , f.  2006 , epjc , 46 , 669 heesen , v. , beck , r. , krause , m. , & dettmar , r .- j .
2009a , , 494 , 563 heesen , v. , krause , m. , beck , r. , & dettmar , r .- j .  2009b , , 506 , 1123 hnig , s.  f. , prieto , m.  a. , & beckert , t.  2008 , , 485 , 33 kamenetzky , j. , glenn , j. , maloney , p.  r. , et al .
2011 , , 731 , 83 karachentsev , i.  d. , grebel , e.  k. , sharina , m.  e. , et al .
2003 , , 404 , 93 karachentsev , i.  d. , tully , r.  b. , dolphin , a. , et al .
2007 , , 133 , 504 kelner , s.  r. , aharonian , f.  a. , & bugayov , v.  v.  2006 , , 74 , 034018 kornei , k.  a. , & mccrady , n.  2009 , , 697 , 1180 krips , m. , martn , s. , eckart , a. , et al .
2011 , , 736 , 37 lacki , b.  c. , & beck , r.  2013 , , 430 , 3171 lacki , b.  c. , thompson , t.  a. , & quataert , e.  2010 , , 717 , 1 lacki , b.  c. , thompson , t.  a. , quataert , e. , loeb , a. , & waxman , e.  2011 , , 734 , 107 lehmer , b.  d. , wik , d.  r. , hornschemeier , a.  e. , et al .
2013 , , 771 , 134 lenain , j .-
p . , ricci , c. , trler , m. , dorner , d. , & walter , r.  2010 , , 524 , 72 lenc , e. , & tingay , s.  j.  2006 , , 132 , 1333 mitsuishi , i. , yamasaki , n.  y. , & takei , y.  2011 , , 742 , l31 mller - snchez , f. , gonzlez - martn , o. , fernndez - ontiveros , j.  a. , acosta - pulido , j.  a. , & prieto , m.  a.  2010 , , 716 , 1166 muxlow , t.  w.  b. , pedlar , a. , holloway , a.  j. , gallimore , j.  f. , & antonucci , r.  r.  j.  1996 , , 278 , 854 paglione , t.  a.  d. , & abrahams , r.  d.  2012 , , 755 , 106 paglione , t.  a.  d. , marscher , a.  p. , jackson , j.  m. , & bertsch , d.  l.  1996 , , 460 , 295 peng , r. , zhou , s. , whiteoak , j.  b. , lo , k.  y. , & sutton , e.  c.  1996 , , 470 , 821 persic , m. , rephaeli , y. , & arieli , y.  2008 , , 486 , 143 ricci , r. , prandoni , i. , gruppioni , c. , sault , r.  j. , & de zotti , g.  2006 , , 445 , 465 rybicki , g.  b. , & lightman , a.  p.  1979
, radiative processes in astrophysics ( new york : wiley ) sakamoto , k. , mao , r .- q . ,
matsushita , s. , et al .
2011 , , 735 , 19 schinnerer , e. , eckart , a. , tacconi , l.  j. , genzel , r. , & downes , d.  2000 , , 533 , 850 schlickeiser , r.  2002 , cosmic ray astrophysics ( new york : springer ) schlickeiser , r. , vainio , r. , bttcher , m. , et al .
2002 , , 393 , 69 schulz , h. , & wegner , g.  1992 , , 266 , 167 scoville , n.  z. , matthews , k. , carico , d.  p. , & sanders , d.  b.  1988 , , 327 , l61 spinoglio , l. , pereira - santaella , m. , busquet , g. , et al .
2012 , , 758 , 108 stecker , f.  w.  1979 , , 228 , 919 strickland , d.  k. , heckman , t.  m. , colbert , e.  j.  m. , hoopes , c.  g. , & weaver , k.  a.  2004 , , 151 , 193 strickland , d.  k. , heckman , t.  m. , weaver , k.  a. , dahlem , m.  2000 , , 120 , 2965 storchi - bergmann , t. , riffel , r.  a. , riffel , r. , et al .
2012 , , 755 , 87 telesco , c.  m. , & harper , d.  a.  1980 , , 235 , 392 thompson , t.  a. , quataert , e. , & waxman , e.  2007 , , 654 , 219 thronson , h.  a. , jr . ,
hereld , m. , majewski , s. , et al .
1989 , , 343 , 158 torres , d.  f.  2004 , , 617 , 966 ulrich , m .- h .  1978 , , 219 , 424 ulvestad , j.  s. , & antonucci , r.  r.  j.  1997 , , 488 , 621 van buren , d. , & greenhouse , m.  a.  1994 , , 431 , 640 waller , w.  h. , kleinmann , s.  g. , ricker , g.  r.  1988 , , 95 , 1057 wei , a. , kovcs , a. , gsten , r. , et al .
2008 , , 490 , 77 westmoquette , m.  s. , smith , l.  j. , & gallagher , j.  s. , iii  2011 , , 414 , 3719 williams , p.  k.  g. , & bower , g.  c.  2010 , , 710 , 1462 wilson , a.  s. , & ulvestad , j.  s.  1983 , , 275 , 8 yoast - hull , t.  m. , everett , j.  e. , gallagher , j.  s. , iii , & zweibel , e.  g.  2013 , , 768 , 53 ( yegz )