arXiv.org > astro-ph: sn

Electron Energization in Quasi-Parallel Shocks: Test-Particle-Electrons in a Proton Driven Turbulence

Hanusch, Adrian; Liseykina, Tatyana; Malkov, Mikhail

In situ observations of energetic particles at the Earth's bow-shock attainable by the satellite missions have long created an opinion that electrons are most efficiently accelerated in a quasi-perpendicular shock geometry. However, shocks that deemed responsible for the production of cosmic ray electrons and their radiation from sources such as supernova remnants are much more powerful and larger than the Earth's bow-shock. Their remote observations suggest that electrons are accelerated very efficiently in the quasi-parallel shocks as well. In this paper we investigate the possibility that protons accelerated to high energies create sufficient wave turbulence necessary for the electron preheating and subsequent injection into the diffusive shock acceleration in a quasi-parallel shock geometry. An additional test-particle-electron population, meant as a low-density addition to the electron core-distribution the hybrid simulation operates on, is introduced. Its purpose is to investigate how these electrons are energized by the "hybrid" electromagnetic field. The reduced spatial dimensionality allowed us to dramatically increase the number of macro-ions per numerical cell and achieve the converged results for the velocity distributions of test-electrons. We discuss the electron preheating mechanisms which can make a significant part of thermal electrons accessible to the ion-driven waves observed in hybrid simulations. We find that the precursor wave field supplied by ions has considerable potential to preheat the electrons before they are shocked at the subshock. Our results indicate that a downstream thermal equilibration of the hot test-electrons and protons does not occur. Instead, the resulting electron-to-proton temperature ratio is a decreasing function of the shock Mach number, $M_\mathrm{A},$ tending to a saturation at high $M_\mathrm{A}.$

Convection and Spin-Up During Common Envelope Evolution: The Formation of Short-Period Double White Dwarfs

Wilson, E. C.; Nordhaus, J.

The formation channels and predicted populations of double-white dwarfs (DWDs) are important because a subset will evolve to be gravitational-wave sources and/or progenitors of Type Ia supernovae. Given the observed population of short-period DWDs, we calculate the outcomes of common envelope evolution when convective effects are included. For each observed white dwarf in a DWD system, we identify all progenitor stars with an equivalent proto-WD core mass from a comprehensive suite of stellar evolution models. With the second observed white dwarf as the companion, we calculate the conditions under which convection can accommodate the energy released as the orbit decays, including (if necessary), how much the envelope must spin-up during the common envelope phase. The predicted post-CE final separations closely track the observed DWD orbital parameter space, further strengthening the view that convection is a key ingredient in common envelope evolution.

Combined analysis of Planck and SPTPol data favors the early dark energy models

Chudaykin, Anton; Gorbunov, Dmitry; Nedelko, Nikita

We study the implications of the Planck temperature power spectrum at low multipoles, $\ell<1000$, and SPTPol data. We show that this combination predicts consistent lensing-induced smoothing of acoustic peaks within $\Lambda$CDM cosmology and yields the robust predictions of the cosmological parameters. Combining only the Planck large-scale temperature data and the SPTPol polarization and lensing measurements within $\Lambda$CDM model we found substantially lower values of linear matter density perturbation $\sigma_8$ which bring the late-time parameter $S_8=\sigma_8\sqrt{\Omega_m/0.3}=0.763\pm0.022$ into accordance with galaxy clustering and weak lensing measurements. It also raises up the Hubble constant $H_0=69.68\pm1.00{\rm \,\,km\,s^{-1}Mpc^{-1}}$ that reduces the Hubble tension to the $2.5\sigma$ level. We examine the residual tension in the Early Dark Energy (EDE) model which produces the brief energy injection prior to recombination. We implement both the background and perturbation evolutions of the scalar field which potential scales as $V(\phi)\propto \phi^{2n}$. Including cosmic shear measurements (KiDS, VIKING-450, DES) and local distance-ladder data (SH0ES) to the combined fit we found that EDE completely alleviates the Hubble tension while not degradating the fit to large-scale structure data. The EDE scenario significantly improves the goodness-of-fit by $2.9\sigma$ in comparison with the concordance $\Lambda$CDM model. The account for the intermediate-redshift data (the supernova dataset and baryon acoustic oscillation data) fits perfectly to our parameter predictions and indicates the preference of EDE over $\Lambda$CDM at $3\sigma$.

Discovery of a hot ultramassive rapidly rotating DBA White Dwarf

We report the discovery of a nearby massive white dwarf with He-H atmosphere. The white dwarf is located at a distance of $74.5\pm0.9$ pc. Its radius, mass, effective temperature, H/He ratio and age are $R=2500\pm100$ km, $M=1.33\pm0.01$ $\rm M_{\odot}$, $T_{\mathrm{eff}}=31200\pm 1200$ K, ${\mathrm{H/He}}\sim0.1$ and $330\pm40$ Myr, respectively. The observed spectrum is redshifted by $V_{\mathrm{r}}=+240\pm15$ km s$^{-1},$ which is mostly attributed to the gravitational redshift. The white dwarf shows a regular stable photometric variability with amplitude $\Delta g\approx0.06^{\mathrm{m}}$ and period $P=353.456$ s suggesting rapid rotation. This massive, hot and rapidly rotating white dwarf is likely to originate from the merging of close binary white dwarf system that avoided explosion in a thermonuclear type Ia supernova at the Carboniferous Period of the Earth history.

Consistency of cubic Galileon cosmology: Model-independent bounds from Pantheon data

Brahma, Suddhasattwa; Hossain, Md. Wali

We revisit the cosmological dynamics of the cubic Galileon model in light of the recently proposed model-independent analyses of the Pantheon supernova data. At the background level, it is shown to be compatible with data and preferred over standard quintessence models. Furthermore, the model is shown to be consistent with the trans-Planckian censorship conjecture (as well as other Swampland conjectures). At the level of linear perturbations, falsifiability of the model is established from the behaviour of the growth index at large redshift.

A cosmic Zevatron based on cyclotron auto-resonance

Salamin, Yousef I.; Wen, Meng; Keitel, Christoph H.

A Zevatron is an accelerator scheme envisaged to accelerate particles to ZeV energies (1 ZeV = $10^{21}$ eV). Schemes, most notably the internal shock model, have been proposed to explain the acceleration of ultra-high-energy cosmic-ray (UHECR) particles that have been sporadically detected reaching Earth since 1962. Here, the cyclotron auto-resonance acceleration (CARA) mechanism is tailored and used to demonstrate possible acceleration of particles ejected as a result of violent astrophysical processes such as the merger of a binary system or a supernova explosion. Such events result in emission of highly-energetic particles and ultra-intense beamed radiation. In the simultaneous presence of a super-strong magnetic field, the condition for cyclotron auto-resonance may be met. Thus CARA can act like a {\it booster} for particles pre-accelerated inside their progenitor by shock waves, possibly among other means. As examples, it is shown that nuclei of hydrogen, helium, and iron-56, may reach ZeV energies by CARA, under which conditions the particles, while gyrating around the lines of an ultra-strong magnetic field, also surf on the waves of a super-intense radiation field. When radiation-reaction is taken into account, it is shown that the ZeV energy gained by a particle can fall off by less than an order-of-magnitude if the resonance condition is missed by roughly less than 20\%.

Luminous supernovae associated with ultra-long gamma-ray bursts from hydrogen-free progenitors extended by pulsational pair-instability

Moriya, Takashi J.; Marchant, Pablo; Blinnikov, Sergei I.

We show that the luminous supernovae associated with ultra-long gamma-ray bursts can be caused by the slow adiabatic cooling of hydrogen-free progenitors extended by pulsational pair-instability. In the accompanying paper (Marchant & Moriya 2020), we have shown that some rapidly-rotating hydrogen-free gamma-ray burst progenitors that experience pulsational pair-instability can keep an extended structure caused by pulsational pair-instability until the core collapse. Such progenitors have large radii exceeding 10 Rsun and they sometimes reach beyond 1000 Rsun at the time of the core collapse. They are, therefore, promising progenitors of ultra-long gamma-ray bursts. We here perform the light-curve modeling of the explosions of one extended hydrogen-free progenitor with a radius of 1962 Rsun. Thanks to the large progenitor radius, the ejecta experience slow adiabatic cooling and they become rapidly evolving (<~ 10 days) luminous (>~ 1e43 erg/s) supernovae even without the energy input from the 56Ni nuclear decay. They also have fast photospheric velocity above 10,000 km/s and hot photospheric temperature above 10,000 K at around the peak luminosity. We find that SN 2011kl, which is associated with the ultra-long gamma-ray burst GRB 111209A, can be explained by the adiabatic cooling phase of the extended progenitor. The ultra-long gamma-ray burst progenitors proposed in Marchant & Moriya 2020 can explain both the ultra-long gamma-ray burst duration and the accompanying supernova properties. When the gamma-ray burst jet is off-axis or choked, the luminous supernovae could be observed as fast blue optical transients without accompanying gamma-ray bursts.

The impact of stellar rotation on the black hole mass-gap from pair-instability supernovae

Marchant, Pablo; Moriya, Takashi

Models of pair-instability supernovae (PISNe) predict a gap in black hole (BH) masses between $\sim 45M_\odot-120M_\odot$, which is referred to as the upper BH mass-gap. With the advent of gravitational-wave astrophysics it has become possible to test this prediction, and there is an important associated effort to understand what theoretical uncertainties modify the boundaries of this gap. In this work we study the impact of rotation on the hydrodynamics of PISNe, which leave no compact remnant, as well as the evolution of pulsational-PISNe (PPISNe), which undergo thermonuclear eruptions before forming a compact object. We perform simulations of non-rotating and rapidly-rotating stripped helium stars in a metal poor environment $(Z_\odot/50)$ in order to resolve the lower edge of the upper mass-gap. We find that the outcome of our simulations is dependent on the efficiency of angular momentum transport, with models that include efficient coupling through the Spruit-Tayler dynamo shifting the lower edge of the mass-gap upwards by $\sim 4\%$, while simulations that do not include this effect shift it upwards by $\sim 15\%$. From this, we expect the lower edge of the upper mass-gap to be dependent on BH spin, which can be tested as the number of observed BH mergers increases. Moreover, we show that stars undergoing PPISNe have extended envelopes ($R\sim 10-1000~R_\odot$) at iron-core collapse, making them promising progenitors for ultra-long gamma-ray bursts.

Plasma Diagnostics of the Supernova Remnant N132D Using Deep XMM-Newton Observations with the Reflection Grating Spectrometer

We present XMM-Newton observations of N132D, the X-ray brightest supernova remnant (SNR) in the Large Magellanic Cloud (LMC), using the Reflection Grating Spectrometer (RGS) that enables high-resolution spectroscopy in the soft X-ray band. A dozen emission lines from L-shell transitions of various elements at intermediate charge states are newly detected in the RGS data integrating the ~200-ks on-axis observations. The 0.3-2.0-keV spectra require at least three components of thermal plasmas with different electron temperatures and indicate clear evidence of non-equilibrium ionization (NEI). Our detailed spectral diagnostics further reveal that the forbidden-to-resonance line ratios of O VII and Ne IX are both higher than expected for typical NEI plasmas. This enhancement could be attributed to either resonance scattering or emission induced by charge exchange in addition to a possible contribution from the superposition of multiple temperature components, although the lack of spatial information prevents us from concluding which is most likely.

Does GW190425 require an alternative formation pathway than a fast-merging channel?

Safarzadeh, Mohammadtaher; Ramirez-Ruiz, Enrico; Berger, Edo

The LIGO/Virgo Scientific Collaboration (LVC) recently announced the detection of a compact object binary merger, GW190425, with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, and individual component masses in the range of about 1.1 to 2.5 $M_{\odot}$. If the constituent compact objects are neutron stars, then the total mass is five standard deviations higher than the mean of $2.66\pm 0.12$ M$_{\odot}$ for Galactic binary neutron stars (BNSs). LVC suggests such massive BNS systems are born from a fast-merging channel, and therefore, their non-detection in the Galaxy to be due to a selection effect. However, we are unable to reconcile the inferred formation efficiency from the reported merger rate, $\mathcal{R}_{\rm GW190425}=460^{+1050}_{-390}$ yr$^{-1}$ Gpc$^{-3}$, with predictions from our own study for fast-merging BNS systems. Moreover, the comparable merger rates of GW190425, and GW170817 are possibly in tension with our results for two reasons: (i) more massive systems are expected to have a lower formation rate, and (ii) fast merging channels should constitute $\lesssim 10\%$ of the total BNS systems if case BB unstable mass transfer is permitted to take place as a formation pathway. We argue that to account for the high merger rate of GW190425 as a BNS system requires: (i) revisiting our understanding of NS formation in supernova explosions, or (ii): that more massive NSs need to be preferentially born with either very weak or very high magnetic fields so that they would be undetectable in the radio surveys. Perhaps massive NSs detected in NS-white dwarf binaries are our clues to the formation path of GW190425 systems.

A Search for the modern counterparts of the Far Eastern guest stars 369 CE, 386 CE, and 393 CE

Hoffmann, Susanne M; Vogt, Nikolaus

In this study, we apply our previously developed method to investigate ancient transient sightings in order to derive consequences for modern astrophysical problems.We present case studies of three observations of so called 'guest stars' in the 4th century CE which lasted several months each. These three observations had been suggested and discussed as possible supernovae but slow novae are also viable alternatives. Our careful re-interpretation of the historical texts and the currently known objects in the given fields shed new light on this topic. In particular, for the two events in 386 and 393 CE we suggest possible supernova identifications, while in all three cases there are interesting candidates for past classical or recurrent nova eruptions among known cataclysmic variables (CVs) and/or symbiotic stars. The transient of 369, we suggest to explain as classical and possibly recurrent nova instead of a supernova. The most plausible candidates are BZ Cam, a CV with a possible nova shell, or CQ Dra, a naked eye multiple system perhaps able to permit an overwhelmingly bright nova with day-time visibility.

Explosion Energies for Core-collapse Supernovae I: Analytic, Spherically Symmetric Solutions

Gogilashvili, Mariam; Murphy, Jeremiah W.; Mabanta, Quintin

