Jump to: Tests | Visualizers | Files | Wiki

EAM_Dynamo_MendelevSordeletKramer_2007_CuZr__MO_120596890176_005

Interatomic potential for Copper (Cu), Zirconium (Zr).
Use this Potential

Title
A single sentence description.
Finnis-Sinclair potential (LAMMPS cubic hermite tabulation) for the Cu-Zr system developed by Mendelev, Sordelet and Kramer (2007) v005
Description
A short description of the Model describing its key features including for example: type of model (pair potential, 3-body potential, EAM, etc.), modeled species (Ac, Ag, ..., Zr), intended purpose, origin, and so on.
We propose a method of using atomistic computer simulations to obtain partial pair correlation functions from wide angle diffraction experiments with metallic liquids and their glasses. In this method, a model is first created using a semiempirical interatomic potential and then an additional atomic force is added to improve the agreement with experimental diffraction data. To illustrate this approach, the structure of an amorphous Cu64.5Zr35.5 alloy is highlighted, where we present the results for the semiempirical many-body potential and fitting to x-ray diffraction data. While only x-ray diffraction data were used in the present work, the method can be easily adapted to the case when there are also data from neutron diffraction or even in combination. Moreover, this method can be employed in the case of multicomponent systems when the data of several diffraction experiments can be combined. Note that the version published in [M.I. Mendelev, M.J. Kramer, R.T. Ott, D.J. Sordelet, D. Yagodin and P. Popel, Phil. Mag. 89, 967-987 (2009).] is better for simulation of liquid/glass structure.
Species
The supported atomic species.
Cu, Zr
Disclaimer
A statement of applicability provided by the contributor, informing users of the intended use of this KIM Item.
None
Content Origin http://www.ctcms.nist.gov/potentials/Cu.html
Contributor mendelev
Maintainer mendelev
Author Mikhail I. Mendelev
Publication Year 2018
Item Citation

This Model originally published in [1] is archived in OpenKIM [2-5].

[1] Mendelev MI, Sordelet DJ, Kramer MJ. Using atomistic computer simulations to analyze x-ray diffraction data from metallic glasses. Journal of Applied Physics. 2007;102(4):043501. doi:10.1063/1.2769157

[2] Mendelev MI. Finnis-Sinclair potential (LAMMPS cubic hermite tabulation) for the Cu-Zr system developed by Mendelev, Sordelet and Kramer (2007) v005. OpenKIM; 2018. doi:10.25950/8556dbaf

[3] Elliott RS. EAM Model Driver for tabulated potentials with cubic Hermite spline interpolation as used in LAMMPS v005. OpenKIM; 2018. doi:10.25950/68defa36

[4] Tadmor EB, Elliott RS, Sethna JP, Miller RE, Becker CA. The potential of atomistic simulations and the Knowledgebase of Interatomic Models. JOM. 2011;63(7):17. doi:10.1007/s11837-011-0102-6

[5] Elliott RS, Tadmor EB. Knowledgebase of Interatomic Models (KIM) Application Programming Interface (API). OpenKIM; 2011. doi:10.25950/ff8f563a

Click here to download the above citation in BibTeX format.
Funding Not available
Short KIM ID
The unique KIM identifier code.
MO_120596890176_005
Extended KIM ID
The long form of the KIM ID including a human readable prefix (100 characters max), two underscores, and the Short KIM ID. Extended KIM IDs can only contain alpha-numeric characters (letters and digits) and underscores and must begin with a letter.
EAM_Dynamo_MendelevSordeletKramer_2007_CuZr__MO_120596890176_005
DOI 10.25950/8556dbaf
https://doi.org/10.25950/8556dbaf
https://search.datacite.org/works/10.25950/8556dbaf
KIM Item Type
Specifies whether this is a Portable Model (software implementation of an interatomic model); Portable Model with parameter file (parameter file to be read in by a Model Driver); Model Driver (software implementation of an interatomic model that reads in parameters).
Portable Model using Model Driver EAM_Dynamo__MD_120291908751_005
DriverEAM_Dynamo__MD_120291908751_005
KIM API Version2.0
Potential Type eam
Programming Language(s)
The programming languages used in the code and the percentage of the code written in each one. "N/A" means "not applicable" and refers to model parameterizations which only include parameter tables and have no programming language.
N/A
Previous Version EAM_Dynamo_MendelevSordeletKramer_2007_CuZr__MO_120596890176_004

Verification Check Dashboard

(Click here to learn more about Verification Checks)

