Showing 1–50 of 122 results for author: Peres, G

Unraveling the Effects of Dense Medium on a Near to Bohm-Limit Acceleration in Kepler's SNR

Authors: Vincenzo Sapienza, Marco Miceli, Oleh Petruk, Aya Bamba, Salvatore Orlando, Fabrizio Bocchino, Giovanni Peres

Abstract: The maximum energy of electrons accelerated by supernova remnants (SNR) is typically limited by radiative losses. In this scenario, the synchrotron cooling time scale is equal to the acceleration time scale. On the other hand, the low propagation speed of a shock in a dense medium is expected to result in an extended acceleration time scale, thus inducing a decrease in the maximum electron energy… ▽ More The maximum energy of electrons accelerated by supernova remnants (SNR) is typically limited by radiative losses. In this scenario, the synchrotron cooling time scale is equal to the acceleration time scale. On the other hand, the low propagation speed of a shock in a dense medium is expected to result in an extended acceleration time scale, thus inducing a decrease in the maximum electron energy for a given SNR age and in the X-ray nonthermal flux. The young Kepler's SNR shows an enhanced efficiency of the acceleration process, which is close to the Bohm limit in the north of its shell, where the shock is slowed down by a dense circumstellar medium. Conversely, in the south, where no interaction with a dense medium is evident and the shock speed is high, the acceleration proceeds with a higher Bohm factor. To investigate this scenario, we studied the temporal evolution of the non-thermal emission, taking advantage of two Chandra X-ray observations of Kepler's SNR (performed in 2006 and 2014). We analyzed the spectra of different filaments both in the north and south of the shell, and measured their proper motion. We found a region with low shock velocity where we measured a significant decrease in flux from 2006 to 2014. This could be the first evidence of fading synchrotron emission in Kepler's SNR. This result suggests that under a certain threshold of shock speed the acceleration process could exit the loss-limited regime. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 8 pages, 2 figures, 2 tables

Report number: PoS(ICRC2023)843

Spatial and dynamical structure of the NGC 2264 star-forming region

Authors: E. Flaccomio, G. Micela, G. Peres, S. Sciortino, E. Salvaggio, L. Prisinzano, M. G. Guarcello, L. Venuti, R. Bonito, I. Pillitteri

Abstract: The processes that lead from molecular clouds to protostars, stellar associations, and the evolution of both, are active research topics. The target of this study, NGC2264, is a benchmark star-forming region (SFR) in which these issues can be profitably studied. We assembled a new extensive sample of NGC 2264 members. To this end we used new X-ray data obtained with the XMM-Newton telescope, GAI… ▽ More The processes that lead from molecular clouds to protostars, stellar associations, and the evolution of both, are active research topics. The target of this study, NGC2264, is a benchmark star-forming region (SFR) in which these issues can be profitably studied. We assembled a new extensive sample of NGC 2264 members. To this end we used new X-ray data obtained with the XMM-Newton telescope, GAIA eDR3 data, and an extensive collection of public and published catalogs. Following a previous suggestion that the SFR might extend significantly beyond the better studied areas, our search covers a wide 2.5x2.5 degrees region in the sky. Our catalog comprises more than 2200 candidate members, a ~100% increase over previous determinations. We analyze their spatial distribution and define new substructures. Using GAIA parallaxes we estimate a new average distance to NGC2264 of 722+/-2pc and suggest that the embedded Spokes subregion is ~20pc farther away within the molecular cloud. A complex dynamics is unveiled by the available proper motions and radial velocities: we observe signs of global expansion and rotation. At the same time, we observe the collapse and coalescence of two substructures in a region where active star formation is taking place. The fraction of stars with disks and of those undergoing circumstellar accretion varies significantly across the field, suggesting that star formation has been occurring for several million years. A particularly low disk fraction around the O VII star S Mon might be attributed to external disk photoevaporation or to an older age of the stars in the region. NGC2264 is not dynamically relaxed and its present configuration is the result of multiple dynamical processes. The cloud has been forming stars for several million years and we identify the process that is likely responsible for the ongoing formation activity. △ Less

Submitted 2 December, 2022; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: 24 pages, 18 figures, accepted for publication on Astronomy & Astrophysics

Journal ref: A&A 670, A37 (2023)

The supernova remnant SN 1006 as a Galactic particle accelerator

Authors: Roberta Giuffrida, Marco Miceli, Damiano Caprioli, Anne Decourchelle, Jacco Vink, Salvatore Orlando, Fabrizio Bocchino, Emanuele Greco, Giovanni Peres

Abstract: The origin of cosmic rays is a pivotal open issue of high-energy astrophysics. Supernova remnants are strong candidates to be the Galactic factory of cosmic rays, their blast waves being powerful particle accelerators. However, supernova remnants can power the observed flux of cosmic rays only if they transfer a significant fraction of their kinetic energy to the accelerated particles, but conclus… ▽ More The origin of cosmic rays is a pivotal open issue of high-energy astrophysics. Supernova remnants are strong candidates to be the Galactic factory of cosmic rays, their blast waves being powerful particle accelerators. However, supernova remnants can power the observed flux of cosmic rays only if they transfer a significant fraction of their kinetic energy to the accelerated particles, but conclusive evidence for such efficient acceleration is still lacking. In this scenario, the shock energy channeled to cosmic rays should induce a higher post-shock density than that predicted by standard shock conditions. Here we show this effect, and probe its dependence on the orientation of the ambient magnetic field, by analyzing deep X-ray observations of the Galactic remnant of SN 1006. By comparing our results with state-of-the-art models, we conclude that SN 1006 is an efficient source of cosmic rays and obtain an observational support for the quasi-parallel acceleration mechanism. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Published in Nature Communications. The published version of the paper, including supplementary material, is freely available online here https://doi.org/10.1038/s41467-022-32781-4 or as a PDF here https://www.nature.com/articles/s41467-022-32781-4.pdf

A spatially resolved study of hard X-ray emission in Kepler's SNR: indications of different regimes of particle acceleration

Authors: Vincenzo Sapienza, Marco Miceli, Aya Bamba, Satoru Katsuda, Tsutomu Nagayoshi, Yukikatsu Terada, Fabrizio Bocchino, Salvatore Orlando, Giovanni Peres

Abstract: Synchrotron X-ray emission in young supernova remnants (SNRs) is a powerful diagnostic tool to study the population of high energy electrons accelerated at the shock front and the acceleration process. We performed a spatially resolved spectral analysis of NuSTAR and XMM-Newton observations of the young Kepler's SNR, aiming to study in detail its non-thermal emission in hard X-rays. We selected a… ▽ More Synchrotron X-ray emission in young supernova remnants (SNRs) is a powerful diagnostic tool to study the population of high energy electrons accelerated at the shock front and the acceleration process. We performed a spatially resolved spectral analysis of NuSTAR and XMM-Newton observations of the young Kepler's SNR, aiming to study in detail its non-thermal emission in hard X-rays. We selected a set of regions all around the rim of the shell and extracted the corresponding spectra. The spectra were analyzed by adopting a model of synchrotron radiation in the loss-limited regime, to constrain the dependence of the cutoff energy of the synchrotron radiation on the shock velocity. We identify two different regimes of particle acceleration, characterized by different Bohm factors. In the north, where the shock interacts with a dense circumstellar medium (CSM), we found a more efficient acceleration than in the south, where the shock velocity is higher and there are no signs of shock interaction with dense CSM. Our results suggest an enhanced efficiency of the acceleration process in regions where the shock-CSM interaction generates an amplified and turbulent magnetic field. By combining hard X-ray spectra with radio and $γ-$ray observations of Kepler's SNR, we modelled the spectral energy distribution. In the light of our results we propose that the observed $γ-$ray emission is mainly hadronic, and originates in the northern part of the shell. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 16 pages, 7 figures, 5 tables, Accepted for publication in ApJ

Additional evidence for a pulsar wind nebula in the heart of SN 1987A from multi-epoch X-ray data and MHD modeling

Authors: Emanuele Greco, Marco Miceli, Salvatore Orlando, Barbara Olmi, Fabrizio Bocchino, Shigehiro Nagataki, Lei Sun, Jacco Vink, Vincenzo Sapienza, Masaomi Ono, Akira Dohi, Giovanni Peres

Abstract: Since the day of its explosion, supernova (SN) 1987A has been closely monitored to study its evolution and to detect its central compact relic. In fact, the formation of a neutron star is strongly supported by the detection of neutrinos from the SN. However, besides the detection in the Atacama Large Millimeter/submillimeter Array (ALMA) data of a feature that is compatible with the emission arisi… ▽ More Since the day of its explosion, supernova (SN) 1987A has been closely monitored to study its evolution and to detect its central compact relic. In fact, the formation of a neutron star is strongly supported by the detection of neutrinos from the SN. However, besides the detection in the Atacama Large Millimeter/submillimeter Array (ALMA) data of a feature that is compatible with the emission arising from a proto-pulsar wind nebula (PWN), the only hint for the existence of such elusive compact object is provided by the detection of hard emission in NuSTAR data up to ~ 20 keV. We report on the simultaneous analysis of multi-epoch observations of SN 1987A performed with Chandra, XMM-Newton and NuSTAR. We also compare the observations with a state-of-the-art 3D magnetohydrodynamic (MHD) simulation of SN 1987A. A heavily absorbed power-law, consistent with the emission from a PWN embedded in the heart of SN 1987A, is needed to properly describe the high-energy part of the observed spectra. The spectral parameters of the best-fit power-law are in agreement with the previous estimate, and exclude diffusive shock acceleration as a possible mechanism responsible for the observed non-thermal emission. The information extracted from our analysis are used to infer the physical characteristics of the pulsar and the broad-band emission of its nebula, in agreement with the ALMA data. Analysis of the synthetic spectra also show that, in the near future, the main contribution to Fe K emission line will originate in the outermost shocked ejecta of SN 1987A. △ Less

Submitted 14 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ

Evidence for past interaction with an asymmetric circumstellar shell in the young SNR Cassiopeia A

Authors: S. Orlando, A. Wongwathanarat, H. -T. Janka, M. Miceli, S. Nagataki, M. Ono, F. Bocchino, J. Vink, D. Milisavljevic, D. J. Patnaude, G. Peres

Abstract: Observations of the SNR Cassiopeia A (Cas A) show asymmetries in the reverse shock that cannot be explained by models describing a remnant expanding through a spherically symmetric wind of the progenitor star. We investigate whether a past interaction of Cas A with an asymmetric circumstellar shell can account for the observed asymmetries. We performed 3D MHD simulations that describe the remnant… ▽ More Observations of the SNR Cassiopeia A (Cas A) show asymmetries in the reverse shock that cannot be explained by models describing a remnant expanding through a spherically symmetric wind of the progenitor star. We investigate whether a past interaction of Cas A with an asymmetric circumstellar shell can account for the observed asymmetries. We performed 3D MHD simulations that describe the remnant evolution from the SN to its interaction with a circumstellar shell. The initial conditions are provided by a 3D neutrino-driven SN model whose morphology resembles Cas A. We explored the parameter space of the shell, searching for a set of parameters able to produce reverse shock asymmetries at the age of 350 years analogous to those observed in Cas A. The interaction of the remnant with the shell can produce asymmetries resembling those observed in the reverse shock if the shell was asymmetric with the densest portion in the nearside to the northwest (NW). The reverse shock shows the following asymmetries at the age of Cas A: i) it moves inward in the observer frame in the NW region, while it moves outward in other regions; ii) the geometric center of the reverse shock is offset to the NW from the geometric center of the forward shock; iii) the reverse shock in the NW region has enhanced nonthermal emission because, there, the ejecta enter the reverse shock with a higher velocity (between 4000 and 7000 km/s) than in other regions (below 2000 km/s). The asymmetries observed in the reverse shock of Cas A can be interpreted as signatures of the interaction of the remnant with an asymmetric circumstellar shell that occurred between 180 and 240 years after the SN event. We suggest that the shell was, most likely, the result of a massive eruption from the progenitor star that occurred between $10^4$ and $10^5$ years prior to core-collapse. We estimate a total mass of the shell of the order 2 Msun. △ Less

