Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NEUTRON 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Neutron-Antineutron Transition Search at DUNE

Neutron-antineutron transition is a baryon number violating process with ΔB=2, providing a unique insight into potential explanations of the baryon asymmetry in our universe, especially in the context of post-sphaleron baryogenesis. Studies have been conducted across various neutron-rich environments, including free neutron sources, neutron stars, and bound neutrons in large underground neutrino detectors, among which large neutrino experiments have historically obtained the most stringent constraints. The forthcoming Deep Underground Neutrino Experiment (DUNE) will offer especially strong prospects for this search, enabled by the high spatial resolution of its liquid argon time projection chamber. This poster presents a recent study on neutron-antineutron transition at DUNE using machine-learning based event classification.

Yu, Kaiwen [UC, Irvine (main)]↗
DOE OSTI2024

The Evolution of Inclined Binary Black Holes in the Disks of Active Galactic Nuclei

The accretion disks that fuel active galactic nuclei (AGNs) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational-wave events involving black holes in the pair-instability mass gap, and may facilitate electromagnetic counterparts to black hole binary mergers. However, many orders of magnitude separate the typical length scales of binary formation and those on which gravitational waves can drive binary inspirals, making binary mergers inside the disk uncertain. Previous hydrodynamical simulations of binaries have either been restricted to two dimensions entirely, or focused on binaries aligned with the midplane of the disk. Herein we present the first three-dimensional, high-resolution, local-shearing-box, inviscid hydrodynamical simulations of disk-embedded binaries over a range of orbital inclinations. We find that retrograde binaries can shrink up to 4 times as quickly as prograde binaries, and that all binaries not perfectly aligned (or anti-aligned) with the AGN disk are driven into alignment. An important consequence of this is that initially retrograde binaries will traverse the inclinations where von Zeipel–Lidov–Kozai oscillations can drive binary eccentricities to large values, potentially facilitating mergers. We also find that interactions with the AGN disk may excite eccentricities in retrograde binaries and cause the orbits of embedded binaries to precess.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Search for Neutron-Antineutron Transition in Large Neutrino Detectors

Neutron-antineutron transition is a baryon number violation process with $\Delta B=2$, providing a unique insight into the potential explanation of the baryon asymmetry in our university, especially in the context of post-sphaleron baryogenesis. Studies have been conducted in various neutron-rich environment including free neutron sources, neutron stars, and bound neutrons in large neutrino detectors, among which the large neutrino experiments obtained the most stringent constraints.

Wan, Linyan [Fermilab] (ORCID:0000000155246137)↗
DOE OSTI2023

Constraining Inputs to Realistic Kilonova Simulations through Comparison to Observed r -process Abundances

Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with the electromagnetic counterpart to GW170817 predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming that the observed source is prototypical, this inferred abundance pattern in turn must match r -process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun, as a practical prototype for the potentially universal abundance signature from neutron star mergers. We characterize the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a dynamical ejecta composition calculated using the FRDM2012 nuclear mass and FRLDM fission models with extremely neutron-rich ejecta ($Y_e$ = 0.035) along with moderately neutron-rich ($Y_e$ = 0.27) wind ejecta composition yields a wind-to-dynamical mass ratio of $M_w$/$M_d$ = 0.47, which best matches the observed AT2017gfo kilonova light curves while also producing the bestmatching abundance of neutron capture elements in the solar system, though, allowing for systematics, the ratio may be as high as of order unity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Properties of Hot Nuclear Matter

A fully quantitative description of the equilibrium and dynamical properties of hot nuclear matter will be needed for the interpretation of the available and forthcoming astrophysical data, providing information on the post-merger phase of a neutron star coalescence. We discuss the results of a recently developed theoretical model, based on a phenomenological nuclear Hamiltonian including two- and three-nucleon potentials, to study the temperature dependence of average and single-particle properties of nuclear matter relevant to astrophysical applications. The potential of the proposed approach for describing dissipative processes leading to the appearance of bulk viscosity in neutron star matter is also outlined.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Simulations of Radiation Hydrodynamics in Strongly Coupled Plasmas (Final Report)

Radiation propagates throughout the Universe and is often the only information we have from distant astronomical objects. Additionally, radiation plays a key role in systems ranging from the interior of stars to the establishment of shocks at the edges of forming galaxies to neutron star mergers. Here radiative processes, absorption, emission, and scattering, are fundamental to transport of momentum and energy, as well as, the overall evolution of these and other astrophysical systems. In some cases, radiation-dominant systems are also sufficiently dense and ionized, such that, the ions are strongly coupled (ie electrostatic energy >> thermal energy or Γ >> 1), which will affect radiative processes and the transport of radiation. During this project, we developed a computational ray-tracing algorithm to study radiation transport through a simulated neutron star capsule implosion under conditions relevant to Omega-60 laser experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Radius-to-frequency Mapping and FRB Frequency Drifts

