SEARCH · Engineering Papers

Results for “traces”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Rubidium and potassium isotopic variations in chondrites and Mars: Accretion signatures and planetary overprints

As moderately volatile elements, isotopes of Rb and K can trace volatilization processes in planetary bodies. Rubidium isotopic data are however very scarce, especially for non-carbonaceous meteorites. Here, in this study, we report combined Rb and K isotopic data (δ 87/85 Rb and δ 41/39 Κ) for 7 ordinary, 6 enstatite, and 4 Martian meteorite falls to understand the causes for the variations in volatile abundances and isotopic compositions. Bulk Rb and K isotopic compositions of planetary bodies are estimated to be (Table 1): Mars +0.10 ± 0.03 ‰ for Rb and -0.26 ± 0.05 ‰ for K, bulk OCs $-0.12^{+0.15}_{-0.24}$ ‰ for Rb and $-0.72^{+0.28}_{-0.41}$ ‰ for K, bulk ECs $-0.02^{+0.29}_{-0.26}$ ‰ for Rb and $-0.33^{+0.37}_{-0.23}$ ‰ for K. The bulk K isotopic compositions of subgroup OCs are estimated to be $-0.72^{+0.26}_{-0.55}$ ‰ for H chondrites, $-0.71^{+0.23}_{-0.39}$ ‰ for L chondrites, and $-0.77^{+0.63}_{-0.30}$ ‰ for LL chondrites. A broad correlation between the Rb and K isotopic compositions of planetary bodies is observed. The correlation follows a slope that is consistent with kinetic evaporation and condensation processes, suggesting volatility-controlled mass-dependent isotope fractionation (as opposed to nucleosynthetic anomalies). Individual ordinary and enstatite chondrites show large Rb and K isotopic variations (-1.02 to +0.29 ‰ for Rb and -0.91 to -0.15 ‰ for K). Samples of lower metamorphic grades display correlated elemental and isotopic fractionations between Rb and K, while samples of higher metamorphic grades show great scatter, suggesting that chondrite parent-body processes have decoupled the two elements and their isotopes at the sample scale. Several processes could have contributed to the observed isotopic variations of Rb and K, including (i) chondrule “nugget effect”, (ii) volatilization during parent-body thermal metamorphism (heat-induced vaporization and gas transport within parent bodies), (iii) thermal diffusion during parent-body metamorphism, and (iv) impact/shock heating. Quantitative modeling of the first two processes suggests that neither of them could produce isotopic variations large enough to explain the observed isotopic variations. Volatilization during parent-body thermal metamorphism [the scenario (ii)], which has been commonly invoked to explain the isotopic variations of volatile elements, is gas transport-limited and its effect on isotopic fractionations of moderately volatile elements should be negligible. Modeling of diffusion processes suggests that (iii) could produce K isotopic variation comparable to the observed variation. The large isotopic variations in non-carbonaceous meteorites are thus most likely due to diffusive redistribution of K and Rb during metamorphism and/or shock-induced heating and vaporization.

58 GEOSCIENCES↗

DOE OSTI2023

Influence of clay mineral weathering on green rust formation at iron-reducing conditions

