SEARCH · Engineering Papers

Results for “Bulk behavior”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Accelerated Discovery of Solar Thermochemical Hydrogen Production Materials via High-Throughput Computational and Experimental Methods

In this project, combinatorial synthesis and testing methods were combined with high-throughput materials theory calculations to greatly accelerate the discovery of thermodynamically suitable candidates for green hydrogen production via a two-stage solar thermochemical water splitting (STCH) process. Over the course of the project, more than 8000 quinary and higher oxide compositions were computationally screened for STCH viability, and detailed stability calculations were performed for more than 30 of the most promising identified compositional archetypes. As a result, three new STCH capable compositional families were discovered and experimentally verified. The first, Ce x Sr 2-x MnO 4 (CSM), represents the first known Ruddlesden-Popper compound to show STCH activity, and thus demonstrates that perovskite-related structures may hold promise for this application. The second family, Sr 1-x Ce x MnO 3 (SCM), is the simple perovskite sister-analog to CSM. Sr 0.7 Ce 0.3 MnO 3 (SCM30), a member of this compositional family, was found to produce the highest hydrogen yields of any compound tested in this project, exceeding the end of project milestone target of > 150 μmol H 2 /gram oxide at a reduction temperature of 1350 °C, although only at steam-to-hydrogen ratios greater than 1000:1. Finally, we proved that a third novel Sr-and Mn-containing family, Sr 1-x Ca x Ti 1-y Mn y O 3 (SCTM), which was identified by Materials Project tools, also splits water. The behavior of the SCTM system was found to be similar to the previously discovered Sr 1-x La x Al 1-y Mn y O 3 (SLMA) family, albeit with lower H 2 yields. Across the three thrusts of the project (computational, combinatorial, and bulk testing), five journal articles were published. As part of Program End Analysis and Data Dissemination, relevant data used for the publications was uploaded to the HydroGEN Data Hub for public access, and in certain cases, results were added to public materials databases.

08 HYDROGEN↗

NASA NTRS1989

Neutron and high speed photogrammetric arcjet diagnosis

Two methods for real time internal diagnostics of arcjet engines are described. One method uses cold, thermal, or epithermal neutrons. Cold neutrons are used to detect the presence and location of hydrogenous propellants. Thermal neutrons are used to delineate the edge contours of anode and cathode surfaces and to measure stress/strain. Epithermal neutrons are used to measure temperatures on arcjet surfaces, bulk material temperatures, and point temperatures in bulk materials. It is found that this method, with an exposure time of 10 min, produces at temperature accuracy for W or Re of + or - 2.5 C. The other method uses visible-light high-speed photogrammetry to obtain images of the transient behavior of the arc during start-up and to relate this behavior to electrial supply characteristics such as voltage, current, and ripple.

Stewart, P. A. E.↗

DOE OSTI2025

Unravelling the radiation-induced redox chemistry of plutonium ions in aqueous solution

Plutonium plays a critical role in nuclear fuel cycle technologies, but our understanding of its fundamental radiation-induced redox chemistry is limited. Changes in oxidation states affect the speciation and transport of plutonium ions in solution. For example, solvent extraction techniques used to separate and recover plutonium from used nuclear fuel rely on the selective formation, maintenance, and complexation of specific plutonium oxidation states. However, radiolytically generated radicals, ions, and molecules can drive the oxidation state distribution of plutonium ions far from equilibrium, ultimately changing the physical and chemical properties of the bulk system. These radiation-induced processes are inevitable due to the ionizing radiation fields generated by the radioactive decay of plutonium and its daughter nuclides. Therefore, mechanistically understanding how plutonium's various oxidation states respond to ionizing radiation is essential for predicting its behavior in solution. Here, we present significant advances in our understanding of radiation-induced plutonium redox chemistry by using time-resolved (electron pulse) and dose accumulation (alpha and gamma) irradiation techniques, along with quantitative multiscale modeling methods.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2024

Impact of grain boundary and surface diffusion on predicted fission gas bubble behavior and release in UO 2 fuel

