Neutron stars in astrophysical binary systems represent exciting sources for multimessenger astrophysics. A potential source of electromagnetic transients from compact binary systems is the neutron star ocean, the external fluid layer encasing a neutron star. We present a groundwork study into tidal waves in neutron star oceans and their consequences. Specifically, we investigate how oscillation modes in neutron star oceans can be tidally excited during compact binary inspirals and parabolic encounters. We find that neutron star oceans can sustain tidal waves with frequencies between 0.01 and 20 Hz. Our results suggest that tidally resonant neutron star ocean waves may serve as a neverbefore studied source of precursor electromagnetic emission prior to neutron star–black hole and binary neutron star mergers. If accompanied by electromagnetic flares, tidally resonant neutron star ocean waves, whose energy budget can reach 1046 erg, may serve as early warning signs (≳1 min before merger) for compact binary mergers. Similarly, excited ocean tidal waves will coincide with neutron star parabolic encounters. Depending on the neutron star ocean model and a flare emission scenario, tidally resonant ocean flares may be detectable by Fermi and Nuclear Spectroscopic Telescope Array (NuSTAR) out to ≳100 Mpc with detection rates as high as ∼7 yr−1 for binary neutron stars and ∼0.6 yr−1 for neutron star–black hole binaries. Observations of emission from neutron star ocean tidal waves along with gravitational waves will provide insight into the equation of state at the neutron star surface, the composition of neutron star oceans and crusts, and neutron star geophysics.
 NSFPAR ID:
 10330979
 Date Published:
 Journal Name:
 Classical and Quantum Gravity
 Volume:
 38
 Issue:
 13
 ISSN:
 02649381
 Page Range / eLocation ID:
 135008
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

ABSTRACT 
Abstract Gravitationalwave observations of binary neutron star mergers provide valuable information about neutron star structure and the equation of state of dense nuclear matter. Numerous methods have been proposed to analyze the population of observed neutron stars, and previous work has demonstrated the necessity of jointly fitting the astrophysical distribution and the equation of state in order to accurately constrain the equation of state. In this work, we introduce a new framework to simultaneously infer the distribution of binary neutron star masses and the nuclear equation of state using Gaussian mixture model density estimates, which mitigates some of the limitations previously used methods suffer from. Using our method, we reproduce previous projections for the expected precision of our joint mass distribution and equationofstate inference with tens of observations. We also show that mismodeling the equation of state can bias our inference of the neutron star mass distribution. While we focus on neutron star masses and matter effects, our method is widely applicable to population inference problems.

ABSTRACT Accurate measurements of the masses of neutron stars are necessary to test binary evolution models, and to constrain the neutron star equation of state. In pulsar binaries with no measurable postKeplerian parameters, this requires an accurate estimate of the binary system’s inclination and the radial velocity of the companion star by other means than pulsar timing. In this paper, we present the results of a new method for measuring this radial velocity using the binary synthesis code Icarus. This method relies on constructing a model spectrum of a tidally distorted, irradiated star as viewed for a given binary configuration. This method is applied to optical spectra of the newly discovered black widow PSR J1555–2908. By modeling the optical spectroscopy alongside optical photometry, we find that the radial velocity of the companion star is 397 ± 4 km s−1 (errors quoted at 95 per cent confidence interval), as well as a binary inclination of >75°. Combined with γray pulsation timing information, this gives a neutron star mass of 1.67$^{+0.15}_{0.09}$ M⊙ and a companion mass of 0.060$^{+0.005}_{0.003}$ M⊙, placing PSR J1555–2908 at the observed upper limit of what is considered a black widow system.

Summary We present a spatially varying Robin interface condition for solving fluid‐structure interaction problems involving incompressible fluid flows and nonuniform flexible structures. Recent studies have shown that for uniform structures with constant material and geometric properties, a constant one‐parameter Robin interface condition can improve the stability and accuracy of partitioned numerical solution procedures. In this work, we generalize the parameter to a spatially varying function that depends on the structure's local material and geometric properties, without varying the exact solution of the coupled fluid‐structure system. We present an algorithm to implement the Robin interface condition in an embedded boundary method for coupling a projection‐based incompressible viscous flow solver with a nonlinear finite element structural solver. We demonstrate the numerical effects of the spatially varying Robin interface condition using two example problems: a simplified model problem featuring a nonuniform Euler‐Bernoulli beam interacting with an inviscid flow and a generalized Turek‐Hron problem featuring a nonuniform, highly flexible beam interacting with a viscous laminar flow. Both cases show that a spatially varying Robin interface condition can clearly improve numerical accuracy (by up to two orders of magnitude in one instance) for the same computational cost. Using the second example problem, we also demonstrate and compare two models for determining the local value of the combination function in the Robin interface condition.

null (Ed.)Binary neutron star mergers provide a unique probe of the densematter equation of state (EoS) across a wide range of parameter space, from the zerotemperature EoS during the inspiral to the hightemperature EoS following the merger. In this paper, we implement a new model for calculating parametrized finitetemperature EoS effects into numerical relativity simulations. This "M* model" is based on a twoparameter approximation of the particle effective mass and includes the leadingorder effects of degeneracy in the thermal pressure and energy. We test our numerical implementation by performing evolutions of rotating single stars with zero and nonzero temperature gradients, as well as evolutions of binary neutron star mergers. We find that our new finitetemperature EoS implementation can support stable stars over many dynamical timescales. We also perform a first parameter study to explore the role of the M* parameters in binary neutron star merger simulations. All simulations start from identical initial data with identical cold EoSs, and differ only in the thermal part of the EoS. We find that both the thermal profile of the remnant and the postmerger gravitational wave signal depend on the choice of M* parameters, but that the total merger ejecta depends only weakly on the finitetemperature part of the EoS across a wide range of parameters. Our simulations provide a first step toward understanding how the finitetemperature properties of dense matter may affect future observations of binary neutron star mergers.more » « less