XENON1T recently reported an excess of low-energy electron recoil events that may be attributable to either new physics or to the radioactive decay of tritium. It is likely that hydrogen is not be effectively removed by the hot zirconium getters deployed in the detector. Cosmogenic activation of the xenon underground is found to be insufficient to describe the observed excess, although gases diffusing out of detector materials from cosmogenic activation on surface may contribute. Changes in the operation of gas purification systems for XENON1T and other liquid nobel gas detectors could both confirm the tritium hypothesis and remove it from the detector.

This article presents an analysis and the resulting limits on light dark matter inelastically scattering off of electrons, and on dark photon and axion-like particle absorption, using a second-generation SuperCDMS high-voltage eV-resolution detector. The 0.93 gram Si detector achieved a 3 eV phonon energy resolution; for a detector bias of 100 V, this corresponds to a charge resolution of 3% of a single electron-hole pair. The energy spectrum is reported from a blind analysis with 1.2 gram-days of exposure acquired in an above-ground laboratory. With charge carrier trapping and impact ionization effects incorporated into the dark matter signal models, the dark matter-electron cross section $\bar{\sigma}_{e}$ is constrained for dark matter masses from 0.5--$10^{4} $MeV$/c^{2}$; in the mass range from 1.2--50 eV$/c^{2}$ the dark photon kinetic mixing parameter $\varepsilon$ and the axioelectric coupling constant $g_{ae}$ are constrained. The minimum 90% confidence-level upper limits within the above mentioned mass ranges are $\bar{\sigma}_{e}\,=\,8.7\times10^{-34}$ cm$^{2}$, $\varepsilon\,=\,3.3\times10^{-14}$, and $g_{ae}\,=\,1.0\times10^{-9}$.

Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $\tau^+ \tau^-$ ($b\bar{b}$) with a local density of $\rho_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $\rho_{\mathrm{DM}}$ and are independent of the DM velocity distribution.

The size and structure of the dusty circumnuclear torus in active galactic nuclei (AGN) can be investigated by analyzing the temporal response of the torus's infrared (IR) dust emission to variations in the AGN ultraviolet/optical luminosity. This method, reverberation mapping, is applicable over a wide redshift range, but the IR response is sensitive to several poorly constrained variables relating to the dust distribution and its illumination, complicating the interpretation of measured reverberation lags. We have used an enhanced version of our torus reverberation mapping code (TORMAC) to conduct a comprehensive exploration of the torus response functions at selected wavelengths, for the standard interstellar medium grain composition. The shapes of the response functions vary widely over the parameter range covered by our models, with the largest variations occurring at shorter wavelengths ($\leq 4.5\,\mu$m). The reverberation lag, quantified as the response-weighted delay (RWD), is most affected by the radial depth of the torus, the steepness of the radial cloud distribution, the degree of anisotropy of the AGN radiation field, and the volume filling factor. Nevertheless, we find that the RWD provides a reasonably robust estimate, to within a factor of $\sim 3$, of the luminosity-weighted torus radius, confirming the basic assumption underlying reverberation mapping. However, overall, the models predict radii at $2.2\,\mu$m that are a typically factor of $\sim 2$ larger than those derived from K-band reverberation mapping. This is likely an indication that the innermost region of the torus is populated by clouds dominated by large graphite grains.

In their recent preprint [arXiv:2002.06937], Kurinsky, Baxter, Kahn, and Krnjaic assume an unphysical ionization yield for plasmon excitations in order to claim a possible dark matter signal. Their proposed signal is not possible based on known physics, but their proposed detection method warrants further investigation.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

Future direct searches for low-mass dark matter particles with germanium detectors, such as SuperCDMS SNOLAB, are expected to be limited by backgrounds from radioactive isotopes activated by cosmogenic radiation inside the germanium. There are limited experimental data available to constrain production rates and a large spread of theoretical predictions. We examine the calculation of expected production rates, and analyze data from the second run of the CDMS low ionization threshold experiment (CDMSlite) to estimate the rates for several isotopes. We model the measured CDMSlite spectrum and fit for contributions from tritium and other isotopes. Using the knowledge of the detector history, these results are converted to cosmogenic production rates at sea level. The production rates in atoms/(kg$\cdot$day) are 74$\pm$9 for $^3$H, 1.5$\pm$0.7 for $^{55}$Fe, 17$\pm$5 for $^{65}$Zn, and 30$\pm$18 for $^{68}$Ge.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

We present the first limits on inelastic electron-scattering dark matter and dark photon absorption using a prototype SuperCDMS detector having a charge resolution of 0.1 electron-hole pairs (CDMS HVeV, a 0.93 gram CDMS HV device). These electron-recoil limits significantly improve experimental constraints on dark matter particles with masses as low as 1 MeV/$\mathrm{c^2}$. We demonstrate a sensitivity to dark photons competitive with other leading approaches but using substantially less exposure (0.49 gram days). These results demonstrate the scientific potential of phonon-mediated semiconductor detectors that are sensitive to single electronic excitations.

We present the molecular gas morphology and kinematics of seven nearby Seyfert galaxies obtained from our 230~GHz ALMA observations. The CO J=2-1 kinematics within the inner $\sim30$" ($\lesssim9$~kpc) reveals rotation patterns that have been explored using the Bertola rotation model and a modified version of the Kinemetry package. The latter algorithm reveals various deviations from pure circular rotation in the inner kiloparsec of all seven galaxies, including kinematic twists, decoupled and counter-rotating cores. A comparison of the global molecular gas and stellar kinematics show overall agreement in the position angle of the major axis and the systemic velocity, but larger discrepancies in the disc inclination. The residual maps obtained with both the methods shows the presence of non-circular motions in most of the galaxies. Despite its importance, a detailed interpretation of the physics responsible for non-circular motions will be discussed in a forthcoming work.

The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) searches for interactions between dark matter particles and germanium nuclei in cryogenic detectors. The experiment has achieved a low energy threshold with improved sensitivity to low-mass (<10 GeV/c$^2$) dark matter particles. We present an analysis of the final CDMSlite data set, taken with a different detector than was used for the two previous CDMSlite data sets. This analysis includes a data "salting" method to protect against bias, improved noise discrimination, background modeling, and the use of profile likelihood methods to search for a dark matter signal in the presence of backgrounds. We achieve an energy threshold of 70 eV and significantly improve the sensitivity for dark matter particles with masses between 2.5 and 10 GeV/c$^2$ compared to previous analyses. We set an upper limit on the dark matter-nucleon scattering cross section in germanium of 5.4$\times$10$^{-42}$ cm$^2$ at 5 GeV/c$^2$, a factor of $\sim$2.5 improvement over the previous CDMSlite result.

The Super Cryogenic Dark Matter Search experiment (SuperCDMS) at the Soudan Underground Laboratory studied energy loss associated with Frenkel defect formation in germanium crystals at mK temperatures using in situ $^{210}$Pb sources. We examine the spectrum of $^{206}$Pb nuclear recoils near its expected 103 keV endpoint energy and determine an energy loss of $\left(6.08\pm0.18\right)$ %, which we attribute to defect formation. From this result and using TRIM simulations, we extract the first experimentally determined average displacement threshold energy of $\left(19.7^{+0.6}_{-0.5}\right)$ eV for germanium. This has implications for the analysis thresholds of future germanium-based dark matter searches.

We present a multiwavelength study of the OH megamaser galaxy (OHMG) IRAS17526+3253, based on new Gemini Multi-Object Spectrograph Integral Field Unit (GMOS/IFU) observations, Hubble Space Telescope F814W and H$\alpha$+[N{\sc ii}] images, and archival 2MASS and 1.49GHz VLA data. The HST images clearly reveal a mid-to-advanced stage major merger whose northwestern and southeastern nuclei have a projected separation of $\sim$8.5kpc. Our HST/H$\alpha$+[N{\sc ii}] image shows regions of ongoing star-formation across the envelope on $\sim$10kpc scales, which are aligned with radio features, supporting the interpretation that the radio emission originates from star-forming regions. The measured H$\alpha$ luminosities imply that the unobscured star-formation rate is $\sim$10-30\,M$_{\odot}$yr$^{-1}$. The GMOS/IFU data reveal two structures in northwestern separated by 850\,pc and by a discontinuity in the velocity field of $\sim$~200~km~s$^{-1}$. We associate the blue-shifted and red-shifted components with, respectively, the distorted disk of northwestern and tidal debris, possibly a tail originating in southeastern. Star-formation is the main ionization source in both components, which have SFRs of $\sim$2.6-7.9\,M$_{\odot}$yr$^{-1}$ and $\sim$1.5-4.5\,M$_{\odot}$yr$^{-1}$, respectively. Fainter line emission bordering these main components is consistent with shock ionization at a velocity $\sim$200~km~s$^{-1}$ and may be the result of an interaction between the tidal tail and the northwestern galaxy's disk. IRAS17526+3253 is one of only a few systems known to host both luminous OH and H$_{2}$O masers. The velocities of the OH and H$_{2}$O maser lines suggest that they are associated with the northwestern and southeastern galaxies, respectively.

Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total background rate as the previous 1167-kg-day exposure at 3.3 keV. This increased exposure is enabled in part by a new optical tracking analysis to better identify events near detector walls, permitting a larger fiducial volume. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 2.5 $\times$ 10$^{-41}$ cm$^2$ for a 25 GeV WIMP, and improve on previous PICO results for 3-5 GeV WIMPs by an order of magnitude.

The Cryogenic Dark Matter Search (CDMS II) experiment aims to detect dark matter particles that elastically scatter from nuclei in semiconductor detectors. The resulting nuclear-recoil energy depositions are detected by ionization and phonon sensors. Neutrons produce a similar spectrum of low-energy nuclear recoils in such detectors, while most other backgrounds produce electron recoils. The absolute energy scale for nuclear recoils is necessary to interpret results correctly. The energy scale can be determined in CDMS II silicon detectors using neutrons incident from a broad-spectrum $^{252}$Cf source, taking advantage of a prominent resonance in the neutron elastic scattering cross section of silicon at a recoil (neutron) energy near 20 (182) keV. Results indicate that the phonon collection efficiency for nuclear recoils is $4.8^{+0.7}_{-0.9}$% lower than for electron recoils of the same energy. Comparisons of the ionization signals for nuclear recoils to those measured previously by other groups at higher electric fields indicate that the ionization collection efficiency for CDMS II silicon detectors operated at $\sim$4 V/cm is consistent with 100% for nuclear recoils below 20 keV and gradually decreases for larger energies to $\sim$75% at 100 keV. The impact of these measurements on previously published CDMS II silicon results is small.

We have obtained observations of the OH Megamaser galaxy IRAS03056+0234 using Gemini Multi-Object Spectrograph (GMOS) Integral Field Unit (IFU), Very Large Array (VLA) and Hubble Space Telescope (HST). The HST data reveals spiral arms containing knots of emission associated to star forming regions. The GMOS-IFU data cover the spectral range of 4500 to 7500 \AA\ at a velocity resolution of 90 km s$^{-1}$ and spatial resolution of 506 pc. The emission-line flux distributions reveal a ring of star forming regions with radius of 786 pc centred at the nucleus of the galaxy, with an ionized gas mass of 1.2$\times$ 10$^{8}$M$_{\odot}$, an ionizing photon luminosity of log Q[H$^{+}$]=53.8 and a star formation rate of 4.9 M$_{\odot}$ yr$^{-1}$. The emission-line ratios and radio emission suggest that the gas at the nuclear region is excited by both starburst activity and an active galactic nucleus. The gas velocity fields are partially reproduced by rotation in the galactic plane, but show, in addition, excess redshifts to the east of the nucleus, consistent with gas inflows towards the nucleus, with velocity of $\sim$45 km s$^{-1}$ and a mass inflow rate of $\sim$7.7$\times$10$^{-3}$ M$_{\odot}$ yr$^{-1}$.

