Skip to Content.
Sympa Menu

forum - Re: Re: [abinit-forum] Response  to electric fields,  strain at other "qpt"&nbsp ;other than Gamma?

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: Re: [abinit-forum] Response  to electric fields,  strain at other "qpt"&nbsp ;other than Gamma?


Chronological Thread 
  • From: "duan_semi" <abinit_duan@126.com>
  • To: "xavier gonze" <gonze@pcpm.ucl.ac.be>
  • Cc: forum@abinit.org
  • Subject: Re: Re:&nbsp;[abinit-forum]&nbsp;Response &nbsp;to&nbsp;electric&nbsp;fields,&nbsp; strain&nbsp;at&nbsp;other&nbsp;"qpt"&nbsp ;other&nbsp;than&nbsp;Gamma?
  • Date: Wed, 25 Oct 2006 10:02:31 +0800 (CST)

Dear Xavier Gonze,
          Thank you very much for you warmhearded help. I have calculated the strain and homogeneous electric perturbation successfully when "qpt" is set to (0, 0, 0). My question is whether I can do the same calculaton when "qpt" is set to (0.5, 0.0, 0.0)?
Best regards,
WD Duan
 
 
 


-----原始邮件-----
发件人:"Xavier Gonze"
发送时间:2006-10-25 00:08:00
收件人:forum@abinit.org
抄送:(无)
主题:Re: [abinit-forum] Response to electric fields, strain at other "qpt" other than Gamma?


Dear WD Duan,

The strain and homogeneous electric field perturbations have  
(naturally) zero wavevector,
as mentioned in the ERROR message that you get.
So, please, use qpt 0.0 0.0 0.0 to compute electric field and strain  
response functions.
There are lots of examples in the corresponding tutorials. E.g. for  
ABINITv5 see the entry tutorial page
http://www.abinit.org/Infos_v5.2/tutorial/welcome.html .

Xavier Gonze

On 20 Oct 2006, at 17:59, abinit_duan@126.com wrote:

> Dear Horacio Wagner Leite Alves,
>        Firstly, thanks a lot for your help. I have calculated the  
> phonon spectrum successfully when "qpt" is set to (0.5  0.0  0.0),  
> but at the same condition when I compute the response function to  
> "electric field" and "strain",some thing is wrong as following:
>   "getcut: wavevector=  0.0000  0.0000  0.0000  ngfft=  18  18  18
>          ecut(hartree)=      6.000   => boxcut(ratio)=   2.07835
>  ==>  initialize data related to q vector <==
>  respfn : ERROR -
>   The treatment of non-zero wavevector q is restricted to phonons.   
> However, the input normalized qpt is    5.000000E-01    0.000000E 
> +00    0.000000E+00,
>   while rfelfd= 2, and rfstrs= 0.
>   Action : change qpt, or rfelfd or rfstrs in the input file."
> I do not know how to deal with it. So I wonder if we could  
> calculate the response functions to "electric field" and "strain",  
> if it can be calculated, how? Would you mind sending me some  
> example of this kind?
> Thanks in advance!
> WD Duan
>








金 秋 美 女 靓 妆 大 比 拼 ( 图 ) 绝 对 经 典 !
周 末 和 女 友 去 buy 的 好 东 东 ( 组 图 )


Archive powered by MHonArc 2.6.16.

Top of Page