Cosmic rays (CRs) are a plausible mechanism for launching winds of cool material from the discs of star-forming galaxies. However, there is no consensus on what types of galaxies likely host CR-driven winds, or what role these winds might play in regulating galaxies' star formation rates. Using a detailed treatment of the transport and losses of hadronic CRs developed in the previous paper in this series, here we develop a semi-analytic model that allows us to assess the viability of using CRs to launch cool winds from galactic discs. In particular, we determine the critical CR fluxes -- and corresponding star formation rate surface densities -- above which hydrostatic equilibrium within a given galaxy is precluded because CRs drive the gas off in a wind or otherwise render it unstable. We show that, for star-forming galaxies with lower gas surface densities typical of the Galaxy and local dwarfs, the locus of this CR stability curve patrols the high side of the observed distribution of galaxies in the Kennicutt-Schmidt parameter space of star formation rate versus gas surface density. However, hadronic losses render CRs unable to drive winds in galaxies with higher surface densities. Our results show that quiescent, low surface density galaxies like the Milky Way are poised on the cusp of instability, such that small changes to ISM parameters can lead to the launching of CR-driven outflows, and we suggest that, as a result, CR feedback sets an ultimate limit to the star formation efficiency of most modern galaxies.

ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of $M_V \approx -20$ mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am declines extremely rapidly for a Type II SN, with a photospheric-phase decline rate of $\sim6.0$ mag (100d)$^{-1}$. Owing to the weakening of HI and the appearance of HeI in its later phases, ASASSN-18am is spectroscopically a Type IIb SN with a partially stripped envelope. However, its photometric and spectroscopic evolution show significant differences from typical SNe IIb. Using a radiative diffusion model, we find that the light curve requires a high synthesised $^{56}$Ni mass of $M_{\rm Ni} \sim0.4\rm M_\odot$ and ejecta with high kinetic energy $E_{\rm kin}= (5-9)\times10^{51} $ erg. Introducing a magnetar central engine still requires $M_{\rm Ni} \sim0.3\, M_\odot$ and $E_{\rm kin}= 3\times10^{51} $ erg. The high $^{56}$Ni mass is consistent with strong iron-group nebular lines in its spectra, which are also similar to several SNe Ic-BL with high $^{56}$Ni yields. The earliest spectrum shows "flash ionization" features, from which we estimate a mass-loss rate of $ \dot{M}\approx 2\times10^{-4} \, \rm M_\odot~yr^{-1} $. This wind density is too low to power the luminous light curve by ejecta-CSM interaction. We measure expansion velocities as high as $ 17,000 \rm~km~s^{-1}$ for H$\alpha$, which is remarkably high compared to other SNe II. We estimate an oxygen core mass of 1.7-3.1 $\rm M_\odot$ using the [OI] luminosity measured from a nebular-phase spectrum, implying a progenitor with a zero-age main sequence mass of 19-24 $\rm M_\odot$.

In the Milky Way, cosmic rays (CRs) are dynamically important in the interstellar medium, contribute to hydrostatic balance, and may help regulate star formation. However, we know far less about the importance of CRs in galaxies whose gas content or star formation rate differ significantly from those of the Milky Way. Here we construct self-consistent models for hadronic CR transport, losses, and contribution to pressure balance as a function of galaxy properties, covering a broad range of parameters from dwarfs to extreme starbursts. While the CR energy density increases from $\sim 1$ eV cm$^{-3}$ to $\sim 1$ keV cm$^{-3}$ over the range from sub-Milky Way dwarfs to bright starbursts, strong hadronic losses render CRs increasingly unimportant dynamically as the star formation rate surface density increases. In Milky Way-like systems, CR pressure is typically comparable to turbulent gas and magnetic pressure at the galactic midplane, but the ratio of CR pressure to gas pressure drops to $\sim 10^{-3}$ in dense starbursts. Galaxies also become increasingly CR calorimetric and gamma-ray bright in this limit. The degree of calorimetry at fixed galaxy properties is sensitive to the assumed model for CR transport, and in particular to the time CRs spend interacting with neutral ISM, where they undergo strong streaming losses. We also find that in some regimes of parameter space hydrostatic equilibrium discs cannot exist, and in Paper II of this series we use this result to derive a critical surface in the plane of star formation surface density and gas surface density beyond which CRs may drive large-scale galactic winds.

The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) $V-$band light curves for sources brighter than V$\lesssim17$ mag across the whole sky. We produced V-band light curves for a total of ${\sim}61.5$ million sources and systematically searched these sources for variability. We identified ${\sim} 426,000$ variables, including ${\sim} 219,000$ new discoveries. Most (${\sim}74\%$) of our discoveries are in the Southern hemisphere. Here we use spectroscopic information from LAMOST, GALAH, RAVE, and APOGEE to study the physical and chemical properties of these variables. We find that metal-poor eclipsing binaries have orbital periods that are shorter than metal-rich systems at fixed temperature. We identified rotational variables on the main-sequence, red giant branch and the red clump. A substantial fraction (${\gtrsim}80\%$) of the rotating giants have large $v_{\rm rot}$ or large NUV excesses also indicative of fast rotation. The rotational variables have unusual abundances suggestive of analysis problems. Semi-regular variables tend to be lower metallicity ($\rm [Fe/H]{\sim}-0.5$) than most giant stars. We find that the APOGEE DR16 temperatures of oxygen-rich semi-regular variables are strongly correlated with the $W_{RP}-W_{JK}$ color index for $\rm T_{eff}\lesssim3800$ K. Using abundance measurements from APOGEE DR16, we find evidence for Mg, O and N enrichment in the semi-regular variables. We find that the Aluminum abundances of the semi-regular variables are strongly correlated with the pulsation period, where the variables with $\rm P\gtrsim 60$ days are significantly depleted in Al.

We utilize deep Chandra X-ray Observatory imaging and spectra of M82, the prototype of a starbursting galaxy with a multiphase wind, to map the hot plasma properties along the minor axis of the galaxy. We extract spectra from 11 regions up to 2.5 kpc from the starbursting midplane and model the data as a multi-temperature, optically thin thermal plasma with contributions from a non-thermal (power-law) component and from charge exchange (CX). We examine the gradients in best-fit parameters, including the intrinsic column density, plasma temperature, metal abundances, and number density of the hot gas as a function of distance from the M82 nucleus. We find that the temperatures and number densities of the warm-hot and hot plasma peak at the starbursting ridge and decreases along the minor axis. The temperature and density profiles are inconsistent with spherical adiabatic expansion of a super-heated wind and suggest mass loading and mixing of the hot phase with colder material. Non-thermal emission is detected in all of the regions considered, and CX comprises 7-22% of the total absorption-corrected, broad-band (0.5-7 keV) X-ray flux. We show that the abundances of O, Ne, Mg, and Fe are roughly constant across the regions considered, while Si and S peak within 500 pc of the central starburst. These findings support a direct connection between the M82 superwind and the warm-hot, metal-rich circumgalactic medium (CGM).

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

We present observations of ASASSN-19dj, a nearby tidal disruption event (TDE) discovered in the post-starburst galaxy KUG 0810+227 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of d $\simeq$ 98 Mpc. We observed ASASSN-19dj from $-$21 to 392 days relative to peak UV/optical emission using high-cadence, multi-wavelength spectroscopy and photometry. From the ASAS-SN $g$-band data, we determine that the TDE began to brighten on 2019 February 6.8 and for the first 25 days the rise was consistent with a flux $\propto$ $t^2$ power-law. ASASSN-19dj peaked in the UV/optical on 2019 March 6.5 (MJD = 58548.5) at a bolometric luminosity of $L = (6.2 \pm 0.2) \times 10^{44} \text{ erg s}^{-1}$. Initially remaining roughly constant in X-rays and slowly fading in the UV/optical, the X-ray flux increased by over an order of magnitude $\sim$225 days after peak, resulting from the expansion of the X-ray emitting surface. The late-time X-ray emission is well-fit by a blackbody with an effective radius of $\sim 1 \times 10^{12} \text{ cm}$ and a temperature of $\sim 6 \times 10^{5} \text{ K}$. Analysis of Catalina Real-Time Transient Survey images reveals a nuclear outburst roughly 14.5 years earlier with a smooth decline and a luminosity of $L_V$ $\geq$ $1.4 \times 10^{43}$ erg s$^{-1}$, although the nature of the flare is unknown. ASASSN-19dj occurred in the most extreme post-starburst galaxy yet to host a TDE, with Lick H$\delta_{A}$ = $7.67 \pm 0.17$ \AA.

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

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.

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

We present the fourth of the Cholla Galactic OutfLow Simulations suite (CGOLS). Using a physically-motivated prescription for clustered supernova feedback, we successfully drive a multiphase outflow from a disk galaxy. The high resolution ($< 5\,\mathrm{pc}$) across a relatively large domain ($20\,\mathrm{kpc}$) allows us to capture the hydrodynamic mixing and dynamical interactions between the hot and cool ($T \sim 10^4\,\mathrm{K}$) phases in the outflow, which in turn leads to direct evidence of a qualitatively new mechanism for cool gas acceleration in galactic winds. We show that mixing of momentum from the hot phase to the cool phase accelerates the cool gas to $800\,\mathrm{km}\,\mathrm{s}^{-1}$ on kpc scales, with properties inconsistent with the physical models of ram pressure acceleration or with bulk cooling from the hot phase. The mixing process also affects the hot phase, modifying its radial profiles of temperature, density, and velocity from the expectations of radial supersonic flow. This mechanism provides a physical explanation for the high velocity, blue shifted, low ionization absorption lines often observed in the spectra of starburst and high redshift galaxies.

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$\alpha$ 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$\alpha$ 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$\alpha$ luminosity remains approximately constant after its initial detection at phase +37 d, and that the H$\alpha$ velocity evolution does not trace that of the Fe~III$~\lambda4660$ emission. These suggest that the H$\alpha$ 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$\alpha$ 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$\alpha$ from CSM interaction to be discovered in an early-type galaxy.

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}$.

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 (A sub-type) EW binaries from the shorter period, late-type (W sub-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.

Although initially thought to be promising for production of the r-process nuclei, standard models of neutrino-heated winds from proto-neutron stars (PNSs) do not reach the requisite neutron-to-seed ratio for production of the lanthanides and actinides. However, the abundance distribution created by the r-, rp-, or $\nu p$-processes in PNS winds depends sensitively on the entropy and dynamical expansion timescale of the flow, which may be strongly affected by high magnetic fields. Here, we present results from magnetohydrodynamic simulations of non-rotating neutrino-heated PNS winds with strong dipole magnetic fields from $10^{14}-10^{16}$ G, and assess their role in altering the conditions for nucleosynthesis. The strong field forms a closed zone and helmet streamer configuration at the equator, with episodic dynamical mass ejections in toroidal plasmoids. We find dramatically enhanced entropy in these regions and conditions favorable for third-peak r-process nucleosynthesis if the wind is neutron-rich. If instead the wind is proton-rich, the conditions will affect the abundances from the $\nu p$-process. We quantify the distribution of ejected matter in entropy and dynamical expansion timescale, and the critical magnetic field strength required to affect the entropy. For $B\gtrsim10^{15}$ G, we find that $\gtrsim10^{-6}$ M$_{\odot}$ and up to $\sim10^{-5}$ M$_{\odot}$ of high entropy material is ejected per highly-magnetized neutron star birth in the wind phase, providing a mechanism for prompt heavy element enrichment of the universe. Former binary companions identified within (magnetar-hosting) supernova remnants, the remnants themselves, and runaway stars may exhibit overabundances. We provide a comparison with a semi-analytic model of plasmoid eruption and discuss implications and extensions.

Galactic winds from star-forming galaxies play at key role in the evolution of galaxies and the inter-galactic medium. They transport metals out of galaxies, chemically-enriching the inter-galactic medium and modifying the chemical evolution of galaxies. They affect the surrounding inter-stellar and circum-galactic media, thereby influencing the growth of galaxies through gas accretion and star-formation. In this contribution we first summarize the physical mechanisms by which the momentum and energy output from a population of massive stars and associated supernovae can drive galactic winds. We use the proto-typical example of M82 to illustrate the multiphase nature of galactic winds. We then describe how the basic properties of galactic winds are derived from the data, and summarize how the properties of galactic winds vary systematically with the properties of the galaxies that launch them. We conclude with a brief discussion of the broad implications of galactic winds.

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\sigma$ 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\sigma$ 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 $\Delta m_{15} (B)$.

We characterize an all-sky catalog of ${\sim} 8,400$ $\delta$ Scuti variables in ASAS-SN, which includes ${\sim} 3,300$ new discoveries. Using distances from 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 $\delta$ Scuti stars discovered by ASAS-SN. From this sample, we identified two new $\delta$ 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.

We present the discovery of PS18kh, a tidal disruption event (TDE) discovered at the center of SDSS J075654.53+341543.6 ($d\simeq322$ Mpc) by the Pan-STARRS Survey for Transients. Our dataset includes pre-discovery survey data from Pan-STARRS, the All-Sky Automated Survey for Supernovae (ASAS-SN), and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as well as high-cadence, multi-wavelength follow-up data from ground-based telescopes and Swift, spanning from 56 days before peak light until 75 days after. The optical/UV emission from PS18kh is well-fit as a blackbody with temperatures ranging from $T\simeq12000$ K to $T\simeq25000$ K and it peaked at a luminosity of $L\simeq8.8\times10^{43}$ ergs s$^{-1}$. PS18kh radiated $E=(3.45\pm0.22)\times10^{50}$ ergs over the period of observation, with $(1.42\pm0.20)\times10^{50}$ ergs being released during the rise to peak. Spectra of PS18kh show a changing, boxy/double-peaked H$\alpha$ emission feature, which becomes more prominent over time. We use models of non-axisymmetric accretion disks to describe the profile of the H$\alpha$ line and its evolution. We find that at early times the high accretion rate leads the disk to emit a wind which modifies the shape of the line profile and makes it bell-shaped. At late times, the wind becomes optically thin, allowing the non-axisymmetric perturbations to show up in the line profile. The line-emitting portion of the disk extends from $r_{\rm in}\sim60r_{\rm g}$ to an outer radius of $r_{\rm out}\sim1400r_{\rm g}$ and the perturbations can be represented either as an eccentricity in the outer rings of the disk or as a spiral arm in the inner disk.

We present the discovery and early evolution of ASASSN-19bt, a tidal disruption event (TDE) discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $d\simeq115$ Mpc and the first TDE to be detected by TESS. As the TDE is located in the TESS Continuous Viewing Zone, our dataset includes 30-minute cadence observations starting on 2018 July 25, and we precisely measure that the TDE begins to brighten $\sim8.3$ days before its discovery. Our dataset also includes 18 epochs of Swift UVOT and XRT observations, 2 epochs of XMM-Newton observations, 13 spectroscopic observations, and ground data from the Las Cumbres Observatory telescope network, spanning from 32 days before peak through 37 days after peak. ASASSN-19bt thus has the most detailed pre-peak dataset for any TDE. The TESS light curve indicates that the transient began to brighten on 2019 January 21.6 and that for the first 15 days its rise was consistent with a flux $\propto t^2$ power-law model. The optical/UV emission is well-fit by a blackbody SED, and ASASSN-19bt exhibits an early spike in its luminosity and temperature roughly 32 rest-frame days before peak and spanning up to 14 days that has not been seen in other TDEs, possibly because UV observations were not triggered early enough to detect it. It peaked on 2019 March 04.9 at a luminosity of $L\simeq1.3\times10^{44}$ ergs s$^{-1}$ and radiated $E\simeq3.2\times10^{50}$ ergs during the 41-day rise to peak. X-ray observations after peak indicate a softening of the hard X-ray emission prior to peak, reminiscent of the hard/soft states in X-ray binaries.

In this catalog we compile information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) as well as all other bright ($m_{peak}\leq17$), spectroscopically confirmed supernovae found in 2017, totaling 308 supernovae. We also present UV through near-IR magnitudes gathered from public databases of all host galaxies for the supernovae in the sample. We perform statistical analyses of our full bright supernova sample, which now contains 949 supernovae discovered since 2014 May 1, including supernovae from our previous catalogs. This is the fourth of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team, and this work presents updated data and measurements, including light curves, redshifts, classifications, and host galaxy identifications, that supersede information contained in any previous publications.

Cosmic rays (CRs) and magnetic fields may be dynamically important in driving large-scale galactic outflows from rapidly star-forming galaxies. We construct two-dimensional axisymmetric models of the local starburst and super-wind galaxy M82 using the CR propagation code GALPROP. Using prescribed gas density and magnetic field distributions, wind profiles, CR injection rates, and stellar radiation fields, we simultaneously fit both the integrated gamma-ray emission and the spatially-resolved multi-frequency radio emission extended along M82's minor axis. We explore the resulting constraints on the gas density, magnetic field strength, CR energy density, and the assumed CR advection profile. In accord with earlier one-zone studies, we generically find low central CR pressures, strong secondary electron/positron production, and an important role for relativistic bremsstrahlung losses in shaping the synchrotron spectrum. We find that the relatively low central CR density produces CR pressure gradients that are weak compared to gravity, strongly limiting the role of CRs in driving M82's fast and mass-loaded galactic outflow. Our models require strong magnetic fields and advection speeds of order ~1000 km/s on kpc scales along the minor axis in order to reproduce the extended radio emission. Degeneracies between the controlling physical parameters of the model and caveats to these findings are discussed.

The isolated binary evolution model for merging neutron stars (NSs) involves processes such as mass transfer, common-envelope evolution, and natal kicks, all of which are poorly understood. Also, the predicted NS-NS merger rates are typically lower than the rates inferred from the LIGO GW170817 event. Here, we investigate merger rates of NS and black hole (BH)-NS binaries in hierarchical triple-star systems. In such systems, the tertiary can induce Lidov-Kozai (LK) oscillations in the inner binary, accelerating its coalescence, and potentially enhancing compact object merger rates. However, since compact objects originate from massive stars, the prior evolution should also be taken into account. Natal kicks, in particular, could significantly reduce the rates by unbinding the tertiary before it can affect the inner binary through LK evolution. We carry out simulations of massive triples taking into account stellar evolution starting from the main sequence, secular and tidal evolution, and the effects of supernovae. For large NS birth kicks (\sigma_k = 265 km/s), we find that the triple NS-NS merger rate (several hundred Gpc^{-3} yr^{-1}) is lower by a factor of ~2-3 than the binary rate, but for no kicks (\sigma_k = 0 km/s), the triple rate (several thousand Gpc^{-3} yr^{-1}) is comparable to the binary rate. Our results indicate that a significant fraction of NS-NS mergers could originate from triples if a substantial portion of the NS population is born with low kick velocities, as indicated by other work. However, uncertainties and open questions remain because of our simplifying assumption of dynamical decoupling after inner binary interaction has been triggered.

After the identification of the gamma-ray blazar TXS 0506+056 as the first compelling IceCube neutrino source candidate, we perform a systematic analysis of all high-energy neutrino events satisfying the IceCube realtime trigger criteria. We find one additional known gamma-ray source, the blazar GB6 J1040+0617, in spatial coincidence with a neutrino in this sample. The chance probability of this coincidence is 30% after trial correction. For the first time, we present a systematic study of the gamma-ray flux, spectral and optical variability, and multi-wavelength behavior of GB6 J1040+0617 and compare it to TXS 0506+056. We find that TXS 0506+056 shows strong flux variability in the Fermi-LAT gamma-ray band, being in an active state around the arrival of IceCube-170922A, but in a low state during the archival IceCube neutrino flare in 2014/15. In both cases the spectral shape is statistically compatible ($\leq 2\sigma$) with the average spectrum showing no indication of a significant relative increase of a high-energy component. While the association of GB6 J1040+0617 with the neutrino is consistent with background expectations, the source appears to be a plausible neutrino source candidate based on its energetics and multi-wavelength features, namely a bright optical flare and modestly increased gamma-ray activity. Finding one or two neutrinos originating from gamma-ray blazars in the given sample of high-energy neutrinos is consistent with previously derived limits of neutrino emission from gamma-ray blazars, indicating the sources of the majority of cosmic high-energy neutrinos remain unknown.

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$ $\delta$-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.

We report comprehensive multi-wavelength observations of a peculiar Type Ia-like supernova ("SN Ia-pec") ASASSN-15pz. ASASSN-15pz is a spectroscopic "twin" of SN 2009dc, a so-called "Super-Chandrasekhar-mass" SN, throughout its evolution, but it has a peak luminosity M_B,peak = -19.69 +/- 0.12 mag that is \approx 0.6 mag dimmer and comparable to the SN 1991T sub-class of SNe Ia at the luminous end of the normal width-luminosity relation. The synthesized Ni56 mass of M_Ni56 = 1.13 +/- 0.14 M_sun is also substantially less than that found for several 2009dc-like SNe. Previous well-studied 2009dc-like SNe have generally suffered from large and uncertain amounts of host-galaxy extinction, which is negligible for ASASSN-15pz. Based on the color of ASASSN-15pz, we estimate a host extinction for SN 2009dc of E(B-V)_host=0.12 mag and confirm its high luminosity (M_B, peak[2009dc] \approx -20.3 mag). The 2009dc-like SN population, which represents ~1% of SNe Ia, exhibits a range of peak luminosities, and do not fit onto the tight width-luminosity relation. Their optical light curves also show significant diversity of late-time (>~ 50 days) decline rates. The nebular-phase spectra provide powerful diagnostics to identify the 2009dc-like events as a distinct class of SNe Ia. We suggest referring to these sources using the phenomenology-based "2009dc-like SN Ia-pec" instead of "Super-Chandrasekhar SN Ia," which is based on an uncertain theoretical interpretation.

Type Ia supernovae (SNe Ia) could arise from mergers of carbon-oxygen white dwarfs (WDs) triggered by Lidov-Kozai (LK) oscillations in hierarchical triple-star systems. However, predicted merger rates are several orders of magnitude lower than the observed SNe Ia rate. The low predicted rates can be attributed in part to the fact that many potential WD-merger progenitor systems, with high mutual orbital inclination, merge or interact before the WD stage. Recently, evidence was found for the existence of natal kicks imparted on WDs with a typical magnitude of 0.75 km/s. In triples, kicks change the mutual inclination and in general increase the outer orbit eccentricity, bringing the triple into an active LK regime at late stages and avoiding the issue of pre-WD merger or interaction. Stars passing by the triple can result in similar effects. However, both processes can also disrupt the triple. In this paper, we quantitatively investigate the impact of WD kicks and flybys on the rate of WD mergers using detailed simulations. We find that WD kicks and flybys combine to increase the predicted WD merger rates by a factor of ~2.5, resulting in a time-integrated rate of ~1.1e-4 Msun^-1. Despite the significant boost, the predicted rates are still more than one order of magnitude below the observed rate of ~1e-3 Msun^-1. However, many systematic uncertainties still remain in our calculations, in particular the potential contributions from tighter triples, dynamically unstable systems, unbound systems due to WD kicks, and quadruple systems.

We explore the synergy between photometric and spectroscopic surveys by searching for periodic variable stars among the targets observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) using photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN). We identified 1924 periodic variables among more than 258000 APOGEE targets; 465 are new discoveries. We homogeneously classified 430 eclipsing and ellipsoidal binaries, 139 classical pulsators (Cepheids, RR Lyrae and delta Scuti), 719 long period variables (pulsating red giants) and 636 rotational variables. The search was performed using both visual inspection and machine learning techniques. The light curves were also modeled with the damped random walk stochastic process. We find that the median [Fe/H] of variable objects is lower by 0.3 dex than that of the overall APOGEE sample. Eclipsing binaries and ellipsoidal variables are shifted to a lower median [Fe/H] by 0.2 dex. Eclipsing binaries and rotational variables exhibit significantly broader spectral lines than the rest of the sample. We make ASAS-SN light curves for all the APOGEE stars publicly available and provide parameters for the variable objects.

We report the discovery of the first likely black hole in a non-interacting binary system with a field red giant. By combining radial velocity measurements from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) with photometric variability data from the All-Sky Automated Survey for Supernovae (ASAS-SN), we identified the bright rapidly-rotating giant 2MASS J05215658+4359220 as a binary system with a massive unseen companion. Subsequent radial velocity measurements reveal a system with an orbital period of 83 days and near-zero eccentricity. The photometric variability period of the giant is consistent with the orbital period, indicative of star spots and tidal synchronization. Constraints on the giant's mass and radius from its luminosity, surface gravity, and temperature imply an unseen companion with mass of $3.3^{+2.8}_{-0.7}$ M$_\odot$, indicating a low-mass black hole or an exceedingly massive neutron star. Measurement of the astrometric binary motion by {\it Gaia} will further characterize the system. This discovery demonstrates the potential of massive spectroscopic surveys like APOGEE and all-sky, high-cadence photometric surveys like ASAS-SN to revolutionize our understanding of the compact object mass function, and to test theories of binary star evolution and the supernova mechanism.

We present a detailed study of the 2017 eruption of the classical nova ASASSN-17pf (LMCN 2017-11a), which is located in the Large Magellanic Cloud, including data from AAVSO, ASAS-SN, SALT, SMARTS, SOAR, and the Neil Gehrels \textit{Swift} Observatory. The optical light-curve is characterized by multiple maxima (flares) on top of a slowly evolving light-curve (with a decline time, $t_2>$ 100 d). The maxima correlate with the appearance of new absorption line systems in the optical spectra characterized by increasing radial velocities. We suggest that this is evidence of multiple episodes of mass-ejection with increasing expansion velocities. The line profiles in the optical spectra indicate very low expansion velocities (FWHM $\sim$ 190 km s$^{-1}$), making this nova one of the slowest expanding ever observed, consistent with the slowly evolving light-curve. The evolution of the colors and spectral energy distribution show evidence of decreasing temperatures and increasing effective radii for the pseudo-photosphere during each maximum. The optical and infrared light-curves are consistent with dust formation 125 days post-discovery. We speculate that novae showing several optical maxima have multiple mass-ejection episodes leading to shocks that may drive $\gamma$-ray emission and dust formation.

