应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助成功】Dmol3计算Al2Fe2的磁性问题，重金答谢
【调剂】北京石油化工学院2024年16个专业接受调剂
141/1返回列表
查看: 3189  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wang0912302

铁虫 (正式写手)

[交流] 【求助成功】Dmol3计算Al2Fe2的磁性问题，重金答谢

各位大虾好，小弟想用Dmol3计算Al2Fe2的磁性，构型已经在高斯里面优化了是最稳定的。那么，我计算时的设置为：Task:Energy
选择Spin unrestricted，Multiplicity是我指定的与高斯计算的自旋多重度是一样的数值。Properties里面选择了Density of states, Electron density, Electrostatics, Orbital, Population analysis。因为是我新手很多不懂，所以计算的时候几乎把所有能选的项目都选择了，一下是输出文件：
Warning: at least C1 SYM file should be there for optgeom
Aluminum     nbas=  1, z=  13, nrfn=  8, rcut=  9.07, e_ref= -0.011956 Ha
                                  rcore=  0.00  zval=  13.00  13.00
   n=1  L=0  occ= 2.00 e=     -55.136017Ha     -1500.3280eV
   n=2  L=0  occ= 2.00 e=      -3.915734Ha      -106.5526eV
   n=2  L=1  occ= 6.00 e=      -2.544879Ha       -69.2497eV
   n=3  L=0  occ= 2.00 e=      -0.273787Ha        -7.4501eV
   n=3  L=1  occ= 1.00 e=      -0.085494Ha        -2.3264eV
   n=3  L=0  occ= 0.00 e=      -0.877651Ha       -23.8821eV
   n=3  L=1  occ= 0.00 e=      -0.638658Ha       -17.3788eV
   n=3  L=2  occ= 0.00 e=      -0.888455Ha       -24.1761eV
Iron         nbas=  2, z=  26, nrfn= 10, rcut=  9.07, e_ref= -0.134172 Ha
                                  rcore=  0.00  zval=  26.00  26.00
   n=1  L=0  occ= 2.00 e=    -254.216093Ha     -6917.5748eV
   n=2  L=0  occ= 2.00 e=     -29.555879Ha      -804.2567eV
   n=2  L=1  occ= 6.00 e=     -25.542745Ha      -695.0538eV
   n=3  L=0  occ= 2.00 e=      -3.351442Ha       -91.1974eV
   n=3  L=1  occ= 6.00 e=      -2.178342Ha       -59.2757eV
   n=3  L=2  occ= 6.00 e=      -0.286130Ha        -7.7860eV
   n=4  L=0  occ= 2.00 e=      -0.189647Ha        -5.1605eV
   n=3  L=2  occ= 0.00 e=      -1.455356Ha       -39.6023eV
   n=4  L=0  occ= 0.00 e=      -0.906806Ha       -24.6755eV
   n=4  L=1  occ= 0.00 e=      -0.660055Ha       -17.9610eV

Symmetry orbitals C1                           
    n  norb    representation
    1    84        a   
  total number of valence orbitals:     84


molecule charge=      0.0   active electron number=      78.0
(without charge=             78.0)
spin polarization=   6.00
electron temperature=     0.050_Ha      1.36_eV    15789._K


Integration in 3d (fixed mesh): ipa=    6  f3= F
npri,iomax,iomin, thres,  rmaxp,    sp
   0    6    1   0.00001   12.00    1.20
Integration points, checksum  0      9712     19468      78.000065

real array elements, matrices vectors etc:       4.0 MB
integer arrays                           :       0.3 MB
min recommended for all-incl workspace   :       6.9 MB
Total memory allocated for arrays         :       8.3 MB

Default spin density on atom Al(  1), L= 1 was used on start: M_s=  1.0
    1  1.000 start spin density on Al(   1), L= 1 M_s=  1.0


            Total Energy           Binding E       Cnvgnce     Time   Iter
