Hello, everyone<br><br> I have a question regarding how to set up lattice in pwscf. Usually I like present lattice in the following format ( aka vasp format)<br><br>lattice:<br> a0 (in angstrom)<br><br> a11 a12 a13<br>
a21 a22 a23<br> a31 a32 a33<br><br>and atomic coordinates in the fractional (crystal) format. So, I guessed that in pwscf it should be something like that:<br><br>ibrav = 0 <br>celldm(1) = 1.88973 (corresponds 1 Ang)
<br><br>CELL_PARAMETERS cubic<br> a11 a12 a13<br>
a21 a22 a23<br>
a31 a32 a33<br>
ATOMIC_POSITIONS {crystal}<br>.....<br>Using this representation for my system (slab) I didn't get convergence at all, and after a long careful examination, I found that this is due to wrong geometry representation. To confirm this fact, I took the simple system from example25 in espresso distribution and changed lattice like:
<br><br>from <br> ibrav= 0, celldm(1)=8.19<br> CELL_PARAMETERS<br> 0.50 0.50 1.00 <br> 0.50 1.00 0.50<br> 1.00 0.50 0.50<br>ATOMIC_POSITIONS {crystal}<br><br>to<br> ibrav= 0, celldm(1)=1.00 ,<br> CELL_PARAMETERS
<br> 4.95 4.95 8.19<br> 4.95 8.19 4.95 <br> 8.19 4.95 4.95<br><br><br>I thought these are equivalent representation, but the code give me different cell volumes and I obtained the different total energy, completely different magnetic moments and long convergence in the second case. So, I suggest that this representation (case 2) is not correct. Am I right?
<br><br> Best regards,<br> JK<br>