Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “single atom alloys”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Synthesis and Transport Properties of the Family of Zintl Phases Ca 3 RESb 3 (RE = La–Nd, Sm, Gd–Tm, Lu): Exploring the Roles of Crystallographic Disorder and Core 4f Electrons for Enhancing Thermoelectric Performance

Zintl phases with complex crystal structures have been studied as promising candidate-materials for thermoelectric (TE) applications. Here, we report the syntheses of the family of rare-earth metal Zintl phases with the general formula Ca 4–x RE x Sb 3 (x ≈ 1; RE = La–Nd, Sm, Gd–Tm, Lu). The structural elucidation is based on refinements of single-crystal X-ray diffraction data for 12 unique chemical compositions. The cubic structure is confirmed as belonging to the anti-Th 3 P 4 structure type (space group I4¯3d, no. 220, Z = 4), where the Ca and RE atoms share the same atomic site with ca. 75% and 25% occupancies, respectively. Such crystallographic disordering of divalent Ca and trivalent RE atoms in the structure provides a pathway to intricate bonding. The latter, together with the presence of heavy elements such as Sb and the lanthanides, are expected to enhance the scattering probability of phonons, thereby leading to as low thermal conductivity κ as that of the ordered RE 4 Sb 3 . The drive of the hypothetical parent compound Ca 4 Sb 3 to be stabilized by alloying with rare-earth metals can be understood following the Zintl-Klemm concept, as the resultant formula may be rationalized as (Ca 2+ ) 3 RE 3+ (Sb 3– ) 3 , indicating the realization of closed-shell electronic configurations for all elements. This notion is confirmed by electronic structure calculations, which reveal narrow bandgaps E g = 0.77 and 0.53 eV for Ca 3 LaSb 3 and Ca 3 LuSb 3 , respectively. Additionally, the incorporation of RE atoms into the structure drives the phase into a state of a degenerate semiconductor with dominant hole charge carriers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tetragonal structure and uniaxial magnetic anisotropy in the arc-melted (Ce,Zr) 2 (Fe,M) 17 (M = Mo, W, Co) alloys

The tetragonal ThMn 12 -type structure is stabilized in (Ce 1-x Zr x ) 2 T 16 M (x = 0.2 – 0.3; T = Fe or Fe/Co, M = Mo or W) arc-melted alloys. Approximately 5 at% of Mo or W admixture is sufficient to transform the hexagonal Th 2 Ni 17 -type structure of (Ce 1-x Zr x ) 2 Fe 17 into the tetragonal ThMn 12 -type structure of nearly single-phase (Ce 1-x Zr x ) 2 Fe 16 M and/or (Ce 1-x Zr x ) 2 Fe 15 CoM bulk alloys. X-ray Rietveld refinements reveal that Zr and Mo (W) substitute different sites in the tetragonal crystal structure. Zirconium preferentially replaces Ce, whereas Mo (W) substitutes Fe. At room temperature, the tetragonal phases exhibit strong ferromagnetism and a uniaxial magneto-crystalline anisotropy with anisotropy fields of 15 – 17 kOe. These materials possess room-temperature saturation magnetizations of 95 – 105 emu/g and Curie temperature of 420 – 505 K. At some concentrations, concurrent Zr and Mo (W) site occupancy facilitates 1:12 structure formation in bulk alloys with minimal presence of non-magnetic element in the Fe sublattice, thus securing the highest known magnetic moment per Fe atom (∼ 1.37 µ B ) in this type of materials. The intrinsic magnetic characteristics, as well as the absence of critical rare earths (RE), make these compounds interesting for development as low-cost permanent magnets. In conclusion, small Co additions improve the Curie temperature, especially if combined with W.

36 MATERIALS SCIENCE↗
DOE OSTI2023

TaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-vanadium (Tanu-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the element V. For Ta, p semi-core states as valence were chosen. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name case-N where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. and Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15-50 (1996) (4) Kresse, G. and Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-tantalum (Nb-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb. For Ta, p semi-core states as valence were chosen. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name case-N where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. and Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15-50 (1996) (4) Kresse, G. and Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

Non-Equilibrium Synthesis of Semiconductor Alloys for Energy Applications

This project explores the synthesis under equilibrium/non-equilibrium conditions of semiconductor alloys critical to SSL to impose top-down influence on the formation of atomic-scale microstructures. A plasma-enhanced MOCVD (PE-MOCVD) technique is employed, in which plasma generated within the low-pressure reactor environment introduces additional energy that is delivered to the gas-phase reactants and growth surface, allowing relatively low substrate temperatures to be maintained, rapidly switching the crystal-growth environment into a non-equilibrium state. The project coordinates on-going work at National Renewable Energy Laboratory (NREL) that addresses order/disorder Al x In 1-x P (AlInP) confinement structures, which are predicted to surpass the power-conversion efficiency of current materials used for red/amber LEDs as components of so-called color-mixed LEDs, which integrate multiple LEDs in a single device for improved color rendering. Ordering, phase separation, and strain relaxation are central to this technology. During growth, each spontaneous process could be suppressed or enhanced on demand by application of plasma as an additional handle. This technique could have broader applications to other energy technologies, such as photovoltaics.

14 SOLAR ENERGY↗
DOE OSTI2024

Epitaxial (Al x Ga 1-x-y In y ) 2 O 3 Alloys Lattice Matched to Monoclinic Ga 2 O 3 Substrates

We have epitaxially stabilized a series of monoclinic (Al x Ga 1-x-y In y ) 2 O 3 alloys by careful choice of molecular beam epitaxy growth conditions, which balance alloy growth with suboxide desorption. The films are pseudomorphic to (010) β-Ga 2 O 3 substrates at thicknesses up to 150 nm with compositions ranging from (Al 0.01 Ga 0.83 In 0.16 ) 2 O 3 to (Al 0.24 Ga 0.75 In 0.03 ) 2 O 3 . The absorption edge shifts from approximately 4.62-5.14 eV with coincidently increasing Al and decreasing In mole fractions. J-V measurements reveal an increase in resistivity over four orders of magnitude with a maximum value of 4.2 x 10 5 Ω-cm for (Al 0.17 Ga 0.76 In 0.07 ) 2 O 3 , which has nearly identical lattice parameters (both in-plane and out-of-plane) to the underlying β-Ga 2 O 3 . Scanning transmission electron microscopy of this sample reveals a mostly uniform and single crystalline film, though we identify areas of non-uniform In incorporation and some γ-phase inclusions. This work demonstrates the feasibility of thick layers lattice-matched to β-Ga 2 O 3 with increased bandgap compared to phase-separation limited (Al,Ga) 2 O 3 . These alloys can enable higher bandgap epitaxial dielectrics and high sheet charge density transistors by increasing the conduction band offset with respect to β-Ga 2 O 3 .

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

CO 2 Adsorption and Hydrogenation on Inverse InO x /Cu(111) Catalysts: Active Role of the Oxide–Metal Interface

The direct conversion of carbon dioxide (CO 2 ) into methanol via hydrogenation is essential for industrial applications. Recent studies on catalysts that contain an inverse oxide/metal configuration have shown very good catalytic performance for the CO 2 hydrogenation to methanol process. Here, in this study, we investigated the behavior of indium oxide-Cu(111) interfaces under pure CO 2 and CO 2 /H 2 mixtures using synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS). Initially, a single layer of copper oxide (Cu x O) was grown on the Cu(111) surface by controlled oxidation. On this surface, indium was deposited at room temperature. Oxygen atoms transferred from Cu x O/Cu(111) to the indium metal upon deposition, forming In-O-Cu bonds and active interfaces. Although Cu(111) is not very active for the binding and activation of CO 2 , the formed InO x -Cu(111) interfaces had no problem adsorbing and dissociating the molecule at room temperature. Reaction of CO 2 with H 2 on InO x -Cu(111) yielded surface-bound H 3 CO, CO 2 δ− , CO 3 , and CH x species that are typical intermediates in the production of methanol and other oxygenates. The InO x -Cu(111) interface underwent dynamic chemical changes under reaction conditions, forming In-Cu alloys at low indium coverages (< 0.05 monolayer), while at higher indium coverages a mixture of an In-Cu and InO x was detected in XPS. These findings indicate that InO x /In-Cu interfaces can play a key role in processes aimed at the trapping and valorization of CO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase instabilities in austenitic steels during particle bombardment at high and low dose rates

Disruption of phase stability by energetic particle bombardment is a major challenge in designing advanced radiation-tolerant alloys and ion beam processing of nanocomposites. Particularly, ballistic dissolution susceptibility of different solute nanocluster species in alloys is poorly understood. Here, low dose rate neutron irradiations were conducted on a Fe-Cr-Ni based austenitic steel in the BOR-60 reactor (9.4 × 10 -7 dpa/s, 318 °C) followed by accelerated dose rate ion irradiations at multiple temperatures ( 10 -3 dpa/s, 380 – 420 °C). Using atom probe tomography, the stability of radiation-enhanced Cu-rich and radiation-induced Ni-Si-Mn-rich nanoclusters was evaluated. During neutron irradiation, Cu-rich clusters nucleated with their core concentrations progressively increasing with dose, while Ni-Si-Mn-rich clusters formed and evolved into G-phase precipitates. Ion irradiations dramatically altered the nanoclusters. Cu-rich clusters were ballistically dissolved, but Ni-Si-Mn-rich clusters remained stable and coarsened with dose at 400 and 420 °C, highlighting that different nanocluster species in a single microstructure can have innately distinct ballistic dissolution susceptibilities. Solute-specific recoil rates were incorporated into the Heinig precipitate stability model, which shows that in addition to radiation-enhanced diffusion, recovery from ballistic dissolution depends on solute concentration gradient near cluster interfaces. Finally, the combined experimental-modeling study quantified the critical temperatures and damage rates where ballistic dissolution dominates for each cluster species.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗
NASA NTRS2024

Molecular Dynamics Simulations of Microstructural Effects on Austenite-Martensite Interfaces in NiTi

Formation and migration of austenite-martensite interfaces plays the key role in reversible martensitic transformations of shape memory alloys (SMAs). How these interfaces interact with the SMA microstructure is a primary determining factor in important functional properties such as hysteresis and transformation span. As such, successful microstructural engineering of SMAs requires in-depth knowledge of interface behavior. The rapid nature of martensitic transformations makes experimental observations of moving austenite-martensite interfaces challenging. Molecular dynamics (MD) simulation is a unique tool which can probe the atomic-scale details of austenite-martensite interfaces as they migrate through different microstructures. However, in focusing on the entire transformation process, including the nucleation of new phases, MD studies are usually performed so far from equilibrium that their relevance to experiment is questionable. Here, we demonstrate new MD techniques to generate energetically preferred austenite-martensite interfaces in NiTi under near-equilibrium conditions. The interfaces are semi-coherent, exhibiting a series of structural disconnections, and they can migrate rapidly through single crystals under only small thermodynamic driving forces. In contrast, when interfaces migrate in polycrystals, their motion is impeded by thermoelastic effects as well as changes in orientation relationships at grain boundaries. Microstructures which accumulate large amounts of elastic energy tend to release some fraction through irreversible, hysteresis-inducing mechanisms. We demonstrate that engineering microstructures with less constraints is a viable strategy to produce SMAs with reduced hysteresis and transformation span. Similar thermoelastic and hysteresis-inducing mechanisms also arise when austenite-martensite interfaces encounter precipitates and can be controlled by tuning characteristics of the precipitate distribution.

Gabriel Plummer↗
NASA NTRSDate not supplied

Molecular Dynamics Simulations of Microstructural Effects on Austenite-Martensite Interfaces in NiTi

The formation and migration of austenite-martensite interfaces plays the key role in reversible martensitic transformations of shape memory alloys (SMAs). How these interfaces interact with the SMA microstructure is a primary determining factor in important functional properties such as hysteresis and transformation span. As such, successful microstructural engineering of SMAs requires in-depth knowledge of interface behavior. The rapid nature of martensitic transformations makes experimental observations of moving austenite-martensite interfaces challenging. Molecular dynamics (MD) simulation is a unique tool which can probe the atomic-scale details of austenite-martensite interfaces as they migrate through different microstructures. However, in focusing on the entire transformation process, including the nucleation of new phases, MD studies are usually performed so far from equilibrium that their relevance to experiment is questionable. Here, we demonstrate new MD techniques to generate energetically preferred austenite-martensite interfaces in NiTi under near-equilibrium conditions. The interfaces are semi-coherent, exhibiting a series of structural disconnections, and they can migrate rapidly through single crystals under only small thermodynamic driving forces. In contrast, when interfaces migrate in polycrystals, their motion is impeded by thermoelastic effects as well as changes in orientation relationships at grain boundaries. Microstructures which accumulate large amounts of elastic energy tend to release some fraction through irreversible, hysteresis-inducing mechanisms. We demonstrate that engineering microstructures with less constraints is a viable strategy to produce SMAs with reduced hysteresis and transformation span. Similar thermoelastic and hysteresis-inducing mechanisms also arise when austenite-martensite interfaces encounter precipitates and can be controlled by tuning characteristics of the precipitate distribution.

Gabriel Plummer↗
NASA NTRSDate not supplied

Molecular Dynamics Simulations of Microstructural Effects on Austenite-Martensite Interfaces in NiTi

Formation and migration of austenite-martensite interfaces plays the key role in reversible martensitic transformations of shape memory alloys (SMAs). How these interfaces interact with the SMA microstructure is a primary determining factor in important functional properties such as hysteresis and transformation span. As such, successful microstructural engineering of SMAs requires in-depth knowledge of interface behavior. The rapid nature of martensitic transformations makes experimental observations of moving austenite-martensite interfaces challenging. Molecular dynamics (MD) simulation is a unique tool which can probe the atomic-scale details of austenite-martensite interfaces as they migrate through different microstructures. However, in focusing on the entire transformation process, including the nucleation of new phases, MD studies are usually performed so far from equilibrium that their relevance to experiment is questionable. Here, we demonstrate new MD techniques to generate energetically preferred austenite-martensite interfaces in NiTi under near-equilibrium conditions. The interfaces are semi-coherent, exhibiting a series of structural disconnections, and they can migrate rapidly through single crystals under only small thermodynamic driving forces. In contrast, when interfaces migrate in polycrystals, their motion is impeded by thermoelastic effects as well as changes in orientation relationships at grain boundaries. Microstructures which accumulate large amounts of elastic energy tend to release some fraction through irreversible, hysteresis-inducing mechanisms. We demonstrate that engineering microstructures with less constraints is a viable strategy to produce SMAs with reduced hysteresis and transformation span. Similar thermoelastic and hysteresis-inducing mechanisms also arise when austenite-martensite interfaces encounter precipitates and can be controlled by tuning characteristics of the precipitate distribution.

Gabriel Plummer↗
DOE OSTI2024

Creation of Multi-Principal Element Alloy NiCoCr Nanostructures via Nanosecond Laser-Induced Dewetting

The multi-principal element alloy nanoparticles (MPEA NPs), a new class of nanomaterials, present a highly rewarding opportunity to explore new or vastly different functional properties than the traditional mono/bi/multimetallic nanostructures due to their unique characteristics of atomic-level homogeneous mixing of constituent elements in the nanoconfinements. Here, in this paper, the successful creation of NiCoCr nanoparticles, a well-known MPEA system is reported, using ultrafast nanosecond laser-induced dewetting of alloy thin films. Nanoparticle formation occurs by spontaneously breaking the energetically unstable thin films in a melt state under laser-induced hydrodynamic instability and subsequently accumulating in a droplet shape via surface energy minimization. While NiCoCr alloy shows a stark contrast in physical properties compared to individual metallic constituents, i.e., Ni, Co, and Cr, yet the transient nature of the laser-driven process facilitates a homogeneous distribution of the constituents (Ni, Co, and Cr) in the nanoparticles. Using high-resolution chemical analysis and scanning nanodiffraction, the environmental stability and grain arrangement in the nanoparticles are further investigated. Thermal transport simulations reveal that the ultrashort (≈100 ns) melt-state lifetime of NiCoCr during the dewetting event helps retain the constituent elements in a single-phase solid solution with homogenous distribution and opens the pathway to create the unique MPEA nanoparticles with laser-induced dewetting process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MX precipitate behavior in an irradiated advanced Fe-9Cr steel: Self-ion irradiation effects on phase stability

In an effort to optimize Fe-9Cr reduced activation ferritic/martensitic (RAFM) steels and to inform the design and operation of fusion reactors, this work represents the first in a series of cohesive studies dedicated to the evolution of MX-TiC precipitates under accelerated single and dual ion irradiations. This study investigates CNA9, a simplified Fe-9Cr RAFM steel featuring initial MX-TiC precipitate densities of (2.3±0.3)×10²¹ m⁻³. This material was subjected to single self-ion irradiation at damage levels ranging from 1 to 100 displacements per atom (dpa) over a temperature range of 300 to 600°C, with a nominal dose rate of 7×10⁻⁴ dpa/s. Irradiation-induced coarsening was observed, as evidenced by statistically significant increases in mean diameter sizes, at 15 dpa at both 500°C and 600°C, whereas no coarsening was noted at 300°C or 400°C. Further, complete dissolution of precipitates occurred at damage levels of 50 and 100 dpa across the two temperatures tested (300°C and 500°C) while no significant changes were observed at any doses below 15 dpa at 500°C. Experimentally parameterized recoil resolution modeling suggests that the observed radiation stability of MX-TiC precipitates is intricately linked to diffusional changes of solutes resulting from the co-evolution of microstructural features within the experiments. The findings align with current theoretical perspectives on radiation-induced precipitate stability in complex alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Larger than uncorrelated vacancy diffusion contributions in chemically disordered crystalline materials

Arising from a variational approach to compute diffusion coefficients, we introduce “superkinetic kinosons”, which are contributions to the diffusion flux of mobile defects and atomic species by single jump mechanisms that can exceed their uncorrelated values. In vacancy-mediated diffusion in crystalline random and high entropy materials, these contributions primarily arise from intermittent jumps of slow moving atoms that remove the vacancy out of correlation traps and enable a quantification of such effects on the diffusion coefficients. They can be significant contributors to diffusion, even when faster moving atoms have infinite percolation networks. Furthermore, their existence and importance provide fundamental insights about underlying aspects of vacancy diffusion such as formation of localized correlation traps and deviation of vacancy diffusion from percolation behavior.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Constraints on the Nature and Distribution of Iridium Host Phases at the Cretaceous-Tertiary Boundary: Implications for Projectile Identity and dispersal on impact

Among Cretaceous-Tertiary boundary sites worldwide, variations in the concentrations and ratios of elements commonly enriched in meteorites complicate traditional geochemical attempts at impactor identification. Yet they may provide constraints on the physical and chemical processes associated with large-body disruption and dispersal, as well as with diagenesis of projectile components. To this end, we continue our efforts to identify the mineral host-phases of projectile-derived elements, particularly for Ir, and to document their partitioning between crater deposits and ejecta resulting from the Chicxulub basin-forming impact. Building on earlier work, we used INAA to measure Ir concentrations in successively smaller splits of finely powdered impact melt breccia from the Chicxulub Crater in Mexico (sample Y6Nl9-R(b)), and K/T boundary fish clay from Stevns Klint, Denmark (sample FC-1, split from 40 kg of homogenized material intended as an analytical standard). Results for the Chicxulub sample show a heterogeneous Ir distribution and document that at least five discrete Ir-bearing host phases were isolated in subsequent splits, having Ir masses equivalent to pure Ir spheres from about 0.8 to about 3.5 mm in diameter. Three of these are within a sufficiently reduced mass of powder to warrant searching for them using backscattered electron microscopy. In contrast, successively smaller splits of the Stevns Klint fish clay show no statistically significant deviation from the reported value of 32 +/- 2 ng/g Ir, suggesting a uniform Ir host-phase distribution. For the smallest split obtained thus far (100 +/- 40 ng/g Ir), a pure Ir sphere of equivalent Ir mass would be <0.05 min in diameter. (n.b. Although homogenizing and sieving of FC-1 to <75 min obviously obscured variations in stratigraphic distribution, it is unlikely to have affected the size-frequency distribution of Ir host phases.) We previously identified micrometer-scale Ir host phases by electron microscopy in melt-rock samples from two widely separated drill holes at the Chicxulub Basin, including a replicate split of Y6-NI9-R. One is an aggregate of subhedral Ir metal grains enclosed in silicate, in which no other Pt group elements (PGE) were detected. A second particle with twice the mass as the first, concentrated predominantly in a single grain, is associated with minor concentrations of Os, Ru, and Pt, and with adhering particles of corundum and perovskite. A third Ir-rich particle, with a greater apparent Os concentration, was identified before being lost as a result of charging under the electron beam. In addition to demonstrating the preservation of projectile components within the Chicxulub Crater, analogous phase associations in Ca- and Al-rich inclusions (CAI) from C2 and C3 chondrites suggest to us that these melt-rock Ir host phases are relics from a carbonaceous chondrite K/T boundary impactor Although the obviously low Ru/Ir ratios of the Chicxulub Ir host phases are qualitatively consistent with suggested PGE fractionation with distance during condensation in an ejecta cloud, it seems difficult to explain the accumulation of the about 3 x 10(exp 11) Ir atoms required to form a about 10(exp -10) g nugget of pure Ir metal within a jet of vaporized projectile expanding at 1-4 km/s, or to effectively exclude or remove commonly alloyed PGE and siderophile elements by fractionation processes resulting from condensation, oxidation, sulfidization, exsolution, or autometamorphism during cooling of the melt. We do not dismiss the importance of these processes entirely; on the contrary, other geochemical and mineralogical aspects of the melt rocks require them, and condensation from the expanding ejecta cloud appears to best explain the primary Ir host-phase distribution in the fish clay, as well as the high Ir concentrations associated with spinel-bearing spheroids at the K/T boundary in the Pacific Ocean . If the "relict" hypothesis is correct, micronuggets of other PGEs and alloys, not detected by our INAA screening, should also occur in the melt rocks. Possibly, the discrete host phases with lesser Ir masses are such alloys with subordinate Ir, rather than simply smaller, predominantly Ir-bearing particles. A CAI source for the relics would be consistent with either a comet or an asteroid K/T impact at Chicxulub. (Additional information contained in the original.)

Schuraytz, B. C.↗
DOE OSTI2020

Unprecedented non-hysteretic superelasticity of [001]-oriented NiCoFeGa single crystals

Superelasticity associated with the martensitic transformation has found a broad range of engineering applications. However, the intrinsic hysteresis and temperature sensitivity of the first-order phase transformation significantly hinder the usage of smart metallic components in many critical areas. In this paper, we report a large superelasticity up to 15.2% strain in [001]-oriented NiCoFeGa single crystals, exhibiting non-hysteretic mechanical responses, a small temperature dependence and high-energy-storage capability and cyclic stability over a wide temperature and composition range. In situ synchrotron X-ray diffraction measurements show that the superelasticity is correlated with a stress-induced continuous variation of lattice parameter accompanied by structural fluctuation. Neutron diffraction and electron microscopy observations reveal an unprecedented microstructure consisting of atomic-level entanglement of ordered and disordered crystal structures, which can be manipulated to tune the superelasticity. The discovery of the large elasticity related to the entangled structure paves the way for exploiting elastic strain engineering and development of related functional materials. NiCoFeGa single crystals exhibit large non-hysteretic superelasticity over broad temperature and composition ranges. It is attributed to the continuous phase transition with applied stress, which is related to the fluctuation of entangled ordered and disordered crystal structures.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Lattice-Matched Semiconductor Layers on Single Crystalline Sapphire Substrate

SiGe is an important semiconductor alloy for high-speed field effect transistors (FETs), high-temperature thermoelectric devices, photovoltaic solar cells, and photon detectors. The growth of SiGe layer is difficult because SiGe alloys have different lattice constants from those of the common Si wafers, which leads to a high density of defects, including dislocations, micro-twins, cracks, and delaminations. This innovation utilizes newly developed rhombohedral epitaxy of cubic semiconductors on trigonal substrates in order to solve the lattice mismatch problem of SiGe by using trigonal single crystals like sapphire (Al2O3) as substrate to give a unique growth-orientation to the SiGe layer, which is automatically controlled at the interface upon sapphire (0001). This technology is different from previous silicon on insulator (SOI) or SGOI (SiGe on insulator) technologies that use amorphous SiO2 as the growth plane. A cubic semiconductor crystal is a special case of a rhombohedron with the inter-planar angle, alpha = 90 deg. With a mathematical transformation, all rhombohedrons can be described by trigonal crystal lattice structures. Therefore, all cubic lattice constants and crystal planes (hkl) s can be transformed into those of trigonal crystal parameters. These unique alignments enable a new opportunity of perfect lattice matching conditions, which can eliminate misfit dislocations. Previously, these atomic alignments were thought to be impossible or very difficult. With the invention of a new x-ray diffraction measurement method here, growth of cubic semiconductors on trigonal crystals became possible. This epitaxy and lattice-matching condition can be applied not only to SiGe (111)/sapphire (0001) substrate relations, but also to other crystal structures and other materials, including similar crystal structures which have pointgroup rotational symmetries by 120 because the cubic (111) direction has 120 rotational symmetry. The use of slightly miscut (less than plus or minus 10 deg.) sapphire (0001) substrate can be used to improve epitaxial relationships better by providing attractive atomic steps in the epitaxial process.

Choi, Sang↗