[Pw_forum] Pd pseudopotential

2011-08-08 Thread Christopher Heard
Hello all,

I want to do some structure relaxations on Pd particles, and I noticed all the 
available pseudopotentials have electronic configuration d9s1.  I would have 
thought d10 is more likely, and am worried as to whether I will therefore have 
to treat all calculations as openshell, with either a fixed tot_magnetization, 
or worse, with full magnetic convergence with a guess at starting_magnetization.

Both will slow me down a lot, and I'd rather not have to.  
Does anyone know, firstly if indeed I will have to do these magnetic 
calculations, and secondly, why d9s1 should be more stable than closed shell?

Thanks,
Chris 

---
PhD student,
University of Birmingham, UK
--- 


[Pw_forum] Slow calculation of (5, 5) CNT transmission and complex bands.

2011-08-08 Thread chengyu yang
Hi all.
As an update I calculated the bandstructure assuming that the negative rho
has little effect on the final result, and then I got a bandstructure with a
certain lines in it. I am hoping if anyone can tell me how to get all the k
point for the nscf input, so I can get a more smooth bandstructure finally.
Thank you.

Best regards.

Chengyu Yang
MMAE, University of Central Florida
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110808/12326d11/attachment.htm
 


[Pw_forum] Co Pseudopotential

2011-08-08 Thread jia chen
Hi,

If I understand you correctly, there is electron transfer between the
surface and the Co atom which makes the spin polarization of Co atom to be
1. Not very surprising to me.

On Mon, Aug 8, 2011 at 7:16 PM, Izaak Williamson  wrote:

> Thank you Heather and Jia for your replies. For a single Co atom, I do get
> the total magnetization of 3 but when Co is on a surface, its spin
> polarization comes out to be 2, which is S=1. The all-electron calculations
> were also done using pbe functional and on the same surface and with the
> same coverage of Co.
>
> Heather, what do you mean by locally? I do calculate the spin polarization
> from pdos. The whole system does have an even number of electrons but why
> does it result in S=1? Any explanation would be appreciated.
>
>
> --
> Izaak Williamson
> Research Assistant
> Physics Department
> Boise State University
>
>
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>


-- 
Jia Chen
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110808/b2c4bb69/attachment.htm
 


[Pw_forum] vdw-df2

2011-08-08 Thread Robin H
I think you can use revPBE to produce the PPs for using,and in your
calculation you can set the input_dft='vDW-DF',that's what the example
mean.Could you tell me where you find the information of the example?
Thanks.
Robin H

2011/7/28 DONG Rui 

> Hi, all
>
> I have something not clear about the vdw-df implementation.
>
> It is mentioned in the example that we'd better use revised pbe when using
> the vdw-df, but "revPBE" itself is another xc-functional? What do we do now?
>
> Do we write the input_dft="revPBE+vdW-DF2"?
>
> or we do "revPBE" first and use the result to restart a "vdW-DF2"
> calculation?
>
> Thank you very much for help!
>
> --
> Sincerely,
>
> Rui DONG
>
> ___
> 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/20110808/a0acbb89/attachment.htm
 


[Pw_forum] Co Pseudopotential

2011-08-08 Thread heather kulik
Total magnetization as printed in the code is equivalent to nelup-neldw. You 
cannot get an odd total magnetization if you have an even number of electrons. 

The total magnetization is only the same as the magnetization on cobalt if 
cobalt is the only atom in your simulation cell with a magnetic moment. 

If, as you have indicated, you have created a surface with an odd number of 
electrons, it must also have a net spin or, alternatively, charge transfer may 
occur between Co and the surface so that the surface has no net spin. Either 
way, mulliken population analysis, projections onto a localized basis, or some 
alternative approach is required to approximately identify the magnetic moment 
on each atom.

If you provide the reference you're comparing to, it might be clearer what you 
are trying to compare against.

Hope that helps,

Heather Kulik
Stanford University
On Monday, August 8, 2011 at 4:16 PM, Izaak Williamson wrote: 
> Thank you Heather and Jia for your replies. For a single Co atom, I do get 
> the total magnetization of 3 but when Co is on a surface, its spin 
> polarization comes out to be 2, which is S=1. The all-electron calculations 
> were also done using pbe functional and on the same surface and with the same 
> coverage of Co.
> 
> Heather, what do you mean by locally? I do calculate the spin polarization 
> from pdos. The whole system does have an even number of electrons but why 
> does it result in S=1? Any explanation would be appreciated.
> 
> -- 
> Izaak Williamson
> Research Assistant
> Physics Department
> Boise State University
> 
> ___
> 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/20110808/a3ebde1e/attachment.htm
 


