We present new numerical relativity results of neutron star mergers with chirp mass $1.188M_\odot$ and mass ratios $q=1.67$ and $q=1.8$ using finite-temperature equations of state (EOS), approximate neutrino transport and a subgrid model for magnetohydrodynamics-induced turbulent viscosity. The EOS are compatible with nuclear and astrophysical constraints and include a new microphysical model derived from ab-initio calculations based on the Brueckner-Hartree-Fock approach. We report for the first time evidence for accretion-induced prompt collapse in high-mass-ratio mergers, in which the tidal disruption of the companion and its accretion onto the primary star determine prompt black hole formation. As a result of the tidal disruption, an accretion disc of neutron-rich and cold matter forms with baryon masses ${\sim}0.15M_\odot$, and it is significantly heavier than the remnant discs in equal-masses prompt collapse mergers. Massive dynamical ejecta of order ${\sim}0.01M_\odot$ also originate from the tidal disruption. They are neutron rich and expand from the orbital plane with a crescent-like geometry. Consequently, bright, red and temporally extended kilonova emission is predicted from these mergers. Our results show that prompt black hole mergers can power bright electromagnetic counterparts for high-mass-ratio binaries, and that the binary mass ratio can be in principle constrained from multimessenger observations.

Thanks to recent measurements of tidal deformability and radius, the nuclear equation of state and structure of neutron stars are now better understood. Here, we show that through resonant tidal excitations in a binary inspiral, the neutron crust generically undergoes elastic-to-plastic transition, which leads to crust heating and eventually meltdown. This process could induce $\sim \mathcal{O}(0.1)$ phase shift in the gravitational waveform. Detecting the timing and induced phase shift of this crust meltdown will shed light on the crust structure, such as the core-crust transition density, which previous measurements are insensitive to. A direct search using GW170817 data has not found this signal, possibly due to limited signal-to-noise ratio. We predict that such signal may be observable with Advanced LIGO Plus and more likely with third-generation gravitational-wave detectors such as the Einstein Telescope and Cosmic Explorer.

The AT2017gfo kilonova counterpart of the binary neutron star merger event GW170817 was characterized by an early-time bright peak in optical and UV bands. Such blue kilonova is commonly interpreted as a signature of weak $r$-process nucleosynthesis in a fast expanding wind whose origin is currently debated. Numerical-relativity simulations with microphysical equations of state, approximate neutrino transport, and turbulent viscosity reveal a new mechanism that can power the blue kilonova. Spiral density waves in the remnant generate a characteristic wind of mass ${\sim}10^{-2}~M_{\odot}$ and velocity ${\sim}0.2$c. The ejected material has electron fraction mostly distributed above $0.25$ being partially reprocessed by hydrodynamic shocks in the expanding arms. The combination of dynamical ejecta and spiral-wave wind can account for solar system abundances of $r$-process elements and early-time observed light curves.

Observations of gravitational waves from inspiralling neutron star binaries---such as GW170817---can be used to constrain the nuclear equation of state by placing bounds on stellar tidal deformability. For slowly rotating neutron stars, the response to a weak quadrupolar tidal field is characterized by four internal-structure-dependent constants called "Love numbers." The tidal Love numbers $k_2^\text{el}$ and $k_2^\text{mag}$ measure the tides raised by the gravitoelectric and gravitomagnetic components of the applied field, and the rotational-tidal Love numbers $\mathfrak{f}^\text{o}$ and $\mathfrak{k}^\text{o}$ measure those raised by couplings between the applied field and the neutron star spin. In this work we compute these four Love numbers for perfect fluid neutron stars with realistic equations of state. We discover (nearly) equation-of-state independent relations between the rotational-tidal Love numbers and the moment of inertia, thereby extending the scope of I-Love-Q universality. We find that similar relations hold among the tidal and rotational-tidal Love numbers. These relations extend the applications of I-Love universality in gravitational-wave astronomy. As our findings differ from those reported in the literature, we derive general formulas for the rotational-tidal Love numbers in post-Newtonian theory and confirm numerically that they agree with our general-relativistic computations in the weak-field limit.

This paper considers the yielding response of a neutron star crust to smooth, unbalanced Maxwell stresses imposed at the core-crust boundary, and the coupling of the dynamic crust to the external magnetic field. Stress buildup and yielding in a magnetar crust is a global phenomenon: an elastic distortion radiating from one plastically deforming zone is shown to dramatically increase the creep rate in distant zones. Runaway creep to dynamical rates is shown to be possible, being enhanced by in situ heating and suppressed by thermal conduction and shearing of an embedded magnetic field. A global and time-dependent model of elastic, plastic, magnetic, and thermal evolution is developed. Fault-like structures develop naturally, and a range of outburst timescales is observed. Transient events with time profiles similar to giant magnetar flares (millisecond rise, $\sim$ 0.1 s duration, and decaying power-law tails) result from runaway creep that starts in localized sub-km-sized patches and spreads across the crust. A one-dimensional model of stress relaxation in the vertically stratified crust shows that a modest increase in applied stress allows embedded magnetic shear to escape the star over $\sim$ 3-10 ms, dissipating greater energy if the exterior field is already sheared. Several such zones coupled to each other naturally yield a burst of duration $\sim$ 0.1 s, as is observed over a wide range of burst energies. The collective interaction of many plastic zones forces an overstability of global elastic modes of the crust, consistent with QPO activity extending over $\gtrsim$ 100 s. Giant flares probably involve sudden meltdown in localized zones, with high-frequency ($\gg$ 100 Hz) QPOs corresponding to standing Alfv\'en waves within these zones.

Scalar-tensor theories of gravity are extensions of General Relativity (GR) including an extra, nonminimally coupled scalar degree of freedom. A wide class of these theories, albeit indistinguishable from GR in the weak field regime, predicts a radically different phenomenology for neutron stars, due to a nonperturbative, strong-field effect referred to as spontaneous scalarization. This effect is known to occur in theories where the effective linear coupling $\beta_0$ between the scalar and matter fields is sufficiently negative, i.e. $\beta_0 \lesssim -4.35$, and has been strongly constrained by pulsar timing observations. In the test-field approximation, spontaneous scalarization manifests itself as a tachyonic-like instability. Recently, it was argued that, in theories where $\beta_0>0$, a similar instability would be triggered by sufficiently compact neutron stars obeying realistic equations of state. In this work we investigate the endstate of this instability for some representative coupling functions with $\beta_0>0$. This is done both through an energy balance analysis of the existing equilibrium configurations, and by numerically determining the nonlinear Cauchy development of unstable initial data. We find that, contrary to the $\beta_0<0$ case, the final state of the instability is highly sensitive to the details of the coupling function, varying from gravitational collapse to spontaneous scalarization. In particular, we show, for the first time, that spontaneous scalarization can happen in theories with $\beta_0>0$, which could give rise to novel astrophysical tests of the theory of gravity.

We suggest that the redshift of photons traveling from past to future null infinity through a collapsing object could provide an observational signature capable of differentiating between the formation of a globally naked singularity and the formation of an event horizon. Supporting evidence for this idea is drawn from the analysis of photons with zero angular momentum through the center of a collapsing spherical dust cloud. We show that the frequency shift as a function of proper time with respect to stationary observers has distinct features depending on whether the object collapses to a black hole or a naked singularity.