We present the results from the timing and spectral study of Mrk\,421 based mainly on the \emph{Swift} data in the X-ray energy range obtained during the time interval 2015\,December--2018\,April. The most extreme X-ray flaring activity on the long-term, daily and intraday timescales was observed during the 2-month period which started in 2017\,December when the 0.3--10\,keV flux exceeded a level of 5$\times$10$^{-9}$erg\,cm$^{-2}$s$^{-1}$, recorded only twice previously. While the TeV-band and X-ray variabilities mostly were correlated, the source often varied in a complex manner in the MeV--GeV and radio--UV energy ranges, indicating that the multifrequency emission of Mrk\,421 could not be always generated in a single zone. The longer-term flares at X-rays and $\gamma$-rays showed a lognormal character, possibly indicating a variability imprint of the accretion disk onto the jet. A vast majority of the 0.3--10\,keV spectra were consistent with the log-parabolic model, showing relatively low spectral curvature and correlations between the different spectral parameters, predicted in the case of the first and second-order Fermi processes. The position of the synchrotron spectral energy distribution (SED) peak showed an extreme variability on diverse timescales between the energies $E_{\rm p}$$<$0.1\,keV and $E_{\rm p}$$>$15\,keV, with 15\% of the spectra peaking at hard X-rays and was related to the peak height as $S_{\rm p}$$\varpropto$$E^{\alpha}_{\rm p}$ with $\alpha$$\sim$0.6, which is expected for the transition from Kraichnan-type turbulence into the \textquotedblleft hard-sphere\textquotedblright~ one. The 0.3--300\,GeV spectra showed the features of the hadronic contribution, jet-star interaction and upscatter in the Klein-Nishina regime in different time intervals.

We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

The University of Michigan 26-m paraboloid was dedicated to obtaining linear polarization and total flux density observations of blazars from the mid-1960s until June 2012 providing an unprecedented record tracking centimeter-band variability over decades at 14.5, 8.0, and 4.8 GHz for both targeted objects and members of flux-limited samples. In the mid-1970s through the mid-1980s, and during the last decade of the program, observations were additionally obtained of circular polarization for a small sample of radio-bright (S>5Jy), active sources. Key program results include evidence supporting class-dependent differences in the magnetic field geometry of BL Lac and QSO jets, identification of linear polarization changes temporally associated with flux outbursts supporting a shock-in-jet scenario, and determination of the spectral evolution of the Stokes V amplitude and polarity for testing proposed models. Recent radiative transfer modeling during large flares supports a jet scenario with a kinetically-dominated, relativistic flow at parsec scales with embedded turbulent magnetic fields and dynamically-weak ordered components which may be helical; the circular polarization observations are consistent with linear-to-circular mode conversion within this turbulent jet environment.

Since the mid-1980s the shock-in-jet model has been the preferred paradigm to explain radio-band flaring in blazar jets. We describe our radiative transfer model incorporating relativistically-propagating shocks, and illustrate how the 4.8, 8, and 14.5 GHz linear polarization and total flux density data from the University of Michigan monitoring program, in combination with the model, constrain jet flow conditions and shock attributes. Results from strong Fermi-era flares in 4 blazars with widely-ranging properties are presented. Additionally, to investigate jet evolution on decadal time scales we analyze 3 outbursts in OT 081 spanning nearly 3 decades and find intrinsic changes attributable to flow changes at a common spatial location, or, alternatively, to a change in the jet segment viewed. The model's success in reproducing these data supports a scenario in which relativistic shocks compress a plasma with an embedded passive, initially-turbulent magnetic field, with additional ordered magnetic field components, one of which may be helical.

As part of a program to identify the physical conditions in the jets of gamma-ray-flaring blazars detected by Fermi, including the role of shocks in the production of high-energy flaring, we obtained 4 years of 3-frequency, centimeter-band total flux density and linear polarization monitoring observations of the radio-bright blazar S5 0716+714 with the University of Michigan 26-m paraboloid. Light curves constructed from these data exhibit a series of rapid, high-amplitude, centimeter-band total flux density outbursts, and changes in the linear polarization consistent with the passage of shocks during the gamma-ray flaring. The observed spectral evolution of the radio-band flares, in combination with radiative transfer simulations incorporating propagating shocks, was used to constrain the shock and jet flow conditions in the parsec-scale regions of the jet. Eight forward-moving, transverse shocks with unusually-strong shock compression factors, a very fast Lorentz factor of the shocks of 77, a bulk Lorentz factor of the flow of 20, a viewing angle of 12 degrees, and an intrinsic opening angle of the radio jet of 5.2 degrees were identified.

We analyze the shock-in-jet models for the gamma-ray flaring blazars 0420-014, OJ 287, and 1156+295 presented in Aller et al. (2014, Paper I), quantifying how well the modeling constrains internal properties of the flow (low energy spectral cutoff, partition between random and ordered magnetic field), the flow dynamics (quiescent flow speed and orientation), and the number and strength of the shocks responsible for radio-band flaring. We conclude that well-sampled, multifrequency polarized flux light curves are crucial for defining source properties. We argue for few, if any, low energy particles in these flows, suggesting no entrainment and efficient energization of jet material, and for approximate energy equipartition between the random and ordered magnetic field components, suggesting that ordered field is built by non-trivial dynamo action from the random component, or that the latter arises from a jet instability that preserves the larger-scale, ordered flow. We present evidence that the difference between orphan radio-band (no gamma-ray counterpart) and non-orphan flares is due to more complex shock interactions in the latter case.

This paper examines the radiative transfer of synchrotron radiation in the presence of a magnetic field configuration resulting from the compression of a highly disordered magnetic field. It is shown that, provided Faraday rotation and circular polarization can be neglected, the radiative transfer equations for synchrotron radiation separate for this configuration, and the intensities and polarization values for sources that are uniform on large scales can be found straightforwardly in the case where opacity is significant. Although the emission and absorption coefficients must, in general, be obtained numerically, the process is much simpler than a full numerical solution to the transfer equations. Some illustrative results are given and an interesting effect, whereby the polarization increases while the magnetic field distribution becomes less strongly confined to the plane of compression, is discussed. The results are of importance for the interpretation of polarization near the edges of lobes in radio galaxies and of bright features in the parsec--scale jets of AGN, where such magnetic field configurations are believed to exist. (Note: The original ApJ version of this paper contained two errata, which are corrected in this version of the paper. First, Fig. 2 was plotted as a mirror image of the correct version, reflected about the line $\log_{10}(\nu/ \nu_0)=0$. Second, as stated, Equation A16 is valid for $\gamma=2$, not $\gamma=3$.)

To investigate parsec-scale jet flow conditions during GeV gamma-ray flares detected by the Fermi Large Angle Telescope, we obtained centimeter-band total flux density and linear polarization monitoring observations from 2009.5 through 2012.5 with the 26-meter Michigan radio telescope for a sample of core-dominated blazars. We use these data to constrain radiative transfer simulations incorporating propagating shocks oriented at an arbitrary angle to the flow direction in order to set limits on the jet flow and shock parameters during flares temporally associated with gamma-ray flares in 0420-014, OJ 287, and 1156+295; these AGN exhibited the expected signature of shocks in the linear polarization data. Both the number of shocks comprising an individual radio outburst (3-4) and the range of the compression ratios of the individual shocks (0.5-0.8) are similar in all three sources; the shocks are found to be forward-moving with respect to the flow. While simulations incorporating transverse shocks provide good fits for 0420-014 and 1156+295, oblique shocks are required for modeling the OJ~287 outburst, and an unusually-low value of the low energy cutoff of the radiating particles' energy distribution is also identified. Our derived viewing angles and shock speeds are consistent with independent VLBA results. While a random component dominates the jet magnetic field, as evidenced by the low fractional linear polarization, to reproduce the observed spectral character requires that a significant fraction of the magnetic field energy is in an ordered axial component. Both straight and low pitch angle helical field lines are viable scenarios.

As part of a program to investigate jet flow conditions during GeV gamma-ray flares detected by Fermi, we are using UMRAO multi-frequency, centimeter-band total flux density and linear polarization monitoring observations to constrain radiative transfer models incorporating propagating shocks orientated at an arbitrary angle to the flow direction. We describe the characteristics of the model, illustrate how the data are used to constrain the models, and present results for three program sources with diverse characteristics: PKS 0420-01, OJ 287, and 1156+295. The modeling of the observed spectral behavior yields information on the sense, strength and orientation of the shocks producing the radio-band flaring; on the energy distribution of the radiating particles; and on the observer's viewing angle with respect to the jet independent of VLBI data. We present evidence that, while a random component dominates the jet magnetic field, a distinguishing feature of those radio events with an associated gamma-ray flare is the presence of a weak but non-negligible ordered magnetic field component along the jet axis.

