Why is there no low-temperature phase transition in NaOH?

When Quantum Fluctuations Meet Structural Instabilities: The Isotope- and Pressure-Induced Phase Transition in the Quantum Paraelectric NaOH

Anhydrous sodium hydroxide, a common and structurally simple compound, shows spectacular isotope effects: NaOD undergoes a first-order transition, which is absent in NaOH. By combining ab initio electronic structure calculations with Feynman path integrals, we show that NaOH is an unusual example of a quantum paraelectric: zero-point quantum fluctuations stretch the weak hydrogen bonds (HBs) into a region where they are unstable and break. By strengthening the HBs via isotope substitution or applied pressure, the system can be driven to a broken-symmetry antiferroelectric phase. In passing, we provide a simple quantitative criterion for HB breaking in layered crystals and show that nuclear quantum effects are crucial in paraelectric to ferroelectric transitions in hydrogen-bonded hydroxides.

Thermal and Nuclear Quantum Effects at the Antiferroelectric to Paraelectric Phase Transition in KOH and KOD Crystals

Crystalline KOH undergoes an antiferroelectric (AFE) proton ordering phase transition at low temperatures, which results in a monoclinic bilayer structure held together by a network of weak hydrogen bonds (HBs). The Curie temperature shifts up when the compound is deuterated, an effect that classical MD is not able to catch. For deeper insights into the transition mechanism, we carry out ab initio MD simulations of KOH and KOD crystals by including quantum effects on the nuclei through Feynman path integrals. The geometric isotope effect and the evolution of the lattice parameters with temperature agree with the experimental data, while the purely classical description is not appropriate. Our results show that deuteration strengthens the HBs in the low-T AFE ordered phase. The transition is characterized by the flipping of OH/OD groups along a bending mode. Above the transition, the system is driven into a dynamical disordered paraelectric phase.

Anharmonic spectral features via trajectory-based quantum dynamics: A perturbative analysis of the interplay between dynamics and sampling

The performance of different approximate algorithms for computing anharmonic features in vibrational spectra is analyzed and compared on model and more realistic systems that present relevant nuclear quantum effects. The methods considered combine approximate sampling of the quantum thermal distribution with classical time propagation and include Matsubara dynamics, path integral dynamics approaches, linearized initial value representation, and the recently introduced adaptive quantum thermal bath. A perturbative analysis of these different methods enables us to account for the observed numerical performance on prototypes for overtones and combination bands and to draw qualitatively correct trends for the numerical results obtained for Fermi resonances. Our results prove that the unequal performances of these approaches often derive from the method employed to sample initial conditions and not, as usually assumed, from the lack of coherence in the time propagation. Furthermore, as confirmed by the analysis reported in Benson and Althorpe, J. Chem. Phys. 130, 194510 (2021), we demonstrate, both via the perturbative approach and numerically, that path integral dynamics methods fail to reproduce the intensities of these anharmonic features and follow purely classical trends with respect to their temperature behavior. Finally, the remarkably accurate performance of the adaptive quantum thermal bath approach is documented and motivated.

Pressure-induced localisation of the hydrogen-bond network in KOH-VI

Using a combination of ab initio crystal structure prediction and neutron diffraction techniques, we have solved the full structure of KOH-VI at 7 GPa. Rather than being orthorhombic and proton-ordered as had previously be proposed, we find that this high-pressure phase of potassium hydroxide is tetragonal (space group I4/mmm) and proton disordered. It has an unusual hydrogen bond topology, where the hydroxyl groups form isolated hydrogen-bonded square planar (OH)4 units. This structure is stable above 6.5 GPa and, despite being macroscopically proton-disordered, local ice rules enforce microscopic order of the hydrogen bonds. We suggest the use of this novel type of structure to study concerted proton tunneling in the solid state, while the topology of the hydrogen bond network could conceivably be exploited in data storage applications based solely on the manipulations of hydrogen bonds. The unusual localisation of the hydrogen bond network under applied pressure is found to be favored by a more compact packing of the constituents in a distorted cesium chloride structure.