By performing Monte Carlo simulations of the evolution of binary primordial black hole (PBH) systems, we estimate the effect of distant encounters with single PBHs upon the coalescence time and merger rate of binary PBHs. We find that, for models where PBHs compose a large fraction of dark matter, $\fPBH\sim 1$, the expected fractional change in coalescence time is negligible, of order $10^{-6}$ for most binaries. For models with significantly lower PBH abundances, $\fPBH\ll 1$, we find that the average change in binary lifetime due to encounters can be as large as $\mathcal{O}(10^{-2})$, with a small number of binaries experiencing an order unity change in lifetime. In the absence of encounters, we also compare the use of an analytic approximation for the coalescence time to numerically evolving the binary system, finding that the analytic approximation results in an order $10\%$ error in the coalescence time. However, when these effects are taken into consideration, there is a negligible change to the calculated merger rate, placing previous constraints on the PBH abundance arising from observed gravitational wave signals from merging binary black holes on a more secure footing.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

We consider the application of peaks theory to the calculation of the number density of peaks relevant for primordial black hole (PBH) formation. For PBHs, the final mass is related to the amplitude and scale of the perturbation from which it forms, where the scale is defined as the scale at which the compaction function peaks. We therefore extend peaks theory to calculate not only the abundance of peaks of a given amplitude, but peaks of a given amplitude and scale. A simple fitting formula is given in the high-peak limit relevant for PBH formation. We also adapt the calculation to use a Gaussian smoothing function, ensuring convergence regardless of the choice of power spectrum.

We present a novel analysis technique for liquid xenon time projection chambers that allows for a lower threshold by relying on events with a prompt scintillation signal consisting of single detected photons. The energy threshold of the LUX dark matter experiment is primarily determined by the smallest scintillation response detectable, which previously required a 2-fold coincidence signal in its photomultiplier arrays, enforced in data analysis. The technique presented here exploits the double photoelectron emission effect observed in some photomultiplier models at vacuum ultraviolet wavelengths. We demonstrate this analysis using an electron recoil calibration dataset and place new constraints on the spin-independent scattering cross section of weakly interacting massive particles (WIMPs) down to 2.5 GeV/c$^2$ WIMP mass using the 2013 LUX dataset. This new technique is promising to enhance light WIMP and astrophysical neutrino searches in next-generation liquid xenon experiments.

In this paper, the criterion used to determine whether a density perturbation will collapse to form a primordial black hole (PBH) is re-examined, in respect of its use to determine the abundance of PBHs. There is particular focus on which parameter to use, the time at which the abundance should be calculated, and the use of different smoothing functions. It is concluded that, with the tools currently available, the smoothed density contrast should be used rather than the peak value, and should be calculated from the time-independent component of the density contrast in the super-horizon limit (long before perturbations enter the horizon) rather than at horizon crossing. For the first time the effect of the choice of smoothing function upon the formation criterion is calculated, and, for a given abundance of PBHs, it is found that the uncertainty in the amplitude of the power spectrum due to this is $\mathcal{O}(10\%)$, an order of magnitude smaller than previous calculations suggest. The relation between the formation criterion stated in terms of the density contrast and the curvature perturbation $\mathcal{R}$ is also discussed.

The formation and abundance of primordial black holes (PBHs) arising from the curvature perturbation $\zeta$ is studied. The non-linear relation between $\zeta$ and the density contrast $\delta$ means that, even when $\zeta$ has an exactly Gaussian distribution, significant non-Gaussianities affecting PBH formation must be considered. Numerical simulations are used to investigate the critical value and the mass of PBHs which form, and peaks theory is used to calculate the mass fraction of the universe collapsing to form PBHs at the time of formation. A formalism to calculate the total present day PBH abundance and mass function is also derived. It is found that the abundance of PBHs is very sensitive to the non-linear effects, and that the power spectrum $\mathcal{P}_\zeta$ must be a factor of $\sim2$ larger to produce the same number of PBHs as the linear model (where the exact value depends on the critical value for a region to collapse and form a PBH). This also means that the derived constraints on the small-scale power spectrum from constraints on the abundance of PBHs are weaker by the same factor.

If the primordial curvature perturbation followed a Gaussian distribution, primordial black holes (PBHs) will be Poisson distributed with no additional clustering. We consider local non-Gaussianity and its impact on the initial PBH clustering and mass function due to mode coupling between long and short wavelength modes. We show that even a small amount of non-Gaussianity results in a significant enhancement on the PBH initial clustering and subsequent merger rate and that the PBH mass function shifts to higher mass PBHs. However, as the clustering becomes strong, the local number density of PBHs becomes large, leading to a large theoretical uncertainty in the merger rate.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately $0.7M_\odot$, with the formation rate increasing by at least two orders of magnitude due to the softening of the equation of state at this time, with a range of approximately $0.3M_\odot<M<1.4M_\odot$ at full width half-maximum. PBH formation is increased by a smaller amount for PBHs with masses spanning a large range, $10^{-3}M_\odot<M_{\rm PBH}<10^{3}M_\odot$, which includes the masses of the BHs that LIGO detected. The most significant source of uncertainty in the number of PBHs formed is now due to unknowns in the formation process, rather than from the phase transition. A near scale-invariant density power spectrum tuned to generate a population with mass and merger rate consistent with that detected by LIGO should also produce a much larger energy density of PBHs with solar mass. The existence of BHs below the Chandresekhar mass limit would be a smoking gun for a primordial origin and they could arguably constitute a significant fraction of the cold dark matter density. They also pose a challenge to inflationary model building which seek to produce the LIGO BHs without overproducing lighter PBHs.

