Showing 1–28 of 28 results for author: Shingles, L

ZTF SN Ia DR2: The secondary maximum in Type Ia supernovae

Authors: M. Deckers, K. Maguire, L. Shingles, G. Dimitriadis, M. Rigault, M. Smith, A. Goobar, J. Nordin, J. Johansson, M. Amenouche, U. Burgaz, S. Dhawan, M. Ginolin, L. Harvey, W. D. Kenworthy, Y. -L. Kim, R. R. Laher, N. Luo, S. R. Kulkarni, F. J. Masci, T. E. Müller-Bravo, P. E. Nugent, N. Pletskova, J. Purdum, B. Racine , et al. (2 additional authors not shown)

Abstract: Type Ia supernova (SN Ia) light curves have a secondary maximum that exists in the $r$, $i$, and near-infrared filters. The secondary maximum is relatively weak in the $r$ band, but holds the advantage that it is accessible, even at high redshift. We used Gaussian Process fitting to parameterise the light curves of 893 SNe Ia from the Zwicky Transient Facility's (ZTF) second data release (DR2), an… ▽ More Type Ia supernova (SN Ia) light curves have a secondary maximum that exists in the $r$, $i$, and near-infrared filters. The secondary maximum is relatively weak in the $r$ band, but holds the advantage that it is accessible, even at high redshift. We used Gaussian Process fitting to parameterise the light curves of 893 SNe Ia from the Zwicky Transient Facility's (ZTF) second data release (DR2), and we were able to extract information about the timing and strength of the secondary maximum. We found $>5σ$ correlations between the light curve decline rate ($Δm_{15}(g)$) and the timing and strength of the secondary maximum in the $r$ band. Whilst the timing of the secondary maximum in the $i$ band also correlates with $Δm_{15}(g)$, the strength of the secondary maximum in the $i$ band shows significant scatter as a function of $Δm_{15}(g)$. We found that the transparency timescales of 97 per cent of our sample are consistent with double detonation models, and that SNe Ia with small transparency timescales ($<$ 32 d) reside predominantly in locally red environments. We measured the total ejected mass for the normal SNe Ia in our sample using two methods, and both were consistent with medians of $1.3\ \pm \ 0.3$ and $1.2\ \pm\ 0.2$ solar masses. We find that the strength of the secondary maximum is a better standardisation parameter than the SALT light curve stretch ($x_1$). Finally, we identified a spectral feature in the $r$ band as Fe II, which strengthens during the onset of the secondary maximum. The same feature begins to strengthen at $<$ 3 d post maximum light in 91bg-like SNe. Finally, the correlation between $x_1$ and the strength of the secondary maximum was best fit with a broken line, with a split at $x_1^0\ =\ -0.5\ \pm\ 0.2$, suggestive of the existence of two populations of SNe Ia. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 18 pages, 14 figures

Primary and secondary source of energy in the superluminous supernova 2018ibb

Authors: Alexandra Kozyreva, Luke Shingles, Petr Baklanov, Alexey Mironov, Fabian R. N. Schneider

Abstract: We examine the pair-instability origin of superluminous supernova 2018ibb. As the base model, we use a non-rotating stellar model with an initial mass of 250 Msun at about 1/15 solar metallicity. We consider three versions of the model as input for radiative transfer simulations done with the STELLA and ARTIS codes: with 25 Msun of 56Ni, 34 Msun of 56Ni, and a chemically mixed case with 34 Msun of… ▽ More We examine the pair-instability origin of superluminous supernova 2018ibb. As the base model, we use a non-rotating stellar model with an initial mass of 250 Msun at about 1/15 solar metallicity. We consider three versions of the model as input for radiative transfer simulations done with the STELLA and ARTIS codes: with 25 Msun of 56Ni, 34 Msun of 56Ni, and a chemically mixed case with 34 Msun of 56Ni. We present light curves and spectra in comparison to the observed data of SN 2018ibb, and conclude that the pair-instability supernova model with 34 Msun of 56Ni explains broad-band light curves reasonably well between -100 and 250 days around the peak. Our synthetic spectra have many similarities with the observed spectra. The luminosity excess in the light curves and the blue-flux excess in the spectra can be explained by an additional energy source, which may be interaction of the SN ejecta with circumstellar matter. We discuss possible mechanisms of the origin of the circumstellar matter being ejected in the decades before the pair-instability explosion. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 7 pages, 4 figures, accepted for publication in A & A

Including a Luminous Central Remnant in Radiative Transfer Simulations for Type Iax Supernovae

Authors: F. P. Callan, S. A. Sim, C. E. Collins, L. J. Shingles, F. Lach, F. K. Roepke, R. Pakmor, M. Kromer, S. Srivastav

Abstract: Type Iax supernovae (SNe Iax) are proposed to arise from deflagrations of Chandrasekhar mass white dwarfs (WDs). Previous deflagration simulations have achieved good agreement with the light curves and spectra of intermediate-luminosity and bright SNe Iax. However, the model light curves decline too quickly after peak, particularly in red optical and near-infrared (NIR) bands. Deflagration models… ▽ More Type Iax supernovae (SNe Iax) are proposed to arise from deflagrations of Chandrasekhar mass white dwarfs (WDs). Previous deflagration simulations have achieved good agreement with the light curves and spectra of intermediate-luminosity and bright SNe Iax. However, the model light curves decline too quickly after peak, particularly in red optical and near-infrared (NIR) bands. Deflagration models with a variety of ignition configurations do not fully unbind the WD, leaving a remnant polluted with $^{56}\mathrm{Ni}$. Emission from such a remnant may contribute to the luminosity of SNe Iax. Here we investigate the impact of adding a central energy source, assuming instantaneous powering by $^{56}\mathrm{Ni}$ decay in the remnant, in radiative transfer calculations of deflagration models. Including the remnant contribution improves agreement with the light curves of SNe Iax, particularly due to the slower post-maximum decline of the models. Spectroscopic agreement is also improved, with intermediate-luminosity and faint models showing greatest improvement. We adopt the full remnant $^{56}\mathrm{Ni}$ mass predicted for bright models, but good agreement with intermediate-luminosity and faint SNe Iax is only possible for remnant $^{56}\mathrm{Ni}$ masses significantly lower than those predicted. This may indicate that some of the $^{56}\mathrm{Ni}$ decay energy in the remnant does not contribute to the radiative luminosity but instead drives mass ejection, or that escape of energy from the remnant is significantly delayed. Future work should investigate the structure of remnants predicted by deflagration models and the potential roles of winds and delayed energy escape, as well as extend radiative transfer simulations to late times. △ Less

