skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Understanding Chemical Bonding in Alloys and the Representation in Atomistic Simulations

Journal Article · · Journal of Physical Chemistry. C
 [1];  [2];  [3]; ORCiD logo [4]; ORCiD logo [2]; ORCiD logo [1]
  1. Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
  2. Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, United States
  3. Univ. of California, Los Angeles, CA (United States). Department of Physics and Astronomy
  4. Univ. of California, Los Angeles, CA (United States). Department of Materials Science and Engineering

Alloys are widely used in catalysts and structural materials. The nature of chemical bonding and the origin of alloy formation energies, defect energies, and interfacial properties have not been well understood to date but are critical to material performance. In this contribution, we explain the polar nature of chemical bonding and an implementation in classical and reactive atomistic simulations to understand such properties more quantitatively. Electronegativity differences between metal atoms lead to polar bonding, and exothermic alloy formation energies are related to charge transfer between the different elements. These differences can be quantified by atomic charges using pairwise charge increments, determined by matching the computed alloy formation energy to experimentally measured alloy formation energies using pair potentials for the pure metals. The polar character of alloys is comparable to organic molecules and partially ionic minerals, for example, AlNi and AlNi3 alloys assume significant atomic charges of ±0.40e and +0.60e/–0.20e, respectively. The subsequent analysis of defect sites and defect energies using force-field-based calculations shows excellent agreement with calculations using density functional theory and embedded atom models (EAM). The formation of vacancy and antisite defects is characterized by a redistribution of charge in the first shell of neighbor atoms in the classical models whereby electroneutrality is maintained and charge increments correlate with differences in electronegativity. The proposed atomic charges represent internal dipole and multipole moments, consistent with existing definitions for organic and inorganic compounds and with the extended Born model (Heinz, H.; Suter, U. W. J. Phys. Chem. B2004,108 (47), 18341–18352). The method can be applied to any alloy and has a reproducibility of ±10%. In contrast, quantum mechanical charge schemes remain associated with deviations exceeding ±100%. The atomic charges for alloys provide a simple initial measure for the internal electronic structure, surface adsorption of molecules, and reactivity in catalysis and corrosion. The models are compatible with the Interface force field (IFF), CHARMM, AMBER, OPLS-AA, PCFF, CVFF, and GROMOS for reliable atomistic simulations of alloys and their interfaces with minerals and electrolytes from the nanometer scale to the micrometer scale.

Research Organization:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Argonne National Lab. (ANL), Argonne, IL (United States)
Sponsoring Organization:
USDOE Office of Science (SC)
DOE Contract Number:
AC02-06CH11357
OSTI ID:
1565666
Journal Information:
Journal of Physical Chemistry. C, Vol. 122, Issue 26; ISSN 1932-7447
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English

Similar Records

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials
Journal Article · Fri Jan 29 00:00:00 EST 2021 · npj Computational Materials · OSTI ID:1565666

Interpretable molecular models for molybdenum disulfide and insight into selective peptide recognition
Journal Article · Wed Aug 26 00:00:00 EDT 2020 · Chemical Science · OSTI ID:1565666

Force Field for Mica-Type Silicates and Dynamics of Octadecylammonium Chains Grafted to Montmorillonite
Journal Article · Sat Jan 01 00:00:00 EST 2005 · Chem. Mater. · OSTI ID:1565666