[Pw_forum] Atomic rho

2008-10-14 Thread Stefano de Gironcoli
In reciprocal space the sum of atomic densities is

sum_rho_at(G) = \sum_ityp S_ityp(G) rho_atomic_ityp(G)

where rho_atomic_ityp(G) is the Fourier trasform of the atomic charge  
density for atomic type ityp and S_ityp(G) is the strcuture factor for  
that atomic type

S_ityp(G) = \sum_{na such that ityp_of_na=ityp) esp(-i G tau_na) / Omega

If you want a single atom charge density you simply need to use the  
espression above with

S_ityp(G) =  esp(-i G tau_na_you_want) / Omega

this is actually a periodic charge density corresponding to atom  
na_you_want plus its
periodic images ... so the cell must be large enough that these images  
do not overlap.

Hope this helps,

  stefano


Quoting Leonardo Matheus :

> Thanks for your reply Axel.
>
>> well, for practical reasons it is more efficient to loop over
>> atom types and atoms per atom type. per atom-type parameters
>> only need to be computed once this way.
>
> I understand why the calculation is like this, but still couldn't find
> this "atoms per atom type" loop.
>
>> i suggest you first have a look into textbook material
>> explaining the implementation of electronic structure
>> calculation with plane wave pseudopotentials (e.g. the
>> books by jorge kohnoff or richard martin). and then
>> read the marx / hutter review which also gives an overview
>> on how the implementation can be done (on the example of
>> the CPMD code, which is not that different as far as the
>> basic principles are concerned). the CPMD code actually
>> does have an implementation of hirshfeld charges determination,
>> but only for norm-conserving pseudopotentials, IIRC, so you
>> can look up how the implementation principles.
>
> I've already read the book by Kohanoff, and now I'm waiting for the
> licence to download the CPMD code, to see how Hirshfeld charges are
> implemented there, thanks for the hint. ABINIT code also has Hirshfeld
> charges implemented, I'm looking there too.
>
>> for implementing it into QE you'll have to figure out
>> how to handle ultra-soft psps contribution to the atomic
>> density on top of everything else.
>
> I'm also using the SaX code, which accepts just NC PPs, so I haven't
> thought deeply on this, but I would first try to adapt the "addusdens"
> subroutine to the individual atomic densities (at least, this is what
> I understood the subroutine does).
>
> Actually, I've based most of my code on the "voronoy" program, so what
> I really need now is just to obtain the individual atomic rho, the
> rest (I think) is mostly done.
>
>> good luck,
>>   axel.
>
> Thanks
>
> Leonardo Matheus
> Doctorate student - Physics Institute, University of S?o Paulo, Brazil
> Nanomol group - http://nanomol.if.usp.br
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>




   SISSA Webmail https://webmail.sissa.it/
   Powered by Horde http://www.horde.org/




[Pw_forum] Getting DOS at Fermi Energy

2008-10-14 Thread Paolo Giannozzi
Paul M. Grant wrote:

> BTW, have you had a chance to look
> at the k-points file for aluminum I sent?

it looks like a uniform unshifted 12x12x12 grid in a simple
cubic lattice. You should get it (or an equivalent one)
directly from pw.x using options nosym=.true., noinv=.true.
(the latter is undocumented: removes q=>-q symmetry) and
K_POINT {automatic}
12 12 12 0 0 0
or else, from kpoints.x, answering 'T' to the question
   write all k? [f] >>
I am not aware of any limitation to the simple cubic grid,
but this should be asked to the authors of epsilon.x

Paolo
-- 
Paolo Giannozzi, Democritos and University of Udine, Italy


[Pw_forum] Atomic rho

2008-10-14 Thread Leonardo Matheus
Thanks for your reply Axel.

> well, for practical reasons it is more efficient to loop over
> atom types and atoms per atom type. per atom-type parameters
> only need to be computed once this way.

I understand why the calculation is like this, but still couldn't find
this "atoms per atom type" loop.

