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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (3828)
>论文投稿 (344)
>虫友互识 (273)
>文献求助 (266)
>学术会议 (220)
>导师招生 (154)
>基金申请 (121)
>考博 (113)
>硕博家园 (111)
>考研 (69)
>博后之家 (60)
>教师之家 (55)
>论文道贺祈福 (52)
>公派出国 (39)
>休闲灌水 (35)
>绿色求助(高悬赏) (30)

返回列表

lionz

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5.2
红花: 1
帖子: 36
在线: 24.2小时
虫号: 637622
注册: 2008-10-27
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】Siesta 原子弛豫很慢已有1人参与

做练习，联系Ni 7层 slab的原子弛豫；放开表面前三层Z方向的坐标；用了16个进程，花了9个小时。请大家帮忙那个参数有问题，input文件如下，谢谢了：
SystemName    NiSurface_LDA  # Descriptive name of the system
SystemLabel    NiSurface          # Short name for naming files

NumberOfSpecies 1
NumberOfAtoms 28
# ------------------------------------------------------------------------------
# Output options
WriteForces    True
WriteCoorXmol    True
WriteMullikenPop  1
WriteCoorStep    True

#-------------------------------------------------------------------------------
# Species and atoms
%block ChemicalSpeciesLabel
1 28 Ni
%endblock ChemicalSpeciesLabel

# Basis
User.Basis       True

#-------------------------------------------------------------------------------
#Structure & K-Sample
LatticeConstant 1.00000  Ang

%block LatticeVectors
4.25386 0.000    0.000
0.000    4.91194 0.000
0.000    0.000    52.0317
%endblock LatticeVectors

%block kgrid_Monkhorst_Pack
  11  0 0  0.0
  0 11  0  0.0
  0 0 1  0.0
%endblock kgrid_Monkhorst_Pack
#----------------------------------------------------------------------------------
# SCF options
MaxSCFIterations    200       # Maximum number of SCF iter
DM.MixingWeight    0.01          # New DM amount for next SCF cycle
#DM.Tolerance       1.d-4       # Tolerance in maximum difference
                                 # between input and output DM
DM.UseSaveDM       true       # to use continuation files
DM.NumberPulay       6

SolutionMethod       diagon       # OrderN or Diagon
Diag.ParallelOverK    True
ElectronicTemperature  25 meV    # Temp. for Fermi smearing

#xc function& spin
xc.functional       LDA          # Exchange-correlation functional
xc.authors          CA          # Exchange-correlation version
SpinPolarized       true       # Logical parameters are: yes or no
MeshCutoff          200. Ry       # Mesh cutoff. real space mesh
#-----------------------------------------------------------------------------------
# MD options

MD.TypeOfRun          cg          # Type of dynamics:
MD.NumCGsteps       60          # Number of CG steps for
                                 # coordinate optimization
#MD.MaxCGDispl       0.1 Ang    # Maximum atomic displacement
                                 # in one CG step (Bohr)
MD.MaxForceTol       0.04 eV/Ang  # Tolerance in the maximum
                                 # atomic force (Ry/Bohr)
#------------------------------------------------------------------------------------
# Atomic coordinates
AtomicCoordinatesFormat    Ang #Fractional
AtomicCoordinatesAndAtomicSpecies < InputNiSurface.xyz

#Relaxation
%block GeometryConstraints
position 1,2,3,6,8,9,11,12,15,16,17,20,22,23,25,26
position 4,5,7,10,13,14,18,19,21,24,27,28  1.0 0.0 0.0
position 4,5,7,10,13,14,18,19,21,24,27,28  0.0 1.0 0.0
%endblock GeometryConstraints

回复此楼

» 猜你喜欢

1楼 2010-07-23 23:08:46

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

hongjw04

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 84
散金: 20
红花: 3
帖子: 58
在线: 40.7小时
虫号: 301365
注册: 2006-11-26
性别: GG
专业: 凝聚态物性 II ：电子结构

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):谢谢回复 2010-07-24 07:52:40

你可以把DM.MixingWeight调大一些， electric Temperture 也可以调大一些，再试试～

赞一下(2人)