Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate computing”

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 253 records · Page 14

Machine learning magnetism classifiers from atomic coordinates

The determination of magnetic structure poses a long-standing challenge in condensed matter physics and materials science. Experimental techniques such as neutron diffraction are resource-limited and require complex structure refinement protocols, while computational approaches such as first-principles density functional theory (DFT) need additional semi-empirical correction, and reliable prediction is still largely limited to collinear magnetism. Here, we present a machine learning model that aims to classify the magnetic structure by inputting atomic coordinates containing transition metal and rare earth elements. By building a Euclidean equivariant neural network that preserves the crystallographic symmetry, the magnetic structure (ferromagnetic, antiferromagnetic, and nonmagnetic) and magnetic propagation vector (zero or non-zero) can be predicted with an average accuracy of 77.8% and 73.6%. In particular, a 91% accuracy is reached when predicting no magnetic ordering even if the structure contains magneticelement(s). Ourworkrepresents onestepforwardtosolvingthegrand challenge of full magnetic structure determination.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Videos, photos, and AI-derived grain size data associated with “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization” under review. This data package includes five data types: 1) raw photos and videos from drone survey and walking smartphone surveys; 2) images derived from raw videos; 3) manual labeling of reference scales; 4) metadata for all images and photo resolution derived from artificial intelligence (AI) models or manual labels, 5) grain size data obtained from AI models for all photos, 6) metadata and grain size data after quality control, 7) summaries of sample efficiency for all data, and 8) computational fluid dynamics (CFD) data used to support hydro-biogeochemical (HBGC) parameter estimation. Such data is used to 1) demonstrate significant improvements in accuracy, efficiency, and quality control for grain size data collection with the help of AI models, 2) study the spatial heterogeneity of grain size and observation reproducibility based on tens of thousands of data points generated by the AI models, and 3) evaluate the impacts of grain size heterogeneity on key HBGC parameters across sediment-to-reach and hourly-to-yearly scales. In particular, the data package contains 116 folders and 179696 files. The files include 41 videos in .mov format, 64047 photos in .jpg format, 13541 video-derived photos in .png format, 12747 segmentation mask data in .tif format, 12747 segmentation data in .json format, 24771 .csv files that with metadata and grain size for each individual photo as well as water depth and velocity data from CFD and observation, 51791 .txt files of raw AI predicted labels, and 11 flight record data in .srt format. The summary for all metadata and grain size statistics information is included in “Scales_V3_NG.csv” and “Statistics_V3_NG.csv”. The summary for data that pass data quality control (QC) level 0-2 is included in “QCStatistics_V3_NG.csv”. The QC level 0 represents photos whose photo resolution is positive, excluding photos that miss reference scale. The QC level 1 means reference scale circularity uncertainty is less than 5% for smartphone images while representing photo resolution is larger than 0.44 mm/pixel for drone images. The QC level 2 means excluding photos whose grain number is less than 100, a minimum number of grains recommended by classic literature. The summary for each video’s name, length, frame rates, survey area, grain number, survey efficiency, etc. can be found in “QCSummary_V3_NG.csv”. The summary for site name, GPS coordinates, and number of images at each site can be found in “SitesSummary_V3_*.csv” files. Overall computational efficiency summary is reported in Table 4 of accompanying manuscript. Additionally, the nitrate concentration data used in this work was downloaded from an existing dataset published on ESS-DIVE (Boat-Dragged Sensor Hanford Reach.csv; Conner A. et al., 2020). We thank the United States Forest Service, Washington Department of Fish and Wildlife, Washington Department of Natural Resources, Cowiche Canyon Conservatory, Port of Benton, and the Confederated Tribes and Bands of the Yakama Nation for access to field locations where the data were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate data collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

Coordinate transformation and construction of finite element mesh in a diverted tokamak geometry

