Search

Recent studies suggest that dark matter (DM) spikes around intermediate-mass black holes could cause observable dephasing in gravitational wave (GW) signals from Intermediate Mass Ratio Inspirals (IMRIs). Previous research primarily used non-self-consistent analytic methods to estimate the impact of DM spikes on eccentric IMRIs. Our study provides the first self-consistent treatment of this phenomenon using $N$-body simulations, incorporating Post-Newtonian effects up to the 2.5 order for accurate and robust results. Contrary to prior works, which posited that the cumulative effect of two-body encounters (dynamical friction; DF) is the primary mechanism for energy dissipation, we reveal that a three-body effect (slingshot mechanism) plays a more significant role in driving the binary system's energy loss and consequent orbital shrinkage. We find that binaries counter-rotating with respect to the DM spike merge faster, while co-rotating binaries merge slower, contrary to expectations from the DF theory. Using Fokker-Planck methods, we also assess the presence and detectability of spikes in realistic environments. When interacting with surrounding materials, DM spikes can have shallower slopes and lower densities than previously considered, leading to smaller signals and lower detection prospects via dephasing. Our results suggest that `deshifting' rather than dephasing might be a more optimistic signature, as it is more robust even in low-density environments., Comment: 21 pages, 20 figures. Accepted in MNRAS. The code used for simulations is available here: https://github.com/dipto4/falcon_dm

We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssim\Delta r<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Compared with the full AGN population, duals are characterized by a low $M_\text{BH}/M_*$ ratio, a high specific star-formation rate (sSFR) of $\sim 1\,\text{Gyr}^{-1}$, and a high Eddington ratio ($\sim 0.05$, double that of single AGN). The dual AGN are formed in major galaxy mergers (typically involving $M_\text{halo}<10^{13}\,M_\odot$), with BHs that have similar masses. At small separations (when their host galaxies are in the late phase of the merger) duals become $2\sim8$ times brighter (albeit more obscured) than at larger separations. $80\%$ of these bright, close duals merge in the simulation within $\sim500\,\text{Myrs}$. Notably, the initially less-massive BH in duals frequently becomes the brighter AGN during the galaxy merger. In offset AGN, the active BH is typically $\gtrsim 10$ times more massive than its non-active counterpart and than most BHs in duals. Offsets are predominantly formed in minor galaxy mergers with the active BH residing in the center of massive halos ($ M_\text{ halo}\sim 10^{13-14}\,M_\odot$). In these deep potentials, gas stripping is common and the secondary quickly deactivates. The stripping also leads to inefficient orbital decay amongst offsets, which stall at $\Delta r\sim5\,\text{kpc}$ for a few hundred Myrs., Comment: 17 Pages, 18 Figures; comments welcome!

We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to $1.5\;\text{ckpc}/h$. We calculate initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above $0.7$ for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below ${\sim 200\,\text{pc}}$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($\sim 500\,\text{Myrs}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries ($M_\text{BH}<2M_\text{seed}$). As a result, only $\lesssim20\%$ of seed MBH pairs merge at $z>3$ after considering both unresolved DF evolution and binary hardening. These $z>3$ seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $>10^9\,M_\odot$. With the higher initial eccentricity prediction from Astrid, we estimate an expected merger rate of $0.3-0.7$ per year from the $z>3$ MBH population. This is a factor of $\sim 7$ higher than the prediction using the circular orbit assumption. The LISA events are expected at a similar rate, and comprise $\gtrsim 60\%$ seed-seed mergers, $\sim 30\%$ involving only one seed-mass MBH, and $\sim 10\%$ mergers of non-seed MBHs., Comment: 17 pages, 13 Figures; comments are welcome

We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030$+$0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of $\gamma = 5$ and a reference frequency of $f_\text{yr} = 1 \text{yr}^{-1}$. Among the upper limits for eight general families of metric theories of gravity, we find the values of $A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}$ and $A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}$ for the family of metric spacetime theories that contain both TT and ST modes., Comment: 24 pages, 18 figures, 3 appendices. Please send any comments/questions to Nima Laal (laaln@oregonstate.edu)

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals., Comment: 13 pages, 4 figures. v2: updated to match published version

