Source author record

Nadia L. Zakamska

Nadia L. Zakamska appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

52works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

52 published item(s)

preprint2022arXiv

Enhancement of Double-Close-Binary Quadruples

Double-close-binary quadruples (2+2 systems) are hierarchical systems of four stars where two short-period binary systems move around their common center of mass on a wider orbit. Using Gaia Early Data Release 3, we search for comoving pairs where both components are eclipsing binaries. We present eight 2+2 quadruple systems with inner orbital periods of $<$ 0.4 days and with outer separations of $\gtrsim 1000$ AU. All but one system are newly discovered by this work, and we catalog their orbital information measured from their light curves. We find that the occurrence rate of 2+2 quadruples is 7.3 $\pm$ 2.6 times higher than what is expected from random pairings of field stars. At most a factor of $\sim$2 enhancement may be explained by the age and metallicity dependence of the eclipsing binary fraction in the field stellar population. The remaining factor of $\sim$3 represents a genuine enhancement of the production of short-period binaries in wide-separation ($>10^3$ AU) pairs, suggesting a close binary formation channel that may be enhanced by the presence of wide companions.

preprint2022arXiv

Probing the Inner Circumgalactic Medium and Quasar Illumination around the Reddest `Extremely Red Quasar' (ERQ)

Dusty quasars might be in a young stage of galaxy evolution with prominent quasar feedback. A recently discovered population of luminous, extremely red quasars at $z\sim$~2--4 has extreme spectral properties related to exceptionally powerful quasar-driven outflows. We present Keck/KCWI observations of the reddest known ERQ, at $z=$\,2.3184, with extremely fast [\ion{O}{III}]~$λ$5007 outflow at $\sim$6000~km~s$^{-1}$. The Ly$α$ halo spans $\sim$100~kpc. The halo is kinematically quiet, with velocity dispersion $\sim$300~km~s$^{-1}$ and no broadening above the dark matter circular velocity down to the spatial resolution $\sim$6~kpc from the quasar. We detect spatially-resolved \ion{He}{II}~$λ$1640 and \ion{C}{IV}~$λ$1549 emissions with kinematics similar to the Ly$α$ halo and a narrow component in the [\ion{O}{III}]~$λ$5007. Quasar reddening acts as a coronagraph allowing views of the innermost halo. A narrow Ly$α$ spike in the quasar spectrum is inner halo emission, confirming the broad \ion{C}{IV}~$λ$1549 in the unresolved quasar is blueshifted by $2240$~km~s$^{-1}$ relative to the halo frame. We propose the inner halo is dominated by moderate-speed outflow driven in the past and the outer halo dominated by inflow. The high central concentration of the halo and the symmetric morphology of the inner region are consistent with the ERQ being in earlier evolutionary stage than blue quasars. The \ion{He}{II}~$λ$1640/Ly$α$ ratio of the inner halo and the asymmetry level of the overall halo are dissimilar to Type~II quasars, suggesting unique physical conditions for this ERQ that are beyond orientation differences from other quasar populations. We find no evidence of mechanical quasar feedback in the Ly$α$-emitting halo.

preprint2022arXiv

The eccentricity distribution of wide binaries and their individual measurements

Eccentricity of wide binaries is difficult to measure due to their long orbital periods. With Gaia's high-precision astrometric measurements, eccentricity of a wide binary can be constrained by the angle between the separation vector and the relative velocity vector (the $v$-$r$ angle). In this paper, by using the $v$-$r$ angles of wide binaries in Gaia Early Data Release 3, we develop a Bayesian approach to measure the eccentricity distribution as a function of binary separations. Furthermore, we infer the eccentricities of individual wide binaries and make them publicly available. Our results show that the eccentricity distribution of wide binaries at $10^2$ AU is close to uniform and becomes superthermal at $>10^{3}$ AU, suggesting two formation mechanisms dominating at different separation regimes. The close binary formation, most likely disk fragmentation, results in a uniform eccentricity distribution at $<10^{2}$ AU. The wide binary formation that leads to highly eccentric wide binaries at $>10^{3}$ AU may be turbulent fragmentation and/or the dynamical unfolding of compact triples. With Gaia, measuring eccentricities is now possible for a large number of wide binaries, opening a new window to understanding binary formation and evolution.

preprint2022arXiv

The SNIa Runaway LP 398-9: Detection of Circumstellar Material and Surface Rotation

A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back $\approx 10^5$ years to a known supernova remnant. Here we report the discovery of carbon-rich circumstellar material around LP 398-9, revealed by a strong infrared excess and analyzed with follow-up spectroscopy. The circumstellar material is most plausibly composed of inflated layers from the star itself, mechanically and radioactively heated by the past companion's supernova. We also detect a 15.4 hr periodic signal in the UV and optical light curves of LP 398-9, which we interpret as surface rotation. The rotation rate is consistent with theoretical predictions from this supernova mechanism, and the brightness variations could originate from surface inhomogeneity deposited by the supernova itself. Our observations strengthen the case for this double-degenerate SNIa progenitor channel, and motivate the search for more runaway SNIa donors.

preprint2022arXiv

Wide binaries from the H3 survey: the thick disk and halo have similar wide binary fractions

Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled $\sim$150,000 spectra for thick-disk and halo stars to date. We identify 800 high-confidence (a contamination rate of 4%) wide binaries and two resolved triples, with binary separations mostly between $10^3$-$10^5$ AU and a lowest [Fe/H] of $-2.7$. Based on their Galactic kinematics, 33 of them are halo wide binaries, and most of those are associated with the accreted Gaia-Sausage-Enceladus galaxy. The wide binary fraction in the thick disk decreases toward the low metallicity end, consistent with the previous findings for the thin disk. Our key finding is that the halo wide binary fraction is consistent with the thick-disk stars at a fixed [Fe/H]. There is no significant dependence of the wide binary fraction on the $α$-captured abundance. Therefore, the wide binary fraction is mainly determined by the iron abundance, not their disk or halo origin nor the $α$-captured abundance. Our results suggest that the formation environments play a major role for the wide binary fraction, instead of other processes like radial migration that only apply to disk stars.

preprint2021arXiv

Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA): Hubble Space Telescope Discovers Double Quasars

Dual supermassive black holes (SMBHs) at $\sim$kpc scales are the progenitor population of SMBH mergers and play an important role in understanding the pairing and dynamical evolution of massive black holes in galaxy mergers. Because of the stringent resolution requirement and the apparent rareness of these small-separation pairs, there are scarce observational constraints on this population, with few confirmed dual SMBHs at $<10$kpc separations at $z>1$. Here we present results from a pilot search for kpc-scale dual quasars selected with Gaia Data release 2 (DR2) astrometry and followed up with Hubble Space Telescope (HST) Wide Field Camera 3 dual-band (F475W and F814W) snapshot imaging. Our targets are quasars primarily selected with the varstrometry technique, i.e., light centroid jitter caused by asynchronous variability from both members in an unresolved quasar pair, supplemented by sub-arcsec pairs already resolved by Gaia DR2. We find an overall high fraction of HST-resolved pairs among the varstrometry-selected quasars (unresolved in Gaia DR2), $\sim 30-50\%$, increasing toward high redshift ($\sim 60-80\%$ at $z>1.5$). We discuss the nature of the 43 resolved sub-arcsec pairs based on HST and supplementary data. A substantial fraction ($\sim 40\%$) of these pairs are likely physical quasar pairs or gravitationally lensed quasars. We also discover a triple quasar candidate and a quadruply lensed quasar, which is among the smallest-separation quadruple lenses. These results provide important guidelines to improve varstrometry selection and follow-up confirmation of $\sim$kpc-scale dual SMBHs at high redshift.

preprint2021arXiv

X-ray analysis of SDSS J165202.60+172852.4, an obscured quasar with outflows at peak galaxy formation epoch

We report on deep XMM-Newton and NuSTAR observations of the high redshift, z=2.94, extremely red quasar (ERQ), SDSS J165202.60+172852.4, with known galactic ionized outflows detected via spatially-resolved [OIII] emission lines. X-ray observations allow us to directly probe the accretion disk luminosity and the geometry and scale of the circumnuclear obscuration. We fit the spectra from the XMM-Newton/EPIC and NuSTAR detectors with a physically motivated torus model and constrain the source to exhibit a near Compton-thick column density of $N_H=(1.02^{+0.76}_{-0.41})\times10^{24}\textrm{ cm}^{-2}$, a near edge-on geometry with the line-of-sight inclination angle of $θ_i=85^{\circ}$, and a scattering fraction of $f_{sc}\sim 3$ %. The absorption-corrected, intrinsic 2-10 keV X-ray luminosity of $L_{\textrm{2-10}}=(1.4^{+1}_{-1}) \times10^{45} \textrm{ erg s}^{-1}$ reveals a powerful quasar that is not intrinsically X-ray weak, consistent with observed trends in other ERQs. We also estimate the physical properties of the obscuration, although highly uncertain: the warm ionized scattering density of $n_e \sim 7.5\times(10^2-10^3)\textrm{ cm}^{-3}$ and the obscuration mass of $M_{obsc} \sim 1.7\times(10^4-10^6) M_{\odot}$. As previously suggested with shallower X-ray observations, optical and infrared selection of ERQ has proved effective in finding obscured quasars with powerful outflow signatures. Our observations provide an in-depth view into the X-ray properties of ERQs and support the conclusions of severely photon-limited studies of obscured quasar populations at high redshifts.

preprint2020arXiv

A Gravitational Redshift Measurement of the White Dwarf Mass-Radius Relation

The mass-radius relation of white dwarfs is largely determined by the equation of state of degenerate electrons, which causes the stellar radius to decrease as mass increases. Here we observationally measure this relation using the gravitational redshift effect, a prediction of general relativity that depends on the ratio between stellar mass and radius. Using observations of over three thousand white dwarfs from the Sloan Digital Sky Survey and the Gaia space observatory, we derive apparent radial velocities from absorption lines, stellar radii from photometry and parallaxes, and surface gravities by fitting atmospheric models to spectra. By averaging the apparent radial velocities of white dwarfs with similar radii and, independently, surface gravities, we cancel out random Doppler shifts and measure the underlying gravitational redshift. Using these results, we empirically measure the white dwarf mass-radius relation across a wide range of stellar masses. Our results are consistent with leading theoretical models, and our methods could be used with future observations to empirically constrain white dwarf core composition and evolution.

preprint2020arXiv

Computational Tools for the Spectroscopic Analysis of White Dwarfs

The spectroscopic features of white dwarfs are formed in the thin upper layer of their stellar photosphere. These features carry information about the white dwarf's surface temperature, surface gravity, and chemical composition (hereafter 'labels'). Existing methods to determine these labels rely on complex ab-initio theoretical models which are not always publicly available. Here we present two techniques to determine atmospheric labels from white dwarf spectra: a generative fitting pipeline that interpolates theoretical spectra with artificial neural networks, and a random forest regression model using parameters derived from absorption line features. We test and compare our methods using a large catalog of white dwarfs from the Sloan Digital Sky Survey (SDSS), achieving the same accuracy and negligible bias compared to previous studies. We package our techniques into an open-source Python module 'wdtools' that provides a computationally inexpensive way to determine stellar labels from white dwarf spectra observed from any facility. We will actively develop and update our tool as more theoretical models become publicly available. We discuss applications of our tool in its present form to identify interesting outlier white dwarf systems including those with magnetic fields, helium-rich atmospheres, and double-degenerate binaries.

preprint2020arXiv

Formation of Orion Fingers

"Orion fingers" are a system of dozens of bowshocks, with the wings of shocks pointing to a common system of origin, which is centered on a dynamically disintegrating system of several massive stars. The shock heads propagate with velocities of up to 300-400 km/s, but the formation and physical properties of the "bullets" leading the shocks are not known. Here we summarize two possible scenarios for the formation of the "bullets" and the resulting bowshocks ("fingers"). In the first scenario, bullets are self-gravitating, Jupiter-mass objects which were formed rapidly and then ejected during the strong dynamical interactions of massive stars and their disks. This scenario naturally explains the similar timescales for the outflow of bullets and for the dynamical interaction of the massive stars, but has some difficulty explaining the observed high velocities of the bullets. In the second scenario, bullets are formed via hydrodynamic instabilities in a massive, infrared-driven wind, naturally explaining the high velocities and the morphology of outflow, but the bullets are not required to be self-gravitating. The processes that created the Orion fingers are likely not unique to this particular star-forming region and may result in free-floating, high-velocity, core-less planets.

preprint2020arXiv

Ionised gas outflow signatures in SDSS-IV MaNGA active galactic nuclei

The prevalence of outflow and feedback signatures in active galactic nuclei (AGN) is a major unresolved question which large integral field unit (IFU) surveys now allow to address. In this paper, we present a kinematic analysis of the ionised gas in 2778 galaxies at z $\sim$ 0.05 observed by SDSS-IV MaNGA. Specifically, we measure the kinematics of the [OIII] λ5007Å emission line in each spatial element and fit multiple Gaussian components to account for possible non-gravitational motions of gas. Comparing the kinematics of the ionised gas between 308 MaNGA-selected AGN that have been previously identified through emission line diagnostics and sources not classified as AGN, we find that while 25% of MaNGA-selected AGN show [OIII] components with emission line widths of $>$ 500 km/s in more than 10% of their spaxels, only 7% of MaNGA non-AGN show a similar signature. Even the AGN that do not show nuclear AGN photoionisation signatures and that were only identified as AGN based on their larger scale photoionisation signatures show similar kinematic characteristics. In addition to obscuration, another possibility is that outflow and mechanical feedback signatures are longer lived than the AGN itself. Our measurements demonstrate that high velocity gas is more prevalent in AGN compared to non-AGN and that outflow and feedback signatures in low-luminosity, low-redshift AGN may so far have been underestimated. We show that higher luminosity MaNGA-selected AGN are able to drive larger scale outflows than lower luminosity AGN. But estimates of the kinetic coupling efficiencies are $<$ 1% and suggest that the feedback signatures probed in this paper are unlikely to have a significant impact on the AGN host galaxies. However, continuous energy injection may still heat a fraction of the cool gas and delay or suppress star formation in individual galaxies even when the AGN is weak.

preprint2020arXiv

Ionized and hot molecular outflows in the inner 500 pc of NGC1275

The role of feedback from Active Galactic Nuclei (AGN) in the evolution of galaxies is still not not fully understood, mostly due to the lack of observational constraints in the multi-phase gas kinematics on the ten to hundred parsec scales. We have used the Gemini Near-infrared Integral Field Spectrograph (NIFS) to map the molecular and ionized gas kinematics in the inner 900$\times$900 pc$^2$ of the Seyfert galaxy NGC1275 at a spatial resolution of $\sim$70 pc. From the fitting of the CO absorption bandheads in the K-band, we derive a stellar velocity dispersion of $265\pm26$ km s$^{-1}$, which implies a black hole mass of $M_{\rm SMBH}=1.1^{+0.9}_{-0.5}\times10^9$ M$_\odot$. We find hot ($T\gtrsim1000$ K) molecular and ionized outflows with velocities of up to 2 000 km s$^{-1}$ and mass outflow rates of $2.7\times10^{-2} {\rm M_\odot}$ yr$^{-1}$ and $1.6 {\rm M_\odot}$ yr$^{-1}$, respectively, in each of these gas phases. The kinetic power of the ionized outflows corresponds to only 0.05 per cent of the luminosity of the AGN of NGC 1275, indicating that they are not powerful enough to provide significant AGN feedback, but may be effective in redistributing the gas in the central region of the galaxy. The AGN driven outflows seem to be responsible for the shocks necessary to produce the observed H$_2$ and [Fe II] line emission.

preprint2020arXiv

Spatially-resolved UV diagnostics of AGN feedback: radiation pressure dominates in a prototypical quasar-driven superwind

Galactic-scale winds driven by active galactic nuclei (AGN) are often invoked to suppress star formation in galaxy evolution models, but the mechanisms driving these outflows are hotly debated. Two key AGN feedback models are (1) radiation pressure accelerating cool gas and (2) a hot outflowing wind entraining the ISM. Highly ionized emission-line diagnostics represent a powerful means of differentiating these scenarios because of their sensitivity to the expected compression of the ISM clouds by the hot wind. Here, we report the first spatially resolved UV emission spectroscopy of a prototypical (radio-quiet) quasar-driven superwind around the obscured quasar SDSSJ1356+1026 at z=0.123. We observe ratios of OVI/CIV, NV/CIV, and CIV/HeII that are remarkably similar for outflowing gas clouds <100 pc and ~10 kpc from the nucleus. Such similarity is expected for clouds with AGN radiation pressure dominated dynamics. Comparing the observed line emission to models of clouds in balance with radiation pressure and/or a hot wind, we rule out the presence of a dynamically important hot wind and constrain the ratio of hot gas pressure to radiation pressure to P_hot/P_rad<0.25 both at <100 pc and ~10 kpc from the nucleus. Moreover, the predictions of the radiation pressure confined cloud models that best fit observed UV line ratios are consistent with the observed diffuse X-ray spectrum. These results indicate that this AGN superwind is driven by radiation pressure or was driven by a hot wind that has since dissipated despite on-going AGN activity.

preprint2020arXiv

Very wide companion fraction from Gaia DR2: a weak or no enhancement for hot jupiter hosts, and a strong enhancement for contact binaries

There is an ongoing debate on whether hot jupiter hosts are more likely to be found in wide binaries with separations of $\gtrsim 100$ AU. In this paper, we search for comoving, very wide companions with separations of $10^3-10^4$ AU for hot jupiter hosts and main-sequence contact binaries in Gaia DR2, and compare the very wide companion fractions with their object-by-object-matched field star samples. We find that $11.9\pm 2.5$% of hot jupiter hosts and $14.1\pm 1.0$% of contact binaries have companions at separations of $10^3-10^4$ AU. While the very wide companion fraction of hot jupiter hosts is a factor of $1.9\pm0.5$ larger than their matched field star sample, it is consistent, within $\sim1σ$, with that of matched field stars if the matching is only with field stars without close companions (within $\sim50$ AU) as is the case for hot jupiter hosts. The very wide companion fraction of contact binaries is a factor of $3.1\pm0.5$ larger than their matched field star sample, suggesting that the formation and evolution of contact binaries are either tied to or correlated with the presence of wide companions. In contrast, the weak enhancement of very wide companion fraction for hot jupiter hosts implies that the formation of hot jupiters is not as sensitive to those environment properties. Our results also hint that the occurrence rates of dual hot jupiter hosts and dual contact binaries may be higher than the expected values from random pairing of field stars, which may be due to their underlying metallicity and age dependence.

preprint2019arXiv

Hubble Space Telescope Observations of Mira Variables in the Type Ia Supernova Host NGC 1559: An Alternative Candle to Measure the Hubble Constant

We present year-long, near-infrared Hubble Space Telescope WFC3 observations used to search for Mira variables in NGC 1559, the host galaxy of the Type Ia supernova (SN Ia) 2005df. This is the first dedicated search for Miras, highly-evolved low-mass stars, in a SN Ia host and subsequently the first calibration of the SN Ia luminosity using Miras in a role historically played by Cepheids. We identify a sample of 115 O-rich Miras with P < 400 days based on their light curve properties. We find that the scatter in the Mira Period-Luminosity Relation (PLR) is comparable to Cepheid PLRs seen in SN Ia supernova host galaxies. Using a sample of O-rich Miras discovered in NGC 4258 with HST F160W and its maser distance, we measure a distance modulus for NGC 1559 of mu1559 = 31.41 +/- 0.050 (statistical) +/- 0.060 (systematic) mag. Based on the light curve of the normal, well-observed, low-reddening SN 2005df, we obtain a measurement of the fiducial SN Ia absolute magnitude of MB0 = -19.27 +/- 0.13 mag. With the Hubble diagram of SNe Ia we find H0 = 72.7 +/- 4.6 kms-1 Mpc-1. Combining the calibration from the NGC 4258 megamaser and the Large Magellanic Cloud detached eclipsing binaries gives a best value of H0 = 73.3 +/- 4.0 km s-1 Mpc-1. This result is within 1-sigma of the Hubble constant derived using Cepheids and multiple calibrating SNe Ia. This is the first of four expected calibrations of the SN Ia luminosity from Miras which should reduce the error in H0 via Miras to ~3%. In light of the present Hubble tension and JWST, Miras have utility in the extragalactic distance scale to check Cepheid distances or calibrate nearby SNe in early-type host galaxies that would be unlikely targets for Cepheid searches.

preprint2016arXiv

Constraining the dynamical importance of hot gas and radiation pressure in quasar outflows using emission line ratios

Quasar feedback models often predict an expanding hot gas bubble which drives a galaxy-scale outflow. In many circumstances this hot gas radiates inefficiently and is therefore difficult to observe directly. We present an indirect method to detect the presence of a hot bubble using hydrostatic photoionization calculations of the cold (~10^4 K) line-emitting gas. We compare our calculations with observations of the broad line region, the inner face of the torus, the narrow line region (NLR), and the extended NLR, and thus constrain the hot gas pressure at distances 0.1 pc -- 10 kpc from the center. We find that emission line ratios observed in the average quasar spectrum are consistent with radiation-pressure-dominated models on all scales. On scales <40 pc a dynamically significant hot gas pressure is ruled out, while on larger scales the hot gas pressure cannot exceed six times the local radiation pressure. In individual quasars, ~25% of quasars exhibit NLR ratios that are inconsistent with radiation-pressure-dominated models, though in these objects the hot gas pressure is also unlikely to exceed the radiation pressure by an order of magnitude or more. The derived upper limits on the hot gas pressure imply that the instantaneous gas pressure force acting on galaxy-scale outflows falls short of the time-averaged force needed to explain the large momentum fluxes \dot{p} >> L_AGN/c inferred for galaxy-scale outflows. This apparent discrepancy can be reconciled if optical quasars previously experienced a buried, fully-obscured phase during which the hot gas bubble was more effectively confined and during which galactic wind acceleration occurred.

preprint2016arXiv

Discovery of extreme [OIII]5007A outflows in high-redshift red quasars

Black hole feedback is now a standard component of galaxy formation models. These models predict that the impact of black hole activity on its host galaxy likely peaked at z=2-3, the epoch of strongest star formation activity and black hole accretion activity in the Universe. We used XShooter on the Very Large Telescope to measure rest-frame optical spectra of four z~2.5 extremely red quasars with infrared luminosities ~10^47 erg/sec. We present the discovery of very broad (full width at half max= 2600-5000 km/sec), strongly blue-shifted (by up to 1500 km/sec) [OIII]5007A emission lines in these objects. In a large sample of obscured and red quasars, [OIII] kinematics are positively correlated with infrared luminosity, and the four objects in our sample are on the extreme end both in [OIII] kinematics and infrared luminosity. We estimate that ~3% of the bolometric luminosity in these objects is being converted into the kinetic power of the observed wind. Photo-ionization estimates suggest that the [OIII] emission might be extended on a few kpc scales, which would suggest that the extreme outflow is affecting the entire host galaxy of the quasar. These sources may be the signposts of the most extreme form of quasar feedback at the peak epoch of galaxy formation, and may represent an active "blow-out" phase of quasar evolution.

preprint2016arXiv

Do galaxy global relationships emerge from local ones? The SDSS IV MaNGA surface mass density-metallicity relation

We present the stellar surface mass density {\it vs.} gas metallicity ($Σ_*-Z$) relation for more than 500,000 spatially-resolved star-forming resolution elements (spaxels) from a sample of 653 disk galaxies included in the SDSS IV MaNGA survey. We find a tight relation between these local properties, with higher metallicities as the surface density increases. This relation extends over three orders of magnitude in the surface mass density and a factor of four in metallicity. We show that this local relationship can simultaneously reproduce two well-known properties of disk galaxies: their global mass-metallicity relationship {\it and} their radial metallicity gradients. We also find that the $Σ_* - Z$ relation is largely independent of the galaxy's total stellar mass and specific star-formation rate (sSFR), except at low stellar mass and high sSFR. These results suggest that in the present-day universe local properties play a key role in determining the gas-phase metallicity in typical disk galaxies.

preprint2016arXiv

Evidence for the Thermal Sunyaev-Zel'dovich Effect Associated with Quasar Feedback

Using a radio-quiet subsample of the Sloan Digital Sky Survey spectroscopic quasar catalogue, spanning redshifts 0.5-3.5, we derive the mean millimetre and far-infrared quasar spectral energy distributions (SEDs) via a stacking analysis of Atacama Cosmology Telescope and Herschel-Spectral and Photometric Imaging REceiver data. We constrain the form of the far-infrared emission and find 3$σ$-4$σ$ evidence for the thermal Sunyaev-Zel'dovich (SZ) effect, characteristic of a hot ionized gas component with thermal energy $(6.2 \pm 1.7)\times 10^{60}$ erg. This amount of thermal energy is greater than expected assuming only hot gas in virial equilibrium with the dark matter haloes of $(1-5)\times 10^{12}h^{-1}$M$_\odot$ that these systems are expected to occupy, though the highest quasar mass estimates found in the literature could explain a large fraction of this energy. Our measurements are consistent with quasars depositing up to $(14.5 \pm 3.3)~τ_8^{-1}$ per cent of their radiative energy into their circumgalactic environment if their typical period of quasar activity is $τ_8\times~10^8$ yr. For high quasar host masses, $\sim10^{13}h^{-1}$M$_\odot$, this percentage will be reduced. Furthermore, the uncertainty on this percentage is only statistical and additional systematic uncertainties enter at the 40 per cent level. The SEDs are dust dominated in all bands and we consider various models for dust emission. While sufficiently complex dust models can obviate the SZ effect, the SZ interpretation remains favoured at the 3$σ$-4$σ$ level for most models.

preprint2016arXiv

Evidence of suppression of star formation by quasar-driven winds in gas-rich host galaxies at z<1?

Feedback from active galactic nuclei (AGN) is widely considered to be the main driver in regulating the growth of massive galaxies through heating or driving gas out of the galaxy, preventing further increase in stellar mass. Observational proof for this scenario has, however, been scarce. We have assembled a sample of 132 radio-quiet type-2 and red AGN at 0.1<z<1. We measure the kinematics of the AGN-ionized gas, the host galaxies' stellar masses and star formation rates and investigate the relationships between AGN luminosities, specific star formation rates (sSFR) and outflow strengths W90 -- the 90\% velocity width of the [OIII]5007 line power and a proxy for the AGN-driven outflow speed. Outflow strength is independent of sSFR for AGN selected on their mid-IR luminosity, in agreement with previous work demonstrating that star formation is not sufficient to produce the observed ionized gas outflows which have to be powered by AGN activity. More importantly, we find a negative correlation between W90 and sSFR in the AGN hosts with the highest SFRs, i.e., with the highest gas content, where presumably the coupling of the AGN-driven wind to the gas is strongest. This implies that AGN with strong outflow signatures are hosted in galaxies that are more `quenched' than galaxies with weaker outflow signatures. Despite the galaxies' high SFRs, we demonstrate that the outflows are not star-formation driven but indeed due to AGN-powering. This observation is consistent with the AGN having a net suppression, `negative' impact, through feedback on the galaxies' star formation history.

preprint2016arXiv

Extremely Red Quasars from SDSS, BOSS and WISE: Classification of Optical Spectra

Quasars with extremely red infrared-to-optical colours are an interesting population that can test ideas about quasar evolution as well as orientation, obscuration and geometric effects in the so-called AGN unified model. To identify such a population we match the quasar catalogues of the Sloan Digital Sky Survey (SDSS), the Baryon Oscillation Spectroscopic Survey (BOSS) to the Wide-Field Infrared Survey Explorer (WISE) to identify quasars with extremely high infrared-to-optical ratios. We identify 65 objects with r(AB)-W4(Vega)>14 mag (i.e., F_nu(22um)/F_nu(r) > ~1000). This sample spans a redshift range of 0.28<z<4.36 and has a bimodal distribution, with peaks at z~0.8 and z~2.5. It includes three z>2.6 objects that are detected in the W4-band but not W1 or W2 (i.e., W1W2-dropouts). The SDSS/BOSS spectra show that the majority of the objects are reddened Type 1 quasars, Type 2 quasars (both at low and high redshift) or objects with deep low-ionization broad absorption lines (BALs) that suppress the observed r-band flux. In addition, we identify a class of Type 1 permitted broad-emission line objects at z~2-3 which are characterized by emission line rest-frame equivalent widths (REWs) of >~150Ang , much larger than those of typical quasars. In particular, 55% (45%) of the non-BAL Type 1s with measurable CIV in our sample have REW(CIV) > 100 (150)Ang, compared to only 5.8% (1.3%) for non-BAL quasars in BOSS. These objects often also have unusual line ratios, such as very high NV/Ly-alpha ratios. These large REWs might be caused by suppressed continuum emission analogous to Type 2 quasars; however, there is no obvious mechanism in standard Unified Models to suppress the continuum without also obscuring the broad emission lines.

preprint2016arXiv

Sensitive Radio Survey of Obscured Quasar Candidates

We study the radio properties of moderately obscured quasars over a range of redshifts to understand the role of radio activity in accretion using the Jansky Very Large Array (JVLA) at 6.0GHz and 1.4GHz. Our z~2.5 sample consists of optically-selected obscured quasar candidates, all of which are radio-quiet, with typical radio luminosities of $νL_ν$[1.4 GHz] < $10^{40}$ erg s$^{-1}$. Only a single source is individually detected in our deep (rms~10 $μ$Jy) exposures. This population would not be identified by radio-based selection methods used for distinguishing dusty star-forming galaxies and obscured active nuclei. In our pilot A-array study of z~0.5 radio-quiet quasars, we spatially resolve four of five objects on scales ~ 5 kpc and find they have steep spectral indices. Therefore, radio emission in these sources could be due to jet-driven or radiatively driven bubbles interacting with interstellar material on the scale of the host galaxy. Finally, we also study the population of ~ 200 faint (~40 $μ$Jy - 40 mJy) radio sources observed over ~ 120 arcmin$^2$ of our data. 60% of these detections are matched in the SDSS and/or WISE and are, in roughly equal shares, active nuclei at a broad range of redshifts, passive galaxies with no other signs of nuclear activity and IR-bright but optically faint sources. Spectroscopically or photometrically confirmed star-forming galaxies constitute only a small minority of the matches. Such sensitive radio surveys allow us to address important questions of AGN evolution and evaluate the AGN contribution to the radio-quiet sky.

preprint2016arXiv

SpIES: The Spitzer IRAC Equatorial Survey

We describe the first data release from the Spitzer-IRAC Equatorial Survey (SpIES); a large-area survey of 115 deg^2 in the Equatorial SDSS Stripe 82 field using Spitzer during its 'warm' mission phase. SpIES was designed to probe sufficient volume to perform measurements of quasar clustering and the luminosity function at z > 3 to test various models for "feedback" from active galactic nuclei (AGN). Additionally, the wide range of available multi-wavelength, multi-epoch ancillary data enables SpIES to identify both high-redshift (z > 5) quasars as well as obscured quasars missed by optical surveys. SpIES achieves 5σ depths of 6.13 μJy (21.93 AB magnitude) and 5.75 μJy (22.0 AB magnitude) at 3.6 and 4.5 microns, respectively - depths significantly fainter than WISE. We show that the SpIES survey recovers a much larger fraction of spectroscopically-confirmed quasars (98%) in Stripe 82 than are recovered by WISE (55%). This depth is especially powerful at high-redshift (z > 3.5), where SpIES recovers 94% of confirmed quasars, whereas WISE only recovers 25%. Here we define the SpIES survey parameters and describe the image processing, source extraction, and catalog production methods used to analyze the SpIES data. In addition to this survey paper, we release 234 images created by the SpIES team and three detection catalogs: a 3.6 μm-only detection catalog containing 6.1 million sources, a 4.5 μm-only detection catalog containing 6.5 million sources, and a dual-band detection catalog containing 5.4 million sources.

preprint2016arXiv

Towards a comprehensive picture of powerful quasars, their host galaxies and quasar winds at z ~ 0.5

Luminous type-2 quasars in which the glow from the central black hole is obscured by dust are ideal targets for studying their host galaxies and the quasars' effect on galaxy evolution. Such feedback appears ubiquitous in luminous obscured quasars where high velocity ionized nebulae have been found. We present rest-frame yellow-band (~5000 Angstroms) observations using the Hubble Space Telescope for a sample of 20 luminous quasar host galaxies at 0.2 < z < 0.6 selected from the Sloan Digital Sky Survey. For the first time, we combine host galaxy observations with geometric measurements of quasar illumination using blue-band HST observations and [OIII] integral field unit observations probing the quasar winds. The HST images reveal bright merger signatures in about half the galaxies; a significantly higher fraction than in comparison inactive ellipticals. We show that the host galaxies are primarily bulge-dominated, with masses close to M*, but belong to < 30% of elliptical galaxies that are highly star-forming at z ~ 0.5. Ionized gas signatures are uncorrelated with faint stellar disks (if present), confirming that the ionized gas is not concentrated in a disk. Scattering cones and [OIII] ionized gas velocity field are aligned with the forward scattering cones being co-spatial with the blue-shifted side of the velocity field, suggesting the high velocity gas is indeed photo-ionized by the quasar. Based on the host galaxies' high star-formation rates and bright merger signatures, we suggest that this low-redshift outbreak of luminous quasar activity is triggered by recent minor mergers. Combining these novel observations, we present new quasar unification tests, which are in agreement with expectations of the orientation-based unification model for quasars.

preprint2015arXiv

Giant scattering cones in obscured quasars

We analyze Hubble Space Telescope observations of scattering regions in 20 luminous obscured quasars at $0.24<z<0.65$ (11 new observations and 9 archival ones) observed at rest-frame $\sim 3000$Å. We find spectacular $5-10$ kpc-scale scattering regions in almost all cases. The median scattering efficiency at this wavelength (the ratio of observed to estimated intrinsic flux) is 2.3\%, and 73\% of the observed flux at this wavelength is due to scattered light, which if unaccounted for may strongly bias estimates of quasar hosts' star formation rates. Modeling these regions as illuminated dusty cones, we estimate the radial density distributions of the interstellar medium as well as the geometric properties of circumnuclear quasar obscuration -- inclinations and covering factors. Small derived opening angles (median half-angle and standard deviation 27\dg$\pm$9\dg) are inconsistent with a 1:1 type 1 / type 2 ratio. We suggest that quasar obscuration is patchy and that the observer has a $\sim 40\%$ chance of seeing a type 1 source even through the obscuration. We estimate median density profile of the scattering medium to be $n_{\rm H}=0.04-0.5$ $(1{\rm kpc}/r)^2$ cm$^{-3}$, depending on the method. Quasars in our sample likely exhibit galaxy-wide winds, but if these consist of optically thick clouds then only a small fraction of the wind mass ($\la 10\%$) contributes to scattering.

preprint2015arXiv

Merger-driven Fueling of Active Galactic Nuclei: Six Dual and Offset Active Galactic Nuclei Discovered with Chandra and Hubble Space Telescope Observations

Dual active galactic nuclei (AGNs) and offset AGNs are kpc-scale separation supermassive black holes pairs created during galaxy mergers, where both or one of the black holes are AGNs, respectively. These dual and offset AGNs are valuable probes of the link between mergers and AGNs but are challenging to identify. Here we present Chandra/ACIS observations of 12 optically-selected dual AGN candidates at z < 0.34, where we use the X-rays to identify AGNs. We also present HST/WFC3 observations of 10 of these candidates, which reveal any stellar bulges accompanying the AGNs. We discover a dual AGN system with separation of 2.2 kpc, where the two stellar bulges have coincident [O III] and X-ray sources. This system is an extremely minor merger (460:1) that may include a dwarf galaxy hosting an intermediate mass black hole. We also find six single AGNs, and five systems that are either dual or offset AGNs with separations < 10 kpc. Four of the six dual AGNs and dual/offset AGNs are in ongoing major mergers, and these AGNs are 10 times more luminous, on average, than the single AGNs in our sample. This hints that major mergers may preferentially trigger higher luminosity AGNs. Further, we find that confirmed dual AGNs have hard X-ray luminosities that are half of those of single AGNs at fixed [O III] luminosity, on average. This could be explained by high densities of gas funneled to galaxy centers during mergers, and emphasizes the need for deeper X-ray observations of dual AGN candidates.

preprint2015arXiv

Similarity of ionized gas nebulae around unobscured and obscured quasars

Quasar feedback is suspected to play a key role in the evolution of massive galaxies, by removing or reheating gas in quasar host galaxies and thus limiting the amount of star formation. In this paper we continue our investigation of quasar-driven winds on galaxy-wide scales. We conduct Gemini Integral Field Unit spectroscopy of a sample of luminous unobscured (type 1) quasars, to determine the morphology and kinematics of ionized gas around these objects, predominantly via observations of the [O III]5007 emission line. We find that ionized gas nebulae extend out to ~13 kpc from the quasar, that they are smooth and round, and that their kinematics are inconsistent with gas in dynamical equilibrium with the host galaxy. The observed morphological and kinematic properties are strikingly similar to those of ionized gas around obscured (type 2) quasars with matched [O III] luminosity, with marginal evidence that nebulae around unobscured quasars are slightly more compact. Therefore in samples of obscured and unobscured quasars carefully matched in [O III] luminosity we find support for the standard geometry-based unification model of active galactic nuclei, in that the intrinsic properties of quasars, of their hosts and of their ionized gas appear to be very similar. Given the apparent ubiquity of extended ionized regions, we are forced to conclude that either the quasar is at least partially illuminating pre-existing gas or that both samples of quasars are seen during advanced stages of quasar feedback. In the latter case, we may be biased by our [O III]-based selection against quasars in the early "blow-out" phase, for example due to dust obscuration.

preprint2015arXiv

Star formation in quasar hosts and the origin of radio emission in radio-quiet quasars

Radio emission from radio-quiet quasars may be due to star formation in the quasar host galaxy, to a jet launched by the supermassive black hole, or to relativistic particles accelerated in a wide-angle radiatively-driven outflow. In this paper we examine whether radio emission from radio-quiet quasars is a byproduct of star formation in their hosts. To this end we use infrared spectroscopy and photometry from Spitzer and Herschel to estimate or place upper limits on star formation rates in hosts of ~300 obscured and unobscured quasars at z<1. We find that low-ionization forbidden emission lines such as [NeII] and [NeIII] are likely dominated by quasar ionization and do not provide reliable star formation diagnostics in quasar hosts, while PAH emission features may be suppressed due to the destruction of PAH molecules by the quasar radiation field. While the bolometric luminosities of our sources are dominated by the quasars, the 160 micron fluxes are likely dominated by star formation, but they too should be used with caution. We estimate median star formation rates to be 6-29 Msun/year, with obscured quasars at the high end of this range. This star formation rate is insufficient to explain the observed radio emission from quasars by an order of magnitude, with log(L_radio, observed/L_radio, SF)=0.6-1.3 depending on quasar type and star formation estimator. Although radio-quiet quasars in our sample lie close to the 8-1000 micron infrared / radio correlation characteristic of the star-forming galaxies, both their infrared emission and their radio emission are dominated by the quasar activity, not by the host galaxy.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

preprint2014arXiv

A Glimpse at Quasar Host Galaxy Far-UV Emission, Using DLAs as Natural Coronagraphs

In merger-driven models of massive galaxy evolution, the luminous quasar phase is expected to be accompanied by vigorous star formation in quasar host galaxies. In this paper, we use high column density Damped Lyman Alpha (DLA) systems along quasar sight lines as natural coronagraphs to directly study the far-UV (FUV) radiation from the host galaxies of luminous background quasars. We have stacked the spectra of $\sim$2,000 DLA systems (N_HI>10^{20.6} cm^{-2}) with a median absorption redshift < z > = 2.6 selected from quasars observed in the SDSS-III Baryon Oscillation Spectroscopic Survey. We detect residual flux in the dark troughs of the composite DLA spectra. The level of this residual flux significantly exceeds systematic errors in the SDSS fiber sky subtraction; furthermore, the residual flux is strongly correlated with the continuum luminosity of the background quasar, while uncorrelated with DLA column density or metallicity. We conclude that the flux could be associated with the average FUV radiation from the background quasar host galaxies (with medium redshift < z > = 3.1) that is not blocked by the intervening DLA. Assuming all of the detected flux originates from quasar hosts, for the highest quasar luminosity bin (< L >= 2.5x 10^{13} L_sun), the host galaxy has a FUV intensity of 1.5 +/- 0.2 x 10^{40} erg s^{-1} A^{-1}; this corresponds to an unobscured UV star formation rate of 9 M_sun/yr.

preprint2014arXiv

ALMA Observations of a Candidate Molecular Outflow in an Obscured Quasar

We present Atacama Large Millimeter/Submillimeter Array (ALMA) CO (1-0) and CO (3-2) observations of SDSS J135646.10+102609.0, an obscured quasar and ultra-luminous infrared galaxy (ULIRG) with two merging nuclei and a known 20-kpc-scale ionized outflow. The total molecular gas mass is M_{mol} ~ 9^{+19}_{-6} x 10^8 Msun, mostly distributed in a compact rotating disk at the primary nucleus (M_{mol} ~ 3 x 10^8 Msun) and an extended tidal arm (M_{mol} ~ 5 x 10^8 Msun). The tidal arm is one of the most massive molecular tidal features known; we suggest that it is due to the lower chance of shock dissociation in this elliptical/disk galaxy merger. In the spatially resolved CO (3-2) data, we find a compact (r ~ 0.3 kpc) high velocity (v ~ 500 km/s) red-shifted feature in addition to the rotation at the N nucleus. We propose a molecular outflow as the most likely explanation for the high velocity gas. The outflowing mass of M_{mol} ~ 7 x 10^7 Msun and the short dynamical time of t_{dyn} ~ 0.6 Myr yield a very high outflow rate of \dot{M}_{mol} ~ 350 Msun/yr and can deplete the gas in a million years. We find a low star formation rate (< 16 Msun/yr from the molecular content and < 21 Msun/yr from the far-infrared spectral energy distribution decomposition) that is inadequate to supply the kinetic luminosity of the outflow (\dot{E} ~ 3 x 10^43 erg/s). Therefore, the active galactic nucleus, with a bolometric luminosity of 10^46 erg/s, likely powers the outflow. The momentum boost rate of the outflow (\dot{p}/(Lbol/c) ~ 3) is lower than typical molecular outflows associated with AGN, which may be related to its compactness. The molecular and ionized outflows are likely two distinct bursts induced by episodic AGN activity that varies on a time scale of 10^7 yr.

preprint2014arXiv

Extended X-ray Emission From a Quasar-Driven Superbubble

We present observations of extended, 20-kpc scale soft X-ray gas around a luminous obscured quasar hosted by an ultra-luminous infrared galaxy caught in the midst of a major merger. The extended X-ray emission is well fit as a thermal gas with a temperature of kT ~ 280 eV and a luminosity of L_X ~ 10^42 erg/s and is spatially coincident with a known ionized gas outflow. Based on the X-ray luminosity, a factor of ~10 fainter than the [OIII] emission, we conclude that the X-ray emission is either dominated by photoionization, or by shocked emission from cloud surfaces in a hot quasar-driven wind.

preprint2014arXiv

Gemini Long-slit Observations of Luminous Obscured Quasars: Further Evidence for an Upper Limit on the Size of the Narrow-Line Region

We examine the spatial extent of the narrow-line regions (NLRs) of a sample of 30 luminous obscured quasars at $0.4 < z < 0.7$ observed with spatially resolved Gemini-N GMOS long-slit spectroscopy. Using the [OIII]$\lambda5007$ emission feature, we estimate the size of the NLR using a cosmology-independent measurement: the radius where the surface brightness falls to 10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We then explore the effects of atmospheric seeing on NLR size measurements and conclude that direct measurements of the NLR size from observed profiles are too large by 0.1 - 0.2 dex on average, as compared to measurements made to best-fit Sérsic or Voigt profiles convolved with the seeing. These data, which span a full order of magnitude in IR luminosity ($\log{(L_{8 μ\mathrm{m}} / \mathrm{erg\, s}^{-1})} = 44.4 - 45.4$) also provide strong evidence that there is a flattening of the relationship between NLR size and AGN luminosity at a seeing-corrected size of $\sim 7$ kpc. The objects in this sample have high luminosities which place them in a previously under-explored portion of the size-luminosity relationship. These results support the existence of a maximal size of the narrow-line region around luminous quasars; beyond this size either there is not enough gas, or the gas is over-ionized and does not produce enough [OIII]$\lambda5007$ emission.

preprint2014arXiv

Quasar feedback and the origin of radio emission in radio-quiet quasars

We conduct kinematic analysis of the SDSS spectra of 568 obscured luminous quasars, with the emphasis on the kinematic structure of the [OIII]5007 emission line. [OIII] emission tends to show blueshifts and blue excess, which indicates that at least part of the narrow-line gas is undergoing an organized outflow. The velocity width containing 90% of line power ranges from 370 to 4780 km/sec, suggesting outflow velocities up to 2000 km/sec. The velocity width of the [OIII] emission is positively correlated with the radio luminosity among the radio-quiet quasars. We propose that radio emission in radio-quiet quasars is due to relativistic particles accelerated in the shocks within the quasar-driven outflows; star formation in quasar hosts is insufficient to explain the observed radio emission. The median radio luminosity of the sample of nu L_nu[1.4GHz] = 10^40 erg/sec suggests a median kinetic luminosity of the quasar-driven wind of L_wind=3x10^44 erg/sec, or about 4% of the estimated median bolometric luminosity L_bol=8x10^45 erg/sec. Furthermore, the velocity width of [OIII] is positively correlated with mid-infrared luminosity, which suggests that outflows are ultimately driven by the radiative output of the quasar. As the outflow velocity increases, some emission lines characteristic of shocks in quasi-neutral medium increase as well, which we take as further evidence of quasar-driven winds propagating into the interstellar medium of the host galaxy. None of the kinematic components show correlations with the stellar velocity dispersions of the host galaxies, so there is no evidence that any of the gas in the narrow-line region of quasars is in dynamical equilibrium with the host galaxy. Quasar feedback appears to operate above the threshold luminosity of L_bol=3x10^45 erg/sec.

preprint2014arXiv

Warm molecular hydrogen in outflows from Ultraluminous Infrared Galaxies

Ultraluminous infrared galaxies (ULIRGs) show on average three times more emission in the rotational transitions of molecular hydrogen than expected based on their star formation rates. Using Spitzer archival data we investigate the origin of excess warm H_2 emission in 115 ULIRGs of the IRAS 1 Jy sample. We find a strong correlation between H_2 and [FeII] line luminosities, suggesting that excess H_2 is produced in shocks propagating within neutral or partially ionized medium. This view is supported by the correlations between H_2 and optical line ratios diagnostic of such shocks. The galaxies powered by star formation and those powered by active nuclei follow the same relationship between H_2 and [FeII], with emission line width being the major difference between these classes (about 500 and 1000 km/sec, respectively). We conclude that excess H_2 emission is produced as the supernovae and active nuclei drive outflows into the neutral interstellar medium of the ULIRGs. A weak positive correlation between H_2 and the length of the tidal tails indicates that these outflows are more likely to be encountered in more advanced mergers, but there is no evidence for excess H_2 produced as a result of the collision shocks during the final coalescence.

preprint2013arXiv

Candidate Type II Quasars at 2 < z < 4.3 in the Sloan Digital Sky Survey III

At low redshifts, dust-obscured quasars often have strong yet narrow permitted lines in the rest-frame optical and ultraviolet, excited by the central active nucleus, earning the designation Type II quasars. We present a sample of 145 candidate Type II quasars at redshifts between 2 and 4.3, encompassing the epoch at which quasar activity peaked in the universe. These objects, selected from the quasar sample of the Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey III, are characterized by weak continuum in the rest-frame ultraviolet (typical continuum magnitude of i \approx 22) and strong lines of CIV and Ly α, with Full Width at Half Maximum less than 2000 kms-1. The continuum magnitudes correspond to an absolute magnitude of -23 or brighter at redshift 3, too bright to be due exclusively to the host galaxies of these objects. Roughly one third of the objects are detected in the shorter-wavelength bands of the WISE survey; the spectral energy distributions (SEDs) of these objects appear to be intermediate between classic Type I and Type II quasars seen at lower redshift. Five objects are detected at rest frame 6μm by Spitzer, implying bolometric luminosities of several times 10^46 erg s-1. We have obtained polarization measurements for two objects; they are roughly 3% polarized. We suggest that these objects are luminous quasars, with modest dust extinction (A_V ~ 0.5 mag), whose ultraviolet continuum also includes a substantial scattering contribution. Alternatively, the line of sight to the central engines of these objects may be partially obscured by optically thick material.

preprint2013arXiv

Observations of Feedback from Radio-Quiet Quasars - II. Kinematics of Ionized Gas Nebulae

The prevalence and energetics of quasar feedback is a major unresolved problem in galaxy formation theory. In this paper, we present Gemini Integral Field Unit observations of ionized gas around eleven luminous, obscured, radio-quiet quasars at z~0.5 out to ~15 kpc from the quasar; specifically, we measure the kinematics and morphology of [O III]5007 emission. The round morphologies of the nebulae and the large line-of-sight velocity widths (with velocities containing 80% of the emission as high as 1000 km/s combined with relatively small velocity difference across them (from 90 to 520 km/s) point toward wide-angle quasi-spherical outflows. We use the observed velocity widths to estimate a median outflow velocity of 760 km/s, similar to or above the escape velocities from the host galaxies. The line-of-sight velocity dispersion declines slightly toward outer parts of the nebulae (by 3% per kpc on average). The majority of nebulae show blueshifted excesses in their line profiles across most of their extents, signifying gas outflows. For the median outflow velocity, we find a kinetic energy flow between 4x10^{44} and 3x10^{45} erg/s and mass outflow rate between 2000 and 20000 Msun/yr. These values are large enough for the observed quasar winds to have a significant impact on their host galaxies. The median rate of converting bolometric luminosity to kinetic energy of ionized gas clouds is ~2%. We report four new candidates for "super-bubbles" -- outflows that may have broken out of the denser regions of the host galaxy.

preprint2013arXiv

Observations of Feedback from Radio-Quiet Quasars: I. Extents and Morphologies of Ionized Gas Nebulae

Black hole feedback -- the strong interaction between the energy output of supermassive black holes and their surrounding environments -- is routinely invoked to explain the absence of overly luminous galaxies, the black hole vs. bulge correlations and the similarity of black hole accretion and star formation histories. Yet direct probes of this process in action are scarce and limited to small samples of active nuclei. We present Gemini IFU observations of the distribution of ionized gas around luminous, obscured, radio-quiet (RQ) quasars at z~0.5. We detect extended ionized gas nebulae via [O III]5007 emission in every case, with a mean diameter of 28 kpc. These nebulae are nearly perfectly round. The regular morphologies of nebulae around RQ quasars are in striking contrast with lumpy or elongated nebulae seen around radio galaxies at low and high redshifts. We present the uniformly measured size-luminosity relationship of [O III] nebulae around Seyfert 2 galaxies and type 2 quasars spanning 6 orders of magnitude in luminosity and confirm the flat slope of the correlation (R ~ L^{0.25+/-0.02}). We find a universal behavior of the [O III]/H-beta ratio in our entire RQ quasar sample: it persists at a constant value (~10) in the central regions, until reaching a "break" isophotal radius ranging from 4 to 11 kpc where it starts to decrease. We propose a model of clumpy nebulae in which clouds that produce line emission transition from being ionization-bounded at small distances from the quasar to being matter-bounded in the outer parts of the nebula, which qualitatively explains the observed line ratio and surface brightness profiles. It is striking that we see such smooth and round large-scale gas nebulosities in this sample, which are inconsistent with illuminated merger debris and which we suggest may be the signature of accretion energy from the nucleus reaching gas at large scales.

preprint2013arXiv

SALT Long-slit Spectroscopy of Luminous Obscured Quasars: An Upper Limit on the Size of the Narrow-Line Region?

We present spatially resolved long-slit spectroscopy from the Southern African Large Telescope (SALT) to examine the spatial extent of the narrow-line regions (NLRs) of a sample of 8 luminous obscured quasars at 0.10 < z < 0.43. Our results are consistent with an observed shallow slope in the relationship between NLR size and L_[OIII], which has been interpreted to indicate that NLR size is limited by the density and ionization state of the NLR gas rather than the availability of ionizing photons. We also explore how the NLR size scales with a more direct measure of instantaneous AGN power using mid-IR photometry from WISE, which probes warm to hot dust near the central black hole and so, unlike [OIII], does not depend on the properties of the NLR. Using our results as well as samples from the literature, we obtain a power-law relationship between NLR size and L_8micron that is significantly steeper than that observed for NLR size and L_[OIII]. We find that the size of the NLR goes approximately as L^(1/2)_8micron, as expected from the simple scenario of constant-density clouds illuminated by a central ionizing source. We further see tentative evidence for a flattening of the relationship between NLR size and L_8micron at the high luminosity end, and propose that we are seeing a limiting NLR size of 10 - 20 kpc, beyond which the availability of gas to ionize becomes too low. We find that L_[OIII] ~ L_8micron^(1.4), consistent with a picture in which the L_[OIII] is dependent on the volume of the NLR. These results indicate that high-luminosity quasars have a strong effect in ionizing the available gas in a galaxy.

preprint2011arXiv

A spectacular outflow in an obscured quasar

SDSS J1356+1026 is a pair of interacting galaxies at redshift z=0.123 that hosts a luminous obscured quasar in its northern nucleus. Here we present two long-slit Magellan LDSS-3 spectra that reveal a pair of symmetric ~10 kpc-size outflows emerging from this nucleus, with observed expansion velocities of ~250 km/s in projection. We present a kinematic model of these outflows and argue that the deprojected physical velocities of expansion are likely ~1000 km/s and that the kinetic energy of the expanding shells is likely 10^44-10^45 erg/s, with an absolute minimum of >10^42 erg/s. Although a radio counterpart is detected at 1.4GHz, it is faint enough that the quasar is considered to be radio-quiet by all standard criteria, and there is no evidence of extended emission due to radio lobes, whether aged or continuously powered by an ongoing jet. We argue that the likely level of star formation is probably insufficient to power the observed energetic outflow and that SDSS J1356+1026 makes a strong case for radio-quiet quasar feedback. In further support of this hypothesis, polarimetric observations show that the direction of quasar illumination is coincident with the direction of the outflow.

preprint2011arXiv

Feedback in Luminous Obscured Quasars

We use spatially resolved long-slit spectroscopy from Magellan to investigate the extent, kinematics, and ionization structure in the narrow-line regions of 15 luminous, obscured quasars with z<0.5. Increasing the dynamic range in luminosity by an order of magnitude, as well as improving the depth of existing observations by a similar factor, we revisit relations between narrow-line region size and the luminosity and linewidth of the narrow emission lines. We find a slope of 0.22 +/- 0.04 for the power-law relationship between size and luminosity, suggesting that the nebulae are limited by availability of gas to ionize at these luminosities. In fact, we find that the active galactic nucleus is effectively ionizing the interstellar medium over the full extent of the host galaxy. Broad (~300-1000 km/s) linewidths across the galaxies reveal that the gas is kinematically disturbed. Furthermore, the rotation curves and velocity dispersions of the ionized gas remain constant out to large distances, in striking contrast to normal and starburst galaxies. We argue that the gas in the entire host galaxy is significantly disturbed by the central active galactic nucleus. While only ~10^7-10^8 M_sun worth of gas are directly observed to be leaving the host galaxies at or above their escape velocities, these estimates are likely lower limits because of the biases in both mass and outflow velocity measurements and may in fact be in accord with expectations from recent feedback models. Additionally, we report the discovery of two dual obscured quasars, one of which is blowing a large-scale (~10 kpc) bubble of ionized gas into the intergalactic medium.

preprint2010arXiv

H_2 emission arises outside photodissociation regions in ultra-luminous infrared galaxies

Ultra-luminous infrared galaxies are among the most luminous objects in the local universe and are thought to be powered by intense star formation. It has been shown that in these objects the rotational spectral lines of molecular hydrogen observed at mid-infrared wavelengths are not affected by dust obscuration, leaving unresolved the source of excitation of this emission. Here I report an analysis of archival Spitzer Space Telescope data on ultra-luminous infrared galaxies and demonstrate that star formation regions are buried inside optically thick clouds of gas and dust, so that dust obscuration affects star-formation indicators but not molecular hydrogen. I thereby establish that the emission of H_2 is not co-spatial with the buried starburst activity and originates outside the obscured regions. This is rather surprising in light of the standard view that H_2 emission is directly associated with star-formation activity. Instead, I propose that H_2 emission in these objects traces shocks in the surrounding material, which are in turn excited by interactions with nearby galaxies, and that powerful large-scale shocks cooling by means of H_2 emission may be much more common than previously thought. In the early universe, a boost in H_2 emission by this process may speed up the cooling of matter as it collapsed to form the first stars and galaxies and would make these first structures more readily observable.

preprint2010arXiv

Observational biases in determining extrasolar planet eccentricities in single-planet systems

We investigate potential biases in the measurements of exoplanet orbital parameters obtained from radial velocity observations for single-planet systems. We create a mock catalog of radial velocity data, choosing input planet masses, periods, and observing patterns from actual radial velocity surveys and varying input eccentricities. We apply Markov Chain Monte Carlo (MCMC) simulations and compare the resulting orbital parameters to the input values. We find that a combination of the effective signal-to-noise ratio of the data, the maximal gap in phase coverage, and the total number of periods covered by observations is a good predictor of the quality of derived orbit parameters. As eccentricity is positive definite, we find that eccentricities of planets on nearly circular orbits are preferentially overestimated, with typical bias of 1-2 times the median eccentricity uncertainty in a survey (e.g., 0.04 in the Butler et al. 2006 catalog). When performing population analysis, we recommend using the mode of the marginalized posterior eccentricity distribution to minimize potential biases. While the Butler et al. (2006) catalog reports eccentricities below 0.05 for just 17% of single-planet systems, we estimate that the true fraction of e < 0.05 orbits is about f(0.05)=38\pm 9%. For planets with P > 10 days, we find f(0.05)=28\pm 8% versus 10% from Butler et al. (2006). These planets either never acquired a large eccentricity or were circularized following any significant eccentricity excitation.

preprint2010arXiv

Space density of optically-selected type 2 quasars

Type 2 quasars are luminous active galactic nuclei (AGN) whose central regions are obscured by large amounts of gas and dust. In this paper, we present a catalog of type 2 quasars from the Sloan Digital Sky Survey (SDSS), selected based on their optical emission lines. The catalog contains 887 objects with redshifts z < 0.83; this is six times larger than the previous version and is by far the largest sample of type 2 quasars in the literature. We derive the [OIII]5008 luminosity function for 10^8.3 Lsun < L[OIII] < 10^10 Lsun (corresponding to intrinsic luminosities up to M[2400A]-28 mag or bolometric luminosities up to 4x10^47 erg/sec). This luminosity function provides strong lower limits to the actual space density of obscured quasars, due to our selection criteria, the details of the spectroscopic target selection, as well as other effects. We derive the equivalent luminosity function for the complete sample of type 1 (unobscured) quasars; then, we determine the ratio of type 2/type 1 quasar number densities. Our best data constrain this ratio to be at least 1.5:1 for 10^8.3 Lsun < L[OIII] < 10^9.5 Lsun at z < 0.3, and at least 1.2:1 for L[OIII]=10^10 Lsun at 0.3 < z < 0.83. Type 2 quasars are at least as abundant as type 1 quasars in the relatively nearby Universe (z < 0.8) for the highest luminosities.

preprint2009arXiv

Eight-Dimensional Mid-Infrared/Optical Bayesian Quasar Selection

We explore the multidimensional, multiwavelength selection of quasars from mid-IR (MIR) plus optical data, specifically from Spitzer-IRAC and the Sloan Digital Sky Survey (SDSS). We apply modern statistical techniques to combined Spitzer MIR and SDSS optical data, allowing up to 8-D color selection of quasars. Using a Bayesian selection method, we catalog 5546 quasar candidates to an 8.0 um depth of 56 uJy over an area of ~24 sq. deg; ~70% of these candidates are not identified by applying the same Bayesian algorithm to 4-color SDSS optical data alone. Our selection recovers 97.7% of known type 1 quasars in this area and greatly improves the effectiveness of identifying 3.5<z<5 quasars. Even using only the two shortest wavelength IRAC bandpasses, it is possible to use our Bayesian techniques to select quasars with 97% completeness and as little as 10% contamination. This sample has a photometric redshift accuracy of 93.6% (Delta Z +/-0.3), remaining roughly constant when the two reddest MIR bands are excluded. While our methods are designed to find type 1 (unobscured) quasars, as many as 1200 of the objects are type 2 (obscured) quasar candidates. Coupling deep optical imaging data with deep mid-IR data could enable selection of quasars in significant numbers past the peak of the quasar luminosity function (QLF) to at least z~4. Such a sample would constrain the shape of the QLF and enable quasar clustering studies over the largest range of redshift and luminosity to date, yielding significant gains in our understanding of quasars and the evolution of galaxies.

preprint2009arXiv

Host Galaxies of Luminous Type 2 Quasars at z ~ 0.5

We present deep Gemini GMOS optical spectroscopy of nine luminous quasars at redshifts z ~ 0.5, drawn from the SDSS type 2 quasar sample. Our targets were selected to have high intrinsic luminosities (M_V < -26 mag) as indicated by the [O III] 5007 A emission-line luminosity (L_[O III]). Our sample has a median black hole mass of ~ 10^8.8 M_sun inferred assuming the local M_BH-sigma_* relation and a median Eddington ratio of ~ 0.7, using stellar velocity dispersions sigma_* measured from the G band. We estimate the contamination of the stellar continuum from scattered quasar light based on the strength of broad H-beta, and provide an empirical calibration of the contamination as a function of L_[O III]; the scattered light fraction is ~ 30% of L_5100 for objects with L_[O III] = 10^9.5 L_sun. Population synthesis indicates that young post-starburst populations (< 0.1 Gyr) are prevalent in luminous type 2 quasars, in addition to a relatively old population (> 1 Gyr) which dominates the stellar mass. Broad emission complexes around He II 4686 A with luminosities up to 10^8.3 L_sun are unambiguously detected in three out of the nine targets, indicative of Wolf-Rayet populations. Population synthesis shows that ~ 5-Myr post-starburst populations contribute substantially to the luminosities (> 50% of L_5100) of all three objects with Wolf-Rayet detections. We find two objects with double cores and four with close companions. Our results may suggest that luminous type 2 quasars trace an early stage of galaxy interaction, perhaps responsible for both the quasar and the starburst activity.

preprint2008arXiv

Hot self-similar relativistic MHD flows

We consider axisymmetric relativistic jets with a toroidal magnetic field and an ultrarelativistic equation of state, with the goal of studying the lateral structure of jets whose pressure is matched to the pressure of the medium through which they propagate. We find all self-similar steady-state solutions of the relativistic MHD equations for this setup. One of the solutions is the case of a parabolic jet being accelerated by the pressure gradient as it propagates through a medium with pressure declining as p(z)\propto z^{-2}. As the jet material expands due to internal pressure gradients, it runs into the ambient medium resulting in a pile-up of material along the jet boundary, while the magnetic field acts to produce a magnetic pinch along the axis of the jet. Such jets can be in a lateral pressure equilibrium only if their opening angle θ_j at distance z is smaller than about 1/γ, where γis the characteristic bulk Lorentz-factor at this distance; otherwise, different parts of the jet cannot maintain causal contact. We construct maps of optically thin synchrotron emission from our self-similar models. We suggest that the boundary pile-up may be the reason for the limb-brightening of the sub-parsec jet of M87. We find that if the synchrotron emissivity falls with the distance from the jet axis, the polarization fraction rises toward the edge, as seen in 3C273 and Mkn501. Projection effects and the emissivity pattern of the jet have a strong effect on the observed polarization signal, so the interpretation of the polarization data in terms of the geometry of magnetic fields is rather uncertain. For example, jets with toroidal magnetic fields display the `spine-sheath' polarization angle pattern seen in some BL Lac objects.

