Dear Forum,<br><br>I have some questions on setting the parameter set (ecfixed,q2sigma,qcurtz) in cp.x of espresso-4.2.1.<br>I use SiO2 of example18 in source code.<br><br>I found that when I adjust (ecfixed,q2sigma,qcurtz) the average volume at equilibrium state varies.<br>
How do think of this volume disagreement? <br>
What is the criteria to select the best (ecfixed,q2sigma,qcurtz)  value?<br><br><br><br>I summarized my results in detail following.  <br><br><br>In this code, the original value is  (ecfixed=16,q2sigma=2,qcutz=150).<br>
I did adjustment with all other parameter fixed as the same to check the pressure, volume, and total energy at equilibrium state of 10Kbar-10K. <br>I adjust ecfixed first, then q2sigma, then qcutz.<br><br>1. adjusting ecfixed.<br>
ecfixed (with q2sigma=2,qcutz=150)        average  volume              volume vibarion<br>17.1                                                                1340                     small <br>17.0                                                                1350                     small                                                                 <br>
16.9                                                                1360                     small <br>16.8                                                                1370                     small <br>16                                                                   1330                     larger  <br>
19                                                                   1390                     larger  <br>All above have similar pressure vibration.<br>The larger ecfixed, the lower total energy.<br><br>How do think of this volume discrepancy? Is this satisfying with inaccuracy?<br>
What is the criteria to select the best ecfixed value?<br>Now, without anyother base, personally, I think that ecfixed=16.8 has small volume vibration and it has the largest volume among those with small volume vibration.<br>
Maybe ecfixed=16.8 it is. <br><br>2. adjust q2sigma<br>q2sigma (with ecfixed=16.8,qcutz=150)           volume                 volume vibration<br>
1.6                                                                1420                           larger<br>
1.8                                                                1345
                           larger                                      <br>
2.0                                                                1370                           small<br>2.2                                                                1380                           larger<br>
2.4                                                                1370                           larger<br>
2.6                                                                1330                           larger<br>All above have similar pressure vibration.<br>
The smaller q2sigma, the lower total energy.<br><br>I think that q2sigma=2 is the best because it has smallest volume vibration. <br>Right?<br><br>3. adjusting qcutz<br>
qcutz (with ecfixed=16.8,q2sigma=2)           volume                 volume vibration<br>

50                                                                1310                           small<br>

80                                                                1335
                           small                                      <br>

120                                                              1360                           small<br>
150                                                              1370                           small<br>

180                                                              1380                           small<br>

400                                                              1370                           small<br>
All above have similar pressure vibration.<br>

The smaller qcutz, the lower total energy.<br><br>They have similar results except the volume discrepancy. <br>So what is the criteria to select ? <br><br><br><br>In sum, regarding the volume discrepancy, how to select the best value for  (ecfixed,q2sigma,qcurtz) ?<br>
<br><br>Thanks for you attention.<br>Any comment is grateful.<br><br><br>WANG Riping <br>2010.1.7 <br><br> <br>-- <br><div>******************************************************************************<br>WANG Riping<br>Ph.D student, </div>

<div>Institute for Study of the Earth's Interior,Okayama University,<br>827 Yamada, Misasa, Tottori-ken 682-0193, Japan<br>Tel: +81-858-43-3739(Office), 1215(Inst)</div>
<div>E-mail: <a href="mailto:wang.riping.81@gmail.com" target="_blank">wang.riping.81@gmail.com</a><br>******************************************************************************<br><br><br></div><br>