Submitted 19 April, 2024; v1 submitted 22 March, 2024; originally announced March 2024.

Comments: 17 pages, 6 figures. Lightcurves and spectra available at https://hesma.h-its.org

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 530, Issue 2, May 2024, Pages 1457 to 1473

GW190425: Pan-STARRS and ATLAS coverage of the skymap and limits on optical emission associated with FRB190425

Authors: S. J. Smartt, M. Nicholl, S. Srivastav, M. E. Huber, K. C. Chambers, K. W. Smith, D. R. Young, M. D. Fulton, J. L. Tonry, C. W. Stubbs, L. Denneau, A. J. Cooper, A. Aamer, J. P. Anderson, A. Andersson, J. Bulger, T. -W Chen, P. Clark, T. de Boer, H. Gao, J. H. Gillanders, A. Lawrence, C. C. Lin, T. B. Lowe, E. A. Magnier , et al. (10 additional authors not shown)

Abstract: GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverag… ▽ More GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverage of the skymap beginning within 1 hour and 3 hours respectively of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 190425 has been suggested, given their spatial and temporal coincidence. The smaller sky localisation area of FRB 190425 and its dispersion measure have led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 +/- 10 Mpc. Our optical imaging covered the galaxy 6.0 hrs after GW190425 was detected and 3.5 hrs after the FRB 190425. No optical emission was detected and further imaging at +1.2 and +13.2 days also revealed no emission. If the FRB 190425 and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 190425 from a BNS merger involves a supramassive magnetised neutron star spinning down by dipole emission on the timescale of hours. We show that magnetar enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB-GW link for this source. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS, 20th Sept 2023, 9 pages

Towards inferring the geometry of kilonovae

Authors: Christine E. Collins, Luke J. Shingles, Andreas Bauswein, Stuart A. Sim, Theodoros Soultanis, Vimal Vijayan, Andreas Floers, Oliver Just, Gerrit Leck, Georgios Lioutas, Gabriel Martínez-Pinedo, Albert Sneppen, Darach Watson, Zewei Xiong

Abstract: Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hy… ▽ More Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4-7 per cent). △ Less

Submitted 23 February, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted by MNRAS

Helium as a signature of the double detonation in Type Ia supernovae

Authors: Christine E. Collins, Stuart A. Sim, Luke. J. Shingles, Sabrina Gronow, Friedrich K. Roepke, Ruediger Pakmor, Ivo R. Seitenzahl, Markus Kromer

Abstract: The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer… ▽ More The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer simulations have not been able to fully address whether helium spectral features would form. This is because helium can not be sufficiently excited to form spectral features by thermal processes, but can be excited by collisions with non-thermal electrons, which most studies have neglected. We carry out a full non-local thermodynamic equilibrium (non-LTE) radiative transfer simulation for an instance of a double detonation explosion model, and include a non-thermal treatment of fast electrons. We find a clear He I λ 10830 feature which is strongest in the first few days after explosion and becomes weaker with time. Initially this feature is blended with the Mg II λ 10927 feature but over time separates to form a secondary feature to the blue wing of the Mg II λ 10927 feature. We compare our simulation to observations of iPTF13ebh, which showed a similar feature to the blue wing of the Mg II λ 10927 feature, previously identified as C I. Our simulation shows a good match to the evolution of this feature and we identify it as high velocity He I λ 10830. This suggests that He I λ 10830 could be a signature of the double detonation scenario. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 7 pages, accepted by MNRAS

Self-consistent 3D radiative transfer for kilonovae: directional spectra from merger simulations

Authors: Luke J. Shingles, Christine E. Collins, Vimal Vijayan, Andreas Flörs, Oliver Just, Gerrit Leck, Zewei Xiong, Andreas Bauswein, Gabriel Martínez-Pinedo, Stuart A. Sim

Abstract: We present three-dimensional radiative transfer calculations for the ejecta from a neutron star merger that include line-by-line opacities for tens of millions of bound-bound transitions, composition from an r-process nuclear network, and time-dependent thermalization of decay products from individual $α$ and $β^-$ decay reactions. In contrast to expansion opacities and other wavelength-binned tre… ▽ More We present three-dimensional radiative transfer calculations for the ejecta from a neutron star merger that include line-by-line opacities for tens of millions of bound-bound transitions, composition from an r-process nuclear network, and time-dependent thermalization of decay products from individual $α$ and $β^-$ decay reactions. In contrast to expansion opacities and other wavelength-binned treatments, a line-by-line treatment enables us include fluorescence effects and associate spectral features with the emitting and absorbing lines of individual elements. We find variations in the synthetic observables with both the polar and azimuthal viewing angles. The spectra exhibit blended features with strong interactions by Ce III, Sr II, Y II, and Zr II that vary with time and viewing direction. We demonstrate the importance of wavelength-calibration of atomic data using a model with calibrated Sr, Y, and Zr data, and find major differences in the resulting spectra, including a better agreement with AT2017gfo. The synthetic spectra for near-polar inclination show a feature at around 8000 A, similar to AT2017gfo. However, they evolve on a more rapid timescale, likely due to the low ejecta mass (0.005 M$_\odot$) as we take into account only the early ejecta. The comparatively featureless spectra for equatorial observers gives a tentative prediction that future observations of edge-on kilonovae will appear substantially different from AT2017gfo. We also show that 1D models obtained by spherically averaging the 3D ejecta lead to dramatically different direction-integrated luminosities and spectra compared to full 3D calculations. △ Less

Submitted 1 September, 2023; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: 12 pages, 5 figures. Accepted by ApJL

Journal ref: 2023, ApJL, 954, L41

1100 days in the life of the supernova 2018ibb -- The best pair-instability supernova candidate, to date

