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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【解决】讨论，siesta中作HOMO 与LUMO图的方法
131/1返回列表
查看: 3481  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

minmin_0082003

金虫 (正式写手)

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

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

[ Last edited by minmin_0082003 on 2011-3-2 at 09:20 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 siesta问题

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-02-24 14:30:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★
minmin_0082003(金币+10): 2011-02-24 18:19:44
zzy870720z(金币+1): 谢谢交流 2011-02-25 12:07:08
引用回帖:
Originally posted by minmin_0082003 at 2011-02-24 07:30:31:
第一种是用denchar， 但是存在问题是denchar会给出你所有的波函数，需要一个个查看
第二种是用LDOS， 计算在Fermi level 附近的能量的电子分布，转换成*。XSF 查看，但是不能确定具体的LUMO和HOMO的具体位置
大 ...

谁说denchar会给出所有的波函数？你可以在输入文件中设置要计算那些轨道的。提前数好了，你的HOMO,LUMO是哪几条，不就OK了！
请看一下WaveFuncKPoints的用法。
赞一下(1人) 回复此楼
2楼2011-02-24 15:26:57
已阅   回复此楼   关注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的回帖

zhangguangping

木虫 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+3): 谢谢交流。 2011-02-24 19:52:03
引用回帖:
Originally posted by minmin_0082003 at 2011-02-24 11:19:27:
你说的这个定义，我是已经定义过了
不过也有可能是数错了
我首先是做一下优化，看的*.EIG 里面的能量，因为我只要Gamma点，所以没有考虑其他的，我认为的HOMO是小于费米能级的所有数中最大的一个，而HOMO是大 ...

用你这个方法也行。不过怎么可能HOMO是19，LUMO是23呢？计算一下你的体系里面有多少电子，按照你体系的自旋简并情况，就知道HOMO，LUMO的位置了。
我不清楚你是不是你的格式不对的原因。
%block WaveFuncKPoints
0.000 0.000 0.000 19 23 # Gamma wavefuncs 3 4
%endblock WaveFuncKPoints
我写的话，我会这样写：
%block WaveFuncKPoints
0.000 0.000 0.000 from 19 to 23    # Gamma wavefuncs 3 4
%endblock WaveFuncKPoints

加上一个from to
手册上有这个，所以要按照手册来好些！我不清楚你的那个写法管不管用，起码我没用过。我这样写的，都可以的。
赞一下(2人) 回复此楼
4楼2011-02-24 19:24:17
已阅   回复此楼   关注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

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

还有，如果你不只是找homo跟Lumo的情况下怎么办呢，比如你想找个O2的π轨道，排除一个个筛选的办法
赞一下 回复此楼
6楼2011-02-25 09:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指点 2011-03-01 14:11:30
引用回帖:
Originally posted by minmin_0082003 at 2011-02-25 02:49:47:
还有，如果你不只是找homo跟Lumo的情况下怎么办呢，比如你想找个O2的π轨道，排除一个个筛选的办法

WriteWaveFunctions        T                 
WaveFuncKPointsScale     ReciprocalLatticeVectors
%block WaveFuncKPoints
0.000 0.000 0.000 from    1 to 14     
%endblock WaveFuncKPoints
这个是在siesta计算的fdf文件中设置的。然后计算的结果再用denchar处理。denchar有自己的输入文件。WriteWaveFunctions        T 你是不是忘记了写这一个选项啊。
具体的关于波函数的设置请参照使用手册3.0-b来说是60页。也就是6.13小节。里面有很细节的讲述。另外注意波函数的输出只有用diagon的时候才可以用。OderN是不行的。
赞一下(2人) 回复此楼
7楼2011-03-01 13:41:59
已阅   回复此楼   关注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的回帖

zhangguangping

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
minmin_0082003(金币+10): 谢谢 2011-03-01 14:52:52
引用回帖:
Originally posted by minmin_0082003 at 2011-03-01 07:23:49:
估计问题出在这里吧
WriteWaveFunctions        T
与%block WaveFuncKPoints
0.000 0.000 0.000 from    1 to 14     
%endblock WaveFuncKPoints
一起设置的情况下，会出现所有的波函数。
不设置WriteW ...

这个和版本没有关系。
建议最好WriteWaveFunctions和COOP.write混合用，如果你只想得到一部分波函数的话。
赞一下(1人) 回复此楼
9楼2011-03-01 14:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)
★
zzy870720z(金币+1): 如果问题解决，请帮忙吧标题求助修改为求助成功，谢谢合作 2011-03-01 15:51:10
引用回帖:
Originally posted by zhangguangping at 2011-03-01 14:27:52:
这个和版本没有关系。
建议最好WriteWaveFunctions和COOP.write混合用，如果你只想得到一部分波函数的话。

解决了，不用COOP.Write 就可以得到波函数了，灰常感谢
赞一下(1人) 回复此楼
10楼2011-03-01 14:52:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yashumimi

银虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1): 鼓励交流 2011-03-07 16:52:46
我用的第二种方法 在费米面附近取0。1看电荷密度分布。。。
赞一下(2人) 回复此楼
11楼2011-03-07 15:35:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangtian612

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
2楼: Originally posted by zhangguangping at 2011-02-24 15:26:57
谁说denchar会给出所有的波函数？你可以在输入文件中设置要计算那些轨道的。提前数好了，你的HOMO,LUMO是哪几条，不就OK了！
请看一下WaveFuncKPoints的用法。...

想请问一下这个几条,是怎么数的,比较外行,求详解
赞一下 回复此楼
12楼2013-11-11 16:52:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
12楼: Originally posted by zhangtian612 at 2013-11-11 09:52:02
想请问一下这个几条,是怎么数的,比较外行,求详解...

比如体系200个价电子，如果自旋简并。那么第100条轨道就是HOMO了。101条就是LUMO了。
赞一下(1人) 回复此楼
13楼2013-11-11 17:09:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minmin_0082003 的主题更新
131/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 明天应该可查了！？ +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