Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relativistic 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 55 records · Page 3

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

79 ASTRONOMY AND ASTROPHYSICS↗

X-Ray Spectra from General Relativistic Radiation Magnetohydrodynamic Simulations of Thin Disks

We compare X-ray emission from several general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin black hole accretion disks with different accretion rates and spins. The simulations were performed using the M1 closure scheme, resolved with 12 frequency (energy) bins logarithmically spaced from 5 × 10 −3 to 5 × 10 3 ​​​ keV. We apply a general relativistic Monte Carlo transport code to postprocess the simulation data with greater fidelity in frequency resolution and Compton scattering treatment. Despite the relatively few energy bins and Kompaneets approximation to Compton scattering utilized in the M1 method, we find generally good agreement between the methods. Both produce prominent thermal profiles with peaks around 2–2.5 keV, where agreement is particularly strong and representative of the soft state. Both also find weaker (lower-luminosity) thermally sourced emission extending out to 100 keV due to the hotter innermost regions of the disks. Inverse Compton scattering becomes increasingly effective at hardening spectral outputs with increasing black hole spin, and becomes the dominant mechanism for photons that escape with energies between 10 to several hundred keV. At very high rates of spin, the radiation flux in this upscattered component becomes comparable to the thermal flux, a phenomenon typically associated with intermediate states. Beyond 10 4 keV, we observe faint, free–free emission from hot, optically thin coronal regions developing near the horizon, common to both spinning and nonspinning black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

Large-scale Evolution of Seconds-long Relativistic Jets from Black Hole–Neutron Star Mergers

We present the first numerical simulations that track the evolution of a black hole–neutron star (BH–NS) merger from premerger to r ≳ 10 11 cm. The disk that forms after a merger of mass ratio q = 2 ejects massive disk winds (3–5 × 10 –2 M ⊙ ). We introduce various postmerger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off as L j ~ t –2 . All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For postmerger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetization σ0 > 100 retain significant magnetization (σ $\gg$ 1) at r > 10 10 cm, emphasizing the importance of magnetic processes in the prompt emission. The jet–wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗

Turbulence modelling in neutron star merger simulations

Observations of neutron star mergers have the potential to unveil detailed physics of matter and gravity in regimes inaccessible by other experiments. Quantitative comparisons to theory and parameter estimation require nonlinear numerical simulations. However, the detailed physics of energy and momentum transfer between different scales, and the formation and interaction of small scale structures, which can be probed by detectors, are not captured by current simulations. This is where turbulence enters neutron star modelling. This review will outline the theory and current status of turbulence modelling for relativistic neutron star merger simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetized Outflows from Short-lived Neutron Star Merger Remnants Can Produce a Blue Kilonova

Abstract We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino–matter interactions and is well suited to studying the resulting nucleosynthesis and kilonova emission. A magnetized wind is driven from the remnant and ejects neutron-rich material at a quasi-steady-state rate of 0.8 × 10 −1 M ⊙ s −1 . We find that the ejecta in our simulations underproduce r -process abundances beyond the second r -process peak. For sufficiently long-lived remnants, these outflows alone can produce blue kilonovae, including the blue kilonova component observed for AT2017gfo.

Astronomy & Astrophysics↗

The Role of Strong Gravity and the Nuclear Equation of State on Neutron-star Common-envelope Accretion

Common-envelope evolution is important in the formation of neutron star binaries within the isolated binary formation channel. As a neutron star inspirals within the envelope of a primary massive star, it accretes and spins up. Because neutron stars are in the strong-gravity regime, they have a substantial relativistic mass deficit, i.e., their gravitational mass is less than their baryonic mass. This effect causes some fraction of the accreted baryonic mass to convert into neutron star binding energy. The relativistic mass deficit also depends on the nuclear equation of state, since more compact neutron stars will have larger binding energies. We model the mass growth and spin-up of neutron stars inspiraling within common-envelope environments and quantify how different initial binary conditions and hadronic equations of state affect the post-common-envelope neutron star's mass and spin. From these models, we find that neutron star mass growth is suppressed by ≈15%–30%. We also find that for a given amount of accreted baryonic mass, more compact neutron stars will spin-up faster while gaining less gravitational mass, and vice versa. Here, this work demonstrates that a neutron star's strong gravity and nuclear microphysics plays a role in neutron-star-common-envelope evolution, in addition to the macroscopic astrophysics of the envelope. Strong gravity and the nuclear equation of state may thus affect both the population properties of neutron star binaries and the cosmic double neutron star merger rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Detection of X-ray emission from the PSR 1259-63/SS 2883 binary system