Submitted 12 August, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 20 pages, 11 figures, 2 Tables, 1 interactive graphics. Accepted for publication on A&A, high quality figures will be available on the published paper

Journal ref: A&A 666, A2 (2022)

Negative and Positive Feedback from a Supernova Remnant with SHREC: A detailed Study of the Shocked Gas in IC443

Authors: G. Cosentino, I. Jiménez-Serra, J. C. Tan, J. D. Henshaw, A. T. Barnes, C. -Y. Law, S. Zeng, F. Fontani, P. Caselli, S. Viti, S. Zahorecz, F. Rico-Villas, A. Megías, M. Miceli, S. Orlando, S. Ustamujic, E. Greco, G. Peres, F. Bocchino, R. Fedriani, P. Gorai, L. Testi, J. Martín-Pintado

Abstract: Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is… ▽ More Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions. △ Less

Submitted 9 January, 2022; originally announced January 2022.

Comments: Accepted for Publication in MNRAS; 11 pages, 7 figures

Correcting the effect of stellar spots on ARIEL transmission spectra II.The limb darkening effect

Authors: G. Cracchiolo, G. Micela, G. Morello, G. Peres

Abstract: This paper is part of an effort to correct the transmission spectra of a transiting planet orbiting an active star. In Paper I (Cracchiolo et al. 2020) we have demonstrated a methodology to minimize the potential bias induced by unocculted star spots on the transmission spectrum, assuming a spot model parameterized by filling factor and temperature. In this work we introduce the limb darkening eff… ▽ More This paper is part of an effort to correct the transmission spectra of a transiting planet orbiting an active star. In Paper I (Cracchiolo et al. 2020) we have demonstrated a methodology to minimize the potential bias induced by unocculted star spots on the transmission spectrum, assuming a spot model parameterized by filling factor and temperature. In this work we introduce the limb darkening effect, therefore the position of the spot in the stellar disk and the impact parameter of the transiting planet now play a key role. The method is tested on simulations of planetary transits of three representative kinds of planetary systems, at ARIEL resolution. We find that a realistic treatment of the limb darkening is required to reliably estimate both the spots parameters and the transmission spectrum of the transiting planet. Furthermore, we show that the influence of the spots onthe retrieval of the planetary transmission spectrum is significant for spots close to the center of the star, covering a fraction greater than 0.05 and with a temperature contrast greater than 500 K, and that for these cases our method can confidently extract the transmission spectrum and the impact parameter of the transiting planet for both cases of occulted and not occulted spots, provided that we have an accurate characterization of the stellar parameters and a reliable simulator of the instrument performances. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Comments: 14 pages, 9 figures

ACM Class: F.2.2; I.2.7

Journal ref: MNRAS, 2021

Modeling the Remnants of Core-collapse Supernovae from Luminous Blue Variable stars

Authors: S. Ustamujic, S. Orlando, M. Miceli, F. Bocchino, M. Limongi, A. Chieffi, C. Trigilio, G. Umana, F. Bufano, A. Ingallinera, G. Peres

Abstract: LBVs are massive evolved stars that suffer sporadic and violent mass-loss events. They have been proposed as the progenitors of some core-collapse SNe, but this idea is still debated due to the lack of direct evidence. Since SNRs can carry in their morphology the fingerprints of the progenitor stars as well as of the inhomogeneous CSM sculpted by the progenitors, the study of SNRs from LBVs could… ▽ More LBVs are massive evolved stars that suffer sporadic and violent mass-loss events. They have been proposed as the progenitors of some core-collapse SNe, but this idea is still debated due to the lack of direct evidence. Since SNRs can carry in their morphology the fingerprints of the progenitor stars as well as of the inhomogeneous CSM sculpted by the progenitors, the study of SNRs from LBVs could help to place core-collapse SNe in context with the evolution of massive stars. We investigate the physical, chemical and morphological properties of the remnants of SNe originating from LBVs, in order to search for signatures, revealing the nature of the progenitors, in the ejecta distribution and morphology of the remnants. As a template of LBVs, we considered the actual LBV candidate Gal 026.47+0.02. We selected a grid of models, which describe the evolution of a massive star with properties consistent with those of Gal 026.47+0.02 and its final fate as core-collapse SN. We developed a 3D HD model that follows the post-explosion evolution of the ejecta from the breakout of the shock wave at the stellar surface to the interaction of the SNR with a CSM characterized by two dense nested toroidal shells, parametrized in agreement with multi-wavelength observations of Gal 026.47+0.02. Our models show a strong interaction of the blast wave with the CSM which determines an important slowdown of the expansion of the ejecta in the equatorial plane where the two shells lay, determining a high degree of asymmetry in the remnant. After 10000 years of evolution the ejecta show an elongated shape forming a broad jet-like structure caused by the interaction with the shells and oriented along the axis of the toroidal shells. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 18 pages, 12 figures; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 654, A167 (2021)

Searching for solar KDAR with DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti, M. P. Andrews , et al. (1157 additional authors not shown)

Abstract: The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search.… ▽ More The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search. In this work, we evaluate the proposed KDAR neutrino search strategies by realistically modeling both neutrino-nucleus interactions and the response of DUNE. We find that, although reconstruction of the neutrino energy and direction is difficult with current techniques in the relevant energy range, the superb energy resolution, angular resolution, and particle identification offered by DUNE can still permit great signal/background discrimination. Moreover, there are non-standard scenarios in which searches at DUNE for KDAR in the Sun can probe dark matter interactions. △ Less

Submitted 26 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 19 pages, 13 figures

Report number: FERMILAB-PUB-21-322-LBNF-ND

Journal ref: JCAP10(2021)065

Ariel: Enabling planetary science across light-years

Authors: Giovanna Tinetti, Paul Eccleston, Carole Haswell, Pierre-Olivier Lagage, Jérémy Leconte, Theresa Lüftinger, Giusi Micela, Michel Min, Göran Pilbratt, Ludovic Puig, Mark Swain, Leonardo Testi, Diego Turrini, Bart Vandenbussche, Maria Rosa Zapatero Osorio, Anna Aret, Jean-Philippe Beaulieu, Lars Buchhave, Martin Ferus, Matt Griffin, Manuel Guedel, Paul Hartogh, Pedro Machado, Giuseppe Malaguti, Enric Pallé , et al. (293 additional authors not shown)

Abstract: Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths.… ▽ More Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

Report number: ESA/SCI(2020)1

X-ray emitting structures in the Vela SNR: ejecta anisotropies and progenitor stellar wind residuals

Authors: Vincenzo Sapienza, Marco Miceli, Giovanni Peres, Fabrizio Bocchino, Salvatore Orlando, Emanuele Greco, Jorge Ariel Combi, Federico García, Manami Sasaki

Abstract: The Vela supernova remnant (SNR) shows several ejecta fragments protruding beyond the forward shock (shrapnel). Recent studies have revealed high Si abundance in two shrapnel (A and G), located in opposite directions with respect to the SNR center. This suggests the possible existence of a Si-rich jet-counterjet structure. We analyzed an XMM-Newton observation of a bright clump, behind shrapnel G,… ▽ More The Vela supernova remnant (SNR) shows several ejecta fragments protruding beyond the forward shock (shrapnel). Recent studies have revealed high Si abundance in two shrapnel (A and G), located in opposite directions with respect to the SNR center. This suggests the possible existence of a Si-rich jet-counterjet structure. We analyzed an XMM-Newton observation of a bright clump, behind shrapnel G, which lies along the direction connecting A and G. The aim is to study the physical and chemical properties of this clump to ascertain whether it is part of this putative jet-like structure. We produced background-corrected and adaptively-smoothed count-rate images and median photon energy maps, and performed a spatially resolved spectral analysis. We identified two structures with different physical properties. The first one is remarkably elongated along the direction connecting A and G. Its X-ray spectrum is much softer than that of the other two shrapnel, to the point of hindering the determination of the Si abundance, however its physical and chemical properties are consistent with those of shrapnel A and G. The second structure, running along the southeast-northwest direction, has a higher temperature and appears like a thin filament. By analyzing the ROSAT data, we have found that this filament is part of a very large and coherent structure that we identified in the western rim of the shell. We obtained a thorough description of the tail of Shrapnel G. In addition we discovered a coherent and very extended feature that we interpret as a signature of an earlier interaction of the remnant with the stellar wind of its progenitor star. The peculiar Ne/O ratio we found in the wind residual may be suggestive of a Wolf-Rayet progenitor for Vela SNR, though further analysis is required to address this point. △ Less

Submitted 23 March, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 10 pages, 13 figures, Accepted to A&A

Journal ref: A&A 649, A56 (2021)

Indication of a Pulsar Wind Nebula in the hard X-ray emission from SN 1987A

Authors: Emanuele Greco, Marco Miceli, Salvatore Orlando, Barbara Olmi, Fabrizio Bocchino, Shigehiro Nagataki, Masaomi Ono, Akira Dohi, Giovanni Peres

Abstract: Since the day of its explosion, SN 1987A (SN87A) was closely monitored with the aim to study its evolution and to detect its central compact relic. The detection of neutrinos from the supernova strongly supports the formation of a neutron star (NS). However, the constant and fruitless search for this object has led to different hypotheses on its nature. Up to date, the detection in the ALMA data o… ▽ More Since the day of its explosion, SN 1987A (SN87A) was closely monitored with the aim to study its evolution and to detect its central compact relic. The detection of neutrinos from the supernova strongly supports the formation of a neutron star (NS). However, the constant and fruitless search for this object has led to different hypotheses on its nature. Up to date, the detection in the ALMA data of a feature somehow compatible with the emission arising from a proto Pulsar Wind Nebula (PWN) is the only hint of the existence of such elusive compact object. Here we tackle this 33-years old issue by analyzing archived observations of SN87A performed Chandra and NuSTAR in different years. We firmly detect nonthermal emission in the $10-20$ kev energy band, due to synchrotron radiation. The possible physical mechanism powering such emission is twofold: diffusive shock acceleration (DSA) or emission arising from an absorbed PWN. By relating a state-of-the-art magneto-hydrodynamic simulation of SN87A to the actual data, we reconstruct the absorption pattern of the PWN embedded in the remnant and surrounded by cold ejecta. We found that, even though the DSA scenario cannot be firmly excluded, the most likely scenario that well explains the data is the PWN emission. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: Accepted for publication in ApJ Letters

Linking core-collapse supernova explosions to supernova remnants through 3D MHD modeling

Authors: S. Orlando, A. Wongwathanarat, H. -T. Janka, M. Miceli, M. Ono, S. Nagataki, F. Bocchino, G. Peres

