English

Testing Lennard-Jones Clusters for Optimality

Atomic and Molecular Clusters 2024-09-04 v3 Chemical Physics

Abstract

This note advertises a simple necessary condition for optimality that any list Nvx(N)N\mapsto v^x(N) of computer-generated putative lowest average pair energies vx(N)v^x(N) of clusters that consist of NN monomers has to satisfy, whenever the monomers interact with each other through pair forces satisfying Newton's ``actio equals re-actio.'' These can be quite complicated, as for instance in the TIP5P model with five-site potential for a rigid tetrahedral-shaped H2_2O monomer of water, or as simple as the Lennard-Jones single-site potential for the center of an atomic monomer (which is also used for one site of the H2_2O monomer in the TIP5P model, that in addition has four peripheral sites with Coulomb potentials). The empirical usefulness of the necessary condition is demonstrated by testing a list of publicly available Lennard-Jones cluster data that have been pooled from 17 sources, covering the interval 2N16102\leq N\leq 1610 without gaps. The data point for N=447N=447 failed this test, meaning the listed 447447-particle Lennard-Jones cluster energy was not optimal. To implement this test for optimality in search algorithms for putatively optimal configurations is an easy task. Publishing only the data that pass the test would increase the odds that these are actually optimal, without guaranteeing it, though.

Keywords

Cite

@article{arxiv.2305.10600,
  title  = {Testing Lennard-Jones Clusters for Optimality},
  author = {Michael K. -H. Kiessling},
  journal= {arXiv preprint arXiv:2305.10600},
  year   = {2024}
}

Comments

12 pages, final revised and enlarged version based on comments received from referee. Accepted for publication in J. Chem. Phys

R2 v1 2026-06-28T10:37:40.847Z