Nonpulsed but variable X-ray emission has been detected from the binary system containing the radio pulsar PSR 1259-63 during two pointed ROSAT observations, taken 5 months apart. This 47.7 ms radio pulsar is in a highly eccentric (epsilon approximately 0.85) binary system with the 10-15 solar mass Be star SS 2883. It is the first radio pulsar found to be in a binary system with a massive main-sequence companion; it is also the most highly eccentric binary system known to contain a neutron star. The level of X-ray flux detected in the ROSAT observations has increased with orbital phase by a factor of at least 10 between 1992 February and 1993 February. The X-ray flux is significantly greater than expected from the Be star's corona and seems likely to originate either from low-level stellar wind accretion onto the neutron star or from the shock between the stellar wind and the relativistic pulsar wind. The system may be the progenitor of the more slowly rotating Be X-ray binary pulsar systems.

Cominsky, Lynn↗

LISA Mission Concept Study, Laser Interferometer Space Antenna for the Detection and Observation of Gravitational Waves

This document presents the results of a design feasibility study for LISA (Laser Interferometer Space Antenna). The goal of LISA is to detect and study low-frequency astrophysical gravitational radiation from strongly relativistic regions. Astrophysical sources potentially visible to LISA include extra-galactic massive black hole binaries at cosmological distances, binary systems composed of a compact star and a massive black hole, galactic neutron star-black hole binaries, and background radiation from the Big Bang. The LISA mission will comprise three spacecraft located five million kilometers apart forming an equilateral triangle in an Earth-trailing orbit. Fluctuations in separation between shielded test masses located within each spacecraft will be determined by optical interferometry which determines the phase shift of laser light transmitted between the test masses.

Folkner, W. M.↗

GR-Athena++: Puncture Evolutions on Vertex-centered Oct-tree Adaptive Mesh Refinement

Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ~1.2 × 10 4 CPUs and excellent weak scaling is shown up to ~10 5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.

79 ASTRONOMY AND ASTROPHYSICS↗

Acceleration by pulsar winds in binary systems

In the absence of accretion torques, a pulsar in a binary system will spin down due to electromagnetic dipole radiation, and the spin-down power will drive a wind of relativistic electron-position pairs. Winds from pulsars with short periods will prevent any subsequent accretion but may be confined by the companion star atmosphere, wind, or magnetosphere to form a standing shock. This paper investigates the possibility of particle acceleration at such a pulsar wind shock and the production of VHE and UHE gamma-rays from interactions of accelerated protons in the companion star's wind or atmosphere. It is found that, in close binaries containing active pulsars, protons will be shock accelerated to a maximum energy dependent on the pulsar spin-down luminosity. If a significant fraction of the spin-down power goes into particle acceleration, these systems should be sources of VHE and possibly UHE gamma-rays. The application of the pulsar wind model to binary sources such as Cygnus X-3 is discussed, as well as the possibility of observing VHE gamma-rays from known binary radio pulsar systems.

Harding, Alice K.↗

Acceleration by pulsar winds in binary systems

In the absence of accretion torques, a pulsar in a binary system will spin down due to electromagnetic dipole radiation and the spin-down power will drive a wind of relativistic electron-positron pairs. Winds from pulsars with short periods will prevent any subsequent accretion but may be confined by the companion star atmosphere, wind, or magnetosphere to form a standing shock. The authors investigate the possibility of particle acceleration at such a pulsar wind shock and the production of very high energy (VHE) and ultra high energy (UHE) gamma rays from interactions of accelerated protons in the companion star's wind or atmosphere. They find that in close binaries containing active pulsars, protons will be shock accelerated to a maximum energy dependent on the pulsar spin-down luminosity. If a significant fraction of the spin-down power goes into particle acceleration, these systems should be sources of VHE and possibly UHE gamma rays. The authors discuss the application of the pulsar wind model to binary sources such as Cygnus X-3, as well as the possibility of observing VHE gamma-rays from known binary radio pulsar systems.

Harding, Alice K.↗

Relativistic Shapiro Delay Measurements of an Extremely Massive Millisecond Pulsar

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing—the process of accounting for every rotation of a pulsar over long time periods—can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system1. One of these, called relativistic Shapiro delay2, can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-yr data set with recent orbitalphase- specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be 2.14+0:10-0:09 M⊙ (68.3% credibility interval; the 95.4% credibility interval is 2.14+0:20-0:18 M⊙). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

H T Cromartie↗

Identification of the nebula G70.7 + 1.2 as a bow shock powered by a pulsar/Be-star binary

