Source author record

Anthony J. Stone

Anthony J. Stone appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2016arXiv

Ab initio atom-atom potentials using CamCASP: Many-body potentials for the pyridine dimer

In Part I of this two-part investigation we described a methodology for the development of robust, analytic, many-body atom-atom potentials for small organic molecules from first principles and demonstrated how the CamCASP program can be used to derive the damped, distributed multipole models for pyridine. Here we demonstrate how the theoretical ideas for the short-range models described in Part I, which are implemented in the CamCASP suite of programs, can be used to develop a series of many-body potentials for the pyridine system. Even the simplest of these potentials exhibit r.m.s. errors of only about 0.6kJ mol-1 for the low-energy pyridine dimers, significantly surpassing the best empirical potentials. Our best model is shown to support eight stable minima, four of which have not been reported in the literature before. Further, the functional form can be made systematically more elaborate so as to improve the accuracy without a significant increase in the human-time spent in their generation. We investigate the effects of anisotropy, rank of multipoles, and choice of polarizability and dispersion models.

preprint2016arXiv

Ab initio atom-atom potentials using CamCASP: Theory and application to multipole models for the pyridine dimer

Creating accurate, analytic atom--atom potentials for small organic molecules from first principles can be a time-consuming and computationally intensive task, particularly if we also require them to include explicit polarization terms, which are essential in many systems. In this first part of a two-part investigation, we describe how the {\sc CamCASP} suite of programs can be used to generate such potentials using some of the most accurate electronic structure methods currently applicable. In particular, we derive the long-range terms from monomer properties, and determine the short-range anisotropy parameters by a novel and robust method based on the iterated stockholder atom approach. We use the techniques described here to develop distributed multipole models for the electrostatic, polarization and dispersion interactions in the pyridine dimer. In the second part of this work we will apply these methods to develop a series of many-body potentials for the pyridine system.

preprint2016arXiv

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

Short-range repulsion within inter-molecular force fields is conventionally described by either Lennard-Jones (${A}/{r^{12}}$) or Born-Mayer ($A\exp(-Br)$) 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. Finally, 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.

preprint2014arXiv

Anomalous non-additive dispersion interactions in systems of three one-dimensional wires

The non-additive dispersion contribution to the binding energy of three one-dimensional (1D) wires is investigated using wires modelled by (i) chains of hydrogen atoms and (ii) homogeneous electron gases. We demonstrate that the non-additive dispersion contribution to the binding energy is significantly enhanced compared with that expected from Axilrod-Teller-Muto-type triple-dipole summations and follows a different power-law decay with separation. The triwire non-additive dispersion for 1D electron gases scales according to the power law $d^{-β}$, where $d$ is the wire separation, with exponents $β(r_s)$ smaller than 3 and slightly increasing with $r_s$ from 2.4 at $r_s = 1$ to 2.9 at $r_s=10$, where $r_s$ is the density parameter of the 1D electron gas. This is in good agreement with the exponent $β=3$ suggested by the leading-order charge-flow contribution to the triwire non-additivity, and is a significantly slower decay than the $\sim d^{-7}$ behaviour that would be expected from triple-dipole summations.

preprint2014arXiv

Distributed multipoles from a robust basis-space implementation of the iterated stockholder atoms procedure

The recently developed iterated stockholder atoms (ISA) approach of Lillestolen and Wheatley (Chem. Commun. {\bf 2008}, 5909 (2008)) offers a powerful method for defining atoms in a molecule. However, the real-space algorithm is known to converge very slowly, if at all. Here we present a robust, basis-space algorithm of the ISA method and demonstrate its applicability on a variety of systems. We show that this algorithm exhibits rapid convergence (taking around 10--80 iterations) with the number of iterations needed being unrelated to the system size or basis set used. Further, we show that the multipole moments calculated using this basis-space ISA method are as good as, or better than those obtained from Stone's distributed multipole analysis (J. Chem. Theory Comput. {\bf 1}, 1128 (2005) ), exhibiting better convergence properties and resulting in better behaved penetration energies. This can have significant consequences in the development of intermolecular interaction models.

preprint2010arXiv

Dispersion interactions between semiconducting wires

The dispersion energy between extended molecular chains (or equivalently infinite wires) with non-zero band gaps is generally assumed to be expressible as a pair-wise sum of atom-atom terms which decay as $R^{-6}$. Using a model system of two parallel wires with a variable band gap, we show that this is not the case. The dispersion interaction scales as $z^{-5}$ for large interwire separations $z$, as expected for an insulator, but as the band gap decreases the interaction is greatly enhanced; while at shorter (but non-overlapping) separations it approaches a power-law scaling given by $z^{-2}$, \emph{i.e.} the dispersion interaction expected between \emph{metallic} wires. We demonstrate that these effects can be understood from the increasing length scale of the plasmon modes (charge fluctuations), and their increasing contribution to the molecular dipole polarizability and the dispersion interaction, as the band gaps are reduced. This result calls into question methods which invoke locality assumptions in deriving dispersion interactions between extended small-gap systems.