preprint2008arXiv

Mid-infrared spectra of optically selected type 2 quasars

Type 2 quasars are luminous Active Galactic Nuclei whose central engines are seen through large amounts of gas and dust. We present Spitzer spectra of twelve type 2 quasars selected on the basis of their optical emission line properties. Within this sample, we find a surprising diversity of spectra, from those that are featureless to those showing strong PAH emission, deep silicate absorption at 10 micron, hydrocarbon absorption, high-ionization emission lines and H_2 rotational emission lines. About half of the objects in the sample are likely Compton-thick, including the two with the deepest Si absorption. The median star-formation luminosity of the objects in our sample measured from the strength of the PAH features is 5x10^11 L_sun, much higher than for field galaxies or for any other AGN sample, but similar to other samples of type 2 quasars. This suggests an evolutionary link between obscured quasars and peak star formation activity in the host galaxy. Despite the high level of star formation, the bolometric output is dominated by the quasar in all cases. For a given strength of 10 micron Si absorption, ULIRGs are significantly colder than are type 2 quasars (their F_nu[14.5 micron]/F_nu[27.5 micron] ratio is 0.5 dex lower), perhaps reflecting different obscuration geometries in these sources. We find that the appearance of the 10 micron feature (i.e., whether it shows in emission or in absorption) is well-correlated with the optical classification in type 1 and type 2 quasars, contrary to some models of clumpy obscuration. Furthermore, this correlation is significantly stronger in quasars (L_bol>10^45 erg/s) than it is in Seyfert galaxies (L_bol<<10^45 erg/s).

