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

登录注册

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

24小时热门版块排行榜

返回列表

joyce79928

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 44
红花: 1
帖子: 64
在线: 23.9小时
虫号: 2823689
注册: 2013-11-24
专业: 材料冶金过程工程

[求助] 结果与文献不相符，是哪里出问题?

我想要做出这篇'Origin Of Current-Induced Forces In An Atomic Gold Wire: A First Principles Study'文献中的结果，力的结果中我得到了相同的趋势，但力的值却与文献相差不少。这会是哪里出了问题?以下是我的输入档:
任何意见都可以提出!!谢谢

Elec:
SystemName  Elec-Au
SystemLabel Elec-Au

# Atomic coordinates
NumberOfAtoms 018
AtomicCoordinatesFormat Ang
%block AtomicCoordinatesAndAtomicSpecies
0.00000000 1.77345822 0.00000000 1
3.07172169 1.77345822 0.00000000 1
6.14344338 1.77345822 0.00000000 1
1.53586084 4.43364706 0.00000000 1
3.07172169 7.09383590 0.00000000 1
4.60758253 4.43364706 0.00000000 1
6.14344338 7.09383590 0.00000000 1
7.67930422 4.43364706 0.00000000 1
9.21516507 7.09383590 0.00000000 1
-1.53586084 0.88672961 2.50127158 1
1.53586084 0.88672961 2.50127158 1
4.60758253 0.88672961 2.50127158 1
0.00000000 3.54691845 2.50127158 1
1.53586084 6.20710728 2.50127158 1
3.07172169 3.54691845 2.50127158 1
4.60758253 6.20710728 2.50127158 1
6.14344338 3.54691845 2.50127158 1
7.67930422 6.20710728 2.50127158 1
%endblock AtomicCoordinatesAndAtomicSpecies

# UNIT CELL
LatticeConstant    1.0 Ang
%block LatticeVectors
9.21529607 0.00000000 0.00000000
4.60764854 7.98068051 0.00000000
0.00000000 0.00000000  4.99155000
%endblock LatticeVectors

# K-points

%block kgrid_Monkhorst_Pack
1 0 0 0.0
0 1 0 0.0
0 0 10 0.5
%endblock kgrid_Monkhorst_Pack

# Chemical species

NumberOfSpecies 1
%block ChemicalSpeciesLabel
  1  79 Au
%endblock ChemicalSpeciesLabel

# Basis set variables
PAO.BasisType split
PAO.BasisSize SZP
PAO.energyshift  0.005  Ry
PAO.BasisType split

# General variables

ElectronicTemperature  300 K
SolutionMethod diagon
MeshCutoff          200. Ry
xc.functional       LDA          # Exchange-correlation functional
xc.authors          CA
SpinPolarized false
      # Exchange-correlation version

# SCF variables

DM.MixSCF1 T
MaxSCFIterations    100          # Maximum number of SCF iter
DM.MixH             T
DM.MixingWeight    0.03       # New DM amount for next SCF cycle
DM.Tolerance       1.d-4       # Tolerance in maximum difference
DM.UseSaveDM       true       # to use continuation files
DM.NumberPulay       6
OccupationFunction    MP

# MD variables

MD.FinalTimeStep 1
MD.TypeOfRun CG
MD.NumCGsteps    000
MD.UseSaveXV    .true.

# Output variables

WriteMullikenPop             1
WriteBands                   .false.
SaveRho                      .false.
SaveDeltaRho                   .false.
SaveHS                         .false.
SaveElectrostaticPotential    no
SaveTotalPotential             no
WriteCoorXmol                .true.
WriteMDXmol                   .true.
WriteMDhistory                .false.
WriteEigenvalues             yes

scatter:
SystemName  BDT-Au-1K
SystemLabel BDT-Au-1K

TS.SaveHS .true.
TS.SaveLead .true.

