Source author record

Matthew A. Malkan

Matthew A. Malkan 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

68works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

68 published item(s)

preprint2026arXiv

A census of star-formation and gas mass tracers in two lensed $z \sim 4$ dusty star-forming galaxies

We present new and archival Atacama Large Millimeter/submillimeter Array (ALMA) observations of two strongly lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope survey, SPT0418-47 $(z = 4.225)$ and SPT2147-50 $(z = 3.760)$. We study the [C II], CO(7-6), [C I](2-1), and, in SPT0418-47, $p$-H$_2$O emission, which along with the underlying continuum (rest-frame 160 $μ$m and 380 $μ$m) are routinely used as tracers of gas mass and/or star-formation rate (SFR). We perform a pixel-by-pixel analysis of both sources in the image plane to study the resolved Kennicutt-Schmidt relation, finding generally good agreement between the slopes of the SFR versus gas mass surface density using the different tracers. Using lens modeling methods, we find that the dust emission is more compact than the line emission in both sources, with CO(7-6) and [C I](2-1) similar in extent and [C II] the most extended, reminiscent of recent findings of extended [C II] spatial distributions in galaxies at similar cosmic epochs. We develop the [C I](2-1) / CO(7-6) flux density ratio as an observable proxy for gas depletion timescale ($τ_{\rm dep}$), which can be applied to large samples of DSFGs, in lieu of more detailed inferences of this timescale which require analysis of observations at multiple wavelengths. Furthermore, the extended [C II] emission in both sources, compared to the total continuum and line emission, suggests that [C II], used in recent years as a molecular gas mass and SFR tracer in high-$z$ galaxies, may not always be a suitable tracer of these physical quantities.

preprint2026arXiv

A precessing jet from an active galactic nucleus drives gas outflow from a disk galaxy

To reproduce observed galaxy properties, cosmological simulations require that massive galaxies experience feedback from active galactic nuclei, which regulates star formation within those galaxies. However, the energetics and timescales of these feedback processes are poorly constrained. We combine optical, infrared, sub-millimeter and radio observations of the active galaxy VV 340a, hosting a low-power jet launched from a supermassive black hole at its center. We find that the jet undergoes precession, with a period of (8.2 $\pm~$5.5) $\times~$10$^5$ years, and drives an outflow of gas at a rate of 19.4 $\pm~$7.9 solar masses per year. The jet shocks the gas, producing highly ionized plasma extending several kiloparsecs from the nucleus. The outflow ejects sufficient gas from the galaxy to influence its star formation rate.

preprint2023arXiv

New Identifications and Multi-wavelength Properties of Extragalactic Fermi Gamma-Ray Sources in the SPT-SZ Survey Field

The fourth Fermi Large Area Telescope (LAT) catalog (4FGL) contains 5064 $γ$-ray sources detected at high significance, but 26% of them still lack associations at other wavelengths. The SPT-SZ survey, conducted between 2008 and 2011 with the South Pole Telescope (SPT), covers 2500 $\mathrm{deg^2}$ of the Southern sky in three millimeter-wavelength (mm) bands and was used to construct a catalog of nearly 5000 emissive sources. In this study, we introduce a new cross-matching scheme to search for multi-wavelength counterparts of extragalactic $γ$-ray sources using a mm catalog. We apply a Poissonian probability to evaluate the rate of spurious false associations and compare the multi-wavelength associations from the radio, mm, near-infrared, and X-ray with 4FGL $γ$-ray sources. In the SPT-SZ survey field, 85% of 4FGL sources are associated with mm counterparts. These mm sources include 94% of previously associated 4FGL sources and 56% of previously unassociated 4FGL sources. The latter group contains 40 4FGL sources for which SPT has provided the first identified counterparts. Nearly all of the SPT-associated 4FGL sources can be described as flat-spectrum radio quasars or blazars. We find that the mm band is the most efficient wavelength for detecting $γ$-ray blazars when considering both completeness and purity. We also demonstrate that the mm band correlates better to the $γ$-ray band than the radio or X-ray bands. With the next generation of CMB experiments, this technique can be extended to greater sensitivities and more sky area to further complete the identifications of the remaining unknown $γ$-ray blazars.

preprint2022arXiv

Beyond the Local Volume. II. Population Scaleheights and Ages of Ultracool Dwarfs in Deep HST/WFC3 Parallel Fields

Ultracool dwarfs represent a significant proportion of stars in the Milky Way,and deep samples of these sources have the potential to constrain the formation history and evolution of low-mass objects in the Galaxy. Until recently, spectral samples have been limited to the local volume (d<100 pc). Here, we analyze a sample of 164 spectroscopically-characterized ultracool dwarfs identified by Aganze et al. (2022) in the Hubble Space Telescope WFC3 Infrared Spectroscopic Parallel (WISP) Survey and 3D-HST. We model the observed luminosity function using population simulations to place constraints on scaleheights, vertical velocity dispersions and population ages as a function of spectral type. Our star counts are consistent with a power-law mass function and constant star formation history for ultracool dwarfs, with vertical scaleheights 249$_{-61}^{+48}$ pc for late M dwarfs, 153$_{-30}^{+56}$ pc for L dwarfs, and 175$_{-56}^{+149}$ pc for T dwarfs. Using spatial and velocity dispersion relations, these scaleheights correspond to disk population ages of 3.6$_{-1.0}^{+0.8}$ for late M dwarfs, 2.1$_{-0.5}^{+0.9}$ Gyr for L dwarfs, and 2.4$_{-0.8}^{+2.4}$ Gyr for T dwarfs, which are consistent with prior simulations that predict that L-type dwarfs are on average a younger and less dispersed population. There is an additional 1-2 Gyr systematic uncertainty on these ages due to variances in age-velocity relations. We use our population simulations to predict the UCD yield in the JWST PASSAGES survey, a similar and deeper survey to WISPS and 3D-HST, and find that it will produce a comparably-sized UCD sample, albeit dominated by thick disk and halo sources.

preprint2022arXiv

Chaotic and Clumpy Galaxy Formation in an Extremely Massive Reionization-Era Halo

The SPT0311-58 system at z=6.900 is an extremely massive structure within the reionization epoch, and offers a chance to understand the formation of galaxies in an extreme peak in the primordial density field. We present 70mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and CII 158um emission in the central pair of galaxies and reach physical resolution ~100-350pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800Myr later at z~4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6<z<8 galaxies identified in rest-UV/optical deep field surveys, but with star formation rates elevated by ~3-5x. Internally, the clumps themselves bear close resemblance to greatly scaled-up versions of virialized cloud-scale structures identified in low-redshift galaxies. Our observations are qualitatively similar to the chaotic and clumpy assembly within massive halos seen in simulations of high-redshift galaxies.

preprint2022arXiv

First Census of Gas-phase Metallicity Gradients of Star-forming Galaxies in Overdense Environments at Cosmic Noon

We report the first spatially resolved measurements of gas-phase metallicity radial gradients in star-forming galaxies in overdense environments at $z\gtrsim2$. The spectroscopic data are acquired by the \mg\ survey, a Hubble Space Telescope (HST) cycle-28 medium program. This program is obtaining 45 orbits of WFC3/IR grism spectroscopy in the density peak regions of three massive galaxy protoclusters (BOSS 1244, BOSS 1542 and BOSS 1441) at $z=2-3$. Our sample in the BOSS 1244 field consists of 20 galaxies with stellar-mass ranging from $10^{9.0}$ to $10^{10.3}$ \Msun\ , star formation rate (SFR) from 10 to 240 \Msun\,yr$^{-1}$, and global gas-phase metallicity (\oh) from 8.2 to 8.6. At $1σ$ confidence level, 2/20 galaxies in our sample show positive (inverted) gradients -- the relative abundance of oxygen increasing with galactocentric radius, opposite the usual trend. Furthermore, 1/20 shows negative gradients and 17/20 are consistent with flat gradients. This high fraction of flat/inverted gradients is uncommon in simulations and previous observations conducted in blank fields at similar redshifts. To understand this, we investigate the correlations among various observed properties of our sample galaxies. We find an anticorrelation between metallicity gradient and global metallicity of our galaxies residing in extreme overdensities, and a marked deficiency of metallicity in our massive galaxies as compared to their coeval field counterparts. We conclude that the cold-mode gas accretion plays an active role in shaping the chemical evolution of galaxies in the protocluster environments, diluting their central chemical abundance, and flattening/inverting their metallicity gradients.

preprint2022arXiv

GOALS-JWST: Tracing AGN Feedback on the Star-Forming ISM in NGC 7469

We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study AGN feedback in the local Universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ~100 pc scales. The starburst ring exhibits prominent Polycyclic Aromatic Hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ~30%, and a total star formation rate of $\rm 10 - 30 \ M_\odot$/yr derived from fine structure and recombination emission lines. Using pure rotational lines of H2, we detect 1.2$\times$10$^{7} \rm \ M_\odot$ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker towards the central source, where larger and possibly more ionized grains dominate the emission. However, the bulk of the dust and molecular gas in the ring appears unaffected by the ionizing radiation or the outflowing wind from the AGN. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the co-evolution of supermassive black holes and their hosts.

preprint2022arXiv

Multi-Phase ISM in the z = 5.7 Hyperluminous Starburst SPT0346-52

SPT0346-52 (z=5.7) is the most intensely star-forming galaxy discovered by the South Pole Telescope, with Sigma_SFR ~ 4200 Msol yr^-1 kpc^-2. In this paper, we expand on previous spatially-resolved studies, using ALMA observations of dust continuum, [NII]205 micron, [CII]158 micron, [OI]146 micron, and undetected [NII]122 micron and [OI]63 micron emission to study the multi-phase interstellar medium (ISM) in SPT0346-52. We use pixelated, visibility-based lens modeling to reconstruct the source-plane emission. We also model the source-plane emission using the photoionization code CLOUDY and find a supersolar metallicity system. We calculate T_dust = 48.3 K and lambda_peak = 80 micron, and see line deficits in all five lines. The ionized gas is less dense than comparable galaxies, with n_e < 32 cm^-3, while ~20% of the [CII]158 emission originates from the ionized phase of the ISM. We also calculate the masses of several phases of the ISM. We find that molecular gas dominates the mass of the ISM in SPT0346-52, with the molecular gas mass ~4x higher than the neutral atomic gas mass and ~100x higher than the ionized gas mass.

preprint2022arXiv

The Lick AGN Monitoring Project 2016: Dynamical Modeling of Velocity-Resolved H\b{eta} Lags in Luminous Seyfert Galaxies

We have modeled the velocity-resolved reverberation response of the H\b{eta} broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitioring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the H\b{eta} BLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as log10(FWHM/σ), on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad H\b{eta} emission line and the Eddington ratio, when using the root-mean-square spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends.

preprint2021arXiv

SOFIA observations of far-IR fine-structure lines in galaxies to measure metallicity

