<div dir="ltr"><div>Hi All,</div><div><br></div><div>I am trying to calculate the Equation of State curve for Na2S unit cell using CP2K. <br></div><div><br></div><div>To start off the calculations, I am performing a cell optimization of a single unit cell of Na2S. Please find attached the input for cell optimization. After optimization I obtain cell parameters to be a=b=c=4.575 as compared to (cell parameters from VASP of same unit cell a=b=c=4.646). The lattice parameter values are fairly close to one another. However, the fractional coordinates change significantly in the final optimized unit cell from the starting configuration. <br></div><div><br></div><div>How do I resolve this issue? I am currently using GPW method with a very high cut off of 1300 Ry to account for hard basis sets for Na. I have also tried using GAPW method with a low cutoff of 550 Ry but with GAPW method I obtain the charge density on grid value to be ~ 0.0000005694 as compared to charge density on grid of 0.0000000003 for GPW methods. Also find attached input file for GAPW method. With GAPW method I get the lattice parameters to be a=b=c=4.642. I would like to know if I should proceed with GAPW method or GPW methods to calculate the equation of state curves for Na2S?<br></div><div><br></div><div>Regards,</div><div>Aniruddha M Dive</div><div>PhD Candidate</div><div>Washington State University<br></div></div>