Source author record

Luis C. Ho

Luis C. Ho 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

210works
6topics
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

210 published item(s)

preprint2026arXiv

Asymmetric torus variability in active galactic nuclei driven by global brightening and dimming

Temporal asymmetry in the flux variability of active galactic nuclei (AGNs) offers key insights into the physical mechanisms driving AGN variability. In this study, we investigated the variability of the torus by analyzing temporal asymmetry in the mid-infrared (MIR) continuum. We compared ensemble structure functions between the brightening and dimming phases for AGNs at $0.15<z<0.4$, using monitoring data in the optical from the Zwicky Transient Facility and in the MIR from the Near-Earth Object Wide-field Infrared Survey Explorer. We found that AGNs with bluer optical-to-MIR colors exhibit positive temporal asymmetry in the MIR, indicating that their variability amplitude is larger when brightening. Conversely, those with redder colors show negative asymmetry, exhibiting larger variability amplitude when decaying. However, there is no significant temporal asymmetry in the $g$-band variability driven by the accretion disk, suggesting that the temporal asymmetry in the MIR continuum primarily originates from intrinsic processes in the torus instead of the reflection of the ultraviolet-optical variability from the accretion disk. Analysis of the composite light curves revealed that AGNs with bluer optical-to-MIR colors tend to brighten gradually in the MIR, leading to the observed temporal asymmetry. This finding suggests that hot-dust-rich AGNs evolve with a gradual decline in hot dust emission, while hot-dust-poor AGNs are associated with a steady increase.

preprint2026arXiv

Composite spectrum of Little Red Dot from a standard inner disk and an unstable outer disk

James Webb Space Telescope (JWST) has revealed a new class of high-redshift, very red, compact broad-line sources, termed as "little red dots" (LRDs). The physical mechanism driving these properties remains elusive. We construct spectral energy distributions (SEDs) with spectroscopic redshift for 28 LRDs and find they exhibit V-shaped SEDs with a common break frequency of $ν_{\rm b}\simeq10^{14.96\pm0.06}$ Hz. We propose that the unique SEDs can be well explained by the combination of an inner standard disk and an outer gravitationally unstable accretion disk with Toomre parameter $Q\sim1$, where the outer disk has a temperature of $\sim2000-4000 K$ and mainly radiates in near-infrared to optical wavebands. The composite spectrum from this model naturally explains the V-shaped continuum and reproduces intrinsically luminous infrared-optical emission without requiring extreme dust extinction or unusual stellar populations. Even considering possible dense gas around the disk to account for pronounced Balmer breaks in some LRDs, the intrinsic optical-UV emission is only suppressed by factors of $\lesssim2-3$, which suggests that most LRDs are sub-Eddington and intrinsically weak. These results provide new insights into early-phase black hole growth and galaxy evolution.

preprint2026arXiv

Galactic bars are already mature at Cosmic Noon: bar strength and flatness at z ~ 1.5

In this work, we explore the nature of $z>1$ galactic bars. Once thought to be highly transient, our results demonstrate otherwise. Our sample consists of nine massive ($>10^{10.5}\,\rm M_{\odot}$) star-forming barred-spiral galaxies at $z_{\rm spec} \sim 1.5$. Using rest-frame near-IR (F444W) JWST/NIRCam imaging, we apply ellipse fitting along with 1D and 2D morphological modeling to directly measure bar properties. We find that five galaxies host flat surface brightness profiles (bar Sérsic index $<0.4$), indicative of highly evolved, "mature" bars. By contrast, only two galaxies show exponential profiles, characteristic of young bars, and these are also shorter in absolute length than the flat bars. We therefore conclude that a large fraction of bars at this epoch have already matured, thereby indicating the presence of well-settled disks required to facilitate bar formation and sustained evolution well before $z\sim1.5$. To assess the gravitational impact of the bars, we calculate the maximum transverse-to-radial force ratio ($Q_{b}$). We find that $Q_{b}$ values are comparable to, or weaker than, those of bars in the local Universe, Seven of the nine bars show only a marginal increase in strength with maturity (from exponential to flat bars). Contrarily however, the remaining two bars are flat, but have the lowest $Q_{b}$ values in our sample. We hence propose that the mature bars at $z\sim 1.5$ may experience phases of weakening due to rapid gas inflows and/or minor mergers. In conclusion, our work sheds light on the rapidly evolving nature of high-z bars and paves the way for larger statistical studies.

preprint2026arXiv

JWST Insights into Narrow-line Little Red Dots

James Webb Space Telescope (JWST) has revealed a population of red and compact objects with a unique V-shape SED at z >= 4 known as Little Red Dots (LRDs). Most of the LRDs with existing spectral observations exhibit broad Balmer lines and are thus likely to host active galactic nuclei (AGNs). Here we present a study of LRDs with no broad H-alpha component. Our sample consists of five LRDs at z~5 with H-alpha line widths of about 250 km/s. They are selected from 32 LRDs that have NIRSpec high- or medium-resolution grating spectra covering H-alpha. During our construction of the sample, we find that approximately 20 percent of the LRD candidates previously selected do not show red continuum emission but resemble the V-shape spectra due to strong line emission. Compared to normal star-forming galaxies, narrow-line LRDs tend to have relatively higher H-alpha line widths and luminosities. If these LRDs are dominated by galaxies, our SED modeling suggests that they are dusty, compact star-forming galaxies with high stellar masses and star formation rates (SFRs). Alternatively, if their SEDs are produced by AGNs, the inferred central black hole masses (MBH) are in the range of 10^5 to 10^6 solar masses, placing them at the low-mass end of the AGN population. They may represent an early stage of super-Eddington growth, where the black holes have yet to accumulate significant masses. With large measurement uncertainties, these black holes appear slightly overmassive relative to the local MBH-Mstar relation, but consistent or undermassive with respect to the MBH-sigma and MBH-Mdyn relations. We further find that nearly half of the high-redshift broad-line AGNs exhibit V-shape SEDs. (abridged)

preprint2026arXiv

The Dichotomy in the Nuclear and Host Galaxy Properties of High-redshift Quasars

The early growth of high-redshift quasars and their host galaxies raises critical questions about their cosmic evolution. We exploit the angular resolution and sensitivity of NIRCam to investigate the host galaxies of 31 quasars at $4\lesssim z\lesssim7$ drawn from multiple JWST surveys. Using a new multi-band forward-modeling code (\textsc{GalfitS}) that incorporates physically motivated priors, we securely detect and quantify the host emission in 30 objects, while simultaneously characterizing the nuclear spectral energy distribution. The host galaxies of high-redshift quasars are $\sim 0.3$~dex more compact than star-forming galaxies of comparable mass. A striking dichotomy emerges: luminous ``blue'' quasars ($L_{5100}\gtrsim10^{45}\,{\rm erg\,s^{-1}}$) reside in bulge-dominated galaxies ($n \approx 5$) and exhibit a narrow range of ultraviolet nuclear slopes (median $β_{\rm UV} \approx -1.4$), while fainter ``red'' quasars inhabit disk-like hosts ($n\approx 1$) and display a broad range of slopes ($β_{\rm UV}\approx-2$ to 4). These two populations differ markedly in their black hole-to-stellar mass ratios, with high-luminosity quasars showing $M_{\mathrm{BH}}/M_\ast = 1.2\%$ compared to $4.7\%$ for lower luminosity sources, placing them collectively $\sim$0.6~dex above the local $M_{\mathrm{BH}}-M_\ast$ relation. This offset likely reflects rapid black hole growth in early gas-rich environments, where feedback from the active galactic nucleus becomes effective only after substantial gas depletion. Our findings suggest that the observed dichotomy, whether due to intrinsic spectral differences or dust extinction, fundamentally shapes the coevolution of supermassive black holes and their host galaxies in the early Universe.

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2023arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: The Black Hole Mass$-$Stellar Mass Relations at $0.2\lesssim z\lesssim 0.8$

We measure the correlation between black-hole mass $M_{\rm BH}$ and host stellar mass $M_*$ for a sample of 38 broad-line quasars at $0.2\lesssim z\lesssim 0.8$ (median redshift $z_{\rm med}=0.5$). The black-hole masses are derived from a dedicated reverberation mapping program for distant quasars, and the stellar masses are estimated from two-band optical+IR HST imaging. Most of these quasars are well centered within $\lesssim 1$kpc from the host galaxy centroid, with only a few cases in merging/disturbed systems showing larger spatial offsets. Our sample spans two orders of magnitude in stellar mass ($\sim 10^9-10^{11}\,M_\odot$) and black-hole mass ($\sim 10^7-10^9\,M_\odot$), and reveals a significant correlation between the two quantities. We find a best-fit intrinsic (i.e., selection effects corrected) $M_{\rm BH}-M_{\rm *,host}$ relation of $\log (M_{\rm BH}/M_{\rm \odot})=7.01_{-0.33}^{+0.23} + 1.74_{-0.64}^{+0.64}\log (M_{\rm *,host}/10^{10}M_{\rm \odot})$, with an intrinsic scatter of $0.47_{-0.17}^{+0.24}$dex. Decomposing our quasar hosts into bulges and disks, there is a similar $M_{\rm BH}-M_{\rm *,bulge}$ relation with a slightly larger scatter, likely caused by systematic uncertainties in the bulge-disk decomposition. The $M_{\rm BH}-M_{\rm *,host}$ relation at $z_{\rm med}=0.5$ is similar to that in local quiescent galaxies, with negligible evolution over the redshift range probed by our sample. With direct black-hole masses from reverberation mapping and a large dynamical range of the sample, selection biases do not appear to affect our conclusions significantly. Our results, along with other samples in the literature, suggest that the locally-measured black-hole mass$-$host stellar mass relation is already in place at $z\sim 1$.

preprint2022arXiv

A Candidate of a Least-Massive Black Hole at the First 1.1 Billion Years of the Universe

We report a candidate of a low-luminosity active galactic nucleus (AGN) at z=5 that was selected from the first near-infrared images of the JWST CEERS project. This source, named CEERS-AGN-z5-1 at absolute 1450 A magnitude M1450=-19.5 +/- 0.3, was found via a visual selection of compact sources from a catalog of Lyman break galaxies at z>4, taking advantage of the superb spatial resolution of the JWST/NIRCam images. The 20 photometric data available from CFHT, HST, Spitzer, and JWST suggest that the continuum shape of this source is reminiscent of that for an unobscured AGN, and there is a clear color excess in the filters where the redshifted Hbeta+[OIII] and Halpha are covered. The estimated line luminosity is L_Hbeta+[OIII] =10^43.0 erg s-1 and L_Halpha =10^42.9 erg s-1 with the corresponding rest-frame equivalent width EW_{Hbeta+[OIII]} =1100 A and EW_Halpha =1600 A, respectively. Our SED fitting analysis favors the scenario that this object is either a strong broad-line emitter or even a super-Eddington accreting black hole (BH), although a possibility of an extremely young galaxy with moderate dust attenuation is not completely ruled out. The bolometric luminosity, L_bol=2.5 +/- 0.3 \times 10^44 erg s-1, is consistent with those of z<0.35 broad-line AGNs with M_BH = 10^6 M_sun accreting at the Eddington limit. This new AGN population at the first 1.1 billion years of the universe may close the gap between the observed BH mass range at high redshift and that of BH seeds. Spectroscopic confirmation is awaited to secure the redshift and its AGN nature.

preprint2022arXiv

A quasar shedding its dust cocoon at redshift 2

We present the first near-IR spectroscopy and joint analyses of multi-wavelength observations for SDSS J082747.14+425241.1, a dust-reddened, weak broad emission-line quasar (WLQ) undergoing a remarkable broad absorption line (BAL) transformation. The systemic redshift is more precisely measured to be $z=2.070\pm0.001$ using H$β$ compared to $z=2.040\pm0.003$ using \mgii\ from the literature, signifying an extreme \mgii\ blueshift of $2140\pm530$ \kms\ relative to H$β$. Using the H$β$-based single-epoch scaling relation with a systematic uncertainty of 0.3 dex, its black hole (BH) mass and Eddington ratio are estimated to be $M_{\rm BH}\sim6.1\times10^8M_\odot$ and $λ_{\rm Edd}\sim0.71$, indicative of being in a rapidly accreting phase. Our investigations confirm the WLQ nature and the LoBAL$\rightarrow$HiBAL transformation, along with a factor of 2 increase in the \mgii+\feii\ emission strength and a decrease of 0.1 in $E(B-V)$ over two decades. The kinetic power of this LoBAL wind at $R\sim$15 pc from its BH is estimated to be $\sim$43\% of the Eddington luminosity, sufficient for quasar feedback upon its host galaxy albeit with an order-of-magnitude uncertainty. This quasar provides a clear example of the long-sought scenario where LoBAL quasars are surrounded by dust cocoons, and wide-angle nuclear winds play a key role in the transition for red quasars evolving into the commonly seen blue quasars.

preprint2022arXiv

Accretion Disk Size Measurements of Active Galactic Nuclei Monitored by the Zwicky Transient Facility

We compile a sample of 92 active galactic nuclei (AGNs) at z<0.75 with $gri$ photometric light curves from the archival data of the Zwicky Transient Facility and measure the accretion disk sizes via continuum reverberation mapping. We employ Monte Carlo simulation tests to assess the influences of data sampling and broad emission lines and select out the sample with adequately high sampling cadences (3 days apart in average) and minimum contaminations of broad emission lines. The inter-band time delays of individual AGNs are calculated using the interpolated cross-correlation function and then these delays are fitted with a generalized accretion disk model, in which inter-band time delays are a power function of wavelength, black hole mass, and luminosity. A Markov-chain Monte Carlo method is adopted to determine the best parameter values. Overall the inter-band time delays can be fitted with the $τ\ \propto λ^{4/3}$ relation as predicted from a steady-state, optically thick, geometrically thin accretion disk, however, the yielded disk size is systematically larger than expected, although the ratio of the measured to theoretical disk sizes depend on using the emissivity -- or responsivity -- weighted disk radius. These results are broadly consistent with previous studies, all together raising a puzzle about the "standard" accretion disk model.

preprint2022arXiv

Accretion Disk Sizes from Continuum Reverberation Mapping of AGN Selected from the ZTF Survey

We present the accretion disk size estimates for a sample of 19 active galactic nuclei (AGN) using the optical $g$, $r$, and $i$ band light curves obtained from the Zwicky Transient Facility (ZTF) survey. All the AGN have reliable supermassive black hole (SMBH) mass estimates based on previous reverberation mapping measurements. The multi-band light curves are cross-correlated, and the reverberation lag is estimated using the Interpolated Cross-Correlation Function (ICCF) method and the Bayesian method using the {\sc javelin} code. As expected from the disk reprocessing arguments, the $g-r$ band lags are shorter than the $g-i$ band lags for this sample. The interband lags for all, but 5 sources, are larger than the sizes predicted from the standard Shakura Sunyaev (SS) analytical model. We fit the light curves directly using a thin disk model implemented through the {\sc javelin} code to get the accretion disk sizes. The disk sizes obtained using this model are on an average 3.9 times larger than the prediction based on the SS disk model. We find a weak correlation between the disk sizes and the known physical parameters, namely, the luminosity and the SMBH mass. In the near future, a large sample of AGN covering broader ranges of luminosity and SMBH mass from large photometric surveys would be helpful in a better understanding of the structure and physics of the accretion disk.

preprint2022arXiv

An Elusive Population of Massive Disk Galaxies Hosting Double-lobed Radio-loud AGNs

It is commonly accepted that radio-loud active galactic nuclei are hosted exclusively by giant elliptical galaxies. We analyze high-resolution optical Hubble Space Telescope images of a sample of radio galaxies with extended double-lobed structures associated with disk-like optical counterparts. After systematically evaluating the probability of chance alignment between the radio lobes and the optical counterparts, we obtain a sample of 18 objects likely to have genuine associations. The host galaxies have unambiguous late-type morphologies, including spiral arms, large-scale dust lanes among the edge-on systems, and exceptionally weak bulges, as judged by the low global concentrations, small global Sérsic indices, and low bulge-to-total light ratios (median $B/T = 0.13$). With a median Sérsic index of 1.4 and low effective surface brightnesses, the bulges are consistent with being pseudo bulges. The majority of the hosts have unusually large stellar masses (median $M_* = 1.3\times 10^{11}\, M_\odot$) and red optical colors (median $g-r = 0.69\,$mag), consistent with massive, quiescent galaxies on the red sequence. We suggest that black hole mass (stellar mass) plays a fundamental role in launching large-scale radio jets, and that the rarity of extended radio lobes in late-type galaxies is the consequence of the steep stellar mass function at the high-mass end. The disk radio galaxies have mostly Fanaroff-Riley type II morphologies yet lower radio power than sources of a similar type traditionally hosted by ellipticals. The radio jets show no preferential alignment with the minor axis of the galactic bulge or disk, apart from a possible mild tendency for alignment among the most disk-dominated systems.

preprint2022arXiv

Asymmetric star formation triggered by gas inflow in a barred lenticular galaxy PGC 34107

Comparing to the inactive and gas-poor normal lenticular galaxies (S0s) in the local universe, we study a barred star-forming S0 galaxy, PGC 34107, which has been observed by the Centro Astronómico Hispano Alemán (CAHA) 3.5-m telescope and the Northern Extended Millimeter Array (NOEMA). The spatially resolved ionized gas and molecular gas traced by $^{12}$CO(1-0), hereafter CO(1-0), show the similar distribution and kinematics to the stellar component with an off-center star-forming region, $\sim$380 pc away from the center. The resolved kinematics of molecular CO(1-0) emission reveals that there is a blueshifted (redshifted) velocity component on the receding (approaching) side of the galaxy along the stellar bar. This might provide a plausible evidence of non-circular motion, such as the bar-induced molecular gas inflow. The velocity of molecular gas inflow decreases with approaching towards the peak of the off-center star formation in the north, which might be associated with the inner Lindblad resonance (ILR). In addition to CO(1-0), we also detect the isotopic line of $^{13}$CO(1-0). Most $\rm Hα$, CO(1-0) and $^{13}$CO(1-0) emissions are concentrated on this northern star-forming region. We find that PGC 34107 follows the local stellar mass-metallicity relation, star-forming main sequence, and the Kennicutt-Schmidt law. The resolved and integrated molecular gas main sequence suggest that there is a higher gas fraction in the galaxy central region, which supports a scenario that the bar-induced gas reservoir provides the raw material, and subsequently triggers the central star formation.

preprint2022arXiv

Centrally Concentrated H I Distribution Enhances Star Formation in Galaxies

We use a sample of 13,511 nearby galaxies from the ALFALFA and SDSS spectroscopic surveys to study the relation between the spatial distribution of H I 21 cm emission and star formation rate (SFR). We introduce a new non-parametric quantity $K$, measured from the curve-of-growth of the line, to describe the shape of the integrated H I profile. The value of $K$ increases from double-horned to single-peaked profiles, depending on projection effects and the spatial and velocity distribution of the gas. Using carefully chosen samples to control for the competing factors that influence the integrated line profile, we argue that useful inferences can be made on the spatial distribution of the gas. We find that galaxies with a high value of $K$ tend to have more centrally concentrated H I distribution within the optical disk of the galaxy at fixed conditions, and that larger values of $K$ are associated with higher levels of total and central SFR. The results suggest that the global concentration of H I plays an important role in facilitating the conversion of neutral atomic hydrogen to molecular hydrogen gas, which, in turn, affects the star formation activity throughout the optical disk. Our sample is biased against quiescent galaxies, and thus the conclusions may not hold for galaxies with low SFR or low H I content.

preprint2022arXiv

Cold Gas in Massive Galaxies as A Critical Test of Black Hole Feedback Models

Black hole feedback has been widely implemented as the key recipe to quench star formation in massive galaxies in modern semi-analytic models and hydrodynamical simulations. As the theoretical details surrounding the accretion and feedback of black holes continue to be refined, various feedback models have been implemented across simulations, with notable differences in their outcomes. Yet, most of these simulations have successfully reproduced some observations, such as stellar mass function and star formation rate density in the local Universe. We use the recent observation on the change of neutral hydrogen gas mass (including both ${\rm H_2}$ and ${\rm HI}$) with star formation rate of massive central disc galaxies as a critical constraint of black hole feedback models across several simulations. We find that the predictions of IllustrisTNG agree with the observations much better than the other models tested in this work. This favors IllustrisTNG's treatment of active galactic nuclei - where kinetic winds are driven by black holes at low accretion rates - as more plausible amongst those we test. In turn, this also indirectly supports the idea that the massive central disc galaxy population in the local Universe was likely quenched by AGN feedback.

preprint2022arXiv

Definitive upper bound on the negligible contribution of quasars to cosmic reionization

Cosmic (hydrogen) reionization marks one of the major phase transitions of the universe at redshift z >= 6. During this epoch, hydrogen atoms in the intergalactic medium (IGM) were ionized by Lyman continuum (LyC) photons. However, it remains challenging to identify the major sources of the LyC photons responsible for reionization. In particular, individual contributions of quasars (or active galactic nuclei, AGNs) and galaxies are still under debate. Here we construct the far-ultraviolet (far-UV) luminosity function for type 1 quasars at z >= 6 that spans 10 magnitudes (-19 < M_UV < -29), conclusively showing that quasars made a negligible contribution to reionization. We mainly search for quasars in the low-luminosity range of M_UV > -23 mag that is critical to determine quasars' total LyC photon production but has been barely explored previously. We find that the quasar population can only provide less than 7% (95% confidence level) of the total photons needed to keep the universe ionized at z = 6.0 - 6.6. Our result suggests that galaxies, presumably low-luminosity star-forming systems, are the major sources of hydrogen reionization.

preprint2022arXiv

Evidence for quasar fast outflows being accelerated at the scale of tens of parsecs

Quasar outflows may play a crucial role in regulating the host galaxy, although the spatial scale of quasar outflows remain a major enigma, with their acceleration mechanism poorly understood. The kinematic information of outflow is the key to understanding its origin and acceleration mechanism. Here, we report the galactocentric distances of different outflow components for both a sample and an individual quasar. We find that the outflow distance increases with velocity, with a typical value from several parsecs to more than one hundred parsecs, providing direct evidence for an acceleration happening at a scale of the order of 10 parsecs. These outflows carry ~1% of the total quasar energy, while their kinematics are consistent with a dust driven model with a launching radius comparable to the scale of a dusty torus, indicating that the coupling between dust and quasar radiation may produce powerful feedback that is crucial to galaxy evolution.

preprint2022arXiv

Ionized Outflows in Nearby Quasars are Poorly Coupled to their Host Galaxies

We analyze Multi Unit Spectroscopic Explorer observations of nine low-redshift (z < 0.1) Palomar-Green quasar host galaxies to investigate the spatial distribution and kinematics of the warm, ionized interstellar medium, with the goal of searching for and constraining the efficiency of active galactic nucleus (AGN) feedback. After separating the bright AGN from the starlight and nebular emission, we use pixel-wise, kpc-scale diagnostics to determine the underlying excitation mechanism of the line emission, and we measure the kinematics of the narrow-line region (NLR) to estimate the physical properties of the ionized outflows. The radial size of the NLR correlates with the AGN luminosity, reaching scales of $\sim 5\,$kpc and beyond. The geometry of the NLR is well-represented by a projected biconical structure, suggesting that the AGN radiation preferably escapes through the ionization cone. We find enhanced velocity dispersions ($\sim 100\,$km$\,$s$^{-1}$) traced by the H$α$ emission line in localized zones within the ionization cones. Interpreting these kinematic features as signatures of interaction between an AGN-driven ionized gas outflow and the host galaxy interstellar medium, we derive mass outflow rates of $\sim 0.008-1.6\, M_\odot \,$yr$^{-1}$ and kinetic injection rates of $\sim 10^{39}-10^{42} \,$erg$\,$s$^{-1}$, which yield extremely low coupling efficiencies of $\lesssim 10^{-3}$. These findings add to the growing body of recent observational evidence that AGN feedback is highly ineffective in the host galaxies of nearby AGNs.

preprint2022arXiv

Is there a sub-parsec-scale jet base in the nearby dwarf galaxy NGC 4395?

NGC 4395 is a dwarf galaxy at a distance of about 4.3 Mpc (scale: ~0.021 pc mas$^{-1}$). It hosts an intermediate-mass black hole (IMBH) with a mass between ~10$^4$ and ~10$^5$ solar masses. The early radio observations of NGC 4395 with the very long baseline interferometry (VLBI) network, High Sensitivity Array (HSA), at 1.4 GHz in 2005 showed that its nucleus has a sub-mJy outflow-like feature (E) extending over 15 mas. To probe the possibility of the feature E as a continuous jet with a base physically coupled with the accretion disc, we performed deep VLBI observations with the European VLBI Network (EVN) at 5 GHz, and analysed the archival data obtained with the HSA at 1.4 GHz in 2008, NSF's Karl G. Jansky Very Large Array (VLA) at 12-18 GHz and the Atacama Large Millimetre/submillimetre Array (ALMA) at 237 GHz. The feature E displays more diffuse structure in the HSA image of 2008 and has no compact substructure detected in the EVN image. Together with the optically thin steep spectrum and the extremely large angular offset (about 220 mas) from the accurate optical Gaia position, we explain the feature E as nuclear shocks likely formed by the IMBH's episodic ejection or wide-angle outflow. The VLA and ALMA observations find a sub-mJy pc-scale diffuse feature, possibly tracing a thermal free-free emission region near the IMBH. There is no detection of a jet base at the IMBH position in the VLBI maps. The non-detections give an extremely low luminosity of <=4.7 x 10$^{33}$ erg s$^{-1}$ at 5 GHz and indicate no evidence of a disc-jet coupling on sub-pc scales.

preprint2022arXiv

Massive Galaxy Mergers Have Distinctive Global HI Profiles

The global 21 cm HI emission-line profile of a galaxy encodes valuable information on the spatial distribution and kinematics of the neutral atomic gas. Galaxy interactions significantly influence the HI disk and imprint observable features on the integrated HI line profile. In this work, we study the neutral atomic gas properties of galaxy mergers selected from the Great Observatories All-sky LIRG Survey. The HI spectra come from new observations with the Five-hundred-meter Aperture Spherical Telescope and from a collection of archival data. We quantify the HI profile of the mergers with a newly developed method that uses the curve-of-growth of the line profile. Using a control sample of non-merger galaxies carefully selected to match the stellar mass of the merger sample, we show that mergers have a larger proportion of single-peaked HI profiles, as well as a greater tendency for the HI central velocity to deviate from the systemic optical velocity of the galaxy. By contrast, the HI profiles of mergers are not significantly more asymmetric than those of non-mergers.

preprint2022arXiv

Mid-infrared Variability of Low-redshift Active Galactic Nuclei : Constraints on a Hot Dust Component with a Variable Covering Factor

We utilize mid-infrared multi-epoch data from Wide-field Infrared Survey Explorer (WISE) over a $\sim10$-year period in W1 (3.4$μ$m) and W2 (4.6$μ$m) bands to investigate the structure of dusty torus in low-redshift ($0.15<z\le 0.4$) active galactic nuclei (AGN). We calculate a Spearman correlation coefficient ($r_{12}$) between W1 magnitude and W1$-$W2 color based on the light curve in individual objects. Interestingly, $r_{12}$ spans a broad range from $-1$ to $1$, and is detected to be correlated with mean W1$-$W2 color and AGN bolometric luminosity, in the sense that objects with blue W1$-$W2 color and low AGN luminosity tend to become redder (bluer) with increasing (decreasing) W1 brightness in the light curve (i.e. $r_{12} < 0$), although the correlation of $r_{12}$ with the bolometric luminosity is relatively weak. The fit for the spectral energy distribution reveals a significant contribution from the host galaxy in the W1 and W2 bands. However, the dependencies of $r_{12}$ on W1$-$W2 color and AGN luminosity still persist even after careful elimination of the host light contribution. We propose that this result can be explained if covering factor of the hot dust component decreases as AGN luminosity increases.

preprint2022arXiv

Monitoring AGNs with H$β$ Asymmetry. III. Long-term Reverberation Mapping Results of 15 Palomar-Green Quasars

In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire light curves and the measurements from individual observing seasons are provided. Reverberation mapping of 9 of our targets has been attempted for the first time, while the results for 6 others can be compared with previous campaigns. We measure the H$β$ time lags over periods of years and estimate their black hole masses. The long duration of the campaign enables us to investigate their broad line region (BLR) geometry and kinematics for different years by using velocity-resolved lags, which demonstrate signatures of diverse BLR geometry and kinematics. The BLR geometry and kinematics of individual objects are discussed. In this sample, the BLR kinematics of Keplerian/virialized motion and inflow is more common than outflow.

preprint2022arXiv

Relation between Black Hole Mass and Bulge Luminosity in Hard X-ray selected Type 1 AGNs

Using $I$-band images of 35 nearby ($z<0.1$) type 1 active galactic nuclei (AGNs) obtained with Hubble Space Telescope, selected from the 70-month Swift-BAT X-ray source catalog, we investigate the photometric properties of the host galaxies. With a careful treatment of the point-spread function (PSF) model and imaging decomposition, we robustly measure the $I$-band brightness and the effective radius of bulges in our sample. Along with black hole (BH) mass estimates from single-epoch spectroscopic data, we present the relation between BH mass and $I$-band bulge luminosity ($M_{\rm BH}-M_{I,\rm bul}$ relation) of our sample AGNs. We find that our sample lies offset from the $M_{\rm BH}-M_{I,\rm bul}$ relation of inactive galaxies by 0.4 dex, i.e., at a given bulge luminosity, the BH mass of our sample is systematically smaller than that of inactive galaxies. We also demonstrate that the zero point offset in the $M_{\rm BH}-M_{I,\rm bul}$ relation with respect to inactive galaxies is correlated with the Eddington ratio. Based on the Kormendy relation, we find that the mean surface brightness of ellipticals and classical bulges in our sample is comparable to that of normal galaxies, revealing that bulge brightness is not enhanced in our sample. As a result, we conclude that the deviation in the $M_{\rm BH}-M_{I,\rm bul}$ relation from inactive galaxies is possibly because the scaling factor in the virial BH mass estimator depends on the Eddington ratio.

preprint2022arXiv

Signature of Supersonic Turbulence in Galaxy Clusters Revealed by AGN-driven H$α$ Filaments

The hot intracluster medium (ICM) is thought to be quiescent with low observed velocity dispersions. Surface brightness fluctuations of the ICM also suggest that its turbulence is subsonic with a Kolmogorov scaling relation, indicating that the viscosity is suppressed and the kinetic energy cascades to small scales unscathed. However, recent observations of the cold gas filaments in galaxy clusters find that the scaling relations are steeper than that of the hot plasma, signaling kinetic energy losses and the presence of supersonic flows. In this work we use high-resolution simulations to explore the turbulent velocity structure of the cold filaments at the cores of galaxy clusters. Our results indicate that supersonic turbulent structures can be "frozen" in the cold gas that cools and fragments out of a fast, $10^7$ K outflow driven by the central active galactic nucleus (AGN), when the radiative cooling time is shorter than the dynamical sound-crossing time. After the cold gas formation, however, the slope of the velocity structure function (VSF) flattens significantly over short, 10 Myr timescales. The lack of flattened VSF in observations of H$α$ filaments indicates that the H$α$-emitting phase is short-lived for the cold gas in galaxy clusters. On the other hand, the ubiquity of supersonic turbulence revealed by observed filaments strongly suggests that supersonic outflows are an integral part of AGN-ICM interaction, and that AGN activity plays a crucial role at driving turbulence in galaxy clusters.

preprint2022arXiv

Statistical Analysis of H I Profile Asymmetry and Shape for Nearby Galaxies

We present a uniform analysis of the integrated profile of the H I emission line of 29,958 nearby ($z < 0.06$) galaxies extracted from the ALFALFA 21 cm survey. We apply the curve-of-growth technique to derive a database of spectral parameters and robust estimates of their associated uncertainties. Besides the central velocity and total flux, the main catalog provides new measures of line width, profile asymmetry, and profile shape. For a subsample of 13,511 galaxies with optical properties available from the Sloan Digital Sky Survey, we compute inclination angle-corrected line widths, rotation velocities empirically calibrated from spatially resolved observations, and dynamical masses based on H I sizes estimated from the H I mass. To facilitate subsequent scientific applications of the database, we also compile a number of ancillary physical properties of the galaxies, including their optical morphology, stellar mass, and various diagnostics of star formation activity. We use the homogeneous catalog of H I parameters to examine the statistical properties of profile asymmetry and shape. Across the full sample, which covers a wide range of stellar masses and environments, statistically significant H I profile asymmetry is detected in $\sim 20\%$ of the galaxy population. The global H I profiles are $35.2 \pm 0.3\%$ single-peaked, $26.9 \pm 0.3\%$ flat-topped, and $37.9 \pm 0.3\%$ double-horned. At a given inclination angle, double-horned profiles are preferentially associated with galaxies of higher stellar mass or optical concentration, while galaxies of lower mass or concentration tend to have single-peaked profiles.

preprint2022arXiv

Strong spiral arms drive secular growth of pseudo bulges in disk galaxies

Spiral-driven instabilities may drive gas inflow to enhance central star formation in disk galaxies. We investigate this hypothesis using the Sloan Digital Sky Survey (SDSS) in a sample of 2779 nearby unbarred star-forming main-sequence spiral galaxies. The strength of spiral arms is quantified by their average Fourier amplitude relative to the axisymmetric disk. The star formation properties in the central 1--3\,kpc region were derived from the SDSS spectra. We show that galaxies with stronger spiral arms not only tend to have more intense central specific star formation rate (sSFR), larger Balmer absorption line index, and lower 4000-Å break strength, but also have enhanced central sSFR relative to sSFR measured for the whole galaxy. This link is independent of redshift, stellar mass, surface density, and concentration. There is a lack of evidence for strong spiral arms being associated with a significant fraction of starburst or post-starburst galaxies, implying that the spiral-induced central star formation is likely continuous rather than bursty. We also show that stronger spiral arms tend to have an increasing fraction of pseudo bulges, a relatively unchanged fraction of star-forming classical bulges, and a decreasing fraction of quenched classical bulges. Moreover, the concentration of galaxies hosting pseudo bulges mildly increases with stronger spiral arms, implying that spirals help pseudo bulges grow. The connection between spirals and the bulge type is partly attributed to the suppression of spirals by classical bulges and partly to the enhanced central star formation driven by spirals. We explain our results in a picture where spiral arms transport cold gas inward to trigger continuous central star formation, which facilitates the build-up of pseudo bulges. Spiral arms thus play a role in the secular evolution of disk galaxies.

preprint2022arXiv

The Age of Discovery with the James Webb: Excavating the Spectral Signatures of the First Massive Black Holes

The James Webb Space Telescope (JWST) will open a new window of the most distant universe and unveil the early growth of supermassive black holes (BHs) in the first galaxies. In preparation for deep JWST imaging surveys, it is crucial to understand the color selection of high-redshift accreting seed BHs. We model the spectral energy distribution of super-Eddington accreting BHs with millions of solar masses in metal-poor galaxies at $z\gtrsim 8$, applying post-process line transfer calculations to radiation hydrodynamical simulation results. Ten kilosecond exposures with the NIRCam and MIRI broad-band filters are sufficient to detect the radiation flux from the seed BHs with bolometric luminosities of $L_{\rm bol}\simeq 10^{45}~{\rm erg~s}^{-1}$. While the continuum colors are similar to those of typical low-$z$ quasars, strong H$α$ line emission with a rest-frame equivalent width ${\rm EW}_{\rm rest}\simeq 1300~Å$ is so prominent that the line flux affects the broad-band colors significantly. The unique colors, for instance F356W$-$F560W $\gtrsim 1$ at $7<z<8$ and F444W$-$F770W $\gtrsim 1$ at $9<z<12$, provide robust criteria for photometric selection of the rapidly growing seed BHs. Moreover, NIRSpec observations of low-ionization emission lines can test whether the BH is fed via a dense accretion disk at super-Eddington rates.

preprint2022arXiv

The buildup of strongly barred galaxies in the TNG100 simulation

We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stellar mass, in agreement with observational results. Strongly barred galaxies exhibit, overall, lower gas-to-stellar mass ratios compared to unbarred galaxies. The majority of barred galaxies are quenched (sSFR $\sim10^{-11.7} $yr$^{-1}$), whereas unbarred galaxies continue to be active (sSFR $\sim10^{-10.3}$yr$^{-1}$) on the main sequence of star-forming galaxies. We explore the evolution of strongly barred and unbarred galaxies to investigate their formation and quenching histories. We find that strong bars form between $0.5< z< 1.5$, with more massive galaxies hosting older bars. Strong bars form in galaxies with an early-established prominent disc component, undergoing periods of enhanced star formation and black hole accretion, possibly assisted by cosmological inflows. Unbarred galaxies, on the other hand, assemble most of their mass and disc component at late times. The nuclear region of strongly barred galaxies quenches shortly after bar formation, while unbarred galaxies remain active across time. Our findings are indicative of bar quenching, possibly assisted by nuclear feedback processes. We conclude that the cosmological environment, together with small scale feedback processes, determine the chances of a galaxy to form a bar and to rapidly quench its central region.

preprint2022arXiv

The Host Galaxy and Rapidly Evolving Broad-line Region in the Changing-look Active Galactic Nucleus 1ES 1927+654

Changing-look active galactic nuclei (AGNs) present an important laboratory to understand the origin and physical properties of the broad-line region (BLR). We investigate follow-up optical spectroscopy spanning $\sim 500$ days after the outburst of the changing-look AGN 1ES\,1927+654. The emission lines displayed dramatic, systematic variations in intensity, velocity width, velocity shift, and symmetry. Analysis of optical spectra and multi-band images indicate that the host galaxy contains a pseudobulge and a total stellar mass of $3.56_{-0.35}^{+0.38} \times 10^{9}\, M_\odot$. Enhanced continuum radiation from the outburst produced an accretion disk wind, which condensed into BLR clouds in the region above and below the temporary eccentric disk. Broad Balmer lines emerged $\sim 100$ days after the outburst, together with an unexpected, additional component of narrow-line emission. The newly formed BLR clouds then traveled along a similar eccentric orbit ($e \approx 0.6$). The Balmer decrement of the BLR increased by a factor of $\sim 4-5$ as a result of secular changes in cloud density. The drop in density at late times allowed the production of \hei\ and \heii\ emission. The mass of the black hole cannot be derived from the broad emission lines because the BLR is not virialized. Instead, we use the stellar properties of the host galaxy to estimate $M_\mathrm{BH} = 1.38_{-0.66}^{+1.25} \times 10^{6}\, M_\odot$. The nucleus reached near or above its Eddington limit during the peak of the outburst. We discuss the nature of the changing-look AGN 1ES\,1927+654 in the context of other tidal disruption events.

preprint2022arXiv

The initial conditions for young massive cluster formation in the Galactic Centre: convergence of large-scale gas flows

Young massive clusters (YMCs) are compact ($\lesssim$1 pc), high-mass (>10${}^4$ M${}_{\odot}$) stellar systems of significant scientific interest. Due to their rarity and rapid formation, we have very few examples of YMC progenitor gas clouds before star formation has begun. As a result, the initial conditions required for YMC formation are uncertain. We present high-resolution (0.13$^{\prime\prime}$, $\sim$1000 au) ALMA observations and Mopra single-dish data, showing that Galactic Centre dust ridge `Cloud d' (G0.412$+$0.052, mass$\sim 7.6 \times 10^4$ M$_{\odot}$, radius$\sim 3.2$ pc) has the potential to become an Arches-like YMC (10$^4$ M$_{\odot}$, r$\sim$1 pc), but is not yet forming stars. This would mean it is the youngest known pre-star forming massive cluster and therefore could be an ideal laboratory for studying the initial conditions of YMC formation. We find 96 sources in the dust continuum, with masses $\lesssim$3 M$_{\odot}$ and radii of $\sim$10${}^3$ au. The source masses and separations are more consistent with thermal rather than turbulent fragmentation. It is not possible to unambiguously determine the dynamical state of most of the sources, as the uncertainty on virial parameter estimates is large. We find evidence for large-scale ($\sim$1 pc) converging gas flows, which could cause the cloud to grow rapidly, gaining 10$^4$ M$_{\odot}$ within 10$^5$ yr. The highest density gas is found at the convergent point of the large-scale flows. We expect this cloud to form many high-mass stars, but find no high-mass starless cores. If the sources represent the initial conditions for star formation, the resulting IMF will be bottom-heavy.

preprint2022arXiv

The Paschen Jump as a Diagnostic of the Diffuse Nebular Continuum Emission in Active Galactic Nuclei

Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 Å is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope STIS spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multi-component spectral models over the range $6800-9700$ Å and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from $\sim$10% to 50% at 8000 Å, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe II emission, high-order Balmer line emission, and other components.

preprint2022arXiv

The Photon Ring in M87*

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

preprint2022arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: UV-Optical Accretion Disk Measurements with Hubble Space Telescope

We present accretion-disk structure measurements from UV-optical reverberation mapping observations of a sample of eight quasars at 0.24<z<0.85. Ultraviolet photometry comes from two cycles of Hubble Space Telescope monitoring, accompanied by multi-band optical monitoring by the Las Cumbres Observatory network and Liverpool Telescopes. The targets were selected from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project sample with reliable black-hole mass measurements from Hbeta reverberation mapping results. We measure significant lags between the UV and various optical griz bands using JAVELIN and CREAM methods. We use the significant lag results from both methods to fit the accretion-disk structure using a Markov chain Monte Carlo approach. We study the accretion disk as a function of disk normalization, temperature scaling, and efficiency. We find direct evidence for diffuse nebular emission from Balmer and FeII lines over discrete wavelength ranges. We also find that our best-fit disk color profile is broadly consistent with the Shakura \& Sunyaev disk model. We compare our UV-optical lags to the disk sizes inferred from optical-optical lags of the same quasars and find that our results are consistent with these quasars being drawn from a limited high-lag subset of the broader population. Our results are therefore broadly consistent with models that suggest longer disk lags in a subset of quasars, for example, due to a nonzero size of the ionizing corona and/or magnetic heating contributing to the disk response.

preprint2022arXiv

The Twin Radio galaxy TRG J104454+354055

We report observations of a bright dumb-bell system of galaxies ($z =0.162$) with the upgraded Giant Metrewave Radio Telescope ($u$GMRT), which show that each member of this gravitationally bound pair of galaxies hosts bipolar radio jets extended on 100 kiloparsec scales. Only two cases of such radio morphology have been reported previously, both being dumb-bell systems, as well. The famous first example, 3C 75, was discovered 4 decades ago, and the second case was discovered 3 decades ago. This implies that such `Twin-Radio-Galaxies' (TRGs) are an exceedingly rare phenomenon. As in the case of its two senior cousins, the bi-polar radio jets of the present TRG (J104454+354055) exhibit strong wiggles and are edge-darkened (Fanaroff-Riley class I). However, there are important differences, too. For instance, the jets in the present TRG do not merge and, moreover, show no signs of distortion due to an external crosswind. This makes the present TRG a much neater laboratory for studying the physics of (sideway) colliding jets of relativistic plasma. This TRG has a Wide-Angle-Tail (WAT) neighbour hosted by another bright galaxy belonging to the same group, which appears to be moving towards the TRG.

preprint2021arXiv

"A more probable explanation" is still impossible to explain GN-z11-flash: in response to Steinhardt et al. (arXiv:2101.12738)

In Jiang et al. (2020), we reported a possible bright flash (hereafter GN-z11-flash) from a galaxy GN-z11 at z ~ 11. Recently, Steinhardt et al. (2021; arXiv:2101.12738) found 27 images with transient signals in Keck MOSFIRE archival data and claimed that GN-z11-flash was more likely from a moving object in our Solar system. We show that the Steinhardt et al.'s definition of the chance probability and their methodology of finding GN-z11-flash-like transients are problematic in several aspects. In particular, none of their transients is analogous to GN-z11-flash, and none of them is positionally coincident with a known object in their imaging data. In Jiang et al., we performed a comprehensive analysis of the origin of GN-z11-flash and ruled out, to the best of our knowledge, the possibility of known man-made objects or moving objects in the Solar system, based on all available information and our current understanding of these objects. Steinhardt et al. did not use such information and did not analyse the GN-z11-flash event itself. The majority of their transients are apparently low-Earth orbit satellites or aircrafts. Therefore, their analysis can neither prove nor disprove our results. Finally, we present a method to estimate the chance probability of finding GN-z11-flash-like transients in archival data. Based on this method and the archival data used by Steinhardt et al., we obtain a loose upper limit of the probability that actually support the original results of Jiang et al. (2020).

preprint2021arXiv

A New Iron Emission Template for Active Galactic Nuclei. I. Optical Template for the H$β$ region

We present a new empirical template for iron emission in active galactic nuclei (AGN) covering the $4000-5600$ A range. The new template is based on a spectrum of the narrow-line Seyfert 1 galaxy Mrk 493 obtained with the Hubble Space Telescope. In comparison with the canonical iron template object I Zw 1, Mrk 493 has narrower broad-line widths, lower reddening, and a less extreme Eddington ratio, making it a superior choice for template construction. We carried out a multicomponent spectral decomposition to produce a template incorporating all permitted and forbidden lines of Fe II identified in the Mrk 493 spectrum over this wavelength range, as well as lines from Ti II, Ni II, and Cr II. We tested the template by fitting it to AGN spectra spanning a broad range of iron emission properties, and we present a detailed comparison with fits using other widely used monolithic and multi-component iron emission templates. The new template generally provides the best fit (lowest $χ^2$) compared to other widely used monolithic empirical templates. In addition, the new template yields more accurate spectral measurements including a significantly better match of the derived Balmer line profiles (H$β$, H$γ$, H$δ$), in contrast with results obtained using the other templates. Our comparison tests show that the choice of iron template can introduce a systematic bias in measurements of the H$β$ line width, which consequently impacts single-epoch black hole mass estimates by $\sim0.1$ dex on average and possibly up to $\sim0.3-0.5$ dex individually.

preprint2021arXiv

Evidence for the connection between star formation rate and evolutionary phases of quasars

Both theory and observations suggest that outflows driven by an active central supermassive black hole (SMBH) has a feedback effect on shaping the global properties of the host galaxy. However, whether feedback from the outflow is effective, and if so, whether it is positive or negative, has long been controversial. Here, using the latest catalog from the Sloan Digital Sky Survey (SDSS), we use the flux ratio of the [O II] to [Ne V] emission lines as a proxy to compare the star formation rate (SFR) in the hosts of quasars with different types of broad absorption lines (BALs): low-ionization (Lo)BAL, high-ionization (Hi)BAL, and non-BAL. We find that SFR decreases from LoBAL to HiBAL quasars, and then increases from HiBAL to non-BAL quasars. Assuming that the sequence of LoBAL to HiBAL to non-BAL represents evolution, our results are consistent with a quenching and subsequent rebound of star formation in quasar host galaxies. This phenomenon can be explained that the SFR is suppressed by the outflow, which then rebounds once the outflow disappears as the quasars evolve from HiBALs to non-BALs. Our result suggests that the quasar outflow has a negative global feedback on galaxy evolution.

preprint2021arXiv

Metallicity in Quasar Broad Line Regions at Redshift $\sim$ 6

Broad line regions (BLRs) in high-redshift quasars provide crucial information of chemical enrichment in the early universe. Here we present a study of BLR metallicities in 33 quasars at redshift $5.7<z<6.4$. Using the near-IR spectra of the quasars obtained from the Gemini telescope, we measure their rest-frame UV emission line flux and calculate flux ratios. We then estimate BLR metallicities with empirical calibrations based on photoionization models. The inferred median metallicity of our sample is a few times the solar value, indicating that the BLR gas had been highly metal-enriched at $z\sim6$. We compare our sample with a low-redshift quasar sample with similar luminosities and find no evidence of redshift evolution in quasar BLR metallicities. This is consistent with previous studies. The Fe II$/$Mg II flux ratio, a proxy for the Fe$/α$ element abundance ratio, shows no redshift evolution as well, further supporting rapid nuclear star formation at $z\sim6$. We also find that the black hole mass-BLR metallicity relation at $z\sim6$ is consistent with the relation measured at $2<z<5$, suggesting that our results are not biased by a selection effect due to this relation.

preprint2021arXiv

On the Mass Loading of AGN-Driven Outflows in Elliptical Galaxies and Clusters

Outflows driven by active galactic nuclei (AGN) are an important channel for accreting supermassive black holes (SMBHs) to interact with their host galaxies and clusters. Properties of the outflows are however poorly constrained due to the lack of kinetically resolved data of the hot plasma that permeates the circumgalactic and intracluster space. In this work, we use a single parameter, outflow-to-accretion mass-loading factor $m=\dot{M}_{\rm out}/\dot{M}_{\rm BH}$, to characterize the outflows that mediate the interaction between SMBHs and their hosts. By modeling both M87 and Perseus, and comparing the simulated thermal profiles with the X-ray observations of these two systems, we demonstrate that $m$ can be constrained between $200-500$. This parameter corresponds to a bulk flow speed between $4,000-7,000\,{\rm km\,s}^{-1}$ at around 1 kpc, and a thermalized outflow temperature between $10^{8.7}-10^{9}\,{\rm K}$. Our results indicate that the dominant outflow speeds in giant elliptical galaxies and clusters are much lower than in the close vicinity of the SMBH, signaling an efficient coupling with and deceleration by the surrounding medium on length scales below 1 kpc. Consequently, AGNs may be efficient at launching outflows $\sim10$ times more massive than previously uncovered by measurements of cold, obscuring material. We also examine the mass and velocity distribution of the cold gas, which ultimately forms a rotationally supported disk in simulated clusters. The rarity of such disks in observations indicates that further investigations are needed to understand the evolution of the cold gas after it forms.

preprint2021arXiv

SAGAN -- III: New insights into giant radio quasars

Giant radio quasars (GRQs) are radio-loud active galactic nuclei (AGNs), propelling megaparsec-scale jets. In order to understand GRQs and their properties, we have compiled all known GRQs ("the GRQ catalogue"), and a subset of small (size <700 kpc) radio quasars (SRQs) from the literature. In this process, we have found 10 new FR-II GRQs, in the redshift range of 0.66 < z < 1.72, which we include in the GRQ catalogue. Using the above samples, we have carried out a systematic comparative study of GRQs and SRQs, using optical and radio data. Our results show that the GRQs and SRQs statistically have similar spectral index and black hole mass distributions. However, SRQs have higher radio core power, core dominance factor, total radio power, jet kinetic power and Eddington ratio compared to GRQs. On the other hand, when compared to giant radio galaxies (GRGs), GRQs have higher black hole mass and Eddington ratio. The high core dominance factor of SRQs is an indicator of them lying closer to the line of sight than GRQs. We also find a correlation of the accretion disc luminosity with the radio core and jet power of GRQs, which provides evidence for disc-jet coupling. Lastly, we find the distributions of Eddington ratios of GRGs and GRQs to be bi-modal, similar to that found in small radio galaxies (SRGs) and SRQs, which indicate that size is not strongly dependent on the accretion state. Using all of these, we provide a basic model for the growth of SRQs to GRQs.

preprint2021arXiv

Star formation in 'the Brick': ALMA reveals an active proto-cluster in the Galactic centre cloud G0.253+0.016

G0.253+0.016, aka 'the Brick', is one of the most massive (> 10^5 Msun) and dense (> 10^4 cm-3) molecular clouds in the Milky Way's Central Molecular Zone. Previous observations have detected tentative signs of active star formation, most notably a water maser that is associated with a dust continuum source. We present ALMA Band 6 observations with an angular resolution of 0.13" (1000 AU) towards this 'maser core', and report unambiguous evidence of active star formation within G0.253+0.016. We detect a population of eighteen continuum sources (median mass ~ 2 Msun), nine of which are driving bi-polar molecular outflows as seen via SiO (5-4) emission. At the location of the water maser, we find evidence for a protostellar binary/multiple with multi-directional outflow emission. Despite the high density of G0.253+0.016, we find no evidence for high-mass protostars in our ALMA field. The observed sources are instead consistent with a cluster of low-to-intermediate-mass protostars. However, the measured outflow properties are consistent with those expected for intermediate-to-high-mass star formation. We conclude that the sources are young and rapidly accreting, and may potentially form intermediate and high-mass stars in the future. The masses and projected spatial distribution of the cores are generally consistent with thermal fragmentation, suggesting that the large-scale turbulence and strong magnetic field in the cloud do not dominate on these scales, and that star formation on the scale of individual protostars is similar to that in Galactic disc environments.

preprint2021arXiv

The Disk Veiling Effect of the Black Hole Low-Mass X-ray Binary A0620-00

The optical light curves of quiescent black hole low-mass X-ray binaries often exhibit significant non-ellipsoidal variabilities, showing the photospheric radiation of the companion star is veiled by other source of optical emission. Assessing this "veiling" effect is critical to the black hole mass measurement. Here in this work, we carry out a strictly simultaneous spectroscopic and photometric campaign on the prototype of black hole low-mass X-ray binary A0620-00. We find that for each observation epoch, the extra optical flux beyond a pure ellipsoidal modulation is positively correlated with the fraction of veiling emission, indicating the accretion disk contributes most of the non-ellipsoidal variations. Meanwhile, we also obtain a K2V spectral classification of the companion, as well as the measurements of the companion's rotational velocity $v \sin i = 83.8\pm1.9$ km s$^{-1}$ and the mass ratio between the companion and the black hole $q=0.063\pm0.004$.

preprint2021arXiv

The Ionization and Destruction of Polycyclic Aromatic Hydrocarbons in Powerful Quasars

We reanalyze the mid-infrared (5-40 $μ$m) Spitzer spectra of 86 low-redshift ($z < 0.5$) Palomar-Green quasars to investigate the nature of polycyclic aromatic hydrocarbon (PAH) emission and its utility as a star formation rate (SFR) indicator for the host galaxies of luminous active galactic nuclei (AGNs). We decompose the spectra with our recently developed template-fitting technique to measure PAH fluxes and upper limits, which we interpret using mock spectra that simulate the effects of AGN dilution. While luminous quasars can severely dilute and affect the detectability of emission lines, PAHs are intrinsically weak in some sources that are otherwise gas-rich and vigorously forming stars, conclusively demonstrating that powerful AGNs destroy PAH molecules. Comparing PAH-based SFRs with independent SFRs derived from the mid-infrared fine-structure neon lines and the total infrared luminosity reveals that PAHs can trace star formation activity in quasars with bolometric luminosities $\lesssim 10^{46}\, \rm erg\,s^{-1}$, but increasingly underestimate the SFR for more powerful quasars, typically by $\sim 0.5$ dex. Relative to star-forming galaxies and low-luminosity AGNs, quasars have a comparable PAH 11.3 $μ$m/7.7 $μ$m ratio but characteristically lower ratios of 6.2 $μ$m/7.7 $μ$m, 8.6 $μ$m/7.7 $μ$m, and 11.3 $μ$m/17.0 $μ$m. We suggest that these trends indicate that powerful AGNs preferentially destroy small grains and enhance the PAH ionization fraction.

preprint2021arXiv

The Relation between Morphological Asymmetry and Nuclear Activity in Low-redshift Galaxies

The morphology of galaxies reflects their assembly history and ongoing dynamical perturbations from the environment. Analyzing i-band images from the Pan-STARRS1 Survey, we study the optical morphological asymmetry of the host galaxies of a large, well-defined sample of nearby active galactic nuclei (AGNs) to investigate the role of mergers and interactions in triggering nuclear activity. The AGNs, comprising 245 type 1 and 4514 type 2 objects, are compared with 4537 star-forming galaxies matched in redshift and stellar mass. We develop a comprehensive masking strategy to isolate the emission of the target from foreground stars and other contaminating sources, all the while retaining projected companions of comparable brightness that may be major mergers. Among three variants of nonparametric indices, both the popular CAS asymmetry parameter and the outer asymmetry parameter (A_outer) yield robust measures of morphological distortion for star-forming galaxies and type 2 AGNs, while only A_outer is effective for type 1 AGNs. The shape asymmetry, by comparison, is affected more adversely by background noise. Asymmetry indices > 0.4 effectively trace systems that are candidate ongoing mergers. Contrary to theoretical expectations, galaxy interactions and mergers are not the main drivers of nuclear activity, at least not in our sample of low-redshift, relatively low-luminosity AGNs, whose host galaxies are significantly less asymmetric than the control sample of star-forming galaxies. Moreover, type 2 AGNs are morphologically indistinguishable from their type 1 counterparts. The level of AGN activity does not correlate with asymmetry, not even among the major merger candidates. As a by-product, we find, consistent with previous studies, that the average asymmetry of star-forming galaxies increases above the main sequence, although not all major mergers exhibit enhanced star formation.

preprint2021arXiv

What is Important? Morphological Asymmetries are Useful Predictors of Star Formation Rates of Star-forming Galaxies in SDSS Stripe 82

Morphology and structure of galaxies reflect their star formation and assembly histories. We use the framework of mutual information ($\mathrm{MI}$) to quantify interdependence among several structural variables and to rank them according to their relevance for predicting specific star formation rate (SSFR) by comparing the $\mathrm{MI}$ of the predictor variables with SSFR and penalizing variables that are redundant. We apply this framework to study $\sim 3,700$ face-on star-forming galaxies (SFGs) with varying degrees of bulge dominance and central concentration and with stellar mass $M_\star \approx 10^9\,M_\odot - 5\times 10^{11}\,M_\odot$ at redshift $z = 0.02-0.12$. We use the Sloan Digital Sky Survey (SDSS) Stripe 82 deep $i$-band imaging data, which improve measurements of asymmetry and bulge dominance indicators. We find that star-forming galaxies are a multi-parameter family. In addition to $M_\star$, asymmetry emerges as the most powerful predictor of SSFR residuals of SFGs, followed by bulge prominence/concentration. Star-forming galaxies with higher asymmetry and stronger bulges have higher SSFR at a given $M_\star$. The asymmetry reflects both irregular spiral arms and lopsidedness in seemingly isolated SFGs and structural perturbations by galaxy interactions or mergers.

preprint2020arXiv

AGN Feedback and Star Formation of Quasar Host Galaxies: Insights from the Molecular Gas

The molecular gas serves as a key probe of the complex interplay between black hole accretion and star formation in the host galaxies of active galactic nuclei (AGNs). We use CO(2-1) observations from a new ALMA survey, in conjunction with literature measurements, to investigate the molecular gas properties of a representative sample of 40 z<0.3 Palomar-Green quasars, the largest and most sensitive study of molecular gas emission to date for nearby quasars. We find that the AGN luminosity correlates with both the CO luminosity and black hole mass, suggesting that AGN activity is loosely coupled to the cold gas reservoir of the host. The observed strong correlation between host galaxy total infrared luminosity and AGN luminosity arises from their common dependence on the molecular gas. We argue that the total infrared luminosity, at least for low-redshift quasars, can be used to derive reliable star formation rates for the host galaxy. The host galaxies of low-redshift quasars have molecular gas content similar to that of star-forming galaxies of comparable stellar mass. Moreover, they share similar gas kinematics, as evidenced by their CO Tully-Fisher relation and the absence of detectable molecular outflows down to sensitive limits. There is no sign that AGN feedback quenches star formation for the quasars in our sample. On the contrary, the abundant gas supply forms stars prodigiously, at a rate that places most of them above the star-forming main sequence and with an efficiency that rivals that of starburst systems.

preprint2020arXiv

Black Hole Accretion Correlates with Star Formation Rate and Star Formation Efficiency in Nearby Luminous Type 1 Active Galaxies

We investigate the relationship between black hole accretion and star formation in a sample of 453 $z\approx0.3$ type 1 active galactic nuclei (AGNs). We use available CO observations to demonstrate that the combination of nebular dust extinction and metallicity provides reliable estimates of the molecular gas mass even for the host galaxies of type 1 AGNs. Consistent with other similar but significantly smaller samples, we reaffirm the notion that powerful AGNs have comparable gas content as nearby star-forming galaxies and that AGN feedback does not deplete the host of cold gas instantaneously. We demonstrate that while the strong correlation between star formation rate and black hole accretion rate is in part driven by the mutual dependence of these parameters on molecular gas mass, the star formation rate and black hole accretion rate are still weakly correlated after removing the dependence of star formation rate on molecular gas mass. This, together with a positive correlation between star formation efficiency and black hole accretion rate, may be interpreted as evidence for positive AGN feedback.

preprint2020arXiv

Black Hole Mass Measurements of Radio Galaxies NGC 315 and NGC 4261 Using ALMA CO Observations

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 5 and Cycle 6 observations of CO(2$-$1) and CO(3$-$2) emission at 0.2''$-$0.3'' resolution in two radio-bright, brightest group/cluster early-type galaxies, NGC 315 and NGC 4261. The data resolve CO emission that extends within their black hole (BH) spheres of influence ($r_\mathrm{g}$), tracing regular Keplerian rotation down to just tens of parsecs from the BHs. The projected molecular gas speeds in the highly inclined ($i>60^\circ$) disks rises at least 500 km s$^{-1}$ near their galaxy centers. We fit dynamical models of thin-disk rotation directly to the ALMA data cubes, and account for the extended stellar mass distributions by constructing galaxy surface brightness profiles corrected for a range of plausible dust extinction values. The best-fit models yield $(M_\mathrm{BH}/10^9\,M_\odot)=2.08\pm0.01(\mathrm{stat})^{+0.32}_{-0.14}(\mathrm{sys})$ for NGC 315 and $(M_\mathrm{BH}/10^9\,M_\odot)=1.67\pm0.10(\mathrm{stat})^{+0.39}_{-0.24}(\mathrm{sys})$ for NGC 4261, the latter of which is larger than previous estimates by a factor of $\sim$3. The BH masses are broadly consistent with the relations between BH masses and host galaxy properties. These are among the first ALMA observations to map dynamically cold gas kinematics well within the BH-dominated regions of radio galaxies, resolving the respective $r_\mathrm{g}$ by factors of $\sim$5$-$10. The observations demonstrate ALMA's ability to precisely measure BH masses in active galaxies, which will enable more confident probes of accretion physics for the most massive galaxies.

preprint2020arXiv

Broad-line Region of the Quasar PG 2130+099 from a Two-Year Reverberation Mapping Campaign with High Cadence

As one of the most interesting Seyfert 1 galaxies, PG 2130+099 has been the target of several reverberation mapping (RM) campaigns over the years. However, its measured broad H$β$ line responses have been inconsistent, with time lags of $\sim$200 days, $\sim$25 days, and $\sim$10 days being reported for different epochs while its optical luminosity changed no more than 40%. To investigate this issue, we conducted a new RM-campaign with homogenous and high cadence (about $\sim$3 days) for two years during 2017--2019 to measure the kinematics and structure of the ionized gas. We successfully detected time lags of broad H$β$, He II, He I, and Fe II lines with respect to the varying 5100Åcontinuum, revealing a stratified structure that is likely virialized with Keplerian kinematics in the first year of observations, but an inflow kinematics of the broad-line region from the second year. With a central black hole mass of $0.97_{-0.18}^{+0.15}\times 10^7~M_{\odot}$, PG 2130+099 has an accretion rate of $10^{2.1\pm0.5}L_{\rm Edd}c^{-2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is speed of light, implying that it is a super-Eddington accretor and likely possesses a slim, rather than thin, accretion disk. The fast changes of the ionization structures of the three broad lines remain puzzling.

preprint2020arXiv

Catching butterflies in the sky: Extended catalog of winged or X-shaped radio sources from the latest FIRST data release

We present a catalog of 290 "winged" or X-shaped radio galaxies (XRGs) extracted from the latest (2014 December 17) data release of the "Very Large Array Faint Images of the Radio Sky at Twenty centimeter." We have combined these radio images with their counterparts in the TIFR GMRT sky survey at 150 MHz, in an attempt to identify any low surface brightness radio emission present in these sources. This has enabled us to assemble a sample of 106 "strong" XRG candidates and 184 "probable" XRG candidates whose XRG designation needs to be verified by further observations. The present sample of 290 XRG candidates is almost twice as large as the number of XRGs currently known. Twenty-five of our 290 XRG candidates (9 "strong" and 16 "probable") are identified as quasars. Double-peaked narrow emission lines are seen in the optical spectra of three of the XRG candidates (two "strong" and one "probable"). Nearly 90% of the sample is located in the FR II domain of the Owen-Ledlow diagram. A few of the strong XRG candidates have a rather flat radio spectrum (spectral index alpha flatter than -0.3) between 150 MHz and 1.4 GHz, or between 1.4 and 5 GHz. Since this is not expected for lobe-dominated extragalactic radio sources (like nearly all known XRGs), these sources are particularly suited for follow-up radio imaging and near-simultaneous measurement of the radio spectrum.

preprint2020arXiv

Chandra Survey of Nearby Galaxies: an Extended Catalog

She et al. (Paper I) assembled a catalog of nearby galaxies observed with the Chandra X-ray observatory, by cross-matching galaxies in the NASA Extragalactic Database (NED) within 50 Mpc and the Chandra archive. That sample has enabled searches of low-mass black holes associated with late-type, bulgeless galaxies and studies of the accretion physics related to low-luminosity active galactic nuclei (LLAGNs). Using a similar approach, here we construct an extended catalog up to 150 Mpc and make a cross-correlation with a catalog of nearby galaxy groups. The new catalog consists of 1,964 galaxies, out of which 1,692 have a redshift independent distance, 1,557 are listed in the galaxy group catalog with group properties available, and 782 are identified to be X-ray AGN candidates. Compared with the AGN sample in Paper I, the new sample is 2.5 times larger in size (782 vs. 314), with ~80% of the new members having an Eddington ratio less than 10^-4. We confirm that the conclusions based on the previous sample remain. With the new sample, we compare AGN fractions between early-type and late-type galaxies, and between central and satellite galaxies in groups, and find no significant difference. This suggests that the secular process is not the dominant mechanism feeding AGNs in the local universe.

preprint2020arXiv

CMZoom: Survey Overview and First Data Release

We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular hydrogen column density of 10^23 cm^-2 at a resolution of ~3" (0.1 pc). In this paper, we focus on the 1.3 mm dust continuum and its data release, but also describe CMZoom spectral line data which will be released in a forthcoming publication. While CMZoom detected many regions with rich and complex substructure, its key result is an overall deficit in compact substructures on 0.1 - 2 pc scales (the compact dense gas fraction: CDGF). In comparison with clouds in the Galactic disk, the CDGF in the CMZ is substantially lower, despite having much higher average column densities. CMZ clouds with high CDGFs are well-known sites of active star formation. The inability of most gas in the CMZ to form compact substructures is likely responsible for the dearth of star formation in the CMZ, surprising considering its high density. The factors responsible for the low CDGF are not yet understood but are plausibly due to the extreme environment of the CMZ, having far-reaching ramifications for our understanding of the star formation process across the cosmos.

preprint2020arXiv

Correlation of structure and stellar properties of galaxies in Stripe 82

Establishing a correlation (or lack thereof) between the bimodal colour distribution of galaxies and their structural parameters is crucial to understand the origin of bimodality. To achieve that, we have performed 2D mass-based structural decomposition (bulge+disc) of all disc galaxies (total$=$1263) in the Herschel imaging area of the Stripe 82 region using $K_s$ band images from the VICS82 survey. The scaling relations thus derived are found to reflect the internal kinematics and are employed in combination to select an indubitable set of classical and pseudo bulge hosting disc galaxies. The rest of the galaxies ($<20\%$) are marked as discs with "ambiguous" bulges. Pseudo and classical bulge disc galaxies exhibit clear bimodality in terms of all stellar parameters ($M_*$, sSFR, $r-K_s$). All pseudo bulge disc galaxies are blue and star-forming and all classical bulge disc galaxies are red and quiescent with less than $5\%$ digressions. Ambiguous bulge disc galaxies are intermittent to pseudo and classical bulge disc galaxies in the distribution of all structural and stellar parameters. $Δ$$\langleμ_{eb}\rangle$ - based on the placement of bulges on the Kormendy relation - is found to be the most efficient single structural indicator of both bulge type and stellar activity. The placement of ambiguous bulge disc galaxies on scaling relations and fundamental plane, in addition to their peculiar stellar properties suggest that they are dominantly a part of the green valley.

preprint2020arXiv

Discovery of a [CI]-faint, CO-bright Galaxy: ALMA Observations of the Merging Galaxy NGC 6052

We report sensitive [\ion{C}{1}]~$^3P_1$--$^3P_0$ and $^{12}$CO~$J$=4--3 observations of the nearby merging galaxy NGC 6052 using the Morita (Atacama Compact) Array of ALMA. We detect $^{12}$CO~$J$=4--3 toward the northern part of NGC 6052, but [\ion{C}{1}]~$^3P_1$--$^3P_0$ is not detected with a [\ion{C}{1}]~$^3P_1$--$^3P_0$ to $^{12}$CO~$J$=4--3 line luminosity ratio of$~\lesssim0.07$. According to models of photodissociation regions, the unusual weakness of [\ion{C}{1}]~$^3P_1$--$^3P_0$ relative to $^{12}$CO~$J$=4--3 can be explained if the interstellar medium has a hydrogen density larger than $10^5\,{\rm cm}^{-3}$, conditions that might arise naturally in the ongoing merging process in NGC 6052. Its [\ion{C}{1}]~$^3P_1$--$^3P_0$ emission is also weaker than expected given the molecular gas mass inferred from previous measurements of $^{12}$CO~$J$=1--0 and $^{12}$CO~$J$=2--1. This suggests that [\ion{C}{1}]~$^3P_1$--$^3P_0$ may not be a reliable tracer of molecular gas mass in this galaxy. NGC 6052 is a unique laboratory to investigate how the merger process impacts the molecular gas distribution.

preprint2020arXiv

Discovery of a damped Ly$α$ galaxy at z $\sim$ 3 towards the quasar SDSS J011852+040644

We report the detection of the host galaxy of a damped Ly$α$ system (DLA) with log N(HI) $ [\rm cm^{-2}]$ = $21.0 \pm 0.10$ at $z \approx 3.0091$ towards the background quasar SDSS J011852+040644 using the Palomar Cosmic Web Imager (PCWI) at the Hale (P200) telescope. We detect Ly$α$ emission in the dark core of the DLA trough at a 3.3$σ$ confidence level, with Ly$α$ luminosity of $L_{\rm Lyα}$ $\rm = (3.8 \pm 0.8) \times 10^{42}\ erg\ s^{-1}$, corresponding to a star formation rate of $\gtrsim 2\ \rm M_{\odot}\ yr^{-1}$ (considering a lower limit on Ly$α$ escape fraction $f_{esc}^{Lyα} \sim 2\%$) as typical for Lyman break galaxies at these redshifts. The Ly$α$ emission is blueshifted with respect to the systemic redshift derived from metal absorption lines by $281 \pm 43$ km/s. The associated galaxy is at very small impact parameter of $\lesssim 12 \rm\ kpc$ from the background quasar, which is in line with the observed anticorrelation between column density and impact parameter in spectroscopic searches tracing the large-scale environments of DLA host galaxies.

preprint2020arXiv

From haloes to galaxies -- II. The fundamental relations in star formation and quenching

Star formation and quenching are two of the most important processes in galaxy formation and evolution. We explore in the local Universe the interrelationships among key integrated galaxy properties, including stellar mass $M_*$, star formation rate (SFR), specific SFR (sSFR), molecular gas mass $M_{\rm H_2}$, star formation efficiency (SFE) of the molecular gas and molecular gas to stellar mass ratio $μ$. We aim to identify the most fundamental scaling relations among these key galaxy properties and their interrelationships. We show the integrated $M_{\rm H_2}$-SFR, SFR-$M_*$ and $M_{\rm H_2}$-$M_*$ relation can be simply transformed from the $μ$-sSFR, SFE-$μ$ and SFE-sSFR relation, respectively. The transformation, in principle, can increase or decrease the scatter of each relation. Interestingly, we find the latter three relations all have significantly smaller scatter than the former three corresponding relations. We show the probability to achieve the observed small scatter by accident is extremely close to zero. This suggests that the smaller scatters of the latter three relations are driven by a more fundamental physical connection among these quantities. We then show the large scatters in the former relations are due to their systematic dependence on other galaxy properties, and on star formation and quenching process. We propose the sSFR-$μ$-SFE relation as the Fundamental Formation Relation (FFR), which governs the star formation and quenching process, and provides a simple framework to study galaxy evolution. Other scaling relations, including integrated Kennicutt-Schmidt law, star-forming main sequence and molecular gas main sequence, can all be derived from the FFR.

preprint2020arXiv

Gemini Multi-Object Spectrograph Integral Field Unit Spectroscopy of the Double-peaked Broad Emission Line of a Red Active Galactic Nucleus

Galaxy mergers are expected to produce multiple supermassive black holes (SMBHs) in close-separation, but the detection of such SMBHs has been difficult. 2MASS J165939.7$+$183436 is a red active galactic nucleus (AGN) that is a prospective merging SMBH candidate owing to its merging features in Hubble Space Telescope imaging and double-peaked broad emission lines (BELs). Herein, we report a Gemini Multi-Object Spectrograph Integral Field Unit observation of a double-peaked broad H$α$ line of 2MASS J165939.7$+$183436. Furthermore, we confirm the existence of two BEL peaks that are kinematically separated by 3000\,$\rm km\,s^{-1}$, with the SMBH of each BEL component weighing at $10^{8.92\pm0.06}\,M_{\rm \odot}$ and $10^{7.13\pm0.06}\,M_{\rm \odot}$, if they arise from independent BELs near the two SMBHs. The BEL components were not separated at $>0\farcs1$; however, under several plausible assumptions regarding the fitting of each spaxel, the two components are found to be spatially separated at $0\farcs085$ ($\sim250$\,pc). Different assumptions for the fitting can lead to a null ($< 0\farcs05$) or a larger spatial separation ($\sim0\farcs15$). Given the uncertainty regarding the spatial separation, various models, such as the disk emitter and multiple SMBH models, are viable solutions to explain the double BEL components. These results will promote future research for finding more multiple SMBH systems in red AGNs, and higher-resolution imaging validates these different models.

preprint2020arXiv

Global Spiral Density Wave Modes in Protoplanetary Disks: Morphology of Spiral Arms

We analyze two-armed global spiral density wave modes generated by gravitational instability in razor-thin, non-viscous, self-gravitating protoplanetary disks to understand the dependence of spiral arm morphology (pitch angle $α$ and amplitude) on various disk conditions. The morphologies of the resulting spiral density wave modes closely resemble observations. Their pitch angles and pattern speeds are insensitive to the boundary conditions adopted. Gaussian disks exhibit more tightly wound spirals (smaller pitch angle) than power law disks under the same conditions. We find that at a fixed disk-to-star mass ratio ($M_d/M_*$), pitch angle increases with average Toomre's stability parameter ($\overline Q$) or average disk aspect ratio ($\overline h$). For a given $\overline Q$, density wave modes with higher $M_d/M_*$ have larger pitch angles, while the behavior reverses for a given $\overline h$. The interdependence between pitch angle and disk properties can be roughly approximated by $α\propto c_s^2/M_d$, where $c_s$ is the sound speed. Our gravitational instability-excited spiral density waves can be distinguished from planet-launched spirals: (1) massive cool disks have spiral pitch angle falling with radius, while low-mass hot disks have spiral pitch angle rising with radius; (2) the profile of spiral amplitude presents several dips and bumps. We propose that gravitational instability-excited density waves can serve as an alternative scenario to explain the observed spiral arms in self-gravitating protoplanetary disks.

preprint2020arXiv

Kinematic decomposition of IllustrisTNG disk galaxies: morphology and relation with morphological structures

We recently developed an automated method, auto-GMM to decompose simulated galaxies. It extracts kinematic structures in an accurate, efficient, and unsupervised way. We use auto-GMM to study the stellar kinematic structures of disk galaxies from the TNG100 run of IllustrisTNG. We identify four to five structures that are commonly present among the diverse galaxy population. Structures having strong to moderate rotation are defined as cold and warm disks, respectively. Spheroidal structures dominated by random motions are classified as bulges or stellar halos, depending on how tightly bound they are. Disky bulges are structures that have moderate rotation but compact morphology. Across all disky galaxies and accounting for the stellar mass within 3 half-mass radii, the kinematic spheroidal structures, obtained by summing up stars of bulges and halos, contribute ~45% of the total stellar mass, while the disky structures constitute 55%. This study also provides important insights about the relationship between kinematically and morphologically derived galactic structures. Comparing the morphology of kinematic structures with that of traditional bulge+disk decomposition, we conclude: (1) the morphologically decomposed bulges are composite structures comprised of a slowly rotating bulge, an inner halo, and a disky bulge; (2) kinematically disky bulges, akin to what are commonly called pseudo bulges in observations, are compact disk-like components that have rotation similar to warm disks; (3) halos contribute almost 30% of the surface density of the outer part of morphological disks when viewed face-on; and (4) both cold and warm disks are often truncated in central regions.

preprint2020arXiv

KMTNet Nearby Galaxy Survey II. Searching for Dwarf Galaxies in Deep and Wide-field Images of the NGC 1291 system

We present newly discovered dwarf galaxy candidates in deep and wide-field images of NGC 1291 obtained with the Korea Microlensing Telescope Network. We identify 15 dwarf galaxy candidates by visual inspection. Using imaging simulations, we demonstrate that the completeness rate of our detection is greater than 70% for the central surface brightness value of $μ_{0,R}\lesssim26$ mag arcsec$^{-2}$ and for magnitudes $M_R\lesssim-10$ mag. The structural and photometric properties of the dwarf galaxy candidates appear to be broadly consistent with those of ordinary dwarf galaxies in nearby groups and clusters, with $μ_{0,R}\sim$ 22.5 to 26.5 mag arcsec$^{-2}$ and effective radii of 200 pc to 1 kpc. The dwarf galaxy candidates show a concentration towards NGC 1291 and tend to be redder the closer they are to the center, possibly indicating that they are associated with NGC 1291. The dwarf candidates presented in this paper appear to be bluer than those in denser environments, revealing that the quenching of star formation in dwarf galaxies is susceptible to the environment, while the morphology shaping is not.

preprint2020arXiv

On the Determination of Rotation Velocity and Dynamical Mass of Galaxies Based on Integrated H I Spectra

The integrated 21 cm H I emission profile of a galaxy encodes valuable information on the kinematics, spatial distribution, and dynamical state of its cold interstellar medium. The line width, in particular, reflects the rotation velocity of the galaxy, which, in combination with a size scale, can be used to constrain the dynamical mass of the system. We introduce a new method based on the concept of the curve of growth to derive a set of robust parameters to characterize the line width, asymmetry, and concentration of the integrated H I spectra. We use mock spectra to evaluate the performance of our method, to estimate realistic systematic uncertainties for the proposed parameters, and to correct the line widths for the effects of instrumental resolution and turbulence broadening. Using a large sample of nearby galaxies with available spatially resolved kinematics, we demonstrate that the newly defined line widths can predict the rotational velocities of galaxies to within an accuracy of $\lesssim 30$ km s$^{-1}$. We use the calibrated line widths, in conjunction with the empirical relation between the size and mass of H I disks, to formulate a prescription for estimating the dynamical mass within the H I-emitting region of gas-rich galaxies. Our formalism yields dynamical masses accurate to $\sim 0.3$ dex based solely on quantities that can be derived efficiently and robustly from current and future extragalactic H I surveys. We further extend the dynamical mass calibration to the scale of the dark matter halo.

preprint2020arXiv

SCUBA-2 Ultra Deep Imaging Eao Survey (Studies) III: Multi-wavelength properties, luminosity functions and preliminary source catalog of 450-$μ$m-selected galaxies

We construct a SCUBA-2 450-$μ$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$σ$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $μ$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$μ$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$β_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$μ$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

preprint2020arXiv

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) IV: Spatial clustering and halo masses of 450-$μ$m-selected sub-millimeter galaxies

We analyze an extremely deep 450-$μ$m image ($1σ=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$μ$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$μ$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$μ$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$μ$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$μ$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$μ$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

preprint2020arXiv

Stellar Properties of the Host Galaxy of an Ultraluminous X-ray Source in NGC 5252

An ultraluminous X-ray source (ULX) in NGC 5252 has been known as a strong candidate for an off-nuclear intermediate-mass black hole. We present near-infrared imaging data of the ULX obtained with the William Herschel Telescope. Using this data we estimate a stellar mass associated with the ULX of $\approx 10^{7.9\pm0.1}M_{\rm \odot}$ , suggesting that it could be (the remnant of) a dwarf galaxy that is in the process of merging with NGC 5252. Based on a correlation between the mass of the central black hole (BH) and host galaxy, the ULX is powered by a $10^5M_{\rm \odot}$ black hole. Alternatively, if the BH mass is $\approx 10^6M_{\rm \odot}$ or larger, the host galaxy of the ULX must have been heavily stripped during the merger. The ULX $K_s$-band luminosity is two orders of magnitude smaller than that expected from an ordinary active galactic nucleus with the observed [O III] luminosity, which also suggests the ULX lacks a dusty torus. We discuss how these findings provide suggestive evidence that the ULX is hosting an intermediate-mass black hole.

preprint2020arXiv

Strong Mg II and Fe II Absorbers at 2.2 < z < 6.0

We present a study of strong intervening absorption systems in the near-IR spectra of 31 luminous quasars at $z>5.7$. The quasar spectra were obtained with {\it Gemini} GNIRS that provide continuous wavelength coverage from $\sim$0.9 to $\sim$2.5 $μ$m. We detect 32 strong Mg II doublet absorbers with rest-frame equivalent width $W_r$ ($\lambda2796$) $>1.0$ Åat $2.2 < z < 6.0$. Each Mg II absorber is confirmed by at least two associated Fe II absorption lines in the rest-frame wavelength range of $\sim 1600-2600$ Å. We find that the comoving line density ($dN/dX$) of the strong Fe II-bearing Mg II absorbers decreases towards higher redshift at $z>3$, consistent with previous studies. Compared with strong Mg II absorbers detected in damped Ly$α$ systems at 2 $<z<$ 4, our absorbers are potentially less saturated and show much larger rest-frame velocity widths. This suggests that the gas traced by our absorbers are potentially affected by galactic superwinds. We analyze the {\it Hubble Space Telescope} near-IR images of the quasars and identify possible associated galaxies for our strong absorbers. There are a maximum of two galaxy candidates found within 5" radius of each absorber. The median F105W-band magnitude of these galaxy candidates is 24.8 mag, which is fainter than the $L^*$ galaxy luminosity at $z\sim$ 4. By using our observed $dN/dX$ of strong Mg II absorbers and galaxy candidates median luminosity, we suggest that at high redshift, strong Mg II absorbers tend to have a more disturbed environment but smaller halo size than that at $z <$ 1.

preprint2020arXiv

Supermassive black holes with high accretion rates in active galactic nuclei. XI. Accretion disk reverberation mapping of Mrk 142

We performed an intensive accretion disk reverberation mapping campaign on the high accretion rate active galactic nucleus Mrk 142 in early 2019. Mrk 142 was monitored with the Neil Gehrels Swift Observatory for 4 months in X-rays and 6 UV/optical filters. Ground-based photometric monitoring was obtained from the Las Cumbres Observatory, Liverpool Telescope and Dan Zowada Memorial Observatory in ugriz filters and the Yunnan Astronomical Observatory in V. Mrk 142 was highly variable throughout, displaying correlated variability across all wavelengths. We measure significant time lags between the different wavelength light curves, finding that through the UV and optical the wavelength-dependent lags, $τ(λ)$, generally follow the relation $τ(λ) \propto λ^{4/3}$, as expected for the $T\propto R^{-3/4}$ profile of a steady-state optically-thick, geometrically-thin accretion disk, though can also be fit by $τ(λ) \propto λ^{2}$, as expected for a slim disk. The exceptions are the u and U band, where an excess lag is observed, as has been observed in other AGN and attributed to continuum emission arising in the broad-line region. Furthermore, we perform a flux-flux analysis to separate the constant and variable components of the spectral energy distribution, finding that the flux-dependence of the variable component is consistent with the $f_ν\proptoν^{1/3}$ spectrum expected for a geometrically-thin accretion disk. Moreover, the X-ray to UV lag is significantly offset from an extrapolation of the UV/optical trend, with the X-rays showing a poorer correlation with the UV than the UV does with the optical. The magnitude of the UV/optical lags is consistent with a highly super-Eddington accretion rate.

preprint2020arXiv

The Carnegie-Irvine Galaxy Survey. IX. Classification of Bulge Types and Statistical Properties of Pseudo Bulges

We study the statistical properties of 320 bulges of disk galaxies in the Carnegie-Irvine Galaxy Survey, using robust structural parameters of galaxies derived from image fitting. We apply the Kormendy relation to classify classical and pseudo bulges and characterize bulge dichotomy with respect to bulge structural properties and physical properties of host galaxies. We confirm previous findings that pseudo bulges on average have smaller Sérsic indices, smaller bulge-to-total ratios, and fainter surface brightnesses when compared with classical bulges. Our sizable sample statistically shows that pseudo bulges are more intrinsically flattened than classical bulges. Pseudo bulges are most frequent (incidence $\gtrsim 80\%$) in late-type spirals (later than Sc). Our measurements support the picture in which pseudo bulges arose from star formation induced by inflowing gas, while classical bulges were born out of violent processes such as mergers and coalescence of clumps. We reveal differences with the literature that warrant attention: (1) the bimodal distribution of Sérsic indices presented by previous studies is not reproduced in our study; (2) classical and pseudo bulges have similar relative bulge sizes; and (3) the pseudo bulge fraction is considerably smaller in early-type disks compared with previous studies based on one-dimensional surface brightness profile fitting. We attribute the above differences to our improved image quality, more robust bulge-to-disk decomposition technique, and different classification criteria applied. Moreover, we find that barred galaxies do not host more pseudo bulges or more prominent pseudo bulges than unbarred galaxies. Various implications of these findings are discussed.

preprint2020arXiv

The Correlation between Black Hole Mass and Stellar Mass for Classical Bulges and the Cores of Ellipticals

The correlation between black hole mass and the stellar mass of the bulge of the host galaxy has attracted much attention ever since its discovery. While traditional investigations of this correlation have treated elliptical galaxies as single, monolithic spheroids, the recent realization that massive elliptical galaxies have undergone significant late-time ($z \lesssim 2$) dissipationless assembly since their initially dense "red nugget" phase strongly suggests that black holes in present-day ellipticals should be associated only with their cores and not with their extended envelopes. We perform two-dimensional image decomposition of Two Micron All Sky Survey $K_s$-band images to derive the stellar mass of the cores of 35 nearby ellipticals with reliably measured black hole masses. We revisit the relation between black hole mass and bulge stellar mass by combining classical bulges with the cores of ellipticals. The new relation exhibits nearly identical slope ($M_{\bullet} \propto M_{\text{core}}^{1.2}$) as the conventional relation but a factor of $\sim2$ higher normalization and moderately larger intrinsic scatter (0.4 dex). At a core mass of $10^{11}\,M_{\odot}$, $M_{\bullet}/M_{\text{core}} = 0.9\%$, but it rises to $M_{\bullet}/M_{\text{core}}=1.5\%$ for the most massive cores with mass $10^{12}\,M_{\odot}$. Fast and slow rotator ellipticals follow the same correlation. The $M_{\bullet}-M_{\text{core}}$ relation provides a revised benchmark for studies of black hole-galaxy coevolution in the high-redshift Universe.

preprint2020arXiv

The Formation History of Subhalos and the Evolution of Satellite Galaxies

Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the subhalos before they are accreted by massive host halos. We find that satellite galaxies in young subhalos (low $a_{\rm nf}$) are less massive and more gas rich, and have stronger star formation and a higher fraction of ex situ stellar mass than satellites in old subhalos (high $a_{\rm nf}$). Furthermore, these low $a_{\rm nf}$ satellites require longer timescales to be quenched as a population than the high $a_{\rm nf}$ counterparts. We find very different merger histories between satellites in fast accretion (FA, $a_{\rm nf}<1.3$) and slow accretion (SA, $a_{\rm nf}>1.3$) subhalos. For FA satellites, the galaxy merger frequency dramatically increases just after accretion, which enhances the star formation at accretion. While, for SA satellites, the mergers occur smoothly and continuously across the accretion time. Moreover, mergers with FA satellites happen mainly after accretion, while a contrary trend is found for SA satellites. Our results provide insight into the evolution and star formation quenching of the satellite population.

preprint2020arXiv

The MALATANG Survey: Dense Gas and Star Formation from High Transition HCN and HCO+ maps of NGC253

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $Σ_\text{stellar}$, and SFE$_\text{dense}$ versus $Σ_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $Σ_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars.

preprint2020arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Photometric g and i Light Curves

The Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) program monitors 849 active galactic nuclei (AGN) both spectroscopically and photometrically. The photometric observations used in this work span over four years and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory's Bok telescope and the CFHT telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN's variability. We find that the Welch-Stetson J-index provides a useful characterization of AGN variability and can be used to select AGNs for further study.

preprint2020arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: The H$β$ Radius-Luminosity Relation

Results from a few decades of reverberation mapping (RM) studies have revealed a correlation between the radius of the broad-line emitting region (BLR) and the continuum luminosity of active galactic nuclei. This "radius-luminosity" relation enables survey-scale black-hole mass estimates across cosmic time, using relatively inexpensive single-epoch spectroscopy, rather than intensive RM time monitoring. However, recent results from newer reverberation mapping campaigns challenge this widely used paradigm, reporting quasar BLR sizes that differ significantly from the previously established radius-luminosity relation. Using simulations of the radius--luminosity relation with the observational parameters of SDSS-RM, we find that this difference is not likely due to observational biases. Instead, it appears that previous RM samples were biased to a subset of quasar properties, and the broader parameter space occupied by the SDSS-RM quasar sample has a genuinely wider range of BLR sizes. We examine the correlation between the deviations from the radius-luminosity relation and several quasar parameters; the most significant correlations indicate that the deviations depend on UV/optical SED and the relative amount of ionizing radiation. Our results indicate that single-epoch black-hole mass estimates that do not account for the diversity of quasars in the radius-luminosity relation could be overestimated by an average of ~0.3 dex.

preprint2020arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: The M_BH-Host Relations at 0.2<z<0.6 from Reverberation Mapping and Hubble Space Telescope Imaging

We present the results of a pilot Hubble Space Telescope (HST) imaging study of the host galaxies of ten quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Probing more than an order of magnitude in BH and stellar masses, our sample is the first statistical sample to study the BH-host correlations beyond z>0.3 with reliable BH masses from reverberation mapping rather than from single-epoch spectroscopy. We perform image decomposition in two HST bands (UVIS-F606W and IR-F110W) to measure host colors and estimate stellar masses using empirical relations between broad-band colors and the mass-to-light ratio. The stellar masses of our targets are mostly dominated by a bulge component. The BH masses and stellar masses of our sample broadly follow the same correlations found for local RM AGN and quiescent bulge-dominant galaxies, with no strong evidence of evolution in the M_BH-M_*,bulge relation to z~0.6. We further compare the host light fraction from HST imaging decomposition to that estimated from spectral decomposition. We found a good correlation between the host fractions derived with both methods. However, the host fraction derived from spectral decomposition is systematically smaller than that from imaging decomposition by ~30%, indicating different systematics in both approaches. This study paves the way for upcoming more ambitious host galaxy studies of quasars with direct RM-based BH masses at high redshift.

preprint2020arXiv

Universal transition diagram from dormant to actively accreting supermassive black holes

The vast majority of supermassive black holes (SMBHs) in the local universe exhibit levels of activity much lower than those expected from gas supplying rates onto the galactic nuclei, and only a small fraction of silent SMBHs can turn into active galactic nuclei. Revisiting observational data of very nearby SMBHs whose gravitational spheres of influence are spatially reached by the Chandra X-ray satellite, we find that the level of BH activity drastically increases from the quiescent phase when the inflow rate outside of the BH influence radius is higher than 0.1% of the Eddington accretion rate. We also show that the relation between the nuclear luminosity and gas accretion rate from the BH influence radius measured from X-ray observations is well described by the universal state transition of accreting SMBHs, as predicted by recent hydrodynamical simulations with radiative cooling and BH feedback. After the state transition, young massive stars should form naturally in the nucleus, as observed in the case of the nearest SMBH, Sagittarius A$^\ast$, which is currently quiescent but was recently active.

preprint2019arXiv

A New Channel of Bulge Formation via The Destruction of Short Bars

Short (inner) bars of sub-kiloparsec radius have been hypothesized to be an important mechanism for driving gas inflows to small scales, thus feeding central black holes. Recent numerical simulations have shown that the growth of central black holes in galaxies can destroy short bars, when the black hole reaches a mass of $\sim 0.1\%$ of the total stellar mass of the galaxy. We study $N$-body simulations of galaxies with single and double bars to track the long-term evolution of the central stellar mass distribution. We find that the destruction of the short bar contributes significantly to the growth of the bulge. The final bulge mass is roughly equal to the sum of the masses of the initial pseudo bulge and short bar. The initially boxy/peanut-shaped bulge of Sérsic index $n\lesssim1$ is transformed into a more massive, compact structure that bears many similarities to a classical bulge, in terms of its morphology ($n \approx 2$), kinematics (dispersion-dominated, isotropic), and location on standard scaling relations (Kormendy relation, mass-size relation, and correlations between black hole mass and bulge stellar mass and velocity dispersion). Our proposed channel for forming classical bulges relies solely on the destruction of short bars without any reliance on mergers. We suggest that some of the less massive, less compact classical bulges were formed in this manner.

