Changeset eb69cce in sasmodels for sasmodels/models/lamellarCaille.py
- Timestamp:
- Nov 30, 2015 9:18:41 PM (8 years ago)
- Branches:
- master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- d18f8a8
- Parents:
- d138d43
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/lamellarCaille.py
r3e428ec reb69cce 3 3 This model provides the scattering intensity, $I(q) = P(q) S(q)$, for a 4 4 lamellar phase where a random distribution in solution are assumed. 5 Here a Caille $S(Q)$ is used for the lamellar stacks. 5 Here a Caille $S(q)$ is used for the lamellar stacks. 6 7 Definition 8 ---------- 6 9 7 10 The scattering intensity $I(q)$ is 8 11 9 .. math: 12 .. math:: 10 13 11 14 I(q) = 2\pi \frac{P(q)S(q)}{\delta q^2} … … 13 16 The form factor is 14 17 15 .. math: 18 .. math:: 16 19 17 20 P(q) = \frac{2\Delta\rho^2}{q^2}\left(1-\cos q\delta \right) … … 19 22 and the structure factor is 20 23 21 .. math: 24 .. math:: 22 25 23 26 S(q) = 1 + 2 \sum_1^{N-1}\left(1-\frac{n}{N}\right) … … 26 29 where 27 30 28 .. math: 31 .. math:: 29 32 30 \begin{eqnarray} 31 \alpha(n) &=& \frac{\eta_{cp}}{4\pi^2} \left(\ln(\pi n)+\gamma_E\right) \\ 32 \gamma_E &=& 0.5772156649 && \text{Euler's constant} \\ 33 \eta_{cp} &=& \frac{q_o^2k_B T}{8\pi\sqrt{K\overline{B}}} && \text{Caille constant} 34 \end{eqnarray} 33 %\begin{align*} % isn't working with pdflatex 34 \begin{array}{rll} 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 & 42 \end{array} 35 43 36 44 Here $d$ = (repeat) spacing, $\delta$ = bilayer thickness, … … 54 62 q = \sqrt{q_x^2 + q_y^2} 55 63 56 The returned value is in units of |cm^-1|, on absolute scale. 64 .. figure:: img/lamellarCaille_1d.jpg 57 65 58 .. image:: img/lamellarCaille_1d.jpg 66 1D plot using the default values (w/6000 data point). 59 67 60 *Figure. 1D plot using the default values (w/6000 data point).* 61 62 Our model uses the form factor calculations as implemented in a c library 63 provided by the NIST Center for Neutron Research (Kline, 2006). 64 65 REFERENCE 66 --------- 68 References 69 ---------- 67 70 68 71 F Nallet, R Laversanne, and D Roux, J. Phys. II France, 3, (1993) 487-502
Note: See TracChangeset
for help on using the changeset viewer.