The binary neutron star merger responsible for the gravitational wave event, GW170817, strengthened the merger association with short-duration gamma-ray bursts (SGRBs) following the detection of the SGRB counterpart, GRB 170817A. Here we consider the constraints on a population of local gamma-ray bursts with moderately short duration ($T_{90}<4 \textrm{ s}$) and within $d < 200 \textrm{ Mpc}$, that may have originated from similar compact binary mergers. Using well localised gamma-ray bursts from $\sim14.5 \textrm{ years}$ of Swift/Burst Alert Telescope monitoring, we find no events with high likelihood of being in this distance range, and place an upper limit for the all-sky rate of such events of $<4\ \textrm{y}^{-1}$. For Fermi/Gamma-ray Burst Monitor (GBM) and CGRO/Burst And Transient Source Experiment (BATSE) detected bursts, where the localisation has considerably larger uncertainties, we cross-correlated with 2MASS Redshift Survey galaxies at $d<100 \textrm{ Mpc}$, obtaining a weaker constraint of $< 12\ \textrm{y}^{-1}$. A separate correlation search from the GBM and BATSE bursts for giant flares originating from soft gamma-ray repeaters in nearby galaxies ($d < 11 \textrm{ Mpc} $ ) yields an upper limit of $<3\ \textrm{y}^{-1}$.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Following the faint gamma-ray burst, GRB 170817A, coincident with a gravitational wave-detected binary neutron star merger at $d\sim40$ Mpc, we consider the constraints on a local population of faint short duration GRBs (defined here broadly as $T_{90}<4$ s). We review proposed low-redshift short-GRBs and consider statistical limits on a $d\lessapprox200$ Mpc population using Swift/Burst Alert Telescope (BAT), Fermi/Gamma-ray Burst Monitor (GBM), and Compton Gamma-Ray Observatory (CGRO) Burst and Transient Source Experiment (BATSE) GRBs. Swift/BAT short-GRBs give an upper limit for the all-sky rate of $<4$ y$^{-1}$ at $d<200$ Mpc, corresponding to $<5$% of SGRBs. Cross-correlation of selected CGRO/BATSE and Fermi/GBM GRBs with $d<100$ Mpc galaxy positions returns a weaker constraint of $\lessapprox12\ {\rm y^{-1}}$. A separate search for correlations due to SGR giant flares in nearby ($d<11$ Mpc) galaxies finds an upper limit of $<3\ {\rm y^{-1}}$. Our analysis suggests that GRB 170817A-like events are likely to be rare in existing SGRB catalogues. The best candidate for an analogue remains GRB 050906, where the Swift/BAT location was consistent with the galaxy IC0327 at $d\approx132$ Mpc. If binary neutron star merger rates are at the high end of current estimates, then our results imply that at most a few percent will be accompanied by detectable gamma-ray flashes in the forthcoming LIGO/Virgo science runs.

We calculate the electromagnetic signal of a gamma-ray flare coming from the surface of a neutron star shortly before merger with a black hole companion. Using a new version of the Monte Carlo radiation transport code Pandurata that incorporates dynamic spacetimes, we integrate photon geodesics from the neutron star surface until they reach a distant observer or are captured by the black hole. The gamma-ray light curve is modulated by a number of relativistic effects, including Doppler beaming and gravitational lensing. Because the photons originate from the inspiraling neutron star, the light curve closely resembles the corresponding gravitational waveform: a chirp signal characterized by a steadily increasing frequency and amplitude. We propose to search for these electromagnetic chirps using matched filtering algorithms similar to those used in LIGO data analysis.

We examine the effects of higher-order multipole contributions of rotating neutron star (NS) spacetimes on the propagation of corrugation (c-)modes within a thin accretion disk. We find that the Lense-Thirring precession frequency, which determines the propagation region of the low-frequency fundamental corrugation modes, can experience a turnover allowing for c-modes to become self-trapped for sufficiently high dimensionless spin $j$ and quadrupole rotational deformability $\alpha$. If such self-trapping c-modes can be detected, e.g. through phase-resolved spectroscopy of the iron line for a high-spin low-mass accreting neutron star, this could potentially constrain the spin-induced NS quadrupole and the NS equation of state.

Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative interactions. In this Letter, we develop the "slimplectic" integrator, a new type of numerical integrator that shares many of the benefits of traditional symplectic integrators yet is applicable to general nonconservative systems. We utilize a fixed time-step variational integrator formalism applied to the principle of stationary nonconservative action developed in Galley, 2013; Galley, Tsang & Stein, 2014. As a result, the generalized momenta and energy (Noether current) evolutions are well-tracked. We discuss several example systems, including damped harmonic oscillators, Poynting-Robertson drag, and gravitational radiation reaction, by utilizing our new publicly available code to demonstrate the slimplectic integrator algorithm. Slimplectic integrators are well-suited for integrations of systems where nonconservative effects play an important role in the long-term dynamical evolution. As such they are particularly appropriate for cosmological or celestial N-body dynamics problems where nonconservative interactions, e.g. gas interactions or dissipative tides, can play an important role.

