Source author record

C. S. Kochanek

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

148works
9topics
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

148 published item(s)

preprint2023arXiv

Optical/$γ$-ray blazar flare correlations: understanding the high-energy emission process using ASAS-SN and Fermi light curves

Using blazar light curves from the optical All-Sky Automated Survey for Supernovae (ASAS-SN) and the $γ$-ray \textit{Fermi}-LAT telescope, we performed the most extensive statistical correlation study between both bands, using a sample of 1,180 blazars. This is almost an order of magnitude larger than other recent studies. Blazars represent more than 98\% of the AGNs detected by \textit{Fermi}-LAT and are the brightest $γ$-ray sources in the extragalactic sky. They are essential for studying the physical properties of astrophysical jets from central black holes. However, their $γ$-ray flare mechanism is not fully understood. Multi-wavelength correlations help constrain the dominant mechanisms of blazar variability. We search for temporal relationships between optical and $γ$-ray bands. Using a Bayesian Block Decomposition, we detect 1414 optical and 510 $γ$-ray flares, we find a strong correlation between both bands. Among all the flares, we find 321 correlated flares from 133 blazars, and derive an average rest-frame time delay of only 1.1$_{-8.5}^{+7.1}$ days, with no difference between the flat-spectrum radio quasars, BL Lacertae-like objects or low, intermediate, and high-synchrotron peaked blazar classes. Our time-delay limit rules out the hadronic proton-synchrotron model as the driver for non-orphan flares and suggests a leptonic single-zone model. Limiting our search to well-defined light curves and removing 976 potential but unclear ``orphan'' flares, we find 191 (13\%) and 115 (22\%) clear ``orphan'' optical and $γ$-ray flares. The presence of ``orphan'' flares in both bands challenges the standard one-zone blazar flare leptonic model and suggests multi-zone synchrotron sites or a hadronic model for some blazars.

preprint2022arXiv

A Rapid Ionization Change in the Nebular-Phase Spectra of the Type Ia SN 2011fe

We present three new spectra of the nearby Type Ia supernova (SN Ia) 2011fe covering $\approx 480-850~$days after maximum light and show that the ejecta undergoes a rapid ionization shift at $\sim 500~$days after explosion. The prominent [FeIII] emission lines at $\approx 4600~$Å are replaced with the permitted FeI+FeII blends at $\sim 4400~$Å and $\sim 5400~$Å. The $\approx 7300~$Å feature, which is produced by [FeII]+[NiII] at $\lesssim 400~$days after explosion, is replaced by broad ($\approx \pm 15\,000~\rm{km}~\rm{s}^{-1}$) symmetric [CaII] emission. Models predict this ionization transition occurring $\sim 100$ days later than what is observed, which we attribute to clumping in the ejecta. Finally, we use the nebular-phase spectra to test several proposed progenitor scenarios for SN2011fe. Non-detections of H and He exclude nearby non-degenerate companions, [OI] non-detections disfavor the violent merger of two white dwarfs, and the symmetric emission-line profiles favor a symmetric explosion.

preprint2022arXiv

An X-ray View of the Ambiguous Nuclear Transient AT2019pev

AT2019pev is a nuclear transient in a narrow-line Seyfert 1 galaxy at $z=0.096$. The archival ultraviolet, optical and infrared data showed features of both tidal disruption events (TDEs) and active galactic nuclei (AGNs), and its nature is not fully understood. We present detailed X-ray observations of AT2019pev taken with Swift, Chandra and NICER over 173 days of its evolution since the first Swift XRT epoch. The X-ray luminosity increases by a factor of five in five days from the first Swift XRT epoch to the lightcurve peak. The lightcurve decays by a factor of ten over $\sim$75 days and then flattens with a weak re-brightening trend at late times. The X-ray spectra show a "harder-when-brighter" trend before peak and a "harder-when-fainter" trend after peak, which may indicate a transition of accretion states. The archival ground-based optical observations show similar time evolution as the X-ray lightcurves. Beyond the seasonal limit of the ground-based observations, the Gaia lightcurve is rising toward an equally bright or brighter peak 223 days after the optical discovery. Combining our X-ray analysis and archival multi-wavelength data, AT2019pev more closely resembles an AGN transient.

preprint2022arXiv

The "Giraffe": Discovery of a stripped red giant in an interacting binary with a ${\sim}2~M_\odot$ lower giant

We report the discovery of a stripped giant + lower giant binary, 2M04123153+6738486 (2M0412), identified during a search for non-interacting compact object-star binaries. 2M0412 is an evolved ($T_{\rm eff, giant}\simeq4000$ K), luminous ($L_{\rm giant}\simeq150~L_\odot$) red giant in a circular $P=81.2$ day binary. 2M0412 is a known variable star previously classified as a semi-regular variable. The cross-correlation functions of follow-up Keck/HIRES and LBT/PEPSI spectra show an RV-variable second component with implied mass ratio $q=M_{\rm giant}/M_{\rm comp}\simeq0.20\pm0.01$. The ASAS-SN, ATLAS, TESS and ZTF light curves show that the giant is a Roche lobe filling ellipsoidal variable with an inclination of $49.4^\circ{}\pm{0.3^{\circ}}$, and a giant mass of $M_{\rm giant}=0.38\pm0.01~ M_\odot$ for a distance of $\simeq3.7$ kpc. The mass of the giant indicates that its envelope has been stripped. The giant companion on the lower red giant branch has a mass of $M_{\rm comp}=1.91\pm0.03~M_\odot$ with $T_{\rm eff, comp}\simeq5000$ K, $L_{\rm comp}\simeq60~L_\odot$ and $R_{\rm comp}\simeq11~R_\odot$. We also identify an orbital phase dependent, broad $\rm Hα$ emission line which could indicate ongoing accretion from the stripped red giant onto the companion.

preprint2022arXiv

The ASAS-SN Catalog of Variable Stars X: Discovery of 116,000 New Variable Stars Using g-band Photometry

The All-Sky Automated Survey for Supernovae (ASAS-SN) is the first optical survey to monitor the entire sky, currently with a cadence of $\lesssim 24$ hours down to $g \lesssim 18.5$ mag. ASAS-SN has routinely operated since 2013, collecting $\sim$ 2,000 to over 7,500 epochs of $V$ and $g-$band observations per field to date. This work illustrates the first analysis of ASAS-SN's newer, deeper, higher cadence $g-$band data. From an input source list of ${\sim}55$ million isolated sources with $g<18$~mag, we identified $1.5\times10^6$ variable star candidates using a random forest classifier trained on features derived from $\textit{Gaia}$, 2MASS, and AllWISE. Using ASAS-SN $g-$band light curves, and an updated random forest classifier augmented with data from Citizen ASAS-SN, we classified the candidate variables into 8 broad variability types. We present a catalog of ${\sim}116,000$ new variable stars with high classification probabilities, including ${\sim}111,000$ periodic variables and ${\sim}5,000$ irregular variables. We also recovered ${\sim}263,000$ known variable stars.

preprint2022arXiv

The search for failed supernovae with the Large Binocular Telescope: a new candidate and the failed SN fraction with 11 yr of data

We present updated results of the Large Binocular Telescope Search for Failed Supernovae. This search monitors luminous stars in 27 nearby galaxies with a current baseline of 11~yr of data. We re-discover the failed supernova (SN) candidate N6946-BH1 as well as a new candidate, M101-OC1. M101-OC1 is a blue supergiant that rapidly disappears in optical wavelengths with no evidence for significant obscuration by warm dust. While we consider other options, a good explanation for the fading of M101-OC1 is a failed SN, but follow-up observations are needed to confirm this. Assuming only one clearly detected failed SN, we find a failed SN fraction $f = 0.16^{+0.23}_{-0.12}$ at 90 per~cent confidence. We also report on a collection of stars that show slow ($\sim$decade), large amplitude ($ΔL/L > 3$) luminosity changes.

preprint2022arXiv

To TDE or not to TDE: The luminous transient ASASSN-18jd with TDE-like and AGN-like qualities

We present the discovery of ASASSN-18jd (AT 2018bcb), a luminous optical/UV/X-ray transient located in the nucleus of the galaxy 2MASX J22434289--1659083 at $z=0.1192$. Over the year after discovery, Swift UVOT photometry shows the UV SED of the transient to be well modeled by a slowly shrinking blackbody with temperature $T \sim 2.5 \times 10^{4} \rm ~K$, a maximum observed luminosity of $L_\text{max} = 4.5^{+0.6}_{-0.3} \times 10^{44} \rm ~erg ~s^{-1}$, and a radiated energy of $E = 9.6^{+1.1}_{-0.6} \times 10^{51} \rm ~erg$. X-ray data from Swift XRT and XMM-Newton show a transient, variable X-ray flux with blackbody and power-law components that fade by nearly an order of magnitude over the following year. Optical spectra show strong, roughly constant broad Balmer emission as well as transient features attributable to He II, N III-V, O III, and coronal Fe. While ASASSN-18jd shares similarities with Tidal Disruption Events (TDEs), it is also similar to the newly-discovered nuclear transients seen in quiescent galaxies and faint Active Galactic Nuclei (AGNs).

preprint2022arXiv

Using AGN lightcurves to map accretion disc temperature fluctuations

We introduce a new model for understanding AGN continuum variability. We start from a Shakura--Sunyaev thin accretion disc with a steady-state radial temperature profile $T(R)$ and assume that the variable flux is due to axisymmetric temperature perturbations $δT(R,t)$. After linearizing the equations, we fit UV-optical AGN lightcurves to determine $δT(R,t)$ for a sample of seven AGNs. We see a diversity of $|δT/T| \sim 0.1$ fluctuation patterns which are not dominated by outgoing waves traveling at the speed of light as expected for the "lamppost" model used to interpret disc reverberation mapping studies. Rather, the most common pattern resembles slow ($v \ll c$) ingoing waves. An explanation for our findings is that these ingoing waves trigger central temperature fluctuations that act as a lamppost, producing lower amplitude temperature fluctuations moving outwards at the speed of light. The lightcurves are dominated by the lamppost signal -- even though the temperature fluctuations are dominated by other structures with similar variability time-scales -- because the discs exponentially smooth the contributions from the slower moving ($v \ll c$) fluctuations to the observed lightcurves. This leads to lightcurves that closely resemble the expectations for a lamppost model but with the slow variability time-scales of the ingoing waves. This also implies that longer time-scale variability signals will increasingly diverge from lamppost models because the smoothing of slower-moving waves steadily decreases as their period or spatial wavelength increases.

preprint2021arXiv

An Extreme Amplitude, Massive Heartbeat System in the LMC Characterized Using ASAS-SN and TESS

Using ASAS-SN data, we find that the bright ($V{\sim}13.5$ mag) variable star MACHO 80.7443.1718 (ASASSN-V J052624.38-684705.6) is the most extreme heartbeat star yet discovered. This massive binary, consisting of at least one early B-type star, has an orbital period of $P_{\rm ASAS-SN}=32.83627\pm0.00846\,{\rm d},$ and is located towards the LH58 OB complex in the LMC. Both the ASAS-SN and TESS light curves show extreme brightness variations of ${\sim}40\%$ at periastron and variations of ${\sim}10\%$ due to tidally excited oscillations outside periastron. We fit an analytical model of the variability caused by the tidal distortions at pericenter to find orbital parameters of $ω=-61.4$ deg, $i=44.8$ deg and $e=0.566$. We also present a frequency analysis to identify the pulsation frequencies corresponding to the tidally excited oscillations.

preprint2021arXiv

Citizen ASAS-SN Data Release I: Variable Star Classification Using Citizen Science

We present the first results from Citizen ASAS-SN, a citizen science project for the All-Sky Automated Survey for Supernovae (ASAS-SN) hosted on the Zooniverse platform. Citizen ASAS-SN utilizes the newer, deeper, higher cadence ASAS-SN $g$-band data and tasks volunteers to classify periodic variable star candidates based on their phased light curves. We started from 40,640 new variable candidates from an input list of ${\sim} 7.4$ million stars with $δ< -60^\circ$ and the volunteers identified 10,420 new discoveries which they classified as 4,234 pulsating variables, 3,132 rotational variables, 2,923 eclipsing binaries, and 131 variables flagged as Unknown. They classified known variable stars with an accuracy of 89% for pulsating variables, 81% for eclipsing binaries, and 49% for rotational variables. We examine user performance, agreement between users, and compare the citizen science classifications with our machine learning classifier updated for the $g$-band light curves. In general, user activity correlates with higher classification accuracy and higher user agreement. We used the user's "Junk" classifications to develop an effective machine learning classifier to separate real from false variables, and there is a clear path for using this "Junk" training set to significantly improve our primary machine learning classifier. We also illustrate the value of Citizen ASAS-SN for identifying unusual variables with several examples.

preprint2021arXiv

Citizen ASAS-SN: Citizen Science with The All-Sky Automated Survey for SuperNovae (ASAS-SN)

We present "Citizen ASAS-SN", a citizen science project hosted on the Zooniverse platform which utilizes data from the All-Sky Automated Survey for SuperNovae (ASAS-SN). Volunteers are presented with ASAS-SN $g$-band light curves of variable star candidates. The classification workflow allows volunteers to classify these sources into major variable groups, while also allowing for the identification of unique variable stars for additional follow-up.

preprint2021arXiv

Discovery of a Highly Eccentric, Chromospherically Active Binary: ASASSN-V J192114.84+624950.8

As part of an All-Sky Automated Survey for SuperNovae (ASAS-SN) search for sources with large flux decrements, we discovered a transient where the quiescent, stellar source, ASASSN-V J192114.84+624950.8, rapidly decreased in flux by $\sim55\%$ ($\sim0.9$ mag) in the g-band. The \textit{TESS} light curve revealed that the source is a highly eccentric, eclipsing binary. Fits to the light curve using \textsc{phoebe} find the binary orbit to have $e=0.79$, $P_{\rm orb}=18.462~\text{days}$, and $i=88.6^{\circ}$ and the ratios of the stellar radii and temperatures to be $R_2/R_1 = 0.71$ and $T_{e,2}/T_{e,1} = 0.82$. Both stars are chromospherically active, allowing us to determine their rotational periods of $P_1=1.52$ days and $P_2=1.79$ days, respectively. A LBT/MODS spectrum shows that the primary is a late-G or early-K type dwarf. Fits to the SED show that the luminosities and temperatures of the two stars are $L_1 = 0.48~L_{\sun}$, $T_1= 5050~K$, $L_2 = 0.12~L_{\sun}$, and $T_{2} = 4190~K$. We conclude that ASASSN-V J192114.84+624950.8 consists of two chromospherically active, rotational variable stars in a highly elliptical eclipsing orbit.

preprint2021arXiv

The ASAS-SN Catalog of Variable Stars II: Uniform Classification of 412,000 Known Variables

The variable stars in the VSX catalog are derived from a multitude of inhomogeneous data sources and classification tools. This inhomogeneity complicates our understanding of variable star types, statistics, and properties, and it directly affects attempts to build training sets for current (and next) generation all-sky, time-domain surveys. We homogeneously analyze the ASAS-SN V-band light curves of ${\sim}412,000$ variables from the VSX catalog. The variables are classified using an updated random forest classifier with an $F_1$ score of 99.4\% and refinement criteria for individual classifications. We have derived periods for ${\sim}52,000$ variables in the VSX catalog that lack a period, and have reclassified ${\sim} 17,000$ sources into new broad variability groups with high confidence. We have also reclassified ${\sim} 94,000$ known variables with miscellaneous/generic classifications. The light curves, classifications, and a range of properties of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). We also include the V-band light curves for a set of ${\sim}4,000$ rare variables and transient sources, including cataclysmic variables, symbiotic binaries and flare stars.

preprint2021arXiv

The ASAS-SN Catalog of Variable Stars III: Variables in the Southern TESS Continuous Viewing Zone

