Source author record

E. Baron

E. Baron appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.HE astro-ph astro-ph.CO astro-ph.IM astro-ph.EP astro-ph.GA Computational Engineering, Finance, and Science physics.ao-ph physics.flu-dyn physics.space-ph

Catalog footprint

What is connected

60works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A Speed Bump: SN 2021aefx Shows that Doppler Shift Alone can Explain Early-Excess Blue Flux in Some Type Ia Supernovae

We present early-time photometric and spectroscopic observations of the Type Ia Supernova (SN Ia) 2021aefx. The early time u-band light curve shows an excess flux when compared to normal SNe Ia. We suggest that the early-excess blue flux may be due to a rapid change in spectral velocity in the first few days post explosion, produced by the emission of the Ca II H&K feature passing from the u to the B bands on the time scale of a few days. This effect could be dominant for all SNe Ia which have broad absorption features and early-time velocities over 25,000 km/s. It is likely to be one of the main causes of early-excess u-band flux in SNe Ia which have early-time high-velocities. This effect may also be dominant in the UV filters, as well as in places where the SN spectral energy distribution is quickly rising to longer wavelengths. The rapid change in velocity can only produce a monotonic change (in flux-space) in the u-band. For objects which explode at lower velocities, and have a more structured shape in the early-excess emission, there must also be an additional parameter producing the early-time diversity. More early time observations, in particular early spectra, are required to determine how prominent this effect is within SNe Ia.

preprint2022arXiv

A Tale of Two Type Ia Supernovae: The fast-declining siblings SNe 2015bo and 1997cn

We present optical and near-infrared photometric and spectroscopic observations of the fast-declining Type Ia Supernova (SN) 2015bo. SN 2015bo is under-luminous (M$_B$ = -17.67 $\pm$ 0.15 mag) and has a fast-evolving light curve ($Δ\mathrm{m}_{15}\mathrm{(B)}$ = 1.91 $\pm$ 0.01 mag and $s_{BV}$ = 0.48 $\pm$ 0.01). It has a unique morphology in the $V-r$ color curve, where it is bluer than all other SNe in the comparison sample. A $^{56}$Ni mass of 0.17 $\pm$ 0.03 $M_{\odot}$ was derived from the peak bolometric luminosity, which is consistent with its location on the luminosity-width relation. Spectroscopically, SN 2015bo is a Cool SN in the Branch classification scheme. The velocity evolution measured from spectral features is consistent with 1991bg-like SNe. SN 2015bo has a SN twin (similar spectra) and sibling (same host galaxy), SN 1997cn. Distance moduli of $μ$ = 34.33 $\pm$ 0.01 (stat) $\pm$0.11 (sys) mag and $μ$ = 34.34 $\pm$ 0.04 (stat) $\pm$ 0.12 (sys) mag were derived for SN 2015bo and SN 1997cn, respectively. These distances are consistent at the 0.06-$σ$ level with each other, and are also consistent with distances derived using surface-brightness fluctuations and redshift-corrected cosmology. This suggests that fast-declining SNe could be accurate distance indicators which should not be excluded from future cosmological analyses.

preprint2022arXiv

SN 2019va: A Type IIP Supernova with Large Influence of Nickel-56 Decay on the Plateau-phase Light Curve

We present multi-band photometric and spectroscopic observations of the type II supernova, (SN) 2019va, which shows an unusually flat plateau-phase evolution in its V-band light curve. Its pseudo-bolometric light curve even shows a weak brightening towards the end of the plateau phase. These uncommon features are related to the influence of 56Ni decay on the light curve during the plateau phase, when the SN emission is usually dominated by cooling of the envelope. The inferred 56Ni mass of SN 2019va is 0.088+/-0.018 solar mass, which is significantly larger than most SNe II. To estimate the influence of 56Ni decay on the plateau-phase light curve, we calculate the ratio (dubbed as eta_Ni) between the integrated time-weighted energy from 56Ni decay and that from envelope cooling within the plateau phase, obtaining a value of 0.8 for SN 2019va, which is the second largest value among SNe II that have been measured. After removing the influence of 56Ni decay on the plateau-phase light curve, we found that the progenitor/explosion parameters derived for SN 2019va are more reasonable. In addition, SN 2019va is found to have weaker metal lines in its spectra compared to other SNe IIP at similar epochs, implying a low-metallicity progenitor, which is consistent with the metal-poor environment inferred from the host-galaxy spectrum. We further discuss the possible reasons that might lead to SN 2019va-like events.

preprint2021arXiv

A 3D radiative transfer framework: XII. Many-core, vector and GPU methods

3D detailed radiative transfer is computationally taxing, since the solution of the radiative transfer equation involves traversing the six dimensional phase space of the 3D domain. With modern supercomputers the hardware available for wallclock speedup is rapidly changing, mostly in response to requirements to minimize the cost of electrical power. Given the variety of modern computing architectures, we aim to develop and adapt algorithms for different computing architectures to improve performance on a wide variety of platforms. We implemented the main time consuming kernels for solving 3D radiative transfer problems for vastly different computing architectures using MPI, OpenMP, OpenACC and vector algorithms. Adapted algorithms lead to massively improved speed for all architectures, making extremely large model calculations easily feasible. These calculations would have previously been considered impossible or prohibitively expensive. Efficient use of modern computing devices is entirely feasible, but unfortunately requires the implementation of specialized algorithms for them.

preprint2021arXiv

Carnegie Supernova Project-II: Near-infrared Spectroscopy of Stripped-Envelope Core-Collapse Supernovae

We present 75 near-infrared (NIR; 0.8$-$2.5 $μ$m) spectra of 34 stripped-envelope core-collapse supernovae (SESNe) obtained by the Carnegie Supernova Project-II (CSP-II), encompassing optical spectroscopic Types IIb, Ib, Ic, and Ic-BL. The spectra range in phase from pre-maximum to 80 days past maximum. This unique data set constitutes the largest NIR spectroscopic sample of SESNe to date. NIR spectroscopy provides observables with additional information that is not available in the optical. Specifically, the NIR contains the resonance lines of He I and allows a more detailed look at whether Type Ic supernovae are completely stripped of their outer He layer. The NIR spectra of SESNe have broad similarities, but closer examination through statistical means reveals a strong dichotomy between NIR "He-rich" and "He-poor" SNe. These NIR subgroups correspond almost perfectly to the optical IIb/Ib and Ic/Ic-BL types, respectively. The largest difference between the two groups is observed in the 2 $μ$m region, near the He I $λ$2.0581 $μ$m line. The division between the two groups is not an arbitrary one along a continuous sequence. Early spectra of He-rich SESNe show much stronger He I $λ$2.0581 $μ$m absorption compared to the He-poor group, but with a wide range of profile shapes. The same line also provides evidence for trace amounts of He in half of our SNe in the He-poor group.

preprint2020arXiv

Carnegie Supernova Project II: The slowest rising Type Ia supernova LSQ14fmg and clues to the origin of super-Chandrasekhar/03fg-like events

The Type Ia supernova (SN Ia) LSQ14fmg exhibits exaggerated properties which may help to reveal the origin of the "super-Chandrasekhar" (or 03fg-like) group. The optical spectrum is typical of a 03fg-like SN Ia, but the light curves are unlike those of any SNe Ia observed. The light curves of LSQ14fmg rise extremely slowly. At -23 rest-frame days relative to B-band maximum, LSQ14fmg is already brighter than $M_V$=-19 mag before host extinction correction. The observed color curves show a flat evolution from the earliest observation to approximately one week after maximum. The near-infrared light curves peak brighter than -20.5 mag in the J and H bands, far more luminous than any 03fg-like SNe Ia with near-infrared observations. At one month past maximum, the optical light curves decline rapidly. The early, slow rise and flat color evolution are interpreted to result from an additional excess flux from a power source other than the radioactive decay of the synthesized $^{56}Ni$. The excess flux matches the interaction with a typical superwind of an asymptotic giant branch (AGB) star in density structure, mass-loss rate, and duration. The rapid decline starting at around one month past B-band maximum may be an indication of rapid cooling by active carbon monoxide (CO) formation, which requires a low temperature and high density environment. These peculiarities point to an AGB progenitor near the end of its evolution and the core degenerate scenario as the likely explosion mechanism for LSQ14fmg.

preprint2020arXiv

Carnegie Supernova Project-II: A new method to photometrically identify sub-types of extreme Type Ia Supernovae

We present a new method to photometrically delineate between various sub-types of type Ia supernovae (SNe Ia). Using the color-stretch parameters, $s_{BV}$ or $s_{gr}$, and the time of i-band primary maximum relative to the B-band or g-band maximum it is demonstrated that 2003fg-like, 1991bg-like, and 2002cx-like SNe Ia can readily be identified. In the cases of these extreme SNe Ia, their primary i-band maximum occurs after the time of the B or g band maxima. We suggest that the timing of the i-band maximum can reveal the physical state of the SN Ia explosion as it traces: i) the speed of the recombination front of iron group elements in the ejecta, ii) the temperature evolution and rate of adiabatic cooling in the ejecta and, iii) the presence of interaction with a stellar envelope. This photometric sub-typing can be used in conjunction with other SNe analysis, such as the Branch diagram, to examine the physics and diversity of SNe Ia. The results here can also be used to screen out non-Ia SNe from cosmological samples that do not have complete spectroscopic typing. Finally, as future surveys like LSST create large databases of light curves of many objects this photometric identification can be used to readily identify and study the rates and bulk properties of peculiar SNe Ia.

