We present an analysis of the core linear polarization properties of 387 parsec-scale active galactic nuclei (AGN) jets. Using 15 GHz VLBA data, we revisit the conclusions of the first paper in this series with multi-epoch measurements and more detailed analysis of a larger AGN sample which spans a broader range of synchrotron peak frequencies. Each AGN has been observed for at least five epochs between 1996 and 2017. We find that BL Lac objects have core electric vector position angles (EVPAs) which tend towards alignment with the local jet direction; compared to flat spectrum radio quasars (FSRQs), their EVPAs are also less variable over time. The AGN cores which are most fractionally polarized and least variable in polarization have EVPAs that are closely aligned with the local jet direction; they also have low variability in EVPA. These results support the popular model of a standing transverse shock at the base of the jet which collimates the jet magnetic field perpendicular to the jet direction, increasing the fractional polarization and leading to greater polarization stability over time. High-synchrotron-peaked (HSP) BL Lac objects form a low luminosity, low fractional polarization population. The five narrow-line Seyfert 1 galaxies in our sample have low fractional polarization and large EVPA-jet misalignments. Although AGN detected at gamma-rays are thought to be more Doppler boosted than non-detected AGN, we find no significant differences in fractional polarization based on detection by Fermi-LAT; the gamma-loud AGN are, however, more variable in core EVPAs.

We used 15 GHz multi-epoch Very Long Baseline Array (VLBA) polarization sensitive observations of 484 sources within a time interval 1996--2016 from the MOJAVE program, and also from the NRAO data archive. We have analyzed the linear polarization characteristics of the compact core features and regions downstream, and their changes along and across the parsec-scale active galactic nuclei (AGN) jets. We detected a significant increase of fractional polarization with distance from the radio core along the jet as well as towards the jet edges. Compared to quasars, BL Lacs have a higher degree of polarization and exhibit more stable electric vector position angles (EVPAs) in their core features and a better alignment of the EVPAs with the local jet direction. The latter is accompanied by a higher degree of linear polarization, suggesting that compact bright jet features might be strong transverse shocks, which enhance magnetic field regularity by compression.

We present 5321 milliarcsecond-resolution total intensity and linear polarization maps of 437 active galactic nuclei (AGNs) obtained with the VLBA at 15 GHz as part of the MOJAVE survey, and also from the NRAO data archive. The former is a long-term program to study the structure and evolution of powerful parsec-scale outflows associated with AGNs. The targeted AGNs are drawn from several flux-limited radio and gamma-ray samples, and all have correlated VLBA flux densities greater than about 50 mJy at 15 GHz. Approximately 80% of these AGNs are associated with gamma-ray sources detected by the Fermi LAT instrument. The vast majority were observed with the VLBA on 5 to 15 occasions between 1996 January 19 and 2016 December 26, at intervals ranging from a month to several years, with the most typical sampling interval being six months. A detailed analysis of the linear and circular polarization evolution of these AGN jets are presented in other papers in this series.

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.