| 4 | 1/1 | 返回列表 |
| 查看: 721 | 回復(fù): 3 | ||
ice_burg新蟲 (初入文壇)
|
[求助]
Abinit結(jié)構(gòu)優(yōu)化時(shí)使用自動(dòng)優(yōu)化給出的不是最優(yōu)解? 已有1人參與
|
|
先謝過 初學(xué)abinit 真心問題多。。。 按照論壇里面臺(tái)灣人的教程 把經(jīng)驗(yàn)值acell c/a帶入optcell=0 ionmov=3 得到一組xred 把xred帶入做optcell=2 ionmov=3 計(jì)算 給輸入文件: ------------------- 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有點(diǎn)奇怪。。。但是還是驗(yàn)證下 都保持不變 給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 ![]() ![]() 不知道怎么個(gè)情況了 一般來說 程序?qū)ぶ挡皇菓?yīng)該比手工靠譜+精確嗎?。。。另外 這個(gè)13.25Bohr*0.407*2的結(jié)果倒是和經(jīng)驗(yàn)值很接近呢。。。 |
新蟲 (初入文壇)
銅蟲 (小有名氣)
銅蟲 (小有名氣)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +9 | Koxui 2026-03-28 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 291求調(diào)劑 +7 | Y-cap 2026-03-29 | 8/400 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|