<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Manoj. the forces you are getting seem really small to me.<div>If you enforce cubic symmetry, forces vanish exactly. If you do not, the code itself has to "figure out" that the atoms are indeed at their proper equilibrium positions, and this can only be achieved within the given overall numerical accuracy. One thing that may explicitly (although very slightly) break cubic symmetry in your supercell is tha FFT grid. Not that it makes any physical sense, but if you really want to overkill zero fores, you may want to play a little bit with the kinetic-energy cutoff and/or FFT grid. SB</div><div><br><div><div>On Oct 7, 2009, at 11:23 PM, Manoj Srivastava wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear PWSCF users, <br> I have a question regarding relaxation of low symmetry lattices. I did<br>relaxation of 5 atoms FCC unit cell of copper. This is pure copper, no<br>interface, no defect or anything. The unit cell is just not conventional,<br>but a different way of choosing bulk FCC. For the default value of<br>convergence the calculation finishes at the end of 1st scf cycle, but I<br>see different forces when I describe nosym=.TRUE. or cell_parameters<br>{cubic}. Following is the input file.<br><br>&control<br> calculation='relax'<br> pseudo_dir = '/home/manoj/espresso-4.0.4/pseudo',<br> outdir='./',<br> prefix='lcu',<br> wf_collect=.TRUE.<br> /<br> &system<br> ibrav =0,<br> celldm(1)=6.81650937063832<br> nat= 5,<br> ntyp= 1,<br> ecutwfc =50.0,<br> occupations='smearing',<br> smearing='gaussian',<br> degauss=0.02,<br> ecutrho=500<br> /<br>&electrons<br> conv_thr = 1.0e-8<br> mixing_beta=0.7<br>/<br>&ions<br>ion_dynamics='bfgs'<br>pot_extrapolation='second_order'<br>wfc_extrapolation='second_order'<br>/<br>ATOMIC_SPECIES<br> Cu 63.55 Cu.pz-d-rrkjus.UPF<br>ATOMIC_POSITIONS<br> Cu 0.0 0.0 0.0<br> Cu 0.5 0.5 0.0<br> Cu 0.0 1.0 0.0<br> Cu 1.0 1.0 0.0<br> Cu 0.5 1.5 0.0<br><br>K_POINTS (automatic)<br>8 8 8 0 0 0<br>CELL_PARAMETERS {cubic}<br> 1.5 0.5 0.0<br>-0.5 1.5 0.0<br> 0.0 0.5 0.5<br><br>Result of relaxation : The covergence is achieved with follwing forces on<br>the atoms. <br>Forces acting on atoms (Ry/au):<br><br> Forces acting on atoms (Ry/au):<br><br> atom 1 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 2 type 1 force = 0.00000567 0.00000185 -0.00000172<br> atom 3 type 1 force = 0.00001155 0.00000377 -0.00000372<br> atom 4 type 1 force = -0.00001155 -0.00000377 0.00000372<br> atom 5 type 1 force = -0.00000567 -0.00000185 0.00000172<br><br><br>I realized that I should not have chosen cubic symmetry group, as my unit<br>cell is not cubic. So, as described in the input file parameters, for the<br>low symmetry lattices, I specifed nosym=.TRUE., and again the relaxation<br>stopped after 1 scf cycle, acheiving convergence. But the forces on the<br>atoms are - <br><br>Forces acting on atoms (Ry/au):<br><br> atom 1 type 1 force = 0.00001355 0.00000524 -0.00000516<br> atom 2 type 1 force = -0.00001493 -0.00000395 0.00000395<br> atom 3 type 1 force = 0.00000253 -0.00000043 0.00000045<br> atom 4 type 1 force = 0.00000820 0.00000230 -0.00000239<br> atom 5 type 1 force = -0.00000935 -0.00000316 0.00000315<br><br>So, still forces are .001 ev/a.u, which i think is higher for a perfect<br>FCC lattice. I can specify stricter convergence of force, but that would<br>move atoms around and break symmetry. So, as long as the lattice parameter<br>is correct, I think for perfect FCC lattice, forces on atoms should be<br>smaller than what they are presently. I have checked the lattice constant<br>for bulk FCC with a different calculation, so lattice constant is not an<br>issue. Another possibility is symmetry group, as you can see with the same<br>input variables, changing symmetry group changes results.I was wondering<br>if there is a way i can choose correct symmetry group for this low<br>symmetry lattice which might give me smaller forces? Are there any other<br>parameters I need to specify? On a broad view I am curious to know how to<br>do relaxation of low symmetry lattices.<br><br>Regards, <br>Manoj Srivastava<br>University of Florida<br>Gainesville, USA<br><br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></div></blockquote></div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">---</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Stefano Baroni - SISSA</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">&</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">DEMOCRITOS National Simulation Center - Trieste</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://stefano.baroni.me">http://stefano.baroni.me</a> [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><br style="font-size: 12px; "></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><span class="Apple-style-span" style="font-size: 14px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">La morale est une logique de l'action comme la logique est une morale de la pensée - Jean Piaget</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br></span></font></div></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Please, if possible, don't</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">send me MS Word or PowerPoint attachments</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Why? See:</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://www.gnu.org/philosophy/no-word-attachments.html">http://www.gnu.org/philosophy/no-word-attachments.html</a></span></font></div><div><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br class="webkit-block-placeholder"></span></font></div></div></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"> </div><br></div></body></html>