Abstract: Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of machine learning IAPs (ML-IAPs) based on four local environment descriptors—atom-centered symmetry functions (ACSF), smooth overlap of atomic positions (SOAP), the spectral neighbor analysis potential (SNAP) bispectrum components, and moment tensors—using a diverse data set generated using high-throughput density functional theory (DFT) calculations. The data set comprising bcc (Li, Mo) and fcc (Cu, Ni) metals and diamond group IV semiconductors (Si, Ge) is chosen to span a range of crystal structures and bonding. All descriptors studied show excellent performance in predicting energies and forces far surpassing that of classical IAPs, as well as predicting properties such as elastic constants and phonon dispersion curves. We observe a general trade-off between accuracy and the degrees of freedom of each model and, consequently, computational cost. We will discuss these trade-offs in the context of model selection for molecular dynamics and other applications.

Abstract: Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of machine learning IAPs (ML-IAPs) based on four local environment descriptors—atom-centered symmetry functions (ACSF), smooth overlap of atomic positions (SOAP), the spectral neighbor analysis potential (SNAP) bispectrum components, and moment tensors—using a diverse data set generated using high-throughput density functional theory (DFT) calculations. The data set comprising bcc (Li, Mo) and fcc (Cu, Ni) metals and diamond group IV semiconductors (Si, Ge) is chosen to span a range of crystal structures and bonding. All descriptors studied show excellent performance in predicting energies and forces far surpassing that of classical IAPs, as well as predicting properties such as elastic constants and phonon dispersion curves. We observe a general trade-off between accuracy and the degrees of freedom of each model and, consequently, computational cost. We will discuss these trade-offs in the context of model selection for molecular dynamics and other applications.

Abstract: Availability of a reliable interatomic potential is one of the major challenges in utilizing molecular dynamics (MD) for simulations of metals at near the melting temperatures and melting point (MP). Here, we propose a novel approach to address this challenge in the concept of modified-embedded-atom (MEAM) interatomic potential; also, we apply the approach on iron, nickel, copper, and aluminum as case studies. We propose adding experimentally available high temperature elastic constants and MP of the element to the list of typical low temperature properties used for the development of MD interatomic potential parameters. We show that the proposed approach results in a reasonable agreement between the MD calculations of melting properties such as latent heat, expansion in melting, liquid structure factor, and solid-liquid interface stiffness and their experimental/computational counterparts. Then, we present the physical properties of mentioned elements near melting temperatures using the new MEAM parameters. We observe that the behavior of elastic constants, heat capacity and thermal linear expansion coefficient at room temperature compared to MP follows an empirical linear relation (α±β × MP) for transition metals. Furthermore, a linear relation between the tetragonal shear modulus and the enthalpy change from room temperature to MP is observed for face-centered cubic materials.

Abstract: In recent years, efficient interatomic potentials approaching the accuracy of density functional theory (DFT) calculations have been developed using rigorous atomic descriptors satisfying strict invariances, for example, for translation, rotation, permutation of homonuclear atoms, among others. In this paper, we generalize the spectral neighbor analysis potential (SNAP) model to bcc-fcc binary alloy systems. We demonstrate that machine-learned SNAP models can yield significant improvements even over the well-established high-performing embedded atom method (EAM) and modified EAM potentials for fcc Cu and Ni. We also report on the development of a SNAP model for the fcc Ni-bcc Mo binary system by machine learning a carefully constructed large computed data set of elemental and intermetallic compounds. We demonstrate that this binary Ni-Mo SNAP model can achieve excellent agreement with experiments in the prediction of a Ni-Mo phase diagram as well as near-DFT accuracy in the prediction of many key properties, such as elastic constants, formation energies, melting points, etc., across the entire binary composition range. In contrast, the existing Ni-Mo EAM has significant errors in the prediction of the phase diagram and completely fails in binary compounds. This paper provides a systematic model development process for multicomponent alloy systems, including an efficient procedure to optimize the hyperparameters in the model fitting, and paves the way for long-time large-scale simulations of such systems.

Abstract: Primary radiation damage formation in solid materials typically involves collisions between atoms that have up to a few hundred keV of kinetic energy. During these collisions, the distance between two colliding atoms can approach 0.05 nm. At such small atomic separations, force fields fitted to equilibrium properties tend to significantly underestimate the potential energy of the colliding dimer. To enable molecular dynamics simulations of high-energy collisions, it is common practice to use a screened Coulomb force field to describe the interactions and to smoothly join this to the equilibrium force field at a suitable interatomic spacing. However, there is no accepted standard method for choosing the parameters used in the joining process, and our results prove that defect production is sensitive to how the force fields are linked. A new procedure is presented that involves the use of ab initio calculations to determine the magnitude and spatial dependence of the pair interactions at intermediate distances, along with systematic criteria for choosing the joining parameters. Results are presented for the case of nickel, which demonstrate the use and validity of the procedure.

Abstract: The two-phase solid–liquid coexisting structures of Ni, Cu, and Al are studied by molecular dynamics (MD) simulations using the second nearest-neighbor (2NN) modified-embedded atom method (MEAM) potential. For this purpose, the existing 2NN-MEAM parameters for Ni and Cu were modified to make them suitable for the MD simulations of the problems related to the two-phase solid–liquid coexistence of these elements. Using these potentials, we compare calculated low-temperature properties of Ni, Cu, and Al, such as elastic constants, structural energy differences, vacancy formation energy, stacking fault energies, surface energies, specific heat and thermal expansion coefficient with experimental data. The solid–liquid coexistence approach is utilized to accurately calculate the melting points of Ni, Cu, and Al. The MD calculations of the expansion in melting, latent heat and the liquid structure factor are also compared with experimental data. In addition, the solid–liquid interface free energy and surface anisotropy of the elements are determined from the interface fluctuations, and the predictions are compared to the experimental and computational data in the literature.

Abstract: A new interatomic potential for the Ni–Zr system is presented. This potential was developed specifically to match experimental scattering data from Ni, Zr and NiZr2 liquids. Both ab initio and published thermodynamic data were used to optimise the potential to study the liquid and amorphous structure of the NiZr2 alloy. This potential has the C16 phase, being more stable than C11b phase in the NiZr2 alloy, consistent with experiments. The potential leads to the correct glass structure in the molecular dynamics simulation and, therefore, can be used to study the liquid–glass transformation in the NiZr2 alloy.

Abstract: Recent molecular dynamics simulations of the growth of [Ni0.8Fe0.2/Au] multilayers have revealed the formation of misfit-strain-reducing dislocation structures very similar to those observed experimentally. Here we report similar simulations showing the formation of edge dislocations near the interfaces of vapor-deposited (111) [NiFe/CoFe/Cu] multilayers. Unlike misfit dislocations that accommodate lattice mismatch, the dislocation structures observed here increase the mismatch strain energy. Stop-action observations of the dynamically evolving atomic structures indicate that during deposition on the (111) surface of a fcc lattice, adatoms may occupy either fcc sites or hcp sites. This results in the random formation of fcc and hcp domains, with dislocations at the domain boundaries. These dislocations enable atoms to undergo a shift from fcc to hcp sites, or vice versa. These shifts lead to missing atoms, and therefore a later deposited layer can have missing planes compared to a previously deposited layer. This dislocation formation mechanism can create tensile stress in fcc films. The probability that such dislocations are formed was found to quickly diminish under energetic deposition conditions.

Abstract: Modified embedded atom method (MEAM) potentials for fcc elements Cu, Ag, Au, Ni, Pd, Pt, Al, and Pb have been newly developed using the original first nearest-neighbor MEAM and the recently developed second nearest-neighbor MEAM formalisms. It was found that the original MEAM potentials for fcc elements show some critical shortcomings such as structural instability and incorrect surface reconstructions on (100), (110), and/or (111) surfaces. The newly developed MEAM potentials solve most of the problems and describe the bulk properties (elastic constants, structural energy differences), point defect properties (vacancy and interstitial formation energy and formation volume, activation energy of vacancy diffusion), planar defect properties (stacking fault energy, surface energy, surface relaxation and reconstruction), and thermal properties (thermal expansion coefficients, specific heat, melting point, heat of melting) of the fcc elements considered, in good agreement with relevant experimental information. It has been shown that in the MEAM the degree of many-body screening (C_min) is an important material property and that structural stability at finite temperatures should be included as a checkpoint during development of semiempirical potentials.

Abstract: We demonstrate an approach to the development of many-body interatomic potentials for monoatomic metals with improved accuracy and reliability. The functional form of the potentials is that of the embedded-atom method, but the interesting features are as follows: (1) The database used for the development of a potential includes both experimental data and a large set of energies of different alternative crystalline structures of the material generated by ab initio calculations. We introduce a rescaling of interatomic distances in an attempt to improve the compatibility between experimental and ab initio data. (2) The optimum parametrization of the potential for the given database is obtained by alternating the fitting and testing steps. The testing step includes a comparison between the ab initio structural energies and those predicted by the potential. This strategy allows us to achieve the best accuracy of fitting within the intrinsic limitations of the potential model. Using this approach we develop reliable interatomic potentials for Al and Ni. The potentials accurately reproduce basic equilibrium properties of these metals, the elastic constants, the phonon-dispersion curves, the vacancy formation and migration energies, the stacking fault energies, and the surface energies. They also predict the right relative stability of different alternative structures with coordination numbers ranging from 12 to 4. The potentials are expected to be easily transferable to different local environments encountered in atomistic simulations of lattice defects.

Abstract: A detailed derivation of the simplest form of the effective medium theory for bonding in metallic systems is presented, and parameters for the fcc metals Ni, Pd, Pt, Cu, Ag and Au are given. The derivation of parameters is discussed in detail to show how new parameterizations can be made. The method and the parameterization is tested for a number of surface and bulk problems. In particular we present calculations of the energetics of metal atoms deposited on metal surfaces. The calculated energies include heats of adsorption, energies of overlayers, both pseudomorphic and relaxed, as well as energies of atoms alloyed into the first surface layer.

