Stephens Peakshapes for TOF NPD

jmas5

I've been working on some combined SXRD/NPD refinements with data from Diamond (I-11) and WISH (ISIS). I have always found the Stephens peakshape macros to be incredibly helpful for XRD, but have not been able to use these for my NPD. I believe the issue is due to the TOF (microsecond) data in NPD, as opposed to the 2theta data in XRD, as the macro calls "th":

macro Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta)
{
prm mhkl = H⁴ s400 + K⁴ s040 + L⁴ s004 +
H² K² s220 + H² L² s202 + K² L² s022 +
H K² L s121 +
H L³ s103 + H³ L s301;

prm pp = D_spacing² * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) Tan(Th) + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta Tan(Th) + 0.0001;
}

We are presently trying to edit this macro to be compatible with our TOF data, but have not had any success. We are effectively trying to combine the TOF_x_axis_calibration macro with Stephens.

As WISH is multibank (i.e. each bank has its own set of t0, t1, t2 e.t.c), which are fit as part of the standard refinement (NaCaAlF) (and then t1/t2 against the sample in later stages), it isn't an option for us to just convert the data from TOF to 2theta, then read that in. (We're not entirely certain of what the implications of this are, and whether this would result in a high degree of correlation between those refined parameters, and those exclusive to this macro...)

This is the current amalgamation of macros:

macro TOF_Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta, t0, t0v, t1, t1v, t2, t2v)
{

prm mhkl = H⁴ s400 + K⁴ s040 + L⁴ s004 +
H² K² s220 + H² L² s202 + K² L² s022 +
H K² L s121 +
H L³ s103 + H³ L s301;

prm pp = D_spacing² * Sqrt(Max(mhkl,0)) / 1000;

pk_xo = CeV(t0,t0v) + CeV(t1,t1v) D_spacing + CeV(t2,t2v) D_spacing²;

'note flightpath, t0, t1, t2 defined locally for each bank's xdd
prm !mytheta = ArcSin(pk_xo/505.56flightpathD_spacing);

gauss_fwhm = 1.8/3.1415927 pp (1-eta) Tan(mytheta) + 0.0001;
low_fwhm = 1.8/3.1415927 pp eta Tan(mytheta) + 0.0001;
}

Which is giving the following error:

"Cannot find match for macro TOF_Stephens_monoclinic
Number of arguements 20"

I also tried this with the whole TOF_x_axis_calibration macro in there, but got the same error.

Has anybody tried this before/can explain the errors I'm making above? Apologies if this is an obvious fix, I am not as familiar with the syntax of TOPAS as I would like to be!

Thanks in advance,
James

jmas5

Have realised that the formatting above automatically italicises anything inside *'s. Where above random lines are in italic, that represents: asterisk xyz asterisk

rowlesmr

I'm not commenting your amalgamation, as I'm not qualified.

How are you calling your TOF macro? It has 16 arguments. The error

"Cannot find match for macro TOF_Stephens_monoclinic
Number of arguements 20"

says that you're calling it in your code with 20 arguments, and TOPAS can't find a macro called "TOF_Stephens_monoclinic" which takes 20 arguments.

johnsoevans

In tof you can just use D_spacing which topas will know for each hkl instead of tan Th. You don’t need to send over calibration constants. This is covered (more or less) in the last section of chapter 4 in the topas/Rietveld text book.

jmas5

Hi Both,

Thanks for your help here. I have modified the original stephens macro to replace tan Th with D_spacing as suggested:

macro TOF_Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta)
{
prm mhkl = H4 s400 + K4 s040 + L4 s004 +
H2 K2 s220 + H2 L2 s202 + K2 L2 s022 +
H K2 L s121 +
H L3 s103 + H3 L s301;

prm pp = D_spacing2 * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
}

Which I am calling as such in the refinement (note for the XRD I call "Stephens_monoclinic" in this same way, with no issues):

