| 查看: 1922 | 回复: 1 | ||||||
043114076木虫 (正式写手)
|
[交流]
PW test 学习 已有1人参与
|
|
PW test 学习 这些例子只是提供计算方法的测试，参数的选择并不严格。 atom: 氧原子单点能计算 • open shell 和 close shell 比较。注意总能量和band gap • nspin 关键字定义是否考虑自旋极化 • occupancies 定义电子排布 • pz(LDA) 与 pbe(GGA) 比较 berry: PbTiO3的极化计算 • berry phase理论参考： ∘ 文献 http://prb.aps.org/abstract/PRB/v47/i3/p1651_1 ∘ blog：http://liuxialei580.blog.163.com ... 034092009559494150/ ∘ PW的文件：bp_c_phase.f90 • 计算的体系为铁电材料 ∘ http://zh.wikipedia.org/wiki/%E9%93%81%E7%94%B5%E6%80%A7 • 计算过程 ∘ scf 得到电子密度。berry.in ∘ nscf 计算极化，对于较低对称性的体系需要计算各个方向的极化，然后加和。berry.in1 和 berry.in2。二者的差别在nbnd，结果没有差别。 ∘ 参数： ‣ lberry ‣ gdir：direction of the k-point strings in reciprocal space. ‣ nppstr:the number of k-points to be calculated along each symmetry-reduced string cluster：孤立体系的计算。 • N，NH4+，H2O。适合与带电，带磁，偶极体系？ • 关键字：assume_isolated • 文献： ∘ makov-payne：Periodic boundary conditions in ab initio calculations，Phys. Rev. B 51, 4014–4022 (1995) ∘ ESM：First-principles calculations of charged surfaces and interfaces: A plane-wave nonrepeated slab approach，Phys. Rev. B 73, 115407 (2006) ∘ martyna-tuckerman：A reciprocal space based method for treating long range interactions in ab-initio and force-field-based calculation in clusters, J.Chem.Phys. 110, 2810 (1999) ∘ dcc: Electrostatic periodic boundary conditions and real-space corrections, Phys.Rev.B 77, 115139 (2008) dipole:CO on Ni 偶极场矫正 • 文献：Dipole correction for surface supercell calculations, L. Bengtsson, PRB 59, 12301 (1999) • 关键字：tefield，dipfield,edir, eamp, emaxpos, eopreg ∘ 例子中eamp为0，不知为何这样设置。emaxpos和eopreg设置有什么标准吗？ ∘ 计算过程中，貌似amp在变化。 • 能量中多了一项：electric field correction electric：加有限的电场使用Berry's phase 方法 • Si • 关键字：lelfield，nppstr，gdir，nberrycyc，efield • k 点不能使用automatic 关键字 • 使用了 ∘ startingwfc='file', ∘ startingpot='file', ∘ 看来需要首先在没有电场的情况下运行scf（？） lda+U： • FeO • lda_plus_u，Hubbard_U ∘ 注意单位是eV • 这种方法只能使用于个别的原子。可以修改set_hubbard_l.f90 和 tabd.f90 增加自己的原子类型 lsda: • Ni • nspin, starting_magnetization, occupations, tot_magnetization ∘ 使用starting_magnetization规定的自旋在计算的过程中会变化，而tot_magnetization则固定。而且两者的标准是不一样。前者规定全自旋向上为1，全自旋向下为-1.后者是自旋向上和向下电子数的差。 ∘ 两种方法计算的结果不同，starting_magnetization 能量更低但是自旋为0.73；tot_magnetization自旋为2，与电子排布相同。 ∘ 注意：规定starting_magnetization 的例子nbnd都不够。 md: • Si • 关键字：dt，nstep, noncolin: • Fe • 关键字：noncolin, constrained_magnetization，angle， starting_magnetization spinorbit: • Pt • 关键字：lspinorb，noncolin， • 适用范围？？？ vdw: van der Waals 相互作用。DFT-D 和vdwDF方法 • 注意如果是新安装的QE，第一次运行vdw1.in会通不过，原因是没有vdW_kernel_table, 运行命令generate_vdW_kernel_table.x 生成vdW_kernel_table放到赝势的目录 |
回复此楼» 收录本帖的淘帖专辑推荐
 科研技巧 | 第一性原理计算经验 | pwscf |
» 猜你喜欢
材料与化工304求B区调剂
已经有7人回复
-
0703化学336分求调剂
已经有4人回复
-
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
-
化学工程321分求调剂
已经有12人回复
-
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有4人回复
-
材料专硕326求调剂
已经有5人回复
-
东南大学364求调剂
已经有5人回复
-
国自科面上基金字体
已经有7人回复
-
药学383 求调剂
已经有4人回复
-
286求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- PWscf 安装运行问题
已经有14人回复
- Problem - nscf - c_bands: eigenvalues not converged
已经有20人回复
- Dalton2011并行编译详记
已经有33人回复
- ESPRESSO4.3.1 install problem
已经有11人回复
- vasp计算结果处理态密度要用到的程序
已经有18人回复
- 【pwzhou个人文集】Molpro 2010.1 PL20 并行编译
已经有5人回复
- 【求助】含时密度泛函计算金团簇的激发态输出出错
已经有8人回复
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
2楼2012-07-20 20:37:58
20