Source author record

Todd A. Thompson

Todd A. Thompson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.GA astro-ph.SR astro-ph.CO astro-ph hep-ph gr-qc astro-ph.EP astro-ph.IM

Catalog footprint

What is connected

91works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2023arXiv

X-ray Properties of NGC 253's Starburst-Driven Outflow

We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central $\sim$250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and non-spherical outflow geometry. Our derived temperatures and densities yield few-Myr cooling times in the nuclear region, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates on the mass outflow rates which are $2.8\:M_{\odot}/\rm{yr}$ in the northern outflow and $3.2\:M_{\odot}/\rm{yr}$ in the southern outflow. Additionally, we detect emission from charge exchange and find it has a significant contribution ($20-42$\%) to the total broad-band ($0.5-7$ keV) X-ray emission in the central and southern regions of the outflow.

preprint2022arXiv

ALMA 200 pc imaging of a z~7 quasar reveals a compact, disk-like host galaxy

We present 0".035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x10^9 M_sun, about a factor two higher than that of the central supermassive black hole. Consequently we do not resolve the black hole's sphere of influence, and find no kinematic signature of the central supermassive black hole. Kinematic modeling of the [CII] line shows that the dynamical mass at large radii is consistent with the gas mass, leaving little room for a significant mass contribution by stars and/or dark matter. The Toomre-Q parameter is less than unity throughout the disk, and thus is conducive to star formation, consistent with the high infrared luminosity of the system. The dust in the central region is optically thick, at a temperature >132 K. Using standard scaling relations of dust heating by star formation, this implies an unprecedented high star formation rate density of >10^4 M_sun / yr / kpc^2. Such a high number can still be explained with the Eddington limit for star formation under certain assumptions, but could also imply that the central supermassive black hole contributes to the heating of the dust in the central 110 pc.

preprint2022arXiv

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

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

preprint2022arXiv

Galactic winds and bubbles from nuclear starburst rings

Galactic outflows from local starburst galaxies typically exhibit a layered geometry, with cool $10^4\,$K flow sheathing a hotter $10^7\,$K, cylindrically-collimated, X-ray emitting plasma. Here, we argue that winds driven by energy-injection in a ring-like geometry can produce this distinctive large-scale multi-phase morphology. The ring configuration is motivated by the observation that massive young star clusters are often distributed in a ring at the host galaxy's inner Lindblad resonance, where larger-scale spiral arm structure terminates. We present parameterized three-dimensional radiative hydrodynamical simulations that follow the emergence and dynamics of energy-driven hot winds from starburst rings. In this Letter, we show that the flow shocks on itself within the inner ring hole, maintaining high $10^7$\,K temperatures, whilst flows that emerge from the wind-driving ring unobstructed can undergo rapid bulk cooling down to $10^4\,$K, producing a fast hot bi-conical outflow enclosed by a sheath of cooler nearly co-moving material without ram-pressure acceleration. The hot flow is collimated along the ring axis, even in the absence of pressure confinement from a galactic disk or magnetic fields. In the early stages of expansion, the emerging wind forms a bubble-like shape reminiscent of the Milky Way's eROSITA and Fermi bubbles and can reach velocities usually associated with AGN-driven winds. We discuss the physics of the ring configuration, the conditions for radiative bulk cooling, and the implications for future X-ray observations.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

Cosmic rays across the star-forming galaxy sequence. I: Cosmic ray pressures and calorimetry

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.

preprint2021arXiv

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

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

preprint2021arXiv

Return of the TEDI: revisiting the Triple Evolution Dynamical Instability channel in triple stars