> i suggest you first have a look into textbook material
> explaining the implementation of electronic structure
> calculation with plane wave pseudopotentials (e.g. the
> books by jorge kohnoff or richard martin). and then
> read the marx / hutter review which also gives an overview
> on how the implementation can be done (on the example of
> the CPMD code, which is not that different as far as the
> basic principles are concerned). the CPMD code actually
> does have an implementation of hirshfeld charges determination,
> but only for norm-conserving pseudopotentials, IIRC, so you
> can look up how the implementation principles.

I've already read the book by Kohanoff, and now I'm waiting for the
licence to download the CPMD code, to see how Hirshfeld charges are
implemented there, thanks for the hint. ABINIT code also has Hirshfeld
charges implemented, I'm looking there too.

> for implementing it into QE you'll have to figure out
> how to handle ultra-soft psps contribution to the atomic
> density on top of everything else.

I'm also using the SaX code, which accepts just NC PPs, so I haven't
thought deeply on this, but I would first try to adapt the "addusdens"
subroutine to the individual atomic densities (at least, this is what
I understood the subroutine does).

Actually, I've based most of my code on the "voronoy" program, so what
I really need now is just to obtain the individual atomic rho, the
rest (I think) is mostly done.

> good luck,
>   axel.

Thanks

Leonardo Matheus
Doctorate student - Physics Institute, University of S?o Paulo, Brazil
Nanomol group - http://nanomol.if.usp.br


[Pw_forum] scf-convergence

2008-10-14 Thread Eyvaz Isaev
Dear Asa,

A few comments on your input file.

1. UsXCrysDen have shown that, presumably, what you study is not a monowire, it 
looks quite strange. Be sure that you have correct structure.

2. You use UltraSoft pseudopotentials, so ecutwfc should be around 10-12 times 
of ecutrho.

3. Again, you mean you have body centered orthorhombic structure, but looking 
at the structure (by XCrysDen) there is no cell center.

Bests,
Eyvaz.


---
Prof. Eyvaz Isaev, 
Theoretical Physics Department, Moscow State Institute of Steel & Alloys, 
Russia, 
Department of Physics, Chemistry, and Biology (IFM), Linkoping University, 
Sweden 
Condensed Matter Theory Group, Uppsala University, Sweden 
Eyvaz.Isaev at fysik.uu.se, isaev at ifm.liu.se, eyvaz_isaev at yahoo.com


--- On Mon, 10/13/08, asa aravindh  wrote:

> From: asa aravindh 
> Subject: [Pw_forum] scf-convergence
> To: pw_forum at pwscf.org
> Date: Monday, October 13, 2008, 3:28 PM
> Dear all..
> I am trying to do scf calculation of an fe/co nanowire
> along the bcc (110) direction.
> However the calculation is not getting converged..after
> even 100 iterations.. the scf accuarcy is not reached and
> suddenly the program stops ...
> can anybody pls tell me if anything is wrong with my input
> file...
> the input file is attached
> any reply would be greatly appreciated
> regards,
> asa
> 
> 
>   Add more friends to your messenger and enjoy! Go to
> http://messenger.yahoo.com/invite/___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum


  


[Pw_forum] possible bug in scatter_forw.f90

2008-10-14 Thread Manoj Srivastava
Dear Alexader, 

 Thank you very much for your reply. I was out of town and could not see
the message till yesterday. I have some questions on it, it would be great
if you could answer them.
On Fri, 10 Oct 2008, Alexander wrote:

> Dear Manoj
> 
> On Friday 10 October 2008 05:20, Manoj Srivastava wrote:
> > Any follow ups on the message below? Alexander?
> >
> > Regards,
> > Manoj
> >
> > On Wed, 8 Oct 2008, Manoj Srivastava wrote:
> > > Dear All,
> > >  Is any body familiar with the scatter_forw.f90 subroutine of PWCOND? I
> > > think there is a bug in this subroutine at the place where you calculate
> > > intw2, which is z integration of nonlocal wavefunction with beta
> > > function. I have looked up the necessary formulae in the Choi & Ihm's
> > > paper (PRB 59, 2267, Jan 1999). In the paper, nonlocal wavefunction has 3
> > > terms, one of which contains beta function, say W. The other two parts
> > > dont have any W in them, rather they are just plane wave solution.
> > > (Please have a look at equation 24 and 26 of the paper ). So, when we are
> > > doing z integration of nonlocal wavefunction with beta function W, we
> > > should have three terms, one of which should contain two W, but rest
> > > should just have one W. On the other hand in the code all three terms
> > > contain two beta functions!! (please have a look at line 228 of
> > > scatter_forw.f90). I am wondering if my understanding is right?
> 
> Your understanding is wrong. On the line 228 you add to the integral ONLY
> contribution with two beta functions (from the 1st term in the nonlocal 
> function, see Eq. 24) but on the line 393 the contribution from the 2nd term 
> is added too. There is no 3rd term since every step you rotate your local 
> solutions \psi_n in such a way that b_{\lambda \alpha lm} vanish (see the 
> paragraph after Eq. 37).
  I believe that the equations 40, 41 and 42 implemented in the code are
