<font size="2">Dear PWscf developers and users,<br></font><div class="gmail_quote"><font size="2"><br>I am currently studying an epitaxial thin-film system. I would like to find the optimized the structure with the lattice parameters a and b being constrained, namely, finding the optimized lattice parameter c (and the internal structure).<br>

<br>I am aware that this procedure can be done manually. But that would cost lots of CPU hours.<br><br>I am considering using cell_dofree, and I am wondering how this option works. Suppose I put<br></font><font size="2">cell_dofree= 'z' under the namelist &CELL. How should I specify the target pressure?<br>

<br>For a thin film with constrained a, b and optimized c, the stress tensor should look like<br><br>sigma_xx = sigma_yy = a non-vanishing value.<br>sigma_xy = sigma_yx = also a non-vanishing value<br>sigma_zz = sigma_xz = sigma_yz = 0<br>

<br>But when specifying press= 0, I think (maybe I am wrong) the program is searching for a structure with <br><br>sigma_ij=0, <br>i,j=x,y,z.<br><br>How should I specify the target pressure (stress tensor) if I perform vc-relax with dofree= 'z'? <br>

<br>Thank you in advance for your response and hint. <br><br>Han</font><br><br clear="all"></div>-- <br>*************************************************************************<br>Han Hsu, Ph.D.<br>Postdoctoral Associate<br>
Department of Chemical Engineering & Materials Science<br>University of Minnesota<br><a href="http://www.cems.umn.edu/%7Ehsuhan" target="_blank">http://www.cems.umn.edu/~hsuhan</a><br>*************************************************************************<br>