Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAxisymmetric Radiative Transfer Models of Kilonovae
80
0
0.0
(
0
)
Authors
Oleg Korobkin
Ask ChatGPT about the research
No Arabic abstract
The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g. in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide curtaining, which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent {P Cygni} features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.
rate research
Read More
The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergo rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a ``kilonova, which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ~hours-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.
The precise location of the water ice condensation front (snow line) in the protosolar nebula has been a debate for a long time. Its importance stems from the expected substantial jump in the abundance of solids beyond the snow line, which is conducive to planet formation, and from the higher stickiness in collisions of ice-coated dust grains, which may help the process of coagulation of dust and the formation of planetesimals. In an optically thin nebula, the location of the snow line is easily calculated to be around 3 AU. However, in its first 5 to 10 million years, the solar nebula was optically thick, implying a smaller snow line radius due to shielding from direct sunlight, but also a larger radius because of viscous heating. Several models have attempted to treat these opposing effects. However, until recently treatments beyond an approximate 1+1D radiative transfer were unfeasible. We revisit the problem with a fully self-consistent 3D treatment in an axisymmetric disk model, including a density-dependent treatment of the dust and ice sublimation. We find that the location of the snow line is very sensitive to the opacities of the dust grains and the mass accretion rate of the disk. We show that previous approximate treatments are quite efficient at determining the location of the snow line if the energy budget is locally dominated by viscous accretion. We derive an analytic estimate of the location of the snow line that compares very well with results from this and previous studies. Using solar abundances of the elements we compute the abundance of dust and ice and find that the expected jump in solid surface density at the snow line is smaller than previously assumed. We further show that in the inner few AU the refractory species are partly evaporated, leading to a significantly smaller solid state surface density in the regions where the rocky planets were formed.
There are two classes of viable progenitors for normal Type Ia supernovae (SNe Ia): systems in which a white dwarf explodes at the Chandrasekhar mass ($M_{ch}$), and systems in which a white dwarf explodes below the Chandrasekhar mass (sub-$M_{ch}$). It is not clear which of these channels is dominant; observations and light curve modeling have provided evidence for both. Here we use an extensive grid of 4500 time-dependent, multiwavelength radiation transport simulations to show that the sub-$M_{ch}$ model can reproduce the entirety of the width-luminosity relation (WLR), while the $M_{ch}$ model can only produce the brighter events $(0.8 < Delta M_{15}(B) < 1.55)$, implying that fast-declining SNe Ia come from sub-$M_{ch}$ explosions. We do not assume a particular theoretical paradigm for the progenitor or explosion mechanism, but instead construct parameterized models that vary the mass, kinetic energy, and compositional structure of the ejecta, thereby realizing a broad range of possible outcomes of white dwarf explosions. We provide fitting functions based on our large grid of detailed simulations that map observable properties of SNe Ia such as peak brightness and light curve width to physical parameters such as $^{56}mathrm{Ni}$ and total ejected mass. These can be used to estimate the physical properties of observed SNe Ia.
Accreting supermassive black holes are sources of polarized radiation that propagates through highly curved spacetime before reaching the observer. In order to help interpret observations of such polarized emission, accurate and efficient numerical schemes for polarized radiative transfer in curved spacetime are needed. In this manuscript we extend our publicly available radiative transfer code RAPTOR to include polarization. We provide a brief review of different codes and methods for covariant polarized radiative transfer available in the literature and existing codes, and present an efficient new scheme. For the spacetime-propagation aspect of the computation, we develop a compact, Lorentz-invariant representation of a polarized ray. For the plasma-propagation aspect of the computation, we perform a formal analysis of the stiffness of the polarized radiative-transfer equation with respect to our explicit integrator, and develop a hybrid integration scheme that switches to an implicit integrator in case of stiffness, in order to solve the equation with optimal speed and accuracy for all possible values of the local optical/Faraday thickness of the plasma. We perform a comprehensive code verification by solving a number of well-known test problems using RAPTOR and comparing its output to exact solutions. We also demonstrate convergence with existing polarized radiative-transfer codes in the context of complex astrophysical problems. RAPTOR is capable of performing polarized radiative transfer in arbitrary, highly curved spacetimes. This capability is crucial for interpreting polarized observations of accreting black holes, which can yield information about the magnetic-field configuration in such accretion flows. The efficient formalism implemented in RAPTOR is computationally light and conceptually simple. The code is publicly available.
We use a Monte Carlo radiative transfer model (MCRTM) to simulate the UBVRI light curves, images and linear polarization of a light echo from supernova SN$~$1987A in the Large Magellanic Cloud (LMC) using various dust cloud shapes, sizes, and optical properties. We compare the theoretical simulations to the observations of AT2019xis, a light echo detected at a large angular distance (4.05$^{}$) from SN$~$1987A. We estimate the size and optical thickness of the dust cloud based on the simulation results and the observations of Optical Gravitational Lensing Experiment (OGLE-IV) Transient Detection System (OTDS) I-band light curve. The mass of the dust cloud is calculated using the estimated size, optical thickness and extinction coefficient. If the dust cloud is assumed to correspond to a gas-to-dust ratio of 300, the total mass of the dust cloud is approximately 7.8-9.3 $M_{odot}$. Based on these theoretical models, we show that the morphological shapes of the light echoes in the wavelength range in or shorter than the U-band to be very different from those in the longer wavelength bands, and the difference carries important information on the early UV radiation of SN$~$1987A.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
