Showing 1–20 of 20 results for author: Tolstov, A

A Model for Fast Blue Optical Transient AT2018Cow: Circumstellar Interaction of a Pulsational Pair-instability Supernova

Authors: Shing-Chi Leung, Sergei Blinnikov, Ken'ichi Nomoto, Petr Baklanov, Elena Sorokina, Alexey Tolstov

Abstract: The Fast Blue Optical Transient (FBOT) ATLAS18qqn (AT2018cow) has a light curve as bright as superluminous supernovae but rises and falls much faster. We model this light curve by circumstellar interaction of a pulsational pair-instability (PPI) supernova (SN) model based on our PPISN models studied in previous work. We focus on the 42 $M_\odot$ He star (core of a 80 $M_{\odot}$ star) which has ci… ▽ More The Fast Blue Optical Transient (FBOT) ATLAS18qqn (AT2018cow) has a light curve as bright as superluminous supernovae but rises and falls much faster. We model this light curve by circumstellar interaction of a pulsational pair-instability (PPI) supernova (SN) model based on our PPISN models studied in previous work. We focus on the 42 $M_\odot$ He star (core of a 80 $M_{\odot}$ star) which has circumstellar matter of mass 0.50 $M_\odot$. With the parameterized mass cut and the kinetic energy of explosion $E$, we perform hydrodynamical calculations of nucleosynthesis and optical light curves of PPISN models. The optical light curve of the first $\sim$ 20 days of AT2018cow is well-reproduced by the shock heating of circumstellar matter for the $42 ~M_{\odot}$ He star with $E = 5 \times 10^{51}$ erg. After day 20, the light curve is reproduced by the radioactive decay of 0.6 $M_\odot$ $^{56}$Co, which is a decay product of $^{56}$Ni in the explosion. We also examine how the light curve shape depends on the various model parameters, such as CSM structure and composition. We also discuss (1) other possible energy sources and their constraints, (2) origin of observed high-energy radiation, and (3) how our result depends on the radiative transfer codes. Based on our successful model for AT2018cow and the model for SLSN with the CSM mass as large as $20 ~M_\odot)$, we propose the working hypothesis that PPISN produces SLSNe if CSM is massive enough and FBOTs if CSM is less than $\sim 1 ~M_\odot$. △ Less

Submitted 12 March, 2021; v1 submitted 26 August, 2020; originally announced August 2020.

Comments: 17 pages, 31 figures, published in Astrophysical Journal, submitted at 7 February 2020, accepted at 18 September 2020, published at 3 November 2020

Journal ref: Astrophysical Journal 903, 66 (2020)

Three-dimensional Models of Core-collapse Supernovae From Low-mass Progenitors With Implications for Crab

Authors: G. Stockinger, H. -Th. Janka, D. Kresse, T. Melson, T. Ertl, M. Gabler, A. Gessner, A. Wongwathanarat, A. Tolstov, S. -C. Leung, K. Nomoto, A. Heger

Abstract: We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse… ▽ More We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth. △ Less

Submitted 10 June, 2020; v1 submitted 5 May, 2020; originally announced May 2020.

Comments: 47 pages, 27 figures, 6 tables; minor revisions, accepted by MNRAS

A rapidly declining transient discovered with Subaru/Hyper Suprime-Cam

Authors: Nozomu Tominaga, Tomoki Morokuma, Masaomi Tanaka, Naoki Yasuda, Hisanori Furusawa, Masayuki Tanaka, Ji-an Jiang, Alexey Tolstov, Sergei Blinnikov, Mamoru Doi, Ikuru Iwata, Hanindyo Kuncarayakti, Takashi J. Moriya, Tohru Nagao, Ken'ichi Nomoto, Junichi Noumaru, Tadafumi Takata

