Next Article in Journal
Active Bionanocomposite Coating Quality Assessments of Some Cucumber Properties with Some Diverse Applications during Storage Condition by Chitosan, Nano Titanium Oxide Crystals, and Sodium Tripolyphosphate
Next Article in Special Issue
Interplay between Habit Plane and Orientation Relationship in an Electron Backscatter Diffraction Analysis: Using the Example of η′-Al8Fe3 in η-Al5Fe2
Previous Article in Journal
Precipitant-Free Crystallization of Lysozyme and Glucose Isomerase by Drying
Previous Article in Special Issue
Crystal Structure, Microstructure and Electronic Properties of a Newly Discovered Ternary Phase in the Al-Cr-Sc System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Correspondence Theory and Its Application to NiTi Shape Memory Alloys

EPFL (Ecole Polytechnique Fédérale de Lausanne), LMTM, PX Group Chair, 2000 Neuchatel, Switzerland
Crystals 2022, 12(2), 130; https://doi.org/10.3390/cryst12020130
Submission received: 23 November 2021 / Revised: 11 January 2022 / Accepted: 13 January 2022 / Published: 18 January 2022
(This article belongs to the Special Issue Crystallography of Structural Phase Transformations (Volume II))

Abstract

:
Martensite crystallography is usually described by the phenomenological theory of martensite crystallography (PTMC). This theory relies on stretch matrices and compatibility equations, but it does not give a global view on the structures of variants, and it masks the relative roles of the symmetries and metrics. Here, we propose an alternative theory called correspondence theory (CT) based on correspondences and symmetries. The compatibility twins between the martensite variants are inherited by correspondence from the symmetry elements of austenite. We show that, for the B2 to B19′ transformation, there is a one-to-one relation between the specific misorientations and the specific inter-correspondences between the variants. For each type of misorientation, the twin of its junction plane can be predicted without calculating the stretch matrices, as in PTMC. The rational elements of the twins do not depend on the metrics; all the transformation twins are thus “generic”. We also introduce the concept of a weak plane that permits to explain the junction planes for polar pairs of variants for which the PTMC compatibility equations cannot be solved. The predictions are validated by comparison with experimental Transmission Kikuchi Diffraction (TKD) maps.

1. Introduction

Shape memory alloys (SMA) such as nickel-titanium (NiTi) are widely used in stents, springs, actuators, sensors and connectors. Their astonishing mechanical properties directly result from the crystallography of the austenite (A) to martensite (M) phase transformation [1,2,3,4]. Depending on the composition, the alloy is martensitic or austenitic at room temperature, which dramatically affects its mechanical behavior. The shape memory effect occurs in martensitic alloys and is explained as follows. When a hot austenitic straight bar is cooled down to room temperature, the austenite grains are transformed into martensite variants, and the bar remains straight. The bar is soft at room temperature and can be bent at low force; it keeps its deformation when the stress is released, with nearly no elastic return. This unusual plastic behavior comes from the fact that the deformation is not mediated by dislocations but by variant reorientation under stress. When the bent bar is heated up again, the martensite variants are transformed back into their parent austenite grains, and the bar returns to its initial straight shape. The superelasticity effect is different and occurs in metastable austenitic alloys. When the bar is bent at room temperature, the austenite is first deformed elastically and becomes unstable; the austenite grains are progressively transformed into martensite variants that are well-oriented in the stress field. When the stress is released, the bar comes back to its initial austenite state and to its straight shape in a way that looks like elasticity but with a larger amplitude.
Before the 1960s, optical microscopy was the unique technique to visualize the martensitic domains, but no direct link could be made between martensite morphology and crystallography. Everything changed with Transmission Electron Microscopy, which has permitted to study simultaneously the shapes, the internal twins, the disorientations and the junction planes of martensite at nm–μm scales [5,6,7,8,9]. High-Resolution TEM (HRTEM) reveals fine details of the microstructure at the atomic level [10,11,12,13,14,15]. Despite its undisputable interest, one should recognize that TEM is however limited to small fields of view and does not permit to get statistical information on a large number of martensitic domains because of the difficulty of aligning them along zone axes. Digital Image Correlation (DIC) strain mapping with optical microscopy images can cover large dimensions (cm–mm); when associated with Scanning Electron Microscopy (SEM) images, it can reach spatial resolutions of a few μm [16,17], and when coupled with in situ X-ray diffraction, it can give precious information on the shear systems implied in the phase transformation [18,19,20]. However, these techniques do not allow studying locally the morphologies, the orientations and the twins of martensite.
Phase and orientation maps can nowadays be acquired by Electron Back-Scatter Diffraction (EBSD) on large surfaces (up to some mm2) with nanometric resolution thanks to field emission gun scanning (FEG) SEM and to the recent development of Complementary Metal Oxide Semiconductor (CMOS) cameras. The speed can reach 3000 pixels/s, and very large maps > 1 million pixels are easily acquired with a very good angular resolution (<0.1°) and spatial resolution (<5 nm). In addition, Transmission Kikuchi Diffraction (TKD) can be attempted as a complementary technique to improve the spatial resolution down to 2 nm. TKD works exactly as EBSD, except that the sample is no longer a bulk specimen tilted at 70° but a TEM lamella positioned close to the SEM pole piece and tilted at −20°, such that the Kikuchi patterns are collected from the bottom surface of the specimen [21,22,23]. The significant progresses made in SEM with the possibility to use high currents without deteriorating the spatial resolution (1 nm at 20 nA and 30 kV) render this technique very attractive. The first EBSD and TKD maps of B19′ martensite variants in NiTi alloys were recently presented in Ref. [24]. Many important results could be obtained. It was shown that the microstructure is made of large B19′ laths visible in EBSD and of sub-μm twins visible by TKD. There are a predominant orientation relationship (OR) and many others ORs close to it with a continuum between them. The dominant OR between austenite and B19′ martensite was noted: AQ; it is also the “natural” OR, i.e., that for which the dense planes and dense directions of austenite and martensite are parallel. The habit plane of the large martensitic laths predicted from the individual lattice distortion associated with the OR AQ is (1 1 ¯ 2)B2  (10 1 ¯ )B19′, which is in good agreement with literature (see Section 2) and with the EBSD maps [24]. The additional ORs were interpreted as “closing-gap” ORs required to maintain the compatibility between variants. A general picture explaining the variants, junction planes, additional ORs and continuums of ORs was proposed. We also made the hypothesis that, for two distortion variants in contact, there is at least one parent symmetry element of the parent phase that is preserved by correspondence and that the variants remain linked together by this symmetry element. The rotation gradients between the natural OR and the closing-gap ORs permit to maintain locally this parent symmetry element. For a parent mirror symmetry, the variants remain linked by a type I twin through the mirror plane, and for 180° rotation axes, they remain linked by a type II twin around the rotation direction. For variants linked by a parent symmetry that is a non-two-fold operation, the twin cannot be of type I or II, but we proposed that it could be a “weak twin”, a concept that will be explained in Section 4.4. It was concluded without giving details that it should be possible to predict the junction planes by using the inter-correspondence operators.
The aim of the study is to explain in detail the principles of the theory that will be called “correspondence theory”, give the details of the calculations and compare the results to the phenomenological theory of martensite crystallography (PTMC) and to experimental TKD maps. The paper is quite long, because it is self-consistent, but some parts related to the PTMC can be easily skipped by readers already familiar with this theory.

2. Toward a Change of Paradigm to Understand the Crystallography of Martensite?

2.1. The Phenomenological Theory of Martensite Crystallography

The crystallographic features of martensite (habit planes, junction planes and orientations) are usually discussed with the PTMC. This theory was born from optical observations of the twinned structures of AuCd alloys in the 1940s, but its development in the 1950s was mainly driven by the will to explain the irrational habit planes of martensite in steels. It was clear that martensitic transformations are non-diffusive and share many characteristics with deformation twinning. The straight boundaries or the midribs of martensite in steels made researchers think that martensite necessarily results from an invariant plane strain (IPS) deformation. An IPS is a simple shear combined with an expansion perpendicular to the shear plane. It is also called “shape strain”, because it is the strain associated with the morphology of the martensite product. However, IPS alone was not sufficient to build the theory. In steels, the face-centered cubic (fcc) lattice of austenite cannot be directly transformed by an IPS into the body-centered cubic (bcc) or body-centered tetragonal (bct) lattice of martensite. The PTMC solves this issue by assuming that the martensite products (laths or lenticles) are actually “composites”. Two versions of this idea were proposed nearly simultaneously by two groups of researchers: Bowles and Mackenzie [25,26] and Wechsler, Liebermann and Read [27]. In the first version, the martensite structure is supposed to be constituted of one variant and regular defects that do not change the lattice. These defects are dislocations or mechanical twins, and their strain is called “lattice invariant strain” (LIS) or “lattice invariant deformation” (LID). In the second version of the theory, the martensite product is supposed to be constituted by two variants linked by a twin relation. Both versions assume that the martensite product is inhomogeneous and constituted of parts linked together by a simple shear; they are mathematically equivalent [28]. Another cornerstone of the PTMC is the existence of a stretch distortion, which, in steels, is the well-known Bain distortion that links the fcc and bcc structures [29]. All these crystallographic elements (IPS, LID and Bain) are conciliated thanks to a free rotation, such that the Bain stretch combined with this rotation gives the lattice distortion of one variant, which, when combined with the LID, gives the IPS. There are numerous books about PTMC, and among the most didactic ones on the first version of PTMC are those written by Bhadeshia [30] and Wayman [31]. Bowles and Mackenzie’s version of PTMC with twinning for LIS has been used for a long time for SMA. It is thus necessary to come back on the theory of twinning.
Twinning has a long history [32]. In the 1880s, Mügge [33] found the general equations that describe the orientation relationship between a crystal and its twin. He introduced the four geometrical elements of twinning: the plane K1 and the direction η1 parallel to K1 and their conjugates K2 and η2 parallel to K2. It is important to note that K1 and η2 are rational (i.e., the indices of the plane and direction are small integers), and K2 and η1 are, in general, irrational. There are two types of twins. Type I twins are characterized by a shear plane K1 and a shear direction along η1 and type II twins by a shear plane K2 and a shear direction along η2. Both modes are conjugate and have the same shear amplitude s. When all the twin elements K1, η1, K2 and η2 are rational, the twin is called “compound”. For type I twins, K1 is as a mirror plane for the twin edifice (it is also the interface between the two individual crystals). For type II twins, the direction η2 is a 180° rotation symmetry of the edifice. The mathematical approach initiated by Mügge was continued in 1950–1970 by Kihô [34], Jaswon and Dove [35] and Bilby, Bevis and Crocker [36,37]. Among the excellent reviews about twinning, one can cite Cahn [38] and Christian and Mahajan [39].
PTMC has used for a long time in the twin systems calculated by Bilby, Bevis and Crocker’s method, but unfortunately, the list of predicted twins can be quite long, and there is no absolute criterion to choose the twins that should solve the PTMC equations. Therefore, instead of guessing the twin system that should act as LID, researchers have used TEM experimental observations to introduce them into the equations. Otsuka et al. observed by TEM that (11 1 ¯ )B19′ was a twinning plane between some variants [5]. Other TEM studies showed other junction planes, such as the (001)B19′, (100)B19′ and (011)B19′ planes [40,41]. These twins were identified as type I twins. Knowles and Smith [6] were the first to report a <011>B19′ type II twin in NiTi alloys. For these twins, the expected junction plane is close to (0.72,1, 1 ¯ )B2 ≈ (3,4, 4 ¯ )B2. This plane was later confirmed by other TEM and HRTEM studies [11,42]; the irrational interface is made of edges and ledges on lower indices planes that remain the subject of theoretical studies [43]. The outputs of the PTMC are the OR between austenite and martensite and the habit plane between the two phases. The habit plane is the plane parallel to the shape of the “composite” lath; it is the interface plane with austenite. It should not be confused with the junction plane that is the interface between the variants that constitute the “composite” lath. Miyazaki et al. [44] could experimentally identify two groups of habit planes, one around (5,6,14)B2 and one around (8,9,14)B2, with a large scattering of ±4°. They compared them to the PTMC predictions by considering the two LIDs already assessed by Knowles and Smith [6], i.e., the (11 1 ¯ )B19′ type I twins and the <011>B19′ type II twins, and they found that the one closest to the experience was (0.215, 0.405, 0.888)B2. A few years later, Miyazaki, Otsuka and Wayman [45] measured the habit planes close to (0.39,0.48,0.78)B2, which is the intermediate between the two previously reported habit planes.
The fact that the junction (twin) planes had to be picked up in a theoretical list or chosen among those observed by TEM has left researchers unsatisfied. Fortunately, Wechsler, Liebermann and Read’s version of PTMC evolved at the turn of the 1990s with the works of Ball, James, Bhattacharya, Pitteri and Zanzotto [2,46,47,48,49,50,51]. They established a mathematical formulation in which the twin is not anymore an arbitrary choice but an output of the calculations required to maintain the compatibility between the variants. The main idea is that the distortion matrices of two variants expressed by polar decomposition by F i and F j are “compatible” at the junction plane i/j if and only if they are rank-1 connected, which means that there is a plane of normal n and a direction a in this plane such that F i F j = a n . The compatibility criterion is thus equivalent to a simple shear, as for deformation twinning, but this shear is no longer an input. The method that solves rank-1 equations uses polar decompositions of F i = Q i   U i and F j = Q j   U j , where the matrices U are “Bain” stretches and the matrices Q are “free” rotations. The equation can be written Q i j   U i U j = a n , where Q i j = Q j 1 Q i and a = Q j 1 a . It is solved by calculating the eigenvalues λ 1 , λ 2 and λ 3 and eigenvectors e 1 , e 2 and e 3 of the matrix F j t   F i t   F i   F j 1 = U j t   U i t   U i   U j 1 . Some solutions exist for a and n if and only if λ 1 1 , λ 2 = 1 , λ 3 1 , and they are expressed as linear functions of e 1 and e 3 , with coefficients that are square roots of fractional functions of λ 1 , λ 2 and λ 3 . The habit plane of a “composite” martensite product made of the two variants i and j is calculated by assuming that it is the plane of the IPS; its equation is λ F i + ( 1 λ ) F j = I + b n , where λ a real number between 0 and 1 that represents the volume faction of each variant. This austenite–martensite compatibility equation is solved thanks to another intermediate matrix (the details are skipped here), and the solutions λ , b and n are expressed as functions of the eigenvalues and eigenvectors of that matrix.
This “modern” and mathematized form of PTMC was applied to NiTi alloys by Hane and Shield in 1999 [49]. The calculations were made for all the pairs of variants; the results could be grouped according to six sets of pairs. For four sets, the possible twins were found on prior (100)B2 or (110)B2 for type I twins or along the <100>B2 or <110>B2 axes for type II twins, which, when expressed in martensite coordinates, give the junction planes ( 011 ) B 19 and ( 11 1 ¯ ) B 19 and the irrational plane close to ( 34 4 ¯ ) B 19 in agreement with the “older” version of PTMC. No solution could be found, however, for the two sets of pairs of variants linked by 90° and 120° rotations. It will be shown that, actually, even for these sets of pairs, junction planes exist, and they can be predicted with the new concept of an “axial weak twin”. By a classification of the different sets of pairs quite similar to Hane and Shield’s one, Pitteri and Zanzotto noticed that two-fold symmetry twins (180° rotations or reflections) did not depend on the exact values of the stretch matrices U i ; these twins were called “generic” [47]. These authors also came to conclude that the pairs of variants linked by the 90° and 120° rotations are “non generic”. It will be shown in Section 6 that, actually, all the twins are generic for their rational elements, i.e., the twin plane K1 for type I twins, the twin direction η 2 for type II twins and the axial direction for the axial weak twins. For the last two decades, most of the experimental observations on SMAs have been interpreted with the modern PTMC, such as the “hearing-bone” assemblies made by two groups of two B19′ variants linked by ( 001 ) B 19 junction planes [10] or the “hexangular” assemblies made of three groups of two variants [7,8,50]. One of the recent advances of the PTMC lead to the design of new alloys with specific compositions for which the metrics of the parent and daughter phases follow extra “super-compatibility” conditions [51,52].
The modern PTMC is quite simple in its principles, since it relies on compatibility conditions between the variants, but the calculations are quite tricky because of the numerous cases to be treated and the equations to be solved. Some problems also remain. For example, it was shown that the habit planes of the composite martensite products (i,j) and (k,l) are not compatible. This issue was already studied and mathematically solved by Bhattacharya by introducing extra rotation matrices J i j k l to obtain a compatibility at the (i,j)/(k,l) interfaces [2]. However, the same problem could be repeated at larger scales by considering two larger assemblies (i,j) + (k,l) and (m,n) + (o,p), etc. The modern PTMC relies on the free rotations Q i j , J i j k l , etc., but these additional rotations depend on the pairs or quadruplets of variants that are considered, which renders irrelevant the notion of austenite–martensite OR that was however judged to be of prime importance in the earlier versions of PTMC. The PTMC also remains mute on atomic displacements, which has forced some researchers to develop alternative approaches based on molecular dynamics to get a better understanding of the rearrangement of the atoms during the austenite–martensite transformation [53,54].

2.2. The Hidden Algebraic Structure of Variants

