Yes, the enthalpy is decreasing the enthalpy is decreasing. Below the values of the  enthaply.<br><br>Best regards,<br>André Luis Martinotto<br>
<br>Andre Martinotto<br>Email: almartinotto at <a href="http://gmail.com/" target="_blank">gmail.com</a><br>Computing Department<br>Universidade de Caxias do Sul<br>Caxias do Sul - RS, Brazil<br><br>     enthalpy old            =   -1018.4377035275 Ry<br>
     enthalpy new            =   -1018.4494900797 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4494900797 Ry<br>     enthalpy new            =   -1018.4568547798 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>
     enthalpy old            =   -1018.4568547798 Ry<br>     enthalpy new            =   -1018.4624673172 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4624673172 Ry<br>     enthalpy new            =   -1018.4675637420 Ry<br>
     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4675637420 Ry<br>     enthalpy new            =   -1018.4745374801 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4745374801 Ry<br>
     enthalpy new            =   -1018.4787658712 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4787658712 Ry<br>     enthalpy new            =   -1018.4839453395 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>
     enthalpy old            =   -1018.4839453395 Ry<br>     enthalpy new            =   -1018.4906317150 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4906317150 Ry<br>     enthalpy new            =   -1018.4927767110 Ry<br>
     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4927767110 Ry<br>     enthalpy new            =   -1018.4939707446 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4939707446 Ry<br>
     enthalpy new            =   -1018.4946012005 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4946012005 Ry<br>     enthalpy new            =   -1018.4951369286 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>
     enthalpy old            =   -1018.4951369286 Ry<br>     enthalpy new            =   -1018.4955690712 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4955690712 Ry<br>     enthalpy new            =   -1018.4961210855 Ry<br>
     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4961210855 Ry<br>     enthalpy new            =   -1018.4964627786 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4964627786 Ry<br>
     enthalpy new            =   -1018.4971749436 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4971749436 Ry<br>     enthalpy new            =   -1018.4985789257 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>
     enthalpy old            =   -1018.4985789257 Ry<br>     enthalpy new            =   -1018.4996031008 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.4996031008 Ry<br>     enthalpy new            =   -1018.5008017865 Ry<br>
     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5008017865 Ry<br>     enthalpy new            =   -1018.5036232518 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5036232518 Ry<br>
     enthalpy new            =   -1018.5070595112 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5070595112 Ry<br>     enthalpy new            =   -1018.5117630729 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>
     enthalpy old            =   -1018.5117630729 Ry<br>     enthalpy new            =   -1018.5170591719 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5170591719 Ry<br>     enthalpy new            =   -1018.5222913399 Ry<br>
     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5222913399 Ry<br>     enthalpy new            =   -1018.5282241124 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5282241124 Ry<br>
     enthalpy new            =   -1018.5342597827 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5342597827 Ry<br>     enthalpy new            =   -1018.5406966898 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>
     enthalpy old            =   -1018.5406966898 Ry<br>     enthalpy new            =   -1018.5469217329 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5469217329 Ry<br>     enthalpy new            =   -1018.5530869862 Ry<br>
     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5530869862 Ry<br>     enthalpy new            =   -1018.5599277199 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5599277199 Ry<br>
     enthalpy new            =   -1018.5668277117 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>     enthalpy old            =   -1018.5668277117 Ry<br>     enthalpy new            =   -1018.5740583719 Ry<br>     CASE: enthalpy_new &lt; enthalpy_old<br>
     enthalpy old            =   -1018.5740583719 Ry<br>     enthalpy new            =   -1018.5816796313 Ry<br><br><span style="visibility: visible;" id="main"><span style="visibility: visible;" id="search"><em></em></span></span><br>
