Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial correlation”

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

Moiré-driven topological electronic crystals in twisted graphene

In a dilute two-dimensional electron gas, Coulomb interactions can stabilize the formation of a Wigner crystal. Although Wigner crystals are topologically trivial, it has been predicted that electrons in a partially filled band can break continuous translational symmetry and time-reversal symmetry spontaneously, resulting in a type of topological electron crystal known as an anomalous Hall crystal. Here we report signatures of a generalized version of the anomalous Hall crystal in twisted bilayer–trilayer graphene, whose formation is driven by the moiré potential. The crystal forms at a band filling of one electron per four moiré unit cells (ν = 1/4) and quadruples the unit-cell area, coinciding with an integer quantum anomalous Hall effect. The Chern number of the state is exceptionally tunable, and it can be switched reversibly between +1 and −1 by electric and magnetic fields. Several other topological electronic crystals arise in a modest magnetic field, originating from ν = 1/3, 1/2, 2/3 and 3/2. As a result, the quantum geometry of the interaction-modified bands is likely to be very different from that of the original parent band, which enables possible future discoveries of correlation-driven topological phenomena.

Electronic properties and materials↗
DOE OSTI2025

Restructuring of Benzimidazole-Based Copper Complexes during Electrochemical CO 2 Reduction

Copper-based molecular catalysts have been widely used for the electrochemical reduction of carbon dioxide (ERCO 2 ). However, establishing a correlation between structure and catalytic performance has been challenging due to the dynamic nature of these materials under ERCO 2 working conditions. Here, we report on the restructuring of the copper complexes [Cu(bzimpy)Cl 2 ] and [Cu(pyrben) 2 (NO 3 )]NO 3 during ERCO 2 in an acidic medium, evidencing its effects on the observed activity, selectivity, and stability. In situ X-ray absorption spectroscopy (XAS) suggested that, at sufficiently negative potentials, the Cu 2+ centers of both complexes are reduced to undercoordinated Cu 0 species. Pb underpotential deposition of post-ERCO 2 samples indicated that such Cu 0 species also present a dominant contribution of the (100) domain, which is in line with the high Faradaic efficiency toward C 2 products (∼45%) shown by both complexes at −1.24 V vs RHE. By combining online electrochemical mass spectrometry, in situ XAS, electron paramagnetic resonance (EPR), and quantitative product analysis, we found that [Cu(bzimpy)Cl 2 ] exhibits a partially reversible restructuring, regenerating the main complex moiety, [Cu(bzimpy)] 2+ , and its geometry.

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

Extended quantum anomalous Hall states in graphene/hBN moiré superlattices

Electrons in topological flat bands can form new topological states driven by correlation effects. The pentalayer rhombohedral graphene/hexagonal boron nitride (hBN) moiré superlattice was shown to host fractional quantum anomalous Hall effect (FQAHE) at approximately 400 mK, triggering discussions around the underlying mechanism and role of moiré effects. In particular, new electron crystal states with non-trivial topology have been proposed. Here, in this paper, we report electrical transport measurements in rhombohedral pentalayer and tetralayer graphene/hBN moiré superlattices at electronic temperatures down to below 40 mK. We observed two more fractional quantum anomalous Hall (FQAH) states and smaller R xx values in pentalayer devices than those previously reported. In the new tetralayer device, we observed FQAHE at moiré filling factors v = 3/5 and 2/3. With a small current at the base temperature, we observed a new extended quantum anomalous Hall (EQAH) state and magnetic hysteresis, where R xy = h/e 2 and vanishing R xx spans a wide range of v from 0.5 to 1.3. At increased temperature or current, EQAH states disappear and partially transition into the FQAH liquid. Furthermore, we observed displacement field-induced quantum phase transitions from the EQAH states to the Fermi liquid, FQAH liquid and the likely composite Fermi liquid. Our observations established a new topological phase of electrons with quantized Hall resistance at zero magnetic field and enriched the emergent quantum phenomena in materials with topological flat bands.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Real-Frequency Response Functions at Finite Temperature

