[Wien] lapw2 aborted during SCF using 3rd structure in force minimization
I was running force minimization of a 56 atomic cell ( an oxide material) by constraining the coordinates of elements other than oxygen as that is the only variable parameter (other coordinates are fixed obeying space group symmetry) in wien2k 11.1 I was using NEW1 minimization method following user guide -- In case of constraint calculations one should use NEW1 You did not read the UG correctly: a) When positions are fixed by symmetry, this is NOT a constrain at all for the minimizers. b) The UG says: when using linear constrains, use NEW1. Linear constrains are specified in case.constrain and restrict something like x+y=0.5 or similar. About your error: you have not found anything, so how should we find anything ? Two structure files having different atomic coordinates were already created. while when running SCF calculation using the 3rd struct file, SCF was aborted by mentioning the line at the 4th iteration: lapw2 -c -up (15:01:51) Abort 3.348u 0.320s 0:01.61 227.3%0+0k 0+20376io 0pf+0w error: command /home/shamik81/WIEN2k/lapw2c uplapw2.def failed stop error Another unusual thing is that the energy and charge convergence in different iterations of this SCF cycle went on as below: EC: .41 CC: .0007728 EC: .00600493 CC. .0028021 EC: .002993535 CC. .0028017 There was no error file created and I haven't found any error indicated in case.output2 (up/dn) file and case.scf2up file was empty as lapw2 got aborted when it was just started. Charge distance in the above SCF cycle went on as following: :DIS : CHARGE DISTANCE ( 0.0007728 for atom 48 spin 2) 0.0002381 :DIS : CHARGE DISTANCE ( 0.0028021 for atom 29 spin 2) 0.0019903 :DIS : CHARGE DISTANCE ( 0.0028017 for atom 29 spin 2) 0.0019911 We are not able to figure out what the problem is. Any response in this regard will be very helpful for us. with regards, -- Shamik Chakrabarti Senior Research Fellow Dept. of Physics Meteorology Material Processing Solid State Ionics Lab IIT Kharagpur Kharagpur 721302 INDIA ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- P.Blaha -- Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna Phone: +43-1-58801-165300 FAX: +43-1-58801-165982 Email: blaha at theochem.tuwien.ac.atWWW: http://info.tuwien.ac.at/theochem/ --
[Wien] lapw2 aborted during SCF using 3rd structure in force minimization
Dear Prof. Peter Blaha, Thank you very much for your reply and clarification. I will rerun the calculation again without using any constraints. Thanks once again. with best regards, On Fri, Aug 24, 2012 at 12:59 PM, Peter Blaha pblaha at theochem.tuwien.ac.atwrote: I was running force minimization of a 56 atomic cell ( an oxide material) by constraining the coordinates of elements other than oxygen as that is the only variable parameter (other coordinates are fixed obeying space group symmetry) in wien2k 11.1 I was using NEW1 minimization method following user guide -- In case of constraint calculations one should use NEW1 You did not read the UG correctly: a) When positions are fixed by symmetry, this is NOT a constrain at all for the minimizers. b) The UG says: when using linear constrains, use NEW1. Linear constrains are specified in case.constrain and restrict something like x+y=0.5 or similar. About your error: you have not found anything, so how should we find anything ? Two structure files having different atomic coordinates were already created. while when running SCF calculation using the 3rd struct file, SCF was aborted by mentioning the line at the 4th iteration: lapw2 -c -up (15:01:51) Abort 3.348u 0.320s 0:01.61 227.3%0+0k 0+20376io 0pf+0w error: command /home/shamik81/WIEN2k/lapw2c uplapw2.def failed stop error Another unusual thing is that the energy and charge convergence in different iterations of this SCF cycle went on as below: EC: .41 CC: .0007728 EC: .00600493 CC. .0028021 EC: .002993535 CC. .0028017 There was no error file created and I haven't found any error indicated in case.output2 (up/dn) file and case.scf2up file was empty as lapw2 got aborted when it was just started. Charge distance in the above SCF cycle went on as following: :DIS : CHARGE DISTANCE ( 0.0007728 for atom 48 spin 2) 0.0002381 :DIS : CHARGE DISTANCE ( 0.0028021 for atom 29 spin 2) 0.0019903 :DIS : CHARGE DISTANCE ( 0.0028017 for atom 29 spin 2) 0.0019911 We are not able to figure out what the problem is. Any response in this regard will be very helpful for us. with regards, -- Shamik Chakrabarti Senior Research Fellow Dept. of Physics Meteorology Material Processing Solid State Ionics Lab IIT Kharagpur Kharagpur 721302 INDIA __**_ Wien mailing list Wien at zeus.theochem.tuwien.ac.**at Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.**ac.at/mailman/listinfo/wienhttp://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- P.Blaha --**--** -- Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna Phone: +43-1-58801-165300 FAX: +43-1-58801-165982 Email: blaha at theochem.tuwien.ac.atWWW: http://info.tuwien.ac.at/** theochem/ http://info.tuwien.ac.at/theochem/ --**--** -- __**_ Wien mailing list Wien at zeus.theochem.tuwien.ac.**at Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.**ac.at/mailman/listinfo/wienhttp://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Shamik Chakrabarti Senior Research Fellow Dept. of Physics Meteorology Material Processing Solid State Ionics Lab IIT Kharagpur Kharagpur 721302 INDIA -- next part -- An HTML attachment was scrubbed... URL: http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20120824/a99db4ab/attachment.htm
[Wien] lapw2 aborted during SCF using 3rd structure in force minimization
Dear Prof. Laurence Marks, I have constrained some fixed positions for keeping them intact due to the following reason: (i) I am doing a calculation which is a primitive supercell generated from the original unit cell having a particular space group. (ii) Similar type of calculation on a different material (by generating primitive space group) showed earlier that if the original structure has a coordinate say (0.5, 0, 0) after force minimization with the primitive space group the optimized coordinates eventually turned out to be (0.4995, 0.9998, 0.9997) and hence we have again came back to our original position (0.5, 0, 0). Hence, optimizing these fixed coordinates seems to be an wastage of computational time while only varying the variable coordinates (of the primitive space group) may give the equilibrium structure in a relatively small time. As you have mentioned, that I have not give enough information for finding out the reason behind the crash, Sir, can you please suggest me which information should I give?as already I have given the following information: (i) There is no error file (ii) case.scf2up is empty (as lapw2 -c -up crashed) (iii) I am not able to figure out any error messages from case.output2 Also, the initial positions for the structure has been entered carefully. The coordinates are determined from Rietveld analysis of experimental XRD which is also confronted with earlier findings. Looking forward to your response. with best regards, On Thu, Aug 23, 2012 at 6:19 PM, Laurence Marks L-marks at northwestern.eduwrote: You have misunderstood the UG. Wien2k knows about fixed positions due to symmetry and so long as you have entered your positions with precision, e.g. 0. not 0. you should not impose additional constraints yourself. In addition there is no reason to use NEW1 if you do have constraints as both PORT and MSR1a obey them. However, unless you have a very good reason I do not recommend using constraints. As to why your calculation crashed, no idea as you have not provided enough information for anyone to guess. I suspect that you have chosen inappropriate parameters and/or have a mistake in your input. Most oxides are centrosymmetric and a common mistake is to get a lower symmetry by not entering the initial positions with enough accuracy. On Thu, Aug 23, 2012 at 5:39 AM, shamik chakrabarti shamikphy at gmail.com wrote: Dear wien2k users, I was running force minimization of a 56 atomic cell ( an oxide material) by constraining the coordinates of elements other than oxygen as that is the only variable parameter (other coordinates are fixed obeying space group symmetry) in wien2k 11.1 I was using NEW1 minimization method following user guide -- In case of constraint calculations one should use NEW1 Two structure files having different atomic coordinates were already created. while when running SCF calculation using the 3rd struct file, SCF was aborted by mentioning the line at the 4th iteration: lapw2 -c -up (15:01:51) Abort 3.348u 0.320s 0:01.61 227.3% 0+0k 0+20376io 0pf+0w error: command /home/shamik81/WIEN2k/lapw2c uplapw2.def failed stop error Another unusual thing is that the energy and charge convergence in different iterations of this SCF cycle went on as below: EC: .41 CC: .0007728 EC: .00600493 CC. .0028021 EC: .002993535 CC. .0028017 There was no error file created and I haven't found any error indicated in case.output2 (up/dn) file and case.scf2up file was empty as lapw2 got aborted when it was just started. Charge distance in the above SCF cycle went on as following: :DIS : CHARGE DISTANCE ( 0.0007728 for atom 48 spin 2) 0.0002381 :DIS : CHARGE DISTANCE ( 0.0028021 for atom 29 spin 2) 0.0019903 :DIS : CHARGE DISTANCE ( 0.0028017 for atom 29 spin 2) 0.0019911 We are not able to figure out what the problem is. Any response in this regard will be very helpful for us. with regards, -- Shamik Chakrabarti Senior Research Fellow Dept. of Physics Meteorology Material Processing Solid State Ionics Lab IIT Kharagpur Kharagpur 721302 INDIA -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu 1-847-491-3996 Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Shamik Chakrabarti Senior Research Fellow Dept. of Physics Meteorology Material Processing Solid State Ionics Lab IIT Kharagpur Kharagpur 721302 INDIA -- next part -- An HTML attachment was scrubbed... URL: