Source author record

Christopher S. Kochanek

Christopher S. Kochanek 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

92works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

92 published item(s)

preprint2022arXiv

ASAS-SN follow-up of IceCube high-energy neutrino alerts

We report on the search for optical counterparts to IceCube neutrino alerts released between April 2016 and August 2021 with the All-Sky Automated Survey for SuperNovae (ASAS-SN). Despite the discovery of a diffuse astrophysical high-energy neutrino flux in 2013, the source of those neutrinos remains largely unknown. Since 2016, IceCube has published likely-astrophysical neutrinos as public realtime alerts. Through a combination of normal survey and triggered target-of-opportunity observations, ASAS-SN obtained images within 1 hour of the neutrino detection for 20% (11) of all observable IceCube alerts and within one day for another 57% (32). For all observable alerts, we obtained images within at least two weeks from the neutrino alert. ASAS-SN provides the only optical follow-up for about 17% of IceCube's neutrino alerts. We recover the two previously claimed counterparts to neutrino alerts, the flaring-blazar TXS 0506+056 and the tidal disruption event AT2019dsg. We investigate the light curves of previously-detected transients in the alert footprints, but do not identify any further candidate neutrino sources. We also analysed the optical light curves of Fermi 4FGL sources coincident with high-energy neutrino alerts, but do not identify any contemporaneous flaring activity. Finally, we derive constraints on the luminosity functions of neutrino sources for a range of assumed evolution models.

preprint2022arXiv

Investigating the Nature of the Luminous Ambiguous Nuclear Transient ASASSN-17jz

We present observations of the extremely luminous but ambiguous nuclear transient (ANT) ASASSN-17jz, spanning roughly 1200 days of the object's evolution. ASASSN-17jz was discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) in the galaxy SDSS J171955.84+414049.4 on UT 2017 July 27 at a redshift of $z=0.1641$. The transient peaked at an absolute $B$-band magnitude of $M_{B,{\rm peak}}=-22.81$, corresponding to a bolometric luminosity of $L_{\rm bol,peak}=8.3\times10^{44}$~erg~s$^{-1}$, and exhibited late-time ultraviolet emission that was still ongoing in our latest observations. Integrating the full light curve gives a total emitted energy of $E_{\rm tot}=(1.36\pm0.08)\times10^{52}$~erg, with $(0.80\pm0.02)\times10^{52}$~erg of this emitted within 200 days of peak light. This late-time ultraviolet emission is accompanied by increasing X-ray emission that becomes softer as it brightens. ASASSN-17jz exhibited a large number of spectral emission lines most commonly seen in active galactic nuclei (AGNs) with little evidence of evolution. It also showed transient Balmer features which became fainter and broader over time, and are still being detected $>1000$ days after peak brightness. We consider various physical scenarios for the origin of the transient, including supernovae (SNe), tidal disruption events (TDEs), AGN outbursts, and ANTs. We find that the most likely explanation is that ASASSN-17jz was an SN~IIn occurring in or near the disk of an existing AGN, and that the late-time emission is caused by the AGN transitioning to a more active state.

preprint2022arXiv

Unveiling the Nature of SN 2011fh: a Young and Massive Star Gives Rise to a Luminous SN 2009ip-like Event

In recent years, many Type IIn supernovae have been found to share striking similarities with the peculiar SN 2009ip, whose true nature is still under debate. Here, we present 10 years of observations of SN 2011fh, an interacting transient with spectroscopic and photometric similarities to SN 2009ip. SN 2011fh had a M$_r \sim -16$ mag brightening event, followed by a brighter M$_r \sim -18$ mag luminous outburst in August 2011. The spectra of SN 2011fh are dominated by narrow to intermediate Balmer emission lines throughout its evolution, with P Cygni profiles indicating fast-moving material at $\sim 6400 \ \textrm{km s}^{-1}$. HST/WFC3 observations from October of 2016 revealed a bright source with M$_{F814W} \approx -13.3$ mag, indicating that we are seeing the ongoing interaction of the ejecta with the circumstellar material or that the star might be going through an eruptive phase five years after the luminous outburst of 2011. Using HST photometry of the stellar cluster around SN 2011fh, we estimated an age of $\sim 4.5$ Myr for the progenitor, which implies a stellar mass of $\sim 60$ M$_\odot$, using single-star evolution models, or a mass range of $35 - 80 \ \textrm{M}_\odot$, considering a binary system. We also show that the progenitor of SN 2011fh exceeded the classical Eddington limit by a large factor in the months preceding the luminous outburst of 2011, suggesting strong super-Eddington winds as a possible mechanism for the observed mass-loss. These findings favor an energetic outburst in a young and massive star, possibly a luminous blue variable.

preprint2022arXiv

Variability Selected Active Galactic Nuclei from ASAS-SN Survey: Constraining the Low Luminosity AGN Population

Low luminosity active galactic nuclei (LLAGN) probe accretion physics in the low Eddington regime and can provide additional clues about galaxy evolution. AGN variability is ubiquitous and thus provides a reliable tool for finding AGN. We analyze the All-Sky Automated Survey for SuperNovae light curves of 1218 galaxies with $g<14$ mag and Sloan Digital Sky Survey spectra in search of AGN. We find 37 objects that are both variable and have AGN-like structure functions, which is about 3% of the sample. The majority of the variability selected AGN are LLAGN with Eddington ratios ranging from $10^{-4}$ to $10^{-2}$. We thus estimate the fraction of LLAGN in the population of galaxies as 2% down to a median Eddington ratio of $2\times10^{-3}$. Combining the BPT line ratio diagnostics and the broad-line AGN, up to $\sim$60% of the AGN candidates are confirmed spectroscopically. The BPT diagnostics also classified 10-30% of the candidates as star forming galaxies rather than AGN.

preprint2021arXiv

Optical Confirmation of X-ray selected Galaxy clusters from the Swift AGN and Cluster survey with MDM and Pan-STARRS Data (Paper III)

To understand structure formation in the universe and impose stronger constraints on the cluster mass function and cosmological models, it is important to have large galaxy cluster catalogs. The Swift AGN and Cluster Survey is a serendipitous X-ray survey aimed at building a large statistically selected X-ray cluster catalog with 442 cluster candidates in its first release. Our initial SDSS follow-up study confirmed $50\%$ of clusters in the SDSS footprint as z $<$ 0.5 clusters. Here, we present further optical follow-up analysis of 248 (out of 442) cluster candidates from the Swift cluster catalog using multi-band imaging from the MDM $2.4m$ telescope and the Pan-STARRS survey. We report the optical confirmation of 55 clusters with $> 3σ$ galaxy overdensities and detectable red sequences in the color-magnitude space. The majority of these confirmed clusters have redshifts z $<$ 0.6. The remaining candidates are potentially higher redshift clusters that are excellent targets for infrared observations. We report the X-ray luminosity and the optical richness for these confirmed clusters. We also discuss the distinction between X-ray and optical observables for the detected and non-detected cluster candidates.

preprint2021arXiv

The Rapid X-ray and UV Evolution of ASASSN-14ko

ASASSN-14ko is a recently discovered periodically flaring transient at the center of the AGN ESO 253-G003 with a slowly decreasing period. Here we show that the flares originate from the northern, brighter nucleus in this dual-AGN, post-merger system. The light curves for the two flares that occurred in May 2020 and September 2020 are nearly identical over all wavelengths. For both events, Swift observations showed that the UV and optical wavelengths brightened in unison. The effective temperature of the UV/optical emission rises and falls with the increase and subsequent decline in the luminosity. The X-ray flux, in contrast, first rapidly drops over $\sim$2.6 days, rises for $\sim$5.8 days, drops again over $\sim$4.3 days and then recovers. The X-ray spectral evolution of the two flares differ, however. During the May 2020 peak the spectrum softened with increases in the X-ray luminosity, while we observed the reverse for the September 2020 peak.

preprint2020arXiv

Direct evidence for shock-powered optical emission in a nova

Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

preprint2020arXiv

Examining a Peak-Luminosity/Decline-Rate Relationship for Tidal Disruption Events

We compare the luminosity, radius, and temperature evolution of the UV/optical blackbodies for 21 well-observed tidal disruption events (TDEs), 8 of which were discovered by the All-Sky Automated Survey for Supernovae. We find that the blackbody radii generally increase prior to peak and slowly decline at late times. The blackbody temperature evolution is generally flat, with a few objects showing small-scale variations. The bolometric UV/optical luminosities generally evolve smoothly and flatten out at late times. Finally, we find an apparent correlation between the peak luminosity and the decline rate of TDEs. This relationship is strongest when comparing the peak luminosity to its decline over 40 days. A linear fit yields $\log_{10}( L_{peak}) = (44.1^{+0.1}_{-0.1}) + (1.6^{+0.4}_{-0.2})(ΔL_{40} + 0.5)$ in cgs, where $ΔL_{40} = \log_{10}(L_{40} / L_{peak})$.

preprint2020arXiv

Progenitor, Precursor and Evolution of the Dusty Remnant of the Stellar Merger M31-LRN-2015

M31-2015-LRN is a likely stellar merger discovered in the Andromeda Galaxy in 2015. We present new optical to mid-infrared photometry and optical spectroscopy for this event. Archival data shows that the source started to brighten $\sim$2 years before the nova event. During this precursor phase, the source brightened by $\sim$3 mag. The lightcurve at 6 and 1.5 months before the main outburst may show periodicity, with periods of 16$\pm$0.3 and 28.1$\pm$1.4 days respectively. This complex emission may be explained by runaway mass loss from the system after the binary undergoes Roche-lobe overflow, leading the system to coalesce in tens of orbital periods. While the progenitor spectral energy distribution shows no evidence of pre-existing warm dust in system, the remnant forms an optically thick dust shell at $\sim$4 months after the outburst peak. The optical depth of the shell increases dramatically after 1.5 years, suggesting the existence of shocks that enhance the dust formation process. We propose that the merger remnant is likely an inflated giant obscured by a cooling shell of gas with mass $\sim0.2$ M$_{\odot}$ ejected at the onset of the common envelope phase.

preprint2020arXiv

Response to Comment on "A Non-Interacting Low-Mass Black Hole -- Giant Star Binary System"

van den Heuvel & Tauris argue that if the red giant star in the system 2MASS J05215658+4359220 has a mass of 1 solar mass (M$_\odot$), then its unseen companion could be a binary composed of two 0.9 M$_\odot$ stars, making a triple system. We contend that the existing data are most consistent with a giant of mass $3.2^{+1.0}_{-1.0}$ M$_\odot$, implying a black hole companion of $3.3^{+2.8}_{-0.7}$ M$_\odot$.

preprint2020arXiv

The Rise and Fall of ASASSN-18pg: Following a TDE from Early To Late Times

We present nearly 500 days of observations of the tidal disruption event ASASSN-18pg, spanning from 54 days before peak light to 441 days after peak light. Our dataset includes X-ray, UV, and optical photometry, optical spectroscopy, radio observations, and the first published spectropolarimetric observations of a TDE. ASASSN-18pg was discovered on 2018 July 11 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $d=78.6$ Mpc, and with a peak UV magnitude of $m\simeq14$ it is both one of the nearest and brightest TDEs discovered to-date. The photometric data allow us to track both the rise to peak and the long-term evolution of the TDE. ASASSN-18pg peaked at a luminosity of $L\simeq2.2\times10^{44}$ erg s$^{-1}$, and its late-time evolution is shallower than a flux $\propto t^{-5/3}$ power-law model, similar to what has been seen in other TDEs. ASASSN-18pg exhibited Balmer lines and spectroscopic features consistent with Bowen fluorescence prior to peak which remained detectable for roughly 225 days after peak. Analysis of the two-component H$α$ profile indicates that, if they are the result of reprocessing of emission from the accretion disk, the different spectroscopic lines may be coming from regions between $\sim10$ and $\sim60$ light-days from the black hole. No X-ray emission is detected from the TDE and there is no evidence of a jet or strong outflow detected in the radio. Our spectropolarimetric observations give no strong evidence for significant asphericity in the emission region, with the emission region having an axis ratio of at least $\sim0.65$.

preprint2019arXiv

A Catalog of M-dwarf Flares with ASAS-SN

We analyzed the light curves of 1376 early-to-late, nearby M dwarfs to search for white-light flares using photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN). We identified 480 M dwarfs with at least one potential flare employing a simple statistical algorithm that searches for sudden increases in $V$-band flux. After more detailed evaluation, we identified 62 individual flares on 62 stars. The event amplitudes range from $0.12 <ΔV < 2.04$ mag. Using classical-flare models, we place lower limits on the flare energies and obtain $V$-band energies spanning $2.0\times10^{30} \lesssim E_{V} \lesssim 6.9\times10^{35}$ erg. The fraction of flaring stars increases with spectral type, and most flaring stars show moderate to strong H$α$ emission. Additionally, we find that 14 of the 62 flaring stars are rotational variables, and they have shorter rotation periods and stronger H$α$ emission than non-flaring rotational variable M dwarfs.

