Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential transformation”

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

Spurious solar-wind effects on acceleration noise in LISA Pathfinder

Spurious solar-wind effects are a potential noise source in future Laser Interferometer Space Antenna (LISA) measurements. One noise coupling mechanism is constrained by estimating solar-wind effects on acceleration noise in LISA Pathfinder (LPF). While LISA is designed for drag-free differential measurement, predicting the realistic impact both bounds the operational environment and assesses whether LISA could provide serendipitous space-weather observations. Data from NASA's Advanced Composition Explorer (ACE), situated at the L1 Lagrange point, serves as a reliable source of solar-wind data. The data sets are compared over the 114 d time period from 1 March 2016 to 23 June 2016. This period gives the longest readily-available open data set, without interference from other commissioning activities. To evaluate space weather effects, the data from both satellites are formatted, gap-filled/interpolated, and fast-Fourier transformed for amplitude spectral density and coherence comparisons. Solar wind effects are not seen in a coherence plot between LPF and ACE; modest coherence in the planned LISA observational frequency band can be attributed to chance. This result indicates that measurable correlation due to solar-wind acceleration noise over 3 month timescales will be a negligible noise source. LISA is unlikely to inform solar wind measurements routinely. Another source of noise from the Sun, solar radiation pressure, is estimated to impart greater acceleration noise, but has yet to be analyzed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

GrainGNN: A dynamic graph neural network for predicting 3D grain microstructure

We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.

36 MATERIALS SCIENCE↗
DOE OSTI2020

MEUMAPPS (Microstructure Evolution Using Massively Parallel Phase-field Simulations)

The software “MEUMAPPS” is a h high-performance computing code used to simulate microstructure evolution associated with diffusional solid-state transformations in structural alloys. Understanding microstructure evolution during thermo-mechanical processing of structural alloys is the first step towards designing and processing alloys for specific applications by meeting property requirements demanded by the application. In this instance, the code is an important component in predicting processing-structure linkages during additive manufacturing of structural alloys as a function of processing parameters and alloy composition used in an additive manufacturing process. The code provides a detailed three-dimensional distribution of different types of phases / constituents and their morphologies, as well as the compositions of these phases that make up the microstructure. The microstructure is simulated by solving the governing partial differential equations using a Fourier Spectral Method.

Radhakrishnan, Balasubram↗
DOE OSTI2020

Reprint of "Selective conversion of CO into ethanol on Cu(511) surface reconstructed from Cu(pc): Operando studies by electrochemical scanning tunneling microscopy, mass spectrometry, quartz crystal nanobalance, and infrared spectroscopy"

A polycrystalline copper, surface-terminated by a well-defined (511)-oriented facet, was electrochemically generated by a series of step-wise surface reconstruction and iterations of mild oxidative-reductive processes in 0.1 M KOH. The electrochemical reduction of CO on the resultant stepped surface was investigated by four surface-sensitive operando methodologies: electrochemical scanning tunneling microscopy (STM), electrochemical quartz crystal nanobalance (EQCN), differential electrochemical mass spectrometry (DEMS), and polarization-modulation infrared spectroscopy (PMIRS). The stepped surface catalyzed the facile conversion of CO into ethanol, the exclusive alcohol product at a low overpotential of -1.06 V (SHE) or - 0.3 V (RHE). Additionally, the chemisorption of CO was found to be a necessary prelude to ethanol production; i.e. the surface coverages, rather than solution concentrations, of CO and its surface-bound intermediates primarily dictate the reaction rates (current densities). Contrary to the expected predominance of undercoordinated step-site reactivity over the coordination chemistry of vicinal surfaces, vibrational spectroscopic evidence reveals the involvement of terrace-bound CO adsorbates during the multi-atomic transformations associated with the production of ethanol.

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

In-vivo Raman microspectroscopy reveals differential nitrate concentration in different developmental zones in Arabidopsis roots