The PTMC was born from the observations of twins in martensite in AuCd alloys. Their presence is essential to Wechsler, Liebermann and Read’s version of PTMC and its modern mathematized form. The author’s research on martensite started with a study of low-carbon steels. Martensite is made of intricate bcc laths that do not contain twins, which makes modern PTMC nearly irrelevant. An alternative approach was required. By using group theory, it can be shown that the orientation variants are simple cosets built from the group of symmetries common to the parent and daughter phases called the “intersection group”. The specific misorientations between the variants are isomorph to double cosets also built on the intersection group. The double cosets were called “operators”. The use of cosets and double cosets was initially proposed by Janovec for ferroelectrics [55,56]. The number of variants and operators is given by Lagrange and Burnside’s formulas, respectively. The algebraic structure of the variants and their operators is not a group, as often claimed in the literature, but a groupoid [57]. A computer program called GenOVa for the “generation of orientational variants” was written in Python to algebraically calculate the variants and the groupoid composition table; it also simulates the electron diffraction patterns and the pole figures of a parent crystal with its variants [58]. The results are used in another computer program called ARPGE that reconstructs the prior parent grains from experimental EBSD maps of martensite [59]. Thanks to ARPGE, odd continuous features in the pole figures of the bcc martensite variants in the prior austenite grains in low-carbon steels were observed. They could not be clarified by the usual plasticity or by the PTMC [60]. Two-step and one-step atomistic models were then proposed to explain these experimental results [61,62]. A hard-sphere approximation was introduced to obtain the analytical equations of atomic displacements during the fcc → bcc lattice distortion [62]. The same approach was applied to the transformations in fcc-bcc-hcp systems [63], assuming that, for each transformation between these phases, there is a “natural” OR for which the dense directions and planes are parallel. This natural OR is the Kurdjumov-Sachs OR for fcc-bcc, Burgers OR for bcc-hcp and Shoji-Nishiyama OR for fcc-hcp transformations. In these works, the concept of simple shear was replaced by that of “angular distortion”. A simple explanation could be proposed for the {225} habit planes of twinned martensite in high-carbon steels [64]. The results were the same as those expected by PTMC if one assumes that the values of the lattice parameters are those of an ideal hard-sphere packing and not the exact ones. Simulation movies showing at the atomistic displacements during the fcc-bcc transformation could be proposed for the first time. The variant selection of martensite in steels was also quantified with the interaction work between the lattice distortion and the external stress [65,66]. Mathematically, the interaction work is the product W = Γ i j   i j , where Γ is the matrix of the external stress field, and is the deformation matrix calculated from the distortion matrix F by = F I . The simple use of the distortion matrix F without the complex machinery of PTMC was successful to understand the variant selection in AuCu red gold alloys [67,68]. The knowledge of distortion matrix F associated with the natural OR is also sufficient to predict the habit planes of martensite in NiTi alloys [24].
The approach developed over the last decade has allowed us to explain many features of martensite, but it does not treat, in a general way, the twins and junction planes observed in shape memory alloys. The aim of the paper is thus to propose a method to calculate the twins and the junction planes. We will show that there is no need to consider the stretch distortion matrices, such as in PTMC; the knowledge of the correspondence and symmetries is sufficient to predict the nature of the twins. The theory will be called “correspondence theory” (CT). Contrarily to the distortion or stretch matrices, the correspondence matrix does not explicitly contain information on the metrics. If one considers a crystallographic direction (or plane) of the parent phase, the correspondence tells in which direction (or plane) of the martensite phase this direction is transformed. The correspondence matrices are constituted of integers or simple rational numbers (Section 5). The theory is built on the algebraic structure of variants and their operators (Section 3.3). The sets calculated by Hane and Shield [49] are actually double cosets, and the rational element of the twins is directly given by the parent symmetries in the double cosets. The junction planes, the closing-gap ORs and the orientation gradients will be predicted only from the inter-correspondence matrices. The twins will be shown to be identical those calculated by the PTMC, but their meaning is clearer, and our approach permits to determine which of the type I or type II is favored, whereas the PTMC remains mute on this point, because both have the same shear amplitude. The concept of an axial weak twin will be introduced for the pairs of variants linked by the 90° and 120° rotations. In the second part of the paper, the predictions will be compared to the experimental TKD maps.

3. Crystallography for Phase Transformations

3.1. Directions and Planes

Some elements of crystallography are briefly recalled, and the notations are introduced. The vectors are written in bold lower cases and the matrices in bold capital letters. The exceptions will be for the symmetry matrices, noted g when they are elements of a point group and, more specifically, m for a reflection symmetry and R for a rotational symmetry. A set of matrices is noted in bold capital letters and a point group in double-struck letters, for example, G . A vector d of the direct space is in a column, and a vector p* of the reciprocal space is in a line. The same reciprocal vector is simply written p when it is in a column, i.e., p * = p t with the symbol “t” in the superscript meaning “transpose of”. A scalar product is calculated by multiplying term-by-term the coordinates of a vector of the reciprocal space and a vector of the direct space and summing them, for example, p * d = p t d = p i d i , with Einstein’s convention. The dyadic product of the vectors d p = ( d . p t ) is the matrix d i p j . The dyadic product notation has all the properties of the matrix product, for example, ( d   p t )   u = ( p t   u ) , which is the vector d multiplied by the scalar product p t   u .
We recall that, if a distortion matrix F acts on the vectors u of the direct space, d = F   d , the same distortion acts on the plane p by p = F *   p , with F * = F t (inverse of the transpose of F). It can be checked that any direction d lying on the plane p remains on the plane after distortion, since p t   d = ( F t   p ) t   F   d = p t d = 0 . A plane p unrotated by F is an eigenvector of F * . A plane is “globally invariant” by a distortion F when p = F *   p = p . If, in addition, all the vectors in the planes are invariant by F, i.e., u   such   that   p t u = 0 ,   F   u = u , the plane is said to be “fully invariant”. For any crystal, a crystallographic basis c = ( a ,   b ,   c ) formed by the usual crystallographic vectors can be defined. At the basis, c can be associated to a 3 × 3 unit cell matrix B c = [ a ,   b ,   c ] by writing the coordinates of the vectors a ,   b   and   c in columns in a unit orthonormal reference frame. The metric of the crystal is defined by the metric tensor
= B c t   B c = [ a 2 b t   a c t   a a t   b b 2 c t   b a t   c b t   c c 2 ]
The metric tensor is the coordinate transformation matrix between the reciprocal space and the direct space = [ c * c ] . It has the properties to be symmetric = t and * = [ c c * ] = 1 . The scalar product between the vectors u and v of the direct space is determined by expressing one of the two vectors in the reciprocal space, thanks to the metric tensor, ( u · v ) = u t v . The norm d of a vector d of the direct space and the norm p * of a vector p of the reciprocal space are, respectively, given by d = d t d and p * = p t * p . The notation d ˜ applied to a direct vector means that d is normalized by d , and the notation p ˜ applied to a reciprocal vector means that p is normalized by p * , i.e., d ˜ = d d * and p ˜ = p p * . The inter-reticular distance d h k l between the planes p of Miller indices p = (h,k,l) is d h k l = 1 p * . The unit normal direction in the direct space of a plane p is n given by n = * p ˜ . It can be verified that n t n = n t   p ˜ = 1 .
It is possible to introduce orthonormal bases linked to the crystallographic basis c by following the instructions given in the crystallography textbooks, but their choice is actually arbitrary. If we note # = ( x , y , z ) as one of these bases, the link with c is given by a matrix called a “structure tensor” S = [ # c ] . It can be easily proven that S t   S = , whatever the instruction that is chosen to calculate S .
We note u A , a crystallographic direction (rational indices) of the austenite lattice (A). This direction can be expressed in other bases, for example, in the crystallographic basis c M of martensite (M); in that case, it is specified u / c M A or just u / M A . We will sometimes use a short notation u A for a direction of an austenite crystal expressed in the crystallographic basis of the same austenite crystal, explicitly u A = u / A A .

3.2. The Transformation Matrices

This section briefly recalls the three types of transformation matrices associated with a displacive transformation from a parent austenite crystal (A) to a daughter martensite crystal (M). More details can be found in Ref. [69]. The lattice distortion is assumed to be linear; it takes the form of an active matrix F A . Any direction u A is transformed by distortion into a new direction u A = F A   u A . The distortion matrix F c A is usually expressed in the usual crystallographic basis of the parent phase c A = ( a A , b A , c A ) ; it is given by F c A = [ c A c A ] = ( a A , b A , c A ) by writing the coordinates of the three vectors a A , b A   and   c A in the column in the basis c A . The “c” index, meaning the crystallographic basis, is not systematically mentioned; the distortion matrix is often simply noted as F A . Instead of using the crystallographic basis, one may prefer using an orthonormal basis # linked to c A by the structure tensor S A . In this basis, the distortion matrix is F # A = S A   F c A   S A 1 . It can be written by polar decomposition as the product of a rotation matrix and a pure stretch (i.e., symmetric) matrix F # A = Q # A   U # A . The PTMC mainly uses the stretch matrices U # A in the calculations.
The misorientation between the austenite crystal and one of the martensite variants M is given by the coordinate transformation matrix T A M = [ c A c A ] . To make it shorter, this matrix is simply called an orientation matrix. It is a passive matrix that changes the coordinates of a fixed vector u between the parent and daughter bases as follows: u / A = T A M   u / M . By using the structure tensor, is possible to replace T A M by a rotation matrix [69].
The correspondence matrix C M A gives the coordinates in the daughter basis c M of the images by distortion of the parent basis vectors a A ,   b A   and   c A . Explicitly, C M A = ( a / M A ,   b / M A , c / M A   ) . Any direction u becomes, after lattice distortion, a direction u that, when written in the crystallographic martensite basis, is   u / M = C M A   u / A . Since a crystallographic direction of austenite becomes a crystallographic direction of martensite, the correspondence matrix is made up of simple integers or rational numbers. It can be shown that   T A M = ( T M A ) 1 and   C A M = ( C M A ) 1 . The three transformation matrices are linked by the equation C M A = T M A   F A [69].
To simplify the notation, we will sometimes write F for F A , T for T A M and C for C A M . Let us explain again their physical meaning. The distortion matrix F encodes the way the lattice is distorted. One could experimentally measure some of its components by observing before and after the transformation the deviation of a scratch at the surface of a polished surface using optical or electron microscopy or the formation of a relief using interferometric microscopy or by measuring the strain field around an isolated martensite variant using EBSD or advanced X-ray microdiffraction techniques. The orientation matrix T encodes the orientation of the martensite, but its knowledge tells nothing about the displacive or diffusive character of the transformation. Experimentally, it can be deduced from the Euler angles measured in the EBSD maps. The correspondence matrix C contains the most important information on the mechanism, because it tells how the crystallographic directions and, thus, the atomic bonding along these directions are transformed. To our knowledge, there is no experimental method in metallurgy to measure the correspondence; it is generally “guessed” by considering the crystallographic structures of the parent and martensite, as Bain did in 1924 for the fcc to bcc transformation [29] (actually, Bain did not explicitly write the correspondence; this was done by Jaswon and Wheeler in 1948 [70]). Note that the correspondence matrix does not depend on the exact value of the metrics of the phases. It is important to keep in mind that the transformation matrices F , T and C are written in the basis c and that they “work” on the directions, i.e., in the direct space. To know how they act on the planes, one must express them in the reciprocal space. They are simply F * = F t , T * = T t   and   C * = C t .

3.3. The Algebraic Structure of the Variants with Their Operators

This section is concisely written to give a global picture of the structure of variants without bogging down the reader with the details of group theory. The algorithms to build simple and double cosets are explained in simple terms in Appendix A.
The variants can be defined by coset decomposition. Let us explain it for the orientations. For an orientation matrix T A M of a specific variant M = M 1 , it can be shown that the symmetries that are common to austenite and martensite form a set T A = G A   T A M G M ( T A M ) 1 , which is the intersection between the point group of austenite G A and the point group of martensite G M . Geometrically, T A constitutes the parent and daughter symmetry elements that are in coincidence. It is a subgroup of G A called the “intersection group”. The orientation variants are defined by the left cosets M i = g i A T A , and their orientations are T A M i = g i A   T A   T A M , as detailed in Ref. [57]. It is implicitly assumed that g 1 A is the identity. The matrices that belong to g 1 A T A define the orientation variant M 1 , the matrices that belong to g 2 A T A with g 2 A T A define the orientation variant M 2 , etc. The number of orientation variants of martensite is given by Lagrange’s formula [57]: N T M = | G A | | T A | . Each orientation variant can be represented by one matrix arbitrarily chosen in its coset. It is often written that the number of variants is | G A | | G M | , but this formula is, in general, incorrect, because the type of variants (correspondence, orientation, distortion and stretch) is not specified and because the intersection group is not necessarily an isomorph to the martensite point group [69].
The specific pairs of variants ( M i , M j ) classified according to their misorientations are given by the sets of matrices T M i M j = T M i A   T A M j = ( T A M ) 1   T A   g k A   T A   T A M , where g k A = ( g i A ) 1 g j A . We call “orientation operator” a double coset, O k = T A   g k A   T A , and the associated set of misorientations is ( T A M ) 1 O k T A M . The number of operators is given by Burnside’s formula [57]. Each type of misorientation can be represented by one matrix arbitrarily chosen in its double coset.
The same definition is used to define the correspondence variants and the correspondence operators, and the only difference is that C A M should be used in place of T A M , as detailed in Ref. [69]. The different inter-correspondences C M i M j are given by the double cosets ( C A M ) 1   C A   g k A   C A   C A M . Each type of inter-correspondence can be represented by one matrix arbitrarily chosen in its double coset. When applied to a direction of M j , C M i M j gives a direction of M i , such that both directions are inherited from the same parent direction.
The distortion variants F i A are not defined as for the orientation and correspondence variants. The distortion variants are the distinct matrices g i A   F A   ( g i A ) 1 . They result from the conjugacy action of G A on F A . The number of distortion variants is given by the orbit-stabilizer theorem. It is N F M = | G A | | F A | , where F A is a subgroup of G A called the stabilizer of F A . There is another way to define the distortion variants that is closer to that used for the orientation and correspondence variants. It is based on the intersection group F A = G A   F A   G A F A 1 and coset decomposition. These variants are actually “distorted-shape” variants, but the distortion variants and the distorted-shape variants are identical when F A is close to identity, as in shape memory alloys [69].
The different correspondence, orientation and distortion matrices of a variant M i are still linked by the relation F i A = T A M i ( C A M i ) 1 . The reader can check it with F i A = g i A   F A   ( g i A ) 1 , T A M i = g i A   T A M i and C A M i = g i A   C A M .

4. The Main Principles of the Correspondence Theory

4.1. Compatibility by Symmetry Preservation

In the previous section, the orientation variants were defined by assuming a unique “natural” OR or, equivalently, since the correspondence is known, a unique “natural” distortion matrix. No attention has been paid yet in to the compatibility. Compatibility is required, because experience has shown that there are no holes or cracks between the variants. As introduced in Ref. [24], the compatibility between the variants induces new “closing-gap” ORs and continuous rotation fields. Let us explain it with the simple 2D example of a (square → parallelogram) transformation. It is assumed that the natural OR is that for which the a-axes of the two phases are parallel (Figure 1a). There are four orientation variants, four correspondence variants and four distortion variants. This result is quite obvious here, but it can be proven rigorously by coset decomposition (Appendix B). The distortion variants are, for the moment, clearly incompatible; there is no direction or plane that can be preserved between two variants, F i A and F j A . The PTMC builds the compatibility equations and determines the twins by using the stretch part of these matrices; the rotations extracted by polar decomposition are considered as “free” parameters. This method is a priori sensitive to the metrics of the phases, and it is not clear why some twins are generic and others are not. Actually, the PTMC compatibility equations are not required to determine the twins. Indeed, one can notice that the compatibility can be obtained along the mirror planes of the austenite with the help of continuous rotation fields. If the variants are linked by the vertical mirror symmetry of austenite, such as the variants M 1 and M 3 , the mirror plane p A = m y A can be maintained invariant thanks to a small rotation field (Figure 1b). If the variants are linked by a diagonal mirror symmetry of austenite, such as the variants M 1 and M 2 , the mirror plane p A = m x y A can be maintained invariant thanks to another rotation field (Figure 1c). The main reason is that, although the distortions F i A and F j A associated with the natural OR transform the mirror plane p A into distinct planes ( F i A ) * p A ( F j A ) * p A , this plane is transformed similarly by the correspondence variants C A M i g i A   C A   C A M and C A M j g j A   C A   C A M , i.e., ( C M i A ) * p A = ± ( C M j A ) * p A . The sign ± depends actually on whether the symmetry is a 180° rotation or a reflection, but this makes no difference in this 2D case because of the centro-symmetry of the lattices. This example shows that, for a pair of variants linked by a parent mirror symmetry, two opposite rotation fields are sufficient to maintain locally the compatibility at this plane, and this plane becomes the junction plane between the variants. These rotation fields have small amplitudes, because the lattice distortion is quite close to identity in shape memory alloys.
The fact that the compatibilities directly result from the parent symmetry elements allows us to incorporate them into the algebraic structure of the variants and operators recalled in Section 3.3. In most of the cases, as in the 2D example of Figure 1 or with NiTi alloys (see next sections), the orientation double cosets (misorientations) and the correspondence double cosets (inter-correspondences) are the same; this allows a one-to-one relation between the misorientations and the junction planes, as it will be explained in Section 4.2 and Section 4.3.

4.2. Junction Planes for Variants Linked by a Parent Reflection Symmetry

