Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pair distribution functions”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2020

Grid-Following Inverters and Synchronous Condensers: A Grid-Forming Pair?: Preprint

As power systems across the globe continue to achieve higher instantaneous penetrations of power electronic converterinterfaced renewable sources, the stability of these power systems is challenged because of the removal of synchronous machines, from which stability is traditionally derived. Although technological solutions to these stability challenges are on the horizon, such as the use of grid-forming inverters, they are not yet widely applied to larger power systems, which presents operational challenges for power systems achieving these high instantaneous penetrations today. An interim solution using existing technologies is to pair synchronous condensers with grid-following inverters, which might prolong the stability of an operating power system while synchronous generators are turned off during periods of high renewable energy availability. This work examines the transient stability of such a solution using PSCAD simulations of a two-bus system coupled to a variedlength transmission line, with a synchronous condenser at one bus and a grid-following inverter with grid-support functionality at the other. The system is exposed to load step, balanced fault, and unbalanced fault perturbations. It was found that this simple system returned to a steady state for transmission line lengths up to 125 km after 10% load steps, or a variety of fault types.

14 SOLAR ENERGY↗
DOE OSTI2024

Enhancing Cloud Cybersecurity: Prescriptive Controls for Operational Technology

This whitepaper provides strategic insights and recommendations into security cloud-based solutions for electric utilities, encompassing operational technology (OT), virtual power plants (VPP), distributed energy resources (DERs), applications, networks, and data storage as they transition to and leverage cloud infrastructure through managed service providers (MSPs) and cloud service providers (CSPs). Principles derived from established frameworks serve as a foundation for best practices across cybersecurity projects and remove the constraints of settling on a single framework. For organizations that prefer not to integrate a specific framework altogether, elements of the proposed approach could be adopted or tailored to best fit defined requirements and expected functionalities. The Cirrus assessment, a utility cloud feasibility tool, and the roadmap it provides serve as a precursor to this paper, which seeks to be a valuable resource for defining next steps following cloud technology integration feasibility appraisal. With its comprehensive approach to adoption, the Cirrus framework offers strategic guidance on responsibly preparing for or deploying a utility cloud solution. The previously published whitepaper, “Use Case-Informed Framework for Utility Cloud Migration,” details the guiding strategy, research, and deployment of cloud solutions within electric and interconnected grid systems. Before implementing the controls suggested in this document, it is recommended that stakeholders complete Cirrus's cloud integration assessment and pair the results with their unique cybersecurity controls to form a comprehensive cloud-based utility cybersecurity plan. The Cirrus outcome will consider a series of future architectures for the grid before and after the energy transition and evaluate the arguments for and against cloud applications for each electric and interconnected grid layer. This document is a companion to the original whitepaper, "Use Case-Informed Framework for Utility Cloud Migration" to further identify and recommend security controls based on Cirrus’s cloud integration assessment output. The following whitepaper outlines the cybersecurity controls that secure cloud-service models pertinent to the electric sector using the predefined categories identify, protect, detect, and respond and recover. The objective is to outline prescriptive security controls based on the type of architecture and data stored in the cloud. The focus includes dissecting the shared responsibility model and elucidating what on-premises Infrastructure as a Service (IaaS), Platform as a Service (PaaS), and Software as a Service (SaaS) entail. A pivotal consideration in this context is allocating responsibility for foundational cybersecurity aspects—having used Cirrus for the cloud integration assessment. The ensuing controls detailed herein also represent a checklist of controls necessary for a secure cloud transition, equipping utilities with the knowledge to navigate this digital transformation with confidence and strategic foresight in a safe and responsible manner.

42 ENGINEERING↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2020

Longitudinal Wobbling Motion in Au 187

The rare phenomenon of nuclear wobbling motion has been investigated in the nucleus 187Au. A longitudinal wobbling-bands pair has been identified and clearly distinguished from the associated signature-partner band on the basis of angular distribution measurements. Theoretical calculations in the framework of the particle rotor model are found to agree well with the experimental observations. This is the first experimental evidence for longitudinal wobbling bands where the expected signature partner band has also been identified, and establishes this exotic collective mode as a general phenomenon over the nuclear chart.

150 ≤ A ≤ 189↗
DOE OSTI2026

