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

Title: Beyond Born-Mayer: Improved models for short-range repulsion in ab initio force fields

Abstract

Short-range repulsion within inter-molecular force fields is conventionally described by either Lennard-Jones or Born-Mayer forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of inter-molecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, and robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Lastly, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.

Authors:
 [1];  [2];  [3];  [1]
  1. Univ. of Wisconsin-Madison, Madison, WI (United States)
  2. Queen Mary Univ. of London, London (United Kingdom)
  3. Univ. of Cambridge, Cambridge (United Kingdom)
Publication Date:
Research Org.:
Univ. of Wisconsin-Madison, Madison, WI (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
1299438
Grant/Contract Number:
SC0014059
Resource Type:
Journal Article: Published Article
Journal Name:
Journal of Chemical Theory and Computation
Additional Journal Information:
Journal Name: Journal of Chemical Theory and Computation; Journal ID: ISSN 1549-9618
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
74 ATOMIC AND MOLECULAR PHYSICS

Citation Formats

Van Vleet, Mary J., Misquitta, Alston J., Stone, Anthony J., and Schmidt, Jordan R. Beyond Born-Mayer: Improved models for short-range repulsion in ab initio force fields. United States: N. p., 2016. Web. doi:10.1021/acs.jctc.6b00209.
Van Vleet, Mary J., Misquitta, Alston J., Stone, Anthony J., & Schmidt, Jordan R. Beyond Born-Mayer: Improved models for short-range repulsion in ab initio force fields. United States. doi:10.1021/acs.jctc.6b00209.
Van Vleet, Mary J., Misquitta, Alston J., Stone, Anthony J., and Schmidt, Jordan R. 2016. "Beyond Born-Mayer: Improved models for short-range repulsion in ab initio force fields". United States. doi:10.1021/acs.jctc.6b00209.
@article{osti_1299438,
title = {Beyond Born-Mayer: Improved models for short-range repulsion in ab initio force fields},
author = {Van Vleet, Mary J. and Misquitta, Alston J. and Stone, Anthony J. and Schmidt, Jordan R.},
abstractNote = {Short-range repulsion within inter-molecular force fields is conventionally described by either Lennard-Jones or Born-Mayer forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of inter-molecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, and robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Lastly, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.},
doi = {10.1021/acs.jctc.6b00209},
journal = {Journal of Chemical Theory and Computation},
number = ,
volume = ,
place = {United States},
year = 2016,
month = 6
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1021/acs.jctc.6b00209

Citation Metrics:
Cited by: 4works
Citation information provided by
Web of Science

Save / Share:
  • Short-range repulsion within inter-molecular force fields is conventionally described by either Lennard-Jones or Born-Mayer forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of inter-molecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, andmore » robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Lastly, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.« less
  • When a liquid is cooled well below its melting temperature at a rate that exceeds the critical cooling rate R{sub c}, the crystalline state is bypassed and a metastable, amorphous glassy state forms instead. R{sub c} (or the corresponding critical casting thickness d{sub c}) characterizes the glass-forming ability (GFA) of each material. While silica is an excellent glass-former with small R{sub c} < 10{sup −2} K/s, pure metals and most alloys are typically poor glass-formers with large R{sub c} > 10{sup 10} K/s. Only in the past thirty years have bulk metallic glasses (BMGs) been identified with R{sub c} approachingmore » that for silica. Recent simulations have shown that simple, hard-sphere models are able to identify the atomic size ratio and number fraction regime where BMGs exist with critical cooling rates more than 13 orders of magnitude smaller than those for pure metals. However, there are a number of other features of interatomic potentials beyond hard-core interactions. How do these other features affect the glass-forming ability of BMGs? In this manuscript, we perform molecular dynamics simulations to determine how variations in the softness and non-additivity of the repulsive core and form of the interatomic pair potential at intermediate distances affect the GFA of binary alloys. These variations in the interatomic pair potential allow us to introduce geometric frustration and change the crystal phases that compete with glass formation. We also investigate the effect of tuning the strength of the many-body interactions from zero to the full embedded atom model on the GFA for pure metals. We then employ the full embedded atom model for binary BMGs and show that hard-core interactions play the dominant role in setting the GFA of alloys, while other features of the interatomic potential only change the GFA by one to two orders of magnitude. Despite their perturbative effect, understanding the detailed form of the intermetallic potential is important for designing BMGs with cm or greater casting thickness.« less
  • We introduce new generalized (reverting to the original) and extended (not reverting to the original) 4-parameter forms of the (B-2) Potential Energy Function (PEF) of Wang etal. (L.-P. Wang, J. Chen and T. van Voorhis, J. Chem. Theor. Comp. 9, 452 (2013)), which is itself a modification of the Buckingham exponential-6 PEF. The new forms have a tunable, singularity-free short-range repulsion and an adjustable long-range attraction. They produce fits to high quality ab initio data for the X–(H2O), X=F, Cl, Br, I and M+(H2O), M=Li, Na, K, Rb, Cs dimers that are between 1 and 2 orders of magnitude bettermore » than the original 3-parameter (B-2) and modified Buckingham exponential-6 PEFs. They are also slightly better than the 4-parameter generalized Buckingham exponential-6(gBe-6) and of comparable quality with the 4-parameter extended Morse (eM) PEFs introduced recently by us.« less