Abstract: This paper addresses the energy associated with the trapping of hydrogen to defects in a nickel lattice. Several dislocations and grain boundaries which occur in nickel are studied. The dislocations include an edge, a screw, and a Lomer dislocation in the locked configuration, i.e. a Lomer-Cottrell lock (LCL). For both the edge and screw dislocations, the maximum trap site energy is approximately 0.1 eV occurring in the region where the lattice is in tension approximately 3-4 angstroms from the dislocation core. For the Lomer-Cottrell lock, the maximum binding energy is 0.33 eV and is located at the core of the a/6(110) dislocation. Several low-index coincident site lattice grain boundaries are investigated, specifically the Sigma 3(112), Sigma 9(221) and Sigma 11(113) tilt boundaries. The boundaries all show a maximum binding energy of approximately 0.25 eV at the tilt boundary. Relaxation of the boundary structures produces an asymmetric atomic structure for both the Sigma 3 and Sigma 9 boundaries and a symmetric structure for the Sigma 11 tilt boundary. The results of this study can be compared to recent experimental studies showing that the activation energy for hydrogen-initiated failure is approximately 0.3-0.4 eV in the Fe-based superalloy IN903. From the results of this comparison it can be concluded that the embrittlement process is likely associated with the trapping of hydrogen to grain boundaries and Lomer-Cottrell locks.

Abstract: The activation energies for self-diffusion of transition metals (Au, Ag, Cu, Ni, Pd, Pt) have been calculated with the Embedded Atom Method (EAM); the results agree well with available experimental data for both mono-vacancy and di-vacancy mechanisms. The EAM was also used to calculate activation energies for vacancy migration near dilute impurities. These energies determine the atomic jump frequencies of the classic "five-frequency formula," which yields the diffusion rates of impurities by a mono-vacancy mechanism. These calculations were found to agree fairly well with experiment and with Neumann and Hirschwald's "Tm" model.

Abstract: Using the approach of Finnis and Sinclair, N-body potentials for copper, silver, gold and nickel have been constructed. The total energy is regarded as consisting of a pair-potential part and a many body cohesive part. Both these parts are functions of the atomic separations only and are represented by cubic splines, fitted to various bulk properties. For the noble metals, the pair-potentials were fitted at short range to pressure-volume relationships calculated by Christensen and Heine so that interactions at separations smaller than that of the first-nearest neighbours can be treated in this scheme. Using these potentials, point defects, surfaces (including the surface reconstructions) and grain boundaries have been studied and satisfactory agreement with available experimental data has been found.

Abstract: A consistent set of embedding functions and pair interactions for use with the embedded-atom method [M.S. Daw and M. I. Baskes, Phys. Rev. B 29, 6443 (1984)] have been determined empirically to describe the fcc metals Cu, Ag, Au, Ni, Pd, and Pt as well as alloys containing these metals. The functions are determined empirically by fitting to the sublimation energy, equilibrium lattice constant, elastic constants, and vacancy-formation energies of the pure metals and the heats of solution of the binary alloys. The validity of the functions is tested by computing a wide range of properties: the formation volume and migration energy of vacancies, the formation energy, formation volume, and migration energy of divacancies and self-interstitials, the surface energy and geometries of the low-index surfaces of the pure metals, and the segregation energy of substitutional impurities to (100) surfaces.

Abstract: The Morse parameters were calculated using experimental values for the energy of vaporization, the lattice constant, and the compressibility. The equation of state and the elastic constants which were computed using the Morse parameters, agreed with experiment for both face-centered and body-centered cubic metals. All stability conditions were also satisfied for both the face-centered and the body-centered metals. This shows that the Morse function can be applied validly to problems involving any type of deformation of the cubic metals.

Abstract: Recent molecular dynamics simulations of the growth of [Ni0.8Fe0.2/Au] multilayers have revealed the formation of misfit-strain-reducing dislocation structures very similar to those observed experimentally. Here we report similar simulations showing the formation of edge dislocations near the interfaces of vapor-deposited (111) [NiFe/CoFe/Cu] multilayers. Unlike misfit dislocations that accommodate lattice mismatch, the dislocation structures observed here increase the mismatch strain energy. Stop-action observations of the dynamically evolving atomic structures indicate that during deposition on the (111) surface of a fcc lattice, adatoms may occupy either fcc sites or hcp sites. This results in the random formation of fcc and hcp domains, with dislocations at the domain boundaries. These dislocations enable atoms to undergo a shift from fcc to hcp sites, or vice versa. These shifts lead to missing atoms, and therefore a later deposited layer can have missing planes compared to a previously deposited layer. This dislocation formation mechanism can create tensile stress in fcc films. The probability that such dislocations are formed was found to quickly diminish under energetic deposition conditions.

Abstract: A detailed derivation of the simplest form of the effective medium theory for bonding in metallic systems is presented, and parameters for the fcc metals Ni, Pd, Pt, Cu, Ag and Au are given. The derivation of parameters is discussed in detail to show how new parameterizations can be made. The method and the parameterization is tested for a number of surface and bulk problems. In particular we present calculations of the energetics of metal atoms deposited on metal surfaces. The calculated energies include heats of adsorption, energies of overlayers, both pseudomorphic and relaxed, as well as energies of atoms alloyed into the first surface layer.

Abstract: The activation energies for self-diffusion of transition metals (Au, Ag, Cu, Ni, Pd, Pt) have been calculated with the Embedded Atom Method (EAM); the results agree well with available experimental data for both mono-vacancy and di-vacancy mechanisms. The EAM was also used to calculate activation energies for vacancy migration near dilute impurities. These energies determine the atomic jump frequencies of the classic "five-frequency formula," which yields the diffusion rates of impurities by a mono-vacancy mechanism. These calculations were found to agree fairly well with experiment and with Neumann and Hirschwald's "Tm" model.

Abstract: A consistent set of embedding functions and pair interactions for use with the embedded-atom method [M.S. Daw and M. I. Baskes, Phys. Rev. B 29, 6443 (1984)] have been determined empirically to describe the fcc metals Cu, Ag, Au, Ni, Pd, and Pt as well as alloys containing these metals. The functions are determined empirically by fitting to the sublimation energy, equilibrium lattice constant, elastic constants, and vacancy-formation energies of the pure metals and the heats of solution of the binary alloys. The validity of the functions is tested by computing a wide range of properties: the formation volume and migration energy of vacancies, the formation energy, formation volume, and migration energy of divacancies and self-interstitials, the surface energy and geometries of the low-index surfaces of the pure metals, and the segregation energy of substitutional impurities to (100) surfaces.

Abstract: An Ag–Ni semi-empirical potential was developed to simulate the segregation of Ni solutes at Ag grain boundaries (GBs). The potential combines a new Ag potential fitted to correctly reproduce the stable and unstable stacking fault energies in this metal and the existing Ni potential from Mendelev et al (2012 Phil. Mag. 92 4454–69). The Ag–Ni cross potential functions were fitted to ab initio data on the liquid structure of the Ag₈₀Ni₂₀ alloy to properly incorporate the Ag–Ni interaction at small atomic separations, and to the Ni segregation energies at different sites within a high-energy Σ9 (221) symmetric tilt GB. By deploying this potential with hybrid Monte Carlo/molecular dynamics simulations, it was found that heterogeneous segregation and clustering of Ni atoms at GBs and twin boundary defects occur at low Ni concentrations, 1 and 2 at%. This behavior is profoundly different from the homogeneous interfacial dispersion generally observed for the Cu segregation in Ag. A GB transformation to amorphous intergranular films was found to prevail at higher Ni concentrations (10 at%). The developed potential opens new opportunities for studying the selective segregation behavior of Ni solutes in interface-hardened Ag metals and its effect on plasticity.

Abstract: Complex concentrated alloys (CCAs) are a new generation of metallic alloys composed of three or more principal elements with physical and mechanical properties that can be tuned by adjusting their compositions. The extensive compositional workspace of CCAs makes it impractical to perform a comprehensive search for a specific material property using experimental measurements. The use of computational methods can rapidly narrow down the search span, improving the efficiency of the design process. We carried out a high-throughput parameterization of modified embedded atom method (MEAM) interatomic potentials for combinations of Cu, Ti, Ni, Cr, Co, Al, Fe, and Mn using a genetic algorithm. Unary systems were parameterized based on DFT calculations and experimental results. MEAM potentials for 28 binary and 56 ternary combinations of the elements were parameterized to DFT results that were carried out with semi-automated frameworks. Specific attention was made to reproduce properties that impact compositional segregation, material strength, and mechanics.

Abstract: A set of embedded atom model (EAM) interatomic potentials was developed to represent highly idealized face-centered cubic (FCC) mixtures of Fe–Ni–Cr–Co–Al at near-equiatomic compositions. Potential functions for the transition metals and their crossed interactions are taken from our previous work for Fe–Ni–Cr–Co–Cu [D. Farkas and A. Caro: J. Mater. Res. 33 (19), 3218–3225, 2018], while cross-pair interactions involving Al were developed using a mix of the component pair functions fitted to known intermetallic properties. The resulting heats of mixing of all binary equiatomic random FCC mixtures not containing Al is low, but significant short-range ordering appears in those containing Al, driven by a large atomic size difference. The potentials are utilized to predict the relative stability of FCC quinary mixtures, as well as ordered L1₂ and B2 phases as a function of Al content. These predictions are in qualitative agreement with experiments. This interatomic potential set is developed to resemble but not model precisely the properties of this complex system, aiming at providing a tool to explore the consequences of the addition of a large size-misfit component into a high entropy mixture that develops multiphase microstructures.

Abstract: Interatomic potentials for the Ni-Co binary and Ni-Al-Co ternary systems have been developed on the basis of the second nearest-neighbor modified embedded-atom method (2NN MEAM) formalism. The potentials describe structural, thermodynamic, deformation and defect properties of solid solution phases or compound phases in reasonable agreements with experiments or first-principles calculations. The results demonstrate the transferability of the potentials and their applicability to large-scale atomistic simulations to investigate the effect of an alloying element, cobalt, on various microstructural factors related to mechanical properties of Ni-based superalloys on an atomic scale.

