
We review the reports of possible yearlong quasiperiodicities 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 quasiperiodic 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 quasiperiodic BL Lac objects is larger than that of
quasiperiodic quasars. If the BL Lac quasiperiodicities 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 quasiperiodicities 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 BondiHoyleLyttleton accretion rate, since the stellar
envelope is inhomogeneous, but may still be able to accrete at hypercritical
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 lowmass Xray binaries like Sco X1, 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 precommon envelope
electromagnetic observations.

We present $\texttt{ENIGMA}$, a time domain, inspiralmergerringdown
waveform model that describes nonspinning binary black holes systems that
evolve on moderately eccentric orbits. The inspiral evolution is described
using a consistent combination of postNewtonian theory, selfforce and black
hole perturbation theory. Assuming eccentric binaries that circularize prior to
coalescence, we smoothly match the eccentric inspiral with a standalone,
quasicircular merger, which is constructed using machine learning algorithms
that are trained with quasicircular numerical relativity waveforms. We show
that $\texttt{ENIGMA}$ reproduces with excellent accuracy the dynamics of
quasicircular 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 massratios 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, nonspinning binary black hole mergers. We use
$\texttt{ENIGMA}$ to show that GW150914, GW151226, GW170104, GW170814 and
GW170608 can be effectively recovered with spinning, quasicircular 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
quasicircular 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 higherorder
waveform multipoles in gravitational wave searches of eccentric binary black
hole mergers.

We have extended Cosmos++, a multidimensional 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 entropybased artificial viscosity and slope limiting
procedures for regularization of shocks. High order multistage forward Euler
and strong stability preserving RungeKutta 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
(cellbycell) 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
pulsedpowerdriven 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 fewcmthick stagnation layer, as observed in both
fastframing camera images and multichord interferometer data, consistent with
collisional shock formation [E. C. Merritt et al., Phys. Rev. Lett. {\bf 111},
085003 (2013)].