Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility characterization”

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 523 records · Page 29

Characterization of cathode keeper wear by surface layer activation

In this study, the erosion rates of the discharge cathode keeper in a 30 cm NSTAR configuration ion thruster were measured using a technique known as Surface Layer Activation (SLA). This diagnostic technique involves producing a radioactive tracer in a given surface by bombardment with high energy ions. The decrease in activity of the tracer material may be monitored as the surface is subjected to wear processes and correlated to a depth calibration curve, yielding the eroded depth. Analysis of the activities was achieved through a gamma spectroscopy system. The primary objectives of this investigation were to reproduce erosion data observed in previous wear studies in order to validate the technique, and to determine the effect of different engine operating parameters on erosion rate. The erosion profile at the TH 15 (23 kw) setting observed during the 8200 hour Life Demonstration Test (LDT) was reproduced. The maximum keeper erosion rate at this setting was determined to be 0.085 pm/hr. Testing at the TH 8 (1.4 kw) setting demonstrated lower erosion rates than TH 15, along with a different wear profile. Varying the keeper voltage was shown to have a significant effect on the erosion, with a positive bias with respect to cathode potential decreasing the erosion rate significantly. Accurate measurements were achieved after operating times of only 40 to 70 hours, a significant improvement over other erosion diagnostic methods.

surface layer activation↗

An Optothermal Raman Model Accounting for Convection, Radiation, and Substrate Effects for Evaluation of the Thermal Conductivity of Two‐Dimensional Coatings

Thermal conductivity of two‐dimensional (2D) materials is important to characterize as this will govern thermal transport physics in the various thermal and energy applications benefitting from the unique properties of these materials. Here, analytical and numerical techniques are presented to determine the thermal conductivity and interfacial conductance of a 2D material coating using confocal micro‐Raman spectroscopy as a noncontact diagnostic. While several methods have been proposed for calculating the values of thermal conductivity, interfacial conductance, and spatial temperature profile, these models often do not consider effects of radiation, convection, and substrate resistance on the temperature and heat flux profiles. We present a model to calculate the thermal conductivity and interfacial conductance which accounts for convection, radiation, and substrate effects to characterize a variety of 2D material coatings, which is demonstrated using large area graphene transferred onto copper and nickel substrates. Convection and radiation effects are found to have a negligible effect on the temperature profile of supported coatings, while the substrate effects have a considerable effect and therefore impact the quantification of thermal conductivity and interfacial conductance in addition to limiting the spatial resolution of the technique. These findings will allow for more accurate and reproducible extraction of the thermal conductivity of graphene and other 2D material coatings when employing noncontact optothermal Raman spectroscopy methods.

2D coatings↗

Evaluating E3SM Global Storm‐Resolving Model Simulations of Deep Convection: Insights From DP‐SCREAM During TRACER

Global Storm-Resolving Models (GSRMs) are becoming increasingly vital for advancing climate modeling and improving the prediction of extreme weather events. Houston, a coastal region frequently affected by deep convective storms, offers an ideal setting to evaluate the ability of GSRMs to simulate deep convection. This study assesses the performance of the Doubly Periodic Simple Cloud-Resolving E3SM (Energy Exascale Earth System Model) Atmosphere Model (DP-SCREAM) using observations from the TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign. DP-SCREAM effectively reproduces the diurnal cycles of clouds and precipitation, demonstrating much greater skill than the E3SM single column model. The DP-SCREAM is demonstrated to be applicable to coastal regions, partially due to the forcing data sets already capturing the influence of breezes. DP-SCREAM also replicates biases persistent in the global version of SCREAM: the underrepresentation of boundary layer shallow clouds, a lack of mid-level congestus clouds, and the popcorn convection, characterized by small and disorganized convective cells generating the strongest precipitation. To investigate these issues, two sensitivity experiments were conducted: increasing the mixing length and scaling up the buoyancy flux within the Simplified Higher Order Closure scheme. Increasing the mixing length improved mid-level congestus representation and reduced unrealistic early morning fog occurrence. Enhancing buoyancy flux only marginally improved the bias of underproduced big convective cells. In conclusion, an additional resolution sensitivity test at 0.5 km grid spacing demonstrated that a refined horizontal resolution alone is insufficient to resolve these biases.

