### Re: Anisotropic peak broadening with TOPAS

```
Dear Topas users,

with a (probably clumsy) solution to achieve my
goal (see my original post far down below).

For those who are interested, I'll explain it
in relatively high detail, but before I do so,
I want to make some statements to prevent
unnecessary discussions.

The reasoning behind my doing is the following:
Investigating a large series of similar but
somehow variable samples, my goal is to derive
numerical parameters for each sample from its
powder XRD. Using these parameters, I can compare
and group samples, e.g. by making statements
like sample A and B are similar (or dissimilar)
with respect to this parameter. Thus, the primary
task is to parametrize the XRD results. Ideally,
such parameters would have some physical meaning,
like lattice parameters, crystallite size etc.
However, this does not necessarily mean that I
would interpret or trust parameters like e.g.
LVol-IB on an absolute scale!!! After all, it
is mainly relative trends I'm interested in.

LVol-IB is is one of the parameters I get and
tabulate if the peak broadening can be successfully
[For details on LVol-IB, see Topas (v3) Users
Manual, sections 3.4.1 and 3.4.2)]
If, however, the peak broadening is clearly
anisotropic, applying the isotropic model gives
inferior fit results. LVol-IB is still calculated,
but more or less meaningless.
Thus, I wanted an anisotropic fit model that BOTH
(a) yields a satisfactory fit AND (b) still
delivers parameters with a similar meaning as
the isotropic LVol-IB.

Applying a spherical harmonics function satisfied
condition (a), but not (b) (maybe just due to my
lack in mathematical insight).

Applying Peter Stephens' code (but modified for
size broadening) met condition (a) and brought
me halfway to reach condition (b). As I did not
find a way of teaching (coding) Topas to do
all calculations I wanted in the launch mode,
I developed a workaround to reach (b).

Now in detail:
The modified Stephens code I use looks like
this:
prm s400  29.52196`_2.88202 min 0
prm s040  40.52357`_4.10160 min 0
prm s004  6631.09739`_227.63909 min 0
prm s220  54.23582`_13.82762
prm s202  1454.83518`_489.04664
prm s022  5423.10499`_765.48349
prm mhkl = H^4 s400 + K^4 s040 + L^4 s004 + H^2 K^2 s220 + H^2
L^2 s202 + K^2 L^2 s022;

lor_fwhm = (D_spacing^2 * Sqrt(Max(0,mhkl)) / 1) / Cos(Th);
Compared to Peters original code, I have changed
the strain dependence * Tan(Th) into the size
dependence / Cos(Th) and re-arranged the remaining
terms in that line of code.
[Peter has mentioned that from the fundamental
theoretical point of view, spherical harmonics
might be better justified then his formula
for the case of SIZE broadening. Anyway,
it works for me from the practical point of
view, thus I'll use it.]
This rearrangement emphasizes the analogy between
the isotropic case c / Cos(Th) (where c is valid
for ALL peaks) and the anisotropic one, where
c is replaced by the hkl dependent term
(D_spacing^2 * Sqrt(Max(0,mhkl)) / 1).
Thus, I freely interpret this term as some
sort of c(hkl), which I will use for some
specific values of hkl to derive hkl dependent
analogues of LVol-IB.
The first step of this calculation I managed
to code for Topas based on Peters equations,
but for specifc hkl values:
prm ch00 = (Lpa^2 * Sqrt(s400) / 1); :  0.24382`_0.01190
prm c0k0 = (Lpb^2 * Sqrt(s040) / 1); :  0.45185`_0.02287
prm c00l = (Lpc^2 * Sqrt(s004) / 1); :  0.13114`_0.00225
As you can see, the c(hkl) term becomes
very simple for the cases hkl = 100, 010
and 001. For these reflections, the
D_spacing can be replaced with the
corresponding reserved variables for the
lattice parameters.

Now, I want to calculate something like
LVol-IB(hkl) from these c(hkl) values.
I still don't know how to code this in launch
mode, so I developed a workaround and switch
back to GUI mode.
In the GUI, I do the following:
1) I fix the lattice parameters to the
values obtained after running the anisotropic
fit in launch mode.
2) I deactivate the Cry Size L parameter on
the Structure tab and add a Lorentzian type
Convolutions tab instead. Both do the same
thing, but only the latter has the input format
I need.
3) I paste one of the c(hkl) values from
the *.out file of the anisotropic refinement,
e.g. the value calculated for ch00, into the
Value box of that convolution and keep
it fixed. Then I start the refinement.
Since lattice parameters and peak shape
are fixed, there is not much to be refined.
Not surprisingly, I obtain an ugly looking
fit in which only few reflections are
adequately described (e.g. the h00 series in
this case).
4) I use the Capture option to grab the
calculated pattern which is then inserted as
a new range. I export this range as ```

### Re: cif files for austenite

```
Dear Antonio,

