Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BH”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

BSE versus StarTrack: Implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N-body evolution program NBODY7. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. All the above newly implemented schemes are also shared with the standalone versions of SSE and BSE. Furthermore, all these demonstrations are performed with both the updated standalone BSE and the updated NBODY7/BSE. When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing. Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: I. Multi-coordinate implementation and testing

General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole (BH) shadows from the Galactic Center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as BHs and neutron stars. To do so, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. Here, we provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the E INSTEIN T OOLKIT . The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of BH magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Black hole and galaxy co-evolution in radio-loud active galactic nuclei at z ~ 0.3–4

The relation between the mass of the supermassive black hole (SMBH) in the center of galaxies and their bulge mass or central velocity dispersion is well known. This suggests a coevolution between the SMBHs and their galaxy hosts. Our aim is to study this relation, specifically, for radio loud galaxies, and as a function of redshift z. We selected a sample of 42 radio galaxies and active galactic nuclei (AGN) with broad emission lines and spectroscopic redshifts between z = 0.3 - 4 by cross-matching the low radio frequency sources from Very Large Array (VLA) FIRST with spectroscopically confirmed galaxies from wide-field surveys, including Sloan Digital Sky Survey (SDSS) DR14 ugriz and Dark Energy Survey (DES) DR2 grzY in the optical, Wield Infrared Survey Explorer (WISE), and the Galaxy And Mass Assembly (GAMA) spectroscopic survey. We characterized the stellar mass (M $\star$ ), star formation, and black hole properties (mass of the central SMBH, Eddington ratio η, and jet power, Q jet ). The relation between SMBH mass, M $\star$ , η, and z is placed into context by comparing them with scaling relations (M BH –M $\star$ , M BH /M $\star$ –z, M BH –Q jet , and Q jet –η) from the literature. On the basis of a multiwavelength spectral energy distribution modeling, our radio sources are broadly consistent with being on the star-forming main sequence. They have sub-Eddington accretion rates, η ≃ 1% on average, as typically found in type I AGN, while higher accretion rates favor more powerful jets to be launched by the central engine. We find overmassive SMBHs in (17 ± 5)% of our radio sources, similarly to previous studies on nearby early-type galaxies. Altogether, an evolutionary scenario in which radio-mode AGN feedback regulates the accretion onto the SMBHs and the stellar mass assembly of the radio sources is discussed, which may explain the observed phenomenology. This pilot study represents a benchmark for future studies using wide-field surveys such as those with Euclid and the Vera Rubin Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The miniJPAS survey: AGN and host galaxy coevolution of X-ray-selected sources

Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass (MBH) and the stellar mass of their host galaxies (M * ). They even show similar histories across cosmic times of their differential terms: the star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being understood. We selected an X-ray sample of active galactic nuclei (AGN) up to z = 2.5 in the miniJPAS footprint. Their X-ray to infrared spectral energy distributions (SEDs) have been modeled with the CIGALE code, constraining the emission to 68 bands, from which 54 are the narrow filters from the miniJPAS survey. For a final sample of 308 galaxies, we derived their physical properties, such as their M * , SFR, star formation history (SFH), and the luminosity produced by the accretion process of the central BH (L AGN ). For a subsample of 113 sources, we also fit their optical spectra to obtain the gas velocity dispersion from the broad emission lines and estimated the M BH . We calculated the BHAR in physical units depending on two radiative efficiency regimes. We find that the Eddington ratios (λ Edd ) and its popular proxy (L X /M * ) have a difference of 0.6 dex, on average, and a KS test indicates that they come from different distributions. Our sources exhibit a considerable scatter on the MBH - M * scaling relation, which can explain the difference between λ Edd and its proxy. We also modeled three evolution scenarios for each source to recover the integral properties at z = 0. Using the SFR and BHAR, we show a notable diminution in the scattering between MBH - M * . For the last scenario, we considered the SFH and a simple energy budget for the AGN accretion, and we retrieved a relation similar to the calibrations known for the local Universe. Our study covers ~1 deg 2 in the sky and is sensitive to biases in luminosity. Nevertheless, we show that, for bright sources, the link between the differential values (SFR and BHAR) and their decoupling based on an energy limit is the key that leads to the local M BH - M * scaling relation. In the future, we plan to extend this methodology to a thousand degrees of the sky using JPAS with an X-ray selection from eROSITA, to obtain an unbiased distribution of BHAR and Eddington ratios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Galaxy mergers in eagle do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN

ABSTRACT We investigate the connection between galaxy–galaxy mergers and enhanced black hole (BH) growth using the cosmological hydrodynamical eagle simulation. We do this via three methods of analysis, investigating: the merger fraction of AGN, the AGN fraction of merging systems, and the AGN fraction of galaxies with close companions. In each case, we find an increased abundance of AGN within merging systems relative to control samples of inactive or isolated galaxies (by up to a factor of ≈3 depending on the analysis method used), confirming that mergers are enhancing BH accretion rates for at least a subset of the galaxy population. The greatest excess of AGN triggered via a merger are found in lower mass (M* ∼ 1010 M⊙) gas rich (fgas > 0.2) central galaxies with lower mass BHs (MBH ∼ 107 M⊙) at lower redshifts (z < 1). We find no enhancement of AGN triggered via mergers in more massive galaxies (M* ≳ 1011 M⊙). The enhancement of AGN is not uniform throughout the phases of a merger, and instead peaks within the early remnants of merging systems (typically lagging ≈300 Myr post-coalescence of the two galaxies at z = 0.5). We argue that neither major ($M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \ge \frac{1}{4}$) nor minor mergers ($\frac{1}{10} \lt M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \lt \frac{1}{4}$) are statistically relevant for enhancing BH masses globally. Whilst at all redshifts the galaxies experiencing a merger have accretion rates that are on average 2–3 times that of isolated galaxies, the majority of mass that is accreted on to BHs occurs outside the periods of a merger. We compute that on average no more than 15 per cent of a BHs final day mass comes from the enhanced accretion rates triggered via a merger.

McAlpine, Stuart↗
DOE OSTI2020

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first ~1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass haloes. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant, and scaling with stellar mass) from both star formation ( |$\langle f_{\rm esc}^{\rm sf} \rangle$| ) and AGN ( |$f_{\rm esc}^{\rm bh}$| ) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ( |$M_*\lt 10^9 \, {\rm \rm M_\odot }$| ) galaxies at z > 6 with the AGN contribution (driven by |$M_{bh}\gt 10^6 \, {\rm \rm M_\odot }$| black holes in |$M_* \gtrsim 10^9\, {\rm \rm M_\odot}$| galaxies) dominating at lower redshifts; (ii) AGN only contribute |$10-25{{\ \rm per\ cent}}$| to the cumulative ionizing emissivity by z = 4 for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of |$\langle f_{\rm esc}^{\rm sf} \rangle$| is shallower than |$f_{\rm esc}^{\rm bh}$| , at z < 7 a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at z ~ 6.8, they take-over at stellar masses that are a tenth of the knee mass by z = 4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dwarf AGNs from Optical Variability for the Origins of Seeds (DAVOS): insights from the dark energy survey deep fields

