abstract: We investigate the consequences of a model of the assembly and growth of massive black holes (MBHs) from primordial seeds, remnants of the first generation of stars in a hierarchical structure formation scenario. Our model traces the build-up of MBHs from an early epoch, and follows the merger history of dark matter halos and their associated holes via Monte Carlo realizations of the merger hierarchy from early times to the present time in a $\Lambda$CDM cosmology. The sequence of minor and major mergers experienced by galactic halos in their hierarchical growth affects the merger history of MBHs embedded in their nuclei. So, if the formation route for the assembly of SMBHs dates back to the early universe, a large number of BH interactions is inevitable. Binary black holes coalescence timescales can be long enough for a third BH to fall in and interact with the central binary. These BH triple interactions lead typically to the final expulsion of one of the three bodies and to the recoil of the binary. Also, asymmetric emission of gravitational waves in the last stages of the black hole merging can give a recoil velocity to the centre of mass of the coalescing binary. This scenario leads to the prediction of a population of intermediate mass BHs (IMBHs) wandering in galaxy halos at the present epoch. We compute the luminosity distribution produced by these IMBHs accreting from their circumstellar medium. We find that in a Milky Way-sized galaxy they are unable to account for sources with luminosities $\ga 10^{39}$ erg/s unless they carry a baryonic remnant from which they are able to accrete for a long time. We also find that, for typical spiral galaxies, the bright end of the point source distribution correlates with the mass of the galaxy, and the most luminous sources are expected to be found in the disk.
author: Marta Volonteri$^{1,2}$[^1] and Rosalba Perna$^{3}$
$^{1}$Department of Astronomy & Astrophysics, University of California, Santa Cruz, CA 95064, USA
$^{2}$Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
$^{3}$Department of Astrophysical and Planetary Sciences, University of Colorado at Boulder, 440 UCB, Boulder, CO 80309, USA
title: Dynamical evolution of intermediate mass black holes and their observable signatures in the nearby Universe

# Introduction

Observations of the centres of nearby galaxies indicate that most host a supermassive black hole (Magorrian et al. 1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000); supermassive black holes (SMBHs) seem to be unequivocally linked to galactic bulges, their masses scaling with the bulge luminosity and stellar velocity dispersion. In cold dark matter (CDM) 'bottom-up' cosmogonies, mergers of galaxies are a central part of the galaxy formation process in the hierarchical structure formation picture, with major mergers (i.e. those between comparable-mass systems) expected to result in the formation of elliptical galaxies (Barnes 1988).

Volonteri, Haardt & Madau 2003 (VHM03) and Volonteri, Madau, & Haardt 2003 (VMH03) have assessed a model for the assembly of SMBHS at the centre of galaxies that trace their hierarchical build-up, assuming that the first 'seed' black holes (BHs) had intermediate masses, $m_\bullet\approx 150\,\,{\rm M_\odot}$, and formed in (mini)halos collapsing at $z\sim 20$ from high-$\sigma$ density fluctuations. These pregalactic holes evolve in a hierarchical fashion, following the merger history of their host halos. During a merger event BHs approach each other owing to dynamical friction, and form a binary system. Stellar dynamical processes drive the binary to harden and eventually coalesce.

The lifetime of BH binaries can be long enough (Begelman, Blandford & Rees; Quinlan & Hernquist 1997; Milosavljevic & Merritt 2001; Yu 2002) that following another galactic merger a third BH can fall in and disturb the evolution of the central system (VHM03). The three BHs are likely to undergo a complicated resonance scattering interaction, leading to the final expulsion of one of the three bodies ('gravitational slingshot') and to the recoil of the binary. Any slingshot, in addition, modifies the binding energy of the binary, typically creating more tightly bound systems (VHM03).

Another interesting gravitational interaction between black holes happens during the last stage of coalescence, when the leading physical process for the binary evolution becomes the emission of gravitational waves. If the system is not symmetric (e.g. BHs have unequal masses or spins) there would be a recoil due to the non-zero net linear momentum carried away by gravitational waves in the coalescence ('gravitational rocket'). The coalesced binary would then be displaced from the galactic centre, leaving straightaway the host halo or sinking back to the centre owing to dynamical friction.

Historically, the possibility of these gravitational interactions between BHs has been envisaged even before the CDM-hierarchical picture has become popular. Slingshots have been proposed, as a contrast to the unified -twin beam- theory (e.g. Blandford & Rees 1974) to explain the nature of double radio sources (Saslaw, Valtonen & Aarseth 1974; Mikkola & Valtonen 1990; Valtonen & Heinämäki 2000): ejected BHs can produce extended radio lobes by interacting with the intergalactic medium in the areas of the lobes.

Slingshots and rockets were both involved when trying to explain the nature of galactic dark halos as composed of BHs (Hut & Rees 1992; Xu & Ostriker 1994). Both processes basically give BHs a recoil velocity, that can even exceed the escape velocity from the host halo and spread BHs outside galactic nuclei (VHM03). So, the outlined scheme predicts, along nuclear SMBHs hosted in galaxy bulges, a number of *wandering BHs* that are the result of rockets, slingshots and mergers with a dynamical friction timescale longer than the Hubble time, so at $z=0$ the BHs are still on their way to the galactic centre. Note that most of the wandering BHs masses are in an intermediate range between stellar mass BHs and supermassive ones.

In this paper, we improve the modelling of the dynamical interactions between galactic halos and between BHs and their host halos. In particular, we include the lengthening of the dynamical friction timescale of infalling satellites due to tidal interactions with the potential of the primary halo, and follow the orbital decay of every single satellite, determining its position within the primary at any given time. We also trace step-by-step the trajectories of BHs ejected from galaxy centres following an energetic BH interaction; this enables us to keep track of the position and velocity of all wandering IMBHs as a function of cosmic time.

The population of wandering IMBHs in local galaxies that our model predicts will have some observational consequences. In particular, due to accretion from their surrounding interstellar medium (ISM), these BHs will appear as bright X-ray sources. We compute here the luminosity distribution expected from such sources, and relate it to current observations. In particular, we discuss whether and under what conditions these IMBHs can make up a substantial fraction of a population of ultra-luminous sources (ULXs), characterised by luminosities $L_X\ga 10^{39}$ erg/s, which has been recently identified in large *Chandra* surveys (e.g. Colbert et al. 2004). The properties of this population as a whole (see Fabbiano 1988 for a review) hint that we might be actually observing an heterogeneous population, where super-Eddington luminosities could for example be due to a photon-bubble instability (Begelman 2002) or could be apparent due to strong beaming (King et al. 2001). A contribution to the high-end tail of the X-ray point source luminosity distribution could also be due to young millisecond pulsars (Perna & Stella 2004). However, accreting IMBHs constitute an obvious candidate for the ULXs, as it has been often suggested (e.g. Colbert & Mushotzky 1999; Miller & Hamilton 2002; Cropper et al. 2004). Part of our study here is therefore a detailed investigation of the present-day observational properties of the population of IMBHs formed in mini-halos at high redshifts and evolved during the process of galaxy formation.

The paper is organised as follows: in §2, we summarise the main features of our cosmological scenario of BH growth; the dynamical evolution of SMBH binaries is discussed in §3. The dynamical evolution of satellites BHs is studied in §4, while §5 describes the accretion model. §6 is dedicated to the observational consequences of these IMBHs in present-day galaxies and §7 discusses our results in the context of previous work in the field, and finally, in §8, we summarise our work.

# Assembly and mergers of SMBHs

We briefly summarise here the main features of our scenario for the hierarchical growth of SMBHs in a $\Lambda$CDM cosmology (see VHM03 for a thorough discussion). In our model, pregalactic 'seed' holes form with intermediate masses ($m_\bullet=150\,\,{\rm M_\odot}$) in (mini)halos collapsing at $z=20$ from rare 3.5-$\sigma$ peaks of the primordial density field (VHM03; Madau & Rees 2001). The assumed 'bias' assures that almost all halos above $10^{11}\,\,{\rm M_\odot}$ actually host a BH at all epochs.

The merger history of dark matter halos and associated BHs is followed by cosmological Monte Carlo realizations of the merger hierarchy from early times until the present in a $\Lambda$CDM cosmology. The dynamical evolution of the BH population is followed in detail with a semi-analytical technique. In this paper, we focus on the properties of two sets of halos with $M=2\times10^{12}\,\,{\rm M_\odot}$ and $M=10^{13}\,\,{\rm M_\odot}$ respectively, creating 30 merger trees of the former and 10 of the latter and averaging the results. For the smaller halo, we will also consider a case in which seed holes are more numerous and populate the 3-$\sigma$ peaks instead.

In order to follow the history of halos and massive black holes, we adopt a two-component model for galaxy halos. The dark matter is distributed according to a NFW profile (Navarro et al. 1997). During the merger of two halo$+$BH systems of comparable masses, dynamical friction against the dark matter background drags in the satellite hole towards the centre of the newly merged system, leading to the formation of a bound BH binary. When two halos of mass $M$ and $M_s$ merge, the 'satellite' (less massive) progenitor (mass $M_s$) is assumed to sink to the centre of the more massive pre-existing system on the dynamical friction (against the dark matter background) timescale, which depends on the orbital parameters of the infalling satellite, which we take from van den Bosch et al. (1999).

