Rapid generation of optimal generalized Monkhorst-Pack grids

Abstract Computational modeling of the properties of crystalline materials has become an increasingly important aspect of materials research, consuming hundreds of millions of CPU-hours at scientific computing centers around the world each year, if not more. A routine operation in such calculations is the evaluation of integrals over the Brillouin zone. We have previously demonstrated that performing such integrals using generalized Monkhorst-Pack k-point grids can roughly double the speed of these calculations relative to the widely-used traditional Monkhorst-Pack grids. However the generation of optimal generalized Monkhorst-Pack grids is not implemented in most software packages due to the computational cost and difficulty of identifying the best grids. To address this problem, we present new algorithms that allow rapid generation of optimal generalized Monkhorst-Pack grids on the fly. We demonstrate that the grids generated by these algorithms are on average significantly more efficient than those generated using existing algorithms across a range of grid densities. For grids that correspond to a real-space supercell with at least 50 A between lattice points, which is sufficient to converge density functional theory calculations within 1 meV/atom for nearly all materials, our algorithm finds optimized grids in an average of 0.19 s on a single processing core. To facilitate the widespread adoption of this approach, we present new open-source tools including a library designed for integration with third-party software packages.

[1]  H. Monkhorst,et al.  SPECIAL POINTS FOR BRILLOUIN-ZONE INTEGRATIONS , 1976 .

[2]  Tim Mueller,et al.  Mechanisms for hydrogen evolution on transition metal phosphide catalysts and a comparison to Pt(111). , 2019, Physical chemistry chemical physics : PCCP.

[4]  Gus L. W. Hart,et al.  Generating derivative structures: Algorithm and applications , 2008, 0804.3544.

[5]  Tim Mueller,et al.  Computationally generated maps of surface structures and catalytic activities for alloy phase diagrams , 2019, Proceedings of the National Academy of Sciences.

[6]  Timothy Keith Mueller.  Computational studies of hydrogen storage materials and the development of related methods , 2007 .

[7]  Soler,et al.  Optimal meshes for integrals in real- and reciprocal-space unit cells. , 1992, Physical review. B, Condensed matter.

[8]  Gus L. W. Hart,et al.  Efficiency of Generalized Regular k-point grids , 2018, Computational Materials Science.

[9]  Carmelo Giacovazzo,et al.  Fundamentals of Crystallography , 2002 .

[10]  Xiaoqing Pan,et al.  Differential Surface Elemental Distribution Leads to Significantly Enhanced Stability of PtNi-Based ORR Catalysts , 2019 .

[11]  Gus L. W. Hart,et al.  Generalized regular k-point grid generation on the fly , 2019, Computational Materials Science.

[12]  K. Livi,et al.  Low-Overpotential Electroreduction of Carbon Monoxide Using Copper Nanowires , 2017 .

[13]  Tim Mueller,et al.  Efficient generation of generalized Monkhorst-Pack grids through the use of informatics , 2016 .

[14]  Kamal Choudhary,et al.  Convergence and machine learning predictions of Monkhorst-Pack k-points and plane-wave cut-off in high-throughput DFT calculations. , 2018, Computational materials science.

[15]  T. Puig,et al.  Local strain-driven migration of oxygen vacancies to apical sites in YBa2Cu3O7-x. , 2020, Nanoscale.

[16]  Conrad W. Rosenbrock,et al.  Machine-learned multi-system surrogate models for materials prediction , 2018, npj Computational Materials.

[17]  S. Froyen,et al.  Brillouin-zone integration by Fourier quadrature: Special points for superlattice and supercell calculations. , 1989, Physical review. B, Condensed matter.

[18]  Adarsh Balasubramanian,et al.  Fast, accurate, and transferable many-body interatomic potentials by symbolic regression , 2019, npj Computational Materials.

[19]  Lattice-constant and band-gap tuning in wurtzite and zincblende BInGaN alloys , 2019, Journal of Applied Physics.

[20]  Tim Mueller,et al.  Improved Prediction of Nanoalloy Structures by the Explicit Inclusion of Adsorbates in Cluster Expansions , 2018, The Journal of Physical Chemistry C.

[21]  Gus L. W. Hart,et al.  A robust algorithm for $k$-point grid generation and symmetry reduction. , 2018 .

[22]  M. Tinkham.  Group Theory and Quantum Mechanics , 1964 .

[23]  K. Livi,et al.  Mechanistic Insights for Low-Overpotential Electroreduction of CO2 to CO on Copper Nanowires , 2017 .

[24]  M. Dupuis,et al.  Synthesis and Anisotropic Electrocatalytic Activity of Covellite Nanoplatelets with Fixed Thickness and Tunable Diameter. , 2018, ACS applied materials & interfaces.

[25]  J. Neugebauer,et al.  The Basics of Electronic Structure Theory for Periodic Systems , 2019, Front. Chem..

[26]  Thomas J. Kempa,et al.  Substrate-directed synthesis of MoS2 nanocrystals with tunable dimensionality and optical properties , 2019, Nature Nanotechnology.

[27]  Tim Mueller,et al.  Ensemble Effect in Bimetallic Electrocatalysts for CO2 Reduction. , 2019, Journal of the American Chemical Society.

[28]  L. Williams,et al.  BInGaN alloys nearly lattice-matched to GaN for high-power high-efficiency visible LEDs , 2017, 1902.02692.

[29]  Yanli Wang,et al.  Tunable Electronic Structures of Hydrogenated Zigzag and Armchair Dumbbell Silicene Nanosheets: A Computational Study , 2018, The Journal of Physical Chemistry C.