
<html>

<head>

<style>

.hmmessage P

{

margin:0px;

padding:0px

}

body.hmmessage

{

font-size: 10pt;

font-family:Verdana

}

</style>

</head>

<body class='hmmessage'>

Dear Tim,<BR>here attached my grep on volume.<BR>

cheers<BR>

Hania Djani ait aissa<BR>

Algeria<BR>> Date: Sat, 24 Jan 2009 10:36:57 -0800<BR>> From: thmason24@yahoo.com<BR>> To: pw_forum@pwscf.org<BR>> CC: akohlmey@cmm.chem.upenn.edu<BR>> Subject: Re: [Pw_forum] vc-relax<BR>> <BR>> Axel,<BR>> <BR>> thanks for the insight. That is outstanding information. I expect many people are running into this problem (I've seen the force oscillation energy decrease phenomenom posted many times). I am using what I thought were conservative ecuts (40 ryd) but perhaps not big enough for the stress tensor when using USPP. This characteristic would be a great thing to warn people about in the wiki vc-relax section. I'd offer to do it myself but I'd feel better leaving it to an expert. I'll report back if this improves my calculation but with the higher ecuts it won't be for a while. Is the issue confined to USPP? would you recommend NCPP for VC-relax as a work around?<BR>> <BR>> <BR>> <BR>> <BR>> Tim Mason<BR>> University of Missouri St. Louis<BR>> <BR>> <BR>> <BR>> <BR>> <BR>> <BR>> On Sat, Jan 24, 2009 at 12:12 PM, Axel Kohlmeyer <akohlmey@cmm.chem.upenn.edu> wrote:<BR>> <BR>> On Sat, 24 Jan 2009, Timothy Mason wrote:<BR>> <BR>> TM> Dear Hania,<BR>> TM><BR>> <BR>> TM> as someone else struggling with VC-relax, I'd like to see a grep of<BR>> TM> your cell volumes. (just grep on volume). I find that in these<BR>> TM> cases where the forces oscillate but the engergies decrease (all<BR>> TM> cases for me) the volumes continue to expand about linearly with<BR>> TM> time. I'm not sure what to make of it. It even happens starting<BR>> <BR>> good point. that could very likely mean, that your stress<BR>> tensor is not sufficiently converged wrt. to the wavefunction<BR>> or density cutoff. one of the most problematic issues with<BR>> vc-relax calculations is that the stress tensor need a much<BR>> higher wavefunction cutoff (and very much highe density cutoff<BR>> in case of ultra-soft pseudos) to be well converged. to make<BR>> matters worse, in a vc-relax calculations, you do not keep<BR>> those cutoffs constant, but the number of plane waves. so for<BR>> an expanding cell you effectively decrease the cutoffs, and<BR>> if you stress tensor is not well converged, it will become<BR>> quickly not converged. there are some tricks to make this<BR>> effect less of a problem, but the immediate solution is to<BR>> start at a lattice constant that is definitely too large,<BR>> i.e. the cell will never expand beyond it. you still have<BR>> consider pulay stress, but that is relatively minor concern.<BR>> <BR>> TM> with experimental structures often increasing them by 50 percent<BR>> TM> before convergence thresholds are met.<BR>> <BR>> a 50% increase sounds unphysical. have you tried restarting<BR>> from the final configuration (resets the "effective" basis<BR>> set size) and see if the results are consistent?<BR>> <BR>> cheers,<BR>> axel.<BR>> <BR>> <BR>> _______________________________________________<BR>> Pw_forum mailing list<BR>> Pw_forum@pwscf.org<BR>> http://www.democritos.it/mailman/listinfo/pw_forum<BR><BR><br /><hr />Découvrez toutes les possibilités de communication <a href='http://www.microsoft.com/windows/windowslive/default.aspx' target='_new'>avec vos proches</a></body>

</html>