Skip to content

Metal cluster (M2–M60, M = Au, Cu, Ni, Pt) formation as investigated using the reactive force field

Evidence and attribution

Authority of statements

Prose sections below (Summary, Methods, Findings, etc.) are curated summaries of the publication identified by doi, title, and pdf_path in the front matter above. They are not new primary claims by this wiki.

For definitive numerical values, reaction schemes, and interpretations, use the peer-reviewed article (and optional records under normalized/papers/ when present)—not this page alone.

Summary

Small-to-medium metal clusters M\(_N\) with N = 2–60 for Au, Cu, Ni, Pt are generated using ReaxFF within an atom-by-atom addition workflow at 298 K, automated by a wrapper program for high-throughput ReaxFF jobs. Structures and energies are compared to DFT and experimental references where available to validate the force fields and protocol. The study reports that ReaxFF is a reasonably accurate tool for structures and, to some extent, energies of bare clusters, while flagging discrepancies for very small (2–6 atom) Au and Pt subsets as motivation for further FF development. A key structural observation is that lowest-energy medium (20–60 atom) clusters under this growth protocol often do not reach highly symmetric geometries, plausibly because stable subunits from earlier steps resist rearrangement—so kinetic accessibility dominates over global symmetry minima.

Methods

1 — MD application (atomistic dynamics)

ReaxFF cluster generation treats M\(_2\)–M\(_{60}\) for M = Au, Cu, Ni, Pt using an atom-by-atom addition workflow at 298 K with Berendsen NVT thermostat control, automated by a wrapper program that iterates “add one atom → relax” (normalized/extracts/201215-33-29-tex-output-2-ctn03-2041_p1-2.txt + article abstract).

  • Engine / code: ReaxFF reactive modeling with a documented \(E_\mathrm{total}=E_\mathrm{valence}+E_\mathrm{Coulomb}+E_\mathrm{vdW}\) split and valence terms enumerated as in Eq. (2) (paper; extract references Chenoweth et al. SI). N/A — a specific MD engine (e.g., LAMMPS) is not named on the indexed excerpt pages.
  • System size & composition: Small-to-medium gas-phase clusters from dimers up to 60 atoms for each metal; gold includes an anionic series [Au\(_2\)]⁻–[Au\(_{60}\)]⁻ in addition to neutrals (abstract, extract).
  • Boundaries / periodicity: Non-periodic gas-phase clusters are implied by the cluster workflow (no bulk PBC claim on the indexed excerpt pages).
  • Ensemble: NVT with Berendsen thermostat (indexed excerpt).
  • Timestep / duration / stages: N/A — not stated on the indexed excerpt pages.
  • Thermostat: Berendsen (N/A — damping/time constant not recovered from the indexed excerpt pages—verify pdf_path).
  • Barostat: N/A — NVT protocol; no NPT control described on the indexed excerpt pages.
  • Temperature: 298 K (indexed excerpt).
  • Pressure / stress: N/A — not stated on the indexed excerpt pages.
  • Electric field: N/A — not stated on the indexed excerpt pages.
  • Replica / enhanced sampling: N/A — not stated on the indexed excerpt pages.

2 — Force-field training

N/A — this work applies published ReaxFF metal parameterizations and focuses on cluster construction + validation rather than refitting terms (abstract framing, extract).

3 — Static QM / DFT-only

DFT (and experiment where available) are used as references to assess cluster structures/energies (abstract); N/A — detailed DFT functional/basis/k-mesh tables are not on the indexed excerpt pages—verify pdf_path.

Analysis metrics: the paper defines energy per atom \(E_N/N\) and relative binding energy increments \(d(E_N/N)=E_N/N-E_{N+1}/(N+1)\) (Eq. (3), paper; referenced in wiki abstract).

Findings

The abstract positions this as the first reported sweep of ReaxFF cluster generation over N = 2–60 for these metals under the stated atom-by-atom protocol. Validation against DFT and experiment is described as generally good for structures and partially for energies, with explicit caveats motivating further FF work for very small (2–6 atom) Au and Pt subsets. For medium (20–60 atom) clusters, lowest-energy structures under this growth pathway often do not reach highly symmetric geometries; the authors attribute this to kinetic trapping by stable subunits formed in earlier steps and note that high-symmetry minima may not lie on the minimum-energy growth pathway.

Limitations

  • Atom-by-atom growth may miss global minima accessible from global optimization or statistical sampling of initial structures (explicitly noted in the article).
  • Filename in manifest encodes TeX timestamps; ignore as bibliographic noise—use JCTN citation above.

Relevance to group

ReaxFF metal parameterization and workflow tooling relevant to high-throughput nanoparticle exploration adjacent to group interests.

Citations and evidence anchors

  • DOI: 10.1166/jctn.2012.2041
  • Text-aligned pointer: normalized/extracts/201215-33-29-tex-output-2-ctn03-2041_p1-2.txt
  • Supporting structures: companion SI PDF paper:2011si-vgomzi-metalclust-venue-paper