<div dir="ltr">Hello,<div><br></div><div>I am a graduate assistant at Virginia Commonwealth University and I am using Quantum Espresso to perform CP molecular dynamics calculations. I would like to run a simple set up of a single Carbon Diamond crystal structure. My ultimate research goal is to generate equations of state for warm dense matter, starting with Carbon in a diamond structure. It is my understanding that there are multiple ways to input this system into Quantum Espresso but I am not getting the same results with each method. I want to make sure I have an intimate understanding of the cp.x executable. This is why I want to calculate the same system using different inputs. (I am using QE 7.0) I have tried many different input methods with the main differences being the number of atoms that are explicitly defined. (2,8,& 18) </div><div><br></div><div>It is my understanding that if you use input parameter ibrav=2 where 2 atoms are explicitly defined under "ATOMIC_POSITIONS", this is the carbon diamond primitive cell. Based on QE related literature, cp.x will automatically fill in the remaining atoms to complete the carbon diamond crystal structure. I am skeptical that cp.x is actually doing this. Most of the QE related literature pertained to pw.x so perhaps each executable defines the crystal structure differently? I have looked into documentation for both cp.x and pw.x to no avail. I would appreciate anyone who can help clarify this input for me. </div><div><br></div><div>If you use input parameter ibrav=0 then you should be able to explicitly define all 18 atoms that comprise the conventional carbon diamond lattice structure. However, if I try to use ibrav=0 and explicitly define 18 atomic positions I get infinite values as outputs. My next attempt was to use ibrav=0 and explicitly define 8 atomic positions. The 8 atom run gave reasonable values but the pressure was orders of magnitude smaller than the 2 atom run. This is concerning to me and it has become clear that I am misunderstanding the cp.x input. I have scoured QE forums and documentation to no avail. </div><div><br></div><div>At this point it would seem the most accurate method is the 8 atom run since I can use these positions to generate a supercell that has the periodic carbon diamond structure, albeit with vacancies in some of the corners and sides. I assume if I make the supercell large enough the errors due to imperfections will be negligible. Also, no real structure is perfectly periodic so we should still be able to use the 8 atom unit cell to construct at usable supercell for modeling purposes. Does this assumption sound reasonable?</div><div><br></div><div>To summarize by questions:</div><div>What is the most accurate way to input the carbon diamond crystal structure in cp.x?</div><div>Why does ibrav=0 with 18 atoms explicitly defined not work?</div><div>Is my understanding of how to input the primitive cell in cp.x correct?</div><div>Is my assumption regarding the "8 atom generated supercell" reasonable?</div><div><br></div><div>I understand this is quite a bit to ask but I am at a loss for what to do next. I would appreciate any help that can be offered. I have provided example input files for each type of input described. (see google drive link)</div><div><br></div><div>Sincerely,</div><div>Juniper Savchick</div><div>Graduate Assistant</div><div>Virginia Commonwealth University</div><div>BS Physics, minor Mathematics</div><div>Pursuing a PhD in Mechanical and Nuclear Engineering</div><br><div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:#f5f5f5;padding:5px;color:#222;font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid #ddd;line-height:1"><a href="https://drive.google.com/drive/folders/1BLY8nVULkYzAWleKz-n6nwjFaKHc_U4m" target="_blank" style="display:inline-block;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0px;border:none;width:100%" aria-label="QE inputs"><img style="vertical-align: bottom; border: none;" src="https://drive-thirdparty.googleusercontent.com/16/type/application/vnd.google-apps.folder+shared"> <span dir="ltr" style="color:#15c;text-decoration:none;vertical-align:bottom">QE inputs</span></a></div><div> <br></div><div><br></div></div>