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

登录注册

应《网络安全法》要求，自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

回复此楼

» 猜你喜欢

化学专业申博已经有4人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
评审有感已经有27人回复
论文撤稿了已经有7人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有4人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有5人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有4人回复
河北省自然科学基金已经有9人回复

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 的主题更新

返回列表

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

[考博] 化学专业申博 +3	赵子羊 2026-05-23	4/200	2026-05-24 18:10 by 工大学长
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	1rx34o113h 2026-05-23	3/150	2026-05-24 17:41 by 0i3mu4vkjz
[基金申请] 评审有感 +16	popular289 2026-05-18	27/1350	2026-05-24 17:34 by hhs666
[教师之家] 论文撤稿了 +4	bjvtcliu 2026-05-24	7/350	2026-05-24 17:29 by bjvtcliu
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4	hvkbtfonbv 2026-05-23	4/200	2026-05-24 17:21 by 75ui6h7z2t
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	hvkbtfonbv 2026-05-23	3/150	2026-05-24 17:10 by 75ui6h7z2t
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	hvkbtfonbv 2026-05-23	3/150	2026-05-24 17:01 by 75ui6h7z2t
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	a2tycdlnq1 2026-05-23	5/250	2026-05-24 16:21 by hhx1yx9evi
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	a2tycdlnq1 2026-05-23	4/200	2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6	Peterchao 2026-05-18	9/450	2026-05-24 16:02 by 130067131
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！ +4	zju2000 2026-05-21	5/250	2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4	ZXW0202 2026-05-22	9/450	2026-05-24 08:47 by bjvtcliu
[基金申请] 青B发送上会通知了吗 +5	chemBioBro 2026-05-22	7/350	2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3	焦晓明 2026-05-21	3/150	2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助，期刊 +4	希冀，有书读 2026-05-20	8/400	2026-05-22 10:16 by 希冀，有书读
[文学芳草园] 献血感触 +7	呀呀好傻 2026-05-19	13/650	2026-05-21 20:15 by 呀呀好傻
[基金申请] 国自然上会要求 +7	无名者登山 2026-05-18	11/550	2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9	fummck 2026-05-20	10/500	2026-05-21 14:17 by draco1987
[有机交流] 反应很差，大量原料没有反应 5+3	Mr.Zot 2026-05-19	8/400	2026-05-20 22:19 by Equinoxhua
[考博] 博士申请 +5	星…… 2026-05-18	6/300	2026-05-18 23:49 by 糊糊涂涂好