<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Dear Quantum Espresso users<br><br>I have done a 'vc-relax' calculation to find the cell parameters of graphite. The selected k-points are (4*4*4) and automatic.<br><br>obtained cell parameters:<br> a = 2.469,<br> b = 2.469,<br> c = 7.069,<br><br>then I used these values to 'relax' the atoms in cell (k-points were 4*4*4 too) . Finally I used the obtained cell parameters to find the proper k-points for graphite.<br>I performed some 'scf' calculations with different k-points and in whole of them the cell parameters and positions were what is obtained from the relax part with 4*4*4 k-points.<br><br><u><b>'scf' energy with different k-points:</b></u><br><br>k-points energy <br>
<meta http-equiv="CONTENT-TYPE" content="text/html; charset=utf-8">
<title></title>
<meta name="GENERATOR" content="OpenOffice.org 3.1 (Unix)">
<style type="text/css">
<!--
@page { margin: 0.79in }
P { margin-bottom: 0.08in }
-->
</style>
<p style="margin-bottom: 0in;">(4*4*4): -45.5732791969 Ry
</p>
<p style="margin-bottom: 0in;">(5*5*5): -45.57637408 Ry<br></p>
<p style="margin-bottom: 0in;">(6*6*6): -45.57592043 Ry<br></p>
<p style="margin-bottom: 0in;">(7*7*7): -45.57420434 Ry<br></p>
<p style="margin-bottom: 0in;">(8*8*8): -45.57474553 Ry<br></p>
<p style="margin-bottom: 0in;">(9*9*9): -45.57500776 Ry<br></p>
<p style="margin-bottom: 0in;">(10*10*10): -45.57454272 Ry<br></p>
<p style="margin-bottom: 0in;">(11*11*11): -45.57464463 Ry<br></p>
<p style="margin-bottom: 0in;">(12*12*12): -45.57480118 Ry<br></p>
<p style="margin-bottom: 0in;">(13*13*13): -45.57462800 Ry<br></p>
<p style="margin-bottom: 0in;">(14*14*14): -45.57463705 Ry<br></p>
<p style="margin-bottom: 0in;">(15*15*15): -45.57473135 Ry<br></p>
<p style="margin-bottom: 0in;">(16*16*16): -45.57465384 Ry</p><p style="margin-bottom: 0in;"><br></p><p style="margin-bottom: 0in;">I appreciate if one help me to find the answers to these questions:</p><p style="margin-bottom: 0in;"><br></p><p style="margin-bottom: 0in;">1. Is the whole procedure plausible ?</p><p style="margin-bottom: 0in;">2. why a energy convergence is not occurred after selecting these high values of k-points?</p><p style="margin-bottom: 0in;">3. How can I find the proper k-points safely?</p><p style="margin-bottom: 0in;"><br></p><p style="margin-bottom: 0in;"><b><u>the input file for 'scf' calculation:</u></b></p><p style="margin-bottom: 0in;"><br></p>&CONTROL<br><p style="margin-bottom: 0in;"> calculation = "scf",<br> pseudo_dir = "/home/koa/soft/qe4.2/espresso-4.2/pseudo",<br> outdir = "/home/koa/tmp",<br> etot_conv_thr= 1.0D-4,<br> forc_conv_thr= 1.0D-3,<br> nstep=150,<br> /<br>&SYSTEM<br> ibrav = 4,<br> a = 2.469,<br> b = 2.469,<br> c = 7.069,<br> cosab = -0.5,<br> cosac = 1.0,<br> cosbc = 1.0,<br> nat = 4,<br> ntyp = 1,<br> ecutwfc = 30.D0,<br> ecutrho = 200.D0<br> occupations = 'smearing'<br> smearing ='mp', <br> degauss = 0.01,<br> nspin = 2,<br> starting_magnetization(1)= 1.0,<br> /<br>&ELECTRONS<br> conv_thr = 1.D-6,<br> mixing_beta = 0.7D0,<br> diagonalization = "david",<br>/<br>&IONS<br> ion_dynamics="bfgs"<br>/<br>ATOMIC_SPECIES<br>C 12.0107 C.pbe-van_ak.UPF<br>ATOMIC_POSITIONS {angstrom}<br>C 0.000000000 0.000000000 0.000000000<br>C 0.000000000 1.425534302 0.000000000<br>C 0.000000000 0.000000000 3.534553118<br>C 1.234548920 0.712767151 3.534553118
<br>K_POINTS {automatic}<br>16 16 16 1 1 1<br><br><b><u>the vc-relax input: </u></b></p><p style="margin-bottom: 0in;"><br></p><p style="margin-bottom: 0in;">&CONTROL<br> calculation = "vc-relax",<br> pseudo_dir = "/home/koa/soft/qe4.2/espresso-4.2/pseudo",<br> outdir = "/home/koa/tmp",<br> etot_conv_thr= 1.0D-4,<br> forc_conv_thr= 1.0D-3,<br> dt=70,<br> nstep=150,<br> /<br>&SYSTEM<br> ibrav = 4,<br> a = 2.456,<br> b = 2.456,<br> c = 6.696,<br> cosab = -0.5,<br> cosac = 1.0,<br> cosbc = 1.0,<br> nat = 4,<br> ntyp = 1,<br> ecutwfc = 30.D0,<br> ecutrho = 200.D0<br> occupations = 'smearing'<br> smearing ='mp', <br> degauss = 0.01,<br> nspin = 2,<br> starting_magnetization(1)= 1.0,<br> /<br>&ELECTRONS<br> conv_thr = 1.D-6,<br> mixing_beta = 0.7D0,<br> diagonalization = "david",<br>/<br>&IONS<br> ion_dynamics="bfgs"<br>/<br>&CELL<br>cell_dynamics = 'damp-w',<br>press = 0.0 ,<br>/<br>ATOMIC_SPECIES<br>C 12.0107 C.pbe-van_ak.UPF<br>ATOMIC_POSITIONS {angstrom}<br> C 0.00000000 0.00000000 0.00000000
<br> C 0.00000000 1.41797226 0.00000000
<br> C 0.00000000 0.00000000 3.34800000
<br> C 1.22800000 0.70898613 3.34800000
<br>K_POINTS {automatic}<br>4 4 4 1 1 1<br></p><p style="margin-bottom: 0in;"><br></p><p style="margin-bottom: 0in;">very best wishes<br></p><p style="margin-bottom: 0in;"><br>Many Thanks</p><p style="margin-bottom: 0in;">Masoud Nahali</p><p style="margin-bottom: 0in;">Sharif University of Technology<br></p><p style="margin-bottom: 0in;"><br></p>
<br /><hr />Hotmail: Trusted email with Microsoft’s powerful SPAM protection. <a href='https://signup.live.com/signup.aspx?id=60969' target='_new'>Sign up now.</a></body>
</html>