Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “log structured”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of structure and volcanic stratigraphy on groundwater and surface water flow: Hat Creek basin, California, USA

Abstract Hydrogeologic systems in the southern Cascade Range in California (USA) develop in volcanic rocks where morphology, stratigraphy, extensional structures, and attendant basin geometry play a central role in groundwater flow paths, groundwater/surface-water interactions, and spring discharge locations. High-volume springs (greater than 3 m 3 /s) flow from basin-filling (<800 ka) volcanic rocks in the Hat Creek and Fall River tributaries and contribute approximately half of the average annual flow of the Pit River, the largest tributary to Shasta Lake. A hydrogeologic conceptual framework is constructed for the Hat Creek basin combining new geologic mapping, water-well lithologic logs, a database of active faults, LiDAR mapping of faults and volcanic landforms, streamflow measurements and airborne thermal infrared remote sensing of stream temperature. These data are used to integrate the geologic structure and the volcanic and volcaniclastic stratigraphy to create a three-dimensional interpretation of the hydrogeology in the basin. Two large streamflow gains from focused groundwater discharge near Big Spring and north of Sugarloaf Peak result from geologic barriers that restrict lateral groundwater flow and force water into Hat Creek. The inferred groundwater-flow barriers divide the aquifer system into at least three leaky compartments. The two downstream compartments lose streamflow in the upstream reaches (immediately downstream of the groundwater-flow barriers) and gain in downstream reaches with the greatest inflows immediately upstream of the barriers.

Geology↗
DOE OSTI2024

Shallow Rate-Redox Potential Scaling in Aqueous Molecular Oxygen Reduction Electrocatalysis Across a Family of Iron Macrocycles

Rate-overpotential scaling relationships have been employed widely to understand trends in oxygen reduction reaction (ORR) electrocatalysis by dissolved metal macrocycles in organic electrolytes. Similar scaling relationships remain unknown for surface-adsorbed ORR electrocatalysts in the acidic aqueous environments germane to proton-exchange membrane (PEM) fuel cells. Herein, we examine ORR catalysis in aqueous perchloric acid media for a structurally diverse array of iron macrocycle complexes adsorbed on Vulcan carbon black. The macrocycles encompass Fe– N 4 , Fe–N 2 N' 2 and Fe–N x C 4-x motifs bearing pyrrolic, pyridinic, and N-heterocyclic carbene (NHC) moieties in the primary ligation sphere, giving rise to a 670 mV range in Fe(III/II) redox potentials, E Fe(III/II) . Experimental Tafel data in the micropolarization regime were extrapolated to the E Fe(III/II) to furnish estimated per-site-normalized current density (j per-site ) values that span ~4.6 orders of magnitude across the family of compounds. Despite the structural diversity of this family of compounds, extrapolated j per-site values correlate with the Fe(III/II) redox potentials in a roughly log-linear fashion with a shallow scaling factor of approximately 145 mV/decade. Further, these findings highlight that negative shifts in E Fe(III/II) lead to diminishing returns in catalytic rate promotion and suggest that changes to the primary ligating environment in a macrocycle are insufficient to break fundamental rate-potential scaling relationships in aqueous ORR catalysis. Together these studies motivate the further development of higher-potential iron complexes that employ motifs beyond the equatorial ligation plane to enhance ORR catalysis.

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

Predicting Volume of Distribution in Humans: Performance of In Silico Methods for a Large Set of Structurally Diverse Clinical Compounds

Volume of distribution at steady state (V D,ss ) is one of the key pharmacokinetic parameters estimated during the drug discovery process. Despite considerable efforts to predict V D,ss , accuracy and choice of prediction methods remain a challenge, with evaluations constrained to a small set (<150) of compounds. To address these issues, a series of in silico methods for predicting human V D,ss directly from structure were evaluated using a large set of clinical compounds. Machine learning (ML) models were built to predict V D,ss directly and to predict input parameters required for mechanistic and empirical V D,ss predictions. In addition, log D, fraction unbound in plasma (fup), and blood-to-plasma partition ratio (BPR) were measured on 254 compounds to estimate the impact of measured data on predictive performance of mechanistic models. Furthermore, the impact of novel methodologies such as measuring partition (Kp) in adipocytes and myocytes (n = 189) on V D,ss predictions was also investigated. In predicting V D,ss directly from chemical structures, both mechanistic and empirical scaling using a combination of predicted rat and dog V D,ss demonstrated comparable performance (62%–71% within 3-fold). The direct ML model outperformed other in silico methods (75% within 3-fold, r 2 = 0.5, AAFE = 2.2) when built from a larger data set. Scaling to human from predicted V D,ss of either rat or dog yielded poor results (<47% within 3-fold). Measured fup and BPR improved performance of mechanistic V D,ss predictions significantly (81% within 3-fold, r 2 = 0.6, AAFE = 2.0). Adipocyte intracellular Kp showed good correlation to the V D,ss but was limited in estimating the compounds with low V D,ss .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