preprint2006arXiv

Quasars Probing Quasars I: Optically Thick Absorbers Near Luminous Quasars

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We search 149 moderate resolution background quasar spectra, from Gemini, Keck, the MMT, and the SDSS to survey Lyman Limit Systems (LLSs) and Damped Ly-alpha systems (DLAs) in the vicinity of 1.8 < z < 4.0 luminous foreground quasars. A sample of 27 new quasar-absorber pairs is uncovered with column densities, 17.2 < log (N_HI/cm^2) < 20.9, and transverse (proper) distances of 22 kpc/h < R < 1.7 Mpc/h, from the foreground quasars. If they emit isotropically, the implied ionizing photon fluxes are a factor of ~ 5-8000 times larger than the ambient extragalactic UV background over this range of distances. The observed probability of intercepting an absorber is very high for small separations: six out of eight projected sightlines with transverse separations R < 150 kpc/h have an absorber coincident with the foreground quasar, of which four have log N_HI > 10^19. The covering factor of log N_HI > 10^19 absorbers is thus ~ 50 % (4/8) on these small scales, whereas < 2% would have been expected at random. There are many cosmological applications of these new sightlines: they provide laboratories for studying fluorescent Ly-alpha recombination radiation from LLSs, constrain the environments, emission geometry, and radiative histories of quasars, and shed light on the physical nature of LLSs and DLAs.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: I. Sample Selection

