Changeset 5f3c534 in sasmodels for sasmodels/models/squarewell.py
- Timestamp:
- Mar 27, 2019 12:11:45 PM (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/squarewell.py
r0507e09 r5f3c534 1 1 # Note: model title and parameter table are inserted automatically 2 2 r""" 3 This calculates the interparticle structure factor for a square well fluid 4 spherical particles. The mean spherical approximation (MSA) closure was 5 used for this calculation, and is not the most appropriate closure for 6 a n attractive interparticle potential. This solution has been compared7 t o Monte Carlo simulations for a square well fluid, showing this calculation8 to be limited in applicability to well depths $\epsilon < 1.5$ kT and 9 volume fractions $\phi < 0.08$.3 Calculates the interparticle structure factor for a hard sphere fluid 4 with a narrow, attractive, square well potential. **The Mean Spherical 5 Approximation (MSA) closure relationship is used, but it is not the most 6 appropriate closure for an attractive interparticle potential.** However, 7 the solution has been compared to Monte Carlo simulations for a square 8 well fluid and these show the MSA calculation to be limited to well 9 depths $\epsilon < 1.5$ kT and volume fractions $\phi < 0.08$. 10 10 11 11 Positive well depths correspond to an attractive potential well. Negative 12 12 well depths correspond to a potential "shoulder", which may or may not be 13 physically reasonable. The stickyhardsphere model may be a better choice in 14 some circumstances. Computed values may behave badly at extremely small $qR$. 13 physically reasonable. The :ref:`stickyhardsphere` model may be a better 14 choice in some circumstances. 15 16 Computed values may behave badly at extremely small $qR$. 17 18 .. note:: 19 20 Earlier versions of SasView did not incorporate the so-called 21 $\beta(q)$ ("beta") correction [2] for polydispersity and non-sphericity. 22 This is only available in SasView versions 4.2.2 and higher. 15 23 16 24 The well width $(\lambda)$ is defined as multiples of the particle diameter … … 18 26 19 27 The interaction potential is: 20 21 .. image:: img/squarewell.png22 28 23 29 .. math:: … … 29 35 \end{cases} 30 36 31 where $r$ is the distance from the center of thesphere of a radius $R$.37 where $r$ is the distance from the center of a sphere of a radius $R$. 32 38 33 In sasview the effective radius may be calculated from the parameters39 In SasView the effective radius may be calculated from the parameters 34 40 used in the form factor $P(q)$ that this $S(q)$ is combined with. 35 41 … … 46 52 .. [#] R V Sharma, K C Sharma, *Physica*, 89A (1977) 213 47 53 54 .. [#] M Kotlarchyk and S-H Chen, *J. Chem. Phys.*, 79 (1983) 2461-2469 55 48 56 Source 49 57 ------ … … 56 64 * **Author:** 57 65 * **Last Modified by:** 58 * **Last Reviewed by:** 66 * **Last Reviewed by:** Steve King **Date:** March 27, 2019 59 67 * **Source added by :** Steve King **Date:** March 25, 2019 60 68 """ … … 64 72 65 73 name = "squarewell" 66 title = "Square well structure factor , with MSAclosure"74 title = "Square well structure factor with Mean Spherical Approximation closure" 67 75 description = """\ 68 76 [Square well structure factor, with MSA closure] 69 Interparticle structure factor S(Q) for a hard sphere fluid with70 a narrow attractive well. Fits are prone to deliver non-physical71 parameters, use with care and read the references in the full manual.72 In sasview the effective radius will be calculated from the73 parameters used in P(Q).77 Interparticle structure factor S(Q) for a hard sphere fluid 78 with a narrow attractive well. Fits are prone to deliver non- 79 physical parameters; use with care and read the references in 80 the model documentation.The "beta(q)" correction is available 81 in versions 4.2.2 and higher. 74 82 """ 75 83 category = "structure-factor"
Note: See TracChangeset
for help on using the changeset viewer.