The observed binary black hole (BBH) mergers indicate a large Galactic progenitor population continuously evolving from large orbital separations and low gravitational wave (GW) frequencies to the final merger phase. We investigate the equilibrium distribution of binary black holes in the Galaxy. Given the observed BBH merger rate, we contrast the expected number of systems radiating in the low-frequency $0.1-10\,$mHz GW band under two assumptions: (1) that all merging systems originate from near-circular orbits, as may be indicative of isolated binary evolution, and (2) that all merging systems originate at very high eccentricity, as predicted by models of dynamically-formed BBHs and triple and quadruple systems undergoing Lidov-Kozai eccentricity oscillations. We show that the equilibrium number of systems expected at every frequency is higher in the eccentric case (2) than in the circular case (1) by a factor of $\simeq 2-15$. This follows from the fact that eccentric systems spend more time than circular systems radiating in the low-frequency GW bands. The GW emission comes in pulses at periastron separated by the orbital period, which may be days to years. For a LISA-like sensitivity curve, we show that if eccentric systems contribute significantly to the observed merger rate, then $\simeq 10$ eccentric systems should be seen in the Galaxy.

We report distinctly double-peakedH-alpha and H-beta emission lines in the late-time, nebular-phase spectra (>~200 d) of the otherwise normal at early phases (<~ 100 d) Type IIP supernova ASASSN-16at (SN 2016X). Such distinctly double-peaked nebular Balmer lines have never been observed for a Type II SN. The nebular-phase Balmer emission is driven by the radioactive Co56 decay, so the observed line-profile bifurcation suggests a strong bipolarity in the Ni56 distribution or in the line-forming region of the inner ejecta. The strongly bifurcated blue- and red-shifted peaks are separated by ~3x10^3 km/s and are roughly symmetrically positioned with respect to the host-galaxy rest frame, implying that the inner ejecta are composed of two almost detached blobs. The red peak progressively weakens relative to the blue peak, and disappears in the 740 d spectrum. One possible reason for the line-ratio evolution is increasing differential extinction from continuous formation of dust within the envelope, which is also supported by the near-infrared flux excess that develops after ~100 d.