AmeriFlux US-Cst Crossett Experimental Forest

This is the AmeriFlux version of the carbon flux data for the site US-Cst Crossett Experimental Forest. Site Description - The study takes place in the Crossett Experimental Forest (CEF, 32°2′ N, 91°57′ W), which was established in 1934 with the objective of developing effective management protocols for loblolly and shortleaf pine forests in the Gulf Coastal Plain. Mean annual temperature and precipitation are 17.6 °C and 1,410 mm, respectively. Most of the soils on and near the CEF are silt loams (primarily Glossaquic Fragiudalfs). The tower is situated in the the old-growth pine management sector of the CEF, which emphasizes management for old-growth-like conditions, with the goal of emulating the structural and compositional attributes of historic forests from this region using a combination of timber harvesting, prescribed fire, and other tools. This includes periodic commercial logging of intermediate-sized loblolly pines to support further restoration treatments, such as prescribed fire, hardwood and exotic species competition control, and the increase of plant species more compatible with fire (e.g., shortleaf pine and understory grasses)

Novick, Kim↗
DOE OSTI2024

AmeriFlux FLUXNET-1F US-Cst Crossett Experimental Forest

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Cst Crossett Experimental Forest. This is the FLUXNET version of the carbon flux data for the site US-Cst Crossett Experimental Forest produced by applying the standard ONEFlux (1F) software. Site Description - The study takes place in the Crossett Experimental Forest (CEF, 32°2′ N, 91°57′ W), which was established in 1934 with the objective of developing effective management protocols for loblolly and shortleaf pine forests in the Gulf Coastal Plain. Mean annual temperature and precipitation are 17.6 °C and 1,410 mm, respectively. Most of the soils on and near the CEF are silt loams (primarily Glossaquic Fragiudalfs). The tower is situated in the the old-growth pine management sector of the CEF, which emphasizes management for old-growth-like conditions, with the goal of emulating the structural and compositional attributes of historic forests from this region using a combination of timber harvesting, prescribed fire, and other tools. This includes periodic commercial logging of intermediate-sized loblolly pines to support further restoration treatments, such as prescribed fire, hardwood and exotic species competition control, and the increase of plant species more compatible with fire (e.g., shortleaf pine and understory grasses)

Novick, Kim↗
DOE OSTI2026

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gigaton commercial-scale carbon storage and mineralization potential in stacked Columbia River basalt reservoirs

This work presents a detailed supercritical CO 2 storage resource estimation for the stacked basalt reservoirs in the Grande Ronde Basalt of the Columbia River Basalt Group in eastern Washington and Oregon. The assessment aims to derisk the commercialization potential of geologic carbon storage in basalt by leveraging both structural and mineralization trapping of CO 2 in basalt. The structural closures formed by anticlinal ridges and synclinal valleys in Yakima Fold Belt are excellent physical traps to accommodate injected supercritical CO 2 . Rigorous hydraulic testing, well logs and simulation results from the Wallula Basalt Pilot #1 well showed the occurrence of 17 suitable permeable injection zones (up to 2,496 mD) intercalated with dense seals (~2.6E-10 mD) in the Grand Ronde Basalt. In addition, geochemical studies showed fast reactions between supercritical CO 2 and dissolved basalt minerals to form stable carbonates. In conclusion, our calculation indicates up to 40 gigatons (P90) of mineralization storage resources exist in the Grande Ronde Basalt reservoirs.

58 GEOSCIENCES↗
DOE OSTI2025

Quantum chemically calculated Abraham parameters for quantifying and predicting polymer hydrophobicity

The leakage and accumulation of plastic in the environment is a significant and growing problem with numerous detrimental impacts and has led to a push toward the design and development of more environmentally benign materials. To this end, we have developed a quantum chemistry-based model for predicting the mobility of polymer materials from molecular structure. Hydrophobicity is used as a surrogate for mobility given that hydrophobic interactions drive much of the partitioning of contaminants in and out of various environmentally relevant compartments. To model polymer hydrophobicity, we adjusted a previously developed Quantum Chemically Calculated Abraham Parameter model to calculate Abraham parameters of small molecules from molecular structure information. The resulting model predicted the octanol-water partition coefficient (K OW ) of polymer repeating units with a root mean square error (RMSE) of 0.48 (log scale). Additionally, the hydrophobicity of high molecular weight polymer materials was captured through solubility parameters and Nile red staining experiments from the literature and predicted with RMSEs of 1.21 (J/cc) 0.5 and 3.42 nm, respectively. Finally, to test the environmental applicability of the model, the relative adsorption capacity of three polymers was predicted and used to unify sorption isotherms across multiple sorbates and polymer sorbents.

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

A strong blend in the morning: studying the circumgalactic medium before cosmic noon with strong, blended Ly α forest systems

ABSTRACT We study of the properties of a new class of circumgalactic medium absorbers identified in the Ly α forest: ‘Strong, Blended Lyman-α’ (or SBLA) absorption systems. We study SBLAs at 2.4 < z < 3.1 in SDSS-IV/eBOSS spectra by their strong extended Ly α absorption complexes covering 138 $\, \, {\rm km}\, {\rm s}^{-1}$ with an integrated $\log (N_{\rm H\, {\small I}}/\mathrm{cm}^{-2}) =16.04$$\substack{+0.05 \\ -0.06}$ and Doppler parameter b = 18.1$\substack{+0.7 \\ -0.4}$$\, \, {\rm km}\, {\rm s}^{-1}$. Clustering with the Ly α forest provides a large-scale structure bias of b = 2.34 ± 0.06 and halo mass estimate of $M_h \approx 10^{12}\, h^{-1}\, {\rm M_{\odot }}$ for our SBLA sample. We measure the ensemble mean column densities of 22 metal features in the SBLA composite spectrum and find that no single-population multiphase model for them is viable. We therefore explore the underlying SBLA population by forward modelling the SBLA absorption distribution. Based on covariance measurements and favoured populations we find that ≈25 per cent of our SBLAs have stronger metals. Using silicon only we find that our strong metal SBLAs trace gas with a log (nH/cm−3) > −2.40 for T = 103.5 K and show gas clumping on <210 parsec scales. We fit multiphase models to this strong subpopulation and find a low ionization phase with nH = 1 cm−3, T = 103.5 K, and [X/H] = 0.8, an intermediate ionization phase with log (nH/cm−3) = −3.05, T = 103.5 K and [X/H] = −0.8, and a poorly constrained higher ionization phase. We find that the low ionization phase favours cold, dense super-solar metallicity gas with a clumping scale of just 0.009 parsecs.

Morrison, Sean (ORCID:0000000267702627)↗
DOE OSTI2024

Synthesizing realistic sand assemblies with denoising diffusion in latent space

Abstract The shapes and morphological features of grains in sand assemblies have far‐reaching implications in many engineering applications, such as geotechnical engineering, computer animations, petroleum engineering, and concentrated solar power. Yet, our understanding of the influence of grain geometries on macroscopic response is often only qualitative, due to the limited availability of high‐quality 3D grain geometry data. In this paper, we introduce a denoising diffusion algorithm that uses a set of point clouds collected from the surface of individual sand grains to generate grains in the latent space. By employing a point cloud autoencoder, the three‐dimensional point cloud structures of sand grains are first encoded into a lower‐dimensional latent space. A generative denoising diffusion probabilistic model is trained to produce synthetic sand that maximizes the log‐likelihood of the generated samples belonging to the original data distribution measured by a Kullback‐Leibler divergence. Numerical experiments suggest that the proposed method is capable of generating realistic grains with morphology, shapes and sizes consistent with the training data inferred from an F50 sand database. We then use a rigid contact dynamic simulator to pour the synthetic sand in a confined volume to form granular assemblies in a static equilibrium state with targeted distribution properties. To ensure third‐party validation, 50,000 synthetic sand grains and the 1542 real synchrotron microcomputed tomography (SMT) scans of the F50 sand, as well as the granular assemblies composed of synthetic sand grains are made available in an open‐source repository.

Vlassis, Nikolaos N.↗
DOE OSTI2025

Classical eikonal from Magnus expansion

