Engineering topics

Pan, Kuo-Chuan

Publications and source records attributed to Pan, Kuo-Chuan.

13 records

DOE OSTI2022

Evolution of Main-sequence-like Surviving Companions in Type Ia Supernova Remnants

Recent theoretical and numerical studies of Type Ia supernovae (SNe Ia) explosions within the single-degenerate scenario suggest that the nondegenerate companions could survive during the supernova impact and could be detectable in nearby supernova remnants. However, observational efforts show less promising evidence of the existence of surviving companions from the standard single-degenerate channels. The spin-up/spin-down models are possible mechanisms to explain the nondetection of surviving companions. In these models, the spin-up phase could increase the critical mass for explosion, leading to a super-Chandrasekhar-mass explosion, and the spin-down phase could lead to extra mass loss and angular momentum redistribution. Since the spin-down timescale for the delayed explosion of a rotating white dwarf is unclear, in this paper we explore a vast parameter space of main-sequence-like surviving companions via two-dimensional hydrodynamic simulations of supernova impact and the subsequent stellar evolution of surviving companions. Tight universal relations to describe the mass-stripping effect, supernova kick, and depth of supernova heating are provided. Our results suggest that the not-yet-detected surviving companions from observations of nearby SN Ia remnants might favor low-mass companions, short binary separation, or stronger supernova explosion energies than the standard single-degenerate channels.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Stellar Mass Black Hole Formation and Multimessenger Signals from Three-dimensional Rotating Core-collapse Supernova Simulations

