BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Sean M. Couch

Sean M. Couch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-phastro-ph.SRastro-ph.GAastro-ph.IMnucl-thphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2026arXiv

thornado+Flash-X: A Hybrid DG-IMEX and Finite-Volume Framework for Neutrino-Radiation Hydrodynamics in Core-Collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement (AMR)-based multiphysics simulation framework Flash-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $O(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino-matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino-electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved with a hybrid DG-FV representation of the fluid variables and operator-split evolution in Flash-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within Flash-X. thornado's neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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}).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Antesonic Condition for the Explosion of Core-Collapse Supernovae II: Rotation and Turbulence

In the problem of steady free-fall onto a standing shockwave around acentral mass, the &#34;antesonic&#34; condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\leq3/16$, where $c_T$ is the isothermal sound speed in the subsonic post-shock flow, and $v_\mathrm{esc}$ is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and theshock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models ofaccretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a $\sim\!5$ per-cent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by $\sim\!20-40$ per-cent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multi-dimensional simulations explode more easily than their one-dimensional counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Impact of Three Dimensions in Simulations of Neutrino-Driven Core-Collapse Supernova Explosions

We present 1D, 2D, and 3D hydrodynamical simulations of core-collapse supernovae including a parameterized neutrino heating and cooling scheme in order to investigate the critical core neutrino luminosity (L_crit) required for explosion. In contrast to some previous works, we find that 3D simulations explode later than 2D simulations, and that L_crit at fixed mass accretion rate is somewhat higher in 3D than in 2D. We find, however, that in 2D L_crit increases as the numerical resolution of the simulation increases. In contrast to some previous works, we argue that the average entropy of the gain region is in fact not a good indicator of explosion but is rather a reflection of the greater mass in the gain region in 2D. We compare our simulations to semi-analytic explosion criteria and examine the nature of the convective motions in 2D and 3D. We discuss the balance between neutrino-driven-buoyancy and drag forces. In particular, we show that the drag force will be proportional to a buoyant plume&#39;s surface area while the buoyant force is proportional to a plume&#39;s volume and, therefore, plumes with greater volume-to-surface area ratios will rise more quickly. We show that buoyant plumes in 2D are inherently larger, with greater volume-to-surface area ratios, than plumes in 3D. In the scenario that the supernova shock expansion is dominated by neutrino-driven buoyancy, this balance between buoyancy and drag forces may explain why 3D simulations explode later than 2D simulations and why L_crit increases with resolution. Finally, we provide a comparison of our results with other calculations in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Dependence of the Neutrino Mechanism of Core-Collapse Supernovae on the Equation of State

We study the dependence of the delayed neutrino-heating mechanism for core-collapse supernovae on the equation of state. Using a simplified treatment of the neutrino physics with a parameterized neutrino luminosity, we explore the relationship between explosion time, mass accretion rate, and neutrino luminosity for a 15 Msun progenitor in 1D and 2D. We test three different equations of state commonly used in core-collapse simulations: the models of Lattimer & Swesty (1991) with incompressibility of 180 MeV and 220 MeV, and the model of Shen et al. (1998), in order of increasing stiffness. We find that for a given neutrino luminosity the time after bounce until explosion increases with the stiffness of the equation of state: the Lattimer & Swesty EOS explode more easily than that of Shen et al. We find this holds in both 1D and 2D, while for all models explosions are obtained more easily in 2D than in 1D. We also discuss the relevance of approximate instability criteria to realistic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Identification of Faint Chandra X-ray Sources in the Core-Collapsed Globular Cluster NGC 6397: Evidence for a Bimodal Cataclysmic Variable Population

We have searched for optical identifications for 79 Chandra X-ray sources that lie within the half-mass radius of the nearby, core-collapsed globular cluster NGC 6397, using deep Hubble Space Telescope Advanced Camera for Surveys Wide Field Channel imaging in H-alpha, R, and B. Photometry of these images allows us to classify candidate counterparts based on color-magnitude diagram location. In addition to recovering nine previously detected cataclysmic variables (CVs), we have identified six additional faint CV candidates, a total of 42 active binaries (ABs), two millisecond pulsars (MSPs), one candidate active galactic nucleus, and one candidate interacting galaxy pair. Of the 79 sources, 69 have a plausible optical counterpart. The 15 likely and possible CVs in NGC 6397 mostly fall into two groups: a brighter group of six for which the optical emission is dominated by contributions from the secondary and accretion disk, and a fainter group of seven for which the white dwarf dominates the optical emission. There are two possible transitional objects that lie between these groups. The faintest CVs likely lie near the minimum of the CV period distribution, where an accumulation is expected. The spatial distribution of the brighter CVs is much more centrally concentrated than those of the fainter CVs and the active binaries. This may represent the result of an evolutionary process in which CVs are produced by dynamical interactions, such as exchange reactions, near the cluster center and are scattered to larger orbital radii, over their lifetimes, as they age and become fainter.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Accretion onto Seed Black Holes in the First Galaxies

The validity of the hypothesis that the massive black holes in high redshift quasars grew from stellar-sized &#34;seeds&#34; is contingent on a seed&#39;s ability to double its mass every few ten million years. This requires that the seed accrete at approximately the Eddington-limited rate. In the specific case of radiatively efficient quasiradial accretion in a metal-poor protogalactic medium, for which the Bondi accretion rate is often prescribed in cosmological simulations of massive black hole formation, we examine the effects of the radiation emitted near the black hole&#39;s event horizon on the structure of the surrounding gas flow. We find that the radiation pressure from photoionization significantly reduces the steady-state accretion rate and renders the quasiradial accretion flow unsteady and inefficient. The time-averaged accretion rates are a small fraction of the Eddington-limited accretion rate for Thomson scattering. The pressure of Ly-alpha photons trapped near the HII region surrounding the black hole may further attenuate the inflow. These results suggest that an alternative to quasiradial, radiatively efficient Bondi-like accretion should be sought to explain the rapid growth of quasar-progenitor seed black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Aspherical Core-Collapse Supernovae in Red Supergiants Powered by Nonrelativistic Jets

We explore the observational characteristics of jet-driven supernovae by simulating bipolar-jet-driven explosions in a red supergiant progenitor. We present results of four models in which we hold the injected kinetic energy at a constant $10^{51}$ ergs across all jet models but vary the specific characteristics of the jets to explore the influence of the nature of jets on the structure of the supernova ejecta. We evolve the explosions past shock-breakout and into quasi-homologous expansion of the supernova envelope into a red supergiant wind. The oppositely-directed, nickel-rich jets give a large-scale asymmetry that may account for the non-spherical excitation and substructure of spectral lines such as H$α$ and He I 10830Å. Jets with a large fraction of kinetic to thermal energy punch through the progenitor envelope and give rise to explosions that would be observed to be asymmetric from the earliest epochs, inconsistent with spectropolarimetric measurements of Type II supernovae. Jets with higher thermal energy fractions result in explosions that are roughly spherical at large radii but are significantly elongated at smaller radii, deep inside the ejecta, in agreement with the polarimetric observations. We present shock breakout light curves that indicate that strongly aspherical shock breakouts are incompatible with recent {\it GALEX} observations of shock breakout from red supergiant stars. Comparison with observations indicates that jets must deposit their kinetic energy efficiently throughout the ejecta while in the hydrogen envelope. Thermal energy-dominated jets satisfy this criterion and yield many of the observational characteristics of Type II supernovae.

0 citations0 reviewsNo code linkedOpen access