Triple-star systems exhibit a phenomenon known as the Triple Evolution Dynamical Instability (TEDI), in which mass loss in evolving triples triggers short-term dynamical instabilities, potentially leading to collisions of stars, exchanges, and ejections. Previous work has shown that the TEDI is an important pathway to head-on stellar collisions in the Galaxy, significantly exceeding the rate of collisions due to random encounters in globular clusters. Here, we revisit the TEDI evolutionary pathway using state-of-the-art population synthesis methods that self-consistently take into account stellar evolution and binary interactions, as well as gravitational dynamics and perturbations from passing stars in the field. We find Galactic TEDI-induced collision rates on the order of 1e-4/yr, consistent with previous studies which were based on more simplified methods. The majority of TEDI-induced collisions involve main sequence stars, potentially producing blue straggler stars. Collisions are also possible involving more evolved stars, potentially producing eccentric post-common-envelope systems, and white dwarfs collisions leading to Type Ia supernovae (although the latter with a negligible contribution to the Galactic rate). In our simulations, the TEDI is not only triggered by adiabatic wind mass loss, but also by Roche lobe overflow in the inner binary: when the donor star becomes less massive than the accretor, the inner binary orbit widens, triggering triple dynamical instability. We find that collision rates are increased by ~17% when fly-bys in the field are taken into account. In addition to collisions, we find that the TEDI produces ~1e-4/yr of unbound stars, although none with escape speeds in excess of 1e3 km/s.

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

Cosmic rays across the star-forming galaxy sequence. II: Stability limits and the onset of cosmic ray-driven outflows

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.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

The Antesonic Condition for the Explosion of Core-Collapse Supernovae II: Rotation and Turbulence

In the problem of steady free-fall onto a standing shockwave around acentral mass, the "antesonic" condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\leq3/16$, where $c_T$ is the isothermal sound speed in the subsonic post-shock flow, and $v_\mathrm{esc}$ is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and theshock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models ofaccretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a $\sim\!5$ per-cent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by $\sim\!20-40$ per-cent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multi-dimensional simulations explode more easily than their one-dimensional counterparts.

preprint2020arXiv

The Physical Nature of Starburst-Driven Galactic Outflows

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.

preprint2020arXiv

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

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

preprint2019arXiv

Cosmic Rays and Magnetic Fields in the Core and Halo of the Starburst M82: Implications for Galactic Wind Physics

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.

preprint2019arXiv

Nebular Spectra of 111 Type Ia Supernovae Disfavor Single Degenerate Progenitors

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

preprint2016arXiv

Constraining Gamma-Ray Emission from Luminous Infrared Galaxies with Fermi-LAT; Tentative Detection of Arp 220

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.

preprint2015arXiv

An Origin for Multi-Phase Gas in Galactic Winds and Halos

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.

preprint2015arXiv

Dynamical formation & scattering of hierarchical triples: Cross sections, Kozai-Lidov oscillations, and collisions

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.

preprint2015arXiv

Dynamics of Dusty Radiation Pressure Driven Shells and Clouds: Fast Outflows from Galaxies, Star Clusters, Massive Stars, and AGN

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.

preprint2015arXiv

Short GMC lifetimes: an observational estimate with the PdBI Arcsecond Whirlpool Survey (PAWS)

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.

preprint2015arXiv

Sub-Eddington Star-Forming Regions are Super-Eddington: Momentum Driven Outflows from Supersonic Turbulence

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.

preprint2015arXiv

The Gender Breakdown of the Applicant Pool for Tenure-Track Faculty Positions at a Sample of North American Research Astronomy Programs

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.

preprint2015arXiv

The Landscape of the Neutrino Mechanism of Core-Collapse Supernovae: Neutron Star and Black Hole Mass Functions, Explosion Energies and Nickel Yields

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.

preprint2014arXiv

