<div>Dear all,</div><div>I am trying to optimize a single layer of Al (111) in the (XY free) plane by constraining the Z component. Al (111) is stretched a = ~ 3 Ang.</div><div>Using a pure PW + PBE scheme, the system converges towards a lattice as reported in the attached file Al_111.pdb</div><div>Using cp2k I cannot reproduce the same result because the system reorganizes itself into smaller metal domains, breaking the original Al (111) lattice.</div><div>I tried with several XCs, but I could not reproduce the "metallic state" obtained as in the PW scheme, which is necessary to treat a more complex system.</div><div>Therefore, I don't really care about the correct behavior of this simple system, I want try to find a way to reproduce PW results using cp2k.</div><div>Do you have any suggestions on how to bias the result of cp2k towards the result obtained with PW-PBE scheme?</div><div><br></div><div>Thank you,</div><div><br></div><div>Luca</div><div><br></div><div><br></div>