The fate of the infalling halo has been investigated recently (Taffoni et al. 2003) in a NFW profile. Dynamical friction appears to be very efficient for 'major' mergers with (total) mass ratio of the progenitors, $P=M_s/M>0.1$. In each major merger the more massive hole accretes at the Eddington rate a gas mass that scales with the fifth power of the circular velocity of the host halo: the normalisation is fixed a posteriori in order to reproduce the observed local $m_{\rm BH}-\sigma_*$ relation (Ferrarese 2002) and the observed luminosity function of optically-selected quasars in the redshift range $1<z<5$.

Satellites of intermediate mass ($0.01 M_{\rm h}<M_{\rm
s,0}<0.1 M_{\rm h}$) suffer severe mass losses by the tidal perturbations induced by the gravitational field of the primary halo. Tidal perturbations cause the progressive stripping of the satellite and, as the magnitude of the frictional drag depends on the mass of the satellite, this progressive mass loss increases the decay time. After a Hubble time satellites have typical masses $\sim 1-10\%$ of their initial mass. The lightest satellites ($P<0.01$) seem to be almost unaffected by orbital decay, so they survive and keep orbiting on rather circular, peripheral orbits.

We adopt the relations suggested by Taffoni et al. for the orbital decay of satellites in a NFW halo. Figure shows the evolution of the dynamical friction timescales with redshift and the dependence on the mass ratio for halos both hosting massive black holes.

# Black hole binary dynamical evolution

The subsequent evolution of a BH binary is determined by the initial central stellar distribution. The binary will initially shrink by dynamical friction from distant stars acting on each BH individually. But as the binary separation decays, the effectiveness of dynamical friction slowly declines because distant encounters perturb only the binary centre of mass but not its semi-major axis. The BH pair then hardens via three-body interactions, i.e., by capturing and ejecting at much higher velocities the stars passing by within a distance of order the semi-major axis of the binary. The merger timescale is computed adopting a simple semi-analytical scheme that qualitatively reproduces the evolution observed in N-body simulations (Merritt 2000; VHM03, VMH03). The binary evolution is slowed down due to the declining stellar density, with a hardening time that becomes increasingly long as the binary shrinks. If the hardening continues sufficiently far, gravitational radiation losses finally take over, and the two BHs coalesce in less than a Hubble time.

We have adopted a simplified description of the BH coalescence process, ignoring for instance the depopulation of the loss cone (e.g. Yu 2002) and the loss of orbital angular momentum to a gaseous disc (Armitage & Natarajan 2002, Escala et al. 2004).

## Gravitational Rocket

When the binary semi-major axis is such that the timescale for hardening by stellar scatterings equals the timescale for shrinking by emission of gravitational waves, the latter process takes over. At this point, if the holes have unequal masses, or spin, the binary recoils due to the non-zero net linear momentum carried away by gravitational waves in the coalescence ('gravitational rocket').

The recoil velocity has still large uncertainties, and we adopt here the results from recent calculations by Favata et al. (2004). They estimate the recoil velocity as a function of the binary mass ratio $q=m_2/m_1$ (with BH masses $m_1\ge m_2$) and the spin parameter during final plunge and coalescence. We adopt here the expression they suggest for estimating 'upper limits' for the recoil, converting it into estimates of the binaries recoil velocity, $v_{\rm CM}$, following the procedure described in Merritt et al. (2004).

If the recoil velocity of the coalescing binary exceeds the escape speed $v_{\rm esc}$, the holes will leave the galaxy altogether. If instead $v_{\rm CM}<v_{\rm esc}$, we assume that the binary is kicked out to a radius $r_{\rm max}$ such that the total energy is conserved, i.e. we solve for $r_{\rm max}$ the equation $(m_1+m_2)\,v_{\rm CM}^2=2(|\phi_0|-|\phi(r_{\rm max})|)$, where $\phi_0$ is the galactic potential at the initial location of the binary and $\phi(r_{\rm max})$ is the potential at the apocentre.

## Gravitational slingshot

The dynamical evolution of SMBH binaries may be disturbed by a third incoming BH, if another major merger takes place before the pre-existing binary has had time to coalesce or a wandering BH can reach the centre owing to dynamical friction (e.g. Hut & Rees 1992; Xu & Ostriker 1994). If the incoming hole reaches the sphere of influence of the central binary, the three BHs are likely to undergo a complicated resonance scattering interaction, leading to the final expulsion of one, typically the lightest, of the three bodies (gravitational slingshot).

Typically an encounter with an intruder of mass $m_{\rm int}$ smaller than both binary members leads to a scattering event, where the binary recoils by momentum conservation and the incoming lighter BH is ejected from the galaxy nucleus. By contrast, when the intruder is more massive than one or both binary components, the probability of an exchange is extremely high: the incoming hole becomes the member of a new binary, and the lightest BH of the original pair gets ejected (Hills & Fullerton 1980).

All triple interactions are followed along the merger tree, adopting the same scheme described in VHM03.

In all cases we have estimated the recoil velocities conserving the binary energy and momentum and estimated the apocentre $r_{\rm max}$ solving the equation $(m_1+m_2)\,v_{\rm CM}^2=2(|\phi_0|-|\phi(r_{\rm max})|)$. Dynamical friction timescales have been calculated by the Chandrasekhar dynamical friction timescale for circular orbits, so they are presumably upper limits to these timescales (thus giving an upper limit on the number of wandering BHs), as the ejection probably occurs on eccentric orbits.

In figure the distribution of BH recoil velocities is plotted in units of the stellar velocity dispersion, ranging from $5 \,{\rm km\,s^{-1}}$ to $200 \,{\rm km\,s^{-1}}$ depending on redshift and halo mass. As a rule of thumb, if the recoil velocity is smaller than the velocity dispersion, the BH remains in the proximity of the centre and the orbital decay is fast (Madau & Quataert 2004). Figure shows that the ejection of all three BHs is therefore not a likely outcome under the assumptions of our BH growth model (see VHM03 for a thorough discussion).

Clearly, most of the single holes either escape their host halos or are slung to the periphery of the galaxy with consequently long dynamical friction timescales. Most of the binaries recoil instead within the core and fall back to the centre soon afterwards, with $t_{\rm df}<0.01\,$Gyr.

Note that triple interactions are rare events, in the whole history of a MW-sized halo with its satellites, on average only $\sim 3$ slingshots happen, in the case of seed holes in 3.5-$\sigma$ peaks.

# Dynamical evolution of off-nuclear black holes

The assumptions underlying our scenario lead to the prediction of a population of massive BHs *wandering* in galaxy halos, due to the processes described above.

During the evolution and mass growth of DM halos, we track the position of all wandering BHs, following the modifications of the galactic potential due to mergers. In case of minor mergers, the potential varies slightly, so the wandering BHs orbits evolve smoothly, while in case of major mergers the potential of both interacting systems change abruptly, causing larger modifications to the BH orbits. Although numerical simulations with resolution high enough to follow the dynamical evolution of intermediate mass BHs in galactic mergers are not available, we can rely on simulations studying the mixing of particles in major mergers (White 1980, Nipoti et al 2003), assuming that during a merger a BH follows the evolution of other particles. On average particles which were at a given radius (normalised at the half mass radius) in one of the galaxies before the merger, are found at a similar radius (normalised to the new half mass radius) in the merger remnant.

To track the orbital evolution of wandering BHs during galactic mergers we adopt the following scheme:

-   MINOR MERGERS: we let the wandering BHs in the more massive progenitor unperturbed, while the wandering BHs in the satellite galaxy are assumed to decay with their host.

-   MAJOR MERGERS: each wandering BH outside the core is placed in the new halo at a distance conserving the same ratio with the virial radius as it had in its parent system. If a wandering BH is instead inside the core, it is assumed to reach the centre of the newly formed halo within the same merging timescale of the halos (see section 2).

# Detecting Intermediate-Mass Black Holes in local galaxies

## Naked wandering BHs