preprint2019arXiv

Active Galactic Nuclei with Ultra-fast Outflows Monitoring Project: The Broad-line Region of Mrk 79 as a Disk Wind

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$β~\lambda4861$, H$γ~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$β~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$γ~\lambda4340$, H$β~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

preprint2019arXiv

An ALMA CO(2-1) Survey of Nearby Palomar-Green Quasars

The properties of the molecular gas can shed light on the physical conditions of quasar host galaxies and the effect of feedback from accreting supermassive black holes. We present a new CO(2-1) survey of 23 z<0.1 Palomar-Green quasars conducted with the Atacama Large Millimeter/submillimeter Array. CO emission was successfully detected in 91% (21/23) of the objects, from which we derive CO luminosities, molecular gas masses, and velocity line widths. Together with CO(1-0) measurements in the literature for 32 quasars (detection rate 53%), there are 15 quasars with both CO(1-0) and CO(2-1) measurements and in total 40 sources with CO measurements. We find that the line ratio R_21 = L'_CO(2-1)/L'_CO(1-0) is subthermal, broadly consistent with nearby galaxies and other quasars previously studied. No clear correlation is found between R_21 and the intensity of the interstellar radiation field or the luminosity of the active nucleus. As with the general galaxy population, quasar host galaxies exhibit a strong, tight, linear L_IR-L_CO relation, with a normalization consistent with that of starburst systems. We investigate the molecular-to-total gas mass fraction with the aid of total gas masses inferred from dust masses previously derived from infrared observations. Although the scatter is considerable, the current data do not suggest that the CO-to-H_2 conversion factor of quasar host galaxies significantly differs from that of normal star-forming galaxies.

preprint2019arXiv

Kinematic signatures of reverberation mapping of close binaries of supermassive black holes in active galactic nuclei. II. Atlas of two-dimensional transfer functions

Most large galaxies harbor supermassive black holes (SMBHs) in their centers, and galaxies merge. Consequently, binary SMBHs should be common in galactic nuclei. However, close binaries of SMBH (CB-SMBHs) with sub-parsec separation cannot be imaged directly using current facilities. Some indirect signatures, such as periodic signals in light curves and double peaks in emission-line profile, have been used to find CB-SMBH candidates, but ambiguities still exist and no definitive conclusions can be made. We have recently proposed a new method focusing on kinematic signatures that can be derived from reverberation mapping of CB-SMBHs, one that offers a promising avenue to address this important problem. In this paper, we calculated models for a wide range of parameters, but BLRs of two BHs are close but still not merged. The purpose of this supplementary paper is to provide an atlas of two-dimensional transfer functions of CB-SMBHs with a wide range of orbital and geometrical parameters to aid more efficient identification of CB-SMBH candidates in reverberation mapping data.

preprint2019arXiv

X-shaped Radio Galaxies: Optical Properties, Large-scale Environment and Relationship to Radio Structure

In order to find clues to the origin of the "winged" or "X-shaped" radio galaxies (XRGs) we investigate here the parent galaxies of a large sample of 106 XRGs for optical-radio axes alignment, interstellar medium, black hole mass, and large-scale environment. For 41 of the XRGs it was possible to determine the optical major axis and the primary radio axis and the strong tendency for the two axes to be fairly close is confirmed. However, several counter-examples were also found and these could challenge the widely discussed backflow diversion model for the origin of the radio wings. Comparison with a well-defined large sample of normal FR II radio galaxies has revealed that: (i) XRGs possess slightly less massive central black holes than the normal radio galaxies (average masses being log$M_{\rm BH} \sim$ 8.81 $M_{\odot}$ and 9.07 $M_{\odot}$, respectively); (ii) a much higher fraction of XRGs ($\sim$ 80%) exhibits red mid-IR colors ($W2 - W3 > 1.5$), indicating a population of young stars and/or an enhanced dust mass, probably due to relatively recent galaxy merger(s). A comparison of the large-scale environment (i.e., within $\sim$ 1 Mpc) shows that both XRGs and FRII radio galaxies inhabit similarly poor galaxy clustering environments (medium richness being 8.94 and 11.87, respectively). Overall, the origin of XRGs seems difficult to reconcile with a single dominant physical mechanism and competing mechanisms seem prevalent.

preprint2018arXiv

On the Connection Between Spiral Arm Pitch Angle and Galaxy Properties

We measure the pitch angle ($φ$) of spiral arms in a sample of 79 galaxies to perform a systematic study of the dependence of $φ$ on galaxy morphology, mass, and kinematics to investigate the physical origin of spiral arms. We find that $φ$ decreases (arms are more tightly wound), albeit with significant scatter, in galaxies with earlier Hubble type, more prominent bulges, higher concentration, and larger total galaxy stellar mass ($M_*^{\rm gal}$). For a given concentration, galaxies with larger stellar masses tend to have tighter spiral arms, and vice versa. We also find that $φ$ obeys a tight inverse correlation with central stellar velocity dispersion for $σ_c$$\gtrsim$$100$ km s$^{-1}$, whereas $φ$ remains approximately constant for $σ_c\lesssim100$ km s$^{-1}$. We demonstrate that the $φ$-$σ_c$ and $φ$-$M_*^{\rm gal}$ relations are projections of a more fundamental three-dimensional $φ-σ_c-M_*^{\rm gal}$ relation, such that pitch angle is determined by $σ_c$ for massive galaxies but by $M_*^{\rm gal}$ for less massive galaxies. Contrary to previous studies, we find that $φ$ correlates only loosely with the galaxy's shear rate. For a given shear rate, spirals generated from $N$-body simulations exhibit much higher $φ$ than observed, suggesting that galactic disks are dynamically cooler (Toomre's $Q \approx 1.2$). Instead, the measured pitch angles show a much stronger relation with morphology of the rotation curve of the central region, such that galaxies with centrally peaked rotation curves have tight arms, while those with slow-rising rotation curves have looser arms. These behaviors are qualitatively consistent with predictions of density wave theory.

preprint2016arXiv

Chandra X-ray and Hubble Space Telescope Imaging of Optically Selected Kiloparsec-Scale Binary Active Galactic Nuclei II: Host Galaxy Morphology and AGN Activity

Binary active galactic nuclei (AGNs) provide clues to how gas-rich mergers trigger and fuel AGNs and how supermassive black hole (SMBH) pairs evolve in a gas-rich environment. While significant effort has been invested in their identification, the detailed properties of binary AGNs and their host galaxies are still poorly constrained. In a companion paper, we examined the nature of ionizing sources in the double nuclei of four kpc-scale binary AGNs with redshifts between 0.1~0.2. Here, we present their host galaxy morphology based on F336W (U-band) and F105W (Y-band) images taken by the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. Our targets have double-peaked narrow emission lines and were confirmed to host binary AGNs with follow up observations. We find that kpc-scale binary AGNs occur in galaxy mergers with diverse morphological types. There are three major mergers with intermediate morphologies and a minor merger with a dominant disk component. We estimate the masses of the SMBHs from their host bulge stellar masses and obtain Eddington ratios for each AGN. Compared with a representative control sample drawn at the same redshift and stellar mass, the AGN luminosities and Eddington ratios of our binary AGNs are similar to those of single AGNs. The U-Y color maps indicate that clumpy star forming regions could significantly affect the X-ray detection of binary AGNs, e.g., the hardness ratio. Considering the weak X-ray emission in AGNs triggered in merger systems, we suggest that samples of X-ray selected AGNs may be biased against gas-rich mergers.

preprint2016arXiv

Improving the Flux Calibration in Reverberation Mapping by Spectral Fitting: Application to the Seyfert Galaxy MCG--6-30-15

We propose a method for the flux calibration of reverberation mapping spectra based on accurate measurement of [O III] $λ5007$ emission by spectral fitting. The method can achieve better accuracy than the traditional method of van Groningen & Wanders (1992), allowing reverberation mapping measurements for object with variability amplitudes as low as $\sim$ 5%. As a demonstration, we reanalyze the data of the Seyfert 1 galaxy MCG--6-30-15 taken from the 2008 campaign of the Lick AGN Monitoring Project, which previously failed to obtain a time lag for this weakly variable object owing to a relatively large flux calibration uncertainty. We detect a statistically significant rest-frame time lag of $6.38_{-2.69}^{+3.07}$ days between the H$β$ and $V$-band light curves. Combining this lag with FWHM(H$β$) = $1933\pm81$ $\rm km~s^{-1}$ and a virial coefficient of $f$ = 0.7, we derive a virial black hole mass of $3.26_{-1.40}^{+1.59}\times10^6$ $M_{\odot}$, which agrees well with previous estimates by other methods.

preprint2016arXiv

Measurement of the Black Hole Mass in NGC 1332 from ALMA Observations at 0.044 Arcsecond Resolution

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3 observations of CO(2-1) emission from the circumnuclear disk in the E/S0 galaxy NGC 1332 at 0.044" resolution. The disk exhibits regular rotational kinematics and central high-velocity emission (+/-500 km/s) consistent with the presence of a compact central mass. We construct models for a thin, dynamically cold disk in the gravitational potential of the host galaxy and black hole, and fit the beam-smeared model line profiles directly to the ALMA data cube. Model fits successfully reproduce the disk kinematics out to r=200 pc. Fitting models just to spatial pixels within projected r=50 pc of the nucleus (two times larger than the black hole's gravitational radius of influence), we find M_BH=6.64(-0.63,+0.65)*10^8 solar masses. This observation demonstrates ALMA's powerful capability to determine the masses of supermassive black holes by resolving gas kinematics on small angular scales in galaxy nuclei.

preprint2016arXiv

NGC5252: a pair of radio-emitting active galactic nuclei?

The X-ray source CXO J133815.6+043255 has counterparts in the UV, optical, and radio bands. Based on the multi-band investigations, it has been recently proposed by Kim et al. (2015) as a rarely-seen off-nucleus ultraluminous X-ray (ULX) source with a black hole mass of >= 10^4 solar mass in the nearby Seyfert galaxy NGC 5252. To explore its radio properties at very high angular resolution, we performed very long-baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.7 GHz. We find that the radio counterpart is remarkably compact among the known ULXs. It does not show a resolved structure with a resolution of a few milliarcsecond (mas), and the total recovered flux density is comparable to that measured in earlier sub-arcsecond-resolution images. The compact radio structure, the relatively flat spectrum, and the high radio luminosity are consistent with a weakly accreting supermassive black hole in a low-luminosity active galactic nucleus. The nucleus of NGC 5252 itself has similar radio properties. We argue that the system represents a relatively rare pair of active galactic nuclei, where both components emit in the radio.

preprint2016arXiv

On the limits of measuring the bulge and disk properties of local and high-redshift massive galaxies

A considerable fraction of the massive quiescent galaxies at \emph{z} $\approx$ 2, which are known to be much more compact than galaxies of comparable mass today, appear to have a disk. How well can we measure the bulge and disk properties of these systems? We simulate two-component model galaxies in order to systematically quantify the effects of non-homology in structures and the methods employed. We employ empirical scaling relations to produce realistic-looking local galaxies with a uniform and wide range of bulge-to-total ratios ($B/T$), and then rescale them to mimic the signal-to-noise ratios and sizes of observed galaxies at \emph{z} $\approx$ 2. This provides the most complete set of simulations to date for which we can examine the robustness of two-component decomposition of compact disk galaxies at different $B/T$. We confirm that the size of these massive, compact galaxies can be measured robustly using a single Sérsic fit. We can measure $B/T$ accurately without imposing any constraints on the light profile shape of the bulge, but, due to the small angular sizes of bulges at high redshift, their detailed properties can only be recovered for galaxies with $B/T$ \gax\ 0.2. The disk component, by contrast, can be measured with little difficulty.

preprint2016arXiv

Reverberation Mapping of the Broad-line Region in NGC 5548: Evidence for Radiation Pressure?

NGC 5548 is the best-observed reverberation-mapped active galactic nucleus with long-term, intensive monitoring. Here we report results from a new observational campaign between January and July, 2015. We measure the centroid time lag of the broad H$β$ emission line with respect to the 5100 Åcontinuum and obtain $τ_{\rm cent} = 7.20^{+1.33}_{-0.35}$ days in the rest frame. This yields a black hole mass of $M_{\bullet}=8.71^{+3.21}_{-2.61} $x$ 10^{7}M_{\odot}$ using a broad H$β$ line dispersion of $3124\pm302$ km s$^{-1}$ and a virial factor of $f_{_{\rm BLR}}=6.3\pm1.5$ for the broad-line region (BLR), consistent with the mass measurements from previous H$β$ campaigns. The high-quality data allow us to construct a velocity-binned delay map for the broad H$β$ line, which shows a symmetric response pattern around the line center, a plausible kinematic signature of virialized motion of the BLR. Combining all the available measurements of H$β$ time lags and the associated mean 5100 Å luminosities over 18 campaigns between 1989 and 2015, we find that the H$β$ BLR size varies with the mean optical luminosity, but, interestingly, with a possible delay of $2.35_{-1.25}^{+3.47}$ yrs. This delay coincides with the typical BLR dynamical timescale of NGC 5548, indicating that the BLR undergoes dynamical changes, possibly driven by radiation pressure.

preprint2016arXiv

Spectroscopic Indication for a Centi-parsec Supermassive Black Hole Binary in the Galactic Center of NGC 5548

As a natural consequence of cosmological hierarchical structure formation, sub-parsec supermassive black hole binaries (SMBHBs) should be common in galaxies but thus far have eluded spectroscopic identification. Based on four decades of optical spectroscopic monitoring, we report that the nucleus of NGC 5548, a nearby Seyfert galaxy long suspected to have experienced a major merger about one billion years ago, exhibits long-term variability with a period of 14 years in the optical continuum and broad Hbeta emission line. Remarkably, the double-peaked profile of Hbeta shows systematic velocity changes with a similar period. These pieces of observations plausibly indicate that a SMBHB resides in the center of NGC 5548. The complex, secular variations in the line profiles can be explained by orbital motion of a binary with equal mass and a semi-major axis of 22 light-days (corresponding to 18 milli-parsec). At a distance of 75 Mpc, NGC 5548 is one of the nearest sub-parsec SMBHB candidates that offers an ideal laboratory for gravitational wave detection.

preprint2016arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. V. A New Size-Luminosity Scaling Relation for the Broad-Line Region

This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$β$ lags ($τ_{_{\rm Hβ}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm Hβ}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$β$ lags tightly depends on accretion rates, $τ_{_{\rm Hβ}}/τ_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $τ_{_{R-L}}$ is the H$β$ lag from the normal $R_{_{\rm Hβ}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm Hβ}}=α_1\ell_{44}^{β_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-γ_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 Å continuum luminosity, and coefficients of $α_1=(29.6_{-2.8}^{+2.7})$ lt-d, $β_1=0.56_{-0.03}^{+0.03}$, $γ_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.

preprint2016arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. VI. Velocity-resolved Reverberation Mapping of H$β$ Line

In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$β$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$β$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.

preprint2016arXiv

The Carnegie-Irvine Galaxy Survey. IV. A Method to Determine the Average Mass Ratio of Mergers That Built Massive Elliptical Galaxies

Many recent observations and numerical simulations suggest that nearby massive, early-type galaxies were formed through a "two-phase" process. In the proposed second phase, the extended stellar envelope was accumulated through many dry mergers. However, details of the past merger history of present-day ellipticals, such as the typical merger mass ratio, are difficult to constrain observationally. Within the context and assumptions of the two-phase formation scenario, we propose a straightforward method, using photometric data alone, to estimate the average mass ratio of mergers that contributed to the build-up of massive elliptical galaxies. We study a sample of nearby massive elliptical galaxies selected from the Carnegie-Irvine Galaxy Survey, using two-dimensional analysis to decompose their light distribution into an inner, denser component plus an extended, outer envelope, each having a different optical color. The combination of these two substructures accurately recovers the negative color gradient exhibited by the galaxy as whole. The color difference between the two components (<Δ(B-V)> ~ 0.10 mag; <Δ(B-R)> ~ 0.14 mag), based on the slope of the M_stellar-color relation for nearby early-type galaxies, can be translated into an estimate of the average mass ratio of the mergers. The rough estimate, 1:5 to 1:10, is consistent with the expectation of the two-phase formation scenario, suggesting that minor mergers were largely responsible for building up to the outer stellar envelope of present-day massive ellipticals. With the help of accurate photometry, large sample size, and more choices of colors promised by ongoing and future surveys, the approach proposed here can reveal more insights into the growth of massive galaxies during the last few Gyr.

preprint2016arXiv

The Fundamental Plane of the Broad-line Region in Active Galactic Nuclei

Broad emission lines in active galactic nuclei (AGNs) mainly arise from gas photoionized by continuum radiation from an accretion disk around a central black hole. The shape of the broad-line profile, described by ${\cal D}_{_{\rm Hβ}}={\rm FWHM}/σ_{_{\rm Hβ}}$, the ratio of full width at half maximum to the dispersion of broad H$β$, reflects the dynamics of the broad-line region (BLR) and correlates with the dimensionless accretion rate ($\dot{\mathscr{M}}$) or Eddington ratio ($L_{\rm bol}/L_{\rm Edd}$). At the same time, $\dot{\mathscr{M}}$ and $L_{\rm bol}/L_{\rm Edd}$ correlate with ${\cal R}_{\rm Fe}$, the ratio of optical Fe II to H$β$ line flux emission. Assembling all AGNs with reverberation mapping measurements of broad H$β$, both from the literature and from new observations reported here, we find a strong bivariate correlation of the form $\log(\dot{\mathscr{M}},L_{\rm bol}/L_{\rm Edd})=α+β{\cal D}_{_{\rm Hβ}}+γ{\cal R}_{\rm Fe},$ where $α=(2.47,0.31)$, $β=-(1.59,0.82)$ and $γ=(1.34,0.80)$. We refer to this as the fundamental plane of the BLR. We apply the plane to a sample of $z < 0.8$ quasars to demonstrate the prevalence of super-Eddington accreting AGNs are quite common at low redshifts.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: An Investigation of Biases in CIV Emission-Line Properties

We investigate the dependence on data quality of quasar properties measured from the CIV emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch and coadded spectra, focusing on the CIV emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ~2-4), we find increasing systematic offsets with decreasing S/N. The systematic difference (measurement uncertainty) in our lowest S/N (<5) subsample between the single-epoch and coadded spectrum (i) CIV equivalent width is 17A (31A), (ii) centroid wavelength is <1A (2A), and fractional velocity widths, Delta V/V, characterized by (iii) the line dispersion, sigma_l, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1-sigma measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the CIV full-width at half maximum (FWHM) velocity width and the CIV profile characterized by FWHM/sigma_l are only smaller than the statistical uncertainties when S/N>10, although offsets in lower S/N spectra exceed the statistical uncertainties by only a factor of ~1.5. Characterizing the CIV line profile by the kurtosis is the least robust property investigated, as the median systematic coadded--single-epoch measurement differences are larger than the statistical uncertainties for all S/N subsamples.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Biases in z>1.46 Redshifts due to Quasar Diversity

We use the coadded spectra of 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project observations of 482 quasars with z>1.46 to highlight systematic biases in the SDSS- and BOSS-pipeline redshifts due to the natural diversity of quasar properties. We investigate the characteristics of this bias by comparing the BOSS-pipeline redshifts to an estimate from the centroid of HeII 1640. HeII has a low equivalent width but is often well-defined in high-S/N spectra, does not suffer from self-absorption, and has a narrow component that, when present (the case for about half of our sources), produces a redshift estimate that, on average, is consistent with that determined from [OII] to within 1-sigma of the quadrature sum of the HeII and [OII] centroid measurement uncertainties. The large redshift differences of ~1000 km/s, on average, between the BOSS-pipeline and HeII-centroid redshifts suggest there are significant biases in a portion of BOSS quasar redshift measurements. Adopting the HeII-based redshifts shows that CIV does not exhibit a ubiquitous blueshift for all quasars, given the precision probed by our measurements. Instead, we find a distribution of CIV centroid blueshifts across our sample, with a dynamic range that (i) is wider than that previously reported for this line, and (ii) spans CIV centroids from those consistent with the systemic redshift to those with significant blueshifts of thousands of kilometers per second. These results have significant implications for measurement and use of high-redshift quasar properties and redshifts and studies based thereon.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Velocity Shifts of Quasar Emission Lines

Quasar emission lines are often shifted from the systemic velocity due to various dynamical and radiative processes in the line-emitting region. The level of these velocity shifts depends both on the line species and on quasar properties. We study velocity shifts for the line peaks of various narrow and broad quasar emission lines relative to systemic using a sample of 849 quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. The coadded (from 32 epochs) spectra of individual quasars have sufficient signal-to-noise ratio (SNR) to measure stellar absorption lines to provide reliable systemic velocity estimates, as well as weak narrow emission lines. The sample also covers a large dynamic range in quasar luminosity (~2 dex), allowing us to explore potential luminosity dependence of the velocity shifts. We derive average line peak velocity shifts as a function of quasar luminosity for different lines, and quantify their intrinsic scatter. We further quantify how well the peak velocity can be measured for various lines as a function of continuum SNR, and demonstrate there is no systematic bias in the line peak measurements when the spectral quality is degraded to as low as SNR~3 per SDSS pixel. Based on the observed line shifts, we provide empirical guidelines on redshift estimation from [OII]3728, [OIII]5008, [NeV]3426, MgII, CIII], HeII1640, broad Hbeta, CIV, and SiIV, which are calibrated to provide unbiased systemic redshifts in the mean, but with increasing intrinsic uncertainties of 46, 56, 119, 205, 233, 242, 400, 415, and 477 km/s, in addition to the measurement uncertainties. These more realistic redshift uncertainties are generally much larger than the formal uncertainties reported by the redshift pipelines for spectroscopic quasar surveys, and demonstrate the infeasibility of measuring quasar redshifts to better than ~200 km/s with only broad lines.

preprint2015arXiv

A Dynamical Study of the Black Hole X-ray Binary Nova Muscae 1991

We present a dynamical study of the Galactic black hole binary system Nova Muscae 1991 (GS/GRS 1124-683). We utilize 72 high resolution Magellan Echellette (MagE) spectra and 72 strictly simultaneous V-band photometric observations; the simultaneity is a unique and crucial feature of this dynamical study. The data were taken on two consecutive nights and cover the full 10.4-hour orbital cycle. The radial velocities of the secondary star are determined by cross-correlating the object spectra with the best-match template spectrum obtained using the same instrument configuration. Based on our independent analysis of five orders of the echellette spectrum, the semi-amplitude of the radial velocity of the secondary is measured to be K_2 = 406.8+/-2.7 km/s, which is consistent with previous work, while the uncertainty is reduced by a factor of 3. The corresponding mass function is f(M) = 3.02+/-0.06 M_\odot. We have also obtained an accurate measurement of the rotational broadening of the stellar absorption lines (v sin i = 85.0+/-2.6 km/s) and hence the mass ratio of the system q = 0.079+/-0.007. Finally, we have measured the spectrum of the non-stellar component of emission that veils the spectrum of the secondary. In a future paper, we will use our veiling-corrected spectrum of the secondary and accurate values of K_2 and q to model multi-color light curves and determine the systemic inclination and the mass of the black hole.

preprint2015arXiv

A Luminous X-ray Flare From The Nucleus of The Dormant Bulgeless Spiral Galaxy NGC 247

NGC 247 is a nearby late-type bulgeless spiral galaxy that contains an inactive nucleus. We report a serendipitous discovery of an X-ray flare from the galaxy center with a luminosity up to 2*10^39 erg/s in the 0.3-10 keV band with XMM-Newton. A Chandra observation confirms that the new X-ray source is spatially coincident with the galaxy nucleus. The XMM-Newton data revealed a hard power-law spectrum with a spectral break near 3-4 keV, no pulsations on timescales longer than 150 ms, and a flat power spectrum consistent with Poisson noise from 1 mHz to nearly 10 Hz. Follow-up observations with Swift detected a second flux peak followed by a luminosity drop by factor of almost 20. The spectral and temporal behaviors of the nuclear source are well consistent with the scenario that the flare was due to an outburst of a low-mass X-ray binary that contains a stellar-mass black hole emitting near its Eddington limit at the peak. However, it cannot be ruled out that the sudden brightening in the nucleus was due to accretion onto a possible low-mass nuclear black hole, fed by a tidally disrupted star or a gas cloud; the MAXI observations limit the peak luminosity of the flare to less than ~10^43 erg/s, suggesting that it is either a low mass black hole or an inefficient tidal disruption event (TDE).

preprint2015arXiv

A Revised Calibration of the Virial Mass Estimator for Black Holes in Active Galaxies Based on Single-epoch H$β$ Spectra

The masses of supermassive black holes in broad-line active galactic nuclei (AGNs) can be measured through reverberation mapping, but this method currently cannot be applied to very large samples or to high-redshift AGNs. As a practical alternative, one can devise empirical scaling relations, based on the correlation between broad-line region size and AGN luminosity and the relation between black hole mass and bulge stellar velocity dispersion, to estimate the virial masses of black holes from single-epoch spectroscopy. We present a revised calibration of the black hole mass estimator for the commonly used H$β$ emission line. Our new calibration takes into account the recent determination of the virial coefficient for pseudo and classical bulges.

preprint2015arXiv

An Off-Nucleus Nonstellar Black Hole in the Seyfert Galaxy NGC 5252

We report the discovery of a ultraluminous X-ray source (ULX; CXO J133815.6+043255) in NGC 5252. This ULX is an off-nuclear point-source, which is 22$^{\prime\prime}$ away from the center of NGC 5252, and has an X-ray luminosity of 1.5 $\times$ $10^{40}$erg s$^{-1}$. It is one of the rare examples of ULX, which exhibits clear counterparts in radio, optical, UV bands. Follow-up optical spectrum of the ULX shows strong emission lines. The redshift of [O III] emission line coincides with the systematic velocity of NGC 5252, suggesting the ULX is gravitationally bound to NGC 5252. The flux of [O III] appears to be correlated with both X-ray and radio luminosity in the same manner as ordinary AGNs, indicating that the [O III] emission is intrinsically associated with the ULX. Based on the multiwavelength data, we argue that the ULX is unlikely to be a background AGN. A more likely option is an accreting BH with a black hole mass of $\geq 10^4M_\odot$, which might be a stripped remnant of a merging dwarf galaxy.

preprint2015arXiv

An ultra-luminous quasar at z=5.363 with a ten billion solar mass black hole and a Metal-Rich DLA at z~5

We report the discovery of an ultra-luminous quasar J030642.51+185315.8 (hereafter J0306+1853) at redshift 5.363, which hosts a super-massive black hole (SMBH) with $M_{BH} = (1.07 \pm 0.27) \times10^{10}~M_\odot$. With an absolute magnitude $M_{1450}=-28.92$ and bolometric luminosity $L_{bol}\sim3.4\times10^{14} L_{\odot}$, J0306+1853 is one of the most luminous objects in the early Universe. It is not likely to be a beamed source based on its small flux variability, low radio loudness and normal broad emission lines. In addition, a $z=4.986$ Damped Ly$α$ system (DLA) with $\rm [M/H]=-1.3\pm0.1$, among the most metal rich DLAs at $z \gtrsim 5$, is detected in the absorption spectrum of this quasar. This ultra-luminous quasar puts strong constraint on the bright-end of quasar luminosity function and massive-end of black hole mass function. It will provide a unique laboratory to the study of BH growth and the co-evolution between BH and host galaxy with multi-wavelength follow-up observations. The future high resolution spectra will give more insights to the DLA and other absorption systems along the line-of-sight of J0306+1853.

preprint2015arXiv

Correlation between Galaxy Mergers and Luminous AGN

It is not yet clear what triggers the activity of active galactic nuclei (AGNs), but galaxy merging has been suspected to be one of the main mechanisms fuelling the activity. Using deep optical images taken at various ground-based telescopes, we investigate the fraction of galaxy mergers in 39 luminous AGNs (M$_{R}\, \lesssim$ -22.6 mag) at $z \leq$ 0.3 (a median redshift of 0.155), of which the host galaxies are generally considered as early-type galaxies. Through visual inspection of the images, we find that 17 of 39 AGN host galaxies (43.6%) show the evidence for current or past mergers like tidal tails, shells, and disturbed morphology. In order to see if this fraction is abnormally high, we also examined the merging fraction of normal early-type galaxies in the Sloan Digital Sky Survey (SDSS) Strip 82 data (a median redshift of 0.04), of which the surface-brightness limit is comparable to our imaging data. To correct for the effects related to the redshift difference of the two samples, we performed an image simulation by putting a bright point source as an artificial AGN in the images of SDSS early-type galaxies and placing them onto the redshifts of AGNs. The merging fraction in this realistic sample of simulated AGNs is only $\sim 5 - 15\%$ ($1/4$ to $1/8$ of that of real AGNs). Our result strongly suggests that luminous AGN activity is associated with galaxy merging.

preprint2015arXiv

Determining the Locations of Dust Sources in FeLoBAL Quasars

We conduct a spectroscopic search of quasars observed by the Sloan Digital Sky Survey (SDSS) with broad absorption line (BAL) troughs due to Mg II and troughs due to Fe II that simultaneously exhibit strong Balmer narrow emission lines (NELs). We find that in a redshift range of 0.4 less than or equal to z less than or equal to 0.9 approximately 23 of the 70 Mg II BALs and 4 of a subset of 15 Fe II BALs exhibit strong Balmer emission. We also find significant fractions of Mg II BALs (approximately 23%) and those Mg II BALs with Fe II troughs (approximately 27%) have strong continuum reddening, E(B - V) greater than or equal to 0.1. From measurements of the Balmer decrement in three objects, we find similarly significant reddening of the NEL region in three of the four objects; the NELs in the fourth object are not measurable. We also include one object in this study not taken from the SDSS sample that shows Fe II absorption and strong narrow emission, but due to measurement uncertainty and low continuum reddening the comparison is consistent but inconclusive. We find a trend in both the Mg II and Fe II BAL samples between the NEL reddening and continuum reddening. Because the narrow line reddening is consistent with the continuum reddening in every object in the two SDSS samples, it suggests that the reddening sources in these objects likely exist at larger radial distances than the narrow line regions from the central nucleus.

preprint2015arXiv

Differences between CO- and calcium triplet-derived velocity dispersions in spiral galaxies: evidence for central star formation?

We examine the stellar velocity dispersions (sigma) of a sample of 48 galaxies, 35 of which are spirals, from the Palomar nearby galaxy survey. It is known that for ultra-luminous infrared galaxies (ULIRGs) and merger remnants thesigma derived from the near-infrared CO band-heads is smaller than that measured from optical lines, while no discrepancy between these measurements is found for early-type galaxies. No such studies are available for spiral galaxies - the subject of this paper. We used cross-dispersed spectroscopic data obtained with the Gemini Near-Infrared Spectrograph (GNIRS), with spectral coverage from 0.85 to 2.5um, to obtain sigma measurements from the 2.29 $μ$m CO band-heads (sigma_{CO}), and the 0.85 um calcium triplet (sigma_{CaT}). For the spiral galaxies in the sample, we found that sigma_{CO} is smaller than sigma_{CaT}, with a mean fractional difference of 14.3%. The best fit to the data is given by sigma_{opt} = (46.0+/-18.1) + (0.85+/-0.12)sigma_{CO}. This "sigma discrepancy" may be related to the presence of warm dust, as suggested by a slight correlation between the discrepancy and the infrared luminosity. This is consistent with studies that have found no sigma-discrepancy in dust-poor early-type galaxies, and a much larger discrepancy in dusty merger remnants and ULIRGs. That sigma_{CO}$ is lower than sigma_{opt} may also indicate the presence of a dynamically cold stellar population component. This would agree with the spatial correspondence between low sigma_{CO} and young/intermediate-age stellar populations that has been observed in spatially-resolved spectroscopy of a handful of galaxies.