Recent multi-dimensional simulations of core-collapse supernovae are producing successful explosions and explosion-energy predictions. In general, the explosion-energy evolution is monotonic and relatively smooth, suggesting a possible analytic solution. We derive analytic solutions for the expansion of the gain region under the following assumptions: spherical symmetry, one-zone shell, and powered by neutrinos and $\alpha$ particle recombination. We consider two hypotheses: I) explosion energy is powered by neutrinos and $\alpha$ recombination, II) explosion energy is powered by neutrinos alone. Under these assumptions, we derive the fundamental dimensionless parameters and analytic scalings. For the neutrino-only hypothesis (II), the asymptotic explosion energy scales as $E_{\infty} \approx 1.5 M_g v_0^2 \eta^{2/3}$, where $M_g$ is the gain mass, $v_0$ is the free-fall velocity at the shock, and $\eta$ is a ratio of the heating and dynamical time scales. Including both neutrinos and recombination (hypothesis I), the asymptotic explosion energy is $E_{\infty} \approx M_g v_0^2 (1.5\eta^{2/3} + \beta f(\rho_0))$, where $\beta$ is the dimensionless recombination parameter. We use Bayesian inference to fit these analytic models to simulations. Both hypotheses fit the simulations of the lowest progenitor masses that tend to explode spherically. The fits do not prefer hypothesis I or II; however, prior investigations suggest that $\alpha$ recombination is important. As expected, neither hypothesis fits the higher-mass simulations that exhibit aspherical explosions. In summary, this explosion-energy theory is consistent with the spherical explosions of low progenitor masses; the inconsistency with higher progenitor-mass simulations suggests that a theory for them must include aspherical dynamics.

Pulsational pair-instability supernovae in very close binaries

Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these events have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor ($Z_\odot/10$) single helium stars using the \texttt{MESA} software instrument. For PPISN progenitors with pre-pulse masses $>50M_\odot$ we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than $100R_\odot$ for periods of $10^2-10^4\;$~yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by $\sim 30\%$, induces eccentricities below the threshold of detectability of the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.

Model independent measurement on the cosmic curvature

Yang, Yingjie; Gong, Yungui

Inflation predicts that the Universe is spatially flat. The Planck 2018 measurements of the cosmic microwave background anisotropy favour a spatially closed universe at more than 2$\sigma$ confidence level. We use model independent methods to study the issue of cosmic curvature. The method reconstructs the Hubble parameter $H(z)$ from cosmic chronometers data with the Gaussian process method. The distance modulus is then calculated with the reconstructed function $H(z)$ and fitted by type Ia supernovae data. Combining the cosmic chronometers and type Ia supernovae data, we obtain $\Omega_{k0}h^2=0.102\pm 0.066$ which shows a spatially open universe at more than 1$\sigma$ confidence level. By adding the redshift space distortions data to the type Ia supernovae data, we obtain $\Omega_{k0}h^2=0.117^{+0.058}_{-0.045}$ and no interesting deviation from $\Lambda$CDM model is found. The almost $2\sigma$ evidence of a spatially open universe is independent of cosmological model and gravitational theory.

Stellar-mass black holes in young massive and open stellar clusters IV: updated stellar-evolutionary and black hole spin models and comparisons with the LIGO-Virgo O1/O2 merger-event data

Banerjee, Sambaran

I present a set of long-term, direct, relativistic many-body computations of model dense stellar clusters with up-to-date stellar-evolutionary, supernova (SN), and remnant natal-kick models, including pair instability and pulsation pair instability supernova (PSN and PPSN), using an updated version of NBODY7 N-body simulation program. The N-body model also includes stellar evolution-based natal spins of BHs and treatments of binary black hole (BBH) mergers based on numerical relativity. These, for the first time in a direct N-body simulation, allow for second-generation BBH mergers. The set of 65 evolutionary models have initial masses $10^4M_\odot-10^5M_\odot$, sizes 1 pc-3 pc, metallicity $0.0001-0.02$, with the massive stars in primordial binaries and they represent young massive clusters (YMC) and moderately massive open clusters (OC). Such models produce dynamically-paired BBH mergers that agree well with the observed masses, mass ratios, effective spin parameters, and final spins of the LVC O1/O2 merger events, provided BHs are born with low or no spin but spin up after undergoing a BBH merger or matter accretion onto it. In particular, the distinctly higher mass, effective spin parameter, and final spin of GW170729 merger event is naturally reproduced, as also the mass asymmetry of the O3 event GW190412. The computed models produce massive, $\sim100M_\odot$ BBH mergers with primary mass within the 'PSN gap' and also yield mergers involving remnants in the 'mass gap'. They also suggest that YMCs and OCs produce persistent, Local-Universe GW sources detectable by LISA. Such clusters are also capable of producing eccentric LIGO-Virgo mergers.

Photometric observations of two type Ic-BL Supernovae: 2016coi and 2018ebt

Tsvetkov, D. Yu.; Pavlyuk, N. N.; Echeistov, V. A.

CCD BVRI photometry is presented for the type Ic-BL Supernovae 2016coi and 2018ebt. The shape of the light curves for both objects is typical for this class of SNe. A change of brightness decline rate at about 100 days after maximum light is observed for SN 2016coi. The light curves of SN 2018ebt are best matched by those of SN Ic-BL 2002ap. For SN 2018ebt, we derived the dates and magnitudes of maximum light, rates of decline at late stages of evolution and presented evidence for extinction in the host galaxy being negligible.

Global fits of axion-like particles to XENON1T and astrophysical data

The excess of electron recoil events seen by the XENON1T experiment has been interpreted as a potential signal of axion-like particles (ALPs), either produced in the Sun, or constituting part of the dark matter halo of the Milky Way. It has also been explained as a consequence of trace amounts of tritium in the experiment. We consider the evidence for the solar and dark-matter ALP hypotheses from the combination of XENON1T data and multiple astrophysical probes, including horizontal branch stars, red giants, and white dwarfs. We briefly address the influence of ALP decays and supernova cooling. While the different datasets are in clear tension for the case of solar ALPs, all measurements can be simultaneously accommodated for the case of a sub-dominant fraction of dark-matter ALPs. Nevertheless, this solution requires the tuning of several a priori unknown parameters, such that for our choices of priors a Bayesian analysis shows no strong preference for the ALP interpretation of the XENON1T excess over the background hypothesis.

Distribution and kinematics of 26Al in the Galactic disc

Fujimoto, Yusuke; Krumholz, Mark R.; Inutsuka, Shu-ichiro

26Al is a short-lived radioactive isotope thought to be injected into the interstellar medium (ISM) by massive stellar winds and supernovae. However, all-sky maps of 26Al emission show a distribution with a much larger scale height and faster rotation speed than either massive stars or the cold ISM. We investigate the origin of this discrepancy using an N-body+hydrodynamics simulation of a Milky-Way-like galaxy, self-consistently including self-gravity, star formation, stellar feedback, and 26Al production. We find no evidence that the Milky Way's spiral structure explains the 26Al anomaly. Stars and the 26Al bubbles they produce form along spiral arms, but, because our simulation produces material arms that arise spontaneously rather than propagating arms forced by an external potential, star formation occurs at arm centres rather than leading edges. As a result, we find a scale height and rotation speed for 26Al similar to that of the cold ISM. However, we also show that a synthetic 26Al emission map produced for a possible Solar position at the edge of a large 26Al bubble recovers many of the major qualitative features of the observed 26Al sky. This suggests that the observed anomalous 26Al distribution is the product of foreground emission from the 26Al produced by a nearby, recent supernova.

An Active-to-Sterile Neutrino Transition Dipole Moment and the XENON1T Excess

Shoemaker, Ian M.; Tsai, Yu-Dai; Wyenberg, Jason

In this short letter, we find that a magnetic transition dipole moment between tau and sterile neutrinos can account for the XENON1T excess events. Unlike the ordinary neutrino dipole moment, the introduction of the new sterile mass scale allows for astrophysical bounds to be suppressed. Interestingly, the best-fit regions that are compatible with the SN1987A imply either boron-8 or CNO neutrinos as the source flux. We find that sterile neutrinos of either $\sim$ 260 keV or in the $\sim$(500 - 800) keV mass range are capable of evading astrophysical constraints while being able to successfully explain the XENON1T event rate. The sterile neutrino in the best fit parameter space may have significant effects on big bang nucleosynthesis (BBN). We show the region in which a low reheating temperature of the Universe may allow the BBN constraints to be alleviated.

Constraining early-time dust formation in core-collapse supernovae

Priestley, F. D.; Bevan, A.; Barlow, M. J.; De Looze, I.

There is currently a severe discrepancy between theoretical models of dust formation in core-collapse supernovae (CCSNe), which predict $\gtrsim 0.01$ M$_\odot$ of ejecta dust forming within $\sim 1000$ days, and observations at these epochs, which infer much lower masses. We demonstrate that, in the optically thin case, these low dust masses are robust despite significant observational and model uncertainties. For a sample of 11 well-observed CCSNe, no plausible model reaches carbon dust masses above $10^{-4}$ M$_\odot$, or silicate masses above $\sim 10^{-3}$ M$_\odot$. Optically thick models can accommodate larger dust masses, but the dust must be clumped and have a low ($<0.1$) covering fraction to avoid conflict with data at optical wavelengths. These values are insufficient to reproduce the observed infrared fluxes, and the required covering fraction varies not only between SNe but between epochs for the same object. The difficulty in reconciling large dust masses with early-time observations of CCSNe, combined with well-established detections of comparably large dust masses in supernova remnants, suggests that a mechanism for late-time dust formation is necessary.

LAMOST Medium-Resolution Spectral Survey of Galactic Nebula (LAMOST MRS-N): An overview of Scientific goals and Survey plan

Since Sep. 2018, LAMOST started the medium-resolution (R $\sim 7500$) spectral survey (MRS). We proposed the spectral survey of Galactic nebulae, including \hii\ regions, HH objects, supernova remnants, planetary nebulae and the special stars with MRS (LAMOST MRS-N). LAMOST MRS-N covers about 1700 square degrees of the northern Galactic plane within $40^\circ < l < 215^\circ$ and $-5^\circ < b < 5^\circ$. In this 5-year survey, we plan to observe about 500 thousand nebulae spectra. According to the commissioning observations, the nebulae spectra can provide precise radial velocity with uncertainty less than 1 km s$^{-1}$. These high-precision spectral data are of great significance to our understanding of star formation and evolution.

Preparing Nuclear Astrophysics for Exascale

Astrophysical explosions such as supernovae are fascinating events that require sophisticated algorithms and substantial computational power to model. Castro and MAESTROeX are nuclear astrophysics codes that simulate thermonuclear fusion in the context of supernovae and X-ray bursts. Examining these nuclear burning processes using high resolution simulations is critical for understanding how these astrophysical explosions occur. In this paper we describe the changes that have been made to these codes to transform them from standard MPI + OpenMP codes targeted at petascale CPU-based systems into a form compatible with the pre-exascale systems now online and the exascale systems coming soon. We then discuss what new science is possible to run on systems such as Summit and Perlmutter that could not have been achieved on the previous generation of supercomputers.

Core-collapse supernova neutrino emission and detection informed by state-of-the-art three-dimensional numerical models

Nagakura, Hiroki; Burrows, Adam; Vartanyan, David; Radice, David

Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star (PNS) convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral Standing Accretion Shock Instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate ($\gtrsim 50 \%$), but that the time-integrated signal has a relatively modest asymmetry ($\lesssim 20 \%$). Both features are associated with the lepton-number emission self-sustained asymmetry (LESA) and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavors of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of $\sim$20\% if the distance to the CCSN is $\lesssim 6$ kpc.

Charged-particle branching ratios above the neutron threshold in $^{19}$F: constraining $^{15}$N production in core-collapse supernovae

Spatially-correlated overabundances of $^{15}$N and $^{18}$O observed in some low-density graphite meteoritic grains have been connected to nucleosynthesis taking place in the helium-burning shell during core-collapse supernovae. Two of the reactions which have been identified as important to the final abundances of $^{15}$N and $^{18}$O are $^{18}$F($n,\alpha$)$^{15}$N and $^{18}$F($n,p$)$^{18}$O. The relative strengths of the $^{18}$F($n,\alpha$)$^{15}$N and $^{18}$F($n,p$)$^{18}$O reactions depend on the relative $\alpha_0$ and $p_0$ decays from states above the neutron threshold in $^{19}$F in addition to other properties. Experimental data on the charged-particle decays from these highly excited states are lacking or inconsistent. Two experiments were performed using proton inelastic scattering from LiF targets and magnetic spectrographs. The first experiment used the high-resolution Q3D spectrograph at Munich to constrain properties of levels in $^{19}$F. A second experiment using the Orsay Split-Pole spectrograph and an array of silicon detectors was performed in order to measure the charged-particle decays of neutron-unbound levels in $^{19}$F. A number of levels in $^{19}$F have been identified along with their corresponding charged-particle decays. The first state above the neutron threshold which has an observed proton-decay branch to the ground state of $^{18}$O lies 68 keV above the neutron threshold while the $\alpha$-particle decays from the neutron-unbound levels are generally observed to be much stronger than the proton decays. Neutron-unbound levels in $^{19}$F are observed to decay predominantly by $\alpha$-particle emission, supporting the role of $^{18}$F($n,\alpha$)$^{15}$N in the production of $^{15}$N in the helium-burning shell of supernovae. Improved resonant-scattering reaction data are required in order to be able to determine the reaction rates accurately.

Exploring the Lower Mass Gap and Unequal Mass Regime in Compact Binary Evolution

Zevin, Michael; Spera, Mario; Berry, Christopher P. L.; Kalogera, Vicky

On August 14, 2019, the LIGO and Virgo detectors observed GW190814, a gravitational-wave signal originating from the merger of a $\simeq 23 M_\odot$ black hole with a $\simeq 2.6 M_\odot$ compact object. GW190814's compact-binary source is atypical both in its highly asymmetric masses and in its lower-mass component lying between the heaviest known neutron star and lightest known black hole in a compact-object binary. If formed through isolated binary evolution, the mass of the secondary is indicative of its mass at birth. We examine the formation of such systems through isolated binary evolution across a suite of assumptions encapsulating many physical uncertainties in massive-star binary evolution. We update how mass loss is implemented for the neutronization process during the collapse of the proto-compact object to eliminate artificial gaps in the mass spectrum at the transition between neutron stars and black holes. We find it challenging for population modeling to match the empirical rate of GW190814-like systems whilst simultaneously being consistent with the rates of other compact binary populations inferred by gravitational-wave observations. Nonetheless, the formation of GW190814-like systems at any measurable rate requires a supernova engine model that acts on longer timescales such that the proto-compact object can undergo substantial accretion immediately prior to explosion, hinting that if GW190814 is the result of massive-star binary evolution, the mass gap between neutron stars and black holes may be narrower or nonexistent.