Green rusts (GR) are important drivers for trace metal and nutrient cycling in suboxic environments. We investigated whether green rusts would incorporate aluminum (Al) or other elements from naturally-formed clay minerals containing easily-weatherable clay minerals (e.g. mica, interlayered clays). We isolated the clay minerals from a Matapeake silt loam soil by removal of silt and sand, organic matter, and reducible oxides to study mechanisms of interaction between Fe(II) and soil-sourced clay minerals. We conducted batch Fe(II) sorption experiments at multiple near-neutral pHs (6.5–7.5) and reaction times (2h–365 days). Mineral transformations were characterized by selective extractions, X-ray diffraction (XRD), and Fe X-ray absorption spectroscopy (XAS) analyzed by shell-fitting and linear combination fitting (LCF) with natural and synthetic standards. Clay mineral fraction contained a mixture of quartz, kaolinite, interlayered vermiculite, mica, and chlorite with significant structural Fe (2.6% wt). Uptake of Fe(II) increased with pH and kinetics were rapid until 5 days, followed by slow continuous Fe(II) uptake. Citrate-bicarbonate desorption kinetics from Fe(II) sorbed clay released more Al and silicon (Si) compared with unreacted soil clay fraction whereas magnesium (Mg) and potassium (K) were unaffected. Citrate-bicarbonate extracted Fe contained more Fe(II) than an ideal GR with an Fe(II)/Fe(III) molar ratio of 5.50. Analysis of the Fe EXAFS by both LCF and shell fitting was best modeled as a combination of Fe(III)-clay reduction to Fe(II) and precipitation of GR and Fe(II)-Al LDH. After 7 days of Fe(II) sorption, LCF identified 55.2% total Fe in clay, 33.4% GR(Cl) and 11.4% Fe(II)-Al LDH. Further, these results provide novel evidence of Fe(II)-Al LDHs precipitating on naturally-formed soil clay minerals as a minor phase to GR. The geochemical implications are that GRs formed in soils and sediments should be considered to have Al and Si as well as Mg substitutions affecting their structure and reactivity.

58 GEOSCIENCES↗

DOE OSTI2023

Molybdenum in basalt-hosted seafloor hydrothermal systems: Experimental, theoretical, and field sampling approaches

Here, seafloor hydrothermal vents represent potential sources of Mo and other biologically relevant transition metals to the global ocean, complementing continental runoff. Here, we use a combination of experimental, theoretical, and field-sampling approaches to investigate the behavior of Mo in basalt-hosted seafloor hydrothermal systems to provide insight into the processes controlling Mo concentrations in hydrothermal fluids and to derive estimates of vent fluid Mo concentrations and fluxes. Results of this study demonstrate that reaction fluids generated from 350°C, 500 bar hydrothermal basalt alteration experiments contain 775–801 nmol/kg Mo and are thus comparable to a recently collected time series of natural seafloor vent fluids that contained 200–220 nmol/kg Mo at 302°C and 29–30 nmol/kg at 281–282°C (Evans et al., 2023). Synchrotron-based analyses of experimentally altered basalt produced in this study and additional natural samples of altered oceanic crust originally collected from Pito Deep Rift reveal the presence of Mo-rich particles consistent with trace molybdenite. Comparisons of Mo:Cu ratios in natural vent fluids and near-vent sediment trap samples from Main Endeavour Field indicate that vent fluid Mo is readily incorporated into buoyant plume particles and advected out of the near-vent field, analogous to previous mass balance studies of Cu in this region. Thermodynamic calculations of molybdenite solubility in the context of mineral-buffered hydrothermal fluids and comparisons with natural and experimental hydrothermal fluids suggest that high-temperature vent fluids contain 30–1500 nmol/kg Mo. While a minor component of the modern Mo budget, hydrothermal Mo fluxes are estimated to have constituted 0.3–200× the contemporaneous continental weathering fluxes prior to the ~2.4 Ga ago “Great Oxidation Event” and widespread oxidative continental weathering. Overall, identification of hydrothermal vents as a source of Mo-rich plume particles with potential for dispersal into the wider marine environment has significant implications for hypotheses regarding the co-evolution of Life and Earth’s environments, specifically the form and availability of Mo in anoxic Archean-Eon oceans, where Mo-dependent enzymatic pathways are thought to have emerged and subsequently evolved.

58 GEOSCIENCES↗

DOE OSTI2023

Experimental determinations of carbon and hydrogen isotope fractionations and methane clumped isotope compositions associated with ethane pyrolysis from 550 to 600°C

