Non-Covalent Interactions in Molecular Crystals: Exploring the Accuracy of the Exchange-Hole Dipole Moment Model with Local Orbitals

• Published in: Journal of chemical theory and computation 2018-11, Vol.14 (11), p.5715-5724
• Main Authors: LeBlanc, Luc M, Weatherby, Joseph A, Otero-de-la-Roza, Alberto, Johnson, Erin R
• Format: Article
• Language: English
• Subjects: Accuracy; Crystal lattices; Dimers; Dipole moments; Error correction; Exchanging; Geometric accuracy; Mathematical models; Model accuracy; Orbitals; Phosphorene; Tetragonal zirconia polycrystals
• ISSN: 1549-9618; 1549-9626
• DOI: 10.1021/acs.jctc.8b00797

We present the first implementation of the exchange-hole dipole moment (XDM) model in combination with a numerical finite-support local orbital method (the SIESTA method) for the modeling of non-covalent interactions in periodic solids. The XDM model is parametrized for both the B86bPBE and PBE functionals using double-ζ- and triple-ζ-quality basis sets (DZP and TZP). The use of finite-support local orbitals is shown to have minimal impact on the computed dispersion coefficients for van der Waals molecular dimers and small molecular solids. However, the quality of the basis set affects the accuracy of calculated dimer binding energies and molecular-crystal lattice energies quite significantly; the size of the counterpoise correction indicates that this is caused by basis-set incompleteness error. In the case of the DZP basis set, its performance for weakly bound gas-phase dimers is similar to that of a double-ζ Gaussian basis set without diffuse functions. The new XDM implementation was tested on graphite and phosphorene exfoliation, and on the X23 benchmark set of molecular-crystal lattice energies. Our results indicate that lattice energies similar to plane-wave calculations can be obtained only if the counterpoise correction is applied. Alternatively, the calculated equilibrium geometries are reasonably close to the plane-wave equivalents, and composite approaches in which a single-point plane-wave calculation is used at the XDM/DZP equilibrium geometry yield good accuracy at a significantly lower computational cost.