source: sasmodels/sasmodels/models/lamellarCailleHG.py @ d18f8a8

core_shell_microgelscostrafo411magnetic_modelrelease_v0.94release_v0.95ticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change on this file since d18f8a8 was d18f8a8, checked in by Paul Kienzle <pkienzle@…>, 8 years ago

fix multiline equation alignment

  • Property mode set to 100644
File size: 4.6 KB
Line 
1# Note: model title and parameter table are inserted automatically
2r"""
3This model provides the scattering intensity, $I(q) = P(q)S(q)$, for a lamellar
4phase where a random distribution in solution are assumed. Here a Caille $S(q)$
5is used for the lamellar stacks.
6
7The scattering intensity $I(q)$ is
8
9.. math::
10
11    I(q) = 2 \pi \frac{P(q)S(q)}{\delta q^2}
12
13
14The form factor $P(q)$ is
15
16.. math::
17
18        P(q) = \frac{4}{q^2}\big\{
19        \Delta\rho_H \left[\sin[q(\delta_H + \delta_T)] - \sin(q\delta_T)\right]
20            + \Delta\rho_T\sin(q\delta_T)\big\}^2
21
22and the structure factor $S(q)$ is
23
24.. math::
25
26    S(q) = 1 + 2 \sum_1^{N-1}\left(1-\frac{n}{N}\right)
27        \cos(qdn)\exp\left(-\frac{2q^2d^2\alpha(n)}{2}\right)
28
29where
30
31.. math::
32    :nowrap:
33
34    \begin{align*}
35    \alpha(n) &= \frac{\eta_{cp}}{4\pi^2} \left(\ln(\pi n)+\gamma_E\right)
36              &&  \\
37    \gamma_E  &= 0.5772156649
38              && \text{Euler's constant} \\
39    \eta_{cp} &= \frac{q_o^2k_B T}{8\pi\sqrt{K\overline{B}}}
40              && \text{Caille constant}
41    \end{align*}
42
43
44$\delta_T$ is the tail length (or *tail_length*), $\delta_H$ is the head
45thickness (or *head_length*), $\Delta\rho_H$ is SLD(headgroup) - SLD(solvent),
46and $\Delta\rho_T$ is SLD(tail) - SLD(headgroup). Here $d$ is (repeat) spacing,
47$K$ is smectic bending elasticity, $B$ is compression modulus, and $N$ is the
48number of lamellar plates (*Nlayers*).
49
50NB: **When the Caille parameter is greater than approximately 0.8 to 1.0, the
51assumptions of the model are incorrect.**  And due to a complication of the
52model function, users are responsible for making sure that all the assumptions
53are handled accurately (see the original reference below for more details).
54
55Non-integer numbers of stacks are calculated as a linear combination of
56results for the next lower and higher values.
57
58The 2D scattering intensity is calculated in the same way as 1D, where
59the $q$ vector is defined as
60
61.. math::
62
63    q = \sqrt{q_x^2 + q_y^2}
64
65.. figure:: img/lamellarCailleHG_1d.jpg
66
67    1D plot using the default values (w/6000 data point).
68
69References
70----------
71
72F Nallet, R Laversanne, and D Roux, J. Phys. II France, 3, (1993) 487-502
73
74also in J. Phys. Chem. B, 105, (2001) 11081-11088
75"""
76from numpy import inf
77
78name = "lamellarCailleHG"
79title = "Random lamellar sheet with Caille structure factor"
80description = """\
81    [Random lamellar phase with Caille  structure factor]
82        randomly oriented stacks of infinite sheets
83        with Caille S(Q), having polydisperse spacing.
84        layer thickness =(H+T+T+H) = 2(Head+Tail)
85        sld = Tail scattering length density
86        sld_head = Head scattering length density
87        sld_solvent = solvent scattering length density
88        background = incoherent background
89        scale = scale factor
90"""
91category = "shape:lamellae"
92
93parameters = [
94              #   [ "name", "units", default, [lower, upper], "type",
95              #     "description" ],
96              [ "tail_length", "Ang",  10, [0, inf], "volume",
97                "Tail thickness" ],
98              [ "head_length", "Ang",  2, [0, inf], "volume",
99                "head thickness" ],
100              [ "Nlayers", "",  30, [0, inf], "",
101                "Number of layers" ],
102              [ "spacing", "Ang", 40., [0.0,inf], "volume",
103                "d-spacing of Caille S(Q)" ],
104              [ "Caille_parameter", "", 0.001, [0.0,0.8], "",
105                "Caille parameter" ],
106              [ "sld", "1e-6/Ang^2", 0.4, [-inf,inf], "",
107                "Tail scattering length density" ],
108              [ "head_sld", "1e-6/Ang^2", 2.0, [-inf,inf], "",
109                "Head scattering length density" ],
110              [ "solvent_sld", "1e-6/Ang^2", 6, [-inf,inf], "",
111                "Solvent scattering length density" ],
112    ]
113
114source = [ "lamellarCailleHG_kernel.c"]
115
116# No volume normalization despite having a volume parameter
117# This should perhaps be volume normalized?
118form_volume = """
119    return 1.0;
120    """
121
122Iqxy = """
123    return Iq(sqrt(qx*qx+qy*qy), IQ_PARAMETERS);
124    """
125
126# ER defaults to 0.0
127# VR defaults to 1.0
128
129demo = dict(
130            scale=1, background=0,
131            Nlayers=20,
132            spacing=200., Caille_parameter=0.05,
133            tail_length=15,head_length=10,
134            #sld=-1, head_sld=4.0, solvent_sld=6.0,
135            sld=-1, head_sld=4.1, solvent_sld=6.0,
136            tail_length_pd= 0.1, tail_length_pd_n=20,
137            head_length_pd= 0.05, head_length_pd_n=30,
138            spacing_pd= 0.2, spacing_pd_n=40
139           )
140
141oldname = 'LamellarPSHGModel'
142oldpars = dict(tail_length='deltaT',head_length='deltaH',Nlayers='n_plates',Caille_parameter='caille', sld='sld_tail', head_sld='sld_head',solvent_sld='sld_solvent')
Note: See TracBrowser for help on using the repository browser.