We present new and archival SOFIA FIFI-LS far-IR spectroscopic observations of 25 local galaxies of either the [OIII]52um and/or the [NIII]57um lines. Including other 31 galaxies from Herschel-PACS, we discuss a local sample of 47 galaxies, including HII region, luminous IR, low-metallicity dwarf and Seyfert galaxies. Analyzing the mid- to far-IR fine-structure lines of this sample, we assess the metallicity and compare with the optical spectroscopy estimates. Using the IR, we find a similar O/H--N/O relation to that known in the optical. As opposite, we find systematically lower N/O IR abundances when compared to the optical determinations, especially at high values of N/O (log(N/O) > -0.8). We explore various hypotheses to account for this difference: (i) difference in ionization structure traced by optical (O+, N+ regions) versus IR lines (O++, N++ regions); (ii) contamination of diffuse ionized gas affecting the optical lines used to compute the N/O abundance; (iii) dust obscuration affecting the optical-based determinations. However, we have not found any correlation of the Delta(N/O)= (N/O)_OPT-(N/O)_IR with either ionization, or electron density, or optical extinction. We speculatively suggest that accretion of metal-poor gas from the circumgalactic medium could provide an explanation for this difference, because the rapid decrease of total abundances during infall is followed by a N/O ratio decrease due to primary production of young - possibly embedded - massive stars, are preferentially traced by the IR diagnostics, while optical diagnostics would better trace the secondary production, when both N/O and O/H abundance ratios will increase.

preprint2021arXiv

The Lick AGN Monitoring Project 2016: Velocity-Resolved Hβ Lags in Luminous Seyfert Galaxies

We carried out spectroscopic monitoring of 21 low-redshift Seyfert 1 galaxies using the Kast double spectrograph on the 3-m Shane telescope at Lick Observatory from April 2016 to May 2017. Targeting active galactic nuclei (AGN) with luminosities of λLλ (5100 Å) = 10^44 erg/s and predicted Hβ lags of 20-30 days or black hole masses of 10^7-10^8.5 Msun, our campaign probes luminosity-dependent trends in broad-line region (BLR) structure and dynamics as well as to improve calibrations for single-epoch estimates of quasar black hole masses. Here we present the first results from the campaign, including Hβ emission-line light curves, integrated Hβ lag times (8-30 days) measured against V-band continuum light curves, velocity-resolved reverberation lags, line widths of the broad Hβ components, and virial black hole mass estimates (10^7.1-10^8.1 Msun). Our results add significantly to the number of existing velocity-resolved lag measurements and reveal a diversity of BLR gas kinematics at moderately high AGN luminosities. AGN continuum luminosity appears not to be correlated with the type of kinematics that its BLR gas may exhibit. Follow-up direct modeling of this dataset will elucidate the detailed kinematics and provide robust dynamical black hole masses for several objects in this sample.

preprint2021arXiv

The mass-metallicity relation at cosmic noon in overdense environments: first results from the MAMMOTH-Grism HST slitless spectroscopic survey

The MAMMOTH-Grism slitless spectroscopic survey is a Hubble Space Telescope (HST) cycle-28 medium program, which is obtaining 45 orbits of WFC3/IR grism spectroscopy in the density peak regions of three massive galaxy protoclusters at $z=2-3$ discovered using the MAMMOTH technique. We introduce this survey by presenting the first measurement of the mass-metallicity relation (MZR) at high redshift in overdense environments via grism spectroscopy. From the completed MAMMOTH-Grism observations in the field of the BOSS1244 protocluster at $z=2.24\pm0.02$, We secure a sample of 36 protocluster member galaxies at $z\sim2.24$, showing strong nebular emission lines ([O III], H$β$ and [O II]) in their G141 spectra. Using the multi-wavelength broad-band deep imaging from HST and ground-based telescopes, we measure their stellar masses in the range of $[10^{9},10^{10.4}]M_\odot$, instantaneous star formation rates (SFR) from 10 to 240$M_\odot yr^{-1}$, and global gas-phase metallicities [$\frac{1}{3}$,1] of solar. Compared with similarly selected field galaxy sample at the same redshift, our galaxies show on average increased SFRs by $\sim$0.06dex and $\sim$0.18dex at $\sim$10$^{10.1}M_\odot$ and $\sim$10$^{9.8}M_\odot$, respectively. Using the stacked spectra of our sample galaxies, we derive the MZR in the BOSS1244 protocluster core as $12+\log({\rm O/H})=(0.136\pm0.018)\times\log(M_\ast/M_\odot)+(7.082\pm0.175)$, showing significantly shallower slope than that in the field. This shallow MZR slope is likely caused by the combined effects of efficient recycling of feedback-driven winds and cold-mode gas accretion in protocluster environments. The former effect helps low-mass galaxies residing in overdensities retain their metal production, whereas the latter effect dilutes the metal content of high-mass galaxies, making them more metal poor than their coeval field counterparts.

preprint2020arXiv

A Census of Sub-kiloparsec Resolution Metallicity Gradients in Star-forming Galaxies at Cosmic Noon from HST Slitless Spectroscopy

We present hitherto the largest sample of gas-phase metallicity radial gradients measured at sub-kiloparsec resolution in star-forming galaxies in the redshift range of $z\in[1.2, 2.3]$. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 76 galaxies with stellar mass ranging from 10$^7$ to 10$^{10}$ $M_\odot$, instantaneous star-formation rate in the range of [1, 100] $M_\odot$/yr, and global metallicity [$\frac{1}{12}$, 2] solar. At 2-$σ$ confidence level, 15/76 galaxies in our sample show negative radial gradients, whereas 7/76 show inverted gradients. Combining ours and all other metallicity gradients obtained at similar resolution currently available in the literature, we measure a negative mass dependence of $Δ\log({\rm O/H})/Δr~ [\mathrm{dex~kpc^{-1}}] = \left(-0.020\pm0.007\right) + \left(-0.016\pm0.008\right) \log(M_\ast/10^{9.4} M_\odot)$ with the intrinsic scatter being $σ=0.060\pm0.006$ over four orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemo-structural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star-formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of $M_\ast\lesssim10^9 M_\odot$.

preprint2020arXiv

Photometric Redshifts in the North Ecliptic Pole Wide Field based on a Deep Optical Survey with Hyper Suprime-Cam

The $AKARI$ space infrared telescope has performed near- to mid-infrared (MIR) observations on the North Ecliptic Pole Wide (NEPW) field (5.4 deg$^2$) for about one year. $AKARI$ took advantage of its continuous nine photometric bands, compared with NASA's $Spitzer$ and WISE space telescopes, which had only four filters with a wide gap in the MIR. The $AKARI$ NEPW field lacked deep and homogeneous optical data, limiting the use of nearly half of the IR sources for extra-galactic studies owing to the absence of photometric redshifts (photo-zs). To remedy this, we have recently obtained deep optical imaging over the NEPW field with 5 bands ($g$, $r$, $i$, $z$, and $Y$) of the Hyper Suprime-Camera (HSC) on the Subaru 8m telescope. We optically identify AKARI-IR sources along with supplementary $Spitzer$ and WISE data as well as pre-existing optical data. In this work, we derive new photo-zs using a $χ^2$ template-fitting method code ($Le$ $Phare$) and reliable photometry from 26 selected filters including HSC, $AKARI$, CFHT, Maidanak, KPNO, $Spitzer$ and WISE data. We take 2026 spectroscopic redshifts (spec-z) from all available spectroscopic surveys over the NEPW to calibrate and assess the accuracy of the photo-zs. At z < 1.5, we achieve a weighted photo-z dispersion of $σ_{Δ{z/(1+z)}}$ = 0.053 with $η$ = 11.3% catastrophic errors.

preprint2020arXiv

Search for Optically Dark Infrared Galaxies without Counterparts of Subaru Hyper Suprime-Cam in the AKARI North Ecliptic Pole Wide Survey Field

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $μ$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $μ$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

preprint2019arXiv

Discovery of Protoclusters at z~3.7 & 4.9: Embedded in Primordial Superclusters

We have carried out follow-up spectroscopy on three overdense regions of $g$- and $r$-dropout galaxies in the Canada-France-Hawaii Telescope Legacy Survey Deep Fields, finding two new protoclusters at $z=4.898$, 3.721 and a possible protocluster at $z=3.834$. The $z=3.721$ protocluster overlaps with a previously identified protocluster at $z=3.675$. The redshift separation between these two protoclusters is $Δz=0.05$, which is slightly larger than the size of typical protoclusters. Therefore, if they are not the progenitors of a $>10^{15}\,\mathrm{M_\odot}$ halo, they would grow into closely-located independent halos like a supercluster. The other protocluster at $z=4.898$ is also surrounded by smaller galaxy groups. These systems including protoclusters and neighboring groups are regarded as the early phase of superclusters. We quantify the spatial distribution of member galaxies of the protoclusters at $z=3.675$ and 3.721 by fitting triaxial ellipsoids, finding a tentative difference: one has a pancake-like shape while the other is filamentary. This could indicate that these two protoclusters are in different stages of formation. We investigate the relation between redshift and the velocity dispersion of protoclusters, including other protoclusters from the literature, in order to compare their dynamical states. Although there is no significant systematic trend in the velocity dispersions of protoclusters with redshift, the distribution is skewed to higher velocity dispersion over the redshift range of $z=2\mathrm{-}6$. This could be interpreted as two phases of cluster formation, one dominated by the steady accretion of galaxies, and the other by the merging between group-size halos, perhaps depending on the surrounding large-scale environments.

preprint2019arXiv

The mass relations between supermassive black holes and their host galaxies at 1<z<2 with HST-WFC3

Correlations between the mass of a supermassive black hole and the properties of its host galaxy (e.g., total stellar mass (M*), luminosity (Lhost)) suggest an evolutionary connection. A powerful test of a co-evolution scenario is to measure the relations MBH-Lhost and MBH-M* at high redshift and compare with local estimates. For this purpose, we acquired HST imaging with WFC3 of 32 X-ray-selected broad-line AGN at 1.2<z<1.7 in deep survey fields. By applying state-of-the-art tools to decompose the HST images including available ACS data, we measured the host galaxy luminosity and stellar mass along with other properties through the 2D model fitting. The black hole mass was determined using the broad Halpha line, detected in the near-infrared with Subaru/FMOS, which potentially minimizes systematic effects using other indicators. We find that the observed ratio of MBH to total M* is 2.7 times larger at z~1.5 than in the local universe, while the scatter is equivalent between the two epochs. A non-evolving mass ratio is consistent with the data at the 2-3 sigma confidence level when accounting for selection effects and their uncertainties. The relationship between MBH-Lhost paints a similar picture. Therefore, our results cannot distinguish whether SMBHs and their total M* and Lhost proceed in lockstep or whether the growth of the former somewhat overshoots the latter, given the uncertainties. Based on a statistical estimate of the bulge-to-total mass fraction, the ratio MBH/M* is offset from the local value by a factor of ~7 which is significant even accounting for selection effects. Taken together, these observations are consistent with a scenario in which stellar mass is subsequently transferred from an angular momentum supported component of the galaxy to the pressure supported one through secular processes or minor mergers at a faster rate than mass accretion onto the SMBH.

preprint2019arXiv

The Seoul National University AGN monitoring project. II. BLR size and black hole mass of two AGNs