preprint2015arXiv

Dissecting the Power Sources of Low-Luminosity Emission-Line Galaxy Nuclei via Comparison of HST-STIS and Ground-Based Spectra

Using a sample of ~100 nearby line-emitting galaxy nuclei, we have built the currently definitive atlas of spectroscopic measurements of H_alpha and neighboring emission lines at subarcsecond scales. We employ these data in a quantitative comparison of the nebular emission in Hubble Space Telescope (HST) and ground-based apertures, which offer an order-of-magnitude difference in contrast, and provide new statistical constraints on the degree to which Transition Objects and low-ionization nuclear emission-line regions (LINERs) are powered by an accreting black hole at <10 pc. We show that while the small-aperture observations clearly resolve the nebular emission, the aperture dependence in the line ratios is generally weak, and this can be explained by gradients in the density of the line-emitting gas: the higher densities in the more nuclear regions potentially flatten the excitation gradients, suppressing the forbidden emission. The Transition Objects show a threefold increase in the incidence of broad H_alpha emission in the high-resolution data, as well as the strongest density gradients, supporting the composite model for these systems as accreting sources surrounded by star-forming activity. The narrow-line LINERs appear to be the weaker counterparts of the Type 1 LINERs, where the low accretion rates cause the disappearance of the broad-line component. The enhanced sensitivity of the HST observations reveals a 30% increase in the incidence of accretion-powered systems at z~0. A comparison of the strength of the broad-line emission detected at different epochs implies potential broad-line variability on a decade-long timescale, with at least a factor of three in amplitude.

preprint2015arXiv

Globular Cluster Populations: Results Including S$^4$G Late-Type Galaxies

Using 3.6 and 4.5$μ$m images of 73 late-type, edge-on galaxies from the S$^4$G survey, we compare the richness of the globular cluster populations of these galaxies to those of early type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M$_*$, specifically $\log({\rm M}_*/{\rm M}_\odot) < 10$, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per $10^9\ {\rm M}_\odot$ of galaxy stellar mass, T$_{\rm N}$, on M$_*$. For $8.5 < \log ({\rm M}_*/{\rm M}_\odot) < 10.5$ we find the relationship can be satisfactorily described as T$_{\rm N} = ({\rm M}_*/10^{6.7})^{-0.56}$ when M$_*$ is expressed in solar masses. The functional form of the relationship is only weakly constrained and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early-types, do not show the tendency for low mass galaxies to split into two T$_{\rm N}$ families. Using these results and a galaxy stellar mass function from the literature, we calculate that in a volume limited, local Universe sample, clusters are most likely to be found around fairly massive galaxies (M$_* \sim 10^{10.8}$ M$_\odot$) and present a fitting function for the volume number density of clusters as a function of parent galaxy stellar mass. We find no correlation between T$_{\rm N}$ and large-scale environment, but do find a tendency for galaxies of fixed M$_*$ to have larger T$_{\rm N}$ if they have converted a larger proportion of their baryons into stars.

preprint2015arXiv

H-alpha kinematics of S4G spiral galaxies-II. Data description and non-circular motions

We present a kinematical study of 29 spiral galaxies included in the Spitzer Survey of Stellar Structure in Galaxies, using Halpha Fabry-Perot data obtained with the Galaxy Halpha Fabry-Perot System instrument at the William Herschel Telescope in La Palma, complemented with images in the R-band and in Halpha. The primary goal is to study the evolution and properties of the main structural components of galaxies through the kinematical analysis of the FP data, complemented with studies of morphology, star formation and mass distribution. In this paper we describe how the FP data have been obtained, processed and analysed. We present the resulting moment maps, rotation curves, velocity model maps and residual maps. Images are available in FITS format through the NASA/IPAC Extragalactic Database and the Centre de Données Stellaires. With these data products we study the non-circular motions, in particular those found along the bars and spiral arms. The data indicate that the amplitude of the non-circular motions created by the bar does not correlate with the bar strength indicators. The amplitude of those non-circular motions in the spiral arms does not correlate with either arm class or star formation rate along the spiral arms. This implies that the presence and the magnitude of the streaming motions in the arms is a local phenomenon.

preprint2015arXiv

Herschel Survey of the Palomar-Green QSOs at Low Redshift

We investigate the global cold dust properties of 85 nearby (z < 0.5) QSOs, chosen from the Palomar-Green sample of optically luminous quasars. We determine their infrared spectral energy distributions and estimate their rest-frame luminosities by combining Herschel data from 70 to 500 microns with near-infrared and mid-infrared measurements from the Two Micron All Sky Survey (2MASS) and the Wide-Field Infrared Survey Explorer (WISE). In most sources the far-infrared (FIR) emission can be attributed to thermally heated dust. Single temperature modified black body fits to the FIR photometry give an average dust temperature for the sample of 33~K, with a standard deviation of 8~K, and an average dust mass of 7E6 Solar Masses with a standard deviation of 9E6 Solar Masses. Estimates of star-formation that are based on the FIR continuum emission correlate with those based on the 11.3 microns PAH feature, however, the star-formation rates estimated from the FIR continuum are higher than those estimated from the 11.3 microns PAH emission. We attribute this result to a variety of factors including the possible destruction of the PAHs and that, in some sources, a fraction of the FIR originates from dust heated by the active galactic nucleus and by old stars.

preprint2015arXiv

Measuring the mass of the central black hole in the bulgeless galaxy NGC 4395 from gas dynamical modeling

NGC 4395 is a bulgeless spiral galaxy, harboring one of the nearest known type 1 Seyfert nuclei. Although there is no consensus on the mass of its central engine, several estimates suggest it to be one of the lightest massive black holes (MBHs) known. We present the first direct dynamical measurement of the mass of this MBH from a combination of two-dimensional gas kinematic data, obtained with the adaptive optics assisted near infrared integral field spectrograph Gemini/NIFS, and high-resolution multiband photometric data from Hubble Space Telescope's Wide Field Camera 3 (HST/WFC3). We use the photometric data to model the shape and stellar mass-to-light ratio (M/L) of the nuclear star cluster. From the Gemini/NIFS observations, we derive the kinematics of warm molecular hydrogen gas as traced by emission through the H$_2$ 1--0 S(1) transition. These kinematics show a clear rotational signal, with a position angle orthogonal to NGC 4395's radio jet. Our best fitting tilted ring models of the kinematics of the molecular hydrogen gas contain a black hole with mass $M=4_{-3}^{+8}\times 10^5$ M$_\odot$ (3$σ$ uncertainties) embedded in a nuclear star cluster of mass $M=2 \times 10^6$ M$_\odot$. Our black hole mass measurement is in excellent agreement with the reverberation mapping mass estimate of Peterson et al. (2005), but shows some tension with other mass measurement methods based on accretion signals.

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

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. III. Detection of Fe II Reverberation in Nine Narrow-Line Seyfert 1 Galaxies

This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 2012--2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six of the newly detected lags are indistinguishable from the Hbeta lags measured in the same sources. Two are significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we find a Fe II radius--luminosity relationship similar to the one for Hbeta, although our sample by itself shows no clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for the majority of the newly reported objects, is indistinguishable from the Hbeta-emitting gas. We also present a tentative correlation between the lag and intensity of Fe II and Hbeta and comment on its possible origin.

preprint2015arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. IV. H$β$ Time Lags and Implications for Super-Eddington Accretion

We have completed two years of photometric and spectroscopic monitoring of a large number of active galactic nuclei (AGNs) with very high accretion rates. In this paper, we report on the result of the second phase of the campaign, during 2013--2014, and the measurements of five new H$β$ time lags out of eight monitored AGNs. All five objects were identified as super-Eddington accreting massive black holes (SEAMBHs). The highest measured accretion rates for the objects in this campaign are $\dot{\mathscr{M}}\gtrsim 200$, where $\dot{\mathscr{M}}= \dot{M}_{\bullet}/L_{\rm Edd}c^{-2}$, $\dot{M}_{\bullet}$ is the mass accretion rates, $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. We find that the H$β$ time lags in SEAMBHs are significantly shorter than those measured in sub-Eddington AGNs, and the deviations increase with increasing accretion rates. Thus, the relationship between broad-line region size ($R_{_{\rm Hβ}}$) and optical luminosity at 5100Å, $R_{_{\rm Hβ}}-L_{5100}$, requires accretion rate as an additional parameter. We propose that much of the effect may be due to the strong anisotropy of the emitted slim-disk radiation. Scaling $R_{_{\rm Hβ}}$ by the gravitational radius of the black hole, we define a new radius-mass parameter ($Y$) and show that it saturates at a critical accretion rate of $\dot{\mathscr{M}}_c=6\sim 30$, indicating a transition from thin to slim accretion disk and a saturated luminosity of the slim disks. The parameter $Y$ is a very useful probe for understanding the various types of accretion onto massive black holes. We briefly comment on implications to the general population of super-Eddington AGNs in the universe and applications to cosmology.

preprint2015arXiv

The Odd Offset between the Galactic Disk and Its Bar in NGC 3906

We use mid-infrared 3.6 and 4.5microns imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ~720 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter sub structure in the galaxy's halo or from a recent interaction with a fast moving neighbor which remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar / kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Ensemble Spectroscopic Variability of Quasar Broad Emission Lines

We explore the variability of quasars in the MgII and Hbeta broad emission lines and UV/optical continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic epochs from SDSS-RM over $6$ months, containing 357 quasars with MgII and 41 quasars with Hbeta . On longer timescales, the study is also supplemented with two-epoch data from SDSS-I/II. The SDSS-I/II data include an additional $2854$ quasars with MgII and 572 quasars with Hbeta. The MgII emission line is significantly variable ($Δf/f$ 10% on 100-day timescales), a necessary prerequisite for its use for reverberation mapping studies. The data also confirm that continuum variability increases with timescale and decreases with luminosity, and the continuum light curves are consistent with a damped random-walk model on rest-frame timescales of $\gtrsim 5$ days. We compare the emission-line and continuum variability to investigate the structure of the broad-line region. Broad-line variability shows a shallower increase with timescale compared to the continuum emission, demonstrating that the broad-line transfer function is not a $δ$-function. Hbeta is more variable than MgII (roughly by a factor of $1.5$), suggesting different excitation mechanisms, optical depths and/or geometrical configuration for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM project have important consequences for future studies of reverberation mapping and black hole mass estimation of $1<z<2$ quasars.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: First Broad-line Hbeta and MgII Lags at z>~0.3 from six-Month Spectroscopy

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: No Evidence for Evolution in the M-sigma Relation to z~1

We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1<z<1 (46 at z>0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal-to-noise ratio coadded spectra (average S/N~30 per 69 km/s pixel) from SDSS-RM allowed decomposition of the host and quasar spectra, and measurement of the host stellar velocity dispersions and black hole (BH) masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity (L5100=10^(43.2-44.7) erg/s) lead to the first clear detection of a correlation between SE virial BH mass and host stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly flatter than the local relation, suggesting that there are selection biases in high-z luminosity-threshold quasar samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution of the M-sigma relation free of caveats by comparing different samples/analyses at disjoint redshifts. We do not observe evolution of the M-sigma relation in our sample, up to z~1, but there is an indication that the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that a constant M-sigma relation is favored to z~1. Our results highlight the scientific potential of deep coadded spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and galaxies at earlier times.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Post-Starburst Signatures in Quasar Host Galaxies at z < 1

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years, which was subsequently quenched or suppressed. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.

preprint2015arXiv

The Stellar Spectral Features of Nearby Galaxies in the Near-Infrared: Tracers of Thermally-Pulsing Asymptotic Giant Branch Stars?

We analyze the stellar absorption features in high signal-to-noise ratio near-infrared (NIR) spectra of the nuclear region of 12 nearby galaxies, mostly spirals. The features detected in some or all of the galaxies in this sample are the TiO (0.843 $μ$m\ and 0.886 $μ$m), VO (1.048 $μ$m), CN (1.1 $μ$m\ and 1.4 $μ$m), H$\rm _2$O (1.4 $μ$m\ and 1.9 $μ$m) and CO (1.6 $μ$m\ and 2.3 $μ$m) bands. The C$\rm _2$ (1.17 $μ$m\ and 1.76 $μ$m) bands are generally weak or absent, although C$\rm _2$ (1.76 $μ$m) may be weakly present in the mean galaxy spectrum. A deep feature near 0.93 $μ$m, likely caused by CN, TiO and/or ZrO, is also detected in all objects. Fitting a combination of stellar spectra to the mean spectrum shows that the absorption features are produced by evolved stars: cool giants and supergiant stars in the early- or thermally-pulsing asymptotic giant branch (E-AGB or TP-AGB) phases. The high luminosity of TP-AGB stars, and the appearance of VO and ZrO features in the data, suggest that TP-AGB stars dominate these spectral features. However, a contribution from other evolved stars is also likely. Comparison with evolutionary population synthesis models shows that models based on empirical libraries that predict relatively strong NIR features provide a more accurate description of the data. However, none of the models tested accurately reproduces all of the features observed in the spectra. To do so, the models will need to not only improve the treatment of TP-AGB stars, but also include good quality spectra of red giant and E-AGB stars. The uninterrupted wavelength coverage, high S/N, and quantity of features we present here will provide a benchmark for the next generation of models aiming to explain and predict the NIR properties of galaxies.

preprint2015arXiv

The Structure of Nuclear Star Clusters in Nearby Late-type Spiral Galaxies from Hubble Space Telescope Wide Field Camera 3 Imaging

We obtained Hubble Space Telescope/Wide Field Camera 3 imaging of a sample of ten of the nearest and brightest nuclear clusters residing in late-type spiral galaxies, in seven bands that span the near-ultraviolet to the near-infrared. Structural properties of the clusters were measured by fitting two-dimensional surface brightness profiles to the images using GALFIT. The clusters exhibit a wide range of structural properties. For six of the ten clusters in our sample, we find changes in the effective radius with wavelength, suggesting radially varying stellar populations. In four of the objects, the effective radius increases with wavelength, indicating the presence of a younger population which is more concentrated than the bulk of the stars in the cluster. However, we find a general decrease in effective radius with wavelength in two of the objects in our sample, which may indicate extended, circumnuclear star formation. We also find a general trend of increasing roundness of the clusters at longer wavelengths, as well as a correlation between the axis ratios of the NCs and their host galaxies. These observations indicate that blue disks aligned with the host galaxy plane are a common feature of nuclear clusters in late-type galaxies, but are difficult to detect in galaxies that are close to face-on. In color-color diagrams spanning the near-UV through the near-IR, most of the clusters lie far from single-burst evolutionary tracks, showing evidence for multi-age populations. Most of the clusters have integrated colors consistent with a mix of an old population (> 1 Gyr) and a young population (~100-300 Myr). The wide wavelength coverage of our data provides a sensitivity to populations with a mix of ages that would not be possible to achieve with imaging in optical bands only.

preprint2015arXiv

Tomography of X-ray Nova Muscae 1991: Evidence for ongoing mass transfer and stream-disc overflow

We present a spectroscopic analysis of the black hole binary Nova Muscae 1991 in quiescence using data obtained in 2009 with MagE on the Magellan Clay telescope and in 2010 with IMACS on the Magellan Baade telescope at the Las Campanas Observatory. Emission from the disc is observed in H alpha, H beta and Ca II (8662 A). A prominent hotspot is observed in the Doppler maps of all three emission lines. The existence of this spot establishes ongoing mass transfer from the donor star in 2009-2010 and, given its absence in the 1993-1995 observations, demonstrates the presence of a variable hotspot in the system. We find the radial distance to the hotspot from the black hole to be consistent with the circularization radius. Our tomograms are suggestive of stream-disc overflow in the system. We also detect possible Ca II (8662 A) absorption from the donor star.

preprint2014arXiv

A Bayesian Approach to Estimate the Size and Structure of the Broad-line Region in Active Galactic Nuclei Using Reverberation Mapping Data

This is the first paper in a series devoted to systematic study of the size and structure of the broad-line region (BLR) in active galactic nuclei (AGNs) using reverberation mapping (RM) data. We employ a recently developed Bayesian approach that statistically describes the variabibility as a damped random walk process and delineates the BLR structure using a flexible disk geometry that can account for a variety of shapes, including disks, rings, shells, and spheres. We allow for the possibility that the line emission may respond non-linearly to the continuum, and we detrend the light curves when there is clear evidence for secular variation. We use a Markov Chain Monte Carlo implementation based on Bayesian statistics to recover the parameters and uncertainties for the BLR model. The corresponding transfer function is obtained self-consistently. We tentatively constrain the virial factor used to estimate black hole masses; more accurate determinations will have to await velocity-resolved RM data. Application of our method to RM data with Hbeta monitoring for about 40 objects shows that the assumed BLR geometry can reproduce quite well the observed emission-line fluxes from the continuum light curves. We find that the Hbeta BLR sizes obtained from our method are on average ~20% larger than those derived from the traditional cross-correlation method. Nevertheless, we still find a tight BLR size-luminosity relation with a slope of alpha=0.55\pm0.03 and an intrinsic scatter of ~0.18 dex. In particular, we demonstrate that our approach yields appropriate BLR sizes for some objects (such as Mrk 142 and PG 2130+099) where traditional methods previously encountered difficulties.

preprint2014arXiv

A New Approach to Constrain Black Hole Spins in Active Galaxies Using Optical Reverberation Mapping

A tight relation between the size of the broad-line region (BLR) and optical luminosity has been established in about 50 active galactic nuclei studied through reverberation mapping of the broad Hbeta emission line. The R_blr-L relation arises from simple photoionization considerations. Using a general relativistic model of an optically thick, geometrically thin accretion disk, we show that the ionizing luminosity jointly depends on black hole mass, accretion rate, and spin. The non-monotonic relation between the ionizing and optical luminosity gives rise to a complicated relation between the BLR size and the optical luminosity. We show that the reverberation lag of Hbeta to the varying continuum depends very sensitively to black hole spin. For retrograde spins, the disk is so cold that there is a deficit of ionizing photons in the BLR, resulting in shrinkage of the hydrogen ionization front with increasing optical luminosity, and hence shortened Hbeta lags. This effect is specially striking for luminous quasars undergoing retrograde accretion, manifesting in strong deviations from the canonical R_blr-L relation. This could lead to a method to estimate black hole spins of quasars and to study their cosmic evolution. At the same time, the small scatter of the observed R_blr-L relation for the current sample of reverberation-mapped active galaxies implies that the majority of these sources have rapidly spinning black holes.

preprint2014arXiv

Correlation between the photon index and X-ray luminosity of black hole X-ray binaries and active galactic nuclei: observations and interpretation

We investigate the observed correlation between the 2--10 keV X-ray luminosity (in unit of the Eddington luminosity; $l_X \equiv L_X/L_{Edd}$) and the photon index ($Γ$) of the X-ray spectrum for both black hole X-ray binaries (BHBs) and active galactic nuclei (AGNs). We construct a large sample, with $10^{-9} < l_X < 10^{-1}$. We find that $Γ$ is positively and negatively correlated with $l_X$ when $l_X > 10^{-3}$ and $10^{-6.5} < l_X < 10^{-3}$ respectively, while $Γ$ is nearly a constant when $l_X < 10^{-6.5}$. We explain the above correlation in the framework of a coupled hot accretion flow -- jet model. The radio emission always come from the jet while the X-ray emission comes from the accretion flow and jet when $l_X$ is above and below $10^{-6.5}$, respectively. More specifically, we assume that with the increase of mass accretion rate, the hot accretion flow develops into a clumpy and further a disc -- corona two-phase structure because of thermal instability. We argue that such kind of two-phase accretion flow can explain the observed positive correlation.

preprint2014arXiv

Globular Cluster Populations: First Results from S$^4$G Early-Type Galaxies

Using 3.6$μ$m images of 97 early-type galaxies, we develop and verify methodology to measure globular cluster populations from the S$^4$G survey images. We find that 1) the ratio, T$_{\rm N}$, of the number of clusters, N$_{\rm CL}$, to parent galaxy stellar mass, M$_*$, rises weakly with M$_*$ for early-type galaxies with M$_* > 10^{10}$ M$_\odot$ when we calculate galaxy masses using a universal stellar initial mass function (IMF), but that the dependence of T$_{\rm N}$ on M$_*$ is removed entirely once we correct for the recently uncovered systematic variation of IMF with M$_*$, and 2) for M$_* < 10^{10}$ M$_\odot$ there is no trend between N$_{\rm CL}$ and M$_*$, the scatter in T$_{\rm N}$ is significantly larger (approaching 2 orders of magnitude), and there is evidence to support a previous, independent suggestion of two families of galaxies. The behavior of N$_{\rm CL}$ in the lower mass systems is more difficult to measure because these systems are inherently cluster poor, but our results may add to previous evidence that large variations in cluster formation and destruction efficiencies are to be found among low mass galaxies. The average fraction of stellar mass in clusters is $\sim$ 0.0014 for M$_* > 10^{10}$ M$_\odot$ and can be as large as $\sim 0.02$ for less massive galaxies. These are the first results from the S$^4$G sample of galaxies, and will be enhanced by the sample of early-type galaxies now being added to S$^4$G and complemented by the study of later type galaxies within S$^4$G.

preprint2014arXiv

How Robust Are the Size Measurements of High-redshift Compact Galaxies?

Massive quiescent galaxies at $z \approx 2$ are apparently much more compact than galaxies of comparable mass today. How robust are these size measurements? We perform comprehensive simulations to determine possible biases and uncertainties in fitting single-component light distributions to real galaxies. In particular, we examine the robustness of the measurements of the luminosity, size, and other structural parameters. We devise simulations with increasing realism to systematically disentangle effects due to the technique (specifically using GALFIT) and the intrinsic structures of the galaxies. By accurately capturing the detailed substructures of nearby elliptical galaxies and then rescaling their sizes and signal-to-noise to mimic galaxies at different redshifts, we confirm that the massive quiescent galaxies at $z \approx 2$ are significantly more compact intrinsically than their local counterparts. Their observed compactness is not a result of missing faint outer light due to systematic errors in modeling. In fact, we find that fitting multi-component galaxies with a single Sérsic profile, the procedure most commonly adopted in the literature, biases the inferred sizes higher by up to 10% - 20%, which accentuates the amount of size evolution required. If the sky estimation has been done robustly and the model for the point-spread function is fairly accurate, GALFIT can retrieve the properties of single-component galaxies over a wide range of signal-to-noise ratios without introducing any systematic errors.

preprint2014arXiv

Inter-comparison of Radio-Loudness Criteria for Type 1 AGNs in the XMM-COSMOS Survey

Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24μm}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500Å})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(νL_ν(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24μm}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria.

preprint2014arXiv

IRAC Near-Infrared Features in the Outer Parts of S4G Galaxies

We present a catalogue and images of visually detected features, such as asymmetries, extensions, warps, shells, tidal tails, polar rings, and obvious signs of mergers or interactions, in the faint outer regions (at and outside of R_25) of nearby galaxies. This catalogue can be used in future quantitative studies that examine galaxy evolution due to internal and external factors. We are able to reliably detect outer region features down to a brightness level of 0.03 MJy/sr per pixel at 3.6 microns in the Spitzer Survey of Stellar Structure in Galaxies (S4G). We also tabulate companion galaxies. We find asymmetries in the outer isophotes in 22+/-1 per cent of the sample. The asymmetry fraction does not correlate with galaxy classification as an interacting galaxy or merger remnant, or with the presence of companions. We also compare the detected features to similar features in galaxies taken from cosmological zoom re-simulations. The simulated images have a higher fraction (33 per cent) of outer disc asymmetries, which may be due to selection effects and an uncertain star formation threshold in the models. The asymmetries may have either an internal (e.g., lopsidedness due to dark halo asymmetry) or external origin.

preprint2014arXiv

Medium Resolution Near-Infrared Spectra of the Host Galaxies of Nearby Quasars

We present medium resolution near-infrared host galaxy spectra of low redshift quasars, PG 0844 + 349 (z=0.064), PG 1226 + 023 (z=0.158), and PG 1426+015 (z=0.086). The observations were done by using the Infrared Camera and Spectrograph (IRCS) at the Subaru 8.2 m telescope. The full width at half maximum of the point spread function was about 0.3 arcsec by operations of an adaptive optics system, which can effectively resolve the quasar spectra from the host galaxy spectra. We spent up to several hours per target and developed data reduction methods to reduce the systematic noises of the telluric emissions and absorptions. From the obtained spectra, we identified absorption features of Mg I (1.503 um), Si I (1.589 um) and CO (6-3) (1.619 um), and measured the velocity dispersions of PG 0844 + 349 to be 132+/-110 km s-1 and PG 1426 + 015 to be 264+/-215 km s-1. By using an M_BH-sigma relation of elliptical galaxies, we derived the black hole (BH) mass of PG 0844+349, log(M_BH/M_SUN) = 7.7+/-5.5 and PG 1426+015, log(M_BH/M_SUN) = 9.0+/-7.5. These values are consistent with the BH mass values from broad emission lines with an assumption of a virial factor of 5.5.

preprint2014arXiv

Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: II. The conversion from light to mass

We present a new approach for estimating the 3.6 micron stellar mass-to-light ratio in terms of the [3.6]-[4.5] colors of old stellar populations. Our approach avoids several of the largest sources of uncertainty in existing techniques. By focusing on mid-IR wavelengths, we gain a virtually dust extinction-free tracer of the old stars, avoiding the need to adopt a dust model to correctly interpret optical or optical/NIR colors normally leveraged to assign M/L. By calibrating a new relation between NIR and mid-IR colors of GLIMPSE giant stars we also avoid discrepancies in model predictions for the [3.6]-[4.5] colors of old stellar populations due to uncertainties in molecular line opacities. We find that the [3.6]-[4.5] color, which is driven primarily by metallicity, provides a tight constraint on M/L_3.6, which varies intrinsically less than at optical wavelengths. The uncertainty on M/L_3.6 of ~0.07 dex due to unconstrained age variations marks a significant improvement on existing techniques for estimating the stellar M/L with shorter wavelength data. A single M/L_3.6=0.6 (assuming a Chabrier IMF), independent of [3.6]-[4.5] color, is also feasible as it can be applied simultaneously to old, metal-rich and young, metal-poor populations, and still with comparable (or better) accuracy (~0.1 dex) as alternatives. We expect our M/L_3.6 to be optimal for mapping the stellar mass distributions in S4G galaxies, for which we have developed an Independent Component Analysis technique to first isolate the old stellar light at 3.6 micron from non-stellar emission (e.g. hot dust and the 3.3 PAH feature). Our estimate can also be used to determine the fractional contribution of non-stellar emission to global (rest-frame) 3.6 micron fluxes, e.g. in WISE imaging, and establishes a reliable basis for exploring variations in the stellar IMF.

preprint2014arXiv

SDSS J013127.34$-$032100.1: A newly discovered radio-loud quasar at $z=5.18$ with extremely high luminosity

Only very few z>5 quasars discovered to date are radio-loud, with a radio-to-optical flux ratio (radio-loudness parameter) higher than 10. Here we report the discovery of an optically luminous radio-loud quasar, SDSS J013127.34-032100.1 (J0131-0321 in short), at z=5.18+-0.01 using the Lijiang 2.4m and Magellan telescopes. J0131-0321 has a spectral energy distribution consistent with that of radio-loud quasars. With an i-band magnitude of 18.47 and radio flux density of 33 mJy, its radio-loudness parameter is ~100. The optical and near-infrared spectra taken by Magellan enable us to estimate its bolometric luminosity to be L_bol ~ 1.1E48 erg/s, approximately 4.5 times greater than that of the most distant quasar known to date. The black hole mass of J0131-0321 is estimated to be 2.7E9 solar masses, with an uncertainty up to 0.4 dex. Detailed physical properties of this high-redshift, radio-loud, potentially super-Eddington quasar can be probed in the future with more dedicated and intensive follow-up observations using multi-wavelength facilities.

preprint2014arXiv

The Baryonic Tully-Fisher Relationship for S$^4$G Galaxies and the "Condensed" Baryon Fraction of Galaxies

We combine data from the Spitzer Survey for Stellar Structure in Galaxies (S$^4$G), a recently calibrated empirical stellar mass estimator from Eskew et al., and an extensive database of HI spectral line profiles to examine the baryonic Tully-Fisher (BTF) relation. We find 1) that the BTF has lower scatter than the classic Tully-Fisher (TF) relation and is better described as a linear relationship, confirming similar previous results, 2) that the inclusion of a radial scale in the BTF decreases the scatter but only modestly, as seen previously for the TF relation, and 3) that the slope of the BTF, which we find to be $3.5\pm 0.2$ ($Δ$ log $M_{baryon}/Δ$ log $v_c$), implies that on average a nearly constant fraction ($\sim 0.4$) of all baryons expected to be in a halo are "condensed" onto the central region of rotationally supported galaxies. The condensed baryon fraction, $M_{baryon}/M_{total}$, is, to our measurement precision, nearly independent of galaxy circular velocity (our sample spans circular velocities, $v_c$, between 60 and 250 km s$^{-1}$, but is extended to $v_c\sim 10$ km s$^{-1}$ using data from the literature). The observed galaxy-to-galaxy scatter in this fraction is generally $\le$ a factor of 2 despite fairly liberal selection criteria. These results imply that cooling and heating processes, such as cold vs. hot accretion, mass loss due to stellar winds, and AGN driven feedback, to the degree that they affect the global galactic properties involved in the BTF, are independent of halo mass for galaxies with $10 < v_c < 250$ km/s and typically introduce no more than a factor of two range in the resulting $M_{baryon}/M_{total}$. Recent simulations by Aumer et al. of a small sample of disk galaxies are in excellent agreement with our data, suggesting that current simulations are capable of reproducing the global properties of individual disk galaxies.

preprint2014arXiv

The Black Hole Mass Scale of Classical and Pseudo Bulges in Active Galaxies

The mass estimator used to calculate black hole (BH) masses in broad-line active galactic nuclei (AGNs) relies on a virial coefficient (the "$f$ factor") that is determined by comparing reverberation-mapped (RM) AGNs with measured bulge stellar velocity dispersions against the $M_{\rm BH}-σ_*$ relation of inactive galaxies. It has recently been recognized that only classical bulges and ellipticals obey a tight $M_{\rm BH}-σ_*$ relation; pseudobulges have a different zero point and much larger scatter. Motivated by these developments, we reevaluate the $f$ factor for RM AGNs with available $σ_*$ measurements, updated H$β$ RM lags, and new bulge classifications based on detailed decomposition of high-resolution ground-based and space-based images. Separate calibrations are provided for the two bulge types, whose virial coefficients differ by a factor of $\sim 2$: $f=6.3\pm1.5$ for classical bulges and ellipticals and $f = 3.2\pm0.7$ for pseudobulges. The structure and kinematics of the broad-line region, at least as crudely encoded in the $f$ factor, seems to related to the large-scale properties or formation history of the bulge. Lastly, we investigate the bulge stellar masses of the RM AGNs, show evidence for recent star formation in the AGN hosts that correlates with Eddington ratio, and discuss the potential utility of the $M_{\rm BH}-M_{\rm bulge}$ relation as a more promising alternative to the conventionally used $M_{\rm BH}-σ_*$ relation for future refinement of the virial mass estimator for AGNs.

