We consider the cosmological models based on Palatini f(R)-theory for the function f(R)=aR-2bR^2-3c/R, which, when only dust visible matter is considered, is called dune cosmology in view of the shape of the function f(R(a)) (being a the scale factor). We discuss about the meaning of solving the model, and interpret it according to Ehlers-Pirani-Schild framework as defining a Weyl geometry on spacetime. Accordingly, we extend the definitions of luminosity distance, proper distance, and redshift to Weyl geometries and fit the values of parameters to SNIa data. Since the theoretical prediction is model-dependent, we argue that the it is affected by an extra choice, namely a model for atomic clocks, which, in principle, produces observable effects. To the best of our knowledge, these effects have not being considered in the literature before.

Cluster weak lensing is a sensitive probe of cosmology, particularly the amplitude of matter clustering $\sigma_8$ and matter density parameter $\Omega_m$. The main nuisance parameter in a cluster weak lensing cosmological analysis is the scatter between the true halo mass and the relevant cluster observable, denoted $\sigma_{\ln Mc}$. We show that combining the cluster weak lensing observable $\Delta \Sigma$ with the projected cluster-galaxy cross-correlation function $w_{p,cg}$ and galaxy auto-correlation function $w_{p,gg}$ can break the degeneracy between $\sigma_8$ and $\sigma_{\ln Mc}$ to achieve tight, percent-level constraints on $\sigma_8$. Using a grid of cosmological N-body simulations, we compute derivatives of $\Delta \Sigma$, $w_{p,cg}$, and $w_{p,gg}$ with respect to $\sigma_8$, $\Omega_m$, $\sigma_{\ln Mc}$ and halo occupation distribution (HOD) parameters describing the galaxy population. We also compute covariance matrices motivated by the properties of the Dark Energy Suvery (DES) cluster and weak lensing survey and the BOSS CMASS galaxy redshift survey. For our fiducial scenario combining $\Delta \Sigma$, $w_{p,cg}$, and $w_{p,gg}$ measured over $0.3-30.0 \; h^{-1} \; \mathrm{Mpc}$, for clusters at $z=0.35-0.55$ above a mass threshold $M_c\approx 2\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$, we forecast a $1.4\%$ constraint on $\sigma_8$ while marginalizing over $\sigma_{\ln Mc}$ and all HOD parameters. Reducing the mass threshold to $1\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$ and adding a $z=0.15-0.35$ redshift bin sharpens this constraint to $0.8\%$. The small scale $(r_p < 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``mass function'' and large scale $(r_p > 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``halo-mass cross-correlation'' regimes of $\Delta \Sigma$ have comparable constraining power, allowing internal consistency tests from such an analysis.

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 .

(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\'{o}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$\sigma$ 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 $\delta<+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.

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

The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy bias. Using a grid of cosmological N-body simulations, we construct a Taylor-expansion emulator that predicts the galaxy autocorrelation $\xi_{\text{gg}}(r)$ and galaxy-matter cross-correlation $\xi_{\text{gm}}(r)$ as a function of $\sigma_8$, $\Omega_m$, and halo occupation distribution (HOD) parameters, which are allowed to vary with large scale environment to represent possible effects of galaxy assembly bias. We present forecasts for a fiducial case that corresponds to BOSS LOWZ galaxy clustering and SDSS-depth weak lensing (effective source density $\sim 0.3$ arcmin$^{-2}$). Using tangential shear and projected correlation function measurements over $0.5 \leq r_p \leq 30 \, h^{-1}$ Mpc yields a 1.8% constraint on the parameter combination $\sigma_8\Omega_m^{0.58}$, a factor of two better than a constraint that excludes non-linear scales ($r_p > 2 \, h^{-1}$ Mpc, $4 \, h^{-1}$ Mpc for $\gamma_t,w_p$). Much of this improvement comes from the non-linear clustering information, which breaks degeneracies among HOD parameters that would otherwise degrade the inference of matter clustering from GGL. Increasing the effective source density to $3$ arcmin$^{-2}$ sharpens the constraint on $\sigma_8\Omega_m^{0.58}$ by a further factor of two. With robust modeling into the non-linear regime, low-redshift measurements of matter clustering at the 1-percent level with clustering+GGL alone are well within reach of current data sets such as those provided by the Dark Energy Survey.

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

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

The baryon acoustic oscillation (BAO) feature in the clustering of matter in the universe serves as a robust standard ruler and hence can be used to map the expansion history of the universe. We use high force resolution simulations to analyze the effects of galaxy bias on the measurements of the BAO signal. We apply a variety of Halo Occupation Distributions (HODs) and produce biased mass tracers to mimic different galaxy populations. We investigate whether galaxy bias changes the non-linear shifts on the acoustic scale relative to the underlying dark matter distribution presented by Seo et al (2009). For the less biased HOD models (b < 3), we do not detect any shift in the acoustic scale relative to the no-bias case, typically 0.10% \pm 0.10%. However, the most biased HOD models (b > 3) show a shift at moderate significance (0.79% \pm 0.31% for the most extreme case). We test the one-step reconstruction technique introduced by Eisenstein et al. (2007) in the case of realistic galaxy bias and shot noise. The reconstruction scheme increases the correlation between the initial and final (z = 1) density fields achieving an equivalent level of correlation at nearly twice the wavenumber after reconstruction. Reconstruction reduces the shifts and errors on the shifts. We find that after reconstruction the shifts from the galaxy cases and the dark matter case are consistent with each other and with no shift. The 1-sigma systematic errors on the distance measurements inferred from our BAO measurements with various HODs after reconstruction are about 0.07% - 0.15%.

We measure shifts of the acoustic scale due to nonlinear growth and redshift distortions to a high precision using a very large volume of high-force-resolution simulations. We compare results from various sets of simulations that differ in their force, volume, and mass resolution. We find a consistency within 1.5-sigma for shift values from different simulations and derive shift alpha(z) -1 = (0.300\pm 0.015)% [D(z)/D(0)]^{2} using our fiducial set. We find a strong correlation with a non-unity slope between shifts in real space and in redshift space and a weak correlation between the initial redshift and low redshift. Density-field reconstruction not only removes the mean shifts and reduces errors on the mean, but also tightens the correlations: after reconstruction, we recover a slope of near unity for the correlation between the real and redshift space and restore a strong correlation between the low and the initial redshifts. We derive propagators and mode-coupling terms from our N-body simulations and compared with Zeldovich approximation and the shifts measured from the chi^2 fitting, respectively. We interpret the propagator and the mode-coupling term of a nonlinear density field in the context of an average and a dispersion of its complex Fourier coefficients relative to those of the linear density field; from these two terms, we derive a signal-to-noise ratio of the acoustic peak measurement. We attempt to improve our reconstruction method by implementing 2LPT and iterative operations: we obtain little improvement. The Fisher matrix estimates of uncertainty in the acoustic scale is tested using 5000 (Gpc/h)^3 of cosmological PM simulations from Takahashi et al. (2009). (abridged)

We introduce a new statistic omega_l for measuring and analyzing large-scale structure and particularly the baryon acoustic oscillations. omega_l is a band-filtered, configuration space statistic that is easily implemented and has advantages over the traditional power spectrum and correlation function estimators. Unlike these estimators, omega_l can localize most of the acoustic information into a single dip at the acoustic scale while also avoiding sensitivity to the poorly constrained large scale power (i.e., the integral constraint) through the use of a localized and compensated filter. It is also sensitive to anisotropic clustering through pair counting and does not require any binning. We measure the shift in the acoustic peak due to nonlinear effects using the monopole omega_0 derived from subsampled dark matter catalogues as well as from mock galaxy catalogues created via halo occupation distribution (HOD) modeling. All of these are drawn from 44 realizations of 1024^3 particle dark matter simulations in a 1h^{-1}Gpc box at z=1. We compare these shifts with those obtained from the power spectrum and conclude that the results agree. This indicates that any distance measurements obtained from omega_0 and P(k) will be consistent with each other. We also show that it is possible to extract the same amount of acoustic information using either omega_0 or P(k) from equal volume surveys.

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.

Since type Ia Supernovae (SNe) explode in galaxies, they can, in principle, be used as the same tracer of the large-scale structure as their hosts to measure baryon acoustic oscillations (BAOs). To realize this, one must obtain a dense integrated sampling of SNe over a large fraction of the sky, which may only be achievable photometrically with future projects such as the Large Synoptic Survey Telescope. The advantage of SN BAOs is that SNe have more uniform luminosities and more accurate photometric redshifts than galaxies, but the disadvantage is that they are transitory and hard to obtain in large number at high redshift. We find that a half-sky photometric SN survey to redshift z = 0.8 is able to measure the baryon signature in the SN spatial power spectrum. Although dark energy constraints from SN BAOs are weak, they can significantly improve the results from SN luminosity distances of the same data, and the combination of the two is no longer sensitive to cosmic microwave background priors.

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.

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.

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.

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.

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

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.

We discuss the issue of toy model building for the dark energy component of the universe. Specifically, we consider two generic toy models recently proposed as alternatives to quintessence models, known as Cardassian expansion and the Chaplygin gas. We show that the former is enteriely equivalent to a class of quintessence models. We determine the observational constraints on the latter, coming from recent supernovae results and from the shape of the matter power spectrum. As expected, these restrict the model to a behaviour that closely matches that of a standard cosmological constant $\Lambda$.

The median statistic has recently been discussed by Gott \textit{et al.} as a more reliable alternative to the standard $\chi^2$ likelihood analysis, in the sense of requiring fewer assumptions about the data and being almost as constraining. We apply this statistic to the currently available combined dataset of 92 distant type Ia supernovae, and also to a mock SNAP-class dataset. We find that the performances of the modified median and $\chi^2$ statistics are comparable, particularly in the latter case. We further extend the work of Gott \textit{et al.} by modifying the median statistic to account for the number and size of sequences of consecutive points above or below the median. We also comment on how the performance of the statistic depends on the choice of free parameters that one is estimating.

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.

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.

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.

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.

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.

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.

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.

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.

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

We generate theoretical albedo and reflection spectra for a full range of extrasolar giant planet (EGP) models, from Jovian to 51-Pegasi class objects. Our albedo modeling utilizes the latest atomic and molecular cross sections, a Mie theory treatment of extinction by condensates, a variety of particle size distributions, and an extension of the Feautrier radiative transfer method which allows for a general treatment of the scattering phase function. We find that due to qualitative similarities in the compositions and spectra of objects within each of five broad effective temperature ranges, it is natural to establish five representative EGP albedo classes: a ``Jovian'' class (T$_{\rm eff} \lesssim 150$ K; Class I) with tropospheric ammonia clouds, a ``water cloud'' class (T$_{\rm eff} \sim 250$ K; Class II) primarily affected by condensed H$_2$O, a ``clear'' class (T$_{\rm eff} \gtrsim 350$ K; Class III) which lacks clouds, and two high-temperature classes: Class IV (900 K $\lesssim$ T$_{\rm{eff}}$ $\lesssim$ 1500 K) for which alkali metal absorption predominates, and Class V (T$_{\rm{eff}}$ $\gtrsim$ 1500 K and/or low surface gravity ($\lesssim$ 10$^3$ cm s$^{-2}$)) for which a high silicate layer shields a significant fraction of the incident radiation from alkali metal and molecular absorption. The resonance lines of sodium and potassium are expected to be salient features in the reflection spectra of Class III, IV, and V objects. We derive Bond albedos and effective temperatures for the full set of known EGPs and explore the possible effects of non-equilibrium condensed products of photolysis above or within principal cloud decks. As in Jupiter, such species can lower the UV/blue albedo substantially, even if present in relatively small mixing ratios.

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

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.

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.

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.