版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

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

[交流] 【解决】讨论，siesta中作HOMO 与LUMO图的方法

第一种是用denchar，但是存在问题是denchar会给出你所有的波函数，需要一个个查看
第二种是用LDOS，计算在Fermi level 附近的能量的电子分布，转换成*。XSF 查看，但是不能确定具体的LUMO和HOMO的具体位置
大家都用哪种？

[ Last edited by minmin_0082003 on 2011-3-2 at 09:20 ]

回复此楼

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2011-02-24 14:30:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

Originally posted by zhangguangping at 2011-02-24 15:26:57:
谁说denchar会给出所有的波函数？你可以在输入文件中设置要计算那些轨道的。提前数好了，你的HOMO,LUMO是哪几条，不就OK了！
请看一下WaveFuncKPoints的用法。

你说的这个定义，我是已经定义过了
不过也有可能是数错了
我首先是做一下优化，看的*.EIG 里面的能量，因为我只要Gamma点，所以没有考虑其他的，我认为的HOMO是小于费米能级的所有数中最大的一个，而HOMO是大于费米能级中最小的一个数。然后我就开始数这两个数的位置，假设这两个数排在第19位与23位，那我该这样定义
%block WaveFuncKPoints
0.000 0.000 0.000 19 23 # Gamma wavefuncs 3 4
%endblock WaveFuncKPoints
对吧？嘿嘿

赞一下

回复此楼

3楼2011-02-24 18:19:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

Originally posted by zhangguangping at 2011-02-24 19:24:17:
用你这个方法也行。不过怎么可能HOMO是19，LUMO是23呢？计算一下你的体系里面有多少电子，按照你体系的自旋简并情况，就知道HOMO，LUMO的位置了。
我不清楚你是不是你的格式不对的原因。
%block WaveFuncKPoi ...

貌似还是不对，一个波函数都没了，如果有时间帮我看一眼，如果我找出错误了，会立刻贴在后面的
SystemName H2O       # Descriptive name of the system
SystemLabel  H2O    # Short name for naming files

NumberOfAtoms          3       # Number of atoms
NumberOfSpecies       2          # Number of species

%block Chemical_Species_Label
  1 8 O-wc
  2 1 H_wc_lya
%endblock Chemical_Species_Label

# Lattice, coordinates, k-sampling

AtomicCoordinatesFormat  Ang

%block AtomicCoordinatesAndAtomicSpecies
3.50499483 4.43907340 4.99999972 1    1  O-wc
4.48480041 4.48209231 4.99999999 2    2  H_wc_lya
3.22063236 5.37726685 4.99999999 2    3  H_wc_lya
%endblock AtomicCoordinatesAndAtomicSpecies

LatticeConstant    1.0 Ang
%block LatticeVectors
10.0000000 0.00000000 0.00000000
0.00000000 10.00000 0.00000000
0.00000000 0.00000 10.0
%endblock LatticeVectors

%block WaveFuncKPoints
0.000 0.000 0.000 from 4 to 5
%endblock WaveFuncKPoints

PAO.EnergyShift       10 meV

# Spin

SpinPolarized       .true.

# DFT, Grid, SCF

XC.functional          GGA       # Exchange-correlation functional type
XC.authors             WC       # Particular parametrization of xc func
MeshCutoff             200. Ry    # Equivalent planewave cutoff for the grid
MaxSCFIterations       200       # Maximum number of SCF iterations per step
DM.MixingWeight       0.05       # New DM amount for next SCF cycle
DM.Tolerance          1.d-4    # Tolerance in maximum difference
                                 # between input and output DM
DM.NumberPulay       3          # Number of SCF steps between pulay mixing

Denchar.TypeOfRun 3D
Denchar.PlotCharge .true.
#Denchar.PlotWaveFunction .true.
Denchar.PlaneGeneration threePoints
%block Denchar.Coor3Points
3.50499483 4.43907340 4.99999972
4.48480041 4.48209231 4.99999999
3.22063236 5.37726685 4.99999999
%endblock Denchar.Coor3Points
Denchar.DensityUnits Ele/Ang**3
Denchar.CoorUnits Ang
Denchar.MinX       -20.0  bohr    # Minimum coordinate of the window in X-dir
Denchar.MaxX       30.0  bohr    # Maximum coordinate of the window in X-dir

Denchar.MinY       -20.0  bohr    # Minimum coordinate of the window in Y-dir
Denchar.MaxY       20.0  bohr    # Maximum coordinate of the window in Y-dir

Denchar.MinZ       -20.0  bohr    # Minimum coordinate of the window in Z-dir
Denchar.MaxZ       30.0  bohr    # Maximum coordinate of the window in Z-dir

Denchar.NumberPointsX    80       # Number of points in X-axis
Denchar.NumberPointsY    80       # Number of points in Y-axis
Denchar.NumberPointsZ    80       # Number of points in Z-axis

# Eigenvalue problem: order-N or diagonalization

SolutionMethod       diagon    # OrderN or Diagon
ElectronicTemperature 5 K       # Temp. for Fermi smearing

# Molecular dynamics and relaxations

MD.TypeOfRun          CG       # Type of dynamics:
MD.NumCGsteps          800

# Output options

WriteEigenvalues       .true.
WriteMullikenPop       1          # Write Mulliken Population Analysis
WriteCoorXmol          .true.
COOP.write             .true.
# Options for saving/reading information

SaveRho                .true.    # Write valence pseudocharge at the mesh
SaveDeltaRho          .false.    # Write RHOscf-RHOatm at the mesh
SaveElectrostaticPotential .false. # Write the total elect. pot. at the mesh
SaveTotalPotential    .false.    # Write the total pot. at the mesh
WriteSiestaDim       .true.    # Write minimum dim to siesta.h and stop
WriteDenchar          .true.    # Write information for DENCHAR

赞一下

回复此楼

5楼2011-02-25 09:47:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

还有，如果你不只是找homo跟Lumo的情况下怎么办呢，比如你想找个O2的π轨道，排除一个个筛选的办法

赞一下

回复此楼

6楼2011-02-25 09:49:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

Originally posted by zhangguangping at 2011-03-01 13:41:59:
WriteWaveFunctions       T
WaveFuncKPointsScale    ReciprocalLatticeVectors
%block WaveFuncKPoints
0.000 0.000 0.000 from 1 to 14
%endblock WaveFuncKPoints
这个是在 ...

估计问题出在这里吧
WriteWaveFunctions       T
与%block WaveFuncKPoints
0.000 0.000 0.000 from 1 to 14
%endblock WaveFuncKPoints
一起设置的情况下，会出现所有的波函数。
不设置WriteWaveFunctions       T 只设置与%block WaveFuncKPoints 没有波函数

不知道为什么？难道是code的原因，我用的是siesta3.0-rc，去试试siesta3.0-b 先

赞一下

回复此楼

8楼2011-03-01 14:23:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