as far as I know, austenite is just a solid solution
of a few mole-% of carbon in gamma-iron.
Thus, you essentially need the structure of gamma-Fe,
which is very simple, as it is an fcc metal
(space group Fm-3m, a = 3.6468 A, just one
Fe atom located on 0,0,0).
I found one ICSD entry for austenite.
However, the dissolved carbon is placed sharing
the Fe site in this structure, which I doubt makes
sense from the chemical point of view.
I would bet the carbon atoms are actually located
Anyway, as we are talking about a few MOLE-%
of carbon here, and carbon is much lighter than
iron, I think neglecting the carbon will not
make much of a difference (both for the accuracy
of the calculated relative intensities and
the quantification in wt-%.
Thus, I would just use gamma-Fe as a replacement
for austenite here.

Cheers,
Frank

--
Frank Girgsdies
Department of Inorganic Chemistry
Fritz Haber Institute (Max Planck Society)
--

antonio josé wrote:

Dear All.

I trying analyze phases in a steel sample using   X-Ray technique. The
analysis of X-Ray pattern shows that the main phases are ferrite and
austenite. The next step would be quantify these phases by Rietveld
method using TOPAS 3.0. I am looking for cif files for austenite and I
didn’t find it in WEB free database or in ICSD database.

Do anyone has this cif file or could tell me where I can find?

Thanks a lot,

Antonio José

[EMAIL PROTECTED]

```

### RE: Anisotropic peak broadening with TOPAS

```Hi Frank

I'm not 100% sure what you have been doing but I think the copying and
pasting to get LVol values for particular sets of hkls can be maybe done in
a more simple manner as shown below. If it's not clear then contact me off
the list.

prm lor_h00 300 min .3
prm lor_0k0 300 min .3
prm lor_00l 300 min .3
prm lor_hkl 300 min .3
prm gauss_h00 300 min .3
prm gauss_0k0 300 min .3
prm gauss_00l 300 min .3
prm gauss_hkl 300 min .3

prm = 1 / IB_from_CS(gauss_h00, lor_h00); : 0 ' This is LVol
prm = 0.89 / Voigt_FWHM_from_CS(gauss_h00, lor_h00); : 0 ' This is
LVol_FWHM

lor_fwhm =
(0.1 Rad Lam / Cos(Th)) /
IF And(K == 0, L == 0) THEN
lor_h00
ELSE IF And(H == 0, L == 0) THEN
lor_0k0
ELSE IF And(H == 0, K == 0) THEN
lor_00L
ELSE
lor_hkl
ENDIF
ENDIF
ENDIF
;

gauss_fwhm =
(0.1 Rad Lam / Cos(Th)) /
IF And(K == 0, L == 0) THEN
gauss_h00
ELSE IF And(H == 0, L == 0) THEN
gauss_0k0
ELSE IF And(H == 0, K == 0) THEN
gauss_00L
ELSE
gauss_hkl
ENDIF
ENDIF
ENDIF
;

Cheers
Alan

-Original Message-
From: Frank Girgsdies [mailto:[EMAIL PROTECTED]
Sent: Friday, 31 October 2008 11:29 PM
To: Frank Girgsdies; Rietveld_l@ill.fr
Subject: Re: Anisotropic peak broadening with TOPAS

Dear Topas users,

with a (probably clumsy) solution to achieve my
goal (see my original post far down below).

For those who are interested, I'll explain it
in relatively high detail, but before I do so,
I want to make some statements to prevent
unnecessary discussions.

The reasoning behind my doing is the following:
Investigating a large series of similar but
somehow variable samples, my goal is to derive
numerical parameters for each sample from its
powder XRD. Using these parameters, I can compare
and group samples, e.g. by making statements
like sample A and B are similar (or dissimilar)
with respect to this parameter. Thus, the primary
task is to parametrize the XRD results. Ideally,
such parameters would have some physical meaning,
like lattice parameters, crystallite size etc.
However, this does not necessarily mean that I
would interpret or trust parameters like e.g.
LVol-IB on an absolute scale!!! After all, it
is mainly relative trends I'm interested in.

LVol-IB is is one of the parameters I get and
tabulate if the peak broadening can be successfully
[For details on LVol-IB, see Topas (v3) Users
Manual, sections 3.4.1 and 3.4.2)]
If, however, the peak broadening is clearly
anisotropic, applying the isotropic model gives
inferior fit results. LVol-IB is still calculated,
but more or less meaningless.
Thus, I wanted an anisotropic fit model that BOTH
(a) yields a satisfactory fit AND (b) still
delivers parameters with a similar meaning as
the isotropic LVol-IB.

Applying a spherical harmonics function satisfied
condition (a), but not (b) (maybe just due to my
lack in mathematical insight).

Applying Peter Stephens' code (but modified for
size broadening) met condition (a) and brought
me halfway to reach condition (b). As I did not
find a way of teaching (coding) Topas to do
all calculations I wanted in the launch mode,
I developed a workaround to reach (b).

Now in detail:
The modified Stephens code I use looks like
this:
prm s400  29.52196`_2.88202 min 0
prm s040  40.52357`_4.10160 min 0
prm s004  6631.09739`_227.63909 min 0
prm s220  54.23582`_13.82762
prm s202  1454.83518`_489.04664
prm s022  5423.10499`_765.48349
prm mhkl = H^4 s400 + K^4 s040 + L^4 s004 + H^2 K^2 s220 +
H^2
L^2 s202 + K^2 L^2 s022;
lor_fwhm = (D_spacing^2 * Sqrt(Max(0,mhkl)) / 1) /
Cos(Th);
Compared to Peters original code, I have changed
the strain dependence * Tan(Th) into the size
dependence / Cos(Th) and re-arranged the remaining
terms in that line of code.
[Peter has mentioned that from the fundamental
theoretical point of view, spherical harmonics
might be better justified then his formula
for the case of SIZE broadening. Anyway,
it works for me from the practical point of
view, thus I'll use it.]
This rearrangement emphasizes the analogy between
the isotropic case c / Cos(Th) (where c is valid
for ALL peaks) and the anisotropic one, where
c is replaced by the hkl dependent term
(D_spacing^2 * Sqrt(Max(0,mhkl)) / 1).
Thus, I freely interpret this term as some
sort of c(hkl), which I will use for some
specific values of hkl to derive hkl dependent
analogues of LVol-IB.
The first step of this calculation I managed
to code for Topas based on Peters equations,
but for specifc hkl values:
prm ch00 = (Lpa^2 * Sqrt(s400) / 1); :  0.24382`_0.01190
```