Impact of inter-correlated initial binary parameters on double black hole and neutron star mergers

The distributions of the initial main-sequence binary parameters are one of the key ingredients in obtaining evolutionary predictions for compact binary (BH-BH / BH-NS / NS-NS) merger rates. Until now, such calculations were done under the assumption that initial binary parameter distributions were independent. Here, we implement empirically derived inter-correlated distributions of initial binary parameters primary mass (M1), mass ratio (q), orbital period (P), and eccentricity (e). Unexpectedly, the introduction of inter-correlated initial binary parameters leads to only a small decrease in the predicted merger rates by a factor of 2 $-$ 3 relative to the previously used non-correlated initial distributions. The formation of compact object mergers in the isolated classical binary evolution favors initial binaries with stars of comparable masses (q = 0.5 $-$ 1) at intermediate orbital periods (log P (days) = 2 $-$ 4). New distributions slightly shift the mass ratios towards smaller values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods only negligibly increase the number of progenitors. Additionally, we discuss the uncertainty of merger rate predictions associated with possible variations of the massive-star initial mass function (IMF). We argue that evolutionary calculations should be normalized to a star formation rate (SFR) that is obtained from the observed amount of UV light at wavelength 1500{\AA} (SFR indicator). In this case, contrary to recent reports, the uncertainty of the IMF does not affect the rates by more than a factor of 2. Any change to the IMF slope for massive stars requires a change of SFR in a way that counteracts the impact of IMF variations on the merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play.Comment: accepted for publication in A&

Stream network analysis and geomorphic flood plain mapping from orbital and suborbital remote sensing imagery application to flood hazard studies in central Texas

The author has identified the following significant results. Development of a quantitative hydrogeomorphic approach to flood hazard evaluation was hindered by (1) problems of resolution and definition of the morphometric parameters which have hydrologic significance, and (2) mechanical difficulties in creating the necessary volume of data for meaningful analysis. Measures of network resolution such as drainage density and basin Shreve magnitude indicated that large scale topographic maps offered greater resolution than small scale suborbital imagery and orbital imagery. The disparity in network resolution capabilities between orbital and suborbital imagery formats depends on factors such as rock type, vegetation, and land use. The problem of morphometric data analysis was approached by developing a computer-assisted method for network analysis. The system allows rapid identification of network properties which can then be related to measures of flood response

Explaining LIGO's observations via isolated binary evolution with natal kicks

We compare binary evolution models with different assumptions about black-hole natal kicks to the first gravitational-wave observations performed by the LIGO detectors. Our comparisons attempt to reconcile merger rate, masses, spins, and spin-orbit misalignments of all current observations with state-of-the-art formation scenarios of binary black holes formed in isolation. We estimate that black holes (BHs) should receive natal kicks at birth of the order of $\sigma\simeq 200$ (50) km/s if tidal processes do (not) realign stellar spins. Our estimate is driven by two simple factors. The natal kick dispersion $\sigma$ is bounded from above because large kicks disrupt too many binaries (reducing the merger rate below the observed value). Conversely, the natal kick distribution is bounded from below because modest kicks are needed to produce a range of spin-orbit misalignments. A distribution of misalignments increases our models' compatibility with LIGO's observations, if all BHs are likely to have natal spins. Unlike related work which adopts a concrete BH natal spin prescription, we explore a range of possible BH natal spin distributions. Within the context of our models, for all of the choices of $\sigma$ used here and within the context of one simple fiducial parameterized spin distribution, observations favor low BH natal spin.Comment: 19 pages, 14 figures, as published in PR

Double Compact Objects III: Gravitational Wave Detection Rates

The unprecedented range of second-generation gravitational-wave (GW) observatories calls for refining the predictions of potential sources and detection rates. The coalescence of double compact objects (DCOs)---i.e., neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS), and black hole-black hole (BH-BH) binary systems---is the most promising source of GWs for these detectors. We compute detection rates of coalescing DCOs in second-generation GW detectors using the latest models for their cosmological evolution, and implementing inspiral-merger-ringdown (IMR) gravitational waveform models in our signal-to-noise ratio calculations. We find that: (1) the inclusion of the merger/ringdown portion of the signal does not significantly affect rates for NS-NS and BH-NS systems, but it boosts rates by a factor $\sim 1.5$ for BH-BH systems; (2) in almost all of our models BH-BH systems yield by far the largest rates, followed by NS-NS and BH-NS systems, respectively, and (3) a majority of the detectable BH-BH systems were formed in the early Universe in low-metallicity environments. We make predictions for the distributions of detected binaries and discuss what the first GW detections will teach us about the astrophysics underlying binary formation and evolution.Comment: published in ApJ, 19 pages, 11 figure

Regge Calculus in Teleparallel Gravity

