Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binary 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 73 records · Page 4

The role of Cr concentration and temperature on cavity swelling with co-injected helium in dual-ion irradiated Fe and Fe-Cr alloys

The level of chromium plays an essential role in irradiation tolerance of Fe-Cr ferritic alloys. However, conflicting results have been reported regarding the dependence of cavity swelling under irradiation on Cr level and temperature. Here, we have performed a comprehensive set of simultaneous dual-ion (Ni + He) irradiations to high dose (~30 displacements per atom, dpa) at 400–550 °C on a series of ultra-high purity Fe and Fe-Cr binary alloys (3–14 wt.%Cr). Helium co-implantation rates of 0.1 and 10 appm He/dpa were selected to examine He synergistic effects relevant for fission and fusion reactor conditions, respectively. Cavities were observed in all irradiated samples by transmission electron microscopy. The results show that higher He implantation rate causes a shift in the swelling peak to higher temperatures in both Fe and Fe-Cr alloys. When assuming smaller cavities as biased sinks, the non-monotonic nature of the cavity swelling behavior is related to the ratio of biased to unbiased point defect sink strengths. Cr-enriched precipitates were observed in Fe-14Cr irradiated at 400 °C by atom probe tomography. Our analysis suggests the formation of Cr-enriched precipitates could suppress cavity swelling for Fe-Cr alloys with Cr content above 10 wt%.

36 MATERIALS SCIENCE↗

A laser ultrasonics-based approach for rapid screening of high entropy alloys

The primary objective of this seed project is to develop a laser ultrasonics-based characterization methodology for rapid metallic materials design and discovery via in situ determination of phases, microstructures and elastic properties with respect to temperature. Ultrasonic waves are strongly affected by material microstructure, and therefore, serve as a facile means to probe elastic properties, phase content and their size distributions and volume fractions. In this study, a laser ultrasonic technique will be used to systematically study the evolution of properties in a set of interrelated simple binary alloys and high entropy alloys (HEA). Phase transformations and microstructural changes inferred from the ultrasonic signals will be correlated with electron microscopy data and predictions using calculations of phase diagrams (CALPHAD). The non-contact and non-destructive ultrasonic testing approach developed in this study could overcome several limitations associated with current material characterization methods for materials discovery. It is expected that the products of the proposed work will be utilized to evaluate novel graded composition HEAs currently being developed at the Idaho National Laboratory (INL) using advanced manufacturing methods based on direct energy deposition and spark plasma sintering processes, and contribute to accelerate the discovery of HEAs.

36 MATERIALS SCIENCE↗

Deformation mechanisms in single crystal Ni-based concentrated solid solution alloys by nanoindentation

Nanoindentation is a critical technique to probe mechanical properties at the micrometer and sub-micrometer scales, accompanied by challenges from indentation size effect, pile-up/sink-in effect, and strain rate sensitivity. In this work, different nanoindentation techniques have been employed to explore Ni-based concentrated solid solution alloys (CSAs) with the addition of 3d transition metal elements including Co, Cr, Mn, and Fe, including unique single-crystal Ni, NiCo, NiFe, Ni 80 Cr 20 , and NiCoFeCr samples with (100) surfaces. A procedure of nanoindentation tests and data analysis/correction have been developed, and a data set of hardness, elastic modulus, strain rate sensitivity, and activation volume for Ni-based CSAs are provided, including the less explored binary alloys such as Ni 80 Cr 20 and Ni 80 Mn 20 . The results show that the type of alloying elements is more critical than the number of elements in strengthening: Co does not provide strengthening in NiCo, while Cr, Mn, and Fe are effective strengthening elements. Cr is the most effective among all the 3d transition metal elements. Furthermore, atomic-level lattice distortion is responsible for the strengthening and the role of stacking fault energy is insignificant in Ni-based CSAs at room temperature. In summary, nanoindentation shows increasing promise as a reliable and fast tool to provide comprehensive mechanical information for new alloy design and development.

36 MATERIALS SCIENCE↗

VZr_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 vanadium-zirconium (V-Zr). 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 V and Zr. 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.

36 MATERIALS SCIENCE↗

TaTi_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-titanium (Ta-Ti). 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 elements Hf and Nb. 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 Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

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↗

Diffusion bonding of tungsten-vanadium-zirconium using vacuum hot pressing for the development of a low decay heat cladding solution for tungsten spallation targets

