Dear Stefan Cottenier Sir,<div><br></div><div>                                      Thank you very much for your reply. Yes you are right about your guess and Sir your answer really pulled us out from lots of worries. Thank you very much Sir.</div>
<div><br></div><div>with best regards,</div><div>Shamik Chakrabarti<br><br><div class="gmail_quote">On Tue, Oct 26, 2010 at 5:21 PM, Stefaan Cottenier <span dir="ltr">&lt;<a href="mailto:Stefaan.Cottenier@ugent.be">Stefaan.Cottenier@ugent.be</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><br>
My guess might be wrong, but it looks like you call two atoms &#39;bonded/unbonded&#39; if there is/isn&#39;t a connecting rod between them in a 3D display? If so, then you shouldn&#39;t worry. This is just a way of drawing (determined by the interatomic distance and the Z-value of the atoms involved). It has nothing to do with real chemical bonding. It is a display option that can be given as input in most (all) crystal viewers.<br>
<font color="#888888">
<br>
Stefaan</font><div class="im"><br>
<br>
<br>
On 26/10/2010 11:06, shamik chakrabarti wrote:<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Dear Peter Blaha Sir and wien2k users,<br>
<br>
                               We have implemented your suggestion and<br>
was able to remove the space group tolerence problem. But still I have a<br>
question regarding the structure. In the energy and force optimized<br>
structure there are some atomic networks and atoms which are not bonded<br>
with each other. The networks and atoms are remain disconnected with the<br>
main atomic network. If we repeat the unit cell along x, y and z<br>
direction, some networks get bonded but still there remain some free<br>
networks and atoms. Thise happened as we replaced C atom in A2BCO4 by D<br>
atom. Due to this replacement lattice parameter enhances but the space<br>
group symmetry remains same. Hence the atoms which remain bonded in<br>
A2BCO4 separated much from each other and get unbonded in A2BDO4. *Now<br>
my question is this types of unbonded atomic structure is really<br>
feasible although the energy and force are minimum in that structure??....*<br>
<br>
In some experimentally generated pure phase structure also we have seen<br>
such unbonded atoms and atomic network. As they are already<br>
experimentally prepared they are feasible although in the unit cell some<br>
atoms can be shown to be remaind unbonded. *Then how can we explain the<br>
presence of the unbonded atomic networks in this experimentally prepared<br>
materials??...*<br>
*<br>
*<br>
*with best regards,*<br>
*Shamik Chakrabarti<br>
*<br>
On Sat, Oct 23, 2010 at 12:21 PM, shamik chakrabarti<br></div><div class="im">
&lt;<a href="mailto:shamikphy@gmail.com" target="_blank">shamikphy@gmail.com</a> &lt;mailto:<a href="mailto:shamikphy@gmail.com" target="_blank">shamikphy@gmail.com</a>&gt;&gt; wrote:<br>
<br>
    Dear Peter Blaha Sir,<br>
<br>
                                 Thank you very much for your reply. I<br>
    have understand what you have said and will try to implement the<br>
    same. Yes it should solve the problem. Thank you sir, thank you very<br>
    much!!....<br>
<br>
    with best regards,<br>
    Shamik Chakrabarti<br>
<br>
<br>
    On Sat, Oct 23, 2010 at 12:17 PM, Peter Blaha<br></div><div><div></div><div class="h5">
    &lt;<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a> &lt;mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>&gt;&gt;<br>

    wrote:<br>
<br>
        Then it simply means, that after min_lapw your atoms arrived at<br>
        positions, which made them nearly equivalent, so that sgroup<br>
        does not know within its accuracy limits if two atoms are<br>
        equivalent or not.<br>
<br>
        Eventually &quot;correct&quot; positions manually, eg. setting atoms at<br>
        0.50002 to 0.50000  and so on.<br>
<br>
        Am 23.10.2010 08:42, schrieb shamik chakrabarti:<br>
<br>
            Dear Peter Blaha Sir,<br>
                                           I am using wien2k_08. No it<br>
            is not an<br>
            hexagonal lattice. Actually the space group is unknown. We<br>
            have taken a<br>
            structure A2BCO4 which has space group Pnma and then<br>
            *replace atom C<br>
            with atom D *considering *primitive space group*. We have<br>
            used the<br>
            atomic coordinates and lattice constants of A2BCO4. But we<br>
            have made all<br>
            the atoms inequivalent in the new material A2BDO4 by considering<br>
            primitive space group. It was done as we *do not want to<br>
            include the<br>
            symmetry constraints* for the atoms of this new material as<br>
            A2BDO4 may<br>
            have some different structure than A2BCO4. We want to see if<br>
            A2BDO4<br>
            structure is feasible or not and infact DFT has found a<br>
            structure with<br>
            minimum energy and minimum force on the atoms.<br>