Authors: Steve Schulze, Claes Fransson, Alexandra Kozyreva, Ting-Wan Chen, Ofer Yaron, Anders Jerkstrand, Avishay Gal-Yam, Jesper Sollerman, Lin Yan, Tuomas Kangas, Giorgos Leloudas, Conor M. B. Omand, Stephen J. Smartt, Yi Yang, Matt Nicholl, Nikhil Sarin, Yuhan Yao, Thomas G. Brink, Amir Sharon, Andrea Rossi, Ping Chen, Zhihao Chen, Aleksandar Cikota, Kishalay De, Andrew J. Drake , et al. (41 additional authors not shown)

Abstract: Abridged - Stars with ZAMS masses between 140 and $260 M_\odot$ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN2018ibb is a H-poor SLS… ▽ More Abridged - Stars with ZAMS masses between 140 and $260 M_\odot$ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN2018ibb is a H-poor SLSN at $z=0.166$ that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the NIR with 2-10m class telescopes. SN2018ibb radiated $>3\times10^{51} \rm erg$ during its evolution, and its bolometric light curve reached $>2\times10^{44} \rm erg\,s^{-1}$ at peak. The long-lasting rise of $>93$ rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source ($^{56}$Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions for their photometric and spectroscopic properties. SN2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25-44 $M_\odot$ of freshly nucleosynthesised $^{56}$Ni, pointing to the explosion of a metal-poor star with a He-core mass of 120-130 $M_\odot$ at the time of death. This interpretation is also supported by the tentative detection of [Co II]$λ$1.025$μ$m, which has never been observed in any other PISN candidate or SLSN before. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN2018ibb by far the best candidate for being a PISN, to date. △ Less

Submitted 24 November, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: Accepted in A&A, the revised version includes a PISN rate estimate and an additional test with PISN models. 47 pages, main text 41 pages, 38 figures, 16 Tables

Photometric study of the late-time near-infrared plateau in Type Ia supernovae

Authors: M. Deckers, O. Graur, K. Maguire, L. Shingles, S. J. Brennan, J. P. Anderson, J. Burke, T. -W. Chen, L. Galbany, M. J. P. Grayling, C. P. Gutiérrez, L. Harvey, D. Hiramatsu, D. A. Howell, C. Inserra, T. Killestein, C. McCully, T. E. Müller-Bravo, M. Nicholl, M. Newsome, E. Padilla Gonzalez, C. Pellegrino, G. Terreran, J. H. Terwel, M. Toy , et al. (1 additional authors not shown)

Abstract: We present an in-depth study of the late-time near-infrared plateau in Type Ia supernovae (SNe Ia), which occurs between 70-500 d. We double the existing sample of SNe Ia observed during the late-time near-infrared plateau with new observations taken with the Hubble Space Telescope, Gemini, New Technology Telescope, the 3.5m Calar Alto Telescope, and the Nordic Optical Telescope. Our sample consis… ▽ More We present an in-depth study of the late-time near-infrared plateau in Type Ia supernovae (SNe Ia), which occurs between 70-500 d. We double the existing sample of SNe Ia observed during the late-time near-infrared plateau with new observations taken with the Hubble Space Telescope, Gemini, New Technology Telescope, the 3.5m Calar Alto Telescope, and the Nordic Optical Telescope. Our sample consists of 24 nearby SNe Ia at redshift < 0.025. We are able to confirm that no plateau exists in the Ks band for most normal SNe Ia. SNe Ia with broader optical light curves at peak tend to have a higher average brightness on the plateau in J and H, most likely due to a shallower decline in the preceding 100 d. SNe Ia that are more luminous at peak also show a steeper decline during the plateau phase in H. We compare our data to state-of-the-art radiative transfer models of nebular SNe Ia in the near-infrared. We find good agreement with the sub-Mch model that has reduced non-thermal ionisation rates, but no physical justification for reducing these rates has yet been proposed. An analysis of the spectral evolution during the plateau demonstrates that the ratio of [Fe II] to [Fe III] contribution in a near-infrared filter determines the light curve evolution in said filter. We find that overluminous SNe decline slower during the plateau than expected from the trend seen for normal SNe Ia △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 17 pages, 8 figures, Accepted for publication in MNRAS

Phase Curves of Kuiper Belt Objects, Centaurs, and Jupiter Family Comets from the ATLAS Survey

Authors: Matthew M. Dobson, Megan E. Schwamb, Susan D. Benecchi, Anne J. Verbiscer, Alan Fitzsimmons, Luke J. Shingles, Larry Denneau, A. N. Heinze, Ken W. Smith, John L. Tonry, Henry Weiland, David. R. Young

Abstract: The Kuiper belt objects, the Centaurs, and the Jupiter-family comets form an evolutionary continuum of small outer Solar System objects, and their study allows us to gain insight into the history and evolution of the Solar System. Broadband photometry can be used to measure their phase curves, allowing a first-order probe into the surface properties of these objects, though limited telescope time… ▽ More The Kuiper belt objects, the Centaurs, and the Jupiter-family comets form an evolutionary continuum of small outer Solar System objects, and their study allows us to gain insight into the history and evolution of the Solar System. Broadband photometry can be used to measure their phase curves, allowing a first-order probe into the surface properties of these objects, though limited telescope time makes measuring accurate phase curves difficult. We make use of serendipitous broadband photometry from the long-baseline, high-cadence ATLAS survey to measure the phase curves for a sample of 18 Kuiper belt objects, Centaurs, and Jupiter-family comets with unprecedentedly large datasets. We find phase curves with previously reported negative slopes become positive with increased data and are thus due to insufficient sampling of the phase curve profile, and not a real physical effect. We search for correlations between phase curve parameters, finding no strong correlations between any parameter pair, consistent with the findings of previous studies. We search for instances of cometary activity in our sample, finding a previously reported outburst by Echeclus and a new epoch of increased activity by Chiron. Applying the main belt asteroid HG1G2 phase curve model to three Jupiter-family comets in our sample with large phase angle spans, we find their slope parameters imply surfaces more consistent with those of carbonaceous main belt asteroids than silicaceous ones. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: 42 pages, 24 figures, accepted for publication to the Planetary Science Journal. Supplementary materials will be available in the online journal version

The Birth of a Relativistic Jet Following the Disruption of a Star by a Cosmological Black Hole

Authors: Dheeraj R. Pasham, Matteo Lucchini, Tanmoy Laskar, Benjamin P. Gompertz, Shubham Srivastav, Matt Nicholl, Stephen J. Smartt, James C. A. Miller-Jones, Kate D. Alexander, Rob Fender, Graham P. Smith, Michael D. Fulton, Gulab Dewangan, Keith Gendreau, Eric R. Coughlin, Lauren Rhodes, Assaf Horesh, Sjoert van Velzen, Itai Sfaradi, Muryel Guolo, N. Castro Segura, Aysha Aamer, Joseph P. Anderson, Iair Arcavi, Sean J. Brennan , et al. (41 additional authors not shown)

Abstract: A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to under… ▽ More A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to understand the radiative mechanisms operating in super-Eddington jets. Here, we present multi-wavelength (X-ray, UV, optical, and radio) observations of the optically discovered transient \target at $z=1.193$. Its unusual X-ray properties, including a peak observed luminosity of $\gtrsim$10$^{48}$ erg s$^{-1}$, systematic variability on timescales as short as 1000 seconds, and overall duration lasting more than 30 days in the rest-frame are traits associated with relativistic TDEs. The X-ray to radio spectral energy distributions spanning 5-50 days after discovery can be explained as synchrotron emission from a relativistic jet (radio), synchrotron self-Compton (X-rays), and thermal emission similar to that seen in low-redshift TDEs (UV/optical). Our modeling implies a beamed, highly relativistic jet akin to blazars but requires extreme matter-domination, i.e, high ratio of electron-to-magnetic field energy densities in the jet, and challenges our theoretical understanding of jets. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: To appear in Nature Astronomy on 30th November 2022. Also see here for an animation explaining the result: https://youtu.be/MQHdSbxuznY

The luminous type Ia supernova 2022ilv and its early excess emission

Authors: Shubham Srivastav, S. J. Smartt, M. E. Huber, G. Dimitriadis, K. C. Chambers, Michael D. Fulton, Thomas Moore, F. P. Callan, James H. Gillanders, K. Maguire, M. Nicholl, Luke J. Shingles, S. A. Sim, K. W. Smith, J. P. Anderson, Thomas de Boer, Ting-Wan Chen, Hua Gao, D. R. Young

Abstract: We present observations and analysis of the host-less and luminous type Ia supernova 2022ilv, illustrating it is part of the 2003fg-like family, often referred to as super-Chandrasekhar (Ia-SC) explosions. The ATLAS light curve shows evidence of a short-lived, pulse-like early excess, similar to that detected in another luminous type Ia supernova (SN 2020hvf). The light curve is broad and the earl… ▽ More We present observations and analysis of the host-less and luminous type Ia supernova 2022ilv, illustrating it is part of the 2003fg-like family, often referred to as super-Chandrasekhar (Ia-SC) explosions. The ATLAS light curve shows evidence of a short-lived, pulse-like early excess, similar to that detected in another luminous type Ia supernova (SN 2020hvf). The light curve is broad and the early spectra are remarkably similar to SN 2009dc. Adopting a redshift of $z=0.026 \pm 0.005$ for SN 2022ilv based on spectral matching, our model light curve requires a large $^{56}$Ni mass in the range $0.7-1.5$ M$_{\odot}$, and a large ejecta mass in the range $1.6-2.3$ M$_{\odot}$. The early excess can be explained by fast-moving SN ejecta interacting with a thin, dense shell of circumstellar material close to the progenitor ($\sim 10^{13}$ cm), a few hours after the explosion. This may be realised in a double-degenerate scenario, wherein a white dwarf merger is preceded by ejection of a small amount ($\sim 10^{-3}-10^{-2}$ M$_{\odot}$) of hydrogen and helium-poor tidally stripped material. A deep pre-explosion Pan-STARRS1 stack indicates no host galaxy to a limiting magnitude of $r \sim 24.5$. This implies a surprisingly faint limit for any host of $M_r \gtrsim -11$, providing further evidence that these types of explosion occur predominantly in low-metallicity environments. △ Less

Submitted 22 January, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: Accepted to ApJL after minor revision

StaNdaRT: A repository of standardized test models and outputs for supernova radiative transfer

Authors: Stéphane Blondin, Sergei Blinnikov, Fionntan P. Callan, Christine E. Collins, Luc Dessart, Wesley Even, Andreas Flörs, Andrew G. Fullard, D. John Hillier, Anders Jerkstrand, Daniel Kasen, Boaz Katz, Wolfgang Kerzendorf, Alexandra Kozyreva, Jack O'Brien, Ezequiel A. Pássaro, Nathaniel Roth, Ken J. Shen, Luke Shingles, Stuart A. Sim, Jaladh Singhal, Isaac G. Smith, Elena Sorokina, Victor P. Utrobin, Christian Vogl , et al. (4 additional authors not shown)

Abstract: We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M… ▽ More We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M$_\odot$) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time, are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe in detail the test models, radiative-transfer codes and output formats and provide access to the repository. We present example results of several key diagnostic features. △ Less

