Home > Scripts > Models > Specialized > sf_sticky_hard_spheres.m



y = sf_sticky_hard_spheres(p, x, [y]) : Sticky Hard Sphere structure factor [Baxter/Menon]


function y=sf_sticky_hard_spheres(varargin)


 y = sf_sticky_hard_spheres(p, x, [y]) : Sticky Hard Sphere structure factor [Baxter/Menon]

   iFunc/sf_sticky_hard_spheres Sticky Hard Sphere structure factor 
          (for e.g. spheres in dielectric liquids)
     The 'x' wave-vector/momentum axis is usually in nm-1 or Angs-1.
     Typical values for parameters are R=50 Angs, rho=0.04, tau=0.15
     The Sticky Hard Sphere model converges to the Hard Sphere model [PY] with
       increasing tau (e.g. tau > 10).
     The model returns the S(q) structure factor.

     See: S.V.G. Menon, C. Manobar, and K. Srinivasa Rao. J. Chem. Phys., 95 (1991) 9186
          Extracted from sasfit/sasfit_sq/sasfit_sq_StickyHardSphere.c
     I. Bressler, et al, Journal of Applied Crystallography, 2015, 48 (5), 1587-1598

 input:  p: sticky hard sphere model parameters (double)
            p = [ R=hard_sphere_radius rho=Volume_Fraction tau=stickiness ]
          or 'guess'
         x: wave-vector/momentum axis (double, e.g. nm-1 or Angs-1)
         y: when values are given and p='guess', a guess of the parameters is performed (double)
 output: y: model value
 ex:     y=sf_sticky_hard_spheres([10 0.2 0.15], 0:0.01:1); or plot(sf_sticky_hard_spheres,[10 0.2 0.15], 0:0.01:1)

 Version: Nov. 26, 2018
 See also iFunc, iFunc/fits, iFunc/plot
 (c) E.Farhi, ILL. License: EUPL.


This function calls: This function is called by:
Generated on Mon 26-Nov-2018 15:08:42 by m2html © 2005. iFit (c) E.Farhi/ILL EUPL 1.1