We analyze the 476 SN Ia host galaxies from the All-Sky Automated Survey for Supernova (ASAS-SN) Bright Supernova Catalogs to determine the observed relative Type Ia supernova (SN) rates as a function of luminosity and host galaxy properties. We find that the luminosity distribution of the SNe Ia in our sample is reasonably well described by a Schechter function with a faint-end slope $\alpha \approx 1.5$ and a knee $M_{\star} \approx -18.0$. Our specific SN Ia rates are consistent with previous results but extend to far lower host galaxy masses. We find an overall rate that scales as $(M_{\star}/10^{10} M_{\odot})^{\alpha}$ with $\alpha \approx -0.5$. This shows that the specific SN Ia rate continues rising towards lower masses even in galaxies as small as $\log(M_{\star} / M_{\odot}) \lesssim 7.0$, where it is enhanced by a factor of $\sim10-20$ relative to host galaxies with stellar masses $\sim10^{10}M_{\odot}$. We find no strong dependence of the specific SN Ia rate on the star formation activity of the host galaxies, but additional observations are required to improve the constraints on the star formation rates.

Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. \citet{Pejcha2012} found that an `antesonic condition' exists for time-steady spherically symmetric models, potentially giving a physical explanation for the critical curve observed in simulations. In this paper, we extend that analysis to time-dependent, spherically symmetric polytropic models. We verify the critical antesonic condition in our simulations, showing that models exceeding it drive transonic winds whereas models below it exhibit steady accretion. In addition, we find that (1) high spatial resolution is needed for accurate determination of the antesonic ratio and shock radius at the critical curve, and that low resolution simulations systematically underpredict these quantities, making explosion more difficult at lower resolution; (2) there is an important physical connection between the critical mass accretion rate at explosion and the mass loss rate of the post-explosion wind: the two are directly proportional at criticality, implying that, at criticality, the wind kinetic power is tied directly to the accretion power; (3) the value of the post-shock adiabatic index $\Gamma$ has a large effect on the length and time scales of the post-bounce evolution of the explosion larger values of $\Gamma$ result in a longer transition from the accretion to wind phases.

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.

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 $\rm P_{orb}=32.83\pm0.05$ 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 $\omega=53.5$ deg, $i=45.7$ deg and $e=0.58$. We also carried out a frequency analysis to identify the pulsation frequencies corresponding to the tidally excited oscillations.

We report the discovery of 19 new R Coronae Borealis (RCB) star and DY Per candidates with light curves from the All-Sky Automated Survey for Supernovae (ASAS-SN). We examined both an existing set of 1602 near/mid-IR selected candidates and an additional 2615 candidates selected to have near/mid-IR SEDs consistent with those of known R Coronae Borealis stars. We visually inspected the light curves for the characteristic variability of these systems.

We present optical photometry and spectroscopy of the Type II supernova ASASSN-14jb, together with VLT MUSE IFU observations of its host galaxy and a nebular-phase spectrum. This SN, in the nearby galaxy ESO 467-G051 ($z=0.006$), was discovered and followed-up by the All Sky Automated Survey for SuperNovae (ASAS-SN). We present LCOGTN $BVgri$ and $Swift$ $w2m1w1ubv$ optical and near-UV/optical light curves and several optical spectra in the early photospheric phase. ASASSN-14jb exploded $\sim 2$ kpc above the star-forming disk of ESO 467-G051, an edge-on disk galaxy. The large projected distance from the disk and non-detection of any H II region in a 1.4 kpc radius in projection are in conflict with the standard environment of core-collapse supernova progenitors and suggests the possible scenario that the progenitor received a kick in a binary interaction. Our analysis yields $D=25\pm 1$ Mpc, $M(^{56}Ni) = 0.0210 \pm 0.0025$ M$_\odot$, $E_{exp}\approx 0.25 \times 10^{51}$ ergs, $M_{ej}\approx 6$ M$_\odot$, and progenitor $R_* = 580 \pm 28$ R$_\odot$, which seems to be consistent with the sub-Solar metallicity of $0.3 \pm 0.1$ Z$_\odot$ derived from the SN spectrum. The nebular spectrum constrains $M_{prog} \simeq 10-12 $ M$_{\odot}$. We discuss the low oxygen abundance of the host galaxy derived from the MUSE data, $12+\log{(O/H)} = 8.27^{+0.16}_{-0.20}$ (O3N2 strong line method), and compare it with the supernova spectra, which is also consistent with a sub-Solar metallicity progenitor. We finally discuss the possible scenarios for the unusual environment for ASASSN-14jb and conclude that either the in-situ star formation or a runaway would imply a low mass progenitor, agreeing with our estimate from the supernova nebular spectrum. We show that the detailed study of the environment can agree with the stronger constrains of the transient observations. (abridged)

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

We employ VLTI GRAVITY to resolve, for the first time, the two images generated by a gravitational microlens. The measurements of the image separation \theta_{-,+}=3.78 +/- 0.05 mas, and hence the Einstein radius \theta_E =1.87 +/- 0.03 mas, are precise. This demonstrates the robustness of the method, provided that the source is bright enough for GRAVITY (K <~ 10.5) and the image separation is of order or larger than the fringe spacing. When \theta_E is combined with a measurement of the "microlens parallax" \pi_E, the two will together yield the lens mass and lens-source relative parallax and proper motion. Because the source parallax and proper motion are well measured by Gaia, this means that the lens characteristics will be fully determined, whether or not it proves to be luminous. This method can be a powerful probe of dark, isolated objects, which are otherwise quite difficult to identify, much less characterize. Our measurement contradicts Einstein's (1936) prediction that "the luminous circle [i.e., microlensed image] cannot be distinguished" from a star.

NGC 253 hosts the nearest nuclear starburst. Previous observations show a region rich in molecular gas, with dense clouds associated with recent star formation. We used ALMA to image the 350 GHz dust continuum and molecular line emission from this region at 2 pc resolution. Our observations reveal ~14 bright, compact (~2-3 pc FWHM) knots of dust emission. Most of these sources are likely to be forming super star clusters (SSCs) based on their inferred dynamical and gas masses, association with 36 GHz radio continuum emission, and coincidence with line emission tracing dense, excited gas. One source coincides with a known SSC, but the rest remain invisible in Hubble near-infrared (IR) imaging. Our observations imply that gas still constitutes a large fraction of the overall mass in these sources. Their high brightness temperature at 350 GHz also implies a large optical depth near the peak of the IR spectral energy distribution. As a result, these sources may have large IR photospheres and the IR radiation force likely exceeds L/c. Still, their moderate observed velocity dispersions suggest that feedback from radiation, winds, and supernovae are not yet disrupting most sources. This mode of star formation appears to produce a large fraction of stars in the burst. We argue for a scenario in which this phase lasts ~1 Myr, after which the clusters shed their natal cocoons but continue to produce ionizing photons. The strong feedback that drives the observed cold gas and X-ray outflows likely occurs after the clusters emerge from this early phase.

We present the discovery of ASASSN-18ey (MAXI J1820+070), a new black hole low-mass X-ray binary discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN). A week after ASAS-SN discovered ASASSN-18ey as an optical transient, it was detected as an X-ray transient by MAXI/GCS. Here, we analyze ASAS-SN and Asteroid Terrestrial-impact Last Alert System (ATLAS) pre-outburst optical light curves, finding evidence of intrinsic variability for several years prior to the outburst. While there was no long-term rise leading to outburst, as has been seen in several other systems, the start of the outburst in the optical preceded that in the X-rays by $7.20\pm0.97~\rm days$. We analyze the spectroscopic evolution of ASASSN-18ey from pre-maximum to $> 100~\rm days$ post-maximum. The spectra of ASASSN-18ey exhibit broad, asymmetric, double-peaked H$\alpha$ emission. The Bowen blend ($\lambda \approx 4650$\AA) in the post-maximum spectra shows highly variable double-peaked profiles, likely arising from irradiation of the companion by the accretion disk, typical of low-mass X-ray binaries. The optical and X-ray luminosities of ASASSN-18ey are consistent with black hole low-mass X-ray binaries, both in outburst and quiescence.

The large columns of dusty gas enshrouding and fuelling star-formation in young, massive stellar clusters may render such systems optically thick to radiation well into the infrared. This raises the prospect that both "direct" radiation pressure produced by absorption of photons leaving stellar surfaces and "indirect" radiation pressure from photons absorbed and then re-emitted by dust grains may be important sources of feedback in such systems. Here we evaluate this possibility by deriving the conditions under which a spheroidal, self-gravitating, mixed gas-star cloud can avoid catastrophic disruption by the combined effects of direct and indirect radiation pressure. We show that radiation pressure sets a maximum star cluster formation efficiency of $\epsilon_{\rm max} \sim 0.9$ at a (very large) gas surface density of $\sim 10^5 M_\odot$ pc$^{-2} (Z_\odot/Z) \simeq 20$ g cm$^{-2} (Z_\odot/Z)$, but that gas clouds above this limit undergo significant radiation-driven expansion during star formation, leading to a maximum stellar surface density very near this value for all star clusters. Data on the central surface mass density of compact stellar systems, while sparse and partly confused by dynamical effects, are broadly consistent with the existence of a metallicity-dependent upper-limit comparable to this value. Our results imply that this limit may preclude the formation of the progenitors of intermediate-mass black holes for systems with $Z \gtrsim 0.2 Z_\odot$.

We report the identification of the bright (V${\sim}13.3$ mag) star FY Sct as a long period detached eclipsing binary using a combined ASAS-SN and ASAS light curve spanning 2000-2018. The orbital period is P${\sim}2.57$ years and the primary eclipse lasts ${\sim}73$ days. The eclipse profile is suggestive of a disk eclipsing binary rather than a stellar component. We also detect ${\sim}0.4$ mag pulsations with a period of P${\sim78}$ d. The next eclipse begins on September, 28, 2018. Further photometric and spectroscopic observations are encouraged, particularly when the system is in eclipse.

The All-Sky Automated Survey for Supernovae (ASAS-SN) is the first optical survey to routinely monitor the whole sky with a cadence of $\sim2-3$ days down to V$\lesssim17$ mag. ASAS-SN has monitored the whole sky since 2014, collecting $\sim100-500$ epochs of observations per field. The V-band light curves for candidate variables identified during the search for supernovae are classified using a random forest classifier and visually verified. We present a catalog of 66,179 bright, new variable stars discovered during our search for supernovae, including 27,479 periodic variables and 38,700 irregular variables. V-band light curves for the ASAS-SN variables are available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). The database will begin to include the light curves of known variable stars in the near future along with the results for a systematic, all-sky variability survey.

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}44,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.

Successful models of the low redshift circumgalactic medium (CGM) must account for (1) a large amount of gas, (2) relatively slow gas velocities, (3) a high degree of metal enrichment, (4) the similar absorption properties around both star-forming and passive galaxies, and (5) the observationally inferred temperature and densities of the CGM gas. We show that galactic wind-driven bubbles can account for these observed properties. We develop a model describing the motion of bubbles driven by a hot, fast galactic wind characteristic of supernova energy injection. The bubble size grows slowly to hundreds of kiloparsecs over $5-10$ Gyr. For high star formation rates or high wind mass loading $\dot M_w/\dot M_\star$, the free-flowing hot wind, the shocked hot wind in the interior of the bubble, and the swept-up halo gas within the bubble shell can all radiatively cool and contribute to low-ionization state metal line absorption. We verify that if the free-flowing wind cools, the shocked wind does as well. We find effective mass loading factors of $(M_w+M_\mathrm{swept})/M_\star\sim5-12$ as the bubbles sweep into the CGM. We predict cool gas masses, velocities, column densities, metal content, and absorption line velocities and linewidths of the bubble for a range of parameter choices. This picture can reproduce many of the COS-Halos and Keeney et al. (2017) observations of low-ionization state metal absorption lines around both star-forming and passive galaxies.

We present the ATLAS discovery and initial analysis of the first 18 days of the unusual transient event, ATLAS18qqn/AT2018cow. It is characterized by a high peak luminosity ($\sim$1.7 $\times$ 10$^{44}$ erg s$^{-1}$), rapidly evolving light curves ($>$5 mag rise to peak in $\sim$3.5 days), and hot blackbody spectra, peaking at $\sim$27000 K that are relatively featureless and unchanging over the first two weeks. The bolometric light curve cannot be powered by radioactive decay under realistic assumptions. The detection of high-energy emission may suggest a central engine as the powering source. Using a magnetar model, we estimated an ejected mass of $0.1-0.4$ \msol, which lies between that of low-energy core-collapse events and the kilonova, AT2017gfo. The spectra cooled rapidly from 27000 to 15000 K in just over 2 weeks but remained smooth and featureless. Broad and shallow emission lines appear after about 20 days, and we tentatively identify them as He I although they would be redshifted from their rest wavelengths. We rule out that there are any features in the spectra due to intermediate mass elements up to and including the Fe-group. The presence of r-process elements cannot be ruled out. If these lines are due to He, then we suggest a low-mass star with residual He as a potential progenitor. Alternatively, models of magnetars formed in neutron-star mergers give plausible matches to the data.

Galactic outflows commonly contain multiphase gas, and its physical origin requires explanation. Using the CGOLS (Cholla Galactic OutfLow Simulations) suite of high-resolution isolated galaxy models, we demonstrate the viability of rapid radiative cooling as a source of fast-moving ($v \sim 1000$ km/s), cool ($10^4$ K) gas observed in absorption line studies of outflows around some star-forming galaxies. By varying the mass-loading and geometry of the simulated winds, we identify a region of parameter space that leads to cool gas in outflows. In particular, when using an analytically-motivated central feedback model, we find that cooling flows can be produced with reasonable mass-loading rates ($\dot{M}_{wind} / \dot{M}_{SFR} \sim 0.5$), provided the star formation rate surface density is high. When a more realistic clustered feedback model is applied, destruction of high density clouds near the disk and interactions between different outflow regions indicate that lower mass-loading rates of the hot gas within the feedback region may still produce multiphase outflows. These results suggest an origin for fast-moving cool gas in outflows that does not rely on directly accelerating cool gas from the interstellar medium. These cooling flows may additionally provide an explanation for the multiphase gas ubiquitously observed in the halos of star-forming galaxies at low redshift.

Galaxy formation simulations demonstrate that cosmic-ray (CR) feedback may be important in the launching of galactic-scale winds. CR protons dominate the bulk of the CR population, yet most observational constraints of CR feedback come from synchrotron emission of CR electrons. In this paper, we present an analysis of 105 months of Fermi Gamma-ray Space Telescope observations of the Small Magellanic Cloud (SMC), with the aim of exploring CR feedback and transport in an external galaxy. We produce maps of the 2-300 GeV emission and detect statistically significant, extended emission along the Bar and the Wing, where active star formation is occurring. Gamma-ray emission is not detected above 13 GeV, and we set stringent upper-limits on the flux above this energy. We find the best fit to the gamma-ray spectrum is a single-component model with a power-law of index $\Gamma=-2.11\pm0.06\pm0.06$ and an exponential cutoff energy of $E_{\rm c} =13.1\pm5.1\pm1.6$ GeV. We assess the relative contribution of pulsars and CRs to the emission, and we find that pulsars may produce up to 14$^{+4}_{-2}$% of the flux above 100 MeV. Thus, we attribute most of the gamma-ray emission (based on its spectrum and morphology) to CR interactions with the ISM. We show that the gamma-ray emissivity of the SMC is five times smaller than that of the Milky Way and that the SMC is far below the calorimetric limit, where all CR protons experience pion losses. We interpret these findings as evidence that CRs are escaping the SMC via advection and diffusion.

We report the identification of a new, long period detached eclipsing binary by the All-Sky Automated Survey for Supernovae (ASAS-SN). Combining V-band photometry from ASAS-SN, ASAS and KELT, we constructed a light curve spanning 2006-2018. ASASSN-V J192543.72+402619.0 is an eccentric detached eclipsing binary with an orbital period of $P\sim7.3$ yr, a minimum orbital eccentricity of $e_{min}\sim0.2$ and a primary eclipse lasting $\sim160$ d. We also detect a secondary eclipse that is $\sim 0.6$ mag shallower. Multi-band photometry suggests that the more luminous component of the binary is an M-giant. Further photometric and spectroscopic observations are encouraged.

The importance of radiation pressure feedback in galaxy formation has been extensively debated over the last decade. The regime of greatest uncertainty is in the most actively star-forming galaxies, where large dust columns can potentially produce a dust-reprocessed infrared radiation field with enough pressure to drive turbulence or eject material. Here we derive the conditions under which a self-gravitating, mixed gas-star disc can remain hydrostatic despite trapped radiation pressure. Consistently taking into account the self-gravity of the medium, the star- and dust-to-gas ratios, and the effects of turbulent motions not driven by radiation, we show that galaxies can achieve a maximum Eddington-limited star formation rate per unit area $\dot{\Sigma}_{\rm *,crit} \sim 10^3 M_{\odot}$ pc$^{-2}$ Myr$^{-1}$, corresponding to a critical flux of $F_{\rm *,crit} \sim 10^{13} L_{\odot}$ kpc$^{-2}$ similar to previous estimates; higher fluxes eject mass in bulk, halting further star formation. Conversely, we show that in galaxies below this limit, our one-dimensional models imply simple vertical hydrostatic equilibrium and that radiation pressure is ineffective at driving turbulence or ejecting matter. Because the vast majority of star-forming galaxies lie below the maximum limit for typical dust-to-gas ratios, we conclude that infrared radiation pressure is likely unimportant for all but the most extreme systems on galaxy-wide scales. Thus, while radiation pressure does not explain the Kennicutt-Schmidt relation, it does impose an upper truncation on it. Our predicted truncation is in good agreement with the highest observed gas and star formation rate surface densities found both locally and at high redshift.

We report and characterize a white-light superflare on a previously undiscovered M dwarf detected by the ASAS-SN survey. Employing various color-magnitude and color-spectral type relationships, we estimate several stellar parameters, including the quiescent V-band magnitude, from which we derive a flare amplitude of $\Delta V \sim 10$. We determine an r-band absolute magnitude of $M_{r} = 11.4$, consistent with a mid-M dwarf, and an approximate distance to the source of $2.2$ kpc. Using classical-flare models, we infer a flare energy of $E_{V} \simeq (4.1\pm 2.2)\times 10^{36}$ ergs, making this one of the strongest flares documented on an M dwarf.

[Abridged]We investigate the long-term secular dynamics and Lidov-Kozai(LK) eccentricity oscillations of quadruple systems composed of two binaries at quadrupole and octupole order in the perturbing Hamiltonian. We show that the fraction of systems reaching high eccentricities is enhanced relative to triple systems, over a broader range of parameter space. We show that this fraction grows with time, unlike triple systems evolved at quadrupole order. This is fundamentally because with their additional degrees of freedom, quadruple systems do not have a maximal set of commuting constants of the motion, even in secular theory at quadrupole order. We discuss these results in the context of star-star and white dwarf-white dwarf(WD) binaries, with emphasis on WD-WD mergers and collisions relevant to the Type Ia supernova problem. For star-star systems, we find that more than 30% of systems reach high eccentricity within a Hubble time, potentially forming triple systems via stellar mergers or close binaries. For WD-WD systems, taking into account general relativistic and tidal precession and dissipation, we show that the merger rate is enhanced in quadruple systems relative to triple systems by a factor of 3.5-10, and that the long-term evolution of quadruple systems leads to a delay-time distribution ~1/t for mergers and collisions. In GW-driven mergers of compact objects, we classify the mergers by their evolutionary patterns in phase space and identify a regime in ~8% of orbital shrinking mergers, where eccentricity oscillations occur on the general relativistic precession timescale, rather than the much longer LK timescale. Finally, we generalize previous treatments of oscillations in the inner binary eccentricity to eccentric mutual orbits. We assess the merger rate in quadruple and triple systems and the implications for their viability as progenitors of stellar mergers and SNe Ia.

The neutron star (NS) merger GW170817 was followed over several days by optical-wavelength ("blue") kilonova (KN) emission likely powered by the radioactive decay of light r-process nuclei synthesized by ejecta with a low neutron abundance (electron fraction Ye ~ 0.25-0.35). While the composition and high velocities of the blue KN ejecta are consistent with shock-heated dynamical material, the large quantity is in tension with the results of numerical simulations. We propose an alternative ejecta source: the neutrino-heated, magnetically-accelerated wind from the strongly-magnetized hypermassive NS (HMNS) remnant. A rapidly-spinning HMNS with an ordered surface magnetic field of strength B ~ 1-3e14 G and lifetime t_rem ~ 0.1-1 s can simultaneously explain the velocity, total mass, and electron fraction of the blue KN ejecta. The inferred HMNS lifetime is close to its Alfven crossing time, suggesting global magnetic torques could be responsible for bringing the HMNS into solid body rotation and instigating its gravitational collapse. Different origins for the KN ejecta may be distinguished by their predictions for the emission in the first hours after the merger, when the luminosity is enhanced by heating from internal shocks; the latter are likely generic to any temporally-extended ejecta source (e.g. magnetar or accretion disk wind) and are not unique to the emergence of a relativistic jet. The same shocks could mix and homogenizes the composition to a low but non-zero lanthanide mass fraction, X_La ~ 1e-3, as advocated by some authors, but only if the mixing occurs after neutrons are consumed in the r-process on a timescale >~ 1 s.

We present a catalog of 4880 newly identified RR Lyrae variable stars (4433 RRab, 446 RRc and 1 RRd) found during the search for supernovae by the All-Sky Automated Survey for Supernovae (ASAS-SN). The light curves, classified using the "Upsilon" random forest classifier and visually verified, are available through the ASAS-SN variable star database (https://asas-sn.osu.edu/variables). The database also contains light curves for an expanding number of additional new variables (~50,000 at present) and will begin to include the light curves of known variable stars in the near future.

We use the multi-epoch radial velocities acquired by the APOGEE survey to perform a large scale statistical study of stellar multiplicity for field stars in the Milky Way, spanning the evolutionary phases between the main sequence and the red clump. We show that the distribution of maximum radial velocity shifts (\drvm) for APOGEE targets is a strong function of \logg, with main sequence stars showing \drvm\ as high as $\sim$300 \kms, and steadily dropping down to $\sim$30 \kms\ for \logg$\sim$0, as stars climb up the Red Giant Branch (RGB). Red clump stars show a distribution of \drvm\ values comparable to that of stars at the tip of the RGB, implying they have similar multiplicity characteristics. The observed attrition of high \drvm\ systems in the RGB is consistent with a lognormal period distribution in the main sequence and a multiplicity fraction of 0.35, which is truncated at an increasing period as stars become physically larger and undergo mass transfer after Roche Lobe Overflow during H shell burning. The \drvm\ distributions also show that the multiplicity characteristics of field stars are metallicity dependent, with metal-poor ([Fe/H]$\lesssim-0.5$) stars having a multiplicity fraction a factor 2-3 higher than metal-rich ([Fe/H]$\gtrsim0.0$) stars. This has profound implications for the formation rates of interacting binaries observed by astronomical transient surveys and gravitational wave detectors, as well as the habitability of circumbinary planets.

The radiation force on dust grains may be dynamically important in driving turbulence and outflows in rapidly star-forming galaxies. Recent studies focus on the highly optically-thick limit relevant to the densest ultra-luminous galaxies and super star clusters, where reprocessed infrared photons provide the dominant source of electromagnetic momentum. However, even among starburst galaxies, the great majority instead lie in the so-called "single-scattering" limit, where the system is optically-thick to the incident starlight, but optically-thin to the re-radiated infrared. In this paper we present a stability analysis and multidimensional radiation-hydrodynamic simulations exploring the stability and dynamics of isothermal dusty gas columns in this regime. We describe our algorithm for full angle-dependent radiation transport based on the discontinuous Galerkin finite element method. For a range of near-Eddington fluxes, we show that the medium is unstable, producing convective-like motions in a turbulent atmosphere with a scale height significantly inflated compared to the gas pressure scale height and mass-weighted turbulent energy densities of $\sim 0.01-0.1$ of the midplane radiation energy density, corresponding to mass-weighted velocity dispersions of Mach number $\sim 0.5-2$. Extrapolation of our results to optical depths of $10^3$ implies maximum turbulent Mach numbers of $\sim20$. Comparing our results to galaxy-averaged observations, and subject to the approximations of our calculations, we find that radiation pressure does not contribute significantly to the effective supersonic pressure support in star-forming disks, which in general are substantially sub-Eddington. We further examine the time-averaged vertical density profiles in dynamical equilibrium and comment on implications for radiation-pressure-driven galactic winds.

Hydrogen-poor superluminous supernovae (SLSNe-I) have been predominantly found in low-metallicity, star-forming dwarf galaxies. Here we identify Gaia17biu/SN 2017egm as an SLSN-I occurring in a "normal" spiral galaxy (NGC 3191) in terms of stellar mass (several times 10^10 M_sun) and metallicity (roughly Solar). At redshift z=0.031, Gaia17biu is also the lowest redshift SLSN-I to date, and the absence of a larger population of SLSNe-I in dwarf galaxies of similar redshift suggests that metallicity is likely less important to the production of SLSNe-I than previously believed. With the smallest distance and highest apparent brightness for an SLSN-I, we are able to study Gaia17biu in unprecedented detail. Its pre-peak near-ultraviolet to optical color is similar to that of Gaia16apd and among the bluest observed for an SLSN-I while its peak luminosity (M_g = -21 mag) is substantially lower than Gaia16apd. Thanks to the high signal-to-noise ratios of our spectra, we identify several new spectroscopic features that may help to probe the properties of these enigmatic explosions. We detect polarization at the ~0.5% level that is not strongly dependent on wavelength, suggesting a modest, global departure from spherical symmetry. In addition, we put the tightest upper limit yet on the radio luminosity of an SLSN-I with <5.4x10^26 erg/s/Hz (at 10 GHz), which is almost a factor of 40 better than previous upper limits and one of the few measured at an early stage in the evolution of an SLSN-I. This limit largely rules out an association of this SLSNe-I with known populations of gamma-ray burst (GRB) like central engines.

We present the first spatially and spectrally resolved image of the molecular outflow in the western nucleus of Arp\,220. The outflow, seen in HCN~(1--0) by ALMA, is compact and collimated, with extension $\lesssim$ 120\,pc. Bipolar morphology emerges along the minor axis of the disk, with redshifted and blueshifted components reaching maximum inclination-corrected velocity of $\sim\,\pm$\,840\,km\,s$^{-1}$. The outflow is also seen in CO and continuum emission, the latter implying that it carries significant dust. We estimate a total mass in the outflow of $\geqslant$\,10$^{6}$\,M$_{\odot}$, a dynamical time of $\sim$\,10$^{5}$\,yr, and mass outflow rates of $\geqslant55$\,M$_{\odot}$\,yr$^{-1}$ and $\geqslant\,15$\,M$_{\odot}$\,yr$^{-1}$ for the northern and southern lobes, respectively. Possible driving mechanisms include supernovae energy and momentum transfer, radiation pressure feedback, and a central AGN. The latter could explain the collimated morphology of the HCN outflow, however we need more complex theoretical models, including contribution from supernovae and AGN, to pinpoint the driving mechanism of this outflow.

We present observations of the rise and peak of the Type IIn supernova SN 2017hcc/ATLAS17lsn obtained by the All-Sky Automated Survey for Supernovae (ASAS-SN) and Swift UVOT. The light curve of SN 2017hcc/ATLAS17lsn peaks at $V\simeq 13.7$ mag, which from the estimated redshift of the host galaxy ($z=0.0168$, $D\simeq 73$ Mpc) implies an absolute peak magnitude $M_{V,peak} \simeq -20.7$ mag. The near-UV to optical spectral energy distribution of SN 2017hcc/ATLAS17lsn from Swift UVOT is consistent with a hot, but cooling blackbody with $\rm T_{bb}\simeq 16500$ K on Oct. 28.4 and $\rm T_{bb} \simeq 11700$ K on Nov. 19.6. The estimated peak bolometric luminosity $L_{bol, peak}\simeq 1.3\times 10^{44}$ erg s$^{-1}$ makes SN2017hcc/ATLAS17lsn one of the most luminous Type IIn supernovae studied to date. From the bolometric light curve we constrain the risetime to be $\sim 27$ days and the total radiated energy of the event to date is $4\times 10^{50}$ erg.

This catalog summarizes information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) and all other bright ($m_{peak}\leq17$), spectroscopically confirmed supernovae discovered in 2016. We then gather the near-IR through UV magnitudes of all host galaxies and the offsets of the supernovae from the centers of their hosts from public databases. We illustrate the results using a sample that now totals 668 supernovae discovered since 2014 May 1, including the supernovae from our previous catalogs, with type distributions closely matching those of the ideal magnitude limited sample from Li et al. (2011). This is the third of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team.

Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large Area Telescope have shown that many classical novae are accompanied by gigaelectronvolt gamma-ray emission. This emission likely originates from strong shocks, providing new insights into the properties of nova outflows and allowing them to be used as laboratories to study the unknown efficiency of particle acceleration in shocks. Here we report gamma-ray and optical observations of the Milky Way nova ASASSN-16ma, which is among the brightest novae ever detected in gamma-rays. The gamma-ray and optical light curves show a remarkable correlation, implying that the majority of the optical light comes from reprocessed emission from shocks rather than the white dwarf. The ratio of gamma-ray to optical flux in ASASSN-16ma directly constrains the acceleration efficiency of non-thermal particles to be ~0.005, favouring hadronic models for the gamma-ray emission. The need to accelerate particles up to energies exceeding 100 gigaelectronvolts provides compelling evidence for magnetic field amplification in the shocks.

Winds arising from galaxies, star clusters, and active galactic nuclei are crucial players in star and galaxy formation, but it has proven remarkably difficult to use observations of them to determine physical properties of interest, particularly mass fluxes. Much of the difficulty stems from a lack of a theory that links a physically-realistic model for winds' density, velocity, and covering factors to calculations of light emission and absorption. In this paper we provide such a model. We consider a wind launched from a turbulent region with a range of column densities, derive the differential acceleration of gas as a function of column density, and use this result to compute winds' absorption profiles, emission profiles, and emission intensity maps in both optically thin and optically thick species. The model is sufficiently simple that all required computations can be done analytically up to straightforward numerical integrals, rendering it suitable for the problem of deriving physical parameters by fitting models to observed data. We show that our model produces realistic absorption and emission profiles for some example cases, and argue that the most promising methods of deducing mass fluxes are based on combinations of absorption lines of different optical depths, or on combining absorption with measurements of molecular line emission. In the second paper in this series, we expand on these ideas by introducing a set of observational diagnostics that are significantly more robust that those commonly in use, and that can be used to obtain improved estimates of wind properties.

Following work by W\"unsch and collaborators, we investigate a self-enrichment scenario for second generation star formation in globular clusters wherein wind material from first generation massive stars rapidly radiatively cools. Radiative energy loss allows retention of fast winds within the central regions of clusters, where it fuels star formation. Secondary star formation occurs in $\sim3-5$ Myr, before supernovae, producing uniform iron abundances in both populations. We derive the critical criteria for radiative cooling of massive star winds and the second generation mass as a function of cluster mass, radius, and metallicity. We derive a critical condition on $M/R$, above which second generation star formation can occur. We speculate that above this threshold the strong decrease in the cluster wind energy and momentum allows ambient gas to remain from the cluster formation process. We reproduce large observed second generation fractions of $\sim30-80\%$ if wind material mixes with ambient gas. Importantly, the mass of ambient gas required is only of order the first generation's stellar mass. Second generation helium enrichment $\Delta Y$ is inversely proportional to mass fraction in the second generation; a large second generation can form with $\Delta Y\sim0.001-0.02$, while a small second generation can reach $\Delta Y\sim0.16$. Like other self-enrichment models for the second generation, we are not able to simultaneously account for both the full range of the Na-O anticorrelation and the second generation fraction.

We compare the structure of molecular gas at $40$ pc resolution to the ability of gas to form stars across the disk of the spiral galaxy M51. We break the PAWS survey into $370$ pc and $1.1$ kpc resolution elements, and within each we estimate the molecular gas depletion time ($\tau_{\rm Dep}^{\rm mol}$), the star formation efficiency per free fall time ($\epsilon_{\rm ff}$), and the mass-weighted cloud-scale (40 pc) properties of the molecular gas: surface density, $\Sigma$, line width, $\sigma$, and $b\equiv\Sigma/\sigma^2\propto\alpha_{\rm vir}^{-1}$, a parameter that traces the boundedness of the gas. We show that the cloud-scale surface density appears to be a reasonable proxy for mean volume density. Applying this, we find a typical star formation efficiency per free-fall time, $\epsilon_{ff} \left( \left< \Sigma_{40pc} \right> \right) \sim 0.3{-}0.36\%$, lower than adopted in many models and found for local clouds. More, the efficiency per free fall time anti-correlates with both $\Sigma$ and $\sigma$, in some tension with turbulent star formation models. The best predictor of the rate of star formation per unit gas mass in our analysis is $b \equiv \Sigma / \sigma^2$, tracing the strength of self gravity, with $\tau_{\rm Dep}^{\rm mol} \propto b^{-0.9}$. The sense of the correlation is that gas with stronger self-gravity (higher $b$) forms stars at a higher rate (low $\tau_{\rm Dep}^{\rm mol}$). The different regions of the galaxy mostly overlap in $\tau_{\rm Dep}^{\rm mol}$ as a function of $b$, so that low $b$ explains the surprisingly high $\tau_{\rm Dep}^{\rm mol}$ found towards the inner spiral arms found by by Meidt et al. (2013).

The All-Sky Automated Survey for Supernovae (ASAS-SN) is working towards imaging the entire visible sky every night to a depth of V~17 mag. The present data covers the sky and spans ~2-5~years with ~100-400 epochs of observation. The data should contain some ~1 million variable sources, and the ultimate goal is to have a database of these observations publicly accessible. We describe here a first step, a simple but unprecedented web interface https://asas-sn.osu.edu/ that provides an up to date aperture photometry light curve for any user-selected sky coordinate. Because the light curves are produced in real time, this web tool is relatively slow and can only be used for small samples of objects. However, it also imposes no selection bias on the part of the ASAS-SN team, allowing the user to obtain a light curve for any point on the celestial sphere. We present the tool, describe its capabilities, limitations, and known issues, and provide a few illustrative examples.

We present Very Large Array observations of the 33 GHz radio continuum emission from 22 local ultraluminous and luminous infrared (IR) galaxies (U/LIRGs). These observations have spatial (angular) resolutions of 30--720 pc (0.07"-0.67") in a part of the spectrum that is likely to be optically thin. This allows us to estimate the size of the energetically dominant regions. We find half-light radii from 30 pc to 1.7 kpc. The 33 GHz flux density correlates well with the IR emission, and we take these sizes as indicative of the size of the region that produces most of the energy. Combining our 33 GHz sizes with unresolved measurements, we estimate the IR luminosity and star formation rate per area, and the molecular gas surface and volume densities. These quantities span a wide range (4 dex) and include some of the highest values measured for any galaxy (e.g., $\mathrm{\Sigma_{SFR}^{33GHz} \leq 10^{4.1} M_{\odot} yr^{-1} kpc^{-2}}$). At least $13$ sources appear Compton thick ($\mathrm{N_{H}^{33GHz} \geq 10^{24} cm^{-2}}$). Consistent with previous work, contrasting these data with observations of normal disk galaxies suggests a nonlinear and likely multi-valued relation between SFR and molecular gas surface density, though this result depends on the adopted CO-to-H$_{2}$ conversion factor and the assumption that our 33 GHz sizes apply to the gas. 11 sources appear to exceed the luminosity surface density predicted for starbursts supported by radiation pressure and supernovae feedback, however we note the need for more detailed observations of the inner disk structure. U/LIRGs with higher surface brightness exhibit stronger [{\sc Cii}] 158$\mu$m deficits, consistent with the suggestion that high energy densities drive this phenomenon.

Multiple observational lines of evidence support a connection between hydrogen-poor superluminous supernovae (SLSNe) and long duration gamma-ray bursts (GRBs). Both events require a powerful central energy source, usually attributed to a millisecond magnetar or an accreting black hole. The GRB-SLSN link raises several theoretical questions: What distinguishes the engines responsible for these different phenomena? Can a single engine power both a GRB and a luminous SN in the same event? We propose a new unifying model for magnetar thermalization and jet formation: misalignment between the rotation (${\bf \Omega}$) and magnetic dipole (${\bf \mu}$) axes thermalizes a fraction of the spindown power by reconnection in the striped equatorial wind, providing a guaranteed source of "thermal" emission to power the supernova. The remaining un-thermalized power energizes a relativistic jet. In this picture, the GRB-SLSN dichotomy is directly linked to ${\bf \Omega \cdot \mu}$. We extend earlier work to show that even weak relativistic jets of luminosity $\sim10^{46}$ erg s$^{-1}$ can escape the expanding SN ejecta hours after the explosion, implying that escaping relativistic jets may accompany many SLSNe. We calculate the observational signature of these jets. We show that they may produce transient UV cocoon emission lasting a few hours when the jet breaks out of the ejecta surface. A longer-lived optical/UV signal may originate from a mildly-relativistic wind driven from the interface between the jet and the ejecta walls. This provides a new explanation for the secondary early-time maximum observed in some SLSNe light curves, such as LSQ14bdq. This scenario also predicts a population of GRB from on-axis jets with extremely long durations, potentially similar to the population of "jetted tidal disruption events", in coincidence with a small subset of SLSNe.

We investigate the properties of Type IIP supernovae that are dominantly powered by the rotational kinetic energy of the newly born neutron star. While the spin-down of a magnetar has previously been proposed as a viable energy source in the context of super-luminous supernovae, we show that a similar mechanism could produce both normal and peculiar Type IIP supernova light curves from red supergiant progenitors for a range of initial spin periods and equivalent dipole magnetic field strengths. Although the formation channel for such magnetars in a typical red supergiant progenitor is unknown, it is tantalizing that this proof of concept model is capable of producing ordinary Type IIP lightcurve properties, perhaps implying that rotation rate and magnetic field strength may play important roles in some ordinary looking Type IIP supernova explosions.

Efficient thermalization of overlapping supernovae within star-forming galaxies may produce a supernova-heated fluid that drives galactic winds. For fiducial assumptions about the timescale for cloud shredding from high-resolution simulations (which neglect magnetic fields) we show that cool clouds with temperature from $T_{c}\sim 10^{2}-10^{4}$ K seen in emission and absorption in galactic winds cannot be accelerated to observed velocities by the ram pressure of a hot wind. Taking into account both the radial structure of the hot flow and gravity, we show that this conclusion holds over a wide range of galaxy, cloud, and hot wind properties. This finding calls into question the prevailing picture whereby the cool atomic gas seen in galactic winds is entrained and accelerated by the hot flow. Given these difficulties with ram pressure acceleration, we discuss alternative models for the origin of high velocity cool gas outflows. Another possibility is that magnetic fields in cool clouds are sufficiently important that they prolong the cloud's life. For $T_{c}=10^{3}$\,K and $10^{4}$\,K clouds, we show that if conductive evaporation can be neglected, the cloud shredding timescale must be $\sim15$ and 5 times longer, respectively, than the values from hydrodynamical simulations in order for cool cloud velocities to reach those seen in observations.

Given its peak luminosity and early-time spectra, ASASSN-15lh was classified as the most luminous supernova (SN) ever discovered (Dong et al. 2016). Here we report a UV rebrightening of ASASSN-15lh observed with Swift during our follow-up campaign. The rebrightening began at $t \simeq$ 90 days (observer frame) after the primary peak and was followed by a $\sim 120$-day long plateau in the bolometric luminosity, before starting to fade again at $t\simeq 210$ days. ASASSN-15lh rebrightened in the Swift UV bands by $\Delta m_{UVW2} \simeq -1.75$ mag, $\Delta m_{UVM2} \simeq -1.25$ mag, and $\Delta m_{UVW1} \simeq -0.8$ mag, but did not rebrighten in the optical bands. Throughout its initial decline, subsequent rebrightening, and renewed decline, the spectra did not show evidence of interactions between the ejecta and circumstellar medium (CSM) such as narrow emission lines. There are hints of weak H$\alpha$ emission at late-times, but Margutti et al. (2016) have shown that it is narrow line emission consistent with star formation in the host nucleus. By fitting a blackbody we find that during the rebrightening the effective photospheric temperature increased from $T_{BB} \simeq 11000$ K to $T_{BB} \simeq 18000$ K. Over the $\sim$ 550 days since its detection, ASASSN-15lh has radiated $\sim 1.7 - 1.9 \times 10^{52}$ ergs. Although its physical nature remains uncertain, the evolution of ASASSN-15lh's photospheric radius, its radiated energy, and the implied event rate, are all more similar to those of H-poor superluminous supernovae (SLSNe-I) than to tidal disruption events (TDEs).

This manuscript presents information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) during 2015, its second full year of operations. The same information is presented for bright ($m_V\leq17$), spectroscopically confirmed supernovae discovered by other sources in 2015. As with the first ASAS-SN bright supernova catalog, we also present redshifts and near-UV through IR magnitudes for all supernova host galaxies in both samples. Combined with our previous catalog, this work comprises a complete catalog of 455 supernovae from multiple professional and amateur sources, allowing for population studies that were previously impossible. This is the second of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team.

We present 90 mas (37 pc) resolution ALMA imaging of Arp 220 in the CO (1-0) line and continuum at $\lambda = 2.6$ mm. The internal gas distribution and kinematics of both galactic nuclei are well-resolved for the first time. In the West nucleus, the major gas and dust emission extends out to 0.2\arcsec radius (74 pc); the central resolution element shows a strong peak in the dust emission but a factor 3 dip in the CO line emission. In this nucleus, the dust is apparently optically thick ($\tau_{\rm 2.6mm} \sim1$) at $\lambda = 2.6$ mm with a dust brightness temperature $\sim147$ K. The column of ISM at this nucleus is $\rm N_{H2} \geq 2\times10^{26}$ cm$^{-2}$, corresponding to $\sim$900 gr cm$^{-2}$. The East nucleus is more elongated with radial extent 0.3\arcsec or $\sim111$ pc. The derived kinematics of the nuclear disks provide a good fit to the line profiles, yielding the emissivity distributions, the rotation curves and velocity dispersions. In the West nucleus, there is evidence of a central Keplerian component requiring a central mass of $8\times10^8$ \msun. The intrinsic widths of the emission lines are $\Delta \rm v (FWHM)$ = 250 (West) and 120 (East) \kms. Given the very short dissipation timescales for turbulence ($\lesssim10^5$ yrs), we suggest that the line widths may be due to semi-coherent motions within the nuclear disks. The symmetry of the nuclear disk structures is impressive -- implying the merger timescale is significantly longer than the rotation period of the disks.

We use idealized 3D hydrodynamic simulations to study the dynamics and thermal structure of the circumgalactic medium (CGM). Our simulations quantify the role of cooling, stellar feedback driven galactic winds and cosmological gas accretion in setting the properties of the CGM in dark matter haloes ranging from $10^{11}$ to $10^{12}$ M$_\odot$. Our simulations support a conceptual picture in which the properties of the CGM, and the key physics governing it, change markedly near a critical halo mass of M$_{\rm crit} \approx 10^{11.5}$ M$_\odot$. As in calculations without stellar feedback, above M$_{\rm crit}$ halo gas is supported by thermal pressure created in the virial shock. The thermal properties at small radii are regulated by feedback triggered when $t_{\rm cool}/t_{\rm ff}\lesssim10$ in the hot gas. Below M$_{\rm crit}$, however, there is no thermally supported halo and self-regulation at $t_{\rm cool}/t_{\rm ff}\sim10$ does not apply. Instead, the gas is out of hydrostatic equilibrium and largely supported against gravity by bulk flows (turbulence and coherent inflow/outflow) arising from the interaction between cosmological gas inflow and outflowing galactic winds. In these lower mass haloes, the phase structure depends sensitively on the outflows' energy per unit mass and mass-loading, which may allow measurements of the CGM thermal state to constrain the nature of galactic winds. Our simulations account for some of the properties of the multiphase halo gas inferred from quasar absorption line observations, including the presence of significant mass at a wide range of temperatures, and the characteristic OVI and CIV column densities and kinematics. However, we underpredict the neutral hydrogen content of the $z\sim0$ CGM.

We explore heavy element nucleosynthesis in neutrino-driven winds from rapidly-rotating, strongly magnetized proto-neutron stars for which the magnetic dipole is aligned with the rotation axis, and the field is assumed to be a static force-free configuration. We process the proto-magnetar wind trajectories calculated by Vlasov et al 2014 through the r-process nuclear reaction network SkyNet using contemporary models for the evolution of the wind electron fraction during the proto-neutron star cooling phase. Although we do not find a successful second or third peak r-process for any rotation period P, we show that proto-magnetars with P around 1-5 ms produce heavy element abundance distributions that extend to higher nuclear mass number than from otherwise equivalent spherical winds (with the mass fractions of some elements enhanced by factors of 100-1000). The heaviest elements are synthesized by outflows emerging along flux tubes which graze the closed zone and pass near the equatorial plane outside the light cylinder. Due to dependence of the nucleosynthesis pattern on the magnetic field strength and rotation rate of the proto-neutron star, natural variations in these quantities between core collapse events could contribute to the observed diversity of the abundances of weak r-process nuclei in metal-poor stars. Further diversity, including possibly even a successful third-peak r-process, could be achieved for misaligned rotators with non-zero magnetic inclination with respect to the rotation axis. If proto-magnetars are central engines for GRBs, their relativistic jets should contain a high mass fraction of heavy nuclei of characteristic average mass number A of order 100, providing a possible source for ultra-high energy cosmic rays comprised of heavy nuclei with an energy spectrum that extends beyond the nominal GZK cut-off for protons or iron nuclei.

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 ($\Delta 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 $\mu = 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}}$.

Star-forming galaxies produce gamma-rays primarily via pion production, resulting from inelastic collisions between cosmic ray protons and the interstellar medium (ISM). The dense ISM and high star formation rates of luminous and ultra-luminous infrared galaxies (LIRGs and ULIRGs) imply that they should be strong gamma-ray emitters, but so far only two LIRGs have been detected. Theoretical models for their emission depend on the unknown fraction of cosmic ray protons that escape these galaxies before interacting. We analyze Fermi-LAT data for 82 of the brightest IRAS LIRGs and ULIRGs. We examine each system individually and carry out a stacking analysis to constrain their gamma-ray fluxes. We report the detection of the nearest ULIRG Arp 220 (~4.6sigma). We observe a gamma-ray flux (0.8--100 GeV) of 2.4e-10 phot cm^-2 s^-1 with photon index of 2.23 (8.2e10^41 ergs s^-1 at 77 Mpc) We also derive upper limits for the stacked LIRGs and ULIRGs. The gamma-ray luminosity of Arp~220 and the stacked upper limits agree with calorimetric predictions for dense star-forming galaxies. With the detection of Arp 220, we extend the gamma-ray--IR luminosity correlation to the high luminosity regime with log(L_[0.1-100 GeV]) = 1.25log(L_[8-1000]) + 26.7 as well as the gamma-ray--radio continuum luminosity correlation with log(L_[0.1-100 GeV]) = 1.22log(L_[1.4 GHz]) + 13.3. The current survey of Fermi-LAT is on the verge of detecting more LIRGs/ULIRGs in the local universe, and we expect even more detections with deeper Fermi-LAT observations or the next generation of gamma-ray detectors.

The transport of gas towards the centre of galaxies is critical for black hole feeding and, indirectly, it can control active galactic nucleus (AGN) feedback. We have quantified the molecular gas inflow in the central R<1kpc of M51 to be 1 Msun/yr, using a new gravitational torque map and the molecular gas traced by the PdBI Arcsecond Whirlpool Survey (PAWS). The nuclear stellar bar is responsible for this gas inflow. We have also used torque profiles to estimate the location of dynamical resonances, suggesting a corotation for the bar at R~20", and a corotation for the spiral at R~100". We demonstrate how important it is to correct 3.6um images for dust emission in order to compute gravitational torques, and we carefully examine further sources of uncertainty. Our observational measurement of gas inflow can be compared with nuclear molecular outflow rates and provide useful constraints for numerical simulations.

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.

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.

Dynamical scattering of binaries and triple systems of stars, planets, and compact objects may produce highly inclined triple systems subject to Kozai-Lidov (KL) oscillations, potentially leading to collisions, mergers, Type Ia supernovae, and other phenomena. We present the results of more than 400 million gravitational scattering experiments of binary-binary, triple-single, and triple-binary scattering. We compute the cross sections for all possible outcomes and explore their dependencies on incoming velocity, mass, semi-major axis, and eccentricity, including analytic fits and discussion of the velocity dependence. For the production of new triple systems by scattering we find that compact triples are preferred, with ratios of outer to inner semi-major axes of ~few--100, flat or quasi-thermal eccentricity distributions, and flat distributions in cosine of the mutual inclination. Dynamically formed triples are thus subject to strong KL oscillations, the "eccentric Kozai mechanism," and non-secular effects. For single and binary flyby encounters with triple systems, we compute the cumulative cross section for changes to the mutual inclination, eccentricity, and semi-major axis ratio. We apply these results to scattering events in the field, open clusters, and globular clusters, and explore the implications for Type Ia supernovae via collisions and mergers, stellar collisions, and the lifetime and dynamical isolation of triple systems undergoing KL oscillations. An Appendix provides an analysis of the velocity dependence of the collision cross section in binary-single scattering.

We show that the turbulent gas in the star-forming regions of galaxies is unstable to wind formation via momentum deposition by radiation pressure or other momentum sources like supernova explosions, even if the system is below the average Eddington limit. This conclusion follows from the fact that the critical momentum injection rate per unit mass for unbinding gas from a self-gravitating system is proportional to the gas surface density and that a turbulent medium presents a broad distribution of column densities to the sources. For an average Eddington ratio of <Gamma>~0.1 and for turbulent Mach numbers >30, we find that ~1% of the gas is ejected per dynamical timescale at velocities larger than the local escape velocity. Because of the lognormal shape of the surface density distribution, the mass loss rate is highly sensitive to the average Eddington ratio, reaching 20-40% of the gas mass per dynamical time for <Gamma>=1. Using this model we find a large scatter in the mass-loading factor for star-forming galaxies, ranging from 0.001-10, but with significant uncertainties. Implications for the efficiency of star formation in giant molecular clouds are highlighted. For radiation pressure feedback alone, we find an increasing star formation efficiency as a function of initial gas surface density. Uncertainties are discussed.

The physical origin of high velocity cool gas seen in galactic winds remains unknown. Following Wang (1995), we argue that radiative cooling in initially hot thermally-driven outflows can produce fast neutral atomic and photoionized cool gas. The inevitability of adiabatic cooling from the flow's initial 10^7-10^8K temperature and the shape of the cooling function for T<10^7K imply that outflows with hot gas mass-loss rate relative to star formation rate of beta=Mdot_hot/Mdot_star > 0.5 cool radiatively on scales ranging from the size of the energy injection region to tens of kpc. We highlight the beta and star formation rate surface density dependence of the column density, emission measure, radiative efficiency, and velocity. At r_cool, the gas produces X-ray and then UV/optical line emission with a total power bounded by 10^{-2} L_star if the flow is powered by steady-state star formation with luminosity L_star. The wind is thermally unstable at r_cool, potentially leading to a multi-phase medium. Cooled winds decelerate significantly in the extended gravitational potential of galaxies. The cool gas precipitated from hot outflows may explain its prevalence in galactic halos. We forward a picture of winds whereby cool clouds are initially accelerated by the ram pressure of the hot flow, but are rapidly shredded by hydrodynamical instabilities, thereby increasing beta, seeding radiative and thermal instability, and cool gas rebirth. If the cooled wind shocks as it sweeps up the circumgalactic medium, its cooling time is short, thus depositing cool gas far out into the halo. Finally, conduction can dominate energy transport in low-beta hot winds, leading to flatter temperature profiles than otherwise expected, potentially consistent with X-ray observations of some starbursts.

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.

We describe and execute a novel approach to observationally estimate the lifetimes of giant molecular clouds (GMCs). We focus on the cloud population between the two main spiral arms in M51 (the inter-arm region) where cloud destruction via shear and star formation feedback dominates over formation processes. By monitoring the change in GMC number densities and properties from one side of the inter-arm to the other, we estimate the lifetime as a fraction of the inter-arm travel time. We find that GMC lifetimes in M51's inter-arm are finite and short, 20 to 30 Myr. Such short lifetimes suggest that cloud evolution is influenced by environment, in which processes can disrupt GMCs after a few free-fall times. Over most of the region under investigation shear appears to regulate the lifetime. As the shear timescale increases with galactocentric radius, we expect cloud destruction to switch primarily to star formation feedback at larger radii. We identify a transition from shear- to feedback-dominated disruption through a change in the behavior of the GMC number density. The signature suggests that shear is more efficient at completely dispersing clouds, whereas feedback transforms the population, e.g. by fragmenting high mass clouds into lower mass pieces. Compared to the characteristic timescale for molecular hydrogen in M51, our short lifetimes suggest that gas can remain molecular while clouds disperse and reassemble. We propose that galaxy dynamics regulates the cycling of molecular material from diffuse to bound (and ultimately star-forming) objects, contributing to long observed molecular depletion times in normal disk galaxies. We also speculate that, in more extreme environments such as elliptical galaxies and concentrated galaxy centers, star formation can be suppressed when the shear timescale becomes so short that some clouds can not survive to collapse and form stars.

The demographics of the field of Astronomy, and the gender balance in particular, is an important active area of investigation. A piece of information missing from the discussion is the gender breakdown of the applicant pool for faculty positions. For a sample of 35 tenure-track faculty positions at 25 research universities advertised over the last few years in astronomy and astrophysics, I find that the ratio of female applicants to the total number of applicants is ~0.2, with little dispersion and with no strong dependence on the total number of applicants. Some discussion is provided in the context of the fraction of women at the graduate student, postdoctoral researcher, and assistant professor levels, but strong conclusions are not possible given the limitations of the study. Current and future faculty search committees will likely be interested to compare their numbers to this distribution to decide whether or not they could be doing more to attract an applicant pool that is representative of the community.

It is typically assumed that radiation pressure driven winds are accelerated to an asymptotic velocity of V ~ v_esc, where v_esc is the escape velocity from the central source. We note that this is not the case for dusty shells and clouds. Instead, if the shell or cloud is initially optically-thick to the UV emission from the source of luminosity L, then there is a significant boost in V that reflects the integral of the momentum absorbed as it is accelerated. For shells reaching a generalized Eddington limit, we show that V ~ (4R_UV L/M_sh c)^1/2, in both point-mass and isothermal-sphere potentials, where R_UV is the radius where the shell becomes optically-thin to UV photons, and M_sh is the mass of the shell. The asymptotic velocity significantly exceeds v_esc for typical parameters, and can explain the ~1000-2000km/s outflows observed from rapidly star-forming galaxies and active galactic nuclei if the surrounding halo has low gas density. Similarly fast outflows from massive stars can be accelerated on few - 10^3 yr timescales. These results carry over to clouds that subtend only a small fraction of the solid angle from the source of radiation and that expand as a consequence of their internal sound speed. We further consider the dynamics of shells that sweep up a dense circumstellar or circumgalactic medium. We calculate the "momentum ratio" Mdot v/(L/c) in the shell limit and show that it can only significantly exceed ~2 if the effective optical depth of the shell to re-radiated FIR photons is much larger than unity. We discuss simple prescriptions for the properties of galactic outflows for use in large-scale cosmological simulations. We also briefly discuss applications to the dusty ejection episodes of massive stars, the disruption of giant molecular clouds, and AGN.

If the neutrino luminosity from the proto-neutron star formed during a massive star core collapse exceeds a critical threshold, a supernova (SN) results. Using spherical quasi-static evolutionary sequences for hundreds of progenitors over a range of metallicities, we study how the explosion threshold maps onto observables, including the fraction of successful explosions, the neutron star (NS) and black hole (BH) mass functions, the explosion energies (E_SN) and nickel yields (M_Ni), and their mutual correlations. Successful explosions are intertwined with failures in a complex pattern that is not simply related to initial progenitor mass or compactness. We predict that progenitors with initial masses of 15 +/- 1, 19 +/- 1, and 21-26 M_Sun are most likely to form BHs, that the BH formation probability is non-zero at solar-metallicity and increases significantly at low metallicity, and that low luminosity, low Ni-yield SNe come from progenitors close to success/failure interfaces. We qualitatively reproduce the observed E_SN-M_Ni correlation, we predict a correlation between the mean and width of the NS mass and E_SN distributions, and that the means of the NS and BH mass distributions are correlated. We show that the observed mean NS mass of ~1.33 M_Sun implies that the successful explosion fraction is higher than 0.35. Overall, we show that the neutrino mechanism can in principle explain the observed properties of SNe and their compact objects. We argue that the rugged landscape of progenitors and outcomes mandates that SN theory should focus on reproducing the wide ranging distributions of observed SN properties.

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emission surrounding both nuclei and conclude that is mostly synchrotron in nature. The spatial distributions of radio emission in both nuclei are well described by exponential profiles. These have deconvolved half-light radii of 51 and 35 pc for the eastern and western nuclei, and they match the number density profile of radio supernovae observed with very long baseline interferometry. This similarity might be due to the fast cooling of cosmic rays electrons caused by the presence of a strong (~ mG) magnetic field in this system. We estimate high luminosity surface densities of $\mathrm{\Sigma_{IR} \sim 4.2^{+1.6}_{-0.7} \times 10^{13}}$ (east) and $\mathrm{\sim 9.7^{+3.7}_{-2.4} \times 10^{13}~(west)~L_{\odot}~kpc^{-2}}$, and star formation rate surface densities of $\mathrm{\Sigma_{SFR} \sim 10^{3.7\pm0.1}}$ (east) and $\mathrm{\sim 10^{4.1\pm0.1}~(west)~M_{\odot}~yr^{-1}~kpc^{-2}}$. These values, especially for the western nucleus are, to our knowledge, the highest luminosity and star formation rate surface densities measured for any star-forming system. Despite these high values, the nuclei lie below the dusty Eddington limit in which radiation pressure is balanced only by self-gravity. The small measured sizes also imply that the nuclei of Arp 220 are only transparent in the frequency range ~ 5 to 350 GHz. Our results offer no clear evidence that an active galactic nucleus dominates the emission from either nucleus at 33 GHz.

The accretion induced collapse (AIC) of a white dwarf to a neutron star has long been suggested as a natural theoretical outcome in stellar evolution, but there has never been a direct detection of such an event. This is not surprising since the small amount of radioactive nickel synthesized ($\sim10^{-3}\,M_\odot$) implies a relatively dim optical transient. Here we argue that a particularly strong signature of an AIC would occur for an oxygen-neon-magnesium (ONeMg) white dwarf accreting from a star that is experiencing Roche-lobe overflow as it becomes a red giant. In such cases, the $\sim10^{50}\,{\rm erg}$ explosion from the AIC collides with and shock-heats the surface of the extended companion, creating an X-ray flash lasting $\sim1\,{\rm hr}$ followed by an optical signature that peaks at an absolute magnitude of $\sim -16$ to $-18$ and lasts for a few days to a week. These events would be especially striking in old stellar environments where hydrogen-rich supernova-like, transients would not normally be expected. Although the rate of such events is not currently known, we describe observing strategies that could be utilized with high cadence surveys that should either detect these events or place strong constraints on their rates.

We calculate the steady-state properties of neutrino-driven winds from strongly magnetized, rotating proto-neutron stars (`proto-magnetars') under the assumption that the outflow geometry is set by the force-free magnetic field of an aligned dipole. Our goal is to assess proto-magnetars as sites of r-process nucleosynthesis and gamma-ray burst engines. One dimensional solutions calculated along flux tubes corresponding to different polar field lines are stitched together to determine the global properties of the flow at a given neutrino luminosity and rotation period. Proto-magnetars with rotation periods of P~2-5 ms are shown to produce outflows more favorable for the production of third-peak r-process nuclei due to their much shorter expansion times through the seed nucleus formation region, yet only moderately lower entropies, as compared to normal spherical PNS winds. Proto-magnetars with moderately rapid birth periods P~3-5 ms may thus represent a promising Galactic r-process site which is compatible with a variety of other observations, including the recent discovery of possible magnetar-powered supernovae in metal poor galaxies. We also confirm previous results that the outflows from proto-magnetars with P~1-2 ms can achieve maximum Lorentz factors Gamma ~ 100-1000 in the range necessary to power gamma-ray bursts (GRBs). The implications of GRB jets with a heavy nuclei-dominated composition as sources of ultra-high energy cosmic rays are also addressed.