# ATOMIC POSITIONS
# Atomic coordinates
NumberOfAtoms 84
AtomicCoordinatesFormat Ang
%block AtomicCoordinatesAndAtomicSpecies
0.00000000 0.00000000 0.00000000 1
3.07172169 0.00000000 0.00000000 1
6.14344338 0.00000000 0.00000000 1
1.53586084 2.66018884 0.00000000 1
3.07172169 5.32037767 0.00000000 1
4.60758253 2.66018884 0.00000000 1
6.14344338 5.32037767 0.00000000 1
7.67930422 2.66018884 0.00000000 1
9.21516507 5.32037767 0.00000000 1
0.00000000 1.77345822 2.50127458 1
3.07172169 1.77345822 2.50127458 1
6.14344338 1.77345822 2.50127458 1
1.53586084 4.43364706 2.50127458 1
3.07172169 7.09383590 2.50127458 1
4.60758253 4.43364706 2.50127458 1
6.14344338 7.09383590 2.50127458 1
7.67930422 4.43364706 2.50127458 1
9.21516507 7.09383590 2.50127458 1
-1.53586084 0.88672961 5.00254615 1
1.53586084 0.88672961 5.00254615 1
4.60758253 0.88672961 5.00254615 1
0.00000000 3.54691845 5.00254615 1
1.53586084 6.20710728 5.00254615 1
3.07172169 3.54691845 5.00254615 1
4.60758253 6.20710728 5.00254615 1
6.14344338 3.54691845 5.00254615 1
7.67930422 6.20710728 5.00254615 1
0.00000100 0.00000000 7.49282457 1
3.07172169 0.00000000 7.49282457 1
6.14344338 0.00000000 7.49282457 1
1.53586084 2.66018884 7.49282557 1
3.07172169 5.32037767 7.49282557 1
4.60758253 2.66018884 7.49282557 1
6.14344238 5.32037767 7.49282557 1
7.67930422 2.66018884 7.49282557 1
9.21516507 5.32037767 7.49282557 1
0.00000000 1.77345822 10.00842135 1
3.07172169 1.77345822 10.00842135 1
6.14344338 1.77345822 10.00842135 1
1.53586184 4.43364706 10.00842135 1
3.07172269 7.09383590 10.00842135 1
4.60758353 4.43364706 10.00842135 1
6.14344438 7.09383590 10.00842135 1
7.67930522 4.43364706 10.00842135 1
9.21516507 7.09383590 10.00842135 1
3.07172234 3.54691955 12.42011351 1
3.07172280 3.54691895 15.36850129 1
3.07172241 3.54691961 18.31689056 1
0.00000000 1.77345822 20.72857382 1
3.07172169 1.77345822 20.72857382 1
6.14344338 1.77345822 20.72857382 1
1.53586184 4.43364706 20.72857382 1
3.07172269 7.09383590 20.72857382 1
4.60758353 4.43364706 20.72857382 1
6.14344438 7.09383590 20.72857382 1
7.67930522 4.43364706 20.72857382 1
9.21516507 7.09383590 20.72857382 1
0.00000100 0.00000000 23.24416960 1
3.07172169 0.00000000 23.24416960 1
6.14344338 0.00000000 23.24416960 1
1.53586084 2.66018884 23.24416960 1
3.07172169 5.32037767 23.24416960 1
4.60758253 2.66018884 23.24416960 1
6.14344238 5.32037767 23.24416960 1
7.67930422 2.66018884 23.24416960 1
9.21516507 5.32037767 23.24416960 1
-1.53586084 0.88672961 25.73444802 1
1.53586084 0.88672961 25.73444802 1
4.60758253 0.88672961 25.73444802 1
0.00000000 3.54691845 25.73444802 1
1.53586084 6.20710728 25.73444802 1
3.07172169 3.54691845 25.73444802 1
4.60758253 6.20710728 25.73444802 1
6.14344338 3.54691845 25.73444802 1
7.67930422 6.20710728 25.73444802 1
0.00000000 1.77345822 28.23571960 1
3.07172169 1.77345822 28.23571960 1
6.14344338 1.77345822 28.23571960 1
1.53586084 4.43364706 28.23571960 1
3.07172169 7.09383590 28.23571960 1
4.60758253 4.43364706 28.23571960 1
6.14344338 7.09383590 28.23571960 1
7.67930422 4.43364706 28.23571960 1
9.21516507 7.09383590 28.23571960 1
%endblock AtomicCoordinatesAndAtomicSpecies