preprint2019arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. X: The Tenth Year (2017)

Continuing the project described by Kato et al. (2009, PASJ, 61, S395, arXiv/0905.1757), we collected times of superhump maxima for 102 SU UMa-type dwarf novae observed mainly during the 2017 season and characterized these objects. WZ Sge-type stars identified in this study are PT And, ASASSN-17ei, ASASSN-17el, ASASSN-17es, ASASSN-17fn, ASASSN-17fz, ASASSN-17hw, ASASSN-17kd, ASASSN-17la, PNV J20205397+2508145 and TCP J00332502-3518565. We obtained new mass ratios for 7 objects using growing superhumps (stage A). ASASSN-17gf is an EI Psc-type object below the period minimum. CRTS J080941.3+171528 and DDE 51 are objects in the period gap and both showed long-lasting phase of stage A superhumps. We also summarized the recent advances in understanding of SU UMa-type and WZ Sge-type dwarf novae.

preprint2019arXiv

The Case for Strangulation in Low-Mass Hosts: DDO 113

We investigate the case for environmental quenching of the Fornax-mass satellite DDO 113, which lies only 9 kpc in projection from its host, the Large-Magellanic-Cloud-mass galaxy NGC 4214. DDO 113 was quenched about 1 Gyr ago and is virtually gas-free, while analogs in the field are predominantly star-forming and gas-rich. We use deep imaging obtained with the Large Binocular Telescope to show that DDO 113 exhibits no evidence of tidal disruption to a surface brightness of $μ_V\sim29$ mag $\text{arcsec}^{-2}$, based on both unresolved emission and resolved stars. Mass-analogs of DDO 113 in Illustris-1 with similar hosts, small projected separations, and no significant tidal stripping first fell into their host halo 2--6 Gyr ago, showing that tidal features (or lack thereof) can be used to constrain infall times in systems where there are few other constraints on the orbit of the satellite. With the infall time setting the clock for environmental quenching mechanisms, we investigate the plausibility of several such mechanisms. We find that strangulation, the cessation of cold gas inflows, is likely the dominant quenching mechanism for DDO 113, requiring a time-averaged mass-loading factor of $η=6-11$ for star-formation-driven outflows that is consistent with theoretical and observational constraints. Motivated by recent numerical work, we connect DDO 113's strangulation to the presence of a cool circumgalactic medium (CGM) around NGC 4214. This discovery shows that the CGM of low-mass galaxies can affect their satellites significantly and motivates further work on understanding the baryon cycle in low-mass galaxies.

preprint2016arXiv

An Extreme Analogue of $ε$ Aurigae: An M-giant Eclipsed Every 69 Years by a Large Opaque Disk Surrounding a Small Hot Source

We present TYC 2505-672-1 as a newly discovered and remarkable eclipsing system comprising an M-type red giant that undergoes a ~3.45 year long, near-total eclipse (depth of ~4.5 mag) with a very long period of ~69.1 yr. TYC 2505-672-1 is now the longest-period eclipsing binary system yet discovered, more than twice as long as that of the currently longest-period system, $ε$ Aurigae. We show from analysis of the light curve including both our own data and historical data spanning more than 120 yr and from modeling of the spectral energy distribution, both before and during eclipse, that the red giant primary is orbited by a moderately hot source (T$_{eff}$~8000 K) that is itself surrounded by an extended, opaque circumstellar disk. From the measured ratio of luminosities, the radius of the hot companion must be in the range 0.1-0.5 Rsun (depending on the assumed radius of the red giant primary), which is an order of magnitude smaller than that for a main sequence A star and 1-2 orders of magnitude larger than that for a white dwarf. The companion is therefore most likely a "stripped red giant" subdwarf-B type star destined to become a He white dwarf. It is however somewhat cooler than most sdB stars, implying a very low mass for this "pre-He-WD" star. The opaque disk surrounding this hot source may be a remnant of the stripping of its former hydrogen envelope. However, it is puzzling how this object became stripped, given that it is at present so distant (orbital semi-major axis of ~24 AU) from the current red giant primary star. Extrapolating from our calculated ephemeris, the next eclipse should begin in early UT 2080 April and end in mid UT 2083 September (eclipse center UT 2081 December 24). This system is poised to become an exemplar of a very rare class of systems, even more extreme in several respects than the well studied archetype $ε$ Aurigae.

preprint2016arXiv

DM Ori: A Young Star Occulted by a Disturbance in its Protoplanetary Disk

In some planet formation theories, protoplanets grow gravitationally within a young star's protoplanetary disk, a signature of which may be a localized disturbance in the disk's radial and/or vertical structure. Using time-series photometric observations by the Kilodegree Extremely Little Telescope South (KELT-South) project and the All-Sky Automated Survey for SuperNovae (ASAS-SN), combined with archival observations, we present the discovery of two extended dimming events of the young star, DM Ori. This young system faded by $\sim$1.5 mag from 2000 March to 2002 August and then again in 2013 January until 2014 September (depth $\sim$1.7 mag). We constrain the duration of the 2000-2002 dimming to be $<$ 860 days, and the event in 2013-2014 to be $<$ 585 days, separated by $\sim$12.5 years. A model of the spectral energy distribution (SED) indicates a large infrared excess consistent with an extensive circumstellar disk. Using basic kinematic arguments, we propose that DM Ori is likely being periodically occulted by a feature (possibly a warp or perturbation) in its circumstellar disk. In this scenario, the occulting feature is located $>$6 AU from the host star, moving at $\sim$14.6 km/s, and is $\sim$4.9 AU in width. This localized structure may indicate a disturbance such as may be caused by a protoplanet early in its formation.

preprint2016arXiv

Quasar Variability in the Mid-Infrared

The Decadal IRAC Bootes Survey is a mid-IR variability survey of the ~9 sq. deg. of the NDWFS Bootes Field and extends the time baseline of its predecessor, the Spitzer Deep, Wide-Field Survey (SDWFS), from 4 to 10 years. The Spitzer Space Telescope visited the field five times between 2004 and 2014 at 3.6 and 4.5 microns. We provide the difference image analysis photometry for a half a million mostly extragalactic sources. In mid-IR color-color plane, sources with quasar colors constitute the largest variability class (75%), 16% of the variable objects have stellar colors and the remaining 9% have the colors of galaxies. Adding the fifth epoch doubles the number of variable active galactic nuclei (AGNs) for the same false positive rates as in SDWFS, or increases the number of sources by 20% while decreasing the false positive rates by factors of 2-3 for the same variability amplitude. We quantify the ensemble mid-IR variability of ~1500 spectroscopically confirmed AGNs using single power-law structure functions (SFs), which we find to be steeper (index $γ=0.45$) than in the optical ($γ=0.3$), leading to much lower amplitudes at short time-lags. This provides evidence for large emission regions, smoothing out any fast UV/optical variations, as the origin of infrared quasar variability. The mid-IR AGN SF slope $γ$ seems to be uncorrelated with both the luminosity and rest-frame wavelength, while the amplitude shows an anti-correlation with the luminosity and a correlation with the rest-frame wavelength.

preprint2016arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. VIII: The Eighth Year (2015-2016)

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 128 SU UMa-type dwarf novae observed mainly during the 2015-2016 season and characterized these objects. The data have improved the distribution of orbital periods, the relation between the orbital period and the variation of superhumps, the relation between period variations and the rebrightening type in WZ Sge-type objects. Coupled with new measurements of mass ratios using growing stages of superhumps, we now have a clearer and statistically greatly improved evolutionary path near the terminal stage of evolution of cataclysmic variables. Three objects (V452 Cas, KK Tel, ASASSN-15cl) appear to have slowly growing superhumps, which is proposed to reflect the slow growth of the 3:1 resonance near the stability border. ASASSN-15sl, ASASSN-15ux, SDSS J074859.55+312512.6 and CRTS J200331.3-284941 are newly identified eclipsing SU UMa-type (or WZ Sge-type) dwarf novae. ASASSN-15cy has a short (~0.050 d) superhump period and appears to belong to EI Psc-type objects with compact secondaries having an evolved core. ASASSN-15gn, ASASSN-15hn, ASASSN-15kh and ASASSN-16bu are candidate period bouncers with superhump periods longer than 0.06 d. We have newly obtained superhump periods for 79 objects and 13 orbital periods, including periods from early superhumps. In order that the future observations will be more astrophysically beneficial and rewarding to observers, we propose guidelines how to organize observations of various superoutbursts.

preprint2016arXiv

The BOSS Emission-Line Lens Survey. III. : Strong Lensing of Ly$α$ Emitters by Individual Galaxies

We introduce the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$α$ EmitteR sYstems (BELLS GALLERY) Survey, which is a Hubble Space Telescope program to image a sample of galaxy-scale strong gravitational lens candidate systems with high-redshift Ly$α$ emitters (LAEs) as the background sources. The goal of the BELLS GALLERY Survey is to illuminate dark substructures in galaxy-scale halos by exploiting the small-scale clumpiness of rest-frame far-UV emission in lensed LAEs, and to thereby constrain the slope and normalization of the substructure-mass function. In this paper, we describe in detail the spectroscopic strong-lens selection technique, which is based on methods adopted in the previous Sloan Lens ACS (SLACS) Survey, BELLS, and SLACS for the Masses Survey. We present the BELLS GALLERY sample of the 21 highest-quality galaxy--LAE candidates selected from $\approx 1.4 \times 10^6$ galaxy spectra in the BOSS of the Sloan Digital Sky Survey III. These systems consist of massive galaxies at redshifts of approximately 0.5 strongly lensing LAEs at redshifts from 2--3. The compact nature of LAEs makes them an ideal probe of dark substructures, with a substructure-mass sensitivity that is unprecedented in other optical strong-lens samples. The magnification effect from lensing will also reveal the structure of LAEs below 100 pc scales, providing a detailed look at the sites of the most concentrated unobscured star formation in the universe. The source code used for candidate selection is available for download as a part of this release.

preprint2016arXiv

The BOSS Emission-Line Lens Survey. IV. : Smooth Lens Models for the BELLS GALLERY Sample

We present \textsl{Hubble Space Telescope} (\textsl{HST}) F606W-band imaging observations of 21 galaxy-Ly$α$ emitter lens candidates in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$α$ EmitteR sYstems (BELLS GALLERY) survey. 17 systems are confirmed to be definite lenses with unambiguous evidence of multiple imaging. The lenses are primarily massive early-type galaxies (ETGs) at redshifts of approximately $0.55$, while the lensed sources are Ly$α$ emitters (LAEs) at redshifts from 2 to 3. Although the \textsl{HST} imaging data are well fit by smooth lens models consisting of singular isothermal ellipsoids in an external shear field, a thorough exploration of dark substructures in the lens galaxies is required. The Einstein radii of the BELLS GALLERY lenses are on average $60\%$ larger than those of the BELLS lenses because of the much higher source redshifts which will allow a detailed investigation of the radius evolution of the mass profile in ETGs. With the aid of the average $\sim 13 \times$ lensing magnification, the LAEs are resolved to comprise individual star-forming knots of a wide range of properties with characteristic sizes from less than 100 pc to several kpc, rest-frame far UV apparent AB magnitudes from 29.6 to 24.2, and typical projected separations of 500 pc to 2 kpc.

preprint2016arXiv

The Eruption of the Candidate Young Star ASASSN-15qi

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

preprint2016arXiv

The SDSS-III BOSS quasar lens survey: discovery of thirteen gravitationally lensed quasars