Abstract: We perform a high-cadence transient survey with Subaru Hyper Suprime-Cam (HSC), which we call the Subaru HSC survey Optimized for Optical Transients (SHOOT). We conduct HSC imaging observations with time intervals of about one hour on two successive nights, and spectroscopic and photometric follow-up observations. A rapidly declining blue transient SHOOT14di at $z=0.4229$ is found in observations… ▽ More We perform a high-cadence transient survey with Subaru Hyper Suprime-Cam (HSC), which we call the Subaru HSC survey Optimized for Optical Transients (SHOOT). We conduct HSC imaging observations with time intervals of about one hour on two successive nights, and spectroscopic and photometric follow-up observations. A rapidly declining blue transient SHOOT14di at $z=0.4229$ is found in observations on two successive nights with an image subtraction technique. The rate of brightness change is $+1.28^{+0.40}_{-0.27}~{\rm mag~day^{-1}}$ ($+1.83^{+0.57}_{-0.39}~{\rm mag~day^{-1}}$) in the observer (rest) frame and the rest-frame color between $3400$ and $4400~\unicode[.8,0]{x212B}$ is $M_{\rm 3400\unicode[.8,0]{x212B}}-M_{\rm 4400\unicode[.8,0]{x212B}}=-0.4$. The nature of the object is investigated by comparing its peak luminosity, decline rate, and color with those of transients and variables previously observed, and those of theoretical models. None of the transients or variables share the same properties as SHOOT14di. Comparisons with theoretical models demonstrate that, while the emission from the cooling envelope of a Type IIb supernova shows a slower decline rate than SHOOT14di, and the explosion of a red supergiant star with a dense circumstellar wind shows a redder color than SHOOT14di, the shock breakout at the stellar surface of the explosion of a $25M_{\odot}$ red supergiant star with a small explosion energy of $\leq0.4\times10^{51}$ erg reproduces the multicolor light curve of SHOOT14di. This discovery shows that a high-cadence, multicolor optical transient survey at intervals of about one hour, and continuous and immediate follow-up observations, is important for studies of normal core-collapse supernovae at high redshifts. △ Less

Submitted 10 October, 2019; v1 submitted 7 September, 2019; originally announced September 2019.

Comments: 17 pages, 11 figures. Accepted for publication in the Astrophysical Journal

Neutrino self-interaction and MSW effects on the supernova neutrino-process

Authors: Heamin Ko, Myung-Ki Cheoun, Eunja Ha, Motohiko Kusakabe, Takehito Hayakawa, Hirokazu Sasaki, Toshitaka Kajino, Masa-aki Hashimoto, Masaomi Ono, Mark D. Usang, Satoshi Chiba, Ko Nakamura, Alexey Tolstov, Ken'ichi Nomoto, Toshihiko Kawano, Grant J. Mathews

Abstract: We calculate the abundances of $^{7}$Li, $^{11}$B, $^{92}$Nb, $^{98}$Tc, $^{138}$La, and $^{180}$Ta produced by neutrino $(ν)$ induced reactions in a core-collapse supernova explosion. We consider the modification by $ν$ self-interaction ($ν$-SI) near the neutrinosphere and the Mikheyev-Smirnov-Wolfenstein effect in outer layers for time-dependent neutrino energy spectra. Abundances of $^{7}$Li an… ▽ More We calculate the abundances of $^{7}$Li, $^{11}$B, $^{92}$Nb, $^{98}$Tc, $^{138}$La, and $^{180}$Ta produced by neutrino $(ν)$ induced reactions in a core-collapse supernova explosion. We consider the modification by $ν$ self-interaction ($ν$-SI) near the neutrinosphere and the Mikheyev-Smirnov-Wolfenstein effect in outer layers for time-dependent neutrino energy spectra. Abundances of $^{7}$Li and heavy isotopes $^{92}$Nb, $^{98}$Tc and $^{138}$La are reduced by a factor of 1.5-2.0 by the $ν$-SI. In contrast, $^{11}$B is relatively insensitive to the $ν$-SI. We find that the abundance ratio of heavy to light nucleus, $^{138}$La/$^{11}$B, is sensitive to the neutrino mass hierarchy, and the normal mass hierarchy is more likely to be consistent with the solar abundances. △ Less

Submitted 23 November, 2019; v1 submitted 28 February, 2019; originally announced March 2019.

Journal ref: Astrophys. J. Lett. 891, L24 (2023)

Catching Element Formation In The Act