# UNIT CELL
LatticeConstant    1.0 Ang
%block LatticeVectors
9.21529607 0.00000000 0.00000000
4.60764854 7.98068051 0.00000000
0.00000000 0.00000000 30.73699118
%endblock LatticeVectors

# K-points

%block kgrid_Monkhorst_Pack
4 0 0 0.0
0 4 0 0.0
0 0 4 0.0
%endblock kgrid_Monkhorst_Pack

# Chemical species

NumberOfSpecies 1
%block ChemicalSpeciesLabel
  1  79 Au
%endblock ChemicalSpeciesLabel

# Basis set variables

PAO.BasisType split
PAO.BasisSize SZP
PAO.BasisType                split          # Type of basis ('nones', 'nonodes', 'split')
PAO.energyshift 0.005  Ry

# General variables

ElectronicTemperature  300 K

SolutionMethod transiesta
MeshCutoff          200. Ry
xc.functional       LDA          # Exchange-correlation functional
xc.authors          CA
SpinPolarized false
      # Exchange-correlation version

# SCF variables

DM.MixSCF1 T
MaxSCFIterations    100          # Maximum number of SCF iter
DM.MixH             T
DM.MixingWeight    0.03       # New DM amount for next SCF cycle
DM.Tolerance       1.d-4       # Tolerance in maximum difference
DM.UseSaveDM       true       # to use continuation files
DM.NumberPulay       5
OccupationFunction    MP

# MD variables

MD.FinalTimeStep 1
MD.TypeOfRun CG
MD.NumCGsteps    000
#%block GeometryConstraints
#position from 1 to 36
#position from 41 to 85
#%endblock GeometryConstraints
MD.UseSaveXV    .true.
UseSaveData .true.

# Output variables

WriteMullikenPop             1
WriteBands                   .false.
SaveRho                      .false.
SaveDeltaRho                   .false.
SaveHS                         .true.
SaveElectrostaticPotential    no
SaveTotalPotential             no
WriteCoorXmol                .true.
WriteMDXmol                   .true.
WriteMDhistory                .false.
WriteEigenvalues             yes

TS.Voltage 2.0 eV
TS.biasContour.NumPoints  100
TS.biasContour.Eta 0.001 eV
SaveTotalPotential T
SaveRho T
SaveElectrostaticPotential T

#------------------------------------
# ELECTRODES
# -----------------------------------
# LEFT ELECTRODE
TS.HSFileLeft '../../elec/4x4/Elec-Au.TSHS'
TS.NumUsedAtomsLeft  18
TS.BufferAtomsLeft  9
# RIGHT ELECTRODE
TS.HSFileRight  '../../elecCBA/4x4/Elec-Au.TSHS'
TS.NumUsedAtomsRight  18

Prof. Brandbyge's result.png

my_result.png

回复此楼

» 猜你喜欢

恳请有学校收留已经有6人回复
求调剂推荐已经有8人回复
0854求调剂已经有14人回复
300求调剂已经有9人回复
291求调剂已经有7人回复
289 分105500药学专硕求调剂(找B区学校) 已经有4人回复
初试324 中药学一志愿天中医求调剂已经有3人回复
药学求调剂已经有14人回复
327求调剂已经有27人回复
急需调剂已经有5人回复

1楼 2014-11-16 20:19:34

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

