Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional priors”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Design and Optimization of Processes for Recovering Rare Earth Elements from End-of-Life Hard Disk Drives

In this conference paper, we propose a superstructure-based approach to finding the optimal pathways for recovering rare earth elements in their commercialized rare earth oxide form from end-of-life HDDs. The proposed superstructure was modeled as a MILP optimization problem, selecting the net present value as the objective function. Whenever possible, costing data taken from the literature was used to inform this mode. However, due to the novelty of this research area data were often not available thus requiring the generation of flowsheets that were implemented in Aspen Plus. To establish the base case optimal result, projections for the number of EOL HDDs in the U.S. available for recycling and estimates of the projected rare earth oxide prices over the lifetime of the plant were used to inform the model. The model was then expanded to include the recycling of EOL HDDs generated prior to the beginning of plant production (period ranging from 2006 through 2024).

Laliwala, Chris↗
DOE OSTI2022

Psilocybin: crystal structure solutions enable phase analysis of prior art and recently patented examples

Psilocybin {systematic name: 3-[2-(dimethylamino)ethyl]-1 H -indol-4-yl dihydrogen phosphate} is a zwitterionic tryptamine natural product found in numerous species of fungi known for their psychoactive properties. Following its structural elucidation and chemical synthesis in 1959, purified synthetic psilocybin has been evaluated in clinical trials and has shown promise in the treatment of various mental health disorders. In a recent process-scale crystallization investigation, three crystalline forms of psilocybin were repeatedly observed: Hydrate A, Polymorph A, and Polymorph B. The crystal structure for Hydrate A was solved previously by single-crystal X-ray diffraction. This article presents new crystal structure solutions for the two anhydrates, Polymorphs A and B, based on Rietveld refinement using laboratory and synchrotron X-ray diffraction data, and density functional theory (DFT) calculations. Utilizing the three solved structures, an investigation was conducted via Rietveld method (RM) based quantitative phase analysis (QPA) to estimate the contribution of the three different forms in powder X-ray diffraction (PXRD) patterns provided by different sources of bulk psilocybin produced between 1963 and 2021. Over the last 57 years, each of these samples quantitatively reflect one or more of the hydrate and anhydrate polymorphs. In addition to quantitatively evaluating the composition of each sample, this article evaluates correlations between the crystal forms present, corresponding process methods, sample age, and storage conditions. Furthermore, revision is recommended on characterizations in recently granted patents that include descriptions of crystalline psilocybin inappropriately reported as a single-phase `isostructural variant.' Rietveld refinement demonstrated that the claimed material was composed of approximately 81% Polymorph A and 19% Polymorph B, both of which have been identified in historical samples. In this article, we show conclusively that all published data can be explained in terms of three well-defined forms of psilocybin and that no additional forms are needed to explain the diffraction patterns.

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

Full core LOCA safety analysis for a PWR containing high burnup fuel

For economic reasons, the US nuclear industry is renewing efforts to build a technical basis to extend peak rod average burnup limits above the current regulatory burnup limit of 62 GWd/MTU. The primary driver of these efforts is to economically increase pressurized water reactor (PWR) cycle lengths to 24 months, reduce the number of fresh fuel assemblies, increase time online, thereby reducing number of outages and their associated cost, higher fuel utilization, and possibly reduce core design constraints. In order for US nuclear utilities to leverage these economic efficiencies, the US Nuclear Regulatory Commission (NRC) will likely require nuclear power plants (NPPs) to analyze a number of potential operational occurrences and their potential consequences with each new core design prior to resuming normal operation. Potential operational occurrences can be divided into three primary regimes: (1) normal operation, (2) anticipated operation occurrences (AOOs), and (3) design basis accidents (DBAs). Normal plant operation is an operating regime in which the plant operates within specified operational limits until the end of the cycle, whereas AOOs are events that result in the NPP deviating outside the normal operating regime. A key attribute of an AOO is that the occurrence should be expected. However, by definition, the occurrence of an AOO does not result in significant impact to critical safety functions. The last potential operational occurrence is a DBA. From the fuel performance point-of-view, DBAs can be subdivided into two bounding categories: (1) loss of coolant accidents (LOCAs), and (2) reactivity insertion accidents (RIAs). Unlike AOOs, DBAs may result in fuel rod failure. The NRC imposes fundamental acceptance criteria to minimize radiological consequences to the public and onsite staff. Furthermore, safety criteria are typically linked to the fulfillment of other acceptance criteria related to reactor safety equipment designed to mitigate DBAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Importance of enforcing Hund’s rules in density functional theory calculations of rare earth magnetocrystalline anisotropy

Density functional theory (DFT) and its extensions, such as DFT+U and DFT+dynamical mean-field theory, are invaluable for studying magnetic properties in solids. However, rare-earth (R) materials remain challenging due to self-interaction errors and the lack of proper orbital polarization. We show how the orbital dependence of self-interaction error contradicts Hund’s rules and plagues magnetocrystalline anisotropy (MA) calculations, and how analyzing DFT states that respect Hund’s rules can mitigate this issue. We benchmark MA in RCo 5 , R 2 Fe 14 B, and RFe 12 , extending prior work on RMn 6 Sn 6 , achieving excellent agreement with experiments. Additionally, we illustrate a semi-analytical perturbation approach that treats crystal fields as a perturbation in the large spin-orbit coupling limit. Using Gd-4f crystal-field splitting, this method provides a microscopic understanding of MA and enables rapid screening of high-MA materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Equilibrium-gated pattern formation: How molecular dissociation thermodynamics drive emergent behavior in dissipative polymeric systems

Emergent patterns in biological systems arise through dissipative processes that balance reaction and transport phenomena, producing highly functional properties from self-regulating mechanisms. Synthetic fabrication, by contrast, often relies on user-controlled, multistep methods that lack the self-organizing capabilities of natural systems. Inspired by nature, we sought chemical systems that integrate strongly coupled reaction and transport phenomena, identifying frontal ring-opening metathesis polymerization (FROMP) as a method capable of creating diverse forms and functions through reactive processing. By employing discrete molecular initiators, FROMP allows precise control of key reaction steps—inhibition, initiation, and propagation. Using an integrated computational and experimental framework, we uncover how near-equilibrium inhibition dynamics, coupled with far-from-equilibrium reaction kinetics, drive pattern formation in frontally polymerized synthetic materials. We propose the concept of equilibrium-gated pattern formation, demonstrating how initiator chemistry can be tuned to achieve programmable macroscale properties. Our study reveals a surprising insight: Emergent behavior in FROMP systems arises from the inhibition-dominated regime of resin composition, expanding prior observations that such behavior is confined to a narrow compositional space near the boundary between front quenching and uniform front propagation. We identify a broader compositional window, far from the quenching regime, where emergent behavior reliably manifests. This expanded design space significantly enhances the operational flexibility of reactive systems and their capacity for self-organization. Furthermore, these insights provide a roadmap for designing bioinspired materials with self-organizing capabilities, unlocking possibilities in synthetic manufacturing.

Architected Polymers↗
DOE OSTI2020

Tracer particle transport dynamics in the diffusive sandpile cellular automaton

The confinement properties of the diffusive running sandpile are characterized by tracking the motion of a population of marked grains of sand. It is found that, as the relative strength of the avalanching to the diffusive transport channel is varied, a point is reached at which the particle global confinement time and the probability density functions of the jump-sizes and waiting-times of the tracked grains experience a sudden change, thus revealing a dynamical transition, that is consistent with previous studies (Newman DE et al., Phys Rev Lett 2002;88(20):204304). Across this transition, the sandpile moves from a regime characterized by self-similarity and memory, where avalanches of all possible sizes dominate transport across the system, to another regime where transport is taken over by near system-size, quasi-periodic avalanches. Values for the fractional transport exponents that quantify effective transport across the sandpile prior to the transition are also obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Predeployment progress of the Canister Deposition Field Demonstration

This report updates the high-level test plan for evaluating surface deposition on three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUTECH Horizontal Modular Storage (NUHOMS) Advanced Horizontal Storage Modules (AHSMs) from Orano (formerly Transnuclear Inc.) and provides a summary of the surface sampling activities that have been conducted to date. The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or requirements. One goal of this testing is to collect defensible and detailed dust deposition measurements from the surface of dry storage canisters in a marine coastal environment to guide chloride-induced stress corrosion cracking (CISCC) research. Another goal is to provide data for the validation of computational fluid dynamics (CFD) based deposition modeling. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the decay heat and thermal hydraulic environment. Test and heater design is supported by detailed CFD modeling. Instrumentation throughout the canister, storage module, and environment will provide extensive information about the thermal-hydraulic behavior of horizontal dry cask storage systems. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer detailed quantification and composition of the deposited particulates from a realistic storage environment. Discussions of a potential host site for the Canister Deposition Field Demonstration (CDFD) are ongoing. Until a host site is chosen, testing of key CDFD hardware components including the heater assemblies, power skid, and remote data acquisition system will continue. Functional testing of the finalized heater assemblies and test apparatus started this fiscal year. These initial heater tests have shown the assemblies are performing within design specifications. Staged surface sampling of a mockup of a canister outside the AHSM on a transfer skid was also performed. Refinements to the sampling procedures and techniques were captured from observation of these activities and lessons-learned debriefs. These updated sampling procedures and techniques are planned to be tested again in the field using the mockup in order to assure personnel are using the most accurate and repeatable methods possible prior to deployment for actual CDFD testing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2022

In Situ 29 Si solid-state NMR study of grafting of organoalkoxysilanes to mesoporous silica nanoparticles

We report that since the catalytic activity and the stability of silica-bound organometallic complexes are affected by their interactions with hydroxyl groups on the surface, isolated hydroxyls are often created prior to the introduction of catalytic species. Here, we investigate a method to remove the indigenous hydroxyls and create new isolated hydroxyls by grafting organo-trimethoxysilane (R-TMS) to generate a silicon T 2 site, (≡SiO-) 2 SiR(-OH). We used in situ 29 Si solid-state NMR experiments to monitor the evolution of T n sites, (≡SiO - ) n SiR(-OH) 3-n (n = 1, 2, 3). The study indicates that i) the grafting proceeds in a consecutive manner as T 1 → T 2 → T 3 , and ii) the kinetics depend on the type of functional groups in the silane. However, the rates of T 1 formation and T 2 → T 3 conversion are also controlled to a significant extent by the entropy loss associated to the initial silane binding and the spatial arrangement of surface hydroxyls, respectively. The grafting of R-TMS with a basic functional group leads to a lower concentration of T 1 sites. The nucleophilicity of the functional group facilitates the grafting process by lowering the enthalpy barrier, while the T 1 formation rate is more influenced by the entropy barrier than the T 1 → T 2 conversion rate. Thus, the basic functional group promotes the T 1 → T 2 conversion more than the T 1 formation, resulting in a lower concentration of T 1 sites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2020

Cross-field electron diffusion due to the coupling of drift-driven microinstabilities

In this paper, the nonlinear interaction between kinetic instabilities driven by multiple ion beams and magnetized electrons is investigated. Electron diffusion across magnetic field lines is enhanced by the coupling of plasma instabilities. Here, a two-dimensional collisionless particle-in-cell simulation is performed accounting for singly and doubly charged ions in a cross-field configuration. Consistent with prior linear kinetic theory analysis and observations from coherent Thomson scattering experiments, the present simulations identify an ion-ion two-stream instability due to multiply charged ions (flowing in the direction parallel to the applied electric field) which coexists with the electron cyclotron drift instability (propagating perpendicular to the applied electric field and parallel to the ExB drift). Small-scale fluctuations due to the coupling of these naturally driven kinetic modes are found to be a mechanism that can enhance cross-field electron transport and contribute to the broadening of the ion velocity distribution functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Ultraviolet to far-infrared dielectric function of n -doped cadmium oxide thin films

Spectroscopic ellipsometry and Fourier transform infrared spectroscopy were applied to extract the ultraviolet to far-infrared (150 – 33333 cm -1 ) complex dielectric functions of high-quality, sputtered indium-doped cadmium oxide (In:CdO) thin crystalline films on MgO substrates possessing carrier densities (N d ) ranging from 1.1 x 10 19 cm -1 to 4.1 x 10 20 cm -1 . A multiple oscillator fit model was used to identify and analyze the three major contributors to the dielectric function and their dependence on doping density: interband transitions in the visible, free-carrier excitations (Drude response) in the near- to far-infrared, and IR-active optic phonons in the far-infrared. More specifically, values pertinent to the complex dielectric function such as the optical band gap (E g ), are shown here to be dependent upon carrier density, increasing from approximately 2.5–3 eV, while the high-frequency permittivity (ε ∞ ) decreases from 5.6 to 5.1 with increasing carrier density. The plasma frequency (ω p ) scales as $\sqrt{N_d}$, resulting in ω p values occurring within the mid- to near-IR, and the effective mass (m *) was also observed to exhibit doping density-dependent changes, reaching a minimum of 0.11 m o in unintentionally doped films (1.1 x 10 19 cm -1 ). Good quantitative agreement with prior work on polycrystalline, higher-doped CdO films is also demonstrated, illustrating the generality of the results. The analysis presented here in this paper will aid in predictive calculations for CdO-based next-generation nanophotonic and optoelectronic devices, while also providing an underlying physical description of the key properties dictating the dielectric response in this atypical semiconductor system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Restoration of Geographically Isolated Wetlands: An Amphibian-Centric Review of Methods and Effectiveness

Geographically isolated wetlands provide a critical habitat for pond-breeding amphibians, a taxa of broad conservation concern. Global wetland loss and degradation has made restoration essential for amphibian conservation. Restoration goals typically include recovering the wetlands’ physiochemical, hydrological, and ecological functions. However, for pond-breeding amphibians, successful restoration should also result in sustained populations, which is difficult to assess and infrequently reported. In this paper, we review the available evidence that restoration of geographically isolated wetlands promotes pond-breeding amphibian occupancy and population persistence. We provide an overview of restoration practices addressing hydrology, vegetation, and ecological processes within these unique environments and across spatial scales. We then summarize the evidence, and discuss the limitations, for evaluating successful restoration within the context of amphibian conservation across these categories. Finally, we provide recommendations for researchers and practitioners to leverage prior successes and establish systematic data collection and dissemination. Moving restoration of wetlands for amphibian conservation forward will require more robust data collection and reporting.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Ka-Band ARM Zenith Radar Corrections (KAZRCOR, KAZRCFRCOR) Value-Added Products

The Ka-Band Atmospheric Radiation Measurement (ARM) Zenith Radar Corrections (KAZRCOR) and KA-Band ARM Zenith Radar CF-Radial, Corrected (KAZRCFRCOR) value-added products (VAPs) perform several corrections to the ingested KAZR moments and also create a significant detection mask for each radar mode. The VAPs compute gaseous attenuation as a function of time and radial distance from the radar antenna, based on ambient meteorological observations, and correct observed reflectivities for that effect. Mean Doppler velocities are dealiased to correct velocities whose magnitudes exceed the radar’s Nyquist velocity. Input KAZR data fields are passed through to the KAZRCOR or KAZRCFRCOR output files, in their native time and range coordinates. Complementary corrected reflectivity and velocity fields are provided, along with a mask of significant detections and a number of data quality flags. This report covers the KAZRCOR VAP as applied to the original KAZR radars and the upgraded KAZR2 radars. Originally, prior to late 2019, there were two separate code bases for the different radar versions. Following the harmonization of KAZR and KAZR2 data formats in 2019, only a single code base is required. The new combined KAZR and KAZR2 code base is called the KAZRCFRCOR VAP. The ‘cfr portion of the VAP name refers to the Radial Climate and Forecasting data format standards that are used in these data sets. Throughout this report, references to ‘KAZRCOR’ should be taken to apply to the ‘KAZRCFRCOR’ VAP as well unless there is an explicit statement to the contrary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Performance of Embedded Sensors in 3D Printed SiC

This report summarizes recent advances in embedding sensors in 3D printed silicon carbide (SiC) ceramic components under the Transformational Challenge Reactor (TCR) program. The additive manufacturing technologies developed under this program will enable fabrication of complex structures with embedded fuels and sensors. The sensors will be capable of characterizing fuel performance using spatially distributed measurements at the most strategic locations that would be otherwise inaccessible using conventional manufacturing techniques. While previous programmatic updates describe initial concepts for embedding sensors, materials selection, and initial characterization of the embedded sensors, the technology requires further demonstration, and quality-significant procedures must be established before the embedding technology is ready for adoption by industry. To this end, this report describes the most effective techniques that have been used to embed functional sensors in 3D printed components using a combination of binder-jet additive manufacturing and chemical vapor infiltration (CVI). A detailed procedure describes each step in the process and is available upon request. Molybdenum (Mo)-sheathed thermocouples have been successfully embedded in complex SiC components, and temperatures were monitored in situ during the embedding process. Post-embedding measurements showed no significant hysteresis, and characterization of the interface revealed qualitatively strong bonding around the entire perimeter of the sensor sheath. Distributed fiber-optic temperature sensors were able to briefly measure temperature profiles during CVI, but they ultimately failed prior to completion of the CVI run. The failure appears to be related to the fiber coating at temperatures close to 1,000°C. Future work will focus on irradiation testing of embedded thermocouples and other sheathed electrical sensors, as well as the identification of fiber-optic sensor coatings that can survive CVI.

42 ENGINEERING↗
DOE OSTI2022

Atomic step disorder on polycrystalline surfaces leads to spatially inhomogeneous work functions

Structural disorder causes materials’ surface electronic properties, e.g., work function (ϕ), to vary spatially, yet it is challenging to prove exact causal relationships to underlying ensemble disorder, e.g., roughness or granularity. For polycrystalline Pt, nanoscale resolution photoemission threshold mapping reveals a spatially varying ϕ=5.70±0.03 eV over a distribution of (111) vicinal grain surfaces prepared by sputter deposition and annealing. With regard to field emission and related phenomena, e.g., vacuum arc initiation, a salient feature of the ϕ distribution is that it is skewed with a long tail to values down to 5.4 eV, i.e., far below the mean, which is exponentially impactful to field emission via the Fowler–Nordheim relation. We show that the ϕ spatial variation and distribution can be explained by ensemble variations of granular tilts and surface slopes via a Smoluchowski smoothing model wherein local ϕ variations result from spatially varying densities of electric dipole moments, intrinsic to atomic steps, that locally modify ϕ. Atomic step-terrace structure is confirmed with scanning tunneling microscopy (STM) at several locations on our surfaces, and prior works showed STM evidence for atomic step dipoles at various metal surfaces. From our model, we find an atomic step edge dipole μ=0.12 D/edge atom, which is comparable to values reported in studies that utilized other methods and materials. Our results elucidate a connection between macroscopic ϕ and the nanostructure that may contribute to the spread of reported ϕ for Pt and other surfaces and may be useful toward more complete descriptions of polycrystalline metals in the models of field emission and other related vacuum electronics phenomena, e.g., arc initiation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Inverse design of pore wall chemistry and topology through active learning of surface group interactions

Design of next-generation membranes requires a nanoscopic understanding of the effect of biologically inspired heterogeneous surface chemistries and topologies (roughness) on local water and solute behavior. In particular, the rejection of small, neutral solutes, such as boric acid, poses a heretofore unsolved challenge. In prior work, a computational inverse design technique using an evolutionary optimization successfully uncovered new surface design strategies for optimized transport of water over solutes in smooth, model pores consisting of two surface chemistries. However, extending such an approach to more complex (and realistic) scenarios involving many surface chemistries as well as surface roughness is challenging due to the expanded design space. In this work, we develop a new approach that uses active learning to optimize in a reduced feature space of surface group interactions, finding parameters that lead to their assembly into ordered, optimal patterns. This approach rapidly identifies novel surface functionalizations that maximize the difference in water and boric acid transport through the nanopore. Moreover, we find that the roughness of the nanopore wall, independent of its chemistry, can be leveraged to enhance transport selectivity: oscillations in the pore wall diameter optimally inhibit boric acid transport by creating energetic wells from which the solute must escape to transport down the pore. Furthermore, this proof-of-concept demonstrates the potential for active learning strategies, in concert with molecular simulations, to rapidly navigate complex design spaces of aqueous interfaces and is promising as a tool for engineering water-mediated surface interactions for a broad range of applications.

36 MATERIALS SCIENCE↗