in the momentum space. In the code, intw2 has 3 terms. We can tell this by
looking at the subroutine integrals.f90, where int2d is defined. Now, by
looking at the choi& Ihm 's paper, the very first term of psi{alpha,l,m}
has f{lambda,alpha,l,m}(equation26).And, if we substitute this in equation
40, we should just get one term. I am just confuesd on the fact that how
come this one term in paper gets transformed into 3 terms of the code!
Would you mind explaining? Also, the integrals that are defined in the
code are only in the slabs (nz1), so how are you taking the integration in
the region (0,d), which we need according to equations 40, 41 and 42 ?

> 
>  > > I have one more question in the later part of the same subroutine. What
> > > does the lapack subroutine
> > > ZGESV(2*n2d,2*n2d+norb*npol,amat,2*n2d,ipiv,xmat,2*n2d,info)
> > > do?
> > > I tried to look it up and i got the idea that it is trying to solve
> > > amat*x=xmat, with amat and xmat known and x unknown, and at the end of
> > > calculation it stores x in xmat.
> > > so, basically it does--  x=[(amat)^(-1)]*xmat. Am i right? So, it changes
> > > the structure of xmat from what is defined in 'constructs matrices'
> > > part. Is it correct?
> 
> Here, you are right.
> 
> > > Also, afterwards in this code where  it 'rotates integrals' is not very
> > > clear to me.
> > > Could somebody please tell me in little detail, what is going on
> > > here? 
> As I explained before rotation means that every step (at each slab)
> you make a linear combination of local solutions \psi_n with the 
> transformation matrix  h_{nn'} (see the last paragraph on the p. 2270).
> so that the new solutions have the property: 
> b_{\lambda n} = \delta_{\lambda n} and the new nonlocal solutions have
> b_{\lambda \alpha lm}=0. Doing this transformation of local and nonlocal
> solutions you should also perform the corresponding transformations
> of the integrals.
> 
> > > Also, is this subroutine written just on the basis of Choi and Ihm 
> > > paper, or are there more reference to it? If yes, would someone mind
> > > mentioning them?
> Unfortunately, the paper of Choi and Ihm is very detail, and in the 
> PWCOND code we followed it very closely, so no more references
> are needed. In our papers we just extended those ideas to ultrasoft
> pseudopotentials, magnetism, spin-orbit coupling and so on.
>   
> Hope this helped you,
> regards, Alexander 
> 
> > >
> > > Regards,
> > > Manoj Srivastava
> > > Physics Graduate Student
> > > University of Florida,
> > > Gainesville,FL, USA
> >
> > ___
> > Pw_forum mailing list
> > Pw_forum at pwscf.org
> > http://www.democritos.it/mailman/listinfo/pw_forum
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
> 

Once again thank you very much for your help. I appreciate it. 

Regards, 
Manoj Srivastava



[Pw_forum] question about output of SCF

2008-10-14 Thread yumin qian
I have done a scf calculation of ground state Pbnm CaCrO3 crystal , the
total band nbnd=101
there are 152 irreducible k-points  ,but the output file in X.save directory
only contians 76 folders : K1-K00076, and the charge-desity.xml only
contains  to  , what is the relation between charge desity grid,
and the number of K000XX folders.

   I also do not know why we define Main grid and wavefunction grid
respectively ? Since the grid number are determined by the same ecutwfc ,
could someone explain more clearly about the all the grids system (include
FFT grid, and charge density grid in real or G-space) in DFT ?

  Many thanks

-- 
Sincerely Y. M. Qian
Lab.of Condensed Matter Theory and Materials Computation
Institute of Physics
Chinese Academy of Sciences
Tel:  + 8610 8264 9147
E-Mail:yuminqian at gmail.com <E-Mail%3Ayuminqian at gmail.com>
P.O.Box 603   Beijing 100190
China
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20081014/5ba1febb/attachment.htm
 


[Pw_forum] Bug in axsf generation?

2008-10-14 Thread Tone Kokalj
On Fri, 2008-10-10 at 08:23 -0600, Jan Beck wrote:
> I have an axsf file that I have generated with quantum espresso's cppp. I
> have checked with the specification of the axsf file format, and the
> output from quantum espresso is consistent with that.

> WARNING: Atom 2 and atom 6 overlap !!!
>   Atom 6 deleted !!!
> WARNING: Atom 2 and atom 10 overlap !!!
>   Atom 10 deleted !!!

I have looked at your file. It is nothing wrong with its syntax, but the
structure is a complete mass: there are many overlapping atoms (look
into the file). Therefore you receive the atom-deleted warnings!

Either your original cp.x file already contains massy structure, or the
cppp utility has a bug (i never used the latter).

Regards, Tone
-- 
Tone Kokalj 
J. Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia (tel:
+386-1-477-3523 // fax:+386-1-477-3822)



[Pw_forum] k-point grids...and all that

2008-10-14 Thread Deyu Lu
>Finally, on the point of mesh visualization.  I use XCrySDen a lot,
>especially for BZ and charge-potential visualizations.  But I don't know how
>to configure or generate a .xsf file for a 3D mesh of discrete points.

W.r.t. visualization of k-points, if the coordinates of k-points and the 
boundary
of the BZ are known, they can be visualized in VMD with the .xyz file format.
Below is an example of the 10 fcc special k-points (blue spheres) in the 
irreducible wedge (white spheres).

http://picasaweb.google.com/deyu.lu/Misc#5257056499972518002

Best
Deyu Lu


***
Deyu Lu (Ph.D)
190 Chemistry Building
University of California, Davis 
One Shields Avenue
Davis, CA 95616
Office phone: (530) 754-9663
Group Webpage: http://angstrom.ucdavis.edu/
***



  


[Pw_forum] Getting DOS at Fermi Energy

2008-10-14 Thread Paolo Giannozzi
Paul M. Grant wrote:

> Look in the PH directory and not PW or PP?  How obvious!

it is the only part of the entire distribution where the DOS
at E_F is presently used. There is a reason for everything in
Q-E.  Sometimes it is not a good reason, or it is an obsolete
reason, or nobody knows any longer what the reason was, but
there is always a reason...

> However, thanks for pointing to the dos_ef function in elphon.f90 
> (which also needs the external routine w0gauss.f90)

it is in PW/w0gauss.f90 and is used in many other places

> I think my main challenge will be to figure out how to extract
> the arrays et and wk from the nefarious *.save folders.
>  Maybe you can point me to some routines where this task
> is handled...and for spin polarized states as well?

if you start from scratch, it will take a lot of time. I think
the minimum energy path here is one of the following:
- add a call to dos_ef in the main pw.x code, whenever the Fermi
   energy is calculated (PW/weights.f90) ;
- modify any of the codes in PP/ that read the nefarious .save
   folders. They typically rebuild the structure of variables as
   used in pw.x ;
- build upon the "qexml" library, as suggested by Axel, but
   beware: I am not sure it is updated

Paolo
-- 
Paolo Giannozzi, Democritos and University of Udine, Italy


[Pw_forum] question about output of SCF

2008-10-14 Thread Paolo Giannozzi
yumin qian wrote:

> I have done a scf calculation of ground state Pbnm CaCrO3 crystal , the 
> total band nbnd=101 
> there are 152 irreducible k-points  ,but the output file in X.save 
> directory only contians 76 folders : K1-K00076

so there are 76 irreducible k-points, and since you are using
LSDA, the code doubles the number of k-points during the
calculation, so that you have 76 k-points for spin-up states
and 76 for spin-down. At the end, when results are saved,
spin-up and spin-down states are collected into a single
folder.

> I also do not know why we define Main grid and wavefunction grid 
> respectively ? Since the grid number are determined by the same 
> ecutwfc

not for USPP or PAW

> , could someone explain more clearly about the all the grids system 
> (include FFT grid, and charge density grid in real or G-space) in DFT ?

the best explanations I managed to produce are reported in
the 3rd paper listed here:
http://www.quantum-espresso.org/wiki/index.php/Ref:_First-principles_Molecular_Dynamics
You can get a copy here:
http://www.fisica.uniud.it/~giannozz/Papers/JCP_57.pdf

Paolo
-- 
Paolo Giannozzi, Democritos and University of Udine, Italy


[Pw_forum] Getting DOS at Fermi Energy

2008-10-14 Thread Paul M. Grant
..."there is always a reason"...sounds like some sort of creationist mantra
to me!

Thanks for the insight into pw.x...maybe in a few light-years I'll be able
to contribute to the cause.

I'll have a look at method you suggest, but last night I think I've found a
straightforward (relatively) way to hack epsilon.f90 with a dos_ef like
option.  Let you know how it turns out.  BTW, have you had a chance to look
at the k-points file for aluminum I sent?

Paul M. Grant, PhD
Principal, W2AGZ Technologies
Visiting Scholar, Applied Physics, Stanford (2005-2008)
EPRI Science Fellow (Retired)
IBM Research Staff Member Emeritus
w2agz at pacbell.net
http://www.w2agz.com
?
?


-Original Message-
From: pw_forum-bounces at pwscf.org [mailto:pw_forum-boun...@pwscf.org] On
Behalf Of Paolo Giannozzi
Sent: Tuesday, October 14, 2008 12:25 AM
To: PWSCF Forum
Subject: Re: [Pw_forum] Getting DOS at Fermi Energy

Paul M. Grant wrote:

> Look in the PH directory and not PW or PP?  How obvious!

it is the only part of the entire distribution where the DOS
at E_F is presently used. There is a reason for everything in
Q-E.  Sometimes it is not a good reason, or it is an obsolete
reason, or nobody knows any longer what the reason was, but
there is always a reason...

> However, thanks for pointing to the dos_ef function in elphon.f90 
> (which also needs the external routine w0gauss.f90)

it is in PW/w0gauss.f90 and is used in many other places

> I think my main challenge will be to figure out how to extract
> the arrays et and wk from the nefarious *.save folders.
>  Maybe you can point me to some routines where this task
> is handled...and for spin polarized states as well?

if you start from scratch, it will take a lot of time. I think
the minimum energy path here is one of the following:
- add a call to dos_ef in the main pw.x code, whenever the Fermi
   energy is calculated (PW/weights.f90) ;
- modify any of the codes in PP/ that read the nefarious .save
   folders. They typically rebuild the structure of variables as
   used in pw.x ;
- build upon the "qexml" library, as suggested by Axel, but
   beware: I am not sure it is updated

Paolo
-- 
Paolo Giannozzi, Democritos and University of Udine, Italy
___
Pw_forum mailing list
Pw_forum at pwscf.org
http://www.democritos.it/mailman/listinfo/pw_forum



[Pw_forum] Bug in axsf generation?

2008-10-14 Thread Paolo Giannozzi
Jan Beck wrote:

> this is my third attempt at posting this reply, I cannot paste the
> .axsf file or attach it, because its too large. The message bounces
> back to me with a "moderator approval pending" message. I think
> approval is still pending

apparently you didn't notice the "Important notice" here: 
http://www.quantum-espresso.org/wiki/index.php/Introduction#Contacts
Executive summary:
---
Messages "waiting for moderator approval" are automatically
deleted with no further processing (sorry, too much spam)
---
> so I apologize for making you download the
> file from a 3rd party site:
> 
> http://rapidshare.com/files/153722255/data.axsf.html

you shouldn't apologize: this is the correct way to send
large files, in my opinion.

Paolo
-- 
Paolo Giannozzi, Democritos and University of Udine, Italy