<div dir="ltr"><div>Thanks, Giovanni and Jess. I am trying to use Cu (101) and TiO2 (111).  So now I am trying this:<span style="font-size:24px;font-family:Arial;color:rgb(67,67,67);background-color:transparent;font-weight:400;font-style:normal;font-variant:normal;text-decoration:none;vertical-align:baseline" id="docs-internal-guid-73f07742-e63e-9877-2edc-0b44de2d097a"></span><br><br></div><br><div>&CONTROL<br>  calculation='scf',<br>  outdir='.',<br>  prefix='calc',<br>  pseudo_dir = './',<br>  verbosity='high',<br>  tstress=.true.,<br>  tprnfor=.true.,<br>/<br><br>&SYSTEM<br>  ibrav=0,<br>  celldm(1)=30.7407795807d0,<br>  nat=35,<br>  ntyp=3,<br>  ecutwfc=40,<br>  ecutrho=200,<br>  occupations='smearing',<br>  smearing='gaussian',<br>  degauss=0.01<br>/<br><br>&ELECTRONS<br>  diagonalization='david',<br>  conv_thr=1d-08,<br>  mixing_mode='plain',<br>  mixing_beta=0.7,<br>/<br><br>ATOMIC_SPECIES<br>  Cu 63.546000d0 Cu.pw-mt_fhi.UPF<br>  O 15.999400d0  O.pw-mt_fhi.UPF<br>  Ti 47.867000d0 Ti.pw-mt_fhi.UPF<br><br>ATOMIC_POSITIONS {alat}<br>  Cu   0.2785048041d0  -0.2776939344d0   0.6164550154d0<br>  Cu   0.5007270263d0  -0.2776939344d0   0.6164550154d0<br>  Cu   0.7229492485d0  -0.2776939344d0   0.6164550154d0<br>  Cu   0.2785048041d0  -0.4348287746d0   0.6164550154d0<br>  Cu   0.5007270263d0  -0.4348287746d0   0.6164550154d0<br>  Cu   0.7229492485d0  -0.4348287746d0   0.6164550154d0<br>  Cu   0.1673936930d0  -0.1991265143d0   0.6950224355d0<br>  Cu   0.3896159152d0  -0.1991265143d0   0.6950224355d0<br>  Cu   0.6118381374d0  -0.1991265143d0   0.6950224355d0<br>  Cu   0.1673936930d0  -0.3562613545d0   0.6950224355d0<br>  Cu   0.3896159152d0  -0.3562613545d0   0.6950224355d0<br>  Cu   0.6118381374d0  -0.3562613545d0   0.6950224355d0<br>  O    0.3936860281d0  -0.4600585799d0   0.7186400989d0<br>  Ti   0.7529053026d0  -0.4539572050d0   0.7252748967d0<br>  O    0.4454579103d0  -0.2907209900d0   0.7319096945d0<br>  O    0.2148034173d0  -0.4236541955d0   0.7582271981d0<br>  Ti   0.5740226917d0  -0.2604179804d0   0.7648619959d0<br>  O    0.2665752994d0  -0.2543166056d0   0.7714967938d0<br>  O    0.7025874730d0  -0.2301149708d0   0.7978142973d0<br>  Ti   0.3951400808d0  -0.2240135960d0   0.8044490951d0<br>  O    0.7543593552d0  -0.3750470614d0   0.8110838930d0<br>  O    0.5237048621d0  -0.5079802669d0   0.8374013965d0<br>  Ti   0.2162574699d0  -0.3447440518d0   0.8440361944d0<br>  O    0.5754767443d0  -0.3386426770d0   0.8506709922d0<br>  O    0.3448222512d0  -0.3144410423d0   0.8769884957d0<br>  Ti   0.7040415257d0  -0.3083396674d0   0.8836232936d0<br>  O    0.3965941334d0  -0.4593731329d0   0.8902580914d0<br>  O    0.8326063070d0  -0.2780366579d0   0.9165755950d0<br>  Ti   0.5251589148d0  -0.4290701233d0   0.9232103928d0<br>  O    0.2177115225d0  -0.4229687485d0   0.9298451906d0<br>  O    0.6537236961d0  -0.3987671137d0   0.9561626942d0<br>  Ti   0.3462763039d0  -0.3926657389d0   0.9627974920d0<br>  O    0.7054955783d0  -0.2294295238d0   0.9694322898d0<br>  O    0.4748410852d0  -0.3623627293d0   0.9957497934d0<br>  Ti   0.1673936930d0  -0.1991265143d0   1.0023845912d0<br><br>K_POINTS {automatic}<br>  8 8 8 1 1 1<br><br>CELL_PARAMETERS {alat}<br>  1.000000000000d0  0.000000000000d0  0.000000000000d0<br>  0.000000000000d0  -0.707106781187d0  0.000000000000d0<br>  0.000000000000d0  0.000000000000d0  1.618839606570d0<br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 5, 2016 at 2:55 PM, Giovanni Cantele <span dir="ltr"><<a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">giovanni.cantele@spin.cnr.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>The ibrav=0 is usually needed if you want to setup a unit cell that is not a default one in QE (see INPUT_PW.txt, under the description of ibrav).</div><div><br></div><div>If the unit cell you mentioned is just that you need (that is, if you have correctly setup the interface), it is orthorhombic, so in that case you can use the corresponding</div><div>ibrav. However, even with ibrav=0 setting up the vacuum is achieved just by increasing the length of the unit cell vector orthogonal to the interface.</div><div><br></div><div>If instead you question was also about how to build the interface, you should also specify which TiO2 and Cu surfaces you’re interested in.</div><div><br></div><div>Giovanni</div><br><div><blockquote type="cite"><div><div class="h5"><div>On 05 Apr 2016, at 12:09, Sarah Alpine <<a href="mailto:sarahalpine1@gmail.com" target="_blank">sarahalpine1@gmail.com</a>> wrote:</div><br></div></div><div><div><div class="h5"><div dir="ltr">If anyone has any suggestions, I'd really appreciate it. Thanks!<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 30, 2016 at 11:28 AM, Sarah Alpine <span dir="ltr"><<a href="mailto:sarahalpine1@gmail.com" target="_blank">sarahalpine1@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Dear All,<br></div>    I am simulating a Copper-TiO2 interface, and since Copper is FCC and TiO2 is tetragonal, I think I need to use ibrav=0. In that case, how do I set up the vacuum around the interface?<br></div>Thanks,<br></div>Sarah<br><br>&CONTROL<br>  calculation='vc-relax',<br>  outdir='.',<br>  prefix='calc',<br>  pseudo_dir = './',<br>  verbosity='high',<br>  tstress=.true.,<br>  tprnfor=.true.,<br>/<br><br>&SYSTEM<br>  ibrav=0,<br>  celldm(1)=20.4938530538d0,<br>  nat=35,<br>  ntyp=3,<br>  ecutwfc=40,<br>  ecutrho=200,<br>  occupations='smearing',<br>  smearing='gaussian',<br>  degauss=0.01<br>/<br><br>&ELECTRONS<br>  diagonalization='david',<br>  conv_thr=1d-08,<br>  mixing_mode='plain',<br>  mixing_beta=0.7,<br>/<br><br>ATOMIC_SPECIES<br>  Cu 63.546000d0 <a href="http://cu.pw" target="_blank">Cu.pw</a>-mt_fhi.UPF<br>  O 15.999400d0  <a href="http://o.pw" target="_blank">O.pw</a>-mt_fhi.UPF<br>  Ti 47.867000d0 <a href="http://ti.pw" target="_blank">Ti.pw</a>-mt_fhi.UPF<br><br>ATOMIC_POSITIONS {alat}<br>  Cu   0.1666666667d0  -0.1178511302d0   0.1178511302d0<br>...<br><br><br>K_POINTS {automatic}<br>  8 8 8 1 1 1<br><br>CELL_PARAMETERS {alat}<br>  1.000000000000d0  0.000000000000d0  0.000000000000d0<br>  0.000000000000d0  -0.471404520791d0  0.000000000000d0<br>  0.000000000000d0  0.000000000000d0  1.618839606572d0<br><br><br></div>
</blockquote></div><br></div></div></div><span class="">
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org" target="_blank">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></span></div></blockquote></div><br><div>
-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br>e-mail: <a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: <a href="http://people.na.infn.it/~cantele" target="_blank">http://people.na.infn.it/~cantele</a><br>

</div>
<br></div><br>_______________________________________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>