Methane clumped isotope compositions signify the relative natural abundances of rare, doubly substituted isotopic species of methane ( 13 CH 3 D and 12 CH 2 D 2 ) and have emerged as a new isotopic tool to trace the sources, sinks, and lifecycles of methane in the environment. Such measurements can identify equilibration (or reequilibration) temperatures if found to be in isotopic equilibrium or non-equilibrium processes (e.g., kinetically controlled reactions or mixing) if not in isotopic equilibrium. Naturally occurring thermogenic methane—formed by the thermally activated breakdown of larger organic molecules—has been found to have clumped isotope compositions consistent with equilibrium at reasonable gas formation temperatures in some settings and non-equilibrium processes occurring during either formation, migration, storage, or extraction in others. To explore the potential controls on the isotopic composition of thermogenic methane, we conducted isothermal time-series ethane pyrolysis experiments at 550 and 600 °C to measure methane and ethane 13 C/ 12 C and D/H fractionations and methane clumped isotope compositions (resolved 13 CH 3 D and 12 CH 2 D 2 ). We explore the effects of modifying the initial clumped isotope composition of ethane and the addition of water vapor to pyrolysis experiments. We observe that ethane and methane 13 C/ 12 C are controlled by kinetic isotope effects and Rayleigh distillation processes. In contrast, ethane and methane D/H and methane clumped isotope compositions appear to be controlled by a combination of these processes and hydrogen isotope exchange. The hydrogen isotope exchange processes lead to isotopic equilibrium as reaction completion is approached for both D/H (ethane/methane) and methane clumped isotope compositions. Here, we develop a chemical model based on a mass balance approach that accounts for inheritance vs. hydrogen-abstraction formation pathways for singly and doubly substituted isotopologues of ethane and methane that is compared to the experimental data. The model allows the determination of carbon and hydrogen kinetic isotope effects associated with ethane cracking and hydrogen abstraction reactions that, where applicable, we compare to prior theoretical constraints. From the comparison of the model to the experimental data, we infer that the kinetically controlled ethane and methane bulk isotope compositions and methane clumped isotope compositions are controlled by kinetic isotope effects (both primary and secondary) associated with both C–C bond and C–H bond cleavage reactions. Specifically, the methane clumped isotope compositions likely result from a combination of clumped isotope effects associated with ethane breakdown and/or assembly of methane isotopologues (expressed in terms of γ-factor parameters ≠ 1) and combinatorial effects that arise probabilistically. We discuss our experimental results in the context of recent pyrolysis experiments and observations of naturally occurring thermogenic methane. We consider a proposal consistent with observations from nature that the hydrogen isotope exchange reactions that promote equilibration of methane isotopic molecules at or near formation temperature may be facilitated by free radicals generated by pyrolysis reactions. In this framework, isotope exchange effectively ceases when pyrolysis effectively ceases locking in compositions that can be consistent with peak formation temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Antimony stable isotope fractionation during adsorption onto birnessite: A molecular perspective from X-ray absorption spectroscopy and density functional theory

Sorption of antimony (Sb) onto birnessite significantly influences the fate of Sb in oceanic and terrestrial environments and fractionates Sb isotopes. Nevertheless, little is known about Sb isotopic fractionation during its adsorption on birnessite. Here, in this study, we show the value of Δ 123 Sb adsorbed-aqueous increases from −0.398 to −0.332 ‰ in 1 h and then decreases and stabilizes at −0.384 ‰ in 72 h. The enrichment of the light Sb isotope is predominantly due to the distortion of the octahedral symmetry. X-ray absorption spectroscopy results indicate Sb first forms a double-corner-sharing complex on birnessite and then transforms to a double-edge-sharing complex during adsorption. The optimized bond distances for double-corner-sharing (3.37 Å) and double-edge-sharing (2.90 Å) complexes calculated using density functional theory (DFT) fits well with the structure (3.41 and 3.00 Å) revealed by X-ray absorption spectroscopy, respectively. The fractionation derived from reduced partition function ratios calculated using DFT aligns well with the experimental results. Therefore, the variation in Sb isotopic fractionation during adsorption is attributed to the evolving structure of Sb complexes on birnessite. Our results demonstrate the isotopic fractionation of Sb during adsorption on birnessite and provide a molecular-scale understanding of Sb behavior, contributing to the correct reconstruction of the Sb isotope composition of ancient seawater using ferromanganese crusts and nodules, and efforts to trace Sb migration in epigenetic mining environments.