Authors: Chris L. Fryer, Frank Timmes, Aimee L. Hungerford, Aaron Couture, Fred Adams, Wako Aoki, Almudena Arcones, David Arnett, Katie Auchettl, Melina Avila, Carles Badenes, Eddie Baron, Andreas Bauswein, John Beacom, Jeff Blackmon, Stephane Blondin, Peter Bloser, Steve Boggs, Alan Boss, Terri Brandt, Eduardo Bravo, Ed Brown, Peter Brown, Steve Bruenn. Carl Budtz-Jorgensen, Eric Burns , et al. (194 additional authors not shown)

Abstract: Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray… ▽ More Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 14 pages including 3 figures

Report number: LA-UR-18-29748

Observations of SN 2017ein Reveal Shock Breakout Emission and A Massive Progenitor Star for a Type Ic Supernova

Authors: Danfeng Xiang, Xiaofeng Wang, Jun Mo, Lingjun Wang, Stephen Smartt, Morgan Fraser, Shuhrat A. Ehgamberdiev, Davron Mirzaqulov, Jujia Zhang, Tianmeng Zhang, Jozsef Vinko, J. Craig Wheeler, Griffin Hosseinzadeh, D. Andrew Howell, Curtis McCully, James M DerKacy, E. Baron, Peter Brown, Xianfei Zhang, Shaolan Bi, Hao Song, Kaicheng Zhang, A. Rest, Ken'ichi Nomoto, Alexey Tolstov , et al. (1 additional authors not shown)

Abstract: We present optical and ultraviolet observations of nearby type Ic supernova SN 2017ein as well as detailed analysis of its progenitor properties from both the early-time observations and the prediscovery Hubble Space Telescope (HST) images. The optical light curves started from within one day to $\sim$275 days after explosion, and optical spectra range from $\sim$2 days to $\sim$90 days after expl… ▽ More We present optical and ultraviolet observations of nearby type Ic supernova SN 2017ein as well as detailed analysis of its progenitor properties from both the early-time observations and the prediscovery Hubble Space Telescope (HST) images. The optical light curves started from within one day to $\sim$275 days after explosion, and optical spectra range from $\sim$2 days to $\sim$90 days after explosion. Compared to other normal SNe Ic like SN 2007gr and SN 2013ge, \mbox{SN 2017ein} seems to have more prominent C{\footnotesize II} absorption and higher expansion velocities in early phases, suggestive of relatively lower ejecta mass. The earliest photometry obtained for \mbox{SN 2017ein} show indications of shock cooling. The best-fit obtained by including a shock cooling component gives an estimate of the envelope mass as $\sim$0.02 M$_{\odot}$ and stellar radius as 8$\pm$4 R$_{\odot}$. Examining the pre-explosion images taken with the HST WFPC2, we find that the SN position coincides with a luminous and blue point-like source, with an extinction-corrected absolute magnitude of M$_V$$\sim$$-$8.2 mag and M$_I$$\sim$$-$7.7 mag.Comparisons of the observations to the theoretical models indicate that the counterpart source was either a single WR star or a binary with whose members had high initial masses, or a young compact star cluster. To further distinguish between different scenarios requires revisiting the site of the progenitor with HST after the SN fades away. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: 28 pages, 19 figures; accepted for publication in The Astrophysical Journal

Light curve modeling of fast-evolving supernova KSN 2015K: explosion in circumstellar matter of Super-AGB progentor

Authors: Alexey Tolstov, Ken'ichi Nomoto, Elena Sorokina, Sergei Blinnikov, Nozomu Tominaga, Yoshiaki Taniguchi

Abstract: Recent supernova searches revealed a number of fast-evolving luminous transients. We perform radiation hydrodynamic simulations of light curves of several models of supernova from super--asymptotic giant branch (super-AGB, SAGB) stars with low mass envelopes ($M_{\rm env}$ $\sim$ $0.05-1.25$ $M_{\odot}$). The differences in the light curves found among the models are used to link the observed even… ▽ More Recent supernova searches revealed a number of fast-evolving luminous transients. We perform radiation hydrodynamic simulations of light curves of several models of supernova from super--asymptotic giant branch (super-AGB, SAGB) stars with low mass envelopes ($M_{\rm env}$ $\sim$ $0.05-1.25$ $M_{\odot}$). The differences in the light curves found among the models are used to link the observed events to the most appropriate models. In particular, we propose that KSN 2015K is an electron-capture supernova. We assume "optically" thick CSM around SAGB and the circumstellar interaction powers the peak luminosity of the light curve with a short rise time. The faint tail might be influenced by the spin-down luminosity of a newborn Crab-like pulsar. Our fits indicate an ejecta mass of $0.02-0.05$ $M_{\odot}$, a circumstellar medium (CSM) mass of $0.10-0.12$ $M_{\odot}$, a radius of the CSM photosphere $\sim$ 10$^{14}$ cm, a kinetic energy of $\sim$ 3 $\times$ 10$^{50}$ erg, a photospheric velocity $v_{\rm ph} \gtrsim 10,000$ km s$^{-1}$ and a pulsar total spin energy ($2.5-4$) $\times$ 10$^{49}$ erg. △ Less