Building on previous developments [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 99, 035120 (2019); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B101, 125109 (2020); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B102, 045115 (2020), B. Holm and U. von Barth, Phys. Rev. B 57, 2108 (1998), J. Vičičevic and M. Ferrero, Phys. Rev. B 101, 075113 (2020)], we show that the diagrammatic Monte Carlo technique allows us to compute finite-temperature response functions directly on the real-frequency axis within any field-theoretical formulation of the interacting fermion problem. There are no limitations on the type and nature of the system’s action or whether partial summation and self-consistent treatment of certain diagram classes are used. In particular, by eliminating the need for numerical analytic continuation from a Matsubara representation, our scheme allows us to study spectral densities of arbitrary complexity with controlled accuracy in models with frequency-dependent effective interactions. Finally, for illustrative purposes we consider the problem of the plasmon linewidth in a homogeneous electron gas (jellium).

3-dimensional systems↗
DOE OSTI2020

TRANSENSOR: Transformer Real-time Assessment INtelligent System with Embedded Network of Sensors and Optical Readout. Final Report

Utilities across the world are wrestling with evolving market dynamics, population growth and climate change. Distributed Energy Resources (DERs) such as solar photovoltaics, distributed generators and energy storage systems are becoming important parts of the U.S. energy mix. These factors are driving an increasing need for low-cost grid asset monitoring. Aside from being costly, traditional utility monitoring systems are not sufficiently robust and do not provide real-time visibility into the condition of grid assets such as transformers. Lack of accurate real-time measurements on performance has resulted in the use of lagging indicators, such as oil sample analysis. To address this need, an innovative embedded optical sensing technology, Transformer Real-time Assessment Intelligent System (TRANSENSOR) was developed, validated and demonstrated in this project. To date, TRANSENSOR has focused on transformers but is extendable to other grid assets. In addition to the technology being embedded into new transformers during manufacturing, a retrofit configuration that can be installed on existing transformers in the field was also developed. It is anticipated that the ability to retrofit existing transformers will help accelerate adoption of the underlying TRANSENSOR technology. Phase 1 of the project focused on laboratory development and qualification of the technology. Following iterations and exploration of relevant optical sensing modalities and multiplexed configurations of interest for the transformer environment, an effective candidate configuration was agreed upon, down-selected and custom-designed for embedding into General Electric (GE) network transformers. Two new GE 500 kilovolt ampere (kVA), 27 kilovolt (kV) distribution network transformers were built with embedded fiber-optic (FO) sensors and successfully qualified per industry standards at GE’s Shreveport, Louisiana facility during Phase 1 of the project. Following the successful completion of Phase 1, the team proceeded to Phase 2, which focused on a field demonstration of the technology. Over Phase 2, the first new GE transformer built with embedded fiber-optic sensors was installed in an above-ground cage and connected to the grid at Con Edison’s Astoria facility. A second transformer equipped with TRANSENSOR was installed in an underground vault. Additionally, an older (1982 year model) GE transformer in an above-ground cage was retrofitted with fiber- optic sensors and reconnected to the grid at ConEd’s facility. Analysis of data acquired from the sensors showed interesting correlations with transformer loading. Additionally, key events such as transformer low-voltage network connection, primary-side energizing, and a pressure loss event from an oil sampling were detected by the TRANSENSOR system. Online data processing/feature extraction algorithms for the second transformer in the underground vault detected key features and event alerts that were transmitted through a wireless 4G connection. The remote deployment concept showed promising results with data collection running for a total of 8 months across the three (3) GE transformers instrumented for the Phase 2 field demonstration. This sets the stage well for further development and commercialization.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Inferring the Focal Depths of Small Earthquakes in Southern California Using Physics-Based Waveform Features

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ~8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < M L < 4) recorded at local distances (<150 km), for which depth estimation is especially challenging. We find that differential magnitudes (M w /M L –M c ) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡log 10 (M 0 )+3log 10 (f c ), and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. Here, we also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

58 GEOSCIENCES↗
DOE OSTI2025

Effect of Sulfonation Level on the Percolated Morphology and Proton Conductivity of Hydrated Fluorine-Free Copolymers: Experiments and Simulations

Using all-atom molecular dynamics simulations and a variety of experimental methods, we previously reported on a linear polyethylene with pendant phenyl sulfonated groups precisely on every fifth carbon along the backbone. With increasing relative humidity this fluorine-free polymer self-assembled to form nanoscale water channels and exhibited exceptional proton conductivity. Expanding upon those findings, here we explore partially sulfonated random copolymers, referred to as p 5PhSH-Y. Using either acetyl sulfate or sulfuric acid, a wide range of sulfonation levels were prepared ( Y = 34−98%) corresponding to ion-exchange capacities (IEC) of 2.0−4.4 mmol/g. Combining experimental techniques and all-atom molecular dynamics simulations, we study the effect of Y on water uptake, nanoscale morphology, and the proton/water transport properties of p5PhSH- Y . The proton conductivity of p 5PhSH- Y increases with relative humidity and with Y and achieves values in excess of 0.1 S/cm. These high conductivities are attributed to high IEC and welldeveloped nanoscale percolated hydrophilic domains made possible by the flexible backbone. We quantitatively describe the nature of the water channels using the characteristic distance, channel width distribution, the area per sulfonate group at the hydrophilic/ hydrophobic interface, and the fractal dimension. Notably, the channel widths and the areas per sulfonate group are nominally independent of the level of sulfonation, while depending significantly on the level of hydration. The fractal dimension of the water channels correlates strongly with the water diffusion coefficients calculated from the molecular dynamics (MD) simulations. These findings demonstrate that the p 5PhSH- Y hydrocarbon copolymers can be modified to tune properties, particularly proton conductivity.

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

Steering the methanol steam reforming reactivity of intermetallic Cu–In compounds by redox activation: stability vs . formation of an intermetallic compound–oxide interface

To compare the inherent methanol steam reforming properties of intermetallic compounds and a corresponding intermetallic compound–oxide interface, we selected the Cu–In system as a model to correlate the stability limits, self-activation and redox activation properties with the catalytic performance. Three distinct intermetallic Cu–In compounds – Cu 7 In 3 , Cu 2 In and Cu 11 In 9 – were studied both in an untreated and redox-activated state resulting from alternating oxidation–reduction cycles. The stability of all studied intermetallic compounds during methanol steam reforming (MSR) operation is essentially independent of the initial stoichiometry and all accordingly resist substantial structural changes. The inherent activity under batch MSR conditions is highest for Cu 2 In, corroborating the results of a Cu 2 In/In 2 O 3 sample accessed through reactive metal–support interaction. Under flow MSR operation, Cu 7 In 3 displays considerable deactivation, while Cu 2 In and Cu 11 In 9 feature stable performance at simultaneously high CO 2 selectivity. Here, the missing significant self-activation is most evident in the operando thermogravimetric experiments, where no oxidation is detected for any of the intermetallic compounds. In situ X-ray diffraction allowed us to monitor the partial decomposition and redox activation of the Cu–In intermetallic compounds into Cu 0.9 In 0.1 /In 2 O 3 (from Cu 7 In 3 ), Cu 7 In 3 /In 2 O 3 (from Cu 2 In) and Cu 7 In 3 /Cu 0.9 In 0.1 /In 2 O 3 (from Cu 11 In 9 ) interfaces with superior MSR performance compared to the untreated samples. Although the catalytic profiles appear surprisingly similar, the latter interface with the highest indium content exhibits the least deactivation, which we explain by formation of stabilizing In 2 O 3 patches under MSR conditions.

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

Microwave-Assisted Hydrogen Generation from Hydrocarbon-Bearing Reservoir Rocks: Stage-Dependent Thermal Runaway and In-Situ Carbonate Engineering

