[Wien] Error in BoltzTraP with non-centrosymmetric materials
I was calculating many compounds with non-centrosymmetric structure but did not encounter problems. Did you use a very large number of k-points ? the number of points generated by kgen is, indded, higher for non-centrosymmetric structures May be there was some restriction, but I don't remember. Can you check with a more simple non-centrosymmetric example (something with Zincblende structure) and if it works increase the number of k-points in steps (rerun only lapw1c to have the energy files) Last not least: all inputs are correct ? Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag von quot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Sonntag, 12. Februar 2012 07:33 An: A Mailing list for WIEN2k users Betreff: [Wien] Error in BoltzTraP with non-centrosymmetric materials Sorry for having a mistake with the title of the latest mail, Recently I installed BoltzTraP to calculate the transport values of CaIrSi3, but I encountered the following message at the end of case.outputtrans and cannot get the transport properties: == End WIEN interface === Input file read successfully NON-CENTROSYMMETRIC. ADDING i. The case.outputtrans of CoSb3 (an example in tests/) is == End WIEN interface === Input file read successfully Finally! Starting Boltzmann calculation! Calling DOS. I don't know why the program does not work in the case of non-centrosymmetric compounds, but if someone knows how to deal with it, could you tell me? Thanks in advance, Gaku Eguchi -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
[Wien] bug for Cs LO ?
OK, I got it. Thanks! Stefaan On 12/02/2012 21:57, Peter Blaha wrote: One has to make a decision whether a certain state should be treated by an extra LO or as normal valence and the criterium is its energy in an atom. Cs-p states are just on the border. You are absolutely correct: for Cs metal one should add a LO. However, for ionic Cs compounds the 5p semicore state comes even closer to EF (just 5 eV below EF) and putting an extra LO might be dangerous and lead to ghostbands. I do not see an easy optimal solution for this problem and one has to live with some handwork. It seems my question was formulated in a misleading way. It was not about what -in1ef does or does not do, but rather about why the default case.in1 for bulk Cs looks the way it does. I'm only rarely disappointed by the case.in1 which is proposed by default -- bulk Cs is the exception. When inspecting case.scf, it looks pretty obvious that there should be a p-LO (fully occupied p-orbitals only 0.8 Ry below E_Fermi). The small QTL-B warning is consistent with this. And the resulting E(V) curve without p-LO is really noisy. Yet the default case.in1 has only the valence p. That does not look right, or at least: not optimal. Stefaan Am 11.02.2012 11:19, schrieb Stefaan Cottenier: When running a straightforward calculation for bcc-Cs (case.struct pasted underneath) with the -in1ef switch, the following case.in1 is generated: WFFIL EF= 0.5 (WFFIL, WFPRI, ENFIL, SUPWF) 7.50 10 4 (R-MT*K-MAX; MAX L IN WF, V-NMT -.0792 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW) 2 -.0792 0.000 CONT 1 2 -5.08 0.001 STOP 1 0 -1.66 0.002 CONT 1 0 -.0792 0.000 CONT 1 1 -.0792 0.000 CONT 1 K-VECTORS FROM UNIT:4 -9.0 4.5 21 red emin/emax/nband There is only one linearization energy given for p-states (valence 6p). The semi-core 5p is missing, although it is correctly given in case.scf: :BAN7: 7 -0.672053 -0.662419 1. :BAN8: 8 -0.669460 -0.662395 1. :BAN9: 9 -0.669324 -0.662351 1. This leads to a small QTL-B value, which triggers a warning only, yet which leads to intolerable scatter in the E(V) curve: QTL-B VALUE .EQ. 2.38282 in Band of energy -0.66932 ATOM= 1 L= 1 Adding the missing semi-core linearization energy manually does solve the problem. However, the fact that it is not put there automatically probably is due to a small bug? Stefaan bcc-Cs B LATTICE,NONEQUIV.ATOMS: 1 229 Im-3m MODE OF CALC=RELA unit=bohr 11.602049 11.602049 11.602049 90.00 90.00 90.00 ATOM 1: X=0. Y=0. Z=0. MULT= 1 ISPLIT= 2 Cs1 NPT= 781 R0=.1 RMT= 2.5 Z: 55. LOCAL ROT MATRIX: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 0 NUMBER OF SYMMETRY OPERATIONS ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
[Wien] bug for Cs LO ?
By the way, Cs 5p exhibits a nice spin-orbit splitting in all experiments Ciao Gerhard Dr. Gerhard H. Fecher Max Planck Institute for Chemical Physics of Solids 01187 Dresden Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag von quot;Stefaan Cottenier [Stefaan.Cottenier at UGent.be] Gesendet: Montag, 13. Februar 2012 08:59 An: A Mailing list for WIEN2k users Betreff: Re: [Wien] bug for Cs LO ? OK, I got it. Thanks! Stefaan On 12/02/2012 21:57, Peter Blaha wrote: One has to make a decision whether a certain state should be treated by an extra LO or as normal valence and the criterium is its energy in an atom. Cs-p states are just on the border. You are absolutely correct: for Cs metal one should add a LO. However, for ionic Cs compounds the 5p semicore state comes even closer to EF (just 5 eV below EF) and putting an extra LO might be dangerous and lead to ghostbands. I do not see an easy optimal solution for this problem and one has to live with some handwork. It seems my question was formulated in a misleading way. It was not about what -in1ef does or does not do, but rather about why the default case.in1 for bulk Cs looks the way it does. I'm only rarely disappointed by the case.in1 which is proposed by default -- bulk Cs is the exception. When inspecting case.scf, it looks pretty obvious that there should be a p-LO (fully occupied p-orbitals only 0.8 Ry below E_Fermi). The small QTL-B warning is consistent with this. And the resulting E(V) curve without p-LO is really noisy. Yet the default case.in1 has only the valence p. That does not look right, or at least: not optimal. Stefaan Am 11.02.2012 11:19, schrieb Stefaan Cottenier: When running a straightforward calculation for bcc-Cs (case.struct pasted underneath) with the -in1ef switch, the following case.in1 is generated: WFFIL EF= 0.5 (WFFIL, WFPRI, ENFIL, SUPWF) 7.50 10 4 (R-MT*K-MAX; MAX L IN WF, V-NMT -.0792 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW) 2 -.0792 0.000 CONT 1 2 -5.08 0.001 STOP 1 0 -1.66 0.002 CONT 1 0 -.0792 0.000 CONT 1 1 -.0792 0.000 CONT 1 K-VECTORS FROM UNIT:4 -9.0 4.5 21 red emin/emax/nband There is only one linearization energy given for p-states (valence 6p). The semi-core 5p is missing, although it is correctly given in case.scf: :BAN7: 7 -0.672053 -0.662419 1. :BAN8: 8 -0.669460 -0.662395 1. :BAN9: 9 -0.669324 -0.662351 1. This leads to a small QTL-B value, which triggers a warning only, yet which leads to intolerable scatter in the E(V) curve: QTL-B VALUE .EQ. 2.38282 in Band of energy -0.66932 ATOM= 1 L= 1 Adding the missing semi-core linearization energy manually does solve the problem. However, the fact that it is not put there automatically probably is due to a small bug? Stefaan bcc-Cs B LATTICE,NONEQUIV.ATOMS: 1 229 Im-3m MODE OF CALC=RELA unit=bohr 11.602049 11.602049 11.602049 90.00 90.00 90.00 ATOM 1: X=0. Y=0. Z=0. MULT= 1 ISPLIT= 2 Cs1 NPT= 781 R0=.1 RMT= 2.5 Z: 55. LOCAL ROT MATRIX: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 0 NUMBER OF SYMMETRY OPERATIONS ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
[Wien] Error in BoltzTraP with non-centrosymmetric materials
Dear Fecher Gerhard, I appreciate your reply, and it's really helpful for me. I understood the error would not come from the inversion symmetry lacking. I've used k=4000 and I think it's not so large. Along with a manual, I put the .energyso and .struct file into a folder with .intrans file, then executed BoltzTraP. The .intrans file was copied from tests/Bi2Te3 and only the Fermi level was changed. Several files were generated, but the process seemed to stop with the message `NON-CENTROSYMMETRIC. ADDING i'. The .engre file was empty and I felt it's strange. The test run with examples were successful, so the installation would be fine. I might have to change the .intrans file which I didn't notice. I'm so sorry to bother you, but if you have any comments or know howto, could you tell me? Thanks in advance, Gaku Eguchi (12/02/13 15:19), Fecher, Gerhard wrote: I was calculating many compounds with non-centrosymmetric structure but did not encounter problems. Did you use a very large number of k-points ? the number of points generated by kgen is, indded, higher for non-centrosymmetric structures May be there was some restriction, but I don't remember. Can you check with a more simple non-centrosymmetric example (something with Zincblende structure) and if it works increase the number of k-points in steps (rerun only lapw1c to have the energy files) Last not least: all inputs are correct ? Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag vonquot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Sonntag, 12. Februar 2012 07:33 An: A Mailing list for WIEN2k users Betreff: [Wien] Error in BoltzTraP with non-centrosymmetric materials Sorry for having a mistake with the title of the latest mail, Recently I installed BoltzTraP to calculate the transport values of CaIrSi3, but I encountered the following message at the end of case.outputtrans and cannot get the transport properties: == End WIEN interface === Input file read successfully NON-CENTROSYMMETRIC. ADDING i. The case.outputtrans of CoSb3 (an example in tests/) is == End WIEN interface === Input file read successfully Finally! Starting Boltzmann calculation! Calling DOS. I don't know why the program does not work in the case of non-centrosymmetric compounds, but if someone knows how to deal with it, could you tell me? Thanks in advance, Gaku Eguchi -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en
[Wien] Error in BoltzTraP with non-centrosymmetric materials
you have to check what it is doing after the message `NON-CENTROSYMMETRIC. ADDING i'. this is a regular message that it has found a noncentro structure (note: Wien gives itself a simmilar message somewhere during initialisation) did you calculate with SO ? in case, try first without and use the case.energy Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag von quot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Montag, 13. Februar 2012 13:02 An: wien at zeus.theochem.tuwien.ac.at Betreff: Re: [Wien] Error in BoltzTraP with non-centrosymmetric materials Dear Fecher Gerhard, I appreciate your reply, and it's really helpful for me. I understood the error would not come from the inversion symmetry lacking. I've used k=4000 and I think it's not so large. Along with a manual, I put the .energyso and .struct file into a folder with .intrans file, then executed BoltzTraP. The .intrans file was copied from tests/Bi2Te3 and only the Fermi level was changed. Several files were generated, but the process seemed to stop with the message `NON-CENTROSYMMETRIC. ADDING i'. The .engre file was empty and I felt it's strange. The test run with examples were successful, so the installation would be fine. I might have to change the .intrans file which I didn't notice. I'm so sorry to bother you, but if you have any comments or know howto, could you tell me? Thanks in advance, Gaku Eguchi (12/02/13 15:19), Fecher, Gerhard wrote: I was calculating many compounds with non-centrosymmetric structure but did not encounter problems. Did you use a very large number of k-points ? the number of points generated by kgen is, indded, higher for non-centrosymmetric structures May be there was some restriction, but I don't remember. Can you check with a more simple non-centrosymmetric example (something with Zincblende structure) and if it works increase the number of k-points in steps (rerun only lapw1c to have the energy files) Last not least: all inputs are correct ? Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag vonquot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Sonntag, 12. Februar 2012 07:33 An: A Mailing list for WIEN2k users Betreff: [Wien] Error in BoltzTraP with non-centrosymmetric materials Sorry for having a mistake with the title of the latest mail, Recently I installed BoltzTraP to calculate the transport values of CaIrSi3, but I encountered the following message at the end of case.outputtrans and cannot get the transport properties: == End WIEN interface === Input file read successfully NON-CENTROSYMMETRIC. ADDING i. The case.outputtrans of CoSb3 (an example in tests/) is == End WIEN interface === Input file read successfully Finally! Starting Boltzmann calculation! Calling DOS. I don't know why the program does not work in the case of non-centrosymmetric compounds, but if someone knows how to deal with it, could you tell me? Thanks in advance, Gaku Eguchi -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
[Wien] Error in BoltzTraP with non-centrosymmetric materials
Did you switch to CALC in intrans or did you stay with NOCALC ? I guess its clear what the latter means. Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag von quot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Montag, 13. Februar 2012 13:02 An: wien at zeus.theochem.tuwien.ac.at Betreff: Re: [Wien] Error in BoltzTraP with non-centrosymmetric materials Dear Fecher Gerhard, I appreciate your reply, and it's really helpful for me. I understood the error would not come from the inversion symmetry lacking. I've used k=4000 and I think it's not so large. Along with a manual, I put the .energyso and .struct file into a folder with .intrans file, then executed BoltzTraP. The .intrans file was copied from tests/Bi2Te3 and only the Fermi level was changed. Several files were generated, but the process seemed to stop with the message `NON-CENTROSYMMETRIC. ADDING i'. The .engre file was empty and I felt it's strange. The test run with examples were successful, so the installation would be fine. I might have to change the .intrans file which I didn't notice. I'm so sorry to bother you, but if you have any comments or know howto, could you tell me? Thanks in advance, Gaku Eguchi (12/02/13 15:19), Fecher, Gerhard wrote: I was calculating many compounds with non-centrosymmetric structure but did not encounter problems. Did you use a very large number of k-points ? the number of points generated by kgen is, indded, higher for non-centrosymmetric structures May be there was some restriction, but I don't remember. Can you check with a more simple non-centrosymmetric example (something with Zincblende structure) and if it works increase the number of k-points in steps (rerun only lapw1c to have the energy files) Last not least: all inputs are correct ? Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag vonquot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Sonntag, 12. Februar 2012 07:33 An: A Mailing list for WIEN2k users Betreff: [Wien] Error in BoltzTraP with non-centrosymmetric materials Sorry for having a mistake with the title of the latest mail, Recently I installed BoltzTraP to calculate the transport values of CaIrSi3, but I encountered the following message at the end of case.outputtrans and cannot get the transport properties: == End WIEN interface === Input file read successfully NON-CENTROSYMMETRIC. ADDING i. The case.outputtrans of CoSb3 (an example in tests/) is == End WIEN interface === Input file read successfully Finally! Starting Boltzmann calculation! Calling DOS. I don't know why the program does not work in the case of non-centrosymmetric compounds, but if someone knows how to deal with it, could you tell me? Thanks in advance, Gaku Eguchi -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
[Wien] Error in BoltzTraP with non-centrosymmetric materials
/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en -- next part -- An HTML attachment was scrubbed... URL: http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20120213/2a12dfa5/attachment.htm
[Wien] Error in BoltzTraP with non-centrosymmetric materials
Dear Dr. Gerhard. Did you switch to CALC in intrans or did you stay with NOCALC ? This was the cause! It worked after changing NOCALC -- CALC. Thank you very very much. Best, Gaku (12/02/13 21:29), Fecher, Gerhard wrote: Did you switch to CALC in intrans or did you stay with NOCALC ? I guess its clear what the latter means. Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag vonquot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Montag, 13. Februar 2012 13:02 An: wien at zeus.theochem.tuwien.ac.at Betreff: Re: [Wien] Error in BoltzTraP with non-centrosymmetric materials Dear Fecher Gerhard, I appreciate your reply, and it's really helpful for me. I understood the error would not come from the inversion symmetry lacking. I've used k=4000 and I think it's not so large. Along with a manual, I put the .energyso and .struct file into a folder with .intrans file, then executed BoltzTraP. The .intrans file was copied from tests/Bi2Te3 and only the Fermi level was changed. Several files were generated, but the process seemed to stop with the message `NON-CENTROSYMMETRIC. ADDING i'. The .engre file was empty and I felt it's strange. The test run with examples were successful, so the installation would be fine. I might have to change the .intrans file which I didn't notice. I'm so sorry to bother you, but if you have any comments or know howto, could you tell me? Thanks in advance, Gaku Eguchi (12/02/13 15:19), Fecher, Gerhard wrote: I was calculating many compounds with non-centrosymmetric structure but did not encounter problems. Did you use a very large number of k-points ? the number of points generated by kgen is, indded, higher for non-centrosymmetric structures May be there was some restriction, but I don't remember. Can you check with a more simple non-centrosymmetric example (something with Zincblende structure) and if it works increase the number of k-points in steps (rerun only lapw1c to have the energy files) Last not least: all inputs are correct ? Ciao Gerhard Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at zeus.theochem.tuwien.ac.at]quot; im Auftrag vonquot;EGUCHI Gaku [geguchi at scphys.kyoto-u.ac.jp] Gesendet: Sonntag, 12. Februar 2012 07:33 An: A Mailing list for WIEN2k users Betreff: [Wien] Error in BoltzTraP with non-centrosymmetric materials Sorry for having a mistake with the title of the latest mail, Recently I installed BoltzTraP to calculate the transport values of CaIrSi3, but I encountered the following message at the end of case.outputtrans and cannot get the transport properties: == End WIEN interface === Input file read successfully NON-CENTROSYMMETRIC. ADDING i. The case.outputtrans of CoSb3 (an example in tests/) is == End WIEN interface === Input file read successfully Finally! Starting Boltzmann calculation! Calling DOS. I don't know why the program does not work in the case of non-centrosymmetric compounds, but if someone knows how to deal with it, could you tell me? Thanks in advance, Gaku Eguchi -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Gaku Eguchi Department of Physics, Graduate School of Science, Kyoto University Oiwake-cho, Kitashirakawa, Sakyou-ku, Kyoto 606-8502, Japan Laboratory TEL : +81-75-753-3744 E-mail : geguchi at scphys.kyoto-u.ac.jp http://www.ss.scphys.kyoto-u.ac.jp/index.html.en ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
[Wien] lapwso crash
Nilton, That's surprising. Usually, Update 3 works and Update 8 doesn't. I can only guess that it is because you are using the static mkl libraries, whereas the dynamic libraries were likely used in other cases. Anyway, glad to hear that Update 8 is working well for you. Gavin On 2/13/2012 7:12 AM, Nilton wrote: Dear Gavin, actually, the Update 3 don't work in my case, but Update sp1.8.273 works very well. Nilton 2012/2/10 Nilton nilton.dantas at gmail.com mailto:nilton.dantas at gmail.com Dear Gavin, thanks a lot for the clues. I will to try. Nilton 2012/2/10 Gavin Abo gsabo at crimson.ua.edu mailto:gsabo at crimson.ua.edu Dear Nilton, The following steps might be useful: (1) You reported that you are using Update 2 of composer xe 2011. I would suggested trying Update 3, since there seems to be a lot of success with this compiler version. If the error goes away with your compiler settings, then you will know that the cause is coming from the compiler. It seems users are running into SIGSEGV errors during runtime with composer xe 2011 versions. Furthermore, this problem seems particular to systems with Intel Xeon processors (and maybe only those that are x64). Falling back on Update 3 or an ifort 11 usually provided a solution. (2) If you add -g your compiler options, it might give more details on the cause of the SIGSEGV error. (3) The links below list some causes and troubleshooting tips for SIGSEGV errors. http://www.democritos.it/pipermail/pw_forum/2005-March/002272.html http://software.intel.com/en-us/articles/determining-root-cause-of-sigsegv-or-sigbus-errors/ Hope this helps, Gavin On 2/9/2012 2:36 PM, Nilton wrote: Dear Peter, I dont think so. I have used case.inso generated by initso_lapw. See below his content. WFFIL 4 0 0 llmax,ipr,kpot -10 1.5Emin, Emax 0 0 1 h,k,l (direction of magnetization) 0 number of atoms with RLO 0 0 number of atoms without SO, atomnumbers -- Nilton S. Dantas Universidade Estadual de Feira de Santana Departamento de Ci?ncias Exatas ?rea de Inform?tica Av. Transnordestina, S/N, Bairro Novo Horizonte CEP 44036900 - Feira de Santana, Bahia, Brasil Tel./Fax +55 75 31618086 tel:%2B55%2075%2031618086 http://www2.ecomp.uefs.br/ http://www.uefs.br/portal ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at mailto:Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at mailto:Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Nilton S. Dantas Universidade Estadual de Feira de Santana Departamento de Ci?ncias Exatas ?rea de Inform?tica Av. Transnordestina, S/N, Bairro Novo Horizonte CEP 44036900 - Feira de Santana, Bahia, Brasil Tel./Fax +55 75 31618086 tel:%2B55%2075%2031618086 http://www2.ecomp.uefs.br/ http://www.uefs.br/portal -- Nilton S. Dantas Universidade Estadual de Feira de Santana Departamento de Ci?ncias Exatas ?rea de Inform?tica Av. Transnordestina, S/N, Bairro Novo Horizonte CEP 44036900 - Feira de Santana, Bahia, Brasil Tel./Fax +55 75 31618086 http://www2.ecomp.uefs.br/ http://www.uefs.br/portal ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- next part -- An HTML attachment was scrubbed... URL: http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20120213/a28de0a2/attachment.htm
[Wien] Time taken to run a 60 atom/cell program
Dear all, I am new to this Wien2k program. So would you please answer me? Will it take only hours or 1-2 days to run? I know the question might be too trivial, but hopefully those who have experiences can help me out. Thanks, chee -- next part -- An HTML attachment was scrubbed... URL: http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20120213/5440530f/attachment.htm
