24小时热门版块排行榜

返回列表

本帖产生 1 个 1ST强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangguangping

木虫 (著名写手)

1ST强帖: 23
应助: 71 (初中生)
贵宾: 0.031
金币: 1973.1
散金: 9059
红花: 79
帖子: 2681
在线: 2128.8小时
虫号: 529624
注册: 2008-03-20
性别: GG
专业: 原子和分子物理

【答案】应助回帖

franch(1ST强帖+1): 谢谢指点 2011-05-16 20:29:15

引用回帖:

Originally posted by deadbook at 2011-05-16 12:17:33:
这两列数据应该不是自旋吧。因为最后还有一行是Z=1或者z=2，这个应该是和自选相关的，我提取PDOS是用了andrei的代码，不过就是不知道为啥算出来的都是零，相当郁闷啊。

上面的回复不是说明白了吗，z=1或者z=2表示的是这个轨道是第一个zeta.
那两列数据分别代表自旋向上和自旋向下的PDOS.
下面是我随便找了一下体系算的自旋极化的PDOS(红色为批注)
#表示下面的内容是PDOS
2   #表示自旋为2
416 #表示计算中用了416个基函数
#表示下面开始列出能量格点,这个地方有多少点,下面的PDOS也相应有多少点,一一对应
         -3.00000
         -2.93939
         -2.87878
         -2.81818
         -2.75757
         -2.69697
         -2.63636
         -2.57575
         -2.51515
         -2.45454
         -2.39394
         -2.33333
         -2.27272
         -2.21212
         -2.15151
         -2.09091
         -2.03030
         -1.96969
         -1.90909
         -1.84848
         -1.78788
         -1.72727
         -1.66666
         -1.60606
         -1.54545
         -1.48485
         -1.42424
         -1.36363
         -1.30303
         -1.24242
         -1.18182
         -1.12121
         -1.06060
         -1.00000
         -0.93939
         -0.87879
         -0.81818
         -0.75757
         -0.69697
         -0.63636
         -0.57576
         -0.51515
         -0.45454
         -0.39394
         -0.33333
         -0.27273
         -0.21212
         -0.15151
         -0.09091
         -0.03030
         0.03030
         0.09091
         0.15151
         0.21212
         0.27273
         0.33333
         0.39394
         0.45454
         0.51515
         0.57576
         0.63636
         0.69697
         0.75757
         0.81818
         0.87879
         0.93939
         1.00000
         1.06060
         1.12121
         1.18182
         1.24242
         1.30303
         1.36363
         1.42424
         1.48485
         1.54545
         1.60606
         1.66666
         1.72727
         1.78788
         1.84848
         1.90909
         1.96969
         2.03030
         2.09091
         2.15151
         2.21212
         2.27272
         2.33333
         2.39394
         2.45454
         2.51515
         2.57575
         2.63636
         2.69697
         2.75757
         2.81818
         2.87878
         2.93939
         3.00000
   #表示能量点列表完毕
#表示开始列出每一个基函数
index="                      1"       #表示这是第1个基函数
atom_index="                      1" #表示这个基函数是原子序号为1的原子的
species="C"                      #表示第一个原子为C原子
position=" 0.000000 0.000000 -18.831582"  #给出第一个C原子的空间坐标,单位为Bohr
n="                      2" #表示该基函数对应径向函数的主量子数n=2
l="                      0"  #表示该基函数对应角向函数的角量子数l=0
m="                      0" #表示该基函数对应角向函数的磁量子数m=0
z="                      1" #表示该基函数是对应着zeta=1的基函数
>             #基函数信息完毕,总之这个基函数表示C原子的2S轨道价键劈裂基组中的第一个zeta函数
       #开始该基函数的PDOS,分别表示自旋向上和自旋向下的PDOS
0.00082 0.00065
0.00063 0.00064
0.00043 0.00053
0.00029 0.00038
0.00028 0.00027
0.00045 0.00027
0.00080 0.00044
0.00126 0.00078
0.00167 0.00123
0.00186 0.00164
0.00173 0.00184
0.00137 0.00173
0.00094 0.00137
0.00059 0.00094
0.00039 0.00059
0.00029 0.00038
0.00022 0.00028
0.00017 0.00022
0.00011 0.00016
0.00006 0.00011
0.00004 0.00006
0.00003 0.00003
0.00004 0.00002
0.00006 0.00002
0.00007 0.00003
0.00007 0.00005
0.00006 0.00008
0.00003 0.00010
0.00001 0.00011
0.00000 0.00010
  -0.00001 0.00009
  -0.00001 0.00006
  -0.00001 0.00004
  -0.00002 0.00001
  -0.00003  -0.00001
  -0.00003  -0.00002
0.00000  -0.00002
0.00009 0.00002
0.00030 0.00013
0.00067 0.00039
0.00119 0.00084
0.00174 0.00148
0.00213 0.00216
0.00217 0.00263
0.00184 0.00267
0.00131 0.00226
0.00078 0.00160
0.00041 0.00095
0.00022 0.00049
0.00018 0.00024
0.00025 0.00017
0.00038 0.00022
0.00050 0.00033
0.00055 0.00043
0.00051 0.00048
0.00039 0.00045
0.00025 0.00035
0.00014 0.00023
0.00006 0.00013
0.00002 0.00006
0.00001 0.00002
0.00000 0.00001
0.00000 0.00000
0.00000 0.00000
0.00000 0.00000
0.00000 0.00000
0.00000 0.00000
0.00000 0.00000
0.00000 0.00000
0.00001 0.00000
0.00002 0.00001
0.00007 0.00002
0.00023 0.00008
0.00076 0.00029
0.00215 0.00094
0.00518 0.00260
0.01047 0.00609
0.01771 0.01198
0.02503 0.01969
0.02959 0.02705
0.02933 0.03108
0.02459 0.02994
0.01805 0.02442
0.01284 0.01749
0.01066 0.01229
0.01118 0.01034
0.01272 0.01114
0.01366 0.01293
0.01352 0.01411
0.01305 0.01427
0.01324 0.01430
0.01427 0.01533
0.01561 0.01749
0.01742 0.01990
0.02219 0.02234
0.03492 0.02753
0.06085 0.04178
0.10123 0.07209
0.14978 0.11991
0.19316 0.17531
      #该基函数的PDOS信息完毕,下面重复如此格式分别给出各个基函数的PDOS信息

如果想得到某一个原子的PDOS,那么将属于这个原子的基函数上面的PDOS求和即可
还有什么没说明白,如果再看不明白,那我也没招了.