Abstract: The structure and morphology of supernova remnants (SNRs) reflect the properties of the parent supernovae (SNe) and the characteristics of the inhomogeneous environments through which the remnants expand. Linking the morphology of SNRs to anisotropies developed in their parent SNe can be essential to obtain key information on many aspects of the explosion processes associated with SNe. Nowadays, o… ▽ More The structure and morphology of supernova remnants (SNRs) reflect the properties of the parent supernovae (SNe) and the characteristics of the inhomogeneous environments through which the remnants expand. Linking the morphology of SNRs to anisotropies developed in their parent SNe can be essential to obtain key information on many aspects of the explosion processes associated with SNe. Nowadays, our capability to study the SN-SNR connection has been largely improved thanks to multi-dimensional models describing the long-term evolution from the SN to the SNR as well as to observational data of growing quality and quantity across the electromagnetic spectrum which allow to constrain the models. Here we used the numerical resources obtained in the framework of the "Accordo Quadro INAF-CINECA (2017)" together with a CINECA ISCRA Award N.HP10BARP6Y to describe the full evolution of a SNR from the core-collapse to the full-fledged SNR at the age of 2000 years. Our simulations were compared with observations of SNR Cassiopeia A (Cas A) at the age of $\sim 350$~years. Thanks to these simulations we were able to link the physical, chemical and morphological properties of a SNR to the physical processes governing the complex phases of the SN explosion. △ Less

Submitted 25 December, 2020; originally announced December 2020.

Comments: 4 pages, 2 figures, report on the numerical resources obtained in the framework of the "Accordo Quadro INAF-CINECA (2017)". arXiv admin note: text overlap with arXiv:2009.01789

Modeling particle acceleration and non-thermal emission in supernova remnants

Authors: S. Orlando, M. Miceli, S. Ustamujic, A. Tutone, E. Greco, O. Petruk, F. Bocchino, G. Peres

Abstract: According to the most popular model for the origin of cosmic rays (CRs), supernova remnants (SNRs) are the site where CRs are accelerated. Observations across the electromagnetic spectrum support this picture through the detection of non-thermal emission that is compatible with being synchrotron or inverse Compton radiation from high energy electrons, or pion decay due to proton-proton interaction… ▽ More According to the most popular model for the origin of cosmic rays (CRs), supernova remnants (SNRs) are the site where CRs are accelerated. Observations across the electromagnetic spectrum support this picture through the detection of non-thermal emission that is compatible with being synchrotron or inverse Compton radiation from high energy electrons, or pion decay due to proton-proton interactions. These observations of growing quantity and quality promise to unveil many aspects of CRs acceleration and require more and more accurate tools for their interpretation. Here, we show how multi-dimensional MHD models of SNRs, including the effects on shock dynamics due to back-reaction of accelerated CRs and the synthesis of non-thermal emission, turned out to be very useful to investigate the signatures of CRs acceleration and to put constraints on the acceleration mechanism of high energy particles. These models have been used to interpret accurately observations of SNRs in various bands (radio, X-ray and $γ$-ray) and to extract from them key information about CRs acceleration. △ Less

Submitted 24 December, 2020; originally announced December 2020.

Comments: 16 pages, 5 figures, invited review accepted for publication on New Astronomy

Modeling the mixed-morphology supernova remnant IC 443. Origin of its complex morphology and X-ray emission

Authors: S. Ustamujic, S. Orlando, E. Greco, M. Miceli, F. Bocchino, A. Tutone, G. Peres

Abstract: The morphology and the distribution of material observed in SNRs reflect the interaction of the SN blast wave with the ambient environment, the physical processes associated with the SN explosion and the internal structure of the progenitor star. IC 443 is a MM SNR located in a quite complex environment: it interacts with a molecular cloud in the NW and SE areas and with an atomic cloud in the NE.… ▽ More The morphology and the distribution of material observed in SNRs reflect the interaction of the SN blast wave with the ambient environment, the physical processes associated with the SN explosion and the internal structure of the progenitor star. IC 443 is a MM SNR located in a quite complex environment: it interacts with a molecular cloud in the NW and SE areas and with an atomic cloud in the NE. In this work we aim at investigating the origin of the complex morphology and multi-thermal X-ray emission observed in SNR IC 443, through the study of the effect of the inhomogeneous ambient medium in shaping its observed structure, and the exploration of the main parameters characterizing the remnant. We developed a 3D HD model for IC 443, which describes the interaction of the SNR with the environment, parametrized in agreement with the results of the multi-wavelength data analysis. We performed an ample exploration of the parameter space describing the initial blast wave and the environment, and the surrounding clouds. From the simulations, we synthesized the X-ray emission maps and spectra and compared them with actual X-ray data collected by XMM-Newton. Our model explains the origin of the complex X-ray morphology of SNR IC 443 in a natural way, being able to reproduce, for the first time, most of the observed features, including the centrally-peaked X-ray morphology (characteristic of MM SNRs) when considering the origin of the explosion at the position where the PWN CXOU J061705.3+222127 was at the time of the explosion. In the model which best reproduces the observations, the mass of the ejecta and the energy of the explosion are $\sim 7 M_\odot$ and $\sim10^{51}$ erg, respectively. From the exploration of the parameter space, we found that the density of the clouds is $n>300$ cm$^{-3}$ and that the age of SNR IC 443 is $\sim8000$ yr. △ Less

Submitted 4 March, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 16 pages, 12 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 649, A14 (2021)

Correcting the effect of stellar spots on ARIEL transmission spectra

Authors: Gianluca Cracchiolo, Giuseppina Micela, Giovanni Peres

Abstract: The goal of this study is to assess the impact of the stellar spots on the extraction of the planetary transmission spectra observed by ARIEL. We develop a method to model the stellar spectrum of a star in the presence of spots by using the out-of-transit observations. It is based on a chi squared minimization procedure of the out-of-transit spectrum on a grid of stellar spectra with different siz… ▽ More The goal of this study is to assess the impact of the stellar spots on the extraction of the planetary transmission spectra observed by ARIEL. We develop a method to model the stellar spectrum of a star in the presence of spots by using the out-of-transit observations. It is based on a chi squared minimization procedure of the out-of-transit spectrum on a grid of stellar spectra with different sizes and temperatures of the spots. The approach allows us also to study the temporal evolution of the spots when comparing stellar spectra observed at different epochs. We also present a method to correct the transit depth variations due to non-occulted stellar spots and estimate the error we introduce if we apply the same correction to crossings over the stellar spots. The method is tested on three types of stellar targets that ARIEL will observe in its 4-years mission lifetime. In all the explored cases, the approach allows to reliably recover the spot parameters (size and temperature) from out-of-transit observations and, for non-occulted spots, to confidently recover the planetary atmosphere transmission spectrum within the noise level(with average uncertainty of at most 3.3% of the planetary signal). Conversely, we find systematic biases in the inferred planetary spectra due to the occulted spots, with measurable effects for the brightest targets especially for more contrasted spots. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Comments: 16 pages, 13 figures

Journal ref: Monthly Notices of the Royal Astronomical Society (2020)

The fully developed remnant of a neutrino-driven supernova: Evolution of ejecta structure and asymmetries in SNR Cassiopeia A

Authors: S. Orlando, A. Wongwathanarat, H. -T. Janka, M. Miceli, M. Ono, S. Nagataki, F. Bocchino, G. Peres

Abstract: Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant f… ▽ More Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant formation. The simulations cover 2000 years of expansion and include all physical processes relevant to describe the complexities in the SN evolution and the subsequent interaction of the stellar debris with the wind of the progenitor star. The interaction of large-scale asymmetries left from the earliest phases of the explosion with the reverse shock produces, at the age of $\approx 350$~years, an ejecta structure and a remnant morphology which are remarkably similar to those observed in Cas A. Small-scale structures in the large-scale Fe-rich plumes created during the initial stages of the SN, combined with HD instabilities that develop after the passage of the reverse shock, naturally produce a pattern of ring- and crown-like structures of shocked ejecta. The consequence is a spatial inversion of the ejecta layers with Si-rich ejecta being physically interior to Fe-rich ejecta. The full-fledged remnant shows voids and cavities in the innermost unshocked ejecta resulting from the expansion of Fe-rich plumes and their inflation due to the decay of radioactive species. The asymmetric distributions of $^{44}$Ti and $^{56}$Fe and their abundance ratio are both compatible with those inferred from high-energy observations of Chandra and NuSTAR. The main asymmetries observed in the ejecta distribution of Cas A can be explained by the interaction of the reverse shock with the large-scale asymmetries that developed from stochastic processes that originate during the first seconds of the SN blast. △ Less

Submitted 1 December, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 33 pages, 29 figures, 1 Table, 4 movies, 3 interactive graphics, accepted for publication on A&A, high quality figures will be available on the published paper

Journal ref: A&A 645, A66 (2021)

3D modeling from the onset of the SN to the full-fledged SNR: Role of an initial ejecta anisotropy on matter mixing

Authors: Antonio Tutone, Salvatore Orlando, Marco Miceli, Sabina Ustamujic, Masaomi Ono, Shigehiro Nagataki, Gilles Ferrand, Emanuele Greco, Giovanni Peres, Donald C. Warren, Fabrizio Bocchino

