Adjusting the melting point of a model system via Gibbs-Duhem integration: application to a model of Aluminum

Abstract

Model interaction potentials for real materials are generally optimized with respect to only those experimental properties that are easily evaluated as mechanical averages (e.g., elastic constants (at T=0 K), static lattice energies and liquid structure). For such potentials, agreement with experiment for the non-mechanical properties, such as the melting point, is not guaranteed and such values can deviate significantly from experiment. We present a method for re-parameterizing any model interaction potential of a real material to adjust its melting temperature to a value that is closer to its experimental melting temperature. This is done without significantly affecting the mechanical properties for which the potential was modeled. This method is an application of Gibbs-Duhem integration [D. Kofke, Mol. Phys.78, 1331 (1993)]. As a test we apply the method to an embedded atom model of aluminum [J. Mei and J.W. Davenport, Phys. Rev. B 46, 21 (1992)] for which the melting temperature for the thermodynamic limit is 826.4 +/- 1.3K - somewhat below the experimental value of 933K. After re-parameterization, the melting temperature of the modified potential is found to be 931.5K +/- 1.5K.

9 pages, 5 figures, 4 tables

References in corpus (1)

• Symplectic algorithm for constant-pressure molecular dynamics using a Nose-Poincare thermostat

Cited by in corpus (11)

• Determination of phase diagrams via computer simulation: Methodology and applications to water, electrolytes and proteins
• Development of an Interatomic Potential for the Simulation of Defects, Plasticity and Phase Transformations in Titanium
• Entropy based fingerprint for local crystalline order
• Enhancing entropy and enthalpy fluctuations to drive crystallization in atomistic simulations
• Absence of superheating for ice Ih with a free surface : a new method of determining the melting point of different water models
• The atomistic representation of first strain-gradient elastic tensors
• Molecular Dynamics Simulation of Ga Penetration along Sigma 5 Symmetric Tilt Grain Boundaries in an Al Bicrystal
• Automated determination of grain boundary energy and potential-dependence using the OpenKIM framework
• Faceting transition in aluminum as a grain boundary phase transition
• Incorporating Gibbs free energy into interatomic potential fitting
• Generalized correction to embedded-atom potentials for modeling equilibrium and non-equilibrium properties of metals