0% found this document useful (0 votes)
17 views4 pages

Thermo Stat

This document contains 5 problems related to statistical thermodynamics of complex liquids. Problem 1 asks to calculate scattering length densities of various molecules. Problem 2 derives the form factor of a core-shell particle. Problem 3 derives the form factor of a prolate spheroid. Problem 4 derives the normalized structure factor of a cylindrical particle and explores asymptotic limits. Problem 5 (optional) derives the scattering intensity of a Gaussian chain.

Uploaded by

Hagere Ethiopia
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
17 views4 pages

Thermo Stat

This document contains 5 problems related to statistical thermodynamics of complex liquids. Problem 1 asks to calculate scattering length densities of various molecules. Problem 2 derives the form factor of a core-shell particle. Problem 3 derives the form factor of a prolate spheroid. Problem 4 derives the normalized structure factor of a cylindrical particle and explores asymptotic limits. Problem 5 (optional) derives the scattering intensity of a Gaussian chain.

Uploaded by

Hagere Ethiopia
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

1

8.575J , 10.44J, 22.52J Statistical Thermodynamics of Complex Liquids


(Spring 2004) Problem Set 1 (Prof. Chen) Due March 18.

1. Given the following table for the scattering lengths of common elements:

Isotope Hydrogen Deuterium Carbon Nitrogen Oxygen


bcoh (10-12 cm) -0.37423 0.6674 0.66484 0.936 0.5805

estimate the molecular volume (from the molecular weight and density) and then calculate the scattering
length densities of the following molecules, in unit of1010 cm-2: H2O, D2O, Octane, Deuterated octane,
and Pluronic P-84, a tri-block co-polymer, [(PEO)19 (PPO)43 (PEO)19], where PEO = -(CH2)2O-, having
a molecular volume 72.4 Å3, and PPO = -(CH2)3O-, having a molecular volume 95.4 Å3 .

2. Show that the form factor of a spherical particle with an internal core of radius R1 and a scattering
length density (sld) ρ1, surrounded by a shell with an outer radius R2 and sld ρ2, immersed in a solvent
of sld ρs, is given by:

4 3 ⎡3 j1(QR1 ) ⎤ 4 3 ⎡ 3j (QR ) ⎤
F2−shell (Q) = πR1 (ρ1 − ρ2 ) ⎢ ⎥ + πR 2 (ρ2 − ρs )⎢ 1 2

3 ⎣ QR 1 ⎦ 3 ⎣ QR 2 ⎦

Use this result to calculate and plot the normalized particle structure factor P(Q) of a co-polymer
micelle having an inner core radius R1 and an outer radius R2. In a core-shell model of the micelle [Y.C.
Liu et al, Phys. Rev. E 54, 1698 (1996)], the inner and outer radii can be determined from the
aggregation (N) and hydration (H) numbers of the micelle. Plot the P(Q) for the case of N = 63 and H
= 290.

3. Show that the form factor of a randomly oriented prolate spheroid, with semi-major and minor axes
of a and b, is given by:

4π 2 ⎡ 3j (u) ⎤
Fellipsoid (Q) = ab (∆ρ)∫0 dµ ⎢ 1 ⎥
1

3 ⎣ u ⎦ (2)
u = Q a µ + b (1 − µ )
2 2 2 2

where (∆ρ) is the contrast between the particle and the solvent, µ the cosine of the angle between the
major axis of the spheroid and the Q-vector.
2

4. Derive a normalized particle structure factor P(Q) of a uniform cylindrical particle of radius R and
r
length L. Assuming that the particle is randomly oriented with respect to the Q vector.
(A) Show that:
1
2 2⎡ 2 ⎤
2
1 vv
1 ⎡ sin QLµ / 2 ⎤ 2J (QR 1 − µ )
⎥ ⎢⎢ 1 ⎥ . (3)
iQ ⋅ r 3
P(Q) =< ∫ e d r >= ∫−1 dµ ⎢
1
Vp Vp 2 ⎣ QLµ / 2 ⎦ ⎣ QR 1 − µ 2 ⎥⎦

v
In Eq.3, V p denotes the volume of the particle, r the position vector of an arbitrary point in the
interior of the particle, and µ the cosine of the angle between the axis of the cylinder and the Q-
vector. The bracket means that we are considering an average over random orientations of the
particle.
(B) Show that for a long and thin cylinder, one has asymptotic formulae:

