[Wien] optimization doku

Jack London yshzhang88 at yahoo.com
Mon Nov 15 19:26:05 CET 2004


Dear Peter and L. D. Marks

I download the "Optimization notes" and read it.
Thinking about my current work testing RKMAX in
vicinal surface by PORT method, I have some questions
on some files (.min_hess, .minrestart, case.tmpM and
case.finM). 

Based on your "note", I should copy .min_hess to
.minrestart then delete the old case.tmpM and
case.finM files before I change RKMAX. But in my
working script file, I didn't do that.

The following script is a part of a series script, at
first it calculat RKMAX=4.8, then ONLY clean the
directory and copy 113.in1_20(RKMAX=4.92) to 113.in1,
after that, copy 1 file to 113.in2 (to change the
"FOR" to "TOT") and "run_lapw -p -i 99 -cc 0.0001
-in1new 5" (without this commander, the :DIS can't
coverge well only using min_lapw), finally, copy 2 to
113.in2 (to change the "TOT" to "FOR") and run
"min_lapw" to get new structure and total energy. 

rm -fr 113-19.*
rm -fr 113_*
rm -fr 113.scf_mini
cp 113.in1_20 113.in1
cp 1 113.in2
run_lapw -p -i 99 -cc 0.0001 -in1new 5
clean -s
save 113-20
cp 2 113.in2
min_lapw -j "run_lapw -p -i 99 -I -in1new 3 -fc 0.5"
clean -s
save 113_in20
cp *.* /home/zhang/work/vicinal/113/hollow111/113/20/.

In my case, 
1), I didn't delete the old case.tmpM and case.finM
before changing RKMAX. 
Do you think it will affect my result? 

2), I didn't copy .min_hess to .minrestart. (min_lapw
will do it for me). But I don't know if the .min_hess
is still valid after I converged charge density at
first. 

In one word, is my script file correct to get
reasonable result? 

> http://www.wien2k.at/reg_user/textbooks/

Thank you very much!

zhang


		
__________________________________ 
Do you Yahoo!? 
Check out the new Yahoo! Front Page. 
www.yahoo.com 
 




More information about the Wien mailing list