Submitted 9 June, 2019; v1 submitted 22 November, 2018; originally announced November 2018.

Comments: 5 pages, 5 figures, Accepted for publication in the Astrophysical Journal Letters on 06 Jun 2019

Type Ib/Ic supernovae: effect of nickel mixing on the early-time color evolution and implications for the progenitors

Authors: Sung-Chul Yoon, Wonseok Chun, Alexey Tolstov, Sergey Blinnikov, Luc Dessart

Abstract: We investigate the effect of mixing of radioactive nickel ($^{56}$Ni) on the early-time color evolution of Type Ib and Ic supernovae (SNe Ib/Ic) using multi-group radiation hydrodynamics simulations. We consider both helium-rich and helium-poor progenitors. Mixing of $^{56}$Ni is parameterized using a Gaussian distribution function. We find that the early-time color evolution with a weak $^{56}$Ni… ▽ More We investigate the effect of mixing of radioactive nickel ($^{56}$Ni) on the early-time color evolution of Type Ib and Ic supernovae (SNe Ib/Ic) using multi-group radiation hydrodynamics simulations. We consider both helium-rich and helium-poor progenitors. Mixing of $^{56}$Ni is parameterized using a Gaussian distribution function. We find that the early-time color evolution with a weak $^{56}$Ni mixing is characterized by three different phases: initial rapid reddening, blueward evolution due to the delayed effect of $^{56}$Ni heating, and redward evolution thereafter until the transition to the nebular phase. With a strong $^{56}$Ni mixing, the second phase disappears. We compare our models with the early-time color evolution of several SNe Ib/Ic (SN1999ex, SN 2008D, SN 2009jf, iPTF13bvn, SN 1994I, SN 2007gr, SN 2013ge, and 2017ein) and find signatures of relatively weak and strong $^{56}$Ni mixing for SNe Ib and SNe Ic, respectively. This suggests that SNe Ib progenitors are distinct from SN Ic progenitors in terms of helium content and that $^{56}$Ni mixing is generally stronger in the carbon-oxygen core and weaker in the helium-rich envelope. We conclude that the early-time color evolution is a powerful probe of $^{56}$Ni mixing in SNe Ib/Ic. △ Less

Submitted 7 October, 2018; originally announced October 2018.

Comments: 20 pages, 12 figures, 3 tables, submitted

Ultraviolet light curves of Gaia16apd in superluminous supernova models

Authors: Alexey Tolstov, Andrey Zhiglo, Ken'ichi Nomoto, Elena Sorokina, Alexandra Kozyreva, Sergei Blinnikov

Abstract: Observations of Gaia16apd revealed extremely luminous ultraviolet emission among superluminous supernovae (SLSNe). Using radiation hydrodynamics simulations we perform a comparison of UV light curves, color temperatures and photospheric velocities between the most popular SLSN models: pair-instability supernova, magnetar and interaction with circumstellar medium. We find that the interaction model… ▽ More Observations of Gaia16apd revealed extremely luminous ultraviolet emission among superluminous supernovae (SLSNe). Using radiation hydrodynamics simulations we perform a comparison of UV light curves, color temperatures and photospheric velocities between the most popular SLSN models: pair-instability supernova, magnetar and interaction with circumstellar medium. We find that the interaction model is the most promising to explain the extreme UV luminosity of Gaia16apd. The differences in late-time UV emission and in color evolution found between the models can be used to link an observed SLSN event to the most appropriate model. Observations at UV wavelengths can be used to clarify the nature of SLSNe and more attention should be paid to them in future follow-up observations. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: 5 pages, 6 figures, accepted for publication in the Astrophysical Journal Letters on 18 Jul 2017

