| 查看: 897 | 回复: 2 | |||
newton3915铁杆木虫 (正式写手)
|
[求助]
siesta 优化石墨结构
|
|
请教各位，小弟初学siesta，尝试着优化石墨的结构，发现程序运行到如下所示没多少就停了，不知什么问题 siesta: ******************** Simulation parameters **************************** siesta: siesta: The following are some of the parameters of the simulation. siesta: A complete list of the parameters used, including default values, siesta: can be found in file out.fdf siesta: redata: Non-Collinear-spin run = F redata: SpinPolarized (Up/Down) run = F redata: Number of spin components = 1 redata: Long output = F redata: Number of Atomic Species = 1 redata: Charge density info will appear in .RHO file --------------------------------------------------------------------------------------------------------------------------- 下面是我的fdf文件，不知道哪里出错了，希望大家点拨一下。不胜感激！ # General System Descriptors SystemName graphite geoopt SystemLabel graphite_opt NumberOfAtoms 4 NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 6 C %endblock ChemicalSpeciesLabel PAO.BasisSize DZP PAO.EnergyShift 0.02 Ry # Structure and K-sampling LatticeConstant 2.46 Ang %block LatticeVectors 1.000 0.000 0.000 -0.500 0.866 0.000 0.000 0.000 2.725 %endblock LatticeVectors AtomicCoordinatesFormat Fractional AtomicCoorFormatOut Ang %block AtomicCoordinatesAndAtomicSpecies 0.000 0.000 0.000 1 0.333 0.667 0.000 1 0.000 0.000 0.500 1 0.667 0.333 0.500 1 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 12 0 0 0.5 0 12 0 0.5 0 0 4 0.5 %endblock kgrid_Monkhorst_Pack # DFT XC.functional LDA XC.authors CA SpinPolarized false MaxSCFIterations 50 DM.MixingWeight 0.25 DM.NumberPulay 3 DM.Tolerance 1.d-4 MeshCutoff 100 Ry SolutionMethod diagon ElectronicTemperature 300.0 K # MD and Relaxations MD.TypeOfRun cg MD.VariableCell true MD.MaxForceTol 0.04 eV/Ang MD.MaxStressTol 1.0 GPa MD.NumCGsteps 200 MD.MaxCGDispl 0.2 Bohr # Output options WriteCoorInitial true WriteCoorStep true WriteForces true WriteKpoints false WriteEigenvalues false Writekbands false WriteBands false WriteMullikenPop false WriteWaveFunction false WriteCoorXmol false WriteCoorCerius false WriteMDCoorXmol false WriteMDhistory false WriteMDXmol false WriteDM true # Options for saving/reading information DM.UseSaveDM false MD.UseSaveXV false MD.UseSaveCG false SaveRho false SaveDeltaRho false SaveElectrostaticPotential false SaveTotalPotential false |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有13人回复
-
关于代码变化问题，想知道的进来
已经有20人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp磁性材料结构优化问题
已经有16人回复
- 静态计算，加U计算不收敛！当U>=7时，就不收敛了，当U小于7时很快收敛！
已经有8人回复
- 请教如何加快siesta的收敛速度
已经有9人回复
- 双层石墨烯的间距多大啊 优化后 谢谢
已经有7人回复
- 晶格弛豫与结构弛豫
已经有18人回复
- 请教，关于Dmol3，如何在优化好的结构中计算能量和频率
已经有4人回复
- VASP NEB 不收敛问题
已经有25人回复
- VASP 计算石墨烯与石墨烯带的区别
已经有11人回复
- Siesta从官网下载赝势计算出错
已经有7人回复
- VASP计算表面吸附 但是能量总是不收敛/。。。
已经有9人回复
- 用Dmol3结构优化时不收敛
已经有5人回复
- 被积函数很快收敛到0，为什么积分不收敛。
已经有3人回复
- 声子谱图（速度自相关函数的傅里叶变换之类的问题）
已经有9人回复
- siesta结构优化收敛标准
已经有5人回复
- 结构优化 siesta 和 vasp
已经有7人回复
- siesta能像vasp一样对于优化后的结构进行振动频率分析吗？
已经有5人回复
- siesta-3.1 计算石墨烯条带不收敛，不知是不是我的输入有问题。
已经有24人回复
- SIESTA计算能级不准怎么解决
已经有10人回复
- 【求助】石墨烯的MS建模
已经有13人回复
- 【求助】单层石墨烯建模后，K点坐标选取?
已经有14人回复
- 【求助】二维材料k点的设置
已经有11人回复
- 【求助】请问用siesta算输运怎么建模啊 我是搞石墨烯
已经有7人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
newton3915
铁杆木虫 (正式写手)
|
顶一下 [ 发自小木虫客户端 ] |
2楼2014-11-24 15:28:40
newton3915
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 4742.3
- 散金: 40
- 红花: 4
- 帖子: 418
- 在线: 391.3小时
- 虫号: 609189
- 注册: 2008-09-22
- 专业: 凝聚态物性 II ：电子结构
3楼2014-11-24 22:59:35
20