We further develop a recently introduced variational principle of stationary action for problems in nonconservative classical mechanics and extend it to classical field theories. The variational calculus used is consistent with an initial value formulation of physical problems and allows for time-irreversible processes, such as dissipation, to be included at the level of the action. In this formalism, the equations of motion are generated by extremizing a nonconservative action $\mathcal{S}$, which is a functional of a doubled set of degrees of freedom. The corresponding nonconservative Lagrangian contains a potential $K$ which generates nonconservative forces and interactions. Such a nonconservative potential can arise in several ways, including from an open system interacting with inaccessible degrees of freedom or from integrating out or coarse-graining a subset of variables in closed systems. We generalize Noether's theorem to show how Noether currents are modified and no longer conserved when $K$ is non-vanishing. Consequently, the nonconservative aspects of a physical system are derived solely from $K$. We show how to use the formalism with examples of nonconservative actions for discrete systems including forced damped harmonic oscillators, radiation reaction on an accelerated charge, and RLC circuits. We present examples for nonconservative classical field theories. Our approach naturally allows for irreversible thermodynamic processes to be included in an unconstrained variational principle. We present the nonconservative action for a Navier-Stokes fluid including the effects of viscous dissipation and heat diffusion, as well as an action that generates the Maxwell model for viscoelastic materials, which can be easily generalized to more realistic rheological models. We show that the nonconservative action can be derived as the classical limit of a more complete quantum theory.

We show that the first order (non co-orbital) corotation torques are significantly modified by entropy gradients in a non-barotropic protoplanetary disk. Such non-barotropic torques can dramatically alter the balance that, for barotropic cases, results in the net eccentricity damping for giant gap-clearing planets embedded in the disk. We demonstrate that stellar illumination can heat the gap enough for the planet's orbital eccentricity to instead be excited. We also discuss the "Eccentricity Valley" noted in the known exoplanet population, where low-metallicity stars have a deficit of eccentric planets between $\sim 0.1$ and $\sim 1$ AU compared to metal-rich systems (Dawson & Murray-Clay 2013). We show that this feature in the planet distribution may be due to the self-shadowing of the disk by a rim located at the dust sublimation radius $\sim 0.1$ AU, which is known to exist for several T Tauri systems. In the shadowed region between $\sim 0.1$ and $\sim 1$ AU lack of gap insolation allows disk interactions to damp eccentricity. Outside such shadowed regions stellar illumination can heat the planetary gaps and drive eccentricity growth for giant planets. We suggest that the self-shadowing does not arise at higher metallicity due to the increased optical depth of the gas interior to the dust sublimation radius.

I derive a fully analytic expression for the linear corotation torque to first order in eccentricity for planets in non-barotropic protoplanetary disks, taking into account the effect of disk entropy gradients. This torque formula is applicable to both the co-orbital corotation torque and the non co-orbital corotation torques -- for planets in orbits with non-zero eccentricity -- in disks where the thermal diffusivity and viscosity are sufficient to maintain linearity of these interactions. While the co-orbital corotation torque is important for migration of planets in Type I migration, the non co-orbital corotation torque plays an important role in the eccentricity evolution of giant planets that have opened gaps in the disk. The presence of an entropy gradient in the disk can significantly modify the corotation torque in both these cases.

Using a fast semi-analytic raytracing code, we study the variability of relativistically broadened iron lines due to discoseismic oscillations concentrated in the inner-most regions of accretion discs around black holes. The corrugation mode, or c-mode, is of particular interest as its natural frequency corresponds well to the ~0.1-15Hz range observed for low-frequency quasi-periodic oscillations (LFQPOs) for lower spins. Comparison of the oscillation phase dependent variability and QPO-phase stacked iron line observations will allow such discoseismic models to be confirmed or ruled out as a source of particular LFQPOs. The spectral range and frequency of the variability of the iron line due to corrugation modes can also potentially be used to constrain the black hole spin if observed with sufficient temporal and spectral resolution.