Abstract: The aim of this work is to bridge the gap between CC SNe and their remnants by investigating how post-explosion anisotropies in the ejecta influence the structure and chemical properties of the remnant at later times. We performed three-dimensional magneto-hydrodynamical simulations starting soon after the SN event and following the evolution of the system in the circumstellar medium (consisting o… ▽ More The aim of this work is to bridge the gap between CC SNe and their remnants by investigating how post-explosion anisotropies in the ejecta influence the structure and chemical properties of the remnant at later times. We performed three-dimensional magneto-hydrodynamical simulations starting soon after the SN event and following the evolution of the system in the circumstellar medium (consisting of the wind of the stellar progenitor), for 5000 years, obtaining the physical scenario of a SNR. Here we focused the analysis on the case of a progenitor red supergiant of 19.8 M_sun. We also investigated how a post-explosion large-scale anisotropy in the SN affects the ejecta distribution and the matter mixing of heavy elements in the remnant, during the first 5000 years of evolution. In the case of a spherically symmetric SN explosion without large-scale anisotropies, the remnant roughly keeps memory of the original onion-like layering of ejecta soon after the SN event. Nevertheless, as the reverse shock hits the ejecta, the element distribution departs from a homologous expansion, because of the slowing down of the outermost ejecta layers due to interaction with the reverse shock. In the case of a large-scale anisotropy developed after the SN, we found that the chemical stratification in the ejecta can be strongly modified and the original onion-like layering is not preserved. The anisotropy may cause spatial inversion of ejecta layers, for instance leading to Fe/Si-rich ejecta outside the O shell, and may determine the formation of Fe/Si-rich jet-like features that may protrude the remnant outline. The level of matter mixing and the properties of the jet-like feature are sensitive to the initial physical (density and velocity) and geometrical (size and position) initial characteristics of the anisotropy. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Comments: 16 pages, 14 figures, accepted to A&A

Report number: RIKEN-iTHEMS-Report-20

Journal ref: A&A 642, A67 (2020)

Supernova Neutrino Burst Detection with the Deep Underground Neutrino Experiment

Authors: DUNE collaboration, B. Abi, R. Acciarri, M. A. Acero, G. Adamov, D. Adams, M. Adinolfi, Z. Ahmad, J. Ahmed, T. Alion, S. Alonso Monsalve, C. Alt, J. Anderson, C. Andreopoulos, M. P. Andrews, F. Andrianala, S. Andringa, A. Ankowski, M. Antonova, S. Antusch, A. Aranda-Fernandez, A. Ariga, L. O. Arnold, M. A. Arroyave, J. Asaadi , et al. (949 additional authors not shown)

Abstract: The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The gen… ▽ More The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The general capabilities of DUNE for neutrino detection in the relevant few- to few-tens-of-MeV neutrino energy range will be described. As an example, DUNE's ability to constrain the $ν_e$ spectral parameters of the neutrino burst will be considered. △ Less

Submitted 29 May, 2021; v1 submitted 15 August, 2020; originally announced August 2020.

Comments: 29 pages, 17 figures; paper based on DUNE Technical Design Report. arXiv admin note: substantial text overlap with arXiv:2002.03005

Report number: FERMILAB-PUB-20-380-LBNF

Unveiling pure-metal ejecta X-ray emission in supernova remnants through their radiative recombination continuum

Authors: Emanuele Greco, Jacco Vink, Marco Miceli, Salvatore Orlando, Vladimir Domcek, Ping Zhou, Fabrizio Bocchino, Giovanni Peres

Abstract: Spectral analysis of X-ray emission from ejecta in supernova remnants (SNRs) is hampered by the low spectral resolution of CCD cameras, which creates a degeneracy between the best-fit values of abundances and emission measure. The combined contribution of shocked ambient medium and ejecta to the X-ray emission complicates the determination of the ejecta mass and chemical composition, leading to bi… ▽ More Spectral analysis of X-ray emission from ejecta in supernova remnants (SNRs) is hampered by the low spectral resolution of CCD cameras, which creates a degeneracy between the best-fit values of abundances and emission measure. The combined contribution of shocked ambient medium and ejecta to the X-ray emission complicates the determination of the ejecta mass and chemical composition, leading to big uncertainties in mass estimates and it can introduce a bias in the comparison between the observed ejecta composition and the yields predicted by explosive nucleosynthesis. We explore the capabilities of present and future spectral instruments with the aim of identifying a spectral feature which may allow us to discriminate between metal-rich and pure-metal plasmas in X-ray spectra of SNRs. We studied the behavior of the most common X-ray emission processes of an optically thin plasma in the high-abundance regime. We investigated spectral features of bremsstrahlung, radiative recombination continua (RRC) and line emission, by exploring a wide range of chemical abundances, temperatures and ionization parameters. We synthesized X-ray spectra from a 3D hydrodynamic (HD) simulation of Cas A, by using the response matrix from the Chandra/ACIS-S CCD detector and that of the XRISM/Resolve X-ray calorimeter. We found that a bright RRC shows up when the plasma is made of pure-metal ejecta, and a high spectral resolution is needed to identify this ejecta signature. We verified the applicability of our novel diagnostic tool and we propose a promising target for the future detection of such spectral feature: the southeastern Fe-rich clump of Cas A. While there is no way to unambiguously reveal pure-metal ejecta emission with CCD detectors, X-ray calorimeters will be able to pinpoint the presence of pure-metal RRC and to recover correctly absolute mass and the chemical composition of the ejecta. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Comments: 13 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 638, A101 (2020)

Hydrodynamic simulations unravel the progenitor-supernova-remnant connection in SN 1987A

Authors: S. Orlando, M. Ono, S. Nagataki, M. Miceli, H. Umeda, G. Ferrand, F. Bocchino, O. Petruk, G. Peres, K. Takahashi, T. Yoshida

Abstract: (Abridged) We aim at linking the dynamical and radiative properties of the remnant of SN 1987A to the geometrical and physical characteristics of the parent aspherical SN explosion and to the internal structure of its progenitor star. We performed 3D hydrodynamic simulations which describe the long-term evolution of SN 1987A from the onset of the SN to the full-fledged remnant at the age of 50 yea… ▽ More (Abridged) We aim at linking the dynamical and radiative properties of the remnant of SN 1987A to the geometrical and physical characteristics of the parent aspherical SN explosion and to the internal structure of its progenitor star. We performed 3D hydrodynamic simulations which describe the long-term evolution of SN 1987A from the onset of the SN to the full-fledged remnant at the age of 50 years, accounting for the pre-SN structure of the progenitor star. The simulations include all physical processes relevant for the complex phases of SN evolution and for the interaction of the SNR with the highly inhomogeneous ambient environment around SN 1987A. From the simulations, we synthesize observables to be compared with observations. By comparing the model results with observations, we constrained the initial SN anisotropy causing Doppler shifts observed in emission lines of heavy elements from ejecta, and leading to the remnant evolution observed in the X-ray band in the last 30 years. In particular, we found that the high mixing of ejecta unveiled by high redshifts and broadenings of [FeII] and $^{44}$Ti lines require a highly asymmetric SN explosion channeling a significant fraction of energy along an axis almost lying in the plane of the central equatorial ring around SN 1987A, roughly along the line-of-sight but with an offset of 40 deg, with the lobe propagating away from the observer slightly more energetic than the other. We found unambiguously that the observed distribution of ejecta and the dynamical and radiative properties of the SNR can be best reproduced if the structure of the progenitor star was that of a blue supergiant resulted from the merging of two massive stars. △ Less

Submitted 6 December, 2019; originally announced December 2019.

Comments: Accepted for publication in Astronomy & Astrophysics, 20 pages, 13 figures, 5 tables

Journal ref: A&A 636, A22 (2020)

Effects of Radiation in Accretion Regions of Classical T Tauri Stars: Pre-heating of accretion column in non-LTE regime

Authors: Salvatore Colombo, Laurent Ibgui, Salvatore Orlando, Rafael Rodriguez, Guadalupe Espinosa, Matthias González, Chantal Stehlé, Lionel de Sá, Costanza Argiroffi, Rosaria Bonito, Giovanni Peres

Abstract: Models and observations indicate that the impact of matter accreting onto the surface of young stars produces regions at the base of accretion columns, in which optically thin and thick plasma components coexist. Thus an accurate description of these impacts requires to account for the effects of absorption and emission of radiation. We study the effects of radiation emerging from shock-heated pla… ▽ More Models and observations indicate that the impact of matter accreting onto the surface of young stars produces regions at the base of accretion columns, in which optically thin and thick plasma components coexist. Thus an accurate description of these impacts requires to account for the effects of absorption and emission of radiation. We study the effects of radiation emerging from shock-heated plasma in impact regions on the structure of the pre-shock downfalling material. We investigate if a significant absorption of radiation occurs and if it leads to a pre-shock heating of the accreting gas. We developed a radiation hydrodynamics model describing an accretion column impacting onto the surface of a Classical T Tauri Star. The model takes into account the stellar gravity, the thermal conduction, and the effects of both radiative losses and absorption of radiation by matter in the non local thermodynamic equilibrium regime. After the impact, a hot slab of post-shock plasma develops at the base of the accretion column. Part of radiation emerging from the slab is absorbed by the pre-shock accreting material. As a result, the pre-shock accretion column gradually heats up to temperatures of $10^5$ K, forming a radiative precursor of the shock. The precursor has a thermal structure with the hottest part at $T \approx 10^5$ K, with size comparable to that of the hot slab, above the post-shock region. At larger distances the temperature gradually decreases to $T \approx 10^4$ K. △ Less

Submitted 20 August, 2019; v1 submitted 19 August, 2019; originally announced August 2019.

Comments: Letter Accepted for publication in A&A

Journal ref: A&A 629, L9 (2019)

Non-LTE radiation hydrodynamics in PLUTO

Authors: Salvatore Colombo, Laurent Ibgui, Salvatore Orlando, Raphael Rodriguez, Guadalupe Espinosa, Matthias González, Chantal Stehlé, Giovanni Peres

Abstract: Modeling the dynamics of most astrophysical structures requires an adequate description of the radiation-matter interaction. Several numerical (magneto)hydrodynamics codes were upgraded with a radiation module to fulfill this request. However, those among them that use either the flux-limited diffusion (FLD) or the M1 radiation moment approaches are restricted to the local thermodynamic equilibriu… ▽ More Modeling the dynamics of most astrophysical structures requires an adequate description of the radiation-matter interaction. Several numerical (magneto)hydrodynamics codes were upgraded with a radiation module to fulfill this request. However, those among them that use either the flux-limited diffusion (FLD) or the M1 radiation moment approaches are restricted to the local thermodynamic equilibrium (LTE). This assumption may be not valid in some astrophysical cases. We present an upgraded version of the LTE radiation-hydrodynamics module implemented in the PLUTO code, originally developed by Kolb et al. (2013), which we have extended to handle non-LTE regimes. Starting from the general frequency-integrated comoving-frame equations of radiation hydrodynamics (RHD), we have justified all the assumptions made to obtain the non-LTE equations actually implemented in the module, under the FLD approximation. An operator-split method is employed, with two substeps: the hydrodynamic part is solved with an explicit method by the solvers already available in PLUTO, the non-LTE radiation diffusion and energy exchange part is solved with an implicit method. The module is implemented in the PLUTO environment. It uses databases of radiative quantities that can be provided independently by the user: the radiative power loss, the Planck and Rosseland mean opacities. Our implementation has been validated through different tests, in particular radiative shock tests. The agreement with the semi-analytical solutions (when available) is good, with a maximum error of 7%. Moreover, we have proved that non-LTE approach is of paramount importance to properly model accretion shock structures. Our radiation FLD module represents a step toward the general non-LTE RHD modeling. The module is available, under request, for the community. △ Less

Submitted 10 July, 2019; originally announced July 2019.

Comments: 15 pages, 4 Figures, Accepted for publication in Section 15: Numerical methods and codes of Astronomy and Astrophysics. The official date of acceptance is 05/07/2019

A stellar flare-coronal mass ejection event revealed by X-ray plasma motions

Authors: C. Argiroffi, F. Reale, J. J. Drake, A. Ciaravella, P. Testa, R. Bonito, M. Miceli, S. Orlando, G. Peres

Abstract: Coronal mass ejections (CMEs), often associated with flares, are the most powerful magnetic phenomena occurring on the Sun. Stars show magnetic activity levels up to 10^4 times higher, and CME effects on stellar physics and circumstellar environments are predicted to be significant. However, stellar CMEs remain observationally unexplored. Using time-resolved high-resolution X-ray spectroscopy of a… ▽ More Coronal mass ejections (CMEs), often associated with flares, are the most powerful magnetic phenomena occurring on the Sun. Stars show magnetic activity levels up to 10^4 times higher, and CME effects on stellar physics and circumstellar environments are predicted to be significant. However, stellar CMEs remain observationally unexplored. Using time-resolved high-resolution X-ray spectroscopy of a stellar flare on the active star HR 9024 observed with Chandra/HETGS, we distinctly detected Doppler shifts in S XVI, Si XIV, and Mg XII lines that indicate upward and downward motions of hot plasmas (~10-25 MK) within the flaring loop, with velocity v~100-400 km/s, in agreement with a model of flaring magnetic tube. Most notably, we also detected a later blueshift in the O VIII line which reveals an upward motion, with v=90+/-30 km/s, of cool plasma (~4 MK), that we ascribe to a CME coupled to the flare. From this evidence we were able to derive a CME mass of 1x10^21 g and a CME kinetic energy of 5x10^34 erg. These values provide clues in the extrapolation of the solar case to higher activity levels, suggesting that CMEs could indeed be a major cause of mass and angular momentum loss. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Comments: Letter published on Nature Astronomy, see https://rdcu.be/bEAta

A new view of the corona of classical T Tauri stars: Effects of flaring activity in circumstellar disks

Authors: Salvatore Colombo, Salvatore Orlando, Giovanni Peres, Fabio Reale, Costanza Argiroffi, Rosaria Bonito, Laurent Ibgui, Chantal Stehlé

Abstract: Classical T Tauri stars (CTTSs) are young low-mass stellar objects accreting mass from their circumstellar disks. They are characterized by high levels of coronal activity as revealed by X-ray observations. This activity may affect the disk stability and the circumstellar environment. Here we investigate if an intense coronal activity due to flares occurring close to the accretion disk may perturb… ▽ More Classical T Tauri stars (CTTSs) are young low-mass stellar objects accreting mass from their circumstellar disks. They are characterized by high levels of coronal activity as revealed by X-ray observations. This activity may affect the disk stability and the circumstellar environment. Here we investigate if an intense coronal activity due to flares occurring close to the accretion disk may perturb the inner disk stability, disrupt the inner part of the disk and, possibly, trigger accretion phenomena with rates comparable with those observed. We model a magnetized protostar surrounded by an accretion disk through 3D magnetohydrodinamic simulations. We explore cases characterized by a dipole plus an octupole stellar magnetic field configuration and different density of the disk or by different levels of flaring activity. As a result of the simulated intense flaring activity, we observe the formation of several loops that link the star to the disk; all these loops build up a hot extended corona with an X-ray luminosity comparable with typical values observed in CTTSs. The intense flaring activity close to the disk can strongly perturb the disk stability. The flares trigger overpressure waves which travel through the disk and modify its configuration. Accretion funnels may be triggered by the flaring activity, thus contributing to the mass accretion rate of the star. Accretion rates synthesized from the simulations are in a range between 10^-10 and 10^-9M_sun yr^-1 The accretion columns can be perturbed by the flares and they can interact with each other, possibly merging together in larger streams. As a result, the accretion pattern can be rather complex: the streams are highly inhomogeneous, with a complex density structure, and clumped △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 12 pages, 9 figures

Journal ref: A&A 624, A50 (2019)

Collisionless shock heating of heavy ions in SN 1987A

Authors: Marco Miceli, Salvatore Orlando, David N. Burrows, Kari A. Frank, Costanza Argiroffi, Fabio Reale, Giovanni Peres, Oleh Petruk, Fabrizio Bocchino

Abstract: Astrophysical shocks at all scales, from those in the heliosphere up to the cosmological shock waves, are typically "collisionless", because the thickness of their jump region is much shorter than the collisional mean free path. Across these jumps, electrons, protons, and ions are expected to be heated at different temperatures. Supernova remnants (SNRs) are ideal targets to study collisionless pr… ▽ More Astrophysical shocks at all scales, from those in the heliosphere up to the cosmological shock waves, are typically "collisionless", because the thickness of their jump region is much shorter than the collisional mean free path. Across these jumps, electrons, protons, and ions are expected to be heated at different temperatures. Supernova remnants (SNRs) are ideal targets to study collisionless processes because of their bright post-shock emission and fast shocks. Although optical observations of Balmer-dominated shocks in young SNRs showed that the post-shock proton temperature is higher than the electron temperature, the actual dependence of the post-shock temperature on the particle mass is still widely debated. We tackle this longstanding issue through the analysis of deep multi-epoch and high-resolution observations of the youngest nearby supernova remnant, SN 1987A, made with the Chandra X-ray telescope. We introduce a novel data analysis method by studying the observed spectra in close comparison with a dedicated full 3-D hydrodynamic simulation. The simulation is able to reproduce self-consistently the whole broadening of the spectral lines of many ions altogether. We can therefore measure the post shock temperature of protons and selected ions through comparison of the model with observations. We have obtained information about the heating processes in collisional shocks by finding that the ion to proton temperature ratio is always significantly higher than one and increases linearly with the ion mass for a wide range of masses and shock parameters. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: The final version of this manuscript (including Methods) is published in Nature Astronomy and is publicly available at this link: https://rdcu.be/bhO7Z

3D MHD modeling of the expanding remnant of SN 1987A. Role of magnetic field and non-thermal radio emission

Authors: S. Orlando, M. Miceli, O. Petruk, M. Ono, S. Nagataki, M. A. Aloy, P. Mimica, S. -H. Lee, F. Bocchino, G. Peres, M. Guarrasi

Abstract: Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface… ▽ More Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface to the expansion of its remnant through the surrounding nebula by 3D MHD simulations. The model considers the radiative cooling, the deviations from equilibrium of ionization, the deviation from temperature-equilibration between electrons and ions, and a plausible configuration of the pre-SN ambient magnetic field. From the simulations, we synthesize both thermal X-ray and non-thermal radio emission and compare the model results with observations. Results. The presence of an ambient magnetic field reduces the erosion and fragmentation of the dense equatorial ring after the impact of the SN blast wave. As a result, the ring survives the passage of the blast, at least, during the time covered by the simulations (40 years). Our model is able to reproduce the morphology and lightcurves of SN 1987A in both X-ray and radio bands. The model reproduces the observed radio emission if the flux originating from the reverse shock is heavily suppressed. In this case, the radio emission originates mostly from the forward shock traveling through the H II region. Possible mechanisms for the suppression of emission from the reverse shock are investigated. We find that synchrotron self-absorption and free-free absorption have negligible effects on the emission during the interaction with the nebula. We suggest that the emission from the reverse shock at radio frequencies might be limited by highly magnetized ejecta. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: 15 pages, 9 Figures, accepted for publication on A&A, high quality figures will be available on the published paper

Journal ref: A&A 622, A73 (2019)

Discovery of a jet-like structure with overionized plasma in the SNR IC443

Authors: Emanuele Greco, Marco Miceli, Salvatore Orlando, Giovanni Peres, Eleonora Troja, Fabrizio Bocchino

Abstract: IC443 is a supernova remnant located in a quite complex environment since it interacts with nearby clouds. Indications for the presence of overionized plasma have been found though the possible physical causes of overionization are still debated. Moreover, because of its peculiar position and proper motion, it is not clear if the pulsar wind nebula (PWN) within the remnant is the relic of the IC44… ▽ More IC443 is a supernova remnant located in a quite complex environment since it interacts with nearby clouds. Indications for the presence of overionized plasma have been found though the possible physical causes of overionization are still debated. Moreover, because of its peculiar position and proper motion, it is not clear if the pulsar wind nebula (PWN) within the remnant is the relic of the IC443 progenitor star or just a rambling one seen in projection on the remnant. Here we address the study of IC443 plasma in order to clarify the relationship PWN-remnant, the presence of overionization and the origin of the latter. We analyzed two \emph{XMM-Newton} observations producing background-subtracted, vignetting-corrected and mosaicked images in two different energy bands and we performed a spatially resolved spectral analysis of the X-ray emission. We identified an elongated (jet-like) structure with Mg-rich plasma in overionization. The head of the jet is interacting with a molecular cloud and the jet is aligned with the position of the PWN at the instant of the supernova explosion. Interestingly, the direction of the jet of ejecta is somehow consistent with the direction of the PWN jet. Our discovery of a jet of ejecta in IC443 enlarge the sample of core-collapse SNRs with collimated ejecta structures. IC443's jet is the first one which shows overionized plasma, possibly associated with the adiabatic expansion of ejecta. The match between the jet's direction and the original position of the PWN strongly supports the association between the neutron star and IC443. △ Less

Submitted 17 April, 2018; originally announced April 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 615, A157 (2018)

Redshifted X-rays from the material accreting onto TW Hya: evidence of a low-latitude accretion spot

Authors: C. Argiroffi, J. J. Drake, R. Bonito, S. Orlando, G. Peres, M. Miceli

Abstract: High resolution spectroscopy, providing constraints on plasma motions and temperatures, is a powerful means to investigate the structure of accretion streams in CTTS. In particular, the accretion shock region, where the accreting material is heated to temperatures of a few MK as it continues its inward bulk motion, can be probed by X-ray spectroscopy. To attempt to detect for the first time the mo… ▽ More High resolution spectroscopy, providing constraints on plasma motions and temperatures, is a powerful means to investigate the structure of accretion streams in CTTS. In particular, the accretion shock region, where the accreting material is heated to temperatures of a few MK as it continues its inward bulk motion, can be probed by X-ray spectroscopy. To attempt to detect for the first time the motion of this X-ray-emitting post-shock material, we searched for a Doppler shift in the deep Chandra/HETGS observation of the CTTS TW Hya. This test should unveil the nature of this X-ray emitting plasma component in CTTS, and constrain the accretion stream geometry. We searched for a Doppler shift in the X-ray emission from TW Hya with two different methods, by measuring the position of a selected sample of emission lines, and by fitting the whole TW Hya X-ray spectrum, allowing the line-of-sight velocity to vary. We found that the plasma at T~2-4 MK has a line-of-sight velocity of 38.3+/-5.1 km/s with respect to the stellar photosphere. This result definitively confirms that this X-ray-emitting material originates in the post-shock region, at the base of the accretion stream, and not in coronal structures. The comparison of the observed velocity along the line of sight, 38.3+/-5.1 km/s, with the inferred intrinsic velocity of the post shock of TW Hya, v_post~110-120 km/s, indicates that the footpoints of the accretion streams on TW Hya are located at low latitudes on the stellar surface. Our results indicate that complex magnetic field geometries, such as that of TW Hya, permit low-latitude accretion spots. Moreover, since on TW Hya the redshift of the soft X-ray emission is very similar to that of the narrow component of the CIV resonance doublet at 1550 Ang, as found by Ardila et al. (2013), then the plasma at 2-4 MK and that at 0.1 MK likely originate in the same post-shock regions. △ Less

Submitted 4 September, 2017; v1 submitted 10 August, 2017; originally announced August 2017.

Comments: Accepted for publication in Astronomy & Astrophysics; 2nd version after language editor corrections; 16 pages, 8 figures, 6 tables

Accretion disk coronae of Intermediate Polar Cataclysmic Variables - 3D MagnetoHydro-Dynamic modeling and thermal X-ray emission

Authors: Enrico Barbera, Salvatore Orlando, Giovanni Peres

Abstract: IPCVs contain a magnetic, rotating white dwarf surrounded by a magnetically truncated accretion disk. To explain their strong flickering X-ray emission, accretion has been successfully taken into account. Nevertheless, observations suggest that accretion phenomena could not be the only process behind it. An intense flaring activity occurring on the surface of the disk may generate a corona, contri… ▽ More IPCVs contain a magnetic, rotating white dwarf surrounded by a magnetically truncated accretion disk. To explain their strong flickering X-ray emission, accretion has been successfully taken into account. Nevertheless, observations suggest that accretion phenomena could not be the only process behind it. An intense flaring activity occurring on the surface of the disk may generate a corona, contribute to the thermal X-ray emission and influence the system stability. Our purposes are: investigating the formation of an extended corona above the accretion disk, due to an intense flaring activity occurring on the disk surface; studying its effects on the disk and stellar magnetosphere; assessing its contribution to the observed X-ray flux. We have developed a 3D MHD model of a IPCV. The model takes into account gravity, disk viscosity, thermal conduction, radiative losses and coronal flare heating. To perform a parameter space exploration, several system conditions have been considered, with different magnetic field intensity and disk density values. From the results of the evolution of the model, we have synthesized the thermal X-ray emission. The simulations show the formation of an extended corona, linking disk and star. The flaring activity is capable of strongly influencing the disk configuration and its stability, effectively deforming the magnetic field lines. Hot plasma evaporation phenomena occur in the layer immediately above the disk. The flaring activity gives rise to a thermal X-ray emission in both the [0.1-2.0] keV and the [2.0-10] keV bands. An intense coronal activity occurring on the disk surface of an IPCV can affect the structure of the disk depending noticeably on the density of the disk and the magnetic field of the central object. Moreover, the synthesis of the thermal X-ray fluxes shows that this flaring activity may contribute to the observed thermal X-ray emission. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Journal ref: A&A 600, A105 (2017)

Hydrodynamic Modeling of Accretion Impacts in Classical T Tauri Stars: Radiative Heating of the Pre-shock Plasma

Authors: G. Costa, S. Orlando, G. Peres, C. Argiroffi, R. Bonito

Abstract: Context. It is generally accepted that, in Classical T Tauri Stars, the plasma from the circumstellar disc accretes onto the stellar surface with free fall velocity, and the impact generates a shock. The impact region is expected to contribute to emission in different spectral bands; many studies have confirmed that the X-rays arise from the post-shock plasma but, otherwise, there are no studies i… ▽ More Context. It is generally accepted that, in Classical T Tauri Stars, the plasma from the circumstellar disc accretes onto the stellar surface with free fall velocity, and the impact generates a shock. The impact region is expected to contribute to emission in different spectral bands; many studies have confirmed that the X-rays arise from the post-shock plasma but, otherwise, there are no studies in the literature investigating the origin of the observed UV emission which is apparently correlated to accretion. Aims. We investigated the effect of radiative heating of the infalling material by the post-shock plasma at the base of the accretion stream with the aim to identify in which region a significant part of the UV emission originates. Methods. We developed a 1D hydrodynamic model describing the impact of an accretion stream onto the stellar surface; the model takes into account the gravity, the radiative cooling of an optically thin plasma, the thermal conduction, and the heating due to absorption of X-ray radiation. The latter term represents the heating of the infalling plasma due to the absorption of X-rays emitted from the post-shock region. Results. We found that the radiative heating of the pre-shock plasma plays a non-negligible role in the accretion phenomenon. In particular, the dense and cold plasma of the pre-shock accretion column is gradually heated up to few $10^5$\;K due to irradiation of X-rays arising from the shocked plasma at the impact region. This heating mechanism does not affect significantly the dynamics of the post-shock plasma. On the other hand, a region of radiatively heated gas (that we consider a precursor) forms in the unshocked accretion column and contributes significantly to UV emission. Our model naturally reproduces the luminosity of UV emission lines correlated to accretion and shows that most of the UV emission originates from the precursor. △ Less

Submitted 5 September, 2016; originally announced September 2016.

Journal ref: A&A 597, A1 (2017)

X-Raying the Dark Side of Venus - Scatter from Venus Magnetotail?

Authors: M. Afshari, G. Peres, P. R. Jibben, A. Petralia, F. Reale, M. Weber

Abstract: This work analyzes the X-ray, EUV and UV emission apparently coming from the Earth-facing (dark) side of Venus as observed with Hinode/XRT and SDO/AIA during a transit across the solar disk occurred in 2012. We have measured significant X-Ray, EUV and UV flux from Venus dark side. As a check we have also analyzed a Mercury transit across the solar disk, observed with Hinode/XRT in 2006. We have us… ▽ More This work analyzes the X-ray, EUV and UV emission apparently coming from the Earth-facing (dark) side of Venus as observed with Hinode/XRT and SDO/AIA during a transit across the solar disk occurred in 2012. We have measured significant X-Ray, EUV and UV flux from Venus dark side. As a check we have also analyzed a Mercury transit across the solar disk, observed with Hinode/XRT in 2006. We have used the latest version of the Hinode/XRT Point Spread Function (PSF) to deconvolve Venus and Mercury X-ray images, in order to remove possible instrumental scattering. Even after deconvolution, the flux from Venus shadow remains significant while in the case of Mercury it becomes negligible. Since stray-light contamination affects the XRT Ti-poly filter data from the Venus transit in 2012, we performed the same analysis with XRT Al-mesh filter data, which is not affected by the light leak. Even the Al-mesh filter data show residual flux. We have also found significant EUV (304 A, 193 A, 335 A) and UV (1700 A) flux in Venus shadow, as measured with SDO/AIA. The EUV emission from Venus dark side is reduced when appropriate deconvolution methods are applied; the emission remains significant, however. The light curves of the average flux of the shadow in the X-ray, EUV, and UV bands appear different as Venus crosses the solar disk, but in any of them the flux is, at any time, approximately proportional to the average flux in a ring surrounding Venus, and therefore proportional to the average flux of the solar regions around Venus obscuring disk line of sight. The proportionality factor depends on the band. This phenomenon has no clear origin; we suggest it may be due to scatter occurring in the very long magnetotail of Venus. △ Less

Submitted 22 July, 2016; originally announced July 2016.

Comments: This paper has been accepted in The Astrophysical Journal

3D MHD modeling of twisted coronal loops

Authors: F. Reale, S. Orlando, M. Guarrasi, A. Mignone, G. Peres, A. W. Hood, E. R. Priest

Abstract: We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tu… ▽ More We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tube, in the solar atmosphere extending from the high-beta chromosphere to the low-beta corona through the steep transition region. The magnetic field expands from the chromosphere to the corona. The maximum resolution is ~30 km. We obtain an overall evolution typical of loop models and realistic loop emission in the EUV and X-ray bands. The plasma confined in the flux tube is heated to active region temperatures (~3 MK) after ~2/3 hr. Upflows from the chromosphere up to ~100 km/s fill the core of the flux tube to densities above 10^9 cm^-3. More heating is released in the low corona than the high corona and is finely structured both in space and time. △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: 29 pages, 14 figures, accepted for publication, movies available upon request to the first author

Impacts of fragmented accretion streams onto Classical T Tauri Stars: UV and X-ray emission lines

Authors: Salvatore Colombo, Salvatore Orlando, Giovanni Peres, Costanza Argiroffi, Fabio Reale

Abstract: Context. The accretion process in Classical T Tauri Stars (CTTSs) can be studied through the analysis of some UV and X-ray emission lines which trace hot gas flows and act as diagnostics of the post-shock downfalling plasma. In the UV band, where higher spectral resolution is available, these lines are characterized by rather complex profiles whose origin is still not clear. Aims. We investigate… ▽ More Context. The accretion process in Classical T Tauri Stars (CTTSs) can be studied through the analysis of some UV and X-ray emission lines which trace hot gas flows and act as diagnostics of the post-shock downfalling plasma. In the UV band, where higher spectral resolution is available, these lines are characterized by rather complex profiles whose origin is still not clear. Aims. We investigate the origin of UV and X-ray emission at impact regions of density structured (fragmented) accretion streams.We study if and how the stream fragmentation and the resulting structure of the post-shock region determine the observed profiles of UV and X-ray emission lines. Methods. We model the impact of an accretion stream consisting of a series of dense blobs onto the chromosphere of a CTTS through 2D MHD simulations. We explore different levels of stream fragmentation and accretion rates. From the model results, we synthesize C IV (1550 Å) and OVIII (18.97 Å) line profiles. Results. The impacts of accreting blobs onto the stellar chromosphere produce reverse shocks propagating through the blobs and shocked upflows. These upflows, in turn, hit and shock the subsequent downfalling fragments. As a result, several plasma components differing for the downfalling velocity, density, and temperature are present altoghether. The profiles of C IV doublet are characterized by two main components: one narrow and redshifted to speed $\approx$ 50 km s$^{-1}$ and the other broader and consisting of subcomponents with redshift to speed in the range 200 $\approx$ 400 km s$^{-1}$. The profiles of OVIII lines appear more symmetric than C IV and are redshifted to speed $\approx$ 150 km s$^{-1}$. Conclusions. Our model predicts profiles of C IV line remarkably similar to those observed and explains their origin in a natural way as due to stream fragmentation. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Comments: 11 pages, 10 figures

Modeling the shock-cloud interaction in SN 1006: unveiling the origin of nonthermal X-ray and gamma-ray emission

Authors: Marco Miceli, Salvatore Orlando, Victor Pereira, Fabio Acero, Satoru Katsuda, Anne Decourchelle, Frank P. Winkler, Rosaria Bonito, Fabio Reale, Giovanni Peres, Jiangtao Li, Gloria Dubner

Abstract: The supernova remnant SN 1006 is a source of high-energy particles and its southwestern limb is interacting with a dense ambient cloud, thus being a promising region for gamma-ray hadronic emission. We aim at describing the physics and the nonthermal emission associated with the shock-cloud interaction to derive the physical parameters of the cloud (poorly constrained by the data analysis), to asc… ▽ More The supernova remnant SN 1006 is a source of high-energy particles and its southwestern limb is interacting with a dense ambient cloud, thus being a promising region for gamma-ray hadronic emission. We aim at describing the physics and the nonthermal emission associated with the shock-cloud interaction to derive the physical parameters of the cloud (poorly constrained by the data analysis), to ascertain the origin of the observed spatial variations in the spectral properties of the X-ray synchrotron emission, and to predict spectral and morphological features of the resulting gamma-ray emission. We performed 3-D magnetohydrodynamic simulations modeling the evolution of SN 1006 and its interaction with the ambient cloud, and explored different model setups. By applying the REMLIGHT code on the model results, we synthesized the synchrotron X-ray emission, and compared it with actual observations, to constrain the parameters of the model. We also synthesized the leptonic and hadronic gamma-ray emission from the models, deriving constraints on the energy content of the hadrons accelerated at the southwestern limb. We found that the impact of the SN 1006 shock front with a uniform cloud with density 0.5 cm^-3 can explain the observed morphology, the azimuthal variations of the cutoff frequency of the X-ray synchrotron emission, and the shock proper motion in the interaction region. Our results show that the current upper limit for the total hadronic energy in the southwestern limb is 2.5e49 erg. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 593, A26 (2016)

EUV flickering of solar coronal loops: a new diagnostic of coronal heating

Authors: E. Tajfirouze, F. Reale, G. Peres, P. Testa

Abstract: A previous work of ours found the best agreement between EUV light curves observed in an active region core (with evidence of super-hot plasma) and those predicted from a model with a random combination of many pulse-heated strands with a power-law energy distribution. We extend that work by including spatially resolved strand modeling and by studying the evolution of emission along the loops in t… ▽ More A previous work of ours found the best agreement between EUV light curves observed in an active region core (with evidence of super-hot plasma) and those predicted from a model with a random combination of many pulse-heated strands with a power-law energy distribution. We extend that work by including spatially resolved strand modeling and by studying the evolution of emission along the loops in the EUV 94 A and 335 A channels of the Atmospheric Imaging Assembly on-board the Solar Dynamics Observatory. Using the best parameters of the previous work as the input of the present one, we find that the amplitude of the random fluctuations driven by the random heat pulses increases from the bottom to the top of the loop in the 94 A channel and, viceversa, from the top to the bottom in the 335 A channel. This prediction is confirmed by the observation of a set of aligned neighbouring pixels along a bright arc of an active region core. Maps of pixel fluctuations may therefore provide easy diagnostics of nano-flaring regions. △ Less

Submitted 15 January, 2016; originally announced January 2016.

Hidden Interactions of Sterile Neutrinos As a Probe For New Physics

Authors: Zahra Tabrizi, O. L. G. Peres

Abstract: Recent results from neutrino experiments show evidence for light sterile neutrinos which do not have any Standard Model interactions. In this work we study the hidden interaction of sterile neutrinos with an "MeV scale" gauge boson (the $ν_s$HI model) with mass $M_X$ and leptonic coupling $g^\prime_l$. By performing an analysis on the $ν_s$HI model using the data of the MINOS neutrino experiment w… ▽ More Recent results from neutrino experiments show evidence for light sterile neutrinos which do not have any Standard Model interactions. In this work we study the hidden interaction of sterile neutrinos with an "MeV scale" gauge boson (the $ν_s$HI model) with mass $M_X$ and leptonic coupling $g^\prime_l$. By performing an analysis on the $ν_s$HI model using the data of the MINOS neutrino experiment we find that the values above $G_X/G_F=92.4$ are excluded by more than $2σ$ C.L., where $G_F$ is the Fermi constant and $G_X$ is the field strength of the $ν_s$HI model. Using this model we can also probe other new physics scenarios. We find that the region allowed by the $(g-2)_μ$ discrepancy is entirely ruled out for $M_X\lesssim 100$ MeV. Finally, the secret interaction of sterile neutrinos has been to solve a conflict between the sterile neutrinos and cosmology. It is shown here that such an interaction is excluded by MINOS for $g^\prime_s> 1.6\times10^{-2}$. This exclusion, however, does depend on the value of $g_l^{\prime}$. △ Less

Submitted 17 February, 2016; v1 submitted 23 July, 2015; originally announced July 2015.

Comments: 11 pages, 2 figures, Improved version with new figures and further details. Accepted in PRD

Journal ref: Phys. Rev. D 93, 053003 (2016)

Neutrino Decay and Solar Neutrino Seasonal Effect

Authors: R. Picoreti, M. M. Guzzo, P. C. de Holanda, O. L. G. Peres

Abstract: We consider the possibility of solar neutrino decay as a sub-leading effect on their propagation between production and detection. Using current oscillation data, we set a new lower bound to the $ν_2$ neutrino lifetime at $τ_2\, /\, m_2 \geq 7.2 \times 10^{-4}\,\,\hbox{s}\,.\,\hbox{eV}^{-1}$ at $99\%\,$C.L.. Also, we show how seasonal variations in the solar neutrino data can give interesting addi… ▽ More We consider the possibility of solar neutrino decay as a sub-leading effect on their propagation between production and detection. Using current oscillation data, we set a new lower bound to the $ν_2$ neutrino lifetime at $τ_2\, /\, m_2 \geq 7.2 \times 10^{-4}\,\,\hbox{s}\,.\,\hbox{eV}^{-1}$ at $99\%\,$C.L.. Also, we show how seasonal variations in the solar neutrino data can give interesting additional information about neutrino lifetime. △ Less

Submitted 26 June, 2015; originally announced June 2015.

UV variability and accretion dynamics in the young open cluster NGC 2264

Authors: Laura Venuti, Jerome Bouvier, Jonathan Irwin, John Stauffer, Lynne Hillenbrand, Luisa Rebull, Ann Marie Cody, Silvia Alencar, Giuseppina Micela, Ettore Flaccomio, Giovanni Peres

Abstract: We explore UV and optical variability signatures for several hundred members of NGC 2264 (3 Myr). We performed simultaneous u- and r-band monitoring over two full weeks with CFHT/MegaCam. About 750 young stars are probed; 40% of them are accreting. Statistically distinct variability properties are observed for accreting and non-accreting cluster members. The accretors exhibit a significantly highe… ▽ More We explore UV and optical variability signatures for several hundred members of NGC 2264 (3 Myr). We performed simultaneous u- and r-band monitoring over two full weeks with CFHT/MegaCam. About 750 young stars are probed; 40% of them are accreting. Statistically distinct variability properties are observed for accreting and non-accreting cluster members. The accretors exhibit a significantly higher level of variability than the non-accretors, especially in the UV. The amount of u-band variability correlates statistically with UV excess in disk-bearing objects, which suggests that accretion and star-disk interaction are the main sources of variability. Cool magnetic spots, several hundred degrees colder than the photosphere and covering from 5 to 30% of the stellar surface, appear to be the leading factor of variability for the non-accreting stars. In contrast, accretion spots, a few thousand degrees hotter than the photosphere and covering a few percent of the stellar surface, best reproduce the variability of accreting objects. The color behavior is also found to be different between accreting and non-accreting stars. Typical variability amplitudes for accreting members rapidly increase from r to u, which indicates a much stronger contrast at short wavelengths; a lower color dependence in the amplitudes is instead measured for diskless stars. We find that u-band variability on hour timescales is typically about 10% of the peak-to-peak variability on day timescales, while longer term (years) variability is consistent with amplitudes measured over weeks. We conclude that for both accreting and non-accreting stars, the mid-term rotational modulation by spots is the leading timescale for a variability of up to several years. In turn, this suggests that the accretion process is essentially stable over years, although it exhibits low-level shorter term variations in single accretion events. △ Less

Submitted 23 June, 2015; originally announced June 2015.

Comments: 14 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 581, A66 (2015)

Magnetohydrodynamic modeling of the accretion shocks in classical T Tauri stars: the role of local absorption on the X-ray emission

Authors: R. Bonito, S. Orlando, C. Argiroffi, M. Miceli, G. Peres, T. Matsakos, C. Stehle, L. Ibgui

Abstract: We investigate the properties of X-ray emission from accretion shocks in classical T Tauri stars (CTTSs), generated where the infalling material impacts the stellar surface. Both observations and models of the accretion process reveal several aspects that are still unclear: the observed X-ray luminosity in accretion shocks is below the predicted value, and the density versus temperature structure… ▽ More We investigate the properties of X-ray emission from accretion shocks in classical T Tauri stars (CTTSs), generated where the infalling material impacts the stellar surface. Both observations and models of the accretion process reveal several aspects that are still unclear: the observed X-ray luminosity in accretion shocks is below the predicted value, and the density versus temperature structure of the shocked plasma, with increasing densities at higher temperature, deduced from the observations, is at odds with that proposed in the current picture of accretion shocks. To address these open issues we investigate whether a correct treatment of the local absorption by the surrounding medium is crucial to explain the observations. To this end, we describe the impact of an accretion stream on a CTTS by considering a magnetohydrodynamic model. From the model results we synthesize the X-ray emission from the accretion shock by producing maps and spectra. We perform density and temperature diagnostics on the synthetic spectra, and we directly compare the results with the observations. Our model shows that the X-ray fluxes inferred from the emerging spectra are lower than expected because of the complex local absorption by the optically thick material of the chromosphere and of the unperturbed stream. Moreover, our model including the effects of local absorption explains in a natural way the apparently puzzling pattern of density versus temperature observed in the X-ray emission from accretion shocks. △ Less

Submitted 10 October, 2014; originally announced October 2014.

Comments: Accepted for publication in Astrophysical Journal Letters; 5 pages, 4 figures

Mapping accretion and its variability in the young open cluster NGC 2264: a study based on u-band photometry

Authors: Laura Venuti, Jérôme Bouvier, Ettore Flaccomio, Silvia H. P. Alencar, Jonathan Irwin, John R. Stauffer, Ann Marie Cody, Paula S. Teixeira, Alana P. Sousa, Giuseppina Micela, Jean-Charles Cuillandre, Giovanni Peres

Abstract: We aim at characterizing the accretion properties of several hundred members of the star-forming cluster NGC 2264 (3 Myr). We performed a deep u,g,r,i mapping and a simultaneous u+r monitoring of the region with CFHT/MegaCam in order to directly probe the accretion process from UV excess measurements. Photometric properties and stellar parameters are determined homogeneously for about 750 monitore… ▽ More We aim at characterizing the accretion properties of several hundred members of the star-forming cluster NGC 2264 (3 Myr). We performed a deep u,g,r,i mapping and a simultaneous u+r monitoring of the region with CFHT/MegaCam in order to directly probe the accretion process from UV excess measurements. Photometric properties and stellar parameters are determined homogeneously for about 750 monitored young objects, spanning the mass range 0.1-2 Mo. About 40% are classical (accreting) T Tauri stars, based on various diagnostics (H_alpha, UV and IR excesses). The remaining non-accreting members define the (photospheric+chromospheric) reference UV emission level over which flux excess is detected and measured. We revise the membership status of cluster members based on UV accretion signatures and report a new population of 50 CTTS candidates. A large range of UV excess is measured for the CTTS population, varying from a few 0.1 to 3 mag. We convert these values to accretion luminosities and obtain mass accretion rates ranging from 1e-10 to 1e-7 Mo/yr. Taking into account a mass-dependent detection threshold for weakly accreting objects, we find a >6sigma correlation between mass accretion rate and stellar mass. A power-law fit, properly accounting for upper limits, yields M_acc $\propto$ M^{1.4+/-0.3}. At any given stellar mass, we find a large spread of accretion rates, extending over about 2 orders of magnitude. The monitoring of the UV excess on a timescale of a couple of weeks indicates that its variability typically amounts to 0.5 dex, much smaller than the observed spread. We suggest that a non-negligible age spread across the cluster may effectively contribute to the observed spread in accretion rates at a given mass. In addition, different accretion mechanisms (like, e.g., short-lived accretion bursts vs. more stable funnel-flow accretion) may be associated to different M_acc regimes. △ Less

Submitted 2 August, 2014; originally announced August 2014.

Comments: 24 pages, 21 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 570, A82 (2014)

Constraining the Violation of Equivalence Principle with IceCube Atmospheric Neutrino Data

Authors: A. Esmaili, D. R. Gratieri, M. M. Guzzo, P. C. de Holanda, O. L. G. Peres, G. A. Valdiviesso

Abstract: The recent high-statistics high-energy atmospheric neutrino data collected by IceCube open a new window to probe new physics scenarios that are suppressed in lower energy neutrino experiments. In this paper we analyze the IceCube atmospheric neutrino data to constrain the Violation of Equivalence Principle (VEP) in the framework of three neutrinos with non-universal gravitational couplings. In thi… ▽ More The recent high-statistics high-energy atmospheric neutrino data collected by IceCube open a new window to probe new physics scenarios that are suppressed in lower energy neutrino experiments. In this paper we analyze the IceCube atmospheric neutrino data to constrain the Violation of Equivalence Principle (VEP) in the framework of three neutrinos with non-universal gravitational couplings. In this scenario the effect of VEP on neutrino oscillation probabilities can be parametrized by two parameters $Δγ_{21}\equiv γ_2-γ_1$ and $Δγ_{31}\equiv γ_3-γ_1$, where $γ_i$'s denote the coupling of neutrino mass eigenstates to gravitational field. By analyzing the latest muon-tracks data sets of IceCube-40 and IceCube-79, besides providing the 2D allowed regions in $(φΔγ_{21},φΔγ_{31})$ plane, we obtain the upper limits $|φΔγ_{21}| < 9.1\times 10^{-27}$ (at 90\% C.L.) which improves the previous limit by $\sim4$ orders of magnitude and $|φΔγ_{31}| \lesssim 6\times 10^{-27}$ (at 90\% C.L.) which improves the current limit by $\sim1$ order of magnitude. Also we discuss in detail and analytically the effect of VEP on neutrino oscillation probabilities. △ Less

Submitted 21 August, 2014; v1 submitted 14 April, 2014; originally announced April 2014.

Comments: 15 pages, 5 figures, v2: references added, matches the version published in PRD

Journal ref: Phys. Rev. D 89, 113003 (2014)

A new look at Spitzer primary transit observations of the exoplanet HD189733b

Authors: Giuseppe Morello, Ingo P. Waldmann, Giovanna Tinetti, Giovanni Peres, Giuseppina Micela, Ian D. Howarth

Abstract: Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$μ$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datas… ▽ More Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$μ$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datasets reported discrepant results, hence the necessity of the reanalyses. The method we used here is based on the Independent Component Analysis (ICA) statistical technique, which ensures a high degree of objectivity. The use of ICA to detrend single photometric observations in a self-consistent way is novel in the literature. The advantage of our reanalyses over previous work is that we do not have to make any assumptions on the structure of the unknown instrumental systematics. Such "admission of ignorance" may result in larger error bars than reported in the literature, up to a factor $1.6$. This is a worthwhile trade-off for much higher objectivity, necessary for trustworthy claims. Our main results are (1) improved and robust values of orbital and stellar parameters, (2) new measurements of the transit depths at 3.6$μ$m, (3) consistency between the parameters estimated from the two observations, (4) repeatability of the measurement within the photometric level of $\sim 2 \times 10^{-4}$ in the IR, (5) no evidence of stellar variability at the same photometric level within 1 year. △ Less

Submitted 12 March, 2014; originally announced March 2014.

Comments: 43 pages, 18 figures

Journal ref: Morello et al. 2014, ApJ, 786, 22

Impact of sterile neutrinos on the early time flux from a galactic supernova

Authors: Arman Esmaili, O. L. G. Peres, Pasquale Dario Serpico

Abstract: We study the impact of the existence of an eV-mass scale sterile neutrino---with parameters in the ballpark of what is required to fit the laboratory anomalies---on the early time profile of the electron neutrino and antineutrino fluxes associated to a core-collapse supernova (SN). In particular, we focus on the universal feature of neutronization burst expected in the first tens of ms of the sign… ▽ More We study the impact of the existence of an eV-mass scale sterile neutrino---with parameters in the ballpark of what is required to fit the laboratory anomalies---on the early time profile of the electron neutrino and antineutrino fluxes associated to a core-collapse supernova (SN). In particular, we focus on the universal feature of neutronization burst expected in the first tens of ms of the signal: provided that a detector with sufficient sensitivity is available, it is well-known that in the 3 neutrino framework the detection of the neutronization burst in neutrino channel would signal inverted mass hierarchy. This conclusion is dramatically altered in the presence of a sterile neutrino: we study here both analytically and numerically the region in parameter space where this characteristic signal disappears, mimicking normal hierarchy expectations. Conversely, the detection of a peak consistent with expectations for inverted mass hierarchy would exclude the existence of a sterile state over a much wider parameter space than what required by laboratory anomalies fits, or even probed by detectors coming on-line in the near future. Additionally, we show the peculiar alteration in the energy-time double differential flux, with a delayed peak appearing for kinematical reasons, which might offer a remarkable signature in case of favorable parameters and for a high statistics detection of a Galactic SN. We also comment on additional potentially interesting effects in the electron antineutrino channel, if more than one angle in the active-sterile sector is non-vanishing. As an ancillary result that we derived in the technical resolution of the equations, in an appendix we report the Cayley-Hamilton formalism for the evolution of a four neutrino system in matter, generalizing existing results in the literature. △ Less

Submitted 21 August, 2014; v1 submitted 6 February, 2014; originally announced February 2014.

Comments: 23 pages, 7 figures; v2: improved discussions, figures and references added, a typo in one of the equations and figures corrected, conclusions unchanged; matches the published version in PRD

Journal ref: Phys. Rev. D 90, 033013 (2014)

Mass Accretion Processes in Young Stellar Objects: Role of Intense Flaring Activity

Authors: S. Orlando, F. Reale, G. Peres, A. Mignone

Abstract: According to the magnetospheric accretion scenario, young low-mass stars are surrounded by circumstellar disks which they interact with through accretion of mass. The accretion builds up the star to its final mass and is also believed to power the mass outflows, which may in turn have a significant role in removing the excess angular momentum from the star-disk system. Although the process of mass… ▽ More According to the magnetospheric accretion scenario, young low-mass stars are surrounded by circumstellar disks which they interact with through accretion of mass. The accretion builds up the star to its final mass and is also believed to power the mass outflows, which may in turn have a significant role in removing the excess angular momentum from the star-disk system. Although the process of mass accretion is a critical aspect of star formation, some of its mechanisms are still to be fully understood. On the other hand, strong flaring activity is a common feature of young stellar objects (YSOs). In the Sun, such events give rise to perturbations of the interplanetary medium. Similar but more energetic phenomena occur in YSOs and may influence the circumstellar environment. In fact, a recent study has shown that an intense flaring activity close to the disk may strongly perturb the stability of circumstellar disks, thus inducing mass accretion episodes (Orlando et al. 2011). Here we review the main results obtained in the field and the future perspectives. △ Less

Submitted 2 January, 2014; originally announced January 2014.

Comments: 4 pages, 2 Figures; accepted for publication on Acta Polytechnica (Proceedings of the Frascati Workshop 2013)

XIPE: the X-ray Imaging Polarimetry Explorer

Authors: Paolo Soffitta, Xavier Barcons, Ronaldo Bellazzini, João Braga, Enrico Costa, George W. Fraser, Szymon Gburek, Juhani Huovelin, Giorgio Matt, Mark Pearce, Juri Poutanen, Victor Reglero, Andrea Santangelo, Rashid A. Sunyaev, Gianpiero Tagliaferri, Martin Weisskopf, Roberto Aloisio, Elena Amato, Primo Attiná, Magnus Axelsson, Luca Baldini, Stefano Basso, Stefano Bianchi, Pasquale Blasi, Johan Bregeon , et al. (74 additional authors not shown)

Abstract: X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and wh… ▽ More X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017 but not selected. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus and two additional GPDs filled with pressurized Ar-DME facing the sun. The Minimum Detectable Polarization is 14 % at 1 mCrab in 10E5 s (2-10 keV) and 0.6 % for an X10 class flare. The Half Energy Width, measured at PANTER X-ray test facility (MPE, Germany) with JET-X optics is 24 arcsec. XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil). △ Less

Submitted 26 September, 2013; originally announced September 2013.

Comments: 49 pages, 14 figures, 6 tables. Paper published in Experimental Astronomy http://link.springer.com/journal/10686

Journal ref: Experimental Astronomy (Aug. 2013)

Neutrinos from collapsars

Authors: Florencia L. Vieyro, Gustavo E. Romero, Orlando L. G. Peres

Abstract: Long gamma-ray bursts (GRBs) are associated with the gravitational collapse of very massive stars. The central engine of a GRB can collimate relativistic jets that propagate inside the stellar envelope. The shock waves produced when the jet disrupts the stellar surface are capable of accelerating particles up to very high energies. If the jet has hadronic content, neutrinos will be produced via ch… ▽ More Long gamma-ray bursts (GRBs) are associated with the gravitational collapse of very massive stars. The central engine of a GRB can collimate relativistic jets that propagate inside the stellar envelope. The shock waves produced when the jet disrupts the stellar surface are capable of accelerating particles up to very high energies. If the jet has hadronic content, neutrinos will be produced via charged pion decays. The main goal of this work is to estimate the neutrino emission produced in the region close to the surface of the star, taking pion and muon cooling into account, along with subtle effects arising from neutrino production in a highly magnetized medium. We estimate the maximum energies of the different kinds of particles and solve the coupled transport equations for each species. Once the particle distributions are known, we calculate the intensity of neutrinos. We study the different effects on the neutrinos that can change the relative weight of different flavors. In particular, we consider the effects of neutrino oscillations, and of neutrino spin precession caused by strong magnetic fields. The expected neutrino signals from the shocks in the uncorking regions of Population III events is very weak, but the neutrino signal produced by Wolf-Rayet GRBs with $z<0.5$ is not far from the level of the atmospheric background. The IceCube experiment does not have the sensitivity to detect neutrinos from the implosion of the earliest stars, but a number of high-energy neutrinos may be detected from nearby long GRBs. The cumulative signal should be detectable over several years ($\sim 10$ yr) of integration with the full 86-string configuration. △ Less

