关于LDA与GGA的问题。【纯新手求助】

didi_dudu

万里云 发表于 2015-6-26 10:48
输入里有个tstress，指定输出应力。

还是先找篇文章，看看别人怎么算的吧。

好的谢谢大神！我正在找这方面的东西看，还有个问题。现在quantum-espresso有没有这种功能呢，“固定晶胞参数中的某一部分，然后优化其他的部分参数”。

万里云

didi_dudu 发表于 2015-6-26 13:31
好的谢谢大神！我正在找这方面的东西看，还有个问题。现在quantum-espresso有没有这种功能呢，“固定晶胞 ...

有。设置calculation = 'vc-relax'，然后在CELL下面有个cell_dofree参数，详细看帮助文档。

didi_dudu

万里云 发表于 2015-6-26 16:23
有。设置calculation = 'vc-relax'，然后在CELL下面有个cell_dofree参数，详细看帮助文档。

太感谢了！！！！大神让我给你充点话费或者神马会员之类的表示下感谢吧

万里云

didi_dudu 发表于 2015-6-29 11:28
太感谢了！！！！大神让我给你充点话费或者神马会员之类的表示下感谢吧

这个就不必了。

didi_dudu

万里云 发表于 2015-6-29 14:21
这个就不必了。

那我就捐给论坛啦

didi_dudu

万里云 发表于 2015-6-29 14:21
这个就不必了。

大神我又来问问题了，我想算AlN的自旋极化。我尝试用PHonon里面example04中的输入来计算。
=================================
vibrations of ch4
&inputph
  tr2_ph=4.0d-17,
  prefix='ch4',
  outdir='/home/wangzf/software/espresso-5.1.2/tempdir',
  amass(1)=1.d0,
  amass(2)=12.d0,
  asr=.true.,
  zue=.true.,
  epsil=.true.,
  trans=.true.,
  fildyn='ch4.dyn.xml',
/
0.0 0.0 0.0
========================
输出来的数据是polarizability，单位是 Polarizability (a.u.)^3       ，然后我在文献中查到的是polarization，单位是C/m^2。然后我在phonon的Doc中搜索polarization关键词，发现主要出现在developer_man里面。是不是espresso程序还有一个developer版本里面具有更多的可能有点不成熟的计算功能呢？之前在文献中看到的born effective charge也主要出现在developer_man里面。 然后现在的espresso可以计算polarization么？谢谢！！！！！！
============================
PS，能不能给个建议哪本书有讲解polarization和polarizability的区别，以及effective charge和born effective charge的区别？ 我知道这种基础问题应该自己解决不过我这方面的知识实在太小白了。

万里云

didi_dudu 发表于 2015-7-1 09:17
大神我又来问问题了，我想算AlN的自旋极化。我尝试用PHonon里面example04中的输入来计算。
============ ...

developer_man是给写扩展程序的人看的，QE倒没有developer版本。

我没有算过声子，因此不知道这几个术语的区别。

算自旋极化不需要算声子啊。

didi_dudu

万里云 发表于 2015-7-1 09:37
developer_man是给写扩展程序的人看的，QE倒没有developer版本。

我没有算过声子，因此不知道这几个术 ...

额 打的时候手误了。想算的是极化。因为我想算AlN的压电系数。看文献是需要算出polarization和stress。spin-polarization和polarization是有区别的吧？ ---------------------
ps： 问个小白问题。输出中的 total energy和 totalall-electron energy 指的都是什么呢？ 谢谢！！！！！

万里云

didi_dudu 发表于 2015-7-1 10:04
额 打的时候手误了。想算的是极化。因为我想算AlN的压电系数。看文献是需要算出polarization和stress。sp ...

spin-polarization是磁相互作用，polarization是电现象，本质上是不同的事物。

total energy是总能，即Kohn-Sham粒子的动能、离子势能、Hartree相互作用能和交换关联能之和。

totall-electron energy没见过。

didi_dudu

万里云 发表于 2015-7-1 11:39
spin-polarization是磁相互作用，polarization是电现象，本质上是不同的事物。

total energy是总能， ...

大神我又来了。。我发现PWscf中可以用关键词
!#   * lberry (.TRUE. or .FALSE.)                                             #!
!#     Tells PWSCF that a Berry phase calcultion is desired.  
来算polarization。 你有没有用过这个关键词呢？ 这个polarization和通常所说的spontaneous polarization是不是一回事？

万里云

didi_dudu 发表于 2015-7-2 07:26
大神我又来了。。我发现PWscf中可以用关键词
!#   * lberry (.TRUE. or .FALSE.)                       ...

我没算过极化。

可以看安装目录下面的PW/examples/example04，这个就是用Berry相算极化的。

didi_dudu

万里云 发表于 2015-7-2 09:15
我没算过极化。

可以看安装目录下面的PW/examples/example04，这个就是用Berry相算极化的。

大神我又双叒叕来了，QE计算时候有没有要设置内存什么的？我现在算一个52个原子的晶胞，发现取k点6 6 6 cutoff 60就算不动了。。。

万里云

didi_dudu 发表于 2015-7-22 21:46
大神我又双叒叕来了，QE计算时候有没有要设置内存什么的？我现在算一个52个原子的晶胞，发现取k点6 6 6 c ...

Electrons下面有个电子对角化方式，将david改成cg，会省点内存。

didi_dudu

万里云 发表于 2015-7-22 22:17
Electrons下面有个电子对角化方式，将david改成cg，会省点内存。

大神我又双叒叕来了。。请问下QE在计算时候有没有固定某些原子然后优化其他原子的功能呢？你知不知道设置的关键词神马的。谢谢！！

万里云

didi_dudu 发表于 2015-8-13 14:10
大神我又双叒叕来了。。请问下QE在计算时候有没有固定某些原子然后优化其他原子的功能呢？你知不知道设置 ...

有，不过我现在想不起来了。

你看下input_pw.html这个文件，在CARD部分。