Massive Torus Research Articles

Supernovalike explosions of massive rotating stars from disks surrounding a black hole

We perform a new general-relativistic viscous-radiation hydrodynamics simulation for supernovalike explosions associated with stellar core collapse of rotating massive stars to a system of a black hole and a massive torus, paying particular attention to large-mass progenitor stars with the zero-age main-sequence mass of MZAMS=20, 35, and 45M⊙ of Aguilera-Dena [ ]. Assuming that a black hole is formed in a short timescale after the onset of the stellar collapse, the new simulations are started from initial data of a spinning black hole and infalling matter that self-consistently satisfy the constraint equations of general relativity. It is found that, with a reasonable size of the viscous parameter, the supernovalike explosion is driven by the viscous heating effect in the torus around the black hole, irrespective of the progenitor mass. The typical explosion energy and ejecta mass for the large-mass cases (MZAMS=35 and 45M⊙) are ∼1052 erg and ∼5M⊙, respectively, with Ni56 mass larger than 0.15M⊙. These are consistent with the observational data of stripped-envelope and high-energy supernovae such as broad-lined type Ic supernovae. This indicates that rotating stellar collapses of massive stars to a black hole surrounded by a massive torus can be a central engine for high-energy supernovae. By artificially varying the angular velocity of the initial data, we explore the dependence of the explosion energy and ejecta mass on the initial angular momentum and find that the large explosion energy ∼1052 erg and large Ni56 mass ≥0.15M⊙ are possible only when a large-mass compact torus with mass ≳1M⊙ is formed. Published by the American Physical Society 2024

Contact Mechanics between Torus and Pit Border for Developing Air-Seeding Seal in Aspirated Bordered Pits

Bordered pits in conifers have been recognized as a significant evolutionary characteristic that served to impede the spread of embolisms between tracheids. Nevertheless, there was a lack of comprehensive understanding regarding the mechanical properties of the torus and the pit border in relation to the formation of contact seals in aspirated pits. A solid mechanics model was developed to study aspirated bordered pits, incorporating the elastic deformation of the torus–margo structure and the contact behavior between the torus and the pit border. Ten pit samples were reconstructed using transmission electron microscopy (TEM) and scanning electron microscopy (SEM) pictures in order to develop mechanical models for pits in the stems of Platycladus orientalis. Due to the limited contact area on the torus surface, the contact pressure between the torus and the pit border exceeded the air pressure in aspirated bordered pits. The external force and the duration required to seal pits decreased when the pit diameter increased and the pit depth decreased. The augmentation of the torus and margo mass necessitated a lengthier duration for the closure of cavities. The findings presented in this study offered theoretical support for the capillary-seeding hypothesis. The aspiration in bordered pits mechanically depended on the structural and material characteristics of the torus and margo.

Self-Consistent Picture of the Mass Ejection from a One Second Long Binary Neutron Star Merger Leaving a Short-Lived Remnant in a General-Relativistic Neutrino-Radiation Magnetohydrodynamic Simulation.

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on the Japanese supercomputer Fugaku using about 85million CPU hours with 20 736 CPUs. We consider an asymmetric binary neutron star merger with masses of 1.2M_{⊙} and 1.5M_{⊙} and a "soft" equation of state SFHo. It results in a short-lived remnant with the lifetime of ≈0.017 s, and subsequent massive torus formation with the mass of ≈0.05M_{⊙} after the remnant collapses to a black hole. For the first time, we find that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the postmerger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity in a single simulation and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

Neutrino absorption and other physics dependencies in neutrino-cooled black hole accretion discs