Precise measurement of the tt¯ production cross-section and lepton differential distributions in eμ dilepton events from s=13TeVpp collisions with the ATLAS detector

The inclusive top quark pair (tt¯$$t\bar{t}$$) cross-section σtt¯$$\sigma _{t\bar{t}}$$ has been measured in proton–proton collisions at s=13TeV$$\sqrt{s}=13\,\text {TeV}$$, using 140fb-1$$140\,{\text {fb}^{-1}} $$ of data collected by the ATLAS experiment at the Large Hadron Collider. Using events with an opposite-charge eμ$$e\mu $$ pair and b-tagged jets, the cross-section is measured to be: σtt¯=829.3±1.3(stat)±8.0(syst)±7.3(lumi)±1.9(beam)pb,$$\begin{aligned} \sigma _{t\bar{t}} & = 829.3\pm 1.3\,\mathrm {(stat)}\ \pm 8.0\,\mathrm {(syst)}\ \pm 7.3\,\mathrm {(lumi)}\ \\ & \quad \pm 1.9\,\mathrm {(beam)}\,\textrm{pb}, \end{aligned}$$where the uncertainties reflect the limited size of the data sample, experimental and theoretical systematic effects, the integrated luminosity, and the proton beam energy, giving a total uncertainty of 1.3%. The result is used to determine the top quark pole mass via the dependence of the predicted cross-section on mtpole$${m_{t}^\textrm{pole}}$$, giving mtpole=172.8-1.7+1.5$${m_{t}^\textrm{pole}}=172.8^{+1.5}_{-1.7}$$ GeV$$\text {GeV}$$. The same event sample is used to measure absolute and normalised differential cross-sections for the tt¯→eμνν¯bb¯$$t\bar{t} \rightarrow e\mu u \bar{ u }b\bar{b} $$ process as a function of single-lepton and dilepton kinematic variables. Complementary measurements of eμbb¯$$e\mu b\bar{b} $$ production, treating both tt¯$$t\bar{t}$$ and Wt events as signal, are also provided. Both sets of differential cross-sections are compared to the predictions of various Monte Carlo event generators, demonstrating that the state-of-the-art generators Powheg MiNNLO and Powhegbb4l describe the data better than Powheghvq.

Aad, G↗
DOE OSTI2020

Grid-Following Inverters and Synchronous Condensers: A Grid-Forming Pair?

As power systems across the globe continue to achieve higher instantaneous penetrations of power electronic converter-interfaced renewable sources, the stability of these power systems is challenged because of the removal of synchronous machines, from which stability is traditionally derived. Although technological solutions to these stability challenges are on the horizon, such as the use of grid-forming inverters, they are not yet widely applied to larger power systems, which presents operational challenges for power systems achieving these high instantaneous penetrations today. A potential interim solution using existing technologies is to pair synchronous condensers with grid-following inverters, which might prolong the stability of an operating power system while synchronous generators are turned off during periods of high renewable energy availability. This work examines the transient stability of such a solution using PSCAD simulations of a two-bus system coupled to a varied-length transmission line, with a synchronous condenser at one bus and a grid-following inverter with grid-support functionality at the other. The system is exposed to load step, balanced fault, and unbalanced fault perturbations. It was found that this simple system returned to a steady state for transmission line lengths up to 125 km after 10% load steps, or a variety of fault types.

grid-following inverters↗
DOE OSTI2021

Hyperheavy spherical and toroidal nuclei: The role of shell structure

The properties of toroidal hyperheavy even-even nuclei and the role of toroidal shell structure are extensively studied within covariant density functional theory. The general trends in the evolution of toroidal shapes in the Z ≈ 130–180 region of the nuclear chart are established for the first time. These nuclei are stable with respect to breathing deformations. The most compact fat toroidal nuclei are located in the Z ≈ 136, N ≈ 206 region of the nuclear chart, but thin toroidal nuclei become dominant with increasing proton number and toward proton and neutron drip lines. The roles of toroidal shell structure, its regularity, supershell structure, and shell gaps as well as the role of different groups of the pairs of the orbitals in its formation are investigated in detail. The lowest in energy solutions at axial symmetry are characterized either by large shell gaps or low density of the single-particle states in the vicinity of the Fermi level in at least one of the subsystems (proton or neutron). Related quantum shell effects are expected to act against the instabilities in breathing and sausage deformations for these subsystems. The investigation with a large set of covariant energy density functionals reveals that substantial proton Z = 154 and 186 and neutron N = 228, 308, and 406 spherical shell gaps exist in all functionals. The nuclei in the vicinity of the combination of these particle numbers form the islands of stability of spherical hyperheavy nuclei. Furthermore, the study suggests that the N = 210 toroidal shell gap plays a substantial role in the stabilization of fat toroidal nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Energy dependence of coherent photonuclear production of J/ψ mesons in ultra-peripheral Pb-Pb collisions at $\sqrt{{\textrm{s}}_{\textrm{NN}}}$ = 5.02 TeV