The optical and near-IR emission from some classes of supernovae (SNe), including Type IIn and possibly some super-luminous SNe, is likely powered by a collision between the SN ejecta and dense circumstellar material (CSM). We argue that for a range of CSM masses and their radii, a collisionless shock can form, allowing for efficient cosmic-ray (CR) acceleration. We show that pp collisions between these newly accelerated CRs and the CSM leads to not only gamma rays but also secondary electrons and positrons that radiate synchrotron photons in the high-frequency radio bands. Our estimates imply that various facilities including the Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA) may observe such SNe at Gpc distances by followup observations in months-to-years, although the detectability strongly depends on the CSM density as well as observed frequency. Detecting this signal would give us a unique probe of CR acceleration at early times, and even non-detections can put interesting limits on the possibility of CR ion acceleration. Following our previous work, we also show that GeV gamma rays can escape from the system without severe attenuation, encouraging point-source and stacking analyses with Fermi. We provide recipes for diagnosing interaction-powered SN scenario with multi-messenger (neutrino and gamma-ray) observations.

Kozai-Lidov (KL) oscillations can accelerate compact object mergers via gravitational wave (GW) radiation by driving the inner binaries of hierarchical triples to high eccentricities. We perform direct three-body integrations of high mass ratio compact object triple systems using Fewbody including post-Newtonian terms. We find that the inner binary undergoes rapid eccentricity oscillations (REOs) on the timescale of the outer orbital period which drive it to higher eccentricities than secular theory would otherwise predict, resulting in substantially reduced merger times. For a uniform distribution of tertiary eccentricity ($e_2$), ~40% of systems merge within ~1-2 eccentric KL timescales whereas secular theory predicts that only ~20% of such systems merge that rapidly. This discrepancy becomes especially pronounced at low $e_2$, with secular theory overpredicting the merger time by many orders of magnitude. We show that a non-negligible fraction of systems have eccentricity > 0.8 when they merge, in contrast to predictions from secular theory. Our results are applicable to high mass ratio triple systems containing black holes or neutron stars. In objects in which tidal effects are important, such as white dwarfs, stars, and planets, REOs can reduce the tidal circularization timescale by an order of magnitude and bring the components of the inner binary into closer orbits than would be possible in the secular approximation.

Galactic superwinds may be driven by very hot outflows generated by overlapping supernovae within the host galaxy. We use the Chevalier & Clegg (CC85) wind model and the observed correlation between X-ray luminosities of galaxies and their SFRs to constrain the mass loss rates (\dot{M}_hot) across a wide range of star formation rates (SFRs), from dwarf starbursts to ultra-luminous infrared galaxies. We show that for fixed thermalization efficiency and mass loading rate, the X-ray luminosity of the hot wind scales as L_X ~ SFR^2, significantly steeper than is observed for star-forming galaxies: L_X ~ SFR. Using this difference we constrain the mass-loading and thermalization efficiency of hot galactic winds. For reasonable values of the thermalization efficiency (<~ 1) and for SFR >~ 10 M_sun/yr we find that \dot{M}_hot/SFR <~ 1, significantly lower than required by integrated constraints on the efficiency of stellar feedback in galaxies, and potentially too low to explain observations of winds from rapidly star-forming galaxies. In addition, we highlight the fact that heavily mass-loaded winds cannot be described by the adiabatic CC85 model because they become strongly radiative.

The kinematic complexity and the favorable position of M51 on the sky make this galaxy an ideal target to test different theories of spiral arm dynamics. Taking advantage of the new high resolution PdBI Arcsecond Whirlpool Survey (PAWS) data, we undertake a detailed kinematic study of M51 to characterize and quantify the origin and nature of the non-circular motions. Using a tilted-ring analysis supported by several other archival datasets we update the estimation of M51's position angle (PA=(173 +/- 3) deg) and inclination (i=(22 +/- 5) deg). Harmonic decomposition of the high resolution (40 pc) CO velocity field shows the first kinematic evidence of an m=3 wave in the inner disk of M51 with a corotation at R(CR,m=3)=1.1 +/- 0.1 kpc and a pattern speed of Omega_p(m=3) = 140 km/(s kpc). This mode seems to be excited by the nuclear bar, while the beat frequencies generated by the coupling between the m=3 mode and the main spiral structure confirm its density-wave nature. We observe also a signature of an m=1 mode that is likely responsible for the lopsidedness of M51 at small and large radii. We provide a simple method to estimate the radial variation of the amplitude of the spiral perturbation (Vsp) attributed to the different modes. The main spiral arm structure has <Vsp>=50-70 km/s, while the streaming velocity associated with the m=1 and m=3 modes is, in general, 2 times lower. Our joint analysis of HI and CO velocity fields at low and high spatial resolution reveals that the atomic and molecular gas phases respond differently to the spiral perturbation due to their different vertical distribution and emission morphology.

Using data from the PdBI Arcsecond Whirlpool Survey (PAWS), we have generated the largest extragalactic Giant Molecular Cloud (GMC) catalog to date, containing 1,507 individual objects. GMCs in the inner M51 disk account for only 54% of the total 12CO(1-0) luminosity of the survey, but on average they exhibit physical properties similar to Galactic GMCs. We do not find a strong correlation between the GMC size and velocity dispersion, and a simple virial analysis suggests that 30% of GMCs in M51 are unbound. We have analyzed the GMC properties within seven dynamically-motivated galactic environments, finding that GMCs in the spiral arms and in the central region are brighter and have higher velocity dispersions than inter-arm clouds. Globally, the GMC mass distribution does not follow a simple power law shape. Instead, we find that the shape of the mass distribution varies with galactic environment: the distribution is steeper in inter-arm region than in the spiral arms, and exhibits a sharp truncation at high masses for the nuclear bar region. We propose that the observed environmental variations in the GMC properties and mass distributions are a consequence of the combined action of large-scale dynamical processes and feedback from high mass star formation. We describe some challenges of using existing GMC identification techniques for decomposing the 12CO(1-0) emission in molecule-rich environments, such as M51's inner disk.

Rotating and magnetized protoneutron stars (PNSs) may drive relativistic magneto-centrifugally accelerated winds as they cool immediately after core collapse. The wind fluid near the star is composed of neutrons and protons, and the neutrons become relativistic while collisionally coupled with the ions. Here, we argue that the neutrons in the flow eventually undergo inelastic collisions around the termination shock inside the stellar material, producing ~0.1-1 GeV neutrinos, without relying on cosmic-ray acceleration mechanisms. Even higher-energy neutrinos may be produced via particle acceleration mechanisms. We show that PINGU and Hyper-Kamiokande can detect such neutrinos from nearby core-collapse supernovae, by reducing the atmospheric neutrino background via coincident detection of MeV neutrinos or gravitational waves and optical observations. Detection of these GeV and/or higher-energy neutrinos would provide important clues to the physics of magnetic acceleration, nucleosynthesis, the relation between supernovae and gamma-ray bursts, and the properties of newly born neutron stars.

To study the role of protosellar jets and outflows in the time evolution of the parent cores and the protostars, the astronomical community needs a large enough data base of infrared images of protostars at the highest spatial resolution possible, to reveal the details of their morphology. Spitzer provides unprecedented sensitivity in the infrared to study both the jet and outflow features, however its spatial resolution is limited by its 0.85m mirror. Here we use a high resolution deconvolution algorithm, "HiRes", to improve the visualization of spatial morphology by enhancing resolution (to sub-arcsecond levels in the IRAC bands) and removing the contaminating sidelobes from bright sources in a sample of 89 protostellar objects. These reprocessed images are useful to detect: (i) wide angle outflow seen in scattered light; (ii) morphological details of H2 emission in jets and bow shocks; and (iii) compact features in MIPS 24 micron images as protostar/ disk and atomic/ionic line emissions associated with the jets. The HiRes fits image data of such a large homogeneous sample presented here will be useful to the community in studying these protostellar objects. To illustrate the utility of this HiRes sample, we show how the opening angle of the wide angle outflows in 31 sources, all observed in the HiRes processed Spitzer images, correlates with age. Our data suggest a power law fit to opening angle versus age with an exponent of ~0.32 and 0.02, respectively for ages less than 8000 yr and greater than 8000 yr.

White dwarf-white dwarf (WD-WD) mergers may lead to type Ia supernovae events. Thompson (2011) suggested that many such binaries are produced in hierarchical triple systems. The tertiary induces eccentricity oscillations in the inner binary via the Kozai-Lidov mechanism, driving the binary to high eccentricities, and significantly reducing the gravitational wave merger timescale (T_GW) over a broad range of parameter space. Here, we investigate the role of tidal forces in these systems. We show that tidal effects are important in the regime of moderately high initial relative inclination between the inner binary and the outer tertiary. For 85 < i_0 < 90 degrees (prograde) and 97 < i_0 < 102 degrees (retrograde), tides combine with GW radiation to dramatically decrease T_GW. In the regime of high inclinations between 91 < i_0 < 96 degrees, the inner binary likely suffers a direct collision, as in the work of Katz & Dong (2012) and tidal effects do not play an important role.

The energy output (per logarithmic interval of particle energies) of Cosmic Rays (CRs) with energies $10{\rm GeV}\lesssim\varepsilon_p\lesssim100{\rm GeV}$ is $\sim 10^{47}\rm erg$ per solar mass of star$-$formation, based on the CR production rate in the Milky Way and in starburst galaxies, implying a generation rate of $\varepsilon_p^2Q\sim 10^{45}\rm erg~Mpc^{-3}~yr^{-1}$ in the local universe. It is only $\sim 10$ times larger than the output, $\varepsilon_p^2 Q=0.5\pm0.2\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$, of Ultra High Energy CRs (UHECRs) at energies $10^{10.5}{\rm GeV}<\varepsilon_p<10^{12}\rm GeV$ (obtained assuming they are mostly protons), which in turn is comparable to the lower limit of $\varepsilon_p^2 Q\ge 0.5\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$ of high energy CRs with $10^6{\rm GeV}\lesssim\varepsilon_p\lesssim 10^{8}\rm GeV$ implied by the saturation of the Waxman-Bahcall bound by the neutrino excess recently discovered by IceCube. These similarities are consistent with a flat production spectrum, $\varepsilon_p^2 Q\sim \text{const}$ for CRs at all observed energies. If a flat production spectrum is generated by our galaxy, the observed CR flux in the range $10^{6.5}-10^{9.5}\rm GeV$, above the "knee", is suppressed compared to lower energies due to propagation effects rather than acceleration upper limits. As suggested by Parizot and Aublin, the most exciting possibility is that cosmic rays at all energies are emitted from a single type of (unknown) sources, which can not be supernova remnants.

We compare the properties of giant molecular clouds (GMCs) in M51 identified by the Plateau de Bure Interferometer Whirlpool Arcsecond Survey (PAWS) with GMCs identified in wide-field, high resolution surveys of CO emission in M33 and the Large Magellanic Cloud (LMC). We find that GMCs in M51 are larger, brighter and have higher velocity dispersions relative to their size than equivalent structures in M33 and the LMC. These differences imply that there are genuine variations in the average mass surface density of the different GMC populations. To explain this, we propose that the pressure in the interstellar medium surrounding the GMCs plays a role in regulating their density and velocity dispersion. We find no evidence for a correlation between size and linewidth in any of M51, M33 or the LMC when the CO emission is decomposed into GMCs, although moderately robust correlations are apparent when regions of contiguous CO emission (with no size limitation) are used. Our work demonstrates that observational bias remains an important obstacle to the identification and study of extragalactic GMC populations using CO emission, especially in molecule-rich galactic environments.

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

The PdBI (Plateau de Bure Interferometer) Arcsecond Whirlpool Survey (PAWS) has mapped the molecular gas in the central ~9kpc of M51 in its 12CO(1-0) line emission at cloud-scale resolution of ~40pc using both IRAM telescopes. We utilize this dataset to quantitatively characterize the relation of molecular gas (or CO emission) to other tracers of the interstellar medium (ISM), star formation and stellar populations of varying ages. Using 2-dimensional maps, a polar cross-correlation technique and pixel-by-pixel diagrams, we find: (a) that (as expected) the distribution of the molecular gas can be linked to different components of the gravitational potential, (b) evidence for a physical link between CO line emission and radio continuum that seems not to be caused by massive stars, but rather depend on the gas density, (c) a close spatial relation between the PAH and molecular gas emission, but no predictive power of PAH emission for the molecular gas mass,(d) that the I-H color map is an excellent predictor of the distribution (and to a lesser degree the brightness) of CO emission, and (e) that the impact of massive (UV-intense) young star-forming regions on the bulk of the molecular gas in central ~9kpc can not be significant due to a complex spatial relation between molecular gas and star-forming regions that ranges from co-spatial to spatially offset to absent. The last point, in particular, highlights the importance of galactic environment -- and thus the underlying gravitational potential -- for the distribution of molecular gas and star formation.

We study the orbital evolution of hierarchical quadruple systems composed of two binaries on a long mutual orbit, where each binary acts as a Kozai-Lidov (KL) perturber on the other. We find that the coupling between the two binaries qualitatively changes the behavior of their KL cycles. The binaries can experience coherent eccentricity oscillations as well as excursions to very high eccentricity that occur over a much larger fraction of the parameter space than in triple systems. For a ratio of outer to inner semi-major axes of 10 to 20, about 30 to 50% of equal-mass quadruples reach eccentricity 1-e < 10^{-3} in one of the binaries. This is about 4 to 12 times more than for triples with equivalent parameters. Orbital "flips" and collisions without previous tidal interaction are similarly enhanced in quadruples relative to triples. We argue that the frequency of evolutionary paths influenced by KL cycles is comparable in the triple and quadruple populations even though field quadruples are a factor of ~5 less frequent than triples. Additionally, quadruples might be a non-negligible source of triples and provide fundamentally new evolutionary outcomes involving close binaries, mergers, collisions, and associated transients, which occur without any fine tuning of parameters. Finally, we study the perturbations to a planetary orbit due to a distant binary and we find that the fraction of orbital flips is a factor of 3 to 4 higher than for the corresponding triple system given our fiducial parameters with implications for hot Jupiters and star-planet collisions.

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