We report the discovery of 13 confirmed two-image quasar lenses from a systematic search for gravitationally lensed quasars in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We adopted a methodology similar to that used in the SDSS Quasar Lens Search (SQLS). In addition to the confirmed lenses, we report 11 quasar pairs with small angular separations ($\lesssim$2") confirmed from our spectroscopy, which are either projected pairs, physical binaries, or possibly quasar lens systems whose lens galaxies have not yet been detected. The newly discovered quasar lens system, SDSS J1452+4224 at zs$\approx$4.8 is one of the highest redshift multiply imaged quasars found to date. Furthermore, we have over 50 good lens candidates yet to be followed up. Owing to the heterogeneous selection of BOSS quasars, the lens sample presented here does not have a well-defined selection function.

preprint2016arXiv

The Swift AGN and Cluster Survey. II. Cluster Confirmation with SDSS Data

We study 203 (of 442) Swift AGN and Cluster Survey extended X-ray sources located in the SDSS DR8 footprint to search for galaxy over-densities in three dimensional space using SDSS galaxy photometric redshifts and positions near the Swift cluster candidates. We find 104 Swift clusters with a >3sigma galaxy over-density. The remaining targets are potentially located at higher redshifts and require deeper optical follow-up observations for confirmation as galaxy clusters. We present a series of cluster properties including the redshift, BCG magnitude, BCG-to-X-ray center offset, optical richness, and X-ray luminosity. We also detect red sequences in ~85% of the 104 confirmed clusters. The X-ray luminosity and optical richness for the SDSS confirmed Swift clusters are correlated and follow previously established relations. The distribution of the separations between the X-ray centroids and the most likely BCG is also consistent with expectation. We compare the observed redshift distribution of the sample with a theoretical model, and find that our sample is complete for z <~ 0.3 and is still 80% complete up to z ~= 0.4, consistent with the SDSS survey depth. These analysis results suggest that our Swift cluster selection algorithm has yielded a statistically well-defined cluster sample for further studying cluster evolution and cosmology. We also match our SDSS confirmed Swift clusters to existing cluster catalogs, and find 42, 23 and 1 matches in optical, X-ray and SZ catalogs, respectively, so the majority of these clusters are new detections.

preprint2016arXiv

The TDE ASASSN-14li and its host resolved at parsec scales with the EVN

We report European Very Long Baseline Interferometry Network (EVN) radio continuum observations of ASASSN-14li, one of the best studied tidal disruption events (TDEs) to date. At 1.7 GHz with ~12x6mas resolution, the emission is unresolved. At 5.0 GHz with ~3x2mas resolution, the radio emission shows an extended structure that can be modeled with two components: a core-like component and a fainter, possibly elongated source 4.3mas (~2pc) away. Our observations are not conclusive as to the nature of the components, but three scenarios are discussed. One possibility is a core-jet/outflow morphology, thus making of ASASSN-14li the first TDE jet/outflow directly imaged. For this case, the projected separation between the two components can only be explained by superluminal motion, rather than the lower velocities inferred from low-resolution radio observations. However, typical fast moving jets have brightness temperatures ~5 orders of magnitude higher than we find, thus making this scenario less likely. The second possibility is that we are imaging a non-relativistic jet from past AGN/TDE activity. In this case a past TDE is preferred given that the spatial extension and radio luminosity of the elongated component are consistent with the theoretical predictions for a TDE outflow. Alternatively, the two sources could indicate the presence of a binary black hole, which would then naturally explain the enhanced TDE rates of post-starburst galaxies. Future EVN observations will help us to distinguish between these scenarios.

preprint2015arXiv

A Connection between Obscuration and Star Formation in Luminous Quasars

We present a measurement of the star formation properties of a uniform sample of mid-IR selected, unobscured and obscured quasars (QSO1s and QSO2s) in the Boötes survey region. We use an spectral energy distribution (SED) analysis for photometric data spanning optical to far-IR wavelengths to decompose AGN and host galaxy components. We find that when compared to a matched sample of QSO1s, the QSO2s have higher far-IR detection fractions, far-IR fluxes and infrared star formation luminosities ($L_{\rm IR}^{\rm SF}$) by a factor of $\sim2$. Correspondingly, we show that the AGN obscured fraction rises from 0.3 to 0.7 between $4-40\times10^{11}L_\odot$. We also find evidence associating the absorption in the X-ray emission with the presence of far-IR emitting dust. Overall, these results are consistent with galaxy evolution models in which quasar obscurations can be associated with a dust-enshrouded starburst galaxies.

preprint2015arXiv

Rest-frame UV single-epoch black hole mass estimates of low-luminosity AGN at intermediate redshifts

The ability to accurately derive black hole (BH) masses at progressively higher redshifts and over a wide range of continuum luminosities has become indispensable in the era of large-area extragalactic spectroscopic surveys. In this paper we present an extension of existing comparisons between rest-frame UV and optical virial BH mass estimators to intermediate redshifts and luminosities comparable to the local H$β$ reverberation mapped active galactic nuclei (AGN). We focus on the MgII, CIV, and CIII] broad emission lines and compare them to both H$α$ and H$β$. We use newly acquired near-infrared spectra from the FMOS instrument on the Subaru telescope for 89 broad-lined AGN at redshifts between 0.3 and 3.5, complemented by data from the AGES survey. We employ two different prescriptions for measuring the emission line widths and compare the results. We confirm that MgII shows a tight correlation with H$α$ and H$β$, with a scatter of ~0.25 dex. The CIV and CIII] estimators, while showing larger scatter, are viable virial mass estimators after accounting for a trend with the UV-to-optical luminosity ratio. We find an intrinsic scatter of ~0.37 dex between Balmer and carbon virial estimators by combining our dataset with previous high redshift measurements. This updated comparison spans a total of 3 decades in BH mass. We calculate a virial factor for CIV/CIII] logf(CIV/CIII])=0.87 with an estimated systematic uncertainty of ~0.4 dex and find excellent agreement between the local reverberation mapped AGN sample and our high-z sample.

preprint2015arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. VII: The Seventh Year (2014-2015)

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 102 SU UMa-type dwarf novae observed mainly during the 2014-2015 season and characterized these objects. Our project has greatly improved the statistics of the distribution of orbital periods, which is a good approximation of the distribution of cataclysmic variables at the terminal evolutionary stage, and confirmed the presence of a period minimum at a period of 0.053 d and a period spike just above this period. The number density monotonically decreased toward the longer period and there was no strong indication of a period gap. We detected possible negative superhumps in Z Cha. It is possible that normal outbursts are also suppressed by the presence of a disk tilt in this system. There was no indication of enhanced orbital humps just preceding the superoutburst, and this result favors the thermal-tidal disk instability as the origin of superoutbursts. We detected superhumps in three AM CVn-type dwarf novae. Our observations and recent other detections suggest that 8% of objects showing dwarf nova-type outbursts are AM CVn-type objects. AM CVn-type objects and EI Psc-type object may be more abundant than previously recognized. OT J213806, a WZ Sge-type object, exhibited a remarkably different feature between the 2010 and 2014 superoutbursts. Although the 2014 superoutburst was much fainter the plateau phase was shorter than the 2010 one, the course of the rebrightening phase was similar. This object indicates that the O-C diagrams of superhumps can be indeed variable at least in WZ Sge-type objects. Four deeply eclipsing SU UMa-type dwarf novae (ASASSN-13cx, ASASSN-14ag, ASASSN-15bu, NSV 4618) were identified. We studied long-term trends in supercycles in MM Hya and CY UMa and found systematic variations of supercycles of ~20%.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

The SWIFT AGN and Cluster Survey I: Number Counts of AGN and Galaxy Clusters

The Swift AGN and Cluster Survey (SACS) uses 125 deg^2 of Swift XRT serendipitous fields with variable depths surrounding gamma-ray bursts to provide a medium depth (4e-15 erg/s/cm^2) and area survey filling the gap between deep, narrow Chandra/XMM-Newton surveys and wide, shallow ROSAT surveys. Here we present a catalog of 22,563 point sources and 442 extended sources and examine the number counts of the AGN and galaxy cluster populations. SACS provides excellent constraints on the AGN number counts at the bright end with negligible uncertainties due to cosmic variance, and these constraints are consistent with previous measurements. We use Wise mid-infrared (MIR) colors to classify the sources. For AGN we can roughly separate the point sources into MIR-red and MIR-blue AGN, finding roughly equal numbers of each type in the soft X-ray band (0.5-2 keV), but fewer MIR-blue sources in the hard X-ray band (2-8 keV). The cluster number counts, with 5% uncertainties from cosmic variance, are also consistent with previous surveys but span a much larger continuous flux range. Deep optical or IR follow-up observations of this cluster sample will significantly increase the number of higher redshift (z > 0.5) X-ray-selected clusters.

preprint2014arXiv

A UV to Mid-IR Study of AGN Selection

We classify the spectral energy distributions (SEDs) of 431,038 sources in the 9 sq. deg Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). There are up to 17 bands of data available per source, including ultraviolet (GALEX), optical (NDWFS), near-IR (NEWFIRM), and mid-infrared (IRAC/MIPS) data, as well as spectroscopic redshifts for ~20,000 objects, primarily from the AGN and Galaxy Evolution Survey (AGES). We fit galaxy, AGN, stellar, and brown dwarf templates to the observed SEDs, which yield spectral classes for the Galactic sources and photometric redshifts and galaxy/AGN luminosities for the extragalactic sources. The photometric redshift precision of the galaxy and AGN samples are sigma/(1+z)=0.040 and sigma/(1+z)=0.169, respectively, with the worst 5% outliers excluded. Based on the reduced chi-squared of the SED fit for each SED model, we are able to distinguish between Galactic and extragalactic sources for sources brighter than I=23.5. We compare the SED fits for a galaxy-only model and a galaxy+AGN model. Using known X-ray and spectroscopic AGN samples, we confirm that SED fitting can be successfully used as a method to identify large populations of AGN, including spatially resolved AGN with significant contributions from the host galaxy and objects with the emission line ratios of "composite" spectra. We also use our results to compare to the X-ray, mid-IR, optical color and emission line ratio selection techniques. For an F-ratio threshold of F>10 we find 16,266 AGN candidates brighter than I=23.5 and a surface density of ~1900 AGN per deg^2.

preprint2014arXiv

Discovery and Observations of ASASSN-13db, an EX Lupi-Type Accretion Event on a Low-Mass T Tauri Star

We discuss ASASSN-13db, an EX Lupi-type ("EXor") accretion event on the young stellar object (YSO) SDSS J051011.01$-$032826.2 (hereafter SDSSJ0510) discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN). Using archival photometric data of SDSSJ0510 we construct a pre-outburst spectral energy distribution (SED) and find that it is consistent with a low-mass class II YSO near the Orion star forming region ($d \sim 420$ pc). We present follow-up photometric and spectroscopic observations of the source after the $ΔV \sim-$5.4 magnitude outburst that began in September 2013 and ended in early 2014. These data indicate an increase in temperature and luminosity consistent with an accretion rate of $\sim10^{-7}$ $\rm{M}_\odot$ yr$^{-1}$, three or more orders of magnitude greater than in quiescence. Spectroscopic observations show a forest of narrow emission lines dominated by neutral metallic lines from Fe I and some low-ionization lines. The properties of ASASSN-13db are similar to those of the EXor prototype EX Lupi during its strongest observed outburst in late 2008.

preprint2014arXiv

New Limits On Gamma-Ray Emission From Galaxy Clusters

Galaxy clusters are predicted to produce gamma-rays through cosmic ray interactions and/or dark matter annihilation, potentially detectable by the Fermi Large Area Telescope (Fermi-LAT). We present a new, independent stacking analysis of Fermi-LAT photon count maps using the 78 richest nearby clusters (z<0.12) from the Two Micron All-Sky Survey (2MASS) cluster catalog. We obtain the lowest limit on the photon flux to date, 2.3e-11 ph/s/cm^2 (95% confidence) per cluster in the 0.8-100 GeV band, which corresponds to a luminosity limit of 3.5e44 ph/s. We also constrain the emission limits in a range of narrower energy bands. Scaling to recent cosmic ray acceleration and gamma-ray emission models, we find that cosmic rays represent a negligible contribution to the intra-cluster energy density and gas pressure.

preprint2014arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. VI: The Sixth Year (2013-2014)

Continuing the project described by Kato et al. (2009, PASJ, 61, S395, arXiv:0905.1757), we collected times of superhump maxima for 56 SU UMa-type dwarf novae mainly observed during the 2013-2014 season and characterized these objects. We detected negative superhumps in VW Hyi and indicated that the low number of normal outbursts in some supercycle can be interpreted as a result of the disk tilt. This finding, combined with the Kepler observation of V1504 Cyg and V344 Lyr, suggests that the disk tilt is responsible for modulating the outburst pattern in SU UMa-type dwarf novae. We also studied the deeply eclipsing WZ Sge-type dwarf nova MASTER OT J005740.99+443101.5 and found evidence of a sharp eclipse during the phase of early superhumps. The profile can be reproduced by a combination of the eclipse of the axisymmetric disk and the uneclipsed light source of early superhumps. This finding confirms the lack of evince of a greatly enhanced hot spot during the early stage of WZ Sge-type outburst. We detected growing (stage A) superhumps in MN Dra and give a suggestion that some of SU UMa-type dwarf novae situated near the critical condition of tidal instability may show long-lasting stage A superhumps. The large negative period derivatives reported in such systems can be understood a result of the combination of stage A and B superhumps. The WZ Sge-type dwarf novae AL Com and ASASSN-13ck showed a long-lasting (plateau-type) rebrightening. In the early phase of the rebrightening, both objects showed a precursor-like outburst, suggesting that the long-lasting rebrightening is triggered by a precursor outburst.

preprint2014arXiv

The Optical, Ultraviolet, and X-ray Structure of the Quasar HE 0435-1223

