Sqw_Bosify: apply the 'Bose' factor (detailed balance) to a classical data set.
The initial data set should obey S(q,w) = S(q,-w), i.e. be 'classical'.
The S(q,w) is a dynamic structure factor aka scattering function.
Sqw_Bosify(Sqw, 0) makes a check of the input data set and returns.
input:
s: Sqw data set (classical, symmetric in energy, no T Bose factor)
e.g. 2D data set with w as 1st axis (rows, meV), q as 2nd axis (Angs-1).
T: when given, Temperature to use for Bose. When not given, the Temperature
is searched in the object. The temperature is in [K]. 1 meV=11.605 K.
type: 'Schofield' or 'harmonic' or 'standard' (default)
conventions:
omega = Ei-Ef = energy lost by the neutron, given in [meV]
omega > 0, neutron looses energy, can not be higher than Ei (Stokes)
omega < 0, neutron gains energy, anti-Stokes
Egelstaff, Eq (9.25) p189
S(q,-w) = exp(-hw/kT) S(q,w)
S(q,w) = exp( hw/kT) S(q,-w)
S(q,w) = Q(w) S*(q,w) with S*=classical limit
for omega > 0, S(q,w) > S(q,-w)
The semi-classical correction, Q, aka 'quantum' correction factor,
can be selected from the optional 'type' argument:
Q = exp(hw_kT/2) 'Schofield' or 'Boltzmann'
Q = hw_kT./(1-exp(-hw_kT)) 'harmonic' or 'Bader'
Q = 2./(1+exp(-hw_kT)) 'standard' or 'Frommhold' (default)
The 'Boltzmann' correction leads to a divergence of the S(q,w) for e.g. w above
few 100 meV. The 'harmonic' correction provides a reasonable correction but does
not fully avoid the divergence at large energies.
Bose factor: n(w) = 1./(exp(w*11.605/T) -1) ~ exp(-w*11.605/T)
w in [meV], T in [K]
References:
B. Hehr, http://www.lib.ncsu.edu/resolver/1840.16/7422 PhD manuscript (2010).
S. A. Egorov, K. F. Everitt and J. L. Skinner. J. Phys. Chem., 103, 9494 (1999).
P. Schofield. Phys. Rev. Lett., 4, 239 (1960).
J. S. Bader and B. J. Berne. J. Chem. Phys., 100, 8359 (1994).
T. D. Hone and G. A. Voth. J. Chem. Phys., 121, 6412 (2004).
L. Frommhold. Collision-induced absorption in gases, 1 st ed., Cambridge
Monographs on Atomic, Molecular, and Chemical Physics, Vol. 2,
Cambridge Univ. Press: London (1993).
Example: s = Sqw_Bosify(s, 300);
See also: Sqw_deBosify, Sqw_symmetrize, Sqw_dynamic_range, Sqw_scatt_xs
(c) E.Farhi, ILL. License: EUPL.