Abstract: Ni–Al–Co is a promising system for ferromagnetic shape memory applications. This paper reports on the development of a ternary embedded-atom potential for this system by fitting to experimental and first-principles data. Reasonably good agreement is achieved for physical properties between values predicted by the potential and values known from experiment and/or first-principles calculations. The potential reproduces basic features of the martensitic phase transformation from the B2-ordered high-temperature phase to a tetragonal CuAu-ordered low-temperature phase. The compositional and temperature ranges of this transformation and the martensite microstructure predicted by the potential compare well with existing experimental data. These results indicate that the proposed potential can be used for simulations of the shape memory effect in the Ni–Al–Co system.

Abstract: Classical effective potentials are indispensable for any large-scale atomistic simulations, and the relevance of simulation results crucially depends on the quality of the potentials used. For complex alloys such as quasicrystals, however, realistic effective potentials are almost non-existent. We report here our efforts to develop effective potentials especially for quasicrystalline alloy systems. We use the so-called force-matching method, in which the potential parameters are adapted so as to reproduce the forces and energies optimally in a set of suitably chosen reference configurations. These reference data are calculated with ab-initio methods. As a first application, embedded-atom method potentials for decagonal Al–Ni–Co, icosahedral Ca–Cd, and both icosahedral and decagonal Mg–Zn quasicrystals have been constructed. The influence of the potential range and degree of specialization on the accuracy and other properties is discussed and compared.

Abstract: Classical effective potentials are indispensable for any large-scale atomistic simulations, and the relevance of simulation results crucially depends on the quality of the potentials used. For complex alloys such as quasicrystals, however, realistic effective potentials are almost non-existent. We report here our efforts to develop effective potentials especially for quasicrystalline alloy systems. We use the so-called force-matching method, in which the potential parameters are adapted so as to reproduce the forces and energies optimally in a set of suitably chosen reference configurations. These reference data are calculated with ab-initio methods. As a first application, embedded-atom method potentials for decagonal Al–Ni–Co, icosahedral Ca–Cd, and both icosahedral and decagonal Mg–Zn quasicrystals have been constructed. The influence of the potential range and degree of specialization on the accuracy and other properties is discussed and compared.

Abstract: Ordered L1₂ γ′ Ni₃Al intermetallics are essential strengthening components to maintain the high strength of Ni-based superalloys and recently developed high entropy alloys at elevated temperatures. Under service conditions, structural disorder is usually encountered in these intermetallics, resulting in significant loss of their strengthening functionality. Thus, retaining the degree of order of the L1₂ intermetallics is vital for their long-term reliability and serviceability. In this study, atomistic simulations and rate equation analysis are employed to highlight a notable enhancement in the reordering ability of L1₂ intermetallics by incorporating multiple principal elements. Specifically, we examine the effects of Co and Ti addition on the irradiation-induced disordering and kinetic reordering process of L1₂ intermetallics. Our results reveal that the incorporation of Ti in the Al sublattice can maintain comparable disordering resistance as Ni₃Al. Better yet, the introduction of Ti or Co fosters vacancy migration, which accelerates the kinetic reordering rate during the defect diffusion stage. A synergistic effect of Ti and Co in promoting kinetic reordering is also observed. Our work thus suggests a promising approach for designing irradiation-resistant multicomponent intermetallics, which can retain a high degree of structural order under irradiation with chemical disorder contributed by desirable compositions.

Abstract: This paper addresses the energy associated with the trapping of hydrogen to defects in a nickel lattice. Several dislocations and grain boundaries which occur in nickel are studied. The dislocations include an edge, a screw, and a Lomer dislocation in the locked configuration, i.e. a Lomer-Cottrell lock (LCL). For both the edge and screw dislocations, the maximum trap site energy is approximately 0.1 eV occurring in the region where the lattice is in tension approximately 3-4 angstroms from the dislocation core. For the Lomer-Cottrell lock, the maximum binding energy is 0.33 eV and is located at the core of the a/6(110) dislocation. Several low-index coincident site lattice grain boundaries are investigated, specifically the Sigma 3(112), Sigma 9(221) and Sigma 11(113) tilt boundaries. The boundaries all show a maximum binding energy of approximately 0.25 eV at the tilt boundary. Relaxation of the boundary structures produces an asymmetric atomic structure for both the Sigma 3 and Sigma 9 boundaries and a symmetric structure for the Sigma 11 tilt boundary. The results of this study can be compared to recent experimental studies showing that the activation energy for hydrogen-initiated failure is approximately 0.3-0.4 eV in the Fe-based superalloy IN903. From the results of this comparison it can be concluded that the embrittlement process is likely associated with the trapping of hydrogen to grain boundaries and Lomer-Cottrell locks.

Abstract: Second nearest neighbor modified embedded-atom method (2NN-MEAM) interatomic potentials are developed for binary aluminum (Al) alloys applicable from room temperature to the melting point. The binary alloys studied in this work are Al-Cu, Al-Fe and Al-Ni. Sensitivity and uncertainty analyses are performed on potential parameters based on the perturbation approach. The outcome of the sensitivity analysis shows that some of the MEAM parameters interdependently influence all MEAM model outputs, allowing for the definition of an ordered calibration procedure to target specific MEAM outputs. Using these 2NN-MEAM interatomic potentials, molecular dynamics (MD) simulations are performed to calculate low and high-temperature properties, such as the formation energies of stable phases and unstable intermetallics, lattice parameters, elastic constants, thermal expansion coefficients, enthalpy of formation of solids, liquid mixing enthalpy, and liquidus temperatures at a wide range of compositions. The computed data are compared with the available first principle calculations and experimental data, showing high accuracy of the 2NN-MEAM interatomic potentials. In addition, the liquidus temperature of the Al binary alloys is compared to the phase diagrams determined by the CALPHAD method.

Abstract: An interatomic potential for the Ni–Al system is presented within the third-generation charge optimized many-body (COMB3) formalism. The potential has been optimized for Ni₃Al, or the γ' phase in Ni-based superalloys. The formation energies predicted for other Ni–Al phases are in reasonable agreement with first-principles results. The potential further predicts good mechanical properties for Ni₃Al, which includes the values of the complex stacking fault (CSF) and the anti-phase boundary (APB) energies for the (1 1 1) and (1 0 0) planes. It is also used to investigate dislocation propagation across the Ni₃Al (1 1 0)–Ni (1 1 0) interface, and the results are consistent with simulation results reported in the literature. The potential is further used in combination with a recent COMB3 potential for Al₂O₃ to investigate the Ni₃Al (1 1 1)–Al₂O₃ (0 0 0 1) interface, which has not been modeled previously at the classical atomistic level due to the lack of a reactive potential to describe both Ni₃Al and Al₂O₃ as well as interactions between them. The calculated work of adhesion for this interface is predicted to be 1.85 J m⁻², which is in agreement with available experimental data. The predicted interlayer distance is further consistent with the available first-principles results for Ni (1 1 1)–Al₂O₃ (0 0 0 1).

Abstract: We construct an interatomic potential for the Ni-Al system within the embedded-atom method formalism. The potential is based on previously developed accurate potentials for pure Ni and Al. The cross-interactions are fitted to experimental cohesive energy, lattice parameter and elastic constants of B2-NiAl, as well as to ab initio formation energies of several real or imaginary intermetallic compounds with different crystal structures and chemical compositions. The potential accurately reproduces a variety of physical properties of the NiAl and Ni3Al phases, and shows reasonable agreement with experimental and ab initio data for phase stability across the Ni-Al phase diagram. Most of the properties reproduced by the new potential were not involved in the fitting process, which demonstrates its excellent transferability. Advantages and certain weaknesses of the new potential in comparison with other existing potentials are discussed in detail. The potential is expected to be especially suitable for simulations of heterophase interfaces and mechanical behavior of Ni-Al alloys.

Abstract: Complex equilibria and phase transformations involving diffusion can now be calculated quickly and efficiently. Detailed examples are given for cases which involve varying degrees of non-equilibrium and therefore time-dependence. Despite very good agreement between such calculations and experimental results, many potential end-users are still not convinced that such techniques could be usefully applied to their own specific problems. Friendly graphic interface versions of calculating software are now generally available, so the authors conclude that the most likely source of the reluctance to use such tools lies in the formulation of relevant questions and the interpretation of the results. Although the potential impact of such tools was foreseen many years ago [M. Hillert, Calculation of phase equilibria, in: Conference on Phase Transformations, 1968], few changes in the relevant teaching curricula have taken into account the availability and power of such techniques.
This paper has therefore been designed not only as a collection of interesting problems, but also highlights the critical steps needed to achieve a solution. Each example includes a presentation of the "real" problem, any simplifications that are needed for its solution, the adopted thermodynamic formulation, and a critical evaluation of the results. The availability of such examples should facilitate changes in subject matter that will both make it easier for the next generation of students to use these tools, and at the same time reduce the time and effort currently needed to solve such problems by less efficient methods.
The first set of detailed examples includes the deoxidation of steel by aluminum; heat balance calculations associated with ladle additions to steel; the determination of conditions that avoid undesirable inclusions; the role of methane in sintering atmospheres; interface control during the physical vapour deposition of cemented carbide; oxidation of γ-TiAl materials; and simulation of the thermolysis of metallorganic precursors for Si-C-N ceramics and interface reaction of yttrium silicates with SiC-coated C/C-SiC composites for heat shield applications.
A second set of examples, more dependent on competitive nucleation and growth, includes segregation and carburization in multicomponent steels and features a series of sophisticated simulatons using DICTRA software.
Interfacial and strain energies become increasingly important in defining phase nucleation and morphology in such problems, but relatively little information is available compared to free energy and diffusion databases. The final section therefore demonstrates how computational thermodynamics, semi-empirical atomistic approaches and first-principles calculations are being used to aid filling this gap in our knowledge.

