The main SciRate homepage is down (when not logged in). We are working to fix it. See https://github.com/scirate/scirate/issues/337 for updates.

# Astrophysics of Galaxies (astro-ph.GA)

• We use observations of the neutral atomic hydrogen (HI) 21-cm emission line to study the spatial distribution of the HI gas in a 80$\degree\times~$90$\degree$ region of the Galaxy halo. The HI column densities in the range of \mbox3--11$\times$10$^{20}$ cm$^{-2}$ have been estimated for some of the studied regions. In our map---obtained with a spectral sensitivity of $\sim$2 K---we do not detect any HI 21-cm emission line above 2$\sigma$ at Galactic latitudes higher than $\sim$46$\degree$. This report summarizes our contribution presented at the conference on the origin and evolution of barionic Galaxy halos.
• The traditional knowledge of the mechanisms that caused the formation and evolution of early-type galaxies (ETG) in a hierarchical universe was challenged by the unexpected finding by ATLAS3D that 86% of the ETGs show signs of a fast-rotating disk. This implies a common origin of most spiral galaxies, followed by a quenching phase, while only a minority of the most massive systems are slow rotators and were likely to be the products of merger events. Our aim is to improve our knowledge on the content and distribution of ionized hydrogen and their usage to form stars in a representative sample of ETGs for which the kinematics and detailed morphological classification were known from ATLAS3D. Using narrow-band filters centered on the redshifted Halpha line along with a broad-band (r-Gunn) filter to recover the stellar continuum, we observed or collected existing imaging observations for 147 ETGs (including members of the Virgo cluster) that are representative of the whole ATLAS3D survey. Fifty-five ETGs (37%) were detected in the Halpha line above our detection threshold, (Halpha E.W. <= -1 AA), and 21 harbor a strong source (Halpha E.W. <=-5 AA). The strong Halpha emitters appear associated with low-mass (M 10^10 M_odot) S0 galaxies that contain conspicuous stellar and gaseous discs. These harbor significant star formation at their interior, including their nuclei. The weak Halpha emitters are almost one order of magnitude more massive, contain gas-poor discs and harbor an AGN at their centers. Their emissivity is dominated by [NII] and does not imply star formation. The 92 undetected ETGs are gas-free systems that lack a disc and exhibit passive spectra even in their nuclei. These pieces of evidence reinforce the conclusion that the evolution of ETGs followed the secular channel for the less massive systems and the dry merging channel for the most massive galaxies.
• Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.
• Based on CO(2-1) public data, we study the monoxide oxygen gas excitation conditions and the magnetic field strength of four spiral galaxies. For the galaxy outskirts, we found kinetic temperatures in the range of $\lesssim$35--38 K, CO column densities $\lesssim 10^{15}$--$10^{16}$ cm$^{-2}$, and H$_2$ masses $\lesssim 4\times 10^6$--$6\times10^8$ M$_\odot$. An H$_2$ density $\lesssim 10^3$ cm$^{-3}$ is suitable to explain the 2$\sigma$ upper limits of the CO(2-1) line intensity. We constrain the magnetic field strength for our sample of spiral galaxies and their outskirts by using their masses and H$_2$ densities to evaluate a simplified magneto-hydrodynamic equation. Our estimations provide values for the magnetic field strength on the order of $\lesssim$6--31 $\upmu$G.
• The Galactic object MWC 137 was suggested to belong to the group of B[e] supergiants. However, with its large-scale optical bipolar ring nebula and the high velocity jet and knots, it is a rather atypical representative of this class. We performed multi-wavelength observations spreading from the optical to the radio regime. Based on optical imaging and long-slit spectroscopic data we found that the northern parts of the large-scale nebula are predominantly blue-shifted, while the southern regions appear mostly red-shifted. We developed a geometrical model consisting of two double-cones. While various observational features can be approximated with such a scenario, the observed velocity pattern is more complex. Using near-infrared integral-field unit spectroscopy we studied the hot molecular gas in the close vicinity of the star. The emission from the hot CO gas arises in a small-scale disk revolving around the star on Keplerian orbits. While the disk itself cannot be spatially resolved, its emission is reflected by dust arranged in arc-like structures and clumps surrounding MWC 137 on small scales. In the radio regime we mapped the cold molecular gas in the outskirts of the optical nebula. We found that large amounts of cool molecular gas and warm dust embrace the optical nebula in the east, south and west. No cold gas or dust were detected in the north and north-western regions. Despite the new insights on the nebula kinematics gained from our studies, the real formation scenario of the large-scale nebula remains an open issue.
• Shell galaxies are believed to form through the collision of a dwarf galaxy with an elliptical galaxy. Shell structures and kinematics have been noticed to be independent tools to measure the gravitational potential of the shell galaxies. In this work, as it is missing in the literature, we compare theoretically the formation of shells in Type I shell galaxies in different gravity theories, including Newtonian plus dark halo gravity, and two non-Newtonian gravity models, MOG and MOND in identical initial systems. We investigate the effect of dynamical friction, which by slowing down the dwarf galaxy in the dark halo models, limits the range of shell radii to small values. Under the same initial conditions, shells appear on a smaller time-scale and over a smaller range of distances in the presence of dark matter compared to the corresponding non-Newtonian gravity models. If galaxies are embedded in dark matter halo, then the merging time may be too rapid to allow multi-generation shell formation as required by observed systems due to the large dynamical friction effect. Starting from the same initial state, in the dark halo model the observation of small bright shells should be accompanied by large faint ones, while in the absence of dark matter the next shell generation patterns iterate with a specific time delay. The first shell generation pattern shows a degeneracy with the age of the shells and different theories, but the relative distance of the shells and the shell expansion velocity can break this degeneracy.
• The aim of this work is to analyse the radio properties of the massive and dynamical disturbed clusters Abell 1451 and Zwcl 0634.1+4750, especially focusing on the possible presence of diffuse emission. We present new GMRT 320 MHz and JVLA 1.5 GHz observations of these two clusters. We found that both Abell 1451 and Zwcl 0634.1+4750 host a radio halo with a typical spectrum ($\alpha\sim1-1.3$). Similarly to a few other cases reported in the recent literature, these radio halos are significantly fainter in radio luminosity with respect to the current radio power-mass correlations and they are smaller than classical giant radio halos. These underluminous sources might contribute to shed light on the complex mechanisms of formation and evolution of radio halos. Furthermore, we detected a candidate radio relic at large distance from the cluster center in Abell 1451 and a peculiar head tail radio galaxy in Zwcl 0634.1+4750, which might be interacting with a shock front.
• We present near-IR spectroscopy of 22 luminous low-ionization broad absorption line quasars (LoBAL QSOs) at redshift 1.3<z<2.5, with 12 objects at z~1.5 and 10 at z~2.3. The spectra cover the rest-frame H$\alpha$ and H$\beta$ line regions, allowing us to obtain robust black hole mass estimates based on the broad H$\alpha$ line. We use these data, augmented by a lower redshift sample from the SDSS, to test the proposed youth scenario for LoBALs, which suggests LoBALs to constitute an early short lived evolutionary stage of quasar activity, by probing for any difference in their masses, Eddington ratios or rest-frame optical spectroscopic properties compared to normal quasars. In addition we construct the UV to mid-IR spectral energy distributions (SED) for the LoBAL sample and a matched non-BAL quasar sample. We do not find any statistically significant difference between LoBAL QSOs and non-BAL QSOs in their black hole mass or Eddington ratio distributions. The mean UV to mid-IR SED of the LoBAL QSOs is consistent with non-BAL QSOs, apart from their stronger reddening. At z>1 there is no clear difference in their optical emission line properties. We do not see particularly weak [OIII] nor strong FeII emission. The LoBAL QSOs do not show a stronger prevalence of ionized gas outflows as traced by the [OIII] line, compared to normal QSOs of similar luminosity. We conclude that the optical-MIR properties of LoBAL QSOs are consistent with the general quasar population and do not support them to constitute a special phase of AGN evolution.
• The existence of binary supermassive black holes (SBHs) is predicted by models of hierarchical galaxy formation. To date, only a single binary SBH has been imaged, at a projected separation of 7.3 parsecs. Here we report the detection of a candidate dual SBH with projected separation of 0.35 pc in the gas-rich interacting spiral galaxy NGC 7674 (Mrk 533). This peculiar Seyfert galaxy possesses a $\sim$0.7 kpc Z-shaped radio jet; the leading model for the formation of such sources postulates the presence of an uncoalesced binary SBH created during the infall of a satellite galaxy. Using very long baseline interferometry (VLBI), we imaged the central region of Mrk 533 at radio frequencies of 2, 5, 8 and 15 GHz. Two, possibly inverted-spectrum radio cores were detected at 15 GHz only; the 8-15 GHz spectral indices of the two cores are $\ge-0.33$ and $\ge-0.38$ ($\pm 30\%$), consistent with accreting SBHs. We derive a jet speed $\sim0.28c$ from multi-epoch parsec-scale data of the hotspot region, and a source age $\ge8.2\times10^3$ yrs.
• We present a multi-wavelength data analysis of the extended star-forming region S235 (hereafter E-S235), where two molecular clouds are present. In E-S235, using the 12CO (1-0) and 13CO (1-0) line data, a molecular cloud linked with the site "S235main" is traced in a velocity range [-24, -18] km s^-1, while the other one containing the sites S235A, S235B, and S235C (hereafter "S235ABC") is depicted in a velocity range [-18, -13] km s^-1. In the velocity space, these two clouds are separated by ~4 km s^-1, and are interconnected by a lower intensity intermediate velocity emission, tracing a broad bridge feature. In the velocity channel maps, a possible complementary molecular pair at [-21, -20] km s^-1 and [-16, -15] km s^-1 is also evident. The sites, "S235ABC", East 1, and South-West are spatially seen in the interface of two clouds. Together, these observed features are consistent with the predictions of numerical models of the cloud-cloud collision (CCC) process, favoring the onset of the CCC in E-S235 about 0.5 Myr ago. Deep UKIDSS near-infrared photometric analysis of point-like sources reveals significant clustering of young stellar populations toward the sites located at the junction, and the "S235main". The sites, "S235ABC" harbor young compact H II regions having dynamical ages of ~0.06-0.22 Myr, and these sites (including South-West and East 1) also contain dust clumps (having M_clump ~40 to 635 M_sun). Our observational findings suggest that the star formation activities (including massive stars) appear to be influenced by the CCC mechanism at the junction.
• The Large Ultraviolet / Optical / Infrared Surveyor (LUVOIR) is one of four large mission concepts currently undergoing community study for consideration by the 2020 Astronomy and Astrophysics Decadal Survey. The LUVOIR Ultraviolet Multi-Object Spectrograph, LUMOS, is being designed to support all of the UV science requirements of LUVOIR, from exoplanet host star characterization to tomography of circumgalactic halos to water plumes on outer solar system satellites. LUMOS offers point source and multi-object spectroscopy across the UV bandpass, with multiple resolution modes to support different science goals. The instrument will provide low (R = 8,000-18,000) and medium (R = 30,000-65,000) resolution modes across the far-ultraviolet (FUV: 100-200 nm) and near-ultraviolet (NUV: 200-400 nm) windows, and a very low resolution mode (R = 500) for spectroscopic investigations of extremely faint objects in the FUV. Imaging spectroscopy will be accomplished over a 3 x 1.6 arcminute field-of-view by employing holographically-ruled diffraction gratings to control optical aberrations, microshutter arrays (MSA), advanced optical coatings for high-throughput in the FUV, and next generation large-format photon-counting detectors. The spectroscopic capabilities of LUMOS are augmented by an FUV imaging channel (100-200nm, 13 milliarcsecond angular resolution, 2 x 2 arcminute field-of-view) that will employ a complement of narrow and medium-band filters. We present an overview of LUMOS' observing modes and estimated performance curves for effective area, spectral resolution, and imaging performance. Example "LUMOS 100-hour Highlights" observing programs are presented to demonstrate the potential power of LUVOIR's ultraviolet spectroscopic capabilities.
• Within the standard hierarchical structure formation scenario, Milky Way-mass dark matter haloes have hundreds of dark matter subhaloes with mass $\gtrsim 10^8 \, {\rm M_{\odot}}$. Over the lifetime of a galactic disc a fraction of these may pass close to the central region and interact with the disc. We extract the properties of subhaloes, such as their mass and trajectories, from a realistic cosmological simulation to study their potential effect on stellar discs. We find that massive subhalo impacts can generate disc heating, rings, bars, warps, lopsidedness as wells as spiral structures in the disc. Specifically, strong counter-rotating single-armed spiral structures form each time a massive subhalo passes through the disc. Such single-armed spirals wind up relatively quickly (over $1-2$ Gyrs) and are generally followed by co-rotating two-armed spiral structures that both develop and wind up more slowly. In our simulations self-gravity in the disc is not very strong and these spiral structures are found to be kinematic density waves. We demonstrate that there is a clear link between each spiral mode in the disc and a given subhalo that caused it, and by changing the mass of the subhalo we can modulate the strength of the spirals. Furthermore, we find that the majority of subhaloes interact with the disc impulsively, such that the strength of spirals generated by subhaloes is proportional to the total torque they exert. We conclude that only a handful of encounters with massive subhaloes is sufficient for re-generating and sustaining spiral structures in discs over their entire lifetime.
• We present the results of a new study of the relationship between infrared excess (IRX), UV spectral slope (beta) and stellar mass at redshifts 2<z<3, based on a deep Atacama Large Millimeter Array (ALMA) 1.3-mm continuum mosaic of the Hubble Ultra Deep Field (HUDF). Excluding the most heavily-obscured sources, we use a stacking analysis to show that z~2.5 star-forming galaxies in the mass range 9.25 <= log(M/Msun) <= 10.75 are fully consistent with the IRX-beta relation expected for a relatively grey attenuation curve, similar to the commonly adopted Calzetti law. Based on a large, mass complete, sample of 2 <= z <= 3 star-forming galaxies drawn from multiple surveys, we proceed to derive a new empirical relationship between beta and stellar mass, making it possible to predict UV attenuation (A_1600) and IRX as a function of stellar mass, for any assumed attenuation law. Once again, we find that z~2.5 star-forming galaxies follow A_1600-mass and IRX-mass relations consistent with a relatively grey attenuation law, and find no compelling evidence that star-forming galaxies at this epoch follow a reddening law as steep as the Small Magellanic Cloud (SMC) extinction curve. In fact, we use a simple simulation to demonstrate that previous determinations of the IRX-beta relation may have been biased toward low values of IRX at red values of beta, mimicking the signature expected for an SMC-like dust law. We show that this provides a plausible mechanism for reconciling apparently contradictory results in the literature and that, based on typical measurement uncertainties, stellar mass provides a cleaner prediction of UV attenuation than beta. Although the situation at lower stellar masses remains uncertain, we conclude that for 2<z<3 star-forming galaxies with log(M/Msun) >= 9.75, both the IRX-beta and IRX-mass relations are well described by a Calzetti-like attenuation law.
• Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used to study the young star cluster (YSC) population of this interacting system. Our newly extracted source catalogue contains 2834 cluster candidates, morphologically classified to be compact and uniform in colour, for which ages, masses and extinction are derived. In this first work we study the main properties of the YSC population of the whole galaxy, considering a mass-limited sample. Both luminosity and mass functions follow a power law shape with slope -2, but at high luminosities and masses a dearth of sources is observed. The analysis of the mass function suggests that it is best fitted by a Schechter function with slope -2 and a truncation mass at $1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm this result and link the shape of the luminosity function to the presence of a truncation in the mass function. A mass limited age function analysis, between 10 and 200 Myr, suggests that the cluster population is undergoing only moderate disruption. We observe little variation in the shape of the mass function at masses above $1\times10^4\ M_\odot$, over this age range. The fraction of star formation happening in the form of bound clusters in M51 is $\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over the whole range from 1 to 200 Myr.
• We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).
• We utilize the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to search for extended Lyman-Alpha emission around the z~6.6 QSO J0305-3150. After carefully subtracting the point-spread-function, we reach a nominal 5-sigma surface brightness limit of SB = 1.9x10$^{-18}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, collapsing 5 wavelength slices centered at the expected location of the redshifted Lyman-Alpha emission (i.e. at 9256 Ang.). Current data suggest the presence (5-sigma, accounting for systematics) of a Lyman-Alpha nebula that extends for 9 kpc around the QSO. This emission is displaced and redshifted by 155 km/s with respect to the location of the QSO host galaxy traced by the [CII] emission line. The total luminosity is L = 3.0x10$^{42}$ erg/s. Our analysis suggests that this emission is unlikely to rise from optically thick clouds illuminated by the ionizing radiation of the QSO. It is more plausible that the Lyman-Alpha emission is due to fluorescence of the highly ionized optically thin gas. This scenario implies a high hydrogen volume density of n$_H$ ~ 6 cm$^{-3}$. In addition, we detect a Lyman-Alpha emitter (LAE) in the immediate vicinity of the QSO: i.e., with a projected separation of 12.5 kpc and a line-of-sight velocity difference of 560 km/s. The luminosity of the LAE is L = 2.1x10$^{42}$ erg/s and its inferred star-formation-rate is SFR ~ 1.3 M$_\odot$/yr. The probability of finding such a close LAE is one order of magnitude above the expectations based on the QSO-galaxy cross-correlation function. This discovery is in agreement with a scenario where dissipative interactions favour the rapid build-up of super-massive black holes at early Cosmic times.
• Inspiraling massive black-hole binaries (MBHBs) forming in the aftermath of galaxy mergers are expected to be the loudest gravitational-wave (GW) sources relevant for pulsar-timing arrays (PTAs) at nHz frequencies. The incoherent overlap of signals from a cosmic population of MBHBs gives rise to a stochastic GW background (GWB) with characteristic strain around $h_c\sim10^{-15}$ at a reference frequency of 1 yr$^{-1}$, although uncertainties around this value are large. Current PTAs are piercing into the GW amplitude range predicted by state-of-the-art MBHB-population models, but no detection has been reported so far. To assess the future success prospects of PTA experiments, it is therefore important to estimate the minimum GWB level consistent with our current understanding of the formation and evolution of galaxies and massive black holes (MBHs). To this purpose, we couple a state-of-the-art semianalytic model of galaxy evolution and an extensive study of the statistical outcome of triple MBH interactions. We show that even in the most pessimistic scenario where all MBHBs stall before entering the GW-dominated regime, triple interactions resulting from subsequent galaxy mergers inevitably drive a considerable fraction of the MBHB population to coalescence. In the nHz frequency range relevant for PTA, the resulting GWB is only a factor of 2-to-3 suppressed compared to a fiducial model where binaries are allowed to merge over Gyr timescales after their host galaxies merge. Coupled with current estimates of the expected GWB amplitude range, our findings suggest that the minimum GWB from cosmic MBHBs is unlikely to be lower than $h_c\sim10^{-16}$ (at $f = 1$ yr$^{-1}$), well within the expected sensitivity of projected PTAs based on future observations with FAST, MeerKAT and SKA.
• We present a flux-limited sample of $z\sim0.3$ Ly$\alpha$ emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism spectroscopic data. The published GALEX $z\sim0.3$ LAE sample is pre-selected from continuum-bright objects and thus is biased against high equivalent width (EW) LAEs. We remove this continuum pre-selection and compute the EW distribution and the luminosity function of the Ly$\alpha$ emission line directly from our sample. We examine the evolution of these quantities from $z\sim0.3$ to $2.2$ and find that the EW distribution shows little evidence for evolution over this redshift range. As shown by previous studies, the Ly$\alpha$ luminosity density from star-forming galaxies declines rapidly with declining redshift. However, we find that the decline in Ly$\alpha$ luminosity density from $z=2.2$ to $z=0.3$ may simply mirror the decline seen in the H$\alpha$ luminosity density from $z=2.2$ to $z=0.4$, implying little change in the volumetric Ly$\alpha$ escape fraction. Finally, we show that the observed Ly$\alpha$ luminosity density from AGNs is comparable to the observed Ly$\alpha$ luminosity density from star-forming galaxies at $z=0.3$. We suggest that this significant contribution from AGNs to the total observed Ly$\alpha$ luminosity density persists out to $z\sim2.2$.
• Massive black hole binaries (MBHBs) are expected to form at the centre of merging galaxies during the hierarchical assembly of the cosmic large scale structure, and are therefore expected to be the loudest sources of gravitational waves (GWs) in the frequency window from nHz to tens of mHz. However, because of the dearth of relevant energy exchanges with background stars and gas, many of these MBHBs may stall at separations too large for GW emission to drive them to coalescence in less than a Hubble time. Triple MBH systems are then bound to form after a further galaxy merger, triggering a complex and rich dynamics that can eventually lead to MBH coalescence. Here we report on the results of a large set of numerical simulations performed with the code presented in Bonetti et al. (2016) where MBH triplets are set in spherical stellar potentials and MBH dynamics is followed through 2.5 post-Newtonian order in the equations of motion. We characterise each simulated system by the mass of the heavier MBH, the inner and outer mass ratios, the initial eccentricities of the inner and outer binaries, and the relative inclination, running a total of about 15k simulations. From our full suite of simulated systems we find that a fraction 20-30% of the MBH binaries that would otherwise stall are led to coalesce within a Hubble time. The corresponding coalescence timescale has a log-normal distribution, with a mean value around 250 Myr, while the eccentricity close to the plunge, albeit small, is non-negligible (~0.1). We construct and discuss marginalised probability distributions of the main parameters involved and, in a companion paper, we will use the results presented here to forecast the contribution of MBH triplets to the GW signal in the nHz regime probed by PTA experiments. In a follow-up paper, we will perform a similar exercise for MBHBs in the mHz regime targeted by LISA. [Abridged]