Exponential Data Fitting and Its Applications

Type Ia supernova light curves are characterized by a rapid rise from zero luminosity to a peak value, followed by a slower quasi-exponential decline. The rise and peak last for a few days, while the decline persists for many months. It is widely believed that the decline is powered by the radioactive decay chain 56Ni → 56Co → 56Fe, but the rates of decline in luminosity do not exactly match the decay rates of Ni and Co. In 1976, Rust, Leventhal, and McCall [19] presented evidence that the declining part of the light curve is well modelled by a linear combination of two exponentials whose decay rates were proportional to, but not exactly equal to, the decay rates for Ni and Co. The proposed reason for the lack of agreement between the rates was that the radioactive decays take place in the interior of a white dwarf star, at densities much higher than any encountered in a terrestrial environment, and that these higher densities accelerate the two decays by the same factor. This paper revisits this model, demonstrating that a variant of it provides excellent fits to observed luminosity data from 6 supernovae.

[1]  J. V. Hise Light-decay curve of the supernova in IC 4182 , 1974 .

[2]  Almerico Murli,et al.  An implementation of a Fourier series method for the numerical inversion of the Laplace transform , 1999, TOMS.

[3]  Taka Sakurai,et al.  Numerical inversion for Laplace transforms of functions with discontinuities , 2004, Advances in Applied Probability.

[4]  David Branch,et al.  Type Ia Supernovae as Standard Candles , 1993 .

[5]  H. A. Watts,et al.  DEPAC - design of a user oriented package of ODE solvers , 1980 .

[6]  B. Rust Use of supernovae light curves for testing the expansion hypothesis and other cosmological relations , 1974 .

[7]  BVRI Light Curves for 29 Type Ia Supernovae , 1996, astro-ph/9609064.

[8]  H. Gudbjartsson,et al.  The rician distribution of noisy mri data , 1995, Magnetic resonance in medicine.

[9]  J. Abate,et al.  Multi‐precision Laplace transform inversion , 2004 .

[10]  J. A. C. Weideman,et al.  Algorithms for Parameter Selection in the Weeks Method for Inverting the Laplace Transform , 1999, SIAM J. Sci. Comput..

[11]  M. Paluszny,et al.  A Multi-strategy Method for MRI Segmentation , 2009 .

[12]  Koenraad Van Leemput,et al.  A unifying framework for partial volume segmentation of brain MR images , 2003, IEEE Transactions on Medical Imaging.

[13]  B. Davies,et al.  Numerical Inversion of the Laplace Transform: A Survey and Comparison of Methods , 1979 .

[14]  Ward Whitt,et al.  On the Laguerre Method for Numerically Inverting Laplace Transforms , 1996, INFORMS J. Comput..

[15]  A. Talbot The Accurate Numerical Inversion of Laplace Transforms , 1979 .

[16]  Neil A. Thacker,et al.  The Importance of Partial Voluming in Multi-dimensional Medical Image Segmentation , 2001, MICCAI.

[17]  Joseph P. Hornak,et al.  From impractical to practical: solving an MRI problem using parallelism , 2005 .

[18]  B. Rust,et al.  Evidence for a radioactive decay hypothesis for supernova luminosity , 1976, Nature.

[19]  S. Provencher CONTIN: A general purpose constrained regularization program for inverting noisy linear algebraic and integral equations , 1984 .

[20]  A. M. Cohen Numerical Methods for Laplace Transform Inversion , 2007 .

[21]  A. Pastorello,et al.  SN 2003du: 480 days in the life of a normal type Ia supernova , 2007, 0704.1244.

[22]  Brian Davies,et al.  Integral transforms and their applications , 1978 .

[23]  Gene H. Golub,et al.  The differentiation of pseudo-inverses and non-linear least squares problems whose variables separate , 1972, Milestones in Matrix Computation.

[24]  M. Phillips,et al.  The Absolute Magnitudes of Type IA Supernovae , 1993 .

[25]  Salvatore Cuomo,et al.  Computation of the inverse Laplace transform based on a collocation method which uses only real values , 2007 .

[26]  S. Provencher A constrained regularization method for inverting data represented by linear algebraic or integral equations , 1982 .

[27]  M. Paluszny,et al.  A Quasi-Analytical Method for Relaxation Rate Distribution Determination of T2-Weighted MRI in Brain , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[28]  Ward Whitt,et al.  The Fourier-series method for inverting transforms of probability distributions , 1992, Queueing Syst. Theory Appl..