We aim to map the distribution and kinematics of molecular and ionized gas in a sample of active galaxies, to quantify the nuclear inflows and outflows. Here, we analyze the nuclear kinematics of NGC 1566 via ALMA observations of the CO J:2-1 emission at 24 pc spatial and $\sim$2.6 km s$^{-1}$ spectral resolution, and Gemini-GMOS/IFU observations of ionized gas emission lines and stellar absorption lines at similar spatial resolution, and 123 km s$^{-1}$ of intrinsic spectral resolution. The morphology and kinematics of stellar, molecular (CO) and ionized ([N II]) emission lines are compared to the expectations from rotation, outflows, and streaming inflows. While both ionized and molecular gas show rotation signatures, there are significant non-circular motions in the innermost 200 pc and along spiral arms in the central kpc (CO). The nucleus shows a double-peaked CO profile (Full Width at Zero Intensity of 200 km s$^{-1}$), and prominent ($\sim$80 km s$^{-1}$) blue and redshifted lobes are found along the minor axis in the inner arcseconds. Perturbations by the large-scale bar can qualitatively explain all features in the observed velocity field. We thus favour the presence of a molecular outflow in the disk with true velocities of $\sim$180 km s$^{-1}$ in the nucleus and decelerating to 0 by $\sim$72 pc. The implied molecular outflow rate is $5.6~[M_{o}yr^{-1}]$, with this gas accumulating in the nuclear 2 arcsec arms. The ionized gas kinematics support an interpretation of a similar, but more spherical, outflow in the inner 100 pc, with no signs of deceleration. There is some evidence of streaming inflows of $\sim$50 km s$^{-1}$ along specific spiral arms, and the estimated molecular mass inflow rate, $\sim0.1~[M_{o}yr^{-1}]$, is significantly larger than the SMBH accretion rate ($\dot{m}=4.8\times10^{-5}~[M_{o}yr^{-1}]$).

We present two-dimensional stellar and gaseous kinematics of the inner 0.7 $\times$ 1.2 kpc$^{2}$ of the Seyfert galaxy ESO 362-G18, derived from optical spectra obtained with the GMOS/IFU on the Gemini South telescope at a spatial resolution of $\approx$170 pc and spectral resolution of 36 km s$^{-1}$. ESO 362-G18 is a strongly perturbed galaxy of morphological type Sa or S0/a, with a minor merger approaching along the NE direction. Previous studies have shown that the [OIII] emission shows a fan-shaped extension of $\approx$ 10\arcsec\ to the SE. We detect the [OIII] doublet, [NII] and H${\alpha}$ emission lines throughout our field of view. The stellar kinematics is dominated by circular motions in the galaxy plane, with a kinematic position angle of $\approx$137$^{\circ}$. The gas kinematics is also dominated by rotation, with kinematic position angles ranging from 122$^{\circ}$ to 139$^{\circ}$. A double-Gaussian fit to the [OIII]$\lambda$5007 and H${\alpha}$ lines, which have the highest signal to noise ratios of the emission lines, reveal two kinematic components: (1) a component at lower radial velocities which we interpret as gas rotating in the galactic disk; and (2) a component with line of sight velocities 100-250 km s$^{-1}$ higher than the systemic velocity, interpreted as originating in the outflowing gas within the AGN ionization cone. We estimate a mass outflow rate of 7.4 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$ in the SE ionization cone (this rate doubles if we assume a biconical configuration), and a mass accretion rate on the supermassive black hole (SMBH) of 2.2 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$. The total ionized gas mass within $\sim$84 pc of the nucleus is 3.3 $\times$ 10$^{5}$ M$_{\odot}$; infall velocities of $\sim$34 km s$^{-1}$ in this gas would be required to feed both the outflow and SMBH accretion.

We present two-dimensional maps of emission-line fluxes and kinematics, as well as of the stellar kinematics of the central few kpc of five bright nearby Seyfert galaxies -- Mrk\,6, Mrk\,79, Mrk\,348, Mrk\,607 and Mrk\,1058 -- obtained from observations with the Gemini Multi-Object Spectrograph (GMOS) Integral Field Unit (IFU) on the Gemini North Telescope. The data cover the inner 3\farcs5$\times$5\farcs0 -- corresponding to physical scales in the range 0.6$\times$0.9 to 1.5$\times$2.2\,kpc$^2$ -- at a spatial resolution ranging from 110 to 280 pc with a spectral coverage of 4300 -- 7100\,\AA\ and velocity resolution of $\approx$ 90\,km\,s$^{-1}$. The gas excitation is Seyfert like everywhere but show excitation, but show excitation gradients that are correlated with the gas kinematics, reddening and/or the gas density. The gas kinematics show in all cases two components: a rotation one similar to that observed in the stellar velocity field, and an outflow component. In the case of Mrk607, the gas is counter-rotating relative to the stars. Enhanced gas velocity dispersion is observed in association to the outflows according to two patterns: at the locations of the highest outflow velocities along the ionization axis or perpendicularly to it in a strip centered at the nucleus that we attribute to an equatorial outflow. Bipolar outflows are observed in Mrk\,348 and Mrk\,79, while in Mrk\,1058 only the blueshifted part is clearly observed, while in the cases of Mrk\,6 and Mrk\,607 the geometry of the outflow needs further constraints from modeling to be presented in a forthcoming study, where the mass flow rate and powers will also be obtained.

The SuperCDMS experiment is designed to directly detect weakly interacting massive particles (WIMPs) that may constitute the dark matter in our Galaxy. During its operation at the Soudan Underground Laboratory, germanium detectors were run in the CDMSlite mode to gather data sets with sensitivity specifically for WIMPs with masses ${<}$10 GeV/$c^2$. In this mode, a higher detector-bias voltage is applied to amplify the phonon signals produced by drifting charges. This paper presents studies of the experimental noise and its effect on the achievable energy threshold, which is demonstrated to be as low as 56 eV$_{\text{ee}}$ (electron equivalent energy). The detector-biasing configuration is described in detail, with analysis corrections for voltage variations to the level of a few percent. Detailed studies of the electric-field geometry, and the resulting successful development of a fiducial parameter, eliminate poorly measured events, yielding an energy resolution ranging from ${\sim}$9 eV$_{\text{ee}}$ at 0 keV to 101 eV$_{\text{ee}}$ at ${\sim}$10 eV$_{\text{ee}}$. New results are derived for astrophysical uncertainties relevant to the WIMP-search limits, specifically examining how they are affected by variations in the most probable WIMP velocity and the Galactic escape velocity. These variations become more important for WIMP masses below 10 GeV/$c^2$. Finally, new limits on spin-dependent low-mass WIMP-nucleon interactions are derived, with new parameter space excluded for WIMP masses $\lesssim$3 GeV/$c^2$

We present Gemini Multi-Object Spectrograph (GMOS) Integral field Unit (IFU), Very Large Array (VLA) and Hubble Space Telescope (HST) observations of the OH Megamaser (OHM) galaxy IRASF23199+0123. Our observations show that this system is an interacting pair, with two OHM sources associated to the eastern (IRAS23199E) member. The two members of the pair present somewhat extended radio emission at 3 and 20~cm, with flux peaks at each nucleus. The GMOS-IFU observations cover the inner $\sim$6kpc of IRAS23199E at a spatial resolution of 2.3~kpc. The GMOS-IFU flux distributions in H$\alpha$ and [NII]$\lambda$6583 are similar to that of an HST [NII]+H$\alpha$ narrow-band image, being more extended along the northeast-southwest direction, as also observed in the continuum HST F814W image. The GMOS-IFU H$\alpha$ flux map of IRAS23199E shows three extranuclear knots attributed to star-forming complexes. We have discovered a Seyfert 1 nucleus in this galaxy, as its nuclear spectrum shows an unresolved broad (FWHM$\approx$2170 kms$^{-1}$) double-peaked H$\alpha$ component, from which we derive a black hole mass of M$_{BH}$= 3.8$^{+0.3}_{-0.2}\times 10^{6}$M$_{\odot}$. The gas kinematics shows low velocity dispersions ($\sigma$) and low [NII]/H$\alpha$ ratios for the star-forming complexes and higher $\sigma$ and [NII]/H$\alpha$ surrounding the radio emission region, supporting interaction between the radio-plasma and ambient gas. The two OH masers detected in IRASF23199E are observed in the vicinity of these enhanced $\sigma$ regions, supporting their association with the active nucleus and its interaction with the surrounding gas. The gas velocity field can be partially reproduced by rotation in a disk, with residuals along the north-south direction being tentatively attributed to emission from the front walls of a bipolar outflow.

The obscuring circumnuclear torus of dusty molecular gas is one of the major components of active galactic nuclei (AGN). The torus can be studied by analyzing the time response of its infrared (IR) dust emission to variations in the AGN continuum luminosity, a technique known as reverberation mapping. The IR response is the convolution of the AGN ultraviolet/optical light curve with a transfer function that contains information about the size, geometry, and structure of the torus. Here, we describe a new computer model that simulates the reverberation response of a clumpy torus. Given an input optical light curve, the code computes the emission of a 3D ensemble of dust clouds as a function of time at selected IR wavelengths, taking into account light travel delays. We present simulated dust emission responses at 3.6, 4.5, and 30 $\mu$m that explore the effects of various geometrical and structural properties, dust cloud orientation, and anisotropy of the illuminating radiation field. We also briefly explore the effects of cloud shadowing (clouds are shielded from the AGN continuum source). Example synthetic light curves have also been generated, using the observed optical light curve of the Seyfert 1 galaxy NGC 6418 as the input. The torus response is strongly wavelength-dependent, due to the gradient in cloud surface temperature within the torus, and because the cloud emission is strongly anisotropic at shorter wavelengths. Anisotropic illumination of the torus also significantly modifies the torus response, reducing the lag between the IR and optical variations.

We use optical spectra from the inner 1.8 $\times$ 2.5kpc$^2$ of the Seyfert 2 galaxy NGC1358, obtained with the GMOS integral field spectrograph on the Gemini South telescope at a spatial resolution of $\approx$ 165pc, to assess the feeding and feedback processes in this nearby active galaxy. Five gaseous kinematical components are observed in the emission line profiles. One of the components is present in the entire field-of-view and we interpret it as due to gas rotating in the disk of the galaxy. Three of the remaining components we interpret as associated to active galactic nucleus (AGN) feedback: a compact unresolved outflow in the inner 1 arcsec and two gas clouds observed at opposite sides of the nucleus, which we propose have been ejected in a previous AGN burst. The disk component velocity field is strongly disturbed by a large scale bar. The subtraction of a velocity model combining both rotation and bar flows reveals three kinematic nuclear spiral arms: two in inflow and one in outflow. We estimate the mass inflow rate in the inner 180pc obtaining $\dot{M}_{in}$ $\approx$ 1.5 $\times 10^{-2}$M$_{\odot}$yr$^{-1}$, about 160 times larger than the accretion rate necessary to power this AGN.

New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results.

