• Citation: A. Nichol, and G.J. Ackland (2016), "Property trends in simple metals: An empirical potential approach", Physical Review B, 93(18), 184101. DOI: 10.1103/physrevb.93.184101.
    Abstract: We demonstrate that the melting points and other thermodynamic quantities of the alkali metals can be calculated based on static crystalline properties. To do this we derive analytic interatomic potentials for the alkali metals fitted precisely to cohesive and vacancy energies, elastic moduli, the lattice parameter, and crystal stability. These potentials are then used to calculate melting points by simulating the equilibration of solid and liquid samples in thermal contact at ambient pressure. With the exception of lithium, remarkably good agreement is found with experimental values. The instability of the bcc structure in Li and Na at low temperatures is also reproduced and, unusually, is not due to a soft T1N phonon mode. No forces or finite-temperature properties are included in the fit, so this demonstrates a surprisingly high level of intrinsic transferability in the simple potentials. Currently, there are few potentials available for the alkali metals, so in addition to demonstrating trends in behavior, we expect that the potentials will be of broad general use.

    Notes: G.J. Ackland noted that lattice parameters, elastic constants and cohesive energies were used in the fitting process, so the values produced by this conversion should match known values. He noted that bcc crystal structure should be stable and produce a melting temperature of 312 K. Publication information was updated on 12 Oct. 2017. Prior publication listing for this potential was Han, S., Zepeda-Ruiz, L. A., Ackland, G. J., Car, R., and Srolovitz, D. J. (2003). Interatomic potential for vanadium suitable for radiation damage simulations. Journal of Applied Physics, 93(6), 3328. DOI: 10.1063/1.1555275

    Related Models:
  • Moldy FS (2016--Nichol-A--Rb--MOLDY--ipr1)
    Notes: The parameters in Rb.moldy were obtained from http://homepages.ed.ac.uk/graeme/moldy/moldy.html and posted with the permission of G.J. Ackland.
    File(s):
  • LAMMPS pair_style eam/fs (2016--Nichol-A--Rb--LAMMPS--ipr1)
    See Computed Properties
    Notes: This conversion was performed by G.J. Ackland and submitted on 8 Dec. 2015.
    File(s): superseded


  • LAMMPS pair_style eam/fs (2016--Nichol-A--Rb--LAMMPS--ipr2)
    See Computed Properties
    Notes: A new conversion to LAMMPS performed by G.J. Ackland was submitted on 10 Oct. 2017. The previous setfl version above had a spurious oscillation period in the tabulated r*phi function that influenced measurements, most notably static elastic constant evaluations.
    File(s):
  • See Computed Properties
    Notes: Listing found at https://openkim.org. This KIM potential is based on the files from 2016--Nichol-A--Rb--LAMMPS--ipr2.
    Link(s):
Date Created: October 5, 2010 | Last updated: November 20, 2024