-
The $\texttt{Abacus}$ Cosmological $N$-body Code
Authors:
Lehman H. Garrison,
Daniel J. Eisenstein,
Douglas Ferrer,
Nina A. Maksimova,
Philip A. Pinto
Abstract:
We present $\texttt{Abacus}$, a fast and accurate cosmological $N$-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct $1/r^2$ summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per sec…
▽ More
We present $\texttt{Abacus}$, a fast and accurate cosmological $N$-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct $1/r^2$ summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of $10^{-5}$. We express the simulation time step as an event-driven "pipeline", incorporating asynchronous events such as completion of co-processor work, Input/Output, and network communication. $\texttt{Abacus}$ has been used to produce the largest suite of $N$-body simulations to date, the $\texttt{AbacusSummit}$ suite of 60 trillion particles (Maksimova et al., 2021), incorporating on-the-fly halo finding. $\texttt{Abacus}$ enables the production of mock catalogs of the volume and resolution required by the coming generation of cosmological surveys.
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
A High-Fidelity Realization of the Euclid Code Comparison $N$-body Simulation with Abacus
Authors:
Lehman H. Garrison,
Daniel J. Eisenstein,
Philip A. Pinto
Abstract:
We present a high-fidelity realization of the cosmological $N$-body simulation from the Schneider et al. (2016) code comparison project. The simulation was performed with our Abacus $N$-body code, which offers high force accuracy, high performance, and minimal particle integration errors. The simulation consists of $2048^3$ particles in a $500\ h^{-1}\mathrm{Mpc}$ box, for a particle mass of…
▽ More
We present a high-fidelity realization of the cosmological $N$-body simulation from the Schneider et al. (2016) code comparison project. The simulation was performed with our Abacus $N$-body code, which offers high force accuracy, high performance, and minimal particle integration errors. The simulation consists of $2048^3$ particles in a $500\ h^{-1}\mathrm{Mpc}$ box, for a particle mass of $1.2\times 10^9\ h^{-1}\mathrm{M}_\odot$ with $10\ h^{-1}\mathrm{kpc}$ spline softening. Abacus executed 1052 global time steps to $z=0$ in 107 hours on one dual-Xeon, dual-GPU node, for a mean rate of 23 million particles per second per step. We find Abacus is in good agreement with Ramses and Pkdgrav3 and less so with Gadget3. We validate our choice of time step by halving the step size and find sub-percent differences in the power spectrum and 2PCF at nearly all measured scales, with $<0.3\%$ errors at $k<10\ \mathrm{Mpc}^{-1}h$. On large scales, Abacus reproduces linear theory better than $0.01\%$. Simulation snapshots are available at http://nbody.rc.fas.harvard.edu/public/S2016 .
△ Less
Submitted 6 March, 2019; v1 submitted 5 October, 2018;
originally announced October 2018.
-
The Abacus Cosmos: A Suite of Cosmological N-body Simulations
Authors:
Lehman H. Garrison,
Daniel J. Eisenstein,
Douglas Ferrer,
Jeremy L. Tinker,
Philip A. Pinto,
David H. Weinberg
Abstract:
We present a public data release of halo catalogs from a suite of 125 cosmological $N$-body simulations from the Abacus project. The simulations span 40 $w$CDM cosmologies centered on the Planck 2015 cosmology at two mass resolutions, $4\times 10^{10}\;h^{-1}M_\odot$ and $1\times 10^{10}\;h^{-1}M_\odot$, in $1.1\;h^{-1}\mathrm{Gpc}$ and $720\;h^{-1}\mathrm{Mpc}$ boxes, respectively. The boxes are…
▽ More
We present a public data release of halo catalogs from a suite of 125 cosmological $N$-body simulations from the Abacus project. The simulations span 40 $w$CDM cosmologies centered on the Planck 2015 cosmology at two mass resolutions, $4\times 10^{10}\;h^{-1}M_\odot$ and $1\times 10^{10}\;h^{-1}M_\odot$, in $1.1\;h^{-1}\mathrm{Gpc}$ and $720\;h^{-1}\mathrm{Mpc}$ boxes, respectively. The boxes are phase-matched to suppress sample variance and isolate cosmology dependence. Additional volume is available via 16 boxes of fixed cosmology and varied phase; a few boxes of single-parameter excursions from Planck 2015 are also provided. Catalogs spanning $z=1.5$ to $0.1$ are available for friends-of-friends and Rockstar halo finders and include particle subsamples. All data products are available at https://lgarrison.github.io/AbacusCosmos
△ Less
Submitted 23 April, 2018; v1 submitted 15 December, 2017;
originally announced December 2017.
-
Improving Initial Conditions for Cosmological $N$-Body Simulations
Authors:
Lehman H. Garrison,
Daniel J. Eisenstein,
Douglas Ferrer,
Marc V. Metchnik,
Philip A. Pinto
Abstract:
In cosmological $N$-body simulations, the representation of dark matter as discrete "macroparticles" suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near $k_{\rm Nyquist}$. Marcos et al. demonstrate that this is due to sparse sampling of modes near $k_{\rm Nyquist}$ and that the often-assumed continuum growing modes are not proper growing…
▽ More
In cosmological $N$-body simulations, the representation of dark matter as discrete "macroparticles" suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near $k_{\rm Nyquist}$. Marcos et al. demonstrate that this is due to sparse sampling of modes near $k_{\rm Nyquist}$ and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate $N$-body code Abacus to implement 2LPT in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1% accuracy in the power spectrum down to $k_{\rm Nyquist}$, versus $k_{\rm Nyquist}/4$ without rescaling or $k_{\rm Nyquist}/13$ without 2LPT, relative to an oversampled reference simulation. We anticipate that our 2LPT will be useful for large simulations where FFTs are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogs. Code to generate initial conditions is available at https://github.com/lgarrison/zeldovich-PLT
△ Less
Submitted 3 August, 2016; v1 submitted 8 May, 2016;
originally announced May 2016.
-
LSST Science Book, Version 2.0
Authors:
LSST Science Collaboration,
Paul A. Abell,
Julius Allison,
Scott F. Anderson,
John R. Andrew,
J. Roger P. Angel,
Lee Armus,
David Arnett,
S. J. Asztalos,
Tim S. Axelrod,
Stephen Bailey,
D. R. Ballantyne,
Justin R. Bankert,
Wayne A. Barkhouse,
Jeffrey D. Barr,
L. Felipe Barrientos,
Aaron J. Barth,
James G. Bartlett,
Andrew C. Becker,
Jacek Becla,
Timothy C. Beers,
Joseph P. Bernstein,
Rahul Biswas,
Michael R. Blanton,
Joshua S. Bloom
, et al. (223 additional authors not shown)
Abstract:
A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south…
▽ More
A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.
△ Less
Submitted 1 December, 2009;
originally announced December 2009.
-
LSST: from Science Drivers to Reference Design and Anticipated Data Products
Authors:
Željko Ivezić,
Steven M. Kahn,
J. Anthony Tyson,
Bob Abel,
Emily Acosta,
Robyn Allsman,
David Alonso,
Yusra AlSayyad,
Scott F. Anderson,
John Andrew,
James Roger P. Angel,
George Z. Angeli,
Reza Ansari,
Pierre Antilogus,
Constanza Araujo,
Robert Armstrong,
Kirk T. Arndt,
Pierre Astier,
Éric Aubourg,
Nicole Auza,
Tim S. Axelrod,
Deborah J. Bard,
Jeff D. Barr,
Aurelian Barrau,
James G. Bartlett
, et al. (288 additional authors not shown)
Abstract:
(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the…
▽ More
(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pachón in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$σ$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $δ<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.
△ Less
Submitted 23 May, 2018; v1 submitted 15 May, 2008;
originally announced May 2008.
-
The Carnegie Supernova Project: The Low-Redshift Survey
Authors:
Mario Hamuy,
Gastón Folatelli,
Nidia I. Morrell,
Mark M. Phillips,
Nicholas B. Suntzeff,
S. E. Persson,
Miguel Roth,
Sergio Gonzalez,
Wojtek Krzeminski,
Carlos Contreras,
Wendy L. Freedman,
D. C. Murphy,
Barry F. Madore,
P. Wyatt,
José Maza,
Alexei V. Filippenko,
Weidong Li,
P. A. Pinto
Abstract:
Supernovae are essential to understanding the chemical evolution of the Universe. Type Ia supernovae also provide the most powerful observational tool currently available for studying the expansion history of the Universe and the nature of dark energy. Our basic knowledge of supernovae comes from the study of their photometric and spectroscopic properties. However, the presently available data s…
▽ More
Supernovae are essential to understanding the chemical evolution of the Universe. Type Ia supernovae also provide the most powerful observational tool currently available for studying the expansion history of the Universe and the nature of dark energy. Our basic knowledge of supernovae comes from the study of their photometric and spectroscopic properties. However, the presently available data sets of optical and near-infrared light curves of supernovae are rather small and/or heterogeneous, and employ photometric systems that are poorly characterized. Similarly, there are relatively few supernovae whose spectral evolution has been well sampled, both in wavelength and phase, with precise spectrophotometric observations. The low-redshift portion of the Carnegie Supernova Project (CSP) seeks to remedy this situation by providing photometry and spectrophotometry of a large sample of supernovae taken on telescope/filter/detector systems that are well understood and well characterized. During a five-year program which began in September 2004, we expect to obtain high-precision u'g'r'i'BVYJHKs light curves and optical spectrophotometry for about 250 supernovae of all types. In this paper we provide a detailed description of the CSP survey observing and data reduction methodology. In addition, we present preliminary photometry and spectra obtained for a few representative supernovae during the first observing campaign.
△ Less
Submitted 1 December, 2005;
originally announced December 2005.
-
Probing Galaxy Formation with He II Cooling Lines
Authors:
Yujin Yang,
Ann I. Zabludoff,
Romeel Davé,
Daniel J. Eisenstein,
Philip A. Pinto,
Neal Katz,
David H. Weinberg,
Elizabeth J. Barton
Abstract:
Using high resolution cosmological simulations, we study hydrogen and helium gravitational cooling radiation. We focus on the HeII cooling lines, which arise from gas with a different temperature history (T_max ~ 10^5K) than HI line emitting gas. We examine whether three major atomic cooling lines, HI 1216A, HeII 1640A and HeII 304A are observable, finding that HI 1216A and HeII 1640A cooling em…
▽ More
Using high resolution cosmological simulations, we study hydrogen and helium gravitational cooling radiation. We focus on the HeII cooling lines, which arise from gas with a different temperature history (T_max ~ 10^5K) than HI line emitting gas. We examine whether three major atomic cooling lines, HI 1216A, HeII 1640A and HeII 304A are observable, finding that HI 1216A and HeII 1640A cooling emission at z=2-3 are potentially detectable with deep narrow band (R>100) imaging and/or spectroscopy from the ground. While the expected strength of HI 1216A cooling emission depends strongly on the treatment of the self-shielded phase of the IGM in the simulations, our predictions for the HeII 1640A line are more robust because the HeII 1640A emissivity is negligible below T~10^4.5 K and less sensitive to the UV background. Although HeII 1640A cooling emission is fainter than HI 1216A by at least a factor of 10 and, unlike HI 1216A, might not be resolved spatially with current observational facilities, it is more suitable to study gas accretion in the galaxy formation process because it is optically thin and less contaminated by the recombination lines from star-forming galaxies. The HeII 1640A line can be used to distinguish among mechanisms for powering the so-called "Lyman alpha blobs" -- including gravitational cooling radiation, photoionization by stellar populations, and starburst-driven superwinds -- because (1) HeII 1640A emission is limited to very low metallicity (log(Z/Z_sun) < -5.3) and Population III stars, and (2) the blob's kinematics are probed unambiguously through the HeII 1640A line width, which, for cooling radiation, is narrower (sigma < 400 km/s) than typical wind speeds.
△ Less
Submitted 31 August, 2005;
originally announced September 2005.
-
Optical and Infrared Photometry of the Type Ia Supernovae 1991T, 1991bg, 1999ek, 2001bt, 2001cn, 2001cz, and 2002bo
Authors:
Kevin Krisciunas,
Nicholas B. Suntzeff,
Mark M. Phillips,
Pablo Candia,
Jose Luis Prieto,
R. Antezana,
R. Chassagne,
H. -W. Chen,
M. Dickinson,
P. R. Eisenhardt,
J. Espinoza,
P. M. Garnavich,
D. Gonzalez,
T. E. Harrison,
M. Hamuy,
V. D. Ivanov,
W. Krzeminski,
C. Kulesa,
P. McCarthy,
A. Moro-Martin,
C. Muena,
A. Noriega-Crespo,
S. E. Persson,
P. A. Pinto,
M. Roth
, et al. (6 additional authors not shown)
Abstract:
We present optical and/or infrared photometry of the Type Ia supernovae SN 1991T, SN 1991bg, SN 1999ek, SN 2001bt, SN 2001cn, SN 2001cz, and SN 2002bo. All but one of these supernovae have decline rate parameters Delta m_15(B) close to the median value of 1.1 for the whole class of Type Ia supernovae. The addition of these supernovae to the relationship between the near-infrared absolute magnitu…
▽ More
We present optical and/or infrared photometry of the Type Ia supernovae SN 1991T, SN 1991bg, SN 1999ek, SN 2001bt, SN 2001cn, SN 2001cz, and SN 2002bo. All but one of these supernovae have decline rate parameters Delta m_15(B) close to the median value of 1.1 for the whole class of Type Ia supernovae. The addition of these supernovae to the relationship between the near-infrared absolute magnitudes and Delta m_15(B) strengthens the previous relationships we have found, in that the maximum light absolute magnitudes are essentially independent of the decline rate parameter. (SN 1991bg, the prototype of the subclass of fast declining Type Ia supernovae, is a special case.) The dispersion in the Hubble diagram in JHK is only ~0.15 mag. The near-infrared properties of Type Ia supernovae continue to be excellent measures of the luminosity distances to the supernova host galaxies, due to the need for only small corrections from the epoch of observation to maximum light, low dispersion in absolute magnitudes at maximum light, and the minimal reddening effects in the near-infrared.
△ Less
Submitted 17 December, 2004; v1 submitted 1 September, 2004;
originally announced September 2004.
-
Unified 1-D Simulations of Gamma-Ray Line Emission from Type Ia Supernovae
Authors:
P. A. Milne,
A. L. Hungerford,
C. L. Fryer,
T. M. Evans,
T. J. Urbatsch,
S. E. Boggs,
J. Isern,
E. Bravo,
A. Hirschmann,
S. Kumagai,
P. A. Pinto,
L. -S. The
Abstract:
The light curves of Type Ia Supernovae (SN Ia) are powered by gamma-rays emitted by the decay of radioactive elements such as $^{56}$Ni and its decay products. These gamma-rays are downscattered,absorbed, and eventually reprocessed into the optical emission which makes up the bulk of all supernova observations. Detection of the gamma-rays that escape the expanding star provide the only direct me…
▽ More
The light curves of Type Ia Supernovae (SN Ia) are powered by gamma-rays emitted by the decay of radioactive elements such as $^{56}$Ni and its decay products. These gamma-rays are downscattered,absorbed, and eventually reprocessed into the optical emission which makes up the bulk of all supernova observations. Detection of the gamma-rays that escape the expanding star provide the only direct means to study this power source for SN Ia light curves. Unfortunately, disagreements between calculations for the gamma-ray lines have made it difficult to interpret any gamma-ray observations.
Here we present a detailed comparison of the major gamma-ray line transport codes for a series of 1-dimensional Ia models. Discrepancies in past results were due to errors in the codes, and the corrected versions of the seven different codes yield very similar results. This convergence of the simulation results allows us to infer more reliable information from the current set of gamma-ray observations of SNe Ia. The observations of SNe 1986G, 1991T and 1998bu are consistent with explosion models based on their classification: sub-luminous, super-luminous and normally-luminous respectively.
△ Less
Submitted 7 June, 2004;
originally announced June 2004.
-
Evidence for Core Collapse in the Type Ib/c SN 1999ex
Authors:
Mario Hamuy,
Maximilian Stritzinger,
M. M. Phillips,
Nicholas B. Suntzeff,
Jose Maza,
Philip A. Pinto
Abstract:
We present optical and infrared spectra of SN 1999ex that reveal a clear example of an intermediate Type Ib/c case. This suggests a continuous spectroscopic sequence between Type Ib and Type Ic supernovae. We report UBVRIz photometric observations of SN 1999ex that started only one day after explosion, which permitted us to witness an elusive transient cooling phase that lasted 4 days. The initi…
▽ More
We present optical and infrared spectra of SN 1999ex that reveal a clear example of an intermediate Type Ib/c case. This suggests a continuous spectroscopic sequence between Type Ib and Type Ic supernovae. We report UBVRIz photometric observations of SN 1999ex that started only one day after explosion, which permitted us to witness an elusive transient cooling phase that lasted 4 days. The initial cooling and subsequent radioactive heating produced a dip in the lightcurve which is consistent with explosion models involving core collapse of evolved massive helium stars, and inconsistent with lightcurves resulting from the thermonuclear runaway of compact white dwarfs.
△ Less
Submitted 16 December, 2002;
originally announced December 2002.
-
Shock Breakout in Core-Collapse Supernovae and its Neutrino Signature
Authors:
Todd A. Thompson,
Adam Burrows,
Philip A. Pinto
Abstract:
(Abridged) We present results from dynamical models of core-collapse supernovae in one spatial dimension, employing a newly-developed Boltzmann neutrino radiation transport algorithm, coupled to Lagrangean hydrodynamics and a consistent high-density nuclear equation of state. We focus on shock breakout and its neutrino signature and follow the dynamical evolution of the cores of 11 M_sun, 15 M_s…
▽ More
(Abridged) We present results from dynamical models of core-collapse supernovae in one spatial dimension, employing a newly-developed Boltzmann neutrino radiation transport algorithm, coupled to Lagrangean hydrodynamics and a consistent high-density nuclear equation of state. We focus on shock breakout and its neutrino signature and follow the dynamical evolution of the cores of 11 M_sun, 15 M_sun, and 20 M_sun progenitors through collapse and the first 250 milliseconds after bounce. We examine the effects on the emergent neutrino spectra, light curves, and mix of species of artificial opacity changes, the number of energy groups, the weak magnetism/recoil corrections, nucleon-nucleon bremsstrahlung, neutrino-electron scattering, and the compressibility of nuclear matter. Furthermore, we present the first high-resolution look at the angular distribution of the neutrino radiation field both in the semi-transparent regime and at large radii and explore the accuracy with which our tangent-ray method tracks the free propagation of a pulse of radiation in a near vacuum. Finally, we fold the emergent neutrino spectra with the efficiencies and detection processes for a selection of modern underground neutrino observatories and argue that the prompt electron-neutrino breakout burst from the next galactic supernova is in principle observable and usefully diagnostic of fundamental collapse/supernova behavior. Though we are not in this study focusing on the supernova mechanism per se, our simulations support the theoretical conclusion (already reached by others) that spherical (1D) supernovae do not explode when good physics and transport methods are employed.
△ Less
Submitted 2 April, 2003; v1 submitted 10 November, 2002;
originally announced November 2002.
-
Optical and Infrared Spectroscopy of SN 1999ee and SN 1999ex
Authors:
Mario Hamuy,
Jose Maza,
Philip A. Pinto,
M. M. Phillips,
Nicholas B. Suntzeff,
R. D. Blum,
K. A. G. Olsen,
David J. Pinfield,
Valentin D. Ivanov,
T. Augusteijn,
S. Brillant,
M. Chadid,
J. -G. Cuby,
V. Doublier,
O. R. Hainaut,
E. Le Floc'h,
C. Lidman,
Monika G. Petr-Gotzens,
E. Pompei,
L. Vanzi
Abstract:
We report optical and infrared spectroscopic observations of the Type Ia SN 1999ee and the Type Ib/c SN 1999ex, both of which were hosted by the galaxy IC 5179. For SN 1999ee we obtained a continuous sequence with an unprecedented wavelength and temporal coverage beginning 9 days before maximum light and extending through day 42. Before maximum light SN 1999ee displayed a normal spectrum with a…
▽ More
We report optical and infrared spectroscopic observations of the Type Ia SN 1999ee and the Type Ib/c SN 1999ex, both of which were hosted by the galaxy IC 5179. For SN 1999ee we obtained a continuous sequence with an unprecedented wavelength and temporal coverage beginning 9 days before maximum light and extending through day 42. Before maximum light SN 1999ee displayed a normal spectrum with a strong Si II 6355 absorption, thus showing that not all slow-declining SNe are spectroscopically peculiar at these evolutionary phases. A comparative study of the infrared spectra of SN 1999ee and other Type Ia supernovae shows that there is a remarkable homogeneity among the Branch-normal SNe Ia during their first 60 days of evolution. SN 1991bg-like objects, on the other hand, display spectroscopic peculiarities at infrared wavelengths. SN 1999ex was characterized by the lack of hydrogen lines, weak optical He I lines, and strong He I 10830,20581, thus providing an example of an intermediate case between pure Ib and Ic supernovae. We conclude therefore that SN 1999ex provides first clear evidence for a link between the Ib and Ic classes and that there is a continuous spectroscopic sequence ranging from the He deficient SNe Ic to the SNe Ib which are characterized by strong optical He I lines.
△ Less
Submitted 27 March, 2002;
originally announced March 2002.
-
Type II Supernovae as Standardized Candles
Authors:
Mario Hamuy,
Philip A. Pinto
Abstract:
We present evidence for a correlation between expansion velocities of the ejecta of Type II plateau supernovae and their bolometric luminosities during the plateau phase. This correlation permits one to standardize the candles and decrease the scatter in the Hubble diagram from ~1 mag to a level of 0.4 and 0.3 mag in the V and I bands, respectively. When we restrict the sample to the eight objec…
▽ More
We present evidence for a correlation between expansion velocities of the ejecta of Type II plateau supernovae and their bolometric luminosities during the plateau phase. This correlation permits one to standardize the candles and decrease the scatter in the Hubble diagram from ~1 mag to a level of 0.4 and 0.3 mag in the V and I bands, respectively. When we restrict the sample to the eight objects which are well in the Hubble flow (cz > 3,000 km/s) the scatter drops even further to only 0.2 mag (or 9% in distance), which is comparable to the precision yielded by Type Ia supernovae and far better than the ``expanding photosphere method'' applied to Type II supernovae. Using SN 1987A to calibrate the Hubble diagrams we get Ho=55+/-12.
△ Less
Submitted 17 January, 2002;
originally announced January 2002.
-
The Distance to SN 1999em from the Expanding Photosphere Method
Authors:
M. Hamuy,
P. A. Pinto,
J. Maza,
N. B. Suntzeff,
M. M. Phillips,
R. G. Eastman,
R. C. Smith,
C. J. Corbally,
D. Burstein,
Y. Li,
V. Ivanov,
A. Moro-Martin,
L. G. Strolger,
R. E. de Souza,
S. dos Anjos,
E. M. Green,
T. E. Pickering,
L. Gonzalez,
R. Antezana,
M. Wischnjewsky,
G. Galaz,
M. Roth,
S. E. Persson,
W. L. Freedman,
R. A. Schommer
Abstract:
We present optical and IR spectroscopy of the first two months of evolution of the Type II SN 1999em. We combine these data with high-quality optical/IR photometry beginning only three days after shock breakout, in order to study the performance of the ``Expanding Photosphere Method'' (EPM) in the determination of distances. With this purpose we develop a technique to measure accurate photospher…
▽ More
We present optical and IR spectroscopy of the first two months of evolution of the Type II SN 1999em. We combine these data with high-quality optical/IR photometry beginning only three days after shock breakout, in order to study the performance of the ``Expanding Photosphere Method'' (EPM) in the determination of distances. With this purpose we develop a technique to measure accurate photospheric velocities by cross-correlating observed and model spectra. The application of this technique to SN 1999em shows that we can reach an average uncertainty of 11% in velocity from an individual spectrum. Our analysis shows that EPM is quite robust to the effects of dust. In particular, the distances derived from the VI filters change by only 7% when the adopted visual extinction in the host galaxy is varied by 0.45 mag. The superb time sampling of the BVIZJHK light-curves of SN 1999em permits us to study the internal consistency of EPM and test the dilution factors computed from atmosphere models for Type II plateau supernovae. We find that, in the first week since explosion, the EPM distances are up to 50% lower than the average, possibly due the presence of circumstellar material. Over the following 65 days, on the other hand, our tests lend strong credence to the atmosphere models, and confirm previous claims that EPM can produce consistent distances without having to craft specific models to each supernova. This is particularly true for the VI filters which yield distances with an internal consistency of 4%. From the whole set of BVIZJHK photometry, we obtain an average distance of 7.5+/-0.5 Mpc, where the quoted uncertainty (7%) is a conservative estimate of the internal precision of the method obtained from the analysis of the first 70 days of the supernova evolution.
△ Less
Submitted 30 April, 2001;
originally announced May 2001.
-
Theory of Extrasolar Giant Planet Transits
Authors:
W. B. Hubbard,
J. J. Fortney,
J. I. Lunine,
A. Burrows,
D. Sudarsky,
P. A. Pinto
Abstract:
We present a synthesis of physical effects influencing the observed lightcurve of an extrasolar giant planet (EGP) transiting its host star. The synthesis includes a treatment of Rayleigh scattering, cloud scattering, refraction, and molecular absorption of starlight in the EGP atmosphere. Of these effects, molecular absorption dominates in determining the transit-derived radius $R$. Using a gen…
▽ More
We present a synthesis of physical effects influencing the observed lightcurve of an extrasolar giant planet (EGP) transiting its host star. The synthesis includes a treatment of Rayleigh scattering, cloud scattering, refraction, and molecular absorption of starlight in the EGP atmosphere. Of these effects, molecular absorption dominates in determining the transit-derived radius $R$. Using a generic model for the atmosphere of EGP HD209458b, we perform a fit to the best available transit lightcurve data, and infer that this planet has a radius at a pressure of 1 bar, $R_1$, equal to 94430 km. We predict that $R$ will be a function of wavelength of observation, with a robust prediction of variations of $\pm 1$% at infrared wavelengths where H$_2$O opacity in the high EGP atmosphere dominates.
△ Less
Submitted 9 June, 2001; v1 submitted 2 January, 2001;
originally announced January 2001.
-
A Test for the Nature of the Type Ia Supernova Explosion Mechanism
Authors:
Philip A. Pinto,
Ronald G. Eastman,
Tamara Rogers
Abstract:
Currently popular models for Type Ia supernovae (SNe Ia) fall into two general classes. The first comprises explosions of nearly pure carbon/oxygen (C/O) white dwarfs at the Chandrasekhar limit which ignite near their centers. The second consists of lower-mass C/O cores which are ignited by the detonation of an accreted surface helium layer. Explosions of the latter type produce copious Fe, Co a…
▽ More
Currently popular models for Type Ia supernovae (SNe Ia) fall into two general classes. The first comprises explosions of nearly pure carbon/oxygen (C/O) white dwarfs at the Chandrasekhar limit which ignite near their centers. The second consists of lower-mass C/O cores which are ignited by the detonation of an accreted surface helium layer. Explosions of the latter type produce copious Fe, Co and Ni K-alpha emission from 56Ni and 56Co decay in the detonated surface layers, emission which is much weaker from Chandrasekhar-mass models. The presence of this emission provides a simple and unambiguous discriminant between these two models for SNe Ia. Both mechanisms may produce 0.1-0.6 solar masses of 56Ni, making them bright gamma-ray line emitters. The time to maximum brightness of 56Ni decay lines is distinctly shorter in the sub-Chandrasekhar mass class of model (approximately 15 days) than in the Chandrasekhar mass model (approximately 30 days), making gamma-ray line evolution another direct test of the explosion mechanism. It should just be possible to detect K-shell emission from a sub-Chandrasekhar explosion from SNe Ia as far away as the Virgo cluster with the XMM Observatory. A 1 to 2 square meter X-ray telescope such as the proposed Con-X Observatory could observe K-alpha emission from sub-Chandrasekhar mass SNe Ia in the Virgo cluster, providing not just a detection, but high-accuracy flux and kinematic information.
△ Less
Submitted 21 August, 2000;
originally announced August 2000.
-
The Type Ia Supernova Width-Luminosity Relation
Authors:
Philip A. Pinto,
Ronald G. Eastman
Abstract:
Observations of type Ia supernovae at high redshifts have become an important tool for studying the geometry of the universe. The relation between the duration of the peak phase of a type Ia supernova's lightcurve and its luminosity (broader is brighter) forms the cornerstone of this measurement, yet it is a purely empirical relation. In this paper we show that the relation is a natural conseque…
▽ More
Observations of type Ia supernovae at high redshifts have become an important tool for studying the geometry of the universe. The relation between the duration of the peak phase of a type Ia supernova's lightcurve and its luminosity (broader is brighter) forms the cornerstone of this measurement, yet it is a purely empirical relation. In this paper we show that the relation is a natural consequence of the radiation transport in type Ia supernovae and suggest constraints on the nature of the explosions which arise from our interpretation of the observed relation. The principle parameter underlying the relation is the mass of radioactive 56Ni produced in the explosion. The relation is shown to be relatively insensitive to most other parameters in Chandrasekhar-mass explosions. Cosmological results are thus unlikely to suffer from systematic effects stemming from evolution in the explosions' progenitors.
△ Less
Submitted 13 June, 2000;
originally announced June 2000.
-
On the Relation Between Peak Luminosity and Parent Population of Type Ia Supernovae: A New Tool for Probing the Ages of Distant Galaxies
Authors:
Valentin D. Ivanov,
Mario Hamuy,
Philip A. Pinto
Abstract:
We study the properties of Type Ia Supernovae (SNe Ia) as functions of the radial distance from their host galaxy centers. Using a sample of 62 SNe Ia with reliable luminosity, reddening, and decline rate determinations, we find no significant radial gradients of SNe Ia peak absolute magnitudes or decline rates in elliptical+S0 galaxies, suggesting that the diversity of SN properties is not rela…
▽ More
We study the properties of Type Ia Supernovae (SNe Ia) as functions of the radial distance from their host galaxy centers. Using a sample of 62 SNe Ia with reliable luminosity, reddening, and decline rate determinations, we find no significant radial gradients of SNe Ia peak absolute magnitudes or decline rates in elliptical+S0 galaxies, suggesting that the diversity of SN properties is not related to the metallicity of their progenitors. We do find that the range in brightness and light curve width of supernovae in spiral galaxies extends to brighter, broader values. These results are interpreted as support for an age, but not metallicity, related origin of the diversity in SNe Ia. If confirmed with a larger and more accurate sample of data, the age-luminosity relation would offer a new and powerful tool to probe the ages and age gradients of stellar populations in galaxies at redshift as high as $z\sim1-2$. The absence of significant radial gradients in the peak $\rm (B-V)_0$ and $\rm (V-I)_0$ colors of SNe Ia supports the redding correction method of Phillips et al (1999). We find no radial gradient in residuals from the SN Ia luminosity-width relation, suggesting that the relation is not affected by properties of the progenitor populations and supporting the reliability of cosmological results based upon the use of SNe Ia as distance indicators.
△ Less
Submitted 3 June, 2000;
originally announced June 2000.
-
A Search for Environmental Effects on Type Ia Supernovae
Authors:
Mario Hamuy,
S. C. Trager,
Philip A. Pinto,
M. M. Phillips,
R. A. Schommer,
Valentin Ivanov,
Nicholas B. Suntzeff
Abstract:
We use integrated colors and B and V absolute magnitudes of Type Ia supernova (SN) host galaxies in order to search for environmental effects on the SN optical properties. With the new sample of 44 SNe we confirm the conclusion by Hamuy et al. (1996a) that bright events occur preferentially in young stellar environments. We find also that the brightest SNe occur in the least luminous galaxies, a…
▽ More
We use integrated colors and B and V absolute magnitudes of Type Ia supernova (SN) host galaxies in order to search for environmental effects on the SN optical properties. With the new sample of 44 SNe we confirm the conclusion by Hamuy et al. (1996a) that bright events occur preferentially in young stellar environments. We find also that the brightest SNe occur in the least luminous galaxies, a possible indication that metal-poorer neighbourhoods produce the more luminous events. The interpretation of these results is made difficult, however, due to the fact that galaxies with younger stellar populations are also lower in luminosity. In an attempt to remove this ambiguity we use models for the line strengths in the absorption spectrum of five early-type galaxies, in order to estimate metallicities and ages of the SN host galaxies. With the addition of abundance estimates from nebular analysis of the emission spectra of three spiral galaxies, we find possible further evidence that luminous SNe are produced in metal-poor neighborhoods. Further spectroscopic observations of the SN host galaxies will be necessary to test these results and assist in disentangling the age/metallicity effects on Type Ia SNe.
△ Less
Submitted 9 May, 2000;
originally announced May 2000.
-
A New Algorithm for Supernova Neutrino Transport and Some Applications
Authors:
A. Burrows,
T. Young,
P. A. Pinto,
R. Eastman,
T. Thompson
Abstract:
We have developed an implicit, multi-group, time-dependent, spherical neutrino transport code based on the Feautrier variables, the tangent-ray method, and accelerated ${\bf Λ}$ iteration. The code achieves high angular resolution, is good to O($v/c$), is equivalent to a Boltzmann solver (without gravitational redshifts), and solves the transport equation at all optical depths with precision. In…
▽ More
We have developed an implicit, multi-group, time-dependent, spherical neutrino transport code based on the Feautrier variables, the tangent-ray method, and accelerated ${\bf Λ}$ iteration. The code achieves high angular resolution, is good to O($v/c$), is equivalent to a Boltzmann solver (without gravitational redshifts), and solves the transport equation at all optical depths with precision. In this paper, we present our formulation of the relevant numerics and microphysics and explore protoneutron star atmospheres for snapshot post-bounce models. Our major focus is on spectra, neutrino-matter heating rates, Eddington factors, angular distributions, and phase-space occupancies. In addition, we investigate the influence on neutrino spectra and heating of final-state electron blocking, stimulated absorption, velocity terms in the transport equation, neutrino-nucleon scattering asymmetry, and weak magnetism and recoil effects. Furthermore, we compare the emergent spectra and heating rates obtained using full transport with those obtained using representative flux-limited transport formulations to gauge their accuracy and viability. Finally, we derive useful formulae for the neutrino source strength due to nucleon-nucleon bremsstrahlung and determine bremsstrahlung's influence on the emergent $ν_μ$ and $ν_τ$ neutrino spectra.
△ Less
Submitted 27 March, 2000; v1 submitted 11 May, 1999;
originally announced May 1999.
-
Selection Effects, Biases, and Constraints in the Calan/Tololo Supernova Survey
Authors:
Mario Hamuy,
Philip A. Pinto
Abstract:
We use Monte Carlo simulations of the Calan/Tololo photographic supernova survey to show that a simple model of the survey's selection effects accounts for the observed distributions of recession velocity, apparent magnitude, angular offset, and projected radial distance between the supernova and the host galaxy nucleus for this sample of Type Ia supernovae (SNe Ia). The model includes biases du…
▽ More
We use Monte Carlo simulations of the Calan/Tololo photographic supernova survey to show that a simple model of the survey's selection effects accounts for the observed distributions of recession velocity, apparent magnitude, angular offset, and projected radial distance between the supernova and the host galaxy nucleus for this sample of Type Ia supernovae (SNe Ia). The model includes biases due to the flux-limited nature of the survey, the different light curve morphologies displayed by different SNe Ia, and the difficulty of finding events projected near the central regions of the host galaxies. From these simulations we estimate the bias in the zero-point and slope of the absolute magnitude-decline rate relation used in SNe Ia distance measurements. For an assumed intrinsic scatter of 0.15 mag about this relation, these selection effects decrease the zero-point by 0.04 mag. The slope of the relation is not significantly biased. We conclude that despite selection effects in the survey, the shape and zero-point of the relation determined from the Calan/Tololo sample are quite reliable. We estimate the degree of incompleteness of the survey as a function of decline rate and estimate a corrected luminosity function for SNe Ia in which the frequency of SNe appears to increase with decline rate (the fainter SNe are more common). Finally, we compute the integrated detection efficiency of the survey in order to infer the rate of SNe Ia from the 31 events found. For a value of Ho=65 km/sec/Mpc we obtain a SN Ia rate of 0.21(+0.30)(-0.13) SNu. This is in good agreement with the value 0.16+/-0.05 SNu recently determined by Capellaro et al. (1997).
△ Less
Submitted 3 December, 1998;
originally announced December 1998.
-
Iron Opacity and the Pulsar of Supernova 1987A
Authors:
C. L. Fryer,
S. A. Colgate,
P. A. Pinto
Abstract:
Neutron stars formed in Type II supernovae are likely to be initially obscured by late-time fallback. Although much of the late-time fallback is quickly accreted via neutrino cooling, some material remains on the neutron star, forming an atmosphere which slowly accretes through photon emission. In this paper, we derive structure equations of the fallback atmosphere and present results of one-dim…
▽ More
Neutron stars formed in Type II supernovae are likely to be initially obscured by late-time fallback. Although much of the late-time fallback is quickly accreted via neutrino cooling, some material remains on the neutron star, forming an atmosphere which slowly accretes through photon emission. In this paper, we derive structure equations of the fallback atmosphere and present results of one-dimensional simulations of that fallback. The atmosphere remaining after neutrino cooling becomes unimportant (less than the Compton Eddington limit) is only a fraction of the total mass accreted (10^-8 of the accreted mass or 10^-9 solar masses.) Recombined iron dominates the opacity in the outer regions leading to an opacity 1000-10,000 times higher than that of electron scattering alone. The resultant photon emission of the remnant atmosphere is limited to 1/1000th the Compton Eddington Luminosity. The late-time evolution of this system leads to the formation of a photon-driven wind from the accretion of the inner portion of the atmosphere, leaving, for most cases, a bare neutron star on timescales shorter than a year. The degenerate remnant of 1987a may not be a black hole. Instead, the fallback material may have already accreted or blown off in the accretion-driven wind. If the neutron star has either a low magnetic field or a low rotational spin frequency, we would not expect to see the neutron star remnant of 1987a.
△ Less
Submitted 27 August, 1998;
originally announced August 1998.
-
Infrared and Optical Spectroscopy of Type Ia Supernovae in the Nebular Phase
Authors:
E. J. C. Bowers,
W. P. S. Meikle,
T. R. Geballe,
N. A. Walton,
P. A. Pinto,
V. S. Dhillon,
S. B. Howell,
M. K. Harrop-Allin
Abstract:
We present near-infrared (NIR) spectra for Type Ia supernovae at epochs of 13 to 338 days after maximum blue light. Some contemporary optical spectra are also shown. All the NIR spectra exhibit considerable structure throughout the J-, H- and K-bands. In particular they exhibit a flux `deficit' in the J-band which persists as late as 175 days. This is responsible for the well-known red J-H colou…
▽ More
We present near-infrared (NIR) spectra for Type Ia supernovae at epochs of 13 to 338 days after maximum blue light. Some contemporary optical spectra are also shown. All the NIR spectra exhibit considerable structure throughout the J-, H- and K-bands. In particular they exhibit a flux `deficit' in the J-band which persists as late as 175 days. This is responsible for the well-known red J-H colour. To identify the emission features and test the $^{56}$Ni hypothesis for the explosion and subsequent light curve, we compare the NIR and optical nebular-phase data with a simple non-LTE nebular spectral model. We find that many of the spectral features are due to iron-group elements and that the J-band deficit is due to a lack of emission lines from species which dominate the rest of the IR/optical spectrum. Nevertheless, some emission is unaccounted for, possibly due to inaccuracies in the cobalt atomic data. For some supernovae, blueshifts of 1000--3000 km/s are seen in infrared and optical features at 3 months. We suggest this is due to clumping in the ejecta. The evolution of the cobalt/iron mass ratio indicates that $^{56}$Co-decay dominates the abundances of these elements. The absolute masses of iron-group elements which we derive support the basic thermonuclear explosion scenario for Type Ia supernovae. A core-collapse origin is less consistent with our data.
△ Less
Submitted 10 July, 1997;
originally announced July 1997.
-
The Physics of Type Ia Supernova Lightcurves: I. Opacity and Diffusion
Authors:
Philip A. Pinto,
Ronald G. Eastman
Abstract:
We examine the basic physics of type Ia supernova (SNe Ia) light curves with a view toward interpreting the relations between peak luminosity, peak width, and late-time slope in terms of the properties of the underlying explosion models. We develop an analytic solution of the radiation transport problem in the dynamic diffusion regime and show that it reproduces bolometric light curves produced…
▽ More
We examine the basic physics of type Ia supernova (SNe Ia) light curves with a view toward interpreting the relations between peak luminosity, peak width, and late-time slope in terms of the properties of the underlying explosion models. We develop an analytic solution of the radiation transport problem in the dynamic diffusion regime and show that it reproduces bolometric light curves produced by more detailed calculations under the assumption of a constant extinction coefficient. This model is used to derive the thermal conditions in the interior of SNe Ia, to address the issue of time dependence in the modeling of SNe Ia light curves and spectra and to show that these are intrinsically time-dependent phenomena. The analytic model is then used to study the sensitivity of light curves to various properties of supernova explosions. We show that the dominant opacity arises from line transitions and discuss the nature of line opacities in expanding media and the appropriateness of various mean opacities used in light curve calculations.
△ Less
Submitted 22 November, 1996;
originally announced November 1996.