In a classical scattering problem, the classical eikonal is defined as the generator of the canonical transformation that maps in-states to out-states. It can be regarded as the classical limit of the log of the quantum S-matrix. In a classical analog of the Born approximation in quantum mechanics, the classical eikonal admits an expansion in oriented tree graphs, where oriented edges denote retarded/advanced worldline propagators. The Magnus expansion, which takes the log of a time-ordered exponential integral, offers an efficient method to compute the coefficients of the tree graphs to all orders. We exploit a Hopf algebra structure behind the Magnus expansion to develop a fast algorithm which can compute the tree coefficients up to the 12th order (over half a million trees) in less than an hour. In a relativistic setting, our methods can be applied to the post-Minkowskian (PM) expansion for gravitational binaries in the worldline formalism. We demonstrate the methods by computing the 3PM eikonal and find agreement with previous results based on amplitude methods. Importantly, the Magnus expansion yields a finite eikonal, while the naïve eikonal based on the time-symmetric propagator is infrared-divergent from 3PM on.

Black Holes↗
DOE OSTI2022

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

Spatially resolved properties of early-type group-dominant galaxies with MUSE: gas content, ionization mechanisms, and metallicity gradients

ABSTRACT With the goal of a thorough investigation of the ionized gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). These data allowed us to study the spatially resolved warm gas properties, including the morphology of the ionized gas, EW(H α), and kinematics as well as the gas-phase metallicity (12 + log(O/H)) of these systems. In order to distinguish between different ionization mechanisms, we used the emission-line ratios [O iii]/H β and [N ii]/H α in the BPT diagrams and EW(H α). We find that the ionization sources in our sample have variable impacts at different radii; central regions are more influenced by low-luminosity active galactic nuclei, while extended regions of low-ionization nuclear emission-line region-like emission are ionized by other mechanisms with post-asymptotic giant branch stars photoionization likely contributing significantly. We classified our sample into three H α +[N ii] emission morphology types. We calculate the gas-phase metallicity assuming several methods and ionization sources. In general, 12 + log(O/H) decreases with radius from the centre for all galaxies, independently of nebular morphology type, indicating a metallicity gradient in the abundance profiles. Interestingly, the more extended filamentary structures and all extranuclear star-forming regions present shallow metallicity gradients. Within the uncertainties these extended structures can be considered chemically homogeneous. We suggest that group-dominant galaxies in our sample likely acquired their cold gas in the past as a consequence of one or more mechanisms, e.g. gas-clouds or satellite mergers/accretion and/or cooling flows that contribute to the growth of the ionized gas structures.

Lagos, P. (ORCID:0000000223218657)↗
DOE OSTI2023

The effect of composition, temperature and pressure on the oxidation state and coordination environment of copper in silicate melts

Copper is a redox variable element that may occur as Cu 0 , Cu + , and Cu 2+ in the Earth’s crust. The oxidation state will affect its partitioning between coexisting minerals, melts and fluids and hence its behaviour in magmatic processes. Copper bearing silicate glasses were quenched from melts with 23 synthetic compositions (19 CaO-Na 2 O-MgO-Al 2 O 3 -SiO 2 (CNMAS), two “granites” containing K 2 O ± H 2 O, and Fe-bearing “MORB” and “andesite”) equilibrated at oxygen fugacities (fO 2 ), expressed in log units relative to the fayalite-magnetite-quartz (FMQ) buffer, ranging from −0.7 to 14, temperatures from 900 to 1500 °C and pressures from 0 to 2.5 GPa. Cu K-edge X-ray absorption near edge structure (XANES) spectra were recorded from the glasses and a pre-edge feature in the XANES spectra was found to scale with the proportion of Cu + . Cu + /ΣCu (where ΣCu = Cu + + Cu 2+ ) was quantified by fitting the intensity of the pre-edge feature as a function of fO 2 to the thermodynamically expected relationship. Cu + /ΣCu was found to only weakly depend on melt composition, with more basic melts (e.g., basalts rather than granites) preferentially stabilising Cu 2+ . Increasing temperature stabilises Cu + , while increasing pressure had little effect on Cu + /ΣCu in CNMAS melts but preferentially stabilised Cu 2+ in granite melts. Cu + /ΣCu can be predicted in silicate melts by the empirical equation: log(Cu 2+ /Cu + ) = 0.25(ΔFMQ + 8.58 − 25050/T + 940P/T – 0.02P) − 4.73 + 5400/T + 1.99Λ + (280P – 90P2) where T is temperature in K, P is pressure in GPa and is the optical basicity of the composition. The effects of fO 2 , melt composition, temperature and pressure on Cu + /ΣCu indicate that Cu + will be the dominant oxidation state in terrestrial silicate melts (e.g., Cu + /ΣCu = 99% in an andesitic melt at 900 °C, 1 GPa and ΔFMQ = 1). The electron exchange reaction Cu 2+ + Fe 2+ → Cu + + Fe 3+ occurs on cooling and, given the abundance of Fe in natural melts, the oxidation state of Cu in natural glasses is unlikely to correspond to that in the original melt.