Adsorption↗

DOE OSTI2025

Temperature dependence of uranium and thorium partitioning in igneous zircons

Zircon is a key mineral in geochronology because of its chemical and physical durability and tendency to incorporate radioactive trace elements such as U and Th. Quantifying the partitioning of the actinide elements is critical to constrain initial non-secular equilibrium amounts of 230 Th in zircon. An excess or deficit of 206 Pb will be produced from such an initial excess/deficit of 230 Th from the secular equilibrium condition, which influences the calculated 206 Pb/ 238 U age (Schärer, 1984, Mattinson, 1973). However, there is no standard way to calculate Th/U partitioning ratios when applying age corrections to young igneous zircon, making uncertainties hard to estimate. To better understand U and Th partitioning of zircon/melt that is the cornerstone of the secular disequilibrium correction, zircon was synthesized in one-atmosphere experiments using basaltic andesite, andesite, and rhyolite starting materials, doped with Zr, U, and Th. Different experimental temperatures and oxygen fugacity conditions (ΔQFM−4 to ΔQFM+4) were explored to examine their effects on U and Th partitioning. In addition, we specifically quantify the effects of sector zoning, fractional crystallization, and melt composition on U and Th partitioning. By combining experimental and natural zircon data, we find that temperature has the primary control on the partitioning of U and Th in the zircon and calibrate an inverse relationship between these partition coefficients and crystallization temperatures. The calibrated equation can be applied to the 230 Th correction for an improvement in the accuracy of Th-corrected 206 Pb/ 238 U dates when the zircon crystallization temperature is known.

206Pb/238U dating↗

DOE OSTI2025

Spatial evidence of cryptic methane cycling and methylotrophic metabolisms along a land–ocean transect in salt marsh sediment

Methylotrophic methanogenesis in the sulfate-rich zone of coastal and marine sediments couples with anaerobic oxidation of methane (AOM), forming the cryptic methane cycle. This study provides evidence of cryptic methane cycling in the sulfate-rich zone across a land–ocean transect of four stations–two brackish, one marine, and one hypersaline–within the Carpinteria Salt Marsh Reserve (CSMR), southern California, USA. Samples from the top 20 cm of sediment from the transect were analyzed through geochemical and molecular (16S rRNA) techniques, in-vitro methanogenesis incubations, and radiotracer incubations utilizing 35 S-SO 4 , 14 C-mono-methylamine, and 14 C-CH 4 . Sediment methane concentrations were consistently low (3 to 28 µM) at all stations, except for the marine station, where methane increased with depth reaching 665 µM. Methanogenesis from mono-methylamine was detected throughout the sediment at all stations with estimated CH 4 production rates in the sub-nanomolar to nanomolar range per cm 3 sediment and day. 16S rRNA analysis identified methanogenic archaea (Methanosarcinaceae, Methanomassiliicoccales, and Methanonatronarchaeacea) capable of producing methane from methylamines in sediment where methylotrophic methanogenesis was found to be active. Metabolomic analysis of porewater showed mono-methylamine was mostly undetectable (<3 µM) or present in trace amounts (<10 µM) suggesting rapid metabolic turnover. In-vitro methanogenesis incubations of natural sediment showed no linear methane buildup, suggesting a process limiting methane emissions. AOM activity, measured with 14 C-CH 4 , overlapped with methanogenesis from mono-methylamine activity at all stations, with rates ranging from 0.03 to 19.4 nmol cm −3 d −1 . Geochemical porewater analysis showed the CSMR sediments are rich in sulfate and iron. Porewater sulfate concentrations (9–91 mM) were non-limiting across the transect, supporting sulfate reduction activity (1.5–2,506 nmol cm −3 d −1 ). Porewater sulfide and iron (II) profiles indicated that the sediment transitioned from a predominantly iron-reducing environment at the two brackish stations to a predominantly sulfate-reducing environment at the marine and hypersaline stations, which coincided with the presence of phyla (Desulfobacterota) involved in these processes. AOM activity overlapped with sulfate reduction and porewater iron (II) concentrations suggesting that AOM is likely coupled to sulfate and possibly iron reduction at all stations. However, 16S rRNA analysis identified anaerobic methanotrophs (ANME-2) only at the marine and hypersaline stations while putative methanogens were found in sediment across all stations. In one sediment horizon at the marine station, methanogen families (Methanosarcinaceae, Methanosaetaceae, Methanomassiliicoccales, and Methanoregulaceae) and ANME 2a,2b, and 2c groups were found together. Collectively, our data suggest that at the brackish stations methanogens alone may be involved in cryptic methane cycling, while at the marine and hypersaline stations both groups may be involved in the process. Differences in rate constants from incubations with 14 C-labeled methane and mono-methylamine suggest a non-methanogenic process oxidizing mono-methylamine to inorganic carbon, likely mediated by sulfate-reducing bacteria. Understanding the potential competition of sulfate reducers with methanogens for mono-methylamine needs further investigation as it might be another important process responsible for low methane emissions in salt marshes.

