Engineering Papers⌕ Search

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.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Exploring the connection between compact object mergers and fast X-ray transients

Context. The connection between compact object mergers and some extragalactic fast X-ray transients (FXRTs) has long been hypothesised but never ultimately established. Aims. In this work, we investigate two FXRTs, the LEIA X-ray Transient LXT 240402A and the Einstein Probe EP 250207b, whose precise positions lie close to nearby (z ≲ 0.1) quiescent galaxies with a negligible probability of chance coincidence, identifying them as particularly promising cases of merger-driven explosions in the local Universe. Methods. We used Chandra to derive accurate localisations for both events and secure otherwise ambiguous associations with their optical counterparts. Deep optical and near-infrared observations with VLT, GTC, and LBT were performed to characterise the surrounding environment and search for kilonova emission, the hallmark of neutron star mergers. Complementary early-time X-ray monitoring with Swift and Einstein Probe was used to constrain the non-thermal afterglow. Results. We find that both FXRTs remain compatible with a compact binary merger progenitor, which produced low-mass ejecta and kilonova emission subdominant to the afterglow. However, alternative explanations such as a distant (z ≳ 1) core-collapse supernova cannot be conclusively ruled out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Horizon radiation reaction forces

Using Effective Field Theory (EFT) methods, we compute the effects of horizon dissipation on the gravitational interactions of relativistic binary black hole systems. We assume that the dynamics is perturbative, i.e it admits an expansion in powers of Newton’s constant (post-Minkowskian, or PM, approximation). As applications, we compute corrections to the scattering angle in a black hole collision due to dissipative effects to leading PM order, as well as the post-Newtonian (PN) corrections to the equations of motion of binary black holes in non-relativistic orbits, which represents the leading order finite size effect in the equations of motion. The methods developed here are also applicable to the case of more general compact objects, eg. neutron stars, where the magnitude of the dissipative effects depends on non-gravitational physics (e.g, the equation of state for nuclear matter).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Periodic Gamma-Ray Modulation of the Blazar PG 1553+113 Confirmed by Fermi-LAT and Multiwavelength Observations

A 2.1 yr periodic oscillation of the gamma-ray flux from the blazar PG 1553+113 has previously been tentatively identified in ∼7 yr of data from the Fermi Large Area Telescope. After 15 yr of Fermi sky-survey observations, doubling the total time range, we report >7 cycle gamma-ray modulation with an estimated significance of 4σ against stochastic red noise. Independent determinations of oscillation period and phase in the earlier and the new data are in close agreement (chance probability <0.01). Pulse timing over the full light curve is also consistent with a coherent periodicity. Multiwavelength new data from Swift X-Ray Telescope, Burst Alert Telescope, and UVOT, and from KAIT, Catalina Sky Survey, All-Sky Automated Survey for Supernovae, and Owens Valley Radio Observatory ground-based observatories as well as archival Rossi X-Ray Timing Explorer satellite-All Sky Monitor data, published optical data of Tuorla, and optical historical Harvard plates data are included in our work. Optical and radio light curves show clear correlations with the gamma-ray modulation, possibly with a nonconstant time lag for the radio flux. We interpret the gamma-ray periodicity as possibly arising from a pulsational accretion flow in a sub-parsec binary supermassive black hole system of elevated mass ratio, with orbital modulation of the supplied material and energy in the jet. Other astrophysical scenarios introduced include instabilities, disk and jet precession, rotation or nutation, and perturbations by massive stars or intermediate-mass black holes in polar orbit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Is GW170817 a Multimessenger Neutron Star-Primordial Black Hole Merger?

We investigate the possibility of the gravitational-wave event GW170817 being a light, solar-mass black hole (BH) - neutron star (NS) merger. We explore two exotic scenarios involving primordial black holes (PBH) that could produce such an event, taking into account available observational information on NGC 4993. First, we entertain the possibility of dynamical NS-PBH binary formation where a solar-mass PBH and a NS form a binary through gravitational interaction. We find that while dynamical NS-PBH formation could account for the GW170817 event, the rate is highly dependent on unknown density contrast factors and could potentially be affected by galaxy mergers. We also find that PBH-PBH binaries would likely have a larger merger rate, assuming the density contrast boost factor of an order similar to the NS-PBH case. These exotic merger formations could provide new channels to account for the volumetric rate of compact-object mergers reported by LIGO/Virgo. Secondly, we consider the case where one of the NS's in a binary NS system is imploded by a microscopic PBH. We find that the predicted rate for NS implosion into a BH is very small, at least for the specific environment of NGC 4993. We point out that similar existing (e.g. GW190425 and GW190814) and future observations will shed additional light on these scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A helium-burning white dwarf binary as a supersoft X-ray source