preprint2014arXiv

The diversity of quasars unified by accretion and orientation

Quasars are rapidly accreting supermassive black holes at the center of massive galaxies. They display a broad range of properties across all wavelengths, reflecting the diversity in the physical conditions of the regions close to the central engine. These properties, however, are not random, but form well-defined trends. The dominant trend is known as Eigenvector 1, where many properties correlate with the strength of optical iron and [OIII] emission. The main physical driver of Eigenvector 1 has long been suspected to be the quasar luminosity normalized by the mass of the hole (the Eddington ratio), an important quantity of the black hole accretion process. But a definitive proof has been missing. Here we report an analysis of archival data that reveals that Eddington ratio indeed drives Eigenvector 1. We also find that orientation plays a significant role in determining the observed kinematics of the gas, implying a flattened, disklike geometry for the fast-moving clouds close to the hole. Our results show that most of the diversity of quasar phenomenology can be unified with two simple quantities, Eddington ratio and orientation.

preprint2014arXiv

The Mass Profile and Shape of Bars in the Spitzer Survey of Stellar Structure in Galaxies (S4G): Search for an Age Indicator for Bars

We have measured the radial light profiles and global shapes of bars using two-dimensional 3.6 $μm $ image decompositions for 144 face-on barred galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G). The bar surface brightness profile is correlated with the stellar mass and bulge-to-total (B/T) ratio of their host galaxies. Bars in massive and bulge-dominated galaxies (B/T$>$0.2) show a flat profile, while bars in less massive, disk-dominated galaxies (B/T$\sim$0) show an exponential, disk-like profile with a wider spread in the radial profile than in the bulge-dominated galaxies. The global two-dimensional shapes of bars, however, are rectangular/boxy, independent of the bulge or disk properties. We speculate that because bars are formed out of disk, bars initially have an exponential (disk-like) profile which evolves over time, trapping more stars into the boxy bar orbits. This leads bars to become stronger and have flatter profiles. The narrow spread of bar radial profiles in more massive disks suggests that these bars formed earlier (z$>$1), while the disk-like profiles and a larger spread in the radial profile in less massive systems imply a later and more gradual evolution, consistent with the cosmological evolution of bars inferred from observational studies. Therefore, we expect that the flatness of the bar profile can be used as a dynamical age indicator of the bar to measure the time elapsed since the bar formation. We argue that cosmic gas accretion is required to explain our results on bar profile and the presence of gas within the bar region.

preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

preprint2013arXiv

A comparative analysis of virial black-hole mass estimates of moderate-luminosity active galactic nuclei using Subaru/FMOS

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by <log(L_3000/L_Halpha)> = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2.

preprint2013arXiv

A Novel Approach to Constrain the Mass Ratio of Minor Mergers in Elliptical Galaxies: Application to NGC 4889, the Brightest Cluster Galaxy in Coma

Minor mergers are thought to be important for the build-up and structural evolution of massive elliptical galaxies. In this work, we report the discovery of a system of four shell features in NGC 4889, one of the brightest members of the Coma cluster, using optical images taken with the Hubble Space Telescope and the Sloan Digital Sky Survey. The shells are well aligned with the major axis of the host and are likely to have been formed by the accretion of a small satellite galaxy. We have performed a detailed two-dimensional photometric decomposition of NGC 4889 and of the many overlapping nearby galaxies in its vicinity. This comprehensive model allows us not only to firmly detect the low-surface brightness shells, but, crucially, also to accurately measure their luminosities and colors. The shells are bluer than the underlying stars at the same radius in the main galaxy. We make use of the colors of the shells and the color-magnitude relation of the Coma cluster to infer the luminosity (or mass) of the progenitor galaxy. The shells in NGC 4889 appear to have been produced by the minor merger of a moderate-luminosity (M_I ~ -18.7 mag) disk (S0 or spiral) galaxy with a luminosity (mass) ratio of ~ 90:1 with respect to the primary galaxy. The novel methodology presented in this work can be exploited to decode the fossil record imprinted in the photometric substructure of other nearby early-type galaxies.

preprint2013arXiv

A Physical Link Between Jet Formation and Hot Plasma in Active Galactic Nuclei

Recent observations suggest that in black hole X-ray binaries jet/outflow formation is related to the hot plasma in the vicinity of the black hole, either in the form of an advection-dominated accretion flow at low accretion rates or in a disk corona at high accretion rates. We test the viability of this scenario for supermassive black holes using two samples of active galactic nuclei distinguished by the presence (radio-strong) and absence (radio-weak) of well-collimated, relativistic jets. Each is centered on a narrow range of black hole mass but spans a very broad range of Eddington ratios, effectively simulating, in a statistical manner, the behavior of a single black hole evolving across a wide spread in accretion states. Unlike the relationship between the radio and optical luminosity, which shows an abruptly break between high- and low-luminosity sources at an Eddington ratio of ~1, the radio emission-a measure of the jet power-varies continuously with the hard X-ray (2-10 keV) luminosity, roughly as L_R \propto L_X^(0.6-0.75). This relation, which holds for both radio-weak and radio-strong active galaxies, is similar to the one seen in X-ray binaries. Jet/outflow formation appears to be closely linked to the conditions that give rise to the hot, optically thin coronal emission associated with accretion flows, both in the regime of low and high accretion rates.

preprint2013arXiv

A Quasar-Galaxy Mixing Diagram: Quasar Spectral Energy Distribution Shapes in the Optical to Near-Infrared

We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1μm to 3000Å and from 1μm to 3μm in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the mixing diagram, we find that a combination of the Elvis et al. (1994, hereafter E94) quasar SED with various contributions from galaxy emission and some dust reddening is remarkably effective in describing the SED shape from 0.3-3μm for large ranges of redshift, luminosity, black hole mass and Eddington ratio of type 1 AGN. In particular, the location in the mixing diagram of the highest luminosity AGN is very close (within 1σ) to that of the E94 SED. The mixing diagram can also be used to estimate the host galaxy fraction and reddening in quasar. We also show examples of some outliers which might be AGN in different evolutionary stages compared to the majority of AGN in the quasar-host galaxy co-evolution cycle.

preprint2013arXiv

Coevolution (Or Not) of Supermassive Black Holes and Host Galaxies

We review the observed demographics and inferred evolution of supermassive black holes (BHs) found by dynamical modeling of spatially resolved kinematics. Most influential was the discovery of a tight correlation between BH mass and the velocity dispersion of the host-galaxy bulge. It and other correlations led to the belief that BHs and bulges coevolve by regulating each other's growth. New results are now replacing this simple story with a richer and more plausible picture in which BHs correlate differently with different galaxy components. BHs are found in pure-disk galaxies, so classical (elliptical-galaxy-like) bulges are not necessary to grow BHs. But BHs do not correlate with galaxy disks. And any correlations with disk-grown pseudobulges or halo dark matter are so weak as to imply no close coevolution. We suggest that there are four regimes of BH feedback. 1- Local, stochastic feeding of small BHs in mainly bulgeless galaxies involves too little energy to result in coevolution. 2- Global feeding in major, wet galaxy mergers grows giant BHs in short, quasar-like "AGN" events whose feedback does affect galaxies. This makes classical bulges and coreless-rotating ellipticals. 3- At the highest BH masses, maintenance-mode feedback into X-ray gas has the negative effect of helping to keep baryons locked up in hot gas. This happens in giant, core-nonrotating ellipticals. They inherit coevolution magic from smaller progenitors. 4- Independent of any feedback physics, the averaging that results from successive mergers helps to engineer tight BH correlations.

preprint2013arXiv

Coevolution (Or Not) of Supermassive Black Holes and Host Galaxies: Supplemental Material

This is the Supplemental Material to Kormendy and Ho 2013, ARAA, 51, 511 (arXiv:1304.7762). Section S1 summarizes indirect methods that are used to estimate black hole (BH) masses for galaxies with active nuclei (AGNs). Section S2 lists the observational criteria that are used to classify classical and pseudo bulges. The (pseudo)bulge classifications used in the main paper are not based on physical interpretation; rather, they are based on these observational criteria. Section S3 supplements the BH database in Section 5 of the main paper and Section S4 here. It discusses corrections to galaxy and BH parameters, most importantly to 2MASS K-band apparent magnitudes. It presents evidence that corrections are needed because 2MASS misses light at large radii when the images of galaxies subtend large angles on the sky or have shallow outer brightness gradients. Section S4 reproduces essentially verbatim the first part of Section 5 in the main paper, the BH database. It includes the list of BH and host-galaxy properties (Tables 2 and 3). Its most important purpose is to provide all of the notes on individual objects.

preprint2013arXiv

Embedded Star Formation in S4G Galaxy Dust Lanes

Star-forming regions that are visible at 3.6 microns and Halpha but not in the u,g,r,i,z bands of the Sloan Digital Sky survey (SDSS), are measured in five nearby spiral galaxies to find extinctions averaging ~3.8 mag and stellar masses averaging ~5x10^4 Msun. These regions are apparently young star complexes embedded in dark filamentary shock fronts connected with spiral arms. The associated cloud masses are ~10^7 Msun. The conditions required to make such complexes are explored, including gravitational instabilities in spiral shocked gas and compression of incident clouds. We find that instabilities are too slow for a complete collapse of the observed spiral filaments, but they could lead to star formation in the denser parts. Compression of incident clouds can produce a faster collapse but has difficulty explaining the semi-regular spacing of some regions along the arms. If gravitational instabilities are involved, then the condensations have the local Jeans mass. Also in this case, the near-simultaneous appearance of equally spaced complexes suggests that the dust lanes, and perhaps the arms too, are relatively young.

preprint2013arXiv

Evidence for AGN-driven Outflows in Young Radio Quasars

We present near-infrared spectra of young radio quasars [P(1.4GHz) ~ 26-27 W/Hz] selected from the Wide-Field Infrared Survey Explorer. The detected objects have typical redshifts of z ~ 1.6-2.5 and bolometric luminosities ~ 10^47 erg/s. Based on the intensity ratios of narrow emission lines, we find that these objects are mainly powered by active galactic nuclei (AGNs), although star formation contribution cannot be completely ruled out. The host galaxies experience moderate levels of extinction, A(V) ~ 0-1.3 mag. The observed [O III] luminosities and rest-frame J-band magnitudes constrain the black hole masses to lie in the range ~ 10^8.9-10^9.7 solar mass. From the empirical correlation between black hole mass and host galaxy mass, we infer stellar masses of ~ 10^11.3-10^12.2 solar mass. The [O III] line is exceptionally broad, with full width at half maximum ~1300 to 2100 km/s, significantly larger than that of ordinary distant quasars. We argue that these large line widths can be explained by jet-induced outflows, as predicted by theoretical models of AGN feedback.

preprint2013arXiv

Evolution of Broad-line Emission from Active Galactic Nuclei

Apart from viewing-dependent obscuration, intrinsic broad-line emission from active galactic nuclei (AGNs) follows an evolutionary sequence: Type $1 \to 1.2/1.5 \to 1.8/1.9 \to 2$ as the accretion rate onto the central black hole is decreasing. This spectral evolution is controlled, at least in part, by the parameter $L_{\rm bol}/M^{2/3}$, where $L_{\rm bol}$ is the AGN bolometric luminosity and $M$ is the black hole mass. Both this dependence and the double-peaked profiles that emerge along the sequence arise naturally in the disk-wind scenario for the AGN broad-line region.

preprint2013arXiv

Fossil Evidence for the Two-phase Formation of Elliptical Galaxies

Massive early-type galaxies have undergone dramatic structural evolution over the last 10 Gyr. A companion paper shows that nearby elliptical galaxies with M*>1.3x10^{11} M_sun generically contain three photometric subcomponents: a compact inner component with effective radius Re<1 kpc, an intermediate-scale middle component with Re~2.5 kpc, and an extended outer envelope with $R_e \approx 10$ kpc. Here we attempt to relate these substructures with the properties of early-type galaxies observed at higher redshifts. We find that a hypothetical structure formed from combining the inner plus the middle components of local ellipticals follows a strikingly tight stellar mass-size relation, one that resembles the distribution of early-type galaxies at z~1.5. Outside of the central kpc, the median stellar mass surface density profiles of this composite structure agree closest with those of massive galaxies that have similar cumulative number density at 1.5<z<2.0 within the uncertainty. We propose that the central substructures in nearby ellipticals are the evolutionary descendants of the "red nuggets" formed under highly dissipative ("wet") conditions at high redshifts, as envisioned in the initial stages of the two-phase formation scenario recently advocated for massive galaxies. Subsequent accretion, plausibly through dissipationless ("dry") minor mergers, builds the outer regions of the galaxy identified as the outer envelope in our decomposition. The large scatter exhibited by this component on the stellar mass-size plane testifies to the stochastic nature of the accretion events.

preprint2013arXiv

Hα Kinematics of S4G spiral galaxies-I. NGC 864

We present a study of the kinematics of the isolated spiral galaxy NGC 864, using Hα Fabry-Perot data obtained with the GHαFaS instrument at the William Herschel Telescope in La Palma, complemented with images at 3.6 μm, in the R band and in Hα filter, and integral field spectroscopic data. The resulting data cubes and velocity maps allow the study of the kinematics of the galaxy, including in-depth investigations of the rotation curve, velocity moment maps, velocity residual maps, gradient maps and position-velocity diagrams. We find asymmetries in the velocity field in the bar zone, caused by non-circular motions, probably in response to the potential of the bar. We also find a flat-profile bar, in agreement with the strong bar, with the grand design spiral pattern, and with the gap between the ends of the bar and the start of the spiral arms. We quantify the rate of massive star formation, which is concentrated in the two spiral arms.

preprint2013arXiv

On the Origin of Lopsidedness in Galaxies as Determined from the Spitzer Survey of Stellar Structure in Galaxies (S4G)

We study the m=1 distortions (lopsidedness) in the stellar components of 167 nearby galaxies that span a wide range of morphologies and luminosities. We confirm the previous findings of 1) a high incidence of lopsidedness in the stellar distributions, 2) increasing lopsidedness as a function of radius out to at least 3.5 exponential scale lengths, and 3) greater lopsidedness, over these radii, for galaxies of later type and lower surface brightness. Additionally, the magnitude of the lopsidedness 1) correlates with the character of the spiral arms (stronger arm patterns occur in galaxies with less lopsidedness), 2) is not correlated with the presence or absence of a bar, or the strength of the bar when one is present, 3) is inversely correlated to the stellar mass fraction, f_*, within one radial scale length, and 4) correlates directly with f_* measured within the radial range over which we measure lopsidedness. We interpret these findings to mean that lopsidedness is a generic feature of galaxies and does not, generally, depend on a rare event, such as a direct accretion of a satellite galaxy onto the disk of the parent galaxy. While lopsidedness may be caused by several phenomena, moderate lopsidedness (<A_1>_i + <A_1>_o)/2 < 0.3) is likely to reflect halo asymmetries to which the disk responds or a gravitationally self-generated mode . We hypothesize that the magnitude of the stellar response depends both on how centrally concentrated the stars are with respect to the dark matter and whether there are enough stars in the region of the lopsidedness that self-gravity is dynamically important.

preprint2013arXiv

Spectral Energy Distributions of Type 1 AGN in XMM-COSMOS Survey II - Shape Evolution

The mid-infrared to ultraviolet (0.1 -- 10 $μm$) spectral energy distribution (SED) shapes of 407 X-ray-selected radio-quiet type 1 AGN in the wide-field ``Cosmic Evolution Survey" (COSMOS) have been studied for signs of evolution. For a sub-sample of 200 radio-quiet quasars with black hole mass estimates and host galaxy corrections, we studied their mean SEDs as a function of a broad range of redshift, bolometric luminosity, black hole mass and Eddington ratio, and compared them with the Elvis et al. (1994, E94) type 1 AGN mean SED. We found that the mean SEDs in each bin are closely similar to each other, showing no statistical significant evidence of dependence on any of the analyzed parameters. We also measured the SED dispersion as a function of these four parameters, and found no significant dependencies. The dispersion of the XMM-COSMOS SEDs is generally larger than E94 SED dispersion in the ultraviolet, which might be due to the broader ``window function'' for COSMOS quasars, and their X-ray based selection.

preprint2013arXiv

The M87 Black Hole Mass from Gas-dynamical Models of Space Telescope Imaging Spectrograph Observations

The supermassive black hole of M87 is one of the most massive black holes known and has been the subject of several stellar and gas-dynamical mass measurements; however the most recent revision to the stellar-dynamical black hole mass measurement is a factor of about two larger than the previous gas-dynamical determinations. Here, we apply comprehensive gas-dynamical models that include the propagation of emission-line profiles through the telescope and spectrograph optics to new Space Telescope Imaging Spectrograph observations from the Hubble Space Telescope. Unlike the previous gas-dynamical studies of M87, we map out the complete kinematic structure of the emission-line disk within about 40 pc from the nucleus, and find that a small amount of velocity dispersion internal to the gas disk is required to match the observed line widths. We examine a scenario in which the intrinsic velocity dispersion provides dynamical support to the disk, and determine that the inferred black hole mass increases by only 6%. Incorporating this effect into the error budget, we ultimately measure a mass of M_BH = (3.5^{+0.9}_{-0.7}) x 10^9 M_sun (68% confidence). Our gas-dynamical black hole mass continues to differ from the most recent stellar-dynamical mass by a factor of two, underscoring the need for carrying out more cross-checks between the two main black hole mass measurement methods.

preprint2013arXiv

UM 625 Revisited: Multiwavelength Study of A Seyfert 1 Galaxy with a Low-mass Black Hole

UM 625, previously identified as a narrow-line active galactic nucleus (AGN), actually exhibits broad \ha\ and \hb\ lines whose width and luminosity indicate a low black hole mass of $1.6 \times 10^6$ \msun. We present a detailed multiwavelength study of the nuclear and host galaxy properties of UM 625. Analysis of \chandra\ and \xmm\ observations suggests that this system contains a heavily absorbed and intrinsically X-ray weak ($\aox=-1.72$) nucleus. Although not strong enough to qualify as radio-loud, UM 625 does belong to a minority of low-mass AGNs detected in the radio. The broad-band spectral energy distribution constrains the bolometric luminosity to $\lbol\approx(0.5-3)\times10^{43}$ \lum\ and $\lratio\approx0.02-0.15$. A comprehensive analysis of Sloan Digital Sky Survey and {\it Hubble Space Telescope}\ images shows that UM 625 is a nearly face-on S0 galaxy with a prominent, relatively blue pseudobulge (\sersic\ index $n = 1.60$) that accounts for $\sim$60% of the total light in the $R$ band. The extended disk is featureless, but the central $\sim150-400$ pc contains a conspicuous semi-ring of bright, blue star-forming knots, whose integrated ultraviolet luminosity suggests a star formation rate of $\sim$0.3 \msun yr$^{-1}$. The mass of the central black hole roughly agrees with the value predicted from its bulge velocity dispersion but is significantly lower than that expected from its bulge luminosity.

preprint2013arXiv

Unveiling the Structure of Barred Galaxies at 3.6 $μ{\rm m}$ with the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G): I. Disk Breaks

We have performed two-dimensional multicomponent decomposition of 144 local barred spiral galaxies using 3.6 $μ{\rm m}$ images from the Spitzer Survey of Stellar Structure in Galaxies. Our model fit includes up to four components (bulge, disk, bar, and a point source) and, most importantly, takes into account disk breaks. We find that ignoring the disk break and using a single disk scale length in the model fit for Type II (down-bending) disk galaxies can lead to differences of 40% in the disk scale length, 10% in bulge-to-total luminosity ratio (B/T), and 25% in bar-to-total luminosity ratios. We find that for galaxies with B/T $\geq$ 0.1, the break radius to bar radius, $r_{\rm br}/R_{\rm bar}$, varies between 1 and 3, but as a function of B/T the ratio remains roughly constant. This suggests that in bulge-dominated galaxies the disk break is likely related to the outer Lindblad Resonance (OLR) of the bar, and thus moves outwards as the bar grows. For galaxies with small bulges, B/T $<$ 0.1, $r_{\rm br}/R_{\rm bar}$ spans a wide range from 1 to 6. This suggests that the mechanism that produces the break in these galaxies may be different from that in galaxies with more massive bulges. Consistent with previous studies, we conclude that disk breaks in galaxies with small bulges may originate from bar resonances that may be also coupled with the spiral arms, or be related to star formation thresholds.

preprint2013arXiv

X-ray nuclear activity in S4G barred galaxies: No link between bar strength and co-occurrent supermassive black hole fueling

Stellar bars can lead to gas inflow toward the center of a galaxy and stimulate nuclear star formation. However, there is no compelling evidence on whether they also feed a central supermassive black hole: by measuring the fractions of barred active and inactive galaxies, previous studies have yielded conflicting results. In this paper, we aim to understand the lack of observational evidence for bar-driven active galactic nucleus (AGN) activity by studying a sample of 41 nearby (d < 35 Mpc) barred galaxies from the Spitzer Survey for Stellar Structure in Galaxies. We use Chandra observations to measure nuclear 2--10 keV X-ray luminosities and estimate Eddington ratios, together with Spitzer 3.6um imaging to quantify the strength of the stellar bar in two independent ways: (1) from its structure, as traced by its ellipticity and boxiness, and (2) from its gravitational torque Q_b, taken as the maximum ratio of the tangential force to the mean background radial force. In this way, rather than discretizing the presence of both stellar bars and nuclear activity, we are able to account for the continuum of bar strengths and degrees of AGN activity. We find nuclear X-ray sources in 31 out of 41 galaxies with median X-ray luminosity and Eddington ratio of L_X=4.3x10^{38} erg/s and L_bol/L_Edd=6.9x10^{-6} respectively, consistent with low-luminosity AGN activity. Including upper limits for those galaxies without nuclear detections, we find no significant correlation between any of the bar strength indicators and the degree of nuclear activity, irrespective of galaxy luminosity, stellar mass, Hubble type, or bulge size. Strong bars do not favor brighter or more efficient nuclear activity, implying that at least for the low-luminosity regime, supermassive black hole fueling is not closely connected to large scale features.

preprint2012arXiv

A unified picture of breaks and truncations in spiral galaxies from SDSS and S^{4}G imaging

The mechanism causing breaks in the radial surface brightness distribution of spiral galaxies is not yet well known. Despite theoretical efforts, there is not a unique explanation for these features and the observational results are not conclusive. In an attempt to address this problem, we have selected a sample of 34 highly inclined spiral galaxies present both in the Sloan Digital Sky Survey and in the Spitzer Survey of Stellar Structure in Galaxies. We have measured the surface brightness profiles in the five Sloan optical bands and in the 3.6$μm$ Spitzer band. We have also calculated the color and stellar surface mass density profiles using the available photometric information, finding two differentiated features: an innermost break radius at distances of $\sim 8 \pm 1$ kpc [$0.77 \pm 0.06$ $R_{25}$] and a second characteristic radius, or truncation radius, close to the outermost optical extent ($\sim 14 \pm 2$ kpc [$1.09 \pm 0.05$ $R_{25}$]) of the galaxy. We propose in this work that the breaks might be a phenomena related to a threshold in the star formation, while truncations are more likely a real drop in the stellar mass density of the disk associated with the maximum angular momentum of the stars.

preprint2012arXiv

A Uniformly Selected Sample of Low-mass Black Holes in Seyfert 1 Galaxies

We have conducted a systematic search of low-mass black holes (BHs) in active galactic nuclei (AGNs) with broad Halpha emission lines, aiming at building a homogeneous sample that is more complete than previous ones for fainter, less highly accreting sources. For this purpose, we developed a set of elaborate, automated selection procedures and applied it uniformly to the Fourth Data Release of the Sloan Digital Sky Survey. Special attention is given to AGN--galaxy spectral decomposition and emission-line deblending. We define a sample of 309 type 1 AGNs with BH masses in the range $8 \times 10^4$--$2 \times 10^6$ \msun (with a median of $1.2 \times 10^6$ solar mass), using the virial mass estimator based on the broad Halpha line. About half of our sample of low-mass BHs differs from that of Greene & Ho, with 61 of them discovered here for the first time. Our new sample picks up more AGNs with low accretion rates: the Eddington ratios of the present sample range from $<~0.01$ to ~1, with 30% below 0.1. This suggests that a significant fraction of low-mass BHs in the local Universe are accreting at low rates. The host galaxies of the low-mass BHs have luminosities similar to those of $L^*$ field galaxies, optical colors of Sbc spirals, and stellar spectral features consistent with a continuous star formation history with a mean stellar age of less than 1 Gyr.

preprint2012arXiv

Berkeley Supernova Ia Program I: Observations, Data Reduction, and Spectroscopic Sample of 582 Low-Redshift Type Ia Supernovae

In this first paper in a series we present 1298 low-redshift (z\leq0.2) optical spectra of 582 Type Ia supernovae (SNe Ia) observed from 1989 through 2008 as part of the Berkeley SN Ia Program (BSNIP). 584 spectra of 199 SNe Ia have well-calibrated light curves with measured distance moduli, and many of the spectra have been corrected for host-galaxy contamination. Most of the data were obtained using the Kast double spectrograph mounted on the Shane 3 m telescope at Lick Observatory and have a typical wavelength range of 3300-10,400 Ang., roughly twice as wide as spectra from most previously published datasets. We present our observing and reduction procedures, and we describe the resulting SN Database (SNDB), which will be an online, public, searchable database containing all of our fully reduced spectra and companion photometry. In addition, we discuss our spectral classification scheme (using the SuperNova IDentification code, SNID; Blondin & Tonry 2007), utilising our newly constructed set of SNID spectral templates. These templates allow us to accurately classify our entire dataset, and by doing so we are able to reclassify a handful of objects as bona fide SNe Ia and a few other objects as members of some of the peculiar SN Ia subtypes. In fact, our dataset includes spectra of nearly 90 spectroscopically peculiar SNe Ia. We also present spectroscopic host-galaxy redshifts of some SNe Ia where these values were previously unknown. [Abridged]

preprint2012arXiv

Breaks in thin and thick discs of edge-on galaxies imaged in the Spitzer Survey of Stellar Structure in Galaxies (S4G)

Breaks in the radial luminosity profiles of galaxies have been until now mostly studied averaged over discs. Here we study separately breaks in thin and thick discs in 70 edge-on galaxies using imaging from the Spitzer Survey of Stellar Structure in Galaxies. We built luminosity profiles of the thin and the thick discs parallel to midplanes and we found that thin discs often truncate (77%). Thick discs truncate less often (31%), but when they do, their break radius is comparable with that in the thin disc. This suggests either two different truncation mechanisms - one of dynamical origin affecting both discs simultaneously and another one only affecting the thin disc - or a single mechanism that creates a truncation in one disc or in both depending on some galaxy property. Thin discs apparently antitruncate in around 40% of galaxies. However, in many cases, these antitruncations are an artifact caused by the superposition of a thin disc and a thick disc with the latter having a longer scale length. We estimate the real thin disc antitruncation fraction to be less than 15%. We found that the ratio of the thick and thin stellar disc mass is roughly constant (0.2<M_T/M_t<0.7) for circular velocities v_c>120 km/s, but becomes much larger at smaller velocities. We hypothesize that this is due to a combination of a high efficiency of supernova feedback and a slower dynamical evolution in lower-mass galaxies causing stellar thin discs to be younger and less massive than in higher-mass galaxies.

preprint2012arXiv

Chandra Observations of the Nuclear Star Cluster and Ultraluminous X-ray Sources in NGC 2139

We report Chandra observations of the Scd galaxy NGC 2139, which is known to host a recently formed (10^7.6 yrs) nuclear star cluster. The star cluster is undetected in X-rays, with an upper bound on 0.5-7 keV luminosity of L_X < 7.1 x 10^37 erg/s. This bound implies a bolometric accretion luminosity <0.3 percent of the Eddington luminosity for a black hole with the mass (approximately 3400 M_sun) expected from extrapolation of the M-sigma relation. The lack of X-ray emission indicates that a black hole, if present, is not undergoing significant accretion at the current time. While the central cluster is undetected, the data reveal a substantial population of bright X-ray point sources elsewhere in this galaxy, with eight qualifying as ultraluminous X-ray sources with L_X > 10^39 erg/s. We use archival Hubble Space Telescope images to identify candidate optical counterparts for seven Chandra sources, which in most cases have optical luminosities and spatial profiles consistent with star clusters. Compared with other galaxies, the number of luminous X-ray sources in NGC 2139 is larger by a factor of 4 - 10 than expected based on its present star formation rate and stellar mass. This finding can be understood if NGC 2139 has concluded a burst of star formation in the recent past, and suggests that this galaxy could be important for testing the use of X-ray source populations as a chronometer of star formation history.

preprint2012arXiv

Cosmological Evolution of Supermassive Black Holes. II. Evidence for Downsizing of Spin Evolution

The spin is an important but poorly constrained parameter for describing supermassive black holes (SMBHs). Using the continuity equation of SMBH number density, we explicitly obtain the mass-dependent cosmological evolution of the radiative efficiency for accretion, which serves as a proxy for SMBH spin. Our calculations make use of the SMBH mass function of active and inactive galaxies (derived in the first paper of this series), the bolometric luminosity function of active galactic nuclei (AGNs), corrected for the contribution from Compton-thick sources, and the observed Eddington ratio distribution. We find that the radiative efficiency generally increases with increasing black hole mass at high redshifts (z>~1), roughly as η\propto M_bh^0.5, while the trend reverses at lower redshifts, such that the highest efficiencies are attained by the lowest mass black holes. Black holes with M_bh>~10^8.5M_sun maintain radiative efficiencies as high as η~0.3-0.4 at high redshifts, near the maximum for rapidly spinning systems, but their efficiencies drop dramatically (by an order of magnitude) by z~0. The pattern for lower mass holes is somewhat more complicated but qualitatively similar. Assuming that the standard accretion disk model applies, we suggest that the accretion history of SMBHs and their accompanying spins evolve in two distinct regimes: an early phase of prolonged accretion, plausibly driven by major mergers, during which the black hole spins up, then switching to a period of random, episodic accretion, governed by minor mergers and internal secular processes, during which the hole spins down. The transition epoch depends on mass, mirroring other evidence for "cosmic downsizing" in the AGN population; it occurs at z~2 for high-mass black holes, and somewhat later, at z~1, for lower-mass systems.

preprint2012arXiv

Early type galaxies with tidal debris and their scaling relations in the Spitzer Survey of Stellar Structure in Galaxies (S4G)

