Topo 440
Topo 440
Topo 440
exe):
Gougoutas’ original paper on the specification of topotactic alignments appeared in 1972, and presented
an elegant, straightforward methodology for specifying the orientation of mother and daughter phases.1 A later
methodology, and an Appendix containing the mathematical relationships was published by Etter in 1976.2
These two papers are essential reading for any scientist interested in establishing, by experiment, the alignment
of a mother phase and one or more daughter phases arising from a crystal-to-crystal solid-state transformation.
Despite their importance, the papers have been cited only 15 and 11 times, respectively. As we returned to the
field of topotactic transformations in 2010, we sought to use modern diffractometry to solve the problem.
It is important here to repeat the definition of topotaxy so that an important point may be addressed. In
topotaxy, a single crystal of a starting material is converted into a pseudomorph2 containing one or more
products in a definite crystallographic orientation; the conversion takes place throughout the entire volume of
the crystal.3 Thus, topotaxy pertains to either a single crystal-to-single crystal process (a solitary daughter phase
is produced) OR a process during which the mother phase is converted into more than one daughter (twinned
phases are produced). Thus, in order to keep the discussion general, we will term the processes described in the
manuscript as “crystal-to-crystal”. In fact, in the present work, all the transformations studied involve twinned
product phases. The important point is that, for topotactic transformations, the diffraction properties of the
daughter phase allow the orientation(s) of the mother and daughter phases(s) to be established experimentally,
In the 1970’s the alignment of mother and daughter phases in molecular crystals was established through
examination of aligned Weissenberg or precession photographs.1,2,4-6 The analysis led to the following
interlattice relationship, which is important in interpreting the topotaxy observed in the experiment (symbols
bi = φij aj
where φij is the topotactic transformation matrix, while bi and aj are the direct basis set vectors for the product
and reactant lattices, respectively. We refer the interested reader to Etter’s paper for other elements of the
theory, and an interesting historical account of how φij was derived from observations on X-ray photographs.3
In the present day, our approach first centered on the use of the original CAD-4 command set, as
implemented in the CAD-4PC software,7 which may still be run on a modern PC! For example, after entering
the orientation matrix (taken from a *.p4p file obtaind using our Bruker-Nonius Kappa-Apex2 instrument) for a
mother phase, we could use the HB command (hkl to ‘bisecting setting’ angles: θ, φ, χ) to get the setting angles
of, for example, 100, 010, and 001 for the mother phase. Then, after cooling our crystal, and observing
appearance of the daughter unit cell, we enter the orientation matrix of the daughter crystal into CAD4-PC, and
use the BH command (‘bisecting setting’ to hkl), along with the original 100, 010, and 001 angles settings from
the mother phase. The resultant values show the correspondence between mother and daughter vectors before
and after the transformation. While this worked beautifully, we imagined that the use of an older, though
elegant, program would discourage general adoption by the community. Further, we had not yet solved the
The orientation matrix for a crystal mounted on a diffractometer has been part of standard printed output
for many decades, and now appears in one or more files, either in one or both of *.cif or *.p4p files. Thus, we
decided to write a program that, in its initial incarnation, would read *.p4p files each of a mother and daughter
phase, automatically calculate angles between all nine combinations of 100, 010 and 001 mother vs. daughter
vectors, and further, allow (a) the calculation of the corresponding values for h, k, and l for a daughter phase
referenced to the mother phase; (b) the calculation of the corresponding values for h, k, and l for a mother phase
referenced to the daughter phase; and (c) the angle between a specific h, k, l combination in the mother phase
and any chosen h, k, l combination in the daughter phase. Finally, the program would provide the opportunity
to perform the calculations either in reciprocal space or direct space. After reading in two *.p4p files, the
program automatically calculates the diffractometer coordinates (cmx, cmy, cmz) for an axial hkl set (in reciprocal
or direct space) of the mother crystal, using Rm, the orientation matrix for the mother crystal:
cmx hm
cmy = Rm km
cmz l
m
After getting the coordinates, we then use the inverse of the orientation matrix for the daughter crystal, R-1d , to
obtain the correspondence between the indices of the mother crystal referenced to the daughter unit cell:
hd cmx
−1
kd = R d cmy
l cmz
d
With the above information in hand, one can obtain direction cosines, and calculate angles between mother &
daughter vectors as well. As we wrote the program in g95 Fortran 95,8 it became obvious that calculation of the
φij = R-1d Rm
As the source in our 60Co gamma-irradiation apparatus reached its eighth half-life, it became apparent that the
exciting work we’d carried out on radical reactions in the solid state10,11 could no longer be carried out at
Brandeis, and alternative sources were not available to us. We thus decided to return to studies of crystal-to-
crystal, thermally-induced reactions and/or phase transformations. At the outset, we faced a dilemma: in the
1970’s the alignment of reactant and product (‘mother’ and ‘daughter’) phases was established using X-ray
precession photographs, taken before and after the transformation. While a modern diffractometer produces no
shortage of images, one does not work with aligned samples, and thus the images cannot be used directly to
establish relative orientations. However, using the approach described by Gougoutas et al.1,4 and by Etter,2 it is
straightforward to carry out identical analyses to those done in the 1970’s. The technique is simple and its
application provides an opportunity to experimentally establish the orientation relationship between a mother
phase and one or more daughter phases or twin orientations. In a recent contribution we describe the
methodology, simple mathematics, and a program (available from the authors) that simplifies the calculation
steps and provides all necessary information to establish orientation relationships and provide an accurate
method to visually overlay mother and daughter packing diagrams.12 Taking a lead from Peggy Etter’s
presentation style, the paper combines the scientific results with a tutorial approach.2,12
1. Gougoutas, J. Z. Isr. J. Chem. 1972, 10, 395-407. In “conservative twinning”, the symmetry of the
daughter phase conserves the symmetry of the mother phase via twinning; the concept is fully illustrated
2. Etter, M. C. J. Am. Chem. Soc. 1976, 98, 5331-5339. Etter defines a pseudomorphic transformation as
“one in which the morphology of the reactant crystal does not change during reaction. In a single crystal
to single crystal pseudomorphic transformation, the diffraction pattern of the reacting crystal changes
but the dimensions and shape of the crystal itself does not (color and density changes often occur also).”
3. Dent Glasser, L. S.; Glasser, F. P.; Taylor, H. F. W. Quart. Rev. 1962, 16, 343-60.
4. (a) Gougoutas, J. Z.; Clardy, J. C. J. Solid State Chem. 1972, 4, 230-242; (b) Gougoutas, J. Z.;
Lessinger, L. J. Solid State Chem. 1973, 7, 175-185; (c) Gougoutas, J. Z.; Lessinger, L. J. Solid State
Chem. 1974, 9, 155-164; (d) Gougoutas, J. Z.; Lessinger, L. J. Solid State Chem. 1975, 12, 51-62; (e)
Gougoutas, J. Z.; Naae, D. G. J. Solid State Chem. 1976, 16, 271-281; (f) Gougoutas, J. Z., Chang, K.
H.; Etter, M. C. J. Solid State Chem. 1976, 16, 283-291; (g) Etter, M. C. J. Solid State Chem. 1976, 16,
399-411; (h) Etter, M. C. J. Am. Chem. Soc. 1976, 98, 5326-5331; (i); (j) Gougoutas, J. Z. J. Am. Chem.
Soc. 1977, 99, 127-132; (k) Gougoutas, J. Z.; Naae, D. G. J. Phys. Chem. 1978, 82, 393-401; (l)
Gougoutas, J. Z.; Johnson, J. J. Am. Chem. Soc. 1978, 100, 5816-5820; (m) Gougoutas, J. Z. J. Am.
6. Cheng, K.; Foxman, B. M.; Gersten, S. W. Mol. Cryst. Liq. Cryst. 1979, 52, 77-82; Foxman, B. M.;
7. Schagen, J. D.; Straver, L.; van Meurs, F.; Williams, G. CAD4 Version 5.0; Enraf-Nonius, Delft, The
http://www.nonius.nl/cad4/manuals/user/chapter02.html and
http://www.nonius.nl/cad4/manuals/user/chapter08.html .
8. Vaught, A. G95; Mesa, Arizona, 2006; see also http://www.g95.org/index.shtml .
9. Experience to date suggests that, from experiment to experiment, and compound to compound,
variations in the observed values of interplanar/interaxial angles may range from as little as 0.3º up to
ca. 2º. In part, this depends on (a) the mechanical stability of the mounted crystal; (b) the
quality/”singleness” of the crystal at the temperatures measured; and (c) the number of data used to
10. Hickey, M. B.; Schlam, R. F.; Guo, C.; Cho, T. H.; Snider, B. B.; Foxman, B. M. CrystEngComm 2011,
13, 3146-3155.
11. Shang, W.; Hickey, M. B.; Enkelmann, V.; Snider, B. B.; Foxman, B. M. CrystEngComm 2011, 13,
4339-4350.
12. Posner, S. R.; Lorson, L. C.; Gell, A. R.; Foxman, B. M. Cryst. Growth Des., 2015, 15, 3407–3416.
ANNOTATED INSTRUCTIONS FOR TOPO v4.40
Below we provide annotated instructions using TOPO. These pages, as well as the details provided in
reference 12, provide a complete picture of the extensive options available in the paper, as well as guidelines on
the preparation of animations that illustrate the molecular motion which occurs between the initial and final
stages of a transformation. In the following example, the transformation between the RT phase and one of the
four twin components of the LT phase of ferrocenium hexafluorophosphate is explored. The downloadable .zip
file at http://www.xray.chem.brandeis.edu contains the TOPO executable, input files and output for checking
program use. Do note that a correct experiment involves a crystal transformation carried out in a single
diffractometer session: the crystal may not be removed between data collection of Form I and the generation of
other phases!!!
---------------------------------------------------------------------------------------------------------------------------------
C:\Users\foxman1\Desktop\topo12\twinlaws>topo440
[ENTER Y OR N]: n
[DEFAULT(CR) = P4P] :
[1 = PARENT MATRIX...NO TWIN]: 2 (I’ll look at twin orientation 2 for the daughter; rerun fr 1, 3, and 4 later)
ENTER FILE NAME FOR OUTPUT [DEFAULT=TOPO.OUT]: topo440.out
Reciprocal: Part IR -- AUTOMATIC GENERATION OF ANGLES This provides an automatic look at some of
BETWEEN PRIMARY R. L. VECTORS the main intervector angles; things can jump
M = MOTHER, D= DAUGHTER out at us here!
HM KM LM HD KD LD ANGLE
1. 0. 0. 0. 0. 1. 93.646
0. 1. 0. 1. 0. 0. 90.261
0. 0. 1. 1. 0. 0. 97.204
0. 0. 1. 0. 1. 0. 89.405
Reciprocal: Part IIR(a) -- GENERATE DAUGHTER INDICES THAT CORRESPOND Let’s have a look at what the
TO INPUT MOTHER INDICES
mother 100 r.l. vector
NOTE THAT DAUGHTER INDICES MAY BE FRACTIONAL! corresponds to in the daughter
MOTHER DAUGHTER
H K L H K L
_______________________________________________________________
1.000 0.000 0.000 0.982 -0.002 -0.255 (mom’s 100 is daughter ‘s 4 0 -1)
MOTHER DAUGHTER
H K L H K L
_______________________________________________________________
M = MOTHER, D= DAUGHTER
ANOTHER? [Y OR N]: n
M = MOTHER, D= DAUGHTER
1. 0. 0. 1. 0. 0. 0.517
1. 0. 0. 0. 1. 0. 90.163
Very similar to the above input
1. 0. 0. 0. 0. 1. 97.207 and output, except that all this is
now in DIRECT SPACE.
0. 1. 0. 1. 0. 0. 90.262
0. 1. 0. 0. 1. 0. 179.189
0. 1. 0. 0. 0. 1. 90.555
0. 0. 1. 1. 0. 0. 93.643
0. 0. 1. 0. 1. 0. 89.200
0. 0. 1. 0. 0. 1. 169.575
TYPE<CR> to CONTINUE:
MOTHER DAUGHTER
[H K L] [H K L]
_______________________________________________________________
M = MOTHER, D= DAUGHTER
HM KM LM HD KD LD ANGLE
ANOTHER? [Y OR N]: n
(H K L) [H K L]
_______________________________________________________________ Now we check the obliquity for this
TLQS twin, a reasonable 3.3°.
0. 0. 1. 1. 0. 4. 3.248