In this work, we quantify the impact of grain boundary (GB) and surface diffusion on fission gas bubble evolution and fission gas release in UO 2 nuclear fuel using simulations with a hybrid phase field/cluster dynamics model. Here, we begin with a comprehensive literature review of uranium vacancy and xenon atom diffusivity in UO 2 through the bulk, along GBs, and along surfaces. In our model we represent fast GB and surface diffusion using a heterogeneous diffusivity that is a function of the order parameters that represent bubbles and grains. We find that the GB diffusivity directly impacts the rate of gas release via GB transport, and that the GB diffusivity is likely below 104 times the lower value from Olander and van Uffelen. We also find that the surface diffusivity impacts bubble coalescence and mobility, and that the bubble surface diffusivity is likely below 10 -4 times the value from Zhou and Olander.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Overscreening, Coion-Dominated Electroosmosis, and Electric Field Strength Mediated Flow Reversal in Polyelectrolyte Brush Functionalized Nanochannels

Controlling the direction and strength of nanofluidic electrohydrodyanmic transport in the presence of an externally applied electric field is extremely important in a number of nanotechnological applications. Here, we employ all-atom molecular dynamics simulations to discover the possibility of changing the direction of electroosmotic (EOS) liquid flows by merely changing the electric field strength in a nanochannel functionalized with polyelectrolyte (PE) brushes. In exploring this, we have uncovered three facets of nanoconfined PE brush behavior and resulting EOS transport. First, we identify the onset of an overscreening effect: such overscreening refers to the presence of more counterions (Na + ) within the brush layer than needed to neutralize the negative brush charges. Accordingly, as a consequence of the overscreening, in the bulk liquid outside the brush layer, there is a greater number of co-ions (Cl – ) than counterions in the presence of an added salt (NaCl). Second, this specific ion distribution ensures that the overall EOS flow is along the direction of motion of the co-ions. Such co-ion-dictated EOS transport directly contradicts the notion that EOS flow is always dictated by the motion of the counterions. Finally, for large-enough electric fields, the brush height reduces significantly, causing some of the excess overscreening-inducing counterions to squeeze out of the PE brush layer into the brush-free bulk. As a result, the overscreening effect disappears and the number of co-ions and counterions outside the PE brush layer become similar. Despite that there is an EOS transport, this EOS transport, unlike the standard EOS transport that occurs due to the imbalance of the co-ions and counterions, occurs since a larger residence time of the water molecules in the first solvation shell of the counterions (Na + ) ensures a water transport in the direction of motion of the counterions. In this work, the net effect is the reversal of the direction of the EOS transport by merely changing the strength of the electric field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Anomalous thickness-dependent electrical conductivity in van der Waals layered transition metal halide, Nb 3 Cl 8

Understanding the electronic transport properties of layered, van der Waals transition metal halides (TMHs) and chalcogenides is a highly active research topic today. Of particular interest is the evolution of those properties with changing thickness as the 2D limit is approached. Here, we present the electrical conductivity of exfoliated single crystals of the TMH, cluster magnet, Nb 3 Cl 8 , over a wide range of thicknesses both with and without hexagonal boron nitride (hBN) encapsulation. The conductivity is found to increase by more than three orders of magnitude when the thickness is decreased from 280 µm to 5 nm, at 300 K. At low temperatures and below ~50 nm, the conductance becomes thickness independent, implying surface conduction is dominating. Temperature dependent conductivity measurements indicate Nb 3 Cl 8 is an insulator, however, the effective activation energy decreases from a bulk value of 310 meV to 140 meV by 5 nm. X-ray photoelectron spectroscopy (XPS) shows mild surface oxidation in devices without hBN capping, however, no significant difference in transport is observed when compared to the capped devices, implying the thickness dependent transport behavior is intrinsic to the material. A conduction mechanism comprised of a higher conductivity surface channel in parallel with a lower conductivity interlayer channel is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Investigation of the High-Pressure Behaviors of Amblygonite by Single-Crystal X-ray Diffraction, Raman Spectroscopy, and DFT Calculations

