Showing 1–50 of 91 results for author: O'Connor, E

Unveiling the Nature of Gravitational-Wave Emission in Core-collapse Supernovae with Perturbative Analysis

Authors: Shuai Zha, Oliver Eggenberger Andersen, Evan P. O'Connor

Abstract: Gravitational waves (GWs) can provide crucial information about the central engines of core-collapse supernovae (CCSNe). In order to unveil the nature of GW emission in CCSNe, we apply perturbative analyses with the same underlying equations as simulations to diagnose oscillations of the proto-neutron star (PNS) during $\sim$1 s postbounce. In the pseudo-Newtonian case, we find that radial profile… ▽ More Gravitational waves (GWs) can provide crucial information about the central engines of core-collapse supernovae (CCSNe). In order to unveil the nature of GW emission in CCSNe, we apply perturbative analyses with the same underlying equations as simulations to diagnose oscillations of the proto-neutron star (PNS) during $\sim$1 s postbounce. In the pseudo-Newtonian case, we find that radial profiles of GW emission match well between the perturbative analysis with $l=2$ and simulations inside the PNS at \emph{any} frequency and time. This confirms that the GW emission of CCSNe arises from the global PNS oscillations in the perturbative regime. Based on this, we solve for the discrete eigenmodes with a free PNS surface and tentatively identify a set of $g$ modes and the $f$ mode contributing to the peak GW emission. We also offer a possible explanation for the power gap in the GW spectrum found in simulations that lies at the frequency with vanishing cumulative emission of the PNS. Our results enhance the predictive power of perturbative analyses in the GW signals of CCSNe. △ Less

Submitted 5 April, 2024; v1 submitted 4 March, 2024; originally announced March 2024.

Comments: 12 pages, 9 figures, accepted by PRD

Gray Two-moment Neutrino Transport: Comprehensive Tests and Improvements for Supernova Simulations

Authors: Haakon Andresen, Evan P. O'Connor, Oliver Eggenberger Andersen, Sean M. Couch

Abstract: In this work we extended an energy-integrated neutrino transport method to facilitate efficient, yet precise, modeling of compact astrophysical objects. We particularly focus on core-collapse supernovae. We implemented a gray neutrino-transport framework from the literature into FLASH and performed a detailed evaluation of its accuracy in core-collapse supernova simulations. Based on comparisons w… ▽ More In this work we extended an energy-integrated neutrino transport method to facilitate efficient, yet precise, modeling of compact astrophysical objects. We particularly focus on core-collapse supernovae. We implemented a gray neutrino-transport framework from the literature into FLASH and performed a detailed evaluation of its accuracy in core-collapse supernova simulations. Based on comparisons with results from simulations using energy-dependent neutrino transport, we incorporated several improvements to the original scheme. Our analysis shows that our gray neutrino transport method successfully reproduces key aspects from more complex energy-dependent transport across a variety of progenitors and equations of state. We find both qualitative and reasonable quantitative agreement with multi-group M1 transport simulations. However, the gray scheme tends to slightly favor shock revival. In terms of gravitational wave and neutrino signals, there is a good alignment with the energy-dependent transport, although we find 15-30% discrepancies in the average energy and luminosity of heavy-lepton neutrinos. Simulations using the gray transport are around four times faster than those using energy-dependent transport. △ Less

Submitted 1 July, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: 18 pages, 10 figures, published in A&A

Journal ref: A&A 687, A55 (2024)

On the pollution of white dwarfs by exo-Oort cloud comets

Authors: Christopher E. O'Connor, Dong Lai, Darryl Z. Seligman

Abstract: A large fraction of white dwarfs (WDs) have metal-polluted atmospheres, which are produced by accreting material from remnant planetary systems. The composition of the accreted debris broadly resembles that of rocky Solar System objects. Volatile-enriched debris with compositions similar to long-period comets (LPCs) is rarely observed. We attempt to reconcile this dearth of volatiles with the prem… ▽ More A large fraction of white dwarfs (WDs) have metal-polluted atmospheres, which are produced by accreting material from remnant planetary systems. The composition of the accreted debris broadly resembles that of rocky Solar System objects. Volatile-enriched debris with compositions similar to long-period comets (LPCs) is rarely observed. We attempt to reconcile this dearth of volatiles with the premise that exo-Oort clouds (XOCs) occur around a large fraction of planet-hosting stars. We estimate the comet accretion rate from an XOC analytically, adapting the 'loss cone' theory of LPC delivery in the Solar System. We investigate the dynamical evolution of an XOC during late stellar evolution. Using numerical simulations, we show that 1 to 30 per cent of XOC objects remain bound after anisotropic stellar mass loss imparting a WD natal kick of $\sim$1 km/s. We also characterize the surviving comets' distribution function. Surviving planets orbiting a WD can prevent the accretion of XOC comets by the star. A planet's 'dynamical barrier' is effective at preventing comet accretion if the energy kick imparted by the planet exceeds the comet's orbital binding energy. By modifying the loss cone theory, we calculate the amount by which a planet reduces the WD's accretion rate. We suggest that the scarcity of volatile-enriched debris in polluted WDs is caused by an unseen population of 10-100 AU scale giant planets acting as barriers to incoming LPCs. Finally, we constrain the amount of volatiles delivered to a planet in the habitable zone of an old, cool WD. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: 18 pages, 12 figures; submitted to MNRAS. Comments welcome!

Giant planet engulfment by evolved giant stars: light curves, asteroseismology, and survivability

Authors: Christopher E. O'Connor, Lars Bildsten, Matteo Cantiello, Dong Lai

Abstract: About ten percent of Sun-like ($1$-$2 M_\odot$) stars will engulf a $1$-$10 M_{\rm J}$ planet as they expand during the red giant branch (RGB) or asymptotic giant branch (AGB) phase of their evolution. Once engulfed, these planets experience a strong drag force in the star's convective envelope and spiral inward, depositing energy and angular momentum. For these mass ratios, the inspiral takes… ▽ More About ten percent of Sun-like ($1$-$2 M_\odot$) stars will engulf a $1$-$10 M_{\rm J}$ planet as they expand during the red giant branch (RGB) or asymptotic giant branch (AGB) phase of their evolution. Once engulfed, these planets experience a strong drag force in the star's convective envelope and spiral inward, depositing energy and angular momentum. For these mass ratios, the inspiral takes $\sim 10$-$10^{2}$ years ($\sim 10^{2}$-$10^{3}$ orbits); the planet undergoes tidal disruption at a radius of $\sim R_\odot$. We use the Modules for Experiments in Stellar Astrophysics (MESA) software instrument to track the stellar response to the energy deposition while simultaneously evolving the planetary orbit. For RGB stars, as well as AGB stars with $M_{\rm p} \lesssim 5 M_{\rm J}$ planets, the star responds quasistatically but still brightens measurably on a timescale of years. In addition, asteroseismic indicators, such as the frequency spacing or rotational splitting, differ before and after engulfment. For AGB stars, engulfment of a $M_{\rm p} \gtrsim 5 M_{\rm J}$ planet drives supersonic expansion of the envelope, causing a bright, red, dusty eruption similar to a "luminous red nova." Based on the peak luminosity, color, duration, and expected rate of these events, we suggest that engulfment events on the AGB could be a significant fraction of low-luminosity red novae in the Galaxy. We do not find conditions where the envelope is ejected prior to the planet's tidal disruption, complicating the interpretation of short-period giant planets orbiting white dwarfs as survivors of common-envelope evolution. △ Less

Submitted 13 June, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: 24 pages, 11 figures, 1 table. Accepted to ApJ

The Force Explosion Condition is Consistent with Spherically Symmetric CCSN Explosions

Authors: Mariam Gogilashvili, Jeremiah W. Murphy, Evan P. O'Connor

