<div dir="ltr"><pre>Kondrin, thank you so much for your quick response. But I want to discuss <br>some more things. I need only bulk modulus not shear modulus. Instead I need deformation potential also. <br>So your suggestion of calculating the bulk modulus from the energy difference <br>between the optimized structure at normal and high pressure can resolve only the bulk modulus problem.<br></pre><pre>But can you suggest me something about deformation potential. For my bulk modulus calculation, initially <br> I have used 'vc-relaxd' structure so that I know that my structure is optimized.<br></pre><pre>Then after deforming the lattice I ran only 'scf' calculations. should I still 'relax' after deforming every<br></pre><pre>time. Also can you suggest me some reference relating elastic constants to bulk modulus for <br></pre><pre>triclinic structure.<br></pre></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 15, 2015 at 4:53 PM, VARSHA RANI <span dir="ltr"><<a href="mailto:varsha24_sps@jnu.ac.in" target="_blank">varsha24_sps@jnu.ac.in</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear all,</div><div><br></div><div>I am trying to calculate the bulk modulus of single crystal Pentacene which has a triclinic unit cell with 72 atoms.  Starting from an optimized structure, with ibrav = 14 (corresponding to triclinc structure), I performed a series of SCF calculations. I changed one of the lattice constants by modifying celldm(1), celldm(2) and celldm(3) appropriately. I then used ev.x to extract the bulk modulus (with the non-cubic option). I also tried to calculate it directly from the parabolic fitting of energy vs volume.</div><div><br></div><div>From both approaches, I am getting a bulk modulus ~ 350 GPa, which seems unrealistic. I have verified that I am able to get the right value of the Bulk modulus for Silicon using both approaches. </div><div><br></div><div>I was wondering if there is a different way to go about calculating bulk modulus values for triclinic structure, or if there is something wrong about the way I have set up my input files. I have attached my input files for two different configurations, the equilibrium lattice constant and the strained value where celldm(1) = a, has been changed by 5%.</div><div><br></div><div>I would really appreciate some help with this.</div><div><br></div><div>Thanks,</div><div> Varsha Rani </div><div><br></div><div>PhD student</div><div>School of Physical Sciences</div><div>Jawaharlal Nehru University,</div><div>New Delhi, India</div><div><br></div></div>
</blockquote></div><br></div>