Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accretion”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

How Does Feedback Affect the Star Formation Histories of Galaxies?

Star formation in galaxies is regulated by the interplay of a range of processes that shape the multiphase gas in the interstellar and circumgalactic media. Using the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) suite of cosmological simulations, we study the effects of varying feedback and cosmology on the average star formation histories (SFHs) of galaxies at z ∼ 0 across the IllustrisTNG, SIMBA, and ASTRID galaxy formation models. We find that galaxy SFHs in all three models are sensitive to changes in stellar feedback, which affect the efficiency of baryon cycling and the rates at which central black holes grow, whereas the effects of varying active galactic nucleus (AGN) feedback depend on model-specific implementations of black hole seeding, accretion, and feedback. We also find strong interaction terms that couple stellar and AGN feedback, usually by regulating the amount of gas available for the central black hole to accrete. Using a double power law to describe the average SFHs, we derive a general set of equations relating the shape of the SFHs to physical quantities like baryon fraction and black hole mass across all three models. We find that a single set of equations (albeit with different coefficients) can describe the SFHs across all three CAMELS models, with cosmology dominating the SFH at early times, followed by halo accretion, and feedback and baryon cycling at late times. Galaxy SFHs provide a novel, complementary probe to constrain cosmology and feedback, and can connect the observational constraints from current and upcoming galaxy surveys with the physical mechanisms responsible for regulating galaxy growth and quenching.