The second supernova that forms double-neutron-star systems is expected to occur in a progenitor that is ultra-stripped due to binary interactions. Thus, the secondary neutron star's mass as well as the post-supernova binary's orbital parameters will depend on the nature of the collapsing progenitor core. Since neutron stars are in the strong-gravity regime, their binding energy makes up a significant fraction of their total mass-energy budget. The second neutron star's binding energy may thus provide a unique insight as to whether its progenitor was a low-mass iron core or an oxygen-neon-magnesium core. I obtain relations for the baryonic mass and binding energy incorporating both a hadronic equation-of-state catalog as well as recent multi-messenger neutron-star observations. With these relations, I obtain the first tight constraints on the baryonic mass and binding energy of three neutron stars that are thought to have formed from an ultra-stripped progenitor. With these constraints, I test if each neutron star is consistent with forming from an ONeMg core that undergoes an electron-capture supernova. From these tests, I find that this scenario can be ruled out for one of three neutron stars. Neutron-star properties and the dense-matter equation of state can thus help distinguish binary formation scenarios., Comment: 7 pages, 5 figures

Common-envelope evolution is important in the formation of neutron star binaries within the isolated binary formation channel. As a neutron star inspirals within the envelope of a primary massive star, it accretes and spins up. Because neutron stars are in the strong-gravity regime, they have a substantial relativistic mass deficit, i.e., their gravitational mass is less than their baryonic mass. This effect causes some fraction of the accreted baryonic mass to convert into neutron star binding energy. The relativistic mass deficit also depends on the nuclear equation of state, since more compact neutron stars will have larger binding energies. We model the mass growth and spin-up of neutron stars inspiraling within common-envelope environments and quantify how different initial binary conditions and hadronic equations of state affect the post-common-envelope neutron star's mass and spin. From these models, we find that neutron star mass growth is suppressed by $\approx 15-30\%$. We also find that for a given amount of accreted baryonic mass, more compact neutron stars will spin-up faster while gaining less gravitational mass, and vice versa. This work demonstrates that a neutron star's strong gravity and nuclear microphysics plays a role in neutron-star-common-envelope evolution, in addition to the macroscopic astrophysics of the envelope. Strong gravity and the nuclear equation of state may thus affect both the population properties of neutron star binaries and the cosmic double neutron star merger rate., Comment: 11 pages, 6 figures, revised version accepted to ApJL