54 ENVIRONMENTAL SCIENCES↗

Let it rip: The mechanics of self-bisection in asexual planarians determines their population reproductive strategies

Asexual freshwater planarians reproduce by transverse bisection (binary fission) into two pieces. This process produces a head and a tail, which fully regenerate within 1-2 weeks. How planarians split into two offspring - using only their musculature and substrate traction – is a challenging biomechanics problem. We found that three different species, Dugesia japonica, Girardia tigrina and Schmidtea mediterranea, have evolved three different mechanical solutions to self-bisect. Using time lapse imaging of the fission process, we quantitatively characterize the main steps of division in the three species and extract the distinct and shared key features. Across the three species, planarians actively alter their body shape, regulate substrate traction, and use their muscles to generate tensile stresses large enough to overcome the ultimate tensile strength of the tissue. Moreover, we show that how each planarian species divides dictates how resources are split among its offspring. This ultimately determines offspring survival and reproductive success. Furthermore, heterospecific differences in the mechanics of self-bisection of individual worms explain the observed differences in the population reproductive strategies of different planarian species.

59 BASIC BIOLOGICAL SCIENCES↗

Multicode benchmark on simulated Ti K-edge x-ray absorption spectra of Ti-O compounds

X-ray absorption spectroscopy (XAS) is an element-specific materials characterization technique that is sensitive to structural and electronic properties. First-principles simulated XAS has been widely used as a powerful tool to interpret experimental spectra and draw physical insights. Recently, there has also been growing interest in building computational XAS databases to enable data analytics and machine learning applications. However, there are non-trivial differences among commonly used XAS simulation codes, both in underlying theoretical formalism and in technical implementation. Reliable and reproducible computational XAS databases require systematic benchmark studies. In this work, we benchmarked Ti K-edge XAS simulations of ten representative Ti-O binary compounds, which we refer to as the Ti-O-10 dataset, using three state-of-the-art codes: XSPECTRA, OCEAN and exciting. We systematically studied the convergence behavior with respect to the input parameters and developed a workflow to automate and standardize the calculations to ensure converged spectra. Our benchmark comparison considers a 35 eV spectral range starting from the K-edge onset, representative of widely used near-edge spectra. Quantitative comparison over this range is based on Spearman’s rank correlation score (r sp ). Our results show: (1) the two Bethe-Salpeter equation (BSE) codes (OCEAN and exciting) have excellent agreement with an average $r$ $sp$ of 0.998; (2) good agreement is obtained between the core-hole potential code (XSPECTRA) and BSE codes (OCEAN and exciting) with an average $r$ $sp$ of 0.990, and this smaller $r$ $sp$ reflects the noticeable differences in the main edge spectral shape that can be primarily attributed to the difference in the strength of the screened core-hole potential; (3) simulations from both methods overall reproduce well the main experimental spectral features of rutile and anatase, and the different treatments of the screened core-hole potential have visible impact on pre-edge intensities and the peak ratio of the main edge; (4) there exist moderate differences in the relative edge alignment of the three codes with a standard deviation of about 0.2 eV, which arise from multiple contributions including the frozen core approximation, final state effects, and different approximations used for the self-energy correction. In conclusion, our benchmark study provides important standards for first-principles XAS simulations with broad impact in data-driven XAS analysis.

36 MATERIALS SCIENCE↗

Development of high Tc (greater than 100 K) Bi, Tl and Y-based materials as superconducting circuit elements