High-Resolution Radio Continuum Measurements of the Nuclear Disks of Arp 220

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{Σ_{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{Σ_{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.

preprint2014arXiv

Hot Galactic Winds Constrained by the X-Ray Luminosities of Galaxies

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.

preprint2014arXiv

Neutrino-heated winds from rotating proto-magnetars

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.

preprint2014arXiv

Probing Cosmic-Ray Ion Acceleration with Radio-Submm and Gamma-Ray Emission from Interaction-Powered Supernovae

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.

preprint2014arXiv

Quasithermal Neutrinos from Rotating Protoneutron Stars Born during Core Collapse of Massive Stars

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.

preprint2014arXiv

Rapid Eccentricity Oscillations and the Mergers of Compact Objects in Hierarchical Triples

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.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Environmental Dependence of Giant Molecular Cloud Properties in M51

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.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Multi-phase cold gas kinematic of M51

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.

preprint2014arXiv

The Signature of Single-Degenerate Accretion Induced Collapse

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.

preprint2013arXiv

A Comparative Study of Giant Molecular Clouds in M51, M33 and the Large Magellanic Cloud

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.

preprint2013arXiv

Diffuse Hard X-ray Emission in Starburst Galaxies as Synchrotron from Very High Energy Electrons

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

preprint2013arXiv

Dust Scattering and the Radiation Pressure Force in the M82 Superwind

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.

preprint2013arXiv

Effects of stellar rotation on star formation rates and comparison to core-collapse supernova rates

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.

preprint2013arXiv

Gas Kinematics on GMC scales in M51 with PAWS: cloud stabilization through dynamical pressure

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

preprint2013arXiv

Greatly Enhanced Eccentricity Oscillations in Quadruple Systems Composed of Two Binaries: Implications for Stars, Planets, and Transients

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.

preprint2013arXiv

Modeling The Molecular Composition in an AGN Disk

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.

preprint2013arXiv

Numerical Simulations of Radiatively-Driven Dusty Winds

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

preprint2013arXiv

On WD-WD Mergers in Triple Systems: The Role of Kozai Resonance with Tidal Friction

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.

preprint2013arXiv

Probability Distribution Functions OF 12CO(J = 1-0) Brightness and Integrated Intensity in M51: The PAWS View

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)

preprint2013arXiv

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

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

preprint2013arXiv

The Impact of Metallicity on the Rate of Type Ia Supernovae

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

preprint2013arXiv

The Mass-Loss Induced Eccentric Kozai Mechanism: A New Channel for the Production of Close Compact Object-Stellar Binaries

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.

preprint2013arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy

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.

preprint2012arXiv

A Quadruple or Triple Origin For Tycho B and SN 1572

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.

preprint2012arXiv

Direct Numerical Simulation of Radiation Pressure-Driven Turbulence and Winds in Star Clusters and Galactic Disks

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

preprint2012arXiv

Effect of Collective Neutrino Oscillations on the Neutrino Mechanism of Core-Collapse Supernovae

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.

preprint2012arXiv

Measurement of the Mass and Stellar Population Distribution in M82 with the LBT

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 μ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γ, 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.

preprint2012arXiv

Radiation Pressure Driven Galactic Winds from Self-Gravitating Discs

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

preprint2012arXiv

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

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

preprint2012arXiv

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from MARVELS II: A Short-period Companion Orbiting an F Star with Evidence of a Stellar Tertiary And Significant Mutual Inclination

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 σ, 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.

preprint2011arXiv

Accelerating Compact Object Mergers in Triple Systems with the Kozai Resonance: A Mechanism for "Prompt" Type Ia Supernovae, Gamma-Ray Bursts, and Other Exotica

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.

preprint2011arXiv

New Class of High-Energy Transients from Crashes of Supernova Ejecta with Massive Circumstellar Material Shells

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.

preprint2011arXiv

On The GeV & TeV Detections of the Starburst Galaxies M82 & NGC 253

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.

preprint2011arXiv

Radiation Pressure Feedback in Galaxies

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.

preprint2011arXiv

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

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

preprint2011arXiv

The Physics of the Neutrino Mechanism of Core-Collapse Supernovae

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

preprint2011arXiv

Upper Bound on the Energy of Particles and Their Secondary Neutrinos

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.

preprint2010arXiv

Assessing Radiation Pressure as a Feedback Mechanism in Star-Forming Galaxies

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.

preprint2010arXiv

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

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

preprint2010arXiv

Proto-Magnetars as GRB Central Engines: Uncertainties, Limitations, and Particulars

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.

preprint2010arXiv

Radiation pressure from massive star clusters as a launching mechanism for super-galactic winds

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.

preprint2010arXiv

The Physics of the FIR-Radio Correlation: I. Calorimetry, Conspiracy, and Implications

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

preprint2010arXiv

The Physics of the FIR-Radio Correlation: II. Synchrotron Emission as a Star-Formation Tracer in High-Redshift Galaxies

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.

preprint2009arXiv

A Maximum Stellar Surface Density in Dense Stellar Systems

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.

preprint2009arXiv

Radio Emission from Supernova Remnants: Implications for Post-Shock Magnetic Field Amplification and the Magnetic Fields of 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.

preprint2009arXiv

Slow Diffusive Gravitational Instability Before Decoupling

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)

preprint2009arXiv

The Disruption of Giant Molecular Clouds by Radiation Pressure and the Efficiency of Star Formation in Galaxies

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.

preprint2007arXiv

Mass Transfer in Close, Rapidly Accreting Protobinaries: An Origin for Massive Twins?

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.

preprint2007arXiv

Proto-Neutron Star Winds with Magnetic Fields and Rotation

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.

preprint2006arXiv

Aspects of Neutrino Production in Supernovae

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.

preprint2006arXiv

The Ionization State of Sodium in Galactic Winds

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.

preprint2006arXiv

The Spin Periods and Rotational Profiles of Neutron Stars at Birth

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.

preprint2005arXiv

Millisecond Proto-Magnetars & Gamma Ray Bursts

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.

preprint2005arXiv

Radiation Pressure Supported Starburst Disks and AGN Fueling

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.

preprint2004arXiv

Magnetar Spindown, Hyper-Energetic Supernovae, and Gamma Ray Bursts

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.

preprint2003arXiv

Protoneutron Star Winds

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.

preprint2001arXiv

The Physics of Protoneutron Star Winds: Implications for r-Process Nucleosynthesis

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 ($τ_ρ\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.

Todd A. Thompson

What is connected

Connect this record

See the researcher in context

Building this map preview

91 published item(s)

X-ray Properties of NGC 253's Starburst-Driven Outflow

ALMA 200 pc imaging of a z~7 quasar reveals a compact, disk-like host galaxy

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

Galactic winds and bubbles from nuclear starburst rings

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

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

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

Cosmic rays across the star-forming galaxy sequence. I: Cosmic ray pressures and calorimetry

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

Return of the TEDI: revisiting the Triple Evolution Dynamical Instability channel in triple stars

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

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

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

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

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

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

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

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

Cosmic rays across the star-forming galaxy sequence. II: Stability limits and the onset of cosmic ray-driven outflows

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

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

The Antesonic Condition for the Explosion of Core-Collapse Supernovae II: Rotation and Turbulence

The Physical Nature of Starburst-Driven Galactic Outflows

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

Cosmic Rays and Magnetic Fields in the Core and Halo of the Starburst M82: Implications for Galactic Wind Physics

Nebular Spectra of 111 Type Ia Supernovae Disfavor Single Degenerate Progenitors

Constraining Gamma-Ray Emission from Luminous Infrared Galaxies with Fermi-LAT; Tentative Detection of Arp 220

An Origin for Multi-Phase Gas in Galactic Winds and Halos

Dynamical formation & scattering of hierarchical triples: Cross sections, Kozai-Lidov oscillations, and collisions

Dynamics of Dusty Radiation Pressure Driven Shells and Clouds: Fast Outflows from Galaxies, Star Clusters, Massive Stars, and AGN

Short GMC lifetimes: an observational estimate with the PdBI Arcsecond Whirlpool Survey (PAWS)

Sub-Eddington Star-Forming Regions are Super-Eddington: Momentum Driven Outflows from Supersonic Turbulence

The Gender Breakdown of the Applicant Pool for Tenure-Track Faculty Positions at a Sample of North American Research Astronomy Programs

The Landscape of the Neutrino Mechanism of Core-Collapse Supernovae: Neutron Star and Black Hole Mass Functions, Explosion Energies and Nickel Yields

High-Resolution Radio Continuum Measurements of the Nuclear Disks of Arp 220

Hot Galactic Winds Constrained by the X-Ray Luminosities of Galaxies

Neutrino-heated winds from rotating proto-magnetars

Probing Cosmic-Ray Ion Acceleration with Radio-Submm and Gamma-Ray Emission from Interaction-Powered Supernovae

Quasithermal Neutrinos from Rotating Protoneutron Stars Born during Core Collapse of Massive Stars

Rapid Eccentricity Oscillations and the Mergers of Compact Objects in Hierarchical Triples

The PdBI Arcsecond Whirlpool Survey (PAWS): Environmental Dependence of Giant Molecular Cloud Properties in M51

The PdBI Arcsecond Whirlpool Survey (PAWS): Multi-phase cold gas kinematic of M51

The Signature of Single-Degenerate Accretion Induced Collapse

A Comparative Study of Giant Molecular Clouds in M51, M33 and the Large Magellanic Cloud

Diffuse Hard X-ray Emission in Starburst Galaxies as Synchrotron from Very High Energy Electrons

Dust Scattering and the Radiation Pressure Force in the M82 Superwind

Effects of stellar rotation on star formation rates and comparison to core-collapse supernova rates

Gas Kinematics on GMC scales in M51 with PAWS: cloud stabilization through dynamical pressure

Greatly Enhanced Eccentricity Oscillations in Quadruple Systems Composed of Two Binaries: Implications for Stars, Planets, and Transients

Modeling The Molecular Composition in an AGN Disk

Numerical Simulations of Radiatively-Driven Dusty Winds

On WD-WD Mergers in Triple Systems: The Role of Kozai Resonance with Tidal Friction

Probability Distribution Functions OF 12CO(J = 1-0) Brightness and Integrated Intensity in M51: The PAWS View

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

The Impact of Metallicity on the Rate of Type Ia Supernovae

The Mass-Loss Induced Eccentric Kozai Mechanism: A New Channel for the Production of Close Compact Object-Stellar Binaries

The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy

A Quadruple or Triple Origin For Tycho B and SN 1572

Direct Numerical Simulation of Radiation Pressure-Driven Turbulence and Winds in Star Clusters and Galactic Disks

Effect of Collective Neutrino Oscillations on the Neutrino Mechanism of Core-Collapse Supernovae

Measurement of the Mass and Stellar Population Distribution in M82 with the LBT

Radiation Pressure Driven Galactic Winds from Self-Gravitating Discs

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

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from MARVELS II: A Short-period Companion Orbiting an F Star with Evidence of a Stellar Tertiary And Significant Mutual Inclination

Accelerating Compact Object Mergers in Triple Systems with the Kozai Resonance: A Mechanism for "Prompt" Type Ia Supernovae, Gamma-Ray Bursts, and Other Exotica

New Class of High-Energy Transients from Crashes of Supernova Ejecta with Massive Circumstellar Material Shells

On The GeV & TeV Detections of the Starburst Galaxies M82 & NGC 253

Radiation Pressure Feedback in Galaxies

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

The Physics of the Neutrino Mechanism of Core-Collapse Supernovae

Upper Bound on the Energy of Particles and Their Secondary Neutrinos

Assessing Radiation Pressure as a Feedback Mechanism in Star-Forming Galaxies

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

Proto-Magnetars as GRB Central Engines: Uncertainties, Limitations, and Particulars