Eccentric grazing envelope evolution towards type IIb supernova progenitors

Shishkin, Dmitry; Soker, Noam

We simulate the evolution of eccentric binary systems in the frame of the grazing envelope evolution (GEE) channel for the formation of Type IIb supernovae (SNe IIb), and find that extra mass removal by jets increases the parameter space for the formation of SNe IIb in this channel. To explore the role of eccentricity and the extra mass removal by jets we use the stellar evolutionary code MESA binary. The initial primary and secondary masses are M1i=15Mo and M2i=2.5Mo. We examine initial semi-major axes of 600-1000Ro, and eccentricities of e=0-0.9. Both Roche lobe overflow (RLOF) and mass removal by jets, followed by a wind, leave a hydrogen mass in the exploding star of M(H,f)=0.05Mo, compatible with a SN IIb progenitor. When the initial orbit is not circular the final orbit might have a very high eccentricity. In many cases, with and without the extra mass removal by jets, the system can enter a common envelope evolution (CEE) phase, and then gets out from it. For some ranges of eccentricities the jets are more efficient in preventing the CEE. Despite the large uncertainties, extra mass removal by jets substantially increases the likelihood of the system to get out from a CEE. This strengthens earlier conclusions for circular orbits. In some cases RLOF alone, without mass removal by jets, can form SN IIb progenitors. We estimate that the extra mass removal by jets in the GEE channel increases the number of progenitors relative to that by RLOF alone by about a factor of two.

Gravitational-wave Signature of a First-order Quantum Chromodynamics Phase Transition in Core-Collapse Supernovae

Zha, Shuai; O'Connor, Evan P.; Chu, Ming-chung; Lin, Lap-Ming; Couch, Sean M.

A first-order quantum chromodynamics (QCD) phase transition (PT) may take place in the protocompact star (PCS) produced by a core-collapse supernova (CCSN). In this work, we study the consequences of such a PT in a non-rotating CCSN with axisymmetric hydrodynamic simulations. We find that the PT leads to the collapse of the PCS and results in a loud burst of gravitational waves (GWs). The amplitude of this GW burst is $\sim30$ times larger than the post-bounce GW signal normally found for non-rotating CCSN. It shows a broad peak at high frequencies ($\sim2500-4000$ Hz) in the spectrum, has a duration of $\lesssim5 {\rm ms}$, and carries $\sim3$ orders of magnitude more energy than the other episodes. Also, the peak frequency of the PCS oscillation increases dramatically after the PT-induced collapse. In addition to a second neutrino burst, the GW signal, if detected by the ground-based GW detectors, is decisive evidence of the first-order QCD PT inside CCSNe and provides key information about the structure and dynamics of the PCS.

The SiTian project

Liu, JiFeng; Soria, Roberto; Wu, Xue-Feng; Wu, Hong; Shang, Zhaohui

SiTian is an ambitious ground-based all-sky optical monitoring project, developed by the Chinese Academy of Sciences. The concept is an integrated network of dozens of 1-m-class telescopes deployed partly in China and partly at various other sites around the world. The main science goals are the detection, identification and monitoring of optical transients (such as gravitational wave events, fast radio bursts, supernovae) on the largely unknown timescales of less than 1 day; SiTian will also provide a treasure trove of data for studies of AGN, quasars, variable stars, planets, asteroids, and microlensing events. To achieve those goals, SiTian will scan at least 10,000 square deg of sky every 30 min, down to a detection limit of $V \approx 21$ mag. The scans will produce simultaneous light-curves in 3 optical bands. In addition, SiTian will include at least three 4-m telescopes specifically allocated for follow-up spectroscopy of the most interesting targets. We plan to complete the installation of 72 telescopes by 2030 and start full scientific operations in 2032.

The Role of Magnetic Fields in Neutrino-Driven Supernovae

We study the impact of a small-scale dynamo in core-collapse supernovae using a 3D neutrino magnetohydrodynamics simulation of a $15 M_\odot$ progenitor. The weak seed field is amplified exponentially in the gain region once neutrino-driven convection develops, and remains dominated by small-scale structures. About $250\, \mathrm{ms}$ after bounce, the field energy in the gain region reaches $\mathord{\sim} 50\%$ of kinetic equipartition. This supports the development of a neutrino-driven explosion with modest global anisotropy, which does not occur in a corresponding model without magnetic fields. Our results suggest that magnetic fields may play a beneficial subsidiary role in neutrino-driven supernovae even without rapid progenitor rotation. Further investigation into the nature of magnetohydrodynamic turbulence in the supernova core is required.

Probing the Seesaw Mechanism with Cosmological Data

Rodrigues, J. G.; Benetti, Micol; Campista, Marcela; Alcaniz, Jailson

We investigate cosmological consequences of an inflationary model which incorporates a generic seesaw extension (types I and II) of the Standard Model of Particle Physics. A non-minimal coupling between the inflaton field and the Ricci scalar is considered as well as radiative corrections at one loop order. This connection between the inflationary dynamics with neutrino physics results in a predictive model whose observational viability is investigated in light of the current cosmic microwave background data, baryon acoustic oscillation observations and type Ia supernovae measurements. Our results show that the non-minimal coupled seesaw potential provides a good description of the observational data when radiative corrections are positive. Such result favours the type II seesaw mechanism over type I and may be an indication for physics beyond the Standard Model.

On the gamma-ray emission of W44 and its surroundings

Peron, Giada; Aharonian, Felix; Casanova, Sabrina; Zanin, Roberta; Romoli, Carlo

We present the analysis of 9.7 years Fermi-LAT data of the middle-aged supernova remnant W44 and the massive molecular gas complex that surrounds it. We derived a high-quality spectral energy distribution of gamma-radiation of the shell over three decades. The very hard spectrum below 1 GeV supports the earlier claims regarding the hadronic origin of radiation. We also confirm the presence of two extended $\gamma$-ray structures located at two opposite edges of the remnant along its major axis. Based on the high-resolution gas maps, we demonstrate that the gamma-ray structures are caused by the enhanced cosmic-ray density rather than the gradient of the gas distribution. We argue that the revealed cosmic-ray "clouds" suggest an anisotropic character of the escape of high-energy particles from the shell along the magnetic field of the remnant.

SNR G39.2-0.3, an Hadronic Cosmic Rays Accelerator

Recent results obtained with gamma-ray satellites have established supernova remnants as accelerators of GeV hadronic cosmic rays. In such processes, CRs accelerated in SNR shocks interact with particles from gas clouds in their surrounding. In particular, the rich medium in which core-collapse SNRs explode provides a large target density to boost hadronic gamma-rays. SNR G39.2-0.3 is one of the brightest SNR in infrared wavelengths, and its broad multi-wavelength coverage allows detailed modeling of its radiation from radio to high energies. We reanalyzed the Fermi-LAT data on this region and compare it with new radio observations from the MWISP survey. The modeling of the spectral energy distribution from radio to GeV energies favors an hadronic origin of the gamma-ray emission and constrains the SNR magnetic field to be at least ~100 uG. Despite the large magnetic field, the present acceleration of protons seems to be limited to ~10 GeV, which points to a drastic slow down of the shock velocity due to the dense wall traced by the CO observations, surrounding the remnant. Further investigation of the gamma-ray spectral shape points to a dynamically old remnant subjected to severe escape of CRs and a decrease of acceleration efficiency. The low-energy peak of the gamma-ray spectrum also suggests that that the composition of accelerated particles might be enriched by heavy nuclei which is certainly expected for a core-collapse SNR. Alternatively, the contribution of the compressed pre-existing Galactic cosmic rays is discussed, which is, however, found to not likely be the dominant process for gamma-ray production.

Kinematic distinction of the two subpopulations of X-ray pulsars

Pri?egen, M.

The population of Be/X-ray binaries shows strong evidence of bimodality, especially in the spin period of neutron stars. Several physical mechanisms may produce this bimodality. The most favored candidate mechanisms are two distinct supernova channels or different accretion modes of the neutron stars in Be/X-ray binaries. Investigating the kinematics of these systems may provide some additional insight into the physics of this bimodality. If the two Be/X-ray binary subpopulations arise from two distinct supernova types, then the two subpopulations should have different peculiar (systemic) velocities. This can be tested either directly, by measuring the velocity of the system, or indirectly, by measuring the position of the system with respect to its birthplace. Using the most recent Gaia dataset and the newest catalogs of Small Magellanic Cloud (SMC) star clusters, we analyzed the tangential peculiar velocities of Be/X-ray binaries in the Galaxy and the positions of Be/X-ray binaries in the SMC. We used the distance of the system from the nearest young star cluster as a proxy to the tangential velocity of the system. We applied statistical testing to investigate whether the two subpopulations that are divided by the spin of the neutron star are also kinematically distinct. There is evidence that the two subpopulations are indeed kinematically distinct. However, the tangential peculiar velocities of the two subpopulations are the reverse from what is expected from the distinct supernova channel hypothesis. We find some marginal evidence ($p \approx$~0.005) that the Galactic Be/X-ray binaries from the short-spin subpopulation have systematically higher peculiar velocities than the systems from the long-spin subpopulation. The same effect, but weaker, is also recovered for the SMC Be/X-ray binaries for all considered cluster catalogs.

Detailed abundances in a sample of very metal poor stars

