Residual Entropy of Ice: A Study Based on Transfer Matrices

Abstract

Residual entropy of ice systems has long been a significant and intriguing issue in condensed matter physics and statistical mechanics. The exact solutions for the residual entropy of realistic three-dimensional ice systems remain unknown. In this study, we focus on two typical realistic ice systems, namely the hexagonal ice (ice Ih) and cubic ice (ice Ic). We present a transfer matrix description of the number of ice-ruled configurations for these two systems. First, a transfer matrix $M$ is constructed for ice Ic, where each element is the number of ice-ruled configurations of a hexagonal monolayer under certain condition. The product of $M$ and its transpose corresponds to a bilayer unit in ice Ih lattice, therefore is exactly a transfer matrix for ice Ih. Making use of this, we simply show that the residual entropy of ice Ih is not less than that of ice Ic in the thermodynamic limit, which was first proved by Onsager in 1960s. Furthermore, we find an alternative transfer matrix $M^{\prime }$ for ice Ih, which is based on a monolayer periodic unit. Some interesting properties of $M$, $M{M}^T$ and $M^{\prime }$ are illustrated, specifically the summation of all elements, the element in the first row and first column, and the trace. Each property is equivalent with the residual entropy of a two-dimensional ice model. Our work rediscovers the relationship between the residual entropies of ice Ih and ice Ic, and provides an effective description for various two-dimensional ice models.

30 pages, 7 figures

References in corpus (25)

• An efficient, multiple range random walk algorithm to calculate the density of states
• Matrix Product States, Projected Entangled Pair States, and variational renormalization group methods for quantum spin systems
• Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram
• Ising models of quantum frustration
• Fast Algorithm to Calculate Density of States
• Why spin ice obeys the ice rules
• A new macroscopically degenerate ground state in the spin ice compound Dy2Ti2O7 under a magnetic field
• Magnetization process of spin ice in a [111] magnetic field
• Dimers on two-dimensional lattices
• Residual entropies for three-dimensional frustrated spin systems with tensor networks
• Residual Entropy of Ordinary Ice from Multicanonical Simulations
• Corrections to Pauling residual entropy and single tetrahedron based approximations for the pyrochlore lattice Ising antiferromagnet
• Configurational entropy of hydrogen-disordered ice polymorphs
• Configurational entropy of ice from thermodynamic integration
• Thermodynamics of the classical spin-ice model with nearest neighbour interactions using the Wang-Landau algorithm
• Entropy of the diluted antiferromagnetic Ising models on the frustrated lattices using the Wang-Landau method
• Large peaks in the entropy of the diluted nearest-neighbor spin-ice model on the pyrochlore lattice in a [111] magnetic field
• Monte Carlo study of degenerate groundstates and residual entropy in a frustrated honeycomb lattice Ising model
• A Transfer Matrix study of the staggered BCSOS model
• The odd eight-vertex model
• Residual entropy of the dilute Ising chain in a magnetic field
• Polarization Plateaus in Hexagonal Water Ice $I_h$
• Calculation of the Residual Entropy of Ice Ih by Monte Carlo simulation with the Combination of the Replica-Exchange Wang-Landau algorithm and Multicanonical Replica-Exchange Method
• Residual Entropy of a Two-dimensional Ising Model with Crossing and Four-spin Interactions
• Exact Results for the Residual Entropy of Ice Hexagonal Monolayer

Cited by in corpus (1)

• Equivalence of residual entropy of hexagonal and cubic ices from tensor network methods