A New, Faint Population of X-ray Transients

Authors: F. E. Bauer, E. Treister, K. Schawinski, S. Schulze, B. Luo, D. M. Alexander, W. N. Brandt, A. Comastri, F. Forster, R. Gilli, D. A. Kann, K. Maeda, K. Nomoto, M. Paolillo, P. Ranalli, D. P. Schneider, O. Shemmer, M. Tanaka, A. Tolstov, N. Tominaga, P. Tozzi, C. Vignali, J. Wang, Y. Xue, G. Yang

Abstract: We report on the detection of a remarkable new fast high-energy transient found in the Chandra Deep Field-South, robustly associated with a faint ($m_{\rm R}=27.5$ mag, $z_{\rm ph}$$\sim$2.2) host in the CANDELS survey. The X-ray event is comprised of 115$^{+12}_{-11}$ net 0.3-7.0 keV counts, with a light curve characterised by a $\approx$100 s rise time, a peak 0.3-10 keV flux of $\approx$5… ▽ More We report on the detection of a remarkable new fast high-energy transient found in the Chandra Deep Field-South, robustly associated with a faint ($m_{\rm R}=27.5$ mag, $z_{\rm ph}$$\sim$2.2) host in the CANDELS survey. The X-ray event is comprised of 115$^{+12}_{-11}$ net 0.3-7.0 keV counts, with a light curve characterised by a $\approx$100 s rise time, a peak 0.3-10 keV flux of $\approx$5$\times$10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, and a power-law decay time slope of $-1.53\pm0.27$. The average spectral slope is $Γ=1.43^{+0.23}_{-0.13}$, with no clear spectral variations. The \hbox{X-ray} and multi-wavelength properties effectively rule out the vast majority of previously observed high-energy transients. A few theoretical possibilities remain: an "orphan" X-ray afterglow from an off-axis short-duration Gamma-ray Burst (GRB) with weak optical emission; a low-luminosity GRB at high redshift with no prompt emission below $\sim$20 keV rest-frame; or a highly beamed Tidal Disruption Event (TDE) involving an intermediate-mass black hole and a white dwarf with little variability. However, none of the above scenarios can completely explain all observed properties. Although large uncertainties exist, the implied rate of such events is comparable to those of orphan and low-luminosity GRBs as well as rare TDEs, implying the discovery of an untapped regime for a known transient class, or a new type of variable phenomena whose nature remains to be determined. △ Less

Submitted 14 February, 2017; originally announced February 2017.

Comments: 18 pages, 7 figures, 4 tables. Accepted for publication in MNRAS

Pulsational Pair-Instability Model for Superluminous Supernova PTF12dam: Interaction and Radioactive Decay

Authors: Alexey Tolstov, Ken'ichi Nomoto, Sergei Blinnikov, Elena Sorokina, Robert Quimby, Petr Baklanov

Abstract: Being a superluminous supernova (SLSN), PTF12dam can be explained by a $^{56}$Ni-powered model, a magnetar-powered model or an interaction model. We propose that PTF12dam is a pulsational pair instability supernova, where the outer envelope of a progenitor is ejected during the pulsations. Thus, it is powered by double energy source: radioactive decay of $^{56}$Ni and a radiative shock in a dense… ▽ More Being a superluminous supernova (SLSN), PTF12dam can be explained by a $^{56}$Ni-powered model, a magnetar-powered model or an interaction model. We propose that PTF12dam is a pulsational pair instability supernova, where the outer envelope of a progenitor is ejected during the pulsations. Thus, it is powered by double energy source: radioactive decay of $^{56}$Ni and a radiative shock in a dense circumstellar medium. To describe multicolor light curves and spectra we use radiation hydrodynamics calculations of STELLA code. We found that light curves are well described in the model with 40M$_{\odot}$ ejecta and 20-40M$_{\odot}$ circumstellar medium. The ejected $^{56}$Ni mass is about 6M$_{\odot}$ which results from explosive nucleosynthesis with large explosion energy (2-3)$\cdot$10$^{52}$ ergs. In comparison with alternative scenarios of pair-instability supernova and magnetar-powered supernova, in interaction model all the observed main photometric characteristics are well reproduced: multicolor light curves, color temperatures, and photospheric velocities. △ Less

