| 6 | 1/1 | 返回列表 |
| 查看: 1266 | 回復(fù): 5 | |||
enola捐助貴賓 (正式寫手)
|
[交流]
【分享】pwscf模擬水分子 已有5人參與
|
|
看見有人在MS版面上提到水分子的模擬,今天有空自己做了下,可是出來結(jié)果不對(輸入時,特意把H-O-H鍵角是105 deg,但pw優(yōu)化后成了97.5deg )不知道是不是水中的H得考慮自旋阿,還是其他有忽略的地方,輸入文件如下,請大家批評阿: &CONTROL calculation = 'relax', prefix = 'H2O', pseudo_dir = '/home/enola/pwscf/espresso-3.2.3/pseudo', outdir = '/home/enola/tmp', tprnfor=.ture forc_conv_thr=1.0D-4 nstep=100 / &SYSTEM ibrav = 0, nat = 3, ntyp = 2, ecutwfc = 35.D0, ecutrho = 150.D0, / &ELECTRONS conv_thr = 1.D-7, mixing_beta = 0.7D0, / &IONS ion_dynamics='bfgs', pot_extrapolation = 'second_order', wfc_extrapolation = 'second_order', / CELL_PARAMETERS cubic 20.0 0.0 0.0 0.0 20.0 0.0 0.0 0.0 20.0 ATOMIC_SPECIES H 1.00d0 H.fpmd.UPF O 16.0d0 O.BLYP.UPF ATOMIC_POSITIONS (bohr) O 0.0099 0.0099 0.0000 H 1.4325 -0.2243 -0.0001 H -0.2243 1.8325 0.0002 K_POINTS {Gamma} |

鐵桿木蟲 (著名寫手)
金蟲 (著名寫手)
|
贗勢用的是模守恒吧,我想ecutrho = 150.D0應(yīng)該差不多了。但是ecutwfc太小,有H在這里,你必須用大一些。建議至少60Ry. 此外,我想很可能是初始構(gòu)型的問題。我對水分子進行過無數(shù)次優(yōu)化,都是對的。 你可以考慮用別的贗勢試一下。 版本太低了啊。呵呵 我總覺得整個輸入文件缺點什么。 [ Last edited by goldenfisher on 2010-6-5 at 19:49 ] |
捐助貴賓 (正式寫手)

|
我用pwscf做出了單體水 也做出了二聚體的水 結(jié)構(gòu)也很好包括氫鍵,1.97Angstrom H-O為0.976Angstrom 104.5degree. 氫鍵能0.229eV是很好的還計算了diamond-phase Ice 也很好和實驗符合的很好,給你一個建議,就是換一下勢函數(shù)。我的問題是怎么加VdW修正?除了pwscf中的例子 要像cpmd中的VdW correction! thank you regards! |

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +6 | 嘉年新程 2026-03-15 | 6/300 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|