Polarization spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multidimensional supernova explosion models. The approach utilizes 'virtual-packets' that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is not only vital for calculating synthetic spectropolarimetry (since the degree of polarization is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealized test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarization from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarization spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

Selecting superluminous supernovae in faint galaxies from the first year of the Pan-STARRS1 Medium Deep Survey

The Pan-STARRS1 (PS1) survey has obtained imaging in five bands (grizyP1) over 10 Medium Deep Survey (MDS) fields covering a total of 70 square degrees. This paper describes the search for apparently hostless supernovae (SNe) within the first year of PS1 MDS data with an aim of discovering superluminous supernovae (SLSNe). A total of 249 hostless transients were discovered down to a limiting magnitude of MAB ∼ 23.5, of which 76 were classified as Type Ia supernovae (SNe Ia). There were 57 SNe with complete light curves that are likely core-collapse SNe (CCSNe) or type Ic SLSNe and 12 of these have had spectra taken. Of these 12 hostless, non-Type Ia SNe, 7 were SLSNe of type Ic at redshifts between 0.5 and 1.4. This illustrates that the discovery rate of type Ic SLSNe can be maximized by concentrating on hostless transients and removing normal SNe Ia. We present data for two possible SLSNe; PS1-10pm (z = 1.206) and PS1-10ahf (z = 1.1), and estimate the rate of type Ic SLSNe to be between 3+3−2×10−53+3−2×10−53+3−2×10−5 and 8+2−1×10−58+2−1×10−58+2−1×10−5 that of the CCSN rate within 0.3 ≤ z ≤ 1.4 by applying a Monte Carlo technique. The rate of slowly evolving, type Ic SLSNe (such as SN2007bi) is estimated as a factor of 10 lower than this range.

Systematic Uncertainties Associated with the Cosmological Analysis of the First PAN-STARRS1 Type Ia Supernova Sample

We probe the systematic uncertainties from the 113 Type Ia supernovae (SN Ia) in the Pan-STARRS1 (PS1) sample along with 197 SN Ia from a combination of low-redshift surveys. The companion paper by Rest et al. describes the photometric measurements and cosmological inferences from the PS1 sample. The largest systematic uncertainty stems from the photometric calibration of the PS1 and low-z samples. We increase the sample of observed Calspec standards from 7 to 10 used to define the PS1 calibration system. The PS1 and SDSS-II calibration systems are compared and discrepancies up to ∼0.02 mag are recovered. We find uncertainties in the proper way to treat intrinsic colors and reddening produce differences in the recovered value of w up to 3%. We estimate masses of host galaxies of PS1 supernovae and detect an insignificant difference in distance residuals of the full sample of 0.037 ± 0.031 mag for host galaxies with high and low masses. Assuming flatness and including systematic uncertainties in our analysis of only SNe measurements, we find w = -1.120^+0.360_-0.206(Stat)^+0.269_-0.291(Sys). With additional constraints from Baryon acoustic oscillation, cosmic microwave background (CMB) (Planck) and H₀ measurements, we find w = -1.166^+0.072_-0.069 and Ω_m = 0.280^+0.013_-0.012 (statistical and systematic errors added in quadrature). The significance of the inconsistency with w = -1 depends on whether we use Planck or Wilkinson Microwave Anisotropy Probe measurements of the CMB: w_{BAO+H0+SN+WMAP} = -1.124^+0.083_-0.065.

Transition redshift: new kinematic constraints from supernovae

The transition redshift (deceleration/acceleration) is discussed by expanding the deceleration parameter to first order around its present value. A detailed study is carried out by considering two different parametrizations, q = q(0) + q(1)z and q = q(0) + q(1)z(1 + z)(-1), and the associated free parameters (q(0), q(1)) are constrained by three different supernovae (SNe) samples. A previous analysis by Riess et al. using the first expansion is slightly improved and confirmed in light of their recent data (Gold07 sample). However, by fitting the model with the Supernova Legacy Survey (SNLS) type Ia sample, we find that the best fit to the redshift transition is z(t) = 0.61, instead of z(t) = 0.46 as derived by the High-z Supernovae Search (HZSNS) team. This result based in the SNLS sample is also in good agreement with the sample of Davis et al., z(t) = 0.60(-0.11)(+0.28) (1 sigma). Such results are in line with some independent analyses and accommodate more easily the concordance flat model (Lambda CDM). For both parametrizations, the three SNe Ia samples considered favour recent acceleration and past deceleration with a high degree of statistical confidence level. All the kinematic results presented here depend neither on the validity of general relativity nor on the matter-energy contents of the Universe.

Type Ia Supernova Cosmology : Quantitative Spectral Analysis

Type Ia supernovae have been successfully used as standardized candles to study the expansion history of the Universe. In the past few years, these studies led to the exciting result of an accelerated expansion caused by the repelling action of some sort of dark energy. This result has been confirmed by measurements of cosmic microwave background radiation, the large-scale structure, and the dynamics of galaxy clusters. The combination of all these experiments points to a “concordance model” of the Universe with flat large-scale geometry and a dominant component of dark energy. However, there are several points related to supernova measurements which need careful analysis in order to doubtlessly establish the validity of the concordance model. As the amount and quality of data increases, the need of controlling possible systematic effects which may bias the results becomes crucial. Also important is the improvement of our knowledge of the physics of supernovae events to assure and possibly refine their calibration as standardized candle. This thesis addresses some of those issues through the quantitative analysis of supernova spectra. The stress is put on a careful treatment of the data and on the definition of spectral measurement methods. The comparison of measurements for a large set of spectra from nearby supernovae is used to study the homogeneity and to search for spectral parameters which may further refine the calibration of the standardized candle. One such parameter is found to reduce the dispersion in the distance estimation of a sample of supernovae to below 6%, a precision which is comparable with the current lightcurve-based calibration, and is obtained in an independent manner. Finally, the comparison of spectral measurements from nearby and distant objects is used to test the possibility of evolution with cosmic time of the intrinsic brightness of type Ia supernovae.

Cosmic Explosions: Observations Of Infant Hydrogen-Free Supernovae Towards An Understanding Of Their Parent Systems

Radiation in the first days of supernova explosions contains rich information about physical properties of the exploding stars. In the past three years, I used the intermediate Palomar Transient Factory to conduct one-day cadence surveys, in order to systematically search for infant supernovae. I show that the one-day cadences in these surveys were strictly controlled, that the realtime image subtraction pipeline managed to deliver transient candidates within ten minutes of images being taken, and that we were able to undertake follow-up observations with a variety of telescopes within hours of transients being discovered. So far iPTF has discovered over a hundred supernovae within a few days of explosions, forty-nine of which were spectroscopically classified within twenty-four hours of discovery.

Our observations of infant Type Ia supernovae provide evidence for both the single-degenerate and double-degenerate progenitor channels. On the one hand, a low-velocity Type Ia supernova iPTF14atg revealed a strong ultraviolet pulse within four days of its explosion. I show that the pulse is consistent with the expected emission produced by collision between the supernova ejecta and a companion star, providing direct evidence for the single degenerate channel. By comparing the distinct early-phase light curves of iPTF14atg to an otherwise similar event iPTF14dpk, I show that the viewing angle dependence of the supernova-companion collision signature is probably responsible to the difference of the early light curves. I also show evidence for a dark period between the supernova explosion and the first light of the radioactively-powered light curve. On the other hand, a peculiar Type Ia supernova iPTF13asv revealed strong near-UV emission and absence of iron in the spectra within the first two weeks of explosion, suggesting a stratified ejecta structure with iron group elements confined to the slow-moving part of the ejecta. With its total ejecta mass estimated to exceed the Chandrasekhar limit, I show that the stratification and large mass of the ejecta favor the double-degenerate channel.

In a separate approach, iPTF found the first progenitor system of a Type Ib supernova iPTF13bvn in the pre-explosion HST archival mages. Independently, I used the early-phase optical observations of this supernova to constrain its progenitor radius to be no larger than several solar radii. I also used its early radio detections to derive a mass loss rate of 3e-5 solar mass per year for the progenitor right before the supernova explosion. These constraints on the physical properties of the iPTF13bvn progenitor provide a comprehensive data set to test Type Ib supernova theories. A recent HST revisit to the iPTF13bvn site two years after the supernova explosion has confirmed the progenitor system.

Moving forward, the next frontier in this area is to extend these single-object analyses to a large sample of infant supernovae. The upcoming Zwicky Transient Facility with its fast survey speed, which is expected to find one infant supernova every night, is well positioned to carry out this task.

Acceleration of the Cosmological Expansion as an Effect of Inhomogeneities

The cosmological observations of light from type Ia supernovae, the cosmic microwave background and the galaxy distribution seem to indicate that the expansion of the universe has accelerated during the latter half of its age. Within standard cosmology, this is ascribed to dark energy, a uniform fluid with large negative pressure that gives rise to repulsive gravity but also entails serious theoretical problems. Understanding the physical origin of the perceived accelerated expansion has been described as one of the greatest challenges in theoretical physics today. In this thesis, we discuss the possibility that, instead of dark energy, the acceleration would be caused by an effect of the nonlinear structure formation on light, ignored in the standard cosmology. A physical interpretation of the effect goes as follows: due to the clustering of the initially smooth matter with time as filaments of opaque galaxies, the regions where the detectable light travels get emptier and emptier relative to the average. As the developing voids begin to expand the faster the lower their matter density becomes, the expansion can then accelerate along our line of sight without local acceleration, potentially obviating the need for the mysterious dark energy. In addition to offering a natural physical interpretation to the acceleration, we have further shown that an inhomogeneous model is able to match the main cosmological observations without dark energy, resulting in a concordant picture of the universe with 90% dark matter, 10% baryonic matter and 15 billion years as the age of the universe. The model also provides a smart solution to the coincidence problem: if induced by the voids, the onset of the perceived acceleration naturally coincides with the formation of the voids. Additional future tests include quantitative predictions for angular deviations and a theoretical derivation of the model to reduce the required phenomenology. A spin-off of the research is a physical classification of the cosmic inhomogeneities according to how they could induce accelerated expansion along our line of sight. We have identified three physically distinct mechanisms: global acceleration due to spatial variations in the expansion rate, faster local expansion rate due to a large local void and biased light propagation through voids that expand faster than the average. A general conclusion is that the physical properties crucial to account for the perceived acceleration are the growth of the inhomogeneities and the inhomogeneities in the expansion rate. The existence of these properties in the real universe is supported by both observational data and theoretical calculations. However, better data and more sophisticated theoretical models are required to vindicate or disprove the conjecture that the inhomogeneities are responsible for the acceleration.

Measuring the Cosmic Microwave Background : Topics along the analysis pipeline

The first quarter of the 20th century witnessed a rebirth of cosmology, study of our Universe, as a field of scientific research with testable theoretical predictions. The amount of available cosmological data grew slowly from a few galaxy redshift measurements, rotation curves and local light element abundances into the first detection of the cos- mic microwave background (CMB) in 1965. By the turn of the century the amount of data exploded incorporating fields of new, exciting cosmological observables such as lensing, Lyman alpha forests, type Ia supernovae, baryon acoustic oscillations and Sunyaev-Zeldovich regions to name a few. -- CMB, the ubiquitous afterglow of the Big Bang, carries with it a wealth of cosmological information. Unfortunately, that information, delicate intensity variations, turned out hard to extract from the overall temperature. Since the first detection, it took nearly 30 years before first evidence of fluctuations on the microwave background were presented. At present, high precision cosmology is solidly based on precise measurements of the CMB anisotropy making it possible to pinpoint cosmological parameters to one-in-a-hundred level precision. The progress has made it possible to build and test models of the Universe that differ in the way the cosmos evolved some fraction of the first second since the Big Bang. -- This thesis is concerned with the high precision CMB observations. It presents three selected topics along a CMB experiment analysis pipeline. Map-making and residual noise estimation are studied using an approach called destriping. The studied approximate methods are invaluable for the large datasets of any modern CMB experiment and will undoubtedly become even more so when the next generation of experiments reach the operational stage. -- We begin with a brief overview of cosmological observations and describe the general relativistic perturbation theory. Next we discuss the map-making problem of a CMB experiment and the characterization of residual noise present in the maps. In the end, the use of modern cosmological data is presented in the study of an extended cosmological model, the correlated isocurvature fluctuations. Current available data is shown to indicate that future experiments are certainly needed to provide more information on these extra degrees of freedom. Any solid evidence of the isocurvature modes would have a considerable impact due to their power in model selection.

Strongly magnetized cold degenerate electron gas: Mass-radius relation of the magnetized white dwarf

We consider a relativistic, degenerate electron gas at zero temperature under the influence of a strong, uniform, static magnetic field, neglecting any form of interactions. Since the density of states for the electrons changes due to the presence of the magnetic field (which gives rise to Landau quantization), the corresponding equation of state also gets modified. In order to investigate the effect of very strong magnetic field, we focus only on systems in which a maximum of either one, two, or three Landau level(s) is/are occupied. This is important since, if a very large number of Landau levels are filled, it implies a very low magnetic field strength which yields back Chandrasekhar's celebrated nonmagnetic results. The maximum number of occupied Landau levels is fixed by the correct choice of two parameters, namely, the magnetic field strength and the maximum Fermi energy of the system. We study the equations of state of these one-level, two-level, and three-level systems and compare them by taking three different maximum Fermi energies. We also find the effect of the strong magnetic field on the mass-radius relation of the underlying star composed of the gas stated above. We obtain an exciting result that it is possible to have an electron-degenerate static star, namely, magnetized white dwarfs, with a mass significantly greater than the Chandrasekhar limit in the range 2.3-2.6M(circle dot), provided it has an appropriate magnetic field strength and central density. In fact, recent observations of peculiar type Ia supernovae-SN 2006gz, SN 2007if, SN 2009dc, SN 2003fg-seem to suggest super-Chandrasekhar-mass white dwarfs with masses up to 2.4-2.8M(circle dot) as their most likely progenitors. Interestingly, our results seem to lie within these observational limits.

Reconstructing the properties of dark energy using standard sirens

Future space-based gravity wave (GW) experiments such as the Big Bang Observatory (BBO), with their excellent projected, one sigma angular resolution, will measure the luminosity distance to a large number of GW sources to high precision, and the redshift of the single galaxies in the narrow solid angles towards the sources will provide the redshifts of the gravity wave sources. One sigma BBO beams contain the actual source in only 68% of the cases; the beams that do not contain the source may contain a spurious single galaxy, leading to misidentification. To increase the probability of the source falling within the beam, larger beams have to be considered, decreasing the chances of finding single galaxies in the beams. Saini et al. T.D. Saini, S.K. Sethi, and V. Sahni, Phys. Rev. D 81, 103009 (2010)] argued, largely analytically, that identifying even a small number of GW source galaxies furnishes a rough distance-redshift relation, which could be used to further resolve sources that have multiple objects in the angular beam. In this work we further develop this idea by introducing a self-calibrating iterative scheme which works in conjunction with Monte Carlo simulations to determine the luminosity distance to GW sources with progressively greater accuracy. This iterative scheme allows one to determine the equation of state of dark energy to within an accuracy of a few percent for a gravity wave experiment possessing a beam width an order of magnitude larger than BBO (and therefore having a far poorer angular resolution). This is achieved with no prior information about the nature of dark energy from other data sets such as type Ia supernovae, baryon acoustic oscillations, cosmic microwave background, etc. DOI:10.1103/PhysRevD.87.083001

