2D TDS intensity distribution without interpolation:ΒΆ

The script make_TDS_intensityPlaneClone calculates x-ray and neutron TDS intensities using only the eigenvectors and eigenvalues provided py the input file. The following steps must have been performed beforehand:

  • Symmetrysation
  • Fourier Interpolation
  • Debye-Waller factors for the temperature of choice
  • The usage is

    make_TDS_IntensityPlaneClone QW_castep_filename input_filename castep_filename 

    the file castep_filename can be either the name of the original castep output our the associated hdf5 file. In any case the associated hdf5 file must exist already.The file QW_castep_filename must contain eigenvectors and eigenvalues for the q-values of interest. A list of q-values can be created by the script make_TDS_IntensityPlane, eigenvectors and eigenvalues on these q-points must be calculated by an external program (e.g. CASTEP)

  • The input_file must set the variables:

    • Mandatories
      • Nfour_interp
      • Temperature
      • Saturation
      • lowerLimit
      • bottom_meV
      • Nqs
      • DQ
    • Definition of Plane through 3 points * redA, redB, redC * redCenter (optional)
    • Definition of plane via normal and point * redNormal * redX_idea * redCenter
    • If NEUTRONCALC==0 ( default )
      • Lambda
      • CohB
      • NeutronE
    • Optional
      • branchWeight
      • energyWindow
      • Eigscal
      • UniqueIon
      • LINEAR
      • COLOR
  • Interactions

    • right mouse click on graph plots intensity profile and displays position in reciprocal space

The input variables are documented with docstrings below

make_TDS_IntensityPlaneClone.APPLYTIMEREVERSAL = 1

this is one by default. Must be coherent with previous steps

make_TDS_IntensityPlaneClone.Nfour_interp = 4

The number of points in each direction of the 3D reciprocal grid. The hdf5 file must contain a previous pretreatement done with the same parameter.

make_TDS_IntensityPlaneClone.Temperature = 100
The hdf5 file must include DW factors calculated ALSO at this temperature.
The temperature at which DW factors have been calculated.
Units are Kelvin
make_TDS_IntensityPlaneClone.Lambda = 0.2

For X-ray scattering: the wavelenght in Angstroems.

make_TDS_IntensityPlaneClone.NEUTRONCALC = 0

This activate neutron scattering calculations.

make_TDS_IntensityPlaneClone.CohB = None

When neutroncalc is on, this must be a dictionary: for each atom name the Coherent lenght.

make_TDS_IntensityPlaneClone.NeutronE = 0

When neutroncalc is on, the neutron energy in meV

make_TDS_IntensityPlaneClone.lowerLimit = 1e-06

To correct the dynamical range, when displaying, the intensity is clipped to prevent it going below this value

make_TDS_IntensityPlaneClone.bottom_meV = 0

Calculated eigenvalues are clipped to this value before use.

make_TDS_IntensityPlaneClone.Nqs = 100

Number of q-points along one axis

make_TDS_IntensityPlaneClone.DQ = 0.08

q-resolution in 1/Ang

make_TDS_IntensityPlaneClone.redA = [1.0, 1.0, 0.0]
  1. q-point defining plane
make_TDS_IntensityPlaneClone.redB = [0.0, 0.0, 0.0]
  1. q-point defining plane
make_TDS_IntensityPlaneClone.redC = [0.0, 0.0, 1.0]
  1. q-point defining plane
make_TDS_IntensityPlaneClone.redCenter = [0.0, 0.0, 0.0]

center of plane

make_TDS_IntensityPlaneClone.redNormal = [0.0, 0.0, 1.0]

Normal to plane

make_TDS_IntensityPlaneClone.redX_idea = [1.0, -1.0, 0.0]

Orientation of plane

make_TDS_IntensityPlaneClone.branchWeight = None

a list of weights : one per branch.

make_TDS_IntensityPlaneClone.energyWindow = None

a list of two numbers : minimum and maximum in meV

make_TDS_IntensityPlaneClone.Eigscal = 0

EigScal==1 Intensity are calculated WITH only eigenvector scalar products ; ==2 Intensity are calculated WITH only eigenvector scalar products PLUS Mass factor

make_TDS_IntensityPlaneClone.UniqueIon = -1

If >=0. Selects one ion. All the other will be silent.

make_TDS_IntensityPlaneClone.RemoveBose = 0

RemoveBose=1 desactivates the Bose statistics.

make_TDS_IntensityPlaneClone.LINEAR = 1

LINEAR = 1 plotting intensity in linear scale, linear = 0 in log scale

make_TDS_IntensityPlaneClone.COLOR = 1

Intensity in color or grayscale