-
Simulating short GRB jets in late binary neutron star merger environments
Authors:
Matteo Pais,
Tsvi Piran,
Kenta Kiuchi,
Masaru Shibata
Abstract:
The electromagnetic emission and the afterglow observations of the binary neutron star merger event GW 170817A confirmed the association of the merger with a short gamma-ray burst (sGRB) harboring a narrow ($5$°-$10$°) and powerful ($10^{49}$-$10^{50}~$erg) jet. Using the 1~second-long neutrino-radiation-GR-MHD simulation of coalescing neutron stars of Kiuchi et al. (2023) and following the semi-a…
▽ More
The electromagnetic emission and the afterglow observations of the binary neutron star merger event GW 170817A confirmed the association of the merger with a short gamma-ray burst (sGRB) harboring a narrow ($5$°-$10$°) and powerful ($10^{49}$-$10^{50}~$erg) jet. Using the 1~second-long neutrino-radiation-GR-MHD simulation of coalescing neutron stars of Kiuchi et al. (2023) and following the semi-analytical estimates of Pais et al. (2023), we inject a narrow, powerful, unmagnetized jet into the post-merger phase. We explore different opening angles, luminosities, central engine durations, and times after the merger. We explore early ($0.1~$s following the merger) and late ($1~$s) jet launches; the latter is consistent with the time delay of $\approx 1.74~$s observed between GW 170817 and GRB 170817A. We demonstrate that the semi-analytical estimates correctly predict the jets' breakout and collimation conditions. When comparing our synthetic afterglow light curves to the observed radio data of GW170807, we find a good agreement for a $3 \times 10^{49}$ ergs jet launched late with an opening angle in the range $\simeq 5$°-$7$°.
△ Less
Submitted 26 July, 2024;
originally announced July 2024.
-
The collimation of relativistic jets in post-neutron star binary merger simulations
Authors:
Matteo Pais,
Tsvi Piran,
Yuri Lyubarsky,
Kenta Kiuchi,
Masaru Shibata
Abstract:
The gravitational waves from the binary neutron star merger GW170817 were accompanied by a multi-wavelength electromagnetic counterpart, which confirms the association of the merger with a short gamma-ray burst (sGRB). The afterglow observations implied that the event was accompanied by a narrow, $\sim 5~$deg, and powerful, $\sim 10^{50}$ erg, jet. We study the propagation of a Poynting flux-domin…
▽ More
The gravitational waves from the binary neutron star merger GW170817 were accompanied by a multi-wavelength electromagnetic counterpart, which confirms the association of the merger with a short gamma-ray burst (sGRB). The afterglow observations implied that the event was accompanied by a narrow, $\sim 5~$deg, and powerful, $\sim 10^{50}$ erg, jet. We study the propagation of a Poynting flux-dominated jet within the merger ejecta (kinematic, neutrino-driven and MRI turbulence-driven) of a neutrino-radiation-GR-MHD simulation of two coalescing neutron stars. We find that the presence of a post-merger low-density/low-pressure polar cavity, that arose due to angular momentum conservation, is crucial to let the jet break out. At the same time the ejecta collimates the jet to a narrow opening angle. The collimated jet has a narrow opening angle of $\sim 4$-$7$ deg and an energy of $10^{49}$-$10^{50}~$erg, in line with the observations of GW170817 and other sGRBs.
△ Less
Submitted 13 March, 2023; v1 submitted 16 November, 2022;
originally announced November 2022.
-
The velocity distribution of outflows driven by choked jets in stellar envelopes
Authors:
Matteo Pais,
Tsvi Piran,
Ehud Nakar
Abstract:
Many stripped envelope supernovae (SNe) present a signature of high-velocity material responsible for broad absorption lines in the observed spectrum. These include SNe that are associated with long gamma-ray bursts (LGRBs) and low-luminosity GRBs (llGRBs), and SNe that are not associated with GRBs. Recently it was suggested that this high velocity material originates from a cocoon that is driven…
▽ More
Many stripped envelope supernovae (SNe) present a signature of high-velocity material responsible for broad absorption lines in the observed spectrum. These include SNe that are associated with long gamma-ray bursts (LGRBs) and low-luminosity GRBs (llGRBs), and SNe that are not associated with GRBs. Recently it was suggested that this high velocity material originates from a cocoon that is driven by a relativistic jet. In LGRBs this jet breaks out successfully from the stellar envelope, while in llGRBs and SNe that are not associated with GRBs the jet is choked. Here we use numerical simulations to explore the velocity distribution of an outflow that is driven by a choked jet and its dependence on the jet and progenitor properties. We find that in all cases where the jet is not choked too deep within the star, the outflow carries a roughly constant amount of energy per logarithmic scale of proper velocity over a wide range of velocities, which depends mostly on the cocoon volume at the time of its breakout. This is a universal property of jets driven outflows, which does not exist in outflows of spherically symmetric explosions or when the jets are choked very deep within the star. We therefore conclude that jets that are choked (not too deep) provide a natural explanation to the fast material seen in the early spectra of stripped envelope SNe that are not associated with LGRBs and that properties of this material could reveal information on the otherwise hidden jets.
△ Less
Submitted 8 December, 2022; v1 submitted 30 August, 2022;
originally announced August 2022.
-
Simulating TeV gamma-ray morphologies of shell-type supernova remnants
Authors:
Matteo Pais,
Christoph Pfrommer
Abstract:
Supernova remnant (SNR) shocks provide favourable sites of cosmic ray (CR) proton acceleration if the local magnetic field direction is quasi-parallel to the shock normal. Using the moving-mesh magneto-hydrodynamical (MHD) code AREPO we present a suite of SNR simulations with CR acceleration in the Sedov-Taylor phase that combine different magnetic field topologies, density distributions with grad…
▽ More
Supernova remnant (SNR) shocks provide favourable sites of cosmic ray (CR) proton acceleration if the local magnetic field direction is quasi-parallel to the shock normal. Using the moving-mesh magneto-hydrodynamical (MHD) code AREPO we present a suite of SNR simulations with CR acceleration in the Sedov-Taylor phase that combine different magnetic field topologies, density distributions with gradients and large-scale fluctuations, and -- for our core-collapse SNRs -- a multi-phase interstellar medium with dense clumps with a contrast of $10^4$. Assuming the hadronic gamma-ray emission model for the TeV gamma-ray emission, we find that large-amplitude density fluctuations of $δρ/ρ_0\gtrsim75$ per cent are required to strongly modulate the gamma-ray emissivity in a straw man's model in which the acceleration efficiency is independent of magnetic obliquity. However, this causes strong corrugations of the shock surface that are ruled out by gamma-ray observations. By contrast, magnetic obliquity-dependent acceleration can easily explain the observed variance in gamma-ray morphologies ranging from SN1006 (with a homogeneous magnetic field) to Vela Junior and RX J1713 (with a turbulent field) in a single model that derives from plasma particle-in-cell simulations. Our best-fit model for SN1006 has a large-scale density gradient of $ \nabla{n}\simeq 0.0034~\mathrm{cm}^{-3}~\mathrm{pc}^{-1} $ pointing from south-west to north-east and a magnetic inclination with the plane of the sky of $\lesssim10^\circ$. Our best-fit model for Vela Junior and RX J1713 adopts a combination of turbulent magnetic field and dense clumps to explain their TeV gamma-ray morphologies and moderate shock corrugations.
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
Evolution and observational signatures of the cosmic ray electron spectrum in SN 1006
Authors:
Georg Winner,
Christoph Pfrommer,
Philipp Girichidis,
Maria Werhahn,
Matteo Pais
Abstract:
Supernova remnants (SNRs) are believed to be the source of Galactic cosmic rays (CRs). SNR shocks accelerate CR protons and electrons which reveal key insights into the non-thermal physics by means of their synchrotron and $γ$-ray emission. The remnant SN 1006 is an ideal particle acceleration laboratory because it is observed across all electromagnetic wavelengths from radio to $γ$-rays. We perfo…
▽ More
Supernova remnants (SNRs) are believed to be the source of Galactic cosmic rays (CRs). SNR shocks accelerate CR protons and electrons which reveal key insights into the non-thermal physics by means of their synchrotron and $γ$-ray emission. The remnant SN 1006 is an ideal particle acceleration laboratory because it is observed across all electromagnetic wavelengths from radio to $γ$-rays. We perform three-dimensional (3D) magnetohydrodynamics (MHD) simulations where we include CR protons and follow the CR electron spectrum. By matching the observed morphology and non-thermal spectrum of SN 1006 in radio, X-rays and $γ$-rays, we gain new insight into CR electron acceleration and magnetic field amplification. 1. We show that a mixed leptonic-hadronic model is responsible for the $γ$-ray radiation: while leptonic inverse-Compton emission and hadronic pion-decay emission contribute equally at GeV energies observed by Fermi, TeV energies observed by imaging air Cherenkov telescopes are hadronically dominated. 2. We show that quasi-parallel acceleration (i.e., when the shock propagates at a narrow angle to the upstream magnetic field) is preferred for CR electrons and that the electron acceleration efficiency of radio-emitting GeV electrons at quasi-perpendicular shocks is suppressed at least by a factor ten. This precludes extrapolation of current one-dimensional plasma particle-in-cell simulations of shock acceleration to realistic SNR conditions. 3. To match the radial emission profiles and the $γ$-ray spectrum, we require a volume-filling, turbulently amplified magnetic field and that the Bell-amplified magnetic field is damped in the immediate post-shock region. Our work connects micro-scale plasma physics simulations to the scale of SNRs.
△ Less
Submitted 5 December, 2020; v1 submitted 11 June, 2020;
originally announced June 2020.
-
Constraining the coherence scale of the interstellar magnetic field using TeV gamma-ray observations of supernova remnants
Authors:
Matteo Pais,
Christoph Pfrommer,
Kristian Ehlert,
Maria Werhahn,
Georg Winner
Abstract:
Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ra…
▽ More
Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ray morphology of SNRs such as Vela Jr. and SN1006 solely by varying the magnetic morphology. This implies that gamma-ray bright regions result from quasi-parallel shocks (i.e., when the shock propagates at a narrow angle to the upstream magnetic field), which are known to efficiently accelerate CR protons, and that gamma-ray dark regions point to quasi-perpendicular shock configurations. Comparison of the simulated gamma-ray morphology to observations allows us to constrain the magnetic coherence scale $λ_B$ around Vela Jr. and SN1006 to $λ_B \simeq 13_{-4.3}^{+13}$ pc and $λ_B >200_{-10}^{+80}$ pc, respectively, where the ambient magnetic field of SN1006 is consistent with being largely homogeneous. We find consistent pure hadronic and mixed hadronic-leptonic models that both reproduce the multi-frequency spectra from the radio to TeV gamma rays and match the observed gamma-ray morphology. Finally, to capture the propagation of a SNR shock in a clumpy interstellar medium, we study the interaction of a shock with a dense cloud with numerical simulations and analytics. We construct an analytical gamma-ray model for a core collapse SNR propagating through a structured interstellar medium, and show that the gamma-ray luminosity is only biased by 30% for realistic parameters.
△ Less
Submitted 9 June, 2020; v1 submitted 8 May, 2018;
originally announced May 2018.
-
The effect of cosmic-ray acceleration on supernova blast wave dynamics
Authors:
M. Pais,
C. Pfrommer,
K. Ehlert,
R. Pakmor
Abstract:
Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh c…
▽ More
Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code {\sc arepo} for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.
△ Less
Submitted 30 April, 2018;
originally announced May 2018.
-
TRAGALDABAS. First results on cosmic ray studies and their relation with the solar activity, the Earth magnetic field and the atmospheric properties
Authors:
J. A. Garzón,
J. Collazo,
J. Cuenca-García,
D. Garcí a Castro,
J. Otero,
M. Yermo,
J. J. Blanco,
T. Kurtukian,
A. Morozova,
M. A. Pais,
A. Blanco,
P. Fonte,
L. Lopes,
G. Kornakov,
H. Álvarez-Pol,
P. Cabanelas,
A. Pazos,
M. Seco,
I. Riádigos,
V. Pérez Muñuzuri,
A. Gómez-Tato,
J. C. Mouriño,
P. Rey,
J. Taboada
Abstract:
Cosmic rays originating from extraterrestrial sources are permanently arriving at Earth atmosphere, where they produce up to billions of secondary particles. The analysis of the secondary particles reaching to the surface of the Earth may provide a very valuable information about the Sun activity, changes in the geomagnetic field and the atmosphere, among others. In this article, we present the fi…
▽ More
Cosmic rays originating from extraterrestrial sources are permanently arriving at Earth atmosphere, where they produce up to billions of secondary particles. The analysis of the secondary particles reaching to the surface of the Earth may provide a very valuable information about the Sun activity, changes in the geomagnetic field and the atmosphere, among others. In this article, we present the first preliminary results of the analysis of the cosmic rays measured with a high resolution tracking detector, TRAGALDABAS, located at the Univ. of Santiago de Compostela, in Spain.
△ Less
Submitted 25 January, 2017;
originally announced January 2017.
-
Temporal Evolution of Sunspot Areas and Estimation of Related Plasma Flows
Authors:
R. Gafeira,
C. C. Fonte,
M. A. Pais,
J. Fernandes
Abstract:
The increased amount of information provided by ongoing missions such as the Solar Dynamics Observatory (SDO) represents a great challenge for the understanding of basic questions such as the internal structure of sunspots and how they evolve with time. Here, we contribute with the exploitation of new data, to provide a better understanding of the separate growth and decay of sunspots, umbra and p…
▽ More
The increased amount of information provided by ongoing missions such as the Solar Dynamics Observatory (SDO) represents a great challenge for the understanding of basic questions such as the internal structure of sunspots and how they evolve with time. Here, we contribute with the exploitation of new data, to provide a better understanding of the separate growth and decay of sunspots, umbra and penumbra. Using fuzzy sets to compute separately the areas of sunspot umbra and penumbra, the growth and decay rates for active regions NOAA 11117, NOAA 11428, NOAA 11429, and NOAA 11430 are computed, from the analysis of intensitygrams obtained by the Helioseismic and Magnetic Imager onboard SDO. A simplified numerical model is proposed for the decay phase, whereby an empirical irrotational and uniformly convergent horizontal velocity field interacting with an axially symmetric and height invariant magnetic field, reproduces the large-scale features of the much more complex convection observed inside sunspots.
△ Less
Submitted 27 October, 2013;
originally announced October 2013.
-
Variation rate of sunspot area
Authors:
R. Gafeira,
C. C. Fonte,
M. A. Pais,
J. Fernandes
Abstract:
The emergence of the magnetic field through the photosphere has multiple manifestations and sunspots are the most prominent examples of this. One of the most relevant sunspot properties, to study both its structure and evolution, is the sunspot area: either total, umbra or penumbra area. Recently Schlichenmaier et al. (2010) studied the evolution of the active region (AR) NOAA 11024 concluding tha…
▽ More
The emergence of the magnetic field through the photosphere has multiple manifestations and sunspots are the most prominent examples of this. One of the most relevant sunspot properties, to study both its structure and evolution, is the sunspot area: either total, umbra or penumbra area. Recently Schlichenmaier et al. (2010) studied the evolution of the active region (AR) NOAA 11024 concluding that during the penumbra formation the umbra area remains constant and that the increase of the total sunspot area is caused exclusively by the penumbra growth. In this presentation the Schlichenmaier's conclusion is firstly tested, investigating the evolution of four different ARs. Hundreds of Intensitygram images from the Helioseismic and Magnetic Imager (HMI) images are used, obtained by the Solar Dynamics Observatory, in order to describe the area evolution of the above ARs and estimate the increase and decrease rates for umbra and penumbra areas, separately. A simple magnetohydrodynamic model is then tentatively used in a first approximation to explain the observed results.
△ Less
Submitted 29 November, 2012;
originally announced November 2012.
-
Long-time temperature variations in Portugal over the last 140 years and the effect of the solar activity
Authors:
A. L. Morozova,
P. Thejll,
M. A. Pais
Abstract:
We present the analysis of temperature variations in Portugal for 140 years (from 1865 to 2005). The two stations with the longest time series of temperature measurements (Lisbon and Coimbra) have been used to study the dependence of the portuguese climate variations on the changes of both global circulation and solar activity. Monthly averaged temperature series have been analyzed together with m…
▽ More
We present the analysis of temperature variations in Portugal for 140 years (from 1865 to 2005). The two stations with the longest time series of temperature measurements (Lisbon and Coimbra) have been used to study the dependence of the portuguese climate variations on the changes of both global circulation and solar activity. Monthly averaged temperature series have been analyzed together with monthly North-Atlantic Oscillation index data, sunspot numbers and cosmic ray flux intensity. Different statistical methods (the correlation analysis and the multiple regression analysis) were used. Our results show that the temperature in Portugal depends not only on the atmospheric variations in the studied region but also on the variations of the solar parameters such as sunspot numbers and the cosmic rays flux intensity. Also, the dependence of temperature on solar parameters is strong during the cold season (November to February), while much weaker during the warm season. For some months, strong correlations between the temperature variations of the current month and the North-Atlantic Oscillation index values of the previous month have been found. The correlation between climatic and solar parameters shows up best on the decadal and decadal-to-centennial timescale. It is found that the temperature correlates positively with the sunspot numbers and negatively with the CR flux intensity throughout the year. Besides, the absolute values of the correlation coefficients between the temperature and the CR are higher than those between the temperature and the sunspot numbers. Our results are consistent with some of the proposed mechanisms that relate solar activity to Earth climate and could be explained through the effect of the solar UV radiation and stratosphere-troposphere coupling or/and through the effect of the CR particles on clouds and stratospheric and tropospheric conditions.
△ Less
Submitted 3 September, 2010;
originally announced September 2010.