In the present study, we extensively explored the high-pressure behaviors and vibrational properties of amblygonite LiAlPO 4 F with elevated pressures up to 34.3 GPa based on single-crystal X-ray diffraction measurements, Raman spectroscopy, and DFT calculations. The compressibility and elastic properties of amblygonite were determined first. Specifically, the obtained isothermal bulk modulus of LiAlPO 4 F is 128(4) GPa and the triclinic phase exhibited anisotropic compression with axial compressibility β c > β a > β b with a ratio of 1.11:1.00:1.20. The Raman spectra showed no indication of phase transformation and were used to obtained mode Grüneisen parameters. The average Grüneisen parameter for PO 4 tetrahedral sites was smaller than for the LiO 4 F sites. Our results provide new insights into the phase stability and elastic properties of lithium-fluorite granites at extreme conditions.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Bulk Formation of Metallic Glasses and Amorphous Silicon from the Melt

By using metallic glass compositions with a high relative glass transition temperature, such as Pd40Ni40P20, homogeneous nucleation also becomes negligible. Large (5g) masses of this alloys were obtained using a molten B2O3 flux. Presently, bulk glass formation in iron based glasses is being investigated. It is expected that if an undercooling of about 250K can be achieved in a Ge or Si melt, formation of the amorphous semiconductor phase (rather than the crystal) may be kinetically favored. The volumetric behavior of undercooled liquid Ga droplet dispersion is investigated by dilatometry. A theoretical model (both analytical and numerical) was developed for transient nucleation in glass forming melts. The model, originally designed for isothermal conditions, was extended to continuous quenching. It is being applied to glass formation in various metallic and oxide systems. A further refinement will be the inclusion of diffusion controlled interfacial rearrangements governing the growth of the crystal embryos.

Spaepen, F.↗

NASA NTRS1996

Nucleation and Crystallization of Globular Proteins: What we Know and What is Missing

Recently. much progress has been made in understanding the nucleation and crystallization of globular proteins, including the formation of compositional and structural crystal defects, Insight into the interactions of (screened) protein macro-ions in solution, obtained from light scattering, small angle X-ray scattering and osmotic pressure studies. can guide the search for crystallization conditions. These studies show that the nucleation of globular proteins is governed by the same principles as that of small molecules. However, failure to account for direct and indirect (hydrodynamic) protein interactions in the solutions results in unrealistic aggregation scenarios. Microscopic studies of numerous proteins reveal that crystals grow by the attachment of growth units through the same layer-spreading mechanisms as inorganic crystals. Investigations of the growth kinetics of hen-egg-white lysozyme (HEWL) reveal non-steady behavior under steady external conditions. Long-term variations in growth rates are due to changes in step-originating dislocation groups. Fluctuations on a shorter timescale reflect the non-linear dynamics of layer growth that results from the interplay between interfacial kinetics and bulk transport. Systematic gel electrophoretic analyses suggest that most HEWL crystallization studies have been performed with material containing other proteins at percent levels. Yet, sub-percent levels of protein impurities impede growth step propagation and play a role in the formation of structural/compositional inhomogeneities. In crystal growth from highly purified HEWL solutions, however, such inhomogeneities are much weaker and form only in response to unusually large changes in growth conditions. Equally important for connecting growth conditions to crystal perfection and diffraction resolution are recent advances in structural characterization through high-resolution Bragg reflection profiling and X-ray topography.

Rosenberger, F.↗

NASA NTRS2008

Evolution of Indarch (EH4 Chondrite) at 1 GPa and High Temperature

The chondritic meteorites are materials that are as old as the solar system itself characterized by variations in bulk chemical and oxidation state, and have long been considered possible building blocks that accreted to form the terrestrial inner planets. Enstatite chondrites contain nearly FeO free enstatite, silicon-rich kamacite and various sulfides indicating formation under highly reducing conditions. These materials could have participated in the formation of the Earth. However, "fingerprinting" of meteoritic materials has shown that no known meteoritic class corresponds to a hypothetical bulk Earth composition in every aspect. To derive constraints on early accretion and differentiation processes and possibly resolve the debate on the formation of the Earth, it is required to study experimentally a variety of chondritic materials and investigate their melting relations and elemental partitioning behavior at variable pressure (P), temperature (T) and oxygen fugacities (fO2). Variations in fO2 can indeed change chemical features and phase equilibria dramatically. The P-T phase diagrams of peridotites and carbonaceous chondrites have been extensively studied experimentally up to pressures and temperatures corresponding to the transition zone and lower mantle. Even though partial melting experiments have been conducted at ambient pressure on the enstatite chondrite Indarch, enstatite meteorites have never been experimentally investigated at high PT. The following investigation focuses on the effect of the fO2 on the phase relations of Indarch, an EH4 chondrite.

