Light Curves and Spectra from Kilonova Models [Poster]
Kilonovae are the radioactively powered optical/IR transients following a binary neutron star or neutron star black hole merger.
79 ASTRONOMY AND ASTROPHYSICS↗
SEARCH · Engineering Papers
Results for “binary stars”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Kilonova Detection
The first observation of gravitational waves (GW) from a binary neutron star collision, known as GW170817, was detected by Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo GW observatories. This event marked a profound discovery of the kilonova, which occurs when two neutron stars collide, resulting in an energetically explosive release of gamma-ray bursts (GRBs) and emitting bright electromagnetic radiation. This project focuses on assessing how well kilonovae will be detected in the Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory. By using Kasen models with varying ejecta mass, ejecta velocity, and lanthanide fraction, we analyze these 329 models and create kilonova light curves to assess the number of visible kilonovae for each band, based on the minimum magnitude at various redshifts.
Khalid, Anna↗
The All-sky Medium Energy Gamma-ray Observatory eXplorer (AMEGO-X) Mission Concept
The All-sky Medium Energy Gamma-ray Observatory eXplorer (AMEGO-X) is designed to identify and characterize gamma rays from extreme explosions and accelerators. The main science themes include: supermassive black holes and their connections to neutrinos and cosmic rays; binary neutron star mergers and the relativistic jets they produce; cosmic ray particle acceleration sources including Galactic supernovae; and continuous monitoring of other astrophysical events and sources over the full sky in this important energy range. AMEGO-X will probe the medium energy gamma-ray band using a single instrument with sensitivity up to an order of magnitude greater than previous telescopes in the energy range 100 keV to 1 GeV that can be only realized in space. During its three-year baseline mission, AMEGO-X will observe nearly the entire sky every two orbits, building up a sensitive all-sky map of gamma-ray sources and emission. AMEGO-X was submitted in the recent 2021 NASA MIDEX Announcement of Opportunity.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
High deuteron and neutron yields from the interaction of a petawatt laser with a cryogenic deuterium jet
A compact high-flux, short-pulse neutron source would have applications from nuclear astrophysics to cancer therapy. Laser-driven neutron sources can achieve fluxes much higher than spallation and reactor neutron sources by reducing the volume and time in which the neutron-producing reactions occur by orders of magnitude. We report progress towards an efficient laser-driven neutron source in experiments with a cryogenic deuterium jet on the Texas Petawatt laser. Neutrons were produced both by laser-accelerated multi-MeV deuterons colliding with Be and mixed metallic catchers and by d (d,n) 3 He fusion reactions within the jet. We observed deuteron yields of 10 13 /shot in quasi-Maxwellian distributions carrying ~8-10% of the input laser energy. We obtained neutron yields greater than 10 10 /shot and found indications of a deuteron-deuteron fusion neutron source with high peak flux (>10 22 cm -2 s -1 ). The estimated fusion neutron yield in our experiment is one order of magnitude higher than any previous laser-induced dd fusion reaction. Though many technical challenges will have to be overcome to convert this proof-of-principle experiment into a consistent ultra-high flux neutron source, the neutron fluxes achieved here suggest laser-driven neutron sources can support laboratory study of the rapid neutron-capture process, which is otherwise thought to occur only in astrophysical sites such as core-collapse supernova, and binary neutron star mergers.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Gamma Ray Source Localization for Time Projection Chamber Telescopes Using Convolutional Neural Networks
Diverse phenomena such as positron annihilation in the Milky Way, merging binary neutron stars, and dark matter can be better understood by studying their gamma ray emission. Despite their importance, MeV gamma rays have been poorly explored at sensitivities that would allow for deeper insight into the nature of the gamma emitting objects. In response, a liquid argon time projection chamber (TPC) gamma ray instrument concept called GammaTPC has been proposed and promises exploration of the entire sky with a large field of view, large effective area, and high polarization sensitivity. Optimizing the pointing capability of this instrument is crucial and can be accomplished by leveraging convolutional neural networks to reconstruct electron recoil paths from Compton scattering events within the detector. In this investigation, we develop a machine learning model architecture to accommodate a large data set of high fidelity simulated electron tracks and reconstruct paths. We create two model architectures: one to predict the electron recoil track origin and one for the initial scattering direction. We find that these models predict the true origin and direction with extremely high accuracy, thereby optimizing the observatory’s estimates of the sky location of gamma ray sources.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics
The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.
79 ASTRONOMY AND ASTROPHYSICS↗
Spectroscopic Orbits of 11 Nearby, Mid-to-late M-dwarf Binaries
We present the spectroscopic orbits of 11 nearby, mid-to-late M dwarf binary systems in a variety of configurations: 2 single-lined binaries (SB1s), 7 double-lined binaries (SB2s), 1 double-lined triple (ST2), and 1 triple-lined triple (ST3). Eight of these orbits are the first published for these systems, while five are newly identified multiples. We obtained multi-epoch, high-resolution spectra with the TRES instrument on the 1.5 m Tillinghast Reflector at the Fred Lawrence Whipple Observatory located on Mt. Hopkins in AZ. Using the TiO molecular bands at 7065–7165 Å, we calculated radial velocities for these systems, from which we derived their orbits. We find LHS 1817 to have in a 7 hr period a companion that is likely a white dwarf, due to the ellipsoidal modulation we see in our MEarth-North light-curve data. We find G 123-45 and LTT 11586 to host companions with minimum masses of 41 M {sub Jup} and 44 M {sub Jup} with orbital periods of 35 and 15 days, respectively. We find 2MA 0930+0227 to have a rapidly rotating stellar companion in a 917 day orbital period. GJ 268, GJ 1029, LP 734-34, GJ 1182, G 258-17, and LTT 7077are SB2s with stellar companions with orbital periods of 10, 96, 34, 154, 5, and 84 days; LP 655-43 is an ST3 with one companion in an 18 day orbital period and an outer component in a longer undetermined period. In addition, we present radial velocities for both components of L 870-44AB and for the outer components of LTT 11586 and LP 655-43.
79 ASTRONOMY AND ASTROPHYSICS↗
Photometric Study and Absolute Parameter Estimation of Six Totally Eclipsing Contact Binaries
High precision CCD observations of six totally eclipsing contact binaries were presented and analyzed. It is found that only one target is an A-type contact binary (V429 Cam), while the others are W-type contact ones. By analyzing the times of light minima, we discovered that two of them exhibit secular period increase while three manifest long-term period decrease. For V1033 Her, a cyclic variation superimposed on the long-term increase was discovered. By comparing the Gaia distances with those calculated by the absolute parameters of 173 contact binaries, we found that the Gaia distance can be applied to estimate the absolute parameters for most contact binaries. The absolute parameters of our six targets were estimated using their Gaia distances. The evolutionary status of contact binaries was studied, and we found that the A- and W-subtype contact binaries may have different formation channels. The relationship between the spectroscopic and photometric mass ratios for 101 contact binaries was presented. It was discovered that the photometric mass ratios are in good agreement with the spectroscopic ones for almost all of the totally eclipsing systems, which is corresponding to the results derived by Pribulla et al. and Terrell & Wilson.
47 OTHER INSTRUMENTATION↗
Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows
Here we present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature thereby integrating over the thermal seed photon flux from the disk surface. The method either assumes the ion and electron temperatures are equal (1T) or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. This new approach leads to radiative efficiency values sime65% above Novikov–Thorne, with a larger fraction of total cooling in the corona as compared to simulations performed using the original target-temperature cooling function. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach, which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).
79 ASTRONOMY AND ASTROPHYSICS↗
Modeling Kilonova Light Curves: Dependence on Nuclear Inputs
The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing β-decay, α-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.
79 ASTRONOMY AND ASTROPHYSICS↗
Observation of the Gamma-Ray Binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS Telescopes
The results of gamma-ray observations of the binary system HESS J0632 + 057 collected during 450 hr over 15 yr, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the very high-energy gamma-ray fluxes with a period of 316.7 ± 4.4 days is reported, consistent with the period of 317.3 ± 0.7 days obtained with a refined analysis of X-ray data. The analysis of data from four orbital cycles with dense observational coverage reveals short-timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over a timescale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance but cannot find any correlation of optical H$α$ parameters with fluxes at X-ray or gamma-ray energies in simultaneous observations. The key finding is that the emission of HESS J0632 + 057 in the X-ray and gamma-ray energy bands is highly variable on different timescales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system and discussed in the context of results reported on other gamma-ray binary systems.
79 ASTRONOMY AND ASTROPHYSICS↗
Spin and Accretion Rate Dependence of Black Hole X-Ray Spectra
We present a survey of how the spectral features of black hole X-ray binary systems depend on spin, accretion rate, viewing angle, and Fe abundance when predicted on the basis of first-principles physical calculations. The power-law component hardens with increasing spin. The thermal component strengthens with increasing accretion rate. The Compton bump is enhanced by higher accretion rate and lower spin. The Fe Kα equivalent width grows sublinearly with Fe abundance. Strikingly, the Kα profile is more sensitive to accretion rate than to spin because its radial surface brightness profile is relatively flat, and higher accretion rate extends the production region to smaller radii. The overall radiative efficiency is at least 30%–100% greater than as predicted by the Novikov–Thorne model.
79 ASTRONOMY AND ASTROPHYSICS↗
Hubble Space Telescope Observations of GW170817: Complete Light Curves and the Properties of the Galaxy Merger of NGC 4993
We present the complete set of Hubble Space Telescope imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 yr post-merger, we reanalyze the full light curve of AT 2017gfo across 12 bands from 5 to 1273 rest-frame days after merger. We obtain four new detections of the short γ-ray burst 170817A afterglow from 109 to 170 rest-frame days post-merger. These detections are consistent with the previously observed β = -0.6 spectral index in the afterglow light curve with no evidence for spectral evolution. We also analyze our limits in the context of kilonova afterglow or IR dust echo emission but find that our limits are not constraining for these models. We use the new data to construct deep optical and IR stacks, reaching limits of M = -6.3 to -4.6 mag, to analyze the local environment around AT 2017gfo and low surface brightness features in its host galaxy NGC 4993. We rule out the presence of any globular cluster at the position of AT 2017gfo to 2.3 × 10$^{4}$L$_{⊙}$, including those with the reddest V - H colors. Finally, we analyze the substructure of NGC 4993 in deep residual imaging and find shell features that extend up to 71."8 (14.2 kpc) from NGC 4993. The shells have a cumulative stellar mass of 6.3 × 10$^{8}$M$_{⊙}$, roughly 2% of NGC 4993, and mass-weighted ages of >3 Gyr. We conclude that it was unlikely that the GW170817 progenitor system formed in the galaxy merger.
79 ASTRONOMY AND ASTROPHYSICS↗
Angular-momentum Transport in Proto-neutron Stars and the Fate of Neutron Star Merger Remnants
Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (viscosity) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants where dΩ/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2 M TOV for forming a stable solid-body rotating NS remnant (where M TOV is the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS
79 ASTRONOMY AND ASTROPHYSICS↗
Contracting and Expanding Binary Black Holes in 3D Low-mass AGN Disks: The Importance of Separation
LIGO/Virgo has detected several binary black hole (BBH) merger events that may have originated in the accretion disks of active galactic nuclei (AGNs). These events require individual black hole masses that fall within the pair instability supernova mass gap, and therefore these black holes may have been grown from hierarchical mergers. AGN disks are a prime environment for hierarchical mergers and thus a potential location for the progenitors of BBH gravitational wave events. Understanding how a BBH embedded in an AGN disk interacts with the surrounding environment is thus crucial for determining whether this interaction can lead to its merger. However, there are few high-fidelity simulations of this process, and almost all are two-dimensional. We present the results from 3D, high-resolution, local shearing-box simulations of an embedded BBH interacting with an AGN disk. In these first simulations of their kind, we focus on determining the mass accretion rate and the orbital evolution rate at different BBH separations. We find that circular, equal-mass BBHs with separations greater than 10% of their Hill radius contract while accreting at a super-Eddington rate. At smaller separations, however, our 3D simulations find that BBHs expand their orbits. This result suggests that it may be difficult for an AGN disk to push a BBH to merger, but we discuss several mechanisms, including MHD turbulence and radiative and mechanical feedback, which could alleviate this difficulty.
79 ASTRONOMY AND ASTROPHYSICS↗
The Origin and Evolution of the Normal Type Ia SN 2018aoz with Infant-phase Reddening and Excess Emission
SN 2018aoz is a Type Ia SN with a B-band plateau and excess emission in infant-phase light curves ≲1 day after the first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN 2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show that the SN is intermediate between two subtypes of normal Type Ia: core normal and broad line. The excess emission may be attributable to the radioactive decay of surface iron-peak elements as well as the interaction of ejecta with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on Hα and He ι favor a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O ι] and He ι disfavors a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN 2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe ιι] and [Ni ιι]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1D models are incompatible with the infant-phase excess emission, $B_{max}$--$V_{max}$ color, and weak strength of nebular-phase [Ca ιι]. Although the explosion processes of SN 2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear to be normal after ~1 day.
79 ASTRONOMY AND ASTROPHYSICS↗
Deep Search for Gamma-Ray Emission from the Accreting X-Ray Pulsar 1A 0535+262
Binary systems are a well-established subclass of gamma-ray sources. The high mass X-ray binary pulsar 1A 0535+262 has been considered to be a possible gamma-ray emitter for a long time, although former gamma-ray searches using the Fermi Large Area Telescope (LAT) and VERITAS data resulted in upper limits only. We aim at a deep search for gamma-ray emission and pulsations from 1A 0535+262 using more than 13 yr of Fermi-LAT data. The analysis was performed for both the whole Fermi-LAT data set, as well as for the X-ray outbursts that 1A 0535+262 has experienced since the launch of Fermi. Various X-ray observations have been used to generate the ephemeris for the pulsation search. We also investigate the long-term gamma-ray flux variability and perform orbital phase-resolved analysis for the outbursts. We did not detect any steady or pulsed gamma-ray emission from 1A 0535+262 during the whole Fermi-LAT mission span or its X-ray outbursts. We thus derived the deepest gamma-ray luminosity upper limits to date at the 95% confidence level to be around (2.3–4.7) × 10 32 erg s -1 depending on different spectral indices assumed, which results in a ratio of L γ to L X (2–150 keV) being (1.9–3.9) × 10 -6 .
79 ASTRONOMY AND ASTROPHYSICS↗
Proposed Lunar Measurements of r-Process Radioisotopes to Distinguish the Origin of Deep-sea 244 Pu
Abstract 244 Pu has recently been discovered in deep-sea deposits spanning the past 10 Myr, a period that includes two 60 Fe pulses from nearby supernovae. 244 Pu is among the heaviest r -process products, and we consider whether it was created in supernovae, which is disfavored by nucleosynthesis simulations, or in an earlier kilonova event that seeded the nearby interstellar medium with 244 Pu that was subsequently swept up by the supernova debris. We discuss how these possibilities can be probed by measuring 244 Pu and other r -process radioisotopes such as 129 I and 182 Hf, both in lunar regolith samples returned to Earth by missions such as Chang’e and Artemis, and in deep-sea deposits.
79 ASTRONOMY AND ASTROPHYSICS↗