 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 454 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問(wèn)題，作者清風(fēng)一郎將贈(zèng)送您 10 個(gè)金幣 | |||
清風(fēng)一郎銅蟲(chóng) (初入文壇)
|
[求助]
QE計(jì)算分子晶體的拉曼光譜，運(yùn)行ph.in報(bào)錯(cuò) 已有1人參與
|
||
|
各位老師好！ 我在用quantumespresso計(jì)算一個(gè)分子晶體的拉曼光譜，在運(yùn)行ph.in文件時(shí)，運(yùn)行一段時(shí)間后，遇到如下報(bào)錯(cuò)，請(qǐng)問(wèn)可能是什么原因引起的呢？ error in routine d2ionq_dispd3 (1): the hessian file: ./tmp/cbp-alpha.hess is missing. 以下為ph.in的輸入文件 &inputph tr2_ph = 1.0d-14 prefix = 'cbp-alpha' outdir = './tmp' fildyn = 'cbp-alpha' ldisp = .false. epsil = .true. trans = .true. qplot = .false. nq1 = 1, nq2 = 1, nq3 = 1 / 0.0 0.0 0.0 以下為上一步scf計(jì)算的pw.in輸入文件 &control title='cbp-alpha-scf', calculation='scf', prefix = 'cbp-alpha' pseudo_dir='./', outdir='./tmp', verbosity='high', / &system ibrav= 0, nat= 124, ntyp= 3, occupations = 'fixed', ecutwfc = 80, ecutrho = 800, vdw_corr = 'grimme-d3', / &electrons conv_thr = 1.0d-10 mixing_beta = 0.5d0 / &ions / &cell / atomic_species c 12.0107 c.pbe-n-rrkjus_psl.1.0.0.upf h 1.00794 h.pbe-rrkjus_psl.1.0.0.upf n 14.00674 n.pbe-n-rrkjus_psl.1.0.0.upf cell_parameters (angstrom) 8.459549526 0.000000000 -0.013055016 0.000000000 15.662496878 0.000000000 -3.735405692 0.000000000 9.404073167 atomic_positions (crystal) c 0.9391414775 0.6511901948 0.0001336347 c 0.3784976529 0.9006437189 0.0155178630 h 0.7757116114 0.2124087035 0.0187583199 h 1.0019841321 0.5897111019 0.0306536179 c 0.9206504794 0.8052563194 0.0383402022 c 0.5506940455 0.0355539285 0.0435877664 c 0.7189543912 0.1652396522 0.0665723707 h 0.6325579656 0.4887538700 0.0548401428 h 0.8993127218 0.9906024777 0.0795285774 h 0.6023302268 0.3313104779 0.0765492092 c 0.0024835333 0.7250014214 0.0795842995 c 0.2261448565 0.1888536358 0.0848644218 h 0.3951836956 0.9024270494 0.1273945662 c 0.2063954886 0.3409479726 0.1226546038 h 0.2897961331 0.1277071056 0.1175317552 h 0.3012245132 0.6118530530 0.1146420810 c 0.0201530199 0.8672945206 0.1411095687 c 0.7292089644 0.4535095030 0.1395458206 h 0.2897961331 0.3722928944 0.6175317552 h 0.3012245132 0.8881469470 0.6146420810 c 0.0201530199 0.6327054794 0.6411095687 c 0.7292089644 0.0464904970 0.6395458206 c 0.0060365193 0.5444443856 0.6545881225 c 0.7123141847 0.1346406675 0.6496738235 c 0.2885512487 0.2379342436 0.6651907325 h 0.2546885844 0.1027339140 0.6882598666 c 0.5811894619 0.4575321578 0.6897656561 n 0.1472689194 0.7637812090 0.7040610279 c 0.7493496897 0.3299254200 0.7119607479 h 0.5258479272 0.5034909736 0.7409130304 c 0.3216422077 0.8918055436 0.7271974022 h 0.8773951383 0.9330155564 0.7239177372 c 0.1596653066 0.6772655983 0.7411508941 c 0.8685822085 0.0017679097 0.7353366102 h 0.4011250486 0.2417865996 0.7628008744 h 0.5988749514 0.7417865996 0.7371991256 c 0.1314177915 0.5017679097 0.7646633898 c 0.8403346934 0.1772655983 0.7588491059 h 0.1226048617 0.4330155564 0.7760822628 c 0.6783577923 0.3918055436 0.7728025978 h 0.4741520728 0.0034909736 0.7590869696 c 0.2506503103 0.8299254200 0.7880392521 n 0.8527310806 0.2637812090 0.7959389721 c 0.4188105381 0.9575321578 0.8102343439 h 0.7453114156 0.6027339140 0.8117401334 c 0.7114487513 0.7379342436 0.8348092675 c 0.2876858153 0.6346406675 0.8503261765 c 0.9939634807 0.0444443856 0.8454118775 c 0.2707910356 0.5464904970 0.8604541794 c 0.9798469801 0.1327054794 0.8588904313 h 0.6987754868 0.3881469470 0.8853579190 h 0.7102038669 0.8722928944 0.8824682448 c 0.7936045114 0.6590520274 0.8773453962 h 0.6048163044 0.0975729506 0.8726054338 c 0.7738551435 0.8111463642 0.9151355782 c 0.9975164667 0.2749985786 0.9204157005 h 0.3976697732 0.6686895221 0.9234507908 h 0.1006872782 0.0093975223 0.9204714226 h 0.3674420344 0.5112461300 0.9451598572 c 0.2810456088 0.8347603478 0.9334276293 c 0.4493059545 0.9644460715 0.9564122336 c 0.0793495206 0.1947436806 0.9616597978 h -0.0019841321 0.4102888981 0.9693463821 h 0.2242883886 0.7875912965 0.9812416801 c 0.6215023471 0.0993562811 0.9844821370 c 0.0608585225 0.3488098052 0.9998663653 k_points automatic 1 1 1 0 0 0 |
回復(fù)此樓» 猜你喜歡
264求調(diào)劑
已經(jīng)有7人回復(fù)
-
085600 材料與化工 298
已經(jīng)有8人回復(fù)
-
266求調(diào)劑
已經(jīng)有5人回復(fù)
-
0854總分272
已經(jīng)有6人回復(fù)
-
EST拒稿重投
已經(jīng)有3人回復(fù)
-
歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考
已經(jīng)有11人回復(fù)
-
0703化學(xué)學(xué)碩 288分
已經(jīng)有4人回復(fù)
-
接收調(diào)劑
已經(jīng)有23人回復(fù)
-
化學(xué) 0703求調(diào)劑 總分293 一志愿211
已經(jīng)有3人回復(fù)
-
304分材料專碩求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
HER12025
銅蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 17
- 散金: 33
- 帖子: 124
- 在線: 38.1小時(shí)
- 蟲(chóng)號(hào): 36722173
- 注冊(cè): 2025-12-26
- 專業(yè): 仿生材料
2樓2026-01-12 18:04:51
rlafite
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 189 (高中生)
- 金幣: 9870.1
- 紅花: 28
- 沙發(fā): 1
- 帖子: 886
- 在線: 339.2小時(shí)
- 蟲(chóng)號(hào): 4563202
- 注冊(cè): 2016-04-02
- 性別: GG
- 專業(yè): 同步輻射技術(shù)及其應(yīng)用
【答案】應(yīng)助回帖
|
你在 ph.x 中啟用了 Grimme-D3 色散修正 (vdw_corr='grimme-d3')， 但 當(dāng)前的 DFPT (ph.x) 并沒(méi)有生成或找到 D3 所需的 Hessian 文件， 因此在計(jì)算拉曼（需要二階導(dǎo)數(shù)）時(shí)直接報(bào)錯(cuò)退出。 vdw_corr = 'grimme-d3' 這是 經(jīng)驗(yàn)色散項(xiàng)，不是嚴(yán)格的 Kohn–Sham 泛函的一部分。 拉曼強(qiáng)度需要：力常數(shù)矩陣（Hessian）/ 極化率對(duì)應(yīng)變的導(dǎo)數(shù) 而當(dāng)你開(kāi)啟 D3 時(shí)：D3 對(duì)能量有貢獻(xiàn) / 對(duì)力有貢獻(xiàn) / 對(duì) Hessian 也有貢獻(xiàn) 問(wèn)題在這里👇 QE 的 ph.x 并不會(huì)自動(dòng)為 D3 構(gòu)造 Hessian |
3樓2026-01-13 02:25:29
rlafite
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 189 (高中生)
- 金幣: 9870.1
- 紅花: 28
- 沙發(fā): 1
- 帖子: 886
- 在線: 339.2小時(shí)
- 蟲(chóng)號(hào): 4563202
- 注冊(cè): 2016-04-02
- 性別: GG
- 專業(yè): 同步輻射技術(shù)及其應(yīng)用
【答案】應(yīng)助回帖
|
聲子 / 拉曼計(jì)算時(shí)關(guān)掉 D3 做法: SCF 結(jié)構(gòu)優(yōu)化階段：用 vdw_corr='grimme-d3' 最終 SCF + ph.x 階段：不用 D3 實(shí)際操作 1️⃣ 再跑一個(gè) 不含 D3 的 SCF vdw_corr = 'none' 保持：樣的結(jié)構(gòu), 同樣的 cutoff, 同樣的 prefix / outdir（或新的也可以） 2️⃣ 用這個(gè) SCF 接 ph.x 這是目前文獻(xiàn)中最常見(jiàn)、審稿人完全接受的做法 拉曼頻率對(duì) vdW 修正通常不敏感; vdW 更重要的是幾何結(jié)構(gòu) |
4樓2026-01-13 02:29:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 清風(fēng)一郎 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[論文投稿]
EST拒稿重投
 5+3
|
15102603076 2026-03-02 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
|
[考研] 085600材料與化工調(diào)劑 280分 +10 | yyqqhh 2026-03-03 | 10/500 |
|
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
|
[考研] 080500材料科學(xué)與工程 +4 | 202114020319 2026-03-03 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[論文投稿]
通訊作者寫(xiě)誰(shuí)，問(wèn)題是你意想不到的問(wèn)題
15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
|
[考研] 282求調(diào)劑 +4 | 2103240126 2026-03-02 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[考研] 0856材料與化工，270求調(diào)劑 +11 | YXCT 2026-03-01 | 13/650 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
|
[考研] 085600材料工程一志愿中科大總分312求調(diào)劑 +9 | 吃宵夜1 2026-02-28 | 11/550 |
|
|
|
[考研] 化學(xué)，材料，環(huán)境類求調(diào)劑 +7 | 考研版棒棒 2026-03-02 | 7/350 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|