We build a model of radius-to-frequency mapping in magnetospheres of neutron stars and apply it to frequency drifts observed in fast radio bursts (FRBs). We assume that an emission patch propagates along the dipolar magnetic field lines, producing coherent emission with frequency, direction, and polarization defined by the local magnetic field. The observed temporal evolution of the frequency depends on the relativistic effects of time contraction and the curvature of the magnetic field lines. The model generically produces linear scaling of the drift rate, $\dot{\omega}$ ∝-${\omega}$, matching both numerically and parametrically the rates observed in FBRs; a more complicated behavior of $\dot{\omega}$ is also possible. Fast rotating magnetospheres produce higher drifts rates for similar viewing parameters than the slowly rotating ones. In the case of repeaters, the same source may show variable drift patterns depending on the observing phase. We expect rotational of polarization position angle through a burst, though by smaller amount than in radio pulsars. All of these findings compare favorably with properties of FBRs, strengthening their possible loci in the magnetospheres of neutron stars.

Magnetars↗
DOE OSTI2023

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. In conclusion, our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bayesian Exploration of Phenomenological EoS of Neutron/Hybrid Stars with Recent Observations

The description of the stellar interior of compact stars remains as a big challenge for the nuclear astrophysics community. The consolidated knowledge is restricted to density regions around the saturation of hadronic matter ρ 0 =2.8 × 10 14 g cm -3 , regimes where our nuclear models are successfully applied. As one moves towards higher densities and extreme conditions up to the quark/gluons deconfinement, little can be said about the microphysics of the equation of state (EoS). Here, we employ a Markov Chain Monte Carlo (MCMC) strategy to access the variability at high density regions of polytropic piecewise models for neutron star (NS) EoS or possible hybrid stars, i.e., a NS with a small quark-matter core. With a fixed description of the hadronic matter for low density, below the nuclear saturation density, we explore a variety of models for the high density regimes leading to stellar masses near to 2.5 M ⊙ , in accordance with the observations of massive pulsars. The models are constrained, including the observation of the merger of neutrons stars from VIRGO-LIGO and with the pulsar observed by NICER. In addition, we also discuss the possibility of the use of a Bayesian power regression model with heteroscedastic error. The set of EoS from the Laser Interferometer Gravitational-Wave Observatory (LIGO) was used as input and treated as the data set for the testing case.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Constraining the Dense Matter Equation of State with Joint Analysis of NICER and LIGO/Virgo Measurements

The Neutron Star Interior Composition Explorer collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars (NSs). In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of NSs from the gravitational wave signal of the compact binary merger GW170817. Further, we compare the constraints on the EOS to those set by the recent measurement of a 2.14 Me pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a NS–black hole merger, which yields weaker constraints on the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Broad Grid of 2D Kilonova Emission Models

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M ⊙ and 0.05 to 0.3c, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (Y e ) ejecta with a robust r-process composition and either a spherical or lobed high-Y e ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamics of Laterally Propagating Flames in X-Ray Bursts. II. Realistic Burning and Rotation

In this work, we continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and thermal structure on the flame. We find that at lower rotation rates the flame becomes harder to ignite, whereas at higher rotation rates the nuclear burning is enhanced by increased confinement from the Coriolis force and the flame propagates steadily. At higher crustal temperatures, the flame moves more quickly and accelerates as it propagates through the atmosphere. If the temperature is too high, instead of a flame propagating across the surface the entire atmosphere burns uniformly. Our findings could have implications for the relationship between observed burst rise times and neutron star rotation and accretion rates. All of the software used for these simulations is freely available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Gravitational Wave Eigenfrequencies from Neutrino-driven Core-collapse Supernovae

Abstract Core-collapse supernovae (CCSNe) are predicted to produce gravitational waves (GWs) that may be detectable by Advanced LIGO/Virgo. These GW signals carry information from the heart of these cataclysmic events, where matter reaches nuclear densities. Recent studies have shown that it may be possible to infer the properties of the proto-neutron star (PNS) via GWs generated by hydrodynamic perturbations of the PNS. However, we lack a comprehensive understanding of how these relationships may change with the properties of CCSNe. In this work, we build a self-consistent suite of over 1000 exploding CCSNe from a grid of progenitor masses and metallicities combined with six different nuclear equations of state (EOS). Performing a linear perturbation analysis on each model, we compute the resonant GW frequencies of the PNS, and we motivate a time-agnostic method for identifying characteristic frequencies of the dominant GW emission. From this, we identify two characteristic frequencies, of the early- and late-time signal, that measure the surface gravity of the cold remnant neutron star, and simultaneously constrain the hot nuclear EOS. However, we find that the details of the CCSN model, such as the treatment of gravity or the neutrino transport, and whether it explodes, noticeably change the magnitude and evolution of the PNS eigenfrequencies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Implementation of Monte Carlo Transport in the General Relativistic SpEC Code

Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the last is a subdominant source of error at the accuracy reached by current simulations and for the interactions currently included in our code. We also provide tests of the most important features of this code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