The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) light curves for sources brighter than V$\lesssim17$ mag across the whole sky. The Transiting Exoplanet Survey Satellite (TESS) has started to produce high-quality light curves with a baseline of at least 27 days, eventually for most of the sky. The combination of ASAS-SN and TESS light curves probes both long and short term variability in great detail, especially towards the TESS continuous viewing zones (CVZ) at the ecliptic poles. We have produced ${\sim}1.3$ million V-band light curves covering a total of ${\sim}1000 \, \rm deg^2$ towards the southern TESS CVZ and have systematically searched these sources for variability. We have identified ${\sim} 11,700$ variables, including ${\sim} 7,000$ new discoveries. The light curves and characteristics of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). We also introduce an online resource to obtain pre-computed ASAS-SN V-band light curves (https://asas-sn.osu.edu/photometry) starting with the light curves of the ${\sim}1.3$ million sources studied in this work. This effort will be extended to provide ASAS-SN light curves for ${\sim}50\;$million sources over the entire sky.

preprint2021arXiv

The ASAS-SN Catalog of Variable Stars V: Variables in the Southern Hemisphere

The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) light curves for sources brighter than V$\lesssim17$ mag across the whole sky. As part of our effort to characterize the variability of all the stellar sources visible in ASAS-SN, we have produced ${\sim}30.1$ million V-band light curves for sources in the southern hemisphere using the APASS DR9 catalog as our input source list. We have systematically searched these sources for variability using a pipeline based on random forest classifiers. We have identified ${\sim} 220,000$ variables, including ${\sim} 88,300$ new discoveries. In particular, we have discovered ${\sim}48,000$ red pulsating variables, ${\sim}23,000$ eclipsing binaries, ${\sim}2,200$ $δ$-Scuti variables and ${\sim}10,200$ rotational variables. The light curves and characteristics of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). The pre-computed ASAS-SN V-band light curves for all the ${\sim}30.1$ million sources are available through the ASAS-SN photometry database (https://asas-sn.osu.edu/photometry). This effort will be extended to provide ASAS-SN light curves for sources in the northern hemisphere and for V$\lesssim17$ mag sources across the whole sky that are not included in APASS DR9.

preprint2021arXiv

The ASAS-SN Catalog of Variable Stars VI: An All-Sky Sample of $δ$ Scuti Stars

We characterize an all-sky catalog of ${\sim} 8,400$ $δ$ Scuti variables in ASAS-SN, which includes ${\sim} 3,300$ new discoveries. Using distances from \textit{Gaia} DR2, we derive period-luminosity relationships (PLRs) for both the fundamental mode and overtone pulsators in the $W_{JK}$, $V$, Gaia DR2 $G$, $J$, $H$, $K_s$ and $W_1$ bands. We find that the overtone pulsators have a dominant overtone mode, with many sources pulsating in the second overtone or higher order modes. The fundamental mode pulsators have metallicity dependent periods, with $\log_{10} (\rm P){\sim}-1.1$ for $\rm [Fe/H]<-0.3$ and $\log_{10} (\rm P){\sim}-0.9$ for $\rm [Fe/H]>0$ which leads to a period-dependant scale height. Stars with $\rm P>0.100\,d$ are predominantly located close to the Galactic disk ($\rm |Z|<0.5\,kpc$). The median period at a scale height of $\rm Z{\sim}0\, kpc$ also increases with the Galactocentric radius $\rm R$, from $\log_{10} (\rm P){\sim}-0.94$ for sources with $\rm R>9\, kpc$ to $\log_{10} (\rm P){\sim}-0.85$ for sources with $\rm R<7\, kpc$, which is indicative of a radial metallicity gradient. To illustrate potential applications of this all-sky catalog, we obtained 30 min cadence, image subtraction TESS light curves for a sample of 10 fundamental mode and 10 overtone $δ$ Scuti stars discovered by ASAS-SN. From this sample, we identified two new $δ$ Scuti eclipsing binaries, ASASSN-V J071855.62$-$434247.3 and ASASSN-V J170344.20$-$615941.2 with short orbital periods of $\rm P_{orb}=2.6096$ d and $\rm P_{orb}=2.5347$ d respectively.

preprint2021arXiv

The ASAS-SN Catalog of Variable Stars VII: Contact Binaries are Different Above and Below the Kraft Break

We characterize ${\sim} 71,200$ W UMa type (EW) contact binaries, including ${\sim} 12,600$ new discoveries, using ASAS-SN $V$-band all-sky light curves along with archival data from Gaia, 2MASS, AllWISE, LAMOST, GALAH, RAVE, and APOGEE. There is a clean break in the EW period-luminosity relation at $\rm \log (\rm P/d){\simeq}-0.30$, separating the longer period early-type EW binaries from the shorter period, late-type systems. The two populations are even more cleanly separated in the space of period and effective temperature, by $\rm T_{eff}=6710\,K-1760\,K\,\log(P/0.5\,d)$. Early-type and late-type EW binaries follow opposite trends in $\rm T_{eff}$ with orbital period. For longer periods, early-type EW binaries are cooler, while late-type systems are hotter. We derive period-luminosity relationships (PLRs) in the $W_{JK}$, $V$, Gaia DR2 $G$, $J$, $H$, $K_s$ and $W_1$ bands for the late-type and early-type EW binaries separated both by period and effective temperature, and by period alone. The dichotomy of contact binaries is almost certainly related to the Kraft break and the related changes in envelope structure, winds and angular momentum loss.

preprint2021arXiv

The Progenitor of the Vela Pulsar

With Gaia parallaxes it is possible to study the stellar populations associated with individual Galactic supernova remnants (SNR) to estimate the mass of the exploding star. Here we analyze the luminous stars near the Vela pulsar and SNR to find that its progenitor was probably (>90%) low mass (8.1-10.3Msun). The presence of the O star gamma2 Vel a little over 100 pc from Vela is the primary ambiguity, as including it in the analysis volume significantly increases the probability (to 5%) of higher mass (>20Msun) progenitors. However, to be a high mass star associated with gamma2 Vel's star cluster at birth, the progenitor would have to a runaway star from an unbound binary with an unusually high velocity. The primary impediment to analyzing large numbers of Galactic SNRs in this manner is the lack of accurate distances. This can likely be solved by searching for absorption lines from the SNR in stars as a function of distance, a method which yielded a distance to Vela in agreement with the direct pulsar parallax. If Vela was a 10Msun supernova in an external galaxy, the 50 pc search region used in extragalactic studies would contain only ~10% of the stars formed in a 50~pc region around the progenitor at birth and ~90% of the stars in the search region would have been elsewhere.

preprint2020arXiv

A mildly relativistic outflow from the energetic, fast-rising blue optical transient CSS161010 in a dwarf galaxy

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_ν\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $Γβc\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($Γβc\sim0.1c$) and closer to ZTF18abvkwla ($Γβc\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

preprint2020arXiv

ASASSN-18tb: A Most Unusual Type Ia Supernova Observed by TESS and SALT

We present photometric and spectroscopic observations of the unusual Type Ia supernova ASASSN-18tb, including a series of SALT spectra obtained over the course of nearly six months and the first observations of a supernova by the Transiting Exoplanet Survey Satellite (TESS). We confirm a previous observation by Kollmeier et al. (2019) showing that ASASSN-18tb is the first relatively normal Type Ia supernova to exhibit clear broad ($\sim1000$ km s$^{-1}$) H$α$ emission in its nebular phase spectra. We find that this event is best explained as a sub-Chandrasekhar mass explosion with $M_{Ni} \approx 0.3\; \rm{M}_\odot$. Despite the strong H$α$ signature at late times, we find that the early rise of the supernova shows no evidence for deviations from a single-component power-law and is best fit with a moderately shallow power-law of index $1.69\pm0.04$. We find that the H$α$ luminosity remains approximately constant after its initial detection at phase +37 d, and that the H$α$ velocity evolution does not trace that of the Fe~III$~\lambda4660$ emission. These suggest that the H$α$ emission arises from circumstellar medium (CSM) rather than swept up material from a non-degenerate companion. However, ASASSN-18tb is strikingly different from other known CSM-interacting Type Ia supernovae in a number of significant ways. Those objects typically show an H$α$ luminosity two orders of magnitude higher than what is seen in ASASSN-18tb, pushing them away from the empirical light-curve relations that define "normal" Type Ia supernovae. Conversely, ASASSN-18tb exhibits a fairly typical light curve and luminosity for an underluminous or transitional SN Ia, with $M_R \approx -18.1$ mag. Moreover, ASASSN-18tb is the only SN Ia showing H$α$ from CSM interaction to be discovered in an early-type galaxy.

preprint2020arXiv

Beyond Gaia: Asteroseismic Distances of M giants using Ground-Based Transient Surveys

Evolved stars near the tip of the red giant branch (TRGB) show solar-like oscillations with periods spanning hours to months and amplitudes ranging from $\sim$1 mmag to $\sim$100 mmag. The systematic detection of the resulting photometric variations with ground-based telescopes would enable the application of asteroseismology to a much larger and more distant sample of stars than is currently accessible with space-based telescopes such as \textit{Kepler} or the ongoing Transiting Exoplanet Survey Satellite (\textit{TESS}) mission. We present an asteroseismic analysis of 493 M giants using data from two ground-based surveys: the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the All-Sky Automated Survey for Supernovae (ASAS-SN). By comparing the extracted frequencies with constraints from \textit{Kepler}, the Sloan Digital Sky Survey Apache Point Observatory Galaxy Evolution Experiment (APOGEE), and Gaia we demonstrate that ground-based transient surveys allow accurate distance measurements to oscillating M giants with a precision of $\sim$15$\%$. Using stellar population synthesis models we predict that ATLAS and ASAS-SN can provide asteroseismic distances to $\sim$2$\times$10$^{6}$ galactic M giants out to typical distances of $20-50 \; \rm{kpc}$, vastly improving the reach of Gaia and providing critical constraints for Galactic archaeology and galactic dynamics.

preprint2020arXiv

Deep Late-Time Observations of the Supernova Impostors SN 1954J and SN 1961V

SN 1954J in NGC 2403 and SN 1961V in NGC 1058 were two luminous transients whose definitive classification as either non-terminal eruptions or supernovae remains elusive. A critical question is whether a surviving star can be significantly obscured by dust formed from material ejected during the transient. We use three lines of argument to show that the candidate surviving stars are not significantly optically extincted ($τ\lesssim 1$) by dust formed in the transients. First, we use SED fits to new HST optical and near-IR photometry. Second, neither source is becoming brighter as required by absorption from an expanding shell of ejected material. Third, the ejecta masses implied by the H$α$ luminosities are too low to produce significant dust absorption. The latter two arguments hold independent of the dust properties. The H$α$ fluxes should also be declining with time as $t^{-3}$, and this seems not to be observed. As a result, it seems unlikely that recently formed dust can be responsible for the present faintness of the sources compared to their progenitors, although this can be verified with \textit{JWST}. This leaves three possibilities: 1) the survivors were misidentified; 2) they are intrinsically less luminous; 3) SN 1954J and SN 1961V were true supernovae.

preprint2020arXiv

Early Time Light Curves of Type Ia Supernovae Observed with TESS

We present early time light curves of Type Ia supernovae observed in the first six sectors of TESS data. Ten of these supernovae were discovered by ASAS-SN, seven by ATLAS, six by ZTF, and one by \textit{Gaia}. For nine SNe with sufficient dynamic range ($>$3.0 mag from detection to peak), we fit power law models and search for signatures of companion stars. We find a diversity of early time light curve shapes, although most of our sources are consistent with fireball models where the flux increases $\propto t^2$. Three SN display a flatter rise with flux $\propto t$. We do not find any evidence for additional structure such as multiple power law components in the early rising light curves. For assumptions about the SN properties and the observer viewing angle, and further assuming that companion stars would be in Roche-lobe overflow, we place limits on the radii of companions for six SNe with complete coverage of the early time light curves. The upper limits are $\lesssim$\,32 R$_\odot$ for these six supernovae, $\lesssim$\,20 R$_\odot$ for five of these six, and $\lesssim$\,4 R$_\odot$ for two of these six. The small sample size does not constrain occurrence rates of single degenerate Type Ia SN progenitors, but we expect that TESS observed enough SNe in its primary mission (26 sectors) to inform this measurement. We also show that TESS is capable of detecting emission from a 1 \rsun\ companion for a Type Ia SN within 50 Mpc, and may do so after about six years.

preprint2020arXiv

On Reverberation Mapping Lag Uncertainties

We broadly explore the effects of systematic errors on reverberation mapping lag uncertainty estimates from {\tt JAVELIN} and the interpolated cross-correlation function (ICCF) method. We focus on simulated lightcurves from random realizations of the lightcurves of five intensively monitored AGNs. Both methods generally work well even in the presence of systematic errors, although {\tt JAVELIN} generally provides better error estimates. Poorly estimated lightcurve uncertainties have less effect on the ICCF method because, unlike {\tt JAVELIN}, it does not explicitly assume Gaussian statistics. Neither method is sensitive to changes in the stochastic process driving the continuum or the transfer function relating the line lightcurve to the continuum. The only systematic error we considered that causes significant problems is if the line lightcurve is not a smoothed and shifted version of the continuum lightcurve but instead contains some additional sources of variability.

preprint2020arXiv

On The Red Supergiant Problem

We examine the problem of estimating the mass range corresponding to the observed red supergiant (RSG) progenitors of Type IIP supernovae. Using Monte Carlo simulations designed to reproduce the properties of the observations, we find that the approach of Davies & Beasor (2018) significantly overestimates the maximum mass, yielding an upper limit of Mh/Msun=20.5+/-2.6 for an input population with Mh/Msun=18. Our preferred Bayesian approach does better, with Mh/Msun=18.6+/-2.1 for the same input populations, but also tends to overestimate Mh. For the actual progenitor sample and a Salpeter initial mass function we find Mh/Msun=19.01-2.04+4.04 for the Eldridge et al. (2004) mass-luminosity relation used by Smartt et al. (2009) and Davies & Beasor (2018), and Mh/Msun=21.28_-2.28+4.52 for the Sukhbold et al. (2018) mass-luminosity relation. Based on the Monte Carlo simulations, we estimate that these are overestimated by 3.3+/-0.8Mh. The red supergiant problem remains.

preprint2020arXiv

Over-constrained Models of Time Delay Lenses Redux: How the Angular Tail Wags the Radial Dog

The two properties of the radial mass distribution of a gravitational lens that are well-constrained by Einstein rings are the Einstein radius R_E and xi2 = R_E alpha''(R_E)/(1-kappa_E), where alpha''(R_E) and kappa_E are the second derivative of the deflection profile and the convergence at R_E. However, if there is a tight mathematical relationship between the radial mass profile and the angular structure, as is true of ellipsoids, an Einstein ring can appear to strongly distinguish radial mass distributions with the same xi2. This problem is beautifully illustrated by the ellipsoidal models in Millon et al. (2019). When using Einstein rings to constrain the radial mass distribution, the angular structure of the models must contain all the degrees of freedom expected in nature (e.g., external shear, different ellipticities for the stars and the dark matter, modest deviations from elliptical structure, modest twists of the axes, modest ellipticity gradients, etc.) that work to decouple the radial and angular structure of the gravity. Models of Einstein rings with too few angular degrees of freedom will lead to strongly biased likelihood distinctions between radial mass distributions and very precise but inaccurate estimates of H0 based on gravitational lens time delays.

preprint2020arXiv

Quasar Accretion Disk Sizes from Continuum Reverberation Mapping in the DES Standard Star Fields

Measurements of the physical properties of accretion disks in active galactic nuclei are important for better understanding the growth and evolution of supermassive black holes. We present the accretion disk sizes of 22 quasars from continuum reverberation mapping with data from the Dark Energy Survey (DES) standard star fields and the supernova C fields. We construct continuum lightcurves with the \textit{griz} photometry that span five seasons of DES observations. These data sample the time variability of the quasars with a cadence as short as one day, which corresponds to a rest frame cadence that is a factor of a few higher than most previous work. We derive time lags between bands with both JAVELIN and the interpolated cross-correlation function method, and fit for accretion disk sizes using the JAVELIN Thin Disk model. These new measurements include disks around black holes with masses as small as $\sim10^7$ $M_{\odot}$, which have equivalent sizes at 2500Å\, as small as $\sim 0.1$ light days in the rest frame. We find that most objects have accretion disk sizes consistent with the prediction of the standard thin disk model when we take disk variability into account. We have also simulated the expected yield of accretion disk measurements under various observational scenarios for the Large Synoptic Survey Telescope Deep Drilling Fields. We find that the number of disk measurements would increase significantly if the default cadence is changed from three days to two days or one day.

preprint2020arXiv

Signatures of Bimodality in Nebular Phase Type Ia Supernova Spectra

One observational prediction for Type Ia supernova (SNe Ia) explosions produced through white dwarf-white dwarf collisions is the presence of bimodal velocity distributions for the $^{56}$Ni decay products, although this signature can also be produced by an off-center ignition in a delayed detonation explosion. These bimodal velocity distributions can manifest as double-peaked or flat-topped spectral features in late-time spectroscopic observations for favorable viewing angles. We present nebular-phase spectroscopic observations of 17 SNe Ia obtained with the Large Binocular Telescope (LBT). Combining these observations with an extensive search of publicly available archival data, we collect a total sample of 48 SNe Ia and classify them based on whether they show compelling evidence for bimodal velocity profiles in three features associated with $^{56}$Ni decay products: the [Fe II] and [Fe III] feature at $\sim5300$ Å, the [Co III] $\lambda5891$ feature, and the [Co III] and [Fe II] feature at $\sim6600$ Å. We identify 9 bimodal SNe in our sample, and we find that these SNe have average peak $M_V$ about 0.3 mag fainter than those which do not. This is consistent with theoretical predictions for explosions created by nearly head-on collisions of white dwarfs due to viewing angle effects and $^{56}$Ni yields.

preprint2020arXiv

Space Telescope and Optical Reverberation Mapping Project. IX. Velocity-Delay Maps for Broad Emission Lines in NGC 5548

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

preprint2020arXiv

The ASAS-SN Catalog of Variable Stars VIII: "Dipper" Stars in the Lupus Star-Forming Region

Some young stellar objects such as T Tauri-like "dipper" stars vary due to transient partial occultation by circumstellar dust, and observations of this phenomenon inform us of conditions in the planet-forming zones close to these stars. Although many dipper stars have been identified with space missions such as $Kepler$/$K2$, ground-based telescopes offer longer term and multi-wavelength perspectives. We identified 11 dipper stars in the Lupus star forming region in data from the All-Sky Automated Survey for SuperNovae (ASAS-SN), and further characterized these using observations by the Las Cumbres Global Observatory Telescope (LCOGT) and the Transiting Exoplanet Survey Satellite $TESS$, as well as archival data from other missions. Dipper stars were identified from a catalog of nearby young stars and selected based on the statistical significance, asymmetry, and quasi-periodicity or aperiodicity of variability in their ASAS-SN light curves. All 11 stars lie above or red-ward of the zero-age main sequence and have infrared excesses indicating the presence of full circumstellar disks. We obtain reddening-extinction relations for the variability of 7 stars using our combined ASAS-SN-$TESS$ and LCOGT photometry. In all cases the slopes are below the ISM value, suggesting larger grains, and we find a tentative relation between the slope (grain size) and the $\text{K}_\text{s}-[22\:μ\text{m}]$ infrared color regarded as a proxy for disk evolutionary state.

preprint2020arXiv

The Highly Luminous Type Ibn Supernova ASASSN-14ms

We present photometric and spectroscopic follow-up observations of the highly luminous Type Ibn supernova ASASSN-14ms, which was discovered on UT 2014-12-26.61 at $m_V \sim 16.5$. With a peak absolute $V$-band magnitude brighter than $-20.5$, a peak bolometric luminosity of $1.7 \times 10^{44}$ ergs s$^{-1}$, and a total radiated energy of $2.1 \times 10^{50}$ ergs, ASASSN-14ms is one of the most luminous Type Ibn supernovae yet discovered. In simple models, the most likely power source for this event is a combination of the radioactive decay of $^{56}$Ni and $^{56}$Co at late times and the interaction of supernova ejecta with the progenitor's circumstellar medium at early times, although we cannot rule out the possibility of a magnetar-powered light curve. The presence of a dense circumstellar medium is indicated by the intermediate-width He I features in the spectra. The faint ($m_g \sim 21.6$) host galaxy SDSS J130408.52+521846.4 has an oxygen abundance below $12+\log(O/H) \lesssim 8.3$, a stellar mass of $M_* \sim 2.6 \times 10^8 M_{\odot}$, and a star formation rate of $\textrm{SFR} \sim 0.02$ $M_{\odot}$ yr$^{-1}$.

preprint2020arXiv

The Most Rapidly Declining Type I Supernova 2019bkc/ATLAS19dqr

We report observations of the hydrogen-deficient supernova (SN) 2019bkc/ATLAS19dqr. With B- and r-band decline between peak and 10 days post peak of Delta m_10(B)=5.24+/-0.07 mag and Delta m_10(r)=3.85+/-0.10$ mag, respectively, SN 2019bkc is the most rapidly declining SN I discovered so far. While its closest matches are the rapidly declining SN 2005ek and SN 2010X, the light curves and spectra of SN 2019bkc show some unprecedented characteristics. SN 2019bkc appears "hostless," with no identifiable host galaxy near its location, although it may be associated with the galaxy cluster MKW1 at z = 0.02. We evaluate a number of existing models of fast-evolving SNe, and we find that none of them can satisfactorily explain all aspects of SN 2019bkc observations.

preprint2020arXiv

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

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

preprint2019arXiv

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived from multidimensional numerical simulations and phenomenological extreme emission models. For neutrino-driven explosions the distance at which we reach 50% detection efficiency is approaching 5 kpc, and for magnetorotationally-driven explosions is up to 54 kpc. However, waveforms for extreme emission models are detectable up to 28 Mpc. For the first time, the gravitational-wave data enabled us to exclude part of the parameter spaces of two extreme emission models with confidence up to 83%, limited by coincident data coverage. Besides, using ad hoc harmonic signals windowed with Gaussian envelopes we constrained the gravitational-wave energy emitted during core-collapse at the levels of $4.27\times 10^{-4}\,M_\odot c^2$ and $1.28\times 10^{-1}\,M_\odot c^2$ for emissions at 235 Hz and 1304 Hz respectively. These constraints are two orders of magnitude more stringent than previously derived in the corresponding analysis using initial LIGO, initial Virgo and GEO 600 data.

preprint2019arXiv

Nebular Spectra of 111 Type Ia Supernovae Disfavor Single Degenerate Progenitors

We place statistical constraints on Type Ia supernova (SN Ia) progenitors using 227 nebular phase spectra of 111 SNe Ia. We find no evidence of stripped companion emission in any of the nebular phase spectra. Upper limits are placed on the amount of mass that could go undetected in each spectrum using recent hydrodynamic simulations. With these null detections, we place an observational $3σ$ upper limit on the fraction of SNe Ia that are produced through the classical H-rich non-degenerate companion scenario of < 5.5%. Additionally, we set a tentative $3σ$ upper limit on He star progenitor scenarios of < 6.4%, although further theoretical modelling is required. These limits refer to our most representative sample including normal, 91bg-like, 91T-like, and "Super Chandrasekhar" \sne but excluding SNe Iax and SNe Ia-CSM. As part of our analysis, we also derive a Nebular Phase Phillips Relation, which approximates the brightness of a SN Ia from $150-500$~days after maximum using the peak magnitude and decline rate parameter $Δm_{15} (B)$.

preprint2019arXiv

Over-constrained Gravitational Lens Models and the Hubble Constant

It is well known that measurements of H0 from gravitational lens time delays scale as H0~1-k_E where k_E is the mean convergence at the Einstein radius R_E but that all available lens data other than the delays provide no direct constraints on k_E. The properties of the radial mass distribution constrained by lens data are R_E and the dimensionless quantity x=R_E a''(R_E)/(1-k_E)$ where a''(R_E) is the second derivative of the deflection profile at R_E. Lens models with too few degrees of freedom, like power law models with densities ~r^(-n), have a one-to-one correspondence between x and k_E (for a power law model, x=2(n-2) and k_E=(3-n)/2=(2-x)/4). This means that highly constrained lens models with few parameters quickly lead to very precise but inaccurate estimates of k_E and hence H0. Based on experiments with a broad range of plausible dark matter halo models, it is unlikely that any current estimates of H0 from gravitational lens time delays are more accurate than ~10%, regardless of the reported precision.

preprint2019arXiv

The New EXor Outburst of ESO-H-alpha~99 observed by Gaia ATLAS and TESS

We report photometry and spectroscopy of the outburst of the young stellar object ESO-Halpha 99. The outburst was first noticed in Gaia alert Gaia18dvc and later by ATLAS. We have established the outburst light curve with archival ATLAS ``Orange'' filter photometry, Gaia data, new V-band photometry, and J, H, and K_s photometry from IRIS and UKIRT. The brightness has fluctuated several times near the light curve maximum. The TESS satellite observed ESO-Halpha 99 with high cadence during one of these minor minima and found brightness fluctuations on timescales of days and hours. Imaging with UKIRT shows the outline of an outflow cavity, and we find one knot of H_2 1-0 S(1) emission, now named MHO 1520, on the symmetry axis of this nebula, indicating recent collimated outflow activity from ESO-Halpha 99. Its pre-outburst SED shows a flat FIR spectrum, confirming its early evolutionary state and its similarity to other deeply embedded objects in the broader EXor class. The pre-outburst luminosity is 34 L_sun, a much higher luminosity than typical EXors, indicating that ESO-Halpha 99 may be a star of intermediate mass. Infrared and optical spectroscopy show a rich emission line spectrum, including HI lines, strong red CaII emission, as well as infrared CO bandhead emission, all characteristic EXors in the broadest sense. Comparison of the present spectra with an optical spectrum obtained in 1993, presumably in the quiescent state of the object, shows that during the present outburst the continuum component of the spectrum has increased notably more than the emission lines. The Halpha equivalent width during the outburst is down to one half of its 1993 level, and shock-excited emission lines are much less prominent.

