\(\renewcommand\AA{\unicode{x212B}}\)

Units

What are units?

Units are a set of small classes in Mantid that define a unit of measure, and the conversions between various units.

The Unit Factory is a Dynamic Factory that creates and hands out instances of Mantid Unit objects.

Available TOF Convertible units

The following units are available in the default Mantid distribution. These units are TOF convertible.

Name

ID (as known by Unit Factory)

Unit

Relevant equation

Time of flight

TOF

\(\mu s\)

TOF

Wavelength

Wavelength

\(\mathrm{\AA}\)

\(\lambda = \frac{h}{p} = \frac{h \times \mathrm{tof}}{m_N \times L_{tot}}\) (see below)

Energy

Energy

\(meV\)

\(E = \frac{1}{2} mv^2 = \frac{m_N}{2} \left ( \frac{L_{tot}}{\mathrm{tof}} \right )^2\)

Energy in wavenumber

Energy_inWavenumber

\(cm^{-1}\)

\(8.06554465 \times E\)

Momentum (k)

Momentum

\(\mathrm{\AA}^{-1}\)

\(k = \frac{2 \pi }{\lambda}=\frac{2 \pi \times m_N \times L_{tot}}{h \times \mathrm{tof}}\)

d-spacing

dSpacing

\(\mathrm{\AA}\)

\(TOF = DIFA \, d^2 + DIFC d + TZERO\) (see below)

Momentum transfer (Q)

MomentumTransfer

\(\mathrm{\AA}^{-1}\)

\(Q = 2 \, k \, sin \, \theta = \frac{4 \pi sin \theta}{\lambda}\)

Momentum transfer squared (\(Q^2\))

QSquared

\(\mathrm{\AA}^{-2}\)

\(Q^2 \frac{}{}\)

Energy transfer

DeltaE

\(meV\)

\(\Delta E = E_{i}-\frac{1}{2}m_N \left ( \frac{L_2}{\mathrm{tof}-L_1\sqrt{\frac{m_N}{2E_i}}} \right )^2\)

Energy transfer in wavenumber

DeltaE_inWavenumber

\(cm^{-1}\)

\(8.06554465 \times \Delta E\)

Spin Echo Length

SpinEchoLength

\(nm\)

\(constant \times \lambda^2\)
The constant is supplied in eFixed

Spin Echo Time

SpinEchoTime

\(ns\)

\(constant \times \lambda^3\)
The constant is supplied in eFixed

d-spacingPerpendicular

dSpacingPerpendicular

\(\mathrm{\AA}\)

\(d_{\perp} = \sqrt{\lambda^2 - 2\log\cos\theta}\)

Where \(L_1\) and \(L_2\) are sample to the source and sample to detector distances respectively, \(L_{tot} = L_1+L_2\) and \(E_i\) is the energy of neutrons leaving the source. \(\theta\) here is the Bragg scattering angle (e.g. half of the \(\theta\)-angle used in spherical coordinate system directed along Mantid z-axis)

Note on Wavelength: If the emode property in ConvertUnits is specified as inelastic Direct/Indirect (inelastic) then the conversion to wavelength will take into account the fixed initial/final energy respectively. Units conversion into elastic momentum transfer (MomentumTransfer) will throw in elastic mode (emode=0) on inelastic workspace (when energy transfer is specified along x-axis)

Note on d-spacing: The coefficients DIFA, DIFC and TZERO may be obtained via calibration of a TOF diffraction instrument. In the absence of a calibration, DIFA=TZERO=0 and the default value of DIFC is:

\(DIFC = 10^{-4} \frac{m_N}{h} (L_1 + L_2) 2 \sin(\theta)\)

where the scaling factor adjusts for the fact that DIFC is required in units of \(\mu s\) per \(\mathrm{\AA}\).

d-spacingPerpendicular is a unit invented in J. Appl. Cryst. (2015) 48, pp. 1627–1636 for 2D Rietveld refinement of angular and wavelength-dispersive neutron time-of-flight powder diffraction data. Together with the d-Spacing \(d\), d-SpacingPerpendicular \(d_{\perp}\) forms a new orthogonal coordinate system.

Available non-TOF Convertible units

The following units are available in the default Mantid distribution. These units cannot be converted to or from TOF.

Name

ID (as known by Unit Factory)

Unit

Description

Empty

No unit

An empty label

t

Time

\(s\)

An independent unit of time not related to energy or TOF

Scattering angle

Degrees

\(degrees\)

Degrees is a measurement of angular position

Temperature

Temperature

\(K\)

Temperature in Kelvin

Working with Units in Python

Accessing units on workspaces

Units on MatrixWorkspaces are accessed via the Axis.

ws = CreateSampleWorkspace()
for i in range(ws.axes()):
    axis = ws.getAxis(i)
    print("Axis {0} is a {1}{2}{3}".format(i,
                                           "Spectrum Axis" if axis.isSpectra() else "",
                                           "Text Axis" if axis.isText() else "",
                                           "Numeric Axis" if axis.isNumeric() else ""))

    unit = axis.getUnit()
    print("\t caption:{0}".format(unit.caption()))
    print("\t symbol:{0}".format(unit.symbol()))

Output:

Axis 0 is a Numeric Axis
   caption:Time-of-flight
   symbol:microsecond
Axis 1 is a Spectrum Axis
   caption:Spectrum
   symbol:

Setting the axisLabel to a Label of your choice

ws = CreateSampleWorkspace()
axis = ws.getAxis(1)
# Create a new axis
axis.setUnit("Label").setLabel('Temperature', 'K')

unit = axis.getUnit()
print("New caption:{0}".format(unit.caption()))
print("New symbol:{0}".format(unit.symbol()))

Output:

New caption:Temperature
New symbol:K

Adding new units

Writing and adding a new unit is relatively straightforward. Instructions will appear here in due course. In the meantime if a unit that you require is missing, then please contact the development team and we will add it to the default Mantid library.

Category: Concepts