Iyer, Kartheik G. [Columbia Univ., New York, NY (U↗
DOE OSTI2023

Collapsar Black Holes Are Likely Born Slowly Spinning

Collapsing stars constitute the main black hole (BH) formation channel, and are occasionally associated with the launch of relativistic jets that power γ-ray bursts (GRBs). Thus, collapsars offer an opportunity to infer the natal (before spin-up/down by accretion) BH spin directly from observations. We show that once the BH saturates with a large-scale magnetic flux, the jet power is dictated by the BH spin and mass accretion rate. Core-collapse simulations by Halevi et al. and GRB observations favor stellar density profiles that yield an accretion rate of $\dot m\approx 10^{-2} M_{\odot}s^{-1}$, weakly dependent on time. This leaves the spin as the main factor that governs the jet power. By comparing the jet power to characteristic GRB luminosities, we find that the majority of BHs associated with jets are likely born slowly spinning with a dimensionless spin of a ≃ 0.2, or a ≃ 0.5 for wobbling jets, with the main uncertainty originating in the unknown γ-ray radiative efficiency. This result could be applied to the entire core-collapse BH population, unless an anticorrelation between the stellar magnetic field and angular momentum is present. In a companion paper, Jacquemin-Ide et al., we show that regardless of the natal spin, the extraction of BH rotational energy leads to spin-down to a ≲ 0.2, consistent with gravitational-wave observations. We verify our results by performing the first 3D general-relativistic magnetohydrodynamic simulations of collapsar jets with characteristic GRB energies, powered by slowly spinning BHs. We find that jets of typical GRB power struggle to escape from the star, providing the first numerical indication that many jets fail to generate a GRB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

First Sagittarius A* Event Horizon Telescope Results. VIII. Physical Interpretation of the Polarized Ring

Abstract In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

XRISM Spectroscopy of the Fe Kα Emission Line in the Seyfert Active Galactic Nucleus NGC 4151 Reveals the Disk, Broad-line Region, and Torus

Abstract We present an analysis of the first two XRISM/Resolve spectra of the well-known Seyfert-1.5 active galactic nucleus (AGN) in NGC 4151, obtained in 2023 December. Our work focuses on the nature of the narrow Fe K α emission line at 6.4 keV, the strongest and most common X-ray line observed in AGN. The total line is found to consist of three components. Even the narrowest component of the line is resolved with evident Fe K α ,1 (6.404 keV) and K α ,2 (6.391 keV) contributions in a 2:1 flux ratio, fully consistent with neutral gas with negligible bulk velocity. Subject to the limitations of our models, the narrowest and intermediate-width components are consistent with emission from optically thin gas, suggesting that they arise in a disk atmosphere and/or wind. Modeling the three line components in terms of Keplerian broadening, they are readily associated with (1) the inner wall of the “torus,” (2) the innermost optical “broad-line region” (or “X-ray BLR”), and (3) a region with a radius of r ≃ 100 GM / c 2 that may signal a warp in the accretion disk. Viable alternative explanations of the broadest component include a fast-wind component and/or scattering; however, we find evidence of variability in the narrow Fe K α line complex on timescales consistent with small radii. The best-fit models are statistically superior to simple Voigt functions, but when fit with Voigt profiles the time-averaged lines are consistent with a projected velocity broadening of FWHM = 1600 − 200 + 400 km s − 1 . Overall, the resolution and sensitivity of XRISM show that the narrow Fe K line in AGN is an effective probe of all key parts of the accretion flow, as it is currently understood. We discuss the implications of these findings for our understanding of AGN accretion, future studies with XRISM, and X-ray-based black hole mass measurements.

Audard, Marc (ORCID:000000034721034X)↗
DOE OSTI2022

Reduced Atmospheres of Post-impact Worlds: The Early Earth

Impacts may have had a significant effect on the atmospheric chemistry of the early Earth. Reduced phases in the impactor (e.g., metallic iron) can reduce the planet's H 2 O inventory to produce massive atmospheres rich in H 2 . While previous studies have focused on the interactions between the impactor and atmosphere in such scenarios, we investigate two further effects: (1) the distribution of the impactor's iron inventory during impact between the target interior, target atmosphere, and escaping the target; and (2) interactions between the post-impact atmosphere and the impact-generated melt phase. We find that these two effects can potentially counterbalance each other, with the melt–atmosphere interactions acting to restore reducing power to the atmosphere that was initially accreted by the melt phase. For a ~10 22 kg impactor, when the iron accreted by the melt phase is fully available to reduce this melt, we find an equilibrium atmosphere with H 2 column density ~10 4 moles cm –2 (pH 2 ~ 120 bars, X H2 ~ 0.77), consistent with previous estimates. However, when the iron is not available to reduce the melt (e.g., sinking out in large diameter blobs), we find significantly less H 2 (7 × 10 2 – 5 × 10 3 moles cm –2 , pH 2 ≲ 60 bars, X H2 ≲ 0.41). These lower H 2 abundances are sufficiently high that species important to prebiotic chemistry can form (e.g., NH 3 , HCN), but sufficiently low that the greenhouse heating effects associated with highly reducing atmospheres, which are problematic to such chemistry, are suppressed. The manner in which iron is accreted by the impact-generated melt phase is critical in determining the reducing power of the atmosphere and resolidified melt pool in the aftermath of impact.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Origin of the ring ellipticity in the black hole images of M87*

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of τ = 0.08−0.02+0.03 with a position angle of ξ = 50.1−7.6+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87* emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.Key words: accretion, accretion disks / black hole physics / gravitation / galaxies: active** These authors have contributed equally to this work.*** NASA Hubble Fellowship Program, Einstein Fellow.† Deceased.

Dahale, Rohan [Granada U., Theor. Phys. Astrophys.↗
DOE OSTI2023

Evaluation of aerosol–cloud interactions in E3SM using a Lagrangian framework

A Lagrangian framework is used to evaluate aerosol–cloud interactions in the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 1 (E3SMv1) for measurements taken at Graciosa Island in the Azores where a U.S. Department of Energy Atmosphere Radiation Measurement (ARM) site is located. This framework uses direct measurements of cloud condensation nuclei (CCN) concentration (instead of relying on satellite retrievals of aerosol optical depth) and incorporates a suite of ground-based ARM measurements, satellite retrievals, and meteorological reanalysis products that when applied to over a 1500 trajectories provides key insights into the evolution of low-level clouds and aerosol radiative forcing that is not feasible from a traditional Eulerian analysis framework. Significantly lower concentrations (40 %) of surface CCN concentration are measured when precipitation rates in 48 h back trajectories average above 1.2 mm d -1 in the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG) product. The depletion of CCN concentration when precipitation rates are elevated is nearly twice as large in the ARM observations compared to E3SMv1 simulations. The model CCN concentration bias remains significant despite modifying the autoconversion and accretion rates in warm clouds. As the clouds in trajectories associated with larger surface-based CCN concentration advect away from Graciosa Island, they maintain higher values of droplet number concentrations (N d ) over multiple days in observations and E3SM simulations compared to trajectories that start with lower CCN concentrations. The response remains robust even after controlling for meteorological factors such as lower troposphere stability, the degree of cloud coupling with the surface, and island wake effects. E3SMv1 simulates a multi-day aerosol effect on clouds and a Twomey radiative effect that is within 30 % of the ARM and satellite observations. However, the mean cloud droplet concentration is more than 2–3 times larger than in the observations. While Twomey radiative effects are similar amongst autoconversion and accretion sensitivity experiments, the liquid water path and cloud fraction adjustments are positive when using a regression model as opposed to negative when using the present-day minus pre-industrial aerosol emissions approach. This result suggests that tuning the autoconversion and accretion alone is unlikely to produce the desired aerosol susceptibilities in E3SMv1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

eROSITA Final Equatorial-Depth Survey (eFEDS) eFEDS X-ray view of WERGS radio galaxies selected by the Subaru/HSC and VLA/FIRST survey

We constructed the eROSITA X-ray catalog of radio galaxies discovered by the WERGS survey, made by cross-matching the wide area Subaru/Hyper Suprime-Cam (HSC) optical survey and VLA/FIRST 1.4 GHz radio survey. We report finding 393 eROSITA detected radio galaxies in the 0.5-2 keV band in the eFEDS field covering 140 deg 2 . Thanks to the wide and medium depth eFEDS X-ray survey down to f 0.5 - 2 keV = 6.5 × 10 -15 erg s -1 cm -2 , the sample contains the rare and most X-ray luminous radio galaxies above the knee of the X-ray luminosity function, spanning 44 < log(L 0.5-2 keV (abs,corr) /erg s -1 ) < 46.5 at 1 < z < 4. The sample also contains the sources around and below the knee for the sources 41 < log(L 0.5-2 keV (abs,corr) /erg s -1 ) < 45 at z < 1. Based on the X-ray properties obtained by the spectral fitting, 37 sources show obscured active galactic nucleus (AGN) signatures with log(N H /cm -2 ) > 22. These obscured and radio AGN reside in 0.4 < z < 3.2, indicating that they are obscured counterparts of the radio-loud quasar, which were missed in the previous optical quasar surveys. By combining radio and X-ray luminosities, we also investigated the jet production efficiency η jet = η rad P jet /L AGN,bol by utilizing the jet power of P jet . We find that there are 14 sources with extremely high jet production efficiency at η jet ≈ 1. This high η jet value might be a result of the decreased radiation efficiency of η rad < 0.1, due to the low accretion rate for those sources, and/or of the boosting due to the decline of L AGN,bol by a factor of 10-100 by keeping P jet constant in the previous Myr, indicating the experience of the AGN feedback. Finally, inferring the BH masses from the stellar mass, we find that X-ray luminous sources show the excess of the radio emission with respect to the value estimated from the fundamental plane. This radio emission excess cannot be explained by the Doppler boosting alone, and therefore the disk–jet connection of X-ray luminous eFEDS-WERGS is fundamentally different from the conventional fundamental plane which mainly covers the low-accretion regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The antesonic condition for the explosion of core-collapse supernovae – II. Rotation and turbulence

In the problem of steady free fall on to a standing shockwave around a central mass, the ‘antesonic’ condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\le 3/16$, where cT is the isothermal sound speed in the subsonic post-shock flow, and vesc is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and the shock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models of accretion 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 ~5 percent 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 ~20 to 40 percent. 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 multidimensional simulations explode more easily than their one-dimensional counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The SOUX AGN sample: optical/UV/X-ray SEDs and the nature of the disc

We use the SOUX sample of ~700 active galactic nucleus (AGN) to form average optical-ultraviolet (UV)-X-rays spectral energy distributions (SEDs) on a two-dimensional (2D) grid of M BH and L 2500 . We compare these with the predictions of a new AGN SED model, qsosed, which includes prescriptions for both hot and warm Comptonization regions as well as an outer standard disc. This predicts the overall SED fairly well for 7.5 < log(M BH /M ⊙ ) < 9.0 over a wide range in L/L Edd , but at higher masses the outer disc spectra in the model are far too cool to match the data. We create optical-UV composites from the entire Sloan Digital Sky Survey sample and use these to show that the mismatch is due to there being no significant change in spectral shape of the optical-UV continuum across several decades of M BH at constant luminosity. We show for the first time that this cannot be matched by standard disc models with high black hole spin. These apparently fit, but are not self-consistent as they do not include the General Relativistic effects for the emission to reach the observer. At high spin, increased gravitational redshift compensates for almost all of the higher temperature emission from the smaller inner disc radii. The data do not match the predictions made by any current accretion flow model. Either the disc is completely covered by a warm Comptonization layer whose properties change systematically with L/L Edd , or the accretion flow structure is fundamentally different to that of the standard disc models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

XRISM high-resolution X-ray spectroscopy of Cygnus X-1: Highly ionized iron absorption structures

We present the first high-resolution X-ray spectral analysis of Cygnus X-1 using XRISM. The observation wa3s carried out from 2024 April 7 to 10, covering the orbital phase range 0.65–0.17 during its low/hard state. Taking advantage of the exceptional energy resolution of the Resolve instrument, we examined highly ionized iron absorption lines and characterized the ionization states, column densities, and line-of-sight velocities of the absorbing plasma. Spectral analysis revealed an ionization parameter of $\xi \sim 3$, column densities of a few $\times 10^{21}$ cm$^{-2}$, and a blueshifted velocity of $\sim$100 km s$^{-1}$. The observation was divided into two phases: before and after orbital phase $\phi _{\rm {orb}} = 0.9$, corresponding to non-dipping and dipping intervals. While only weak absorption features were present before $\phi _{\rm {orb}} = 0.9$, strong absorption by He-like and H-like Fe appeared during the dipping phase. We measured equivalent widths of 2.3, 0.4, and 1.2 eV for He-like Fe K$\alpha$ and H-like Ly$\alpha _1$ and Ly$\alpha _2$, respectively—demonstrating the capability of XRISM Resolve to securely detect narrow absorption features of only a few eV. These measurements trace the motion of the absorbing material and offer insight into the kinematics and spatial distribution of the wind in the vicinity of the black hole. These findings enhance our understanding of wind-fed accretion in Cygnus X-1 and highlight the importance of continued high-resolution X-ray observations to further constrain the physical properties of winds and accretion flows in high-mass X-ray binaries.

Astronomy and AstroPhysics↗
DOE OSTI2026

XRISM high-resolution X-ray spectroscopy of Cygnus X-1: Orbital and short-term variability of iron absorption

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe K band on second timescales. Using data from the Performance Verification phase in 2024 April, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser, and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve’s capability to probe plasma environments in high-mass X-ray binaries.

Astronomy and AstroPhysics↗
DOE OSTI2026

Changing-look Active Galactic Nuclei Powered by Disk Tearing

Changing-look active galactic nuclei (CLAGN) feature order-of-magnitude variability in both the continuum and broad line luminosities on months-to-years-long timescales, and are currently unexplained. Simulations have demonstrated that rotating black holes (BHs) sometimes tear apart tilted accretion disks. These tearing events violently restructure the disk on timescales much shorter than a viscous timescale, hinting at a connection to CLAGN. Here, we show that disk tearing can power changing-look events. We report synthetic observations of an extremely high-resolution, three-dimensional, general-relativistic magnetohydrodynamic simulation of a geometrically thin, tilted accretion disk around a rapidly rotating, 10 8 M ⊙ BH. We perform ray-tracing calculations that follow the thermal disk emission to both a line-of-sight camera and to a distribution of cameras in a prescribed torus-like broad-line region (BLR). The continuum photoionizes the BLR, and we calculate the resulting spectrum. Both the continuum and line luminosities undergo order-of-magnitude swings on months-to-years-long timescales. We find shorter, week-long variability driven by the geometric precession of the inner disk and an intraday quasiperiodic oscillation driven by radial breathing of the inner disk. When the torn disk precesses, it causes asymmetric illumination of the BLR, driving time-evolving red-to-blue asymmetries of the broad emission lines that may be a smoking gun for disk tearing. We also make predictions for future photometric observations from ULTRASAT and Vera Rubin Observatory, given a Shakura–Sunyaev–like effective temperature profile, both of which may play an important role in detecting future changing-look events. Future work should consider nonthermal emission, especially from the time-variable corona.

accretion↗
DOE OSTI2020

The Non-carbonaceous–Carbonaceous Meteorite Dichotomy

The isotopic dichotomy between non-carbonaceous (NC) and carbonaceous (CC) meteorites indicates that meteorite parent bodies derive from two genetically distinct reservoirs, which presumably were located inside (NC) and outside (CC) the orbit of Jupiter and remained isolated from each other for the first few million years of the solar system. Here we review the discovery of the NC–CC dichotomy and its implications for understanding the early history of the solar system, including the formation of Jupiter, the dynamics of terrestrial planet formation, and the origin and nature of Earth’s building blocks. The isotopic difference between the NC and CC reservoirs is probably inherited from the solar system’s parental molecular cloud and has been maintained through the rapid formation of Jupiter that prevented significant exchange of material from inside (NC) and outside (CC) its orbit. The growth and/or migration of Jupiter resulted in inward scattering of CC bodies, which accounts for the co-occurrence of NC and CC bodies in the present-day asteroid belt and the delivery of presumably volatile-rich CC bodies to the growing terrestrial planets. Earth’s primitive mantle, at least for siderophile elements like Mo, has a mixed NC–CC composition, indicating that Earth accreted CC bodies during the final stages of its growth, perhaps through the Moon-forming giant impactor. The late-stage accretion of CC bodies to Earth is sufficient to account for the entire budget of Earth’s water and highly volatile species.

58 GEOSCIENCES↗
DOE OSTI2022

A hierarchical stochastic modeling approach for representing point bar geometries and petrophysical property variations

The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.

58 GEOSCIENCES↗
DOE OSTI2025

The persistent shadow of the supermassive black hole of M87. II. Model comparisons and theoretical interpretations

The Event Horizon Telescope (EHT) observation of M87∗ in 2018 has revealed a ring with a diameter that is consistent with the 2017 observation. The brightest part of the ring is shifted to the southwest from the southeast. In this paper, we provide theoretical interpretations for the multi-epoch EHT observations for M87∗ by comparing a new general relativistic magnetohydrodynamics model image library with the EHT observations for M87∗ in both 2017 and 2018. The model images include aligned and tilted accretion with parameterized thermal and nonthermal synchrotron emission properties. The 2018 observation again shows that the spin vector of the M87∗ supermassive black hole is pointed away from Earth. A shift of the brightest part of the ring during the multi-epoch observations can naturally be explained by the turbulent nature of black hole accretion, which is supported by the fact that the more turbulent retrograde models can explain the multi-epoch observations better than the prograde models. The EHT data are inconsistent with the tilted models in our model image library. Assuming that the black hole spin axis and its large-scale jet direction are roughly aligned, we expect the brightest part of the ring to be most commonly observed 90 deg clockwise from the forward jet. This prediction can be statistically tested through future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗
DOE OSTI2020

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