Jump to: Tests | Visualizers | Files | Wiki

Sim_LAMMPS_BOP_ZhouWardFoster_2016_AlCu__SM_566399258279_000

Interatomic potential for Aluminum (Al), Copper (Cu).
Use this Potential

Title
A single sentence description.
 LAMMPS BOP potential for the Al-Cu system developed by Zhou, Ward, and Foster (2016) v000
Description Al-rich Al1-xCux alloys are important structural materials in the aerospace industry due to their high strength to density ratio. They are also emerging materials for hydrogen containing structures due to their potentially high resistance to hydrogen embrittlement. To enable accurate simulations of the mechanical behavior of Al1-xCux alloys that can guide material improvement, we have developed a high-fidelity analytical bond-order potential (BOP) for the Al-Cu system (the code is publically available in molecular dynamics package LAMMPS). The formalism of the potential is derived from quantum mechanical theories, and the parameters are optimized in an iteration fashion. The iterations begin by fitting properties of a variety of elemental and compound configurations (with coordination varying from 1 to 12) including small clusters, bulk lattices, defects, and surfaces. Following the fitting process, crystalline growth of important equilibrium phases is checked through molecular dynamics simulations of vapor deposition. It is demonstrated that this Al-Cu bond-order potential has unique advantages relative to existing literature potentials in reproducing structural and property tends from experiments and quantum-mechanical calculations, and providing good descriptions of melting temperature, defect characteristics, and surface energies. Most importantly, this BOP is the only potential currently available capable of capturing the Al-rich end of the Al-Cu phase diagram. This capability is rigorously verified by the potential's ability to capture the crystalline growth of the ground-state structures for elemental Al and Cu, as well as, the theta and theta' phases of the Al2Cu compound in vapor deposition simulations.
Species
The supported atomic species.
Al, Cu
Disclaimer
A statement of applicability provided by the contributor, informing users of the intended use of this KIM Item.
None
Content Origin LAMMPS package 22-Sep-2017
Contributor Ronald E. Miller
Maintainer Ronald E. Miller
Published on KIM 2019
How to Cite Click here to download this citation in BibTeX format.
Funding Not available
Short KIM ID
The unique KIM identifier code.
SM_566399258279_000
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.
Sim_LAMMPS_BOP_ZhouWardFoster_2016_AlCu__SM_566399258279_000
DOI 10.25950/9e693ed2
https://doi.org/10.25950/9e693ed2
https://commons.datacite.org/doi.org/10.25950/9e693ed2
KIM Item TypeSimulator Model
KIM API Version2.1
Simulator Name
The name of the simulator as defined in kimspec.edn.
LAMMPS
Potential Type bop
Simulator Potential bop

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
F 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
F 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
N/A 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

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: Al
Species: Cu


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: Cu
Species: Al


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: Al


Dislocation Core Energies

This graph shows the dislocation core energy of a cubic crystal at zero temperature and pressure for a specific set of dislocation core cutoff radii. After obtaining the total energy of the system from conjugate gradient minimizations, non-singular, isotropic and anisotropic elasticity are applied to obtain the dislocation core energy for each of these supercells with different dipole distances. Graphs are generated for each species supported by the model.

(No matching species)

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: Al
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: Cu
Species: Al


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: Al
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: Al
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: Al


Cubic Crystal Basic Properties Table

Species: Al

Species: Cu



Tests


CohesiveEnergyVsLatticeConstant__TD_554653289799_003
Cohesive energy versus lattice constant curve for monoatomic cubic lattices v003

Creators:
Contributor: karls
Publication Year: 2019
DOI: https://doi.org/10.25950/64cb38c5

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 multiplied 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)
Cohesive energy versus lattice constant curve for bcc Al v003 view 2687
Cohesive energy versus lattice constant curve for bcc Cu v003 view 3423
Cohesive energy versus lattice constant curve for diamond Al v003 view 2623
Cohesive energy versus lattice constant curve for diamond Cu v003 view 2879
Cohesive energy versus lattice constant curve for fcc Al v003 view 2751
Cohesive energy versus lattice constant curve for fcc Cu v003 view 3359
Cohesive energy versus lattice constant curve for sc Al v003 view 2687
Cohesive energy versus lattice constant curve for sc Cu v003 view 2847

ElasticConstantsCubic__TD_011862047401_006
Elastic constants for cubic crystals at zero temperature and pressure v006

Creators: Junhao Li and Ellad Tadmor
Contributor: tadmor
Publication Year: 2019
DOI: https://doi.org/10.25950/5853fb8f

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 multiplied 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)
Elastic constants for bcc Al at zero temperature v006 view 5566
Elastic constants for bcc Cu at zero temperature v006 view 10652
Elastic constants for diamond Al at zero temperature v001 view 41426
Elastic constants for diamond Cu at zero temperature v001 view 29110
Elastic constants for fcc Al at zero temperature v006 view 12028
Elastic constants for fcc Cu at zero temperature v006 view 16090
Elastic constants for sc Al at zero temperature v006 view 4702
Elastic constants for sc Cu at zero temperature v006 view 9373

ElasticConstantsHexagonal__TD_612503193866_004
Elastic constants for hexagonal crystals at zero temperature v004

Creators: Junhao Li
Contributor: jl2922
Publication Year: 2019
DOI: https://doi.org/10.25950/d794c746

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 multiplied 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)
Elastic constants for hcp Al at zero temperature v004 view 5030
Elastic constants for hcp Cu at zero temperature v004 view 9041

GrainBoundaryCubicCrystalSymmetricTiltRelaxedEnergyVsAngle__TD_410381120771_002
Relaxed energy as a function of tilt angle for a symmetric tilt grain boundary within a cubic crystal v002

Creators: Brandon Runnels
Contributor: brunnels
Publication Year: 2019
DOI: https://doi.org/10.25950/4723cee7

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 multiplied 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)
Relaxed energy as a function of tilt angle for a 100 symmetric tilt grain boundary in fcc Al v002 view 42185062
Relaxed energy as a function of tilt angle for a 100 symmetric tilt grain boundary in fcc Cu v000 view 71448848
Relaxed energy as a function of tilt angle for a 110 symmetric tilt grain boundary in fcc Al v000 view 129792702
Relaxed energy as a function of tilt angle for a 110 symmetric tilt grain boundary in fcc Cu v000 view 208566786
Relaxed energy as a function of tilt angle for a 111 symmetric tilt grain boundary in fcc Al v000 view 78561917
Relaxed energy as a function of tilt angle for a 111 symmetric tilt grain boundary in fcc Cu v000 view 109895708
Relaxed energy as a function of tilt angle for a 112 symmetric tilt grain boundary in fcc Al v000 view 310905542
Relaxed energy as a function of tilt angle for a 112 symmetric tilt grain boundary in fcc Cu v000 view 428574223

LatticeConstantCubicEnergy__TD_475411767977_007
Equilibrium lattice constant and cohesive energy of a cubic lattice at zero temperature and pressure v007

Creators: Daniel S. Karls and Junhao Li
Contributor: karls
Publication Year: 2019
DOI: https://doi.org/10.25950/2765e3bf

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 multiplied 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)
Equilibrium zero-temperature lattice constant for bcc Al v007 view 3615
Equilibrium zero-temperature lattice constant for bcc Cu v007 view 6206
Equilibrium zero-temperature lattice constant for diamond Al v007 view 17818
Equilibrium zero-temperature lattice constant for diamond Cu v007 view 17914
Equilibrium zero-temperature lattice constant for fcc Al v007 view 13052
Equilibrium zero-temperature lattice constant for fcc Cu v007 view 14267
Equilibrium zero-temperature lattice constant for sc Al v007 view 6526
Equilibrium zero-temperature lattice constant for sc Cu v007 view 6302

LatticeConstantHexagonalEnergy__TD_942334626465_005
Equilibrium lattice constants for hexagonal bulk structures at zero temperature and pressure v005

Creators: Daniel S. Karls and Junhao Li
Contributor: karls
Publication Year: 2019
DOI: https://doi.org/10.25950/c339ca32

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 multiplied 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)
Equilibrium lattice constants for hcp Al v005 view 3961290
Equilibrium lattice constants for hcp Cu v005 view 3773266

LinearThermalExpansionCoeffCubic__TD_522633393614_001
Linear thermal expansion coefficient of cubic crystal structures v001

Creators: Mingjian Wen
Contributor: mjwen
Publication Year: 2019
DOI: https://doi.org/10.25950/fc69d82d

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 multiplied 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)
Linear thermal expansion coefficient of fcc Al at 293.15 K under a pressure of 0 MPa v001 view 104465294
Linear thermal expansion coefficient of fcc Cu at 293.15 K under a pressure of 0 MPa v001 view 238052285

PhononDispersionCurve__TD_530195868545_004
Phonon dispersion relations for an fcc lattice v004

Creators: Matt Bierbaum
Contributor: mattbierbaum
Publication Year: 2019
DOI: https://doi.org/10.25950/64f4999b

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 multiplied 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)
Phonon dispersion relations for fcc Al v004 view 62411
Phonon dispersion relations for fcc Cu v004 view 62411

StackingFaultFccCrystal__TD_228501831190_002
Stacking and twinning fault energies of an fcc lattice at zero temperature and pressure v002

Creators:
Contributor: SubrahmanyamPattamatta
Publication Year: 2019
DOI: https://doi.org/10.25950/b4cfaf9a

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 multiplied 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)
Stacking and twinning fault energies for fcc Al v002 view 126724871
Stacking and twinning fault energies for fcc Cu v002 view 330022708

SurfaceEnergyCubicCrystalBrokenBondFit__TD_955413365818_004
High-symmetry surface energies in cubic lattices and broken bond model v004

Creators: Matt Bierbaum
Contributor: mattbierbaum
Publication Year: 2019
DOI: https://doi.org/10.25950/6c43a4e6

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 multiplied 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)
Broken-bond fit of high-symmetry surface energies in fcc Al v004 view 345321
Broken-bond fit of high-symmetry surface energies in fcc Cu v004 view 612941


Errors

StackingFaultFccCrystal__TD_228501831190_001

No Driver
Verification Check Error Categories Link to Error page
ForcesNumerDeriv__VC_710586816390_002 other view


Files

AlCu.bop.table
CMakeLists.txt
LICENSE
kimprovenance.edn
kimspec.edn
smspec.edn

Download

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

Wiki

Wiki is ready to accept new content.
Login to edit Wiki content