Changeset 872
- Timestamp:
- Mar 22, 2013 10:24:31 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/help/Charge Flipping - sucrose.htm
r870 r872 24 24 <o:Author>Von Dreele</o:Author> 25 25 <o:LastAuthor>Von Dreele</o:LastAuthor> 26 <o:Revision>1 3</o:Revision>27 <o:TotalTime>1 342</o:TotalTime>26 <o:Revision>14</o:Revision> 27 <o:TotalTime>1423</o:TotalTime> 28 28 <o:Created>2013-03-15T18:32:00Z</o:Created> 29 <o:LastSaved>2013-03-2 1T16:12:00Z</o:LastSaved>29 <o:LastSaved>2013-03-22T15:24:00Z</o:LastSaved> 30 30 <o:Pages>17</o:Pages> 31 <o:Words>33 95</o:Words>32 <o:Characters>19 358</o:Characters>31 <o:Words>3367</o:Words> 32 <o:Characters>19197</o:Characters> 33 33 <o:Company>Argonne National Laboratory</o:Company> 34 <o:Lines>1 61</o:Lines>34 <o:Lines>159</o:Lines> 35 35 <o:Paragraphs>45</o:Paragraphs> 36 <o:CharactersWithSpaces>22 708</o:CharactersWithSpaces>36 <o:CharactersWithSpaces>22519</o:CharactersWithSpaces> 37 37 <o:Version>14.00</o:Version> 38 38 </o:DocumentProperties> … … 953 953 fitting-indexing exercise</a>. If you didnt do that exercise, do it now.</p> 954 954 955 <h2> <span class=GramE>Step 1.</span>Setup for Pawley refinement</h2>955 <h2>Step 1: Setup for Pawley refinement</h2> 956 956 957 957 <p class=MsoNormal>There are a number of steps that must be done in preparation … … 969 969 Symbol;mso-ascii-font-family:"Times New Roman";mso-hansi-font-family:"Times New Roman"; 970 970 mso-char-type:symbol;mso-symbol-font-family:Symbol'><span style='mso-char-type: 971 symbol;mso-symbol-font-family:Symbol'>Q</span></span>. It is an area relatively 972 clear of peaks and makes a good place to end the calculations with a reflection 973 set of ~1Å resolution as needed for charge flipping. Enter <b style='mso-bidi-font-weight: 974 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 975 minor-latin;mso-hansi-theme-font:minor-latin'>23.55</span></b> for <span 976 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 977 "Calibri","sans-serif";mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 978 minor-latin'>Tmax</span></b></span> as the upper limit and make a note of the 979 exact d-spacing (1.0125Å) using the mouse cursor on the plot.</p> 971 symbol;mso-symbol-font-family:Symbol'>Q</span></span>. There is an area 972 relatively clear of peaks that makes a good place to end the calculations with 973 a reflection set of <span class=SpellE>d<sub>min</sub></span> ~1Å resolution needed 974 for charge flipping. Enter <b style='mso-bidi-font-weight:normal'><span 975 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 976 mso-hansi-theme-font:minor-latin'>23.55</span></b> for <span class=SpellE><b 977 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 978 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Tmax</span></b></span> 979 as the upper limit and make a note of the exact d-spacing (1.0125Å) using the 980 mouse cursor on the plot.</p> 980 981 981 982 <p class=MsoListParagraph style='text-indent:-.25in;mso-list:l6 level1 lfo2'><![if !supportLists]><span … … 1021 1022 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 1022 1023 minor-latin;mso-hansi-theme-font:minor-latin'>Do Pawley refinement</span></b> 1023 box and enter the d-spacing (1.0125) into the <b style='mso-bidi-font-weight: 1024 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 1025 minor-latin;mso-hansi-theme-font:minor-latin'>Pawley <span class=SpellE>dmin</span></span></b> 1026 box.</p> 1024 box and enter the d-spacing (<b style='mso-bidi-font-weight:normal'><span 1025 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1026 mso-hansi-theme-font:minor-latin'>1.0125</span></b>) into the <b 1027 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1028 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Pawley <span 1029 class=SpellE>dmin</span></span></b> box.</p> 1027 1030 1028 1031 <p class=MsoNormal>When done the <b style='mso-bidi-font-weight:normal'><span … … 1093 1096 surfaces</p> 1094 1097 1095 <p class=MsoNormal><span class=GramE>a nd</span>preferred orientation1098 <p class=MsoNormal><span class=GramE>as</span> well as preferred orientation 1096 1099 correction curves.</p> 1097 1100 … … 1119 1122 you didnt skip a step as it might not work correctly if you did.</p> 1120 1123 1121 <h2>Step 2 Pawley refinement</h2>1124 <h2>Step 2: Pawley refinement</h2> 1122 1125 1123 1126 <p class=MsoNormal>To start the refinement, do <b style='mso-bidi-font-weight: … … 1197 1200 1198 1201 <p class=MsoNormal>What you see is that the background has a substantial 1199 diffuse scattering feature at ~5deg as well as some peak shape issues (which 1200 may be a consequence of the poor background fit). Also notice the negative <b 1201 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1202 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Y</span></b> 1203 value in <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1204 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Instrument 1205 parameters</span></b>; the problem is that the default <span class=SpellE>mustrain</span> 1206 value (1000.0) is much too high for this material and wed like the fit to be 1202 diffuse scattering feature at ~5deg and there are some peak shape issues (which 1203 may be a consequence of the poor background fit). Wed like the fit to be 1207 1204 better before trying charge flipping. First select <span class=GramE><b 1208 1205 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; … … 1263 1260 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 1264 1261 minor-latin;mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b> until 1265 the refinement converges. After a set of refinement cycles go to the Pawley reflections 1262 the refinement converges. After a set of refinement cycles go to the <b 1263 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1264 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Pawley reflections</span></b> 1266 1265 tab and check the intensities. If there are many negatives amongst the highest 1267 1266 angle reflections do <b style='mso-bidi-font-weight:normal'><span 1268 1267 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1269 1268 mso-hansi-theme-font:minor-latin'>Pawley update</span></b>, this replaces all 1270 the values with those in the extracted from the profile a ndlisted in the <b1269 the values with those in the extracted from the profile as listed in the <b 1271 1270 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1272 1271 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Reflection 1273 list</span></b> for the PWDR data set. Next in the <b style='mso-bidi-font-weight: 1274 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 1275 minor-latin;mso-hansi-theme-font:minor-latin'>General</span></b> tab set the <b 1272 list</span></b> for the PWDR data set. This makes all the negative intensities 1273 positive. Next in the <b style='mso-bidi-font-weight:normal'><span 1274 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1275 mso-hansi-theme-font:minor-latin'>General</span></b> tab set the <b 1276 1276 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1277 1277 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Pawley neg. … … 1293 1293 <p class=MsoNormal>Lets see if this is good enough for charge flipping.</p> 1294 1294 1295 <h2>Step 3 Charge flipping</h2>1295 <h2>Step 3: Charge flipping</h2> 1296 1296 1297 1297 <p class=MsoNormal>A few steps are needed to set up for a charge flipping run. … … 1317 1317 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1318 1318 mso-hansi-theme-font:minor-latin'>Reflection Lists</span></b> under the PWDR 1319 data set; these should match the Pawley reflection set.</p> 1319 data set; these should match the <b style='mso-bidi-font-weight:normal'><span 1320 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1321 mso-hansi-theme-font:minor-latin'>Pawley reflection</span></b> set.</p> 1320 1322 1321 1323 <p class=MsoListParagraph style='margin-left:.75in;text-indent:-.25in; … … 1351 1353 "Calibri","sans-serif";mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 1352 1354 minor-latin'>10-12</span></b> for equal atom problems and larger for heavy atom 1353 ones. <span class=GramE>For sucrose, set this to <b style='mso-bidi-font-weight:1354 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 1355 m inor-latin;mso-hansi-theme-font:minor-latin'>12.</span></b></span> When ready1356 s elect <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif";1355 ones. For sucrose set it to <b style='mso-bidi-font-weight:normal'><span 1356 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1357 mso-hansi-theme-font:minor-latin'>12.</span></b> When ready select <b 1358 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1357 1359 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Compute/Charge 1358 1360 flipping</span></b>; a progress bar dialog will appear that tracks the residual 1359 1361 found in each charge flipping cycle. It will quickly drop to some level 1360 (42-44%) and then perhaps drop again into the 20-25% range after a considerable1361 length of time for a successful charge flipping run. You should save the1362 (42-44%) and then for a successful charge flipping run perhaps drop again into 1363 the 20-25% range after a considerable length of time. You should save the 1362 1364 project if you get an apparently good solution. Many times it will not drop 1363 1365 below ~40%; just try <b style='mso-bidi-font-weight:normal'><span … … 1401 1403 22YiTra2kYMu1l1tQD30/bPm3wwYN0x18gb45AdQl1tp5j/sFB2T0FQ7R0nTNEV3j6o9feQzro1i 1402 1404 OWA14Fm+Q8a1a8+Bvu/d/dMb2JY5uiPbhG/ktn4cqGU/er3pcvwCAAD//wMAUEsDBBQABgAIAAAA 1403 IQ COfwQ4JQQAANgQAAAOAAAAZHJzL2Uyb0RvYy54bWzsWFuP4jYUfq/U/2DlsRKThIQA0TCrKZdR1405 IQDQSJmzIgQAANIQAAAOAAAAZHJzL2Uyb0RvYy54bWzsWFuP4jYUfq/U/2DlsRKThIQA0TCrKZdR 1404 1406 pW270tIfYBznok3s1DaE2VX/e48vgcAUtdqVqqqCB7Dj43P5zvH5HB7fHZsaHaiQFWcLL3wIPEQZ 1405 1407 4VnFioX323YzmnlIKswyXHNGF94rld67p++/e+zalI55yeuMCgRKmEy7duGVSrWp70tS0gbLB95S … … 1410 1412 Rcq12xdOomhit4WRHmpfcaodcE47J3WEUHnyDK78NnA/lrilJmdSA9eDO+7BfRYERVMLrZHocZVD 1411 1413 UAcr2j8J2P81nDojFwDG4Ty0GPUYzsM4tlCY0RAJnJK9VC+UAyo4xYf3UsEyVHEGIztwNbGFJORN 1412 DWfjBx /FCeqQVhb0ByjrxaCITmIBKm9IARonqXl4U1k0EAObt7RBdCdtAUoC+KBkMokSi0Nxcg4q1413 4rYglMcpbFz2SJAjc1DACGHd1ba2bFsudQFqYKA4twZ00AFiGriz9MQUeS8NkWvpyBWkkba7nBkB 1414 Leu6WQkPQbPa2XBarLR32ogeom7huVSg8jyUSvBPcPYCY73hB7rlZovSTseTKTRP44rJIHhwFiH7 1415 XUV+pJ+HG6bBfBoAfOcT2hpFvWV4Pg7iZO6O55tFZ9DW3qWBmg0NDTy72NNL9b/WwkC6j6MXIDWX 1416 1NrTMP0NXHlVA0EwoKj/OF6QKRuOG5hK0AU0OMeMbyAcc5BrpusjnAcTWwiS11WmVzUeUhS7ZS3Q 1417 AUPsm40+N64sL8SAhFhmtJUUZ2s3Vriq7Ris16bmoeG4wtStx7Dfl3kwX8/Ws3gUj5P1KA5Wq9Hz 1418 Z hmPkk04nayi1XK5Cv/QmIdxWlZZRpn2rmfiMP5nzdjdCSyHnrj4IoqrYDfweRusf+mG4QyIpf+11419 YPcN2VLHjmev0JwNq8GJgqsQ0FvJxWcPdXCtWHjy9z0W1EP1TwyoxRwYuIeYiSlfD4nhym64ghkB 1420 VQuPKGgBdrJU9vayb0VVlGArNBXL+DOwbF7p/g2kLFPrl5sAv/1bRAc9294iDNHNNMbaHaDCbyM6 1421 i3AYzHSZ2mJ2vJdoxht0pjvvWaK4896A1G7yxZ337rx3573/M++Zd1R4fTYs7l719fv5cG548vyH 1422 xNOfAAAA//8DAFBLAwQUAAYACAAAACEAVz1x3t8AAAALAQAADwAAAGRycy9kb3ducmV2LnhtbEyP 1423 zU7DMBCE70i8g7VI3Kjz1zakcSqE6BlRuHBzYzexiNep7TaBp2c50duM9tPsTL2d7cAu2gfjUEC6 1424 SIBpbJ0y2An4eN89lMBClKjk4FAL+NYBts3tTS0r5SZ805d97BiFYKikgD7GseI8tL22MizcqJFu 1425 R+etjGR9x5WXE4XbgWdJsuJWGqQPvRz1c6/br/3ZCkCfrU7TMIXXstgVL+Wn+clPRoj7u/lpAyzq 1426 Of7D8FefqkNDnQ7ujCqwQUCRp2tCSSxT2kDEMl8/AjuQyIoceFPz6w3NLwAAAP//AwBQSwECLQAU 1427 AAYACAAAACEAtoM4kv4AAADhAQAAEwAAAAAAAAAAAAAAAAAAAAAAW0NvbnRlbnRfVHlwZXNdLnht 1428 bFBLAQItABQABgAIAAAAIQA4/SH/1gAAAJQBAAALAAAAAAAAAAAAAAAAAC8BAABfcmVscy8ucmVs 1429 c1BLAQItABQABgAIAAAAIQCOfwQ4JQQAANgQAAAOAAAAAAAAAAAAAAAAAC4CAABkcnMvZTJvRG9j 1430 LnhtbFBLAQItABQABgAIAAAAIQBXPXHe3wAAAAsBAAAPAAAAAAAAAAAAAAAAAH8GAABkcnMvZG93 1431 bnJldi54bWxQSwUGAAAAAAQABADzAAAAiwcAAAAA1414 DWfjBx8FqENaV9Cfn6yXghoaSJU3pACMk1R4U1c0kArQLV0Q2klXgJIAPiiZTKLEglCcXINyuC0I 1415 tXGKGZc9DOTIHA4wQli3tK2t2ZZLXX0aFajMrUEcdICYRu0sPTEV3ktD3Fo6ctVopO0uZ0ZAv7ru 1416 VMJD0Kl2NpwWK+2dNqKHqFt4LhGoPA+lEvwTHLzAWG/4gW652aK00/FkCp3TuGLyBx6cRch+V5Ef 1417 6efhhmkwnwYA3/l4tkZRbxmej4M4mbuz+WbRGbSFd2mgZkNDA88u9vRS/a+1MJDu4+gFSM0ltfY0 1418 TH8DV17VwA4M+Ok/jhdkyobjBqYSdAENDjHjGwjHnOKa6foI58HEFoLkdZXpVY2HFMVuWQt0wBD7 1419 ZqPPjSvLCzFgIJYZbSXF2dqNFa5qOwbrtal56DauMHXfMdT3ZR7M17P1LB7F42Q9ioPVavS8Wcaj 1420 ZBNOJ6totVyuwj805mGcllWWUaa962k4jP9ZJ3YXAkugJyK+iOIq2A183gbrX7phCANi6X8t2H03 1421 tryx49krdGZDaXCi4B4E3FZy8dlDHdwpFp78fY8F9VD9EwNeMQcGLiFmYsrXQ2K4shuuYEZA1cIj 1422 ClqAnSyVvbrsW1EVJdgKTcUy/gwUm1e6eQMjy9T65SZAbv8Wy0HLtlcIw3IzjbF2B3jw21jOIhwG 1423 M12mtpgd6SWa7gad6U56hr57Pr6TXnqTLO6kdye9O+n9n0nPvJ3Ci7OhcPeSr9/Mh3NDkue/Ip7+ 1424 BAAA//8DAFBLAwQUAAYACAAAACEAVz1x3t8AAAALAQAADwAAAGRycy9kb3ducmV2LnhtbEyPzU7D 1425 MBCE70i8g7VI3Kjz1zakcSqE6BlRuHBzYzexiNep7TaBp2c50duM9tPsTL2d7cAu2gfjUEC6SIBp 1426 bJ0y2An4eN89lMBClKjk4FAL+NYBts3tTS0r5SZ805d97BiFYKikgD7GseI8tL22MizcqJFuR+et 1427 jGR9x5WXE4XbgWdJsuJWGqQPvRz1c6/br/3ZCkCfrU7TMIXXstgVL+Wn+clPRoj7u/lpAyzqOf7D 1428 8FefqkNDnQ7ujCqwQUCRp2tCSSxT2kDEMl8/AjuQyIoceFPz6w3NLwAAAP//AwBQSwECLQAUAAYA 1429 CAAAACEAtoM4kv4AAADhAQAAEwAAAAAAAAAAAAAAAAAAAAAAW0NvbnRlbnRfVHlwZXNdLnhtbFBL 1430 AQItABQABgAIAAAAIQA4/SH/1gAAAJQBAAALAAAAAAAAAAAAAAAAAC8BAABfcmVscy8ucmVsc1BL 1431 AQItABQABgAIAAAAIQDQSJmzIgQAANIQAAAOAAAAAAAAAAAAAAAAAC4CAABkcnMvZTJvRG9jLnht 1432 bFBLAQItABQABgAIAAAAIQBXPXHe3wAAAAsBAAAPAAAAAAAAAAAAAAAAAHwGAABkcnMvZG93bnJl 1433 di54bWxQSwUGAAAAAAQABADzAAAAiAcAAAAA 1432 1434 "> 1433 1435 <v:shape id="Arc_x0020_37" o:spid="_x0000_s1027" style='position:absolute; … … 1456 1458 " path="m457200,nsc709705,,914400,204695,914400,457200r-457200,l457200,xem457200,nfc709705,,914400,204695,914400,457200e" 1457 1459 filled="f" strokecolor="red" strokeweight="1.5pt"> 1458 <v:path arrowok="t" o:connecttype="custom" o:connectlocs="0,0; 1,0"1460 <v:path arrowok="t" o:connecttype="custom" o:connectlocs="0,0;0,0" 1459 1461 o:connectangles="0,0"/> 1460 1462 </v:shape><v:shape id="Arc_x0020_38" o:spid="_x0000_s1028" style='position:absolute; … … 1484 1486 " path="m457200,nsc709705,,914400,204695,914400,457200r-457200,l457200,xem457200,nfc709705,,914400,204695,914400,457200e" 1485 1487 filled="f" strokecolor="red" strokeweight="1.5pt"> 1486 <v:path arrowok="t" o:connecttype="custom" o:connectlocs="0,0; 1,0"1488 <v:path arrowok="t" o:connecttype="custom" o:connectlocs="0,0;0,0" 1487 1489 o:connectangles="0,0"/> 1488 1490 </v:shape></v:group><![endif]--><![if !vml]><span style='mso-ignore:vglayout; … … 1500 1502 <p class=MsoNormal>The result looks very promising! Charge flipping operates 1501 1503 without regard to symmetry and GSAS-II attempts to locate the map with respect 1502 to the symmetry elements from the computed reflection phases. NB: charge1503 flipping also will give a solution without regard to handedness . This solution1504 may be either the correct one for sucrose or its reverse.Errors in these1504 to the symmetry elements from the computed reflection phases. (NB: charge 1505 flipping also will give a solution without regard to handedness; this solution 1506 may be either the correct one for sucrose or its inverse.) Errors in these 1505 1507 phases can disrupt this process so that the map and peak positions are offset 1506 1508 from their true positions. To discover this and fix it, rotate the drawing … … 1538 1540 src="Charge%20Flipping%20-%20sucrose_files/image015.jpg" v:shapes="Picture_x0020_12"><![endif]></span></p> 1539 1541 1540 <h2>Step 4 Obtain solution from charge flip result</h2>1542 <h2>Step 4: Obtain solution from charge flip result</h2> 1541 1543 1542 1544 <p class=MsoNormal>The next step is to extract the unique peak positions from 1543 the map and make atoms out of them. The map peak list is sorted with the1544 highest peak first; you can change the sorting to be along any of the column1545 the map and make atoms out of them. The map peak list is initially sorted by 1546 peak magnitude; you can change the sorting to be along any of the column 1545 1547 headings in the table by a single click on the desired heading. The <span 1546 1548 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: … … 1652 1654 src="Charge%20Flipping%20-%20sucrose_files/image022.png" v:shapes="Picture_x0020_17"><![endif]></span></p> 1653 1655 1654 <p class=MsoNormal>Notice that in this case the atom list is neatly sorted out into1655 two groups with a clear gap between M79 and M64. The higher ones are likely to 1656 be the 11 O-atoms and the rest are the 12 C-atoms. There may be a bit of 1657 ambiguity in the middle, but just assign the top 11 as O and the rest as C and 1658 see what happens. You can do this by selecting a block of atoms, go to <b1656 <p class=MsoNormal>Notice that in this case the atom list is neatly sorted out 1657 into two groups with a clear gap between M79 and M64. The higher ones are 1658 likely to be the 11 O-atoms and the rest are the 12 C-atoms. There may be a bit 1659 of ambiguity in the middle, but just assign the top 11 as O and the rest as C 1660 and see what happens. You can do this by selecting a block of atoms, go to <b 1659 1661 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1660 1662 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Edit/Modify … … 1678 1680 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1679 1681 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Draw atoms</span></b> 1680 tab. Pick the appropriate column and make the selection from the dialog1681 box.<span style='mso-spacerun:yes'> </span></p>1682 tab. Double click the appropriate column heading and make the selection from 1683 the dialog box.<span style='mso-spacerun:yes'> </span></p> 1682 1684 1683 1685 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 1688 1690 src="Charge%20Flipping%20-%20sucrose_files/image025.png" v:shapes="Picture_x0020_19"><![endif]></span></p> 1689 1691 1690 <p class=MsoNormal>Interestingly the atoms identified exactly matched with the1691 magnitudes of the charge flip map peaks with O at the top and C at the bottom.</p>1692 1693 1692 <p class=MsoNormal>Now one can safely remove the map peaks (do <b 1694 1693 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; … … 1748 1747 they will probably be outside the unit cell box. Do this on the drawing with <b 1749 1748 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1750 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Shift</span></b> 1751 left & right mouse buttons, they will show in green. </p> 1749 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Shift left</span></b> 1750 for the first one & <b style='mso-bidi-font-weight:normal'><span 1751 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1752 mso-hansi-theme-font:minor-latin'>Shift right</span></b> for <span 1753 class=SpellE>the</span> rest, they will show in green. </p> 1752 1754 1753 1755 <p class=MsoListParagraph style='margin-left:.75in;text-indent:-.25in; … … 1819 1821 src="Charge%20Flipping%20-%20sucrose_files/image027.jpg" v:shapes="Picture_x0020_20"><![endif]></span></p> 1820 1822 1821 <p class=MsoNormal>Save the project file; the next step is to <span 1822 class=GramE>do<span style='mso-spacerun:yes'> </span>the</span> Rietveld 1823 refinement on the completed structure.</p> 1824 1825 <h2>Step 5 Complete the structure</h2> 1823 <p class=MsoNormal>Save the project file; the next step is to complete the 1824 structure by finding H-atoms and then do the final Rietveld refinements.</p> 1825 1826 <h2>Step 5: Complete the structure</h2> 1826 1827 1827 1828 <p class=MsoNormal>To complete the structure you will need to perform a … … 1873 1874 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 1874 1875 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b> from the main 1875 GSAS-II data tree menu; all that was refined is scale and background. I got <span1876 GSAS-II data tree menu; only scale and background are refined. I got <span 1876 1877 class=GramE>a</span> <span class=SpellE>Rwp</span> ~20%. The plot shows clear 1877 1878 intensity differences; the atom positions obviously need refinement.</p> … … 1892 1893 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1893 1894 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Atoms</span></b> 1894 tab. Note the warning at the bottom of the window; the origin in space group P 1895 2<sub>1</sub> has an arbitrary location along the y-axis. Then double click the 1896 <bstyle='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif";1895 tab. Note the warning at the bottom of the window; the origin in space group P2<sub>1</sub> 1896 has an arbitrary location along the y-axis. Then double click the <b 1897 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 1897 1898 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>refine</span></b> 1898 1899 column heading, a dialog box will appear, select the <b style='mso-bidi-font-weight: … … 2004 2005 the <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 2005 2006 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Reflection 2006 set from</span></b> item (the only choice). The Fourier calculation will use the2007 structure factors shown in the <b style='mso-bidi-font-weight:normal'><span2007 set from</span></b> item (the only choice). The Fourier calculation will use 2008 the structure factors shown in the <b style='mso-bidi-font-weight:normal'><span 2008 2009 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2009 2010 mso-hansi-theme-font:minor-latin'>Reflection list</span></b> item for this PWDR … … 2049 2050 position. Recall that the refinement was just done without any H-atoms; it will 2050 2051 make the best fit compensating for the missing H-atoms by displacing the C 2051 & O atoms in the direction of the missing H-atoms. </p>2052 & O atoms in the direction of the missing H-atoms. </p> 2052 2053 2053 2054 <p class=MsoListParagraph style='margin-left:.75in;text-indent:-.25in; … … 2062 2063 mso-list:l4 level1 lfo14'><![if !supportLists]><span style='mso-fareast-font-family: 2063 2064 "Times New Roman"'><span style='mso-list:Ignore'>2.<span style='font:7.0pt "Times New Roman"'> 2064 </span></span></span><![endif]>Select a C or O atom swith possible H-atom2065 </span></span></span><![endif]>Select a C or O atom with possible H-atom 2065 2066 density nearby. It will appear green. Do <b style='mso-bidi-font-weight:normal'><span 2066 2067 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; … … 2102 2103 recognizable H-atom positions you can find; the goal is to find all 22 of them 2103 2104 (you probably will find less than that, but you could guess many of the rest). 2104 I found all but 5 of them.</p>2105 I easily found all but 5 of them.</p> 2105 2106 2106 2107 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2130 2131 mso-hansi-theme-font:minor-latin'>General</span></b> tab, adjust the <b 2131 2132 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 2132 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Contour level</span></b> 2133 and <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 2134 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Map radius</span></b> 2135 as desired. Look for the remaining H-atoms; you may have to guess positions for 2136 some of them using known stereochemistry rules. My model looked like</p> 2133 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Contour 2134 level</span></b> and <b style='mso-bidi-font-weight:normal'><span 2135 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2136 mso-hansi-theme-font:minor-latin'>Map radius</span></b> as desired. Look for 2137 the remaining H-atoms; you may have to guess positions for some of them using 2138 known stereochemistry rules. My model looked like</p> 2137 2139 2138 2140 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2203 2205 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; the <span 2204 2206 class=SpellE>Rwp</span> will start out high but should converge to where it was 2205 before this change. I got 8.08% and the Data tab looked like</p> 2207 before this change. I got 8.08% and the <b style='mso-bidi-font-weight:normal'><span 2208 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2209 mso-hansi-theme-font:minor-latin'>Data tab</span></b> looked like</p> 2206 2210 2207 2211 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2212 2216 src="Charge%20Flipping%20-%20sucrose_files/image043.png" v:shapes="Picture_x0020_42"><![endif]></span></p> 2213 2217 2214 <p class=MsoNormal> This no change in model yet but next select <b2215 style=' mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif";2216 mso- ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Generalized</span></b>2217 for the <span class=SpellE><b style='mso-bidi-font-weight:normal'><span 2218 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2219 m so-hansi-theme-font:minor-latin'>Mustrain</span></b></span><b2220 style=' mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif";2221 mso- ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'> model</span></b>;2222 this is the Stephens phenomenological model for <span class=SpellE>microstrain</span> 2223 which has 9terms for monoclinic structures. Check all the boxes for it.</p>2218 <p class=MsoNormal>Next select <b style='mso-bidi-font-weight:normal'><span 2219 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2220 mso-hansi-theme-font:minor-latin'>Generalized</span></b> for the <span 2221 class=SpellE><b style='mso-bidi-font-weight:normal'><span style='font-family: 2222 "Calibri","sans-serif";mso-ascii-theme-font:minor-latin;mso-hansi-theme-font: 2223 minor-latin'>Mustrain</span></b></span><b style='mso-bidi-font-weight:normal'><span 2224 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2225 mso-hansi-theme-font:minor-latin'> model</span></b>; this is the Stephens 2226 phenomenological model for <span class=SpellE>microstrain</span> which has 9 2227 terms for monoclinic structures. Check all the boxes for it.</p> 2224 2228 2225 2229 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2234 2238 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>; the <span 2235 2239 class=SpellE>Rwp</span> will start out high but should converge lower than 2236 where it was before this change. I got 7.65% and the Data tab looked like. If 2237 you check the <span class=SpellE>Mustrain</span> button in the Select plot type 2238 box, you will see the interesting <span class=SpellE>mustrain</span> surface 2239 plot for ground sucrose.</p> 2240 where it was before this change. I got 7.65% and the <b style='mso-bidi-font-weight: 2241 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 2242 minor-latin;mso-hansi-theme-font:minor-latin'>Data tab</span></b> looked like. 2243 If you check the <span class=SpellE><b style='mso-bidi-font-weight:normal'><span 2244 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2245 mso-hansi-theme-font:minor-latin'>Mustrain</span></b></span> button in the <b 2246 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 2247 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Select plot 2248 type</span></b> box, you will see the interesting <span class=SpellE>mustrain</span> 2249 surface plot for ground sucrose.</p> 2240 2250 2241 2251 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2341 2351 normal'><span style='font-family:"Calibri","sans-serif";mso-ascii-theme-font: 2342 2352 minor-latin;mso-hansi-theme-font:minor-latin'>109.54</span></b>. This is the 2343 tetrahedral angle so just press OK. 36 restraints will be listed; they are 2344 listed in order of the B atom. O-atoms should have one angle apiece; C-atoms 2345 will have more. </p> 2353 tetrahedral angle so just press <b style='mso-bidi-font-weight:normal'><span 2354 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2355 mso-hansi-theme-font:minor-latin'>OK</span></b>. 36 restraints will be listed; 2356 they are listed in order of the B atom. O-atoms should have one angle apiece; 2357 C-atoms will have more. </p> 2346 2358 2347 2359 <p class=MsoNormal><span style='mso-no-proof:yes'><!--[if gte vml 1]><v:shape … … 2367 2379 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2368 2380 mso-hansi-theme-font:minor-latin'>X-coordinates</span></b>; this will set those 2369 for all H-atoms. Now do <span style='font-family:"Calibri","sans-serif"; 2370 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Calculate/Refine</span>, 2371 the console will list progress on the Penalty function (i.e. the restraints) 2372 and the <span class=SpellE>Rwp</span> drops in my case to ~7.05%. One can increase 2373 the <b style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 2381 for all H-atoms. Now do <b style='mso-bidi-font-weight:normal'><span 2382 style='font-family:"Calibri","sans-serif";mso-ascii-theme-font:minor-latin; 2383 mso-hansi-theme-font:minor-latin'>Calculate/Refine</span></b>, the console will 2384 list progress on the Penalty function (i.e. the restraints) and the <span 2385 class=SpellE>Rwp</span> drops in my case to ~7.05%. One can increase the <b 2386 style='mso-bidi-font-weight:normal'><span style='font-family:"Calibri","sans-serif"; 2374 2387 mso-ascii-theme-font:minor-latin;mso-hansi-theme-font:minor-latin'>Restraint 2375 2388 weight factors</span></b> to tighten them up at some expense in <span
Note: See TracChangeset
for help on using the changeset viewer.