[abridged] Radiation pressure on dust grains may be an important mechanism in driving winds in a wide variety of astrophysical systems. However, the efficiency of the coupling between the radiation field and the dusty gas is poorly understood in environments characterized by high optical depths. We present a series of idealized numerical experiments, performed with the radiation-hydrodynamic code ORION, in which we study the dynamics of such winds and quantify their properties. We find that, after wind acceleration begins, radiation Rayleigh-Taylor instability forces the gas into a configuration that reduces the rate of momentum transfer from the radiation field to the gas by a factor ~ 10 - 100 compared to an estimate based on the optical depth at the base of the atmosphere; instead, the rate of momentum transfer from a driving radiation field of luminosity L to the gas is roughly L/c multiplied by one plus half the optical depth evaluated using the photospheric temperature, which is far smaller than the optical depth one would obtain using the interior temperature. When we apply our results to conditions appropriate to ULIRGs and star clusters, we find that the asymptotic wind momentum flux from such objects should not significantly exceed that carried by the direct radiation field, L/c. This result constrains the expected mass loss rates from systems that exceed the Eddington limit to be of order the so-called "single-scattering" limit, and not significantly higher. We present an approximate fitting formula for the rate of momentum transfer from radiation to dusty gas through which it passes, which is suitable for implementation in sub-grid models of galaxy formation. Finally, we provide a first map of the column density distribution of gas in a radiatively-driven wind as a function of velocity, and velocity dispersion.

[Abdriged] The origin of the diffuse hard X-ray (2 - 10 keV) emission from starburst galaxies is a long-standing problem. We suggest that synchrotron emission of 10 - 100 TeV electrons and positrons (e+/-) can contribute to this emission, because starbursts have strong magnetic fields. We consider three sources of e+/- at these energies: (1) primary electrons directly accelerated by supernova remnants; (2) pionic secondary e+/- created by inelastic collisions between CR protons and gas nuclei in the dense ISMs of starbursts; (3) pair e+/- produced between the interactions between 10 - 100 TeV gamma-rays and the intense far-infrared (FIR) radiation fields of starbursts. We create one-zone steady-state models of the CR population in the Galactic Center (R <= 112 pc), NGC 253, M82, and Arp 220's nuclei, assuming a power law injection spectrum for electrons and protons. We compare these models to extant radio and GeV and TeV gamma-ray data for these starbursts, and calculate the diffuse synchrotron X-ray and Inverse Compton (IC) luminosities of these starbursts. If the primary electron spectrum extends to ~PeV energies and has a proton/electron injection ratio similar to the Galactic value, we find that synchrotron contributes 2 - 20% of their unresolved, diffuse hard X-ray emission. Inverse Compton emission is likewise a minority of the unresolved X-ray emission in these starbursts, from 0.1% in the Galactic Center to 10% in Arp 220's nuclei. We also model generic starbursts, including submillimeter galaxies, in the context of the FIR--X-ray relation, finding that up to 2% in the densest starbursts with our fiducial assumptions. Neutrino and TeV gamma-ray data can further constrain the synchrotron X-ray emission of starbursts. Our models do not constrain hard synchrotron X-ray emission from any additional hard components of primary e+/- from sources like pulsars in starbursts.

We investigate star formation rate (SFR) calibrations in light of recent developments in the modeling of stellar rotation. Using new published non-rotating and rotating stellar tracks, we study the integrated properties of synthetic stellar populations and find that the UV to SFR calibration for the rotating stellar population is 30% smaller than for the non-rotating stellar population, and 40% smaller for the Halpha to SFR calibration. These reductions translate to smaller SFR estimates made from observed UV and Halpha luminosities. Using the UV and Halpha fluxes of a sample of ~300 local galaxies, we derive a total (i.e., sky-coverage corrected) SFR within 11 Mpc of 120-170 Msun/yr and 80-130 Msun/yr for the non-rotating and rotating estimators, respectively. Independently, the number of core-collapse supernovae discovered in the same volume requires a total SFR of 270^{+110}_{-80} Msun/yr, suggesting a mild tension with the SFR estimates made with rotating calibrations. More generally, when compared with the estimated SFR, the local supernova discoveries strongly constrain any physical effects that might increase the energy output of massive stars, including, but not limited to, stellar rotation. The cosmic SFR and cosmic supernova rate data on the other hand show the opposite trend, with the cosmic SFR almost double that inferred from cosmic supernovae, constraining a significant decrease in the energy output of massive stars. Together, these lines of evidence suggest that the true SFR calibration factors cannot be too far from their canonical values.

Radiation pressure on dust grains may be an important physical mechanism driving galaxy-wide superwinds in rapidly star-forming galaxies. We calculate the combined dust and gas Eddington ratio (Gamma) for the archetypal superwind of M82. By combining archival GALEX data, a standard dust model, Monte Carlo dust scattering calculations, and the Herschel map of the dust surface density distribution, the observed FUV/NUV surface brightness in the outflow constrains both the total UV luminosity escaping from the starburst along its minor axis (L_*, UV) and the flux-mean opacity, thus allowing a calculation of Gamma. We find that L_(*, UV) ~ 1-6*10^42 ergs s^-1, ~2-12 times greater than the UV luminosity observed from our line of sight. On a scale of 1-3 kpc above the plane of M82, we find that Gamma ~ 0.01 - 0.06. On smaller scales (~0.25-0.5 kpc), where the enclosed mass decreases, our calculation of L_(*, UV) implies that Gamma ~ 0.1 with factor of few uncertainties. Within the starburst itself, we estimate the single-scattering Eddington ratio to be of order unity. Thus, although radiation pressure is weak compared to gravity on kpc scales above the plane of M82, it may yet be important in launching the observed outflow. We discuss the primary uncertainties in our calculation, the sensitivity of Gamma to the dust grain size distribution, and the time evolution of the wind following M82's recent starburst episodes.

We use the high spatial and spectral resolution of the PAWS CO(1-0) survey of the inner 9 kpc of the iconic spiral galaxy M51 to examine the effect of gas streaming motions on the star-forming properties of individual GMCs. We compare our view of gas flows in M51 -- which arise due to departures from axi-symmetry in the gravitational potential (i.e. the nuclear bar and spiral arms) -- with the global pattern of star formation as traced by Halpha and 24\mu m emission. We find that the dynamical environment of GMCs strongly affects their ability to form stars, in the sense that GMCs situated in regions with large streaming motions can be stabilized, while similarly massive GMCs in regions without streaming go on to efficiently form stars. We argue that this is the result of reduced surface pressure felt by clouds embedded in an ambient medium undergoing large streaming motions, which prevents collapse. Indeed, the variation in gas depletion time expected based on the observed streaming motions throughout the disk of M51 quantitatively agrees with the variation in observed gas depletion time scale. The example of M51 shows that streaming motions, triggered by gravitational instabilities in the form of bars and spiral arms, can alter the star formation law; this can explain the variation in gas depletion time among galaxies with different masses and morphologies. In particular, we can explain the long gas depletion times in spiral galaxies compared to dwarf galaxies and starbursts. We suggest that adding a dynamical pressure term to the canonical free-fall time produces a single star formation law that can be applied to all star-forming regions and galaxies, across cosmic time.

We present the data of the Plateau de Bure Arcsecond Whirlpool Survey (PAWS), a high spatial and spectral resolution 12CO(1-0) line survey of the inner 10 x 6 kpc of the M51 system, and the first wide-field imaging of molecular gas in a star-forming spiral galaxy with resolution matched to the typical size of Giant Molecular Clouds (40 pc). We describe the observation, reduction, and combination of the Plateau de Bure Interferometer (PdBI) and IRAM-30m "short spacing" data. The final data cube attains 1.1"-resolution over the 270" x 170" field of view, with sensitivity to all spatial scales from the combination of PdBI and IRAM-30m data, and brightness sensitivity of 0.4 K (1 sigma) in each 5 km/s-wide channel map. We find a CO-luminosity of 9.10^8 K km/s pc^2, corresponding to a molecular gas mass of 4.10^9 Msol for a standard CO-to-H2 conversion factor. Unexpectedly, we find that a large fraction, (50 +/- 10)%, of this emission arises mostly from spatial scales larger than 36" or about 1.3 kpc. Through a series of tests, we demonstrate that this extended emission does not result from a processing artifact. We discuss its origin in light of the stellar component, the 12CO/13CO ratio, and the difference between the kinematics and structure of the PdBI-only and hybrid synthesis (PdBI + IRAM-30m) images. The extended emission is consistent with a thick, diffuse disk of molecular gas with a typical scale height of about 200 pc, substructured in unresolved filaments which fills about 0.1% of the volume.

We analyse the distribution of CO brightness temperature and integrated intensity in M51 at ~40 pc resolution using new CO data from the Plateau de Bure Arcsecond Whirlpool Survey (PAWS). We present probability distribution functions (PDFs) of the CO emission within the PAWS field, which covers the inner 11 x 7 kpc of M51. We find variations in the shape of CO PDFs within different M51 environments, and between M51 and M33 and the Large Magellanic Cloud (LMC). Globally, the PDFs for the inner disk of M51 can be represented by narrow lognormal functions that cover 1 to 2 orders of magnitude in CO brightness and integrated intensity. The PDFs for M33 and the LMC are narrower and peak at lower CO intensities. However, the CO PDFs for different dynamical environments within the PAWS field depart from the shape of the global distribution. The PDFs for the interarm region are approximately lognormal, but in the spiral arms and central region of M51, they exhibit diverse shapes with a significant excess of bright CO emission. The observed environmental dependence of the shape of the CO PDFs is qualitatively consistent with changes that would be expected if molecular gas in the spiral arms has a larger range of average densities, gas temperatures and velocity fluctuations, though further work is required to disentangle the importance of large-scale dynamical effects versus star formation feedback in regulating these properties. We show that the shape of the CO PDFs for different M51 environments is only weakly related to global properties of the CO emission, but is strongly correlated with some properties of the local giant molecular cloud (GMC) and young stellar cluster populations. For galaxies with strong spiral structure such as M51, our results indicate that galactic-scale dynamical processes play a significant role in the formation and evolution of GMCs and stellar clusters.(abridged)

Over a broad range of initial inclinations and eccentricities an appreciable fraction of hierarchical triple star systems with similar masses are essentially unaffected by the Kozai-Lidov mechanism (KM) until the primary in the central binary evolves into a compact object. Once it does, it may be much less massive than the other components in the ternary, enabling the "eccentric Kozai mechanism (EKM):" the mutual inclination between the inner and outer binary can flip signs driving the inner binary to very high eccentricity, leading to a close binary or collision. We demonstrate this "Mass-loss Induced Eccentric Kozai" (MIEK) mechanism by considering an example system and defining an ad-hoc minimal separation between the inner two members at which tidal affects become important. For fixed initial masses and semi-major axes, but uniform distributions of eccentricity and cosine of the mutual inclination, ~10% of systems interact tidally or collide while the primary is on the MS due to the KM or EKM. Those affected by the EKM are not captured by earlier quadrupole-order secular calculations. We show that fully ~30% of systems interact tidally or collide for the first time as the primary swells to AU scales, mostly as a result of the KM. Finally, ~2% of systems interact tidally or collide for the first time after the primary sheds most of its mass and becomes a WD, mostly as a result of the MIEK mechanism. These findings motivate a more detailed study of mass-loss in triple systems and the formation of close NS/WD-MS and NS/WD-NS/WD binaries without an initial common envelope phase.

We use a high-temperature chemical network to derive the molecular abundances in axisymmetric accretion disk models around active galactic nuclei (AGNs) within 100 pc using simple radial and vertical density and temperature distributions motivated by more detailed physical models. We explore the effects of X-ray irradiation and cosmic ray ionization on the spatial distribution of the molecular abundances of CO, CN, CS, HCN, HCO+, HC3N, C2H, and c-C3H2 using a variety of plausible disk structures. These simple models have molecular regions with a layer of X-ray dominated regions, a midplane without the strong influence of X-rays, and a high-temperature region in the inner portion with moderate X-ray flux where families of polyynes (C$_{\rm n}$H$_{2}$) and cyanopolyynes can be enhanced. For the high midplane density disks we explore, we find that cosmic rays produced by supernovae do not significantly affect the regions unless the star formation efficiency significantly exceeds that of the Milky Way. We highlight molecular abundance observations and ratios that may distinguish among theoretical models of the density distribution in AGN disks. Finally, we assess the importance of the shock crossing time and the accretion time relative to the formation time for various chemical species. Vertical column densities are tabulated for a number of molecular species at both the characteristic shock crossing time and steady state. Although we do not attempt to fit any particular system or set of observations, we discuss our models and results in the context of the nearby AGN NGC 1068.

I show that the temperature of nuclear star clusters, starburst clusters in M82, compact high-z galaxies, and some globular clusters of the Galaxy likely exceeded the ice line temperature (T_Ice ~ 150-170 K) during formation for a time comparable to the planet formation timescale. The protoplanetary disks within these systems will thus not have an ice line, decreasing the total material available for building protoplanetary embryos, inhibiting the formation of gas- and ice-giants if they form by core accretion, and prohibiting habitability. Planet formation by gravitational instability is similarly suppressed because Toomre's Q > 1 in all but the most massive disks. I discuss these results in the context of the observed lack of planets in 47 Tuc. I predict that a similar search for planets in the globular cluster NGC 6366 ([Fe/H] = -0.82) should yield detections, whereas (counterintuitively) the relatively metal-rich globular clusters NGC 6440, 6441, and 6388 should be devoid of giant planets. The characteristic stellar surface density above which T_Ice is exceeded in star clusters is ~6 x 10^3 M_sun/pc^2 f_{dg, MW}^{-1/2}, where f_{dg, MW} is the dust-to-gas ratio of the embedding material, normalized to the Milky Way value. Simple estimates suggest that ~5 - 50% of the stars in the universe formed in an environment exceeding this surface density. Caveats and uncertainties are detailed.

Kerzendorf et al. (2012) recently reported the startling discovery of a metal-poor ([Fe/H]=-1 +/- 0.4) A-type star near the center of the Tycho supernova remnant. We propose two possible explanations. In the first, Tycho B is a blue straggler, formed from the merger of a close K- or G-type binary system, which was previously in a quadruple system with the binary that produced SN 1572. Both binaries were likely brought to tidal contact by Kozai-Lidov oscillations acting in concert with tidal friction. Analogous progenitor systems may include CzeV343, VW LMi, and KIC 4247791. In the second, Tycho B is the surviving tertiary component of a triple system, which was also likely affected by Kozai-Lidov oscillations. Rates are briefly discussed. Problems with each evolutionary scenario are presented. Finally, a chance alignment between Tycho B and the supernova remnant is not excluded.

[abridged] The pressure exerted by the radiation of young stars may be an important feedback mechanism in forming star clusters and the disks of starburst galaxies. However, there is great uncertainty in how efficiently radiation couples to matter in these high optical depth environments. In particular, it is unclear what levels of turbulence the radiation can produce, and whether the infrared radiation trapped by the dust opacity can give rise to heavily mass-loaded winds. In this paper we report a series of two-dimensional flux-limited diffusion radiation-hydrodynamics calculations performed with the code ORION in which we drive strong radiation fluxes through columns of dusty matter confined by gravity. We consider both systems where the radiation flux is sub-Eddington throughout the gas column, and where it is super-Eddington at the midplane but sub-Eddington in the atmosphere. In the latter, we find that the radiation-matter interaction gives rise to radiation-driven Rayleigh-Taylor instability, which drives supersonic turbulence at a level sufficient to fully explain the turbulence seen in Galactic protocluster gas clouds, and to make a non-trivial contribution to the turbulence observed in starburst galaxy disks. However, the instability also produces a channel structure in which the radiation-matter interaction is reduced because the radiation field is not fully trapped. For astrophysical parameters relevant to forming star clusters and starburst galaxies, we find that this effect reduces the net momentum deposition rate in the dusty gas by a factor of ~2-6 compared to simple analytic estimates, and that in steady state the Eddington ratio reaches unity and there are no strong winds. We provide an approximation formula, appropriate for implementation in analytic models and non-radiative simulations, for the force exerted by the infrared radiation field in this regime.

We present a K-band spectroscopic study of the stellar and gas kinematics, mass distribution, and stellar populations of the archetypical starburst galaxy M82. Our results are based on a single spectrum at a position angle of 67.5 degrees through the K-band nucleus. We used the CO stellar absorption band head at 2.29 {\mu}m (CO_2.29) to measure the rotation curve out to nearly 4 kpc radius on both the eastern and western sides of the galaxy. Our data show that the rotation curve is flat from 1 - 4 kpc. This stands in sharp contrast to some previous studies, which have interpreted H I and CO emission-line position-velocity diagrams as evidence for a declining rotation curve. The kinematics of the Br\gamma, H_2, and He I emission lines are consistent with, although characterized by slightly higher velocities than, the stellar kinematics. We derived M82's mass distribution from our stellar kinematic measurements and estimate its total dynamical mass is ~10^10 Msun. We measured the equivalent width of CO_2.29 (W_2.29) as a function of distance from the center of the galaxy to investigate the spatial extent of the red supergiant (RSG) population. The variation in W_2.29 with radius clearly shows that RSGs dominate the light inside 500 pc radius. M82's superwind is likely launched from this region, where we estimate the enclosed mass is <= 2 x 10^9 Msun.

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 \sigma, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

In the seconds after collapse of a massive star, the newborn proto-neutron star (PNS) radiates neutrinos of all flavors. The absorption of electron-type neutrinos below the radius of the stalled shockwave may drive explosions (the "neutrino mechanism"). Because the heating rate is proportional to the square of neutrino energy, flavor conversion of mu and tau neutrinos to electron-type neutrinos via collective neutrino oscillations (CnuO) may in principle increase the heating rate and drive explosions. In order to assess the potential importance of CnuO for the shock revival, we solve the steady-state boundary value problem of spherically-symmetric accretion between the PNS surface (r_nu) and the shock (r_S), including a scheme for flavor conversion via CnuO. For a given r_nu, PNS mass (M), accretion rate (Mdot), and assumed values of the neutrino energies from the PNS, we calculate the critical neutrino luminosity above which accretion is impossible and explosion results. We show that CnuO can decrease the critical luminosity by a factor of at most ~1.5, but only if the flavor conversion is fully completed inside r_S and if there is no matter suppression. The magnitude of the effect depends on the model parameters (M, Mdot, and r_nu) through the shock radius and the physical scale for flavor conversion. We quantify these dependencies and find that CnuO could lower the critical luminosity only for small M and Mdot, and large r_nu. However, for these parameter values CnuO are suppressed due to matter effects. By quantifying the importance of CnuO and matter suppression at the critical neutrino luminosity for explosion, we show in agreement with previous studies that CnuO are unlikely to affect the neutrino mechanism of core-collapse supernovae significantly.

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

(Abridged) We study large-scale winds driven from uniformly bright self-gravitating discs radiating near the Eddington limit. We show that the ratio of the radiation pressure force to the gravitational force increases with height above the disc surface to a maximum of twice the value of the ratio at the disc surface. Thus, uniformly bright self-gravitating discs radiating at the Eddington limit are fundamentally unstable to driving large-scale winds. These results contrast with the spherically symmetric case, where super-Eddington luminosities are required for wind formation. We apply this theory to galactic winds from rapidly star-forming galaxies that approach the Eddington limit for dust. For hydrodynamically coupled gas and dust, we find that the asymptotic velocity of the wind is v_\infty ~ 1.5 v_rot and that v_\infty SFR^{0.36}, where v_rot is the disc rotation velocity and SFR is the star formation rate, both of which are in agreement with observations. However, these results of the model neglect the gravitational potential of the surrounding dark matter halo and an old passive stellar bulge or extended disc, which act to decrease v_\infty. A more realistic treatment shows that the flow can either be unbound, or bound, forming a "fountain flow" with a typical turning timescale of t_turn ~ 0.1-1 Gyr. We provide quantitative criteria and scaling relations for assessing whether or not a rapidly star-forming galaxy of given properties can drive unbound flows via the mechanism described in this paper. Importantly, we note that because t_turn is longer than the star formation timescale in the rapidly star-forming galaxies and ULIRGs for which our theory is most applicable, if rapidly star-forming galaxies are selected as such, they may be observed to have strong outflows, even though their winds are eventually bound on large scales.

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.

(Abridged) Neutrino heating may drive core-collapse supernova explosions. Although it is known that the stalled accretion shock turns into explosion when the neutrino luminosity from the collapsed core exceeds a critical value (L_crit) (the "neutrino mechanism"), the physics of L_crit, as well as its dependence on the properties of the proto-neutron star (PNS) and changes to the microphysics has never been systematically explored. We solve the one-dimensional steady-state accretion problem between the PNS surface and the accretion shock. We quantify the deep connection between the solution space of steady-state accretion flows with bounding shocks and the neutrino mechanism. We show that there is a maximum, critical sound speed above which it is impossible to maintain accretion with a standoff shock, because the shock jump conditions cannot be satisfied. The physics of this critical sound speed is general and does not depend on a specific heating mechanism. For the simple model of pressure-less free-fall onto a shock bounding an isothermal accretion flow with sound speed c_T, we show that if c_T^2/v_escape^2 > 3/16 explosion results. We generalize this result to the more complete supernova problem, showing explicitly that the same physics determines L_crit. We find that the critical condition for explosion can be written as c_S^2/v_escape^2 = 0.19, where c_S is the adiabatic sound speed. This "antesonic" condition describes L_crit over a broad range in accretion rate and microphysics. We show that the addition of the accretion luminosity (L_acc) reduces L_crit non-trivially. As in previous work, we find that L_crit is always significantly higher than the maximum possible value of L_acc. Finally, we provide evidence that the reduction in L_crit seen in recent multi-dimensional simulations results from a reduction in the efficiency of cooling, rather than an increase in the heating rate.

Approximately 1/4-1/2 of short duration Gamma-Ray Bursts are followed by variable X-ray emission lasting ~ 100 s with a fluence comparable or exceeding that of the initial burst itself. The long duration and significant energy of this `extended emission'(EE) poses a major challenge to the standard binary neutron star merger model. metzger08 recently proposed that the EE is powered by the spin-down of a strongly magnetized neutron star. However, the effects of surrounding material on the magnetar outflow have not yet been considered. Here we present time-dependent axisymmetric relativistic magnetohydrodynamic simulations of the interaction of the relativistic proto-magnetar wind with a surrounding 10^-1 10^-3 M_\odot envelope, which represents material ejected during the merger; the supernova following AIC; or via outflows from the accretion disk. The collision between the relativistic magnetar wind and the expanding ejecta produces a magnetized nebula inside the ejecta. A strong toroidal magnetic field builds up in the nebula, which drives a bipolar jet out through the ejecta, similar to the magnetar model developed in the case of long duration GRBs. We quantify the `break-out' time and opening angle of the jet theta_j as a function of the wind energy flux dot{E} and ejecta mass M_ej. We show that dot{E} and theta_j are inversely correlated, such that the beaming-corrected (isotropic) luminosity of the jet is primarily a function of M_ej. Both variability arguments, and the lower limit on the power of magnetar outflows capable of producing bright emission, suggest that the true opening angle of the magnetar jet must be relatively large. The model thus predicts a class of events for which the EE is observable with no associated short GRB. These may appear as long-duration GRBs or X-Ray Flashes, which may be detected by future all-sky X-ray survey missions.

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

A new class of core-collapse supernovae (SNe) has been discovered in recent years by optical/infrared surveys; these SNe suggest the presence of one or more extremely dense (~10^5-10^11 cm^-3) shells of circumstellar material (CSM) on 10^2-10^4 AU scales. We consider the collisions of the SN ejecta with these massive CSM shells as potential cosmic-ray (CR) accelerators. If ~10% of the SN energy goes into CRs, multi-TeV neutrinos and/or GeV-TeV gamma rays almost simultaneous with the optical/infrared light curves are detectable for SNe at <20-30 Mpc. A new type of coordinated multi-messenger search for such transients of duration ~1 - 10 months is required; these may give important clues to the physical origin of such SNe and to CR acceleration mechanisms.

White dwarf-white dwarf (WD-WD) and neutron star-neutron star (NS-NS) mergers may produce Type Ia supernovae and gamma-ray bursts (GRBs), respectively. A general problem is how to produce binaries with semi-major axes small enough to in significantly less than the Hubble time (t_H), and thus accommodate the observation that these events closely follow episodes of star formation. I explore the possibility that such systems are not binaries at all, but actually coeval, or dynamical formed, triple systems. The tertiary induces Kozai oscillations in the inner binary, driving it to high eccentricity, and reducing its gravitational wave (GW) merger timescale. This effect significantly increases the allowed range of binary period P such that the merger time is t_merge < t_H. In principle, Chandrasekhar-mass binaries with P ~ 300 days can merge in < t_H if they contain a prograde solar-mass tertiary at high enough inclination. For retrograde tertiaries, the maximum P such that t_merge < t_H is yet larger. In contrast, P < 0.3 days is required in the absence of a tertiary. I discuss implications of these findings for the production of transients formed via compact object binary mergers. Based on the statistics of solar-type binaries, I argue that many such binaries should be in triple systems affected by the Kozai resonance. If true, expectations for the mHz GW signal from individual sources, the diffuse background, and the foreground for GW experiments like LISA are modified. This work motivates future studies of triples systems of A, B,and O stars, and new types of searches for WD-WD binaries in triple systems.