ABSTRACT Black hole (BH) accretion discs formed in compact-object mergers or collapsars may be major sites of the rapid-neutron-capture (r-)process, but the conditions determining the electron fraction (Ye) remain uncertain given the complexity of neutrino transfer and angular-momentum transport. After discussing relevant weak-interaction regimes, we study the role of neutrino absorption for shaping Ye using an extensive set of simulations performed with two-moment neutrino transport and again without neutrino absorption. We vary the torus mass, BH mass and spin, and examine the impact of rest-mass and weak-magnetism corrections in the neutrino rates. We also test the dependence on the angular-momentum transport treatment by comparing axisymmetric models using the standard α-viscosity with viscous models assuming constant viscous length-scales (lt) and 3D magnetohydrodynamic (MHD) simulations. Finally, we discuss the nucleosynthesis yields and basic kilonova properties. We find that absorption pushes Ye towards ∼0.5 outside the torus, while inside increasing the equilibrium value $Y_\mathrm{ e}^{\mathrm{eq}}$ by ∼0.05–0.2. Correspondingly, a substantial ejecta fraction is pushed above Ye = 0.25, leading to a reduced lanthanide fraction and a brighter, earlier, and bluer kilonova than without absorption. More compact tori with higher neutrino optical depth, τ, tend to have lower $Y_\mathrm{ e}^{\mathrm{eq}}$ up to τ ∼ 1–10, above which absorption becomes strong enough to reverse this trend. Disc ejecta are less (more) neutron rich when employing an lt = const. viscosity (MHD treatment). The solar-like abundance pattern found for our MHD model marginally supports collapsar discs as major r-process sites, although a strong r-process may be limited to phases of high mass-infall rates, $\dot{M}\, \, \raise0.14em\rm{\gt }\lower0.28em\rm{\sim }\, \, 2\times 10^{-2}$ M⊙ s−1.

Ultra-delayed Neutrino-driven Explosion of Rotating Massive-star Collapse

Abstract Long-term neutrino-radiation hydrodynamics simulations in full general relativity are performed for the collapse of rotating massive stars that are evolved from He-stars with initial masses of 20 and 32 M ⊙. It is shown that if the collapsing stellar core has sufficient angular momentum, the rotationally supported proto-neutron star (PNS) survives for seconds accompanying the formation of a massive torus of mass larger than 1 M ⊙. Subsequent mass accretion onto the central region produces a massive and compact central object, and eventually enhances the neutrino luminosity beyond 1053 erg s−1, resulting in a very delayed neutrino-driven explosion, in particular toward the polar direction. The kinetic energy of the explosion can be appreciably higher than 1052 erg for a massive progenitor star and compatible with that of energetic supernovae like broad-line type-Ic supernovae. By the subsequent accretion, the massive PNS collapses eventually into a rapidly spinning black hole, which could be a central engine for gamma-ray bursts if a massive torus surrounds it.

Viscous evolution of a massive disk surrounding stellar-mass black holes in full general relativity

Long-term viscous neutrino-radiation hydrodynamics simulations in full general relativity are performed for a massive disk surrounding spinning stellar-mass black holes with mass $M_{\rm BH}=4$, $6$, and $10M_\odot$ and initial dimensionless spin $\chi \approx 0.8$. The initial disk is chosen to have mass $M_{\rm disk}\approx 0.1$ or $3M_\odot$ as plausible models of the remnants for the merger of black hole-neutron star binaries or the stellar core collapse from a rapidly rotating progenitor, respectively. For $M_{\rm disk} \approx 0.1M_\odot$ with the outer disk edge initially located at $r_{\rm out} \sim 200$ km, we find that $15$%-$20$% of $M_{\rm disk}$ is ejected and the average electron fraction of the ejecta is $\langle Y_e \rangle = 0.30$-$0.35$ as found in the previous study. For $M_{\rm disk} \approx 3M_\odot$, we find that $\approx 10$%-$20$% of $M_{\rm disk}$ is ejected for $r_{\rm out}\approx 200$-$1000$ km. In addition, $\langle Y_e \rangle$ of the ejecta can be enhanced to be $\gtrsim 0.4$ because the electron fraction is increased significantly during the long-term viscous expansion of the disk with high neutrino luminosity until the mass ejection sets in. Our results suggest that not heavy $r$-process elements but light trans-iron elements would be synthesized in the matter ejected from a massive torus surrounding stellar-mass black holes. We also find that the outcomes of the viscous evolution for the high-mass disk case is composed of a rapidly spinning black hole surrounded by a torus with a narrow funnel, which appears to be suitable for generating gamma-ray bursts.

Axisymmetric hydrodynamics in numerical relativity using a multipatch method

We describe a method of implementing the axisymmetric evolution of general-relativistic hydrodynamics and magnetohydrodynamics through modification of a multipatch grid scheme. In order to ease the computational requirements required to evolve the post-merger phase of systems involving binary compact massive objects in numerical relativity, it is often beneficial to take advantage of these system’s tendency to rapidly settle into states that are nearly axisymmetric, allowing for 2D evolution of secular timescales. We implement this scheme in the spectral Einstein code and show the results of application of this method to four test systems including viscosity, magnetic fields, and neutrino radiation transport. Our results show that this method can be used to quickly allow already existing 3D infrastructure that makes use of local coordinate system transformations to be made to run in axisymmetric 2D with the flexible grid creation capabilities of multipatch methods. Our code tests include a simple model of a binary neutron star postmerger remnant, for which we confirm the formation of a massive torus which is a promising source of post-merger ejecta.