Active galactic nuclei (AGNs) show a correlation between the size of the broad line region (BLR) and the monochromatic continuum luminosity at 5100 Å, allowing black hole mass estimation based on single-epoch spectra. However, the validity of the correlation is yet to be clearly tested for high-luminosity AGNs. We present the first reverberation-mapping results of the Seoul National University AGN monitoring program (SAMP), which is designed to focus on luminous AGNs for probing the high end of the size-luminosity relation. We report time lag measurements of two AGNs, namely, 2MASS J10261389+5237510 and SDSS J161911.24+501109.2, using the light curves obtained over a $\sim$1000 day period with an average cadence of $\sim$10 and $\sim$20 days, respectively for photometry and spectroscopy monitoring. Based on a cross-correlation analysis and H$β$ line width measurements, we determine the H$β$ lag as $41.8^{+4.9}_{-6.0}$ and $52.6^{+17.6}_{-14.7}$ days in the observed-frame, and black hole mass as $3.65^{+0.49}_{-0.57} \times 10^7 M_{\odot}$ and $23.02^{+7.81}_{-6.56} \times 10^7 M_{\odot}$, respectively for 2MASS J1026 and SDSS J1619.

preprint2016arXiv

A Systematic Survey of Protoclusters at $z\sim3\mathrm{-}6$ in the CFHTLS Deep Fields

We present the discovery of three protoclusters at $z\sim3\mathrm{-}4$ with spectroscopic confirmation in the Canada-France-Hawaii Telescope (CFHT) Legacy Survey Deep Fields. In these fields, we investigate the large-scale projected sky distribution of $z\sim3\mathrm{-}6$ Lyman break galaxies and identify 21 protocluster candidates from regions that are overdense at more than $4σ$ overdensity significance. Based on cosmological simulations, it is expected that more than $76\%$ of these candidates will evolve into a galaxy cluster of at least a halo mass of $10^{14}\,\mathrm{M_\odot}$ at $z=0$. We perform follow-up spectroscopy for eight of the candidates using Subaru/FOCAS, KeckII/DEIMOS, and Gemini-N/GMOS. In total we target 462 dropout candidates and obtain 138 spectroscopic redshifts. We confirm three real protoclusters at $z=3\mathrm{-}4$ with more than five members spectroscopically identified, and find one to be an incidental overdense region by mere chance alignment. The other four candidate regions at $z\sim5\mathrm{-}6$ require more spectroscopic follow-up in order to be conclusive. A $z=3.67$ protocluster, which has eleven spectroscopically confirmed members, shows a remarkable core-like structure composed of a central small region ($<0.5\,\mathrm{physical\>Mpc}$) and an outskirts region ($\sim1.0\,\mathrm{physical\>Mpc}$). The Ly$α$ equivalent widths of members of the protocluster are significantly smaller than those of field galaxies at the same redshift while there is no difference in the UV luminosity distributions. These results imply that some environmental effects start operating as early as at $z\sim4$ along with the growth of the protocluster structure.

preprint2016arXiv

An Empirical Determination of the Intergalactic Background Light from UV to FIR Wavelengths Using FIR Deep Galaxy Surveys and the Gamma-ray Opacity of the Universe

We have previously calculated the intergalactic background light (IBL) as a function of redshift in the far ultraviolet to near infrared range, based purely on data from deep galaxy surveys. Here we utilize similar methods to determine the mid- and far infrared IBL out to a wavelength of 850 microns. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. By also including the effect of the 2.7 K cosmic background photons, we determine 68% confidence upper and lower limits on the opacity of the universe to gamma-rays up to PeV energies. Our direct results on the IBL are consistent with those from complimentary gamma-ray analyses using observations from the Fermi $γ$-ray space telescope and the H.E.S.S. air Cherenkov telescope. Thus, we find no evidence of previously suggested processes for the modification of gamma-ray spectra other than that of absorption by pair production alone.

preprint2016arXiv

Be Stars in the Open Cluster NGC 6830

We report the discovery of 2 new Be stars, and re-identify one known Be star in the open cluster NGC 6830. Eleven H-alpha emitters were discovered using the H-alpha imaging photometry of the Palomar Transient Factory Survey. Stellar membership of the candidates was verified with photometric and kinematic information using 2MASS data and proper motions. The spectroscopic confirmation was carried out by using the Shane 3-m telescope at Lick observatory. Based on their spectral types, three H-alpha emitters were confirmed as Be stars with H-alpha equivalent widths > -10 Angstrom. Two objects were also observed by the new spectrograph SED-Machine on the Palomar 60 inch Telescope. The SED-Machine results show strong H-alpha emission lines, which are consistent with the results of the Lick observations. The high efficiency of the SED-Machine can provide rapid observations for Be stars in a comprehensive survey in the future.

preprint2016arXiv

Broad Hbeta Emission-Line Variability in a Sample of 102 Local Active Galaxies

A sample of 102 local (0.02 < z < 0.1) Seyfert galaxies with black hole masses MBH > 10^7 M_sun was selected from the Sloan Digital Sky Survey (SDSS) and observed using the Keck 10-m telescope to study the scaling relations between MBH and host galaxy properties. We study profile changes of the broad Hbeta emission line within the ~3-9 year time-frame between the two sets of spectra. The variability of the broad Hbeta emission line is of particular interest, not only since it is used to estimate MBH, but also since its strength and width is used to classify Seyfert galaxies into different types. At least some form of broad-line variability (in either width or flux) is observed in the majority (~66%) of the objects, resulting in a Seyfert-type change for ~38% of the objects, likely driven by variable accretion and/or obscuration. The broad Hbeta line virtually disappears in 3/102 (~3%) extreme cases. We discuss potential causes for these changing-look AGNs. While similar dramatic transitions have previously been reported in the literature, either on a case-by-case basis or in larger samples focusing on quasars at higher redshifts, our study provides statistical information on the frequency of H$β$ line variability in a sample of low-redshift Seyfert galaxies.

preprint2016arXiv

The Metal Abundances across Cosmic Time ($\mathcal{MACT}$) Survey. I. Optical Spectroscopy in the Subaru Deep Field

Deep rest-frame optical spectroscopy is critical for characterizing and understanding the physical conditions and properties of the ionized gas in galaxies. Here, we present a new spectroscopic survey called "Metal Abundances across Cosmic Time" or $\mathcal{MACT}$, which will obtain rest-frame optical spectra for $\sim$3000 emission-line galaxies. This paper describes the optical spectroscopy that has been conducted with MMT/Hectospec and Keck/DEIMOS for $\approx$1900 $z=0.1-1$ emission-line galaxies selected from our narrowband and intermediate-band imaging in the Subaru Deep Field. In addition, we present a sample of 164 galaxies for which we have measured the weak [OIII]$λ$4363 line (66 with at least 3$σ$ detections and 98 with significant upper limits). This nebular emission line determines the gas-phase metallicity by measuring the electron temperature of the ionized gas. This paper presents the optical spectra, emission-line measurements, interstellar properties (e.g., metallicity, gas density), and stellar properties (e.g., star formation rates, stellar mass). Paper II of the $\mathcal{MACT}$ survey (Ly et al.) presents the first results on the stellar mass--gas metallicity relation at $z\lesssim1$ using the sample with [OIII]$λ$4363 measurements.

preprint2016arXiv

The Metal Abundances across Cosmic Time ($\mathcal{MACT}$) Survey. II. Evolution of the Mass-Metallicity Relation over 8 Billion Years, using [OIII]$λ$4363Å-based Metallicities

We present the first results from MMT and Keck spectroscopy for a large sample of $0.1\leq z\leq1$ emission-line galaxies selected from our narrow-band imaging in the Subaru Deep Field. We measured the weak [OIII]$λ$4363 emission line for 164 galaxies (66 with at least 3$σ$ detections, and 98 with significant upper limits). The strength of this line is set by the electron temperature for the ionized gas. Because the gas temperature is regulated by the metal content, the gas-phase oxygen abundance is inversely correlated with [OIII]$λ$4363 line strength. Our temperature-based metallicity study is the first to span $\approx$8 Gyr of cosmic time and $\approx$3 dex in stellar mass for low-mass galaxies, $\log{\left(M_{\rm star}/M_{\rm sun}\right)}\approx6.0-9.0$. Using extensive multi-wavelength photometry, we measure the evolution of the stellar mass--gas metallicity relation and its dependence on dust-corrected star formation rate (SFR). The latter is obtained from high signal-to-noise Balmer emission-line measurements. Our mass-metallicity relation is consistent with Andrews & Martini at $z\leq0.3$, and evolves toward lower abundances at a given stellar mass, $\log{({\rm O/H})}\propto(1+z)^{-2.32^{+0.52}_{-0.26}}$. We find that galaxies with lower metallicities have higher SFRs at a given stellar mass and redshift, although the scatter is large ($\approx$0.3 dex), and the trend is weaker than seen in local studies. We also compare our mass--metallicity relation against predictions from high-resolution galaxy formation simulations, and find good agreement with models that adopt energy- and momentum-driven stellar feedback. We have identified 16 extremely metal-poor galaxies with abundances less than a tenth of solar; our most metal-poor galaxy at $z\approx0.84$ is similar to I Zw 18.

preprint2016arXiv

The Optical Variability of SDSS Quasars from Multi-epoch Spectroscopy. III. A Sudden UV Cutoff in Quasar SDSS J2317+0005

We have collected near-infrared to X-ray data of 20 multi-epoch heavily reddened SDSS quasars to investigate the physical mechanism of reddening. Of these, J2317+0005 is found to be a UV cutoff quasar. Its continuum, which usually appears normal, decreases by a factor 3.5 at 3000Å, compared to its more typical bright state during an interval of 23 days. During this sudden continuum cut-off, the broad emission line fluxes do not change, perhaps due to the large size of the Broad Line Region (BLR), r > 23 / (1+z) days. The UV continuum may have suffered a dramatic drop out. However, there are some difficulties with this explanation. Another possibility is that the intrinsic continuum did not change, but was temporarily blocked out, at least towards our line of sight. As indicated by X-ray observations, the continuum rapidly recovers after 42 days. A comparison of the bright state and dim states would imply an eclipse by a dusty cloud with a reddening curve having a remarkably sharp rise shortward of 3500Å. Under the assumption of being eclipsed by a Keplerian dusty cloud, we characterized the cloud size with our observations, however, which is a little smaller than the 3000Å continuum-emitting size inferred from accretion disk models. Therefore, we speculate this is due to a rapid outflow or inflow with a dusty cloud passing through our line-of-sight to the center.

preprint2015arXiv

A Local Baseline of the Black Hole Mass Scaling Relations for Active Galaxies. III. The BH mass - $σ$ relation

We create a baseline of the black hole (BH) mass (MBH) - stellar-velocity dispersion (sigma) relation for active galaxies, using a sample of 66 local (0.02<z<0.09) Seyfert-1 galaxies, selected from the Sloan Digital Sky Survey (SDSS). Analysis of SDSS images yields AGN luminosities free of host-galaxy contamination and morphological classification. 51/66 galaxies have spiral morphology. 28 bulges have Sersic index n<2 and are considered candidate pseudo bulges, with eight being definite pseudo bulges based on multiple classification criteria met. Only 4/66 galaxies show sign of interaction/merging. High signal-to-noise ratio Keck spectra provide the width of the broad Hbeta emission line free of FeII emission and stellar absorption. AGN luminosity and Hbeta line widths are used to estimate MBH. The Keck-based spatially-resolved kinematics is used to determine stellar-velocity dispersion within the spheroid effective radius. We find that sigma can vary on average by up to 40% across definitions commonly used in the literature, emphasizing the importance of using self-consistent definitions in comparisons and evolutionary studies. The MBH-sigma relation for our Seyfert-1 galaxies has the same intercept and scatter as that of reverberation-mapped AGNs as well as quiescent galaxies, consistent with the hypothesis that our single epoch MBH estimator and sample selection do not introduce significant biases. Barred galaxies, merging galaxies, and those hosting pseudo bulges do not represent outliers in the MBH-sigma relation. This is in contrast with previous work, although no firm conclusion can be drawn due to the small sample size and limited resolution of the SDSS images.