We used the Gemini Multi-Object Spectrograph Integral Field Unit to map the gas distribution, excitation and kinematics within the inner kiloparsec of four nearby low-luminosity active galaxies: NGC3982, NGC4501, NGC2787 and NGC4450. The observations cover the spectral range 5600-7000{\AA} at a velocity resolution of 120km/s and spatial resolution ranging from 50 to 70pc at the galaxies. Extended emission in H{\alpha}, [NII]{\lambda}{\lambda}6548,6583, [SII]{\lambda}{\lambda}6716,6730 over most of the field-of-view is observed for all galaxies, while only NGC3982 shows [OI]{\lambda}6300 extended emission. The H{\alpha} equivalent widths combined with the [NII]/H{\alpha} line ratios reveal that NGC3982 and NGC4450 harbor Seyfert nuclei surrounded by regions with LINER excitation, while NGC2787 and NGC4501 harbor LINER nuclei. NGC3982 shows a partial ring of recent star-formation at 500pc from the nucleus, while in NGC4501 a region at 500pc west of the nucleus shows LINER excitation but has been interpreted as an aging HII region with the gas excitation dominated by shocks from supernovae. The line-of-sight velocity field of the gas shows a rotation pattern for all galaxies, with deviations from pure disk rotation observed in NGC3982, NGC4501 and NGC4450. For NGC4501 and NGC4450, many of these deviations are spatially coincident with dust structures seen in optical continuum images, leading to the interpretation that the deviations are due to shocks in the gas traced by the dust. A speculation is that these shocks lead to loss of angular momentum, allowing the gas to be transferred inwards to feed the AGN. In the case of NGC2787, instead of deviations in the rotation field, we see a misalignment of 40{^\circ} between the orientation of the line of nodes of the gas rotation and the photometric major axis of the galaxy. Evidence of compact nuclear outflows are seen in NGC4501 and NGC4450.

Proposed dark matter detectors with eV-scale sensitivities will detect a large background of atomic (nuclear) recoils from coherent photon scattering. This background climbs steeply below $\sim10$~eVnr, far exceeding the declining rate of low-energy Compton recoils. The upcoming generation of dark matter detectors will not be limited by this background, but further development of eV-scale and sub-eV detectors will require the use of low-$Z$ target materials, such as helium, to avoid a large rate of coherent photon scattering, or highly efficient methods to reject photons when they scatter.

We report a measurement of the ionization efficiency of silicon nuclei recoiling with sub-keV kinetic energy in the bulk silicon of a charge-coupled device (CCD). Nuclear recoils are produced by low-energy neutrons ($<$24 keV) from a $^{124}$Sb-$^{9}$Be photoneutron source, and their ionization signal is measured down to 60 eV electron equivalent. This energy range, previously unexplored, is relevant for the detection of low-mass dark matter particles. The measured efficiency is found to deviate from the extrapolation to low energies of the Lindhard model. This measurement also demonstrates the sensitivity to nuclear recoils of CCDs employed by DAMIC, a dark matter direct detection experiment located in the SNOLAB underground laboratory.

SuperCDMS SNOLAB will be a next-generation experiment aimed at directly detecting low-mass (< 10 GeV/c$^2$) particles that may constitute dark matter by using cryogenic detectors of two types (HV and iZIP) and two target materials (germanium and silicon). The experiment is being designed with an initial sensitivity to nuclear recoil cross sections ~ 1 x 10$^{-43}$ cm$^2$ for a dark matter particle mass of 1 GeV/c$^2$, and with capacity to continue exploration to both smaller masses and better sensitivities. The phonon sensitivity of the HV detectors will be sufficient to detect nuclear recoils from sub-GeV dark matter. A detailed calibration of the detector response to low energy recoils will be needed to optimize running conditions of the HV detectors and to interpret their data for dark matter searches. Low-activity shielding, and the depth of SNOLAB, will reduce most backgrounds, but cosmogenically produced $^{3}$H and naturally occurring $^{32}$Si will be present in the detectors at some level. Even if these backgrounds are x10 higher than expected, the science reach of the HV detectors would be over three orders of magnitude beyond current results for a dark matter mass of 1 GeV/c$^2$. The iZIP detectors are relatively insensitive to variations in detector response and backgrounds, and will provide better sensitivity for dark matter particle masses (> 5 GeV/c$^2$). The mix of detector types (HV and iZIP), and targets (germanium and silicon), planned for the experiment, as well as flexibility in how the detectors are operated, will allow us to maximize the low-mass reach, and understand the backgrounds that the experiment will encounter. Upgrades to the experiment, perhaps with a variety of ultra-low-background cryogenic detectors, will extend dark matter sensitivity down to the "neutrino floor", where coherent scatters of solar neutrinos become a limiting background.

This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

We employ an $^{88}$Y/Be photoneutron source to derive the quenching factor for neutron-induced nuclear recoils in germanium, probing recoil energies from a few hundred eV$_{nr}$ to 8.5keV$_{nr}$. A comprehensive Monte Carlo simulation of our setup is compared to experimental data employing a Lindhard model with a free electronic energy loss $k$ and an adiabatic correction for sub-keV$_{nr}$ nuclear recoils. The best fit $k=0.179\pm 0.001$ obtained using a Monte Carlo Markov Chain (MCMC) ensemble sampler is in good agreement with previous measurements, confirming the adequacy of the Lindhard model to describe the stopping of few-keV ions in germanium crystals at a temperature of $\sim$77 K. This value of $k$ corresponds to a quenching factor of 13.7 % to 25.3 % for nuclear recoil energies between 0.3 keV$_{nr}$ and 8.5 keV$_{nr}$, respectively.

The $^9$Be$(\gamma,n)^8$Be reaction is enhanced by a near threshold $1/2^+$ state. Contradictions between existing measurements of this reaction cross-section affect calculations of astrophysical r-process yields, dark matter detector calibrations, and the theory of the nuclear structure of $^9$Be. Select well-documented radioisotope $^9$Be$(\gamma,n)$ source yield measurements have been reanalyzed, providing a set of high-accuracy independently measured cross sections. A Breit-Wigner fit of these corrected measurements yields $E_R=1738.8\pm1.9$ keV, $\Gamma_\gamma=0.771\pm0.021$ eV, and $\Gamma_n=268\pm15$ keV for the $1/2^+$ state. A virtual $1/2^+$ state is excluded with 99.3\% confidence.

We observed the nuclear region of the galaxy NGC 1365 with the integral field unit of the Gemini Multi Object Spectrograph mounted on the GEMINI-South telescope. The field of view covers $13^{\prime\prime} \times 6^{\prime\prime}$ ($1173 \times 541$ pc$^{2}$) centered on the nucleus, at a spatial resolution of $52$ pc. The spectral coverage extends from $5600$ \AA\ to $7000$ \AA, at a spectral resolution $R=1918$. NGC 1365 hosts a Seyfert 1.8 nucleus, and exhibits a prominent bar extending out to $100^{\prime\prime}$ (9 kpc) from the nucleus. The field of view lies within the inner Lindblad resonance. Within this region, we found that the kinematics of the ionized gas (as traced by [OI], [NII], H$\alpha$, and [SII]) is consistent with rotation in the large-scale plane of the galaxy. While rotation dominates the kinematics, there is also evidence for a fan-shaped outflow, as found in other studies based on the [OIII] emission lines. Although evidence for gas inflowing along nuclear spirals has been found in a few barred galaxies, we find no obvious signs of such features in the inner kiloparsec of NGC 1365. However, the emission lines exhibit a puzzling asymmetry that could originate from gas which is slower than the gas responsible for the bulk of the narrow-line emission. We speculate that it could be tracing gas which lost angular momentum, and is slowly migrating from the inner Lindblad resonance towards the nucleus of the galaxy.

New data are reported from a second run of the 2-liter PICO-2L C$_3$F$_8$ bubble chamber with a total exposure of 129$\,$kg-days at a thermodynamic threshold energy of 3.3$\,$keV. These data show that measures taken to control particulate contamination in the superheated fluid resulted in the absence of the anomalous background events observed in the first run of this bubble chamber. One single nuclear-recoil event was observed in the data, consistent both with the predicted background rate from neutrons and with the observed rate of unambiguous multiple-bubble neutron scattering events. The chamber exhibits the same excellent electron-recoil and alpha decay rejection as was previously reported. These data provide the most stringent direct detection constraints on weakly interacting massive particle (WIMP)-proton spin-dependent scattering to date for WIMP masses $<$ 50$\,$GeV/c$^2$.

New data are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 36.8 kg of CF$_3$I and located in the SNOLAB underground laboratory. PICO-60 is the largest bubble chamber to search for dark matter to date. With an analyzed exposure of 92.8 livedays, PICO-60 exhibits the same excellent background rejection observed in smaller bubble chambers. Alpha decays in PICO-60 exhibit frequency-dependent acoustic calorimetry, similar but not identical to that reported recently in a C$_3$F$_8$ bubble chamber. PICO-60 also observes a large population of unknown background events, exhibiting acoustic, spatial, and timing behaviors inconsistent with those expected from a dark matter signal. These behaviors allow for analysis cuts to remove all background events while retaining $48.2\%$ of the exposure. Stringent limits on weakly interacting massive particles interacting via spin-dependent proton and spin-independent processes are set, and most interpretations of the DAMA/LIBRA modulation signal as dark matter interacting with iodine nuclei are ruled out.

We present optical integral field spectroscopy of the inner $2.5 \times 3.4$ kpc$^2$ of the broad-line radio galaxy Pictor A, at a spatial resolution of $\approx 400$ pc. Line emission is observed over the whole field-of-view, being strongest at the nucleus and in an elongated linear feature (ELF) crossing the nucleus from the south-west to the north-east along PA $\sim 70^\circ$. Although the broad double-peaked H$\alpha$ line and the [OI]6300/H$\alpha$ and [SII]6717+31/H$\alpha$ ratios are typical of AGNs, the [NII]6584/H$\alpha$ ratio (0.15 - 0.25) is unusually low. We suggest that this is due to the unusually low metallicity of the gas. Centroid velocity maps show mostly blueshifts to the south and redshifts to the north of the nucleus, but the velocity field is not well fitted by a rotation model. Velocity dispersions are low (< 100 km s$^{-1}$) along the ELF, ruling out a jet-cloud interaction as the origin of this structure. The ELF shows both blueshifts and redshifts in channel maps, suggesting that it is close to the plane of the sky. The ELF is evidently photoionized by the AGN, but its kinematics and inferred low metallicity suggest that this structure may have originated in a past merger event with another galaxy. We suggest that the gas acquired in this interaction may be feeding the ELF.

Several studies indicate that radio-loud (RL) Active Galactic Nuclei (AGN) are produced only by the most massive black holes (BH), $M_{\rm BH} \sim 10^8$-$10^{10} M_\odot$. This idea has been challenged by the discovery of RL Narrow Line Seyfert 1 (RL NLSy1), having estimated masses of $M_{\rm BH}$$\sim$$10^6$-$10^7$ M$_\odot$. However, these low $M_{\rm BH}$ estimates might be due to projection effects. Spectropolarimetry allows us to test this possibility by looking at RL NLSy1s under a different perspective, i.e., from the viewing angle of the scattering material. We here report the results of a pilot study of VLT spectropolarimetric observations of the RL NLSy1 PKS 2004-447. Its polarization properties are remarkably well reproduced by models in which the scattering occurs in an equatorial structure surrounding its broad line region, seen close to face-on. In particular, we detect a polarized H$\alpha$ line with a width of $\sim$ 9,000 km s$^{-1}$, $\sim 6$ times broader than the width seen in direct light. This corresponds to a revised estimate of $M_{\rm BH}$$\sim$$6\times10^8$ M$_\odot$, well within the typical range of RL AGN. The double-peaked polarized broad H$\alpha$ profile of the target suggests that the rare combination of the orientation effects and a broad line region dominated by the rotation might account for this class of objects, casting doubts on the virial estimates of BH mass for type-I AGN.