Submitted 15 April, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A. 27 pages, 12 figures (v4: updated to match published version). The ejecta models and output files from the simulations are available at https://github.com/sn-rad-trans/data1

Journal ref: A&A 668, A163 (2022)

3D radiative transfer kilonova modelling for binary neutron star merger simulations

Authors: Christine E. Collins, Andreas Bauswein, Stuart A. Sim, Vimal Vijayan, Gabriel Martínez-Pinedo, Oliver Just, Luke J. Shingles, Markus Kromer

Abstract: The detection of GW170817 and the accompanying electromagnetic counterpart, AT2017gfo, have provided an important set of observational constraints for theoretical models of neutron star mergers, nucleosynthesis, and radiative transfer for kilonovae. We apply the 3D Monte Carlo radiative transfer code ARTIS to produce synthetic light curves of the dynamical ejecta from a neutron star merger, which… ▽ More The detection of GW170817 and the accompanying electromagnetic counterpart, AT2017gfo, have provided an important set of observational constraints for theoretical models of neutron star mergers, nucleosynthesis, and radiative transfer for kilonovae. We apply the 3D Monte Carlo radiative transfer code ARTIS to produce synthetic light curves of the dynamical ejecta from a neutron star merger, which has been modelled with 3D smooth-particle hydrodynamics (SPH) and included neutrino interactions. Nucleosynthesis calculations provide the energy released from radioactive decays of r-process nuclei, and radiation transport is performed using grey opacities given as functions of the electron fraction. We present line-of-sight dependent bolometric light curves, and find the emission along polar lines of sight to be up to a factor of ~2 brighter than along equatorial lines of sight. Instead of a distinct emission peak, our bolometric light curve exhibits a monotonic decline, characterised by a shoulder at the time when the bulk ejecta becomes optically thin. We show approximate band light curves based on radiation temperatures and compare these to the observations of AT2017gfo. We find that the rapidly declining temperatures lead to a blue to red colour evolution similar to that shown by AT2017gfo. We also investigate the impact of an additional, spherically symmetric secular ejecta component, and we find that the early light curve remains nearly unaffected, while after about 1 day the emission is strongly enhanced and dominated by the secular ejecta, leading to the shift of the shoulder from 1-2 to 6-10 days. △ Less