ABSTRACT We present a sample of 706, z < 1.5 active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg2. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 d, a 6-yr baseline, and a single-epoch imaging depth of up to g ∼ 24.5. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast }\lt 10^{9.5}\, {\rm M}_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability time-scales (∼ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole (BH) masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalogue, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower BH masses in deep fields, which may be more ‘pristine’ analogues of supermassive BH seeds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Long-term 3D MHD simulations of black hole accretion discs formed in neutron star mergers

ABSTRACT We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in flash 4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disc simulations that vary the initial magnetic field geometry and disc compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modelling the BH’s gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively inefficient state in the disc. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disc compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively inefficient state at late times. The electron fraction distribution of the disc outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding ${\sim} 50{{\ \rm per\ cent}}$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

Fahlman, Steven (ORCID:0000000348759940)↗
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 OSTI2023

Secular outflows from 3D MHD hypermassive neutron star accretion disc systems

Magnetized hypermassive neutron stars (HMNSs) have been proposed as a way for neutron star mergers to produce high electron fraction, high-velocity ejecta, as required by kilonova models to explain the observed light curve of GW170817. The HMNS drives outflows through neutrino energy deposition and mechanical oscillations, and raises the electron fraction of outflows through neutrino interactions before collapsing to a black hole (BH). Here, in this study, we perform 3D numerical simulations of HMNS–torus systems in ideal magnetohydrodynamics, using a leakage/absorption scheme for neutrino transport, the nuclear APR equation of state, and Newtonian self-gravity, with a pseudo-Newtonian potential added after BH formation. Due to the uncertainty in the HMNS collapse time, we choose two different parametrized times to induce collapse. We also explore two initial magnetic field geometries in the torus, and evolve the systems until the outflows diminish significantly ($\sim\!\! 1\!\! - \!\!2\ \mathrm{s}$). We find bluer, faster outflows as compared to equivalent BH–torus systems, producing M ∼ 10 −3 M ⊙ of ejecta with Y e ≥ 0.25 and v ≥ 0.25c by the simulation end. Approximately half the outflows are launched in disc winds at times $t\lesssim 500 \ \mathrm{ms}$, with a broad distribution of electron fractions and velocities, depending on the initial condition. The remaining outflows are thermally driven, characterized by lower velocities and electron fractions. Nucleosynthesis with tracer particles shows patterns resembling solar abundances in all models. Although outflows from our simulations do not match those inferred from two-component modelling of the GW170817 kilonova, self-consistent multidimensional detailed kilonova models are required to determine whether our outflows can power the blue kilonova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

GW190521 from the Merger of Ultradwarf Galaxies

We present an alternative formation scenario for the gravitational wave event GW190521 that can be explained as the merger of central black holes (BHs) from two ultradwarf galaxies of stellar mass ∼105–106 M⊙, which had themselves previously undergone a merger. The GW190521 components’ masses of 85-14+21 M⊙ and 66-18+17 M⊙ challenge standard stellar evolution models, as they fall in the so-called mass gap. We demonstrate that the merger history of ultradwarf galaxies at high redshifts (1≲z≲2) matches well the LIGO-Virgo inferred merger rate for BHs within the mass range of the GW190521 components, resulting in a likely time delay of ≲4 Gyr considering the redshift of this event. We further demonstrate that the predicted timescales are consistent with expectations for central BH mergers, although with large uncertainties due to the lack of high-resolution simulations in low-mass dwarf galaxies. Our findings show that this BH production and merging channel is viable and extremely interesting as a new way to explore galaxies’ BH seeds and galaxy formation. We recommend this scenario be investigated in detail with simulations and observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Materials Data on RbB10H9 by Materials Project

Rb(BH)7HB2BH crystallizes in the orthorhombic Pna2_1 space group. The structure is zero-dimensional and consists of eight boranediylradical molecules, eight HB2 clusters, and four Rb(BH)7 clusters. In each HB2 cluster, there are two inequivalent B+0.60- sites. In the first B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.32 Å. In the second B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.32 Å. H+0.56+ is bonded in a water-like geometry to two B+0.60- atoms. In each Rb(BH)7 cluster, there are two inequivalent Rb1+ sites. In the first Rb1+ site, Rb1+ is bonded in a 8-coordinate geometry to eight H+0.56+ atoms. There are a spread of Rb–H bond distances ranging from 2.85–3.16 Å. In the second Rb1+ site, Rb1+ is bonded in a 7-coordinate geometry to seven H+0.56+ atoms. There are a spread of Rb–H bond distances ranging from 2.91–3.20 Å. There are fourteen inequivalent B+0.60- sites. In the first B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.19 Å. In the second B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the third B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the fourth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the fifth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the sixth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the seventh B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the eighth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the ninth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the tenth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the eleventh B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the twelfth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the thirteenth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.20 Å. In the fourteenth B+0.60- site, B+0.60- is bonded in a distorted single-bond geometry to one H+0.56+ atom. The B–H bond length is 1.19 Å. There are fourteen inequivalent H+0.56+ sites. In the first H+0.56+ site, H+0.56+ is bonded in a distorted single-bond geometry to one Rb1+ and one B+0.60- atom. In the second H+0.56+ site, H+0.56+ is bonded in a distorted single-bond geometry to one Rb1+ and one B+0.60- atom. In the third H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to two Rb1+ and one B+0.60- atom. In the fourth H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to one Rb1+ and one B+0.60- atom. In the fifth H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to one Rb1+ and one B+0.60- atom. In the sixth H+0.56+ site, H+0.56+ is bonded in a distorted single-bond geometry to one Rb1+ and one B+0.60- atom. In the seventh H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to one Rb1+ and one B+0.60- atom. In the eighth H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to one Rb1+ and one B+0.60- atom. In the ninth H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to one Rb1+ and one B+0.60- atom. In the tenth H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to one Rb1+ and one B+0.60- atom. In the eleventh H+0.56+ site, H+0.56+ is bonded in a distorted single-bond geometry to one Rb1+ and one B+0.60- atom. In the twelfth H+0.56+ site, H+0.56+ is bonded in a distorted single-bond geometry to one Rb1+ and one B+0.60- atom. In the thirteenth H+0.56+ site, H+0.56+ is bonded in a distorted single-bond geometry to one Rb1+ and one B+0.60- atom. In the fourteenth H+0.56+ site, H+0.56+ is bonded in a single-bond geometry to one Rb1+ and one B+0.60- atom.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cavern Leaching Monitoring CY18 And CY19

The U.S. Strategic Petroleum Reserve (SPR) is a crude oil storage system run by the U.S. Department of Energy (DOE). The reserve consists of 60 active storage caverns spread across four sites in Louisiana and Texas, near the Gulf of Mexico. Beginning in 2016, the SPR began executing U.S. congressionally mandated oil sales. The configuration of the reserve, with a total capacity of greater than 700 MMB, requires raw water to be used instead of saturated brine for oil withdrawals such as for sales. All sales will produce leaching within the caverns used for oil delivery. Thirty-six caverns had a combined total of over 29 MMB of water injected from CY18-CY19 for mandatory sales. Leaching effects were monitored in these caverns to understand how the sales operations may impact the long-term integrity of the caverns. While frequent sonars are the best way to monitor changes in cavern shape, they can be resource intensive for the number of caverns involved in sales and exchanges. An intermediate option is to model the leaching effects and see if any concerning features develop. The leaching effects were modeled here using the Sandia Solution Mining Code (SANSMIC). The results indicate that leaching induced features are not of concern in the majority of the caverns, 32 of 36. Four caverns, BH-107, BH-108, BH-114 and WH-114 have features that may grow with additional leaching and should be monitored as leaching continues in those caverns. Six caverns had post sale sonars which were compared with SANSMIC results. SANSMIC was able to capture the leaching well. A deviation in the SANSMIC and sonar cavern shapes was observed near the cavern floor in caverns with significant floor rise, a process not captured by SANSMIC. These results suggest SANSMIC is a useful tool for monitoring changes in cavern shape due to leaching effects related to sales and exchanges.

02 PETROLEUM↗
DOE OSTI2021

Strategic Petroleum Reserve Cavern Leaching Monitoring CY20

The U.S. Strategic Petroleum Reserve is a crude oil storage system run by the U.S. Department of Energy. The reserve consists of 60 active storage caverns spread across four sites in Louisiana and Texas, near the Gulf of Mexico. Beginning in 2016, the SPR began executing U.S. congressionally mandated oil sales. The configuration of the reserve, with a total capacity of greater than 700 MMB, requires raw water to be used instead of saturated brine for oil withdrawals such as for sales. All sales will produce leaching within the caverns used for oil delivery. Twenty-five caverns had a combined total of over 39 MMB of water injected in CY 20 as part of the Exchange for Storage program; oil was withdrawn in the same manner as for congressionally mandated sales. Leaching effects were monitored in these caverns to understand how the oil withdrawals may impact the long-term integrity of the caverns. While frequent sonars are the best way to monitor changes in cavern shape, they can be resource intensive for the number of caverns involved in sales and exchanges. An intermediate option is to model the leaching effects and see if any concerning features develop. The leaching effects were modeled here using the Sandia Solution Mining Code (SANSMIC) . The results indicate that leaching induced features are not of concern in the majority of the caverns, 19 of 25. Six caverns, BH-107, BH-113, BH-114, BM-4, BM-106, and WH-114 have features that may grow with additional leaching and should be monitored as leaching continues in those caverns. Ten caverns had post sale sonars that were compared with SANSMIC results. SANSMIC was able to capture the leaching well , particularly the formation of shelves and flares. A deviation in the SANSMIC and sonar cavern shapes was observed near the cavern floor in caverns with significant floor rise, a process not captured by SANSMIC. These results suggest SANSMIC is a useful tool for monitoring changes in cavern shape due to leaching effects related to sales and exchanges.

02 PETROLEUM↗
DOE OSTI2020

The Ubiquity of AGN Winds in Seyfert 1 Galaxies at Low Redshift

The optical spectra of 3896 Seyfert 1 (Sy1) galaxies detected with WISE at z < 0.4 were analyzed for evidence of outflows. In 37% of the Sy1s in our sample, the outflows appear as broad, blue-shifted, spectrally resolved components of the [O iii]λ5007 Å emission line, with a mean maximum velocity V {sub max} ∼ 1014 km s{sup −1} that is consistent with AGN winds. For each Sy1, we deduced the black hole (BH) mass, bolometric luminosity, Eddington ratio, and power-law index of the continuum, which we compared with the star formation rate (SFR) and host morphology. Having separated our sample in two spectroscopic subgroups—Sy1s with only broad Balmer lines (Sy1B) and those with both narrow and broad (Sy1N) lines—and distinguishing those that show as outflow (Sy1Bw and Sy1Nw), we report the following differences: (1) the BH mass is systematically higher and the power-law steeper in the Sy1B–Sy1Bw than in the Sy1N–Sy1Nw; (2) V {sub max} is higher in the Sy1Bw than in the Sy1Nw, correlated in both groups with the BH mass and bolometric luminosity; (3) the Eddington ratio and SFR are higher in the Sy1 with outflows; and (4) the specific star formation rates (sSFRs) of the Sy1s are normal for their morphology and mass, typical of early-type spiral galaxies in the green valley, far from the quenched regime. From these results, we conclude that AGN winds in Sy1s are triggered by higher accretion rates and probably radiatively launched, and there is no clear evidence of an effect on the star formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Equation of State and Progenitor Dependence of Stellar-mass Black Hole Formation

The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS, it will become gravitationally unstable and further collapse into a black hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of pre-supernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-specific entropy phase space ${M}_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and is almost EOS independent, apart from the final end point. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum gravitational mass supported by a hot NS with constant specific entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Quasar Shedding Its Dust Cocoon at Redshift 2

We present the first near-IR spectroscopy and joint analyses of multiwavelength observations for SDSS J082747.14+425241.1, a dust-reddened, weak broad emission-line quasar (WLQ) undergoing a remarkable broad-absorption line (BAL) transformation. The systemic redshift is more precisely measured to be z = 2.070 ± 0.001 using Hβ compared to z = 2.040 ± 0.003 using Mg ιι from the literature, signifying an extreme Mg ιι blueshift of 2140 ± 530 km s –1 relative to Hβ. Using the Hβ-based single-epoch scaling relation with a systematic uncertainty of 0.3 dex, its black hole (BH) mass and Eddington ratio are estimated to be M BH ~ 6.1 × 10 8 M ⊙ and λ Edd ~ 0.71, indicative of being in a rapidly accreting phase. Our investigations confirm the WLQ nature and the LoBAL → HiBAL transformation, along with a factor of 2 increase in the Mg ιι+Fe ιι emission strength and a decrease of 0.1 in E(B – V) over two decades. The kinetic power of this LoBAL wind at R ~ 15 pc from its BH is estimated to be ~43% of the Eddington luminosity, sufficient for quasar feedback upon its host galaxy albeit with an order-of-magnitude uncertainty. This quasar provides a clear example of the long-sought scenario where LoBAL quasars are surrounded by dust cocoons, and wide-angle nuclear winds play a key role in the transition of red quasars evolving into the commonly seen blue quasars.

79 ASTRONOMY AND ASTROPHYSICS↗