BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daniel Kasen

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.COastro-ph.SRastro-ph.GAastro-phgr-qcnucl-thastro-ph.IMhep-phphysics.flu-dyn

Catalog footprint

What is connected

60works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

60 published item(s)

preprint2022arXiv

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they undergo core-collapse supernovae (CCSNe). The resulting dense circumstellar medium (CSM) may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomena were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2-4 km/s convective motions seed the intrinsic convective instability associated with the high entropy material deep in the envelope, enabling gas from deep within the envelope to escape, increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by 10x along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Carbon/Oxygen-dominated Atmosphere Days After Explosion for the "Super-Chandrasekhar" Type Ia SN 2020esm

Seeing pristine material from the donor star in a Type Ia supernova (SN Ia) explosion can reveal the nature of the binary system. In this paper, we present photometric and spectroscopic observations of SN 2020esm, one of the best-studied SNe of the class of "super-Chandrasekhar" SNe Ia (SC SNe Ia), with data obtained $-12$ to +360 days relative to peak brightness, obtained from a variety of ground- and space-based telescopes. Initially misclassified as a Type II supernova, SN 2020esm peaked at $M_{B} = -19.9$ mag, declined slowly ($Δm_{15}(B) = 0.92$ mag), and had particularly blue UV and optical colors at early times. Photometrically and spectroscopically, SN 2020esm evolved similarly to other SC SNe Ia, showing the usual low ejecta velocities, weak intermediate mass elements (IMEs), and the enhanced fading at late times, but its early spectra are unique. Our first few spectra (corresponding to a phase of $\gtrsim$10~days before peak) reveal a nearly-pure carbon/oxygen atmosphere during the first days after explosion. This composition can only be produced by pristine material, relatively unaffected by nuclear burning. The lack of H and He may further indicate that SN 2020esm is the outcome of the merger of two carbon/oxygen white dwarfs (WDs). Modeling its bolometric light curve, we find a $^{56}$Ni mass of $1.23^{+0.14}_{-0.14}$ M$_{\odot}$ and an ejecta mass of $1.75^{+0.32}_{-0.20}$ M$_{\odot}$, in excess of the Chandrasekhar mass. Finally, we discuss possible progenitor systems and explosion mechanisms of SN 2020esm and, in general, the SC SNe Ia class.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nebular Models of Sub-Chandrasekhar Mass Type Ia Supernovae: Clues to the Origin of Ca-rich Transients

We use non-local thermal equilibrium radiative transport modeling to examine observational signatures of sub-Chandrasekhar mass double detonation explosions in the nebular phase. Results range from spectra that look like typical and subluminous Type Ia supernovae (SNe) for higher mass progenitors to spectra that look like Ca-rich transients for lower mass progenitors. This ignition mechanism produces an inherent relationship between emission features and the progenitor mass as the ratio of the nebular [Ca II]/[Fe III] emission lines increases with decreasing white dwarf mass. Examining the [Ca II]/[Fe III] nebular line ratio in a sample of observed SNe we find further evidence for the two distinct classes of SNe Ia identified in Polin et al. by their relationship between Si II velocity and B-band magnitude, both at time of peak brightness. This suggests that SNe Ia arise from more than one progenitor channel, and provides an empirical method for classifying events based on their physical origin. Furthermore, we provide insight to the mysterious origin of Ca-rich transients. Low-mass double detonation models with only a small mass fraction of Ca (1%) produce nebular spectra that cool primarily through forbidden [Ca II] emission.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparing Moment-Based and Monte Carlo Methods of Radiation Transport Modeling for Type II-Plateau Supernova Light Curves

Time-dependent electromagnetic signatures from core-collapse supernovae are the result of detailed transport of the shock-deposited and radioactively-powered radiation through the stellar ejecta. Due to the complexity of the underlying radiative processes, considerable approximations are made to simplify key aspects of the radiation transport problem. We present a systematic comparison of the moment-based radiation hydrodynamical code STELLA and the Monte Carlo radiation transport code Sedona in the 1D modeling of Type II-Plateau supernovae. Based on explosion models generated from the Modules for Experiments in Stellar Astrophysics (MESA) instrument, we find remarkable agreements in the modeled light curves and the ejecta structure thermal evolution, affirming the fidelity of both radiation transport modeling approaches. The radiative moments computed directly by the Monte Carlo scheme in Sedona also verify the accuracy of the moment-based scheme. We find that the coarse resolutions of the opacity tables and the numerical approximations in STELLA have insignificant impact on the resulting bolometric light curves, making it an efficient tool for the specific task of optical light curve modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Observable signatures of dark photons from supernovae

A dark photon is a well-motivated new particle which, as a component of an associated dark sector, could explain dark matter. One strong limit on dark photons arises from excessive cooling of supernovae. We point out that even at couplings where too few dark photons are produced in supernovae to violate the cooling bound, they can be observed directly through their decays. Supernovae produce dark photons which decay to positrons, giving a signal in the 511 keV annihilation line observed by SPI/INTEGRAL. Further, prompt gamma-ray emission by these decaying dark photons gives a signal for gamma-ray telescopes. Existing GRS observations of SN1987a already constrain this, and a future nearby SN could provide a detection. Finally, dark photon decays from extragalactic SN would produce a diffuse flux of gamma rays observable by detectors such as SMM and HEAO-1. Together these observations can probe dark photon couplings several orders of magnitude beyond current constraints for masses of roughly 1 - 100 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Observing the earliest moments of supernovae using strong gravitational lenses

We determine the viability of exploiting lensing time delays to observe strongly gravitationally lensed supernovae (gLSNe) from first light. Assuming a plausible discovery strategy, the Legacy Survey of Space and Time (LSST) and the Zwicky Transient Facility (ZTF) will discover $\sim$ 110 and $\sim$ 1 systems per year before the supernova (SN) explosion in the final image respectively. Systems will be identified $11.7^{+29.8}_{-9.3}$ days before the final explosion. We then explore the possibility of performing early-time observations for Type IIP and Type Ia SNe in LSST-discovered systems. Using a simulated Type IIP explosion, we predict that the shock breakout in one trailing image per year will peak at $\lesssim$ 24.1 mag ($\lesssim$ 23.3) in the $B$-band ($F218W$), however evolving over a timescale of $\sim$ 30 minutes. Using an analytic model of Type Ia companion interaction, we find that in the $B$-band we should observe at least one shock cooling emission event per year that peaks at $\lesssim$ 26.3 mag ($\lesssim$ 29.6) assuming all Type Ia gLSNe have a 1 M$_\odot$ red giant (main sequence) companion. We perform Bayesian analysis to investigate how well deep observations with 1 hour exposures on the European Extremely Large Telescope would discriminate between Type Ia progenitor populations. We find that if all Type Ia SNe evolved from the double-degenerate channel, then observations of the lack of early blue flux in 10 (50) trailing images would rule out more than 27% (19%) of the population having 1 M$_\odot$ main sequence companions at 95% confidence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Radiative Emission Mechanisms of Tidal Disruption Events

We describe how the various outcomes of stellar tidal disruption give rise to observable radiation. We separately consider the cases where gas circularizes rapidly into an accretion disc, as well as the case when shocked debris streams provide the observable emission without having fully circularized. For the rapid circularization case, we describe how outflows, absorption by reprocessing layers, and Comptonization can cause the observed radiation to depart from that of a bare disc, possibly giving rise to the observed optical/UV emission along with soft X-rays from the disc. If, instead, most of the debris follows highly eccentric orbits for a significant time, many properties of the observed optical/UV emission can be explained by the scale of those eccentric orbits and the shocks embedded in the debris flow near orbital apocenter. In this picture, soft X-ray emission at early times results from the smaller amount of debris mass deflected into a compact accretion disc by weak shocks near the stellar pericenter. A general proposal for the near-constancy of the ultraviolet/optical color temperatures is provided, by linking it to incomplete thermalization of radiation in the atmosphere of the emitting region. We also briefly discuss the radio signals from the interaction of unbound debris and jets with the black hole environment.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SN 2019ehk: A Double-Peaked Ca-rich Transient with Luminous X-ray Emission and Shock-Ionized Spectral Features

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$α$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Effect of Jet-Ejecta Interaction on the Viewing Angle Dependence of Kilonova Light Curves

The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet-ejecta interaction as initial conditions for multi-dimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at 1.5h after merger. We find that on this timescale, jet shock heating does not affect the kilonova emission. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near ($\lesssim 30 °$) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Ultraviolet Spectrum of the Tidal Disruption Flare ASASSN-14li

