<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif;font-size:small;color:#000000">Thank you for the clarification. <br clear="all"></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="trebuchet ms, sans-serif"><br></font></div><div dir="ltr"><font face="trebuchet ms, sans-serif">Regards,</font><div><font face="trebuchet ms, sans-serif"><br></font></div><div><b><font face="trebuchet ms, sans-serif">Chaman Gupta</font></b><div><font face="trebuchet ms, sans-serif">Graduate Research Student, <i><a href="http://depts.washington.edu/nrglab/" target="_blank">Novosselov Research Group</a>, UW</i><br></font></div><div><font face="trebuchet ms, sans-serif">Graduate Research Student, <i><a href="https://depts.washington.edu/pzlab/wordpress/" target="_blank">Pauzauskie Research Group</a>, UW</i></font></div><div><font face="trebuchet ms, sans-serif">University of Washington Seattle, <i>M.S. in Mechanical Engineering</i></font></div><div><font face="trebuchet ms, sans-serif">IIT Kharagpur, <i>B.Tech in Metallurgical and Materials Engineering</i></font></div><div><font face="trebuchet ms, sans-serif"><b><a href="https://www.linkedin.com/in/professorchamangupta" target="_blank">Linked In</a> | <a href="mailto:chaman@uw.edu" target="_blank">Email</a> | </b></font><span style="font-family:"trebuchet ms",sans-serif"><b>P</b>: +1 (206) 383-3514</span></div></div></div></div></div></div></div></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 20, 2020 at 5:50 AM M.J. Hutcheon <<a href="mailto:mjh261@cam.ac.uk">mjh261@cam.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="font-size:10pt">
<p>Dear Chaman,</p>
<p>When I say final geometry, I mean a,b,c,alpha,beta,gamma and all of the fractional coordinates of the atoms: enough to specify the structure in it's entirety. These can all change with pressure, as well as the unit cell shrinking, atoms sometimes move around when the crystal is pressurised.</p>
<blockquote type="cite" style="padding:0px 0.4em;border-left:2px solid rgb(16,16,255);margin:0px">
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">Or is it all accounted for when I just change the cell geometry (a,b,c and alpha,beta,gamma) to calculate the new band structure at a different pressure?</div>
</blockquote>
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"> </div>
<div dir="ltr">This is right; you do a geometry optimisation to determine a,b,c alpha,beta,gamma and the fractional atom coordinates at the new pressure. The geometry is then all you need to determine the bandstructure.</div>
<div dir="ltr"> </div>
<div dir="ltr">Best,</div>
<div dir="ltr"> </div>
<div dir="ltr">Michael</div>
<div dir="ltr"> </div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p id="gmail-m_-8627903021149689257reply-intro">On 2020-05-20 13:05, Chaman Gupta wrote:</p>
<blockquote type="cite" style="padding:0px 0.4em;border-left:2px solid rgb(16,16,255);margin:0px">
<div id="gmail-m_-8627903021149689257replybody1">
<div dir="ltr">
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">Great explanation. That clears up so many of my questions, and the in-general how to approach my problem.</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">When you say feed the final geometry, you are referring to the a,b,c and alpha,beta,gamma values right?</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)"> </div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">A fundamental question:</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)"> </div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">When I say that I have a structure and I apply external pressure to it, then I wanna see how the band structure changes. So the band structure is basically the electronic energy levels (valence and conduction). Now apart from the pressure changing my cell geometry, shouldn't there be an effect of this external pressure on these energies? Or is it all accounted for when I just change the cell geometry (a,b,c and alpha,beta,gamma) to calculate the new band structure at a different pressure?</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"><span style="font-family:"trebuchet ms",sans-serif"> </span></div>
<div dir="ltr"><span style="font-family:"trebuchet ms",sans-serif">Regards,</span>
<div><span style="font-family:"trebuchet ms",sans-serif"> </span></div>
<div><strong><span style="font-family:"trebuchet ms",sans-serif">Chaman Gupta</span></strong>
<div><span style="font-family:"trebuchet ms",sans-serif">Graduate Research Student, <em><a href="http://depts.washington.edu/nrglab/" rel="noopener noreferrer" target="_blank">Novosselov Research Group</a>, UW</em><br></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">Graduate Research Student, <em><a href="https://depts.washington.edu/pzlab/wordpress/" rel="noopener noreferrer" target="_blank">Pauzauskie Research Group</a>, UW</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">University of Washington Seattle, <em>M.S. in Mechanical Engineering</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">IIT Kharagpur, <em>B.Tech in Metallurgical and Materials Engineering</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif"><strong><a href="https://www.linkedin.com/in/professorchamangupta" rel="noopener noreferrer" target="_blank">Linked In</a> | <a href="mailto:chaman@uw.edu" rel="noreferrer" target="_blank">Email</a> | </strong></span><span style="font-family:"trebuchet ms",sans-serif"><strong>P</strong>: +1 (206) 383-3514</span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<div>
<div dir="ltr">On Wed, May 20, 2020 at 4:43 AM M.J. Hutcheon <<a href="mailto:mjh261@cam.ac.uk" rel="noreferrer" target="_blank">mjh261@cam.ac.uk</a>> wrote:</div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div style="font-size:10pt">
<p>Dear Chaman,</p>
<p>Firstly, the pressures in QE are in kbar, not GPa (I don't know if this is different in the GUI version), see <a href="https://www.quantum-espresso.org/Doc/INPUT_PW.html" rel="noopener noreferrer" target="_blank">https://www.quantum-espresso.org/Doc/INPUT_PW.html</a>. Because the cell changes shrinks (and possibly the atoms move) as you increase the pressure, you should do a geometry optimisation at each new pressure, before calculating anything. The final geometry at each pressure can then be fed into further calculations, which now no longer need to know about the external pressure, as they only depend on the geometry (this is why the scf calculation will not accept a pressure). Anything that QE can calculate can be calculated as a function of pressure by this method.</p>
<p>Best,</p>
<p>Michael</p>
<p id="gmail-m_-8627903021149689257v1gmail-m_-8147440817993983812reply-intro">On 2020-05-20 12:22, Chaman Gupta wrote:</p>
<blockquote style="padding:0px 0.4em;border-left:2px solid rgb(16,16,255);margin:0px">
<div id="gmail-m_-8627903021149689257v1gmail-m_-8147440817993983812replybody1">
<div dir="ltr">
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">Thanks a lot Michael for the explanation. Diving a bit more into my specific situation, here's my problem statement:</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)"> </div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">I have a 'Si simple hexagonal phase (Phase V)' which is stable in the pressure range of 17 to 30 GPa (approx.). I want to generate the band structure for this Si at different Pressure values.</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">*Note: My cell dimensions will vary with change in Pressure*</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)"> </div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">1. I have the cif file for this phase. Should I in this case run a geometry optimization and look at the HOMO LUMO energy difference to compute the band gap,</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">     or, Should I run a scf calculation and look at the HOMO LUMO energy difference to compute the band gap.</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)"> </div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">2. When I am running a 'scf' calculation, and I try to enter (P = 19.5 GPa):</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">    <span style="font-family:Arial,Helvetica,sans-serif;color:rgb(34,34,34)">&CELL</span></div>
    cell_dofree    = "all"<br>    cell_dynamics  = "bfgs"<br>    press          =  1.95000e+02<br>    press_conv_thr =  5.00000e-01