El-CID: a filter for gravitational-wave electromagnetic counterpart identification

ABSTRACT As gravitational-wave (GW) interferometers become more sensitive and probe ever more distant reaches, the number of detected binary neutron star mergers will increase. However, detecting more events farther away with GWs does not guarantee corresponding increase in the number of electromagnetic counterparts of these events. Current and upcoming wide-field surveys that participate in GW follow-up operations will have to contend with distinguishing the kilonova (KN) from the ever increasing number of transients they detect, many of which will be consistent with the GW sky-localization. We have developed a novel tool based on a temporal convolutional neural network architecture, trained on sparse early-time photometry and contextual information for Electromagnetic Counterpart Identification (El-CID). The overarching goal for El-CID is to slice through list of new transient candidates that are consistent with the GW sky localization, and determine which sources are consistent with KNe, allowing limited target-of-opportunity resources to be used judiciously. In addition to verifying the performance of our algorithm on an extensive testing sample, we validate it on AT2017gfo – the only EM counterpart of a binary neutron star merger discovered to date – and AT2019npv – a supernova that was initially suspected as a counterpart of the GW event, GW190814, but was later ruled out after further analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2026

Direct inference of nuclear equation-of-state parameters from gravitational-wave observations

The observation of neutron star mergers with gravitational waves (GWs) has provided a new method to constrain the dense-matter equation of state (EOS) and to better understand its nuclear physics. However, inferring nuclear microphysics from GW observations necessitates the sampling of EOS model parameters that serve as input for each EOS used during the GW data analysis. The sampling of the EOS parameters requires solving the Tolman–Oppenheimer–Volkoff (TOV) equations a large number of times—a process that slows down each likelihood evaluation in the analysis on the order of a few seconds. Here, we employ emulators for the TOV equations built using multilayer perceptron neural networks to enable direct inference of nuclear EOS parameters from GW strain data. Our emulators allow us to rapidly solve the TOV equations, taking in EOS parameters and outputting the associated tidal deformability of a neutron star in only a few tens of milliseconds. We implement these emulators in PyCBC to directly infer the EOS parameters using the event GW170817, providing posteriors on these parameters informed solely by GWs. We benchmark these runs against analyses performed using the full TOV solver and find that the emulators achieve speed ups of nearly two orders of magnitude, with negligible differences in the recovered posteriors. Additionally, we constrain the slope and curvature of the symmetry energy at the 90% upper credible interval to be $L$ sym ≲ 106 MeV and $K$ sym ≲ 26 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Magnetars and axion-like particles: probes with the hard X-ray spectrum

Quiescent hard X-ray and soft gamma-ray emission from neutron stars constitute a promising frontier to explore axion-like-particles (ALPs). ALP production in the core peaks at energies of a few keV to a few hundreds of keV; subsequently, the ALPs escape and convert to photons in the magnetosphere. The emissivity goes as ~T 6 while the conversion probability is enhanced for large magnetic fields, making magnetars, with their high core temperatures and strong magnetic fields, ideal targets for probing ALPs. We compute the energy spectrum of photons resulting from conversion of ALPs in the magnetosphere and then compare it against hard X-ray data from NuSTAR, INTEGRAL, and XMM-Newton, for a set of eight magnetars for which such data exists. Upper limits are placed on the product of the ALP-nucleon and ALP-photon couplings. For the production in the core, we perform a calculation of the ALP emissivity in degenerate nuclear matter modeled by a relativistic mean field theory. The reduction of the emissivity due to improvements to the one-pion exchange approximation is incorporated, as is the suppression of the emissivity due to proton superfluidity in the neutron star core. A range of core temperatures is considered, corresponding to different models of the steady heat transfer from the core to the stellar surface. For the subsequent conversion, we solve the coupled differential equations mixing ALPs and photons in the magnetosphere. The conversion occurs due to a competition between the dipolar magnetic field and the photon refractive index induced by the external magnetic field. Semi-analytic expressions are provided alongside the full numerical results. Furthermore, we also present an analysis of the uncertainty on the axion limits we derive due to the uncertainties in the magnetar masses, nuclear matter equation of state, and the proton superfluid critical temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