Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gravitational collapse”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Occurrence of Gravitational Collapse in the Accreting Neutron Stars of Binary-driven Hypernovae

The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon–oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. This article aims to determine the conditions under which a black hole (BH) forms from NS collapse induced by the accretion and the impact on the GRB's observational properties and taxonomy. We perform three-dimensional, smoothed particle hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state, with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS in the binary to the critical mass for gravitational collapse into a BH and its occurrence time, t col . We include a nonzero angular momentum of the NSs and find that t col ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about 5 minute orbital period), promptly form a BH, and release ≳10 52 erg of energy. They form NS–BH binaries with tens of kiloyears merger timescales by gravitational-wave emission. BdHNe II and III do not form BHs, and release ~10 50 –10 52 erg and ≲10 50 erg of energy, respectively. They form NS–NS binaries with a range of merger timescales larger than for NS–BH binaries. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a nonrotating NS. Magnetic braking by a 10 13 G field can delay BH formation, leading to BH–BH or NS–BH with tens of kiloyears merger timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Unveiling horizons in quantum critical collapse

Critical gravitational collapse offers a unique window into regimes of arbitrarily high curvature, culminating in a naked singularity arising from smooth initial data — thus providing a dynamical counterexample to weak cosmic censorship. Near the critical regime, quantum effects from the collapsing matter are expected to intervene before full quantum gravity resolves the singularity. Despite its fundamental significance, a self-consistent treatment has so far remained elusive. In this work, we perform a one-loop semiclassical analysis using the robust anomaly-based method in the canonical setup of Einstein gravity minimally coupled to a free, massless scalar field. Focusing on explicitly solvable near-critical solutions in both 2 + 1 and 3 + 1 dimensions, we analytically solve the semiclassical Einstein equations and obtain controlled, quantitative results for several long-standing questions within the dominant s-wave sector. We find that regularity uniquely selects a Boulware-like quantum state, encoding genuine vacuum polarization effects from the collapsing matter. Remarkably, the resulting quantum corrections manifest as a growing mode. Horizon-tracing analyses, incorporating both classical and quantum modes, reveal the emergence of a finite mass gap, signaling a phase transition from classical Type II to quantum-modified Type I behavior, thereby providing a quantum enforcement of the weak cosmic censorship. The most non-trivial aspect of our analysis involves dealing with non-conformal matter fields in explicitly time-dependent critical spacetimes. Along the way, we uncover intriguing and previously underexplored features of quantum field theory in curved spacetime.

2D Gravity↗
DOE OSTI2023

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

How do axisymmetric black holes grow monopole and dipole hair?

Here, we study the dynamical formation of scalar monopole and dipole hair in scalar Gauss-Bonnet theory and dynamical Chern-Simons theory. We prove that the spherically symmetric mode of the dipole hair is completely determined by the product of the mass of the spacetime and the value of the monopole hair. We then show that the dynamics of the ℓ = 1 mode of the dipole hair is intimately tied to the appearance of the event horizon during axisymmetric collapse, which results in the radiation of certain modes that could have been divergent in the future of the collapse. We confirm these analytical predictions by simulating the gravitational collapse of a rapidly rotating neutron star in the decoupling limit, both in scalar Gauss-Bonnet and dynamical Chern-Simons theory. Our results, combined with those provide a clear physical picture of the dynamics of scalar monopole and dipole radiation in axisymmetric and spherical gravitational collapse in these theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Equation-of-state Dependence of Gravitational Waves in Core-collapse Supernovae

Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars 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 CCSN simulations, subsequently impacting the GW emission. In this study, 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*, which is crucial in setting the thermal dependence of the EOS, and (2) the isoscalar incompressibility modulus, K sat . While K 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* 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 ~1250 Hz; we investigate its origin and report its EOS dependence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Embedding condition describing shearing, dissipative collapse

