Stimulating Simulations

Dr. Stuart Wright, Senior Scientist, EDAX

It has been an interesting experience to build our OIM Matrix™ software package. As you may know, OIM Matrix is partially a front-end user interface to the EMsoft package developed by Professor Marc De Graef’s group at Carnegie Mellon University to make it convenient to use within the framework of OIM Analysis™. I have learned a lot in the process and am grateful for Marc’s patience with my many questions. Will Lenthe recently joined the EBSD group at EDAX. Will worked as a Post-Doc in Marc’s group, and his additional insights have been invaluable as we are striving to build the second generation of OIM Matrix. It will be easier to use, more robust, and provide some significant speed gains.

While our initial focus for OIM Matrix was on helping users improve the indexing of EBSD patterns from difficult-to-index materials, I’ve been surprised by how useful it has been for testing our software. It has also helped us in developing some of our new features. Having well-simulated patterns for known orientations and EBSD/SEM geometries is very helpful.

I used OIM Matrix for a study on feldspars. According to Wikipedia:

“Feldspars are a group of rock-forming aluminum tectosilicate minerals containing sodium, calcium, potassium, or barium. The most common members of the feldspar group are the plagioclase (sodium-calcium) feldspars and the alkali (potassium-sodium) feldspars. Feldspars make up about 60% of the Earth’s crust and 41% of the Earth’s continental crust by weight.”

Given that feldspars are relatively common, we are frequently asked to help index them. They are difficult, as a poster at the 2019 Quantitative Microanalysis (QMA) conference detailed [1]. I thought it might be interesting to see what we could learn about the limits of EBSD in characterizing these materials. I won’t give you all that we learned in that little study, but what I thought was an interesting snapshot. Figure 1 shows a phase diagram for the feldspar group of minerals.

Figure 1. Phase diagram for the feldspar group.

To start, I looked in the American Mineralogist Crystal Structure Database (AMCSD) for all the relevant entries I could find. There are a lot of variants. Here is a table:

Table 1. Number of entries in AMCSD for each feldspar.

I enjoy seeing pattern simulation results, but producing 149 master patterns [2] would take more patience than I have (each master pattern calculation can take several hours for these low-symmetry materials). So, I selected one entry for each mineral type. I tried to find one that seemed most representative of all the other entries in the set. After calculating the eight master patterns, I simulated one individual pattern at the same orientation for each mineral, as shown in Figure 2. Note that they are all similar, with the most deviation coming from the anorthite and sanidine end members of the series.

Figure 2. Patterns were simulated at Euler angles of (30°, 30°, 30°) for each feldspar.

To quantify the differences, I calculated the normalized dot-products [3] for all pattern pairs to get the following table. A value of “1” indicates the patterns are identical. As expected by the initial observation, the biggest difference is the sanidine to albite pair of patterns.

Table 2. Normalized dot products.

Of course, the next step would be to see how this holds up to experimental patterns and dictionary indexing [4]. I hope to eventually do this with samples Professor Rudy Wenk of Stanford University kindly gave me. Rudy has been one of the major contributors to the entries in the AMCSD for feldspars.

There was one more virtual experiment I thought would be interesting. I wanted to ascertain how much the chemical species in the feldspar series influenced the patterns. To do this, I created an average structure instead of using the lattice parameters for each feldspar. I then populated these structures with atoms to maintain the chemical composition ratios specified for each series. A master pattern for each ideal structure was calculated. Three hundred forty patterns were simulated uniformly, covering orientation space with a spacing of approximately 30° between orientations. The average normalized dot products were calculated for each pattern against the albite pattern at the same orientation. Figure 3 shows the results.

Figure 3. The normalized dot product of simulated patterns for idealized structures against the albite simulated patterns.

Clearly, the dot products are all very near 1, indicating that the differences in the simulated patterns due to chemical composition are small for these chemical species. This suggests that coupling EBSD with EDS is critical when trying to differentiate the different feldspar minerals. While this small study has not changed the world of feldspar indexing, it has, at least, been a stimulating study of simulating for me.

[1] B Schneider, and J Fournelle (2019) “Using Quantitative and Qualitative Analysis to Confirm Phase Identities for Large Area EBSD Mapping of Geological Thin Sections” Poster at Microanalysis Society Topical Conference: Quantitative Microanalysis, University of Minnesota, Minneapolis MN, June 2019.

[2] PG Callahan, and M De Graef (2013) “Dynamical electron backscatter diffraction patterns. Part I: Pattern simulations” Microscopy and Microanalysis, 19, 1255-1265.

[3] S Singh, and M De Graef (2016) “Orientation sampling for dictionary-based diffraction pattern indexing methods” Modelling and Simulation in Materials Science and Engineering, 24, 085013.

[4] K Marquardt, M De Graef, S Singh, H Marquardt, A Rosenthal, and S Koizuimi (2019) “Quantitative electron backscatter diffraction (EBSD) data analyses using the dictionary indexing (DI) approach: Overcoming indexing difficulties on geological materials” American Mineralogist: Journal of Earth and Planetary Materials, 102, 1843-1855.

Leave a ReplyCancel reply