Tidal debris around galaxies can yield important clues on their evolution. We have identified tidal debris in 11 early type galaxies (T \leq 0) from a sample of 65 early types drawn from the Spitzer Survey of Stellar Structure in Galaxies (S4G). The tidal debris includes features such as shells, ripples and tidal tails. A variety of techniques, including two-dimensional decomposition of galactic structures, was used to quantify the residual tidal features. The tidal debris contributes ~3 - 10% to the total 3.6 μm luminosity of the host galaxy. Structural parameters of the galaxies were estimated using two-dimensional profile fitting. We investigate the locations of galaxies with tidal debris in the Fundamental Plane and Kormendy relation. We find that galaxies with tidal debris lie within the scatter of early type galaxies without tidal features. Assuming that the tidal debris is indicative of recent gravitational interaction or merger, this suggests that these galaxies have either undergone minor merging events so that the overall structural properties of the galaxies are not significantly altered, or they have undergone a major merging events but already have experienced sufficient relaxation and phase-mixing so that their structural properties become similar to those of the non-interacting early type galaxies.

preprint2012arXiv

Physical Properties of the Narrow-Line Region of Low-Mass Active Galaxies

We present spectroscopic observations of 27 active galactic nuclei (AGN) with some of the lowest black hole (BH) masses known. We use the high spectral resolution and small aperture of our Keck data, taken with the Echellette Spectrograph and Imager, to isolate the narrow-line regions (NLRs) of these low-mass BHs. We investigate their emission-line properties and compare them with those of AGN with higher-mass black holes. While we are unable to determine absolute metallicities, some of our objects plausibly represent examples of the low-metallicity AGN described by Groves et al. (2006), based on their [N II]/H_alpha ratios and their consistency with the Kewley & Ellison (2008) mass-metallicity relation. We find tentative evidence for steeper far-UV spectral slopes in lower-mass systems. Overall, NLR emission lines in these low-mass AGN exhibit trends similar to those seen in AGN with higher-mass BHs, such as increasing blueshifts and broadening with increasing ionization potential. Additionally, we see evidence of an intermediate line region whose intensity correlates with L/L_Edd, as seen in higher-mass AGN. We highlight the interesting trend that, at least in these low-mass AGN, the [O III] equivalent width (EW) is highest in symmetric NLR lines with no blue wing. This trend of increasing [O III] EW with line symmetry could be explained by a high covering factor of lower ionization gas in the NLR. In general, low-mass AGN preserve many well-known trends in the structure of the NLR, while exhibiting steeper ionizing continuum slopes and somewhat lower gas-phase metallicities.

preprint2012arXiv

Rapid Infrared Variability of Three Radio-loud Narrow-line Seyfert 1 Galaxies: A View from the Wide-field Infrared Survey Explorer

Using newly released data from the Wide-field Infrared Survey Explorer, we report the discovery of rapid infrared variability in three radio-loud narrow-line Seyfert 1 galaxies (NLS1s) selected from the 23 sources in the sample of Yuan et al. (2008). J0849+5108 and J0948+0022 clearly show intraday variability, while J1505+0326 has a longer measurable time scale within 180 days. Their variability amplitudes, corrected for measurement errors, are $\sim 0.1-0.2$ mag. The detection of intraday variability restricts the size of the infrared-emitting region to $\sim 10^{-3}$ pc, significantly smaller than the scale of the torus but consistent with the base of a jet. The three variable sources are exceptionally radio-loud, have the highest radio brightness temperature among the whole sample, and all show detected $γ$-ray emission in Fermi/LAT observations. Their spectral energy distributions resemble those of low-energy-peaked blazars, with a synchrotron peak around infrared wavelengths. This result strongly confirms the view that at least some radio-loud NLS1s are blazars with a relativistic jet close to our line of sight. The beamed synchrotron emission from the jet contributes significantly to and probably dominates the spectra in the infrared and even optical bands.

preprint2012arXiv

Simultaneous Ultraviolet and Optical Emission-line Profiles of Quasars: Implications for Black Hole Mass Determination

The X-shooter instrument on the VLT was used to obtain spectra of seven moderate-redshift quasars simultaneously covering the spectral range 3000 Ang to 2.5 microns. At z ~ 1.5, most of the prominent broad emission lines in the ultraviolet to optical region are captured in their rest frame. We use this unique dataset, which mitigates complications from source variability, to intercompare the line profiles of C IV 1549, C III] 1909, Mg II 2800, and Halpha and evaluate their implications for black hole mass estimation. We confirm that Mg II and the Balmer lines share similar kinematics and that they deliver mutually consistent black hole mass estimates with minimal internal scatter (< 0.1 dex) using the latest virial mass estimators. Although no virial mass formalism has yet been calibrated for C III], this line does not appear promising for such an application because of the large spread of its velocity width compared to lines of both higher and lower ionization; part of the discrepancy may be due to the difficulty of deblending C III] from its neighboring lines. The situation for C IV is complex and, because of the limited statistics of our small sample, inconclusive. On the one hand, slightly more than half of our sample (4/7) have C IV line widths that correlate reasonably well with Halpha line widths, and their respective black hole mass estimates agree to within ~0.15 dex. The rest, on the other hand, exhibit exceptionally broad C IV profiles that overestimate virial masses by factors of 2-5 compared to Halpha. As C IV is widely used to study black hole demographics at high redshifts, we urgently need to revisit our analysis with a larger sample.

preprint2012arXiv

The Carnegie-Irvine Galaxy Survey. III. The Three-Component Structure of Nearby Elliptical Galaxies

Motivated by recent developments in our understanding of the formation and evolution of massive galaxies, we explore the detailed photometric structure of a representative sample of 94 bright, nearby elliptical galaxies, using high-quality optical images from the Carnegie-Irvine Galaxy Survey. The sample spans a range of environments and stellar masses, from M* = 10^{10.2} to 10^{12.0} solar mass. We exploit the unique capabilities of two-dimensional image decomposition to explore the possibility that local elliptical galaxies may contain photometrically distinct substructure that can shed light on their evolutionary history. Compared with the traditional one-dimensional approach, these two-dimensional models are capable of consistently recovering the surface brightness distribution and the systematic radial variation of geometric information at the same time. Contrary to conventional perception, we find that the global light distribution of the majority (>75%) of elliptical galaxies is not well described by a single Sersic function. Instead, we propose that local elliptical galaxies generically contain three subcomponents: a compact (R_e < 1 kpc) inner component with luminosity fraction f ~ 0.1-0.15; an intermediate-scale (R_e ~ 2.5 kpc) middle component with f ~ 0.2-0.25; and a dominant (f = 0.6), extended (R_e ~ 10 kpc) outer envelope. All subcomponents have average Sersic indices n ~ 1-2, significantly lower than the values typically obtained from single-component fits. The individual subcomponents follow well-defined photometric scaling relations and the stellar mass-size relation. We discuss the physical nature of the substructures and their implications for the formation of massive elliptical galaxies.

preprint2012arXiv

The Low-mass, Highly Accreting Black Hole Associated with the Active Galactic Nucleus 2XMM J123103.2+110648

Optical spectra and images taken with the Baade 6.5 meter Magellan telescope confirm that 2XMM J123103.2+110648, a highly variable X-ray source with an unusually soft spectrum, is indeed associated with a type 2 (narrow-line) active nucleus at a redshift of z = 0.11871. The absence of broad Halpha or Hbeta emission in an otherwise X-ray unabsorbed source suggests that it intrinsically lacks a broad-line region. If, as in other active galaxies, the ionized gas and stars in J1231+1106 are in approximate virial equilibrium, and the black hole mass versus stellar velocity dispersion relation holds, the exceptionally small velocity dispersion of 33.5 km/s for [O III] 5007 implies that the black hole mass is approximately 10^5 solar masses, among the lowest ever detected. Such a low black hole mass is consistent with the general characteristics of the host, a small, low-luminosity, low-mass disk galaxy. We estimate the Eddington ratio of the black hole to be > 0.5, in good agreement with expectations based on the X-ray properties of the source.

preprint2012arXiv

The S4G perspective on circumstellar dust extinction of AGB stars in M100

We examine the effect of circumstellar dust extinction on the near-IR contribution of asymptotic giant branch (AGB) stars in intermediate-age clusters throughout the disk of M100. For our sample of 17 AGB-dominated clusters we extract optical-to-mid-IR SEDs and find that NIR brightness is coupled to the mid-IR dust emission in such a way that a significant reduction of AGB light, of up to 1 mag in K-band, follows from extinction by the dust shell formed during this stage. Since the dust optical depth varies with AGB chemistry (C-rich or O-rich), our results suggest that the contribution of AGB stars to the flux from their host clusters will be closely linked to the metallicity and the progenitor mass of the AGB star, to which dust chemistry and mass-loss rate are sensitive. Our sample of clusters--each the analogue of a ~1 Gyr old post-starburst galaxy--has implications within the context of mass and age estimation via SED modelling at high z: we find that the average ~0.5 mag extinction estimated here may be sufficient to reduce the AGB contribution in (rest-frame) K-band from ~70%, as predicted in the latest generation of synthesis models, to ~35%. Our technique for selecting AGB-dominated clusters in nearby galaxies promises to be effective for discriminating the uncertainties associated with AGB stars in intermediate-age populations that plague age and mass estimation in high-z galaxies.

preprint2012arXiv

Two Populations of Molecular Clouds in the Antennae Galaxies