preprint2015arXiv

Multiple Images of a Highly Magnified Supernova Formed by an Early-Type Cluster Galaxy Lens

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

preprint2015arXiv

Physical conditions of the interstellar medium in star-forming galaxies at z~1.5

We present results from Subaru/FMOS near-infrared (NIR) spectroscopy of 118 star-forming galaxies at $z\sim1.5$ in the Subaru Deep Field. These galaxies are selected as [OII]$λ$3727 emitters at $z\approx$ 1.47 and 1.62 from narrow-band imaging. We detect H$α$ emission line in 115 galaxies, [OIII]$λ$5007 emission line in 45 galaxies, and H$β$, [NII]$λ$6584, and [SII]$λλ$6716,6731 in 13, 16, and 6 galaxies, respectively. Including the [OII] emission line, we use the six strong nebular emission lines in the individual and composite rest-frame optical spectra to investigate physical conditions of the interstellar medium in star-forming galaxies at $z\sim$1.5. We find a tight correlation between H$α$ and [OII], which suggests that [OII] can be a good star formation rate (SFR) indicator for galaxies at $z\sim1.5$. The line ratios of H$α$/[OII] are consistent with those of local galaxies. We also find that [OII] emitters have strong [OIII] emission lines. The [OIII]/[OII] ratios are larger than normal star-forming galaxies in the local Universe, suggesting a higher ionization parameter. Less massive galaxies have larger [OIII]/[OII] ratios. With evidence that the electron density is consistent with local galaxies, the high ionization of galaxies at high redshifts may be attributed to a harder radiation field by a young stellar population and/or an increase in the number of ionizing photons from each massive star.

preprint2014arXiv

An Empirical Determination of the Intergalactic Background Light Using NIR Deep Galaxy Survey Data out to 5 microns and the Gamma-ray Opacity of the Universe

We extend our previous model-independent determination of the intergalactic background light (IBL), based purely on galaxy survey data, out to a wavelength of 5 microns. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. We further determine a 68% confidence upper and lower limit on the opacity of the universe to gamma-rays up to energies of 1.6/(1+z) TeV. A comparison of our lower limit redshift-dependent opacity curves to the opacity limits derived from the results of both ground-based air Cherenkov telescope and Fermi-LAT observations of PKS 1424+240 allows us to place a new upper limit on the redshift of this source, independent of IBL modeling.

preprint2014arXiv

Cosmic Evolution of Black Holes and Spheroids. V. The Relation Between Black Hole Mass and Host Galaxy Luminosity for a Sample of 79 Active Galaxies

We investigate the cosmic evolution of the black hole (BH) mass -- bulge luminosity relation using a sample of 52 active galaxies at $z \sim 0.36$ and $z \sim 0.57$ in the BH mass range of $10^{7.4-9.1} M_{\odot}$. By consistently applying multi-component spectral and structural decomposition to high-quality Keck spectra and high-resolution HST images, BH masses ($M_{\rm BH}$) are estimated using the H$β$ broad emission line combined with the 5100 Å nuclear luminosity, and bulge luminosities ($L_{\rm bul}$) are derived from surface photometry. Comparing the resulting $M_{\rm BH}-L_{\rm bul}$ relation to local active galaxies and taking into account selection effects, we find evolution of the form $M_{\rm BH} / L_{\rm bul} \propto (1+z)^γ$ with $γ=1.8\pm0.7$, consistent with BH growth preceding that of the host galaxies. Including an additional sample of 27 active galaxies with $0.5<z<1.9$ taken from the literature and measured in a consistent way, we obtain $γ=0.9\pm0.7$ for the $M_{\rm BH}-L_{\rm bul}$ relation and $γ=0.4\pm0.5$ for the $M_{\rm BH}$--total host galaxy luminosity ($L_{\rm host}$) relation. The results strengthen the findings from our previous studies and provide additional evidence for host-galaxy bulge growth being dominated by disk-to-bulge transformation via minor mergers and/or disk instabilities.

preprint2014arXiv

Far-infrared line spectra of Seyfert galaxies from the Herschel-PACS Spectrometer

We present spectroscopic observations of FIR fine-structure lines of 26 Seyfert galaxies obtained with the Herschel-PACS spectrometer. These observations are complemented by spectroscopy with Spitzer-IRS and Herschel-SPIRE. The ratios of the OIII, NII, SIII and NeV lines have been used to determine electron densities in the ionised gas regions. The CI lines, observed with SPIRE, have been used to measure the densities in the neutral gas, while the OI lines provide a measure of the gas temperature, at densities below 10000 cm-3. Using the OI145/63um and SIII33/18um line ratios we find an anti-correlation of the temperature with the gas density. Using various fine-structure line ratios, we find that density stratification is common in these active galaxies. On average, the electron densities increase with the ionisation potential of the ions producing the NII, SIII and NeV emission. The infrared emission lines arise partly in the Narrow Line Region (NLR) photoionised by the AGN central engine, partly in HII regions photo ionised by hot stars and partly in neutral gas in photo-dissociated regions (PDRs). We attempt to separate the contributions to the line emission produced in these different regions by comparing our emission line ratios to empirical and theoretical values. In particular, we tried to separate the contribution of AGN and star formation by using a combination of Spitzer and Herschel lines, and we found that, besides the well known mid-IR line ratios, the mixed mid-IR/far-IR line ratio of OIII88um/OIV26um can reliably discriminate the two emission regimes, while the far-IR line ratio of CII157um/OI63um is only able to mildly separate the two regimes. By comparing the observed CII157um/NII205um ratio with photoionisation models, we also found that most of the CII emission in the galaxies we examined is due to PDRs.

preprint2014arXiv

Luminosity and Redshift Dependence of the Covering Factor of AGNs viewed with WISE and SDSS

In this work, we investigate the dependence of the covering factor (CF) of active galactic nuclei (AGNs) (i) on the mid-infrared (MIR) luminosity and (ii) on the redshift. We constructed 12- and 22-micron luminosity functions (LFs) at 0.006 < z < 0.3 using the Wide-field Infrared Survey Explorer} (WISE) data. Combining the WISE catalog with the Sloan Digital Sky Survey (SDSS) spectroscopic data, we selected 223,982 galaxies at 12 micron and 25,721 galaxies at 22 micron for spectroscopic classification. We then identified 16,355 AGNs at 12 micron and 4,683 AGNs at 22 micron by their optical emission lines and cataloged classifications in the SDSS. Following that, we estimated the CF as the fraction of type 2 AGN in all AGNs whose MIR emissions are dominated by the active nucleus (not their host galaxies) based on their MIR colors. We found that (i) the CF decreased with increasing MIR luminosity, regardless of the choice of type 2 AGN classification criteria, and (ii) the CF did not change significantly with the redshift for z < 0.2. Furthermore, we carried out various tests to determine the influence of selection bias and confirmed similar dependences exist even when taking these uncertainties into account. The luminosity dependence of the CF can be explained by the receding torus model, but the "modified" receding torus model gives a slightly better fit, as suggested by Simpson.

preprint2013arXiv

"Direct" Gas-phase Metallicities, Stellar Properties, and Local Environments of Emission-line Galaxies at Redshift below 0.90

Using deep narrow-band (NB) imaging and optical spectroscopy from the Keck telescope and MMT, we identify a sample of 20 emission-line galaxies (ELGs) at z=0.065-0.90 where the weak auroral emission line, [OIII]4363, is detected at >3σ. These detections allow us to determine the gas-phase metallicity using the "direct'' method. With electron temperature measurements and dust attenuation corrections from Balmer decrements, we find that 4 of these low-mass galaxies are extremely metal-poor with 12+log(O/H) <= 7.65 or one-tenth solar. Our most metal-deficient galaxy has 12+log(O/H) = 7.24^{+0.45}_{-0.30} (95% confidence), similar to some of the lowest metallicity galaxies identified in the local universe. We find that our galaxies are all undergoing significant star formation with average specific star formation rate (SFR) of (100 Myr)^{-1}, and that they have high central SFR surface densities (average of 0.5 Msun/yr/kpc^2. In addition, more than two-thirds of our galaxies have between one and four nearby companions within a projected radius of 100 kpc, which we find is an excess among star-forming galaxies at z=0.4-0.85. We also find that the gas-phase metallicities for a given stellar mass and SFR lie systematically below the local M-Z-(SFR) relation by \approx0.2 dex (2σ significance). These results are partly due to selection effects, since galaxies with strong star formation and low metallicity are more likely to yield [OIII]4363 detections. Finally, the observed higher ionization parameter and electron density suggest that they are lower redshift analogs to typical z>1 galaxies.

preprint2013arXiv

Discovery and Early Multi-Wavelength Measurements of the Energetic Type Ic Supernova PTF12gzk: A Massive-Star Explosion in a Dwarf Host Galaxy

We present the discovery and extensive early-time observations of the Type Ic supernova (SN) PTF12gzk. Our finely sampled light curves show a rise of 0.8mag within 2.5hr. Power-law fits [f(t)\sim(t-t_0)^n] to these data constrain the explosion date to within one day. We cannot rule out the expected quadratic fireball model, but higher values of n are possible as well for larger areas in the fit parameter space. Our bolometric light curve and a dense spectral sequence are used to estimate the physical parameters of the exploding star and of the explosion. We show that the photometric evolution of PTF12gzk is slower than that of most SNe Ic, and its high ejecta velocities (~30,000km/s four days after explosion) are closer to the observed velocities of broad-lined SNe Ic associated with gamma-ray bursts (GRBs) than to the observed velocities in normal Type Ic SNe. The high velocities are sustained through the SN early evolution, and are similar to those of GRB-SNe when the SN reach peak magnitude. By comparison with the spectroscopically similar SN 2004aw, we suggest that the observed properties of PTF12gzk indicate an initial progenitor mass of 25-35 solar mass and a large (5-10E51 erg) kinetic energy, close to the regime of GRB-SN properties. The host-galaxy characteristics are consistent with GRB-SN hosts, and not with normal SN Ic hosts as well, yet this SN does not show the broad lines over extended periods of time that are typical of broad-line Type Ic SNe.

preprint2013arXiv

Fueling AGN-I: How the Global Characteristics of the Central Kiloparsec of Seyferts differ from Quiescent Galaxies

In a matched sample of Seyfert and quiescent galaxies we simultaneously probe the stellar and molecular gas kinematics from 1 kpc down to 50 pc with the aim of identifying the dynamical processes dictating black hole accretion rates. This first paper compares the global characteristics of a sample of ten galaxies. We find several differences within a radius of 500 pc that are correlated with Active Galactic Nucleus (AGN) activity. The Seyferts have: (1) a more centrally concentrated nuclear stellar surface brightness with lower stellar luminosities beyond a radius of 100 pc, (2) a lower stellar velocity dispersion within a radius of 200 pc, (3) elevated H2 luminosity out to a radius of at least 250 pc, and (4) more centrally concentrated H2 surface brightness profiles. These observed differences can be interpreted as evidence for Seyfert galaxies having a dynamically cold (in comparison to the bulge) nuclear structure composed of a significant gas reservoir and a relatively young stellar population. This structure is undetected (and possibly does not exist) in quiescent galaxies. The presence of such a nuclear structure in Seyfert galaxies provides evidence for inflow of the surrounding interstellar medium since the nuclear stellar population requires a supply of gas from which to form. The fueling of a Seyfert AGN is thus associated with the formation of a dynamically cold component of gas and stars on scales of hundreds of parsecs.