In this thesis, the effect of non-Gaussianity upon the abundance of primordial black holes (PBHs), and the implications of such an effect are considered. It is shown that even small non-Gaussianity parameters can have a large effect on the constraints that can be placed on the primordial curvature perturbation power spectrum - which can become stronger or weaker by an order of magnitude. The effects of super-horizon curvature perturbation modes at the time of PBH formation are considered, and it is shown that these have little effect on the formation of a PBH, but can have an indirect effect on the abundance of PBHs due to modal coupling to horizon-scale modes in the presence of non-Gaussianity. By taking into account the effect of modal coupling to CMB-scale modes, many models can be ruled out as a mechanism to produce enough PBHs to constitute dark matter.

We report an absolute calibration of the ionization yields($\textit{Q$_y$})$ and fluctuations for electronic recoil events in liquid xenon at discrete energies between 186 eV and 33.2 keV. The average electric field applied across the liquid xenon target is 180 V/cm. The data are obtained using low energy $^{127}$Xe electron capture decay events from the 95.0-day first run from LUX (WS2013) in search of Weakly Interacting Massive Particles (WIMPs). The sequence of gamma-ray and X-ray cascades associated with $^{127}$I de-excitations produces clearly identified 2-vertex events in the LUX detector. We observe the K- (binding energy, 33.2 keV), L- (5.2 keV), M- (1.1 keV), and N- (186 eV) shell cascade events and verify that the relative ratio of observed events for each shell agrees with calculations. The N-shell cascade analysis includes single extracted electron (SE) events and represents the lowest-energy electronic recoil $\textit{in situ}$ measurements that have been explored in liquid xenon.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

The first searches for axions and axion-like particles with the Large Underground Xenon (LUX) experiment are presented. Under the assumption of an axio-electric interaction in xenon, the coupling constant between axions and electrons, gAe is tested, using data collected in 2013 with an exposure totalling 95 live-days $\times$ 118 kg. A double-sided, profile likelihood ratio statistic test excludes gAe larger than 3.5 $\times$ 10$^{-12}$ (90% C.L.) for solar axions. Assuming the DFSZ theoretical description, the upper limit in coupling corresponds to an upper limit on axion mass of 0.12 eV/c$^{2}$, while for the KSVZ description masses above 36.6 eV/c$^{2}$ are excluded. For galactic axion-like particles, values of gAe larger than 4.2 $\times$ 10$^{-13}$ are excluded for particle masses in the range 1-16 keV/c$^{2}$. These are the most stringent constraints to date for these interactions.

We present experimental constraints on the spin-dependent WIMP-nucleon elastic cross sections from the total 129.5 kg-year exposure acquired by the Large Underground Xenon experiment (LUX), operating at the Sanford Underground Research Facility in Lead, South Dakota (USA). A profile likelihood ratio analysis allows 90% CL upper limits to be set on the WIMP-neutron (WIMP-proton) cross section of $\sigma_n$ = 1.6$\times 10^{-41}$ cm$^{2}$ ($\sigma_p$ = 5$\times 10^{-40}$ cm$^{2}$) at 35 GeV$c^{-2}$, almost a sixfold improvement over the previous LUX spin-dependent results. The spin-dependent WIMP-neutron limit is the most sensitive constraint to date.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

We study the significance of mergers in the quenching of star formation in galaxies at z~1 by examining their color-mass distributions for different morphology types. We perform two-dimensional light profile fits to GOODS iz images of ~5000 galaxies and X-ray selected active galactic nucleus (AGN) hosts in the CANDELS/GOODS-north and south fields in the redshift range 0.7<z<1.3. Distinguishing between bulge-dominated and disk-dominated morphologies, we find that disks and spheroids have distinct color-mass distributions, in agreement with studies at z~0. The smooth distribution across colors for the disk galaxies corresponds to a slow exhaustion of gas, with no fast quenching event. Meanwhile, blue spheroids most likely come from major mergers of star-forming disk galaxies, and the dearth of spheroids at intermediate green colors is suggestive of rapid quenching. The distribution of moderate luminosity X-ray AGN hosts is even across colors, in contrast, and we find similar numbers and distributions among the two morphology types with no apparent dependence on Eddington ratio. The high fraction of bulge-dominated galaxies that host an AGN in the blue cloud and green valley is consistent with the scenario in which the AGN is triggered after a major merger, and the host galaxy then quickly evolves into the green valley. This suggests AGN feedback may play a role in the quenching of star formation in the minority of galaxies that undergo major mergers.

We report constraints on spin-independent weakly interacting massive particle (WIMP)-nucleon scattering using a 3.35e4 kg-day exposure of the Large Underground Xenon (LUX) experiment. A dual-phase xenon time projection chamber with 250 kg of active mass is operated at the Sanford Underground Research Facility under Lead, South Dakota (USA). With roughly fourfold improvement in sensitivity for high WIMP masses relative to our previous results, this search yields no evidence of WIMP nuclear recoils. At a WIMP mass of 50 GeV/c^2, WIMP-nucleon spin-independent cross sections above 2.2e-46 cm^2 are excluded at the 90% confidence level. When combined with the previously reported LUX exposure, this exclusion strengthens to 1.1e-46 cm^2 at 50 GeV/c^2.

