Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complexity bounding”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

2023 American Society for Mass Spectrometry (ASMS) 71st Annual Conference on Mass Spectrometry and Allied Topics

Introduction (120 words max) Understanding metal-cluster chemistry occuring at solvent boundaries in the aqueous and organic phases has applications in environments from cellular processes to nuclear fuel reprocessing. Transport of metal ions at the boundaries from aqueous to organic phases involves forming a metal-ligand complex, and wherever the initial metal coordination environment is significantly different from the final one, the metal transitions through a series of transient species in passing from one phase to another. Here an investigation of the role of coordination in the chemistry of the transient species using gas-phase measurements that are free of solvent effects to better understand the binding of complexes of metals with triphenylphosphine chalcogenide ligands, examining metal-ligand homo- and hetero-dimers to better understand transient species. Methods (120 word max) Mass spectrometry and collision induced dissociation (CID) experiments were performed with a Bruker (Billerica, MA, USA) micrOTOF-Q II quadrupole time-of-flight mass spectrometer (QTOF) and Bruker amaZon speed ETD (ion trap). High resolution/high mass accuracy spectra were generated using the QTOF. External calibration was performed with Agilent (Santa Clara, CA, USA) ESI-L Low Concentration tuning mix. Both mass spectrometers were equipped with either the electrospray ionization source or nanospray sources. Metal samples were prepared between 40 – 60 uM of the metal-ligand complex in 25% water and 75% acetonitrile. Metal-ion clusters were isolated and subjected to collision induced dissociation. Density functional theory calculations were performed. Preliminary Data or Plenary Speakers Abstract (300 words max) Metal ion clusters with triphenylphosphine chalcogenide ligands were observed for group I metals with triphenylphosphine chalcogenide samples in the mass spectrum upon electrospray ionization. For each metal ligand complex of interest, the parent ion was isolated and collision induced dissociation fragmentation spectra were acquired. We observed clusters of group I metal with triphenylphosphine oxide, triphenyl phosphine sulfide, and triphenylphosphine selenide, with homodimers and heterodimer formation. In samples where the ligands were mixed, we observed mixed sodium ligand clusters at varying amounts. These mixed ligand clusters were fragmented. Metal clusters of mixed ligand dimers containing triphenylphosphine oxide showed preferential loss of the other ligand, either triphenylphosphine sulfide or triphenyl selenide. In samples with mixed triphenylphosphine sulfide and triphenylphosphine selenium ligands, sodium bound similarly between the ligands, and losses were more evenly split, showing loss ratio upon CID with losses of triphenylphosphine sulfide:triphenylphosphine selenide 43:57 ratio observed on CID. These results suggest that the oxide binds significantly more strongly than either the selenium or sulfur triphenylphosphine ligand, and the sulfur and selenium ligands are more evenly bound. Calculations were performed using density functional theory to calculate likely structures and bond energies between the group I metal and the ligands. Novel Aspect Novel analysis of sodium bound dimers with chalcogenide triphenylphosphine ligands were investigated using mass spectrometry and theoretical calculations.

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

Selective Detection of Diglycolamide Diastereomers By Metal-Ligand Cluster Formation