Grade Name Category Brief Description Full Results Aux File(s)
P vc-species-supported-as-stated mandatory
The model supports all species it claims to support; see full description.
Results Files
N/A vc-periodicity-support mandatory
Periodic boundary conditions are handled correctly; see full description.
Results Files
P vc-permutation-symmetry mandatory
Total energy and forces are unchanged when swapping atoms of the same species; see full description.
Results Files
B vc-forces-numerical-derivative consistency
Forces computed by the model agree with numerical derivatives of the energy; see full description.
Results Files
P vc-dimer-continuity-c1 informational
The energy versus separation relation of a pair of atoms is C1 continuous (i.e. the function and its first derivative are continuous); see full description.
Results Files
P vc-objectivity informational
Total energy is unchanged and forces transform correctly under rigid-body translation and rotation; see full description.
Results Files
P vc-inversion-symmetry informational
Total energy is unchanged and forces change sign when inverting a configuration through the origin; see full description.
Results Files
P vc-memory-leak informational
The model code does not have memory leaks (i.e. it releases all allocated memory at the end); see full description.
Results Files
P vc-thread-safe mandatory
The model returns the same energy and forces when computed in serial and when using parallel threads for a set of configurations. Note that this is not a guarantee of thread safety; see full description.
Results Files
P vc-unit-conversion mandatory
The model is able to correctly convert its energy and/or forces to different unit sets; see full description.
Results Files

Visualizers (in-page)


BCC Lattice Constant

This bar chart plot shows the mono-atomic body-centered cubic (bcc) lattice constant predicted by the current model (shown in the unique color) compared with the predictions for all other models in the OpenKIM Repository that support the species. The vertical bars show the average and standard deviation (one sigma) bounds for all model predictions. Graphs are generated for each species supported by the model.

Species: Cu
Species: Zr


Cohesive Energy Graph

This graph shows the cohesive energy versus volume-per-atom for the current mode for four mono-atomic cubic phases (body-centered cubic (bcc), face-centered cubic (fcc), simple cubic (sc), and diamond). The curve with the lowest minimum is the ground state of the crystal if stable. (The crystal structure is enforced in these calculations, so the phase may not be stable.) Graphs are generated for each species supported by the model.

Species: Zr
Species: Cu


Diamond Lattice Constant

This bar chart plot shows the mono-atomic face-centered diamond lattice constant predicted by the current model (shown in the unique color) compared with the predictions for all other models in the OpenKIM Repository that support the species. The vertical bars show the average and standard deviation (one sigma) bounds for all model predictions. Graphs are generated for each species supported by the model.

Species: Cu
Species: Zr


FCC Elastic Constants

This bar chart plot shows the mono-atomic face-centered cubic (fcc) elastic constants predicted by the current model (shown in blue) compared with the predictions for all other models in the OpenKIM Repository that support the species. The vertical bars show the average and standard deviation (one sigma) bounds for all model predictions. Graphs are generated for each species supported by the model.

Species: Zr
Species: Cu


FCC Lattice Constant

This bar chart plot shows the mono-atomic face-centered cubic (fcc) lattice constant predicted by the current model (shown in red) compared with the predictions for all other models in the OpenKIM Repository that support the species. The vertical bars show the average and standard deviation (one sigma) bounds for all model predictions. Graphs are generated for each species supported by the model.

Species: Zr
Species: Cu


FCC Stacking Fault Energies

This bar chart plot shows the intrinsic and extrinsic stacking fault energies as well as the unstable stacking and unstable twinning energies for face-centered cubic (fcc) predicted by the current model (shown in blue) compared with the predictions for all other models in the OpenKIM Repository that support the species. The vertical bars show the average and standard deviation (one sigma) bounds for all model predictions. Graphs are generated for each species supported by the model.

Species: Cu


FCC Surface Energies

This bar chart plot shows the mono-atomic face-centered cubic (fcc) relaxed surface energies predicted by the current model (shown in blue) compared with the predictions for all other models in the OpenKIM Repository that support the species. The vertical bars show the average and standard deviation (one sigma) bounds for all model predictions. Graphs are generated for each species supported by the model.

Species: Cu


SC Lattice Constant

This bar chart plot shows the mono-atomic simple cubic (sc) lattice constant predicted by the current model (shown in the unique color) compared with the predictions for all other models in the OpenKIM Repository that support the species. The vertical bars show the average and standard deviation (one sigma) bounds for all model predictions. Graphs are generated for each species supported by the model.

Species: Cu
Species: Zr


Cubic Crystal Basic Properties Table

Species: Cu

Species: Zr



Tests