<br>
            with best regards,<br>
            Shamik Chakrabarti<br>
            On Sat, Oct 23, 2010 at 11:48 AM, Peter Blaha<br>
            &lt;<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a><br>
            &lt;mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>&gt;<br>
            &lt;mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a><br>
            &lt;mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>&gt;&gt;&gt; wrote:<br>
<br>
                Are you using WIEN2k_10 ???  It should not happen anymore.<br>
<br>
                Hexagonal lattice ???  min_lapw could have produced atomic<br>
                positions, which are not completely equivalent anymore (like<br>
                0.3333333 and 0.666666) because of rounding errors.<br>
<br>
                Check your atomic positions manually.<br>
<br>
                Am 23.10.2010 06:22, schrieb shamik chakrabarti:<br>
<br>
                    Dear wien2k users,<br>
<br>
                                                 we have done structure<br>
            optimization<br>
                    (both<br>
                    volume and force) of a A2BCO4 type material. After<br>
            optimization the<br>
                    maximum force on an atom is *2.036 mRy/a.u.* After<br>
            optimization<br>
                    if we<br>
                    view the structure with xcrysden we can see that<br>
            some atoms remain<br>
                    unbonded i.e. one or two free atom and the rest of<br>
            the structure is<br>
                    connected through bond. Now my questions are:<br>
<br>
                               (1) Can such a structure be stable<br>
            although DFT has<br>
                    able to<br>
                    found a solution in which the energy of the unit<br>
            cell and forces<br>
                    on the<br>
                    atoms are minimum?<br>
<br>
                               (2) When I take this optimized struct<br>
            file and tried to<br>
                    initilize and run a SCF with this struct file, after<br>
            *x sgroup*<br>
                    command<br>
                    it will show :<br>
                    *Accuracy problem. Please run with different<br>
            tolerance (x sgroup<br>
                    -settol<br>
                    .00000100)*<br>
                    *<br>
                    *<br>
                                   Now if we ignore this problem or we<br>
            run the command x<br>
                    sgroup -settol 0.00000100, irrespective of that all<br>
            the atoms got<br>
                    deleted after x symmetry command!!<br>
<br>
                                   So, my question is whether this<br>
            optimized struct<br>
                    file is<br>
                    not feasible or is their any way to remove this<br>
            accuracy problem?<br>
<br>
                    any response will be very helpful for us. Thanking you,<br>
<br>
                    with regards,<br>
<br>
<br>
                    --<br>
                    Shamik Chakrabarti<br>
                    Research Scholar<br>
                    Dept. of Physics &amp; Meteorology<br>
                    Material Processing &amp; Solid State Ionics Lab<br>
                    IIT Kharagpur<br>
                    Kharagpur 721302<br>
                    INDIA<br>
<br>
<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>
            &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>&gt;<br>
            &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
            &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>&gt;&gt;<br>
<br>
            <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br>
<br>
                --<br>
                Peter Blaha<br>
                Inst.Materials Chemistry<br>
                TU Vienna<br>
                Getreidemarkt 9<br>
                A-1060 Vienna<br>
                Austria<br>
                +43-1-5880115671<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>
            &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>&gt;<br>
            &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
            &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>&gt;&gt;<br>
<br>
            <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br>
<br>
<br>
<br>
            --<br>
            Shamik Chakrabarti<br>
            Research Scholar<br>
            Dept. of Physics &amp; Meteorology<br>
            Material Processing &amp; Solid State Ionics Lab<br>
            IIT Kharagpur<br>
            Kharagpur 721302<br>
            INDIA<br>
<br>
<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>
            &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>&gt;<br>
            <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br>
<br>
        --<br>
        Peter Blaha<br>
        Inst.Materials Chemistry<br>
        TU Vienna<br>
        Getreidemarkt 9<br>
        A-1060 Vienna<br>
        Austria<br>
        +43-1-5880115671<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>
        &lt;mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>&gt;<br>
        <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br>
<br>
<br>
<br>
    --<br>
    Shamik Chakrabarti<br>
    Research Scholar<br>
    Dept. of Physics &amp; Meteorology<br>
    Material Processing &amp; Solid State Ionics Lab<br>
    IIT Kharagpur<br>
    Kharagpur 721302<br>
    INDIA<br>
<br>
<br>
<br>
<br>
--<br>
Shamik Chakrabarti<br>
Research Scholar<br>
Dept. of Physics &amp; Meteorology<br>
Material Processing &amp; Solid State Ionics Lab<br>
IIT Kharagpur<br>
Kharagpur 721302<br>
INDIA<br>
<br>
<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" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
</div></div></blockquote><div><div></div><div class="h5">
<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" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Shamik Chakrabarti<br>Research Scholar<br>Dept. of Physics &amp; Meteorology<br>Material Processing &amp; Solid State Ionics Lab<br>IIT Kharagpur<br>Kharagpur 721302<br>
INDIA<br>
</div>