Unevolved metal poor stars are the witness of the early evolution of the Galaxy. The determination of their detailed chemical composition is an important tool to understand the chemical history of our Galaxy. The study of their chemical composition can also be used to constrain the nucleosynthesis of the first generation of supernovae that enriched the interstellar medium. The aim is to observe a sample of extremely metal poor stars (EMP stars) candidates selected from SDSS DR12 release and determine their chemical composition. We obtained high resolution spectra of a sample of five stars using HDS on Subaru telescope and used standard 1D models to compute the abundances. The stars we analysed have a metallicity [Fe/H] between -3.50 dex and -4.25 dex . We confirm that the five metal poor candidates selected from low resolution spectra are very metal poor. We present, the discovery of a new ultra metal-poor star (UMP star) with a metallicity of [Fe/H]= -4.25 dex (SDSS~J1050032.34$-$241009.7). We measured in this star an upper limit of lithium ( log(Li/H) <= 2.0. We found that the 4 most metal poor stars of our sample have a lower lithium abundance than the Spite plateau lithium value. We obtain upper limits for carbon in the sample of stars. None of them belong to the high carbon band. We measured abundances of Mg and Ca in most of the stars and found three new alpha-poor stars.

Determining the Hubble Constant without the Sound Horizon Scale: Measurements from CMB Lensing

Baxter, Eric J.; Sherwin, Blake D.

Measurements of the Hubble constant, $H_0$, from the cosmic distance ladder are currently in tension with the value inferred from Planck observations of the CMB and other high redshift datasets if a flat $\Lambda$CDM cosmological model is assumed. One of the few promising theoretical resolutions of this tension is to invoke new physics that changes the sound horizon scale in the early universe; this can bring CMB and BAO constraints on $H_0$ into better agreement with local measurements. In this paper, we discuss how a measurement of the Hubble constant can be made from the CMB without using information from the sound horizon scale, $r_s$. In particular, we show how measurements of the CMB lensing power spectrum can be used to place interesting constraints on $H_0$ when combined with measurements of either supernovae or galaxy weak lensing, which constrain the matter density parameter. The constraints arise from the sensitivity of the CMB lensing power spectrum to the horizon scale at matter-radiation equality (in projection); this scale could have a different dependence on new physics than the sound horizon. From an analysis of current CMB lensing data from Planck and Pantheon supernovae with conservative external priors, we derive an $r_s$-independent constraint of $H_0 = 73.5\pm 5.3$ km/s/Mpc. Forecasts for future CMB surveys indicate that improving constraints beyond an error of $\sigma(H_0) = 3$ km/s/Mpc will be difficult with CMB lensing, although applying similar methods to the galaxy power spectrum may allow for further improvements.

Systematic investigation of the effect of 56Ni mixing in the early photospheric velocity evolution of stripped-envelope supernovae

Mixing of 56Ni, whose nuclear decay energy is a major luminosity source in stripped-envelope supernovae, is known to affect the observational properties of stripped-envelope supernovae such as light-curve and color evolution. Here we systematically investigate the effect of 56Ni mixing on the photospheric velocity evolution in stripped-envelope supernovae. We show that 56Ni mixing significantly affects the early photospheric velocity evolution. The photospheric velocity, which is often used to constrain the ejecta mass and explosion energy, significantly varies by just changing the degree of 56Ni mixing. In addition, the models with a small degree of 56Ni mixing show a flattening in the early photospheric velocity evolution, while the fully mixed models show a monotonic decrease. The velocity flattening appears in both helium and carbon+oxygen progenitor explosions with a variety of ejecta mass, explosion energy, and 56Ni mass. Some stripped-envelope supernovae with early photospheric velocity information do show such a flattening. We find that Type Ib SN 2007Y, which has early photospheric velocity information, has a signature of a moderate degree of 56Ni mixing in the photospheric velocity evolution and about a half of the ejecta is mixed in it. The immediate spectroscopic follow-up observations of stripped-envelope supernovae shortly after the explosion providing the early photospheric evolution give an important clue to constrain 56Ni mixing in the ejecta.

A new class of Roche lobe-filling hot subdwarf binaries

We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD.

Constraints and cosmography of $\Lambda$CDM in presence of viscosity

In this work, we study two scenarios of the Universe filled by a perfect fluid following the traditional dark energy and a viscous fluid as dark matter. In this sense, we explore the simplest case for the viscosity in the Eckart formalism, a constant, and then, a polynomial function of the redshift. We constrain the phase-space of the model parameters by performing a Bayesian analysis based on Markov Chain Monte Carlo method and using the latest data of the Hubble parameter (OHD), Type Ia Supernovae (SNIa) and Strong Lensing Systems. The first two samples cover the region $0.01

Quantum Black Hole Seismology II: Applications to Astrophysical Black Holes

Oshita, Naritaka; Tsuna, Daichi; Afshordi, Niayesh

With the advent of gravitational wave astronomy, searching for gravitational wave echoes from black holes (BHs) is becoming an interesting probe of their quantum nature near their horizons. Newborn BHs may be strong emitters of echoes, as they accompany large perturbations in the surrounding spacetime upon formation. Utilizing the Quantum Black Hole Seismology framework (Oshita et al. 2020), we study the expected echoes upon BH formation resulting from neutron star mergers and failed supernovae. For BH remnants from neutron star mergers, we evaluate the consistency of these models with the recent claim on the existence of echoes following the neutron star merger event GW170817. We find that the claimed echoes in GW170817, if real, suggest that overtones contribute a significant amount of energy in the ringdown of the remnant BH. We finally discuss the detectability of echoes from failed supernovae by second and third-generation gravitational wave detectors, and find that current (future) detectors constrain physical reflectivity models for events occurring within a few Mpc (a few x 10 Mpc). Detecting such echo signals may significantly constrain the maximum mass and equation of state of neutron stars.

Intermediate Luminosity Red Transients by Black Holes Born from Erupting Massive Stars

We consider black hole formation in failed supernovae when a dense circumstellar medium (CSM) is present around the massive star progenitor. By utilizing radiation hydrodynamical simulations, we calculate the mass ejection of blue supergiants and Wolf-Rayet stars in the collapsing phase and the radiative shock occurring between the ejecta and the ambient CSM. We find that the resultant emission is redder and dimmer than normal supernovae (bolometric luminosity of $\sim 10^{40}-10^{41}\ {\rm erg\ s^{-1}}$, effective temperature of $\sim 5\times 10^3$ K, and timescale of 10-100 days) and shows a characteristic power-law decay, which may comprise a fraction of intermediate luminosity red transients (ILRTs) including AT 2017be. In addition to searching for the progenitor star in the archival data, we encourage X-ray follow-up observations of such ILRTs $\sim$ 1-10 yr after the collapse, targeting the fallback accretion disk.

The low rate of Galactic pevatrons

Cristofari, P.; Blasi, P.; Amato, E.

Although supernova remnants remain the main suspects as sources of Galactic cosmic rays up to the knee, the supernova paradigm still has many loose ends. The weakest point in this construction is the possibility that individual supernova remnants can accelerate particles to the rigidity of the knee, $\sim 10^{6}$ GV. This scenario heavily relies upon the possibility to excite current driven non-resonant hybrid modes while the shock is still at the beginning of the Sedov phase. These modes can enhance the rate of particle scattering thereby leading to potentially very-high maximum energies. Here we calculate the spectrum of particles released into the interstellar medium from the remnants of different types of supernovae. We find that only the remnants of very powerful, rare core-collapse supernova explosions can accelerate light elements such as hydrogen and helium nuclei, to the knee rigidity, and that the local spectrum of cosmic rays directly constrains the rate of such events, if they are also source of PeV cosmic rays. On the other hand, for remnants of typical core-collapse supernova explosions, the Sedov phase is reached at late times, when the maximum energy is too low and the spectrum at very-high energies is very steep, being mostly produced during the ejecta dominated phase. For typical thermonuclear explosions, resulting in type Ia supernovae, we confirm previous findings that these objects can only produce cosmic rays up to $\lesssim 10^{5}$ GeV. The implications for the overall cosmic ray spectrum observed at the Earth and for the detection of PeVatrons by future gamma-ray observatories are discussed.

Production of secondary particles in heavy nuclei interactions in supernova remnants

Depending on their type, supernova remnants may have ejecta material with high abundance of heavy elements such as carbon or oxygen. In addition, core-collapse supernovae explode in the wind material of their progenitor star that may also have a high abundance of heavy elements. Hadronic collisions in these enriched media spawn the production of gamma rays, neutrinos, and secondary electrons with spectra that cannot be scaled from those calculated for pp collisions, potentially leading to erroneous results. We used Monte-Carlo event generators to calculate the differential production rate of particles such as gamma rays, neutrinos, and secondary electrons for H, He, C, and O nuclei as projectiles and as target material. The cross sections and the multiplicity spectra are separately computed for each of the 16 combinations of projectile and target. We describe characteristic effects of heavy nuclei in the shape and normalization of the spectra of various particles produced.

Binary population synthesis with probabilistic remnant mass and kick prescriptions

We report on the impact of a probabilistic prescription for compact remnant masses and spins on massive binary population synthesis. We find that this prescription populates the putative mass gap between neutron stars and black holes with low-mass black holes. However, evolutionary effects reduce the number of X-ray binary candidates with low-mass black holes, consistent with the dearth of such systems in the observed sample. We further find that this prescription is consistent with the formation of heavier binary neutron stars such as GW190425, but over-predicts the masses of Galactic double neutron stars. The revised natal kicks, particularly increased ultra-stripped supernova kicks, do not directly explain the observed Galactic double neutron star orbital period--eccentricity distribution. Finally, this prescription allows for the formation of systems similar to the recently discovered extreme mass ratio binary GW190814, but only if we allow for the survival of binaries in which the common envelope is initiated by a donor crossing the Hertzsprung gap, contrary to our standard model.

Discovery of recombining plasma associated with the candidate supernova remnant G189.6+3.3 with Suzaku

Yamauchi, Shigeo; Oya, Moe; Nobukawa, Kumiko K.; Pannuti, Thomas G.

We present the results of an X-ray spectral analysis of the northeast region of the candidate supernova remnant G189.6+3.3 with Suzaku. K-shell lines from highly ionized Ne, Mg, Si, and S were detected in the spectrum for the first time. In addition, a radiative recombining continuum (RRC) from He-like Si was clearly seen near 2.5 keV. This detection of an RRC reveals for the first time that G189.6+3.3 possesses an X-ray-emitting recombining plasma (RP). The extracted X-ray spectrum in the 0.6-10.0 keV energy band is well fitted with a model consisting of a collisional ionization equilibrium plasma component (associated with the interstellar medium) and an RP component (associated with the ejecta). The spectral feature shows that G189.6+3.3 is most likely to be a middle-aged SNR with an RP.

Strongly lensed SN Refsdal: refining time delays based on the supernova explosion models

Baklanov, Petr; Lyskova, Natalia; Blinnikov, Sergei; Nomoto, Ken'ichi

We explore the unique supernova (SN) "Refsdal" - the first discovered gravitationally lensed SN with multiple images. A large magnification provided by the galactic-scale lens, augmented by the cluster lens, gave us the opportunity to perform a detailed modelling of a distant SN at $z \simeq 1.5$. % We present results of radiation hydrodynamics modelling of \snref. According to our calculations, the \snref progenitor is likely to be a more massive and energetic version of \sna, i.e. a blue supergiant star with a radius of ${\rm R} = 50 R_{\odot}$ before an explosion and a total mass of ${\rm M}= 26 M_{\odot}$. The total energy release is ${\rm E}=5\times10^{51}$ erg. Reconstruction of SN light curves allowed us to obtain time delays and magnifications for the images S2-S4 relative to S1 with higher accuracy than previous estimates of Rodney et al. We estimate the Hubble constant $H_0 = 68.6^{+13.6}_{-9.7}$ km s$^{-1}$ Mpc$^{-1}$ via re-scaling the time delays predicted by different lens models to match the values obtained in this work. With more photometric data on the fifth image SX, we will be able to further refine the time delay and magnification estimates for SX and obtain competitive constraints on $H_0$.

The Role of Stellar Feedback in the Chemical Evolution of a Low Mass Dwarf Galaxy

Emerick, Andrew; Bryan, Greg L.; Mac Low, Mordecai-Mark

We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust grains, ionization and associated heating, Lyman-Werner (LW) dissociation of H$_2$, and radiation pressure; and winds from massive main sequence (neglecting their energy input) and asymptotic giant branch (AGB) stars. Radiative transfer is done by ray tracing. We consider the effects each of these processes have on regulating the star formation rate, global properties, multi-phase interstellar medium (ISM), and driving of galactic winds. We follow individual metal species from distinct nucleosynthetic enrichment channels (AGB winds, massive star stellar winds, core collapse and Type Ia supernovae) and pay particular attention to how these feedback processes regulate metal mixing in the ISM, the metal content of outflows, and the stellar abundance patterns in our galaxy. We find that---for a low-metallicity, low-mass dwarf galaxy ---stellar radiation, particularly ionizing radiation and LW radiation, are important sources of stellar feedback whose effects dominate over photoelectric heating and HI radiation pressure. However, feedback is coupled non-linearly, and the inclusion or exclusion of each process produces non-negligible effects. We find strong variations with: the star formation history; the ejection fractions of metals, mass, and energy; and the distribution of elements from different nucleosynthetic sources in both the gas and stars.

The Spectacular Ultraviolet Flash From the Peculiar Type Ia Supernova 2019yvq

Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

Lifting the core-collapse supernova bounds on keV-mass sterile neutrinos

Suliga, Anna M.; Tamborra, Irene; Wu, Meng-Ru

We explore the energy and entropy transport as well as the lepton number variation induced from the mixing between electron and sterile neutrinos with keV mass in the supernova core. We develop a radial- and time-dependent treatment of the sterile-electron neutrino mixing, by including ordinary matter effects, reconversions between sterile and electron antineutrinos, as well as the collisional production of sterile particles for the first time. The dynamical feedback due to the production of sterile particles on the composition and thermodynamic properties of the core only leads to major implications for the supernova physics for large mixing angles ($\sin^2 2 \theta > 10^{-10}$). Our findings suggest that a self-consistent appraisal of the electron-sterile conversion physics in the supernova core would relax the bounds on the sterile neutrino mixing parameters reported in the literature for mixing angles smaller than $10^{-6}$, leaving the parameter space of the mass and mixing angles of sterile neutrinos relevant to dark matter searches unconstrained by supernovae.

Early recombination as a solution to the $H_0$ tension

Sekiguchi, Toyokazu; Takahashi, Tomo

We show that the $H_0$ tension can be resolved by making recombination earlier, keeping the fit to cosmic microwave background (CMB) data almost intact. We provide a suite of general necessary conditions to give a good fit to CMB data while realizing a high value of $H_0$ suggested by local measurements. As a concrete example for a successful scenario with early recombination, we demonstrate that a model with time-varying $m_e$ can indeed satisfy all the conditions. We further show that such a model can also be well fitted to low-$z$ distance measurements of baryon acoustic oscillation (BAO) and type-Ia supernovae (SNeIa) with a simple extension of the model. Time-varying $m_e$ in the framework of $\Omega_k\Lambda$CDM is found to be a sufficient and excellent example as a solution to the $H_0$ tension, yielding $H_0=72.3_{-2.8} ^{+2.7}\,$km/sec/Mpc from the combination of CMB, BAO and SNeIa data even without incorporating any direct local $H_0$ measurements. Apart from the $H_0$ tension, this model is also favored from the viewpoint of the CMB lensing anomaly.

The Carnegie Supernova Project II. Observations of SN 2014ab reveal a 2010jl-like SN IIn with pre-existing dust?

We present optical and near-infrared photometry and spectroscopy of the Type IIn supernova (SN) 2014ab, obtained by the Carnegie Supernova Project II (CSP-II) and initiated immediately after its optical discovery. We also present mid-infrared photometry obtained by the Wide-field Infrared Survey Explorer (WISE) satellite extending from 56 days prior to the optical discovery to over 1600 days. The light curve of SN 2014ab evolves slowly, while the spectra exhibit strong emission features produced from the interaction between rapidly expanding ejecta and dense circumstellar matter. The light curve and spectral properties are very similar to those of SN 2010jl. The estimated mass-loss rate of the progenitor of SN 2014ab is of the order of 0.1 Msun/yr under the assumption of spherically symmetric circumstellar matter and steady mass loss. Although the mid-infrared luminosity increases due to emission from dust, which is characterized by a blackbody temperature close to the dust evaporation temperature (~ 2000 K), no clear signatures of in situ dust formation within the cold dense shell located behind the forward shock are observed in SN 2014ab in early phases. Mid-infrared emission of SN 2014ab may originate from pre-existing dust located within dense circumstellar matter that is heated by the SN shock or shock-driven radiation. Finally, for the benefit of the community, we also present in an Appendix five near-infrared spectra of SN 2010jl obtained between 450 to 1300 days post discovery.

Testing Modified Gravity theory (MOG) with Type Ia Supernovae, Cosmic Chronometers and Baryon Acoustic Oscillations

We analyse the MOdified Gravity (MOG) theory, proposed by Moffat, in a cosmological context. We use data from Type Ia Supernovae (SNe Ia), Baryon Acoustic Oscillations (BAO) and Cosmic Chronometers (CC) to test MOG predictions. For this, we perform $\chi^2$ tests considering fixed values of $H_0$ and $V_G$, the self-interaction potential of one of the scalar fields in the theory. Our results show that the MOG theory is in agreement with all data sets for some particular values of $H_0$ and $V_G$, being the BAO data set the most powerful tool to test MOG predictions, due to its constraining power.

Element Stratification in the Middle-Aged Type Ia Supernova Remnant G344.7-0.1

Despite their importance, a detailed understanding of Type Ia supernovae (SNe Ia) remains elusive. X-ray measurements of the element distributions in supernova remnants (SNRs) offer important clues for understanding the explosion and nucleosynthesis mechanisms for SNe Ia. However, it is challenging to observe the entire ejecta mass in X-rays for young SNRs, because the central ejecta may not have been heated by the reverse shock yet. Here we present over 200 kilosecond Chandra observations of the Type Ia SNR G344.7-0.1, whose age is old enough for the reverse shock to have reached the SNR center, providing an opportunity to investigate the distribution of the entire ejecta mass. We reveal a clear stratification of heavy elements with a centrally peaked distribution of the Fe ejecta surrounded by intermediate-mass elements (IMEs: Si, S, Ar Ca) with an arc-like structure. The centroid energy of the Fe K emission is marginally lower in the central Fe-rich region than in the outer IME-rich regions, suggesting that the Fe ejecta were shock-heated more recently. These results are consistent with the prediction for standard SN Ia models, where the heavier elements are synthesized in the interior of an exploding white dwarf. We find, however, that the peak location of the Fe K emission is slightly offset to the west with respect to the geometric center of the SNR. This apparent asymmetry is likely due to the inhomogeneous density distribution of the ambient medium, consistent with our radio observations of the ambient molecular and neutral gas.

On the origin of GW190425

The LIGO/Virgo collaborations recently announced the detection of a likely binary neutron star merger, GW190425. The total mass of GW190425 is significantly larger than the masses of Galactic double neutron stars known through radio astronomy. This suggests that GW190425 formed differently from Galactic double neutron stars. We hypothesize that GW190425 formed via unstable "case BB" mass transfer. According to this hypothesis, the progenitor of GW190425 was a binary consisting of a neutron star and a ${\sim} 4-5 {M_\odot}$ helium star, which underwent a common-envelope process. Following the supernovae of the helium star core, a tight, eccentric, double neutron star was formed, which merged in ${\lesssim}$ 10 Myr. The helium star progenitor may explain the unusually large mass of GW190425, while the short time to merger may explain why we do not see similar systems in radio. In order to test this hypothesis, we measure the eccentricity of GW190425 using publicly available LIGO/Virgo data. We constrain the eccentricity at 10 Hz to be $e \leq 0.007$ with $90\%$ confidence. This result provides no evidence for or against the unstable mass transfer scenario because the binary is likely to have circularized to $e\lesssim10^{-4}$ by the time it entered the LIGO/Virgo band. Future detectors operating in lower frequency bands will enable us to discern the formation channel of mergers similar to GW190425 using eccentricity measurements.

Defying the laws of Gravity I: model-independent reconstruction of the Universe expansion from growth data

L'Huillier, Benjamin; Shafieloo, Arman; Polarski, David; Starobinsky, Alexei A.

Using redshift space distortion data, we perform model-independent reconstructions of the growth history of matter inhomogeneity in the expanding Universe using two methods: crossing statistics and Gaussian processes. We then reconstruct the corresponding history of the Universe background expansion and fit it to type Ia supernovae data, putting constraints on $(\Omega_\mathrm{m,0},\sigma_{8,0})$. The results obtained are consistent with the concordance flat-$\Lambda$CDM model and General Relativity as the gravity theory given the current quality of the inhomogeneity growth data.

Neutrinos from the cosmic noon: a probe of the cosmic star formation history

Riya; Rentala, Vikram

Multiple astrophysical probes of the cosmic star formation history yield widely different inferences of this rate at redshifts z > 1. While all probes seem to indicate a period of peak star formation known as the cosmic noon between 1.5 < z < 3, the detailed inferences from these probes are in disagreement. In particular, the magnitude of the peak star formation rate density indicated by H-alpha data is higher by a factor of ~ 4 compared to the magnitude of the peak indicated by UV/IR data. In this work, we explore the potential of future measurements of the diffuse supernova neutrino background at the Hyper-Kamiokande (HK) experiment to resolve the discrepancy and help pin down the magnitude of the peak cosmic star formation rate. We find that, depending upon the cosmic core-collapse supernova neutrino spectrum, HK loaded with 0.1% Gadolinium by mass has the potential to discriminate between the different star formation histories with between 1.6-20 years of data collection.

Boosting MonteCarlo sampling with a non-Gaussian fit

We propose a new method, called MonteCarlo Posterior Fit, to boost the MonteCarlo sampling of likelihood (posterior) functions. The idea is to approximate the posterior function by an analytical multidimensional non-Gaussian fit. The many free parameters of this fit can be obtained by a smaller sampling than is needed to derive the full numerical posterior. In the examples that we consider, based on supernovae and cosmic microwave background data, we find that one needs an order of magnitude smaller sampling than in the standard algorithms to achieve comparable precision. This method can be applied to a variety of situations and is expected to significantly improve the performance of the MonteCarlo routines in all the cases in which sampling is very time-consuming. Finally, it can also be applied to Fisher matrix forecasts, and can help solve various limitations of the standard approach.

A Search for Supernova Light Echoes in NGC 6946 with SITELLE

Radica, M. C.; Welch, D.; Rousseau-Nepton, L.

We present the analysis of four hours of spectroscopic observations of NGC 6946 with the SITELLE Imaging Fourier Transform Spectrometer on the Canada-France-Hawaii Telescope, acquired to search for supernova light echoes from its ten modern supernovae. We develop a novel spectroscopic search method: identifying negatively sloped continua in the narrow-band SN3 filter as candidate highly-broadened P-Cygni profiles in the H$\alpha$ line, which would be characteristic of the spectra of supernovae ejecta. We test our methodology by looking for light echoes from any of the ten supernovae observed in NGC 6946 in the past 100 years. We find no evidence of light echoes above the survey surface brightness limit of 1$\times$10$^{-15}$erg/s/cm$^2$/arcsec$^2$.

Neutron Star Extreme Matter Observatory: A kilohertz-band gravitational-wave detector in the global network

Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

Accurate and precision Cosmology with redshift unknown gravitational wave sources

Gravitational waves can provide an accurate measurement of the luminosity distance to the source, but cannot provide the source redshift unless the degeneracy between mass and redshift can be broken. This makes it essential to infer the redshift of the source independently to measure the expansion history of the Universe. We show that by exploiting the clustering scale of the gravitational wave sources with galaxies of known redshift, we can infer the expansion history from redshift unknown gravitational wave sources. Using gravitational wave sources with unknown redshift that are detectable from the network of gravitational wave detectors with Advanced LIGO design sensitivity, we will be able to obtain accurate and precise measurements of the local Hubble constant, the expansion history of the universe, and the gravitational wave bias parameter, which captures the distribution of gravitational wave sources with respect to the redshift tracer distribution. This technique is not only limited to the low redshift gravitational wave sources, but will be also applicable to the high redshift gravitational wave sources detectable from Laser Interferometer Space Antenna (LISA), Cosmic Explorer (CE), and Einstein Telescope (ET). Moreover, this method is also going to be applicable to samples of supernovae and fast radio bursts with unknown or photometric redshifts.

INTEGRAL discovery of a burst with associated radio emission from the magnetar SGR 1935+2154

We report on INTEGRAL observations of the soft $\gamma$-ray repeater SGR 1935+2154 performed between 2020 April 28 and May 3. Several short bursts with fluence of $\sim10^{-7}-10^{-6}$ erg cm$^{-2}$ were detected by the IBIS instrument in the 20-200 keV range. The burst with the hardest spectrum, discovered and localized in real time by the INTEGRAL Burst Alert System, was spatially and temporally coincident with a short and very bright radio burst detected by the CHIME and STARE2 radio telescopes at 400-800 MHz and 1.4 GHz, respectively. Its lightcurve shows three narrow peaks separated by $\sim$29 ms time intervals, superimposed on a broad pulse lasting $\sim$0.6 s. The brightest peak had a delay of 6.5$\pm$1.0 ms with respect to the 1.4 GHz radio pulse (that coincides with the second and brightest component seen at lower frequencies). The burst spectrum, an exponentially cut-off power law with photon index $\Gamma=0.7_{-0.2}^{+0.4}$ and peak energy $E_p=65\pm5$ keV, is harder than those of the bursts usually observed from this and other magnetars. By the analysis of an expanding dust scattering ring seen in X-rays with the {\it Neil Gehrels Swift Observatory} XRT instrument, we derived a distance of 4.4$_{-1.3}^{+2.8}$ kpc for SGR 1935+2154, independent of its possible association with the supernova remnant G57.2+0.8. At this distance, the burst 20-200 keV fluence of $(6.1\pm 0.3)\times10^{-7}$ erg cm$^{-2}$ corresponds to an isotropic emitted energy of $\sim1.4\times10^{39}$ erg. This is the first burst with a radio counterpart observed from a soft $\gamma$-ray repeater and it strongly supports models based on magnetars that have been proposed for extragalactic fast radio bursts.

Radiative-transfer modeling of nebular-phase type II supernovae. Dependencies on progenitor and explosion properties

Dessart, Luc; Hillier, D. John

Nebular phase spectra of core-collapse supernovae (SNe) provide critical and unique information on the progenitor massive star and its explosion. We present a set of 1-D steady-state non-local thermodynamic equilibrium radiative transfer calculations of type II SNe at 300d after explosion. Guided by results for a large set of stellar evolution simulations, we craft ejecta models for type II SNe from the explosion of a 12, 15, 20, and 25Msun star. The ejecta density structure and kinetic energy, the 56Ni mass, and the level of chemical mixing are parametrized. Our model spectra are sensitive to the adopted line Doppler width, a phenomenon we associate with the overlap of FeII and OI lines with Lyalpha and Lybeta. Our spectra show a strong sensitivity to 56Ni mixing since it determines where decay power is absorbed. Even at 300d after explosion, the H-rich layers reprocess the radiation from the inner metal rich layers. In a given progenitor model, variations in 56Ni mass and distribution impact the ejecta ionization, which can modulate the strength of all lines. Such ionization shifts can quench CaII line emission. In our set of models, the OI6300 doublet strength is the most robust signature of progenitor mass. However, we emphasize that convective shell merging in the progenitor massive star interior can pollute the O-rich shell with Ca, which will weaken the OI6300 doublet flux in the resulting nebular SN II spectrum. This process may occur in Nature, with a greater occurrence in higher mass progenitors, and may explain in part the preponderance of progenitor masses below 17Msun inferred from nebular spectra.

Evaluating the Calibration of SN Ia Anchor Datasets with a Bayesian Hierarchical Model

Inter-survey calibration remains an important systematic uncertainty in cosmological studies using type Ia supernova (SNe Ia). Ideally, each survey would measure its system throughputs, for instance with bandpass measurements combined with observations of well-characterized spectrophotometric standard stars; however, many important nearby-SN surveys have not done this. We recalibrate these surveys by tying their tertiary survey stars to Pan-STARRS1 g, r, and i, and SDSS/CSP u. This improves upon previous recalibration efforts by taking the spatially variable zeropoints of each telescope/camera into account, and applying improved color transformations in the surveys' natural instrumental photometric systems. Our analysis uses a global hierarchical model of the data which produces a covariance matrix of magnitude offsets and bandpass shifts, quantifying and reducing the systematic uncertainties in the calibration. We call our method CROSS-CALIBration with a Uniform Reanalysis (X-CALIBUR). This approach gains not only from a sophisticated analysis, but also from simply tying our calibration to more color calibrators, rather than just the one color calibrator (BD+17 4708) as many previous efforts have done. The results presented here have the potential to help understand and improve calibration uncertainties upcoming SN Ia cosmological analyses.

New Early Dark Energy

Niedermann, Florian; Sloth, Martin S.

New measurements of the expansion rate of the universe have plunged the standard model of cosmology into a severe crisis. In this letter we propose a simple resolution to the problem. We propose that a first order phase transition in a dark sector in the early universe, before recombination, can resolve the problem. This will lead to a short phase of a New Early Dark Energy (NEDE) component and can explain the observations. Fitting our model to measurements of the Cosmic Microwave Background, Baryonic Acoustic Oscillations, and supernovae yields a significant improvement of the best-fit compared with the standard cosmological model without NEDE. We find the mean value of the present Hubble parameter in the NEDE model to be $H_0=71.4 \pm 1.0 ~\textrm{km}\, \textrm{s}^{-1}\, \textrm{Mpc}^{-1}$ ($68\, \%$ C.L.).

Resolving the Hubble Tension with New Early Dark Energy

Niedermann, Florian; Sloth, Martin S.

New Early Dark Energy (NEDE) is a component of vacuum energy at the electron volt scale, which decays in a first order phase transition shortly before recombination [arXiv:1910.10739]. The NEDE component has the potential to resolve the tension between recent local measurements of the expansion rate of the universe using supernovae (SN) data, and the expansion rate inferred from the early universe through measurements of the Cosmic Microwave Background (CMB) when assuming $\Lambda$CDM. We discuss in depth the two-scalar field model of the NEDE phase transition including the process of bubble percolation, collision and coalescence. We also estimate the gravitational wave signal produced during the collision phase and argue that it can be searched for using pulsar timing arrays. In a second step, we construct an effective cosmological model, which describes the phase transition as an instantaneous process, and derive the covariant equations that match perturbations across the transition surface. Fitting the cosmological model to CMB, Baryonic Acoustic Oscillations and SN data, we report $H_0 = 69.6^{+1.0}_{-1.3} \, \textrm{km}\, \textrm{s}^{-1}\, \textrm{Mpc}^{-1}$ $(68 \%$ C.L.) without the local measurement of the Hubble parameter, bringing the tension down to $2.5\, \sigma$. Including the local input, we find $H_0 = 71.4 \pm 1.0 \, \textrm{km}\, \textrm{s}^{-1}\, \textrm{Mpc}^{-1}$ $(68 \%$ C.L.) and strong evidence for a non-vanishing NEDE component with a $\simeq 4\, \sigma$ significance.

Nickel clumps and light curve bumps: Light curves and spectra of models with macroscopic $^{56}$Ni clumps in the outer ejecta

Magee, M. R.; Maguire, K.

An excess of flux in the early light curves of type Ia supernovae has been observed in a few cases. Multiple scenarios have been proposed to explain this. Recently, it has been shown that for at least one object (SN$~$2018oh) the excess emission observed could be the result of a large-scale clump of $^{56}$Ni in the outer ejecta of $\sim$0.03$~M_{\rm{\odot}}$. We present model light curves and spectra for ejecta profiles containing $^{56}$Ni clumps of varying masses (0.01, 0.02, 0.03, and 0.04$~M_{\rm{\odot}}$) and shapes. We find that even for our lowest mass $^{56}$Ni clump, an increase of $\gtrsim$2 magnitudes is produced in the bolometric light curve at one day after explosion, relative to models without a $^{56}$Ni clump. We show that the colour evolution of models with a $^{56}$Ni clump differs significantly from those without, and shows a colour inversion similar to some double-detonation explosion models. Furthermore, spectra of our $^{56}$Ni clump models show that strong suppression of flux between $\sim$3$~$700 $-$ 4$~$000$~$Ang. close to maximum light appears to be a generic feature for this class of model. Comparing our models to observations of SNe$~$2017cbv and 2018oh, we show that a $^{56}$Ni clump of 0.02 $-$ 0.04$~M_{\rm{\odot}}$ can match shapes of the early light curve bumps, however the colour and spectral evolution are in disagreement. This would indicate that an alternative origin for the flux excess is necessary. In addition, based on existing explosion scenarios, producing a large-scale, macroscopic $^{56}$Ni clump in the outer ejecta as required to match the light curve shape, without the presence of additional short-lived radioactive material, may prove challenging. Given that only a small amount of $^{56}$Ni in the outer ejecta is required to produce a bump in the light curve, such a large clump in the outer ejecta must be rare, if it were to occur at all.

The lowest frequency Fast Radio Bursts: Sardinia Radio Telescope detection of the periodic FRB 180916 at 328 MHz

We report on the lowest-frequency detection to date of three bursts from the fast radio burst FRB 180916, observed at 328 MHz with the Sardinia Radio Telescope (SRT). The SRT observed the periodic repeater FRB 180916 for five days from 2020 February 20 to 24 during a time interval of active radio bursting, and detected the three bursts during the first hour of observations; no more bursts were detected during the remaining ~ 30 hours. Simultaneous SRT observations at 1548 MHz did not detect any bursts. Burst fluences are in the range 37 to 13 Jy ms. No relevant scattering is observed for these bursts. We also present the results of the multi-wavelength campaign we performed on FRB 180916, during the five days of the active window. Simultaneously with the SRT observations, others with different time spans were performed with the Northern Cross at 408 MHz, with XMM-Newton, NICER, INTEGRAL, AGILE, and with the TNG and two optical telescopes in Asiago, which are equipped with fast photometers. XMM-Newton obtained data simultaneously with the three bursts detected by the SRT, and determined a luminosity upper limit in the 0.3-10 keV energy range of ~$10^{45}$ erg/s for the burst emission. AGILE obtained data simultaneously with the first burst and determined a fluence upper limit in the MeV range for millisecond timescales of $ 10^{-8}$ erg cm$^{-2}$.Our results show that absorption from the circumburst medium does not significantly affect the emission from FRB 180916, thus limiting the possible presence of a superluminous supernova around the source, and indicate that a cutoff for the bursting mechanism, if present, must be at lower frequencies. Our multi-wavelength campaign sensitively constrains the broadband emission from FRB 180916, and provides the best limits so far for the electromagnetic response to the radio bursting of this remarkable source of fast radio bursts.

Supernova Triggers for End-Devonian Extinctions?

The Late Devonian was a protracted period of low speciation resulting in biodiversity decline, culminating in extinction events near the Devonian-Carboniferous boundary. Recent evidence indicates that the final extinction event may have coincided with a dramatic drop in stratospheric ozone, possibly due to a global temperature rise. Here we study an alternative possible cause for the postulated ozone drop: a nearby supernova explosion that could inflict damage by accelerating cosmic rays that can deliver ionizing radiation for up to $\sim 100$ kyr. We therefore propose that end-Devonian extinction was triggered by one or more supernova explosions at $\sim 20 \ \rm pc$, somewhat beyond the ``kill distance'' that would have precipitated a full mass extinction. Nearby supernovae are likely due to core-collapses of massive stars in clusters in the thin Galactic disk in which the Sun resides. Detecting any of the long-lived radioisotopes \sm146, \u235 or \pu244 in one or more end-Devonian extinction strata would confirm a supernova origin, point to the core-collapse explosion of a massive star, and probe supernova nucleosythesis. Other possible tests of the supernova hypothesis are discussed.

Dalek -- a deep-learning emulator for TARDIS

Supernova spectral time series contain a wealth of information about the progenitor and explosion process of these energetic events. The modeling of these data requires the exploration of very high dimensional posterior probabilities with expensive radiative transfer codes. Even modest parametrizations of supernovae contain more than ten parameters and a detailed exploration demands at least several million function evaluations. Physically realistic models require at least tens of CPU minutes per evaluation putting a detailed reconstruction of the explosion out of reach of traditional methodology. The advent of widely available libraries for the training of neural networks combined with their ability to approximate almost arbitrary functions with high precision allows for a new approach to this problem. Instead of evaluating the radiative transfer model itself, one can build a neural network proxy trained on the simulations but evaluating orders of magnitude faster. Such a framework is called an emulator or surrogate model. In this work, we present an emulator for the TARDIS supernova radiative transfer code applied to Type Ia supernova spectra. We show that we can train an emulator for this problem given a modest training set of a hundred thousand spectra (easily calculable on modern supercomputers). The results show an accuracy on the percent level (that are dominated by the Monte Carlo nature of TARDIS and not the emulator) with a speedup of several orders of magnitude. This method has a much broader set of applications and is not limited to the presented problem.

Constraints on electron-scattering interpretation of XENON1T excess

Chigusa, So; Endo, Motoi; Kohri, Kazunori

Recently, the XENON1T experiment has observed an excess in the electronic recoil data in the recoil energy range of $1$-$7$ keV. One of the most favored new physics interpretations is electron scattering with a boosted particle with a velocity of $\sim 0.1$ and a mass of $\gtrsim 0.1\,\mathrm{MeV}$. If such a particle has a strong interaction with electrons, it may affect the standard scenario of cosmology or be observed at low-threshold direct detection experiments. We study various constraints, mainly focusing on those from the big-bang nucleosynthesis, supernova cooling, and direct detection experiments. We discuss the implication of these constraints on electron-scattering interpretation of the XENON1T excess.

Ultralight Scalar Decay and the Hubble Tension

Gonzalez, Mark; Hertzberg, Mark P.; Rompineve, Fabrizio

We examine whether the Hubble tension, the mismatch between early and late measurements of $H_0$, can be alleviated by ultralight scalar fields in the early universe, and we assess its plausibility within UV physics. Since their energy density needs to rapidly redshift away, we explore decays to massless fields around the era of matter-radiation equality. We highlight a concrete implementation of ultralight pseudo-scalars, axions, that decay to an abelian dark sector. This scenario circumvents major problems of other popular realizations of early universe scalar models in that it uses a regular scalar potential that is quadratic around the minimum, instead of the extreme fine-tuning of many existing models. The idea is that the scalar is initially frozen in its potential until $H\sim m$, then efficient energy transfer from the scalar to the massless field can occur shortly after the beginning of oscillations due to resonance. We introduce an effective fluid model which captures the transition from the frozen scalar phase to the radiation dark sector phase. We perform a fit to a combined Planck 2018, BAO, SH$_0$ES and Pantheon supernovae dataset and find that the model gives $H_0=69.69_{-0.83}^{+0.81}$ km/s/Mpc with $\Delta\chi^2 \approx -8$ compared to $\Lambda$CDM; while inclusions of other data sets may worsen the fit. Importantly, we find that large values of the coupling between fields is required for sufficiently rapid decay: For axion-gauge field models $\phi F\tilde{F}/\Lambda$ it requires $\Lambda\lesssim f/80$, where $2\pi f$ is the field range. We find related conclusions for scalar-scalar models $\sim\phi\,\chi^2$ and for models that utilize perturbative decays. We conclude that these sorts of ultralight scalar models that purport to alleviate the Hubble tension, while being reasonable effective field theories, require features that are difficult to embed within UV physics.

The broad-line type Ic SN 2020bvc: signatures of an off-axis gamma-ray burst afterglow

Long-duration gamma-ray bursts (GRBs) are almost unequivocally associated with very energetic, broad-lined supernovae (SNe) of Type Ic-BL. While the gamma-ray emission is emitted in narrow jets, the SN emits radiation isotropically. Therefore, some SN Ic-BL not associated with GRBs have been hypothesized to arise from events with inner engines such as off-axis GRBs or choked jets. Here we present observations of the nearby ($d = 120$ Mpc) SN 2020bvc (ASAS-SN 20bs) which support this scenario. \textit{Swift} UVOT observations reveal an early decline (up to two days after explosion) while optical spectra classify it as a SN Ic-BL with very high expansion velocities ($\approx$ 70,000 km/s), similar to that found for the jet-cocoon emission in SN 2017iuk associated with GRB 171205A. Moreover, \textit{Swift} X-Ray Telescope and \textit{CXO} X-ray Observatory detected X-ray emission only three days after the SN and decaying onwards, which can be ascribed to an afterglow component. Cocoon and X-ray emission are both signatures of jet-powered GRBs. In the case of SN 2020bvc, we find that the jet is off axis (by $\approx$ 23 degrees), as also indicated by the lack of early ($\approx 1$ day) X-ray emission which explains why no coincident GRB was detected promptly or in archival data. These observations suggest that SN 2020bvc is the first orphan GRB detected through its associated SN emission.

Observable signatures of dark photons from supernovae

A dark photon is a well-motivated new particle which, as a component of an associated dark sector, could explain dark matter. One strong limit on dark photons arises from excessive cooling of supernovae. We point out that even at couplings where too few dark photons are produced in supernovae to violate the cooling bound, they can be observed directly through their decays. Supernovae produce dark photons which decay to positrons, giving a signal in the 511 keV annihilation line observed by SPI/INTEGRAL. Further, prompt gamma-ray emission by these decaying dark photons gives a signal for gamma-ray telescopes. Existing GRS observations of SN1987a already constrain this, and a future nearby SN could provide a detection. Finally, dark photon decays from extragalactic SN would produce a diffuse flux of gamma rays observable by detectors such as SMM and HEAO-1. Together these observations can probe dark photon couplings several orders of magnitude beyond current constraints for masses of roughly 1 - 100 MeV.

Statistical characterization and classification of astronomical transients with Machine Learning in the era of the Vera Rubin Survey Telescope

Vicedomini, M.; Brescia, M.; Cavuoti, S.; Longo, G.; Riccio, G.

Astronomy has entered the multi-messenger data era and Machine Learning has found widespread use in a large variety of applications. The exploitation of synoptic (multi-band and multi-epoch) surveys, like LSST (Large Synoptic Survey Telescope), requires an extensive use of automatic methods for data processing and interpretation. With data volumes in the petabyte domain, the discrimination of time-critical information has already exceeded the capabilities of human operators and crowds of scientists have extreme difficulty to manage such amounts of data in multi-dimensional domains. This work is focused on an analysis of critical aspects related to the approach, based on Machine Learning, to variable sky sources classification, with special care to the various types of Supernovae, one of the most important subjects of Time Domain Astronomy, due to their crucial role in Cosmology. The work is based on a test campaign performed on simulated data. The classification was carried out by comparing the performances among several Machine Learning algorithms on statistical parameters extracted from the light curves. The results make in evidence some critical aspects related to the data quality and their parameter space characterization, propaedeutic to the preparation of processing machinery for the real data exploitation in the incoming decade.

On the Distance of SGR 1935+2154 Associated with FRB 200428 and Hosted in SNR G57.2+0.8

Zhong, S. Q.; Dai, Z. G.; Zhang, H. M.; Deng, C. M.

Owing to the detection of an extremely bright fast radio burst (FRB) 200428 associated with a hard X-ray counterpart from the magnetar soft gamma-ray repeater (SGR) 1935+2154, the distance of SGR 1935+2154 potentially hosted in the supernova remnant (SNR) G57.2+0.8 can be revisited. Under the assumption that the SGR and the SNR are physically related, in this Letter, by investigating the dispersion measure (DM) of the FRB contributed by the foreground medium of our Galaxy and the local environments and combining with other observational constraints, we find that the distance of SGR 1935+2154 turns out to be $9.0\pm2.5\,$kpc and the SNR radius falls into $10$ to $18\,$pc since the local DM contribution is as low as $0-18\,$pc cm$^{-3}$. These results are basically consistent with the previous studies. In addition, an estimate for the Faraday rotation measure of the SGR and SNR is also carried out.

Local starburst conditions and formation of GRB 980425 / SN 1998bw within a collisional ring

We present the first spatially resolved study of molecular gas in the vicinity of a Gamma Ray Burst, using CO(2-1) emission line observations with the Atacama Large Millimetre Array (ALMA) at ~50 pc scales. The host galaxy of GRB 980425 contains a ring of high column density HI gas which is likely to have formed due to a collision between the GRB host and its companion galaxy, within which the GRB is located. We detect eleven molecular gas clumps in the galaxy, seven of which are within the gas ring. The clump closest to the GRB position is at a projected separation of ~280 pc. Although it is plausible that the GRB progenitor was ejected from clusters formed in this clump, we argue that the in situ formation of the GRB progenitor is the most likely scenario. We measure the molecular gas masses of the clumps and find them to be sufficient for forming massive star clusters. The molecular gas depletion times of the clumps show a variation of ~2 dex, comparable with the large variation in depletion times found in starburst galaxies in the nearby Universe. This demonstrates the presence of starburst modes of star formation on local scales in the galaxy, even while the galaxy as a whole cannot be categorised as a starburst based on its global properties. Our findings suggest that the progenitor of GRB 9802425 was originated in a young massive star cluster formed in the starburst mode of star formation.

Higgs boson emerging from the dark

We propose a new non-thermal mechanism of dark matter production based on vacuum misalignment. A global $X$-charge asymmetry is generated at high temperatures, under which both the will-be Higgs and the dark matter are charged. At lower energies, the vacuum changes alignment and breaks the $U(1)_X$, leading to the emergence of the Higgs and of a fraction of charge asymmetry stored in the stable dark matter relic. This mechanism can be present in a wide variety of models based on vacuum misalignment, and we demonstrate it in a composite Higgs template model, where all the necessary ingredients are naturally present. A light pseudo-scalar $\eta$ is always predicted, with interesting implications for cosmology, future supernova observations and exotic $Z \to \gamma \eta$ decays.

Electron acceleration in non-relativistic quasi-perpendicular collisionless shocks

Xu, Rui; Spitkovsky, Anatoly; Caprioli, Damiano

We study diffusive shock acceleration (DSA) of electrons in non-relativistic quasi-perpendicular shocks using self-consistent one-dimensional particle-in-cell (PIC) simulations. By exploring the parameter space of sonic and Alfv\'{e}nic Mach numbers we find that high Mach number quasi-perpendicular shocks can efficiently accelerate electrons to power-law downstream spectra with slopes consistent with DSA prediction. Electrons are reflected by magnetic mirroring at the shock and drive non-resonant waves in the upstream. Reflected electrons are trapped between the shock front and upstream waves and undergo multiple cycles of shock drift acceleration before the injection into DSA. Strong current-driven waves also temporarily change the shock obliquity and cause mild proton pre-acceleration even in quasi-perpendicular shocks, which otherwise do not accelerate protons. These results can be used to understand nonthermal emission in supernova remnants and intracluster medium in galaxy clusters.

The possible disappearance of a massive star in the low metallicity galaxy PHL 293B

We investigate a suspected very massive star in one of the most metal-poor dwarf galaxies, PHL~293B. Excitingly, we find the sudden disappearance of the stellar signatures from our 2019 spectra, in particular the broad H lines with P~Cygni profiles that have been associated with a massive luminous blue variable (LBV) star. Such features are absent from our spectra obtained in 2019 with the ESPRESSO and X-shooter instruments of the ESO's VLT. We compute radiative transfer models using CMFGEN that fit the observed spectrum of the LBV and are consistent with ground-based and archival Hubble Space Telescope photometry. Our models show that during 2001--2011 the LBV had a luminosity $L_* = 2.5-3.5 \times 10^6 ~L_{\odot}$, a mass-loss rate $\dot{M} = 0.005-0.020 ~M_{\odot}$~yr$^{-1}$, a wind velocity of 1000~km~s$^{-1}$, and effective and stellar temperatures of $T_\mathrm{eff} = 6000-6800$~K and $T_\mathrm{*}=9500-15000$~K. These stellar properties indicate an eruptive state. We consider two main hypotheses for the absence of the broad emission components from the spectra obtained since 2011. One possibility is that we are seeing the end of an LBV eruption of a surviving star, with a mild drop in luminosity, a shift to hotter effective temperatures, and some dust obscuration. Alternatively, the LBV could have collapsed to a massive black hole without the production of a bright supernova.

The Supernova Remnant Population of NGC6946 as Observed in [Fe II] 1.644 $\mu$m with HST

Long, Knox S.; Blair, William P.; Winkler, P. Frank; Lacey, Christina K.

NGC6946 is a high star formation rate face-on spiral galaxy that has hosted ten supernovae since 1917. Not surprisingly, a large number of supernova remnants and candidates have been identified either as optical nebulae with high [S II]:H$\alpha$ line ratios (147) or as compact non-thermal radio sources (35). However, there are only seven overlaps between these two samples. Here, we apply [Fe II] 1.644 $\mu$m emission as a new diagnostic to search for supernova remnants in an attempt to resolve this discrepancy. [Fe II] is expected to be relatively strong in the radiative shocks of supernova remnants and almost absent in HII regions. It is less susceptible to the effects of absorption along the line of sight than the optical lines normally used to identify remnants. Using data from the WFC3 camera on HST}, we identify 132 [Fe II] emission nebulae in NGC6946 as likely supernova remnants. Of these, 54 align with previously known optical supernova remnants. The remaining 78 objects are new; of these 44 are visible in new HST imagery in H$\alpha$ and [S II]. This brings the total number of supernova remnant candidates (from optical and/or IR data) in NGC6946 to 225. A total of 14 coincidences with radio supernova remnant candidates (out of 30 in our search area) are found in this expanded list. The identification of so many new remnant candidates validates the use of [Fe II] imagery for finding remnants, and suggests that previous remnant searches in other galaxies may be far from complete.