We have compiled a large sample of low-redshift active galactic nuclei (AGN) identified via their emission line characteristics from the spectroscopic data of the Sloan Digital Sky Survey. Since emission lines are often contaminated by stellar absorption lines, we developed an objective and efficient method of subtracting the stellar continuum from every galaxy spectrum before making emission line measurements. The distribution of the measured H$α$ Full Width at Half Maxima values of emission line galaxies is strongly bimodal, with two populations separated at about 1,200km s$^{-1}$. This feature provides a natural separation between narrow-line and broad-line AGN. The narrow-line AGN are identified using standard emission line ratio diagnostic diagrams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 galaxy spectra selected over 1151 square degrees. This sample is used in a companion paper to determine the emission-line luminosity function of AGN.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: II. Emission-Line Luminosity Function

The emission line luminosity function of active galactic nuclei (AGN) is measured from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky Survey within a redshift range of $0<z<0.15$. The $\Ha$ and [OIII]$λ5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$α$ and the shapes are well fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming a universal conversion from emission line strength to continuum luminosity, the inferred B band magnitude luminosity function is comparable both to the AGN luminosity function of previous studies and to the low redshift quasar luminosity function derived from the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of all galaxies and about $6\times 10^{-3}$ of the total light of galaxies in the $r$-band comes from the nuclear activity. The numbers of Seyfert 1s and Seyfert 2s are comparable at low luminosity, while at high luminosity, Seyfert 1s outnumber Seyfert 2s by a factor of 2-4. In making the luminosity function measurements, we assumed that the nuclear luminosity is independent of the host galaxy luminosity, an assumption we test {\it a posteriori}, and show to be consistent with the data. Given the relationship between black hole mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host luminosity suggests that the main variable that determines the AGN luminosity is the Eddington ratio, not the black hole mass. This appears to be different from luminous quasars, which are most likely to be shining near the Eddington limit.

preprint2003arXiv

Models of Galaxy Clusters with Thermal Conduction

We present a simple model of hot gas in galaxy clusters, assuming hydrostatic equilibrium and energy balance between radiative cooling and thermal conduction. For five clusters, A1795, A1835, A2199, A2390 and RXJ1347.5-1145, the model gives a good description of the observed radial profiles of electron density and temperature, provided we take the thermal conductivity $κ$ to be about 30% of the Spitzer conductivity. Since the required $κ$ is consistent with the recent theoretical estimate of Narayan & Medvedev (2001) for a turbulent magnetized plasma, we consider a conduction-based equilibrium model to be viable for these clusters. We further show that the hot gas is thermally stable because of the presence of conduction. For five other clusters, A2052, A2597, Hydra A, Ser 159-03 and 3C295, the model requires unphysically large values of $κ$ to fit the data. These clusters must have some additional source of heat, most likely an active galactic nucleus since all the clusters have strong radio galaxies at their centers. We suggest that thermal conduction, though not dominant in these clusters, may nevertheless play a significant role by preventing the gas from becoming thermally unstable.