preprint2017arXiv

Where Is the Flux Going? The Long-Term Photometric Variability of Boyajian's Star

We present ~800 days of photometric monitoring of Boyajian's Star (KIC 8462852) from the All-Sky Automated Survey for Supernovae (ASAS-SN) and ~4000 days of monitoring from the All Sky Automated Survey (ASAS). We show that from 2015 to the present the brightness of Boyajian's Star has steadily decreased at a rate of 6.3 +/- 1.4 mmag yr^-1, such that the star is now 1.5% fainter than it was in February 2015. Moreover, the longer time baseline afforded by ASAS suggests that Boyajian's Star has also undergone two brightening episodes in the past 11 years, rather than only exhibiting a monotonic decline. We analyze a sample of ~1000 comparison stars of similar brightness located in the same ASAS-SN field and demonstrate that the recent fading is significant at >99.4% confidence. The 2015-2017 dimming rate is consistent with that measured with Kepler data for the time period from 2009 to 2013. This long-term variability is difficult to explain with any of the physical models for the star's behavior proposed to date.

preprint2016arXiv

Application of Stochastic Modeling to Analysis of Photometric Reverberation Mapping Data

We use both simulated and real quasar light curves to explore modeling photometric reverberation-mapping (RM) data as a stochastic process. We do this using modifications to our previously developed RM method based on modeling quasar variability as a damped random walk. We consider the feasibility of one and two-band photometric RM and compare the results with those fr om spectroscopic RM. We find that our method for two-band photometric RM can be competitive with spectroscopic RM only for strong (large equivalent width) lines like Halpha and Hbeta, and that the one-band method is also feasible, but requires very high-precision photometry. We fail to robustly detect Halpha lags in single-band quasar light curves from OGLE-III and IV despite the outstanding cadence and time span of the data, on account of photometric uncertainties in the range 0.02-0.04 mag. Simulations suggest that success could be achieved if the photometric uncertainties were of order 0.01 mag. Single-band RM for all lines and two-band RM for lower equivalent width lines are likely only feasible for statistical estimate s of mean lags for large samples of AGN of similar properties (e.g., luminosity) rather than for individual quasars. Our approach is directly applicable to the time-domain programs within ongoing and future wide-field imaging surveys, and could provide robust lag measurements for an unprecedented number of systems.

preprint2016arXiv

ASASSN-15lh: A Highly Super-Luminous Supernova

We report the discovery of ASASSN-15lh (SN 2015L), which we interpret as the most luminous supernova yet found. At redshift z = 0.2326, ASASSN-15lh reached an absolute magnitude of M_{u,AB} = -23.5+/-0.1 and bolometric luminosity L_bol = (2.2+/-0.2)x 10^45 ergs s^-1, which is more than twice as luminous as any previously known supernova. It has several major features characteristic of the hydrogen-poor super-luminous supernovae (SLSNe-I), whose energy sources and progenitors are currently poorly understood. In contrast to most previously known SLSNe-I that reside in star-forming dwarf galaxies, ASASSN-15lh appears to be hosted by a luminous galaxy (M_K ~ -25.5) with little star formation. In the 4 months since first detection, ASASSN-15lh radiated (1.1+/- 0.2)x10^52 ergs, challenging the magnetar model for its engine.

preprint2016arXiv

ASASSN-15oi: A Rapidly Evolving, Luminous Tidal Disruption Event at 216 Mpc

We present ground-based and Swift photometric and spectroscopic observations of the tidal disruption event (TDE) ASASSN-15oi, discovered at the center of 2MASX J20390918-3045201 ($d\simeq216$ Mpc) by the All-Sky Automated Survey for SuperNovae (ASAS-SN). The source peaked at a bolometric luminosity of $L\simeq1.3\times10^{44}$ ergs s$^{-1}$ and radiated a total energy of $E\simeq6.6\times10^{50}$ ergs over the first $\sim3.5$ months of observations. The early optical/UV emission of the source can be fit by a blackbody with temperature increasing from $T\sim2\times10^4$ K to $T\sim4\times10^4$ K while the luminosity declines from $L\simeq1.3\times10^{44}$ ergs s$^{-1}$ to $L\simeq2.3\times10^{43}$ ergs s$^{-1}$, requiring the photosphere to be shrinking rapidly. The optical/UV luminosity decline during this period is most consistent with an exponential decline, $L\propto e^{-(t-t_0)/τ}$, with $τ\simeq46.5$ days for $t_0\simeq57241.6$ (MJD), while a power-law decline of $L\propto (t-t_0)^{-α}$ with $t_0\simeq57212.3$ and $α=1.62$ provides a moderately worse fit. ASASSN-15oi also exhibits roughly constant soft X-ray emission that is significantly weaker than the optical/UV emission. Spectra of the source show broad helium emission lines and strong blue continuum emission in early epochs, although these features fade rapidly and are not present $\sim3$ months after discovery. The early spectroscopic features and color evolution of ASASSN-15oi are consistent with a TDE, but the rapid spectral evolution is unique among optically-selected TDEs.

preprint2016arXiv

ASASSN-16ae: A Powerful White-Light Flare on an Early-L Dwarf

We report the discovery and classification of SDSS~J053341.43+001434.1 (SDSS0533), an early-L dwarf first discovered during a powerful $ΔV < -11$ magnitude flare observed as part of the ASAS-SN survey. Optical and infrared spectroscopy indicate a spectral type of L0 with strong H$α$ emission and a blue NIR spectral slope. Combining the photometric distance, proper motion, and radial velocity of SDSS0533 yields three-dimensional velocities of $(U,V,W)=(14\pm13,-35\pm14,-94\pm22)$~km~s$^{-1}$, indicating that it is most likely part of the thick disk population and probably old. The three detections of SDSS0533 obtained during the flare are consistent with a total $V$-band flare energy of at least $4.9\times10^{33}$~ergs (corresponding to a total thermal energy of at least $E_{\rm tot}>3.7\times10^{34}$~erg), placing it among the strongest detected M dwarf flares. The presence of this powerful flare on an old L0 dwarf may indicate that stellar-type magnetic activity persists down to the end of the main sequence and on older ML transition dwarfs.

preprint2016arXiv

Discovery and Observations of the Unusually Luminous Type-Defying II-P/II-L Supernova ASASSN-13co

We present photometric and spectroscopic observations of ASASSN-13co, an unusually luminous Type II supernova and the first core-collapse supernova discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN). First detection of the supernova was on UT 2013 August 29 and the data presented span roughly 3.5 months after discovery. We use the recently developed model from Pejcha & Prieto (2015) to model the multi-band light curves of ASASSN-13co and derive the bolometric luminosity curve. We compare ASASSN-13co to other Type II supernovae to show that it was unusually luminous for a Type II supernova and that it exhibited an atypical light curve shape that does not cleanly match that of either a standard Type II-L or Type II-P supernova.

preprint2016arXiv

Hello Darkness My Old Friend: The Fading of the Nearby TDE ASASSN-14ae

We present late-time optical spectroscopy taken with the Large Binocular Telescope's Multi-Object Double Spectrograph, an improved ASAS-SN pre-discovery non-detection, and late-time SWIFT observations of the nearby ($d=193$ Mpc, $z=0.0436$) tidal disruption event (TDE) ASASSN-14ae. Our observations span from $\sim$20 days before to $\sim$750 days after discovery. The proximity of ASASSN-14ae allows us to study the optical evolution of the flare and the transition to a host dominated state with exceptionally high precision. We measure very weak H$α$ emission 300 days after discovery ($L_{\rm Hα} \simeq 4\times 10^{39}$ ergs s$^{-1}$) and the most stringent upper limit to date on the H$α$ luminosity $\sim$750 days after discovery ($L_{\rm Hα} \lesssim 10^{39}$ ergs s$^{-1}$), suggesting that the optical emission arising from a TDE can vanish on a timescale as short as 1 year. Our results have important implications for both spectroscopic detection of TDE candidates at late times, as well as the nature of TDE host galaxies themselves.

preprint2016arXiv

MUSE Reveals a Recent Merger in the Post-starburst Host Galaxy of the TDE ASASSN-14li

We present MUSE integral field spectroscopic observations of the host galaxy (PGC 043234) of one of the closest ($z=0.0206$, $D\simeq 90$ Mpc) and best-studied tidal disruption events (TDE), ASASSN-14li. The MUSE integral field data reveal asymmetric and filamentary structures that extend up to $\gtrsim 10$ kpc from the post-starburst host galaxy of ASASSN-14li. The structures are traced only through the strong nebular [O III] $λ$5007, [N II] $λ$6584, and H$α$ emission lines. The total off nuclear [O III] $λ$5007 luminosity is luminosity is $4.7\times 10^{39}$ erg s$^{-1}$ and the ionized H mass is $\rm \sim 10^4(500/n_e)\,M_{\odot}$. Based on the BPT diagram, the nebular emission can be driven by either AGN photoionization or shock excitation, with AGN photoionization favored given the narrow intrinsic line widths. The emission line ratios and spatial distribution strongly resemble ionization nebulae around fading AGNs such as IC 2497 (Hanny's Voorwerp) and ionization "cones" around Seyfert 2 nuclei. The morphology of the emission line filaments strongly suggest that PGC 043234 is a recent merger, which likely triggered a strong starburst and AGN activity leading to the post-starburst spectral signatures and the extended nebular emission line features we see today. We briefly discuss the implications of these observations in the context of the strongly enhanced TDE rates observed in post-starburst galaxies and their connection to enhanced theoretical TDE rates produced by supermassive black-hole binaries.

preprint2016arXiv

Observations of the Lensed Quasar Q2237+0305 with CanariCam at GTC

We present new mid-IR observations of the quadruply lensed quasar Q2237+0305 taken with CanariCam on the Gran Telescopio Canarias. Mid-IR emission by hot dust, unlike the optical and near-IR emission from the accretion disk, is unaffected by the interstellar medium (extinction/scattering) or stellar microlensing. We compare these "true" ratios to the (stellar) microlensed flux ratios observed in the optical/near-IR to constrain the structure of the quasar accretion disk. We find a half-light radius of $R_{1/2}=3.4_{-2.1}^{+5.3}\sqrt{\langle M \rangle/0.3\,\rm{M_{\odot}}}$ light-days at $λ_{rest}=1736$ Å, and an exponent for the temperature profile $R \propto λ^{p}$ of $p=0.79\pm0.55$, where $p=4/3$ for a standard thin-disk model. If we assume that the differences in the mid-IR flux ratios measured over the years are due to microlensing variability, we find a lower limit for the size of the mid-IR-emitting region of $R_{1/2} \gtrsim 200\,\sqrt{\langle M \rangle/0.3\,\rm{M_{\odot}}}$ light-days. We also test for the presence of substructure/satellites by comparing the observed mid-IR flux ratios with those predicted from smooth lens models. We can explain the differences if the surface density fraction in satellites near the lensed images is $α= 0.033_{-0.019}^{+0.046}$ for a singular isothermal ellipsoid plus external shear mass model or $α= 0.013_{-0.008}^{+0.019}$ for a mass model combining ellipsoidal NFW and de Vaucouleurs profiles in an external shear.

preprint2016arXiv

Six Months of Multi-Wavelength Follow-up of the Tidal Disruption Candidate ASASSN-14li and Implied TDE Rates from ASAS-SN

We present ground-based and Swift photometric and spectroscopic observations of the candidate tidal disruption event (TDE) ASASSN-14li, found at the center of PGC 043234 ($d\simeq90$ Mpc) by the All-Sky Automated Survey for SuperNovae (ASAS-SN). The source had a peak bolometric luminosity of $L\simeq10^{44}$ ergs s$^{-1}$ and a total integrated energy of $E\simeq7\times10^{50}$ ergs radiated over the $\sim6$ months of observations presented. The UV/optical emission of the source is well-fit by a blackbody with roughly constant temperature of $T\sim35,000$ K, while the luminosity declines by roughly a factor of 16 over this time. The optical/UV luminosity decline is broadly consistent with an exponential decline, $L\propto e^{-t/t_0}$, with $t_0\simeq60$ days. ASASSN-14li also exhibits soft X-ray emission comparable in luminosity to the optical and UV emission but declining at a slower rate, and the X-ray emission now dominates. Spectra of the source show broad Balmer and helium lines in emission as well as strong blue continuum emission at all epochs. We use the discoveries of ASASSN-14li and ASASSN-14ae to estimate the TDE rate implied by ASAS-SN, finding an average rate of $r \simeq 4.1 \times 10^{-5}~{\rm yr}^{-1}$ per galaxy with a 90% confidence interval of $(2.2 - 17.0) \times 10^{-5}~{\rm yr}^{-1}$ per galaxy. ASAS-SN found roughly 1 TDE for every 70 Type Ia supernovae in 2014, a rate that is much higher than that of other surveys.

preprint2016arXiv

SN 2015bn: a detailed multi-wavelength view of a nearby superluminous supernova

We present observations of SN 2015bn (= PS15ae = CSS141223-113342+004332 = MLS150211-113342+004333), a Type I superluminous supernova (SLSN) at redshift $z=0.1136$. As well as being one of the closest SLSNe I yet discovered, it is intrinsically brighter ($M_U\approx-23.1$) and in a fainter galaxy ($M_B\approx-16.0$) than other SLSNe at $z\sim0.1$. We used this opportunity to collect the most extensive dataset for any SLSN I to date, including densely-sampled spectroscopy and photometry, from the UV to the NIR, spanning $-$50 to +250 days from optical maximum. SN 2015bn fades slowly, but exhibits surprising undulations in the light curve on a timescale of 30-50 days, especially in the UV. The spectrum shows extraordinarily slow evolution except for a rapid transformation between +7 and +20-30 days. No narrow emission lines from slow-moving material are observed at any phase. We derive physical properties including the bolometric luminosity, and find slow velocity evolution and non-monotonic temperature and radial evolution. A deep radio limit rules out a healthy off-axis gamma-ray burst, and places constraints on the pre-explosion mass loss. The data can be consistently explained by a $\gtrsim10\,{\rm M}_\odot$ stripped progenitor exploding with $\sim 10^{51}\,$erg kinetic energy, forming a magnetar with a spin-down timescale of $\sim20$ days (thus avoiding a gamma-ray burst) that reheats the ejecta and drives ionization fronts. The most likely alternative scenario -- interaction with $\sim20\,{\rm M}_\odot$ of dense, inhomogeneous circumstellar material -- can be tested with continuing radio follow-up.

preprint2016arXiv

Space Telescope and Optical Reverberation Mapping Project. III. Optical Continuum Emission and Broad-Band Time Delays in NGC 5548

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,Å to the $z$ band ($\sim\!9160$\,Å). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$λ1640$ and $λ4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $τ\propto λ^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

preprint2016arXiv

Space Telescope and Optical Reverberation Mapping Project. IV. Anomalous behavior of the broad ultraviolet emission lines in NGC 5548

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

preprint2016arXiv

Structure of the Accretion Disk in the Lensed Quasar Q2237+0305 from Multi-Epoch and Multi-Wavelength Narrow Band Photometry

We present estimates for the size and the logarithmic slope of the disk temperature profile of the lensed quasar Q2237+0305, independent of the component velocities. These estimates are based on six epochs of multi-wavelength narrowband images from the Nordic Optical Telescope. For each pair of lensed images and each photometric band, we determine the microlensing amplitude and chromaticity using pre-existing mid-IR photometry to define the baseline for no microlensing magnification. A statistical comparison of the combined microlensing data (6 epochs $\times$ 5 narrow bands $\times$ 6 image pairs) with simulations based on microlensing magnification maps gives Bayesian estimates for the half-light radius of $R_{1/2}=8.5^{+7.5}_{-4.0}\sqrt{ \langle M \rangle/0.3\, M_\odot}$ light-days, and $p=0.95\pm0.33$ for the exponent of the logarithmic temperature profile $T\propto R^{ -1/p}$. This size estimate is in good agreement with most recent studies. Other works based on the study of single microlensing events predict smaller sizes, but could be statistically biased by focusing on high-magnification events.

preprint2016arXiv

The ASAS-SN Bright Supernova Catalog $-$ I. 2013$-$2014

We present basic statistics for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) during its first year-and-a-half of operations, spanning 2013 and 2014. We also present the same information for all other bright ($m_V\leq17$), spectroscopically confirmed supernovae discovered from 2014 May 1 through the end of 2014, providing a comparison to the ASAS-SN sample starting from the point where ASAS-SN became operational in both hemispheres. In addition, we present collected redshifts and near-UV through IR magnitudes, where available, for all host galaxies of the bright supernovae in both samples. This work represents a comprehensive catalog of bright supernovae and their hosts from multiple professional and amateur sources, allowing for population studies that were not previously possible because the all-sky emphasis of ASAS-SN redresses many previously existing biases. In particular, ASAS-SN systematically finds bright supernovae closer to the centers of host galaxies than either other professional surveys or amateurs, a remarkable result given ASAS-SN's poorer angular resolution. This is the first of a series of yearly papers on bright supernovae and their hosts that will be released by the ASAS-SN team.

preprint2016arXiv

The Young and Bright Type Ia Supernova ASASSN-14lp: Discovery, Early-Time Observations, First-Light Time, Distance to NGC 4666, and Progenitor Constraints

On 2014 Dec. 9.61, the All-Sky Automated Survey for SuperNovae (ASAS-SN or "Assassin") discovered ASASSN-14lp just $\sim2$ days after first light using a global array of 14-cm diameter telescopes. ASASSN-14lp went on to become a bright supernova ($V = 11.94$ mag), second only to SN 2014J for the year. We present prediscovery photometry (with a detection less than a day after first light) and ultraviolet through near-infrared photometric and spectroscopic data covering the rise and fall of ASASSN-14lp for more than 100 days. We find that ASASSN-14lp had a broad light curve ($Δm_{15}(B) = 0.80 \pm 0.05$), a $B$-band maximum at $2457015.82 \pm 0.03$, a rise time of $16.94^{+ 0.11 }_{- 0.10 }$ days, and moderate host--galaxy extinction ($E(B-V)_{\textrm{host}} = 0.33 \pm 0.06$). Using ASASSN-14lp we derive a distance modulus for NGC 4666 of $μ= 30.8 \pm 0.2$ corresponding to a distance of $14.7 \pm 1.5$ Mpc. However, adding ASASSN-14lp to the calibrating sample of Type Ia supernovae still requires an independent distance to the host galaxy. Finally, using our early-time photometric and spectroscopic observations, we rule out red giant secondaries and, assuming a favorable viewing angle and explosion time, any non-degenerate companion larger than $0.34 R_{\textrm{sun}}$.

preprint2016arXiv

Tidal Disruption Event (TDE) Demographics

We survey the properties of stars destroyed in TDEs as a function of BH mass, stellar mass and evolutionary state, star formation history and redshift. For Mbh<10^7Msun, the typical TDE is due to a M*~0.3Msun M-dwarf, although the mass function is relatively flat for $M*<Msun. The contribution from older main sequence stars and sub-giants is small but not negligible. From Mbh~10^7.5-10^8.5Msun, the balance rapidly shifts to higher mass stars and a larger contribution from evolved stars, and is ultimately dominated by evolved stars at higher BH masses. The star formation history has little effect until the rates are dominated by evolved stars. TDE rates should decline very rapidly towards higher redshifts. The volumetric rate of TDEs is very high because the BH mass function diverges for low masses. However, any emission mechanism which is largely Eddington-limited for low BH masses suppresses this divergence in any observed sample and leads to TDE samples dominated by Mbh~10^6.0-10^7.5Msun BHs with roughly Eddington peak accretion rates. The typical fall back time is relatively long, with 16% having Tfb<10^(-1) years (37 days), and 84% having longer time scales. Many residual rate discrepancies can be explained if surveys are biased against TDEs with these longer Tfb, which seems very plausible if Tfb has any relation to the transient rise time. For almost any BH mass function, systematic searches for fainter, faster time scale TDEs in smaller galaxies, and longer time scale TDEs in more massive galaxies are likely to be rewarded.

preprint2015arXiv

Abundance Anomalies In Tidal Disruption Events

The ~10% of tidal disruption events (TDEs) due to stars more massive than the Sun should show abundance anomalies due to stellar evolution in helium, carbon and nitrogen, but not oxygen. Helium is always enhanced, but only by up to ~25% on average because it becomes inaccessible once it is sequestered in the high density core as the star leaves the main sequence. However, portions of the debris associated with the disrupted core of a main sequence star can be enhanced in helium by factors of 2-3 for debris at a common orbital period. These helium abundance variations may be a contributor to the observed diversity of hydrogen and helium line strengths in TDEs. A still more striking anomaly is the rapid enhancement of nitrogen and the depletion of carbon due to the CNO cycle -- stars more massive than the Sun quickly show an increase in their average N/C ratio by factors of 3-10. Because low mass stars evolve slowly and high mass stars are rare, TDEs showing high N/C will almost all be due to 1-2Msun stars disrupted on the main sequence. Like helium, portions of the debris will show still larger changes in C and N, and the anomalies decline as the star leaves the main sequence. The enhanced [N/C] abundance ratio of these TDEs provides the first natural explanation for the rare, nitrogen rich quasars and also explains the strong nitrogen emission seen in ultraviolet spectra of ASASSN-14li.

