Geometry optimization of zirconium sulfophenylphosphonate layers by molecular simulation methods