Changeset b0c4271 in sasmodels for sasmodels/models/adsorbed_layer.py
- Timestamp:
- Oct 5, 2016 7:40:41 PM (8 years ago)
- Branches:
- master, core_shell_microgels, costrafo411, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- 63e9cd8
- Parents:
- bb2d187
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/adsorbed_layer.py
r52ec91e rb0c4271 1 #adsorbed_layer model 2 #conversion of Core2ndMomentModel.py 3 #converted by Steve King, Mar 2016 1 r""" 2 Definition 3 ---------- 4 4 5 r"""6 5 This model describes the scattering from a layer of surfactant or polymer 7 6 adsorbed on large, smooth, notionally spherical particles under the conditions … … 19 18 the mean of the density distribution (ie, the distance of the centre-of-mass 20 19 of the distribution from the interface), $\sigma = \sqrt{t^2/12}$. 21 22 Definition23 ----------24 20 25 21 .. math:: … … 44 40 calculation is exact). 45 41 42 The code for this model is based originally on a a fortran implementation by 43 Steve King at ISIS in the SANDRA package c. 1990. 44 46 45 References 47 46 ---------- 48 47 49 S King, P Griffiths, J Hone, and T Cosgrove, 50 *SANS from Adsorbed Polymer Layers*, *Macromol. Symp.*, 190 (2002) 33-42. 48 .. [#] S King, P Griffiths, J Hone, and T Cosgrove, *SANS from Adsorbed Polymer 49 Layers*, *Macromol. Symp.*, 190 (2002) 33-42. 50 51 Authorship and Verification 52 ---------------------------- 53 54 * **Author:** Jae-Hi Cho **Date:** pre 2010 55 * **Last Modified by:** Paul Kienzle **Date:** April 14, 2016 56 * **Last Reviewed by:** Steve King **Date:** March 18, 2016 51 57 """ 52 58 … … 80 86 def Iq(q, second_moment, adsorbed_amount, density_shell, radius, 81 87 volfraction, sld_shell, sld_solvent): 82 # pylint: disable = missing-docstring83 #deltarhosqrd = (sld_shell - sld_solvent) * (sld_shell - sld_solvent)84 #numerator = 6.0 * pi * volfraction * (adsorbed_amount * adsorbed_amount)85 #denominator = (q * q) * (density_shell * density_shell) * radius86 #eterm = exp(-1.0 * (q * q) * (second_moment * second_moment))87 ##scale by 10^-2 for units conversion to cm^-188 #inten = 1.0e-02 * deltarhosqrd * ((numerator / denominator) * eterm)89 88 with errstate(divide='ignore'): 90 89 aa = ((sld_shell - sld_solvent)/density_shell * adsorbed_amount) / q
Note: See TracChangeset
for help on using the changeset viewer.