对照文献合成不出，请求大神指导已经有12人回复
Latex参考文献不同期刊的要求写的格式不一样，怎么办？已经有5人回复
Gaussian新手计算HOMO-LUMO的数值跟文献所给出数值不一样已经有11人回复
Endnote插入参考文献后在word中输出和预览不一致已经有4人回复
两篇文献合成MOF-235方法不一样？已经有13人回复
有一个文献找不到出处！！急！！已经有11人回复
论文的结果和讨论部分引用参考文献这样好不好？已经有29人回复
查不到doi 的文献怎么办？已经有6人回复
一个简单的反应但是按照文献做不出来，求做过类似的大神指教已经有12人回复
新人求助：Elsevier投稿格式中的问题，figure caption，参考文献等，不胜感激啊！已经有22人回复
本人查阅文献，遇到一个问题，请教哪位给我解答哈，在下将不胜感激！已经有18人回复
Elsevier 在线投稿能看到生成的pdf，但是参考文献检查结果不正确，请问是为什么呢已经有6人回复
投的JMC，参考文献的格式严重出错，很丢脸，不知道对后面有影响吗？已经有20人回复
【求助】空间群与文献中的不一样，怎么办？已经有21人回复
【讨论】H-BEA分子筛优化后结果和文献不一样，是怎么回i事？高手指点一下已经有4人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 joyce79928 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 327求调剂 +25	Xxjc1107. 2026-04-13	27/1350	2026-04-15 23:22 by Equinoxhua
[考研] 271求调剂 +32	2261744733 2026-04-11	33/1650	2026-04-15 22:03 by noqvsozv
[考研] 284求调剂 +21	让我上岸吧阿西 2026-04-09	21/1050	2026-04-15 21:53 by noqvsozv
[考研] 085404 22408 309分求调剂 +9	lzmk 2026-04-14	10/500	2026-04-15 20:02 by 学员JpLReM
[考研] 本科211，报考085601-310分 +16	ararak 2026-04-13	16/800	2026-04-14 14:55 by Delta2012
[考研] 求调剂 +12	何气正 2026-04-13	13/650	2026-04-14 14:47 by zs92450
[考研] 085600材料与化工349分求调剂 +16	李木子啊哈哈 2026-04-12	17/850	2026-04-14 09:11 by fenglj492
[考研] 305求调剂 +8	玛卡巴卡boom 2026-04-11	8/400	2026-04-14 09:04 by pengliang8036
[考研] 0856专硕求调剂希望是a区院校 +24	好好休息好不好 2026-04-09	27/1350	2026-04-13 22:22 by pies112
[考研] 材料复试求调剂 +24	xhhdjdjsjks 2026-04-09	24/1200	2026-04-13 15:49 by 幸免 ..
[考研] 求调剂288 +7	ioodiiij 2026-04-10	9/450	2026-04-13 08:33 by Hayaay
[考研] 调剂 +10	只叙离别辞 2026-04-09	12/600	2026-04-11 20:57 by 逆水乘风
[考研] 303求调剂 +14	SereinQ 2026-04-10	15/750	2026-04-11 20:43 by 蓝云思雨
[考研] 283求调剂 +22	那个噜子 2026-04-09	22/1100	2026-04-11 10:41 by 逆水乘风
[考研] 一志愿东北大学控制工程085406数二英二385，求调剂 +8	Ezra_Zhang 2026-04-09	8/400	2026-04-11 09:15 by 猪会飞
[考研] 263能源动力专硕求调剂 +3	加大号饭盒袋 2026-04-10	3/150	2026-04-10 22:23 by 286640313
[考研] 284求调剂 +9	让我上岸吧阿西 2026-04-09	11/550	2026-04-10 19:18 by 靖jing
[考研] 一志愿沪9，326生物学求相关专业调剂 +4	刘墨墨 2026-04-09	4/200	2026-04-10 12:07 by pengliang8036
[考研] 344求调剂 +7	丶风雪夜归人丶 2026-04-09	7/350	2026-04-10 12:05 by pengliang8036
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3	1青科大。 2026-04-09	5/250	2026-04-10 09:58 by 翩翩一书生