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

Abstract

The long-standing question of whether the residual entropy of hexagonal ice ($S_h$) equals that of cubic ice ($S_c$) remains unresolved despite decades of research on ice-type models. While analytical studies have established the inequality $S_h\ge S_c$, numerical investigations suggest that the two values are very close. In this work, we revisit this problem using high-precision tensor-network methods. In Monte Carlo approaches the residual entropy cannot be directly obtained by sampling the ground-state degeneracy space, however, the tensor-network framework enables an explicit encoding of the "ice rule'' into local tensors, and then the residual entropy is transformed into finding the largest eigenvalue of a transfer operator in the form of a projected entangled-pair operator, which allows high-accuracy numerical evaluation. Meanwhile, we propose a new perspective based on analyzing the normality of the transfer operator, and demonstrate that if the operator is normal, the equality $S_h=S_c$ follows directly. Then the variational tensor network methods are employed to numerically verify this normality. Finally both residual entropies are directly computed by using our recently developed split corner transfer matrix renormalization group algorithm, providing a rigorous evidence supporting the equality between $S_h$ and $S_c$.

10 pages, 9 figures, three tables

References in corpus (19)

• Matrix Product States and Projected Entangled Pair States: Concepts, Symmetries, and Theorems
• Ice structures, patterns, and processes: A view across the ice-fields
• Advances in the experimental exploration of water's phase diagram
• Magnetization process of spin ice in a [111] magnetic field
• Fermionic projected entangled-pair states and topological phases
• Configurational entropy of hydrogen-disordered ice polymorphs
• A scaling hypothesis for projected entangled-pair states
• Solving frustrated Ising models using tensor networks
• Partial lifting of degeneracy in the $J_1-J_2-J_3$ Ising antiferromagnet on the kagome lattice
• Artificial out-of-plane Ising antiferromagnet on the kagome lattice with very small further neighbour couplings
• Large peaks in the entropy of the diluted nearest-neighbor spin-ice model on the pyrochlore lattice in a [111] magnetic field
• Unified tensor network theory for frustrated classical spin models in two dimensions
• Matrix product state fixed points of non-Hermitian transfer matrices
• 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
• Variationally optimizing infinite projected entangled-pair states at large bond dimensions: A split corner transfer matrix renormalization group approach
• Exact Results for the Residual Entropy of Ice Hexagonal Monolayer
• Quantum Anharmonic Effects in Hydrogen-Bond Symmetrization of High-Pressure Ice
• Efficient optimization of variational tensor-network approach to three-dimensional statistical systems
• Residual Entropy of Ice: A Study Based on Transfer Matrices