Experimental work on this project over the last four years has resulted in establishing processing and characterization techniques for producing both the Bi-based and Tl-based superconductors in their high temperature (2223) forms. In the bulk, dry pressed form, maximum critical temperatures (Tc) of 108.2 K and 117.8 K, respectively, were measured. Results have further shown that the Bi and Tl-based superconducting materials in bulk form are noticeably different from the Y-based 123 material in that superconductivity is considerably harder to achieve, maintain, and reproduce. This is due primarily to the difficulty in obtaining the higher Tc phase in pure form since it commonly co-exists with other undesirable, lower Tc phases. In particular, it has been found that long processing times for calcining and firing (20 - 200 hrs.) and close control of temperatures which are very near the melting point are required in order to obtain higher proportions of the desirable, high Tc (2223) phase. Thus far, the BSCCO bulk materials has been prepared in uniaxially pressed, hot pressed, and tapecast form. The uniaxially pressed material has been synthesized by the mixed oxide, coprecipitation, and melt quenching processes. The tapecast and hot pressed materials have been prepared via the mixed oxide process. In addition, thick films of BSCCO (2223 phase) have been prepared by screen printing on to yttria and magnesia stabilized zirconia with only moderate success; i.e., superconductivity was achieved in these thick films, but the highest Tc obtained in these films was 89.0 K. The Tc's of the bulk hot pressed, tapecast, and screen printed thick film materials were found to be 108.2, 102.4, and 89.0 K, respectively.

Haertling, Gene↗

Evaluation of a Reference-Free Collision Cross Section Calibration Strategy for Proteomics Using SLIM-Based High-Resolution Ion Mobility Spectrometry–Mass Spectrometry

Ion mobility spectrometry (IMS) is a gas-phase analytical technique that separates ions with different sizes and shapes and is compatible with mass spectrometry (MS) to provide an additional separation dimension. The rapid nature of the IMS separation combined with the high sensitivity of MS-based detection and the ability to derive structural information on analytes in the form of the property collision cross section (CCS) makes IMS particularly well-suited for characterizing complex samples in -omics applications. In such applications, the quality of CCS from IMS measurements is critical to confident annotation of the detected components in the complex -omics samples. However, most IMS instrumentation in mainstream use requires calibration to calculate CCS from measured arrival times, with the most notable exception being drift tube IMS measurements using multifield methods. The strategy for calibrating CCS values, particularly selection of appropriate calibrants, has important implications for CCS accuracy, reproducibility, and transferability between laboratories. The conventional approach to CCS calibration involves explicitly defining calibrants ahead of data acquisition and crucially relies upon availability of reference CCS values. In this work, we present a novel reference-free approach to CCS calibration which leverages trends among putatively identified features and computational CCS prediction to conduct calibrations post-data acquisition and without relying on explicitly defined calibrants. We demonstrated the utility of this reference-free CCS calibration strategy for proteomics application using high-resolution structures for lossless ion manipulations (SLIM)-based IMS-MS. In conclusion, we first validated the accuracy of CCS values using a set of synthetic peptides and then demonstrated using a complex peptide sample from cell lysate.

59 BASIC BIOLOGICAL SCIENCES↗

Systematics in asteroseismic modelling: application of a correlated noise model for oscillation frequencies

ABSTRACT The detailed modelling of stellar oscillations is a powerful approach to characterizing stars. However, poor treatment of systematics in theoretical models leads to misinterpretations of stars. Here, we propose a more principled statistical treatment for the systematics to be applied to fitting individual mode frequencies with a typical stellar model grid. We introduce a correlated noise model based on a Gaussian process (GP) kernel to describe the systematics given that mode frequency systematics are expected to be highly correlated. We show that tuning the GP kernel can reproduce general features of frequency variations for changing model input physics and fundamental parameters. Fits with the correlated noise model better recover stellar parameters than traditional methods that either ignore the systematics or treat them as uncorrelated noise.

Li, Tanda (ORCID:0000000163962563)↗

Failure Analysis of a New Polyamide-Based Fluoropolymer-Free Backsheet After Combined-Accelerated Stress Testing

