[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [ccp4bb]: anisotropic ellipsoids
*** For details on how to be removed from this list visit the ***
*** CCP4 home page http://www.dl.ac.uk/CCP/CCP4/main.html ***
Norbert Straeter wrote:
>
> *** For details on how to be removed from this list visit the ***
> *** CCP4 home page http://www.dl.ac.uk/CCP/CCP4/main.html ***
>
> According to many textbooks the first three of the thermal parameters
> U11 U22 U33 U12 U13 and U23 describe the displacements along the
> perpendicular principal axis of the ellipsoid and the latter three give
> the orientation of the principal axes with respect to the unit cell axes.
> However, I can't find anywhere how U12 U13 and U23 (apparently as
> direction cosini) exactly describe the orientation of the ellipsoid, say
> in a cartesian system.
>
> Any hint is appreciated (but don't suggest to try to follow the ortep
> code)...
>
> Thanks, Norbert
Is this any help? Taken from $CLIBS/rwbrook.f
C PDB files contain cartesian anisotropic temperature factors as
orthogonal Us.
C Stored as U11 U22 U33 U12 U13 U23
C The anisotropic temperature factors can be input/output to this
routine
C as orthogonal or as crystallographic Us.
C
C Shelx defines Uf to calculate temperature factor as:
C T(aniso_Uf) = exp (-2PI**2 ( (h*ast)**2 Uf_11 + (k*bst)**2 Uf_22 + ...
C + 2hk*ast*bst*Uf_12 +..)
C
C Note: Uo_ji == Uo_ij and Uf_ji == Uf_ij.
C
C [Uo_ij] listed on ANISOU card satisfy the relationship:
C [Uo_ij] = [RFu]-1 [Uf_ij] {[RFu]-1}T C
C where [Rfu] is the normalised [Rf] matrix read from the SCALEi
cards.
C see code. [ROu] == [RFu]-1
C
C T(aniso_Uo) = U(11)*H**2 + U(22)*K**2 + 2*U(12)*H*K + ...
C where H,K,L are orthogonal reciprocal lattice indecies. ( EJD: I
think????)
C
C Biso = 8*PI**2 (Uo_11 + Uo_22 + Uo_33) / 3.0
C
C [Uf(symm_j)] = [Symm_j] [Uf] [Symm_j]T
C
--
------------------------------------------------------------------
Eleanor J.Dodson, Chemistry Department, University of York, U.K.
Tel: Home +44 (1904) 42 44 49, work: +44 (1904) 43 25 65
------------------------------------------------------------------