Planetary nebulae with UVIT

The high excitation planetary nebula, NGC 6302, has been imaged in two far-ultraviolet (FUV) filters, F169M (Sapphire; λeff: 1608 Å) and F172M (Silica; λeff: 1717 Å) and two near-UV (NUV) filters, N219M (B15; λeff: 2196 Å) and N279N (N2; λeff: 2792 Å) with the Ultra Violet Imaging Telescope (UVIT). The FUV F169M image shows faint emission lobes that extend to about 5 arcmin on either side of the central source. Faint orthogonal collimated jet-like structures are present on either side of the FUV lobes through the central source. These structures are not present in the two NUV filters or in the FUV F172M filter. Optical and infrared (IR) images of NGC 6302 show bright emission bipolar lobes in the east-west direction with a massive torus of molecular gas and dust seen as a dark lane in the north-south direction. The FUV lobes are much more extended and oriented at a position angle of 113°. They and the jet-like structures might be remnants of an earlier evolutionary phase, prior to the dramatic explosive event that triggered the Hubble type bipolar flows approximately 2200 years ago. The source of the FUV lobe and jet emission is not known, but is likely due to fluorescent emission from H2 molecules. The cause of the difference in orientation of optical and FUV lobes is not clear and, we speculate, could be related to two binary interactions.

Long-term GRMHD simulations of neutron star merger accretion discs: implications for electromagnetic counterparts

We investigate the long-term evolution of black hole accretion disks formed in neutron star mergers. These disks expel matter that contributes to an $r$-process kilonova, and can produce relativistic jets powering short gamma-ray bursts. Here we report the results of a three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulation of such a disk which is evolved for long enough ($\sim 9$s, or $\sim 6\times 10^5 r_{\rm g}/c$) to achieve completion of mass ejection far from the disk. Our model starts with a poloidal field, and fully resolves the most unstable mode of the magnetorotational instability. We parameterize the dominant microphysics and neutrino cooling effects, and compare with axisymmetric hydrodynamic models with shear viscosity. The GRMHD model ejects mass in two ways: a prompt MHD-mediated outflow and a late-time, thermally-driven wind once the disk becomes advective. The total amount of unbound mass ejected ($0.013M_\odot$, or $\simeq 40\%$ of the initial torus mass) is twice as much as in hydrodynamic models, with higher average velocity ($0.1c$) and a broad electron fraction distribution with a lower average value ($0.16$). Scaling the ejected fractions to a disk mass of $\sim 0.1M_\odot$ can account for the red kilonova from GW170817 but underpredicts the blue component. About $\sim 10^{-3}M_\odot$ of material should undergo neutron freezout and could produce a bright kilonova precursor in the first few hours after the merger. With our idealized initial magnetic field configuration, we obtain a robust jet and sufficient ejecta with Lorentz factor $\sim 1-10$ to (over)produce the non-thermal emission from GW1708107.

Expansion of the Potential of a Homogeneous Circular Torus in Terms of Geometrical Parameter

The spatial potential of a homogeneous gravitating (or statically charged) circular torus is expanded in terms of geometrical parameter q. The first term of the expansion (coefficient of zero-power parameter q) coincides with the potential of a thin ring the mass of which is equal to the mass of torus and the radius of which is equal to the radius of the axial line. It is shown that the coefficients of terms with odd powers are zeros. Even terms of the expansion of the potential are analytically represented in terms of standard complete elliptic integrals. The corresponding series describes the potential of torus in entire space including internal region. The proposed method for representation of the potential makes it possible to obtain gravitational energy of a homogeneous circular torus. The method is used to calculate masses of two rings of the Chariklo asteroid. The mass of the internal ring is Mr1 ≈ 9.8 × 1018 g and the ratio of such a mass to the mass of asteroid is $$\frac{{{M_{r1}}}}{{{M_0}}}$$ ≈ 0.001. The results for the external ring are Mr2 ≈ 1018 g and $$\frac{{{M_{r1}}}}{{{M_0}}}$$ ≈ 10–4.

Three-Dimensional General-Relativistic Magnetohydrodynamic Simulations of Remnant Accretion Disks from Neutron Star Mergers: Outflows and r-Process Nucleosynthesis.