Ef        -3006.078159Ha        -0.9945740Ha      4.57E-01     0.0m      1
Ef        -3006.124976Ha        -1.0413906Ha      4.01E-01     0.0m      2
Ef        -3006.016702Ha        -0.9331170Ha      1.60E-01     0.0m      3
Ef        -3006.006076Ha        -0.9224911Ha      9.17E-02     0.0m      4
Ef        -3006.006748Ha        -0.9231631Ha      5.41E-02     0.0m      5
Ef        -3006.006584Ha        -0.9229989Ha      3.34E-02     0.0m      6
Ef        -3006.006836Ha        -0.9232513Ha      2.25E-02     0.0m      7
Ef        -3006.007033Ha        -0.9234480Ha      1.62E-02     0.0m      8
Ef        -3006.007096Ha        -0.9235108Ha      1.16E-02     0.0m      9
Ef        -3006.007063Ha        -0.9234785Ha      7.51E-03     0.0m     10
Ef        -3006.007049Ha        -0.9234636Ha      4.41E-03     0.0m     11
Ef        -3006.007056Ha        -0.9234711Ha      8.83E-04     0.0m     12
Ef        -3006.007055Ha        -0.9234705Ha      6.51E-04     0.0m     13
Ef        -3006.007054Ha        -0.9234690Ha      3.00E-04     0.0m     14
Ef        -3006.007054Ha        -0.9234687Ha      4.52E-04     0.0m     15
Ef        -3006.007053Ha        -0.9234685Ha      2.24E-04     0.0m     16
Ef        -3006.007053Ha        -0.9234683Ha      1.26E-04     0.0m     17
Ef        -3006.007053Ha        -0.9234683Ha      1.02E-04     0.0m     18
Ef        -3006.007053Ha        -0.9234684Ha      5.52E-05     0.0m     19
Ef        -3006.007053Ha        -0.9234684Ha      4.14E-05     0.0m     20
Ef        -3006.007053Ha        -0.9234685Ha      2.43E-05     0.0m     21
Ef        -3006.007053Ha        -0.9234685Ha      1.47E-05     0.0m     22
Ef        -3006.007053Ha        -0.9234685Ha      7.11E-06     0.0m     23
Ef        -3006.007053Ha        -0.9234685Ha      4.69E-06     0.0m     24
Ef        -3006.007053Ha        -0.9234684Ha      2.69E-06     0.0m     25
Ef        -3006.007053Ha        -0.9234684Ha      1.19E-06     0.0m     26
Ef        -3006.007053Ha        -0.9234684Ha      8.11E-07     0.0m     27
Message: SCF converged


