24小时热门版块排行榜

返回列表

043114076

木虫 (正式写手)

应助: 49 (小学生)
金币: 3364.9
散金: 52
红花: 23
帖子: 788
在线: 369.5小时
虫号: 570592
注册: 2008-06-07
性别: GG
专业: 理论和计算化学

[交流] 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放到赝势的目录

回复此楼