Interacting dark energy in the early 2020s: a promising solution to the $H_0$ and cosmic shear tensions

Di Valentino, Eleonora; Melchiorri, Alessandro; Mena, Olga; Vagnozzi, Sunny

We examine interactions between dark matter and dark energy in light of the latest cosmological observations, focusing on a specific model with coupling proportional to the dark energy density. Our data includes Cosmic Microwave Background (CMB) measurements from the Planck 2018 legacy data release, late-time measurements of the expansion history from Baryon Acoustic Oscillations (BAO) and Supernovae Type Ia (SNeIa), galaxy clustering and cosmic shear measurements from the Dark Energy Survey Year 1 results, and the 2019 local distance ladder measurement of the Hubble constant $H_0$ from the Hubble Space Telescope. Considering Planck data both in combination with BAO or SNeIa data reduces the $H_0$ tension to a level which could possibly be compatible with a statistical fluctuation. The very same model also significantly reduces the $\Omega_{\rm m}-\sigma_8$ tension between CMB and cosmic shear measurements. Interactions between the dark sectors of our Universe remain therefore a promising joint solution to these persisting cosmological tensions.

SN 2019ehk: A Double-Peaked Ca-rich Transient with Luminous X-ray Emission and Shock-Ionized Spectral Features

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

On the cosmological performance of photometric classified supernovae with machine learning