The merger of binary neutron stars, or of a neutron star and a stellar-mass black hole, can result in the formation of a massive rotating torus around a spinning black hole. In addition to providing collimating media for γ-ray burst jets, unbound outflows from these disks are an important source of mass ejection and rapid neutron capture (r-process) nucleosynthesis. We present the first three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations of neutrino-cooled accretion disks in neutron star mergers, including a realistic equation of state valid at low densities and temperatures, self-consistent evolution of the electron fraction, and neutrino cooling through an approximate leakage scheme. After initial magnetic field amplification by magnetic winding, we witness the vigorous onset of turbulence driven by the magnetorotational instability (MRI). The disk quickly reaches a balance between heating from MRI-driven turbulence and neutrino cooling, which regulates the midplane electron fraction to a low equilibrium value Y_{e}≈0.1. Over the 380-ms duration of the simulation, we find that a fraction ≈20% of the initial torus mass is unbound in powerful outflows with asymptotic velocities v≈0.1c and electron fractions Y_{e}≈0.1-0.25. Postprocessing the outflows through a nuclear reaction network shows the production of a robust second- and third-peak r process. Though broadly consistent with the results of previous axisymmetric hydrodynamical simulations, extrapolation of our results to late times suggests that the total ejecta mass from GRMHD disks is significantly higher. Our results provide strong evidence that postmerger disk outflows are an important site for the r process.

General relativistic viscous hydrodynamics of differentially rotating neutron stars

Employing a simplified version of the Israel-Stewart formalism for general-relativistic shear-viscous hydrodynamics, we perform axisymmetric general-relativistic simulations for a rotating neutron star surrounded by a massive torus, which can be formed from differentially rotating stars. We show that with our choice of a shear-viscous hydrodynamics formalism, the simulations can be stably performed for a long time scale. We also demonstrate that with a possibly high shear-viscous coefficient, not only viscous angular momentum transport works but also an outflow could be driven from a hot envelope around the neutron star for a time scale $\ensuremath{\gtrsim}100\text{ }\text{ }\mathrm{ms}$ with the ejecta mass $\ensuremath{\gtrsim}{10}^{\ensuremath{-}2}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$, which is comparable to the typical mass for dynamical ejecta of binary neutron-star mergers. This suggests that massive neutron stars surrounded by a massive torus, which are typical outcomes formed after the merger of binary neutron stars, could be the dominant source for providing neutron-rich ejecta, if the effective shear viscosity is sufficiently high, i.e., if the viscous $\ensuremath{\alpha}$ parameter is $\ensuremath{\gtrsim}{10}^{\ensuremath{-}2}$. The present numerical result indicates the importance of a future high-resolution magnetohydrodynamics simulation that is the unique approach to clarify the viscous effect in the merger remnants of binary neutron stars by the first-principle manner.

Probing the active galactic nucleus unified model torus properties in Seyfert galaxies

We studied the physical parameters of a sample comprising of all Spitzer/IRS public spectra of Seyfert galaxies in the mid-infrared (5.2-38$\mu$m range) under the active galactic nuclei (AGN) unified model. We compare the observed spectra with $\sim10^6$ CLUMPY model spectral energy distributions, which consider a torus composed of dusty clouds. We find a slight difference in the distribution of line-of-sight inclination angle, $i$, requiring larger angles for Seyfert 2 (Sy2) and a broader distribution for Seyfert 1 (Sy1). We found small differences in the torus angular width, $\sigma$, indicating that Sy1 may host a slightly narrower torus than Sy2. The torus thickness, together with the bolometric luminosities derived, suggest a very compact torus up to $\sim$6 pc from the central AGN. The number of clouds along the equatorial plane, $N$, as well the index of the radial profile, $q$, are nearly the same for both types. These results imply that the torus cloud distribution is nearly the same for type 1 and type 2 objects. The torus mass is almost the same for both types of activity, with values in the range of $M_{tor}\sim$10$^{4}-$10$^{7}\rm M_{\odot}$. The main difference appears to be related to the clouds' intrinsic properties: type 2 sources present higher optical depths $\tau_V$. The results presented here reinforce the suggestion that the classification of a galaxy may depend also on the intrinsic properties of the torus clouds rather than simply on their inclination. This is in contradiction with the simple geometric idea of the unification model.

On extreme transient events from rotating black holes and their gravitational wave emission