Fermi Energy    -0.140616 Ha    -3.826 eV  up spin and    -0.190269 Ha    -5.177 eV  -kTS_e= -0.9201246 Ha down spin
  HOMO is orbital number    80
  LUMO is orbital number    77

    state                         eigenvalue        occupation
                              (au)          (ev)

    1   +   1   a       -254.117580     -6914.894     1.000
    2   +   2   a       -254.117580     -6914.894     1.000
    3   -   1   a       -254.116048     -6914.852     1.000
    4   -   2   a       -254.116048     -6914.852     1.000
    5   +   3   a        -55.131957     -1500.218     1.000
    6   +   4   a        -55.131957     -1500.218     1.000
    7   -   3   a        -55.130267     -1500.172     1.000
    8   -   4   a        -55.130267     -1500.172     1.000
    9   +   5   a        -29.468259      -801.872     1.000
   10   +   6   a        -29.468259      -801.872     1.000
   11   -   5   a        -29.430870      -800.855     1.000
   12   -   6   a        -29.430870      -800.855     1.000
   13   +   7   a        -25.451774      -692.578     1.000
   14   +   8   a        -25.451774      -692.578     1.000
   15   +   9   a        -25.451684      -692.576     1.000
   16   +  10   a        -25.451684      -692.576     1.000
   17   +  11   a        -25.451509      -692.571     1.000
   18   +  12   a        -25.451509      -692.571     1.000
   19   -   7   a        -25.423192      -691.801     1.000
   20   -   8   a        -25.423191      -691.801     1.000
   21   -   9   a        -25.422536      -691.783     1.000
   22   -  10   a        -25.422536      -691.783     1.000
   23   -  11   a        -25.422418      -691.779     1.000
   24   -  12   a        -25.422418      -691.779     1.000
   25   +  13   a         -3.915975      -106.559     1.000
   26   +  14   a         -3.915906      -106.557     1.000
   27   -  13   a         -3.913725      -106.498     1.000
   28   -  14   a         -3.913658      -106.496     1.000
   29   +  15   a         -3.295833       -89.684     1.000
   30   +  16   a         -3.292801       -89.602     1.000
   31   -  15   a         -3.221021       -87.648     1.000
   32   -  16   a         -3.217869       -87.563     1.000
   33   +  17   a         -2.545402       -69.264     1.000
   34   +  18   a         -2.545141       -69.257     1.000
   35   +  19   a         -2.545085       -69.255     1.000
   36   +  20   a         -2.544915       -69.251     1.000
   37   +  21   a         -2.544559       -69.241     1.000
   38   +  22   a         -2.544355       -69.235     1.000
   39   -  17   a         -2.542683       -69.190     1.000
   40   -  18   a         -2.542450       -69.184     1.000
   41   -  19   a         -2.542396       -69.182     1.000
   42   -  20   a         -2.542219       -69.177     1.000
   43   -  21   a         -2.541894       -69.168     1.000
   44   -  22   a         -2.541679       -69.163     1.000
   45   +  23   a         -2.134305       -58.077     1.000
   46   +  24   a         -2.123901       -57.794     1.000
   47   +  25   a         -2.123414       -57.781     1.000
   48   +  26   a         -2.122072       -57.745     1.000
   49   +  27   a         -2.121422       -57.727     1.000
   50   +  28   a         -2.116878       -57.603     1.000
   51   -  23   a         -2.062410       -56.121     1.000
   52   -  24   a         -2.050448       -55.796     1.000
   53   -  25   a         -2.049404       -55.767     1.000
   54   -  26   a         -2.048690       -55.748     1.000
   55   -  27   a         -2.047297       -55.710     1.000
   56   -  28   a         -2.043976       -55.619     1.000
   57   +  29   a         -0.401951       -10.938     0.995
   58   -  29   a         -0.380948       -10.366     0.978
   59   +  30   a         -0.307777        -8.375     0.966
   60   +  31   a         -0.297454        -8.094     0.958
   61   +  32   a         -0.277851        -7.561     0.940
   62   -  30   a         -0.271922        -7.399     0.837
   63   +  33   a         -0.267208        -7.271     0.926
   64   +  34   a         -0.257613        -7.010     0.912
   65   +  35   a         -0.244705        -6.659     0.889
   66   +  36   a         -0.243617        -6.629     0.887
   67   -  31   a         -0.243407        -6.623     0.743
   68   +  37   a         -0.230961        -6.285     0.859
   69   +  38   a         -0.224737        -6.115     0.843
   70   -  32   a         -0.221462        -6.026     0.651
   71   +  39   a         -0.221037        -6.015     0.833
   72   +  40   a         -0.219823        -5.982     0.830
   73   +  41   a         -0.216930        -5.903     0.821
   74   -  33   a         -0.209681        -5.706     0.596
   75   -  34   a         -0.201166        -5.474     0.554
   76   -  35   a         -0.197735        -5.381     0.537
   77   -  36   a         -0.181703        -4.944     0.457
   78   -  37   a         -0.180422        -4.910     0.451
   79   -  38   a         -0.176186        -4.794     0.430
   80   +  42   a         -0.173947        -4.733     0.661
   81   -  39   a         -0.169027        -4.599     0.395
   82   -  40   a         -0.153748        -4.184     0.325
   83   -  41   a         -0.152297        -4.144     0.319
   84   +  43   a         -0.135948        -3.699     0.477
   85   -  42   a         -0.128411        -3.494     0.225
   86   -  43   a         -0.109720        -2.986     0.166
   87   +  44   a         -0.101116        -2.752     0.312
   88   +  45   a         -0.098070        -2.669     0.299
   89   +  46   a         -0.077878        -2.119     0.222
   90   +  47   a         -0.077626        -2.112     0.221
   91   -  44   a         -0.075736        -2.061     0.092
   92   -  45   a         -0.072207        -1.965     0.086
   93   -  46   a         -0.054861        -1.493     0.062
   94   -  47   a         -0.051752        -1.408     0.059
   95   +  48   a         -0.009766        -0.266     0.068
   96   +  49   a          0.013139         0.358     0.044
   97   -  48   a          0.019757         0.538     0.015
   98   -  49   a          0.031766         0.864     0.012
   99   +  50   a          0.095389         2.596     0.009
  100   +  51   a          0.100933         2.747     0.008
  101   -  50   a          0.116502         3.170     0.002
  102   +  52   a          0.122835         3.343     0.005
  103   -  51   a          0.125217         3.407     0.002
  104   +  53   a          0.126953         3.455     0.005
  105   +  54   a          0.140418         3.821     0.004
  106   -  52   a          0.147902         4.025     0.001
  107   -  53   a          0.157117         4.275     0.001
  108   -  54   a          0.159671         4.345     0.001
  109   +  55   a          0.182202         4.958     0.002
  110   +  56   a          0.198994         5.415     0.001
  111   -  55   a          0.203699         5.543     0.000
  113   +  57   a          0.225792         6.144     0.001
  114   +  58   a          0.226573         6.165     0.001
  115   +  59   a          0.228264         6.211     0.001
  116   +  60   a          0.239299         6.512     0.001
  118   +  61   a          0.247413         6.732     0.000