We have used single-dish centimeter-band, multi-frequency linear polarization and total flux density monitoring observations obtained by UMRAO at 14.5, 8, and 4.8 GHz during gamma-ray flaring to constrain propagating shock-in-jet models; this procedure allows us to probe jet conditions at or near the presumed parsec-scale gamma-ray emission site within the relativistic jet. Results are presented for the blazar 0420-014 during the gamma-ray flare detected by the LAT which peaked in late January 2010. Our work identifies the passage of three forward-moving transverse shocks during the radio-band flaring, a shock Lorentz factor of 12, and an observer's viewing angle of 2 degrees and sets limits on the energy distribution of the radiating particles.

We present centimeter-band total flux density and linear polarization light curves illustrating the signature of shocks during radio band outbursts associated in time with gamma-ray flares detected by the Fermi LAT. The general characteristics of the spectral evolution during these events is well-explained by new radiative transfer simulations incorporating propagating oblique shocks and assuming an initially turbulent magnetic field. This finding supports the idea that oblique shocks in the jet are a viable explanation for activity from the radio to the gamma-ray band in at least some gamma-ray flares.

The `shock in jet' model for cm-waveband blazar variability is revisited, allowing for arbitrary shock orientation with respect to the jet flow direction, and both random and ordered magnetic field. It is shown that oblique shocks can explain events with swings in polarization position angle much less than the 90 deg. associated with transverse structures, while retaining the general characteristics of outbursts, including spectral behavior and level of peak percentage polarization. Models dominated by a force-free, minimum energy magnetic field configuration (essentially helical) display a shallow rise in percentage polarization and frequency dependent swing in polarization position angle not in agreement with the results of single-dish monitoring observations, implying that the field is predominantly random in the quiescent state. Outbursts well-explained by the `shock in jet' model are present during gamma-ray flaring in several sources, supporting the idea that shock events are responsible for activity from the radio to gamma-ray bands.

The properties of blazar variability in the radio band are studied using the unique combination of temporal resolution from single dish monitoring and spatial resolution from VLBA imaging; such measurements, now available in all four Stokes parameters, together with theoretical simulations, identify the origin of radio band variability and probe the characteristics of the radio jet where the broadband blazar emission originates. Outbursts in total flux density and linear polarization in the optical-to-radio bands are attributed to shocks propagating within the jet spine, in part based on limited modeling invoking transverse shocks; new radiative transfer simulations allowing for shocks at arbitrary angle to the flow direction confirm this picture by reproducing the observed centimeter-band variations observed more generally, and are of current interest since these shocks may play a role in the gamma-ray flaring detected by Fermi. Recent UMRAO multifrequency Stokes V studies of bright blazars identify the spectral variability properties of circular polarization for the first time and demonstrate that polarity flips are relatively common. All-Stokes data are consistent with the production of circular polarization by linear-to-circular mode conversion in a region that is at least partially self-absorbed. Detailed analysis of single-epoch, multifrequency, all-Stokes VLBA observations of 3C 279 support this physical picture and are best explained by emission from an electron-proton plasma.

In the 1990s a comparison of sparse EGRET measurements with single-dish flux density monitoring from the Metsahovi and UMRAO programs established a temporal connection between the onset of flaring at radio band and the occurrence of gamma-ray activity. Correlations between the emergence of new VLBI components from the core, flares in linearly polarized radio flux, and gamma-ray activity in bright EGRET-detected blazars supported a picture in which the gamma-ray and the radio band emission arises in the same shocked region of the jet, with the high energy emission produced via inverse Compton scattering by the synchrotron-emitting electrons in the jet. Quantitative tests of this scenario, however, were hampered by insufficient temporal sampling of the data and the simple nature of the models adopted. The extensive data from Fermi coupled with the wealth of well-sampled radio band data from old as well as new programs such as the F-GAMMA project now permit statistical studies for large numbers of sources, including weak HBLs, and detailed analyses of individual highly-active class members. I summarize progress in understanding the origin of the gamma-ray emission using these new measurements. I focus on three areas: attempts to isolate the physical site of the high energy emission using time delay information; investigation of the emission process using the characteristics of the variability; and quantitative tests of the shock model picture using high-time-sampled multifrequency linear polarization data, VLBP imaging, and new models of propagating oblique relativistic shocks incorporating detailed radiative transfer calculations.

We describe and present initial results from a Fermi cycle 2 program designed to monitor the behavior of the centimeter-band linear polarization and total flux density emitted by gamma-ray-bright blazars during flaring. The goal of the program is to identify changes in the magnetic field structure in the radio jet associated with gamma-ray flaring and ultimately to test whether gamma-ray flaring is associated with the onset of shocks in the radio jet. Light curves illustrating radio band variability patterns are shown for sample program sources.

In this letter, we describe results of new high-resolution axisymmetric relativistic MHD simulations of Pulsar Wind Nebulae. The simulations reveal strong breakdown of the equatorial symmetry and highly variable structure of the pulsar wind termination shock. The synthetic synchrotron maps, constructed using a new more accurate approach, show striking similarity with the well known images of the Crab Nebula obtained by Chandra, and the Hubble Space Telescope. In addition to the jet-torus structure, these maps reproduce the Crab's famous moving wisps whose speed and rateof production agree with the observations. The variability is then analyzed using various statistical methods, including the method of structure function and wavelet transform. The results point towards the quasi-periodic behaviour with the periods of 1.5-3yr and MHD turbulence on scales below 1yr. The full account of this study will be presented in a follow up paper.

We have undertaken the simulation of hydrodynamic flows with bulk Lorentz factors in the range 10^2--10^6. We discuss the application of an existing relativistic, hydrodynamic primitive-variable recovery algorithm to a study of pulsar winds, and, in particular, the refinement made to admit such ultra-relativistic flows. We show that an iterative quartic root finder breaks down for Lorentz factors above 10^2 and employ an analytic root finder as a solution. We find that the former, which is known to be robust for Lorentz factors up to at least 50, offers a 24% speed advantage. We demonstrate the existence of a simple diagnostic allowing for a hybrid primitives recovery algorithm that includes an automatic, real-time toggle between the iterative and analytical methods. We further determine the accuracy of the iterative and hybrid algorithms for a comprehensive selection of input parameters and demonstrate the latter's capability to elucidate the internal structure of ultra-relativistic plasmas. In particular, we discuss simulations showing that the interaction of a light, ultra-relativistic pulsar wind with a slow, dense ambient medium can give rise to asymmetry reminiscent of the Guitar nebula leading to the formation of a relativistic backflow harboring a series of internal shockwaves. The shockwaves provide thermalized energy that is available for the continued inflation of the PWN bubble. In turn, the bubble enhances the asymmetry, thereby providing positive feedback to the backflow.

Quasi-periodic signals have yielded important constraints on the masses of black holes in galactic X-ray binaries, and here we extend this to active galactic nuclei (AGNs). We employ a wavelet technique to analyze 19 observations of 10 AGNs obtained with the XMM-Newton EPIC pn camera. We report the detection of a candidate 3.3 ks quasi period in 3C 273. If this period represents an orbital timescale originating near a last stable orbit of 3Rs, it implies a central black hole mass of 7.3x10^6 Msun. For a maximally rotating black hole with a last stable orbit of 0.6Rs, a central black hole mass of 8.1x10^7 Msun is implied. Both of these estimates are substantially lower than previous reverberation-mapping results, which place the central black hole mass of 3C 273 at about 2.35 x10^8 Msun. Assuming that this reverberation mass is correct, the X-ray quasi period would be caused by a higher order oscillatory mode of the accretion disk.

In this letter we present the results of a wavelet analysis of the radio light curve of the BL Lac Object PKS B0048-097 from the University of Michigan Radio Astronomy Observatory monitoring program at 8GHz during twenty-five years, from 1979 to 2004. The results show a remarkable periodicity of 450-470 days in the early 1980s that changed to a ~585 day periodicity in the late 1980s to early 1990s. A less pronounced ~400 day periodicity is found after \~1995. Very-long-baseline interferometry imaging at 15GHz shows dramatic structural changes in the usually unresolved source between two epochs, 1995.57 and 2002.38. The pronounced northward directed jet seen in the 2002 image differs by more than 90 degrees in direction from the source structure found in earlier epochs. These findings make PKS B0048-097 a primary target for multi-wavelength observations and intensive radio monitoring to decipher the blazar-variability phenomenon.

Evidence that the magnetic fields of extragalactic jets have a significant fraction of their energy in a random component is briefly summarized, and a detailed model of evolving, jet polarization structures is constructed, based on this picture. The evolving magnetic field structure of an oblique shock complex that forms in a relativistic jet simulation is explored by using velocity data from the hydrodynamical simulation to advect an initially random magnetic field distribution. Radiative transfer calculations reveal that emission from a propagating region of magnetic field, `ordered' by the shock, and lying approximately transverse to the flow direction, merges with that from an evolving sheared region at the flow periphery. If such a flow were barely resolved, observation would suggest evolution from a somewhat oblique, to a more longitudinal, magnetic field structure with respect to the flow axis, while higher resolution observations would infer a component following a non-linear trajectory, and with a magnetic field orientation that rotates during evolution. This result highlights the ambiguity in interpreting VLBP data, and illustrates the importance of simulations in providing a framework for proper interpretation of such data.

We introduce an algorithm for applying a cross-wavelet transform to analysis of quasiperiodic variations in a time-series, and introduce significance tests for the technique. We apply a continuous wavelet transform and the cross-wavelet algorithm to the Pearson-Readhead VLBI survey sources using data obtained from the University of Michigan 26-m parabloid at observing frequencies of 14.5, 8.0, and 4.8 GHz. Thirty of the sixty-two sources were chosen to have sufficient data for analysis, having at least 100 data points for a given time-series. Of these thirty sources, a little more than half exhibited evidence for quasiperiodic behavior in at least one observing frequency, with a mean characteristic period of 2.4 yr and standard deviation of 1.3 yr. We find that out of the thirty sources, there were about four time scales for every ten time series, and about half of those sources showing quasiperiodic behavior repeated the behavior in at least one other observing frequency.

Using UMRAO centimeter-band total flux density and linear polarization monitoring observations of the complete Pearson-Readhead extragalactic source sample obtained between August 1984 and March 2001, we identify the range of variability in extragalactic objects as functions of optical and radio morphological classification, and relate total flux density variations to structural changes in published coeval VLBI maps in selected objects. We have identified preferred orientations of the electric vector of the polarized emission at 14.5 and 4.8 GHz in roughly half of the objects, and compared these with orientations of the flow direction indicated by VLBI morphology. When comparing the distributions of the orientation offsets for the BL Lacs and for the QSOs, we find differences in both range and mean value, in support of intrinsic class differences. In the shock-in-jet scenario, we attribute this to the allowed range of obliquities of shocks developing in the flow relative to the flow direction. The fact that we find longterm stability in EVPA over many events implies that a dominant magnetic field orientation persists; in the core-dominated objects, with small contribution from the underlying quiescent jet, this plausibly suggests that the magnetic field has a longterm memory, with subsequent shock events exhibiting similar EVPA orientation, or, alternatively, the presence of a standing shock in the core. We have looked for systematic, monotonic changes in EVPA which might be expected in the emission from a precessing jet, a model currently invoked for some AGNs; none were identified. Further, we carried out a Scargle periodogram analysis of the total flux density observations, but found no strong evidence for periodicity in any of the sample sources.

We interpret the linear polarization structures observed in extragalactic radio sources, even those oriented at oblique angles to the jet flows, to be due to oblique, relativistic shock fronts in the emitting regions. Many sources exhibit indications of such oblique structures, and the goal of this investigation is to test this hypothesis quantitatively. A selected group of ten highly-variable extragalactic sources were observed with the VLBA at 15 and 43 GHz, on nine epochs spanning a 30-month period; five of these objects were also observed at 8.0 and 22 GHz. The integrated total flux densities and linear polarizations of the selected objects were also observed several times a month at 4.8, 8.0 and 14.5 GHz with the University of Michigan 26-meter telescope. All objects exhibited variability with several exhibiting more than one independent outburst during the period. We show the evolution of the polarized components with time and discuss the relativistic shock parameters required to match the observed polarization structures. Even cases where the magnetic field is apparently oriented along the jet flow can be fit by oblique shock models when relativistic aberration effects are included.

Monitoring results are presented for the Stangellini 1 Jy GHz-peaked-spectrum source sample, illustrating that several members exhibit variability in total flux and/or linear polarization over timescales of order a decade. The variability occurs while the spectrum, based on the integrated fluxes, remains steep and characteristic of a transparent source. Total flux variability is unexpected in view of recent VLBI observations indicating no or hidden cores in several members. However, both the variability, and the detection of circular polarization in one class member, argue for the presence of opacity in some portion of the radio jet.

Relativistic jets in Galactic superluminals and extragalactic AGN may be surrounded by a wind near to the central engine. Theoretical analysis and numerical simulation reveal considerable stabilization of relativistic jet flow by a wind to helical and higher order asymmetric modes of jet distortion. When velocities are measured in the source (inlet) frame, reduction in the absolute velocity difference between jet and wind, $\Delta v = v_{jet} - v_{wnd}$, provides stabilization in addition to stabilization provided by a high jet Lorentz factor, but a high Lorentz factor wind is not needed to stabilize a high Lorentz factor jet. However, the fundamental pinch mode is not similarly affected and knots with spacing a few times the jet radius are anticipated to develop in such flows. Thus, we identify a mechanism that can suppress large scale asymmetric structures while allowing axisymmetric structures to develop. Relativistic jets surrounded by outflowing winds will be more stable than if surrounded by a stationary or backflowing external medium. Knotty structures along a straight jet like that in 3C 175 could be triggered by pinching of an initially low Mach number jet surrounded by a suitable wind. As the jet enters the radio lobe, suppression of any surrounding outflow or backflow associated with the high pressure lobe triggers exponential growth of helical twisting.

We describe a new numerical 3D relativistic hydrodynamical code, the results of validation tests, and a comparison with earlier, 2D studies. The 3D code has been used to study the deflection and precession of relativistic flows. We find that even quite fast jets (gamma~10) can be significantly influenced by impinging on an oblique density gradient, exhibiting a rotation of the Mach disk in the jet's head. The flow is bent via a potentially strong, oblique internal shock that arises due to asymmetric perturbation of the flow by its cocoon. In extreme cases this cocoon can form a marginally relativistic flow orthogonal to the jet, leading to large scale dynamics quite unlike that normally associated with astrophysical jets. Exploration of a gamma=5 flow subject to a large amplitude precession (semi-angle 11.25dg) shows that it retains its integrity, with modest reduction in Lorentz factor and momentum flux, for almost 50 jet-radii, but thereafter, the collimated flow is disrupted. The flow is approximately ballistic, with velocity vectors not aligned with the local jet `wall'. We consider simple estimators of the flow emissivity in each case and conclude that a) while the oblique internal shocks which mediate a small change in the direction of the deflected flows have little impact on the global dynamics, significantly enhanced flow emission (by a factor of 2-3) may be associated with such regions; and b) the convolution of rest frame emissivity and Doppler boost in the case of the precessed jet invariably leads to a core-jet-like structure, but that intensity fluctuations in the jet cannot be uniquely associated with either change in internal conditions or Doppler boost alone, but in general are a combination of both factors.

We present analysis of the flux and polarization variability of parsec scale radio jets from a dual-frequency, six-epoch, VLBA polarization experiment monitoring 12 blazars. The observations were made at 15 and 22 GHz at bimonthly intervals over 1996. Here we analyze the flux, fractional polarization, and polarization position angle behavior of core regions and jet features, considering both the linear trends of these quantities with time and more rapid fluctuations about the linear trends. The dual frequency nature of the observations allows us to examine spectral evolution, to separate Faraday effects from changes in magnetic field order, and also to deduce empirical estimates for the uncertainties in measuring properties of VLBI jet features (abridged).

Three dimensional numerical simulations of the response of a Lorentz factor 2.5 relativistic jet to precession at three different frequencies have been performed. Low, moderate and high precession frequencies have been chosen relative to the maximally unstable frequency predicted by a Kelvin-Helmholtz stability analysis. Transverse motion and velocity decreases as the precession frequency increases. Although helical displacement of the jet decreases in amplitude as the precession frequency increases, a helical shock is generated in the medium external to the jet at all precession frequencies. Complex pressure and velocity structure inside the jet is shown to be produced by a combination of the helical surface and first body modes predicted by a normal mode analysis of the relativistic hydrodynamic equations. The surface and first body mode have different wave speed and wavelength, are launched in phase by the periodic precession, and exhibit beat patterns in synthetic emission images. Wave (pattern) speeds range from $0.41c$ to $0.86c$ but beat patterns remain stationary. Thus, we find a mechanism that can produce differentially moving and stationary features in the jet.

We perform fully relativistic hydrodynamic simulations of the deceleration and lateral expansion of a relativistic jet as it expands into an ambient medium. The hydrodynamic calculations use a 2D adaptive mesh refinement (AMR) code, which provides adequate resolution of the thin shell of matter behind the shock. We find that the sideways propagation is different than predicted by simple analytic models. The physical conditions at the sides of the jet are found to be significantly different than at the front of the jet, and most of the emission occurs within the initial opening angle of the jet. The light curves, as seen by observers at different viewing angles with respect to the jet axis, are then calculated assuming synchrotron emission. For an observer along the jet axis, we find a sharp achromatic `jet break' in the light curve at frequencies above the typical synchrotron frequency, at $t_{jet}\approx 5.8(E_{52}/n_1)^{1/3}(\theta_0/0.2)^{8/3}$ days, while the temporal decay index $\alpha$ ($F_{\nu}\propto t^{\alpha}$) after the break is steeper than $-p$ ($\alpha=-2.85$ for $p=2.5$). At larger viewing angles $t_{jet}$ increases and the jet break becomes smoother.

We review recent relativistic hydrodynamic simulations of jets, and their interpretation in terms of the results from linear stability analysis. These studies show that, interpreted naively, the distribution of synchrotron intensity will in general be a poor guide to the physical state (density and pressure) of the underlying flow, and that even if the physical state can be inferred, it, in turn, may prove to be a poor guide to the source dynamics, in terms of the transport of energy and momentum from the central engine. However, we demonstrate that an interplay of simulation and linear stability analysis provides a powerful tool for elucidating the nature and character of structures that jets may sustain. From such studies we can explain the complex behavior of observed jets, which manifest both stationary and propagating structures, without recourse to ad hoc macroscopic disturbances. This provides a framework for the interpretation of multi-epoch total intensity data wherein an understanding of the character of individual flow features will allow the effects of physical state and dynamics to be deconvolved.

We present proper motions obtained from a dual frequency, six-epoch, VLBA polarization experiment monitoring a sample of 12 blazars. The observations were made at 15 GHz and 22 GHz at bi-monthly intervals over 1996. Ten of the eleven sources for which proper motion could be reliably determined are superluminal. Only J2005+77 has no superluminal components. Three sources (OJ287, J1224+21, and J1512-09) show motion faster than 10h^{-1}c, requiring $\gamma_{pattern}$ of at least 10h^{-1}c (H_0 = 100h km/s/Mpc). We compare our results to those in the literature and find motions outside the previously observed range for four sources. While some jet components exhibit significant non-radial motion, most motion is radial. In at least two sources there are components moving radially at significantly different structural position angles. In five of six sources (3C120, J1224+21, 3C273, 3C279, J1512-09, and J1927+73) that have multiple components with measurable proper motion, the innermost component is significantly slower than the others, suggesting that acceleration occurs in the jet. In the motions of individual components we observe at least one decelerating motion and two ``bending'' accelerations which tend to align their motions with larger scale structure. We also discuss in detail our procedures for obtaining robust kinematical results from multi-frequency VLBI data spanning several epochs.

The amplitude of jet distortions and accompanying pressure and velocity fluctuations resulting from Kelvin-Helmholtz instability of three dimensional relativistic jets are explored. The effect of instability on jets as they accelerate from sub- to super-Alfv\'enic speeds is explored and a relativistic stabilization mechanism for trans-Alfv\'enic jets is proposed. The level to which asymmetric instabilities on supermagnetosonic relativistic jets will grow is predicted theoretically and a Doppler boosted ``apparent'' emissivity is computed. Effects due to helically twisted filamentary structure produced by asymmetric modes of instability should be readily observable on relativistic jets.

We compare results from a relativistic and a nonrelativistic set of 2D axisymmetric jet simulations. For a set of five relativistic simulations that either increase the Lorentz factor or decrease the adiabatic index we compute nonrelativistic simulations with equal useful power or thrust. We examine these simulations for morphological and dynamical differences, focusing on the velocity field, the width of the cocoon, the age of the jets, and the internal structure of the jet itself. The primary result of these comparisons is that the velocity field of nonrelativistic jet simulations cannot be scaled up to give the spatial distribution of Lorentz factors seen in relativistic simulations. Since the local Lorentz factor plays a major role in determining the total intensity for parsec scale extragalactic jets, this suggests that a nonrelativistic simulation cannot yield the proper intensity distribution for a relativistic jet. Another general result is that each relativistic jet and its nonrelativistic equivalents have similar ages (in dynamical time units, = R/a_a, where R is the initial radius of a cylindrical jet and a_a is the sound speed in the ambient medium). In addition to these comparisons, we have completed four new relativistic simulations to investigate the effect of varying thermal pressure on relativistic jets. The simulations generally confirm that faster (larger Lorentz factor) and colder jets are more stable, with smaller amplitude and longer wavelength internal variations. The apparent stability of these jets does not follow from linear normal mode analysis, which suggests that there are available growing Kelvin-Helmholtz modes. (Abridged.)

Results from a long-term program to quantify the range of behavior of the cm-wavelength total flux and linear polarization variability properties of a sample of 41 radio-loud BL Lac objects using weekly to tri-monthly observations with the University of Michigan 26-m telescope operating at 14.5, 8.0, and 4.8 GHz are presented; these observations are used to identify class-dependent differences between these BL Lacs and QSOs in the Pearson-Readhead sample. The BL Lacs are found to be more highly variable in total flux density than the QSOs, exhibiting changes that are often nearly-simultaneous and of comparable amplitude at 14.5 and 4.8 GHz in contrast to the behavior in the QSOs and supporting the existence of class-dependent differences in opacity within the parsec-scale jet flows. Structure function analyses of the flux observations quantify that a characteristic timescale is identifiable in only 1/3 of the BL Lacs. The time-averaged fractional linear polarizations are only on the order of a few percent and are consistent with the presence of tangled magnetic fields within the emitting regions. In many sources a preferred long-term orientation of the EVPA is present; when compared with the VLBI structural axis, no preferred position angle difference is identified. The polarized flux typically exhibits variability with timescales of months to a few years and shows the signature of a propagating shock during several resolved outbursts. The observations indicate that the source emission is predominately due to evolving source components and support the occurrence of more frequent shock formation in BL Lac parsec-scale flows than in QSO jets. The differences in variability behavior and polarization between BL Lacs and QSOs can be explained by differences in jet stability.

We use a continuous wavelet transform to analyze more than two decades of data for the BL~Lac object OJ~287 acquired as part of the UMRAO variability program. We find clear evidence for a persistent modulation of the total flux and polarization with period $\sim 1.66$ years, and for another signal that dominates activity in the 1980s with period $\sim 1.12$ years. The relationship between these two variations can be understood in terms of a `shock-in-jet' model, in which the longer time scale periodicity is associated with an otherwise quiescent jet, and the shorter time scale activity is associated with the passage of a shock; the different periodicities of these two components may reflect different internal conditions of the two flow domains, leading to different wave speeds, or different contractions of a single underlying periodicity, due to the different Doppler factors of the two flow components. We suggest that the modulation arises from a wave driven by some asymmetric disturbance close to the central engine. The periodic behavior in polarization exhibits excursions in $U$ which correspond to a direction $\sim 45^{\circ}$ from the VLBI jet axis. This behavior is not explained by the random walk in the $Q$-$U$ plane which is expected from models in which a pattern of randomly aligned magnetic field elements propagate across the visible portion of the flow, and suggests a small amplitude, cyclic variation in the flow direction in that part of the flow that dominates cm-wavelength emission.

A series of simulated maps showing the appearance in total intensity of flows computed using a recently developed relativistic hydrodynamic code (Duncan \& Hughes 1994: ApJ, 436, L119) are presented. The radiation transfer calculations were performed by assuming the flow is permeated by a magnetic field and fast particle distribution in energy equipartition, with energy density proportional to the hydrodynamic energy density (i.e., pressure). We find that relativistic flows subject to strong perturbations exhibit a density structure consisting of a series of nested bow shocks, and that this structure is evident in the intensity maps for large viewing angles. However, for viewing angles $<30^{\circ}$, differential Doppler boosting leads to a series of axial knots of emission, similar to the pattern exhibited by many VLBI sources. The appearance of VLBI knots is determined primarily by the Doppler boosting of parts of a more extended flow. To study the evolution of a perturbed jet, a time series of maps was produced and an integrated flux light curve created. The light curve shows features characteristic of a radio loud AGN: small amplitude variations and a large outburst. We find that in the absence of perturbations, jets with a modest Lorentz factor ($\sim 5$) exhibit complex intensity maps, while faster jets (Lorentz factor $\sim 10$) are largely featureless. We also study the appearance of kiloparsec jet-counterjet pairs by producing simulated maps at relatively large viewing angles; we conclude that observed hot spot emission is more likely to be associated with the Mach disk than with the outer, bow shock.

We describe a method for the numerical solution of the relativistic Euler equations which we have found to be both robust and efficient, and which has enabled us to simulate relativistic jets. The technique employs a solver of the Godunov-type, with approximate solution of the local Riemann problems, applied to laboratory frame variables. Lorentz transformations provide the rest frame quantities needed for the estimation of wave speeds, etc. This is applied within the framework of an adaptive mesh refinement algorithm, allowing us to perform high-resolution, 2-D simulations with modest computing resources. We present the results of nonrelativistic, and relativistic ($\gamma=5$ and $10$) runs, for adiabatic indices of $5/3$ and $4/3$. We find the same gross morphology in all cases, but the relativistic runs exhibit little instability and less well-defined structure internal to the jet: this might explain the difference between (relatively slow) BL~Lacs and (faster) QSOs. We find that the choice of adiabatic index makes a small but discernible difference to the structure of the shocked jet and ambient media.