Anaerobic oxidation of methane↗

DOE OSTI2022

Lithological controls on soil geochemistry and clay mineralogy across Spodosols in the coastal temperate rainforest of southeast Alaska

The interaction of similar soil forming factors with different parent materials determines soil chemical weathering, influencing soil processes and properties. In the coastal temperate rainforests of southeast Alaska, Spodosols is the dominant soil order in well-drained parent materials regardless of the lithology. Yet, the role of lithology on chemical weathering, base cation depletion, and mineral transformation during pedogenesis remains elusive. Here, we established a lithosequence comprised of soils derived from tonalite, phyllite, slate, and metavolcanic rocks to test the hypothesis that despite the soils in southeast Alaska presenting similar taxonomy and morphology, the influence of lithology in chemical weathering intensity and mineralogy can be detected. We evaluated physicochemical properties, clay mineralogy, Fe-Al oxides, chemical weathering, and elemental mass balance on eleven Spodosols sampled along the lithosequence. We also propose a new weathering index, the Weathering Index for Spodosols (WISP), suitable to evaluate the weathering degree of soils with Al mobility and base-rich parent materials. This index fills a gap in traditional weathering indices that uses Al as an immobile element and/or does not evaluate the leaching of main base cations. We found the pedons across the lithosequence expressed similar physicochemical properties, with predominantly thick profiles, presenting andic properties, relatively thick spodic horizons, and thin E horizons. Podzolization imposed similar mineral transformation trends, mainly dissolution of chlorite, mica alteration to interstratified mica-vermiculite, and formation of smectite in E horizons. Kaolinite was detected in the soils from all lithologies, but only in trace amounts in phyllite soils. Fe oxides depth distribution was similar but with a higher concentration in the slate (44.5 ± 4.2 g/kg), followed by metavolcanic (31.6 ± 4.0 g/kg), phyllite (27.9 ± 3.8 g/kg), and tonalite soils (19.9 ± 5.7 g/kg) and predominantly as organometallic complexes (extracted by Na-pyrophosphate). The stronger weathering of chlorite and mica in the E horizons was reflected in larger losses of Mg and K, while Ca and Na were more depleted in the spodic horizons, suggesting plagioclase weathering. Slate soils were the most depleted in base cations, followed by phyllite, tonalite, and metavolcanic soils. The more depleted status of slate soils was supported by the higher WISP values (WISP = 62), followed by phyllite (WISP = 55), tonalite (WISP = 53), and metavolcanic soils (WISP = 45). Our results demonstrated that lithology controlled elemental depletion intensities, Fe oxide concentrations, and partially the clay fraction mineralogy in addition to podzolization acting as a dominant pedogenic process across all sites. These findings advance our understanding of the role of lithology on soil mineralogy/geochemistry that impacts critical soil functions, such as the soil carbon cycle, elemental fluxes, organo-mineral interactions, and solid-solution reactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Magmatic, hydrothermal and ore element transfer processes of the southeastern Archean Superior Province implied from electrical resistivity structure