We conduct a multiwavelength continuum variability study of the Seyfert 1 galaxy NGC 5548 to investigate the temperature structure of its accretion disk. The 19 overlapping continuum light curves (1158 to 9157 angstroms) combine simultaneous HST , Swift , and ground-based observations over a 180 day period from 2014 January to July. Light-curve variability is interpreted as the reverberation response of the accretion disk to irradiation by a central time-varying point source. Our model yields the disk inclination, i, temperature T1 at 1 light day from the black hole, and a temperature-radius slope, alpha. We also infer the driving light curve and find that it correlates poorly with both the hard and soft X-ray light curves, suggesting that the X-rays alone may not drive the ultraviolet and optical variability over the observing period. We also decompose the light curves into bright, faint, and mean accretion-disk spectra. These spectra lie below that expected for a standard blackbody accretion disk accreting at L/LEdd = 0.1

The Large Underground Xenon (LUX) experiment is a dual-phase liquid xenon time projection chamber (TPC) operating at the Sanford Underground Research Facility in Lead, South Dakota. A calibration of nuclear recoils in liquid xenon was performed $\textit{in situ}$ in the LUX detector using a collimated beam of mono-energetic 2.45 MeV neutrons produced by a deuterium-deuterium (D-D) fusion source. The nuclear recoil energy from the first neutron scatter in the TPC was reconstructed using the measured scattering angle defined by double-scatter neutron events within the active xenon volume. We measured the absolute charge ($Q_{y}$) and light ($L_{y}$) yields at an average electric field of 180 V/cm for nuclear recoil energies spanning 0.7 to 74 keV and 1.1 to 74 keV, respectively. This calibration of the nuclear recoil signal yields will permit the further refinement of liquid xenon nuclear recoil signal models and, importantly for dark matter searches, clearly demonstrates measured ionization and scintillation signals in this medium at recoil energies down to $\mathcal{O}$(1 keV).

We present constraints on weakly interacting massive particles (WIMP)-nucleus scattering from the 2013 data of the Large Underground Xenon dark matter experiment, including $1.4\times10^{4}\;\mathrm{kg\; day}$ of search exposure. This new analysis incorporates several advances: single-photon calibration at the scintillation wavelength, improved event-reconstruction algorithms, a revised background model including events originating on the detector walls in an enlarged fiducial volume, and new calibrations from decays of an injected tritium $\beta$ source and from kinematically constrained nuclear recoils down to 1.1 keV. Sensitivity, especially to low-mass WIMPs, is enhanced compared to our previous results which modeled the signal only above a 3 keV minimum energy. Under standard dark matter halo assumptions and in the mass range above 4 $\mathrm{GeV}\,c^{-2}$, these new results give the most stringent direct limits on the spin-independent WIMP-nucleon cross section. The 90% C.L. upper limit has a minimum of 0.6 zb at 33 $\mathrm{GeV}\,c^{-2}$ WIMP mass.

We present measurements of the electron-recoil (ER) response of the LUX dark matter detector based upon 170,000 highly pure and spatially-uniform tritium decays. We reconstruct the tritium energy spectrum using the combined energy model and find good agreement with expectations. We report the average charge and light yields of ER events in liquid xenon at 180 V/cm and 105 V/cm and compare the results to the NEST model. We also measure the mean charge recombination fraction and its fluctuations, and we investigate the location and width of the LUX ER band. These results provide input to a re-analysis of the LUX Run3 WIMP search.

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

Primordial black holes represent a unique probe to constrain the early universe on small scales - providing the only constraints on the primordial power spectrum on the majority of scales. However, these constraints are strongly dependent on even small amounts of non-Gaussianity, which is unconstrained on scales significantly smaller than those visible in the CMB. This paper goes beyond previous considerations to consider the effects of a bispectrum of the equilateral, orthogonal and local shapes with arbitrary magnitude upon the abundance of primordial black holes. Non-Gaussian density maps of the early universe are generated from a given bispectrum and used to place constraints on the small scale power spectrum. When small, we show that the skewness provides an accurate estimate for how the constraint depends on non-Gaussianity, independently of the shape of the bispectrum. We show that the orthogonal template of non-Gaussianity has an order of magnitude weaker effect on the constraints than the local and equilateral templates.

The design and performance of the LUX-ZEPLIN (LZ) detector is described as of March 2015 in this Conceptual Design Report. LZ is a second-generation dark-matter detector with the potential for unprecedented sensitivity to weakly interacting massive particles (WIMPs) of masses from a few GeV/c2 to hundreds of TeV/c2. With total liquid xenon mass of about 10 tonnes, LZ will be the most sensitive experiment for WIMPs in this mass region by the end of the decade. This report describes in detail the design of the LZ technical systems. Expected backgrounds are quantified and the performance of the experiment is presented. The LZ detector will be located at the Sanford Underground Research Facility in South Dakota. The organization of the LZ Project and a summary of the expected cost and current schedule are given.

We introduce the first phase of the Kepler-Swift Active Galaxies and Stars survey (KSwAGS), a simultaneous X-ray and UV survey of ~6 square degrees of the Kepler field using the Swift XRT and UVOT. We detect 93 unique X-ray sources with S/N>3 with the XRT, of which 60 have observed UV counterparts. We use the Kepler Input Catalog (KIC) to obtain the optical counterparts of these sources, and construct the X-ray to optical flux ratio as a first approximation of the classification of the source. The survey produces a mixture of stellar sources, extragalactic sources, and sources which we are not able to classify with certainty. We have obtained optical spectra for thirty of these targets, and are conducting an ongoing observing campaign to fully identify the sample. For sources classified as stellar or AGN with certainty, we construct SEDs using the 2MASS, UBV and GALEX data supplied for their optical counterparts by the KIC, and show that the SEDs differ qualitatively between the source types, and so can offer a method of classification in absence of a spectrum. Future papers in this series will analyze the timing properties of the stars and AGN in our sample separately. Our survey provides the first X-ray and UV data for a number of known variable stellar sources, as well as a large number of new X-ray detections in this well-studied portion of the sky. The KSwAGS survey is currently ongoing in the K2 ecliptic plane fields.