Highlights: • Perturbative approach to the time-dependent Karmarkar condition in shearing collapse. • Analysis of the (in)stability of the collapsing fluid using the TOV equation. • Study of the temperature profiles using a causal heat transport equation. In this work we employ the time-dependent Karmarkar condition to model dissipative gravitational collapse in the presence of shear which leads to a final static configuration. We believe that this is a first attempt at using the time-dependent Karmarkar condition describing a radiating star with shear in which the collapsing core dissipates energy in the form of a radial heat flux. The final static configuration which is described by an embedding class I spacetime was obtained by Hansraj and Moodly who showed that this solution satisfies all the physical requirements for a static compact object. As the collapse proceeds and the heat flux tends to zero the dynamical core evolves into a static compact core with the exterior spacetime being described by the vacuum Schwarzschild solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Protoplanetary Disk Rings as Sites for Planetesimal Formation

Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic-jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI itself requires some initial concentration of dust. Here we present the first 3D simulations of planetesimal formation in the presence of a pressure bump modeled specifically after those seen by Atacama Large Millimeter/submillimeter Array. We place a pressure bump at the center of a large 3D shearing box, along with an initial solid-to-gas ratio of Z = 0.01, and we include both particle back-reaction and particle self-gravity. We consider millimeter-sized and centimeter-sized particles separately. For simulations with centimeter-sized particles, we find that even a small pressure bump leads to the formation of planetesimals via the SI; a pressure bump does not need to fully halt radial particle drift for the SI to become efficient. Furthermore, pure gravitational collapse via concentration in pressure bumps (such as would occur at sufficiently high concentrations and without the SI) is not responsible for planetesimal formation. For millimeter-sized particles, we find tentative evidence that planetesimal formation does not occur. If this result is confirmed at higher resolution, it could put strong constraints on where planetesimals can form. Ultimately, our results show that for centimeter-sized particles planetesimal formation in pressure bumps is extremely robust.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The Fuzzball Paradigm