Copper↗
DOE OSTI2024

A high internal heat flux and large core in a warm Neptune exoplanet

Interactions between exoplanetary atmospheres and internal properties have long been proposed to be drivers of the inflation mechanisms of gaseous planets and apparent atmospheric chemical disequilibrium conditions. However, transmission spectra of exoplanets have been limited in their ability to observationally confirm these theories owing to the limited wavelength coverage of the Hubble Space Telescope (HST) and inferences of single molecules, mostly H 2 O. Here, in this work, we present the panchromatic transmission spectrum of the approximately 750 K, low-density, Neptune-sized exoplanet WASP-107b using a combination of HST Wide Field Camera 3 (WFC3) and JWST Near-Infrared Camera (NIRCam) and Mid-Infrared Instrument (MIRI). From this spectrum, we detect spectroscopic features resulting from H 2 O (21σ), CH 4 (5σ), CO (7σ), CO 2 (29σ), SO 2 (9σ) and NH 3 (6σ). The presence of these molecules enables constraints on the atmospheric metal enrichment (M/H is 10–18× solar), vertical mixing strength (log 10 K zz = 8.4–9.0 cm 2 s -1 ) and internal temperature (>345 K). The high internal temperature is suggestive of tidally driven inflation4 acting on a Neptune-like internal structure, which can naturally explain the large radius and low density of the planet. These findings suggest that eccentricity-driven tidal heating is a critical process governing atmospheric chemistry and interior-structure inferences for most of the cool (<1,000 K) super-Earth-to-Saturn-mass exoplanet population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Cosmic Ultraviolet Baryon Survey (CUBS) – III. Physical properties and elemental abundances of Lyman-limit systems at z < 1

ABSTRACT We present a systematic investigation of physical conditions and elemental abundances in four optically thick Lyman-limit systems (LLSs) at z = 0.36–0.6 discovered within the Cosmic Ultraviolet Baryon Survey (CUBS). Because intervening LLSs at z < 1 suppress far-UV (ultraviolet) light from background QSOs, an unbiased search of these absorbers requires a near-UV-selected QSO sample, as achieved by CUBS. CUBS LLSs exhibit multicomponent kinematic structure and a complex mix of multiphase gas, with associated metal transitions from multiple ionization states such as C ii, C iii, N iii, Mg ii, Si ii, Si iii, O ii, O iii, O vi, and Fe ii absorption that span several hundred km s−1 in line-of-sight velocity. Specifically, higher column density components (log N(H i)/cm−2≳ 16) in all four absorbers comprise dynamically cool gas with $\langle T \rangle =(2\pm 1) \times 10^4\,$K and modest non-thermal broadening of $\langle b_\mathrm{nt} \rangle =5\pm 3\,$km s−1. The high quality of the QSO absorption spectra allows us to infer the physical conditions of the gas, using a detailed ionization modelling that takes into account the resolved component structures of H i and metal transitions. The range of inferred gas densities indicates that these absorbers consist of spatially compact clouds with a median line-of-sight thickness of $160^{+140}_{-50}$ pc. While obtaining robust metallicity constraints for the low density, highly ionized phase remains challenging due to the uncertain $N\mathrm{(H\, {\small I})}$, we demonstrate that the cool-phase gas in LLSs has a median metallicity of $\mathrm{[\alpha /H]_{1/2}}=-0.7^{+0.1}_{-0.2}$, with a 16–84 percentile range of [α/H] = (−1.3, −0.1). Furthermore, the wide range of inferred elemental abundance ratios ([C/α], [N/α], and [Fe/α]) indicate a diversity of chemical enrichment histories. Combining the absorption data with deep galaxy survey data characterizing the galaxy environment of these absorbers, we discuss the physical connection between star-forming regions in galaxies and diffuse gas associated with optically thick absorption systems in the z < 1 circumgalactic medium.

Zahedy, Fakhri S. (ORCID:0000000178692551)↗