Changeset 5151
- Timestamp:
- Jan 21, 2022 2:51:28 PM (21 months ago)
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
Tutorials/SeqRefine/SequentialTutorial.htm
r5039 r5151 2056 2056 <P></P> 2057 2057 <p class=MsoNormal> 2058 For the initial fitting of datasets over a range of conditions, we2059 will want to start from the fitted values we have fit, but only refine2060 a minimal number of parameters until the scale factor, background and2061 lattice parameters have been optimized for each dataset and more2062 variables can be added back.2063 2058 <P></P> 2064 2059 <p class=MsoNormal>There is one significant way that sequential … … 2070 2065 2071 2066 <P></P> 2072 <p class=MsoNormal><I> 2067 <p class=MsoNormal style='PADDING-LEFT:0.25in'><I> 2068 More details: 2073 2069 Note that GSAS-II shows the user unit cell parameters in the 2074 2070 traditional a, b, c, α, β & γ (distances in Å, angles in … … 2095 2091 </I> 2096 2092 2097 2098 2093 <H3>1. Remove refinement of sensitive parameters (lattice, phase 2099 2094 fractions, size/μstrain, background peak and sample displacement) 2100 2095 and add refinement of the <tt>D<I>ij</i></tt> terms.</H3> 2096 2097 <p class=MsoNormal style='PADDING-LEFT:0.25in'> 2098 For the initial fitting of datasets over a range of conditions, we 2099 will want to start from the fitted values we have fit, but until 2100 we have reasonably good values for the scale factor, background and 2101 lattice parameters, we do not want to refine any thing else.</P> 2102 2101 2103 2102 2104 <OL style='PADDING-LEFT:0.75in' type="A"> … … 2460 2462 src="SequentialTutorial_files/image075.png" v:shapes="Picture_x0020_77"><![endif]></span></p> 2461 2463 2462 <P></P><LI> 2464 <P></P><LI> 2465 <B><I>Important:</I></B> 2463 2466 Check the "<b>Copy results to next histogram</b>?" checkbox. When this 2464 2467 is selected, the parameters obtained after fitting each histogram 2465 2468 will be used as the starting point for the next fit, which is 2466 2469 essential here since the lattice parameters change considerably. If 2467 this were not selected, then refinement of each histogramwould2468 start with the current parameters -- that were obtained from fitting 2469 the first (~7 K) dataset. 2470 The first few refinements will go well, but when reaching datasets2470 this were not selected, then refinement for <I>each histogram</I> would 2471 start with the current parameters for that histogram -- 2472 these were obtained from fitting the first (~7 K) dataset. If we did that 2473 the first few refinements will go well, but when reaching datasets 2471 2474 where the temperature 2472 2475 was sufficiently different from 7 K, the starting lattice parameters may be 2473 2476 too far off and the refinements will fail. Even if that were not 2474 true, since the temperatures increase in the series, the previous 2475 histogram will have the closest set of lattice parameters and will 2476 thus refine most quickly. </p> 2477 true, since the temperatures are fairly close between datasets in the 2478 series, 2479 each fit will start with a fairly close to correct set of lattice 2480 parameters and thus fewer cycles of refinement will be needed to get 2481 to the correct values. 2482 </p> 2477 2483 2478 2484 <P></P><LI>Also on the Controls window, change the <B>Max cycles</B> … … 2493 2499 once the histograms were selected for a sequential fit in the previous 2494 2500 step, the menu command was relabeled from Calculate/Sequential 2495 refine.) The refinement may take a few minutes as each of the 17 2501 refine.)</P><P></P> 2502 <p class=MsoListParagraphCxSpMiddle <p 2503 style='PADDING-LEFT:0.25in'><I>Note that a window 2504 saying "There were constraint warnings..." will be generated. 2505 This is because we have specified a constraint on the phase 2506 fractions in <A href="#I.3">step I.3</A>, but we are no longer varying 2507 them, so this constraint must be ignored. 2508 Also note the discussion on <A href="#VI.1.F">sequential fit 2509 constraint modes</A> in step VI, below. 2510 This alert is shown to provide a chance to review 2511 possible problems with constraints and not continue with the fit, 2512 when this indicates that something has been left out from the varied 2513 parameters. However, in this case, this phase fraction constraint 2514 will be fine to ignore since we do not want to refine phase fractions 2515 at this stage. This message will stop being shown when phase fractions 2516 are refined again (or if the constraint would be removed). Press 2517 "Yes" to continue with this initial sequential fit, as this 2518 warning can be ignored. 2519 <P></P></I> 2520 <p class=MsoListParagraphCxSpMiddle> 2521 The refinement may take a few minutes as each of the 17 2496 2522 datasets are fit in turn.</p> 2497 2498 2499 2523 <h3>3. Accept the sequential refinement results.</h3> 2500 2524 … … 2639 2663 2640 2664 <P></P><LI> 2665 <A name="VI.1.F"></A> 2641 2666 Select the <B>Constraints</B> data tree item and note that the 2642 2667 second line of controls that appears for sequential fits. … … 2644 2669 2645 2670 <P></P> 2646 <p class=MsoListParagraphCxSpMiddle> 2671 </OL> 2672 <p class=MsoListParagraphCxSpMiddle style='PADDING-LEFT:0.25in'> 2647 2673 <I>Note: 2648 2674 Refinement of phase … … 2651 2677 As of version 5038 of GSAS-II three modes were added for constraints 2652 2678 in sequential refinements. 2653 <UL style='PADDING-LEFT: 0.75in'>2679 <UL style='PADDING-LEFT:1.0in'> 2654 2680 <LI>In </I>"Ignore unless hist=*"<I> mode (which was 2655 2681 the only mode prior to this version), any constraints with a … … 2667 2693 </UL> 2668 2694 </I> 2669 </OL>2670 2695 2671 2696 <h3>2. Perform the refinement.</h3> -
trunk/GSASIIdataGUI.py
r5139 r5151 5610 5610 else: 5611 5611 l = str(len(allwarnings)) + ' histograms' 5612 msg = 'There were c hanges to constraints needed in {}. See console output or visitConstraints data tree item for more details. Continue with refinement?'.format(l)5613 dlg = wx.MessageDialog(self,msg,caption='Constraint changes',style=wx.YES_NO)5612 msg = 'There were constraint warnings generated in {}. See console output or visit the Constraints data tree item for more details. Continue with refinement?'.format(l) 5613 dlg = wx.MessageDialog(self,msg,caption='Constraint warnings',style=wx.YES_NO) 5614 5614 dlg.CenterOnParent() 5615 5615 result = wx.ID_NO
Note: See TracChangeset
for help on using the changeset viewer.