Dear David <br>You set "starting_magnetization(3)=0.5", but you only have two elements. Therefore , the starting_magnetization of Fe is actually zero. <br>Best wishes.<br>Pang Rui<br><div><sign signid="99"><div style="color:#909090;font-family:Arial Narrow;font-size:12px"><br><br><br><br>------------------</div><div style="font-size:14px;font-family:Verdana;color:#000;"><div><div class="logo" style="width:305px;height:35px;line-height:35px;margin:20px 0 0 0;"><img src="https://exmail.qq.com/cgi-bin/viewfile?type=logo&domain=sustc.edu.cn"></div><div class="c_detail" style="margin:10px 0 0 0;"><h4 class="name" style="margin:0;font-size:14px;font-weight:bold;line-height:28px;zoom:1;">庞瑞(PANG Rui)</h4><p class="position" style="margin:0;line-height:22px;color:#a0a0a0;"></p><p class="department" style="margin:0;line-height:22px;color:#a0a0a0;">South University of Science and Technology of China/Department of Physics</p><p class="department" style="margin:0;line-height:22px;color:#a0a0a0;">No.1088,Xueyuan Road, Shenzhen,Guangdong<br></p><p class="phone" style="margin:0;line-height:22px;color:#a0a0a0;"></p></div></div></div></sign></div><div> </div><div><includetail><div> </div><div> </div><div style="font:Verdana normal 14px;color:#000;"><div style="FONT-SIZE: 12px;FONT-FAMILY: Arial Narrow;padding:2px 0 2px 0;">------------------ Original ------------------</div><div style="FONT-SIZE: 12px;background:#efefef;padding:8px;"><div id="menu_sender"><b>From: </b> "David Foster"<davidfoster751@yahoo.com>;</div><div><b>Date: </b> Fri, Jul 17, 2015 08:44 PM</div><div><b>To: </b> "pw_forum"<pw_forum@pwscf.org>; <wbr></div><div></div><div><b>Subject: </b> [Pw_forum] Magnetization on the Fe atom doped on the graphene</div></div><div> </div>Dear Users<br><br>I have doped one Fe atom in the supercell of graphene, and used "starting_magnetization" keyword (due to the electronic configuration of Fe which is 4s2 3d6).<br><br>Here it is my input:<br><br>=====<br>&CONTROL<br>                       title = 'graph44'<br>                 calculation = 'vc-relax'<br>                restart_mode = 'from_scratch'<br>                      outdir = './graph44_relax'<br>                  pseudo_dir = './'<br>                      prefix = 'graph44'<br>                     disk_io = 'default'<br>                   verbosity = 'default'<br>                   etot_conv_thr=1.0D-6<br>                   forc_conv_thr=1.0D-2<br>                   nstep=1000<br>                   tstress=.true.<br>                   tprnfor=.true.<br>/<br> &SYSTEM<br>                       ibrav = 4<br>                       nat = 32<br>                       celldm(1)=18.783876326<br>                       celldm(3)=2.004008<br>                          ntyp = 2<br>                     ecutwfc = 75<br>                     ecutrho = 500<br>                       starting_magnetization(1)=0.5<br>                       starting_magnetization(3)=0.5<br>                      nspin=2<br>                      occupations='smearing'<br>                      degauss=0.02<br>                      smearing='mv'<br>                       nbnd=160<br>                       <br>/<br> &ELECTRONS<br>            electron_maxstep = 1000<br>                    conv_thr = 1.0D-7<br>                 mixing_mode = 'plain'<br>                 mixing_beta = 0.5<br>                 mixing_ndim = 15<br>             diagonalization = 'david'<br>/<br> &IONS<br>            ion_dynamics = 'bfgs'<br>/<br> &CELL<br>            cell_dynamics = 'bfgs'<br>            cell_dofree='2Dxy'<br>/<br>ATOMIC_SPECIES<br>   C   12.0107         C.pbe-n-rrkjus_psl.1.0.0.UPF<br>   Fe  55.845          Fe.pbe-spn-rrkjus_psl.1.0.0.UPF<br><br>ATOMIC_POSITIONS crystal<br>  C   0.0827257631375200   0.1654515262750400   0.0000000000000000 1 1 0<br>  C   0.1654515262077270   0.0827257630794374   0.0000000000000000 1 1 0<br>  C   0.3309030525741389   0.1654515262750400   0.0000000000000000 1 1 0<br>  C   0.4136288157449489   0.0827257630794371   0.0000000000000000 1 1 0<br>  C   0.5790803421113629   0.1654515262750400   0.0000000000000000 1 1 0<br>  C   0.6618061051815698   0.0827257630794374   0.0000000000000000 1 1 0<br>  C   0.8272576315479819   0.1654515262750410   0.0000000000000000 1 1 0<br>  C   0.9099833946181889   0.0827257630794374   0.0000000000000000 1 1 0<br>  C   0.0827257632048355   0.4136288157456879   0.0000000000000000 1 1 0<br>  C   0.1654515263331270   0.3309030526662500   0.0000000000000000 1 1 0<br>  C   0.3309030525408510   0.4136288157456870   0.0000000000000000 1 1 0<br>  C   0.4136288157697440   0.3309030526662500   0.0000000000000000 1 1 0<br>  C   0.5790803420780759   0.4136288157456880   0.0000000000000000 1 1 0<br>  C   0.6618061053069669   0.3309030526662500   0.0000000000000000 1 1 0<br>  C   0.8272576315146950   0.4136288157456880   0.0000000000000000 1 1 0<br>  C   0.9099833947435888   0.3309030526662499   0.0000000000000000 1 1 0<br>  C   0.0827257632296303   0.6618061053324988   0.0000000000000000 1 1 0<br>  C   0.1654515264004400   0.5790803421368960   0.0000000000000000 1 1 0<br>  C   0.3309030526662500   0.6618061053325008   0.0000000000000000 1 1 0<br>  C   0.4136288157364559   0.5790803421368959   0.0000000000000000 1 1 0<br>  C   0.5790803421028690   0.6618061053324990   0.0000000000000000 1 1 0<br>  C   0.8272576316400939   0.6618061053324990   0.0000000000000000 1 1 0<br>  C   0.9099833947102990   0.5790803421368960   0.0000000000000000 1 1 0<br>  C   0.0827257631382572   0.9099833946869810   0.0000000000000000 1 1 0<br>  C   0.1654515263090670   0.8272576314913760   0.0000000000000000 1 1 0<br>  C   0.3309030526754809   0.9099833946869810   0.0000000000000000 1 1 0<br>  C   0.4136288157456880   0.8272576314913770   0.0000000000000000 1 1 0<br>  C   0.5790803420114971   0.9099833946869810   0.0000000000000000 1 1 0<br>  C   0.6618061051823049   0.8272576314913748   0.0000000000000000 1 1 0<br>  C   0.8272576315487209   0.9099833946869799   0.0000000000000000 1 1 0<br>  C   0.9099833947195301   0.8272576314913750   0.0000000000000000 1 1 0<br> Fe   0.6618061052736799   0.5790803421368970   0.0000000000000000 1 1 0<br>K_POINTS automatic<br>2 2 1 0 0 0<br>====================<br><br><br>and got following scf-converged magnetization in the first step of geometry optimization:<br>========<br><br>Program PWSCF v.5.1.1 starts on 17Jul2015 at 13:41:18 <br><br>     This program is part of the open-source Quantum ESPRESSO suite<br>     for quantum simulation of materials; please cite<br>         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);<br>          URL http://www.quantum-espresso.org", <br>     in publications or presentations arising from this work. More details at<br>     http://www.quantum-espresso.org/quote<br><br>     Parallel version (MPI), running on    12 processors<br>     R & G space division:  proc/nbgrp/npool/nimage =      12<br>     Reading input from graph44.in<br><br>     Current dimensions of program PWSCF are:<br>     Max number of different atomic species (ntypx) = 10<br>     Max number of k-points (npk) =  40000<br>     Max angular momentum in pseudopotentials (lmaxx) =  3<br><br>     Subspace diagonalization in iterative solution of the eigenvalue problem:<br>     scalapack distributed-memory algorithm (size of sub-group:  2*  2 procs)<br><br> <br>     Parallelization info<br>     --------------------<br>     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW<br>     Min        1011     606    159               180960    84085   11235<br>     Max        1013     607    160               180966    84119   11239<br>     Sum       12145    7279   1915              2171547  1009209  134837<br> <br>     Generating pointlists ...<br>     new r_m :   0.0591 (alat units)  1.1102 (a.u.) for type    1<br>     new r_m :   0.0591 (alat units)  1.1102 (a.u.) for type    2<br><br>     Title: <br>     graph44                                                                    <br><br><br>     bravais-lattice index     =            4<br>     lattice parameter (alat)  =      18.7839  a.u.<br>     unit-cell volume          =   11502.3279 (a.u.)^3<br>     number of atoms/cell      =           32<br>     number of atomic types    =            2<br>     number of electrons       =       140.00<br>     number of Kohn-Sham states=          160<br>     kinetic-energy cutoff     =      75.0000  Ry<br>     charge density cutoff     =     500.0000  Ry<br>     convergence threshold     =      1.0E-07<br>     mixing beta               =       0.5000<br>     number of iterations used =           15  plain     mixing<br>     Exchange-correlation      = PBE ( 1  4  3  4 0 0)<br>     nstep                     =         1000<br><br><br>     celldm(1)=  18.783876  celldm(2)=   0.000000  celldm(3)=   2.004008<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   2.004008 )  <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.499000 )  <br><br><br>     PseudoPot. # 1 for  C read from file:<br>     ./C.pbe-n-rrkjus_psl.1.0.0.UPF<br>     MD5 check sum: d965a6b284613baf0982652bb1fc1f03<br>     Pseudo is Ultrasoft + core correction, Zval =  4.0<br>     Generated using "atomic" code by A. Dal Corso  v.5.1.1<br>     Using radial grid of 1073 points,  4 beta functions with: <br>                l(1) =   0<br>                l(2) =   0<br>                l(3) =   1<br>                l(4) =   1<br>     Q(r) pseudized with 0 coefficients <br><br><br>     PseudoPot. # 2 for Fe read from file:<br>     ./Fe.pbe-spn-rrkjus_psl.1.0.0.UPF<br>     MD5 check sum: 12c5fd6419f2a80ca8ce5aff3429efb6<br>     Pseudo is Ultrasoft + core correction, Zval = 16.0<br>     Generated using "atomic" code by A. Dal Corso  v.5.1.1<br>     Using radial grid of 1191 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>        C              4.00    12.01070      C( 1.00)<br>        Fe            16.00    55.84500     Fe( 1.00)<br><br>     Starting magnetic structure <br>     atomic species   magnetization<br>        C            0.500<br>        Fe           0.000<br><br>      2 Sym. Ops. (no inversion) found<br><br><br><br>   Cartesian axes<br><br>     site n.     atom                  positions (alat units)<br>         1           C   tau(   1) = (   0.0000000   0.1432852   0.0000000  )<br>         2           C   tau(   2) = (   0.1240886   0.0716426   0.0000000  )<br>         3           C   tau(   3) = (   0.2481773   0.1432852   0.0000000  )<br>         4           C   tau(   4) = (   0.3722659   0.0716426   0.0000000  )<br>         5           C   tau(   5) = (   0.4963546   0.1432852   0.0000000  )<br>         6           C   tau(   6) = (   0.6204432   0.0716426   0.0000000  )<br>         7           C   tau(   7) = (   0.7445319   0.1432852   0.0000000  )<br>         8           C   tau(   8) = (   0.8686205   0.0716426   0.0000000  )<br>         9           C   tau(   9) = (  -0.1240886   0.3582131   0.0000000  )<br>        10           C   tau(  10) = (   0.0000000   0.2865704   0.0000000  )<br>        11           C   tau(  11) = (   0.1240886   0.3582131   0.0000000  )<br>        12           C   tau(  12) = (   0.2481773   0.2865704   0.0000000  )<br>        13           C   tau(  13) = (   0.3722659   0.3582131   0.0000000  )<br>        14           C   tau(  14) = (   0.4963546   0.2865704   0.0000000  )<br>        15           C   tau(  15) = (   0.6204432   0.3582131   0.0000000  )<br>        16           C   tau(  16) = (   0.7445319   0.2865704   0.0000000  )<br>        17           C   tau(  17) = (  -0.2481773   0.5731409   0.0000000  )<br>        18           C   tau(  18) = (  -0.1240886   0.5014983   0.0000000  )<br>        19           C   tau(  19) = (  -0.0000000   0.5731409   0.0000000  )<br>        20           C   tau(  20) = (   0.1240886   0.5014983   0.0000000  )<br>        21           C   tau(  21) = (   0.2481773   0.5731409   0.0000000  )<br>        22           C   tau(  22) = (   0.4963546   0.5731409   0.0000000  )<br>        23           C   tau(  23) = (   0.6204432   0.5014983   0.0000000  )<br>        24           C   tau(  24) = (  -0.3722659   0.7880687   0.0000000  )<br>        25           C   tau(  25) = (  -0.2481773   0.7164261   0.0000000  )<br>        26           C   tau(  26) = (  -0.1240886   0.7880687   0.0000000  )<br>        27           C   tau(  27) = (  -0.0000000   0.7164261   0.0000000  )<br>        28           C   tau(  28) = (   0.1240886   0.7880687   0.0000000  )<br>        29           C   tau(  29) = (   0.2481773   0.7164261   0.0000000  )<br>        30           C   tau(  30) = (   0.3722659   0.7880687   0.0000000  )<br>        31           C   tau(  31) = (   0.4963546   0.7164261   0.0000000  )<br>        32           Fe  tau(  32) = (   0.3722659   0.5014983   0.0000000  )<br><br>     number of k points=     8  Marzari-Vanderbilt smearing, width (Ry)=  0.0200<br>                       cart. coord. in units 2pi/alat<br>        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.2500000<br>        k(    2) = (   0.0000000  -0.5773503   0.0000000), wk =   0.2500000<br>        k(    3) = (   0.5000000  -0.2886751   0.0000000), wk =   0.2500000<br>        k(    4) = (  -0.5000000  -0.2886751   0.0000000), wk =   0.2500000<br>        k(    5) = (   0.0000000   0.0000000   0.0000000), wk =   0.2500000<br>        k(    6) = (   0.0000000  -0.5773503   0.0000000), wk =   0.2500000<br>        k(    7) = (   0.5000000  -0.2886751   0.0000000), wk =   0.2500000<br>        k(    8) = (  -0.5000000  -0.2886751   0.0000000), wk =   0.2500000<br><br>     Dense  grid:  2171547 G-vectors     FFT dimensions: ( 135, 135, 270)<br><br>     Smooth grid:  1009209 G-vectors     FFT dimensions: ( 108, 108, 216)<br><br>     Largest allocated arrays     est. size (Mb)     dimensions<br>        Kohn-Sham Wavefunctions        25.76 Mb     (   10552,  160)<br>        NL pseudopotentials            42.83 Mb     (   10552,  266)<br>        Each V/rho on FFT grid         12.79 Mb     (  419175,   2)<br>        Each G-vector array             1.38 Mb     (  180962)<br>        G-vector shells                 1.38 Mb     (  180962)<br>     Largest temporary arrays     est. size (Mb)     dimensions<br>        Auxiliary wavefunctions       103.05 Mb     (   10552,  640)<br>        Each subspace H/S matrix        1.56 Mb     (     320,  320)<br>        Each <psi_i|beta_j> matrix      0.65 Mb     (     266,  160)<br>        Arrays for rho mixing          95.94 Mb     (  419175,   15)<br><br>     Check: negative/imaginary core charge=   -0.000001    0.000000<br><br>     Initial potential from superposition of free atoms<br>     Check: negative starting charge=(component1):   -0.000923<br>     Check: negative starting charge=(component2):   -0.000347<br><br>     starting charge  139.99827, renormalised to  140.00000<br><br>     negative rho (up, down):  9.232E-04 3.468E-04<br>     Starting wfc are  134 randomized atomic wfcs +   26 random wfc<br><br>     total cpu time spent up to now is       85.0 secs<br><br>     per-process dynamical memory:   427.7 Mb<br><br>     Self-consistent Calculation<br><br>     iteration #  1     ecut=    75.00 Ry     beta=0.50<br>     Davidson diagonalization with overlap<br>     ethr =  1.00E-02,  avg # of iterations =  7.8<br><br>     negative rho (up, down):  8.979E-04 5.715E-04<br> <br>     Magnetic moment per site:<br>     atom:    1    charge:    1.6553    magn:    0.4208    constr:    0.0000<br>     atom:    2    charge:    1.6564    magn:    0.4234    constr:    0.0000<br>     atom:    3    charge:    1.6879    magn:    0.4233    constr:    0.0000<br>     atom:    4    charge:    1.6476    magn:    0.4301    constr:    0.0000<br>     atom:    5    charge:    1.6901    magn:    0.4219    constr:    0.0000<br>     atom:    6    charge:    1.6540    magn:    0.4268    constr:    0.0000<br>     atom:    7    charge:    1.6913    magn:    0.4182    constr:    0.0000<br>     atom:    8    charge:    1.6145    magn:    0.4197    constr:    0.0000<br>     atom:    9    charge:    1.6534    magn:    0.4214    constr:    0.0000<br>     atom:   10    charge:    1.6897    magn:    0.4144    constr:    0.0000<br>     atom:   11    charge:    1.6909    magn:    0.4336    constr:    0.0000<br>     atom:   12    charge:    1.6847    magn:    0.4380    constr:    0.0000<br>     atom:   13    charge:    1.8452    magn:    0.4414    constr:    0.0000<br>     atom:   14    charge:    1.6811    magn:    0.4409    constr:    0.0000<br>     atom:   15    charge:    1.6933    magn:    0.4303    constr:    0.0000<br>     atom:   16    charge:    1.6533    magn:    0.4169    constr:    0.0000<br>     atom:   17    charge:    1.6536    magn:    0.4270    constr:    0.0000<br>     atom:   18    charge:    1.6926    magn:    0.4206    constr:    0.0000<br>     atom:   19    charge:    1.6876    magn:    0.4265    constr:    0.0000<br>     atom:   20    charge:    1.6853    magn:    0.4360    constr:    0.0000<br>     atom:   21    charge:    1.8448    magn:    0.4443    constr:    0.0000<br>     atom:   22    charge:    1.8439    magn:    0.4420    constr:    0.0000<br>     atom:   23    charge:    1.6483    magn:    0.4396    constr:    0.0000<br>     atom:   24    charge:    1.6135    magn:    0.4298    constr:    0.0000<br>     atom:   25    charge:    1.6880    magn:    0.4233    constr:    0.0000<br>     atom:   26    charge:    1.6520    magn:    0.4145    constr:    0.0000<br>     atom:   27    charge:    1.6883    magn:    0.4244    constr:    0.0000<br>     atom:   28    charge:    1.6521    magn:    0.4257    constr:    0.0000<br>     atom:   29    charge:    1.6838    magn:    0.4368    constr:    0.0000<br>     atom:   30    charge:    1.6563    magn:    0.4340    constr:    0.0000<br>     atom:   31    charge:    1.6501    magn:    0.4374    constr:    0.0000<br>     atom:   32    charge:   11.4685    magn:   -0.0036    constr:    0.0000<br><br>     total cpu time spent up to now is      378.9 secs<br><br>     total energy              =    -630.04480293 Ry<br>     Harris-Foulkes estimate   =    -629.83142751 Ry<br>     estimated scf accuracy    <       5.43017607 Ry<br><br>     total magnetization       =     4.29 Bohr mag/cell<br>     absolute magnetization    =     4.85 Bohr mag/cell<br><br>     iteration #  2     ecut=    75.00 Ry     beta=0.50<br>     Davidson diagonalization with overlap<br>     ethr =  3.88E-03,  avg # of iterations =  5.8<br><br>     negative rho (up, down):  2.473E-03 2.683E-03<br> <br>.<br>.<br>.<br>.<br>.<br>total cpu time spent up to now is     3943.6 secs<br><br>     total energy              =    -632.02402317 Ry<br>     Harris-Foulkes estimate   =    -632.02402327 Ry<br>     estimated scf accuracy    <       0.00000019 Ry<br><br>     total magnetization       =     0.01 Bohr mag/cell<br>     absolute magnetization    =     0.01 Bohr mag/cell<br><br><br>iteration # 25     ecut=    75.00 Ry     beta=0.50<br>     Davidson diagonalization with overlap<br>     ethr =  1.34E-10,  avg # of iterations =  2.4<br><br>     negative rho (up, down):  8.601E-04 8.602E-04<br> <br>     Magnetic moment per site:<br>     atom:    1    charge:    1.6815    magn:   -0.0001    constr:    0.0000<br>     atom:    2    charge:    1.6811    magn:    0.0001    constr:    0.0000<br>     atom:    3    charge:    1.7051    magn:    0.0000    constr:    0.0000<br>     atom:    4    charge:    1.6858    magn:   -0.0000    constr:    0.0000<br>     atom:    5    charge:    1.7065    magn:    0.0000    constr:    0.0000<br>     atom:    6    charge:    1.6828    magn:    0.0001    constr:    0.0000<br>     atom:    7    charge:    1.7076    magn:   -0.0001    constr:    0.0000<br>     atom:    8    charge:    1.6525    magn:    0.0001    constr:    0.0000<br>     atom:    9    charge:    1.6832    magn:   -0.0001    constr:    0.0000<br>     atom:   10    charge:    1.7073    magn:    0.0001    constr:    0.0000<br>     atom:   11    charge:    1.7140    magn:   -0.0001    constr:    0.0000<br>     atom:   12    charge:    1.7107    magn:    0.0001    constr:    0.0000<br>     atom:   13    charge:    1.8439    magn:   -0.0001    constr:    0.0000<br>     atom:   14    charge:    1.7097    magn:    0.0000    constr:    0.0000<br>     atom:   15    charge:    1.7123    magn:   -0.0000    constr:    0.0000<br>     atom:   16    charge:    1.6812    magn:    0.0001    constr:    0.0000<br>     atom:   17    charge:    1.6792    magn:   -0.0000    constr:    0.0000<br>     atom:   18    charge:    1.7101    magn:    0.0001    constr:    0.0000<br>     atom:   19    charge:    1.7064    magn:   -0.0000    constr:    0.0000<br>     atom:   20    charge:    1.7108    magn:    0.0001    constr:    0.0000<br>     atom:   21    charge:    1.8440    magn:   -0.0001    constr:    0.0000<br>     atom:   22    charge:    1.8405    magn:   -0.0001    constr:    0.0000<br>     atom:   23    charge:    1.6823    magn:    0.0001    constr:    0.0000<br>     atom:   24    charge:    1.6504    magn:   -0.0000    constr:    0.0000<br>     atom:   25    charge:    1.7100    magn:    0.0001    constr:    0.0000<br>     atom:   26    charge:    1.6796    magn:    0.0000    constr:    0.0000<br>     atom:   27    charge:    1.7138    magn:    0.0000    constr:    0.0000<br>     atom:   28    charge:    1.6791    magn:   -0.0000    constr:    0.0000<br>     atom:   29    charge:    1.7077    magn:    0.0001    constr:    0.0000<br>     atom:   30    charge:    1.6835    magn:   -0.0001    constr:    0.0000<br>     atom:   31    charge:    1.6800    magn:    0.0001    constr:    0.0000<br>     atom:   32    charge:   11.1637    magn:    0.0012    constr:    0.0000<br><br>     total cpu time spent up to now is     4107.5 secs<br><br>     End of self-consistent calculation<br><br> ------ SPIN UP ------------<br><br><br>          k = 0.0000 0.0000 0.0000 (126127 PWs)   bands (ev):<br><br>===================================================================<br><br><br>we know that for a single electron |μS| = √3 μB, or approximately 1.73 Bohr magnetons. <br><br>For my supercell, I have only one Fe, so, I expected to get magnetization greater than one (for Fe atom, there is 4 unpaired electron). But, I got 0.01 Bohr-magneton. Any help?<br><br><br>Regards<br><br>David Foster<br><br>Ph.D. Student of Chemistry<br><br>_______________________________________________<br>Pw_forum mailing list<br>Pw_forum@pwscf.org<br>http://pwscf.org/mailman/listinfo/pw_forum</div><!--<![endif]--></includetail></div>