<div dir="ltr">Dear Prof. Peter,<div><br></div><div>                    I have sent the struct file after 4 successive optimizations (a) volume, (b) c/a, (3) b/a, (4) force.  That's why the value is a bit different than the reported one in the mail. It is true that different values of c/a do not change the value of constant volume. But my question is the difference in conventional unit cell volume is coming between the cases when we are obtaining unit cell volume from volume optimization and in the next step while optimized c/a is calculated by keeping the original volume constant. </div><div><br></div><div>Looking forward to listen from you.</div><div><br></div><div>with regards,  . </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 12 Jun 2020 at 22:08, Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Just to contribute to the complicated algebra:<br>
<br>
With your struct file the volume (:VOL in scf) is:<br>
:VOL  : UNIT CELL VOLUME =    1072.33603<br>
<br>
This is different from what you wrote in your email ??? !!!!<br>
 > coming as 1073.2425 bohr^3 which get doubles when we consider ....<br>
<br>
If I do<br>
x optimize<br>
option 3    (not 2, which is for tetragonal and hexagonal lattices!!!!)<br>
and 0 and 2% change:<br>
<br>
  ca-test_coa____0.00.struct<br>
<br>
  11.264663 11.699300 16.273574 90.000000<br>
  ca-test_coa____2.00.struct<br>
<br>
  11.190551 11.622329 16.489838 90.000000<br>
