Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic models”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Enhanced liquid metal wetting on oxide surfaces via patterned particles

Liquid metals typically do not intrinsically wet oxide surfaces due to chemical dissimilarity, preventing the formation of well-bonded interfaces in many applications. This paper investigates how the wetting and spreading of liquid silver (Ag) droplets on oxide surfaces can be enhanced via patterned nickel (Ni) particles using multiscale modeling. First, a force field (FF) for an Ag/yttria-stabilized-zirconia system was developed by fitting the interfacial binding energy to density-functional theory (DFT) calculations. Then, molecular dynamics (MD) simulations using these FF parameters revealed that the good wetting of Ag on Ni caused Ag droplets to be attracted to the Ni particles, promoting spreading. Thus, the optimal hexagonal Ni pattern was shown to enhance the wetting area by 224% compared to a bare oxide surface. Second, an analytical model predicting the maximum possible enhancement in wetting area as a function of the intrinsic wetting angle, size of the Ag droplet, and size of the Ni particles was derived and validated with MD simulations, to generalize this method to other substrates. Here, both the analytical model and MD results suggested that the ability of Ni patterns to enhance the Ag wetting area becomes more significant as the wettability of substrate becomes worse.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reversing inverse Hall-Petch and direct computation of Hall-Petch coefficients

A 2D bicrystal atomistic model of dislocation transmission through a Σ11<101>{131} symmetric-tilt grain boundary reveals details of Hall-Petch breakdown in the single dislocation regime in Al, Ni, and Cu. Partly based on a previous study, this research determines the stress required for a single dislocation to be transmitted through the Σ11 boundary and finds that single dislocation transmission typically deviates from Hall-Petch behavior because the leading partial dislocation becomes trapped in the boundary and emits glissile grain boundary disconnections that cause boundary sliding deformation at stresses well below boundary transmission stresses. Thus, mechanisms of inverse Hall-Petch, namely grain boundary shear and sliding, are shown to operate in lieu of grain boundary transmission in the single dislocation regime. However, by controlling the applied stresses, inverse Hall-Petch can be reversed and typical Hall- Petch grain boundary transmission regained. This, in turn, allows the direct computation of Hall-Petch coefficients. Further, a second dislocation on the same slip plane preserves typical Hall-Petch behavior for Al and Ni and does not lead to boundary sliding events, thus implying a critical grain size for Hall-Petch breakdown based on the pileup model.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Permeation of niobium through grain boundaries in copper

Low mutual solubility is expected to avert interdiffusion in layered composites of phase-separating metals. However, we show that Nb diffuses through polycrystalline Cu—despite the minimal bulk solubility of these elements—due to short circuit transport of Nb along certain grain boundaries in Cu. Atomistic modeling demonstrates that Nb-permeable Cu grain boundaries exhibit negative enthalpy of mixing of Nb, resulting in an enthalpically-stabilized solution of highly mobile Nb atoms that easily diffuse through the boundary. By contrast, Mo, which also has minimal solubility with Cu, has positive enthalpy of mixing at Nb-permeable Cu grain boundaries and does not diffuse through them. Furthermore, our findings suggest material selection and grain boundary engineering as pathways for designing improved diffusion barriers and thermally stable laminate composites.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Permeation of niobium through grain boundaries in copper

Low mutual solubility is expected to avert interdiffusion in layered composites of phase-separating metals. However, we show that Nb diffuses through polycrystalline Cu—despite the minimal bulk solubility of these elements—due to short circuit transport of Nb along certain grain boundaries in Cu. Atomistic modeling demonstrates that Nb-permeable Cu grain boundaries exhibit negative enthalpy of mixing of Nb, resulting in an enthalpically-stabilized solution of highly mobile Nb atoms that easily diffuse through the boundary. By contrast, Mo, which also has minimal solubility with Cu, has positive enthalpy of mixing at Nb-permeable Cu grain boundaries and does not diffuse through them. Furthermore, our findings suggest material selection and grain boundary engineering as pathways for designing improved diffusion barriers and thermally stable laminate composites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lpp positions peptidoglycan at the AcrA-TolC interface in the AcrAB-TolC multidrug efflux pump