The gravitational-wave (GW) detection of GW190521 has provided new insights on the mass distribution of black holes and new constraints for astrophysical formation channels. With independent claims of GW190521 having significant pre-merger eccentricity, we investigate what this implies for GW190521-like binaries that form dynamically. The Laser Interferometer Space Antenna (LISA) will also be sensitive to GW190521-like binaries if they are circular from an isolated formation channel. We show, however, that GW190521-like binaries that form dynamically may skip the LISA band entirely. To this end, we simulate GW190521 analogues that dynamically form via post-Newtonian binary-single scattering. From these scattering experiments, we find that GW190521-like binaries may enter the LIGO-Virgo band with significant eccentricity as suggested by recent studies, though well below an eccentricity of $e_{\rm 10Hz} \lesssim 0.7$. Eccentric GW190521-like binaries further motivate the astrophysical science case for a decihertz GW observatory, such as the kilometer-scale version of the Midband Atomic Gravitational-wave Interferometric Sensor (MAGIS). Pre-merger observations of GW190521-like binaries with such a decihertz GW detector would be able to constrain the eccentricity of GW190521-like binaries to greater precision than with just LIGO-Virgo alone. These eccentricity constraints would also provide additional insights into the possible environments that GW190521-like binaries form in., Comment: 8 pages, 4 figures, revised version accepted to ApJL

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature., Comment: 25 pages, 14 figures, 5 tables, 3 appendices. Published in The Astrophysical Journal Letters. Please send any comments/questions to Joseph Simon (joe.simon@nanograv.org). Jupyter notebook tutorials and some MCMC chain files are available at https://github.com/nanograv/12p5yr_stochastic_analysis

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607$\pm$7 days) at 99.95% significance (with a global p-value of ~$10^{-3}$ accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67-002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ~20-yr time baseline spans ~4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302-102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna., Comment: Accepted to MNRAS, 17 pages, 9 figures

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg$^2$ overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-year long multi-color ($griz$) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to $r{\sim}$23.5 mag) powered by less massive black holes (with masses $\gtrsim10^{8.5}M_{\odot}$) at high redshift for the first time. We find five candidates with significant (at $>$99.74% single-frequency significance in at least two bands with a global p-value of $\sim$7$\times10^{-4}$--3$\times10^{-3}$ accounting for the look-elsewhere effect) periodicity with observed periods of $\sim$3--5 years (i.e., 1--2 years in rest frame) having $\sim$4--6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of ${\sim}0.8^{+0.5}_{-0.3}$% or ${\sim}1.1^{+0.7}_{-0.5}$ quasar per deg$^2$. Our detection rate is 4--80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. We discuss implications on the future direct detection of low-frequency gravitational waves. Continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins., Comment: 21 pages, 8 figures, MNRAS accepted

The Neutron Star Interior Composition Explorer (NICER) recently measured the mass and equatorial radius of the isolated neutron star PSR J0030+0451. We use these measurements to infer the moment of inertia, the quadrupole moment, and the surface eccentricity of an isolated neutron star for the first time, using relations between these quantities that are insensitive to the unknown equation of state of supranuclear matter. We also use these results to forecast the moment of inertia of neutron star A in the double pulsar binary J0737-3039, a quantity anticipated to be directly measured in the coming decade with radio observations. Combining this information with the measurement of the tidal Love number with LIGO/Virgo observations, we propose and implement the first theory-agnostic and equation-of-state-insensitive test of general relativity. Specializing these constraints to a particular modified theory, we find that consistency with general relativity places the most stringent constraint on gravitational parity violation to date, surpassing all other previously reported bounds by seven orders of magnitude and opens the path for future a test of general relativity with multimessenger neutron star observations., Comment: 12 pages, 7 figures, 5 tables. (v2) Major revision; expanded supplemental material. (v3) Matches version published in Phys. Rev. Lett

Some fraction of compact binaries that merge within a Hubble time may have formed from two massive stars in isolation. For this isolated-binary formation channel, binaries need to survive two supernova (SN) explosions in addition to surviving common-envelope evolution. For the SN explosions, both the mass loss and natal kicks change the orbital characteristics, producing either a bound or unbound binary. We show that gravitational waves (GWs) may be produced not only from the core-collapse SN process, but also from the SN mass loss and SN natal kick during the pre-SN to post-SN binary transition. We model the dynamical evolution of a binary at the time of the second SN explosion with an equation of motion that accounts for the finite timescales of the SN mass loss and the SN natal kick. From the dynamical evolution of the binary, we calculate the GW burst signals associated with the SN natal kicks. We find that such GW bursts may be of interest to future mid-band GW detectors like DECIGO. We also find that the energy radiated away from the GWs emitted due to the SN mass loss and natal kick may be a significant fraction, ${\gtrsim}10\%$, of the post-SN binary's orbital energy. For unbound post-SN binaries, the energy radiated away in GWs tends to be higher than that of bound binaries., Comment: 7 pages, 4 figures, submitted to MNRAS, in revision