TOF_Stephens_monoclinic
(
eta, 0.5,
s400, 200,
s040, 200,
s004, 200,
s220, 200,
s202, 200,
s022, 200,
s121, 200,
s103, 200,
s301, 200
)

However, I'm still receiving the same error message:

"Cannot find match for macro TOF_Stephens_monoclinic
Number of arguements 20"

As suggested above, it is saying that it cannot find this macro, which takes 20 arguments (corresponding to the 10 parameter names and 10 parameter values in the called macro).

I'm defining this, as I would any temporary macro, at the top of the file (as opposed to in the .inc). This works just fine for the other macros I have defined this way, such as the below:

macro Out_Tick_XRD_Q(file)
{
phase_out file load out_record out_fmt out_eqn
{
" %11.5f 200\n" = 2 *Pi / D_spacing;
}
}

Is there any reason why TOPAS would not be recognising this user defined macro? I have tried changing the name in case there is anything special about it, but I receive the same error, only with the new name. Functionally, the only differences between the macro I am trying to use (unsuccessfully), and the default stephens_monoclinic, is the new name, and the replacement of tan Th with D_spacing.

Checking the stephens_monoclinic macro in our topas.inc (V7) actually shows it is a bit different to that on the topas website:

macro Stephens_monoclinic
(
etac, etav,
s400c,s400v,
s040c,s040v,
s004c,s004v,
s220c,s220v,
s202c,s202v,
s022c,s022v,
s121c,s121v,
s103c,s103v,
s301c,s301v
)
{
#m_argu s400c If_Prm_Eqn_Rpt(s400c, s400v,)
#m_argu s040c If_Prm_Eqn_Rpt(s040c, s040v,)
#m_argu s004c If_Prm_Eqn_Rpt(s004c, s004v,)
#m_argu s220c If_Prm_Eqn_Rpt(s220c, s220v,)
#m_argu s202c If_Prm_Eqn_Rpt(s202c, s202v,)
#m_argu s022c If_Prm_Eqn_Rpt(s022c, s022v,)
#m_argu s121c If_Prm_Eqn_Rpt(s121c, s121v,)
#m_argu s103c If_Prm_Eqn_Rpt(s103c, s103v,)
#m_argu s301c If_Prm_Eqn_Rpt(s301c, s301v,)
Stephens_lor_gauss(
etac, etav,
CeV(s400c,s400v) H⁴ +
CeV(s040c,s040v) K⁴ +
CeV(s004c,s004v) L⁴ +
CeV(s220c,s220v) H² K² +
CeV(s202c,s202v) H² L² +
CeV(s022c,s022v) K² L² +
CeV(s121c,s121v) H K² L +
CeV(s103c,s103v) H L³ +
CeV(s301c,s301v) H³ L
)
}

I am a little confused, as this seems to entirely remove reference to the tan Th, which we thought was causing the TOF incompatibility to begin with. Apologies for the wall of text, if you have any suggestions for why this might still not be working, any advice greatly appreciated!

Thanks again for your time and help!

rowlesmr

jmas5 Checking the stephens_monoclinic macro in our topas.inc (V7) actually shows it is a bit different to that on the topas website:

macro Stephens_monoclinic
...
}

I am a little confused, as this seems to entirely remove reference to the tan Th, which we thought was causing the TOF incompatibility to begin with.

It just abstracts away the Tan(Th) to the Stephens_lor_gauss macro, which is

macro Stephens_lor_gauss(etac, etav, mhkl)
{
   #m_argu etac  If_Prm_Eqn_Rpt(etac, etav, min 0 max 1)   
               local #m_unique pp_ = D_spacing^2 * Sqrt(Max(mhkl,0)) Tan(Th) 0.0018/3.1415927;
               gauss_fwhm =  pp_ (1-CeV(etac, etav));
               lor_fwhm   =  pp_ CeV(etac, etav);
}

johnsoevans