CohesiveEnergyVsLatticeConstant__TD_554653289799_003
This Test Driver uses LAMMPS to compute the cohesive energy of a given monoatomic cubic lattice (fcc, bcc, sc, or diamond) at a variety of lattice spacings. The lattice spacings range from a_min (=a_min_frac*a_0) to a_max (=a_max_frac*a_0) where a_0, a_min_frac, and a_max_frac are read from stdin (a_0 is typically approximately equal to the equilibrium lattice constant). The precise scaling and number of lattice spacings sampled between a_min and a_0 (a_0 and a_max) is specified by two additional parameters passed from stdin: N_lower and samplespacing_lower (N_upper and samplespacing_upper). Please see README.txt for further details.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
CohesiveEnergyVsLatticeConstant_bcc_Cu__TE_864632638496_003 view 4223
CohesiveEnergyVsLatticeConstant_bcc_Zr__TE_783403151694_003 view 3935
CohesiveEnergyVsLatticeConstant_diamond_Cu__TE_596332570306_003 view 4319
CohesiveEnergyVsLatticeConstant_diamond_Zr__TE_742267498137_003 view 4574
CohesiveEnergyVsLatticeConstant_fcc_Cu__TE_311348891940_003 view 4255
CohesiveEnergyVsLatticeConstant_fcc_Zr__TE_241660333240_003 view 4255
CohesiveEnergyVsLatticeConstant_sc_Cu__TE_767437873249_003 view 4478
CohesiveEnergyVsLatticeConstant_sc_Zr__TE_943773936272_003 view 4382
ElasticConstantsCubic__TD_011862047401_006
Computes the cubic elastic constants for some common crystal types (fcc, bcc, sc, diamond) by calculating the hessian of the energy density with respect to strain. An estimate of the error associated with the numerical differentiation performed is reported.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
ElasticConstantsCubic_bcc_Cu__TE_091603841600_006 view 2047
ElasticConstantsCubic_bcc_Zr__TE_286034503723_006 view 6398
ElasticConstantsCubic_diamond_Cu__TE_330878926469_001 view 6366
ElasticConstantsCubic_diamond_Zr__TE_820660699174_001 view 11580
ElasticConstantsCubic_fcc_Cu__TE_188557531340_006 view 2367
ElasticConstantsCubic_fcc_Zr__TE_026250508553_006 view 2015
ElasticConstantsCubic_sc_Cu__TE_319353354686_006 view 1887
ElasticConstantsCubic_sc_Zr__TE_103738020637_006 view 6526
ElasticConstantsHexagonal__TD_612503193866_004
Computes the elastic constants for hcp crystals by calculating the hessian of the energy density with respect to strain. An estimate of the error associated with the numerical differentiation performed is reported.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
ElasticConstantsHexagonal_hcp_Cu__TE_198002759922_004 view 2260
ElasticConstantsHexagonal_hcp_Zr__TE_924023808720_004 view 1942
GrainBoundaryCubicCrystalSymmetricTiltRelaxedEnergyVsAngle__TD_410381120771_002
Computes grain boundary energy for a range of tilt angles given a crystal structure, tilt axis, and material.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
GrainBoundaryCubicCrystalSymmetricTiltRelaxedEnergyVsAngle_fcc100_Cu__TE_529988253259_000 view 5526549
GrainBoundaryCubicCrystalSymmetricTiltRelaxedEnergyVsAngle_fcc110_Cu__TE_708214008908_000 view 16787241
GrainBoundaryCubicCrystalSymmetricTiltRelaxedEnergyVsAngle_fcc111_Cu__TE_603516505525_000 view 8291766
GrainBoundaryCubicCrystalSymmetricTiltRelaxedEnergyVsAngle_fcc112_Cu__TE_288691353820_000 view 35002941
LatticeConstantCubicEnergy__TD_475411767977_007
Equilibrium lattice constant and cohesive energy of a cubic lattice at zero temperature and pressure.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
LatticeConstantCubicEnergy_bcc_Cu__TE_873531926707_007 view 2495
LatticeConstantCubicEnergy_bcc_Zr__TE_819253466839_007 view 2399
LatticeConstantCubicEnergy_diamond_Cu__TE_939141232476_007 view 3263
LatticeConstantCubicEnergy_diamond_Zr__TE_184605903050_007 view 4606
LatticeConstantCubicEnergy_fcc_Cu__TE_387272513402_007 view 4223
LatticeConstantCubicEnergy_fcc_Zr__TE_010442444476_007 view 2751
LatticeConstantCubicEnergy_sc_Cu__TE_904717264736_007 view 2335
LatticeConstantCubicEnergy_sc_Zr__TE_107850120912_007 view 2335
LatticeConstantHexagonalEnergy__TD_942334626465_005
Calculates lattice constant of hexagonal bulk structures at zero temperature and pressure by using simplex minimization to minimize the potential energy.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
LatticeConstantHexagonalEnergy_hcp_Cu__TE_344176839725_005 view 53389
LatticeConstantHexagonalEnergy_hcp_Zr__TE_888140777754_005 view 23177
LinearThermalExpansionCoeffCubic__TD_522633393614_001
This Test Driver uses LAMMPS to compute the linear thermal expansion coefficient at a finite temperature under a given pressure for a cubic lattice (fcc, bcc, sc, diamond) of a single given species.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
LinearThermalExpansionCoeffCubic_fcc_Cu__TE_335019190158_001 view 13562216
PhononDispersionCurve__TD_530195868545_004
Calculates the phonon dispersion relations for fcc lattices and records the results as curves.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
PhononDispersionCurve_fcc_Cu__TE_575177044018_004 view 52846
StackingFaultFccCrystal__TD_228501831190_002
Intrinsic and extrinsic stacking fault energies, unstable stacking fault energy, unstable twinning energy, stacking fault energy as a function of fractional displacement, and gamma surface for a monoatomic FCC lattice at zero temperature and pressure.
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
StackingFaultFccCrystal_0bar_Cu__TE_090810770014_002 view 13553803
SurfaceEnergyCubicCrystalBrokenBondFit__TD_955413365818_004
Calculates the surface energy of several high symmetry surfaces and produces a broken-bond model fit. In latex form, the fit equations are given by:

