| 查看: 3414 | 回复: 12 | |||||
[交流]
【解决】讨论，siesta中作HOMO 与LUMO图的方法
|
|
第一种是用denchar， 但是存在问题是denchar会给出你所有的波函数，需要一个个查看 第二种是用LDOS， 计算在Fermi level 附近的能量的电子分布，转换成*。XSF 查看，但是不能确定具体的LUMO和HOMO的具体位置 大家都用哪种？ [ Last edited by minmin_0082003 on 2011-3-2 at 09:20 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | siesta问题 |
» 猜你喜欢
面上提前没消息，有中的吗
已经有17人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 循环伏安曲线基础问题--氧化峰，还原峰，氧化电位，还原电位及计算HOMO，LUMO的方法
已经有6人回复
- 【求助】如何用高斯计算一个分子中的某一部分HOMO和LUMO？
已经有7人回复
- 【求助】用gaussian view 做HOMO,LUMO图
已经有4人回复
- 【求助】Gaussian计算HOMO-LUMO问题
已经有17人回复
- 【求助】请问，如何在MS中分析HOMO与LUMO，谢谢您的帮助
已经有12人回复
» 抢金币啦！回帖就可以得到:
-
国自然发榜在即：此刻的你，是什么心情？
+2/950
-
今年面上没有消息，想来是没了，散个金攒攒人品，为明年祈福
+1/600
-
2027年博士生（硕博连读/直博）吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/472
-
坐标北京，92男，诚征结婚对象
+1/166
-
威格科技张德龙：二十年深耕，铸就国产惰性气氛设备新高度
+1/88
-
具有超高耐腐蚀性的铝合金阳极氧化后处理工艺
+1/79
-
南京大学医工交叉方向薛璐璐团队招聘脂质纳米材料与RNA递送方向博士后1-2名
+1/77
-
有机合成高级研究员
+1/30
-
中国特种设备检测研究院招聘1名高温涂层方向博士后
+2/26
-
招收博士后（神经退行性疾病、认知障碍疾病）
+1/22
-
高弯曲强度玻璃制品的特殊加工工艺
+1/6
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘（长期有效）
+5/5
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/5
-
北京师范大学（珠海校区）赵亚然副教授招收2027年电催化推免硕士/博士研究生
+1/5
-
北理工集成电路杰青团队 | 诚招科助理
+1/5
-
多本Scopus收录国际英文期刊征稿，全免 APC！含医学期刊
+1/4
-
【2026 PhD/RA】半导体器件和存内计算 （港科广）
+1/3
-
（Scopus期刊）Biofunctional Materials《生物功能性材料》招募青年编委及征稿
+1/3
-
宁波诺丁汉大学新能源与环境材料技术研究中心招收2027年春季入学博士生
+1/2
-
国自然发榜在即：此刻的你，是什么心情？
+1/1
2楼2011-02-24 15:26:57
3楼2011-02-24 18:19:27
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+3): 谢谢交流。 2011-02-24 19:52:03
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+3): 谢谢交流。 2011-02-24 19:52:03
|
用你这个方法也行。不过怎么可能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 手册上有这个，所以要按照手册来好些！我不清楚你的那个写法管不管用，起码我没用过。我这样写的，都可以的。 |
4楼2011-02-24 19:24:17
|
貌似还是不对，一个波函数都没了，如果有时间帮我看一眼，如果我找出错误了，会立刻贴在后面的 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
6楼2011-02-25 09:49:47
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指点 2011-03-01 14:11:30
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指点 2011-03-01 14:11:30
|
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是不行的。 |
7楼2011-03-01 13:41:59
8楼2011-03-01 14:23:49
9楼2011-03-01 14:27:52
10楼2011-03-01 14:52:15
11楼2011-03-07 15:35:21
12楼2013-11-11 16:52:02
13楼2013-11-11 17:09:11