Although magmatic- and metamorphic- derived fluids are widely recognized mineralizing agents, the role of crustal architecture in defining source and sink zones within the middle to lower crust and upper mantle of ancient orogens remains enigmatic. The globally largest and best-preserved Archean greenstone belts lie in the Superior Province, Canada. They provide an ideal location to investigate the influence of igneous construction and subsequent syn-deformational plutonism and metamorphism on the localization of metal-rich melts and fluids throughout the crustal column. Integration of three-dimensional magnetotelluric modelling and seismic reflection sections across the Abitibi subprovince reveals details of a 'whole-of-crust' magmatic and hydrothermal system. East-west low resistivity structures broadly underlie the surface traces of the major deformation zones that are host to significant gold endowment (>200 Moz), while mid-crustal cross trends suggest mineralized fluids flowed along 'pipes' within fault planes. Most low resistivity structures are inferred to represent domains containing interconnected zones of graphite and/or sulfide. Furthermore, these delineate relict mantle source/transit domains and crustal pathways enriched by the flow of magmas or metamorphic fluids genetically related to a late-stage pulse of ore-bearing magmatism, possibly as a result of slab break-off or delamination. Using the combined MT and seismic data, we develop a 3-D crustal-scale model which highlights how evolving orogenic architecture-controlled mass transfer and metallogenic processes developed.

58 GEOSCIENCES↗

DOE OSTI2022

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A numerical study of the spectral radiative properties of packed bed with mixed bauxite and silica spheres

Bauxite and silica particles have gained increasing attention for applications in the field of concentrated solar power. In this work, a Monte Carlo ray-tracing simulation is performed to predict the radiative properties (absorptance, reflectance, and transmittance) of packed beds with mixed bauxite and silica spherical particles at wavelengths of 0.5, 2.6, 9.0, and 9.35 µm. These wavelengths are representative for the visible, near-infrared, and the mid-infrared regions that are important for solar and thermal radiation, respectively. A repeating unit-column approach is used to mathematically represent the particle bed. The effects of particle mixing ratio, volume fraction, and wavelength on the predicted radiative properties are examined. The obtained radiative properties are inputted to an inverse method to retrieve the effective absorption and scattering coefficients as well as the scattering albedo, which may be later used in a continuous-scale radiative heat transfer analysis. Furthermore, the independent scattering model is used to obtain the absorption coefficient and scattering albedo based on the absorption and scattering cross sections predicted by a Monte Carlo algorithm for a single particle. It is shown that the independent scattering model underpredicts the scattering coefficient for opaque particles but overpredicts the scattering coefficient for semitransparent particles for sufficiently high particle volume fractions. Here, the radiative properties calculated from the independent scattering model are compared to the full Monte Carlo simulation of the particle bed to examine the influence of particle mixing on dependent scattering.

42 ENGINEERING↗

DOE OSTI2025

Accurate and efficient predictions of keyhole dynamics in laser materials processing using machine learning-aided simulations

The keyhole phenomenon has been widely observed in laser materials processing, including laser welding, remelting, cladding, drilling, and additive manufacturing. Keyhole-induced defects, primarily pores, dramatically affect the performance of final products, impeding the broad use of these laser-based technologies. The formation of these pores is typically associated with the dynamic behavior of the keyhole. So far, the accurate characterization and prediction of keyhole features, particularly keyhole depth, as a function of time, has been a challenging task. In situ characterization of keyhole dynamic behavior using the synchrotron X-ray technique is informative but complicated and expensive. Current simulations are generally hindered by their poor accuracy and generalization abilities in predicting keyhole depths due to the lack of accurate laser absorptance data. In this study, we develop a machine learning-aided simulation method that accurately predicts keyhole dynamics, especially in keyhole depth fluctuations, over a wide range of processing parameters. In two case studies involving titanium and aluminum alloys, we achieve keyhole depth prediction with a mean absolute percentage error of 10 %, surpassing those simulated using the ray-tracing method with an error margin of 30 %, while also reducing computational time. This exceptional fidelity and efficiency empower our model to serve as a cost-effective alternative to synchrotron experiments. Our machine learning-aided simulation method is affordable and readily deployable for a large variety of materials, opening new doors to eliminate or reduce defects for a wide range of laser materials processing techniques.