As our first model of IMBHs in present-day galaxies, we assume that all BHs are *naked*, i.e. that the BHs left over by minor mergers have been completely stripped of their envelope. We embed all wandering BHs in a galactic disc/bulge/halo model (Springel & White 1999), and determine the corresponding density of the medium in the surrounding of the BHs at their respective positions. We assume that a fraction $m_{\rm d}$ of the initial halo mass is into a thin stellar disc with an exponential surface density, $$\rho_{\rm d}(R,z) = \frac{\Sigma_0}{2z_0} \exp\left(-\frac{R}{R_{\rm d}}\right) {\rm
sech}^2\left(\frac{z}{z_0}\right) .$$ Here $M_{\rm d}=m_{\rm d} M$ is the total mass of the disc and $R_{\rm d}$ is its scale radius. We adopt $m_{\rm d}=0.05$, while the disc scale length is determined through the fitting formula provided in Shen et al. (2003) for the Petrosian half-light (or mass) radius: $$R_{\rm d} = 0.17 M_{\rm d}^{0.14}  (1+M_{\rm d}/3.98\times10^{10})^{0.25}.$$ Most spiral galaxies seem to be consistent with a constant vertical scale length with a value of $z_0\simeq 0.2 R_{\rm d}$. Bulges are modelled with a spherical Hernquist (1990) profile of the form $$\rho_{\rm b}(r)=
\frac{M_{\rm b}}{2\pi}\frac{r_{\rm b}}{r(r_{\rm b}+r)^3} .$$ In analogy to the treatment of the disc, we assume that the bulge mass is a fraction $m_{\rm b}=0.02$ of the halo mass and that the bulge scale radius $r_{\rm b}$ is a fraction $f_{\rm b}=0.1$ of that of the disc, i.e. $r_{\rm b}=f_{\rm
b} R_{\rm d}$. We have then assumed that the gas traces the baryons, the gas fraction (of the total baryonic mass) is $20\%$ in the disk and $4\%$ in bulge and halo (Bell & de Jong 2000, Fukugita, Hogan & Peebles 1998).

To model an elliptical galaxy (in the case $M=10^{13}\,\,{\rm M_\odot}$) we have considered the Hernquist profile only, assuming that the baryonic content of the galaxy is a fraction $0.1$ of the halo mass. For the elliptical, we have then assumed that the gas fraction is $1\%$(Mathews & Brighenti 2003).

## Accretion model

In this section, we discuss how we estimate the luminosities of off-nuclear BHs resulting from the process of accretion from the surrounding medium.

Let $\rho$ be the density of the medium in the surrounding of a BH, and let $c_s$ be its sound speed. The accretion rate onto a BH of mass $M_{\rm BH}$ can be estimated using the Bondi-Hoyle formula (Bondi & Hoyle 1944), $$\dot{M}_{\rm Bondi}=\frac{\lambda \,4\pi\, G^2\,
M_{\rm BH}^2\,\rho}{(v^2+c_s^2)^{3/2}}\; ,
\label{eq:dotm}$$ where $v$ is the velocity of the BH with respect to the medium, and $\lambda$ is a parameter on the order of 1. The BH radial velocities are computed by integrating the Chandrasekar formula (Binney & Tremaine 1987) as a function of the initial conditions, while the sound speed velocity, $c_s\sim 10(T/10^4{\rm K})^{1/2}$ km/s, depends on the local temperature at the BH location, which is discussed later in this section.

The bolometric luminosity of the BH can be written as $$L_{\rm bol}=\eta \dot{M} c^2\;,
\label{eq:lbol}$$ where $\eta$ represents the fraction of the accreted mass that is radiated away. The nature of the accretion process, and the consequent value of $\eta$ in BHs is rather uncertain. Active Galactic nuclei (AGNs) accrete through accretion discs with an high efficiency ($\eta\sim 0.1$). Their luminosities are close to their Eddington values. At the other end of the luminosity function, are the SMBHs at the centres of our own and nearby galaxies, whose luminosities can be as low as $\sim 10^{-9}-10^{-8}$ of their Eddington values (e.g. Loewenstein et al. 2001). Low luminosities could be due either to low efficiencies with $\dot{M}\sim \dot{M}_{\rm Bondi}$, or to relatively high efficiencies with $\dot{M}\ll \dot{M}_{\rm Bondi}$ (or also to a combination of both low accretion rates and low efficiencies). Both scenarios have been suggested in the literature. Narayan & Yi (1995a) found that, when the accretion rate decreases below some critical value $\dot{m}\sim 0.1\alpha^2$ ($\alpha$ being the viscosity parameter in the disc), the thin disc switches to an Advection Dominated Accretion Flow, characterised by a low radiative efficiency, due to the inability of the electrons in the plasma to cool (see also Rees et al. 1982; Begelman & Celotti 2004). On the other hand, an accretion rate at values below the Bondi rate could be due to the presence of winds (or outflows; Narayan & Yi 1995b; Blandford & Begelman 1999). This scenario appears to be supported by observations of both accreting black holes (Loewenstein et al. 2001), as well as old, accreting neutron stars (Perna et al. 2003).

In this work, when considering a model for the X-ray luminosity from accretion onto the IMBHs similar to that in SMBHs, we will adopt an empirical approach, and simply assume that the X-ray luminosity is a fraction $\eta_{\rm X}$ of the accretion luminosity, with $\eta_{\rm X}$ calibrated on X-ray observations of SMBHs. This approach is similar to the one adopted by Agol & Kamionkowski (2002), in their calculations of the X-ray luminosity from isolated, stellar-mass BHs in the Milky Way (they parameterised their results in terms of the X-ray accretion efficiency).

As a second model, we consider an optically thick, geometrically thin accretion disc, as in the standard Shakura-Sunyaev (1973, SS) formulation. Radiation is emitted locally with a blackbody spectrum of temperature $$T(r) = \left(\frac{3GM_{\rm BH}\dot{M}}{8\pi r^3\sigma}\right)^{1/4}\left[
1-\left(\frac{r_{\rm in}}{r}\right)^{1/2}\right]^{1/4}\;.
\label{temp}$$ The emitted spectrum (per unit frequency) from the whole surface of the disc is given by $$L_\nu= 2\pi \frac{h\nu}{c^2}\int_{r_{\rm in}}^{r_{\rm max}}\frac{r\, dr}
{e^{h\nu/k T(r)}-1}\;.
\label{eq:bb}$$ We take $R_{\rm in}=3 R_S$ and $R_{\rm max}=1000 R_S$, where $r_s=2GM_{\rm BH}/c^2$ is the Schwarzschild radius of the BH; this yields an efficiency $\eta\sim0.06$. As described in the previous sections, our simulations allow us to identify whether the BHs are in the halo, bulge or disc of the galaxy. For those in the halo and bulge, we assume that the temperature is the virial temperature of the accreting interstellar medium (ISM). This is $T_{\rm vir}=5.8\times 10^5$ K for the galaxy of mass $2\times 10^{12}\,{\rm M_\odot}$, and $T_{\rm vir}=2.3\times 10^{6}$ K for the galaxy of mass $10^{13}\,{\rm M_\odot}$. In our Galaxy, we know that the disc is composed of several phases (e.g. Bland-Hawthorn & Reynolds 2000), from cold ($T\sim 100$ K), to warm ($T\sim 8000$ K), to hot ($T\sim 10^6$ K). Modelling the various phases is beyond the scope of this paper essentially because, while some information on these phases is available for our Galaxy, no such studies have been possible for other galaxies. Moreover, due to the process of radiative feedback (Blaes, Warren & Madau 1995), the gas in the vicinity of the BH will be heated by the ionising radiation produced by accretion. The amount of heating and cooling depends on a number of factors, such as the density of the medium, its metallicity, the velocity of the BH, the accretion efficiency. For a medium made of H and He only, the ratio between the total heating rate to the initial energy density of the gas at the accretion radius is $\sim 2\times 10^5\eta v_{40}^{-2}T_4^{-1}$ (Agol & Kamionkowski 2002), where $T_4=T_410^4$ K, and $v=v_{40}40$ km/s is the BH velocity. For BHs with velocities smaller than a few tens of km/s, heating is important even for very low accretion efficiencies. Once H and He are ionized, the opacity drops substantially, and the heating rate becomes much longer than the accretion rate, leaving the gas at a temperature of around $2\times 10^4$ K (Agol & Kamionkowski 2002; Blaes, Warren & Madau 1995). We will take this to be the temperature of the accreting gas in the disc, assuming that, if there is a hot phase, this does not generally dominate. Note that, in those cases where the velocity of the BH is sufficiently high that dynamical heating is not important (and therefore the gas is actually cooler than assumed), the velocity of the BH dominates over that of the sound speed in the expression for the accretion rate, hence making the precise value of the temperature of the gas unimportant for our calculations.

For any BH, the density of the accreting material is determined at its location in the galaxy according to the model described in §6.1. For satellites, we also consider the possibility that they might accrete from a remnant core as described in §5.4.

# Results and implications for ULXs

At a given time, a galaxy not only contains off-centre BHs slingshot or rocketed, but can also have BHs, possibly embedded in the remnant of their original host, still on their way to the centre following a minor merger.

Figure shows the evolution with redshift of the mass function of wandering BHs in a $M=2\times10^{12}\,\,{\rm M_\odot}$ halo, assuming that seed holes form in 3-$\sigma$ density peaks at z=20. Figure shows instead the mass function of wandering BHs in a $M=2\times10^{12}\,\,{\rm M_\odot}$ halo at $z=0$, for the case of seed holes in 3.5-$\sigma$ peaks and in 3-$\sigma$ peaks. Note that most of the wandering BH masses are in an intermediate range between stellar mass BHs and supermassive ones. The upper end of the mass distribution of IMBHs, as well as their total number, is an increasing function of the galaxy mass. This can be seen by comparing Figure 4 (top panel) and Figure , which shows the mass function of wandering BHs in a $M=10^{13}\,\,{\rm M_\odot}$ halo at $z=0$, for the case of seed holes in 3.5-$\sigma$ peaks.

Figures and show the corresponding radial distributions of these BHs. Typically, BHs in satellite remnants dwell in the outskirts of galaxies, while those caused by slingshots or rockets have a broader distribution and can be found even in the inner regions. Interestingly, there is not a correlation between BH masses and radial distance (Figure ). Naively, more massive BHs would be expected in inner regions, due to the inverse dependence of orbital decay timescale on the mass. On the other hand, the shortness of the orbital decay timescale for the more massive BHs implies that very few of these are expected in the inner regions, because once they sink to the inner regions, they quickly sink further to the very center and merge with the central black hole or are re-ejected. Larger holes, however, appear only at later times, so that the time elapsed from their initial interaction with the galaxy is short. Figure shows also the density of gas surrounding the same BHs.

Resulting luminosities are shown in Figures 9, 10 and 11. For the bolometric luminosities, we show results for an efficiency $\eta=0.1$ (more typical of AGNs) and for $\eta=10^{-5}$ (more representative of the dormant SMBHs in local galaxies), while X-ray luminosities are shown for the disc model and for an X-ray efficiency[^2] $\eta_{\rm X}=10^{-6}$.

For our galaxy model with mass $M=2\times 10^{12}M_\odot$, and accretion from the local ISM, we find (see Figure 9) that there are no sources with bolometric luminosity larger than $\sim 10^{39}$ erg/s, even for a large efficiency of $0.1$. Clearly, these sources cannot contribute to the ultra-luminous X-ray sources, with $L_X\ga 10^{39}$ erg/s, observed in local galaxies. We find the X-ray luminosity to be at most a few $\times 10^{38}$ erg/s (see Figure 10), if accretion proceeds at high efficiency through an SS disc. We need to stress that, although we find some very massive slung black holes at typical ISM densities (see Fig. 8), they have rather large velocities, and this suppresses the accretion rate and prevents these sources from exhibiting very large X-ray luminosities It should also be noted that, although the number of BHs is largest in the halo and the smallest in the disc, the bright end of the luminosity function is dominated by the BHs in the disc when accretion is solely from the ISM. This is due to the typically larger densities and lower temperatures of the disc. If the accretion mode of these BHs were similar to the one in the Galactic Centre, the bright end of X-ray luminosity would barely reach the $\sim 10^{34}$ erg/s level, several orders of magnitude below the Eddington limit for a solar mass BH.

In our model in which local IMBHs are formed through mergers of smaller BHs during the galaxy assembling process, bigger galaxies are more likely to harbour more massive wandering BHs (compare Figures 3 and 4). Given the scaling of the accretion rate with $M_{\rm BH}^2$, this implies that, the larger the galaxy mass, the more luminous will be the sources associated with accreting IMBHs. This is shown in Figure 11, where the bolometric luminosity function is computed for a galaxy of $M=10^{13}M_\odot$. Both a spiral and an elliptical galaxy models are considered. If the accretion efficiency is high ($\eta\sim 0.1$), a spiral galaxy has a 10% probability of hosting a source with $L_x\ga 10^{40}$ erg/s. This is however not the case for an elliptical galaxy, where the lower ISM densities suppress the accretion emission. No ULXs due to IMBHs are expected even in a very massive elliptical galaxy and for a high accretion efficiency.

## "Dressed" wandering BHs

All the previous results have been produced under the assumption that the BHs left over by minor mergers have been completely stripped of their envelope and accrete from their surrounding ISM (*naked* BHs scenario). However, it is possible that, during the merger process, BHs are not completely stripped of their surrounding core of material. Following tidal stripping, a bound core can remain, enveloping the central BH. Mayer & Wadsley (2003) have performed high-resolution N-Body/SPH simulations that follow the evolution of small disc satellites falling into the potential of a much larger halo. Tidal shocks remove most of the outer dark and baryonic mass. In addition, ram pressure affects the gaseous component, stripping it even more severely. After a Hubble time, about $30\%$ of the stars have been stripped, and $80-90$% of the gas disc has been removed, leaving a system resembling a gas-poor dwarf spheroidal. A reasonable value for the gas density in satellite remnants, is the typical gas density $\rho\sim 10^{-3}{\rm cm}^{-3}$, found in dwarf galaxies of the local group (Grebel, Gallagher & Harbeck 2003). The density scatter is large, varying by about 3 orders of magnitude from one system to another. An optimistic upper limit for the gas densities was adopted by Islam et al. (2003b, c). Their estimate assumes that all the baryons in the satellite remnants are under the form of gas, and that the baryon content equals the cosmological baryon fraction. In this case, for a standard $\Lambda$CDM cosmology the gas density in the remnants is of order $\rho\sim 8\times 10^{-2}{\rm cm}^{-3}$.

Given the uncertainties involved in the estimates of the density of these possible remnant cores, we prefer to just give some indicative results within the context of this scenario. In particular, we find that about 10% of the galaxies with mass $\sim 2\times 10^{12}M_\odot$ will have a source with bolometric luminosity of $L_{\rm bol}\ga 1.5\times
10^{39}(\eta/0.1)(\rho/0.1\,{\rm cm}^{-3})$ erg/s for the 3$\sigma$ peaks, and $L_{\rm bol}\ga 2\times 10^{41}(\eta/0.1)(\rho/0.1\,{\rm cm}^{-3})$ for the 3.5$\sigma$ peaks. It should also be noted that, if the BHs do accrete from a remnant core, then the luminosity function of the BHs will be independent of their location, and therefore the most luminous ones are likely to be found in the halo, where their number is the largest.

# Comparison with previous work

Our results are much more pessimistic than previous investigations involving IMBHs remnants of the first stars (Islam, Taylor & Silk 2003a,b,c). Although the hierarchical framework is similar this paper and the Islam et al. papers rely on different assumptions about the abundance of primordial seeds and about the dynamical and mass growth history of MBHs. First, their preferred model assumed seed holes, of $m_\bullet= 260\,\,{\rm M_\odot}$, in 3-$\sigma$ peaks at $z=24.6$, leading to a number of seeds almost 3 orders of magnitude larger than that in our reference model (3.5-$\sigma$ peaks at $z=20$). Second, their merger efficiency is much larger, since they assume that any MBH coming within $1\%$ of the virial radius of the host halo actually merges with the central BH. Finally, they ignore the mass growth of MBHs due to accretion of gas, in a luminous phase as quasars, thus requiring a large number of MBH mergers in order to reproduce the observed local $m_{\rm BH}-\sigma_*$ relation. In our model, instead, it is gas accretion, following galactic major mergers, which determines the mass of SMBHs today, with mergers playing a secondary role, thus requiring a much smaller number of seeds to reproduce the observed SMBHs masses.

As a result, our predictions for the number of X-ray sources in a galaxy of given mass generally fall below the predictions of Islam et al. (2003c), when similar estimates for the density and temperature of the core remnant are adopted[^3].

# Summary and discussion

In this paper we have studied the dynamical evolution with cosmic time of a population of IMBHs wandering in galactic halos, in the context of popular hierarchical structure formation theories. We follow the merger history of dark matter halos and associated SMBHs via cosmological Monte Carlo realizations of the merger hierarchy from the end of the dark ages to the present in a $\Lambda$CDM cosmology. Massive black holes form as the end-product of the first generation of stars and then evolve in a hierarchical fashion following the fate of their host halos. Along the merger history of dark matter halos, the BHs embedded in their inner regions sink to the centre of coalescing larger and larger galactic structures owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, and form a binary system. Stellar dynamical processes drive the binary to harden and eventually coalesce. The dynamical evolution of SMBH binaries may be disturbed by a third incoming BH, if another major merger takes place before the pre-existing binary has had time to coalesce. The three BHs are likely to undergo a complicated resonance scattering interaction, leading to the recoil of the binary and of the single BH, thus creating a population of wandering BHs. Other coalescing binary BHs become wandering due to the non-zero net linear momentum carried away by gravitational waves if a binary is asymmetric. BHs retained within the potential well of the host halo are then bound to slowly travel back towards the galactic centre under the effect of dynamical friction. The population of wandering BHs comprises also another class of holes, those originally hosted by halos which are small compared to the primary halo. In this case the wandering holes can, or not, be still surrounded by the remnants of their old host.

Our simulations have allowed us to estimate the number of black holes, their mass distribution and position within the galaxy. The total number of wandering IMBHs correlates with the mass of the galaxy at $z=0$, and the BHs have typically larger masses. As a result, the number of point sources due to IMBHs is also expected to correlate with the mass of the galaxy, and so the luminosity of the brightest sources. To zeroth order, no correlation is expected in our scenario of IMBHs formation and the star formation rate (SFR) in the galaxy at $z=0$, as the mass function and positions of the BHs within the galaxy are essentially determined by its past merger history. However, if a galaxy has a high SFR, then it is likely to also have a larger filling fraction of dense molecular clouds, which would significantly enhance the accretion rate if a BH happened to wander inside them. In the scenario recently proposed by Krolik (2004), BHs can be formed from primordial abundance gas or from stellar mergers [^4], and these mechanisms would make them to preferentially populate the disk of galaxies, therefore making them brighter, especially if they move for most of the time inside dense regions.