Santos, Marcelo Vargas dos; Quartin, Miguel; Reis, Ribamar R. R.

The efficient classification of different types of supernova is one of the most important problems for observational cosmology. However, spectroscopic confirmation of most objects in upcoming photometric surveys, such as the The Rubin Observatory Legacy Survey of Space and Time (LSST), will be unfeasible. The development of automated classification processes based on photometry has thus become crucial. In this paper we investigate the performance of machine learning (ML) classification on the final cosmological constraints using simulated lightcurves from The Supernova Photometric Classification Challenge, released in 2010. We study the use of different feature sets for the lightcurves and many different ML pipelines based on either decision tree ensembles or automated search processes. To construct the final catalogs we propose a threshold selection method, by employing a Bias-Variance tradeoff. This is a very robust and efficient way to minimize the Mean Squared Error. With this method we were able to get very strong cosmological constraints, which allowed us to keep ~75% of the total information in the type Ia SNe when using the SALT2 feature set and ~33% for the other cases (based on either the Newling model or on standard wavelet decomposition).

The Carnegie Supernova Project-I: Correlation Between Type Ia Supernovae and Their Host Galaxies from Optical to Near-Infrared Bands

We present optical and near-infrared ($ugriYJH$) photometry of host galaxies of Type Ia supernovae (SN~Ia) observed by the \textit{Carnegie Supernova Project-I}. We determine host galaxy stellar masses and, for the first time, study their correlation with SN~Ia standardized luminosity across optical and near-infrared ($uBgVriYJH$) bands. In the individual bands, we find that SNe~Ia are more luminous in more massive hosts with luminosity offsets ranging between $-0.07 \pm0.03$ mag to $-0.15\pm0.04$ mag after light-curve standardization. The slope of the SN~Ia Hubble residual-host mass relation is negative across all $uBgVriYJH$ bands with values ranging between $-0.036\pm 0.025$ mag/dex to $-0.097\pm 0.027$ mag/dex -- implying that SNe~Ia in more massive galaxies are brighter than expected. The near-constant observed correlations across optical and near-infrared bands indicate that dust may not play a significant role in the observed luminosity offset--host mass correlation. We measure projected separations between SNe~Ia and their host centers, and find that SNe~Ia that explode beyond a projected 10 kpc have a $\rm 30\% \ to \ 50\%$ reduction of the dispersion in Hubble residuals across all bands -- making them a more uniform subset of SNe~Ia. Dust in host galaxies, peculiar velocities of nearby SN~Ia, or a combination of both may drive this result as the color excesses of SNe~Ia beyond 10 kpc are found to be generally lower than those interior, but there is also a diminishing trend of the dispersion as we exclude nearby events. We do not find that SN~Ia average luminosity varies significantly when they are grouped in various host morphological types. Host galaxy data from this work will be useful, in conjunction with future high-redshift samples, in constraining cosmological parameters.