preprint2015arXiv

Almost Gone: SN~2008S and NGC~300 2008OT-1 are Fainter than their Progenitors

We present late-time Hubble and Spitzer Space Telescope imaging of SN 2008S and NGC 300 2008OT-1, the prototypes of a common class of stellar transients whose true nature is debated. Both objects are still fading and are now >15 times fainter than the progenitors in the mid-IR and are undetected in the optical and near-IR. Data from the Large Binocular Telescope and Magellan show that neither source has been variable in the optical since fading in 2010. We present models of surviving sources obscured by dusty shells or winds and find that extreme dust models are needed for surviving stars to be successfully hidden by dust. Explaining these transients as supernovae explosions, such as the electron capture supernovae believed to be associated with extreme AGB stars, seems an equally viable solution. Though SN 2008S is not detected in Chandra X-Ray Observatory data taken in 2012, the flux limits allow the fading IR source to be powered solely by the shock interaction of ejecta with the circumstellar medium if the shock velocity at the time of the observation was >20% slower than estimated from emission line widths while the transient was still optically bright. Continued SST monitoring and 10-20 micron observations with JWST can resolve any remaining ambiguities.

preprint2015arXiv

Discovery of Five Candidate Analogs for $η\,$Carinae in Nearby Galaxies

The late-stage evolution of very massive stars such as $η$ Carinae may be dominated by episodic mass ejections which may later lead to Type II superluminous supernova (SLSN-II; e.g., SN 2006gy). However, as long as $η$ Car is one of a kind, it is nearly impossible to quantitatively evaluate these possibilities. Here we announce the discovery of five objects in the nearby ($\sim4-8$ Mpc) massive star-forming galaxies M51, M83, M101 and NGC6946 that have optical through mid-IR photometric properties consistent with the hitherto unique $η$ Car. The Spitzer mid-IR spectral energy distributions of these $L_{bol}\simeq3-8\times10^{6} L_\odot$ objects rise steeply in the $3.6-8 μ$m bands, then turn over between $8$ and $24 μ$m, indicating the presence of warm ($\sim400-600$ K) circumstellar dust. Their optical counterparts in HST images are $\sim1.5-2$ dex fainter than their mid-IR peaks and require the presence of $\sim5-10 M_\odot$ of obscuring material. Our finding implies that the rate of $η$ Car-like events is a fraction $f=0.094$ ($0.040 < f < 0.21$ at $90\%$ confidence) of the core-collapse supernova (ccSN) rate. If there is only one eruption mechanism and SLSN-II are due to ccSN occurring inside these dense shells, then the ejection mechanism is likely associated with the onset of carbon burning ($\sim 10^3 - 10^4$ years) which is also consistent with the apparent ages of massive Galactic shells.

preprint2015arXiv

LOSS's First Supernova? New Limits on the "Impostor" SN 1997bs

We present new, late-time Hubble Space Telescope and Spitzer Space Telescope observations of the archetypal supernova (SN) impostor SN 1997bs. We show that SN 1997bs remains much fainter than its progenitor, posing a challenge for the canonical picture of late-time obscuration by dust forming in a shell ejected during the transient. The possibility that the star survived cloaked behind a dusty, steady wind is also disfavored. The simplest explanation is that SN 1997bs was a subluminous Type IIn SN, although it is currently impossible to rule out the possibility that the star survived either behind an obscuring dusty shell $\gtrsim 1 M_{\odot}$ or with a significantly decreased intrinsic luminosity. Continued monitoring of the system is the most promising means of ultimately resolving the question.

preprint2015arXiv

Massive stars exploding in a He-rich circumstellar medium - VII. The metamorphosis of ASASSN-15ed from a narrow line Type Ibn to a normal Type Ib Supernova

We present the results of the spectroscopic and photometric monitoring campaign of ASASSN-15ed. The transient was discovered quite young by the All Sky Automated Survey for SuperNovae (ASAS-SN). Amateur astronomers allowed us to sample the photometric SN evolution around maximum light, which we estimate to have occurred on JD = 2457087.4 +- 0.6 in the r-band. Its apparent r-band magnitude at maximum was r = 16.91 +- 0.10, providing an absolute magnitude M(r) ~ -20.04 +- 0.20, which is slightly more luminous than the typical magnitudes estimated for Type Ibn SNe. The post-peak evolution was well monitored, and the decline rate (being in most bands around 0.1 mag/d during the first 25 d after maximum) is marginally slower than the average decline rates of SNe Ibn during the same time interval. The object was initially classified as a Type Ibn SN because early-time spectra were characterized by a blue continuum with superimposed narrow P-Cygni lines of He I, suggesting the presence of a slowly moving (1200-1500 km/s), He-rich circumstellar medium. Later on, broad P-Cygni He I lines became prominent. The inferred velocities, as measured from the minimum of the broad absorption components, were between 6000 and 7000 km/s. As we attribute these broad features to the SN ejecta, this is the first time we have observed the transition of a Type Ibn SN to a Type Ib SN.

preprint2015arXiv

Optical observations of the luminous Type IIn Supernova 2010jl for over 900 days

The luminous Type IIn Supernova (SN) 2010jl shows strong evidence for the interaction of the SN ejecta with dense circumstellar material (CSM). We present observations of SN 2010jl for $t \sim 900$ d after its earliest detection, including a sequence of optical spectra ranging from $t = 55$ to $909$ d. We also supplement our late time spectra and the photometric measurements in the literature with an additional epoch of new, late time $BVRI$ photometry. Combining available photometric and spectroscopic data, we derive a semi-bolometric optical light curve and calculate a total radiated energy in the optical for SN 2010jl of $\sim 3.5\times10^{50}$ erg. We also examine the evolution of the H$α$ emission line profile in detail and find evidence for asymmetry in the profile for $t \gtrsim 775$ d that is not easily explained by any of the proposed scenarios for this fascinating event. Finally, we discuss the interpretations from the literature of the optical and near-infrared light curves, and propose that the most likely explanation of their evolution is the formation of new dust in the dense, pre-existing CSM wind after $\sim 300$ d.

preprint2015arXiv

Probing the dark matter radial profile in lens galaxies and the size of X-ray emitting region in quasars with microlensing

We use X-ray and optical microlensing measurements to study the shape of the dark matter density profile in the lens galaxies and the size of the (soft) X-ray emission region. We show that single epoch X-ray microlensing is sensitive to the source size. Our results, in good agreement with previous estimates, show that the size of the X-ray emission region scales roughly linearly with the black hole mass, with a half light radius of $R_{1/2}\simeq(24\pm14) r_g$ where $r_g=GM_{BH}/c^2$. This corresponds to a size of $\log(R_{1/2}/cm)=15.6^{+0.3}_{-0.3}$ or $\sim$ 1 light day for a black hole mass of $M_{BH}=10^9 M_\sun$. We simultaneously estimated the fraction of the local surface mass density in stars, finding that the stellar mass fraction is $α=0.20\pm0.05$ at an average radius of $\sim 1.9 R_{e}$, where $R_e$ is the effective radius of the lens. This stellar mass fraction is insensitive to the X-ray source size and in excellent agreement with our earlier results based on optical data. By combining X-ray and optical microlensing data, we can divide this larger sample into two radial bins. We find that the surface mass density in the form of stars is $α=0.31\pm0.15$ and $α=0.13\pm0.05$ at $(1.3\pm0.3) R_{e}$ and $(2.3\pm0.3) R_{e}$, respectively, in good agreement with expectations and some previous results.

preprint2015arXiv

Simulations of the OzDES AGN Reverberation Mapping Project

As part of the OzDES spectroscopic survey we are carrying out a large scale reverberation mapping study of $\sim$500 quasars over five years in the 30 deg$^2$ area of the Dark Energy Survey (DES) supernova fields. These quasars have redshifts ranging up to 4 and have apparent AB magnitudes between $16.8<r<22.5$ mag. The aim of the survey is to measure time lags between fluctuations in the quasar continuum and broad emission line fluxes of individual objects in order to measure black hole masses for a broad range of AGN and constrain the radius-luminosity ($R-L$) relationship. Here we investigate the expected efficiency of the OzDES reverberation mapping campaign and its possible extensions. We expect to recover lags for $\sim$35-45\% of the quasars. AGN with shorter lags and greater variability are more likely to yield a lag, and objects with lags $\lesssim$6 months or $\sim$1 year are expected be recovered the most accurately. The baseline OzDES reverberation mapping campaign is predicted to produce an unbiased measurement of the $R-L$ relationship parameters for H$β$, Mg II $λ$2798, and C IV $λ$1549. However, extending the baseline survey by either increasing the spectroscopic cadence, extending the survey season, or improving the emission line flux measurement accuracy will significantly improve the $R-L$ parameter constraints for all broad emission lines.

preprint2015arXiv

Space Telescope and Optical Reverberation Mapping Project. I. Ultraviolet Observations of the Seyfert 1 Galaxy NGC 5548 with the Cosmic Origins Spectrograph on Hubble Space Telescope

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis.

preprint2015arXiv

Space Telescope and Optical Reverberation Mapping Project. II. Swift and HST Reverberation Mapping of the Accretion Disk of NGC 5548

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a τpropto λ^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

preprint2015arXiv

Spitzer Point Source Catalogs of $\sim300,000$ Stars in Seven Nearby Galaxies

We present Spitzer IRAC $3.6-8$um and MIPS $24$um point-source catalogs for seven galaxies: NGC$6822$, M$33$, NGC$300$, NGC$2403$, M$81$, NGC$0247$, and NGC$7793$. The catalogs contain a total of $\sim300,000$ sources and were created by dual-band selection of sources with $>3σ$ detections at both $3.6$um and $4.5$um. The source lists become significantly incomplete near $m_{3.6}=m_{4.5}\simeq18$. We complement the $3.6$um and $4.5$um fluxes with $5.8$um, $8.0$um and $24$um fluxes or $3σ$ upper limits using a combination of PSF and aperture photometry. This catalog is a resource as an archive for studying mid-infrared transients and for planning observations with the James Webb Space Telescope.

preprint2015arXiv

The Cepheid distance to the maser-host galaxy NGC 4258: Studying systematics with the Large Binocular Telescope

We identify and phase a sample of 81 Cepheids in the maser-host galaxy NGC 4258 using the Large Binocular Telescope (LBT), and obtain calibrated mean magnitudes in up to 4 filters for a subset of 43 Cepheids using archival HST data. We employ 3 models to study the systematic effects of extinction, the assumed extinction law, and metallicity on the Cepheid distance to NGC 4258. We find a correction to the Cepheid colors consistent with a grayer extinction law in NGC 4258 compared to the Milky Way ($R_V =4.9$), although we believe this is indicative of other systematic effects. If we combine our Cepheid sample with previously known Cepheids, we find a significant metallicity adjustment to the distance modulus of $γ_1 = -0.61 \pm 0.21$ mag/dex, for the Zaritsky et al. (1994) metallicity scale, as well as a weak trend of Cepheid colors with metallicity. Conclusions about the absolute effect of metallicity on Cepheid mean magnitudes appear to be limited by the available data on the metallicity gradient in NGC 4258, but our Cepheid data require at least some metallicity adjustment to make the Cepheid distance consistent with independent distances to the LMC and NGC 4258. From our ensemble of models and the geometric maser distance of NGC 4258 ($μ_{N4258} = 29.40 \pm 0.06$ mag), we estimate $μ_{LMC} = 18.57 \pm 0.14$ mag ($51.82 \pm 3.23$ kpc).

preprint2015arXiv

The Search for Failed Supernovae with The Large Binocular Telescope: First Candidates

We are monitoring 27 galaxies within 10 Mpc using the Large Binocular Telescope to search for failed supernovae (SNe), massive stars that collapse to form a black hole without a SN explosion. We present the results from the first 4 years of survey data, during which these galaxies were observed to produce 3 successful core-collapse SNe. We search for stars that have "vanished" over the course of our survey, by examining all stars showing a decrease in luminosity of $ΔνL_ν \ge 10^4L_{\odot}$ from the first to the last observation. We also search for the low luminosity, long duration transients predicted by \citet{lovegrove2013} for failed explosions of red supergiants. After analyzing the first 4 years of data in this first direct search for failed SNe, we are left with one candidate requiring further study. This candidate has an estimated mass of 18-25$M_{\odot}$, a mass range likely associated with failed SNe and, if real, implies that failed SN represent a median fraction of $f \simeq 0.30$ of core-collapses, with symmetric 90% confidence limits of $0.07 \le f \le 0.62$. If follow up data eliminate this candidate, we find an upper limit on the fraction of core collapses leading to a failed SN of $f<0.40$ at 90% confidence. As the duration of the survey continues to increase, it will begin to constrain the $f \simeq$ 10-30% failure rates needed to explain the deficit of massive SN progenitors and the observed black hole mass function.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Rapid CIV Broad Absorption Line Variability

We report the discovery of rapid variations of a high-velocity CIV broad absorption line trough in the quasar SDSS J141007.74+541203.3. This object was intensively observed in 2014 as a part of the Sloan Digital Sky Survey Reverberation Mapping Project, during which 32 epochs of spectroscopy were obtained with the Baryon Oscillation Spectroscopic Survey spectrograph. We observe significant (>4sigma) variability in the equivalent width of the broad (~4000 km/s wide) CIV trough on rest-frame timescales as short as 1.20 days (~29 hours), the shortest broad absorption line variability timescale yet reported. The equivalent width varied by ~10% on these short timescales, and by about a factor of two over the duration of the campaign. We evaluate several potential causes of the variability, concluding that the most likely cause is a rapid response to changes in the incident ionizing continuum. If the outflow is at a radius where the recombination rate is higher than the ionization rate, the timescale of variability places a lower limit on the density of the absorbing gas of n_e > 3.9 x 10^5 cm^-3. The broad absorption line variability characteristics of this quasar are consistent with those observed in previous studies of quasars, indicating that such short-term variability may in fact be common and thus can be used to learn about outflow characteristics and contributions to quasar/host-galaxy feedback scenarios.

preprint2015arXiv

Total eclipse of the heart: The AM CVn Gaia14aae / ASSASN-14cn

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

preprint2014arXiv

ASASSN-14ae: A Tidal Disruption Event at 200 Mpc

ASASSN-14ae is a candidate tidal disruption event (TDE) found at the center of SDSS J110840.11+340552.2 ($d\simeq200$~Mpc) by the All-Sky Automated Survey for Supernovae (ASAS-SN). We present ground-based and Swift follow-up photometric and spectroscopic observations of the source, finding that the transient had a peak luminosity of $L\simeq8\times10^{43}$~erg~s$^{-1}$ and a total integrated energy of $E\simeq1.7\times10^{50}$ ergs radiated over the $\sim5$ months of observations presented. The blackbody temperature of the transient remains roughly constant at $T\sim20,000$~K while the luminosity declines by nearly 1.5 orders of magnitude during this time, a drop that is most consistent with an exponential, $L\propto e^{-t/t_0}$ with $t_0\simeq39$~days. The source has broad Balmer lines in emission at all epochs as well as a broad He II feature emerging in later epochs. We compare the color and spectral evolution to both supernovae and normal AGN to show that {\name} does not resemble either type of object and conclude that a TDE is the most likely explanation for our observations. At $z=0.0436$, ASASSN-14ae is the lowest-redshift TDE candidate discovered at optical/UV wavelengths to date, and we estimate that ASAS-SN may discover $0.1 - 3$ of these events every year in the future.

preprint2014arXiv

Constraints on Core Collapse from the Black Hole Mass Function

We model the observed black hole mass function under the assumption that black hole formation is controlled by the compactness of the stellar core at the time of collapse. Low compactness stars are more likely to explode as supernovae and produce neutron stars, while high compactness stars are more likely to be failed supernovae that produce black holes with the mass of the helium core of the star. Using three sequences of stellar models and marginalizing over a model for the completeness of the black hole mass function, we find that the compactness xi(2.5) above which 50% of core collapses produce black holes is xi(2.5)=0.24 (0.15 < xi(2.5) < 0.37) at 90% confidence). While models with a sharp transition between successful and failed explosions are always the most likely, the width of the transition between the minimum compactness for black hole formation and the compactness above which all core collapses produce black holes is not well constrained. The models also predict that f=0.18 (0.09 < f < 0.39) of core collapses fail assuming a minimum mass for core collapse of 8Msun. We tested four other criteria for black hole formation based on xi(2.0) and xi(3.0), the compactnesses at enclosed masses of 2.0 or 3.0 rather than 2.5Msun, the mass of the iron core, and the mass inside the oxygen burning shell. We found that xi(2.0) works as well as xi(2.5), while the compactness xi(3.0) works significantly worse, as does using the iron core mass or the mass enclosed by the oxygen burning shell. As expected from the high compactness of 20-25Msun stars, black hole formation in this mass range provides a natural explanation of the red supergiant problem.

preprint2014arXiv

Dark Matter Fraction in Lens Galaxies: New Estimates from Microlensing

We present a joint estimate of the stellar/dark matter mass fraction in lens galaxies and the average size of the accretion disk of lensed quasars from microlensing measurements of 27 quasar image pairs seen through 19 lens galaxies. The Bayesian estimate for the fraction of the surface mass density in the form of stars is $α=0.21\pm0.14$ near the Einstein radius of the lenses ($\sim 1 - 2$ effective radii). The estimate for the average accretion disk size is $R_{1/2}=7.9^{+3.8}_{-2.6}\sqrt{M/0.3M_\sun}$ light days. The fraction of mass in stars at these radii is significantly larger than previous estimates from microlensing studies assuming quasars were point-like. The corresponding local dark matter fraction of 79\% is in good agreement with other estimates based on strong lensing or kinematics. The size of the accretion disk inferred in the present study is slightly larger than previous estimates.

preprint2014arXiv

Dust Formation By Failed Supernovae

We consider dust formation during the ejection of the hydrogen envelope of a red supergiant during a failed supernova (SN) creating a black hole. While the dense, slow moving ejecta are very efficient at forming dust, only the very last phases of the predicted visual transient will be obscured. The net grain production consists of ~0.01 solar masses of very large grains (10 to 1000 microns). This means that failed SNe could be the source of the very large extrasolar dust grains identified by Ulysses, Galileo and radar studies of meteoroid re-entry trails rather than their coming from an ejection process associated with protoplanetary or other disks.

preprint2014arXiv

Dust Formation in the Presence of Photons I: Evaporation Rates for Small Dust Grains

The temperature of newly forming dust is controlled by the radiation field. As dust forms around stars, stellar transients, quasars or supernovae, the grains must grow through a regime where they are stochastically heated by individual photons. Since evaporation rates increase exponentially with temperature while cooling times decrease only as a power law, the evaporation rates for these small grains are dominated by the temperature spikes. We calculate effective evaporation temperatures for a broad range of input spectra that are encapsulated in a series of simple interpolation formulae for both graphitic and silicate grains. These can be easily used to first determine if dust formation is possible and then to estimate the radius or time at which it commences for a broad range of radiation environments. With these additional physical effects, very small grains may form earlier than in standard models of AGB winds. Even for very high mass loss rates, the hottest stars that can form dust are G and F stars particularly in the case of silicate dusts. For hotter stars, the higher fluxes of ultraviolet photons prevent dust formation. Thus, episodic dust formation by OH/IR stars and LBVs is primarily driven by fluctuations in their apparent temperatures rather than changes in luminosity or mass loss rates.

preprint2014arXiv

Finding $η$ Car Analogs in Nearby Galaxies Using Spitzer: II. Identification of An Emerging Class of Extragalactic Self-Obscured Stars

Understanding the late-stage evolution of the most massive stars such as $η$ Carinae is challenging because no true analogs of $η$ Car have been clearly identified in the Milky Way or other galaxies. In Khan et. al. (2013), we utilized Spitzer IRAC images of $7$ nearby ($\lesssim4$ Mpc) galaxies to search for such analogs, and found $34$ candidates with flat or red mid-IR spectral energy distributions. Here, in Paper II, we present our characterization of these candidates using multi-wavelength data from the optical through the far-IR. Our search detected no true analogs of $η$ Car, which implies an eruption rate that is a fraction $0.01\lesssim F \lesssim 0.19$ of the ccSN rate. This is roughly consistent with each $M_{ZAMS} \gtrsim 70M_\odot$ star undergoing $1$ or $2$ outbursts in its lifetime. However, we do identify a significant population of $18$ lower luminosity $\left(\log(L/L_\odot)\simeq5.5-6.0\right)$ dusty stars. Stars enter this phase at a rate that is fraction $0.09 \lesssim F \lesssim 0.55$ of the ccSN rate, and this is consistent with all $25 < M_{ZAMS} < 60M_\odot$ stars undergoing an obscured phase at most lasting a few thousand years once or twice. These phases constitute a negligible fraction of post-main sequence lifetimes of massive stars, which implies that these events are likely to be associated with special periods in the evolution of the stars. The mass of the obscuring material is of order $\sim M_\odot$, and we simply do not find enough heavily obscured stars for theses phases to represent more than a modest fraction ($\sim 10\%$ not $\sim 50\%$) of the total mass lost by these stars. In the long term, the sources that we identified will be prime candidates for detailed physical analysis with JWST.

preprint2014arXiv

Reverberation Mapping of the Seyfert 1 Galaxy NGC 7469

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

preprint2014arXiv

Stellar Mergers Are Common

The observed Galactic rate of stellar mergers or the initiation of common envelope phases brighter than M_V=-3 (M_I=-4) is of order 0.5 (0.3)/year with 90% confidence statistical uncertainties of 0.24-1.1 (0.14-0.65) and factor of 2 systematic uncertainties. The (peak) luminosity function is roughly dN/dL L^(-1.4+/-0.3), so the rates for events more luminous than V1309 Sco (M_V=-7 mag) or V838Mon (M_V=-10 mag) are lower at r~0.1/year and 0.03/year, respectively. The peak luminosity is a steep function of progenitor mass, L M^(2-3). This very roughly parallels the scaling of luminosity with mass on the main sequence, but the transients are ~2000-4000 times more luminous at peak. Combining these, the mass function of the progenitors, dN/dM M^(-2.0+/-0.8), is consistent with the initial mass function, albeit with broad uncertainties. These observational results are also broadly consistent with the estimates of binary population synthesis models. While extragalactic variability surveys can better define the rates and properties of the high luminosity events, systematic, moderate depth (I>16 mag) surveys of the Galactic plane are needed to characterize the low luminosity events. The existing Galactic samples are only ~20% complete and Galactic surveys are (at best) reaching a typical magnitude limit of <13 mag.

preprint2014arXiv

The Average Size and Temperature Profile of Quasar Accretion Disks

We use multi-wavelength microlensing measurements of a sample of 10 image pairs from 8 lensed quasars to study the structure of their accretion disks. By using spectroscopy or narrow band photometry we have been able to remove contamination from the weakly microlensed broad emission lines, extinction and any uncertainties in the large-scale macro magnification of the lens model. We determine a maximum likelihood estimate for the exponent of the size versus wavelength scaling ($r_s\propto λ^p$ corresponding to a disk temperature profile of $T\propto r^{-1/p}$) of $p=0.75^{+0.2}_{-0.2}$, and a Bayesian estimate of $p=0.8\pm0.2$, which are significantly smaller than the prediction of thin disk theory ($p=4/3$). We have also obtained a maximum likelihood estimate for the average quasar accretion disk size of $r_s=4.5^{+1.5}_{-1.2} $ lt-day at a rest frame wavelength of $λ= 1026~{\mathrm Å}$ for microlenses with a mean mass of $M=1 M_\sun$, in agreement with previous results, and larger than expected from thin disk theory.

preprint2014arXiv

The Clustering and Halo Masses of Star Forming Galaxies at z<1

We present clustering measurements and halo masses of star forming galaxies at 0.2 < z < 1.0. After excluding AGN, we construct a sample of 22553 24 μm sources selected from 8.42 deg^2 of the Spitzer MIPS AGN and Galaxy Evolution Survey of Boötes. Mid-infrared imaging allows us to observe galaxies with the highest star formation rates (SFRs), less biased by dust obscuration afflicting the optical bands. We find that the galaxies with the highest SFRs have optical colors which are redder than typical blue cloud galaxies, with many residing within the green valley. At z > 0.4 our sample is dominated by luminous infrared galaxies (LIRGs, L_TIR > 10^11 Lsun) and is comprised entirely of LIRGs and ultra-luminous infrared galaxies (ULIRGs, L_TIR > 10^12 Lsun) at z > 0.6. We observe weak clustering of r_0 = 3-6 Mpc/h for almost all of our star forming samples. We find that the clustering and halo mass depend on L_TIR at all redshifts, where galaxies with higher L_TIR (hence higher SFRs) have stronger clustering. Galaxies with the highest SFRs at each redshift typically reside within dark matter halos of M_halo ~ 10^12.9 Msun/h. This is consistent with a transitional halo mass, above which star formation is largely truncated, although we cannot exclude that ULIRGs reside within higher mass halos. By modeling the clustering evolution of halos, we connect our star forming galaxy samples to their local descendants. Most star forming galaxies at z < 1.0 are the progenitors of L < 2.5L* blue galaxies in the local universe, but star forming galaxies with the highest SFRs (L_TIR >10^11.7 Lsun) at 0.6<z<1.0 are the progenitors of early-type galaxies in denser group environments.

preprint2014arXiv

The Man Behind the Curtain: X-rays Drive the UV through NIR Variability in the 2013 AGN Outburst in NGC 2617

After the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered a significant brightening of the inner region of NGC 2617, we began a ~70 day photometric and spectroscopic monitoring campaign from the X-ray through near-infrared (NIR) wavelengths. We report that NGC 2617 went through a dramatic outburst, during which its X-ray flux increased by over an order of magnitude followed by an increase of its optical/ultraviolet (UV) continuum flux by almost an order of magnitude. NGC 2617, classified as a Seyfert 1.8 galaxy in 2003, is now a Seyfert 1 due to the appearance of broad optical emission lines and a continuum blue bump. Such "changing look Active Galactic Nuclei (AGN)" are rare and provide us with important insights about AGN physics. Based on the Hbeta line width and the radius-luminosity relation, we estimate the mass of central black hole to be (4 +/- 1) x 10^7 M_sun. When we cross-correlate the light curves, we find that the disk emission lags the X-rays, with the lag becoming longer as we move from the UV (2-3 days) to the NIR (6-9 days). Also, the NIR is more heavily temporally smoothed than the UV. This can largely be explained by a simple model of a thermally emitting thin disk around a black hole of the estimated mass that is illuminated by the observed, variable X-ray fluxes.

preprint2013arXiv

Characterizing a Dramatic $ΔV\sim-9$ Flare on an Ultracool Dwarf Found by the ASAS-SN Survey

We analyze a $ΔV\sim-9$ magnitude flare on the newly identified M8 dwarf SDSS J022116.84+194020.4 (hereafter SDSSJ0221) detected as part of the All-Sky Automated Survey for Supernovae (ASAS-SN). Using infrared and optical spectra, we confirm that SDSSJ0221 is a relatively nearby (d$\sim$76 pc) M8 dwarf with strong quiescent H$α$ emission. Based on kinematics and the absence of features consistent with low-gravity (young) ultracool dwarfs, we place a lower limit of 200 Myr on the age of SDSSJ0221. When modeled with a simple, classical flare light-curve, this flare is consistent with a total $U$-band flare energy $E_U\sim$ 10$^{34}$ erg, confirming that the most dramatic flares are not limited to warmer, more massive stars. Scaled to include a rough estimate of the emission line contribution to the $V$ band, we estimate a blackbody filling factor of $\sim$$10-30\%$ during the flare peak and $\sim$$0.5-1.6\%$ during the flare decay phase. These filling factors correspond to flare areas that are an order of magnitude larger than those measured for most mid-M dwarf flares.

preprint2013arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XIII: Time delays and 9-yr optical monitoring of the lensed quasar RX J1131-1231

We present the results from nine years of optically monitoring the gravitationally lensed z=0.658 quasar RX J1131-1231. The R-band light curves of the four individual images of the quasar were obtained using deconvolution photometry for a total of 707 epochs. Several sharp quasar variability features strongly constrain the time delays between the quasar images. Using three different numerical techniques, we measure these delays for all possible pairs of quasar images while always processing the four light curves simultaneously. For all three methods, the delays between the three close images A, B, and C are compatible with being 0, while we measure the delay of image D to be 91 days, with a fractional uncertainty of 1.5% (1 sigma), including systematic errors. Our analysis of random and systematic errors accounts in a realistic way for the observed quasar variability, fluctuating microlensing magnification over a broad range of temporal scales, noise properties, and seasonal gaps. Finally, we find that our time-delay measurement methods yield compatible results when applied to subsets of the data.

preprint2013arXiv

Dark energy with gravitational lens time delays

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

preprint2013arXiv

Failed Supernovae Explain the Compact Remnant Mass Function

One explanation for the absence of higher mass red supergiants (16.5 Msun < M < 25Msun) as the progenitors of Type IIP supernovae (SNe) is that they die in failed SNe creating black holes. Simulations show that such failed SNe still eject their hydrogen envelopes in a weak transient, leaving a black hole with the mass of the star's helium core (5-8Msun). Here we show that this naturally explains the typical masses of observed black holes and the gap between neutron star and black hole masses without any fine-tuning of the SN mechanism beyond having it fail in a mass range where many progenitor models have density structures that make the explosions more likely to fail. There is no difficulty including this ~20% population of failed SNe in any accounting of SN types over the progenitor mass function. And, other than patience, there is no observational barrier to either detecting these black hole formation events or limiting their rates to be well below this prediction.

preprint2013arXiv

Is Quasar Optical Variability a Damped Random Walk?

The damped random walk (DRW) model is increasingly used to model the variability in quasar optical light curves, but it is still uncertain whether the DRW model provides an adequate description of quasar optical variability across all time scales. Using a sample of OGLE quasar light curves, we consider four modifications to the DRW model by introducing additional parameters into the covariance function to search for deviations from the DRW model on both short and long time scales. We find good agreement with the DRW model on time scales that are well sampled by the data (from a month to a few years), possibly with some intrinsic scatter in the additional parameters, but this conclusion depends on the statistical test employed and is sensitive to whether the estimates of the photometric errors are correct to within ~10%. On very short time scales (below a few months), we see some evidence of the existence of a cutoff below which the correlation is stronger than the DRW model, echoing the recent finding of Mushotzky et al. (2011) using quasar light curves from Kepler. On very long time scales (> a few years), the light curves do not constrain models well, but are consistent with the DRW model.

preprint2013arXiv

Microlensing of Quasar Broad Emission Lines: Constraints on Broad Line Region Size

We measure the differential microlensing of the broad emission lines between 18 quasar image pairs in 16 gravitational lenses. We find that high ionization lines such as CIV are more strongly microlensed than low ionization lines, indicating that the high ionization line emission regions are more compact. If we statistically model the distribution of microlensing magnifications, we obtain estimates for the broad line region radius of 24 (-15/+22) and 55 (-35/+150) light-days (90% confidence) for the high and low ionization lines, respectively. When the sample is divided attending to quasar luminosity, we find that the line emission regions of more luminous quasars are larger, with a slope consistent with the expected scaling from photoionization models. Our estimates also agree well with the results from local reveberation mapping studies.

preprint2013arXiv

Microlensing of Quasar UV Iron Emission

We measure the differential microlensing of the UV Fe II and Fe III emission line blends between 14 quasar image pairs in 13 gravitational lenses. We find that the UV iron emission is strongly microlensed in 4 cases with amplitudes comparable to that of the continuum. Statistically modeling the magnifications we infer a typical size of r ~ 4*sqrt(M/Msun) light-days for the Fe line emitting regions which is comparable to the size of the region generating the UV continuum (3 ~ 7 light-days). This may indicate that a significant part of the UV Fe II and Fe III emission originates in the quasar accretion disk.

preprint2013arXiv

Observing the Next Galactic Supernova

We model the distance, extinction, and magnitude probability distributions of a successful Galactic core-collapse supernova (ccSN), its shock breakout radiation, and its massive star progenitor. We find, at very high probability (~100%), that the next Galactic SN will easily be detectable in the near-IR and that near-IR photometry of the progenitor star very likely (~92%) already exists in the 2MASS survey. Most ccSNe (~98%) will be easily observed in the optical, but a significant fraction (~43%) will lack observations of the progenitor due to a combination of survey sensitivity and confusion. If neutrino detection experiments can quickly disseminate a likely position (~3 deg), we show that a modestly priced IR camera system can probably detect the shock breakout radiation pulse even in daytime (~64% for the cheapest design). Neutrino experiments should seriously consider adding such systems, both for their scientific return and as an added and internal layer of protection against false triggers. We find that shock breakouts from failed ccSNe of red supergiants may be more observable than those of successful SNe. We review the process by which neutrinos from a Galactic ccSN would be detected and announced. We provide new information on the EGADS system and its potential for providing instant neutrino alerts. We also discuss the distance, extinction, and magnitude probability distributions for the next Galactic Type Ia SN. Based on our modeled observability, we find a Galactic ccSN rate of 3.2 (+7.3/-2.6) per century and a Galactic Type Ia SN rate of 1.4 (+1.4/-0.8) per century for a total Galactic SN rate of 4.6 (+7.4/-2.7) per century is needed to account for the SNe observed over the last millennium.

preprint2013arXiv

The Size of the Narrow-Line Emitting Region in the Seyfert 1 Galaxy NGC 5548 from Emission-Line Variability

The narrow [O III] 4959, 5007 emission-line fluxes in the spectrum of the well-studied Seyfert 1 galaxy NGC 5548 are shown to vary with time. From this we show that the narrow line-emitting region has a radius of only 1-3 pc and is denser (n ~ 10^5 cm^{-3}) than previously supposed. The [O III] line width is consistent with virial motions at this radius given previous determinations of the black hole mass.Since the [O III] emission-line flux is usually assumed to be constant and is therefore used to calibrate spectroscopic monitoring data, the variability has ramifications for the long-term secular variations of continuum and emission-line fluxes, though it has no effect on shorter-term reverberation studies. We present corrected optical continuum and broad Hbeta emission-line light curves for the period 1988 to 2008.

preprint2013arXiv

The Structure of the X-ray and Optical Emitting Regions of the Lensed Quasar Q 2237+0305

We use gravitational microlensing to determine the size of the X-ray and optical emission regions of the quadruple lens system Q 2237+0305. The optical half-light radius, log(R_{1/2,V}/cm)=16.41\pm0.18 (at lambda_{rest}=2018 Å), is significantly larger than the observed soft, log(R_{1/2,soft}/cm)=15.76^{+0.41}_{-0.34} (1.1-3.5 keV in the rest frame), and hard, log(R_{1/2,hard}/cm)=15.46^{+0.34}_{-0.29} (3.5-21.5 keV in the rest frame), band X-ray emission. There is a weak evidence that the hard component is more compact than the soft, with log(R_{1/2,soft}/R_{1/2,hard}) \sim 0.30^{+0.53}_{-0.45}. This wavelength-dependent structure agrees with recent results found in other lens systems using microlensing techniques, and favors geometries in which the corona is concentrated near the inner edge of the accretion disk. While the available measurements are limited, the size of the X-ray emission region appears to be roughly proportional to the mass of the central black hole.

preprint2013arXiv

Time Delay and Accretion Disk Size Measurements in the Lensed Quasar SBS 0909+532 from Multiwavelength Microlensing Analysis

We present three complete seasons and two half-seasons of SDSS r-band photometry of the gravitationally lensed quasar SBS 0909+532 from the U.S. Naval Observatory, as well as two seasons each of SDSS g-band and r-band monitoring from the Liverpool Robotic Telescope. Using Monte Carlo simulations to simultaneously measure the system's time delay and model the r-band microlensing variability, we confirm and significantly refine the precision of the system's time delay to Δt_{AB} = 50^{+2}_{-4} days, where the stated uncertainties represent the bounds of the formal 1σ confidence interval. There may be a conflict between the time delay measurement and a lens consisting of a single galaxy. While models based on the Hubble Space Telescope astrometry and a relatively compact stellar distribution can reproduce the observed delay, the models have somewhat less dark matter than we would typically expect. We also carry out a joint analysis of the microlensing variability in the r- and g-bands to constrain the size of the quasar's continuum source at these wavelengths, obtaining log[(r_{s,r}/cm) [cos{i}/0.5]^{1/2}] = 15.3 \pm 0.3 and log[(r_{s,g}/cm) [cos{i}/0.5]^{1/2}] = 14.8 \pm 0.9, respectively. Our current results do not formally constrain the temperature profile of the accretion disk but are consistent with the expectations of standard thin disk theory.

preprint2013arXiv

Type Ia Single Degenerate Survivors Must Be Overluminous

In the single-degenerate (SD) channel of a Type Ia supernovae (SN Ia) explosion, a main-sequence (MS) donor star survives the explosion but it is stripped of mass and shock heated. An essentially unavoidable consequence of mass loss during the explosion is that the companion must have an overextended envelope after the explosion. While this has been noted previously, it has not been strongly emphasized as an inevitable consequence. We calculate the future evolution of the companion by injecting 2-6 10^47 ergs into the stellar evolution model of a 1 Msun donor star based on the post-explosion progenitors seen in simulations. We find that, due to the Kelvin-Helmholtz collapse of the envelope, the companion must become significantly more luminous (10 - 10^3 Lsun) for a long period of time (10^3 - 10^4 years). The lack of such a luminous "leftover" star in the LMC supernova remnant SNR 0609-67.5 provides another piece of evidence against the SD scenario. We also show that none of the stars proposed as the survivors of the Tycho supernova, including Tycho G, could plausibly be the donor star. Additionally, luminous donors closer than ~10 Mpc should be observable with the Hubble Space Telescope starting ~2 years post-peak. Such systems include SN 1937C, SN 1972E, SN 1986G, and SN 2011fe. Thus, the SD channel is already ruled out for at least two nearby SNe Ia and can be easily tested for a number of additional ones. We also discuss similar implications for the companions of core-collapse SNe.

preprint2012arXiv

A Robust Determination of the size of quasar accretion disks using gravitational microlensing

Using microlensing measurements from a sample of 27 image-pairs of 19 lensed quasars we determine a maximum likelihood estimate for the accretion disk size of an {{\em}average} quasar of $r_s=4.0^{+2.4}_{-3.1} $ light days at rest frame $<λ>=1736$Å for microlenses with a mean mass of $<M>=0.3M_\odot$. This value, in good agreement with previous results from smaller samples, is roughly a factor of 5 greater than the predictions of the standard thin disk model. The individual size estimates for the 19 quasars in our sample are also in excellent agreement with the results of the joint maximum likelihood analysis.

preprint2012arXiv

A Study of Cepheids in M81 with the Large Binocular Telescope (Efficiently Calibrated with HST)

We identify and phase a sample of 107 Cepheids with 10<P/days<100 in M81 using the LBT and calibrate their BVI mean magnitudes with archival HST data. The use of a ground-based telescope to identify and phase the Cepheids and HST only for the final calibration reduces the demand on HST by nearly an order of magnitude and yields Period-Luminosity (PL) relations with dispersions comparable to the best LMC samples. We fit the sample using the OGLE-II LMC PL relations and are unable to find a self-consistent distance for different band combinations or radial locations within M81. We can do so after adding a radial dependence to the PL zero point that corresponds to a luminosity dependence on metallicity of g_mu=-0.56+/-0.36 mag/dex. We find marginal evidence for a shift in color as a function of metallicity, distinguishable from the effects of extinction, of g_2=+0.07+/-0.03 mag/dex. We find a distance modulus for M81, relative to the LMC, of mu(M81-LMC)=9.39+/-0.14 mag, including uncertainties due to the metallicity corrections. This corresponds to a distance to M81 of 3.6+/-0.2 Mpc, assuming a LMC distance modulus of 18.41 mag. We carry out a joint analysis of M81 and NGC4258 Cepheids and simultaneously solve for the distance of M81 relative to NGC4258 and the metallicity corrections. Given the current data, the uncertainties of such joint fits are dominated by the relative metallicities and the abundance gradients rather than by measurement errors of the Cepheid magnitudes or colors. We find mu(M81-LMC)=9.40 (-0.11/+0.15) mag, mu(N4258-LMC)=11.08 (-0.17/+0.21) mag and mu(N4258-M81)=1.68+/-0.08 mag and joint metallicity corrections of g_mu=-0.62 (-0.35/+0.31) mag/dex and g_2=0.01+/-0.01 mag/dex. Quantitative analyses of Cepheid distances must take into account both the metallicity dependencies of the Cepheids and the uncertainties in the abundance estimates. (ABRIDGED)

preprint2012arXiv

Dust to Dust: 3 years in the Evolution of the Unusual SN 2008S

We obtained late-time optical and near-IR imaging of SN 2008S with the Large Binocular Telescope (LBT), near-IR data with the Hubble Space Telescope (HST), and mid-IR data with the Spitzer Space Telescope (SST). We find that (1) it is again invisible at optical (UBVR) wavelengths to magnitude limits of approximately 25 mag, (2) while detected in the near-IR (H) at approximately 24.8 mag, it is fading rapidly, and (3) it is still brighter than the progenitor at 3.6 and 4.5 microns in the mid-IR with a slow, steady decline. The IR detections in December 2010 are consistent with dust emission at a blackbody temperature of T ~ 640 K and a total luminosity of L ~ 200000 Lsun, much higher than the L ~ 40000 Lsun luminosity of the obscured progenitor star. The local environment also shows no evidence for massive (M >= 10 Msun) stars in the vicinity of the transient, consistent with the progenitor being a massive AGB star.