We present two-dimensional gaseous kinematics of the inner 1.2 $\times$ 1.8 kpc$^2$ of the Seyfert 2 galaxy NGC3081, from optical spectra (5600--7000\r{A}) obtained with the GMOS integral field spectrograph on the Gemini North telescope at a spatial resolution of $\approx$ 100pc. We have identified two-components in the line emitting gas. A narrower component (FWHM $\approx$ 60-100km s$^{-1}$), which appears to be gas in the galaxy disk, and which shows a distorted rotation pattern, is observed over the whole field of view. A broader component (FWHM $\approx$150-250 km s$^{-1}$) is present in the inner $\approx$ 2arcsec (200pc) and shows blueshifts and redshifts in the near and far sides of the galaxy, respectively, consistent with a bipolar outflow. Assuming this to be the case, we estimate that the mass outflow rate in ionized gas ($\dot{M}_{out}$) is between 1.9 $\times 10^{-3}$M$_{\odot}$ yr$^{-1}$ and 6.9 $\times 10^{-3}$M$_{\odot}$ yr$^{-1}$. The subtraction of a rotation model from the narrower component velocity field reveals a pattern of excess blueshifts of $\approx$ 50km s$^{-1}$ in the far side of the galaxy and similar excess redshifts in the near side, which are cospatial with a previously known nuclear bar. We interpret these residuals as due to gas following non-circular orbits in the barred potential. Under the assumption that these motions may lead to gas inflows, we estimate an upper limit for the mass inflow rate in ionized gas of $\phi$ $\approx$ 1.3 $\times 10^{-2}$M$_{\odot}$ yr$^{-1}$.

New data are reported from the operation of a 2-liter C$_3$F$_8$ bubble chamber in the 2100 meter deep SNOLAB underground laboratory, with a total exposure of 211.5 kg-days at four different recoil energy thresholds ranging from 3.2 keV to 8.1 keV. These data show that C3F8 provides excellent electron recoil and alpha rejection capabilities at very low thresholds, including the first observation of a dependence of acoustic signal on alpha energy. Twelve single nuclear recoil event candidates were observed during the run. The candidate events exhibit timing characteristics that are not consistent with the hypothesis of a uniform time distribution, and no evidence for a dark matter signal is claimed. These data provide the most sensitive direct detection constraints on WIMP-proton spin-dependent scattering to date, with significant sensitivity at low WIMP masses for spin-independent WIMP-nucleon scattering.

We present optical integral field spectroscopy of the circum-nuclear gas of the Seyfert 2 galaxy NGC 1386. The data cover the central 7$^{\prime\prime} \times 9^{\prime\prime}$ (530 $\times$ 680 pc) at a spatial resolution of 0.9" (68 pc), and the spectral range 5700-7000 \AA\ at a resolution of 66 km s$^{-1}$. The line emission is dominated by a bright central component, with two lobes extending $\approx$ 3$^{\prime\prime}$ north and south of the nucleus. We identify three main kinematic components. The first has low velocity dispersion ($\bar \sigma \approx $ 90 km s$^{-1}$), extends over the whole field-of-view, and has a velocity field consistent with gas rotating in the galaxy disk. We interpret the lobes as resulting from photoionization of disk gas in regions where the AGN radiation cones intercept the disk. The second has higher velocity dispersion ($\bar \sigma \approx$ 200 km s$^{-1}$) and is observed in the inner 150 pc around the continuum peak. This component is double peaked, with redshifted and blueshifted components separated by $\approx$ 500 km s$^{-1}$. Together with previous HST imaging, these features suggest the presence of a bipolar outflow for which we estimate a mass outflow rate of $\mathrm{\dot M} \gtrsim $ 0.1 M$_{\odot}$ yr$^{-1}$. The third component is revealed by velocity residuals associated with enhanced velocity dispersion and suggests that outflow and/or rotation is occurring approximately in the equatorial plane of the torus. A second system of velocity residuals may indicate the presence of streaming motions along dusty spirals in the disk.

We present results from a fifteen-month campaign of high-cadence (~ 3 days) mid-infrared Spitzer and optical (B and V ) monitoring of the Seyfert 1 galaxy NGC 6418, with the objective of determining the characteristic size of the dusty torus in this active galactic nucleus (AGN). We find that the 3.6 $\mu$m and 4.5 $\mu$m flux variations lag behind those of the optical continuum by $37.2^{+2.4}_{-2.2}$ days and $47.1^{+3.1}_{-3.1}$ days, respectively. We report a cross-correlation time lag between the 4.5 $\mu$m and 3.6 $\mu$m flux of $13.9^{+0.5}_{-0.1}$ days. The lags indicate that the dust emitting at 3.6 $\mu$m and 4.5 $\mu$m is located at a distance of approximately 1 light-month (~ 0.03 pc) from the source of the AGN UV-optical continuum. The reverberation radii are consistent with the inferred lower limit to the sublimation radius for pure graphite grains at 1800 K, but smaller by a factor of ~ 2 than the corresponding lower limit for silicate grains; this is similar to what has been found for near-infrared (K-band) lags in other AGN. The 3.6 and 4.5 $\mu$m reverberation radii fall above the K-band $\tau \propto L^{0.5}$ size-luminosity relationship by factors $\lesssim 2.7$ and $\lesssim 3.4$, respectively, while the 4.5 $\mu$m reverberation radius is only 27% larger than the 3.6 $\mu$m radius. This is broadly consistent with clumpy torus models, in which individual optically thick clouds emit strongly over a broad wavelength range.

A bubble chamber has been developed to be used as an active target system for low energy nuclear astrophysics experiments. Adopting ideas from dark matter detection with superheated liquids, a detector system compatible with gamma-ray beams has been developed. This detector alleviates some of the limitations encountered in standard measurements of the minute cross sections of interest to stellar environments. While the astrophysically relevant nuclear reaction processes at hydrostatic burning temperatures are dominated by radiative captures, in this experimental scheme we measure the time-reversed processes. Such photodisintegrations allow us to compute the radiative capture cross sections when transitions to excited states of the reaction products are negligible. Due to the transformation of phase space, the photodisintegration cross sections are up to two orders of magnitude higher. The main advantage of the new target-detector system is a density several orders of magnitude higher than conventional gas targets. Also, the detector is virtually insensitive to the gamma-ray beam itself, thus allowing us to detect only the products of the nuclear reaction of interest. The development and the operation as well as the advantages and disadvantages of the bubble chamber are discussed.

We present a multiwavelength study of the OH Megamaser galaxy (OHMG) IRAS16399-0937, based on new HST/ACS F814W and H$\alpha$+[NII] images and archive data from HST, 2MASS, Spitzer, Herschel and the VLA. This system has a double nucleus, whose northern (IRAS16399N) and southern (IRAS16399S) components have a projected separation of $\sim$ 6'' (3.4 kpc) and have previously been identified based on optical spectra as a Low Ionization Nuclear Emission Line Region (LINER) and starburst nucleus, respectively. The nuclei are embedded in a tidally distorted common envelope, in which star formation is mostly heavily obscured. The infrared spectrum is dominated by strong polycyclic aromatic hydrocarbon (PAH), but deep silicate and molecular absorption features are also present, and are strongest in the IRAS16399N nucleus. The 0.435 - 500$\mu$m SED was fitted with a model including stellar, ISM and AGN torus components using our new MCMC code, clumpyDREAM. The results indicate that the IRAS16399N contains an AGN (L$_{bol} \sim 10^{44}$ ergs/s) deeply embedded in a quasi-spherical distribution of optically-thick clumps with a covering fraction $\approx1$. We suggest that these clumps are the source of the OHM emission in IRAS16399-0937. The high torus covering fraction precludes AGN-photoionization as the origin of the LINER spectrum, however, the spectrum is consistent with shocks (v $\sim100-200$ km s$^{-1}$). We infer that the $\sim10^8$ M$_{\odot}$ black-hole in IRAS16399N is accreting at a small fraction ($\sim1$%) of its Eddington rate. The low accretion-rate and modest nuclear SFRs suggest that while the gas-rich major merger forming the IRAS16399-0937 system has triggered widespread star formation, the massive gas inflows expected from merger simulations have not yet fully developed.

The coalescence of a binary black hole can be accompanied by a large gravitational recoil due to anisotropic emission of gravitational waves. A recoiling supermassive black hole (SBH) can subsequently undergo long-lived oscillations in the potential well of its host galaxy, suggesting that offset SBHs may be common in the cores of massive ellipticals. We have analyzed HST archival images of 14 nearby core ellipticals, finding evidence for small ($\lesssim 10$ pc) displacements between the AGN (locating the SBH) and the center of the galaxy (the mean photocenter) in 10 of them. Excluding objects that may be affected by large-scale isophotal asymmetries, we consider six galaxies to have detected displacements, including M87, where a displacement was previously reported by Batcheldor et al. 2010. In individual objects, these displacements can be attributed to residual gravitational recoil oscillations following a major or minor merger within the last few Gyr. For plausible merger rates, however, there is a high probability of larger displacements than those observed, if SBH coalescence took place in these galaxies. Remarkably, the AGN-photocenter displacements are approximately aligned with the radio source axis in four of the six galaxies with displacements, including three of the four having relatively powerful kpc-scale jets. This suggests intrinsic asymmetries in radio jet power as a possible displacement mechanism, although approximate alignments are also expected for gravitational recoil. Orbital motion in SBH binaries and interactions with massive perturbers can produce the observed displacement amplitudes but do not offer a ready explanation for the alignments.

It remains uncertain which continuum and emission line diagnostics best indicate the bolometric powers of active galactic nuclei (AGNs), especially given the attenuation caused by the circumnuclear material and the possible contamination by components related to star formation. Here we use mid-IR spectra along with multiwavelength data to investigate the merit of various diagnostics of AGN radiative power, including the mid-IR [Ne III] lambda25.89 mum and [O IV] lambda25.89 mum fine-structure lines, the optical [O III] lambda5007 forbidden line, and mid-IR 24 mum, 5 GHz radio, and X-ray continuum emission, for complete samples of 46 2Jy radio galaxies (0.05 < z < 0.7) and 17 3CRR FRII radio galaxies (z < 0.1). We find that the mid-IR [O IV] line is the most reliable indicator of AGN power for powerful radio-loud AGNs. By assuming that the [O IV] is emitted isotropically, and comparing the [O III] and 24 mum luminosities of the broad- and narrow-line AGNs in our samples at fixed [O IV] luminosity, we show that the [O III] and 24 mum emission are both mildly attenuated in the narrow-line compared to the broad-line objects by a factor of #2. However, despite this attenuation, the [O III] and 24 mum luminosities are better AGN power indicators for our sample than either the 5 GHz radio or the X-ray continuum luminosities. We also detect the mid-IR 9.7 mum silicate feature in the spectra of many objects but not ubiquitously: at least 40% of the sample shows no clear evidence for these features. We conclude that, for the majority of powerful radio galaxies, the mid-IR lines are powered by AGN photoionization.

