Simple refinements of Brillouin zone integration

Calculations of thermodynamic properties of crystals by means of quasi-harm onic lattice dynamics require numerical integrations over the Brillouin zon e, using successively finer grids to achieve convergence ro the required pr ecision; but for complex crystals convergence may be uneconomically slow. A model for orthorhombic polyethylene is used to show how convergence may be improved (1) at low temperatures by taking successively finer grids close to the origin of reciprocal space, and (2) at all temperatures by using a t hree-dimensional Simpson's rule.