Feature request: Slip transmission in ORTools

[hgtdwdvjiu8876543-debug]

Hi @AzdiarGazder
I hope you are doing well. I have recently started using MTEX for crystallography-related research and have become familiar with the basic workflow. I believe crystallographic approaches can be highly valuable for analysing the strength–toughness balance of materials; however, I have encountered some practical difficulties in reproducing a published figure. I would greatly appreciate your guidance.

Specifically, I am trying to reproduce Fig.10(a-c) from the following paper: https://doi.org/10.1016/j.jmrt.2025.11.047

[AzdiarGazder]

Hi @hgtdwdvjiu8876543-debug,

Fig. 10 in the paper you reference uses the (1-m'_max) parameter calculated using the STABiX - Slip Transmission Matlab toolbox for all variant pair boundaries formed with V1, (i..e.- V1 - V1-24), for the ideal K–S OR.

What is missing from the paper is:

1. Exactly which slip systems were used to do this calculation on either side of the variant boundary pair, and
2. The justification as to why it was only used for the ideal and rational K-S OR case.

I suspect the answers are:

1. Every slip system combination was tried for a given variant boundary pair, and the one with the highest m'_max value was selected.
   If this is the method used, it lacks a fundamental basis - as deformation accommodation via activation of a slip system within a grain heavily depends on the type and direction of the externally imposed deformation and the orientation of the grain itself.
2. The study is only based on an ideal rational K-S OR case - not the irrational ORs that operate in reality.

Statistics from slip transmission is definitely a useful feature to add to ORTools.
But I'm afraid there is no simple working solution to your question just yet.
It will take some time for me to think of a workable solution.

Warm regards,
Azdi

[AzdiarGazder]

Hi @hgtdwdvjiu8876543-debug ,

The issue was closed because it was resolved.

If you read through the entire message thread you'll see the following message - pasted below - from me.

Hope this helps! 🙂
Have a good weekend.

Warm regards,
Azdi

Hi @wanghaijian217,

A new function computeGrainPairs.m has been added to the ORTools Developer to help you calculate grain pairs based on variant, packet, Bain group id OR any other type of user-specified id. The advantage of the last option is that you can now return grain pair information based on equivalent ids or any other type of id that you define.

Please see Example 10 for a detailed demonstration on how to use this new function.

Hope this helps with your analysis.

Warm regards,
Azdi

[hgtdwdvjiu8876543-debug]

Hi @AzdiarGazder

Following your tutorial, I was able to complete the full workflow (FCC→BCC reconstruction and KS variant assignment) and successfully run it on my own EBSD dataset.

However, when I visualize KS variant-pair boundaries and the corresponding density/frequency histogram using plotMap_KSvariantPairs (with the parentGrainId option), the output only shows a limited set of V1-related variant pairs:
V1–V2, V1–V3 (V5), V1–V4, and V1–V6.

In contrast, the following V1-related pairs do not appear in either the boundary map or the histogram:
V1–V7, V1–V8, V1–V9 (V19), V1–V10 (V14), V1–V11 (V13), V1–V12 (V20), V1–V15 (V23), V1–V16, V1–V17, V1–V18 (V22), V1–V21, and V1–V24.

For reference, below are example outputs of what I am expecting to obtain:

Warm regards,
shuai

[AzdiarGazder]

Hi @hgtdwdvjiu8876543-debug ,

Again - please click on this link to see Example 10.
Please refer the following section between lines 162 and 233 starting with:

%% OPTION 2:  Groups of variant id for child grain pair analysis
...
%%

It contains a detailed demonstration on how to get this graph.
Hope this helps.

Warm regards,
Azdi

[hgtdwdvjiu8876543-debug]

Hi @AzdiarGazder

Thanks a lot for the quick reply and for pointing me to Example 10.

Sorry if my previous message was unclear. My understanding is:

Option 2 gives the KDE/density plots for the 16 child grain variant pairs, but it does not plot the corresponding boundaries (the “plot the boundaries” result I’m trying to reproduce).

The boundary-style plots I want are the ones produced in Option 3 in Example 10. However, in the current script it seems Option 3 only demonstrates a subset of pairs (e.g., V1–V2, V1–V3(V5), V1–V6, V1–V4 / 3 figures).

Thanks again for your help!

Warm regards,
shuai

[hgtdwdvjiu8876543-debug]

Hi @AzdiarGazder @frankNiessen

Sorry to bother you again. I’m writing to follow up on the discussion regarding “Slip transmission in ORTools” (Issue #99). About two weeks ago, you mentioned that there might not be a simple solution and that you would need some time to think of a feasible approach.     
May I ask if there has been any progress or any preliminary direction you would recommend? Even a rough idea of the next steps (e.g., what assumptions or additional inputs might be needed) would be extremely helpful for my current work.     

If it would help, I’m happy to provide more details about my specific use case, share example datasets, or test any candidate implementation on my side.     

Thank you very much for your time and guidance.     

Warm regards,
shuai