The viability of novel coextruded, fluoropolymer-free backsheets for photovoltaic (PV) modules has been questioned as a result of a large number of early-life backsheet failures in PV installations containing one of the earliest co-extruded polyamide (PA)-based backsheet to reach the market, “AAA.” New PV reliability testing protocols have been recently developed and applied to backsheets to reproduce failures observed in the field and evaluate the durability of novel backsheet materials and designs prior to commercialization. A new co-extruded PA-based backsheet was tested using combined-accelerated stress testing (C-AST) and demonstrated a greater lifetime than AAA, and some other fluoropolymer-based backsheets such as polyvinylidene fluoride. The improved PA-based backsheet also eventually failed by through-thickness cracking. Using surface and bulk material characterization techniques, we performed a comprehensive study of material properties before and after the stress testing. Aging of the backsheet resulted in an increase of surface roughness by erosion of the outer PA layer. However the failure is more likely related to an increase in crystallinity of the polyolefin core layer reducing the backsheet tearing energy. The analysis can ultimately inform on the specific weaknesses of the materials so that the manufacturer can improve the backsheet design to extend its lifetime.

14 SOLAR ENERGY↗

Impacts of long-range transport of aerosols on marine-boundary-layer clouds in the eastern North Atlantic

Vertical profiles of aerosols are inadequately observed and poorly represented in climate models, contributing to the current large uncertainty associated with aerosol–cloud interactions. The US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) aircraft field campaign near the Azores islands provided ample observations of vertical distributions of aerosol and cloud properties. Here we utilize the in situ aircraft measurements from the ACE-ENA and ground-based remote-sensing data along with an aerosol-aware Weather Research and Forecast (WRF) model to characterize the aerosols due to long-range transport over a remote region and to assess their possible influence on marine-boundary-layer (MBL) clouds. The vertical profiles of aerosol and cloud properties measured via aircraft during the ACE-ENA campaign provide detailed information revealing the physical contact between transported aerosols and MBL clouds. The European Centre for Medium-Range Weather Forecasts Copernicus Atmosphere Monitoring Service (ECMWF-CAMS) aerosol reanalysis data can reproduce the key features of aerosol vertical profiles in the remote region. The cloud-resolving WRF sensitivity experiments with distinctive aerosol profiles suggest that the transported aerosols and MBL cloud interactions (ACIs) require not only aerosol plumes to get close to the marine-boundary-layer top but also large cloud top height variations. Based on those criteria, the observations show that the occurrence of ACIs involving the transport of aerosol over the eastern North Atlantic (ENA) is about 62 % in summer. For the case with noticeable long-range-transport aerosol effects on MBL clouds, the susceptibilities of droplet effective radius and liquid water content are -0.11 and +0.14, respectively. When varying by a similar magnitude, aerosols originating from the boundary layer exert larger microphysical influence on MBL clouds than those entrained from the free troposphere.

54 ENVIRONMENTAL SCIENCES↗

Airborne and Shipborne Polarimetric Measurements over Open Ocean and Coastal Waters: Intercomparisons and Implications for Spaceborne Observations

Comprehensive polarimetric closure is demonstrated using observations from two in-situ polarimeters and Vector Radiative Transfer (VRT) modeling. During the Ship-Aircraft Bio-Optical Research (SABOR) campaign, the novel CCNY HyperSAS-POL polarimeter was mounted on the bow of the R/V Endeavor and acquired hyperspectral measurements from just above the surface of the ocean, while the NASA GISS Research Scanning Polarimeter was deployed onboard the NASA LaRC's King Air UC-12B aircraft. State-of-the-art, ancillary measurements were used to characterize the atmospheric and marine contributions in the VRT model, including those of the High Spectral Resolution Lidar (HSRL), the AErosol RObotic NETwork for Ocean Color (AERONET-OC), a profiling WETLabs ac-9 spectrometer and the Multi-spectral Volume Scattering Meter (MVSM). An open-ocean and a coastal scene are analyzed, both affected by complex aerosol conditions. In each of the two cases, it is found that the model is able to accurately reproduce the Stokes components measured simultaneously by each polarimeter at different geometries and viewing altitudes. These results are mostly encouraging, considering the different deployment strategies of RSP and HyperSAS-POL, which imply very different sensitivities to the atmospheric and ocean contributions, and open new opportunities in above-water polarimetric measurements. Furthermore, the signal originating from each scene was propagated to the top of the atmosphere to explore the sensitivity of polarimetric spaceborne observations to changes in the water type. As expected, adding polarization as a measurement capability benefits the detection of such changes, reinforcing the merits of the full-Stokes treatment in modeling the impact of atmospheric and oceanic constituents on remote sensing observations.

