Unfortunately, the calculation of forces is not yet implemented <br>with the fully rotational DFT+U+J scheme (lda_plus_u_kind=1).<br>This works of course for a simplified DFT+U scheme of Matteo<br>(lda_plus_u_kind=0). Probably it is not that difficult but <br>
I did not look yet ...<br>regards,<br>Alexander <br><br><div class="gmail_quote">2012/11/10 nuttachai jutong <span dir="ltr"><<a href="mailto:nuttachai.jutong@gmail.com" target="_blank">nuttachai.jutong@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Matteo<br>
Here is my in put file for relaxation. Actually It work if I use only<br>
'lda_plus_u = .TRUE.' , I mean only LDA+U without J.<br>
<br>
 &control<br>
    calculation='vc-relax'<br>
    restart_mode='from_scratch',<br>
    pseudo_dir = '/home/ltmp/jutongnu/espresso-5.0/PW/EuO-LDAU_Espresso/Pseudo',<br>
    outdir='/home/ltmp/jutongnu/espresso-5.0/PW/EuO-LDAU_Espresso/Vary_Uf_Fix_all/EuO_Uf_7.1_Jf_0.77_cut70_k8'<br>
    prefix='EuO_Uf_7.1_Jf_0.77_cut70_k8',<br>
    etot_conv_thr= 1.0D-4,<br>
    forc_conv_thr= 1.0D-3,<br>
    wf_collect=.true.<br>
    verbosity = 'high'<br>
 /<br>
 &system<br>
    ibrav=  0, a= 1,<br>
    nat=  2, ntyp= 2,<br>
    ecutwfc = 70.0, ecutrho = 630,<br>
    nspin=2, starting_magnetization(1)=1,<br>
    occupations='smearing', smearing='fd', degauss=0.003,<br>
    nosym  =.FALSE.<br>
    lda_plus_u = .TRUE.,  lda_plus_u_kind= 1,<br>
    Hubbard_U(1)=7.1 , Hubbard_J(1,1)= 0, Hubbard_J(2,1)= 0,<br>
Hubbard_J(3,1)= 0.77,<br>
    Hubbard_U(2)=4.6 , Hubbard_J(1,2)=1.2,<br>
 /<br>
 &electrons<br>
    conv_thr = 1.0d-10<br>
    mixing_beta = 0.3<br>
    electron_maxstep=1000,<br>
 /<br>
&IONS<br>
 ion_dynamics="bfgs"<br>
/<br>
&CELL<br>
cell_dynamics = 'damp-w',<br>
press = 0.0 ,<br>
/<br>
<br>
ATOMIC_SPECIES<br>
 Eu  151.964  063-Eu-ca-nspd-hanghui.uspp.format.UPF<br>
 O   16       O.pz-rrkjus.UPF<br>
CELL_PARAMETERS (angstroms)<br>
 0.000    2.57195  2.57195<br>
 2.57195  0.000    2.57195<br>
 2.57195  2.57195  0.000<br>
ATOMIC_POSITIONS (angstroms)<br>
Eu   0.00000000    0.00000000    0.00000000   1 1 1<br>
O    2.57195000    2.57195000    2.57195000   1 1 1<br>
K_POINTS automatic<br>
8 8 8 0 0 0<br>
<br>
Best regard,<br>
 Nuttachai<br>
<br>
><br>
> Message: 8<br>
> Date: Sat, 10 Nov 2012 04:37:47 -0600<br>
> From: Matteo Cococcioni <<a href="mailto:matteo@umn.edu">matteo@umn.edu</a>><br>
> Subject: Re: [Pw_forum] Relaxation with LDA+U+J<br>
> To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br>
> Message-ID:<br>
>       <<a href="mailto:CAMZASgFhmaNmTDTq99ck4X82A_Z1iZUoRLn-ZdQxQB7QfMvyWw@mail.gmail.com">CAMZASgFhmaNmTDTq99ck4X82A_Z1iZUoRLn-ZdQxQB7QfMvyWw@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset=ISO-8859-1<br>
><br>
> Dear Nuttachai,<br>
><br>
> could you please post the entire input you used for this calculation?<br>
><br>
> from the error message you are getting it seems that you are asking<br>
> for something not yet implemented. However it is difficult to guess<br>
> what that is without seeing your input.<br>
><br>
> Matteo<br>
><br>
><br>
><br>
> On Fri, Nov 9, 2012 at 8:24 PM, nuttachai jutong<br>
> <<a href="mailto:nuttachai.jutong@gmail.com">nuttachai.jutong@gmail.com</a>> wrote:<br>
>> Dear all,<br>
>>  I have tried to calculate a gemoetry relaxation with LDA+U+J, however<br>
>> I got the error<br>
>> ' Error in routine force_hub (1):<br>
>>       forces in full LDA+U scheme are not yet implemented'<br>
>><br>
>> I performed calculation the electronics band structure with LDA+U+J, by QE<br>
>> 5.0.<br>
>> 'lda_plus_u = .TRUE.,  lda_plus_u_kind= 1,'<br>
>><br>
>> Everything go very well, however I got error when I try to calculate<br>
>> the gemoetry relaxation by<br>
>><br>
>> 'calculation  = "vc-relax'.<br>
>><br>
>> and I set<br>
>><br>
>> &IONS<br>
>>  ion_dynamics="bfgs"<br>
>> /<br>
>> &CELL<br>
>> cell_dynamics = 'damp-w',<br>
>> press = 0.0 ,<br>
>> /<br>
>><br>
>>  Dose someone please suggest me How do I solve this problem ?<br>
>> Is it a bug ? or I have to set some flag ?<br>
>><br>
>> Best regard,<br>
>>  Nuttachai<br>
>> _______________________________________________<br>
>> Pw_forum mailing list<br>
>> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
>> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
><br>
><br>
> --<br>
> Matteo Cococcioni<br>
> Department of Chemical Engineering and Materials Science,<br>
> University of Minnesota<br>
> 421 Washington Av. SE<br>
> Minneapolis, MN 55455<br>
> Tel. <a href="tel:%2B1%20612%20624%209056" value="+16126249056">+1 612 624 9056</a>    Fax <a href="tel:%2B1%20612%20626%207246" value="+16126267246">+1 612 626 7246</a><br>
><br>
><br>
> ------------------------------<br>
><br>
> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
> End of Pw_forum Digest, Vol 65, Issue 23<br>
> ****************************************<br>
><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br>