Abstract: A new embedded-atom potential has been developed for Ni₃Al by fitting to experimental and first-principles data. The potential describes lattice properties of Ni₃Al, point defects, planar faults, as well as the γ and γ′ fields on the Ni–Al phase diagram. The potential is applied to calculate the energies of coherent Ni/Ni₃Al interphase boundaries with three different crystallographic orientations. Depending on the orientation, the interface energy varies between 12 and 46 mJ/m2. Coherent γ/γ′ interfaces existing at high temperatures are shown to be more diffuse and are likely to have a lower energy than Ni/Ni₃Al interfaces.

Abstract: An embedded-atom potential has been constructed for the intermetallic compound B2−NiAl by fitting to both experimental properties and ab initio data. The ab initio data have been generated in the form of energy-volume relations for a number of alternative structures of NiAl and Ni3Al, as well as for Ni and Al. The potential accurately reproduces the basic lattice properties of B2−NiAl, planar faults, and point-defect characteristics. It also reproduces the energetics and stability of all other structures included in the fit. The potential is applied to calculate equilibrium point-defect concentrations in B2−NiAl as functions of temperature and composition near the stoichiometry. In contrast to previous calculations, the defect formation entropies arising from atomic vibrations are included in our calculation within the quasiharmonic approximation. Such entropies tend to increase the concentrations of thermal point defects in B2−NiAl at high temperatures, but the atomic disorder mechanism remains triple-defect type up to the melting point.

Abstract: An interatomic potential for the Ni–Al system is presented within the third-generation charge optimized many-body (COMB3) formalism. The potential has been optimized for Ni₃Al, or the γ' phase in Ni-based superalloys. The formation energies predicted for other Ni–Al phases are in reasonable agreement with first-principles results. The potential further predicts good mechanical properties for Ni₃Al, which includes the values of the complex stacking fault (CSF) and the anti-phase boundary (APB) energies for the (1 1 1) and (1 0 0) planes. It is also used to investigate dislocation propagation across the Ni₃Al (1 1 0)–Ni (1 1 0) interface, and the results are consistent with simulation results reported in the literature. The potential is further used in combination with a recent COMB3 potential for Al₂O₃ to investigate the Ni₃Al (1 1 1)–Al₂O₃ (0 0 0 1) interface, which has not been modeled previously at the classical atomistic level due to the lack of a reactive potential to describe both Ni₃Al and Al₂O₃ as well as interactions between them. The calculated work of adhesion for this interface is predicted to be 1.85 J m⁻², which is in agreement with available experimental data. The predicted interlayer distance is further consistent with the available first-principles results for Ni (1 1 1)–Al₂O₃ (0 0 0 1).

Abstract: Interatomic potentials for the Ni-Ti and Ni-Al-Ti systems have been developed based on the second nearest-neighbor modified embedded-atom method (2NN-MEAM) formalism. The Ni-Ti binary potential reproduces fundamental materials properties (structural, elastic, thermodynamic, and thermal stability) of alloy systems in reasonable agreement with experiments, first-principles calculations and thermodynamic calculations. Atomistic simulations using the Ni-Al-Ti ternary potential validate that the potential can be applied successfully to atomic-scale investigations to clarify the effects of titanium on important materials phenomena (site preference in γ', γ-γ' phase transition, and segregation on grain boundaries) in Ni-Al-Ti ternary superalloys.

Abstract: A computational approach combining dispersion-corrected density functional theory (DFT) and classical molecular dynamics is employed to characterize the geometrical and thermomechanical properties of a recently proposed 2D transition metal dihalide NiCl₂. The characterization is performed using a classical interatomic force field whose parameters are determined and verified through the comparison with the results of DFT calculations. The developed force field is used to study the mechanical response, thermal stability, melting and solidification of a NiCl₂ monolayer on the atomistic level of detail. The 2D NiCl₂ sheet is found to be thermally stable at temperatures below its melting point of ~695 K. At higher temperatures, several subsequent structural transformations of NiCl₂ are observed, namely a transition into a porous 2D sheet and a 1D nanowire. The MD simulations of NiCl₂ cooling show that the molten NiCl₂ system solidifies into an amorphous porous 2D structure at T ~450 K. The resulting structure has lower cohesive energy with respect to the initial 2D sheet. The computational methodology presented through the case study of NiCl₂ can also be utilized to study the properties of other novel 2D materials, including recently synthesized NiO₂, NiS₂, and NiSe₂.

Abstract: High Entropy Alloys (HEA) attract attention as possible radiation resistant materials, a feature observed in some experiments that has been attributed to several unique properties of HEA, in particular to the disorder-induced reduced thermal conductivity and to the peculiar defect properties originating from the chemical complexity. To explore the origin of such behavior we study the early stages (less than 0.1 ns), of radiation damage response of a HEA using molecular dynamics simulations of collision cascades induced by primary knock-on atoms (PKA) with 10, 20 and 40 keV, at room temperature, on an idealized model equiatomic quinary fcc FeNiCrCoCu alloy, the corresponding "Average Atom" (AA) material, and on pure Ni. We include accurate corrections to describe short-range atomic interactions during the cascade. In all cases the average number of defects in the HEA is lower than for pure Ni, which has been previously used to help claiming that HEA is radiation resistant. However, simulated defect evolution during primary damage, including the number of surviving Frenkel Pairs, and the defect cluster size distributions are nearly the same in all cases, within our statistical uncertainty. The number of surviving FP in the alloy is predicted fairly well by analytical models of defect production in pure materials. All of this indicates that the origin of radiation resistance in HEAs as observed in experiments may not be related to a reduction in primary damage due to chemical disorder, but is probably caused by longer-time defect evolution.

Abstract: A set of embedded atom method model interatomic potentials is presented to represent a high-entropy alloy with five components. The set is developed to resemble but not model precisely face-centered cubic (fcc) near-equiatomic mixtures of Fe–Ni–Cr–Co–Cu. The individual components have atomic sizes deviating up to 3%. With the heats of mixing of all binary equiatomic random fcc mixtures being less than 0.7 kJ/mol and the corresponding value for the quinary being −0.0002 kJ/mol, the potentials predict the random equiatomic fcc quinary mixture to be stable with respect to phase separation or ordering and with respect to bcc and hcp random mixtures. The details of lattice distortion, strain, and stress states in this phase are reported. The standard deviation in the individual nearest neighbor bond lengths was found to be in the range of 2%. Most importantly, individual atoms in the alloy were found to be under atomic strains up to 0.5%, corresponding to individual atomic stresses up to several GPa.

Abstract: The single-phase equiatomic CoCrFeMnNi alloy is a random solid solution of five elements on the face-centered cubic lattice, whose pure constituents crystallize in very different structures and exhibit diverse magnetic properties. Due to the randomness of the alloy, 80% of nearest neighbor bonds are between unlike elements and thus the details of bonding in pure structures are less important. The elastic moduli of this alloy give rise to small Cauchy pressure C₁₂ − C₄₄, which suggests that the dominant part of bonding may be described by a simple pair potential. We test this hypothesis by developing a long-range Lennard-Jones potential in which the equilibrium crystal structures of pure constituents are taken as reference. The standard mixing rules for regular solutions are then adopted to obtain parameters for bonds between unlike elements in the quinary system. The transferability of this potential to quaternary CoCrFeNi, ternary CoCrNi, and binary FeNi alloys is investigated and the predictions compared with experiments and density functional theory calculations. By sampling over a large number of random configurations, we investigate the effect of compositional randomness on misfit volumes, energies of point defects and stacking faults, and the dislocation friction stresses experienced by moving edge and screw dislocations.

Abstract: Although high-entropy alloys (HEAs) are attracting interest, the physical metallurgical mechanisms related to their properties have mostly not been clarified, and this limits wider industrial applications, in addition to the high alloy costs. We clarify the physical metallurgical reasons for the materials phenomena (sluggish diffusion and micro-twining at cryogenic temperatures) and investigate the effect of individual elements on solid solution hardening for the equiatomic CoCrFeMnNi HEA based on atomistic simulations (Monte Carlo, molecular dynamics and molecular statics). A significant number of stable vacant lattice sites with high migration energy barriers exists and is thought to cause the sluggish diffusion. We predict that the hexagonal close-packed (hcp) structure is more stable than the face-centered cubic (fcc) structure at 0 K, which we propose as the fundamental reason for the micro-twinning at cryogenic temperatures. The alloying effect on the critical resolved shear stress (CRSS) is well predicted by the atomistic simulation, used for a design of non-equiatomic fcc HEAs with improved strength, and is experimentally verified. This study demonstrates the applicability of the proposed atomistic approach combined with a thermodynamic calculation technique to a computational design of advanced HEAs.

Abstract: High entropy alloys (HEA) composition-structure relationships are crucial for guiding their design and applications. Here, we use a combined experimental and molecular dynamics (MD) approach to investigate phase formation during physical vapor deposition (PVD) of CoFeNiTi_x and CrFeNiTi_x HEA thin films. We vary titanium molar ratio from 0 to 1 to understand the role of a larger element in the alloy mixture. The experiments show that a high titanium content favors amorphous phase formation in the samples produced by magnetron co-sputtering. In contrast, a low titanium content results in the formation of a face-centered cubic (FCC) structure in both HEA families. This effect of titanium content on the stability of the amorphous and FCC phases is reproduced in PVD MD simulations. The threshold titanium molar ratio is identified to be ~0.53 and ~0.16 for the CoFeNiTi_x and CrFeNiTi_x films in the experiments, and ~0.53 and ~0.53 in the MD simulations. In addition, the atomistic modeling allows for energy versus volume calculations with increasing titanium content, which demonstrate the stabilization of the amorphous phase with respect to crystalline structures. To isolate the effect of atomic sizes, additional simulations are performed using an average-atom model, which disregards differences in atomic radii while preserving the average properties of the alloy. In these simulations, the energetic stability of the amorphous phase disappears. The combined experimental and simulation results demonstrate that the formation of the amorphous phase in HEA thin films generated by PVD is directly caused by the atomic size difference.