preprint2012arXiv

Finding η Car Analogs in Nearby Galaxies Using Spitzer: I. Candidate Selection

The late-stage evolution of the most massive stars such as η Carinae is controlled by the effects of mass loss, which may be dominated by poorly understood eruptive mass ejections. Understanding this population is challenging because no true analogs of η Car have been clearly identified in the Milky Way or other galaxies. We utilize Spitzer IRAC images of 7 nearby (=<4 Mpc) galaxies to search for such analogs. We find 34 candidates with a flat or rising mid-IR spectral energy distributions towards longer mid-infrared wavelengths that emit 10^5 L_sun in the IRAC bands (3.6 to 8.0 micron) and are not known to be background sources. Based on our estimates for the expected number of background sources, we expect that follow-up observations will show that most of these candidates are not dust enshrouded massive stars, with an expectation of only 6+-6 surviving candidates. Since we would detect true analogs of η Car for roughly 200 years post-eruption, this implies that the rate of eruptions like η Car is less than the ccSN rate. It is possible, however, that every M > 40 M_sun star undergoes such eruptions given our initial results. In Paper II we will characterize the candidates through further analysis and follow-up observations, and there is no barrier to increasing the galaxy sample by an order of magnitude.

preprint2012arXiv

On Absorption by Circumstellar Dust, With the Progenitor of SN2012aw as a Case Study

We use the progenitor of SN2012aw to illustrate the consequences of modeling circumstellar dust using Galactic (interstellar) extinction laws that (1) ignore dust emission in the near-IR and beyond; (2) average over dust compositions, and (3) mis-characterize the optical/UV absorption by assuming that scattered photons are lost to the observer. The primary consequences for the progenitor of SN2012aw are that both the luminosity and the absorption are significantly over-estimated. In particular, the stellar luminosity is most likely in the range 10^4.8 < L/Lsun < 10^5.0 and the star was not extremely massive for a Type IIP progenitor, with M < 15Msun. Given the properties of the circumstellar dust and the early X-ray/radio detections of SN2012aw, the star was probably obscured by an on-going wind with Mdot ~ 10^-5.5 to 10^-5.0 Msun/year at the time of the explosion, roughly consistent with the expected mass loss rates for a star of its temperature (T_* ~ 3600K) and luminosity. In the spirit of Galactic extinction laws, we supply simple interpolation formulas for circumstellar extinction by dusty graphitic and silicate shells as a function of wavelength (>0.3 micron) and total (absorption plus scattering) V-band optical depth (tau < 20). These do not include the contributions of dust emission, but provide a simple, physical alternative to incorrectly using interstellar extinction laws.

preprint2012arXiv

Revealing the Structure of an Accretion Disk Through Energy Dependent X-ray Microlensing

We present results from monitoring observations of the gravitationally lensed quasar RX J1131-1231 performed with the Chandra X-ray Observatory. The X-ray observations were planned with relatively long exposures that allowed a search for energy-dependent microlensing in the soft (0.2-2 keV) and hard (2-10 keV) light curves of the images of RX J1131-1231. We detect significant microlensing in the X-ray light-curves of images A and D, and energy-dependent microlensing of image D. The magnification of the soft band appears to be larger than that in the hard band by a factor of ~ 1.3 when image D becomes more magnified. This can be explained by the difference between a compact, softer-spectrum corona that is producing a more extended, harder spectrum reflection component off the disk. This is supported by the evolution of the fluorescent iron line in image D over three consecutive time-averaged phases of the light curve. In the first period, a Fe line at E = 6.36(-0.16,+0.13) keV is detected (at > 99% confidence). In the second period, two Fe lines are detected, one at E = 5.47(-0.08,+0.06) keV (detected at > 99% confidence) and another at E = 6.02(-0.07,+0.09) keV (marginally detected at > 90% confidence), and in the third period, a broadened Fe line at 6.42(-0.15,+0.19) keV is detected (at > 99% confidence). This evolution of the Fe line profile during the microlensing event is consistent with the line distortion expected when a caustic passes over the inner disk where the shape of the fluorescent Fe line is distorted by General Relativistic and Doppler effects.

preprint2012arXiv

Reverberation Mapping Results for Five Seyfert 1 Galaxies

We present the results from a detailed analysis of photometric and spectrophotometric data on five Seyfert 1 galaxies observed as a part of a recent reverberation mapping program. The data were collected at several observatories over a 140-day span beginning in 2010 August and ending in 2011 January. We obtained high sampling-rate light curves for Mrk 335, Mrk 1501, 3C120, Mrk 6, and PG2130+099, from which we have measured the time lag between variations in the 5100 Angstrom continuum and the H-beta broad emission line. We then used these measurements to calculate the mass of the supermassive black hole at the center of each of these galaxies. Our new measurements substantially improve previous measurements of MBH and the size of the broad line-emitting region for four sources and add a measurement for one new object. Our new measurements are consistent with photoionization physics regulating the location of the broad line region in active galactic nuclei.

preprint2012arXiv

SN2002bu -- Another SN2008S-like Transient

We observed SN2002bu in the near-IR with the Hubble Space Telescope, the mid-IR with the Spitzer Space Telescope and in X-rays with Swift 10 years after the explosion. If the faint L_H\sim100 Lsun HST near-IR source at the transient position is the near-IR counterpart of SN2002bu, then the source has dramatically faded between 2004 and 2012, from L\sim10^6.0 Lsun to L\sim10^4.5 Lsun. It is still heavily obscured, tau_V\sim5 in graphitic dust models, with almost all the energy radiated in the mid-IR. The radius of the dust emission is increasing as R\simt^(0.7+/-0.4) and the optical depth is dropping as tau_V\simt^(-1.3+/-0.4). The evolution expected for an expanding shell of material, tau_V\sim1/t^2, is ruled out at approximately 2 sigma while the tau_V\simt^(-0.8) to t^(-1) optical depth scaling for a shock passing through a pre-existing wind is consistent with the data. If the near-IR source is a chance superposition, the present day source can be moderately more luminous, significantly more obscured and evolving more slowly. While we failed to detect X-ray emission, the X-ray flux limits are consistent with the present day emissions being powered by an expanding shock wave. SN2002bu is clearly a member of the SN2008S class of transients, but continued monitoring of the evolution of the spectral energy distribution is needed to conclusively determine the nature of the transient.

preprint2012arXiv

The Clustering of Extremely Red Objects

We measure the clustering of Extremely Red Objects (EROs) in ~8 deg^2 of the NOAO Deep Wide Field Survey Boötes field in order to establish robust links between ERO z~1.2 and local galaxy z<0.1 populations. Three different color selection criteria from the literature are analyzed to assess the consequences of using different criteria for selecting EROs. Specifically, our samples are (R-K_s)>5.0 (28,724 galaxies), (I-K_s)>4.0 (22,451 galaxies) and (I-[3.6])>5.0 (64,370 galaxies). Magnitude-limited samples show the correlation length (r_0) to increase for more luminous EROs, implying a correlation with stellar mass. We can separate star-forming and passive ERO populations using the (K_s-[24]) and ([3.6]-[24]) colors to K_s=18.4 and [3.6]=17.5, respectively. Star-forming and passive EROs in magnitude limited samples have different clustering properties and host dark halo masses, and cannot be simply understood as a single population. Based on the clustering, we find that bright passive EROs are the likely progenitors of >4L^* elliptical galaxies. Bright EROs with ongoing star formation were found to occupy denser environments than star-forming galaxies in the local Universe, making these the likely progenitors of >L^* local ellipticals. This suggests that the progenitors of massive >4L^* local ellipticals had stopped forming stars by z>1.2, but that the progenitors of less massive ellipticals (down to L^*) can still show significant star formation at this epoch.

preprint2012arXiv

The Importance of Broad Emission-Line Widths in Single Epoch Black Hole Mass Estimates

Estimates of the mass of super-massive black holes (BHs) in distant active galactic nuclei (AGNs) can be obtained efficiently only through single-epoch spectra, using a combination of their broad emission-line widths and continuum luminosities. Yet the reliability and accuracy of the method, and the resulting mass estimates, M_BH, remain uncertain. A recent study by Croom using a sample of SDSS, 2QZ and 2SLAQ quasars suggests that line widths contribute little information about the BH mass in these single-epoch estimates and can be replaced by a constant value without significant loss of accuracy. In this Letter, we use a sample of nearby reverberation-mapped AGNs to show that this conclusion is not universally applicable. We use the bulge luminosity (L_Bulge) of these local objects to test how well the known M_BH - L_Bulge correlation is recovered when using randomly assigned line widths instead of the measured ones to estimate M_BH. We find that line widths provide significant information about M_BH, and that for this sample, the line width information is just as significant as that provided by the continuum luminosities. We discuss the effects of observational biases upon the analysis of Croom and suggest that the results can probably be explained as a bias of flux-limited, shallow quasar samples.

preprint2012arXiv

The Mid-IR Contribution Of Dust Enshrouded Stars In Six Nearby Galaxies

We measure the integrated contributions of dusty AGB stars and other luminous red mid-IR sources to the mid-IR luminosities of 6 galaxies (M81, NGC 2403, NGC 300, M33 and the Magellanic Clouds). We find the dusty AGB stars whose mid-IR fluxes are dominated by dust rather than photospheric emission contribute from 0.6% (M81) to 5.6% (SMC) of the 3.6 micron flux and 1.0% (M81) to 10.1% (SMC) of the 4.5 micron flux. We find a trend of decreasing AGB contribution with increasing galaxy metallicity, luminosity and mass and decreasing SSFR. However, these galaxy properties are strongly correlated in our sample and the simplest explanation of the trend is galaxy metallicity. Bright, red sources other than dusty AGB stars represent a smaller fraction of the luminosity, ~1.2% at 3.6 microns, however their dust is likely cooler and their contributions are likely larger at longer wavelengths. Excluding the SMC, the contribution from these red sources correlates with the specific star formation rate as we would expect for massive stars. In total, after correcting for dust emission at other wavelengths, the dust around AGB stars radiates 0.1-0.8% of the bolometric luminosities of the galaxies. Thus, hot dust emission from AGB and other luminous dusty stars represent a small fraction of the total luminosities of the galaxies but a significant fraction of their mid-IR emissions.

preprint2012arXiv

The Structure of the Broad Line Region in AGN: I. Reconstructed Velocity-Delay Maps

We present velocity-resolved reverberation results for five active galactic nuclei. We recovered velocity-delay maps using the maximum-entropy method for four objects: Mrk 335, Mrk 1501, 3C120, and PG2130+099. For the fifth, Mrk 6, we were only able to measure mean time delays in different velocity bins of the Hbeta emission line. The four velocity-delay maps show unique dynamical signatures for each object. For 3C120, the Balmer lines show kinematic signatures consistent with both an inclined disk and infalling gas, but the HeII 4686 emission line is suggestive only of inflow. The Balmer lines in Mrk 335, Mrk 1501, and PG 2130+099 show signs of infalling gas, but the HeII emission in Mrk 335 is consistent with an inclined disk. We also see tentative evidence of combined virial motion and infalling gas from the velocity-binned analysis of Mrk 6. The maps for 3C120 and Mrk 335 are two of the most clearly defined velocity-delay maps to date. These maps constitute a large increase in the number of objects for which we have resolved velocity-delay maps and provide evidence supporting the reliability of reverberation-based black hole mass measurements.

preprint2012arXiv

The Unusual Temporal and Spectral Evolution of the Type IIn Supernova 2011ht

We present very early UV to optical photometric and spectroscopic observations of the peculiar Type IIn supernova (SN) 2011ht in UGC 5460. The UV observations of the rise to peak are only the second ever recorded for a Type IIn SN and are by far the most complete. The SN, first classified as a SN impostor, slowly rose to a peak of M_V \sim -17 in \sim55 days. In contrast to the \sim2 magnitude increase in the v-band light curve from the first observation until peak, the UV flux increased by >7 magnitudes. The optical spectra are dominated by strong, Balmer emission with narrow peaks (FWHM\sim600 km/s), very broad asymmetric wings (FWHM\sim4200 km/s), and blue shifted absorption (\sim300 km/s) superposed on a strong blue continuum. The UV spectra are dominated by FeII, MgII, SiII, and SiIII absorption lines broadened by \sim1500 km/s. Merged X-ray observations reveal a L_(0.2-10)=(1.0+/-0.2)x10^(39) erg/s. Some properties of SN 2011ht are similar to SN impostors, while others are comparable to Type IIn SNe. Early spectra showed features typical of luminous blue variables at maximum and during giant eruptions. However, the broad emission profiles coupled with the strong UV flux have not been observed in previous SN impostors. The absolute magnitude and energetics (~2.5x10^(49) ergs in the first 112 days) are reminiscent of normal Type IIn SN, but the spectra are of a dense wind. We suggest that the mechanism for creating this unusual profile could be a shock interacting with a shell of material that was ejected a year before the discovery of the SN.

preprint2012arXiv

Unmasking the Supernova Impostors

(ABRIDGED) The canonical picture of a supernova impostor is a -11 < M_V < -14 optical transient from a massive (M > 40Msun) star during which the star ejects a dense shell of material. Dust formed in the ejecta then obscures the star. In this picture, the geometric expansion of the shell leads to clear predictions for the evolution of the optical depths and hence the evolution of the optical through mid-IR emissions. Here we review the theory of this standard model and then examine the impostors SN1954J, SN1997bs, SN1999bw, SN2000ch, SN2001ac, SN2002bu, SN2002kg and SN2003gm, as well as the potential archetype eta Carinae. SN1999bw, SN2000ch, SN2001ac, SN2002bu and SN2003gm all show mid-IR emission indicative of dust, and the luminosities of SN1999bw, SN2001ac, SN2002bu and SN2003gm are dominated by dust emission. The properties of these sources are broadly inconsistent with the predictions of the canonical model. There are probably two classes of sources. In one class (eta Carinae, SN1954J, SN1997bs, and (maybe) SN2003gm), the optical transient is a signal that the star is entering a phase with very high mass loss rates that must last far longer than the visual transient. The second class (SN1999bw, SN2001ac, SN2002bu and (maybe) SN2003gm) has the different physics of SN2008S and the 2008 NGC300 transient, where they are obscured by dust re-forming in a pre-existing wind after it was destroyed by an explosive transient. There are no cases where the source at late times is significantly fainter than the progenitor star. All these dusty transients are occurring in relatively low mass (M < 25Msun) stars rather than high mass (M > 40Msun) stars radiating near the Eddington limit like eta Carinae. The durations and energetics of these transients cannot be properly characterized without near/mid-IR observations.

preprint2011arXiv

A New Microlensing Event in the Doubly-Imaged Quasar Q0957+561

We present evidence for ultraviolet/optical microlensing in the gravitationally lensed quasar Q0957+561. We combine new measurements from our optical monitoring campaign at the United States Naval Observatory, Flagstaff (USNO) with measurements from the literature and find that the time-delay-corrected r-band flux ratio m_A - m_B has increased by ~0.1 magnitudes over a period of five years beginning in the fall of 2005. We apply our Monte Carlo microlensing analysis procedure to the composite light curves, obtaining a measurement of the optical accretion disk size, log {(r_s/cm)[cos(i)/0.5]^{1/2}} = 16.2^{+0.5}_{-0.6}, that is consistent with the quasar accretion disk size - black hole mass relation.

preprint2011arXiv

A Reverberation Lag for the High-Ionization Component of the Broad Line Region in the Narrow-Line Seyfert 1 Mrk 335

We present the first results from a detailed analysis of photometric and spectrophotometric data on the narrow-line Seyfert 1 galaxy Mrk 335, collected over a 120-day span in the fall of 2010. From these data we measure the lag in the He II 4686 broad emission line relative to the optical continuum to be 2.7 \pm 0.6 days and the lag in the Hβ4861 broad emission line to be 13.9 \pm 0.9 days. Combined with the line width, the He II lag yields a black hole mass, MBH = (2.6 \pm 0.8)\times 10^7 Msun. This measurement is consistent with measurements made using the Hβ4861 line, suggesting that the He II emission originates in the same region as Hβ, but at a much smaller radius. This constitutes the first robust lag measurement for a high-ionization line in a narrow-line Seyfert 1 galaxy.

preprint2011arXiv

AGES: The AGN and Galaxy Evolution Survey

The AGN and Galaxy Evolution Survey (AGES) is a redshift survey covering, in its standard fields, 7.7 square degrees of the Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). The final sample consists of 23745 redshifts. There are well-defined galaxy samples in ten bands (the Bw, R, I, J, K, IRAC 3.6, 4.5, 5.8 and 8.0 micron and MIPS 24 micron bands) to a limiting magnitude of I<20 mag for spectroscopy. For these galaxies, we obtained 18163 redshifts from a sample of 35200 galaxies, where random sparse sampling was used to define statistically complete sub-samples in all ten photometric bands. The median galaxy redshift is 0.31, and 90% of the redshifts are in the range 0.085<z<0.66. AGN were selected as radio, X-ray, IRAC mid-IR and MIPS 24 micron sources to fainter limiting magnitudes (I<22.5 mag for point sources). Redshifts were obtained for 4764 quasars and galaxies with AGN signatures, with 2926, 1718, 605, 119 and 13 above redshifts of 0.5, 1, 2, 3 and 4, respectively. We detail all the AGES selection procedures and present the complete spectroscopic redshift catalogs, spectra, and spectral energy distribution decompositions. The photometric redshift estimates are for all sources in the AGES samples.

preprint2011arXiv

An Alternative Approach To Measuring Reverberation Lags in Active Galactic Nuclei

Motivated by recent progress in the statistical modeling of quasar variability, we develop a new approach to measuring emission-line reverberation lags to estimate the size of broad-line regions (BLRs) in active galactic nuclei. Assuming that all emission-line light curves are scaled, smoothed, and displaced versions of the continuum, this alternative approach fits the light curves directly using a damped random walk model and aligns them to recover the time lag and its statistical confidence limits. We introduce the mathematical formalism of this approach and demonstrate its ability to cope with some of the problems for traditional methods, such as irregular sampling, correlated errors, and seasonal gaps. We redetermine the lags for 87 emission lines in 31 quasars and reassess the BLR size--luminosity relationship using 60 H-beta lags. We confirm the general results from the traditional cross-correlation methods, with a few exceptions. Our method, however, also supports a broad range of extensions. In particular, it can simultaneously fit multiple lines and continuum light curves which improves the lag estimate for the lines and provides estimates of the error correlations between them. Determining these correlations is of particular importance for interpreting emission-line velocity--delay maps. We can also include parameters for luminosity-dependent lags or line responses. We use this to detect the scaling of the BLR size with continuum luminosity in NGC 5548.

preprint2011arXiv

Black Hole Mass Estimates Based on CIV are Consistent with Those Based on the Balmer Lines

Using a sample of high-redshift lensed quasars from the CASTLES project with observed-frame ultraviolet or optical and near-infrared spectra, we have searched for possible biases between supermassive black hole (BH) mass estimates based on the CIV, Halpha and Hbeta broad emission lines. Our sample is based upon that of Greene, Peng & Ludwig, expanded with new near-IR spectroscopic observations, consistently analyzed high S/N optical spectra, and consistent continuum luminosity estimates at 5100A. We find that BH mass estimates based on the FWHM of CIV show a systematic offset with respect to those obtained from the line dispersion, sigma_l, of the same emission line, but not with those obtained from the FWHM of Halpha and Hbeta. The magnitude of the offset depends on the treatment of the HeII and FeII emission blended with CIV, but there is little scatter for any fixed measurement prescription. While we otherwise find no systematic offsets between CIV and Balmer line mass estimates, we do find that the residuals between them are strongly correlated with the ratio of the UV and optical continuum luminosities. Removing this dependency reduces the scatter between the UV- and optical-based BH mass estimates by a factor of approximately 2, from roughly 0.35 to 0.18 dex. The dispersion is smallest when comparing the CIV sigma_l mass estimate, after removing the offset from the FWHM estimates, and either Balmer line mass estimate. The correlation with the continuum slope is likely due to a combination of reddening, host contamination and object-dependent SED shapes. When we add additional heterogeneous measurements from the literature, the results are unchanged.

preprint2011arXiv

Discovery of Variability of the Progenitor of SN 2011dh in M51 Using the Large Binocular Telescope

We show that the candidate progenitor of the core-collapse SN 2011dh in M51 (8 Mpc away) was fading by 0.039 +- 0.006 mag/year during the three years prior to the supernova, and that this level of variability is moderately unusual for other similar stars in M 51. While there are uncertainties about whether the true progenitor was a blue companion to this candidate, the result illustrates that there are no technical challenges to obtaining fairly high precision light curves of supernova progenitors using ground based observations of nearby (<10 Mpc) galaxies with wide field cameras on 8m-class telescopes. While other sources of variability may dominate, it is even possible to reach into the range of evolution rates required by the quasi-static evolution of the stellar envelope. For M 81, where we have many more epochs and a slightly longer time baseline, our formal 3 sigma sensitivity to slow changes is presently 3 millimag/year for a M_V ~= -8 mag star. In short, there is no observational barrier to determining whether the variability properties of stars in their last phases of evolution (post Carbon ignition) are different from earlier phases.

