版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

joyce79928

銅蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 44
紅花: 1
帖子: 64
在線: 23.9小時
蟲號: 2823689
注冊: 2013-11-24
專業(yè): 材料冶金過程工程

[求助] 結果與文獻不相符，是哪里出問題?

我想要做出這篇'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

回復此樓

» 猜你喜歡

311（085601）求調(diào)劑已經(jīng)有4人回復
305求調(diào)劑已經(jīng)有4人回復
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑已經(jīng)有7人回復
一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑已經(jīng)有12人回復
289求調(diào)劑已經(jīng)有11人回復
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有7人回復
生物學學碩，一志愿湖南大學，初試成績338 已經(jīng)有6人回復
321求調(diào)劑已經(jīng)有7人回復
343求調(diào)劑已經(jīng)有4人回復
本科新能源科學與工程，一志愿華理能動285求調(diào)劑已經(jīng)有6人回復

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

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

對照文獻合成不出，請求大神指導已經(jīng)有12人回復
Latex參考文獻不同期刊的要求寫的格式不一樣，怎么辦？已經(jīng)有5人回復
Gaussian新手計算HOMO-LUMO的數(shù)值跟文獻所給出數(shù)值不一樣已經(jīng)有11人回復
Endnote插入?yún)⒖嘉墨I后在word中輸出和預覽不一致已經(jīng)有4人回復
兩篇文獻合成MOF-235方法不一樣？已經(jīng)有13人回復
有一個文獻找不到出處��！急！！已經(jīng)有11人回復
論文的結果和討論部分引用參考文獻這樣好不好？已經(jīng)有29人回復
查不到doi 的文獻怎么辦？已經(jīng)有6人回復
一個簡單的反應但是按照文獻做不出來，求做過類似的大神指教已經(jīng)有12人回復
新人求助：Elsevier投稿格式中的問題，figure caption，參考文獻等，不勝感激��！已經(jīng)有22人回復
本人查閱文獻，遇到一個問題，請教哪位給我解答哈，在下將不勝感激！已經(jīng)有18人回復
Elsevier 在線投稿能看到生成的pdf，但是參考文獻檢查結果不正確，請問是為什么呢已經(jīng)有6人回復
投的JMC，參考文獻的格式嚴重出錯，很丟臉，不知道對后面有影響嗎？已經(jīng)有20人回復
【求助】空間群與文獻中的不一樣，怎么辦？已經(jīng)有21人回復
【討論】H-BEA分子篩優(yōu)化后結果和文獻不一樣，是怎么回i事？高手指點一下已經(jīng)有4人回復

點擊這里搜索更多相關資源

科研從小木蟲開始，人人為我，我為人人

相關版塊跳轉(zhuǎn) 我要訂閱樓主 joyce79928 的主題更新