The isotope-effect in the phase transition of KDP: New insights from ab initio path-integral simulations

TL;DR

This study uses ab initio path-integral simulations to explore how quantum effects influence isotope-dependent phase transitions in KDP, revealing distinct mechanisms for protons and deuterons.

Contribution

It provides new insights into isotope effects in KDP by distinguishing quantum delocalization from vibration-assisted hopping using advanced simulations.

Findings

Different phase transition mechanisms for protons and deuterons

Qualitative differences in paraelectric phase behavior between isotopes

Agreement with experimental neutron scattering data

Abstract

We investigate the quantum-mechanical localization of protonated and deterated isotopes in the symmetric low-barrier hydrogen-bonds of potassium dihydrogen phosphate (KDP) crystals in the paraelectric phase. The spatial density distributions of these hydrogen atoms are suspected to be responsible for the surprisingly large isotope effect observed for the ferroelectric phase transition in KDP. We employ ab initio path integral molecular dynamics simulations to obtain the nuclear real-space and momentum-space densities n(R) and n(k) of protons and deuterons, which are compared to experimental Neutron Compton Scattering data. Our results suggest a qualitative difference in the nature of the paraelectic phase in KDP between the two isotopes. We are able to discriminate between real quantum delocalization and vibration-assisted hopping and thus provide evidence for two distinct mechanisms of the ferroelectric phase transition in this class of materials.