preprint2011arXiv

Dusty Explosions from Dusty Progenitors: The Physics of SN2008S and the 2008 NGC300-OT

SN2008S and the 2008 NGC300-OT were explosive transients of stars self-obscured by very dense, dusty stellar winds. An explosive transient with an un-observed shock break-out luminosity of order 10^10 Lsun is required to render the transients little obscured and visible in the optical at their peaks. Such a large break-out luminosity then implies that the progenitor stars were cool, red supergiants, most probably ~9 Msun extreme AGB (EAGB) stars. As the shocks generated by the explosions propagate outward through the dense wind, they produce a shock luminosity in soft X-rays that powers the long-lived luminosity of the transients. Unlike typical cases of transients exploding into a surrounding circumstellar medium, the progenitor winds in these systems are optically thick to soft X-rays, easily absorb radio emission and rapidly reform dust destroyed by the peak luminosity of the transients. As a result, X-rays are absorbed by the gas and the energy is ultimately radiated by the reformed dust. Three years post-peak, both systems are still significantly more luminous than their progenitor stars, but they are again fully shrouded by the re-formed dust and only visible in the mid-IR. The high luminosity and heavy obscuration may make it difficult to determine the survival of the progenitor stars for ~10 years. However, our model indicates that SN2008S, but not the NGC300-OT, should now be a detectable X-ray source. SN2008S has a higher estimated shock velocity and a lower density wind, so the X-rays begin to escape at a much earlier phase.

preprint2011arXiv

Object-X: The Brightest Mid-IR Point Source in M33

We discuss the nature of the brightest mid-IR point source (which we dub Object X) in the nearby galaxy M33. Although multi-wavelength data on this object have existed in the literature for some time, it has not previously been recognized as the most luminous mid-IR object in M33 because it is entirely unremarkable in both optical and near-IR light. In the Local Group Galaxies Survey, Object X is a faint red source visible in VRI and H-alpha but not U or B. It was easily seen at JHK_s in the 2MASS survey. It is the brightest point source in all four Spitzer IRAC bands and is also visible in the MIPS 24-micron band. Its bolometric luminosity is 5x10^5 L_sun. The source is optically variable on short time scales (tens of days) and is also slightly variable in the mid-IR, indicating that it is a star. Archival photographic plates (from 1949 and 1991) show no optical source, so the star has been obscured for at least half a century. Its properties are similar to those of the Galactic OH/IR star IRC+10420 which has a complex dusty circumstellar structure resulting from episodic low velocity mass ejections. We propose that Object X is a M>30 M_sun evolved star obscured in its own dust ejected during episodic mass loss events over at least half a century. It may emerge from its current ultra-short evolutionary phase as a hotter post-RSG star analogous to M33 Var A. The existence and rarity of such objects can be an important probe of a very brief yet eventful stellar evolutionary phase.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

The Astrophysical Implications of Dust Formation During The Eruptions of Hot, Massive Stars

Dust formation in the winds of hot stars is inextricably linked to the classic eruptive state of luminous blue variables (LBVs) because it requires very high mass loss rates, Mdot>10^(-2.5) Msun/year, for grains to grow and for the non-dust optical depth of the wind to shield the dust formation region from the true stellar photosphere. Thus, dusty shells around hot stars trace the history of "great" eruptions, and the statistics of such shells in the Galaxy indicate that these eruptions are likely the dominant mass loss mechanism for evolved, M>40Msun stars. Dust formation at such high Mdot also explains why very large grains (amax>1 micron) are frequently found in these shells, since amax \propto Mdot. The statistics of these shells (numbers, ages, masses, and grain properties such as amax) provide an archaeological record of this mass loss process. In particular, the velocities v, transient durations (where known) and ejected masses M of the Galactic shells and the supernova "impostors" proposed as their extragalactic counterparts are very different. While much of the difference is a selection effect created by shell lifetimes ~1/(v sqrt(M)), more complete Galactic and extragalactic surveys are needed to demonstrate that the two phenomena share a common origin given that their observed properties are essentially disjoint. If even small fractions (1%) of SNe show interactions with such dense shells of ejecta, as is currently believed, then the driving mechanism of the eruptions must be associated with the very final phases of stellar evolution, suggestive of some underlying nuclear burning instability.

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The Magellanic Quasars Survey. II. Confirmation of 144 New Active Galactic Nuclei Behind the Southern Edge of the Large Magellanic Cloud

We quadruple the number of quasars known behind the Large Magellanic Cloud (LMC) from 55 (42 in the LMC fields of the third phase of the Optical Gravitational Lensing Experiment (OGLE)) to 200 by spectroscopically confirming 169 (144 new) quasars from a sample of 845 observed candidates in four ~3 deg^2 Anglo-Australian Telescope/AAOmega fields south of the LMC center. The candidates were selected based on their Spitzer mid-infrared colors, X-ray emission, and/or optical variability properties in the database of the OGLE microlensing survey. The contaminating sources can be divided into 115 young stellar objects (YSOs), 17 planetary nebulae (PNe), 39 Be and 24 blue stars, 68 red stars, and 12 objects classed as either YSO/PN or blue star/YSO. There are also 402 targets with either featureless spectra or too low signal-to-noise ratio for source classification. Our quasar sample is 50% (30%) complete at I = 18.6 mag (19.3 mag). The newly discovered active galactic nuclei (AGNs) provide many additional reference points for proper motion studies of the LMC, and the sample includes 10 bright AGNs (I < 18 mag) potentially suitable for absorption line studies. Their primary use, however, is for detailed studies of quasar variability, as they all have long-term, high cadence, continuously growing light curves from the microlensing surveys of the LMC. Completing the existing Magellanic Quasars Survey fields in the LMC and Small Magellanic Cloud should yield a sample of ~700 well-monitored AGNs, and expanding it to the larger regions covered by the OGLE-IV survey should yield a sample of ~3600 AGNs.

preprint2011arXiv

The Mid-IR and X-ray Selected QSO Luminosity Function

We present the J-band luminosity function of 1838 mid-infrared and X-ray selected AGNs in the redshift range 0<z<5.85. These luminosity functions are constructed by combining the deep multi-wavelength broad-band observations from the UV to the mid-IR of the NDWFS Bootes field with the X-ray observations of the XBootes survey and the spectroscopic observations of the same field by AGES. Our sample is primarily composed of IRAC-selected AGNs, targeted using modifications of the Stern et al.(2005) criteria, complemented by MIPS 24 microns and X-ray selected AGNs to alleviate the biases of IRAC mid-IR selection against z~4.5 quasars and AGNs faint with respect to their hosts. This sample provides an accurate link between low and high redshift AGN luminosity functions and does not suffer from the usual incompleteness of optical samples at z~3. We find that the space density of the brightest quasars strongly decreases from z=3 to z=0, while the space density of faint quasars is at least flat, and possibly increasing, over the same redshift range. At z>3 we observe a decrease in the space density of quasars of all brightnesses. We model the luminosity function by a double power-law and find that its evolution cannot be described by either pure luminosity or pure density evolution, but must be a combination of both. Our best-fit model has bright and faint power-law indices consistent with the low redshift measurements based on the 2QZ and 2SLAQ surveys and it generally agrees with the number of bright quasars predicted by other LFs at all redshifts. If we construct the QSO luminosity function using only the IRAC-selected AGNs, we find that the biases inherent to this selection method significantly modify the behavior of phi*(z) only for z<1 and have no significant impact upon the characteristic magnitude M*_J(z).

preprint2010arXiv

Census of Self-Obscured Massive Stars in Nearby Galaxies with Spitzer: Implications for Understanding the Progenitors of SN 2008S-Like Transients

A new link in the causal mapping between massive stars and potentially fatal explosive transients opened with the 2008 discovery of the dust-obscured progenitors of the luminous outbursts in NGC 6946 and NGC 300. Here we carry out a systematic mid-IR photometric search for massive, luminous, self-obscured stars in four nearby galaxies: M33, NGC 300, M81, and NGC 6946. For detection, we use only the 3.6 micron and 4.5 micron IRAC bands, as these can still be used for multi-epoch Spitzer surveys of nearby galaxies (=<10 Mpc). We combine familiar PSF and aperture-photometry with an innovative application of image subtraction to catalog the self-obscured massive stars in these galaxies. In particular, we verify that stars analogous to the progenitors of the NGC 6946 (SN 2008S) and NGC 300 transients are truly rare in all four galaxies: their number may be as low as ~1 per galaxy at any given moment. This result empirically supports the idea that the dust-enshrouded phase is a very short-lived phenomenon in the lives of many massive stars and that these objects constitute a natural extension of the AGB sequence. We also provide mid-IR catalogs of sources in NGC 300, M81, and NGC 6946.

preprint2010arXiv

Low Resolution Spectral Templates for AGNs and Galaxies

We present a set of low resolution empirical SED templates for AGNs and galaxies in the wavelength range from 0.03 to 30 microns. These templates form a non-negative basis of the color space of such objects and have been derived from a combination 14448 galaxies and 5347 likely AGNs in the NDWFS Bootes field. We briefly describe how the templates are derived and discuss some applications of them. In particular, we discuss biases in commonly used AGN mid-IR color selection criteria and the expected distribution of sources in the current WISE satellite mission.

preprint2010arXiv

Low Resolution Spectral Templates For AGNs and Galaxies From 0.03 -- 30 microns

We present a set of low resolution empirical SED templates for AGNs and galaxies in the wavelength range from 0.03 to 30 microns based on the multi-wavelength photometric observations of the NOAO Deep-Wide Field Survey Bootes field and the spectroscopic observations of the AGN and Galaxy Evolution Survey. Our training sample is comprised of 14448 galaxies in the redshift range 0<~z<~1 and 5347 likely AGNs in the range 0<~z<~5.58. We use our templates to determine photometric redshifts for galaxies and AGNs. While they are relatively accurate for galaxies, their accuracies for AGNs are a strong function of the luminosity ratio between the AGN and galaxy components. Somewhat surprisingly, the relative luminosities of the AGN and its host are well determined even when the photometric redshift is significantly in error. We also use our templates to study the mid-IR AGN selection criteria developed by Stern et al.(2005) and Lacy et al.(2004). We find that the Stern et al.(2005) criteria suffers from significant incompleteness when there is a strong host galaxy component and at z =~ 4.5, when the broad Halpha emission line is redshifted into the [3.6] band, but that it is little contaminated by low and intermediate redshift galaxies. The Lacy et al.(2004) criterion is not affected by incompleteness at z =~ 4.5 and is somewhat less affected by strong galaxy host components, but is heavily contaminated by low redshift star forming galaxies. Finally, we use our templates to predict the color-color distribution of sources in the upcoming WISE mission and define a color criterion to select AGNs analogous to those developed for IRAC photometry. We estimate that in between 640,000 and 1,700,000 AGNs will be identified by these criteria, but will have serious completeness problems for z >~ 3.4.

preprint2010arXiv

Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk

We model the time variability of ~9,000 spectroscopically confirmed quasars in SDSS Stripe 82 as a damped random walk. Using 2.7 million photometric measurements collected over 10 years, we confirm the results of Kelly et al. (2009) and Kozłowski et al. (2010) that this model can explain quasar light curves at an impressive fidelity level (0.01-0.02 mag). The damped random walk model provides a simple, fast [O(N) for N data points], and powerful statistical description of quasar light curves by a characteristic time scale (tau) and an asymptotic rms variability on long time scales (SF_inf). We searched for correlations between these two variability parameters and physical parameters such as luminosity and black hole mass, and rest-frame wavelength. We find that tau increases with increasing wavelength with a power law index of 0.17, remains nearly constant with redshift and luminosity, and increases with increasing black hole mass with power law index of 0.21+/-0.07. The amplitude of variability is anti-correlated with the Eddington ratio, which suggests a scenario where optical fluctuations are tied to variations in the accretion rate. The radio-loudest quasars have systematically larger variability amplitudes by about 30%, when corrected for the other observed trends, while the distribution of their characteristic time scale is indistinguishable from that of the full sample. We do not detect any statistically robust differences in the characteristic time scale and variability amplitude between the full sample and the small subsample of quasars detected by ROSAT. Our results provide a simple quantitative framework for generating mock quasar light curves, such as currently used in LSST image simulations. (abridged)

preprint2010arXiv

Quasar Selection Based on Photometric Variability

We develop a method for separating quasars from other variable point sources using SDSS Stripe 82 light curve data for ~10,000 variable objects. To statistically describe quasar variability, we use a damped random walk model parametrized by a damping time scale, tau, and an asymptotic amplitude (structure function), SF_inf. With the aid of an SDSS spectroscopically confirmed quasar sample, we demonstrate that variability selection in typical extragalactic fields with low stellar density can deliver complete samples with reasonable purity (or efficiency, E). Compared to a selection method based solely on the slope of the structure function, the inclusion of the tau information boosts E from 60% to 75% while maintaining a highly complete sample (98%) even in the absence of color information. For a completeness of C=90%, E is boosted from 80% to 85%. Conversely, C improves from 90% to 97% while maintaining E=80% when imposing a lower limit on tau. With the aid of color selection, the purity can be further boosted to 96%, with C= 93%. Hence, selection methods based on variability will play an important role in the selection of quasars with data provided by upcoming large sky surveys, such as Pan-STARRS and the Large Synoptic Survey Telescope (LSST). For a typical (simulated) LSST cadence over 10 years and a photometric accuracy of 0.03 mag (achieved at i~22), C is expected to be 88% for a simple sample selection criterion of tau>100 days. In summary, given an adequate survey cadence, photometric variability provides an even better method than color selection for separating quasars from stars.

preprint2010arXiv

SDWFS-MT-1: A Self-Obscured Luminous Supernova at z~0.2

We report the discovery of a six-month-long mid-infrared transient, SDWFS-MT-1 (aka SN 2007va), in the Spitzer Deep, Wide-Field Survey of the NOAO Deep Wide-Field Survey Bootes field. The transient, located in a z=0.19 low luminosity (M_[4.5]~-18.6 mag, L/L_MilkyWay~0.01) metal-poor (12+log(O/H)~7.8) irregular galaxy, peaked at a mid-infrared absolute magnitude of M_[4.5]~-24.2 in the 4.5 micron Spitzer/IRAC band and emitted a total energy of at least 10^51 ergs. The optical emission was likely fainter than the mid-infrared, although our constraints on the optical emission are poor because the transient peaked when the source was "behind" the Sun. The Spitzer data are consistent with emission by a modified black body with a temperature of ~1350 K. We rule out a number of scenarios for the origin of the transient such as a Galactic star, AGN activity, GRB, tidal disruption of a star by a black hole and gravitational lensing. The most plausible scenario is a supernova exploding inside a massive, optically thick circumstellar medium, composed of multiple shells of previously ejected material. If the proposed scenario is correct, then a significant fraction (~10%) of the most luminous supernova may be self-enshrouded by dust not only before but also after the supernova occurs. The spectral energy distribution of the progenitor of such a supernova would be a slightly cooler version of eta Carina, peaking at 20-30 microns.

preprint2010arXiv

The end of the white dwarf cooling sequence in M 67

In this paper, we present for the first time a proper-motion-selected white dwarf (WD) sample of the old Galactic open cluster M 67, down to the bottom of the WD cooling sequence (CS). The color-magnitude diagram is based on data collected with the LBC-Blue camera at the prime focus of LBT. As first epoch data, we used CFHT-archive images collected 10 years before LBC data. We measured proper motions of all the identified sources. Proper motions are then used to separate foreground and background objects from the cluster stars, including WDs. Finally, the field-object cleaned WD CS in the V vs. B-I color-magnitude diagram is compared with the models. We confirm that the age derived from the location of the bottom of the WD CS is consistent with the turn off age.

preprint2010arXiv

The SN 2008S Progenitor Star: Gone or Again Self-Obscured?

We obtained late-time optical and near-IR imaging of SN 2008S with the Large Binocular Telescope (LBT). We find that (1) it is again invisible at optical (UBVR) wavelengths to magnitude limits of approximately 25 mag, and (2) while detected in the near-IR (HK) at approximately 20 mag, it is fading rapidly. The near-IR detections in March and May 2010 are consistent with dust emission at a blackbody temperature of T ~ 900 K and a total luminosity of L ~ 40000 Lsun, comparable to the luminosity of the obscured progenitor star. If it is a supernova, the near-IR emission is likely due to shock heated dust since the elapsed time from peak is too long to support a near-IR dust echo and the decline in luminosity is shallower than the 56Co slope. If it is reprocessed emission from a surviving progenitor, a dust photosphere must have reestablished itself closer to the star than before the transient (~40 AU rather than 150 AU), unless there is a second, cooler dust component that dominates at mid-IR wavelengths. The continued rapid fading at roughly constant temperature favors transient emission, but the SED peaks in the mid-IR and future Spitzer observations will be needed to close the case.

preprint2010arXiv

The Supernova Impostor Impostor SN 1961V: Spitzer Shows That Zwicky Was Right (Again)

SN 1961V, one of Zwicky's defining Type V supernovae (SN), was a peculiar transient in NGC 1058 that has variously been categorized as either a true core collapse SN leaving a black hole (BH) or neutron star (NS) remnant, or an eruption of a luminous blue variable (LBV) star. The former case is suggested by its association with a decaying non-thermal radio source, while the latter is suggested by its peculiar transient light curve and its low initial expansion velocities. The crucial difference is that the star survives a transient eruption but not an SN. All stars identified as possible survivors are significantly fainter, L_opt ~ 10^5 Lsun, than the L_opt ~ 3 10^6 Lsun progenitor star at optical wavelengths. While this can be explained by dust absorption in a shell of material ejected during the transient, the survivor must then be present as a L_IR ~ 3 10^6 Lsun mid-infrared source. Using archival Spitzer observations of the region, we show that such a luminous mid-IR source is not present. The brightest source of dust emission is only L_IR ~ 10^5 Lsun and does not correspond to the previously identified candidates for the surviving star. The dust cannot be made sufficiently distant and cold to avoid detection unless the ejection energy, mass and velocity scales are those of a SN or greater. We conclude that SN 1961V was a peculiar, but real, supernova. Its peculiarities are probably due to enhanced mass loss just prior to the SN, followed by the interactions of the SN blast wave with this ejecta. This adds to the evidence that there is a population of SN progenitors that have major mass loss episodes shortly before core collapse. The progenitor is a low metallicity, ~1/3 solar, high mass, M_ZAMS > 80 Msun, star, which means either that BH formation can be accompanied by an SN or that surprisingly high mass stars can form a NS.

preprint2009arXiv

Go Long, Go Deep: Finding Optical Jet Breaks for Swift-Era GRBs with the LBT

Using the 8.4m Large Binocular Telescope, we observed six GRB afterglows from 2.8 hours to 30.8 days after the burst triggers to systematically probe the late time behaviors of afterglows including jet breaks, flares, and supernova bumps. We detected five afterglows with Sloan r' magnitudes ranging from 23.0-26.3 mag. The depth of our observations allows us to extend the temporal baseline for measuring jet breaks by another decade in time scale. We detected two jet breaks and a third candidate, all of which are not detectable without deep, late time optical observations. In the other three cases, we do not detect the jet breaks either because of contamination from the host galaxy light, the presence of a supernova bump, or the intrinsic faintness of the optical afterglow. This suggests that the basic picture that GRBs are collimated is still valid and that the apparent lack of Swift jet breaks is due to poorly sampled afterglow light curves, particularly at late times.

preprint2009arXiv

Mid-Infrared Galaxy Luminosity Functions from the AGN and Galaxy Evolution Survey

We present galaxy luminosity functions at 3.6, 4.5, 5.8, and 8.0 micron measured by combining photometry from the IRAC Shallow Survey with redshifts from the AGN and Galaxy Evolution Survey of the NOAO Deep Wide-Field Survey Bootes field. The well-defined IRAC samples contain 3800-5800 galaxies for the 3.6-8.0 micron bands with spectroscopic redshifts and z < 0.6. We obtained relatively complete luminosity functions in the local redshift bin of z < 0.2 for all four IRAC channels that are well fit by Schechter functions. We found significant evolution in the luminosity functions for all four IRAC channels that can be fit as an evolution in M* with redshift, ΔM* = Qz. While we measured Q=1.2\pm0.4 and 1.1\pm0.4 in the 3.6 and 4.5 micron bands consistent with the predictions from a passively evolving population, we obtained Q=1.8\pm1.1 in the 8.0 micron band consistent with other evolving star formation rate estimates. We compared our LFs with the predictions of semi-analytical galaxy formation and found the best agreement at 3.6 and 4.5 micron, rough agreement at 8.0 micron, and a large mismatch at 5.8 micron. These models also predicted a comparable Q value to our luminosity functions at 8.0 micron, but predicted smaller values at 3.6 and 4.5 micron. We also measured the luminosity functions separately for early and late-type galaxies. While the luminosity functions of late-type galaxies resemble those for the total population, the luminosity functions of early-type galaxies in the 3.6 and 4.5 micron bands indicate deviations from the passive evolution model, especially from the measured flat luminosity density evolution. Combining our estimates with other measurements in the literature, we found (53\pm18)% of the present stellar mass of early-type galaxies has been assembled at z=0.7.

