Jump to: Tests | Visualizers | Files | Wiki

Sim_LAMMPS_BOP_ZhouWardFoster_2015_CuH__SM_404135993060_000

Title
A single sentence description.
LAMMPS BOP potential for the Cu-H system developed by Zhou et al. (2015) v000
Description Despite extensive studies in the past, deterioration of mechanical properties due to hydrogen environment exposure remains a serious problem for structural materials. More effective improvement of a material's resilience requires advanced computational methods to elucidate the fundamental mechanisms of the hydrogen effects. To enable accurate molecular dynamics (MD) studies of the hydrogen effects on metals, we have developed a high-fidelity analytical bond-order potential (BOP) for the copper-hydrogen binary system as a representative case. This potential is available through the publically available MD code LAMMPS. The potential parameters are optimized using an iterative process. First, the potential is fitted to static and reactive properties of a variety of elemental and binary configurations including small clusters and bulk lattices (with coordination varying from 1 to 12). Then the potential is put through a series of rigorous MD simulation tests (e.g., vapor deposition and solidification) that involve chaotic initial configurations. It is demonstrated that this Cu-H BOP not only gives structural and property trends close to those seen in experiments and quantum mechanical calculations, but also predicts the correct phase transformations and chemical reactions in direct MD simulations. The correct structural evolution from chaotic initial states strongly verifies the transferability of the potential. A highly transferable potential is the reason that a well-parameterized analytical BOP can enable MD simulations of metal-hydrogen interactions to reach a fidelity level not achieved in the past.
Species
The supported atomic species.
Cu, H
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 ronmiller
Maintainer ronmiller
Author Ronald E. Miller
Publication Year 2019
Item Citation

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

[1] Zhou XW, Ward DK, Foster M, Zimmerman JA. An analytical bond-order potential for the copper–hydrogen binary system. Journal of Materials Science [Internet]. 2015Jan;50(7):2859–75. Available from: https://doi.org/10.1007/s10853-015-8848-9 doi:10.1007/s10853-015-8848-9

[2] Miller RE. LAMMPS BOP potential for the Cu-H system developed by Zhou et al. (2015) v000. OpenKIM; 2019. doi:10.25950/845f3cc8

[3] 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

[4] 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.
Short KIM ID
The unique KIM identifier code.
SM_404135993060_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_2015_CuH__SM_404135993060_000
DOI 10.25950/845f3cc8
https://doi.org/10.25950/845f3cc8
https://search.datacite.org/works/10.25950/845f3cc8
KIM Item TypeSimulator Model
KIM API Version2.1
Simulator Name
The name of the simulator as defined in kimspec.edn.
LAMMPS
Potential Type 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: Cu
Species: H

Click on any thumbnail to get a full size image.



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

Click on any thumbnail to get a full size image.



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

Click on any thumbnail to get a full size image.



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

Click on any thumbnail to get a full size image.



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

Click on any thumbnail to get a full size image.



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

Click on any thumbnail to get a full size image.



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

Click on any thumbnail to get a full size image.



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

Click on any thumbnail to get a full size image.



Cubic Crystal Basic Properties Table

Species: Cu

Species: H



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 2975
CohesiveEnergyVsLatticeConstant_diamond_Cu__TE_596332570306_003 view 2815
CohesiveEnergyVsLatticeConstant_fcc_Cu__TE_311348891940_003 view 3167
CohesiveEnergyVsLatticeConstant_sc_Cu__TE_767437873249_003 view 2815
ElasticConstantsCubic__TD_011862047401_005
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_H__TE_632848943374_005 view 21304
ElasticConstantsCubic_diamond_H__TE_744263053723_000 view 38695
ElasticConstantsCubic_fcc_H__TE_627409417266_005 view 27272
ElasticConstantsCubic_sc_H__TE_648951261715_005 view 21850
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 7837
ElasticConstantsCubic_diamond_Cu__TE_330878926469_001 view 28758
ElasticConstantsCubic_fcc_Cu__TE_188557531340_006 view 13243
ElasticConstantsCubic_sc_Cu__TE_319353354686_006 view 8957
ElasticConstantsHexagonal__TD_612503193866_003
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_003 view 9609
ElasticConstantsHexagonal_hcp_H__TE_758877215376_003 view 3934
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 6110
LatticeConstantCubicEnergy_bcc_H__TE_166241045885_007 view 21561
LatticeConstantCubicEnergy_diamond_Cu__TE_939141232476_007 view 16858
LatticeConstantCubicEnergy_diamond_H__TE_257661677950_007 view 23992
LatticeConstantCubicEnergy_fcc_Cu__TE_387272513402_007 view 12700
LatticeConstantCubicEnergy_fcc_H__TE_384479542888_007 view 15099
LatticeConstantCubicEnergy_sc_Cu__TE_904717264736_007 view 6302
LatticeConstantCubicEnergy_sc_H__TE_478794314457_007 view 12092
LatticeConstantHexagonalEnergy__TD_942334626465_004
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_004 view 91799
LatticeConstantHexagonalEnergy_hcp_H__TE_510142503120_004 view 250246
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 64490
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 283948887
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 573403




Wiki

Wiki is ready to accept new content.

Login to edit Wiki content