<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt">
<div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">
<div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">Hi, <br><br>Use attached ev.f90 instead. I think output file&nbsp; is quite self-explanatory.<br>All you need is just answering correctly questions requested by the program (ev.x).<br><br>Besides, you need to learn about the methods to calculate ground state parameters, including bulk modulus. See, at least, VLab lecture given by&nbsp; R. Wentzcovitch , <br><span><a target="_blank" href="http://www.vlab.msi.umn.edu/events/download/vlab_lectures/renata/lecture3.pdf">http://www.vlab.msi.umn.edu/events/download/vlab_lectures/renata/lecture3.pdf</a></span><br><br>I also recommend to search this site to find very good lectures by&nbsp; another authors.<br><br>I&nbsp; remember, last year there was a PRB paper concerning applications of different methods <br>(<font size="2">Phys. Rev. B
                      78,
                      214108
                      (2008) </font>).<br><br>Bests,<br>Eyvaz<br>-------------------------------------------------------------------<br>Prof. Eyvaz Isaev, <br>Department of Physics, Chemistry, and Biology (IFM), Linkoping University, Sweden <br>Theoretical Physics Department, Moscow State Institute of Steel &amp; Alloys, Russia, <br>isaev@ifm.liu.se, eyvaz_isaev@yahoo.com<div><br></div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font face="Tahoma" size="2"><hr size="1"><b><span style="font-weight: bold;">From:</span></b> Tram Bui
 &lt;trambui@u.boisestate.edu&gt;<br><b><span style="font-weight: bold;">To:</span></b> PWSCF Forum &lt;pw_forum@pwscf.org&gt;<br><b><span style="font-weight: bold;">Sent:</span></b> Tue, January 4, 2011 11:06:44 PM<br><b><span style="font-weight: bold;">Subject:</span></b> [Pw_forum] Bulk Modulus Calculation<br></font><br>
Hi All,<br>&nbsp;&nbsp;&nbsp;&nbsp; Hope everyone had great holidays, I had a question before regarding the Bulk Modulus calculation (ev.f90) But I'm having hard time understand&nbsp; the file and don't know how to use the ev.f to help me calculate the bulk modulus if I have an input file is like below and the structure is zinc blend (b-SiC). Would you help me if you get a chance? here is my file and you can run it in xcyrsden to see the structure.<br>
&nbsp;<br>&amp;CONTROL<br>&nbsp;&nbsp; calculation='scf'<br>&nbsp;&nbsp; restart_mode='from_scratch',<br>&nbsp;&nbsp; prefix='Silicon',<br>&nbsp;&nbsp; tprnfor=.true.<br>&nbsp;&nbsp; pseudo_dir='/home/trambui/QE/espresso-4.2.1/pseudo',<br>&nbsp;&nbsp; outdir='./tmp7/',<br>
/<br>&amp;SYSTEM<br>&nbsp;&nbsp; ibrav=2, celldm(1)=8.33, nat=2, ntyp=2,<br>&nbsp;&nbsp; ecutwfc= 30,<br>/<br>&amp;ELECTRONS<br>&nbsp;&nbsp; diagonalization= 'david',<br>&nbsp;&nbsp; mixing_mode= 'plain',<br>&nbsp;&nbsp; mixing_beta= 0.5,<br>&nbsp;&nbsp; conv_thr= 1.0d-7,<br>
/<br>ATOMIC_SPECIES<br>Si 28.086 Si.pbe-n-van.UPF<br>&nbsp;C 12.011 C.pbe-van_ak.UPF<br>ATOMIC_POSITIONS<br>Si 0.00 0.00 0.00<br>&nbsp;C 0.25 0.25 0.25<br>K_POINTS automatic<br>10 10 10 0 0 0<br><br>I really appreciate your helps!<br clear="all">
<br>Tram Bui<br><br>B.S. Materials Science &amp; Engineering<br><a rel="nofollow" ymailto="mailto:trambui@u.boisestate.edu" target="_blank" href="mailto:trambui@u.boisestate.edu">trambui@u.boisestate.edu</a><br><br>
</div></div>
</div></div>
</div><br>

      </body></html>