Source author record

Emanuele Sobacchi

Emanuele Sobacchi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2026arXiv

Temporal variability of polarization in blazars

We investigate the temporal variability of polarization of synchrotron radiation from blazar jets. Multiwavelength observations revealed high-amplitude rotations of the electric vector position angle (EVPA), both in the optical and in the X-rays. More often, the polarization degree and the EVPA show a seemingly erratic variability. To interpret these observations, we present a geometric and deterministic model in which off-axis, compact emitting features (i.e.,~blobs) propagate along the jet with the local velocity of the flow. The dynamics of the blobs is determined by the jet electromagnetic fields, which are calculated self-consistently using an analytical model of magnetically dominated outflows. The jet is axisymmetric, and its electromagnetic fields do not have a turbulent component. We show that the observed polarization is sensitive to the initial spatial configurations of the blobs. For the same jet structure, we observe several remarkably complex polarization patterns, including large EVPA rotations of $180^{\circ}$ or more in both directions and more erratic fluctuations. Simultaneous high-amplitude variations of the polarization degree and the EVPA can coincide with peaks of the observed luminosity. However, seemingly uncorrelated variations are also possible. We discuss the feasibility of constraining the particle acceleration mechanism from multifrequency polarimetric observations.

preprint2024arXiv

Nuclear rings are the inner edge of a gap around the Lindblad Resonance

Gaseous nuclear rings are large-scale coherent structures commonly found at the centres of barred galaxies. We propose that they are an accumulation of gas at the inner edge of an extensive gap that forms around the Inner Lindblad Resonance (ILR). The gap initially opens because the bar potential excites strong trailing waves near the ILR, which remove angular momentum from the gas disc and transport the gas inwards. The gap then widens because the bar potential continuously excites trailing waves at the inner edge of the gap, which remove further angular momentum, moving the edge further inwards until it stops at a distance of several wavelengths from the ILR. The gas accumulating at the inner edge of the gap forms the nuclear ring. The speed at which the gap edge moves and its final distance from the ILR strongly depend on the sound speed, explaining the puzzling dependence of the nuclear ring radius on the sound speed in simulations.

preprint2022arXiv

Filamentation of Fast Radio Bursts in magnetar winds

Magnetars are the most promising progenitors of Fast Radio Bursts (FRBs). Strong radio waves propagating through the magnetar wind are subject to non-linear effects, including modulation/filamentation instabilities. We derive the dispersion relation for modulations of strong waves propagating in magnetically-dominated pair plasmas focusing on dimensionless strength parameters $a_0\lesssim 1$, and discuss implications for FRBs. As an effect of the instability, the FRB radiation intensity develops sheets perpendicular to the direction of the wind magnetic field. When the FRB front expands outside the radius where the instability ends, the radiation sheets are scattered due to diffraction. The FRB scattering timescale depends on the properties of the magnetar wind. In a cold wind, the typical scattering timescale is $τ_{\rm sc}\sim{\rm\; μs-ms}$ at the frequency $ν\sim 1{\rm\; GHz}$. The scattering timescale increases at low frequencies, with the scaling $τ_{\rm sc}\proptoν^{-2}$. The frequency-dependent broadening of the brightest pulse of FRB 181112 is consistent with this scaling. From the scattering timescale of the pulse, one can estimate that the wind Lorentz factor is larger than a few tens. In a warm wind, the scattering timescale can approach $τ_{\rm sc}\sim{\rm\; ns}$. Then scattering produces a frequency modulation of the observed intensity with a large bandwidth, $Δν\sim 1/τ_{\rm sc}\gtrsim 100{\rm\; MHz}$. Broadband frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind.

preprint2022arXiv

The Physical Origin and the Properties of Arm Spurs/Feathers in Local Simulations of the Wiggle Instability