A coordinate transformation technique between straight magnetic field line coordinate system (Ψ, θ) and Cartesian coordinate system (R, Z) is presented employing a Solov'ev solution of the Grad-Shafranov equation. Employing the equilibrium solution, the poloidal magnetic flux Ψ(R, Z) of a diverted tokamak, magnetic field line equation is solved computationally to find curves of constant poloidal angle θ, which provides us with explicit relations R = R(Ψ, θ) and Z = Z(Ψ, θ). Correspondingly, conversion from one coordinate to the other along particle trajectories in the vicinity of separatrix is demonstrated. Based on the magnetic structure, a finite element mesh is generated in a diverted tokamak geometry to solve Poisson's equation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

S4PST: Sustainability for Programming Systems and Tools: May Workshop Report

The US Department of Energy (DOE) Exascale Computing Project (ECP) has fostered and strengthened the use of modern software engineering practices for developing applications and libraries, and this effort has resulted in the coordinated and interoperable E4S1 and xSDK2 ecosystems. Although this approach is cost-effective, it relies on robust programming systems and tools (PST) as the underlying foundation for our HPC software. At present, our primary PST stack consists of traditional high-performance computing (HPC) languages, namely Fortran, C, C++, and the popular Python language for data analysis and AI workflows. These languages support various programming frameworks and run-time abstractions that enable parallelism and concurrency across multiple node architectures and thousands of nodes through a variety of interconnect systems. However, to accommodate users’ diverse needs, certain aspects of the HPC ecosystem are delegated to vendor-specific or third-party implementations that extend beyond a particular scientific domain. This broader scope results in a multitude of specifications and variations, which leads to a complex orchestration of many-ecosystems. Unfortunately, this complexity in the ecosystem imposes additional overhead costs on consumers during the latter stages of the development cycle. In addition to the software ecosystem challenge, the upcoming conclusion of the ECP by December 2023 has raised significant concerns within the HPC programming systems community, from both the economic and social perspectives. The ECP has implemented a management structure for software development and funding decisions across all ECP participants by following a conventional hierarchical and centralized approach. However, this structure has prompted certain considerations within the community, particularly in anticipation of the Software Sustainability initiative by the DOE’s Advanced Scientific Computing Research Program (ASCR). For the success of this new initiative, it is of utmost importance to secure consistent funding and foster close engagement with researchers and core developers of existing programming-system products. This collaboration is vital to maintaining the critical capabilities of the current software during the transition phase while proactively adapting to future technology and workforce trends. The community recognizes the significance of adapting to emerging trends and is aware of the inherent fragility of the HPC software ecosystem, particularly in relation to programming systems that cater to all users. The ability to adapt and evolve is essential to staying relevant and effectively addressing these technical, economic, and social challenges. The S4PST team, which represents one of the six ASCR Software Sustainability seedling projects, is dedicated to tackling these challenges through community-based approaches that go beyond the scope of the DOE. This involves collaboration between national laboratories with academia, non-DOE institutions, hardware and system vendors, and international partners. By fostering these partnerships, we aim to create a robust and sustainable HPC software ecosystem that can effectively meet the needs of the community. This new community effort, driven by the eight DOE labs, will take on the responsibility of guiding funding decisions for programming-systems development and maintenance with transparency and consistency across all decisions. Additionally, the team will offer common technical services to the programming systems community, irrespective of their funding situations, and facilitate community-wide incubation to proactively nurture the software ecosystem. By actively engaging with stakeholders and employing a collaborative approach, we can collectively shape the future of programming systems and ensure a robust and thriving HPC software landscape. On May 11–12, 2023, the S4PST team conducted its inaugural kick-off workshop at the Innovative Computing Laboratory (ICL) in the University of Tennessee, Knoxville, hosted by Hartwig Anzt. The workshop encompassed various sessions dedicated to presentations and discussions, with the aim of comprehending the team members’ perspectives on the vision of software sustainability. Additionally, the workshop aimed to identify the technical, economic, and social requirements for sustaining the programming-systems community in the field of HPC. This report provides a summary of the S4PST effort by highlighting five major thrust areas discussed during the workshop: (i) community, (ii) technical support, (iii) training and diversity, (iv) verification, validation and correctness, and (v) emerging technologies. It also encompasses an overview of the presentations and discussions held throughout the event, our views and potential synergies with other seedling efforts, along with the outcomes and key takeaways from our initial discussions.

97 MATHEMATICS AND COMPUTING↗

Solution and Active Site Speciation Drive Selectivity for Electrocatalytic Reactive Carbon Capture in Diethanolamine over Ni–N–C Catalysts

Direct conversion of captured forms of carbon, or reactive carbon capture (RCC), presents an opportunity to reduce the energy intensity and cost of direct CO 2 utilization from dilute sources. While amine-based sorbents effectively capture CO 2 , their use for RCC presents numerous challenges with typical pure metal catalysts used for electrochemical CO 2 reduction (CO 2 R). Here, using both theory and experiments, we find that Ni–N–C single atom catalysts are effective for RCC conversion to CO using a diethanolamine sorbent, in contrast to pure metal catalysts. Computational analysis reveals that RCC can proceed directly through direct reduction of the sorbent-CO 2 adduct or indirectly by C–N bond breaking facilitating CO 2 adsorption and subsequent reduction. We find that the latter mechanism is most prevalent at low overpotentials where we experimentally observe RCC selectivity. We also find experimentally that the rate of CO production for RCC with Ni–N–C catalysts can exceed pure bicarbonate solutions at intermediate sorbent concentration (0.1–0.5 M DEA) under dilute (10–25%) streams of CO 2 at low overpotentials. The coordination environment of Ni sites and the solution speciation influence their RCC activity, with changes in protonation to coordinating N/C atoms resulting in changing the RCC mechanism and consequent activity. In situ X-ray absorption spectroscopy and computational analysis reveal restructuring under RCC conditions due to hydrogen coadsorption with DEA that limits the stability of Ni–N–C catalysts. This work highlights the importance of carefully controlling the catalyst and solution environment to achieve active and stable RCC electrocatalysis.

Chemistry↗

Four-Coordinate Fe N 2 and Imido Complexes Supported by a Hemilabile NNC Heteroscorpionate Ligand

Inspired by mechanistic proposals for N 2 reduction at the nitrogenase FeMo cofactor, we report herein a new, strongly σ-donating heteroscorpionate ligand featuring two weak-field pyrazoles and an alkyl donor. This ligand supports four-coordinate Fe(I)-N 2 , Fe(II)-Cl, and Fe(III)-imido complexes, which we have characterized using a variety of spectroscopic and computational methods. Structural and quantum mechanical analysis reveal the nature of the Fe–C bonds to be essentially invariant between the complexes, with conversion between the (formally) low-valent Fe-N 2 and high-valent Fe-imido complexes mediated by pyrazole hemilability. This presents a useful strategy for substrate reduction at such low-coordinate centers and suggests a mechanism by which FeMoco might accommodate the binding of both π-acidic and π-basic nitrogenous substrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Analysis of Bis(trifluoromethylsulfonyl)imide Interactions with Metal Cations Through a Chemical Informatics Approach

Nominally weakly coordinating anions are useful for modulating the solubility and chemical properties of metal complexes, but identification and analysis of the systematics of the interactions of anions with cationic metal complexes has not received the attention it deserves. Here, a chemical informatics approach is demonstrated for identifying and quantitatively analyzing the ways that the bis(trifluoromethylsulfonyl)imide anion (TFSI) can interact with metal-containing species. An open access computer program (PyCIFTer) was developed to facilitate large-scale structural analysis of TFSI-containing species by utilization of experimental atomic coordinate data from single-crystal X-ray diffraction (XRD) studies obtained from the Cambridge Structural Database (CSD). PyCIFTer establishes a three-dimensional vector space from the raw atomic coordinates, generating acyclic, undirected graphs that are used to rapidly analyze the structural properties (bond lengths and angles) of TFSI in individual structures in sequential/batch fashion. The structures are sorted by PyCIFTer into groups based on pre-set and chemically sensible criteria, affording a comprehensive and systematic view of TFSI structural chemistry. This approach avoids tedious one-at-a-time interrogation of structures, a prospect unreasonable in this case, and many others of contemporary chemical relevance; there were over 1500 structures in the CSD containing TFSI as of November 2024. The results demonstrate that TFSI only rarely binds to cations in the solid state, favoring the formation of species in which TFSI is found in cations’ outer coordination spheres. The prospect of applying PyCIFTer to other moieties is also discussed. PyCIFTer is also schematically compared to the commercial CSD Python application programming interface (API). Taken together, this work demonstrates the usefulness of modular workflows for sequential/batch analysis of structural data from XRD, an approach that appears poised to accelerate the translation of legacy structural results into new chemical insights and hypotheses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strong Electronic and Magnetic Coupling in M 4 (M = Ni, Cu) Clusters via Direct Orbital Interactions between Low-Coordinate Metal Centers

