<font size="4">Dear PWscf developers and users,<br><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</font><br><font size="4">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>-- <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>