Significantly super-Chandrasekhar mass-limit of white dwarfs in noncommutative geometry

Kalita, Surajit; Mukhopadhyay, Banibrata; Govindarajan, T. R.

Chandrasekhar made the startling discovery about nine decades back that the mass of compact object white dwarf has a limiting value, once nuclear fusion reactions stop therein. This is the Chandrasekhar mass-limit which is $\sim1.4M_\odot$ for a non-rotating non-magnetized white dwarf. On approaching this limiting mass, a white dwarf is believed to spark off with an explosion called type Ia supernova, which is considered to be a standard candle. However, observations of several over-luminous, peculiar type Ia supernovae indicate that the Chandrasekhar mass-limit to be significantly larger. By considering noncommutativity among the components of position and momentum variables, hence uncertainty in their measurements, at the quantum scales, we show that the mass of white dwarfs could be significantly super-Chandrasekhar and thereby arrive at a new mass-limit $\sim 2.6M_\odot$, explaining a possible origin of over-luminous peculiar type Ia supernovae. The idea of noncommutativity, apart from the Heisenberg's uncertainty principle, is there for quite sometime, without any observational proof however. Our finding offers a plausible astrophysical evidence of noncommutativity, arguing for a possible second standard candle, which has many far-reaching implications.

Implications of a possible TeV break in the cosmic-ray electron and positron flux

Ding, Yu-Chen; Li, Nan; Wei, Chun-Cheng; Wu, Yue-Liang; Zhou, Yu-Feng

A TeV spectral break in the total flux of cosmic-ray electrons and positrons (CREs) at which the spectral power index softens from $\sim 3$ to $\sim 4$ has been observed by H.E.S.S. and recently confirmed by DAMPE with a high significance of $6.6~\sigma$. Such an observation is apparently inconsistent with the data from other experiments such as Fermi-LAT, AMS-02 and CALET. We perform a global analysis to the latest CRE data including Fermi-LAT, AMS-02, CALET, DAMPE and H.E.S.S. with energy scale uncertainties taken into account to improve the consistency between the data sets. The fit result strongly favors the existence of the break at $\sim 1$ TeV with an even higher statistical significance of $13.3~\sigma$. In view of the tentative CRE break, we revisit a number of models of nearby sources, such as a single generic Pulsar Wind Nebula (PWN), known multiple PWNe from the ATNF catalog, and their combinations with either an additional Dark Matter (DM) component or a Supernova Remnant (SNR). We show that the CRE break at $\sim 1$ TeV, together with the known CR positron excess points towards the possibility that the nearby sources should be highly charge asymmetric. Among the models under consideration, the one with a PWN plus SNR is most favored by the current data. The favoured distance and age of the PWN and SNR sources are both within $0.6$ kpc and around $10^{5}$ yr respectively. Possible candidate sources include PSR J0954-5430, Vela and Monogem ring, etc. We find that for the models under consideration, the additional DM component is either unnecessary, or predicts too much photons in tension with the H.E.S.S. data of $\gamma$-rays from the direction of the Galactic Center. We also show that the current measurement of the anisotropies in the arrival direction of the CRE can be useful in determining the property of the sources.

Cobalt and copper abundances in 56 Galactic bulge red giants

The Milky Way bulge is an important tracer of the early formation and chemical enrichment of the Galaxy. The abundances of different iron-peak elements in field bulge stars can give information on the nucleosynthesis processes that took place in the earliest supernovae. Cobalt (Z=27) and copper (Z=29) are particularly interesting.We aim to identify the nucleosynthesis processes responsible for the formation of the iron-peak elements Co and Cu. Methods. We derived abundances of the iron-peak elements cobalt and copper in 56 bulge giants, 13 of which were red clump stars. High-resolution spectra were obtained using FLAMES-UVES at the ESO Very Large Telescope by our group in 2000-2002, which appears to be the highest quality sample of high-resolution data on bulge red giants obtained in the literature to date. Over the years we have derived the abundances of C, N, O, Na, Al, Mg; the iron-group elements Mn and Zn; and neutron-capture elements. In the present work we derive abundances of the iron-peak elements cobalt and copper. We also compute chemodynamical evolution models to interpret the observed behaviour of these elements as a function of iron. The sample stars show mean values of [Co/Fe]~0.0 at all metallicities, and [Cu/Fe]~0.0 for [Fe/H]>-0.8 and decreasing towards lower metallicities with a behaviour of a secondary element. We conclude that [Co/Fe] varies in lockstep with [Fe/H], which indicates that it should be produced in the alpha-rich freezeout mechanism in massive stars. Instead [Cu/Fe] follows the behaviour of a secondary element towards lower metallicities, indicating its production in the weak s-process nucleosynthesis in He-burning and later stages. The chemodynamical models presented here confirm the behaviour of these two elements (i.e. [Co/Fe] vs. [Fe/H]~constant and [Cu/Fe] decreasing with decreasing metallicities).

Formation of mass-gap objects in highly asymmetric mergers

Safarzadeh, Mohammadtaher; Loeb, Abraham

The LIGO/Virgo Collaboration (LVC) recently reported the detection of GW190814, a merger of a $23^{+1.0}_{-0.9}~M_{\odot}$ primary black hole (BH), and a $2.6^{+0.08}_{-0.08}~M_{\odot}$ secondary. The secondary's mass falls into the mass-gap regime, which refers to the scarcity of compact objects in the mass range of 2-5 $M_{\odot}$. The first clue to the formation of the GW190814 lies in the fact that the primary is a very massive BH. We suggest that the secondary was born as a neutron star (NS) where a significant amount of the supernova ejecta mass from its formation remained bound to the binary due to the presence of the massive BH companion. The bound mass forms a circumbinary accretion disk, and its accretion onto the NS created a mass-gap object. In this scenario, LIGO/Virgo will only detect mass-gap objects in binary mergers with an extreme mass ratio. We also predict a correlation between the mass of the secondary and the mass of the primary in such asymmetric mergers. Our model can be tested with future data from the LVC's third-observing run.

Energy Conservation in the thin layer approximation: II. The asymmetric classic case for supernovae remnants

Zaninetti, Lorenzo

Here we extend the conservation of energy in the framework of the thin layer approximation to the asymmetrical case. Four types of interstellar medium are analysed, in which the density follows an inverse square profile, a power law profile, an exponential profile and a toroidal profile. An analytical solution for the radius as a function of time and the polar angle in spherical coordinates is derived in the case of the inverse square profile. The analytical and numerical results are applied to two supernova remnants: SN 1987A and SN 1006. The back reaction due to the radiative losses is evaluated in the case of the inverse square profile for the surrounding medium. Two models for the image formation are presented, which explain the triple ring visible in SN 1987A and the jet feature of SN 1006.

MANTRA: A Machine Learning reference lightcurve dataset for astronomical transient event recognition