We have directly measured the energy threshold and efficiency for bubble nucleation from iodine recoils in a CF3I bubble chamber in the energy range of interest for a dark matter search. These interactions cannot be probed by standard neutron calibration methods, so we develop a new technique by observing the elastic scattering of 12 GeV/c negative pions. The pions are tracked with a silicon pixel telescope and the reconstructed scattering angle provides a measure of the nuclear recoil kinetic energy. The bubble chamber was operated with a nominal threshold of (13.6+-0.6) keV. Interpretation of the results depends on the response to fluorine and carbon recoils, but in general we find agreement with the predictions of the classical bubble nucleation theory. This measurement confirms the applicability of CF3I as a target for spin-independent dark matter interactions and represents a novel technique for calibration of superheated fluid detectors.

Dark matter detectors require calibrations of their energy scale and efficiency to detect nuclear recoils in the $1-50$ keV range. Most calibrations use neutron scattering and require MCNP or Geant4 simulations of neutron propagation through the detector. For most nuclei heavier than $^{16}$O, these simulations' libraries ignore the contribution of resolved resonances to the neutron elastic differential cross-section. For many isotopes and neutron energies of importance to dark matter detection, this invalid assumption can severely distort simulated nuclear recoil spectra. The correct angular distributions can be calculated from the resonance parameters using R-matrix formalism. A set of neutron scattering libraries with high resolution angular distributions for MCNP and Geant4 of $^{19}$F, $^{40}$Ar, $^{50,52}$Cr, $^{56}$Fe, $^{136}$Xe, and $^{206,207,208}$Pb is presented. An MCNPX library for simulating the production of low-energy neutrons in the $^9$Be$(\gamma,n)^8$Be reaction is also presented. Example dark matter detector calibrations are simulated with the new libraries showing how detector sensitivity could be overestimated by factors of two by relying on existing MCNP and Geant4 libraries.

We present two-dimensional gaseous kinematics of the inner 1.1 x 1.6kpc^2 of the Seyfert 2 galaxy NGC2110, from optical spectra obtained with the GMOS integral field spectrograph on the Gemini South telescope at a spatial resolution of 100pc. Gas emission is observed over the whole field-of-view, with complex - and frequently double - emission-line profiles. We have identified four components in the emitting gas, according to their velocity dispersion (sigma), which we refer to as: (1) warm gas disk (sigma = 100-220km/s); (2) cold gas disk (sigma = 60-90km/s); (3) nuclear component (sigma = 220-600km/s); and (4) northern cloud (sigma = 60-80km/s). Both the cold and warm disk components are dominated by rotation and have similar gas densities, but the cold gas disk has lower velocity dispersions and reaches higher rotation velocities. We attribute the warm gas disk to a thick gas layer which encompasses the cold disk as observed in some edge-on spiral galaxies. After subtraction of a rotation model from the cold disk velocity field, we observe excess blueshifts of 50km/s in the far side of the galaxy as well as similar excess redshifts in the near side. These residuals can be interpreted as due to nuclear inflow in the cold gas, with an estimated ionized gas mass inflow rate of 2.2 x 10^(-2)Msun/yr. We have also subtracted a rotating model from the warm disk velocity field and found excess blueshifts of 100km/s to the SW of the nucleus and excess redshifts of 40km/s to the NE, which we attribute to gas disturbed by an interaction with a nuclear spherical outflow. This nuclear outflow is the origin of the nuclear component observed within the inner 300pc and it has a mass outflow rate of 0.9Msun/yr. In a region between 1" and 4" north of the nucleus we find a new low sigma component of ionized gas which we attribute to a high latitude cloud photoionized by the nuclear source.

We present a two-dimensional analysis of the gaseous excitation and kinematics of the inner 2.5 x 1.7 kpc^2 of the LINER/Seyfert 1 galaxy Arp 102B, from optical spectra obtained with the GMOS integral field spectrograph on the Gemini North telescope at a spatial resolution of 250 pc. Emission-line flux maps show the same two-armed nuclear spiral we have discovered in previous observations with the HST-ACS camera. One arm reaches 1 kpc to the east and the other 500 pc to the west, with a 8.4 GHz VLA bent radio jet correlating with the former. The gas density is highest (500 - 900 cm^(-3)) at the nucleus and in the northern border of the east arm, at a region where the radio jet seems to be deflected. Channel maps show blueshifts but also some redshifts at the eastern arm and jet location which can be interpreted as originated in the front and back walls of an outflow pushed by the radio jet, suggesting also that the outflow is launched close to the plane of the sky. We estimate a mass outflow rate along the east arm of 0.26 - 0.32 Msun yr^(-1) (depending on the assumed outflow geometry), which is between one and two orders of magnitude higher than the mass accretion rate to the active nucleus, implying that there is mass-loading of the nuclear outflow from circumnuclear gas. The power of this outflow is 0.06 - 0.3%Lbol. We propose a scenario in which gas has been recently captured by Arp 102B in an interaction with Arp 102A, settling in a disk rotating around the nucleus of Arp 102B and triggering its nuclear activity. A nuclear jet is pushing the circumnuclear gas, giving origin to the nuclear arms. A blueshifted emitting gas knot is observed at 300 pc south-east from the nucleus and can be interpreted as another (more compact) outflow, with a possible counterpart to the north-west.

From an optical spectroscopic survey of 3CR radiogalaxies (RGs) with z<0.3, we discovered three objects characterized by an extremely low level of gas excitation and a large deficit of line emission with respect to RGs of similar radio luminosity. We interpreted these objects as relic active galactic nuclei (AGN), i.e., sources observed after a large drop in their nuclear activity. We here present new spectroscopic observations for these three galaxies and for a group of "candidate" relics. None of the candidates can be convincingly confirmed. From the new data for the three relics, we estimate the density of the line-emitting gas. This enables us to explore the temporal evolution of the line ratios after the AGN "death". The characteristic timescale is the light-crossing time of the emission line region, a few thousand years, too short to correspond to a substantial population of relic RGs. Additional mechanisms of gas ionization, such as "relic shocks" from their past high power phase or stellar sources, should also be considered to account for the spectroscopic properties of the relic RGs. Relic RGs appear to be a mixed bag of sources in different phases of evolution, including AGN recently (~10,000 years ago) quenched, galaxies that have been inactive for at least one million years, and objects caught during the transition from a powerful RG to a low power FRI source.

We have devised a technique for measuring some of the most important nuclear reactions in stars which we expect to provide considerable improvement over previous experiments. Adapting ideas from dark matter search experiments with bubble chambers, we have found that a superheated liquid is sensitive to recoils produced from \gamma-rays photodisintegrating the nuclei of the liquid. The main advantage of the new target-detector system is a gain in yield of six orders of magnitude over conventional gas targets due to the higher mass density of liquids. Also, the detector is practically insensitive to the \gamma-ray beam itself, thus allowing it to detect only the products of the nuclear reaction of interest. The first set of tests of a superheated target with a narrow bandwidth \gamma-ray beam was completed and the results demonstrate the feasibility of the scheme. The new data are successfully described by an R-matrix model using published resonance parameters. With the increase in luminosity of the next generation \gamma-ray beam facilities, the measurement of thermonuclear rates in the stellar Gamow window would become possible.

We present new CRIRES spectroscopic observations of BrGamma in the nuclear region of the Circinus galaxy, obtained with the aim of measuring the black hole (BH) mass with the spectroastrometric technique. The Circinus galaxy is an ideal benchmark for the spectroastrometric technique given its proximity and secure BH measurement obtained with the observation of its nuclear H2O maser disk. The kinematical data have been analyzed both with the classical method based on the analysis of the rotation curves and with the new method developed by us and based on spectroastrometry. The classical method indicates that the gas disk rotates in the gravitational potential of an extended stellar mass distribution and a spatially unresolved mass of (1.7 +- 0.2) 10^7 Msun, concentrated within r < 7 pc. The new method is capable of probing gas rotation at scales which are a factor ~3.5 smaller than those probed by the rotation curve analysis. The dynamical mass spatially unresolved with the spectroastrometric method is a factor ~2 smaller, 7.9 (+1.4 -1.1) 10^6 Msun indicating that spectroastrometry has been able to spatially resolve the nuclear mass distribution down to 2 pc scales. This unresolved mass is still a factor ~4.5 larger than the BH mass measurement obtained with the H2O maser emission indicating that it has not been possible to resolve the sphere of influence of the BH. Based on literature data, this spatially unresolved dynamical mass distribution is likely dominated by molecular gas and it has been tentatively identified with the circum-nuclear torus which prevents a direct view of the central BH in Circinus. This mass distribution, with a size of ~2pc, is similar in shape to that of the star cluster of the Milky Way suggesting that a molecular torus, forming stars at a high rate, might be the earlier evolutionary stage of the nuclear star clusters which are common in late type spirals.

New data are reported from the operation of a 4.0 kg CF$_{3}$I bubble chamber in the 6800-foot-deep SNOLAB underground laboratory. The effectiveness of ultrasound analysis in discriminating alpha-decay background events from single nuclear recoils has been confirmed, with a lower bound of $>$99.3% rejection of alpha-decay events. Twenty single nuclear recoil event candidates and three multiple bubble events were observed during a total exposure of 553 kg-days distributed over three different bubble nucleation thresholds. The effective exposure for single bubble recoil-like events was 437.4 kg-days. A neutron background internal to the apparatus, of known origin, is estimated to account for five single nuclear recoil events and is consistent with the observed rate of multiple bubble events. This observation provides world best direct detection constraints on WIMP-proton spin-dependent scattering for WIMP masses $>$20 GeV/c$^{2}$ and demonstrates significant sensitivity for spin-independent interactions.

(Abridged) We present the results from new 15 ks Chandra-ACIS and 4.9 GHz Very Large Array observations of 13 galaxies hosting low luminosity AGN. This completes the multiwavelength study of a sample of 51 nearby early-type galaxies described in Capetti & Balmaverde (2005, 2006); Balmaverde & Capetti (2006). The aim of the three previous papers was to explore the connection between the host galaxies and AGN activity in a radio-selected sample. We detect nuclear X-ray emission in eight sources and radio emission in all but one (viz., UGC6985). The new VLA observations improve the spatial resolution by a factor of ten: the presence of nuclear radio sources in 12 of the 13 galaxies confirms their AGN nature. As previously indicated, the behavior of the X-ray and radio emission in these sources depends strongly on the form of their optical surface brightness profiles derived from Hubble Space Telescope imaging, i.e., on their classification as "core", "power-law" or "intermediate" galaxies. With more than twice the number of "power-law" and "intermediate" galaxies compared to previous work, we confirm with a much higher statistical significance that these galaxies lie well above the radio-X-ray correlation established in FRI radio galaxies and the low-luminosity "core" galaxies. This result highlights the fact that the "radio-loud/radio-quiet" dichotomy is a function of the host galaxy's optical surface brightness profile. We present radio-optical-X-ray spectral indices for all 51 sample galaxies. Survival statistics point to significant differences in the radio-to-optical and radio-to-X-ray spectral indices between the "core" and "power-law" galaxies (Gehan's Generalized Wilcoxon test probability "p" for the two classes being statistically similar is <10^-5), but not in the optical-to-X-ray spectral indices (p=0.25).