Matteo Ottaviani↗

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗

A forward model-based validation of cardiovascular system identification

We present a theoretical evaluation of a cardiovascular system identification method that we previously developed for the analysis of beat-to-beat fluctuations in noninvasively measured heart rate, arterial blood pressure, and instantaneous lung volume. The method provides a dynamical characterization of the important autonomic and mechanical mechanisms responsible for coupling the fluctuations (inverse modeling). To carry out the evaluation, we developed a computational model of the cardiovascular system capable of generating realistic beat-to-beat variability (forward modeling). We applied the method to data generated from the forward model and compared the resulting estimated dynamics with the actual dynamics of the forward model, which were either precisely known or easily determined. We found that the estimated dynamics corresponded to the actual dynamics and that this correspondence was robust to forward model uncertainty. We also demonstrated the sensitivity of the method in detecting small changes in parameters characterizing autonomic function in the forward model. These results provide confidence in the performance of the cardiovascular system identification method when applied to experimental data.

NASA Discipline Cardiopulmonary↗

Influence of Terminal Carboxyl Groups on the Structure and Reactivity of Functionalized m-Carboranethiolate Self-Assembled Monolayers

The structure and function of self-assembled monolayers (SAMs) at the nanoscale are determined by the steric and electronic effects of their building blocks. Carboranethiol molecules form pristine monolayers that provide tunable two-dimensional systems to probe lateral and interfacial interactions. Additional ω-functionality, such as carboxyl groups, can be introduced to change the properties of the exposed surfaces. Here, two geometrically similar isomeric m-carborane analogues of m-mercaptobenzoic acid, 1-COOH-7-SH-1,7-C 2 B 10 H 10 and racem-1-COOH-9-SH-1,7-C 2 B 10 H 10 , are characterized and their SAMs on Au{111} are examined. The latter isomer belongs to the rare group of chiral cage molecules and becomes, to our knowledge, the first example assembled on Au{111}. Although different in symmetry, molecules of both isomers assemble into similar hexagonal surface patterns. The nearest-neighbor spacing of 8.4 ± 0.4 Å is larger than that of non-carboxylated isomers, consistent with the increased steric demands of the carboxyl groups. Computational modeling reproduced this spacing and suggests a tilt relative to the surface normal. However, tilt domains are not observed experimentally, suggesting the presence of strong lateral interactions. Analyses of the influence of the functional groups through the pseudo-aromatic m-carborane skeleton showed that the thiol group attached to either carbon or boron atoms increases the carboxyl group acidity in solution. In contrast, the acidity of the exposed carboxyl group in the SAMs decreases upon surface attachment; computational analyses suggest that the driving force of this shift is the dielectric of the environment in the monolayer as a result of confined intermolecular interactions, proximity to the Au surface, and partial desolvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗

Sub-pilot-scale Production of High-Value Products from U.S. Coals