New data are presented here that firmly establish the nonthermal nature of the radio emission from the enigmatic radio and optical nebula G70.7 + 1.2. H-alpha and forbidden O I Fabry-Perot observations are used to argue that the extended optical emission from the nebula arises from a bow shock powered by a mass-losing luminous star moving supersonically through dense gas. The strong nonthermal radio emission from the object is then explained as the shocked relativistic wind from a pulsar, which is proposed here as a companion to the Be star. The coincidence of the optical and radio emission requires the pulsar and stellar winds to be mixed together. The system has a large overall velocity of about 60 km/s which is inexplicable in all other models but which is typical of binary pulsars.

Kulkarni, S. R.↗

Radiation dynamics in X-ray binaries. I - Type 1 bursts. II - Type 2 bursts. III - Extremely compact objects

Equations describing the evolution of a thin, axisymmetric, viscous, relativistic, irradiated accretion disk are presented, as well as numerical solutions of these equations in the case where irradiation results from a thermonuclear flash on the surface of the accreting neutron star. These calculations verify the notion that the radiation torque induces a substantial increase in the accretion rate, during a type 1 X-ray burst, and provide insight into the factors which influence the dynamical response of the disk. A new model for the source XBT 1730-335, the rapid burster, is presented. Temporal and spectral properties are calculated. The rapid burster is found to be a nonmagnetic, 'critically compact', slowly rotating neutron star in a highly eccentric binary system with a period of 6 mo. The spectral modifications which arise from the scattering of photons by accretion disks around nonmagnetic neutron stars are calculated. The 'black hole candidates' are interpreted as extremely compact stars.

Walker, Mark A.↗

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗

Oscillations of red dwarfs in evolved low-mass binaries with neutron stars

We investigate a novel aspect of a problem related to the properties of low-mass binaries (LMBs) with millisecond pulsars: the pulsations of the red dwarf (donor) companion of the neutron star (NS). The illumination of the donor star by the pulsar's high-energy nonthermal radiation and relativistic wind may substantially affect its structure. We present a quantitative analysis of the oscillation spectrum of a red dwarf which has evolved in an LMB and has undergone the stage of evaporation. We calculate the p- and g-modes for red dwarfs with masses in the interval (0.2-0.6) stellar mass. For comparison, similar calculations are presented for zero age main-sequence (ZAMS) stars of the same masses. For less massive donor stars (approximately 0.2 stellar mass) the oscillation spectrum becomes quantitatively different from that of their ZAMS counterparts. The differnce is due to the fact that a ZAMS star of 0.2 stellar mass is fully convective, while the donor star in an LMB is expected to be far from thermal equilibrium and not fully convective. As a result, in contrast to a low-mass ZAMS star, a red dwarf of the same mass in an LMB allows the existence of g-modes. We also consider tidally forced g-modes, and perform a linear analysis of these oscillations for different degrees of nonsynchronism between the orbital and spin rotation of the red dwarf component. We demonstrate the existence of a series of reasonances for the low-order g-modes which may occur in LMBs at a late stage of their evolution. We discuss the possibility that these oscillations may trigger Roche lobe overflow and sudden mass loss by the donor star. Further implications of this effect for gamma- and X-ray burst phenomena are outlined.

Sarna, Marek J.↗

Effects of Landau quantization on neutrino emission and absorption

Some neutron stars known as magnetars possess very strong magnetic fields, with surface fields as large as 10 15 G and internal fields that are possibly stronger. Recent observations of the radio pulsar GLEAM-X J1627 suggest it may have a surface field as strong as 10 16 G. In the presence of a strong magnetic field, Leinson and Perez showed and Baiko and Yakovlev confirmed that the energy levels of electrons and protons are quantized and the direct Urca process allows neutron stars to cool rapidly, even at low density. For the case of magnetic fields 𝐵 ≳ 10 16 G, we find features in the emissivity due to energy quantization that are not captured by the frequently employed quasiclassical approximation where energy levels are treated as nearly continuous. Resonances can result in amplification of the neutrino emissivity at specific densities compared to a calculation that neglects quantization, particularly at low temperature. These effects are not important for the thermal evolution of an entire neutron star, but may be relevant for phenomena that depend on behavior at specific densities. Here, we present a fully relativistic calculation of the direct Urca rate in a strong magnetic field using the standard V-A weak Lagrangian incorporating mean field nuclear effects and discuss approaches to the numerical challenge the modified wave functions present and a new semianalytic approximation. These tools are also applicable to calculating neutrino opacities in strong magnetic fields in the ejecta of binary neutron star mergers. We calculate the opacities for neutrinos capturing on free nucleons at subsaturation densities and temperatures exceeding an MeV. We find an enhancement to capture processes of the lowest energy neutrinos by an order of magnitude or more due to suppression of electron Pauli blocking in the case of capture on neutrons, and from the effect of the nucleon magnetic moments in the case of capture on protons.

astrophysical electromagnetic fields↗