I think you're just missing a #m_argu for (eta,etav) which is in your macro call.

jmas5

Hi Both, thanks again for your help and patience with me on this!

johnsoevans For John's comment, I think that's correct, but for some reason that macro works just fine for the xray xdds and has never had an issue. I think possibly because it is defined in the "Stephens_lor_gauss" macro that it calls? For consistency, I have tried making the below macro, where I have replaced references to tan Th with D_spacing:

Macro:

macro TOF_Stephens_lor_gauss(etac, etav, mhkl)
{
#m_argu etac If_Prm_Eqn_Rpt(etac, etav, min 0 max 1)
local #m_unique pp_ = D_spacing² * Sqrt(Max(mhkl,0)) D_spacing 0.0018/3.1415927;
gauss_fwhm = pp_ (1-CeV(etac, etav));
lor_fwhm = pp_ CeV(etac, etav);
}

macro TOF_Stephens_monoclinic
(
etac, etav,
s400c,s400v,
s040c,s040v,
s004c,s004v,
s220c,s220v,
s202c,s202v,
s022c,s022v,
s121c,s121v,
s103c,s103v,
s301c,s301v
)
{
#m_argu s400c If_Prm_Eqn_Rpt(s400c, s400v,)
#m_argu s040c If_Prm_Eqn_Rpt(s040c, s040v,)
#m_argu s004c If_Prm_Eqn_Rpt(s004c, s004v,)
#m_argu s220c If_Prm_Eqn_Rpt(s220c, s220v,)
#m_argu s202c If_Prm_Eqn_Rpt(s202c, s202v,)
#m_argu s022c If_Prm_Eqn_Rpt(s022c, s022v,)
#m_argu s121c If_Prm_Eqn_Rpt(s121c, s121v,)
#m_argu s103c If_Prm_Eqn_Rpt(s103c, s103v,)
#m_argu s301c If_Prm_Eqn_Rpt(s301c, s301v,)
TOF_Stephens_lor_gauss(
etac, etav,
CeV(s400c,s400v) H⁴ +
CeV(s040c,s040v) K⁴ +
CeV(s004c,s004v) L⁴ +
CeV(s220c,s220v) H² K² +
CeV(s202c,s202v) H² L² +
CeV(s022c,s022v) K² L² +
CeV(s121c,s121v) H K² L +
CeV(s103c,s103v) H L³ +
CeV(s301c,s301v) H³ L
)
}

This yields the below error, when called like this (with the 20 arguments):

Call:

TOF_Stephens_monoclinic
(
etac, 0.5,
s400c, 200,
s040c, 200,
s004c, 200,
s220c, 200,
s202c, 200,
s022c, 200,
s121c, 200,
s103c, 200,
s301c, 200
)

Error:

" Uninitialized_Variable in equation: m672c9d7a_1b"

rowlesmr Taking Matt's approach for the below macro, is also resulting in a different error currently:

Macro:

prm pp = D_spacing2 * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
}

Call:

prm eta 0.5
prm s400 200
prm s040 200
prm s004 200
prm s220 200
prm s202 200
prm s022 200
prm s121 200
prm s103 200
prm s301 200
TOF_Stephens_monoclinic(eta, s400, s040, s004, s220, s202, s022, s121, s103, s301)

Error:
"Uninitialized_Variable in equation: K2"

I wonder if the error is the same between them both, just that it is named in the second macro and not in the first? Thanks again for your help!

rowlesmr

So you have defined a macro

macro TOF_Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta)
{
prm mhkl = H4 s400 + K4 s040 + L4 s004 +
H2 K2 s220 + H2 L2 s202 + K2 L2 s022 +
H K2 L s121 +
H L3 s103 + H3 L s301;

prm pp = D_spacing2 * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
}

which you are then calling in your input file as

TOF_Stephens_monoclinic(eta, 0.5, s400, 200, s040, 200, s004, 200, s220, 200, s202, 200, s022, 200, s121, 200, s103, 200, s301, 200)

