BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jeremiah W. Murphy

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.HE

Catalog footprint

What is connected

11works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

11 published item(s)

preprint2021arXiv

Gaia EDR3 confirms that Westerlund 1 is closer and older than previously thought

Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previously assumed. Compared to DR2, the EDR3 individual parallax uncertainties for each star decreased by 30\%. However, the aggregate parallax uncertainty for the cluster remained the same. This suggests that the uncertainty is dominated by systematics, which is possibly due to crowding, motions within the cluster, or motions due to binary orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring the parallax of Westerlund 1 from Gaia DR2

Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using {\it Gaia} Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The {\it Gaia} DR2 parallaxes reduce the uncertainty from a factor of 3 to 18\% and rules out the most often quoted value of 5 kpc with 99\% confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M$_{\odot}$; the new {\it Gaia} DR2 distance shifts this down to about 22 M$_{\odot}$. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Convection-Aided Explosions in One-Dimensional Core-Collapse Supernova Simulations I: Technique and Validation

Most one-dimensional core-collapse simulations fail to explode, yet multi-dimensional simulations often explode. A dominant multi-dimensional effect aiding explosion is neutrino-driven convection. We incorporate a convection model in approximate one-dimensional core-collapse supernova (CCSN) simulations. This is the 1D+ method. This convection model lowers the neutrino luminosity required for explosion by 30%, similar to the reduction observed in multi-dimensional simulations. The model is based upon the global turbulence model of Mabanta & Murphy (2018) and models the mean-field turbulent flow of neutrino-driven convection. In this preliminary investigation, we use simple neutrino heating and cooling algorithms to compare the critical condition in the 1D+ simulations with the critical condition observed in two-dimensional simulations. Qualitatively, the critical conditions in the 1D+ and the two-dimensional simulations are similar. The assumptions in the convection model affect the radial profiles of density, entropy, and temperature, and comparisons with the profiles of three dimensional simulations will help to calibrate these assumptions. These 1D+ simulations are consistent with the profiles and explosion conditions of equivalent two-dimensional CCSN simulations but are ~100 times faster, and the 1D+ prescription has the potential to be ~100,000 faster than three-dimensional CCSN simulations. The 1D+ technique will be ideally suited to test the explodability of thousands of progenitor models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Supernova Progenitor Mass Distributions of M31 and M33: Further Evidence for an Upper Mass Limit

Using Hubble Space Telescope (HST) photometry to measure star formation histories, we age-date the stellar populations surrounding supernova remnants (SNRs) in M31 and M33. We then apply stellar evolution models to the ages to infer the corresponding masses for their supernova progenitor stars. We analyze 33 M33 SNR progenitors and 29 M31 SNR progenitors in this work. We then combine these measurements with 53 previously published M31 SNR progenitor measurements to bring our total number of progenitor mass estimates to 115. To quantify the mass distributions, we fit power laws of the form $dN/dM \propto M^{-α}$. Our new, larger sample of M31 progenitors follows a distribution with $α= 4.4\pm 0.4$, and the M33 sample follows a distribution with $α= 3.8^{+0.4}_{-0.5}$. Thus both samples are consistent within the uncertainties, and the full sample across both galaxies gives $α= 4.2\pm 0.3$. Both the individual and full distributions display a paucity of massive stars when compared to a Salpeter initial mass function (IMF), which we would expect to observe if all massive stars exploded as SN that leave behind observable SNR. If we instead fix $α= 2.35$ and treat the maximum mass as a free parameter, we find $M_{max} \sim 35-45M_{sun}$, indicative of a potential maximum cutoff mass for SN production. Our results suggest that either SNR surveys are biased against finding objects in the youngest (<10 Myr old) regions, or the highest mass stars do not produce SNe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Dominance of Neutrino-Driven Convection in Core-Collapse Supernovae

Multi-dimensional instabilities have become an important ingredient in core-collapse supernova (CCSN) theory. Therefore, it is necessary to understand the driving mechanism of the dominant instability. We compare our parameterized three-dimensional CCSN simulations with other buoyancy-driven simulations and propose scaling relations for neutrino-driven convection. Through these comparisons, we infer that buoyancy-driven convection dominates post-shock turbulence in our simulations. In support of this inference, we present four major results. First, the convective fluxes and kinetic energies in the neutrino-heated region are consistent with expectations of buoyancy-driven convection. Second, the convective flux is positive where buoyancy actively drives convection, and the radial and tangential components of the kinetic energy are in rough equipartition (i.e. K_r ~ K_θ + K_ϕ). Both results are natural consequences of buoyancy-driven convection, and are commonly observed in simulations of convection. Third, buoyant driving is balanced by turbulent dissipation. Fourth, the convective luminosity and turbulent dissipation scale with the driving neutrino power. In all, these four results suggest that in neutrino-driven explosions, the multi-dimensional motions are consistent with neutrino-driven convection.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Investigation into the Character of Pre-Explosion Core-Collapse Supernova Shock Motion