求助成功！！

[ Last edited by zzy870720z on 2011-2-21 at 08:15 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2011-01-20 22:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang0912302

铁虫 (正式写手)

续

Orbital occupation is:
   28  a  ( 1) alpha and  28  a  ( 1) beta, and
    14.00 alpha elec in  28 orbitals plus    8.00  beta elec in  24 orbitals
  Total number electrons:              78.0

Ef        -3006.007053Ha        -0.9234684Ha      9.22E-07     0.0m     28

df  binding energy      -0.9234684Ha         -25.12887eV        -579.497kcal/mol

Energy components:
   Sum of atomic energies =    -3005.0835850Ha
                  Kinetic =        3.0395867Ha
            Electrostatic =       -4.4807708Ha
     Exchange-correlation =        0.2254589Ha
        Spin polarization =        0.2922567Ha

Df  binding energy extrapolated to T=0K      -0.4634061 Ha       -12.60993 eV

+++ Entering Properties Section +++



Charge partitioning by Hirshfeld method: (  0.0001    6.0000)
  Al    1 charge   spin                    0.0121    0.6098
  Al    2 charge   spin                    0.0121    0.6098
  Fe    3 charge   spin                   -0.0121    2.3902
  Fe    4 charge   spin                   -0.0121    2.3902



Density and Electric Field Gradient at nucleus [au]

  Al    1  density at nucleus       1391.831628
          efg:  xx yy zz /  yz xz xy     0.4166   -0.1656   -0.2510
                                        -0.6291    0.0000    0.0000
          principal axis components     -0.8388    0.4166    0.4222
          largest principal component    0.8388
          asymmetry parameter            0.0066
  Al    2  density at nucleus       1391.831628
          efg:  xx yy zz /  yz xz xy     0.4166   -0.1656   -0.2510
                                         0.6291    0.0000    0.0000
          principal axis components     -0.8388    0.4166    0.4222
          largest principal component    0.8388
          asymmetry parameter            0.0066
  Fe    3  density at nucleus      11860.900581
          efg:  xx yy zz /  yz xz xy    -0.7158    0.4693    0.2466
                                         0.0000    0.2302    0.0000
          principal axis components     -0.7680    0.2988    0.4693
          largest principal component    0.7680
          asymmetry parameter            0.2220
  Fe    4  density at nucleus      11860.900580
          efg:  xx yy zz /  yz xz xy    -0.7158    0.4693    0.2466
                                         0.0000   -0.2302    0.0000
          principal axis components     -0.7680    0.2988    0.4693
          largest principal component    0.7680
          asymmetry parameter            0.2220



Spin density at nucleus and Asymmetric hyperfine tensor [au]

  Al    1  spin density at nucleus     0.394085
          aht:  xx yy zz /  yz xz xy     0.0243   -0.0070   -0.0173
                                        -0.0456    0.0000    0.0000
          principal axis components     -0.0580    0.0243    0.0337
  Al    2  spin density at nucleus     0.394085
          aht:  xx yy zz /  yz xz xy     0.0243   -0.0070   -0.0173
                                         0.0456    0.0000    0.0000
          principal axis components     -0.0580    0.0243    0.0337
  Fe    3  spin density at nucleus     1.197554
          aht:  xx yy zz /  yz xz xy     0.0224   -0.0361    0.0137
                                         0.0000   -0.0346    0.0000
          principal axis components     -0.0361   -0.0168    0.0529
  Fe    4  spin density at nucleus     1.197554
          aht:  xx yy zz /  yz xz xy     0.0224   -0.0361    0.0137
                                         0.0000    0.0346    0.0000
          principal axis components     -0.0361   -0.0168    0.0529


Dipole moment vectors (au):
                       x                y                z      
  electronic:       -0.000002         0.000000        -0.440034
  nuclear:           0.000000         0.000000        -0.000053
  net:              -0.000002         0.000000        -0.440087
dipole magnitude:     0.44009 au        1.1186 debye


Mulliken Population analysis

Population analysis for representation  1    a   
              spin up    down    charge    spin
Al(  1) 10 0   1.000   1.000     2.000   0.000
Al(  1) 20 0   1.000   1.000     2.000   0.000
Al(  1) 30 0   0.816   0.682     1.499   0.134
Al(  1) 30 0  -0.006  -0.004    -0.010  -0.002
Al(  1) 21-1   1.000   1.000     2.000   0.000
Al(  1) 31-1   0.313   0.169     0.482   0.144
Al(  1) 31-1  -0.003  -0.003    -0.007   0.000
Al(  1) 21 0   1.000   1.000     2.000   0.000
Al(  1) 31 0   0.343   0.178     0.521   0.166
Al(  1) 31 0  -0.004  -0.002    -0.006  -0.003
Al(  1) 21 1   1.000   1.000     2.000   0.000
Al(  1) 31 1   0.335   0.167     0.501   0.168
Al(  1) 31 1  -0.004  -0.002    -0.006  -0.002
Al(  1) 32-2   0.014   0.009     0.023   0.006
Al(  1) 32-1   0.017   0.011     0.028   0.006
Al(  1) 32 0   0.015   0.011     0.026   0.004
Al(  1) 32 1   0.011   0.006     0.017   0.005
Al(  1) 32 2   0.014   0.010     0.024   0.004
Al(  2) 10 0   1.000   1.000     2.000   0.000
Al(  2) 20 0   1.000   1.000     2.000   0.000
Al(  2) 30 0   0.816   0.682     1.499   0.134
Al(  2) 30 0  -0.006  -0.004    -0.010  -0.002
Al(  2) 21-1   1.000   1.000     2.000   0.000
Al(  2) 31-1   0.313   0.169     0.482   0.144
Al(  2) 31-1  -0.003  -0.003    -0.007   0.000
Al(  2) 21 0   1.000   1.000     2.000   0.000
Al(  2) 31 0   0.343   0.178     0.521   0.166
Al(  2) 31 0  -0.004  -0.002    -0.006  -0.003
Al(  2) 21 1   1.000   1.000     2.000   0.000
Al(  2) 31 1   0.335   0.167     0.501   0.168
Al(  2) 31 1  -0.004  -0.002    -0.006  -0.002
Al(  2) 32-2   0.014   0.009     0.023   0.006
Al(  2) 32-1   0.017   0.011     0.028   0.006
Al(  2) 32 0   0.015   0.011     0.026   0.004
Al(  2) 32 1   0.011   0.006     0.017   0.005
Al(  2) 32 2   0.014   0.010     0.024   0.004
Fe(  3) 10 0   1.000   1.000     2.000   0.000
Fe(  3) 20 0   1.000   1.000     2.000   0.000
Fe(  3) 30 0   1.000   1.001     2.001  -0.001
Fe(  3) 40 0   0.563   0.338     0.901   0.226
Fe(  3) 40 0   0.000   0.000     0.000   0.001
Fe(  3) 21-1   1.000   1.000     2.000   0.000
Fe(  3) 31-1   1.000   1.001     2.001   0.000
Fe(  3) 41-1   0.079   0.040     0.119   0.039
Fe(  3) 21 0   1.000   1.000     2.000   0.000
Fe(  3) 31 0   1.001   1.001     2.002   0.000
Fe(  3) 41 0   0.080   0.040     0.120   0.040
Fe(  3) 21 1   1.000   1.000     2.000   0.000
Fe(  3) 31 1   1.001   1.001     2.002  -0.001
Fe(  3) 41 1   0.075   0.034     0.109   0.041
Fe(  3) 32-2   0.863   0.453     1.316   0.410
Fe(  3) 32-2   0.005   0.006     0.011  -0.001
Fe(  3) 32-1   0.865   0.452     1.318   0.413
Fe(  3) 32-1   0.004   0.005     0.009  -0.001
Fe(  3) 32 0   0.866   0.458     1.324   0.408
Fe(  3) 32 0   0.004   0.005     0.008  -0.001
Fe(  3) 32 1   0.862   0.455     1.317   0.406
Fe(  3) 32 1   0.005   0.006     0.011  -0.002
Fe(  3) 32 2   0.862   0.469     1.331   0.393
Fe(  3) 32 2   0.002   0.003     0.006  -0.001
Fe(  4) 10 0   1.000   1.000     2.000   0.000
Fe(  4) 20 0   1.000   1.000     2.000   0.000
Fe(  4) 30 0   1.000   1.001     2.001  -0.001
Fe(  4) 40 0   0.563   0.338     0.901   0.226
Fe(  4) 40 0   0.000   0.000     0.000   0.001
Fe(  4) 21-1   1.000   1.000     2.000   0.000
Fe(  4) 31-1   1.000   1.001     2.001   0.000
Fe(  4) 41-1   0.079   0.040     0.119   0.039
Fe(  4) 21 0   1.000   1.000     2.000   0.000
Fe(  4) 31 0   1.001   1.001     2.002   0.000
Fe(  4) 41 0   0.080   0.040     0.120   0.040
Fe(  4) 21 1   1.000   1.000     2.000   0.000
Fe(  4) 31 1   1.001   1.001     2.002  -0.001
Fe(  4) 41 1   0.075   0.034     0.109   0.041
Fe(  4) 32-2   0.863   0.453     1.316   0.410
Fe(  4) 32-2   0.005   0.006     0.011  -0.001
Fe(  4) 32-1   0.865   0.452     1.318   0.413
Fe(  4) 32-1   0.004   0.005     0.009  -0.001
Fe(  4) 32 0   0.866   0.458     1.324   0.408
Fe(  4) 32 0   0.004   0.005     0.008  -0.001
Fe(  4) 32 1   0.862   0.455     1.317   0.406
Fe(  4) 32 1   0.005   0.006     0.011  -0.002
Fe(  4) 32 2   0.862   0.469     1.331   0.393
Fe(  4) 32 2   0.002   0.003     0.006  -0.001


summarized population analysis
              spin up    down    charge    spin
Al(  1) 10     1.000   1.000     2.000   0.000
Al(  1) 20     1.000   1.000     2.000   0.000
Al(  1) 30     0.811   0.678     1.489   0.132
Al(  1) 21     3.000   3.000     6.000   0.000
Al(  1) 31     0.980   0.507     1.487   0.472
Al(  1) 32     0.072   0.046     0.118   0.025
Al(  2) 10     1.000   1.000     2.000   0.000
Al(  2) 20     1.000   1.000     2.000   0.000
Al(  2) 30     0.811   0.678     1.489   0.132
Al(  2) 21     3.000   3.000     6.000   0.000
Al(  2) 31     0.980   0.507     1.487   0.472
Al(  2) 32     0.072   0.046     0.118   0.025
Fe(  3) 10     1.000   1.000     2.000   0.000
Fe(  3) 20     1.000   1.000     2.000   0.000
Fe(  3) 30     1.000   1.001     2.001  -0.001
Fe(  3) 40     0.564   0.338     0.901   0.226
Fe(  3) 21     3.000   3.000     6.000   0.000
Fe(  3) 31     3.002   3.003     6.005  -0.001
Fe(  3) 41     0.234   0.114     0.348   0.120
Fe(  3) 32     4.338   2.313     6.651   2.025
Fe(  4) 10     1.000   1.000     2.000   0.000
Fe(  4) 20     1.000   1.000     2.000   0.000
Fe(  4) 30     1.000   1.001     2.001  -0.001
Fe(  4) 40     0.564   0.338     0.901   0.226
Fe(  4) 21     3.000   3.000     6.000   0.000
Fe(  4) 31     3.002   3.003     6.005  -0.001
Fe(  4) 41     0.234   0.114     0.348   0.120
Fe(  4) 32     4.338   2.313     6.651   2.025

Mulliken atomic charges:
          charge    spin
Al(  1)  -0.094   0.630
Al(  2)  -0.094   0.630
Fe(  3)   0.094   2.370
Fe(  4)   0.094   2.370

Mulliken bond orders
    Al  1    Al  2    0.1593
    Al  1    Fe  3    0.1010
    Al  2    Fe  3    0.1010
    Al  1    Fe  4    0.1010
    Al  2    Fe  4    0.1010
    Fe  3    Fe  4    0.1871

Mayer bond orders
    Al  1    Al  2    0.4455
    Al  1    Fe  3    0.3473
    Al  2    Fe  3    0.3473
    Al  1    Fe  4    0.3473
    Al  2    Fe  4    0.3473
    Fe  3    Fe  4    0.5545

Mayer total valence
  Al    1    1.2496
  Al    2    1.2496
  Fe    3    2.1450
  Fe    4    2.1450

Mayer free valence
  Al    1    0.1096
  Al    2    0.1096
  Fe    3    0.8960
  Fe    4    0.8960


Plotting output:
  property:                                                file name:
  DMol3 total electron density                             Al2Fe2_density.grd            
  DMol3 deformation electron density                       Al2Fe2_deformation_density.grd
  DMol3 spin density                                       Al2Fe2_spin_density.grd      
  DMol3 electrostatic potential Ha*electron(-1)            Al2Fe2_potential.grd         
  orbital   80  42 +a     -0.1739  0.66   Al2Fe2_homo.grd               
  orbital   77  36 -a     -0.1817  0.46   Al2Fe2_lumo.grd               

grid specifications: I_dim, Origin, n_intervals to corner:
    3   -7.6534   -7.9368   -6.8030   52   -7.6534   -7.9368    7.9368
   57   -7.6534    8.2203   -6.8030   54    7.6534   -7.9368   -6.8030
我现在困惑的是找不到具体的每个原子的磁矩大小？难道就是各个spin up 和spin down之和吗？我还需要各个原子的s,p,d轨道的磁矩大小，我看了很多用Dmol3做的paper都给出了原子具体的s,p,d轨道的磁矩大小数值，请问如何从输出文件里面读取呢？？
小弟万分感谢，相助者，重金答谢！小弟感激不尽！！
赞一下 回复此楼
2楼2011-01-20 22:29:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
★
zzy870720z(金币+1): 谢谢建议 2011-01-22 18:56:17
wang0912302(金币+10): 2011-01-25 15:28:26
建议楼住不要用 Dmol3 算磁性 ＶＡＳＰ不是更好吗
赞一下(1人) 回复此楼
3楼2011-01-22 17:46:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang0912302

铁虫 (正式写手)

我算的是cluster

★
zzy870720z(金币+1): 如果问题解决，请帮忙把标题求助修改为求助成功，修改后归还金币，谢谢合作 2011-02-03 18:01:37
引用回帖:
Originally posted by jghe at 2011-01-22 17:46:51:
建议楼住不要用 Dmol3 算磁性 ＶＡＳＰ不是更好吗

我算的是cluster，涉及非周期性的，用VASP的话要建立晶包，必须考虑两个层之间的影响，同时VASP算出来的磁性是所有s,p,d轨道的磁性之和，没有具体到每个轨道的磁性大小。当然，用Dmol3算出来的也是大致的spin up spin down，磁性大小直接于此相关，最近JCP发的几篇cluster的都是用Dmol3算的。
谢谢啊！我的问题已经解决了！黑嘿嘿嘿嘿
赞一下(1人) 回复此楼
4楼2011-01-25 15:31:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chdx0618

木虫 (小有名气)

建议

★
zzy870720z(金币+1): 鼓励交流 2011-02-03 18:01:50
建议以后楼主把生成的文件上传！
赞一下(1人) 回复此楼
5楼2011-02-01 20:20:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yasistar

禁虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
本帖内容被屏蔽
6楼2011-04-14 09:15:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdsdlgh

银虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
489017楼: Originally posted by wang0912302 at 2011-01-20 22:29:31
Orbital occupation is:
   28  a  ( 1) alpha and  28  a  ( 1) beta, and
    14.00 alpha elec in  28 orbitals plus    8.00  beta elec in  24 orbitals
  Total number electrons:              78.0

...

楼主怎么看的。倒是说说解决办法啊
赞一下 回复此楼
7楼2012-10-29 19:31:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangml_2000

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主问题解决了，恭喜恭喜！
现在我也遇到同样的问题，能不能把解决方案详细说明一下，本人感激不尽
赞一下 回复此楼
8楼2013-03-01 11:47:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangml_2000

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主怎么不回话呢？
有问题解决了是好事，但是大家都在问你问题了，忙得没有看见吗？
赞一下 回复此楼
9楼2013-03-05 09:25:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang0912302

铁虫 (正式写手)
引用回帖:
1631530楼: Originally posted by zhangml_2000 at 2013-03-01 11:47:08
楼主问题解决了，恭喜恭喜！
现在我也遇到同样的问题，能不能把解决方案详细说明一下，本人感激不尽

对不起，我最近忙着签证的事情，，，这个东西是很久很久之前做的东西，我没印象了，等我回头忙好了，找找之前的数据和资料在给你详细解答吧，抱歉，，，，
赞一下 回复此楼
10楼2013-03-05 19:50:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang0912302

铁虫 (正式写手)
引用回帖:
1631798楼: Originally posted by zhangml_2000 at 2013-03-05 09:25:04
楼主怎么不回话呢？
有问题解决了是好事，但是大家都在问你问题了，忙得没有看见吗？

对不起，我最近忙着签证的事情，，，这个东西是很久很久之前做的东西，我没印象了，等我回头忙好了，找找之前的数据和资料在给你详细解答吧，抱歉，，，，
赞一下 回复此楼
11楼2013-03-05 19:50:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangml_2000

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
哦，原来是这样！你忙完了帮忙看看！谢谢
赞一下 回复此楼
12楼2013-03-06 08:32:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_shuangyan

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
4楼: Originally posted by wang0912302 at 2011-01-25 15:31:58
我算的是cluster，涉及非周期性的，用VASP的话要建立晶包，必须考虑两个层之间的影响，同时VASP算出来的磁性是所有s,p,d轨道的磁性之和，没有具体到每个轨道的磁性大小。当然，用Dmol3算出来的也是大致的spin up  ...

向你打听一下，用castep计算超胞的磁性时，每个原子都要加U吗？
赞一下 回复此楼
13楼2013-05-06 20:48:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

醒神星期一

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主，这么多年过去了，这个问题还是困扰着很多人，如果你还活跃在小木虫，能不能麻烦你透露一下具体的方法？
赞一下 回复此楼
14楼2018-07-28 08:35:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wang0912302 的主题更新
141/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求 +3 zju2000 2024-04-18 3/150 2024-04-20 00:48 by cococat951
[硕博家园] 哎 +5 s150535912 2024-04-18 6/300 2024-04-19 20:33 by s150535912
[考博] 申请24博士 材料/化工/环境 +4 满目_星辰 2024-04-17 4/200 2024-04-19 20:10 by 前行的道路
[找工作] 事业单位还是大学好？ +17 青萍之沫 2024-04-16 18/900 2024-04-19 17:45 by charles-c
[教师之家] 美国教授开除中国研究生：我就不该录取你！ +7 苏东坡二世 2024-04-15 7/350 2024-04-19 16:20 by feng6531
[基金申请] 院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度 +7 zju2000 2024-04-16 10/500 2024-04-19 11:10 by dwuab
[论文投稿] 国家重点研发计划论文可以发MDPI里的期刊吗 +4 lin_2018 2024-04-17 4/200 2024-04-19 08:46 by TsolmonZ
[考研] 267求调剂 +3 工科材料267 2024-04-15 3/150 2024-04-18 19:45 by lature00
[博后之家] 博后换方向可行吗？ +3 越越不暴躁 2024-04-15 3/150 2024-04-18 10:58 by ciompman
[考研] 一志愿北化085600 328求线上调剂 +10 Mrshark 2024-04-13 10/500 2024-04-18 08:30 by 刘国宁
[考研] 274求调剂 +7 四个大字，你滴� 2024-04-14 10/500 2024-04-17 22:38 by 华北刘兵
[论文投稿] MDPI初审7天，是不是被拒 20+6 applegirl5762 2024-04-15 8/400 2024-04-17 08:49 by bobvan
[基金申请] 请问教育部人文社科难度大吗 +9 锦衣卫寒战 2024-04-15 14/700 2024-04-16 23:10 by 锦衣卫寒战
[考研] 294求调剂 +3 694062003 2024-04-15 4/200 2024-04-16 15:01 by 邹邹哈哈
[考研] 329求调剂 +6 Kaylawander 2024-04-13 7/350 2024-04-16 12:00 by 风来花开1
[考研] 273求调剂 +5 Late婉安 2024-04-15 7/350 2024-04-15 13:01 by Late婉安
[考研] 专硕调剂招生 +3 电致发光 2024-04-15 4/200 2024-04-15 07:34 by ashorewmj
[考研] 290求调剂 +3 杨yhr 2024-04-14 5/250 2024-04-14 21:50 by coco1981
[考研] 086000生医没考数学304分，一志愿西农求调剂 +3 王有田有学上 2024-04-13 5/250 2024-04-13 16:18 by coco1981
[考博] 现在申博很晚了吗 +3 熙风慕寒 2024-04-13 3/150 2024-04-13 11:08 by Nanamiwww
信息提示
关闭
请填处理意见
关闭