next up previous contents index
Next: Bulk Solvent Mask Up: Crystallographic Diffraction Data Previous: Example: Expand a

Partial Structure Factors

  Suppose that one wants to freeze a major portion of the stucture during refinement. Under certain conditions, one can save CPU time by storing the calculated structure factors of the frozen part of the molecule and then adding them to the calculated structure factors of the free part of the molecule. The computation of R values, crystallographic target functions, or electron density maps always refers to (cf. Chapter 12). Note, however, that is not added to when performing XREFin DO operations (Section 12.5).

In the following example, the structure factors of residues 1 and 2 are computed and written along with to the file ``partial.data":

 xrefin
   ...
   reflection @fobs.fob end
   selection=( resid 1:2 ) 
   update 
   do (FPART=FCALC)   
   write reflection fobs fpart output=partial.data end        
 end

In the next example, the partial structure factors are read, the structure factors of the remaining residues are computed, and R for whole structure is computed:

 xrefin
   ...
   reflection reset @partial.data end
   selection=( not ( resid 1:2 ) ) 
   update
   print rfactor
 end

In addition, one should fix the frozen atoms during minimization and dynamics by including the statement

constraints fix=( resid 1:2 ) end
before running minimization or dynamics. One should use ``method=direct" instead of ``method=FFT" if the frozen part contains more than 80% of the molecule.



Sat Mar 11 09:37:37 PST 1995