Super star clusters --- extremely massive clusters found predominately in starburst environments --- are essential building blocks in the formation of galaxies and thought to dominate star formation in the high-redshift universe. However, the transformation from molecular gas into these ultra-compact star clusters is not well understood. To study this process, we used the Submillimeter Array and the Plateau de Bure Interferometer to obtain high angular resolution (~1.5" or 160 pc) images of the Antennae overlap region in CO(2--1) to search for the molecular progenitors of the super star clusters. We resolve the molecular gas distribution into a large number of clouds, extending the differential cloud mass function down to a 5σcompleteness limit of 3.8x10^5 M_sun. We identify a distinct break in the mass function around log M_mol/M_sun ~ 6.5, which separates the molecular clouds into two distinct populations. The smaller, less massive clouds reside in more quiescent areas in the region, while the larger, more massive clouds cluster around regions of intense star formation. A broken power-law fit to the mass function yields slopes of α= -1.39+/-0.10 and α= -1.44+/-0.14 for the low- and high-mass cloud population, well-matched to the mass function found for super star clusters in the Antennae galaxies. We find large velocity gradients and velocity dispersions at the locations of intense star formation, suggestive of compressive shocks. It is likely that these environmental factors contribute to the formation of the observed massive molecular clouds and super star clusters in the Antennae galaxies.

preprint2012arXiv

Two-Component Structure of the Hbeta Broad-Line Region in Quasars. I. Evidence from Spectral Principal Component Analysis

We report on a spectral principal component analysis (SPCA) of a sample of 816 quasars, selected to have small Fe II velocity shifts with spectral coverage in the rest wavelength range 3500--5500 Å. The sample is explicitly designed to mitigate spurious effects on SPCA induced by Fe II velocity shifts. We improve the algorithm of SPCA in the literature and introduce a new quantity, \emph{the fractional-contribution spectrum}, that effectively identifies the emission features encoded in each eigenspectrum. The first eigenspectrum clearly records the power-law continuum and very broad Balmer emission lines. Narrow emission lines dominate the second eigenspectrum. The third eigenspectrum represents the Fe II emission and a component of the Balmer lines with kinematically similar intermediate velocity widths. Correlations between the weights of the eigenspectra and parametric measurements of line strength and continuum slope confirm the above interpretation for the eigenspectra. Monte Carlo simulations demonstrate the validity of our method to recognize cross talk in SPCA and firmly rule out a single-component model for broad Hbeta. We also present the results of SPCA for four other samples that contain quasars in bins of larger Fe II velocity shift; similar eigenspectra are obtained. We propose that the Hbeta-emitting region has two kinematically distinct components: one with very large velocities whose strength correlates with the continuum shape, and another with more modest, intermediate velocities that is closely coupled to the gas that gives rise to Fe II emission.

preprint2012arXiv

X-ray Properties of Intermediate-mass Black Holes in Active Galaxies. III. Spectral Energy Distribution and Possible Evidence for Intrinsically X-ray-weak AGNs

We present a systematic X-ray study, the third in a series, of 49 active galactic nuclei with intermediate-mass black holes (IMBH; ~10^5-10^6 M_sun) using Chandra observations. We detect 42 out of 49 targets with a 0.5-2 keV X-ray luminosity 10^41-10^43 erg/s. We perform spectral fitting for the 10 objects with enough counts (>200), and they are all well fit by a simple power-law model modified by Galactic absorption, with no sign of significant intrinsic absorption. While we cannot fit the X-ray spectral slope directly for the rest of the sample, we estimate it from the hardness ratio and find a range of photon indices consistent with those seen in more luminous and massive objects. The X-ray-to-optical spectral slope (alphaox) of our IMBH sample is systematically flatter than in active galaxies with more massive black holes, consistent with the well-known correlation between alphaox and UV luminosity. Thanks to the wide dynamic range of our sample, we find evidence that alphaox increases with decreasing M_BH as expected from accretion disk models, where the UV emission systematically decreases as M_BH decreases and the disk temperature increases. We also find a long tail toward low alphaox values. While some of these sources may be obscured, given the high L_bol/L_Eddington values in the sample, we argue that some may be intrinsically X-ray-weak, perhaps owing to a rare state that radiates very little coronal emission.

preprint2011arXiv

AEGIS: Demographics of X-ray and Optically Selected AGNs

We develop a new diagnostic method to classify galaxies into AGN hosts, star-forming galaxies, and absorption-dominated galaxies by combining the [O III]/Hbeta ratio with rest-frame U-B color. This can be used to robustly select AGNs in galaxy samples at intermediate redshifts (z<1). We compare the result of this optical AGN selection with X-ray selection using a sample of 3150 galaxies with 0.3<z<0.8 and I_AB<22, selected from the DEEP2 Galaxy Redshift Survey and the All-wavelength Extended Groth Strip International Survey (AEGIS). Among the 146 X-ray sources in this sample, 58% are classified optically as emission-line AGNs, the rest as star-forming galaxies or absorption-dominated galaxies. The latter are also known as "X-ray bright, optically normal galaxies" (XBONGs). Analysis of the relationship between optical emission lines and X-ray properties shows that the completeness of optical AGN selection suffers from dependence on the star formation rate and the quality of observed spectra. It also shows that XBONGs do not appear to be a physically distinct population from other X-ray detected, emission-line AGNs. On the other hand, X-ray AGN selection also has strong bias. About 2/3 of all emission-line AGNs at L_bol>10^44 erg/s in our sample are not detected in our 200 ks Chandra images, most likely due to moderate or heavy absorption by gas near the AGN. The 2--7 keV detection rate of Seyfert 2s at z~0.6 suggests that their column density distribution and Compton-thick fraction are similar to that of local Seyferts. Multiple sample selection techniques are needed to obtain as complete a sample as possible.

preprint2011arXiv

Black Hole Mass and Bulge Luminosity for Low-mass Black Holes

We study the scaling between bulge magnitude and central black hole (BH) mass in galaxies with virial BH masses < 10^6 solar mass. Based on careful image decomposition of a snapshot Hubble Space Telescope I-band survey, we found that these BHs are found predominantly in galaxies with pseudobulges. Here we show that the \mbulge\ relation for the pseudobulges at low mass is significantly different from classical bulges with BH masses >10^7 solar mass. Specfically, bulges span a much wider range of bulge luminosity, and on average the luminosity is larger, at fixed black hole mass. The trend holds both for the active galaxies from Bentz et al. and the inactive sample of Gultekin et al. and cannot be explained by differences in stellar populations, as it persists when we use dynamical bulge masses. Put another way, the ratio between bulge and BH mass is much larger than $\sim 1000$ for our sample. This is consistent with recent suggestions that black hole mass does not scale with the pseudobulge luminosity. The low-mass scaling relations appear to flatten, consistent with predictions from Volonteri & Natarajan for massive seed BHs.

preprint2011arXiv

Cosmological evolution of supermassive black holes. I. mass function at 0<z<2

We present the mass function of supermassive black holes (SMBHs) over the redshift range z=0-2, using the latest deep luminosity and mass functions of field galaxies to constrain the masses of their spheroids, which we relate to SMBH mass through the empirical correlation between SMBH and spheroid mass (the M_bh-M_sph relation). In ddition to luminosity fading of the stellar content of the spheroids, we carefully consider the variation of the bulge-to-total luminosity ratio of the galaxy populations and the M_bh/M_sph ratio, which, according to numerous recent studies, evolves rapidly with redshift. The SMBH mass functions derived from the galaxy luminosity and mass functions show very good agreement, both in shape and in normalization. The resultant SMBH mass function and integrated mass density for the local epoch (z~0) match well those derived independently by other studies. Consistent with other evidence for cosmic downsizing, the upper end of the mass function remains roughly constant since $z\approx2$, while the space density of lower mass black holes undergoes strong evolution. We carefully assess the impact of various sources of uncertainties on our calculations. A companion paper uses the mass function derived in this work to determine the radiative efficiency of black hole accretion and constraints that can be imposed on the cosmological evolution of black hole spin.

preprint2011arXiv

Exploring the Low-Mass End of The M-Sigma Relation with Active Galaxies

We present new measurements of stellar velocity dispersions, using spectra obtained with the Keck Echellette Spectrograph and Imager (ESI) and the Magellan Echellette (MagE), for 76 Seyfert 1 galaxies from the recent catalogue of Greene & Ho. These objects were selected from the Sloan Digital Sky Survey (SDSS) to have estimated black hole (BH) masses below 2\times10^6 M\odot. Combining our results with previous ESI observations of similar objects, we obtain an expanded sample of 93 galaxies and examine the relation between BH mass and velocity dispersion (the M-Sigma relation) for active galaxies with low BH masses. The low-mass active galaxies tend to follow the extrapolation of the M-Sigma relation of inactive galaxies. Including results for active galaxies of higher BH mass from the literature, we find a zero pointα= 7.68\pm0.08 and slope of β= 3.32\pm0.22 for the M-Sigma relation [in the form log Mbh=α+βlog(σ*/200 km s-1)], with intrinsic scatter of 0.46\pm0.03 dex. This result is consistent, within the uncertainties, with the slope of the M-Sigma relation for reverberation-mapped active galaxies with BH masses from 10^6 to 10^9 M\odot. For the subset of our sample having morphological information from Hubble Space Telescope images, we examine the slope of the M-Sigma relation separately for subsamples of barred and unbarred host galaxies, and find no significant evidence for a difference in slope. We do find a mild offset between low-inclination and high-inclination disk galaxies, such that more highly inclined galaxies tend to have larger σ* at a given value of BH mass, presumably due to the contribution of disk rotation within the spectroscopic aperture. We also find that the velocity dispersion of the ionized gas trace the stellar velocity dispersion well for this large sample of low-mass Seyfert 1 galaxies.

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.

preprint2011arXiv

Medium Resolution Spectral Library of Late-Type Stellar Templates in Near-Infrared Band

We present medium resolution (R = 5000-6000) spectra in the near-infrared band, 1.4-1.8 microns, for template stars in G, K, and M types observed by the echelle spectrometer, IRCS, at the SUBARU 8.2 m telescope. The identification of lines is based on the spectra of Arcturus (K2 III) in the literature. We measured the equivalent of widths and compared our results to those of Meyer et al. (1998). We conclude that our spectral resolution (R = 6000) data can investigate more accurately the properties of lines in stellar spectra. The library of the template stellar spectra in ASCII format are available for download on the World Wide Web.

preprint2011arXiv

Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: I. Correcting for contamination by PAH, hot dust, and intermediate age stars

With the aim of constructing accurate 2D maps of the stellar mass distribution in nearby galaxies from S4G 3.6 and 4.5 micron images, we report on the separation of the light from old stars from the emission contributed by contaminants (e.g. hot dust and the 3.3 micron PAH feature). Results for a small sample of six disk galaxies (NGC 1566, NGC 2976, NGC 3031, NGC 3184, NGC 4321, and NGC 5194) with a range of morphological properties, dust contents and star formation histories are presented to demonstrate our approach. We use an Independent Component Analysis (ICA) technique designed to separate statistically independent source distributions, maximizing the distinction in the [3.6]-[4.5] colors of the sources. The technique also removes emission from intermediate-age evolved red objects with a low mass-to-light ratio, such as asymptotic giant branch (AGB) and red supergiant (RSG) stars, revealing maps of the underlying old distribution of light with [3.6]-[4.5] colors consistent with the colors of K and M giants. Contaminants are identified via comparison to the non-stellar emission imaged at 8 microns, which is dominated by the broad PAH feature. Using the measured 3.6/8 micron ratio to select the individual contaminants, we find that hot dust and PAH together contribute between ~5-15% to the integrated light at 3.6 microns, while light from regions dominated by intermediate-age stars accounts for only 1-5%. Locally, however, the contribution from either contaminant can reach much higher levels; dust contributes on average 22% to the emission in star-forming regions throughout the sample, while intermediate age-stars contribute upwards of 50% in localized knots. The removal of these contaminants with ICA leaves maps of the old stellar disk that retain a high degree of structural information and are ideally suited for tracing the stellar mass, as will be the focus in a companion paper.

preprint2011arXiv

The Carnegie-Irvine Galaxy Survey. I. Overview and Atlas of Optical Images

The Carnegie-Irvine Galaxy Survey (CGS) is a long-term program to investigate the photometric and spectroscopic properties of a statistically complete sample of 605 bright (B_T < 12.9 mag), southern (Dec. < 0) galaxies using the facilities at Las Campanas Observatory. This paper, the first in a series, outlines the scientific motivation of CGS, defines the sample, and describes the technical aspects of the optical broadband (BVRI) imaging component of the survey, including details of the observing program, data reduction procedures, and calibration strategy. The overall quality of the images is quite high, in terms of resolution (median seeing 1"), field of view (8.9' X 8.9'), and depth (median limiting surface brightness 27.5, 26.9, 26.4, and 25.3 mag/arcsec2 in the B, V, R, and I bands, respectively). We prepare a digital image atlas showing several different renditions of the data, including three-color composites, star-cleaned images, stacked images to enhance faint features, structure maps to highlight small-scale features, and color index maps suitable for studying the spatial variation of stellar content and dust. In anticipation of upcoming science analyses, we tabulate an extensive set of global properties for the galaxy sample. These include optical isophotal and photometric parameters derived from CGS itself, as well as published information on multiwavelength (ultraviolet, U-band, near-infrared, far-infrared) photometry, internal kinematics (central stellar velocity dispersions, disk rotational velocities), environment (distance to nearest neighbor, tidal parameter, group or cluster membership), and H I content. The digital images and science-level data products will be made publicly accessible to the community.

preprint2011arXiv

The Carnegie-Irvine Galaxy Survey. II. Isophotal Analysis

The Carnegie-Irvine Galaxy Survey (CGS) is a comprehensive investigation of the physical properties of a complete, representative sample of 605 bright (B_T <= 12.9 mag) galaxies in the southern hemisphere. This contribution describes the isophotal analysis of the broadband (BVRI) optical imaging component of the project. We pay close attention to sky subtraction, which is particularly challenging for some of the large galaxies in our sample. Extensive crosschecks with internal and external data confirm that our calibration and sky subtraction techniques are robust with respect to the quoted measurement uncertainties. We present a uniform catalog of one-dimensional radial profiles of surface brightness and geometric parameters, as well as integrated colors and color gradients. Composite profiles highlight the tremendous diversity of brightness distributions found in disk galaxies and their dependence on Hubble type. A significant fraction of S0 and spiral galaxies exhibit non-exponential profiles in their outer regions. We perform Fourier decomposition of the isophotes to quantify non-axisymmetric deviations in the light distribution. We use the geometric parameters, in conjunction with the amplitude and phase of the m=2 Fourier mode, to identify bars and quantify their size and strength. Spiral arm strengths are characterized using the m=2 Fourier profiles and structure maps. Finally, we utilize the information encoded in the m=1 Fourier profiles to measure disk lopsidedness. The databases assembled here and in Paper I lay the foundation for forthcoming scientific applications of CGS.

preprint2011arXiv

The Host Galaxies of Low-mass Black Holes

Using HST observations of 147 host galaxies of low-mass black holes (BHs), we systematically study the structures and scaling relations of these active galaxies. Our sample is selected to have central BHs with virial masses ~10^5-10^6 solar mass. The host galaxies have total I-band magnitudes of -23.2<M_I<-18.8 mag and bulge magnitudes of -22.9<M_I<-16.1 mag. Detailed bulge-disk-bar decompositions with GALFIT show that 93% of the galaxies have extended disks, 39% have bars and 5% have no bulges at all at the limits of our observations. Based on the Sersic index and bulge-to-total ratio, we conclude that the majority of the galaxies with disks are likely to contain pseudobulges and very few of these low-mass BHs live in classical bulges. The fundamental plane of our sample is offset from classical bulges and ellipticals in a way that is consistent with the scaling relations of pseudobulges. The sample has smaller velocity dispersion at fixed luminosity in the Faber-Jackson plane, compared with classical bulges and elliptical galaxies. The galaxies without disks are structurally more similar to spheroidals than to classical bulges according to their positions in the fundamental plane, especially the Faber-Jackson projection. Overall, we suggest that BHs with mass < 10^6 solar mass live in galaxies that have evolved secularly over the majority of their history. A classical bulge is not a prerequisite to host a black hole.

preprint2011arXiv

The unusual vertical mass distribution of NGC 4013 seen through the Spitzer Survey of Stellar Structure in Galaxies (S4G)

NGC 4013 is a nearby Sb edge-on galaxy known for its "prodigious" HI warp and its "giant" tidal stream. Previous work on this unusual object shows that it cannot be fitted satisfactorily by a canonical thin+thick disk structure. We have produced a new decomposition of NGC 4013, considering three stellar flattened components (thin+thick disk plus an extra and more extended component) and one gaseous disk. All four components are considered to be gravitationally coupled and isothermal. To do so, we have used the 3.6micron images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We find evidence for NGC4013 indeed having a thin and a thick disk and an extra flattened component. This smooth and extended component (scaleheight z_EC kpc) could be interpreted as a thick disk or as a squashed ellipsoidal halo and contains ~20% of the total mass of all three stellar components. We argue it is unlikely to be related to the ongoing merger or due to the off-plane stars from a warp in the other two disk components. Instead, we favor a scenario in which the thick disk and the extended component were formed in a two-stage process, in which an initially thick disk has been dynamically heated by a merger soon enough in the galaxy history to have a new thick disk formed within it.

preprint2011arXiv

Thick disks of edge-on galaxies seen through the Spitzer Survey of Stellar Structure in Galaxies (S4G): Lair of missing baryons?

Most, if not all, disk galaxies have a thin (classical) disk and a thick disk. In most models thick disks are thought to be a necessary consequence of the disk formation and/or evolution of the galaxy. We present the results of a study of the thick disk properties in a sample of carefully selected edge-on galaxies with types ranging from T=3 to T=8. We fitted one-dimensional luminosity profiles with physically motivated functions - the solutions of two stellar and one gaseous isothermal coupled disks in equilibrium - which are likely to yield more accurate results than other functions used in previous studies. The images used for the fits come from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We found that thick disks are on average more massive than previously reported, mostly due to the selected fitting function. Typically, the thin and the thick disk have similar masses. We also found that thick disks do not flare significantly within the observed range in galactocentric radii and that the ratio of thick to thin disk scaleheights is higher for galaxies of earlier types. Our results tend to favor an in situ origin for most of the stars in the thick disk. In addition the thick disk may contain a significant amount of stars coming from satellites accreted after the initial build-up of the galaxy and an extra fraction of stars coming from the secular heating of the thin disk by its own overdensities. Assigning thick disk light to the thin disk component may lead to an underestimate of the overall stellar mass in galaxies, because of different mass to light ratios in the two disk components. On the basis of our new results, we estimate that disk stellar masses are between 10% and 50% higher than previously thought and we suggest that thick disks are a reservoir of "local missing baryons".

preprint2011arXiv

What controls the Fe II strength in active galactic nuclei?

We used a large, homogeneous sample of 4178 z <= 0.8 Seyfert 1 galaxies and QSOs selected from the Sloan Digital Sky Survey to investigate the strength of Fe II emission and its correlation with other emission lines and physical parameters of active galactic nuclei. We find that the strongest correlations of almost all the emission-line intensity ratios and equivalent widths (EWs) are with the Eddington ratio (L/L_{Edd}), rather than with the continuum luminosity at 5100Å (L_{5100}) or black hole mass (M_{BH}); the only exception is the EW of ultraviolet Fe II emission, which does not correlate at all with broad-line width, L_{5100}, M_{BH}, or L/L_{Edd}. By contrast, the intensity ratios of both the ultraviolet and optical Fe II emission to Mg II λ2800 correlate quite strongly with L/L_{Edd}. Interestingly, among all the emission lines in the near-UV and optical, the EW of narrow optical Fe II emission has the strongest correlation with L/L_{Edd}. We suggest that the variation of the emission-line strength in active galaxies is regulated by L/L_{Edd} because it governs the global distribution of the column density of the clouds gravitationally bound in the line-emitting region, as well as its overall gas supply. The systematic dependence on L/L_{Edd} must be corrected when using the FeII/MgII intensity ratio as a measure of the Fe/Mg abundance ratio to study the history of chemical evolution in QSO environments.

preprint2010arXiv

Detailed Decomposition of Galaxy Images. II. Beyond Axisymmetric Models

We present a two-dimensional (2-D) fitting algorithm (GALFIT, Version 3) with new capabilities to study the structural components of galaxies and other astronomical objects in digital images. Our technique improves on previous 2-D fitting algorithms by allowing for irregular, curved, logarithmic and power-law spirals, ring and truncated shapes in otherwise traditional parametric functions like the Sersic, Moffat, King, Ferrer, etc., profiles. One can mix and match these new shape features freely, with or without constraints, apply them to an arbitrary number of model components and of numerous profile types, so as to produce realistic-looking galaxy model images. Yet, despite the potential for extreme complexity, the meaning of the key parameters like the Sersic index, effective radius or luminosity remain intuitive and essentially unchanged. The new features have an interesting potential for use to quantify the degree of asymmetry of galaxies, to quantify low surface brightness tidal features beneath and beyond luminous galaxies, to allow more realistic decompositions of galaxy subcomponents in the presence of strong rings and spiral arms, and to enable ways to gauge the uncertainties when decomposing galaxy subcomponents. We illustrate these new features by way of several case studies that display various levels of complexity.

preprint2010arXiv

The Prevalence of Narrow Optical Fe II Emission Lines in Type 1 Active Galactic Nuclei

From detailed spectral analysis of a large sample of low-redshift active galactic nuclei (AGNs) selected from the Sloan Digital Sky Survey, we demonstrate---statistically for the first time---that narrow optical Fe II emission lines, both permitted and forbidden, are prevalent in type 1 AGNs. Remarkably, these optical lines are completely absent in type 2 AGNs, across a wide luminosity range, from Seyfert 2 galaxies to type 2 quasars. We suggest that the narrow FeII-emitting gas is confined to a disk-like geometry in the innermost regions of the narrow-line region on physical scales smaller than the obscuring torus.

preprint2010arXiv

The Radio Properties of Type 2 Quasars

This paper presents the first high-resolution and high-sensitivity study of the radio properties of optically selected type~2 quasars. We used the Very Large Array at 8.4 GHz to observe 59 sources drawn from the Sloan Digital Sky Survey sample of Zakamska et al. (2003). The detection rate of our survey is 59% (35/59), comparable to the detection rate in FIRST at 1.4 GHz. Ongoing star formation, although present, contributes negligible radio emission at the current sensitivity limit. Comparing the radio powers with the [O III] $λ$5007 luminosities, we find that roughly 15%$\pm$5% of the sample can be considered radio-loud. Intriguingly, the vast majority of the detected sources in our sample fall in a region intermediate between those traditionally occupied by radio-loud and radio-quiet quasars. Moreover, most of these "radio-intermediate" sources tend to have flat or inverted radio spectra, which we speculate may be caused by free-free absorption by ionized gas in the narrow-line region. The incidence of flat-spectrum sources in type 2 quasars appears to be much higher than in type 1 quasars, in apparent violation of the simple orientation-based unified model for active galaxies.

preprint2010arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S^4G)

The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.

preprint2010arXiv

The thick disk in the galaxy NGC 4244 from S4G imaging

If thick disks are ubiquitous and a natural product of disk galaxy formation and/or evolution processes, all undisturbed galaxies which have evolved during a significant fraction of a Hubble time should have a thick disk. The late-type spiral galaxy NGC 4244 has been reported as the only nearby edge-on galaxy without a confirmed thick disk. Using data from the Spitzer Survey of Stellar Structure in Galaxies (S4G) we have identified signs of two disk components in this galaxy. The asymmetries between the light profiles on both sides of the mid-plane of NGC 4244 can be explained by a combination of the galaxy not being perfectly edge-on and a certain degree of opacity of the thin disk. We argue that the subtlety of the thick disk is a consequence of either a limited secular evolution in NGC 4244, a small fraction of stellar material in the fragments which built the galaxy, or a high amount of gaseous accretion after the formation of the galaxy.

preprint2009arXiv

A Quintet of Black Hole Mass Determinations

We report five new measurements of central black hole masses based on STIS and WFPC2 observations with the Hubble Space Telescope and on axisymmetric, three-integral, Schwarzschild orbit-library kinematic models. We selected a sample of galaxies within a narrow range in velocity dispersion that cover a range of galaxy parameters (including Hubble type and core/power-law surface density profile) where we expected to be able to resolve the galaxy's sphere of influence based on the predicted value of the black hole mass from the M-sigma relation. We find masses in units of 10^8 solar masses for the following galaxies: NGC 3585, M_BH = 3.4 (+1.5, -0.6); NGC 3607, M_BH = 1.2 (+0.4, -0.4); NGC 4026, M_BH = 2.1 (+0.7, -0.4); and NGC 5576, M_BH = 1.8 (+0.3, -0.4), all significantly excluding M_BH = 0. For NGC 3945, M_BH = 0.09 (+0.17, -0.21), which is significantly below predictions from M-sigma and M-L relations and consistent with M_BH = 0, though the presence of a double bar in this galaxy may present problems for our axisymmetric code.

preprint2009arXiv

Estimating Black Hole Masses in Active Galactic Nuclei Using the MgII 2800 Emission Line

We investigate the relationship between the linewidths of broad Mg II \lambda2800 and Hbeta in active galactic nuclei (AGNs) to refine them as tools to estimate black hole (BH) masses. We perform a detailed spectral analysis of a large sample of AGNs at intermediate redshifts selected from the Sloan Digital Sky Survey, along with a smaller sample of archival ultraviolet spectra for nearby sources monitored with reverberation mapping. Careful attention is devoted to accurate spectral decomposition, especially in the treatment of narrow-line blending and Fe II contamination. We show that, contrary to popular belief, the velocity width of Mg II tends to be smaller than that of Hbeta, suggesting that the two species are not cospatial in the broad-line region. Using these findings and recently updated BH mass measurements from reverberation mapping, we present a new calibration of the empirical prescriptions for estimating virial BH masses for AGNs using the broad Mg II and Hbeta lines. We show that the BH masses derived from our new formalisms show subtle but important differences compared to some of the mass estimators currently used in the literature.

preprint2009arXiv

Magellan Spectroscopy of Low-Redshift Active Galactic Nuclei

We present an atlas of moderate-resolution (R ~ 1200-1600) optical spectra of 94 low-redshift (z < 0.5) active galactic nuclei taken with the Magellan 6.5 m Clay Telescope. The spectra mostly cover the rest-frame region 3600-6000 Ang. All the objects have preexisting Hubble Space Telescope imaging, and they were chosen as part of an ongoing program to investigate the relationship between black hole mass and their host galaxy properties. A significant fraction of the sample has no previous quantitative spectroscopic measurements in the literature. We perform spectral decomposition of the spectra and present detailed fits and basic measurements of several commonly used broad and narrow emission lines, including [O II] 3727, He II 4686, Hbeta, and [O III] 4959, 5007. Eight of the objects are narrow-line sources that were previously misclassified as broad-line (type 1) Seyfert galaxies; of these, five appear not to be accretion-powered.

preprint2009arXiv

On the Disappearance of the Broad-Line Region in Low-Luminosity Agns

The disk-wind scenario for the broad-line region (BLR) and toroidal obscuration in active galactic nuclei predicts the disappearance of the BLR at low luminosities. In accordance with the model predictions, data from a nearly complete sample of nearby AGNs show that the BLR disappears at luminosities lower than $5\times\E{39} (M/10^7\Mo)^{2/3}$ erg s$^{-1)$, where $M$ is the black hole mass. The radiative efficiency of accretion onto the black hole is $\la \E{-3}$ for these sources, indicating that their accretion is advection-dominated.

preprint2009arXiv

The M-sigma and M-L Relations in Galactic Bulges and Determinations of their Intrinsic Scatter

We derive improved versions of the relations between supermassive black hole mass (M_BH) and host-galaxy bulge velocity dispersion (sigma) and luminosity (L) (the M-sigma and M-L relations), based on 49 M_BH measurements and 19 upper limits. Particular attention is paid to recovery of the intrinsic scatter (epsilon_0) in both relations. We find log(M_BH / M_sun) = alpha + beta * log(sigma / 200 km/s) with (alpha, beta, epsilon_0) = (8.12 +/- 0.08, 4.24 +/- 0.41, 0.44 +/- 0.06) for all galaxies and (alpha, beta, epsilon_0) = (8.23 +/- 0.08, 3.96 +/- 0.42, 0.31 +/- 0.06) for ellipticals. The results for ellipticals are consistent with previous studies, but the intrinsic scatter recovered for spirals is significantly larger. The scatter inferred reinforces the need for its consideration when calculating local black hole mass function based on the M-sigma relation, and further implies that there may be substantial selection bias in studies of the evolution of the M-sigma relation. We estimate the M-L relationship as log(M_BH / M_sun) = alpha + beta * log(L_V / 10^11 L_sun,V) of (alpha, beta, epsilon_0) = (8.95 +/- 0.11, 1.11 +/- 0.18, 0.38 +/- 0.09); using only early-type galaxies. These results appear to be insensitive to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We show that culling the sample according to the resolution of the black hole's sphere of influence biases the relations to larger mean masses, larger slopes, and incorrect intrinsic residuals.

preprint2009arXiv

X-ray Properties of Intermediate-Mass Black Holes in Active Galaxies. II. X-ray-Bright Accretion and Possible Evidence for Slim Disks

We present X-ray properties of optically-selected intermediate-mass (~10^5--10^6 M_Sun) black holes (BHs) in active galaxies (AGNs), using data from the Chandra X-Ray Observatory. Our observations are a continuation of a pilot study by Greene & Ho (2007). Of the 8 objects observed, 5 are detected with X-ray luminosities in the range L_0.5-2 keV = 10^41--10^43 erg s^-1, consistent with the previously observed sample. Objects with enough counts to extract a spectrum are well fit by an absorbed power law. We continue to find a range of soft photon indices 1 < Γ_s < 2.7, where N(E) \propto E^-Γ_s, consistent with previous AGN studies, but generally flatter than other narrow-line Seyfert 1 active nuclei (NLS1s). The soft photon index correlates strongly with X-ray luminosity and Eddington ratio, but does not depend on BH mass. There is no justification for the inclusion of any additional components, such as a soft excess, although this may be a function of the relative inefficiency of detecting counts above 2 keV in these relatively shallow observations. As a whole, the X-ray-to-optical spectral slope α_ox is flatter than in more massive systems, even other NLS1s. Only X-ray-selected NLS1s with very high Eddington ratios share a similar α_ox. This is suggestive of a physical change in the accretion structure at low masses and at very high accretion rates, possibly due to the onset of slim disks. Although the detailed physical explanation for the X-ray loudness of these intermediate-mass BHs is not certain, it is very striking that targets selected on the basis of optical properties should be so distinctly offset in their broader spectral energy distributions.

preprint2008arXiv

A Stellar Dynamical Measurement of the Black Hole Mass in the Maser Galaxy NGC 4258

We determine the mass of the black hole at the center of the spiral galaxy NGC 4258 by constructing axisymmetric dynamical models of the galaxy. These models are constrained by high spatial resolution imaging and long-slit spectroscopy of the nuclear region obtained with the {\em Hubble Space Telescope}, complemented by ground-based observations extending to larger radii. Our best mass estimate is $\MBH = (3.3 \pm 0.2) \times 10^7 \MSun $ for a distance of 7.28 Mpc (statistical errors only). This is within 15% of $ (3.82\pm 0.01) \times 10^7 \MSun$, the mass determined from the kinematics of water masers (rescaled to the same distance) assuming they are in Keplerian rotation in a warped disk. The construction of accurate dynamical models of NGC 4258 is somewhat compromised by an unresolved active nucleus and color gradients, the latter caused by variations in the stellar population and/or obscuring dust. These problems are not present in the $\sim 30$ other black hole mass determinations from stellar dynamics that have been published by us and other groups; thus, the relatively close agreement between the stellar dynamical mass and the maser mass in NGC 4258 enhances our confidence in the black hole masses determined in other galaxies from stellar dynamics using similar methods and data of comparable quality.

preprint2008arXiv

Nuclear Activity in Nearby Galaxies

A significant fraction of nearby galaxies show evidence of weak nuclear activity unrelated to normal stellar processes. Recent high-resolution, multiwavelength observations indicate that the bulk of this activity derives from black hole accretion with a wide range of accretion rates. The low accretion rates that typify most low-luminosity active galactic nuclei induce significant modifications to their central engine. The broad-line region and obscuring torus disappear in some of the faintest sources, and the optically thick accretion disk transforms into a three-component structure consisting of an inner radiatively inefficient accretion flow, a truncated outer thin disk, and a jet or outflow. The local census of nuclear activity supports the notion that most, perhaps all, bulges host a central supermassive black hole, although the existence of active nuclei in at least some late-type galaxies suggests that a classical bulge is not a prerequisite to seed a nuclear black hole.

preprint2007arXiv

Fe K Line Profile in Low-redshift Quasars: Average Shape and Eddington Ratio Dependence

We analyze X-ray spectra of 43 Palomar-Green quasars observed with {\it XMM-Newton} in order to investigate their mean Fe K line profile and its dependence on physical properties. The continuum spectra of 39 objects are well reproduced by a model consisting of a power law and a blackbody modified by Galactic absorption. The spectra of the remaining four objects require an additional power-law component absorbed with a column density of $\sim 10^{23} {\rm cm}^{-2}$. We fit the entire sample simultaneously to derive average Fe line parameters by assuming a common Fe line shape. The Fe line is relatively narrow ($σ=0.36$ keV), with a center energy of 6.48 keV and a mean equivalent width (EW) of 248 eV. By combining black hole masses estimated from the virial method and bolometric luminosities derived from full spectral energy distributions, we examine the dependence of the Fe K line profile on Eddington ratio. As the Eddington ratio increases, the line becomes systematically stronger (EW = 130 to 280 eV), broader ($σ=0.1$ to 0.7 keV), and peaks at higher energies (6.4 to 6.8 keV). This result suggests that the accretion rate onto the black hole directly influences the geometrical structure and ionization state of the accretion disk.

preprint2007arXiv

The Mass Function of Active Black Holes in the Local Universe

We present the first measurement of the black hole (BH) mass function for broad-line active galaxies in the local Universe. Using the ~9000 broad-line active galaxies from the Fourth Data Release of the Sloan Digital Sky Survey, we construct a broad-line luminosity function that agrees very well with the local soft X-ray luminosity function. Using standard virial relations, we then convert observed broad-line luminosities and widths into BH masses. A mass function constructed in this way has the unique capability to probe the mass region <10^6 M_sun, which, while insignificant in terms of total BH mass density, nevertheless may place important constraints on the mass distribution of seed BHs in the early Universe. The characteristic local active BH has a mass of ~10^7 M_sun radiating at 10% of the Eddington rate. The active fraction is a strong function of BH mass; at both higher and lower masses the active mass function falls more steeply than one would infer from the distribution of bulge luminosity. The deficit of local massive radiating BHs is a well-known phenomenon, while we present the first robust measurement of a decline in the space density of active BHs at low mass.

preprint2007arXiv

The Masses of Nuclear Black Holes in Luminous Elliptical Galaxies and Implications for the Space Density of the Most Massive Black Holes

Black hole masses predicted from the Mbh-sigma relationship conflict with those predicted from the Mbh-L relationship for the most luminous galaxies, such as brightest cluster galaxies (BCGs). This is because stellar velocity dispersion, sigma, increases only weakly with L for BCGs and other giant ellipticals. The Mbh-L relationship predicts that the most luminous BCGs may have Mbh approaching 10^{10}M_sol, while the M-sigma relationship always predicts Mbh<3X10^9M_sol. We argue that the Mbh-L relationship is a plausible or even preferred description for BCGs and other galaxies of similar luminosity. If cores in central stellar density are formed by binary BHs, then the inner-core cusp radius, r_gamma, may be an independent witness of Mbh. Using structural parameters derived from a large sample of early-type galaxies observed by HST, we argue that L is superior to sigma as an indicator of r_gamma in luminous galaxies. The observed r_gamma-Mbh relationship for 11 core galaxies with measured Mbh appears to be consistent with the Mbh-L relationship for BCGs. BCGs have large cores appropriate for their large luminosities that may be difficult to generate with the modest BH masses inferred from the Mbh-sigma relationship. Mbh~L may be expected to hold for BCGs, if they were formed in dissipationless mergers, which should preserve ratio of BH to stellar mass. This picture appears to be consistent with the slow increase in sigma with L and the more rapid increase in effective radii with L seen in BCGs. If BCGs have large BHs commensurate with their luminosities, then the local BH mass function for Mbh>3X10^9M_sol may be nearly an order of magnitude richer than that inferred from the Mbh-sigma relationship. The volume density of QSOs at earlier epochs may favor the predictions from the Mbh-L relationship.

preprint2007arXiv

The radio properties of type II quasars

Quasars (of type I) are the luminous analogs of type I Seyfert galaxies. Within the framework of unified models of active galaxies, the population of quasars of type II recently discovered with the Sloan Digital Sky Survey are the luminous analogs of type II Seyfert galaxies. Since our knowledge and understanding of the radio properties of these type II quasars are very limited, we have performed Very Large Array observations for a sub-sample of such sources. Our detection rate of 61% is consistent with the detection rate for other AGN samples. We do not find a correlation between radio and [O III] $λ$5007 luminosities for these sources. Although the distribution of spectral indices is similar to that of the 3C sources, the lack of dependence of radio luminosity on [O III] $λ$5007 luminosity suggests that not all sources in the sample are genuine AGNs.

preprint2006arXiv

Constraints on the Star Formation Rate in Active Galaxies

The [O II] 3727 emission line is often used as an indicator of star formation rate in extragalactic surveys, and it can be an equally effective tracer of star formation in systems containing luminous active galactic nuclei (AGNs). In order to investigate the ongoing star formation rate of the host galaxies of AGNs, we measured the strength of [O II] and other optical emission lines from a large sample (~ 3600) of broad-line (Type 1) AGNs selected from the Sloan Digital Sky Survey. We performed a set of photoionization calculations to help evaluate the relative contribution of stellar and nonstellar photoionization to the observed strength of [O II]. Consistent with the recent study of Ho (2005), we find that the observed [O II] emission can be explained entirely by photoionization from the AGN itself, with little or no additional contribution from HII regions. This indicates that the host galaxies of Type 1 AGNs experience very modest star formation concurrent with the optically active phase of the nucleus. By contrast, we show that the sample of ``Type 2'' quasars selected from the Sloan Digital Sky Survey does exhibit substantially stronger HII emission consistent with an elevated level of star formation, a result that presents a challenge to the simplest form of the AGN unification model.

preprint2006arXiv

HST/STIS Spectra of Nuclear Star Clusters in Spiral Galaxies: Dependence of Age and Mass On Hubble Type

(Abridged) We study the nuclear star clusters in spiral galaxies of various Hubble types using spectra obtained with STIS on-board HST. We observed the nuclear clusters in 40 galaxies, selected from two previous HST/WFPC2 imaging surveys. The spectra provide a better separation of cluster light from underlying galaxy light than is possible with ground-based spectra. To infer the star formation history, metallicity and dust extinction, we fit weighted superpositions of single-age stellar population templates to the spectra. The luminosity-weighted age ranges from 10 Myrs to 10 Gyrs. The stellar populations of NCs are generally best fit as a mixture of populations of different ages. This indicates that NCs did not form in a single event, but instead they had additional star formation long after the oldest stars formed. On average, the sample clusters in late-type spirals have a younger luminosity-weighted mean age than those in early-type spirals (log(age/yr) = 8.37+/-0.25 vs. 9.23+/-0.21). The average cluster masses are smaller in late-type spirals than in early-type spirals (log(M/Msun) = 6.25+/-0.21 vs. 7.63+/-0.24), and exceed the masses typical of globular clusters. The cluster mass correlates strongly with both the Hubble type of the host galaxy and the luminosity of its bulge. The latter correlation has the same slope as the well-known correlation between supermassive black hole mass and bulge luminosity. The properties of both nuclear clusters and black holes are therefore intimately connected to the properties of the host galaxy.

preprint2006arXiv

Stellar populations in the nuclei of late-type spiral galaxies

(Abridged) As part of an ongoing effort to study the stellar nuclei of very late-type, bulge-less spirals, we present results from a high-resolution spectroscopic survey of nine such nuclear star clusters, undertaken with VLT/UVES. We fit the spectra with population synthesis models and measure Lick-type indices to determine mean luminosity-weighted ages, which range from 4.1*10^7 to 1.1*10^10 years and are insensitive to assumed metallicity or internal extinction. The average metallicity of nuclear clusters in late-type spirals is slightly sub-solar (<Z> = 0.015) but shows significant scatter. The fits also show that the nuclear cluster spectra are best described by a mix of several generations of stars. This is supported by the fact that only models with composite stellar populations yield mass-to-light ratios that match those obtained from dynamical measurements. The last star formation episode was on average 34 Myr ago, while all clusters experienced some star formation in the last 100 Myr. We thus conclude that the nuclear clusters undergo repeated episodes of star formation. The robustness with respect to possible contamination from the underlying galaxy disk is demonstrated by comparison to spectra obtained with HST/STIS. Combining these results with those from Walcher et al. (2005), we have thus shown that the stellar nuclei of these bulge-less galaxies are massive and dense star clusters that form stars recurrently until the present day. This unique set of properties is likely due to the central location of these clusters in their host galaxies.

preprint2006arXiv

The Centers of Early-Type Galaxies with HST. VI. Bimodal Central Surface Brightness Profiles

We combine the results from several HST investigations of the central structure of early-type galaxies to generate a large sample of parameterized surface photometry. The studies included were those that used the "Nuker law" to characterize the inner light distributions of the galaxies. The sample comprises WFPC1 and WFPC2 V band observations published earlier by our group, R band WFPC2 photometry of Rest et al., NICMOS H band photometry by Ravindranath et al. and Quillen et al., and the BCG WFPC2 I band photometry of Laine et al. The distribution of the logarithmic slopes of the central profiles strongly affirms that the central structure of elliptical galaxies with Mv < -19 is bimodal, based on both parametric and non-parametric analysis. At the HST resolution limit, most galaxies are either power-law systems, which have steep cusps in surface brightness, or core systems, which have shallow cusps interior to a steeper envelope brightness distribution. A rapid transition between the two forms occurs over the luminosity range -22 < Mv < -20, with cores dominating at the highest luminosities, and power-laws at the lowest. There are a few "intermediate" systems that have both cusp slopes and total luminosities that fall within the core/power-law transition, but they are rare and do not fill in the overall bimodal distribution of cusp slopes. These results are inconsistent with the Ferrarese et al. Virgo Cluster Survey (VCS) analysis. However, using galaxies common to the VCS samples, we demonstrate that the VCS models of the cusps are either a poor match to the observations or consist of forms fitted to the galaxy envelopes and extrapolated inward to the HST resolution limit.

preprint2006arXiv

The Host Galaxy of the Quasar HE 0450-2958

A recent study suggests that the quasar HE 0450-2958 is hosted by a galaxy substantially fainter than that inferred from the correlation between black hole mass and bulge luminosity. As this result has significant bearings on galaxy and black hole evolution, we revisit the issue by performing an independent analysis of the data, using a two-dimensional image fitting technique. We indeed find no evidence of a host galaxy either, but, due to the brightness of the quasar and uncertainties in the point-spread function, the limits are fairly weak. To derive an upper limit on the host galaxy luminosity, we perform simulations to deblend the quasar from the host under conditions similar to those actually observed. We find that the host galaxy has an absolute magnitude upper limit of -20 < M_V < -21, in good agreement with the previous determination. Since this limit is consistent with the value predicted from the current best estimate of the black hole mass, there is no compelling evidence that the quasar HE 0450-2958 has an abnormally underluminous host galaxy. We also show that, contrary to previous claims, the companion galaxy to HE 0450-2958 should not be be regarded as an ultraluminous infrared galaxy.

preprint2003arXiv

A Low-Mass Central Black Hole in the Bulgeless Seyfert 1 Galaxy NGC 4395

NGC 4395 is one of the least luminous and nearest known type 1 Seyfert galaxies, and it also lacks a bulge. We present an HST I-band image of its nuclear region and Keck high-resolution (8 km/s) echelle spectra containing the Ca II near-infrared triplet. In addition to the unresolved point source, there is a nuclear star cluster of size r ~ 3.9 pc; the upper limit on its velocity dispersion is only 30 km/s. We thus derive an upper limit of 6.2x10^6 solar masses for the mass of the compact nucleus. Based on the amount of spatially resolved light in the HST image, a sizable fraction of this is likely to reside in stars. Hence, this estimate sets a stringent upper limit on the mass of the central black hole. We argue, from other lines of evidence, that the true mass of the black hole is likely to be 10^4-10^5 solar masses. Although the black hole is much less massive than those thought to exist in classical active galactic nuclei, its accretion rate of L_bol/L_Edd ~ 2x10^-2 to 2x10^-3 is consistent with the mass-luminosity relation obeyed by classical AGNs. This may explain why NGC 4395 has a high-excitation (Seyfert) emission-line spectrum; active galaxies having low-ionization spectra seem to accrete at significantly lower rates. NGC 4395, a pure disk galaxy, demonstrates that supermassive black holes are not associated exclusively with bulges.

preprint2002arXiv

Detailed Structural Decomposition of Galaxy Images

We present a two-dimensional (2-D) fitting algorithm (GALFIT) designed to extract structural components from galaxy images, with emphasis on closely modeling light profiles of spatially well-resolved, nearby galaxies observed with the Hubble Space Telescope. Our algorithm improves on previous techniques in two areas, by being able to simultaneously fit a galaxy with an arbitrary number of components, and with optimization in computation speed, suited for working on large galaxy images. We use 2-D models such as the ``Nuker'' law, the Sersic (de Vaucouleurs) profile, an exponential disk, and Gaussian or Moffat functions. The azimuthal shapes are generalized ellipses that can fit disky and boxy components. Many galaxies with complex isophotes, ellipticity changes, and position-angle twists can be modeled accurately in 2-D. When examined in detail, we find that even simple-looking galaxies generally require at least three components to be modeled accurately, rather than the one or two components more often employed. We illustrate this by way of 7 case studies, which include regular and barred spiral galaxies, highly disky lenticular galaxies, and elliptical galaxies displaying various levels of complexities. A useful extension of this algorithm is to accurately extract nuclear point sources in galaxies. We compare 2-D and 1-D extraction techniques on simulated images of galaxies having nuclear slopes with different degrees of cuspiness, and we then illustrate the application of the program to several examples of nearby galaxies with weak nuclei.

preprint2001arXiv

Central Structural Parameters of Early-Type Galaxies as Viewed with HST/NICMOS

We present surface photometry for the central regions of a sample of 33 early-type (E, S0, and S0/a) galaxies observed at 1.6 microns (H band) using the Hubble Space Telescope (HST). We employ a new technique of two-dimensional fitting to extract quantitative parameters for the bulge light distribution and nuclear point sources, taking into consideration the effects of the point-spread function. Parameterizing the bulge profile with a ``Nuker'' law, we confirm that the central surface-brightness distributions largely fall into two categories, each of which correlates with the global properties of the galaxies. ``Core'' galaxies tend to be luminous ellipticals with boxy or pure elliptical isophotes, whereas ``power-law'' galaxies are preferentially lower luminosity systems with disky isophotes. Unlike most previous studies, however, we do not find a clear gap in the distribution of inner cusp slopes; several objects have inner cusp slopes (0.3 < gamma < 0.5) which straddle the regimes conventionally defined for core and power-law type galaxies. The nature of these intermediate objects is unclear. We draw attention to two objects in the sample which appear to be promising cases of galaxies with isothermal cores that are not the brightest members of a cluster. Unresolved nuclear point sources are found in about 50% of the sample galaxies, roughly independent of profile type, with magnitudes in the range m^{nuc}_H = 12.8 to 17.4 mag, which correspond to M_H^{nuc} = -12.8 to -18.4 mag. (Abridged)

preprint2000arXiv

A Relationship Between Nuclear Black Hole Mass and Galaxy Velocity Dispersion

We describe a correlation between the mass M_BH of a galaxy's central black hole and the luminosity-weighted line-of-sight velocity dispersion sigma_e within the half-light radius. The result is based on a sample of 26 galaxies, including 13 galaxies with new determinations of black hole masses from Hubble Space Telescope measurements of stellar kinematics. The best-fit correlation is M_BH = 1.2 (+-0.2) x 10^8 M_sun (sigma_e/200 km/s)^(3.75 (+-0.3))over almost three orders of magnitude in M_BH; the scatter in M_BH at fixed sigma_e is only 0.30 dex and most of this is due to observational errors. The M_BH-sigma_e relation is of interest not only for its strong predictive power but also because it implies that central black hole mass is constrained by and closely related to properties of the host galaxy's bulge.

preprint1999arXiv

Discovery of Radio Outbursts in the Active Nucleus of M81

The low-luminosity active galactic nucleus of M81 has been monitored at centimeter wavelengths since early 1993 as a by-product of radio programs to study the radio emission from Supernova 1993J. The extensive data sets reveal that the nucleus experienced several radio outbursts during the monitoring period. At 2 and 3.6 cm, the main outburst occurred roughly in the beginning of 1993 September and lasted for approximately three months; at longer wavelengths, the maximum flux density decreases, and the onset of the burst is delayed. These characteristics qualitatively resemble the standard model for adiabatically expanding radio sources, although certain discrepancies between the observations and the theoretical predictions suggest that the model is too simplistic. In addition to the large-amplitude, prolonged variations, we also detected milder changes in the flux density at 3.6 cm and possibly at 6 cm on short (less than 1 day) timescales. We discuss a possible association between the radio activity and an optical flare observed during the period that the nucleus was monitored at radio wavelengths.

preprint1999arXiv

Radio Cores in Low-Luminosity AGN: ADAFs or Jets?

We have surveyed two large samples of nearby low-luminosity AGN with the VLA to search for flat-spectrum radio cores, similar to Sgr A* in the Galactic Center. Roughly one third of all galaxies are detected (roughly one half if HII transition objects are excluded from the sample), many of which have compact radio cores. Follow-up observations with the VLBA have confirmed that these cores are non-thermal in origin, with lower limits for the brightness temperatures around ~10^8 K. The brightest of these are resolved into linear structures. The radio spectral indices of the cores are quite flat (alpha~0), with no evidence for the highly inverted radio cores predicted in the ADAF model. Spectrum and morphology of the compact radio emission is typical for radio jets seen also in more luminous AGN. The emission-line luminosity seems to be correlated with the radio core flux. Together with the VLBI observations this suggests that optical and radio emission in at least half the low-luminosity Seyferts and LINERs are black hole powered. We find only a weak correlation between bulge luminosity and radio flux and an apparently different efficiency between elliptical and spiral galaxies for producing radio emission at a given optical luminosity.

preprint1997arXiv

A Search for "Dwarf" Seyfert Nuclei. III. Spectroscopic Parameters and Properties of the Host Galaxies

We have completed an optical spectroscopic survey of the nuclear regions (r < 200 pc) of a large sample of nearby galaxies. Although the main objectives of the survey are to search for low-luminosity active galactic nuclei and to quantify their luminosity function, the database can be used for a variety of other purposes. This paper presents measurements of the spectroscopic parameters for the 418 emission-line nuclei, along with a compilation of the global properties of all 486 galaxies in the survey. Stellar absorption generally poses a serious obstacle to obtaining accurate measurement of emission lines in nearby galactic nuclei. We describe a procedure for removing the starlight from the observed spectra in an efficient and objective manner. The main parameters of the emission lines (intensity ratios, fluxes, profile widths, equivalent widths) are measured and tabulated, as are several stellar absorption-line and continuum indices useful for studying the stellar population. Using standard nebular diagnostics, we determine the probable ionization mechanisms of the emission-line objects. The resulting spectral classifications provide extensive information on the demographics of emission-line nuclei in the local universe. This new catalog contains over 200 objects showing spectroscopic evidence for recent star formation and an equally large number of active galactic nuclei, including 46 which show broad H-alpha emission. These samples will serve as the basis of future studies of nuclear activity in nearby galaxies.

preprint1997arXiv

A Search for "Dwarf" Seyfert Nuclei. IV. Nuclei with Broad H-alpha Emission

We present the results of an optical spectroscopic survey designed to search for low-luminosity, "dwarf" Seyfert nuclei in a magnitude-limited sample of 486 bright, northern galaxies. Moderate-resolution spectra of exceptionally high quality were obtained in part to detect broad H-alpha emission, similar in character to, but much weaker than, the broad permitted lines that define type 1 Seyfert nuclei. One of the primary goals of the survey is to better quantify the faint end of the luminosity function of active galactic nuclei. This paper describes the subset of nuclei showing definite or probable evidence of broad H-alpha emission. We outline the procedures for determining the presence of this elusive spectral feature, steps for its quantitative measurement, and the associated systematic errors. Of the 211 emission-line galaxies classified as having Seyfert or LINER nuclei in our survey, the broad H-alpha line was detected with confidence in 34 objects, and with less certainty in another 12. Most of the detections are reported for the first time, and the detection rate represents a lower limit to the true incidence of active nuclei harboring a broad emission-line region. These statistics imply that broad-lined active nuclei are much more common than previously believed: they exist in at least 20% of all galaxies spectroscopically classified as "active" and in more than 10% of all luminous galaxies at the current epoch.