ABSTRACT
The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These ...observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy r-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound–bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion–opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and thus eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion–opacity formalism. The agreement between the line-binned and expansion–opacity results is explained as arising from the similarity in their opacities in the limit of low optical depth, where radiation transport is important in the ejecta. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semirelativistic approaches.
Abstract
Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The ...characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of two-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1
M
⊙
and 0.05 to 0.3
c
, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (
Y
e
) ejecta with a robust
r
-process composition and either a spherical or lobed high-
Y
e
ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code
SuperNu
using a full suite of lanthanide and fourth-row element opacities. We examine the trends of these models with parameter variation, show how they can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.
ABSTRACT
We report on our observing campaign of the compact binary merger GW190814, detected by the Advanced LIGO and Advanced Virgo detectors on 2019 August 14. This signal has the best localization ...of any observed gravitational wave (GW) source, with a 90 per cent probability area of 18.5 deg2, and an estimated distance of ≈240 Mpc. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering 88 per cent of the probability area down to a limiting magnitude of w = 19.9 AB. Nearby galaxies within the high probability region were targeted with the Lowell Discovery Telescope (LDT), whereas promising candidate counterparts were characterized through multicolour photometry with the Reionization and Transients InfraRed (RATIR) and spectroscopy with the Gran Telescopio de Canarias (GTC). We use our optical and near-infrared limits in conjunction with the upper limits obtained by the community to constrain the possible electromagnetic counterparts associated with the merger. A gamma-ray burst seen along its jet’s axis is disfavoured by the multiwavelength data set, whereas the presence of a burst seen at larger viewing angles is not well constrained. Although our observations are not sensitive to a kilonova similar to AT2017gfo, we can rule out high-mass (>0.1 M⊙) fast-moving (mean velocity ≥0.3c) wind ejecta for a possible kilonova associated with this merger.
Abstract
Recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code ...to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a “flashlight” effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of
ζ
≈ 10
−16
s
−1
; however, within
r
≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.
The collapsar engine behind long-duration gamma-ray bursts extracts the energy released from the rapid accretion of a collapsing star onto a stellar mass black hole. In a collapsing star, this black ...hole can form in two ways: the direct collapse of the stellar core into a black hole and the delayed collapse of a black hole caused by fallback in a weak supernova explosion. In the case of a delayed-collapse black hole, the strong collapsar-driven explosion overtakes the weak supernova explosion before shock breakout, and it is very difficult to distinguish this black hole formation scenario from the direct-collapse scenario. However, the delayed-collapse mechanism, with its double explosion, produces explosive nucleosynthetic yields that are very different from those in the direct-collapse scenario. We present one-dimensional studies of the nucleosynthetic yields from both black hole formation scenarios, deriving differences and trends in their nucleosynthetic yields.
Multidimensional simulations of the neutrino-driven mechanism behind core-collapse supernovae have long shown that the explosions from this mechanism would be asymmetric. Recently, detailed ...core-collapse simulations have shown that the explosion may be strongest in a single direction. We present a suite of simulations modeling these "single-lobe" supernova explosions of a 15 M sub( )red supergiant star, focusing on the effect these asymmetries have on the gamma-ray emission and the mixing in the explosion. We discuss how these asymmetries in the explosion mechanism might explain many of the observed "asymmetries" of supernovae, focusing on features of both supernova 1987A and the Cassiopeia A supernova remnant. In particular, we show that single-lobe explosions provide a promising solution to the redshifted iron lines of supernova 1987A. We also show that the extent of mixing for explosive burning products depends sensitively on the angular profile of the velocity asymmetry and may be much more extensive than previously assumed.
Detailed radiative transfer simulations of kilonovae are difficult to apply directly to observations; they only sparsely cover simulation parameters, such as the mass, velocity, morphology, and ...composition of the ejecta. On the other hand, semianalytic models for kilonovae can be evaluated continuously over model parameters, but neglect important physical details which are not incorporated in the simulations, thus introducing systematic bias. Starting with a grid of two-dimensional anisotropic simulations of kilonova light curves covering a wide range of ejecta properties, we apply adaptive learning techniques to iteratively choose new simulations and produce high-fidelity surrogate models for those simulations. These surrogate models allow for continuous evaluation across model parameters while retaining the microphysical details about the ejecta. Using a code formultimessenger inference developed by our group, we demonstrate how to use our interpolated models to infer kilonova parameters. Comparing to inferences using simplified analytic models, we recover different ejecta properties. We discuss the implications of this analysis which is qualitatively consistent with similar previous work using detailed ejecta opacity calculations and which illustrates systematic challenges for kilonova modeling. An associated data and code release provides our interpolated light-curve models, interpolation implementation which can be applied to reproduce our work or extend to new models, and our multimessenger parameter inference engine.
We report on the use of the Los Alamos suite of relativistic atomic physics codes to generate radiative opacities for the modeling of astrophysically relevant plasmas under local thermodynamic ...equilibrium (LTE) conditions. The atomic structure calculations are carried out in fine-structure detail, including full configuration interaction. Three example applications are considered: iron opacities at conditions relevant to the base of the solar convection zone, nickel opacities for the modeling of stellar envelopes, and samarium opacities for the modeling of light curves produced by neutron star mergers. In the first two examples, comparisons are made between opacities that are generated with the fully and semi-relativistic capabilities in the Los Alamos suite of codes. As expected for these highly charged, iron-peak ions, the two methods produce reasonably similar results, providing confidence that the numerical methods have been correctly implemented. However, discrepancies greater than 10% are observed for nickel and investigated in detail. In the final application, the relativistic capability is used in a preliminary investigation of the complicated absorption spectrum associated with cold lanthanide elements.
We retrospectively reviewed the results of operative decompression for peroneal nerve palsy in thirty-one patients who had been managed between 1980 and 1990. All patients had been managed ...non-operatively for at least two months after they had initially been seen. Intraoperatively, we found epineurial fibrosis and bands of fibrous tissue constricting the peroneal nerve at the level of the fibular head and at the proximal origin of the peroneus longus muscle. At a mean of thirty-six months (range, twelve to seventy-two months) postoperatively, thirty (97 per cent) of the thirty-one patients reported subjective and functional improvement and were able to discontinue the use of the ankle-foot orthosis. In contrast, only three of nine patients who had been managed non-operatively reported subjective and functional improvement (p < 0.01). Peroneal nerve palsy does not always resolve spontaneously; if it is left untreated, the loss of dorsiflexion of the ankle and persistent paresthesias can result in severe functional disability. Therefore, if non-operative measures do not lead to improvement within two months, we believe that operative decompression should be considered.
The currently favored model for long-duration gamma-ray bursts (GRBs) invokes explosions from the collapse of a massive star down to a black hole, either directly or through fallback. Those GRBs ...forming via fallback will produce much less radioactive nickel, and hence it has been argued (without any real calculation) that these systems produce dim supernovae. These fallback black hole GRBs have recently been argued as possible progenitors of a newly discovered set of GRBs lacking any associated supernovae. Here we present the first ever radiation-hydrodynamics calculations of the light curves produced in the hypernova explosion by a delayed-fallback gamma-ray burst. We find that the bolometric light curve is dominated by shock-deposited energy, not the decay of radioactive elements. As such, observations of such bursts actually probe the density in the progenitor wind more than it does the production of radioactive nickel.