The multidrug efflux pumps of Gram-negative bacteria are a class of complexes that span the periplasm, coupling both the inner and outer membranes to expel toxic molecules. The best-characterized example of these tripartite pumps is the AcrAB-TolC complex of Escherichia coli. However, how the complex interacts with the peptidoglycan (PG) cell wall, which is anchored to the outer membrane (OM) by Braun’s lipoprotein (Lpp), is still largely unknown. In this work, we present molecular dynamics simulations of a complete, atomistic model of the AcrAB-TolC complex with the inner membrane, OM, and PG layers all present. We find that the PG localizes to the junction of AcrA and TolC, in agreement with recent cryo-tomography data. Free-energy calculations reveal that the positioning of PG is determined by the length and conformation of multiple Lpp copies anchoring it to the OM. The distance between the PG and OM measured in cryo-electron microscopy images of wild-type E. coli also agrees with the simulation-derived spacing. Finally, sequence analysis of AcrA suggests a conserved role for interactions with PG in the assembly and stabilization of efflux pumps, one that may extend to other trans-envelope complexes as well.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

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

Local structure and Ca/Si ratio in C-S-H gels from hydration of blends of tricalcium silicate and silica fume

At the microscale, C-S-H gels from alite, or neat Portland cements, has a Ca/Si ratio close to 1.80. At the nanoscale, C-S-H is described by a defective tobermorite structure which allows a maximum Ca/Si ratio close to 1.40. There is no agreement in the location of the extra 0.40 mol of Ca(OH){sub 2} at the nanoscale. Atomistic modelling studies reported Ca(OH){sub 2} species within the tobermorite interlayer space. Other works point toward a fine intermixing of defective tobermorite and nanoportlandite. Here, we have prepared a series of alite blended with silica fume and studied the pastes by several techniques including synchrotron pair distribution function (PDF). In the employed conditions, the C-S-H gel formed by the pozzolanic reaction has nearly the same local structure than the primary C-S-H gel. Furthermore, differential PDF points toward Ca(OH){sub 2} excess having a local structure compatible with few-layer thick nanoparticles stretched along the c-axis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The intrinsic mechanical properties of hydromagnesite, Mg 5 (CO 3 ) 4 (OH) 2 ·4H 2 O, a key phase of reactive MgO carbonate cement

To potentially enable CO 2 sequestration, reactive MgO carbonate cement is emerging as an alternative binder to Portland cement. Understanding the mechanical properties of its binding phase is critical for understanding the strength development and performing materials design for reactive MgO cement systems; however, the intrinsic mechanical properties of hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ·4H 2 O), a key binding phase, remain unexplored. Here the present study utilized synchrotron-based high-pressure X-ray diffraction to determine the unit cell-scale, intrinsic mechanical properties of hydromagnesite for the first time. Up to hydrostatic loading of 7.7 GPa, the bulk modulus of hydromagnesite was determined as 59 GPa or 71 GPa fitted using the second-order or third-order Birch-Murnaghan equation of state, which we contextualize with binding phases in various cement systems. The experiment results are applicable in materials design of low-carbon concrete and valuable for the validation and calibration of atomistic models.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring structural transitions at grain boundaries in Nb using a generalized embedded atom interatomic potential

The advancement in experimental techniques, like the atom probe tomography and high resolution electron microscopy, is fueling interest in studying structural transformations of grain boundaries in metal and alloys to uncover correlations between mechanical properties and solute or impurity segregation to grain boundaries. Atomistic modeling is an important tool that can pinpoint the intricate dynamics of grain boundary phase transitions, but the lack of accurate interatomic potentials needed to simulate the complex dynamics of grain boundary structural transitions and identify different metastable phases has been the bottleneck. To this end, we use niobium as a model body centered cubic (BCC) metal and develop an interatomic potential to study grain boundary phase transitions. The potential for Nb is based on a generalization of the embedded atomic method potential and has sufficient flexibility to learn complex energy landscapes using a small set of training structures. We systematically test and validate the using data from ab initio density functional theory calculations and experiments. Using this potential, we calculate energies of multiple symmetric-tilt grain boundaries spanning a wide range of misorientation angles. Additionally, we explore different metastable structures of the Σ 27(552) [$1\overline{1}0$] grain boundary and use molecular dynamic simulations to study the coexistence of metastable phases and grain boundary transitions at finite temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sluggish, chemical bias and percolation phenomena in atomic transport by vacancy and interstitial diffusion in Ni-Fe alloys

