Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PULSARS”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nonequilibrium effects on stability of hybrid stars with first-order phase transitions

The stability of hybrid stars with first-order phase transitions as determined by calculating fundamental radial oscillation modes is known to differ from the predictions of the widely used Bardeen-Thorne-Meltzer criterion. We consider the effects of out-of-chemical-equilibrium physics on the radial modes and hence stability of these objects. For a barotropic equation of state, this is done by allowing the adiabatic sound speed to differ from the equilibrium sound speed. We show that doing so extends the stable branches of stellar models, allowing stars with rapid phase transitions to support stable higher-order stellar multiplets similarly to stars with multiple slow phase transitions. We also derive a new junction condition to impose on the oscillation modes at the phase transition. Termed the reactive condition, it is physically motivated, consistent with the generalized junction conditions between two phases, and has the common rapid and slow conditions as limiting cases. Unlike the two common cases, it can only be applied to nonbarotropic stars. Here, we apply this junction condition to hybrid stellar models generated using a two-phase equation of state consisting of nuclear matter with unpaired quark matter at high densities joined by a first-order phase transition and show that, like in the slow limiting case, stars that are classically unstable are stabilized by a finite chemical reaction speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

G objects and primordial black holes

We suggest that "G objects" recently discovered in the Galactic Center may be clouds of gas bound by the gravitational field of stellar-mass black holes produced in the interactions of sublunar primordial black holes with neutron stars. If dark matter is composed of primordial black holes with masses (10 –16 –10 –10 )M ⊙ , these black holes can be captured by neutron stars in the Galactic Center, where the dark matter density is high. After the capture, the neutron star is consumed by the black hole, resulting in a population of (1–2)M ⊙ black holes. Furthermore, these stellar-mass black holes, accompanied by gaseous atmospheres, can account for the observed properties of the G objects, including their resilience to tidal disruption by the supermassive black hole in the Galactic Center while also producing emission consistent with inferred luminosities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detectability of QCD phase transitions in binary neutron star mergers: Bayesian inference with the next generation gravitational wave detectors

We study the detectability of postmerger QCD phase transitions in neutron star binaries with next-generation gravitational-wave detectors Cosmic Explorer and Einstein Telescope. We perform numerical relativity simulations of neutron star mergers with equations of state that include a quark deconfinement phase transition through either a Gibbs or Maxwell construction. These are followed by Bayesian parameter estimation of the associated gravitational-wave signals using the nrpmw waveform model, with priors inferred from the analysis of the inspiral signal. We assess the ability of the model to measure the postmerger peak frequency $f$$^{peak}_{2}$ and identify aspects that should be improved in the model. We show that, even at postmerger signal to noise ratios as low as 10, the model can distinguish (at the 90% level) $f$$^{peak}_{2}$ between binaries with and without a phase transition in most cases. Phase-transition induced deviations in the $f$$^{peak}_{2}$ from the predictions of equation-of-state insensitive relations can also be detected if they exceed 1.6⁢σ. Our results suggest that next-generation gravitational wave detectors can measure phase transition effects in binary neutron star mergers. Furthermore, unless the phase transition is “strong,” disentangling it from other hadronic physics uncertainties will require significant theory improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗
DOE OSTI2024

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Classification of compact objects and model comparison using EOS knowledge

Nuclear theory and experiments, alongside astrophysical observations, constrain the equation of state (EOS) of supranuclear-dense matter. Conversely, knowledge of the EOS allows an improved interpretation of nuclear or astrophysical data. In this article, we use several established constraints on the EOS and the new NICER measurement of PSR J0437-4715 to comment on the nature of the primary companion in GW230529 and the companion of PSR J0514-4002E. We find that, with a probability of ≳84% and ≳68%, respectively, both objects are black holes. These likelihoods increase to above 95% when one uses GW170817’s remnant as an upper limit on the TOV mass. We also demonstrate that the current knowledge of the EOS substantially disfavors high masses and radii for PSR J⁢0030+0451, inferred recently when combining NICER with XMM-Newton background data and using particular hot-spot models. Lastly, we also use our obtained EOS knowledge to comment on measurements of the nuclear symmetry energy, finding that the large value predicted by the PREX-II measurement displays some mild tension with other constraints on the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Asymptotic-state prediction for fast flavor transformation in neutron star mergers

