skip to main content

Title: Residual entropy of ices and clathrates from Monte Carlo simulation

We calculated the residual entropy of ices (Ih, Ic, III, V, VI) and clathrates (I, II, H), assuming the same energy of all configurations satisfying the Bernal–Fowler ice rules. The Metropolis Monte Carlo simulations in the range of temperatures from infinity to a size-dependent threshold were followed by the thermodynamic integration. Convergence of the simulation and the finite-size effects were analyzed using the quasichemical approximation and the Debye–Hückel theory applied to the Bjerrum defects. The leading finite-size error terms, ln N/N, 1/N, and for the two-dimensional square ice model also 1/N{sup 3/2}, were used for an extrapolation to the thermodynamic limit. Finally, we discuss the influence of unequal energies of proton configurations.
Authors:
 [1]
  1. Department of Physical Chemistry, Institute of Chemical Technology, 166 28 Prague 6 (Czech Republic)
Publication Date:
OSTI Identifier:
22304337
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Physics; Journal Volume: 140; Journal Issue: 20; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; CLATHRATES; COMPUTERIZED SIMULATION; CONFIGURATION; ENTROPY; MONTE CARLO METHOD; PROTONS; TWO-DIMENSIONAL CALCULATIONS