In this paper, we report the results of extended atomistic modeling of intrinsic mobility of point defects and associated atomic transport in Ni-Fe model binary alloys. We consider the effects of composition and temperature and present evidence of the sluggish and chemically biased diffusion, and percolation effects occurring in atomic transport via the vacancy and interstitial migration mechanisms. The results are analyzed and discussed in the light of previous studies and some experimental observations. It is demonstrated that the sluggish diffusion, the chemically biased diffusion, and the percolation are interlinked phenomena that are defined by the chemical complexity of particular alloys. Methods for predicting these phenomena in multicomponent alloys are discussed.We report a fundamental understanding of sluggish diffusion, chemically-biased diffusion, as well as percolation phenomena, in Ni-Fe random alloys for vacancy and interstitial atom migration mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Systematic discovery of new nano-scale metastable intermetallic eutectic phases in laser rapid solidified Aluminum-Germanium alloy

Laser surface remelting of as-cast Al-Ge eutectic alloy is shown to produce ultrafine lamellar eutectic morphology with interlamellar spacing refined up to ∼60 nm and composed of FCC Al solid solution and unusual Al x Ge y intermetallic phases that do not form during near-equilibrium solidification. The microstructures are characterized and analyzed using a combination of selected area electron diffraction, high-resolution scanning transmission electron microscopy, energy dispersive X-ray spectroscopy to obtain high-resolution elemental maps, and atomistic modeling using density functional theory followed by atomic-scale image simulation. Depending on the local solidification conditions, the crystallography of the Al x Ge y intermetallic phases in the eutectic microstructure is either monoclinic (C 2/c) or monoclinic (P 2 1 ), with high densities of defects in both cases. This is in sharp contrast to the as-cast alloys that showed nominally pure Al and Ge phases with significant solute partitioning and equilibrium FCC and diamond cubic crystal structures, respectively. Corresponding kinetic phase diagrams are proposed to interpret the evolution of nano-lamellar eutectic morphologies with equilibrium Al and metastable Al x Ge y phases, and to explain increased solid solubility in the Al phases manifested by precipitation of ultrafine clusters of Ge. Furthermore, the reasons for the formation of these metastable eutectics under laser rapid solidification are discussed from the perspective of the competitive growth criterion.

Al-Ge eutectic↗
DOE OSTI2020

Impact of acid site speciation and spatial gradients on zeolite catalysis

This mini-review provides an overview of the current state of acid site control in zeolite catalysts, including methods of synthesis, advanced characterization, and measured effects of acid properties (speciation, concentration, proximity, siting, and spatial distribution) on a variety of commerciallyrelevant reactions. The diversity of aluminum species is described with respect to their location at specific sites in zeolite crystals as well as mesoscopic gradients in elemental composition that give rise to zoned or core-shell architectures. Challenges in the identification of acid siting are highlighted within the context of trial-and-error synthesis methods for a range of aluminosilicate frameworks, which hinder a priori design of catalysts, and limitations in techniques to suitably characterize active sites. Emphasis is also placed on knowledge gaps in zeolite catalysis wherein broad development of structure-performance relationships relies on future advancement of synthesis and analytical methods in parallel with atomistic modeling.

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

A comparative study of the effectiveness of empirical potentials for molecular dynamics simulations of borosilicate glasses

