Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular statics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Direct strain correlations at the single-atom level in three-dimensional core-shell interface structures

Abstract Nanomaterials with core-shell architectures are prominent examples of strain-engineered materials. The lattice mismatch between the core and shell materials can cause strong interface strain, which affects the surface structures. Therefore, surface functional properties such as catalytic activities can be designed by fine-tuning the misfit strain at the interface. To precisely control the core-shell effect, it is essential to understand how the surface and interface strains are related at the atomic scale. Here, we elucidate the surface-interface strain relations by determining the full 3D atomic structure of Pd@Pt core-shell nanoparticles at the single-atom level via atomic electron tomography. Full 3D displacement fields and strain profiles of core-shell nanoparticles were obtained, which revealed a direct correlation between the surface and interface strain. The strain distributions show a strong shape-dependent anisotropy, whose nature was further corroborated by molecular statics simulations. From the observed surface strains, the surface oxygen reduction reaction activities were predicted. These findings give a deep understanding of structure-property relationships in strain-engineerable core-shell systems, which can lead to direct control over the resulting catalytic properties.

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

Charge influence on point defect interactions with xenon bubbles in uranium dioxide

The interaction of xenon (Xe) bubbles with small self-defects in uranium dioxide (UO2) has been studied using molecular statics simulations. The results show that the pressure and size of bubbles have a minimal impact on the heterogeneity of the interaction between a charge neutral Xe bubble and charged point defects, while the local charge distribution around the bubble overwhelmingly determines the interaction strength and critical interaction distance. The charge effect on the interactions between Xe bubbles and defects was further confirmed by assessing the point defect interaction energies with negatively or positively charged Xe bubbles. The Xe density, or pressure, has a much smaller effect on the interaction energy between charged bubbles and small defects at least for low and medium pressure bubbles. The interactions between charged Xe bubbles and point defects basically follow the Coulomb electrostatic interaction law, which is independent of the empirical potentials employed for the Xe–UO2 system in this work.

Yang, L. (ORCID:0000000322166071)↗
DOE OSTI2026

Extended Rice–Thomson analysis and atomistic simulations revealing grain boundary effects on fracture in refractory high-entropy alloys

Significance This work serves to extend the fundamental ductile vs. brittle fracture theory, specifically the Rice–Thomson criterion, by introducing a grain boundary ahead of an initiating crack which propagates at an oblique angle to impinge the boundary. Atomistic fracture simulations on two refractory complex concentrated alloys, the brittle NbMoTaW and the ductile Nb 45 Ta 25 Ti 15 Hf 15 , demonstrate qualitative correspondence with the extended Rice–Thomson criterion and experimental observations. Abstract Understanding how grain boundaries mediate fracture remains a critical challenge in designing ductile, high-performance refractory alloys. Here, we extend the Rice–Thomson criterion to account for the angle between cracks and the impinging grain boundaries (GBs), capturing the competition between intergranular fracture and dislocation-mediated plasticity. Using machine learning interatomic potentials, we performed molecular statics simulations to probe fracture mechanisms in nanocrystalline NbMoTaW and Nb 45 Ta 25 Ti 15 Hf 15 , each with two different grain sizes, revealing trends consistent with experimental observations and the extended Rice model. Comparison with averaged R-curves for bulk samples demonstrates that GBs enhance ductility in Nb 45 Ta 25 Ti 15 Hf 15 in both grain sizes investigated. In contrast, GBs only locally improve fracture resistance in NbMoTaW when cracks are temporarily pinned at GBs inclined at high angles from the crack, but generally promote brittle intergranular fracture. These contrasting behaviors are attributed to differences in GB cohesion, reflecting clear alloying trends that align with ab-initio calculations and trends observed experimentally. Our results bridge classical fracture theory, atomistic simulations, and experimental observations, providing a comprehensive understanding of the fracture mechanisms in nanocrystalline refractory complex concentrated alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accelerated molecular dynamics simulations of dislocation climb in nickel

We report the mechanical behavior of materials operating under high temperatures is strongly influenced by creep mechanisms such as dislocation climb, which is controlled by the diffusion of vacancies. However, atomistic simulations of these mechanisms have traditionally been impractical due to the long time scales required. To overcome these time scale challenges, we use Parallel Trajectory Splicing (ParSplice), an accelerated molecular dynamics method, to simulate dislocation climb in nickel. We focus on modeling the activity of a vacancy near a jog on an edge dislocation in order to observe vacancy pipe diffusion and vacancy absorption at the jog. From rigorously constructed trajectories encompassing more than 2000 vacancy absorption events over a simulation time of more than 4μs at 900 K, a comprehensive sampling of available atomistic mechanisms is collated and analyzed further with molecular statics calculations. We estimate average rates for pipe diffusion and vacancy absorption into the jog using data from the dynamic and static calculations, finding very good agreement. Our results strongly suggest that the dominant mechanism for vacancy absorption by jogs is via biased diffusion to the dislocation core followed by fast pipe diffusion to the jog.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗
DOE OSTI2022

Invariant surface elastic properties in FCC metals and their correlation to bulk properties revealed by machine learning methods

In this work, we present a combination of machine-learned models that predicts the surface elastic properties of general free surfaces in face-centered cubic (FCC) metals. These models are built by combining a semi-analytical method based on atomistic simulations to calculate surface properties with the artificial neural network (ANN) method or the boosted regression tree (BRT) method. The latter is also used to link bulk properties and surface orientation to surface properties. The surface elastic properties are represented by their invariants considering plane elasticity within a polar method. The resulting models are shown to accurately predict the surface elastic properties of seven pure FCC metals (Cu, Ni, Ag, Au, Al, Pd, Pt). The BRT model reveals the correlations between bulk and corresponding surface properties in terms of invariants, which can be used to guide the design of complex nano-sized particles, wires and films. Finally, by expressing the surface excess energy density as a function of surface elastic invariants, fast predictions of surface energy as a function of in-plane deformations can be made from these model constructs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces

The objective of this project is to develop, and integrate, high-performance simulation tools capable of predicting plasma-facing component (PFC) operating lifetime and the impact of the evolving surface morphology of tungsten-based PFCs on plasma contamination, including the dynamic recycling of fuel species and tritium retention, in future magnetic fusion devices. Establishing a fundamental physical understanding and developing predictive capabilities of plasma-surface interactions (PSI) requires simultaneously addressing complex and diverse physics occurring over a wide range of length (Angstroms to meters) and time (femtoseconds to years) scales, as well as integrating extensive physical processes across the plasma–surface interface. This requires development of not only detailed physics models and computational strategies at each scale, but also algorithms and methods to couple them effectively in a way that can be robustly validated. Deploying these tools requires the continued development and coupling of leadership-scale computational codes to describe the boundary plasma and the evolving PFC surface, as well as a host of simulations that bridge disparate scales to address complex physical and computational issues at the plasma–surface interface in multi-component materials systems for magnetic fusion energy development beyond ITER.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Next-Generation Materials Design: Quantum Mechanics and Data-Driven Modeling

The future of materials design is rapidly advancing through the combination of quantum mechanics and data-driven modeling. These approaches integrate quantum principles with advanced data analysis, enabling precise insights into material behavior. This talk will highlight recent progress in using these methods for computational design, particularly in high-entropy alloy catalysts, emphasizing the role of hierarchical machine-learning architectures for accurate predictions. Additionally, I will discuss our work on developing machine learning interatomic potentials (MLPs) for single-element metals, metal oxides, and alloys under extreme conditions, focusing on melting behavior and phase properties at high temperatures and pressures. We have also refined our MLP models to capture dynamic surface interactions, such as CO2 and CO adsorption on MgO, using both static and molecular dynamics simulations. These models maintain high accuracy while significantly reducing computational costs compared to first-principles calculations. By enabling efficient and accurate simulations, this work supports broader community adoption, optimizes datasets for materials discovery, and extends the accessible time, size, and environmental conditions beyond the limits of experiments and traditional simulations.

machine learning↗
DOE OSTI2025

Optimizing Grain Boundary Structures with LAMMPS Using Evolutionary Algorithms

Grain boundary structure optimization is an important part of materials modeling. Current methods for grain boundary structure optimization involve inefficient, time-consuming processes that do not fully explore the interface parameter space. Evolutionary algorithms have recently been demonstrated to be effective at determining both stable and metastable grain boundary interface structures. In this work, we demonstrate the use of GBOpt, a grain boundary structure optimization software designed to use the Large-scale Atomic/Molecular Massively Parallel Simulation (LAMMPS) software to efficiently determine grain boundary structures. We demonstrate that a only a few manipulations, namely atom insertion, atom removal, and relative grain displacement, are sufficient to explore much of the grain boundary structure parameter space. The efficacy of this approach is demonstrated on an FCC Ni system, and a BCC Fe system. The computational cost is compared against the gamma-surface sampling approach to demonstrate performance improvement.

Evolutionary algorithms↗
DOE OSTI2023

Reducing the matrix effect in mass spectral imaging of biofilms using flow-cell culture

The interactions between soil microorganisms and soil minerals play a crucial role in the formation and evolution of minerals and the stability of soil aggregates. Due to the heterogeneity and diversity of the soil environment, the under-standing of the functions of bacterial biofilms in soil minerals at the microscale is limited. A soil mineral-bacterial biofilm system was used as a model in this study, and it was analyzed by time-of-flight secondary ion mass spectrometry (ToF-SIMS) to acquire molecular level information. Static culture in multi-wells and dynamic flow-cell culture in microfluidics of biofilms were investigated. Our results show that more characteristic molecules of biofilms can be observed in SIMS spectra of the flow-cell culture. In contrast, biofilm signature peaks are buried under the mineral components in SIMS spectra in the static culture case. Spectral overlay was used in peak selection prior to performing Principal component analysis (PCA). Comparisons of the PCA results between the static and flow-cell culture show more pronounced molecular features and higher loadings of organic peaks of the dynamic cultured specimens. For example, fatty acids secreted from bacterial biofilm extracellular polymeric substance are likely to be responsible for biofilm dispersal due to mineral treatment up to 48 h. Such findings suggest that the use of microfluidic cells to dynamically culture biofilms be a more suitable method for reducing the matrix effect arisen from the growth medium and minerals as a perturbation fac-tor for improved spectral and multivariate analysis of complex mass spectral data in ToF-SIMS. These results show that the interaction mechanism between biofilms and soil minerals at the molecular level can be better studied using the flow-cell culture and advanced mass spectral imaging techniques like ToF-SIMS.

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

Mechanisms of Al3+ dimerization in alkaline solutions