NEW MASS LIMIT OF WHITE DWARFS

Is the Chandrasekhar mass limit for white dwarfs (WDs) set in stone? Not anymore, recent observations of over-luminous, peculiar type Ia supernovae can be explained if significantly super-Chandrasekhar WDs exist as their progenitors, thus barring them to be used as cosmic distance indicators. However, there is no estimate of a mass limit for these super-Chandrasekhar WD candidates yet. Can they be arbitrarily large? In fact, the answer is no! We arrive at this revelation by exploiting the flux freezing theorem in observed, accreting, magnetized WDs, which brings in Landau quantization of the underlying electron degenerate gas. This essay presents the calculations which pave the way for the ultimate (significantly super-Chandrasekhar) mass limit of WDs, heralding a paradigm shift 80 years after Chandrasekhar's discovery.

Maximum mass of stable magnetized highly super-Chandrasekhar white dwarfs: stable solutions with varying magnetic fields

We address the issue of stability of recently proposed significantly super-Chandrasekhar white dwarfs. We present stable solutions of magnetostatic equilibrium models for super-Chandrasekhar white dwarfs pertaining to various magnetic field profiles. This has been obtained by self-consistently including the effects of the magnetic pressure gradient and total magnetic density in a general relativistic framework. We estimate that the maximum stable mass of magnetized white dwarfs could be more than 3 solar mass. This is very useful to explain peculiar, overluminous type Ia supernovae which do not conform to the traditional Chandrasekhar mass-limit.