Neutrino flavor instabilities appear to be omnipresent in dense astrophysical environments, thus presenting a challenge to large-scale simulations of core-collapse supernovae and neutron star mergers (NSMs). Subgrid models offer a path forward, but require an accurate determination of the local outcome of such conversion phenomena. Focusing on “fast” instabilities, related to the existence of a crossing between neutrino and antineutrino angular distributions, we consider a range of analytical mixing schemes, including a new, fully three-dimensional one, and also introduce a new machine learning (ML) model. We compare the accuracy of these models with the results of several thousands of local dynamical calculations of neutrino evolution from the conditions extracted from classical NSM simulations. Our ML model shows good overall performance, but struggles to generalize to conditions from a NSM simulation not used for training. The multidimensional analytic model performs and generalizes even better, while other analytic models (which assume axisymmetric neutrino distributions) do not have reliably high performances, as they notably fail as expected to account for effects resulting from strong anisotropies. As a result, the ML and analytic subgrid models extensively tested here are both promising, with different computational requirements and sources of systematic errors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impact of dark matter on tidal signatures in neutron star mergers with the Einstein Telescope

If dark matter (DM) accumulates inside neutron stars (NS), it changes their internal structure and causes a shift of the tidal deformability from the value predicted by the dense-matter equation of state (EOS). In principle, this shift could be observable in the gravitational-wave (GW) signal of binary neutron star (BNS) mergers. We investigate the effect of fermionic, noninteracting DM when observing a large number of GW events from DM-admixed BNSs with the precision of the proposed Einstein telescope (ET). Specifically, we study the impact on the recovery of the baryonic EOS and whether DM properties can be constrained. For this purpose, we create event catalogs of BNS mock events with DM fraction up to 1%, from which we reconstruct the posterior uncertainties with the Fisher matrix approach. Using this data, we perform joint Bayesian inference on the baryonic EOS, DM particle mass, and DM particle fraction in each event. Here, our results reveal that when falsely ignoring DM effects, the EOS posterior is biased toward softer EOSs, though the offset is rather small. Further, we find that within our assumptions of our DM model and population, ET will likely not be able to test the presence of DM in BNSs, even when combining many events and adding Cosmic Explorer (CE) to the next-generation detector network. Likewise, the potential constraints on the DM particle mass will remain weak because of degeneracies with the fraction and EOS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Dark matter catalyzed baryon destruction

WIMP-type dark matter may have additional interactions that break baryon number, leading to induced nucleon decays which are subject to direct experimental constraints from proton decay experiments. In this work, we analyze the possibility of continuous baryon destruction, deriving strong limits from the dark matter accumulating inside old neutron stars, as such a process leads to excess heat generation. We construct the simplest particle dark matter model that breaks the baryon and lepton numbers separately but conserves B - L. Virtual exchange by DM particles in this model results in dinucleon decay via nn → $n\overline{v}$ and np → ne⁺ processes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Josephson currents in neutron stars

We demonstrate that the interface between S-wave and P-wave paired superfluids in neutron stars induces a neutron supercurrent, a charge-neutral analog of the Josephson junction effect in electronic superconductors. The proton supercurrent entrainment by the neutron superfluid generates, in addition to the neutral supercurrent, a charged current across the interface. Beyond this stationary effect, the motion of the neutron vortex lines responding to secular changes in the neutron star’s rotation rate induces a time-dependent oscillating Josephson current across this interface when proton flux tubes are dragged along with them. We show that such motion produces radiation from the interface once clusters of proton flux tubes intersect the interface. The power of radiation exceeds by orders of magnitude the Ohmic dissipation of currents in neutron stars. As a result, this effect appears to be phenomenologically significant enough to heat the star and alter its cooling rate during the photon cooling era.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

How fast can protons decay?