Microlensing has proven an effective probe of the structure of the innermost regions of quasars, and an important test of accretion disk models. We present light curves of the lensed quasar HE 0435-1223 in the R band and in the ultraviolet, and consider them together with X-ray light curves in two energy bands that are presented in a companion paper. Using a Bayesian Monte Carlo method, we constrain the size of the accretion disk in the rest-frame near- and far-UV, and constrain for the first time the size of the X-ray emission regions in two X-ray energy bands. The R-band scale size of the accretion disk is about 10^15.23 cm (~23 r_g), slightly smaller than previous estimates, but larger than would be predicted from the quasar flux. In the UV, the source size is weakly constrained, with a strong prior dependence. The UV to R-band size ratio is consistent with the thin disk model prediction, with large error bars. In soft and hard X-rays, the source size is smaller than ~10^14.8 cm (~10 r_g) at 95% confidence. We do not find evidence of structure in the X-ray emission region, as the most likely value for the ratio of the hard X-ray size to the soft X-ray size is unity. Finally, we find that the most likely value for the mean mass of stars in the lens galaxy is ~0.3 M_sun, consistent with other studies.

preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

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

preprint2013arXiv

A Two-Year Time Delay for the Lensed Quasar SDSS J1029+2623

We present 279 epochs of optical monitoring data spanning 5.4 years from 2007 January to 2012 June for the largest image separation (22.6 arcsec) gravitationally lensed quasar, SDSS J1029+2623. We find that image A leads the images B and C by dt_AB = (744+-10) days (90% confidence); the uncertainty includes both statistical uncertainties and systematic differences due to the choice of models. With only a ~1% fractional error, the interpretation of the delay is limited primarily by cosmic variance due to fluctuations in the mean line-of-sight density. We cannot separate the fainter image C from image B, but since image C trails image B by only 2-3 days in all models, the estimate of the time delay between image A and B is little affected by combining the fluxes of images B and C. There is weak evidence for a low level of microlensing, perhaps created by the small galaxy responsible for the flux ratio anomaly in this system. Interpreting the delay depends on better constraining the shape of the gravitational potential using the lensed host galaxy, other lensed arcs and the structure of the X-ray emission.

preprint2013arXiv

CIV Line-Width Anomalies: The Perils of Low S/N Spectra

Comparison of six high-redshift quasar spectra obtained with the Large Binocular Telescope with previous observations from the Sloan Digital Sky Survey shows that failure to correctly identify absorption and other problems with accurate characterization of the CIV emission line profile in low S/N data can severely limit the reliability of single-epoch mass estimates based on the CIV emission line. We combine the analysis of these new high-quality data with a reanalysis of three other samples based on high S/N spectra of the CIV emission line region. We find that a large scatter between the Hbeta- and CIV-based masses remains even for this high S/N sample when using the FWHM to characterize the BLR velocity dispersion and the standard virial assumption to calculate the mass. However, we demonstrate that using high-quality data and the line dispersion to characterize the CIV line width leads to a high level of consistency between CIV- and Hbeta-based masses, with <0.3 dex of observed scatter, and an estimated ~0.2 dex intrinsic scatter, in the mass residuals.

preprint2013arXiv

Detection of Substructure in the Gravitationally Lensed Quasar MG0414+0534 using Mid-Infrared and Radio VLBI Observations

We present 11.2 micron observations of the gravitationally lensed, radio-loud z_s=2.64 quasar MG0414+0534, obtained using the Michelle camera on Gemini North. We find a flux ratio anomaly of A2/A1= 0.93 +/- 0.02 for the quasar images A1 and A2. When combined with the 11.7 micron measurements from Minezaki et al. (2009), the A2/A1 flux ratio is nearly 5-sigma from the expected ratio for a model based on the two visible lens galaxies. The mid-IR flux ratio anomaly can be explained by a satellite (substructure), 0.3" Northeast of image A2, as can the detailed VLBI structures of the jet produced by the quasar. When we combine the mid-IR flux ratios with high-resolution VLBI measurements, we find a best-fit mass between 10^(6.2) and 10^(7.5) M_sol inside the Einstein radius for a satellite substructure modeled as a singular isothermal sphere at the redshift of the main lens (z_l=0.96). We are unable to set an interesting limit on the mass to light ratio due to its proximity to the quasar image A2. While the observations used here were technically difficult, surveys of flux anomalies in gravitational lenses with the James Webb Space Telescope will be simple, fast, and should well constrain the abundance of substructure in dark matter haloes.

preprint2013arXiv

Dust, Gas, and Metallicities of Cosmologically Distant Lens Galaxies

We homogeneously analyzed the \chandra\ X-ray observations of 10 gravitational lenses, HE 0047-1756, QJ 0158-4325, SDSS 0246-0805, HE 0435-1223, SDSS 0924+0219, SDSS 1004+4112, HE 1104-1805, PG 1115+080, Q 1355-2257, and Q 2237+0305, to measure the differential X-ray absorption between images, the metallicity, and the dust-to-gas ratio of the lens galaxies. We detected differential absorption in all lenses except SDSS 0924+0219 and HE 1104-1805. This doubles the sample of dust-to-gas ratio measurements in cosmologically distant lens galaxies. We successfully measured the gas phase metallicity of three lenses, Q 2237+0305, SDSS 1004+4112, and B 1152+199 from the X-ray spectra. Our results suggest a linear correlation between metallicity and dust-to-gas ratio (i.e., a constant metal-to-dust ratio), consistent with what is found for nearby galaxies. We obtain an average dust-to-gas ratio $E(B-V)/N_H=1.17^{+0.41}_{-0.31} \times 10^{-22}\rm mag\,cm^2\,atom^{-1}$ in the lens galaxies, with an intrinsic scatter of $\rm0.3\,dex$. Combining these results with data from GRB afterglows and quasar foreground absorbers, we found a mean dust-to-gas ratio $\mdtg,$ now significantly lower than the average Galactic value, $1.7\,\times 10^{-22}\,\rm mag\, cm^{2}\, atoms^{-1}.$ This suggests evolution of dust-to-gas ratios with redshift and lower average metallicities for the higher redshift galaxies, consistent with current metal and dust evolution models of interstellar medium. The slow evolution in the metal-to-dust ratio with redshift implies very rapid dust formation in high redshift ($z>2$) galaxies.

preprint2013arXiv

Mid-Infrared Selection of Active Galactic Nuclei with the Wide-Field Infrared Survey Explorer. II. Properties of WISE-Selected Active Galactic Nuclei in the NDWFS Boötes Field

Stern et al.(2012) presented a study of WISE selection of AGN in the 2 deg^2 COSMOS field, finding that a simple criterion W1-W2>=0.8 provides a highly reliable and complete AGN sample for W2<15.05, where the W1 and W2 passbands are centered at 3.4 and 4.6 microns, respectively. Here we extend this study using the larger 9 deg^2 NOAO Deep Wide-Field Survey Bootes field which also has considerably deeper WISE observations than the COSMOS field, and find that this simple color-cut significantly loses reliability at fainter fluxes. We define a modified selection criterion combining the W1-W2 color and the W2 magnitude to provide highly reliable or highly complete AGN samples for fainter WISE sources. In particular, we define a color-magnitude cut that finds 130+/-4 deg^-2 AGN candidates for W2<17.11 with 90% reliability. Using the extensive UV through mid-IR broad-band photometry available in this field, we study the spectral energy distributions of WISE AGN candidates. As expected, the WISE AGN selection is biased towards objects where the AGN dominates the bolometric luminosity output, and that it can identify highly obscured AGN. We study the distribution of reddening in the AGN sample and discuss a formalism to account for sample incompleteness based on the step-wise maximum-likelihood method of Efstathiou et al.(1988). The resulting dust obscuration distributions depend strongly on AGN luminosity, consistent with the trend expected for a Simpson (2005) receding torus. At L_AGN~3x10^44 erg/s, 29+/-7% of AGN are observed as Type 1, while at ~4x10^45 erg/s the fraction is 64+/-13%. The distribution of obscuration values suggests that dust in the torus is present as both a diffuse medium and in optically thick clouds.

preprint2013arXiv

Reconciling 56Ni Production in Type Ia Supernovae with Double Degenerate Scenarios

Binary white dwarf (WD) coalescences driven by gravitational waves or collisions in triple systems are potential progenitors of Type Ia supernovae (SNe Ia). We combine the distribution of 56Ni inferred from observations of SNe Ia with the results of both sub-Chandrasekhar detonation models and direct collision calculations to estimate what mass WDs should be exploding in each scenario to reproduce the observations. These WD mass distributions are then compared with the observed Galactic WD mass distribution and Monte Carlo simulations of WD-WD binary populations. For collisions, we find that the average mass of the individual components of the WD-WD binary must be peaked at ~0.75Msun, significantly higher than the average WD mass in binaries or in the field of ~0.55-0.60Msun. Thus, if collisions produce a large fraction of SNe Ia, then a mechanism must exist that favors large mass WDs. On the other hand, in an old stellar population, collisions would naturally result in a class of low luminosity SNe Ia, and we suggest these may be related to 1991bg-like events. For sub-Chandrasekhar detonations, we find that the average mass of the exploding WDs must be peaked at ~1.1Msun. This is interestingly similar to the average sum of the masses in WD-WD binaries, but it is not clear (and should be further explored) whether double degenerate mergers would be sufficiently efficient at synthesizing 56Ni to match the observed yields. If not, then actual ~1.1Msun WDs would be needed for sub-Chandrasekhar detonations. Since such high mass WDs are produced relatively quickly in comparison to the age of the environments where SNe Ia are found, this would require either accretion onto lower mass WDs prior to ignition or a long timescale between formation of the ~1.1Msun WD and ignition (such as set by gravitational wave emission or binary interactions).

preprint2013arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. V: The Fifth Year (2012-2013)

Continuing the project described by Kato et al. (2009a, arXiv:0905.1757), we collected times of superhump maxima for SU UMa-type dwarf novae mainly observed during the 2012-2013 season. We found three objects (V444 Peg, CSS J203937 and MASTER J212624) having strongly positive period derivatives despite the long orbital period (Porb). By using the period of growing stage (stage A) superhumps, we obtained mass ratios for six objects. We characterized nine new WZ Sge-type dwarf novae. We made a pilot survey of the decline rate of slowly fading part of SU UMa-type and WZ Sge-type outbursts. The decline time scale was found to generally follow the expected Porb^(1/4) dependence and WZ Sge-type outbursts also generally follow this trend. There are some objects which show slower decline rates, and we consider these objects good candidates for period bouncers. We also studied unusual behavior in some objects, including BK Lyn which made a transition from an ER UMa-type state to the novalike (standstill) state in 2013 and unusually frequent occurrence of superoutbursts in NY Ser and CR Boo. We applied least absolute shrinkage and selection operator (Lasso) power spectral analysis, which has been proven to be very effective in analyzing the Kepler data, to ground-based photometry of BK Lyn and detected the dramatic disappearance of the signal of negative superhumps in 2013. We suggested that the mass-transfer rates did not vary strongly between the ER UMa-type state and novalike state in BK Lyn, and this transition was less likely caused by a systematic variation of the mass-transfer rate.

preprint2013arXiv

The Impact of Metallicity on the Rate of Type Ia Supernovae

The metallicity of a star strongly affects both its evolution and the properties of the stellar remnant that results from its demise. It is generally accepted that stars with initial masses below ~8 M_sun leave behind white dwarfs and that some sub-population of these lead to Type Ia supernovae. However, it is often tacitly assumed that metallicity has no effect on the rate of SNe Ia. We propose that a consequence of the effects of metallicity is to significantly increase the SN Ia rate in lower-metallicity galaxies, in contrast to previous expectations. This is because lower-metallicity stars leave behind higher-mass white dwarfs, which should be easier to bring to explosion. We first model SN Ia rates in relation to galaxy masses and ages alone, finding that the elevation in the rate of SNe Ia in lower-mass galaxies measured by LOSS is readily explained. However, we then see that models incorporating this effect of metallicity agree just as well. Using the same parameters to estimate the cosmic SN Ia rate, we again find good agreement with data up to z~2. We suggest that this degeneracy warrants more detailed examination of host galaxy metallicities. We discuss additional implications, including for hosts of high-z SNe Ia, the SN Ia delay time distribution, super-Chandrasekhar SNe, and cosmology.

preprint2013arXiv

The Magellanic Quasars Survey. III. Spectroscopic Confirmation of 758 AGNs Behind the Magellanic Clouds