GRMHD formulation of highly super-Chandrasekhar magnetized white dwarfs: stable configurations of non-spherical white dwarfs

The topic of magnetized super-Chandrasekhar white dwarfs is in the limelight, particularly in the last few years, since our proposal of their existence. By full-scale general relativistic magnetohydrodynamic (GRMHD) numerical analysis, we confirm in this work the existence of stable, highly magnetized, significantly super-Chandrasekhar white dwarfs with mass more than 3 solar mass. While a poloidal field geometry renders the white dwarfs oblate, a toroidal field makes them prolate retaining an overall quasi-spherical shape, as speculated in our earlier work. These white dwarfs are expected to serve as the progenitors of over-luminous type Ia supernovae.

GRMHD formulation of highly super-Chandrasekhar rotating magnetized white dwarfs: stable configurations of non-spherical white dwarfs

Here we extend the exploration of significantly super-Chandrasekhar magnetized white dwarfs by numerically computing axisymmetric stationary equilibria of differentially rotating magnetized polytropic compact stars in general relativity (GR), within the ideal magnetohydrodynamic regime. We use a general relativistic magnetohydrodynamic (GRMHD) framework that describes rotating and magnetized axisymmetric white dwarfs, choosing appropriate rotation laws and magnetic field profiles (toroidal and poloidal). The numerical procedure for finding solutions in this framework uses the 3 + 1 formalism of numerical relativity, implemented in the open source XNS code. We construct equilibrium sequences by varying different physical quantities in turn, and highlight the plausible existence of super-Chandrasekhar white dwarfs, with masses in the range of 2-3 solar mass, with central (deep interior) magnetic fields of the order of 10(14) G and differential rotation with surface time periods of about 1-10 s. We note that such white dwarfs are candidates for the progenitors of peculiar, overluminous Type Ia supernovae, to which observational evidence ascribes mass in the range 2.1-2.8 solar mass. We also present some interesting results related to the structure of such white dwarfs, especially the existence of polar hollows in special cases.

Observational constraints on the unified dark matter and dark energy model based on the quark bag model

In this work we investigate if a small fraction of quarks and gluons, which escaped hadronization and survived as a uniformly spread perfect fluid, can play the role of both dark matter and dark energy. This fluid, as developed in [1], is characterized by two main parameters: beta, related to the amount of quarks and gluons which act as dark matter; and gamma, acting as the cosmological constant. We explore the feasibility of this model at cosmological scales using data from type Ia Supernovae (SNeIa), Long Gamma-Ray Bursts (LGRB) and direct observational Hubble data. We find that: (i) in general, beta cannot be constrained by SNeIa data nor by LGRB or H(z) data; (ii) gamma can be constrained quite well by all three data sets, contributing with approximate to 78% to the energy matter content; (iii) when a strong prior on (only) baryonic matter is assumed, the two parameters of the model are constrained successfully. (C) 2014 The Authors. Published by Elsevier B.V.