<div dir="ltr">So I am thinking I want to do vc-relax allowing A,B,C to vary in Base-centered orthorhombic while keeping the volume fixed. <div><br></div><div>I know we can do cell_dofree in the pw.x input file and can modify cell_base.f90 if an option that we want isn't available to us. </div><div><br></div><div>The issue I'm having is that we can do </div><div><br></div><div>in cell_base.f90 I'll keep the cell volume fixed and vary the vectors of the primitive ORCC cell as so..</div><div><br></div><div>1 1 0</div><div>1 1 0 </div><div>0 0 1</div><div><br></div><div>with fix_volume=.true.</div><div><div><div>            </div></div><div>but if I use that iforceh with fix_volume then I will get cell_parameters that will vary like this:</div><div><br></div><div>(CP11 , CP12 , 0,       ),</div><div><div>(CP11 , CP12 , 0,       ),</div></div></div><div>(0       , 0        , CP33,),</div><div><br></div><div><br></div><div>CP11 != CP12 and CP21 != CP22</div><div><br></div><div>so I am thinking I am losing the ORCC primitive cell shape. Ideally I'd like to make it so </div><div><br></div><div><div>CP11 == CP12 and CP21 == CP22</div></div><div><br></div><div>while varying CP11,CP22, and CP33 and keeping a constant volume.</div><div><br></div><div>One option would be to transform the primitive ORCC cell and atomic_positions to it's conventional counterpart and do cell_dofree='xyz' with adding fix_volume=.true. to cell_base.f90.</div><div><br></div><div>I'd rather not do that just because it would take a while to do a vc-relax with 90 atoms (and I'm thinking I want to do this for a few dozen chemical variations of the same system)</div><div><br></div><div>Any suggestions?</div><div><br></div><div>Thanks,</div><div>Andrew</div><div><br></div><div><br></div></div>