The cross section for coherent photonuclear production of J/ψ is presented as a function of the electromagnetic dissociation (EMD) of Pb. The measurement is performed with the ALICE detector in ultra-peripheral Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{{\textrm{s}}_{\textrm{NN}}}$ = 5.02 TeV. Cross sections are presented in five different J/ψ rapidity ranges within |y| < 4, with the J/ψ reconstructed via its dilepton decay channels. In some events the J/ψ is not accompanied by EMD, while other events do produce neutrons from EMD at beam rapidities either in one or the other beam direction, or in both. The cross sections in a given rapidity range and for different configurations of neutrons from EMD allow for the extraction of the energy dependence of this process in the range 17 < W γ Pb,n < 920 GeV, where W γ Pb,n is the centre-of-mass energy per nucleon of the γPb system. This range corresponds to a Bjorken-x interval spanning about three orders of magnitude: 1.1 × 10 —5 < x < 3.3 × 10 —2 . In addition to the ultra-peripheral and photonuclear cross sections, the nuclear suppression factor is obtained. These measurements point to a strong depletion of the gluon distribution in Pb nuclei over a broad, previously unexplored, energy range. These results, together with previous ALICE measurements, provide unprecedented information to probe quantum chromodynamics at high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Differential charge radii: Proton-neutron interaction effects

The analysis of self-consistency and proton-neutron interaction effects in the buildup of differential charge radii were carried out in covariant density functional theoretical calculations without pairing interaction on the example of selected configurations in the Pb isotopes. The proton-neutron interaction of neutron(s) added to the neutron 𝑁 = 126 core and the protons forming the 𝑍 = 82 proton core is responsible for a major contribution to the buildup of differential charge radii. It depends on the products of proton and neutron wave functions and thus on their nodal structure. This interaction leads to a redistribution of single-particle density of occupied proton states which in turn modifies the charge radii. Furthermore, the microscopic origin of this redistribution and its consequences for differential charge radii were investigated for the first time. Self-consistency effects affecting the shape of proton potential, total proton densities and the energies of the single-particle proton states provide only minor contribution to differential charge radii.

190 ≤ A ≤ 219↗
DOE OSTI2024

Helium-4 gravitational form factors: Exchange currents

We evaluate the leading exchange corrections to the helium-4 gravitational form factors (GFFs) to momenta of the order of the nucleon mass. We use both the K-harmonic method with simple pair nucleon potential, and a Jastrow trial function using the Argonne 𝑣 14 potential, to evaluate the helium-4 GFFs. The exchange current contributions include the pair interaction, plus the seagull and the pion exchange interactions, modulo the recoil corrections. To estimate the off-shellness of the pion nucleon coupling in this momenta range, we discuss the results using either the pseudoscalar (PS) or pseudovector (PV) pion-nucleon couplings. When the PV coupling is used, the pair diagram contribution is higher order in the nonrelativistic expansion. The results for the helium-4 A-GFF are comparable to those given by the impulse approximation, especially for the PS coupling using both the K-harmonic method and variational method. The exchange current contributions with the PS coupling for the charge form factor of helium-4, yield better agreement with the existing data over a broad range of momenta, especially when the Argonne 𝑣 14 potential including the D-wave admixture is used.

A ≤ 5↗
DOE OSTI2022

Exclusive production of a pair of high transverse momentum photons in pion-nucleon collisions for extracting generalized parton distributions