Primordial black holes (PBHs) are black holes which may have formed very early on during the radiation dominated era in the early universe. We present here a method by which the large scale perturbations in the density of primordial black holes may be used to place tight constraints on non-gaussianity if PBHs account for dark matter (DM). The presence of local-type non-gaussianity is known to have a significant effect on the abundance of primordial black holes, and modal coupling from the observed CMB scale modes can significantly alter the number density of PBHs that form within different regions of the universe, which appear as DM isocurvature modes. Using the recent \emph{Planck} constraints on isocurvature perturbations, we show that PBHs are excluded as DM candidates for even very small local-type non-gaussianity, $|f_{NL}|\approx0.001$ and remarkably the constraint on $g_{NL}$ is almost as strong. Even small non-gaussianity is excluded if DM is composed of PBHs. If local non-Gaussianity is ever detected on CMB scales, the constraints on the fraction of the universe collapsing into PBHs (which are massive enough to have not yet evaporated) will become much tighter.

The effects of non-gaussianity on the constraints on the primordial curvature perturbation power spectrum from primordial black holes (PBHs) are considered. We extend previous analyses to include the effects of coupling between the modes of the horizon scale at the time the PBH forms and super-horizon modes. We consider terms of up to third order in the Gaussian perturbation. For the weakest constraints on the abundance of PBHs in the early universe (corresponding to a fractional energy density of PBHs of $10^{-5}$ at the time of formation), in the case of gaussian perturbations, constraints on the power spectrum are $\mathcal{P}_{\zeta}<0.05$ but can significantly tighter when even a small amount of non-gaussianity is considered, to $\mathcal{P}_{\zeta}<0.01$, and become approximately $\mathcal{P}_{\zeta}<0.003$ in more special cases. Surprisingly, even when there is negative skew (which naively would suggest fewer areas of high density, leading to weaker constraints), we find that the constraints on the power spectrum become tighter than the purely gaussian case - in strong contrast with previous results. We find that the constraints are highly sensitive to both the non-gaussianity parameters as well as the amplitude of super-horizon perturbations.

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis.

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a \tau propto \lambda^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

We reinspect the calculation for the mass fraction of primordial black holes (PBHs) which are formed from primordial perturbations, finding that performing the calculation using the comoving curvature perturbation $\mathcal{R}_{c}$ in the standard way vastly overestimates the number of PBHs, by many orders of magnitude. This is because PBHs form shortly after horizon entry, meaning modes significantly larger than the PBH are unobservable and should not affect whether a PBH forms or not - this important effect is not taken into account by smoothing the distribution in the standard fashion. We discuss alternative methods and argue that the density contrast, $\Delta$, should be used instead as super-horizon modes are damped by a factor $k^{2}$. We make a comparison between using a Press-Schechter approach and peaks theory, finding that the two are in close agreement in the region of interest. We also investigate the effect of varying the spectral index, and the running of the spectral index, on the abundance of primordial black holes.

We present an analysis of 2.05 $\mu$m Hubble Space Telescope (HST) polarimetric data for a sample of 13 nearby Fanaroff-Riley type II (FRII) 3CR radio sources ($0.03<z<0.11$) that are classified as narrow line radio galaxies (NLRG) at optical wavelengths. We find that the compact cores of the NLRG in our sample are intrinsically highly polarised in the near-IR ($6 < P_{2.05\mu m} < 60$ per cent), with the electric-vector (E-vector) perpendicular to the radio axis in 54 per cent of the sources. The levels of extinction required to produce near-infrared polarisation by the dichroic extinction mechanism are consistent with the measured values reported in Ram\'irez et al. (2014), provided that this mechanism has its maximum efficiency. This consistency suggests that the nuclear polarisation could be due to dichroic extinction. In this case, toroidal magnetic fields that are highly coherent would be required in the circumnuclear tori to align the elongated dust grains responsible for the dichroic extinction. However, it is not entirely possible to rule out other polarisation mechanisms (e.g. scattering, synchrotron emission) with our observations at only one near-IR wavelength. Therefore further polarimetry observations at mid-IR and radio wavelengths will be required to test whether all the near-IR polarisation is due to dichroic extinction.

Primordial black holes (PBHs) can form in the early Universe from the collapse of rare, large density fluctuations. They have never been observed, but this fact is enough to constrain the amplitude of fluctuations on very small scales which cannot be otherwise probed. Because PBHs form only in very rare large fluctuations, the number of PBHs formed is extremely sensitive to changes in the shape of the tail of the fluctuation distribution - which depends on the amount of non-Gaussianity present. We first study how local non-Gaussianity of arbitrary size up to fifth order affects the abundance and constraints from PBHs, finding that they depend strongly on even small amounts of non-Gaussianity and the upper bound on the allowed amplitude of the power spectrum can vary by several orders of magnitude. The sign of the non-linearity parameters (f_{NL}, g_{NL}, etc) are particularly important. We also study the abundance and constraints from PBHs in the curvaton scenario, in which case the complete non-linear probability distribution is known, and find that truncating to any given order (i.e. to order f_{NL} or g_{NL}, etc) does not give accurate results.