We describe the puzzles that arise in the quantum theory of black holes, and explain how they are resolved in string theory. We review how the Bekenstein entropy is obtained through the count of brane bound states. We describe the fuzzball construction of black hole microstates. These states have no horizon and radiate from their surface like a normal body, so there is no information puzzle. We explain how the semiclassical approximation is violated in gravitational collapse even though curvatures are low at the classical horizon. This violation happens because the collapse leads to a stretching of space that is {\it fast}: light does not have time to travel across the collapsing region to establish the `vecro' correlations needed in the quantum gravitational vacuum. These vecro correlations arise from the existence of virtual fuzzball fluctuations in the gravitational vacuum, and are significant because of the large degeneracy of fuzzball states implied by the Bekenstein entropy. It is plausible that similar effects of fast expansion may be responsible for effects like dark energy and the Early Dark Energy postulated to explain the Hubble tension.

FOS: Physical sciences↗
DOE OSTI2025

Connecting GRBs from binary neutron star mergers to nuclear properties of neutron stars

The fate of the binary neutron star (NS) merger remnants hinges sensitively upon the NS equation of state and the threshold mass, M ls , that separates a long-lived from a short-lived NS remnant. The nature of the electromagnetic counterparts is also influenced by the remnant type, particularly in determining whether a gamma-ray burst from a compact binary merger (cbGRB) is of short or long duration. We propose a novel approach to probe Mls by linking it to the estimated observed ratio of long to short cbGRBs. We find that current observations broadly favor a relatively high value for this transition, M ls ≃ 1.3M TOV , for which M TOV ≲ 2.6M ⊙ , consistent with numerical simulations, as also shown here. Furthermore, our results disfavor nuclear physics scenarios that would lead to catastrophic pressure loss at a few times nuclear density and temperatures of tens of MeV, leading to a rapid gravitational collapse of binaries with total mass M ≲ 1.3M TOV . Future individual gravitational wave events with on-axis cbGRBs can further bound Mls.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Solar mass primordial black holes in moduli dominated universe

In this report we explore the prospect of producing primordial black holes around the solar mass region during an early matter domination epoch. The early matter-dominated epoch can arise when a moduli field comes to dominate the energy density of the Universe prior to big bang nucleosynthesis. The absence of radiation pressure during a matter-dominated epoch enhances primordial black hole formation from the gravitational collapse of primordial density fluctuations. In particular, we find that primordial black holes are produced in the 0.1 - 10M ⊙ mass range with a favorable choice of parameters in the theory. However, they cannot explain all of the merger events detected by the LIGO/Virgo gravitational wave search. In such a case, primordial black holes form about 4% of the total dark matter abundance, of which 95% belongs to the LIGO/Virgo consistent mass range. The rest of the dark matter could be in the form of particles that are produced from the decay of the moduli field during reheating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

From capture to collapse: Revisiting black hole formation by fermionic asymmetric dark matter in neutron stars

Fermionic asymmetric dark matter (ADM) can be captured in neutron stars (NSs) via scatterings with the star constituents. The absence of dark matter annihilation due to its asymmetric nature leads to ADM accumulation in the NS core, potentially reaching densities sufficient to exceed the Chandrasekhar limit and trigger its gravitational collapse into a black hole (BH), eventually consuming the NS from within. Therefore, the existence and observation of old neutron stars provide a means to constrain the properties of ADM. We revisit previous constraints on the mass and scattering cross section off neutrons of fermionic ADM across a class of models. We critically examine common simplifying approximations used in the literature to derive these limits. Our analysis includes improved treatments of dark matter capture, thermalization, BH formation, accretion, and evaporation. We find that previous results can be relaxed by a few orders of magnitude once these effects are properly accounted for.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Testing $f(R)$ gravity with scale dependent cosmic void velocity profiles

Here, we study the impact of cosmological scale modifications to general relativity on the dynamics of halos within voids by comparing N-body simulations incorporating Hu-Sawicki f(R) gravity, with |f R0 | = 10 -6 and 10 -5 , to those of Λ CDM. By examining the radial velocity statistics within voids classified based on their size and density-profile, as “rising” (R-type) or “shell” (S-type), we find that halo motions in small R-type voids, with effective radius <15 Mpc/h, reveal distinctive differences between f(R) and Λ CDM cosmologies. To understand this observed effect, we study the linear and nonlinear fifth forces, and develop an iterative algorithm to accurately solve the nonlinear fifth force equation. We use this to characterize the Chameleon screening mechanism in voids and contrast the behavior with that observed in gravitationally collapsed objects. The force analysis underscores how smaller R-type voids exhibit the highest ratios of fifth force to Newtonian force, which source distinguishable differences in the velocity profiles and thereby provide rich environments in which to constrain gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dust stars in the minimal exponential measure model

Here, we report the existence of horizonless compact object solutions supported by dust in the minimal exponential measure (MEMe) model, a theory which modifies the couplings between gravity and matter without introducing dynamical degrees of freedom. For a perfect fluid source, the field equations for the MEMe model can be rewritten as the Einstein field equations sourced by a perfect fluid with a transformed equation of state, which can endow a sufficiently dense cloud of dust with an effective pressure. The resulting dust-supported horizonless compact objects can have masses below ∼10 −11 ⁢𝑀 ⊙ , making them suitable as massive compact halo objects comprising a significant mass fraction for dark matter. A necessary condition for the existence of these compact object solutions is that the single free parameter in the MEMe model is positive valued. Additionally, we find that this positive sign for the parameter can provide a mechanism for suppressing the formation of (primordial) black holes from the gravitational collapse of matter below a certain mass scale.

alternative gravity theories↗
DOE OSTI2021

Replica wormholes for an evaporating 2D black hole

Quantum extremal islands reproduce the unitary Page curve of an evaporating black hole. This has been derived by including replica wormholes in the gravitational path integral, but for the transient, evaporating black holes most relevant to Hawking’s paradox, these wormholes have not been analyzed in any detail. In this paper we study replica wormholes for black holes formed by gravitational collapse in Jackiw-Teitelboim gravity, and confirm that they lead to the island rule for the entropy. The main technical challenge is that replica wormholes rely on a Euclidean path integral, while the quantum extremal islands of an evaporating black hole exist only in Lorentzian signature. Furthermore, the Euclidean equations for the Schwarzian mode are non-local, so it is unclear how to connect to the local, Lorentzian dynamics of an evaporating black hole. We address these issues with Schwinger-Keldysh techniques and show how the non-local equations reduce to the local ‘boundary particle’ description in special cases.

2D gravity↗
DOE OSTI2022

The role of magnetic fields in the stability and fragmentation of filamentary molecular clouds: two case studies at OMC-3 and OMC-4

We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained 214 $\mu$m HAWC + /SOFIA polarization data, and we make use of archival data for the column density and C 18 O (1–0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anticorrelated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientation (HRO) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line of sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavourable for star formation and is consistent with much fewer YSOs than in OMC-3.

79 ASTRONOMY AND ASTROPHYSICS↗