--------------------------------------------------<br>
cp ca-test_coa____2.00.struct ca-test.struct<br>
x lapw0<br>
grep :VOL *scf0<br>
<br>
The volume is still the same.<br>
<br>
<br>
<br>
<br>
Am 12.06.2020 um 17:32 schrieb shamik chakrabarti:<br>
> Another thing is that from volume optimization the unit cell volume is <br>
> coming as 1073.2425 bohr^3 which get doubles when we consider the volume <br>
> of a conventional cell to become 2146.485 bohr^3  while the printed <br>
> volume of  a conventional unit cell after c/a optimization is coming as <br>
> 2058.7690 bohr^3 ; Why is this discrepancy? Is it that the conventional <br>
> cell volume about 1.9183 times the unit cell volume & not as I have <br>
> assumed. I am sending the fully optimized struct (according to my <br>
> calculation) file for your consideration.<br>
> <br>
> with regards,<br>
> <br>
> On Fri, 12 Jun 2020 at 20:46, shamik chakrabarti <br>
> <<a href="mailto:shamik15041981@gmail.com" target="_blank">shamik15041981@gmail.com</a> <mailto:<a href="mailto:shamik15041981@gmail.com" target="_blank">shamik15041981@gmail.com</a>>> wrote:<br>
> <br>
>     Dear Tomas,<br>
>                     Can you please elaborate how to know a', b' and c'<br>
>     (the final lattice parameters) from a unit cell volume optimization<br>
>     with a:b:c = constant (where a, b, c are initial lattice parameters) ?<br>
> <br>
>     with regards,<br>
> <br>
>     On Fri, 12 Jun 2020 at 17:15, Tomas Kana <<a href="mailto:kana@seznam.cz" target="_blank">kana@seznam.cz</a><br>
>     <mailto:<a href="mailto:kana@seznam.cz" target="_blank">kana@seznam.cz</a>>> wrote:<br>
> <br>
>         Dear Shamik,<br>
>         Ok, I undestand your point.<br>
>         But anyway, I think, that the problem comes<br>
>         from your equation.<br>
>         If you want to relax c/a with constant volume and<br>
>         constant b/a then I would prefer these eqations:<br>
>         original c/a = cAoriginal<br>
>         b/a = bAconstant<br>
>         new c/a = cAnew<br>
> <br>
>         original volume a*b*c<br>
>         new volume (the same as original)<br>
>         aNew * bNew * cNew  =<br>
>         aNew * (aNew * bAConstant) * (aNew * cAnew) =<br>
>         aNew^3 * bAConsant*cAnew = original volume<br>
>         You choose your desired cAnew and know0 bAconstant<br>
>         and the original volume.<br>
>         In this way, you obtain the aNew and it is easy to<br>
>         obtain bNew and cNew then. But I still do not understand the<br>
>         parameter x in your eqation.<br>
>         Maybe this could be the source of your error?<br>
>         Hope this helps<br>
>         Tomas<br>
> <br>
>         Dear Tomas,<br>
> <br>
>                           I am sorry to say that you have misunderstood<br>
>         my equation. I have used abcx^3 = a'b'c' where a', b' and c' are<br>
>         the new lattice parameters & not the old a,b,c !!! :)<br>
> <br>
> <br>
>         with regards,<br>
> <br>
> <br>
>         On Fri, 12 Jun 2020 at 16:31, Tomas Kana <<a href="mailto:kana@seznam.cz" target="_blank">kana@seznam.cz</a><br>
>         <mailto:<a href="mailto:kana@seznam.cz" target="_blank">kana@seznam.cz</a>>> wrote:<br>
> <br>
>         Dear Shamik,<br>
> <br>
> <br>
> <br>
>         If you want to keep a constant volume and<br>
> <br>
>         you use  the equation<br>
> <br>
>         abcx^3 = 332.934 = abc then only one solution is possible:<br>
> <br>
>         x = 1.  ax = a, bx = b  cx = c.<br>
> <br>
>         Your equation does not make sense to me....<br>
> <br>
>         With best regards<br>
> <br>
>         Tomas Kana<br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         where a, b, c are the initial lattice parameters. From that<br>
>         equation I have found the value of x & then finally multiply x<br>
>         with each lattice parameters to get new lattice parameters as <br>
>         a* = ax, b*=bx & c*=cx & then put<br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         ---------- Původní e-mail ----------<br>
> <br>
> <br>
> <br>
>         Od: shamik chakrabarti <<a href="mailto:shamik15041981@gmail.com" target="_blank">shamik15041981@gmail.com</a><br>
>         <mailto:<a href="mailto:shamik15041981@gmail.com" target="_blank">shamik15041981@gmail.com</a>>><br>
> <br>
> <br>
> <br>
>         Komu: A Mailing list for WIEN2k users<br>
>         <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>>><br>
> <br>
> <br>
> <br>
>         Datum: 12. 6. 2020 10:05:37<br>
> <br>
> <br>
> <br>
>         Předmět: Re: [Wien] Volume has changed during running<br>
>         optimization of c/a<br>
> <br>
>           with constant volume & b/a<br>
> <br>
> <br>
> <br>
>         No, I have created a new directory for c/a optimization. I have<br>
>         done the followings,<br>
> <br>
> <br>
> <br>
>         abcx^3 = 332.934 where a, b, c are the initial lattice<br>
>         parameters. From that equation I have found the value of x &<br>
>         then finally multiply x with each lattice parameters to get new<br>
>         lattice parameters as  a* = ax, b*=bx & c*=cx & then put<br>
> <br>
> <br>
> <br>
>         a*,<br>
> <br>
> <br>
> <br>
>         b* &  c* as new inputs for c/a optimization with constant volume<br>
>         in a new directory.<br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         with regards,<br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         On Fri, 12 Jun 2020 at 13:23, Peter Blaha<br>
>         <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>>> wrote:<br>
> <br>
> <br>
> <br>
>         How did you do:<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         "putting the lattice parameters as input for c/a optimization"   ???<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         If you continue in the previous directory, do you realize that<br>
>         there is<br>
> <br>
> <br>
> <br>
>         a   case_initial.struct   which, if present, is taken a "base<br>
>         struct file".<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         Am 12.06.2020 um 04:12 schrieb shamik chakrabarti:<br>
> <br>
> <br>
> <br>
>          > Dear wien2k users,<br>
> <br>
> <br>
> <br>
>          ><br>
> <br>
> <br>
> <br>
>          >            I am running wien2k 19.1. I have started c/a<br>
>         optimization<br>
> <br>
> <br>
> <br>
>          > with constant volume and b/a for an ABO4 compound. From volume<br>
> <br>
> <br>
> <br>
>          > optimization with a:b:c=constant I have obtained an unit cell<br>
>         volume of<br>
> <br>
> <br>
> <br>
>          > 332.934 Ang^3. While after calculation of lattice parameters<br>
>         and putting<br>
> <br>
> <br>
> <br>
>          > the lattice parameters as input for c/a optimization with<br>
>         *constant<br>
> <br>
> <br>
> <br>
>          > volume* & b/a, the printed unit cell volume is coming<br>
>         as 312.8004  Ang^3.<br>
> <br>
> <br>
> <br>
>          ><br>
> <br>
> <br>
> <br>
>          > Where I am doing wrong?<br>
> <br>
> <br>
> <br>
>          ><br>
> <br>
> <br>
> <br>
>          > Thanks in advance.<br>
> <br>
> <br>
> <br>
>          ><br>
> <br>
> <br>
> <br>
>          > with regards,<br>
> <br>
> <br>
> <br>
>          > --<br>
> <br>
> <br>
> <br>
>          > Dr. Shamik Chakrabarti<br>
> <br>
> <br>
> <br>
>          > Research Fellow<br>
> <br>
> <br>
> <br>
>          > Department of Physics<br>
> <br>
> <br>
> <br>
>          > Indian Institute of Technology Patna<br>
> <br>
> <br>
> <br>
>          > Bihta-801103<br>
> <br>
> <br>
> <br>
>          > Patna<br>
> <br>
> <br>
> <br>
>          > Bihar, India<br>
> <br>
> <br>
> <br>
>          ><br>
> <br>
> <br>
> <br>
>          > _______________________________________________<br>
> <br>
> <br>
> <br>
>          > Wien mailing list<br>
> <br>
> <br>
> <br>
>          > <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
> <br>
> <br>
> <br>
>          > <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> <br>
> <br>
> <br>
>          > SEARCH the MAILING-LIST at:<br>
>         <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
> <br>
> <br>
>          ><br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         -- <br>
> <br>
> <br>
> <br>
>         Peter Blaha<br>
> <br>
> <br>
> <br>
>         Inst.Materials Chemistry<br>
> <br>
> <br>
> <br>
>         TU Vienna<br>
> <br>
> <br>
> <br>
>         Getreidemarkt 9<br>
> <br>
> <br>
> <br>
>         A-1060 Vienna<br>
> <br>
> <br>
> <br>
>         Austria<br>
> <br>
> <br>
> <br>
>         +43-1-5880115671<br>
> <br>
> <br>
> <br>
>         _______________________________________________<br>
> <br>
> <br>
> <br>
>         Wien mailing list<br>
> <br>
> <br>
> <br>
>         <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
> <br>
> <br>
> <br>
>         <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> <br>
> <br>
> <br>
>         SEARCH the MAILING-LIST at:<br>
>         <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         -- <br>
> <br>
>         Dr. Shamik Chakrabarti<br>
> <br>
>         Research Fellow<br>
> <br>
>         Department of Physics<br>
> <br>
>         Indian Institute of Technology Patna<br>
> <br>
>         Bihta-801103<br>
> <br>
>         Patna<br>
> <br>
>         Bihar, India<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         _______________________________________________<br>
> <br>
> <br>
> <br>
>         Wien mailing list<br>
> <br>
> <br>
> <br>
>         <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
> <br>
> <br>
> <br>
>         <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> <br>
> <br>
> <br>
>         SEARCH the MAILING-LIST at:<br>
>         <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
> <br>
> <br>
>         _______________________________________________<br>
> <br>
>         Wien mailing list<br>
> <br>
>         <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
> <br>
>         <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> <br>
>         SEARCH the MAILING-LIST at:<br>
>         <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         -- <br>
>         Dr. Shamik Chakrabarti<br>
>         Research Fellow<br>
>         Department of Physics<br>
>         Indian Institute of Technology Patna<br>
>         Bihta-801103<br>
>         Patna<br>
>         Bihar, India<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
>         _______________________________________________<br>
> <br>
>         Wien mailing list<br>
> <br>
>         <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
> <br>
>         <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> <br>
>         SEARCH the MAILING-LIST at:<br>
>         <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
>         _______________________________________________<br>
>         Wien mailing list<br>
>         <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>         <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
>         <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
>         SEARCH the MAILING-LIST at:<br>
>         <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
> <br>
> <br>
>     -- <br>
>     Dr. Shamik Chakrabarti<br>
>     Research Fellow<br>
>     Department of Physics<br>
>     Indian Institute of Technology Patna<br>
>     Bihta-801103<br>
>     Patna<br>
>     Bihar, India<br>
> <br>
> <br>
> <br>
> -- <br>
> Dr. Shamik Chakrabarti<br>
> Research Fellow<br>
> Department of Physics<br>
> Indian Institute of Technology Patna<br>
> Bihta-801103<br>
> Patna<br>
> Bihar, India<br>
> <br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.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/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
<br>
-- <br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a>    WIEN2k: <a href="http://www.wien2k.at" rel="noreferrer" target="_blank">http://www.wien2k.at</a><br>
WWW: <br>
<a href="http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------" rel="noreferrer" target="_blank">http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------</a> <br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.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/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:small">Dr. Shamik Chakrabarti</div><div style="font-size:small">Research Fellow </div><div style="font-size:small">Department of Physics</div><div style="font-size:small">Indian Institute of Technology Patna</div><div style="font-size:small">Bihta-801103</div><div style="font-size:small">Patna</div><div style="font-size:small">Bihar, India</div></div></div></div></div></div></div></div>