An optically and geometrically thick torus obscures the central engine of Active Galactic Nuclei (AGN) from some lines of sight. From a magnetohydrodynamical framework, the torus can be considered to be a particular region of clouds surrounding the central engine where the clouds are dusty and optically thick. In this framework, the magnetic field plays an important role in the creation, morphology and evolution of the torus. If the dust grains within the clouds are assumed to be aligned by paramagnetic alignment, then the ratio of the intrinsic polarisation and visual extinction, P(%)/Av, is a function of the magnetic field strength. To estimate the visual extinction through the torus and constrain the polarisation mechanisms in the nucleus of AGN, we developed a polarisation model to fit both the total and polarised flux in a 1.2" (~263pc) aperture of the type 2 AGN, IC5063. The polarisation model is consistent with the nuclear polarisation observed at K being produced by dichroic absorption from aligned dust grains with a visual extinction through the torus of 48$\pm$2 mag. We estimated the intrinsic polarisation arising from dichroic absorption to be P$_{K}^{dic}$=12.5$\pm$2.7%. We consider the physical conditions and environment of the gas and dust for the torus of IC5063. Then, through paramagnetic alignment, we estimate a magnetic field strength in the range of 12-128mG in the NIR emitting regions of the torus of IC5063. Alternatively, we estimate the magnetic field strength in the plane of the sky using the Chandrasekhar-Fermi method. The minimum magnetic field strength in the plane of the sky is estimated to be 13 and 41 mG depending of the conditions within the torus of IC5063. These techniques afford the chance to make a survey of AGN, to investigate the effects of magnetic field strength on the torus, accretion, and interaction to the host galaxy.

The nuclei of Seyfert 1 galaxies exhibit a range of optical polarization characteristics that can be understood in terms of two scattering regions producing orthogonal polarizations: an extended polar scattering region (PSR) and a compact equatorial scattering region (ESR), located within the circum-nuclear torus. Here we present NICMOS 2.0 micron imaging polarimetry of 6 "polar scattered" Seyfert 1 (S1) galaxies, in which the PSR dominates the optical polarization. The unresolved nucleus (<0.58 arcsec) is significantly polarized in only three objects, but 5 of the 6 exhibit polarization in a 0.58 to 1.5 arcsec circum-nuclear annulus. In Fairall 51 and ESO 323-G077, the polarization position angle at 2 microns (theta2m) is consistent with the average for the optical spectrum (thetav), implying that the nuclear polarization is dominated by polar scattering at both wavelengths. The same is probably true for NGC 3227. In both NGC 4593 and Mrk 766, there is a large difference between theta2m and thetav off nucleus, where polar scattering is expected to dominate. This may be due to contamination by interstellar polarization in NGC 4593, but there is no clear explanation in the case of the strongly polarized Mrk 766. Lastly, in Mrk 1239, a large change (~ 60deg) in theta2m between the nucleus and the annulus indicates that the unresolved nucleus and its immediate surroundings have different polarization states at 2 microns, which we attribute to the ESR and PSR, respectively. A further implication is that the source of the scattered 2 micron emission in the unresolved nucleus is the accretion disk, rather than torus hot dust emission.

We used the CLUMPY torus models and a Bayesian approach to fit the infrared spectral energy distributions (SEDs) and ground-based high-angular resolution mid-infrared spectroscopy of 13 nearby Seyfert galaxies. This allowed us to put tight constraints on torus model parameters such as the viewing angle, the radial thickness of the torus Y, the angular size of the cloud distribution sigma_torus, and the average number of clouds along radial equatorial rays N_0. The viewing angle is not the only parameter controlling the classification of a galaxy into a type 1 or a type 2. In principle type 2s could be viewed at any viewing angle as long as there is one cloud along the line of sight. A more relevant quantity for clumpy media is the probability for an AGN photon to escape unabsorbed. In our sample, type 1s have relatively high escape probabilities, while in type 2s, as expected, tend to be low. Our fits also confirmed that the tori of Seyfert galaxies are compact with torus model radii in the range 1-6pc. The scaling of the models to the data also provided the AGN bolometric luminosities, which were found to be in good agreement with estimates from the literature. When we combined our sample of Seyfert galaxies with a sample of PG quasars from the literature to span a range of L_bol(AGN)~10^{43}-10^{47}erg/s, we found plausible evidence of the receding torus. That is, there is a tendency for the torus geometrical covering factor to be lower at high AGN luminosities than at low AGN luminosities. This is because at low AGN luminosities the tori appear to have wider angular sizes and more clouds along radial equatorial rays. We cannot, however rule out the possibility that this is due to contamination by extended dust structures not associated with the dusty torus at low AGN luminosities, since most of these in our sample are hosted in highly inclined galaxies. (Abridged)

We report spectropolarimetric observations of the quasar E1821+643 (z=0.297), which suggest that it may be an example of gravitational recoil due to anisotropic emission of gravitational waves following the merger of a supermassive black hole (SMBH) binary. In total flux, the broad Balmer lines are redshifted by ~1000 km/s relative to the narrow lines and have highly red asymmetric profiles, whereas in polarized flux the broad H_alpha line exhibits a blueshift of similar magnitude and a strong blue asymmetry. We show that these observations are consistent with a scattering model in which the broad-line region has two components, moving with different bulk velocities away from the observer and towards a scattering region at rest in the host galaxy. If the high velocity system is identified as gas bound to the SMBH, this implies that the SMBH is itself moving with a velocity ~2100 km/s relative to the host galaxy. We discuss some implications of the recoil hypothesis and also briefly consider whether our observations can be explained in terms of scattering of broad-line emission originating from the active component of an SMBH binary, or from an outflowing wind.

