Changeset 1745 for Tutorials/MCsimanneal/MCSA in GSAS.htm
- Timestamp:
- Mar 20, 2015 8:49:59 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Tutorials/MCsimanneal/MCSA in GSAS.htm
r1718 r1745 23 23 <o:Author>vondreele</o:Author> 24 24 <o:LastAuthor>Von Dreele</o:LastAuthor> 25 <o:Revision> 2</o:Revision>26 <o:TotalTime> 581</o:TotalTime>25 <o:Revision>3</o:Revision> 26 <o:TotalTime>767</o:TotalTime> 27 27 <o:Created>2015-03-16T21:14:00Z</o:Created> 28 <o:LastSaved>2015-03-1 6T21:14:00Z</o:LastSaved>28 <o:LastSaved>2015-03-19T20:13:00Z</o:LastSaved> 29 29 <o:Pages>1</o:Pages> 30 <o:Words>24 20</o:Words>31 <o:Characters>13 797</o:Characters>30 <o:Words>2432</o:Words> 31 <o:Characters>13865</o:Characters> 32 32 <o:Company>Argonne National Laboratory</o:Company> 33 <o:Lines>11 4</o:Lines>33 <o:Lines>115</o:Lines> 34 34 <o:Paragraphs>32</o:Paragraphs> 35 <o:CharactersWithSpaces>16 185</o:CharactersWithSpaces>35 <o:CharactersWithSpaces>16265</o:CharactersWithSpaces> 36 36 <o:Version>15.00</o:Version> 37 37 </o:DocumentProperties> … … 44 44 <!--[if gte mso 9]><xml> 45 45 <w:WordDocument> 46 <w:Zoom>116</w:Zoom> 47 <w:SpellingState>Clean</w:SpellingState> 48 <w:GrammarState>Clean</w:GrammarState> 46 49 <w:TrackMoves>false</w:TrackMoves> 47 50 <w:TrackFormatting/> … … 644 647 {font-family:"Cambria Math"; 645 648 panose-1:2 4 5 3 5 4 6 3 2 4; 646 mso-font-charset: 1;649 mso-font-charset:0; 647 650 mso-generic-font-family:roman; 648 mso-font-format:other;649 651 mso-font-pitch:variable; 650 mso-font-signature: 0 0 0 0 00;}652 mso-font-signature:-536870145 1107305727 0 0 415 0;} 651 653 @font-face 652 654 {font-family:Calibri; … … 1098 1100 mso-fareast-theme-font:minor-fareast; 1099 1101 mso-bidi-font-family:Calibri;} 1102 span.SpellE 1103 {mso-style-name:""; 1104 mso-spl-e:yes;} 1105 span.GramE 1106 {mso-style-name:""; 1107 mso-gram-e:yes;} 1100 1108 .MsoChpDefault 1101 1109 {mso-style-type:export-only; … … 1104 1112 mso-ansi-font-size:10.0pt; 1105 1113 mso-bidi-font-size:10.0pt; 1106 font-family:"Calibri",sans-serif;1107 1114 mso-ascii-font-family:Calibri; 1108 1115 mso-hansi-font-family:Calibri; … … 1259 1266 <p class=MsoNormal>In these exercises you will use GSAS-II to solve the 1260 1267 structures of 3-aminoquinoline and α-d-lactose monohydrate from powder 1261 diffraction data via Monte Carlo/Simulated Annealing (MC/SA). The data sets were 1262 kindly provided by Peter Stephens (SUNY StonyBrook) and were originally 1263 collected on NSLS beam line X3b1. The MC/SA technique is needed in these cases 1264 because as you will see the data does not extend very far in 2Θ so one 1265 cannot use charge flipping to solve them. The steps for each begin with peak 1266 selection and fitting, then indexing to identify the lattice parameters. A 1267 Pawley refinement is used to obtain a set of structure factors needed for the 1268 MC/SA runs. We will use previously developed models for the two molecules; I 1269 used the freeware program Avogadro (<a href="http://avogadro.cc/wiki/Main_Page">http://avogadro.cc/wiki/Main_Page</a>) 1268 diffraction data via Monte Carlo/Simulated Annealing (MC/SA). The data sets 1269 were kindly provided by Peter Stephens (SUNY <span class=SpellE>StonyBrook</span>) 1270 and were originally collected on NSLS beam line X3b1. The MC/SA technique is 1271 needed in these cases because as you will see the data does not extend very far 1272 in 2Θ so one cannot use charge flipping to solve them. The steps for each 1273 begin with peak selection and fitting, then indexing to identify the lattice 1274 parameters. A Pawley refinement is used to obtain a set of structure factors 1275 needed for the MC/SA runs. We will use previously developed models for the two 1276 molecules; I used the freeware program Avogadro (<a 1277 href="http://avogadro.cc/wiki/Main_Page">http://avogadro.cc/wiki/Main_Page</a>) 1270 1278 to create them.</p> 1271 1279 1272 <p class=MsoNormal>If you have not done so already, <a 1273 href="../../../GSASII/help/Starting%20GSAS.htm">start GSAS-II</a>.</p> 1280 <p class=MsoNormal>If you have not done so already, start GSAS-II.</p> 1274 1281 1275 1282 <h2>Step 1. Import powder data for 3-aminoquinoline</h2> … … 1285 1292 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1286 1293 mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri;mso-bidi-theme-font: 1287 minor-latin'>Exercises/MCSA</span></b><span style='font-family:"Calibri",sans-serif; 1288 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family: 1289 Calibri;mso-bidi-theme-font:minor-latin'> to find the file; you may need to 1290 change the file type to <b style='mso-bidi-font-weight:normal'>All files (*.*)</b> 1291 to find the desired file.</span></p> 1294 minor-latin'>exercises/<span class=SpellE>MCsimanneal</span></span></b><span 1295 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1296 mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri;mso-bidi-theme-font: 1297 minor-latin'> to find the file; you may need to change the file type to <b 1298 style='mso-bidi-font-weight:normal'>All files (*.*)</b> to find the desired 1299 file.</span></p> 1292 1300 1293 1301 <p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo2'><![if !supportLists]><span 1294 1302 style='mso-fareast-font-family:"Times New Roman"'><span style='mso-list:Ignore'>2.<span 1295 1303 style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]>Select 1296 the <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1304 the <span class=SpellE><b style='mso-bidi-font-weight:normal'><span 1305 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1306 mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri;mso-bidi-theme-font: 1307 minor-latin'>Quinoli.gda</span></b></span><span style='font-family:"Calibri",sans-serif; 1297 1308 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family: 1298 Calibri;mso-bidi-theme-font:minor-latin'>Quinoli.gda</span></b><span 1299 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1300 mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri;mso-bidi-theme-font: 1301 minor-latin'> data</span> file in the first dialog and press <b 1302 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1309 Calibri;mso-bidi-theme-font:minor-latin'> data</span> file in the first dialog 1310 and press <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1303 1311 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family: 1304 1312 Calibri;mso-bidi-theme-font:minor-latin'>Open</span></b>. There will be a 1305 Dialog box asking < b style='mso-bidi-font-weight:normal'><span1313 Dialog box asking <span class=GramE><b style='mso-bidi-font-weight:normal'><span 1306 1314 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1307 1315 mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri;mso-bidi-theme-font: 1308 minor-latin'>Is this the file you want?</span></b> Press <b style='mso-bidi-font-weight: 1309 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1310 minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri; 1311 mso-bidi-theme-font:minor-latin'>Yes</span></b> button to proceed.</p> 1316 minor-latin'>Is</span></b></span><b style='mso-bidi-font-weight:normal'><span 1317 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1318 mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri;mso-bidi-theme-font: 1319 minor-latin'> this the file you want?</span></b> Press <span class=GramE><b 1320 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1321 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family: 1322 Calibri;mso-bidi-theme-font:minor-latin'>Yes</span></b></span> button to 1323 proceed.</p> 1312 1324 1313 1325 <p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo2'><![if !supportLists]><span … … 1316 1328 the <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1317 1329 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family: 1318 Calibri;mso-bidi-theme-font:minor-latin'>BNL- I.instprm</span></b> instrument1319 parameter file in the second dialog and press <b style='mso-bidi-font-weight:1330 Calibri;mso-bidi-theme-font:minor-latin'>BNL-<span class=SpellE>I.instprm</span></span></b> 1331 instrument parameter file in the second dialog and press <b style='mso-bidi-font-weight: 1320 1332 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1321 1333 minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family:Calibri; … … 1323 1335 file type to <b style='mso-bidi-font-weight:normal'><span style='font-family: 1324 1336 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1325 minor-latin'>*. instprm</span></b> (GSAS-II instrument parameter file). The1326 p owder pattern will immediately appear.</p>1337 minor-latin'>*.<span class=SpellE>instprm</span></span></b> (GSAS-II instrument 1338 parameter file). The powder pattern will immediately appear.</p> 1327 1339 1328 1340 <p class=MsoListParagraph><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shapetype … … 1346 1358 <v:path o:extrusionok="f" gradientshapeok="t" o:connecttype="rect"/> 1347 1359 <o:lock v:ext="edit" aspectratio="t"/> 1348 </v:shapetype><v:shape id=" _x0000_i1051" type="#_x0000_t75" style='width:525pt;1349 height:450pt;visibility:visible;mso-wrap-style:square'>1360 </v:shapetype><v:shape id="Picture_x0020_28" o:spid="_x0000_i1051" type="#_x0000_t75" 1361 style='width:525pt;height:450pt;visibility:visible;mso-wrap-style:square'> 1350 1362 <v:imagedata src="MCSA%20in%20GSAS_files/image001.png" o:title=""/> 1351 1363 </v:shape><![endif]--><![if !vml]><img border=0 width=700 height=600 1352 src="MCSA%20in%20GSAS_files/image002.gif" v:shapes=" _x0000_i1051"><![endif]></span></p>1364 src="MCSA%20in%20GSAS_files/image002.gif" v:shapes="Picture_x0020_28"><![endif]></span></p> 1353 1365 1354 1366 <p class=MsoListParagraph><o:p> </o:p></p> … … 1375 1387 picked 21 peaks, skipping a few very weak ones. Then use <b style='mso-bidi-font-weight: 1376 1388 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1377 minor-latin;mso-hansi-theme-font:minor-latin'>Peak fitting/LSQ PeakFit </span></b>to1378 fit them. You will need to vary at least <b style='mso-bidi-font-weight:normal'><span 1379 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin;1380 mso- hansi-theme-font:minor-latin'>X</span></b> & <b style='mso-bidi-font-weight:1381 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1382 m inor-latin;mso-hansi-theme-font:minor-latin'>Y</span></b> in <b1383 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif;1389 minor-latin;mso-hansi-theme-font:minor-latin'>Peak fitting/LSQ <span 1390 class=SpellE>PeakFit</span> </span></b>to fit them. You will need to vary at 1391 least <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1392 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>X</span></b> 1393 & <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1394 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Y</span></b> 1395 in <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1384 1396 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Instrument 1385 Parameters</span></b> along with <b style='mso-bidi-font-weight:normal'><span 1386 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1387 mso-hansi-theme-font:minor-latin'>pos</span></b> & <b style='mso-bidi-font-weight: 1388 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1389 minor-latin;mso-hansi-theme-font:minor-latin'>int</span></b> in the <b 1390 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1391 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Peak List</span></b> 1392 to get a reasonable fit suitable for indexing.</p> 1397 Parameters</span></b> along with <span class=SpellE><b style='mso-bidi-font-weight: 1398 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1399 minor-latin;mso-hansi-theme-font:minor-latin'>pos</span></b></span> & <span 1400 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 1401 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1402 minor-latin'>int</span></b></span> in the <b style='mso-bidi-font-weight:normal'><span 1403 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1404 mso-hansi-theme-font:minor-latin'>Peak List</span></b> to get a reasonable fit 1405 suitable for indexing.</p> 1393 1406 1394 1407 <h2>Step 3. Indexing</h2> … … 1402 1415 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1403 1416 minor-latin;mso-hansi-theme-font:minor-latin'>Unit Cells List</span></b> for 1404 indexing. The most likely Bravais lattice for these small organic molecules is 1405 either Orthorhombic-P or Monoclinic-P; you can try both and see what comes out 1406 (Hint: its <b style='mso-bidi-font-weight:normal'><span style='font-family: 1407 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1408 minor-latin'>Orthorhombic-P</span></b>). The correct unit cell should almost 1409 immediately appear with an M20 ~330 with a=7.748, b=7.651, c=12.736, 1410 Vol=755.01. Do <b style='mso-bidi-font-weight:normal'><span style='font-family: 1417 indexing. The most likely <span class=SpellE>Bravais</span> lattice for these 1418 small organic molecules is either Orthorhombic-P or Monoclinic-P; you can try 1419 both and see what comes out (Hint: its <b style='mso-bidi-font-weight:normal'><span 1420 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1421 mso-hansi-theme-font:minor-latin'>Orthorhombic-P</span></b>). The correct unit 1422 cell should almost immediately appear with an M20 ~200 with a=7.748, b=7.651, 1423 c=12.736, Vol=755.01. NB: your solution may have <span class=GramE>a &</span> 1424 b switched. Do <b style='mso-bidi-font-weight:normal'><span style='font-family: 1411 1425 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1412 1426 minor-latin'>Cell Index/Refine/Copy Cell</span></b>; you can then try various … … 1487 1501 box and enter the d-spacing (1.8133) into the <b style='mso-bidi-font-weight: 1488 1502 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1489 minor-latin;mso-hansi-theme-font:minor-latin'>Pawley dmin</span></b> box.</p> 1503 minor-latin;mso-hansi-theme-font:minor-latin'>Pawley <span class=SpellE>dmin</span></span></b> 1504 box.</p> 1490 1505 1491 1506 <p class=MsoNormal>When done the <b style='mso-bidi-font-weight:normal'><span … … 1494 1509 1495 1510 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1496 id=" Picture_x0020_2" o:spid="_x0000_i1050" type="#_x0000_t75" style='width:465.75pt;1497 height:375pt;visibility:visible;mso-wrap-style:square'>1511 id="_x0000_i1050" type="#_x0000_t75" style='width:487.5pt;height:375pt; 1512 visibility:visible;mso-wrap-style:square'> 1498 1513 <v:imagedata src="MCSA%20in%20GSAS_files/image003.png" o:title=""/> 1499 </v:shape><![endif]--><![if !vml]><img border=0 width=6 21height=5001500 src="MCSA%20in%20GSAS_files/image004.gif" v:shapes=" Picture_x0020_2"><![endif]></span></p>1514 </v:shape><![endif]--><![if !vml]><img border=0 width=650 height=500 1515 src="MCSA%20in%20GSAS_files/image004.gif" v:shapes="_x0000_i1050"><![endif]></span></p> 1501 1516 1502 1517 <p class=MsoNormal><o:p> </o:p></p> … … 1520 1535 1521 1536 <p class=MsoNormal>Select either choice; the desired data set will be added to 1522 this phase for analysis. The data tab now shows the new data set; check the <b 1523 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1524 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Show</span></b> 1525 button to see the full information</p> 1526 1527 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1528 id="Picture_x0020_28" o:spid="_x0000_i1048" type="#_x0000_t75" style='width:359.25pt; 1537 this phase for analysis. The data tab now shows the new data set.</p> 1538 1539 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1540 id="Picture_x0020_30" o:spid="_x0000_i1048" type="#_x0000_t75" style='width:359.25pt; 1529 1541 height:375pt;visibility:visible;mso-wrap-style:square'> 1530 1542 <v:imagedata src="MCSA%20in%20GSAS_files/image007.png" o:title=""/> 1531 1543 </v:shape><![endif]--><![if !vml]><img border=0 width=479 height=500 1532 src="MCSA%20in%20GSAS_files/image008.gif" v:shapes="Picture_x0020_ 28"><![endif]></span></p>1544 src="MCSA%20in%20GSAS_files/image008.gif" v:shapes="Picture_x0020_30"><![endif]></span></p> 1533 1545 1534 1546 <p class=MsoNormal><o:p> </o:p></p> 1535 1547 1536 1548 <p class=MsoNormal>This is the location for all the phase dependent parameters 1537 for this histogram. Notice that it includes phase fraction, size & mustrain 1538 as well as preferred orientation, elastic strain and extinction corrections. 1539 The Babinet correction is intended for protein work where a significant region 1540 of the structure is disordered solvent. There are buttons for plotting size and 1541 mustrain surfaces and preferred orientation correction curves.</p> 1549 for this histogram. Notice that it includes phase fraction, size & <span 1550 class=SpellE>mustrain</span> as well as preferred orientation, elastic strain 1551 and extinction corrections. The <span class=SpellE>Babinet</span> correction is 1552 intended for protein work where a significant region of the structure is 1553 disordered solvent. There are buttons for plotting size and <span class=SpellE>mustrain</span> 1554 surfaces and preferred orientation correction curves.</p> 1542 1555 1543 1556 <p class=MsoNormal>Next find the <b style='mso-bidi-font-weight:normal'><span … … 1549 1562 makes the reflection set over the range covered by setting <b style='mso-bidi-font-weight: 1550 1563 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1551 minor-latin;mso-hansi-theme-font:minor-latin'>Pawley dmin</span></b>. The table1552 should list reflections 0-91. Select and check the <b style='mso-bidi-font-weight: 1553 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1554 m inor-latin;mso-hansi-theme-font:minor-latin'>refine</span></b> column using1555 the same technique you used for the peak list. </p>1564 minor-latin;mso-hansi-theme-font:minor-latin'>Pawley <span class=SpellE>dmin</span></span></b>. 1565 The table should list reflections 0-91. Select and check the <b 1566 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1567 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>refine</span></b> 1568 column using the same technique you used for the peak list. </p> 1556 1569 1557 1570 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 1577 1590 mso-hansi-theme-font:minor-latin'>aminoquinoline.bck0.gpx</span></b>. It can be 1578 1591 used to recover from a bad refinement. A progress dialog box will appear 1579 showing the residual as the refinement proceeds. My refinement completed with 1580 Rwp ~12.9%; a dialog box appears asking if you wish to load the result. Press <b 1581 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1582 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>OK</span></b>. 1583 To see the plot select the PWDR line in the GSAS-II data tree (Ive adjusted 1584 the width and height).</p>1585 1586 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1587 id=" Picture_x0020_6" o:spid="_x0000_i1046" type="#_x0000_t75" style='width:525pt;1588 height:450pt;visibility:visible;mso-wrap-style:square'>1592 showing the residual as the refinement proceeds. My refinement completed with <span 1593 class=SpellE>Rwp</span> ~12.9%; a dialog box appears asking if you wish to load 1594 the result. Press <b style='mso-bidi-font-weight:normal'><span 1595 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1596 mso-hansi-theme-font:minor-latin'>OK</span></b>. To see the plot select the 1597 PWDR line in the GSAS-II data tree (Ive adjusted the width and height).</p> 1598 1599 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1600 id="_x0000_i1046" type="#_x0000_t75" style='width:525pt;height:450pt; 1601 visibility:visible;mso-wrap-style:square'> 1589 1602 <v:imagedata src="MCSA%20in%20GSAS_files/image011.png" o:title=""/> 1590 1603 </v:shape><![endif]--><![if !vml]><img border=0 width=700 height=600 1591 src="MCSA%20in%20GSAS_files/image012.gif" v:shapes=" Picture_x0020_6"><![endif]></span></p>1604 src="MCSA%20in%20GSAS_files/image012.gif" v:shapes="_x0000_i1046"><![endif]></span></p> 1592 1605 1593 1606 <p class=MsoNormal><o:p> </o:p></p> … … 1603 1616 to <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1604 1617 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>1</span></b>, 1605 then change <b style='mso-bidi-font-weight:normal'><span style='font-family: 1618 then change <span class=SpellE><b style='mso-bidi-font-weight:normal'><span 1619 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1620 mso-hansi-theme-font:minor-latin'>pos</span></b></span> to <b style='mso-bidi-font-weight: 1621 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1622 minor-latin;mso-hansi-theme-font:minor-latin'>19</span></b>, <span 1623 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 1606 1624 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1607 minor-latin'>pos</span></b> to <b style='mso-bidi-font-weight:normal'><span 1608 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1609 mso-hansi-theme-font:minor-latin'>19</span></b>, <b style='mso-bidi-font-weight: 1610 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1611 minor-latin;mso-hansi-theme-font:minor-latin'>int</span></b> to <b 1612 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1613 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>100000</span></b> 1614 and <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1615 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>sig</span></b> 1616 to <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1625 minor-latin'>int</span></b></span> to <b style='mso-bidi-font-weight:normal'><span 1626 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1627 mso-hansi-theme-font:minor-latin'>100000</span></b> and <b style='mso-bidi-font-weight: 1628 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1629 minor-latin;mso-hansi-theme-font:minor-latin'>sig</span></b> to <b 1630 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1617 1631 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>100000</span></b>; 1618 1632 for the peak only refine <b style='mso-bidi-font-weight:normal'><span … … 1620 1634 mso-hansi-theme-font:minor-latin'>int</span></b>. I also set <b 1621 1635 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1622 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>No. coeff.</span></b> 1623 to <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1636 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>No. <span 1637 class=SpellE>coeff</span>.</span></b> <span class=GramE>to</span> <b 1638 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1624 1639 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>5</span></b>; 1625 1640 the window shows</p> 1626 1641 1627 1642 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1628 id="Picture_x0020_ 7" o:spid="_x0000_i1045" type="#_x0000_t75" style='width:402pt;1643 id="Picture_x0020_33" o:spid="_x0000_i1045" type="#_x0000_t75" style='width:402pt; 1629 1644 height:197.25pt;visibility:visible;mso-wrap-style:square'> 1630 1645 <v:imagedata src="MCSA%20in%20GSAS_files/image013.png" o:title=""/> 1631 1646 </v:shape><![endif]--><![if !vml]><img border=0 width=536 height=263 1632 src="MCSA%20in%20GSAS_files/image014.gif" v:shapes="Picture_x0020_ 7"><![endif]></span></p>1647 src="MCSA%20in%20GSAS_files/image014.gif" v:shapes="Picture_x0020_33"><![endif]></span></p> 1633 1648 1634 1649 <p class=MsoNormal><o:p> </o:p></p> … … 1636 1651 <p class=MsoNormal>Do a <b style='mso-bidi-font-weight:normal'><span 1637 1652 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1638 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; my Rwp fell to 1639 ~8.5%. In <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1640 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Background</span></b> 1641 select <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1642 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>pos</span></b> 1643 & <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1644 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>sig</span></b> 1645 for refinement and in <b style='mso-bidi-font-weight:normal'><span 1646 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1647 mso-hansi-theme-font:minor-latin'>Instrument Parameters</span></b> select <b 1648 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1649 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>X</span></b> 1650 & <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1651 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Y</span></b> 1652 for refinement; then do <b style='mso-bidi-font-weight:normal'><span 1653 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1654 mso-hansi-theme-font:minor-latin'>Calculate /Refine</span></b> again. I got 1655 Rwp~7.5% at convergence and the fit looked like</p> 1653 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; my <span 1654 class=SpellE>Rwp</span> fell to ~8.5%. In <b style='mso-bidi-font-weight:normal'><span 1655 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1656 mso-hansi-theme-font:minor-latin'>Background</span></b> select <span 1657 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 1658 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1659 minor-latin'>pos</span></b></span> & <b style='mso-bidi-font-weight:normal'><span 1660 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1661 mso-hansi-theme-font:minor-latin'>sig</span></b> for refinement and in <b 1662 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1663 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Instrument 1664 Parameters</span></b> select <b style='mso-bidi-font-weight:normal'><span 1665 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1666 mso-hansi-theme-font:minor-latin'>X</span></b> & <b style='mso-bidi-font-weight: 1667 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1668 minor-latin;mso-hansi-theme-font:minor-latin'>Y</span></b> for refinement; then 1669 do <span class=GramE><b style='mso-bidi-font-weight:normal'><span 1670 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1671 mso-hansi-theme-font:minor-latin'>Calculate</span></b></span><b 1672 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1673 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'> /Refine</span></b> 1674 again. I got Rwp~7.5% at convergence and the fit looked like</p> 1656 1675 1657 1676 <p class=MsoNormal><span style='mso-no-proof:yes'>c<!--[if gte vml 1]><v:shape 1658 id=" Picture_x0020_33" o:spid="_x0000_i1044" type="#_x0000_t75" style='width:525pt;1659 height:450pt;visibility:visible;mso-wrap-style:square'>1677 id="_x0000_i1044" type="#_x0000_t75" style='width:525pt;height:450pt; 1678 visibility:visible;mso-wrap-style:square'> 1660 1679 <v:imagedata src="MCSA%20in%20GSAS_files/image015.png" o:title=""/> 1661 1680 </v:shape><![endif]--><![if !vml]><img border=0 width=700 height=600 1662 src="MCSA%20in%20GSAS_files/image016.gif" v:shapes=" Picture_x0020_33"><![endif]></span></p>1681 src="MCSA%20in%20GSAS_files/image016.gif" v:shapes="_x0000_i1044"><![endif]></span></p> 1663 1682 1664 1683 <p class=MsoNormal><o:p> </o:p></p> … … 1695 1714 mso-hansi-theme-font:minor-latin'>Residue rigid bodies</span></b> and then do <b 1696 1715 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1697 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Edit /Import1716 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Edit Body/Import 1698 1717 XYZ</span></b>. Change the file type to <b style='mso-bidi-font-weight:normal'><span 1699 1718 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1700 mso-hansi-theme-font:minor-latin'>XYZ file (*.xyz)</span></b> and select <b 1701 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1702 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>aminoquinoline.xyz</span></b>. 1703 The <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1704 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Rigid bodies</span></b> 1719 mso-hansi-theme-font:minor-latin'>XYZ file (*.xyz)</span></b> and select <span 1720 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 1721 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1722 minor-latin'>aminoquinoline.xyz</span></b></span>. The <b style='mso-bidi-font-weight: 1723 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1724 minor-latin;mso-hansi-theme-font:minor-latin'>Rigid <span class=GramE>bodies</span></span></b> 1705 1725 window will show</p> 1706 1726 … … 1716 1736 <p class=MsoNormal>If you select the <b style='mso-bidi-font-weight:normal'><span 1717 1737 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1718 mso-hansi-theme-font:minor-latin'>Plot?</span></b> box a drawing will appear</p> 1719 1720 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1721 id="Picture_x0020_1" o:spid="_x0000_i1041" type="#_x0000_t75" style='width:286.5pt; 1722 height:234pt;visibility:visible;mso-wrap-style:square'> 1738 mso-hansi-theme-font:minor-latin'>Plot?</span></b> <span class=GramE>box</span> 1739 a drawing will appear</p> 1740 1741 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1742 id="_x0000_i1041" type="#_x0000_t75" style='width:286.5pt;height:234pt; 1743 visibility:visible;mso-wrap-style:square'> 1723 1744 <v:imagedata src="MCSA%20in%20GSAS_files/image021.png" o:title="" croptop="13685f" 1724 1745 cropbottom="15961f" cropleft="36898f" cropright="13504f"/> 1725 1746 </v:shape><![endif]--><![if !vml]><img border=0 width=382 height=312 1726 src="MCSA%20in%20GSAS_files/image022.jpg" v:shapes=" Picture_x0020_1"><![endif]></span></p>1747 src="MCSA%20in%20GSAS_files/image022.jpg" v:shapes="_x0000_i1041"><![endif]></span></p> 1727 1748 1728 1749 <p class=MsoNormal><o:p> </o:p></p> … … 1788 1809 <p class=MsoNormal>You see the four equivalent molecules for space group P2<sub>1</sub>2<sub>1</sub>2<sub>1</sub> 1789 1810 whose centers fall within the unit cell; the reference one is that in the lower 1790 left corner with C1 at the origin. You can enter values in the <b 1791 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1792 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>x,y,z</span></b> 1793 boxes and see how the structure changes. Similarly, you can enter an angle in <b 1794 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1795 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Oa</span></b> 1796 or a vector in the <b style='mso-bidi-font-weight:normal'><span 1797 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1798 mso-hansi-theme-font:minor-latin'>Oi,Oj,Ok</span></b><span 1811 left corner with C1 at the origin. You can enter values in the <span 1812 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 1813 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1814 minor-latin'>x<span class=GramE>,y,z</span></span></b></span> boxes and see how 1815 the structure changes. Similarly, you can enter an angle in <span class=SpellE><b 1816 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1817 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Oa</span></b></span> 1818 or a vector in the <span class=SpellE><b style='mso-bidi-font-weight:normal'><span 1819 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1820 mso-hansi-theme-font:minor-latin'>Oi<span class=GramE>,Oj,Ok</span></span></b></span><span 1799 1821 style='mso-spacerun:yes'> </span>box to see the effects of rotating the 1800 1822 molecule. Each of these parameters has a defined range for the MC/SA … … 1813 1835 boxes and select <b style='mso-bidi-font-weight:normal'><span style='font-family: 1814 1836 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1815 minor-latin'>AV</span></b> from the <b style='mso-bidi-font-weight:normal'><span 1816 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1817 mso-hansi-theme-font:minor-latin'>Oa</span></b>, <b style='mso-bidi-font-weight: 1818 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1819 minor-latin;mso-hansi-theme-font:minor-latin'>Oi,Oj,Ok Vary?</span></b> 1820 pulldown. The MC/SA window should look like</p> 1821 1822 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1823 id="Picture_x0020_39" o:spid="_x0000_i1037" type="#_x0000_t75" style='width:448.5pt; 1824 height:193.5pt;visibility:visible;mso-wrap-style:square'> 1837 minor-latin'>AV</span></b> from the <span class=SpellE><b style='mso-bidi-font-weight: 1838 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1839 minor-latin;mso-hansi-theme-font:minor-latin'>Oa</span></b></span>, <span 1840 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 1841 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1842 minor-latin'>Oi<span class=GramE>,Oj,Ok</span></span></b></span><b 1843 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1844 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'> Vary?</span></b> 1845 <span class=GramE>pulldown</span>. The MC/SA window should look like</p> 1846 1847 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1848 id="_x0000_i1037" type="#_x0000_t75" style='width:448.5pt;height:193.5pt; 1849 visibility:visible;mso-wrap-style:square'> 1825 1850 <v:imagedata src="MCSA%20in%20GSAS_files/image029.png" o:title=""/> 1826 1851 </v:shape><![endif]--><![if !vml]><img border=0 width=598 height=258 1827 src="MCSA%20in%20GSAS_files/image030.gif" v:shapes=" Picture_x0020_39"><![endif]></span><span1852 src="MCSA%20in%20GSAS_files/image030.gif" v:shapes="_x0000_i1037"><![endif]></span><span 1828 1853 style='mso-spacerun:yes'> </span></p> 1829 1854 … … 1845 1870 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1846 1871 mso-hansi-theme-font:minor-latin'>Pawley reflections</span></b> then the 1847 structure factors will be Pawley model values, Fsq(hkl) in Pawley table, and 1848 the covariance matrix developed from the refinement will be used for peak 1849 overlap effects in the MC/SA calculations. Alternatively, if you choose <b 1850 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1851 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>PWDR 1852 Quinoli.gda Bank 1 </span></b>then the structure factors will be those, Fosq in 1872 structure factors will be Pawley model values, <span class=SpellE><span 1873 class=GramE>Fsq</span></span><span class=GramE>(</span><span class=SpellE>hkl</span>) 1874 in Pawley table, and the covariance matrix developed from the refinement will 1875 be used for peak overlap effects in the MC/SA calculations. Alternatively, if 1876 you choose <b style='mso-bidi-font-weight:normal'><span style='font-family: 1877 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1878 minor-latin'>PWDR <span class=SpellE>Quinoli.gda</span> Bank 1 </span></b>then 1879 the structure factors will be those, <span class=SpellE>Fosq</span> in 1853 1880 Reflection table, extracted during the Pawley refinement and the peak FWHMs 1854 1881 will be used for peak overlap effects in the MC/SA calculations. Notice that 1855 the Pawley Fsq(hkl) values are the same as the Fcsq Reflection List values. 1856 Here I chose <b style='mso-bidi-font-weight:normal'><span style='font-family: 1857 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1858 minor-latin'>Pawley reflections</span></b> for the <b style='mso-bidi-font-weight: 1882 the Pawley <span class=SpellE><span class=GramE>Fsq</span></span><span 1883 class=GramE>(</span><span class=SpellE>hkl</span>) values are the same as the <span 1884 class=SpellE>Fcsq</span> Reflection List values. Here I chose <b 1885 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1886 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>PWDR <span 1887 class=SpellE>Quinoli.gda</span> Bank 1 </span></b>for the <b style='mso-bidi-font-weight: 1859 1888 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 1860 1889 minor-latin;mso-hansi-theme-font:minor-latin'>Reflection set</span></b> as it … … 1886 1915 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>No. trials</span></b>; 1887 1916 this is the number of random tries at each temperature in the annealing 1888 schedule. When done the <b style='mso-bidi-font-weight:normal'><span 1889 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1890 mso-hansi-theme-font:minor-latin'>General</span></b> window should look like</p> 1891 1892 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1893 id="Picture_x0020_16" o:spid="_x0000_i1036" type="#_x0000_t75" style='width:7in; 1917 schedule. I also used <b style='mso-bidi-font-weight:normal'><span 1918 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1919 mso-hansi-theme-font:minor-latin'>4 MC/SA runs</span></b> <span class=GramE>When</span> 1920 done the <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1921 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>General</span></b> 1922 window should look like</p> 1923 1924 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1925 id="Picture_x0020_38" o:spid="_x0000_i1036" type="#_x0000_t75" style='width:7in; 1894 1926 height:375pt;visibility:visible;mso-wrap-style:square'> 1895 1927 <v:imagedata src="MCSA%20in%20GSAS_files/image031.png" o:title=""/> 1896 1928 </v:shape><![endif]--><![if !vml]><img border=0 width=672 height=500 1897 src="MCSA%20in%20GSAS_files/image032.gif" v:shapes="Picture_x0020_ 16"><![endif]></span></p>1929 src="MCSA%20in%20GSAS_files/image032.gif" v:shapes="Picture_x0020_38"><![endif]></span></p> 1898 1930 1899 1931 <p class=MsoNormal><o:p> </o:p></p> … … 1913 1945 1914 1946 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1915 id="Picture_x0020_ 19" o:spid="_x0000_i1035" type="#_x0000_t75" style='width:560.25pt;1947 id="Picture_x0020_39" o:spid="_x0000_i1035" type="#_x0000_t75" style='width:564pt; 1916 1948 height:273.75pt;visibility:visible;mso-wrap-style:square'> 1917 1949 <v:imagedata src="MCSA%20in%20GSAS_files/image033.png" o:title=""/> 1918 </v:shape><![endif]--><![if !vml]><img border=0 width=7 47height=3651919 src="MCSA%20in%20GSAS_files/image034.gif" v:shapes="Picture_x0020_ 19"><![endif]></span></p>1950 </v:shape><![endif]--><![if !vml]><img border=0 width=752 height=365 1951 src="MCSA%20in%20GSAS_files/image034.gif" v:shapes="Picture_x0020_39"><![endif]></span></p> 1920 1952 1921 1953 <p class=MsoNormal><o:p> </o:p></p> … … 1925 1957 1926 1958 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1927 id="Picture_x0020_ 20" o:spid="_x0000_i1034" type="#_x0000_t75" style='width:279.75pt;1928 height:2 43.75pt;visibility:visible;mso-wrap-style:square'>1929 <v:imagedata src="MCSA%20in%20GSAS_files/image035.png" o:title="" croptop=" 11773f"1930 cropbottom=" 16530f" cropleft="37715f" cropright="13079f"/>1931 </v:shape><![endif]--><![if !vml]><img border=0 width= 373 height=3251932 src="MCSA%20in%20GSAS_files/image036.jpg" v:shapes="Picture_x0020_ 20"><![endif]></span></p>1959 id="Picture_x0020_1" o:spid="_x0000_i1034" type="#_x0000_t75" style='width:364.5pt; 1960 height:298.5pt;visibility:visible;mso-wrap-style:square'> 1961 <v:imagedata src="MCSA%20in%20GSAS_files/image035.png" o:title="" croptop="8643f" 1962 cropbottom="21002f" cropleft="39996f" cropright="10456f"/> 1963 </v:shape><![endif]--><![if !vml]><img border=0 width=486 height=398 1964 src="MCSA%20in%20GSAS_files/image036.jpg" v:shapes="Picture_x0020_1"><![endif]></span></p> 1933 1965 1934 1966 <p class=MsoNormal><o:p> </o:p></p> … … 1936 1968 <p class=MsoNormal>If you are lucky (like I was!) then the result is clearly a 1937 1969 good solution with a very low Residual (~4.8%). This problem typically gives 1938 suitable solutions with Residuals of 4-5%. If no good solution appears (e.g. 1939 molecules clashing), then you should just rerun MC/SA perhaps using more runs 1940 or more trials. Be sure to set the <b style='mso-bidi-font-weight:normal'><span 1941 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1942 mso-hansi-theme-font:minor-latin'>Keep</span></b> box for any solutions you 1943 want to retain; the others will be cleared before the next MC/SA run starts. 1944 When you think you have a good one, <b style='mso-bidi-font-weight:normal'><span 1945 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1946 mso-hansi-theme-font:minor-latin'>Select</span></b> it; the parameters will be 1947 copied to the appropriate boxes in the upper part of the MC/SA window and the 1948 structure is drawn again. Next we want to refine this solution so go back to 1949 the <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1950 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>General</span></b> 1951 tab.</p> 1970 suitable solutions with Residuals of 4-5%. If your cell indexing resulted in <span 1971 class=GramE>a &</span> b swapped, you should view the structure along b 1972 (green). If no good solution appears (e.g. molecules clashing), then you should 1973 just rerun MC/SA perhaps using more runs or more trials. Be sure to set the <b 1974 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1975 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Keep</span></b> 1976 box for any solutions you want to retain; the others will be cleared before the 1977 next MC/SA run starts. When you think you have a good one, <span class=GramE><b 1978 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 1979 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Select</span></b></span> 1980 it; the parameters will be copied to the appropriate boxes in the upper part of 1981 the MC/SA window and the structure is drawn again. Next we want to refine this 1982 solution so go back to the <b style='mso-bidi-font-weight:normal'><span 1983 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 1984 mso-hansi-theme-font:minor-latin'>General</span></b> tab.</p> 1952 1985 1953 1986 <p class=MsoNormal>MC/SA refinement is achieved by narrowing the search ranges … … 1964 1997 1965 1998 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape 1966 id="Picture_x0020_ 22" o:spid="_x0000_i1033" type="#_x0000_t75" style='width:569.25pt;1967 height:2 85pt;visibility:visible;mso-wrap-style:square'>1999 id="Picture_x0020_40" o:spid="_x0000_i1033" type="#_x0000_t75" style='width:448.5pt; 2000 height:273.75pt;visibility:visible;mso-wrap-style:square'> 1968 2001 <v:imagedata src="MCSA%20in%20GSAS_files/image037.png" o:title=""/> 1969 </v:shape><![endif]--><![if !vml]><img border=0 width= 759 height=3801970 src="MCSA%20in%20GSAS_files/image038.gif" v:shapes="Picture_x0020_ 22"><![endif]></span></p>2002 </v:shape><![endif]--><![if !vml]><img border=0 width=598 height=365 2003 src="MCSA%20in%20GSAS_files/image038.gif" v:shapes="Picture_x0020_40"><![endif]></span></p> 1971 2004 1972 2005 <p class=MsoNormal><o:p> </o:p></p> 1973 2006 1974 <p class=MsoNormal>This is clearly a good result (R ~ 0.4%!) and the structure</p> 2007 <p class=MsoNormal>This is clearly a good result (R ~ <span class=GramE>2.6%</span>!) 2008 and the structure</p> 1975 2009 1976 2010 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 1984 2018 <p class=MsoNormal><o:p> </o:p></p> 1985 2019 1986 <p class=MsoNormal> isidentical to the published one</p>2020 <p class=MsoNormal><span class=GramE>is</span> identical to the published one</p> 1987 2021 1988 2022 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 1994 2028 src="MCSA%20in%20GSAS_files/image042.jpg" v:shapes="Picture_x0020_24"><![endif]></span></p> 1995 2029 1996 <p class=MsoNormal> asdrawn by the Mercury program </p>2030 <p class=MsoNormal><span class=GramE>as</span> drawn by the Mercury program </p> 1997 2031 1998 2032 <p class=MsoNormal>(<a … … 2049 2083 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 2050 2084 minor-latin;mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b> from 2051 the main GSAS-II data tree window; my R<sub>wp</sub> was ~9%. Next check the <b 2052 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2053 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Refine unit 2054 cell</span></b> in the <b style='mso-bidi-font-weight:normal'><span 2055 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2056 mso-hansi-theme-font:minor-latin'>General</span></b> window and in the <b 2057 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2058 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Data</span></b> 2059 window, using the <b style='mso-bidi-font-weight:normal'><span 2060 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2061 mso-hansi-theme-font:minor-latin'>uniaxial</span></b> model for mustrain, check 2062 both <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2063 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Cryst. size</span></b> 2064 & <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2065 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>mustrain</span></b> 2085 the main GSAS-II data tree window; my <span class=SpellE>R<sub>wp</sub></span> 2086 was ~9%. Next check the <b style='mso-bidi-font-weight:normal'><span 2087 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2088 mso-hansi-theme-font:minor-latin'>Refine unit cell</span></b> in the <b 2089 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2090 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>General</span></b> 2091 window and in the <b style='mso-bidi-font-weight:normal'><span 2092 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2093 mso-hansi-theme-font:minor-latin'>Data</span></b> window, using the <b 2094 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2095 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>uniaxial</span></b> 2096 model for <span class=SpellE>mustrain</span>, check both <span class=SpellE><b 2097 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2098 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Cryst</span></b></span><b 2099 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2100 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>. <span 2101 class=GramE>size</span></span></b> & <span class=SpellE><b 2102 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2103 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>mustrain</span></b></span> 2066 2104 boxes and do another <b style='mso-bidi-font-weight:normal'><span 2067 2105 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2068 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; I got an R<sub>wp</sub> 2069 ~8.5% for this. Weve not refined atom positions or thermal parameters; to do 2070 this we need to describe the structure as a rigid body.</p> 2106 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; I got an <span 2107 class=SpellE>R<sub>wp</sub></span> ~8.5% for this. Weve not refined atom 2108 positions or thermal parameters; to do this we need to describe the structure 2109 as a rigid body.</p> 2071 2110 2072 2111 <h2>Step 9. Rigid body model</h2> … … 2115 2154 corresponding atoms in the structure (C(2), C(3) & C(6)). Do these in 2116 2155 order; the drawing will show the shift of molecular origin and then its 2117 reorientation as the three atoms are selected. Press < b style='mso-bidi-font-weight:2118 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 2119 m inor-latin;mso-hansi-theme-font:minor-latin'>Ready</span></b> when done; the2120 structure will be drawn with yellow bonds indicating that it is now a rigid 2121 body and not independent atoms. The <b style='mso-bidi-font-weight:normal'><span 2122 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2123 m so-hansi-theme-font:minor-latin'>RB Models</span></b> window shows the new2124 rigid body</p>2156 reorientation as the three atoms are selected. Press <span class=GramE><b 2157 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2158 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Ready</span></b></span> 2159 when done; the structure will be drawn with yellow bonds indicating that it is 2160 now a rigid body and not independent atoms. The <b style='mso-bidi-font-weight: 2161 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 2162 minor-latin;mso-hansi-theme-font:minor-latin'>RB Models</span></b> window shows 2163 the new rigid body</p> 2125 2164 2126 2165 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2134 2173 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2135 2174 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Refine?</span></b> 2136 box, select <b style='mso-bidi-font-weight:normal'><span style='font-family: 2175 <span class=GramE>box</span>, select <b style='mso-bidi-font-weight:normal'><span 2176 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2177 mso-hansi-theme-font:minor-latin'>AV</span></b> from the pull down, choose <span 2178 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 2137 2179 "Calibri",sans-serif;mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 2138 minor-latin'>AV</span></b> from the pull down, choose <b style='mso-bidi-font-weight: 2139 normal'><span style='font-family:"Calibri",sans-serif;mso-ascii-theme-font: 2140 minor-latin;mso-hansi-theme-font:minor-latin'>Uiso</span></b> from the thermal 2141 parameter model and check its <b style='mso-bidi-font-weight:normal'><span 2142 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2143 mso-hansi-theme-font:minor-latin'>Refine?</span></b> box. Then do <b 2144 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2145 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; 2146 after a couple of refinement runs I got an R<sub>wp</sub> ~7.9% and the <b 2147 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2148 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>RBModels</span></b> 2149 window shows the new parameters</p> 2180 minor-latin'>Uiso</span></b></span> from the thermal parameter model and check 2181 its <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri",sans-serif; 2182 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Refine?</span></b> 2183 <span class=GramE>box</span>. Then do <b style='mso-bidi-font-weight:normal'><span 2184 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2185 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; after a couple 2186 of refinement runs I got an <span class=SpellE>R<sub>wp</sub></span> ~7.9% and 2187 the <span class=SpellE><b style='mso-bidi-font-weight:normal'><span 2188 style='font-family:"Calibri",sans-serif;mso-ascii-theme-font:minor-latin; 2189 mso-hansi-theme-font:minor-latin'>RBModels</span></b></span> window shows the 2190 new parameters</p> 2150 2191 2151 2192 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2156 2197 src="MCSA%20in%20GSAS_files/image052.gif" v:shapes="Picture_x0020_32"><![endif]></span></p> 2157 2198 2158 <p class=MsoNormal> andthe profile shows the fit</p>2199 <p class=MsoNormal><span class=GramE>and</span> the profile shows the fit</p> 2159 2200 2160 2201 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2167 2208 <p class=MsoNormal>This completes the structure analysis for 3-aminoquinoline. 2168 2209 Next is a more complex case involving some internal degrees of freedom in the 2169 molecule of a-d- lactose.</p>2210 molecule of a-d-<span class=GramE>lactose.</span></p> 2170 2211 2171 2212 </div>
Note: See TracChangeset
for help on using the changeset viewer.