E_{FCC} (\vec{n}) = p_1 (4 \left( |x+y| + |x-y| + |x+z| + |x-z| + |z+y| +|z-y|\right)) + p_2 (8 \left( |x| + |y| + |z|\right)) + p_3 (2 ( |x+ 2y + z| + |x+2y-z| + |x-2y + z| + |x-2y-z| + |2x+y+z| + |2x+y-z| +|2x-y+z| +|2x-y-z| +|x+y+2z| +|x+y-2z| +|x-y+2z| +|x-y-2z| ) + c

E_{BCC} (\vec{n}) = p_1 (6 \left( | x+y+z| + |x+y-z| + |-x+y-z| + |x-y+z| \right)) + p_2 (8 \left( |x| + |y| + |z|\right)) + p_3 (4 \left( |x+y| + |x-y| + |x+z| + |x-z| + |z+y| +|z-y|\right)) +c.

In Python, these two fits take the following form:

def BrokenBondFCC(params, index):

import numpy
x, y, z = index
x = x / numpy.sqrt(x**2.+y**2.+z**2.)
y = y / numpy.sqrt(x**2.+y**2.+z**2.)
z = z / numpy.sqrt(x**2.+y**2.+z**2.)

return params[0]*4* (abs(x+y) + abs(x-y) + abs(x+z) + abs(x-z) + abs(z+y) + abs(z-y)) + params[1]*8*(abs(x) + abs(y) + abs(z)) + params[2]*(abs(x+2*y+z) + abs(x+2*y-z) +abs(x-2*y+z) +abs(x-2*y-z) + abs(2*x+y+z) +abs(2*x+y-z) +abs(2*x-y+z) +abs(2*x-y-z) + abs(x+y+2*z) +abs(x+y-2*z) +abs(x-y+2*z) +abs(x-y-2*z))+params[3]

def BrokenBondBCC(params, x, y, z):


import numpy
x, y, z = index
x = x / numpy.sqrt(x**2.+y**2.+z**2.)
y = y / numpy.sqrt(x**2.+y**2.+z**2.)
z = z / numpy.sqrt(x**2.+y**2.+z**2.)

return params[0]*6*(abs(x+y+z) + abs(x-y-z) + abs(x-y+z) + abs(x+y-z)) + params[1]*8*(abs(x) + abs(y) + abs(z)) + params[2]*4* (abs(x+y) + abs(x-y) + abs(x+z) + abs(x-z) + abs(z+y) + abs(z-y)) + params[3]
Test Test Results Link to Test Results page Benchmark time
Usertime muliplied by the Whetstone Benchmark. This number can be used (approximately) to compare the performance of different models independently of the architecture on which the test was run.

Measured in Millions of Whetstone Instructions (MWI)
SurfaceEnergyCubicCrystalBrokenBondFit_fcc_Cu__TE_689904280697_004 view 39538





Download Dependency

This Model requires a Model Driver. Archives for the Model Driver EAM_Dynamo__MD_120291908751_005 appear below.


EAM_Dynamo__MD_120291908751_005.txz Tar+XZ Linux and OS X archive
EAM_Dynamo__MD_120291908751_005.zip Zip Windows archive

Wiki

Wiki is ready to accept new content.

Login to edit Wiki content