Microwave-assisted hydrogen generation from hydrocarbon-bearing reservoir rocks is strongly influenced by mineralogy, methane activation, carbonate reactions, and thermal runaway behavior. This study investigates a new approach in which carbonate phases are generated in-situ through the reaction of internally produced CO2 with Ca(OH)2 under microwave heating conditions. The objective is to evaluate how rock mineralogy, methane injection, and Ca(OH)2 addition influence hydrogen generation, carbon redistribution, and stage-dependent reaction pathways during microwave exposure. Microwave heating experiments were conducted using Permian Basin reservoir rocks under three experimental conditions: rock-only experiments under Ar atmosphere, CH4–Ar experiments without additive, and CH4–Ar experiments containing 5 wt% Ca(OH)2. Methane-assisted experiments were performed under continuous injection of 30 standard cubic centimeters per minute (sccm) CH4 and 30 sccm Ar. Based on thermal runaway behavior, each experiment was divided into three operational stages: Before Thermal Runaway (BR), After Thermal Runaway–Decrease in Microwave Power (ARD), and After Thermal Runaway– Increase in Microwave Power (ARI). Temperature and gas composition were continuously monitored throughout the experiments. The rock-only experiments demonstrated that hydrogen generation can occur intrinsically from hydrocarbon-bearing rocks under microwave heating, even without externally injected methane. However, hydrogen production did not correlate solely with kerogen content, indicating that mineralogy strongly influences hydrogen-generation pathways. Correlation analyses suggested that kerogen decomposition initially generated CH4, CO, and CO2, followed by secondary hydrocarbon reactions associated with H2 and C2 hydrocarbon formation. Methane-assisted experiments substantially increased hydrogen production; however, identical methane injection rates produced significantly different hydrogen yields among the rock samples, confirming that mineralogical composition controls methaneconversion behavior under microwave heating conditions. The addition of Ca(OH)2 significantly altered carbon evolution behavior in a stage-dependent manner. During the BR stage, Ca(OH)2 reduced gas-phase CO2 production, particularly in carbonate-rich rocks, indicating favorable conditions for in-situ carbonation and carbonate deposition prior to extensive thermal decomposition. The suppression of CO2 during BR became more pronounced with increasing carbonate content of the rock system. After thermal runaway, carbonate-containing systems exhibited enhanced hydrogen generation behavior, suggesting that carbonatederived mineral transformations and carbonate-mediated reactions contribute to high-temperature hydrogen-generation pathways. In carbonate-poor rocks, Ca(OH)2 addition enabled simultaneous URTeC 4493775 2 enhancement of hydrogen production and partial suppression of CO2 release during the post-runaway stages. Overall, the results demonstrate that microwave-assisted hydrogen generation is governed by dynamically evolving interactions among kerogen decomposition, methane activation, mineral transformations, carbonate formation/decomposition, and thermal runaway behavior. This work introduces in-situ carbonate engineering with Ca(OH)2 as a strategy for coupling hydrogen generation with partial insitu carbon management under microwave heating conditions.

03 NATURAL GAS↗
DOE OSTI2023

Tensor train continuous time solver for quantum impurity models

The simulation of strongly correlated quantum impurity models is a significant challenge in modern condensed matter physics that has multiple important applications. Thus far, the most successful methods for approaching this challenge involve Monte Carlo techniques that accurately and reliably sample perturbative expansions to any order. However, the cost of obtaining high precision through these methods is high. Recently, tensor train decomposition techniques have been developed as an alternative to Monte Carlo integration. In this study, we apply these techniques to the single-impurity Anderson model at equilibrium by calculating the systematic expansion in power of the hybridization of the impurity with the bath. Furthermore, we demonstrate the performance of the method in a paradigmatic application, examining the first-order phase transition on the infinite-dimensional Bethe lattice, which can be mapped to an impurity model through dynamical mean field theory. Our results indicate that using tensor train decomposition schemes allows the calculation of finite-temperature Green's functions and thermodynamic observables with unprecedented accuracy. The methodology holds promise for future applications to frustrated multiorbital systems, using a combination of partially summed series with other techniques pioneered in diagrammatic and continuous time quantum Monte Carlo.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Influence of soil depth, irrigation, and plant genotype on the soil microbiome, metaphenome, and carbon chemistry

