Skip to content

Commit b90ecfe

Browse files
committed
fix menu name for RMCprofile; add help on restraint weighting
1 parent 017e024 commit b90ecfe

File tree

5 files changed

+24
-12
lines changed

5 files changed

+24
-12
lines changed

MDhelp/docs/commontreeitems.md

Lines changed: 13 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -154,7 +154,7 @@ Also included when sequential refinement is selected is a menu button labeled "S
154154

155155
This window shows the restraints to be used in a refinement for each phase (if more than one). It is organized into several tabbed pages, one page for each type of restraint. Restraints are developed for an individual phase and act as additional observations to be "fitted" during the refinement.
156156

157-
Note that a restraint "pushes" a refinement towards a target value, but does not require that to happen (unlike a constraint). The strength of the "push" is dependent on the weighting factor and the "esd" used on the target value. The ideal use of a restraint is for an aspect of a structure that is not well-determined by the data. In this case, when the weighting of a restraint is dropped, the value of the restrained parameter may move considerably from the target value, but the quality of the fit (noted by the Rietveld plot, R-factor, reduced $\chi^2$, etc.) is not expected to change very much. If one places a restraint that is not consistent with the model (*e.g.* is wrong), then the fit will improve significantly as the weight is decreased. Such feedback is not possible with a constraint.
157+
Note that a restraint "pushes" a refinement towards a target value, but does not require that to happen (unlike a constraint). The strength of the "push" is dependent on the weighting factor and the "esd" used on the target value. The ideal use of a restraint is for an aspect of a structure that is not well-determined by the data. In this case, when the weighting of a restraint is dropped, the value of the restrained parameter may move considerably from the target value, but the quality of the fit (noted by the Rietveld plot, R-factor, reduced $\chi^2$, etc.) is not expected to change very much. If one places a restraint that is not consistent with the model (*e.g.* is wrong), then the fit will improve significantly as the weight is decreased. Such feedback is not possible with a constraint. [See below](#RestraintWeights) for more discussion on restraint weighting.
158158

159159
<H3 style="color:blue;font-size:1.1em">What can I do here?</H3>
160160

@@ -230,6 +230,18 @@ angle, ignore, C9 C10 C11, OK, ignore, ignore, 116.18,
230230
angle, ignore, C11 C12 N15, OK, ignore, ignore, 123.91, 0.02
231231
```
232232

233+
<a name="RestraintWeights"></a>
234+
### Restraint Weights
235+
236+
There is no simple recipe for how to weight restraints, in that it varies with the problem, but one should understand is that restraints are used computationally just like any other data point, where the refinement works towards minimizing the sum of
237+
$[ w * (obs_i - calc_i ) / \sigma_i ]^2$ where w is usually 1 and $obs_i$ might be a powder diffraction intensity but could also be a restraint bond distance and where $\sigma_i$ is the uncertainty on $obs_i$, which for restraints is the "esd" value. Since there are usually $10^3$ to $10^4$ data points and only circa $10^1$ restraints, I will usually initially weight the restraints quite highly (100-10,000) since I want to start the fit with a model that closely matches the expected distances and angles.
238+
239+
As the refinement progresses and I have a good fit, I will then try lowering the restraint weighting. Ideally, I can set the weight on the distances & angles to zero, in which case the restraints will show me a diagnostic of how far the as-fit distances, etc. deviate from my presupposed values, but do not actually change the fit.
240+
241+
A weight factor of 0 may not produce a good result with complex structures and not very complex diffraction patterns. In those cases, there just is not enough information to define all the structural degrees of freedom uniquely. In that case, lowering the weighting will increase the deviations seen in the restraints with only very trivial improvements in the fit to the powder pattern. One needs to reduce the restraint weight, but only to the point where a chemically reasonable structure is still obtained. Usually, a small value for the weight is needed (around 0.1 to 10) for the final refinements, as very little "force" is needed to keep the structure in agreement with what is expected for the class of materials.
242+
243+
There is one important case to be aware of, and that is where the restraints are incorrect for the material being studied. If one tries to restrain, for example, a bond to a value that is wrong for the material (for example a carbon-carbon single bond where in fact a double-bond is present) and the data have enough sensitivity to determine this, then the restraints and the data will be in conflict. Lowering the restraint weighting will produce a significant improvement in the fit to the powder data as the bond distance refines away from the target value. If this occurs the restraints need to be corrected to the valid structure. This is different from the case where there is not enough information to define all the structural degrees of freedom, where the fit to the powder data improves only slightly as the restraint is down-weighted. In that latter case, there are very many structures that produce about the same quality of fit, but restraints will prejudice the fit towards models that are more expected based on what is known on the class of materials.
244+
233245

234246
<a name="Rigidbodies"></a>
235247
<a name="Rigid_bodies"></a>

RMCProfile-I/RMCProfile-I.htm

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1860,7 +1860,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 2. Set general
18601860
minor-latin;mso-bidi-theme-font:minor-latin'>1</span></b><span
18611861
style='mso-fareast-font-family:"Times New Roman"'> minute. At each save time a
18621862
number of files are written by RMCProfile; these can be viewed by using the
1863-
Operations/View command (more about this later – don’t bother trying it now,
1863+
Operations/Plot command (more about this later – don’t bother trying it now,
18641864
there is nothing to see). For the purposes of this tutorial leave them at their
18651865
defaults. <o:p></o:p></span></p>
18661866

@@ -2213,7 +2213,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 7. Viewing
22132213
<p class=MsoNormal><span style='mso-fareast-font-family:"Times New Roman"'>Do </span><b
22142214
style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif;
22152215
mso-ascii-theme-font:minor-latin;mso-fareast-font-family:"Times New Roman";
2216-
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/View</span></b><span
2216+
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/Plot</span></b><span
22172217
style='mso-fareast-font-family:"Times New Roman"'>; a FileDialog showing only
22182218
*.csv files will appear<o:p></o:p></span></p>
22192219

RMCProfile-II/RMCProfile-II.htm

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1873,7 +1873,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 2. Set general
18731873
minor-latin;mso-bidi-theme-font:minor-latin'>1</span></b><span
18741874
style='mso-fareast-font-family:"Times New Roman"'> minute. At each save time a
18751875
number of files are written by RMCProfile; these can be viewed by using the
1876-
Operations/View command (more about this later – don’t bother trying it now,
1876+
Operations/Plot command (more about this later – don’t bother trying it now,
18771877
there is nothing to see). For the purposes of this tutorial leave them at their
18781878
defaults. <o:p></o:p></span></p>
18791879

@@ -2282,7 +2282,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 7. Viewing
22822282
<p class=MsoNormal><span style='mso-fareast-font-family:"Times New Roman"'>Do </span><b
22832283
style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif;
22842284
mso-ascii-theme-font:minor-latin;mso-fareast-font-family:"Times New Roman";
2285-
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/View</span></b><span
2285+
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/Plot</span></b><span
22862286
style='mso-fareast-font-family:"Times New Roman"'> (you can do this any time
22872287
after a save has been done); a FileDialog showing only *.csv files will appear<o:p></o:p></span></p>
22882288

@@ -4113,7 +4113,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 8. Fixing peak
41134113
the result do </span><b style='mso-bidi-font-weight:normal'><span
41144114
style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin;
41154115
mso-fareast-font-family:"Times New Roman";mso-hansi-theme-font:minor-latin;
4116-
mso-bidi-theme-font:minor-latin'>Operations/View</span></b><span
4116+
mso-bidi-theme-font:minor-latin'>Operations/Plot</span></b><span
41174117
style='mso-fareast-font-family:"Times New Roman"'>; the peaks in the </span><b
41184118
style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif;
41194119
mso-ascii-theme-font:minor-latin;mso-fareast-font-family:"Times New Roman";

RMCProfile-III/RMCProfile-III.htm

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1878,7 +1878,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 2. Set general
18781878
minor-latin;mso-bidi-theme-font:minor-latin'>1</span></b><span
18791879
style='mso-fareast-font-family:"Times New Roman"'> minute. At each save time a
18801880
number of files are written by RMCProfile; these can be viewed by using the
1881-
Operations/View command (more about this later – don’t bother trying it now,
1881+
Operations/Plot command (more about this later – don’t bother trying it now,
18821882
there is nothing to see). For the purposes of this tutorial leave them at their
18831883
defaults. <o:p></o:p></span></p>
18841884

@@ -2314,7 +2314,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 7. Viewing
23142314
<p class=MsoNormal><span style='mso-fareast-font-family:"Times New Roman"'>Do </span><b
23152315
style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif;
23162316
mso-ascii-theme-font:minor-latin;mso-fareast-font-family:"Times New Roman";
2317-
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/View</span></b><span
2317+
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/Plot</span></b><span
23182318
style='mso-fareast-font-family:"Times New Roman"'>; a <span class=SpellE>FileDialog</span>
23192319
showing only *.csv files will appear<o:p></o:p></span></p>
23202320

@@ -4712,7 +4712,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 8. Adjustment
47124712
"Times New Roman"'> on the nag note. RMCProfile will run for another 10 min.
47134713
When done do </span><b style='mso-bidi-font-weight:normal'><span
47144714
style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin;
4715-
mso-fareast-font-family:"Times New Roman";mso-hansi-theme-font:minor-latin'>Operations/View</span></b><span
4715+
mso-fareast-font-family:"Times New Roman";mso-hansi-theme-font:minor-latin'>Operations/Plot</span></b><span
47164716
style='mso-fareast-font-family:"Times New Roman"'> and examine the G(R)
47174717
partials plot<o:p></o:p></span></p>
47184718

@@ -4739,7 +4739,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 8. Adjustment
47394739
new setting. When done (another 10 min) do </span><b style='mso-bidi-font-weight:
47404740
normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:
47414741
minor-latin;mso-fareast-font-family:"Times New Roman";mso-hansi-theme-font:
4742-
minor-latin'>Operations/View</span></b><span style='mso-fareast-font-family:
4742+
minor-latin'>Operations/Plot</span></b><span style='mso-fareast-font-family:
47434743
"Times New Roman"'> to see the results.<o:p></o:p></span></p>
47444744

47454745
<p class=MsoNormal><span style='mso-fareast-font-family:"Times New Roman";

RMCProfile-IV/RMCProfile-IV.htm

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1882,7 +1882,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 2. Set general
18821882
minor-latin;mso-bidi-theme-font:minor-latin'>1</span></b><span
18831883
style='mso-fareast-font-family:"Times New Roman"'> minute. At each save time a
18841884
number of files are written by RMCProfile; these can be viewed by using the
1885-
Operations/View command (more about this later – don’t bother trying it now,
1885+
Operations/Plot command (more about this later – don’t bother trying it now,
18861886
there is nothing to see). For the purposes of this tutorial leave them at their
18871887
defaults. <o:p></o:p></span></p>
18881888

@@ -2311,7 +2311,7 @@ <h3><span style='mso-fareast-font-family:"Times New Roman"'>Step 7. Viewing
23112311
<p class=MsoNormal><span style='mso-fareast-font-family:"Times New Roman"'>Do </span><b
23122312
style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif;
23132313
mso-ascii-theme-font:minor-latin;mso-fareast-font-family:"Times New Roman";
2314-
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/View</span></b><span
2314+
mso-hansi-theme-font:minor-latin;mso-bidi-theme-font:minor-latin'>Operations/Plot</span></b><span
23152315
style='mso-fareast-font-family:"Times New Roman"'>; a <span class=SpellE>FileDialog</span>
23162316
showing only *.csv files will appear<o:p></o:p></span></p>
23172317

0 commit comments

Comments
 (0)