Search

Multifrequency polarimetry is emerging as a powerful probe of blazar jets, especially due to the advent of the Imaging X-ray Polarimetry Explorer (IXPE) space observatory. We study the polarization of High-Synchrotron Peaked (HSP) blazars, where both optical and X-ray emission can be attributed to synchrotron radiation from a population of non-thermal electrons. We adopt an axisymmetric stationary force-free jet model, where the electromagnetic fields are determined by the jet shape. When the jet is nearly parabolic, the X-ray polarization degree is $\Pi_{\rm X}\sim 15-50\%$, and the optical polarization degree is $\Pi_{\rm O}\sim 5-25\%$. The polarization degree is strongly chromatic, as $\Pi_{\rm X}/\Pi_{\rm O}\sim 2-9$. The chromaticity is due to the softening of the electron distribution at high energies, and is much stronger than for a uniform magnetic field. The Electric Vector Position Angle (EVPA) is aligned with the projection of the jet axis on the plane of the sky. These results compare very well with multifrequency polarimetric observations of HSP blazars. Instead, when the jet is nearly cylindrical, the polarization degree is large and weakly chromatic (we find $\Pi_{\rm X}\sim 70\%$ and $\Pi_{\rm O}\sim 60\%$, close to the expected values for a uniform magnetic field). The EVPA is perpendicular to the projection of the jet axis on the plane of the sky. Then, a cylindrical geometry is practically ruled out by current observations. The polarization degree and the EVPA may be less sensitive to the specific particle acceleration process (e.g.,~magnetic reconnection or shocks) than previously thought., Comment: 8 pages, 4 figures, accepted by A&A

The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $\sigma_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $\sigma_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $\mu$G. Finally, we discuss the implications of our results within the observational context in the CMZ.

We investigate the impact on the properties of high-energy emitting BL Lac objects of a population of electrons with an anisotropic momentum distribution. We adopt a simple phenomenological description of the anisotropy, in which the most energetic electrons have a small pitch angle and the least energetic electrons are isotropic, as proposed by Sobacchi \& Lyubarsky (2019). We explore (i) a simple model that assumes a phenomenological shape for the electron energy distribution, and (ii) a self-consistent scheme in which the electrons follow a distribution which is the result of the balance between injection and radiative losses (we include the effects of the anisotropy on the synchrotron cooling rate). Considering the BL Lac object Mkn 421 as representative of the entire class, we show that in both cases the emission can be satisfactorily reproduced under equipartition between the magnetic field and the relativistic electrons. This is in better agreement with the idea that jets are launched as Poynting dominated flows with respect to the standard isotropic scenario, which requires both a low magnetization and a low radiative efficiency to reproduce the observed emission. The hard spectrum predicted for the inverse Compton continuum at TeV energies could be used as a potential test of the anisotropic model., Comment: 7 pages, 5 figures, accepted for publication by MNRAS

Theoretical modelling and observations of AGN jets suggest that the non-thermal electrons emitting the observed radiation should (i) carry an amount of energy comparable to the magnetic fields ($U_e\sim U_B$), which is likely the case if the magnetic fields play a dynamically important role in the jet's acceleration process; (ii) cool efficiently in a dynamical time ($t_{\rm cool}\lesssim t_{\rm dyn}$), which is suggested by the fact that a large fraction of the jet's kinetic energy is promptly converted into radiation. These expectations are at odds with the results of the simplest one-zone Self-Synchro-Compton (SSC) model for the Spectral Energy Distribution (SED) of BL Lacs. Indeed, the model predicts $U_e\gg U_B$ and $t_{\rm cool}\gg t_{\rm dyn}$ for most of the objects. Here we closely investigate one of the key assumptions of this model, namely that the momentum distribution of the non-thermal electrons is isotropic. We find that this assumption may be an oversimplification. If the magnetic energy is dissipated via a turbulent MHD cascade, the highest energy electrons may instead retain a small pitch angle. Since the synchrotron emissivity is suppressed when the pitch angle is small, this effect may importantly affect the modelling of the SED. As an illustrative example, we present an anisotropic model for the electron momentum distribution such that $U_e\sim U_B$ and $t_{\rm cool}\lesssim t_{\rm dyn}$ at the same time. Our model manages to simultaneously solve the two problems with one only more free parameter with respect to the usual isotropic one-zone SSC model., Comment: Submitted to MNRAS. Revised after the first referee report

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI., Comment: 24 pages, 8 figures, accepted for publication in PASA (SPICA special issue)

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features., Comment: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue

We study the CO line luminosity ($L_{\rm CO}$), the shape of the CO Spectral Line Energy Distribution (SLED), and the value of the CO-to-$\rm H_2$ conversion factor in galaxies in the Epoch of Reionization (EoR). To this aim, we construct a model that simultaneously takes into account the radiative transfer and the clumpy structure of giant molecular clouds (GMCs) where the CO lines are excited. We then use it to post-process state-of-the-art zoomed, high resolution ($30\, \rm{pc}$), cosmological simulation of a main-sequence ($M_{*}\approx10^{10}\, \rm{M_{\odot}}$, $SFR\approx 100\,\rm{M_{\odot}\, yr^{-1}}$) galaxy, "Alth{\ae}a", at $z\approx6$. We find that the CO emission traces the inner molecular disk ($r\approx 0.5 \,\rm{kpc}$) of Alth{\ae}a with the peak of the CO surface brightness co-located with that of the [CII] 158$\rm \mu m$ emission. Its $L_{\rm CO(1-0)}=10^{4.85}\, \rm{L_{\odot}}$ is comparable to that observed in local galaxies with similar stellar mass. The high ($\Sigma_{gas} \approx 220\, \rm M_{\odot}\, pc^{-2}$) gas surface density in Alth{\ae}a, its large Mach number (\mach$\approx 30$), and the warm kinetic temperature ($T_{k}\approx 45 \, \rm K$) of GMCs yield a CO SLED peaked at the CO(7-6) transition, i.e. at relatively high-$J$, and a CO-to-$\rm H_2$ conversion factor $\alpha_{\rm CO}\approx 1.5 \, \rm M_{\odot} \rm (K\, km\, s^{-1}\, pc^2)^{-1} $ lower than that of the Milky Way. The ALMA observing time required to detect (resolve) at 5$\sigma$ the CO(7-6) line from galaxies similar to Alth{\ae}a is $\approx13$ h ($\approx 38$ h)., Comment: 16 pages, 14 figures, accepted for publication in MNRAS

In the paradigm of magnetic launching of astrophysical jets, instabilities in the MHD flow are a good candidate to convert the Poynting flux into the kinetic energy of the plasma. If the magnetised plasma fills the almost entire space, the jet is unstable to helical perturbations of its body. However, the growth rate of these modes is suppressed when the poloidal component of the magnetic field has a vanishing gradient, which may be the actual case for a realistic configuration. Here we show that, if the magnetised plasma is confined into a limited region by the pressure of some external medium, the velocity shear at the contact surface excites unstable modes which can affect a significant fraction of the jet's body. We find that when the Lorentz factor of the jet is $\Gamma\sim10$ ($\Gamma\sim 100$), these perturbations typically develop after propagating along the jet for tens (hundreds) of jet's radii. Surface modes may therefore play an important role in converting the energy of the jet from the Poynting flux to the kinetic energy of the plasma, particularly in AGN. The scaling of the dispersion relation with (i) the angular velocity of the field lines and (ii) the sound speed in the confining gas is discussed., Comment: 9 pages, accepted by MNRAS

