I General Questions
- [#gen1 My question is not answered here !]
- [#gen2 How can I get news about the evolution (new releases etc..) of Fox ?]
- [#gen3 What is an "Object" ?]
II Global Optimization
- [#glob1 I know 'Simulated Annealing' (SA), why do you recommend to use "Parallel Tempering" (PT) ?]
- [#glob2 My optimization has reached [enter cost here]. Has it converged enough ?]
- [#glob3 Fox does not optimize profile/unit cell parameters, even if I check the parameters to be optimized ! Why ?]
- [#glob4 What parameters can/should I optimize ?]
- [#glob5 What objects should I add in the Optimization object ?]
- [#glob6 How can I remove an object from an Optimization ?]
III Crystal Modelization
- [#cryst1 How can I take into account Special Positions ? What is the Dynamical Occupancy Correction ?]
- [#cryst2 What happens when two atoms are completely overlapping ? Does one disappear ?]
- [#cryst3 How can I ensure that the algorithm keeps the correct formula for my compound ?]
- [#cryst4 How can I use an Antibump criterion ? How can I see the Antibump parameters ? Does it prevent the "merging" of atoms using the Dynamical Occupancy Correction ?]
- [#cryst5 Should I restrict the positionnal parameters (fractionnal coordinates) to [0;1] using the parameter's limits ? Same for angles in [0,2pi]?]
- [#diff1 I copied the powder pattern 2theta zero (shift) from my profile fitting program, but it looks wrong ?]
- [#diff2 How can I change the powder pattern background points ?]
- [#diff3 Is it possible to exclude parts of a powder pattern ? Can I remove an exclude region ?]
- [#diff4 Can I use several data sets to refine a single crystal structure ?]
- [#diff5 You recommend to use only low-angle data, but why discard this information ? High-angle reflections ARE important for structure determination]
- [#diff6 What shall I put in the "linear polarization rate" field ?]
- [#diff7 Do you take into account the magnetic scattering of neutrons ?]
[#diff8 The Rwp and graphs look fine, but the optimized preferred orientation parameters are wrong !]
Anchor(gen1) My question is not answered here !
Send me an e-mail at firstname.lastname@example.org, I generally answer quickly ! And I like to have queries, that's how the software & algorithms get improved !
Anchor(gen2) How can I get news about the evolution (new releases etc..) of Fox ?
To be kept informed about new releases, you are encouraged to subscribe to the low-volume FOX mailing list at http://lists.sourceforge.net/mailman/listinfo/objcryst-foxx (the archives can be found at http://sourceforge.net/mailarchive/forum.php?forum_id=654).
If you intend to use the ObjCryst++ / Fox source code for your own projects, it is highly recommended to subscribe to the ObjCryst++ mailing list at http://lists.sourceforge.net/mailman/listinfo/objcryst-devel, (the archives can be found at http://sourceforge.net/mailarchive/forum.php?forum_id=1689).
Anchor(gen3) What is an "Object" ?
Almost everything (Crystal, PowderPattern, Scatterer, ScatteringPower, OptimizationObj) is considered an object. In Fox/!ObjCryst++ a generic approach is used, so that the algorithms can optimize them without actually knowing what they actually are: the idea is that the algorithm knows that these objects are "Refinable Objects", so that he can get from them (i) a list of parameters and (ii) a choice of Cost Functions. The advantage of such an approach is that if the modelization (parametrization) of some objects changes, the algorithms do not need any modification since they use a generic approach to the objects.
II Global Optimization
Anchor(glob1) I know 'Simulated Annealing' (SA), why do you recommend to use "Parallel Tempering" (PT) ?
In SA, the 'temperature' of the algorithm (think: high 'temperatures' allows for more improbable configurations, allowing to get of local minima) is decreased slowly following a predetermined law. If the decrease is too fast, you can easily be trapped in a local minimum.With PT, all temperatures are optimized in parallel, so that you can get of a local minimum at all times. The algorithm is invariant with time and therefore does not require from the user to choose the number of trial and the temperature schedule.General Idea: keep the default algorithm and optimization schedule & parameters (temperature and displacement), they are tuned to work with any structure or data. In other words, always use Parallel Tempering, unless you are interested in algorithms and want to play with SA.
Anchor(glob2) My optimization has reached [enter cost here]. Has it converged enough ?
Sorry, there is no general rule for just numbers... Nothing (for a rough idea at least) replaces the eye comparison of calculated and observed patterns (use the zoom by left-click & drag !). A Rwp value of 0.15 can be sufficient (for the structure solution to be good) if the background is low, but it can also be crap if the background is high (with a background at 50 000 and the maximum intensity at 60 000, even Rwp=.05 can be very bad) !
Anchor(glob3) Fox does not optimize profile/unit cell parameters, even if I check the parameters to be optimized ! Why ?
- "It's not a bug, it' a feature": Fox automatically fixes all profile and 2theta correction parameters, for several reasons:
- It is a bad idea to use a Global Optimization algorithm derived from Monte-Carlo to search for this kind of parameters.
- Profile parameters, 2theta shift (zero) can be determined before using Fox for the structure solution, using profile fitting (and therefore you should do it that way).
- By using the "integrated Rwp" (iRwp) cost function (rather than R or Rwp), the algorithm is almost insensitive to profile parameters. Personnally, for all structures I work on and for all examples given with Fox, I never use precise/refined profile parameters, but always adjust them "by eye" within Fox, which is sufficient thanks to the "i"Rwp approach. And generally I only use the 'W' term (constant width). You only need to refine these parameters to get a good full-profile agreement and please referees, but for the structure solution (finding approximate positions of atoms) part, only the integrated intensities matter.
Anchor(glob4) What parameters can/should I optimize ?
- You can optimize positionnal parameters (translations, rotations and conformation), Biso's, occupancies, and texture parameters. However, you should begin to optimize positionnal parameters only (translations, rotations and conformation), after setting "expected" values for Biso, occupancies of 1, and without texture parameters. You may activate the optimization of occupancies after a first optimization if you believe it is necessary. Optimizing Biso is probably a bad idea, since at low angle (the only part used for global optimization) data is not sensitive to Biso.
And do not activate texture parameters unless (i) you have a strong feeling it is necessary and (ii) you have not been able to obtain an un-textured data set. Having strong texture will significantly reduce the probability of finding the structure solution.
Anchor(glob5) What objects should I add in the Optimization object ?
- If you are optimizing a Crystal structure with a single Powder Pattern, you should add these two objects. If you have several data sets, you must add all diffraction data objects.
Anchor(glob6) How can I remove an object from an Optimization ?
- The only way to remove an object or cost function is to edit an xml file, and look at the very end: just remove the line corresponding to the object/cost function you want to remove, and voila ! NOTE: do not remove an object if you are still using on of its cost functions !!
III Crystal Modelization
Anchor(cryst1) How can I take into account Special Positions ? What is the Dynamical Occupancy Correction ?
- To take into account special positions, activate the "use dynamical occupancy correction" option in your Crystal structure interface. Do not change the occupancy of any atom. When atoms (of identical type) overlap, their occupancy will be automatically reduced (in the background, you won't see it) so that n atoms overlapping will all have an occupancy of 1/n. This allows to correct any type of overlapping, for an atom on a special position, or atoms being merged. It is strongly recommended to activate the Dynamical Occupancy Correction for all inorganic structures. Generally it will be useless for molecular structures, as special positions or shared atoms between molecules are extremenly rare.
Anchor(cryst2) What happens when two atoms are completely overlapping ? Does one disappear ?
- When (say) two atoms fully overlap (assuming you have activated the Dynamical Occupancy Correction), both their occupancies are set to 0.5, so that the sum of the two atoms look just like one. However both the atoms are still present in the global optimization, and they will be able to move away from each other later if it gives a better solution.
Anchor(cryst3) How can I ensure that the algorithm keeps the correct formula for my compound ?
- You have to put at least as many independent atoms as there are in the final result ; that is, you have to guess approximately how many atoms will be required. And thanks to the Dynamical Occupancy Correction, you can put too many atoms, and excess atoms will automatically be merged. Of course, don't exaggerate the number of excess atoms or the algorithm may converge very slowly.
Anchor(cryst4) How can I use an Antibump criterion ? How can I see the Antibump parameters ? Does it prevent the "merging" of atoms using the Dynamical Occupancy Correction ?
- To use the antibump:
you need to add (Crystal, menu "Parameters -> Set Antibump Distance") the couple of ScatteringPower (atom types) and declare the interatomic distance below which a penalty is used (this should approximately be the expected distance between the two atoms). Do that for all couples of atom types you want to affect with an antibump.
- the cost will automatically be added to the overall cost (as long as you added the Crystal structure tpo the list of objects to be optimized). If the cost is too small with regard to the Chi^2 of your diffraction data, you can use the "scale" factor displayed beside the antibump cost to increase it. Be warned that a too large antibump factor will prevent the optimization to work correctly; therefore, do not mistake the antibump function for an energetic description of your crystal structure.
You cannot see the antibump parameters, unless you edit an xml save file (where each line <AntiBumpDistance ScattPow1="O" ScattPow2="Pb">1.2</!AntiBumpDistance> corresponds to an antibump distance). You can change an antibump distance by re-setting the antibump distance using the Crystal "Parameters -> Set Antibump Distance" menu. To delete it simply remove the corresponding line.
For identical types of atoms, the antibump penalty vanishes to zero when the atoms get near each over (practically, the penalty goes from reaches its maximum when the distance decreases from the given dmin to dmin/2, and then decreases to 0 when the atoms fully merge). So the antibump cost function does not prevent merging. Note that you cannot use different antibump parameters for (say) two different oxygen types: as the antibump evaluation currently uses the atomic number, all oxygens will necessary have the same antibump parameters.
Anchor(cryst5) Should I restrict the positionnal parameters (fractionnal coordinates) to (0;1) using the parameter's limits ? Same for angles in (0,2pi)?
- No, these are "periodic" parameters, and are automatically corrected, e.g. changed to 0.01 if they get to 1.01 (using limits is then a bad idea, as it restricts movements).
IV Diffraction data
Anchor(diff1) I copied the powder pattern 2theta zero (shift) from my profile fitting program, but it looks wrong ?
- Try to change the sign and/or to multiply it by a factor 2, and see if it is better (after each change, right-click in the powder graph window to update the display).
Anchor(diff2) How can I change the powder pattern background points ?
- Simply change the 2-column (2theta - intensity) file you used to import the points, and re-import it. It does not matter if the number of points used has changed.
Anchor(diff3) Is it possible to exclude parts of a powder pattern ? Can I remove an exclude region ?
Yes, in the PowderPattern Object, use the Pattern->"Add 2theta exclude region", and enter the 2theta limits. You can add several ones.
NOTE: the powder patterns in these regions are still computed, so if you want to avoid computing the higher parts of the powder pattern, rather use the "Max sin(theta)/lambda" entry field.
To remove it, you must save in an xml file, and then remove the line declaring the region, which should be something like : <Exclude2Theta>0 5</!Exclude2Theta>
Anchor(diff4) Can I use several data sets to refine a single crystal structure ?
Yes, generally speaking, you can add any number of data sets. You just have to add them all to the GlobalOptimization object, and declare a cost function (the iRwp preferably) for each data set.
Anchor(diff5) You recommend to use only low-angle data, but why discard this information ? High-angle reflections are important for structure determination !
- First of all, a clear distinction should be made between (i) structure solution and (ii) structure refinement, which are both part of a complete structure determination:
- "Structure solution" consists in starting from relatively little information about the structure (one data set, and chemical composition), and one way or another to determine the approximate (0.1-0.5 Angstroem) position of all atoms.
- "Structure refinement" consist in starting from known approximate atoms positions, and to determine these positions, as well as temperature factors, occupancies, with the best precision possible.
Anchor(diff6) What shall I put in the "linear polarization rate" field ?
- The linear polarization rate is around 1.0 for synchrotron data (typically 0.95-0.98), and 0 for laboratory diffractometer (unless you use a monochromator). A precise value is not critical for global optimization, as we mostly use low-angle data, and as errors can be swallowed by the overall B-factor. The calculation of the polarization factor assumes that the polarization plane is horizontal, and the diffracting plane is vertical (i.e. standard equatorial geometry to maximize intensity on a synchrotron). This is ignored for neutron diffraction.
Anchor(diff7) Do you take into account the magnetic scattering of neutrons ?
Anchor(diff8) The Rwp and graphs look fine, but the optimized preferred orientation parameters are wrong !
- With preferred orientation, it is easy to find several close configuration, by inverting the March coefficient and using a perpendicular vector of preferred orientation. This cannot happen if you use the limits on the March coefficient, either to be above (needles) or below (plate) 1.0.