Let us consider two variants M 1 and M 2 joined by a reflection symmetry of austenite m A . For cubic austenite, the mirror planes are {100} and {110}. From the definition, the variant M 1 = M is the reference variant.; the correspondence matrix C M A is given relative to this variant. We note p A as the mirror plane of the symmetry m A between M 1 and M 2 . As explained in the previous section, this plane should be split into two nonparallel planes by the natural distortions of the variants, but it is not so. Indeed, p A is transformed by correspondence into the same martensitic plane p M , i.e., its Millers indices are the same in M 1 and M 2 , and the parallelism between p A and p M is maintained thanks to the additional rotation gradient. Explicitly, the junction plane is
p M = ( C M A ) * p A = ( C A M ) t p A
We note p A C p M , where C means transformed by correspondence into. The plane p M becomes a mirror plane between the variants M 1 = M and M 2 , and M 1 and M 2 are slightly rotated toward each other. The variants are thus linked by a type I twin with K 1 = p M . The “closing-gap” orientation relationship allowed by the rotation gradient is such that p A p M 1 p M 2 . The junction plane is also p A   p M 1   p M 2 .
The nature of the twin is thus completely defined by the parent symmetry and the correspondence matrix. The twin is of type I, and the mirror plane is K 1 = p M given by Equation (1). It is possible to calculate the inter-correspondence matrix of the twin. Indeed, the correspondence variants are left cosets C A M 1 = I   C A M = C A M and C A M 2 = m A   C A M . The inter-correspondence matrix is thus C M 1 M 2 = C M 1 A   C A M 2 = C M A   m A   C A M = ( C A M ) 1   m A   C A M . To simplify the notation, we note it as C int . It is easy to show that ( C int ) 2 = I . The twin is generic, because the rational twin plane K1 given by Equation (1) does not depend on the metric. Only the shear amplitude and the irrational shear direction depend on the metric, as detailed as follows. The shear amplitude s is immediately given by Bevis and Crocker’s formula [37]; its square is
s 2 = T r   ( C int t     C int   1 ) 3
where is the metric tensor of the martensite phase. The shear direction d M 1 M 2 can be calculated from the correspondence matrix without explicitly using the distortion matrix. As shown in Figure 2a, although the plane p A is not tilted during its transformation into p M , i.e., p A p M , its normal n A is sheared into opposite directions for the two variants. It becomes n M 1 = C M 1 A n A for the variant M 1 and n M 2 = C M 2 A n A for the variant M 2 . The shear direction is d M 1 M 2 = n M 2 n M 1 , but this expression is difficult to calculate, because the vectors are written in their own respective basis. Instead of using n M 1 and n M 2 , let us consider n M 1 the normal to the plane p M . It is given by
n M 1 = 1 p M = 1 ( C M A ) * p A
This direction is originated from a parent direction represented in dashed blue in Figure 2a. The knowledge of this vector is not important; we just need to know that it is transformed into a direction of the variant M 2 , n M 2 . The vectors n M 1 and n M 2 are in inter-correspondence, i.e., n M 2 = C int   n M 1 . The shear direction is d M = d M 1 M 2 = n M 2 n M 1 . Here, again, one has to keep in mind that, in this equation, n M 2 is written in c M 2 and n M 1 in c M 1 . Let us write all the terms in c M 2 . The first term is not changed, but n M 1 is changed into its opposite, n M 1 / M 1 = n M 1 / M 2 because n M 1 is normal to the mirror plane. It comes that the shear direction d M 1 M 2 written in c M 2 is d M = C int   n M 1 + n M 1 = ( C int + I )   n M 1 . Since it belongs to the mirror plane between the variants M 1 and M 2 , its coordinates are the same in c M 1 . They are rewritten in c M 1 by changing again n M 1 into n M 1 , and we get
d M = ( C int + I )   n M 1
If we assume that this direction is directly inherited from a direction of the parent phase by correspondence without any rotation, we get d A = C A M   d M and d A d M . The closing-gap OR and the twinning characteristics are now completely known; they are
{   p A   p M 1   p M 2   ( r a t i o n a l   s h e a r   p l a n e   K 1 ) d A   d M 1   d M 2   ( i r r a t i o n a l   s h e a r   d i r e c t i o n   η 1 )
The “fictive” simple shear distortion matrix S between the two variants is given by S = I + s   d ˜   p   ˜ t = I + s   d ˜ p   ˜ , where d = d M and p = p M , both vector being normalized in their own spaces (direction and reciprocal, respectively), i.e., d ˜ 2 = d ˜ t     d ˜ = 1 and p ˜ 2 = p ˜ t   * p ˜ = 1 . The type I twin elements are fully determined; they are K 1 = p M and η 1 = d M . It is, however, important to keep in mind that the mechanism of transformation from a variant into another one is not exactly a simple shear, because both variants are created from the parent phase. It is often assumed that variant reorientation is mediated by a simple shear that reduces one variant to the profit of another variant that is better orientated in the stress field, but actually, nothing proves that a simple shear mechanism is implied in this process. As it will be discussed in Section 6, we think it is more probably a double-step mechanism, M 1 A M 2 .

4.3. Junction Planes for Variants Linked by a Parent 180° Rotation Symmetry

Let us consider now two variants M 1 and M 2 joined by a 180° rotation symmetry of austenite R π A . For cubic austenite, the rotation axes are <100> and <110>. The variant M 1 is the reference variant. We note u A as the rotation axis of R π A . This axis should be split into two nonparallel axes by the natural distortions of M 1 and M 2 , but it is not so, and u A is maintained invariant thanks to an additional rotation gradient. It is transformed by correspondence into the same martensitic direction u M , i.e., its indices are the same in M 1 and M 2 :
u M = C M A u A
We note u A C u M , where C means transformed by correspondence into. The direction u M becomes a 180° rotation axis between the variants that are slightly rotated toward each other. The variants are thus linked by a type II twin with η 2 = u M . The “closing-gap” OR allowed by the rotation gradient is such that u A u M 1 u M 2 . According to the classical twinning theory, the junction should be the irrational plane p = K 2 . Let us explain how it can be calculated without using the distortion matrices. The inter-correspondence matrix is C int = C M 1 M 2 = C M A   R π A   C A M = ( C A M ) 1   R π A   C A M . We have shown in Ref. [71] that type II twins in the direct space are type I twins in the reciprocal space. They can be calculated exactly as type I twins by exchanging the shear plane and the shear direction p d and the matrices of the direct space with those of the reciprocal space * and C C * . Bevis and Crocker’s formula can then be directly applied after these exchanges. The square of the shear value of s * is s * 2 = T r   ( C int * t   *   C int *   * 1 ) 3 . Since C int 1 = C int , it is
s * 2 = T r   ( C int   1   C int t   ) 3
Since the trace does not depend on the order of the matrices in the product, we can immediately check that s * 2 = s 2 . To calculate the shear plane, instead of considering the inter-corresponding planes p M 1 and p M 2 inherited from p A , we directly use p M 1 , the plane normal to the direction u M . It is given by a vector of the reciprocal space:
p M 1 =   u M =   C M A u A
This plane is inherited from a parent plane whose normal vector is represented in dashed blue in Figure 2b. The knowledge of this vector is not important; we just need to know that it is transformed into a reciprocal vector of the variant M 2 , p M 2 . The vectors p M 1 and p M 2 are in inter-correspondence, i.e., p M 2 = C i n t *   p M 1 . The reciprocal shear direction is d M 1 M 2 * = p M 2 p M 1 . Here, again, one has to keep in mind that, in this equation, p M 2 is written in ( c M 2 ) * and p M 1 in ( c M 1 ) * . Let us write all the terms in ( c M 2 ) * . The first term is not changed. The second term p M 1 is not changed either when written into the basis of M 2 , because it represents a plane that is normal to the rotation axis. It becomes that the reciprocal shear direction d M 1 M 2 * written in ( c M 2 ) * is d M 1 M 2 * = C int *   p M 1 p M 1 = ( C int * E )   p M 1 . This vector can be rewritten in ( c M 1 ) * by taking its opposite. This reciprocal vector is the shear plane in the direct space, which is also the junction plane j p M between the two variants. We can write it as
j p M = d M 1 M 2 * = ( C int * I )   p M 1
The type II twin elements are fully determined; they are η 2 = u M and K 2 = j p M .
If we assume that the junction plane is directly inherited from a plane of the parent phase by correspondence, we get j p A = ( C A M ) *   j p M and j p A j p M . The closing-gap OR and the twinning characteristics are now completely known; they are
{   u A   u M 1   u M 2   ( r a t i o n a l   s h e a r   d i r e c t i o n   η 2 ) j p A   j p M 1   j p M 2 ( i r r a t i o n a l   s h e a r   p l a n e   K 2 )
We have assumed, as in the usual twinning theory, that the junction plane for type II twins is the irrational plane K2. Actually, as it will be discussed in the next section, a different assumption is also conceivable: the junction planes of type II twins could be low-index rational planes that are not fully invariant but slightly distorted. We called them “weak planes”. In most of the cases, the calculations show that the weak planes are very close to the K2 planes, and it is difficult from an experimental point of view to distinguish the two hypotheses (fully invariant irrational K2 plane or slightly distorted but rational weak plane). We will see, however, in the next section, that the junction planes associated with the 90° and 120° rotational symmetries cannot be calculated with the usual twining theory, whereas they can be explained and predicted with the concept of a weak plane.

4.4. Junction Planes for Variants Not Linked by a Two-Fold Parent Symmetry

What is the junction plane between two variants that are linked by a symmetry that is not an order of two such as those inherited from the four-fold rotations around the <100> axes and three-fold rotations around the <111> axes of the parent B2 cubic phase? Bevis and Crocker’s formula of twinning cannot be applied anymore, because the correspondence matrix does not verify the property ( C 12 ) 2 = I on which the theory is based. The PTMC compatibility equations based on the stretch matrices have no solution. These twins are called “non-generic non-conventional” by Pitteri and Zanzotto [47,48]. It is generally inferred that there is no junction plane for these specific misorientations [2,49] or that the junctions may exist only for some specific values of the stretch matrix [47]. We disagree on this point. Our experimental and theoretical previous works on martensite in steels made us conclude that the interface plane between austenite and isolated martensite is not necessarily fully invariant [62], and some deformation twins in magnesium with interfaces that have different and non-equivalent Miller indices have been already observed [72]. The idea that the interface (habit plane or junction plane) should be fully invariant has a long history that comes from the concept of simple shear and twinning dislocations, but it is highly questionable [73]. The concepts of weak plane and weak twins directly result from our conviction that an interface plane is not necessarily fully invariant. The details are reported in a separate paper [74]; we just give here the main idea. We call a “weak plane” an interface between two crystals (1 and 2) that (i) has rational Miller indices that are not necessarily the same or equivalent when the plane is indexed in crystal 1 and in crystal 2, i.e., p 1 p 2 but p 1 = ( h 1   k 1   l 1 ) and p 2 = ( h 2   k 2   l 2 ) can be different, and (ii) is such that a slight intraplanar distortion is sufficient to transform one into the other, ( h 1   k 1   l 1 )   ( h 2   k 2   l 2 ) . An axial weak plane is a weak plane that contains a crystallographic direction that is invariant (unrotated and undistorted). A schematic representation of an axial weak plane is shown in Figure 3. The plane is constituted by two rational directions in crystal 1, u 1 and v 1 , and two directions in crystal 2, u 2 and v 2 . In this figure, the twin axis is the invariant direction u 1 = u 2 . The rational directions v 1 and v 2 have close norms and make nearly the same angle with the axis u 1 of the twin. The weak plane is ( u 1 , v 1 ) that is parallel to ( u 2 , v 2 ) ; its Miller indices are the coordinates of the reciprocal vector p 1 = ( u 1 × v 1 ) in c 1 , and of p 2 = ( u 2 × v 2 ) in c 2 . We note the weak plane is p 1 p 2 .
The exact way crystal 1 can be transformed into crystal 2 is not a simple shear, as it was described in Ref. [74]. A computer program that calculates the axial weak planes was written in Python and incorporated in GenOVa [58]. When applied to the deformation twins in magnesium, it “predicts” an axial weak twin associated with a weak plane ( 0001 ) 1 ( 01 1 ¯ 0 ) 2 around the a-axis (rotation angle 90°) that is a classical basal/prismatic interface observed locally along (86°, a) extension twins [75]. It also predicts the unconventional twin with a weak plane ( 3 1 ¯ 2 ¯ 2 ) 1 ( 0 1 ¯ 12 ) 2 around the axis 22 4 ¯ 3 (rotation angle 58°), which was recently observed by EBSD [76]. GenOVa calculates all the characteristics of the weak twins, including the generalized shear value and the distortion, orientation and correspondence matrices [74].
Let us come back now to the case of two variants linked by a parent symmetry that is not two-fold, for example, for the three-fold rotation symmetry R 2 π / 3 A . The rotation axis u A of this symmetry should be tilted by the natural lattice distortion, but we assume that it is actually maintained unrotated thanks to the additional rotation gradient. It is transformed into u M = C M A u A . The possible weak planes and associated weak twins of the axis u M are calculated by GenOVa, and only those whose correspondence matrix are equal to the inter-correspondence matrix between the two variants C int = C M 1 M 2 = ( C A M ) 1   R 2 π / 3 A   C A M are considered. Depending on the tolerances used to determine the directions v 1 and v 2 (Figure 3), different weak planes are predicted. Those associated with the lowest tolerances are ranked in the first positions. It will be shown that, in the case of NiTi alloys, the predicted “weak” junction planes agree well with the experimental TKD maps.
In summary, the correspondence theory assumes that there is a natural distortion and a natural OR and that the compatibility between the variants is obtained by additional closing-gap rotations of small amplitudes. The variants are linked by the symmetries of the parent phase. In the case of austenite reflection, the variants are linked by the rational mirror plane inherited from the austenitic one. The variants are in type I twin relation, and the junction plane is the plane given by Equation (1). In the case of austenite 180° rotation, the variants are linked by the rational rotation axis inherited from the austenitic one by Equation (6). They are in type II twin relation, and the junction plane is the invariant irrational plane K2 or a weak plane close to it. In the case of an austenite symmetry that is not a 180° rotation (e.g., 90° or 120°), the variants are linked by the rotation axis inherited by correspondence from the parent rotation axis; the twins are nonconventional, and we infer that the junctions are “weak planes”.
All the twins can be predicted from C A M , the correspondence matrix, and , the metric tensor of the martensite phase, without calculating the stretch matrices. The rational “shear” plane K1 of type I twins is the correspondence image of the austenite mirror plane; the rational “shear” direction η 2 of the type II twins and of the weak twins is the correspondence image of the austenite rotation axis. They depend only on the correspondence; they are independent of the metrics. Only the irrational values of the twins, i.e., the “shear” direction of the type I twins and the “shear” plane of the type II twins and the weak plane of the weak twins, require the knowledge of . If one considers only the rational elements of the twins, all the twins are actually “generic”.
The variants of the B2 → B19′ transformation in NiTi alloys and their junction planes will be studied in detail in Section 5. We will show that the calculations with CT give exactly the same results as with PTMC, at least for the type I and II twins. The advantages of the CT over PTMC will appear clearly. The variants and the different pairs of variants will be rigorously determined from coset and double coset decomposition from the early beginning. The nature of the twins (type I, type II or weak) will be immediately recognized nearly without numerical calculations. The predictions will be compared to some experimental TKD maps in Section 5.4. It will be the first time that TKD is used to study the junction planes.

5. Application of the Correspondence Theory to the Junction Planes in NiTi Alloys

5.1. The Variants and the Operators

For the calculations, we used the same lattice parameters as reported for NiTi alloys in most of the PTMC studies, i.e., a = 3.015   for the parent B2 “austenite” and a = 2.889   ,   b = 4.120     and   c = 4.622   ,   β = 96.8 ° for the monoclinic B19′ “martensite”. The reader can keep in mind some symmetry rules in the B19′ phase. For example, the planes ( h ¯ k l ) B 19 and ( h k l ¯ ) B 19 are parallel; they are equivalent by symmetry to the plane ( h ¯ k ¯ l ) B 19 , but they are not equivalent to ( h k l ) B 19 . The same rule on the signs of the indices is valid for the directions.
The metric tensor of B19′ calculated from the lattice parameters is
= ( 8.3463 0 1.5810 0 16.9744 0 1.5810 0 21.3628 )
The correspondence matrix according to Otsuka and Ren’s model [3] is
C B 2 B 19 = ( 0 1 1 0 1 1 1 0 0 )
The natural OR is that for which the dense plane and dense direction of two phases are parallel, i.e., (010)B19′//(110)B2 and [101]B19′//[ 1 ¯ 11]B2. This OR was noted as the AQ in Ref. [24]. The distortion and orientation matrices associated with this OR are
F B 2 = ( 1.0251 0.0589 0.1297 0.0589 1.0251 0.1297 0.0016 0.0016 0.9511 ) T B 2 B 19 = ( 0.1297 0.9663 1.0840 0.1297 0.9663 1.0840 0.9511 0 0.0032 )
The plane (1 1 ¯ 2)B2//(10 1 ¯ )B19′ is unrotated by the distortion F B 2 ; it is the habit plane of the large martensite laths mapped in EBSD in Ref. [24]. For simplification, we will note C = C B 2 B 19 , T = T B 2 B 19 and   F = F B 2 . The point group of the parent B2 phase is G B 2 constituted of the 48 symmetry matrices of the cube. The point group of martensite is G B 19 = { I ,   I ¯ ,   R π , [ 010 ] B 19 ,   m ( 010 ) B 19 } .
The intersection groups for the correspondence and orientation variants are C B 2 = G B 2   C   G B 19 C 1 and T B 2 = G B 2   T   G B 19 T 1 , respectively. The intersection group for the distorted-shaped variants is F B 2 = G B 2   F   G B 2 F 1 [69]. The calculations show that these three intersection groups are equal, whatever the relation, C B 2 = T B 2 = F B 2 . It is a subgroup of G B 2 that will be simply noted B 2 . It is constituted of the four symmetries matrices “common” to the daughter and parent phases and expressed in the parent crystallographic basis, B 2 = { I ,   I ¯ , R π , [ 110 ] B 2 , m ( 110 ) B 2 } , where
I = ( 1 0 0 0 1 0 0 0 1 ) ,   I ¯ = ( 1 0 0 0 1 0 0 0 1 ) ,   R π , [ 110 ] B 2 = ( 0 1 0 1 0 0 0 0 1 ) ,   m ( 110 ) B 2 = ( 0 1 0 1 0 0 0 0 1 )
It is the isomorph to the daughter point group G B 19 , but this is not a general rule.
The correspondence and orientation variants are the simple left cosets multiplied at the right by the transformation matrix C and T, respectively. Since, for the B2 → B19′ transformation, the left cosets and double cosets do not depend on the transformation matrix, we will just use the terms “variant” and “operator” without specifying their type (correspondence, orientation or distortion). The number of variants is 48/4 = 12 by Lagrange’s formula. The calculation shows that there are seven operators, which agrees with Burnside’s formula. The variants, operators and their algebraic representations in cosets are summarized in Table 1. Since the intersection groups are the same, there are one-to-one relations between the correspondence, orientation and distortion variants and one-to-one relations between the correspondence, orientation and distortion operators. Any specific misorientation, i.e., a double coset, T 1   B 2   g j B 2   B 2   T , can be associated a specific inter-correspondence relation, i.e., a double-coset C 1 B 2   g j B 2 B 2 C . This will allow us to predict the junction planes between variants by the unique knowledge of their misorientation.

5.2. Prediction of the Junction Planes and Closing-Gap ORs

The 48 symmetry matrices g j B 2 of the parent B2 phase are partitioned in Table 2 into double cosets B 2 g j B 2 B 2 . The table also indicates for each operator the disorientation between the B19′ variants. The disorientation is the rotation in the set of equivalent rotations and roto-inversions T 1   T B 2 g j B 2 T B 2   T that has the minimum rotation angle. It is the rotation that is given by commercial software when drawing a line between two pixels or grains in an EBSD or TKD map.
The neutral operator O 0 leaves each variant invariant. It is associated with the double coset formed by the matrices g j B 2 B 2 . There is thus no junction plane associated with this operator.
The operator O 2 is characterized by a disorientation between the variants that is a rotation of 180° around the axis [001]B19′. The double coset contains two reflection symmetries: m ( 001 ) B 2 and m ( 1 1 ¯ 0 ) B 2 . The associated mirror planes are transformed by correspondence as follows: ( 001 ) B 2 C * ( 100 ) B 19 and ( 1 1 ¯ 0 ) B 2 C * ( 001 ) B 19 . The predicted K1 junction planes are thus the planes p = ( 100 ) B 19 or p = ( 001 ) B 19 . The first correspondence ( 001 ) B 2 C * ( 100 ) B 19 induces a slight rotation that maintains the contact between the ( 001 ) B 2 plane and the ( 100 ) B 19 plane of the variants in the pairs, i.e., ( 001 ) B 2 ( 100 ) B 19 ( 100 ) B 19 . The notation ( 100 ) B 19 ( 100 ) B 19 is short and means that the (100) plane of one variant is parallel to the (100) plane of the other variant. The second correspondence ( 1 1 ¯ 0 ) B 2 C * ( 001 ) B 19 induces another closing-gap OR, for which ( 1 1 ¯ 0 ) B 2 ( 001 ) B 19 ( 001 ) B 19 . The inter-correspondence matrices deduced from the reflections m ( 001 ) B 2 and m ( 1 1 ¯ 0 ) B 2 are
C int B 19 = ( 1 0 0 0 1 0 0 0 1 )   and   C int B 19 = ( 1 0 0 0 1 0 0 0 1 ) ,   respectively .
These matrices allow a direct calculation of the shear amplitude s and shear direction d by using Equations (2)–(4). With the twin on p = ( 100 ) B 19 , we get d = [ 001 ] B 19 and s = 0.2385 . The shear direction d is in correspondence with [ 1 ¯ 10 ] B 2 . With the twin on p = ( 001 ) B 19 , we get d = [ 100 ] B 19 and s = 0.2389 . The shear direction d is in correspondence with [ 001 ] B 2 . Since the shear directions are rational, the two type I twins are also of type II. They are compound twins. The closing-gap OR associated with the mirror planes m ( 001 ) B 2 is ( 100 ) B 19 ( 001 ) B 2 and [ 001 ] B 19 [ 1 ¯ 10 ] B 2 , which is the OR A observed experimentally in Ref. [24], and that associated with m ( 1 1 ¯ 0 ) B 2 is ( 001 ) B 19 ( 1 1 ¯ 0 ) B 2 and [ 100 ] B 19 [ 001 ] B 2 , which is the OR C observed experimentally in Ref. [24].
The operator O 2 also contains the 180° rotation symmetries R π , [ 001 ] B 2 and R π , [ 1 1 ¯ 0 ] B 2 ; the associated twins are the same compound twins as those determined previously.
The operator O 4 is characterized by a misorientation between the variants that is a rotation of 120° around the axis ~[17,0,16]B19′. It contains two reflection symmetries: m ( 01 1 ¯ ) B 2 and m ( 101 ) B 2 . The junction planes predicted by the correspondence theory are thus ( 01 1 ¯ ) B 2 C * ( 1 ¯ 11 ) B 19 and ( 101 ) B 2 C * ( 11 1 ¯ ) B 19 . The two planes ( 1 ¯ 11 ) B 19 and ( 11 1 ¯ ) B 19 are equivalent by the B19′ symmetries. Let us consider only the first correspondence and the associated junction plane p = ( 1 ¯ 11 ) B 19 between the variants. The inter-correspondence matrices deduced from the reflection m ( 01 1 ¯ ) B 2 are
C int B 19 = ( 0 1 1 1 2 1 2 1 2 1 2 1 2 1 2 )
This matrix allows a direct calculation of the shear direction d and amplitude s. From Equations (2)–(4), we get d = ~   [ 7 ,   6 ¯ , 13 ] B 19 and s = 0.3096 . The correspondence ( 01 1 ¯ ) B 2 C * ( 1 ¯ 11 ) B 19 induces a slight rotation that maintains the contact between the ( 01 1 ¯ ) B 2 plane of one variant with the ( 1 ¯ 11 ) B 19 plane of the other variant, i.e., ( 01 1 ¯ ) B 2 ( 1 ¯ 11 ) B 19 ( 1 ¯ 11 ) B 19 . The martensite direction d is in correspondence with the parent direction by d = ~   [ 7 ,   6 ¯ , 13 ] B 19 C ~   [ 19 ¯ , 7 , 7 ] B 2 . It is assumed that this direction is not rotated by the distortion. The closing-gap OR is thus ( 1 ¯ 11 ) B 19   ( 01 1 ¯ ) B 2 and [ 7 ,   6 ¯ , 13 ] B 19 [ 19 ¯ , 7 , 7 ] B 2 .
The operator O 4 also contains the 180° rotation symmetries R π , [ 01 1 ¯ ] B 2 and R π , [ 101 ] B 2 . By correspondence, the rotation axes become [ 01 1 ¯ ] B 2 C [ 2 ¯ 11 ] B 19 and [ 101 ] B 2 C [ 21 1 ¯ ] B 19 . The two directions [ 2 ¯ 11 ] B 19 and [ 21 1 ¯ ] B 19 are equivalent by the B19′ symmetries. Let us consider only the first correspondence and the associated junction direction u = [ 2 ¯ 11 ] B 19 between the variants. The inter-correspondence matrix deduced from the rotation R π , [ 01 1 ¯ ] B 2 and written in the reciprocal space is
( C int B 19 ) * = ( 0 1 2 1 2 1 1 2 1 2 1 1 2 1 2 )
This matrix allows a direct calculation of the junction (shear) plane j p and amplitude s * . From Equations (7)–(9), we get j p   ~   ( 1 , 2 ¯ , 4 ) B 19 and s * = 0.3096 . The martensite plane j p is in correspondence with a parent plan by ( 1 , 2 ¯ , 4 ) B 19 C * ( 3 ¯ , 1 , 1 ) B 2 . We assume that this plane is not tilted by the distortion; the closing-gap OR is thus [ 2 ¯ 11 ] B 19 [ 01 1 ¯ ] B 2 and ( 1 , 2 ¯ , 4 ) B 19 ( 3 ¯ , 1 , 1 ) B 2 .
The operator O 5 is characterized by a misorientation between the variants that is a rotation of 120° around the axis ~ [ 40 3 ¯ ] B 19 . It contains two reflection symmetries: m ( 011 ) B 2 and m ( 1 ¯ 01 ) B 2 . The junction planes predicted by the correspondence theory are thus ( 011 ) B 2 C * ( 111 ) B 19 and ( 1 ¯ 01 ) B 2 C * ( 1 1 ¯ 1 ) B 19 . The two planes ( 111 ) B 19 and ( 1 1 ¯ 1 ) B 19 are equivalent by the B19′ symmetries. Let us consider only the first correspondence and the associated junction plane p = ( 111 ) B 19 between the variants. The correspondence matrices deduced from the reflection m ( 011 ) B 2 are
C int B 19 = ( 0 1 1 1 2 1 2 1 2 1 2 1 2 1 2 )
This matrix allows a direct calculation of the shear direction d and amplitude s. From Equations (2)–(4), we get d = ~ [ 62 ¯ , 21 , 41 ] B 19 and s = 0.1422 . The martensite direction d is in correspondence with a parent direction by [ 62 ¯ , 21 , 41 ] B 19 C [ 10 ¯ , 31 , 31 ¯ ] B 2 . We assume that this direction is not rotated by the distortion. The closing OR is thus ( 111 ) B 19 ( 011 ) B 2 and [ 62 ¯ , 21 , 41 ] B 19 [ 10 ¯ , 31 , 31 ¯ ] B 2 .
The operator O 5 also contains the 180° rotation symmetries R π , [ 011 ] B 2 and R π , [ 1 ¯ 01 ] B 2 . By correspondence, the rotation axes become [ 011 ] B 2 C [ 211 ] B 19 and [ 101 ] B 2 C [ 2 1 ¯ 1 ] B 19 . The two directions [ 211 ] B 19 and [ 2 1 ¯ 1 ] B 19 are equivalent by the B19′ symmetries. Let us consider only the first correspondence and the associated junction direction u = [ 211 ] B 19 between the variants. The inter-correspondence matrix deduced from the rotation R π , [ 01 1 ¯ ] B 2 and written in the reciprocal space is
( C int B 19 ) * = ( 0 1 2 1 2 1 1 2 1 2 1 1 2 1 2 )
From Equations (7)–(9), we get j p   ~ ( 2 ¯ , 1 ,   3 ) B 19 and s * = 0.1422 . The martensite plane j p is in correspondence with a parent plane by ( 2 ¯ , 1 ,   3 ) B 19 C * ( 1 , ¯ 2 , 2 ¯ ) B 2 . The closing-gap OR is thus [ 211 ] B 19   [ 011 ] B 2 and ~ ( 2 ¯ , 1 ,   3 ) B 19 ~ ( 1 , ¯ 2 , 2 ¯ ) B 2 .
The operator O 6 is characterized by a misorientation between the variants that is a rotation of 90° around the axis ~[11,0,1]B19′. It contains two reflection symmetries: m ( 010 ) B 2 and m ( 100 ) B 2 . The junction planes predicted by the correspondence theory are ( 010 ) B 2 C * ( 011 ) B 19 or ( 100 ) B 2 C * ( 01 1 ¯ ) B 19 . The two planes ( 011 ) B 19 and ( 01 1 ¯ ) B 19 are equivalent by the B19′ symmetries. Let us consider only the first correspondence and the associated junction plane ( 011 ) B 19 between the variants. The inter-correspondence matrices deduced from the reflection m ( 010 ) B 2 is
C int B 19 = ( 1 0 0 0 0 1 0 1 0 )
From Equations (2)–(4), we get d = ~ [ 11 ¯ , 7 ¯ , 7 ] B 19 and s = 0.2804 . The martensite direction d is in correspondence with a parent direction by [ 11 ¯ , 7 ¯ , 7 ] B 19 C [ 14 ¯ , 0 ,   11 ¯ ] B 2 . The closing-gap OR is thus ( 011 ) B 19   ( 010 ) B 2 and [ 11 ¯ , 7 ¯ , 7 ] B 19 [ 14 ¯ , 0 ,   11 ¯ ] B 2 .
The operator O 6 also contains the 180° rotation symmetries R π , [ 010 ] B 2 and R π , [ 100 ] B 2 . By correspondence, the rotation axes become [ 010 ] B 2 C [ 011 ] B 19 and [ 101 ] B 2 C [ 01 1 ¯ ] B 19 . The two directions [ 011 ] B 19 and [ 01 1 ¯ ] B 19 are equivalent by the B19′ symmetries. Let us consider only the first correspondence and the associated junction direction u = [ 011 ] B 19 . The correspondence matrix deduced from the rotation R π , [ 010 ] B 2 and written in the reciprocal space is
( C int B 19 ) * = ( 1 0 0 0 0 1 0 1 0 )
From Equations (7)–(9), we get j p   ~   ( 3 , 4 , 4 ¯ ) B 19 and s * = 0.2804 . The martensite plane j p is in correspondence with a parent plan by ( 3 , 4 , 4 ¯ ) B 19 C * ( 4 , 0 , 3 ) B 2 . The closing-gap OR is thus [ 011 ] B 19   [ 010 ] B 2 and ~ ( 3 , 4 , 4 ¯ ) B 19 ~ ( 4 , 0 , 3 ) B 2 .
The operators O 1 and O 3 are complementary polar operators. Their double cosets do not contain any reflection or 180° rotation symmetry. Contrarily to the previous operators, there is no junction plane of type K1 or K2. The usual belief is that that there should not be junction plane at all, but, as explained in Section 4.4, we disagree, because we think that twinning cannot be reduced to a simple shear [73]. Let us explain how the concept of weak plane and weak twin can be used to establish a list of potential junction planes. First, it is considered in the operators O 1 and O 3 how the 90° rotation axes and the 120° rotation axes are transformed by correspondence: [ 100 ] B 2 C [ 011 ] B 19 and [ 111 ] B 2 C [ 110 ] B 19 . Then, GenOVa [74] is used to calculate the axis weak twins along the axes [ 011 ] B 19 and [ 110 ] B 19 . The B19′ lattice parameters are those given in Section 5.1, and the tolerances chosen for the reticular distances and angles are 2% and 2°, respectively (see Ref. [74] for the details). Only the weak twins for which the correspondence matrices are of type C int B 19 are considered. The inter-correspondence matrices C int B 19 belong to the double set C B 19 B 2   g j B 2   C B 2 B 19 , where g j B 2 are the matrices in the double cosets of O 1 and O 3 . The calculations show that there are 16 inter-correspondence matrices; they are of type C int B 19 = ε 1 ( 0 1 ε 2 1 ε 2 ε 3 1 2 ε 3 ε 4 1 2 1 2 ε 3 1 2 ε 2 ε 3 1 2 ε 2 ε 4 1 2 ε 2 ε 3 ε 4 ) with ε i = { 1 1 for i { 1 , 2 , 3 , 4 } .
Let us choose the matrix with ε i = { 1 , 1 , 1 , 1 } , and we get
C int B 19 = ( 0 1 1 1 2 1 2 1 2 1 2 1 2 1 2 )
There is only one weak twin for which the correspondence matrix is C int B 19 . It is the weak twin around the axis [ 011 ] B 19 . The associated weak plane is ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 . More specifically, the intraplanar transformation occurs by the transformation of directions [ 310 ] B 19 C int B 19 [ 1 ¯ 2 1 ¯ ] B 19 and [ 30 1 ¯ ] B 19   C int B 19 [ 1 ¯ 1 ¯ 2 ¯ ] B 19 , and the transformation outside the weak plane is such that [ 100 ] B 19   C int B 19 [ 01 1 ¯ ] B 19 . The generalized strain and twin index are sg = 0.2912 and qg = 2, and the misorientation is a rotation of 84° around [ 011 ] B 19 . With larger tolerances on the reticular distances and angles of 5% and 5°, two other [ 011 ] B 19 weak twins are also “predicted”; they are the ( 1 1 ¯ 1 ) B 19 ( 11 1 ¯ ) B 19 weak twin, for which the disorientation is a rotation of 86° around [ 011 ] B 19 , and the ( 0 1 ¯ 1 ) B 19 ( 100 ) B 19 weak twins, for which the disorientation is a rotation of 85° around [ 011 ] B 19 . The three weak twins have the same generalized strain values and twin index, but the lowest intraplanar distortion is that of the ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 weak twin. For this twin, the direction [ 011 ] B 19 is in correspondence with [ 010 ] B 2 . This parallelism gives one component of the closing-gap OR. The second component depends on which of the planes ( 1 3 ¯ 3 ) B 19 or ( 31 1 ¯ ) B 19 is considered. In the first case, the closing-gap OR is [ 011 ] B 19   [ 010 ] B 2 and ( 1 3 ¯ 3 ) B 19 ( 3 ¯ 01 ) B 2 , and in the second case, it is [ 011 ] B 19   [ 010 ] B 2 and ( 31 1 ¯ ) B 19 ( 103 ) B 2 .
The junction planes predicted by the correspondence theory are summarized in Table 3. Note that the junction planes for the type II twins are given by the irrational invariant planes K2 written by a rational approximate marked by the symbol ~ . As mentioned in Section 4.4, it is conceivable that the junction planes of type II twins are actually weak planes. We used GenOVa to calculate the rational axial weak planes associated with the 180° rotation axes of operators O 4 ,   O 5 and O 6 , and we found the rational planes ( 1 2 ¯ 4 ) B 19 , ( 2 ¯ 13 ) B 19 and ( 34 4 ¯ ) B 19 , respectively. Note that these planes are those used to write the irrational planes K2 by their rational approximates, but here, with the assumption of a weak plane, the junction planes are supposed to be exactly these rational planes. The deviations between the two planes are so low (less than 1°) that it is impossible to discriminate them from an experimental point of view. It can also be added that, for the operator O 5 , two other weak planes were found; they are ( 1 ¯ 11 ) B 19 and ( 1 ¯ 02 ) B 19 , but their intraplanar angular deviations (0.17° and 1.8°, respectively) are higher than for the plane ( 2 ¯ 13 ) B 19 which is only 0.02°.
The junction planes that are observed experimentally are written in bold in Table 3. One can wonder why some operators have a tendency to form type I twin and not type II twin junction planes or the opposite, for some other operators, whereas both type I and type II twins have the same shear amplitude. For example, why is the junction plane of the operator O 4 ( 1 ¯ 11 ) B 19 and not ~ ( 1 2 ¯ 4 ) B 19 , or why is the junctions plane of operator O 5 ~ ( 2 ¯ 13 ) B 19 and not ( 111 ) B 19 ? One can also wonder why the junction plane of operator O 2 is ( 100 ) B 19 and not ( 001 ) B 19 whereas both are conjugate compound twins. Thanks to the CT, an explanation based on closing-gaps ORs can be proposed. These ORs are summarized in the second column of Table 4. For each operator in Table 4, the deviation angle between the closing-gap ORs and the natural OR AQ given by the matrix T B 2 B 19 in Equation (13) is given. This angle is obtained with GenOVa by calculating the lowest disorientation between the 12 variants of OR AQ and the 12 or 24 variants of the closing-gap OR that is considered. It is striking that, for each operator, the type of twin that has the minimum deviation is also that one that is reported in the literature and observed in TKD (Section 5.4). We think this is not a coincidence. The CT indeed assumes the existence of a natural OR; it is thus expected that the twin is that one that makes the variants deviate as little as possible from their natural orientations, i.e., that one for which the orientation gradients and strains are as low as possible.
The disorientation histogram with the distribution of rotations axes associated with all the variants of all the closing-gap ORs of Table 4 is presented in Figure 4. It is quite similar to those already presented Ref. [24], which shows that the closing-gap ORs are quite close to the ORs A, AQ, C, CQ and I determined experimentally from EBSD. This means that the orientation gradients predicted from the CT also agree quite well with the experiments.

5.3. Comparison with PTMC

The twin characteristics and the shear amplitudes determined by the CT are now compared with those calculated by PTMC [2,49]. The result is presented in Table 5. Note that some typos in the Table 5.1 of Ref. [2] have been corrected, such as the error of the sign in the indices of the plane ( 0.7206 ,   1 ,   1 ¯ ) B 19 and inversion of the type II twins and shears between the categories C and D in Ref. [2].
It is clear that the results are the same for the type I and type II twins. Therefore, what is advantage of using CT instead of PTMC? To answer this question, it is important to recall that the PTMC only considers the distortion variants F i or the stretch variants U i , but it does not distinguish the different types of variants, such as the stretch variants, the distortion variants, the orientation variants and the correspondence variants [73]. This choice may be unfortunate, since it leads to assume that the number of “variants” (which ones?) is the number of symmetries of the parent phase divided by the number of symmetries of the daughter phase, which is correct for NiTi but not true in general (see Section 3.3). Besides this problem, PTMC calculations are unnecessarily complex, and they mask the relative roles of the symmetries and metrics. In PTMC, the metrics of the parent and daughter phases are included in each matrix U i and the symmetries are already treated by the identification of the 12 “equivalent” stretch matrices. As explained in Section 2.1, the PTMC needs to calculate the eigenvalues λ 1 , λ 2 and λ 3 and the eigenvectors e 1 , e 2 and e 3 of the matrices U j T U i T U i U j 1 for all of the 12 × 11 = 132 pairs ( U i , U j ) . Bhattacharya reported 192 twinning modes [2]. Hane and Shield (1999) could also show that the 132 pairs can be grouped into six types of pairs of variants that they called S1–S6. The CT finds the same sets, but the calculations are simple and direct. They require only the determination of the intersection group (group of common symmetries), and the decomposition of the parent B2 point group into double cosets. The matrices are simple 3 × 3 matrices made of 0, 1 and −1. In addition, since, for this transformation, the intersection group does not depend on the type of relation (orientation, distortion and correspondence), a one-to-one link between the misorientations, the “compatibility twins” and their junction plane can be established. The number of different pairs of variants found by Hane and Shield is nothing else than the double cosets (seven, if one counts the neutral operator), and this number is given by Burnside’s formula [57]. The roles of symmetries and metrics clearly appear in the CT. The rational component of the compatibility twins does not depend on the metrics. The mirror plane K1 of the type I twins is inherited by the correspondence from the mirror plane of a parent reflection symmetry, and the 180° rotation axis η 2 of the type II twins is inherited by the correspondence from the rotation axis of a parent two-fold rotation symmetry. Only the shear amplitude and the irrational component of the twin (i.e., the direction η 1 for the type I twins or the plane K2 for the type II twin) depend on the metric and, more specifically, on the metric of the martensite phase.
The CT relies on (i) a natural OR, (ii) additional “closing-gap” ORs required for the compatibility between the variants and (iii) a continuum of orientations between the natural OR and the closing-gap ORs. The fact that, for each operator, the twin that is observed experimentally (see Section 5.4) corresponds to one that makes the variants deviate as little as possible from their natural orientations reinforces the theory. The existence of a natural OR has already allowed us to determine the habit planes of the martensite products [24]. The CT includes the idea that the interfaces between the variants are not necessarily fully invariant but can be slightly distorted [73]. This concept of “weak plane” was detailed in Ref. [74]. The “weak plane” hypothesis permits to predict the formation of junctions for variants linked by polar operators ( O 1 and O 3 ), whereas the PTMC is mute on them. The CT predicts that these junction planes are ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 . We also suggest that the junction planes of type II twins could also be rational weak planes. The predicted rational weak planes and the irrational invariant planes are, however, so close in B19′ martensite that it is not possible for the moment to confirm or infirm their existence. Whatever the operator (polar or ambivalent), whatever the twin (type I, type II or weak), the CT shows that the rational elements of the twins are generic, i.e., the mirror plane for the type I twins, the 180° rotation axis for the type II twins and the rotation axis for the weak twin do not depend on the metrics; they are inherited from a parent symmetry element by correspondence and are preserved by the distortion. More generally, the CT gives a geometrical representation of the variants with the rotation gradients and the compatibility twins (Figure 1).

5.4. TKD Observations of the Junction Planes in NiTi Alloys

The predictions made by the CT are now compared with experimental TKD maps. The TEM samples of fully martensitic NiTi samples have been prepared by dual-jet electropolishing. The TKD orientation maps of the B19′ martensite were acquired around the thin edges around the hole of the lamellae. The details of the sample preparation and TKD acquisition were given in Ref. [24]. ARPGE was then used to automatically plot the boundaries between the B19′ variants with seven colors (one per operator), as shown in the map of Figure 5b. The colors chosen for the two complementary polar operators O 1 and O 3 are red and blue, respectively, but the boundaries appear in purple at medium resolution, because the two colors overlap. It can be observed that some junction planes in Figure 5b are colored in purple, which means that the junctions for polar operators do exist, even if the PTMC equations are not solvable for them. Now, the traces of the junction planes predicted by the CT are plotted on the TKD maps with ARPGE. Practically, the “user” chooses the specific operator O i and enters the expected pair of parallel planes, for example, the operator O 2 and the planes ( 100 ) B 19 . ARPGE then considers in the map all the pairs of B19′ laths that are in contact and for which the misorientation is close to that of the operator O 2 , i.e., whose misorientation is a rotation close to ( 180 ° , [ 001 ] B 19 ) with a tolerance of 5°. Then, it plots the traces of ( 100 ) B 19 only if the plane ( 100 ) B 19 of one variant is parallel to the plane ( 100 ) B 19 of the other variant, with a tolerance of 5°. For the polar operators O 1 and O 3 , the user can enter a parallelism condition between nonequivalent planes, for example, ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 .
The analyses of the traces of the junction planes in the TKD map of Figure 5 are presented in Figure 6. For the operator O 2 , the traces of the ( 100 ) B 19 planes predicted for the compound twin are parallel to the boundaries (Figure 6a). For the operator O 4 , the traces of the predicted ( 1 ¯ 11 ) B 19 planes of the type I twin also agree with the boundaries (Figure 6b). For the operator O 5 , the ~ ( 2 ¯ 13 ) B 19 irrational planes of the type II twin or the ( 2 ¯ 13 ) B 19 rational plane of the weak twin also fit well with most of the boundaries (Figure 6c, white rectangles). We noticed, however, that some boundaries are not parallel to ~ ( 2 ¯ 13 ) B 19 or to ( 111 ) B 19 , which is rational plane predicted for the type I twins. We have thus investigated the possibility of other weak twins of axis [ 211 ] B 19 that would be inherited by correspondence from [ 011 ] B 2 , and we found by increasing the tolerances in GenOVa that, besides the ( 2 ¯ 13 ) B 19 , a weak plane ( 1 ¯ 11 ) B 19 could be possible. This plane gives a better agreement (Figure 6c, red rectangle). For the operator O 6 , a good accordance was found with the predicted irrational ~ ( 34 4 ¯ ) B 19 K2 plane or with the rational ( 34 4 ¯ ) B 19 weak planes (Figure 6d). For the polar operators O 1 and O 3 , approximatively half of the grain boundaries are parallel to the predicted ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 weak plane of axis [ 011 ] B 19 (Figure 6e, yellow rectangles). A quarter of the grain boundaries agree with the other weak planes ( 0 1 ¯ 1 ) B 19 ( 100 ) B 19 that are also along the axis [ 011 ] B 19 (Figure 6e, red rectangle). No agreement could be found for the remaining quarter.
The analysis of the traces of the junction planes in another TKD map is given in Appendix C. Globally, a good match between the traces of the predicted junction planes and the B19′ boundaries in the TKD maps supports the CT. To our knowledge, only CT can predict the type of twin, i.e., type I or type II (or weak), and only the concept of weak twin can explain the formation of junction planes for variants linked by polar operators.

6. Discussion

Let us now continue the comparison between PTMC and CT initiated with the example of the NiTi alloys in Section 5. The notion of invariant plane is at the core of the PTMC. Indeed, the PTMC is founded on the idea that the habit plane of the martensite products should be invariant and that the junction plane between two variants in a martensite product should be that of a simple shear. The PTMC calculations use stretch matrices to predict the junction planes between the variants in the pairs; they are quite laborious, because all the possible pairs are explored, and they do not permit to distinguish the relative roles of the symmetries and the metrics. When the equations have no solutions, the PTMC assumes that there are no junction planes, and when the equations can be solved, two solutions are found, one corresponding to a type I twin and the other one to a type II twin. Both have the same shear amplitude, and the PTMC remains mute on which of the two should form. Eventually, some twins appear “generic”, i.e., independent of the metrics, but the reason is blurred by the details of the calculations.
The CT initially sketched in Ref. [24] proposes another approach. The CT is based on the notion of correspondence and uses rigorous algebraic definitions. The variants are simple cosets, and the operators (types of variant pairs) are double cosets. Both are built in the intersection group, which is the subgroup of symmetries “common” to the parent and daughter phases. In general, the intersection group depends on the notion that is considered to define the term “common” (orientation, distortion or correspondence). However, for many transformations, for example, in NiTi alloys, the orientation, distortion and correspondence intersection groups are the same. It is thus possible to establish one-to-one relations between the orientation, distortion and correspondence variants and one-to-one relations between the orientation, distortion and correspondence operators (Table 1), and there is no need to specify the type of variants or operators that are considered. The number of variants is given by Lagrange’s formula and the number of operators by Burnside’s formula. This approach permits to easily identify the different operators and classify them as “ambivalent” or “polar”, as illustrated in Table 2. These notions were initially introduced by Janovec for ferroelectrics in 1972 [55,56]. The ambivalent operators contain two-fold symmetries of austenite, whereas the polar operators do not. By assuming that the symmetry elements of austenite that are preserved by correspondence become the twin elements in the pairs of martensite variants, the CT predicts the twin characteristics without determining the stretch matrices, as in the PTMC. The calculations are fast and only imply simple matrices. The other advantage of the CT is that it is clear from the beginning on which parts of the twins the metric plays a role.
For the ambivalent operators, the rational K1 plane of a type I twin is inherited by correspondence from a parent symmetry reflection by Equation (1), and the rational η 2 direction of a type II twin (and, more generally, of a weak twin) is inherited by correspondence from a parent symmetry rotation axis by Equation (6). For type I twins, the metric influences the shear amplitude by Equation (2) and the irrational direction η 1 by Equation (4). For type II twins, the metric influences the shear amplitude by Equation (7) and the irrational Miller indices of the K2 plane by Equation (8). For weak twins, the metric affects the number of possible weak planes and their rational Miller indices. The junction planes between the variants in the pairs are the K1 plane for the type I twins and can be the irrational K2 planes for the type II twins by Equation (9) or rational axial weak planes. The weak planes and weak twins were introduced in Section 4.4. A weak plane is a rational plane that can be slightly distorted and transformed into another rational plane whose Miller indices may be nonequivalent to those before distortion. An axial weak plane is a weak plane containing a fully invariant rational direction (Figure 3).
For polar operators, there is no fully invariant plane between the variants, but weak planes may exist. As for the type II twins, their axis is inherited by correspondence from the axis of the austenite rotation symmetry; it is thus “generic”, i.e., metric-independent. For example, in NiTi alloys, it was calculated that the junction “weak” planes could be ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 planes (Table 3), and this prediction is in agreement with many grain boundaries between variants linked by polar operators (Section 5.4).
There is a last but important element in the CT; it is the assumption that there exists a “natural” OR between the austenite and martensite. This OR permits to calculate the “natural distortion” and the planes unrotated by this distortion. They are habit plane candidates. This allowed us to predict the (1 1 ¯ 2)B2//(10 1 ¯ )B19′ habit plane of the martensite products in NiTi alloys [24]. It was shown in Section 5.2 that the assumption of a natural OR also permits to predict which of the type I or type II twins will form for each operator (Table 5). This is impossible in the framework of PTMC, because both twins have the same shear amplitude. Eventually, the CT gives a global picture of the intricacy of the variants. The variants are mainly oriented according to the natural OR, and the compatibility between them is allowed by the twins and the rotation gradients (disclinations), as schematized in Figure 1. The CT predictions on NiTi alloys agree well with the experimental TKD maps, as detailed in Section 5.4. Future investigations will be performed on other metallic and ceramic martensitic alloys with the same approach.
Before ending, we would like to recall that the CT is born from doubts about the PTMC that uses (or abuse) simple shears and invariant plane strains. The new CT twin equations written for the type I and type II twins, i.e., Equations (1)–(5) and (6)–(10), respectively, were given for comparison with the PTMC, but it does not mean that martensite is formed by atomic displacements that would follow an invariant plane strain. As explained in Ref. [73], the concept of simple shear is probably not appropriate to describe martensitic transformations. In our point of view, there is no need to consider the sessile or glissile characters of the dislocations in the junction plane to determine whether an interface “moves” or not, because it is the transformation that makes the interface “move” and not the movement of the interface that generates the transformation, as usually believed [76]. The junction plane is “just” a boundary surrounded by a strong elastic field due to the gradients of rotations between the closing-gap ORs and the natural OR, and it “moves” only because of the B19′ → B19′ reorientation under stress and not the opposite. As already introduced in Ref. [24], it is highly plausible that the variant reorientations and the large recoverable plasticity of the martensitic shape memory alloys result from a double transformation: (martensite variant not well orientated in the applied stress field) → (come back to the parent austenite) → (martensite variant well orientated in the applied stress field). In order to check this hypothesis, a straight bar was deformed at room temperature (10 mm long, 2 mm thick) of the same martensitic NiTi alloy as that used in the present paper and in Ref. [24] such that the bar gets a U-shape after deformation. The EBSD map in the deformed regions shows that the martensite is not anymore lenticular as it was in the un-deformed state; it is now nearly equiaxed, as illustrated in Figure 7a. The reconstruction of the parent B2 grains with ARPGE [24] confirms the similarities between the reorientated equiaxed B19′ martensite grains and the prior B2 grains (Figure 7b). The equiaxed morphologies of the B19′ grains formed by reorientation would be difficult to explain with the concept of “moving interfaces”, whereas it is expected by a variant reorientation by a B19′ → B2 → B19′ mechanism. The variant selection by reorientation under stress is also marked by the fact that most of the B19′ grains have the same reddish color in the IPF-X map (X is the tensile direction). The quantitative analysis of such maps will be the subject of a future study.

7. Conclusions

The present paper proposes a theory called “correspondence theory” (CT) to explain the main crystallographic features of the martensitic variants in shape memory alloys. It is an alternative to the PTMC mainly based on the correspondence matrix, on the symmetries of the austenite and martensite phases and on the assumption of a “natural” orientation relationship between the two phases. The CT assumes that that the compatibility between the martensite variants is obtained along austenite symmetry elements that become correspondence twin elements for the variants. The symmetry matrices can be partitioned into double cosets called “operators”. The operators that contain two-fold symmetries (reflections or 180° rotations) are ambivalent, and those that do not are said to be “polar”. The formation of the compatibility twins makes the variants deviate from their “natural” orientations, and the accommodation is obtained by rotation gradients. The variants, the operators and the rational elements of twins (i.e., the twin plane K 1 for the type I twins and the twin direction η 2 for the type II twins) can be calculated without knowing the metrics. Only the irrational elements of the twins (i.e., the twin direction η 1 for the type I twins and the twin planes K 2 for the type II twins) depend on the metric of the martensite phase. For the ambivalent operators, the junction plane is a rational twin plane K 1 of a type I twin inherited from a parent mirror symmetry or an irrational plane K 2 of a type II twin around the axis η 2 or an axial weak plane. For these rotational twins, the rotation axis is inherited from a parent rotational symmetry. A weak plane is a plane that can be slightly distorted and transformed into another rational plane whose Miller indices are not necessarily equivalent to those of the initial plane. By considering the deviation between the closing-gap OR and the natural OR, the CT permits to predict which of the type I or type II twin should form. For the polar operators, the junction plane is an axial weak plane whose axis is inherited from a parent rotational symmetry.
The CT has been applied to the B2 → B19′ martensitic transformation in NiTi alloys. It gives a global understanding of the structure of variants with their main “natural” orientations, habit planes, twins, junction planes, closing-gap ORs and orientation gradients. The predictions were shown to agree well with the traces of the boundaries in the experimental TKD maps.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Roland Logé is acknowledged for his support.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. How to Cook Cosets and Double Cosets

The left cosets g   represent the specific relations between the parent crystal and its variants, such as their orientations (by using the matrix T) or their correspondence (by using the matrix C). A coset decomposition of a group G by a subgroup is quite simple. Consider the set S = G , take as first coset, i.e., c 1 = , and remove all the elements of c 1 that are in S , i.e., replace S S c 1 . Take an element g of the new set S , for example, the first one, and form the new coset c 2 by multiplying g at its right by all the elements of , i.e., c 2 = g   = { g h i ,   h i } , then remove all these elements from S , i.e., and replace S S c 2 . Repeat the operation until there is no element anymore in S . In this paper, the elements of G are 3 × 3 symmetry matrices written in the crystallographic basis, i.e., they contain only 0, 1 and −1. The set of the left cosets g   is noted G / and is called the quotient set. The elements of G are not anymore individual symmetry matrices, as in G or , but sets of matrices. The indices i of the cosets c i are the indices of the variants; they depend on the order of the symmetry elements in G and and are thus arbitrary.
The double cosets   g   represent the specific relations between the variants, such as their misorientations (by using the matrix T) or their inter-correspondence (by using the matrix C). We called them operators. A double coset decomposition of a group G by a subgroup is as simple as for simple cosets. Consider the set S = G , take as first double coset, i.e., d 1 = , and remove all the elements of d 1 that were in S , i.e., replace S S d 1 . Take an element g of the new set S , for example, the first one, and form the new double coset d 2 by multiplying g at its right and left by all the elements of , i.e., d 2 =   g   = { h i g h j ,   h i ,   h j } , then remove all these elements from S , i.e., and replace S S d 2 . Repeat the operation until there is no element anymore in S . The set of double cosets   g   is noted \ G and is called the double quotient set. The elements of \ G are not anymore individual symmetry matrices, as in G or , but sets of matrices. The indices i of the double cosets d i are the indices of the operators; they depend on the order of the symmetry elements in G and and are thus arbitrary. Each operator can also be written as a set of variants (i.e., cosets c i ) written in reference to variant 1 and is thus a set of pairs ( c 1 , c i ) . This allows to establish the composition table of the operators that we called the groupoid composition table [57]. This table is used in ARPGE to reconstruct the prior parent grains from EBSD maps [58,59].

Appendix B. Illustration of the Correspondence Theory with a Simple 2D Example

The CT assumes that the compatibility between two variants i and j is obtained along the prior austenite symmetry element that should be “lost” by the natural distortion but that is preserved by correspondence. Let us explain this idea with the simple 2D example of transformation from a square “austenite” (A, noted here γ) to a parallelogram “martensite” (M, noted here α). The lattice parameter of the square is aγ, and those of the parallelogram are aα, bα and the angle θ . We consider the case where aα = aγ, r = b α   a α ~ 1 and θ ~ 90 ° . We assume that there is a natural OR for which the dense directions are parallel, here [1,0]γ//[1,0]α and (0,1)γ//(0,1)α. The correspondence, orientation and distortion matrices are C = ( 1 0 0 1 ) , T = F = ( 1 r cos θ 0 r   sin θ ) . The γ point group contain eight symmetries that are the identity I ; the inversion I ¯ ; the four mirror symmetries m x γ , m y γ , m x   y ¯   γ and m xy γ and the two four-fold rotations R π / 2 γ   and R π / 2 γ . The parallelogram symmetries are just I and I ¯ . Here, as shown in Table A1, the subgroup of the common symmetries does not depend on the type of relation that is considered (correspondence C, orientation T and distortion F): C γ = T γ = F γ = { I ,   I ¯ } . We note it as γ . The C, T and F variants are all simple cosets based on this subgroup α i = g i γ γ , multiplied at the right by the matrix that encodes the type of relation, C, T or F. According to Lagrange’s formula, there are 8/2 = 4 variants. They are represented in Figure 1a (A = γ and M = α). As for the variants, the operators do not depend on the type of relation, C, T or F. They are expressed by the double cosets O i =   γ g i γ γ multiplied, respectively, at the right and left by the matrix and its inverse that encode the type of relation, C, T or F, as shown in Table A1. Enumerating the double cosets or applying Burnside’s formula shows that there are four operators.
Table A1. Intersection group, variants and operators in the 2D case of a square (γ) to the parallelogram (α) transformation of correspondence, orientation and distortion types of relations.
Table A1. Intersection group, variants and operators in the 2D case of a square (γ) to the parallelogram (α) transformation of correspondence, orientation and distortion types of relations.
CorrespondenceOrientationDistortion
Intersection subgroups
C γ = G γ   C   G α C 1   T γ = G γ   T   G α T 1   F γ = G γ   F   G α F 1  
= { I ,   I ¯ } = { I ,   I ¯ } = { I ,   I ¯ }
Variants (from simple cosets)
g i γ γ C g i γ   γ T g i γ   γ F
Operators (from double cosets)
C 1 C γ   g j γ C γ   C T 1 T γ   g j γ T γ T   F 1 F γ   g j γ F γ   F
For each operator O i , the disorientation is the rotation with the smallest rotation angles among all the equivalent rotations in the set of matrices T 1 T γ   g j γ T γ T . They are given in the second column of Table A2. The double coset O i is written as a list of γ symmetry matrices or geometrical elements or as a list of arrows between the α variants, given in columns 3, 4 and 5 of Table A2. The operator O 0 is the identity for all the variants, as it leaves each of them invariant. The operators O 1 and O 2 are ambivalent, since their double cosets contain two-fold matrices; geometrically, they contain arrows, and their inverses, for example, the operator O 1 contain the arrows ( α 1 α 3 ) and ( α 3 α 1 ) . A polar operator is a double coset containing matrices whose inverses belong to another operator; the two operators are then called “complementary”. The operator O 2 is expected to be polar, because it contains only the four-fold symmetry elements; however, in this example, since ( R π / 2 γ ) 2 = I ¯ and I ¯ belongs to γ , the operator is ambivalent.
Table A2. Operators O i between α martensite parallelograms. The disorientations are the second column. The other columns give the double cosets of γ austenite symmetry matrices, with their corresponding geometrical elements, and their set of arrows from one variant (the source) to another variant (the target). In this example, the variants do not depend on the relation C, T or F:   O i = C γ   g i γ C γ   =   T γ   g i γ T γ   =   F γ   g i γ F γ .
Table A2. Operators O i between α martensite parallelograms. The disorientations are the second column. The other columns give the double cosets of γ austenite symmetry matrices, with their corresponding geometrical elements, and their set of arrows from one variant (the source) to another variant (the target). In this example, the variants do not depend on the relation C, T or F:   O i = C γ   g i γ C γ   =   T γ   g i γ T γ   =   F γ   g i γ F γ .
Disorient .   between   α   Variants γ   Symmetry   Matrices   in   the   Double   Coset   T γ   g i γ T γ With Their
Geometrical
Elements
Operators as Set of Equivalent Arrows
(Source Variant → Target Variant)
O 0 I ( 1 0 0 1 ) ( 1 0 0 1 ) I , I ¯ ( α 1 α 1 ) ,   ( α 2 α 2 ) ,   ( α 3 α 3 ) ,   ( α 4 α 4 )
O 1 R π , ~ [ 0 , 1 ] α ( 1 0 0 1 ) ( 1 0 0 1 ) m x γ ,   m y γ ( α 1 α 3 ) ,   ( α 2 α 4 ) ,   ( α 3 α 1 ) ,   ( α 4 α 2 )
O 2 R π , [ 1 , 1 ] α ( 0 1 1 0 ) ( 0 1 1 0 ) m xy γ ,   m x   y ¯   γ ( α 1 α 2 ) ,   ( α 2 α 1 ) ,   ( α 3 α 4 ) ,   ( α 4 α 3 )
O 3 R π / 2 α ( 0 1 1 0 ) ( 0 1 1 0 ) R π / 2 γ , R π / 2 γ ( α 1 α 4 ) ,   ( α 2 α 3 ) ,   ( α 3 α 2 ) ,   ( α 4 α 1 )
The junction planes can now be predicted by considering the correspondence operators following the idea initiated in our last work [24]. First, we consider the operator O 1 or, equivalently, the pairs ( α 1 , α 3 ) and ( α 2 , α 4 ) in Table A2. The two variants in each pair have two planes that are in correspondence; they are the planes (1,0)α and (0,1)α inherited from the planes m x γ = ( 0 , 1 ) γ and m y γ = ( 1 , 0 ) γ , as shown in Figure 1. The horizontal plane m x γ is a natural junction plane for the variants α 1 and α 3 , because it is invariant by their distortions. More interestingly, the plane m y γ is also in inter-correspondence between variants α 1 and α 3 , but it is rotated clockwise by the distortion of α1 and anticlockwise by the distortion of α3, which means that the planes ( 1 , 0 ) α 1 and ( 1 , 0 ) α 3 are not parallel anymore once the natural distortions are completed. The compatibility between the variants α 1 and α 3 can, however, be maintained along the plane m y γ , thanks to a slight rotation counterclockwise for α 1 and clockwise for α 3 . We make the hypothesis that the two variants remain joined by a continuum of orientation, such that the plane ( 1 , 0 ) α 1 is continuously rotated counterclockwise and ( 1 , 0 ) α 3 is continuously rotated clockwise to come back parallel to m y γ . This induces a continuum of ORs from the natural OR toward a new OR that is ( 1 , 0 ) α 1 r ( 1 , 0 ) α 3 r ( 1 , 0 ) γ . This OR appears as a secondary OR generated from the natural OR to permit compatibility between the variants α1 and α3 along the plane m y γ , as schematized in Figure 1b. We have noted α 1 r and α 3 r in place of α 1 and α 3 because of the slight rotation. The junction plane between the variants α 1 and α 3 is the horizontal plane m x γ , whereas the junction plane between the variants α 1 r and α 3 r is the vertical plane m y γ . There are thus two possible junction planes for the pair of variants ( α 1 , α 3 ). The same reasoning equally applies to the pair ( α 2 , α 4 ).
The operators O 2 are made of the pairs ( α 1 , α 2 ) and ( α 3 , α 4 ). The distortions associated with α 1 and α 2 leave none of the directions invariant. However, here again, one can notice that the variants are linked by the correspondence m x y γ ( 1 ¯ , 1 ) α 1 and m x y γ ( 1 ¯ , 1 ) α 2 , which induces a secondary OR with rotated variants ( 1 ¯ , 1 ) α 1 r ( 1 ¯ , 1 ) α 2 r ( 1 ¯ , 1 ) γ , as shown in Figure 1c. They are also linked by the correspondence m x   y ¯   γ ( 1 , 1 ) α 1 and m x   y ¯   γ ( 1 , 1 ) α 2 , which induces a secondary OR with rotated variants ( 1 , 1 ) α 1 r ( 1 , 1 ) α 2 r ( 1 , 1 ) γ .
The operators O 3 is made of the pairs ( α 1 , α 4 ) and ( α 2 , α 3 ) . It does not contain mirror symmetry or two-fold rotation. There is thus no obvious junction plane. Depending on the metric of the α phase, some weak planes could be foreseen, as is the case for NiTi alloys (see Section 4.4), but here, the calculations show no possible weak twin with a low generalized strain value, which means that there is no expected junction plane between the variants α 1 and α 4 and between α 2 and α 3 .
In the CT, since the “natural” OR remains the “absolute” OR on which the structure keeps its integrity, all the distortion variants can be formed according to the same principle, without requiring specific additional rotation Qij, Jij, Jijkl, etc., as in the PTMC. The global structure made of all the variants is preserved by the “closing-gap” ORs and their associated continuum of ORs. The gradients orientations predicted by the CT seem to be similar to those already observed in the pole figures of EBSD maps of martensite in steels [60,61,62] and deformation twins in magnesium [75].
Similar results could be obtained with the PTMC, but in a more complex way. One would indeed need to calculate the matrices Δ F ( i , j ) = F j T F i T F i   F j 1 = U j T U i T U i U j 1 for all the pairs of distortion variants F i and F j or of stretch variants U i and U j . The notation Δ F ( i , j ) was chosen in place of the usual notation often noted C ( i , j ) to avoid confusion with the inter-correspondence matrices presented in the main text. Let us do the calculations with r = 1.2 and θ = 80 ° .
For the pairs of variants ( α 1 ,   α 3 ) and ( α 2 ,   α 4 ) , the matrix Δ F = ( 1 . 0.35266 0.35266 1.12436 ) . The eigenvalues are λ 1 = 0.70409 and λ 3 = 1.42028 . We can assume that the value λ 2 = 1 is along a third dimension normal to the figure. The related eigenvectors are e 1 = [ 0.76604 , 0.64279 ] and e 3 = [ 0.64279 ,   0.76604 ] . Applying Ball and James’ formula [46] also reported in Bhattacharya’s book ([2] p. 70), it is eventually numerically found that the junction planes are n ^ = ( 0 , 1 ) or n ^ = ( 1 , 0 ) .
For the pairs of variants ( α 1 ,   α 2 ) and ( α 3 ,   α 4 ) , the matrix Δ F = ( 1.4400 0.07758 0.07758 0.69862 ) . The eigenvalues are λ 1 = 0.69059 and λ 3 = 1.44803 . The related eigenvectors are e 1 = [ 0.45677 ,   0.88959 ] and e 3 = [ 0.99468 , 0.10298 ] . Applying Ball and James’ formula, it is eventually numerically found that the junction planes are n ^ = ( 1 , 1 ) or n ^ = ( 1 , 1 ) .
For the pairs of variants ( α 1 ,   α 2 ) and ( α 2 ,   α 3 ) , the matrix Δ F = ( 1.43400 0.43024 0.43024 0.82299 ) . The eigenvalues are λ 1 = 0.60208 and λ 3 = 1.66091 . The related eigenvectors are e 1 = [ 0.45677 ,   0.88959 ] and e 3 = [ 0.88959 ,   0.45677 ] . Applying Ball and James’ formula, it is eventually numerically found that the junction planes are n ^ = ( 0.21683 ,   0.50402 ) or n ^ = ( 0.50402 ,   0.21683 ) .
Consequently, PTMC and the correspondence theory give the same junction planes for the operators O 1 and O 2 . It should be noticed, however, that the PTMC requires intermediate calculations that mask the fact that the junction planes are inherited from parent mirror planes. The other difference is for the operator O 3 . The PTMC predicts non-integer junction planes whose values depend on the metrics of the parent and daughter phases, whereas the correspondence theory does not predict any junction plane. We will see an opposite situation in the case of B2–B19′ transformation in NiTi, where the PTMC predicts no junction plane for the polar operators, whereas the CT predicts weak planes that will be confirmed by TKD.

Appendix C. Additional TKD Map and Junction Planes

We have acquired various TKD maps and analyzed the traces of the junction planes with ARPGE, as explained in Section 5. Here, in Figure A1, we show another TKD map where the grain boundaries between the B19′ variants are colored according to the operators rather than link them. The comparison of the traces with the predictions of the correspondence theory on different zones of the maps is shown in Figure A2. The agreement is satisfying, except for some boundaries between variants linked by the polar operators.
Figure A1. TKD map of a fully martensitic B19′ NiTi alloy with the junction boundaries colored as function of the operators, such as in Figure 5.
Figure A1. TKD map of a fully martensitic B19′ NiTi alloy with the junction boundaries colored as function of the operators, such as in Figure 5.
Crystals 12 00130 g0a1
Figure A2. TKD map of a fully martensitic B19′ NiTi alloy with the junction boundaries colored as function of the operators, such as in Figure 5.
Figure A2. TKD map of a fully martensitic B19′ NiTi alloy with the junction boundaries colored as function of the operators, such as in Figure 5.
Crystals 12 00130 g0a2

References

  1. Buehler, W.J.; Gilfrich, J.W.; Wiley, R.C. Effect of Low-Temperature Phase Changes on the Mechanical Properties of Alloys near Composition TiNi. J. Appl. Phys. 1963, 34, 1473–1477. [Google Scholar] [CrossRef]
  2. Bhattacharya, K. Microstructure of Martensite. In Why It Forms and How It Gives Rise to the Shape-Memory Effect, 1st ed.; Oxford University Press: New York, NY, USA, 2003. [Google Scholar]
  3. Otsuka, K.; Ren, X. Physical metallurgy of Ti-Ni-based shape memory alloys. Prog. Mat. Sci. 2005, 50, 511–678. [Google Scholar] [CrossRef]
  4. Chowdhury, P.; Sehitoglu, H. Deformation physics of shape memory alloys—Fundamentals at atomistic frontier. Prog. Mater. Sci. 2017, 88, 49–88. [Google Scholar] [CrossRef]
  5. Otsuka, K.; Sawamura, T.; Shimizu, K. Crystal structure and internal defects of equiatomic TiNi martensite. Phys. Stat. Sol. 1971, 5, 457–470. [Google Scholar] [CrossRef]
  6. Knowles, K.M.; Smith, D.A. The crystallography of the martensitic transformation in equiatomic nickel-titanium. Acta Metall. 1981, 29, 101–110. [Google Scholar] [CrossRef]
  7. Nishida, M.; Nishiura, T.; Kawano, H.; Inamura, T. Self-accommodation of B19’ martensite in Ti-Ni shape memory alloys-Part I. Morphological and crystallographic studies of the variant selection rules. Philos. Mag. 2012, 92, 2215–2223. [Google Scholar] [CrossRef]
  8. Nishida, M.; Okunishi, E.; Nishiura, T.; Kawano, H.; Inamura, T.; Li, S.; Hara, T. Self-accommodation of B19’ martensite in Ti-Ni shape memory alloys-Part II. Characteristic interface structures between habit plane variants. Philos. Mag. 2012, 92, 2234–2246. [Google Scholar] [CrossRef]
  9. Šittner, P.; Molnárová, O.; Kadeřávek, L.; Tyc, O.; Heller, L. Deformation twinning in martensite affecting functional behavior of NiTi shape memory alloys. Materialia 2020, 9, 100506. [Google Scholar] [CrossRef]
  10. Waitz, T. The self-accommodated morphology of martensite in nanocrystalline NiTi shape memory alloys. Acta Mater. 2005, 53, 2273–2283. [Google Scholar] [CrossRef]
  11. Liu, Y.; Xie, Z.L. Twinning and detwinning of <011> type II twin in shape memory alloy. Acta Mater. 2003, 51, 5529–5543. [Google Scholar] [CrossRef]
  12. Liu, Y.; Xie, Z.L. The rational Nature of Type II Twin in NiTi Shape Memory Alloy. J. Intell. Mater. Struct. 2006, 17, 1083–1090. [Google Scholar] [CrossRef]
  13. Ye, J.; Mishra, R.K.; Pelton, A.R.; Minor, A.M. Direct observation of the NiTi martensitic phase transformation in nanoscale volumes. Acta Mater. 2010, 58, 490–498. [Google Scholar] [CrossRef]
  14. Ii, S.; Yamauchi, K.; Maruhashi, Y.; Nishida, M. Direct evidence of correlation between {20-1}B19′ and {114}B2 deformation twins in Ti–Ni shape memory alloy. Scr. Mater. 2003, 49, 723–727. [Google Scholar] [CrossRef]
  15. Šittner, P.; Sedlák, P.; Seiner, H.; Sedmák, P.; Pilch, J.; Delville, R.; Heller, L.; Kadeřávek, L. On the coupling between martensitic transformation and plasticity in NiTi: Experiments and continuum based modelling. Prog. Mater. Sci. 2018, 98, 249–298. [Google Scholar] [CrossRef]
  16. Laplanche, G.; Birk, T.; Schneider, S.; Frenzel, J.; Eggeler, G. Effect of temperature and texture on the reorientation of martensite variants in NiTi shape memory alloys. Acta Mater. 2017, 127, 143–152. [Google Scholar] [CrossRef]
  17. Xiao, Y.; Zeng, P.; Lei, L.; Zhang, Y. In situ observation on temperature dependence of martensitic transformation and plastic deformation in superelastic NiTi shape memory alloy. Mater. Des. 2017, 134, 111–120. [Google Scholar] [CrossRef]
  18. Hsu, W.-N.; Polatidis, E.; Šmíd, M.; Casati, N.; Van Petegem, S.; Van Swygenhoven, H. Load path change on superelastic NiTi alloys: In situ synchrotron XRD and SEM DIC. Acta Mater. 2018, 144, 874–883. [Google Scholar] [CrossRef]
  19. Hsu, W.-N.; Polatidis, E.; Šmíd, M.; Van Petegem, S.; Casati, N.; Van Swygenhoven, H. Deformation and degradation of superelastic NiTi under multiaxial loading. Acta Mater. 2019, 167, 149–158. [Google Scholar] [CrossRef]
  20. Polatidis, E.; Šmíd, M.; Kuběna, I.; Hsu, W.-N.; Laplanche, G.; Van Swygenhoven, H. Deformation mechanisms in a superelastic NiTi alloy: An in-situ high resolution digital image correlation study. Mater. Des. 2020, 191, 108622. [Google Scholar] [CrossRef]
  21. Keller, R.R.; Geiss, R.H. Transmission EBSD from 10 nm domains in a scanning electron microscope. J. Microsc. 2012, 245, 245–251. [Google Scholar] [CrossRef]
  22. Trimby, T. Orientation mapping of nanostructured materials using transmission Kikuchi diffraction in the scanning electron microscope. Ultramicrosc 2012, 120, 16–24. [Google Scholar] [CrossRef] [PubMed]
  23. Suzuki, S. Feature of Transmission EBSD and its application. JOM 2013, 65, 1254–1263. [Google Scholar] [CrossRef] [Green Version]
  24. Cayron, C. What EBSD and TKD tell us about the crystallography of the martensitic B2-B19′ transformation in NiTi shape memory alloys. Crystals 2020, 10, 562. [Google Scholar] [CrossRef]
  25. Bowles, J.S.; Mackenzie, J.K. The crystallography of martensitic transformations I. Acta Metall. 1954, 2, 129–137. [Google Scholar] [CrossRef]
  26. Bowles, J.S.; Mackenzie, J.K. The crystallography of martensitic transformations II. Acta Metall. 1954, 2, 138–147. [Google Scholar] [CrossRef]
  27. Wechsler, M.S.; Liebermann, D.S.; Read, T.A. On the theory of the formation of martensite. Trans. AIME 1953, 197, 1503–1515. [Google Scholar]
  28. Christian, J. Applications of the phenomenological theories of martensite I–II. J. Inst. Met. 1956, 84, 386–394. [Google Scholar]
  29. Bain, E.C.; Dunkirk, N. The nature of martensite. Trans. AIME 1924, 70, 25–47. [Google Scholar]
  30. Bhadeshia, H.K.D.H. Worked Examples in the Geometry of Crystals, 2nd ed.; The Institute of Metals: London, UK, 1987. [Google Scholar]
  31. Wayman, C.M. Shear transformations and microstructures. Metallography 1975, 8, 105–130, republished in Mater. Charact. 1997, 39, 235–260. [Google Scholar] [CrossRef]
  32. Duparc, O.B.M.H. A review of some elements for the history of mechanical twinning centred on its German origins until Otto Mügge’s K1 and K2 invariant plane notation. J. Mater. Sci. 2016, 52, 4182–4196. [Google Scholar] [CrossRef]
  33. Mügge, O. Ueber Homogene Deformationen (einfache Schiebungen) an den Triklinen Doppelsalzen BaCdCl4.4aq., Neues Jahrbuch für Mineralogie, Geologie und Palaeontologie Beilage-Band. 1889. Available online: https://www.phase-trans.msm.cam.ac.uk/2001/geometry2/Geometry.pdf (accessed on 22 November 2021).
  34. Kihô, H. The Crystallographic Aspect of the Mechanical Twinning in Metals. J. Phys. Soc. Jpn. 1954, 9, 739–747. [Google Scholar] [CrossRef]
  35. Jaswon, M.A.; Dove, D.B. The Crystallography of Deformation Twinning. Acta Cryst. 1960, 13, 232–240. [Google Scholar] [CrossRef] [Green Version]
  36. Bilby, B.A.; Crocker, A.G. The Theory of the Crystallography of Deformation Twinning. Proc. R. Soc. Lond. Ser. A 1965, 288, 240–255. [Google Scholar]
  37. Bevis, M.; Crocker, A.G. Twinning Shears in Lattices. Proc. R. Soc. Lond. Ser. A 1968, 304, 123–134. [Google Scholar]
  38. Cahn, R.W. Twinned Crystals. Adv. Phys. 1954, 3, 363–445. [Google Scholar] [CrossRef]
  39. Christian, J.W.; Mahajan, S. Deformation Twinning. Prog. Mater. Sci. 1995, 39, 1–157. [Google Scholar] [CrossRef]
  40. Gupta, S.P.; Johnson, A.A. Morphology and crystallography of beta’ martensite in TiNi alloys. Trans. Jpn. Inst. Metals. 1973, 14, 292–302. [Google Scholar] [CrossRef] [Green Version]
  41. Sinclair, R. Origin of stacking faults in NiTi martensite. AIP Conf. 1979, 53, 269. [Google Scholar]
  42. Nishida, M.; Yamauchi, K.; Itai, I.; Ohgi, H.; Chiba, A. High resolution electron microscopy studies of twin boundary structure in B19’ martensite in the Ti-Ni shape memory alloy. Acta Metall. Mater. 1995, 43, 1229–1234. [Google Scholar] [CrossRef]
  43. Mohammed, A.S.K.; Sehitoglu, H. Modeling the interface structure of type II boundary in B19’ NiTi from an atomistic and topological standpoint. Acta Mater. 2020, 183, 93–109. [Google Scholar] [CrossRef]
  44. Miyazaki, S.; Kimura, S.; Otsuka, K.; Suzuki, Y. The habit plane and transformation strains associated with martensitic transformation in Ti-Ni single crystals. Scr. Metall. 1984, 18, 883–888. [Google Scholar] [CrossRef]
  45. Miyazaki, S.; Otsuka, K.; Wayman, C.M. The shape memory mechanism associated with the martensitic transformation in Ti-Ni alloys-I. Self-accommodation. Acta Metall. 1989, 37, 1873–1884. [Google Scholar] [CrossRef]
  46. Ball, J.M.; James, R.D. Finite phase mixtures as minimizers of energy. Arch. Ration. Mech. Anal. 1987, 100, 13–52. [Google Scholar] [CrossRef]
  47. Pitteri, M.; Zanzotto, G. Generic and non-generic cubic-to-monoclinic transitions and their twins. Acta Mater. 1998, 46, 225–237. [Google Scholar] [CrossRef]
  48. Pitteri, M.; Zanzotto, G. Continuum Models for Phase Transitions and Twinning in Crystals; Chapman & Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
  49. Hane, K.F.; Shield, T.W. Microstructure in the cubic to monoclinic transition in titanium-nickel shape memory alloys. Acta Mater. 1999, 47, 2603–2617. [Google Scholar] [CrossRef]
  50. Inamura, T.; Nishiura, T.; Kawano, H.; Hosoda, H.; Nishida, M. Self-accommodation of B19’ martensite in Ti-Ni shape memory alloys-Part III. Analysis of habit plane variant clusters by the geometrically nonlinear theory. Philos. Mag. 2012, 92, 2247–2263. [Google Scholar] [CrossRef]
  51. Gu, H.; Bumke, L.; Chluba, C.; Quandt, E.; James, R.D. Phase engineering and supercompatibility of shape memory alloys. Mater. Todays 2018, 21, 265–277. [Google Scholar] [CrossRef]
  52. Song, Y.; Chen, X.; Dabade, V.; Shield, T.W.; James, R.D. Enhanced reversibility and unusual microstructure of a phase-transforming material. Nature 2013, 502, 85–88. [Google Scholar] [CrossRef]
  53. Chen, X.; Liu, T.; Li, R.; Liu, J.; Zhao, Y. Molecular dynamics simulation on the shape memory effect and superelasticity in NiTi shape memory alloy. Comput. Mater. Sci. 2018, 146, 61–69. [Google Scholar] [CrossRef]
  54. Li, B.; Shen, Y.; An, Q. Structural origin of reversible martensitic transformation and reversible twinning in NiTi shape memory alloy. Acta Mater. 2020, 199, 240–252. [Google Scholar] [CrossRef]
  55. Janovec, V. Group analysis of domains and domain pairs. Czech. J. Phys. B 1972, 22, 975–994. [Google Scholar] [CrossRef]
  56. Janovec, V.; Hahn, T.; Klapper, H. International Tables for Crystallography; Authier, A., Ed.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003; Volume D, pp. 377–391. [Google Scholar]
  57. Cayron, C. Groupoid of orientational variants. Acta Cryst. 2006, 62, 21–40. [Google Scholar] [CrossRef] [Green Version]
  58. Cayron, C. GenOVa: A computer program to generate orientational variants. J. Appl. Cryst. 2007, 40, 1179–1182. [Google Scholar] [CrossRef] [Green Version]
  59. Cayron, C. ARPGE: A computer program to automatically reconstruct the parent grains from electron backscatter diffraction data. J. Appl. Cryst. 2007, 40, 1183–1188. [Google Scholar] [CrossRef] [Green Version]
  60. Cayron, C.; Barcelo, F.; de Carlan, Y. The mechanisms of the fcc–bcc martensitic transformation revealed by pole figures. Acta Mater. 2010, 58, 1395–1402. [Google Scholar] [CrossRef]
  61. Cayron, C. One-step model of the face-centred-cubic to body-centred-cubic martensitic transformation. Acta Cryst. 2013, 69, 498–509. [Google Scholar] [CrossRef]
  62. Cayron, C. Continuous atomic displacements and lattice distortion during fcc–bcc martensitic transformation. Acta Mater. 2015, 96, 189–202. [Google Scholar] [CrossRef] [Green Version]
  63. Cayron, C. Angular distortive matrices of phase transitions in the fcc-bcc-hcp system. Acta Mater. 2016, 111, 417–441. [Google Scholar] [CrossRef] [Green Version]
  64. Baur, A.P.; Cayron, C.; Logé, R. {225}γ habit planes in martensitic steels: From the PTMC to a continuous model. Sci. Rep. 2017, 7, 40938. [Google Scholar] [CrossRef] [PubMed]
  65. Baur, A.P.; Cayron, C.; Logé, R. Variant selection in surface martensite. J. Appl. Cryst. 2017, 50, 1646–1652. [Google Scholar] [CrossRef] [Green Version]
  66. Baur, A.P.; Cayron, C.; Logé, R. Variant selection in Fe-2Ni-1.8C under bending. Crystals 2018, 8, 474. [Google Scholar] [CrossRef] [Green Version]
  67. Larcher, M.N.D.; Cayron, C.; Blatter, A.; Soulignac, R.; Logé, R. EBSD Study of Variant Selection during Ordering Phase Transformation in L10-type Red Gold Alloy. J. Appl. Cryst. 2019, 52, 1202–1213. [Google Scholar] [CrossRef]
  68. Larcher, M.N.D.; Cayron, C.; Blatter, A.; Soulignac, R.; Logé, R. The thermally activated distortion with amplification effect and related variant selection in red gold alloys. Acta Mater. 2020, 198, 242–256. [Google Scholar] [CrossRef]
  69. Cayron, C. The transformation matrices (Distortion, orientation, correspondence), their continuous forms and their variants. Acta Cryst. 2019, 75, 411–437. [Google Scholar]
  70. Jaswon, M.A.; Wheeler, J.A. Atomic Displacements in the Austenite-Martensite Transformation. Acta Cryst. 1948, 1, 216–224. [Google Scholar] [CrossRef]
  71. Cayron, C. Complements to Mügge and Friedel’s theory of twinning. Metals 2020, 10, 231. [Google Scholar] [CrossRef] [Green Version]
  72. Cayron, C.; Logé, R. Evidence of new twinning modes in magnesium questioning the shear paradigm. J. Appl. Cryst. 2018, 51, 809–817. [Google Scholar] [CrossRef] [Green Version]
  73. Cayron, C. Shifting the shear paradigm in the crystallographic models of displacive transformations in metals and alloys. Crystals 2018, 8, 181. [Google Scholar] [CrossRef] [Green Version]
  74. Cayron, C. The weak plane and weak twins. in preparation.
  75. Cayron, C. Hard-sphere displacive model of extension twinning in magnesium. Mater. Des. 2017, 119, 361–375. [Google Scholar] [CrossRef] [Green Version]
  76. Christian, J.W. Deformation by Moving Interfaces. Metall. Trans. A 1982, 13, 509–538. [Google Scholar] [CrossRef]
Figure 1. A 2D example of a square austenite (A) → parallelogram martensite (M) transformation. (a) There are four orientation variants, assuming a M / / a A for the “natural” OR. (b) Closing-gap OR ( 1 , 0 ) M 1 / / ( 1 , 0 ) M 3 / / ( 1 , 0 ) A and associated rotation gradients between the variants M 1 and M 3 . This new OR maintains invariant the vertical reflection plane m y A . (c) Closing-gap OR ( 1 ¯ , 1 ) M 1 / / ( 1 ¯ , 1 ) M 2 / / ( 1 ¯ , 1 ) A and associated rotation gradients between the variants M 1 and M 2 . This new OR maintains invariant the diagonal reflection plane m x y A . (d) Global compatibility between the variants mediated by the closing-gap ORs and associated orientation gradients.
Figure 1. A 2D example of a square austenite (A) → parallelogram martensite (M) transformation. (a) There are four orientation variants, assuming a M / / a A for the “natural” OR. (b) Closing-gap OR ( 1 , 0 ) M 1 / / ( 1 , 0 ) M 3 / / ( 1 , 0 ) A and associated rotation gradients between the variants M 1 and M 3 . This new OR maintains invariant the vertical reflection plane m y A . (c) Closing-gap OR ( 1 ¯ , 1 ) M 1 / / ( 1 ¯ , 1 ) M 2 / / ( 1 ¯ , 1 ) A and associated rotation gradients between the variants M 1 and M 2 . This new OR maintains invariant the diagonal reflection plane m x y A . (d) Global compatibility between the variants mediated by the closing-gap ORs and associated orientation gradients.
Crystals 12 00130 g001
Figure 2. Junction plane between variants M 1 and M 2 deduced from the austenite symmetry operation that links them. (a) Reflection on the austenite plane p A , and (b) 180° rotation around the austenite axis u A .
Figure 2. Junction plane between variants M 1 and M 2 deduced from the austenite symmetry operation that links them. (a) Reflection on the austenite plane p A , and (b) 180° rotation around the austenite axis u A .
Crystals 12 00130 g002
Figure 3. Axial weak plane between two crystals: 1 and 2. (a) The rational directions u 1 and u 2 are the same vector; it is the invariant axis of the twin. The rational directions v 1 and v 2 have close norms and make nearly the same angle with the axis of the twin. The weak plane is ( u 1 , v 1 ) parallel to ( u 2 , v 2 ) ; it is the plane of the page; its Miller indices are the coordinates of the reciprocal vector p 1 = ( u 1 × v 1 ) in the crystallographic basis of crystal 1 c 1 , and of p 2 = ( u 2 × v 2 ) in c 2 . (b) The slight intraplanar distortion represented by the black arrows transforms p 1 into p 2 .
Figure 3. Axial weak plane between two crystals: 1 and 2. (a) The rational directions u 1 and u 2 are the same vector; it is the invariant axis of the twin. The rational directions v 1 and v 2 have close norms and make nearly the same angle with the axis of the twin. The weak plane is ( u 1 , v 1 ) parallel to ( u 2 , v 2 ) ; it is the plane of the page; its Miller indices are the coordinates of the reciprocal vector p 1 = ( u 1 × v 1 ) in the crystallographic basis of crystal 1 c 1 , and of p 2 = ( u 2 × v 2 ) in c 2 . (b) The slight intraplanar distortion represented by the black arrows transforms p 1 into p 2 .
Crystals 12 00130 g003
Figure 4. Theoretical disorientation histogram and rotation axes built from all the orientation variants of all the closing-gap ORs presented in Table 4.
Figure 4. Theoretical disorientation histogram and rotation axes built from all the orientation variants of all the closing-gap ORs presented in Table 4.
Crystals 12 00130 g004
Figure 5. TKD map of a B19′ martensitic NiTi alloy. (a) IPF-Y colored map of the B19′ grains. (b) Specific grain boundaries colored in a TKD map according to the operators. The polar operators O 1 and O 3 in red and blue, respectively, overlap on each side of the grain boundaries and merge into a purple color.
Figure 5. TKD map of a B19′ martensitic NiTi alloy. (a) IPF-Y colored map of the B19′ grains. (b) Specific grain boundaries colored in a TKD map according to the operators. The polar operators O 1 and O 3 in red and blue, respectively, overlap on each side of the grain boundaries and merge into a purple color.
Crystals 12 00130 g005
Figure 6. Traces of the junction planes predicted by the CT for the operators (a) O 2 , (b) O 4 , (c) O 5 , (d) O 6 (e) O 1 and O 3 . The comparison for each pair of B19′ variants is visually made by considering whether the traces are parallel to the grain boundaries. Good agreement was found for most of the boundaries that were considered.
Figure 6. Traces of the junction planes predicted by the CT for the operators (a) O 2 , (b) O 4 , (c) O 5 , (d) O 6 (e) O 1 and O 3 . The comparison for each pair of B19′ variants is visually made by considering whether the traces are parallel to the grain boundaries. Good agreement was found for most of the boundaries that were considered.
Crystals 12 00130 g006
Figure 7. B19′ martensite reorientated by bending a martensitic NiTi bar (10 mm long and 2 mm thick) at room temperature into a U-shape. (a) IPF-X EBSD map of the B19′ phase. The map was acquired in the tensile part, with the tensile direction along X (horizontal direction). The initial microstructure before bending was constituted of long lenticular laths (see Figure 4b of Ref. [24]). The reddish color of the IPF-X map and the high density close to 3 ¯ 16 B 19 in the IPF map presented on the right-hand side (Multiple of Uniform Distribution, MUD = 27) show that the B19′ variants are reorientated under the stress. (b) IPF-X map of the parent B2 grains reconstructed with ARPGE.
Figure 7. B19′ martensite reorientated by bending a martensitic NiTi bar (10 mm long and 2 mm thick) at room temperature into a U-shape. (a) IPF-X EBSD map of the B19′ phase. The map was acquired in the tensile part, with the tensile direction along X (horizontal direction). The initial microstructure before bending was constituted of long lenticular laths (see Figure 4b of Ref. [24]). The reddish color of the IPF-X map and the high density close to 3 ¯ 16 B 19 in the IPF map presented on the right-hand side (Multiple of Uniform Distribution, MUD = 27) show that the B19′ variants are reorientated under the stress. (b) IPF-X map of the parent B2 grains reconstructed with ARPGE.
Crystals 12 00130 g007
Table 1. Intersection subgroup, variants and operators of the B2 → B19′ transformation for correspondence, orientation and distortion relations. For this transformation, the intersection subgroup does not depend on the type of relation, C B 2 = T B 2 = F B 2 and is just noted as B 2 .
Table 1. Intersection subgroup, variants and operators of the B2 → B19′ transformation for correspondence, orientation and distortion relations. For this transformation, the intersection subgroup does not depend on the type of relation, C B 2 = T B 2 = F B 2 and is just noted as B 2 .
CorrespondenceOrientationDistortion
Intersection subgroups
C B 2 = G B 2   C   G B 19 C 1   T B 2 = G B 2   T   G B 19 T 1   F B 2 = G B 2   F   G B 2 F 1  
= { I ,   I ¯ ,   R π , [ 110 ] B 2 ,   m ( 110 ) B 2 } = { I ,   I ¯ ,   R π , [ 110 ] B 2 ,   m ( 110 ) B 2 } = { I ,   I ¯ ,   R π , [ 110 ] B 2 ,   m ( 110 ) B 2 }
Variants (from simple cosets)
g i B 2 B 2 C g i B 2   B 2   T g i B 2   B 2   F
Number of variants (Lagrange’s formula)
121212
Operators (from double cosets)
C 1 B 2   g j B 2 B 2   C T 1 B 2   g j B 2   B 2 T   F 1 B 2   g j B 2 B 2   F
Number of operators (Burnside’s formula)
777
Table 2. List of operators for the B2 → B19′ transformation. The operators are numbered from 0 (identity operator) to 6; the order results from the double-coset decomposition made in GenOVa. The operators are classified here into three categories: the neutral operator O 0 ; the ambivalent operators { O 2 ,   O 4 ,   O 5   and   O 6 } and the complementary polar operators { O 1   and   O 3 } . The second column gives the disorientation, i.e., the smallest rotation in the set T 1 T B 2   g j B 2 T B 2   T . The third and fourth columns give the B2 symmetries g j B 2 in the double cosets by their matrices and by symmetry elements, respectively. The elements in bold are two-fold symmetries.
Table 2. List of operators for the B2 → B19′ transformation. The operators are numbered from 0 (identity operator) to 6; the order results from the double-coset decomposition made in GenOVa. The operators are classified here into three categories: the neutral operator O 0 ; the ambivalent operators { O 2 ,   O 4 ,   O 5   and   O 6 } and the complementary polar operators { O 1   and   O 3 } . The second column gives the disorientation, i.e., the smallest rotation in the set T 1 T B 2   g j B 2 T B 2   T . The third and fourth columns give the B2 symmetries g j B 2 in the double cosets by their matrices and by symmetry elements, respectively. The elements in bold are two-fold symmetries.
B2 Symmetries in the Double Cosets
Disorient.MatricesGeometrical Elements
O 0 I ( 1 0 0 0 1 0 0 0 1 ) ( 1 0 0 0 1 0 0 0 1 ) ( 0 1 0 1 0 0 0 0 1 ) ( 0 1 0 1 0 0 0 0 1 ) I , I ¯ , m ( 110 ) B 2 , R π , [ 110 ] B 2
O 2 R π , [ 001 ] B 19 ( 1 0 0 0 1 0 0 0 1 ) ( 0 1 0 1 0 0 0 0 1 ) ( 1 0 0 0 1 0 0 0 1 ) ( 0 1 0 1 0 0 0 0 1 ) m ( 001 ) B 2 , m ( 1 1 ¯ 0 ) B 2 , R π , [ 001 ] B 2 , R π , [ 1 1 ¯ 0 ] B 2
O 4 R 2 π / 3 , ~ [ 17 , 0 , 16 ] B 19 ( 1 0 0 0 0 1 0 1 0 ) ( 0 0 1 0 1 0 1 0 0 )   ( 0 0 1 1 0 0 0 1 0 )   ( 0 1 0 0 0 1 1 0 0 ) m ( 01 1 ¯ ) B 2 , m ( 101 ) B 2 , R ¯ 2 π / 3 , [ 1 11 ¯ ] B 2 , R ¯ 2 π / 3 , [ 1 11 ¯ ] B 2
( 1 0 0 0 0 1 0 1 0 )   ( 0 0 1 0 1 0 1 0 0 ) ( 0 0 1 1 0 0 0 1 0 )   ( 0 1 0 0 0 1 1 0 0 ) R π , [ 01 1 ¯ ] B 2 , R π , [ 101 ] B 2 , R 2 π / 3 , [ 1 11 ¯ ] B 2 , R 2 π / 3 , [ 1 11 ¯ ] B 2
O 5 R 2 π / 3 , ~ [ 40 3 ¯ ] B 19 ( 1 0 0 0 0 1 0 1 0 ) ( 0 0 1 0 1 0 1 0 0 )   ( 0 0 1 1 0 0 0 1 0 )   ( 0 1 0 0 0 1 1 0 0 ) m ( 011 ) B 2 , m ( 1 ¯ 01 ) B 2 , R ¯ 2 π / 3 , [ 1 ¯ 1 1 ¯ ] B 2 , R ¯ 2 π / 3 , [ 1 ¯ 1 1 ¯ ] B 2
( 1 0 0 0 0 1 0 1 0 )   ( 0 0 1 0 1 0 1 0 0 ) ( 0 1 0 0 0 1 1 0 0 ) ( 0 0 1 1 0 0 0 1 0 ) R π , [ 011 ] B 2 , R π , [ 1 ¯ 01 ] B 2 , R 2 π / 3 , [ 1 ¯ 1 1 ¯ ] B 2 , R 2 π / 3 , [ 1 ¯ 1 1 ¯ ] B 2
O 6 R π / 2 , ~ [ 11 , 0 , 1 ] B 19   ( 1 0 0 0 1 0 0 0 1 ) ( 1 0 0 0 1 0 0 0 1 )   ( 0 1 0 1 0 0 0 0 1 )   ( 0 1 0 1 0 0 0 0 1 ) m ( 010 ) B 2 , m ( 100 ) B 2 , R ¯ π / 2 , [ 001 ] B 2 , R ¯ π / 2 , [ 001 ] B 2
( 1 0 0 0 1 0 0 0 1 )   ( 1 0 0 0 1 0 0 0 1 ) ( 0 1 0 1 0 0 0 0 1 )   ( 0 1 0 1 0 0 0 0 1 ) R π , [ 010 ] B 2 ,   R π , [ 100 ] B 2 ,   R π / 2 , [ 001 ] B 2 ,   R π / 2 , [ 001 ] B 2
O 1 R π / 2 , ~ [ 0 , 8 , 7 ] B 19 ( 1 0 0 0 0 1 0 1 0 )   ( 0 0 1 0 1 0 1 0 0 ) ( 0 0 1 1 0 0 0 1 0 )   ( 0 1 0 0 0 1 1 0 0 ) R π / 2 , [ 100 ] B 2 , R π / 2 , [ 010 ] B 2 , R 2 π / 3 , [ 111 ] B 2 , R 2 π / 3 , [ 11 ¯ 1 ] B 2
( 1 0 0 0 0 1 0 1 0 )   ( 0 0 1 0 1 0 1 0 0 )   ( 0 0 1 1 0 0 0 1 0 )   ( 0 1 0 0 0 1 1 0 0 ) R ¯ π / 2 , [ 100 ] B 2 , R ¯ π / 2 , [ 010 ] B 2 , R ¯ 2 π / 3 , [ 111 ] B 2 , R ¯ 2 π / 3 , [ 11 ¯ 1 ] B 2
O 3 R π / 2 , ~ [ 0 , 8 , 7 ] B 19 ( 1 0 0 0 0 1 0 1 0 )   ( 0 0 1 0 1 0 1 0 0 ) ( 0 1 0 0 0 1 1 0 0 )   ( 0 0 1 1 0 0 0 1 0 ) R π / 2 , [ 100 ] B 2 , R π / 2 , [ 010 ] B 2 , R 2 π / 3 , [ 111 ] B 2 , R 2 π / 3 , [ 11 ¯ 1 ] B 2
( 1 0 0 0 0 1 0 1 0 )   ( 0 0 1 0 1 0 1 0 0 ) ( 0 1 0 0 0 1 1 0 0 ) ( 0 0 1 1 0 0 0 1 0 ) R ¯ π / 2 , [ 100 ] B 2 , R ¯ π / 2 , [ 010 ] B 2 , R ¯ 2 π / 3 , [ 111 ] B 2 , R ¯ 2 π / 3 , [ 11 ¯ 1 ] B 2
Table 3. Junction planes and twin characteristics predicted by the inter-correspondence operator O i between the B19′ variants. We indicate, for the type I twins, the B2 reflection plane and the B19′ twin plane (junction plane) and, for the type II twins, the B2 180° axis and the B19′ twin axis. The invariant but irrational plane K2 is written by a rational approximate tagged by a symbol ~ . The arrows mean both “correspondence” and “parallelism”. The junction planes that agree with the observations (Section 5.4) are written in bold.
Table 3. Junction planes and twin characteristics predicted by the inter-correspondence operator O i between the B19′ variants. We indicate, for the type I twins, the B2 reflection plane and the B19′ twin plane (junction plane) and, for the type II twins, the B2 180° axis and the B19′ twin axis. The invariant but irrational plane K2 is written by a rational approximate tagged by a symbol ~ . The arrows mean both “correspondence” and “parallelism”. The junction planes that agree with the observations (Section 5.4) are written in bold.
Reflection Junction Plane
(Type I)
Rotational Junction Plane
(Type II and Weak)
Shear Amplitude
O 2   ( 001 ) B 2   C * ( 100 ) B 19 or ( 1 1 ¯ 0 ) B 2 C * ( 001 ) B 19 s = 0.2389
O 4 ( 01 1 ¯ ) B 2 C * ( 1 ¯ 11 ) B 19 [ 01 1 ¯ ] B 2 C   [ 2 ¯ 11 ] B 19   i n v .   p l a n e   ~ ( 1 2 ¯ 4 ) B 19 w e a k   p l a n e :   ( 1 2 ¯ 4 ) B 19 s = 0.3096
O 5 ( 011 ) B 2 C * ( 111 ) B 19 [ 011 ] B 2 C [ 211 ] B 19 i n v .   p l a n e ~ ( 2 ¯ 13 ) B 19 w e a k   p l a n e s : ( 2 ¯ 13 ) B 19 ,   ( 1 ¯ 11 ) B 19 ,   ( 1 ¯ 02 ) B 19 s = 0.1422
O 6 ( 010 ) B 2 C * ( 011 ) B 19 [ 100 ] B 2 C [ 011 ] B 19 i n v .   p l a n e ~ ( 34 4 ¯ ) B 19 w e a k   p l a n e : ( 34 4 ¯ ) B 19 s = 0.2804
O 1 O 3 no solution [ 100 ] B 2 C [ 011 ] B 19 n o   i n v a r i a n t   p l a n e w e a k   p l a n e s :   ( 1 1 ¯ 1 ) B 19 ( 11 1 ¯ ) B 19   ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 ( 0 1 ¯ 1 ) B 19 ( 100 ) B 19 [ 111 ] B 2 C [ 110 ] B 19 n o   i n v a r i a n t   p l a n e w e a k   p l a n e :   n o   s o l u t i o n sg = 0.2911
Table 4. Closing-gap ORs between two B19′ variants linked by the operator O i predicted by the correspondence theory. In the last column are the deviation angles between the closing-gap OR and the natural OR AQ [24]. For each operator, the minimum angle is written in bold.
Table 4. Closing-gap ORs between two B19′ variants linked by the operator O i predicted by the correspondence theory. In the last column are the deviation angles between the closing-gap OR and the natural OR AQ [24]. For each operator, the minimum angle is written in bold.
Op.ORDeviation from AQ OR
O 2 case   1 :   ( 100 ) B 19     ( 001 ) B 2   and   [ 001 ] B 19     [ 1 ¯ 10 ] B 2 1.43°
case   2 :   ( 001 ) B 19     ( 1 1 ¯ 0 ) B 2   and   [ 100 ] B 19     [ 001 ] B 2 8.23°
O 4 type   I :   ( 1 ¯ 11 ) B 19       ( 01 1 ¯ )B2 and [ 7 ,   6 ¯ , 13 ] B 19   [ 19 ¯ , 7 , 7 ] B 2 2.21°
type   II :   [ 2 ¯ 11 ] B 19     [ 01 1 ¯ ] B 19 and ( 1 , 2 ¯ , 4 ) B 19     ( 3 ¯ , 1 , 1 ) B 2 8.88°
O 5 type   I :   ( 111 ) B 19     ( 011 )B2 and [ 62 ¯ , 21 , 41 ] B 19     [ 10 ¯ , 31 , 31 ¯ ] B 2 8.11°
type   II :   [ 211 ] B 19     [ 011 ] B 2 and ( 2 ¯ , 1 ,   3 ) B 19     ( 1 , ¯ 2 , 2 ¯ ) B 2 2.15°
O 6 type   I :   ( 011 ) B 19     ( 010 )B2 and [ 11 ¯ , 7 ¯ , 7 ] B 19     [ 14 ¯ , 0 ,   11 ¯ ] B 2 7.93°
type   II :   [ 011 ] B 19     [ 010 ] B 2 and ( 3 , 4 , 4 ¯ ) B 19     ( 4 , 0 , 3 ) B 2 3.94°
O 1   and   O 3 case 1: [ 011 ] B 19   [ 010 ] B 2   and   ( 1 3 ¯ 3 ) B 19     ( 3 ¯ 01 ) B 2 7.04°
case 2:   [ 011 ] B 19   [ 010 ] B 2 and ( 31 1 ¯ ) B 19     ( 103 ) B 2 3.45°
Table 5. Comparison between the CT and the PTMC to predict the B19′ junction planes associated with the operators O i . The letters t.I, t.II and C mean “type I”, “type II” and “compound” twins, respectively. In the column “Name in PTMC”, notation with the letter “S” was used by Hane and Shield [51] and notation with a letter (A, B, etc.) by Bhattacharya [2]. For the correspondence theory, the letter W means “weak twin”. The junction planes experimentally observed are written in bold.
Table 5. Comparison between the CT and the PTMC to predict the B19′ junction planes associated with the operators O i . The letters t.I, t.II and C mean “type I”, “type II” and “compound” twins, respectively. In the column “Name in PTMC”, notation with the letter “S” was used by Hane and Shield [51] and notation with a letter (A, B, etc.) by Bhattacharya [2]. For the correspondence theory, the letter W means “weak twin”. The junction planes experimentally observed are written in bold.
Closing-Gap
Disorientation
Junction Planes CTShearJunction Planes PTMCName in PTMCShear
O 2 180 ° , [ 001 ] B 19
(from C1)
C1: ( 100 ) B 19 , C2: ( 001 ) B 19 0.2385C1: ( 100 ) B 19 , C2: ( 001 ) B 19 S1, A0.2385
O 4 116 ° , ~ [ 1 , 0 , 1 ] B 19
(from I)
t.I: ( 1 ¯ 11 ) B 19
t.II: ~   ( 124 ) B 19
W: ( 124 ) B 19
0.3096t.I: ( 1 ¯ 11 ) B 19
t.II: ( 0.2470 ,   0.5060 ,   1 ) B 19
S3, D/C0.3096
O 5 119 ° , ~ [ 5 , 0 , 4 ¯ ] B 19
(from W)
t.I: ( 111 ) B 19
t.II: ~   ( 2 ¯ 13 ) B 19
W: ( 2 ¯ 13 ) B 19
0.1422t.I: ( 111 ) B 19
t.II: ( 0.6686 ¯ ,   0.3376 ,   1 ) B 19
S4, C/D0.1423
O 6 96 ° , ~ [ 13 , 0 , 1 ] B 19
(from W)
t.I: ( 011 ) B 19
t.II: ~ ( 34 4 ¯ ) B 19
W: ( 34 4 ¯ ) B 19
0.2804t.I: ( 011 ) B 19
t.II: ( 0.7206 ,   1 ,   1 ¯ ) B 19
S2, B0.2804
O 1 84 ° , ~ [ 0 , 1 , 1 ] B 19 W: ( 1 1 ¯ 1 ) B 19 ( 11 1 ¯ ) B 19   ( 1 3 ¯ 3 ) B 19 ( 31 1 ¯ ) B 19 ( 0 1 ¯ 1 ) B 19 ( 100 ) B 19 0.2912no solutionS5no value
O 3 84 ° , ~ [ 0 , 1 , 1 ] B 19 S6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cayron, C. The Correspondence Theory and Its Application to NiTi Shape Memory Alloys. Crystals 2022, 12, 130. https://doi.org/10.3390/cryst12020130

AMA Style

Cayron C. The Correspondence Theory and Its Application to NiTi Shape Memory Alloys. Crystals. 2022; 12(2):130. https://doi.org/10.3390/cryst12020130

Chicago/Turabian Style

Cayron, Cyril. 2022. "The Correspondence Theory and Its Application to NiTi Shape Memory Alloys" Crystals 12, no. 2: 130. https://doi.org/10.3390/cryst12020130

APA Style

Cayron, C. (2022). The Correspondence Theory and Its Application to NiTi Shape Memory Alloys. Crystals, 12(2), 130. https://doi.org/10.3390/cryst12020130

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop