We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

Context. Blazars are a subset of active galactic nuclei (AGN) with jets that are oriented along our line of sight. Variability and spectral energy distribution (SED) studies are crucial tools for understanding the physical processes responsible for observed AGN emission. Aims. We report peculiar behaviour in the bright gamma-ray blazar PKS 1424-418 and use its strong variability to reveal information about the particle acceleration and interactions in the jet. Methods. Correlation analysis of the extensive optical coverage by the ATOM telescope and nearly continuous gamma-ray coverage by the Fermi Large Area Telescope is combined with broadband, time-dependent modeling of the SED incorporating supplemental information from radio and X-ray observations of this blazar. Results. We analyse in detail four bright phases at optical-GeV energies. These flares of PKS 1424-418 show high correlation between these energy ranges, with the exception of one large optical flare that coincides with relatively low gamma-ray activity. Although the optical/gamma-ray behaviour of PKS 1424-418 shows variety, the multiwavelength modeling indicates that these differences can largely be explained by changes in the flux and energy spectrum of the electrons in the jet that are radiating. We find that for all flares the SED is adequately represented by a leptonic model that includes inverse Compton emission from external radiation fields with similar parameters. Conclusions. Detailed studies of individual blazars like PKS 1424-418 during periods of enhanced activity in different wavebands are helping us identify underlying patterns in the physical parameters in this class of AGN.

We have used photometry from the Kepler satellite to characterize the variability of four radio-loud active galactic nuclei (AGN) on timescales from years to minutes. The Kepler satellite produced nearly continuous high precision data sets which provided better temporal coverage than possible with ground based observations. We have now accumulated eleven quarters of data, eight of which were reported in our previous paper. In addition to constructing power spectral densities (PSDs) and characterizing the variability of the last three quarters, we have linked together the individual quarters using a multiplicative scaling process, providing data sets spanning ~2.8 years with >98% coverage at a 30 minute sampling rate. We compute PSDs on these connected data sets that yield power law slopes at low frequencies in the approximate range of -1.5 to -2.0, with white noise seen at higher frequencies. These PSDs are similar to those of both the individual quarters and to those of ground based optical observations of other AGN. We also have explored a PSD binning method intended to reduce a bias toward shallow slope fits by evenly distributing the points within the PSDs. This tends to steepen the computed PSD slopes, especially when the low frequencies are relatively poorly fit. We detected flares lasting several days in which the brightness increased by ~15-20% in one object, as well a smaller flare in another. Two AGN showed only small, ~1-2%, fluctuations in brightness.

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent \gamma-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each \gamma-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the \gamma-ray outbursts, although the \gamma-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the \gamma-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the \gamma-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that \gamma-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells.

We report on multiwavelength observations of the blazar PKS 0537-441 (z = 0.896) obtained from microwaves through gamma rays by SMA, REM, ATOM, Swift and Fermi during 2008 August-2010 April. Strong variability has been observed in gamma rays, with two major flaring episodes (2009 July and 2010 March) and a harder-when-brighter behaviour, quite common for FSRQs and low-synchrotron-peaked BL Lacs, in 2010 March. In the same way the SED of the source cannot be modelled by a simple synchrotron self-Compton model, as opposed to many BL Lacs, but the addition of an external Compton component of seed photons from a dust torus is needed. The 230 GHz light curve showed an increase simultaneous with the gamma-ray one, indicating co-spatiality of the mm and gamma-ray emission region likely at large distance from the central engine. The low, average, and high activity SED of the source could be fit changing only the electron distribution parameters, but two breaks in the electron distribution are necessary. The ensuing extra spectral break, located at NIR-optical frequencies, together with that in gamma rays seem to indicate a common origin, most likely due to an intrinsic feature in the underlying electron distribution. An overall correlation between the gamma-ray band with the R-band and K-band has been observed with no significant time lag. On the other hand, when inspecting the light curves on short time scales some differences are evident. In particular, flaring activity has been detected in NIR and optical bands with no evident gamma-ray counterparts in 2009 September and November. Moderate variability has been observed in X-rays with no correlation between flux and photon index. An increase of the detected X-ray flux with no counter part at the other wavelengths has been observed in 2008 October, suggesting once more a complex correlation between the emission at different energy bands.

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164 in 2008. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with ~0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary core and in the superluminal knot, both at >12 parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.

We have used Kepler photometry to characterize variability in four radio-loud active galactic nuclei (three quasars and one object tentatively identified as a Seyfert 1.5 galaxy) on timescales from minutes to months, comparable to the light crossing time of the accretion disk around the central supermassive black hole or the base of the relativistic jet. Kepler's almost continuous observations provide much better temporal coverage than is possible from ground-based observations. We report the first such data analyzed for quasars. We have constructed power spectral densities using 8 Kepler quarters of long-cadence (30-minute) data for three AGN, 6 quarters for one AGN and 2 quarters of short-cadence (1-minute) data for all four AGN. On timescales longer than about 0.2-0.6 day, we find red noise with mean power-law slopes ranging from -1.8 to -1.2, consistent with the variability originating in turbulence either behind a shock or within an accretion disk. Each AGN has a range of red noise slopes which vary slightly by month and quarter of observation. No quasi-periodic oscillations of astrophysical origin were detected. We detected several days-long flares when brightness increased by 3% - 7% in two objects. No flares on timescales of minutes to hours were detected. Our observations imply that the duty cycle for enhanced activity in these radio-loud AGN is small. These well-sampled AGN light curves provide an impetus to develop more detailed models of turbulence in jets and instabilities in accretion disks.

We present multiwavelength data of the blazar 3C 454.3 obtained during an extremely bright outburst from November 2010 through January 2011. These include flux density measurements with the Herschel Space Observatory at five submillimeter-wave and far-infrared bands, the Fermi Large Area Telescope at gamma-ray energies, Swift at X-ray, ultraviolet (UV), and optical frequencies, and the Submillimeter Array at 1.3 mm. From this dataset, we form a series of 52 spectral energy distributions (SEDs) spanning nearly two months that are unprecedented in time coverage and breadth of frequency. Discrete correlation anlaysis of the millimeter, far-infrared, and gamma-ray light curves show that the variations were essentially simultaneous, indicative of co-spatiality of the emission, at these wavebands. In contrast, differences in short-term fluctuations at various wavelengths imply the presence of inhomegeneities in physical conditions across the source. We locate the site of the outburst in the parsec-scale core, whose flux density as measured on 7 mm Very Long Baseline Array images increased by 70 percent during the first five weeks of the outburst. Based on these considerations and guided by the SEDs, we propose a model in which turbulent plasma crosses a conical standing shock in the parsec-scale region of the jet. Here, the high-energy emission in the model is produced by inverse Compton scattering of seed photons supplied by either nonthermal radiation from a Mach disk, thermal emission from hot dust, or (for X-rays) synchrotron radiation from plasma that crosses the standing shock. For the two dates on which we fitted the model SED to the data, the model corresponds very well to the observations at all bands except at X-ray energies, where the spectrum is flatter than observed.

The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and \gamma-ray flux apparently fades in time, likely because of an increasing optical to \gamma-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor \delta. Under the hypothesis that in the period 2008-2011 all the \gamma-ray and optical variability on a one-week timescale were due to changes in \delta, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle \theta\ only, then \theta\ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

We locate the gamma-ray and lower frequency emission in flares of the BL Lac object AO 0235+164 at >12pc in the jet of the source from the central engine. We employ time-dependent multi-spectral-range flux and linear polarization monitoring observations, as well as ultra-high resolution (~0.15 milliarcsecond) imaging of the jet structure at lambda=7mm. The time coincidence in the end of 2008 of the propagation of the brightest superluminal feature detected in AO 0235+164 (Qs) with an extreme multi-spectral-range (gamma-ray to radio) outburst, and an extremely high optical and 7mm (for Qs) polarization degree provides strong evidence supporting that all these events are related. This is confirmed at high significance by probability arguments and Monte-Carlo simulations. These simulations show the unambiguous correlation of the gamma-ray flaring state in the end of 2008 with those in the optical, millimeter, and radio regime, as well as the connection of a prominent X-ray flare in October 2008, and of a series of optical linear polarization peaks, with the set of events in the end of 2008. The observations are interpreted as the propagation of an extended moving perturbation through a re-collimation structure at the end of the jet's acceleration and collimation zone.

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with 0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.

We monitored the flaring blazar 3C 454.3 during 2005 June-July with the Spitzer Infrared Spectrograph (IRS: 15 epochs), Infrared Array Camera (IRAC: 12 epochs) and Multiband Imaging Photometer (MIPS: 2 epochs). We also made Spitzer IRS, IRAC, and MIPS observations from 2006 December-2007 January when the source was in a low state, the latter simultaneous with a single Chandra X-ray observation. In addition, we present optical and sub-mm monitoring data. The 2005-2007 period saw 3 major outbursts. We present evidence that the radio-optical SED actually consists of two variable synchrotron peaks, the primary at IR and the secondary at sub-mm wavelengths. The lag between the optical and sub-mm outbursts may indicate that these two peaks arise from two distinct regions along the jet separated by a distance of 0.07-5 pc. The flux at 5-35 microns varied by a factor of 40 and the IR peak varied in frequency from <1E13 Hz to 4E13 Hz between the highest and lowest states in 2005 and 2006, respectively. Variability was well correlated across the mid-IR band, with no measurable lag. Flares that doubled in flux occurred on a time scale of 3 days. The IR SED peak moved to higher frequency as a flare brightened, then returned to lower frequency as it decayed. The fractional variability amplitude increased with frequency, which we attribute to decreasing synchrotron-self absorption optical depth. Mid-IR flares may signal the re-energization of a shock that runs into inhomogeneities along the pre-existing jet or in the external medium. The synchrotron peak frequencies during each major outburst may depend upon both the distance from the jet apex and the physical conditions in the shocks. Variation of the Doppler parameter along a curved or helical jet is another possibility. Frequency variability of the IR synchrotron peak may have important consequences for the interpretation of the blazar sequence, and the presence of a secondary peak may give insight into jet structure.

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

The past two Decadal Surveys in Astronomy and Astrophysics recommended the completion of a space-based interferometry mission, known today as SIM PlanetQuest, for its unique ability to detect and characterize nearby rocky planets (Bahcall 1991, McKee & Taylor 2001), as well as contributions to a broad range of problems in astrophysics. Numerous committees of the National Research Council as well as NASA Roadmaps have similarly highlighted SIM as the one technology that offers detection and characterization of rocky planets around nearby stars and which is technically ready. To date, SIM remains the only program with the capability of detecting and confirming rocky planets in the habitable zones of nearby solar-type stars. Moreover, SIM measures masses and three-dimensional orbits of habitable planets around nearby stars (within 25 pc); these are the only stars for which follow-up by other techniques is feasible, such as space-based spectroscopy, ground-based interferometry, and of course TPF.

The X-ray emission from quasars such as 3C273 is generally agreed to arise from Compton scattering of low energy seed photons by relativistic electrons in a relativistic jet oriented close to the line of sight. However there are a number of possible models for the origin of the seed photons. In Paper I (McHardy et al 1999) we showed that the X-ray and IR variability from 3C273 was highly correlated in 1997, with the IR flux leading the X-rays by ~0.75 +/- 0.25 days. The strong correlation, and lag, supports the Synchrotron Self-Compton (SSC) model, where the seed photons are synchroton photons from the jet itself. The previous correlation was based on one moderately well sampled flare and another poorly sampled flare, so the possibility of chance correlated variability exists. Here we report on further X-ray and IR observations of 3C273 which confirm the behaviour seen in Paper I. During a 2 week period of observations we see a flare of amplitude ~25%, lasting for ~5 days, showing a high correlation between IR and X-ray variations, with the X-rays lagging by ~1.45+/- 0.15 days. These observations were not scheduled at any special time, implying that the same mechanism - almost certainly SSC - dominates the X-ray emission on most occasions and that the structure of the emission region is similar in most small flares.

We have calculated inverse-Compton Doppler factors for 3C 279 using the collection of VLBI data recently published by us, and the collection of multiwavelength spectra recently published by Hartman et al. From the Doppler factor and superluminal apparent speed, we then calculate the Lorentz factor and angle to the line-of-sight of the parsec-scale relativistic jet. We model the jet components as homogeneous spheres and the VLBI core as an unresolved inhomogeneous conical jet. The conical-jet model can be made to match both the observed X-ray emission and the VLBI properties of the core with a suitable choice of Doppler factor, implying the core makes a significant contribution to the X-ray emission. The parameters of the conical models indicate the jet is particle dominated at the radii that produce significant emission, and is not in equipartition. At the inner radius of the conical jet the magnetic field is of order 0.1 G and the relativistic-particle number density is of order 10 cm^{-3}. When all components are included in the calculation, then on average the core produces about half of the X-rays, with the other half being split between the long-lived component C4 and the brightest inner-jet component. We calculate an average speed and angle to the line-of-sight for the region of the jet interior to 1 mas of v=0.992c (gamma=8) and 4 degrees, and an average speed and angle to the line-of-sight for C4 (at a distance from the core of 3 mas) of v=0.997c (gamma=13) and 2 degrees. These values imply average Doppler factors of delta=12 for the inner jet, and delta=21 for C4.

We have used the VLBA to produce a high dynamic range image of the nucleus of NGC 6251 at 1.6 GHz and snapshot images at 5.0, 8.4, and 15.3 GHz to search for emission from a parsec-scale counterjet. Previous VLBI images at 1.6 GHz have set a lower limit for the jet/counterjet brightness ratio near the core at about 80:1, which is larger than expected given the evidence that the radio axis is fairly close to the plane of the sky. A possible explanation is that the inner few pc of the counterjet is hidden by free-free absorption by ionized gas associated with an accretion disk or torus. This would be consistent with the nearly edge-on appearance of the arcsecond-scale dust disk seen in the center of NGC 6251 by HST. We detect counterjet emission close to the core at 1.6 GHz, but not at the higher frequencies. Given that the optical depth of free-free absorption falls off more rapidly with increasing frequency than the optically thin synchrotron emission from a typical radio jet, this result implies that the absence of a detectable parsec-scale counterjet at high frequencies is not due to free-free absorption unless the density of ionized gas is extremely high and we have misidentified the core at 1.6 GHz. The most likely alternative is a large jet/counterjet brightness ratio caused by relativistic beaming, which in turn requires the inner radio axis to be closer to our line of sight than the orientation of the HST dust disk would suggest.

NGC4261 (3C270) is an elliptical galaxy containing a 300-parsec scale nuclear disk of gas and dust imaged by the Hubble Space Telescope (HST), around a central supermassive black hole. Previous VLBI observations of NGC4261 revealed a gap in emission in the radio counterjet, presumably due to free-free absorption in the inner parsec of the accretion disk. Here we present three 8 GHz VLBA observations of NGC4261 that allow us to monitor the location and depth of the gap, and check for motions in the jet and counterjet. The separation between the brightest peak and the gap is stable, with an upper limit to its motion of 0.01 c, supporting the interpretation of the gap as absorption by an accretion disk rather than an intrinsic jet feature. These observations span a time of order that required for orbiting material in the disk to transit the counterjet, so we are able to search for density changes (clumps) in the disk by monitoring the optical depth of the gap. The optical depth of the gap is stable to within 20% over 5 years at tau = 1.1 +/- 0.1, corresponding to an electron density in the disk that is constant to within 10%. We measure an apparent speed in the jet of 0.52 +/- 0.07 c. An apparent speed could not be measured for the counterjet due to a lack of identifiable features. From the apparent jet speed and the jet to counterjet brightness ratio, we calculate the viewing angle of the jet to be 63 +/- 3 degrees and its intrinsic speed to be 0.46 +/- 0.02 c. From the inclination and position angles of the parsec-scale radio jet and outer HST disk rotation axis we calculate a difference between the parsec-scale radio jet and outer HST disk rotation axis of 12 +/- 2 degrees.

We examine the coincidence of times of high $\gamma$-ray flux and ejections of superluminal components from the core in EGRET blazars based on a VLBA monitoring program at 22 and 43 GHz from November 1993 to July 1997. In 23 cases of $\gamma$-ray flares for which sufficient VLBA data exist, 10 of the flares (in 8 objects) fall within 1$\sigma$ uncertainties of the extrapolated epoch of zero separation from the core of a superluminal radio component. In each of two sources (0528+134 and 1730-130) two successive $\gamma$-ray flares were followed by the appearance of new superluminal components. We carried out statistical simulations which show that if the number of coincidences $\ge$ 10 the radio and $\gamma$-ray events are associated with each other at >99.999% confidence. Our analysis of the observed behavior, including variability of the polarized radio flux, of the sources before, during, and after the $\gamma$-ray flares suggests that the $\gamma$-ray events occur in the superluminal radio knots. This implies that the $\gamma$-ray flares are caused by inverse Compton scattering by relativistic electrons in the parsec-scale regions of the jet rather than closer to the central engine.

The physical conditions in the inner parsec of accretion disks believed to orbit the central black holes in active galactic nuclei can be probed by imaging the absorption (by ionized gas in the disk) of background emission from a radio counterjet. We report high angular resolution VLBI observations of the nearby (about 40 Mpc) radio galaxy NGC 4261 that confirm free-free absorption of radio emission from a counterjet by a geometrically thin, nearly edge-on disk at 1.6, 4.8, and 8.4 GHz. The angular width and depth of the absorption appears to increase with decreasing frequency, as expected. We derive an average electron density of ~10E4 per cc at a disk radius of about 0.2 pc, assuming that the inner disk inclination and opening angles are the same as at larger radii. Pressure balance between the thermal gas and the magnetic field in the disk implies an average field strength of 0.1 milligauss at a radius of 0.2 pc. These are the closest-in free-free absorption measurements to date of the conditions in an extragalactic accretion disk orbiting a black hole with a well-determined mass. If a standard advection-dominated accretion flow exists in the disk center, then the transition between thin and thick disk regions must occur at a radius less than 0.2 pc (4000 Schwarzschild radii).

We present the results of a program to monitor the structure of the radio emission in 42 $\gamma$-ray bright blazars (31 quasars and 11 BL Lac objects) with the VLBA at 43, 22, and occasionally 15 and 8.4 GHz, over the period from November 1993 to July 1997. We determine proper motions in 33 sources and find that the apparent superluminal motions in $\gamma$-ray sources are much faster than for the general population of bright compact radio sources. This follows the strong dependence of the $\gamma$-ray flux on the level of relativistic beaming for both external-radiation Compton and synchrotron self-Compton emission. There is a positive correlation (correlation coefficient $r$=0.45) between the flux density of the VLBI core and the $\gamma$-ray flux and a moderate correlation (partial correlation coefficient $r$=0.31) between $\gamma$-ray apparent luminosity and superluminal velocities of jet components, as expected if the $\gamma$-ray emission originates in a very compact region of the relativistic jet and is highly beamed. In 43% of the sources the jet bends by more than 20$^\circ$ on parsec scales, which is consistent with amplification by projection effects of modest actual changes in position angle. In 27 of the sources in the sample there is at least one non-core component that appears to be stationary during our observations. Different characteristics of stationary features close to and farther from the core lead us to suggest two different classes of stationary components: those within about 2 milliarcseconds (mas) of the core, probably associated with standing hydrodynamical compressions, and those farther down the jet, which tend to be associated with bends in the jet.

We present results of long-term high-frequency VLBI monitoring of the relativistic jet in 3C279, consisting of 18 epochs at 22 GHz from 1991 to 1997 and 10 epochs at 43 GHz from 1995 to 1997. Three major results of this study are: apparent speeds measured for six superluminal components range from 4.8 to 7.5 c (H_{0}=70 km s^{-1} Mpc^{-1}, q_{0}=0.1), variations in the total radio flux are due primarily to changes in the VLBI core flux, and the uniform-sphere brightness temperature of the VLBI core is about 1 x 10^{13} K at 22 GHz after 1995, one of the highest direct estimates of a brightness temperature. If the variability brightness temperature measured for 3C279 by Lahteenmaki & Valtaoja is an actual value and not a lower limit, then the rest-frame brightness temperature of 3C279 is quite high and limited by inverse Compton effects rather than equipartition. The parsec-scale morphology of 3C279 consists of a bright, compact VLBI core, a jet component (C4) that moved from about 2 mas to about 3.5 mas from the core during the course of our monitoring, and an inner jet that extends from the core to a stationary component, C5, at about 1 mas from the core. Components in the inner jet are relatively short-lived, and fade by the time they reach about 1 mas from the core. The components have different speeds and position angles from each other, but these differences do not match the differences predicted by the precession model of Abraham & Carrara. Although VLBI components were born about six months prior to each of the two observed gamma-ray high states, the sparseness of the gamma-ray data prevents a statistical analysis of possible correlations.