Submitted 23 February, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: Accepted by MNRAS

Modelling the ionisation state of Type Ia supernovae in the nebular-phase

Authors: Luke J. Shingles, Andreas Flörs, Stuart A. Sim, Christine E. Collins, Friedrich K. Roepke, Ivo R. Seitenzahl, Ken J. Shen

Abstract: The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a me… ▽ More The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a method to determine the collisional excitation conditions from [Fe II] line ratios independently from the ionisation state and find that it cannot be applied to highly-ionised models due to the influence of recombination cascades on Fe$^{+}$ level populations. When the ionisation state is artificially lowered, the line ratios (and excitation conditions) are too similar to distinguish between explosion scenarios. We investigate changes to the treatment of non-thermal energy deposition as a way to reconcile over-ionised theoretical models with observations and find that a simple work function approximation provides closer agreement with the data for sub-Mch models than a detailed Spencer-Fano treatment with widely-used cross section data. To quantify the magnitude of additional heating processes that would be required to sufficiently reduce ionisation from fast leptons, we artificially boost the rate of energy loss to free electrons. We find that the equivalent of as much as an eight times increase to the plasma loss rate would be needed to reconcile the sub-Mch model with observed spectra. Future studies could distinguish between reductions in the non-thermal ionisation rates and increased recombination rates, such as by clumping. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 13 pages. Accepted by MNRAS

On the fate of the secondary white dwarf in double-degenerate double-detonation Type Ia supernovae

Authors: R. Pakmor, F. P. Callan, C. E. Collins, S. E. de Mink, A. Holas, W. E. Kerzendorf, M. Kromer, P. G. Neunteufel, John T. O'Brien, F. K. Roepke, A. J. Ruiter, I. R. Seitenzahl, Luke J. Shingles, S. A. Sim, S. Taubenberger

Abstract: The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and exp… ▽ More The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a $1.05\,\mathrm{M_\odot}$ and a $0.7\,\mathrm{M_\odot}$ carbon-oxygen white dwarf with $0.03\,\mathrm{M_\odot}$ helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at $v~>~15\,000$\,km\,s$^{-1}$ are indistinguishable. Light curves and spectra are very similar until $\sim~40\,\mathrm{d}$ after explosion and the ejecta are much more spherical than violent merger models. The inner ejecta differ significantly slowing down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by $\sim20$ per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model. △ Less

Submitted 25 October, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 12 pages, 7 figures, accepted to MNRAS, comments welcome

The influence of line opacity treatment in STELLA on supernova light curves

Authors: Alexandra Kozyreva, Luke Shingles, Alexey Mironov, Petr Baklanov, Sergey Blinnikov

Abstract: We systematically explore the effect of the treatment of line opacity on supernova light curves. We find that it is important to consider line opacity for both scattering and absorption (i.e. thermalisation which mimics the effect of fluorescence.) We explore the impact of degree of thermalisation on three major types of supernovae: Type Ia, Type II-peculiar, and Type II-plateau. For that we use r… ▽ More We systematically explore the effect of the treatment of line opacity on supernova light curves. We find that it is important to consider line opacity for both scattering and absorption (i.e. thermalisation which mimics the effect of fluorescence.) We explore the impact of degree of thermalisation on three major types of supernovae: Type Ia, Type II-peculiar, and Type II-plateau. For that we use radiative transfer code STELLA and analyse broad-band light curves in the context of simulations done with the spectral synthesis code ARTIS and in the context a few examples of observed supernovae of each type. We found that the plausible range for the ratio between absorption and scattering in the radiation hydrodynamics code STELLA is (0.8-1):(0.2-0), i.e. the recommended thermalisation parameter is 0.9. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS

Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae

Authors: L. J. Shingles, S. A. Sim, M. Kromer, K. Maguire, M. Bulla, C. Collins, C. P. Ballance, A. S. Michel, C. A. Ramsbottom, F. K. Roepke, I. R. Seitenzahl, N. B. Tyndall

Abstract: We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-the… ▽ More We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 15 pages, 12 figures. Accepted for publication by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 492, Issue 2, February 2020, Pages 2029-2043

A year-long plateau in the late-time near-infrared light curves of Type Ia supernovae

Authors: Or Graur, Kate Maguire, Russell Ryan, Matt Nicholl, Arturo Avelino, Adam G. Riess, Luke Shingles, Ivo R. Seitenzahl, Robert Fisher

Abstract: The light curves of Type Ia supernovae are routinely used to constrain cosmology models. Driven by radioactive decay of 56Ni, the light curves steadily decline over time, but >150 days past explosion, the near-infrared portion is poorly characterized. We report a year-long plateau in the near-infrared light curve at 150-500 days, followed by a second decline phase accompanied by a possible appeara… ▽ More The light curves of Type Ia supernovae are routinely used to constrain cosmology models. Driven by radioactive decay of 56Ni, the light curves steadily decline over time, but >150 days past explosion, the near-infrared portion is poorly characterized. We report a year-long plateau in the near-infrared light curve at 150-500 days, followed by a second decline phase accompanied by a possible appearance of [Fe I] emission lines. This near-infrared plateau contrasts sharply with Type IIP plateaus and requires a new physical mechanism. We suggest a such as masking of the "near-infrared catastrophe," a predicted yet unobserved sharp light-curve decline, by scattering of ultraviolet photons to longer wavelengths. The transition off the plateau could be due to a change in the dominant ionization state of the supernova ejecta. Our results shed new light on the complex radiative transfer processes that take place in Type Ia supernovae and enhance their use as "standard candles." △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: Authors' version of article published in Nature Astronomy. 12 pages, 6 figures, supplementary data. Published (open) version available here: https://rdcu.be/bTq6b

