| 6 | 1/1 | 返回列表 |
| 查看: 1265 | 回復(fù): 5 | |||
enola捐助貴賓 (正式寫手)
|
[交流]
【分享】pwscf模擬水分子 已有5人參與
|
|
看見有人在MS版面上提到水分子的模擬,今天有空自己做了下,可是出來結(jié)果不對(duì)(輸入時(shí),特意把H-O-H鍵角是105 deg,但pw優(yōu)化后成了97.5deg )不知道是不是水中的H得考慮自旋阿,還是其他有忽略的地方,輸入文件如下,請(qǐng)大家批評(píng)阿: &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)型的問題。我對(duì)水分子進(jìn)行過無數(shù)次優(yōu)化,都是對(duì)的。 你可以考慮用別的贗勢試一下。 版本太低了啊。呵呵 我總覺得整個(gè)輸入文件缺點(diǎn)什么。 [ 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是很好的還計(jì)算了diamond-phase Ice 也很好和實(shí)驗(yàn)符合的很好,給你一個(gè)建議,就是換一下勢函數(shù)。我的問題是怎么加VdW修正?除了pwscf中的例子 要像cpmd中的VdW correction! thank you regards! |

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +3 | NIFFFfff 2026-03-20 | 3/150 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 一志愿中國石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|