In our scenario of BH formation in mini-halos and growth through subsequent galaxy merging, we find that, if the wandering IMBHs do not carry any substantial baryonic core with them, then a Milky Way sized galaxy is not likely to harbour sources with $L\ga 10^{39}$ erg/s. On the other hand, about 10% of spiral galaxies with mass $\sim 10^{13}M_\odot$ would likely have a source with luminosity $\ga
10^{40}$ erg/s. The brightest sources are expected to be found in the disk, due to its typical lower temperature and larger density than the bulge and the halo. No ULXs are expected in ellipticals, even of large mass ($\sim 10^{13}M_\odot$), due to their average lower density and the lack of a cooler, dense disk which boosts the emission. At sub-Eddington luminosities, the luminosity function of the naked IMBH population in a typical, MW-size galaxy is well below the observed point source population in *Chandra* surveys (Colbert et al. 2004). This is not surprising, as this point source population in the sub-Eddington range is dominated by high mass X-ray binaries (Grimm et al. 2004).

Our results imply that, in order to produce ULXs from IMBHs in MW-type galaxies (and ellipticals for that matter), it is necessary that the wandering IMBHs are not completely stripped of their baryonic core, and this remnant material feeds the accretion at a rate higher than what would be allowed from the ISM alone. Even allowing an high efficiency of $\sim 0.1$, a density $\ga 0.01 {\rm cm}^{-3}$ is required in order for a few percent of MW-type galaxies to host a source with luminosity $\ga 10^{40}$ erg/s. However, the density cannot be significantly higher than that. In our models, the high tail of the mass distribution of IMBHs extends to $\sim 10^6-10^8 M_\odot$, and the observational lack of sources with luminosities $\ga 10^{44}$ erg/s implies that these IMBHs must accrete at sub-Eddington rates.

In this paper, however, we have investigated the luminosity that IMBHs can attain in normal, quiescent galaxies. Starbursting galaxies, such as the Antennae (Zezas & Fabbiano 2002) or the Cartwheel galaxy (Gao et al 2003) show much larger population of ULXs compared to quiescent galaxies of the same size. In particular, in the Cartwheel galaxy these sources are found in the outer star-forming ring, believed to be created by radially expanding density waves caused by a plunging merger. Note that all these sources are at a distance larger than $\sim 10$ kpc from the centre of the galaxy, so in our scheme for a quiescent spiral, they would be found in the low-density, gas-poor galactic halo. The density waves expanding outward, though, provide a temporary high density environment that may supply fuel for an IMBH to accrete and shine at high luminosities, higher also than those predicted for IMBHs travelling in galactic discs.

# Acknowledgments

We have benefited from discussions with M. Colpi, T. J. Cox, P. Madau and C. Nipoti.

[^1]: E-mail: email@example.com (MV); email@example.com (RP)

[^2]: The SMBH in the Galactic Centre has an X-ray efficiency of $\sim 4\times 10^{-8}-2\times 10^{-6}$ in the 2-10 keV band (Baganoff et al. 2001). For BHs at the other end of the mass spectrum (i.e. the stellar-mass BHs), Agol & Kamionkowski (2002) used a best guess estimate of $\eta_{\rm X}\sim
    10^{-5}$.

[^3]: Note that we consider the possibility of BHs accreting from a baryonic remnant only for the satellites, whose mass distribution has an upper end cutoff which is smaller than that of the sling BHs (see Fig. 3 and 4).

[^4]: For studies of stellar merging see also Portegies Zwart & McMillan (2002) and Gurkan, Rasio & Freitag (2004); note that the formation of IMBHs by runaway stellar merging is possible only in young, compact clusters.


abstract: Recent advances in microelectronics have enabled the realization of Wireless Body Area Networks (WBANs). Increasing the transmission power of WBAN's nodes improves the Signal to Interference plus Noise Ratio (SINR), and hence decreases the bit error probability. However, this increase may impose interference on nodes within the same WBAN or on other nodes of nearby coexisting WBANs, as these WBANs may use similar frequencies. Due to co-channel interference, packet collisions and retransmissions are increased and consequently, the power consumption of the individual WBANs may increase correspondingly. To address this problem, we adopt the approach of two-hop cooperative communication due to its efficiency in power savings. In this paper, we propose a cooperative power control-based algorithm, namely, IMA, for interference mitigation among the individual sensors of a single WBAN. Basically, our approach selects an optimal set of relays from the nodes within each WBAN to mitigate the interference. Thus, IMA selection criterion relies on the best channel, namely, SINR and power conditions to select the set of best relays. The experimental results illustrate that IMA improves the SINR, the power efficiency and extends WBAN lifetime. In addition, the results illustrate that IMA lowers the bit error probability and improves the throughput.
author:

Email: {mohamad.ali; hassine.moungla; email@example.com
title: Energy Aware Competitiveness Power Control in Relay-Assisted Interference Body Networks

# Introduction

The pervasive use of wireless networks, recent developments in micro-electronics and the miniaturization of low-power sensors led to the existence of WBANs. The technological advancements in the wireless communication enabled low-power, intelligent, miniaturized sensor nodes placed in or around the human body to enhance the quality of patient's life. Various applications such as personal health monitoring, ubiquitous healthcare, sports, entertainment, and military have been found to provide high reliable data communication systems as well as to improve the healthcare conditions. Their individual sensors mainly monitor physical activities, actions and vital signs as glucose percentage in blood, heart beats, respiration, body temperature and/or can record electrocardiography (ECG) .

The co-channel interference is challenging due to the highly mobile and resource constrained nature of WBANs. Firstly, such nature makes the allocation of a global coordinator to manage the coexistence problem among coexisting WBANs as well as the application of advanced antenna and power control techniques used in other networks unsuitable for WBANs. Secondly, the stringent resource in WBANs is the energy, which requires a careful calibration of the transmission power of the individual nodes to minimize interference in the network. Thirdly, due to the infeasibility of coordination among the WBANs and the unpredictable movements of their individual sensors, an interference may arise due to the simultaneous transmissions of sensors of different coexisting WBANs using the same channel. For instance, an interference may happen when the simultaneous reception from two or more distinct sensors of different WBANs collide at the same receiving node such as a WBAN's coordinator. Consequently, the interference may affect the communication links and degrade the performance of each individual WBAN . Therefore, interference mitigation is of the utmost importance to improve the reliability of the whole network. To this end, the IEEE standard proposes three mechanisms for co-channel interference mitigation in WBANs, namely, beacon shifting, channel hopping and active superframe interleaving . In this paper, we tackle these issues and contribute the following:

-   *IMA, a cooperative power control-based scheme for interference mitigation among the individual sensors of a single WBAN*

-   *A relay selection criterion that determines an optimal set of relays using the best-known channel and node energy conditions*

The simulation results show that our proposed approach can significantly lower the interference among the individual sensors of an intra-WBAN as well as increase the power savings at both node- and WBAN- levels. Moreover, IMA significantly avoids the intra-WBAN interference and does not require any mutual coordination among the individual WBAN coordinators. The rest of the paper is organized as follows. Section @ sets our work apart from other approaches in the literature. Section @ summarizes the system model and the assumptions. Section @ describes *IMA* in detail. Section @ presents the simulation results. Finally, the paper is concluded in Section @

# Related Work

The problem of intra-WBAN interference has been addressed through cooperation, power control, link adaption, multiple access and resource allocation schemes. Example schemes that pursued the cooperation methodology include , , . Dong et al., proposed a single-relay cooperative scheme where the best relays are eventually selected in a distributed fashion. The scheme is based on MAC request-to-send (RTS) and clear-to-send (CTS) signaling, where a set of potential relays compute individually the required transmission power to participate in the cooperative communication to significantly mitigate the interference. Also, Dong et al., addressed the problem of coexistence of multiple non-coordinated WBANs. A decode-and-forward protocol with two relays and selection combining at the desired coordinator is used, enabling intra-network and inter-network operation, to allocate slots for each link packet transmission to mitigate intra-WBAN interference. Whilst, Feng et al., proposed a prediction-based dynamic relay transmission scheme through which the problem of "when to relay" and "who to relay" is decided in an optimal way. Other approaches pursued power control schemes for intra-WBAN interference mitigation. Dong et al., proposed a joint two-hop relay-assisted cooperative communication integrated with transmit power control for intra-WBAN interference mitigation. This scheme reduces co-channel interference at the WBAN individual nodes and extends the WBAN energy lifetime.

A number of approaches adopted medium access schemes to mitigate intra-WBAN interference include , . Mahapatro et al., proposed a TDMA-based scheme that enables two or three coexisting WBANs to agree on a common TDMA schedule to reduce the intra-WBAN interference and the transmission latency per each individual node. Whereas, Chen et al., proposed a distributed scheme that adopts a polling-based TDMA for traffic coordination of the individual sensors within a WBAN and a carrier sensing mechanism conducted before each beacon transmission to deal with inter-WBAN interference. Meantime, some approaches adopted the link adaption methodology to mitigate the intra-WBAN interference include , , . Yang et al., focused on the performance at the WBAN coordinator which periodically calculates SINR. Based on this calculation, it commands the nodes within its WBAN to select an appropriate scheme (data rate, modulation, duty cycles, etc.) to mitigate the interference. Meantime, Martelli et al., considered a WBAN where the coordinator periodically queries its sensors to transmit data. The network adopts the carrier sense multiple access with collision avoidance (CSMA/CA) and the nodes adopt link adaptation to select the modulation scheme according to the experienced level of interference. Whilst, Domenicali et al., analyzed the performance of a reference WBAN in terms of bit error rate, throughput, and energy lifetime. The study proved the performance of a WBAN can be improved by the adoption of an optimized time hopping code assignment strategy. A strategy to extend the lifetime of the WBAN is also introduced.

Jamthe et al., pursued the medium access methodology and proposed a quality of service based medium access control (MAC) scheduling approach to avoid inter-WBAN interference and introduce a fuzzy inference engine for intra-WBAN scheduling so as to avoid interference within WBANs. Other approaches pursued the resource allocation include , . Liang et al., proposed a distributed interference detection and mitigation scheme through using adaptive channel hopping for intra- and inter-WBAN interference mitigation. Whilst, Movassaghi et al., proposed a dynamic resource allocation scheme for intra- and inter- WBAN interference mitigation among multiple coexisting WBANs through using orthogonal channels for high interfering nodes.

Though, most of the recent works addressed problems related to co-channel interference at the WBAN-level and do not consider the node-level interference. In this paper, we take a step forward, and exploit the benefits of two-hop cooperation, and propose a cooperative power control-based algorithm for interference mitigation among the individual sensors of a single WBAN. Thus, we depend on the relaying and two-hop communication scheme to minimize intra-WBAN interference and increase the power savings at the node- and WBAN-levels.

# System Model and Assumptions

We consider the realistic scenario when multiple beacon-enabled CSMA/CA-based WBANs coexist in a large hall of a hospital. Each WBAN consists of a single coordinator and up to K sensors, each generates its data based on a predefined sampling rate and transmits data at a maximum rate of 250Kb/s within the license-free 2.4 GHz band. Unlike sensors, we assume all coordinators are equipped with unconstrained energy supply. In fact, the time period between any two consecutive beacons depends on the application, the required bit rate, and the transmitted packet size. The smaller the period is, the smaller data packet size is, the more frequent contention for the channel and the smaller the throughput is. Thus, we assume the followings:

-   A WBAN topology based on two-hop communication

-   A CSMA/CA as a medium access scheme

# Intra-WBAN Interference Mitigation Approach - IMA

As pointed out, a co-channel interference may arise due to the collisions amongst the concurrent transmissions made by sensors in different WBANs using the same channel. To address this issue, we exploit the benefits of the two-hop relaying schemes in order to tackle this problem and reduce the probability of collision while enabling autonomous scheduling of the medium access within each WBAN. Therefore, IMA minimizes the intra-WBAN interference and maximizes the energy lifetime. The coordinator of each WBAN periodically broadcasts beacons, through which the nodes synchronize their transmissions according to its clock. Then, after the successful reception of each beacon, the nodes within each WBAN compete for the channel using CSMA/CA. During the contention access period (CAP), each WBAN's node pursues the communication model with its associated coordinator as follows, RTS-CTS-DATA-Acknowledgment (ACK). Received beacons and RTS packets are used in SINR computation and remain valid as long as they happen to be within the coherence time of the channel. IMA involves four phases; Beacon phase, RTS phase, CTS phase and $CTS_{Reception}$ phase.

## Beacon phase

In this phase, the coordinator periodically broadcasts beacons every T milliseconds (ms) to enable the sensors within its WBAN to synchronize. Based on the power contained within the beacon signal, each node computes the corresponding SINR which is defined at a given node as the power received from the desired transmitter divided by the sum of undesired powers received from other interfering nodes plus the power contained within the noise. Basically, it is used to evaluate the channel quality along with any path from a source to the coordinator through relays. Thus, R denotes the set of WBAN's nodes, where each node can successfully decode the beacon and whose *$SINR >= SINR_{Thr}$*. SINR is computed as in **eq.(1)**. Where P is the desired power at receiver, $I_{i}$ is the interference power received from interfering node *i* and $N_{0}$ is an additive white Gaussian noise $$SINR =\cfrac{P}{\sum_{i=1}^{N} I_i + N_0}$$

## RTS phase

After the successful reception of a beacon, All nodes that have pending data (denoted by sources) start a contention for the channel, where each node transmits its short packet ($RTS_{i}$). The transmissions of such RTS packets may introduce interference at some other nodes. Each node that is in the reception mode when sensing these RTS packets, it immediately computes the SINR of the received signal. If SINR is above a threshold, i.e., *$SINR >= SINR_{Thr}$*, this node will be included in the relay candidate set. Thus, each relay candidate has successfully received a beacon and RTS packet can compete for the relay selection process in the next phase. Therefore, the periodic computations evaluate the channel quality from each relay candidate to its corresponding source node. Up to this information, each relay candidate can evaluate the communication path from a particular source node to the coordinator through itself. We define:

-   S denotes the set of all source nodes have pending data to transmit in the current superframe

-   $N_{i}$ denotes the set of neighbors of a source node denoted by $S_{i}$ that successfully decode $RTS_{i}$

## CTS phase

Each relay candidate has the knowledge of the channel condition from itself to the source nodes and to the coordinator. Basically, all nodes that have low SINR values are excluded from the relay selection process to avoid the interference. The main goal is to select a relay for each source node so that the interference is minimized and the power saving is maximized. This way, each relay sets a timer to compute a waiting time value. Then, that relay computes the SNR and SINR values from the signals received from the coordinator and the source node, respectively. Consequently, the smaller the difference between both values (SNR and SINR) at each relay implies the better the candidate is likely to be chosen as the best relay. knowing that the waiting time is inversely proportional to their values difference. This difference follows a threshold margin, hence, each relay candidate within the set $Q_{i}$ of $S_{i}$ transmits its CTS to $S_{i}$, where $Q_{i}$ denotes the set of nodes successfully decode the $RTS_{i}$ and the beacon, $Q_{i} = R \cap N_{i}$. Then we define the beacon validity as:

**if (reception time of RTS - reception time of beacon) $<$ T ms, then it is valid, otherwise, it is not**

## $CTS_{RECEPTION}$ phase

Upon the reception of CTS packets from the relay candidates, each source node waits a period of time and enqueues the last two CTS packets received from the last two distinct relays. Then, it selects the best relay based on the high SINR and residual energy (ER) of the battery reported within each CTS packet. M denotes the set of nodes successfully decode the CTS. **Algorithm ** provides a summary of IMA.

# Performance Evaluation

## Simulation Environment and Setup

The simulation environment is based on the IEEE CSMA/CA specifications, where a small scale fading channel model and a channel coherence time of 500ms are considered for on-body communication . The WBAN is located within a space of 3x3x3$m^{3}$ and consists of a coordinator and up to M=8 sensor nodes. We also used the industrial, scientific and medical (ISM) radio bands for operating frequency of 2.4 GHz. We compared the performance of a WBAN that employs IMA coexisting with other WBANs that do not employ IMA. To have further fair results and to illustrate the outperformance of our algorithm, we generate the same simulation setup in terms of topology, the channel model, the medium access control (MAC) and the physical (PHY) specifications which involve the node transmission power (-10 dBm), receiver sensitivity (-84.7 dBm), noise floor (-102 dBm), an additive white Gaussian noise (AWGN) of mean zero and $\sigma$=6.81 dBm as well as path loss exponent (4.22) . In this work, we focus on three performance metrics; SINR, lifetime and throughput, which have the direct impact on the performance of each individual WBAN.

1.  **WBAN Lifetime** is defined as the sum of residual energies of WBAN nodes at a particular point in time.

2.  **Battery energy residue** is defined as the amount of energy stored within a particular node's battery at a particular point in time.

3.  **Throughput** is defined as the sum of the number of successful packets delivered in a period of time at a particular node.

4.  **Outage probability (OP)** is defined as the probability of SINR being below a given threshold.

$$OP = Pr(SINR \leq SINR_{threshold})$$

## Network lifetime and battery energy residue

All the nodes of a WBAN have the same initial battery energy (150 milliJoules). The energy residue of IMA scheme denoted by RE and that of the competing scheme, namely, without-IM, are compared in **Figure **. As can be seen from the figure, RE of IMA is always higher than that of the competing scheme all the time. Referring to curves in the figure, it is evident to notice that RE values of both schemes decrease linearly with time. The AVGRE results of with-IMA are always larger than the AVGRE results of without-IMA. with-IMA's RE starts at 150 milliJoules and eventually stabilizes at 150. However, without-IMA's RE decreases very slightly to eventually stabilize to 145 milliJoules. Whilst, in without-IMA, the AVGRE starts at 150 milliJoules and decreases faster than with-IMA to stabilize at 130 milliJoules, while the second stabilizes at 120. The improvement in energy savings and the extension of WBAN lifetime is due to some reasons. 1) IMA selects relay nodes of minimal interferences that were supposed to share in the RTS/CTS contention phase. Such relays are now avoided with the employment of IMA algorithm and so they do not transmit their CTS packets accordingly. Moreover, these relays do not receive their respective ACKs as well as do not introduce interferences at other nodes. Consequently, the number of communicated packets in the system as well as the probability of collisions and hence energy consumption. Also, with IMA employment, no need for nodes that receive only RTS to reply with their CTS which decreases packet transmissions are reduced which lowers the energy consumption. 2) IMA selects nodes with higher SINR values which improve the energy per bit, which increases the probability of a successful delivery of packets at some nodes and hence decreases the number of retries for retransmissions.