Using late-time optical and near-infrared spectra to constrain Type Ia supernova explosion properties

Authors: K. Maguire, S. A. Sim, L. Shingles, J. Spyromilio, A. Jerkstrand, M. Sullivan, T. -W. Chen, R. Cartier, G. Dimitriadis, C. Frohmaier, L. Galbany, C. P. Gutiérrez, G. Hosseinzadeh, D. A. Howell, C. Inserra, R. Rudy, J. Sollerman

Abstract: The late-time spectra of Type Ia supernovae (SNe Ia) are powerful probes of the underlying physics of their explosions. We investigate the late-time optical and near-infrared spectra of seven SNe Ia obtained at the VLT with XShooter at $>$200 d after explosion. At these epochs, the inner Fe-rich ejecta can be studied. We use a line-fitting analysis to determine the relative line fluxes, velocity s… ▽ More The late-time spectra of Type Ia supernovae (SNe Ia) are powerful probes of the underlying physics of their explosions. We investigate the late-time optical and near-infrared spectra of seven SNe Ia obtained at the VLT with XShooter at $>$200 d after explosion. At these epochs, the inner Fe-rich ejecta can be studied. We use a line-fitting analysis to determine the relative line fluxes, velocity shifts, and line widths of prominent features contributing to the spectra ([Fe II], [Ni II], and [Co III]). By focussing on [Fe II] and [Ni II] emission lines in the ~7000-7500 Å region of the spectrum, we find that the ratio of stable [Ni II] to mainly radioactively-produced [Fe II] for most SNe Ia in the sample is consistent with Chandrasekhar-mass delayed-detonation explosion models, as well as sub-Chandrasekhar mass explosions that have metallicity values above solar. The mean measured Ni/Fe abundance of our sample is consistent with the solar value. The more highly ionised [Co III] emission lines are found to be more centrally located in the ejecta and have broader lines than the [Fe II] and [Ni II] features. Our analysis also strengthens previous results that SNe Ia with higher Si II velocities at maximum light preferentially display blueshifted [Fe II] 7155 Å lines at late times. Our combined results lead us to speculate that the majority of normal SN Ia explosions produce ejecta distributions that deviate significantly from spherical symmetry. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Comments: 17 pages, 12 figure, accepted for publication in MNRAS

A kilonova as the electromagnetic counterpart to a gravitational-wave source

Authors: S. J. Smartt, T. -W. Chen, A. Jerkstrand, M. Coughlin, E. Kankare, S. A. Sim, M. Fraser, C. Inserra, K. Maguire, K. C. Chambers, M. E. Huber, T. Kruhler, G. Leloudas, M. Magee, L. J. Shingles, K. W. Smith, D. R. Young, J. Tonry, R. Kotak, A. Gal-Yam, J. D. Lyman, D. S. Homan, C. Agliozzo, J. P. Anderson, C. R. Angus C. Ashall , et al. (96 additional authors not shown)

Abstract: Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a r… ▽ More Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements. △ Less

Submitted 17 October, 2017; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: Nature, in press, DOI 10.1038/nature24303. Data files will be made available at http://www.pessto.org

A chemical signature from fast-rotating low-metallicity massive stars: ROA 276 in omega Centauri

Authors: David Yong, John E. Norris, Gary S. Da Costa, Laura M. Stanford, Amanda I. Karakas, Luke J. Shingles, Raphael Hirschi, Marco Pignatari

Abstract: We present a chemical abundance analysis of a metal-poor star, ROA 276, in the stellar system omega Centauri. We confirm that this star has an unusually high [Sr/Ba] abundance ratio. Additionally, ROA 276 exhibits remarkably high abundance ratios, [X/Fe], for all elements from Cu to Mo along with normal abundance ratios for the elements from Ba to Pb. The chemical abundance pattern of ROA 276, rel… ▽ More We present a chemical abundance analysis of a metal-poor star, ROA 276, in the stellar system omega Centauri. We confirm that this star has an unusually high [Sr/Ba] abundance ratio. Additionally, ROA 276 exhibits remarkably high abundance ratios, [X/Fe], for all elements from Cu to Mo along with normal abundance ratios for the elements from Ba to Pb. The chemical abundance pattern of ROA 276, relative to a primordial omega Cen star ROA 46, is best fit by a fast-rotating low-metallicity massive stellar model of 20 Msun, [Fe/H] = -1.8, and an initial rotation 0.4 times the critical value; no other nucleosynthetic source can match the neutron-capture element distribution. ROA 276 arguably offers the most definitive proof to date that fast-rotating massive stars contributed to the production of heavy elements in the early Universe. △ Less

Submitted 22 February, 2017; originally announced February 2017.

Comments: ApJ in press

Evolution and nucleosynthesis of helium-rich asymptotic giant branch models

Authors: Luke J. Shingles, Carolyn L. Doherty, Amanda I. Karakas, Richard J. Stancliffe, John C. Lattanzio, Maria Lugaro

