[Pw_forum] cannot allocate memory error

2011-09-03 Thread Stefano Baroni

On Sep 3, 2011, at 5:36 PM, Elie Moujaes wrote:

> Thanks a lot for the reply ..I was trying to gain some time as I am running 
> other programs on the server at Uni. I guess I will have to do things step by 
> step

that's a very good idea ;-)
take care - SB

> 
> Thanks again
> 
> Elie

---
Stefano Baroni - SISSA  &  DEMOCRITOS National Simulation Center - Trieste
http://stefano.baroni.me [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni 
(skype)

La morale est une logique de l'action comme la logique est une morale de la 
pens?e - Jean Piaget

-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/4bb201c8/attachment.htm
 


[Pw_forum] cannot allocate memory error

2011-09-03 Thread Elie Moujaes

Thanks a lot for the reply ..I was trying to gain some time as I am running 
other programs on the server at Uni. I guess I will have to do things step by 
step
Thanks again
Elie

> From: ttduyle at gmail.com
> Date: Sat, 3 Sep 2011 11:30:00 -0400
> To: pw_forum at pwscf.org
> Subject: Re: [Pw_forum] cannot allocate memory error
> 
> On Sat, Sep 3, 2011 at 11:23 AM, Elie Moujaes
>  wrote:
> > Dear all,
> > I am having a problem when I am running a scf calculation for a 72-atom
> > supercell of Carbon grain boundaries. The calculations start and then after
> > a few minutes stop with the following error:
> > Operating system error: Cannot allocate memory
> > operating system error: cannot allocate memory
> > out of memory
> > Can anyone tell me what is going on? I am not sure if it is really a memory
> > problem or something else?
> Yes. It said that. Your system is quite large for running on a
> personal desktop. Check output to see how much memory you need for
> this system.
> --
> Duy Le
> PhD Candidate
> Department of Physics
> University of Central Florida.
> 
> "Men don't need hand to do things"
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
  
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/f4e56844/attachment-0001.htm
 


[Pw_forum] cannot allocate memory error

2011-09-03 Thread Elie Moujaes

Dear all,
I am having a problem when I am running a scf calculation for a 72-atom 
supercell of Carbon grain boundaries. The calculations start and then after a 
few minutes stop with the following error:
 Traceback: not available, compile with -ftrace=frame or -ftrace=fullOperating 
system error: Cannot allocate memory
Memory allocation failed
on my desktop at work.  I have the espresso-4.2 there. On my laptop at home 
(where I have downloaded the latest espresso-4.3.2 version). Now  I am getting 
the error :
operating system error: cannot allocate memoryout of memory
Can anyone tell me what is going on? I am not sure if it is really a memory 
problem or something else?
Thanks
Elie MoujaesUniversity of NottinghamNG7 2RDUK   
  
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/f915da4d/attachment.htm
 


[Pw_forum] Method of computing static dielectric constant of materials

2011-09-03 Thread Rajan Pandey
I forgot to add the affiliation in my previous post. Correcting myself, and
regret any inconvenience.

-- Forwarded message --
From: Rajan Pandey 
Date: Sat, Sep 3, 2011 at 4:13 PM
Subject: Method of computing static dielectric constant of materials
To: PWSCF Forum 


Dear Quantum Espresso community,

I am trying to compute the static dielectric constant of materials using the
methodology
as implemented in CP code (part of the Quantum Espresso). The method is
explained
in the example30 of Quantum Espresso distribution. The system discussed is
MgO,
Ref: P.Umari and A.Pasquarello, Physical Review Letters, 89, p.157602
(2002).

The example30 mentioned above has following note:

"NOTE: the electronic dipole is defined modulo a factor (2*L=31.824i a.u.,
during the MD simulation the term "ln det S" changes the Riemann
plane, this must be taken into account when addressing the
electronic dipole."

The above statement is given in the context of a cubic supercell with length
L containing 64 atoms of MgO.

I am studying a well known system, SiO2 (alpha quartz) for sanity check. I
am simulating SiO2 in a 1x1x2
trigonal supercell (18 atoms) as well as in an orthorhombic (in alpha quartz
case it will be tetragonal because a = b )
supercell. The ambiguity is that when I use the formula along with the
"NOTE" of example30 in Quantum Espresso,
mentioned above, I get wrong results for static dielectric constant.
However, when I do not use the
"NOTE" mentioned above, and by using the formula described in example30, I
get the simulated value close
to the experimental static dielectric constant of SiO2. This means that the
method works for non-cubic unit cells too.
I am not able to understand the role of the point mentioned in the NOTE
(above) in general, when the formula described
in example30 (reiterated below), seems to work.

The difference d_Eps  between static and high-frequency dielectric constant
is given by:

d_Eps=4*pi*(D2_el + D2_ion - D1_el - d1_ion)/(0.001 a.u. * Omega)

