-
Meeting the Challenges of Modeling Astrophysical Thermonuclear Explosions: Castro, Maestro, and the AMReX Astrophysics Suite
Authors:
M. Zingale,
A. S. Almgren,
M. G. Barrios Sazo,
V. E. Beckner,
J. B. Bell,
B. Friesen,
A. M. Jacobs,
M. P. Katz,
C. M. Malone,
A. J. Nonaka,
D. E. Willcox,
W. Zhang
Abstract:
We describe the AMReX suite of astrophysics codes and their application to modeling problems in stellar astrophysics. Maestro is tuned to efficiently model subsonic convective flows while Castro models the highly compressible flows associated with stellar explosions. Both are built on the block-structured adaptive mesh refinement library AMReX. Together, these codes enable a thorough investigation…
▽ More
We describe the AMReX suite of astrophysics codes and their application to modeling problems in stellar astrophysics. Maestro is tuned to efficiently model subsonic convective flows while Castro models the highly compressible flows associated with stellar explosions. Both are built on the block-structured adaptive mesh refinement library AMReX. Together, these codes enable a thorough investigation of stellar phenomena, including Type Ia supernovae and X-ray bursts. We describe these science applications and the approach we are taking to make these codes performant on current and future many-core and GPU-based architectures.
△ Less
Submitted 16 November, 2017;
originally announced November 2017.
-
Galactos: Computing the Anisotropic 3-Point Correlation Function for 2 Billion Galaxies
Authors:
Brian Friesen,
Md. Mostofa Ali Patwary,
Brian Austin,
Nadathur Satish,
Zachary Slepian,
Narayanan Sundaram,
Deborah Bard,
Daniel J Eisenstein,
Jack Deslippe,
Pradeep Dubey,
Prabhat
Abstract:
The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to bill…
▽ More
The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to billions of galaxies. We present Galactos, a high-performance implementation of a novel, O(N^2) algorithm that uses a load-balanced k-d tree and spherical harmonic expansions to compute the anisotropic 3PCF. Our implementation is optimized for the Intel Xeon Phi architecture, exploiting SIMD parallelism, instruction and thread concurrency, and significant L1 and L2 cache reuse, reaching 39% of peak performance on a single node. Galactos scales to the full Cori system, achieving 9.8PF (peak) and 5.06PF (sustained) across 9636 nodes, making the 3PCF easily computable for all galaxies in the observable universe.
△ Less
Submitted 31 August, 2017;
originally announced September 2017.
-
Optical and ultraviolet spectroscopic analysis of SN 2011fe at late times
Authors:
Brian Friesen,
E. Baron,
Jerod T. Parrent,
R. C. Thomas,
David Branch,
Peter Nugent,
Peter H. Hauschildt,
Ryan J. Foley,
Darryl E. Wright,
Yen-Chen Pan,
Alexei V. Filippenko,
Kelsey I. Clubb,
Jeffrey M. Silverman,
Keiichi Maeda,
Isaac Shivvers,
Patrick L. Kelly,
Daniel P. Cohen,
Armin Rest,
Daniel Kasen
Abstract:
We present optical spectra of the nearby Type Ia supernova SN 2011fe at 100, 205, 311, 349, and 578 days post-maximum light, as well as an ultraviolet spectrum obtained with Hubble Space Telescope at 360 days post-maximum light. We compare these observations with synthetic spectra produced with the radiative transfer code PHOENIX. The day +100 spectrum can be well fit with models which neglect col…
▽ More
We present optical spectra of the nearby Type Ia supernova SN 2011fe at 100, 205, 311, 349, and 578 days post-maximum light, as well as an ultraviolet spectrum obtained with Hubble Space Telescope at 360 days post-maximum light. We compare these observations with synthetic spectra produced with the radiative transfer code PHOENIX. The day +100 spectrum can be well fit with models which neglect collisional and radiative data for forbidden lines. Curiously, including this data and recomputing the fit yields a quite similar spectrum, but with different combinations of lines forming some of the stronger features. At day +205 and later epochs, forbidden lines dominate much of the optical spectrum formation; however, our results indicate that recombination, not collisional excitation, is the most influential physical process driving spectrum formation at these late times. Consequently, our synthetic optical and UV spectra at all epochs presented here are formed almost exclusively through recombination-driven fluorescence. Furthermore, our models suggest that the ultraviolet spectrum even as late as day +360 is optically thick and consists of permitted lines from several iron-peak species. These results indicate that the transition to the "nebular" phase in Type Ia supernovae is complex and highly wavelength-dependent.
△ Less
Submitted 16 July, 2016;
originally announced July 2016.
-
Comparative Analysis of SN 2012dn Optical Spectra: Days -14 to +114
Authors:
J. T. Parrent,
D. A. Howell,
R. A. Fesen,
S. Parker,
F. B. Bianco,
B. Dilday,
D. Sand,
S. Valenti,
J. Vinkó,
P. Berlind,
P. Challis,
D. Milisavljevic,
N. Sanders,
G. H. Marion,
J. C. Wheeler,
P. Brown,
M. L. Calkins,
B. Friesen,
R. Kirshner,
T. Pritchard,
R. Quimby,
P. Roming
Abstract:
SN 2012dn is a super-Chandrasekhar mass candidate in a purportedly normal spiral (SAcd) galaxy, and poses a challenge for theories of type Ia supernova diversity. Here we utilize the fast and highly parameterized spectrum synthesis tool, SYNAPPS, to estimate relative expansion velocities of species inferred from optical spectra obtained with six facilities. As with previous studies of normal SN Ia…
▽ More
SN 2012dn is a super-Chandrasekhar mass candidate in a purportedly normal spiral (SAcd) galaxy, and poses a challenge for theories of type Ia supernova diversity. Here we utilize the fast and highly parameterized spectrum synthesis tool, SYNAPPS, to estimate relative expansion velocities of species inferred from optical spectra obtained with six facilities. As with previous studies of normal SN Ia, we find that both unburned carbon and intermediate mass elements are spatially coincident within the ejecta near and below 14,000 km/s. Although the upper limit on SN 2012dn's peak luminosity is comparable to some of the most luminous normal SN Ia, we find a progenitor mass exceeding ~1.6 Msun is not strongly favored by leading merger models since these models do not accurately predict spectroscopic observations of SN 2012dn and more normal events. In addition, a comparison of light curves and host-galaxy masses for a sample of literature and Palomar Transient Factory SN Ia reveals a diverse distribution of SN Ia subtypes where carbon-rich material remains unburned in some instances. Such events include SN 1991T, 1997br, and 1999aa where trace signatures of C III at optical wavelengths are presumably detected.
△ Less
Submitted 12 March, 2016;
originally announced March 2016.
-
Spectral Models for Early Time SN 2011fe Observations
Authors:
E. Baron,
P. Hoeflich Brian Friesen,
M. Sullivan,
E. Hsiao,
R. S. Ellis,
A. Gal-Yam,
D. A. Howell,
P. E. Nugent,
I. Dominguez,
K. Krisciunas,
M. M. Phillips,
N. Suntzeff,
L. Wang,
R. C. Thomas
Abstract:
We use observed UV through near IR spectra to examine whether SN 2011fe can be understood in the framework of Branch-normal SNe Ia and to examine its individual peculiarities. As a benchmark, we use a delayed-detonation model with a progenitor metallicity of Z_solar/20. We study the sensitivity of features to variations in progenitor metallicity, the outer density profile, and the distribution of…
▽ More
We use observed UV through near IR spectra to examine whether SN 2011fe can be understood in the framework of Branch-normal SNe Ia and to examine its individual peculiarities. As a benchmark, we use a delayed-detonation model with a progenitor metallicity of Z_solar/20. We study the sensitivity of features to variations in progenitor metallicity, the outer density profile, and the distribution of radioactive nickel. The effect of metallicity variations in the progenitor have a relatively small effect on the synthetic spectra. We also find that the abundance stratification of SN 2011fe resembles closely that of a delayed detonation model with a transition density that has been fit to other Branch-normal Type Ia supernovae. At early times, the model photosphere is formed in material with velocities that are too high, indicating that the photosphere recedes too slowly or that SN 2011fe has a lower specific energy in the outer ~0.1 M_sun than does the model. We discuss several explanations for the discrepancies. Finally, we examine variations in both the spectral energy distribution and in the colors due to variations in the progenitor metallicity, which suggests that colors are only weak indicators for the progenitor metallicity, in the particular explosion model that we have studied. We do find that the flux in the U band is significantly higher at maximum light in the solar metallicity model than in the lower metallicity model and the lower metallicity model much better matches the observed spectrum.
△ Less
Submitted 31 August, 2015; v1 submitted 19 August, 2015;
originally announced August 2015.
-
Near-infrared line identification in type Ia supernovae during the transitional phase
Authors:
Brian Friesen,
E. Baron,
John P. Wisniewski,
Jerod T. Parrent,
R. C. Thomas,
Timothy R. Miller,
G. H. Marion
Abstract:
We present near-infrared synthetic spectra of a delayed-detonation hydrodynamical model and compare them to observed spectra of four normal type Ia supernovae ranging from day +56.5 to day +85. This is the epoch during which supernovae are believed to be undergoing the transition from the photospheric phase, where spectra are characterized by line scattering above an optically thick photosphere, t…
▽ More
We present near-infrared synthetic spectra of a delayed-detonation hydrodynamical model and compare them to observed spectra of four normal type Ia supernovae ranging from day +56.5 to day +85. This is the epoch during which supernovae are believed to be undergoing the transition from the photospheric phase, where spectra are characterized by line scattering above an optically thick photosphere, to the nebular phase, where spectra consist of optically thin emission from forbidden lines. We find that most spectral features in the near-infrared can be accounted for by permitted lines of Fe II and Co II. In addition, we find that [Ni II] fits the emission feature near 1.98 μm, suggesting that a substantial mass of 58Ni exists near the center of the ejecta in these objects, arising from nuclear burning at high density. A tentative identification of Mn II at 1.15 μm may support this conclusion as well.
△ Less
Submitted 29 July, 2014;
originally announced July 2014.
-
A review of type Ia supernova spectra
Authors:
J. Parrent,
B. Friesen,
M. Parthasarathy
Abstract:
SN 2011fe was the nearest and best-observed type Ia supernova in a generation, and brought previous incomplete datasets into sharp contrast with the detailed new data. In retrospect, documenting spectroscopic behaviors of type Ia supernovae has been more often limited by sparse and incomplete temporal sampling than by consequences of signal-to-noise ratios, telluric features, or small sample sizes…
▽ More
SN 2011fe was the nearest and best-observed type Ia supernova in a generation, and brought previous incomplete datasets into sharp contrast with the detailed new data. In retrospect, documenting spectroscopic behaviors of type Ia supernovae has been more often limited by sparse and incomplete temporal sampling than by consequences of signal-to-noise ratios, telluric features, or small sample sizes. As a result, type Ia supernovae have been primarily studied insofar as parameters discretized by relative epochs and incomplete temporal snapshots near maximum light. Here we discuss a necessary next step toward consistently modeling and directly measuring spectroscopic observables of type Ia supernova spectra. In addition, we analyze current spectroscopic data in the parameter space defined by empirical metrics, which will be relevant even after progenitors are observed and detailed models are refined.
△ Less
Submitted 25 February, 2014;
originally announced February 2014.
-
Supernova Resonance--scattering Line Profiles in the Absence of a Photosphere
Authors:
Brian Friesen,
E. Baron,
David Branch,
Bin Chen,
Jerod T. Parrent,
R. C. Thomas
Abstract:
In supernova spectroscopy relatively little attention has been given to the properties of optically thick spectral lines in epochs following the photosphere's recession. Most treatments and analyses of post-photospheric optical spectra of supernovae assume that forbidden-line emission comprises most if not all spectral features. However, evidence exists which suggests that some spectra exhibit lin…
▽ More
In supernova spectroscopy relatively little attention has been given to the properties of optically thick spectral lines in epochs following the photosphere's recession. Most treatments and analyses of post-photospheric optical spectra of supernovae assume that forbidden-line emission comprises most if not all spectral features. However, evidence exists which suggests that some spectra exhibit line profiles formed via optically thick resonance-scattering even months or years after the supernova explosion. To explore this possibility we present a geometrical approach to supernova spectrum formation based on the "Elementary Supernova" model, wherein we investigate the characteristics of resonance-scattering in optically thick lines while replacing the photosphere with a transparent central core emitting non-blackbody continuum radiation, akin to the optical continuum provided by decaying 56Co formed during the explosion. We develop the mathematical framework necessary for solving the radiative transfer equation under these conditions, and calculate spectra for both isolated and blended lines. Our comparisons with analogous results from the Elementary Supernova code SYNOW reveal several marked differences in line formation. Most notably, resonance lines in these conditions form P Cygni-like profiles, but the emission peaks and absorption troughs shift redward and blueward, respectively, from the line's rest wavelength by a significant amount, despite the spherically symmetric distribution of the line optical depth in the ejecta. These properties and others that we find in this work could lead to misidentification of lines or misattribution of properties of line-forming material at post-photospheric times in supernova optical spectra.
△ Less
Submitted 16 September, 2012;
originally announced September 2012.
-
Analysis of the Early-Time Optical Spectra of SN 2011fe in M101
Authors:
J. T. Parrent,
D. A. Howell,
B. Friesen,
R. C. Thomas,
R. A. Fesen,
D. Milisavljevic,
F. B. Bianco,
B. Dilday,
P. Nugent,
E. Baron,
I. Arcavi,
S. Ben-Ami,
D. Bersier,
L. Bildsten,
J. Bloom,
Y. Cao,
S. B. Cenko,
A. V. Filippenko,
A. Gal-Yam,
M. M. Kasliwal,
N. Konidaris,
S. R. Kulkarni,
N. M. Law,
D. Levitan,
K. Maguire
, et al. (13 additional authors not shown)
Abstract:
The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture…
▽ More
The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.
△ Less
Submitted 30 May, 2012; v1 submitted 27 May, 2012;
originally announced May 2012.