Despite their practical importance, atomistic modeling of B 2 O 3 containing glasses have been challenging due to the lack of reliable empirical potentials. Fortunately, a few recent developments have shown promises to simulate these glasses where the boron coordination has complex and non-linear dependence on glass composition. This work aims to provide an evaluation of the effectiveness of three recently developed potentials by a systematic study of a series (~20) of sodium borosilicate glasses with constant K (ratio of [SiO 2 ]/[B 2 O 3 ]) =2 and varying R (ratio of [Na 2 O]/[B 2 O 3 ] ranging from 0.1 to 4) values and several sodium boroaluminosilicate glasses. A comparative assessment was established on the basis of the short- and medium-range structure features, such as boron N 4 values, total correlation functions, bond angle distribution, oxygen speciation, and mechanical properties using experimental or well-established models as criteria. Furthermore, this work provides insights on the choices of empirical potentials for MD simulations of borosilicate glasses and gives directions of future potential development and refining.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Identification and correction of temporal and spatial distortions in scanning transmission electron microscopy

Scanning transmission electron microscopy (STEM) has become the technique of choice for quantitative characterization of atomic structure of materials, where the minute displacements of atomic columns from high-symmetry positions can be used to map strain, polarization, octahedra tilts, and other physical and chemical order parameter fields. The latter can be used as inputs into mesoscopic and atomistic models, providing insight into the correlative relationships and generative physics of materials on the atomic level. However, these quantitative applications of STEM necessitate understanding the microscope induced image distortions and developing the pathways to compensate them both as part of a rapid calibration procedure for in situ imaging, and the post-experimental data analysis stage. Here, we explore the spatiotemporal structure of the microscopic distortions in STEM using multivariate analysis of the atomic trajectories in the image stacks. Based on the behavior of principal component analysis (PCA), we develop the Gaussian process (GP)-based regression method for quantification of the distortion function. The limitations of such an approach and possible strategies for implementation as a part of in-line data acquisition in STEM are discussed. Here, the analysis workflow is summarized in a Jupyter notebook that can be used to retrace the analysis and analyze the reader's data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Broad Distribution of Local Polar States Generates Large Electrothermal Properties in Pb-Free Relaxor Ferroelectrics

Electrothermal energy conversion provides attractive solutions for global energy management, such as energy harvesting from waste heat using pyroelectric energy conversion (PEC) and efficient cooling of portable electronics or data servers using the electrocaloric effect. Relaxor ferroelectrics are attractive for electrothermal energy conversion because of their large pyroelectric coefficients over a wide temperature range. Although Pb-based relaxors are well-known, toxicity concerns have mandated the intense search for Pb-free alternatives. Here, we engineered (Ba,Ca)TiO 3 -based relaxors based on a multisite doping strategy, which show promising electrothermal performance, viz. a maximum PEC efficiency of 14% and electrocaloric refrigeration capacity of 115 J/kg. Using local-scale structural analysis, we provide an atomistic model for large electrothermal properties in the newly designed Pb-free ferroelectrics, whereby a temperature-independent continuous distribution of cation displacement directions creates easy pathways for microscopic polarization reorientation. This research provides key structural insight for future atomic-scale engineering of environmentally sustainable ferroelectrics in energy applications.

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

Simulation of Inelastic Neutron Scattering Spectra Directly from Molecular Dynamics Trajectories

Inelastic neutron scattering (INS) is a widely used technique to study atomic and molecular vibrations. With the increasing complexity of materials and thus the INS spectra, being able to simulate the spectra from various atomistic models becomes an essential step and also a major bottleneck for INS data analysis. The conventional approach using density functional theory and lattice dynamics often falls short when the materials of interest are complex (e.g., defective, disordered, heterogeneous, amorphous, large-scale), for which molecular dynamics driven by an interatomic force field is a more common approach. In this paper, we demonstrate a method to directly convert molecular dynamics trajectories into simulated INS spectra, including not only fundamental but also higher order excitations. In this work, the results are compared with data collected on various representative samples from different neutron spectrometers. This development will open great opportunities by providing the key tool to perform in-depth analysis of INS data and to validate and optimize computer models.

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

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