Long duration Gamma-Ray Bursts (GRBs) originate from the core collapse of massive stars, but the identity of the central engine remains elusive. Previous work has shown that rapidly spinning, strongly magnetized proto-neutron stars (`millisecond proto-magnetars') produce outflows with energies, timescales, and magnetizations sigma_0 (maximum Lorentz factor) that are consistent with those required to produce long GRBs. Here we extend this work in order to construct a self-consistent model that directly connects the properties of the central engine to the observed prompt emission. Just after the launch of the supernova shock, a wind heated by neutrinos is driven from the proto-magnetar. The outflow is collimated into a bipolar jet by its interaction with the star. As the magnetar cools, the wind becomes ultra-relativistic and Poynting-flux dominated (sigma_0 >> 1) on a timescale comparable to that required for the jet to clear a cavity through the star. Although the site and mechanism of the prompt emission are debated, we calculate the emission predicted by two models: magnetic dissipation and internal shocks. Our results favor the magnetic dissipation model in part because it predicts a relatively constant `Band' spectral peak energy E_peak with time during the GRB. The jet baryon loading decreases abruptly when the neutron star becomes transparent to neutrinos at t ~ 10-100 seconds. Jets with ultra-high magnetization cannot effectively accelerate and dissipate their energy, suggesting this transition ends the prompt emission and may explain the steep decay phase that follows. We assess several phenomena potentially related to magnetar birth, including low luminosity GRBs, thermal-rich GRBs/X-ray Flashes, very luminous supernovae, and short duration GRBs with extended emission.

The GeV and TeV emission from M82 and NGC 253 observed by Fermi, HESS, and VERITAS constrains the physics of cosmic rays (CRs) in these dense starbursts. We argue that the gamma rays are predominantly hadronic in origin, as expected by previous studies. The measured fluxes imply that pionic losses are efficient for CR protons in both galaxies: we show that a fraction F_cal ~ 0.2 - 0.4 of the energy injected in high energy primary CR protons is lost to inelastic proton-proton collisions (pion production) before escape, producing gamma rays, neutrinos, and secondary electrons and positrons. We discuss the factor ~2 uncertainties in this estimate, including supernova rate and leptonic contributions to the GeV-TeV emission. We argue that gamma-ray data on ULIRGs like Arp 220 can test whether M82 and NGC 253 are truly calorimetric, and we present upper limits on Arp 220 from the Fermi data. We show that the observed ratio of the GeV to GHz fluxes of the starbursts suggests that non-synchrotron cooling processes are important for cooling the CR electron/positron population. We briefly reconsider previous predictions in light of the gamma-ray detections, including the starburst contribution to the gamma-ray background and CR energy densities. Finally, as a guide for future studies, we list the brightest star-forming galaxies on the sky and present updated predictions for their gamma-ray and neutrino fluxes.

We derive an upper limit to the energy of nuclei accelerated via the Fermi mechanism in any relativistic shockwave, driven by any astrophysical engine. This bound is accessible to current and upcoming ultra-high energy neutrino experiments. Detection of a single neutrino with energy above the upper limit would exclude all sites of shock acceleration, and imply physics beyond the Standard Model. We comment on the possibility that relativistic flows launched by supermassive black hole mergers are the source of the observed ultra-high energy cosmic rays.

We evaluate radiation pressure from starlight on dust as a feedback mechanism in star-forming galaxies by comparing the luminosity and flux of star-forming systems to the dust Eddington limit. The linear LFIR--L'HCN correlation provides evidence that galaxies may be regulated by radiation pressure feedback. We show that star-forming galaxies approach but do not dramatically exceed Eddington, but many systems are significantly below Eddington, perhaps due to the "intermittency" of star formation. Better constraints on the dust-to-gas ratio and the CO- and HCN-to-H2 conversion factors are needed to make a definitive assessment of radiation pressure as a feedback mechanism.

Radiation pressure from the absorption and scattering of starlight by dust grains may be an important feedback mechanism in regulating star-forming galaxies. We compile data from the literature on star clusters, star-forming subregions, normal star-forming galaxies, and starbursts to assess the importance of radiation pressure on dust as a feedback mechanism, by comparing the luminosity and flux of these systems to their dust Eddington limit. This exercise motivates a novel interpretation of the Schmidt Law, the LIR-L'CO correlation, and the LIR-L'HCN correlation. In particular, the linear LIR-L'HCN correlation is a natural prediction of radiation pressure regulated star formation. Overall, we find that the Eddington limit sets a hard upper bound to the luminosity of any star-forming region. Importantly, however, many normal star-forming galaxies have luminosities significantly below the Eddington limit. We explore several explanations for this discrepancy, especially the role of "intermittency" in normal spirals - the tendency for only a small number of subregions within a galaxy to be actively forming stars at any moment because of the time-dependence of the feedback process and the luminosity evolution of the stellar population. If radiation pressure regulates star formation in dense gas, then the gas depletion timescale is 6 Myr, in good agreement with observations of the densest starbursts. Finally, we highlight the importance of observational uncertainties - namely, the dust-to-gas ratio and the CO-H2 and HCN-H2 conversion factors - that must be understood before a definitive assessment of radiation pressure as a feedback mechanism in star-forming galaxies.

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.

The millisecond proto-magnetar model for the central engine of long-duration gamma-ray bursts is briefly reviewed. Limitations and uncertainties in the model are highlighted. A short discussion of the maximum energy, maximum duration, radiative efficiency, jet formation mechanism, late-time energy injection, and (non-)association with supernovae of millisecond magnetar-powered GRBs is provided.

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.

We construct one-zone steady-state models of cosmic ray (CR) injection, cooling, and escape over the entire dynamic range of the FIR-radio correlation (FRC), from normal galaxies to starbursts, over the redshift interval 0 <= z <= 10. Normal galaxies with low star-formation rates become radio-faint at high z, because Inverse Compton (IC) losses off the CMB cool CR electrons and positrons rapidly, suppressing their nonthermal radio emission. However, we find that this effect occurs at higher redshifts than previously expected, because escape, bremsstrahlung, ionization, and starlight IC losses act to counter this effect and preserve the radio luminosity of galaxies. The radio dimming of star-forming galaxies at high z is not just a simple competition between magnetic field energy density and the CMB energy density; the CMB must also compete with every other loss process. We predict relations for the critical redshift when radio emission is significantly suppressed compared to the z ~ 0 FRC as a function of star-formation rate per unit area. Additionally, we provide a quantitative explanation for the relative radio brightness of some high-z submillimeter galaxies. We show that at fixed star formation rate surface density, galaxies with larger CR scale heights are radio bright with respect to the FRC, because of weaker bremsstrahlung and ionization losses compared to compact starbursts. We predict that these "puffy starbursts" should have steeper radio spectra than compact galaxies with the same star-formation rate surface density. We find that radio bright submillimeter galaxies alone cannot explain the excess radio emission reported by ARCADE2, but they may significantly enhance the diffuse radio background with respect to a naive application of the z ~ 0 FRC.

(Abridged) The far-infrared (FIR) and radio luminosities of star-forming galaxies are linearly correlated over a very wide range in star formation rate, from normal spirals like the Milky Way to the most intense starbursts. Using one-zone models of cosmic ray (CR) injection, cooling, and escape in star-forming galaxies, we attempt to reproduce the observed FIR-radio correlation over its entire span. We show that ~2% of the kinetic energy from supernova explosions must go into primary CR electrons and that ~10 - 20% must go into primary CR protons. Secondary electrons and positrons are likely comparable to or dominate primary electrons in dense starburst galaxies. We discuss the implications of our models for the magnetic field strengths of starbursts, the detectability of starbursts by Fermi, and cosmic ray feedback. Overall, our models indicate that both CR protons and electrons escape from low surface density galaxies, but lose most of their energy before escaping dense starbursts. The FIR-radio correlation is caused by a combination of the efficient cooling of CR electrons (calorimetry) in starbursts and a conspiracy of several factors. For lower surface density galaxies, the decreasing radio emission caused by CR escape is balanced by the decreasing FIR emission caused by the low effective UV dust opacity. In starbursts, bremsstrahlung, ionization, and Inverse Compton cooling decrease the radio emission, but they are countered by secondary electrons/positrons and the decreasing critical synchrotron frequency, which both increase the radio emission. Our conclusions hold for a broad range of variations on our fiducial model.

Galactic outflows of low ionization, cool gas are ubiquitous in local starburst galaxies, and in the majority of galaxies at high redshift. How these cool outflows arise is still in question. Hot gas from supernovae has long been suspected as the primary driver, but this mechanism suffers from its tendency to destroy the cool gas as the latter is accelerated. We propose a modification of the supernova scenario that overcomes this difficulty. Star formation is observed to take place in clusters; in a given galaxy, the bulk of the star formation is found in the ~20 most massive clusters. We show that, for L* galaxies, the radiation pressure from clusters with M>10^6 M_sun is able to expel the surrounding gas at velocities in excess of the circular velocity of the disk galaxy. This cool gas can travel above the galactic disk in less than 2 Myr, well before any supernovae erupt in the driving cluster. Once above the disk, the cool outflowing gas is exposed to radiation, and supernovae induced hot gas outflows, from other clusters in the disk, which drive it to distances of several tens to hundreds of kpc. Because the radiatively driven clouds grow in size as they travel, and because the hot gas is more dilute at large distance, the clouds are less subject to destruction if they do eventually encounter hot gas. Therefore, unlike wind driven clouds, radiatively driven clouds can survive to distances ~50 kpc. We identify these cluster-driven winds with large-scale galactic outflows. Another implication of our model is that only starburst galaxies, where massive clusters reside, are able to drive winds cold outflows on galactic scales via this mechanism. We find that the critical star formation rates above which large scale cool outflows will be launched to be ~0.1 M_sun/yr/kpc^2, which is in good agreement with observations.

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.

We argue that the recent groundbreaking discovery by Badenes et al. (2009) of a nearby (~50 pc) white dwarf-neutron star (or black hole) binary (SDSS 1257+5428) with a merger timescale ~500 Myr implies that such systems are common; we estimate that there are of order 10^6 in the Galaxy. Although subject to large uncertainties, the nominal derived merger rate is ~5 x 10^-4 per yr in the Milky Way, just ~3-6 and ~20-40 times less than the Type Ia and core-collapse supernova (SN) rates, respectively. This implies that the merger rate is ~0.5-1 x 10^4 per Gpc^3 per yr in the local universe, ~5000-10000 times more than the observed (beaming-uncorrected) long-duration gamma-ray burst (GRB) rate. We estimate the lower limit on the rate in the Galaxy to be >2.5 x 10^-5 per yr at 95% confidence. We briefly discuss the implications of this finding for the census of long- and short-duration GRBs and their progenitors, the frequency of tight binary companions to Type Ib/c SN progenitors, the origin of ultra-high energy cosmic rays (UHECRs), the formation of rapidly rotating neutron stars and ~2-3 M_sun black holes, the census of faint Ia-like SNe, as well as for upcoming and current transient surveys (e.g., LOSS, PTF, LSST), and for high- (LIGO) and low-frequency (LISA) gravitational wave searches.

Radiative diffusion damps acoustic modes at large comoving wavenumber (k) before decoupling (``Silk damping''). In a simple WKB analysis, neglecting moments of the temperature distribution beyond the quadrupole, damping appears in the acoustic mode as a term of order ik^2/(taudot) where taudot is the scattering rate per unit conformal time. Although the Jeans instability is stabilized on scales smaller than the adiabatic Jeans length, I show that the medium is linearly unstable to first order in (1/taudot) to a slow diffusive mode. At large comoving wavenumber, the characteristic growth rate becomes independent of spatial scale and constant: (t_{KH}a)^-1 ~ (128 pi G/9 kappa_T c)(rho_m/rho_b), where "a" is the scale factor, rho_m and rho_b are the matter and baryon energy density, respectively, and kappa_T is the Thomson opacity. This is the characteristic timescale for a fluid parcel to radiate away its thermal energy content at the Eddington limit, analogous to the Kelvin-Helmholz (KH) time for a massive star or the Salpeter time for black hole growth. Although this mode grows at all times prior to decoupling and on scales smaller than the horizon, the growth time is long, about 100 times the age of the universe at decoupling. Thus, it modifies the density and temperature perturbations on small scales only at the percent level. The physics of this mode is already accounted for in the popular codes CMBFAST and CAMB, but is typically neglected in analytic studies of the growth of primordial perturbations. This work clarifies the physics of this instability in the epoch before decoupling, and emphasizes that the universe is formally unstable on scales below the horizon, even in the limit of large taudot. Analogous instabilities at yet earlier epochs are also mentioned. (Abridged)

We present a Spitzer/IRS low-resolution mid-infrared spectrum (5-14 micron) of the luminous transient discovered in the nearby galaxy NGC 300 in May 2008. The spectrum, obtained three months after discovery, shows that the transient is very luminous in the mid-IR. Furthermore, the spectrum shows strong, broad emission features at 8 and 12 micron that are observed in Galactic carbon-rich proto-planetary nebulae. Combining these data with published optical and near-IR photometry obtained at the same epoch, we find that the mid-IR excess traced by the Spitzer spectrum accounts for ~20% of the total energy output. This component can be well explained by emission from ~3x10^-4 Msun of pre-existing progenitor dust at temperature ~400 K. The spectral energy distribution of the transient also shows a near-IR excess that can be explained by emission from newly-formed dust in the ejecta. Alternatively, both the near-IR and mid-IR excesses can together be explained by a single pre-existing geometrically thick dust shell. In light of the new observations obtained with Spitzer, we revisit the analysis of the optical spectra and kinematics, which were compared to the massive yellow-hypergiant IRC+10420 in previous studies. We show that proto-planetary nebulae share many properties with the NGC 300 transient and SN 2008S. We conclude that even though the explosion of a massive star (> 10 Msun) cannot be ruled out, an explosive event on a massive (6-10 Msun) carbon-rich AGB/super-AGB or post-AGB star is consistent with all observations of the transients and their progenitors presented thus far.

We compile observations of the surface mass density profiles of dense stellar systems, including globular clusters in the Milky Way and nearby galaxies, massive star clusters in nearby starbursts, nuclear star clusters in dwarf spheroidals and late-type disks, ultra-compact dwarfs, and galaxy spheroids spanning the range from low-mass cusp bulges and ellipticals to massive core ellipticals. We show that in all cases the maximum stellar surface density attained in the central regions of these systems is similar, Sigma_max ~ 10^11 M_sun/kpc^2 (~20 g/cm^2), despite the fact that the systems span 7 orders of magnitude in total stellar mass M_star, 5 in effective radius R_e, and have a wide range in effective surface density M_star/R_e^2. The surface density limit is reached on a wide variety of physical scales in different systems and is thus not a limit on three-dimensional stellar density. Given the very different formation mechanisms involved in these different classes of objects, we argue that a single piece of physics likely determines Sigma_max. The radiation fields and winds produced by massive stars can have a significant influence on the formation of both star clusters and galaxies, while neither supernovae nor black hole accretion are important in star cluster formation. We thus conclude that feedback from massive stars likely accounts for the observed Sigma_max, plausibly because star formation reaches an Eddington-like flux that regulates the growth of these diverse systems. This suggests that current models of galaxy formation, which focus on feedback from supernovae and active galactic nuclei, are missing a crucial ingredient.

Supermassive black holes (BHs) obey tight scaling relations between their mass and their host galaxy properties such as total stellar mass, velocity dispersion, and potential well depth. This has led to the development of self-regulated models for BH growth, in which feedback from the central BH halts its own growth upon reaching a critical threshold. However, models have also been proposed in which feedback plays no role: so long as a fixed fraction of the host gas supply is accreted, relations like those observed can be reproduced. Here, we argue that the scatter in the observed BH-host correlations, and its run with scale, presents a demanding constraint on any model for these correlations, and that it favors self-regulated models of BH growth. We show that the scatter in the stellar mass fraction within a radius R in observed ellipticals and spheroids increases strongly at small R. At fixed total stellar mass (or host velocity dispersion), on very small scales near the BH radius of influence, there is an order-of-magnitude scatter in the amount of gas that must have entered and formed stars. In short, the BH appears to 'know more' about the global host galaxy potential on large scales than the stars and gas supply on small scales. This is predicted in self-regulated models; however, models where there is no feedback would generically predict order-of-magnitude scatter in the BH-host correlations. Likewise, models in which the BH feedback in the 'bright' mode does not regulate the growth of the BH itself, but sets the stellar mass of the galaxy by inducing star formation or blowing out a mass in gas much larger than the galaxy stellar mass, are difficult to reconcile with the scatter on small scales.

Radiation pressure from the absorption and scattering of starlight by dust grains may be a crucial feedback mechanism in starburst galaxies and the self-gravitating parsec-scale disks that accompany the fueling of active galactic nuclei. I review the case for radiation pressure in both optically-thin and highly optically-thick contexts. I describe the conditions for which Eddington-limited star formation yields a characteristic flux of ~10^{13} L_sun/kpc^2, and I discuss the physical circumstances for which the flux from radiation pressure supported disks is below or above this value. In particular, I describe the young stellar disk on ~0.1 pc scales at the Galactic Center. I argue that its bolometric flux at formation, ~10^{15} L_sun/kpc^2, and the observed stellar mass and scale height imply that the disk may have been radiation pressure supported during formation.

Despite its probable occurrence in Nature, the accretion-induced collapse (AIC) of a white dwarf has not yet been directly observed (or, at least, it has not been unambiguously identified as such). In this contribution we summarize the observational signatures of AIC, emphasizing its possible role as both an optical and high-energy transient.

The progenitors of SN 2008S and the 2008 transient in NGC300 were dust-enshrouded, with extremely red mid-infrared (MIR) colors and relatively low luminosities. The transients were optically faint (-13 < M_V < -15) compared to normal core-collapse supernovae (SNe), and their spectra exhibited narrow emission lines. These events are unique among transient-progenitor pairs and hence constitute a new class. Whether they are true SNe or bright massive-star eruptions, we argue that their rate is ~20% of the SN rate. This fact is remarkable in light of the observation that a very small fraction of all massive stars have the MIR colors of the SN 2008S and NGC300 progenitors, as we show using MIR and optical luminosity, color, and variability properties of massive stars in M33. We find that the fraction of massive stars with colors consistent with these progenitors is 1/10000. In fact, only < 10 similar objects exist in M33 - all of which lie at the luminous red extremum of the AGB sequence. That these transients are relatively common with respect to SNe, while their progenitors are rare compared to the massive star population, implies that the dust-enshrouded phase is a short-lived phase in the lives of many massive stars. This shrouded epoch can occur only in the last ~10^4 yr before explosion, be it death or merely eruption. We discuss the implications of this finding for the evolution and census of ``low-mass'' massive stars (8-12 Msun), and we connect it with theoretical discussions of electron-capture SNe and the explosive birth of white dwarfs. A systematic census with (warm) Spitzer of galaxies in the local universe for analogous progenitors would significantly improve our knowledge of this channel to massive stellar explosions, and potentially to others with obscured progenitors. (Abridged)

Star formation is slow, in the sense that the gas consumption time is much longer than the dynamical time. It is also inefficient; essentially all star formation in local galaxies takes place in giant molecular clouds (GMCs), but the fraction of a GMC converted to stars is very small, ~5%. In the most luminous starbursts, the GMC lifetime is shorter than the main sequence lifetime of even the most massive stars, so that supernovae can play no role in GMC disruption. We investigate the disruption of GMCs across a wide range of galaxies, from normal spirals to the densest starbursts; we take into account the effects of HII gas pressure, shocked stellar winds, protostellar jets, and radiation pressure produced by the absorption and scattering of starlight on dust grains. In the Milky Way, we find that a combination of three mechanisms, jets, HII gas pressure, and radiation pressure, disrupts the clouds. In more rapidly star forming galaxies such as ``clump'' galaxies at high-redshift, ultra-luminous infrared galaxies (ULIRGs) and submillimeter galaxies, radiation pressure dominates natal cloud distribution. We predict the presence of 10-20 clusters with masses ~10^7 Msun in local ULIRGs such as Arp 220 and a similar number of clusters with M_* ~ 10^8 Msun in high redshift clump galaxies; submillimeter galaxies will have even more massive clusters. We find that the mass fraction of a GMC that ends up in stars is an increasing function of the gas surface density of a galaxy, reaching ~35% in the most luminous starbursts. Furthermore, the disruption of bubbles by radiation pressure stirs the interstellar medium to velocities of ~10 km/s in normal galaxies and to ~100 km/s in ULIRGs like Arp 220, consistent with observations. Thus, radiation pressure may play a dominant role in the ISM of star-forming galaxies.

Using observations from the literature, we show that the non-thermal radio luminosity (L) of supernova remnants (SNRs) is a strong function of the average gas surface density (Sigma) of the galaxy in which the remnants reside, from normal spirals to luminous starbursts. We combine a simple theory for electron cooling in SNRs with the observed radio luminosities to estimate the remnant magnetic field strength (B_SNR): the correlation between L and Sigma implies that B_SNR also increases with Sigma. We explore two interpretations of this correlation: (1) B_SNR is generated by post-shock magnetic field amplification, with B_SNR^2 proportional to Sigma and (2) B_SNR results from shock-compression of the ambient ISM magnetic field (B_ISM), with B_ISM being larger in denser galaxies. We find that shock compression is, on average, sufficient to produce the observed radio emission from SNRs in the densest starbursts; amplification of post-shock magnetic fields is not required. By contrast, in normal spirals post-shock field amplification (by a factor of a few - 10) is consistent with the data; we find tentative evidence that both the Alfven speed and the ratio of B_SNR^2 to the post-shock pressure ("epsilon_B") are constant in SNRs from galaxy to galaxy. We discuss observational tests that can be used to distinguish between these two interpretations of the radio luminosities of SNRs. Regardless of which is correct, the radio emission from SNRs provides an upper limit to B_ISM that is independent of the minimum energy assumption. For the densest starbursts, the ISM magnetic energy density is below that required for hydrostatic equilibrium; thus magnetic fields are not dynamically important on the largest scales in starbursts, in contrast with spiral galaxies like our own. This dichotomy may have implications for galactic dynamo theory.

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spindown of a newly formed millisecond, B ~ 10^{15} G, magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum, and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spindown powers (~ 10^{51}-10^{52} ergs/s), the magnetar wind is super-fast at almost all latitudes, while for lower spindown powers (~ 10^{50} erg/s), the wind is sub-fast but still super-Alfvenic. In all cases, the rates at which the neutron star loses mass, angular momentum, and energy are very similar to the corresponding free wind values (<~ 30% differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated (~5-10 deg.) relativistic jet out along the rotation axis of the star. Nearly all of the spindown power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.

Using two Chandra observations we have derived estimates of the dust distribution and distance to the eclipsing high mass X-ray binary (HMXB) Cen X-3 using the energy-resolved dust-scattered X-ray halo. By comparing the observed X-ray halos in 200 eV bands from 2-5 keV to the halo profiles predicted by the Weingartner & Draine interstellar grain model, we find that the vast majority (about 70%) of the dust along the line of sight to the system is located within about 300 pc of the Sun, although the halo measurements are insensitive to dust very close to the source. One of the Chandra observations occurred during an egress from eclipse as the pulsar emerged from behind the mass-donating primary. By comparing model halo light curves during this transition to the halo measurements, a source distance of 5.7 +/- 1.5 kpc (68% confidence level) is estimated, although we find this result depends on the distribution of dust on very small scales. Nevertheless, this value is marginally inconsistent with the commonly accepted distance to Cen X-3 of 8 kpc. We also find that the energy scaling of the scattering optical depth predicted by the Weingartner & Draine interstellar grain model does not accurately represent the results determined by X-ray halo studies of Cen X-3. Relative to the model, there appears to be less scattering at low energies or more scattering at high energies in Cen X-3.

I review aspects of the theory of long-duration gamma-ray burst (GRB) central engines. I focus on the requirements of any model; these include the angular momentum of the progenitor, the power, Lorentz factor, asymmetry, and duration of the flow, and both the association and the non-association with bright supernovae. I compare and contrast the collapsar and millisecond proto-magnetar models in light of these requirements. The ability of the latter model to produce a flow with Lorentz factor ~100 while simultaneously maintaining a kinetic luminosity of ~10^50 ergs/s for a timescale of ~10-100 s is emphasized.

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.

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.

I consider the physics of gravitational instabilities in the presence of dynamically important radiation pressure and gray radiative diffusion, governed by a constant opacity, kappa. For any non-zero radiation diffusion rate on an optically-thick scale, the medium is unstable unless the classical gas-only isothermal Jeans criterion is satisfied. When diffusion is "slow," although the dynamical Jeans instability is stabilized by radiation pressure on scales smaller than the adiabatic Jeans length, on these same spatial scales the medium is unstable to a diffusive mode. In this regime, neglecting gas pressure, the characteristic timescale for growth is independent of spatial scale and given by (3 kappa c_s^2)/(4 pi G c), where c_s is the adiabatic sound speed. This timescale is that required for a fluid parcel to radiate away its thermal energy content at the Eddington limit, the Kelvin-Helmholz timescale for a radiation pressure supported self-gravitating object. In the limit of "rapid" diffusion, radiation does nothing to suppress the Jeans instability and the medium is dynamically unstable unless the gas-only Jeans criterion is satisfied. I connect with treatments of Silk damping in the early universe. I discuss several applications, including photons diffusing in regions of extreme star formation (starburst galaxies & pc-scale AGN disks), and the diffusion of cosmic rays in normal galaxies and galaxy clusters. The former (particularly, starbursts) are "rapidly" diffusing and thus cannot be supported against dynamical instability in the linear regime by radiation pressure alone. The latter are more nearly "slowly" diffusing. I speculate that the turbulence in starbursts may be driven by the dynamical coupling between the radiation field and the self-gravitating gas, perhaps mediated by magnetic fields. (Abridged)

We present five epochs of simultaneous radio (VLA) and X-ray (Chandra) observations of SS 433, to study the relation between the radio and X-ray emission in the arcsecond-scale jets of the source. We detected X-ray emission from the extended jets in only one of the five epochs of observation, indicating that the X-ray reheating mechanism is transient. The reheating does not correlate with the total flux in the core or in the extended radio jets. However, the radio emission in the X-ray reheating regions is enhanced when X-ray emission is present. Deep images of the jets in linear polarization show that outside of the core, the magnetic field in the jets is aligned parallel to the local velocity vector, strengthening the case for the jets to be composed of discrete bullets rather than being continuous flux tubes. We also observed anomalous regions of polarized emission well away from the kinematic trace, confirming the large-scale anisotropy of the magnetic field in the ambient medium surrounding the jets.

The low-mass X-ray binary GS 1826-238 is presently unique for its consistently regular bursting behavior. In previous Rossi X-Ray Timing Explorer (RXTE) measurements between 1997 November and 2002 July, this source exhibited (nearly) limit-cycle bursts with recurrence times that decreased proportionately as the persistent flux increased. Here we report additional measurements of the burst recurrence time by RXTE, Chandra, and XMM-Newton, as well as observations of optical bursts. On a few occasions we measured burst recurrence times which deviated significantly from the earlier flux-recurrence time relationship, and most of these bursts occurred earlier than would be predicted based on the X-ray flux level. The epochs with early bursts were also accompanied by unusual broadband timing signatures, with the entire power spectrum shifting to higher frequencies. Concurrent XMM-Newton observations during one of these occasions, in 2003 April, indicate that an additional soft component may be present in the spectrum containing enough flux (30% of the total) to account for the burst recurrence time discrepancy. A self-consistent interpretation for the increase in soft flux and accompanying timing changes during 2003 April is that accretion disk extends down to smaller radial distances from the source than during the other observing epochs. The RXTE observations since 2003 April show that the spectral and timing properties have nearly returned to the previously established level.

Evidence is growing for a class of gamma-ray bursts (GRBs) characterized by an initial ~0.1-1 s spike of hard radiation followed, after a ~3-10 s lull in emission, by a softer period of extended emission lasting ~10-100 s. In a few well-studied cases, these ``short GRBs with extended emission'' show no evidence for a bright associated supernova (SN). We propose that these events are produced by the formation and early evolution of a highly magnetized, rapidly rotating neutron star (a ``proto-magnetar'') which is formed from the accretion-induced collapse (AIC) of a white dwarf (WD), the merger and collapse of a WD-WD binary, or, perhaps, the merger of a double neutron star binary. The initial emission spike is powered by accretion onto the proto-magnetar from a small disk that is formed during the AIC or merger event. The extended emission is produced by a relativistic wind that extracts the rotational energy of the proto-magnetar on a timescale ~10-100 s. The ~3-10 s delay between the prompt and extended emission is the time required for the newly-formed proto-magnetar to cool sufficiently that the neutrino-heated wind from its surface becomes ultra-relativistic. Because a proto-magnetar ejects little or no Ni56 (< 1e-3 M_sun), these events should not produce a bright SN-like transient. We model the extended emission from GRB060614 using spin-down calculations of a cooling proto-magnetar, finding reasonable agreement with observations for a magnetar with an initial rotation period of ~1 ms and a surface dipole field of ~3e15 G. If GRBs are indeed produced by AIC or WD-WD mergers, they should occur within a mixture of both early and late-type galaxies and should not produce strong gravitational wave emission. An additional consequence of our model is the existence of X-ray flashes unaccompanied by a bright SN.

We model the interaction between the wind from a newly formed rapidly rotating magnetar and the surrounding progenitor. In the first few seconds after core collapse the magnetar inflates a bubble of plasma and magnetic fields behind the supernova shock, which expands asymmetrically because of the pinching effect of the toroidal magnetic field, as in PWNe, even if the host star is spherically symmetric. The degree of asymmetry depends on the ratio of the magnetic energy to the total energy in the bubble. We assume that the wind by newly formed magnetars inflating these bubbles is more magnetized than for PWNe. We show that for a magnetic to total power supplied by the central magnetar $\sim 0.1$ the bubble expands relatively spherically while for values greater than 0.3, most of the pressure in the bubble is exerted close to the rotation axis, driving a collimated outflow out through the host star. This can account for the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). Given that the wind magnetization increases in time, we thus suggest that the magnetar-driven bubble initially expands relatively spherically (enhancing the energy of the associated supernova) while at late times it becomes progressivelymore collimated (producing the GRB). Similar processes may operate in more modestly rotating neutron stars to produce asymmetric supernovae and lower energy transients such as X-ray flashes.

We present time-dependent axisymmetric magnetohydrodynamic simulations of the interaction of a relativistic magnetized wind produced by a proto-magnetar with a surrounding stellar envelope, in the first $\sim 10$ seconds after core collapse. We inject a super-magnetosonic wind with $\dot E = 10^{51}$ ergs s$^{-1}$ into a cavity created by an outgoing supernova shock. A strong toroidal magnetic field builds up in the bubble of plasma and magnetic field that is at first inertially confined by the progenitor star. This drives a jet out along the polar axis of the star, even though the star and the magnetar wind are each spherically symmetric. The jet has the properties needed to produce a long-duration gamma-ray burst (GRB). At $\sim 5$ s after core bounce, the jet has escaped the host star and the Lorentz factor of the material in the jet at large radii $\sim 10^{11}$ cm is similar to that in the magnetar wind near the source. Most of the spindown power of the central magnetar escapes via the relativistic jet. There are fluctuations in the Lorentz factor and energy flux in the jet on $\sim 0.01-0.1$ second timescale. These may contribute to variability in GRB emission (e.g., via internal shocks).

We calculate the structure and neutron content of neutrino-heated MHD winds driven from the surface of newly-formed magnetars (``proto-magnetars'') and from the midplane of hyper-accreting disks, two of the possible central engines for gamma-ray bursts (GRBs) and hyper-energetic supernovae (SNe). Both the surface of proto-magnetars and the midplane of neutrino-cooled accretion flows (NDAFs) are electron degenerate and neutron-rich (neutron-to-proton ratio n/p >> 1). If this substantial free neutron excess is preserved to large radii in ultra-relativistic outflows, several important observational consequences may result. Weak interaction processes, however, can drive n/p to ~1 in the nondegenerate regions that obtain just above the surfaces of NDAFs and proto-magnetars. Our calculations show that mildly relativistic neutron-rich outflows from NDAFs are possible in the presence of a strong poloidal magnetic field. However, we find that neutron-rich winds possess a minimum mass-loss rate that likely precludes simultaneously neutron-rich and ultra-relativistic (Lorentz factor > 100) NDAF winds accompanying a substantial accretion power. In contrast, proto-magnetars are capable of producing neutron-rich long-duration GRB outflows ~10-30 seconds following core bounce for sub-millisecond rotation periods; such outflows would, however, accompany only extremely energetic events, in which the GRB + SN energy budget exceeds ~ 4e52 ergs. Neutron-rich highly relativistic outflows may also be produced during some short-duration GRBs by geometrically thick accretion disks formed from compact object mergers. The implications for r-process nucleosynthesis, optical transients due to non-relativistic neutron-rich winds, and Nickel production in proto-magnetar and NDAF winds are also briefly discussed.

We model the interaction between the wind from a newly formed rapidly rotating magnetar and the surrounding supernova shock and host star. The dynamics is modeled using the two-dimensional, axisymmetric thin-shell equations. In the first ~10-100 seconds after core collapse the magnetar inflates a bubble of plasma and magnetic fields behind the supernova shock. The bubble expands asymmetrically because of the pinching effect of the toroidal magnetic field, just as in the analogous problem of the evolution of pulsar wind nebulae. The degree of asymmetry depends on E_mag/E_tot. The correct value of E_mag/E_tot is uncertain because of uncertainties in the conversion of magnetic energy into kinetic energy at large radii in relativistic winds; we argue, however, that bubbles inflated by newly formed magnetars are likely to be significantly more magnetized than their pulsar counterparts. We show that for a ratio of magnetic to total power supplied by the central magnetar L_mag/L_tot ~ 0.1 the bubble expands relatively spherically. For L_mag/L_tot ~ 0.3, however, most of the pressure in the bubble is exerted close to the rotation axis, driving a collimated outflow out through the host star. This can account for the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). Outflows from magnetars become increasingly magnetically dominated at late times, due to the decrease in neutrino-driven mass loss as the young neutron star cools. We thus suggest that the magnetar-driven bubble initially expands relatively spherically, enhancing the energy of the associated supernova, while at late times it becomes progressively more collimated, producing the GRB.

