<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000">Thank you for your email,<br><br>I have performed convergence studies ( both on 1x1 and 2x2 unit cell) for k-point, wave-function cutoff, smearing. <br><br><b>Here is the input for Ti bulk</b>,<br><br>&control<br>   calculation      = 'scf'<br>   restart_mode     = 'from_scratch'<br>   pseudo_dir       = '/work/shared/s-tih/pseudo/'<br>   prefix           = 'Ti2_deg4_40_8k'<br>   wf_collect       = .true.<br>   tstress          = .true.<br>   tprnfor          = .true.<br>   forc_conv_thr    = 1.0d-5<br>   verbosity        = 'high'<br> /<br>&system<br>   ibrav            = 4<br>   a = 2.950,b=2.950,c=4.81735,cosbc=0,cosac=0,cosab=-0.5<br>   nat              = 2<br>   ntyp             = 1<br>   ecutwfc          = 40<br>   ecutrho          = 320<br>   occupations      = 'smearing',  smearing='mp', degauss=0.04D0<br> /<br>&electrons<br>   diagonalization  = 'cg'<br>   mixing_beta      = 0.3d00<br>   conv_thr         = 1.0d-7<br> <br>/<br>ATOMIC_SPECIES<br>  Ti  47.8670  Ti.pbe-spn-rrkjus_psl.1.0.0.UPF<br><br><br>ATOMIC_POSITIONS (crystal)<br>   Ti     0.666 0.333 0.7500<br>   Ti     0.333 0.666 0.2500<br><br>K_POINTS {automatic}<br>  8  8  5  0  0  0<br><br><b>And the output</b>, </div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000">..................................................................................................................<br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000">    bravais-lattice index     =            4<br>     lattice parameter (alat)  =       5.5747  a.u.<br>     unit-cell volume          =     245.0076 (a.u.)^3<br>     number of atoms/cell      =            2<br>     number of atomic types    =            1<br>     number of electrons       =        24.00<br>     number of Kohn-Sham states=           16<br>     kinetic-energy cutoff     =      40.0000  Ry<br>     charge density cutoff     =     320.0000  Ry<br>     convergence threshold     =      1.0E-07<br>     mixing beta               =       0.3000<br>     number of iterations used =            8  plain     mixing<br>     Exchange-correlation      = PBE ( 1  4  3  4 0 0)<br><br>     celldm(1)=   5.574692  celldm(2)=   1.000000  celldm(3)=   1.633000<br>     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000<br><br>     crystal axes: (cart. coord. in units of alat)<br>               a(1) = (   1.000000   0.000000   0.000000 )<br>               a(2) = (  -0.500000   0.866025   0.000000 )<br>               a(3) = (   0.000000   0.000000   1.633000 )<br><br>     reciprocal axes: (cart. coord. in units 2 pi/alat)<br>               b(1) = (  1.000000  0.577350  0.000000 )<br>               b(2) = (  0.000000  1.154701  0.000000 )<br>               b(3) = (  0.000000  0.000000  0.612370 )<br><br><br>     PseudoPot. # 1 for Ti read from file:<br>     /work/shared/s-tih/pseudo/Ti.pbe-spn-rrkjus_psl.1.0.0.UPF<br>     MD5 check sum: e281089c08e14b8efcf92e44a67ada65<br>     Pseudo is Ultrasoft + core correction, Zval = 12.0<br>     Generated using "atomic" code by A. Dal Corso v.6.2.2<br>     Using radial grid of 1177 points,  6 beta functions with:<br>                l(1) =   0<br>                l(2) =   0<br>                l(3) =   1<br>                l(4) =   1<br>                l(5) =   2<br>                l(6) =   2<br>     Q(r) pseudized with 0 coefficients<br><br><br>     atomic species   valence    mass     pseudopotential<br>        Ti            12.00    47.86700     Ti( 1.00)<br><br>      8 Sym. Ops., with inversion, found ( 6 have fractional translation)<br>..................................................................................................................<br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000">     the Fermi energy is    12.6981 ev<br><br>!    total energy              =    -239.45976063 Ry<br>     Harris-Foulkes estimate   =    -239.45976064 Ry<br>     estimated scf accuracy    <          6.2E-09 Ry<br><br>     The total energy is the sum of the following terms:<br><br>     one-electron contribution =     -84.68427794 Ry<br>     hartree contribution      =      49.69752058 Ry<br>     xc contribution           =     -36.99691545 Ry<br>     ewald contribution        =    -167.48081236 Ry<br>     smearing contrib. (-TS)   =       0.00472453 Ry<br><br>     convergence has been achieved in   7 iterations<br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000"><b>cohesive Energy  = -239.45976063/2 - -119.34098597 = -0.388894345 Ry = -5.29 eV.   </b>(0.44 eV difference from experimental value).<br><br>I am attaching the i/o files also with this, but most of the important details of the calculation are above.<br><br>Thanks for any help.</div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000"><br></div><div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:12.8px"><b style="font-family:"Times New Roman";font-size:medium">------------------------------------------------------------------------------------------------------</b></div><div style="line-height:22px"><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Ajmalghan MUTHALI</span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px"><br></span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Post doctorate researcher</span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Laboratoire ICB </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">UMR 6303 CNRS-Université de Bourgogne </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">9 Avenue Alain Savary, BP 47870 </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">F-21078 DIJON Cedex, France </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Tel: +33-(0)7.69.28.19.91 </span></font></div></div><div style="line-height:22px"><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px">Email : </span><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px"><a href="mailto:ajmalghan.muthali@u-bourgogne.fr" target="_blank">ajmalghan.muthali@u-bourgogne.fr</a></span></font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Nov 12, 2019 at 12:25 PM Giuseppe Mattioli <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
Dear Ajmalghan<br>
Sorry for asking a possibly stupid question, but you are focusing on  <br>
the calculation of isolated Ti, and the error might be contained in  <br>
the calculation of hcp metal Ti... Are you sure that everything is  <br>
correct in that case?<br>
HTH<br>
Giuseppe<br>
<br>
Quoting Ajmal Ghan <<a href="mailto:ajmalghan1991@gmail.com" target="_blank">ajmalghan1991@gmail.com</a>>:<br>
<br>
> Thanks all for the reply,<br>
><br>
> With all the inputs provided here and mail archives, I made some<br>
> significant changes ( fixed magnetization, increased the size of the cell,<br>
> Gamma point calculation, Mixing beta etc...).<br>
><br>
> *Input of isolated Ti atom*:<br>
><br>
> &control<br>
>    calculation      = 'scf'<br>
>    restart_mode     = 'from_scratch'<br>
>    pseudo_dir       = '/work/shared/s-tih/pseudo/'<br>
>    prefix           = 'Tifree_deg1_40_gk_1'<br>
>    wf_collect       = .true.<br>
>    forc_conv_thr    = 1.0d-5<br>
>    verbosity        = 'high'<br>
>  /<br>
> &system<br>
>    ibrav            = 1<br>
>    celldm(1)        = 30<br>
>    nat              = 1<br>
>    ntyp             = 1<br>
>    ecutwfc          = 40<br>
>    ecutrho          = 320<br>
>    nspin            = 2<br>
>    tot_magnetization = 2<br>
>    nosym            = .true<br>
>    nbnd             = 100<br>
>    occupations      = 'fixed'<br>
>  /<br>
> &electrons<br>
>    diagonalization  = 'cg'<br>
>    mixing_beta      = 0.3d00<br>
>    conv_thr         = 1.0d-7<br>
> /<br>
> ATOMIC_SPECIES<br>
>   Ti  47.8670  Ti.pbe-spn-rrkjus_psl.1.0.0.UPF<br>
><br>
><br>
> ATOMIC_POSITIONS (crystal)<br>
>    Ti     0.5 0.5 0.5<br>
><br>
> K_POINTS GAMMA<br>
><br>
><br>
> And for bulk calculations, i used same forc_conv_thr, ecutwfc, &electrons<br>
> parameters. *But I am still getting 5.23 eV as cohesive energy for Ti*. is<br>
> it possible to get the experimental ( 4.85 eV) using DFT calculations as<br>
> reported in some of the journals?<br>
><br>
> The output of isolated Ti atom calculation looks like,<br>
><br>
>     bravais-lattice index     =            1<br>
>      lattice parameter (alat)  =      30.0000  a.u.<br>
>      unit-cell volume          =   27000.0000 (a.u.)^3<br>
>      number of atoms/cell      =            1<br>
>      number of atomic types    =            1<br>
>      number of electrons       =        12.00 (up:   7.00, down:   5.00)<br>
>      number of Kohn-Sham states=          100<br>
>      kinetic-energy cutoff     =      40.0000  Ry<br>
>      charge density cutoff     =     320.0000  Ry<br>
>      convergence threshold     =      1.0E-07<br>
>      mixing beta               =       0.3000<br>
>      number of iterations used =            8  plain     mixing<br>
>      Exchange-correlation      = PBE ( 1  4  3  4 0 0)<br>
><br>
>      celldm(1)=  30.000000  celldm(2)=   0.000000  celldm(3)=   0.000000<br>
>      celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000<br>
><br>
> etc..............................................................................................................<br>
><br>
>      Starting magnetic structure<br>
>      atomic species   magnetization<br>
>         Ti           0.000<br>
><br>
>      No symmetry found<br>
> etc..............................................................................................................<br>
> ..................................................................................................................<br>
><br>
>      iteration # 23     ecut=    40.00 Ry     beta= 0.30<br>
>      CG style diagonalization<br>
>      ethr =  1.01E-09,  avg # of iterations =  3.1<br>
><br>
>      negative rho (up, down):  1.660E-02 1.306E-01<br>
><br>
>      Magnetic moment per site:<br>
>      atom:    1    charge:   11.9999    magn:    2.0000    constr:    0.0000<br>
><br>
>      total cpu time spent up to now is      108.5 secs<br>
><br>
>      End of self-consistent calculation<br>
><br>
>  ------ SPIN UP ------------<br>
><br>
>           k = 0.0000 0.0000 0.0000 ( 57657 PWs)   bands (ev):<br>
><br>
>    -62.1874 -38.4854 -38.4348 -38.4332  -4.4517  -4.4503  -4.4048  -3.5611<br>
>     -3.5610  -3.5153  -1.4572  -1.4565  -1.1485  -0.4725   0.0828   0.0992<br>
>      0.1709   0.1890   0.1914   0.2168   0.5306   0.5437   0.6058   0.6689<br>
>      0.6794   0.6917   0.7097   0.7127   0.8758   0.8811   0.9122   0.9174<br>
>      1.1513   1.1544   1.1806   1.2368   1.4025   1.4198   1.4417   1.4878<br>
>      1.5485   1.5709   1.9392   1.9444   1.9629   2.0040   2.0535   2.0683<br>
>      2.1492   2.1983   2.2084   2.3100   2.3235   2.3454   2.3845   2.4064<br>
>      2.4189   2.4639   2.4693   2.4865   2.4880   2.5054   2.5055   2.5179<br>
>      2.5188   2.5215   2.7061   2.7211   2.7420   2.7869   2.7963   2.8090<br>
>      2.8145   2.8304   2.9938   3.0003   3.0113   3.0175   3.0230   3.0259<br>
>      3.0280   3.0332   3.0529   3.0580   3.0645   3.0690   3.0757   3.0785<br>
>      3.0808   3.0898   3.5616   3.5684   3.6347   3.8452   3.9717   3.9774<br>
>      4.0381   4.0472   4.0540   4.0617<br>
><br>
>      occupation numbers<br>
>      1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000<br>
><br>
>  ------ SPIN DOWN ----------<br>
>           k = 0.0000 0.0000 0.0000 ( 57657 PWs)   bands (ev):<br>
><br>
>    -60.5692 -37.0219 -36.7521 -36.7481  -3.9282  -2.9464  -2.8480  -2.8461<br>
>     -2.0737  -2.0734  -1.3881  -0.9733  -0.9724  -0.4455   0.0983   0.1054<br>
>      0.2139   0.2301   0.2409   0.2476   0.5776   0.5939   0.6083   0.6823<br>
>      0.6922   0.7172   0.7202   0.7223   0.9007   0.9617   0.9637   0.9653<br>
>      1.1886   1.2019   1.2109   1.2572   1.4214   1.4579   1.4604   1.5175<br>
>      1.5802   1.6023   1.9618   1.9874   1.9882   2.0402   2.1348   2.1844<br>
>      2.2059   2.2275   2.2450   2.3219   2.3301   2.3507   2.4059   2.4124<br>
>      2.4170   2.4823   2.4877   2.4881   2.5037   2.5077   2.5266   2.5297<br>
>      2.5305   2.5376   2.7444   2.7669   2.7706   2.8302   2.8379   2.8443<br>
>      2.8558   2.8670   3.0199   3.0266   3.0335   3.0380   3.0397   3.0427<br>
>      3.0462   3.0493   3.0781   3.0795   3.0805   3.0813   3.0847   3.0869<br>
>      3.1009   3.1038   3.6043   3.6940   3.6955   3.9171   4.0363   4.0648<br>
>      4.0726   4.0791   4.0896   4.0980<br>
><br>
>      occupation numbers<br>
>      1.0000   1.0000   1.0000   1.0000   1.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000<br>
>      0.0000   0.0000   0.0000   0.0000<br>
><br>
>      highest occupied, lowest unoccupied level (ev):    -3.9282   -3.5611<br>
><br>
> !    total energy              =    -119.34098597 Ry<br>
>      Harris-Foulkes estimate   =    -119.34098367 Ry<br>
>      estimated scf accuracy    <       0.00000010 Ry<br>
><br>
>      The total energy is the sum of the following terms:<br>
><br>
>      one-electron contribution =    -167.22501663 Ry<br>
>      hartree contribution      =      79.78699354 Ry<br>
>      xc contribution           =     -18.28393498 Ry<br>
>      ewald contribution        =     -13.61902790 Ry<br>
><br>
>      total magnetization       =     2.00 Bohr mag/cell<br>
>      absolute magnetization    =     2.44 Bohr mag/cell<br>
><br>
>      convergence has been achieved in  23 iterations<br>
><br>
><br>
> Waiting for reply. Thanks for all the inputs.<br>
><br>
> Thanks & Regards,<br>
> *------------------------------------------------------------------------------------------------------*<br>
> Ajmalghan MUTHALI<br>
><br>
> Post doctorate researcher<br>
> Laboratoire ICB<br>
> UMR 6303 CNRS-Université de Bourgogne<br>
> 9 Avenue Alain Savary, BP 47870<br>
> F-21078 DIJON Cedex, France<br>
> Tel: +33-(0)7.69.28.19.91<br>
> Email : <a href="mailto:ajmalghan.muthali@u-bourgogne.fr" target="_blank">ajmalghan.muthali@u-bourgogne.fr</a><br>
><br>
><br>
> On Sun, Nov 10, 2019 at 3:52 PM Ari P Seitsonen <<a href="mailto:Ari.P.Seitsonen@iki.fi" target="_blank">Ari.P.Seitsonen@iki.fi</a>><br>
> wrote:<br>
><br>
>><br>
>> Dear Ajmalghan,<br>
>><br>
>>    Some quick comments:<br>
>><br>
>>   - Why do you use k point in the case of an isolated atom?? Well, that<br>
>> should not matter, just that you are wasting computing time<br>
>><br>
>>   - I guess that the spherical symmetry of the atom is broken; thus I<br>
>> would<br>
>> use a non-cube cell, preferably the orthorhombic cell, with slightly<br>
>> different lengths of the basis vectors of the unit cell. Then to break the<br>
>> symmetry, you can use some randomisation of the initial wave functions.<br>
>> And still, the convergence is probably going to be very difficult... You<br>
>> can indeed try to fix the magnetisation; and I would reduce the<br>
>> 'mixing_beta' to something (very) small - in principle already at the<br>
>> first step the electron density should be close to the self-consistent<br>
>> one, bar the loss of sphericality and the randomised wave function<br>
>><br>
>>   - Please remember that the scale of 'starting_magnetisation' is from -1<br>
>> to +1, meaning that all the electrons are spin-polarised either up or<br>
>> down, whereas in your case you only want to polarise the two valence<br>
>> electrons out of the valence of the pseudo potential that seems to be 12<br>
>> electrons<br>
>><br>
>>    Well, Good Luck. :)<br>
>><br>
>>      Greetings from Paris,<br>
>><br>
>>         apsi<br>
>><br>
>><br>
>> -=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-<br>
>>    Ari Paavo Seitsonen / <a href="mailto:Ari.P.Seitsonen@iki.fi" target="_blank">Ari.P.Seitsonen@iki.fi</a> / <a href="http://www.iki.fi/~apsi/" rel="noreferrer" target="_blank">http://www.iki.fi/~apsi/</a><br>
>>      Ecole Normale Supérieure (ENS), Département de Chimie, Paris<br>
>>      Mobile (F) : +33 789 37 24 25    (CH) : +41 79 71 90 935<br>
>><br>
>><br>
>> On Sat, 9 Nov 2019, Ajmal Ghan wrote:<br>
>><br>
>> > Hello all,<br>
>> ><br>
>> > I have gone through all the archived discussion about cohesive energy<br>
>> calculation.<br>
>> ><br>
>> > I am trying to calculate the cohesive energy of Ti which is  4.85eV. But<br>
>> I am getting 5.23 eV which is closer. But I have found from a previous<br>
>> discussion<br>
>> > here (<br>
>> <a href="https://www.mail-archive.com/users@lists.quantum-espresso.org/msg11410.html" rel="noreferrer" target="_blank">https://www.mail-archive.com/users@lists.quantum-espresso.org/msg11410.html</a><br>
>> ) that, the final magnetization should be 2 Bohr mag/ cell at the end<br>
>> > of calculation for Ti and smearing contribution of energy should be<br>
>> 0.0Ry.<br>
>> > But I am getting to get 3.83 Bohr mag/cell and a significant smearing<br>
>> contribution = -0.03295688 Ry.<br>
>> ><br>
>> > Anyone help me to sort this. what should I change in the input?<br>
>> ><br>
>> ><br>
>> > here is my input,<br>
>> ><br>
>> > &control<br>
>> >    calculation      = 'scf'<br>
>> >    restart_mode     = 'from_scratch'<br>
>> >    pseudo_dir       = '/work/shared/s-tih/pseudo/'<br>
>> >    prefix           = 'Tifree_deg1_50_8k_1'<br>
>> >    wf_collect       = .true.<br>
>> >    tstress          = .true.<br>
>> >    tprnfor          = .true.<br>
>> >    forc_conv_thr    = 1.0d-6<br>
>> >    verbosity        = 'high'<br>
>> >  /<br>
>> > &system<br>
>> >    ibrav            = 1<br>
>> >    celldm(1)     = 20<br>
>> >    nat              = 1<br>
>> >    ntyp            = 1<br>
>> >    ecutwfc       = 50      //( I have performed convergence study. But<br>
>> since degauss is reduced to 0.01, i increased ecut)<br>
>> >    ecutrho       = 400<br>
>> >    nspin          = 2<br>
>> >    starting_magnetization(1) = 1        // I think, the final<br>
>> magnetisation should be 2 bohr mag/ cell at the end of calculation.<br>
>> >    nosym            = .true                    /// I hope this is enough<br>
>> to break the symmetry<br>
>> >    nbnd             = 100<br>
>> >    occupations      = 'smearing',  smearing='mp', degauss=0.01D0<br>
>> >  /<br>
>> > &electrons<br>
>> >    diagonalization  = 'cg'<br>
>> >    mixing_beta      = 0.7d00<br>
>> >    conv_thr         = 1.0d-8     // I used even higher convergence since<br>
>> smearing is reduced.<br>
>> > /<br>
>> > ATOMIC_SPECIES<br>
>> >   Ti  47.8670  Ti.pbe-spn-rrkjus_psl.1.0.0.UPF<br>
>> ><br>
>> ><br>
>> > ATOMIC_POSITIONS (crystal)<br>
>> >    Ti     0.5 0.5 0.5<br>
>> ><br>
>> > K_POINTS {automatic}<br>
>> >   8  8  5  0  0  0<br>
>> ><br>
>> ><br>
>> > I don't really know this is the correct way to post a reply. All input<br>
>> is welcome.<br>
>> ><br>
>> > Thanks & Regards,<br>
>> ><br>
>> ------------------------------------------------------------------------------------------------------<br>
>> > Ajmalghan MUTHALI<br>
>> ><br>
>> > Post doctorate researcher<br>
>> > Laboratoire ICB<br>
>> > UMR 6303 CNRS-Université de Bourgogne<br>
>> > 9 Avenue Alain Savary, BP 47870<br>
>> > F-21078 DIJON Cedex, France<br>
>> > Tel: +33-(0)7.69.28.19.91<br>
>> > Email : <a href="mailto:ajmalghan.muthali@u-bourgogne.fr" target="_blank">ajmalghan.muthali@u-bourgogne.fr</a><br>
>> ><br>
>> >_______________________________________________<br>
>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>><br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>