Abstract: Alloying of Ni with Fe or Co has been shown to reduce primary damage production under ion irradiation. Similar results have been obtained from classical molecular dynamics simulations of 1, 10, 20, and 40 keV collision cascades in Ni, NiFe, and NiCo. In all cases, a mix of imperfect stacking fault tetrahedra, faulted loops with a 1/3⟨111⟩ Burgers vector, and glissile interstitial loops with a 1/2⟨110⟩ Burgers vector were formed, along with small sessile point defect complexes and clusters. Primary damage reduction occurs by three mechanisms. First, Ni-Co, Ni-Fe, Co-Co, and Fe-Fe short-distance repulsive interactions are stiffer than Ni-Ni interactions, which lead to a decrease in damage formation during the transition from the supersonic ballistic regime to the sonic regime. This largely controls final defect production. Second, alloying decreases thermal conductivity, leading to a longer thermal spike lifetime. The associated annealing reduces final damage production. These two mechanisms are especially important at cascades energies less than 40 keV. Third, at the higher energies, the production of large defect clusters by subcascades is inhibited in the alloys. A number of challenges and limitations pertaining to predictive atomistic modeling of alloys under high-energy particle irradiation are discussed.

Abstract: Interatomic potentials for the Ni-Co binary and Ni-Al-Co ternary systems have been developed on the basis of the second nearest-neighbor modified embedded-atom method (2NN MEAM) formalism. The potentials describe structural, thermodynamic, deformation and defect properties of solid solution phases or compound phases in reasonable agreements with experiments or first-principles calculations. The results demonstrate the transferability of the potentials and their applicability to large-scale atomistic simulations to investigate the effect of an alloying element, cobalt, on various microstructural factors related to mechanical properties of Ni-based superalloys on an atomic scale.

Abstract: Ni–Al–Co is a promising system for ferromagnetic shape memory applications. This paper reports on the development of a ternary embedded-atom potential for this system by fitting to experimental and first-principles data. Reasonably good agreement is achieved for physical properties between values predicted by the potential and values known from experiment and/or first-principles calculations. The potential reproduces basic features of the martensitic phase transformation from the B2-ordered high-temperature phase to a tetragonal CuAu-ordered low-temperature phase. The compositional and temperature ranges of this transformation and the martensite microstructure predicted by the potential compare well with existing experimental data. These results indicate that the proposed potential can be used for simulations of the shape memory effect in the Ni–Al–Co system.

Abstract: An interatomic potential (termed EAM-21) has been developed with the embedded atomic method (EAM) for CrFeMnNi quaternary HEAs. This potential is based on a previously developed potential for CrFeNi ternary alloys. The parameters to develop the potential were determined by fitting to experimental values, density functional theory (DFT) and thermodynamic calculations, to reproduce the main crystal characteristics, namely: the stability of the fcc phase, elastic constants, and stacking fault energy. Its applicability for the study of plastic deformation mechanisms was checked by calculations of behaviour of a ½<1 1 0>1 1 1 edge dislocation in equiatomic quaternary CrFeMnNi alloy, as well as its less-complex subsystems (ternaries, binaries, and pure metals). The calculations were performed in the domain of temperatures between 0 and 900 K; smooth and stable glide of an edge dislocation and fcc phase stability in this temperature range was confirmed. This study demonstrates the suitability of the EAM-21 potential for the analysis of plasticity mechanisms and mechanical properties of CrFeMnNi HEAs.

Abstract: Fe‐Ni‐Cr stainless‐steels are important structural materials because of their superior strength and corrosion resistance. Atomistic studies of mechanical properties of stainless‐steels, however, have been limited by the lack of high‐fidelity interatomic potentials. Here using density functional theory as a guide, we have developed a new Fe‐Ni‐Cr embedded atom method potential. We demonstrate that our potential enables stable molecular dynamics simulations of stainless‐steel alloys at high temperatures, accurately reproduces the stacking fault energy-known to strongly influence the mode of plastic deformation (e.g., twinning vs. dislocation glide vs. cross‐slip)-of these alloys over a range of compositions, and gives reasonable elastic constants, energies, and volumes for various compositions. The latter are pertinent for determining short‐range order and solute strengthening effects. Our results suggest that our potential is suitable for studying mechanical properties of austenitic and ferritic stainless‐steels which have vast implementation in the scientific and industrial communities. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.

Abstract: The predictive power of a classical molecular dynamics simulation is largely determined by the physical validity of its underlying empirical potential. In the case of high-energy collision cascades, it was recently shown that correctly modeling interactions at short distances is necessary to accurately predict primary damage production. An ab initio based framework is introduced for modifying an existing embedded-atom method FeNiCr potential to handle these short-range interactions. Density functional theory is used to calculate the energetics of two atoms approaching each other, embedded in the alloy, and to calculate the equation of state of the alloy as it is compressed. The pairwise terms and the embedding terms of the potential are modified in accordance with the ab initio results. Using this reparametrized potential, collision cascades are performed in Ni₅₀Fe₅₀, Ni₈₀Cr₂₀ and Ni₃₃Fe₃₃Cr₃₃. The simulations reveal that alloying Ni and NiCr to Fe reduces primary damage production, in agreement with some previous calculations. Alloying Ni and NiFe to Cr does not reduce primary damage production, in contradiction with previous calculations.

Abstract: Under irradiation NiSiPCr clusters are formed in high-Cr ferritic martensitic steels as well as in FeCr model alloys. In the literature little is known about the origin and contribution to the hardening of these clusters. In this work we performed density functional theory (DFT) calculations to study the stability of small substitutional NiCr-vacancy clusters and interstitial configurations in bcc Fe. Based on DFT data and experimental considerations a ternary potential for the ferritic FeNiCr system was developed. The potential was applied to study the thermodynamic stability of NiCr clusters by means of Metropolis Monte Carlo (MMC) simulations. The results of our simulations show that Cr and Ni precipitate as separate fractions and suggest only a limited synergetic effect between Ni and Cr. Therefore our results suggest that the NiCrSiP clusters observed in experiments must be the result of other mechanisms than the synergy of Cr and Ni at thermal equilibrium.

Abstract: A semi-empirical interatomic potential formalism, the second-nearest-neighbor modified embedded-atom method (2NN MEAM), has been applied to obtaining interatomic potentials for the Fe-Ni, Cr-Ni and Fe-Cr-Ni systems using previously developed MEAM potentials of Fe and Ni and a newly revised potential of Cr. The potential parameters were determined by fitting the experimental data on the enthalpy of formation or mixing, lattice parameter and elastic constant. The present potentials generally reproduced the fundamental physical properties of the Fe-Ni and Cr-Ni alloys. The enthalpy of formation or mixing of the disordered phase at finite temperature and the enthalpy of mixing of the liquid phase are reasonable in agreements with experiment data and CALPHAD calculations. The potentials can be combined with already-developed MEAM potentials to describe Fe-Cr-Ni-based multicomponent alloys. Moreover, the average diffusivities in the unary, some binary and ternary alloys were simulated based on present potential. Good agreement is obtained in comparison with experimental data.

Abstract: The degradation of austenitic stainless steels in a radiation environment is a known problem for the in-core components of nuclear light water reactors. For a better understanding of the prevailing mechanisms responsible for the materials' degradation, large-scale atomistic simulations are desirable. In this framework and as a follow-up on Bonny et al (2011 Modelling Simul. Mater. Sci. Eng. 19 085008), we developed an embedded atom method type interatomic potential for the ternary FeNiCr system to model the production and evolution of radiation defects. Special attention has been drawn to the Fe10Ni20Cr alloy, whose properties were ensured to be close to those of 316L austenitic stainless steels. The potential is extensively benchmarked against density functional theory calculations and the potential developed in our earlier work. As a first validation, the potential is used in AKMC simulations to simulate thermal annealing experiments in order to determine the self-diffusion coefficients of the components in FeNiCr alloys around the Fe10Ni20Cr composition. The results from these simulations are consistent with experiments, i.e., DCr > DNi > DFe.

Abstract: Austenitic stainless steels are commonly used materials for in-core components of nuclear light water reactors. In service, such components are exposed to harsh conditions: intense neutron irradiation, mechanical and thermal stresses, and aggressive corrosion environment which all contribute to the components' degradation. For a better understanding of the prevailing mechanisms responsible for the materials degradation, large-scale atomistic simulations are desirable. In this framework we developed an embedded atom method type interatomic potential for the ternary FeNiCr system to model movement of dislocations and their interaction with radiation defects. Special attention has been drawn to the Fe-10Ni-20Cr alloy, whose properties were ensured to be close to those of 316L austenitic stainless steel. In particular, the stacking fault energy and elastic constants are well reproduced. The fcc phase for the Fe–10Ni-20Cr random alloy was proven to be stable in the temperature range 0–900 K and under shear strain up to 5%. For the same alloy the stable glide of screw dislocations and stability of Frank loops was confirmed.

Abstract: Recent irradiation experiments have shown that smaller vacancy clusters are observed in Ni₈₀Pd₂₀ compared to Ni₈₀Fe₂₀. Using atomistic calculations, we find that the vacancy energetics are significantly different between the two alloys. Ni₈₀Pd₂₀ has lower vacancy migration barriers and lower vacancy-vacancy binding energies than Ni₈₀Fe₂₀. The consequence of these energetic differences is observed in molecular dynamics (MD) simulations, where despite higher vacancy diffusivity that would help in cluster formation, significantly reduced vacancy clusters are observed in Ni₈₀Pd₂₀ than Ni₈₀Fe₂₀. Calculations show that binding energy decreases and formation energy increases with increasing Ni-Ni bond lengths, and larger Ni-Ni bond lengths are observed in Ni₈₀Pd₂₀ than Ni₈₀Fe₂₀. Thus, the reduced vacancy-vacancy binding and higher formation energy due to longer Ni-Ni bonds in Ni₈₀Pd₂₀ are possibly the underlying reasons for smaller vacancy clusters in Ni₈₀-Pd₂₀ than Ni₈₀Fe₂₀. This study illustrates the unique effects of alloying elements on defect energetics and microstructural evolution in random alloys.