In the context of the teleparallel equivalent of general relativity, the Weitzenbock manifold is considered as the limit of a suitable sequence of discrete lattices composed of an increasing number of smaller an smaller simplices, where the interior of each simplex (Delaunay lattice) is assumed to be flat. The link lengths between any pair of vertices serve as independent variables, so that torsion turns out to be localized in the two dimensional hypersurfaces (dislocation triangle, or hinge) of the lattice. Assuming that a vector undergoes a dislocation in relation to its initial position as it is parallel transported along the perimeter of the dual lattice (Voronoi polygon), we obtain the discrete analogue of the teleparallel action, as well as the corresponding simplicial vacuum field equations.Comment: Latex, 10 pages, 2 eps figures, to appear in Class. Quant. Gra

Precision Determination of the Mass Function of Dark Matter Halos

The predicted mass function of dark matter halos is essential in connecting observed galaxy cluster counts and models of galaxy clustering to the properties of the primordial density field. We determine the mass function in the concordance $\Lambda$CDM cosmology, as well as its uncertainty, using sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a mass range of over five orders of magnitude. Using the nested volumes and single-halo tests, we find and correct for a systematic error in the friends-of-friends halo-finding algorithm. We find a fitting form and full error covariance for the mass function that successfully describes the simulations' mass function and is well-behaved outside the simulations' resolutions. Estimated forecasts of uncertainty in cosmological parameters from future cluster count surveys have negligible contribution from remaining statistical uncertainties in the central cosmology multiplicity function. There exists a potentially non-negligible cosmological dependence (non-universality) of the halo multiplicity function.Comment: 4 pages, 3 figures, submitted to ApJ

Torsion Degrees of Freedom in the Regge Calculus as Dislocations on the Simplicial Lattice

Using the notion of a general conical defect, the Regge Calculus is generalized by allowing for dislocations on the simplicial lattice in addition to the usual disclinations. Since disclinations and dislocations correspond to curvature and torsion singularities, respectively, the method we propose provides a natural way of discretizing gravitational theories with torsion degrees of freedom like the Einstein-Cartan theory. A discrete version of the Einstein-Cartan action is given and field equations are derived, demanding stationarity of the action with respect to the discrete variables of the theory

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. V. Rising X-ray Emission from an Off-Axis Jet

We report the discovery of rising X-ray emission from the binary neutron star (BNS) merger event GW170817. This is the first detection of X-ray emission from a gravitational-wave source. Observations acquired with the Chandra X-ray Observatory (CXO) at t~2.3 days post merger reveal no significant emission, with L_x<=3.2x10^38 erg/s (isotropic-equivalent). Continued monitoring revealed the presence of an X-ray source that brightened with time, reaching L_x\sim 9x10^39 erg/s at ~15.1 days post merger. We interpret these findings in the context of isotropic and collimated relativistic outflows (both on- and off-axis). We find that the broad-band X-ray to radio observations are consistent with emission from a relativistic jet with kinetic energy E_k~10^49-10^50 erg, viewed off-axis with theta_obs~ 20-40 deg. Our models favor a circumbinary density n~ 0.0001-0.01 cm-3, depending on the value of the microphysical parameter epsilon_B=10^{-4}-10^{-2}. A central-engine origin of the X-ray emission is unlikely. Future X-ray observations at $t\gtrsim 100$ days, when the target will be observable again with the CXO, will provide additional constraints to solve the model degeneracies and test our predictions. Our inferences on theta_obs are testable with gravitational wave information on GW170817 from Advanced LIGO/Virgo on the binary inclination.Comment: 7 Pages, 4 Figures, ApJL, In Press. Keywords: GW170817, LV

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VI. Radio Constraints on a Relativistic Jet and Predictions for Late-Time Emission from the Kilonova Ejecta

We present Very Large Array (VLA) and Atacama Large Millimeter/sub-millimeter Array ALMA radio observations of GW\,170817, the first Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo gravitational wave (GW) event from a binary neutron star merger and the first GW event with an electromagnetic (EM) counterpart. Our data include the first observations following the discovery of the optical transient at both the centimeter ($13.7$ hours post merger) and millimeter ($2.41$ days post merger) bands. We detect faint emission at 6 GHz at 19.47 and 39.23 days after the merger, but not in an earlier observation at 2.46 d. We do not detect cm/mm emission at the position of the optical counterpart at frequencies of 10-97.5 GHz at times ranging from 0.6 to 30 days post merger, ruling out an on-axis short gamma-ray burst (SGRB) for energies $\gtrsim 10^{48}$ erg. For fiducial SGRB parameters, our limits require an observer viewer angle of $\gtrsim 20^{\circ}$. The radio and X-ray data can be jointly explained as the afterglow emission from an SGRB with a jet energy of $\sim 10^{49}-10^{50}$ erg that exploded in a uniform density environment with $n\sim 10^{-4}-10^{-2}$ cm$^{-3}$, viewed at an angle of $\sim 20^{\circ}-40^{\circ}$ from the jet axis. Using the results of our light curve and spectral modeling, in conjunction with the inference of the circumbinary density, we predict the emergence of late-time radio emission from the deceleration of the kilonova (KN) ejecta on a timescale of $\sim 5-10$ years that will remain detectable for decades with next-generation radio facilities, making GW\,170817 a compelling target for long-term radio monitoring.Comment: 8 pages, 4 figures, 1 table. ApJL, in press. Keywords: GW170817, LV