<div dir="ltr">Dear Sir,<br><div><br></div><div>               In the output-order file there are 3 types of outputs. One is for order of fit 2, another is order of fit 3 and the last one is order of fit is 4. My query which output among those three should be considered?</div><div><br></div><div>with regards,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 6 Aug 2024 at 13:37, Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</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">Why runsp_lapw, isn't SnSe2 a semiconductor without magnetic order and a non-spinpolarized calculation is just fine ?<br>
<br>
According to Pearsons Crystal Database from 2022, SnSe2 should crystallize in P -3m1 (no. 164), <br>
in its crystal class -3m (also for 3m, 32) one has 6 independent elastic constants (11 12 13 14 33 44)<br>
(dependent are 22=11, 55=44, 23=13, 24=-14, 56=14,  c66 =( c11-c12)/2))<br>
<br>
If you use a hexagonal structure, you will have only 5 independent constants (11 12 13 33 44)<br>
in both cases hexagonal and trigonal (also crystal classes 3, -3 with seven independent constants) you have c66 =( c11-c12)/2 which is NOT independent<br>
<br>
As Morteza mentions, Se on 2d has a free z-parameter and thus the position needs to be optimized.<br>
Moreover, when the lattice is distorted, some atoms may occupy sites with free parameters even though the original didn't have.<br>
Therefore you need to minimize the forces and optimize all free positions that appear in all structures during the determination of the elastic constants.<br>
Usually one compares the elastic properties with fixed and relaxed positions. The differences in the equation of state p-V-E curves can be attributed to something like an intenal pressure as in real gases cpmpared to ideal ones. <br>
<br>
Did you use the correct structure and space group ?<br>
On the materials explorer project webpage you will find the matrix with 6 independent values for the stiffness or compliance tensors.<br>
<br>
Everything one needs to know is provided in textbooks !<br>
<br>
Ciao<br>
Gerhard<br>
<br>
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>
"I think the problem, to be quite honest with you,<br>
is that you have never actually known what the question is."<br>
<br>
====================================<br>
Dr. Gerhard H. Fecher<br>
Institut of Physics<br>
Johannes Gutenberg - University<br>
55099 Mainz<br>
________________________________________<br>
Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>] im Auftrag von morteza jamal via Wien [<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>]<br>
Gesendet: Montag, 5. August 2024 20:24<br>
An: <a href="mailto:shamik15041981@gmail.com" target="_blank">shamik15041981@gmail.com</a>; <a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a><br>
Cc: morteza jamal<br>
Betreff: Re: [Wien] Query regarding elastic constants of SnSe2<br>
<br>
Dear  Shamik,<br>
As i checked your scf files for C33 calculation, I saw there are forces on your atoms (atom 2) and you must relax your stuct files ( runsp_lapw  -min .....  or  min_lapw " runsp_lapw ....".<br>
<br>
Try to converge  "Energy,  Charge and force"<br>
runsp_lapw  -min ..... -ec 0.0001 -cc ... -fc ...<br>
<br>
<br>
With best,<br>
Morteza<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><span class="gmail_signature_prefix">-- </span><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>