Gaseous substructures such as feathers and spurs dot the landscape of spiral arms in disc galaxies. One of the candidates to explain their formation is the wiggle instability of galactic spiral shocks. We study the wiggle instability using local 2D hydrodynamical isothermal non-self gravitating simulations. We find that: (1) Simulations agree with analytic linear stability analysis only under stringent conditions. They display surprisingly strong non-linear coupling between the different modes, even for small mode amplitudes ($\sim 1\%$). (2) We demonstrate that the wiggle instability originates from a combination of two physically distinct mechanisms: the first is the Kelvin-Helmholtz instability, and the second is the amplification of infinitesimal perturbations from repeated shock passages. These two mechanisms can operate simultaneously, and which mechanism dominates depends on the underlying parameters. (3) We explore the parameter space and study the properties of spurs/feathers generated by the wiggle instability. The wiggle instability is highly sensitive to the underlying parameters. The feather separation decreases, and the growth rate increases, with decreasing sound speed, increasing potential strength and decreasing interarm distance. (4) We compare our simulations with a sample of 20 galaxies in the HST Archival Survey of Spiral Arm Substructure of La Vigne et al. and find that the wiggle instability is able to reproduce the typical range of feather separations seen in observations. It remains unclear how the wiggle instability relates to competing mechanisms for spur/feather formation such as the magneto-jeans instability and the stochastic accumulation of gas due to correlated supernova feedback.

preprint2021arXiv

A fully-kinetic model for orphan gamma-ray flares in blazars

Blazars emit a highly-variable non-thermal spectrum. It is usually assumed that the same non-thermal electrons are responsible for the IR-optical-UV emission (via synchrotron) and the gamma-ray emission (via inverse Compton). Hence, the light curves in the two bands should be correlated. Orphan gamma-ray flares (i.e., lacking a luminous low-frequency counterpart) challenge our theoretical understanding of blazars. By means of large-scale two-dimensional radiative particle-in-cell simulations, we show that orphan gamma-ray flares may be a self-consistent by-product of particle energization in turbulent magnetically-dominated pair plasmas. The energized particles produce the gamma-ray flare by inverse Compton scattering an external radiation field, while the synchrotron luminosity is heavily suppressed since the particles are accelerated nearly along the direction of the local magnetic field. The ratio of inverse Compton to synchrotron luminosity is sensitive to the initial strength of turbulent fluctuations (a larger degree of turbulent fluctuations weakens the anisotropy of the energized particles, thus increasing the synchrotron luminosity). Our results show that the anisotropy of the non-thermal particle population is key to modeling the blazar emission.

preprint2020arXiv