Berthet, S.↗

DOE OSTI2021

The scaling of kinetic and transport behaviors in the solution-phase chemistry of a plasma–liquid interface

The reactions at a plasma–liquid interface often involve species such as the solvated electron or the hydroxyl radical, which initiate the reduction or oxidation of solution-phase reactants (so-called scavengers) or are consumed by their own second-order recombination. Here, the mathematical scaling of the reaction–diffusion equations at the interface is used to obtain a characteristic time that can be used to determine the transition from highly efficient scavenger reduction or oxidation to lower efficiencies due to transport limitations. The characteristic time (tc) is validated using numerical solutions of the reaction–diffusion equations. When the scavenger kinetics are faster than second-order recombination, this characteristic transition time scales proportionally with the scavenger diffusivity (Ds) and the square of the scavenger bulk concentration (SB) and inversely proportional to the electron flux (J) squared; that is, t c = D s S B 2 F 2 /J 2 , where F is Faraday’s constant. However, when the scavenger kinetics are comparable or slower than second-order recombination, this scaling does not hold. Extending this analysis to three-dimensional systems shows that the profile of the electron flux on the surface affects the spatial location where reactions are most effective. Finally, the assessment of the implications of these behaviors for the reactor design highlights how effectively controlling the electron flux and solution transport may be necessary to improve the efficiency of scavenger reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Local Multimodal Electro-Chemical-Structural Characterization of Solid-Electrolyte Grain Boundaries

Typical models of polycrystalline ionic materials treat the grain boundary properties as single valued, without consideration of the full range of values that define the macroscopically measured average. Here a unique experimental platform suitable for local multimodal characterization of individual grain boundaries in bicrystal fibers is reported. A variation of three orders of magnitude in the grain boundary conductivity of ceria is observed, as measured across six individual bicrystals by both alternating current impedance spectroscopy and direct current (D.C.) linear sweep voltammetry. Nonlinear behavior of the D.C. measurements is consistent with resistance due to a space charge effect. Time-of-flight secondary ion beam spectroscopy reveals a correlation between grain boundary resistance and the concentration of impurities Si, Al, and Ca segregated at the grain boundaries, although the bulk concentrations of these impurities are negligible. Electron backscatter diffraction analysis of the crystal orientations suggests a correlation between the misorientation across the grain boundaries and grain boundary resistance. These correlations point towards a grain boundary resistance that arises from impurity-generated space charge effects and variations in impurity concentration and hence resistivity driven by the energetics of impurity segregation to grain boundaries of differing surface energies.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗

DOE OSTI2021

Strong-field physics in three-dimensional topological insulators

We investigate theoretically the strong-field regime of light-matter interactions in the topological-insulator class of quantum materials. In particular, we focus on the process of nonperturbative high-order harmonic generation from the paradigmatic three-dimensional topological insulator bismuth selenide (Bi 2 Se 3 ) subjected to intense midinfrared laser fields. We analyze the contributions from the spin-orbit-coupled bulk states and the topological surface bands separately and reveal a major difference in how their harmonic yields depend on the ellipticity of the laser field. Furthermore, bulk harmonics show a monotonic decrease in their yield as the ellipticity increases, in a manner reminiscent of high harmonic generation in gaseous media. However, the surface contribution exhibits a highly nontrivial dependence, culminating with a maximum for circularly polarized fields. We attribute the observed anomalous behavior to (i) the enhanced amplitude and the circular pattern of the interband dipole and the Berry connections in the vicinity of the Dirac point and (ii) the influence of the higher-order, hexagonal warping terms in the Hamiltonian, which are responsible for the hexagonal deformation of the energy surface at higher momenta. The latter are associated directly with spin-orbit-coupling parameters. Our results thus establish the sensitivity of strong-field-driven high harmonic emission to the topology of the band structure as well as to the manifestations of spin-orbit interaction.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Phonon transport governed by intrinsic scattering in short-period AlN/GaN superlattices