In this report we present an extensive study of tetranuclear transition-metal cluster compounds M 4 (NP t Bu 3 ) 4 and [M 4 (NP t Bu 3 ) 4 ][B(C 6 F 5 ) 4 ] (M = Ni, Cu; t Bu = tert -butyl), which feature low-coordinate metal centers and direct metal-metal orbital overlap. X-ray diffraction, electrochemical, magnetic, spectroscopic, and computational analysis elucidate the nature of the bonding interactions in these clusters and the impact of these interactions on the electronic and magnetic properties. Direct orbital overlap results in strongly coupled, large-spin ground states in the [Ni 4 (NP t Bu 3 ) 4 ] +/0 clusters and fully delocalized, spin-correlated electrons. Correlated electronic structure calculations confirm the presence of ferromagnetic ground states that arise from direct exchange between magnetic orbitals, and, in the case of the neutral cluster, itinerant electron magnetism similar to that in metallic ferromagnets. The cationic nickel cluster also possesses large magnetic anisotropy exemplified by a large, positive axial zero-field splitting parameter of D = +7.95 or +9.2 cm -1 , as determined by magnetometry or electron paramagnetic resonance spectroscopy, respectively. The [Ni 4 (NP t Bu 3 ) 4 ] + cluster is also the first molecule with easy-plane magnetic anisotropy to exhibit zero-field slow magnetic relaxation, and under a small applied field, it exhibits relaxation exclusively through an Orbach mechanism with a spin relaxation barrier of 16 cm -1 . The S = 1 / 2 complex [Cu 4 (NP t Bu 3 ) 4 ] + exhibits slow magnetic relaxation via a Raman process on the millisecond time scale, supporting the presence of slow relaxation via an Orbach process in the nickel analogue. Overall, this work highlights the unique electronic and magnetic properties that can be realized in metal clusters featuring direct metal-metal orbital interactions between low-coordinate metal centers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comprehensive compilations of computation results and validations for neutronics start-up tests at China Experimental Fast Reactor

This paper compiles and analyzes refined results of the coordinated research project (CRP) on Neutronics Benchmark of China Experimental Fast Reactor (CEFR) Start-Up Tests conducted by the International Atomic Energy Agency (IAEA) since 2018. Twenty-eight research organizations participated with various code systems. The China Institute of Atomic Energy (CIAE) provided the benchmark specifications and participants conducted the benchmark through blind and refined phases. This paper presents the benchmark results on six experimental measurements, criticality, control rod worth, temperature reactivity coefficients, sodium void worth, assembly swap reactivities, and foil activations. Except for a few outlier results, the simulation results show good agreement with the measured data within 1-σ experiment uncertainties. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Modeling separation of lanthanides via heterogeneous ligand binding

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln 3+ ) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln 3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln 3+ coordinated to ligands in water and embedded in metal–organic frameworks (MOFs), and ΔΔE for Ln 3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln 3+ in the periodic table. Static DFT calculations with a full Ln 3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. Finally, we also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗

Pressure-Induced Metal-like Transport and Magnetoresistance in a Au 2+ –Au 3+ Halide Perovskite