Abstract: One of the major challenges in Core-collapse Supernova (CCSN) theory is to predict which stars explode and which collapse to black holes. Gogilashvili and Murphy (2022) derived an analytic force explosion condition (FEC) and showed that the FEC is consistent with CCSN simulations that use the light-bulb approximation for neutrino heating and cooling. In this follow-up manuscript, we show that the… ▽ More One of the major challenges in Core-collapse Supernova (CCSN) theory is to predict which stars explode and which collapse to black holes. Gogilashvili and Murphy (2022) derived an analytic force explosion condition (FEC) and showed that the FEC is consistent with CCSN simulations that use the light-bulb approximation for neutrino heating and cooling. In this follow-up manuscript, we show that the FEC is consistent with the explosion condition when using actual neutrino transport in GR1D simulations (O'CONNOR 2015). Since most 1D simulations do not explode, to facilitate this test, we enhance the heating efficiency within the gain region. To compare the analytic FEC and radiation-hydrodynamic simulations, this manuscript also presents a practical translation of the physical parameters. For example: we replace the neutrino power deposited in the gain region, $L_ντ_g$, with the net neutrino heating in the gain region; rather than assuming that $\dot{M}$ is the same everywhere, we calculate $\dot{M}$ within the gain region; and we use the neutrino opacity at the gain radius. With small, yet practical modifications, we show that the FEC predicts the explosion conditions in spherically symmetric CCSN simulations that use neutrino transport. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: Figure 5 most clearly highlights the fidelity and the efficacy of the Force Explosion Condition

Inferring Type II-P Supernova Progenitor Masses from Plateau Luminosities

Authors: Brandon L. Barker, Evan P. O'Connor, Sean M. Couch

Abstract: Connecting observations of core-collapse supernova explosions to the properties of their massive star progenitors is a long-sought, and challenging, goal of supernova science. Recently, Barker et al. (2022) presented bolometric light curves for a landscape of progenitors from spherically symmetric neutrino-driven core-collapse supernova (CCSN) simulations using an effective model. They find a tigh… ▽ More Connecting observations of core-collapse supernova explosions to the properties of their massive star progenitors is a long-sought, and challenging, goal of supernova science. Recently, Barker et al. (2022) presented bolometric light curves for a landscape of progenitors from spherically symmetric neutrino-driven core-collapse supernova (CCSN) simulations using an effective model. They find a tight relationship between the plateau luminosity of the Type II-P CCSN light curve and the terminal iron core mass of the progenitor. Remarkably, this allows us to constrain progenitor properties with photometry alone. We analyze a large observational sample of Type II-P CCSN light curves and estimate a distribution of iron core masses using the relationship of Barker et al 2022. The inferred distribution matches extremely well with the distribution of iron core masses from stellar evolutionary models, and namely, contains high-mass iron cores that suggest contributions from very massive progenitors in the observational data. We use this distribution of iron core masses to infer minimum and maximum mass of progenitors in the observational data. Using Bayesian inference methods to locate optimal initial mass function parameters, we find M$_{\mathrm{min}}=9.8^{+0.37}_{-0.27}$ and M$_{\mathrm{max}}=24.0^{+3.9}_{-1.9}$ solar masses for the observational data. △ Less

Submitted 29 January, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: Accepted to ApJ Letters. 9 pages, 4 figures

Report number: LA-UR-22-31879

A Parameterized Neutrino Emission Model to Study Mass Ejection in Failed Core-collapse Supernovae

Authors: A. S. Schneider, E. O'Connor

Abstract: Some massive stars end their lives as \textit{failed} core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to n… ▽ More Some massive stars end their lives as \textit{failed} core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within $\sim100\,{\rm ms}$ and up to $\sim10^6\,{\rm s}$. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 21 pages, 12 figures, 2 tables

Impact of Rotation on the Multimessenger Signatures of a Hadron-quark Phase Transition in Core-collapse Supernovae

Authors: Shuai Zha, Evan O'Connor

Abstract: We study the impact of rotation on the multimessenger signals of core-collapse supernovae (CCSNe) with the occurrence of a first-order hadron-quark phase transition (HQPT). We simulate CCSNe with the \texttt{FLASH} code starting from a 20~$M_\odot$ progenitor with different rotation rates, and using the RDF equation of state from \textit{Bastian} 2021 that prescribes the HQPT. Rotation is found to… ▽ More We study the impact of rotation on the multimessenger signals of core-collapse supernovae (CCSNe) with the occurrence of a first-order hadron-quark phase transition (HQPT). We simulate CCSNe with the \texttt{FLASH} code starting from a 20~$M_\odot$ progenitor with different rotation rates, and using the RDF equation of state from \textit{Bastian} 2021 that prescribes the HQPT. Rotation is found to delay the onset of the HQPT and the resulting dynamical collapse of the protocompact star (PCS) due to the centrifugal support. All models with the HQPT experience a second bounce shock which leads to a successful explosion. The oblate PCS as deformed by rotation gives rise to strong gravitational-wave (GW) emission around the second bounce with a peak amplitude larger by a factor of $\sim10$ than that around the first bounce. The breakout of the second bounce shock at the neutrinosphere produces a $\barν_e$-rich neutrino burst with a luminosity of serveral 10$^{53}$~erg~s$^{-1}$. In rapidly rotating models the PCS pulsation following the second bounce generates oscillations in the neutrino signal after the burst. In the fastest rotating model with the HQPT, a clear correlation is found between the oscillations in the GW and neutrino signals immediately after the second bounce. In addition, the HQPT-induced collapse leads to a jump in the ratio of rotational kinetic energy to gravitational energy ($β$) of the PCS, for which persistent GW emission may arise due to secular nonaxisymmetric instabilities. △ Less

Submitted 5 January, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: 15 pages, 15 figures, multimessenger data available at https://doi.org/10.5281/zenodo.6965988

Journal ref: Physical Review D 106, 123037 (2022)

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

Neutrino Driven Explosions aided by Axion Cooling in Multidimensional Simulations of Core-Collapse Supernovae

Authors: Aurore Betranhandy, Evan O'Connor

Abstract: In this study, we present the first multidimensional core-collapse supernovae (CCSNe) simulations including QCD axions in order to assess the impact on the CCSN explosion mechanism. We include axions in our simulations through the nucleon-nucleon bremsstrahlung emission channel and as a pure energy-sink term under the assumption that the axions free-stream after being emitted. We perform both sphe… ▽ More In this study, we present the first multidimensional core-collapse supernovae (CCSNe) simulations including QCD axions in order to assess the impact on the CCSN explosion mechanism. We include axions in our simulations through the nucleon-nucleon bremsstrahlung emission channel and as a pure energy-sink term under the assumption that the axions free-stream after being emitted. We perform both spherically symmetric (1D) and axisymmetric (2D) simulations. In 1D, we utilize a parameterized heating scheme to achieve explosions, whereas in 2D we self-consistently realize explosions through the neutrino heating mechanism. Our 2D results for a $20 M_\odot$ progenitor show an impact of the axion emission on the shock behavior and the explosion time when considering values of the Peccei-Quinn energy scale $f_a \leq 2 \times 10^8$ GeV. The strong cooling due to the axion emission accelerates the contraction of the core and leads to more efficient neutrino heating and earlier explosions. For the axion emission formalism utilized, the values of $f_a$ that impact the explosion are close to, but in tension with current limits based on the neutrinos detected from SN1987A. However, given the non-linear behavior of the emission and the multidimensional nature of CCSNe, we suggest that a self-consistent, multidimensional approach to simulating CCSNe, including any late time accretion and cooling, is needed to fully explore the axion bounds from supernovae and the impact on the CCSN explosion mechanism. △ Less

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

Comments: 18 pages, 12 figures

Journal ref: Phys. Rev. D 106, Vol. 106, Iss. 6, 063019, 2022

Detectability of neutrino-signal fluctuations induced by the hadron-quark phase transition in failing core-collapse supernovae

Authors: Zidu Lin, Shuai Zha, Evan P. O'Connor, Andrew W. Steiner

Abstract: We introduce a systematic and quantitative methodology for establishing the presence of neutrino oscillatory signals due to the hadron-quark phase transition (PT) in failing core-collapse supernovae from the observed neutrino event rate in water- or ice-based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test-statistic; it is employed for quantitative ana… ▽ More We introduce a systematic and quantitative methodology for establishing the presence of neutrino oscillatory signals due to the hadron-quark phase transition (PT) in failing core-collapse supernovae from the observed neutrino event rate in water- or ice-based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test-statistic; it is employed for quantitative analysis of neutrino signals without assuming the frequency, amplitude, starting time, and duration of the PT-induced oscillations present in the neutrino events and thus it is suitable for analyzing neutrino signals from a wide variety of numerical simulations. We test the validity of this method by using a core-collapse simulation of a 17 solar-mass star by Zha \emph{et al.} (2021). Based on this model, we further report the presence of a PT-induced oscillations quantitatively for a core-collapse supernovae out to a distance of $\sim 10$ kpc, $\sim 5$ kpc for IceCube and to a distance of $\sim 10$ kpc, $\sim 5$ kpc and $\sim 1$ kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova and the study of hadron-quark phase in the core of core-collapse supernovae. △ Less

Submitted 26 December, 2023; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: We added an appendix discussing the model-dependence of the numerical simulation and its neutrino signals. Accepted by PRD

Comparison of Electron Capture Rates in the N=50 Region using 1D Simulations of Core-collapse Supernovae

Authors: Zac Johnston, Sheldon Wasik, Rachel Titus, MacKenzie L. Warren, Evan P. O'Connor, Remco Zegers, Sean M. Couch

Abstract: Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N=50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implemen… ▽ More Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N=50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parameterized N=50 rates (LMP), and an older independent-particle approximation (IPA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts on the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and IPA, and yet slightly more favorable to explosion than both. △ Less

Submitted 11 November, 2022; v1 submitted 18 February, 2022; originally announced February 2022.

Comments: 14 pages, 8 figures, 1 table. Updated Fig 5, minor corrections, renamed SNA to IPA. Accepted by ApJ

Journal ref: ApJ 939 15 (2022)

Hydrodynamic Simulations of Electron-capture Supernovae: Progenitor and Dimension Dependence

Authors: Shuai Zha, Evan P. O'Connor, Sean M. Couch, Shing-Chi Leung, Ken'ichi Nomoto

Abstract: We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in \texttt{FLASH} with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with 6 progenitor models that, at the time of collapse, span a range of propagati… ▽ More We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in \texttt{FLASH} with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with 6 progenitor models that, at the time of collapse, span a range of propagating flame front radii. For comparison, we also perform a simulation with the same setup using the canonical, spherically-symmetrical progenitor model n8.8. We found that the variations in the progenitor models inherited from simulations of stellar evolution and flame propagation do not significantly alter the global properties of the neutrino-driven ECSN explosion, such as the explosion energy ($\sim1.36$-$1.48\times10^{50}$ erg) and the mass ($\sim0.017$-$0.018M_\odot$) and composition of the ejecta. Due to aspherical perturbations induced by the 2D flame, the ejecta contains a small amount ($\lesssim1.8\times10^{-3}~M_\odot$) of low-$Y_e$ ($0.35<Y_e<0.4$) component. The baryonic mass of the protoneutron star is $\sim1.34~M_\odot$ ($\sim1.357~M_\odot$) with the new (n8.8) progenitor models when simulations end at $\sim400$ ms and the discrepancy is due to updated weak-interaction rates in the progenitor evolutionary simulations. Our results reflect the nature of ECSN progenitors containing a strongly degenerate ONeMg core and suggest a standardized ECSN explosion initialized by ONeMg core collapse. Moreover, we carry out a rudimentary three-dimensional simulation and find that the explosion properties are fairly compatible with the 2D counterpart. Our paper facilitates a more thorough understanding of ECSN explosions following the ONeMg core collapse, though more three-dimensional simulations are still needed. △ Less

Submitted 13 June, 2022; v1 submitted 30 December, 2021; originally announced December 2021.

Comments: 12 pages, 13 figures, progenitor models available at zenodo, https://doi.org/10.5281/zenodo.5748457, Accepted by MNRAS

Secular chaos in white-dwarf planetary systems: origins of metal pollution and short-period planetary companions

Authors: Christopher E. O'Connor, Jean Teyssandier, Dong Lai

Abstract: Secular oscillations in multi-planet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This "secular chaos" readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging… ▽ More Secular oscillations in multi-planet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This "secular chaos" readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging method, with emphasis on the relationship between the planets' properties and the time-scale and efficiency of chaotic diffusion. We find that secular chaos can excite extreme eccentricities on time-scales spanning several orders of magnitude in a given system. We apply our results to the evolution of planetary systems around white dwarfs (WDs), specifically the tidal disruption and high-eccentricity migration of planetesimals and planets. We find that secular chaos in a planetesimal belt driven by large ($\gtrsim 10 M_{\oplus}$), distant ($\gtrsim 10 \, {\rm au}$) planets can sustain metal accretion onto a WD over Gyr time-scales. We constrain the total mass of planetesimals initially present within the chaotic zone by requiring that the predicted mass delivery rate to the Roche limit be consistent with the observed metal accretion rates of WDs with atmospheric pollution throughout the cooling sequence. Based on the occurrence of long-period exoplanets and exo-asteroid belts, we conclude that secular chaos can be a significant (perhaps dominant) channel for polluting solitary WDs. Secular chaos can also produce short-period planets and planetesimals around WDs in concert with various circularization mechanisms. We discuss prospects for detecting exoplanets driving secular chaos around WDs using direct imaging and microlensing. △ Less

Submitted 26 April, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 21 pages, 10 figures, 2 tables; accepted in MNRAS

Gemini/GMOS Transmission Spectroscopy of the Grazing Planet Candidate WD 1856+534 b

Authors: Siyi Xu, Hannah Diamond-Lowe, Ryan J. MacDonald, Andrew Vanderburg, Simon Blouin, P. Dufour, Peter Gao, Laura Kreidberg, S. K. Leggett, Andrew W. Mann, Caroline V. Morley, Andrew W. Stephens, Christopher E. O'Connor, Pa Chia Thao, Nikole K. Lewis

Abstract: WD 1856+534 b is a Jupiter-sized, cool giant planet candidate transiting the white dwarf WD 1856+534. Here, we report an optical transmission spectrum of WD 1856+534 b obtained from ten transits using the Gemini Multi-Object Spectrograph. This system is challenging to observe due to the faintness of the host star and the short transit duration. Nevertheless, our phase-folded white light curve reac… ▽ More WD 1856+534 b is a Jupiter-sized, cool giant planet candidate transiting the white dwarf WD 1856+534. Here, we report an optical transmission spectrum of WD 1856+534 b obtained from ten transits using the Gemini Multi-Object Spectrograph. This system is challenging to observe due to the faintness of the host star and the short transit duration. Nevertheless, our phase-folded white light curve reached a precision of 0.12 %. WD 1856+534 b provides a unique transit configuration compared to other known exoplanets: the planet is $8\times$ larger than its star and occults over half of the stellar disc during mid-transit. Consequently, many standard modeling assumptions do not hold. We introduce the concept of a `limb darkening corrected, time-averaged transmission spectrum' and propose that this is more suitable than $(R_{\mathrm{p}, λ} / R_{\mathrm{s}})^2$ for comparisons to atmospheric models for planets with grazing transits. We also present a modified radiative transfer prescription. Though the transmission spectrum shows no prominent absorption features, it is sufficiently precise to constrain the mass of WD 1856+534 b to be > 0.84 M$_\mathrm{J}$ (to $2 \, σ$ confidence), assuming a clear atmosphere and a Jovian composition. High-altitude cloud decks can allow lower masses. WD 1856+534 b could have formed either as a result of common envelope evolution or migration under the Kozai-Lidov mechanism. Further studies of WD 1856+534 b, alongside new dedicated searches for substellar objects around white dwarfs, will shed further light on the mysteries of post-main sequence planetary systems. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 21 pages, 10 figures, accepted for publication in AJ

EuCAPT White Paper: Opportunities and Challenges for Theoretical Astroparticle Physics in the Next Decade

Authors: R. Alves Batista, M. A. Amin, G. Barenboim, N. Bartolo, D. Baumann, A. Bauswein, E. Bellini, D. Benisty, G. Bertone, P. Blasi, C. G. Böhmer, Ž. Bošnjak, T. Bringmann, C. Burrage, M. Bustamante, J. Calderón Bustillo, C. T. Byrnes, F. Calore, R. Catena, D. G. Cerdeño, S. S. Cerri, M. Chianese, K. Clough, A. Cole, P. Coloma , et al. (112 additional authors not shown)

Abstract: Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, und… ▽ More Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: White paper of the European Consortium for Astroparticle Theory (EuCAPT). 135 authors, 400 endorsers, 133 pages, 1382 references

Effect of the Nuclear Equation of State on Relativistic-Turbulence Induced Core-Collapse Supernovae

Authors: Luca Boccioli, Grant J. Mathews, In-Saeng Suh, Evan P. O'Connor

Abstract: The nuclear equation of state is an important component in the evolution of core collapse supernovae. In this paper we make a survey of various equations of state in the literature and analyze their effect on spherical core-collapse models in which the effects of three-dimensional turbulence is modeled by a general relativistic formulation of Supernova Turbulence in Reduced dimensionality (STIR).… ▽ More The nuclear equation of state is an important component in the evolution of core collapse supernovae. In this paper we make a survey of various equations of state in the literature and analyze their effect on spherical core-collapse models in which the effects of three-dimensional turbulence is modeled by a general relativistic formulation of Supernova Turbulence in Reduced dimensionality (STIR). We show that the viability of the explosion is quite EOS dependent and that it best correlates with the early-time interior entropy density of the proto-neutron star. We check that this result is not progenitor dependent, although low-mass progenitors show different explosion properties, due to the different pre-collapse nuclear composition. Larger central entropies also induce more vigorous proto-neutron-star convection in our one-dimensional turbulence model, as well as a wider convective layer. △ Less

Submitted 17 January, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: Accepted to ApJ

Journal ref: ApJ 926:2 (2022)

Neutrino Echos following Black Hole Formation in Core-Collapse Supernovae

Authors: Samuel Gullin, Evan O'Connor, Jia-Shian Wang, Jeff Tseng

Abstract: During a failed core-collapse supernova, the protoneutron star eventually collapses under its own gravitational field and forms a black hole. This collapse happens quickly, on the dynamical time of the protoneutron star, $\lesssim$0.5 ms. During this collapse, barring any excessive rotation, the entire protoneutron star is accreted into the newly formed black hole. The main source of neutrinos is… ▽ More During a failed core-collapse supernova, the protoneutron star eventually collapses under its own gravitational field and forms a black hole. This collapse happens quickly, on the dynamical time of the protoneutron star, $\lesssim$0.5 ms. During this collapse, barring any excessive rotation, the entire protoneutron star is accreted into the newly formed black hole. The main source of neutrinos is now removed and the signal abruptly shuts off over this formation timescale. However, while the source of neutrinos is turned off, the arrival times at an Earth-based detector will depend on the neutrino path. We show here that a modest amount of neutrinos, emitted just prior to the black hole forming, scatter on the infalling material into our line of sight and arrive after the formation of the black hole, up to 15 ms in our model. This neutrino echo, which we characterize with Monte Carlo simulations and analytic models, has a significantly higher average energy (upwards of $\sim$ 50 MeV) compared to the main neutrino signal, and for the canonical failed supernova explored here, is likely detectable in $\mathcal{O}$(10 kT) supernova neutrino detectors for Galactic failed supernovae. The presence of this signal is important to consider if using black hole formation as a time post for triangulation or the post black hole timing profile for neutrino mass measurements. On its own, it can also be used to characterize or constrain the structure and nature of the accretion flow. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 16 pages, 5 figures, 1 table, 5 neutrino-echo neutrinos, submitted to AAS journals

Non-Radial Neutrino Emission upon Black Hole Formation in Core Collapse Supernovae

Authors: Jia-Shian Wang, Jeff Tseng, Samuel Gullin, Evan P. O'Connor

Abstract: Black hole formation in a core-collapse supernova is expected to lead to a distinctive, abrupt drop in neutrino luminosity due to the engulfment of the main neutrino-producing regions as well as the strong gravitational redshift of those remaining neutrinos which do escape. Previous analyses of the shape of the cut-off have focused on specific trajectories or simplified models of bulk neutrino tra… ▽ More Black hole formation in a core-collapse supernova is expected to lead to a distinctive, abrupt drop in neutrino luminosity due to the engulfment of the main neutrino-producing regions as well as the strong gravitational redshift of those remaining neutrinos which do escape. Previous analyses of the shape of the cut-off have focused on specific trajectories or simplified models of bulk neutrino transport. In this article, we integrate over simple "ballistic" geodesics to investigate potential effects on the cut-off profile of including all neutrino emission angles from a collapsing surface in the Schwarzschild metric, and from a contracting equatorial mass ring in the Kerr metric. We find that the non-radial geodesics contribute to a softening of the cut-off in both cases. In addition, extreme rotation introduces significant changes to the shape of the tail which may be observable in future neutrino detectors, or combinations of detectors. △ Less

Submitted 27 October, 2021; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: 11 pages, 14 figures, minor text and format changes, final version approved for Phys Rev D

SNEWPY: A Data Pipeline from Supernova Simulations to Neutrino Signals

Authors: Amanda L. Baxter, Segev BenZvi, Joahan Castaneda Jaimes, Alexis Coleiro, Marta Colomer Molla, Damien Dornic, Tomer Goldhagen, Anne M. Graf, Spencer Griswold, Alec Habig, Remington Hill, Shunsaku Horiuchi James P. Kneller Rafael F. Lang, Massimiliano Lincetto, Jost Migenda, Ko Nakamura, Evan O'Connor, Andrew Renshaw, Kate Scholberg, Navya Uberoi, Arkin Worlikar

Abstract: Current neutrino detectors will observe hundreds to thousands of neutrinos from a Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is cur… ▽ More Current neutrino detectors will observe hundreds to thousands of neutrinos from a Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is currently a large gap between supernova simulations and the corresponding signals in neutrino detectors, which will make any comparison between theory and observation very difficult. SNEWPY is an open-source software package which bridges this gap. The SNEWPY code can interface with supernova simulation data to generate from the model either a time series of neutrino spectral fluences at Earth, or the total time-integrated spectral fluence. Data from several hundred simulations of core-collapse, thermonuclear, and pair-instability supernovae is included in the package. This output may then be used by an event generator such as sntools or an event rate calculator such as SNOwGLoBES. Additional routines in the SNEWPY package automate the processing of the generated data through the SNOwGLoBES software and collate its output into the observable channels of each detector. In this paper we describe the contents of the package, the physics behind SNEWPY, the organization of the code, and provide examples of how to make use of its capabilities. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: Software available at https://github.com/SNEWS2/snewpy

Equation of State Dependence of Gravitational Waves in Core-Collapse Supernovae

Authors: Oliver Eggenberger Andersen, Shuai Zha, André da Silva Schneider, Aurore Betranhandy, Sean M. Couch, Evan P. O'Connor

Abstract: Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars (NSs) and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto-neutron star (PNS) and the post-bounce dynamics in CCSNe simulations, subsequently impacting the GW emission.… ▽ More Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars (NSs) and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto-neutron star (PNS) and the post-bounce dynamics in CCSNe simulations, subsequently impacting the GW emission. We perform axisymmetric simulations of CCSNe with Skyrme-type EOSs to study how the GW signal and PNS convection zone are impacted by two experimentally accessible EOS parameters, (1) the effective mass of nucleons, $m^\star$, which is crucial in setting the thermal dependence of the EOS, and (2) the isoscalar incompressibility modulus, $K_{\rm{sat}}$. While $K_{\rm{sat}}$ shows little impact, the peak frequency of the GWs has a strong effective mass dependence due to faster contraction of the PNS for higher values of $m^\star$ owing to a decreased thermal pressure. These more compact PNSs also exhibit more neutrino heating which drives earlier explosions and correlates with the GW amplitude via accretion plumes striking the PNS, exciting the oscillations. We investigate the spatial origin of the GWs and show the agreement between a frequency-radial distribution of the GW emission and a perturbation analysis. We do not rule out overshoot from below via PNS convection as another moderately strong excitation mechanism in our simulations. We also study the combined effect of effective mass and rotation. In all our simulations we find evidence for a power gap near $\sim$1250 Hz, we investigate its origin and report its EOS dependence. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: 21 pages, 14 figures

Progenitor Dependence of Hadron-quark Phase Transition in Failing Core-collapse Supernovae

Authors: Shuai Zha, Evan P. O'Connor, André da Silva Schneider

Abstract: We study the consequences of a hadron-quark phase transition (PT) in failing core-collapse supernovae (CCSNe) which give birth to stellar-mass black holes (BH). We perform a suite of neutrino-transport general-relativistic hydrodynamic simulations in spherical symmetry with 21 progenitor models and a hybrid equation of state (EoS) including hadrons and quarks. We find that the effect of the PT on… ▽ More We study the consequences of a hadron-quark phase transition (PT) in failing core-collapse supernovae (CCSNe) which give birth to stellar-mass black holes (BH). We perform a suite of neutrino-transport general-relativistic hydrodynamic simulations in spherical symmetry with 21 progenitor models and a hybrid equation of state (EoS) including hadrons and quarks. We find that the effect of the PT on the CCSN postbounce dynamics is a function of the bounce compactness parameter $ξ_{2.2}$. For $ξ_{2.2}\gtrsim0.24$, the PT leads to a second dynamical collapse of the protocompact star (PCS). While BH formation starts immediately after this second collapse for models with $ξ_{2.2}\gtrsim0.51$, the PCS experiences a second bounce and oscillations for models with $0.24\lesssimξ_{2.2}\lesssim0.51$. These models emit potent oscillatory neutrino signals with a period of $\sim$ms for tens of ms after the second bounce, which can be a strong indicator of the PT in failing CCSNe if detected in the future. However, no shock revival occurs and BH formation inevitably takes place in our spherically-symmetric simulations. Furthermore, via a diagram of mass-specific entropy evolution of the PCS, the progenitor dependence can be understood through the appearance of third-family of compact stars emerging at large entropy induced by the PT. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 12 pages, 10 figures. Accepted by ApJ

Calibration of the Advanced Spectral Leakage scheme for neutron star merger simulations, and extension to smoothed-particle hydrodynamics

Authors: Davide Gizzi, Christoffer Lundman, Evan O'Connor, Stephan Rosswog, Albino Perego

Abstract: We calibrate a neutrino transport approximation, called Advanced Spectral Leakage (ASL), with the purpose of modeling neutrino-driven winds in neutron star mergers. Based on a number of snapshots we gauge the ASL parameters by comparing against both the two-moment (M1) scheme implemented in the FLASH code and the Monte Carlo neutrino code Sedonu. The ASL scheme contains three parameters, the least… ▽ More We calibrate a neutrino transport approximation, called Advanced Spectral Leakage (ASL), with the purpose of modeling neutrino-driven winds in neutron star mergers. Based on a number of snapshots we gauge the ASL parameters by comparing against both the two-moment (M1) scheme implemented in the FLASH code and the Monte Carlo neutrino code Sedonu. The ASL scheme contains three parameters, the least robust of which results to be a blocking parameter for electron neutrinos and anti-neutrinos. The parameter steering the angular distribution of neutrino heating is re-calibrated compared to the earlier work. We also present a new, fast and mesh-free algorithm for calculating spectral optical depths, which, when using Smoothed Particle Hydrodynamics (SPH), makes the neutrino transport completely particle-based. We estimate a speed-up of a factor of 100 in the optical depth calculation when comparing to a grid-based approach. In the suggested calibration we recover luminosities and mean energies within 25%. A comparison of the rates of change of internal energy and electron fraction in the neutrino-driven wind suggests comparable accuracies of ASL and M1, but a higher computational efficiency of the ASL scheme. We estimate that the ratio between the CPU hours spent on the ASL neutrino scheme and those spent on the hydrodynamics is 0.8 per timestep when considering the SPH code MAGMA2 as source code for the Lagrangian hydrodynamics, to be compared with a factor of 10 from the M1 in FLASH. △ Less

Submitted 19 May, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: Minor revision, accepted for publication in MNRAS. 20 pages, 11 figures

General Relativistic Neutrino-Driven Turbulence in One-Dimensional Core-Collapse Supernovae

Authors: Luca Boccioli, Grant J. Mathews, Evan O'Connor

Abstract: Convection and turbulence in core-collapse supernovae (CCSNe) are inherently three-dimensional in nature. However, 3D simulations of CCSNe are computationally demanding. Thus, it is valuable to modify simulations in spherical symmetry to incorporate 3D effects using some parametric model. In this paper, we report on the formulation and implementation of general relativistic (GR) neutrino-driven tu… ▽ More Convection and turbulence in core-collapse supernovae (CCSNe) are inherently three-dimensional in nature. However, 3D simulations of CCSNe are computationally demanding. Thus, it is valuable to modify simulations in spherical symmetry to incorporate 3D effects using some parametric model. In this paper, we report on the formulation and implementation of general relativistic (GR) neutrino-driven turbulent convection in the spherically symmetric core-collapse supernova code \texttt{GR1D}. This is based upon the recently proposed method of Supernova Turbulence in Reduced-dimensionality (\textit{STIR}) in Newtonian simulations from \cite{Couch2020_STIR}. When the parameters of this model are calibrated to 3D simulations, we find that our GR formulation of \textit{STIR} requires larger turbulent eddies to achieve a shock evolution similar to the original \textit{STIR} model. We also find that general relativity may alter the correspondence between progenitor mass and successful vs.~failed explosions. △ Less

Submitted 12 February, 2021; originally announced February 2021.

Comments: 15 pages, 4 figures

Journal ref: ApJ 912:29 (2021)

Connecting the Light Curves of Type IIP Supernovae to the Properties of their Progenitors

Authors: Brandon L. Barker, Chelsea E. Harris, MacKenzie L. Warren, Evan P. O'Connor, Sean M. Couch

Abstract: Observations of core-collapse supernovae (CCSNe) reveal a wealth of information about the dynamics of the supernova ejecta and its composition but very little direct information about the progenitor. Constraining properties of the progenitor and the explosion requires coupling the observations with a theoretical model of the explosion. Here, we begin with the CCSN simulations of Couch et al 2020 A… ▽ More Observations of core-collapse supernovae (CCSNe) reveal a wealth of information about the dynamics of the supernova ejecta and its composition but very little direct information about the progenitor. Constraining properties of the progenitor and the explosion requires coupling the observations with a theoretical model of the explosion. Here, we begin with the CCSN simulations of Couch et al 2020 ApJ 890 127, which use a non-parametric treatment of the neutrino transport while also accounting for turbulence and convection. In this work we use the SuperNova Explosion Code to evolve the CCSN hydrodynamics to later times and compute bolometric light curves. Focusing on SNe IIP, we then (1) directly compare the theoretical STIR explosions to observations and (2) assess how properties of the progenitor's core can be estimated from optical photometry in the plateau phase alone. First, the distribution of plateau luminosities (L$_{50}$) and ejecta velocities achieved by our simulations is similar to the observed distributions. Second, we fit our models to the light curves and velocity evolution of some well-observed SNe. Third, we recover well-known correlations, as well as the difficulty of connecting any one SN property to zero-age main sequence mass. Finally, we show that there is a usable, linear correlation between iron core mass and L$_{50}$ such that optical photometry alone of SNe IIP can give us insights into the cores of massive stars. Illustrating this by application to a few SNe, we find iron core masses of 1.3-1.5 solar masses with typical errors of ~0.05 solar masses. Data are publicly available online (\url{https://doi.org/10.5281/zenodo.6631964}). △ Less

Submitted 10 June, 2022; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 23 pages, Accepted to ApJ. Data available online https://doi.org/10.5281/zenodo.6631964

Measuring the distance and mass of galactic core-collapse supernovae using neutrinos

Authors: Manne Segerlund, Erin O'Sullivan, Evan O'Connor

Abstract: Neutrinos from a Galactic core-collapse supernova will be measured by neutrino detectors minutes to days before an optical signal reaches Earth. We present a novel calculation showing the ability of current and near-future neutrino detectors to make fast predictions of the progenitor distance and place constraints on the zero-age main sequence mass in order to inform the observing strategy for ele… ▽ More Neutrinos from a Galactic core-collapse supernova will be measured by neutrino detectors minutes to days before an optical signal reaches Earth. We present a novel calculation showing the ability of current and near-future neutrino detectors to make fast predictions of the progenitor distance and place constraints on the zero-age main sequence mass in order to inform the observing strategy for electromagnetic follow-up. We show that for typical Galactic supernovae, the distance can be constrained with an uncertainty of $\sim$5\% using IceCube or Hyper-K and, furthermore, the zero-age main sequence mass can be constrained for extremal values of compactness. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 7 pages + 6 pages of supplemental information, submitted to PRL, data and analysis scripts at https://doi.org/10.5281/zenodo.4467495

Supernova Model Discrimination with Hyper-Kamiokande

Authors: Hyper-Kamiokande Collaboration, :, K. Abe, P. Adrich, H. Aihara, R. Akutsu, I. Alekseev, A. Ali, F. Ameli, I. Anghel, L. H. V. Anthony, M. Antonova, A. Araya, Y. Asaoka, Y. Ashida, V. Aushev, F. Ballester, I. Bandac, M. Barbi, G. J. Barker, G. Barr, M. Batkiewicz-Kwasniak, M. Bellato, V. Berardi, M. Bergevin , et al. (478 additional authors not shown)

Abstract: Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-colla… ▽ More Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-collapse supernovae is not yet well understood. Hyper-Kamiokande is a next-generation neutrino detector that will be able to observe the neutrino flux from the next galactic core-collapse supernova in unprecedented detail. We focus on the first 500 ms of the neutrino burst, corresponding to the accretion phase, and use a newly-developed, high-precision supernova event generator to simulate Hyper-Kamiokande's response to five different supernova models. We show that Hyper-Kamiokande will be able to distinguish between these models with high accuracy for a supernova at a distance of up to 100 kpc. Once the next galactic supernova happens, this ability will be a powerful tool for guiding simulations towards a precise reproduction of the explosion mechanism observed in nature. △ Less

Submitted 20 July, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 21 pages, 7 figures. Article based on thesis published as arXiv:2002.01649. v2: added references and some explanations in response to reviewer comments

Journal ref: Astrophys.J. 916 (2021) 15

Determining the Structure of Rotating Massive Stellar Cores with Gravitational Waves

Authors: Michael A. Pajkos, MacKenzie L. Warren, Sean M. Couch, Evan P. O'Connor, Kuo-Chuan Pan

Abstract: The gravitational wave (GW) signal resulting from stellar core collapse encodes a wealth of information about the physical parameters of the progenitor star and the resulting core-collapse supernova (CCSN). We present a novel approach to constrain CCSN progenitor properties at collapse using two of the most detectable parts of the GW signal: the core-bounce signal and evolution of the dominant fre… ▽ More The gravitational wave (GW) signal resulting from stellar core collapse encodes a wealth of information about the physical parameters of the progenitor star and the resulting core-collapse supernova (CCSN). We present a novel approach to constrain CCSN progenitor properties at collapse using two of the most detectable parts of the GW signal: the core-bounce signal and evolution of the dominant frequency mode from the protoneutron star. We focus on the period after core bounce but before explosion and investigate the predictive power of GWs from rotating CCSNe to constrain properties of the progenitor star. We analyze 34 2D and four 3D neutrinoradiation-hydrodynamic simulations of stellar core collapse in progenitors of varied initial mass and rotation rate. Extending previous work, we verify the compactness of the progenitor at collapse to correlate with the early ramp-up slope, and in rotating cases, also with the core angular momentum. Combining this information with the bounce signal, we present a new analysis method to constrain the pre-collapse core compactness of the progenitor. Because these GW features occur less than a second after core bounce, this analysis could allow astronomers to predict electromagnetic properties of a resulting CCSN even before shock breakout. △ Less

Submitted 22 June, 2021; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 21 pages, 11 figures; Published in ApJ

Journal ref: ApJ 914:80 (2021)

SNEWS 2.0: A Next-Generation SuperNova Early Warning System for Multi-messenger Astronomy

Authors: S. Al Kharusi, S. Y. BenZvi, J. S. Bobowski, W. Bonivento, V. Brdar, T. Brunner, E. Caden, M. Clark, A. Coleiro, M. Colomer-Molla, J. I. Crespo-Anadón, A. Depoian, D. Dornic, V. Fischer, D. Franco, W. Fulgione, A. Gallo Rosso, M. Geske, S. Griswold, M. Gromov, D. Haggard, A. Habig, O. Halim, A. Higuera, R. Hill , et al. (46 additional authors not shown)

Abstract: The next core-collapse supernova in the Milky Way or its satellites will represent a once-in-a-generation opportunity to obtain detailed information about the explosion of a star and provide significant scientific insight for a variety of fields because of the extreme conditions found within. Supernovae in our galaxy are not only rare on a human timescale but also happen at unscheduled times, so i… ▽ More The next core-collapse supernova in the Milky Way or its satellites will represent a once-in-a-generation opportunity to obtain detailed information about the explosion of a star and provide significant scientific insight for a variety of fields because of the extreme conditions found within. Supernovae in our galaxy are not only rare on a human timescale but also happen at unscheduled times, so it is crucial to be ready and use all available instruments to capture all possible information from the event. The first indication of a potential stellar explosion will be the arrival of a bright burst of neutrinos. Its observation by multiple detectors worldwide can provide an early warning for the subsequent electromagnetic fireworks, as well as signal to other detectors with significant backgrounds so they can store their recent data. The Supernova Early Warning System (SNEWS) has been operating as a simple coincidence between neutrino experiments in automated mode since 2005. In the current era of multi-messenger astronomy there are new opportunities for SNEWS to optimize sensitivity to science from the next Galactic supernova beyond the simple early alert. This document is the product of a workshop in June 2019 towards design of SNEWS 2.0, an upgraded SNEWS with enhanced capabilities exploiting the unique advantages of prompt neutrino detection to maximize the science gained from such a valuable event. △ Less

Submitted 18 March, 2021; v1 submitted 30 October, 2020; originally announced November 2020.

Comments: 66 pages, 11 figures

Journal ref: New J. Phys. 23 031201 (2021)

Enhanced Lidov-Kozai migration and the formation of the transiting giant planet WD1856+534b

Authors: Christopher E. O'Connor, Bin Liu, Dong Lai

Abstract: We investigate the possible origin of the transiting giant planet WD1856+534b, the first strong exoplanet candidate orbiting a white dwarf, through high-eccentricity migration (HEM) driven by the Lidov-Kozai (LK) effect. The host system's overall architecture is an hierarchical quadruple in the '2+2' configuration, owing to the presence of a tertiary companion system of two M-dwarfs. We show that… ▽ More We investigate the possible origin of the transiting giant planet WD1856+534b, the first strong exoplanet candidate orbiting a white dwarf, through high-eccentricity migration (HEM) driven by the Lidov-Kozai (LK) effect. The host system's overall architecture is an hierarchical quadruple in the '2+2' configuration, owing to the presence of a tertiary companion system of two M-dwarfs. We show that a secular inclination resonance in 2+2 systems can significantly broaden the LK window for extreme eccentricity excitation ($e \gtrsim 0.999$), allowing the giant planet to migrate for a wide range of initial orbital inclinations. Octupole effects can also contribute to the broadening of this 'extreme' LK window. By requiring that perturbations from the companion stars be able to overcome short-range forces and excite the planet's eccentricity to $e \simeq 1$, we obtain an absolute limit of $a_{1} \gtrsim 8 \, {\rm AU} \, (a_{3} / 1500 \, {\rm AU})^{6/7}$ for the planet's semi-major axis just before migration (where $a_{3}$ is the semi-major axis of the 'outer' orbit). We suggest that, to achieve a wide LK window through the 2+2 resonance, WD1856b likely migrated from $30 \, {\rm AU} \lesssim a_{1} \lesssim 60 \, {\rm AU}$, corresponding to $\sim 10$--$20 \, {\rm AU}$ during the host's main-sequence phase. We discuss possible difficulties of all flavours of HEM affecting the occurrence rate of short-period giant planets around white dwarfs. △ Less

Submitted 11 January, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 8 pages, 4 figures, 1 table; MNRAS 501 (2021) 507-514

Impact of neutrino pair-production rates in Core-Collapse Supernovae

Authors: Aurore Betranhandy, Evan O'Connor

Abstract: In this paper, we present a careful study on the impact of neutrino pair-production on core-collapse supernovae via spherically-symmetric, general-relativistic simulations of two different massive star progenitors with energy-dependent neutrino transport. We explore the impact and consequences of both the underlying microphysics and the implementation in the radiation transport algorithms on the s… ▽ More In this paper, we present a careful study on the impact of neutrino pair-production on core-collapse supernovae via spherically-symmetric, general-relativistic simulations of two different massive star progenitors with energy-dependent neutrino transport. We explore the impact and consequences of both the underlying microphysics and the implementation in the radiation transport algorithms on the supernova evolution, neutrino signal properties, and the explosion dynamics. We consider the two dominant neutrino pair-production processes found in supernovae, electron-positron annihilation as well as nucleon-nucleon bremsstrahlung in combination with both a simplified and a complete treatment of the processes in the radiation transport algorithms. We find that the use of the simplified prescription quantitatively impacts the neutrino signal at the 10\% level and potentially the supernova dynamics, as we show for the case of a 9.6M$_\odot$ progenitor. We also show that the choice of nucleon-nucleon bremsstrahlung interaction can also have a quantitative impact on the neutrino signal. A self-consistent treatment with state-of-the-art microphysics is suggested for precision simulations of core collapse, however the simplified treatment explored here is both computationally less demanding and results in a qualitatively similar evolution. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 15 pages, 7 figures, submitted to PRD

Journal ref: Phys. Rev. D 102, 123015 (2020)

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

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

Abstract: 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 amplitud… ▽ More 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. △ Less

Submitted 29 October, 2021; v1 submitted 9 July, 2020; originally announced July 2020.

Comments: 6 pages, 4 figures, with Supplementary Materials. Corrected characteristic strain spectra, Fig. 4 and Fig. S3(b). Conclusion of original version is not affected. See also the PRL erratum

High-eccentricity migration of planetesimals around polluted white dwarfs

Authors: Christopher E. O'Connor, Dong Lai

Abstract: Several white dwarfs with atmospheric metal pollution have been found to host small planetary bodies (planetesimals) orbiting near the tidal disruption radius. We study the physical properties and dynamical origin of these bodies under the hypothesis that they underwent high-eccentricity migration from initial distances of several astronomical units. We examine two plausible mechanisms for orbital… ▽ More Several white dwarfs with atmospheric metal pollution have been found to host small planetary bodies (planetesimals) orbiting near the tidal disruption radius. We study the physical properties and dynamical origin of these bodies under the hypothesis that they underwent high-eccentricity migration from initial distances of several astronomical units. We examine two plausible mechanisms for orbital migration and circularization: tidal friction and ram-pressure drag in a compact disc. For each mechanism, we derive general analytic expressions for the evolution of the orbit that can be rescaled for various situations. We identify the physical parameters that determine whether a planetesimal's orbit can circularize within the appropriate time-scale and constrain these parameters based on the properties of the observed systems. For tidal migration to work, an internal viscosity similar to that of molten rock is required, and this may be naturally produced by tidal heating. For disc migration to operate, a minimal column density of the disc is implied; the inferred total disc mass is consistent with estimates of the total mass of metals accreted by polluted WDs. △ Less

Submitted 8 October, 2020; v1 submitted 12 May, 2020; originally announced May 2020.

Comments: 17 pages, 13 figures, 2 appendices

Journal ref: MNRAS 498 (2020) 4005-4020

Equation of State and Progenitor Dependence of Stellar-Mass Black-Hole Formation

Authors: A. S. Schneider, E. O'Connor, E. Granqvist, A. Betranhandy, S. M. Couch

Abstract: The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS it will become gravitationally unstable and further collapse into a black-hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of presupernova progenitor stars and equations of state (EOSs) of nuclear ma… ▽ More The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS it will become gravitationally unstable and further collapse into a black-hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of presupernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core-collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-entropy phase space $M_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and almost EOS independent, apart from the final endpoint. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum \textit{gravitational} mass supported by a hot NS with constant entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: 22 pages, 13 figures, 3 tables

Constraining properties of the next nearby core-collapse supernova with multi-messenger signals

Authors: MacKenzie L. Warren, Sean M. Couch, Evan P. O'Connor, Viktoriya Morozova

Abstract: With the advent of modern neutrino and gravitational wave detectors, the promise of multi-messenger detections of the next galactic core-collapse supernova has become very real. Such detections will give insight into the core-collapse supernova mechanism, the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detect… ▽ More With the advent of modern neutrino and gravitational wave detectors, the promise of multi-messenger detections of the next galactic core-collapse supernova has become very real. Such detections will give insight into the core-collapse supernova mechanism, the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detections, a thorough understanding of the landscape of possible core-collapse supernova events, and their multi-messenger signals, is needed. We present detailed predictions of neutrino and gravitational wave signals from 1D simulations of stellar core collapse, spanning the landscape of core-collapse progenitors from $9-120\,\mathrm{M}_{\odot}$. In order to achieve explosions in 1D, we use the STIR model, which includes the effects of turbulence and convection in 1D supernova simulations to mimic the 3D explosion mechanism. We study the gravitational wave emission from the 1D simulations using an astroseismology analysis of the proto-neutron star. We find that the neutrino and gravitational wave signals are strongly correlated with the structure of the progenitor star and remnant compact object. Using these correlations, future detections of the first few seconds of neutrino and gravitational wave emission from a galactic core-collapse supernova may be able to provide constraints on stellar evolution independent of pre-explosion imaging and the mass of the compact object remnant prior to fallback accretion. △ Less

Submitted 12 May, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 25 pages, 15 figures, 8 tables, 462 explosions, 2.7e6 neutrinos. Submitted to AAS Journals

Multimessenger Asteroseismology of Core-Collapse Supernovae

Authors: John Ryan Westernacher-Schneider, Evan O'Connor, Erin O'Sullivan, Irene Tamborra, Meng-Ru Wu, Sean M. Couch, Felix Malmenbeck

Abstract: We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutri… ▽ More We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $\sim 420$ Hz, and yet the mode's frequency in our simulations is $\sim 15\%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $\sim\!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $\sim\!5$ kpc using Advanced LIGO at design sensitivity. △ Less

Submitted 8 November, 2019; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: 25 pages, 16 figures. Version accepted to PRD. Error corrected, new appendices and figures added

Report number: 15 December 2019 issue of Physical Review D (Vol. 100, No. 12)

Journal ref: Phys. Rev. D 100, 123009 (2019)

A multi-dimensional implementation of the Advanced Spectral neutrino Leakage scheme

Authors: Davide Gizzi, Evan O'Connor, Stephan Rosswog, Albino Perego, Ruben Cabezón, Lorenzo Nativi

Abstract: We present a new, multi-dimensional implementation of the Advanced Spectral Leakage (ASL) scheme with the purpose of modelling neutrino-matter interactions in neutron star mergers. A major challenge is the neutrino absorption in the semi-transparent regime, which is responsible for driving winds from the merger remnant. The composition of such winds is crucial in the understanding of the electroma… ▽ More We present a new, multi-dimensional implementation of the Advanced Spectral Leakage (ASL) scheme with the purpose of modelling neutrino-matter interactions in neutron star mergers. A major challenge is the neutrino absorption in the semi-transparent regime, which is responsible for driving winds from the merger remnant. The composition of such winds is crucial in the understanding of the electromagnetic emission in the recently observed macronova following GW170817. Compared to the original version, we introduce an optical-depth-dependent flux factor to model the average angle of neutrino propagation, and a modulation that accounts for flux anisotropies in non-spherical geometries. We scrutinise our approach by first comparing the new scheme against the original one for a spherically symmetric core-collapse supernova snapshot, both in 1D and in 3D, and additionally against a two-moment (M1) scheme as implemented in 1D into the code GR1D. The luminosities and mean energies agree to a few percents in most tests. Finally, for the case of a binary merger remnant snapshot we compare the new ASL scheme with the M1 scheme that is implemented in the Eulerian adaptive mesh refinement code FLASH. We find that the neutrino absorption distribution in the semi-transparent regime is overall well reproduced. Both approaches agree to within $\lesssim 15\%$ for the average energies and to better than $\sim 35 \%$ in the total luminosities. △ Less

Submitted 15 October, 2019; v1 submitted 27 June, 2019; originally announced June 2019.

Comments: v.2, minor corrections, accepted for publication in MNRAS

Equation of state effects in the core collapse of a $20$-$M_\odot$ star

Authors: A. S. Schneider, L. F. Roberts, C. D. Ott, E. O'connor

Abstract: Uncertainties in our knowledge of the properties of dense matter near and above nuclear saturation density are among the main sources of variations in multi-messenger signatures predicted for core-collapse supernovae (CCSNe) and the properties of neutron stars (NSs). We construct 97 new finite-temperature equations of state (EOSs) of dense matter that obey current experimental, observational, and… ▽ More Uncertainties in our knowledge of the properties of dense matter near and above nuclear saturation density are among the main sources of variations in multi-messenger signatures predicted for core-collapse supernovae (CCSNe) and the properties of neutron stars (NSs). We construct 97 new finite-temperature equations of state (EOSs) of dense matter that obey current experimental, observational, and theoretical constraints and discuss how systematic variations in the EOS parameters affect the properties of cold nonrotating NSs and the core collapse of a $20\,M_\odot$ progenitor star. The core collapse of the $20\,M_\odot$ progenitor star is simulated in spherical symmetry using the general-relativistic radiation-hydrodynamics code GR1D where neutrino interactions are computed for each EOS using the NuLib library. We conclude that the effective mass of nucleons at densities above nuclear saturation density is the largest source of uncertainty in the CCSN neutrino signal and dynamics even though it plays a subdominant role in most properties of cold NS matter. Meanwhile, changes in other observables affect the properties of cold NSs, while having little effect in CCSNe. To strengthen our conclusions, we perform six octant three-dimensional CCSN simulations varying the effective mass of nucleons at nuclear saturation density. We conclude that neutrino heating and, thus, the likelihood of explosion is significantly increased for EOSs where the effective mass of nucleons at nuclear saturation density is large. △ Less

Submitted 5 June, 2019; originally announced June 2019.

Comments: 21 pages, 13 figures, 2 tables

Journal ref: Phys. Rev. C 100, 055802 (2019)

Characteristics of ephemeral coronal holes

Authors: Andrew R. Inglis, Rachel E. O'Connor, W. Dean Pesnell, Michael S. Kirk, Nishu Karna

Abstract: Small-scale ephemeral coronal holes may be a recurring feature on the solar disk, but have received comparatively little attention. These events are characterized by compact structure and short total lifetimes, substantially less than a solar disk crossing. We present a systematic search for these events, using Atmospheric Imaging Assembly EUV image data from the Solar Dynamics Observatory, coveri… ▽ More Small-scale ephemeral coronal holes may be a recurring feature on the solar disk, but have received comparatively little attention. These events are characterized by compact structure and short total lifetimes, substantially less than a solar disk crossing. We present a systematic search for these events, using Atmospheric Imaging Assembly EUV image data from the Solar Dynamics Observatory, covering the time period 2010 - 2015. Following strict criteria, this search yielded four clear examples of the ephemeral coronal hole phenomenon. The properties of each event are characterized, including their total lifetime, growth and decay rates, and areas. The magnetic properties of these events are also determined using Helioseismic and Magnetic Imager data. Based on these four events, ephemeral coronal holes experience rapid initial growth of up to 3000 Mm2/hr, while the decay phases are typically more gradual. Like conventional coronal holes, the mean magnetic field in each ephemeral coronal hole displays a consistent polarity, with mean magnetic flux densities generally < 10 G. No evidence of a corresponding signature is seen in solar wind data at 1 AU. Further study is needed to determine whether ephemeral coronal holes are under-reported events or a truly rare phenomenon. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: Accepted for publication in the Astrophysical Journal, 9 pages, 9 figures

Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s

Authors: R. J. Foley, K. D. Alexander, I. Andreoni, I. Arcavi, K. Auchettl, J. Barnes, G. Baym, E. C. Bellm, A. M. Beloborodov, N. Blagorodnova, J. P. Blakeslee, P. R. Brady, M. Branchesi, J. S. Brown, N. Butler, M. Cantiello, R. Chornock, D. O. Cook, J. Cooke, D. L. Coppejans, A. Corsi, S. M. Couch, M. W. Coughlin, D. A. Coulter, P. S. Cowperthwaite , et al. (88 additional authors not shown)

Abstract: As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse s… ▽ More As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020

Evolution of the progenitors of SNe 1993J and 2011dh revealed through late time radio and X-ray studies

Authors: Esha Kundu, Peter Lundqvist, Elena Sorokina, Miguel A. Pérez-Torres, Sergei Blinnikov, Evan O'Connor, Mattias Ergon, Poonam Chandra, Barnali Das

Abstract: We perform hydrodynamical simulations of the interaction between supernova (SN) ejecta and circumstellar medium (CSM) for SN 1993J and SN 2011dh, and calculate the radio and X-ray emissions expected from the shocked gas at late epochs ($t$). Considering the ejecta structure from multi-group radiation hydrodynamics simulation, we find that the observed rapid drop in radio and X-ray light curves of… ▽ More We perform hydrodynamical simulations of the interaction between supernova (SN) ejecta and circumstellar medium (CSM) for SN 1993J and SN 2011dh, and calculate the radio and X-ray emissions expected from the shocked gas at late epochs ($t$). Considering the ejecta structure from multi-group radiation hydrodynamics simulation, we find that the observed rapid drop in radio and X-ray light curves of SN 1993J at $t>$3000 days can be due to a change in the mass-loss rate ($\dot M$) around $\sim$6500 years prior to the explosion of the SN. The exact epoch scales inversely with the assumed wind velocity of $v_{\rm w}=10~ km~s^{-1}$. The progenitor of this SN very likely belonged to a binary system, where, during its evolution, the primary had transferred material to the secondary. It is argued in the paper that the change in $\dot M$ can happen because of a change in the mass accretion efficiency ($η$) of the companion star. It is possible that before $\sim6500~(v_{\rm w}/10~km~s^{-1})^{-1}$years prior to the explosion, $η$ was high, thus the CSM was tenuous, which causes the late time downturn in fluxes. In the case of SN 2011dh, the late time evolution is found to be consistent with a wind medium with $\dot M/v_{\rm w}=4\times10^{-6 }~M_{\odot}~ yr^{-1}/10 ~{km ~s^{-1}}$. It is difficult from our analysis to predict whether the progenitor of this SN had a binary companion, however, if future observations show similar decrease in radio and X-ray fluxes, then that would give strong support to a scenario where both SNe had undergone similar kind of binary evolution before explosion. △ Less

Submitted 6 March, 2019; originally announced March 2019.

Comments: 13 pages, 10 figures; Accepted for publication in ApJ

Simulating Turbulence-aided Neutrino-driven Core-collapse Supernova Explosions in One Dimension

Authors: Sean M. Couch, MacKenzie L. Warren, Evan P. O'Connor

Abstract: The core-collapse supernova (CCSN) mechanism is fundamentally three-dimensional with instabilities, convection, and turbulence playing crucial roles in aiding neutrino-driven explosions. Simulations of CCNSe including accurate treatments of neutrino transport and sufficient resolution to capture key instabilities remain amongst the most expensive numerical simulations in astrophysics, prohibiting… ▽ More The core-collapse supernova (CCSN) mechanism is fundamentally three-dimensional with instabilities, convection, and turbulence playing crucial roles in aiding neutrino-driven explosions. Simulations of CCNSe including accurate treatments of neutrino transport and sufficient resolution to capture key instabilities remain amongst the most expensive numerical simulations in astrophysics, prohibiting large parameter studies in 2D and 3D. Studies spanning a large swath of the incredibly varied initial conditions of CCSNe are possible in 1D, though such simulations must be artificially driven to explode. We present a new method for including the most important effects of convection and turbulence in 1D simulations of neutrino-driven CCSNe, called Supernova Turbulence In Reduced-dimensionality, or STIR. Our new approach includes crucial terms resulting from the turbulent and convective motions of the flow. We estimate the strength of convection and turbulence using a modified mixing length theory (MLT) approach introducing a few free parameters to the model which are fit to the results of 3D simulations. For sufficiently large values of the mixing length parameter, turbulence-aided neutrino-driven explosions are obtained. We compare the results of STIR to high-fidelity 3D simulations and perform a parameter study of CCSN explosion using 200 solar-metallicity progenitor models from 9 to 120 $M_\odot$. We find that STIR is a better predictor of which models will explode in multidimensional simulations than other methods of driving explosions in 1D. We also present a preliminary investigation of predicted observable characteristics of the CCSN population from STIR, such as the distributions of explosion energies and remnant masses. △ Less

Submitted 25 November, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 30 pages, 17 figures. Substantially updated following referee report and resubmitted to AAS Journals

Features of Accretion-phase Gravitational-wave Emission from Two-dimensional Rotating Core-collapse Supernovae

Authors: Michael A. Pajkos, Sean M. Couch, Kuo-Chuan Pan, Evan P. O'Connor

Abstract: We explore the influence of progenitor mass and rotation on the gravitational-wave (GW) emission from core-collapse supernovae, during the postbounce, preexplosion, accretion-phase. We present the results from 15 two-dimensional (2D) neutrino radiation-hydrodynamic simulations from initial stellar collapse to $\sim$300 ms after core bounce. We examine the features of the GW signals for four zero-a… ▽ More We explore the influence of progenitor mass and rotation on the gravitational-wave (GW) emission from core-collapse supernovae, during the postbounce, preexplosion, accretion-phase. We present the results from 15 two-dimensional (2D) neutrino radiation-hydrodynamic simulations from initial stellar collapse to $\sim$300 ms after core bounce. We examine the features of the GW signals for four zero-age main sequence (ZAMS) progenitor masses ranging from 12 $M_\odot$ to 60 $M_\odot$ and four core rotation rates from 0 to 3 rad s$^{-1}$. We find that GW strain immediately around core bounce is fairly independent of ZAMS mass and---consistent with previous findings---that it is more heavily dependent on the core angular momentum. At later times, all nonrotating progenitors exhibit loud GW emission, which we attribute to vibrational g-modes of the protoneutron star (PNS) excited by convection in the postshock layer and the standing accretion shock instability (SASI). We find that increasing rotation rates results in muting of the accretion-phase GW signal due to centrifugal effects that inhibit convection in the postshock region, quench the SASI, and slow the rate at which the PNS peak vibrational frequency increases. Additionally, we verify the efficacy of our approximate general relativistic (GR) effective potential treatment of gravity by comparing our core bounce GW strains with the recent 2D GR results of other groups. △ Less

Submitted 29 June, 2019; v1 submitted 25 January, 2019; originally announced January 2019.

Comments: 16 pages, 12 figures; accepted to ApJ

Journal ref: ApJ 878:13 (2019)

FRIB and the GW170817 Kilonova

Authors: A. Aprahamian, R. Surman, A. Frebel, G. C. McLaughlin, A. Arcones, A. B. Balantekin, J. Barnes, Timothy C. Beers, Erika M. Holmbeck, Jinmi Yoon, Maxime Brodeur, T. M. Sprouse, Nicole Vassh, Jolie A. Cizewski, Jason A. Clark, Benoit Cote, Sean M. Couch, M. Eichler, Jonathan Engel, Rana Ezzeddine, George M. Fuller, Samuel A. Giuliani, Robert Grzywacz, Sophia Han, C. J. Horowitz , et al. (23 additional authors not shown)

Abstract: In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experim… ▽ More In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: Proceedings for the FRIB Theory Alliance workshop "FRIB and the GW170817 kilonova", held 16-27 July 2018 at the Facility for Rare Isotope Beams, Michigan State University, East Lansing, MI USA; 44 pages, 14 figures

Report number: LA-UR-18-28022

Exploring Fundamentally Three-dimensional Phenomena in High-fidelity Simulations of Core-collapse Supernovae

Authors: Evan O'Connor, Sean Couch

Abstract: The details of the physical mechanism that drives core-collapse supernovae (CCSNe) remain uncertain. While there is an emerging consensus on the qualitative outcome of detailed CCSN mechanism simulations in 2D, only recently have high-fidelity 3D simulations become possible. Here we present the results of an extensive set of 3D CCSN simulations using high-fidelity multidimensional neutrino transpo… ▽ More The details of the physical mechanism that drives core-collapse supernovae (CCSNe) remain uncertain. While there is an emerging consensus on the qualitative outcome of detailed CCSN mechanism simulations in 2D, only recently have high-fidelity 3D simulations become possible. Here we present the results of an extensive set of 3D CCSN simulations using high-fidelity multidimensional neutrino transport, high-resolution hydrodynamics, and approximate general relativistic gravity. We employ a state-of-the-art 20 solar mass progenitor generated using the Modules for Experiments in Stellar Astrophysics (MESA; Farmer et al. (2016) Paxton et al. (2011, 2013, 2015, 2018) and the SFHo equation of state of Steiner et al. (2013). While none of our 3D CCSN simulations explode within ~500ms after core bounce, we find that the presence of large scale aspherical motion in the Si and O shells surrounding the collapsing iron core aid shock expansion and bring the models closer to the threshold of explosion. We also find some dependence on resolution and geometry (octant vs. full 4$π$). As has been noted in other recent works, we find that the post-shock turbulence plays an important role in determining the overall dynamical evolution of our simulations. We find a strong standing accretion shock instability (SASI) that develops at late times during the shock recession epoch. The SASI aids in transient shock expansion phases, but is not enough to result in shock revival. We also report that for a subset of our simulations, we find conclusive evidence for the LESA first reported in Tamborra et al. (2014), but until now, not confirmed by other simulation codes. Both the progenitor asphericities and the SASI-induced transient shock expansion phases generate transient gravitational waves and neutrino signal modulations via perturbations of the protoneutron star by turbulent motions. (abridged) △ Less

Submitted 19 July, 2018; originally announced July 2018.

Comments: 20 pages, 14 figures, 0 explosions, submitted to The Astrophysical Journal. Comments welcome

Elastic Scattering in General Relativistic Ray Tracing for Neutrinos

Authors: M. Brett Deaton, Evan O'Connor, Y. L. Zhu, Andy Bohn, Jerred Jesse, Francois Foucart, Matthew D. Duez, G. C. McLaughlin

Abstract: We present a covariant ray tracing algorithm for computing high-resolution neutrino distributions in general relativistic numerical spacetimes with hydrodynamical sources. Our formulation treats the very important effect of elastic scattering of neutrinos off of nuclei and nucleons (changing the neutrino's direction but not energy) by incorporating estimates of the background neutrino fields. Back… ▽ More We present a covariant ray tracing algorithm for computing high-resolution neutrino distributions in general relativistic numerical spacetimes with hydrodynamical sources. Our formulation treats the very important effect of elastic scattering of neutrinos off of nuclei and nucleons (changing the neutrino's direction but not energy) by incorporating estimates of the background neutrino fields. Background fields provide information about the spectra and intensities of the neutrinos scattered into each ray. These background fields may be taken from a low-order moment simulation or be ignored, in which case the method reduces to a standard state-of-the-art ray tracing formulation. The method handles radiation in regimes spanning optically thick to optically thin. We test the new code, highlight its strengths and weaknesses, and apply it to a simulation of a neutron star merger to compute neutrino fluxes and spectra, and to demonstrate a neutrino flavor oscillation calculation. In that environment, we find qualitatively different fluxes, spectra, and oscillation behaviors when elastic scattering is included. △ Less

Submitted 26 June, 2018; originally announced June 2018.

Comments: 25 pages, submitted to Physical Review D

Journal ref: Phys. Rev. D 98, 103014 (2018)

The Impact of Different Neutrino Transport Methods on Multidimensional Core-collapse Supernova Simulations

Authors: Kuo-Chuan Pan, Carlos Mattes, Evan P. O'Connor, Sean M. Couch, Albino Perego, Almudena Arcones

Abstract: Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in… ▽ More Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in one- and two-dimensions with three approximate neutrino transport schemes, each implemented in the {\tt FLASH} simulation framework. We compare a two moment M1 scheme with an analytic closure (M1), the Isotropic Diffusion Source Approximation (IDSA), and the Advanced Spectral Leakage (ASL) method. We identify and discuss the advantages and disadvantages of each scheme. For each approximate transport scheme, we use identical grid setups, hydrodynamics, and gravity solvers to investigate the transport effects on supernova shock dynamics and neutrino quantities. We find that the transport scheme has a small effect on the evolution of protoneutron star (PNS) radius, PNS mass, and the mass accretion rate. The neutrino luminosities, mean energies, and shock radii have a $\sim$ 10-20\% quantitative difference but the overall qualitative trends are fairly consistent between all three approximations. We find larger differences in the gain region properties, including the gain region mass and the net heating rate in the gain region, as well as the strength of PNS convection in the core. We investigate the progenitor, nuclear equation of state, and stochastic perturbation dependence of our simulations and find similar magnitudes of impact on key quantities. We also compare the computational expense of the various approximations. △ Less

Submitted 15 October, 2018; v1 submitted 26 June, 2018; originally announced June 2018.

Comments: 28 pages, 11 figures, accepted by J. Phys. G focus issue on core-collapse supernovae

Global Comparison of Core-Collapse Supernova Simulations in Spherical Symmetry

Authors: Evan O'Connor, Robert Bollig, Adam Burrows, Sean Couch, Tobias Fischer, Hans-Thomas Janka, Kei Kotake, Eric J. Lentz, Matthias Liebendörfer, O. E. Bronson Messer, Anthony Mezzacappa, Tomoya Takiwaki, David Vartanyan

Abstract: We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of agreement between various groups studying the core-collapse supernova cen… ▽ More We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of agreement between various groups studying the core-collapse supernova central engine. Second, we desire to form a strong basis for future simulation codes and methods to compare to. Lastly, we want this work to be a stepping stone for future work exploring more complex simulations of core-collapse supernovae, i.e., simulations in multiple dimensions and simulations with modern neutrino and nuclear physics. We compare the early (first ~500ms after core bounce) spherically-symmetric evolution of a 20 solar mass progenitor star from six different core-collapse supernovae codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, F{\sc{ornax}}, GR1D, and PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic methods employed, we find excellent agreement in many critical quantities, including the shock radius evolution and the amount of neutrino heating. Our results provide an excellent starting point from which to extend this comparison to higher dimensions and compare the development of hydrodynamic instabilities that are crucial to the supernova explosion mechanism, such as turbulence and convection. △ Less

Submitted 11 June, 2018; originally announced June 2018.

Comments: 24 pages, 7 figures, J. Phys. G focus issue on core-collapse supernovae. This document was written collaboratively on Authorea, comments welcome at https://www.authorea.com/users/1943/articles/167397-global-comparison-of-core-collapse-supernova-simulations-in-spherical-symmetry

Black hole-neutron star mergers using a survey of finite-temperature equations of state

Authors: Wyatt Brege, Matthew D. Duez, Francois Foucart, M. Brett Deaton, Jesus Caro, Daniel A. Hemberger, Lawrence E. Kidder, Evan O'Connor, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Each of the potential signals from a black hole-neutron star merger should contain an imprint of the neutron star equation of state: gravitational waves via its effect on tidal disruption, the kilonova via its effect on the ejecta, and the gamma ray burst via its effect on the remnant disk. These effects have been studied by numerical simulations and quantified by semi-analytic formulae. However,… ▽ More Each of the potential signals from a black hole-neutron star merger should contain an imprint of the neutron star equation of state: gravitational waves via its effect on tidal disruption, the kilonova via its effect on the ejecta, and the gamma ray burst via its effect on the remnant disk. These effects have been studied by numerical simulations and quantified by semi-analytic formulae. However, most of the simulations on which these formulae are based use equations of state without finite temperature and composition-dependent nuclear physics. In this paper, we simulate black hole-neutron star mergers varying both the neutron star mass and the equation of state, using three finite-temperature nuclear models of varying stiffness. Our simulations largely vindicate formulae for ejecta properties but do not find the expected dependence of disk mass on neutron star compaction. We track the early evolution of the accretion disk, largely driven by shocking and fallback inflow, and do find notable equation of state effects on the structure of this early-time, neutrino-bright disk. △ Less

Submitted 25 April, 2018; originally announced April 2018.

Comments: 10 pages, 11 figures, submitted to Phys. Rev. D

Journal ref: Phys. Rev. D 98, 063009 (2018)

Core-Collapse Supernova Simulations including Neutrino Interactions from the Virial EOS

Authors: Evan O'Connor, C. J. Horowitz, Zidu Lin, Sean Couch

Abstract: Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the superno… ▽ More Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the supernova shock. Accurate neutrino-matter interaction terms are crucial for simulating these explosions. In this proceedings for IAUS 331, SN 1987A, 30 years later, we explore several corrections to the neutrino-nucleon scattering opacity and demonstrate the effect on the dynamics of the core-collapse supernova central engine via two dimensional neutrino-radiation-hydrodynamics simulations. Our results reveal that the explosion properties are sensitive to corrections to the neutral-current scattering cross section at the 10-20% level, but only for densities at or above $\sim 10^{12}$ g cm$^{-3}$ △ Less

Submitted 21 December, 2017; originally announced December 2017.

Comments: 6 pages, 3 figures, appears in Proc. IAU Symposium 331, SN 1987A, 30 years later - Cosmic Rays and Nuclei from Supernovae and Their Aftermaths

Journal ref: Supernova 1987A:30 years later - Cosmic Rays and Nuclei from Supernovae and their aftermaths, Proceedings of the International Astronomical Union, IAU Symposium, Volume 331, pp. 107-112