<div><br>
<div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">The BURAI does not take in this command when I hit update. Whilst, if I enter the same command while doing a geometry optimization it does update my input file.</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">a. Can this command only be included while doing a geometry optimization calculation ('relax' or 'vc-relax')?</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">b. Is there any way to include external pressure while doing the scf calculation?</div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)"> </div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">3. Is it even possible to study (generate) band structures at different external pressures using QE?</div>
<div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)"> </div>
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">Thanks again for the help.</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"><span style="font-family:"trebuchet ms",sans-serif"> </span></div>
<div dir="ltr"><span style="font-family:"trebuchet ms",sans-serif">Regards,</span>
<div><span style="font-family:"trebuchet ms",sans-serif"> </span></div>
<div><strong><span style="font-family:"trebuchet ms",sans-serif">Chaman Gupta</span></strong>
<div><span style="font-family:"trebuchet ms",sans-serif">Graduate Research Student, <em><a href="http://depts.washington.edu/nrglab/" rel="noopener noreferrer" target="_blank">Novosselov Research Group</a>, UW</em><br></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">Graduate Research Student, <em><a href="https://depts.washington.edu/pzlab/wordpress/" rel="noopener noreferrer" target="_blank">Pauzauskie Research Group</a>, UW</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">University of Washington Seattle, <em>M.S. in Mechanical Engineering</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">IIT Kharagpur, <em>B.Tech in Metallurgical and Materials Engineering</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif"><strong><a href="https://www.linkedin.com/in/professorchamangupta" rel="noopener noreferrer" target="_blank">Linked In</a> | <a href="mailto:chaman@uw.edu" rel="noreferrer" target="_blank">Email</a> | </strong></span><span style="font-family:"trebuchet ms",sans-serif"><strong>P</strong>: +1 (206) 383-3514</span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<div>
<div dir="ltr">On Wed, May 20, 2020 at 4:01 AM M.J. Hutcheon <<a href="mailto:mjh261@cam.ac.uk" rel="noreferrer" target="_blank">mjh261@cam.ac.uk</a>> wrote:</div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div style="font-size:10pt">
<p>Dear Chaman,</p>
<p>Essentially, a geometry optimisation moves the atoms around to reduce the energy of the system. The energy of the system is calculated using an scf calculation. The geometry optimisation is therefore a bunch of scf calculations strung together, with a succession of atomic positions approaching the minimum-energy positions.</p>
<p>You should carry out a geometry optimisation if you suspect that your atoms are not in the ground state for your given parameters.</p>
<p>Best,</p>
<p>Michael Hutcheon</p>
<p>TCM group, Cambridge University</p>
<p><br></p>
<p id="gmail-m_-8627903021149689257v1gmail-m_-8147440817993983812v1gmail-m_-3508657660197467182reply-intro">On 2020-05-20 11:37, Chaman Gupta wrote:</p>
<blockquote style="padding:0px 0.4em;border-left:2px solid rgb(16,16,255);margin:0px">
<div id="gmail-m_-8627903021149689257v1gmail-m_-8147440817993983812v1gmail-m_-3508657660197467182replybody1">
<div dir="ltr">
<div style="font-family:"trebuchet ms",sans-serif;font-size:small;color:rgb(0,0,0)">
<div>Hey everyone,</div>
<div> </div>
<div>Might be a very naive question, but I have been very confused with this recently.</div>
<div> </div>
<div>I have a 'Si - simple hexagonal' cif file. I want to see the band structure and/or the band gap (HOMO - LUMO).</div>
<div>1. I ran a 'scf' calculation on BURAI (a GUI for QE). It generated an output file which gave me the Total energy and the energy of HOMO and LUMO.</div>
<div>2. Now separately, I ran a 'geometry optimization' calculation (both 'relax' and 'vc-relax': thus 2 different calculations), again on BURAI, and I got an output file, which gives me Total energy and energy of HOMO and LUMO.</div>
<div> </div>
<div>*I specified nbnd = 8 (no. of bands) for both the above mentioned calulations*</div>
<div> </div>
<div>I am confused as in which one is correct or more reliable. And also when should we use which calculation?</div>
<div> </div>
<div>Help is really appreciated.</div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"><span style="font-family:"trebuchet ms",sans-serif"> </span></div>
<div dir="ltr"><span style="font-family:"trebuchet ms",sans-serif">Regards,</span>
<div><span style="font-family:"trebuchet ms",sans-serif"> </span></div>
<div><strong><span style="font-family:"trebuchet ms",sans-serif">Chaman Gupta</span></strong>
<div><span style="font-family:"trebuchet ms",sans-serif">Graduate Research Student, <em><a href="http://depts.washington.edu/nrglab/" rel="noopener noreferrer" target="_blank">Novosselov Research Group</a>, UW</em><br></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">Graduate Research Student, <em><a href="https://depts.washington.edu/pzlab/wordpress/" rel="noopener noreferrer" target="_blank">Pauzauskie Research Group</a>, UW</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">University of Washington Seattle, <em>M.S. in Mechanical Engineering</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif">IIT Kharagpur, <em>B.Tech in Metallurgical and Materials Engineering</em></span></div>
<div><span style="font-family:"trebuchet ms",sans-serif"><strong><a href="https://www.linkedin.com/in/professorchamangupta" rel="noopener noreferrer" target="_blank">Linked In</a> | <a href="mailto:chaman@uw.edu" rel="noreferrer" target="_blank">Email</a> | </strong></span><span style="font-family:"trebuchet ms",sans-serif"><strong>P</strong>: +1 (206) 383-3514</span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<div style="margin:0px;padding:0px;font-family:monospace">_______________________________________________<br>Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noopener noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noopener noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div>
</blockquote>
<p><br></p>
</div>
</blockquote>
</div>
</div>
</blockquote>
<p><br></p>
</div>
</blockquote>
</div>
</div>
</blockquote>
<p><br></p>

</div>
</blockquote></div>