We investigate the structure of the stalled supernova shock in both 2D and 3D and explore the differences in the effects of neutrino heating and the standing accretion shock instability (SASI). We find that early on the amplitude of the dipolar mode of the shock is factors of 2 to 3 smaller in 3D than in 2D. However, later in both 3D and 2D the monopole and dipole modes start to grow until explosion. Whereas in 2D the (l,m) = (1,0) mode changes sign quasi-periodically, producing the "up-and-down" motion always seen in modern 2D simulations, in 3D this almost never happens. Rather, in 3D when the dipolar mode starts to grow, it grows in magnitude and wanders stochastically in direction until settling before explosion to a particular patch of solid angle. In 2D we find that the amplitude of the dipolar shock deformation separates into two classes. For the first, identified with the SASI and for a wide range of "low" neutrino luminosities, this amplitude remains small and roughly constant. For the other, identified with higher luminosities and neutrino-driven convection, the dipolar amplitude grows sharply. Importantly, it is only for this higher luminosity class that we see neutrino-driven explosions within ~1 second of bounce. Moreover, for the "low" luminosity runs, the power spectra of these dipolar oscillations peak in the 30-50 Hz range associated with advection timescales, while for the high-luminosity runs the power spectra at lower frequencies are significantly more prominent. We associate this enhanced power at lower frequencies with slower convective effects and the secular growth of the dipolar shock amplitude. On the basis of our study, we hypothesize that neutrino-driven buoyant convection should almost always dominate the SASI when the supernova explosion is neutrino-driven.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dimensional Dependence of the Hydrodynamics of Core-Collapse Supernovae

The multidimensional character of the hydrodynamics in core-collapse supernova (CCSN) cores is a key facilitator of explosions. Unfortunately, much of this work has necessarily been performed assuming axisymmetry and it remains unclear whether or not this compromises those results. In this work, we present analyses of simplified two- and three-dimensional CCSN models with the goal of comparing the multidimensional hydrodynamics in setups that differ only in dimension. Not surprisingly, we find many differences between 2D and 3D models. While some differences are subtle and perhaps not crucial to understanding the explosion mechanism, others are quite dramatic and make interpreting 2D CCSN models problematic. In particular, we find that imposing axisymmetry artificially produces excess power at the largest spatial scales, power that has been deemed critical in the success of previous explosion models and has been attributed solely to the standing accretion shock instability. Nevertheless, our 3D models, which have an order of magnitude less power on large scales compared to 2D models, explode earlier. Since we see explosions earlier in 3D than in 2D, the vigorous sloshing associated with the large scale power in 2D models is either not critical in any dimension or the explosion mechanism operates differently in 2D and 3D. Possibly related to the earlier explosions in 3D, we find that about 25% of the accreted material spends more time in the gain region in 3D than in 2D, being exposed to more integrated heating and reaching higher peak entropies, an effect we associate with the differing characters of turbulence in 2D and 3D. Finally, we discuss a simple model for the runaway growth of buoyant bubbles that is able to quantitatively account for the growth of the shock radius and predicts a critical luminosity relation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Supernova Remnant Progenitor Masses in M31

Using HST photometry, we age-date 59 supernova remnants (SNRs) in the spiral galaxy M31 and use these ages to estimate zero-age main sequence masses (MZAMS) for their progenitors. To accomplish this, we create color-magnitude diagrams (CMDs) and use CMD fitting to measure the recent star formation history (SFH) of the regions surrounding cataloged SNR sites. We identify any young coeval population that likely produced the progenitor star and assign an age and uncertainty to that population. Application of stellar evolution models allows us to infer the MZAMS from this age. Because our technique is not contingent on precise location of the progenitor star, it can be applied to the location of any known SNR. We identify significant young SF around 53 of the 59 SNRs and assign progenitor masses to these, representing a factor of 2 increase over currently measured progenitor masses. We consider the remaining 6 SNRs as either probable Type Ia candidates or the result of core-collapse progenitors that have escaped their birth sites. The distribution of recovered progenitor masses is bottom heavy, showing a paucity of the most massive stars. If we assume a single power law distribution, dN/dM proportional to M^alpha, we find a distribution that is steeper than a Salpeter IMF (alpha=-2.35). In particular, we find values of alpha outside the range -2.7 to -4.4 inconsistent with our measured distribution at 95% confidence. If instead we assume a distribution that follows a Salpeter IMF up to some maximum mass, we find that values of M_max greater than 26 Msun are inconsistent with the measured distribution at 95% confidence. In either scenario, the data suggest that some fraction of massive stars may not explode. The result is preliminary and requires more SNRs and further analysis. In addition, we use our distribution to estimate a minimum mass for core collapse between 7.0 and 7.8 Msun.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Global Turbulence Model for Neutrino-Driven Convection in Core-Collapse Supernovae