Is that correct?

Well, there's your problem. You've defined a macro with 10 arguments and are calling it with 20. TOPAS then complains that it can't find a macro definition with 20 arguments because you haven't provided one.

You need to call it with 10 arguments.

prm eta 0.5
prm s400 200
prm s040 200
prm s004 200
prm s220 200
prm s202 200
prm s022 200
prm s121 200
prm s103 200
prm s301 200
TOF_Stephens_monoclinic(eta, s400, s040, s004, s220, s202, s022, s121, s103, s301)

jmas5

(Edit: re-posting with correct formatting)

Hi Both, thanks again for your help and patience with me on this!

Macro:

macro TOF_Stephens_lor_gauss(etac, etav, mhkl)
{
#m_argu etac If_Prm_Eqn_Rpt(etac, etav, min 0 max 1)
local #m_unique pp_ = D_spacing2 * Sqrt(Max(mhkl,0)) D_spacing 0.0018/3.1415927;
gauss_fwhm = pp_ (1-CeV(etac, etav));
lor_fwhm = pp_ CeV(etac, etav);
}

macro TOF_Stephens_monoclinic
(
etac, etav,
s400c,s400v,
s040c,s040v,
s004c,s004v,
s220c,s220v,
s202c,s202v,
s022c,s022v,
s121c,s121v,
s103c,s103v,
s301c,s301v
)
{
#m_argu s400c If_Prm_Eqn_Rpt(s400c, s400v,)
#m_argu s040c If_Prm_Eqn_Rpt(s040c, s040v,)
#m_argu s004c If_Prm_Eqn_Rpt(s004c, s004v,)
#m_argu s220c If_Prm_Eqn_Rpt(s220c, s220v,)
#m_argu s202c If_Prm_Eqn_Rpt(s202c, s202v,)
#m_argu s022c If_Prm_Eqn_Rpt(s022c, s022v,)
#m_argu s121c If_Prm_Eqn_Rpt(s121c, s121v,)
#m_argu s103c If_Prm_Eqn_Rpt(s103c, s103v,)
#m_argu s301c If_Prm_Eqn_Rpt(s301c, s301v,)
TOF_Stephens_lor_gauss(
etac, etav,
CeV(s400c,s400v) H4 +
CeV(s040c,s040v) K4 +
CeV(s004c,s004v) L4 +
CeV(s220c,s220v) H2 K2 +
CeV(s202c,s202v) H2 L2 +
CeV(s022c,s022v) K2 L2 +
CeV(s121c,s121v) H K2 L +
CeV(s103c,s103v) H L3 +
CeV(s301c,s301v) H3 L
)
}

This yields the below error, when called like this (with the 20 arguments):

Call:

TOF_Stephens_monoclinic
(
etac, 0.5,
s400c, 200,
s040c, 200,
s004c, 200,
s220c, 200,
s202c, 200,
s022c, 200,
s121c, 200,
s103c, 200,
s301c, 200
)

Error:

" Uninitialized_Variable in equation: m672c9d7a_1b"

rowlesmr Taking Matt's approach for the below macro, is also resulting in a different error currently:

Macro:

macro TOF_Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta)
{
prm mhkl = H4 s400 + K4 s040 + L4 s004 +
H2 K2 s220 + H2 L2 s202 + K2 L2 s022 +
H K2 L s121 +
H L3 s103 + H3 L s301;

prm pp = D_spacing2 * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
}

Call:

prm eta 0.5
prm s400 200
prm s040 200
prm s004 200
prm s220 200
prm s202 200
prm s022 200
prm s121 200
prm s103 200
prm s301 200
TOF_Stephens_monoclinic(eta, s400, s040, s004, s220, s202, s022, s121, s103, s301)

Error:
"Uninitialized_Variable in equation: K2"

