<div dir="ltr">Hi, <div>I am trying to calculate the binding energies of multiple adsorbates on a tin oxide surface (a 4 layer slab), as well as the zero point energies. When I relax the geometry of the structure (only the top two layers and adsorbates are allowed to relax) using a gamma point calculation and run a hessian analysis on the optimized structure, all calculated frequencies are real. But, when I run a geometry relaxation with k-Points on the optimized structure from a gamma point calculation and run a subsequent hessian calculation also using k-Points, multiple imaginary frequencies are observed. I have verified that the geometries obtained from both gamma point and k-Point calculation are almost the same. Can anyone help me with what is going wrong here ? </div><div><br></div><div>All input and output files are attached. </div><div><br></div><div>ads_gamma_rlxn.inp/.out (gamma point geometry relaxation) </div><div>ads_gamma_hessian.inp/.out (gamma point hessian analysis)</div><div><br></div><div><div>ads_kP_rlxn.inp/.out (k-point geometry relaxation) </div><div>ads_kP_hessian.inp/.out (k-point hessian analysis)</div></div><div><br></div><div>Thanks for your help,</div><div><br></div><div>- Karthik </div><div><br></div><div><br></div><div><br></div></div>