We show that exclusive production of a pair of high transverse momentum photons in pion-nucleon collisions can be systematically studied in QCD factorization approach if the photon’s transverse momentum q T with respect to the colliding pion is much greater than Λ QCD . We demonstrate that the leading power non-perturbative contributions to the scattering amplitudes of this exclusive process are process-independent and can be systematically factorized into universal pion’s distribution amplitudes (DAs) and nucleon’s generalized parton distributions (GPDs), which are convoluted with corresponding infrared safe and perturbatively calculable short-distance hard parts. The correction to this factorized expression is suppressed by powers of 1/q T . We demonstrate quantitatively that this new type of exclusive processes is not only complementary to existing processes for extracting GPDs, but also capable of providing an enhanced sensitivity to the dependence of both DAs and GPDs on the active parton’s momentum fraction x. We also introduce additional, but the same type of exclusive observables to enhance our capability to explore GPDs, in particular, their x-dependence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Investigation and Expression Analysis of R2R3-MYBs and Anthocyanin Biosynthesis-Related Genes during Seed Color Development of Common Bean ( Phaseolus vulgaris )

Anthocyanins are responsible for the coloration of common bean seeds, and their accumulation is positively correlated with the expression level of anthocyanin biosynthetic genes. The MBW (MYB-bHLH-WD40) complex is thought to regulate the expression of these genes, and MYB proteins, which are a key factor in activating anthocyanin pathway genes, have been identified in several plants. This study demonstrated gene structures, chromosomal placements, gene duplications of R2R3-MYBs, miRNAs associated with R2R3-MYBs, and the interaction of these genes with other flavonoid regulatory genes. qRT-PCR was used to investigate the role of specific R2R3-MYBs and flavonoid genes in common bean seed color development. As a result of a comprehensive analysis with the help of in silico tools, we identified 160 R2R3-MYB genes in the common bean genome. We divided these genes into 16 classes on the basis of their intron-exon and motif structures. Except for three, the rest of the common bean R2R3-MYB members were distributed to all chromosomes with different densities, primarily located on chromosomes 3 and 8. We identified a total of 44 duplicated gene pairs dispersed across 11 chromosomes and evolved under purifying selection (Ka/Ks < 1), 19 of which were derived from a whole-genome duplication. Our research uncovered 25 putative repressor PvMYB proteins that contain the EAR motif. Additionally, fifty different cis-regulatory elements regulated by light, stress, and hormone were identified. Within the genome of the common bean, we discovered a total of 36 microRNAs that target a total of 72 R2R3-MYB transcripts. The effect of 16 R2R3-MYB genes and 16 phenylpropanoid pathway genes, selected on the basis of their interaction in the protein-protein interaction map, playing role in the regulation of seed coat color development was evaluated using qRT-PCR in 5 different tissues at different developmental stages. The results revealed that these specific genes have different expression levels during different developmental periods, with higher levels in the pod filling and early pod stages than in the rest of the developmental periods. Furthermore, it was shown that PvTT8 (bHLH), PvTT2 (PvMYB42), PvMYB113, PvTTG1, and PvWD68 genes have effects on the regulation of seed coat color. The findings of this study, which is the first to use whole-genome analysis to identify and characterize the R2R3-MYB genes in common bean, may serve as a reference for future functional research in the legume.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Ab initio electron dynamics in high electric fields: Accurate prediction of velocity-field curves

Electron dynamics in external electric fields governs the behavior of solid-state electronic devices. First-principles calculations enable precise predictions of charge transport in low electric fields. However, studies of high-field electron dynamics remain elusive due to a lack of accurate and broadly applicable methods. In this work, we develop an efficient approach to solve the real-time Boltzmann transport equation with both the electric field term and ab initio electron-phonon collisions. These simulations provide field-dependent electronic distributions in the time domain, allowing us to investigate both transient and steady-state transport in electric fields ranging from low to high (> 10 kV/cm). The broad capabilities of our approach are shown by computing nonequilibrium electron occupations and velocity-field curves in Si, GaAs, and graphene, obtaining results in quantitative agreement with experiment. Our approach sheds light on microscopic details of transport in high electric fields, including the dominant scattering mechanisms and valley occupation dynamics. Our results demonstrate quantitatively accurate calculations of electron dynamics in low to high electric fields, with broad application to power and micro-electronics, optoelectronics, and sensing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Exploring the Nucleon Structure via Deep Electroproduction Processes

