Dear Lachlan,
I'll have a look at you tutorial. The macro stuff in GSAS has been
extremely valuable for the protein work. Setting up one is easy & takes
about 1/2 hr from scratch (macros & starting model ina PDB file). You may
have noticed the extensive set of amino acid based macros now distributed
with GSAS. Now for you questions/comments:
At 03:57 AM 10/31/00 +0000, you wrote:
>There is a new GSAS tutorial on the CCP14 site:
>
> Setting up GSAS Macro Files for Restrained Refinement on relatively
> large Inorganics (Polymeric Inorganics) (using a "created" Ti Fe O
> example - originally just a new published Titanate phase)
>
> http://www.ccp14.ac.uk/solution/gsas/restrained_inorganic.html
>
>This deals with:
> i) atom co-ordinates macro file
>
> ii) atom bond length restraints macro file (generated from Crystals
> with the option of performing DLS - and Scott Belmonte's new
> "coue" conversion software)
>
> iii) dual atom occupancy constraints for the metal sites
>
> iv) total cell contents restraints macro file
>
>Corrections, suggestions, etc appreciated. Especially new macros for
>making this type of thing easier to set up for certain types of
>restraints where it is important to test the robustness from different
>starting positions and refinement strategies.
>
>
>Lachlan.
>
>PS: While GSAS scales very well in terms of increased problem size and
>interface control - I am having some difficulty with getting convergence
>within Genles on a "real" dataset. (this is on Le Bail fitted data where
>only the background and atomic parameters are allowed to refine - starting
>from an idealised starting model guided from TEM)
Convergence for big problems is somewhat different from little ones. You
can now adjust the convergence criterion (it's a log scale!) and invoke a
very heavy duty global damping factor (Marquardt factor) to keep things
under control. The convergence factor probably should be scaled with the
size of the problem. I use a convergence criterion of ~50 for a 1000 atom
protein (>3000 variables). Note the default is 0.01.
>Does anyone have hints and tricks for parameter release conditions for
>large inorganics (that required restraints to keep the Oxygen distances
>reasonable) (~100 atoms in the asymmetric unit - may have to be doubled
>based on reinterpretation of the TEM data). All the metal sites have
>potential dual occupancy with that can be varied. What can you get away
>with and have a good change of convergence with Genles? How far away from
>their true positions can the heavy atoms be to drop in? Estimates of
>a good numbers of cycles to use to let the atoms rattle around?
Another feature of GSAS is the availability of a band matrix approximation
to the full matrix. Assuming that the atoms can be ordered by their
connectivity - easy for a protein, hard for a complex oxide; then a band
matrix can be used. I typically use a 300 variable band for a 3000 variable
problem. That will cover 10-12 amino acid xyz's which will catch all the
near-neighbor interactions captured in the restraints. With this the
convergence is pretty clean & smooth. Some fiddling then is done to
complete the refinement - mostly fixing messed up stereochemistry and then
repeating the refinement. A 1000 atom protein will converge in about 50-100
LS cycles (including fiddling) with the starting positions 1-2A away from
the final ones.
>Are there possibilities for what I believe on the commercial(?) version
>of Shelx 76(?) was a "cascading" refinement(?) - to automatically
>release small sets of atoms - and cycle through automatically?
>(is it possible to interlink expedt and genles under macro control
>to achieve this?)
I have done this for a protein - refine the structure using a sequence of
overlapping blocks doing one cycle of LS for each block. It can be setup
using DOS bat files and GSAS macros. It works but the convergence is slower
than using the band matrix approach.
By the way I've been doing this with powder data but all these facilities
are available for single crystal data as well. I haven't tried refining a
single crystal protein data set yet but it should work very well.
>--
>Lachlan M. D. Cranswick
>
>Collaborative Computational Project No 14 (CCP14)
> for Single Crystal and Powder Diffraction
>Daresbury Laboratory, Warrington, WA4 4AD U.K
>Tel: +44-1925-603703 Fax: +44-1925-603124
>E-mail: [EMAIL PROTECTED] Ext: 3703 Room C14
> http://www.ccp14.ac.uk