Submitted 24 September, 2013; originally announced September 2013.

Comments: 11 pages, 7 figures, accepted for publication in A&A

Radiative accretion shocks along nonuniform stellar magnetic fields in classical T Tauri stars

Authors: S. Orlando, R. Bonito, C. Argiroffi, F. Reale, G. Peres, M. Miceli, T. Matsakos, C. Stehlé, L. Ibgui, L. de Sá, J. P. Chièze, T. Lanz

Abstract: (abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D a… ▽ More (abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D axisymmetric magnetohydrodynamic simulations. Our model takes into account the gravity, the radiative cooling, and the magnetic-field-oriented thermal conduction. RESULTS. The structure, stability, and location of the shocked plasma strongly depend on the configuration and strength of the magnetic field. For weak magnetic fields, a large component of B may develop perpendicular to the stream at the base of the accretion column, limiting the sinking of the shocked plasma into the chromosphere. An envelope of dense and cold chromospheric material may also develop around the shocked column. For strong magnetic fields, the field configuration determines the position of the shock and its stand-off height. If the field is strongly tapered close to the chromosphere, an oblique shock may form well above the stellar surface. In general, a nonuniform magnetic field makes the distribution of emission measure vs. temperature of the shocked plasma lower than in the case of uniform magnetic field. CONCLUSIONS. The initial strength and configuration of the magnetic field in the impact region of the stream are expected to influence the chromospheric absorption and, therefore, the observability of the shock-heated plasma in the X-ray band. The field strength and configuration influence also the energy balance of the shocked plasma, its emission measure at T > 1 MK being lower than expected for a uniform field. The above effects contribute in underestimating the mass accretion rates derived in the X-ray band. △ Less

Submitted 19 September, 2013; originally announced September 2013.

Comments: 11 pages, 11 Figures; accepted for publication on A&A. Version with full resolution images can be found at http://www.astropa.unipa.it/~orlando/PREPRINTS/sorlando_accretion_shocks.pdf

Probing light sterile neutrinos in medium baseline reactor experiments

Authors: Arman Esmaili, Ernesto Kemp, O. L. G. Peres, Zahra Tabrizi

Abstract: Medium baseline reactor experiments (Double Chooz, Daya Bay and RENO) provide a unique opportunity to test the presence of light sterile neutrinos. We analyze the data of these experiments in the search of sterile neutrinos and also test the robustness of theta_13 determination in the presence of sterile neutrinos. We show that existence of a light sterile neutrino state improves the fit to these… ▽ More Medium baseline reactor experiments (Double Chooz, Daya Bay and RENO) provide a unique opportunity to test the presence of light sterile neutrinos. We analyze the data of these experiments in the search of sterile neutrinos and also test the robustness of theta_13 determination in the presence of sterile neutrinos. We show that existence of a light sterile neutrino state improves the fit to these data moderately. We also show that the measured value of theta_{13} by these experiments is reliable even in the presence of sterile neutrinos, and the reliability owes significantly to the Daya Bay and RENO data. From the combined analysis of the data of these experiments we constrain the mixing of a sterile neutrino with mass squared difference in the range of (10^{-3}-10^{-1}) eV^2 to sin^2 2 theta_{14} <0.1 at 95 % C.L.. △ Less

Submitted 30 October, 2013; v1 submitted 28 August, 2013; originally announced August 2013.

Comments: v2: improved discussion, references added, matched the version published in PRD

Journal ref: Phys. Rev. D 88, 073012 (2013)