<div dir="ltr">On Tue, Sep 19, 2017 at 9:37 PM, Amar Singh <span dir="ltr"><<a href="mailto:amarsingh122014@rediffmail.com" target="_blank">amarsingh122014@rediffmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I tried your suggestion, with approximate values of lattice parameters for rhombohedral setting (as I have lattice parameters for only hexagonal settings) [...] if the unit cell given by me has a volume ~ 64 A3 which fits one BaTiO3 unit, how can space group generate two units of BaTiO3 for same volume ??<br></blockquote><div><br></div>maybe your approximate values of the lattice parameters are not well approximated<br><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">​My other problem is that the CIF file which I have (true data reported from X-ray diffraction) has lattice parameters for hexagonal setting<br>​ a = b = 5.5 A, c = 13 A and alpha = beta = 90, gama = 120. <br>​Now for this unit cell, R3c should be able to generate 6 units of BaTiO3. <br>​How can I write my input script to give me this unit cell optimized with 30 atoms (6 + 6 + 18)<br></blockquote><br></div><div class="gmail_quote">There are 1001 different ways to provide structural data to QE. For instance, you can do what I did years ago for trigonal and hexagonal Al2O3 (see below): write down atomic positions in  crystal axis in the two cases, write down the relation between the two set of axis.<br><br></div><div class="gmail_quote">Paolo<br></div><div class="gmail_quote">---<br>alpha-Al2O3:<br>R-3c (space group: 167, hexagonal axis) Z=6<br>positions Al: 12c  (0,0,z) (0,0,-z+1/2) (0,0,-z) (0,0,z+1/2)<br>positions  O: 18c  (x,0,1/4) (0,x,1/4) (-x,-x,1/4) <br>                   (-x,0,3/4) (0,-x,3/4) (x,x,3/4)<br>+ (0,0,0); (2/3,1/3,1/3); (1/3,2/3,2/3) (obverse setting) or<br>+ (0,0,0); (1/3,2/3,1/3); (2/3,1/3,2/3) (reverse setting)<br>z(Al)=0.35228, x(O)=0.3064<br>a=4.758, c=12.99<br><br>R-3c (space group: 167, rhombohedral axis) Z=2<br>positions Al: 4c  (z,z,z) (-z+1/2,-z+1/2,-z+1/2) (-z,-z,-z) (z+1/2,z+1/2,z+1/2)<br>positions  O: 6e  (x,-x+1/2,1/4) (1/4,x,-x+1/2) (-x+1/2,1/4,x)<br>                  (-x,x+1/2,3/4) (3/4,-x,x+1/2) (x+1/2,3/4,-x)<br>z(Al)=0.35228, x(O)=1/4-0.3064=-0.0564<br>a=5.13, alpha=55.6 (cos=0.565)<br><br>Relation between rhombohedral A_r, B_r, C_r and hexagonal A_h, B_h, C_h:<br><br>A_r = 2/3 A_h + 1/3 B_h + 1/3 C_h       A_h = A_r - B_r<br>B_r =-1/3 A_h + 1/3 B_h + 1/3 C_h       B_h = B_r - C_r<br>C_r =-1/3 A_h - 2/3 B_h + 1/3 C_h       C_h = A_r + B_r + C_r<br><br>a_r = sqrt(3a_h^2+c_h^2)/3 ; sin(alpha/2)= 3/(2*sqrt(3+(c_h/a_h)^2)<br><br>Relation between rhombohedral axis with threefold axis around 001:<br><br>a=( s*sqrt(3)/2, -s/2, r)   b=( 0, s, r)   c = (-s*sqrt(3)/2, -s/2, r)<br><br>s = 2a_0/sqrt(3) * sin(alpha/2),  r = sqrt(a_0^2-s^2)<br><br>and around 111: a=(p,q,q), b=(q,p,q), c=(q,q,p)   (p<q)<br><br>r = (p+2q)/sqrt(3) ;  s = (q-p)*sqrt(3/2)<br>q = (r+s/sqrt(2))/sqrt(3) ;  p = (r-s*sqrt(2))/sqrt(3)<br><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div></div>