[Pw_forum] Co Pseudopotential

2011-08-08 Thread jia chen
Hi,

I did lots of calculations by this one, works decently and I did atomic
calculation, total magnetization is 3.0.

On Mon, Aug 8, 2011 at 2:06 PM, Izaak Williamson  wrote:

> Hello everyone,
>
> I am trying to calculate the magnetic spin of Co adsorbed onto a surface.
> The literature that I am comparing my values with says the Co spin is 3/2
> from an all-electron calculation and I am getting S=1. The pseudopotential
> that I am using for Co is Co.pbe-nd-rrkjus.UPF. Is there something about
> this pseudopotential that is giving me an incorrect spin? Should I be using
> a different pseudopotential?
>
> Thanks for any help.
>
> --
> Izaak Williamson
> Research Assistant
> Physics Department
> Boise State University
>
>
>
>
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>


-- 
Jia Chen
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110808/76e3c5ca/attachment.htm
 


[Pw_forum] Co Pseudopotential

2011-08-08 Thread Izaak Williamson
Thank
you Heather and Jia for your replies. For a single Co atom, I do get the
total magnetization of 3 but when Co is on a surface, its spin
polarization comes out to be 2, which is S=1. The all-electron
calculations were also done using pbe functional and on the same
surface and with the same coverage of Co.

Heather, what do you mean by
locally? I do calculate the spin polarization from pdos. The whole
system does have an even number of electrons but why does it result in
S=1? Any explanation would be appreciated.

-- 
Izaak Williamson
Research Assistant
Physics Department
Boise State University
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110808/8eb63a4e/attachment.htm
 


[Pw_forum] Slow calculation of (5, 5) CNT transmission and complex bands.

2011-08-08 Thread chengyu yang
Thank you Dr.Sclauzero and all,

 I've changed the pseudo potential to USPP, set the ecutwfc down to 35,
and reduce theepsproj=1.d-5, and I ran the whole procedure to calculate the
band structure of my (5,5) CNT again.

   This time the procedure runs well and faster  until I got to the bands.x
step, I got " negative rho (up, down):  0.502E-03 0.000E+00" at the
bands.out, and it stopped there. I couldn't proceed further.

My another question is about the K points, I could generate a k point
path via Xcrysden, but I couldn't generate the all the k points I need to
input manually for the 'nscf'  input as in the examples. I defined 50 k
points in Xcrysden, however, the final result was only the few high symmetry
points I defined. I don't know how to get all the k points the 'nscf' needs?
I've tried the kpoints.x, but I don't know how the mesh n1,n2,n3 can be
connected to the k point path by Xcrysden, and how I should set up them.

  Thanks beforehand for anyone who can explain these things to me, I
really appreciate it.

  Best regards.

Chengyu Yang
Material Science, University of Central Florida
USA






2011/8/2 Gabriele Sclauzero 

>
> Il giorno 29/lug/2011, alle ore 21.16, chengyu yang ha scritto:
>
> Dear All,
>  After two hours , my transmission calculation comes out a result.
>
>
> And you couldn't wait for a couple of hours before writing to the forum
> saying that you were in trouble?...
>
> Assuming that your PWscf calculation has been converged, the important
> thing here is to look at the basis set reduction performed by PWcond:
>
>  ngper, shell number = 2809 299
>  ngper, n2d = 2809 472
>
> It looks to me that the n2d could be too small compared to the original
> size of the basis set (ngper). As I wrote in my previous email, you should
> try to decrease epsproj and check that the results (i.e. the complex bands)
> do not change appreciably. The eigenvalues at real k from PWcond need to be
> on the top of the band structure computed with PWscf.
> Also compare with previous theoretical results, I guess that there should
> be several works about the conductance of CNTs.
> cuplot is needed only to plot the scattering states. I guess that you do
> not need them at this stage.
>
> HTH
>
> GS
>
> But I thought it may be not converged and didn't get the final result.
> May anybody have a look at my output file? I attached in this email.
>My input file is like this:
>&inputcond
> prefixl='55'
> tran_file='trans.cnt55',
> ikind=0
> energy0=4.d0
> denergy=-0.1d0
> ewind=3.d0
> epsproj=1.d-4
> delgep = 1.d-9
> cutplot = 3.d0
>  /
> 1
> 0.0 0.0 1.0
> 5
>
>
> Is there somebody that can give me an advice?
> I really appreciate it.
>
> Chengyu Yang
> Materials Science and Engineering,
> University of Central Florida, USA
>
> ___
>
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>
>
> ? Gabriele Sclauzero, EPFL SB ITP CSEA
> *   PH H2 462, Station 3, CH-1015 Lausanne*
>
>
>
>
>
>
>
>
> ___
> 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/20110808/2e06099f/attachment-0001.htm
 