preprint2009arXiv

Stellar Binary Companions to Supernova Progenitors

For typical models of binary statistics, 50-70% of core-collapse supernova (ccSN) progenitors are members of a stellar binary at the time of the explosion. Independent of any consequences of mass transfer, this has observational consequences that can be used to study the binary properties of massive stars. In particular, the secondary companion to the progenitor of a Type Ib/c SN is frequently (~50%) the more optically luminous star since the high effective temperatures of the stripped progenitors make it relatively easy for a lower luminosity, cooler secondary to emit more optical light. Secondaries to the lower mass progenitors of Type II SN will frequently produce excess blue emission relative to the spectral energy distribution of the red primary. Available data constrain the models weakly. Any detected secondaries also provide an independent lower bound on the progenitor mass and, for historical SN, show that it was not a Type Ia event. Bright ccSN secondaries have an unambiguous, post-explosion observational signature - strong, blue-shifted, relatively broad absorption lines created by the developing supernova remnant. These can be used to locate historical SN with bright secondaries, confirm that a source is a secondary, and, potentially, measure abundances of ccSN ejecta. Luminous, hot secondaries will reionize the SNR on time scales of 100-1000 yrs that are faster than reionization by the reverse shock, creating peculiar HII regions due to the high metallicity and velocities of the ejecta.

preprint2009arXiv

The Sizes of the X-ray and Optical Emission Regions of RXJ1131-1231

We use gravitational microlensing of the four images of the z=0.658 quasar RXJ1131-1231 to measure the sizes of the optical and X-ray emission regions of the quasar. The (face-on) scale length of the optical disk at rest frame 400 nm is 1.3 10^15cm, while the half-light radius of the rest frame 0.3-17 keV X-ray emission is 2.3 10^14cm. The formal uncertainties are factors of 1.6 and 2.0, respectively. With the exception of the lower limit on the X-ray size, the results are very stable against any changes in the priors used in the analysis. Based on the Hbeta line-width, we estimate that the black hole mass is ~10^8 Msun, which corresponds to a gravitational radius of r_g~2 10^13 cm. Thus, the X-ray emission is emerging on scales of ~10r_g and the 400 nm emission on scales of ~70 r_g. A standard thin disk of this size should be significantly brighter than observed. Possible solutions are to have a flatter temperature profile or to scatter a large fraction of the optical flux on larger scales after it is emitted. While our calculations were not optimized to constrain the dark matter fraction in the lens galaxy, dark matter dominated models are favored. With well-sampled optical and X-ray light curves over a broad range of frequencies there will be no difficulty in extending our analysis to completely map the structure of the accretion disk as a function of wavelength.

preprint2008arXiv

A Survey About Nothing: Monitoring a Million Supergiants for Failed Supernovae

Extragalactic transient searches have historically been limited to looking for the appearance of new sources such as supernovae. It is now possible to carry out a new kind of survey that will do the opposite, that is, search for the disappearance of massive stars. This will entail the systematic observation of galaxies within a distance of 10 Mpc in order to watch ~10^6 supergiants. Reaching this critical number ensures that something will occur yearly, since these massive stars must end their lives with a core collapse within ~10^6 years. Using deep imaging and image subtraction it is possible to determine the fates of these stars whether they end with a bang (supernova) or a whimper (fall out of sight). Such a survey would place completely new limits on the total rate of all core collapses, which is critical for determining the validity of supernova models. It would also determine the properties of supernova progenitors, better characterize poorly understood optical transients, such as eta Carina-like mass ejections, find and characterize large numbers of Cepheids, luminous blue variables and eclipsing binaries, and allow the discovery of any new phenomena that inhabit this relatively unexplored parameter space.

preprint2008arXiv

Clusters of Galaxies in the First Half of the Universe from the IRAC Shallow Survey

We have identified 335 galaxy cluster and group candidates, 106 of which are at z > 1, using a 4.5 um selected sample of objects from a 7.25 deg^2 region in the Spitzer Infrared Array Camera (IRAC) Shallow Survey. Clusters were identified as 3-dimensional overdensities using a wavelet algorithm, based on photometric redshift probability distributions derived from IRAC and NOAO Deep Wide-Field Survey data. We estimate only ~10% of the detections are spurious. To date 12 of the z > 1 candidates have been confirmed spectroscopically, at redshifts from 1.06 to 1.41. Velocity dispersions of ~750 km/s for two of these argue for total cluster masses well above 10^14 M_sun, as does the mass estimated from the rest frame near infrared stellar luminosity. Although not selected to contain a red sequence, some evidence for red sequences is present in the spectroscopically confirmed clusters, and brighter galaxies are systematically redder than the mean galaxy color in clusters at all redshifts. The mean I - [3.6] color for cluster galaxies up to z ~ 1 is well matched by a passively evolving model in which stars are formed in a 0.1 Gyr burst starting at redshift z_f = 3. At z > 1, a wider range of formation histories is needed, but higher formation redshifts (i.e. z_f > 3) are favored for most clusters.

preprint2008arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses VII. Time delays and the Hubble constant from WFI J2033-4723

Gravitationally lensed quasars can be used to map the mass distribution in lensing galaxies and to estimate the Hubble constant H0 by measuring the time delays between the quasar images. Here we report the measurement of two independent time delays in the quadruply imaged quasar WFI J2033-4723 (z = 1.66). Our data consist of R-band images obtained with the Swiss 1.2 m EULER telescope located at La Silla and with the 1.3 m SMARTS telescope located at Cerro Tololo. The light curves have 218 independent epochs spanning 3 full years of monitoring between March 2004 and May 2007, with a mean temporal sampling of one observation every 4th day. We measure the time delays using three different techniques, and we obtain Dt(B-A) = 35.5 +- 1.4 days (3.8%) and Dt(B-C) = 62.6 +4.1/-2.3 days (+6.5%/-3.7%), where A is a composite of the close, merging image pair. After correcting for the time delays, we find R-band flux ratios of F_A/F_B = 2.88 +- 0.04, F_A/F_C = 3.38 +- 0.06, and F_A1/F_A2 = 1.37 +- 0.05 with no evidence for microlensing variability over a time scale of three years. However, these flux ratios do not agree with those measured in the quasar emission lines, suggesting that longer term microlensing is present. Our estimate of H0 agrees with the concordance value: non-parametric modeling of the lensing galaxy predicts H0 = 67 +13/-10 km s-1 Mpc-1, while the Single Isothermal Sphere model yields H0 = 63 +7/-3 km s-1 Mpc-1 (68% confidence level). More complex lens models using a composite de Vaucouleurs plus NFW galaxy mass profile show twisting of the mass isocontours in the lensing galaxy, as do the non-parametric models. As all models also require a significant external shear, this suggests that the lens is a member of the group of galaxies seen in field of view of WFI J2033-4723.

preprint2008arXiv

Discovery of the Dust-Enshrouded Progenitor of SN 2008S with Spitzer

We report the discovery of the progenitor of the recent type IIn SN 2008S in the nearby galaxy NGC 6946. Surprisingly, it was not found in deep, pre-explosion optical images of its host galaxy taken with the Large Binocular Telescope, but only through examination of archival Spitzer mid-IR data. A source coincident with the SN 2008S position is clearly detected in the 4.5, 5.8, and 8.0 micron IRAC bands, showing no evident variability in the three years prior to the explosion, yet is undetected at 3.6 and 24 micron. The distinct presence of ~440 K dust, along with stringent LBT limits on the optical fluxes, suggests that the progenitor of SN 2008S was engulfed in a shroud of its own dust. The inferred luminosity of 3.5x10^4 Lsun implies a modest mass of ~10 Msun. We conclude that objects like SN 2008S are not exclusively associated with the deaths or outbursts of very massive eta Carinae-like objects. This conclusion holds based solely on the optical flux limits even if our identification of the progenitor with the mid-IR source is incorrect.

preprint2008arXiv

X-ray Microlensing in RXJ1131-1231 and HE1104-1805

We present results from a monitoring campaign performed with the Chandra X-ray Observatory of the gravitationally lensed quasars RX J1131-1231 and HE 1104-1805. We detect significant X-ray variability in all images of both quasars. The flux variability detected in image A of RX J1131-1231 is of particular interest because of its high amplitude (a factor of ~ 20). We interpret it as arising from microlensing since the variability is uncorrelated with that of the other images and the X-ray flux ratios show larger changes than the optical as we would expect for microlensing of the more compact X-ray emission regions. The differences between the X-ray and optical flux ratios of HE 1104-1805 are less dramatic, but there is no significant soft X-ray or dust absorption, implying the presence of X-ray microlensing in this system as well. Combining the X-ray data with the optical light curves we find that the X-ray emitting region of HE 1104-1805 is compact with a half-light radius ~ 6 r_g, where the gravitational radius is r_g = 3.6 x 10^14 cm, thus placing significant constraints on AGN corona models. We also find that the microlensing in HE 1104-1805 favors mass models for the lens galaxy that are dominated by dark matter. Finally, we better characterize the massive foreground cluster near RX J1131-1231, set limits on other sources of extended X-ray emission, and limit the fluxes of any central odd images to be 30-50 (3 sigma) times fainter than the observed images.

preprint2007arXiv

LBT Discovery of a Yellow Supergiant Eclipsing Binary in the Dwarf Galaxy Holmberg IX

In a variability survey of M81 using the Large Binocular Telescope we have discovered a peculiar eclipsing binary (MV ~ -7.1) in the field of the dwarf galaxy Holmberg IX. It has a period of 272 days and the light curve is well-fit by an overcontact model in which both stars are overflowing their Roche lobes. It is composed by two yellow supergiants (V-I ~ 1 mag, T_eff = 4800 K), rather than the far more common red or blue supergiants. Such systems must be rare. While we failed to find any similar systems in the literature, we did, however note a second example. The SMC F0 supergiant R47 is a bright (MV ~ -7.5) periodic variable whose All Sky Automated Survey (ASAS) light curve is well-fit as a contact binary with a 181 day period. We propose that these systems are the progenitors of supernovae like SN 2004et and SN 2006ov, which appeared to have yellow progenitors. The binary interactions (mass transfer, mass loss) limit the size of the supergiant to give it a higher surface temperature than an isolated star at the same core evolutionary stage. We also discuss the possibility of this variable being a long-period Cepheid.

preprint2007arXiv

Microlensing of the Lensed Quasar SDSS0924+0219

We analyze V, I and H band HST images and two seasons of R-band monitoring data for the gravitationally lensed quasar SDSS0924+0219. We clearly see that image D is a point-source image of the quasar at the center of its host galaxy. We can easily track the host galaxy of the quasar close to image D because microlensing has provided a natural coronograph that suppresses the flux of the quasar image by roughly an order of magnitude. We observe low amplitude, uncorrelated variability between the four quasar images due to microlensing, but no correlated variations that could be used to measure a time delay. Monte Carlo models of the microlensing variability provide estimates of the mean stellar mass in the lens galaxy (0.02 Msun < M < 1.0 Msun), the accretion disk size (the disk temperature is 5 x 10^4 K at 3.0 x 10^14 cm < rs < 1.4 x 10^15 cm), and the black hole mass (2.0 x 10^7 Msun < MBH η_{0.1}^{-1/2} (L/LE)^{1/2} < 3.3 x 10^8 Msun), all at 68% confidence. The black hole mass estimate based on microlensing is consistent with an estimate of MBH = 7.3 +- 2.4 x 10^7 Msun from the MgII emission line width. If we extrapolate the best-fitting light curve models into the future, we expect the the flux of images A and B to remain relatively stable and images C and D to brighten. In particular, we estimate that image D has a roughly 12% probability of brightening by a factor of two during the next year and a 45% probability of brightening by an order of magnitude over the next decade.

preprint2007arXiv

Mid-IR Observations and a Revised Time Delay for the Gravitational Lens System Quasar HE 1104-1805

The mid-IR flux ratios F_A/F_B = 2.84 +/- 0.06 of the two images of the gravitationally lensed quasar HE 1104-1805 show no wavelength dependence to within 3% across 3.6-8.0 um, no time dependence over 6 months and agree with the broad emission line flux ratios. This indicates that the mid-IR emission likely comes from scales large enough to be little affected by microlensing and that there is little differential extinction between the images. We measure a revised time-delay between these two images of 152.2 +2.8-3.0 days from R and V-band data covering 1997 to 2006. This time-delay indicates that the lens has an approximately flat rotation curve over scales of 1-2 R_e. We also observed uncorrelated variations of ~0.05 mag/yr which we attribute to microlensing of the optical emission from the accretion disk. The optical colors have also changed significantly in the sense that image A is now redder than image B, rather than bluer as it was in 1993.

preprint2007arXiv

The Quasar Accretion Disk Size - Black Hole Mass Relation

We use the microlensing variability observed for nine gravitationally lensed quasars to show that the accretion disk size at 2500 Angstroms is related to the black hole mass by log(R_2500/cm) = (15.6+-0.2) + (0.54+-0.28)log(M_BH/10^9M_sun). This scaling is consistent with the expectation from thin disk theory (R ~ M_BH^(2/3)), but it implies that black holes radiate with relatively low efficiency, log(eta) = -1.29+-0.44 + log(L/L_E) where eta=L/(Mdot c^2). These sizes are also larger, by a factor of ~3, than the size needed to produce the observed 0.8 micron quasar flux by thermal radiation from a thin disk with the same T ~ R^(-3/4) temperature profile. More sophisticated disk models are clearly required, particularly as our continuing observations improve the precision of the measurements and yield estimates of the scaling with wavelength and accretion rate.

preprint2006arXiv

Probing the Coevolution of Supermassive Black Holes and Galaxies Using Gravitationally Lensed Quasar Hosts

In the present-day universe, supermassive black hole masses (MBH) appear to be strongly correlated with their galaxy's bulge luminosity, among other properties. In this study, we explore the analogous relationship between MBH, derived using the virial method, and the stellar R-band bulge luminosity (Lr) or stellar bulge mass (M*) at epochs of 1 < z < 4.5 using a sample of 31 gravitationally lensed AGNs and 20 non-lensed AGNs. At redshifts z > 1.7 (10--12 Gyrs ago), we find that the observed MBH--Lr relation is nearly the same (to within ~0.3 mag) as it is today. When the observed Lr are corrected for luminosity evolution, this means that the black holes grew in mass faster than their hosts, with the MBH/M* mass ratio being a factor of > 4(+2)(-1) times larger at z > 1.7 than it is today. By the redshift range 1<z<1.7 (8-10 Gyrs ago), the MBH/M* ratio is at most two times higher than today, but it may be consistent with no evolution. Combining the results, we conclude that the ratio MBH/M* rises with look-back time, although it may saturate at ~6 times the local value. Scenarios in which moderately luminous quasar hosts at z>1.7 were fully formed bulges that passively faded to the present epoch are ruled out.

preprint2000arXiv

The Fundamental Plane of Gravitational Lens Galaxies and The Evolution of Early-Type Galaxies in Low Density Environments

Most gravitational lenses are early-type galaxies in relatively low density environments -- a ``field'' rather than a ``cluster'' population. We show that field early-type galaxies with 0 < z < 1, as represented by the lens galaxies, lie on the same fundamental plane as those in rich clusters at similar redshifts. We then use the fundamental plane to measure the combined evolutionary and K-corrections for early-type galaxies in the V, I and H bands. Only for passively evolving stellar populations formed at z > 2 (H_0=65 km/s Mpc, Omega_0=0.3, Lambda_0=0.7) can the lens galaxies be matched to the local fundamental plane. The high formation epoch and the lack of significant differences between the field and cluster populations contradict many current models of the formation history of early-type galaxies. Lens galaxy colors and the fundamental plane provide good photometric redshift estimates with an empirical accuracy of -0.03 +/- 0.11 for the 17 lenses with known redshifts. A mass model dominated by dark matter is more consistent with the data than either an isotropic or radially anisotropic constant M/L mass model, and a radially anisotropic model is better than an isotropic model.

preprint1997arXiv

The geometry of the quadruply imaged quasar PG 1115+080; implications for Ho

Time delay measurements have recently been reported for the lensed quasar PG 1115+080. These measurements can be used to derive Ho, but only if we can constrain the lensing potential. We have applied a recently developed deconvolution technique to analyse sub-arcsecond I band images of PG 1115+080, obtained at the Nordic Optical Telescope (NOT) and the Canada France Hawaii Telescope (CFHT). The high performance of the deconvolution code allows us to derive precise positions and magnitudes for the four lensed images of the quasar, as well as for the lensing galaxy. The new measurement of the galaxy position improves its precision by a factor of 3 and thus strengthens the constraints on the lensing potential. With the new data, a range of models incorporating some of the plausible systematic uncertainties yields Ho = 53 (+10/-7) km/s/mpc.

preprint1997arXiv

The Optical Properties of Gravitational Lens Galaxies as a Probe of Galaxy Structure and Evolution

We combine photometry and lens modeling to study the properties of 17 gravitational lens galaxies between z=0.1 and 1. Most of the lens galaxies are passively evolving early-type galaxies, with a few spirals. The colors, scale lengths, and ellipticities of lens galaxies are similar to those of the general population of early-type galaxies, although there may be a deficit of apparently round lens galaxies produced by the inclination dependence of lensing cross sections. The projected mass distributions are aligned with the projected light distributions to <~ 10 deg, except in the presence of a strong external tidal perturbation, suggesting that dark matter halos have orbits that are significantly modified by interactions with the baryonic component and are not far out of alignment with the stars. Lens galaxies obey image separation/lens luminosity correlations analogous to the Faber-Jackson and Tully-Fisher relations, which are consistent with standard dark matter lens models. The lens galaxy mass-to-light ratios decrease with redshift as d(log M/L_B)/dz = -0.3\pm0.1 (-0.5\pm0.1) for Omega_0=1 (0.1), thus providing direct evidence of passive evolution for a sample of early-type galaxies in low-density environments. The evolution-corrected mass-to-light ratios are generally larger than predicted by constant M/L dynamical models, although there is significant scatter; with improved photometry, lens galaxy mass-to-light ratios would better distinguish between constant M/L and dark matter models. These conclusions are limited primarily by the quality of lens galaxy photometry.

preprint1995arXiv

Deprojection of Axially Symmetric Objects

The deprojection of axisymmetric density distributions is generally indeterminate to within the addition of certain axisymmetric distributions (konus densities) that are invisible in projection. The known class of konus densities is expanded considerably here through the introduction of {\em semikonus functions}. These functions are closed with respect to multiplication in ordinary space, and the real part of an arbitrary polynomial of semikonus functions is a konus function. This property facilitates the construction of semikonus (and konus) functions with tailored properties, such as asymptotic forms. We also develop a simple technique for constructing several classes of konus distributions with arbitrary density profiles in the equatorial plane.

preprint1995arXiv

Gravitational Lenses and the Structure of Galaxies

Nearly singular isothermal mass distributions with small core radii are consistent with stellar dynamics, lens statistics, and lens models as a model for E/S0 galaxies. Models like the de Vaucouleurs model with a constant mass-to-light ratio are not. While the isothermal distributions are probably an oversimplification, E/S0 galaxies (at least in projection) must have significant amounts of dark matter on scales of an effective radius.

preprint1995arXiv

The Mass Distribution of the Lens Galaxy in MG1131+0456

We present the results of modeling 5, 8 and 15 GHz maps of the gravitational lens MG1131+0456 using the LensClean algorithm. Two models for the mass distribution in the lens were fit: a de Vaucouleurs model and a model with the two-dimensional potential phi = b^{2-alpha}(r^2+s^2)^{alpha/2}/2alpha, both in an external shear field. The best fit de Vaucouleurs model has an effective radius of R_e = 0.83+-0.13 arcseconds, larger than that measured at optical wavelengths. The best fit ``alpha model'' has alpha=0.6+-0.2 and a core radius of s = 0.19 +- 0.07 arcseconds. An isothermal model (alpha=1) is inconsistent with the data. The best alpha model fits the 8 GHz map significantly better than the best de Vaucouleurs model. Although it is difficult to distinguish visually the original image and the reconstructions, none of the models is statistically consistent with the data. We believe the primary problem is that an external shear is an inadequate model for the angular structure of the lens galaxy. The finite core radius in the models is required by the structure of the extended ring, so we confirm that the central component in the 8 GHz map is a central, lensed image of the source rather than the lens galaxy. The model lens positions, the position angle of the shear, and the mass interior to the ring are determined very precisely by the lens models. The time delay between the compact components is predicted with an uncertainty we estimate to be at most +- 9%; the formal uncertainty is +- 4%.