Submitted 29 December, 2016; v1 submitted 5 December, 2016; originally announced December 2016.

Comments: 11 pages, 17 figures, accepted for publication in the Astrophysical Journal on 28 Dec 2016

Fast evolving pair-instability supernova models: evolution, explosion, light curves

Authors: Alexandra Kozyreva, Matthew Gilmer, Raphael Hirschi, Carla Frohlich, Sergey Blinnikov, Ryan T. Wollaeger, Ulrich M. Noebauer, Daniel R. van Rossum, Alexander Heger, Wesley P. Even, Roni Waldman, Alexey Tolstov, Emmanouil Chatzopoulos, Elena Sorokina

Abstract: With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before t… ▽ More With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition. △ Less

Submitted 4 October, 2016; originally announced October 2016.

Comments: 14 pages, 1 table, 19 figures. Accepted for publication in MNRAS

Multicolor light curves simulations of Population III core-collapse supernovae: from shock breakout to $^{56}$Co decay

Authors: Alexey Tolstov, Ken'ichi Nomoto, Nozomu Tominaga, Miho Ishigaki, Sergey Blinnikov, Tomoharu Suzuki

Abstract: The properties of the first generation of stars and their supernova (SN) explosions remains unknown due to the lack of their actual observations. Recently many transient surveys are conducted and the feasibility of the detection of supernovae (SNe) of Pop III stars is growing. In this paper we study the multicolor light curves for a number of metal-free core-collapse SN models (25-100 M$_{\odot}$)… ▽ More The properties of the first generation of stars and their supernova (SN) explosions remains unknown due to the lack of their actual observations. Recently many transient surveys are conducted and the feasibility of the detection of supernovae (SNe) of Pop III stars is growing. In this paper we study the multicolor light curves for a number of metal-free core-collapse SN models (25-100 M$_{\odot}$) to provide the indicators for finding and identification of first generation SNe. We use mixing-fallback supernova explosion models which explain the observed abundance patterns of metal poor stars. Numerical calculations of the multicolor light curves are performed using multigroup radiation hydrodynamic code STELLA. The calculated light curves of metal-free SNe are compared with non-zero metallicity models and several observed SNe. We have found that the shock breakout characteristics, the evolution of the photosphere's velocity, the luminosity, the duration and color evolution of the plateau - all the SN phases are helpful to estimate the parameters of SN progenitor: the mass, the radius, the explosion energy and the metallicity. We conclude that the multicolor light curves can be potentially used to identify first generation SNe in the current (Subaru/HSC) and future transient surveys (LSST, JWST). They are also suitable for identification of the low-metallicity SNe in the nearby Universe (PTF, Pan-STARRS, Gaia). △ Less

Submitted 27 February, 2016; v1 submitted 28 December, 2015; originally announced December 2015.

Comments: 15 pages, 19 figures; Accepted for publication in Astrophysical Journal, February 26, 2016

Type I Superluminous Supernovae as Explosions inside Non-Hydrogen Circumstellar Envelopes

Authors: Elena Sorokina, Sergei Blinnikov, Ken'ichi Nomoto, Robert Quimby, Alexey Tolstov