We present a Hubble Space Telescope STIS spectrum of ASASSN-14li, the first rest-frame UV spectrum of a tidal disruption flare (TDF). The underlying continuum is well fit by a blackbody with $T_{\mathrm{UV}} = 3.5 \times 10^{4}$ K, an order of magnitude smaller than the temperature inferred from X-ray spectra (and significantly more precise than previous efforts based on optical and near-UV photometry). Super-imposed on this blue continuum, we detect three classes of features: narrow absorption from the Milky Way (probably a high-velocity cloud), and narrow absorption and broad (FWHM $\approx 2000$-8000 km s$^{-1}$) emission lines at/near the systemic host velocity. The absorption lines are blueshifted with respect to the emission lines by $Δv = -(250$-400) km s$^{-1}$. Due both to this velocity offset and the lack of common low-ionization features (Mg II, Fe II), we argue these arise from the same absorbing material responsible for the low-velocity outflow discovered at X-ray wavelengths. The broad nuclear emission lines display a remarkable abundance pattern: N III], N IV], He II are quite prominent, while the common quasar emission lines of C III] and Mg II are weak or entirely absent. Detailed modeling of this spectrum will help elucidate fundamental questions regarding the nature of the emission process(es) at work in TDFs, while future UV spectroscopy of ASASSN-14li would help to confirm (or refute) the previously proposed connection between TDFs and "N-rich" quasars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Binary Stars Can Provide the "Missing Photons" Needed for Reionization

