Changeset 5f3c534 in sasmodels for sasmodels/models/hardsphere.py
- Timestamp:
- Mar 27, 2019 10:11:45 AM (5 years ago)
- Branches:
- master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- 9947865
- Parents:
- 055ec4f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/hardsphere.py
r0507e09 r5f3c534 1 1 # Note: model title and parameter table are inserted automatically 2 r"""Calculate the interparticle structure factor for monodisperse 2 r""" 3 Calculates the interparticle structure factor for monodisperse 3 4 spherical particles interacting through hard sphere (excluded volume) 4 interactions. 5 May be a reasonable approximation for other shapes of particles that 6 freely rotate, and for moderately polydisperse systems. Though strictly 7 the maths needs to be modified (no \Beta(Q) correction yet in sasview). 5 interactions. This $S(q)$ may also be a reasonable approximation for 6 other particle shapes that freely rotate (but see the note below), 7 and for moderately polydisperse systems. 8 9 .. note:: 10 11 This routine is intended for uncharged particles! For charged 12 particles try using the :ref:`hayter-msa` $S(q)$ instead. 13 14 .. note:: 15 16 Earlier versions of SasView did not incorporate the so-called 17 $\beta(q)$ ("beta") correction [1] for polydispersity and non-sphericity. 18 This is only available in SasView versions 4.2.2 and higher. 8 19 9 20 radius_effective is the effective hard sphere radius. 10 21 volfraction is the volume fraction occupied by the spheres. 11 22 12 In sasview the effective radius may be calculated from the parameters23 In SasView the effective radius may be calculated from the parameters 13 24 used in the form factor $P(q)$ that this $S(q)$ is combined with. 14 25 15 26 For numerical stability the computation uses a Taylor series expansion 16 at very small $qR$, there may be a very minor glitch at the transition point17 in some circumstances.27 at very small $qR$, but there may be a very minor glitch at the 28 transition point in some circumstances. 18 29 19 Th e S(Q) uses the Percus-Yevick closure where the interparticle20 potentialis30 This S(q) uses the Percus-Yevick closure relationship [2] where the 31 interparticle potential $U(r)$ is 21 32 22 33 .. math:: … … 27 38 \end{cases} 28 39 29 where $r$ is the distance from the center of thesphere of a radius $R$.40 where $r$ is the distance from the center of a sphere of a radius $R$. 30 41 31 42 For a 2D plot, the wave transfer is defined as … … 38 49 References 39 50 ---------- 51 52 .. [#] M Kotlarchyk & S-H Chen, *J. Chem. Phys.*, 79 (1983) 2461-2469 40 53 41 54 .. [#] J K Percus, J Yevick, *J. Phys. Rev.*, 110, (1958) 1 … … 63 76 [Hard sphere structure factor, with Percus-Yevick closure] 64 77 Interparticle S(Q) for random, non-interacting spheres. 65 May be a reasonable approximation for other shapes of66 particles that freely rotate, and for moderately polydisperse67 systems. Though strictly the maths needs to be modified -68 which sasview does not do yet.78 May be a reasonable approximation for other particle shapes 79 that freely rotate, and for moderately polydisperse systems 80 . The "beta(q)" correction is available in versions 4.2.2 81 and higher. 69 82 radius_effective is the hard sphere radius 70 83 volfraction is the volume fraction occupied by the spheres.
Note: See TracChangeset
for help on using the changeset viewer.