[Wien] FW: Internal Coordinates minimization for Monoclinic

Ghosh SUDDHASATTWA ssghosh at igcar.gov.in
Mon Jun 28 07:48:29 CEST 2010


Dear Prof.Blaha, 
Thanks for the mail. 
I have done the following 
1. Generate the structure file and run a normal SCF with spin orbit
coupling. (Since I am working on an heavy element with 18 atoms per unit
cell)
2. Now my question is 
We can either use option 7 and generate 15 structures and optimize the
structure with minimum energy. 
By modifying the command in optimize.job 
Run_lapw -cc 0.0001 -ec 0.0001 -p -so 

Or, we can simply do min_lapw for the structure file we have generated
initially. 
After running min_lapw, we get the atomic positions with minimum partial
forces, which gets saved in the structure file. 
Now, using the revised (optimized) atomic coordinates, run a SCF. 
By doing so, can we say that we have the structure with minimum total energy
(by doing only min_lapw and not using option 7) 

I hope I am clear enough now. 
Thank you 
Suddhasattwa 

-----Original Message-----
From: wien-bounces at zeus.theochem.tuwien.ac.at
[mailto:wien-bounces at zeus.theochem.tuwien.ac.at] On Behalf Of Peter Blaha
Sent: Monday, June 28, 2010 11:02 AM
To: A Mailing list for WIEN2k users
Subject: Re: [Wien] FW: Internal Coordinates minimization for Monoclinic

Most likely nobody can answer this question, because nobody understands your
question.
At least I did not understand what you are really doing.

Did you run    x optimize with option 7 ?
Did you modify optimize.job or just run it unmodified ?
Did you put min_lapw instead or run_lapw ?
Did you find the structure with lowest total energy ?

> Coming back to my question, doing a 4-D optimization of a monoclinic 
> case by choosing Option 7, can we say with authority that
> 
> 1. Internal coordinates minimization of monoclinic case (with 18 atoms 
> per unit cell) keeping the lattice parameters same would be enough to 
> know the total energy.
> 

> In any case, after lapw_mini, the optimized coordinates are saved in the 
> case.struct, so running the SCf again would give the optimized total 
> energy.

Why would you run again an scf ? You already have one ?


-- 
-----------------------------------------
Peter Blaha
Inst. Materials Chemistry, TU Vienna
Getreidemarkt 9, A-1060 Vienna, Austria
Tel: +43-1-5880115671
Fax: +43-1-5880115698
email: pblaha at theochem.tuwien.ac.at
-----------------------------------------
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien



More information about the Wien mailing list