Dear Dr. Deyu Lu<br><br>I am looking for a Ti norm-conserving PSP.<br>I was wondering if I could have the generated Ti-PSP by you.<br> <br>Best Regards,<br>Mohaddeseh<br>---------------------------------------------------------<br>
Mohaddeseh Abbasnejad, <br>Room No. 323, Department of Physics, <br>University of Tehran, North Karegar Ave.,<br>Tehran, P.O. Box: 14395-547- IRAN <br>Tel. No.: +98 21 6111 8634 & Fax No.: +98 21 8800 4781<br>Cellphone: +98 917 731 7514<br>
E-Mail: <a href="mailto:m.abbasnejad@gmail.com" target="_blank">m.abbasnejad@gmail.com</a><br>Website: <a href="http://physics.ut.ac.ir/" target="_blank">http://physics.ut.ac.ir</a><br>---------------------------------------------------------<br>
<br><br><div class="gmail_quote">On Tue, Jun 28, 2011 at 5:57 AM, Deyu Lu <span dir="ltr"><<a href="mailto:deyulu@yahoo.com">deyulu@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dear PW developers and users:<br>
I have been doing norm-conserving PSP (PBE) calculations to<br>
study the structure of rutile tio2. During my study, I found the<br>
resulting lattice constant is sensitive to some fine details of<br>
the PSPs I used. Such effect is kind of counter-intuitive. Perhaps<br>
someone with a better understanding can clarify it for me. The story is a bit long, which I summarize below. Input files and some PSPs are attached.<br>
<br>
All the calculations are carried out on a Linux cluster with QE 4.2.1.<br>
The code is compiled with Intel compiler (ifort) and mvapich2-1.2,<br>
and is linked to the mkl library accompanied with the code. Main<br>
results are verified on a Cray XE6 machine (hopper at NERSC) with the same version of QE.<br>
<br>
The Ti PBE PSP was generated from an ionic configuration<br>
3s2 3p6 3d0 4s0, with semi-core included explicitly. A small rcut<br>
(0.9 0.9 1.0 for 3s 3p 3d) was used to ensure a good transferability.<br>
The psp was generated from opium, and the convergence can be<br>
reached at ecutwfn=160 Ry.<br>
<br>
The structure of Rutile TiO2 was calculated using this Ti PSP<br>
and O.pbe-mt.UPF from QE library. The lattice constant 4.587 angstrom<br>
is in excellent agreement with experiment (4.587 and 4.582 angstrom)<br>
but at odds with all the known theoretical results (ultrasoft and paw) in the literature (4.63 - 4.65). The error is about 1.1%. So our good agreement seems to be accidental. In order to obtain a better understanding, I performed the calculations using paw (vasp) and norm-conserving PSP with HGH (abinit library Theor. Chem. Acc. 114, 145 (2005)) which is known to be very accurate. Indeed, the results are consistent with literature.<br>
<br>
The 1% difference was later found to be related to the Oxygen PSP.<br>
Different Oxygen PSPs are used in the following tests.<br>
<br>
1.d: 08-O.PBE.fhi.UPF (abinit PSP library, converted to UPF, lmax=3, local=2)<br>
a=4.633 angstrom<br>
<br>
2.d: O.pbe-mt.UPF (QE PSP library, lmax=2, local=2)<br>
a=4.587 angstrom<br>
<br>
In order to see if the f channel can cause any difference, 08-O.PBE.fhi.UPF is modified (08-O.PBE.fake.UPF) by removing f channel and making corresponding changes in DIJ.<br>
No difference is found as compared to 1.d.<br>
<br>
6.d: 08-O.PBE.fake.UPF (modified version, lmax=2, local=2)<br>
a=4.633 angstrom<br>
<br>
There seems to be a bug in the <PP_HEADER> section of O.pbe-mt.UPF, where Max angular momentum component is set to 1 instead of 2. I manually changed it to 2 (O.pbe-mt.9.UPF), and no difference found in output as compared to 2.d. So this field is not sensitive in the calculation.<br>
<br>
9.d: O.pbe-mt.9.UPF (minor modification of O.pbe-mt.UPF)<br>
a=4.587 angstrom<br>
<br>
In the rest of tests, the "Max angular momentum component" is left unchanged.<br>
<br>
To find out the difference between 08-O.PBE.fake.UPF and O.pbe-mt.UPF, I plot term by term for two PSPs, and the results are shown in the xmgrace plots for mesh, local, and nonlocal parts of the PSP. The two meshes have the same size, with a max difference within 5d-3 a.u (mesh.agr). There is no difference found in local and non-local components of the PSP (local.agr and nonlocal.agr). In the inset of local.agr, we see that even the small wiggle at the origin is id<br>
entical. My conclusion is that the two PSPs are literally the same.<br>
<br>
Further tests were done by replace local PP and Beta sections in O.pbe-mt.UPF one by one by the corresponding section from 08-O.PBE.fake.UPF, and the difference is negligible as compared to 2.d.<br>
<br>
In the end, I replaced </PP_MESH> section in O.pbe-mt.UPF with that from 08-O.PBE.fake.UPF, and kept everything else unchanged, the lattice constant increased from 4.587 to 4.633.<br>
<br>
8.d: O.pbe-mt.8.UPF ( </PP_MESH> section modified)<br>
a=4.633<br>
<br>
In the end, I was quite confused about the fact that a very small change in the mesh (<5d-3 max at 80 a.u.) can cause a difference in the lattice constant by 1%.<br>
<br>
Best<br>
Deyu Lu<br>
<br>
************************************************<br>
<font color="#888888">Deyu Lu<br>
Assistant Physicist, Theory & Computation Group<br>
the Center for Functional Nanomaterials<br>
Building 735, Brookhaven National Lab<br>
Upton, NY, 11973<br>
************************************************</font><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>