<div><DIV style="font-family:Georgia; color:#000000; font-size:13px;">Hello Cam,<DIV><BR></DIV><DIV>you should not set celldm(3) to 0. QE works with 3D periodic models. So if you are interested in graphene sheet, you have to find the c length so that the graphene sheet does not interact with its images in the c direction. c=10-15 angs should do the job but you have to check this.</DIV><DIV><BR></DIV><DIV>BTW, if you intend to relax the structure, the c parameter should be fixed. Use the cell_dofree keyword for this, as you mentioned.</DIV><DIV><BR></DIV><DIV>Best regards</DIV><DIV>Pascal<BR><BR><DIV>Cameron Foss <cjfoss@umass.edu> wrote:<BLOCKQUOTE type="cite" cite="<CANUcA0MuJbBVeAiJ5RTMassdrau0BTe1=ZR8BN2vu+GR3art=A@mail.gmail.com>"><DIV dir="ltr">Hello all,<DIV><BR></DIV><DIV>I have looked through a lot of the forum and have yet to find anything that helps with this question. However, if I am wrong where should I look for the answer to this?</DIV><DIV>
<BR></DIV><DIV>I am trying to run a scf simulation of a 2D graphene sheet of carbon atoms. Graphite is an hexagonal structure and uses ibrav=4, I have tried using this and setting celldm(3)=0, to indicate no height component, but that doesn't work. I have also tried indicating a free lattice (ibrav=0) and setting c=0, and setting the CELL_PARAMETERS card vectors to only having 2D components (ie the z-components are set to 0). However, that too has errors when I try running the pw.x code or loading the file in XCrySDen. </DIV>
<DIV><BR></DIV><DIV>(NOTE: I have been playing with QE-5.1 for only the past few weeks and am using XCrySDen to visualize the generated unit cells.) </DIV><DIV><BR></DIV><DIV>Particularly the error I get from XCrySDen in the terminal is: Num of Atoms =0, Num of Frames =0.</DIV>
<DIV><BR></DIV><DIV>I have also looked at the cell_dofree variable in the CELL card. but setting that to xy does not seem to work either, which I would suspect since I still have to define  a height to the structure.</DIV>
<DIV> </DIV><DIV>Are any of these approaches a viable way of getting just a 2D structure of graphene and I might be missing something? Or should i be going about this in another way? </DIV><DIV><BR></DIV><DIV>Best,</DIV><DIV>
Cam</DIV><DIV><BR></DIV></DIV>
<PRE style="font-size:13px;">_______________________________________________<BR>Pw_forum mailing list<BR><A href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</A><BR><A href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</A></PRE></BLOCKQUOTE></DIV><BR><BR><DIV><DIV style="font-family:Arial; color:#000000; font-size:10px;"><FONT size="2"><FONT size="1"></FONT><BR>>-----------------<BR>Pascal Boulet<BR>Aix-Marseille University <BR>MADIREL Laboratory<BR>Avenue Normandie-Niemen
<BR>13397 Marseille Cedex 20
<BR>Email: pascal.boulet@univ-amu.fr
<BR>Tel. <A href="tel:+33413551810" target="_blank">+33 413 55 18 10</A>
<BR>Fax  <A href="tel:+33413551850" target="_blank">+33 413 55 18 50</A></FONT><BR><BR></DIV></DIV></DIV></DIV></div>