In the quadrupole approximation of General Relativity in the weak-field limit, a time-varying quadrupole moment generates gravitational radiation. Binary orbits are one of the main mechanisms for producing gravitational waves and are the main sources and backgrounds for gravitational-wave detectors across the multi-band spectrum. In this Paper, we introduce additional contributions to the gravitational radiation from close binaries that arise from time-varying masses along with those produced by orbital motion. We derive phase-dependent formulae for these effects in the quadrupolar limit for binary point masses, which reduce to the formulae that Peters and Mathews (1963) derived when the mass of each component is taken to be constant. We show that gravitational radiation from mass variation can be orders of magnitude greater than that of orbital motion., Comment: 6 pages, 1 figure, accepted

We propose a novel methodology to select host galaxy candidates of future pulsar timing array (PTA) detections of resolved gravitational waves (GWs) from massive black hole binaries (MBHBs). The method exploits the physical dependence of the GW amplitude on the MBHB chirp mass and distance to the observer, together with empirical MBH mass-host galaxy correlations, to rank potential host galaxies in the mass-redshift plane. This is coupled to a null-stream based likelihood evaluation of the GW amplitude and sky position in a Bayesian framework that assigns to each galaxy a probability of hosting the MBHB generating the GW signal. We test our algorithm on a set of realistic simulations coupling the likely properties of the first PTA resolved GW signal to synthetic all-sky galaxy maps. For a foreseeable PTA sky-localization precision of 100 squared degrees, we find that the GW source is hosted with 50% (90%) probability within a restricted number of <50 (<500) potential hosts. These figures are orders of magnitude smaller than the total number of galaxies within the PTA sky error-box, enabling extensive electromagnetic follow-up campaigns on a limited number of targets., Comment: 12 pages, 10 figures

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs., Comment: Accepted to The Astronomy and Astrophysics Review (A&ARv)

Blazars are a sub-population of quasars whose jets are nearly aligned with the line-of-sight, which tend to exhibit multi-wavelength variability on a variety of timescales. Quasi-periodic variability on year-like timescales has been detected in a number of bright sources, and has been connected to the orbital motion of a putative massive black hole binary. If this were indeed the case, those blazar binaries would contribute to the nanohertz gravitational-wave stochastic background. We test the binary hypothesis for the blazar population observed by the \textit{Fermi} Gamma-Ray Space Telescope, which consists of BL Lacertae objects and flat-spectrum radio quasars. Using mock populations informed by the luminosity functions for BL Lacertae objects and flat-spectrum radio quasars with redshifts $z \le 2$, we calculate the expected gravitational wave background and compare it to recent pulsar timing array upper limits. The two are consistent only if a fraction $\lesssim 10^{-3}$ of blazars hosts a binary with orbital periods $<5$ years. We therefore conclude that binarity cannot significantly explain year-like quasi-periodicity in blazars., Comment: 5 pages, 4 figures, accepted by MNRAS Letters

We review the reports of possible year-long quasi-periodicities of BL Lac objects in the $\gamma$-ray and optical bands, and present a homogeneous time analysis of the light curves of PKS2155$-$304, PG1553+113, and BL Lac. Based on results from a survey covering the entire Fermi $\gamma$-ray sky we have estimated the fraction of possible quasi-periodic BL Lac objects. We compared the cyclical behaviour in BL Lac objects with that derived from the search of possible optical periodicities in quasars, and find that at z$\lesssim$1 the cosmic density of quasi-periodic BL Lac objects is larger than that of quasi-periodic quasars. If the BL Lac quasi-periodicities were due to a supermassive binary black hole (SBBH) scenario, there could be a tension with the upper limits on the gravitational wave background measured by the pulsar timing array. The argument clearly indicates the difficulties of generally associating quasi-periodicities of BL Lac objects with SBBHs., Comment: In publication on A&A, 6 pages, 4 figure (11 plots). Minor corrections added

