<DIV>Dear Prof. Blaha and Prof. L. Marks, </DIV>  <DIV>Why can't an accurate code like&nbsp;WIEN&nbsp;perform </DIV>  <DIV>conjugate gradient relaxtion by computing </DIV>  <DIV>the full hessian. Is the computation of the </DIV>  <DIV>second derivative of the energy very </DIV>  <DIV>expensive?&nbsp;Relaxation with WIEN is </DIV>  <DIV>far difficult than I thought. PORT keeps </DIV>  <DIV>failing.&nbsp;It is not easy to find a </DIV>  <DIV>set of deltas and etas for NEWT. </DIV>  <DIV>Is there a thumb rule for selecting </DIV>  <DIV>delta and eta (for NEWT/NEW1) </DIV>  <DIV>for systems with entirely </DIV>  <DIV>different&nbsp;masses? e.g. H and Pt.</DIV>  <DIV>&nbsp;</DIV>  <DIV>Thank you.&nbsp;</DIV>