We examine timing noise in both magnetars and regular pulsars, and find that there exists a component of the timing noise with strong magnetic field dependence above B ~ 10^{12.5} G. The dependence of the timing noise floor on the magnetic field is also reflected in the smallest observable glitch size. We find that magnetospheric torque variation cannot explain this component of timing noise. We calculate the moment of inertia of the magnetic field outside of a neutron star and show that this timing noise component may be due to variation of this moment of inertia, and could be evidence of rapid global magnetospheric variability.

We demonstrate that resonant shattering flares can occur during close passages of neutron stars in eccentric or hyperbolic encounters. We provide updated estimates for the rate of close encounters of compact objects in dense stellar environments, which we find are substantially lower than given in previous works. While such occurrences are rare, we show that shattering flares can provide a strong electromagnetic counterpart to the gravitational wave bursts expected from such encounters, allowing triggered searches for these events to occur.

Magnetars are neutron stars showing dramatic X-ray and soft $\gamma$-ray outbursting behaviour that is thought to be powered by intense internal magnetic fields. Like conventional young neutron stars in the form of radio pulsars, magnetars exhibit "glitches" during which angular momentum is believed to be transferred between the solid outer crust and the superfluid component of the inner crust. Hitherto, the several hundred observed glitches in radio pulsars and magnetars have involved a sudden spin-up of the star, due presumably to the interior superfluid rotating faster than the crust. Here we report on X-ray timing observations of the magnetar 1E 2259+586 which we show exhibited a clear "anti-glitch" -- a sudden spin down. We show that this event, like some previous magnetar spin-up glitches, was accompanied by multiple X-ray radiative changes and a significant spin-down rate change. This event, if of origin internal to the star, is unpredicted in models of neutron star spin-down and is suggestive of differential rotation in the neutron star, further supporting the need for a rethinking of glitch theory for all neutron stars.

The physical origin of high-frequency QPOs (HFQPOs) in black-hole X-ray binaries remains an enigma despite many years of detailed observational studies. Although there exists a number of models for HFQPOs, many of these are simply "notions" or "concepts" without actual calculation derived from fluid or disk physics. Future progress requires a combination of numerical simulations and semi-analytic studies to extract physical insights. We review recent works on global oscillation modes in black-hole accretion disks, and explain how, with the help of general relativistic effects, the energy stored in the disk differential rotation can be pumped into global spiral density modes in the disk, making these modes grow to large amplitudes under certain conditions ("corotational instability"). These modes are robust in the presence of disk magnetic fields and turbulence. The computed oscillation mode frequencies are largely consistent with the observed values for HFQPOs in BH X-ray binaries. The approximate 2:3 frequency ratio is also expected from this model. The connection of HFQPOs with other disk properties (such as production of episodic jets) is also discussed.

The resonant excitation of neutron star (NS) modes by tides is investigated as a source of short gamma-ray burst (sGRB) precursors. We find that the driving of a crust-core interface mode can lead to shattering of the NS crust, liberating ~10^46-10^47 erg of energy seconds before the merger of a NS-NS or NS-black hole binary. Such properties are consistent with Swift/BAT detections of sGRB precursors, and we use the timing of the observed precursors to place weak constraints on the crust equation of state. We describe how a larger sample of precursor detections could be used alongside coincident gravitational wave detections of the inspiral by Advanced LIGO class detectors to probe the NS structure. These two types of observations nicely complement one another, since the former constrains the equation of state and structure near the crust-core boundary, while the latter is more sensitive to the core equation of state.

Waves reflected by the inner edge of a protoplanetary disk are shown to significantly modify Type I migration, even allowing the trapping of planets near the inner disk edge for small planets in a range of disk parameters. This may inform the distribution planets close to their central stars, as observed by the Kepler mission.