ABSTRACT Climate change is causing an increase in drought in many soil ecosystems and a loss of soil organic carbon. Calcareous soils may partially mitigate these losses via carbon capture and storage. Here, we aimed to determine how irrigation-supplied soil moisture and perennial plants impact biotic and abiotic soil properties that underpin deep soil carbon chemistry in an unfertilized calcareous soil. Soil was sampled up to 1 m in depth from irrigated and planted field treatments and was analyzed using a suite of omics and chemical analyses. The soil microbial community composition was impacted more by irrigation and plant cover treatments than by soil depth. By contrast, metabolomes, lipidomes, and proteomes differed more with soil depth than treatments. Deep soil (>50 cm) had higher soil pH and calcium concentrations and higher levels of organic acids, bicarbonate, and triacylglycerides. By contrast, surface soil (0–5 cm) had higher concentrations of soil organic matter, organic carbon, oxidizable carbon, and total nitrogen. Surface soils also had higher amounts of sugars, sugar alcohols, phosphocholines, and proteins that reflect osmotic and oxidative stress responses. The lipidome was more responsive to perennial tall wheatgrass treatments compared to the metabolome or proteome, with a striking change in diacylglyceride composition. Permanganate oxidizable carbon was more consistently correlated to metabolites and proteins than soil organic and inorganic carbon and soil organic matter. This study reveals specific compounds that reflect differences in organic, inorganic, and oxidizable soil carbon fractions that are impacted by interactions between irrigation-supplied moisture and plant cover in calcareous soil profiles. IMPORTANCE Carbon is cycled through the air, plants, and belowground environment. Understanding soil carbon cycling in deep soil profiles will be important to mitigate climate change. Soil carbon cycling is impacted by water, plants, and soil microorganisms, in addition to soil mineralogy. Measuring biotic and abiotic soil properties provides a perspective of how soil microorganisms interact with the surrounding chemical environment. This study emphasizes the importance of considering biotic interactions with inorganic and oxidizable soil carbon in addition to total organic carbon in carbonate-containing soils for better informing soil carbon management decisions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Spirometric traits show quantile-dependent heritability, which may contribute to their gene-environment interactions with smoking and pollution

Background: “Quantile-dependent expressivity” refers to a genetic effect that is dependent upon whether the phenotype (e.g., spirometric data) is high or low relative to its population distribution. Forced vital capacity (FVC), forced expiratory volume in 1 second (FEV 1 ), and the FEV 1 /FVC ratio are moderately heritable spirometric traits. The aim of the analyses is to test whether their heritability (h 2 ) is constant over all quantiles of their distribution. Methods: Quantile regression was applied to the mean age, sex, height and smoking-adjusted spirometric data over multiple visits in 9,993 offspring-parent pairs and 1,930 sibships from the Framingham Heart Study to obtain robust estimates of offspring-parent (β OP ), offspring-midparent (β OM ), and full-sib regression slopes (β FS ). Nonparametric significance levels were obtained from 1,000 bootstrap samples. β OP s were used as simple indicators of quantile-specific heritability (i.e.,h 2 = 2β OP /(1+r spouse ), where r spouse was the correlation between spouses). Results: β OP ± standard error (SE) decreased by 0.0009 ± 0.0003 (P = 0.003) with every one-percent increment in the population distribution of FEV 1 /FVC, i.e., β OP ± SE were: 0.182 ± 0.031, 0.152 ± 0.015; 0.136 ± 0.011; 0.121 ± 0.013; and 0.099 ± 0.013 at the 10th, 25th, 50th, 75th, and 90th percentiles of the FEV 1 /FVC distribution, respectively. These correspond to h 2 ± SEs of 0.350 ± 0.060 at the 10th, 0.292 ± 0.029 at the 25th, 0.262 ± 0.020 at the 50th, 0.234 ± 0.025 at the 75th, and 0.191 ± 0.025 at the 90th percentiles of the FEV 1 /FVC ratio. Maximum mid-expiratory flow (MMEF) h 2 ± SEs increased 0.0025 ± 0.0007 (P = 0.0004) with every one-percent increment in its distribution, i.e.: 0.467 ± 0.046, 0.467 ± 0.033, 0.554 ± 0.038, 0.615 ± 0.042, and 0.675 ± 0.060 at the 10th, 25th, 50th, 75th, and 90th percentiles of its distribution. This was due to forced expiratory flow at 75% of FVC (FEF75%), whose quantile-specific h 2 increased an average of 0.0042 ± 0.0008 for every one-percent increment in its distribution. It is speculated that previously reported gene-environment interactions may be partially attributable to quantile-specific h 2 , i.e., greater heritability in individuals with lower FEV 1 /FVC due to smoking or airborne particles exposure vs. nonsmoking, unexposed individuals. Conclusion: Heritabilities of FEV 1 /FVC, MMEF, and FEF75% from quantile-regression of offspring-parent and sibling spirometric data suggest their quantile-dependent expressivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evolution of Charge-Lattice Dynamics across the Kuramoto Synchronization Phase Diagram of Quantum Tunneling Polarons in Cuprate Superconductors