The Magellanic Quasars Survey (MQS) has now increased the number of quasars known behind the Magellanic Clouds by almost an order of magnitude. All survey fields in the Large Magellanic Cloud (LMC) and 70% of those in the Small Magellanic Cloud (SMC) have been observed. The targets were selected from the third phase of the Optical Gravitational Lensing Experiment (OGLE-III) based on their optical variability, mid-IR and/or X-ray properties. We spectroscopically confirmed 758 (565 LMC and 193 SMC) quasars behind the Clouds, of which 94% (527 LMC and 186 SMC) are newly identified. The MQS quasars have long-term (12 years and growing for OGLE), high-cadence light curves, enabling unprecedented variability studies of quasars. The MQS quasars also provide a dense reference grid for measuring both the internal and bulk proper motions of the Clouds, and 50 quasars are bright enough (I<18 mag) for absorption studies of the interstellar/galactic (ISM/IGM) medium of the Clouds.

preprint2013arXiv

The Structure of HE 1104-1805 from Infrared to X-Ray

The gravitationally lensed quasar HE 1104-1805 has been observed at a variety of wavelengths ranging from the mid-infrared to X-ray for nearly 20 years. We combine flux ratios from the literature, including recent Chandra data, with new observations from the SMARTS telescope and HST, and use them to investigate the spatial structure of the central regions using a Bayesian Monte Carlo analysis of the microlensing variability. The wide wavelength coverage allows us to constrain not only the accretion disk half-light radius r_1/2, but the power-law slope ξ of the size-wavelength relation r_1/2 ~ λ^ξ. With a logarithmic prior on the source size, the (observed-frame) R-band half-light radius log(r_1/2/cm) is 16.0+0.3-0.4, and the slope ξ is 1.0+0.30-0.56. We put upper limits on the source size in soft (0.4-1.2 keV) and hard (1.2-8 keV) X-ray bands, finding 95% upper limits on log (r_1/2/cm) of 15.33 in both bands. A linear prior yields somewhat larger sizes, particularly in the X-ray bands. For comparison, the gravitational radius, using a black hole mass estimated using the Hβ line, is log(r_g/cm) = 13.94. We find that the accretion disk is probably close to face-on, with cos i = 1.0 being four times more likely than cos i = 0.5. We also find probability distributions for the mean mass of the stars in the foreground lensing galaxy, the direction of the transverse peculiar velocity of the lens, and the position angle of the projected accretion disk's major axis (if not face-on).

preprint2012arXiv

A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

preprint2012arXiv

A Global Physical Model for Cepheids

We perform a global fit to ~5,000 radial velocity and ~177,000 magnitude measurements in 29 photometric bands covering 0.3 to 8.0 microns distributed among 287 Galactic, LMC, and SMC Cepheids with P > 10 days. We assume that the Cepheid light curves and radial velocities are fully characterized by distance, reddening, and time-dependent radius and temperature variations. We construct phase curves of radius and temperature for periods between 10 and 100 days, which yield light curve templates for all our photometric bands and can be easily generalized to any additional band. With only 4 to 6 parameters per Cepheid, depending on the existence of velocity data and the amount of freedom in the distance, the models have typical rms light and velocity curve residuals of 0.05 mag and 3.5 km/s. The model derives the mean Cepheid spectral energy distribution and its derivative with respect to temperature, which deviate from a black body in agreement with metal-line and molecular opacity effects. We determine a mean reddening law towards the Cepheids in our sample, which is not consistent with standard assumptions in either the optical or near-IR. Based on stellar atmosphere models we predict the biases in distance, reddening, and temperature determinations due to the metallicity and we quantify the metallicity signature expected for our fit residuals. The observed residuals as a function of wavelength show clear differences between the individual galaxies, which are compatible with these predictions. In particular, we find that metal-poor Cepheids appear hotter. Finally, we provide a framework for optimally selecting filters that yield the smallest overall errors in Cepheid parameter determination, or filter combinations for suppressing or enhancing the metallicity effects on distance determinations. We make our templates publicly available.

preprint2012arXiv

A measurement of the intrahalo light fraction with near-infrared background anisotropies

Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.

preprint2012arXiv

Characterizing the Optical Variability of Bright Blazars: Variability-Based Selection of Fermi AGN

We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.

preprint2012arXiv

Further Evidence that Quasar X-Ray Emitting Regions Are Compact: X-Ray and Optical Microlensing in the Lensed Quasar Q J0158-4325

We present four new seasons of optical monitoring data and six epochs of X-ray photometry for the doubly-imaged lensed quasar Q J0158-4325. The high-amplitude, short-period microlensing variability for which this system is known has historically precluded a time delay measurement by conventional methods. We attempt to circumvent this limitation by application of a Monte Carlo microlensing analysis technique, but we are only able to prove that the delay must have the expected sign (image A leads image B). Despite our failure to robustly measure the time delay, we successfully model the microlensing at optical and X-ray wavelengths to find a half light radius for soft X-ray emission log(r_{1/2,X,soft}/cm) = 14.3^{+0.4}_{-0.5}, an upper limit on the half-light radius for hard X-ray emission log(r_{1/2,X,hard}/cm) <= 14.6 and a refined estimate of the inclination-corrected scale radius of the optical R-band (rest frame 3100 Angstrom) continuum emission region of log(r_s/cm) = 15.6+-0.3.

preprint2012arXiv

The BOSS Emission-Line Lens Survey. II. Investigating Mass-Density Profile Evolution in the SLACS+BELLS Strong Gravitational Lens Sample

We present an analysis of the evolution of the central mass-density profile of massive elliptical galaxies from the SLACS and BELLS strong gravitational lens samples over the redshift interval z ~ 0.1-0.6, based on the combination of strong-lensing aperture mass and stellar velocity-dispersion constraints. We find a significant trend towards steeper mass profiles (parameterized by the power-law density model with rho ~ r^[-gamma]) at later cosmic times, with magnitude d<gamma>/dz = -0.60 +/- 0.15. We show that the combined lens-galaxy sample is consistent with a non-evolving distribution of stellar velocity dispersions. Considering possible additional dependence of <gamma> on lens-galaxy stellar mass, effective radius, and Sersic index, we find marginal evidence for shallower mass profiles at higher masses and larger sizes, but with a significance that is sub-dominant to the redshift dependence. Using the results of published Monte Carlo simulations of spectroscopic lens surveys, we verify that our mass-profile evolution result cannot be explained by lensing selection biases as a function of redshift. Interpreted as a true evolutionary signal, our result suggests that major dry mergers involving off-axis trajectories play a significant role in the evolution of the average mass-density structure of massive early-type galaxies over the past 6 Gyr. We also consider an alternative non-evolutionary hypothesis based on variations in the strong-lensing measurement aperture with redshift, which would imply the detection of an "inflection zone" marking the transition between the baryon-dominated and dark-matter halo-dominated regions of the lens galaxies. Further observations of the combined SLACS+BELLS sample can constrain this picture more precisely, and enable a more detailed investigation of the multivariate dependences of galaxy mass structure across cosmic time.

preprint2012arXiv

The Chandra view of the Largest Quasar Lens SDSS J1029+2623

We present results from Chandra observations of the cluster lens SDSS J1029+2623 at z_l=0.58, which is a gravitationally lensed quasar with the largest known image separation. We clearly detect X-ray emission both from the lensing cluster and the three lensed quasar images. The cluster has an X-ray temperature of kT = 8.1 (+2.0, -1.2) keV and bolometric luminosity of L_X = 9.6e44 erg s^-1. Its surface brightness is centered near one of the brightest cluster galaxies, and it is elongated East-West. We identify a subpeak North-West of the main peak, which is suggestive of an ongoing merger. Even so, the X-ray mass inferred from the hydrostatic equilibrium assumption appears to be consistent with the lensing mass from the Einstein radius of the system. We find significant absorption in the soft X-ray spectrum of the faintest quasar image, which can be caused by an intervening material at either the lens or source redshift. The X-ray flux ratios between the quasar images (after correcting for absorption) are in reasonable agreement with those at optical and radio wavelengths, and all the flux ratios are inconsistent with those predicted by simple mass models. This implies that microlensing effect is not significant for this system and dark matter substructure is mainly responsible for the anomalous flux ratios.

preprint2012arXiv

The Galaxy Optical Luminosity Function from the AGN and Galaxy Evolution Survey (AGES)

We present the galaxy optical luminosity function for the redshift range 0.05<z<0.75 from the AGN and Galaxy Evolution Survey (AGES), a spectroscopic survey of 7.6 sq. deg. in the Bootes field of the NOAO Deep Wide-Field Survey. Our statistical sample is comprised of 12,473 galaxies with known redshifts down to I=20.4 (AB). Our results at low redshift are consistent with those from SDSS; at higher redshift, we find strong evidence for evolution in the luminosity function, including differential evolution between blue and red galaxies. We find that the luminosity density evolves as (1+z)^(0.54+/-0.64) for red galaxies and (1+z)^(1.64+/-0.39) for blue galaxies.

preprint2012arXiv

The Hidden Fortress: Structure and substructure of the complex strong lensing cluster SDSS J1029+2623

We present Hubble Space Telescope (HST) ACS and WFC3 observations of SDSS J1029+2623, a three-image quasar lens system produced by a foreground cluster at z=0.584. Our strong lensing analysis reveals 6 additional multiply imaged galaxies. We confirm the complex nature of the mass distribution of the lensing cluster, with a bimodal distribution which deviates from the Chandra X-ray SB distribution. The Einstein radius is estimated to be θ_E=15.2" \pm 0.5" for the quasar redshift of z=2.197. We derive a radial mass distribution from the combination of strong lensing, HST/ACS weak lensing, and Subaru/Suprime-cam weak lensing analysis results, finding a best-fit virial mass of M_vir=(1.5+0.40-0.35) \times 10^14 h^-1 M_sun and a concentration parameter of c_vir=25.7+14.1-7.5. The lensing mass estimate at the outer radius is smaller than the X-ray mass estimate by a factor of ~2. We ascribe this large mass discrepancy to shock heating of the intracluster gas during a merger, which is also suggested by the complex mass and gas distributions and the high value of the concentration parameter. In the HST image, we also identify a probable galaxy, GX, in the vicinity of the faintest quasar image C. The inclusion of GX explains the anomalous flux ratios between the quasar images. The morphology of the highly elongated quasar host galaxy is also well reproduced. The best-fit model suggests large total magnifications of 30 for the quasar and 35 for the quasar host galaxy, and has an AB time delay consistent with the measured value. (Abridged)

preprint2012arXiv

The observed neutron star mass distribution as a probe of the supernova explosion mechanism

The observed distribution of neutron star (NS) masses reflects the physics of core-collapse supernova explosions and the structure of the massive stars that produce them at the end of their evolution. We present a Bayesian analysis that directly compares the NS mass distribution observed in double NS systems to theoretical models of NS formation. We find that models with standard binary mass ratio distributions are strongly preferred over independently picking the masses from the initial mass function, although the strength of the inference depends on whether current assumptions for identifying the remnants of the primary and secondary stars are correct. Second, NS formation models with no mass fallback are favored because they reduce the dispersion in NS masses. The double NS system masses thus directly point to the mass coordinate where the supernova explosion was initiated, making them an excellent probe of the supernova explosion mechanism. If we assume no fallback and simply vary the mass coordinate separating the remnant and the supernova ejecta, we find that for solar metallicity stars the explosion most likely develops at the edge of the iron core at a specific entropy of about 2.8 k_B. The primary limitations of our study are the poor knowledge of the supernova explosion mechanism and the lack of broad range of SN model explosions of LMC to solar metallicity.

preprint2012arXiv

The Sloan Digital Sky Survey Quasar Lens Search. V. Final Catalog from the Seventh Data Release

We present the final statistical sample of lensed quasars from the Sloan Digital Sky Survey (SDSS) Quasar Lens Search (SQLS). The well-defined statistical lens sample consists of 26 lensed quasars brighter than i=19.1 and in the redshift range of 0.6<z<2.2 selected from 50,836 spectroscopically confirmed quasars in the SDSS Data Release 7 (DR7), where we restrict the image separation range to 1"<θ<20" and the i-band magnitude differences in two image lenses to be smaller than 1.25 mag. The SDSS DR7 quasar catalog also contains 36 additional lenses identified with various techniques. In addition to these lensed quasars, we have identified 81 pairs of quasars from follow-up spectroscopy, 26 of which are physically associated binary quasars. The statistical lens sample covers a wide range of image separations, redshifts, and magnitudes, and therefore is suitable for systematic studies of cosmological parameters and surveys of the structure and evolution of galaxies and quasars.

preprint2012arXiv

The Sloan Digital Sky Survey Quasar Lens Search. VI. Constraints on Dark Energy and the Evolution of Massive Galaxies