Abstract Background Nitrate (NO 3 − ) is one of the two major forms of inorganic nitrogen absorbed by plant roots, and the tissue nitrate concentration in roots is considered important for optimizing developmental programs. Technologies to quantify the expression levels of nitrate transporters and assimilating enzymes at the cellular level have improved drastically in the past decade. However, a technological gap remains for detecting nitrate at a high spatial resolution. Using extraction-based methods, it is challenging to reliably estimate nitrate concentration from a small volume of cells (i.e., with high spatial resolution), since targeting a small or specific group of cells is physically difficult. Alternatively, nitrate detection with microelectrodes offers subcellular resolution with high cell specificity, but this method has some limitations on cell accessibility and detection speed. Finally, optical nitrate biosensors have very good ( in-vivo ) sensitivity (below 1 mM) and cellular-level spatial resolution, but require plant transformation, limiting their applicability. In this work, we apply Raman microspectroscopy for high-dynamic range in-vivo mapping of nitrate in different developmental zones of Arabidopsis thaliana roots in-situ . Results As a proof of concept, we have used Raman microspectroscopy for in-vivo mapping of nitrate content in roots of Arabidopsis seedlings grown on agar media with different nitrate concentrations. Our results revealed that the root nitrate concentration increases gradually from the meristematic zone (~ 250 µm from the root cap) to the maturation zone (~ 3 mm from the root cap) in roots grown under typical growth conditions used for Arabidopsis, a trend that has not been previously reported. This trend was observed for plants grown in agar media with different nitrate concentrations (0.5–10 mM). These results were validated through destructive measurement of nitrate concentration. Conclusions We present a methodology based on Raman microspectroscopy for in-vivo label-free mapping of nitrate within small root tissue volumes in Arabidopsis. Measurements are done in-situ without additional sample preparation. Our measurements revealed nitrate concentration changes from lower to higher concentration from tip to mature root tissue. Accumulation of nitrate in the maturation zone tissue shows a saturation behavior. The presented Raman-based approach allows for in-situ non-destructive measurements of Raman-active compounds.

Fernández González, Alma↗
DOE OSTI2023

Multidimensional quantitative phenotypic and molecular analysis reveals neomorphic behaviors of p53 missense mutants

Abstract Mutations in the TP53 tumor suppressor gene occur in >80% of the triple-negative or basal-like breast cancer. To test whether neomorphic functions of specific TP53 missense mutations contribute to phenotypic heterogeneity, we characterized phenotypes of non-transformed MCF10A-derived cell lines expressing the ten most common missense mutant p53 proteins and observed a wide spectrum of phenotypic changes in cell survival, resistance to apoptosis and anoikis, cell migration, invasion and 3D mammosphere architecture. The p53 mutants R248W, R273C, R248Q, and Y220C are the most aggressive while G245S and Y234C are the least, which correlates with survival rates of basal-like breast cancer patients. Interestingly, a crucial amino acid difference at one position—R273C vs. R273H—has drastic changes on cellular phenotype. RNA-Seq and ChIP-Seq analyses show distinct DNA binding properties of different p53 mutants, yielding heterogeneous transcriptomics profiles, and MD simulation provided structural basis of differential DNA binding of different p53 mutants. Integrative statistical and machine-learning-based pathway analysis on gene expression profiles with phenotype vectors across the mutant cell lines identifies quantitative association of multiple pathways including the Hippo/YAP/TAZ pathway with phenotypic aggressiveness. Further, comparative analyses of large transcriptomics datasets on breast cancer cell lines and tumors suggest that dysregulation of the Hippo/YAP/TAZ pathway plays a key role in driving the cellular phenotypes towards basal-like in the presence of more aggressive p53 mutants. Overall, our study describes distinct gain-of-function impacts on protein functions, transcriptional profiles, and cellular behaviors of different p53 missense mutants, which contribute to clinical phenotypic heterogeneity of triple-negative breast tumors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Position-specific carbon isotopes of Murchison amino acids elucidate extraterrestrial abiotic organic synthesis networks

The Murchison meteorite is a well-studied carbonaceous chondrite with relatively high concentrations of amino acids thought to be endogenous to the meteorite, in part because they are characterized by carbon isotope (δ 13 C) values higher than those typical of terrestrial amino acids. Past studies have proposed that extraterrestrial amino acids in the Murchison meteorite could have formed by Strecker synthesis (for α-amino acids), Michael addition (for β-amino acids), or reductive amination, but a lack of constraints have prevented confident discrimination among these possibilities, or assignment of specific formation pathways to each of several specific amino acids. Position-specific carbon isotope analysis differentiates amongst these mechanisms by relating molecular sites to isotopically distinct carbon sources and by constraining isotope effects associated with elementary chemical reactions. Prior measurements of the position-specific carbon isotopic composition of α-alanine from the Murchison CM chondrite demonstrated that alanine’s high δ 13 C VPDB value is attributable to the amine carbon (δ 13 C VPDB = +142 ± 20‰), consistent with Strecker synthesis drawing on 13 C-rich carbonyl groups in precursors (Chimiak et al., 2021). Here, we measured the δ 13 C composition of fragment ions generated by electron impact ionization of derivatized α-alanine, β-alanine, and aspartic acid from Murchison via gas chromatography-Fourier transform mass spectrometry. α-Alanine’s amine carbon yielded δ 13 C VPDB = +109 ± 21‰, which is consistent with the previously measured value and with formation from 13 C-rich precursors. β-Alanine’s amine carbon presents a lower δ 13 C VPDB = +33 ± 24‰, which supports formation from 13 C-rich precursors but potentially via a Michael addition mechanism rather than Strecker synthesis. Aspartic acid’s amine carbon has δ 13 C VPDB = -14 ± 5‰, suggesting synthesis from precursors distinct from those that generated the alanine isomers. Further, these measurements indicate that Murchison amino acids are a mixture of compounds made from different synthesis mechanisms, though some subsets likely drew on the same substrates; this conclusion highlights the complexity of extraterrestrial organic synthesis networks and the potential of emerging methods of isotope ratio analysis to elucidate the details of those networks.

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

