[dc02af0] | 1 | /* Lamellar_ParaCrystal - Pedersen's model |
---|
| 2 | |
---|
| 3 | */ |
---|
[a0168e8] | 4 | double paraCryst_sn(double ww, double qval, double davg, int Nlayers, double an); |
---|
| 5 | double paraCryst_an(double ww, double qval, double davg, int Nlayers); |
---|
[dc02af0] | 6 | |
---|
[a0168e8] | 7 | static double |
---|
| 8 | Iq(double qval, |
---|
| 9 | double th, |
---|
| 10 | double fp_Nlayers, |
---|
| 11 | double davg, |
---|
| 12 | double pd, |
---|
| 13 | double sld, |
---|
| 14 | double solvent_sld) |
---|
[dc02af0] | 15 | { |
---|
| 16 | //get the fractional part of Nlayers, to determine the "mixing" of N's |
---|
[a0168e8] | 17 | int n1 = (int)(fp_Nlayers); //truncate towards zero |
---|
| 18 | int n2 = n1 + 1; |
---|
| 19 | const double xn = (double)n2 - fp_Nlayers; //fractional contribution of n1 |
---|
[dc02af0] | 20 | |
---|
[a0168e8] | 21 | const double ww = exp(-0.5*square(qval*pd*davg)); |
---|
[319ab14] | 22 | |
---|
[dc02af0] | 23 | //calculate the n1 contribution |
---|
[a0168e8] | 24 | double Znq,Snq,an; |
---|
[dc02af0] | 25 | an = paraCryst_an(ww,qval,davg,n1); |
---|
| 26 | Snq = paraCryst_sn(ww,qval,davg,n1,an); |
---|
[a0168e8] | 27 | |
---|
[dc02af0] | 28 | Znq = xn*Snq; |
---|
| 29 | |
---|
| 30 | //calculate the n2 contribution |
---|
| 31 | an = paraCryst_an(ww,qval,davg,n2); |
---|
| 32 | Snq = paraCryst_sn(ww,qval,davg,n2,an); |
---|
[319ab14] | 33 | |
---|
[dc02af0] | 34 | Znq += (1.0-xn)*Snq; |
---|
| 35 | |
---|
| 36 | //and the independent contribution |
---|
| 37 | Znq += (1.0-ww*ww)/(1.0+ww*ww-2.0*ww*cos(qval*davg)); |
---|
| 38 | |
---|
| 39 | //the limit when Nlayers approaches infinity |
---|
| 40 | // Zq = (1-ww^2)/(1+ww^2-2*ww*cos(qval*davg)) |
---|
| 41 | |
---|
[a0168e8] | 42 | const double xi = th/2.0; //use 1/2 the bilayer thickness |
---|
[5467cd8] | 43 | const double Pbil = square(sas_sinx_x(qval*xi)); |
---|
[dc02af0] | 44 | |
---|
[a0168e8] | 45 | const double contr = sld - solvent_sld; |
---|
| 46 | const double inten = 2.0*M_PI*contr*contr*Pbil*Znq/(qval*qval); |
---|
[319ab14] | 47 | //printf("q=%.7e wwm1=%g ww=%.5e an=% 12.5e Snq=% 12.5e Znq=% 12.5e Pbil=% 12.5e\n",qval,wwm1,ww,an,Snq,Znq,Pbil); |
---|
[a0168e8] | 48 | return 1.0e-4*inten; |
---|
[dc02af0] | 49 | } |
---|
| 50 | |
---|
| 51 | // functions for the lamellar paracrystal model |
---|
| 52 | double |
---|
[a0168e8] | 53 | paraCryst_sn(double ww, double qval, double davg, int Nlayers, double an) { |
---|
[dc02af0] | 54 | |
---|
| 55 | double Snq; |
---|
[319ab14] | 56 | |
---|
[4962519] | 57 | Snq = an/( (double)Nlayers*square(1.0+ww*ww-2.0*ww*cos(qval*davg)) ); |
---|
[dc02af0] | 58 | |
---|
[a0168e8] | 59 | return Snq; |
---|
[dc02af0] | 60 | } |
---|
| 61 | |
---|
| 62 | double |
---|
[a0168e8] | 63 | paraCryst_an(double ww, double qval, double davg, int Nlayers) { |
---|
[dc02af0] | 64 | double an; |
---|
| 65 | |
---|
| 66 | an = 4.0*ww*ww - 2.0*(ww*ww*ww+ww)*cos(qval*davg); |
---|
| 67 | an -= 4.0*pow(ww,(Nlayers+2))*cos((double)Nlayers*qval*davg); |
---|
| 68 | an += 2.0*pow(ww,(Nlayers+3))*cos((double)(Nlayers-1)*qval*davg); |
---|
| 69 | an += 2.0*pow(ww,(Nlayers+1))*cos((double)(Nlayers+1)*qval*davg); |
---|
| 70 | |
---|
[a0168e8] | 71 | return an; |
---|
[dc02af0] | 72 | } |
---|
| 73 | |
---|