Core-collapse supernovae (SNe) of Types Ib and Ic arise from hydrogen-stripped stars, while the latter are also stripped of their helium. Both SN types have a similar temporal evolution, suggesting broadly similar progenitors. However, while some Ic SNe are associated with relativistic jets of gamma-ray bursts (GRBs), no GRB has yet been found in association with Type Ib SNe. Here we find that, even if GRB-like central engines operate in both SNe Ib and Ic, different properties of their envelopes may accommodate this potential tension. In particular, we focus on the case of a low-mass, extended envelope surrounding the progenitor star (as produced, for example, by strong mass losses prior to explosion). If the envelopes of Type Ib SNe are sufficiently massive ($M_{\rm ext}\sim(0.3-1)M_\odot$) and extended ($R_{\rm ext}\sim 10^{13}\text{ cm}$), we show that (i) GRB-like jets cannot break out of the star; (ii) the SN light curve is compatible with current observations. Different envelope properties of Type Ib SNe with respect to (at least some) Type Ic SNe may be connected to the presence of a helium layer surrounding the star., Comment: We realised that eq. 4 and 5 are not appropriate if a large amount of energy is deposited directly into the envelope. We are trying to fix this issue, but this is not immediate and can affect the conclusions. Hence, we prefer to withdraw the paper

Long-duration, spectrally-soft Gamma-Ray Bursts (GRBs) are associated with Type Ic Core Collapse (CC) Supernovae (SNe), and thus arise from the death of massive stars. In the collapsar model, the jet launched by the central engine must bore its way out of the progenitor star before it can produce a GRB. Most of these jets do not break out, and are instead "choked" inside the star, as the central-engine activity time, $t_{\rm e}$, is not long enough. Modelling the long-soft GRB duration distribution assuming a power-law distribution for their central-engine activity times, $\propto t_{\rm e}^{-\alpha}$ for $t_{\rm e}>t_{\rm b}$, we find a steep distribution ($\alpha\sim4$) and a typical GRB jet breakout time of $t_{\rm b}\sim 60\text{ s}$ in the star's frame. The latter suggests the presence of a low-density, extended envelope surrounding the progenitor star, similar to that previously inferred for low-luminosity GRBs. Extrapolating the range of validity of this power law below what is directly observable, to $t_{\rm e}

In the paradigm of magnetic acceleration of relativistic jets, one of the key points is identifying a viable mechanism to convert the Poynting flux into the kinetic energy of the plasma beyond equipartition. A promising candidate is the kink instability, which deforms the body of the jet through helical perturbations. Since the detailed structure of real jets is unknown, we explore a large family of cylindrical, force-free equilibria to get robust conclusions. We find that the growth rate of the instability depends primarily on two parameters: (i) the gradient of the poloidal magnetic field; (ii) the Lorentz factor of the perturbation, which is closely related to the velocity of the plasma. We provide a simple fitting formula for the growth rate of the instability. As a tentative application, we use our results to interpret the dynamics of the jet in the nearby active galaxy M87. We show that the kink instability becomes non-linear at a distance from the central black hole comparable to where the jet stops accelerating. Hence (at least for this object), the kink instability of the jet is a good candidate to drive the transition from a Poynting-dominated to a kinetic-energy-dominated flow., Comment: 7 pages, 5 figures, accepted by MNRAS. Minor revision after referee comments

We investigate the impact of sinks of ionizing radiation on the reionization-era 21-cm signal, focusing on 1-point statistics. We consider sinks in both the intergalactic medium and inside galaxies. At a fixed filling factor of HII regions, sinks will have two main effects on the 21-cm morphology: (i) as inhomogeneous absorbers of ionizing photons they result in smaller and more widespread cosmic HII patches; and (ii) as reservoirs of neutral gas they contribute a non-zero 21-cm signal in otherwise ionized regions. Both effects damp the contrast between neutral and ionized patches during reionization, making detection of the epoch of reionization with 21-cm interferometry more challenging. Here we systematically investigate these effects using the latest semi-numerical simulations. We find that sinks dramatically suppress the peak in the redshift evolution of the variance, corresponding to the midpoint of reionization. As previously predicted, skewness changes sign at midpoint, but the fluctuations in the residual HI suppress a late-time rise. Furthermore, large levels of residual HI dramatically alter the evolution of the variance, skewness and power spectrum from that seen at lower levels. In general, the evolution of the large-scale modes provides a better, cleaner, higher signal-to-noise probe of reionization., Comment: Minor edits to agree with MNRAS published version