preprint2013arXiv

The 9 and 18 Micron Luminosity Function of Various Types of Galaxies with AKARI: Implication for the Dust Torus Structure of AGN

We present the 9 and 18 micron luminosity functions (LFs) of galaxies at 0.006 < z < 0.8 (with an average redshift of ~ 0.04) using the AKARI mid-infrared all-sky survey catalog. We selected 243 galaxies at 9 micron and 255 galaxies at 18 micron from the Sloan Digital Sky Survey (SDSS) spectroscopy region. These galaxies were then classified by their optical emission lines, such as the line width of H_alpha or by their emission line ratios of [OIII]/H_beta and [NII]/H_alpha into five types: Type 1 active galactic nuclei (AGN) (Type 1); Type 2 AGN (Type 2); low-ionization narrow emission line galaxies (LINER); galaxies with both star formation and narrow-line AGN activity (composite galaxies); and star-forming galaxies (SF). We found that (i) the number density ratio of Type 2 to Type 1 AGNs is 1.73 +/- 0.36, which is larger than a result obtained from the optical LF and (ii) this ratio decreases with increasing 18 micron luminosity.

preprint2013arXiv

The Lick AGN Monitoring Project 2011: Fe II Reverberation from the Outer Broad-Line Region

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including H-beta and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases the measured lag (t_cen) of Fe II is longer than that of H-beta, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al. (2004), we find t_cen(Fe II)/t_cen(H-beta) = 1.9+-0.6 in NGC 4593 and 1.5+-0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region.

preprint2013arXiv

The Low-Luminosity End of the Radius-Luminosity Relationship for Active Galactic Nuclei

We present an updated and revised analysis of the relationship between the Hbeta broad-line region (BLR) radius and the luminosity of the active galactic nucleus (AGN). Specifically, we have carried out two-dimensional surface brightness decompositions of the host galaxies of 9 new AGNs imaged with the Hubble Space Telescope Wide Field Camera 3. The surface brightness decompositions allow us to create "AGN-free" images of the galaxies, from which we measure the starlight contribution to the optical luminosity measured through the ground-based spectroscopic aperture. We also incorporate 20 new reverberation-mapping measurements of the Hbeta time lag, which is assumed to yield the average Hbeta BLR radius. The final sample includes 41 AGNs covering four orders of magnitude in luminosity. The additions and updates incorporated here primarily affect the low-luminosity end of the R-L relationship. The best fit to the relationship using a Bayesian analysis finds a slope of alpha = 0.533 (+0.035/-0.033), consistent with previous work and with simple photoionization arguments. Only two AGNs appear to be outliers from the relationship, but both of them have monitoring light curves that raise doubt regarding the accuracy of their reported time lags. The scatter around the relationship is found to be 0.19(+/-0.02) dex, but would be decreased to 0.13 dex by the removal of these two suspect measurements. A large fraction of the remaining scatter in the relationship is likely due to the inaccurate distances to the AGN host galaxies. Our results help support the possibility that the R-L relationship could potentially be used to turn the BLRs of AGNs into standardizable candles. This would allow the cosmological expansion of the Universe to be probed by a separate population of objects, and over a larger range of redshifts.

preprint2012arXiv

A Determination of the Intergalactic Redshift Dependent UV-Optical-NIR Photon Density Using Deep Galaxy Survey Data and the Gamma-ray Opacity of the Universe

We calculate the intensity and photon spectrum of the intergalactic background light (IBL) as a function of redshift using an approach based on observational data obtained in many different wavelength bands from local to deep galaxy surveys. This allows us to obtain an empirical determination of the IBL and to quantify its observationally based uncertainties. Using our results on the IBL, we then place 68% confidence upper and lower limits on the opacity of the universe to gamma-rays, free of the theoretical assumptions that were needed for past calculations. We compare our results with measurements of the extragalactic background light and upper limits obtained from observations made by the Fermi Gamma-ray Space Telescope.

preprint2012arXiv

A Ly-alpha Emitter with an Extremely Large Rest-frame Equivalent Width of ~900A at z=6.5: A Candidate of Population III-dominated Galaxy?

We have identified a very interesting Ly-alpha emitter, whose Ly-alpha emission line has an extremely large observed equivalent width of EW_0=436^{+422}_{-149}A, which corresponds to an extraordinarily large intrinsic rest-frame equivalent width of EW_0^{int}=872^{+844}_{-298}A after the average intergalactic absorption correction. The object was spectroscopically confirmed to be a real Ly-alpha emitter by its apparent asymmetric Ly-alpha line profile detected at z=6.538. The continuum emission of the object was definitely detected in our deep z'-band image; thus, its EW_0 was reliably determined. Follow-up deep near-infrared spectroscopy revealed emission lines of neither He II lambda1640 as an apparent signature of Population III, nor C IV lambda1549 as a proof of active nucleus. No detection of short-lived He II lambda1640 line is not necessarily inconsistent with the interpretation that the underlying stellar population of the object is dominated by Population III. We found that the observed extremely large EW_0 of the Ly-alpha emission and the upper limit on the EW_0 of the He II lambda1640 emission can be explained by population synthesis models favoring a very young age less than 2-4Myr and massive metal-poor (Z<10^{-5}) or even metal-free stars. The observed large EW_0 of Ly-alpha is hardly explained by Population I/II synthesis models with Z>10^{-3}. However, we cannot conclusively rule out the possibility that this object is composed of a normal stellar population with a clumpy dust distribution, which could enhance the Ly-alpha EW_0, though its significance is still unclear.

preprint2012arXiv

Discovery of a protocluster at z~6

We report the discovery of a protocluster at z~6 containing at least eight cluster member galaxies with spectroscopic confirmations in the wide-field image of the Subaru Deep Field (SDF). The overdensity of the protocluster is significant at the 6 sigma level, based on the surface number density of i'-dropout galaxies. The overdense region covers ~36 sq. arcmin, and includes 30 i'-dropout galaxies. Follow-up spectroscopy revealed that 15 of these are real z~6 galaxies (5.7 < z < 6.3). Eight of the 15 are clustering in a narrow redshift range centered at z=6.01, corresponding to a seven-fold increase in number density over the average in redshift space. We found no significant difference in the observed properties, such as Ly-alpha luminosities and UV continuum magnitudes, between the eight protocluster members and the seven non-members. The velocity dispersion of the eight protocluster members is 647 km/s, which is about three times higher than that predicted by the standard cold dark matter model. This discrepancy could be attributed to the distinguishing three-dimensional distribution of the eight protocluster members. We discuss two possible explanations for this discrepancy: either the protocluster is already mature, with old galaxies at the center, or it is still immature and composed of three subgroups merging to become a larger cluster. In either case, this concentration of z=6.01 galaxies in the SDF may be one of the first sites of formation of a galaxy cluster in the universe.

preprint2012arXiv

Dust Attenuation and H-alpha Star Formation Rates of z~0.5 Galaxies

Using deep narrow-band and broad-band imaging, we identify 401 z~0.40 and 249 z~0.49 H-alpha line-emitting galaxies in the Subaru Deep Field. Compared to other H-alpha surveys at similar redshifts, our samples are unique since they probe lower H-alpha luminosities, are augmented with multi-wavelength (rest-frame 1000AA--1.5 microns) coverage, and a large fraction (20%) of our samples has already been spectroscopically confirmed. Our spectra allow us to measure the Balmer decrement for nearly 60 galaxies with H-beta detected above 5-sigma. The Balmer decrements indicate an average extinction of A(H-alpha)=0.7^{+1.4}_{-0.7} mag. We find that the Balmer decrement systematically increases with higher H-alpha luminosities and with larger stellar masses, in agreement with previous studies with sparser samples. We find that the SFRs estimated from modeling the spectral energy distribution (SED) is reliable---we derived an "intrinsic" H-alpha luminosity which is then reddened assuming the color excess from SED modeling. The SED-predicted H-alpha luminosity agrees with H-alpha narrow-band measurements over 3 dex (rms of 0.25 dex). We then use the SED SFRs to test different statistically-based dust corrections for H-alpha and find that adopting one magnitude of extinction is inappropriate: galaxies with lower luminosities are less reddened. We find that the luminosity-dependent dust correction of Hopkins et al. yields consistent results over 3 dex (rms of 0.3 dex). Our comparisons are only possible by assuming that stellar reddening is roughly half of nebular reddening. The strong correspondence argue that with SED modeling, we can derive reliable intrinsic SFRs even in the absence of H-alpha measurements at z~0.5.

preprint2012arXiv

Reliable Identifications of AGN from the WISE, 2MASS and Rosat all-sky surveys

We have developed the "S_IX" statistic to identify bright, highly-likely Active Galactic Nucleus (AGN) candidates solely on the basis of WISE, 2MASS and Rosat all-sky survey data. This statistic was optimized with data from the preliminary WISE survey and the SDSS, and tested with Lick 3-m Kast spectroscopy. We find that sources with S_IX < 0 have a <~95% likelihood of being an AGN (defined in this paper as a Seyfert 1, quasar or blazar). This statistic was then applied to the full WISE/2MASS/RASS dataset, including the final WISE data release, to yield the "W2R" sample of 4,316 sources with S_IX < 0. Only 2,209 of these sources are currently in the Veron-Cetty and Veron (VCV) catalog of spectroscopically confirmed AGN, indicating that the W2R sample contains nearly 2,000 new, relatively bright (J <~ 16) AGN. We utilize the W2R sample to quantify biases and incompleteness in the VCV catalog. We find it is highly complete for bright (J < 14), northern AGN, but the completeness drops below 50% for fainter, southern samples and for sources near the Galactic plane. This approach also led to the spectroscopic identification of 10 new AGN in the Kepler field, more than doubling the number of AGN being monitored by Kepler. This has identified ~1 bright AGN every 10 square degrees, permitting construction of AGN samples in any sufficiently large region of sky.

preprint2012arXiv

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad Line Region in Mrk 50

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-σ* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

preprint2012arXiv

The Stellar Population and Star Formation Rates of z~1.5-1.6 [O II] Emitting Galaxies Selected from Narrow-Band Emission-Line Surveys

