Hi Stefano,<br><br>Thanks for the reply! <br><br>I tried your patch and it did not fix the eigenvalues. So I changed it a little bit:<br><br>Original:<br>       <i>IF ( first .AND. istep == 0 .AND. &<br>
                    starting_pot == 'atomic' ) CALL ns_adj() ;  rhoin%ns
    = rho%ns</i><br>
    <br>New one:<br>       <i>IF ( first .AND. istep == 0 .AND. & <br>       starting_pot == 'atomic' ) CALL ns_adj();<br>          <br>       IF ( first .AND. istep == 0 .AND. &<br>       starting_pot == 'atomic' ) rhoin%ns = rho%ns;</i><br>
<br>I tested both cases on the same setup (using the intel compiler), and the second one works out. Maybe it should be updated on cvs?<br><br>Regards,<br><br>Giovani Faccin<br>IFGW/Unicamp - Brazil<br><br><br><br><div class="gmail_quote">
2011/3/29 Stefano de Gironcoli <span dir="ltr"><<a href="mailto:degironc@sissa.it">degironc@sissa.it</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">


  
    
    
  
  <div text="#000000" bgcolor="#ffffff">
    I think this is a bug ...  <br>
    <br>
    The initial value of ns should be overwritten by the input
    request...<br>
    A simple way to obtain this should be to replace, in electrons.f90,
    the line<br>
    <br>
                IF ( first .AND. istep == 0 .AND. &<br>
                    starting_pot == 'atomic' ) CALL ns_adj() <br>
    by<br>
    <br>
              IF ( first .AND. istep == 0 .AND. &<br>
                    starting_pot == 'atomic' ) CALL ns_adj() ;  rhoin%ns
    = rho%ns<br>
    <br>
    I updated the cvs version.<br><font color="#888888">
    <br>
    stefano</font><div><div></div><div class="h5"><br>
    <br>
    <br>
    On 03/29/2011 06:55 PM, Giovani Faccin wrote:
    <blockquote type="cite">
      <pre>Dear Gabriele,

Yes, I tried mixing_fixed_ns. The problem is that this command fixes the
initial occupation during all the run, not the occupation imposed on step 2
with starting_ns eigenvalues. I'm attaching an example input and output file
that shows this happening.

As you can see on the output.txt file, after iteraction 1 the program
inserts the eigenvalues I wanted on atom 1 spin2 and atom 2 spin2. However,
just after this it prints:


RESET ns to initial values (iter <= mixing_fixed_ns)

on every new iteration. Looking at the end of the calculation, you can see
that the final states are exactly the default states from step 0. In other
words, it overwrites the eigenvalues defined by starting_ns eigenvalues,
placing the default initial occupations instead. Perhaps an "if" clause
could be added to the code to manage situations when both mixing_fixed_ns
and starting_ns eigenvalues are defined by the user? Or maybe I'm doing this
all wrong?

Regards,

Giovani Faccin

IFGW/Unicamp-Brazil





2011/3/29 Gabriele Sclauzero <a href="mailto:sclauzer@sissa.it" target="_blank"><sclauzer@sissa.it></a>

</pre>
      <blockquote type="cite">
        <pre>Dear Giovani,

    have you tried with mixing_fixed_ns?

    Could you please provide your affiliation?


Regards,


GS

mixing_fixed_ns
Il giorno 29/mar/2011, alle ore 18.16, Giovani Faccin ha scritto:

Hi all!

Is it possible, in a LDA+U calculation, to fix the starting_ns eigenvalues
for the entire calculation run, instead of just the second interaction?

Specifically, I'd like to be able to track how a certain electronic
configuration's associated energy varies when other parameters (ex: atomic
distances) change. But in order to do that, I need to be able to impose
specific electronic configurations for the system, even if those
configurations are not the global energy mininum.

Can that be done in QE?

Thanks!

--
Giovani


_______________________________________________
Pw_forum mailing list
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>



§ Gabriele Sclauzero, EPFL SB ITP CSEA
*   PH H2 462, Station 3, CH-1015 Lausanne*


_______________________________________________
Pw_forum mailing list
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>


</pre>
      </blockquote>
      <pre></pre>
      <pre><fieldset></fieldset>
_______________________________________________
Pw_forum mailing list
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>
</pre>
    </blockquote>
    <br>
  </div></div></div>

<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>
<br></blockquote></div><br><br clear="all"><br>-- <br>Giovani<br>