preprint2020arXiv

Discovery and Rapid Follow-up Observations of the Unusual Type II SN 2018ivc in NGC 1068

We present the discovery and high-cadence follow-up observations of SN 2018ivc, an unusual Type II supernova that exploded in NGC 1068 (D=10.1 Mpc). The light curve of SN 2018ivc declines piecewise-linearly, changing slope frequently, with four clear slope changes in the first 30 days of evolution. This rapidly changing light curve indicates that interaction between the circumstellar material and ejecta plays a significant role in the evolution. Circumstellar interaction is further supported by a strong X-ray detection. The spectra are rapidly evolving and dominated by hydrogen, helium, and calcium emission lines. We identify a rare high-velocity emission-line feature blueshifted at ~7800 km/s (in Ha, Hb, Pb, Pg, HeI, CaII), which is visible from day 18 until at least day 78 and could be evidence of an asymmetric progenitor or explosion. From the overall similarity between SN 2018ivc and SN 1996al, the \Ha{} equivalent width of its parent HII region, and constraints from pre-explosion archival Hubble Space Telescope images, we find that the progenitor of SN 2018ivc could be as massive as 52 Msun but is more likely <12 Msun. SN 2018ivc demonstrates the importance of the early discovery and rapid follow-up observations of nearby supernovae to study the physics and progenitors of these cosmic explosions.

preprint2020arXiv

SN 2018hti: a nearby superluminous supernova discovered in a metal-poor galaxy

SN 2018hti is a Type I superluminous supernova (SLSN~I) with an absolute $g$-band magnitude of $-22.2$ at maximum brightness, discovered in a metal-poor galaxy at a redshift of 0.0612. We present extensive photometric and spectroscopic observations of this supernova, covering the phases from $\sim -35$ days to more than +340 days from the $r$-band maximum. Combining our $BVgri$-band photometry with {\it Swift} UVOT optical/ultraviolet photometry, we calculated the peak luminosity as $\sim 3.5\times10^{44}$ erg s$^{-1}$. Modeling the observed light curve reveals that the luminosity evolution of SN 2018hti can be produced by an ejecta mass of 5.8 $M_\odot$ and a magnetar with a magnetic field of $B=1.8\times10^{13}$~G having an initial spin period of $P_0=1.8$ ms. Based on such a magnetar-powered scenario and a larger sample, a correlation between the spin of the magnetar and the kinetic energy of the ejecta can be inferred for most SLSNe~I, suggesting a self-consistent scenario. Like for other SLSNe~I, the host galaxy of SN 2018hti is found to be relatively faint ($M_{g} = -17.75$ mag) and of low metallicity ($Z=0.3~Z_\odot$), with a star-formation rate of 0.3 $M_\odot$ yr$^{-1}$. According to simulation results of single-star evolution, SN 2018hti could originate from a massive, metal-poor star with a zero-age main sequence (ZAMS) mass of 25--40 $M_\odot$, or from a less massive rotating star with $M_\mathrm{ZAMS} \approx 16$--25 $M_\odot$. For the case of a binary system, its progenitor could also be a star with $M_\mathrm{ZAMS} \gtrsim 25$ $M_\odot$.

preprint2020arXiv

The Carnegie Supernova Project II. Observations of the intermediate luminosity red transient SNhunt120

We present multi-wavelength observations of two gap transients followed by the Carnegie Supernova Project-II and supplemented with data obtained by a number of different programs. Here in the first of two papers, we focus on the intermediate luminosity red transient (ILRT) designated SNhunt120, while in a companion paper we examine the luminous red novae AT 2014ej. Our data set for SNhunt120 consists of an early optical discovery, estimated to be within 3 days after outburst, the subsequent optical and near-infrared broadband followup extending over a $\sim$2 month period, two visual- and two near-infrared wavelength spectra, and Spitzer Space Telescope observations extending from early ($+$28 d) to late ($+$1155 d) phases. SNhunt120 resembles other ILRTs such as NGC 300-2008-OT and SN 2008S, and like these other ILRTs, SNhunt120 exhibits prevalent mid-infrared emission at both early and late phases. From the comparison of SNhunt120 and other ILRTs to electron-capture supernova simulations, we find that the current models underestimate the explosion kinetic energy and thereby produce synthetic light curves that over-estimate the luminosity. Finally, examination of pre-outburst Hubble Space Telescope images yields no progenitor detection.

preprint2020arXiv

The Carnegie Supernova Project II. Observations of the luminous red nova AT 2014ej

We present optical and near-infrared broadband photometry and optical spectra of AT 2014ej from the the Carnegie Supernova Project-II. These observations are complemented with data from the CHilean Automatic Supernova sEarch, the Public ESO Spectroscopic Survey of Transient Objects, and from the Backyard Observatory Supernova Search. Observational signatures of AT 2014ej reveal that it is similar to other members of the gap-transient subclass known as luminous red novae (LRNe), including the ubiquitous double hump light curve and spectral properties akin to the LRN SN 2017jfs. A medium-dispersion, visual-wavelength spectrum of AT 2014ej taken the Magellan Clay telescope, exhibits a P Cygni H$α$ feature characterized by a blue velocity at zero intensity of $\approx 110$ km s$^{-1}$ and a P Cygni minimum velocity of $\approx70$ km s$^{-1}$, and which we attribute to emission from a circumstellar wind. Inspection of pre-outbust Hubble Space Telescope images yields no conclusive progenitor detection. In comparison with a sample of LRNe from the literature, AT 2014ej lies among the brighter end of the luminosity distribution. Comparison of the ultra-violet, optical, infrared (UVOIR) light curves of well-observed LRNe to common-envelope evolution models from the literature, indicates the models under predict the luminosity of the comparison sample at all phases and also produce inconsistent time-scales of the secondary peak. Future efforts to model LRNe should expand upon the current parameter space explored and therefore may consider more massive systems and a wider range of dynamical timescales.

preprint2019arXiv

Predicting the Extreme Ultraviolet Radiation Environment of Exoplanets Around Low-Mass Stars: GJ 832, GJ 176, GJ 436

Correct estimates of stellar extreme ultraviolet (EUV; 100 - 1170 Å) flux are important for studying the photochemistry and stability of exoplanet atmospheres, as EUV radiation ionizes hydrogen and contributes to the heating, expansion, and potential escape of a planet's upper atmosphere. Contamination from interstellar hydrogen makes observing EUV emission from M stars particularly difficult, and impossible past 100 pc, and necessitates other means to predict the flux in this wavelength regime. We present EUV -- infrared (100 Å- 5.5 $μ$m) synthetic spectra computed with the PHOENIX atmospheric code of three early M dwarf planet hosts: GJ 832 (M1.5 V), GJ 176 (M2.5 V), and GJ 436 (M3.5 V). These one-dimensional semiempirical nonlocal thermodynamic equilibrium models include simple temperature prescriptions for the stellar chromosphere and transition region, from where ultraviolet (UV; 100 - 3008 Å) fluxes originate. We guide our models with Hubble Space Telescope far- and near-UV spectra and discuss the ability to constrain these models using Galaxy Evolution Explorer UV photometry. Our models closely reproduce the observations and predict the unobservable EUV spectrum at a wavelength resolution of < 0.1 Å. The temperature profiles that best reproduce the observations for all three stars are described by nearly the same set of parameters, suggesting that early M type stars may have similar thermal structures in their upper atmospheres. With an impending UV observation gap and the scarcity of observed EUV spectra for stars less luminous and more distant than the Sun, upper-atmosphere models such as these are important for providing realistic spectra across short wavelengths and for advancing our understanding of the effects of radiation on planets orbiting M stars.

preprint2019arXiv

The Carnegie Supernova Project II. Early observations and progenitor constraints of the Type Ib supernova LSQ13abf

Supernova LSQ13abf was discovered soon after explosion by the La Silla-QUEST Survey and followed by the CSP II at optical and near-IR wavelengths. Our analysis indicates LSQ13abf was discovered within two days of explosion and its first 10 days of evolution reveal a B-band light curve with an abrupt drop in luminosity. Contemporaneously, the V-band light curve exhibits a rise towards a first peak and the r- and i-band light curves show no early peak. The early light-curve evolution of LSQ13abf is reminiscent of the post explosion cooling phase observed in the Type Ib SN 2008D, and the similarity between the two objects extends over weeks. Spectroscopically, LSQ13abf resembles SN 2008D with P Cygni He I features that strengthen over time. Spectral energy distributions are constructed from broad-band photometry, and by fitting black-body (BB) functions a UVOIR light curve is constructed, and the underlying BB-temperature and BB-radius profiles are estimated. Explosion parameters are estimated by simultaneously fitting an Arnett model to the UVOIR light curve and the velocity evolution derived from spectral features, and a post-shock breakout cooling model to the first two epochs of the bolometric evolution. This combined model suggests an explosion energy of 1.3x10$^{51}$ ergs, a relatively high ejecta mass of 5.94 M$_{\odot}$, a Ni mass of 0.16 M$_{\odot}$, and a progenitor-star radius of 28.0 R$_{\odot}$. The ejecta mass suggests the origins of LSQ13abf lie with a >25 M$_{\odot}$ ZAMS progenitor and its radius is three and nine times larger than values estimated from the same analysis applied to observations of SNe 2008D and 1999ex, respectively. Alternatively, comparison of hydrodynamical simulations of >20-25 M$_{\odot}$ ZAMS progenitors that evolve to pre-SN envelope masses around 10 M$_{\odot}$ and extended (~100 R$_{\odot}$) envelopes also match the observations of LSQ13abf.