We provide a model for how Kennicutt-Schmidt (KS) laws, which describe the correlation between star formation rate and gas surface or volume density, depend on the molecular line chosen to trace the gas. We show that, for lines that can be excited at low temperatures, the KS law depends on how the line critical density compares to the median density in a galaxy's star-forming molecular clouds. High critical density lines trace regions with similar physical properties across galaxy types, and this produces a linear correlation between line luminosity and star formation rate. Low critical density lines probe regions whose properties vary across galaxies, leading to a star formation rate that varies superlinearly with line luminosity. We show that a simple model in which molecular clouds are treated as isothermal and homogenous can quantitatively reproduce the observed correlations between galactic luminosities in far infrared and in the CO(1->0) and HCN(1->0) lines, and naturally explains why these correlations have different slopes. We predict that IR-line luminosity correlations should change slope for galaxies in which the median density is close to the line critical density. This prediction may be tested by observations of lines such as HCO^+(1->0) with intermediate critical densities, or by HCN(1->0) observations of intensely star-forming high redshift galaxies with very high densities. Recent observations by Gao et al. hint at just such a change in slope. We argue that deviations from linearity in the HCN(1->0)-IR correlation at high luminosity are consistent with the assumption of a constant star formation efficiency.

We begin by reviewing the theory of thermal, neutrino-driven proto-neutron star (PNS) winds. Including the effects of magnetic fields and rotation, we then derive the mass and energy loss from magnetically-driven PNS winds for both relativistic and non-relativistic outflows, including important multi-dimensional considerations. With these simple analytic scalings we argue that proto-magnetars born with ~ millisecond rotation periods produce relativistic winds just a few seconds after core collapse with luminosities, timescales, mass-loading, and internal shock efficiencies favorable for producing long-duration gamma-ray bursts.

Rapidly accreting massive protostars undergo a phase of deuterium shell burning during pre-main sequence evolution that causes them to swell to tenths of an AU in radius. During this phase, those with close binary companions will overflow their Roche lobes and begin transferring mass. Since massive stars frequently have companions at distances well under 1 AU, this process may affect the early evolution of a substantial fraction of massive stars. We use a simple protostellar evolution model to determine the range in accretion rates, mass ratios, and orbital separations for which mass transfer will occur, and we compute approximately the stability and final outcome of the transfer process. We discuss how mass transfer affects the demographics of massive binaries, and show that it provides a natural explanation for the heretofore unexplained population of massive "twins", high mass binaries with mass ratios very close to unity.

We solve the 1D neutrino-heated non-relativistic MHD wind problem for conditions that range from slowly rotating (spin period P > 10 ms) protoneutron stars (PNSs) with surface field strengths typical of radio pulsars (B < 10^13 G), to "proto-magnetars" with B ~ 10^14-10^15 G in their hypothesized rapidly rotating initial states (P ~ 1 ms). We use the simulations of Bucciantini et al. (2006) to map our monopole results onto a more physical dipole geometry and to estimate the spindown of PNSs when their winds are relativistic. We then quantify the effects of rotation and magnetic fields on the mass loss, energy loss, and r-process nucleosynthesis in PNS winds. We describe the evolution of PNS winds through the Kelvin-Helmholtz cooling epoch, emphasizing the transition between (1) thermal neutrino-driven, (2) non-relativistic magnetically-dominated, and (3) relativistic magnetically-dominated outflows. We find that proto-magnetars with P ~ 1 ms and B > 10^15 G drive relativistic winds with luminosities, energies, and Lorentz factors (magnetization sigma ~ 0.1-1000) consistent with those required to produce long duration gamma-ray bursts and hyper-energetic supernovae (SNe). A significant fraction of the rotational energy may be extracted in only a few seconds, sufficiently rapidly to alter the asymptotic energy of the SN remnant, its morphology, and, potentially, its nucleosynthetic yield. Winds from PNSs with more modest rotation periods (2 - 10 ms) and with magnetar-strength fields produce conditions significantly more favorable for the r-process than winds from slowly rotating PNSs. Lastly, we show that energy and momentum deposition by convectively-excited waves further increase the likelihood of successful r-process in PNS winds.

We solve the time-dependent dynamics of axisymmetric, general relativistic MHD winds from rotating neutron stars. The mass loss rate is obtained self consistently as a solution of the MHD equations, subject to a finite thermal pressure at the stellar surface. Conditions are chosen to be representative of the neutrino driven phase in newly born magnetars, which have been considered as a possible engine for GRBs. We compute the angular momentum and energy losses as a function of $\sigma$ and compare them with the analytic expectation from the classical theory of pulsar winds. We observe the convergence to the force-free limit in the energy loss and we study the evolution of the closed zone for increasing magnetization. Results also show that the dipolar magnetic field and the presence of a closed zone do not modify significantly the acceleration and collimation properties of the wind.

Magnetars are a sizable subclass of the neutron star census. Their very high magnetic field strengths are thought to be a consequence of rapid (millisecond) rotation at birth in a successful core-collapse supernova. In their first tens of seconds of existence, magnetars transition from hot, extended ``proto-''magnetars to the cooled and magnetically-dominated objects we identify $\sim10^4$ years later as Soft Gamma-ray Repeaters (SGRs) and Anamolous X-ray Pulsars (AXPs). Millisecond proto-magnetar winds during this cooling phase likewise transition from non-relativistic and thermally-driven to magneto-centrifugally-driven, and finally to relativistic and Poynting-flux dominated. Here we review the basic considerations associated with that transition. In particular, we discuss the spindown of millisecond proto-magnetars throughout the Kelvin-Helmholtz cooling epoch. Because of their large reservoir of rotational energy, their association with supernovae, and the fact that their winds are expected to become highly relativistic in the seconds after their birth, proto-magnetars have been suggested as the central engine of long-duration gamma ray bursts. We discuss some of the issues and outstanding questions in assessing them as such.

With recent and archival Rossi X-Ray Timing Explorer (RXTE) X-ray measurements of the heavily obscured X-ray pulsar EXO 1722-363 (IGR J17252-3616), we carried out a pulse timing analysis to determine the orbital solution for the first time. The binary system is characterized by a_x sin(i) = 101 +/- 3 lt-s and P_orb = 9.7403 +/- 0.0004 days (90% confidence), with the precision of the orbital period being obtained by connecting datasets separated by more than 7 years (272 orbital cycles). The orbit is consistent with circular, and e < 0.19 at the 90% confidence level. The mass function is 11.7 +/- 1.2 M_sun and confirms that this source is a High Mass X-ray Binary (HMXB) system. The orbital period, along with the previously known ~414 s pulse period, places this system in the part of the Corbet diagram populated by supergiant wind accretors. Using previous eclipse time measurements by Corbet et al. and our orbital solution, combined with the assumption that the primary underfills its Roche lobe, we find i > 61 degrees at the 99% confidence level, the radius of the primary is between 21 R_sun and 37 R_sun, and its mass is less than about 22 M_sun. The acceptable range of radius and mass shows that the primary is probably a supergiant of spectral type B0I-B5I. Photometric measurements of its likely counterpart are consistent with the spectral type and luminosity if the distance to the system is between 5.3 kpc and 8.7 kpc. Spectral analysis of the pulsar as a function of orbital phase reveals an evolution of the hydrogen column density suggestive of dense filaments of gas in the downstream wake of the pulsar, with higher levels of absorption seen at orbital phases 0.5-1.0, as well as a variable Fe K_alpha line.

Roughly 80% of Ultraluminous Infrared Galaxies (ULIRGs) show blue shifted absorption in the resonance lines of neutral sodium, indicating that cool winds are common in such objects, as shown by Rupke et al and by Martin. The neutral sodium (NaI) columns indicated by these absorption lines are ~ 10^{13}-3x10^{14}/cm^2, while the bolometric luminosity varies by a factor of only four. We show that the gas in ULIRG outflows is likely to be in photoionization equilibrium. The very small ULIRG sample of Goldader et al. demonstrates that the ratio of ultraviolet flux to far infrared flux varies by a factor $\sim100$ from object to object. While the Goldader sample does not overlap with those of Rupke et al. and Martin, we show that such a large variation in ultraviolet flux will produce a similar variation in the column of neutral sodium for a fixed mass flux and density. However, if the cold gas is in pressure equilibrium with a hot outflow with a mass loss rate similar to the star formation rate, the range of ionization state is significantly smaller. Measurements of the UV flux for objects in the Martin and Rupke et al. catalogs will definitively determine if photoionization effects are responsible for the wide variation seen in the sodium columns. If they are, a determination of the gas density and mass loss rate in the cool winds will follow, with attendant improvements in our understanding of wind driving mechanisms and of the effects of galaxies on their surroundings.

If cosmic ray protons interact with gas at roughly the mean density of the interstellar medium in starburst galaxies, then pion decay in starbursts is likely to contribute significantly to the diffuse extra-galactic background in both gamma-rays and high energy neutrinos. We describe the assumptions that lead to this conclusion and clarify the difference between our estimates and those of Stecker (2006). Detection of a single starburst by GLAST would confirm the significant contribution of starburst galaxies to the extra-galactic neutrino and gamma-ray backgrounds.

Cosmic ray protons interacting with gas at the mean density of the interstellar medium in starburst galaxies lose energy rapidly via inelastic collisions with ambient nuclei. The resulting pions produce secondary electrons and positrons, high-energy neutrinos, and gamma-ray photons. We estimate the cumulative gamma-ray emission from starburst galaxies. We find a total integrated background above 100 MeV of F_gamma ~ 10^{-6} GeV/cm^2/s/sr and a corresponding specific intensity at GeV energies of nuI_nu ~ 10^{-7} GeV/cm^2/s/sr. Starbursts may thus account for a significant fraction of the extra-galactic $\gamma$-ray background. We show that the FIR-radio correlation provides a strong constraint on the gamma-ray emission from starburst galaxies because pions decay into both gamma-rays and radio-emitting electron/positron pairs. We identify several nearby systems where the potential for observing gamma-ray emission is the most favorable (M82, NGC 253, and IC 342), predict their fluxes, and predict a linear FIR-gamma-ray correlation for the densest starbursts. If established, the FIR-gamma-ray correlation would provide strong evidence for the ``calorimeter'' theory of the FIR-radio correlation and would imply that cosmic rays in starburst galaxies interact with gas at approximately the mean density of the interstellar medium (ISM), thereby providing an important constraint on the physics of the ISM in starbursts.

I discuss neutrino production in supernovae (SNe) and the detection of both Galactic core collapse events and the diffuse extra-galactic MeV neutrino background expected from the integrated history of star formation. In particular, I consider what processes might affect our expectations for both. I focus on ``rapid'' rotation, defined as leading to millisecond initial neutron star spin periods. Rotation affects the neutrino luminosity, the average neutrino energy, the duration of the Kelvin-Helmholtz cooling epoch, and the ratios of luminosities and average energies between neutrino species; it can strongly suppresses the anti-electron as well as mu, anti-mu, tau, and anti-tau neutrino fluxes relative to those for the electron neutrinos. As a result, depending on the prevalence of rapid rotation in SN progenitors through cosmic time, this may affect predictions for the MeV neutrino background and the history of nucleosynthetic enrichment. I emphasize connections between the MeV neutrino background and tracers of the star formation rate density at high redshift in other neutrino and photon wavebands.

With the excellent angular resolution of the Chandra X-ray Observatory, it is possible to geometrically determine the distance to variable Galactic sources, based on the phenomenon that scattered radiation appearing in the X-ray halo has to travel along a slightly longer path than the direct, unscattered radiation. By measuring the delayed variability, constraints on the source distance can be obtained if the halo brightness is large enough to dominate the point spread function (PSF) and to provide sufficient statistics. The distance to Cyg X-3, which has a quasi-sinusoidal light curve, has been obtained with this approach by Predehl et al. Here we examine the feasibility of using the delayed signature of type I X-ray bursts as distance indicators. We use simulations of delayed X-ray burst light curves in the halo to find that the optimal annular region and energy band for a distance measurement with a grating observation is roughly 10-50" and 1-5 keV respectively, assuming Chandra's effective area and PSF, uniformly distributed dust, the input spectrum and optical depth to GX 13+1, and the Weingartner & Draine interstellar grain model. We find that the statistics are dominated by Poisson noise rather than systematic uncertainties, e.g., the PSF contribution to the halo. Using Chandra, a distance measurement to such a source at 4 (8) kpc could be made to about 23% (30%) accuracy with a single burst with 68% confidence. By stacking many bursts, a reasonable estimate of systematic errors limit the distance measurement to about 10% accuracy.

With recent and archival Rossi X-Ray Timing Explorer (RXTE) X-ray measurements of the heavily obscured X-ray pulsar IGR J16393-4643, we carried out a pulse timing analysis to determine the orbital parameters. Assuming a circular orbit, we phase-connected data spanning over 1.5 years. The most likely orbital solution has a projected semi-major axis of 43 +- 2 lt-s and an orbital period of 3.6875 +- 0.0006 days. This implies a mass function of 6.5 +- 1.1 M_sun and confirms that this INTEGRAL source is a High Mass X-ray Binary (HMXB) system. By including eccentricity in the orbital model, we find e < 0.25 at the 2 sigma level. The 3.7 day orbital period and the previously known ~910 s pulse period place the system in the region of the Corbet diagram populated by supergiant wind accretors, and the low eccentricity is also consistent with this type of system. Finally, it should be noted that although the 3.7 day solution is the most likely one, we cannot completely rule out two other solutions with orbital periods of 50.2 and 8.1 days.

We solve for the time-dependent dynamics of axisymmetric, general relativistic magnetohydrodynamic winds from rotating neutron stars. The mass loss rate is obtained self-consistently as a solution to the MHD equations, subject to a finite thermal pressure at the stellar surface. We consider both monopole and dipole magnetic field geometries and we explore the parameter regime extending from low magnetization (low-sigma_o), almost thermally-driven winds to high magnetization (high-sigma_o), relativistic Poynting-flux dominated outflows. We compute the angular momentum and rotational energy loss rates as a function of sigma_o and compare with analytic expectations from the classical theory of pulsars and magnetized stellar winds. In the case of the monopole, our high-sigma_o calculations asymptotically approach the analytic force-free limit. If we define the spindown rate in terms of the open magnetic flux, we similarly reproduce the spindown rate from recent force-free calculations of the aligned dipole. However, even for sigma_o as high as ~20, we find that the location of the Y-type point (r_Y), which specifies the radius of the last closed field line in the equatorial plane, is not the radius of the light cylinder R_L = c/omega (R = cylindrical radius), as has previously been assumed in most estimates and force-free calculations. Instead, although the Alfven radius at intermediate latitudes quickly approaches R_L as sigma_o exceeds unity, r_Y remains significantly less than R_L. Because r_Y < R_L, our calculated spindown rates thus exceed the classic ``vacuum dipole'' rate. We discussthe implications of our results for models of rotation-powered pulsars and magnetars, both in their observed states and in their hypothesized rapidly rotating initial state.

We estimate minimum energy magnetic fields (B_min) for a sample of galaxies with measured gas surface densities, spanning from normal spirals to starbursts. We show that the ratio of the minimum energy magnetic pressure to the total pressure in the ISM decreases substantially with increasing surface density; for Arp 220 this ratio is ~10^-4. Therefore, if the minimum energy estimate is applicable, magnetic fields in starbursts are dynamically weak compared to gravity, in contrast to normal spiral galaxies. We argue, however, that rapid cooling of relativistic electrons in starbursts invalidates the minimum energy estimate. We critically assess a number of independent constraints on the magnetic field strength in starbursts. In particular, we argue that the existence of the FIR-radio correlation implies that the synchrotron cooling timescale for cosmic ray electrons is much shorter than their escape time from the galactic disk; this in turn implies that the true magnetic field in starbursts is significantly larger than B_min. The strongest argument against such large fields is that one might expect starbursts to have steep radio spectra indicative of strong synchrotron cooling, which is not observed. We show, however, that ionization and bremsstrahlung losses can flatten the nonthermal spectra of starburst galaxies even in the presence of rapid cooling, providing much better agreement with observed spectra. We further demonstrate that ionization and bremsstrahlung losses are likely to be important in shaping the radio spectra of most starbursts at GHz frequencies, thereby preserving the linearity of the FIR-radio correlation. We thus conclude that magnetic fields in starbursts are significantly larger than B_min. We highlight several observations that can test this conclusion.

