source: sasmodels/sasmodels/models/hollow_rectangular_prism_thin_walls.py @ cf3d0ce

core_shell_microgelsmagnetic_modelticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change on this file since cf3d0ce was 99658f6, checked in by grethevj, 5 years ago

updated ER functions including cylinder excluded volume, to match 4.x

  • Property mode set to 100644
File size: 4.7 KB
RevLine 
[deb7ee0]1# rectangular_prism model
2# Note: model title and parameter table are inserted automatically
3r"""
4Definition
5----------
6
[6e7d7b6]7
8This model provides the form factor, $P(q)$, for a hollow rectangular
9prism with infinitely thin walls. It computes only the 1D scattering, not the 2D.
[deb7ee0]10The 1D scattering intensity for this model is calculated according to the
[6e7d7b6]11equations given by Nayuk and Huber\ [#Nayuk2012]_.
[deb7ee0]12
13Assuming a hollow parallelepiped with infinitely thin walls, edge lengths
[ab2aea8]14$A \le B \le C$ and presenting an orientation with respect to the
15scattering vector given by $\theta$ and $\phi$, where $\theta$ is the angle
16between the $z$ axis and the longest axis of the parallelepiped $C$, and
17$\phi$ is the angle between the scattering vector (lying in the $xy$ plane)
18and the $y$ axis, the form factor is given by
[deb7ee0]19
20.. math::
21
[ab2aea8]22    P(q) = \frac{1}{V^2} \frac{2}{\pi} \int_0^{\frac{\pi}{2}}
23           \int_0^{\frac{\pi}{2}} [A_L(q)+A_T(q)]^2 \sin\theta\,d\theta\,d\phi
[deb7ee0]24
[ab2aea8]25where
[deb7ee0]26
27.. math::
28
[ab2aea8]29    V &= 2AB + 2AC + 2BC \\
30    A_L(q) &=  8 \times \frac{
31            \sin \left( \tfrac{1}{2} q A \sin\phi \sin\theta \right)
32            \sin \left( \tfrac{1}{2} q B \cos\phi \sin\theta \right)
33            \cos \left( \tfrac{1}{2} q C \cos\theta \right)
34        }{q^2 \, \sin^2\theta \, \sin\phi \cos\phi} \\
35    A_T(q) &=  A_F(q) \times
36      \frac{2\,\sin \left( \tfrac{1}{2} q C \cos\theta \right)}{q\,\cos\theta}
[deb7ee0]37
38and
39
40.. math::
[ab2aea8]41
42  A_F(q) =  4 \frac{ \cos \left( \tfrac{1}{2} q A \sin\phi \sin\theta \right)
43                       \sin \left( \tfrac{1}{2} q B \cos\phi \sin\theta \right) }
[deb7ee0]44                     {q \, \cos\phi \, \sin\theta} +
[ab2aea8]45              4 \frac{ \sin \left( \tfrac{1}{2} q A \sin\phi \sin\theta \right)
46                       \cos \left( \tfrac{1}{2} q B \cos\phi \sin\theta \right) }
[deb7ee0]47                     {q \, \sin\phi \, \sin\theta}
48
49The 1D scattering intensity is then calculated as
50
51.. math::
52
[ab2aea8]53  I(q) = \text{scale} \times V \times (\rho_\text{p} - \rho_\text{solvent})^2 \times P(q)
54
[6e7d7b6]55where $V$ is the surface area of the rectangular prism, $\rho_\text{p}$
56is the scattering length density of the parallelepiped, $\rho_\text{solvent}$
57is the scattering length density of the solvent, and (if the data are in
58absolute units) *scale* is related to the total surface area.
[deb7ee0]59
60**The 2D scattering intensity is not computed by this model.**
61
62
63Validation
64----------
65
66Validation of the code was conducted  by qualitatively comparing the output
[6e7d7b6]67of the 1D model to the curves shown in (Nayuk, 2012\ [#Nayuk2012]_).
[deb7ee0]68
[aa2edb2]69
70References
71----------
[deb7ee0]72
[6e7d7b6]73.. [#Nayuk2012] R Nayuk and K Huber, *Z. Phys. Chem.*, 226 (2012) 837-854
[99658f6]74L. Onsager, Ann. New York Acad. Sci. 51, 627-659 (1949).
[6e7d7b6]75
76
77Authorship and Verification
78----------------------------
79
80* **Author:** Miguel Gonzales **Date:** February 26, 2016
81* **Last Modified by:** Paul Kienzle **Date:** October 15, 2016
82* **Last Reviewed by:** Paul Butler **Date:** September 07, 2018
[deb7ee0]83"""
84
[2d81cfe]85import numpy as np
[deb7ee0]86from numpy import pi, inf, sqrt
87
[3d8283b]88name = "hollow_rectangular_prism_thin_walls"
89title = "Hollow rectangular parallelepiped with thin walls."
[deb7ee0]90description = """
[3d8283b]91    I(q)= scale*V*(sld - sld_solvent)^2*P(q)+background
[deb7ee0]92        with P(q) being the form factor corresponding to a hollow rectangular
93        parallelepiped with infinitely thin walls.
94"""
95category = "shape:parallelepiped"
96
97#             ["name", "units", default, [lower, upper], "type","description"],
[42356c8]98parameters = [["sld", "1e-6/Ang^2", 6.3, [-inf, inf], "sld",
[deb7ee0]99               "Parallelepiped scattering length density"],
[42356c8]100              ["sld_solvent", "1e-6/Ang^2", 1, [-inf, inf], "sld",
[deb7ee0]101               "Solvent scattering length density"],
[a807206]102              ["length_a", "Ang", 35, [0, inf], "volume",
[deb7ee0]103               "Shorter side of the parallelepiped"],
104              ["b2a_ratio", "Ang", 1, [0, inf], "volume",
105               "Ratio sides b/a"],
106              ["c2a_ratio", "Ang", 1, [0, inf], "volume",
107               "Ratio sides c/a"],
108             ]
109
[3d8283b]110source = ["lib/gauss76.c", "hollow_rectangular_prism_thin_walls.c"]
[71b751d]111have_Fq = True
[ee60aa7]112effective_radius_type = [
[99658f6]113    "equivalent cylinder excluded volume", "equivalent outer volume sphere", 
114    "half length_a", "half length_b", "half length_c",
[ee60aa7]115    "equivalent outer circular cross-section",
116    "half ab diagonal", "half diagonal",
117    ]
[deb7ee0]118
119
[31df0c9]120def random():
121    a, b, c = 10**np.random.uniform(1, 4.7, size=3)
122    pars = dict(
123        length_a=a,
124        b2a_ratio=b/a,
125        c2a_ratio=c/a,
126    )
127    return pars
128
129
[deb7ee0]130# parameters for demo
131demo = dict(scale=1, background=0,
[ab2aea8]132            sld=6.3, sld_solvent=1.0,
[a807206]133            length_a=35, b2a_ratio=1, c2a_ratio=1,
134            length_a_pd=0.1, length_a_pd_n=10,
[deb7ee0]135            b2a_ratio_pd=0.1, b2a_ratio_pd_n=1,
136            c2a_ratio_pd=0.1, c2a_ratio_pd_n=1)
137
[6dd90c1]138tests = [[{}, 0.2, 0.837719188592],
139         [{}, [0.2], [0.837719188592]],
[deb7ee0]140        ]
Note: See TracBrowser for help on using the repository browser.