Abstract Type Ia supernovae are cosmic distance indicators 1,2 , and the main source of iron in the Universe 3,4 , but their formation paths are still debated. Several dozen supersoft X-ray sources, in which a white dwarf accretes hydrogen-rich matter from a non-degenerate donor star, have been observed 5 and suggested as Type Ia supernovae progenitors 6–9 . However, observational evidence for hydrogen, which is expected to be stripped off the donor star during the supernova explosion 10 , is lacking. Helium-accreting white dwarfs, which would circumvent this problem, have been predicted for more than 30 years (refs. 7,11,12 ), including their appearance as supersoft X-ray sources, but have so far escaped detection. Here we report a supersoft X-ray source with an accretion disk whose optical spectrum is completely dominated by helium, suggesting that the donor star is hydrogen-free. We interpret the luminous and supersoft X-rays as resulting from helium burning near the surface of the accreting white dwarf. The properties of our system provide evidence for extended pathways towards Chandrasekhar-mass explosions based on helium accretion, in particular for stable burning in white dwarfs at lower accretion rates than expected so far. This may allow us to recover the population of the sub-energetic so-called Type Iax supernovae, up to 30% of all Type Ia supernovae 13 , within this scenario.

Science & Technology - Other Topics↗
DOE OSTI2020

GW190814: Gravitational Waves from the Coalescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object