The molecular speciation of aluminum (Al3+) in alkaline solutions is fundamental to its precipitation chemistry within a number of industrial applications that include ore renement and industrial processing of Al wastes. Under these conditions, Al3+ is predominantly Al(OH) {4 , while at high [Al] dimeric species also known to form. To date, the mechanism of dimer formation remains unclear and it is likely infuenced by complex ion ion interactions. In the present work, we investigate the role of ion pairing and a suite of potential dimerization pathways using a three-pronged approach that consists of density functional theory-based molecular dynamics (DFTMD), static DFT calculations (sDFT) and semi-empirical density functional tight binding molecular dynamics (SEMD). Specic cation eects imparted by the background electrolyte cations Na+, Li+ and K+ have been examined. Our simulations predict that when the Al species are ion paired with either cation the formation of the oxo-bridged Al2O(OH) 2{6 is favored with respect to the dihydroxo-bridged Al2(OH) 2{8 , in agreement with previousspectroscopic work. The formation of both dimers rst proceeds by bridging of two monomeric units via one hydroxo ligand, leading to a labile Al2(OH) 2{8 isomer. The eect of contact ion pairing of Li+ and K+ on the dimerization energetics is distinctly more favorable than that of Na+, which may have an impact on further ligomerization.

Pouvreau, Maxime↗
DOE OSTI2023

Electrical Manipulation of Quantum Coherence in a Two-Level Molecular System

We report the manipulation of ultrafast quantum coherence of a two-level single hydrogen molecular system by employing static electric field from the sample bias in a femtosecond terahertz scanning tunneling microscope. A H 2 molecule adsorbed on the polar Cu 2⁢ N surface develops an electric dipole and exhibits a giant Stark effect. An avoided crossing of the quantum state energy levels is derived from the resonant frequency of the single H 2 two levels in a double-well potential. The dephasing time of the initial wave packet can also be changed by applying the electric field. The electrical manipulation for different tunneling gaps in three dimensions allows quantification of the surface electrostatic fields at the atomic scale. Further, our work demonstrated the potential application of molecules as controllable two-level molecular systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamically tunable membrane metasurfaces for infrared spectroscopy and strong light-matter interactions

Mid-infrared spectroscopy enables biochemical sensing by identifying vibrational molecular fingerprints, but it faces limitations in instrumentation portability and analytical sensitivity. Optical metasurfaces with strong mid-infrared photonic resonances provide an attractive solution towards on-chip spectrometry and sensitive molecular detection, yet their static nature hinders their anticipated impact. Here, we introduce and demonstrate dynamically tunable silicon membrane metasurfaces exhibiting high-Q transmissive resonances in the fingerprint region. By harnessing silicon’s thermo-optical properties, we achieve continuous modulation of coupling-induced transparency (CIT) modes that emerge upon the interference of quasi-bound states in the continuum (q-BICs) and surface lattice modes (SLMs). We measure a spectral tuning rate of 0.06 cm −1 K −1 by continuously sweeping the sharp CIT resonances over a 23.5 cm −1 spectral range across a temperature range of 300–700 K. In the current proof‑of‑concept implementation, the dynamic transmission control enables non-contact chemical analysis of polymer films by detecting characteristic absorption bands of polystyrene (1450 and 1492 cm −1 ) and poly(methyl methacrylate) (1730 cm −1 ) without requiring conventional spectrometers. When analyte molecules fill the metasurface-generated photonic cavities, we demonstrate vibrational strong coupling between the poly(methyl methacrylate)’s carbonyl band and the CIT mode, manifested in a Rabi splitting of ~43 cm −1 . Our results establish a new photonic platform that unites spectral precision, strong field enhancement, and reconfigurability, offering diverse potential for compact mid-infrared spectroscopy, molecular sensing, and programmable polaritonic photonics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

A theoretical study of solid solution strengthening in the refractory medium entropy alloy Nb 45 Ta 25 Ti 15 Hf 15

The refractory medium-entropy alloy (RMEA) Nb 45 Ta 25 Ti 15 Hf 15 exhibits exceptional tensile ductility and fracture toughness at ambient temperature, but its engineering applications are limited by a lack of high temperature strength. Using a machine-learning interatomic potential (MLIP) with near-density functional theory (DFT) accuracy, we conducted molecular dynamics (MD) and statics simulations of the behavior of dislocations with both screw and edge characters. We also analyze experimentally measured yield strengths using the Rao-Suzuki model and the Maresca-Curtin model modified to include a temperature-dependent shear modulus and a bulk modulus-dependent misfit volume, thereby uncovering the mechanisms underlying the yielding of this RMEA. Compared with the published experimental yield strength, the models parameterized by the MLIP effectively reproduce the experimental results over a wide temperature range. The models and MD simulations indicate that yielding is governed by screw dislocations, with dipole dragging as the dominant mechanism. In MD simulations, we observed a potential softening mechanism not considered by the Rao-Suzuki screw model: slow migration of interstitial jogs along the dislocation core, which could lead to the annihilation of vacancy and interstitial jog pairs by their combination.

BCC complex concentrated alloys↗
DOE OSTI2025

MP-ALOE: an r2SCAN dataset for universal machine learning interatomic potentials

We present MP-ALOE, a dataset of nearly 1 million DFT calculations using the accurate r2SCAN meta-generalized gradient approximation. Covering 89 elements, MP-ALOE was created using active learning and primarily consists of off-equilibrium structures. We benchmark a machine learning interatomic potential trained on MP-ALOE, and evaluate its performance on a series of benchmarks, including predicting the thermochemical properties of equilibrium structures; predicting forces of far-from-equilibrium structures; maintaining physical soundness under static extreme deformations; and molecular dynamic stability under extreme temperatures and pressures. MP-ALOE shows strong performance on all of these benchmarks and is made public for the broader community to utilize.

Kuner, Matthew C↗