Next-Generation Simulations of The Remarkable Deaths of Massive Stars [Slides]
Abstract not provided.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
Not Available
Globular clusters are of paramount importance for testing theories of stellar evolution and early galaxy formation. Strong evidence for multiple populations of stars in globular clusters derives from observed abundance anomalies. A puzzling example is the recently detected Mg–K anticorrelation in NGC 2419. We perform Monte Carlo nuclear reaction network calculations to constrain the temperature–density conditions that gave rise to the elemental abundances observed in this elusive cluster. We find a correlation between stellar temperature and density values that provide a satisfactory match between simulated and observed abundances in NGC 2419 for all relevant elements (Mg, Si, K, Ca, Sc, Ti, and V). Except at the highest densities (ρ ≳ 10{sup 8} g cm{sup −3}), the acceptable conditions range from ≈100 MK at ≈10{sup 8} g cm{sup −3} to ≈200 MK at ≈10{sup −4} g cm{sup −3}. This result accounts for uncertainties in nuclear reaction rates and variations in the assumed initial composition. We review hydrogen-burning sites and find that low-mass stars, asymptotic giant branch (AGB) stars, massive stars, or supermassive stars cannot account for the observed abundance anomalies in NGC 2419. Super-AGB stars could be viable candidates for the polluter stars if stellar model parameters can be fine-tuned to produce higher temperatures. Novae, involving either CO or ONe white dwarfs, could be interesting polluter candidates, but a current lack of low-metallicity nova models precludes firmer conclusions. We also discuss whether additional constraints for the first-generation polluters can be obtained by future measurements of oxygen, or by evolving models of second-generation low-mass stars with a non-canonical initial composition.
The supernova impostor PSN J09132750+7627410 in NGC 2748 reached a maximum luminosity of ≈−14 mag. It was quickly realized that it was not a true supernova, but another example of a nonterminal giant eruption. PSN J09132750+7627410 is distinguished by multiple P Cygni absorption minima in the Balmer emission lines that correspond to outflow velocities of −400, −1100, and −1600 km s{sup −1}. Multiple outflows have been observed in only a few other objects. In this paper we describe the evolution of the spectrum and the P Cygni profiles for 3 months past maximum, the post-maximum formation of a cool, dense wind, and the identification of a possible progenitor. One of the possible progenitors is an infrared source. Its pre-eruption spectral energy distribution suggests a bolometric luminosity of −8.3 mag and a dust temperature of 780 K. If it is the progenitor, it is above the AGB limit, unlike the intermediate-luminosity red transients. The three P Cygni profiles could be due to ejecta from the current eruption, the wind of the progenitor, or previous mass-loss events. We suggest that they were all formed as part of the same high-mass-loss event and are due to material ejected at different velocities or energies. We also suggest that multiple outflows during giant eruptions may be more common than reported.
In modern cosmology, determining the Hubble constant (H{sub 0}) using a distance ladder to percent level and comparing with the results from the Planck satellite can shed light on the nature of dark energy, physics of the neutrino, and curvature of the universe. Thanks to the endeavor of the SH0ES team, the uncertainty of the H{sub 0} has be dramatically reduced, from 10% to 2.4%, and with the promise of even reaching 1% in the near future. In this regard, it is fundamentally important to investigate the systematics. This is best done using other good independent distance indicators. One promising method is the flux-weighted gravity luminosity relation (FGLR) of the blue supergiants (BSGs). As BSGs are the brightest objects in galaxies, they can probe distances up to 10 Mpc with negligible blending effects. While the FGLR method delivered distance is in good agreement with other distance indicators, it has been shown that this method delivers greater distances in the cases of M33 and NGC 55. Here, we investigate whether the M33 distance estimate of FGLR suffers systematics from stellar variability. Using CFHT M33 monitoring data, we found that 9 out of 22 BSGs showed variability during the course of 500 days, although with amplitudes as small as 0.05 mag. This suggests that stellar variability plays a negligible role in the FGLR distance determination.
In this paper, we describe a program that identifies in the optical spectrum the main parameters of a spectral line, namely the initial and final wavelengths, and the line depth. Moreover, using numerical calculations, it identifies and removes adjacent lines. Next, the program calculates the equivalent width and the FWHM. The software was tested in a sample of 300 lines in two databases of synthetic spectra generated by the CMFGEN and PoWR codes, and 300 lines in observed spectra from the IACOB database, showing a Gaussian distribution of relative errors, from which it is inferred that 80% of the measured lines have errors less than 17% and only 5% of the lines have errors greater than 26%. The program was also run on the entire database of 45,000 CMFGEN and 202 POWR synthetic spectra, generating a library of H i, He i, and He ii lines necessary to feed the FITspec code for the derivation of stellar parameters: effective temperature, surface gravity, and luminosity.
We present the first 4 π– three-dimensional (3D) simulation of the last minutes of oxygen shell burning in an 18 M {sub ⊙} supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a one-dimensional stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ∼0.1 at collapse, and an ℓ = 2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 M {sub ⊙} to 0.56 M {sub ⊙} due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing-length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12% – 24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.
We present the results of a Hubble Space Telescope ACS/HRC FUV, ACS/WFC optical study into the cluster populations of a sample of 22 Luminous Infrared Galaxies in the Great Observatories All-Sky LIRG Survey. Through integrated broadband photometry, we have derived ages and masses for a total of 484 star clusters contained within these systems. This allows us to examine the properties of star clusters found in the extreme environments of LIRGs relative to lower luminosity star-forming galaxies in the local universe. We find that by adopting a Bruzual and Charlot simple stellar population model and Salpeter initial mass function, the age distribution of the clusters declines as dN/dτ=τ{sup −0.9+/−0.3}, consistent with the age distribution derived for the Antennae Galaxies, and interpreted as evidence for rapid cluster disruption occurring in the strong tidal fields of merging galaxies. The large number of 10{sup 6} M{sub ⊙} young clusters identified in the sample also suggests that LIRGs are capable of producing more high-mass clusters than what is observed to date in any lower luminosity star-forming galaxy in the local universe. The observed cluster mass distribution of dN/dM=M{sup −1.95+/−0.11} is consistent with the canonical −2 power law used to describe the underlying initial cluster mass function (ICMF) for a wide range of galactic environments. We interpret this as evidence against mass-dependent cluster disruption, which would flatten the observed CMF relative to the underlying ICMF distribution.
We present an abundance analysis of seven super star clusters in the disk of M83. The near-infrared spectra of these clusters are dominated by red supergiants, and the spectral similarity in the J-band of such stars at uniform metallicity means that the integrated light from the clusters may be analyzed using the same tools as those applied to single stars. Using data from VLT/KMOS, we estimate metallicities for each cluster in the sample. We find that the abundance gradient in the inner regions of M83 is flat, with a central metallicity of [Z]=0.21±0.11 relative to a solar value of Z {sub ⊙} = 0.014, which is in excellent agreement with the results from an analysis of luminous hot stars in the same regions. Compiling this latest study with our other recent work, we construct a mass–metallicity relation for nearby galaxies based entirely on the analysis of RSGs. We find excellent agreement with the other stellar-based technique—that of blue supergiants—as well as with temperature-sensitive (“auroral” or “direct”) H ii-region studies. Of all the H ii-region strong-line calibrations, those that are empirically calibrated to direct-method studies (N2 and O3N2) provide the most consistent results.
This paper is a sequel to our 2015 paper, Kato et al., which calculated the luminosities and spectra of electron-type anti-neutrinos ( ν-bar {sub e}) from the progenitors of core-collapse supernovae. Expecting that the capability to detect electron-type neutrinos (ν{sub e}) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to include all flavors of neutrinos emitted from the pre-bounce phase. We pick up three progenitor models of electron capture supernovae (ECSNe) and iron-core collapse supernovae (FeCCSNe). We find that the number luminosities reach ∼10{sup 57} s{sup –1} and ∼10{sup 53} s{sup –1 }at maximum for ν{sub e} and ν-bar {sub e}, respectively. We also estimate the numbers of detection events at terrestrial neutrino detectors including DUNE, taking flavor oscillations into account and assuming the distance to the progenitors to be 200 pc. It is demonstrated that ν-bar {sub e} from the ECSN progenitor will be undetected at almost all detectors, whereas we will be able to observe ≳15,900 ν{sub e} at DUNE for the inverted mass hierarchy. From the FeCCSN progenitors, the number of ν-bar {sub e} events will be largest for JUNO, 200–900 ν-bar {sub e}, depending on the mass hierarchy, whereas the number of ν{sub e} events at DUNE is ≳2100 for the inverted mass hierarchy. These results imply that the detection of ν-bar {sub e} is useful to distinguish progenitors of FeCCSNe from those of ECSNe, while ν{sub e} will provide us with detailed information on the collapse phase regardless of the type and mass of the progenitor.
We present a new method for neutrino-matter coupling in multi-dimensional radiation-hydrodynamic simulations of core-collapse supernovae (CCSNe) with the full Boltzmann neutrino transport. This development is motivated by the fact that accurate conservation of momentum is required for reliable numerical modelings of CCSN dynamics including a recoil of proto-neutron stars (PNSs). The new method is built on a hybrid approach in which we use the energy-momentum tensor of neutrinos to compute the momentum feedback from neutrino to matter in the optically thick region while we employ the collision term in the optically thin region. In this method we utilize a general relativistic description of radiation-hydrodynamics with angular moments, which allows us to evaluate the momentum feedback from neutrino to matter without inconsistency with our Boltzmann solver. We demonstrate that the new method substantially improves the accuracy of linear momentum conservation in our CCSN simulations under reasonable angular resolutions in momentum space, alleviating the difficulty in giving the diffusion limit precisely with the discrete ordinate ($S_n$) method. It is the first ever demonstration that the PNS kick can be handled directly and properly in multi-dimensional radiation-hydrodynamic simulations with the full Boltzmann neutrino transport.
Abstract We study the effects of additional cooling due to the emission of a dark matter candidate particle, the dark photon, on the final phases of the evolution of a 15 M ⊙ star and resulting modifications of the pre-supernova neutrino signal. For a substantial portion of the dark photon parameter space the extra cooling speeds up Si burning, which results in a reduced number of neutrinos emitted during the last day before core collapse. This reduction can be described by a systematic acceleration of the relevant timescales and the results can be estimated semi-analytically in good agreement with the numerical simulations. Outside the semi-analytic regime we find more complicated effects. In a narrow parameter range, low-mass dark photons lead to an increase in the number of emitted neutrinos because of additional shell-burning episodes that delay core collapse. Furthermore, relatively strong couplings produce a thermonuclear runaway during O burning, which could result in a complete disruption of the star but requires more detailed simulations to determine the outcome. Our results show that pre-supernova neutrino signals are a potential probe of the dark photon parameter space.
We present 3D full star simulations, reaching up to 90% of the total stellar radius, for three 7 M⊙ stars of different ages: zero-age main sequence (ZAMS), mid–main sequence (midMS), and terminal-age main sequence (TAMS). A comparison with several theoretical prescriptions shows that the generation spectra for all three ages are dominated by convective plumes. Two distinct overshooting layers are observed, with most plumes stopped within the layer situated directly above the convective boundary; overshooting to the second, deeper layer becomes progressively more infrequent with increasing stellar age. Internal gravity wave (IGW) propagation is significantly impacted in the midMS and TAMS models as a result of some IGWs getting trapped within their Brunt–Väisälä frequency spikes. A fundamental change in the wave structure across radius is also observed, driven by the effect of density stratification on IGW propagation causing waves to become evanescent within the radiative zone, with older stars being affected more strongly. We find that the steepness of the frequency spectrum at the surface increases from ZAMS to the older models, with older stars also showing more modes in their spectra.
We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.