<div dir="ltr"><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">Dear all,</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">I deal with the cell optimization and the equation of states of the trigonal structures (especially, GeO2). To get the optimized structure, I started from the experimental structure and performed a series of the relax calculations with "relax" keyword to it. During the relax calculations, first, I changed celldm(1)=a and celldm(3)=c/a parameters by keeping c parameter constant and I found the optimized celldm(1), after that, I kept celldm(1) fixed at its optimized value and I started to change only celldm(3) parameter until getting the optimized c parameter. </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">After getting the optimized celldm(1) and celldm(3), I calculated the equation of states of the structure by applying SCF calculation and finally, I used ev.x to calculate the bulk modulus with the non-cubic option, which is ~140 GPa.</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">However, the bulk modulus value I found is so different from the experimental value which is 42 GPa.  </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">When I applied the same optimization process to the cubic structures (Si-diamond and Ge-diamond), I achieved the almost same results with their experimental values. </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">I did not understand why this procedure does not work for the trigonal structures. I wonder, is there anything that I missed? </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">By the way, I also applied vc-relax calculation to the trigonal structure, again, I got the larger bulk modulus than the experimental one!</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">I would greatly appreciate if you could kindly help me with this problem.</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">Best regards,</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">Nadire</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><pre style="font-family:courier,'courier new',monospace;font-size:14px;white-space:pre-wrap;word-wrap:break-word;margin-top:0px;margin-bottom:0px;color:rgb(0,0,0);line-height:19.6px"><br></pre></div>-- <br><div class="gmail_signature"><div dir="ltr">Nadire Nayir<div>Research Assistant</div><div><div>Middle East Technical University</div><div>Physics Department</div></div></div></div>
</div>