We present results from an extensive set of one- and two-dimensional radiation-hydrodynamic simulations of the supernova core collapse, bounce, and postbounce phases, and focus on the protoneutron star (PNS) spin periods and rotational profiles as a function of initial iron core angular velocity, degree of differential rotation, and progenitor mass. For the models considered, we find a roughly linear mapping between initial iron core rotation rate and PNS spin. The results indicate that the magnitude of the precollapse iron core angular velocities is the single most important factor in determining the PNS spin. Differences in progenitor mass and degree of differential rotation lead only to small variations in the PNS rotational period and profile. Based on our calculated PNS spins, at ~ 200-300 milliseconds after bounce, and assuming angular momentum conservation, we estimate final neutron star rotation periods. We find periods of one millisecond and shorter for initial central iron core periods of below ~ 10 s. This is appreciably shorter than what previous studies have predicted and is in disagreement with current observational data from pulsar astronomy. After considering possible spindown mechanisms that could lead to longer periods we conclude that there is no mechanism that can robustly spin down a neutron star from ~ 1 ms periods to the "injection" periods of tens to hundreds of milliseconds observed for young pulsars. Our results indicate that, given current knowledge of the limitations of neutron star spindown mechanisms, precollapse iron cores must rotate with periods around 50-100 seconds to form neutron stars with periods generically near those inferred for the radio pulsar population.

Using simultaneous observations from Chandra and RXTE, we investigated the LMXB GS 1826-238 with the goal of studying its spectral and timing properties. The uninterrupted Chandra observation captured 6 bursts (RXTE saw 3 of the 6), yielding a recurrence time of 3.54 +/- 0.03 hr. Using the proportional counter array on board RXTE, we made a probable detection of 611 Hz burst oscillations in the decaying phases of the bursts with an average rms signal amplitude of 4.8%. The integrated persistent emission spectrum can be described as the dual Comptonization of ~ 0.3 keV soft photons by a plasma with kT_e ~ 20 keV and an optical depth of about 2.6 (interpreted as emission from the accretion disk corona), plus the Comptonization of hotter ~ 0.8 keV seed photons by a ~ 6.8 keV plasma (interpreted as emission from or near the boundary layer). We discovered evidence for a neutral Fe K\alpha emission line, and we found interstellar Fe L_II and Fe L_III absorption features. The burst spectrum can be fit by fixing the disk Comptonization parameters to the persistent emission best-fit values, and adding a blackbody. The blackbody/seed photon temperature at the peak of the burst is ~ 1.8 keV and returns to ~ 0.8 keV over 200 s. The blackbody radius is consistent with R_bb = 10.3-11.7 km assuming a distance of 6 kpc; however, by accounting for the fraction of the surface that is obscured by the disk as a function of binary inclination, we determined the source distance must actually be near 5 kpc in order for the stellar radius to lie within the commonly assumed range of 10-12 km.

We consider the structure of marginally Toomre-stable starburst disks under the assumption that radiation pressure on dust grains provides the dominant vertical support against gravity. This is particularly appropriate when the disk is optically thick to its own IR radiation, as in the central regions of ULIRGs. Because the disk radiates at its Eddington limit, the Schmidt-law for star formation changes in the optically-thick limit, with the star formation rate per unit area scaling as Sigma_g/kappa, where Sigma_g is the gas surface density and kappa is the mean opacity. We show that optically thick starburst disks have a characteristic flux and dust effective temperature of F ~ 10^{13} L_sun/kpc^2 and T_eff ~ 90K, respectively. We compare our predictions with observations and find good agreement. We extend our model from many-hundred parsec scales to sub-parsec scales and address the problem of fueling AGN. We assume that angular momentum transport proceeds via global torques rather than a local viscosity. We account for the radial depletion of gas due to star formation and find a strong bifurcation between two classes of disk models: (1) solutions with a starburst on large scales that consumes all of the gas with little fueling of a central AGN and (2) models with an outer large-scale starburst accompanied by a more compact starburst on 1-10 pc scales and a bright central AGN. The luminosity of the latter models is in many cases dominated by the AGN. We show that the vertical thickness of the starburst disk on pc scales can approach h ~ r, perhaps accounting for the nuclear obscuration in some Type 2 AGN. We also argue that the disk of young stars in the Galactic Center may be the remnant of such a compact nuclear starburst.

In the seconds after core collapse and explosion, a thermal neutrino-driven wind emerges from the cooling, deleptonizing newly-born neutron star. If the neutron star has a large-scale magnetar-strength surface magnetic field and millisecond rotation period, then the wind is driven primarily by magneto-centrifugal slinging, and only secondarily by neutrino interactions. The strong magnetic field forces the wind to corotate with the stellar surface and the neutron star's rotational energy is efficiently extracted. As the neutron star cools, and the wind becomes increasingly magnetically-dominated, the outflow becomes relativistic. Here I review the millisecond magnetar model for long-duration gamma ray bursts and explore some of the basic physics of neutrino-magnetocentrifugal winds. I further speculate on some issues of collimation and geometry in the millisecond magnetar model.

We construct models of core-collapse supernovae in one spatial dimension, including rotation, angular momentum transport, and viscous dissipation employing an alpha-prescription. We compare the evolution of a fiducial 11 M_sun non-rotating progenitor with its evolution including a wide range of imposed initial rotation profiles (1.25<P_0<8 s, where P_0 is the initial, approximately solid-body, rotation period of the iron core). This range of P_0 covers the region of parameter space from where rotation begins to modify the dynamics (P_0~8 s) to where angular velocities at collapse approach Keplerian (P_0~1 s). Assuming strict angular momentum conservation, all models in this range leave behind neutron stars with spin periods <10 ms, shorter than those of most radio pulsars, but similar to those expected theoretically for magnetars at birth. A fraction of the gravitational binding energy of collapse is stored in the free energy of differential rotation. This energy source may be tapped by viscous processes, providing a mechanism for energy deposition that is not strongly coupled to the mass accretion rate through the stalled supernova shock. This effect yields qualitatively new dynamics in models of supernovae. We explore several potential mechanisms for viscosity in the core-collapse environment: neutrino viscosity, turbulent viscosity caused by the magnetorotational instability (MRI), and turbulent viscosity by entropy- and composition-gradient-driven convection. We argue that the MRI is the most effective. We find that for rotation periods in the range P_0<~5 s, and a range of viscous stresses, that the post-bounce dynamics is significantly effected by the inclusion of this extra energy deposition mechanism; in several cases we obtain strong supernova explosions.

We investigate large-scale galactic winds driven by momentum deposition. Momentum injection is provided by (1) radiation pressure produced by the continuum absorption and scattering of UV photons on dust grains and (2) supernovae. UV radiation can be produced by a starburst or AGN activity. We argue that momentum-driven winds are an efficient mechanism for feedback during the formation of galaxies. We show that above a limiting luminosity, momentum deposition from star formation can expel a significant fraction of the gas in a galaxy. The limiting, Eddington-like luminosity is $L_{\rm M}\simeq(4f_g c/G) \sigma^4$, where $\sigma$ is the galaxy velocity dispersion and $f_g$ is the gas fraction. A starburst that attains $L_{\rm M}$ moderates its star formation rate and its luminosity does not increase significantly further. We argue that ellipticals attain this limit during their growth at $z \gtrsim 1$ and that this is the origin of the Faber-Jackson relation. We show that Lyman break galaxies and ultra-luminous infrared galaxies have luminosities near $L_{\rm M}$. Star formation is unlikely to efficiently remove gas from very small scales in galactic nuclei, i.e., scales much smaller than that of a nuclear starburst. This gas is available to fuel a central black hole (BH). We argue that a BH clears gas out of its galactic nucleus when the luminosity of the BH itself reaches $\approx L_{\rm M}$. This shuts off the fuel supply to the BH and may also terminate star formation in the surrounding galaxy. As a result, the BH mass is fixed to be $M_{\rm BH}\simeq (f_g \kappa_{\rm es}/\pi G^2)\sigma^4$, where $\kappa_{\rm es}$ is the electron scattering opacity. This limit is in accord with the observed $M_{\rm BH}-\sigma$ relation. (Abridged)

Neutrino-matter cross sections and interaction rates are central to the core-collapse supernova phenomenon and, very likely, to the viability of the explosion mechanism itself. In this paper, we describe the major neutrino scattering, absorption, and production processes that together influence the outcome of core collapse and the cooling of protoneutron stars. One focus is on energy redistribution and many-body physics, but our major goal is to provide a useful resource for those interested in supernova neutrino microphysics.

The Kelvin-Helmholtz cooling epoch, lasting tens of seconds after the birth of a neutron star in a successful core-collapse supernova, is accompanied by a neutrino-driven wind. For magnetar-strength ($\sim10^{15}$ G) large scale surface magnetic fields, this outflow is magnetically-dominated during the entire cooling epoch.Because the strong magnetic field forces the wind to co-rotate with the protoneutron star,this outflow can significantly effect the neutron star's early angular momentum evolution, as in analogous models of stellar winds (e.g. Weber & Davis 1967). If the rotational energy is large in comparison with the supernova energy and the spindown timescale is short with respect to the time required for the supernova shockwave to traverse the stellar progenitor, the energy extracted may modify the supernova shock dynamics significantly. This effect is capable of producing hyper-energetic supernovae and, in some cases, provides conditions favorable for gamma ray bursts. We estimate spindown timescales for magnetized, rotating protoneutron stars and construct steady-state models of neutrino-magnetocentrifugally driven winds. We find that if magnetars are born rapidly rotating, with initial spin periods ($P$) of $\sim1$ millisecond, that of order $10^{51}-10^{52}$ erg of rotational energy can be extracted in $\sim10$ seconds. If magnetars are born slowly rotating ($P\gtrsim10$ ms) they can spin down to periods of $\sim1$ second on the Kelvin-Helmholtz timescale.

We consider nucleosynthesis in outflows originating from the inner regions of viscous accretion disks formed after the collapse of a rotating massive star. We show that wind-like outflows driven by viscous and neutrino heating can efficiently synthesize Fe-group elements moving at near-relativistic velocities. The mass of 56Ni synthesized and the asymptotic velocities attained in our calculations are in accord with those inferred from observations of SN1998bw and SN2003dh. These steady wind-like outflows are generally proton rich, characterized by only modest entropies, and consequently synthesize essentially nothing heavier than the Fe-group elements. We also discuss bubble-like outflows resulting from rapid energy deposition in localized regions near or in the accretion disk. These intermittent ejecta emerge with low electron fraction and are a promising site for the synthesis of the A=130 r-process peak elements.

We have developed a time-dependent, multi-energy-group, and multi-angle (S$_n$) Boltzmann transport scheme for radiation hydrodynamics simulations, in one and two spatial dimensions. The implicit transport is coupled to both 1D (spherically-symmetric) and 2D (axially-symmetric) versions of the explicit Newtonian hydrodynamics code VULCAN. The 2D variant, VULCAN/2D, can be operated in general structured or unstructured grids and though the code can address many problems in astrophysics it was constructed specifically to study the core-collapse supernova problem. Furthermore, VULCAN/2D can simulate the radiation/hydrodynamic evolution of differentially rotating bodies. We summarize the equations solved and methods incorporated into the algorithm and present results of a time-dependent 2D test calculation. A more complete description of the algorithm is postponed to another paper. We highlight a 2D test run that follows for 22 milliseconds the immediate post-bounce evolution of a collapsed core. We present the relationship between the anisotropies of the overturning matter field and the distribution of the corresponding flux vectors, as a function of energy group. This is the first 2D multi-group, multi-angle, time-dependent radiation/hydro calculation ever performed in core collapse studies. Though the transport module of the code is not gray and does not use flux limiters (however, there is a flux-limited variant of VULCAN/2D), it still does not include energy redistribution and most velocity-dependent terms.

Neutrino-driven winds are thought to accompany the Kelvin-Helmholtz cooling phase of nascent protoneutron stars in the first seconds after a core-collapse supernova. These outflows are a likely candidate as the astrophysical site for rapid neutron-capture nucleosynthesis (the r-process). In this chapter we review the physics of protoneutron star winds and assess their potential as a site for the production of the heavy r-process nuclides. We show that spherical transonic protoneutron star winds do not produce robust $r$-process nucleosynthesis for `canonical' neutron stars with gravitational masses of 1.4 M_sun and coordinate radii of 10 km. We further speculate on and review some aspects of neutrino-driven winds from protoneutron stars with strong magnetic fields.

(Abridged) We present results from dynamical models of core-collapse supernovae in one spatial dimension, employing a newly-developed Boltzmann neutrino radiation transport algorithm, coupled to Lagrangean hydrodynamics and a consistent high-density nuclear equation of state. We focus on shock breakout and its neutrino signature and follow the dynamical evolution of the cores of 11 M_sun, 15 M_sun, and 20 M_sun progenitors through collapse and the first 250 milliseconds after bounce. We examine the effects on the emergent neutrino spectra, light curves, and mix of species of artificial opacity changes, the number of energy groups, the weak magnetism/recoil corrections, nucleon-nucleon bremsstrahlung, neutrino-electron scattering, and the compressibility of nuclear matter. Furthermore, we present the first high-resolution look at the angular distribution of the neutrino radiation field both in the semi-transparent regime and at large radii and explore the accuracy with which our tangent-ray method tracks the free propagation of a pulse of radiation in a near vacuum. Finally, we fold the emergent neutrino spectra with the efficiencies and detection processes for a selection of modern underground neutrino observatories and argue that the prompt electron-neutrino breakout burst from the next galactic supernova is in principle observable and usefully diagnostic of fundamental collapse/supernova behavior. Though we are not in this study focusing on the supernova mechanism per se, our simulations support the theoretical conclusion (already reached by others) that spherical (1D) supernovae do not explode when good physics and transport methods are employed.

Because of their neutron-richness and association with supernovae, post-explosion protoneutron star winds are thought to be a likely astrophysical site for rapid neutron capture nucleosynthesis (the r-process). However, the most recent models of spherical neutrino-driven protoneutron star winds do not produce robust r-process nucleosynthesis for `canonical' neutron stars with a gravitational mass of 1.4M_sun and coordinate radius of 10km. These models fail variously; either the flow entropy is too low, the electron fraction is too high, or the dynamical expansion timescale is too long. To date, no models have included the effects of an ordered dipole magnetic field. We show that a strong magnetic field can trap the outflow in the neutrino heating region, thus leading to much higher matter entropy. We estimate both the trapping timescale and the resulting entropy amplification. For sufficiently large energy deposition rates, the trapped matter emerges dynamically from the region of closed magnetic field lines and escapes to infinity. We find that ordered dipoles with surface fields of greater than 6 times 10^14 Gauss increase the asymptotic entropy sufficiently for robust r-process nucleosynthesis.

Neutrino-matter interaction rates are central to the core collapse phenomenon and, perhaps, to the viability of the mechanism of core-collapse supernova explosions. In this paper we catalog and discuss the major neutrino scattering, absorption, and production processes that together influence the outcome of core collapse and the cooling of protoneutron stars. These are the essential inputs into the codes used to simulate the supernova phenomenon and an understanding of these processes is a prerequisite to continuing progress in supernova theory.

We review the status of the current quest to understand the mechanism of core-collapse supernovae, if neutrino-driven. In the process, we discuss the spherical explosion paradigm and its problems, some results from our new suite of collapse calculations performed using a recently-developed 1D implicit, multi-group, Feautrier/tangent-ray, Boltzmann solver coupled to explicit predictor/corrector hydrodynamics, the basic energetics of supernova explosions, and the promise of multi-D radiation/hydro simulations to explain why the cores of massive stars explode.

We solve the general-relativistic steady-state eigenvalue problem of neutrino-driven protoneutron star winds, which immediately follow core-collapse supernova explosions. We provide velocity, density, temperature, and composition profiles and explore the systematics and structures generic to such a wind for a variety of protoneutron star characteristics. Furthermore, we derive the entropy, dynamical timescale, and neutron-to-seed ratio in the general relativistic framework essential in assessing this site as a candidate for $r$-process nucleosynthesis. Generally, we find that for a given mass outflow rate ($\dot{M}$), the dynamical timescale of the wind is significantly shorter than previously thought. We argue against the existence or viability of a high entropy ($\gtrsim300$ per k$_{B}$ per baryon), long dynamical timescale $r$-process epoch. In support of this conclusion, we model the protoneutron star cooling phase, calculate nucleosynthetic yields in our steady-state profiles, and estimate the integrated mass loss. We find that transonic winds enter a high entropy phase only with very low $\dot{M}$ ($\lesssim1\times10^{-9}$ M$_\odot$ s$^{-1}$) and extremely long dynamical timescale ($\tau_\rho\gtrsim0.5$ seconds). Our results support the possible existence of an early $r$-process epoch at modest entropy ($\sim150$) and very short dynamical timescale, consistent in our calculations with a very massive or very compact protoneutron star that contracts rapidly after the preceding supernova. We explore possible modifications to our models, which might yield significant $r$-process nucleosynthesis generically. Finally, we speculate on the effect of fallback and shocks on both the wind physics and nucleosynthesis.

In preparation for a set of hydrodynamical simulations of core-collapse supernovae and protoneutron star winds, we investigate the rates of production and thermalization of $\nu_\mu$ and $\nu_\tau$ neutrinos in dense nuclear matter. Included are contributions from electron scattering, electron-positron annihilation, nucleon-nucleon bremsstrahlung, and nucleon scattering. We find that nucleon scattering dominates electron scattering as a thermalization process at neutrino energies greater than $\sim$15 MeV. In addition, nucleon-nucleon bremsstrahlung dominates electron-positron annihilation as a production mechanism at low neutrino energies, near and below the $\nu_\mu$ and $\nu_\tau$ neutrinospheres. Furthermore, we have begun a study of steady-state general relativistic protoneutron star winds employing simple neutrino heating and cooling terms. From this analysis we obtain acceleration profiles as well as asymptotic lepton fractions and baryon entropies essential in assessing the wind as a potential site for $r$-process nucleosynthesis.

We have developed an implicit, multi-group, time-dependent, spherical neutrino transport code based on the Feautrier variables, the tangent-ray method, and accelerated ${\bf \Lambda}$ iteration. The code achieves high angular resolution, is good to O($v/c$), is equivalent to a Boltzmann solver (without gravitational redshifts), and solves the transport equation at all optical depths with precision. In this paper, we present our formulation of the relevant numerics and microphysics and explore protoneutron star atmospheres for snapshot post-bounce models. Our major focus is on spectra, neutrino-matter heating rates, Eddington factors, angular distributions, and phase-space occupancies. In addition, we investigate the influence on neutrino spectra and heating of final-state electron blocking, stimulated absorption, velocity terms in the transport equation, neutrino-nucleon scattering asymmetry, and weak magnetism and recoil effects. Furthermore, we compare the emergent spectra and heating rates obtained using full transport with those obtained using representative flux-limited transport formulations to gauge their accuracy and viability. Finally, we derive useful formulae for the neutrino source strength due to nucleon-nucleon bremsstrahlung and determine bremsstrahlung's influence on the emergent $\nu_{\mu}$ and $\nu_{\tau}$ neutrino spectra.

We investigate the rates of production and thermalization of $\nu_\mu$ and $\nu_\tau$ neutrinos at temperatures and densities relevant to core-collapse supernovae and protoneutron stars. Included are contributions from electron scattering, electron-positron annihilation, nucleon-nucleon bremsstrahlung, and nucleon scattering. For the scattering processes, in order to incorporate the full scattering kinematics at arbitrary degeneracy, the structure function formalism developed by Reddy et al. (1998) and Burrows and Sawyer (1998) is employed. Furthermore, we derive formulae for the total and differential rates of nucleon-nucleon bremsstrahlung for arbitrary nucleon degeneracy in asymmetric matter. We find that electron scattering dominates nucleon scattering as a thermalization process at low neutrino energies ($\epsilon_\nu\lesssim 10$ MeV), but that nucleon scattering is always faster than or comparable to electron scattering above $\epsilon_\nu\simeq10$ MeV. In addition, for $\rho\gtrsim 10^{13}$ g cm$^{-3}$, $T\lesssim14$ MeV, and neutrino energies $\lesssim60$ MeV, nucleon-nucleon bremsstrahlung always dominates electron-positron annihilation as a production mechanism for $\nu_\mu$ and $\nu_\tau$ neutrinos.

We show that, in the context of Moore's Law, overall productivity can be increased for large enough computations by `slacking' or waiting for some period of time before purchasing a computer and beginning the calculation.

Previous works have shown that water Cherenkov detectors have superior sensitivity to those of scintillation counters as applied to detecting extensive air showers (EAS). This is in large part due to their much higher sensitivity to EAS photons which are more than five times more numerous than EAS electrons. Large area water Cherenkov detectors can be constructed relatively cheaply and operated reliably. A sparse detector array has been designed which uses these types of detectors to substantially increase the area over which the Milagro Gamma Ray Observatory collects EAS information. Improvements to the Milagro detector's performance characteristics and sensitivity derived from this array and preliminary results from a prototype array currently installed near the Milagro detector will be presented.

The Milagro detector is an air shower array which uses the water Cherenkov technique and is capable of continuously monitoring the sky at energies near 1 TeV. The detector consists of 20000 metric tons of pure water instrumented with 723 photo-multiplier tubes (PMTs). The PMTs are arranged in a two-layer structure on a lattice of 3 m spacing covering 5000 $m^2$ area. The direction of the shower is determined from the relative timing of the PMT signals, necessitating a common time reference and amplitude slewing corrections to improve the time resolution. The calibration system to provide these consists of a pulsed laser driving 30 diffusing light sources deployed in the pond to allow cross-calibration of the PMTs. The system is capable of calibrating times and the pulse-heights from the PMTs using the time-over-threshold technique. The absolute energy scale is provided using single muons passing through the detector. The description of the calibration system of the Milagro detector and its prototype Milagrito will be presented.

The Milagrito water Cherenkov telescope operated for over a year. The most probable gamma-ray energy was ~1 TeV and the trigger rate was as high as 400 Hz. We have developed an efficient technique for searching the entire sky for short duration bursts of TeV photons. Such bursts may result from "traditional" gamma-ray bursts that were not in the field-of-view of any other instruments, the evaporation of primordial black holes, or some as yet undiscovered phenomenon. We have begun to search the Milagrito data set for bursts of duration 10 seconds. Here we will present the technique and the expected results. Final results will be presented at the conference.

Solar Energetic Particles from the 6 November 1997 solar flare/CME(coronal mass ejection) with energies exceeding 10 GeV have been detected by Milagrito, a prototype of the Milagro Gamma Ray Observatory. While particle acceleration beyond 1 GeV at the Sun is well established, few data exist for protons or ions beyond 10 GeV. The Milagro observatory, a ground based water Cherenkov detector designed for observing very high energy gamma ray sources, can also be used to study the Sun. Milagrito, which operated for approximately one year in 1997/98, was sensitive to solar proton and neutron fluxes above ~5- 10 GeV. Milagrito operated in a scaler mode, which was primarily sensitive to muons, low energy photons, and electrons, and the detector operated in a mode sensitive to showers and high zenith angle muons. In its scaler mode, Milagrito registered a rate increase coincident with the 6 November 1997 ground level event observed by Climax and other neutron monitors. A preliminary analysis suggests the presence of >10 GeV particles.

Milagrito, a prototype for the Milagro detector, operated for 15 months in 1997-8 and collected 8.9 billion events. It was the first extensive air shower (EAS) array sensitive to showers intiated by primaries with energy below 1 TeV. The shadows of the sun and moon observed with cosmic rays can be used to study systematic pointing shifts and measure the angular resolution of EAS arrays. Below a few TeV, the paths of cosmic rays coming toward the earth are bent by the helio- and geo-magnetic fields. This is expected to distort and displace the shadows of the sun and the moon. The moon shadow, offset from the nominal (undeflected) position, has been observed with high statistical significance in Milagrito. This can be used to establish energy calibrations, as well as to search for the anti-matter content of the VHE cosmic ray flux. The shadow of the sun has also been observed with high significance.

The Milagro Gamma Ray Observatory is the world's first large-area water Cherenkov detector capable of continuously monitoring the sky at TeV energies. Located in northern New Mexico, Milagro will perform an all sky survey of the Northern Hemisphere at energies between ~250 GeV and 50 TeV. With a high duty cycle, large detector area (~5000 square meters), and a wide field-of-view (~1 sr), Milagro is uniquely capable of searching for transient and DC sources of high-energy gamma-ray emission. Milagro has been operating since February, 1999. The current status of the Milagro Observatory and initial results will be discussed.

The Milagro water Cherenkov detector near Los Alamos, New Mexico, has been operated as a sky monitor at energies of a few TeV between February 1997 and April 1998. Serving as a test run for the full Milagro detector, Milagrito has taken data during the strong and long-lasting 1997 flare of Mrk 501. We present results from the analysis of Mrk 501 and compare the excess and background rates with expectations from the detector simulations.

Observing gamma ray bursts (GRBs) in the TeV energy range can be extremely valuable in providing insight to GRB radiation mechanisms and in constraining source distances. The Milagrito detector was an air shower array which used the water Cherenkov technique to search for TeV sources. Data from this detector was analyzed to look for a TeV component of GRBs coincident with low energy -rays detected by the BATSE instrument on the Compton Gamma Ray Observatory. A sample of 54 BATSE GRBs which were in the field of view of the Milagrito detector during its lifetime (February 1997 to May 1998) was used.

The Milagrito water Cherenkov detector in the Jemez Mountains near Los Alamos, New Mexico took data from February 1997 to April 1998. Milagrito served as a prototype for the larger Milagro detector, which has just begun operations. Milagrito was the first large-aperture gamma-ray detector with sensitivity to gamma rays below 1 TeV. We report here on a search for steady emission from point sources over most of the northern sky using data from Milagrito.