We present self-consistent 3D core-collapse supernova simulations of a 40 ${M}_{\odot }$ progenitor model using the isotropic diffusion source approximation for neutrino transport and an effective general relativistic potential up to `0.9 s postbounce. We consider three different rotational speeds with initial angular velocities of ${{\rm{\Omega }}}_{0}=0$, 0.5, and 1 rad s -1 and investigate the impact of rotation on shock dynamics, black hole (BH) formation, and gravitational wave (GW) signals. The rapidly rotating model undergoes an early explosion at ~250 ms postbounce and shows signs of the low $T/| W| $ instability. We do not find BH formation in this model within ~460 ms postbounce. In contrast, we find BH formation at 776 ms postbounce and 936 ms postbounce for the nonrotating and slowly rotating models, respectively. The slowly rotating model explodes at ~650 ms postbounce, and the subsequent fallback accretion onto the proto–neutron star (PNS) results in BH formation. In addition, the standing accretion shock instability induces rotation of the PNS in the model that started with a nonrotating progenitor. Assuming conservation of specific angular momentum during BH formation, this corresponds to a BH spin parameter of a = J/M = 0.046. However, if no explosion sets in, all the angular momentum will eventually be accreted by the BH, resulting in a nonspinning BH. The successful explosion of the slowly rotating model drastically slows down the accretion onto the PNS, allowing continued cooling and contraction that results in an extremely high GW frequency (f ~ 3000 Hz) at BH formation, while the nonrotating model generates GW signals similar to our corresponding 2D simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Determining the Structure of Rotating Massive Stellar Cores with Gravitational Waves

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 neutrino-radiation-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.

Astronomy & Astrophysics↗

DOE OSTI2019

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

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 ~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 ⊙ to 60 M ⊙ 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. Finally, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

The impact of different neutrino transport methods on multidimensional core-collapse supernova simulations

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 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 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 ~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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Equation of State Dependent Dynamics and Multi-messenger Signals from Stellar-mass Black Hole Formation

We investigate axisymmetric black hole (BH) formation and its gravitational wave (GW) and neutrino signals with self-consistent core-collapse supernova simulations of a non-rotating 40 M ⊙ progenitor star using the isotropic diffusion source approximation for the neutrino transport and a modified gravitational potential for general relativistic effects. We consider four different neutron star (NS) equations of state (EoS): LS220, SFHo, BHBΛ$\phi$, and DD2, and study the impact of the EoS on BH formation dynamics and GW emission. We find that the BH formation time is sensitive to the EoS from 460 to >1300 ms and is delayed in multiple dimensions for ~100–250 ms due to the finite entropy effects. Depending on the EoS, our simulations show the possibility that shock revival can occur along with the collapse of the proto-neutron star (PNS) to a BH. The gravitational waveforms contain four major features that are similar to previous studies but show extreme values: (1) a low-frequency signal (~300–500 Hz) from core-bounce and prompt convection, (2) a strong signal from the PNS g-mode oscillation among other features, (3) a high-frequency signal from the PNS inner-core convection, and (4) signals from the standing accretion shock instability and convection. The peak frequency at the onset of BH formation reaches to ~2.3 kHz. Thus, the characteristic amplitude of a 10 kpc object at peak frequency is detectable but close to the noise threshold of the Advanced LIGO and KAGRA, suggesting that the next-generation GW detector will need to improve the sensitivity at the kHz domain to better observe stellar-mass BH formation from core-collapse supernovae or failed supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Nature of the Diffuse Source and Its Central Point-like Source in SNR 0509–67.5

We examine a diffuse emission region near the center of SNR 0509−67.5 to determine its nature. Within this diffuse region we observe a point-like source that is bright in the near-IR, but is not visible in the B and V bands. We consider an emission line observed at 6766 Å and the possibilities that it is Ly α , H α , and [O ii] λ 3727. We examine the spectral energy distribution (SED) of the source, comprised of Hubble Space Telescope B , V , I , J , and H bands in addition to Spitzer /IRAC 3.6, 4.5, 5.8, and 8 μ m bands. The peak of the SED is consistent with a background galaxy at z ≈ 0.8 ± 0.2 and a possible Balmer jump places the galaxy at z ≈ 0.9 ± 0.3. These SED considerations support the emission line’s identification as [O ii] λ 3727. We conclude that the diffuse source in SNR 0509−67.5 is a background galaxy at z ≈ 0.82. Furthermore, we identify the point-like source superposed near the center of the galaxy as its central bulge. Finally, we find no evidence for a surviving companion star, indicating a double-degenerate origin for SNR 0509−67.5.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

TWO-DIMENSIONAL CORE-COLLAPSE SUPERNOVA SIMULATIONS WITH THE ISOTROPIC DIFFUSION SOURCE APPROXIMATION FOR NEUTRINO TRANSPORT

The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation–hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped-particle and streaming-particle components. Heavy neutrinos are described by a leakage scheme. Unlike the “ray-by-ray” approach in some other multidimensional supernova models, we use cylindrical coordinates and solve the trapped-particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We provide an IDSA verification by performing one-dimensional (1D) and 2D simulations with 15 and 20 M{sub ⊙} progenitors from Woosley et al. and discuss the difference between our IDSA results and those existing in the literature. Additionally, we perform Newtonian 1D and 2D simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 M{sub ⊙} progenitors from Woosley et al. with the HS(DD2) equation of state. General-relativistic effects are neglected. We obtain robust explosions with diagnostic energies E{sub dia} ≳ 0.1–0.5 B (1 B ≡ 10{sup 51} erg) for all considered 2D models within approximately 100–300 ms after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affects the postbounce evolution, e.g., the neglect of neutrino–electron scattering during collapse will lead to a stronger explosion.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

SIMULATIONS OF THE SYMBIOTIC RECURRENT NOVA V407 CYG. I. ACCRETION AND SHOCK EVOLUTIONS

The shock interaction and evolution of nova ejecta with wind from a red giant (RG) star in a symbiotic binary system are investigated via three-dimensional hydrodynamics simulations. We specifically model the 2010 March outburst of the symbiotic recurrent nova V407 Cygni from its quiescent phase to its eruption phase. The circumstellar density enhancement due to wind–white-dwarf interaction is studied in detail. It is found that the density-enhancement efficiency depends on the ratio of the orbital speed to the RG wind speed. Unlike another recurrent nova, RS Ophiuchi, we do not observe a strong disk-like density enhancement, but instead observe an aspherical density distribution with ∼20% higher density in the equatorial plane than at the poles. To model the 2010 outburst, we consider several physical parameters, including the RG mass-loss rate, nova eruption energy, and ejecta mass. A detailed study of the shock interaction and evolution reveals that the interaction of shocks with the RG wind generates strong Rayleigh–Taylor instabilities. In addition, the presence of the companion and circumstellar density enhancement greatly alter the shock evolution during the nova phase. Depending on the model, the ejecta speed after sweeping out most of the circumstellar medium decreases to ∼100–300 km s{sup −1}, which is consistent with the observed extended redward emission in [N ii] lines in 2011 April.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Search for surviving companions in type Ia supernova remnants

The nature of the progenitor systems of type Ia supernovae (SNe Ia) is still unclear. One way to distinguish between the single-degenerate scenario and double-degenerate scenario for their progenitors is to search for the surviving companions (SCs). Using a technique that couples the results from multi-dimensional hydrodynamics simulations with calculations of the structure and evolution of main-sequence- (MS-) and helium-rich SCs, the color and magnitude of MS- and helium-rich SCs are predicted as functions of time. The SC candidates in Galactic type Ia supernova remnants (Ia SNR) and nearby extragalactic Ia SNRs are discussed. We find that the maximum detectable distance of MS SCs (helium-rich SCs) is 0.6-4 Mpc (0.4-16 Mpc), if the apparent magnitude limit is 27 in the absence of extinction, suggesting that the Large and Small Magellanic Clouds and the Andromeda Galaxy are excellent environments in which to search for SCs. However, only five Ia SNRs have been searched for SCs, showing little support for the standard channels in the singe-degenerate scenario. To better understand the progenitors of SNe Ia, we encourage the search for SCs in other nearby Ia SNRs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

EVOLUTION OF POST-IMPACT REMNANT HELIUM STARS IN TYPE Ia SUPERNOVA REMNANTS WITHIN THE SINGLE-DEGENERATE SCENARIO

The progenitor systems of Type Ia supernovae (SNe Ia) are still under debate. Based on recent hydrodynamics simulations, non-degenerate companions in the single-degenerate scenario (SDS) should survive the supernova (SN) impact. One way to distinguish between the SDS and the double-degenerate scenario is to search for the post-impact remnant stars (PIRSs) in SN Ia remnants. Using a technique that combines multi-dimensional hydrodynamics simulations with one-dimensional stellar evolution simulations, we have examined the post-impact evolution of helium-rich binary companions in the SDS. It is found that these helium-rich PIRSs (He PIRSs) dramatically expand and evolve to a luminous phase (L {approx} 10{sup 4} L{sub Sun }) about 10 yr after an SN explosion. Subsequently, they contract and evolve to become hot blue-subdwarf-like (sdO-like) stars by releasing gravitational energy, persisting as sdO-like stars for several million years before evolving to the helium red-giant phase. We therefore predict that a luminous OB-like star should be detectable within {approx}30 yr after the SN explosion. Thereafter, it will shrink and become an sdO-like star in the central regions of SN Ia remnants within star-forming regions for SN Ia progenitors evolved via the helium-star channel in the SDS. These He PIRSs are predicted to be rapidly rotating (v{sub rot} {approx}> 50 km s{sup -1}) and to have high spatial velocities (v{sub linear} {approx}> 500 km s{sup -1}). Furthermore, if SN remnants have diffused away and are not recognizable at a later stage, He PIRSs could be an additional source of single sdO stars and/or hypervelocity stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

EVOLUTION OF POST-IMPACT COMPANION STARS IN SN Ia REMNANTS WITHIN THE SINGLE-DEGENERATE SCENARIO

The nature of the progenitor systems of Type Ia supernovae is still uncertain. One way to distinguish between the single-degenerate scenario and double-degenerate scenario is to search for the post-impact remnant star. To examine the characteristics of the post-impact remnant star, we have carried out three-dimensional hydrodynamic simulations of supernova impacts on main-sequence-like stars. We explore the evolution of the post-impact remnants using the stellar evolution code MESA. We find that the luminosity and radius of the remnant star dramatically increase just after the impact. After the explosion, post-impact companions continue to expand on a progenitor-dependent timescale of {approx}10{sup 2.5}-10{sup 3} years before contracting. It is found that the time evolution of the remnant star is dependent not only on the amount of energy absorbed but also on the depth of the energy deposition. We examine the viability of the candidate star Tycho G as the possible remnant companion in Tycho's supernova by comparing it to the evolved post-impact remnant stars in our simulations. The closest model in our simulations has a similar effective temperature, but the luminosity and radius are twice as large. By examining the angular momentum distribution in our simulations, we find that the surface rotational speed could drop to {approx}10 km s{sup -1} if the specific angular momentum is conserved during the post-impact evolution, implying that Tycho G cannot be completely ruled out because of its low surface rotation speed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

Scalable Algorithms for Distributed-Memory Adaptive Mesh Refinement

This paper presents scalable algorithms and data structures for adaptive mesh refinement computations. We describe a novel mesh restructuring algorithm for adaptive mesh refinement computations that uses a constant number of collectives regardless of the refinement depth. To further increase scalability, we describe a localized hierarchical coordinate-based block indexing scheme in contrast to traditional linear numbering schemes, which incur unnecessary synchronization. In contrast to the existing approaches which take O(P) time and storage per process, our approach takes only constant time and has very small memory footprint. With these optimizations as well as an efficient mapping scheme, our algorithm is scalable and suitable for large, highly-refined meshes. We present strong-scaling experiments up to 2k ranks on Cray XK6, and 32k ranks on IBM Blue Gene/Q.

Langer, Akhil↗