We report the analysis of near-infrared imaging, polarimetric and spectroscopic observations of the powerful radio galaxy 3C433, obtained with the HST and UKIRT telescopes. The high spatial resolution of HST allows us to study the near-nuclear regions of the galaxy (<1 kpc). In line with previous observations, we find that 3C433 has an unresolved core source that is detected in all near-IR bands, but dominates over the host galaxy emission at 2.05 um. Our analysis reveals: (1) the presence of a dust lane aligned close to perpendicular (PA$=70\pm5\degr$) to the inner radio jet axis (PA$=-12\pm2\degr$); (2) a steep slope to the near-IR SED ($\alpha=5.8\pm0.1$; F$_{\nu}\propto\nu^{-\alpha}$); (3) an apparent lack of broad permitted emission lines at near-IR wavelengths, in particular the absence of a broad Pa$\alpha$ emission line; and (4) high intrinsic polarization for the unresolved core nuclear source ($8.6\pm1$ per cent), with an E-vector perpendicular (PA=$83.0\pm 2.3\degr$) to the inner radio jet. Using five independent techniques we determine an extinction to the compact core source in the range 3<A_V<67 mag. An analysis of the long wavelength SED rules out a synchrotron origin for the high near-IR polarization of the compact core source. Therefore, scattering and dichroic extinction are plausible polarizing mechanisms, although in both of these cases the broad permitted lines from the AGN are required to have a width >10^4 km/s (FWHM) to escape detection in our near-IR spectrum. Dichroic extinction is the most likely polarization mechanism because it is consistent with the various available extinction estimates. In this case, a highly ordered, coherent toroidal magnetic field must be present in the obscuring structure close to the nucleus.

We present models representing the scattering of quasar radiation off free electrons and dust grains in geometries that approximate the structure of quasar host galaxies. We show that, for reasonable assumptions, scattering alone can easily produce ratios of nuclear (point source) to extended fluxes comparable to those determined in studies of quasar hosts. This result suggests that scattered quasar light, as well as stellar emission from the host galaxy, contributes significantly to the detected extended flux, leading to uncertainty in the inferred properties of quasar host. A significant contribution from scattered quasar light will lead to overestimates of the luminosity and hence mass of the host galaxy, and may also distort its morphology. Scattering of quasar light within the host galaxy may provide alternative explanations for the apparent peak in host luminosity at z = 2-3; possibly the overall average higher luminosity of radio-loud host galaxies relative to those of radio-quiet quasars (RQQs), and the apparent preference of high-luminosity RQQs for spheroidal rather than disk galaxies.

We have developed a semi-empirical simulation of the extragalactic radio continuum sky suitable for aiding the design of next generation radio interferometers such as the Square Kilometre Array (SKA). The emphasis is on modelling the large-scale cosmological distribution of radio sources rather than the internal details of individual galaxies. Here we provide a description of the simulation to accompany the online release of a catalogue of 320 million simulated radio sources. The simulation covers 20x20 deg^2 - a plausible upper limit to the instantaneous field of view attainable with future (e.g. SKA) aperture array technologies - out to redshift z=20, and down to flux density limits of 10 nJy at 151, 610 MHz, 1.4, 4.86 and 18 GHz. Five distinct source types are included: radio-quiet AGN, radio-loud AGN of the FRI and FRII structural classes, and star-forming galaxies, the latter split into populations of quiescent and starbursting galaxies. In our semi-empirical approach, the simulated sources are drawn from observed (or extrapolated) luminosity functions and grafted onto an underlying dark matter density field with biases which reflect their measured large-scale clustering. A numerical Press-Schechter-style filtering of the density field is used to identify and populate clusters of galaxies. Radio source structures are built from point source and elliptical sub-components, and for FRI and FRII sources an orientation-based unification and beaming model is used to partition flux between the core and extended lobes and hotspots. The simulation output can be post-processed to achieve more complete agreement with observational data in the years ahead, with the aim of using these 'idealised skies' in telescope simulators to optimise the design of the SKA itself (abridged).

It is now widely accepted that most galaxies undergo an active phase, during which a central super-massive black hole generates vast radiant luminosities through the gravitational accretion of gas. Winds launched from a rotating accretion disk surrounding the black hole are thought to play a critical role, allowing the disk to shed angular momentum that would otherwise inhibit accretion. Such winds are capable of depositing large amounts of mechanical energy in the host galaxy and its environs, profoundly affecting its formation and evolution, and perhaps regulating the formation of large-scale cosmological structures in the early Universe. Although there are good theoretical grounds for believing that outflows from active galactic nuclei originate as disk winds, observational verification has proven elusive. Here we show that structures observed in polarized light across the broad H-alpha emission line in the quasar PG 1700+158 originate close to the accretion disk in an electron scattering wind. The wind has large rotational motions (~4,000 km/s), providing direct observational evidence that outflows from active galactic nuclei are launched from the disks. Moreover, the wind rises nearly vertically from the disk, favouring launch mechanisms that impart an initial acceleration perpendicular to the disk plane.