Understanding the internal structure of the nucleon is a fundamental goal of modern physics, which aims at a comprehensive framework describing the internal dynamics of quarks and gluons. Among other structure functions, Generalized Parton Distributions (GPDs) offer a powerful framework for describing the nucleon dynamics by correlating the longitudinal momentum and the transverse position of its internal partons. Such a correlation provides a three-dimensional picture of the nucleon and enables access to fundamental properties, including the internal pressure distributions and the parton’s angular momentum contribution to the nucleon's total spin, thereby playing a central role in resolving the nucleon spin puzzle. At the Thomas Jefferson National Facility (JLab), polarized electron beam experiments allow for probing GPDs through the measurement of hard exclusive processes. Among the cleanest experimental channels, we find the electro-production of a real photon through the Deeply Virtual Compton Scattering (DVCS) mechanism. The first data-taking period of the CLAS12 program, taking place in 2018, allowed for unique DVCS Beam Spin Asymmetry (BSA) measurements in the phase space covered by a 10.6 GeV polarized electron beam impinging on an unpolarized liquid hydrogen target. Although detecting all final-state particles ensures exclusivity of the process, conservation laws indicate that it is not mandatory. I adopt an approach omitting the direct detection of the recoil proton, providing a simplified yet effective event selection strategy that boosts statistics and gives access to a larger phase space sensitive to the underlying GPD dynamics through BSA and cross section measurements. The Double DVCS (DDVCS) process promises a dedicated mapping of GPDs. Contributing to the electro-production of a lepton pair cross-section, the DDVCS reaction extends DVCS by allowing the final-state photon to be virtual, enriching the kinematic phase space and providing unique access to the internal correlations encoded by GPDs. A feasibility study is conducted to assess the potential of future DDVCS measurements at Jefferson Lab and the future Electron-Ion Collider (EIC). While Jefferson Lab will provide DDVCS measurements in the valence region through the SoLID$\mu$ and $\mu$CLAS12 experimental projects, in the long term, the EIC will provide complementary measurements in the sea region, both accessing unprecedented information about GPDs in a phase space region otherwise inaccessible. Taken together, these investigations demonstrate both the current capabilities and future opportunities for probing GPDs through exclusive processes. The experimental analysis of DVCS at CLAS12 provides precise measurements within an established framework, while the phenomenological study of DDVCS opens the door to richer and more comprehensive explorations with future detectors and facilities.

Alvarado, Juan [Université Paris-Saclay: Gif-sur-Y↗
DOE OSTI2023

Utah FORGE: Reports on Northwestern Nevada Well Doublet Drilling and Testing by Fervo Energy

These reports review Fervo Energy's construction of a commercial enhanced geothermal system (EGS). Fervo has qualified full functionality of the system through production testing at commercially relevant operating conditions. The project site is located in a nearfield setting adjacent to an operating geothermal power station in north-central Nevada and is designed to deliver an uplift in high-temperature geothermal flow rates to increase the power capacity at the facility. The project involved drilling a first-of-a-kind EGS horizontal doublet well system consisting of an injection and production well pair within a high-temperature hard-rock geothermal formation. The lithology of the target reservoir is characterized as a mixed metasedimentary and igneous formation, comprised of phyllite, quartzite, diorite, and granodiorite, representative of the geology across the most prospective geothermal areas throughout the western US. This work included a 16-stage plug-and-perforate stimulation treatment which was the first of its kind in a high-temperature mixed metasedimentary and granitic formation in a fully horizontal doublet geothermal well.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Charge Radius of Neutron-Deficient Ni 54 and Symmetry Energy Constraints Using the Difference in Mirror Pair Charge Radii

he nuclear root-mean-square charge radius of Ni 54 was determined with collinear laser spectroscopy to be R ( Ni 54 ) = 3.737 ( 3 ) fm . In conjunction with the known radius of the mirror nucleus Fe 54 , the difference of the charge radii was extracted as Δ R ch = 0.049 ( 4 ) fm . Based on the correlation between Δ R ch and the slope of the symmetry energy at nuclear saturation density ( L ), we deduced 21 ≤ L ≤ 88 MeV . The present result is consistent with the L from the binary neutron star merger GW170817, favoring a soft neutron matter EOS, and barely consistent with the PREX-2 result within 1 σ error bands. Our result indicates the neutron-skin thickness of Ca 48 as 0.15–0.21 fm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