Abstract: A number of Type I (hydrogenless) superluminous supernova (SLSN) events have been discovered recently. However, their nature remains debatable. One of the most promising ideas is the shock-interaction mechanism, but only simplified semi-analytical models have been applied so far. We simulate light curves for several Type I SLSN (SLSN-I) models enshrouded by dense, non-hydrogen circumstellar envelo… ▽ More A number of Type I (hydrogenless) superluminous supernova (SLSN) events have been discovered recently. However, their nature remains debatable. One of the most promising ideas is the shock-interaction mechanism, but only simplified semi-analytical models have been applied so far. We simulate light curves for several Type I SLSN (SLSN-I) models enshrouded by dense, non-hydrogen circumstellar envelopes, using a multi-group radiation hydrodynamics code that predicts not only bolometric, but also multicolor light curves. We demonstrate that the bulk of SLSNe-I including those with relatively narrow light curves like SN 2010gx or broad ones like PTF09cnd can be explained by the interaction of the SN ejecta with he CS envelope, though the range of parameters for these models is rather wide. Moderate explosion energy ($\sim (2 - 4)\cdot 10^{51}$ ergs) is sufficient to explain both narrow and broad SLSN-I light curves, but ejected mass and envelope mass differ for those two cases. Only 5 to 10 $M_\odot$ of non-hydrogen material is needed to reproduce the light curve of SN 2010gx, while the best model for PTF09cnd is very massive: it contains almost $ 50 M_\odot $ in the CS envelope and only $ 5 M_\odot $ in the ejecta. The CS envelope for each case extends from 10 $R_\odot$ to $\sim 10^5R_\odot$ ($7\cdot 10^{15} $ cm), which is about an order of magnitude larger than typical photospheric radii of standard SNe near the maximum light. We briefly discuss possible ways to form such unusual envelopes. △ Less

Submitted 6 September, 2016; v1 submitted 3 October, 2015; originally announced October 2015.

Comments: 16 pages; ApJ, in press; 2 new figures added with new modeled light curves; the text is expanded and corrected according to referee's and editor's comments

Journal ref: ApJ 829, 17 (2016)

Shock wave structure in astrophysical flows with an account of photon transfer

Authors: Alexey Tolstov, Sergey Blinnikov, Shigehiro Nagataki, Ken'ichi Nomoto

Abstract: For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we consider one-dimensional (1D) special relativistic radiation hydrodynamics by solving the Boltzmann equation for radiative transfer. The structure of a radiativ… ▽ More For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we consider one-dimensional (1D) special relativistic radiation hydrodynamics by solving the Boltzmann equation for radiative transfer. The structure of a radiative shock is calculated for a number of shock tube problems, including strong shock waves, and relativistic- and radiation-dominated cases. Calculations are performed using an iterative technique that consistently solves the equations of relativistic hydrodynamics and relativistic comoving radiative transfer. A comparison of radiative transfer solutions with the Eddington approximation and the M1 closure is made. A qualitative analysis of moment equations for radiation is performed and the conditions for the existence of jump discontinuity for non-relativistic cases are investigated numerically. △ Less

Submitted 29 February, 2016; v1 submitted 3 December, 2014; originally announced December 2014.

Comments: 10 pages, 10 figures, 2 tables, Accepted for publication in Astrophysical Journal, July 23, 2015

Spectral and Polarization Properties of Photospheric Emission From Stratified Jets

Authors: Hirotaka Ito, Shigehiro Nagataki, Jin Matsumoto, Shiu-Hang Lee, Alexey Tolstov, Jirong Mao, Maria Dainotti, Akira Mizuta

