<div dir="ltr"><div><div>Thank you very much Prof. Tomas for letting me the nice tool.</div><div><br></div>I tried the OrthOpt script as usual as suggested "source volumeOrtho" in a pre-initialized case directory on a i3 Laptop with Wien2k_17.1. <br></div>It showed me several errors on screen (just copying few error messages while the optimization process is still running):<br><br>------------------------------<br> beforeCompass N says it has been already computed<br>rm: cannot remove '/home/bhamu/wien_work/palmexOrthOptA1.00/arielOrthOpt*': No such file or directory<br>position NE is A = 33.0 B = 0.5 <br>energy NE is <br><br>parent directory: <br>/home/bhamu/wien_work/palmexOrthOptA1.00<br>beforeCompass NE says it will be computed soon<br>rm: cannot remove '/home/bhamu/wien_work/palmexOrthOptA1.00/arielOrthOpt*': No such file or directory<br>position E is A = 33.0 B = 0.0 <br>energy E is <br><br>parent directory: <br>/home/bhamu/wien_work/palmexOrthOptA1.00<br>beforeCompass E says it will be computed soon<br>rm: cannot remove '/home/bhamu/wien_work/palmexOrthOptA1.00/arielOrthOpt*': No such file or directory<br>position SE is A = 33.0 B = -0.5 <br>energy SE is <br><br>parent directory: <br>/home/bhamu/wien_work/palmexOrthOptA1.00<br>beforeCompass SE says it will be computed soon<br>rm: cannot remove '/home/bhamu/wien_work/palmexOrthOptA1.00/arielOrthOpt*': No such file or directory<br>position S is A = 32.5 B = -0.5 <br>energy S is <br><br>beforeCompass S says it has been already computed<br>rm: cannot remove '/home/bhamu/wien_work/palmexOrthOptA1.00/arielOrthOpt*': No such file or directory<br><br>> stop error<br>cleaning /home/bhamu/wien_work/palmexOrthOptA-1.00/arielOrthOptSE/OrthOpt <br>position SW is A = 32.0 B = -0.5 <br>energy SW is <br><br>beforeCompass SW says it has been already computed<br>rm: cannot remove '/home/bhamu/wien_work/palmexOrthOptA1.00/arielOrthOpt*': No such file or directory<br>position W is A = 32.0 B = 0.0 <br>energy W is <br><br>beforeCompass W says it has been already computed<br>rm: cannot remove '/home/bhamu/wien_work/palmexOrthOptA1.00/arielOrthOpt*': No such file or directory<br>in /home/bhamu/wien_work/palmexOrthOptA-1.00/arielOrthOptSE/OrthOpt<br>position NW is A = 32.0 B = 0.5 <br>energy in SE : <br><br>> stop error<br><br><br><br><br><br><br>> stop error<br>cleaning /home/bhamu/wien_work/palmexOrthOptA0.00/arielOrthOptE/OrthOpt <br>in /home/bhamu/wien_work/palmexOrthOptA0.00/arielOrthOptE/OrthOpt<br>energy in E : <br><br>> stop error<br>cleaning /home/bhamu/wien_work/palmexOrthOptA0.00/arielOrthOptSE/OrthOpt <br>in /home/bhamu/wien_work/palmexOrthOptA0.00/arielOrthOptSE/OrthOpt<br>energy in SE : <br><br>> stop error<br>cleaning /home/bhamu/wien_work/palmexOrthOptA0.00/arielOrthOptNE/OrthOpt <br>in /home/bhamu/wien_work/palmexOrthOptA0.00/arielOrthOptNE/OrthOpt<br>energy in NE : <br><br><div class="gmail_extra">----------------------------------------------</div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">In addition to above, it is mentioned in "compute.job" that the script will use 15 core by default. with below statements:</div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_extra"> # If you want some more parallelization then mind that the number of <br> # processors increases as 15, 30, 45, 60 and so on <br>#cat > .machines <<+<br>#1:${HOSTNAME}:1<br>#+ # The + sign must be at the very beginning of the line</div><div class="gmail_extra"><br></div></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">So does it mean that we do not need to modify anything if we want to use only 15 cores in main cluster job file provided that the core must be greater then 15 in cluster job file? <br></div><div class="gmail_extra"><br></div><div class="gmail_extra">Suppose in my cluster job file I specify 16 cores and then the "source volumeOrtho" should use 15 cores. Is it?</div><div class="gmail_extra"><br></div><div class="gmail_extra">What about if we use 30 cores? how to change above lines in "compute.job" to use 30 cores? <br></div><br><div class="gmail_extra"><br></div><div class="gmail_extra">Kind regards</div><div class="gmail_extra">Bhamu<br clear="all"></div><div class="gmail_extra"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px"><br></span></div><div dir="ltr"><span style="font-size:12.8px"><br></span></div><br></div></div></div></div></div></div></div></div><br><div class="gmail_quote">On Wed, Oct 18, 2017 at 2:11 PM, Tomas Kana <span dir="ltr"><<a href="mailto:kana@seznam.cz" target="_blank">kana@seznam.cz</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Bhamu,<br>
There is an alternative way. You may use the<br>
OrthoOpt package available at<br>
<a href="http://susi.theochem.tuwien.ac.at/reg_user/unsupported/OrthoOpt/" rel="noreferrer" target="_blank">http://susi.theochem.tuwien.<wbr>ac.at/reg_user/unsupported/<wbr>OrthoOpt/</a><br>
It proceeds in a similar way Stefaan has suggested in previous<br>
mail copied below. It first fixes the volume and finds<br>
a minimum energy depending on two variables: c/a and b/a.<br>
Then, it moves to the next volume and so on until it reaches<br>
minimum volume.<br>
HTH<br>
Tomas<br>
<br>
No, not really. You have three (unit cell) degrees of freedom: E(a,b,c). Or similarly E(V, b/a, c/a) (with V=a*b*c). The different options in optimize_lapw give you a series of unit cells with either different values of V (and b/a and c/a constant), different values of b/a (and V and c/a constant) and different values of c/a (and V and b/a constant). These are three slices through the energy landscape. By finding the energy minimum along one slice, fix it, then go on with the next slice, etc. … you can gradually approach the absolute energy minimum. It might take several iterations through this scheme if your starting cell is far away from the minimum. Best is to start first with the degree of freedom that affects energy most (usually the volume).<br>
<br>
<br>
<br>
Stefaan<br>
<br>
Předmět: Re: [Wien] optimization of an orthorhombic cell<br>
Dear List,<br>
Could you please help me to get optimize lattice<br>
parameters from the procedure that I followed as<br>
mentioned below:<br>
The procedure is a four-step task what I am following:<br>
The original lattice parameters were: a, b, c=16.763005, 24.803600, 16.320431 vol 6785.64777697697 bohr^3.<br>
1. VARY C/A RATIO with CONSTANT VOLUME and B/A (orthorh lattice)<br>
getting below optimized data:<br>
==============================<wbr>============================<br>
Lowest point on X-axis is = -1.615800<br>
Minimum value of B/A is = 1.455755<br>
Minimum value of ENE is = -339059.104498 Ry<br>
==============================<wbr>============================<br>
Value of A is = 16.85428 bohr ; 8.91890 Ang<br>
Value of B is = 24.53569 bohr ; 12.98373 Ang<br>
Value of C is = 16.40929 bohr ; 8.68342 Ang<br>
==============================<wbr>============================<br>
Initial value of B/A is = 1.479663<br>
Conv. Unit Cell Vol = 6785.7557 bohr^3 ; **********Ang^<br>
2. Volume optimization using option 1 in "x optimize" and then calculating new volume from eplot option.<br>
Query: Step two is giving optimized volume but I do not know how to calculate lattices parameters from a single volume<br>
value for the next step "VARY C/A RATIO with CONSTANT VOLUME and B/A (orthorh lattice)"?<br>
The original cell volume and lattice parameters are mentioned at the beginning.<br>
and the new volume from step two is: 6815.3427 bohr^3<br>
The above same query is for the remaining process:<br>
3. optimization of C/A<br>
4. optimization of V.<br>
*5 At the final step I will do structural<br>
relaxation (min_.....).<br>
Could you please advice me how to proceed?<br>
Kind regards,<br>
Bhamu<br>
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
</blockquote></div><br></div></div>