The super-luminous object ASASSN-15lh (SN2015L) is an extreme event with a total energy $E_{rad}\simeq 1.1\times 10^{52}$ erg in black body radiation on par with its kinetic energy $E_k$ in ejecta and a late time plateau in the UV, that defies a nuclear origin. It likely presents a new explosion mechanism for hydrogen-deprived supernovae. With no radio emission and no H-rich environment we propose to identify $E_{rad}$ with dissipation of a baryon-poor outflow in the optically thick remnant stellar envelope produced by a central engine. By negligible time scales of light crossing and radiative cooling of the envelope, SN2015L's light curve closely tracks the evolution of this engine. We here model its light curve by the evolution of black hole spin, during angular momentum loss in Alv\'en waves to matter at the Inner Most Stable Circular Orbit (ISCO). The duration is determined by $\sigma=M_T/M$ of the torus mass $M_T$ around the black hole of mass $M$: $\sigma\sim 10^{-7}$ and $\sigma\sim 10^{-2}$ for SN2015L and, respectively, a long GRB. The observed electromagnetic radiation herein represents a minor output of the rotational energy $E_{rot}$ of the black hole, while most is radiated unseen in gravitational radiation. This model explains the high-mass slow-spin binary progenitor of GWB150914, as the remnant of two CC-SNe in an intra-day binary of two massive stars. This model rigorously predicts a change in magnitude $\Delta m\simeq 1.15$ in the light curve post-peak, in agreement with the light curve of SN2015L with no fine-tuning.

NEUTRON-STAR MERGER EJECTA AS OBSTACLES TO NEUTRINO-POWERED JETS OF GAMMA-RAY BURSTS

ABSTRACT We present the first special relativistic, axisymmetric hydrodynamic simulations of black hole-torus systems (approximating general relativistic gravity) as remnants of binary-neutron star (NS–NS) and neutron star–black hole (NS–BH) mergers, in which the viscously driven evolution of the accretion torus is followed with self-consistent energy-dependent neutrino transport and the interaction with the cloud of dynamical ejecta expelled during the NS–NS merging is taken into account. The modeled torus masses, BH masses and spins, and the ejecta masses, velocities, and spatial distributions are adopted from relativistic merger simulations. We find that energy deposition by neutrino annihilation can accelerate outflows with initially high Lorentz factors along polar low-density funnels, but only in mergers with extremely low baryon pollution in the polar regions. NS–BH mergers, where polar mass ejection during the merging phase is absent, provide sufficiently baryon-poor environments to enable neutrino-powered, ultrarelativistic jets with terminal Lorentz factors above 100 and considerable dynamical collimation, favoring short gamma-ray bursts (sGRBs), although their typical energies and durations might be too small to explain the majority of events. In the case of NS–NS mergers, however, neutrino emission of the accreting and viscously spreading torus is too short and too weak to yield enough energy for the outflows to break out from the surrounding ejecta shell as highly relativistic jets. We conclude that neutrino annihilation alone cannot power sGRBs from NS–NS mergers.

Going with the flow: using gas clouds to probe the accretion flow feeding Sgr A*

The massive black hole in our galactic center, Sgr A*, accretes only a small fraction of the gas available at its Bondi radius. The physical processes determining this accretion rate remain unknown, partly due to a lack of observational constraints on the gas at distances between ~10 and ~10$^5$ Schwarzschild radii (Rs) from the black hole. Recent infrared observations identify low-mass gas clouds, G1 and G2, moving on highly eccentric, nearly co-planar orbits through the accretion flow around Sgr A*. Although it is not yet clear whether these objects contain embedded stars, their extended gaseous envelopes evolve independently as gas clouds. In this paper we attempt to use these gas clouds to constrain the properties of the accretion flow at ~10$^3$ Rs. Assuming that G1 and G2 follow the same trajectory, we model the small differences in their orbital parameters as evolution resulting from interaction with the background flow. We find evolution consistent with the G-clouds originating in the clockwise disk. Our analysis enables the first unique determination of the rotation axis of the accretion flow: we localize the rotation axis to within 20 degrees, finding an orientation consistent with the parsec-scale jet identified in x-ray observations and with the circumnuclear disk, a massive torus of molecular gas ~1.5 pc from Sgr A*. This suggests that the gas in the accretion flow comes predominantly from the circumnuclear disk, rather than the winds of stars in the young clockwise disk. This result will be tested by the Event Horizon Telescope within the next year. Our model also makes testable predictions for the orbital evolution of G1 and G2, falsifiable on a 5-10 year timescale.

