<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Dear Baroni,<BR>
thank you, but it seems i am not using US pps. However, I will increase ecutrho and see.<BR>
thanks again<BR>
hania<BR><BR><BR>
<BLOCKQUOTE>
<HR>
From: baroni@sissa.it<BR>To: pw_forum@pwscf.org<BR>Date: Sun, 22 Jun 2008 14:32:02 +0200<BR>Subject: Re: [Pw_forum] cell optimization<BR><BR>I may be wrong, but your problem could be due to subtelties of the US PP formalism. Using US PPs, I am afraid that PW calculations are variational with respect to PW expansion coefficients only when the augmentation charge is at convergence. Convergence with respect to the augmentation charge usually is a rather delicate issue. For sure, the cutoff for the charge-density distribution (ecutrho) should be taken larger (sometimes much larger) than the default value for norm-conserving PPs (4 x ecutwfc). In yor calculation, you seem to use just this value, which may turn out to be inadequate. Have you tried to test the convergence of your calcs wrt echutrho (keeping ecutwfc fixed)? SB
<DIV><BR>
<DIV>
<DIV>On Jun 22, 2008, at 11:06 AM, hania djani-ait aissa wrote:</DIV><BR class=EC_Apple-interchange-newline>
<BLOCKQUOTE><SPAN class=EC_Apple-style-span style="WORD-SPACING: 0px; FONT: 14px Helvetica; TEXT-TRANSFORM: none; COLOR: rgb(0,0,0); TEXT-INDENT: 0px; WHITE-SPACE: normal; LETTER-SPACING: normal; BORDER-COLLAPSE: separate">
<DIV class=EC_hmmessage style="FONT-SIZE: 10pt; FONT-FAMILY: Tahoma">Dear Nicola,<BR>I have tried startingwfc='random' and it gave me the same energies values as for startingwfc='atomic' with still the same problem of drop and rise of total energy. I also increase Ecutoff to 40Ry, the energies values are actually more accurate but the problem persist.<BR> <BR>Hania<BR> <BR> <BR>
<DIV style="PADDING-RIGHT: 0px; PADDING-LEFT: 0px; MARGIN-BOTTOM: 0px; PADDING-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px; PADDING-TOP: 0px"><SPAN lang=EN-GB><FONT color=#000000><FONT face="Courier New"><SPAN> </SPAN>&control<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>calculation<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>= 'scf'<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>restart_mode = 'from_scratch'<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>pseudo_dir<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>= '/root/pwscf/pseudo/'<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>outdir<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>= '/root/tmp/'<BR><SPAN> </SPAN>/<BR><SPAN> </SPAN>&system<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>ibrav=7<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>celldm(1)=7.33, celldm(3)=8.54<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>nat=19<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>ntyp=3<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>ecutwfc=25.0<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>ecutrho=100.0<BR><SPAN> </SPAN>/<BR><SPAN> </SPAN>&electrons<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>conv_thr = 1e-6,<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>mixing_beta=0.2,</FONT></FONT></SPAN></DIV>
<DIV style="PADDING-RIGHT: 0px; PADDING-LEFT: 0px; MARGIN-BOTTOM: 0px; PADDING-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px; PADDING-TOP: 0px"><SPAN lang=EN-GB><FONT color=#000000><FONT face="Courier New"><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>startingwfc=’random’<BR><SPAN> </SPAN>/<BR>ATOMIC_SPECIES<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>Bi<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>208.98<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>Bi.pbe-d-mt.UPF<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN></FONT></FONT></SPAN><FONT color=#000000><FONT face="Courier New"><SPAN lang=IT>Ti<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>47.867<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>Ti.pbe-sp-van_ak.UPF<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>15.9994<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>O.pbe-van_ak.UPF<BR></SPAN>ATOMIC_POSITIONS<SPAN class=EC_Apple-converted-space> </SPAN><BR><SPAN> </SPAN>Bi 0.00000000 0.00000000 0.5713404<BR><SPAN> </SPAN>Bi 0.00000000 0.00000000 7.9816596<BR><SPAN> </SPAN>Bi 0.00000000 0.00000000 1.8072489<BR><SPAN> </SPAN>Bi 0.00000000 0.00000000 6.7457511<BR><SPAN> </SPAN></FONT></FONT><FONT color=#000000><FONT face="Courier New"><SPAN lang=IT>Ti 0.00000000 0.00000000 3.1757289<BR><SPAN> </SPAN>Ti 0.00000000 0.00000000 5.3772711<BR><SPAN> </SPAN>Ti 0.00000000 0.00000000 4.2765000<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.50000000 0.00000000 0.0000<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.00000000 0.50000000 0.0000<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.50000000 0.00000000 2.13825<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.50000000 0.00000000 6.41475<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.00000000 0.00000000 3.771873<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.00000000 0.00000000 4.781127<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.00000000 0.00000000 2.7241305<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.00000000 0.00000000 5.8288695<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.50000000 0.00000000 1.0049775<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.50000000 0.00000000 7.5480225<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.00000000 0.50000000 1.0049775<BR><SPAN> </SPAN>O<SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN>0.00000000 0.50000000 7.5480225<BR>K_POINTS {automatic}<BR><SPAN> <SPAN class=EC_Apple-converted-space> </SPAN></SPAN></SPAN>4 4 4 1 1 1<BR> </FONT></FONT></DIV> <BR><BR>
<HR>
Téléchargez le nouveau Windows Live Messenger !<SPAN class=EC_Apple-converted-space> </SPAN><A href="http://get.live.com/messenger/overview" target=_blank>Téléchargez Messenger, c'est gratuit !</A>_______________________________________________<BR>Pw_forum mailing list<BR><A href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</A><BR><A href="http://www.democritos.it/mailman/listinfo/pw_forum" target=_blank>http://www.democritos.it/mailman/listinfo/pw_forum</A><BR></DIV></SPAN></BLOCKQUOTE></DIV><BR></DIV><BR><BR>
<DIV><SPAN class=EC_Apple-style-span style="WORD-SPACING: 0px; FONT: 14px Helvetica; TEXT-TRANSFORM: none; COLOR: rgb(0,0,0); TEXT-INDENT: 0px; WHITE-SPACE: normal; LETTER-SPACING: normal; BORDER-COLLAPSE: separate">
<DIV style="WORD-WRAP: break-word"><SPAN class=EC_Apple-style-span style="WORD-SPACING: 0px; FONT: 14px Helvetica; TEXT-TRANSFORM: none; COLOR: rgb(0,0,0); TEXT-INDENT: 0px; WHITE-SPACE: normal; LETTER-SPACING: normal; BORDER-COLLAPSE: separate">
<DIV style="WORD-WRAP: break-word"><SPAN class=EC_Apple-style-span style="WORD-SPACING: 0px; FONT: 14px Helvetica; TEXT-TRANSFORM: none; COLOR: rgb(0,0,0); TEXT-INDENT: 0px; WHITE-SPACE: normal; LETTER-SPACING: normal; BORDER-COLLAPSE: separate">
<DIV style="WORD-WRAP: break-word"><SPAN class=EC_Apple-style-span style="WORD-SPACING: 0px; FONT: 14px Helvetica; TEXT-TRANSFORM: none; COLOR: rgb(0,0,0); TEXT-INDENT: 0px; WHITE-SPACE: normal; LETTER-SPACING: normal; BORDER-COLLAPSE: separate">
<DIV style="WORD-WRAP: break-word">
<DIV>
<DIV style="MARGIN-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">---</SPAN></FONT></DIV>
<DIV style="MARGIN-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">Stefano Baroni - SISSA</SPAN></FONT><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"> </SPAN></FONT><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">&</SPAN></FONT><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"> </SPAN></FONT><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">DEMOCRITOS National Simulation Center - Trieste</SPAN></FONT></DIV>
<DIV style="MARGIN-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"><A href="http://www.sissa.it/~baroni" target=_blank>http://www.sissa.it/~baroni</A> / [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)</SPAN></FONT></DIV>
<DIV style="MIN-HEIGHT: 14px; MARGIN-BOTTOM: 0px; FONT: 12px Helvetica; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><BR style="FONT-SIZE: 12px"></DIV>
<DIV style="MIN-HEIGHT: 14px; MARGIN-BOTTOM: 0px; FONT: 12px Helvetica; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><SPAN class=EC_Apple-style-span style="FONT-SIZE: 14px">
<DIV style="MARGIN-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">La morale est une logique de l'action comme la logique est une morale de la pensée - Jean Piaget</SPAN></FONT></DIV>
<DIV style="MARGIN-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"><BR></SPAN></FONT></DIV></SPAN></DIV>
<DIV style="MARGIN-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">Please, if possible, don't</SPAN></FONT><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"> </SPAN></FONT><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">send me MS Word or PowerPoint attachments</SPAN></FONT></DIV>
<DIV style="MARGIN-BOTTOM: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px">Why? See:</SPAN></FONT><FONT class=EC_Apple-style-span size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"> </SPAN></FONT><FONT class=EC_Apple-style-span color=#0023e9 size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"><A href="http://www.gnu.org/philosophy/no-word-attachments.html" target=_blank>http://www.gnu.org/philosophy/no-word-attachments.html</A></SPAN></FONT></DIV>
<DIV><FONT class=EC_Apple-style-span color=#0023e9 size=3><SPAN class=EC_Apple-style-span style="FONT-SIZE: 12px"><BR class=EC_webkit-block-placeholder></SPAN></FONT></DIV></DIV></DIV></SPAN><BR class=EC_Apple-interchange-newline></DIV></SPAN><BR class=EC_Apple-interchange-newline></DIV></SPAN><BR class=EC_Apple-interchange-newline></DIV></SPAN><BR class=EC_Apple-interchange-newline></DIV><BR></BLOCKQUOTE><br /><hr />Téléchargez le nouveau Windows Live Messenger ! <a href='http://get.live.com/messenger/overview' target='_new'>Téléchargez Messenger, c'est gratuit !</a></body>
</html>