Hard Synchrotron Spectra from Magnetically Dominated Plasma Turbulence

Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically-dominated ($σ_0\gg1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim σ_0 γ_{th0}$ (here, $γ_{th0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are roughly isotropic. This energy-dependent anisotropy leads to a synchrotron spectral flux $νF_ν\propto ν^s$ that is much harder than for isotropic particles. Remarkably, for $σ_0\gg1$ we find that the solid-angle-averaged spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim δB_{\rm rms0}^2/B_0^2\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Our findings may help explaining the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of Pulsar Wind Nebulae.

preprint2020arXiv

Magnetic energy dissipation and origin of non-thermal spectra in radiatively efficient relativistic sources

The dissipation of turbulent magnetic fields is an appealing scenario to explain the origin of non-thermal particles in high-energy astrophysical sources. However, it has been suggested that the particle distribution may effectively thermalise when the radiative (synchrotron and/or Inverse Compton) losses are severe. Inspired by recent PIC simulations of relativistic turbulence, which show that electrons are impulsively heated in intermittent current sheets by a strong electric field aligned with the local magnetic field, we instead argue that in plasmas where the particle number density is dominated by the pairs (electron-positron and electron-positron-ion plasmas): (i) as an effect of fast cooling and of different injection times, the electron energy distribution is ${\rm d}n_e/{\rm d}γ\proptoγ^{-2}$ for $γ\lesssimγ_{\rm heat}$ (the Lorentz factor $γ_{\rm heat}$ being close to the equipartition value), while the distribution steepens at higher energies; (ii) since the time scales for the turbulent fields to decay and for the photons to escape are of the same order, the magnetic and the radiation energy densities in the dissipation region are comparable; (iii) if the mass energy of the plasma is dominated by the ion component, the pairs with a Lorentz factor smaller than a critical one (of the order of the proton-to-electron mass ratio) become isotropic, while the pitch angle remains small otherwise. The outlined scenario is consistent with the typical conditions required to reproduce the Spectral Energy Distribution of blazars, and allows one to estimate the magnetisation of the emission site. Finally, we show that turbulence within the Crab Nebula may power the observed gamma-ray flares if the pulsar wind is nearly charge-separated at high latitudes.

preprint2016arXiv

A model for periodic blazars

We describe a scenario to explain blazar periodicities with timescales of $\sim$ few years. The scenario is based on a binary super-massive black hole (SMBH) system in which one of the two SMBH carries a jet. We discuss the various mechanisms that can cause the jet to precess and produce corkscrew patterns through space with a scale of $\sim$ few pc. It turns out that the dominant mechanism responsible for the precession is simply the imprint of the jet-carrying SMBH orbital speed on the jet. Gravitational deflection and Lense-Thirring precession (due to the gravitational field of the other SMBH) are second order effects. We complement the scenario with a kinematical jet model which is inspired to the spine-sheath structure observed in M87. One of the main advantages of such a structure is that it allows the peak of the synchrotron emission to scale with frequency according to $νF\propto ν^ξ$ as the viewing angle is changed, where $ξ$ is not necessarily 3 or 4 as in the case of jets with uniform velocity, but can be $ξ\sim 1$. Finally, we apply the model to the source PG1553+113, which has been recently claimed to show a $T_{\rm obs}=(2.18\pm 0.08)\text{ yr}$ periodicity. We are able to reproduce the optical and gamma-ray light curves and multiple synchrotron spectra simultaneously. We also give estimates of the source mass and size.

preprint2016arXiv

Cross-correlation of the cosmic 21-cm signal and Lyman Alpha Emitters during reionization

Interferometry of the cosmic 21-cm signal is set to revolutionize our understanding of the Epoch of Reionization (EoR), eventually providing 3D maps of the early Universe. Initial detections however will be low signal-to-noise, limited by systematics. To confirm a putative 21-cm detection, and check the accuracy of 21-cm data analysis pipelines, it would be very useful to cross-correlate against a genuine cosmological signal. The most promising cosmological signals are wide-field maps of Lyman alpha emitting galaxies (LAEs), expected from the Subaru Hyper-Suprime Cam (HSC) Ultra-Deep field. Here we present estimates of the correlation between LAE maps at z~7 and the 21-cm signal observed by both the Low Frequency Array (LOFAR) and the planned Square Kilometer Array Phase 1 (SKA1). We adopt a systematic approach, varying both: (i) the prescription of assigning LAEs to host halos; and (ii) the large-scale structure of neutral and ionized regions (i.e. EoR morphology). We find that the LAE-21cm cross-correlation is insensitive to (i), thus making it a robust probe of the EoR. A 1000h observation with LOFAR would be sufficient to discriminate at >1 standard deviation a fully ionized Universe from one with a mean neutral fraction of xHI~0.50, using the LAE-21cm cross-correlation function on scales of R~3-10 Mpc. Unlike LOFAR, whose detection of the LAE-21cm cross-correlation is limited by noise, SKA1 is mostly limited by ignorance of the EoR morphology. However, the planned 100h wide-field SKA1-Low survey will be sufficient to discriminate an ionized Universe from one with xHI~0.25, even with maximally pessimistic assumptions.

preprint2016arXiv

The Evolution of 21-cm Structure (EOS): public, large-scale simulations of Cosmic Dawn and Reionization

We introduce the Evolution of 21-cm Structure (EOS) project: providing periodic, public releases of the latest cosmological 21-cm simulations. 21-cm interferometry is set to revolutionize studies of the Cosmic Dawn (CD) and epoch of reionization (EoR), eventually resulting in 3D maps of the first billion years of our Universe. Progress will depend on sophisticated data analysis pipelines, which are in turn tested on large-scale mock observations. Here we present the 2016 EOS data release, consisting of the largest (1.6 Gpc on side with a 1024^3 grid), public 21-cm simulations of the CD and EoR. We include calibrated, sub-grid prescriptions for inhomogeneous recombinations and photo-heating suppression of star formation in small mass galaxies. We present two simulation runs that approximately bracket the contribution from faint unseen galaxies. From these two extremes, we predict that the duration of reionization (defined as a change in the mean neutral fraction from 0.9 to 0.1) should be between 2.7 < Delta z < 5.7. The large-scale 21-cm power during the advanced EoR stages can be different by up to a factor of ~10, depending on the model. This difference has a comparable contribution from: (i) the typical bias of sources; and (ii) a more efficient negative feedback in models with an extended EoR driven by faint galaxies. We also make detectability forecasts. With a 1000h integration, HERA and SKA1-low should achieve a signal-to-noise of ~few-hundreds throughout the EoR/CD, while in the maximally optimistic scenario of perfect foreground cleaning, all instruments should make a statistical detection of the cosmic signal. We also caution that our ability to clean foregrounds determines the relative performance of narrow/deep vs. wide/shallow surveys expected with SKA1. Our 21-cm power spectra, simulation outputs and visualizations are publicly available.

preprint2015arXiv

Some comments on the electrodynamics of binary pulsars

We consider the electrodynamics of in-spiraling binary pulsars, showing that there are two distinct ways in which they may emit radiation. On the one hand, even if the pulsars do not rotate, we show that in vacuo orbital rotation generates magnetic quadrupole emission, which, in the late stages of the binary evolution becomes nearly as effective as magnetic dipole emission by a millisecond pulsar. On the other hand, we show that interactions of the two magnetic fields generate powerful induction electric fields, which cannot be screened by a suitable distribution of charges and currents like they are in isolated pulsars. We compute approximate electromotive forces for this case.

preprint2015arXiv

The clustering of Lyman alpha emitters at z=7: implications for reionization and host halo masses

The Ly$α$ line of high-redshift galaxies has emerged as a powerful probe of both early galaxy evolution and the epoch of reionization (EoR). Motivated by the upcoming wide-field survey with the Subaru Hyper Suprime-Cam (HSC), we study the angular correlation function (ACF) of narrow-band selected, $z\approx7$ Ly$α$ emitting galaxies (LAEs). The clustering of LAEs is determined by both: (i) their typical host halo masses, $\bar{M}_{\rm h}$; and (ii) the absorption due to a patchy EoR, characterized by an average neutral fraction of the intergalactic medium, $\bar{x}_{\rm HI}$. We bracket the allowed LAE ACFs by exploring extreme scenarios for both the intrinsic Ly$α$ emission and the large-scale pattern (i.e. morphology) of cosmic ionized patches in physical EoR models. Independent of the EoR morphology, current $z\approx 7$ ACF measurements constrain $\bar{x}_{\rm HI}\lesssim 0.5$ (1-$σ$). We also find that the low values of the currently-observed ACF imply that LAEs are hosted by relatively small dark matter halos: $\bar{M}_{\rm h} \lesssim 10^{10} M_{\odot}$, with corresponding duty cycles of $\lesssim$ few per cent. These values are over an order of magnitude lower than the analogous ones for color-selected, Lyman break galaxies, suggesting that $z\approx 7$ narrow-band LAEs searches are preferentially selecting young, star-burst galaxies, residing in less massive halos. The upcoming Ultra Deep campaign with the HSC will significantly improve constraints on both the EoR and LAE host halos.

preprint2014arXiv

Inhomogeneous recombinations during cosmic reionization

By depleting the ionizing photon budget available to expand cosmic HII regions, recombining systems (or Lyman limit systems) can have a large impact during (and following) cosmic reionization. Unfortunately, directly resolving such structures in large-scale reionization simulations is computationally impractical. Instead, here we implement a sub-grid prescription for tracking inhomogeneous recombinations in the intergalactic medium. Building on previous work parameterizing photo-heating feedback on star-formation, we present large-scale, semi-numeric reionization simulations which self-consistently track the local (sub-grid) evolution of both sources and sinks of ionizing photons. Our simple, single-parameter model naturally results in both an extended reionization and a modest, slowly-evolving emissivity, consistent with observations. Recombinations are instrumental in slowing the growth of large HII regions, and damping the rapid rise of the ionizing background in the late stages of (and following) reionization. As a result, typical HII regions are smaller by factors of $\sim 2-3$ throughout reionization. The large-scale ($k\lesssim 0.2\text{ Mpc$^{-1}$}$) ionization power spectrum is suppressed by factors of $\gtrsim 2-3$ in the second half of reionization. Therefore properly modeling recombinations is important in interpreting virtually all reionization observables, including upcoming interferometry with the redshifted 21 cm line. Consistent with previous works, we find the clumping factor of ionized gas to be $C_{\rm HII}\sim 4$ at the end of reionization.

preprint2013arXiv

How does radiative feedback from a UV background impact reionization?

An ionizing UV background (UVB) inhibits gas accretion and photo-evaporates gas from the shallow potential wells of small, dwarf galaxies. During cosmological reionization, this effect can result in negative feedback: suppressing star-formation inside HII regions, thus impeding their continued growth. It is difficult to model this process, given the enormous range of scales involved. We tackle this problem using a tiered approach: combining parameterized results from single-halo collapse simulations with large-scale models of reionization. In the resulting reionization models, the ionizing emissivity of galaxies depends on the local values of the reionization redshift and the UVB intensity. We present a physically-motivated analytic expression for the average minimum mass of star-forming galaxies, which can be readily used in modeling galaxy formation. We find that UVB feedback: (i) delays the end stages of reionization by less than 0.5 in redshift; (ii) results in a more uniform distribution of HII regions, peaked on smaller-scales (with large-scale ionization power suppressed by tens of percent); and (iii) suppresses the global photoionization rate per baryon by a factor of < 2 towards the end of reionization. However, the impact is modest, since the hydrodynamic response of the gas to the UVB occurs on a time-scale comparable to reionization. In particular, the popular approach of modeling UVB feedback with an instantaneous transition in the minimum mass of star-forming galaxies, dramatically overestimates its importance. UVB feedback does not significantly affect reionization unless: (i) molecularly-cooled galaxies contribute significantly to reionization; or (ii) internal feedback processes strongly couple with UVB feedback in the early Universe. Since both are considered unlikely, we conclude that there is no significant self-regulation of reionization by UVB feedback.

preprint2013arXiv

The depletion of gas in high-redshift dwarf galaxies from an inhomogeneous reionization

The reionization of the intergalactic medium (IGM) was likely inhomogeneous and extended. By heating the IGM and photo-evaporating gas from the outskirts of galaxies, this process can have a dramatic impact on the growth of structures. Using a suite of spherically-symmetric collapse simulations spanning a large parameter space, we study the impact of an ionizing ultraviolet background (UVB) on the condensation of baryons onto dark matter halos. We present an expression for the halo baryon fraction, which is an explicit function of: (i) halo mass; (ii) UVB intensity; (iii) redshift; (iv) redshift at which the halo was exposed to a UVB. We also present a corresponding expression for the characteristic or critical mass, defined as the halo mass which retains half of its baryons compared to the global value. Since our results are general and physically-motivated, they can be broadly applied to inhomogeneous reionization models.