Random field optimization

Herein we present a new modeling paradigm for optimization that we call random field optimization. Random fields are a powerful modeling abstraction that aims to capture the behavior of random variables that live on infinite-dimensional spaces (e.g., space and time) such as stochastic processes (e.g., time series, Gaussian processes, and Markov processes), random matrices, and random spatial fields. This paradigm involves sophisticated mathematical objects (e.g., stochastic differential equations and space-time kernel functions) and has been widely used in neuroscience, geoscience, physics, civil engineering, and computer graphics. Despite of this, however, random fields have seen limited use in optimization; specifically, existing optimization paradigms that involve uncertainty (e.g., stochastic programming and robust optimization) mostly focus on the use of finite random variables. This trend is rapidly changing with the advent of statistical optimization (e.g., Bayesian optimization) and multi-scale optimization (e.g., integration of molecular sciences and process engineering). Our work extends a recently-proposed abstraction for infinite-dimensional optimization problems by capturing more general uncertainty representations. Moreover, we discuss solution paradigms for this new class of problems based on finite transformations and sampling, and identify open questions and challenges.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Heat capacity of microgram oxide samples by fast scanning calorimetry

Quantitative scanning calorimetry on microgram-sized samples opens a broad, new range of opportunities for studying the thermodynamic properties of quantity-limited materials, including those produced under extreme conditions or found as rare accessory minerals in nature. We calibrated the Mettler Toledo Flash DSC 2+ calorimeter to obtain quantitative heat capacities in the range 200–350 °C, using samples weighing between 2 and 11.5 μg. Our technique is applied to a new set of oxide materials to which it has never been used before, without the need for melting, glass transitions, or phase transformations. Heat capacity data were obtained for silica in the high pressure stishovite (rutile) structure, dense post-stishovite glass, standard fused quartz, and for TiO2 rutile. These heat capacities agree within 5%–15% with the literature values reported for rutile, stishovite, and fused SiO2 glass. The heat capacity of post-stishovite glass, made by heating stishovite to 1000 °C, is a newly reported value. After accurate calibrations, measured heat capacities were then used to calculate masses for samples in the microgram range, a substantial improvement over measurement in conventional microbalances, which have uncertainties approaching 50%–100% for such small samples. Since the typical uncertainty of heat capacities measured on 10–100 mg samples in conventional differential scanning calorimetry is typically 7% (1%–5% with careful work), flash differential scanning calorimetry, using samples a factor of 1000 smaller, increases the uncertainty of heat capacity measurements by a factor of <3, opening the door for meaningful measurements on ultra-small, high-pressure samples and other quantity-limited materials.

Instruments & Instrumentation↗
DOE OSTI2021

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗
DOE OSTI2025

Strain-driven oxygen vacancy ordering in LaNiO 3 thin films revealed by integrated differential phase contrast imaging in scanning transmission electron microscopy

Rare-earth nickelates, such as LaNiO 3 (LNO), exhibit complex electronic properties, with ordered oxygen vacancies (OOV) influencing conductivity and magnetic behavior. We investigate the structural stability of strain-induced OOV phases in LNO thin films grown on SrTiO 3 substrates and the impact of Ruddlesden–Popper (RP) faults. Using high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and integrated differential phase contrast (iDPC) STEM imaging, we conducted atomic-scale structural and compositional analyses of OOV. Geometric phase analysis (GPA) was employed to measure the strain in fault-free and RP fault regions, while density functional theory (DFT) calculations explored different OOV arrangements in the LNO phase. Simulated iDPC-STEM imaging of energy-stabilized structures was performed to correlate with experimental results. Here, our findings reveal superstructure modulation in the chemical composition and atomic-scale lattice structure in LNO, primarily due to the formation of the OOV in Ni–O layers of the LaNiO 2.5 phase. The out-of-plane compressive strain of about 2% stabilizes this phase, reducing the strain, diminishing OOV, and transforming them into LNO.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bark morphological and chemical features are differentially correlated with disease resistance and yield in hybrid poplar taxa