We present the first detailed study of the stellar populations of star-forming galaxies at z~1.5, which are selected by their [O II] emission line, detected in narrow-band surveys. We identified ~1,300 [O II] emitters at z=1.47 and z=1.62 in the Subaru Deep Field with rest-frame EWs above 13Å. Optical and near-infrared spectroscopic observations for ~10% of our samples show that our separation of [O II] from [O III] emission-line galaxies in two-color space is 99% successful. We analyze the multi-wavelength properties of a subset of ~1,200 galaxies with the best photometry. They have average rest-frame EW of 45Å, stellar mass of 3 x 10^9 M_sun, and stellar age of 100 Myr. In addition, our SED fitting and broad-band colors indicate that [O II] emitters span the full range of galaxy populations at z~1.5. We also find that 80% of [O II] emitters are also photometrically classified as "BX/BM" (UV) galaxies and/or the star-forming "BzK" (near-IR) galaxies. Our [O II] emission line survey produces a far more complete, and somewhat deeper sample of z~1.5 galaxies than either the BX/BM or sBzK selection alone. We constructed average SEDs and find that higher [O II] EW galaxies have somewhat bluer continua. SED model-fitting shows that they have on average half the stellar mass of galaxies with lower [O II] EW. The observed [O II] luminosity is well-correlated with the far-UV continuum with a logarithmic slope slightly 0f 0.89\pm0.22. The scatter of the [O II] luminosity against the far-UV continuum suggests that [O II] can be used as a SFR indicator with a reliability of 0.23 dex.

preprint2012arXiv

Unusual Long and Luminous Optical Transient in the Subaru Deep Field

We present observations of SDF-05M05, an unusual optical transient discovered in the Subaru Deep Field (SDF). The duration of the transient is > ~800 d in the observer frame, and the maximum brightness during observation reached approximately 23 mag in the i' and z' bands. The faint host galaxy is clearly identified in all 5 optical bands of the deep SDF images. The photometric redshift of the host yields z~0.6 and the corresponding absolute magnitude at maximum is ~-20. This implies that this event shone with an absolute magnitude brighter than -19 mag for approximately 300 d in the rest frame, which is significantly longer than a typical supernova and ultra-luminous supernova. The total radiated energy during our observation was 1x10^51 erg. The light curves and color evolution are marginally consistent with some of luminous IIn supernova. We suggest that the transient may be a unique and peculiar supernova at intermediate redshift.

preprint2011arXiv

A Census of Star-Forming Galaxies at z = 1-3 in the Subaru Deep Field

Several UV and near-infrared color selection methods have identified galaxies at z = 1-3. Since each method suffers from selection biases, we have applied three leading techniques (Lyman break, BX/BM, and BzK selection) simultaneously in the Subaru Deep Field. This field has reliable (Δz/(1 + z) = 0.02--0.09) photometric redshifts for ~53,000 galaxies from 20 bands (1500Å--2.2μm). The BzK, LBG, and BX/BM samples suffer contamination from z<1 interlopers of 6%, 8%, and 20%, respectively. Around the redshifts where it is most sensitive (z~1.9 for star-forming BzK, z~1.8 for z~2 LBGs, z~1.6 for BM, and z~2.3 for BX), each technique finds 60-80% of the census of the three methods. In addition, each of the color techniques shares 75-96% of its galaxies with another method, which is consistent with previous studies that adopt identical criteria on magnitudes and colors. Combining the three samples gives a comprehensive census that includes ~90% of z-phot = 1-3 galaxies, using standard magnitude limits similar to previous studies. In fact, we find that among z = 1-2.5 galaxies in the color selection census, 81-90% of them can be selected by just combining the BzK selection with one of the UV techniques (z~2 LBG or BX and BM). The average galaxy stellar mass, reddening and SFRs all decrease systematically from the sBzK population to the LBGs, and to the BX/BMs. The combined color selections yield a total cosmic SFR density of 0.18 $\pm$ 0.03 M_sun yr^{-1} Mpc^{-3} for K_AB <= 24. We find that 65% of the star formation is in galaxies with E(B-V) > 0.25 mag, even though they are only one-fourth of the census by number.

preprint2011arXiv

A Local Baseline of the Black Hole Mass Scaling Relations for Active Galaxies. I. Methodology and Results of Pilot Study

We present high-quality Keck/LRIS longslit spectroscopy of a pilot sample of 25 local active galaxies selected from the SDSS (0.02<z<0.1; MBH>10^7 M_sun) to study the relations between black hole mass (MBH) and host-galaxy properties. We determine stellar kinematics of the host galaxy, deriving stellar-velocity dispersion profiles and rotation curves from three spectral regions (including CaH&K, MgIb triplet, and CaII triplet). In addition, we perform surface photometry on SDSS images, using a newly developed code for joint multi-band analysis. BH masses are estimated from the width of the Hbeta emission line and the host-galaxy free 5100A AGN luminosity. Combining results from spectroscopy and imaging allows us to study four MBH scaling relations: MBH-sigma, MBH-L(sph), MBH-M(sph,*), MBH-M(sph,dyn). We find the following results. First, stellar-velocity dispersions determined from aperture spectra (e.g. SDSS fiber spectra or unresolved data from distant galaxies) can be biased, depending on aperture size, AGN contamination, and host-galaxy morphology. However, such a bias cannot explain the offset seen in the MBH-sigma relation at higher redshifts. Second, while the CaT region is the cleanest to determine stellar-velocity dispersions, both the MgIb region, corrected for FeII emission, and the CaHK region, although often swamped by the AGN powerlaw continuum and emission lines, can give results accurate to within a few percent. Third, the MBH scaling relations of our pilot sample agree in slope and scatter with those of other local active and inactive galaxies. In the next papers of the series we will quantify the scaling relations, exploiting the full sample of ~100 objects.

preprint2011arXiv

AKARI Observation of the North Ecliptic Pole (NEP) Supercluster at z = 0.087: mid-infrared view of transition galaxies

We present the mid-infrared (MIR) properties of galaxies within a supercluster in the North Ecliptic Pole region at z?0.087 observed with the AKARI satellite. We use data from the AKARI NEP-Wide (5.4 deg2) IR survey and the CLusters of galaxies EVoLution studies (CLEVL) mission program. We show that near-IR (3 μm)-mid- IR (11 μm) color can be used as an indicator of the specific star formation rate and the presence of intermediate age stellar populations. From the MIR observations, we find that red-sequence galaxies consist not only of passively evolving red early-type galaxies, but also of 1) "weak-SFG" (disk-dominated star-forming galaxies which have star formation rates lower by \sim 4 \times than blue-cloud galaxies), and 2) "intermediate- MXG" (bulge-dominated galaxies showing stronger MIR dust emission than normal red early-type galaxies). Those two populations can be a set of transition galaxies from blue, star-forming, late-type galaxies evolving into red, quiescent, early-type ones. We find that the weak-SFG are predominant at intermediate masses (1010M\odot < M\star < 1010.5M\odot) and are typically found in local densities similar to the outskirts of galaxy clusters. As much as 40% of the supercluster member galaxies in this mass range can be classified as weak-SFGs, but their proportion decreases to < 10% at larger masses (M\star > 1010.5 M\odot) at any galaxy density. The fraction of the intermediate-MXG among red- sequence galaxies at 1010M\odot < M\star < 1011M\odot also decreases as the density and mass increase. In particular, \sim42% of the red-sequence galaxies with early-type morphologies are classified as intermediate-MXG at intermediate densities. These results suggest that the star formation activity is strongly dependent on the stellar mass, but that the morphological transformation is mainly controlled by the environment.

preprint2011arXiv

Completing the Census of Ly-alpha Emitters at the Reionization Epoch

We carried out extended spectroscopic confirmations of Ly-alpha emitters (LAEs) at z=6.5 and 5.7 in the Subaru Deep Field. Now, the total number of spectroscopically confirmed LAEs is 45 and 54 at z=6.5 and 5.7, respectively, and at least 81% (70%) of our photometric candidates at z=6.5 (5.7) have been spectroscopically identified as real LAEs. We made careful measurements of the Ly-alpha luminosity, both photometrically and spectroscopically, to accurately determine the Ly-alpha and rest-UV luminosity functions (LFs). The substantially improved evaluation of the Ly-alpha LF at z=6.5 shows an apparent deficit from z=5.7 at least at the bright end, and a possible decline even at the faint end, though small uncertainties remain. The rest-UV LFs at z=6.5 and 5.7 are in good agreement, at least at the bright end, in clear contrast to the differences seen in the Ly-alpha LF. These results imply an increase in the neutral fraction of the intergalactic medium from z=5.7 to 6.5. The rest-frame equivalent width (EW_0) distribution at z=6.5 seems to be systematically smaller than z=5.7, and it shows an extended tail toward larger EW_0. The bright end of the rest-UV LF can be reproduced from the observed Ly-alpha LF and a reasonable EW_0-UV luminosity relation. Integrating this rest-UV LF provides the first measurement of the contribution of LAEs to the photon budget required for reionization. The derived UV LF suggests that the fractional contribution of LAEs to the photon budget among Lyman break galaxies significantly increases towards faint magnitudes. Low-luminosity LAEs could dominate the ionizing photon budget, though this inference depends strongly on the uncertain faint-end slope of the Ly-alpha LF.

preprint2011arXiv

The Lick AGN Monitoring Project: Recalibrating Single-Epoch Virial Black Hole Mass Estimates

We investigate the calibration and uncertainties of black hole mass estimates based on the single-epoch (SE) method, using homogeneous and high-quality multi-epoch spectra obtained by the Lick Active Galactic Nucleus (AGN) Monitoring Project for 9 local Seyfert 1 galaxies with black hole masses < 10^8 M_sun. By decomposing the spectra into their AGN and stellar components, we study the variability of the single-epoch Hbeta line width (full width at half-maximum intensity, FWHM_Hbeta; or dispersion, sigma_Hbeta) and of the AGN continuum luminosity at 5100A (L_5100). From the distribution of the "virial products" (~ FWHM_Hbeta^2 L_5100^0.5 or sigma_Hbeta^2 L_5100^0.5) measured from SE spectra, we estimate the uncertainty due to the combined variability as ~ 0.05 dex (12%). This is subdominant with respect to the total uncertainty in SE mass estimates, which is dominated by uncertainties in the size-luminosity relation and virial coefficient, and is estimated to be ~ 0.46 dex (factor of ~ 3). By comparing the Hbeta line profile of the SE, mean, and root-mean-square (rms) spectra, we find that the Hbeta line is broader in the mean (and SE) spectra than in the rms spectra by ~ 0.1 dex (25%) for our sample with FWHM_Hbeta < 3000 km/s. This result is at variance with larger mass black holes where the difference is typically found to be much less than 0.1 dex. To correct for this systematic difference of the Hbeta line profile, we introduce a line-width dependent virial factor, resulting in a recalibration of SE black hole mass estimators for low-mass AGNs.

preprint2011arXiv

The Mass of the Black Hole in Arp 151 from Bayesian Modeling of Reverberation Mapping Data

Supermassive black holes are believed to be ubiquitous at the centers of galaxies. Measuring their masses is extremely challenging yet essential for understanding their role in the formation and evolution of cosmic structure. We present a direct measurement of the mass of a black hole in an active galactic nucleus (Arp 151) based on the motion of the gas responsible for the broad emission lines. By analyzing and modeling spectroscopic and photometric time series, we find that the gas is well described by a disk or torus with an average radius of 3.99 +- 1.25 light days and an opening angle of 68.9 (+21.4, -17.2) degrees, viewed at an inclination angle of 67.8 +- 7.8 degrees (that is, closer to face-on than edge-on). The black hole mass is inferred to be 10^(6.51 +- 0.28) solar masses. The method is fully general and can be used to determine the masses of black holes at arbitrary distances, enabling studies of their evolution over cosmic time.

preprint2011arXiv

