Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userProgenitor properties of type II supernovae: fitting to hydrodynamical models using Markov chain Monte Carlo methods
82
0
0.0
(
0
)
Authors
L. Martinez
Ask ChatGPT about the research
No Arabic abstract
The progenitor and explosion properties of type II supernovae (SNe II) are fundamental to understand the evolution of massive stars. Special interest has been given to the range of initial masses of their progenitors, but despite the efforts made, it is still uncertain. Direct imaging of progenitors in pre-explosion images point out an upper initial mass cutoff of $sim$18$M_{odot}$. However, this is in tension with previous studies in which progenitor masses inferred by light curve modelling tend to favour high-mass solutions. Moreover, it has been argued that light curve modelling alone cannot provide a unique solution for the progenitor and explosion properties of SNe II. We develop a robust method which helps us to constrain the physical parameters of SNe II by fitting simultaneously their bolometric light curve and the evolution of the photospheric velocity to hydrodynamical models using statistical inference techniques. Pre-supernova red supergiant models were created using the stellar evolution code MESA, varying the initial progenitor mass. The explosion of these progenitors was then processed through hydrodynamical simulations, where the explosion energy, synthesised nickel mass, and the latters spatial distribution within the ejecta were changed. We compare to observations via Markov chain Monte Carlo methods. We apply this method to a well-studied set of SNe with an observed progenitor in pre-explosion images and compare with results in the literature. Progenitor mass constraints are found to be consistent between our results and those derived by pre-SN imaging and the analysis of late-time spectral modelling. We have developed a robust method to infer progenitor and explosion properties of SN II progenitors which is consistent with other methods in the literature, which suggests that hydrodynamical modelling is able to accurately constrain physical properties of SNe II.
rate research
Read More
Nebular phase spectra of core-collapse supernovae (SNe) provide critical and unique information on the progenitor massive star and its explosion. We present a set of 1-D steady-state non-local thermodynamic equilibrium radiative transfer calculations of type II SNe at 300d after explosion. Guided by results for a large set of stellar evolution simulations, we craft ejecta models for type II SNe from the explosion of a 12, 15, 20, and 25Msun star. The ejecta density structure and kinetic energy, the 56Ni mass, and the level of chemical mixing are parametrized. Our model spectra are sensitive to the adopted line Doppler width, a phenomenon we associate with the overlap of FeII and OI lines with Lyalpha and Lybeta. Our spectra show a strong sensitivity to 56Ni mixing since it determines where decay power is absorbed. Even at 300d after explosion, the H-rich layers reprocess the radiation from the inner metal rich layers. In a given progenitor model, variations in 56Ni mass and distribution impact the ejecta ionization, which can modulate the strength of all lines. Such ionization shifts can quench CaII line emission. In our set of models, the OI6300 doublet strength is the most robust signature of progenitor mass. However, we emphasize that convective shell merging in the progenitor massive star interior can pollute the O-rich shell with Ca, which will weaken the OI6300 doublet flux in the resulting nebular SN II spectrum. This process may occur in Nature, with a greater occurrence in higher mass progenitors, and may explain in part the preponderance of progenitor masses below 17Msun inferred from nebular spectra.
We have a developed a new method for fitting spectral energy distributions (SEDs) to identify and constrain the physical properties of high-redshift (4 < z < 8) galaxies. Our approach uses an implementation of Bayesian based Markov Chain Monte Carlo (PiMC^2) that allows us to compare observations to arbitrarily complex models and to compute 95% credible intervals that provide robust constraints for the model parameters. The work is presented in 2 sections. In the first, we test PiMC^2 using simulated SEDs to not only confirm the recovery of the known inputs but to assess the limitations of the method and identify potential hazards of SED fitting when applied specifically to high redshift (z>4) galaxies. Our tests reveal five critical results: 1) the ability to confidently constrain metallicity, population ages, and Av all require photometric accuracy better than what is currently achievable (i.e. less than a few percent); 2) the ability to confidently constrain stellar masses (within a factor of two) can be achieved without the need for high-precision photometry; 3) the addition of IRAC photometry does not guarantee that tighter constraints of the stellar masses and ages can be defined; 4) different assumptions about the star formation history can lead to significant biases in mass and age estimates; and 5) we are able to constrain stellar age and Av of objects that are both young and relatively dust free. In the second part of the paper we apply PiMC^2 to 17 4<z<8 objects, including the GRAPES Ly alpha sample (4<z<6), supplemented by HST/WFC3 near-IR observations, and several broad band selected z>6 galaxies. Using PiMC^2, we are able to constrain the stellar mass of these objects and in some cases their stellar age and find no evidence that any of these sources formed at a redshift much larger than z_f=8, a time when the Universe was ~ 0.6 Gyr old.
We analyze the properties of 42 rapidly rotating, low metallicity, quasi-chemically homogeneously evolving stellar models in the mass range between 4 and 45 $,mathrm{M}_odot$ at the time of core collapse. Such models were proposed as progenitors for both superluminous supernovae (SLSNe) and long duration gamma-ray bursts (lGRBs), and the Type Ic-BL supernovae (SNe) that are associated with them. Our findings suggest that whether these models produce a magnetar driven SLSN explosion or a near-critically rotating black hole (BH) is not a monotonic function of the initial mass. Rather, their explodability varies non-monotonically depending on the late core evolution, once chemical homogeneity is broken. Using different explodability criteria we find that our models have a clear preference to produce SLSNe at lower masses, and lGRBs at higher masses; but find several exceptions, expecting lGRBs to form from stars as low as 10 $,mathrm{M}_odot$, and SLSNe with progenitors as massive as 30 $,mathrm{M}_odot$. In general, our models reproduce the predicted angular momenta, ejecta masses and magnetic field strengths at core collapse inferred for SLSNe and lGRBs, and suggest significant interaction with their circumstellar medium, particularly for explosions with low ejecta mass.
In this paper, we study the asymptotic variance of sample path averages for inhomogeneous Markov chains that evolve alternatingly according to two different $pi$-reversible Markov transition kernels $P$ and $Q$. More specifically, our main result allows us to compare directly the asymptotic variances of two inhomogeneous Markov chains associated with different kernels $P_i$ and $Q_i$, $iin{0,1}$, as soon as the kernels of each pair $(P_0,P_1)$ and $(Q_0,Q_1)$ can be ordered in the sense of lag-one autocovariance. As an important application, we use this result for comparing different data-augmentation-type Metropolis-Hastings algorithms. In particular, we compare some pseudo-marginal algorithms and propose a novel exact algorithm, referred to as the random refreshment algorithm, which is more efficient, in terms of asymptotic variance, than the Grouped Independence Metropolis-Hastings algorithm and has a computational complexity that does not exceed that of the Monte Carlo Within Metropolis algorithm.
We present JEKYLL, a new code for modelling of supernova (SN) spectra and lightcurves based on Monte-Carlo (MC) techniques for the radiative transfer. The code assumes spherical symmetry, homologous expansion and steady state for the matter, but is otherwise capable of solving the time-dependent radiative transfer problem in non-local-thermodynamic-equilibrium (NLTE). The method used was introduced in a series of papers by Lucy, but the full time-dependent NLTE capabilities of it have never been tested. Here, we have extended the method to include non-thermal excitation and ionization as well as charge-transfer and two-photon processes. Based on earlier work, the non-thermal rates are calculated by solving the Spencer-Fano equation. Using a method previously developed for the SUMO code, macroscopic mixing of the material is taken into account in a statistical sense. In addition, a statistical Markov-chain model is used to sample the emission frequency, and we introduce a method to control the sampling of the radiation field. Except for a description of JEKYLL, we provide comparisons with the ARTIS, SUMO and CMFGEN codes, which show good agreement in the calculated spectra as well as the state of the gas. In particular, the comparison with CMFGEN, which is similar in terms of physics but uses a different technique, shows that the Lucy method does indeed converge in the time-dependent NLTE case. Finally, as an example of the time-dependent NLTE capabilities of JEKYLL, we present a model of a Type IIb SN, taken from a set of models presented and discussed in detail in an accompanying paper. Based on this model we investigate the effects of NLTE, in particular those arising from non-thermal excitation and ionization, and find strong effects even on the bolometric lightcurve. This highlights the need for full NLTE calculations when simulating the spectra and lightcurves of SNe.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