preprint2016arXiv

Supernova 2014J at M82: II. Direct Analysis of Spectra Obtained with Isaac Newton and William Herschel Telescopes

We analyze a time series of optical spectra of SN~2014J from almost two weeks prior to maximum to nearly four months after maximum. We perform our analysis using the SYNOW code, which is well suited to track the distribution of the ions with velocity in the ejecta. We show that almost all of the spectral features during the entire epoch can be identified with permitted transitions of the common ions found in normal SNe Ia in agreement with previous studies. We show that 2014J is a relatively normal SN Ia. At early times the spectral features are dominated by Si II, S II, Mg II, and Ca II. These ions persist to maximum light with the appearance of Na I and Mg I. At later times iron-group elements also appear, as expected in the stratified abundance model of the formation of normal type Ia SNe. We do not find significant spectroscopic evidence for oxygen, until 100 days after maximum light. The +100 day identification of oxygen is tentative, and would imply significant mixing of unburned or only slight processed elements down to a velocity of 6,000 km/s. Our results are in relatively good agreement with other analyses in the IR. We briefly compare SN 2011fe to SN 2014J and conclude that the differences could be due to different central densities at ignition or differences in the C/O ratio of the progenitors.

preprint2015arXiv

Identification of the feature that causes the I-band secondary maximum of a type Ia supernova

We obtained a time series of spectra covering the secondary maximum in the I-band of the bright Type Ia supernova 2014J in M82 with the TIGRE telescope. Comparing the observations with theoretical models calculated with the time dependent extension of the PHOENIX code, we identify the feature that causes the secondary maximum in the I-band light curve. Fe II 3d6(3D)4s-3d6(5D)4p and similar high excitation transitions produce a blended feature at 7500 Å, which causes the rise of the light curve towards the secondary maximum. The series of observed spectra of SN 2014J and archival data of SN 2011fe confirm this conclusion. We further studied the plateau phase of the Rband light curve of SN 2014J and searched for features which contribute to the flux. The theoretical models do not clearly indicate a new feature that may cause the Rband plateau phase. However, Co II features in the range of 6500 - 7000 Å and the Fe II feature of the I-band are clearly seen in the theoretical spectra, but do not appear to provide all of the flux necessary for the R-band plateau.

preprint2015arXiv

Spectral Models for Early Time SN 2011fe Observations

We use observed UV through near IR spectra to examine whether SN 2011fe can be understood in the framework of Branch-normal SNe Ia and to examine its individual peculiarities. As a benchmark, we use a delayed-detonation model with a progenitor metallicity of Z_solar/20. We study the sensitivity of features to variations in progenitor metallicity, the outer density profile, and the distribution of radioactive nickel. The effect of metallicity variations in the progenitor have a relatively small effect on the synthetic spectra. We also find that the abundance stratification of SN 2011fe resembles closely that of a delayed detonation model with a transition density that has been fit to other Branch-normal Type Ia supernovae. At early times, the model photosphere is formed in material with velocities that are too high, indicating that the photosphere recedes too slowly or that SN 2011fe has a lower specific energy in the outer ~0.1 M_sun than does the model. We discuss several explanations for the discrepancies. Finally, we examine variations in both the spectral energy distribution and in the colors due to variations in the progenitor metallicity, which suggests that colors are only weak indicators for the progenitor metallicity, in the particular explosion model that we have studied. We do find that the flux in the U band is significantly higher at maximum light in the solar metallicity model than in the lower metallicity model and the lower metallicity model much better matches the observed spectrum.

preprint2015arXiv

Theoretical Clues to the Ultraviolet Diversity of Type Ia Supernovae

The effect of metallicity on the observed light of Type Ia supernovae (SNe Ia) could lead to systematic errors as the absolute magnitudes of local and distant SNe Ia are compared to measure luminosity distances and determine cosmological parameters. The UV light may be especially sensitive to metallicity, though different modeling methods disagree as to the magnitude, wavelength dependence, and even the sign of the effect. The outer density structure, ^56 Ni, and to a lesser degree asphericity, also impact the UV. We compute synthetic photometry of various metallicity-dependent models and compare to UV/optical photometry from the Swift Ultra-Violet/Optical Telescope. We find that the scatter in the mid-UV to near-UV colors is larger than predicted by changes in metallicity alone and is not consistent with reddening. We demonstrate that a recently employed method to determine relative abundances using UV spectra can be done using UVOT photometry, but we warn that accurate results require an accurate model of the cause of the variations. The abundance of UV photometry now available should provide constraints on models that typically rely on UV spectroscopy for constraining metallicity, density, and other parameters. Nevertheless, UV spectroscopy for a variety of SN explosions is still needed to guide the creation of accurate models. A better understanding of the influences affecting the UV is important for using SNe Ia as cosmological probes, as the UV light may test whether SNe Ia are significantly affected by evolutionary effects.

preprint2014arXiv

A 3D radiative transfer framework: XI. multi-level NLTE

Multi-level non-local thermodynamic equilibrium (NLTE) radiation transfer calculations have become standard throughout the stellar atmospheres community and are applied to all types of stars as well as dynamical systems such as novae and supernovae. Even today spherically symmetric 1D calculations with full physics are computationally intensive. We show that full NLTE calculations can be done with fully 3 dimensional (3D) radiative transfer. With modern computational techniques and current massive parallel computational resources, full detailed solution of the multi-level NLTE problem coupled to the solution of the radiative transfer scattering problem can be solved without sacrificing the micro physics description. We extend the use of a rate operator developed to solve the coupled NLTE problem in spherically symmetric 1D systems. In order to spread memory among processors we have implemented the NLTE/3D module with a hierarchical domain decomposition method that distributes the NLTE levels, radiative rates, and rate operator data over a group of processes so that each process only holds the data for a fraction of the voxels. Each process in a group holds all the relevant data to participate in the solution of the 3DRT problem so that the 3DRT solution is parallelized within a domain decomposition group. We solve a spherically symmetric system in 3D spherical coordinates in order to directly compare our well-tested 1D code to the 3D case. We compare three levels of tests: a) a simple H+He test calculation, b) H+He+CNO+Mg, c) H+He+Fe. The last test is computationally large and shows that realistic astrophysical problems are solvable now, but they do require significant computational resources. With presently available computational resources it is possible to solve the full 3D multi-level problem with the same detailed micro-physics as included in 1D modeling.

preprint2014arXiv

Comprehensive Observations of the Bright and Energetic Type Iax SN 2012Z: Interpretation as a Chandrasekhar Mass White Dwarf Explosion

We present UV through NIR broad-band photometry, and optical and NIR spectroscopy of Type Iax supernova 2012Z. The data set consists of both early and late-time observations, including the first late phase NIR spectrum obtained for a spectroscopically classified SN Iax. Simple model calculations of its bolometric light curve suggest SN 2012Z produced ~0.3 M_sun of (56)Ni, ejected about a Chandrasekhar mass of material, and had an explosion energy of ~10^51 erg, making it one of the brightest and most energetic SN Iax yet observed. The late phase NIR spectrum of SN 2012Z is found to broadly resemble similar epoch spectra of normal SNe Ia; however, like other SNe Iax, corresponding visual-wavelength spectra differ substantially compared to all supernova types. Constraints from the distribution of IMEs, e.g. silicon and magnesium, indicate that the outer ejecta did not experience significant mixing during or after burning, and the late phase NIR line profiles suggests most of the (56)Ni is produced during high density burning. The various observational properties of SN 2012Z are found to be consistent with the theoretical expectations of a Chandrasekhar mass white dwarf progenitor that experiences a pulsational delayed detonation, which produced several tenths of a solar mass of (56)Ni during the deflagration burning phase and little (or no) (56)Ni during the detonation phase. Within this scenario only a moderate amount of Rayleigh-Taylor mixing occurs both during the deflagration and fallback phase of the pulsation, and the layered structure of the IMEs is a product of the subsequent denotation phase. The fact that the SNe Iax population does not follow a tight brightness-decline relation similar to SNe Ia can then be understood in the framework of variable amounts of mixing during pulsational rebound and variable amounts of (56)Ni production during the early subsonic phase of expansion.

preprint2014arXiv

Near-infrared line identification in type Ia supernovae during the transitional phase

We present near-infrared synthetic spectra of a delayed-detonation hydrodynamical model and compare them to observed spectra of four normal type Ia supernovae ranging from day +56.5 to day +85. This is the epoch during which supernovae are believed to be undergoing the transition from the photospheric phase, where spectra are characterized by line scattering above an optically thick photosphere, to the nebular phase, where spectra consist of optically thin emission from forbidden lines. We find that most spectral features in the near-infrared can be accounted for by permitted lines of Fe II and Co II. In addition, we find that [Ni II] fits the emission feature near 1.98 μm, suggesting that a substantial mass of 58Ni exists near the center of the ejecta in these objects, arising from nuclear burning at high density. A tentative identification of Mn II at 1.15 μm may support this conclusion as well.

