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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】已经修改过很多遍，但是算完后表面结构完全扭曲了，请大家看看错误到底在哪里
52/1返回列表
查看: 881  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁

[交流] 【求助】已经修改过很多遍，但是算完后表面结构完全扭曲了，请大家看看错误到底在哪里

输入文件如下，想计算一个有机分子在Au表面的吸附能，只用了一个Au（111）面，6个Au分子。使用VESTA软件建立的
NumberOfAtoms           21

NumberOfSpecies 5

%block ChemicalSpeciesLabel
        1       79      Au
        2       1       H
        3       6       C
        4       7       N
        5       8       O
%endblock ChemicalSpeciesLabel

PAO.BasisType           split
PAO.BasisSize   DZP
PAO.EnergyShift 200 meV

LatticeConstant   4.0782 Ang

%block Latticeparameters
       1  1  1   90. 90. 90.
%endblock

AtomicCoordinatesFormat Ang
# Default value: NotScaledCartesianBohr
#  * Bohr or NotScaledCartesianBohr
#  * Ang or NotScaledCartesianAng
#  * ScaledCartesian
#  * Fractional or ScaledByLatticeVectors

%block AtomicCoordinatesAndAtomicSpecies
       1.134997585     -0.865427958      2.667642299  1
       1.134997585      2.800524042     -0.998314701 1
       4.800949585     -0.865428958     -0.998314701 1
       1.055172585      1.007996042      0.875067299 1
       3.008423585     -0.945253958      0.875067299 1
       3.008397585      1.007969042     -1.078133701 1
      -1.535468291     -1.667381335     -4.616771325 3
      -2.691536984     -2.022071694     -3.918607325 3
      -0.379399598     -1.312690975     -3.918607325 3
       0.520025508     -1.036740752     -4.461780325 2
      -2.691536984     -2.022071694     -2.522277325 3
      -0.379399598     -1.312690975     -2.522277325  3
      -3.590962091     -2.298021917     -1.979104325 2
      -1.535468291     -1.667381335     -1.824113325 3
      -1.535468291     -1.667381335     -0.737766325  2
      -1.535468291     -1.667381335     -5.703118325 2
       0.895620251     -0.921505529     -1.752278409 4
       1.858789288     -0.625998200     -1.170609729 3
      -3.908601386     -2.395475984     -4.653606289  4
      -4.145020167     -3.594731572     -4.796382267  5
      -4.618939486     -1.486692354     -5.082587506  5
%endblock AtomicCoordinatesAndAtomicSpecies

%block GeometryConstraints
cellside z
%endblock GeometryConstraints

OptialCalculation .true.
Optical.EnergyMinimum 0 Ry
optical.energyMaximum 30 Ry
optical.broaden 0.05ev
Slabdipolecorrection .true.

kgrid_cutoff            30 Ang

XC.functional   GGA
XC.authors      PBE

%block kgrid_Monkhorst_pack
5 0 0 0.0
0 5 0 0.0
0 0 1 0.0
%endblock

MeshCutoff              200 Ry
MaxSCFIterations        100
DM.NumberPulay  3
DM.NumberKick           15
SolutionMethod  diagon
MD.TypeOfRun            CG
MD.NumCGsteps     800
MD.MaxCGDispl           0.2 Ang
MD.MaxForceTol  0.02 eV/Ang
PAO.energy shift 5mRy

[ Last edited by minmin_0082003 on 2010-1-18 at 12:54 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2010-01-18 10:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by yyx19840628 at 2010-1-18 19:05:
我不太懂，来学习学习，楼主要得到什么结果？

要得到吸附能，还没算完呢，不知道对不对
赞一下 回复此楼
13楼2010-01-18 19:09:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

漂泊四方

金虫 (小有名气)
★ ★ ★ ★ ★
minmin_0082003(金币+5,VIP+0):非常感谢。可以指点一下修改的方法吗，如果有时间的话 1-18 10:51
感觉问题比较多。

1. 真空层有取么？
2. Au的层数太少，当然固定衬底的方法也有问题。这个是你表面扭曲的原因。
3. 周期性模型的覆盖度的问题。你这个分子比较大，单胞之间的有机分子在你这个单胞大小的时候肯定有作用。这个是你希望的结果么(SAMs)?
赞一下(3人) 回复此楼
我与我纠缠久
2楼2010-01-18 10:45:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by 漂泊四方 at 2010-1-18 10:45:
感觉问题比较多。

1. 真空层有取么？
2. Au的层数太少，当然固定衬底的方法也有问题。这个是你表面扭曲的原因。
3. 周期性模型的覆盖度的问题。你这个分子比较大，单胞之间的有机分子在你这个单胞大小的时候 ...

没有取真空层，VESTA里面设置真空层的那一项，所以就没取， 老板只让我用一层，说这样可以早点出结果。解决的方法是不是增加层数，或者是固定衬底。如果这个是错误的，该怎么修改才正确呢。对于第三个覆盖度的问题。应该是又三种情况吧，但是本人实在是愚钝，文献看不明白atop之类还有覆盖度是1或者1/3又是该怎么设置。请大侠指点
赞一下 回复此楼
3楼2010-01-18 10:50:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)
★ ★ ★ ★ ★
minmin_0082003(金币+3,VIP+0):谢谢，3层的话估计要耗我好几个月才能做完，siesta没有并行。 1-18 11:26
qasd(金币+2,VIP+0):xiexie~ 1-18 19:45
有机分子太大了，表面太小 取4*4的表面都感觉有点小。
至于真空层 是肯定要取的，大概10A就差不多了。
一层怎么模拟？我看到文献中最少也用三层。
赞一下(12人) 回复此楼
相逢是缘！
4楼2010-01-18 11:17:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