Gravitational collapse of plasmas in General Relativity

We provide a covariant derivation of plasma physics coupled to gravitation by utilizing the 3+1 formulation of general relativity, including a discussion of the Lorentz force law. We then reduce the system to the spherically symmetric case and show that all regions of the spacetime can be represented in a single coordinate system, thus revoking the need for junction conditions. We further show that the region exterior to the collapsing region is naturally described by the charged Vaidya spacetime in non-null coordinates.Comment: Talk given at the Spanish Relativity Meeting, Tenerife, September 200

Spherically Symmetric Gravitational Collapse of Perfect Fluids

Formulating a perfect fluid filled spherically symmetric metric utilizing the 3+1 formalism for general relativity, we show that the metric coefficients are completely determined by the mass-energy distribution, and its time rate of change on an initial spacelike hypersurface. Rather than specifying Schwarzschild coordinates for the exterior of the collapsing region, we let the interior dictate the form of the solution in the exterior, and thus both regions are found to be written in one coordinate patch. This not only alleviates the need for complicated matching schemes at the interface, but also finds a new coordinate system for the Schwarzschild spacetime expressed in generalized Painleve-Gullstrand coordinates.Comment: 3 pages, To appear in the proceedings of the eleventh Marcel Grossmann meeting on general relativity (MGXI), 23-29 July, 2006, Berli

Three-dimensional shapelets and an automated classification scheme for dark matter haloes

We extend the two-dimensional Cartesian shapelet formalism to d-dimensions. Concentrating on the three-dimensional case, we derive shapelet-based equations for the mass, centroid, root-mean-square radius, and components of the quadrupole moment and moment of inertia tensors. Using cosmological N-body simulations as an application domain, we show that three-dimensional shapelets can be used to replicate the complex sub-structure of dark matter halos and demonstrate the basis of an automated classification scheme for halo shapes. We investigate the shapelet decomposition process from an algorithmic viewpoint, and consider opportunities for accelerating the computation of shapelet-based representations using graphics processing units (GPUs).Comment: 19 pages, 11 figures, accepted for publication in MNRA

Mapping the Universe Expansion: Enabling percent-level measurements of the Hubble Constant with a single binary neutron-star merger detection

The joint observation of the gravitational-wave and electromagnetic signal from the binary neutron-star merger GW170817 allowed for a new independent measurement of the Hubble constant $H_0$, albeit with an uncertainty of about 15\% at 1$\sigma$. Observations of similar sources with a network of future detectors will allow for more precise measurements of $H_0$. These, however, are currently largely limited by the intrinsic degeneracy between the luminosity distance and the inclination of the source in the gravitational-wave signal. We show that the higher-order modes in gravitational waves can be used to break this degeneracy in astrophysical parameter estimation in both the inspiral and post-merger phases of a neutron star merger. We show that for systems at distances similar to GW170817, this method enables percent-level measurements of $H_0$ with a single detection. This would permit the study of time variations and spatial anisotropies of $H_0$ with unprecedented precision. We investigate how different network configurations affect measurements of $H_0$, and discuss the implications in terms of science drivers for the proposed 2.5- and third-generation gravitational-wave detectors. Finally, we show that the precision of $H_0$ measured with these future observatories will be solely limited by redshift measurements of electromagnetic counterparts

Gravitational collapse of spherically symmetric plasmas in Einstein-Maxwell spacetimes

We utilize a recent formulation of a spherically symmetric spacetime endowed with a general decomposition of the energy momentum tensor [Phys. Rev. D, 75, 024031 (2007)] to derive equations governing spherically symmetric distributions of electromagnetic matter. We show the system reduces to the Reissner-Nordstrom spacetime in general, spherically symmetric coordinates in the vacuum limit. Furthermore, we show reduction to the charged Vaidya spacetime in non-null coordinates when certain equations of states are chosen. A model of gravitational collapse is discussed whereby a charged fluid resides within a boundary of finite radial extent on the initial hypersurface, and is allowed to radiate charged particles. Our formalism allows for the discussion of all regions in this model without the need for complicated matching schemes at the interfaces between successive regions. As further examples we consider the collapse of a thin shell of charged matter onto a Reissner-Nordstrom black hole. Finally, we reduce the entire system of equations to the static case such that we have the equations for hydrostatic equilibrium of a charged fluid.Comment: Accepted for publication in Phys. Rev.

Spherically Symmetric Gravitational Collapse of General Fluids