Abstract: In this paper, we present a new quaternary interatomic potential for the NiFeCrPd system, which is an extension on the previous NiFeCr potential. Density functional theory is used to calculate the quantities to be fitted, with particular focus on the energetics of point defects with solutes, for the potential to be used towards understanding radiation damage properties. The potential thus will enable the modeling of multi-elemental solid solution alloys consisting of up to four elements. To test the potential, we have performed atomistic kinetic Monte Carlo simulations to investigate the effect of configurational entropy on the self-diffusion coefficients. The self-diffusion coefficients are found to increase with chemical complexity, contrary to the common postulation of sluggish diffusion in high entropy alloys (HEAs). In addition, we have performed molecular dynamics simulations to elucidate the effect of Pd on vacancy diffusion and clustering in pure Ni and binary alloys. In agreement with recent irradiation experiments, our simulations show that while large vacancy clusters, such as stacking fault tetrahedra, are formed in pure Ni, Ni-Fe and Ni-Cr systems, negligible vacancy clustering is observed in Ni-Pd systems, indicating a possible effect of Pd in reducing cluster sizes. We suggest that this potential will be useful for studying the defect evolution in multi-component HEAs.

Abstract: A new interatomic potential has been developed for the Ni–Cr system in the angular-dependent potential (ADP) format by fitting the potential parameters to a set of experimental and first-principles data. The ADP potential reproduces a wide range of properties of both elements as well as binary alloys with reasonable accuracy, including thermal and mechanical properties, defects, melting points of Ni and Cr, and the Ni–Cr phase diagram. The potential can be used for atomistic simulations of solidification, mechanical behavior and microstructure of the Ni-based and Cr-based phases as well as two-phase alloys.

Abstract: A semi-empirical interatomic potential formalism, the second-nearest-neighbor modified embedded-atom method (2NN MEAM), has been applied to obtaining interatomic potentials for the Fe-Ni, Cr-Ni and Fe-Cr-Ni systems using previously developed MEAM potentials of Fe and Ni and a newly revised potential of Cr. The potential parameters were determined by fitting the experimental data on the enthalpy of formation or mixing, lattice parameter and elastic constant. The present potentials generally reproduced the fundamental physical properties of the Fe-Ni and Cr-Ni alloys. The enthalpy of formation or mixing of the disordered phase at finite temperature and the enthalpy of mixing of the liquid phase are reasonable in agreements with experiment data and CALPHAD calculations. The potentials can be combined with already-developed MEAM potentials to describe Fe-Cr-Ni-based multicomponent alloys. Moreover, the average diffusivities in the unary, some binary and ternary alloys were simulated based on present potential. Good agreement is obtained in comparison with experimental data.

Abstract: Radiation-induced embrittlement of bainitic steels is the lifetime limiting factor of reactor pressure vessels in existing nuclear light water reactors. The primary mechanism of embrittlement is the obstruction of dislocation motion produced by nanometric defect structures that develop in the bulk of the material due to irradiation. In view of improving the predictive capability of existing models it is necessary to understand better the mechanisms leading to the formation of these defects, amongst which the so-called "late blooming phases". In this work we study the stability of the latter by means of density functional theory (DFT) calculations and Monte Carlo simulations based on a here developed quaternary FeCuNiMn interatomic potential. The potential is based on extensive DFT and experimental data. The reference DFT data on solute–solute interaction reveal that, while Mn–Ni pairs and triplets are unstable, larger clusters are kept together by attractive binding energy. The NiMnCu synergy is found to increase the temperature range of stability of solute atom precipitates in Fe significantly as compared to binary FeNi and FeMn alloys. This allows for thermodynamically stable phases close to reactor temperature, the range of stability being, however, very sensitive to composition.

Abstract: We present molecular dynamics (MD) simulations of radiation damage in Fe nanoparticles (NP) and bimetallic FeCu core–shell nanoparticles (CSNP). The CSNP includes a perfect body-centered cubic (bcc) Fe core coated with a face-centered cubic (fcc) Cu shell. Irradiation with Fe Primary Knock-on Atoms (PKA) with energies between 1 and 7 keV leads to point defects, without clustering beyond divacancies and very few slightly larger vacancy clusters, and without interstitial clusters, unlike what happens in bulk at the same PKA energies. The Fe-Cu interface and shell can act as a defect sink, absorbing radiation-induced damage and, therefore, the final number of defects in the Fe core is significantly lower than in the Fe NP. In addition, the Cu shell substantially diminishes the number of sputtered Fe atoms, acting as a barrier for recoil ejection. Structurally, the Cu shell responds to the stress generated by the collision cascade by creating and destroying stacking faults across the shell width, which could also accommodate further irradiation defects. We compare our MD results to Monte Carlo Binary Collision Approximation (BCA) simulations using the SRIM code, for the irradiation of an amorphous 3-layer thin film with a thickness equal to the CSNP diameter. BCA does not include defect recombination, so the number of Frenkel pairs is significantly higher than in MD, as expected. Sputtering yield (Y) is underestimated by BCA, which is also expected since the simulation is for a thin film at normal incidence. We also compare MD defect production to bulk predictions of the analytic Athermal Recombination Corrected Displacements Per Atom (arc-dpa) model. The number of vacancies in the Fe core is only slightly lower than arc-dpa predictions, but the number of interstitials is reduced by about one order of magnitude compared to vacancies, at 5 keV. According to the radiation resistance found for FeCu CSNP in our simulations, this class of nanomaterial could be suitable for developing new radiation-resistant coatings, nanostructured components, and shields for use in extreme environments, for instance, in nuclear energy and astrophysical applications.

Abstract: In recent years, the development of atomistic models dealing with microstructure evolution and subsequent mechanical property change in reactor pressure vessel steels has been recognised as an important complement to experiments. In this framework, a literature study has shown the necessity of many-body interatomic potentials for multi-component alloys. In this paper, we develop a ternary many-body Fe–Cu–Ni potential for this purpose. As a first validation, we used it to perform a simulated thermal annealing study of the Fe–Cu and Fe–Cu–Ni alloys. Good qualitative agreement with experiments is found, although fully quantitative comparison proved impossible, due to limitations in the used simulation techniques. These limitations are also briefly discussed.

Abstract: In this atomistic study on the copper–nickel system, a new embedded-atom alloy potential between copper and nickel is fitted to experimental data on the mixing enthalpy, taking available potentials for the pure components from literature. The resulting phase boundaries of the new potential are in very good agreement with a recent CALPHAD prediction. Using this new potential, a high angle symmetrical tilt Σ5 and a coherent Σ3 twin grain boundary (GB) are chosen for a systematic investigation of equilibrium GB segregation in the semi-grandcanonical ensemble at temperatures from 400 K to 800 K. Applying thermodynamically accurate integration techniques, the GB formation energies are calculated exactly and as an absolute value for every temperature and composition, which also enables the evaluation of GB excess entropies. The thorough thermodynamic model of GBs developed by Frolov and Mishin is excellently confirmed by the simulations quantitatively, if the impact of both segregation and GB tension on the change in GB formation energy is accounted for. In the case of the Σ3 coherent GB, it turns out that the change in GB formation energy at low temperatures is for the most part attributed to the GB tension, while segregation only has a small influence. This demonstrated effect of GB tensions should also be taken into account in the interpretation of experiments.

Abstract: We have developed a semi-empirical and many-body type model potential using a modified charge density profile for Cu–Ni alloys based on the embedded-atom method (EAM) formalism with an improved optimization technique. The potential is determined by fitting to experimental and first-principles data for Cu, Ni and Cu–Ni binary compounds, such as lattice constants, cohesive energies, bulk modulus, elastic constants, diatomic bond lengths and bond energies. The generated potentials were tested by computing a variety of properties of pure elements and the alloy of Cu, Ni: the melting points, alloy mixing enthalpy, lattice specific heat, equilibrium lattice structures, vacancy formation and interstitial formation energies, and various diffusion barriers on the (100) and (111) surfaces of Cu and Ni.

Abstract: A semi-empirical interatomic potential, the MEAM, has been applied to obtain an interatomic potential for the Cu–Ni system, based on the previously developed potentials for pure Cu and Ni. The procedure for the determination of potential parameter values is presented. It is shown that the potential describes the basic thermodynamic properties and alloy behaviors of the fcc solid solution (enthalpy of mixing, miscibility gap and lattice parameter) in good agreement with CALPHAD calculation and experimental information. It is also shown how the CALPHAD calculation (enthalpy of mixing) can be used for optimization of the interatomic potential parameters.

Abstract: The surface composition of Ni-Cu alloys has been calculated as a function of atomic layer, crystal face, and bulk composition at a temperature of 800 K. The results show that the composition varies nonmonotonically near the surface with the surface layer strongly enriched in Cu while the near-surface layers are enriched in Ni. The calculations use the embedded-atom method [M. S. Daw and M. I. Baskes, Phys. Rev. B 29, 6443 (1984)] in conjunction with Monte Carlo computer simulations. The embedding functions and pair interactions needed to describe Ni-Cu alloys are developed and applied to the calculation of bulk energies, lattice constants, and short-range order. The heats of segregation are computed for the dilute limit, and the composition profile is obtained for the (100), (110), and (111) surfaces for a variety of bulk compositions. The results are found to be in accord with experimental data.

Abstract: Ni is the second most abundant element in the Earth’s core. Yet, its effects on the inner core’s structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe’s melting point at inner core pressures. The melting temperature of Ni is shown to be 700 to 800 K higher than that of Fe at 323 to 360 GPa. hcp, bcc, and liquid phase relations differ for Fe and Ni. Ni can be a bcc stabilizer for Fe at high temperatures and inner core pressures. A small amount of Ni can accelerate Fe’s crystallization at core pressures. These results suggest that Ni may substantially impact the structure and formation process of the solid inner core.