Because of its sensitivity to the instantaneous structure factor, S(Q,t = 0), Extended X-ray Absorption Fine Structure (EXAFS) is a powerful tool for probing the dynamic structure of condensed matter systems in which the charge and lattice dynamics are coupled. When applied to hole-doped cuprate superconductors, EXAFS has revealed the presence of internal quantum tunneling polarons (IQTPs). An IQTP arises in EXAFS as a two-site distribution for certain Cu–O pairs, which is also duplicated in inelastic scattering but not observed in standard diffraction measurements. The Cu–Sr pair distribution has been found to be highly anharmonic and strongly correlated to both the IQTPs and to superconductivity, as, for example, in YSr 2 Cu 2.75 Mo 0.25 O 7.54 (T c =84 K). In order to describe such nontrivial, anharmonic charge-lattice dynamics, we have proposed a model Hamiltonian for a prototype six-atom cluster, in which two Cu-apical-O IQTPs are charge-transfer bridged through Cu atoms by an O atom in the CuO 2 plane and are anharmonically coupled via a Sr atom. By applying an exact diagonalization procedure to this cluster, we have verified that our model indeed produces an intricate interplay between charge and lattice dynamics. Then, by using the Kuramoto model for the synchronization of coupled quantum oscillators, we have found a first-order phase transition for the IQTPs into a synchronized, phase-locked phase. Most importantly, we have shown that this transition results specifically from the anharmonicity. Finally, we have provided a phase diagram showing the onset of the phase-locking of IQTPs as a function of the charge-lattice and anharmonic couplings in our model. We have found that the charge, initially confined to the apical oxygens, is partially pumped into the CuO 2 plane in the synchronized phase, which suggests a possible connection between the synchronized dynamic structure and high-temperature superconductivity (HTSC) in doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

X-ray Photoelectron Spectroscopy Analysis of Nafion-Containing Samples: Pitfalls, Protocols, and Perceptions of Physicochemical Properties

X-ray photoelectron spectroscopy (XPS) is one of the most common techniques used to analyze the surface composition of catalysts and support materials used in polymer electrolyte membrane (PEM) fuel cells and electrolyzers, providing important insights for further improvement of their properties. Characterization of catalyst layers (CLs) is more challenging, which can be at least partially attributed to the instability of ionomer materials such as Nafion during measurements. This work explores the stability of Nafion during XPS measurements, illuminating and addressing Nafion degradation concerns. The extent of Nafion damage as a function of XPS instrumentation, measurement conditions, and sample properties was evaluated across multiple instruments. Results revealed that significant Nafion damage to the ion-conducting sulfonic acid species (>50% loss in sulfur signal) may occur in a relatively short time frame (tens of minutes) depending on the exact nature of the sample and XPS instrument. This motivated the development and validation of a multipoint XPS data acquisition protocol that minimizes Nafion damage, resulting in reliable data acquisition by avoiding significant artifacts from Nafion instability. The developed protocol was then used to analyze both thin film ionomer samples and Pt/C-based CLs. Comparison of PEM fuel cell CLs to Nafion thin films revealed several changes in Nafion spectral features attributed to charge transfer due to interaction with conductive catalyst and support species. This study provides a method to reliably characterize ionomer-containing samples, facilitating fundamental studies of the catalyst-ionomer interface and more applied investigations of structure-processing-performance correlations in PEM fuel cell and electrolyzer CLs.

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

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The LSST AGN Data Challenge: Selection Methods