We present diffraction limited, 10um imaging polarimetry data for the central regions of the archetypal Seyfert AGN, NGC1068. The position angle of polarization is consistent with three dominant polarizing mechanisms. We identify three distinct regions of polarization: (a) north of the nucleus, arising from aligned dust in the NLR, (b) south, east and west of the nucleus, consistent with dust being channeled toward the central engine and (c) a central minimum of polarization consistent with a compact (<22pc) torus. These observations provide continuity between the geometrically and optically thick torus and the host galaxy's nuclear environments. These images represent the first published mid-IR polarimetry from an 8-m class telescope and illustrate the potential of such observations.

We present results from very high signal-to-noise spectropolarimetric observations of the Seyfert 1 galaxy NGC 3783. Position Angle (PA) changes across the Balmer lines show that the scatterer is resolving the Broad-Emission Line Region (BLR). A broad component seen in polarized light and located bluewards from the H$\beta$ line very likely corresponds to HeII$\lambda4686$. The lack of PA changes across this line suggests that the region responsible for this emission appears to the scatterer as unresolved as the continuum source, in agreement with the stratified BLR structure determined from reverberation mapping.

We present detailed scattering models confirming that distinctive variations in polarization across the broad Halpha line, which are observed in a significant fraction of Seyfert 1 galaxies, can be understood in terms of a rotating line-emitting disc surrounded by a co-planar scattering region (the equatorial scattering region). The predicted polarization properties are: averaged over wavelength, the position angle of polarization is aligned with the projected disc rotation axis and hence also with the radio source axis; (ii) the polarization PA rotates across the line profile, reaching equal but opposite (relative to the continuum PA) rotations in the blue and red wings; (iii) the degree of polarization peaks in the line wings and passes through a minimum in the line core. We identify 11 objects which exhibit these features to different degrees. In order to reproduce the large amplitude PA rotations observed in some cases, the scattering region must closely surround the emission disc and the latter must itself be a relatively narrow annulus -- presumably the Halpha-emitting zone of a larger accretion disc. Asymmetries in the polarization spectra may be attributable to several possible causes, including bulk radial infall in the equatorial scattering region, or contamination by polar scattered light. The broad Halpha lines do not, in general, exhibit double-peaked profiles, suggesting that a second Halpha-emitting component of the broad-line region is present, in addition to the disc.

We have identified 12 Seyfert 1 galaxies that exhibit optical polarization spectra similar to those of Seyfert 2 galaxies in which polarized broad-lines are detected. We present new spectropolarimetric observations of 3 of them: Was 45, Mrk 231 and NGC 3227. These objects appear to be polarized as a result of far-field scattering in the polar illumination cones of the circum-nuclear torus. We estimate that they represent between 10 and 30 per cent of the Seyfert 1 population; they are found amongst all the main spectroscopic sub-types, including narrow-line Seyfert 1 galaxies. We have shown elsewhere that Seyfert 1 nuclei more commonly have polarization characteristics that can be attributed to scattering by a compact 'equatorial' scattering region located inside the torus. We propose that both equatorial and polar scattering regions are present in all Seyfert galaxies and argue that the observed range of polarization properties can be broadly understood as an orientation effect. In this scheme, polar-scattered Seyfert 1 galaxies represent the transition between unobscured (the majority of type 1) and obscured (type 2) Seyferts. They are viewed through the upper layers of the torus and are thus subject to moderate extinction (Av ~1-4 mag) sufficient to suppress polarized light from the equatorial scattering region, but not the broad wings of the Balmer lines. The orientation of the polarization position angle relative to the radio source is broadly consistent with the two-component scattering model. More generally, we find that amongst Seyferts 1 galaxies, parallel, perpendicular and intermediate orientations of the polarization PA relative to the radio axis occur roughly in the proportions 2:1:1.

We present optical spectropolarimetry of the nuclei of 36 Seyfert 1 galaxies, obtained with the William Herschel and the Anglo-Australian Telescopes from 1996 to 1999. In 20 of these, the optical emission from the active nucleus is intrinsically polarized. We have measured a significant level of polarization in a further 7 objects but these may be heavily contaminated by Galactic interstellar polarization. The intrinsically polarized Seyfert 1s exhibit a variety of characteristics, with the average polarization ranging from < 0.5 to 5 per cent and many showing variations in both the degree and position angle of polarization across the broad H alpha emission line. We identify a small group of Seyfert 1s that exhibit polarization properties similar to those of Seyfert 2 galaxies in which polarized broad-lines have been discovered. These objects represent direct observational evidence that a Seyfert 2-like far-field polar scattering region is also present in Seyfert 1s. Several other objects have features that can be explained in terms of equatorial scattering of line emission from a rotating disk. We propose that much of the diversity in the polarization properties of Seyfert galaxies can be understood in terms of a model involving both equatorial and polar scattering, the relative importance of the two geometries as sources of polarized light being determined principally by the inclination of the system axis to the line-of-sight.

We report the results of a spectropolarimetric survey of a complete far infrared selected sample of Seyfert 2 galaxies. We have found polarized broad Halpha emission in one new source, NGC5995. In the sample as a whole, there is a clear tendency for galaxies in which we have detected broad Halpha in polarized light to have warm mid--far infrared colours (F(60um)/F(25um)<4), in agreement with our previous results. However, a comparison of the optical, radio and hard x-ray properties of these systems leads us to conclude that this is a secondary consequence of the true mechanism governing our ability to see scattered light from the broad line region. We find a strong trend for galaxies showing such emission to lie above a critical value of the relative luminosity of the active core to the host galaxy (as measured from the [OIII] 5007A equivalent width) which varies as a function of the obscuring column density as measured from hard x-ray observations. The warmth of the infrared colours is then largely due to a combination of the luminosity of the active core, the obscuring column and the relative importance of the host galaxy in powering the far infrared emission, and not solely orientation as we inferred in our previous paper. Our data may also provide an explanation as to why the most highly polarized galaxies, which appear to have tori that are largely edge-on, are also the most luminous and have the most easily detectable scattered broad Halpha.