We present a statistical analysis of the final lens sample from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). The number distribution of a complete subsample of 19 lensed quasars selected from 50,836 source quasars is compared with theoretical expectations, with particular attention to the selection function. Assuming that the velocity function of galaxies does not evolve with redshift, the SQLS sample constrains the cosmological constant to Ω_Λ=0.79^{+0.06}_{-0.07}(stat.)^{+0.06}_{-0.06}(syst.) for a flat universe. The dark energy equation of state is found to be consistent with w=-1 when the SQLS is combined with constraints from baryon acoustic oscillation (BAO) measurements or results from the Wilkinson Microwave Anisotropy Probe (WMAP). We also obtain simultaneous constraints on cosmological parameters and redshift evolution of the galaxy velocity function, finding no evidence for redshift evolution at z<1 in any combinations of constraints. For instance, number density evolution quantified as ν_n=d\lnϕ_*/d\ln(1+z) and the velocity dispersion evolution ν_σ=d\lnσ_*/d\ln(1+z) are constrained to ν_n=1.06^{+1.36}_{-1.39}(stat.)^{+0.33}_{-0.64}(syst.) and ν_σ=-0.05^{+0.19}_{-0.16}(stat.)^{+0.03}_{-0.03}(syst.) respectively when the SQLS result is combined with BAO and WMAP for flat models with a cosmological constant. We find that a significant amount of dark energy is preferred even after fully marginalizing over the galaxy evolution parameters. Thus the statistics of lensed quasars robustly confirm the accelerated cosmic expansion.

preprint2012arXiv

The Structure of 2MASS Galaxy Clusters

We use a sample of galaxies from the Two Micron All Sky Survey (2MASS) Extended Source Catalog to refine a matched filter method of finding galaxy clusters that takes into account each galaxy's position, magnitude, and redshift if available. The matched filter postulates a radial density profile, luminosity function, and line-of-sight velocity distribution for cluster galaxies. We use this method to search for clusters in the galaxy catalog, which is complete to an extinction-corrected K-band magnitude of 13.25 and has spectroscopic redshifts for roughly 40% of the galaxies, including nearly all brighter than K = 11.25. We then use a stacking analysis to determine the average luminosity function, radial distribution, and velocity distribution of cluster galaxies in several richness classes, and use the results to update the parameters of the matched filter before repeating the cluster search. We also investigate the correlations between a cluster's richness and its velocity dispersion and core radius, using these relations to refine priors that are applied during the cluster search process. After the second cluster search iteration, we repeat the stacking analysis. We find a cluster galaxy luminosity function that fits a Schechter form, with parameters M_K* - 5 log h = -23.64\pm0.04 and α= -1.07\pm0.03. We can achieve a slightly better fit to our luminosity function by adding a Gaussian component on the bright end to represent the brightest cluster galaxy (BCG) population. The radial number density profile of galaxies closely matches a projected Navarro-Frenk-White (NFW) profile at intermediate radii, with deviations at small radii due to well-known cluster centering issues and outside the virial radius due to correlated structure. The velocity distributions are Gaussian in shape, with velocity dispersions that correlate strongly with richness.

preprint2012arXiv

X-ray Monitoring of Gravitational Lenses With Chandra

We present \emph{Chandra} monitoring data for six gravitationally lensed quasars: QJ 0158$-$4325, HE 0435$-$1223, HE 1104$-$1805, SDSS 0924+0219, SDSS 1004+4112, and Q 2237+0305. We detect X-ray microlensing variability in all six lenses with high confidence. We detect energy dependent microlensing in HE 0435$-$1223, SDSS 1004+4112, SDSS 0924+0219 and Q 2237+0305. We present a detailed spectral analysis for each lens, and find that simple power-law models plus Gaussian emission lines give good fits to the spectra. We detect intrinsic spectral variability in two epochs of Q 2237+0305. We detect differential absorption between images in four lenses. We also detect the \feka\ emission line in all six lenses, and the Ni XXVII K$α$ line in two images of Q 2237+0305. The rest frame equivalent widths of the \feka\ lines are measured to be 0.4--1.2 keV, significantly higher than those measured in typical active galactic nuclei of similar X-ray luminosities. This suggests that the \feka\ emission region is more compact or centrally concentrated than the continuum emission region.

preprint2011arXiv

A Study of Gravitational Lens Chromaticity using Ground-based Narrow Band Photometry

We present observations of wavelength-dependent flux ratios for four gravitational lens systems (SDSS~J1650+4251, HE~0435$-$1223, FBQ 0951+2635, and Q~0142$-$100) obtained with the Nordic Optical Telescope. The use of narrowband photometry, as well as the excellent seeing conditions during the observations, allows us to study their chromatic behavior. For SDSS~J1650+4251, we determine the extinction curve of the dust in the $z_L=0.58$ lens galaxy and find that the 2175 Å\ feature is absent. In the case of HE~0435$-$1223, we clearly detect chromatic microlensing. This allows us to estimate the wavelength-dependent size of the accretion disk. We find an R-band disk size of $r^{R}_s=13\pm5$ light days for a linear prior on $r^{R}_s$ and $r^{R}_s=7\pm6$ light days for a logarithmic prior. For a power-law size-wavelength scaling of $r_s\proptoλ^{p}$, we were able to constrain the value of the exponent to $p=1.3\pm0.3$ for both $r^{R}_s$ priors, which is in agreement with the temperature profiles of simple thin disk models ($p=4/3$).

preprint2011arXiv

A Study of Gravitational Lens Chromaticity with the Hubble Space Telescope

We report Hubble Space Telescope observations of 6 gravitational lenses with the Advanced Camera for Surveys. We measured the flux ratios between the lensed images in 6 filters from 8140Å to 2200Å. In 3 of the systems, HE0512$-$3329, B1600+434, and H1413+117, we were able to construct UV extinction curves partially overlapping the 2175Å feature and characterize the properties of the dust relative to the Galaxy and the Magellanic Clouds. In HE1104$-$1804 we detect chromatic microlensing and use it to study the physical properties of the quasar accretion disk. For a Gaussian model of the disk $\exp(-r^2/2 r_s^2)$, scaling with wavelength as $r_s \propto λ^p$, we estimate $r_s(\lambda3363)=4^{+4}_{-2}$ ($7\pm 4$) light-days and $p=1.1\pm 0.6$ ($1.0\pm 0.6$) for a logarithmic (linear) prior on $r_s$. The remaining two systems, FBQ0951+2635 and SBS1520+530, yielded no useful estimates of extinction or chromatic microlensing.

preprint2011arXiv

Analyzing the Flux Anomalies of the Large-Separation Lensed Quasar SDSS J1029+2623

Using a high resolution radio image, we successfully resolve the two fold image components B and C of the quasar lens system SDSS J1029+2623. The flux anomalies associated with these two components in the optical regime persist, albeit less strongly, in our radio observations, suggesting that the cluster must be modeled by something more than a single central potential. We argue that placing substructure close to one of the components can account for a flux anomaly with negligible changes in the component positions. Our best fit model has a substructure mass of ~10^8 solar masses up to the mass-sheet degeneracy, located roughly 0.1 arcsecs West and 0.1 arcsecs North of component B. We demonstrate that a positional offset between the centers of the source components can explain the differences between the optical and radio flux ratios.

preprint2011arXiv

Clustering of obscured and unobscured quasars in the Bootes field: Placing rapidly growing black holes in the cosmic web

We present the first measurement of the spatial clustering of mid-infrared selected obscured and unobscured quasars, using a sample in the redshift range 0.7 < z < 1.8 selected from the 9 deg^2 Bootes multiwavelength survey. Recently the Spitzer Space Telescope and X-ray observations have revealed large populations of obscured quasars that have been inferred from models of the X-ray background and supermassive black hole evolution. To date, little is known about obscured quasar clustering, which allows us to measure the masses of their host dark matter halos and explore their role in the cosmic evolution of black holes and galaxies. In this study we use a sample of 806 mid-infrared selected quasars and ~250,000 galaxies to calculate the projected quasar-galaxy cross-correlation function w_p(R). The observed clustering yields characteristic dark matter halo masses of log (M_halo [h^-1 M_sun]) = 12.7^+0.4_-0.6 and 13.3^+0.3_-0.4 for unobscured quasars (QSO-1s) and obscured quasars (Obs-QSOs), respectively. The results for QSO-1s are in excellent agreement with previous measurements for optically-selected quasars, while we conclude that the Obs-QSOs are at least as strongly clustered as the QSO-1s. We test for the effects of photometric redshift errors on the optically-faint Obs-QSOs, and find that our method yields a robust lower limit on the clustering; photo-z errors may cause us to underestimate the clustering amplitude of the Obs-QSOs by at most ~20%. We compare our results to previous studies, and speculate on physical implications of stronger clustering for obscured quasars.

preprint2011arXiv

Discovery of Energy Dependent X-ray Microlensing in Q2237+0305

We present our long term Chandra X-ray monitoring data for the gravitationally lensed quasar Q2237+0305 with 20 epochs spanning 10 years. We easily detect microlensing variability between the images in the full (0.2--8 keV), soft (0.2--2 keV), and hard (2--8 keV) bands at very high confidence. We also detect, for the first time, chromatic microlensing differences between the soft and hard X-ray bands. The hard X-ray band is more strongly microlensed than the soft band, suggesting that the corona above the accretion disk thought to generate the X-rays has a non-uniform electron distribution, in which the hotter and more energetic electrons occupy more compact regions surrounding the black holes. Both the hard and soft X-ray bands are more strongly microlensed than the optical (restframe UV) emission, indicating that the X-ray emission is more compact than the optical, confirming the microlensing results from other lenses.

preprint2011arXiv

The BOSS Emission-Line Lens Survey (BELLS). I. A large spectroscopically selected sample of Lens Galaxies at redshift ~ 0.5

We present a catalog of 25 definite and 11 probable strong galaxy-galaxy gravitational lens systems with redshifts 0.4 \lesssim z \lesssim 0.7, discovered spectroscopically by the presence of higher redshift emission-lines within the Baryon Oscillation Spectroscopic Survey (BOSS) of luminous galaxies, and confirmed with high-resolution Hubble Space Telescope (HST) images of 44 candidates. Our survey extends the methodology of the Sloan Lens ACS Survey (SLACS: Bolton et al. 2006; 2008) to higher redshift. We describe the details of the BOSS spectroscopic candidate detections, our HST Adanced Camera for Surveys (ACS) image processing and analysis methods, and our strong gravitational lens modeling procedure. We report BOSS spectroscopic parameters and ACS photometric parameters for all candidates, and mass-distribution parameters for the best-fit singular isothermal ellipsoid models of definite lenses. Our sample to date was selected using only the first six months of BOSS survey-quality spectroscopic data. The full 5-year BOSS database should produce a sample of several hundred strong galaxy-galaxy lenses and in combination with SLACS lenses at lower redshift, strongly constrain the redshift evolution of the structure of elliptical, bulge-dominated galaxies as a function of luminosity, stellar mass, and rest-frame color, thereby providing a powerful test for competing theories of galaxy formation and evolution.

preprint2011arXiv

The Cosmic Core-collapse Supernova Rate does not match the Massive-Star Formation Rate

We identify a "supernova rate problem": the measured cosmic core-collapse supernova rate is a factor of ~ 2 smaller (with significance ~ 2 sigma) than that predicted from the measured cosmic massive-star formation rate. The comparison is critical for topics from galaxy evolution and enrichment to the abundance of neutron stars and black holes. We systematically explore possible resolutions. The accuracy and precision of the star formation rate data and conversion to the supernova rate are well supported, and proposed changes would have far-reaching consequences. The dominant effect is likely that many supernovae are missed because they are either optically dim (low-luminosity) or dark, whether intrinsically or due to obscuration. We investigate supernovae too dim to have been discovered in cosmic surveys by a detailed study of all supernova discoveries in the local volume. If possible supernova impostors are included, then dim supernovae are common enough by fraction to solve the supernova rate problem. If they are not included, then the rate of dark core collapses is likely substantial. Other alternatives are that there are surprising changes in our understanding of star formation or supernova rates, including that supernovae form differently in small galaxies than in normal galaxies. These possibilities can be distinguished by upcoming supernova surveys, star formation measurements, searches for disappearing massive stars, and measurements of supernova neutrinos.

preprint2011arXiv

The Effect of a Time-Varying Accretion Disk Size on Quasar Microlensing Light Curves

Microlensing perturbations to the magnification of gravitationally lensed quasar images are dependent on the angular size of the quasar. If quasar variability at visible wavelengths is caused by a change in the area of the accretion disk, it will affect the microlensing magnification. We derive the expected signal, assuming that the luminosity scales with some power of the disk area, and estimate its amplitude using simulations. We discuss the prospects for detecting the effect in real-world data and for using it to estimate the logarithmic slope of the luminosity's dependence on disk area. Such an estimate would provide a direct test of the standard thin accretion disk model. We tried fitting six seasons of the light curves of the lensed quasar HE 0435-1223 including this effect as a modification to the Kochanek et al. (2006) approach to estimating time delays. We find a dramatic improvement in the goodness of fit and relatively plausible parameters, but a robust estimate will require a full numerical calculation in order to correctly model the strong correlations between the structure of the microlensing magnification patterns and the magnitude of the effect. We also comment briefly on the effect of this phenomenon for the stability of time delay estimates.