Comprehensive nucleosynthesis analysis for ejecta of compact binary mergers

We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

Investigation of dust properties of the proto-planetary nebula IRAS 18276−1431

We investigate the circumstellar dust properties of the oxygen-rich bipolar proto-planetary nebula (PPN) IRAS 18276−1431 by means of two-dimensional radiative transfer simulations of the circumstellar dust shell. The model geometry is assumed to have a torus and an envelope which consists of a pair of bipolar lobes and a spherical asymptotic giant branch shell. The parameters of the dust and the dust shell are constrained by comparing the spectral energy distribution (SED) and near-infrared intensity and polarization data with the models. The polarization in the envelope reaches 50–60 per cent and is nearly constant in the H and KS bands in the observations. This weak wavelength dependence of the polarization can be reproduced with a grain-size distribution function for the torus: 0.05 μm ≤ a with n(a)∝a− (p = 5.5)exp ( − a/ac = 0.3 μm). The power index p is significantly steeper than that for interstellar dust (p ∼ 3). Similar results have also been found in some other PPNs and suggest that mechanisms that grind down large particles, such as sputtering, may also have acted when the dust particles formed. The spectral opacity index β is found to be 0.6 ± 0.5 from the 760 μm to 2.6 mm fluxes, which is characterized by the dust in the torus. This low value (<2) indicates the presence of large dust grains in the torus. We discuss two possible dust models for the torus. One has a size distribution function of 1.0 ≤ a ≤ amax = 5 000.0 μm with n(a) ∝ a−(p = 2.5) and the other is 1.0 μm ≤ a ≤ amax = 10 000.0 μm with n(a) ∝ a−(p = 3.5). The former has β of 0.633, but we are not able to find reasonable geometry parameters to fit the SED in the infrared. The latter has β of 1.12, but reproduces the SED better over a wide wavelength range. With this dust model, the geometric parameters are estimated as follows: the inner and outer radii are 30 and 1000 au and the torus mass is 3.0 M⊙. Given that the torii are generally not found to be rotating, a large fraction of the torus material is likely to be expanding. Assuming an expansion velocity of 15 km s−1, the torus formation time and mass-loss rate are found to be ∼300 yr and ∼10−2 M⊙ yr−1, respectively.

MODELING THE NUCLEAR INFRARED SPECTRAL ENERGY DISTRIBUTION OF TYPE II ACTIVE GALACTIC NUCLEI

We present results from model fitting to the Spectral Energy Distribution (SED) of a homogeneous sample of Seyfert II galaxies drawn from the $12\mu$m Galaxy Sample. Imaging and nuclear flux measurements are presented in an accompanying paper (Videla et al., 2013). Here we add IRS Spitzer observations to further constrain the SEDs after careful subtraction of a starburst component. We use the library of CLUMPY torus models from Nenkova et al.~(2008ab) and also test the two-phase models recently produced by Stalevski et al.~(2012). We find that photometric and spectroscopic observations in the mid-IR (>5mu) are crucial to properly constrain the best-fit torus models. About half of our sources show clear near-IR excess of their SEDs above the best fit models. This problem can be less severe when using the Stalevski et al.~(2012) models. It is not clear what is the nature of this emission since best fitted black body temperatures are very high (~1700-2500 K) and the Type II classification of our sources would correspond to a small probability to peer directly into the hottest regions of the torus. Crucially, the derived torus parameters when using CLUMPY models are quite robust,, independently of whether the sources require an additional black body component or not. Our findings suggest that tori are characterized by N_0 > 5, sigma > 40, tau < 25, i > 40 degrees, Y < 50 and A_v^los ~ 100-300. From these we can determine that typical torus sizes and masses of 0.1-5.0 pc and 10^{4-6} M_\odot. We find tentative evidence that those nuclei with a detected Hidden Broad Line Regions are characterized by lower levels of extinction than those without one. Finally, we find no correlation between the torus properties and the presence of circumnuclear or more global star-formation.

COMPACT BINARY PROGENITORS OF SHORT GAMMA-RAY BURSTS

In recent years, detailed observations and accurate numerical simulations have provided support to the idea that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms can extract energy from this system and power the SGRBs. Here we connect observations and numerical simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general-relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By assuming a constant efficiency in converting torus mass into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses >~1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some of the SGRBs, the BH would need to have an initial spin of ~0.9, or higher.