We report results from HST/NICMOS 2 micron imaging polarimetry of the central region of Centaurus A. In the vicinity of the nucleus we observe a complex polarization structure which we explain by a combination of scattering of nuclear light and dichroic polarization associated with the dust lane. The scattered nuclear radiation is found in an angular region which extends over ~ 70 degrees and thus it does not originate from a highly collimated beam, but is associated with more omni-directional nuclear illumination. These observations also show the presence of an unresolved, highly polarized (P = 11.1 %) nuclear source whose polarization angle PA = 148.2 degrees is perpendicular to the jet axis. We set an upper limit of 0.04'' (~0.8 pc) to its extent. The observed nuclear polarization is naturally accounted for if we are observing scattered light from an otherwise obscured nucleus provided that both the scattering region and the occulting torus are extremely compact, with an outer radius of less than ~ 1 pc. Alternatively, we might be directly seeing the infrared counterpart of the radio core, similar to those found in other low luminosity radio-galaxies observed with HST. We discuss these results in the framework of the FRI / BL Lac unifying model.

We present near-infrared broad-band polarization images of the nuclear regions of the Circinus galaxy in the J, H and K bands. For the first time the south-eastern reflection cone is detected in polarized light, which is obscured at optical wavelengths behind the galactic disk. This biconical structure is clearly observed in J and H band polarized flux whilst in the K band a more compact structure is detected. Total flux J-K and H-K colour maps reveal a complex colour gradient toward the south-east direction (where the Circinus galactic disk is nearer to us). We find enhanced extinction in an arc shaped structure, at about 200pc from the nucleus, probably part of the star-formation ring. We model the polarized flux images with the scattering and torus model of Young etal, with the same basic input parameters as used by Alexander etal in the spectropolarimetry modelling of Circinus. The best fit to the polarized flux is achieved with a torus radius of ~16pc, and a visual extinction Av, through the torus, to the near-infrared emission regions of >66 mags.

We present near-IR imaging polarimetry observations of the nucleus of Cygnus A, taken with the NICMOS camera of the HST at a wavelength of 2.0 microns.These maps reveal a highly collimated region of polarized emission straddling the nucleus and extending to a radius of 1.2 arcseconds. Remarkably, this feature coincides with one, but only one, limb of the edge-brightened bicone structure seen in the total intensity image. The high degree (P_k ~ 25%) and orientation of the extended polarization feature are consistent with a scattering origin. Most plausibly, the detection of polarization along only one limb of the bicone is a consequence of intrinsic anisotropy of the near-IR continuum within the radiation cones, with the direction of maximum intensity of the near-IR radiation field significantly displaced from the direction of the radio axis. The unresolved nuclear core source is also highly polarized (P_k > 28%), with a position angle close to the perpendicular to the radio axis. Given that this high degree of nuclear polarization can only be explained in terms of dichroic extinction if the dichroic mechanism is unusually efficient in Cygnus A, it is more likely that the nuclear polarization is due to the scattering of nuclear light in an unresolved scattering region close to the AGN. In this case, the flux of the core source in the K-band is dominated by scattered rather than transmitted quasar light, and previous extinction estimates based on K-band photometry of the core substantially underestimate the true nuclear extinction.

We report on X-ray, optical and infrared follow-up observations of GRB 980703. We detect a previously unknown X-ray source in the GRB error box; assuming a power law decline we find for its decay index alpha < -0.91 (3-sigma). We invoke host galaxy extinction to match the observed spectral slope with the slope expected from `fireball' models. We find no evidence for a spectral break in the infrared to X-ray spectral range on 1998 July 4.4, and determine a lower limit of the cooling break frequency: nu_c > 1.3 x 10^{17} Hz. For this epoch we obtain an extinction of A_V = 1.50 +/- 0.11. From the X-ray data we estimate the optical extinction to be A_V = 20.2 +12.3 -7.3, inconsistent with the former value. Our optical spectra confirm the redshift of z = 0.966 found by Djorgovski et al. (1998). We compare the afterglow of GRB 980703 with that of GRB 970508 and find that the fraction of the energy in the magnetic field, epsilon_B < 6x10^{-5}, is much lower in the case of GRB 980703, which is a consequence of the high frequency of the cooling break.

We report HST NICMOS observations of the central region of NGC 5128 at 2.2mic and in Paschen Alpha (PA). The continuum images show extended emission typical of an elliptical galaxy and a strong unresolved central source we identify as the nucleus of the galaxy. Its position is consistent with ground-based IR and radio data, and with the peak of reddening found with WF/PC-1. In PA, we detect a prominent elongated structure, centered on the nucleus, extended by ~2" at a position angle of 33deg, and with a major to minor axis ratio of ~2. We interpret this as an inclined, ~40 parsec diameter, thin nuclear disk of ionized gas rather than a jet-gas cloud interaction. We do see several weaker PA features, some of which may be circumnuclear gas clouds shocked by the X-ray/radio jet. The disk is one of the smallest ever observed at the nucleus of an AGN. It is not perpendicular to the jet, but consistent with being oriented along the major axis of the bulge. If it represents the warped outer portion of an accretion disk around a black hole, we conclude that even on the scale of a few parsecs, the disk is dominated by the galaxy gravitational potential and not directly related to the symmetry axis of the AGN.