Abstract: A semi-empirical interatomic potential formalism, the second-nearest-neighbor modified embedded-atom method (2NN MEAM), has been applied to obtaining interatomic potentials for the Fe-Ni, Cr-Ni and Fe-Cr-Ni systems using previously developed MEAM potentials of Fe and Ni and a newly revised potential of Cr. The potential parameters were determined by fitting the experimental data on the enthalpy of formation or mixing, lattice parameter and elastic constant. The present potentials generally reproduced the fundamental physical properties of the Fe-Ni and Cr-Ni alloys. The enthalpy of formation or mixing of the disordered phase at finite temperature and the enthalpy of mixing of the liquid phase are reasonable in agreements with experiment data and CALPHAD calculations. The potentials can be combined with already-developed MEAM potentials to describe Fe-Cr-Ni-based multicomponent alloys. Moreover, the average diffusivities in the unary, some binary and ternary alloys were simulated based on present potential. Good agreement is obtained in comparison with experimental data.

Abstract: A many-body interatomic potential for the Fe–Ni system is fitted, capable of describing both the ferritic and austenitic phase. The Fe–Ni system exhibits two stable ordered intermetallic phases, namely, L10 FeNi and L12 FeNi3, that are key issues to be tackled when creating a Fe–Ni potential consistent with thermodynamics. A procedure, based on a rigid lattice Ising model and the theory of correlation functions space, is developed to address all the intermetallics that are possible ground states of the system. While controlling the ground states of the system, the mixing enthalpy and defect properties were fitted. Both bcc and fcc defect properties are compared with density functional theory calculations and other potentials found in the literature. Finally, the potential is thermodynamically validated by constructing the alloy phase diagram. It is shown that the experimental phase diagram is reproduced reasonably well and that our potential gives a globally improved description of the Fe–Ni system in the whole concentration range with respect to the potentials found in the literature.

Abstract: First-principles calculations of the energy of various crystal structures of Fe, Ni and ordered Fe–Ni compounds with different stoichiometries have been performed by the linearized augmented plane wave (LAPW) method in the generalized gradient approximation. The most stable compounds are L12–Ni3Fe, L10–FeNi, C11f–Ni2Fe and C11f–Fe2Ni. The L12-Ni3Fe compound has the largest negative formation energy, which is consistent with the experimental Fe–Ni phase diagram. The L10–FeNi compound has also been observed experimentally in meteorite samples as a metastable phase. It is suggested here that the C11f compounds could also form in Fe–Ni alloys at low temperatures. A new semi-empirical interatomic potential has been developed for the Fe–Ni system by fitting to experimental data and the results of the LAPW calculations. Recognizing the significance of the covalent component of bonding in this system, the potential is based on the embedded-atom method (EAM) but additionally includes a bond-angle dependence. In comparison with the existing modified EAM method, our potential form is simpler, extends interactions to several (3–5) coordination shells and replaces the screening procedure by a smooth cutoff of the potential functions. The potential reproduces a variety of properties of Fe and Ni with a reasonable accuracy. It also reproduces all stability trends across the Fe–Ni system established by the LAPW calculations. The potential can be useful in atomistic simulations of the phases of the Fe–Ni system.

Abstract: Hydrogen ingress into a metal is a persistent source of embrittlement. Fracture surfaces are often intergranular, suggesting favorable cleave crack growth along grain boundaries (GBs) as one driver for embrittlement. Here, atomistic simulations are used to investigate the effects of segregated hydrogen on the behavior of cracks along various symmetric tilt grain boundaries in fcc Nickel. An atomistic potential for Ni–H is first recalibrated against new quantum level computations of the energy of H in specific sites within the NiΣ5(120)⟨100⟩ GB. The binding energy of H atoms to various atomic sites in the NiΣ3(111) (twin), NiΣ5(120)⟨100⟩, NiΣ99(557)⟨110⟩, and NiΣ9(221)⟨110⟩ GBs, and to various surfaces created by separating these GBs into two possible fracture surfaces, are computed and used to determine equilibrium H concentrations at bulk H concentrations typical of embrittlement in Ni. Mode I fracture behavior is then studied, examining the influence of H in altering the competition between dislocation emission (crack blunting; “ductile” behavior) and cleavage fracture (“brittle” behavior) for intergranular cracks. Simulation results are compared with theoretical predictions (Griffith theory for cleavage; Rice theory for emission) using the computed surface energies. The deformation behavior at the GBs is, however, generally complex and not as simple as cleavage or emission at a sharp crack tip, which is not unexpected due to the complexity of the GB structures. In cases predicted to emit dislocations from the crack tip, the presence of H atoms reduces the critical load for emission of the dislocations and no cleavage is found. In the cases predicted to cleave, the presence of H atoms reduces the cleavage stress intensity and makes cleavage easier, including NiΣ9(221)⟨110⟩ which emits dislocations in the absence of H. Aside from the one unusual NiΣ9(221)⟨110⟩ case, no tendency is found for H to cause a ductile-to-brittle transformation for cracks along GBs in Ni, either according to theory or simulation for initial equilibrium H segregation and with no, or limited, H diffusion near the newly-created fracture surfaces. The NiΣ3(111) twin boundary does not absorb H at all, suggesting that embrittlement is more difficult in materials with higher fraction of such twin boundaries, as found experimentally. Experimental observations of cleavage-like failure are thus presumably caused by mechanisms involving H diffusion or dynamic crack growth.

Abstract: Second nearest-neighbor modified embedded-atom method (MEAM) interatomic potentials for the Al-H and Ni-H binary systems have been developed on the basis of previously developed MEAM potentials of pure Al, Ni, and H. The potentials can describe various fundamental physical properties of the relevant binary alloys (structural, thermodynamic, defect, and dynamic properties of metastable hydrides or hydrogen in face-centered cubic solid solutions) in good agreement with experiments or first-principles calculations. The applicability of the present potentials to atomic level investigations of dynamic behavior of hydrogen atoms in metal membranes is also discussed.

Abstract: A semi-empirical methodology for predicting the permeability of hydrogen in metallic alloys is proposed by combining an atomistic simulation and a thermodynamic calculation. An atomistic simulation based on a modified embedded-atom method interatomic potential and a CALPHAD-type thermodynamic calculation technique was used to predict the diffusivity and solubility of hydrogen, respectively. The approach was applied to the prediction of the hydrogen permeability in V–Al and V–Ni alloys that are promising for non-Pd hydrogen separation membranes. The predicted permeability of hydrogen decreases, as Al or Ni concentration increases in the alloys. The predicted permeability is in quite good agreement with experimental data available in literature, successfully reproducing the overall trend for the effect of alloying elements, which enables an alloy design of metallic hydrogen permeable membranes.

Abstract: Interatomic potentials for the Co-Cr, Co-Fe, Co-Mn, Cr-Mn and Mn-Ni binary systems have been developed in the framework of the second nearest-neighbor modified embedded-atom method (2NN MEAM) formalism. The potentials describe various fundamental alloy behaviors (structural, elastic and thermodynamic behavior of solution and compound phases), mostly in reasonable agreements with experimental data or first-principles calculations. The potentials can be utilized to complete the interatomic potential for the CoCrFeMnNi alloy and to investigate the atomic scale physical metallurgy of high entropy alloys.

Abstract: In recent years, efficient interatomic potentials approaching the accuracy of density functional theory (DFT) calculations have been developed using rigorous atomic descriptors satisfying strict invariances, for example, for translation, rotation, permutation of homonuclear atoms, among others. In this paper, we generalize the spectral neighbor analysis potential (SNAP) model to bcc-fcc binary alloy systems. We demonstrate that machine-learned SNAP models can yield significant improvements even over the well-established high-performing embedded atom method (EAM) and modified EAM potentials for fcc Cu and Ni. We also report on the development of a SNAP model for the fcc Ni-bcc Mo binary system by machine learning a carefully constructed large computed data set of elemental and intermetallic compounds. We demonstrate that this binary Ni-Mo SNAP model can achieve excellent agreement with experiments in the prediction of a Ni-Mo phase diagram as well as near-DFT accuracy in the prediction of many key properties, such as elastic constants, formation energies, melting points, etc., across the entire binary composition range. In contrast, the existing Ni-Mo EAM has significant errors in the prediction of the phase diagram and completely fails in binary compounds. This paper provides a systematic model development process for multicomponent alloy systems, including an efficient procedure to optimize the hyperparameters in the model fitting, and paves the way for long-time large-scale simulations of such systems.

Abstract: The state-of-the-art experimental and atomistic simulation techniques were utilized to study the structure of the liquid and amorphous Ni62Nb38 alloy. First, the ab initio molecular dynamics (AIMD) simulation was performed at rather high temperature where the time limitations of the AIMD do not prevent to reach the equilibrium liquid structure. A semi-empirical potential of the Finnis-Sinclair (FS) type was developed to almost exactly reproduce the AIMD partial pair correlation functions (PPCFs) in a classical molecular dynamics simulation. This simulation also showed that the FS potential well reproduces the bond angle distributions. The FS potential was then employed to elongate the AIMD PPCFs and determine the total structure factor (TSF) which was found to be in excellent agreement with X-ray TSF obtained within the present study demonstrating the reliability of the AIMD for the simulation of the structure of the liquid Ni–Nb alloys as well as the reliability of the developed FS potential. The glass structure obtained with the developed potential was also found to be in excellent agreement with the X-ray data. The analysis of the structure revealed that a network of the icosahedra clusters centered on Ni atoms is forming during cooling the liquid alloy down to Tg and the Nb Z14, Z15, and Z16 clusters are attached to this network. This network is the main feature of the Ni62Nb38 alloy and further investigations of the properties of this alloy should be based on study of the behavior of this network.