We introduce MANTRA, an annotated dataset of 4869 transient and 71207 non-transient object lightcurves built from the Catalina Real Time Transient Survey. We provide public access to this dataset as a plain text file to facilitate standardized quantitative comparison of astronomical transient event recognition algorithms. Some of the classes included in the dataset are: supernovae, cataclysmic variables, active galactic nuclei, high proper motion stars, blazars and flares. As an example of the tasks that can be performed on the dataset we experiment with multiple data pre-processing methods, feature selection techniques and popular machine learning algorithms (Support Vector Machines, Random Forests and Neural Networks). We assess quantitative performance in two classification tasks: binary (transient/non-transient) and eight-class classification. The best performing algorithm in both tasks is the Random Forest Classifier. It achieves an F1-score of 96.25% in the binary classification and 52.79% in the eight-class classification. For the eight-class classification, non-transients ( 96.83% ) is the class with the highest F1-score, while the lowest corresponds to high-proper-motion stars ( 16.79% ); for supernovae it achieves a value of 54.57% , close to the average across classes. The next release of MANTRA includes images and benchmarks with deep learning models.

Explosive Nucleosynthesis in Sub-Chandrasekhar Mass White Dwarf Models for Type Ia Supernovae: Dependence on Model Parameters

Leung, Shing-Chi; Nomoto, Ken'ichi

Recent observations of Type Ia supernovae (SNe Ia) have shown diversified properties of the explosion strength, light curves and chemical composition. To investigate possible origins of such diversities in SNe Ia, we have presented multi-dimensional hydrodynamical study of explosions and associated nucleosynthesis in the near Chandrasekhar mass carbon-oxygen (CO) white dwarfs (WDs) for a wide range of parameters (Leung and Nomoto 2018 ApJ). In the present paper, we extend our wide parameter survey of models to the explosions of sub-Chandrasekhar mass CO WDs. We take the double detonation model for the explosion mechanism. The model parameters of the survey include the metallicity of $Z = 0 - 5~Z_\odot$, the CO WD mass of $M = 0.90 - 1.20~M_\odot$, and the He envelope mass of $M_{\rm He} = 0.05 - 0.20~M_\odot$. We also study how the initial He detonation configuration, such as spherical, bubble, and ring shapes, triggers the C detonation. For these parameters, we derive the minimum He envelope mass necessary to trigger the C detonation. We then examine how the explosion dynamics and associated nucleosynthesis depend on these parameters, and compare our results with the previous representative models. We compare our nucleosynthesis yields with the unusual abundance patterns of Fe-peak elements and isotopes observed in SNe Ia 2011fe, 2012cg and 2014J, as well as SN Ia remnant 3C 397 to provide constraints on their progenitors and environments. We provide the nucleosynthesis yields table of the sub-Chandrasekhar mass explosions, to discuss their roles in the galactic chemical evolution and archaeology.

Electron Capture Supernovae of Super-AGB Stars: Sensitivity on Input Physics

Leung, Shing-Chi; Nomoto, Ken'ichi; Suzuki, Tomoharu

Stars of $\sim$ 8 - 10 $M_{\odot}$ on the main-sequence form strongly electron-degenerate O+Ne+Mg core and become super-AGB stars. If such an O+Ne+Mg core grows to 1.38 $M_\odot$, electron captures on $^{20}$Ne$(e,\nu_e)^{20}$F$(e,\nu_e)^{20}$O take place and ignite O-Ne deflagration around the center. In this paper, we perform two-dimensional hydrodynamics simulations of the propagation of the O-Ne flame to see whether such a flame induces a collapse of the O+Ne+Mg core due to subsequent electron capture behind the flame or triggers a thermonuclear explosion. We present a series of models to explore how the outcome depends on model parameters for the central density in the range from $10^{9.80}$ to $10^{10.20}$ g cm$^{-3}$, flame structure of both centered and off-centered ignition kernels, special and general relativistic effects, turbulent flame speed formula and the treatments of laminar burning phase. We find that the O+Ne+Mg core obtained from stellar evolutionary models has a high tendency to collapse into a neutron star. We obtain the bifurcation between the electron-capture collapse and thermonuclear explosion. We discuss the implication in nucleosynthesis and the possible observational signals of this class of supernovae.

The role of galactic dynamics in shaping the physical properties of giant molecular clouds in Milky Way-like galaxies

We examine the role of the large-scale galactic-dynamical environment in setting the properties of giant molecular clouds in Milky Way-like galaxies. We perform three high-resolution simulations of Milky Way-like discs with the moving-mesh hydrodynamics code Arepo, yielding a statistical sample of $\sim 80,000$ giant molecular clouds and $\sim 55,000$ HI clouds. We account for the self-gravity of the gas, momentum and thermal energy injection from supernovae and HII regions, mass injection from stellar winds, and the non-equilibrium chemistry of hydrogen, carbon and oxygen. By varying the external gravitational potential, we probe galactic-dynamical environments spanning an order of magnitude in the orbital angular velocity, gravitational stability, mid-plane pressure and the gradient of the galactic rotation curve. The simulated molecular clouds are highly overdense ($\sim 100\times$) and over-pressured ($\sim 25\times$) relative to the ambient interstellar medium. Their gravo-turbulent and star-forming properties are decoupled from the dynamics of the galactic mid-plane, so that the kpc-scale star formation rate surface density is related only to the number of molecular clouds per unit area of the galactic mid-plane. Despite this, the clouds display clear, statistically-significant correlations of their rotational properties with the rates of galactic shearing and gravitational free-fall. We find that galactic rotation and gravitational instability can influence their elongation, angular momenta, and tangential velocity dispersions. The lower pressures and densities of the HI clouds allow for a greater range of significant dynamical correlations, mirroring the rotational properties of the molecular clouds, while also displaying a coupling of their gravitational and turbulent properties to the galactic-dynamical environment.

ASASSN-18am/SN 2018gk : An overluminous Type IIb supernova from a massive progenitor

ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of $M_V \approx -20$ mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am declines extremely rapidly for a Type II SN, with a photospheric-phase decline rate of $\sim6.0$ mag (100d)$^{-1}$. Owing to the weakening of HI and the appearance of HeI in its later phases, ASASSN-18am is spectroscopically a Type IIb SN with a partially stripped envelope. However, its photometric and spectroscopic evolution show significant differences from typical SNe IIb. Using a radiative diffusion model, we find that the light curve requires a high synthesised $^{56}$Ni mass of $M_{\rm Ni} \sim0.4\rm M_\odot$ and ejecta with high kinetic energy $E_{\rm kin}= (5-9)\times10^{51} $ erg. Introducing a magnetar central engine still requires $M_{\rm Ni} \sim0.3\, M_\odot$ and $E_{\rm kin}= 3\times10^{51} $ erg. The high $^{56}$Ni mass is consistent with strong iron-group nebular lines in its spectra, which are also similar to several SNe Ic-BL with high $^{56}$Ni yields. The earliest spectrum shows "flash ionization" features, from which we estimate a mass-loss rate of $ \dot{M}\approx 2\times10^{-4} \, \rm M_\odot~yr^{-1} $. This wind density is too low to power the luminous light curve by ejecta-CSM interaction. We measure expansion velocities as high as $ 17,000 \rm~km~s^{-1}$ for H$\alpha$, which is remarkably high compared to other SNe II. We estimate an oxygen core mass of 1.7-3.1 $\rm M_\odot$ using the [OI] luminosity measured from a nebular-phase spectrum, implying a progenitor with a zero-age main sequence mass of 19-24 $\rm M_\odot$.

A DESGW Search for the Electromagnetic Counterpart to the LIGO/Virgo Gravitational Wave Binary Neutron Star Merger Candidate S190510g

We present the results from a search for the electromagnetic counterpart of the LIGO/Virgo event S190510g using the Dark Energy Camera (DECam). S190510g is a binary neutron star (BNS) merger candidate of moderate significance detected at a distance of 227$\pm$92 Mpc and localized within an area of 31 (1166) square degrees at 50\% (90\%) confidence. While this event was later classified as likely non-astrophysical in nature within 30 hours of the event, our short latency search and discovery pipeline identified 11 counterpart candidates, all of which appear consistent with supernovae following offline analysis and spectroscopy by other instruments. Later reprocessing of the images enabled the recovery of 6 more candidates. Additionally, we implement our candidate selection procedure on simulated kilonovae and supernovae under DECam observing conditions (e.g., seeing, exposure time) with the intent of quantifying our search efficiency and making informed decisions on observing strategy for future similar events. This is the first BNS counterpart search to employ a comprehensive simulation-based efficiency study. We find that using the current follow-up strategy, there would need to be 19 events similar to S190510g for us to have a 99\% chance of detecting an optical counterpart, assuming a GW170817-like kilonova. We further conclude that optimization of observing plans, which should include preference for deeper images over multiple color information, could result in up to a factor of 1.5 reduction in the total number of followups needed for discovery.

Supernova Muons: New Constraints on Z' Bosons, Axions, and ALPs

Croon, Djuna; Elor, Gilly; Leane, Rebecca K.; McDermott, Samuel D.

New light particles produced in supernovae can lead to additional energy loss and a consequent deficit in neutrino production in conflict with the neutrinos observed from Supernova 1987A (SN1987A). Contrary to the majority of previous SN1987A studies, we examine the impact of $Z'$ bosons, axions, and axion-like particles (ALPs) interacting with the muons produced in SN1987A. For the first time, we find constraints on generic $Z'$ bosons coupled to muons, and apply our results to particle models including gauged $L_\mu-L_\tau$ number, $U(1)_{L_\mu-L_\tau}$, and gauged $B-L$ number, $U(1)_{B-L}$. We constrain $Z'$ bosons with masses up to about 250-500 MeV, and down to about $10^{-9}$ in $Z'$-muon coupling. We also extend previous work on axion-muon couplings by examining the importance of loop-level interactions, as well as performing calculations over a wider range of axion masses. We constrain muon-coupled axions from arbitrarily low masses up to about 200-500 MeV, with bounds extending down to axion-muon couplings of approximately $10^{-8}$ GeV$^{-1}$. We conclude that supernovae broadly provide a sensitive probe of new lightly-coupled particles interacting with muons.

Soundness of Dark Energy properties

Di Valentino, Eleonora; Gariazzo, Stefano; Mena, Olga; Vagnozzi, Sunny

Type Ia Supernovae (SNeIa) used as standardizable candles have been instrumental in the discovery of cosmic acceleration, usually attributed to some form of dark energy (DE). Recent studies have raised the issue of whether intrinsic SNeIa luminosities might evolve with redshift. While the evidence for cosmic acceleration is robust to this possible systematic, the question remains of how much the latter can affect the inferred properties of the DE component responsible for cosmic acceleration. This is the question we address in this work. We use SNeIa distance moduli measurements from the Pantheon and JLA samples. We consider models where the DE equation of state is a free parameter, either constant or time-varying, as well as models where DE and dark matter interact, and finally a model-agnostic parametrization of effects due to modified gravity (MG). When SNeIa data are combined with Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements, we find strong degeneracies between parameters governing the SNeIa systematics, the DE parameters, and the Hubble constant $H_0$. These degeneracies significantly broaden the DE parameter uncertainties, in some cases leading to ${\cal O}(\sigma)$ shifts in the central values. However, including low-redshift Baryon Acoustic Oscillation and Cosmic Chronometer measurements, as well as CMB lensing measurements, considerably improves the previous constraints, and the only remaining effect of the examined systematic is a $\lesssim 40\%$ broadening of the uncertainties on the DE parameters. The constraints we derive on the MG parameters are instead basically unaffected by the systematic in question. We therefore confirm the overall soundness of dark energy properties.

A systematic study on the rise time-peak luminosity relation for bright optical transients powered by wind shock breakout

Suzuki, Akihiro; Moriya, Takashi J.; Takiwaki, Tomoya

A number of astrophyical transients originating from stellar explosions are powered by the collision of the ejected material with the circumstellar medium, which efficiently produces thermal radiation via shock dissipation. We investigate how such interaction-powered transients are distributed in the peak bolometric luminosity vs the rise time phase space. Taking the advantage of less time-consuming one-dimensional simulations with spherical symmetry, we calculated more than 500 models with different circumstellar mass and radius, ejecta mass and energy, and chemical compositions. The peak bolometric luminosity, the total radiated energy, and the rise time of the interaction-powered emission are measured for each simulated light curve. We consider how these characteristic quantities are determined as a function of the model parameters and discuss possible implications for the observed populations of (potential) interaction-powered transients, such as type IIn supernovae and fast blue optical transients.

The Proper Motion of the Central Compact Object RX J0822-4300 in the Supernova Remnant Puppis A, Revisited

Mayer, Martin; Becker, Werner; Patnaude, Daniel; Winkler, P. Frank; Kraft, Ralph

We present an improved proper motion measurement of the central compact object RX J0822-4300, located in the supernova remnant Puppis A. By employing a new data set taken in February 2019 by the High Resolution Camera aboard the Chandra X-ray Observatory, we approximately double the available temporal baseline for our analysis to slightly more than 19 years (7000 days). We correct for the astrometric inaccuracy of Chandra using calibrator stars with known optical positions that are detected in all observations. Thereby, we obtain absolute positions of RX J0822-4300 accurate to around $0.1^{\prime\prime}$ and from these a new best estimate for its total proper motion of $\mu_{\rm tot}= (80.4 \pm 7.7)\,\rm{mas\,yr}^{-1}$. For a remnant distance of 2 kpc, this corresponds to a projected kick velocity of $(763 \pm 73)\, \rm{km\,s}^{-1}$ at a position angle of $\phi_0 = (247.8 \pm 4.4)^{\circ}$. The proper motion measurement of RX J0822-4300 is used for discussing the kinematic age of Puppis A.

The Impact of Type Ia Supernovae in Quiescent Galaxies: I. Formation of the Multiphase Interstellar medium

Li, Miao; Li, Yuan; Bryan, Greg L.; Ostriker, Eve C.; Quataert, Eliot

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

The Impact of Type Ia Supernovae in Quiescent Galaxies: II. Energetics and Turbulence

Li, Miao; Li, Yuan; Bryan, Greg L.; Ostriker, Eve C.; Quataert, Eliot

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.