preprint2014arXiv

SNe Ia: Can Chandrasekhar Mass Explosions Reproduce the Observed Zoo?

The question of the nature of the progenitor of Type Ia supernovae (SNe Ia) is important both for our detailed understanding of stellar evolution and for their use as cosmological probes of the dark energy. Much of the basic features of SNe Ia can be understood directly from the nuclear physics, a fact which Gerry would have appreciated. We present an overview of the current observational and theoretical situation and show that it not incompatible with most SNe Ia being the results of thermonuclear explosions near the Chandrasekhar mass.

preprint2013arXiv

Hydrogen and helium in the spectra of Type Ia supernovae

We present predictions for hydrogen and helium emission line luminosities from circumstellar matter around Type Ia supernovae (SNe Ia) using time dependent photoionization modeling. ESO/VLT optical echelle spectra of the SN Ia 2000cx were taken before and up to 70 days after maximum. We detect no hydrogen and helium lines, and place an upper limit on the mass loss rate for the putative wind of less than 1.3EE{-5} solar masses per year, assuming a speed of 10 km/s and solar abundances for the wind. In a helium-enriched case, the best line to constrain the mass loss would be He I 10,830 A. We confirm the details of interstellar Na I and Ca II absorption towards SN 2000cx as discussed by Patat et al., but also find evidence for 6613.56 A Diffuse Interstellar Band (DIB) absorption in the Milky Way. We discuss measurements of the X-ray emission from the interaction between the supernova ejecta and the wind and we re-evaluate observations of SN 1992A obtained 16 days after maximum by Schlegel & Petre. We find an upper limit of 1.3EE{-5} solar masses per year. These results, together with the previous observational work on the normal SNe Ia 1994D and 2001el, disfavour a symbiotic star in the upper mass loss rate regime from being the likely progenitor scenario for these SNe. To constrain hydrogen in late time spectra, we present ESO/VLT and ESO/NTT optical and infrared observations of SNe Ia 1998bu and 2000cx 251-388 days after maximum. We see no hydrogen line emission in SNe 1998bu and 2000cx at these epochs, and we argue from modeling that the mass of such hydrogen-rich gas must be less than 0.03 solar masses for both supernovae. Comparing similar upper limits with recent models of Pan et al., it seems hydrogen-rich donors with a separation of less than 5 times the radius of the donor may be ruled out for the five SNe Ia 1998bu, 2000cx, 2001el, 2005am and 2005cf.

preprint2013arXiv

Optical and Near-IR Observations of the Faint and Fast 2008ha-like Supernova 2010ae

A comprehensive set of optical and near-infrared (NIR) photometry and spectroscopy is presented for the faint and fast 2008ha-like supernova (SN) 2010ae. Contingent on the adopted value of host extinction, SN 2010ae reached a peak brightness of -13.8 > M_V > -15.3 mag, while modeling of the UVOIR light curve suggests it produced 0.003--0.007 M_sun of (56)Ni, ejected 0.30--0.60 M_sun of material, and had an explosion energy of 0.04--0.30x10^51 erg. The values of these explosion parameters are similar to the peculiar SN 2008ha --for which we also present previously unpublished early phase optical and NIR light curves-- and places these two transients at the faint end of the 2002cx-like SN population. Detailed inspection of the post maximum NIR spectroscopic sequence indicates the presence of a multitude of spectral features, which are identified through SYNAPPS modeling to be mainly attributed to Co II. Comparison with a collection of published and unpublished NIR spectra of other 2002cx-like SNe, reveals that a Co II footprint is ubiquitous to this subclass of transients, providing a link to Type Ia SNe. A visual-wavelength spectrum of SN 2010ae obtained at +252 days past maximum shows a striking resemblance to a similar epoch spectrum of SN 2002cx. However, subtle differences in the strength and ratio of calcium emission features, as well as diversity among similar epoch spectra of other 2002cx-like SNe indicates a range of physical conditions of the ejecta, highlighting the heterogeneous nature of this peculiar class of transients.

preprint2012arXiv

A 3D radiative transfer framework IX. Time dependence

Context. Time-dependent, 3D radiation transfer calculations are important for the modeling of a variety of objects, from supernovae and novae to simulations of stellar variability and activity. Furthermore, time-dependent calculations can be used to obtain a 3D radiative equilibrium model structure via relaxation in time. Aims. We extend our 3D radiative transfer framework to include direct time dependence of the radiation field; i.e., the $\partial I/ \partial t$ terms are fully considered in the solution of radiative transfer problems. Methods. We build on the framework that we have described in previous papers in this series and develop a subvoxel method for the $\partial I/\partial t$ terms. Results. We test the implementation by comparing the 3D results to our well tested 1D time dependent radiative transfer code in spherical symmetry. A simple 3D test model is also presented. Conclusions. The 3D time dependent radiative transfer method is now included in our 3D RT framework and in PHOENIX/3D.

preprint2012arXiv

A 3D radiative transfer framework: X. Arbitrary Velocity Fields in the Co-moving Frame

3-D astrophysical atmospheres will have random velocity fields. We seek to combine the methods we have developed for solving the 1-D problem with arbitrary flows to those that we have developed for solving the fully 3-D relativistic radiative transfer problem in the case of monotonic flows. The methods developed in the case of 3-D atmospheres with monotonic flows, solving the fully relativistic problem along curves defined by an affine parameter, are very flexible and can be extended to the case of arbitrary velocity fields in 3-D. Simultaneously, the techniques we developed for treating the 1-D problem with arbitrary velocity fields are easily adapted to the 3-D problem. The algorithm we present allows the solution of 3-D radiative transfer problems that include arbitrary wavelength couplings. We use a quasi-analytic formal solution of the radiative transfer equation that significantly improves the overall computation speed. We show that the approximate lambda operator developed in previous work gives good convergence, even neglecting wavelength coupling. Ng acceleration also gives good results. We present tests that are of similar resolution to what has been presented using Monte-Carlo techniques, thus our methods will be applicable to problems outside of our test setup. Additional domain decomposition parallelization strategies will be explored in future work.

preprint2012arXiv

A Physical Model for SN 2001ay, a normal, bright, extremely slowly declining Type Ia supernova

We present a study of the peculiar Type Ia supernova 2001ay (SN 2001ay). The defining features of its peculiarity are: high velocity, broad lines, and a fast rising light curve, combined with the slowest known rate of decline. It is one magnitude dimmer than would be predicted from its observed value of Delta-m15, and shows broad spectral features. We base our analysis on detailed calculations for the explosion, light curves, and spectra. We demonstrate that consistency is key for both validating the models and probing the underlying physics. We show that this SN can be understood within the physics underlying the Delta-m15 relation, and in the framework of pulsating delayed detonation models originating from a Chandrasekhar mass, white dwarf, but with a progenitor core composed of 80% carbon. We suggest a possible scenario for stellar evolution which leads to such a progenitor. We show that the unusual light curve decline can be understood with the same physics as has been used to understand the Delta-m15 relation for normal SNe Ia. The decline relation can be explained by a combination of the temperature dependence of the opacity and excess or deficit of the peak luminosity, alpha, measured relative to the instantaneous rate of radiative decay energy generation. What differentiates SN 2001ay from normal SNe Ia is a higher explosion energy which leads to a shift of the Ni56 distribution towards higher velocity and alpha < 1. This result is responsible for the fast rise and slow decline. We define a class of SN 2001ay-like SNe Ia, which will show an anti-Phillips relation.

preprint2012arXiv

Analysis of the Early-Time Optical Spectra of SN 2011fe in M101

The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.

preprint2012arXiv

Near-infrared light curves of type Ia supernovae

Aims. With our time-dependent model atmosphere code PHOENIX, our goal is to simulate light curves and spectra of hydrodynamical models of all types of supernovae. In this work, we simulate near-infrared light curves of SNe Ia and confirm the cause of the secondary maximum. Methods. We apply a simple energy solver to compute the evolution of an SN Ia envelope during the free expansion phase. Included in the solver are energy changes due to expansion, the energy deposition of γ-rays and interaction of radiation with the material. Results. We computed theoretical light curves of several SN Ia hydrodynamical models in the I, J, H, and K bands and compared them to the observed SN Ia light curves of SN 1999ee and SN 2002bo. By changing a line scattering parameter in time, we obtained quite reasonable fits to the observed near-infrared light curves. This is a strong hint that detailed NLTE effects in IR lines have to be modeled, which will be a future focus of our work. Conclusions. We found that IR line scattering is very important for the near-infrared SN Ia light curve modeling. In addition, the recombination of Fe III to Fe II and of Co III to Co II is responsible for the secondary maximum in the near-infrared bands. For future work the consideration of NLTE for all lines (including the IR subordinate lines) will be crucial.

preprint2012arXiv

Quantitative Photospheric Spectral Analysis of the Type IIP Supernova 2007od