I wonder if the error is the same between them both, just that it is named in the second macro and not in the first? Thanks again for your help!

rowlesmr

Its a copy/paste forum-formatting issue (at least for the error in my macro)

macro TOF_Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta)
{
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 +
H K^2 L s121 +
H L^3 s103 + H^3 L s301;

prm pp = D_spacing^2 * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
}

You need the sysmbol in there to make a power.

jmas5

Thanks for this! That does seem to have been the issue with that macro. However, I'm now getting a different error (though I think this now means that the macro itself is working). The new fun error is "Negative FWHM encountered", which I have seen before when I accidentally leave these parameters refining with simulated annealing on. I'm wondering if perhaps the starting parameters that the macro for XRD uses will now need to be different for NPD? To confirm, I'm currently using:

Macro:

`macro TOF_Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta)
{
prm mhkl = H⁴ s400 + K⁴ s040 + L⁴ s004 +
H² K² s220 + H² L² s202 + K² L² s022 +
H K² L s121 +
H L³ s103 + H³ L s301;

prm pp = D_spacing² * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
}
`

Input file:
prm eta 0.5 prm s400 200 prm s040 200 prm s004 200 prm s220 200 prm s202 200 prm s022 200 prm s121 200 prm s103 200 prm s301 200 TOF_Stephens_monoclinic(eta, s400, s040, s004, s220, s202, s022, s121, s103, s301)

Setting all the parameters to 1 (instead of 200) allowed the refinement to initialise, but after one run I encountered the negative FWHM error again. I was able to get around this with the below restrictions, however all the parameters do minimise. Is this physically meaningful?

prm eta_b5 0.5 min 0 max 1 prm s400_b5 0.01 min 0 max 1 prm s040_b5 0.01 min 0 max 1 prm s004_b5 0.01 min 0 max 1 prm s220_b5 0.01 min 0 max 1 prm s202_b5 0.01 min 0 max 1 prm s022_b5 0.01 min 0 max 1 prm s121_b5 0.01 min 0 max 1 prm s103_b5 0.01 min 0 max 1 prm s301_b5 0.01 min 0 max 1 TOF_Stephens_monoclinic(eta_b5, s400_b5, s040_b5, s004_b5, s220_b5, s202_b5, s022_b5, s121_b5, s103_b5, s301_b5)

Thanks again for your continued help!

jmas5

Apologies, macro was not formatted as code. Re-posted below:

`macro TOF_Matt_Stephens_monoclinic(s400, s040, s004, s220, s202, s022, s301, s121, s103, eta)
{
prm mhkl = H⁴ s400 + K⁴ s040 + L⁴ s004 +
H² K² s220 + H² L² s202 + K² L² s022 +
H K² L s121 +
H L³ s103 + H³ L s301;

prm pp = D_spacing² * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
} `

Also for reference, below are the refined values:

prm eta_b5 0.0625563.501034_LIMIT_MIN_0 min 0 max 1
prm s400_b5 0.87625_4709.88862_LIMIT_MIN_0 min 0 max 1 prm s040_b5 0.876254764.84687_LIMIT_MIN_0 min 0 max 1
prm s004_b5 0.12625_875.637423_LIMIT_MIN_0 min 0 max 1 prm s220_b5 0.626254261.83885_LIMIT_MIN_0 min 0 max 1
prm s202_b5 0.87625_4728.59797_LIMIT_MIN_0 min 0 max 1 prm s022_b5 0.876255195.35536_LIMIT_MIN_0 min 0 max 1
prm s121_b5 0.87625_4906.25752_LIMIT_MIN_0 min 0 max 1 prm s103_b5 0.37625_2631.16776_LIMIT_MIN_0 min 0 max 1
prm s301_b5 0.87625_2375.1912_LIMIT_MIN_0 min 0 max 1

jmas5

I do not know why it is not formatting correctly, but the "^" signs are in the macro as below

johnsoevans