Introduction (120 words): Differentiation of isomers using mass spectrometry is a challenge, as these species tend to produce ions of identical mass-to-charge ratio. Therefore, separation with mass spectrometry must rely on chemical differences between isomeric species. Separation of stereoisomers is especially challenging, as these species often have very similar chemistry. However, recent investigations of the solvent extraction behavior of diastereomers of diglycolamide molecules have shown that different stereoisomers can have significantly different lanthanide extraction efficiencies. These studies showed that the R,S stereoisomer of 2,2'-oxybis(N,N-dioctylpropanamide efficiently extracts lanthanides, while the S,S stereoisomer does not. Thus, gas-phase metal complexation may provide a pathway for differentiating these species using mass spectrometry. Methods (120 words): Europium, samarium, and holmium clusters with the R,S, and the S,S diastereomers of 2,2'-oxybis(N,N-didecylpropanamide (trivially tetradecyl diglycolamide, or mTDDGA) were generated in gas-phase using the electrospray ionization source of a Bruker (Billerica, MA, USA) micrOTOF-Q II quadrupole time-of-flight mass spectrometer. Solutions of aqueous lanthanide nitrate and diglycolamide ligands in 2-propanol were added to 50%:50% by volume acetonitrile:2-propanol to generate the spray solutions. Ligand concentration was maintained at approximately 6 µM, while lanthanide concentration was maintained at either 6 µM or 30 µM . High resolution, high mass accuracy spectra were recorded using the time-of-flight. Mass accuracy was ensured using external calibration with Agilent (Santa Clara, CA, USA) ESI-L Low Concentration tuning mix. Preliminary data (300 words): Electrospray of solutions containing only R,S-mTDDGA or S,S-m-TDDGA yield abundant peaks at m/z= 721.7, corresponding to [mTDDGA+H]+1, m/z = 743.7, corresponding to [mTDDGA+Na]+1, and m/z = 1464.4, corresponding to the sodium-bound dimer, [Na(mTDDGA)2]+1. Spray of solutions containing 6 µM europium (Eu) with 3 µM R,S-mTDDGA produced an abundant cluster ion with an isotopic envelope centered at m/z=772.0, corresponding to [Eu(R,S-mTDDGA)3]3+. Europium-containing cluster ions with isotopic envelopes apexing at m/z = 1012.3 and 1189.0 correspond to the [Eu(R,S-mTDDGA)4]3+ and [Eu(R,S-mTDDGA)5]3+ cluster ions, respectively. Protonated and sodiated R,S-mTDDGA are present as well, along with the sodium-bound dimer. In contrast, spray of solutions containing 6 µM europium (Eu) with 3 µM S,S-mTDDGA diastereomer yielded only protonated and sodiated ligand ions; no Eu- S,S-mTDDGA ions were present. This mirrors the solution-phase behavior, where the S,S diestereomer extracts lanthanides with significantly less efficiency than the R,S diestereomer. Experiments with holmium and samarium yielded the same results. Increasing the lanthanide concentration to 30 µM resulted in limited formation of lanthanide-S,S-mTDDGA clusters. However, the abundance of the [Ln(S,S-mTDDGA)3]3+ clusters remained low compared to the abundance of the protonated and sodiated ligand. Additionally, the abundance of the [Ln(R,S-mTDDGA)3]3+ cluster ion was significantly higher than that of the protonated and sodiated ligand in experiments with spray solutions containing 30 µM lanthanide with R,S-mTDDA, supporting the hypothesis that R,S-mTDDGA has a higher intrinsic affinity for the lanthanides than S,S-mTDDGA. Thus, the intrinsic complexation with lanthanides appears to offer a route to differentiate diglycolamide diestereomers in the gas-phase. Novelty (20): Stereoisomers can be differentiated in the gas phase by exploitation of differences in metal complexation chemistry.

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

LibraryX: A Framework for Cross-Library-Call Optimization

Scientific applications utilize performance libraries as a software engineering concept: these libraries encapsulate important and well-understood (mathematical) operations, allow for reuse, and are implemented and tuned by experts. Domain scientists then implement complex algorithms based on these domainspecific libraries. While individual library calls are optimized, larger performance gains across sequences of calls—sometimes spanning multiple libraries—are often unrealized, forcing a trade-off between performance and implementation complexity.To overcome this issue, we propose LibraryX, an approach and a system that allows for cross-library-call optimization even when library calls stem from multiple performance libraries. LibraryX annotates library calls with semantic information and optimizes entire directed acyclic graphs (DAGs) of calls dynamically using the SPIRAL code generation system. We demonstrate its effectiveness across a range of memory bound workloads, achieving significant speedups on Nvidia, AMD, and Intel accelerators compared to code using native libraries without cross-call optimization.

Rao, Sanil [Carnegie Mellon University,Department ↗
DOE OSTI2023

Motion Planning Algorithms for Safety and Quantum Computing Efficiency

Motion planning remains a fundamental problem in robotics. Sampling-based algorithms use randomization to allow efficient solutions to this complex problem. As mobile robots and autonomous vehicles become more prevalent in everyday life, motion planning must be applied to increasingly challenging scenarios. Safety has become a paramount concern in motion planning for ensuring robotic applications enrich human lives. To date, many motion planning techniques to increase safety in the face of uncertain and dynamic environments have been developed. This dissertation first addresses distributional safety of Rapidly-Exploring Random Trees (RRT) through our algorithm W-Safe RRT. To acknowledge distributional uncertainty and poor modeling, W-Safe RRT uses the Wasserstein metric to provide a probabilistic bound on the distributional distance between a robot and obstacles. Human-interpretable environmental agent classification allows online safety margin adaptation. We propose and analyze an integrating region method for online classification that increases actor labeling accuracy based on behavioral feature values when compared to state of the art methods. The method performs class assignments based on local maximum likelihood in a created behavioral feature-space, allowing a notion of classification uncertainty. Model-based methods with safety guarantees can quickly become computationally in tractable, especially with multiple agents, higher dimensions, and plentiful unknowns. Sampling based algorithms have been parallelized for computation with multi-core computers and GPUs. We consider the use of quantum algorithms and computers for sampling-based motion planning for the first time. Quantum computing performs operations on superpositions of states and can solve certain problems much more efficiently than classical computers, but introduces previously unseen challenges. With Quantum-RRT, we recast the motion planning problem into a database-search structure and use Quantum Amplitude Amplification to find reachable states in the database with a quadratic performance increase over classical methods. We address two error sources with this method: quantum measurement and quantum oracle errors. We then extend this method to Parallel Quantum-RRT, which uses a manager-worker architecture with multiple parallel quantum workers to increase database search efficiency. We compare algorithm architectures and characterize probabilities of multiple workers finding solutions. Lastly, we test in simulation the quantum algorithms against classical versions in a wide variety of scenarios, concluding that a similar parallelization improvement is to be found in the quantum case as was found in the parallelization of classical RRT.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Soil organic carbon is not just for soil scientists: measurement recommendations for diverse practitioners

Soil organic carbon (SOC) regulates terrestrial ecosystem functioning, provides diverse energy sources for soil microorganisms, governs soil structure, and regulates the availability of organically bound nutrients. Investigators in increasingly diverse disciplines recognize how quantifying SOC attributes can provide insight about ecological states and processes. Today, multiple research networks collect and provide SOC data, and robust, new technologies are available for managing, sharing, and analyzing large data sets. Here, we advocate that the scientific community capitalize on these developments to augment SOC data sets via standardized protocols. We describe why such efforts are important and the breadth of disciplines for which it will be helpful, and outline a tiered approach for standardized sampling of SOC and ancillary variables that ranges from simple to more complex. We target scientists ranging from those with little to no background in soil science to those with more soil-related expertise, and offer examples of the ways in which the resulting data can be organized, shared, and discoverable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Calcitonin Receptor N-Glycosylation Enhances Peptide Hormone Affinity by Controlling Receptor Dynamics

The class B G protein-coupled receptor (GPCR) calcitonin receptor (CTR) is a drug target for osteoporosis and diabetes. N-glycosylation of asparagine 130 in its extracellular domain (ECD) enhances calcitonin hormone affinity with the proximal GlcNAc residue mediating this effect through an unknown mechanism. Here, we present two crystal structures of salmon calcitonin-bound, GlcNAc-bearing CTR ECD at 1.78 and 2.85 Å resolutions and analyze the mechanism of the glycan effect. The N130 GlcNAc does not contact the hormone. Surprisingly, the structures are nearly identical to a structure of hormone-bound, N-glycan-free ECD, which suggested that the GlcNAc might affect CTR dynamics not observed in the static crystallographic snapshots. Hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations revealed that glycosylation stabilized a β-sheet adjacent to the N130 GlcNAc and the N-terminal α-helix near the peptide-binding site while increasing flexibility of the peptide-binding site turret loop. These changes due to N-glycosylation increased the ligand on-rate and decreased its off-rate. The glycan effect extended to RAMP-CTR amylin receptor complexes and was also conserved in the related CGRP receptor. These results reveal that N-glycosylation can modulate GPCR function by altering receptor dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

QM Investigation of Rare Earth Ion Interactions with First Hydration Shell Waters and Protein-Based Coordination Models

Here, conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln 3+ ) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln 3+ interaction with representatives of solvent (water, H 2 O) and protein binding sites (acetate, CH 3 COO – ), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln 3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln 3+ interactions with (H 2 O) 9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu 3+ , which bound less favorably than the preceding element, Yb 3+ . Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln 3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln 3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.

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

RAD51C-XRCC3 structure and cancer patient mutations define DNA replication roles

$RAD51C$ is an enigmatic predisposition gene for breast, ovarian, and prostate cancer. Currently, missing structural and related functional understanding limits patient mutation interpretation to homology-directed repair (HDR) function analysis. Here we report the RAD51C-XRCC3 (CX3) X-ray co-crystal structure with bound ATP analog and define separable RAD51C replication stability roles informed by its three-dimensional structure, assembly, and unappreciated polymerization motif. Mapping of cancer patient mutations as a functional guide confirms ATP-binding matching RAD51 recombinase yet highlights distinct CX3 interfaces. Analyses of CRISPR/Cas9-edited human cells with $RAD51C$ mutations combined with single-molecule, single-cell and biophysics measurements uncover discrete CX3 regions for DNA replication fork protection, restart and reversal, accomplished by separable functions in DNA binding and implied 5’ RAD51 filament capping. Collective findings establish CX3 as a cancer-relevant replication stress response complex, show how HDR-proficient variants could contribute to tumor development, and identify regions to aid functional testing and classification of cancer mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Self-similar mesocrystals form via interface-driven nucleation and assembly

Crystallization by particle attachment (CPA) is a frequently occurring mechanism of colloidal crystallization that results in hierarchical morphologies. CPA has been exploited to create nanomaterials with unusual properties and is implicated in the development of complex mineral textures. Oriented attachment–a form of CPA in which particles align along specific crystallographic directions–produces mesocrystals that diffract as single crystals do, although the constituent particles are still discernible. The conventional view of CPA is that nucleation provides a supply of particles that aggregate via Brownian motion biased by attractive interparticle potentials. However, mesocrystals often exhibit regular morphologies and uniform sizes. Although many crystal systems form mesocrystals and individual attachment events have been directly visualized10, how random attachment events lead to well defined, self-similar morphologies remains unknown, as does the role of surface-bound ligands, which are ubiquitous in nanoparticle systems. Attempts to understand mesocrystal formation are further complicated in many systems by the presence of precursor nanoparticles with a phase distinct from that of the bulk. Some studies propose that such particles convert before attachment, whereas others attribute conversion to the attachment process itself and yet others conclude that transformation occurs after the mesocrystals exceed a characteristic size. Here we investigate mesocrystal formation by iron oxides, which are important colloidal phases in natural environments and classic examples of systems forming ubiquitous precursor phases and undergoing CPA accompanied by phase transformations. Combining in situ transmission electron microscopy (TEM) at 80 degrees Celsius with ‘freeze-and-look’ TEM, we tracked the formation of haematite (Hm) mesocrystals in the presence of oxalate (Ox), which is abundant in soils, where iron oxides are common. In this work, we find that isolated Hm particles rarely appear, but once formed, interfacial gradients at the Ox-covered surfaces drive Hm particles to nucleate repeatedly about two nanometres from the surfaces, to which they then attach, thereby generating mesocrystals. Comparison to natural and synthetic systems suggests that interface-driven pathways are widespread.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding Strong Scaling on GPUs Using Empirical Performance Saturation Size

The roofline model provides a concise overview of the maximum performance capabilities of a given computer system through a combination of peak memory bandwidth and compute performance rates. The increasing complexity of scheduling and cache in recent GPUs, however, has introduced complicated performance variability that is not captured by arithmetic intensity alone. This work examines the effect of problem size and GPU launch configurations on roofline performance for V100, A100, MI100, and MI250X graphics processing units. We introduce an extended roofline model that takes problem size into account, and find that strong scaling on GPUs can be characterized by saturation problem sizes as additional key metrics. Saturation problem sizes break up a plot of GPU performance vs. problem size into three distinct performance regimes– size-limited, cache-bound, and DRAM-bound. With our extended roofline model, we are able to provide a robust view of these performance regimes across recent GPU architectures.

Eberius, David↗
DOE OSTI2020

Exploring the unknown Lambda-neutron interaction

No published \Lambda Λ n scattering data exist. A relativistic heavy-ion experiment has suggested that a \Lambda Λ nn bound state was seen. However, several theoretical analyses have cast serious doubt on the bound-state assertion. Nevertheless, there could exist a three-body \Lambda Λ nn resonance. Such a resonance could be used to constrain the \Lambda Λ n interaction. We discuss \Lambda Λ nn calculations using nn and \Lambda Λ n pairwise interactions of rank-one, separable form that fit effective range parameters of the nn system and those hypothesized for the as yet unobserved \Lambda Λ n system based upon four different \Lambda Λ N potentials. The use of rank-one separable potentials allows one to analytically continue the \Lambda Λ nn Faddeev equations onto the second complex energy plane in search of resonance poles, by examining the eigenvalue spectrum of the kernel of the Faddeev equations. Although each of the potential models predicts a \Lambda Λ nn sub-threshold resonance pole, scaling of the \Lambda Λ n interaction by as little as \sim ∼ 5% does produce a physical resonance. This suggests that one may use photo-(electro-)production of the \Lambda Λ nn system from tritium as a tool to examine the strength of the \Lambda Λ n interaction.

Gibson, Benjamin↗
DOE OSTI2021

Synthetic Evaluation of Infrasonic Multipole Waveform Inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies are able to help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

58 GEOSCIENCES↗
DOE OSTI2022

Deep Learning Based Frequency Stability Assessment in Power Grid with High Renewables

Frequency stability assessment is one critical aspect of power system security assessment. Traditional N-1 screening method is based on the simulations of a few typical daily and seasonal operation scenarios. However, the increasing integration of inverter-based renewables and the retirement of conventional synchronous generators result in decreasing system inertia and growing complexity of system operating conditions. Selecting a few typical operation scenarios cannot cover all operating conditions, and the time-domain simulation of all operation conditions requires tremendous time. This paper proposes a more efficient frequency stability assessment method based on deep learning. The affinity propagation clustering algorithm is used to divide the dataset into different clusters, so the selected dataset for training can cover the diversified operating conditions as much as possible. Also, feature normalization is applied to both the training dataset and testing dataset in order to remove any unnecessary bias. Especially, trained model based on full dataset normalization has bounded error in the prediction. The case study on the reduced 240-bus WECC system demonstrates that the proposed method can predict accurate frequency nadir with limited training dataset. The deep learning model using the revised feature normalization can predict more accurate frequency nadir than that using the traditional feature normalization and has very small maximum prediction error.

affinity propagation↗
DOE OSTI2023

Recognition of centromere‐specific histone Cse4 by the inner kinetochore Okp1‐Ame1 complex

Abstract Successful mitosis depends on the timely establishment of correct chromosomal attachments to microtubules. The kinetochore, a modular multiprotein complex, mediates this connection by recognizing specialized chromatin containing a histone H3 variant called Cse4 in budding yeast and CENP‐A in vertebrates. Structural features of the kinetochore that enable discrimination between Cse4/CENP‐A and H3 have been identified in several species. How and when these contribute to centromere recognition and how they relate to the overall structure of the inner kinetochore are unsettled questions. More generally, this molecular recognition ensures that only one kinetochore is built on each chromatid and that this happens at the right place on the chromatin fiber. We have determined the crystal structure of a Cse4 peptide bound to the essential inner kinetochore Okp1‐Ame1 heterodimer from budding yeast. The structure and related experiments show in detail an essential point of Cse4 contact and provide information about the arrangement of the inner kinetochore.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Selective CO 2 Reduction by Bis(bipyridine)cobalt(II) Catalysts: The Role of Pendant Pyridine as a Proton Acceptor

Electrochemical CO 2 reduction reaction (CO 2 RR) catalyzed by molecular earth-abundant metal catalysts is a promising strategy to convert CO 2 into value-added products. One recent trend in this field has been focusing on the rational design of catalysts by incorporating redox-active ligands and modifying the secondary coordination sphere (SCS) to achieve efficient and selective CO 2 RR. Herein, we report a series of Co bis­(bipyridine) catalysts featuring various dangling groups, such as pyridine, tertiary amine, or butyl, in the secondary coordination sphere (Co-PyMe, Co-Py, Co-PrN, and Co-Bu). Efficient, selective electrocatalytic CO 2 RR was achieved by the complexes after the generation of triply reduced intermediate consisting of a Co I center and a dianionic ligand, producing CO as the major product and trace amount of H 2 . Strong correlations with the identity of dangling groups and turnover frequency (TOF) have been observed, in which Co-PyMe displayed the highest TOF (1086 s –1 in MeCN/H 2 O). Mechanistic studies indicated that the acceleration of CO 2 RR with pyridine-functionalized catalysts were derived from the protonation of pyridine dangling groups which participated as weak acids in the H-bonding network with exogenous proton sources, stabilizing CO 2 -bound intermediates and facilitating proton transfer. In addition, precatalytic CO 2 binding and activation at the third reduction (−2.1 V) was revealed by CV and SEC-IR studies. The resultant doubly reduced CO-bound species acted as a trapping state which inhibited CO 2 RR electrocatalysis. Regeneration of active species was accessed via reductive dissociation of CO at a more negative potential. In conclusion, this study highlights the combined effects of redox-active ligands and pyridine/pyridinium as SCS groups on CO 2 RR catalysis and provides design principles for future development of CO 2 RR catalysts utilizing pyridine/pyridiniums as SCS functional groups to fine-tune the catalytic activity.

CO2 reduction↗
DOE OSTI2022

Spiers Memorial Lecture: Theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our ability to model and predict the kinetics of such reactions has grown in leaps and bounds through the synergy of ever more sophisticated experimental studies with increasingly high level theoretical analyses. Herein, this review begins with a brief historical overview of the progress from the Lindemann mechanism to the master equation, which now provides the cornerstone for most theoretical analyses. The current status of ab initio transition state theory based implementations of the master equation is then reviewed, beginning with a discussion of the energy resolved chemical conversion rates, followed by a review of the pressure dependence of the thermal kinetics. The latter discussion focusses on the collisional energy and angular momentum transfer rates as well as the complexities of non-thermal effects and of multiple-well multiple-channel reactions. The synergy with recent state-of-the-art experiments is used to motivate these discussions. Attempts to automate the calculations are also briefly reviewed. Throughout the discussion, we provide our perspective on current understanding and continuing challenges and opportunities. One key challenge relates to the coupling of sequences of reactions, which frequently leads to deviations from the classic presumption of thermal reactants.

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

Gradient flow based phase-field modeling using separable neural networks

Allen–Cahn equation is a reaction–diffusion equation and is widely used for modeling phase separation. Machine learning methods for solving the Allen–Cahn equation in its strong form suffer from inaccuracies in collocation techniques, errors in computing higher-order spatial derivatives, and the large system size required by the space–time approach. To overcome these challenges, we propose solving the gradient flow of the Ginzburg–Landau free energy functional, which is equivalent to the Allen–Cahn equation, thereby avoiding the second-order spatial derivatives associated with the Allen–Cahn equation. A minimizing movement scheme is employed to solve the gradient flow problem, eliminating the complexities of a space–time approach. We utilize a separable neural network that efficiently represents the phase field through low-rank tensor decomposition. As we use the minimizing movement scheme to numerically solve the gradient flow problem, we thus, refer to the proposed method as the Separable Deep Minimizing Movement (SDMM) method. The evaluation of the functional in the minimizing movement scheme using the Gauss quadrature technique bypasses the inaccuracies associated with collocation techniques traditionally used to solve partial differential equations. A hyperbolic tangent transformation is introduced on the phase field prior to the evaluation of the functional to ensure that it remains strictly bounded within the values of the two phases. For this transformation, theoretical guarantee for energy stability of the minimizing movement scheme is established. Our results suggest that this transformation helps to improve the accuracy and efficiency significantly. The proposed method resolves the challenges faced by state-of-the-art machine learning techniques, outperforming them in both accuracy and efficiency. It is also the first machine learning method to achieve an order of magnitude speed improvement over the finite element method. In addition to its formulation and computational implementation, several case studies illustrate the applicability of the proposed method.

42 ENGINEERING↗
DOE OSTI2023

Allosterically coupled conformational dynamics in solution prepare the sterol transfer protein StarD4 to release its cargo upon interaction with target membranes

Complex mechanisms regulate the cellular distribution of cholesterol, a critical component of eukaryote membranes involved in regulation of membrane protein functions directly and through the physiochemical properties of membranes. StarD4, a member of the steroidogenic acute regulator-related lipid-transfer (StART) domain (StARD)-containing protein family, is a highly efficient sterol-specific transfer protein involved in cholesterol homeostasis. Its mechanism of cargo loading and release remains unknown despite recent insights into the key role of phosphatidylinositol phosphates in modulating its interactions with target membranes. We have used large-scale atomistic Molecular dynamics (MD) simulations to study how the dynamics of cholesterol bound to the StarD4 protein can affect interaction with target membranes, and cargo delivery. We identify the two major cholesterol (CHL) binding modes in the hydrophobic pocket of StarD4, one near S136&S147 (the Ser-mode), and another closer to the putative release gate located near W171, R92&Y117 (the Trp-mode). We show that conformational changes of StarD4 associated directly with the transition between these binding modes facilitate the opening of the gate. To understand the dynamics of this connection we apply a machine-learning algorithm for the detection of rare events in MD trajectories (RED), which reveals the structural motifs involved in the opening of a front gate and a back corridor in the StarD4 structure occurring together with the spontaneous transition of CHL from the Ser-mode of binding to the Trp-mode. Further analysis of MD trajectory data with the information-theory based NbIT method reveals the allosteric network connecting the CHL binding site to the functionally important structural components of the gate and corridor. Mutations of residues in the allosteric network are shown to affect the performance of the allosteric connection. These findings outline an allosteric mechanism which prepares the CHL-bound StarD4 to release and deliver the cargo when it is bound to the target membrane.

59 BASIC BIOLOGICAL SCIENCES↗