We compare and analyze a time series of spectral observations obtained during the first 30 days of evolution of SN 2007od with the non-LTE code PHOENIX. Despite some spectroscopic particularities in the Balmer features, this supernova appears to be a normal Type II, and the fits proposed are generally in good agreement with the observations. As a starting point we have carried out an analysis with the parameterized synthetic spectrum code SYNOW to confirm line identifications and to highlight differences between the results of the two codes. The analysis computed using PHOENIX suggests the presence of a high velocity feature in Hβ and an Hα profile reproduced with a density profile steeper than that of the other elements. We also show a detailed analysis of the ions velocities of the 6 synthetic spectra. The distance is estimated for each epoch with the Spectral-fitting Expanding Atmosphere Method (SEAM). Consistent results are found using all the spectra which give the explosion date of JD 2454403 (29 October, 2007) and a distance modulus μ= 32.2 \pm 0.3.

preprint2012arXiv

Supernova Resonance--scattering Line Profiles in the Absence of a Photosphere

In supernova spectroscopy relatively little attention has been given to the properties of optically thick spectral lines in epochs following the photosphere's recession. Most treatments and analyses of post-photospheric optical spectra of supernovae assume that forbidden-line emission comprises most if not all spectral features. However, evidence exists which suggests that some spectra exhibit line profiles formed via optically thick resonance-scattering even months or years after the supernova explosion. To explore this possibility we present a geometrical approach to supernova spectrum formation based on the "Elementary Supernova" model, wherein we investigate the characteristics of resonance-scattering in optically thick lines while replacing the photosphere with a transparent central core emitting non-blackbody continuum radiation, akin to the optical continuum provided by decaying 56Co formed during the explosion. We develop the mathematical framework necessary for solving the radiative transfer equation under these conditions, and calculate spectra for both isolated and blended lines. Our comparisons with analogous results from the Elementary Supernova code SYNOW reveal several marked differences in line formation. Most notably, resonance lines in these conditions form P Cygni-like profiles, but the emission peaks and absorption troughs shift redward and blueward, respectively, from the line's rest wavelength by a significant amount, despite the spherically symmetric distribution of the line optical depth in the ejecta. These properties and others that we find in this work could lead to misidentification of lines or misattribution of properties of line-forming material at post-photospheric times in supernova optical spectra.

preprint2012arXiv

The bright Type IIP SN 2009bw, showing signs of interaction

We present photometry and spectroscopy of the type IIP supernova 2009bw in UGC 2890 from few days after the outburst to 241 days. The light curve of SN 2009bw during the photospheric phase is similar to that of normal SNe IIP but with brighter peak and plateau (Mmax R = -17.82 mag, Mplateau R = -17.37 mag). The luminosity drop from the photospheric to the nebular phase is one of the fastest ever observed, ~2.2 mag in about 13 days. The radioactive tail of the bolometric light curve indicates that the amount of ejected 56 Ni is \approx 0.022 M\odot. The photospheric spectra reveal high velocity lines of Hα and Hβ until about 105 days after the shock breakout, suggesting a possible early interaction between the SN ejecta and pre-existent circumstellar material, and the presence of CNO elements. By modeling the bolometric light curve, ejecta expansion velocity and photospheric temperature, we estimate a total ejected mass of 8-12M\odot, a kinetic energy of ~0.3 foe and an initial radius of ~ 3.6 - 7 \times 10^13 cm.

preprint2012arXiv

The UV/optical spectra of the Type Ia supernova SN 2010jn: a bright supernova with outer layers rich in iron-group elements

Radiative transfer studies of Type Ia supernovae (SNe Ia) hold the promise of constraining both the time-dependent density profile of the SN ejecta and its stratification by element abundance which, in turn, may discriminate between different explosion mechanisms and progenitor classes. Here we present a detailed analysis of Hubble Space Telescope ultraviolet (UV) and ground-based optical spectra and light curves of the SN Ia SN 2010jn (PTF10ygu). SN 2010jn was discovered by the Palomar Transient Factory (PTF) 15 days before maximum light, allowing us to secure a time-series of four UV spectra at epochs from -11 to +5 days relative to B-band maximum. The photospheric UV spectra are excellent diagnostics of the iron-group abundances in the outer layers of the ejecta, particularly those at very early times. Using the method of 'Abundance Tomography' we have derived iron-group abundances in SN 2010jn with a precision better than in any previously studied SN Ia. Optimum fits to the data can be obtained if burned material is present even at high velocities, including significant mass fractions of iron-group elements. This is consistent with the slow decline rate (or high 'stretch') of the light curve of SN 2010jn, and consistent with the results of delayed-detonation models. Early-phase UV spectra and detailed time-dependent series of further SNe Ia offer a promising probe of the nature of the SN Ia mechanism.

preprint2011arXiv

A 3D radiative transfer framework: XIII. OpenCL implementation

We discuss an implementation of our 3D radiative transfer (3DRT) framework with the OpenCL paradigm for general GPU computing. We implement the kernel for solving the 3DRT problem in Cartesian coordinates with periodic boundary conditions in the horizontal $(x,y)$ plane, including the construction of the nearest neighbor $\Lstar$ and the operator splitting step. We present the results of a small and a large test case and compare the timing of the 3DRT calculations for serial CPUs and various GPUs. The latest available GPUs can lead to significant speedups for both small and large grids compared to serial (single core) computations.

preprint2011arXiv

Comparative Analysis of Peculiar Type Ia 1991bg-like Supernovae Spectra

Spectroscopic analyses of Type Ia supernovae have shown there exist four spectroscopic groups---cools, broad line, shallow silicon, and core normal---defined by the widths of the Si II features at 5972 Angstroms and 6355 Angstroms. 1991bg-likes are classified as "cools". Cools are dim, undergo a rapid decline in luminosity, and produce significantly less nickel than normal Type Ia supernovae. They also have an unusually deep and wide trough in their spectra around 4200 Angstroms and a relatively strong Si II absorption attributed to the line at 5972 Angstroms. We examine the spectra of supernova (SN) 1991bg and the cools SN 1997cn, SN 1999by, and SN 2005bl using the highly parameterized synthetic spectrum code SYNOW, and find general agreement with similar spectroscopic studies. Our analysis reveals that this group of supernovae is fairly homogeneous, with many of the blue spectral features well fit by Fe II. The nature of the spectroscopic commonalities and the variations in the class are discussed. Finally, we examine intermediates such as SN 2004eo and discuss the spectroscopic subgroup distribution of Type Ia supernovae.

preprint2011arXiv

Constraining the Properties of SNe~Ia Progenitors from Light Curves

We present an analysis of high precision V light curves (LC) for 18 local Type Ia Supernovae, SNe Ia, as obtained with the same telescope and setup at the Las Campanas Observatory (LCO). This homogeneity provides an intrinsic accuracy a few hundreds of a magnitude both with respect to individual LCs and between different objects. Based on the Single Degenerate Scenario, SD, we identify patterns which have been predicted by model calculations as signatures of the progenitor and accretion rate which change the explosion energy and the amount of electron capture, respectively. Using these templates as principle components and the overdetermined system of SN pairs, we reconstruct the properties of progenitors and progenitor systems. All LCO SNe Ia follow the brightness decline relation but 2001ay. After subtraction of the two components, the remaining scatter is reduced to 0.01-0.03m. Type SNe Ia seem to originate from progenitors with Main Sequence masses of 3Mo with the exception of two subluminous SNe Ia with < 2Mo. The component analysis indicates a wide range of accretion rates in the progenitor systems closing the gap to accretion induced collapses (AIC). SN1991t-like objects show differences in $dm15$ but no tracers of our secondary parameters. This may point to a different origin such as DD-Scenario or the Pulsating Delayed Detonations. SN2001ay does not follow the decline relation. It can be understood in the framework of C-rich WDs, and this group may produce an anti-Phillips relation. We suggest that this may be a result of a common envelope phase and mixing during central He burning as in SN1987A.

preprint2011arXiv

Theoretical light curves of type Ia supernovae

Aims. We present the first theoretical SN Ia light curves calculated with the time-dependent version of the general purpose model atmosphere code PHOENIX. Our goal is to produce light curves and spectra of hydro models of all types of supernovae. Methods. We extend our model atmosphere code PHOENIX to calculate type Ia supernovae light curves. A simple solver was imple- mented which keeps track of energy conservation in the atmosphere during the free expansion phase. Results. The correct operation of the new additions to PHOENIX were verified in test calculations. Furthermore, we calculated theo- retical light curves and compared them to the observed SN Ia light curves of SN 1999ee and SN 2002bo. We obtained LTE as well as NLTE model light curves. Conclusions. We have verified the correct operation of our extension into the time domain. We have calculated the first SN Ia model light curves using PHOENIX in both LTE and NLTE. For future work the infrared model light curves need to be further investigated.

preprint2010arXiv

Aspects of Multi-Dimensional Modelling of Substellar Atmospheres