The Relation between Black Hole Mass and Host Spheroid Stellar Mass out to z~2

We combine Hubble Space Telescope images from the Great Observatories Origins Deep Survey with archival Very Large Telescope and Keck spectra of a sample of 11 X-ray selected broad-line active galactic nuclei in the redshift range 1<z<2 to study the black hole mass - stellar mass relation out to a lookback time of 10 Gyrs. Stellar masses of the spheroidal component are derived from multi-filter surface photometry. Black hole masses are estimated from the width of the broad MgII emission line and the 3000A nuclear luminosity. Comparing with a uniformly measured local sample and taking into account selection effects, we find evolution in the form M_BH/M_spheroid ~ (1+z)^(1.96+/-0.55), in agreement with our earlier studies based on spheroid luminosity. However, this result is more accurate because it does not require a correction for luminosity evolution and therefore avoids the related and dominant systematic uncertainty. We also measure total stellar masses. Combining our sample with data from the literature, we find M_BH/M_host ~ (1+z)^(1.15+/-0.15), consistent with the hypothesis that black holes (in the range M_BH ~ 10^8-9 M_sun) predate the formation of their host galaxies. Roughly one third of our objects reside in spiral galaxies; none of the host galaxies reveal signs of interaction or major merger activity. Combined with the slower evolution in host stellar masses compared to spheroid stellar masses, our results indicate that secular evolution or minor mergers play a non-negligible role in growing both BHs and spheroids.

preprint2010arXiv

Overdensity of i'-Dropout Galaxies in the Subaru Deep Field: A Candidate Protocluster at z ~ 6

We investigate the sky distribution of z ~ 6 Lyman break galaxies selected as i'-dropouts having i' - z' > 1.45 down to z' < 26.5 in the Subaru Deep Field (SDF). We discover 37 i'-dropouts clustered in a projected comoving 21.6 x 21.6 Mpc^2 region at z = 6, showing a local density excess. Carrying out follow-up spectroscopy, we identify four of them as Lyman-alpha emitters at z = 5.92, 6.01, 6.03 and 6.03 (spread over a distance of 46.6 Mpc). The number density of the cluster itself in SDF is ~ 2.2 x 10^{-7} Mpc^{-3}, smaller than those of protoclusters (i.e., forming galaxy clusters) at z ~ 2-5.7. Also, the structure shows ~4-21 times larger galaxy number density than those of z ~ 6 galaxies in a general field. It has a mass of M ~ 1.5^{+1.8}_{-0.5} x 10^{15}M_sun, comparable to those of z ~ 0-5 protoclusters. Since the contamination of our sample by interlopers is estimated to be quite low, 5.9%, most of the other unconfirmed i'-dropouts in the overdense region can be also z ~ 6 galaxies. Hence, it could be a candidate forming cluster at z ~ 6, representing a progenitor of galaxy clusters seen in the recent-day Universe.

preprint2010arXiv

Spitzer Space Telescope Constraint on the Stellar Mass of a z = 6.96 Lyman Alpha Emitter

We obtained mid-infrared 3.6 and 4.5 micron imaging of a z=6.96 Lyman alpha emitter (LAE) IOK-1 discovered in the Subaru Deep Field, using Spitzer Space Telescope Infrared Array Camera observations. After removal of a nearby bright source, we find that IOK-1 is not significantly detected in any of these infrared bands to m_3.6 ~ 24.00 and m_4.5 ~ 23.54 at 3 sigma. Fitting population synthesis models to the spectral energy distribution consisting of the upper limit fluxes of the optical to infrared non-detection images and fluxes in detection images, we constrain the stellar mass M* of IOK-1. This LAE could have either a mass as low as M* <~ 2-9 x 10^8 Msun for the young age (<~ 10 Myr) and the low dust reddening (A_V ~ 0) or a mass as large as M* <~ 1-4 x 10^{10} Msun for either the old age (> 100 Myr) or the high dust reddening (A_V ~ 1.5). This would be within the range of masses of z ~ 3-6.6 LAEs studied to date, ~ 10^6-10^{10} Msun. Hence, IOK-1 is not a particularly unique galaxy with extremely high mass or low mass but is similar to one of the LAEs seen at the later epochs.

preprint2010arXiv

The Lick AGN Monitoring Project: Alternate Routes to a Broad-line Region Radius

It is now possible to estimate black hole masses across cosmic time, using broad emission lines in active galaxies. This technique informs our views of how galaxies and their central black holes coevolve. Unfortunately, there are many outstanding uncertainties associated with these "virial" mass estimates. One of these comes from using the accretion luminosity to infer a size for the broad-line region. Incorporating the new sample of low-luminosity active galaxies from our recent monitoring campaign at Lick Observatory, we recalibrate the radius-luminosity relation with tracers of the accretion luminosity other than the optical continuum. We find that the radius of the broad-line region scales as the square root of the X-ray and Hbeta luminosities, in agreement with recent optical studies. On the other hand, the scaling appears to be marginally steeper with narrow-line luminosities. This is consistent with a previously observed decrease in the ratio of narrow-line to X-ray luminosity with increasing total luminosity. The radius of the broad-line region correlates most tightly with Hbeta luminosity, while the X-ray and narrow-line relations both have comparable scatter of a factor of two. These correlations provide useful alternative virial BH masses in objects with no detectable optical/UV continuum emission, such as high-redshift galaxies with broad emission lines, radio-loud objects, or local active galaxies with galaxy-dominated continua.

preprint2010arXiv

The Lick AGN Monitoring Project: Reverberation Mapping of Optical Hydrogen and Helium Recombination Lines

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission, which we have previously reported. We present here the light curves for the Halpha, Hgamma, HeII 4686, and HeI 5876 emission lines and the time lags for the emission-line responses relative to changes in the continuum flux. Combining each emission-line time lag with the measured width of the line in the variable part of the spectrum, we determine a virial mass of the central supermassive black hole from several independent emission lines. We find that the masses are generally consistent within the uncertainties. The time-lag response as a function of velocity across the Balmer line profiles is examined for six of the AGNs. Finally we compare several trends seen in the dataset against the predictions from photoionization calculations as presented by Korista & Goad. We confirm several of their predictions, including an increase in responsivity and a decrease in the mean time lag as the excitation and ionization level for the species increases. Further confirmation of photoionization predictions for broad-line gas behavior will require additional monitoring programs for these AGNs while they are in different luminosity states. [abridged]

preprint2010arXiv

The Lick AGN Monitoring Project: The MBH - sigma Relation For Reverberation-Mapped Active Galaxies

(Abridged) To investigate the black hole mass (MBH) vs. stellar velocity dispersion relation of active galaxies, we measured the velocity dispersions of a sample of local Seyfert 1 galaxies, for which we have recently determined MBH using reverberation mapping. For most objects, velocity dispersions were measured from high S/N ratio optical spectra centered on the Ca II triplet region (~8500 A), obtained at the Keck, Palomar, and Lick Observatories. For two objects, in which the Ca II triplet region was contaminated by nuclear emission, the measurement was based on high-quality H-band spectra obtained with the OSIRIS at the Keck-II Telescope. Combining our new measurements with data from the literature, we assemble a sample of 24 active galaxies with stellar velocity dispersions and reverberation-based MBH measurements in the range of black hole mass 10^6<MBH/M_sun<10^9. We use this sample to obtain reverberation mapping constraints on the slope and intrinsic scatter of the M-sigma relation of active galaxies. Assuming a constant virial coefficient f for the reverberation mapping black hole masses, we find a slope b=3.55+-0.60 and the intrinsic scatter sigma_int=0.43+-0.08 dex in the relation log (MBH/M_sun) = a + b log (sigma/200 km/s), which are consistent with those found for quiescent galaxies. We derive an updated value of the virial coefficient f by finding the value which places the reverberation masses in best agreement with the M-sigma relation of quiescent galaxies; using the quiescent M-sigma relation determined by Gultekin et al. we find log f=0.72^{+0.09-0.10} with an intrinsic scatter of 0.44+-0.07 dex. The uncertainty of the virial coefficient remains one of the main sources of the uncertainty in black hole mass determinations using reverberation mapping, and therefore also in single-epoch spectroscopic estimates of MBH in active galaxies.

preprint2009arXiv

Cosmic Evolution of Black Holes and Spheroids. IV. The BH Mass - Spheroid Luminosity Relation

From high-resolution images of 23 Seyfert-1 galaxies at z=0.36 and z=0.57 obtained with the Near Infrared Camera and Multi-Object Spectrometer on board the Hubble Space Telescope (HST), we determine host-galaxy morphology, nuclear luminosity, total host-galaxy luminosity and spheroid luminosity. Keck spectroscopy is used to estimate black hole mass (M_BH). We study the cosmic evolution of the M_BH-spheroid luminosity (L_sph) relation. In combination with our previous work, totaling 40 Seyfert-1 galaxies, the covered range in BH mass is substantially increased, allowing us to determine for the first time intrinsic scatter and correct evolutionary trends for selection effects. We re-analyze archival HST images of 19 local reverberation-mapped active galaxies to match the procedure adopted at intermediate redshift. Correcting spheroid luminosity for passive luminosity evolution and taking into account selection effects, we determine that at fixed present-day V-band spheroid luminosity, M_BH/L_sph \propto (1+z)^(2.8+/-1.2). When including a sample of 44 quasars out to z=4.5 taken from the literature, with luminosity and BH mass corrected to a self-consistent calibration, we extend the BH mass range to over two orders of magnitude, resulting in M_BH/L_sph \propto (1+z)^(1.4+/-0.2). The intrinsic scatter of the relation, assumed constant with redshift, is 0.3+/-0.1 dex (<0.6 dex at 95% CL). The evolutionary trend suggests that BH growth precedes spheroid assembly. Interestingly, the M_BH-total host-galaxy luminosity relation is apparently non-evolving. It hints at either a more fundamental relation or that the spheroid grows by a redistribution of stars. However, the high-z sample does not follow this relation, indicating that major mergers may play the dominant role in growing spheroids above z~1.

preprint2009arXiv

Lyman Break Galaxies at z~1.8-2.8: GALEX/NUV Imaging of the Subaru Deep Field

Abridged: A photometric sample of ~7100 V<25.3 Lyman break galaxies (LBGs) has been selected by combining Subaru/Suprime-Cam BVRci'z' data with deep GALEX/NUV imaging of the Subaru Deep Field. Follow-up spectroscopy confirmed 24 LBGs at 1.5<z<2.7. Among the optical spectra, 12 have Ly-alpha emission with rest-frame equivalent widths of ~5-60AA. The success rate for identifying LBGs as NUV-dropouts at 1.5<z<2.7 is 86%. The rest-frame UV (1700AA) luminosity function (LF) is constructed from the photometric sample with corrections for stellar contamination and z<1.5 interlopers. The LF is 1.7+/-0.1 times higher than those of z~2 BXs and z~3 LBGs. Three explanations were considered, and it is argued that significantly underestimating low-z contamination or effective comoving volume is unlikely: the former would be inconsistent with the spectroscopic sample at 93% confidence, and the second explanation would not resolve the discrepancy. The third scenario is that different photometric selection of the samples yields non-identical galaxy populations, such that some BX galaxies are LBGs and vice versa. This argument is supported by a higher surface density of LBGs at all magnitudes while the redshift distribution of the two populations is nearly identical. This study, when combined with other star-formation rate (SFR) density UV measurements from LBG surveys, indicates that there is a rise in the SFR density: a factor of 3-6 (3-10) increase from z~5 (z~6) to z~2, followed by a decrease to z~0. This result, along with past sub-mm studies that find a peak at z~2 in their redshift distribution, suggest that z~2 is the epoch of peak star-formation. Additional spectroscopy is required to characterize the complete shape of the z~2 LBG UV LF via measurements of contamination and accurate distances.