Abstract Development of the Rubin Observatory Legacy Survey of Space and Time (LSST) includes a series of Data Challenges (DCs) arranged by various LSST Scientific Collaborations that are taking place during the project's preoperational phase. The AGN Science Collaboration Data Challenge (AGNSC-DC) is a partial prototype of the expected LSST data on active galactic nuclei (AGNs), aimed at validating machine learning approaches for AGN selection and characterization in large surveys like LSST. The AGNSC-DC took place in 2021, focusing on accuracy, robustness, and scalability. The training and the blinded data sets were constructed to mimic the future LSST release catalogs using the data from the Sloan Digital Sky Survey Stripe 82 region and the XMM-Newton Large Scale Structure Survey region. Data features were divided into astrometry, photometry, color, morphology, redshift, and class label with the addition of variability features and images. We present the results of four submitted solutions to DCs using both classical and machine learning methods. We systematically test the performance of supervised models (support vector machine, random forest, extreme gradient boosting, artificial neural network, convolutional neural network) and unsupervised ones (deep embedding clustering) when applied to the problem of classifying/clustering sources as stars, galaxies, or AGNs. We obtained classification accuracy of 97.5% for supervised models and clustering accuracy of 96.0% for unsupervised ones and 95.0% with a classic approach for a blinded data set. We find that variability features significantly improve the accuracy of the trained models, and correlation analysis among different bands enables a fast and inexpensive first-order selection of quasar candidates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Harvest Intensity Effects on Carbon Stocks and Biodiversity Are Dependent on Regional Climate in Douglas-Fir Forests of British Columbia

Temperate forests provide crucial ecosystems services as living sinks for atmospheric carbon (C) and repositories of biodiversity. Applying harvesting at intensities that minimize losses offers one means for mitigating global change. However, little is known of overstory retention levels that best conserve ecosystem services in different regional climates, and likewise as climate changes. To quantify the effect of harvest intensity on C stocks and biodiversity, we compared five harvesting intensities (clearcutting, seedtree retention, 30% patch retention, 60% patch retention, and uncut controls) across a climatic aridity gradient that ranged from humid to semi-arid in the Douglas-fir (Pseudotsuga menziesii) forests of British Columbia. We found that increased harvesting intensity reduced total ecosystem, above ground, and live tree C stocks 1 year post-harvest, and the magnitude of these losses were negatively correlated with climatic aridity. In humid forests, total ecosystem C ranged from 50% loss following clearcut harvest, to 30% loss following large patch retention harvest. In arid forests this range was 60 to 8% loss, respectively. Where lower retention harvests are sought, the small patch retention treatment protected both C stocks and biodiversity in the arid forests, whereas the seedtree method performed as well or better in the humid forests. Below ground C stocks declined by an average of 29% after harvesting, with almost all of the loss from the forest floor and none from the mineral soil. Of the secondary pools, standing and coarse deadwood declined in all harvesting treatments regardless of cutting intensity or aridity, while C stocks in fine fuels and stumps increased. The understory plant C pool declined across all harvesting intensities in the humid forests, but increased in arid forests. Shannon’s diversity and richness of tree and bryoid species declined with harvesting intensity, where tree species losses were greatest in the humid forests and bryoid losses greatest in arid forests. Shrub and herb species were unaffected. This study showed that the highest retention level was best at reducing losses in C stocks and biodiversity, and clearcutting the poorest, and while partial retention of canopy trees can reduce losses in these ecosystem services, outcomes will vary with climatic aridity.

54 ENVIRONMENTAL SCIENCES↗