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.

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.

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.

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.

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)].