Changeset a0fee3b in sasmodels


Ignore:
Timestamp:
Oct 6, 2016 1:48:59 PM (5 years ago)
Author:
richardh
Branches:
master, core_shell_microgels, costrafo411, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
Children:
4a2e1e5
Parents:
7ab7bb3
Message:

improved doc and new figures for core_shell_bicelle

Location:
sasmodels/models
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/core_shell_bicelle.py

    rb0c4271 ra0fee3b  
    1212.. figure:: img/core_shell_bicelle_geometry.png 
    1313 
    14     Core shell geometry (Graphic from ref [#Matusmori]_). 
    15     Note however that the model here calculates for rectangular, not curved, rims. 
     14    Schematic cross-section of bicelle. Note however that the model here 
     15    calculates for rectangular, not curved, rims as shown below. 
     16 
     17.. figure:: img/core_shell_bicelle_parameters.png 
     18 
     19   Cross section of cylindrical symmetry model used here. Users will have  
     20   to decide how to distribute "heads" and "tails" between the rim, face  
     21   and core regions in order to estimate appropriate starting parameters. 
    1622 
    1723The output of the 1D scattering intensity function for randomly oriented 
     
    5359title = "Circular cylinder with a core-shell scattering length density profile.." 
    5460description = """ 
    55     P(q,alpha)= scale/Vs*f(q)^(2) + bkg,  where: f(q)= 2(sld_core 
    56     - solvant_sld)* Vc*sin[qLcos(alpha/2)] 
    57     /[qLcos(alpha/2)]*J1(qRsin(alpha)) 
    58     /[qRsin(alpha)]+2(shell_sld-sld_solvent) 
    59     *Vs*sin[q(L+T)cos(alpha/2)][[q(L+T) 
    60     *cos(alpha/2)]*J1(q(R+T)sin(alpha)) 
    61     /q(R+T)sin(alpha)] 
     61    P(q,alpha)= (scale/Vs)*f(q)^(2) + bkg,  where:  
     62    f(q)= Vt(sld_rim - sld_solvent)* sin[qLt.cos(alpha)/2] 
     63    /[qLt.cos(alpha)/2]*J1(qRout.sin(alpha)) 
     64    /[qRout.sin(alpha)]+ 
     65    (sld_core-sld_face)*Vc*sin[qLcos(alpha)/2][[qL 
     66    *cos(alpha)/2]*J1(qRc.sin(alpha)) 
     67    /qRc.sin(alpha)]+ 
     68    (sld_face-sld_rim)*(Vc+Vf)*sin[q(L+2.thick_face). 
     69    cos(alpha)/2][[q(L+2.thick_face)*cos(alpha)/2]* 
     70    J1(qRc.sin(alpha))/qRc.sin(alpha)] 
    6271 
    6372    alpha:is the angle between the axis of 
    6473    the cylinder and the q-vector 
    65     Vs: the volume of the outer shell 
    66     Vc: the volume of the core 
     74    Vt = pi.(Rc + thick_rim)^2.Lt : total volume 
     75    Vc = pi.Rc^2.L :the volume of the core 
     76    Vf = 2.pi.Rc^2.thick_face 
     77    Rc = radius: is the core radius 
    6778    L: the length of the core 
    68     shell_sld: the scattering length density 
    69     of the shell 
     79    Lt = L + 2.thick_face: total length 
     80    Rout = radius + thick_rim 
     81    sld_core, sld_rim, sld_face:scattering length 
     82    densities within the particle 
    7083    sld_solvent: the scattering length density 
    7184    of the solvent 
    7285    bkg: the background 
    73     T: the thickness 
    74     R+T: is the outer radius 
    75     L+2T: The total length of the outershell 
    7686    J1: the first order Bessel function 
    7787    theta: axis_theta of the cylinder 
Note: See TracChangeset for help on using the changeset viewer.