## Throughput

The sum of successful packets received (SoP) at the WBAN coordinator versus time for a WBAN that employs IMA and another that does not employ IMA are compared in **Figure **. As can be clearly seen in the figure, IMA always provides a higher SoP than the competing scheme (another scheme that does not employ IMA) all the time. Referring to the curves in the figure, for instance, at time=0, both results have similar SoP values and hence, as time increases, SoP of both schemes increase linearly in different speeds, i.e, the curve of IMA is faster than that of the another scheme. The difference between their SoPs curves becomes evident at time=3000 seconds, SoP=300 packets successfully received at the WBAN coordinator that does not employ IMA, whereas, SoP=630 packets received at the WBAN coordinator that employs IMA. The improvement of the throughput in IMA is due to the increase in the SINR and the avoidance of interfering nodes. Hence, the achievement of higher SINR increases the energy per bit; the signal becomes less error-prone and more resistive to interferences, also the fraction of bit errors is decreased and which decreases the packet retransmissions as well. Furthermore, the throughput is increased due to the decrease in the number of collisions and packet retransmissions. In addition, IMA algorithm selects best relays with minimal interference level which increases the chances for successful packets delivery.

## Discussion of RTS and CTS

IMA suggests that any node receives both beacon and RTS packets with a precise SINR difference margin contends for the relay selection process. Firstly, these nodes that receive RTS packets only are excluded from the selection process and hence can not reply by their CTS packets which decrease the number of CTS packets. Secondly, IMA avoids nodes that experience high interference (with minimal SINR) even though they overhear both the coordinator (beacons) and the source nodes (RTS packets). These nodes are also excluded from the selection process and hence they do not transmit their CTS packets. Consequently, this lowers the number of nodes compete for the selection process as well as decreases the number of exchanged CTS and their corresponding ACKs. As a result, the total number of RTS/CTS and ACKs are decreased in the WBAN.

# Conclusion

In this paper, we have presented *IMA*, a distributed CSMA-based power control-based for intra-WBAN interference mitigation scheme based on two-hop cooperative communication. *IMA* exploits the benefits of the two-hop scheme to lower the probability of collisions among transmission of sensors within a single WBAN as well as among the different coexisting *WBAN*s. Accordingly, each distinct sensor within the WBAN selects the best relay to retransmits its packet to the coordinator. Compared with other sample algorithm, *IMA* has low complexity and does not require any inter-*WBAN* coordination. Simulation results show that *IMA* outperforms other sample schemes in terms of interference mitigation, energy savings, and throughput.


author:
II. Physikalisches Institut
University of Goettingen
Friedrich-Hund-Platz 1
37077 Göttingen
E-mail:
title: The Electroweak Sudakov approximation in SHERPA

# Introduction

Monte Carlo event generators provide an important bridge between collider experiments and the underlying physics of the SM. They enable low-multiplicity fixed-order calculations to be compared against the busy hadronic environment that characterises modern-day experimental particle physics. This is achieved by dressing the fixed-order calculation with a parton shower and introducing underlying events and hadronisation effects. The current state-of-the art for these predictions involves matching NLO QCD matrix elements to a parton shower. Moving beyond this perturbative accuracy requires either performing an NNLO QCD or an NLO EW calculation, both of which require new technology and often contribute a similar order of magnitude to the correction. There has been a lot of work on both calculating fixed order NNLO QCD corrections such as  and on the matching of NNLO QCD cross sections with a parton shower . However, the work presented in the following considers the extension of Monte Carlo event generators, specifically the SHERPA  event generator, to include NLO EW corrections. There is a particular focus on the EW Sudakov approximation, the origin of which is outlined and the implementation within SHERPA briefly explained. Some initial results are presented for a 14 TeV LHC.

# NLO QCD

To begin a discussion on NLO EW, it is instructive to consider the comparative case of NLO QCD. Calculating the NLO QCD corrections to a process analytically is given by,

$$\label{EQN:NLOana}
  \sigma_{\text{NLO}}^{\text{QCD}} = \int(\text{B}+ \text{V})\text{d}\Phi_\text{B}+ \int\text{R}\text{d}\Phi_\text{R}\,,$$

which cannot be easily transferred to a numerical calculation due to the divergence of the virtual ($\text{V}$) and real ($\text{R}$) integrals on the RHS. The born contribution, $\text{B}$, obviously does not contain any such divergences. Each term in eqn. () is integrated over its appropriate phase-space to give the full NLO cross section, $\sigma_{\text{NLO}}^{\text{QCD}}$. The solution to this problem which is almost universally adopted is to use a subtraction scheme . This subtracts a quantity, $\text{S}$ which exactly matches the divergent behaviour of the integral over the real emission. This must be analytically integrable over the 1 parton sub-space and not introduce any new divergences . Once these conditions are met, the subtraction term can be added to the virtual integral. Once it is appropriately integrated it will then also match exactly the divergent structure of the virtual contribution. With this introduction, eqn () becomes a sum of finite integrals,

$$\sigma_\text{NLO}^{\text{QCD}} = \int(\text{B}+\text{V}+\int\text{S}\text{d}\Phi_1)\text{d}\Phi_\text{B}+\int(\text{R}-\text{S})\text{d}\Phi_\text{R}\,.$$

# NLO EW

## Performing the calculation

There are several differences between an NLO QCD calculation and an NLO EW calculation. One obvious difference is the appearance of masses of the EW gauge bosons, which regulate the divergences from soft and collinear radiation. The real emission of $W^\pm$ and $Z$ bosons decay into other particles and are, theoretically, distinguishable processes from the underlying Born term. This reduces the problem of including real radiation to that of QED, which can be treated in a similar way to the NLO QCD calculation above. Because the real radiation can be, at least to a large extent, classified as a distinct process, the large corrections in the high energy regime, where the mass of the weak bosons becomes negligible, are physically meaningful. This high energy regime is the limit where the Sudakov logarithmic approximation becomes valid.

Another difference introduced in NLO EW calculations is the dependence of the EW bosons on the helicity of the particles involved, particularly clear in $W^\pm$ boson emission. Unlike in QCD, the couplings of the weak bosons are strongly dependent of the helicity of the particle. Furthermore, the exchange of weak particles can change the underlying Born term, for example mixing electrons with neutrinos, and therefore create interferences between previously distinct processes. These differences affect both the NLO EW calculation and the implementation of the EW Sudakov logarithmic approximation.

## EW Sudakov Approximation

A full NLO EW calculation is very computationally intensive, whereas the EW Sudakov approximation does not introduce much overhead and is therefore comparatively cheap. It is also easier to include on top of NLO QCD calculations. This approximation considers only the logarithmic contribution to the correction, and is dominant in the high-energy regime. The EW Sudakov approximation in SHERPA follows an excellent and clear paper by Denner and Pozzorini , which provides a break down of how to implement Sudakov logarithms in a process-independent way. These logarithms, $L$, typically take the form

$$\label{EQN:sud1}
L\sim\log\left(\frac{s}{M_V^2}\right)\,,$$

which captures the high energy behaviour of the NLO EW calculations. In eqn (), $s$ is the centre of mass energy, and $M_V$ is the mass of the relevant weak boson. There is always a non-logarithmic piece which is neglected, and this creates a limit to the accuracy EW Sudakov logarithms can reach, typically $\sim 1\%$. These logarithms are produced by the exchange of soft-collinear weak bosons, which become divergent in the limit of vanishing mass, and are referred to as mass singular diagrams, as discussed in ref. . These diagrams involve either the exchange of an EW boson between 2 external legs (double logarithms) or the emission of a soft or collinear boson from an external leg (single logarithms). These contributing diagrams are shown in fig. . The right-hand diagram in fig.  includes the wave-function renormalisation terms where the boson is reabsorbed by the emitting line. There are also single logarithms from parameter renormalisation, which are not depicted.

It is clear that the double logarithmic terms are the dominant contribution, with the single logarithms providing a subleading contribution. However, while the double logarithms are typically a negative correction to the Born cross section, the subleading contribution is a positive correction. Considering only the leading contribution therefore leads to an overestimation of the size of the Sudakov correction.

