Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRadiative-transfer modeling of supernovae in the nebular-phase. A novel treatment of chemical mixing in spherical symmetry
204
0
0.0
(
0
)
Authors
Luc Dessart
Ask ChatGPT about the research
No Arabic abstract
Supernova (SN) explosions, through the metals they release, play a pivotal role in the chemical evolution of the Universe and the origin of life. Nebular phase spectroscopy constrains such metal yields, for example through forbidden line emission associated with OI, CaII, FeII, or FeIII. Fluid instabilities during the explosion produce a complex 3D ejecta structure, with considerable macroscopic, but no microscopic, mixing of elements. This structure sets a formidable challenge for detailed nonlocal thermodynamic equilibrium radiative transfer modeling, which is generally limited to 1D in grid-based codes. Here, we present a novel and simple method that allows for macroscopic mixing without any microscopic mixing, thereby capturing the essence of mixing in SN explosions. With this new technique, the macroscopically mixed ejecta is built by shuffling in mass space, or equivalently in velocity space, the shells from the unmixed coasting ejecta. The method requires no change to the radiative transfer, but necessitates high spatial resolution to resolve the rapid variation in composition with depth inherent to this shuffled-shell structure. We show results for a few radiative-transfer simulations for a Type II SN explosion from a 15Msun progenitor star. Our simulations capture the strong variations in temperature or ionization between the various shells that are rich in H, He, O, or Si. Because of nonlocal energy deposition, gamma rays permeate through an extended region of the ejecta, making the details of the shell arrangement unimportant. The greater physical consistency of the method delivers spectral properties at nebular times that are more reliable, in particular in terms of individual emission line strengths, which may serve to constrain the SN yields and, for core collapse SNe, the progenitor mass. The method works for all SN types.
rate research
Read More
We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.
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.
The discontinuous Galerkin finite element method (DG-FEM) is successfully applied to treat a broad variety of transport problems numerically. In this work, we use the full capacity of the DG-FEM to solve the radiative transfer equation in spherical symmetry. We present a discontinuous Galerkin method to directly solve the spherically-symmetric radiative transfer equation as a two-dimensional problem. The transport equation in spherical atmospheres is more complicated than in the plane-parallel case due to the appearance of an additional derivative with respect to the polar angle. The DG-FEM formalism allows for the exact integration of arbitrarily complex scattering phase functions, independent of the angular mesh resolution. We show that the discontinuous Galerkin method is able to describe accurately the radiative transfer in extended atmospheres and to capture discontinuities or complex scattering behaviour which might be present in the solution of certain radiative transfer tasks and can, therefore, cause severe numerical problems for other radiative transfer solution methods.
Type II supernovae (SNe) often exhibit a linear polarization, arising from free-electron scattering, with complicated optical signatures, both in the continuum and in lines. Focusing on the early nebular phase, at a SN age of 200d, we conduct a systematic study of the polarization signatures associated with a 56Ni `blob that breaks spherical symmetry. Our ansatz, supported by nonLTE radiative transfer calculations, is that the primary role of such a 56Ni blob is to boost the local density of free electrons, which is otherwise reduced following recombination in SNe II. Using 2D polarized radiation transfer modeling, we explore the influence of such an electron-density enhancement, varying its magnitude N_e_fac, its velocity location V_blob, and its spatial extent. For plausible N_e_fac values of a few tens, a high-velocity blob can deliver a continuum polarization P_cont of 0.5-1.0% at 200d. Our simulations reproduce the analytic scalings for P_cont, and in particular the linear growth with the blob radial optical depth. The most constraining information is, however, carried by polarized line photons. For a high V_blob, the polarized spectrum appears as a replica of the full spectrum, scaled down by a factor 100 to 1000 (i.e., 1/P_cont), and redshifted by an amount V_blob(1-cos(alpha_los)), where alpha_los is the line of sight angle. As V_blob is reduced, the redshift decreases and the replication deteriorates. Lines whose formation region overlap with the blob appear weaker and narrower in the polarized flux. Because of its dependence on inclination (~ sin^2 alpha_los), the polarization preferentially reveals asymmetries in the plane perpendicular to the line of sight. With the adequate choice of electron-density enhancement, some of these results may apply to asymmetric explosions in general, or to the polarization signatures from newly-formed dust in the outer ejecta. [abridged]
The recent identification of the first complex chiral molecule, propylene oxide (PrO) in space opens up a new window to further study the origin of homochirality on the Earth. There are some recent studies to explain the formation of PrO however additional studies on the formation of this species are needed for better understanding. We seek to prepare a complete reaction network to study the formation of propylene oxide in the astrophysically relevant conditions. Based on our results, a detailed radiative transfer modeling has been carried out to propose some more transitions which would potentially be targeted in the millimeter wave domain. Gas-grain chemical network was used to explain the observed abundance of PrO in a cold shell surrounding the high-mass star-forming region of Sgr B2. Quantum chemical calculations were employed to study various reaction parameters and to compute multiple vibrational frequencies of PrO. To model the formation of PrO in the observed region, we considered a dark cloud model. Additionally, we used a model to check the feasibility of forming PrO in the hot core region. Some potential transitions in the millimeter wave domain are predicted which could be useful for the future astronomical detection. Radiative transfer modeling has been utilized to extract the physical condition which might be useful to know the properties of the source in detail. Moreover, vibrational transitions of PrO has been provided which could be very useful for the future detection of PrO by the upcoming James Webb Space Telescope (JWST).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
