<div dir="ltr"><div><span style="font-family:georgia,serif;font-size:large">Dear Sir,</span><br></div><div><span style="font-family:georgia,serif;font-size:large"><br></span></div><div><font face="georgia, serif" size="4">Good afternoon and happy new year to all of you.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">I am new user here so I need a step by step reply.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Before sending this email I have read the posting rules.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">I used Quantum Espresso 6.4.1 compiled with ifort on a 16 core CPU.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Without having any previous experience on any quantum simulation tool, I spent almost six months to learn quantum espresso and for bulk system I feel comfortable.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">I have done some calculations with a binary semiconductor compound which was vc-relxed for ecutwcf (65), ecurho (52) and k-mesh(9 9 3 0 0 0). The convergence  criteria was less then 1meV for these parameters. Pseudo-potentials are taken from PSlibrary (recommended there).</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Now I want to do some calculation of its monolayer.</font></div><div><font face="georgia, serif" size="4">As I could not get any informative tutorials for this kind of calculations so  I am asking here.</font></div><div><font face="georgia, serif" size="4">My doubts are on k-mesh, &CELL and vacuum for the mono-layer.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">1. k-mesh:</font></div><div><font face="georgia, serif" size="4">Que. Should I keep the same k-mesh but reducing the k-point in vacuum direction to 1?</font></div><div><font face="georgia, serif" size="4">i.e. 9 9 1 0 0 0 (9 9 3 0 0 0 was for bulk with 1meV convergence).<br></font></div><div><font face="georgia, serif" size="4">I have done a k-mesh convergence test and  no significant changes observed after 9 9 1 0 0 0 k-grid.</font></div><div><font face="georgia, serif" size="4"></font></div><div><font face="georgia, serif" size="4">My double arises here because materials cloud suggest 2 k-points in vacuum direction.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">2. &CELL</font></div><div><font face="georgia, serif" size="4">I have adopted two choices to do it.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">A. </font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">I have read the mailing list and <a href="https://lists.quantum-espresso.org/pipermail/users/2017-October/039359.html" style="">found a link</a> where it is mentioned that how to do vc-relax of a mono-layer.<br></font></div><div><font face="georgia, serif" size="4">That user used </font></div><div><font face="georgia, serif" size="4">&CELL<br>cell_dofree='2Dshape'<br>/<br></font></div><div><font face="georgia, serif" size="4">to do a vc-relax.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">With these flags (with k-mesh 9 9 1 0 0 0),</font></div><div><font face="georgia, serif" size="4">The relaxed data are below</font></div><div><font face="georgia, serif" size="4">Total force =     0.000008     Total SCF correction =     0.000000<br></font></div><div><font face="georgia, serif" size="4">          total   stress  (Ry/bohr**3)                   (kbar)     P=   -2.81                # Stress<br></font></div><div><font face="georgia, serif" size="4">     Final enthalpy =    -227.5797893621 Ry<br>Begin final coordinates<br>     new unit-cell volume =   1782.72164 a.u.^3 (   264.17209 Ang^3 )<br>     density =      0.85971 g/cm^3<br><br>CELL_PARAMETERS (alat=  1.00000000)<br>   6.675929970   0.000000000   0.000000000     # Please note a,b,c are same what I have supplied in the input file.<br>   0.000000000   6.675929970   0.000000000<br>   0.000000000   0.000000000  40.000000000<br></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">B.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">from pw.x input descriptions, I found "cell_dofree='2Dxy' " do to a vc-relax in x-y plane.</font></div><div><font face="georgia, serif" size="4">So with below flags</font></div><div><font face="georgia, serif" size="4">&CELL<br>cell_dofree='2Dxy'<br>/<br></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">My results are</font></div><div><font face="georgia, serif" size="4">     Total force =     0.000013     Total SCF correction =     0.000002<br></font></div><div><font face="georgia, serif" size="4">               total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.05                # Stress is 0.5 kbar<br></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">     Final enthalpy =    -227.5803704614 Ry<br>Begin final coordinates<br>     new unit-cell volume =   1742.44890 a.u.^3 (   258.20428 Ang^3 )<br>     density =      0.87958 g/cm^3<br><br>CELL_PARAMETERS (alat=  1.00000000)<br>   6.600092607   0.000000000   0.000000000           # Please note here that a and b has reduced from what I given into input file (6.675929970)<br>  -0.000000000   6.600092607   0.000000000<br>   0.000000000   0.000000000  40.000000000<br></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">My concern arises from the  stress and change in lattice parameters from above two cases.</font></div><div><font face="georgia, serif" size="4">In the case-B  where I have adopted stress and Final enthalpy is more stable then the case-A.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Que. Which approach I should use A or B? </font></div><div><font face="georgia, serif" size="4">Que. Is it okay if the cell parameters (a and b) deviates from the bulk case?</font></div><div><font face="georgia, serif" size="4">Que. Or you can advice me how to do vc-relax for  a mono-layer.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">C. For vacuum</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Que. How to decide the vacuum size?</font></div><div><font face="georgia, serif" size="4">I have done a convergence test from 30 bohr to 50 bohr and the total energy first decrease, reach to a lowest value of then starts to increase.</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Below is the Total every for 30 Bohr to 50 Bohr vacuum size</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">-227.40963239 Ry   # 30<br>-227.43879270 Ry   # 31<br>-227.46529446 Ry   # 32<br>-227.48906558 Ry<br>-227.51005151 Ry<br>-227.52821961 Ry<br>-227.54356349 Ry<br>-227.55610674 Ry<br>-227.56589822 Ry<br>-227.57301424 Ry<br>-227.57755075 Ry<br>-227.57962022 Ry                  # with vacuum 41 bohr<br>-227.57934725 Ry<br>-227.57686484 Ry<br>-227.57231253 Ry<br>-227.56583286 Ry<br>-227.55757040 Ry<br>-227.54767074 Ry<br>-227.53627785 Ry<br>-227.52353424 Ry    # 49<br>-227.50957968 Ry    # 50<br></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><span style="font-family:georgia,serif;font-size:large">Any step wise reply will help me a lot.</span><br></div><div><font face="georgia, serif" size="4"><br></font></div><div><span style="font-family:georgia,serif;font-size:large">Looking forward for a supportive response.</span><br></div><div><font face="georgia, serif" size="4"><br></font></div><div><span style="font-family:georgia,serif;font-size:large">thank you very much sir,</span><br></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Best wishes</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Ms. Rekha</font></div><div><span style="font-family:georgia,serif;font-size:large">Ex-PG student,</span></div><div><span style="font-family:georgia,serif;font-size:large">LBS college, Jaipur</span></div><div><span style="font-family:georgia,serif;font-size:large">India</span><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4">Mob.: +11 90-95 790 71 697</font></div><div><font face="georgia, serif" size="4">Email: <a href="mailto:rekha1997jpr@gmail.com">rekha1997jpr@gmail.com</a></font></div><div><font face="georgia, serif" size="4"><br></font></div><div><br></div></div>