Abins dialog.
Table of Contents
| Name | Direction | Type | Default | Description |
|---|---|---|---|---|
| DFTprogram | Input | string | CASTEP | DFT program which was used for a phonon calculation. Allowed values: [‘CASTEP’, ‘CRYSTAL’, ‘DMOL3’, ‘GAUSSIAN’] |
| PhononFile | Input | string | Mandatory | File with the data from a phonon calculation. Allowed values: [‘phonon’, ‘out’, ‘outmol’, ‘log’] |
| ExperimentalFile | Input | string | File with the experimental inelastic spectrum to compare. Allowed values: [‘raw’, ‘dat’] | |
| Temperature | Input | number | 10 | Temperature in K for which dynamical structure factor S should be calculated. |
| Scale | Input | number | 1 | Scale the intensity by the given factor. Default is no scaling. |
| SampleForm | Input | string | Powder | Form of the sample: Powder. Allowed values: [‘Powder’] |
| Instrument | Input | string | TOSCA | Name of an instrument for which analysis should be performed. Allowed values: [‘TOSCA’] |
| Atoms | Input | str list | List of atoms to use to calculate partial S.If left blank, workspaces with S for all types of atoms will be calculated. | |
| SumContributions | Input | boolean | False | Sum the partial dynamical structure factors into a single workspace. |
| ScaleByCrossSection | Input | string | Incoherent | Scale the partial dynamical structure factors by the scattering cross section. Allowed values: [‘Total’, ‘Incoherent’, ‘Coherent’] |
| QuantumOrderEventsNumber | Input | string | 1 | Number of quantum order effects included in the calculation (1 -> FUNDAMENTALS, 2-> first overtone + FUNDAMENTALS + 2nd order combinations, 3-> FUNDAMENTALS + first overtone + second overtone + 2nd order combinations + 3rd order combinations etc...). Allowed values: [‘1’, ‘2’, ‘3’, ‘4’] |
| OutputWorkspace | Output | Workspace | Mandatory | Name to give the output workspace. |
Abins is a plugin for Mantid which allows scientists to compare experimental and theoretical inelastic neutron scattering spectra (INS).
Abins requires a file with the ab-initio phonon data to perform INS analysis. Currently output data from CASTEP (.phonon) and CRYSTAL (.out) DFT programs can be used to perform analysis. Optionally a user can provide an experimental file with measured dynamical structure factor S in order to directly compare theoretical and experimental spectra. A user can produce one dimensional INS spectrum which can be compared against TOSCA and TOSCA-like instruments.
After successfully performed analysis a user obtains a group workspace which stores theoretical and optionally experimental spectra. Currently a user can produce theoretical spectrum for the given type of atom (for example for benzene two types of atoms: C, H) and for each quantum event (up to fourth order). A user can produce a total theoretical spectrum for the given atom which is a sum over all considered quantum events for that atom. A user can also produce a total spectrum for the whole considered system. Dynamical structure factor S is calculated for all atoms in the system. If needed a user can also include in a simulation elevated temperature.
A description about the implemented working equations can be found here.
Abins is in constant development and suggestions how to make it better are very welcome. For any suggestions of enhancements please contact Krzysztof Dymkowski (krzysztof.dymkowski@stfc.ac.uk).
Note
To run these usage examples please first download the usage data, and add these to your path. In MantidPlot this is done using Manage User Directories.
Example - loading CASTEP phonon data :
benzene_wrk = Abins(DFTProgram="CASTEP", PhononFile="benzene.phonon", QuantumOrderEventsNumber="1")
for name in benzene_wrk.getNames():
print name
Output:
benzene_wrk_C_total
benzene_wrk_C
benzene_wrk_H_total
benzene_wrk_H
Example - loading CRYSTAL phonon data :
wrk=Abins(DFTProgram="CRYSTAL", PhononFile="b3lyp.out", QuantumOrderEventsNumber="1")
for name in wrk.getNames():
print name
Output:
wrk_C_total
wrk_C
wrk_H_total
wrk_H
wrk_N_total
wrk_N
wrk_Na_total
wrk_Na
wrk_O_total
wrk_O
Example - calling Abins with more arguments:
wrk_verbose=Abins(DFTprogram="CASTEP", PhononFile="benzene.phonon", ExperimentalFile="benzene_experimental.dat",
Temperature=10, SampleForm="Powder", Instrument="TOSCA", Atoms="H", SumContributions=True,
QuantumOrderEventsNumber="1", ScaleByCrossSection="Incoherent")
for name in wrk_verbose.getNames():
print name
Output:
experimental_wrk
wrk_verbose_total
wrk_verbose_H_total
wrk_verbose_H
Categories: Algorithms | Simulation