We report the observation of a compact binary coalescence involving a 22.2-24.3 M ⊙ black hole and a compact object with a mass of 2.50-2.67 M ⊙ (all measurements quoted at the 90% credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on 2019 August 14 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg 2 at a distance of 241 +41-41 Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves , 0.112+0.0090.008, , and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to ≤0.07. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1-23 Gpc -3 yr -1 for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models of the formation and mass distribution of compact-object binaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

First Resolved Scattered-light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

In this work, we present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light, and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet–disk system. Uniquely, HD 98363 has a wide comoving stellar companion, Wray 15-788, with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one nondetection to one of the brightest disks in the region. With the new results, a total of 15 A/F stars in Sco-Cen have resolved scattered-light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F stars at this important stage in the development of planetary systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A prescription for the asteroseismic surface correction

In asteroseismology, the surface effect refers to a disparity between the observed and the modelled frequencies in stars with solar-like oscillations. It originates from improper modelling of the surface layers. Correcting the surface effect usually requires using functions with free parameters, which are conventionally fitted to the observed frequencies. On the basis that the correction should vary smoothly across the H–R diagram, we parameterize it as a simple function of surface gravity, effective temperature, and metallicity. We determine this function by fitting a wide range of stars. The absolute amount of the surface correction decreases with luminosity, but the ratio between it and νmax increases, suggesting the surface effect is more important for red giants than dwarfs. Applying the prescription can eliminate unrealistic surface correction, which improves parameter estimations with stellar modelling. Using two open clusters, we found a reduction of scatter in the model-derived ages for each star in the same cluster. As an important application, we provide a new revision for the Δν scaling relation that, for the first time, accounts for the surface correction. The values of the correction factor, f Δν , are up to 2 percent smaller than those determined without the surface effect considered, suggesting decreases of up to 4 percent in radii and up to 8 percent in masses when using the asteroseismic scaling relations. This revision brings the asteroseismic properties into an agreement with those determined from eclipsing binaries. The new correction factor and the stellar models with the corrected frequencies are available at https://www.github.com/parallelpro/surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Discovery and characterization of five new eclipsing AM CVn systems

ABSTRACT AM CVn systems are ultra-compact, hydrogen-depleted, and helium-rich, accreting binaries with degenerate or semidegenerate donors. We report the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.5, 53.3, 37.4, and 35.4 min. These systems were discovered by searching for deep eclipses in the Zwicky Transient Facility (ZTF) light curves of white dwarfs selected using Gaia parallaxes. We obtained phase-resolved spectroscopy to confirm that all systems are AM CVn binaries, and we obtained high-speed photometry to confirm the eclipse and characterize the systems. The spectra show double-peaked H e lines but also show metals, including K and Zn, elements that have never been detected in AM CVn systems before. By modelling the high-speed photometry, we measured the mass and radius of the donor star, potentially constraining the evolutionary channel that formed these AM CVn systems. We determined that the average mass of the accreting white dwarf is ≈0.8 M⊙, and that the white dwarfs in long-period systems are hotter than predicted by recently updated theoretical models. The donors have a high entropy and are a factor of ≈2 more massive compared to zero-entropy donors at the same orbital period. The large donor radius is most consistent with H e-star progenitors, although the observed spectral features seem to contradict this. The discovery of five new eclipsing AM CVn systems is consistent with the known observed AM CVn space density and estimated ZTF recovery efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Tidally Excited Modes and δ Scuti Pulsations in the Eclipsing Triple Star IM Persei

IM Per is a triple star system whose eclipsing pair masses and radii are accurate to within 1%. The Transiting Exoplanet Survey Satellite (TESS) light curve of the program target exhibits partial eclipses and multiple oscillations with millimagnitude-level amplitudes. It is found that the oscillations affect eclipse timing measurements. Binary modeling of the high-quality TESS data indicates that the eclipsing components of the triple system are twin dwarfs with parameters of M {sub 2}/M {sub 1} = 0.995, R {sub 2}/R {sub 1} = 0.901, and Δ(T{sub eff,1}–T{sub eff,2}) = 12 K in an eccentric (e = 0.049), detached configuration. The third light of l {sub 3} = 0.054 may mostly come from a G-type tertiary companion. Our predicted parallax of 1.52 ± 0.09 mas is concurrent with the Gaia measurement of 1.52 ± 0.05 mas. Multifrequency analysis of the outside-eclipse residuals reveals 22 significant pulsation signals: 4 in the gravity-mode region (0.03−2.22 day{sup −1}) and 18 in the pressure-mode region (9.19−25.12 day{sup −1}). Of the low frequencies, f {sub 11} and f {sub 14} are orbital harmonics that can be identified as tidally excited modes. The pulsation periods and constants for the high frequencies, and the position in the Cepheid instability strip demonstrate that the pulsating component of IM Per is a δ Sct variable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Four new deeply eclipsing white dwarfs in Zwicky Transient Facility

ABSTRACT We present the results of a search for deeply eclipsing white dwarfs in the Zwicky Transient Facility (ZTF) Data Release 4 (DR4). We identify nine deeply eclipsing white dwarf candidates, four of which we followed up with high-cadence photometry and spectroscopy. Three of these systems show total eclipses in the ZTF data and our follow-up Apache Point Observatory 3.5 m telescope observations. Even though the eclipse duration is consistent with sub-stellar companions, our analysis shows that all four systems contain a white dwarf with low-mass stellar companions of ∼0.1 M⊙. We provide mass and radius constraints for both stars in each system based on our photometric and spectroscopic fitting. Finally, we present a list of 41 additional eclipsing WD+M candidates identified in a preliminary search of ZTF DR7, including 12 previously studied systems. We identify two new candidate short-period, eclipsing, white dwarf–brown dwarf binaries within our sample of 41 WD+M candidates based on Pan-STARRS colours.

Kosakowski, A. (ORCID:0000000298781647)↗
DOE OSTI2021

Total X-Ray Emission from the LMC Observed with HaloSat

We constrain the global parameters of the Large Magellanic Cloud (LMC) by modeling the first soft X-ray (0.4–6.8 keV) observations of the entire LMC produced in a single pointing with moderate spectral resolution by the HaloSat CubeSat. These data are best fit with the sum of two thermal plasma components in collisional ionization equilibrium and a power law. We find cool (0.210 ± 0.014 keV) and hot (0.89 ± 0.14 keV) components. The total X-ray luminosity of the LMC is (1.08 ± 0.14) × 10{sup 39} erg s{sup −1}. X-ray binaries make up a large fraction of the emission with a luminosity of (6.0 ± 0.8) × 10{sup 38} erg s{sup −1}, followed by cool gas from superbubbles, supernovae, and diffuse emission with a luminosity of (3.0 ± 0.3) × 10{sup 38} erg s{sup −1}. The hot gas from star formation contributes the smallest fraction, with a luminosity of (1.9 ± 0.5) × 10{sup 38} erg s{sup −1}. We estimate the total volume of the cool gas to be between (0.2–1.2) × 10{sup 10} pc{sup 3} and the hot gas to be between (1.0–5) × 10{sup 7} pc{sup 3} for filling factors of f = 1 and 0.2. These volumes result in a total thermal energy for the cool gas between (1.4–3) × 10{sup 54} erg for electron densities of 0.017–0.04 cm{sup −3}, and a thermal energy for the hot gas between (1.7–4) × 10{sup 53} erg for electron densities of 0.12–0.3 cm{sup −3}. This yields cooling timescales for the cool and hot gas of (1.5–3) × 10{sup 8} yr and (1.8–6) × 10{sup 7} yr, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

An Isolated Stellar-mass Black Hole Detected through Astrometric Microlensing

We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t E ≃ 270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of 6 yr, reveals a clear relativistic astrometric deflection of the background star's apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of Earth's motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of 7.1 ± 1.3 M ⊙ and a distance of 1.58 ± 0.18 kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion for the BH. The proper motion is offset from the mean motion of Galactic disk stars at similar distances by an amount corresponding to a transverse space velocity of ~45 km s -1 , suggesting that the BH received a "natal kick" from its supernova explosion. Previous mass determinations for stellar-mass BHs have come from radial velocity measurements of Galactic X-ray binaries and from gravitational radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is the first for an isolated stellar-mass BH using any technique.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Numerical relativity simulations of prompt collapse mergers: Threshold mass and phenomenological constraints on neutron star properties after GW170817

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite-temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, binary tidal parameter at the threshold (Λth), maximum mass of nonrotating NSs, and radii of reference mass NSs. We combine the EOS-insensitive relations, phenomenological constraints on NS properties, and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and a 1.6 M ⊙ NS of 9.81 and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability (Λ) of a 1.4 M ⊙ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters—fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We discuss an approach to constrain the upper as well as lower limit of NS maximum mass using future binary NS detections and their identification as prompt or delayed collapse. With future observations, it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Igniting Weak Interactions in Neutron Star Postmerger Accretion Disks

The merger of two neutron stars or a neutron star and a black hole typically results in the formation of a postmerger accretion disk. Outflows from disks may dominate the overall ejecta from mergers and be a major source of r-process nuclei in our universe. We explore the parameter space of such disks and their outflows and r-process yields by performing 3D general-relativistic magnetohydrodynamic simulations with weak interactions and approximate neutrino transport. We discuss the mapping between the initial binary parameters and the parameter space of the resulting disks, chiefly characterized by their initial accretion rate. We demonstrate the existence of an ignition threshold for weak interactions at around ~10 -3 M⊙ s -1 for typical parameters by means of analytic calculations and numerical simulations. We find a degenerate, self-regulated, neutrino-cooled regime above the threshold and an advection-dominated regime below the threshold. Excess heating in the absence of neutrino cooling below the threshold leads to ≳60% of the initial disk mass being ejected in outflows, with typical velocities of ~(0.1–0.2)c, compared to ≲40% at ~(0.1–0.15)c above the threshold. While disks below the threshold show suppressed production of light r-process elements, disks above the threshold can produce the entire range of r-process elements, in good agreement with the observed solar system abundances. Disks below the ignition threshold may produce an overabundance of actinides seen in actinide-boost stars. As gravitational-wave detectors start to sample the neutron star merger parameter space, different disk realizations may be observable via their associated kilonova emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum field theory, worldline theory, and spin magnitude change in orbital evolution

A previous paper [Z. Bern et al., Binary dynamics through the fifth power of spin at 𝑂⁡(𝐺 2 ), Phys. Rev. Lett. 130, 201402 (2023)] identified a puzzle stemming from the amplitudes-based approach to spinning bodies in general relativity: additional Wilson coefficients appear compared to current worldline approaches to conservative dynamics of generic astrophysical objects, including neutron stars. In this paper we clarify the nature of analogous Wilson coefficients in the simpler theory of electrodynamics. We analyze the original field-theory construction, identifying definite-spin states some of which have negative norms, and relating the additional Wilson coefficients in the classical theory to transitions between different quantum spin states. We produce a new version of the theory which also has additional Wilson coefficients, but no negative-norm states. We match, through 𝒪⁡(𝛼 2 ) and 𝒪⁡(𝑆 2 ), the Compton amplitudes of these field theories with those of a modified worldline theory with extra degrees of freedom introduced by releasing the spin supplementary condition. We build an effective two-body Hamiltonian that matches the impulse and spin kick of the modified field theory and of the worldline theory, displaying additional Wilson coefficients compared to standard worldline approaches. The results are then compactly expressed in terms of an eikonal formula. Our key conclusion is that, contrary to standard approaches, while the magnitude of the spin tensor is still conserved, the magnitude of the spin vector can change under conserved Hamiltonian dynamics and this change is governed by the additional Wilson coefficients. For specific values of Wilson coefficients the results are equivalent to those from a definite spin obeying the spin supplementary condition, but for generic values they are physically inequivalent. These results warrant detailed studies of the corresponding issues in general relativity.

classical black holes↗
DOE OSTI2022

Radiation hydrodynamics modelling of kilonovae with SNEC

ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

79 ASTRONOMY AND ASTROPHYSICS↗