In the southeastern United States, the establishment of short-rotation intensively cultured plantations of hybrid poplar has been hindered by its susceptibility to stem cankers. We evaluated the tradeoffs between biomass yield and disease tolerance in hybrid poplar genotypes belonging to P. deltoides × P. maximowiczii (DM), P. deltoides × P. nigra (DN), P. trichocarpa × P. maximowiczii (TM), and P. deltoides × P. deltoides (DD) taxa. We hypothesized that canker resistant genotypes will have thicker bark but bark thickness and biomass yield will be negatively correlated. After two growing seasons, the DD genotypes developed thicker bark compared to the genotypes of other taxa and bark thickness was not correlated with biomass yield in the DD genotypes (R 2 = 0.002). However, in the TM, DM, and DN genotypes, bark thickness was negatively correlated with biomass yield (R 2 = 0.33–0.77). Disease incidence studies revealed that the DM genotypes were most susceptible to canker whereas no disease was detected in DD genotypes. Furthermore, bark analysis conducted by Fourier transform infrared spectroscopy coupled with multivariate analysis showed that that DD genotypes to be chemically separate from the three hybrid genotypes and that bark chemistry was correlated with canker disease incidence. Taken together, these results reveal that it is possible to generate hybrid poplar genotypes with thicker bark, disease resistance, and higher biomass yields. This insight should guide further efforts to develop genetically improved hybrid poplar genotypes, both in terms of biomass yield and disease tolerance, for cultivation in the southeastern United States. Hybrid poplar cultivation in southeastern United States is hindered by its susceptibility to stem cankers. We evaluated tradeoffs between yield and canker disease resistance in various hybrid poplar genotypes. After two growing seasons, the DD genotypes showed disease resistance and developed thicker bark that was chemically distinct from the other genotypes. Bark thickness was not correlated with yield in the DD genotypes but was negatively correlated with yield in the other genotypes. These results will guide the development of hybrid poplar genotypes that are both disease resistant and high yielding for cultivation in the southeastern United States.

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

Numerical fluid dynamics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. II. Entropy production and irreversibility of RG flows

Herein we demonstrate that the reformulation of renormalization group (RG) flow equations as nonlinear heat equations has severe implications on the understanding of RG flows in general. We demonstrate by explicitly constructing an entropy function for a zero-dimensional Z 2 -symmetric model that the dissipative character of generic nonlinear diffusion equations is also hard-coded in the functional RG equation. This renders RG flows manifestly irreversible, revealing the semigroup property of RG transformations on the level of the flow equation itself. Additionally, we argue that the dissipative character of RG flows, its irreversibility and the entropy production during the RG flow may be linked to the existence of a so-called C– / A-function. In total, this introduces an asymmetry in the so-called RG time—in complete analogy to the thermodynamic arrow of time—and allows for an interpretation of infrared actions as equilibrium solutions of dissipative RG flows equations. The impossibility of resolving microphysics from macrophysics is evident in this framework. Furthermore, we directly link the irreversibility and the entropy production in RG flows to an explicit numerical entropy production, which is manifest in diffusive and non-linear partial differential equations (PDEs) and a standard mathematical tool for the analysis of PDEs. Using exactly solvable zero-dimensional Z 2 -symmetric models, we explicitly compute the (numerical) entropy production related to the total variation nonincreasing property of the PDE during RG flows toward the infrared limit. Finally, we discuss generalizations of our findings and relations to the C– / A-theorem as well as how our work may help to construct truncations of RG flow equations in the future, including numerically stable schemes for solving the corresponding PDEs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Targeted genetic manipulation and yeast-like evolutionary genomics in the green alga Auxenochlorella