<div class="gmail_quote">On Fri, Aug 20, 2010 at 2:53 PM, Stefano de Gironcoli <span dir="ltr">&lt;<a href="mailto:degironc@sissa.it">degironc@sissa.it</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
is the enthalpy decreasing during the relaxation ?<br>
<div><div></div><div class="h5">stefano<br>
<br>
Andre Martinotto wrote:<br>
&gt; Dear Stefano,<br>
&gt;<br>
&gt; Sorry.  I sent the wrong input file. The outputs correspond to the file<br>
&gt; below.<br>
&gt;<br>
&gt; Best regards,<br>
&gt; André Luis Martinotto<br>
&gt;<br>
&gt; Andre Martinotto<br>
&gt; Email: almartinotto at <a href="http://gmail.com" target="_blank">gmail.com</a><br>
&gt; Computing Department<br>
&gt; Universidade de Caxias do Sul<br>
&gt; Caxias do Sul - RS, Brazil<br>
&gt;<br>
&gt; &amp;CONTROL<br>
&gt;                        title = CoSb3 ,<br>
&gt;                  calculation = &#39;vc-relax&#39;,<br>
&gt;                 restart_mode = &#39;restart&#39;,<br>
&gt;                   wf_collect = .true. ,<br>
&gt;                       outdir = &#39;/dados/almartin/est2_a&#39; ,<br>
&gt;                       wfcdir = &#39;/dados/almartin/est2_a&#39; ,<br>
&gt;                   pseudo_dir = &#39;/home/u/almartin/cosb3/pseudo&#39; ,<br>
&gt;                       prefix = &#39;est1&#39; ,<br>
&gt;                etot_conv_thr = 1.0D-4 ,<br>
&gt;                forc_conv_thr = 1.0D-3 ,<br>
&gt;                      tstress = .true. ,<br>
&gt;                      tprnfor = .true.<br>
&gt;  /<br>
&gt;  &amp;SYSTEM<br>
&gt;                        ibrav = 1,<br>
&gt;                            A = 8.5012 ,<br>
&gt;                          nat = 32,<br>
&gt;                         ntyp = 2,<br>
&gt;                      ecutwfc = 30. ,<br>
&gt;                      ecutrho = 300 ,<br>
&gt;                      input_dft = &#39;PBE&#39; ,<br>
&gt;                  occupations = &#39;smearing&#39; ,<br>
&gt;                      degauss = 0.01 ,<br>
&gt;                     smearing = &#39;marzari-vanderbilt&#39; ,<br>
&gt;  /<br>
&gt;  &amp;ELECTRONS<br>
&gt;                     conv_thr = 1.0D-9 ,<br>
&gt;                  startingpot = &#39;atomic&#39; ,<br>
&gt;                  startingwfc = &#39;atomic&#39; ,<br>
&gt;                  mixing_mode = &#39;plain&#39; ,<br>
&gt;                  mixing_beta = 0.7D0 ,<br>
&gt;              diagonalization = &#39;david&#39; ,<br>
&gt;  /<br>
&gt;  &amp;IONS<br>
&gt;                 ion_dynamics = &#39;bfgs&#39; ,<br>
&gt;                ion_positions = &#39;default&#39; ,<br>
&gt;  /<br>
&gt;  &amp;CELL<br>
&gt;                cell_dynamics = &#39;bfgs&#39; ,<br>
&gt;                        press = 0.D0 ,<br>
&gt;               press_conv_thr = 0.2D0 ,<br>
&gt;                  cell_factor = 1.5D0 ,<br>
&gt;  /<br>
&gt; ATOMIC_SPECIES<br>
&gt;    Co   58.93300  Co.pbe-nd-rrkjus.UPF<br>
&gt;    Sb  121.76000  Sb.pbe-rrkjus.UPF<br>
&gt; ATOMIC_POSITIONS crystal<br>
&gt; Co      0.25    0.25    0.25    1 1 1<br>
&gt; Co      0.75    0.75    0.25    1 1 1<br>
&gt; Co      0.75    0.25    0.75    1 1 1<br>
&gt; Co      0.25    0.75    0.75    1 1 1<br>
&gt; Co      0.75    0.75    0.75    1 1 1<br>
&gt; Co      0.25    0.25    0.75    1 1 1<br>
&gt; Co      0.25    0.75    0.25    1 1 1<br>
&gt; Co      0.75    0.25    0.25    1 1 1<br>
&gt; Sb      0       0.6646  0.8421  1 1 1<br>
&gt; Sb      0       0.3354  0.8421  1 1 1<br>
&gt; Sb      0       0.6646  0.1579  1 1 1<br>
&gt; Sb      0       0.3354  0.1579  1 1 1<br>
&gt; Sb      0.8421  0       0.6646  1 1 1<br>
&gt; Sb      0.6646  0.8421  0       1 1 1<br>
&gt; Sb      0.8421  0       0.3354  1 1 1<br>
&gt; Sb      0.3354  0.8421  0       1 1 1<br>
&gt; Sb      0.1579  0       0.6646  1 1 1<br>
&gt; Sb      0.6646  0.1579  0       1 1 1<br>
&gt; Sb      0.1579  0       0.3354  1 1 1<br>
&gt; Sb      0.  0.  0.       0 0 0<br>
&gt; Sb      0.5     0.1646  0.3421  1 1 1<br>
&gt; Sb      0.5     0.8354  0.3421  1 1 1<br>
&gt; Sb      0.5     0.1646  0.6579  1 1 1<br>
&gt; Sb      0.5     0.8354  0.6579  1 1 1<br>
&gt; Sb      0.3421  0.5     0.1646  1 1 1<br>
&gt; Sb      0.1646  0.3421  0.5     1 1 1<br>
&gt; Sb      0.3421  0.5     0.8354  1 1 1<br>
&gt; Sb      0.8354  0.3421  0.5     1 1 1<br>
&gt; Sb      0.6579  0.5     0.1646  1 1 1<br>
&gt; Sb      0.1646  0.6579  0.5     1 1 1<br>
&gt; Sb      0.6579  0.5     0.8354  1 1 1<br>
&gt; Sb      0.5  0.5  0.5     0 0 0<br>
&gt; K_POINTS automatic<br>
&gt;   8 8 8 0 0 0<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; On Fri, Aug 20, 2010 at 12:17 PM, Stefano de Gironcoli &lt;<a href="mailto:degironc@sissa.it">degironc@sissa.it</a>&gt;wrote:<br>
&gt;<br>
&gt;<br>
&gt;&gt; you ARE NOT making a vc-relaxation unless you specify<br>
&gt;&gt; calculation=&#39;vc-relax&#39;<br>
&gt;&gt; stefano<br>
&gt;&gt;<br>
&gt;&gt; Andre Martinotto wrote:<br>
&gt;&gt;<br>
&gt;&gt;&gt; Dear all,<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; I&#39;m trying to make a vc-relax optimization of a structure with vacancies,<br>
&gt;&gt;&gt; but I&#39;m having some problems for convergence. I have a problem with<br>
&gt;&gt;&gt;<br>
&gt;&gt; forces<br>
&gt;&gt;<br>
&gt;&gt;&gt; convergence. Below the input file as are the force and stress results.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Someone would have some suggestion to solve this problem?<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Thanks!<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Best regards,<br>
&gt;&gt;&gt; André Luis Martinotto<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Andre Martinotto<br>
&gt;&gt;&gt; Email: almartinotto at <a href="http://gmail.com" target="_blank">gmail.com</a><br>
&gt;&gt;&gt; Computing Department<br>
&gt;&gt;&gt; Universidade de Caxias do Sul<br>
&gt;&gt;&gt; Caxias do Sul - RS, Brazil<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; /------------------------------------- input file<br>
&gt;&gt;&gt; ---------------------------------------/<br>
&gt;&gt;&gt;  &amp;CONTROL<br>
&gt;&gt;&gt;                        title = CoSb3 ,<br>
&gt;&gt;&gt;                  calculation = &#39;relax&#39;,<br>
&gt;&gt;&gt;                 restart_mode = &#39;restart&#39;,<br>
&gt;&gt;&gt;                   wf_collect = .true. ,<br>
&gt;&gt;&gt;                       outdir = &#39;/dados/almartin/est2_a&#39; ,<br>
&gt;&gt;&gt;                       wfcdir = &#39;/dados/almartin/est2_a&#39; ,<br>
&gt;&gt;&gt;                   pseudo_dir = &#39;/home/u/almartin/cosb3/pseudo&#39; ,<br>
&gt;&gt;&gt;                       prefix = &#39;est1&#39; ,<br>
&gt;&gt;&gt;                etot_conv_thr = 1.0D-4 ,<br>
&gt;&gt;&gt;                forc_conv_thr = 1.0D-3 ,<br>
&gt;&gt;&gt;                      tstress = .true. ,<br>
&gt;&gt;&gt;                      tprnfor = .true.<br>
&gt;&gt;&gt;  /<br>
&gt;&gt;&gt;  &amp;SYSTEM<br>
&gt;&gt;&gt;                        ibrav = 1,<br>
&gt;&gt;&gt;                            A = 8.5012 ,<br>
&gt;&gt;&gt;                          nat = 32,<br>
&gt;&gt;&gt;                         ntyp = 2,<br>
&gt;&gt;&gt;                      ecutwfc = 30. ,<br>
&gt;&gt;&gt;                      ecutrho = 300 ,<br>
&gt;&gt;&gt;                      input_dft = &#39;PBE&#39; ,<br>
&gt;&gt;&gt;                  occupations = &#39;smearing&#39; ,<br>
&gt;&gt;&gt;                      degauss = 0.01 ,<br>
&gt;&gt;&gt;                     smearing = &#39;marzari-vanderbilt&#39; ,<br>
&gt;&gt;&gt;  /<br>
&gt;&gt;&gt;  &amp;ELECTRONS<br>
&gt;&gt;&gt;                     conv_thr = 1.0D-9 ,<br>
&gt;&gt;&gt;                  startingpot = &#39;atomic&#39; ,<br>
&gt;&gt;&gt;                  startingwfc = &#39;atomic&#39; ,<br>
&gt;&gt;&gt;                  mixing_mode = &#39;plain&#39; ,<br>
&gt;&gt;&gt;                  mixing_beta = 0.7D0 ,<br>
&gt;&gt;&gt;              diagonalization = &#39;david&#39; ,<br>
&gt;&gt;&gt;  /<br>
&gt;&gt;&gt;  &amp;IONS<br>
&gt;&gt;&gt;                 ion_dynamics = &#39;bfgs&#39; ,<br>
&gt;&gt;&gt;                ion_positions = &#39;default&#39; ,<br>
&gt;&gt;&gt;  /<br>
&gt;&gt;&gt; &amp;CELL<br>
&gt;&gt;&gt;                cell_dynamics = &#39;bfgs&#39; ,<br>
&gt;&gt;&gt;                        press = 0.D0 ,<br>
&gt;&gt;&gt;               press_conv_thr = 0.2D0 ,<br>
&gt;&gt;&gt;                  cell_factor = 1.5D0 ,<br>
&gt;&gt;&gt;  /<br>
&gt;&gt;&gt; ATOMIC_SPECIES<br>
&gt;&gt;&gt;    Co   58.93300  Co.pbe-nd-rrkjus.UPF<br>
&gt;&gt;&gt;    Sb  121.76000  Sb.pbe-rrkjus.UPF<br>
&gt;&gt;&gt; ATOMIC_POSITIONS crystal<br>
&gt;&gt;&gt; Co      0.25    0.25    0.25    1 1 1<br>
&gt;&gt;&gt; Co      0.75    0.75    0.25    1 1 1<br>
&gt;&gt;&gt; Co      0.75    0.25    0.75    1 1 1<br>
&gt;&gt;&gt; Co      0.25    0.75    0.75    1 1 1<br>
&gt;&gt;&gt; Co      0.75    0.75    0.75    1 1 1<br>
&gt;&gt;&gt; Co      0.25    0.25    0.75    1 1 1<br>
&gt;&gt;&gt; Co      0.25    0.75    0.25    1 1 1<br>
&gt;&gt;&gt; Co      0.75    0.25    0.25    1 1 1<br>
&gt;&gt;&gt; Sb      0       0.6646  0.8421  1 1 1<br>
&gt;&gt;&gt; Sb      0       0.3354  0.8421  1 1 1<br>
&gt;&gt;&gt; Sb      0       0.6646  0.1579  1 1 1<br>
&gt;&gt;&gt; Sb      0       0.3354  0.1579  1 1 1<br>
&gt;&gt;&gt; Sb      0.8421  0       0.6646  1 1 1<br>
&gt;&gt;&gt; Sb      0.6646  0.8421  0       1 1 1<br>
&gt;&gt;&gt; Sb      0.8421  0       0.3354  1 1 1<br>
&gt;&gt;&gt; Sb      0.3354  0.8421  0       1 1 1<br>
&gt;&gt;&gt; Sb      0.1579  0       0.6646  1 1 1<br>
&gt;&gt;&gt; Sb      0.6646  0.1579  0       1 1 1<br>
&gt;&gt;&gt; Sb      0.1579  0       0.3354  1 1 1<br>
&gt;&gt;&gt; Sb      0.  0.  0.           0 0 0<br>
&gt;&gt;&gt; Sb      0.5     0.1646  0.3421  1 1 1<br>
&gt;&gt;&gt; Sb      0.5     0.8354  0.3421  1 1 1<br>
&gt;&gt;&gt; Sb      0.5     0.1646  0.6579  1 1 1<br>
&gt;&gt;&gt; Sb      0.5     0.8354  0.6579  1 1 1<br>
&gt;&gt;&gt; Sb      0.3421  0.5     0.1646  1 1 1<br>
&gt;&gt;&gt; Sb      0.1646  0.3421  0.5     1 1 1<br>
&gt;&gt;&gt; Sb      0.3421  0.5     0.8354  1 1 1<br>
&gt;&gt;&gt; Sb      0.8354  0.3421  0.5     1 1 1<br>
&gt;&gt;&gt; Sb      0.6579  0.5     0.1646  1 1 1<br>
&gt;&gt;&gt; Sb      0.1646  0.6579  0.5     1 1 1<br>
&gt;&gt;&gt; Sb      0.6579  0.5     0.8354  1 1 1<br>
&gt;&gt;&gt; Sb      0.5  0.5  0.5         0 0 0<br>
&gt;&gt;&gt; K_POINTS automatic<br>
&gt;&gt;&gt;   8 8 8 0 0 0<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; /------------------------------------- stress results<br>
&gt;&gt;&gt; --------------------------------------/<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;     entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.77<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 9.48<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 11.91<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 4.78<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -6.89<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -2.65<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 1.31<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 5.44<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 2.47<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -2.82<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -1.30<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.08<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.71<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 1.11<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.70<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.16<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -0.77<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -1.30<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -1.45<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -1.23<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -0.56<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.39<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 1.24<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 1.92<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 2.18<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 2.16<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 1.82<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 1.35<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.75<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; 0.30<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -0.03<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -0.37<br>
&gt;&gt;&gt;      entering subroutine stress ...<br>
&gt;&gt;&gt;           total   stress  (Ry/bohr**3)                   (kbar)     P=<br>
&gt;&gt;&gt; -0.72<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; /------------------------------------- force results<br>
&gt;&gt;&gt; --------------------------------------/<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;      Total force =     0.072779     Total SCF correction =     0.000066<br>
&gt;&gt;&gt;      Total force =     0.053372     Total SCF correction =     0.000053<br>
&gt;&gt;&gt;      Total force =     0.044120     Total SCF correction =     0.000023<br>
&gt;&gt;&gt;      Total force =     0.045980     Total SCF correction =     0.000024<br>
&gt;&gt;&gt;      Total force =     0.049875     Total SCF correction =     0.000027<br>
&gt;&gt;&gt;      Total force =     0.043061     Total SCF correction =     0.000055<br>
&gt;&gt;&gt;      Total force =     0.033627     Total SCF correction =     0.000027<br>
&gt;&gt;&gt;      Total force =     0.024939     Total SCF correction =     0.000032<br>
&gt;&gt;&gt;      Total force =     0.022018     Total SCF correction =     0.000009<br>
&gt;&gt;&gt;      Total force =     0.015573     Total SCF correction =     0.000049<br>
&gt;&gt;&gt;      Total force =     0.013346     Total SCF correction =     0.000039<br>
&gt;&gt;&gt;      Total force =     0.010870     Total SCF correction =     0.000017<br>
&gt;&gt;&gt;      Total force =     0.010344     Total SCF correction =     0.000031<br>
&gt;&gt;&gt;      Total force =     0.016943     Total SCF correction =     0.000021<br>
&gt;&gt;&gt;      Total force =     0.015429     Total SCF correction =     0.000031<br>
&gt;&gt;&gt;      Total force =     0.013450     Total SCF correction =     0.000013<br>
&gt;&gt;&gt;      Total force =     0.015799     Total SCF correction =     0.000010<br>
&gt;&gt;&gt;      Total force =     0.023243     Total SCF correction =     0.000032<br>
&gt;&gt;&gt;      Total force =     0.028594     Total SCF correction =     0.000015<br>
&gt;&gt;&gt;      Total force =     0.038065     Total SCF correction =     0.000030<br>
&gt;&gt;&gt;      Total force =     0.044262     Total SCF correction =     0.000073<br>
&gt;&gt;&gt;      Total force =     0.048789     Total SCF correction =     0.000054<br>
&gt;&gt;&gt;      Total force =     0.053071     Total SCF correction =     0.000047<br>
&gt;&gt;&gt;      Total force =     0.057755     Total SCF correction =     0.000030<br>
&gt;&gt;&gt;      Total force =     0.061205     Total SCF correction =     0.000057<br>
&gt;&gt;&gt;      Total force =     0.063732     Total SCF correction =     0.000055<br>
&gt;&gt;&gt;      Total force =     0.065333     Total SCF correction =     0.000069<br>
&gt;&gt;&gt;      Total force =     0.067697     Total SCF correction =     0.000029<br>
&gt;&gt;&gt;      Total force =     0.071297     Total SCF correction =     0.000042<br>
&gt;&gt;&gt;      Total force =     0.075401     Total SCF correction =     0.000020<br>
&gt;&gt;&gt;      Total force =     0.080403     Total SCF correction =     0.000064<br>
&gt;&gt;&gt;      Total force =     0.086012     Total SCF correction =     0.000030<br>
&gt;&gt;&gt;      Total force =     0.091511     Total SCF correction =     0.000033<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; ------------------------------------------------------------------------<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; _______________________________________________<br>
&gt;&gt;&gt; Pw_forum mailing list<br>
&gt;&gt;&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt;&gt;&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; Pw_forum mailing list<br>
&gt;&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt;&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt; ------------------------------------------------------------------------<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Pw_forum mailing list<br>
&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br>