小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
【求助】LDOS 和PDOS的区别
我想问一下LDOS和PDOS的区别，还有就是在ATK中的two-prob system 中怎么计算PDOS
对，谢谢，我也查了ATK手册，的却不能算PDOS，但是好多用ATK计算的关于双电极系统的文章上都分析了PDOS，感觉不知道怎么算的，有用过ATK的希望指教一下
在ATK的论坛上，有一个叫Nordland的人，写了一个脚本，号称可以计算单个原子的PDOS，很多人就是用的这个脚本来算的PDOS。但事实表明，这个脚本算的PDOS是很不可靠的，甚至是错误的。
可以看一下下面的这个帖子：
/forum/index.php?topic=69.0
不否认，不管是PDOS还是LDOS，都可以分解为单个原子的单个轨道的贡献的总和。但在具体分析的时候，对于LDOS，比较不常见看到将其分解为各个轨道的贡献，尽管可以这么做，这么做也是完全合理的，针对特定目的这种分解可能也是必要的。一般来说，或者对于初学者来说，对于LDOS，最根本应该理解的是它是坐标和能量两者的函数，其中已经包括了所有轨道的贡献。其实LDOS有时候还有另外一个等价的词，叫space resolved density of states，更能说明LDOS的含义.
回答的很赞。通过看前辈的这个帖子解决了困扰了我很久的一个问题。多谢。
那怎样得到这样的LDOS呢？
AF-ZGNR-文献.jpg
同问。如何画出态密度图？
PDOS投影到某个原子态密度 英文是不是指partial DOS,某个原子轨道指的是projected DOS？
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
登录小木虫
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
【求助】lev00 运行出错
进行到以下这一步，读取CHGCAR文件时候出现错误，也查不到什么原因，希望大家帮忙给些建议。
Choose the problem to solve:
============== DOS options =======================
&&Cd. Choose dispersion for the DOS smearing
&&TD. Total DOS: you need only -> band.out, brill.dat
&&PD. Projected DOS calculation
============== DENSITY options ===================
&&D. total/partial electron density, electrostatic potential
SD. total or
electrostatic potential
============== ELECTR. potential options =========
CP. Coulomb potential (point,line,plane)
------------>
....... Choose the file to be read in: ......
& &1. total electron density CHGCAR
& &2. partial electron density PARCHG
& &3. electrostatic potential LOCPOT
& &4. quit: do not read any file
------------>
Incorrect item number! Try again!
forrtl: severe (193): Run-Time Check Failure. The variable 'density_$IJK' is being used without being defined
Image& && && && &&&PC& && && && && & Routine& && && && &Line& && &&&Source
lev00& && && && &&&9F2D&&Unknown& && && && && &Unknown&&Unknown
lev00& && && && &&&8A35&&Unknown& && && && && &Unknown&&Unknown
lev00& && && && &&&BA919&&Unknown& && && && && &Unknown&&Unknown
lev00& && && && &&&89FD&&Unknown& && && && && &Unknown&&Unknown
lev00& && && && &&&93C8&&Unknown& && && && && &Unknown&&Unknown
lev00& && && && &&&AD52C&&density_& && && && && && &103&&density.f90
lev00& && && && &&&89DD&&mainmenu_& && && && && && &89&&mainmenu.f90
lev00& && && && &&&59C9&&MAIN__& && && && && && && &93&&lev00.f90
lev00& && && && &&&324C&&Unknown& && && && && &Unknown&&Unknown
libc.so.6& && && & FCC1D994&&Unknown& && && && && &Unknown&&Unknown
lev00& && && && &&&3159&&Unknown& && && && && &Unknown&&Unknown
按照这个方法，编译仍未成功。make命令之后是$ make
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &atoms.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &menu.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &box.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &kpoints.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &mendeleev.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &code.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &param.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &dos_inc.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &explore.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &siesta_eig.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &device.f90
device.f90(336): remark #7712: This variable has not been used.& &
integer iType(N_at),i,j,k,i1,j1,nat,nt,iErr,N_at
---------------------------------------^
device.f90(509): remark #7712: This variable has not been used.& &
subroutine vector3(DIRC)
-------------------^
device.f90(586): remark #7712: This variable has not been used.& &
subroutine centralP(DIRC)
--------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &bastr.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &do_param.f90
do_param.f90(10): remark #7712: This variable has not been used.& &
character Line*200,cha
-------------------^
do_param.f90(12): remark #7712: This variable has not been used.& &
logical& && && && && && && && && &&&:: found=.false.,Err1=.true.,Err
---------------------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &get_param_siesta.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &hat.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &read_siesta_input.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &read_vasp_geom.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &lev00.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &tools.f90
tools.f90(18): remark #7712: This variable has not been used.& &
&&character answer
------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &tools_strings.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &prep_dos.f90
prep_dos.f90(22): remark #7712: This variable has not been used.& &
integer NKP,item,NB,iErr,ngroup,Num_B
-----------------^
prep_dos.f90(24): remark #7712: This variable has not been used.& &
character cha1*1,cha2*2,answer*1,filen(99)*12,cha,line*16
------------------------^
prep_dos.f90(884): remark #7712: This variable has not been used.& &
character line*80,filen*12,filen1*12,cha*12
-------------------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &dos_add.f90
dos_add.f90(149): remark #7712: This variable has not been used.& &
subroutine phys_band_nmb(E,B_min,B_max,Num_B,N_min,N_max)
---------------------------^
dos_add.f90(482): remark #7712: This variable has not been used.& &
real*8 E(0:3),A(0:3),tot,proj,energy,E21,E0,E10,E20,E30,A10,A20,A30,E1,E31
-----------------------------------------^
dos_add.f90(482): remark #7712: This variable has not been used.& &
real*8 E(0:3),A(0:3),tot,proj,energy,E21,E0,E10,E20,E30,A10,A20,A30,E1,E31
--------------------------------------------------------^
dos_add.f90(482): remark #7712: This variable has not been used.& &
real*8 E(0:3),A(0:3),tot,proj,energy,E21,E0,E10,E20,E30,A10,A20,A30,E1,E31
------------------------------------------------------------^
dos_add.f90(482): remark #7712: This variable has not been used.& &
real*8 E(0:3),A(0:3),tot,proj,energy,E21,E0,E10,E20,E30,A10,A20,A30,E1,E31
----------------------------------------------------------------^
dos_add.f90(482): remark #7712: This variable has not been used.& &
real*8 E(0:3),A(0:3),tot,proj,energy,E21,E0,E10,E20,E30,A10,A20,A30,E1,E31
--------------------------------------------------------------------^
dos_add.f90(483): remark #7712: This variable has not been used.& &
real*8 tot0,tot1,A21,A31,proj0,proj1,B0,Disp,fact,tot2,proj2,fact0,B1,fact1
dos_add.f90(483): remark #7712: This variable has not been used.& &
real*8 tot0,tot1,A21,A31,proj0,proj1,B0,Disp,fact,tot2,proj2,fact0,B1,fact1
-------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &mainmenu.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &read_vasp_psi2.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &plotting.f90
plotting.f90(463): remark #7712: This variable has not been used.& &
& && && && && && && && && && && &&&nclasses,nresol_prv,type_prv)
--------------------------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &prep_disp.f90
prep_disp.f90(543): remark #7712: This variable has not been used.& &
& && && && && && &Nbands0,Nbands1,iFlag,step,Dispers,Broad_Band,iErr)
----------------------------------^
prep_disp.f90(563): remark #7712: This variable has not been used.& &
character line*80,filen*12,filen1*12,cha*12
-------------------------------------^
prep_disp.f90(568): remark #7712: This variable has not been used.& &
real*8 B0_min(NBANDS),B0_max(NBANDS),step,Dispers,Broad_Band,W,a
-------------------------------------------------------------^
prep_disp.f90(570): remark #7712: This variable has not been used.& &
integer iErr,Num_B0,lenght,iFlag,iPh,iPh_min,iPh_max,NB0,NB1,NB,NKP,i,j
-----------------------------------------------------^
prep_disp.f90(570): remark #7712: This variable has not been used.& &
integer iErr,Num_B0,lenght,iFlag,iPh,iPh_min,iPh_max,NB0,NB1,NB,NKP,i,j
---------------------------------------------------------^
prep_disp.f90(570): remark #7712: This variable has not been used.& &
integer iErr,Num_B0,lenght,iFlag,iPh,iPh_min,iPh_max,NB0,NB1,NB,NKP,i,j
-------------------------------------------------------------^
prep_disp.f90(570): remark #7712: This variable has not been used.& &
integer iErr,Num_B0,lenght,iFlag,iPh,iPh_min,iPh_max,NB0,NB1,NB,NKP,i,j
----------------------------------------------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &dipole.f90
dipole.f90(1): remark #7712: This variable has not been used.& &
subroutine dipole(grid,totdens,filen,lenght,Dip)
-----------------------^
dipole.f90(12): remark #7712: This variable has not been used.& &
character iask,cha
----------^
dipole.f90(353): remark #7712: This variable has not been used.& &
real*8 GRID(NGX,NGY,NGZ),R(3),x(3),Dip(3)
-----------------------------------^
dipole.f90(355): remark #7712: This variable has not been used.& &
character iask,cha
----------^
dipole.f90(355): remark #7712: This variable has not been used.& &
character iask,cha
---------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &density.f90
density.f90(20): remark #7712: This variable has not been used.& &
complex*16 :: comp0
--------------^
density.f90(23): remark #7712: This variable has not been used.& &
character answer,cha,iask,cha2*2
----------^
density.f90(23): remark #7712: This variable has not been used.& &
character answer,cha,iask,cha2*2
---------------------^
density.f90(217): remark #7712: This variable has not been used.& &
character iask,cha,line*40,cha1,cha2*2
-------------------^
density.f90(217): remark #7712: This variable has not been used.& &
character iask,cha,line*40,cha1,cha2*2
---------------^
density.f90(718): remark #7712: This variable has not been used.& &
subroutine cut_atoms(grid,totdens)
--------------------------^
density.f90(728): remark #7712: This variable has not been used.& &
real*8 GRID(NGX,NGY,NGZ),R(3),x(3),totdens,drad,rCharge1,rCharge2,rCharge3,charg
------------------------------^
density.f90(732): remark #7712: This variable has not been used.& &
character iask,cha,line*40,answer,cha2*2
---------------^
density.f90(1192): remark #7712: This variable has not been used.& &
real*8 :: vect(3,3),corner(3),da(3),ace1,TotCh=0.0,dV,rCharge
------------------------------------^
density.f90(1194): remark #7712: This variable has not been used.& &
character iask,cha,line*40,cha1,cha2*2
---------------^
density.f90(1194): remark #7712: This variable has not been used.& &
character iask,cha,line*40,cha1,cha2*2
---------------------------^
density.f90(1526): remark #7712: This variable has not been used.& &
integer i,j,k,ix,iy,iz
density.f90(1526): remark #7712: This variable has not been used.& &
integer i,j,k,ix,iy,iz
----------^
density.f90(1526): remark #7712: This variable has not been used.& &
integer i,j,k,ix,iy,iz
------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &write_dens.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &lev_coulmb.f90
lev_coulmb.f90(342): remark #7712: This variable has not been used.& &
real*8 :: Q(NIONS),vMad,factor,gEwald,EPSx,bCENTX,bCENTY,bCENTZ,pot
-------------------------------------------^
lev_coulmb.f90(342): remark #7712: This variable has not been used.& &
real*8 :: Q(NIONS),vMad,factor,gEwald,EPSx,bCENTX,bCENTY,bCENTZ,pot
--------------------------------------------------^
lev_coulmb.f90(342): remark #7712: This variable has not been used.& &
real*8 :: Q(NIONS),vMad,factor,gEwald,EPSx,bCENTX,bCENTY,bCENTZ,pot
---------------------------------------------------------^
lev_coulmb.f90(426): remark #7712: This variable has not been used.& &
real*8 :: R(3),Q(NIONS),tiny=0.00001,dzero=0.0,vMad,factor,gEwald,EPSx
------------------------^
lev_coulmb.f90(429): remark #7712: This variable has not been used.& &
integer lenght,iQuit,item,k2,lenght3,iCheck
-----------------------------^
lev_coulmb.f90(618): remark #7712: This variable has not been used.& &
integer lenght,iQuit,item,k3,k2,k1,lenght3,iCheck
--------------------------------^
lev_coulmb.f90(907): remark #7712: This variable has not been used.& &
real*8 :: Q(NIONS),Pnt(3),a(3),dzero=0.0,EPSx,gEwald,Ew
-------------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &hould.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &invers.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &manip_dens.f90
manip_dens.f90(289): remark #7712: This variable has not been used.& &
integer j0,ijk,iz,iy,ix,NPLWV
------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &plot_add.f90
plot_add.f90(11): remark #7712: This variable has not been used.& &
real*8 GRID(NGX,NGY,NGZ),DIRC(3,3),BCELL(3,3),R(3),x(3),dzero
---------------------------------------------------^
plot_add.f90(212): remark #7712: This variable has not been used.& &
real*8 DIRC(3,3),BCELL(3,3),pA(2),pB(2),pC(2),x(3),GRID(NGX,NGY,NGZ),R(3)
----------------------------------------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &read_density.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &simulate.f90
simulate.f90(606): remark #7712: This variable has not been used.& &
integer n1,n2,n3,ijk,j,i
---------------------^
simulate.f90(634): remark #7712: This variable has not been used.& &
integer ngrid(3),ng(3),ifreq,i,j,ijk,ix,iy,iz,j0,iGrid
simulate.f90(634): remark #7712: This variable has not been used.& &
integer ngrid(3),ng(3),ifreq,i,j,ijk,ix,iy,iz,j0,iGrid
-----------------^
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &read_siesta_pdos.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-c&&-g& &ttag.f90
ifort -debug -check -traceback -fpe0 -warn -Vaxlib&&-g&&-o lev00 atoms.o menu.o&&box.o kpoints.o mendeleev.o code.o param.o dos_inc.o explore.o&&siesta_eig.o device.o bastr.o do_param.o get_param_siesta.o hat.o read_siesta_input.o read_vasp_geom.o lev00.o tools.o tools_strings.o prep_dos.o dos_add.o mainmenu.o read_vasp_psi2.o plotting.o prep_disp.o dipole.o density.o write_dens.o lev_coulmb.o hould.o invers.o manip_dens.o plot_add.o read_density.o simulate.o read_siesta_pdos.o ttag.o&&-L/usr/local/intel/Compiler/11.1/icc/mkl/lib/em64t:/usr/local/intel/Compiler/11.1/ifort/lib/intel64:/usr/local/intel/Compiler/11.1/icc/lib/intel64希望能给予进一步的指导。
我把编译器的设置选项修改了一下，
FCOMPL&&= ifort -debug -check -traceback -fpe0 -warn -Vaxlib $(DEBUG)
FCOMPL&&= ifort&&$(DEBUG)
貌似是没提示编译错误……
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
登录小木虫
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
【求助】siesta 计算DOS应该注意什么问题
输入文件如何设置参数，
DOS 与PDOS还有TDOS的主要区别是什么
做DOS是为了解释说明orbital之外，还有其他用途吗？
请大家指教
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
登录小木虫
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
求教：关于PDOS的处理
如题，使用PAW-PBE泛函，在使用LORBIT=11计算PDOS以后，可以用DOS-split来处理态密度吗？
听有人说，DOS-split只能用来处理超软计算出来的态密度？
如果不能用DOS-split处理，那么有什么小程序可以用来处理PDOS呢？
请各位大侠指教！
:cat39::cat39::cat39::cat39::cat39:
前辈，我好想说的不清楚，
对于出现spd轨道那种，是PDOS还是TDOS？
对于px py pz那种，是PDOS还是TODS？
另外，DOS-split只能用来处理超软计算出来的态密度？
请前辈赐教
态密度（DOS）分析就是将总态密度（TDOS）分解为各个原子的各个轨道（PDOS），它的前提假设是体系中的每个电子状态的轨道是由组成它的各个原子的各个轨道按一定比例混合得到的，态密度分析就是其逆过程。你好像是想分析原子轨道（s、px、py、pz和d轨道）是如何相互作用形成分子轨道或晶体轨道，这些都是通过PDOS来分析的。DOS-split只能用来处理超软计算出来的态密度，为什么，这句话我没听说过。
“DOS-split只能用来处理超软计算出来的态密度，为什么，这句话我没听说过。”
这句话也是别人问我的，吓了一跳，我也没听过……
前辈，我是看到文献里有两种态密度分析，一种是s p d，另一种是
s、px、py、pz……，就是不知道这两种的名称各叫什么，好像都叫混了，感觉……
DOS-split只能用来处理超软计算出来的态密度，我也不确认，需要再查证一下，搞清楚出处和理由。
我觉得你指的两种态密度分析，可能是一种是基于原子轨道的，就是s p d，s px py pz...也是属于这种，px py pz是按取向把p轨道分得更细了；另一种是基于成键类型，如σ、π轨道。
建议读下这篇文献
A chemical and theoretical way to look at bonding on surfaces, Roald Hoffmann, 1988, Reviews of Modern Physics. 把一些概念搞清楚。
首先多謝前輩指教
Roald Hoffmann，一直想讀讀他寫的那本書，但是一直也沒倒出空，唉
我前面好像是說錯了，我一直不太清楚LDOS與PDOS的區別，感覺LDOS就是spd那種，PDOS就是s px py pz那種，不知道這麼理解對不對，請前輩指教
不客气，我做这个也没多久。
那篇Review也不用一次全读完，可以一段一段的理解，并且也不是读一遍就马上全明白的，每隔断时间就要重新拿出来看看，可能会有新的理解。
论坛上有PDOS和LDOS的区别的讨论：http://emuch.net/bbs/viewthread.php?tid=2168825&fpage=1，可以参考下
前輩謙虛了:D
這個帖子我看過，前輩，這段話寫的很好：
我认为PDOS指的是投影或者分波态密度，projected or partial density of states. 意思是投影到轨道的态密度，LDOS指的是局域态密度，可以是单原子的，也可以是某个区域所有原子的，取决于你想分析什么性质。同样的，PDOS也同样可以有LDOS，即每个原子的投影态密度，或者所有原子的投影态密度。
另外小木蟲中所有關於PDOS和LDOS的帖子也看過很多，文獻有時候叫不同的說法，給人說糊塗了……
前輩一定看過Hoffmann的那本書吧，我最近買了，瀏覽了一下，感覺相當不錯
可以这么理解，PDOS偏重于原子轨道的组分，LDOS偏重于空间分布。文献中的不一定是对的，需要自己分辨。
多谢前辈指教，明白了不少
请问Hoffmann的书是什么名字？
固体与表面
不客气，你去小木虫找找，有中文版，也有英文版，找不到联系我
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
登录小木虫
打开微信扫一扫
随时随地聊科研