The convergence of the ab-initio many-body expansion for the cohesive energy of solid mercury

A many-body expansion for mercury clusters of the form E = \sum_{i<j}\Delta \epsilon_{ij} + \sum_{i<j<k}\Delta \epsilon_{ijk} + ... \quad, does not converge smoothly with increasing cluster size towards the solid state. Even for smaller cluster sizes (up to n=6), where van der Waals forces still dominate, one observes bad convergence behaviour. For solid mercury the convergence of the many-body expansion can dramatically be improved by an incremental procedure within an embedded cluster approach. Here one adds the coupled cluster many-body electron correlation contributions of the embedded cluster to the bulk HF energy. In this way we obtain a cohesive energy (not corrected for zero-point vibration) of 0.79 eV in perfect agreement with the experimental value.Comment: 10 pages, 3 figures, accepted PR

A simple environment-dependent overlap potential and Cauchy violation in solid argon

We develop an analytic and environment-dependent interatomic potential for the overlap repulsion in solid argon, based on an approximate treatment of the non-orthogonal Tight-Binding theory for the closed-shell systems. The present model can well reproduce the observed elastic properties of solid argon including Cauchy violation at high pressures, yet very simple. A useful and novel analysis is given to show how the elastic properties are related to the environment-dependence incorporated into a generic pairwise potential. The present study has a close link to the broad field of computational materials science, in which the inclusion of environment dependence in short-ranged repulsive part of a potential model is sometimes crucial in predicting the elastic properties correctly.Comment: 10 pages, 3 figure

First Principles Calculation of Elastic Properties of Solid Argon at High Pressures

The density and the elastic stiffness coefficients of fcc solid argon at high pressures from 1 GPa up to 80 GPa are computed by first-principles pseudopotential method with plane-wave basis set and the generalized gradient approximation (GGA). The result is in good agreement with the experimental result recently obtained with the Brillouin spectroscopy by Shimizu et al. [Phys. Rev. Lett. 86, 4568 (2001)]. The Cauchy condition was found to be strongly violated as in the experimental result, indicating large contribution from non-central many-body force. The present result has made it clear that the standard density functional method with periodic boundary conditions can be successfully applied for calculating elastic properties of rare gas solids at high pressures in contrast to those at low pressures where dispersion forces are important.Comment: 4 pages, 5 figures, submitted to PR

First-principles calculations of the phase stability of TiO2

Published versio

Influence of three-body forces and anharmonic effects on the zero-point energy of rare-gas crystals

The stabilization of the fcc structure of the heavier rare-gas crystals is mainly due to the zero-point energy (ZPE) calculated at the harmonic level with two-body contributions only. We evaluate the influence of the anharmonic contributions on the ZPE within the Einstein approximation. For the influence of the three-body contributions we develop an analytic three- body potential fitted to coupled-cluster ab initio data. Both the anharmonic and three-body contributions change the absolute value of the ZPE, but have no influence on the fcc-hcp energy difference