Establishing the temperature and orientation dependence of the threshold displacement energy in ThO2 via molecular dynamics simulations
Not Available
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
SEARCH · Search NASA
Results for “molecular simulation”
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.
Microscopic Origin of Twist-Dependent Electron Transfer Rate in Bilayer Graphene
Using molecular simulation and continuum dielectric theory, we consider how electrochemical kinetics are modulated by the twist angle in bilayer graphene electrodes. By establishing a connection between the twist angle and the screening length of charge carriers within the electrode, we investigate how tunable metallicity modifies the statistics of the electron transfer energy gap. Constant potential molecular simulations show that the activation free energy for electron transfer increases with screening length, leading to a non-monotonic dependence on the twist angle. Here, the twist angle alters the density of states, tuning the number of thermally accessible channels for electron transfer and the reorganization energy by affecting the stability of the vertically excited state through attenuated image charge interactions. Understanding these effects allows us to express the Marcus rate of interfacial electron transfer as a function of the twist angle in a manner consistent with experimental observations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Evaluation of δ-Phase ZrH1.4 to ZrH1.7 Thermal Neutron Scattering Laws Using Ab Initio Molecular Dynamics Simulations
Zirconium hydride is commonly used for next-generation reactor designs due to its excellent hydrogen retention capacity at temperatures below 1000 K. These types of reactors operate at thermal neutron energies and require accurate representation of thermal scattering laws (TSLs) to optimize moderator performance and evaluate the safety indicators for reactor design. In this work, we present an atomic-scale representation of sub-stoichiometric ZrH2−x(0.3≤x≤0.6), which relies on ab initio molecular dynamics (AIMD) in tandem with velocity auto-correlation (VAC) analysis to generate phonon density of states (DOS) for TSL development. The novel NJOY+NCrystal tool, developed by the European Spallation Source community, was utilized to generate the TSL formulations in the A Compact ENDF (ACE) format for its utility in neutron transport software. First, stoichiometric zirconium hydride cross sections were benchmarked with experiments. Then sub-stoichiometric zirconium hydride TSLs were developed. Significant deviations were observed between the new δ-phase ZrH2−x TSLs and the TSLs in the current ENDF release. It was also observed that varying the hydrogen vacancy defect concentration and sites did not cause as significant a change in the TSLs (e.g., ZrH1.4 vs. ZrH1.7) as was caused by the lattice transformation from ϵ- to δ-phase.
42 ENGINEERING↗
Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation
Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.
Boron↗
Exploration of Tertiary Structure in Sequence-Defined Polymers Using Molecular Dynamics Simulations
Peptoids are a class of sequence-defined biomimetic polymers with peptide-like backbones and side chains located on backbone nitrogens rather than alpha carbons. These materials demonstrate a strong ability for precise control of single-chain structure, multiunit self-assembly, and macromolecular assembly through careful tuning of sequence due to the diversity of available side chains, although the driving forces behind these assemblies are often not understood. Prior experimental work has shown that linked 15mer peptoids can mimic the protein helical hairpin structure by leveraging the chirality-inducing nature of bulky side chains and hydrophobicity, but there are still gaps in our understanding of the relationship between sequence, stability, and particular secondary or tertiary structure. Here, we present a molecular dynamics (MD) study on the folding behavior of these polymers into hairpins, discussing the differences in structure from sequences with various characteristics in water and acetonitrile, and then compare the handedness preference of common helical motifs between solvents.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Origin of Heating-Induced Softening and Enthalpic Reinforcement in Elastomeric Nanocomposites
Molecular simulations demonstrate that the enthalpic softening of elastomeric nanocomposites upon heating can arise naturally from a Poisson’s ratio mismatch between elastomer and nanoparticle networks, providing a more parsimonious explanation for this phenomenon than the widely accepted interpretation based on glassy interparticle bridging. Despite a century of use, the mechanism of nanoparticle-driven mechanical reinforcement of elastomers is unresolved. Here, a major hypothesis attributes it to glassy interparticle bridges, supported by an observed inversion of the variation of the modulus E(T) on heating – from entropic stiffening in elastomers to enthalpic softening in nanocomposites. Here, molecular simulations reveal that elastomer enthalpic softening can instead emerge from a competition over the preferred volumes between elastomer and nanoparticulate networks. A theory for this competition accounting for softening of the bulk modulus on heating predicts the simulated E(T) inversion, suggesting that reinforcement is driven by a volume-competition mechanism unique to cocontinuous systems of soft and rigid networks.
Biopolymers↗
Modeling and Optimization of Zeolites for Contaminant Removal from Coal Combustion Impoundment Leachates
Coal Combustion Residual (CCR) landfills in the U.S. can lead to elevated contaminant concentrations in groundwater and runoff, particularly of arsenic, boron, and selenium. Zeolites can be readily synthesized from materials contained within the coal combustion ash impoundments and can be cation-exchanged to enhance adsorption capacity, selectivity, and reactivity. However, the optimization of zeolites becomes a daunting task when accounting for the variety of Si:Al ratios, the species of extra-framework cations present, and the zeolite pore structure. Molecular simulations provide methods to study and guide the design of zeolites for the sorption of contaminants from aqueous solution. In this work, models that predict the sorption of arsenic, boron, and selenium from water using cation-exchanged zeolites. Because of the lack of experimental adsorption data for these species of contaminants, models were parametrized to reproduce the results of first-principles calculations and then used to predict sorption for zeolites for a dataset containing 6000 combinations of zeolites and sorption conditions. Machine learning was used to train a model to predict sorption for materials in this database based on the results of the molecular simulations. Next, a genetic algorithm was used to optimize zeolites for the removal of each contaminant from aqueous solution for individual impoundment sites based on contaminant concentrations reported by the Electric Power Research Institute.
Findley, John↗
Confinement Effects on Proton Transfer in TiO 2 Nanopores from Machine Learning Potential Molecular Dynamics Simulations
Improved understanding of proton transfer in nanopores is critical for a wide range of emerging applications, yet experimentally probing mechanisms and energetics of this process remains a significant challenge. To help reveal details of this process, we developed and applied a machine learning potential derived from first-principles calculations to examine water reactivity and proton transfer in TiO 2 slit-pores. Here, we find that confinement of water within pores smaller than 0.5 nm imposes strong and complex effects on water reactivity and proton transfer. Although the proton transfer mechanism is similar to that at a TiO 2 interface with bulk water, confinement reduces the activation energy of this process, leading to more frequent proton transfer events. This enhanced proton transfer stems from the contraction of oxygen–oxygen distances dictated by the interplay between confinement and hydrophilic interactions. Our simulations also highlight the importance of the surface topology, where faster proton transport is found in the direction where a unique arrangement of surface oxygens enables the formation of an ordered water chain. In a broader context, our study demonstrates that proton transfer in hydrophilic nanopores can be enhanced by controlling pore size, surface chemistry, and topology.
36 MATERIALS SCIENCE↗
RASPA3
RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.
Dubbeldam, David↗
Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations
Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale molecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.
X-ray imaging↗
Investigation of the iodate sorption mechanism by CoAl LDH through experiments and ab initio molecular dynamics simulations
Radioiodine released during the nuclear-fuel cycle constitutes a persistent radiological hazard. In this study, the IO 3 – uptake mechanism of CoAl LDH was resolved by combining pH-controlled sorption experiments, synchrotron XAFS, and DFT-based AIMD simulations. At pH close to 6, approximately 90% of IO 3 – was removed, and the equilibrium distribution coefficient reached about 1.7 × 10 4 mL g –1 . EXAFS analysis indicated an average iodine–oxygen bond length of 1.81 Å and a coordination number near 3, with the fit R-factor equal to 0.002. The simulations faithfully reproduced the experimental spectrum and revealed transient proton hopping events that generated metastable I–O–H species inside the interlayer, thereby confirming nitrate-to-iodate exchange as the controlling capture pathway. Atomic density profiles and radial distribution functions further showed that IO 3 – adopt an end-on orientation perpendicular to the hydroxide sheets, while water molecules mediate proton migration without disturbing the host lattice. In conclusion, the integrated experimental–computational evidence demonstrates that CoAl LDH can rapidly and selectively sequester IO 3 – under near-neutral conditions, offering atomic scale guidance for the rational engineering of layered sorbents for advanced radioactive-waste treatment.
Kang, Jaehyuk [Jeju National Univ. (Korea, Republi↗
Investigation of thermomechanical behaviors of acrylate-based shape memory polymers by molecular dynamics simulation
Not Available
Xi, Jiaxin (ORCID:0000000171363114)↗
Overcoming the Chemical Complexity Bottleneck in on-the-Fly Machine Learned Molecular Dynamics Simulations
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