preprint2009arXiv

Spitzer-IRS high resolution spectroscopic survey of the 12 micron Seyfert galaxies: II. Results for the Complete Dataset

We present the Spitzer 10-37um IRS high resolution (R~600) spectroscopic survey of the Seyfert galaxies of the 12 micron Galaxy Sample. The new spectra of 61 galaxies, with those already published, gives us a total of 91 12micron Seyfert galaxies observed, out of 112. We use an improved AGN classification for Seyfert galaxies: instead of the type 1 and 2 classes, we use the spectropolarimetric data to divide them into "AGN 1" and "AGN 2", where AGN 1's include all broad-line objects, including the Seyfert 2's showing hidden broad lines in polarized light, while AGN 2's contains only Seyferts with no detectable broad lines at all. We present various mid-IR observables and we find that these properties characterize the AGN 1's objects as a single family, with strongly AGN-dominated spectra. In contrast, the AGN 2's can be divided in two groups, the first one with properties similar to the AGN 1's and the second similar to the non-Seyfert galaxies, such as LINERs or starburst galaxies. We computed a semianalytical model to estimate the AGN and the starburst contributions to the mid-IR galaxy emission at 19um. We find that AGN 1 have an AGN contribution >73% and AGN 2 >45% of their total emission at 19um. The detection of [NeV] lines is an almost perfect signature of energy production by an AGN. We present mean spectra of the various AGN categories. We derive the first local luminosity functions for the mid-infrared lines and the PAH feature. No statistical difference is found in the space densities for Seyfert 1's and 2's, nor for the new classes of AGN 1's and 2's. The global output of accretion-powered galactic nuclei in the local universe is derived from the correlation between [NeV] line and the nonstellar IR continuum luminosity.

preprint2009arXiv

The Lick AGN Monitoring Project: Broad-Line Region Radii and Black Hole Masses from Reverberation Mapping of Hbeta

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7 M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission. We present here the light curves for the objects in this sample and the subsequent Hbeta time lags for the nine objects where these measurements were possible. The Hbeta lag time is directly related to the size of the broad-line region, and by combining the lag time with the measured width of the Hbeta emission line in the variable part of the spectrum, we determine the virial mass of the central supermassive black hole in these nine AGNs. The absolute calibration of the black hole masses is based on the normalization derived by Onken et al. We also examine the time lag response as a function of velocity across the Hbeta line profile for six of the AGNs. The analysis of four leads to ambiguous results with relatively flat time lags as a function of velocity. However, SBS 1116+583A exhibits a symmetric time lag response around the line center reminiscent of simple models for circularly orbiting broad-line region (BLR) clouds, and Arp 151 shows an asymmetric profile that is most easily explained by a simple gravitational infall model. Further investigation will be necessary to fully understand the constraints placed on physical models of the BLR by the velocity-resolved response in these objects.

preprint2009arXiv

The Lick AGN Monitoring Project: Photometric Light Curves and Optical Variability Characteristics

The Lick AGN Monitoring Project targeted 13 nearby Seyfert 1 galaxies with the intent of measuring the masses of their central black holes using reverberation mapping. The sample includes 12 galaxies selected to have black holes with masses roughly in the range 10^6-10^7 solar masses, as well as the well-studied AGN NGC 5548. In conjunction with a spectroscopic monitoring campaign, we obtained broad-band B and V images on most nights from 2008 February through 2008 May. The imaging observations were carried out by four telescopes: the 0.76-m Katzman Automatic Imaging Telescope (KAIT), the 2-m Multicolor Active Galactic Nuclei Monitoring (MAGNUM) telescope, the Palomar 60-in (1.5-m) telescope, and the 0.80-m Tenagra II telescope. Having well-sampled light curves over the course of a few months is useful for obtaining the broad-line reverberation lag and black hole mass, and also allows us to examine the characteristics of the continuum variability. In this paper, we discuss the observational methods and the photometric measurements, and present the AGN continuum light curves. We measure various variability characteristics of each of the light curves. We do not detect any evidence for a time lag between the B- and V-band variations, and we do not find significant color variations for the AGNs in our sample.

preprint2008arXiv

Star Formation Rates and Metallicities of K-selected Star Forming Galaxies at z~2

We present spectroscopy of 15 star-forming BzK galaxies (sBzKs) with K(AB)<23 in the Subaru Deep Field, for which Halpha and some other emission lines are detected in 0.9 to 2.3 micron spectra with a resolution of R=500. Using Halpha luminosities, we obtain star formation rates (SFRs), and then specific SFRs (SSFRs) dividing SFRs by stellar masses, which are derived from SED fitting to BVRi'z'K photometry. It is found that sBzKs with higher stellar masses have larger SFRs. A negative correlation is seen between stellar mass and SSFR, which is consistent with the previous results for z~2 galaxies. This implies that a larger growth of stellar mass occurs in less massive galaxies. In addition, gas-phase oxygen abundances, 12+log(O/H), are derived from the ratio of NII(lambda 6584) to Halpha using the N2 index method. We have found a correlation between stellar mass and oxygen abundance in the sense that more massive sBzKs tend to be more metal rich, which is qualitatively consistent with the relation for UV-selected z~2 galaxies. However, the metallicity of the sBzKs is ~0.2 dex higher than that of UV-selected galaxies with similar stellar masses, which is significant considering the small uncertainties. The sBzKs in our sample have redder R-K colors than the UV-selected galaxies. This galaxy color-dependence in the oxygen abundance may be caused by older or dustier galaxies having higher metallicities at z~2.

preprint2007arXiv

A search for Lyman Break Galaxies at z>8 in the NICMOS Parallel Imaging Survey

We have selected 14 J-dropout Lyman Break Galaxy (LBG) candidates with J110 - H160 > 2.5 from the NICMOS Parallel Imaging Survey. This survey consists of 135 square arcminutes of imaging in 228 independent sight lines, reaching average 5 sigma sensitivities of J110 = 25.8 and H160 = 25.6 (AB). Distinguishing these candidates from dust reddened star forming galaxies at z ~ 2-3 is difficult, and will require longer wavelength observations. We consider the likelihood that any J-dropout LBGs exist in this survey, and find that if L*(z=9.5) is significantly brighter than L*(z=6) (a factor of four), then a few J-dropout LBGs are likely. A similar increase in luminosity has been suggested by Eyles et al. and Yan et al., but the magnitude of this increase is uncertain.

preprint2007arXiv

Comparing and calibrating black hole mass estimators for distant active galactic nuclei

Black hole mass is a fundamental property of active galactic nuclei (AGNs). In the distant universe, black hole mass is commonly estimated using the MgII, Hbeta, or Halpha emission line widths and the optical/UV continuum or line luminosities, as proxies for the characteristic velocity and size of the broad-line region. Although they all have a common calibration in the local universe, a number of different recipes are currently used in the literature. It is important to verify the relative accuracy and consistency of the recipes, as systematic changes could mimic evolutionary trends when comparing various samples. At z=0.36, all three lines can be observed at optical wavelengths, providing a unique opportunity to compare different empirical recipes. We use spectra from the Keck Telescope and the Sloan Digital Sky Survey to compare black hole mass estimators for a sample of nineteen AGNs at this redshift. We compare popular recipes available from the literature, finding that mass estimates can differ up to 0.38+-0.05 dex in the mean (or 0.13+-0.05 dex, if the same virial coefficient is adopted). Finally, we provide a set of 30 internally self consistent recipes for determining black hole mass from a variety of observables. The intrinsic scatter between cross-calibrated recipes is in the range 0.1-0.3 dex. This should be considered as a lower limit to the uncertainty of the black hole mass estimators.

preprint2001arXiv

The Distance to SN 1999em in NGC 1637 from the Expanding Photosphere Method

We present 30 optical spectra and 49 photometric epochs sampling the first 517 days after discovery of supernova (SN) 1999em, and derive its distance through the expanding photosphere method (EPM). SN 1999em is shown to be a Type II-plateau (II-P) event, with a photometric plateau lasting until about 100 days after explosion. We identify the dominant ions responsible for most of the absorption features seen in the optical portion of the spectrum during the plateau phase. We find the distance to SN 1999em to be D = 8.2 +/- 0.6 Mpc, with an explosion date of 5.3 +/- 1.4 days before discovery. We examine potential sources of systematic error in EPM-derived distances, and find the most significant to result from uncertainty in the theoretical modeling of the flux distribution emitted by the SN photosphere (i.e., the ``flux dilution factor''). We compare previously derived EPM distances to 5 SNe II in galaxies for which a recently revised Cepheid distance exists from the HST Key Project and find D(Cepheids) / D(EPM) = 0.96 +/- 0.09. Finally, we investigate the possible use of SNe II-P as standard candles and find that for 8 photometrically confirmed SNe II-P with previously derived EPM distances and SN 1999em, the mean plateau absolute brightness is M_V(plateau) = -16.4 +/- 0.6 mag, implying that distances good to ~30% (1-sigma) may be possible without the need for a complete EPM analysis. At M_V(plateau) = -15.9 +/- 0.2 mag, SN 1999em is somewhat fainter than the average SN II-P. The general consistency of absolute SNe II-P brightness during the plateau suggests that the standard candle assumption may allow SNe II-P to be viable cosmological beacons at z > 2.

preprint1994arXiv

The jet-disk symbiosis. II. Interpreting the Radio/UV correlations in quasars

We investigate the correlation between the accretion disk (UV) luminosity and the radio core emission of a quasar sample. In a radio/$L_{\rm disk}$ plot we find the quasars to be separated into four classes: core dominated quasars (CDQ), lobe dominated quasars (LDQ), radio-intermediate quasars (RIQ) and radio weak quasars. In general the radio core emission scales with the disk luminosity, especially in the radio weak quasars. This shows that radio and UV emission have a common energy source and that the difference between radio loud and radio weak is established already on the parsec scale. We investigate the possibility that radio jets are responsible for the radio core emission in radio loud and radio weak quasars. Comparing our data with a simple jet emission model that takes the limits imposed by energy and mass conservation in a coupled jet-disk system into account, we find that radio loud jets carry a total power $Q_{\rm jet}$ that is at least 1/3 of the observed disk luminosity $L_{\rm disk}$. For the electron population one is forced to postulate an efficient process producing a large number of pairs and/or injecting electrons with a distribution with low-energy cut-off around 50 MeV -- secondary pair production in hadronic cascades could be such a process. The bulk Lorentz factor of the jet is limited to a narrow region ($3\laγ_{\rm j}\la10$)). The radio emission of radio weak quasars can be explained with exactly the same parameters for a powerful relativistic jet if secondary pair production, as suggested for radio loud jets, is inhibited. There is evidence that RIQ are the relativistically boosted population of radio weak quasar jets. We can estimate the hidden disk luminosity of FR II galaxies and find that this is consistent with FRII being misdirected quasars.