Simulations of core-collapse supernovae (CCSNe) result in successful explosions once the neutrino luminosity exceeds a critical curve, and recent simulations indicate that turbulence further enables explosion by reducing this critical neutrino luminosity. We propose a theoretical framework to derive this result and take the first steps by deriving the governing mean-field equations. Using Reynolds decomposition, we decompose flow variables into background and turbulent flows and derive self-consistent averaged equations for their evolution. As basic requirements for the CCSN problem, these equations naturally incorporate steady-state accretion, neutrino heating and cooling, non-zero entropy gradients, and turbulence terms associated with buoyant driving, redistribution, and dissipation. Furthermore, analysis of two-dimensional (2D) CCSN simulations validate these Reynolds-averaged equations, and we show that the physics of turbulence entirely accounts for the differences between 1D and 2D CCSN simulations. As a prelude to deriving the reduction in the critical luminosity, we identify the turbulent terms that most influence the conditions for explosion. Generically, turbulence equations require closure models, but these closure models depend upon the macroscopic properties of the flow. To derive a closure model that is appropriate for CCSNe, we cull the literature for relevant closure models and compare each with 2D simulations. These models employ local closure approximations and fail to reproduce the global properties of neutrino-driven turbulence. Motivated by the generic failure of these local models, we propose an original model for turbulence which incorporates global properties of the flow. This global model accurately reproduces the turbulence profiles and evolution of 2D CCSN simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Progenitor Mass of SN 2011dh from Stellar Populations Analysis

Using Hubble Space Telescope (HST) photometry, we characterize the age of the stellar association in the vicinity of supernova (SN) 2011dh and use it to infer the zero-age main sequence mass (M_{ZAMS}) of the progenitor star. We find two distinct and significant star formation events with ages of <6 and 17^{+3}_{-4} Myrs, and the corresponding M_{ZAMS} are >29 and 13^{+2}_{-1} M_{Sun}, respectively. These two bursts represent 18^{+4}_{-9}% (young) and 64^{+10}_{-14}% (old) of the total star formation in the last 50 Myrs. Adopting these fractions as probabilities suggests that the most probable M_{ZAMS} is 13^{+2}_{-1} M_{Sun}. These results are most sensitive to the luminosity function along the well-understood main sequence and are less sensitive to uncertain late-stage stellar evolution. Therefore, they stand even if the progenitor suffered disruptive post-main-sequence evolution (e.g. eruptive mass loss or binary Roche-lobe overflow). Progenitor identification will help to further constrain the appropriate population. Even though pre-explosion images show a yellow supergiant (YSG) at the site of the SN, panchromatic SN light curves suggest a more compact star as the progenitor. In spite of this, our results suggest an association between the YSG and the SN. Not only was the star located at the SN site, but reinforcing an association, the star's bolometric luminosity is consistent with the final evolutionary stage of the 17 Myr old star burst. If the YSG disappears, then M_{ZAMS}=13^{+2}_{-1} M_{Sun}, but if it persists, then our results allow the possibility that the progenitor was an unseen star of >29 M_{Sun}.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Model for Gravitational Wave Emission from Neutrino-Driven Core-Collapse Supernovae

Abridged: Using a suite of progenitor models, neutrino luminosities, and two- dimensional (2D) simulations, we investigate the matter gravitational-wave (GW) emission from postbounce phases of neutrino-driven core-collapse supernovae (CCSNe). The relevant phases are prompt and steady-state convection, the standing accretion shock instability (SASI), and asymmetric explosions. For the stages before explosion, we propose a model for the source of GW emission. Downdrafts of the postshock-convection/SASI region strike the protoneutron star "surface" with large speeds and are decelerated by buoyancy forces. We find that the GW amplitude is set by the magnitude of deceleration and, by extension, the downdraft's speed and the vigor of postshock-convective/SASI motions. However, the characteristic frequencies, which evolve from ~100 Hz after bounce to ~300-400 Hz, are primarily independent of these speeds, but are set by the deceleration timescale, which is in turn set by the buoyancy frequency at the lower boundary of postshock convection. Consequently, the GW characteristic frequencies are dependent upon a combination of core structure attributes, specifically the dense-matter equation of state (EOS) and details that determine the gradients at the boundary, including the accretion-rate history, the EOS at subnuclear densities, and neutrino transport. During explosion, the high frequency signal wanes and is replaced by a strong low frequency, ~10s of Hz, signal that reveals the general morphology of the explosion (i.e. prolate, oblate, or spherical). However, current and near-future GW detectors are sensitive to GW power at frequencies >50 Hz. Therefore, the signature of explosion will be the abrupt reduction of detectable GW emission.

0 citations0 reviewsNo code linkedOpen access