Theoretical arguments and observations suggest that the atmospheres of Brown Dwarfs and planets are very dynamic on chemical and on physical time scales. The modelling of such substellar atmospheres has, hence, been much more demanding than initially anticipated. This Splinter (http://star-www.st-and.ac.uk/~ch80/CS16/MultiDSplinter_CS16.html) has combined new developments in atmosphere modelling, with novel observational techniques, and new challenges arising from planetary and space weather observations.

preprint2010arXiv

Hydrogen Recombination with Multilevel atoms

Hydrogen recombination is one of the most important atomic processes in many astrophysical objects such as Type II supernova (SN~II) atmospheres, the high redshift universe during the cosmological recombination era, and H II regions in the interstellar medium. Accurate predictions of the ionization fraction can be quite different from those given by a simple solution if one takes into account many angular momentum sub-states, non-resonant processes, and calculates the rates of all atomic processes from the solution of the radiative transfer equation instead of using a Planck function under the assumption of thermal equilibrium. We use the general purpose model atmosphere code PHOENIX 1D to compare how the fundamental probabilities such as the photo-ionization probability, the escape probability, and the collisional de-excitation probability are affected by the presence of other metals in the environment, multiple angular momentum sub-states, and non-resonant processes. Our comparisons are based on a model of SN 1999em, a SNe Type II, 20 days after its explosion.

preprint2010arXiv

Searching for Hydrogen in Type Ib Supernovae

We present synthetic spectral fits of the typical Type Ib SN 1999dn and the Hydrogen Rich Ib SN 2000H using the generalized non-local thermodynamic equilibrium stellar atmospheres code \phx. We fit model spectra to five epochs of SN 1999dn ranging from ten days pre-maximum light to 17 days post-maximum light and the two earliest epochs of SN 2000H available, maximum light and six days post-maximum. Our goal is to investigate the possibility of hydrogen in Type Ib Supernovae (SNe Ib), specifically a feature around 6200Å which has previously been attributed to high velocity H-alpha. In earlier work on SN 1999dn we found the most plausible alternative to H-alpha to be a blend of Si II and Fe II lines which can be adjusted to fit by increasing the metallicity. Our models are simple; they assume a powerlaw density profile with radius, homologous expansion, and solar compositions. The helium core is produced by burning 4H --> He in order to conserve nucleon number. For models with hydrogen the outer skin of the model consists of a shell of solar composition. The hydrogen mass of the standard solar composition shell is M_H less than about 0.001 times the mass of the sun in SN 1999dn and M_H less than about 0.2 times the mass of the sun for SN 2000H. Our models fit the observed spectra reasonably well, successfully reproducing most features including the characteristic He I absorptions. The hydrogen feature in SN 1999dn is clear, but much more pronounced in SN 2000H. We discuss a possible evolutionary scenario that accounts for the dichotomy in the hydrogen shell mass between these two supernovae.

preprint2009arXiv

3D Radiative Transfer with PHOENIX

Using the methods of general relativity Lindquist derived the radiative transfer equation that is correct to all orders in v/c. Mihalas developed a method of solution for the important case of monotonic velocity fields with spherically symmetry. We have developed the generalized atmosphere code PHOENIX, which in 1-D has used the framework of Mihalas to solve the radiative transfer equation (RTE) in 1-D moving flows. We describe our recent work including 3-D radiation transfer in PHOENIX and particularly including moving flows exactly using a novel affine method. We briefly discuss quantitative spectroscopy in supernovae.

preprint2009arXiv

A 3D radiative transfer framework: IV. spherical & cylindrical coordinate systems

We extend our framework for 3D radiative transfer calculations with a non-local operator splitting methods along (full) characteristics to spherical and cylindrical coordinate systems. These coordinate systems are better suited to a number of physical problems than Cartesian coordinates. The scattering problem for line transfer is solved via means of an operator splitting (OS) technique. The formal solution is based on a full characteristics method. The approximate $Λ$ operator is constructed considering nearest neighbors exactly. The code is parallelized over both wavelength and solid angle using the MPI library. We present the results of several test cases with different values of the thermalization parameter for the different coordinate systems. The results are directly compared to 1D plane parallel tests. The 3D results agree very well with the well-tested 1D calculations.

preprint2009arXiv

A 3D radiative transfer framework: V. Homologous Flows

Observations and theoretical calculations have shown the importance of non-spherically symmetric structures in supernovae. Thus, the interpretation of observed supernova spectra requires the ability to solve the transfer equation in 3-D moving atmospheres. We present an implementation of the solution of the radiative transfer equation in 3-D homologously expanding atmospheres in spherical coordinates. The implementation is exact to all orders in v/c. We use a new affine method that makes use of the fact that photons travel on straight lines. We compare our results in 3-D for spherically symmetric test problems with high velocity fields and find excellent agreement. Our well-tested 1-D results are based on methods where the momentum directions vary along the characteristic (co-moving momentum directions). Thus, we are able to verify both the analytic framework and its numerical implementation. Additionally, we have been able to test the parallelization over characteristics. Using 512^2 momentum angles we ran the code on 16,384 Opteron processors and achieved excellent scaling. It is now possible to calculate synthetic spectra from realistic 3D hydro simulations. This should open an era of progress in hydro modeling, similar to that that occurred in the 1980s when 1-D models were confronted with synthetic spectra.

preprint2009arXiv

A 3D radiative transfer framework: VI. PHOENIX/3D example applications

We demonstrate the application of our 3D radiative transfer framework in the model atmosphere code PHOENIX/3D for a number of spectrum synthesis calculations for very different conditions. The 3DRT framework discussed in the previous papers of this series was added to our general-purpose model atmosphere code PHOENIX/1D and an extended 3D version PHOENIX/3D was created. The \phxT code is parallelized via the MPI library using a hierarchical domain decomposition and displays very good strong scaling. We present the results of several test cases for widely different atmosphere conditions and compare the 3D calculations with equivalent 1D models to assess the internal accuracy of the 3D modeling. In addition, we show the results for a number of parameterized 3D structures. With presently available computational resources it is possible to solve the full 3D radiative transfer (including scattering) problem with the same micro-physics as included in 1D modeling.

preprint2009arXiv

On the Hydrogen Recombination Time in Type II Supernova Atmospheres

NLTE radiative transfer calculations of differentially expanding supernovae atmospheres are computationally intensive and are almost universally performed in time-independent snapshot mode. The validity of the steady-state approximation in the rate equations has recently been questioned. We calculate the effective recombination time of hydrogen in SN II using our general purpose model atmosphere code PHOENIX. While we find that the recombination time for the conditions of SNe II at early times is increased over the classical value for the case of a simple hydrogen model atom with energy levels corresponding to just the first 2 principle quantum numbers, the classical value of the recombination time is recovered in the case of a multi-level hydrogen atom. We also find that the recombination time at most optical depths is smaller in the case of a multi-level atom than for a simple two-level hydrogen atom. We find that time dependence in the rate equations is important in the early epochs of a supernova's lifetime. The changes due to the time dependent rate equation (at constant input luminosity) are manifested in physical parameters such as the level populations which directly affects the spectra. The H-alpha profile is affected by the time dependent rate equations at early times. At later times time dependence does not significantly modify the level populations and therefore the H-alpha profile is roughly independent of whether the steady-state or time-dependent approach is used.

preprint2009arXiv

Secondary Parameters of Type Ia Supernova Light Curves

High-quality observations of $B$ and $V$ light curves obtained at Las Campanas Observatory for local Type Ia Supernovae (SNe Ia) show clear evidence that SNe Ia with the same brightness decline or stretch may have systematic and independent deviations at times < 5 days before and at times > 30 days after maximum light. This suggests the existence of two independent secondary parameters which control the shape of SN Ia light curves in addition to the brightness decline relation. stretch. The differences are consistent in morphology of the time dependence and size with predictions by models within the delayed detonation scenario. The secondary parameters may reflect two independent physical effects caused by variations in the progenitor and accretion rates, and link the LC variations in shape with the intrinsic, absolute brightness.

preprint2009arXiv

Time-dependent radiative transfer with PHOENIX

Aims. We present first results and tests of a time-dependent extension to the general purpose model atmosphere code PHOENIX. We aim to produce light curves and spectra of hydro models for all types of supernovae. Methods. We extend our model atmosphere code PHOENIX to solve time-dependent non-grey, NLTE, radiative transfer in a special relativistic framework. A simple hydrodynamics solver was implemented to keep track of the energy conservation of the atmosphere during free expansion. Results. The correct operation of the new additions to PHOENIX were verified in test calculations. Conclusions. We have shown the correct operation of our extension to time-dependent radiative transfer and will be able to calculate supernova light curves and spectra in future work.

preprint2008arXiv

Stellar Population Models and Individual Element Abundances II: Stellar Spectra and Integrated Light Models

The first paper in this series explored the effects of altering the chemical mixture of the stellar population on an element by element basis on stellar evolutionary tracks and isochrones to the end of the red giant branch. This paper extends the discussion by incorporating the fully consistent synthetic stellar spectra with those isochrone models in predicting integrated colors, Lick indices, and synthetic spectra. Older populations display element ratio effects in their spectra at higher amplitude than younger populations. In addition, spectral effects in the photospheres of stars tend to dominate over effects from isochrone temperatures and lifetimes, but, further, the isochrone-based effects that are present tend to fall along the age-metallicity degeneracy vector, while the direct stellar spectral effects usually show considerable orthogonality.

preprint2007arXiv

The ACS Survey of Galactic Globular Clusters. II. Stellar Evolution Tracks, Isochrones, Luminosity Functions, and Synthetic Horizontal-Branch Models

The ACS Survey of Galactic Globular Clusters, an HST Treasury Project, will deliver high quality, homogeneous photometry of 65 globular clusters. This paper introduces a new collection of stellar evolution tracks and isochrones suitable for analyzing the ACS Survey data. Stellar evolution models were computed at [Fe/H]= -2.5, -2.0, -1.5, -1.0, -0.5, and 0; [alpha/Fe]= -0.2, 0, 0.2, 0.4, 0.6, and 0.8; and three initial He abundances for masses from 0.1 to 1.8 Msun and ages from 2 to 15 Gyr. Each isochrone spans a wide range in luminosity from Mv~14 up to the tip of the red giant branch. These are complemented by a set of He-burning tracks that extend from the zero age horizontal branch to the onset of thermal pulsations on the asymptotic giant branch. In addition, a set of computer programs are provided that make it possible to interpolate the isochrones in [Fe/H], generate luminosity functions from the isochrones, and create synthetic horizontal branch models. The tracks and isochrones have been converted to the observational plane with two different color-Teff transformations, one synthetic and one semi-empirical, in ground-based B, V, and I, and F606W and F814W for both ACS-WFC and WFPC2 systems. All models and programs presented in this paper are available from http://stellar.dartmouth.edu/~models/

preprint2006arXiv

Nearby Supernova Factory Observations of SN 2006D: On Sporadic Carbon Signatures in Early Type Ia Supernova Spectra

We present four spectra of the Type Ia supernova (SN Ia) 2006D extending from -7 to +13 days with respect to B-band maximum. The spectra include the strongest signature of unburned material at photospheric velocities observed in a SN Ia to date. The earliest spectrum exhibits C II absorption features below 14,000 km/s, including a distinctive C II λ6580 absorption feature. The carbon signatures dissipate as the SN approaches peak brightness. In addition to discussing implications of photospheric-velocity carbon for white dwarf explosion models, we outline some factors that may influence the frequency of its detection before and around peak brightness. Two effects are explored in this regard, including depopulation of the C II optical levels by non-LTE effects, and line-of-sight effects resulting from a clumpy distribution of unburned material with low volume-filling factor.

preprint2002arXiv

Detailed Spectroscopic Analysis of SN 1987A: The Distance to the LMC using the SEAM method

Supernova 1987A remains the most well-studied supernova to date. Observations produced excellent broad-band photometric and spectroscopic coverage over a wide wavelength range at all epochs. We model the observed spectra from Day 1 to Day 81 using a hydrodynamical model. We show that good agreement can be obtained at times up to about 60 days, if we allow for extended nickel mixing. Later than about 60 days the observed Balmer lines become stronger than our models can reproduce. We show that this is likely due to a more complicated distribution of gamma-rays than we allow for in our spherically symmetric calculations. We present synthetic light curves in UBVRIJHK and a synthetic bolometric light curve. Using this broad baseline of detailed spectroscopic models we find a distance modulus mu = 18.5 +/- 0.2 using the SEAM method of determining distances to supernovae. We find that the explosion time agrees with that of the neutrino burst and is constrained at 68% confidence to within +/- 0.9 days. We argue that the weak Balmer lines of our detailed model calculations casts doubt on the accuracy of the purely photometric EPM method. We also suggest that Type IIP supernovae will be most useful as distance indicators at early times due to a variety of effects.

preprint2002arXiv

Detectibility of Hydrogen Mixing in Type Ia Supernova pre-Maximum Spectra

The presence of a small amount of hydrogen is expected in most single degenerate scenarios for producing a Type Ia supernova (SN Ia). While hydrogen may be detected in very early high resolution optical spectra, in early radio spectra, and in X-ray spectra, here we examine the possibility of detecting hydrogen in early low resolution spectra such as those that will be obtained by proposed large scale searches for nearby SNe Ia. We find that definitive detections will require both very early spectra (less than 5 days after explosion) and perhaps slightly higher amounts of hydrogen than are currently predicted to be mixed into the outer layers of SNe Ia. Thus, the non-detection of hydrogen so far does not in and of itself rule out any current progenitor models. Nevertheless, very early spectra of SNe Ia will provide significant clues to the amount of hydrogen present and hence to the nature of the SN Ia progenitor system. Spectral coverage in both the optical and IR will be required to definitively identify hydrogen in low resolution spectra.

preprint2000arXiv

Preliminary Spectral Analysis of the Type II Supernova 1999em

We have calculated fast direct spectral model fits to two early-time spectra of the Type-II plateau SN 1999em, using the SYNOW synthetic spectrum code. The first is an extremely early blue optical spectrum and the second a combined HST and optical spectrum obtained one week later. Spectroscopically this supernova appears to be a normal Type II and these fits are in excellent agreement with the observed spectra. Our direct analysis suggests the presence of enhanced nitrogen. We have further studied these spectra with the full NLTE general model atmosphere code PHOENIX. While we do not find confirmation for enhanced nitrogen (nor do we rule it out), we do require enhanced helium. An even more intriguing possible line identification is complicated Balmer and He I lines, which we show falls naturally out of the detailed calculations with a shallow density gradient. We also show that very early spectra such as those presented here combined with sophisticated spectral modeling allows an independent estimate of the total reddening to the supernova, since when the spectrum is very blue, dereddening leads to changes in the blue flux that cannot be reproduced by altering the ``temperature'' of the emitted radiation. These results are extremely encouraging since they imply that detailed modeling of early spectra can shed light on both the abundances and total extinction of SNe II, the latter improving their utility and reliability as distance indicators.

preprint1998arXiv

Numerical Solution of the Expanding Stellar Atmosphere Problem

In this paper we discuss numerical methods and algorithms for the solution of NLTE stellar atmosphere problems involving expanding atmospheres, e.g., found in novae, supernovae and stellar winds. We show how a scheme of nested iterations can be used to reduce the high dimension of the problem to a number of problems with smaller dimensions. As examples of these sub-problems, we discuss the numerical solution of the radiative transfer equation for relativistically expanding media with spherical symmetry, the solution of the multi-level non-LTE statistical equilibrium problem for extremely large model atoms, and our temperature correction procedure. Although modern iteration schemes are very efficient, parallel algorithms are essential in making large scale calculations feasible, therefore we discuss some parallelization schemes that we have developed.

preprint1997arXiv

NLTE effects of Ti~I in M dwarfs and giants

We present detailed NLTE Ti~I calculations in model atmospheres of cool dwarf and giant stars. A fully self-consistent NLTE treatment for a Ti~I model atom with 395 levels and 5279 primary bound-bound transitions is included, and we discuss the implication of departures from LTE in this atom for the strengths of Ti~I lines and TiO molecular bands in cool star spectra. We show that in the atmospheric parameter range investigated, LTE is a poor approximation to Ti~I line formation, as expected from the low collisional rates in cool stars. The secondary effects of Ti~I overionization on the TiO number density and the TiO molecular opacities, however, are found to be negligible in the molecular line forming region for the relatively small parameter range studied in this paper.

preprint1996arXiv

NLTE Model Atmospheres for M Dwarfs and Giants

The atmospheres of M stars are dominated by a small number of very strong molecular compounds (H$_2$O, TiO, H$_2$, CO, VO). Most of the hydrogen is locked in molecular H$_2$, most of the carbon in CO; and H$_2$O, TiO and VO opacities define a pseudo-continuum covering the entire flux distribution of these stars. The optical ``continuum'' is due to TiO vibrational bands which are often used as temperature indicators for these stars. These may be the depth of the bands relative to the troughs in between them; or the depth of the VO bands; or of the atomic lines relative to the local ``continuum''; or even the strength of the infrared water bands; all of these depend on the strength of the TiO bands and the amount of flux-redistribution to longer wavelengths exerted by them. Departures from LTE of the Ti I atom, and thus the concentration of the important TiO molecule, could, therefore, have severe and measurable consequences on the atmospheric structure and spectra of these stars. In this paper we discuss NLTE effects of Ti I in fully self-consistent models for a few representative M/Brown dwarf and M giant model atmospheres and spectra.

preprint1996arXiv

Non-LTE Treatment of Fe II in Astrophysical Plasmas

We describe our implementation of an extremely detailed model atom of singly ionized iron for NLTE computations in static and moving astrophysical plasmas. Our model atom includes 617 levels, 13675 primary permitted transitions and up to 1.2 million secondary transitions. Our approach guarantees that the total iron opacity is included at the correct wavelength with reasonable memory and CPU requirements. We find that the lines saturate the wavelength space, such that special wavelength points inserted along the detailed profile functions may be replaced with a statistical sampling method. We describe the results of various test calculations for novae and supernovae.

preprint1996arXiv

Synthetic Spectra of Hydrodynamic Models of Type Ia Supernovae

We present detailed NLTE synthetic spectra of hydrodynamic SNe Ia models. We make no assumptions about the form of the spectrum at the inner boundary. We calculate both Chandrasekhar-mass deflagration models and sub-Chandrasekhar ``helium detonators.'' Gamma-ray deposition is handled in a simple, accurate manner. We have parameterized the storage of energy that arises from the time dependent deposition of radioactive decay energy in a reasonable manner, that spans the expected range. We find that the Chandrasekhar-mass deflagration model W7 of Nomoto etal shows good agreement with the observed spectra of SN 1992A and SN 1994D, particularly in the UV, where our models are expected to be most accurate. The sub-Chandrasekhar models do not reproduce the UV deficit observed in normal SNe Ia. They do bear some resemblance to sub-luminous SNe Ia, but the shape of the spectra (i.e. the colors) are opposite to that of the observed ones and the intermediate mass element lines such as Si II, and Ca II are extremely weak, which seems to be a generic difficulty of the models. Although the sub-Chandrasekhar models have a significant helium abundance (unlike Chandrasekhar-mass models), helium lines are not prominent in the spectra near maximum light and thus do not act as a spectral signature for the progenitor.

preprint1996arXiv

The Effects of Fe II NLTE on Nova Atmospheres and Spectra

The atmospheres of novae at early times in their outbursts are very extended, expanding shells with low densities. Models of these atmospheres show that NLTE effects are very important and must be included in realistic calculations. We have, therefore, been improving our atmospheric studies by increasing the number of ions treated in NLTE. One of the most important ions is Fe II which has a complex structure and numerous lines in the observable spectrum. In this paper we investigate NLTE effects for Fe II for a wide variety of parameters. We use a detailed Fe II model atom with 617 level and 13675 primary lines, treated using a rate-operator formalism. We show that the radiative transfer equation in nova atmospheres {\em must} be treated with sophisticated numerical methods and that simple approximations, such as the Sobolev method, {\em cannot} be used because of the large number of overlapping lines in the co-moving frame. Our results show that the formation of the Fe II lines is strongly affected by NLTE effects. For low effective temperatures, $\Teff < 20,000\,$K, the optical Fe II lines are most influenced by NLTE effects, while for higher $\Teff$ the UV lines of Fe II are very strongly affected by NLTE. The departure coefficients are such that Fe II tends to be overionized in NLTE when compared to LTE. Therefore, Fe II-NLTE must be included with sophisticated radiative transfer in nova atmosphere models in order to reliably analyze observed nova spectra.

preprint1995arXiv

Preliminary Spectral Analysis of SN 1994I

We present optical spectra of the Type Ic supernova 1994I in M51 and preliminary non-LTE analysis of the spectra. Our models are not inconsistent with the explosions of C+O cores of massive stars. While we find no direct evidence for helium in the optical spectra, our models cannot rule out small amounts of helium. More than 0.1~\msol\ of helium seems unlikely.

preprint1995arXiv

Statistical Connections Between the Properties of Type Ia Supernovae and the B--V Colors of Their Parent Galaxies, and the Value of H$_0$

Statistical connections between the properties of Type Ia supernovae (SNe Ia) and the B--V colors of their parent galaxies are established. Compared to SNe Ia in blue galaxies [B-V$\la$0.75], SNe Ia in redder galaxies have (1) a wider dispersion in the blueshifts of their Si II $λ$6355 absorption features, ten days after maximum light; (2) more rapidly declining light curves; and (3) lower luminosities. Even when the spectroscopically peculiar, very subluminous SNe Ia such as SN 1991bg are disregarded, SNe Ia in red galaxies are less luminous than those in blue galaxies by about 0.3 magnitudes. When SNe Ia that are thought to have been significantly extinguished by dust in their parent galaxies are disregarded, those in blue galaxies have observational absolute--magnitude dispersions of only $σ_{obs}(M_B)$=0.20 and $σ_{obs}(M_V)$=0.17, which implies that their intrinsic absolute--magnitude dispersions are very small. We use six SNe Ia whose absolute magnitudes have been calibrated by means of Cepheids, which also indicate that the intrinsic absolute--magnitude dispersions of SNe Ia in blue galaxies are very small, to calibrate SNe Ia in blue galaxies and obtain $\rm H_0=57\pm4\ km\ s^{-1}\ Mpc^{-1}$. This value is in excellent agreement with that obtained by Saha et al. (1995b), in spite of the fact that they do not take into account any dependence of SN Ia absolute magnitude on the nature of the parent galaxy. Some implications of the statistical connections between SNe Ia and the colors of their parent galaxies, for identifying the progenitor binary systems of SNe Ia and for using high--redshift SNe Ia to measure q$_0$, are briefly discussed.

E. Baron

What is connected

Connect this record

See the researcher in context

Building this map preview

60 published item(s)

A Speed Bump: SN 2021aefx Shows that Doppler Shift Alone can Explain Early-Excess Blue Flux in Some Type Ia Supernovae

A Tale of Two Type Ia Supernovae: The fast-declining siblings SNe 2015bo and 1997cn

SN 2019va: A Type IIP Supernova with Large Influence of Nickel-56 Decay on the Plateau-phase Light Curve

A 3D radiative transfer framework: XII. Many-core, vector and GPU methods

Carnegie Supernova Project-II: Near-infrared Spectroscopy of Stripped-Envelope Core-Collapse Supernovae

Carnegie Supernova Project II: The slowest rising Type Ia supernova LSQ14fmg and clues to the origin of super-Chandrasekhar/03fg-like events

Carnegie Supernova Project-II: A new method to photometrically identify sub-types of extreme Type Ia Supernovae

Discovery and Rapid Follow-up Observations of the Unusual Type II SN 2018ivc in NGC 1068

SN 2018hti: a nearby superluminous supernova discovered in a metal-poor galaxy

The Carnegie Supernova Project II. Observations of the intermediate luminosity red transient SNhunt120

The Carnegie Supernova Project II. Observations of the luminous red nova AT 2014ej

Predicting the Extreme Ultraviolet Radiation Environment of Exoplanets Around Low-Mass Stars: GJ 832, GJ 176, GJ 436

The Carnegie Supernova Project II. Early observations and progenitor constraints of the Type Ib supernova LSQ13abf

Supernova 2014J at M82: II. Direct Analysis of Spectra Obtained with Isaac Newton and William Herschel Telescopes

Identification of the feature that causes the I-band secondary maximum of a type Ia supernova

Spectral Models for Early Time SN 2011fe Observations

Theoretical Clues to the Ultraviolet Diversity of Type Ia Supernovae

A 3D radiative transfer framework: XI. multi-level NLTE

Comprehensive Observations of the Bright and Energetic Type Iax SN 2012Z: Interpretation as a Chandrasekhar Mass White Dwarf Explosion

Near-infrared line identification in type Ia supernovae during the transitional phase

SNe Ia: Can Chandrasekhar Mass Explosions Reproduce the Observed Zoo?

Hydrogen and helium in the spectra of Type Ia supernovae

Optical and Near-IR Observations of the Faint and Fast 2008ha-like Supernova 2010ae

A 3D radiative transfer framework IX. Time dependence

A 3D radiative transfer framework: X. Arbitrary Velocity Fields in the Co-moving Frame

A Physical Model for SN 2001ay, a normal, bright, extremely slowly declining Type Ia supernova

Analysis of the Early-Time Optical Spectra of SN 2011fe in M101

Near-infrared light curves of type Ia supernovae

Quantitative Photospheric Spectral Analysis of the Type IIP Supernova 2007od

Supernova Resonance--scattering Line Profiles in the Absence of a Photosphere

The bright Type IIP SN 2009bw, showing signs of interaction

The UV/optical spectra of the Type Ia supernova SN 2010jn: a bright supernova with outer layers rich in iron-group elements

A 3D radiative transfer framework: XIII. OpenCL implementation

Comparative Analysis of Peculiar Type Ia 1991bg-like Supernovae Spectra

Constraining the Properties of SNe~Ia Progenitors from Light Curves

Theoretical light curves of type Ia supernovae

Aspects of Multi-Dimensional Modelling of Substellar Atmospheres

Hydrogen Recombination with Multilevel atoms

Searching for Hydrogen in Type Ib Supernovae

3D Radiative Transfer with PHOENIX

A 3D radiative transfer framework: IV. spherical & cylindrical coordinate systems

A 3D radiative transfer framework: V. Homologous Flows

A 3D radiative transfer framework: VI. PHOENIX/3D example applications

On the Hydrogen Recombination Time in Type II Supernova Atmospheres

Secondary Parameters of Type Ia Supernova Light Curves

Time-dependent radiative transfer with PHOENIX

Stellar Population Models and Individual Element Abundances II: Stellar Spectra and Integrated Light Models

The ACS Survey of Galactic Globular Clusters. II. Stellar Evolution Tracks, Isochrones, Luminosity Functions, and Synthetic Horizontal-Branch Models

Nearby Supernova Factory Observations of SN 2006D: On Sporadic Carbon Signatures in Early Type Ia Supernova Spectra

Detailed Spectroscopic Analysis of SN 1987A: The Distance to the LMC using the SEAM method

Detectibility of Hydrogen Mixing in Type Ia Supernova pre-Maximum Spectra

Preliminary Spectral Analysis of the Type II Supernova 1999em

Numerical Solution of the Expanding Stellar Atmosphere Problem

NLTE effects of Ti~I in M dwarfs and giants

NLTE Model Atmospheres for M Dwarfs and Giants

Non-LTE Treatment of Fe II in Astrophysical Plasmas

Synthetic Spectra of Hydrodynamic Models of Type Ia Supernovae

The Effects of Fe II NLTE on Nova Atmospheres and Spectra

Preliminary Spectral Analysis of SN 1994I

Statistical Connections Between the Properties of Type Ia Supernovae and the B--V Colors of Their Parent Galaxies, and the Value of H$_0$