Abstract: There is now strong evidence that some stars have been born with He mass fractions as high as $Y \approx 0.40$ (e.g., in $ω$ Centauri). However, the advanced evolution, chemical yields, and final fates of He-rich stars are largely unexplored. We investigate the consequences of He-enhancement on the evolution and nucleosynthesis of intermediate-mass asymptotic giant branch (AGB) models of 3, 4, 5,… ▽ More There is now strong evidence that some stars have been born with He mass fractions as high as $Y \approx 0.40$ (e.g., in $ω$ Centauri). However, the advanced evolution, chemical yields, and final fates of He-rich stars are largely unexplored. We investigate the consequences of He-enhancement on the evolution and nucleosynthesis of intermediate-mass asymptotic giant branch (AGB) models of 3, 4, 5, and 6 M$_\odot$ with a metallicity of $Z = 0.0006$ ([Fe/H] $\approx -1.4$). We compare models with He-enhanced compositions ($Y=0.30, 0.35, 0.40$) to those with primordial He ($Y=0.24$). We find that the minimum initial mass for C burning and super-AGB stars with CO(Ne) or ONe cores decreases from above our highest mass of 6 M$_\odot$ to $\sim$ 4-5 M$_\odot$ with $Y=0.40$. We also model the production of trans-Fe elements via the slow neutron-capture process (s-process). He-enhancement substantially reduces the third dredge-up efficiency and the stellar yields of s-process elements (e.g., 90% less Ba for 6 M$_\odot$, $Y=0.40$). An exception occurs for 3 M$_\odot$, where the near-doubling in the number of thermal pulses with $Y=0.40$ leads to $\sim$ 50% higher yields of Ba-peak elements and Pb if the $^{13}$C neutron source is included. However, the thinner intershell and increased temperatures at the base of the convective envelope with $Y=0.40$ probably inhibit the $^{13}$C neutron source at this mass. Future chemical evolution models with our yields might explain the evolution of s-process elements among He-rich stars in $ω$ Centauri. △ Less

Submitted 2 July, 2015; originally announced July 2015.

Comments: 21 pages, 16 figures, accepted for publication by MNRAS. Stellar yields included as online data tables

Journal ref: MNRAS, 2015, 452 (3): 2804-2821

Iron and s-elements abundance variations in NGC5286: comparison with anomalous globular clusters and Milky Way satellites

Authors: A. F. Marino, A. P. Milone, A. I. Karakas, L. Casagrande, D. Yong, L. Shingles, G. Da Costa, J. E. Norris, P. B. Stetson, K. Lind, M. Asplund, R. Collet, H. Jerjen, L. Sbordone, A. Aparicio, S. Cassisi

Abstract: We present a high resolution spectroscopic analysis of 62 red giants in the Milky Way globular cluster NGC5286. We have determined abundances of representative light proton-capture, alpha, Fe-peak and neutron-capture element groups, and combined them with photometry of multiple sequences observed along the colour-magnitude diagram. Our principal results are: (i) a broad, bimodal distribution in s-… ▽ More We present a high resolution spectroscopic analysis of 62 red giants in the Milky Way globular cluster NGC5286. We have determined abundances of representative light proton-capture, alpha, Fe-peak and neutron-capture element groups, and combined them with photometry of multiple sequences observed along the colour-magnitude diagram. Our principal results are: (i) a broad, bimodal distribution in s-process element abundance ratios, with two main groups, the s-poor and s-rich groups; (ii) substantial star-to-star Fe variations, with the s-rich stars having higher Fe, e.g. <[Fe/H]>_s-rich - <[Fe/H]>_s-poor ~ 0.2~dex; and (iii) the presence of O-Na-Al (anti-)correlations in both stellar groups. We have defined a new photometric index, c_{BVI}=(B-V)-(V-I), to maximise the separation in the colour-magnitude diagram between the two stellar groups with different Fe and s-element content, and this index is not significantly affected by variations in light elements (such as the O-Na anticorrelation). The variations in the overall metallicity present in NGC5286 add this object to the class of "anomalous" GCs. Furthermore, the chemical abundance pattern of NGC5286 resembles that observed in some of the anomalous GCs, e.g. M22, NGC1851, M2, and the more extreme Omega Centauri, that also show internal variations in s-elements, and in light elements within stars with different Fe and s-elements content. In view of the common variations in s-elements, we propose the term s-Fe-anomalous GCs to describe this sub-class of objects. The similarities in chemical abundance ratios between these objects strongly suggest similar formation and evolution histories, possibly associated with an origin in tidally disrupted dwarf satellites. △ Less

Submitted 26 February, 2015; originally announced February 2015.

Comments: 28 pages, 21 figures, accepted for publication in MNRAS

The s-Process Enrichment of the Globular Clusters M4 and M22

Authors: Luke J. Shingles, Amanda I. Karakas, Raphael Hirschi, Cherie K. Fishlock, David Yong, Gary S. Da Costa, Anna F. Marino

Abstract: We investigate the enrichment in elements produced by the slow neutron-capture process ($s$-process) in the globular clusters M4 (NGC 6121) and M22 (NGC 6656). Stars in M4 have homogeneous abundances of Fe and neutron-capture elements, but the entire cluster is enhanced in $s$-process elements (Sr, Y, Ba, Pb) relative to other clusters with a similar metallicity. In M22, two stellar groups exhibit… ▽ More We investigate the enrichment in elements produced by the slow neutron-capture process ($s$-process) in the globular clusters M4 (NGC 6121) and M22 (NGC 6656). Stars in M4 have homogeneous abundances of Fe and neutron-capture elements, but the entire cluster is enhanced in $s$-process elements (Sr, Y, Ba, Pb) relative to other clusters with a similar metallicity. In M22, two stellar groups exhibit different abundances of Fe and $s$-process elements. By subtracting the mean abundances of $s$-poor from $s$-rich stars, we derive $s$-process residuals or empirical $s$-process distributions for M4 and M22. We find that the $s$-process distribution in M22 is more weighted toward the heavy $s$-peak (Ba, La, Ce) and Pb than M4, which has been enriched mostly with light $s$-peak elements (Sr, Y, Zr). We construct simple chemical evolution models using yields from massive star models that include rotation, which dramatically increases $s$-process production at low metallicity. We show that our massive star models with rotation rates of up to 50\% of the critical (break-up) velocity and changes to the preferred $^{17}$O($α$,$γ$)$^{21}$Ne rate produce insufficient heavy $s$-elements and Pb to match the empirical distributions. For models that incorporate AGB yields, we find that intermediate-mass yields (with a $^{22}$Ne neutron source) alone do not reproduce the light-to-heavy $s$-element ratios for M4 and M22, and that a small contribution from models with a $^{13}$C pocket is required. With our assumption that $^{13}$C pockets form for initial masses below a transition range between 3.0 and 3.5 M$_\odot$, we match the light-to-heavy s-element ratio in the s-process residual of M22 and predict a minimum enrichment timescale of between 240 and 360 Myr. Our predicted value is consistent with the 300 Myr upper limit age difference between the two groups derived from isochrone fitting. △ Less