Auxenochlorella spp. are diploid oleaginous green algae whose streamlined genomes can be readily manipulated by homologous recombination, making them highly amenable to discovery research and bioengineering. Vegetatively diploid organisms experience specific evolutionary phenomena, including allodiploid hybridization, mitotic recombination, loss-of-heterozygosity, and aneuploidy; however, studies of these forces have largely focused on yeasts. Here, we present a telomere-to-telomere phased diploid genome assembly of Auxenochlorella UTEX 250-A (haploid length 22 Mb) and introduce a genetic toolkit for site-specific manipulation of the nuclear genome in multiple strains, featuring several selectable markers, inducible promoters, and fluorescent reporters for protein localization. UTEX 250-A is an allodiploid hybrid of Auxenochlorella protothecoides and Auxenochlorella symbiontica, two species differentiated by extensive chromosomal rearrangements. UTEX 250-A haplotypes are a mosaic of each parental species following mitotic recombination, and two chromosomes are trisomic. Loss-of-heterozygosity events are pervasive across Auxenochlorella and can evolve rapidly in the laboratory. High-quality structural annotation yielded ∼7,500 genes per haplotype. Auxenochlorella have experienced gene family loss and reduction, including core photosynthesis genes, and exhibit periodic adenine and cytosine methylation at promoters and gene bodies, respectively. Approximately 10% of genes, especially those involved in DNA repair and sex, overlap antisense long noncoding RNAs, which may participate in a regulatory mechanism. We demonstrate the utility of Auxenochlorella for fundamental research by knockout of a chlorophyll biosynthesis enzyme, and confirm one trisomy by allele-specific transformation. These results demonstrate the generality of several evolutionary forces associated with vegetative diploidy and provide a foundation for the use of Auxenochlorella as a reference organism.

CHL27↗
DOE OSTI2022

Multi-state Catalysts Modulated by Mechanical Force (Final Report)

The development of more efficient catalytic processes and new approaches to control catalytic activity and selectivity are central to the realization of more selective, atom economic, and energy efficient routes to value added chemicals and polymers. The reactivity and selectivity of a transition metal catalyst is intimately related to the ligand-sphere geometry and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another, resulting in sub-optimal efficiency. Macroscopic mechanical forces are both large, potentially much larger than interatomic forces, and are directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. As such, mechanical force represents a heretofore untapped approach to modulate catalyst geometry, with the potential to reversibly modulate catalyst geometry on the timescale of catalytic turnover or monomer enchainment. This project has addressed the fundamental challenges in material-to-molecule strain coupling associated with the development of a new class of mechanically responsive catalysts (mechanocatalysts) in which active organotransition metal catalysts are strategically embedded in a flexible polymer network such that application of external mechanical force (stretching or deformation) leads to modulation of catalyst geometry, and hence reactivity and selectivity. Our efforts during the tenure of this grant were directed toward the elucidation of force-reactivity relationships of elementary transformations that occur within the first coordination sphere of a transition metal complex employing stiff stilbene photoswitches tethered to a flexible bidentate phosphine ligand derived from MeOBiphep as molecular force probes which provide a range of compressive and extension forces to the coupled transition metal complex depending on the geometry of the stiff stilbene and length of the tethering chains. During the tenure of this grant, we have quantified the rate of C(sp 2 )-C(sp 2 ) reductive elimination from platinum(II) diaryl complexes containing bis(phosphine) force probe ligands as a function of mechanical force; compressive forces decreased the rate of reductive elimination whereas extension forces increased the rate relative to the strain-free MeOBiphep complex with a 3.4-fold change in rate over a ~290 pN range of restoring forces. In a similar manner, we have quantified the rate of oxidative addition of bromobenzene to low-ligated palladium(0) complexes containing force probe ligands as a function of mechanical force; compressive forces increase the rate of oxidative addition, whereas tensile forces decrease the rate with a ~6 fold change in rate across ~340 pN of force applied to the complexes. In both cases, experimental and computational analyses argue strongly against any significant force-induced perturbation of ground state geometry within the first coordination sphere of the reactant complexes. Rather, the force/rate behavior observed for these transformations across these ranges of forces is attributed to the coupling of force to the nuclear motion comprising the reaction coordinates for reductive elimination and oxidative addition. These results together inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle and suggest opportunities to experimentally map geometry changes associated with reactions in transition metal complexes and potential strategies for force-modulated catalysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Low temperature tetragonal polymorph of CaZrF 6

A new tetragonal polymorph of CaZrF 6 can be prepared by high energy ball milling of a CaF 2 /ZrF 4 mixture, followed by heat treatment at 325 °C. This polymorph is thermodynamically stable with respect to the well-known cubic form at low temperatures. However, it readily transforms to the cubic form on heating above ~400 °C. The tetragonal (β) CaZrF 6 is not isostructural with any previously known alkaline earth AZrF 6 phase. Unlike the cubic form, which shows strong negative thermal expansion over a wide temperature range, the tetragonal form displays positive thermal expansion in all directions (100–400 K: α l ~ +17 × 10 –6 K –1 and +13 × 10 –6 K –1 along the a- and c-axes, respectively).

36 MATERIALS SCIENCE↗