We study the effect of corotation resonance on the inertial-acoustic oscillations (p-modes) of black-hole accretion discs. Previous works have shown that for barotropic flows (where the pressure depends only on the density), wave absorption at the corotation resonance can lead to mode growth when the disc vortensity, $\zeta=\kappa^2/(2\Omega\Sigma)$ (where $\Omega, \kappa, \Sigma$ are the rotation rate, radial epicyclic frequency and surface density of the disc, respectively), has a positive gradient at the corotation radius. Here we generalize the analysis of the corotation resonance effect to non-barotropic fluids. We show that the mode instability criterion is modified by the finite radial Brunt-V\"as\"al\"a frequency of the disc. We derive an analytic expression for the reflectivity when a density wave impinges upon the corotation barrier, and calculate the frequencies and growth rates of global p-modes for disc models with various $\alpha$-viscosity parameterizations. We find that for disc fluids with constant adiabatic index $\Gamma$, super-reflection and mode growth depend on the gradient of the effective vortensity, $\zeta_{\rm eff} = \zeta/S^{2/\Gamma}$ (where $S \equiv P/\Sigma^{\Gamma}$ measures the entropy): when $d\zeta_{\rm eff}/dr > 0$ at the corotation radius, wave absorption leads to amplification of the p-mode. Our calculations show that the lowest-order p-modes with azimuthal wave number $m=2, 3, 4,...$ have the largest growth rates, with the frequencies approximately in (but distinct from) the $2:3:4...$ commensurate ratios. We discuss the implications of our results for the high-frequency quasi-periodic oscillations observed in accreting black-hole systems.

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the inner-most regions of accretion discs or boundary layers. Two simple models are considered in this paper: The magnetosphere-disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star-disc model involves a Keplerian disc terminated at the stellar atomosphere with high density and small density scale height. We find that the interface modes at the magnetosphere-disc boundary are generally unstable due to Rayleigh-Taylor and/or Kelvin-Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh-Taylor instability and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin-Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh-Taylor or Kelvin-Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.

We study the global stability of non-axisymmetric p-modes (also called inertial-acoustic modes) trapped in the inner-most regions of accretion discs around black holes. We show that the lowest-order (highest-frequency) p-modes, with frequencies $\omega=(0.5-0.7) m\Omega_{\rm ISCO}$, can be overstable due to general relativistic effects, according to which the radial epicyclic frequency is a non-monotonic function of radius near the black hole. The mode is trapped inside the corotation resonance radius and carries a negative energy. The mode growth arises primarily from wave absorption at the corotation resonance, and the sign of the wave absorption depends on the gradient of the disc vortensity. When the mode frequency is sufficiently high, such that the slope of the vortensity is positive at corotation positive wave energy is absorbed at the resonance, leading to the growth of mode amplitude. We also study how the rapid radial inflow at the inner edge of the disc affects the mode trapping and growth. Our analysis of the behavior of the fluid perturbations in the transonic flow near the ISCO indicates that, while the inflow tends to damp the mode, the damping effect is sufficiently small under some conditions so that net mode growth can still be achieved. We further clarify the role of the Rossby wave instability and show that it does not operate for black hole accretion discs with smooth-varying vortensity profiles. Overstable non-axisymmetric p-modes driven by the corotational instability provide a plausible explanation for the high-frequency (> 100 Hz) quasi-periodic oscillations (HFQPOs) observed from a number of black-hole X-ray binaries in the very high state. The absence of HFQPOs in the soft (thermal) state may result from mode damping due to the radial infall at the ISCO.

Diskoseismic c-modes in accretion discs have been invoked to explain low-frequency variabilities observed in black-hole X-ray binaries. These modes are trapped in the inner-most region of the disc and have frequencies much lower than the rotation frequency at the disc inner radius. We show that because the trapped waves can tunnel through the evanescent barrier to the corotational wave zone, the c-modes are damped due to wave absorption at the corotation resonance. We calculate the corotational damping rates of various c-modes using the WKB approximation. The damping rate varies widely depending on the mode frequency, the black hole spin parameter and the disc sound speed, and is generally much less than 10% of the mode frequency. A sufficiently strong excitation mechanism is needed to overcome this corotational damping and make the mode observable.

In differentially rotating discs with no self-gravity, density waves cannot propagate around the corotation, where the wave pattern rotation speed equals the fluid rotation rate. Waves incident upon the corotation barrier may be super-reflected (commonly referred to as corotation amplifier), but the reflection can be strongly affected by wave absorptions at the corotation resonance/singularity. The sign of the absorption is related to the Rossby wave zone very near the corotation radius. We derive the explicit expressions for the complex reflection and transmission coefficients, taking into account wave absorption at the corotation resonance. We show that for generic discs, this absorption plays a much more important role than wave transmission across the corotation barrier. Depending on the sign of the gradient of the specific vorticity of the disc the corotation resonance can either enhance or diminish the super-reflectivity, and this can be understood in terms of the location of the Rossby wave zone relative to the corotation radius. Our results provide the explicit conditions (in terms of disc thickness, rotation profile and specific vorticity gradient) for which super-reflection can be achieved. Global overstable disc modes may be possible for discs with super-reflection at the corotation barrier.