We present $\texttt{ENIGMA}$, a time domain, inspiral-merger-ringdown waveform model that describes non-spinning binary black holes systems that evolve on moderately eccentric orbits. The inspiral evolution is described using a consistent combination of post-Newtonian theory, self-force and black hole perturbation theory. Assuming eccentric binaries that circularize prior to coalescence, we smoothly match the eccentric inspiral with a stand-alone, quasi-circular merger, which is constructed using machine learning algorithms that are trained with quasi-circular numerical relativity waveforms. We show that $\texttt{ENIGMA}$ reproduces with excellent accuracy the dynamics of quasi-circular compact binaries. We validate $\texttt{ENIGMA}$ using a set of $\texttt{Einstein Toolkit}$ eccentric numerical relativity waveforms, which describe eccentric binary black hole mergers with mass-ratios between $1 \leq q \leq 5.5$, and eccentricities $e_0 \lesssim 0.2$ ten orbits before merger. We use this model to explore in detail the physics that can be extracted with moderately eccentric, non-spinning binary black hole mergers. We use $\texttt{ENIGMA}$ to show that GW150914, GW151226, GW170104, GW170814 and GW170608 can be effectively recovered with spinning, quasi-circular templates if the eccentricity of these events at a gravitational wave frequency of 10Hz satisfies $e_0\leq \{0.175,\, 0.125,\,0.175,\,0.175,\, 0.125\}$, respectively. We show that if these systems have eccentricities $e_0\sim 0.1$ at a gravitational wave frequency of 10Hz, they can be misclassified as quasi-circular binaries due to parameter space degeneracies between eccentricity and spin corrections. Using our catalog of eccentric numerical relativity simulations, we discuss the importance of including higher-order waveform multipoles in gravitational wave searches of eccentric binary black hole mergers., Comment: 19 pages, 10 figures, 1 Appendix. v2: we use numerical relativity simulations to quantify the importance of including higher-order waveform multipoles for the detection of eccentric binary black hole mergers, references added. Accepted to Phys. Rev. D

We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for regularization of shocks. High order multi-stage forward Euler and strong stability preserving Runge-Kutta time integration options complement high order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, though we note an equivalent capability currently also exists in CosmosDG for Newtonian systems., Comment: 21 pages, 14 figures, submitted to The Astrophysical Journal

The common envelope phase is a likely formation channel for close binary systems containing compact objects. Neutron stars in common envelopes accrete at a fraction of the Bondi-Hoyle-Lyttleton accretion rate, since the stellar envelope is inhomogeneous, but may still be able to accrete at hyper-critical rates (though not enough to become black holes). We show that common envelope systems consisting of a neutron star with a massive primary may be gravitational wave sources detectable in the Advanced LIGO band as far away as the Magellanic Clouds. To characterize their evolution, we perform orbital integrations using 1D models of $12 M_\odot$ and $20 M_\odot$ primaries, considering the effects of density gradient on the accretion onto the NS and spin evolution. From the range of possible accretion rates relevant to common envelope evolution, we find that these systems may be louder gravitational wave sources than low-mass X-ray binaries like Sco X-1, which are currently the target of directed searches for continuous GWs. We also find that their strain amplitude signal may allow for novel constraints on the orbital separation and inspiral timescale in common envelopes when combined with pre-common envelope electromagnetic observations., Comment: 20 pages, 9 figures, accepted to ApJ

We report spatially resolved measurements of the oblique merging of two supersonic laboratory plasma jets. The jets are formed and launched by pulsed-power-driven railguns using injected argon, and have electron density $\sim 10^{14}$ cm$^{-3}$, electron temperature $\approx 1.4$ eV, ionization fraction near unity, and velocity $\approx 40$ km/s just prior to merging. The jet merging produces a few-cm-thick stagnation layer, as observed in both fast-framing camera images and multi-chord interferometer data, consistent with collisional shock formation [E. C. Merritt et al., Phys. Rev. Lett. {\bf 111}, 085003 (2013)]., Comment: 11 pages, 9 figures, accepted by Physics of Plasmas (to be published in 2014 May issue)