The Max(mhkl, 0) should stop pp from going negative. Have you tried just putting the min/max on the eta Gaussian/Lorentzian broadening term? That should be the only term left that would allow a negative fwhm (unless you are also refining other peak shape parameters).

[The code formatting on the wiki isn't topas-specific. It's just something built in that "almost works". Sorry for hassle.]

jmas5

johnsoevans

Hi John,

Thanks for your message. Trying to use the standard values of 200 (Also happens for 20, 2) for the sx terms, with only min 0 max 1 (start 0.5) on eta, still results in the negative FWHM error. I can get it start the refinement by setting them all to 1 (or 0.5/0.05...), but I quickly encounter negative FHWM error on the first step. It seems the only way I can make the refinement happy is to give it min 0 max 1 for all parameters, below are my starting point and the refined values. If this seems physically reasonable I am happy to use this, but would like to avoid adding extra terms otherwise. Thanks again for your help!

Macro:

prm pp = D_spacing² * Sqrt(Max(mhkl,0)) / 1000;

gauss_fwhm = 1.8/3.1415927 pp (1-eta) D_spacing + 0.0001;
lor_fwhm = 1.8/3.1415927 pp eta D_spacing + 0.0001;
} `

Call:

prm eta_b5 0.5 min 0 max 1 prm s400_b5 0.05 min 0 max 1 prm s040_b5 0.05 min 0 max 1 prm s004_b5 0.05 min 0 max 1 prm s220_b5 0.05 min 0 max 1 prm s202_b5 0.05 min 0 max 1 prm s022_b5 0.05 min 0 max 1 prm s121_b5 0.05 min 0 max 1 prm s103_b5 0.05 min 0 max 1 prm s301_b5 0.05 min 0 max 1 TOF_Stephens_monoclinic(eta_b5, s400_b5, s040_b5, s004_b5, s220_b5, s202_b5, s022_b5, s121_b5, s103_b5, s301_b5)

Refined values:

prm eta_b5 0.0625561.686111_LIMIT_MIN_0 min 0 max 1
prm s400_b5 0.88125_4732.37809_LIMIT_MIN_0 min 0 max 1 prm s040_b5 0.881254789.89379_LIMIT_MIN_0 min 0 max 1
prm s004_b5 0.00625_562.368449_LIMIT_MIN_0 min 0 max 1 prm s220_b5 0.256253034.88555_LIMIT_MIN_0 min 0 max 1
prm s202_b5 0.88125_4752.97901_LIMIT_MIN_0 min 0 max 1 prm s022_b5 0.881255428.61298_LIMIT_MIN_0 min 0 max 1
prm s121_b5 0.88125_4917.87234_LIMIT_MIN_0 min 0 max 1 prm s103_b5 0.63125_3928.17146_LIMIT_MIN_0 min 0 max 1
prm s301_b5 0.88125_2386.28893_LIMIT_MIN_0 min 0 max 1

rowlesmr

Those large errors on the coefficients are concerning. They show that something not quite right is going on.

Other things to try are wrapping the gauss and lor equations in Abs(...) or 'Max(..., 0)`, but I can't really grep how those equations can be negative given all of the other constraints.

jmas5

rowlesmr
Updated macro:

Input 1:

Resultant output1 :

Result: Very large errors as before

Input 2:

Output 2:

Result: Quite large values (?), no large errors except for eta which is the only param with the min 0 max 1 left

Input 3:

Output 3:

Result: No large errors, I guess because they are now allowed to go to huge values which they seem to want to do. I don't know if this is physically meaningful.

I think, overall, the Abs(x) has eliminated the prior issues of negative FWHM, so that I don't need the min x max y restraints on the values. My question now then is, are the values it is refining to sensible? I have included the XRD stephens monoclinic parameters below for comparison, as they are quite different. My main concern is Eta, which I believed should be between 0 and 1? I don't know if this is different when working with TOF data?

Thanks again for your help!