[Pw_forum] Run example error in qe-gipaw

2011-08-08 Thread Layla Martin-Samos
Dear Pengju, I just try qe-gipaw example without finding any problem. One of
the thing that could give some difficulties is the fact that the out dir is
set to /tmp/ceresoli. Maybe you could try to change outdir to something
else.

bests regards

Layla
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110808/519a68c1/attachment.htm
 


[Pw_forum] QE v4.3.2 failed to read the input file of v4.3

2011-08-08 Thread Layla Martin-Samos
Dear lilonghua could you send please the input file that causes the problem
you mentioned?

thank you very much

Layla

2011/8/6 lilonghua 

> **
> I think the problem is caused by the character length. when the input
> character len<=79, everything is all right!
>
>
> 2011-08-06
> --
>  lilonghua
> --
> ** lilonghua
> *?* 2011-08-06  11:36:23
> ** pw_forum
> *???*
> *???* QE v4.3.2 failed to read the input file of v4.3
>   Dear all,
>
> I have compiled the lastest version of QE. However, pw.x v4.3.2 failed to 
> read the input file of v4.3, even the EXX_example in v4.3.2. the Error 
> message is:
>   " from  read_namelists  : error #17
>   reading namelist system"
>
> I don't find significant changes of variables in document of pwscf v4.3.2. 
> And pwscf v4.3 is ok to perform the EXX_example of v4.3.2 (actually, those 
> examples are same to that of v4.3). So, what's happen to the io of v4.3.2? 
> and how to solve it?
>Best regards,
> Longhua Li
> 2011-08-06
> 
> lilonghua
>
> ___
> 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/20110808/06d03a5b/attachment.htm
 


[Pw_forum] Co Pseudopotential

2011-08-08 Thread heather kulik
Izaak,

If you really are getting S=1 for the total spin in your calculation, then you 
have an even number of electrons in your system. Neutral cobalt has an odd 
number of valence electrons, 3d74s2, (hence the S=3/2) so this likely means 
that you have an odd number of electrons in whatever you chose for your 
surface. 

That's not to say you couldn't run into other problems. If you mean you see S=1 
locally, e.g. from PDOS, you are likely getting charge transfer between the 
atom and the surface, and you may observe something that looks like S=1 on Co. 
I would need more information about your "all electron calculation" reference 
from the literature (i.e. is it Co on the same surface and with the same 
functional?) in order to suggest anything further.

Hope that helps,
Heather Kulik
Stanford University
On Monday, August 8, 2011 at 11:06 AM, Izaak Williamson wrote: 
> Hello everyone,
> 
> I am trying to calculate the magnetic spin of Co adsorbed onto a surface. The 
> literature that I am comparing my values with says the Co spin is 3/2 from an 
> all-electron calculation and I am getting S=1. The pseudopotential that I am 
> using for Co is Co.pbe-nd-rrkjus.UPF. Is there something about this 
> pseudopotential that is giving me an incorrect spin? Should I be using a 
> different pseudopotential?
> 
> Thanks for any help.
> 
> -- 
> Izaak Williamson
> Research Assistant
> Physics Department
> Boise State University
> 
> 
> 
> ___
> 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/20110808/2398fe74/attachment.htm
 


[Pw_forum] Co Pseudopotential

2011-08-08 Thread Izaak Williamson
Hello everyone,


I am trying to calculate the magnetic spin of Co adsorbed onto a surface. The 
literature that I am comparing my values with says the Co spin?is 3/2 from an 
all-electron calculation and I am getting S=1. The pseudopotential that I am 
using for Co is Co.pbe-nd-rrkjus.UPF. Is there something about this 
pseudopotential that is giving me an incorrect spin? Should I be using a 
different pseudopotential?

Thanks for any help.

-- 
 Izaak Williamson
Research Assistant
Physics Department
Boise State University
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110808/d4c5398b/attachment.htm