[Wien] [by mistake I replies in another thread] problem in 2Doptimize.job or in my method ...
Dr. K. C. Bhamu
kcbhamu85 at gmail.com
Wed Jan 11 22:15:47 CET 2017
>
> Thank you Prof. Gerhard and Lyudmila
> Please see my additional information (some additional queries are there):
>
>
> what switches do you usually use when you run a LDA+U calculation ?
>>
> min -j "runsp_c_lapw -orb -I -fc 1.0 -i 90 "
>
>
>> why don't you use them here if you like to do a LDA+U calculation ?
>>
> Now, I corrected it.
>
>>
>> you have a NN error, when does this error occur in general ?
>>
>
> It occurs at the beginning for
> ERROR status in 2D_V-15.0_COA_-9.0 (-15 -10 -5 0 5 10 15 Vs. -9 -6 -3 0
> 3 6 9)
> ERROR status in 2D_V-10.0_COA_-9.0 (-10 -5 0 5 10 vs. -9 -6 3 0 3 6 9)
> ERROR status in 2D_V-10.0_COA_-8.0 (-10 -5 0 5 10 Vs. -8 -6 -4 -2 0 2 4 6
> 8)
> ERROR status in 2D_V-10.0_COA_-6.0 (-10 -5 0 5 10 Vs. -6 -4 -2 0 2 4 6 )
> ERROR status in 2D_V-8.0_COA_-6.0 (-8 -6 -4 -2 0 2 4 6 8 Vs. -6 -4 -2 0 2
> 4 6)
> ERROR status in 2D_V-6.0_COA_-6.0 (-6 -4 -2 0 2 4 6 Vs. -6 -4 -2 0 2 4 6)
>
> The simplest (with less parameters) case was *-6 -4 -2 0 2 4 6" for vol
> and *-6%,-3%,0,3%,6% (also -6 -4 -2 0 2 4 6 ) for c/a where I
> encountered the nn error.
>
> Now, as per suggestions of Dr. Lyudmila, I tried to see in which file the
> problem was. So:
>
> *In the new test, I did not encounter the nn problem for "-6 -4 -2 0 2 4
> 6" volume changes and -4, -2, 0 2 4 for c/a.*
> *Should I take this small value only? **I wanted to optimise it for some
> more parameters but it giving nn error beyod these range.*
>
>
>
>> how much did you reduce the rmt's and what is the nearest neighbour
>> distance for the structure with smallest volume and c/a ratio ?
>>
>
> I reduced 5% rmts and the attached one is with this reduced rmt
> (CuGaO2.struct) and nn distance is;
>
> ATOM 1 Cu ATOM 3 O
> RMT( 1)=1.75000 AND RMT( 3)=1.50000
> SUMS TO 3.25000 LT. NN-DIST= 3.43911
>
> ATOM 2 Ga ATOM 3 O
> RMT( 2)=1.85000 AND RMT( 3)=1.50000
> SUMS TO 3.35000 LT. NN-DIST= 3.73567
>
> ATOM 3 O ATOM 1 Cu
> RMT( 3)=1.50000 AND RMT( 1)=1.75000
> SUMS TO 3.25000 LT. NN-DIST= 3.43911
>
>
> For original struct, the nn distance is below (CuGaO2_orig.struct):
>
> ATOM 1 Cu ATOM 3 O
>> RMT( 1)=1.84000 AND RMT( 3)=1.58000
>> SUMS TO 3.42000 LT. NN-DIST= 3.43911
>>
>> ATOM 2 Ga ATOM 3 O
>> RMT( 2)=1.95000 AND RMT( 3)=1.58000
>> SUMS TO 3.53000 LT. NN-DIST= 3.73567
>>
>> ATOM 3 O ATOM 1 Cu
>> RMT( 3)=1.58000 AND RMT( 1)=1.84000
>> SUMS TO 3.42000 LT. NN-DIST= 3.43911
>>
>
>
>
>
>> what happens if you change the position of one of the atoms ?
>>
>
> How much I change (+/-?) and for which? would you plz suggest, I will
> apply and test it.
>
>
>>
>> is the struct file the original one or the one where the calculation
>> fails ?
>>
>
> CuGaO2.struct is with 5% rmt reduction and CuGaO2_orig.struct is original.
>
>
>>
>> If you change the script then you should know what you are doing, usually
>> it helps to read the manual to find out what the commands mean
>> a basic knowledge on linux commands is also helpful, what means cp ?
>>
>
> I am sorry. cp just copy one file to other and I understood now it from
> Prof. Lyudmila's reply. I should not change it.
>
>
>> Did you intend to run the scf cycles several times for the same structure
>> ? I guess not.
>>
>
> No, I did not. I want to do 2Doptimisation and then want to apply mBJ
> potential with the optimised one.
>
> As Ludmilla told, it doesn't make sense to arbitrarily delete commands
>> from the script by commenting them,
>> in particular, if those commands do not have anything to do with the
>> error you receive.
>>
>
> I realised it now. Thank you for correcting me.
>
>
> Sincerely
> Bhamu
> Ciao
> Gerhard
>
> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:
> "I think the problem, to be quite honest with you,
> is that you have never actually known what the question is."
>
> ====================================
> Dr. Gerhard H. Fecher
> Institut of Inorganic and Analytical Chemistry
> Johannes Gutenberg - University
> 55099 Mainz
> and
> Max Planck Institute for Chemical Physics of Solids
> 01187 Dresden
> ________________________________________
> Von: Wien [wien-bounces at zeus.theochem.tuwien.ac.at] im Auftrag von Dr. K.
> C. Bhamu [kcbhamu85 at gmail.com]
> Gesendet: Dienstag, 10. Januar 2017 18:18
> An: A Mailing list for WIEN2k users
> Betreff: [Wien] problem in 2Doptimize.job or in my method ...
>
> Dear Users
> I am running the attached struct file for 2Doptimization. The structure is
> hexagonal.
> I am running it for PBE+U.
>
> First query:
> Is it ok to use it for PBE+U ( min -j "runsp_lapw -p -i 200 -ec 0.00005
> -cc 0.0005 -fc 0.5 -I") ??
>
> when I submit the 2Doptimize.job my initialised struct file changed to
> the initial one and it gives me "nn error".
> then I tried to comment/uncomment many lines in the script (dstart and
> clmextrapol_lapw) but the problem did not solve. Finally, I uncommitted
> the line "## cp $i.struct CuGaO2.struct" see below in blue with bold. Now
> the job is running well.
>
> Second query:
> My query is whether what I did is fine of I should not change anything in
> the script. If I do not change anything in the script then the nn problem
> occurs.
>
>
> #ana2D_lapw#!/bin/csh -f
> #
> #Modify this script according to your needs
> unalias rm
> #
> # to reuse previous scf runs (without a new scf run) set answscf=y
> # and use the same "savename".
> #
> # When you make modifications (RKmax, k-mesh, XC-potentials) choose:
> # answscf=no, old_savename=savename and a new savename (eg.
> "_pbe_rk8_1000k").
> set answscf=y
> set savename=
> set old_savename=
> #
> set numbvcoa = `head -1 < numbvcoa`
> set numbcoa = $numbvcoa[2]
> #
> rm -f CuGaO2.Vconst*
> if !(-d Vconst$savename) then
> mkdir Vconst$savename
> endif
> if (-e CuGaO2.clmsum && ! -z CuGaO2.clmsum) then
> x dstart -super
> endif
> if (-e CuGaO2.clmup && ! -z CuGaO2.clmup ) then
> x dstart -super -up
> x dstart -super -dn
> endif
>
> set count = 1
> set j = 0
> @ j = $j + 1
> set jj = 1
> foreach i ( \
>
>
> 2D_V__-15.0_COA__9.0 \
> *****************
> *****************
>
> 2D_V_15.0_COA__9.0 \
> )
> echo "*******************************"
> echo $i
> echo "*******************************"
> if ( -e "${i}$savename.scf" ) then
> if ($answscf == "y" ) then
> echo "Reusing old ${i}$savename.scf file."
> else
> echo "Old ${i}$savename.scf file will not be used."
> endif
> if ($answscf == "y" ) goto exitscf
> endif
> ## cp $i.struct CuGaO2.struct changed
> from cp $i.struct CuGaO2.struct to ## cp $i.struct CuGaO2.struct
>
> # if you have a previous optimize-run:
> # cp ${i}$old_savename.struct CuGaO2.struct
> # cp ${i}$old_savename.clmsum CuGaO2.clmsum
> # cp ${i}$old_savename.clmup CuGaO2.clmup
> # cp ${i}$old_savename.clmdn CuGaO2.clmdn
> # if you want to start with dstart:
> # x dstart # -c
> # x dstart -up # -c
> # x dstart -dn # -c
> # recommended option: use charge extrapolation
> clmextrapol_lapw
> if (-e CuGaO2.clmup && \
> ! -z CuGaO2.clmup ) then
> clmextrapol_lapw -up
> clmextrapol_lapw -dn
> endif
>
> #I HIGHLY RECOMMEND TO RELAX YOUR STRUCT FILE AND THEN USE IT.
>
> #Usually low-symmetry compounds have internal coordinates
> # and you must relax them (activate -min).
>
> #run_lapw -ec 0.0001 # -p -it -min -fc 1.0 -cc 0.01
> min -j "runsp_lapw -p -i 200 -ec 0.00005 -cc 0.0005 -fc 0.5 -I"
> # runsp_lapw -ec 0.0001
>
> set stat = $status
> if ($stat) then
> echo "ERROR status in" $i
> exit 1
> endif
>
> save_lapw -f ${i}$savename
>
> exitscf:
> grepline :ENE "${i}$savename.scf" 1 > analysisENE
> grepline :VOL "${i}$savename.scf" 1 > analysisVOL
> set ene=`grep :ENE analysisENE | cut -f2 -d= `
> set vol=`grep :VOL analysisVOL | cut -f2 -d= `
> set a=`head -4 < ${i}$savename.struct | tail -1|cut -c1-9`
> set c=`head -4 < ${i}$savename.struct | tail -1|cut -c22-30`
> if ( $count > $numbcoa ) then
> @ j = $j + 1
> set numbvcoa = `head -$j < numbvcoa | tail -1`
> set numbcoa = $numbvcoa[2]
> set count = 1
> endif
> set VCOA=`echo "$i" | cut -c4-18`
> set coa=`echo " $c / $a " | bc -l`
> echo $c $a $vol $ene $coa $VCOA>>CuGaO2.Vconst$j
> @ count = $count + 1
> @ jj = $jj + 1
>
> end
> cp CuGaO2.Vconst* Vconst$savename
> echo " all CuGaO2.Vconst* files are also saved in Vconst$savename "
> echo " use ana2D_lapw to analyze the results"
> #ana2D_lapw
>
> Sincerely
> Bhamu
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> SEARCH the MAILING-LIST at: http://www.mail-archive.com/
> wien at zeus.theochem.tuwien.ac.at/index.html
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20170112/9319de88/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CuGaO2.struct
Type: application/octet-stream
Size: 3074 bytes
Desc: not available
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20170112/9319de88/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CuGaO2_orig.struct
Type: application/octet-stream
Size: 3074 bytes
Desc: not available
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20170112/9319de88/attachment-0001.obj>
More information about the Wien
mailing list