We present new Spitzer IRS spectroscopy of Cygnus A, one of the most luminous radio sources in the local universe. Data on the inner 20" are combined with new reductions of MIPS and IRAC photometry as well as data from the literature to form a radio through mid-infrared spectral energy distribution (SED). This SED is then modeled as a combination of torus reprocessed active galactic nucleus (AGN) radiation, dust enshrouded starburst, and a synchrotron jet. This combination of physically motivated components successfully reproduces the observed emission over almost 5 dex in frequency. The bolometric AGN luminosity is found to be 10^12 L_\odot (90% of LIR), with a clumpy AGN-heated dust medium extending to \sim130 pc from the supermassive black hole. Evidence is seen for a break or cutoff in the core synchrotron emission. The associated population of relativistic electrons could in principle be responsible for some of the observed X-ray emission though the synchrotron self-Compton mechanism. The SED requires a cool dust component, consistent with dust-reprocessed radiation from ongoing star formation. Star formation contributes at least 6 \times 10^10 L_\odot to the bolometric output of Cygnus A, corresponding to a star formation rate of \sim10 M_\odot yr-1.

We present deep Spitzer/IRS spectra for complete samples of 46 2Jy radio galaxies (0.05<z<0.7) and 19 3CRR FRII radio galaxies (z<0.1), and use the detection of polycyclic aromatic hydrocarbon (PAH) features to examine the incidence of contemporaneous star formation and radio-loud AGN activity. Our analysis reveals PAH features in only a minority (30%) of the objects with good IRS spectra. Using the wealth of complementary data available for the 2Jy and 3CRR samples we make detailed comparisons between a range of star formation diagnostics: optical continuum spectroscopy, mid- to far-IR (MFIR) color, far-IR excess and PAH detection. There is good agreement between the various diagnostic techniques: most candidates identified to have star formation activity on the basis of PAH detection are also identified using at least two of the other techniques. We find that only 35% of the combined 2Jy and 3CRR sample show evidence for recent star formation activity (RSFA) at optical and/or MFIR wavelengths. This result argues strongly against the idea of a close link between starburst and powerful radio-loud AGN activity, reinforcing the view that, although a large fraction of powerful radio galaxies may be triggered in galaxy interactions, only a minority are triggered at the peaks of star formation activity in major, gas-rich mergers. However, we find that compact radio sources (D < 15 kpc) show a significantly higher incidence of RSFA (>75%) than their more extended counterparts (=15 -- 25%). We discuss this result in the context of a possible bias towards the selection of compact radio sources triggered in gas-rich environments.

We measure the black hole mass in the nearby active galaxy Centaurus A (NGC 5128) using a new method based on spectroastrometry of a rotating gas disk. The spectroastrometric approach consists in measuring the photocenter position of emission lines for different velocity channels. In a previous paper we focused on the basic methodology and the advantages of the spectroastrometric approach with a detailed set of simulations demonstrating the possibilities for black hole mass measurements going below the conventional spatial resolution. In this paper we apply the spectroastrometric method to multiple longslit and integral field near infrared spectroscopic observations of Centaurus A. We find that the application of the spectroastrometric method provides results perfectly consistent with the more complex classical method based on rotation curves: the measured BH mass is nearly independent of the observational setup and spatial resolution and the spectroastrometric method allows the gas dynamics to be probed down to spatial scales of ~0.02", i.e. 1/10 of the spatial resolution and ~1/50 of BH sphere of influence radius. The best estimate for the BH mass based on kinematics of the ionized gas is then log(MBH (sin i)^2/M\odot)=7.5 \pm 0.1 which corresponds to MBH = 9.6(+2.5-1.8) \times 10^7 M\odot for an assumed disk inclination of i = 35deg. The complementarity of this method with the classic rotation curve method will allow us to put constraints on the disk inclination which cannot be otherwise derived from spectroastrometry. With the application to Centaurus A, we have shown that spectroastrometry opens up the possibility of probing spatial scales smaller than the spatial resolution, extending the measured MBH range to new domains which are currently not accessible: smaller BHs in the local universe and similar BHs in more distant galaxies.

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 report the discovery of a two-armed mini-spiral structure within the inner kiloparsec of the E0 LINER/Seyfert 1 galaxy Arp102B. The arms are observed in H-alpha emission and located East and West of the nucleus, extending up to about 1 kpc from it. We use narrow-band imaging from the Hubble Space Telescope Advanced Camera for Surveys, in combination with archival VLA radio images at 3.6 and 6 cm to investigate the origin of the nuclear spiral. From the H-alpha luminosity of the spiral, we obtain an ionized gas mass of the order of one million solar masses. One possibility is that the nuclear spiral represents a gas inflow triggered by a recent accretion event which has replenished the accretion disk, giving rise to the double-peaked emission-line profiles characteristic of Arp102B. However, the radio images show a one-sided curved jet which correlates with the eastern spiral arm observed in the H-alpha image. A published milliarcsecond radio image also shows one-sided structure at position angle about 40 degrees, approximately aligned with the inner part of the eastern spiral arm. The absence of a radio counter-part to the western spiral arm is tentatively interpreted as indicating that the jet is relativistic, with an estimated speed of 0.45c. Estimates of the jet kinetic energy and the ionizing luminosity of the active nucleus indicate that both are capable of ionizing the gas along the spiral arms. We conclude that, although the gas in the nuclear region may have originated in an accretion event, the mini-spiral is most likely the result of a jet-cloud interaction rather than an inflowing stream.

From an optical spectroscopic survey of 3CR radio galaxies with z<0.3, we discovered a new spectroscopic class of powerful radio-loud AGN. The defining characteristics of these galaxies are that compared with radio galaxies of similar radio luminosity they have: a [O III]\Hb ratio of ~0.5, indicative of an extremely low level of gas excitation; a large deficit of [O III] emission and radio core power. We interpret these objects as relic AGN, i.e. sources that experienced a large drop in their level of nuclear activity, causing a decrease in their nuclear and line luminosity. This class opens a novel approach to investigating lifetimes and duty cycles of AGN.

We present stellar and gaseous kinematics of the inner 120x250pc^2 of the Liner/Seyfert 1 galaxy M81, from optical spectra obtained with the GMOS integral field spectrograph on the Gemini North telescope at a spatial resolution of 10pc. The stellar velocity field shows circular rotation but deviations are observed close to the minor axis which can be attributed to stellar motions possibly associated to a nuclear bar. The stellar velocity dispersion of the bulge is 162km/s leading to a black hole mass of M_BH=5.5x10^7M_sun based on the M_BH-sigma relationship. The gas kinematics is dominated by non-circular motions and the subtraction of the stellar velocity field reveals blueshifts of ~-100km/s on the far side of the galaxy and a few redshifts on the near side. These characteristics can be interpreted in terms of streaming towards the center if the gas is in the plane. On the basis of the observed velocities and geometry of the flow, we estimate a mass inflow rate in ionized gas of ~4.0x10^-3M_sun/year, which is of the order of the accretion rate necessary to power the LINER nucleus of M81. We have also applied the technique of Principal Component Analysis (PCA) to our data, which reveals the presence of a rotating nuclear gas disk within ~50pc from the nucleus and a compact outflow, approximately perpendicular to the disk. The PCA combined with the observed gas velocity field shows that the nuclear disk is being fed by gas circulating in the galaxy plane. The presence of the outflow is supported by a compact jet seen in radio observations at a similar orientation, as well as by an enhancement of the [OI]\Halpha line ratio, probably resulting from shock excitation of the circumnuclear gas by the radio jet. With these observations we are thus resolving both the feeding -- via the nuclear disk and observed gas inflow, and the feedback -- via the outflow, around the nucleus of M81.

We have observed the composite active galactic nucleus (AGN)-starburst galaxy NGC 6764 with the Very Long Baseline Array at 1.6 and 4.9 GHz. These observations have detected a "core-jet" structure and a possible weak counterjet component at 1.6 GHz. The upper limits to the core and jet (1.6-4.9 GHz) spectral index are 0.6 and 0.3, respectively. Taken together with the high brightness temperature of ~10^7 K for the core region, the radio emission appears to be coming from a synchrotron jet. At a position angle of 25 degrees, the parsec-scale jet seems to be pointing closely toward the western edge of the southern kpc-scale bubble in NGC 6764. A real connection between the parsec- and sub-kpc-scale emission would not only suggest the presence of a curved jet, but also a close link between the AGN jet and the radio bubbles in NGC 6764. We demonstrate that a precessing jet model can explain the radio morphology from parsec- to sub-kpc scales, and the model best-fit parameters of jet speed and orientation are fully consistent with the observed jet-to-counterjet surface brightness ratio. The jet however appears to be disrupted on scales of hundreds of parsecs, possibly due to interaction with and entrainment of the interstellar medium gas, which subsequently leads to the formation of bubbles. The jet energetics in NGC 6764 suggest that it would take 12-21 Myr to inflate the (southern) bubble. This timescale corresponds roughly to the starburst episode that took place in NGC 6764 about 15-50 Myr ago, and could be indicative of a close connection between jet formation and the starburst activity in this galaxy.

We present Spitzer photometric data for a complete sample of 19 low redshift (z<0.1) 3CRR radio galaxies as part of our efforts to understand the origin of the prodigious mid- to far-infrared (MFIR) emission from radio-loud AGN. Our results show a correlation between AGN power (indicated by [OIII] 5007 emission line luminosity) and 24 micron luminosity. This result is consistent with the 24 micron thermal emission originating from warm dust heated directly by AGN illumination. Applying the same correlation test for 70 micron luminosity against [OIII] luminosity we find this relation to suffer from increased scatter compared to that of 24 micron. In line with our results for the higher-radio-frequency-selected 2Jy sample, we are able to show that much of this increased scatter is due to heating by starbursts which boost the far-infrared emission at 70 micron in a minority of objects (17-35%). Overall this study supports previous work indicating AGN illumination as the dominant heating mechanism for MFIR emitting dust in the majority of low to intermediate redshift radio galaxies (0.03<z<0.7), with the advantage of strong statistical evidence. However, we find evidence that the low redshift broad-line objects (z<0.1) are distinct in terms of their positions on the MFIR vs. [OIII] correlations.

Isophotal analysis of M87, using data from the Advanced Camera for Surveys, reveals a projected displacement of 6.8 +/- 0.8 pc (~ 0.1 arcsec) between the nuclear point source (presumed to be the location of the supermassive black hole, SMBH) and the photo-center of the galaxy. The displacement is along a position angle of 307 +/- 17 degrees and is consistent with the jet axis. This suggests the active SMBH in M87 does not currently reside at the galaxy center of mass, but is displaced in the counter-jet direction. Possible explanations for the displacement include orbital motion of an SMBH binary, gravitational perturbations due to massive objects (e.g., globular clusters), acceleration by an asymmetric or intrinsically one-sided jet, and gravitational recoil resulting from the coalescence of an SMBH binary. The displacement direction favors the latter two mechanisms. However, jet asymmetry is only viable, at the observed accretion rate, for a jet age of >0.1 Gyr and if the galaxy restoring force is negligible. This could be the case in the low density core of M87. A moderate recoil ~1 Myr ago might explain the disturbed nature of the nuclear gas disk, could be aligned with the jet axis, and can produce the observed offset. Alternatively, the displacement could be due to residual oscillations resulting from a large recoil that occurred in the aftermath of a major merger any time in the last 1 Gyr.

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.

The mid-far-infrared spectral energy distributions (SEDs) of 83 active galaxies, mostly Seyfert galaxies, selected from the extended 12 micron sample are presented. The data were collected using all three instruments, IRAC, IRS, and MIPS, aboard the Spitzer Space Telescope. The IRS data were obtained in spectral mapping mode, and the photometric data from IRAC and IRS were extracted from matched, 20 arcsec diameter circular apertures. The MIPS data were obtained in SED mode, providing very low resolution spectroscopy (R ~ 20) between ~ 55 and 90 microns in a larger, 20 by 30 arcsec synthetic aperture. We further present the data from a spectral decomposition of the SEDs, including equivalent widths and fluxes of key emission lines; silicate 10 and 18 micron emission and absorption strengths; IRAC magnitudes; and mid-far infrared spectral indices. Finally, we examine the SEDs averaged within optical classifications of activity. We find that the infrared SEDs of Seyfert 1s and Seyfert 2s with hidden broad line regions (HBLR, as revealed by spectropolarimetry or other technique) are qualitatively similar, except that Seyfert 1s show silicate emission and HBLR Seyfert 2s show silicate absorption. The infrared SEDs of other classes with the 12 micron sample, including Seyfert 1.8-1.9, non-HBLR Seyfert 2 (not yet shown to hide a type 1 nucleus), LINER and HII galaxies, appear to be dominated by star-formation, as evidenced by blue IRAC colors, strong PAH emission, and strong far-infrared continuum emission, measured relative to mid-infrared continuum emission.

This is the first in a series of papers in which we study the application of spectroastrometry in the context of gas kinematical studies aimed at measuring the mass of supermassive black holes. The spectroastrometrical method consists in measuring the photocenter of light emission in different wavelength or velocity channels. In particular we explore the potential of spectroastrometry of gas emission lines in galaxy nuclei to constrain the kinematics of rotating gas disks and to measure the mass of putative supermassive black holes. By means of detailed simulations and test cases, we show that the fundamental advantage of spectroastrometry is that it can provide information on the gravitational potential of a galaxy on scales significantly smaller (~ 1/10) than the limit imposed by the spatial resolution of the observations. We then describe a simple method to infer detailed kinematical informations from spectroastrometry in longslit spectra and to measure the mass of nuclear mass concentrations. Such method can be applied straightforwardly to integral field spectra, which do not have the complexities due to a partial spatial covering of the source in the case of longslit spectra.

We present an analysis of Spitzer IRS spectroscopy of 83 active galaxies from the extended 12 micron sample. We find rank correlations between several tracers of star formation which suggest that (1) the PAH feature is a reliable tracer of star formation, (2) there is a significant contribution to the heating of the cool dust by stars, (3) the H$_2$ emission is also primarily excited by star formation. The 55-90 vs. 20-30 spectral index plot is also a diagnostic of the relative contribution of Starburst to AGN. We see there is a large change in spectral index across the sample. Thus, the contribution to the IR spectrum from the AGN and starburst components can be comparable in magnitude but the relative contribution also varies widely across the sample. We find rank correlations between several AGN tracers. We search for correlations between AGN and Starburst tracers and we conclude that the AGN and Starburst tracers are not correlated. This is consistent with our conclusion that the relative strength of the AGN and Starburst components varies widely across the sample. Thus, there is no simple link between AGN fueling and Black Hole Growth and star formation in these galaxies. The distribution of Sil 10 micron and 18 micron strengths is consistent with the clumpy torus models of Sirocky et al. We find a rank correlation between the [NeV] 14 micron line and the 6.7 micron continuum which may be due to an extended component of hot dust. The Sy 2s with a Hidden Broad Line Region (HBLR) have a higher ratio of AGN to Starburst contribution to the SED than Sy 2s without an HBLR. This may contribute to the detection of the HBLR in polarized light. The Sy 2s with an HBLR are more similar to the Sy 1s than they are to the Sy 2s without an HBLR.

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.

The application of the virial theorem to the Broad Line Region of Active Galactic Nuclei allows Black Hole mass estimates for large samples of objects at all redshifts. In a recent paper we showed that ionizing radiation pressure onto BLR clouds affects virial BH mass estimates and we provided empirically calibrated corrections. More recently, a new test of the importance of radiation forces has been proposed: the MBH-sigma relation has been used to estimate MBH for a sample of type-2 AGN and virial relations (with and without radiation pressure) for a sample of type-1 AGN extracted from the same parent population. The observed L/LEdd distribution based on virial BH masses is in good agreement with that based on MBH-sigma only if radiation pressure effects are negligible, otherwise significant discrepancies are observed. In this paper we investigate the effects of intrinsic dispersions associated to the virial relations providing MBH, and we show that they explain the discrepancies between the observed L/LEdd distributions of type-1 and type-2 AGN. These discrepancies in the L/LEdd distributions are present regardless of the general importance of radiation forces, which must be negligible only for a small fraction of sources with large L/LEdd. Average radiation pressure corrections should then be applied in virial MBH estimators until their dependence on observed source physical properties has been fully calibrated. Finally, the comparison between MBH and L/LEdd distributions derived from sigma-based and virial estimators can constrain the variance of BLR physical properties in AGN.

The findings of a nine orbit calibration plan carried out during HST Cycle 15, to fully determine the NICMOS camera 2 (2.0 micron) polarization calibration to high accuracy, are reported. Recently Ueta et al. and Batcheldor et al. have suggested that NICMOS possesses a residual instrumental polarization at a level of 1.2-1.5%. This would completely inhibit the data reduction in a number of GO programs, and hamper the ability of the instrument to perform high accuracy polarimetry. We obtained polarimetric calibration observations of three polarimetric standards at three spacecraft roll angles separated by ~60deg. Combined with archival data, these observations were used to characterize the residual instrumental polarization in order for NICMOS to reach its full potential of accurate imaging polarimetry at p~1%. Using these data, we place an 0.6% upper limit on the instrumental polarization and calculate values of the parallel transmission coefficients that reproduce the ground-based results for the polarimetric standards. The uncertainties associated with the parallel transmission coefficients, a result of the photometric repeatability of the observations, are seen to dominate the accuracy of p and theta. However, the updated coefficients do allow imaging polarimetry of targets with p~1.0% at an accuracy of +/-0.6% and +/-15deg. This work enables a new caliber of science with HST.

The ability of NICMOS to perform high accuracy polarimetry is currently hampered by an uncalibrated residual instrumental polarization at a level of 1.2-1.5%. To better quantify and characterize this residual we obtained observations of three polarimetric standard stars at three separate space-craft roll angles. Combined with archival data, these observations were used to characterize the residual instrumental polarization to enable NICMOS to reach its full polarimetric potential. Using these data, we calculate values of the parallel transmission coefficients that reproduce the ground-based results for the polarimetric standards. The uncertainties associated with the parallel transmission coefficients, a result of the photometric repeatability of the observations, dominate the accuracy of p and theta. However, the new coefficients now enable imaging polarimetry of targets with p~1.0% at an accuracy of +/-0.6% and +/-15 degrees.

The application of the virial theorem provides a tool to estimate supermassive black hole (BH) masses in large samples of active galactic nuclei (AGN) with broad emission lines at all redshifts and luminosities, if the broad line region (BLR) is gravitationally bound. In this paper we discuss the importance of radiation forces on BLR clouds arising from the deposition of momentum by ionizing photons. Such radiation forces counteract gravitational ones and, if not taken into account, BH masses can be severely underestimated. We provide virial relations corrected for the effect of radiation pressure and we discuss their physical meaning and application. If these corrections to virial masses, calibrated with low luminosity objects, are extrapolated to high luminosities then the BLRs of most quasars might be gravitationally unbound. The importance of radiation forces in high luminosity objects must be thoroughly investigated to assess the reliability of quasar BH masses.

The infrared is a key wavelength regime for probing the dusty, obscured nuclear regions of active galaxies. We present results from an infrared study of 87 nearby Seyfert galaxies using the Spitzer Space Telescope and ground-based telescopes. Combining detailed modelling of the 3 - 100 micron spectral energy distributions with mid-IR spectral diagnostics and near-infrared observations, we find broad support for the unified model of AGNs. The IR emission of Seyfert 1s and 2s is consistent with their having the same type of central engine viewed at a different orientation. The nature of the putative torus is becoming clearer; in particular we present evidence that it is likely a clumpy medium. Mid-infrared correlations between tracers of star formation and AGN ionizing luminosity reveal the starburst-AGN connection implied by the black hole/bulge mass relation, however it is not yet clear if this is due to feedback.

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 consider the effect of radiation pressure from ionizing photons on black hole (BH) mass estimates based on the application of the virial theorem to broad emission lines in AGN spectra. BH masses based only on the virial product V^2R and neglecting the effect of radiation pressure can be severely underestimated especially in objects close to the Eddington limit. We provide an empirical calibration of the correction for radiation pressure and we show that it is consistent with a simple physical model in which BLR clouds are optically thick to ionizing radiation and have average column densities of NH~10^23 cm^-2. This value is remarkably similar to what is required in standard BLR photoionization models to explain observed spectra. With the inclusion of radiation pressure the discrepancy between virial BH masses based on single epoch spectra and on reverberation mapping data drops from 0.4 to 0.2 dex rms. The use of single epoch observations as surrogates of reverberation mapping campaigns can thus provide more accurate BH masses than previously thought. Finally, we show that Narrow Line Seyfert 1 (NLS1) galaxies have apparently low BH masses because they are radiating close to their Eddington limit. After the radiation pressure correction, NLS1 galaxies have BH masses similar to other broad line AGNs and follow the same MBH-sigma/L relations as other active and normal galaxies. Radiation forces arising from ionizing photon momentum deposition constitute an important physical effect which must be taken into account when computing virial BH masses.

We report the discovery of gas streaming motions along nuclear spiral arms towards the LINER nucleus of the galaxy NGC 6951. The observations, obtained using the GMOS integral field spectrograph on the Gemini North telescope, yielded maps of the flux distributions and gas kinematics in the Halpha, [NII]6584 and [SII]6717,31 emission lines of the inner 7x5 arcsec^2 of the galaxy. This region includes a circumnuclear star-forming ring with radius 500pc, a nuclear spiral inside the ring and the LINER nucleus. The kinematics of the ionized gas is dominated by rotation, but subtraction of a kinematic model of a rotating exponential disk reveals deviations from circular rotation within the nuclear ring which can be attributed to (1) streaming motions along the nuclear spiral arms and (2) a bipolar outflow which seems to be associated to a nuclear jet. On the basis of the observed streaming velocities and geometry of the spiral arms we estimate a mass inflow rate of ionized gas of 3x10^(-4) Msun/yr, which is of the order of the accretion rate necessary to power the LINER nucleus of NGC 6951. Similar streaming motions towards the nucleus of another galaxy with LINER nucleus -- NGC 1097 -- have been reported by our group in a previous paper. Taken together, these results support a scenario in which nuclear spirals are channels through which matter is transferred from galactic scales to the nuclear region to feed the supermassive black hole.

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 are conducting a large observing program with the Spitzer Space Telescope to determine the mid-to-far infrared spectral energy distributions of a well-defined sample of 87 nearby, 12 micron-selected Seyfert galaxies. In this paper we present the results of IRS low-resolution spectroscopy of a statistically representative subsample of 51 of the galaxies (59%), with an analysis of the continuum shapes and a comparison of the Seyfert types. We find that the spectra clearly divide into groups based on their continuum shapes and spectral features. Some spectral features are clearly related to a starburst contribution to the IR spectrum, while the mechanisms producing observed power-law continuum shapes, attributed to an AGN component, may be dust or non-thermal emission. The infrared spectral types appear to be related to the Seyfert types. Principal component analysis results suggest that the relative contribution of starburst emission may be the dominant cause of variance in the observed spectra. We find that Sy 1's have higher ratios of IR/radio emission than Sy 2's, as predicted by the unified model if the torus is optically thick in the mid-IR. However, smooth-density torus models predict a much larger difference between type 1's and 2's than observed in our sample. Our observations may be consistent with clumpy torus models containing a steep radial distribution of optically thick dense clumps. (Abridged)