Submitted 3 September, 2014; originally announced September 2014.

Comments: 13 pages. Accepted for publication in ApJ

Journal ref: ApJ, 795, 34 (2014)

Iron and neutron-capture element abundance variations in the globular cluster M2 (NGC 7089)

Authors: David Yong, Ian U. Roederer, Frank Grundahl, Gary S. Da Costa, Amanda I. Karakas, John E. Norris, Wako Aoki, Cherie K. Fishlock, A. F. Marino, A. P. Milone, Luke J. Shingles

Abstract: We present CN and CH indices and CaII triplet metallicities for 34 giant stars and chemical abundances for 33 elements in 14 giants in the globular cluster M2. Assuming the program stars are cluster members, our analysis reveals (i) an extreme variation in CN and CH line strengths, (ii) a metallicity dispersion with a dominant peak at [Fe/H] = -1.7 and smaller peaks at -1.5 and -1.0, (iii) star-to… ▽ More We present CN and CH indices and CaII triplet metallicities for 34 giant stars and chemical abundances for 33 elements in 14 giants in the globular cluster M2. Assuming the program stars are cluster members, our analysis reveals (i) an extreme variation in CN and CH line strengths, (ii) a metallicity dispersion with a dominant peak at [Fe/H] = -1.7 and smaller peaks at -1.5 and -1.0, (iii) star-to-star abundance variations and correlations for the light elements O, Na, Al and Si and (iv) a large (and possibly bimodal) distribution in the abundances of all elements produced mainly via the s-process in solar system material. Following Roederer et al. (2011), we define two groups of stars, "r+s" and "r-only", and subtract the average abundances of the latter from the former group to obtain a "s-process residual". This s-process residual is remarkably similar to that found in M22 and in M4 despite the range in metallicity covered by these three systems. With recent studies identifying a double subgiant branch in M2 and a dispersion in Sr and Ba abundances, our spectroscopic analysis confirms that this globular cluster has experienced a complex formation history with similarities to M22, NGC 1851 and omega Centauri. △ Less

Submitted 28 April, 2014; originally announced April 2014.

Comments: Accepted for publication in MNRAS

Augmented Reality in Astrophysics

Authors: Frédéric P. A. Vogt, Luke J. Shingles

Abstract: Augmented Reality consists of merging live images with virtual layers of information. The rapid growth in the popularity of smartphones and tablets over recent years has provided a large base of potential users of Augmented Reality technology, and virtual layers of information can now be attached to a wide variety of physical objects. In this article, we explore the potential of Augmented Reality… ▽ More Augmented Reality consists of merging live images with virtual layers of information. The rapid growth in the popularity of smartphones and tablets over recent years has provided a large base of potential users of Augmented Reality technology, and virtual layers of information can now be attached to a wide variety of physical objects. In this article, we explore the potential of Augmented Reality for astrophysical research with two distinct experiments: (1) Augmented Posters and (2) Augmented Articles. We demonstrate that the emerging technology of Augmented Reality can already be used and implemented without expert knowledge using currently available apps. Our experiments highlight the potential of Augmented Reality to improve the communication of scientific results in the field of astrophysics. We also present feedback gathered from the Australian astrophysics community that reveals evidence of some interest in this technology by astronomers who experimented with Augmented Posters. In addition, we discuss possible future trends for Augmented Reality applications in astrophysics, and explore the current limitations associated with the technology. This Augmented Article, the first of its kind, is designed to allow the reader to directly experiment with this technology. △ Less

Submitted 23 May, 2013; originally announced May 2013.

Comments: 15 pages, 11 figures. Accepted for publication in Ap&SS. The final publication will be available at link.springer.com

Is the Sulphur Anomaly in Planetary Nebulae Caused by the s-Process?

Authors: Luke J. Shingles, Amanda I. Karakas

Abstract: Motivated by unexplained observations of low sulphur abundances in planetary nebulae (PNe) and the PG1159 class of post asymptotic giant branch (AGB) stars, we investigate the possibility that sulphur may be destroyed by nucleosynthetic processes in low-to-intermediate mass stars during stellar evolution. We use a 3 Msun, Z=0.01 evolutionary sequence to examine the consequences of high and low rea… ▽ More Motivated by unexplained observations of low sulphur abundances in planetary nebulae (PNe) and the PG1159 class of post asymptotic giant branch (AGB) stars, we investigate the possibility that sulphur may be destroyed by nucleosynthetic processes in low-to-intermediate mass stars during stellar evolution. We use a 3 Msun, Z=0.01 evolutionary sequence to examine the consequences of high and low reaction rate estimates of neutron captures onto sulphur and neighbouring elements. In addition, we have tested high and low rates for the neutron producing reactions C13(alpha,n)O16 and Ne22(alpha,n)Mg25. We vary the mass width of a partially mixed zone (PMZ), which is responsible for the formation of a C13 pocket and is the site of the C13(alpha,n)O16 neutron source. We test PMZ masses from zero up to an extreme upper limit of the entire He-intershell mass at 10^-2 Msun. We find that the alternative reaction rates and variations to the partially mixed zone have almost no effect on surface sulphur abundances and do not reproduce the anomaly. To understand the effect of initial mass on our conclusions, 1.8 Msun and 6 Msun evolutionary sequences are also tested with similar results for sulphur abundances. We are able to set a constraint on the size of the PMZ, as PMZ sizes that are greater than half of the He-intershell mass (in the 3 Msun model) are excluded by comparison with neon abundances in planetary nebulae. We compare the 1.8 Msun model's intershell abundances with observations of PG1159-035, whose surface abundances are thought to reflect the intershell composition of a progenitor AGB star. We find general agreement between the patterns of F, Ne, Si, P, and Fe abundances and a very large discrepancy for sulphur where our model predicts abundances that are 30-40 times higher than is observed in the star. △ Less

Submitted 27 February, 2013; v1 submitted 27 February, 2013; originally announced February 2013.

Comments: 11 pages. Accepted for publication in MNRAS

Journal ref: MNRAS 431, 2861-2871 (2013)