Investigators from the University of Utah, University of Wyoming and Marshall University pursued a program to study the conversion of raw coal to high-value products of carbon fiber and silicon carbide. Team members also developed an initial framework for a data portal that can incorporate laboratory data on coal processing and product quality, and also work with tools for machine learning for data analysis, data visualization and economic assessment. Experimental R&D efforts focused on the conversion of raw coal to coal tar and other byproducts, and the resulting tar intermediates were upgraded to form anisotropic and isotropic pitch materials. These pitch materials were produced from coal using both thermal (pyrolysis) and chemical (mild solvolysis liquefaction) decomposition of raw coal. Four different coals were studied: Utah bituminous coal (Sufco), Wyoming PRB coal (Black Thunder), Illinois bituminous coal (Illinois #6), and West Virginia bituminous coal (Flying Eagle). Both metallurgical-grade coking coals and lower-grade steam coals were investigated, and controlled secondary gas-phase reactions were used during a two-stage pyrolysis process to induce cracking and condensation reactions among the pyrolytic tar species. This approach successfully improved the performance of the lower grade coals for yielding pitch materials, with properties more consistent with a commercial-grade pitch that had previously demonstrated success for quality carbon fiber production. The use of waste plastic materials was also studied, to help improve physical and chemical characteristics of the intermediate tars and final pitch product; in particular, for lowering the pitch softening point to an acceptable level for melt spinning carbon fiber. Mild solvolysis liquefaction was also used as a method for producing pitch for carbon fiber production. As expected, significantly higher pitch yields were obtained using this approach, and waste plastic materials were also successfully used to reduce pitch softening point to an acceptable level. The plastic materials were also utilized to create a solvent for the mild solvolysis process, and this plastic-derived solvent was shown to provide results consistent with more expensive commercial chemical solvents, and could thus avoid the need for costly recovery and recycle of a liquefaction solvent. Additional experimental R&D focused on the production of silicon carbide (β-SiC) from the residual char byproduct from pitch production, and also on the production of carbon fiber from the anisotropic pitch. SiC was successfully synthesized using a mixture of residual char and sandstone at a ratio of 1:1. Reaction temperature and residence time were optimized and yielded a product purity of 81%. For carbon fiber production, the most successful pitch samples were obtained from the mild solvolysis liquefaction approach, combined with the use of a plastic (HDPE)-derived solvent. Fiber properties improved over time as laboratory fiber production methodologies improved, and final yields of carbon fiber were obtained with a diameter of 12.14 ± 1.10 um, Modulus of 173.73 ± 15.25 GPa, and Tensile Strength of 1.04 ± 0.10 GPa. A proof-of-concept Modern Community Research Data Portal (MCRDP) was developed and deployed for coal and coal-derived pitch characterization, with the full support of (i) remote web-based access, (ii) distributed analysis, (iii) interactive visualization and exploration, (iv) shared and long-term data access, (v) advanced query capabilities and (vi) real-time collaboration. The Coal to Products Data Portal “coaltoproducts.org” provides researchers with space to store and share data within a project, tools for analyzing and understanding data for scientific investigation, and the ability to publish data to the broader community for reproducibility. The portal leverages the Material Commons 2.0 (MC) platform developed by the Center for PRedictive Integrated Structural Materials Science (PRISMS) of the University of Michigan, to achieve long-term longevity of data collections and, more importantly, collaborative science. A number of data visualization tools were also assessed and implemented for interrogating the experimental and modeling data. The machine learning portion of this project analyzed datasets from two different coal conversion processes performed on a diverse set of coal samples from both the coal pyrolysis experiments and the solvent liquefaction experiments. The work was initiated by exploring standard regression models on the pyrolysis data, aiming to understand the impact of sample characteristics and processing conditions on key product metrics. Over the course of the project, the focus expanded to include a variety of machine learning tools, delving into both supervised and unsupervised learning methods. Models tested on the pyrolysis data included linear, ridge, lasso, elastic-net, Gaussian process, random forest regression, and AutoSklearn, and the approach was continually refined to enhance predictive accuracy and model interpretability. Similar techniques were applied to the liquefaction data with an additional focus on feature engineering. Along with mesophase content, additional outputs of interest were the pitch yield, softening point, and QI content. Insights derived from these analyses are crucial in determining the factors influencing the quality and yield of coal-derived products. As the work progressed, the research evolved from foundational model comparisons to analyses of random forests, decision paths, and feature importance scores. A thorough market analysis was performed to examine the prospects of coal-based carbon fibers. The best opportunities for coal come from its lower and more stable price relative to petroleum, particularly for subbituminous coals, which is the primary advantage that a coal refinery may have over a petroleum refinery. Before a commercial CTP production facility can be modeled, however, several things need to be understood regarding the nature of the would-be coal refinery. These include the technology to be deployed, the size of facility, the volume(s) of co-product(s), and the waste and emissions profile of the plant. The volume of co-products and waste may be substantial and will require separate market analysis to ensure viability. In the near-term, the importance of coal tar pitch, in the form of carbon pitch, to the aluminum and steel industries is likely to overshadow the alternative use of this material as an input for carbon fiber. The importance of steel and aluminum in building materials, and the need for carbon materials in their manufacturing, will ensure that demand for these products remains for the long run. In addition, carbon fiber may also be the best substitute for steel and aluminum well into the future. While society will eventually be able to shift production of much of its electricity needs to renewables, it will not be able to shift away from fossil fuels for production of high-strength construction and vehicular materials. Demand for carbon fiber is expected to increase quickly, but the volume of carbon fiber and the amount of coal that would be needed to produce even a sizeable share of this market may still be relatively small compared to current coal production. Thus, other coal-based products like graphene, graphite, carbon foams, resins, and carbon-based building products will play important roles in sustaining coal production as coal-fired power generation continues to decline.

01 COAL, LIGNITE, AND PEAT↗

Automated eddy current analysis of materials

This research effort focused on the use of eddy current techniques for characterizing flaws in graphite-based filament-wound cylindrical structures. A major emphasis was on incorporating artificial intelligence techniques into the signal analysis portion of the inspection process. Developing an eddy current scanning system using a commercial robot for inspecting graphite structures (and others) has been a goal in the overall concept and is essential for the final implementation for expert system interpretation. Manual scans, as performed in the preliminary work here, do not provide sufficiently reproducible eddy current signatures to be easily built into a real time expert system. The expert systems approach to eddy current signal analysis requires that a suitable knowledge base exist in which correct decisions as to the nature of the flaw can be performed. In eddy current or any other expert systems used to analyze signals in real time in a production environment, it is important to simplify computational procedures as much as possible. For that reason, we have chosen to use the measured resistance and reactance values for the preliminary aspects of this work. A simple computation, such as phase angle of the signal, is certainly within the real time processing capability of the computer system. In the work described here, there is a balance between physical measurements and finite element calculations of those measurements. The goal is to evolve into the most cost effective procedures for maintaining the correctness of the knowledge base.

Workman, Gary L.↗

Bayesian Exploration and Surrogate Emulation of Nonlinear Beam-Response Geometry in the LBNF Beamline

Next-generation long-baseline neutrino experiments aim to achieve multi-MW proton beam power while reducing accelerator-induced systematic uncertainties. At Fermilab, the LBNF beamline is designed for 1.2 MW operation with PIP-II and is upgradeable to 2.4 MW. DUNE will probe the three-flavor neutrino paradigm and search for CP violation, requiring precise neutrino-flux normalization and improved control of accelerator-related uncertainties. Within the LBNF beamline, the System for On-Axis Neutrino Detection (SAND) will constrain flux uncertainties using precision near-detector measurements, while the Muon Monitor System (MuMS) will provide beamline diagnostics sensitive to the proton beam, target, and horn configuration. However, the pion phase space relevant for DUNE depends simultaneously on many correlated parameters, including beam centroid, beam width, horn current and alignment, target position, optics shifts, and radiation-induced changes. Consequently, MuMS observables exhibit nonlinear and coupled responses that are difficult to characterize using traditional one-parameter scans. To address this challenge, we are developing a Bayesian Exploration framework coupled to physics-informed surrogate emulators trained on Geant4 beamline simulations. Gaussian-process emulators provide both fast predictions and uncertainty estimates, enabling adaptive selection of new simulation points in beam-parameter space. As an initial demonstration, we construct surrogate emulators for MuMS response observables using a verified simulation campaign spanning proton-beam steering conditions. The emulators reproduce the simulated dependence of MuMS centroid and gradient observables while providing predictive uncertainties, and serve as the foundation for future multidimensional exploration including beam width, horn current, and additional beamline parameters. This work establishes a framework for uncertainty-aware beam monitoring, adaptive simulation campaigns, and rapid beam-response inference for future DUNE operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