No Arabic abstract
In this first research note of a series of two, we conduct optical/UV investigations of the spectropolarimetric signatures emerging from the structure of quasars (Elvis 2000) applied to a purely theoretical, dusty model. We aim to explore the similarities/differences between an absorbing, disk-born outflow and the usual dusty torus that is supposed to hide the internal regions of active galactic nuclei (AGN). Using radiative transfer Monte Carlo simulations, we compute the continuum polarization signatures emerging from the model setup of Elvis (2000). We find that a dust-filled outflow produces very low amount of wavelength-depend polarization degrees, associated with a photon polarization angle perpendicular to the projected symmetry axis of the model. The polarization percentages are ten times lower than what can be produced by a toroidal model, with a maximal polarization degree found for intermediate viewing angle (i.e. when the observers line-of-sight crosses the outflowing material). The structure for quasars unsuccessfully blocks the radiation from the central irradiating source and shows a spectropolarimetric behavior that cannot be conciliated with observations. Either a new set of morphological parameters or different optical thickness must be considered.
Polarization has been detected at early times for all types of supernova, indicating that such systems result from or quickly develop some form of asymmetry. In addition, the detection of strong line polarization in supernovae is suggestive of chemical inhomogeneities (clumps) in the layers above the photosphere, which may reflect hydrodynamical instabilities during the explosion. We have developed a fast, flexible, approximate semi-analytic code for modeling polarized line radiative transfer within 3-D inhomogeneous rapidly-expanding atmospheres. Given a range of model parameters, the code generates random sets of clumps in the expanding ejecta and calculates the emergent line profile and Stokes parameters for each configuration. The ensemble of these configurations represents both the effects of various host geometries and of different viewing angles. We present results for the first part of our survey of model geometries, specifically the effects of the number and size of clumps (and the related effect of filling factor) on the emergent spectrum and Stokes parameters. Our simulations show that random clumpiness can produce line polarization in the range observed in SNe Ia (~1-2%), as well as the Q-U loops that are frequently seen in all SNe. We have also developed a method to connect the results of our simulations to robust observational parameters such as maximum polarization and polarized equivalent width in the line. Our models, in connection with spectropolarimetric observations, can constrain the 3-D structure of supernova ejecta and offer important insight into the SN explosion physics and the nature of their progenitor systems.
Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly-coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of resonant drag instabilities (RDIs) which de-couple the dust and gas dynamics and could qualitatively change the nonlinear outcome of these outflows. We present the first simulations of radiation-dust driven outflows in stratified, inhomogeneous media, including explicit grain dynamics and a realistic spectrum of grain sizes and charge, magnetic fields and Lorentz forces on grains (which dramatically enhance the RDIs), Coulomb and Epstein drag forces, and explicit radiation transport allowing for different grain absorption and scattering properties. In this paper we consider conditions resembling giant molecular clouds (GMCs), HII regions, and distributed starbursts, where optical depths are modest ($lesssim 1$), single-scattering effects dominate radiation-dust coupling, Lorentz forces dominate over drag on grains, and the fastest-growing RDIs are similar, such as magnetosonic and fast-gyro RDIs. These RDIs generically produce strong size-dependent dust clustering, growing nonlinear on timescales that are much shorter than the characteristic times of the outflow. The instabilities produce filamentary and plume-like or horsehead nebular morphologies that are remarkably similar to observed dust structures in GMCs and HII regions. Additionally, in some cases they strongly alter the magnetic field structure and topology relative to filaments. Despite driving strong micro-scale dust clumping which leaves some gas behind, an order-unity fraction of the gas is always efficiently entrained by dust.
We present evidence for the rapid variability of the high velocity iron K-shell absorption in the nearby ($z=0.184$) quasar PDS456. From a recent long Suzaku observation in 2013 ($sim1$Ms effective duration) we find that the the equivalent width of iron K absorption increases by a factor of $sim5$ during the observation, increasing from $<105$eV within the first 100ks of the observation, towards a maximum depth of $sim500$eV near the end. The implied outflow velocity of $sim0.25$c is consistent with that claimed from earlier (2007, 2011) Suzaku observations. The absorption varies on time-scales as short as $sim1$ week. We show that this variability can be equally well attributed to either (i) an increase in column density, plausibly associated with a clumpy time-variable outflow, or (ii) the decreasing ionization of a smooth homogeneous outflow which is in photo-ionization equilibrium with the local photon field. The variability allows a direct measure of absorber location, which is constrained to within $r=200-3500$$rm{r_{g}}$ of the black hole. Even in the most conservative case the kinetic power of the outflow is $gtrsim6%$ of the Eddington luminosity, with a mass outflow rate in excess of $sim40%$ of the Eddington accretion rate. The wind momentum rate is directly equivalent to the Eddington momentum rate which suggests that the flow may have been accelerated by continuum-scattering during an episode of Eddington-limited accretion.
Ultra-fast outflows (UFOs) are seen in many AGN, giving a possible mode for AGN feedback onto the host galaxy. However, the mechanism(s) for the launch and acceleration of these outflows are currently unknown, with UV line driving apparently strongly disfavoured as the material along the line of sight is so highly ionised that it has no UV transitions. We revisit this issue using the Suzaku X-ray data from PDS 456, an AGN with the most powerful UFO seen in the local Universe. We explore conditions in the wind by developing a new 3-D Monte-Carlo code for radiation transport. The code only handles highly ionised ions, but the data show the ionisation state of the wind is high enough that this is appropriate, and this restriction makes it fast enough to explore parameter space. We reproduce the results of earlier work, confirming that the mass loss rate in the wind is around 30% of the inferred inflow rate through the outer disc. We show for the first time that UV line driving is likely to be a major contribution to the wind acceleration. The mass loss rate in the wind matches that predicted from a purely line driven system, and this UV absorption can take place out of the line of sight. Continuum driving should also play a role as the source is close to Eddington. This predicts that the most extreme outflows will be produced from the highest mass accretion rate flows onto high mass black holes, as observed.
We investigate mass ejection from accretion disks formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH-NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disk configurations expected from BH-NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion disks that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower-limit to the fraction of the initial disk mass ejected. We find a nearly linear inverse dependence of this fraction on disk compactness (BH mass over initial disk radius). The dependence is related to the fraction of the disk mass accreted before the outflow is launched, which depends on the disk position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disk mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disk masses. We estimate the radioactive luminosity from the disk outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH-NS parameter space, the disk outflow contribution is well below the kilonova mass upper limits for GW190814.