Current laboratory bounds imply that protons are extremely long-lived. However, this conclusion may not hold for all time and in all of space. We find that the proton lifetime can be ∼ 15 orders of magnitude shorter in the relatively recent past on Earth or at the present time elsewhere in the Milky Way. A number of terrestrial and astrophysical constraints are examined and potential signals are outlined. We also sketch possible models that could lead to spatial or temporal variations in the proton lifetime. A positive signal could be compelling evidence for a new long range force of nature, with important implications for the limitations of fundamental inferences based solely on laboratory measurements. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Cocoon shock breakout emission from binary neutron star mergers

Shock breakout emission is among the first observable signals in a wide variety of astrophysical phenomena, including neutron star (NS) mergers, and it can be the dominant component in low-luminosity short γ-ray bursts, as exemplified by GRB 170817A. In this work, we investigate the cocoon shock breakout emission in NS mergers and how its signal depends on the outermost layers of the ejecta profile, which we derive from general relativistic radiation hydrodynamic simulations. We study the formation of the cocoon as a consequence of a relativistic jet propagating through the ejecta. To explore the influence of the outermost layers of the ejecta on the breakout emission, we explore cases where the ejecta has a sharp cutoff or an extended smooth tail. We find that the shock breakout emission is strongly influenced by the shape of the ejecta outer layers, with extended tails yielding results consistent with the observed properties of GRB 170817A, whereas sharp cutoffs overestimate the radiated energy. Using a Bayesian analysis, we estimate the best fit parameters of the central engine, considering both accreting black hole and magnetized neutron star scenarios. Our findings indicate a slight preference for the scenarios where the engine is a black hole. Furthermore, our work probes the nature of neutron star mergers and highlights the importance of the shape of the ejecta profile in modeling early electromagnetic counterparts to these mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Open-source library for performance-portable neutrino reaction rates: Application to neutron star mergers

A realistic and detailed description of neutrinos in binary neutron star (BNS) mergers is essential to build reliable models of such systems. To this end, we present bns_nurates, a novel open-source numerical library designed for the efficient on-the-fly computation of neutrino interactions, with particular focus on regimes relevant to BNS mergers. bns_nurates targets a higher level of accuracy and realism in the implementation of commonly employed reactions by accounting for relevant microphysics effects on the interactions, such as weak magnetism and mean field effects. It also includes the contributions of inelastic neutrino scattering off electrons and positrons and (inverse) nucleon decays. Finally, it offers a way to reconstruct the neutrino distribution function in the framework of moment-based transport schemes. As a first application, we compute both energy-dependent and energy-integrated neutrino emissivities and opacities for conditions extracted from a BNS merger simulation with m1 transport scheme. We find some qualitative differences in the results when considering the impact of the additional relevant reactions and of microphysics effects. For example, neutrino-electron/positron scattering reactions are important for the energy exchange of heavy-type neutrinos as they do not undergo semileptonic charged-current processes, when μ± are not accounted for. Moreover, weak magnetism and mean field effects can significantly modify the contribution of β processes for electron-type (anti)neutrinos, increasing at the same time the importance of (inverse) neutron decays. Here, the improved treatment for the reaction rates also modifies the conditions at which neutrinos decouple from matter in the system, potentially affecting their emission spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Signature of hadron-quark crossover in binary-neutron-star mergers

Here, we study observational signatures of the hadron-quark crossover in binary-neutron-star mergers by numerical-relativity simulations with various mass configurations. We employ two equations of state (EoSs) for matter consistent with inference from the observational data. In the crossover scenario the EoS is softened in a density realized in binary-neutron-star mergers and is smoothly continued to quark matter. In the phase transition scenario without crossover, the EoS remains stiff and a first-order phase transition takes place in a density out of reach of mergers. A GW170817-like system forms a remnant massive neutron star in both scenarios, and it collapses into a black hole only in the crossover scenario due to the softening while gravitational-wave emission is strong. This difference is clearly reflected in the sudden shutdown of gravitational waves. For a given EoS, the lifetime of the merger remnant is determined primarily by the total mass of the system. Identifying these features in a variety of future events with the next generation of ground-based gravitational-wave detectors will enable us to clarify details of hadron-quark transition. The mass of the accretion disk surrounding the remnant black hole is affected not only by the lifetime of the remnant but also by the mass ratio of the system. Electromagnetic emission associated with the disk outflow will also be useful for detailed investigation of the hadron-quark transition.

QCD phase transitions↗