<div dir="ltr">Hello, I'm currently performing a CELL_OPT minimisation of a slab to model an external surface where there is approximately 14 Å of vacuum between each slab in the z-axis. As I want to maintain the slab model, I'm constraining the z vector whilst allowing the a and b lattice vectors for my slab to relax along with my nuclei. However, I keep obtaining negative pressure values above 100 Bar and my system is unable to converge as the cell keeps wanting to collapse.<div><br></div><div>I've tried,<div>- Changing the minimiser from BFGS to CG</div><div>- Constraining and fixing the cell angels</div><div>- Increasing the EPS_SCF to obtain more accurate forces</div><div>- Varying the multi-grid cutoff energies</div><div>But sadly I've had no success with obtain reasonable pressure deviations. None of the forces are particularly large on the atoms (the largest are on the order of E-5 a.u.) so I'm unsure why my pressure is so negative. Identical slabs can successfully converge with slightly smaller and larger vacuum sizes so I'm unsure why certain vacuum separations are unsuccessful. I'm running CP2K version 6.1 and I've attached my input and the final step.</div><div><br></div><div>Does anyone have any suggestions for me to try?</div></div></div>