Hydrogen-bonded ferroelectrics are modelled by a coupled spin/nonlinear lattice (polarizability) interaction Hamiltonian, where specifically the geometry of the hydrogen bond is included. The model leads to a structural phase transition and describes correctly the isotope effect due to the substitution H/D in hydrogen-bonded systems in terms of bond length changes.
