| 4 | 1/1 | 返回列表 |
| 查看: 720 | 回復(fù): 3 | ||
ice_burg新蟲 (初入文壇)
|
[求助]
Abinit結(jié)構(gòu)優(yōu)化時使用自動優(yōu)化給出的不是最優(yōu)解? 已有1人參與
|
|
先謝過 初學(xué)abinit 真心問題多。。。 按照論壇里面臺灣人的教程 把經(jīng)驗值acell c/a帶入optcell=0 ionmov=3 得到一組xred 把xred帶入做optcell=2 ionmov=3 計算 給輸入文件: ------------------- dilatmx 1.05 tolmxf 1.0d-05 ionmov 3 optcell 2 ntime 20 ecutsm 0.5 getwfk -1 prtden 0 prteig 0 chksymbreak 0 kptopt 1 ngkpt 6 6 6 toldfe 1.0d-6 nstep 30 diemac 12 ecut 100 enunit 2 ntypat 3 natom 7 znucl 49 80 34 typat 1 1 2 3 3 3 3 acell 10.87 10.87 10.87 rprim 0.5 -0.5 1.02 -0.5 -0.5 1.02 -0.5 -0.5 -1.02 xred 0.0000000000E+00 -5.5511151231E-17 5.5511151231E-17 5.0000000000E-01 2.5000000000E-01 -7.5000000000E-01 0.0000000000E+00 5.0000000000E-01 -5.0000000000E-01 5.1408395117E-01 7.8573180120E-02 -1.2857754903E-01 4.8591604883E-01 6.4266150021E-01 -5.9265713129E-01 5.0004368913E-02 8.7142245097E-01 -3.5733849979E-01 -5.0004368913E-02 4.0734286871E-01 -9.2142681988E-01 ------------------------------------------------------ 給出out結(jié)果: acell 1.3257791925E+01 1.3257791925E+01 1.3257791925E+01 rprim 4.0680634832E-01 -4.0680634832E-01 8.1793471007E-01 -4.0680634832E-01 -4.0680634832E-01 8.1793471007E-01 -4.0680634832E-01 -4.0680634832E-01 -8.1793471007E-01 雖然這組rprim有點奇怪。。。但是還是驗證下 都保持不變 給acell+ 3*0.001 acell1 1.3257791925E+01 1.3257791925E+01 1.3257791925E+01 Bohr acell2 1.3267791925E+01 1.3267791925E+01 1.3267791925E+01 Bohr acell3 1.3277791925E+01 1.3277791925E+01 1.3277791925E+01 Bohr acell4 1.3287791925E+01 1.3287791925E+01 1.3287791925E+01 Bohr acell5 1.3297791925E+01 1.3297791925E+01 1.3297791925E+01 Bohr etotal1 -1.3142002925E+02 etotal2 -1.3142007513E+02 etotal3 -1.3142010670E+02 etotal4 -1.3142012419E+02 etotal5 -1.3142012794E+02 ![]() ![]() 不知道怎么個情況了 一般來說 程序?qū)ぶ挡皇菓?yīng)該比手工靠譜+精確嗎?。。。另外 這個13.25Bohr*0.407*2的結(jié)果倒是和經(jīng)驗值很接近呢。。。 |
新蟲 (初入文壇)
銅蟲 (小有名氣)
銅蟲 (小有名氣)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +7 | 辰沐5211314 2026-03-26 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +6 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +3 | 小張zxy張 2026-03-26 | 5/250 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 317求調(diào)劑 +5 | 十閑wx 2026-03-24 | 5/250 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +6 | 小小甜面團 2026-03-25 | 6/300 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|