preprint2011arXiv

The Magellanic Quasars Survey. I. Doubling the Number of Known AGNs Behind the Small Magellanic Cloud

We report the spectroscopic confirmation of 29 new, 12 plausible, and 3 previously known quasars behind the central ~1.5 deg^2 region of the Small Magellanic Cloud. These were identified in a single 2df/AAOmega observation on the Anglo-Australian Telescope of 268 candidates selected primarily based on their mid-IR colors, along with a smaller number of optically variable sources in OGLE-II close to known X-ray sources. The low detection efficiency was partly expected from the high surface density of SMC as compared to the LMC targets and the faintness of many of them (149 with I>20 mag). The expected number of I<20 mag quasars in the field is ~38, and we found 15 (22 with plausible) or 40% (60%). We did not attempt to determine the nature of the remaining sources, although several appear to be new planetary nebulae. The newly discovered AGNs can be used as reference points for future proper motion studies, to study absorption in the interstellar medium, and to study the physics of quasar variability with the existing long-term, highly cadenced OGLE light curves.

preprint2011arXiv

The Microlensing Properties of a Sample of 87 Lensed Quasars

Gravitational microlensing is a powerful tool for probing the physical properties of quasar accretion disks and properties of the lens galaxy such as its dark matter fraction and mean stellar mass. Unfortunately the number of lensed quasars ($\sim 90$) exceeds our monitoring capabilities. Thus, estimating their microlensing properties is important for identifying good microlensing candidates as well as for the expectations of future surveys. In this work we estimate the microlensing properties of a sample of 87 lensed quasars. While the median Einstein radius crossing time scale is 20.6 years, the median source crossing time scale is 7.3 months. Broadly speaking, this means that on $\sim 10$ year timescales roughly half the lenses will be quiescent, with the source in a broad demagnified valley, and roughly half will be active with the source lying in the caustic ridges. We also found that the location of the lens system relative to the CMB dipole has a modest effect on microlensing timescales, and in theory microlensing could be used to confirm the kinematic origin of the dipole. As a corollary of our study we analyzed the accretion rate parameters in a sub-sample of 32 lensed quasars. At fixed black hole mass, it is possible to sample a broad range of luminosities (i.e., Eddington factors) if it becomes feasible to monitor fainter lenses.

preprint2010arXiv

Eight New Quasar Lenses from the Sloan Digital Sky Survey Quasar Lens Search

We report the discovery and confirmation of eight new two-image lensed quasars by the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. The lenses are SDSSJ0904+1512 (image separation θ=1"13, source redshift z_s=1.826), SDSSJ1054+2733 (θ=1"27, z_s=1.452), SDSSJ1055+4628 (θ=1"15, z_s=1.249), SDSSJ1131+1915 (θ=1"46, z_s=2.915), SDSSJ1304+2001 (θ=1"87, z_s=2.175), SDSSJ1349+1227 (θ=3"00, z_s=1.722), SDSSJ1455+1447 (θ=1"73, z_s=1.424), and SDSSJ1620+1203 (θ=2"77, z_s=1.158). Three of them, SDSSJ1055+4628, SDSSJ1455+1447, and SDSSJ1620+1203, satisfy the criteria for constructing our statistical sample for studying the cosmological model. Based on galactic absorption lines of the lens galaxies, we also derive lens redshifts of z_l=0.398 and z_l=0.513 for SDSSJ1620+1203 and the previously discovered lens SDSSJ0746+4403, respectively.

preprint2010arXiv

Microlensing Evidence That A Type 1 Quasar Is Viewed Face On

Using a microlensing analysis of 11-years of OGLE V-band photometry of the four image gravitational lens Q2237+0305, we measure the inclination i of the accretion disk to be cos i > 0.66 at 68% confidence. Very edge on (cos i < 0.39) solutions are ruled out at 95% confidence. We measure the V-band radius of the accretion disk, defined by the radius where the temperature matches the monitoring band photon emission, to be R_V = 5.8 (+3.8) (-2.3) x 10^15 cm assuming a simple thin disk model and including the uncertainties in its inclination. The projected radiating area of the disk remains too large to be consistent with the observed flux for a T proportional to R^-3/4 thin disk temperature profile. There is no strong correlation between the direction of motion (peculiar velocity) of the lens galaxy and the orientation of the disk.

preprint2010arXiv

Mid-Infrared Variability from the Spitzer Deep, Wide-Field Survey

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

preprint2010arXiv

On the Baryon Fractions in Clusters and Groups of Galaxies

We present the baryon fractions of 2MASS groups and clusters as a function of cluster richness using total and gas masses measured from stacked ROSAT X-ray data and stellar masses estimated from the infrared galaxy catalogs. We detect X-ray emission even in the outskirts of clusters, beyond r_200 for richness classes with X-ray temperatures above 1 keV. This enables us to more accurately determine the total gas mass in these groups and clusters. We find that the optically selected groups and clusters have flatter temperature profiles and higher stellar-to-gas mass ratios than the individually studied, X-ray bright clusters. We also find that the stellar mass in poor groups with temperatures below 1 keV is comparable to the gas mass in these systems. Combining these results with individual measurements for clusters, groups, and galaxies from the literature, we find a break in the baryon fraction at ~1 keV. Above this temperature, the baryon fraction scales with temperature as f_b \propto T^0.20\pm0.03. We see significantly smaller baryon fractions below this temperature, and the baryon fraction of poor groups joins smoothly onto that of systems with still shallower potential wells such as normal and dwarf galaxies where the baryon fraction scales with the inferred velocity dispersion as f_b \propto σ^1.6. The small scatter in the baryon fraction at any given potential well depth favors a universal baryon loss mechanism and a preheating model for the baryon loss. The scatter is, however, larger for less massive systems. Finally, we note that although the broken power-law relation can be inferred from data points in the literature alone, the consistency between the baryon fractions for poor groups and massive galaxies inspires us to fit the two categories of objects (galaxies and clusters) with one relation.

preprint2010arXiv

The Evolution of the Star Formation Rate of Galaxies at 0.0 < z < 1.2

We present the 24 micron rest-frame luminosity function (LF) of star-forming galaxies in the redshift range 0.0 < z < 0.6 constructed from 4047 spectroscopic redshifts from the AGN and Galaxy Evolution Survey of 24 micron selected sources in the Bootes field of the NOAO Deep Wide-Field Survey. This sample provides the best available combination of large area (9 deg^2), depth, and statistically complete spectroscopic observations, allowing us to probe the evolution of the 24 micron LF of galaxies at low and intermediate redshifts while minimizing the effects of cosmic variance. In order to use the observed 24 micron luminosity as a tracer for star formation, active galactic nuclei (AGNs) that could contribute significantly at 24 micron are identified and excluded from our star-forming galaxy sample based on their mid-IR spectral energy distributions or the detection of X-ray emission. The evolution of the 24 micron LF of star-forming galaxies for redshifts of z < 0.65 is consistent with a pure luminosity evolution where the characteristic 24 micron luminosity evolves as (1+z)^(3.8+/-0.3). We extend our evolutionary study to encompass 0.0 < z < 1.2 by combining our data with that of the Far-Infrared Deep Extragalactic Legacy Survey. Over this entire redshift range the evolution of the characteristic 24 micron luminosity is described by a slightly shallower power law of (1+z)^(3.4+/-0.2). We find a local star formation rate density of (1.09+/-0.21) x 10^-2 Msun/yr/Mpc^-3, and that it evolves as (1+z)^(3.5+/-0.2) over 0.0 < z < 1.2. These estimates are in good agreement with the rates using optical and UV fluxes corrected for the effects of intrinsic extinction in the observed sources. This agreement confirms that star formation at z <~ 1.2 is robustly traced by 24 micron observations and that it largely occurs in obscured regions of galaxies. (Abridged)

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Lens Search. IV. Statistical Lens Sample from the Fifth Data Release

We present the second report of our systematic search for strongly lensed quasars from the data of the Sloan Digital Sky Survey (SDSS). From extensive follow-up observations of 136 candidate objects, we find 36 lenses in the full sample of 77,429 spectroscopically confirmed quasars in the SDSS Data Release 5. We then define a complete sample of 19 lenses, including 11 from our previous search in the SDSS Data Release 3, from the sample of 36,287 quasars with i<19.1 in the redshift range 0.6<z<2.2, where we require the lenses to have image separations of 1"<θ<20" and i-band magnitude differences between the two images smaller than 1.25 mag. Among the 19 lensed quasars, 3 have quadruple-image configurations, while the remaining 16 show double images. This lens sample constrains the cosmological constant to be Ω_Λ=0.84^{+0.06}_{-0.08}(stat.)^{+0.09}_{-0.07}(syst.) assuming a flat universe, which is in good agreement with other cosmological observations. We also report the discoveries of 7 binary quasars with separations ranging from 1.1" to 16.6", which are identified in the course of our lens survey. This study concludes the construction of our statistical lens sample in the full SDSS-I data set.

preprint2010arXiv

The Transverse Peculiar Velocity of the Q2237+0305 Lens Galaxy and the Mean Mass of Its Stars

Using 11-years of OGLE V-band photometry of Q2237+0305, we measure the transverse velocity of the lens galaxy and the mean mass of its stars. We can do so because, for the first time, we fully include the random motions of the stars in the lens galaxy in the analysis of the light curves. In doing so, we are also able to correctly account for the Earth's parallax motion and the rotation of the lens galaxy, further reducing systematic errors. We measure a lower limit on the transverse speed of the lens galaxy, v_t > 338 km/s (68% confidence) and find a preferred direction to the East. The mean stellar mass estimate including a well-defined velocity prior is 0.12 <= <M/Msun> <= 1.94 at 68% confidence, with a median of 0.52 Msun. We also show for the first time that analyzing subsets of a microlensing light curve, in this case the first and second halves of the OGLE V-band light curve, give mutually consistent physical results.

preprint2009arXiv

Detection of a Companion Lens Galaxy using the Mid-infrared Flux Ratios of the Gravitationally Lensed Quasar H1413+117

We present the first resolved mid-IR (11 micron) observations of the four-image quasar lens H1413+117 using the Michelle camera on Gemini North. All previous observations (optical, near-IR, and radio) of this lens show a "flux anomaly," where the image flux ratios cannot be explained by a simple, central lens galaxy. We attempt to reproduce the mid-IR flux ratios, which are insensitive to extinction and microlensing, by modeling the main lens as a singular isothermal ellipsoid. This model fails to reproduce the flux ratios. However, we can explain the flux ratios simply by adding to the model a nearby galaxy detected in the H-band by HST/NICMOS-NIC2. This perturbing galaxy lies 4.0" from the main lens and it has a critical radius of 0.63" +/- 0.02" which is similar to that of the main lens, as expected from their similar H-band fluxes. More remarkably, this galaxy is not required to obtain a good fit to the system astrometry, so this represents the first clear detection of an object through its effect on the image fluxes of a gravitational lens. This is a parallel to the detections of visible satellites from astrometric anomalies, and provides a proof of the concept of searching for substructure in galaxies using anomalous flux ratios.

preprint2009arXiv

Evolution of the UV Excess in Early-Type Galaxies

We examine the UV emission from luminous early-type galaxies as a function of redshift. We perform a stacking analysis using Galaxy Evolution Explorer (GALEX) images of galaxies in the NOAO Deep Wide Field Survey (NDWFS) Boötes field and examine the evolution in the UV colors of the average galaxy. Our sample, selected to have minimal ongoing star formation based on the optical to mid-IR SEDs of the galaxies, includes 1843 galaxies spanning the redshift range $0.05\leq z\leq0.65$. We find evidence that the strength of the UV excess decreases, on average, with redshift, and our measurements also show moderate disagreement with previous models of the UV excess. Our results show little evolution in the shape of the UV continuum with redshift, consistent either with the binary model for the formation of Extreme Horizontal Branch (EHB) stars or with no evolution in EHB morphology with look-back time. However, the binary formation model predicts that the strength of the UV excess should also be relatively constant, in contradiction with our measured results. Finally, we see no significant influence of a galaxy's environment on the strength of its UV excess.

preprint2009arXiv