Abstract: A new high-precision angular-dependent potential of the Ni-Pd system was obtained by fitting the experimental data and first-principles calculations. Then, the element segregation characteristics and thermal stability of Ni-Pd bimetallic nanoparticles were investigated by Monte Carlo method and molecular dynamics method. The results show that the chemical ordering pattern of Pd_xNi_1-x nanoparticle is the result of the competition of surface energy, strain energy, bond energy and interface energy. When a small amount of Pd atoms are substitutionally doped into the Ni nanoparticle, all the Pd atoms will be segregated on the surface and dispersed. The synergistic effect of Ni atoms and Pd atoms in the surface will improve the catalytic activity and carbon deposition resistance of Pd_xNi_1-x nanoparticle catalyst in methane dry reforming reaction. Increasing the doping amount of Pd atoms will gradually reduce the melting point of Pd_xNi_1-x nanoparticle, thereby reducing its sintering resistance.

Abstract: Palladium (Pd) has attracted attention as one of the major components of noble metal catalysts due to its excellent reactivity to a wide range of catalytic reactions. It is important to predict and control the atomic arrangement in catalysts because their properties are known to be affected by it. Therefore, to enable atomistic simulations for investigating the atomic scale structural evolution, we have developed interatomic potentials for Pd-M (M = Al, Co, Cu, Fe, Mo, Ni, Ti) binary systems based on the second nearest-neighbor modified embedded-atom method formalism. These potentials reproduce various fundamental properties of the alloys (the structural, elastic and thermodynamic properties of compound and solution phases, and order-disorder transition temperature) in reasonable agreements with experimental data, first-principles calculations and CALPHAD assessments. Herein, we propose that these potentials can be applied to the design of robust bimetallic catalysts by predicting the shape and atomic arrangement of Pd bimetallic nanoparticles.

Abstract: Energy transfer between lattice atoms and electrons is an important channel of energy dissipation during displacement cascade evolution in irradiated materials. On the assumption of small atomic displacements, the intensity of this transfer is controlled by the strength of electron–phonon (el–ph) coupling. The el–ph coupling in concentrated Ni-based alloys was calculated using electronic structure results obtained within the coherent potential approximation. It was found that Ni_0.5Fe_0.5, Ni_0.5Co_0.5 and Ni_0.5Pd_0.5 are ordered ferromagnetically, whereas Ni_0.5Cr_0.5 is nonmagnetic. Since the magnetism in these alloys has a Stoner-type origin, the magnetic ordering is accompanied by a decrease of electronic density of states at the Fermi level, which in turn reduces the el–ph coupling. Thus, the el–ph coupling values for all alloys are approximately 50% smaller in the magnetic state than for the same alloy in a nonmagnetic state. As the temperature increases, the calculated coupling initially increases. After passing the Curie temperature, the coupling decreases. The rate of decrease is controlled by the shape of the density of states above the Fermi level. Introducing a two-temperature model based on these parameters in 10 keV molecular dynamics cascade simulation increases defect production by 10–20% in the alloys under consideration.

Abstract: Interatomic potentials for Pt-M (M = Al, Co, Cu, Mo, Ni, Ti, V) binary systems have been developed on the basis of the second nearest-neighbor modified embedded-atom method (2NN MEAM) formalism. The parameters of pure Mo have also been newly developed to solve a problem in the previous 2NN MEAM potential in which the sigma and α-Mn structures become more stable than the bcc structure. The potentials reproduce various materials properties of alloys (structural, thermodynamic and order-disorder transition temperature) in reasonable agreements with relevant experimental data and other calculations. The applicability of the developed potentials to atomistic investigations for the shape and atomic configuration of Pt bimetallic nanoparticles is demonstrated.

Abstract: A new angular-dependent potential (ADP) of Ni-Rh system was obtained by fitting the experimental data and first principle data, and the effectiveness of the potential was tested. Then, the element segregation characteristics and thermal stability of Ni-Rh nanoparticles were studied by Monte Carlo and molecular dynamics. The results show that the chemical ordering pattern of Ni_1-xRh_x nanoparticles is the result of the competition of surface energy, strain energy, interface energy and bond energy. With the increase of x, Rh atoms are preferentially segregated to the surface and dispersed. The concentration of Rh atoms in the surface decreases with the increase of size or temperature. With the increase of x, the melting point of Ni_1-xRh_x nanoparticle first gradually increased, reached the highest near x = 0.1, then gradually decreased, reached the lowest near x = 0.5, and then gradually increased. The above results theoretically explain the reason why doping a small amount of Rh can improve the coking-resistance and sintering-resistance ability of Ni catalyst.

Abstract: Nickel-titanium (NiTi) shape memory alloys (SMA) are widely used, however simulating the martensitic transformation of NiTi from first principles remains challenging. In this work, we developed a neural network interatomic potential (NNIP) for near-equiatomic Ni-Ti system through active-learning based acquisitions of density functional theory (DFT) training data, which achieves state-of-the-art accuracy. Phonon dispersion and potential-of-mean-force calculations of the temperature-dependent free energy have been carried out. This NNIP predicts temperature-induced, stress-induced, and deformation twinning-induced martensitic transformations from atomic simulations, in significant agreement with experiments. The NNIP can directly simulate the superelasticity of NiTi nanowires, providing a tool to guide their design.

Abstract: Nickel-titanium (NiTi) shape memory alloys (SMA) are widely used, however simulating the martensitic transformation of NiTi from first principles remains challenging. In this work, we developed a neural network interatomic potential (NNIP) for near-equiatomic Ni-Ti system through active-learning based acquisitions of density functional theory (DFT) training data, which achieves state-of-the-art accuracy. Phonon dispersion and potential-of-mean-force calculations of the temperature-dependent free energy have been carried out. This NNIP predicts temperature-induced, stress-induced, and deformation twinning-induced martensitic transformations from atomic simulations, in significant agreement with experiments. The NNIP can directly simulate the superelasticity of NiTi nanowires, providing a tool to guide their design.

Abstract: We developed new interatomic potentials, based on the second nearest-neighbor modified embedded-atom method (2NN-MEAM) formalism, for Ti, Ni, and the binary Ti–Ni system. These potentials were fit to melting points, latent heats, the binary phase diagrams for the Ti rich and Ni rich regions, and the liquid phase enthalpy of mixing for binary alloys, therefore they are particularly suited for calculations of crystal-melt (CM) interface thermodynamic and transport properties. The accuracy of the potentials for pure Ti and pure Ni were tested against both 0 K and high temperature properties by comparing various properties obtained from experiments or density functional theory calculations including structural properties, elastic constants, point-defect properties, surface energies, temperatures and enthalpies of phase transformations, and diffusivity and viscosity in the liquid phase. The fitted binary potential for Ti–Ni was also tested against various non-fitted properties at 0 K and high temperatures including lattice parameters, formation energies of different intermetallic compounds, and the temperature dependence of liquid density at various concentrations. The CM interfacial free energies obtained from simulations, based on the newly developed Ti–Ni potential, show that the bcc alloys tend to have smaller anisotropy compared with fcc alloys which is consistent with the finding from the previous studies comparing single component bcc and fcc materials. Moreover, the interfacial free energy and its anisotropy for Ti-2 atom% Ni were also used to parameterize a 2D phase field (PF) model utilized in solidification simulations. The PF simulation predictions of microstructure development during solidification are in good agreement with a geometric model for dendrite primary arm spacing.

Abstract: Interatomic potentials for the Ni-Ti and Ni-Al-Ti systems have been developed based on the second nearest-neighbor modified embedded-atom method (2NN-MEAM) formalism. The Ni-Ti binary potential reproduces fundamental materials properties (structural, elastic, thermodynamic, and thermal stability) of alloy systems in reasonable agreement with experiments, first-principles calculations and thermodynamic calculations. Atomistic simulations using the Ni-Al-Ti ternary potential validate that the potential can be applied successfully to atomic-scale investigations to clarify the effects of titanium on important materials phenomena (site preference in γ', γ-γ' phase transition, and segregation on grain boundaries) in Ni-Al-Ti ternary superalloys.

Abstract: Phase transitions in nickel-titanium shape-memory alloys are investigated by means of atomistic simulations. A second nearest-neighbor modified embedded-atom method interatomic potential for the binary nickel-titanium system is determined by improving the unary descriptions of pure nickel and pure titanium, especially regarding the physical properties at finite temperatures. The resulting potential reproduces accurately the hexagonal-close-packed to body-centered-cubic phase transition in Ti and the martensitic B2−B19′ transformation in equiatomic NiTi. Subsequent large-scale molecular-dynamics simulations validate that the developed potential can be successfully applied for studies on temperature- and stress-induced martensitic phase transitions related to core applications of shape-memory alloys. A simulation of the temperature-induced phase transition provides insights into the effect of sizes and constraints on the formation of nanotwinned martensite structures with multiple domains. A simulation of the stress-induced phase transition of a nanosized pillar indicates a full recovery of the initial structure after the loading and unloading processes, illustrating a superelastic behavior of the target system.

Abstract: We study the properties of NiTi shape-memory nanoparticles coherently embedded in TiV matrices using three-dimensional atomistic simulations based on the modified embedded-atom method. To this end, we develop and present a suitable NiTiV potential for our simulations. Employing this potential, we identify the conditions under which the martensitic phase transformation of such a nanoparticle is triggered—specifically, how these conditions can be tuned by modifying the size of the particle, the composition of the surrounding matrix, or the temperature and strain state of the system. Using these insights, we establish how the transformation temperature of such particles can be influenced and discuss the practical implications in the context of shape-memory strengthened alloys.

Abstract: A semi-empirical methodology for predicting the permeability of hydrogen in metallic alloys is proposed by combining an atomistic simulation and a thermodynamic calculation. An atomistic simulation based on a modified embedded-atom method interatomic potential and a CALPHAD-type thermodynamic calculation technique was used to predict the diffusivity and solubility of hydrogen, respectively. The approach was applied to the prediction of the hydrogen permeability in V–Al and V–Ni alloys that are promising for non-Pd hydrogen separation membranes. The predicted permeability of hydrogen decreases, as Al or Ni concentration increases in the alloys. The predicted permeability is in quite good agreement with experimental data available in literature, successfully reproducing the overall trend for the effect of alloying elements, which enables an alloy design of metallic hydrogen permeable membranes.