Dear Amin!<br><br>Are you trying to calculate the energy of BH3 dimer? If it so then the vacuum space you have used is very less, there will be lot of spurious interaction. Increase the vacuum space such that the minimum distance between the images is greater than 15 angstrom. Please try it.<br>
<br><div class="gmail_quote">On Mon, Aug 29, 2011 at 4:04 AM, Amin Torabi <span dir="ltr"><<a href="mailto:mtorabi@uwo.ca">mtorabi@uwo.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div dir="ltr">Dear QE-users,<div class="im"><br><br>Will you please take a look at my input file and convergence table below:</div></div><div class="im"><div dir="ltr">
<br></div><div dir="ltr">I am using PBE USPP for both hydrogen and boron, but calculating energies using PBEsol. (ecutrho is set to 12* ecutwfc.)<br>
<br></div></div><div dir="ltr">My question is:<br>1. Does it look normal? I was expecting to see 1mRy accuracy at lower ecutwfc values.</div><div><div></div><div class="h5"><div dir="ltr">2. I didn't use "occupations" or "smearing" and I didn't get any error in output files. Do they affect the convergence?</div>
<div dir="ltr">3. I want to do vc-relax calculation at high pressure. Based on my table, what ecutwfc do you suggest?</div><div dir="ltr"><br></div><div dir="ltr">Thanks in advance.</div>
<div dir="ltr"><br><br><br>ecutwfc energy<br></div>30 -38.28149597<br>35 -38.28437812<br>40 -38.28554856<br>45 -38.28658930<br>
50 -38.28772083<br>55 -38.28883141<br>60 -38.28989311<br>65 -38.29086868<br>70 -38.29167889<br>75 -38.29227767<br>
80 -38.29267847<br>85 -38.29291696<br>90 -38.29302834<br>95 -38.29305836<br>100 -38.29305908<br>105 -38.29307545<br>
110 -38.29313486<br>115 -38.29324982<br><br> &control<br> pseudo_dir = **<br> outdir = **<br> prefix = '$a'<br>
tstress = .true<br> tprnfor = .true<br> verbosity = 'low'<br> /<br>&SYSTEM<br> ibrav = 12<br> celldm(1) = 8.31479435<br> celldm(2) = 1.3<br>
celldm(3) = 1.4772727273<br> celldm(4) = -0.260504508642648<br> nat = 16<br> ntyp = 2<br> ecutwfc = $a<br> ecutrho = `expr 12 \* $a`<br>
input_dft = 'pbesol'<br> /<br>&ELECTRONS<br> /<br>ATOMIC_SPECIES<br> B 10.811 B.pbe-n-van_ak.UPF<br> H 1.00794 H.pbe-van_ak.UPF<br>ATOMIC_POSITIONS crystal <br>B 0.002 0.146 0.042<br>
H -0.194 0.166 0.140<br>H 0.196 0.294 -0.005<br>H 0.104 -0.019 0.112<br>B -0.002 -0.146 -0.042<br>H 0.194 -0.166 -0.140<br>H -0.196 -0.294 0.005<br>H -0.104 0.019 -0.112<br>
B 0.498 0.354 0.542<br>H 0.694 0.334 0.640<br>H 0.304 0.206 0.495<br>H 0.396 0.519 0.612<br>B 0.502 0.646 0.458<br>H 0.306 0.666 0.360<br>H 0.696 0.794 0.505<br>
H 0.604 0.481 0.388<br>K_POINTS automatic<br> 8 8 8 0 0 0<br clear="all"><font color="#888888"><br>-- <br><div dir="ltr">Amin Torabi<br>Ph.D. Student<br>Chemistry Department<br>The University <i>of</i> Western Ontario<br>
London, On Canada, N6A 5B7<br>
<span style="font-family: 'Droid Sans',arial,sans-serif; font-size: 13px; border-collapse: collapse; color: rgb(32, 32, 32);">Phone: <a value="+15196612111" style="color: rgb(103, 117, 58);">519-661-2111</a> Ext: 87871</span><div>
<span style="font-family: 'Droid Sans',arial,sans-serif; font-size: 13px; border-collapse: collapse; color: rgb(32, 32, 32);">****************************************</span></div></div>
</font></div></div></div>
</div><div><div></div><div class="h5"><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Amin Torabi<br>Ph.D. Student<br>Chemistry Department<br>The University <i>of</i> Western Ontario<br>London, On Canada, N6A 5B7<br>
<span style="font-family: 'Droid Sans',arial,sans-serif; font-size: 13px; border-collapse: collapse; color: rgb(32, 32, 32);">Phone: <a value="+15196612111" style="color: rgb(103, 117, 58);">519-661-2111</a> Ext: 87871</span><div>
<span style="font-family: 'Droid Sans',arial,sans-serif; font-size: 13px; border-collapse: collapse; color: rgb(32, 32, 32);">****************************************</span></div></div><br>
</div></div></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Regards,<br>MOHNISH,<br>-----------------------------------------------------------------<br>Mohnish Pandey<div>BTech-Mtech, IIT Kanpur<br>Senior Project Associate,<br>
Department of Chemical Engineering,<br>IIT KANPUR, UP, INDIA<br>-----------------------------------------------------------------</div><br>