D1(2)_el = electronic contribution of the dipole at the beginning (end) of
the relaxation.
D1(2)_ion = ionic contribution of the dipole at the beginning (end) of the
relaxation.
Omega = supercell volume

I shall appreciate any comments from community members.

Thanks, and regards,

Rajan

Rajan K. Pandey, Ph.D.

Advisory Research Engineer,
Semiconductor Research & Development Center
India Systems & Technology Engineering Lab
IBM India Pvt. Ltd.
MD3 1F B354
Manyata Embassy  Business Park
Nagawara, Outer Ring Road
Bangalore - 560045, India
Phone: +91-80-28061262
Mobile: +91-9901850981
Email: rajapand at in.ibm.com
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/dc20ec01/attachment.htm
 


[Pw_forum] Method of computing static dielectric constant of materials

2011-09-03 Thread Rajan Pandey
Dear Quantum Espresso community,

I am trying to compute the static dielectric constant of materials using the
methodology
as implemented in CP code (part of the Quantum Espresso). The method is
explained
in the example30 of Quantum Espresso distribution. The system discussed is
MgO,
Ref: P.Umari and A.Pasquarello, Physical Review Letters, 89, p.157602
(2002).

The example30 mentioned above has following note:

"NOTE: the electronic dipole is defined modulo a factor (2*L=31.824i a.u.,
during the MD simulation the term "ln det S" changes the Riemann
plane, this must be taken into account when addressing the
electronic dipole."

The above statement is given in the context of a cubic supercell with length
L containing 64 atoms of MgO.

I am studying a well known system, SiO2 (alpha quartz) for sanity check. I
am simulating SiO2 in a 1x1x2
trigonal supercell (18 atoms) as well as in an orthorhombic (in alpha quartz
case it will be tetragonal because a = b )
supercell. The ambiguity is that when I use the formula along with the
"NOTE" of example30 in Quantum Espresso,
mentioned above, I get wrong results for static dielectric constant.
However, when I do not use the
"NOTE" mentioned above, and by using the formula described in example30, I
get the simulated value close
to the experimental static dielectric constant of SiO2. This means that the
method works for non-cubic unit cells too.
I am not able to understand the role of the point mentioned in the NOTE
(above) in general, when the formula described
in example30 (reiterated below), seems to work.

The difference d_Eps  between static and high-frequency dielectric constant
is given by:

d_Eps=4*pi*(D2_el + D2_ion - D1_el - d1_ion)/(0.001 a.u. * Omega)

D1(2)_el = electronic contribution of the dipole at the beginning (end) of
the relaxation.
D1(2)_ion = ionic contribution of the dipole at the beginning (end) of the
relaxation.
Omega = supercell volume

I shall appreciate any comments from community members.

Thanks, and regards,

Rajan
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/eee0c9e2/attachment.htm
 


[Pw_forum] Wannier90, question about units in seedname_hr.dat

2011-09-03 Thread Julen Ibanez Azpiroz
Dear Prof. Giannozzi

Thank you very much for the answer, I thought QE and W90 shared the mailing
list, I will try in W90's, thanks again,


-- 

?
 Julen Iba?ez Azpiroz
 Materia Kondentsatuaren Fisika Saila
 Zientzia eta Teknologia Fakultatea
 Euskal Herriko Unibertsitatea
 644 Posta Kutxatila, 48080 Bilbo, Spain

 Telefonoa: +34 946015326
 Mail: jibanez008 at ikasle.ehu.es

-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/55d6277d/attachment.htm
 


[Pw_forum] cannot allocate memory error

2011-09-03 Thread Duy Le
On Sat, Sep 3, 2011 at 11:23 AM, Elie Moujaes
 wrote:
> Dear all,
> I am having a problem when I am running a scf calculation for a 72-atom
> supercell of Carbon grain boundaries. The calculations start and then after
> a few minutes stop with the following error:
> Operating system error: Cannot allocate memory
> operating system error: cannot allocate memory
> out of memory
> Can anyone tell me what is going on? I am not sure if it is really a memory
> problem or something else?
Yes. It said that. Your system is quite large for running on a
personal desktop. Check output to see how much memory you need for
this system.
--
Duy Le
PhD Candidate
Department of Physics
University of Central Florida.

"Men don't need hand to do things"


[Pw_forum] Problem with Bi2Se3 SCF converge

2011-09-03 Thread mohnish pandey
Dear Wu !

Can you please paste your complete input file so that we can take a look and
may be able to tell where is the problem.

On Sat, Sep 3, 2011 at 6:36 AM, WF  wrote:

> Dear Mohnish:
>Thanks for our advices. I tried beta = 0.1 and 0.2, and those
> message is postponed to 100+ iteration, but the convergence still cannot be
> achieved.
>
>
> On Fri, Sep 2, 2011 at 11:45 AM, mohnish pandey 
> wrote:
> >  Try reducing "beta"
>
> >On Wed, Aug 31, 2011 at 8:07 AM, WF  wrote:
>
> > Hello everyone,
> >
> >  I am doing calculation of Bi2Se3 with vdW-DF functional, but
> > it really hard to converge. After several iterations it keeps to
> > prompt things
> > like:
> >
> > ** **
> >
> >  iteration # 58 ecut=60.00 Ry beta=0.70
> >
> >  CG style diagonalization
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  2 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  2 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  1 eigenvalues not converged
> >
> >  c_bands:  2 eigenvalues not converged
> >
> >  c_bands:  2 eigenvalues not converged
> >
> >  c_bands:  2 eigenvalues not converged
> >
> >  c_bands:  2 eigenvalues not converged
> >
> >  
> >
> > and will not converge in 1000 iteration.
> >
> > ** **
> >
> >  But the calculation of LDA and PBE can converge, however. All
> > pseudo-potential is generated by ld1.x with same parameters except
> > functional.
> >
> > ** **
> >
> >  I have tried to change some parameter in SCF but it doesn?t
> > seem to be helpful. Is there any suggestion for this problem?
> >
> > ** **
> >
> >  Thanks.
> >
> > ** **
> >
> > F,Wu
> >
> > ** **
> >
> > ___
> > Pw_forum mailing list
> > Pw_forum at pwscf.org
> > http://www.democritos.it/mailman/listinfo/pw_forum
> >
> >
>
> -
> Regards
>
> ---
> F, Wu
> College of Chemistry and Molecular Engineering
> Peking University
> --
>
>
>
>
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>



-- 
Regards,
MOHNISH,
-
Mohnish Pandey
BTech-Mtech, IIT Kanpur
Senior Project Associate,
Department of Chemical Engineering,
IIT KANPUR, UP, INDIA
-
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/6894c488/attachment.htm
 


[Pw_forum] Problem with Bi2Se3 SCF converge

2011-09-03 Thread WF
Dear Mohnish:
Thanks for our advices. I tried beta = 0.1 and 0.2, and those
message is postponed to 100+ iteration, but the convergence still cannot be
achieved.


On Fri, Sep 2, 2011 at 11:45 AM, mohnish pandey 
wrote:
>  Try reducing "beta"

>On Wed, Aug 31, 2011 at 8:07 AM, WF  wrote:

> Hello everyone,
>
>  I am doing calculation of Bi2Se3 with vdW-DF functional, but 
> it really hard to converge. After several iterations it keeps to 
> prompt things
> like:
>
> ** **
>
>  iteration # 58 ecut=60.00 Ry beta=0.70
>
>  CG style diagonalization
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  2 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  2 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  1 eigenvalues not converged
>
>  c_bands:  2 eigenvalues not converged
>
>  c_bands:  2 eigenvalues not converged
>
>  c_bands:  2 eigenvalues not converged
>
>  c_bands:  2 eigenvalues not converged
>
>  
>
> and will not converge in 1000 iteration.
>
> ** **
>
>  But the calculation of LDA and PBE can converge, however. All 
> pseudo-potential is generated by ld1.x with same parameters except
> functional.
>
> ** **
>
>  I have tried to change some parameter in SCF but it doesn?t 
> seem to be helpful. Is there any suggestion for this problem?
>
> ** **
>
>  Thanks.
>
> ** **
>
> F,Wu
>
> ** **
>
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>

-
Regards

---
F, Wu
College of Chemistry and Molecular Engineering
Peking University
--






[Pw_forum] total force and number of k-points

2011-09-03 Thread Elie Moujaes

Professor Giannozzi,
Thanks a lot for your reply.
Elie

> From: giannozz at democritos.it
> Date: Fri, 2 Sep 2011 22:55:25 +0200
> To: pw_forum at pwscf.org
> Subject: Re: [Pw_forum] total force and number of k-points
> 
> 
> On Sep 2, 2011, at 15:54 , Elie Moujaes wrote:
> 
> > I have a simple question: does the total force acting
> > on a system depend on the number of k-points  used?
> 
> not sure what you mean by "total force", but forces on
> atoms depend upon the number of k-points, for sure
> 
> > That is in my scf calculation I am using a grid of 4x4x1
> > for my system whereas the grid must be denser if phonon
> > calculations are taking place I am using 8x8x1. Do I get
> > different total force acting on the system and does it mean
> > I have to relax the crystal twice one for an scf calculation
> > and another for  a phonon calculation?
> 
> in principle a system should be relaxed for the same set of
> parameters (cutoff, k-points) that are used in the subsequent
> phonon calculation. The difference might be small enough to
> be negligible, but this is something that should be checked.
> 
> P.
> ---
> Paolo Giannozzi, Dept of Chemistry&Physics&Environment,
> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
> Phone +39-0432-558216, fax +39-0432-558222
> 
> 
> 
> 
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
  
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110903/fce43158/attachment.htm