1 | #if !defined(StickyHS_h) |
---|
2 | #define StickyHS_h |
---|
3 | |
---|
4 | /** |
---|
5 | * Structure definition for StickyHS_struct (struncture factor) parameters |
---|
6 | */ |
---|
7 | //[PYTHONCLASS] = StickyHSStructure |
---|
8 | //[DISP_PARAMS] = effect_radius |
---|
9 | //[DESCRIPTION] =<text> Structure Factor for interacting particles: . |
---|
10 | // |
---|
11 | // The interaction potential is |
---|
12 | // |
---|
13 | // U(r)= inf , r < 2R |
---|
14 | // = -Uo , 2R < r < 2R + w |
---|
15 | // = 0 , r >= 2R +w |
---|
16 | // |
---|
17 | // R: effective radius of the hardsphere |
---|
18 | // stickiness = [exp(Uo/kT)]/(12*perturb) |
---|
19 | // perturb = w/(w+ 2R) , 0.01 =< w <= 0.1 |
---|
20 | // w: The width of the square well ,w > 0 |
---|
21 | // v: The volume fraction , v > 0 |
---|
22 | // |
---|
23 | // Ref: Menon, S. V. G.,et.al., J. Chem. |
---|
24 | // Phys., 1991, 95(12), 9186-9190. |
---|
25 | // </text> |
---|
26 | //[FIXED]= effect_radius.width |
---|
27 | typedef struct { |
---|
28 | /// effect radius of hardsphere [A] |
---|
29 | // [DEFAULT]=effect_radius=50.0 [A] |
---|
30 | double effect_radius; |
---|
31 | |
---|
32 | /// Volume fraction |
---|
33 | // [DEFAULT]=volfraction= 0.1 |
---|
34 | double volfraction; |
---|
35 | |
---|
36 | /// Perturbation Parameter ( between 0.01 - 0.1) |
---|
37 | // [DEFAULT]=perturb=0.05 |
---|
38 | double perturb; |
---|
39 | |
---|
40 | /// Stickiness |
---|
41 | // [DEFAULT]=stickiness= 0.2 |
---|
42 | double stickiness; |
---|
43 | } StickyHSParameters; |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | /// 1D scattering function |
---|
48 | double StickyHS_analytical_1D(StickyHSParameters *pars, double q); |
---|
49 | |
---|
50 | /// 2D scattering function |
---|
51 | double StickyHS_analytical_2D(StickyHSParameters *pars, double q, double phi); |
---|
52 | double StickyHS_analytical_2DXY(StickyHSParameters *pars, double qx, double qy); |
---|
53 | |
---|
54 | #endif |
---|