The physical conditions in the inner parsec of accretion disks believed to orbit the central black holes in active galactic nuclei can be probed by imaging the absorption of background radio emission by ionized gas in the disk. We report high angular resolution observations of the nearby galaxy NGC 4261 which show evidence for free-free absorption by a thin, nearly edge-on disk at several frequencies. The angular width, and probably the depth, of the absorption appears to increase with decreasing frequency, as expected. Because free-free absorption is much larger at lower frequencies, the longest possible baselines are needed to provide adequate angular resolution; observing at higher frequencies to improve resolution will not help.

We present the first VLBI Space Observatory Programme (VSOP) observations of the gamma-ray blazar 3C 279 at 1.6 and 5 GHz. The combination of the VSOP and VLBA-only images at these two frequencies maps the jet structure on scales from 1 to 100 mas. On small angular scales the structure is dominated by the quasar core and the bright secondary component `C4' located 3 milliarcseconds from the core (at this epoch). On larger angular scales the structure is dominated by a jet extending to the southwest, which at the largest scale seen in these images connects with the smallest scale structure seen in VLA images. We have exploited two of the main strengths of VSOP: the ability to obtain matched-resolution images to ground-based images at higher frequencies and the ability to measure high brightness temperatures. A spectral index map was made by combining the VSOP 1.6 GHz image with a matched-resolution VLBA-only image at 5 GHz from our VSOP observation on the following day. The spectral index map shows the core to have a highly inverted spectrum, with some areas having a spectral index approaching the limiting value for synchrotron self-absorbed radiation of 2.5. Gaussian model fits to the VSOP visibilities revealed high brightness temperatures (>10^{12} K) that are difficult to measure with ground-only arrays. An extensive error analysis was performed on the brightness temperature measurements. Most components did not have measurable brightness temperature upper limits, but lower limits were measured as high as 5x10^{12} K. This lower limit is significantly above both the nominal inverse Compton and equipartition brightness temperature limits. The derived Doppler factor, Lorentz factor, and angle to the line-of-sight in the case of the equipartition limit are at the upper end of the range of expected values for EGRET blazars.

The structure of AGN accretion disks on sub-parsec scales can be probed through free-free absorption of synchrotron emission from the base of symmetric radio jets. We present new VLBA observations of the nearby FR-I radio galaxy NGC 4261 at 22 and 43 GHz, and combine these with previous VLBA observations at 1.6 and 8.4 GHz to map absorption caused by an inner accretion disk. Assuming the disk is geometrically and optically thin and composed of a uniform 10^4 K plasma, the average electron density in the inner 0.1 pc is 10^3 - 10^8 cm^-3. Equating thermal gas pressure and magnetic field strength gives a disk magnetic field of 10^-4 - 10^-2 Gauss at 0.1 pc. The jet opening angle is between 0.3 and 20 degrees during the first 0.2 pc of the jet, and must be less than 5 degrees during the first 0.8 pc. We include an appendix containing expressions for a simple, optically thin, gas pressure dominated accretion disk model which may be applicable to other galaxies in addition to NGC 4261.

The gamma-ray blazar 3C 279 was monitored on a nearly daily basis with IUE, ROSAT and EGRET for three weeks between December 1992 and January 1993. During this period, the blazar was at a historical minimum at all wavelengths. Here we present the UV data obtained during the above multiwavelength campaign. A maximum UV variation of ~50% is detected, while during the same period the X-ray flux varied by no more than 13%. At the lowest UV flux level the average spectrum in the 1230-2700 A interval is unusually flat for this object (<alpha_UV>~1). The flattening could represent the lowest energy tail of the inverse Compton component responsible for the X-ray emission, or could be due to the presence of a thermal component at ~20000 K possibly associated with an accretion disk. The presence of an accretion disk in this blazar object, likely observable only in very low states and otherwise hidden by the beamed, variable synchrotron component, would be consistent with the scenario in which the seed photons for the inverse Compton mechanism producing the gamma-rays are external to the relativistic jet. We further discuss the long term correlation of the UV flux with the X-ray and gamma-ray fluxes obtained at various epochs. All UV archival data are included in the analysis. Both the X- and gamma-ray fluxes are generally well correlated with the UV flux, approximately with square root and quadratic dependences, respectively.

We present thirty VLBI images of the TeV blazar Markarian 421 (1101+384) at fifteen epochs spanning the time range from 1994 to 1997, and at six different frequencies from 2.3 to 43 GHz. The imaged observations include a high-resolution 5 GHz VLBI Space Observatory Programme (VSOP) observation with the HALCA satellite on 1997 November 14; full-track VLBA observations from 1994 April, 1996 November, and 1997 May at frequencies between 5 and 43 GHz; six epochs of VLBA snapshot observations at frequencies between 2 and 15 GHz from Radio Reference Frame studies; and five geodetic VLBI observations at 2 and 8 GHz from the archive of the Washington VLBI Correlator Facility located at the U.S. Naval Observatory. The dense time coverage of the images allows us to unambiguously track components in the parsec-scale jet over the observed time range. We measure the speeds of three inner jet components located between 0.5 and 5 mas from the core (0.3 to 3 pc projected linear distance) to be 0.19 +/- 0.27, 0.30 +/- 0.07, and -0.07 +/- 0.07 c (H_{0}=65 km s^{-1} Mpc^{-1}). If the sole 43 GHz image is excluded, all measured speeds are consistent with no motion. These speeds differ from tentative superluminal speeds measured by Zhang & B\aa\aa th from three epochs of data from the early 1980's. Possible interpretations of these subluminal speeds in terms of the high Doppler factor demanded by the TeV variability of this source are discussed.

The blazar 3C 279, one of the brightest identified extragalactic objects in the gamma-ray sky, underwent a flare of a factor 10 amplitude in gamma-rays towards the end of a 3-week pointing by CGRO, in 1996 January-February. The flare peak represents the highest gamma-ray intensity ever recorded for this object. During the high state, extremely rapid gamma-ray variability was seen. Coordinated multifrequency observations were carried out with RXTE, ASCA, ROSAT and IUE and from many ground-based observatories, covering most accessible wavelengths. The well-sampled, simultaneous RXTE light curve shows an outburst of lower amplitude (factor of ~3) well correlated with the gamma-ray flare without any apparent lag. The optical-UV light curves, which are not well sampled during the high energy flare, exhibit more modest variations (factor of ~2) and a lower degree of correlation. The flux at millimetric wavelengths was near an historical maximum during the gamma-ray flare peak, with suggestion of a correlated decay. We present simultaneous spectral energy distributions of 3C 279 prior to and near to the flare peak. The gamma-rays vary by more than the square of the observed IR-optical flux change, which poses some problems for specific blazar emission models. The synchrotron-self Compton model would require that the largest synchrotron variability occurred in the mostly unobserved sub-mm/far-infrared region. Alternatively, a large variation in the external photon field could occur over a time scale of few days. This occurs naturally in the ``mirror'' model, wherein the flaring region in the jet photoionizes nearby broad-emission-line clouds, which in turn provide soft photons that are Comptonized to gamma-ray energies.

We observed the nuclear region of NGC 4261 (3C270) with the VLBA at two frequencies (1.6 and 8.4 GHz) simultaneously. We find that the position angle of the parsec-scale radio axis agrees, within the errors, with the position angle of the VLA-scale jet. Our observations also reveal basically symmetric radio structures at both 1.6 and 8.4 GHz. Analysis of these images shows that most of the central 10 pc of this source is not significantly affected by free-free absorption, even though HST images show that the nucleus contains a nearly edge-on disk of gas and dust on larger scales. Our highest angular resolution image shows a very narrow gap in emission, which we interpret as an absorption feature, just east of the radio core. This suggests that there may be a small, dense inner accretion disk whose width is less than 0.1 pc. If the inclination of this inner disk is close to that of the larger-scale HST disk it becomes optically thin to 8.4 GHz radiation at a deprojected radius of about 0.8 pc. The brightness of the pc-scale jets falls off very rapidly on both sides of the core, suggesting that the jets are rapidly expanding during the the first several pc of their travel. It appears that there is a small dense inner disk centered on the radio core (the base of the jets; less than 1 pc), a low density bubble filling most of the the inner several pc of the nucleus (within which the radio jets expand rapidly; ~10 pc), and a surrounding cool, higher density region (of which the HST absorption disk is part; tens to hundreds of pc) within which the transverse expansion of the radio jets, as implied by the rate of decrease in jet brightness, is nearly halted.

We present 5 GHz global VLBI observations of the two Fanaroff Riley Type I radio galaxies B2 0836+29 and 3C465 (2335+26). For 3C465 we present also 1.7 GHz and 8.4 GHz global VLBI data. In addition VLA observations were used to obtain arsecond resolution continuum and polarization maps at 5 GHz. Both sources are very asymmetric on the parsec-scale, with a core and a one-sided jet, aligned with the main arcsecond scale jet. We place a limit on the milliarcsecond jet to counterjet brightness ratio B$_{jet}$/B$_{cjet}$ \gtsim 20 and \gtsim 30 for B2 0836+29 and 3C465 respectively. For 3C465 the strong asymmetry holds to the kiloparsec scale. The brightness asymmetry and the ratio between the core radio power and total radio power allow us to constrain the jet velocity close to the core and the orientation of the radio structure with respect to the line of sight. The results suggest that the plasma speed is relativistic on the parsec scale for both sources, i.e. v$_{jet}$ \gtsim 0.75c for B2 0836+29 and v$_{jet}$ \gtsim 0.6c for 3C465. While v$_{jet}$ decreases from the parsec to the kiloparsec scale in B2 0836+29, in 3C465 the very high v$_{jet}$ holds all the way to the kiloparsec-scale {\t bright spot}. Our results are in agreement with the unification scheme suggestion that FR-I radio galaxies are the unbeamed poulation of BL-Lac objects. Furthermore, they reinforce the idea that the central engine in FR-I and FR-II radio galaxies must be qualitatively similar. The different radio morphology could then be due either to an intrinsically different nuclear power, which affects the torus geometry or to different conditions in the region beyond the parsec scale, where a significant deceleration in the FR-I jets occurs.

We present observations at 1.7 and 8.4 GHz of two Compact Steep Spectrum (CSS) sources from a complete sample of low-intermediate power radio galaxies. 3C346 shows an asymmetric structure with a one-sided ``jet'' and ``hot spot''. Present observations suggest that the classification of this source as a CSS is inappropriate, and that it is a common radio galaxy at a small angle to the line of sight. Its properties are in agreement with the predictions of unified schemes models. 4C31.04 shows more complex structure with the possibility of a centrally located flat spectrum core in between two close lobes. We suggest that this source could be a low redshift Compact Symmetric Object.

We present here new VLBI observations of one FR-I radio galaxy (NGC2484) and two Broad Line FR-II radio galaxies (3C109 and 3C382). For 3C109 new VLA maps are also shown. These sources belong to a complete sample of radio galaxies under study for a better knowledge of their structures at parsec resolution. The parsec structure of these 3 objects is very similar: asymmetric emission which we interpret as the core plus a one-sided jet. The parsec scale jet is always on the same side of the main kpc-scale jet. The limit on the jet to counterjet brightness ratio, the ratio of the core radio power to the total radio power and the Synchrotron-self Compton model allow us to derive some constraints on the jet velocity and orientation with respect to the line of sight. From these data and from those published on 2 other sources of our sample, we suggest that parsec scale jets are relativistic in both FR-I and FR-II radio galaxies and that pc scale properties in FR-I and FR-II radio galaxies are very similar despite of the large difference between these two classes of radio galaxies on the kpc scale.