
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">

<HTML><HEAD>

<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>

<META name=GENERATOR content="MSHTML 8.00.6001.19019">

<STYLE></STYLE>

</HEAD>

<BODY bgColor=#ffffff>

<DIV><FONT size=2 face=Arial>If I were you, I would try degauss=0.006, 0.004, 

0.002 with dense enough k points </FONT></DIV>

<DIV><FONT size=2 face=Arial>and see the tendency.</FONT></DIV>

<DIV><FONT size=2 face=Arial></FONT>&nbsp;</DIV>

<DIV><FONT size=2 face=Arial>zhou huiqun</FONT></DIV>

<DIV><FONT size=2 face=Arial>@earth sciences, nanjing university, 

china</FONT></DIV>

<BLOCKQUOTE 

style="BORDER-LEFT: #000000 2px solid; PADDING-LEFT: 5px; PADDING-RIGHT: 0px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px">

  <DIV style="FONT: 10pt arial">----- Original Message ----- </DIV>

  <DIV 

  style="FONT: 10pt arial; BACKGROUND: #e4e4e4; font-color: black"><B>From:</B> 

  <A title=tkaloni@gmail.com href="mailto:tkaloni@gmail.com">Thaneshwor 

  Kaloni</A> </DIV>

  <DIV style="FONT: 10pt arial"><B>To:</B> <A title=pw_forum@pwscf.org 

  href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A> </DIV>

  <DIV style="FONT: 10pt arial"><B>Sent:</B> Monday, March 14, 2011 12:22 

  AM</DIV>

  <DIV style="FONT: 10pt arial"><B>Subject:</B> [Pw_forum] problem with degauss 

  and magnetization</DIV>

  <DIV><BR></DIV>

  <DIV dir=ltr>Dear Quantum Espresso users,<BR><BR><BR>I am trying to get the 

  magnetization by creating one vacancy <BR>in 5x5 supercell of 

  graphene.<BR>With different values of degauss the total magnetization is 

  <BR>changing surprisingly. Is there anything wrong with my input ?.<BR>Could 

  anybody please tell me&nbsp; some references where I<BR>can read detail about 

  smearing and degauss ?.<BR><BR><BR>The total magnetization with&nbsp; 

  different degauss are given as 

  below,<BR><BR><B>degauss&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  total 

  magnetization<BR>&nbsp;0.05&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.00<BR>&nbsp;0.01&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  &nbsp; 

  0.58<BR>&nbsp;0.001&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.88<BR>&nbsp;0.0001&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.88<BR><BR>&nbsp; </B><BR><BR>I am posting my input file 

  below.<BR><BR><BR><BR>&amp;CONTROL<BR>calculation = 'scf' ,<BR>restart_mode = 

  'from_scratch' ,<BR>outdir = 'OUT/',<BR>pseudo_dir = 

  '/home/kaloni/pseudo/',<BR>prefix = 

  'g',<BR>tstress=.true.,<BR>/<BR>&amp;SYSTEM<BR>&nbsp;&nbsp;&nbsp; ibrav = 4, 

  a=12.199359575, b= 12.199359575,<BR>&nbsp;&nbsp;&nbsp; c=19.999668786, cosab= 

  -0.5, cosac =0.0 , cosbc=0.0,<BR>&nbsp;&nbsp;&nbsp; nat=&nbsp; 49, ntyp= 

  1,<BR>&nbsp;&nbsp;&nbsp; ecutwfc = 40,<BR>&nbsp;&nbsp;&nbsp; ecutrho = 

  400,<BR>&nbsp;&nbsp;&nbsp; nspin =2, starting_magnetization = 

  1.0,<BR>&nbsp;&nbsp;&nbsp; occupations='smearing', smearing=' gaussian', 

  degauss=0.05,<BR>&nbsp;&nbsp;&nbsp; nosym = 

  .true.,<BR><BR>/<BR>&nbsp;&amp;ELECTRONS<BR>&nbsp;&nbsp;&nbsp; mixing_beta = 

  0.3<BR>&nbsp;&nbsp;&nbsp; conv_thr = 

  1.0d-9<BR>&nbsp;/<BR>ATOMIC_SPECIES<BR>C&nbsp;&nbsp;&nbsp; 

  12.0107&nbsp;&nbsp;&nbsp;&nbsp; C.pz-rrkjus.UPF<BR>ATOMIC_POSITIONS 

  (crystal)<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.000437574&nbsp; 

  -0.000580938&nbsp;&nbsp; 

  0.000000015<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.000878501&nbsp;&nbsp; 0.200420327&nbsp; 

  -0.000000008<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.001210858&nbsp;&nbsp; 0.399950578&nbsp; 

  -0.000000071<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.000312205&nbsp;&nbsp; 

  0.599509065&nbsp;&nbsp; 0.000000002<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  -0.000716730&nbsp;&nbsp; 0.798976722&nbsp;&nbsp; 

  0.000000025<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.199479935&nbsp; 

  -0.000986699&nbsp; -0.000000031<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.200591611&nbsp;&nbsp; 0.200243356&nbsp;&nbsp; 

  0.000000020<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.204951004&nbsp;&nbsp; 0.399348478&nbsp; 

  -0.000000039<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.200756083&nbsp;&nbsp; 0.599210178&nbsp; 

  -0.000000027<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.199615366&nbsp;&nbsp; 0.799016660&nbsp;&nbsp; 

  0.000000045<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.399828248&nbsp; 

  -0.000558044&nbsp; -0.000000011<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.399613916&nbsp;&nbsp; 0.200154971&nbsp;&nbsp; 

  0.000000014<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.401785094&nbsp;&nbsp; 0.599046302&nbsp; 

  -0.000000047<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.400264884&nbsp;&nbsp; 0.799510584&nbsp;&nbsp; 

  0.000000019<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.600023737&nbsp;&nbsp; 0.000445858&nbsp; 

  -0.000000037<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.599500576&nbsp;&nbsp; 0.200387180&nbsp;&nbsp; 

  0.000000030<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.594152084&nbsp;&nbsp; 0.399266958&nbsp;&nbsp; 

  0.000000290<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.597008286&nbsp;&nbsp; 0.598889370&nbsp;&nbsp; 

  0.000000057<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.600448742&nbsp;&nbsp; 0.800956977&nbsp; 

  -0.000000036<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.800386518&nbsp;&nbsp; 0.000445469&nbsp; 

  -0.000000021<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.800217273&nbsp;&nbsp; 0.200473648&nbsp; 

  -0.000000180<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.798607049&nbsp;&nbsp; 0.399927530&nbsp;&nbsp; 

  0.000000120<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.798270617&nbsp;&nbsp; 0.599140827&nbsp;&nbsp; 

  0.000000073<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.799080030&nbsp;&nbsp; 0.799392409&nbsp; 

  -0.000000024<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.132951308&nbsp;&nbsp; 0.065395469&nbsp;&nbsp; 

  0.000000019<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.133142262&nbsp;&nbsp; 0.266027786&nbsp; 

  -0.000000014<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.135800585&nbsp;&nbsp; 0.466331814&nbsp; 

  -0.000000062<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.133194345&nbsp;&nbsp; 0.665806346&nbsp;&nbsp; 

  0.000000014<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.132633945&nbsp;&nbsp; 0.865364014&nbsp;&nbsp; 

  0.000000018<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.332386006&nbsp;&nbsp; 0.065362305&nbsp; 

  -0.000000012<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.329787525&nbsp;&nbsp; 0.259493844&nbsp;&nbsp; 

  0.000000003<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.336738546&nbsp;&nbsp; 0.468697263&nbsp; 

  -0.000000042<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.333671463&nbsp;&nbsp; 0.665846245&nbsp; 

  -0.000000020<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.333143624&nbsp;&nbsp; 0.865895587&nbsp;&nbsp; 

  0.000000044<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.533078996&nbsp;&nbsp; 0.066642011&nbsp; 

  -0.000000019<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.532824452&nbsp;&nbsp; 0.265959049&nbsp;&nbsp; 

  0.000000110<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.531539380&nbsp;&nbsp; 0.468474465&nbsp;&nbsp; 

  0.000000226<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.534652606&nbsp;&nbsp; 0.669416138&nbsp; 

  -0.000000013<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.533415252&nbsp;&nbsp; 0.867317865&nbsp; 

  -0.000000015<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.733533974&nbsp;&nbsp; 0.067101960&nbsp; 

  -0.000000071<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.733185901&nbsp;&nbsp; 0.266889900&nbsp; 

  -0.000000007<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.730336277&nbsp;&nbsp; 0.466262841&nbsp;&nbsp; 

  0.000000196<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.731928945&nbsp;&nbsp; 0.665704201&nbsp;&nbsp; 

  0.000000020<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.733842050&nbsp;&nbsp; 0.867300592&nbsp; 

  -0.000000017<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.933534511&nbsp;&nbsp; 0.066657775&nbsp;&nbsp; 

  0.000000016<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.933667182&nbsp;&nbsp; 0.266910925&nbsp; 

  -0.000000085<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.933067875&nbsp;&nbsp; 0.466263286&nbsp;&nbsp; 

  0.000000000<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.932454196&nbsp;&nbsp; 0.665766534&nbsp;&nbsp; 

  0.000000004<BR>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

  0.932632612&nbsp;&nbsp; 0.865835315&nbsp;&nbsp; 0.000000008<BR>K_POINTS 

  AUTOMATIC<BR>4 4 1 1 1 1<BR><BR><BR>Thanks <BR><BR>Thaneshwor Prashad 

  Kaloni<BR>&nbsp;<BR></DIV>

  <P>

  <HR>


  <P></P>_______________________________________________<BR>Pw_forum mailing 

  list<BR>Pw_forum@pwscf.org<BR>http://www.democritos.it/mailman/listinfo/pw_forum<BR></BLOCKQUOTE></BODY></HTML>