Search NASA⌕ Search

SEARCH · Search NASA

Results for “edge 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 1,585 records · Page 88

Efficient frequency allocation for superconducting quantum processors using improved optimization techniques

Building on previous research on frequency allocation optimization for superconducting circuit quantum processors, this work incorporates several techniques to improve overall solution quality. Here, we introduce constraints and imposed edgewise differences help to improve the optimization results. We also introduce optimization variables for the orientation of each edge, defined as the direction from the control qubit to the target qubit, to be chosen during optimization. To scale up to larger processors, multimodule designs are employed with various boundary conditions, thereby enhancing the collective yield. These enhancements allow for greater flexibility in processor design by eliminating the need for handpicked orientations. We support the efficient assembly of large processors with dense connectivity by choosing the best boundary conditions. Examples demonstrate that, at low computational cost, this optimization approach finds a frequency configuration for a square chip with over 1000 qubits and over 10% yield at much larger dispersion levels than required by previous approaches.

Zhang, Zewen [Argonne National Laboratory (ANL), A↗

Post-Marcus Theory and Simulation of Interfacial Charge Transfer Dynamics in Organic Semiconducting Materials (Final Report)

The research program developed, validated and applies predictive computational protocols for calculating charge transfer (CT) rates in complex molecular systems, including molecular dyads and triads in liquid solution and solid-state organic semiconducting (OSC) materials. We have established a transformative computational scheme that goes beyond widely used simplifications, to achieve realistic descriptions of CT processes. The approach properly addresses the contribution of molecular environment at ambient conditions to CT processes. Our approach achieves unique insight on investigated CT processes in relevant experimental efforts. The collaborative team included three principal investigators (PIs), with complimentary expertise in classical molecular dynamics simulations and data science (Cheung), state-of-the-art electronic structure calculations (Dunietz), and cutting-edge theory and simulation techniques for modeling energy, charge and coherence transfer dynamics in molecular systems (Geva).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exceptionally Complex Electronic Structures of Lanthanide Oxides and Small Molecules

Lanthanide (Ln) oxide clusters and molecular systems provide a bottom-up look at the electronic structures of the bulk materials because of close parallels in the patterns of Ln 4fN subshell occupancy between the molecular and bulk Ln 2 O 3 size limits. At the same time, these clusters and molecules offer a challenge to the theory community to find appropriate and robust treatments for the 4f N patterns across the Ln series. Anion photoelectron (PE) spectroscopy provides a powerful experimental tool for studying these systems, mapping the energies of the ground and low-lying excited states of the neutral relative to the initial anion state, providing spectroscopic patterns that reflect the Ln 4f N occupancy. In this Account, we review our anion PE spectroscopic and computational studies on a range of small lanthanide molecules and cluster species. The PE spectra of LnO – (Ln = Ce, Pr, Sm, Eu) diatomic molecules show spectroscopic signatures associated with detachment of an electron from what can be described as a diffuse Ln 6s-like orbital. While the spectra of all four diatomics share this common transition, the fine structure in the transition becomes more complex with increasing 4f occupancy. This effect reflects increased coupling between the electrons occupying the corelike 4f and diffuse 6s orbitals with increasing N. Understanding the PE spectra of these diatomics sets the stage for interpreting the spectra of polyatomic molecular and cluster species. In general, the results confirm that the partial 4f N subshell occupancy is largely preserved between molecular and bulk oxides and borides. However, they also suggest that surfaces and edges of bulk materials may support a low-energy, diffuse Ln 6s band, in contrast to bulk interiors, in which the 6s band is destabilized relative to the 5d band. We also identify cases in which the molecular Ln centers have 4f N+1 occupancy rather than bulklike 4f N , which results in weaker Ln–O bonding. Specifically, Sm centers in mixed Ce–Sm oxides or in Sm x O y – (y ≤ x) clusters have this higher 4f N+1 occupancy. The PE spectra of these particular species exhibit a striking increase in the relative intensities of excited-state transitions with decreasing photon energy (resulting in lower photoelectron kinetic energy). Furthermore, this is opposite of what is expected on the basis of the threshold laws that govern photodetachment. We relate this phenomenon to strong electron–neutral interactions unique to these complex electronic structures. The time scale of the interaction, which shakes up the electronic configuration of the neutral, increases with decreasing electron momentum. From a computational standpoint, we point out that special care must be taken when considering Ln cluster and molecular systems toward the center of the Ln series (e.g., Sm, Eu), where treatment of electrons explicitly or using an effective core potential can yield conflicting results on competing subshell occupancies. However, despite the complex electronic structures associated with partially filled 4f N subshells, we demonstrate that inexpensive and tractable calculations yield useful qualitative insight into the general electronic structural features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

Resonant X-ray Diffraction Reveals the Location of Counterions in Doped Organic Mixed Ionic Conductors

Organic mixed ionic–electronic conductors (OMIECs) have the potential to enable diverse new technologies, ranging from novel in situ biosensors to flexible energy storage devices and neuromorphic computing platforms. However, their complex behavior in functional films involving electrolyte-induced swelling, ion ingress, and electrochemical doping inhibits rational material design. Of critical importance is an understanding of the specific location of the ions in the volumetrically doped material, yet this information is not readily available. Here, in this report, we present the use of grazing- incidence resonant X-ray diffraction (RXRD, also known as anomalous diffraction) at S and Cl K-edges to determine the structure of a doped, prototypical, semicrystalline polymer OMIEC based on oligo(ethylene glycol) substitution of regioregular polythiophene. The RXRD measurement provides two key insights. Quantitative analysis of the RXRD allows the determination of the position of the ion relative to the polymer backbone in the crystalline regions. We find that the anion is relatively distant from the backbone, nearer to the lamella mid-plane naively in conflict with expected Coulombic attraction between the ion and the doped polymer polaron. Comparison of RXRD to Cl – fluorescence (total Cl – ) allows determination of the relative order of doping between the crystalline and amorphous regions. We find preferential doping of the crystalline regions. Both insights, the preferential doping of crystals at low potential and the specific location of the counterion with respect to the polymer backbone, are critical to developing a microscopic understanding of transport in OMIECs.

36 MATERIALS SCIENCE↗

Simultaneous use of camera and probe diagnostics to unambiguously identify and study the dynamics of multiple underlying instabilities during the route to plasma turbulence

We use multiple-tip Langmuir probes and fast imaging to unambiguously identify and study the dynamics of underlying instabilities during the controlled route to fully-developed plasma turbulence in a linear magnetized helicon plasma device. Langmuir probes measure radial profiles of electron temperature, plasma density and potential; from which we compute linear growth rates of instabilities, cross-phase between density and potential fluctuations, Reynold's stress, particle flux, vorticity, time-delay estimated velocity, etc. Fast imaging complements the 1D probe measurements by providing temporally and spatially resolved 2D details of plasma structures associated with the instabilities. We find that three radially separated plasma instabilities exist simultaneously. Density gradient driven resistive drift waves propagating in the electron diamagnetic drift direction separate the plasma into an edge region dominated by strong, velocity shear driven Kelvin-Helmholtz instabilities and a central core region which shows coherent Rayleigh-Taylor modes propagating in the ion diamagnetic drift direction. The simultaneous, complementary use of both probes and camera was crucial to identify the instabilities and understand the details of the very rich plasma dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interior region-of-interest reconstruction using a small, nearly piecewise constant subregion

Purpose: To develop a method to reconstruct an interior region-of-interest (ROI) image with sufficient accuracy that uses differentiated backprojection (DBP) projection onto convex sets (POCS) [H. Kudo et al., ''Tiny a priori knowledge solves the interior problem in computed tomography'', Phys. Med. Biol. 53, 2207-2231 (2008)] and a tiny knowledge that there exists a nearly piecewise constant subregion. Methods: The proposed method first employs filtered backprojection to reconstruct an image on which a tiny region P with a small variation in the pixel values is identified inside the ROI. Total variation minimization [H. Yu and G. Wang, ''Compressed sensing based interior tomography'', Phys. Med. Biol. 54, 2791-2805 (2009); W. Han et al., ''A general total variation minimization theorem for compressed sensing based interior tomography'', Int. J. Biomed. Imaging 2009, Article 125871 (2009)] is then employed to obtain pixel values in the subregion P, which serve as a priori knowledge in the next step. Finally, DBP-POCS is performed to reconstruct f(x,y) inside the ROI. Clinical data and the reconstructed image obtained by an x-ray computed tomography system (SOMATOM Definition; Siemens Healthcare) were used to validate the proposed method. The detector covers an object with a diameter of {approx}500 mm. The projection data were truncated either moderately to limit the detector coverage to diameter 350 mm of the object or severely to cover diameter 199 mm. Images were reconstructed using the proposed method. Results: The proposed method provided ROI images with correct pixel values in all areas except near the edge of the ROI. The coefficient of variation, i.e., the root mean square error divided by the mean pixel values, was less than 2.0% or 4.5% with the moderate or severe truncation cases, respectively, except near the boundary of the ROI. Conclusions: The proposed method allows for reconstructing interior ROI images with sufficient accuracy with a tiny knowledge that there exists a nearly piecewise constant subregion.

60 APPLIED LIFE SCIENCES↗

(Towards) DNS of a Laboratory Lean CH4/H2 Low-Swirl Flame Impinging on an Inclined Wall

Due to downsizing trends, flame-wall interaction (FWI) is increasingly prominent in gas turbines (GTs). FWI has direct consequences on flame stabilization and pollutant emissions, but it is not well understood in turbulent flows representative of GTs. We present results from a direct numerical simulation (DNS) of a turbulent CH4/H2 model GT low-swirl laboratory-scale flame interacting with an inclined wall. The results from the laboratory flame include simultaneous measurements of velocity using stereo particle imaging velocimetry and OHxCH2O planar laser induced fluorescence. The adaptive-mesh refinement solver PeleLMeX is used, with 24-species, 105-reaction reduced Aramco chemical kinetics mechanism. The premixed fuel-air mixture consists of hydrogen-enriched methane with 70% hydrogen volume fraction and 0.4 equivalence ratio. The inflow is prescribed to match experimental measurements at the burner exit. Karlovitz and turbulent Reynolds numbers are 300 and 400, respectively. The simulation and experimental results show excellent agreement. The flame features a bowl-shape stabilization, with a corrugated, continuous flame front at the leading edge, followed by fragmented reaction zones downstream. A large diffuse cloud of CH2O is formed downstream of the quenching point. The simulation results indicate that the cloud of CH2O is the result of incomplete methane combustion, with CH2O "leaking" from the locally quenched reaction zones.The DNS provides fine-grain resolution of turbulence-flame-wall interaction that cannot be captured with experimental measurements. With access to the entire solution vector at each cell of the computational domain, the local quenching.

flame-wall interactions↗

A flux-balanced fluid model for collisional plasma edge turbulence: Numerical simulations with different aspect ratios

We investigate the drift wave–zonal flow dynamics in a shearless slab geometry with the new flux-balanced Hasegawa-Wakatani model. As in previous Hasegawa-Wakatani models, we observe a sharp transition from a turbulence dominated regime to a zonal jet dominated regime as we decrease the plasma resistivity. However, unlike previous models, zonal structures are always present in the flux-balanced model, even for high resistivity, and strongly reduce the level of particles and vorticity flux. The more robust zonal jets also have a higher variability than in previous models, which is further enhanced when the computational domain is chosen to be elongated in the radial direction. In these cases, we observe complex multiscale dynamics, with multiple jets interacting with one another, and intermittent bursts. We present a detailed statistical analysis which highlights how the changes in the aspect ratio of the computational domain affect the third-order statistical moments, and thus modify the turbulent dynamics. By changing the aspect ratio and extending either the radial or the binormal direction, the present model study offers a better approximation to mimic turbulence in flux tube simulations for either very low magnetic shear or very high magnetic shear systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ensemble Federated Machine Learning‐Based Cybersecurity Situational Awareness in Microgrid Network

Cyber-physical microgrids are vulnerable to stealthy cybersecurity threats that disguise their actions through the exploitation of system knowledge. Such actions can severely impacts microgrids deployed in defense bases, slowing the response time of military forces during national emergencies. Several machine-learning algorithms have been proposed to detect intrusions in the grid networks; however, these traditional machine-learning algorithms lack data privacy and are subject to several adversarial machine-learning threats. This paper proposes a novel federated machine learning (FML)-based three-model framework to detect and identify stealthy data-integrity attacks while ensuring data privacy in microgrid networks. The proposed architecture uses a variational mode decomposition technique to extract derived features from incoming measurement and control datasets. The extraction of these derived features allows FML models to learn minute variations in data patterns that allow them to perform significantly better than the models trained with generic datasets consisting of raw features. Our experimental results show the efficient performance of the proposed methodology against different types of data integrity attacks while considering primary and secondary controllers in microgrids. Further, the applied FML-integrated random forest ensemble algorithm outperforms the existing generic FML algorithms during noisy and noise-free datasets with prediction latencies of only 91–134 µs per sample within the 0.1 s sampling interval and requires communication bandwidth of around ∼8.25 KB/s at the control center and ∼2.7 KB/s per edge client for communication.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Applications of high-speed dust injection to magnetic fusion [Slides]

It is now an established fact that a significant amount of dust is produced in magnetic fusion devices due to plasma-wall interactions. Dust inventory must be controlled, in particular for the next-generation steady-state fusion machines like ITER, as it can pose significant safety hazards and degrade performance. Safety concerns are due to tritium retention, dust radioactivity, toxicity, and flammability. Performance concerns include high-Z impurities carried by dust to the fusion core that can reduce plasma temperature and may even induce sudden termination of the plasma. We have recognized that dust transport, dust-plasma interactions in magnetic fusion devices can be effectively studied experimentally by injection of dust with known properties into fusion plasmas. Other applications of injected dust include diagnosis of fusion plasmas and edge localized mode (ELM)'s pacing. In diagnostic applications, dust can be regarded as a source of transient neutrals before complete ionization. ELM's pacing is a promising scheme to prevent disruptions and type I ELM's that can cause catastrophic damage to fusion machines. Different implementation schemes are available depending on applications of dust injection. One of the simplest dust injection schemes is through gravitational acceleration of dust in vacuum. Experiments at Los Alamos and Princeton will be described, both of which use piezoelectric shakers to deliver dust to plasma. In Princeton experiments, spherical particles (40 micron) have been dropped in a systematic and reproducible manner using a computer-controlled piezoelectric bending actuator operating at an acoustic (0,2) resonance. The circular actuator was constructed with a 2.5 mm diameter central hole. At resonance ({approx} 2 kHz) an applied sinusoidal voltage has been used to control the flux of particles exiting the hole. A simple screw throttle located {approx}1mm above the hole has been used to set the magnitude of the flux achieved for a given voltage. Particle fluxes ranging from a few tens of particle per second up to thousands of particles per second have been achieved using this simple device. To achieve higher dust injection speed, another key consideration is how to accelerate dust at controlled amount. In addition to gravity, other possible acceleration mechanisms include electrostatic, electromagnetic, gas-dragged, plasma-dragged, and laser-ablation-based acceleration. Features and limitations of the different acceleration methods will be discussed. We will also describe laboratory experiments on dust acceleration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

District-Scale Analysis of Electricity Load and Strategies to Improve Energy Reliability Using Prototype District Models

Projected increases in electricity demand in the U.S. highlight the urgent need for effective load management to ensure grid reliability. As the building sector accounts for approximately 75% of electricity usage, enhancing energy efficiency and flexibility in this sector is crucial. Adopting district-level approaches offers significant advantages over traditional individual building analyses by enabling shared infrastructure and economies of scale. To navigate the data and computational challenges associated with modeling energy at the district level, prototype district models have been proposed as holistic, system-level solutions that capture complex interactions within typical configurations. This study presents these models as a reference tool for analyzing district-scale energy systems across various climate zones in the U.S. Developed with input from stakeholders, these models integrate varied building characteristics, inter-building connections, and energy system interactions. A case study utilizing the Urban Edge prototype district model, implemented on the URBANopt™ platform, evaluates multiple demand scenarios and the impact of distributed energy resources such as fuel-fired backup generators, photovoltaic systems, and batteries. Findings suggest that while new electric systems can significantly reduce annual energy use, they may also elevate peak electricity loads, with a notable 43% increase in heating-dominant climate zone 5B. The optimal backup power solutions vary based on location, influenced by factors such as utility rates and incentives. For example, PV and batteries perform well in high-cost regions like New York City, while diesel backup generators are more suitable for backup needs in climate zone 3A, such as Atlanta. Thus, this research highlights the importance of prototype district models for future district-scale energy planning.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Defect Properties, Anion Ordering, and Photochromic Mechanism in Yttrium Oxyhydride

Yttrium oxyhydride (YHO) undergoes a reversible photochromic transition when exposed to ultraviolet light. However, the mechanism for this transformation is not fully understood, and the structure and precise chemical composition of YHO remain under debate. Here, we use first-principles density functional theory calculations with a hybrid functional to study the structure, chemical stability, and point defect properties of YHO. As experiments have shown, we find that YHO prefers a cubic structure, with H and O anions present in equal concentrations and located on tetrahedral sites. Stoichiometric and ordered YHO is chemically stable, but it has a wide band gap of 5.01 eV, considerably larger than that measured in experiments (2.4–3.8 eV). On the other hand, Y4H10O has a smaller band gap of 2.97 eV and also has a region of chemical stability; thus, the actual material may include some fraction of this H-rich structure. The defect chemistry of YHO is dominated by anionic antisite species (H O and O H ), with hydrogen interstitials (H i ) and vacancies (V H ) also present in reasonably high concentrations. We show that antisite disorder lowers the band gap relative to the perfectly ordered structure, bringing the magnitude of the gap into closer agreement with experiment. Based on our calculations of defect migration and the positions of defect states relative to the band edges, we link the onset of photochromic behavior to the reaction H O – → V O 0 + H i – , which follows photoexcitation of a H O + defect. H i – can subsequently migrate away and be trapped by additional H O + defects, contributing to the persistence of the reaction, while the resultant oxygen vacancy, V O 0 , introduces an occupied defect state that leads to optical absorption at visible wavelengths. Our results can explain reported discrepancies between experimental and computational results for YHO, and they allow us to propose specific atomic-scale processes that can lead to photochromism. In conclusion, understanding these mechanisms is key for unlocking YHO’s application in devices ranging from smart windows and optoelectronics to electrochemical synapses for neural networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-F-BRD-14: Dose Weighted Linear Energy Transfer Analysis of Critical Structures in Proton Therapy of Pediatric Brain Tumor Patients

Purpose: Dose weighted Linear Energy Transfer (LETd) analysis of critical structures may be useful in understanding the side effects of the proton therapy. The objective is to analyze the differences between LETd and dose distributions in brain tumor patients receiving double scattering proton therapy, to quantify LETd variation in critical organs, and to identify beam arrangements contributing to high LETd in critical organs. Methods: Monte Carlo simulations of 9 pediatric brain tumor patients were performed. The treatment plans were reconstructed with the TOPAS Monte Carlo code to calculate LETd and dose. The beam data were reconstructed proximal to the aperture of the double scattering nozzle. The dose and LETd to target and critical organs including brain stem, optic chiasm, lens, optic nerve, pituitary gland, and hypothalamus were computed for each beam. Results: Greater variability in LETd compared to dose was observed in the brainstem for patients with a variety of tumor types including 5 patients with tumors located in the posterior fossa. Approximately 20%–44% brainstem volume received LETd of 5kev/µm or greater from beams within gantry angles 180°±30° for 5 patients treated with a 3 beam arrangement. Critical organs received higher LETd when located in the vicinity of the beam distal edge. Conclusion: This study presents a novel strategy in the evaluation of the proton treatment impact on critical organs. While the dose to critical organs is confined below the required limits, the LETd may have significant variation. Critical organs in the vicinity of beam distal edge receive higher LETd and depended on beam arrangement, e.g. in posterior fossa tumor treatment, brainstem receive higher LETd from posterior-anterior beams. This study shows importance of the LETd analysis of the radiation impact on the critical organs in proton therapy and may be used to explain clinical imaging observations after therapy.

60 APPLIED LIFE SCIENCES↗

Growth of Well-Defined Model Catalysts for Electrochemistry: From Surface Science Studies to Electrocatalytic CO2 Conversion

We combined ultrahigh vacuum (UHV) surface science techniques, electrochemical measurements, and computational modeling to investigate electrocatalytic systems that are crucial components of the carbon management effort, including oxygen evolution reactions (OER) and CO2 reduction reactions (CO2RR). Well-defined Model catalysts were grown on substrates in the UHV chamber and characterized with X-ray photoelectron spectroscopy (XPS) and scanning tunneling microscopy (STM). Selected model electrocatalysts were then tested in electrochemical cells to establish the structure-property relationships in OER and CO2RR. Our results showed that the edge sites of Fe2O3 grown on Au(111) were the most active toward OER and incorporation of Ni at the edge sites (NiFeOx) further boosted their OER activity. We also resolved the size-dependent electrocatalytic CO2-to-CO conversion of the Ag nanoparticle electrocatalysts with average particle diameter between 2 to 6 nm: smaller diameter (< 3 nm) particles favored H2 evolution reaction (HER) due to a high population of Ag edge sites, whereas larger diameter particles favored CO2RR as the population of Ag(100) surface sites grew. We further discovered that electronic interactions between small diameter Ag particles and highly defective carbon supports could break the size-dependent CO2RR selectivity, resulting in highly selective (CO Faradaic Efficiency > 90%) and active Ag nanoparticle electrocatalysts with sizes < 2 nm diameter.

Deng, Xingyi↗

Almost-Rigidity of Frameworks

Here, we extend the mathematical theory of rigidity of frameworks (graphs embedded in d-dimensional space) to consider nonlocal rigidity and flexibility properties. We provide conditions on a framework under which (I) as the framework flexes continuously it must remain inside a small ball, a property we call “almost-rigidity”; (II) any other framework with the same edge lengths must lie outside a much larger ball; (III) if the framework deforms by some given amount, its edge lengths change by a minimum amount; (IV) there is a nearby framework that is prestress stable, and thus rigid. The conditions can be tested efficiently using semidefinite programming. The test is a slight extension of the test for prestress stability of a framework, and gives analytic expressions for the radii of the balls and the edge length changes. Examples illustrate how the theory may be applied in practice, and we provide an algorithm to test for rigidity or almost-rigidity. We briefly discuss how the theory may be applied to tensegrities. © 2020 Wiley Periodicals LLC.

97 MATHEMATICS AND COMPUTING↗

Advances in Clayff Molecular Simulation of Layered and Nanoporous Materials and Their Aqueous Interfaces

As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. Overall, we conclude with a discussion of expectations for future developments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