• Citation: R.S. Elliott, and A. Akerson (2015), "Efficient "universal" shifted Lennard-Jones model for all KIM API supported species".

    Notes: This is the Co interaction from the "Universal" parameterization for the openKIM LennardJones612 model driver.The parameterization uses a shifted cutoff so that all interactions have a continuous energy function at the cutoff radius. This model was automatically fit using Lorentz-Berthelotmixing rules. It reproduces the dimer equilibrium separation (covalent radii) and the bond dissociation energies. It has not been fitted to other physical properties and its ability to model structures other than dimers is unknown. See the README and params files on the KIM model page for more details.

  • See Computed Properties
    Notes: Listing found at https://openkim.org.
    Link(s):
  • Citation: X.W. Zhou, R.A. Johnson, and H.N.G. Wadley (2004), "Misfit-energy-increasing dislocations in vapor-deposited CoFe/NiFe multilayers", Physical Review B, 69(14), 144113. DOI: 10.1103/physrevb.69.144113.
    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.

    Related Models:
  • FORTRAN (2004--Zhou-X-W--Co--FORTRAN--ipr1)
    Notes: These are the original files sent by X.W. Zhou (Sandia National Laboratory) and posted with his permission. C.A. Becker (NIST) modified create.f to include the reference in the generated potential files and the EAM.input file for this composition. These files can be used to generate alloy potentials for Cu, Ag, Au, Ni, Pd, Pt, Al, Pb, Fe, Mo, Ta, W, Mg, Co, Ti, and Zr by editing EAM.input. However, as addressed in the reference, these potentials were not designed for use with metal compounds.
    File(s): superseded


  • LAMMPS pair_style eam/alloy (2004--Zhou-X-W--Co--LAMMPS--ipr1)
    See Computed Properties
    Notes: This file was generated by C.A. Becker (NIST) from create.f and posted with X.W. Zhou's (Sandia National Laboratory) permission.
    File(s): superseded


  • FORTRAN (2004--Zhou-X-W--Co--FORTRAN--ipr2)
    Notes: The file Zhou04_create_v2.f is an updated version of create.f modified by L.M. Hale (NIST) following advice from X.W. Zhou (Sandia National Laboratory). This version removes spurious fluctuations in the tabulated functions of the original potential files caused by single/double precision floating point number conflicts.
    File(s):
  • LAMMPS pair_style eam/alloy (2004--Zhou-X-W--Co--LAMMPS--ipr2)
    See Computed Properties
    Notes: This file was generated by L.M. Hale from Zhou04_create_v2.f on 13 April 2018 and posted with X.W. Zhou's (Sandia National Laboratory) permission. This version corrects an issue with spurious fluctuations in the tabulated functions.
    File(s):
  • See Computed Properties
    Notes: Listing found at https://openkim.org. This KIM potential is based on the files from 2004--Zhou-X-W--Co--LAMMPS--ipr1.
    Link(s):
  • See Computed Properties
    Notes: Listing found at https://openkim.org. This KIM potential is based on the files from 2004--Zhou-X-W--Co--LAMMPS--ipr2.
    Link(s):
 
  • Citation: X.W. Zhou, R.A. Johnson, and H.N.G. Wadley (2004), "Misfit-energy-increasing dislocations in vapor-deposited CoFe/NiFe multilayers", Physical Review B, 69(14), 144113. DOI: 10.1103/physrevb.69.144113.
    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.

    Notes: This is a combined potential that contains all 16 elements from the source reference. It is provided here due to various requests for more elemental combinations often for high entropy simulations. As a caution, note that all of the cross interactions are determined through a universal mixing function and that most elemental systems were not thoroughly explored and tested by the original authors meaning that most binary and higher-order systems may not be well optimized.

  • See Computed Properties
    Notes: This file was generated by Ilia Nikiforov using the Zhou04_create_v2.f FORTRAN code which can be found on the associated elemental listings. The code was slightly modified to increase the tabulation points to 3000 to ensure good interpolations of the embedding energy function for all elements as W has a noticeably larger delta rho than the other elements. Also, the header was fixed to include all 16 element symbol tags.
    File(s):
 
  • Citation: C. Vailhé, and D. Farkas (1997), "Shear faults and dislocation core structures in B2 CoAl", Journal of Materials Research, 12(10), 2559-2570. DOI: 10.1557/jmr.1997.0340.
    Abstract: Interatomic potentials of the embedded atom and embedded defect type were derived for the Co–Al system by empirical fitting to the properties of the B2 CoAl phase. The embedded atom potentials reproduced most of the properties needed, except that, in using this method, the elastic constants cannot be fitted exactly because CoAl has a negative Cauchy pressure. In order to overcome this limitation and fit the elastic constants correctly, angular forces were added using the embedded defect technique. The effects of angular forces to the embedded atom potentials were seen in the elastic constants, particularly C44. Planar fault energies changed up to 30% in the 110 and 112 γ surfaces and the vacancy formation energies were also very sensitive to the non-central forces. Dislocation core structures and Peierls stress values were computed for the 〈100〉 and 〈111〉 dislocations without angular forces. As a general result, the dislocations with a planar core moved for critical stress values below 250 MPa in contrast with the nonplanar cores for which the critical stress values were above 1500 MPa. The easiest dislocations to move were the 1/2〈111〉 edge superpartials, and the overall preferred slip plane was 110. These results were compared with experimental observations in CoAl and previously simulated dislocations in NiAl.

    Related Models:
  • See Computed Properties
    Notes: Listing found at https://openkim.org.
    Link(s):
 
 
 
 
 
 
 
  • Citation: P. Brommer, and F. Gähler (2006), "Effective potentials for quasicrystals fromab-initiodata", Philosophical Magazine, 86(6-8), 753-758. DOI: 10.1080/14786430500333349.
    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.

    Notes: This is for the Potential A model described in the reference

    Related Models:
  • See Computed Properties
    Notes: Listing found at https://openkim.org.
    Link(s):
  • Citation: P. Brommer, and F. Gähler (2006), "Effective potentials for quasicrystals fromab-initiodata", Philosophical Magazine, 86(6-8), 753-758. DOI: 10.1080/14786430500333349.
    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.

    Notes: This is for the Potential B model described in the reference

    Related Models:
  • See Computed Properties
    Notes: Listing found at https://openkim.org.
    Link(s):
 
 
 
 
 
  • Citation: O.R. Deluigi, R.C. Pasianot, F.J. Valencia, A. Caro, D. Farkas, and E.M. Bringa (2021), "Simulations of primary damage in a High Entropy Alloy: Probing enhanced radiation resistance", Acta Materialia, 213, 116951. DOI: 10.1016/j.actamat.2021.116951.
    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.

    Notes: This is a modified version of 2018--Farkas-D-Caro-A--Fe-Ni-Cr-Co-Cu that adds the ZBL correction at shorter interatomic distances making it suitable for radiation studies.

  • See Computed Properties
    Notes: This file was provided by Diana Farkas (Virginia Tech) on May 16, 2021 and posted with her permission.
    File(s):
 
 
  • Citation: R. Gröger, V. Vitek, and A. Dlouhý (2020), "Effective pair potential for random fcc CoCrFeMnNi alloys", Modelling and Simulation in Materials Science and Engineering, 28(7), 075006. DOI: 10.1088/1361-651x/ab7f8b.
    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 C12 − C44, 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.

    Notes: R. Gröger notes that "This is the Mie n-2n potential, where n=6 was found to give the best results - it is the same as the Lennard-Jones 6-12 potential. These potential files contain parameterizations of the Co-Cr-Fe-Mn-Ni system intended for studies of compositionally complex alloys with spatially random distributions of individual elements. Although it was developed primarily for studies of the quinary fcc CoCrFeMnNi system, the paper above demonstrates that it can be used equally well for quaternaries and ternaries. We emphasize that the model ceases to be applicable for binary and unary systems, where most or all first neighbor bonds are between the same elements."

  • See Computed Properties
    Notes: These files were provided by Roman Gröger on Oct 23, 2021 and posted with his permission. The file mie.mod shows how to include these potentials in LAMMPS simulations via "include mie.mod".
    File(s):
 
  • Citation: A. Liang, D.C. Goodelman, A.M. Hodge, D. Farkas, and P.S. Branicio (2023), "CoFeNiTix and CrFeNiTix high entropy alloy thin films microstructure formation", Acta Materialia, 257, 119163. DOI: 10.1016/j.actamat.2023.119163.
    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 CoFeNiTix and CrFeNiTix 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 CoFeNiTix and CrFeNiTix 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.

  • See Computed Properties
    Notes: This file was provided by Diana Farkas on January 12, 2024.
    File(s):
 
 
  • Citation: M.I. Mendelev (2024), "to be published".

    Notes: Mikhail Mendelev notes "This Ni-Co-Cr potential is designed to simulate the plastic deformation in the NiCoCr HEA. The potential correctly reproduces all element ground states and melting temperature as well as hcp-fcc transformation in pure Co. Note that the potential was fit to the PBEsol VASP data, therefore, it systematically underestimates the lattice parameters."

  • See Computed Properties
    Notes: This file was provided by Mikhail Mendelev on July 9, 2024.
    File(s):
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
  • Citation: L.K. Béland, C. Lu, Y.N. Osetskiy, G.D. Samolyuk, A. Caro, L. Wang, and R.E. Stoller (2016), "Features of primary damage by high energy displacement cascades in concentrated Ni-based alloys", Journal of Applied Physics, 119(8), 085901. DOI: 10.1063/1.4942533.
    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.

    Notes: Prof. Beland notes that "The potential takes elemental Ni from 2004--Mishin-Y--Ni-Al and Co from 2012--Purja-Pun-G-P-Mishin-Y--Co and mixes them. We first applied the effective gauge transformation, and then fitted the cross-term as to reproduce the heat of mixing of Ni(x)-Co(1-x). The potential is very soft at short distances. In order to perform collision cascades, it should be overlaid to the ZBL potential, with an outer cutoff of 2.0 Angstroms."

    Related Models:
  • LAMMPS pair_style eam/alloy (2016--Beland-L-K--Ni-Co--LAMMPS--ipr1)
    See Computed Properties
    Notes: This file was provided by Laurent Béland on 7 Nov 2019 and posted with his permission. Note: The EAM potential by itself is very soft at short distances. In order to perform collision cascades, use the hybrid style listed below.
    File(s):
  • LAMMPS pair_style hybrid/overlay zbl eam/alloy (2016--Beland-L-K--Ni-Co--LAMMPS--ipr2)
    See Computed Properties
    Notes: The eam file was provided by Laurent Béland on 7 Nov 2019 and posted with his permission. It is the same eam/alloy file as the above implementation. example.lammps.in provides an example of how to call the potential with the ZBL overlay applied.
    File(s):
 
 
 
 
 
Date Created: October 5, 2010 | Last updated: March 05, 2025