Computational fluid dynamics↗

DOE OSTI2021

Structures of two novel crystal forms of Aspergillus oryzae alpha amylase (taka-amylase)

The structures of Aspergillus oryzae α-amylase were determined in a tetragonal crystal, having one molecule as asymmetric unit, and a monoclinic crystal with two molecules as asymmetric unit. Both crystal forms were obtained from trace contaminants of an old commercial lipase preparation. Structures were determined and refined to 1.65 Å and 1.43 Å resolution respectively. The latter crystal has a non-crystallographic (NCS) twofold axis within the asymmetric unit. Glycosylation at Asn197 is evident, and in the tetragonal crystal can be seen to include three, partially disordered sugar residues following the initial N-acetyl glucosamine (NAG). Superposition of the tetragonal crystal model on the α-amylases from Bacillus subtilis (PDB:1BAG), pig pancreas (PDB:3L2L), and barley (PDB:1AMY), show a high degree of coincidence, particularly for the (β/α)8-barrel domains, and especially within the active site. Using this structural agreement between amylases, we extrapolated the binding model of a six residue, limit dextrin found in pig pancreas α-amylase to the A. oryzae enzyme model, which predicts substrate interacting amino acid residues.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

La 3+ and Y 3+ interactions with the carboxylic acid moiety at the liquid/vapor interface: Identification of binding complexes, charge reversal, and detection limits

Specific interactions of yttrium and lanthanum ions with a fatty acid Langmuir monolayer were investigated using vibrational sum frequency spectroscopy. The trivalent ions were shown to interact with the charged form of the carboxylic acid group from nanomolar concentrations (<300 nM). Analysis of the spectral features from both the symmetric and the asymmetric carboxylate modes reveals the presence of at least three distinct coordination structures linked to specific binding configurations. Although the same species were identified for both La 3+ and Y 3+ , they display a different concentration dependence, highlighting the ion-specificity of the interaction. From the analysis of the response of interfacial water molecules, the reversal of the surface charge, as well as the formation of yttrium hydroxide complexes, were detected upon increasing the amount of salt in solution. The binding interaction and kinetics of absorption are sensitive to the solution pH, showing a distinct ion speciation in the interfacial region when compared to the bulk. Changing the subphase pH or adding a monovalent background electrolyte that promotes deprotonation of the carboxylic acid headgroup could further improve the detection limit of La 3+ and Y 3+ to concentrations < 100 nM. These findings demonstrate that nM concentrations of trace metals contaminants, typically found on monovalent salts, can significantly influence the binding structure and kinetics in Langmuir monolayers.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2021

A novel energy balance approach for a verifiable and accurate solution of radiation extinction in purely absorbing particle clouds

Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.

42 ENGINEERING↗

DOE OSTI2024

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

X-ray topography characterization of gallium nitride substrates for power device development

Gallium nitride substrates grown by the hydride vapor phase epitaxy (HVPE) method using a patterned growth process have been characterized by synchrotron monochromatic beam X-ray topography in the grazing incidence geometry. Images reveal a starkly heterogeneous distribution of dislocations with areas as large as 0 3 mm 2 containing threading dislocation densities below 10 3 cm -2 in between a grid of strain centers with higher threading dislocation densities (> 10 4 cm -2 ). Basal plane dislocation densities in these areas are as low as 104 cm -2 . By comparing the recorded images of dislocations with ray tracing simulations of expected dislocations in GaN, the Burgers vectors of the dislocations have been determined. The distribution of threading screw/mixed dislocations (TSDs/TMDs), threading edge dislocations (TEDs) and basal plane dislocations (BPDs) is discussed with implications for fabrication of power devices.

36 MATERIALS SCIENCE↗