Here we employ density functional theory based phonon transport methods to provide a rigorous understanding of the nature of thermal transport in coherent short-period AlN/GaN superlattices (SLs), with period lengths up to three unit cells of each, and compare these with properties of their bulk constituents. Increasing the period length leads to phonon band folding with frequency gaps and thus smaller phonon velocities and more phonon scattering than in bulk, both of which reduce lattice thermal conductivity (κ). Contrary to expectations, we find that velocity variations among larger-period AlN/GaN SLs play only a minor role in cross-plane κ reductions, while variations in intrinsic phonon scattering are strongly correlated with their transport behaviors. This work provides insights into the microscopic behaviors of technologically relevant nanostructured materials, which are likely applicable to a wider range of SL systems and other nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Vortex Matter in Confined Superconductors and Mesoscopic Hybrid Heterostructures. Final report

This program studied the fundamental features that underlie the behavior of hybrid ferromagnet/superconductor heterostructures both in the regime of magnetic coupling and proximity effect. In ferromagnet/superconductor hybrid structures, the source of superconducting correlations and ferromagnetism are spatially separated. Therefore, the interplay between ferromagnetism and superconductivity results in physical effects that are not easily observed in bulk superconductors and allows to tune the local superconducting and vortex properties, realizing, thus, a new class of materials. In the case of proximity effect the ferromagnet and the superconductor are in contact and leakage of superconducting pairs is allowed into the ferromagnetic layer. Typically, the decay length for the superconducting pair correlations in the ferromagnet is of the order of one nanometer due to the typically large ferromagnetic exchange field. However, a suitable magnetic structure in the ferromagnetic material can be used to align the spins of the superconducting pairs allowing much longer penetration of superconductivity in the magnetic material. These superconducting pairs with aligned spins can carry spin and provide the basis for superconducting spintronics devices. In the case of magnetic coupling the effect of a spatial inhomogeneous stray field produced by the ferromagnet substantially modifies the nucleation of the order parameter in the superconducting layer and can lead to a significant change in the local superconducting properties. Therefore, magnetically coupled planar ferromagnet-superconductor hybrid structures offer new avenues for manipulation of the superconductivity at the nanoscale and convenient means to control vortex dynamics.

36 MATERIALS SCIENCE↗

NASA NTRS2004

JWST ISIM Distortion Analysis Challenge

Very tight distortion requirements are imposed on the JWST's ISM structure due to the sensitivity of the telescope's mirror segment and science instrument positioning. The ISIM structure is a three dimensional truss with asymmetric gusseting and metal fittings. One of the primary challenges for ISIM's analysis team is predicting the thermal distortion of the structure both from the bulk cooldown from ambient to cryo, and the smaller temperature changes within the cryogenic operating environment. As a first cut to estimate thermal distortions, a finite element model of bar elements was created. Elements representing joint areas and metal fittings use effective properties that match the behavior of the stack-up of the composite tube, gusset and adhesive under mechanical and thermal loads. These properties were derived by matching tip deflections of a solid model simplified T-joint. Because of the structure s asymmetric gusseting, this effective property model is a first attempt at predicting rotations that cannot be captured with a smeared CTE approach. In addition to the finite element analysis, several first order calculations have been performed to gauge the feasibility of the material design. Because of the stringent thermal distortion requirements at cryogenic temperatures, a composite tube material with near zero or negative CTE is required. A preliminary hand analysis of the contribution of the various components along the distortion path between FGS and the other instruments, neglecting second order effects were examined. A plot of bounding tube longitudinal and transverse CTEs for thermal stability requirements was generated to help determine the feasibility of meeting these requirements. This analysis is a work in progress en route to a large degree of freedom hi-fidelity FEA model for distortion analysis. Methods of model reduction, such as superelements, are currently being investigated.

Cifie, Emmanuel↗

DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