Host galaxies, clustering, Eddington ratios, and evolution of radio, X-ray, and infrared-selected AGNs

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the radio, X-ray, and infrared. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts from the AGN and Galaxy Evolution Survey (AGES) and data in the radio (WSRT 1.4 GHz), X-rays (Chandra XBootes), and mid-IR (IRAC Shallow Survey). The radio, X-ray, and IR AGN samples show modest overlap, indicating that to the flux limits of the survey, they represent largely distinct classes of AGNs. We derive host galaxy colors and luminosities, as well as Eddington ratios (lambda), for obscured or optically faint AGNs. We also measure the two-point cross-correlation between AGNs and galaxies on scales of 0.3-10 h^-1 Mpc, and derive typical dark matter halo masses. We find that: (1) radio AGNs are mainly found in luminous red galaxies, are strongly clustered (with M_halo ~ 3x10^13 h^-1 M_sun), and have very low lambda <~ 10^-3; (2) X-ray-selected AGNs are preferentially found in galaxies in the "green valley" of color-magnitude space and are clustered similarly to typical AGES galaxies (M_halo ~ 10^13 h^-1 M_sun), with 10^-3 <~ lambda <~ 1; (3) IR AGNs reside in slightly bluer, less luminous galaxies than X-ray AGNs, are weakly clustered (M_halo <~ 10^12 h^-1 M_sun), and have lambda > 10^-2. We interpret these results in terms of a simple model of AGN and galaxy evolution, whereby a "quasar" phase and the growth of the stellar bulge occurs when a galaxy's dark matter halo reaches a critical mass between ~10^12 and 10^13 M_sun. Subsequently, star formation ceases and AGN accretion shifts from radiatively efficient (optical- and IR- bright) to radiatively inefficient (optically-faint, radio-bright) modes.

preprint2009arXiv

Quantifying Quasar Variability As Part of a General Approach To Classifying Continuously Varying Sources

Robust fast methods to classify variable light curves in large sky surveys are becoming increasingly important. While it is relatively straightforward to identify common periodic stars and particular transient events (supernovae, novae, microlensing), there is no equivalent for non-periodic continuously varying sources (quasars, aperiodic stellar variability). In this paper we present a fast method for modeling and classifying such sources. We demonstrate the method using ~ 86,000 variable sources from the OGLE-II survey of the LMC and ~ 2,700 mid-IR selected quasar candidates from the OGLE-III survey of the LMC and SMC. We discuss the location of common variability classes in the parameter space of the model. In particular we show that quasars occupy a distinct region of variability space, providing a simple quantitative approach to the variability selection of quasars.

preprint2009arXiv

The Star Formation and Nuclear Accretion Histories of Normal Galaxies in the AGES Survey

We combine IR, optical and X-ray data from the overlapping, 9.3 square degree NOAO Deep Wide-Field Survey (NDWFS), AGN and Galaxy Evolution Survey (AGES), and XBootes Survey to measure the X-ray evolution of 6146 normal galaxies as a function of absolute optical luminosity, redshift, and spectral type over the largely unexplored redshift range 0.1 < z < 0.5. Because only the closest or brightest of the galaxies are individually detected in X-rays, we use a stacking analysis to determine the mean properties of the sample. Our results suggest that X-ray emission from spectroscopically late-type galaxies is dominated by star formation, while that from early-type galaxies is dominated by a combination of hot gas and AGN emission. We find that the mean star formation and supermassive black hole accretion rate densities evolve like (1+z)^3, in agreement with the trends found for samples of bright, individually detectable starburst galaxies and AGN. Our work also corroborates the results of many previous stacking analyses of faint source populations, with improved statistics.

preprint2008arXiv

Spectroscopic Confirmation of the Fifth Image of SDSS J1004+4112 and Implications for the M_BH-sigma_* Relation at z=0.68

We present the results of deep spectroscopy for the central region of the cluster lens SDSS J1004+4112 with the Subaru telescope. A secure detection of an emission line of the faint blue stellar object (component E) near the center of the brightest cluster galaxy (G1) confirms that it is the central fifth image of the lensed quasar system. In addition, we measure the stellar velocity dispersion of G1 to be sigma_* = 352+-13 km/s. We combine these results to obtain constraints on the mass M_BH of the putative black hole (BH) at the center of the inactive galaxy G1, and hence on the M_BH-sigma_* relation at the lens redshift z_l=0.68. From detailed mass modeling, we place an upper limit on the black hole mass, M_BH < 2.1x10^{10}M_sun at 1-sigma level (<3.1x10^{10}M_sun at 3-sigma), which is consistent with black hole masses expected from the local and redshift-evolved M_BH-sigma_* relations, M_BH~10^{9}-10^{10}M_sun.

preprint2008arXiv

The Origin of the 24-micron Excess in Red Galaxies

Observations with the Spitzer Space Telescope have revealed a population of red-sequence galaxies with a significant excess in their 24-micron emission compared to what is expected from an old stellar population. We identify 900 red galaxies with 0.15<z<0.3 from the AGN and Galaxy Evolution Survey (AGES) selected from the NOAO Deep Wide-Field Survey Bootes field. Using Spitzer/MIPS, we classify 89 (~10%) with 24-micron infrared excess (f24>0.3 mJy). We determine the prevalence of AGN and star-formation activity in all the AGES galaxies using optical line diagnostics and mid-IR color-color criteria. Using the IRAC color-color diagram from the IRAC Shallow Survey, we find that 64% of the 24-micron excess red galaxies are likely to have strong PAH emission features in the 8-micron IRAC band. This fraction is significantly larger than the 5% of red galaxies with f24<0.3 mJy that are estimated to have strong PAH emission, suggesting that the infrared emission is largely due to star-formation processes. Only 15% of the 24-micron excess red galaxies have optical line diagnostics characteristic of star-formation (64% are classified as AGN and 21% are unclassifiable). The difference between the optical and infrared results suggest that both AGN and star-formation activity is occurring simultaneously in many of the 24-micron excess red galaxies. These results should serve as a warning to studies that exclusively use optical line diagnostics to determine the dominant emission mechanism in the infrared and other bands. We find that ~40% of the 24-micron excess red galaxies are edge-on spiral galaxies with high optical extinctions. The remaining sources are likely to be red galaxies whose 24-micron emission comes from a combination of obscured AGN and star-formation activity.

preprint2007arXiv

The X-ray Properties of Optically-Selected Galaxy Clusters

We stacked the X-ray data from the ROSAT All Sky Survey for over 4,000 clusters selected from the 2MASS catalog and divided into five richness classes. We detected excess X-ray emission over background at the center of the stacked images in all five richness bins. The interrelationships between the mass, X-ray temperature and X-ray luminosity of the stacked clusters agree well with those derived from catalogs of X-ray clusters. Poisson variance in the number of galaxies occupying halos of a given mass leads to significant differences between the average richness at fixed mass and the average mass at fixed richness that we can model relatively easily using a simple model of the halo occupation distribution. These statistical effects probably explain recent results in which optically-selected clusters lie on the same X-ray luminosity-temperature relations as local clusters but have lower optical richnesses than observed for local clusters with the same X-ray properties. When we further binned the clusters by redshift, we did not find significant redshift-dependent biases in the sense that the X-ray luminosities for massive clusters of fixed optical richness show little dependence on redshift beyond that expected from the effects of Poisson fluctuations. Our results demonstrate that stacking of RASS data from optically selected clusters can be a powerful test for biases in cluster selection algorithms.

preprint2006arXiv

Black Hole Masses and Eddington Ratios at 0.3<z<4

We study the distribution of Eddington luminosity ratios, L_bol/L_edd, of active galactic nuclei (AGNs) discovered in the AGN and Galaxy Evolution Survey (AGES). We combine H-beta, MgII, and CIV line widths with continuum luminosities to estimate black hole (BH) masses in 407 AGNs, covering the redshift range z~0.3-4 and the bolometric luminosity range L_bol~10^45-10^47 erg/s. The sample consists of X-ray or mid-infrared (24 micron) point sources with optical magnitude R<=21.5 mag and optical emission line spectra characteristic of AGNs. For the range of luminosity and redshift probed by AGES, the distribution of estimated Eddington ratios is well described as log-normal with a peak at L_bol/L_edd ~ 1/4 and a dispersion of 0.3 dex. Since additional sources of scatter are minimal, this dispersion must account for contributions from the scatter between estimated and true BH mass and the scatter between estimated and true bolometric luminosity. Therefore, we conclude that: (1) neither of these sources of error can contribute more than \~0.3 dex rms; and (2) the true Eddington ratios of optically luminous AGNs are even more sharply peaked. Because the mass estimation errors must be smaller than ~0.3 dex, we can also investigate the distribution of Eddington ratios at fixed BH mass. We show for the first time that the distribution of Eddington ratios at fixed BH mass is peaked, and that the dearth of AGNs at a factor ~10 below Eddington is real and not an artifact of sample selection. These results provide strong evidence that supermassive BHs gain most of their mass while radiating close to the Eddington limit, and they suggest that the fueling rates in luminous AGNs are ultimately determined by BH self-regulation of the accretion flow rather than galactic scale dynamical disturbances.

preprint2005arXiv

SDSS J0806+2006 and SDSS J1353+1138: Two New Gravitationally Lensed Quasars from the Sloan Digital Sky Survey

We report the discoveries of two, two-image gravitationally lensed quasars selected from the Sloan Digital Sky Survey: SDSS J0806+2006 at z_s=1.540 and SDSS J1353+1138 at z_s=1.629 with image separations of 1.40" and 1.41" respectively. Spectroscopic and optical/near-infrared imaging follow-up observations show that the quasar images have identical redshifts and possess extended objects between the images that are likely to be lens galaxies at z_l~0.6 in SDSS J0806+2006 and z_l~0.3 in SDSS J1353+1138. The field of SDSS J0806+2006 contains several nearby galaxies that may significantly perturb the system, and SDSS J1353+1138 has an extra component near its Einstein ring that is probably a foreground star. Simple mass models with reasonable parameters reproduce the quasar positions and fluxes of both systems.

preprint2005arXiv

The XBootes Chandra Survey Paper II: The X-ray Source Catalog

We present results from a Chandra survey of the nine square degree Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). This XBootes survey consists of 126 separate contiguous ACIS-I observations each of approximately 5000 seconds in duration. These unique Chandra observations allow us to search for large scale structure and to calculate X-ray source statistics o ver a wide, contiguous field of view with arcsecond angular resolution and uniform coverage. Optical spectroscopic follow-up observations and the rich NDWFS data set will allow us to identify and classify these X-ray selected sources. Using wavelet decomposition, we detect 4642 point sources with n $\ge$ 2 counts. In order to keep our detections ~99% reliable, we limit our list to sources with n $\ge$ 4 counts. The full 0.5--7 keV band n $\ge$ 4 count list has 3293 point sources. In addition to the point sources, 43 extended sources have been detected consistent, with the depth of these observations and the number counts of clusters. We present here the X-ray catalog for the XBootes survey, including source positions, X-ray fluxes, hardness ratios and their uncertainties. We calculate and present the differential number of sources per flux densit y interval, $N(S)$, for the point sources.

preprint2003arXiv

The central image of a gravitationally lensed quasar

A galaxy can act as a gravitational lens, producing multiple images of a background object. Theory predicts there should be an odd number of images but, paradoxically, almost all observed lenses have 2 or 4 images. The missing image should be faint and appear near the galaxy's center. These ``central images'' have long been sought as probes of galactic cores too distant to resolve with ordinary observations. There are five candidates, but in one case the third image is not necessarily a central image, and in the others, the central component might be a foreground source rather than a lensed image. Here we report the most secure identification of a central image, based on radio observations of PMN J1632-0033, one of the latter candidates. Lens models incorporating the central image show that the mass of the lens galaxy's central black hole is less than 2 x 10^8 M_sun, and the galaxy's surface density at the location of the central image is more than 20,000 M_sun per square parsec, in agreement with expectations based on observations of galaxies hundreds of times closer to the Earth.

preprint1999arXiv

Redshifts of the Gravitational Lenses MG1131+0456 and B1938+666

The redshifts of the gravitational lens galaxies in MG1131+0456 and B1938+666 are 0.844 and 0.881 respectively. Both are early-type galaxies lying at the redshifts predicted by assuming that they are early-type galaxies with old stellar populations lying on the fundamental plane. We also find evidence for a foreground group of galaxies at z=0.343 near MG1131+0456. The source redshifts are predicted to be >1.8 in both systems, but they are so red that infrared spectra will be required to determine their redshifts.

preprint1996arXiv

Structural and dynamical uncertainties in modeling axisymmetric elliptical galaxies

Quantitative dynamical models of galaxies require deprojecting the observed surface brightness to determine the luminosity density of the galaxy. Existing deprojection methods for axisymmetric galaxies assume that a unique deprojection exists for any given inclination, even though the projected density is known to be degenerate to the addition of "konus densities" that are invisible in projection. We develop a deprojection method based on linear regularization that can explore the range of luminosity densities statistically consistent with an observed surface brightness distribution. The luminosity density is poorly constrained at modest inclinations (i > ~30 deg), even in the limit of vanishing observational errors. In constant mass-to-light ratio, axisymmetric, two-integral dynamical models, the uncertainties in the luminosity density result in large uncertainties in the meridional plane velocities. However, the projected line-of-sight velocities show variations comparable to current typical observational uncertainties.