[phenixbb] How does phenix.fmodel treat isotropic and anisotropic B?

Pavel Afonine pafonine at lbl.gov
Tue Sep 14 17:01:58 PDT 2010


  Hi Hailiang,

all phenix.fmodel uses from input PDB file are:

- CRYST1 record;
- Coordinates, occupancies and scattering types from ATOM records;
- ADPs defined in ANISOU records (if ANISOU records are available) or it 
takes isotropic equivalent from ATOM records.

That's all you need to compute Fcalc (but not Fmodel), which is:

   Fmodel = scale * exp(AnisoScale) * (Fcalc + k_sol * exp(-b_sol*s^2/4) 
* Fmask)

   where:

   - Fmodel - total model structure factor (complex value)
   - AnisoScale = -ht*A(-1)*b_cart*A(-1)th/4
   - h - column vector with Miller indices
   - A - orthogonalization matrix
   - b_cart - anisotropic scale matrix
   - t and (-1) denotes transposition and inversion operations
   - scale - overall scale factor
   - Fcalc - structure factors calculated from atomic model
   - k_sol and b_sol - Flat Bulk solvent model parameters
   - Fmask - structure factors calculated from bulk solvent mask

Note, the isotropic equivalent B-factor in ATOM records is the mean of 
the trace of the ANISOU matrix divided by 10000 and multiplied by 8*pi^2.

Pavel.


On 9/14/10 4:48 PM, zhangh1 at umbc.edu wrote:
> Hi,
>
> I am doing a phenix.fmodel calculation for a PDB structure containing
> anisotroipc B. If I understand it correctly, the B factor is composed of
> both the isotropic B for ATOM line (Biso) and the anisotropic B from
> ANISOU line (Uani):
>
> Fc=sigma ( fi * exp (-Biso*s*s) * exp (-st*Uani*s) * exp (i*r*s) )
>
> Please correct me if the above equation is wrong. I also noticed that in
> PDB, the Uani (after scaling by 10^-4) is generally much smaller than
> Biso. So I assume that Biso has a much bigger contribution than Uani, the
> latter only implying the vibrational anisotropicity.
>
> However, I did a test calculation by setting all Bio to 0 whereas keeping
> Uani unchanged. Due to the above assumption, I expected the results should
> be very different from the native structure; however, coot told me they
> were just so close...
>
> Maybe I am wrong somewhere...
>
> Best Regards, Hailiang
>
> _______________________________________________
> phenixbb mailing list
> phenixbb at phenix-online.org
> http://phenix-online.org/mailman/listinfo/phenixbb




More information about the phenixbb mailing list