We express Einstein's field equations for a spherically symmetric ball of general fluid such that they are conducive to an initial value problem. We show how the equations reduce to the Vaidya spacetime in a non-null coordinate frame, simply by designating specific equations of state. Furthermore, this reduces to the Schwarzschild spacetime when all matter variables vanish. We then describe the formulation of an initial value problem, whereby a general fluid ball with vacuum exterior is established on an initial spacelike slice. As the system evolves, the fluid ball collapses and emanates null radiation such that a region of Vaidya spacetime develops. Therefore, on any subsequent spacelike slice there exists three regions; general fluid, Vaidya and Schwarzschild, all expressed in a single coordinate patch with two free-boundaries determined by the equations. This implies complicated matching schemes are not required at the interfaces between the regions, instead, one simply requires the matter variables tend to the appropriate equations of state. We also show the reduction of the system of equations to the static cases, and show staticity necessarily implies zero ``heat flux''. Furthermore, the static equations include a generalization of the Tolman-Oppenheimer-Volkoff equations for hydrostatic equilibrium to include anisotropic stresses in general coordinates.Comment: 11 pages, 3 figures, submitted to Phys. Rev.

Pressure gradients, shell crossing singularities and acoustic oscillations - application to inhomogeneous cosmological models

Inhomogeneous cosmological models have recently become a very interesting alternative to standard cosmology. This is because these models are able to fit cosmological observations without the need for dark energy. However, due to inhomogeneity and pressure-less matter content, these models can suffer from shell crossing singularities. These singularities occur when two shell of dust collide with each other leading to infinite values of the density. In this Letter we show that if inhomogeneous pressure is included then these singularities can be prevented from occurring over the period of structure formation. Thus, a simple incorporation of a gradient of pressure allows for more comprehensive studies of inhomogeneous cosmological models and their application to cosmology.Comment: 5 pages, 3 figures. Accepted for publication in MNRA

Linking the rates of neutron star binaries and short gamma-ray bursts

Short gamma-ray bursts are believed to be produced by both binary neutron star (BNS) and neutron star-black hole (NSBH) mergers. We use current estimates for the BNS and NSBH merger rates to calculate the fraction of observable short gamma-ray bursts produced through each channel. This allows us to constrain merger rates of BNS to $\mathcal{R}_{\rm{BNS}}=384^{+431}_{-213}{\rm{Gpc}^{-3} \rm{yr}^{-1}}$ ($90\%$ credible interval), a $16\%$ decrease in the rate uncertainties from the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. Assuming a top-hat emission profile with a large Lorentz factor, we constrain the average opening angle of gamma-ray burst jets produced in BNS mergers to $\approx 15^\circ$. We also measure the fraction of BNS and NSBH mergers that produce an observable short gamma-ray burst to be $0.02^{+0.02}_{-0.01}$ and $0.01 \pm 0.01$, respectively and find that $\gtrsim 40\%$ of BNS mergers launch jets (90\% confidence). We forecast constraints for future gravitational-wave detections given different modelling assumptions, including the possibility that BNS and NSBH jets are different. With $24$ BNS and $55$ NSBH observations, expected within six months of the LIGO-Virgo-KAGRA network operating at design sensitivity, it will be possible to constrain the fraction of BNS and NSBH mergers that launch jets with $10\%$ precision. Within a year of observations, we can determine whether the jets launched in NSBH mergers have a different structure than those launched in BNS mergers and rule out whether $\gtrsim 80\%$ of binary neutron star mergers launch jets. We discuss the implications of future constraints on understanding the physics of short gamma-ray bursts and binary evolution.Comment: Accepted in Physical Review D: 13 pages, 5 figure

Gravitational waves from Sco X-1: A comparison of search methods and prospects for detection with advanced detectors

The low-mass X-ray binary Scorpius X-1 (Sco X-1) is potentially the most luminous source of continuous gravitational-wave radiation for interferometers such as LIGO and Virgo. For low-mass X-ray binaries this radiation would be sustained by active accretion of matter from its binary companion. With the Advanced Detector Era fast approaching, work is underway to develop an array of robust tools for maximizing the science and detection potential of Sco X-1. We describe the plans and progress of a project designed to compare the numerous independent search algorithms currently available. We employ a mock-data challenge in which the search pipelines are tested for their relative proficiencies in parameter estimation, computational efficiency, robust- ness, and most importantly, search sensitivity. The mock-data challenge data contains an ensemble of 50 Scorpius X-1 (Sco X-1) type signals, simulated within a frequency band of 50-1500 Hz. Simulated detector noise was generated assuming the expected best strain sensitivity of Advanced LIGO and Advanced VIRGO ($4 \times 10^{-24}$ Hz$^{-1/2}$). A distribution of signal amplitudes was then chosen so as to allow a useful comparison of search methodologies. A factor of 2 in strain separates the quietest detected signal, at $6.8 \times 10^{-26}$ strain, from the torque-balance limit at a spin frequency of 300 Hz, although this limit could range from $1.2 \times 10^{-25}$ (25 Hz) to $2.2 \times 10^{-26}$ (750 Hz) depending on the unknown frequency of Sco X-1. With future improvements to the search algorithms and using advanced detector data, our expectations for probing below the theoretical torque-balance strain limit are optimistic.Comment: 33 pages, 11 figure