A bit late but I didnt answer in time.
This is one of the cases where I still wheel out ALMN to do the self rotation. It generates ALL symmetry peaks and allows you to select which axis you want to take as the polar axis. (ncode = 1 c*, ncode = 2 a*, ncode = 3 b* and so on)
And the program output also makes an attempt to explain what it is doing!

A script
#!/bin/csh -f

#
almn \
hklin  /y/people/ccp4/projects/mao/maon5_all_scaleit1.mtz \
MAPOUT /y/work/ccp4//almn.map  \
<<eof
SELF 3 25
RESO  10 3.5
TITL  Jean
CRYS file 1 orth 1 flim 1 10000000000000
LABI FILE 1 F=F_natlo
LIMIT 0 180 5 1    !  Beta limit 90 because of symmetry
!                    If in doubt set Beta limit 180..
FIND 5 40
NOPR
MAP
END
eof



Dirk Kostrewa wrote:
Dear CCP4ers,

I've asked you about symmetry in stereographic projections of self-rotation functions, because I have in a monoclinic space group with beta=97 a peak for a NCS 7-fold axis at Phi=83, Psi=90, Kappa=51.4. In this self-rotation function, calculated with GLRF, the monoclinic b-axis is oriented from south-pole to north-pole. Thus, with Psi=90, the 7-fold lies in the a,c-plane, and I was wondering wether it points along the c-axis by some symmetry element, since 83=180-97. Pierre Rizkallah pointed to the fact that the self-rotation function is calculated in Patterson space which has an inversion centre, and Ian Tickle has pointed me to the POLARRFN documentation that discusses some of these symmetries. From this, I can, in my case, construct the following symmetry-equivalent positions:

(1) Original Peak: Phi, Psi, Kappa ---> 83, 90, 51.4 (2) General equivalent position: 180+Phi, 180-Psi, -Kappa ---> 263, 90, -51.4 (3) Crystallographic Dyad || Poles: 180+Phi, Psi, Kappa ---> 263, 90, 51.4 (4) (3) combined with (2): Phi, 180-Psi, -Kappa ---> 83, 90, -51.4

Thus, none of these combinations gives a Phi-angle of 97 degrees. However, meanwhile Liang Tong, the author of GLRF, explained to me, that in his convention, a positive Phi angle goes from +X (=a) to -Z (=-c*), which means that Phi=83 points really along the (-)c-axis (thus confirming my initial hypothesis).

             7-fold
              /
             /
            /
           b-------------- a
          /|
         / |
        /  |
       /   |
      c    c*

So, many thanks to all of you who replied and helped me to solve this puzzle! I would still be interested in a textbook or paper discussing symmetry in stereographic projections, though ...

Best regards,

Dirk.

Am 22.11.2007 um 14:40 schrieb Dirk Kostrewa:

Dear CCP4ers,

does any of you have a good reference describing the symmetry of crystallographic stereographic projections? There is a lot of literature describing rotational symmetry in Eulerian angular space, but I'm not aware of any for polar angles. In particular, I've calculated a self-rotation function for a crystal in space group C2 with a monoclinic beta-angle of 97 degrees in a convention where Phi is the angle from the x-axis (=a-axis) and Psi is the angle from the monoclinic b-axis. I get a beautiful peak for a seven-fold rotation axis at Phi=83 degrees and Psi=90 degrees. I think that the seven-fold NCS-axis should point along the crystallographic c-axis, but then I would expect Phi=97 and Psi=90. Presumably, there must be a symmetry with the term Phi' = 180-Phi, but I would like to _know_ it. I would be grateful for any pointer to a good reference describing symmetry in crystallographic stereographic projections, also for future cases.

Best regards,

Dirk.

*******************************************************
Dirk Kostrewa
Gene Center, A 5.07
Ludwig-Maximilians-University
Feodor-Lynen-Str. 25
81377 Munich
Germany
Phone:     +49-89-2180-76845
Fax:     +49-89-2180-76999
E-mail:    [EMAIL PROTECTED]
*******************************************************




*******************************************************
Dirk Kostrewa
Gene Center, A 5.07
Ludwig-Maximilians-University
Feodor-Lynen-Str. 25
81377 Munich
Germany
Phone:     +49-89-2180-76845
Fax:     +49-89-2180-76999
E-mail:    [EMAIL PROTECTED]
*******************************************************



Reply via email to