Tantalum has been used as cladding material for water-cooled solid tungsten targets at many leading spallation neutron production facilities thanks to its high neutron yield, manageable radiation damage behavior, and excellent corrosion/erosion resistance in radiation environments. However, from a safety hazard perspective, thermal neutron capture of tantalum in spallation environments causes a high specific decay heat in the target volume, which often becomes a limiting factor in increasing the beam power on the target. In this paper, we studied vacuum hot pressing (VHP) parameters to diffusion bond zirconium to tungsten to explore the feasibility of using zirconium alloys as an alternative cladding material to tantalum. Zirconium alloys have long been used as cladding material for early generation solid spallation targets, and nuclear fuel rods. In spallation environments zirconium has significantly lower decay heat with shorter decay time compared to tantalum. The hot isostatic pressing (HIP) of zirconium and tungsten is known to produce limited bonding quality due to the formation of the brittle ZrW2 intermetallic layer. To overcome this problem, placing a vanadium interlayer between tungsten and zirconium has been proposed by exploring parameter space in binary alloy phase diagrams. Under the VHP conditions, 860 ◦C at 70 MPa for 4 hours, Zr-V and V-W showed good diffusion bonding, which demonstrates the feasibility of a single step HIP process to make the zirconium alloy clad tungsten spallation volumes.

Lee, Yong Joong↗

HfNb_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 hafnium-niobium (Hf-Nb). 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 elements Hf and Nb. 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,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗

A stress-sensitive precipitate nucleation model beyond classical nucleation theory

The dynamic evolution of precipitates and second phases dictates the strength and stability of most engineering alloys. By design, or as a consequence of thermo-mechanical aging, engineering metals and alloys often form precipitates of second phases when subjecting to diverse thermal and mechanical loads. Precipitation is governed by several factors, including the alloy’s composition, processing/operating temperature, and stresses — either as a result of external loads or from residual stresses. However, state-of-the-art models for precipitate nucleation (i.e., classical nucleation theory) typically lacks consistent method to capture the effects of externally applied and/or internal stresses on nucleation; thereby severely limiting the applicability of these models to complex materials systems and to representative loading scenarios. Here, in this work, we extend upon classical nucleation theory to account for the effect of stresses on precipitation kinetics and thermodynamics. This is achieved via the use of an Eshelbian micromechanics framework keeping track of (i) the stress build up resulting from second phase formation as a function of mechanical load and, (ii) the effects of dislocations on precipitate formation. This new model is applied to σ precipitate in Fe–Cr binary alloys and M 23 C 6 precipitate in 316H stainless steel (SS). Simulations demonstrate the important role of both the remotely applied loads and dislocation pile ups on precipitate nucleation.

36 MATERIALS SCIENCE↗

Structural heterogeneity in non-crystalline Te x Se1−x thin films

Rapid crystallization behavior of amorphous TexSe1−x thin films limits the use of these alloys as coatings and in optoelectronic devices. Understanding the short- and medium-range ordering of the amorphous structure and the fundamental physics governing the crystallization of the films is crucial. Although the lack of long range crystalline order restricts the characterization of the amorphous films, electron microscopy offers a way to extract information about the nanoscale ordering. In this paper, the local ordering of amorphous TexSe1−x thin films with x=0.22, 0.61, 0.70, 0.90, and 1 grown by thermal evaporation is investigated using radial distribution function (RDF) and fluctuation electron microscopy (FEM) analysis. RDF results show that the nearest-neighbor distances of selenium (Se) and tellurium (Te) in their crystalline structure are preserved, and their bond lengths increase with the addition of Te. Density functional theory (DFT) calculations predict structures with interatomic distances similar to those measured experimentally. Additionally, fluctuations in atomic coordination are analyzed. Medium range order (MRO) analysis obtained from FEM and DFT calculations suggests that there are at least two populations within the chain network structure, which are close to the Se–Se and Te–Te intrachain distances. For the binary alloy with x > 0.61, TexSe1−x, Te–Te like populations increase and Te fragments might form, suggesting that the glass forming ability decreases rapidly.

Sari, Bengisu (ORCID:0000000244217098)↗

First Principles Evaluation of Phase Stability in the In-Sn Binary System

Abstract The In-Sn binary alloy system exhibits several unusual features that challenge crystallographic and thermodynamic expectations. We combine first principles total energy calculation with simple thermodynamic modeling to address two key points. First, we evaluate energies along the Bain path to interpret the discontinuous transition between the phases α-In (Pearson type tI2) and β-In 3 Sn (also Pearson type tI2) that are identical in symmetry. Second, we demonstrate that the solid solution phases β-In 3 Sn and γ-InSn 4 (Pearson type hP1) exist at high temperatures only, and they exhibit eutectoid decompositions at low temperatures.

Chemistry↗

Mapping structural heterogeneity at the nanoscale with scanning nano-structure electron microscopy (SNEM)

Here, in this work, we explore the use of scanning electron diffraction (also known as 4D-STEM) coupled with electron atomic pair distribution function analysis (ePDF) to understand the local order (structure and chemistry) as a function of position in a complex multicomponent system, a hot rolled, Ni-encapsulated, Zr 65 Cu 17.5 Ni 10 Al 7.5 bulk metallic glass (BMG), with a spatial resolution of 3 nm. We show that it is possible to gain insight into the chemistry and chemical clustering/ordering tendency in different regions of the sample, including in the vicinity of nano-scale crystallites that are identified from virtual dark field images and in heavily deformed regions at the edge of the BMG. In addition to simpler analysis, unsupervised machine learning was used to extract partial PDFs from the material, modeled as a quasi-binary alloy, and map them in space. These maps allowed key insights not only into the local average composition, as validated by EELS, but also a unique insight into chemical short-range ordering tendencies in different regions of the sample during formation. The experiments are straightforward and rapid and, unlike spectroscopic measurements, don’t require energy filters on the instrument. We spatially map different quantities of interest (QoI’s), defined as scalars that can be computed directly from positions and widths of ePDF peaks or parameters refined from fits to the patterns. We developed a flexible and rapid data reduction and analysis software framework that allows experimenters to rapidly explore images of the sample on the basis of different QoI’s. The power and flexibility of this approach are explored and described in detail. Because of the fact that we are getting spatially resolved images of the nanoscale structure obtained from ePDFs we call this approach scanning nano-structure electron microscopy (SNEM), and we believe that it will be powerful and useful extension of current 4D-STEM methods.

36 MATERIALS SCIENCE↗

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗

Reconstructing parent microstructures in martensitic and pearlitic Ti-Cu

Evaluating the microstructural evolution of parent phases has long been a challenge in metals with a partial or complete solid state transformation. Most parent microstructure evolution has to be inferred from the product microstructure, with some information permanently lost. Parent microstructure reconstructions relying on orientation relationships and EBSD data are a useful tool to overcome this challenge, though much background knowledge is required to do so. Here, this work introduces a reconstruction process for a eutectoid composition Ti-Cu binary alloy to evaluate the high temperature β-Ti phase from martensitic and pearlitic product microstructures. Martensitic microstructures were accurately reconstructed, enabling elucidation of β-Ti grain size and texture, and the employed script included for scientific reference. Pearlitic microstructures did not reconstruct effectively with the current process. However, new insight into the interfacial crystallographic orientations of pearlitic Ti-Cu microstructures was gleaned via this analysis, and future investigations of interest are discussed accordingly. Ultimately, these findings demonstrate tools widely implemented on conventional titanium and ferrous alloys can also be applied to understand parent microstructures in Ti-Cu and other alloys of emerging interest. This is especially of use as new manufacturing processes for such materials require understanding how parent phase microstructures respond to new, novel material processing such as AM.

36 MATERIALS SCIENCE↗

Inferring topological transitions in pattern-forming processes with self-supervised learning

Abstract The identification of transitions in pattern-forming processes are critical to understand and fabricate microstructurally precise materials in many application domains. While supervised methods can be useful to identify transition regimes, they need labels, which require prior knowledge of order parameters or relevant microstructures describing these transitions. Instead, we develop a self-supervised, neural-network-based approach that does not require predefined labels about microstructure classes to predict process parameters from observed microstructures. We show that assessing the difficulty of solving this inverse problem can be used to uncover microstructural transitions. We demonstrate our approach by automatically discovering microstructural transitions in two distinct pattern-forming processes: the spinodal decomposition of a two-phase mixture and the formation of binary-alloy microstructures during physical vapor deposition of thin films. This approach opens a path forward for discovering unseen or hard-to-discern transitions and ultimately controlling complex pattern-forming processes.

Abram, Marcin↗

Computed entropy spectra for grain boundary segregation in polycrystals

Abstract Grain boundary solute segregation influences most bulk material properties, and understanding solute thermodynamics at grain boundaries is critical for engineering them. However, the vast grain boundary space in polycrystals is challenging to evaluate due to its size, especially for the intrinsically hard-to-compute segregation excess entropy. Here data science methods are used to generate a database of site-wise grain boundary segregation entropy spectra for 155 dilute binary alloys within the harmonic approximation. The spectral framework allows scale bridging between the calculated atomistic site-wise energy-entropy spectra and macroscopic segregation entropy estimates. The results affirm that macroscopic averaging is not sufficient: a spectral treatment of grain boundary segregation is needed to accurately model bulk temperature dependence of grain boundary solute segregation. The calculated spectral entropy database and thermodynamic framework can be applied for both understanding segregation experiments and alloy design exercises, paving the way to a finite-temperature grain boundary genome.

Chemistry↗

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗

Stability and Mobility of Disconnections in Solute Atmospheres: Insights from Interfacial Defect Diagrams

This Letter explores the stability of disconnections (step-dislocation defects) at grain boundaries in binary alloys. We introduce interfacial defect diagrams, derived from atomistic simulations and segregation theory, to predict the stability of disconnections in the temperature-solute concentration phase space and relate it to governing segregation mechanisms. These diagrams reveal multiple stability regimes influenced by solute-induced clustering and pinning effects impacting the thermal migration of disconnections and offering insights into their thermodynamics and kinetic properties.

Crystal structure↗