Empirical constraints on reionization require galactic ionizing photon escape fractions fesc>20%, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ~1-5%. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the BPASS model). Binary mass transfer and mergers enhance the population of massive stars at late times (>3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the ISM, and so efficiently leak out of galaxies. As a result, the time-averaged "effective" escape fraction (ratio of escaped ionizing photons to observed 1500 A photons) increases by factors 4-10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Late-time Spectroscopy of Type Iax Supernovae

We examine the late-time (t > 200 days after peak brightness) spectra of Type Iax supernovae (SNe Iax), a low-luminosity, low-energy class of thermonuclear stellar explosions observationally similar to, but distinct from, Type Ia supernovae. We present new spectra of SN 2014dt, resulting in the most complete published late-time spectral sequence of a SN Iax. At late times, SNe Iax have generally similar spectra, all with a similar continuum shape and strong forbidden-line emission. However, there is also significant diversity where some late-time SN Iax spectra display narrow P-Cygni features and a continuum indicative of a photosphere in addition to strong narrow forbidden lines, while others have no obvious P-Cygni features, strong broad forbidden lines, and weak narrow forbidden lines. Finally, some SNe Iax have spectra intermediate to these two varieties with weak P-Cygni features and broad/narrow forbidden lines of similar strength. We find that SNe Iax with strong broad forbidden lines also tend to be more luminous and have higher-velocity ejecta at peak brightness. We estimate blackbody and kinematic radii of the late-time photosphere, finding the latter an order of magnitude larger than the former. We propose a two-component model that solves this discrepancy and explains the diversity of the late-time spectra of SNe Iax. In this model, the broad forbidden lines originate from the SN ejecta, while the photosphere, P-Cygni lines, and narrow forbidden lines originate from a wind launched from the remnant of the progenitor white dwarf and is driven by the radioactive decay of newly synthesized material left in the remnant. The relative strength of the two components accounts for the diversity of late-time SN Iax spectra. This model also solves the puzzle of a long-lived photosphere and slow late-time decline of SNe Iax. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Synchrotron Magnetic Fields from Rayleigh-Taylor Instability in Supernovae

Synchrotron emission from a supernova necessitates a magnetic field, but it is unknown how strong the relevant magnetic fields are, and what mechanism generates them. In this study, we perform high-resolution numerical gas dynamics calculations to determine the growth of turbulence due to Rayleigh-Taylor instability, and the resulting kinetic energy in turbulent fluctuations, to infer the strength of magnetic fields amplified by this turbulence. We find that Rayleigh-Taylor instability can produce turbulent fluctuations strong enough to amplify magnetic fields to a few percent of equipartition with the thermal energy. This turbulence stays concentrated near the reverse shock, but averaging this magnetic energy throughout the shocked region (weighting by emissivity) sets the magnetic fields at a minimum of 0.3 percent of equipartition. This suggests a minimum effective magnetic field strength ($ε_B > 0.003$) which should be present in all interacting supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Evolution and Fate of Super-Chandrasekhar Mass White Dwarf Merger Remnants

We present stellar evolution calculations of the remnant of the merger of two carbon-oxygen white dwarfs (CO WDs). We focus on cases that have a total mass in excess of the Chandrasekhar mass. After the merger, the remnant manifests as an $L \sim 3 \times 10^4 L_\odot$ source for $\sim 10^4$ yr. A dusty wind may develop, leading these sources to be self-obscured and to appear similar to extreme AGB stars. Roughly $\sim 10$ such objects should exist in the Milky Way and M31 at any time. As found in previous work, off-center carbon fusion is ignited within the merger remnant and propagates inward via a carbon flame, converting the WD to an oxygen-neon (ONe) composition. By following the evolution for longer than previous calculations, we demonstrate that after carbon-burning reaches the center, neutrino-cooled Kelvin-Helmholtz contraction leads to off-center neon ignition in remnants with masses $\ge 1.35 M_\odot$. The resulting neon-oxygen flame converts the core to a silicon WD. Thus, super-Chandrasekhar WD merger remnants do not undergo electron-capture induced collapse as traditionally assumed. Instead, if the remnant mass remains above the Chandrasekhar mass, we expect that it will form a low-mass iron core and collapse to form a neutron star. Remnants that lose sufficient mass will end up as massive, isolated ONe or Si WDs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The X-ray through Optical Fluxes and Line Strengths of Tidal Disruption Events

Observations of luminous flares resulting from the possible tidal disruption of stars by supermassive black holes have raised a number of puzzles. Outstanding questions include the origin of the optical and ultraviolet (UV) flux, the weakness of hydrogen lines in the spectrum, and the occasional simultaneous observation of x-rays. Here we study the emission from tidal disruption events (TDEs) produced as radiation from black hole accretion propagates through an extended, optically thick envelope formed from stellar debris. We analytically describe key physics controlling spectrum formation, and present detailed radiative transfer calculations that model the spectral energy distribution (SED) and optical line strengths of TDEs near peak brightness. The steady-state transfer is coupled to a non local thermodynamic equilibrium treatment of the excitation and ionization states of hydrogen, helium and oxygen (as a representative metal). Our calculations show how an extended envelope can reprocess a fraction of soft x-rays and produce the observed optical fluxes of order 10^43 ergs per second. Variations in the mass or size of the envelope may help explain how the optical flux changes over time with roughly constant color. For high enough accretion luminosities, x-rays can highly ionize the reprocessing region and escape to be observed simultaneously with the optical flux, producing an SED not described by a single blackbody. Due to optical depth effects, hydrogen Balmer line emission is often strongly suppressed relative to helium line emission (with HeII-to-H line ratios of at least 5:1 in some cases) even in the disruption of a solar-composition star. We discuss the implications of our results to understanding the type of stars destroyed in TDEs and the physical processes responsible for producing the observed flares.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Atomic Chemistry In Turbulent Astrophysical Media I: Effect of Atomic Cooling

We carry out direct numerical simulations of turbulent astrophysical media that explicitly track ionizations, recombinations, and species-by-species radiative cooling. The simulations assume solar composition and follows the evolution of hydrogen, helium, carbon, oxygen, sodium, and magnesium, but they do not include the presence of an ionizing background. In this case, the medium reaches a global steady state that is purely a function of the one-dimensional turbulent velocity dispersion, $σ_{\rm 1D},$ and the product of the mean density and the driving scale of turbulence, $n L.$ Our simulations span a grid of models with $σ_{\rm 1D}$ ranging from 6 to 58 km s$^{-1}$ and $n L$ ranging from 10$^{16}$ to 10$^{20}$ cm$^{-2},$ which correspond to turbulent Mach numbers from $M=0.2$ to 10.6. The species abundances are well described by single-temperature estimates whenever $M$ is small, but local equilibrium models can not accurately predict the global equilibrium abundances when $M \gtrsim 1.$ To allow future studies to account for nonequilibrium effects in turbulent media, we gather our results into a series of tables, which we will extend in the future to encompass a wider range of elements, compositions, and ionizing processes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetar Driven Shock Breakout and Double Peaked Supernova Light Curves

The light curves of some luminous supernovae are suspected to be powered by the spindown energy of a rapidly rotating magnetar. Here we describe a possible signature of the central engine: a burst of shock breakout emission occurring several days after the supernova explosion. The energy input from the magnetar inflates a high-pressure bubble that drives a shock through the pre-exploded supernova ejecta. If the magnetar is powerful enough, that shock will near the ejecta surface and become radiative. At the time of shock breakout, the ejecta will have expanded to a large radius (~10^{14} cm) so that the radiation released is at optical/ultraviolet wavelengths (T ~ 20,000 K) and lasts for several days. The luminosity and timescale of this magnetar driven shock breakout are similar to the first peak observed recently in the double-peaked light curve of SNLSQ14BDQ. However, for a large region of model parameter space, the breakout emission is predicted to be dimmer than the diffusive luminosity from direct magnetar heating. A distinct double peaked light curve may therefore only be conspicuous if thermal heating from the magnetar is suppressed at early times. We describe how such a delay in heating may naturally result from inefficient dissipation and thermalization of the pulsar wind magnetic energy. Without such suppression, the breakout may only be noticeable as a small bump or kink in the early luminosity or color evolution, or as a small but abrupt rise in the photospheric velocity. A similar breakout signature may accompany other central engines in supernovae, such as a black hole accreting fallback material.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Monte Carlo Radiation Hydrodynamics with Implicit Methods

We explore the application of Monte Carlo transport methods to solving coupled radiation-hydrodynamics problems. We use a time-dependent, frequency-dependent, 3-dimensional radiation transport code, that is special relativistic and includes some detailed microphysical interactions such as resonant line scattering. We couple the transport code to two different 1-dimensional (non-relativistic) hydrodynamics solvers: a spherical Lagrangian scheme and a Eulerian Godunov solver. The gas-radiation energy coupling is treated implicitly, allowing us to take hydrodyanimcal time-steps that are much longer than the radiative cooling time. We validate the code and assess its performance using a suite of radiation hydrodynamical test problems, including ones in the radiation energy dominated regime. We also develop techniques that reduce the noise of the Monte Carlo estimated radiation force by using the spatial divergence of the radiation pressure tensor. The results suggest that Monte Carlo techniques hold promise for simulating the multi-dimensional radiation hydrodynamics of astrophysical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Super-Eddington Stellar Winds Driven by Near-Surface Energy Deposition

We develop analytic and numerical models of the properties of super-Eddington stellar winds, motivated by phases in stellar evolution when super-Eddington energy deposition (via, e.g., unstable fusion, wave heating, or a binary companion) heats a region near the stellar surface. This appears to occur in luminous blue variables (LBVs), Type IIn supernovae progenitors, classical novae, and X-ray bursts. We show that when the wind kinetic power exceeds Eddington, the photons are trapped and behave like a fluid. Convection does not play a significant role in the wind energy transport. The wind properties depend on the ratio of a characteristic speed in the problem vc ~ (Edot G)^{1/5} (where Edot is the heating rate) to the stellar escape speed near the heating region vesc(r_h). For vc > vesc(r_h) the wind kinetic power at large radii Edot_w ~ Edot. For vc < vesc(r_h), most of the energy is used to unbind the wind material and thus Edot_w < Edot. Multidimensional hydrodynamic simulations without radiation diffusion using FLASH and one-dimensional hydrodynamic simulations with radiation diffusion using MESA are in good agreement with the analytic predictions. The photon luminosity from the wind is itself super-Eddington but in many cases the photon luminosity is likely dominated by `internal shocks' in the wind. We discuss the application of our models to eruptive mass loss from massive stars and argue that the wind models described here can account for the broad properties of LBV outflows and the enhanced mass loss in the years prior to Type IIn core-collapse supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Difficulty of Getting High Escape Fractions of Ionizing Photons from High-redshift Galaxies: a View from the FIRE Cosmological Simulations

We present a series of high-resolution (20-2000 Msun, 0.1-4 pc) cosmological zoom-in simulations at z~6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 10^9-10^11 Msun and rest-frame ultraviolet magnitude Muv = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01%-20%), while the time-averaged fesc over long time-scales generally remains ~5%, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3-10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g., models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that "sub-grid" star formation models, which do not explicitly resolve star formation in dense clouds with n >> 1 cm^-3, will dramatically over-predict fesc.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Diversity of Transients from Magnetar Birth

Strongly-magnetized, rapidly-rotating neutron stars are contenders for the central engines of both long-duration gamma-ray bursts (LGRBs) and hydrogen-poor super-luminous supernovae (SLSNe-I). Models for typical (~minute long) LGRBs invoke magnetars with high dipole magnetic fields (Bd > 1e15 G) and short spin-down times, while models for SLSNe-I invoke neutron stars with weaker fields and longer spin-down times of weeks. Here we identify a transition region in the space of Bd and birth period for which a magnetar can power both a long GRB and a luminous SN. In particular, we show that a 2 ms period magnetar with a spin-down time of ~1e4 s can explain the observations of both the ultra-long GRB 111209 and its associated luminous SN2011kl. For magnetars with longer spin down times, we predict even longer duration (~1e6 s) GRBs and brighter supernovae, a correlation that extends to Swift J2058+05 (commonly interpreted as a tidal disruption event). We further show that previous estimates of the maximum rotational energy of a proto-magnetar were too conservative and energies up to Emax ~1-2e53 erg are possible. The magnetar model can therefore comfortably accommodate the extreme energy requirements recently posed by the most luminous supernova ASASSN-15lh. The high ionization flux from a pulsar wind nebula powering ASASSN-15lh may lead to an "ionization break-out" X-ray burst over the coming months, which would be accompanied by an abrupt change in the optical spectrum. We conclude by briefly contrasting millisecond magnetar and black hole models for SLSNe and ultra-long GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The interplay of disk wind and dynamical ejecta in the aftermath of neutron star - black hole mergers

We explore the evolution of the different ejecta components generated during the merger of a neutron star (NS) and a black hole (BH). Our focus is the interplay between material ejected dynamically during the merger, and the wind launched on a viscous timescale by the remnant accretion disk. These components are expected to contribute to an electromagnetic transient and to produce r-process elements, each with a different signature when considered separately. Here we introduce a two-step approach to investigate their combined evolution, using two- and three-dimensional hydrodynamic simulations. Starting from the output of a merger simulation, we identify each component in the initial condition based on its phase space distribution, and evolve the accretion disk in axisymmetry. The wind blown from this disk is injected into a three-dimensional computational domain where the dynamical ejecta is evolved. We find that the wind can suppress fallback accretion on timescales longer than ~100 ms. Due to self-similar viscous evolution, the disk accretion at late times nevertheless approaches a power-law time dependence $\propto t^{-2.2}$. This can power some late-time GRB engine activity, although the available energy is significantly less than in traditional fallback models. Inclusion of radioactive heating due to the r-process does not significantly affect the fallback accretion rate or the disk wind. We do not find any significant modification to the wind properties at large radius due to interaction with the dynamical ejecta. This is a consequence of the different expansion velocities of the two components.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kilonova Light Curves from the Disk Wind Outflows of Compact Object Mergers

We study the radioactively-powered transients produced by accretion disk winds following a compact object merger. Starting with the outflows generated in two-dimensional hydrodynamical disk models, we use wavelength-dependent radiative transfer calculations to generate synthetic light curves and spectra. We show that the brightness and color of the resulting kilonova transients carry information about the merger physics. In the regions of the wind where neutrino irradiation raises the electron fraction to Ye > 0.25, r-process nucleosynthesis halts before producing high-opacity, complex ions (the lanthanides). The kilonova light curves thus show two distinct components: a brief (~2 day) blue optical transient produced in the outer lanthanide-free ejecta, and a longer (~10 day) infrared transient produced in the inner, lanthanide line-blanketed region. Mergers producing a longer-lived neutron star, or a more rapidly spinning black hole, have stronger neutrino irradiation, generate more lanthanide-free ejecta, and are optically brighter and bluer. At least some optical emission is produced in all disk wind models, which should enhance the detectability of electromagnetic counterparts to gravitational wave sources. However, the presence of even a small amount (10^{-4} Msun) of overlying, neutron-rich dynamical ejecta will act as a "lanthanide-curtain", obscuring the optical wind emission from certain viewing angles. Because the disk outflows have moderate velocities (~10,000 km/s), numerous resolved line features are discernible in the spectra, distinguishing disk winds from fast-moving dynamical ejecta, and offering a potential diagnostic of the detailed composition of freshly produced r-process material.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutron-powered precursors of kilonovae

The merger of binary neutron stars (NSs) ejects a small quantity of neutron rich matter, the radioactive decay of which powers a day to week long thermal transient known as a kilonova. Most of the ejecta remains sufficiently dense during its expansion that all neutrons are captured into nuclei during the r-process. However, recent general relativistic merger simulations by Bauswein and collaborators show that a small fraction of the ejected mass (a few per cent, or ~1e-4 Msun) expands sufficiently rapidly for most neutrons to avoid capture. This matter originates from the shocked-heated interface between the merging NSs. Here we show that the beta-decay of these free neutrons in the outermost ejecta powers a `precursor' to the main kilonova emission, which peaks on a timescale of a few hours following merger at U-band magnitude ~22 (for an assumed distance of 200 Mpc). The high luminosity and blue colors of the neutron precursor render it a potentially important counterpart to the gravitational wave source, that may encode valuable information on the properties of the merging binary (e.g. NS-NS versus NS-black hole) and the NS equation of state. Future work is necessary to assess the robustness of the fast moving ejecta and the survival of free neutrons in the face of neutrino absorptions, although the precursor properties are robust to a moderate amount of leptonization. Our results provide additional motivation for short latency gravitational wave triggers and rapid follow-up searches with sensitive ground based telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Outflows from accretion disks formed in neutron star mergers: effect of black hole spin

The accretion disk that forms after a neutron star merger is a source of neutron-rich ejecta. The ejected material contributes to a radioactively-powered electromagnetic transient, with properties that depend sensitively on the composition of the outflow. Here we investigate how the spin of the black hole remnant influences mass ejection on the thermal and viscous timescales. We carry out two-dimensional, time-dependent hydrodynamic simulations of merger remnant accretion disks including viscous angular momentum transport and approximate neutrino self-irradiation. The gravity of the spinning black hole is included via a pseudo-Newtonian potential. We find that a disk around a spinning black hole ejects more mass, up to a factor of several, relative to the non-spinning case. The enhanced mass loss is due to energy release by accretion occurring deeper in the gravitational potential, raising the disk temperature and hence the rate of viscous heating in regions where neutrino cooling is ineffective. The mean electron fraction of the outflow increases moderately with BH spin due to a highly-irradiated (though not neutrino-driven) wind component. While the bulk of the ejecta is still very neutron-rich, thus generating heavy r-process elements, the leading edge of the wind contains a small amount of Lanthanide-free material. This component can give rise to a ~1 day blue optical `bump' in a kilonova light curve, even in the case of prompt BH formation, which may facilitate its detection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Confronting Simulations of Optically Thick Gas in Massive Halos with Observations at z=2-3

Cosmological hydrodynamic simulations predict the physical state of baryons in the circumgalactic medium (CGM), which can be directly tested via quasar absorption line observations. We use high resolution "zoom-in" simulations of 21 galaxies to characterize the distribution of neutral hydrogen around halos in the mass range M_vir~ 2x10^11 - 4x10^12 M_sun at z~2. We find that both the mass fraction of cool (T <= 3x10^4 K) gas and the covering fraction of optically-thick Lyman limit systems (LLSs) depend only weakly on halo mass, even around the critical value for the formation of stable virial shocks. The covering fraction of LLSs interior to the virial radius varies between f_c~0.05 - 0.2, with significant scatter among halos. Our simulations of massive halos (M_vir >= 10^12 M_sun) underpredict the covering fraction of optically-thick gas observed in the quasar CGM by a large factor. The reason for this discrepancy is unclear, but several possibilities are discussed. In the lower mass halos (M_vir >= 5x10^11 M_sun) hosting star-forming galaxies, the predicted covering factor agrees with observations, however current samples of quasar-galaxy pairs are too small for a conclusive comparison. To overcome this limitation, we propose a new observable: the small-scale auto-correlation function of optically-thick absorbers detected in the foreground of close quasar pairs. We show that this new observable can constrain the underlying dark halos hosting LLSs at z~2-3, as well as the characteristic size and covering factor of the CGM.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effect of a High Opacity on the Light Curves of Radioactively Powered Transients from Compact Object Mergers

The coalescence of compact objects are a promising astrophysical sources of gravitational wave (GW) signals. The ejection of r-process material from such mergers may lead to a radioactively-powered electromagnetic counterpart which, if discovered, would enhance the science return of a GW detection. As very little is known about the optical properties of heavy r-process elements, previous light curve models have adopted opacities similar to those of iron group elements. Here we report that the presence of heavier elements, particularly the lanthanides, increase the ejecta opacity by several orders of magnitude. We include these higher opacities in time dependent, multi-wavelength radiative transport calculations to predict the broadband light curves of one-dimensional models over a range of parameters (ejecta masses from 0.001 to 0.1 solar masses and velocities from 0.1 to 0.3c). We find that the higher opacities lead to much longer duration light curves which can last a week or more. The emission is shifted toward the infrared bands due to strong optical line blanketing, and the colors at later times are representative of a blackbody near the recombination temperature of the lanthanides (T ~ 2500K). We further consider the case in which a second mass outflow, composed of 56Ni, is ejected from a disk wind, and show that the net result is a distinctive two component spectral energy distribution, with a bright optical peak due to 56Ni and an infrared peak due to r-process ejecta. We briefly consider the prospects for detection and identification of these transients.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nebular spectroscopy of the nearby type IIb supernova 2011dh

We present nebular spectra of the nearby Type IIb supernova (SN) 2011dh taken between 201 and 678 days after core collapse. At these late times, SN 2011dh exhibits strong emission lines including a very broad and persistent Hα feature. New models of the nebular spectra confirm that the progenitor of SN 2011dh was a low-mass giant (M ~ 13 - 15 M_sun) that ejected ~0.07 M_sun of 56Ni and ~0.27 M_sun of oxygen at the time of explosion, consistent with the recent disappearance of a candidate yellow supergiant progenitor. We show that light from the SN location is dominated by the fading SN at very late times (~2 yr) and not, for example, by a binary companion or a background source. We present evidence for interaction between the expanding SN blastwave and a circumstellar medium at late times and show that the SN is likely powered by positron deposition by ~1 yr after explosion. We also examine the geometry of the ejecta and show that SN 2011dh's nebular line profiles indicate a globally spherical explosion with aspherical components or clumps.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Opacities and Spectra of the r-process Ejecta from Neutron Star Mergers

Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements by the r-process. The subsequent radioactive decay of the nuclei can power electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying such events is an important goal of future transient surveys, offering new perspectives on the origin of r-process nuclei and the astrophysical sources of gravitational waves. Predictions of the transient light curves and spectra, however, have suffered from the uncertain optical properties of heavy ions. Here we consider the opacity of expanding r-process material and argue that it is dominated by line transitions from those ions with the most complex valence electron structure, namely the lanthanides. For a few representative ions, we run atomic structure models to calculate radiative data for tens of millions of lines. We find that the resulting r-process opacities are orders of magnitude larger than that of ordinary (e.g., iron-rich) supernova ejecta. Radiative transport calculations using these new opacities indicate that the transient emission should be dimmer and redder than previously thought. The spectra appear pseudo-blackbody, with broad absorption features, and peak in the infrared (~1 micron). We discuss uncertainties in the opacities and attempt to quantify their impact on the spectral predictions. The results have important implications for observational strategies to find and study the radioactively powered electromagnetic counterparts to compact object mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rapidly Fading Supernovae from Massive Star Explosions

Transient surveys have recently discovered a class of supernovae (SNe) with extremely rapidly declining light curves. These events are also often relatively faint, especially compared to Type Ia SNe. The common explanation for these events involves a weak explosion, producing a radioactive outflow with small ejected mass and kinetic energy (M ~ 0.1 Msun and E ~ 0.1 B, respectively), perhaps from the detonation of a helium shell on a white dwarf. We argue, in contrast, that these events may be Type Ib/c SNe with typical masses and energies (M ~ 3 Msun, E ~ 1 B), but which ejected very little radioactive material. In our picture, the light curve is powered by the diffusion of thermal energy deposited by the explosion shock wave, and the rapid evolution is due to recombination, which reduces the opacity and results in an "oxygen-plateau" light curve. Using a radiative transfer code, we generate synthetic spectra and light curves and demonstrate that this model can reasonably fit the observations of one event, SN 2010X. Similar models may explain the features of other rapidly evolving SNe such as SN 2002bj and SN 2005ek. SNe such as these require stripped-envelope progenitors with rather large radii (R ~ 20 Rsun), which may originate from a mass loss episode occurring just prior to explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Very Young Type Ia Supernova 2013dy: Discovery, and Strong Carbon Absorption in Early-Time Spectra

The Type Ia supernova (SN Ia) 2013dy in NGC 7250 (d ~ 13.7 Mpc) was discovered by the Lick Observatory Supernova Search. Combined with a prediscovery detection by the Italian Supernova Search Project, we are able to constrain the first-light time of SN 2013dy to be only 0.10 +/- 0.05 d (2.4 +/- 1.2 hr) before the first detection. This makes SN 2013dy the earliest known detection of an SN Ia. We infer an upper limit on the radius of the progenitor star of R_0 < 0.25 R_sun, consistent with that of a white dwarf. The light curve exhibits a broken power law with exponents of 0.88 and then 1.80. A spectrum taken 1.63 d after first light reveals a C II absorption line comparable in strength to Si II. This is the strongest C II feature ever detected in a normal SN Ia, suggesting that the progenitor star had significant unburned material. The C II line in SN 2013dy weakens rapidly and is undetected in a spectrum 7 days later, indicating that C II is detectable for only a very short time in some SNe Ia. SN 2013dy reached a B-band maximum of M_B = -18.72 +/- 0.03 mag ~17.7 d after first light.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Tidal Tail Ejection as a Signature of Type~Ia Supernovae from White Dwarf Mergers

The merger of two white dwarfs may be preceded by the ejection of some mass in "tidal tails", creating a circumstellar medium around the system. We consider the variety of observational signatures from this material, which depend on the lag time between the start of the merger and the ultimate explosion (assuming one occurs) of the system in a Type Ia supernova. If the time lag is fairly short, the interaction of the supernova ejecta with the tails could lead to detectable shock emission at radio, optical, and/or x-ray wavelengths. At somewhat later times, the tails produce relatively broad NaID absorption lines with velocity widths of order the white dwarf escape speed ($\sim 1000$ \kms). That none of these signatures have been detected in normal SNe Ia constrains the lag time to be either very short ($\lesssim 100$ s) or fairly long ($\gtrsim 100$ yr). If the tails have expanded and cooled over timescales $\sim 10^4$ yr, they could be observable through narrow NaID and CaII H&K absorption lines in the spectra, which are seen in some fraction of SNe Ia. Using a combination of 3D and 1D hydrodynamical codes, we model the mass-loss from tidal interactions in binary systems, and the subsequent interactions with the interstellar medium, which produce a slow-moving, dense shell of gas. We synthesize NaID line profiles by ray-casting through this shell, and show that in some circumstances tidal tails could be responsible for narrow absorptions similar to those observed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Type Ia Supernovae from Merging White Dwarfs I. Prompt Detonations

Merging white dwarfs are a possible progenitor of Type Ia supernovae (SNe Ia). While it is not entirely clear if and when an explosion is triggered in such systems, numerical models suggest that a detonation might be initiated before the stars have coalesced to form a single compact object. Here we study such "peri-merger" detonations by means of numerical simulations, modeling the disruption and nucleosynthesis of the stars until the ejecta reach the coasting phase. Synthetic light curves and spectra are generated for comparison with observations. Three models are considered with primary masses 0.96 Msun, 1.06 Msun, and 1.20 Msun. Of these, the 0.96 Msun dwarf merging with an 0.81 Msun companion, with a Ni56 yield of 0.58 Msun, is the most promising candidate for reproducing common SNe Ia. The more massive mergers produce unusually luminous SNe Ia with peak luminosities approaching those attributed to "super-Chandrasekhar" mass SNe Ia. While the synthetic light curves and spectra of some of the models resemble observed SNe Ia, the significant asymmetry of the ejecta leads to large orientation effects. The peak bolometric luminosity varies by more than a factor of 2 with the viewing angle, and the velocities of the spectral absorption features are lower when observed from angles where the light curve is brightest. The largest orientation effects are seen in the ultraviolet, where the flux varies by more than an order of magnitude. Despite the large variation with viewing angle, the set of three models roughly obeys a width-luminosity relation, with the brighter light curves declining more slowly in the B-band. Spectral features due to unburned carbon from the secondary star are also seen in some cases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Type Ia Supernovae from Merging White Dwarfs II) Post-Merger Detonations

Merging carbon-oxygen (CO) white dwarfs are a promising progenitor system for Type Ia supernovae (SN Ia), but the underlying physics and timing of the detonation are still debated. If an explosion occurs after the secondary star is fully disrupted, the exploding primary will expand into a dense CO medium that may still have a disk-like structure. This interaction will decelerate and distort the ejecta. Here we carry out multi-dimensional simulations of ``tamped" SN Ia models, using both particle and grid-based codes to study the merger and explosion dynamics, and a radiative transfer code to calculate synthetic spectra and light curves. We find that post-merger explosions exhibit an hourglass-shaped asymmetry, leading to strong variations in the light curves with viewing angle. The two most important factors affecting the outcome are the scale-height of the disk, which depends sensitively on the binary mass ratio, and the total ${}^{56}$Ni yield, which is governed by the central density of the remnant core. The synthetic broadband light curves rise and decline very slowly, and the spectra generally look peculiar, with weak features from intermediate mass elements but relatively strong carbon absorption. We also consider the effects of the viscous evolution of the remnant, and show that a longer time delay between merger and explosion probably leads to larger ${}^{56}$Ni yields and more symmetrical remnants. We discuss the relevance of this class of aspherical ``tamped" SN Ia for explaining the class of ``super-Chandrasekhar'' SN Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

3-D Radiative Transfer Calculations of Radiation Feedback from Massive Black Holes: Outflow of Mass from the Dusty "Torus"

Observational and theoretical arguments suggest that the momentum carried in mass outflows from AGN can reach several times L / c, corresponding to outflow rates of hundreds of solar masses per year. Radiation pressure on lines alone may not be sufficient to provide this momentum deposition, and the transfer of reprocessed IR radiation in dusty nuclear gas has been postulated to provide the extra enhancement. The efficacy of this mechanism, however, will be sensitive to multi-dimensional effects such as the tendency for the reprocessed radiation to preferentially escape along sight-lines of lower column density. We use Monte Carlo radiative transfer calculations to determine the radiation force on dusty gas residing within approximately 10 parsecs from an accreting super-massive black hole. We calculate the net rate of momentum deposition in the surrounding gas and estimate the mass-loss rate in the resulting outflow as a function of solid angle for different black hole luminosities, sightline-averaged column densities, clumping parameters, and opening angles of the dusty gas. We find that these dust-driven winds carry momentum fluxes of 1-5 times L / c and correspond to mass-loss rates of 10-100 solar masses per year for a 10^8 solar mass black hole radiating at or near its Eddington limit. These results help to explain the origin of high velocity molecular and atomic outflows in local ULIRGs, and can inform numerical simulations of galaxy evolution including AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pair-Instability Supernovae at the Epoch of Reionization

Pristine stars with masses between ~140 and 260 M_sun are theoretically predicted to die as pair-instability supernovae. These very massive progenitors could come from Pop III stars in the early universe. We model the light curves and spectra of pair-instability supernovae over a range of masses and envelope structures. At redshifts of reionization z >= 6, we calculate the rates and detectability of pair-instability and core collapse supernovae, and show that with the James Webb Space Telescope, it is possible to determine the contribution of Pop III and Pop II stars toward reionization by constraining the stellar initial mass function at that epoch using these supernovae. We also find the rates of Type Ia supernovae, and show that they are not rare during reionization, and can be used to probe the mass function at 4-8 M_sun. If the budget of ionizing photons was dominated by contributions from top-heavy Pop III stars, we predict that the bright end of the galaxy luminosity function will be contaminated by pair-instability supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pair-Instability Supernovae via Collision Runaway in Young Dense Star Clusters

Stars with helium cores between ~64 and 133 M_sun are theoretically predicted to die as pair-instability supernovae. This requires very massive progenitors, which are theoretically prohibited for Pop II/I stars within the Galactic stellar mass limit due to mass loss via line-driven winds. However, the runaway collision of stars in a dense, young star cluster could create a merged star with sufficient mass to end its life as a pair-instability supernova, even with enhanced mass loss at non-zero metallicity. We show that the predicted rate from this mechanism is consistent with the inferred volumetric rate of roughly ~2x10^-9 Mpc^-3 yr^-1 of the two observed pair-instability supernovae, SN 2007bi and PTF 10nmn, neither of which have metal-free host galaxies. Contrary to prior literature, only pair-instability supernovae at low redshifts z<2 will be observable with the Large Synoptic Survey Telescope (LSST). We estimate the telescope will observe ~10^2 such events per year that originate from the collisional runaway mergers in clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Supernova Light Curves Powered by Fallback Accretion

Some fraction of the material ejected in a core collapse supernova explosion may remain bound to the compact remnant, and eventually turn around and fall back. We show that the late time (> days) power associated with the accretion of this "fallback" material may significantly affect the optical light curve, in some cases producing super-luminous or otherwise peculiar supernovae. We use spherically symmetric hydrodynamical models to estimate the accretion rate at late times for a range of progenitor masses and radii and explosion energies. The accretion rate onto the proto-neutron star or black hole decreases as Mdot ~ t^-5/3 at late times, but its normalization can be significantly enhanced at low explosion energies, in very massive stars, or if a strong reverse shock wave forms at the helium/hydrogen interface in the progenitor. If the resulting super-Eddington accretion drives an outflow which thermalizes in the outgoing ejecta, the supernova debris will be re-energized at a time when photons can diffuse out efficiently. The resulting light curves are different and more diverse than previous fallback supernova models which ignored the input of accretion power and produced short-lived, dim transients. The possible outcomes when fallback accretion power is significant include super-luminous (> 10^44 ergs / s) Type II events of both short and long durations, as well as luminous Type I events from compact stars that may have experienced significant mass loss. Accretion power may unbind the remaining infalling material, causing a sudden decrease in the brightness of some long duration Type II events. This scenario may be relevant for explaining some of the recently discovered classes of peculiar and rare supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Long-Term Evolution of Double White Dwarf Mergers

In this paper, we present a model for the long-term evolution of the merger of two unequal mass C/O white dwarfs (WDs). After the dynamical phase of the merger, magnetic stresses rapidly redistribute angular momentum, leading to nearly solid-body rotation on a viscous timescale of 1e-4 to 1 yr, long before significant cooling can occur. Due to heating during the dynamical and viscous phases, the less massive WD is transformed into a hot, slowly rotating, and radially extended envelope supported by thermal pressure. Following the viscous phase of evolution, the maximum temperature near the envelope base may already be high enough to begin off-center convective carbon-burning. If not, Kelvin-Helmholtz contraction of the inner region of the envelope on a thermal timescale of 1e3-1e4 yr compresses the base of the envelope, again yielding off-center burning. As a result, the long-term evolution of the merger remnant is similar to that seen in previous calculations: the burning shell diffuses inwards over ~1e4 yr, eventually yielding a high-mass O/Ne WD or a collapse to a neutron star. During the cooling and shell-burning phases, the merger remnant radiates near the Eddington limit. Given the double WD merger rate of a few per 1000 yr, tens of these ~1e38 erg/s sources should exist in a Milky Way-type galaxy. While the end result is similar to that of previous studies, the physical picture and the dynamical state of the matter in our model differ from previous work. Furthermore, remaining uncertainties related to the convective structure near the photosphere and mass loss during the thermal evolution may significantly affect our conclusions. Thus, future work within the context of the physical model presented here is required to better address the eventual fate of double WD mergers, including those for which one or both of the components is a He WD.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Compact Degenerate Primary-Star Progenitor of SN 2011fe

While a white dwarf is, from a theoretical perspective, the most plausible primary star in Type Ia supernova (SN Ia), many other candidates have not been formally ruled out. Shock energy deposited in the envelope of any exploding primary contributes to the early SN brightness and, since this radiation energy is degraded by expansion after the explosion, the diffusive luminosity depends on the initial primary radius. We present a new non-detection limit of the nearby SN Ia 2011fe, obtained what appears to be just 4 hours after explosion, allowing us to directly constrain the initial primary radius, R_p. Coupled with the non-detection of a quiescent X-ray counterpart and the inferred synthesized Ni mass, we show that R_p <~ 0.02 R_sun (a factor of 5 smaller than previously inferred), that the average density of the primary must be rho_p > 10,000 gm cm^{-3}, and that the effective temperature must be less than a few x 10^5 K. This rules out hydrogen burning main sequence stars and giants. Constructing the helium-burning main sequence and carbon-burning main sequence, we find such objects are also excluded. By process of elimination, we find that only degeneracy-supported compact objects---WDs and neutron stars---are viable as the primary star of SN 2011fe. With few caveats, we also restrict the companion (secondary) star radius to R_ c <~ 0.1 R_sun, excluding Roche-Lobe overflowing red giant and main-sequence companions to high significance.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Absorption line systems in simulated galaxies fed by cold streams

Hydro cosmological simulations reveal that massive galaxies at high redshift are fed by long narrow streams of merging galaxies and a smoother component of cold gas. We post-process seven high-resolution simulated galaxies with radiative transfer to study the absorption characteristics of the gas in galaxies and streams, in comparison with the statistics of observed absorption-line systems. We find that much of the stream gas is ionized by UV radiation from background and local stellar sources, but still optically thick (N_HI > 10^17 cm^-2) so that the streams appear as Lyman-limit systems (LLSs). At z>3, the fraction of neutral gas in streams becomes non-negligible, giving rise to damped Lyman-alpha (DLA) absorbers as well. The gas in the central and incoming galaxies remains mostly neutral, responsible for DLAs. Within one (two) virial radii, the covering factor of optically thick gas is <25% (10%) for LLSs and <5% (1%) for DLAs, slowly declining with time following the universal expansion. Nevertheless, galaxies and their cold streams in the studied mass range, M_vir = 10^10 - 10^12 Msun, account for >30% of the observed absorbers in the foreground of quasars, the rest possibly arising from smaller galaxies or the intergalactic medium. The mean metallicity in the streams is ~1% solar, much lower than in the galaxies. The simulated galaxies reproduce the Lyalpha-absorption equivalent widths observed around Lyman-break galaxies, but they severely underpredict the equivalent widths in metal lines, suggesting that the latter may arise from outflows. We conclude that the observed metal-poor LLSs are likely detections of the predicted cold streams. Revised analysis of the observed LLSs kinematics and simulations with more massive outflows in conjunction with the inflows may enable a clearer distinction between the signatures of the various gas modes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pair Instability Supernovae: Light Curves, Spectra, and Shock Breakout

For the initial mass range (140 < M < 260 Msun) stars die in a thermonuclear runaway triggered by the pair-production instability. The supernovae they make can be remarkably energetic (up to ~10^53 ergs) and synthesize considerable amounts of radioactive isotopes. Here we model the evolution, explosion, and observational signatures of representative pair-instability supernovae (PI SNe) spanning a range of initial masses and envelope structures. The predicted light curves last for hundreds of days and range in luminosity, from very dim to extremely bright, L ~ 10^44 ergs/s. The most massive events are bright enough to be seen at high redshift, but the extended light curve duration (~1 year) -- prolonged by cosmological time-dilation -- may make it difficult to detect them as transients. An alternative approach may be to search for the brief and luminous outbreak occurring when the explosion shock wave reaches the stellar surface. Using a multi-wavelength radiation-hydrodynamics code we calculate that, in the rest-frame, the shock breakout transients of PI SNe reach luminosities of 10^45-10^46 ergs/s, peak at wavelengths ~30-170 Angstroms, and last for several hours. We explore the detectability of PI SNe emission at high redshift, and discuss how observations of the light curves, spectra, and breakout emission can be used to constrain the mass, radius, and metallicity of the progenitor.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Peculiar Type II Supernovae from Blue Supergiants

The vast majority of Type II supernovae (SNe) are produced by red supergiants (RSGs), but SN 1987A revealed that blue supergiants (BSGs) can produce members of this class as well, albeit with some peculiar properties. This best studied event revolutionized our understanding of SNe, and linking it to the bulk of Type II events is essential. We present here optical photometry and spectroscopy gathered for SN 2000cb, which is clearly not a standard Type II SN and yet is not a SN 1987A analog. The light curve of SN 2000cb is reminiscent of that of SN 1987A in shape, with a slow rise to a late optical peak, but on substantially different time scales. Spectroscopically, SN 2000cb resembles a normal SN II but with ejecta velocities that far exceed those measured for SN 1987A or normal SNe II, above 18000 km/s for H-alpha at early times. The red colours, high velocities, late photometric peak, and our modeling of this object all point toward a scenario involving the high-energy explosion of a small-radius star, most likely a BSG, producing 0.1 solar masses of Ni-56. Adding a similar object to the sample, SN 2005ci, we derive a rate of about 2% of the core-collapse rate for this loosely defined class of BSG explosions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simple Models of Metal-Line Absorption and Emission from Cool Gas Outflows

We analyze the absorption and emission-line profiles produced by a set of simple, cool gas wind models motivated by galactic-scale outflow observations. We implement monte carlo radiative transfer techniques that track the propagation of scattered and fluorescent photons to generate 1D spectra and 2D spectral images. We focus on the MgII 2796,28303 doublet and FeII UV1 multiplet at ~2600A, but the results are applicable to other transitions that trace outflows (e.g. NaI, Lya, SiII). By design, the resonance transitions show blue-shifted absorption but one also predicts strong resonance and fine-structure line-emission at roughly the systemic velocity. This line-emission `fills-in' the absorption reducing the equivalent width by up to 50%, shift the absorption-lin centroid by tens of km/s, and reduce the effective opacity near systemic. Analysis of cool gas outflows that ignores this line-emission may incorrectly infer that the gas is partially covered, measure asignificantly lower peak optical depth, and/or conclude that gas at systemic velocity is absent. Because the FeII lines are connected by optically-thin transitions to fine-structure levels, their profiles more closely reproduce the intrinsic opacity of the wind. Together these results naturally explain the absorption and emission-line characteristics observed for star-forming galaxies at z<1. We also study a scenario promoted to describe the outflows of z~3 Lyman break galaxies and find prfiles inconsistent with the observations due to scattered photon emission. Although line-emission complicates the analysis of absorption-line profiles, the surface brightness profiles offer a unique means of assessing the morphology and size of galactic-scale winds. Furthermore, the kinematics and line-ratios offer powerful diagnostics of outflows, motivating deep, spatially-extended spectroscopic observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sub-Chandrasekhar Mass Models For Type Ia Supernovae

For carbon-oxygen white dwarfs accreting hydrogen or helium at rates in the range ~1-10 x 10^(-8) Msun/y, a variety of explosive outcomes is possible well before the star reaches the Chandrasekhar mass. These outcomes are surveyed for a range of white dwarf masses (0.7 - 1.1 Msun), accretion rates (1 - 7 x 10^(-8) Msun/y), and initial white dwarf temperatures (0.01 and 1 Lsun). The results are particularly sensitive to the convection that goes on during the last few minutes before the explosion. Unless this convection maintains a shallow temperature gradient, and unless the density is sufficiently high, the accreted helium does not detonate. Below a critical helium ignition density, which we estimate to be 5 - 10 x 10^5 g cm^(-3), either helium novae or helium deflagrations result. The hydrodynamics, nucleosynthesis, light curves, and spectra of a representative sample of detonating and deflagrating models are explored. Some can be quite faint indeed, powered at peak for a few days by the decay of 48Cr and 48V. Only the hottest, most massive white dwarfs considered with the smallest helium layers, show reasonable agreement with the light curves and spectra of common Type Ia supernovae. For the other models, especially those involving lighter white dwarfs, the helium shell mass exceeds 0.05 Msun and the mass of the 56Ni that is synthesized exceeds 0.01 Msun. These explosions do not look like ordinary Type Ia supernovae, or any other frequently observed transient.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Supernova 2011fe from an Exploding Carbon-Oxygen White Dwarf Star

Type Ia supernovae (SNe Ia) have been used empirically as standardized candles to reveal the accelerating universe even though fundamental details, such as the nature of the progenitor system and how the star explodes, remained a mystery. There is consensus that a white dwarf star explodes after accreting matter in a binary system, but the secondary could be anything from a main sequence star to a red giant, or even another white dwarf. The uncertainty stems from the fact that no recent SN Ia has been discovered close enough to detect the stars before explosion. Here we report early observations of SN 2011fe (PTF11kly) in M101 at a distance of 6.4 Mpc, the closest SN Ia in the past 25 years. We find that the exploding star was likely a carbon-oxygen white dwarf, and from the lack of an early shock we conclude that the companion was most likely a main sequence star. Early spectroscopy shows high-velocity oxygen that varies on a time scale of hours and extensive mixing of newly synthesized intermediate mass elements in the outermost layers of the supernova. A companion paper uses pre-explosion images to rule out luminous red giants and most helium stars as companions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Most Slowly Declining Type Ia Supernova 2001ay

We present optical and near-infrared photometry, as well as ground-based optical spectra and Hubble Space Telescope ultraviolet spectra, of the Type Ia supernova (SN) 2001ay. At maximum light the Si II and Mg II lines indicated expansion velocities of 14,000 km/sec, while Si III and S II showed velocities of 9,000 km/sec There is also evidence for some unburned carbon at 12,000 km/sec. SN 2001ay exhibited a decline-rate parameter Delta m_15(B) = 0.68 \pm 0.05 mag; this and the B-band photometry at t > +25 d past maximum make it the most slowly declining Type Ia SN yet discovered. Three of four super-Chandrasekhar-mass candidates have decline rates almost as slow as this. After correction for Galactic and host-galaxy extinction, SN 2001ay had M_B = -19.19 and M_V = -19.17 mag at maximum light; thus, it was not overluminous in optical bands. In near-infrared bands it was overluminous only at the 2-sigma level at most. For a rise time of 18 d (explosion to bolometric maximum) the implied Ni-56 yield was (0.58 \pm 0.15)/alpha M_Sun, with alpha = L_max/E_Ni probably in the range 1.0 to 1.2. The Ni-56 yield is comparable to that of many Type Ia supernovae. The "normal" Ni-56 yield and the typical peak optical brightness suggest that the very broad optical light curve is explained by the trapping of the gamma rays in the inner regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Low-ionization Line Emission from Starburst Galaxies: A New Probe of Galactic-Scale Outflows

We study the kinematically narrow, low-ionization line emission from a bright, starburst galaxy at z = 0.69 using slit spectroscopy obtained with Keck/LRIS. The spectrum reveals strong absorption in MgII and FeII resonance transitions with Doppler shifts of -200 to -300 km/s, indicating a cool gas outflow. Emission in MgII near and redward of systemic velocity, in concert with the observed absorption, yields a P Cygni-like line profile similar to those observed in the Ly alpha transition in Lyman Break Galaxies. Further, the MgII emission is spatially resolved, and extends significantly beyond the emission from stars and HII regions within the galaxy. Assuming the emission has a simple, symmetric surface brightness profile, we find that the gas extends to distances > ~7 kpc. We also detect several narrow FeII* fine-structure lines in emission near the systemic velocity, arising from energy levels which are radiatively excited directly from the ground state. We suggest that the MgII and FeII* emission is generated by photon scattering in the observed outflow, and emphasize that this emission is a generic prediction of outflows. These observations provide the first direct constraints on the minimum spatial extent and morphology of the wind from a distant galaxy. Estimates of these parameters are crucial for understanding the impact of outflows in driving galaxy evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measuring Ejecta Velocity Improves Type Ia Supernova Distances

We use a sample of 121 spectroscopically normal Type Ia supernovae (SNe Ia) to show that their intrinsic color is correlated with their ejecta velocity, as measured from the blueshift of the Si II 6355 feature near maximum brightness, v_Si. The SN Ia sample was originally used by Wang et al. (2009) to show that the relationship between color excess and peak magnitude, which in the absence of intrinsic color differences describes a reddening law, was different for two subsamples split by v_Si (defined as "Normal" and "High-Velocity"). We verify this result, but find that the two subsamples have the same reddening law when extremely reddened events (E(B-V) > 0.35 mag) are excluded. We also show that (1) the High-Velocity subsample is offset by ~0.06 mag to the red from the Normal subsample in the (B_max - V_max) - M_V plane, (2) the B_max - V_max cumulative distribution functions of the two subsamples have nearly identical shapes, but the High-Velocity subsample is offset by ~0.07 mag to the red in B_max - V_max, and (3) the bluest High-Velocity SNe Ia are ~0.10 mag redder than the bluest Normal SNe Ia. Together, this evidence indicates a difference in intrinsic color for the subsamples. Accounting for this intrinsic color difference reduces the scatter in Hubble residuals from 0.190 mag to 0.130 mag for SNe Ia with A_V < 0.7 mag. The scatter can be further reduced to 0.109 mag by exclusively using SNe Ia from the Normal subsample. Additionally, this result can at least partially explain the anomalously low values of R_V found in large SN Ia samples. We explain the correlation between ejecta velocity and color as increased line blanketing in the High-Velocity SNe Ia, causing them to become redder. We discuss some implications of this result, and stress the importance of spectroscopy for future SN Ia cosmology surveys, with particular focus on the design of WFIRST.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Single or Double Degenerate Progenitors? Searching for Shock Emission in the SDSS-II Type Ia Supernovae

From the set of nearly 500 spectroscopically confirmed type~Ia supernovae and around 10,000 unconfirmed candidates from SDSS-II, we select a subset of 108 confirmed SNe Ia with well-observed early-time light curves to search for signatures from shock interaction of the supernova with a companion star. No evidence for shock emission is seen; however, the cadence and photometric noise could hide a weak shock signal. We simulate shocked light curves using SN Ia templates and a simple, Gaussian shock model to emulate the noise properties of the SDSS-II sample and estimate the detectability of the shock interaction signal as a function of shock amplitude, shock width, and shock fraction. We find no direct evidence for shock interaction in the rest-frame $B$-band, but place an upper limit on the shock amplitude at 9% of supernova peak flux ($M_B > -16.6$ mag). If the single degenerate channel dominates type~Ia progenitors, this result constrains the companion stars to be less than about 6 $M_{\odot}$ on the main sequence, and strongly disfavors red giant companions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Rise and Fall of Type Ia Supernova Light Curves in the SDSS-II Supernova Survey

We analyze the rise and fall times of type Ia supernova (SN Ia) light curves discovered by the SDSS-II Supernova Survey. From a set of 391 light curves k-corrected to the rest frame B and V bands, we find a smaller dispersion in the rising portion of the light curve compared to the decline. This is in qualitative agreement with computer models which predict that variations in radioactive nickel yield have less impact on the rise than on the spread of the decline rates. The differences we find in the rise and fall properties suggest that a single 'stretch' correction to the light curve phase does not properly model the range of SN Ia light curve shapes. We select a subset of 105 light curves well-observed in both rise and fall portions of the light curves and develop a '2-stretch' fit algorithm which estimates the rise and fall times independently. We find the average time from explosion to B-band peak brightness is 17.38 +/- 0.17 days. Our average rise time is shorter than the 19.5 days found in previous studies; this reflects both the different light curve template used and the application of the 2-stretch algorithm. We find that slow declining events tend to have fast rise times, but that the distribution of rise minus fall time is broad and single-peaked. This distribution is in contrast to the bimodality in this parameter that was first suggested by Strovink (2007) from an analysis of a small set of local SNe Ia. We divide the SDSS-II sample in half based on the rise minus fall value, tr-tf <= 2 days and tr-tf>2 days, to search for differences in their host galaxy properties and Hubble residuals; we find no difference in host galaxy properties or Hubble residuals in our sample.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermonuclear .Ia Supernovae from Helium Shell Detonations: Explosion Models and Observables

During the early evolution of an AM CVn system, helium is accreted onto the surface of a white dwarf under conditions suitable for unstable thermonuclear ignition. The turbulent motions induced by the convective burning phase in the He envelope become strong enough to influence the propagation of burning fronts and may result in the onset of a detonation. Such an outcome would yield radioactive isotopes and a faint rapidly rising thermonuclear ".Ia" supernova. In this paper, we present hydrodynamic explosion models and observable outcomes of these He shell detonations for a range of initial core and envelope masses. The peak UVOIR bolometric luminosities range by a factor of 10 (from 5e41 - 5e42 erg/s), and the R-band peak varies from M_R,peak = -15 to -18. The rise times in all bands are very rapid (<10 d), but the decline rate is slower in the red than the blue due to a secondary near-IR brightening. The nucleosynthesis primarily yields heavy alpha-chain elements (40Ca through 56Ni) and unburnt He. Thus, the spectra around peak light lack signs of intermediate mass elements and are dominated by CaII and TiII features, with the caveat that our radiative transfer code does not include the non-thermal effects necessary to produce He features.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Optical Transients from the Unbound Debris of Tidal Disruption

In the tidal disruption of a star by a black hole, roughly half of the stellar mass becomes bound and falls into the black hole, while the other half is ejected at high velocity. Several previous studies have considered the emission resulting from the accretion of bound material; we consider the possibility that the unbound debris may also radiate once it has expanded and become transparent. We show that the gradual energy input from hydrogen recombination compensates for adiabatic loses over significant expansion factors. The opacity also drops dramatically with recombination, and the internal energy can be radiated by means of a cooling-transparency wave propagating from the surface layers inward. The result is a brief optical transient occurring ~1 week after disruption and lasting 3-5 days with peak luminosities of 10^40-10^42 ergs/s, depending on the mass of the disrupted star. These recombination powered transients should accompany the x-ray/ultraviolet flare from the accretion of bound material, and so may be a useful signature for discriminating tidal disruption events, especially for lower and intermediate mass black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Seeing the Collision of a Supernova with its Companion Star

The progenitors of Type Ia and some core collapse supernovae are thought to be stars in binary systems, but little observational evidence exists to confirm the hypothesis. We suggest that the collision of the supernova ejecta with its companion star should produce detectable emission in the hours and days following the explosion. The interaction occurs at distances ~10^11-10^13 cm and shocks the impacting supernova debris, dissipating kinetic energy and re-heating the gas. Initially, some radiation may escape promptly through the evacuated region of the shadowcone, producing a bright X-ray (0.1-2 keV) burst lasting minutes to hours with luminosity L ~ 10^44 ergs/s. Continuing radiative diffusion from deeper layers of shock heated ejecta produces a longer lasting optical/UV emission which exceeds the radioactively powered luminosity of the supernova for the first few days after the explosion. These signatures are prominent for viewing angles looking down upon the shocked region, or about 10% of the time. The properties of the emission provide a straightforward measure of the separation distance between the stars and hence (assuming Roche lobe overflow) the companion's radius. Current optical and UV data sets likely already constrain red giant companions. By systematically acquiring early time data for many supernovae, it should eventually be possible to empirically determine how the parameters of the progenitor system influence the outcome of the explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Supernova Light Curves Powered by Young Magnetars

We show that energy deposited into an expanding supernova remnant by a highly magnetic (B ~ 5 x 10^14 G) neutron star spinning at an initial period of P ~ 2-20 ms can substantially brighten the light curve. For magnetars with parameters in this range, the rotational energy is released on a timescale of days to weeks, which is comparable to the effective diffusion time through the supernova remnant. The late time energy injection can then be radiated without suffering overwhelming adiabatic expansion losses. The magnetar input also produces a central bubble which sweeps ejecta into an internal dense shell, resulting in a prolonged period of nearly constant photospheric velocity in the observed spectra. We derive analytic expressions for the light curve rise time and peak luminosity as a function of B, P and the properties of the supernova ejecta that allow for direct inferences about the underlying magnetar in bright supernovae. We perform numerical radiation hydrodynamical calculations of a few specific instances and compare the resulting light curves to observed events. Magnetar activity is likely to impact more than a few percent of all core collapse supernovae, and may naturally explain some of the brightest events ever seen (e.g., SN 2005ap and SN 2008es) at L > 10^44 ergs/s.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Diversity of Type Ia Supernovae from Broken Symmetries

Type Ia supernovae result when carbon-oxygen white dwarfs in binary systems accrete mass from companion stars, reach a critical mass, and explode. The near uniformity of their light curves makes these supernovae good standard candles for measuring cosmic expansion, but a correction must be applied to account for the fact that the brighter supernovae have broader light curves. One-dimensional modelling, with a certain choice of parameters, can reproduce this general trend in the width-luminosity relation, but the processes of ignition and detonation have recently been shown to be intrinsically asymmetric. Here we report on multi-dimensional modelling of the explosion physics and radiative transfer that reveals that the breaking of spherical symmetry is a critical factor in determining both the width luminosity relation and the observed scatter about it. The deviation from sphericity can also explain the finite polarization detected in the light from some supernovae. The slope and normalization of the width-luminosity relation has a weak dependence on certain properties of the white dwarf progenitor, in particular the trace abundances of elements other than carbon and oxygen. Failing to correct for this effect could lead to systematic overestimates of up to 2% in the distance to remote supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Three-Dimensional Simulations of Tidally Disrupted Solar-Type Stars and the Observational Signatures of Shock Breakout

We describe a three-dimensional simulation of a $1 M_{\odot}$ solar-type star approaching a $10^{6} M_{\odot}$ black hole on a parabolic orbit with a pericenter distance well within the tidal radius. While falling towards the black hole, the star is not only stretched along the orbital direction but even more severely compressed at right angles to the orbit. The overbearing degree of compression achieved shortly after pericenter leads to the production of strong shocks which largely homogenize the temperature profile of the star, resulting in surface temperatures comparable to the initial temperature of the star's core. This phenomenon, which precedes the fallback accretion phase, gives rise to a unique double-peaked X-ray signature that, if detected, may be one of the few observable diagnostics of how stars behave under the influence of strong gravitational fields. If $\sim 10^{6} M_{\odot}$ black holes were prevalent in small or even dwarf galaxies, the nearest of such flares may be detectable by EXIST from no further away than the Virgo Cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Type II Supernovae: Model Light Curves and Standard Candle Relationships

A survey of Type II supernovae explosion models has been carried out to determine how their light curves and spectra vary with their mass, metallicity, and explosion energy. The presupernova models are taken from a recent survey of massive stellar evolution at solar metallicity supplemented by new calculations at subsolar metallicity. Explosions are simulated by the motion of a piston near the edge of the iron core and the resulting light curves and spectra are calculated using full multi-wavelength radiation transport. Formulae are developed that describe approximately how the model observables (light curve luminosity and duration) scale with the progenitor mass, explosion energy, and radioactive nucleosynthesis. Comparison with observational data shows that the explosion energy of typical supernovae (as measured by kinetic energy at infinity) varies by nearly an order of magnitude -- from 0.5 to 4.0 x 10^51 ergs, with a typical value of ~0.9 x 10^51 ergs. Despite the large variation, the models exhibit a tight relationship between luminosity and expansion velocity, similar to that previously employed empirically to make SNe IIP standardized candles. This relation is explained by the simple behavior of hydrogen recombination in the supernova envelope, but we find a sensitivity to progenitor metallicity and mass that could lead to systematic errors. Additional correlations between light curve luminosity, duration, and color might enable the use of SNe IIP to obtain distances accurate to ~20% using only photometric data.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Detonating Failed Deflagration Model of Thermonuclear Supernovae II. Comparison to Observations

We develop and demonstrate the methodology of testing multi-dimensional supernova models against observations by studying the properties of one example of the detonation from failed deflagration (DFD) explosion model of thermonuclear supernovae. Using time-dependent multi-dimensional radiative transfer calculations, we generate the synthetic broadband optical light curves, near-infrared light curves, color evolution curves, full spectral time-series, and spectropolarization of the model, as seen from various viewing angles. All model observables are critically evaluated against examples of well-observed, standard Type Ia supernovae (SNe Ia). We explore the consequences of the intrinsic model asphericity by studying the dependence of the model emission on viewing angle, and by quantifying the resulting dispersion in (and internal correlations between) various model observables. These statistical properties of the model are also evaluated against those of the available observational sample of SNe Ia. On the whole, the DFD model shows good agreement with a broad range of SN Ia observations. Certain deficiencies are also apparent, and point to further developments within the basic theoretical framework. We also identify several intriguing orientation effects in the model which suggest ways in which the asphericity of SNe Ia may contribute to their photometric and spectroscopic diversity and, conversely, how the relative homogeneity of SNe Ia constrains the degree of asymmetry allowable in the models. The comprehensive methodology adopted in this work proves an essential component of developing and validating theoretical supernova models, and helps motivate and clearly define future directions in both the modeling and the observation of SNe Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Type Ia Supernova Light Curves

The diversity of Type Ia supernova (SN Ia) photometry is explored using a grid of 130 one-dimensional models. It is shown that the observable properties of SNe Ia resulting from Chandrasekhar-mass explosions are chiefly determined by their final composition and some measure of ``mixing'' in the explosion. A grid of final compositions is explored including essentially all combinations of 56Ni, stable ``iron'', and intermediate mass elements that result in an unbound white dwarf. Light curves (and in some cases spectra) are calculated for each model using two different approaches to the radiation transport problem. Within the resulting templates are models that provide good photometric matches to essentially the entire range of observed SNe Ia. On the whole, the grid of models spans a wide range in B-band peak magnitudes and decline rates, and does not obey a Phillips relation. In particular, models with the same mass of 56Ni show large variations in their light curve decline rates. We identify the physical parameters responsible for this dispersion, and consider physically motivated ``cuts'' of the models that agree better with the Phillips relation. For example, models that produce a constant total mass of burned material of 1.1 +/- Msun do give a crude Phillips relation, albeit with much scatter. The scatter is further reduced if one restricts that set to models that make 0.1 to 0.3 Msun of stable iron and nickel isotopes, and then mix the ejecta strongly between the center and 0.8 Msun. We conclude that the supernovae that occur most frequently in nature are highly constrained by the Phillips relation and that a large part of the currently observed scatter in the relation is likely a consequence of the intrinsic diversity of these objects.

0 citations0 reviewsNo code linkedOpen access