We present the results of IRS low-resolution spectroscopy of 51 Seyfert galaxies, part of a large Spitzer observing program to determine the mid-to-far infrared spectral energy distributions of a well-defined sample of 87 nearby, 12 micron-selected Seyferts. We find that the spectra clearly divide into groups based on their continuum shapes and spectral features. The infrared spectral types appear to be related to the Seyfert types. Some features are clearly related to a starburst contribution to the IR spectrum, while the observed power-law continuum shapes, attributed to the AGN, may be dust or non-thermal emission. Principal component analysis results suggest that the relative contribution of starburst emission is the dominant cause of variance in the spectra. We find that the Sy 2's show on average stronger starburst contributions than the Sy 1's.

We have used GMOS-IFU and high resolution HST-ACS observations to map, in unprecedented detail, the gas velocity field and structure within the 0.7 kpc circumnuclear ring of the SBb LINER/Seyfert 1 galaxy NGC 1097. We find clear evidence of radial streaming motions associated with spiral structures leading to the unresolved (<3.5 parsecs) nucleus, which we interpret as part of the fueling chain by which gas is transported to the nuclear starburst and supermassive black hole.

The value of accurately knowing the absolute calibration of the polarizing elements in the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) becomes especially important when conducting studies which require measuring degrees of polarization of close to 1% in the near infrared. We present a comprehensive study of all previously observed polarimetric standards using the NIC2 camera on NICMOS. Considering both pre- and post-NICMOS Cooling System observations we find variations in the polarimetry consistent with the effects of sub-pixel mis-alignments and the point spread function. We also measure non-zero results from unpolarized standards indicating an instrumental polarization of p ~ 1.2%, theta ~ 88degrees. The lack of polarized and unpolarized standard stars with which to perform a comprehensive calibration study means we cannot be confident that the current calibration will be effective for a number of recent large NICMOS GO programs. Further observations of polarimetric standards are needed in order to fully characterize the behavior of NICMOS at around p=1%.

Magain et al. (2005) argued that the host galaxy of the quasar in HE0450-2958 is substantially under-luminous given the likely mass of its nuclear black hole. Using kinematical information from the spectra of the quasar and the companion galaxy, an ultra-luminous infrared galaxy, we test the hypothesis that the black hole powering the quasar was ejected from the companion galaxy during a merger. We find that the ejection model can be securely ruled out, since the kick velocity required to remove the black hole from the galaxy is greater than about 500 km/s, inconsistent with the presence of narrow emission line gas at the same redshift as the quasar nucleus. We also show that the quasar in HE0450-2958 has the spectral characteristics of a narrow-line Seyfert 1 galaxy and calculate a mass for its black hole that is roughly an order of magnitude smaller than estimated by Magain et al. The predicted luminosity of the host galaxy is then consistent with the upper limits inferred by those authors.

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 subarcsecond resolution mid infrared images of NGC 4151 at 10.8 micron and 18.2 micron. These images were taken with the University of Florida mid-IR camera/spectrometer OSCIR at the Gemini North 8-m telescope. We resolve emission at both 10.8 micron and 18.2 micron extending ~ 3.5" across at a P.A. of ~ 60 degrees. This coincides with the the narrow line region of NGC 4151 as observed in [OIII] by the Hubble Space Telescope. The most likely explanation for this extended mid-IR emission is dust in the narrow line region heated by a central engine. We find no extended emission associated with the proposed torus and place an upper limit on its mid-IR size of less than or equal to ~ 35 pc.

We analyze published reverberation mapping data for three Seyfert galaxies (NGC 3227, NGC 3516, and NGC 4593) to refine the mass estimate for the supermassive black hole in the center of each object. Treatment of the data in a manner more consistent with other large compilations of such masses allows us to more securely compare our results to wider samples of data, e.g., in the investigation of the M_bh-sigma relationship for active and quiescent galaxies.

We describe the results of a detailed study of the polarization properties of the broad H-alpha emission line in Type 1 Seyfert nuclei. Our analysis of these data points to a model in which the broad Balmer lines are emitted by a rotating disk, and are scattered in two main regions - one co-planar with the disk and within the circum-nuclear torus, the other, the polar scattering region, outside the torus but aligned with its axis. The relative importance of the two sources of polarized light is largely determined by the inclination of the system axis to the line-of-sight.

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 present the results of HST narrow band imaging of eleven Compact Steep Spectrum (CSS) radio sources. Five of them (3C 48, 3C 147, 3C303.1, 3C 277.1 and 4C 12.50) were part of a dedicated ``pointed'' program of deep line imaging ([OIII]). For six additional sources (3C 49, 3C 93.1, 3C 138, 3C 268.3, 3C305.1 and 3C343.1) ``snapshot'' images ([OIII] or [OII]) were taken from the HST archive. In all but one of the targets (3C 49) line emission has been detected and only in one case (3C 138) is unresolved. In four out five of the sources with deep observations, the line emission extends well beyond the size of the radio source but along the radio axis. Structures of similar surface brightness would have not been seen in the snapshot images. These emission line structures extend to scales of 10 to 30 kpc and cover a projected angle, when seen from the nucleus, of 30 to 110 degrees, indicating that the nuclear illumination is anisotropic. Photon counting arguments also support this interpretation. In six objects the radio emission extends over more than 1 arcsec. In these cases the line emission has an elongated structure, linking the nucleus to the radio-lobes, possibly tracing the path of the invisible radio jets. Nevertheless the emission line morphologies do not show the bow shocks at the extremities of the radio lobes expected if they are sources whose expansion is frustrated by a dense external medium. Our data favour the alternative model in which CSSs are the young phase of the large size radio sources. When ``pointed'' pure continuum images are available, there appears to be no alignment between radio and continuum emission which contradicts previous suggestions based on broad-band HST imaging.

The photoionization models for the narrow emission line regions of powerful radio galaxies have yet to be tested in depth. To this end, we present high-quality long-slit spectroscopy of the powerful double-nucleus radio galaxy 3C 321. The data have good enough spatial resolution to be able to trace the variation in emission-line properties on kpc scales. Continuum modelling and subtraction enables the faint emission line fluxes to be measured in several regions across the emission line nebula. We plot diagnostic line-ratio diagrams and compare them with the predictions of various photoionization models, finding that the data is best fit by models which assume a mixture of optically thin and thick clouds illuminated by a power-law continuum. The emission line kinematics, line ratios and deduced physical conditions show remarkably little variation across the source. We deduce a mean electron density of 400 +/- 120 cm-3 and a mean temperature of 11500 +/- 1500 K. Assuming a single population of optically thick line-emitting clouds, we calculate a mean photoionization parameter of (1.1 +/- 0.5) x 10e-2 and hence a photoionizing photon luminosity of Q ~ 10e55 -- 10e56 photon/s/sr. This indicates a central engine as luminous as that of the powerful quasar 3C 273, yet there is no evidence for such an energetically prolific central engine at either far-infrared or radio wavelengths. We therefore conclude that the mixed-media models, which give Q ~ 5 x 10e53 -- 5 x 10e54, represent a more likely scenario. As a by-product of the continuum subtraction we infer that young stellar populations account for ~ 0.4% of the visible stellar mass in the galaxy, and that these populations are spatially extended.

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 present sensitive, high angular resolution (6" x 5") 21-cm observations of the neutral hydrogen in the nearby barred Seyfert galaxy, NGC4151. These HI observations, obtained using the VLA in B-configuration, are the highest resolution to date of this galaxy, and reveal hitherto unprecedented detail in the distribution and kinematics of the HI on sub-kiloparsec scales. A complete analysis and discussion of the HI data are presented and the global properties of the galaxy are related to the bar dynamics presented in Paper I.

The nature of the optical--radio correlations for powerful radio galaxies is investigated using spectroscopic observations of a complete sample of southern 2Jy radio sources. In line with previous work, we find that significant correlations exist between the luminosities of the [OIII]5007, [OII]3727 and Hbeta emission lines and the radio luminosity. However, our observations are not easily reconciled with the idea that these correlations are caused by the increase in the power of the photoionizing quasar as the jet power increases, with average ISM properties not changing appreciably with redshift or radio power: not only do we find that the scatter in the L_[OIII] vs. L_radio correlation is significantly larger than in L_[OII]} vs. L_radio and L_Hbeta vs. L_radio correlations, but the ionization state deduced from the emission lines does not increase with radio power as predicted by the simple, constant ISM, photionization model. We conclude that: (a) there exists a considerable range in the quasar ionizing luminosity at a given redshift; and (b) that the mean density of the emission line clouds is larger in the high redshift/high power radio sources. The latter density enhancement may either be a consequence of the increased importance of jet-cloud interactions or, alternatively, due to a higher pressure in the confining hot ISM, in the high redshift objects. Deep spectra show that many of the sources in our sample are broad line radio galaxies (BLRG). The fact that the BLRG are observed out the redshift limit of the survey, overlapping in redshift with the quasars, argues against the idea that BLRG are simply the low radio power counterparts of high power, high redshift quasars.

We have obtained HST FOC f/48 long-slit spectroscopy of the inner 4" of the Narrow Line Region of NGC 1068 between 3500-5400A with a spectral resolution of 1.78A/pixel. At a spatial scale of 0.0287" per pixel these data provide an order of magnitude improvement in resolution over previous ground based spectra and allow us to trace the interaction between the radio jet and the gas in the NLR. Our results show that, within +/-0.5" of the radio-jet the emission lines are kinematically disturbed and split into two components whose velocity separation is 1500 km/sec. The filaments associated with the radio lobe also show a redshifted kinematic disturbance of the order of 300 km/sec which probably is a consequence of the expansion of the radio plasma. Furthermore, the material enveloping the radio-jet is in a much higher ionization state than that of the surrounding NLR gas. The highest excitation is coincident with the jet axis where emission in the coronal line of [FeVII] 3769A is detected and the HeII 4686A is strong but where [OII] 3727A is depressed. This large localized increase in ionization on the jet axis is accompanied by the presence of an excess continuum. Because the electron density is substantially larger in the jet compared to the surrounding NLR, these results can only be explained if there is a more intense ionizing continuum associated with the jet. This can be accomplished in a variety of ways which include an intrinsically anisotropic nuclear radiation field, a reduced gas covering factor or the presence of a local ionization source. The morphology, kinematics and, possibly, the ionization structure of the NLR in the vicinity of the jet of NGC 1068 are a direct consequence of the interaction with the radio outflow.

We have studied the core morphology of the ultraluminous infrared galaxy Mrk 273 by combining a high-resolution adaptive optics near-infrared image with an optical image from the Hubble Space Telescope and interferometric radio continuum data, all at spatial resolutions of 150 mas or better. The near-infrared image reveals that the nucleus has two main components, both of which have radio counterparts. The strongest component (N) shows very similar extended structure in the radio and near-infrared. It has a flat radio spectrum and is resolved into a double-lobed structure (Ne; Nw), with a separation of 90\pm5 mas (70 parsec). A similar structure is detected in the near-infrared. We identify this component as the location of the active nucleus. The second component (SW), strong in the near-infrared but relatively weak in the radio, is located $\sim1$ arcsecond to the southwest. We interpret this as an obscured starburst region associated with the merger. The radio continuum images show a third, strong, component (SE) which has previously been interpreted as a second nucleus. However, it shows no associated optical or near-infrared emission, suggesting that it is in fact a background source.