π ⎡ 2J 1 (QR) ⎤ π
2 1 2 2
− Q R
P(Q) →
QL>2 π QL ⎢⎣ QR ⎥⎦ → QL e
QR <1
4
(4)

(C) Show that for a flat particle (a lamellar) of QR>>1,

⎡ sin QL / 2 ⎤
2 1
2 2 − Q L
2 2

P(Q) → 2 2
QR>>1 Q R
⎢⎣ QL / 2 ⎥⎦ → 2 2e
QL<1 Q R
12
(5)

where L is the thickness of the flat plate.

(D) From Eq.4 and 5, one can conclude that for a long rod a ln[QI(Q)] vs Q2 plot, and for a flat disk,
a ln[Q2I(Q)] vs Q2 plot, will result in a straight line at large Q with slopes proportional to R2 /4 and
L2/12 respectively. Explore additional system parameters you can extract from the intercept at Q =
0.
(E) In polymer literature, another approximate formula is often used. It is the limit when R goes to
zero, the so called "stiff thin rod" limit. Show that:

⎡ 2⎤
sin QLµ /
2
⎡ sin xµ ⎤ 2 1 2 x sin u ⎛ sin x ⎞ 2
⎥ = ∫0 dµ ⎢ ⎥ = ∫0 du − ⎜ ⎟ . (6)
1 1
P(Q) = ∫0 dµ ⎢
R→0 ⎣ QLµ / 2 ⎦ ⎣ xµ ⎦ x u ⎝ x ⎠

Explore graphically the difference between approximations of the last equation and Eq. 3
3

5. (This problem is added for your interest only. Your answer is optional. In case you solve it,
you will get a bonus points of 20/100)
Scattering intensity of a Gaussian chain. Consider a flexible polymer chain of a contour length L
=Na, where N is the no. of segments and a is the Kuhn length. If the chain makes a random walk in
space, then the mean square end-to-end distance is R = Na . For this chain, the distribution of
2 2

distances between two links (i,j) is Gaussian, namely,


2
⎡ 3 ⎤ ⎡ −3Rij2 ⎤
3/ 2

P(Rij )dRij = ⎢ ⎥ 4 πR exp ⎢2


⎥ dRij , where Rij2 = i − j a 2 . (7)
⎢⎣ 2 π Rij2 ⎥⎦ ⎢⎣ 2 Rij ⎥⎦
ij 2

In order to calculate the normalized particle structure factor of such chain. (a) Start from the
definition:

N N r r r N N r r

∑∑ ∑∑
1 iQ⋅( R i − R j ) 1 iQ⋅ Rij
P (Q) = 2 e = 2 e (8)
N i =1 j =1
N i=1 j =1
gaussian

so that we can evaluate the Gaussian average by integrating the exponential phase factor using the
distribution function given by Eq.7. (b) Show first that :

N N N N

∑∑ ∑ ∑ exp[− 6Q a
1 1 2 1 1
P (Q) = 2 exp[− Rij ] = 2 i− j
2 2
(9)
N i=1 j =1
6 N i =1 j =1

(c) Prove a theorem: For an arbitrary function f(x),


N N N −1

∑ ∑ f (|i − j|) = Nf (0) + 2 ∑ (N − n) f (n)


i =1 j =1 n=1
(10)

(d) Use the theorem to show that the sum in Eq.9 can be evaluated as:

1 ⎡ (N − 1)α − Nα 2 + α N +1 ⎤
P (Q) = 2 ⎢N + 2 ⎥ (11)
N ⎣ (1 − α )2 ⎦

1 1
where α = exp( Q 2a 2 ) ≈ 1 − Q 2 a 2 , because Qa is much smaller than unity in practice.
6 6
(e) Show that in the limit N → ∞ , P (Q) approach the Debye function

2 1 2 2 1 2 2
P (Q) = (x −1 + e ) , where x = Q R = Q a N
−x
2
(12)
x 6 6
(f) Discuss the small and large Q behavior of Eq.12. In particular, show that:
4

1 1 1
lim = + Q2 R 2 (13)
x→ ∞ P(Q) 2 12

You might also like