The Cs 4 Au II Au III 2 Cl 12 perovskite (1), featuring AuCl 4 trimers separated by vacancies, enables the first high-pressure study of Au 2+/3+ mixed-valence. Our computational analysis of the gold frontier orbitals suggests that the Au 2+ →Au 3+ intervalence charge transfer (IVCT) occurs across the vacancies. Computational structures indicate that these vacancies rapidly shrink with pressure and the Au 2+ and Au 3+ coordination spheres become very similar at the phase transition to nearly cubic symmetry at ca. 15 GPa─enabling facile IVCT. Although the activation energy of conductivity of 0.73(4) meV and far-infrared absorption indicate a small but nonzero bandgap, ambient thermal energy drives the IVCT, affording metallic properties: prominent infrared reflectivity and transport values of 10 2 S·cm –1 . This prompted us to perform the first high-pressure studies of magnetoresistance (MR) and Hall effect in halide perovskites. At 16 GPa, the MR increases by 9.3% at 2 K and 9 T; this value is maintained up to 27 GPa, when a local distortion drives electronic localization. By globally fitting the MR and Hall resistance to a two-carrier model we quantify how the carrier densities and mobilities evolve with pressure. Thus, metal-like transport and MR in 1 is driven by a pressure-induced transition from localized to partially delocalized mixed-valence.

Deschene, Christina R. [Stanford University, CA (U↗

Deciphering the Olefin Isomerization-Polymerization Paradox of Palladium(II) Diimine Catalysts: Discovery of Simultaneous and Independent Pathways of Olefin Isomerization and Living Polymerization

This work elucidates a long-standing unexplained paradox commonly observed within the polymerization of α-olefin using palladium (Pd)(II)–diimine catalysts, in which isomerization and living polymerization of α-olefins are both observed. With a classical mechanistic understanding of these complexes, this behavior is often dismissed and interpreted as experimental error. Herein, we present a comprehensive mechanistic investigation into this phenomenon that supports the existence of a novel mechanistic pathway for Pd(II)–diimine complexes. Part one of the mechanistic study lays the foundation of the proposed mechanism, in which neutral Pd(II)–diimine complexes were found to exhibit a moderate to good catalytic activity for olefin isomerization of α-olefins despite the established notion that catalyst activation is required. Extensive experimental and computational studies reveal the possibility of a partial dissociation of the diimine ligand, which frees up one coordination site and enables coordination–insertion. This finding is significant as the coexistence of two reactive coordination sites at the palladium center becomes a valid proposal for the activated cationic Pd(II)–diimine complexes. In part two, we examined and validated the simultaneously observed α-olefin isomerization and living polymerization using the cationic Pd(II)–diimine catalyst, which supports the presence of two independent reaction pathways of isomerization and polymerization, respectively. Furthermore, the addition of a strong Lewis acid, such as AlCl 3 , accelerates the ligand dissociation and the consequential isomerization as it weakens the palladium–nitrogen bond through competitive binding. In part three, Lewis acid-triggered olefin isomerization-polymerization is employed to prepare living olefinic block copolymers and further synthesize novel polyolefin-polar block copolymers with unique architectures, distinct levels of branching, crystallinity, and polar functionality in a one-pot manner.

Catalysts↗

Electrocatalytic synthesis of heterocycles from biomass-derived furfuryl alcohols

It is very attractive yet underexplored to synthesize heterocyclic moieties pertaining to biologically active molecules from biomass-based starting compounds. Herein, we report an electrocatalytic Achmatowicz reaction for the synthesis of hydropyranones from furfuryl alcohols, which can be readily produced from biomass-derived and industrially available furfural. Taking advantage of photo-induced polymerization of a bipyridyl ligand, we demonstrate the facile preparation of a heterogenized nickel electrocatalyst, which effectively drives the Achmatowicz reaction electrochemically. A suite of characterization techniques and density functional theory computations were performed to aid the understanding of the reaction mechanism. It is rationalized that the unsaturated coordination sphere of nickel sites in our electrocatalyst plays an important role at low applied potential, not only allowing the intimate interaction between the nickel center and furfuryl alcohol but also enabling the transfer of hydroxide from nickel to the bound furfuryl alcohol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CO electroreduction on single-atom copper

Electroreduction of carbon dioxide (CO 2 ) or carbon monoxide (CO) toward C 2+ hydrocarbons such as ethylene, ethanol, acetate and propanol represents a promising approach toward carbon-negative electrosynthesis of chemicals. Fundamental understanding of the carbon–carbon (C-C) coupling mechanisms in these electrocatalytic processes is the key to the design and development of electrochemical systems at high energy and carbon conversion efficiencies. Here, we report the investigation of CO electreduction on single-atom copper (Cu) electrocatalysts. Atomically dispersed Cu is coordinated on a carbon nitride substrate to form high-density copper–nitrogen moieties. Chemisorption, electrocatalytic, and computational studies are combined to probe the catalytic mechanisms. Unlike the Langmuir-Hinshelwood mechanism known for copper metal surfaces, the confinement of CO adsorption on the single-copper-atom sites enables an Eley-Rideal type of C-C coupling between adsorbed (*CO) and gaseous [CO(g)] carbon moxide molecules. The isolated Cu sites also selectively stabilize the key reaction intermediates determining the bifurcation of reaction pathways toward different C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

U.S. Department of Energy: National Virtual Biotechnology Laboratory (Technical Report)

With funding from the CARES Act, the U.S Department of Energy (DOE) established the National Virtual Biotechnology Laboratory (NVBL) in March 2020 to address key challenges associated with the COVID-19 crisis. The NVBL brought together the broad scientific and technical expertise and resources of DOE’s 17 national laboratories to help tackle medical supply shortages, discover potential drugs to fight the virus, develop and validate COVID-19 testing methods, model disease spread and impact across the nation, and understand virus transport in buildings and the environment. National laboratory resources leveraged for this effort include a suite of world-leading user facilities broadly available to the research community, such as light and neutron sources, nanoscale science research centers, sequencing and biocharacterization facilities, and high-performance computing facilities. As part of the NVBL framework, DOE rapidly assembled five project teams to (1) identify new targets for medical therapeutics; (2) develop innovations in testing capabilities; (3) provide epidemiological and logistical support; (4) understand viral fate and transport in the environment; and (5) address supply chain bottlenecks by harnessing extensive additive manufacturing capabilities. Each research team was charged with defining high-impact projects that could be completed in a 6-month sprint while coordinating their developments with academia, other government agencies, and the private sector. Within months, NVBL teams used DOE’s high-performance computers and light and neutron sources to identify promising candidates for antibodies and antivirals that universities and drug companies are now evaluating. NVBL researchers also developed new diagnostic targets and sample collection approaches, and supported efforts by the U.S. Food and Drug Administration, Centers for Disease Control and Prevention, and U.S. Department of Defense to establish national guidelines used in administering millions of tests. Researchers used artificial intelligence and high-performance computing to produce near-real-time data analysis to forecast disease transmission, stress on public health infrastructure, and economic impact, which supported decision-makers at the local, state, and national levels. To minimize virus uptake and protect human health, NVBL teams studied how to control indoor virus movement. Researchers also produced innovations in materials and advanced manufacturing that mitigated shortages in test kits and personal protective equipment, creating nearly 1,000 new jobs. Through its NVBL framework, DOE has contributed significantly to the nation’s COVID response, demonstrating in only a few months the critical impact of its national laboratories. NVBL’s accomplishments demonstrate not only the powerful resource represented by DOE’s national laboratories working together to meet national needs, but also the effectiveness of the integrated NVBL framework for rapidly responding to emergencies with research and development solutions. Going forward, the NVBL is poised to apply the unique capabilities and expertise of the national laboratory complex to future national and international emergencies, both natural and engineered. Through this framework, DOE will continue to be an integral component of agency-wide efforts to prepare for and respond to biorisks and other crises. This technical report describes the goals, progress, and results of NVBL’s five project teams—Molecular Design for COVID-19 Therapeutics, COVID-19 Testing, Epidemiological Modeling, Viral Fate and Transport, and Materials and Manufacturing of Critical Supplies—and lists each team’s publications and research output.

42 ENGINEERING↗

Assessing models of force-dependent unbinding rates via infrequent metadynamics

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent “bonds” being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. Furthermore, this work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