In the one-loop approximation, the EW Sudakov logarithms introduce corrections, summed over external legs $i$ and $j$ and EW bosons with typical mass $M$,

$$\label{EQN:sudfin}
   L=\frac{\alpha}{4\pi}\left[A\log^{2}\left(\frac{(p_i+p_j)^2}{M^2}\right)
   +B\log\left(\frac{(p_i+p_j)^2}{M^2}\right)\right]\,,$$

where $p_i$ denotes the momentum of leg $i$.

## EW Sudakov results

Within the SHERPA framework, the EW Sudakov approximation is included as a K-factor that is applied at each phase-space point. As is implied by eqn (), it depends only on the final state legs, and iterates over all possible exchanges of EW bosons. All bosons are assumed to have a mass equal to the $W^\pm$ boson mass. This introduces only a small logarithmic correction for the $Z$ boson, which can be neglected to the order considered here. The mass difference between the photon and the $W^\pm$ boson introduces large logarithms, but these largely cancel against real photon radiation. Therefore, the assumption that the EW bosons all have equal mass does not have a significant impact on the correction. The EW Sudakov approximation only affects the hard process and can be easily employed with the parton shower or applied to an NLO QCD computation. However, the implementation does rely on the COMIX  matrix element generator.

The results shown here are the first results with the implementation of NLO EW Sudakov logarithms within SHERPA. The calculations are performed for a 14 TeV LHC collider. The left-hand side of fig.  shows the effect of including NLO EW Sudakov corrections at 14 TeV in off-shell $W^\pm$+jets production. It is clear that as the $p_T$ of the leading jet increases, the relative correction from the EW Sudakov logarithms becomes larger and more negative. This reaches almost $40\%$ at 1 TeV. The behaviour is similar on the right-hand side of fig. , throughout the $p_T$ spectrum, which shows the same distribution but for on-shell production of the $W^\pm$ boson.

## Full NLO EW

Although the NLO EW Sudakov approximation is comparatively quick and easy to include on top of NLO QCD computations, for some studies, either to a much higher precision or outside of the high energy regime, a full NLO EW computation must be employed. This faces the same challenges and subtleties that the NLO EW Sudakov approach dealt with, alongside ambiguities in process definition. It must be decided, for example, what counts as a photon emitted in the NLO EW calculation and what is simply radiation from a jet. Also, the difference between NLO QCD corrections to EW processes and NLO EW corrections to QCD processes must be defined in order to avoid double counting.

Within the SHERPA event generator, there is currently a working interface to the OpenLoops  loop provider for the NLO EW virtual amplitude and QED subtraction handled within SHERPA. There are already publications on NLO EW corrections, including multijet merging, to $V$+jets , both with an on-shell $W^\pm$ boson and including off-shell effects, however the code is not yet public. There is also ongoing effort to implement an NLO EW interface between SHERPA and Recola , as another one-loop provider for both NLO QCD and NLO EW virtual amplitudes.

# Conclusions

Improving the perturbative accuracy of the hard interaction in Monte Carlo event simulation now involves the calculation of NLO EW contributions, which are often of a comparable size to NNLO QCD. This includes several new challenges, and the full computation is quite time-intensive. EW Sudakov logarithms provide a simple way for the dominant behaviour of the NLO EW calculation to be taken into account without the computational overhead. It is also trivial to include on top of QCD corrections, unlike the full calculation where the interference terms must be carefully considered to avoid double counting. This is an important complementary approach to be implemented alongside the full calculation. There is also a lot of promising progress in the automated evaluation of the full NLO EW correction.


abstract: Our aim is to test potential solar prominence plasma diagnostics as obtained with the new solar capability of the *Atacama Large Millimeter / submillimeter Array* (ALMA). We investigate the thermal and plasma diagnostic potential of ALMA for solar prominences through the computation of brightness temperatures at ALMA wavelengths. The brightness temperature, for a chosen line of sight, is calculated using densities of hydrogen and helium obtained from a radiative transfer code under non local thermodynamic equilibrium (NLTE) conditions, as well as the input internal parameters of the prominence model in consideration. Two distinct sets of prominence models were used: isothermal-isobaric fine-structure threads, and large-scale structures with radially increasing temperature distributions representing the prominence-to-corona transition region. We compute brightness temperatures over the range of wavelengths in which ALMA is capable of observing (0.32 – 9.6mm), however we particularly focus on the bands available to solar observers in ALMA cycles 4 and 5, namely 2.6 – 3.6mm (Band 3) and 1.1 – 1.4mm (Band 6). We show how the computed brightness temperatures and optical thicknesses in our models vary with the plasma parameters (temperature and pressure) and the wavelength of observation. We then study how ALMA observables such as the ratio of brightness temperatures at two frequencies can be used to estimate the optical thickness and the emission measure for isothermal and non-isothermal prominences. From this study we conclude that, for both sets of models, ALMA presents a strong thermal diagnostic capability, provided that the interpretation of observations is supported by the use of non-LTE simulation results.
address: SUPA, School of Physics & Astronomy, University of Glasgow, Glasgow G12 8QQ, Scotland
author:   ;
bibliography: paper.bib
title: Solar prominence modelling and plasma diagnostics at ALMA wavelengths

```latex
\begin{article}

\begin{opening}

\title{Solar prominence modelling and plasma diagnostics at ALMA wavelengths}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Authors Names
%
\author[addressref={gla},corref,email@example.com}]{\inits{A.S.}\fnm{Andrew}~\lnm{Rodger} \orcid{0000-0003-0385-4581}}

\author[addressref={gla},email@example.com}]{\inits{N.}\fnm{Nicolas}~\lnm{Labrosse} \orcid{0000-0002-4638-157X}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Runningheads
%
%\runningauthor{}
%\runningtitle{}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Affilations
%% id shold be the same with \author addressref value.
\address[id={gla}]{SUPA, School of Physics \& Astronomy, University of Glasgow, Glasgow G12 8QQ, Scotland}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Abstract
\begin{abstract}
Our aim is to test potential solar prominence plasma diagnostics as obtained with the new solar capability of the \emph{Atacama Large Millimeter / submillimeter Array} (ALMA).
   We investigate the thermal and plasma diagnostic potential of ALMA for solar prominences through the computation of brightness temperatures at ALMA wavelengths.
   The brightness temperature, for a chosen line of sight, is calculated using densities of hydrogen and helium obtained from a radiative transfer code under non local thermodynamic equilibrium (NLTE) conditions, as well as the input internal parameters of the prominence model in consideration.
   Two distinct sets of prominence models were used: isothermal-isobaric fine-structure threads, and large-scale structures with radially increasing temperature distributions representing the prominence-to-corona transition region.
   We compute brightness temperatures over the range of wavelengths in which ALMA is capable of observing (0.32 -- 9.6mm), however we particularly focus on the bands available to solar observers in ALMA cycles 4 and 5, namely 2.6 -- 3.6mm (Band 3) and 1.1 -- 1.4mm (Band 6).
   We show how the computed brightness temperatures and optical thicknesses in our models vary with the plasma parameters (temperature and pressure) and the wavelength of observation. We then study how ALMA observables such as the ratio of brightness temperatures at two frequencies can be used to estimate the optical thickness and the emission measure for isothermal and non-isothermal prominences.
   From this study we conclude that, for both sets of models, ALMA  presents a strong thermal diagnostic capability, provided that the interpretation of observations is supported by the use of non-LTE simulation results.
\end{abstract}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Keywords
%
\keywords{Corona, Radio Emission; Prominences, Models; Radio Emission,  Quiet; Spectrum, Continuum}

\end{opening}
%-------------------------------------------------

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Sections
%
% \section{}%\label{s:?}

%% Figure
%
% \begin{figure}
% \centerline{\includegraphics[width=0.5\textwidth,clip=]{<fig.eps>}}
% \caption{}%\label{fig:?}
% \end{figure}

%% Table
%
% \begin{table}
% \caption{}%\label{tbl:?}
% \begin{tabular}{}
% \hline
% \multicolumn{2}{c}{<>}
% <data>
% \hline
% \end{tabular}
% \end{table}

\input{introduction.tex}

\input{modelling.tex}

\input{plasma_diagnostics.tex}

\input{discussion.tex}

\input{conclusions.tex}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Appendix
%
% \appendix

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Acknowledgements
%
\begin{acks}
  AR acknowledges support from a STFC studentship.
  NL acknowledges support from STFC grant  ST/L000741/1.
  This research made use of NASA's Astrophysics Data System.

\end{acks}

%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Bibliography
%
% Using BibTeX
%
\bibliographystyle{spr-mp-sola}
\bibliography{paper}
%
% Without BibTeX
% \begin{thebibliography}{}
% \bibitem[\protect\citeauthoryear{Author}{Year}]{key}
%   <bibliographical entry>
%
% \bibitem[\protect\citeauthoryear{}{}]{}
%
%
% \end{thebibliography}

\end{article}
```