Abstract: We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the je… ▽ More We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($θ$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%. △ Less

Submitted 24 May, 2014; originally announced May 2014.

Comments: 21 pages, 12 figures, accepted for publication in ApJ

Matter mixing in aspherical core-collapse supernovae: a search for possible conditions for conveying $^{56}$Ni into high velocity regions

Authors: Masaomi Ono, Shigehiro Nagataki, Hirotaka Ito, Shiu-Hang Lee, Jirong Mao, Masa-aki Hashimoto, Alexey Tolstov

Abstract: We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarif… ▽ More We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarify the key conditions for reproducing such high velocity of $^{56}$Ni, we revisit matter mixing in aspherical core-collapse supernova explosions. Explosions are initiated artificially by injecting thermal and kinetic energies around the interface between the iron core and the silicon-rich layer. Perturbations of 5% or 30% amplitude in the radial velocities are introduced at several points in time. We found that no high velocity $^{56}$Ni can be obtained if we consider bipolar explosions with perturbations (5% amplitude) of pre-supernova origins. If large perturbations (30% amplitude) are introduced or exist due to some unknown mechanism in a later phase just before the shock wave reaches the hydrogen envelope, $^{56}$Ni with a velocity of 3,000 km s$^{-1}$ can be obtained. Aspherical explosions that are asymmetric across the equatorial plane with clumpy structures in the initial shock waves are investigated. We found that the clump sizes affect the penetration of $^{56}$Ni. Finally, we report that an aspherical explosion model that is asymmetric across the equatorial plane with multiple perturbations of pre-supernova origins can cause the penetration of $^{56}$Ni clumps into fast moving matter of 3,000 km s$^{-1}$. We show that both aspherical explosion with clumpy structures and perturbations of pre-supernova origins may be necessary to reproduce the observed high velocity of $^{56}$Ni. To confirm this, more robust three-dimensional simulations are required. △ Less

Submitted 18 June, 2013; originally announced June 2013.

Comments: 29 pages, 28 figures, 2 tables, submitted in ApJ

Coupling of matter and radiation at supernova shock breakout

Authors: A. G. Tolstov, S. I. Blinnikov, D. K. Nadyozhin

Abstract: Some features of the physics of radiation-dominated shock waves are discussed with emphasis on the peculiarities which are important for correct numerical modeling of shock breakouts in supernova. With account of those peculiarities, a number of models for different supernova types is constructed based on multigroup radiation transfer coupled to hydrodynamics. We describe the implementation of a n… ▽ More Some features of the physics of radiation-dominated shock waves are discussed with emphasis on the peculiarities which are important for correct numerical modeling of shock breakouts in supernova. With account of those peculiarities, a number of models for different supernova types is constructed based on multigroup radiation transfer coupled to hydrodynamics. We describe the implementation of a new algorithm RADA, designed for modeling photon transfer at extremely-relativistic motions of matter, into our older code STELLA. The results of numerical simulations of light curves, and continuum spectra are presented. The influence of effects of photon scattering on electrons, of thermalization depth and of special relativity in transfer equation is considered. Some cases are demonstrated, when the appearance of hard X-ray emission is possible at the shock breakout. The necessary refinements in numerical algorithms for radiative transfer and hydrodynamics are pointed out. Prospects for using the results of numerical simulation to analyze and interpret available and future data from space observatories are discussed. △ Less

Submitted 15 December, 2012; originally announced December 2012.

Comments: Accepted for publication in MNRAS 2012 December 06, 20 pages, 30 figures

The influence of jet geometry on light curves and spectra of GRB afterglows

Authors: A. G. Tolstov

Abstract: We have performed detailed calculations of spectra and light curves of GRB afterglows assuming that the observed GRBs can have a jet geometry. The calculations are based on an expanding relativistic shock GRB afterglow model where the afterglow is the result of synchrotron radiation of relativistic electrons with power-law energy distribution at the front of external shock being decelerated in a c… ▽ More We have performed detailed calculations of spectra and light curves of GRB afterglows assuming that the observed GRBs can have a jet geometry. The calculations are based on an expanding relativistic shock GRB afterglow model where the afterglow is the result of synchrotron radiation of relativistic electrons with power-law energy distribution at the front of external shock being decelerated in a circumstellar medium. To determine the intensity on the radiation surface we solve numerically the full time-, angle-, and frequency-dependent special relativistic transfer equation in the comoving frame using the method of long characteristics. △ Less

Submitted 28 June, 2012; originally announced June 2012.

Spectra and Light Curves of GRB Afterglows

Authors: A. G. Tolstov, S. I. Blinnikov

Abstract: We performed accurate numerical calculations of angle-, time-, and frequency-dependent radiative transfer for the relativistic motion of matter in gamma-ray burst (GRB) models. Our technique for solving the transfer equation, which is based on the method of characteristics, can be applied to the motion of matter with a Lorentz factor up to 1000. The effect of synchrotron self-absorption is taken… ▽ More We performed accurate numerical calculations of angle-, time-, and frequency-dependent radiative transfer for the relativistic motion of matter in gamma-ray burst (GRB) models. Our technique for solving the transfer equation, which is based on the method of characteristics, can be applied to the motion of matter with a Lorentz factor up to 1000. The effect of synchrotron self-absorption is taken into account. We computed the spectra and light curves from electrons with a power-law energy distribution in an expanding relativistic shock and compare them with available analytic estimates. The behavior of the optical afterglows from GRB 990510 and GRB 000301c is discussed qualitatively. △ Less

Submitted 3 February, 2004; originally announced February 2004.

Comments: 8 pages, 7 figures

Journal ref: Astron.Lett. 29 (2003) 353-360