I have tried to run a geometry optimisation on a 2D water slab using VDW-DF exchange functional, dipole correction, etc., etc.
I find that by starting from different initial states, I can end up with two equivalent supercells, one of which has a rectangular basis, and the other has a parallelogram as base (angle about 60 degrees). The total energies after minimisation agree to 7 decimal places. So far so good. The in-plane stresses are very different, however: the rectangular one is 2.7 GPa, and the parallelogram one is 2.2 GPa (both in plane stresses are equal as one would expect). Is there any reason why that should be the case? Anyone aware of any pitfalls with the calculation of stresses in the vdW part of the code? I am using QE 5.1.2. Niels --- Prof. Niels R. Walet Phone: +44(0)1613063693 School of Physics and Astronomy Fax: +44(0)1613064303 The University of Manchester Mobile: +44(0)7516622121 Manchester, M13 9PL, UK room 7.7, Schuster Building email: [email protected] web: http://www.theory.physics.manchester.ac.uk/~mccsnrw<../../owa/redir.aspx?C=8f9098c1943846dcba17536b8618a38a&URL=http%3a%2f%2fwww.theory.physics.manchester.ac.uk%2f%7emccsnrw>
_______________________________________________ Pw_forum mailing list [email protected] http://pwscf.org/mailman/listinfo/pw_forum
