IPhO 2023 — Задача 1. Characterization of Soil Colloids (10 points)

Автор: Olympiads XYZ · транскрипция на официалните материали

Проверена срещу оригинала на 13.9.2026 от същия модел, който я е транскрибирал (без независима проверка)

Съдържание

УсловиеРешение

Theory — Q1-1, English (Official): Characterization of Soil Colloids (10 points) — International Physics Olympiad 2023 Tokyo Japan · 10 т.

Условие

Colloidal science is useful to characterize soil particles because many of them can be regarded as colloidal particles of micrometer size. For example, Brownian motion (random motion of colloidal particles) can be used to measure particle sizes.

Part A. Motions of colloidal particles (1.6 points)

We analyze the one-dimensional Brownian motion of a colloidal particle with mass MM. The equation of motion for its velocity v(t)v(t) reads:

Mv˙=γv(t)+F(t)+Fext(t),(1)M\dot{v} = -\gamma v(t) + F(t) + F_{\mathrm{ext}}(t), \qquad (1)

where γ\gamma is the friction coefficient, F(t)F(t) is a force due to random collisions with water molecules, and Fext(t)F_{\mathrm{ext}}(t) is an external force. In Part A, we assume Fext(t)=0F_{\mathrm{ext}}(t) = 0.

In the following, you may use τ\tau in your answers.

Part B. Effective equation of motion (1.8 points)

Results so far imply that particle velocities v(t)v(t) and v(t)v(t') may be regarded as uncorrelated random quantities if ttτ|t - t'| \gg \tau. On this basis, we introduce a theoretical model to approximately describe the one-dimensional Brownian motion, where the velocity changes randomly at each time interval δ (τ)\delta\ (\gg \tau), i.e.,

v(t)=vn(tn1<ttn),(2)v(t) = v_n \quad (t_{n-1} < t \leq t_n), \qquad (2)

with tn=nδ (n=0,1,2,)t_n = n\delta\ (n = 0, 1, 2, \cdots) and a random quantity vnv_n. It satisfies

vn=0,vnvm={C(n=m),0(nm),(3)\langle v_n \rangle = 0, \quad \langle v_n v_m \rangle = \begin{cases} C & (n = m), \\ 0 & (n \neq m), \end{cases} \qquad (3)

with a parameter CC depending on δ\delta. Here X\langle X \rangle indicates the expectation value of XX. That is, if you draw random numbers XX infinite times, the mean will be X\langle X \rangle.

Now we consider the particle displacement Δx(t)=x(t)x(0)\Delta x(t) = x(t) - x(0) for t=Nδt = N\delta with an integer NN.

Part C. Electrophoresis (2.7 points)

Here we discuss electrophoresis, i.e., transport of charged particles by an electric field. Suspension of colloidal particles with mass MM and charge Q (>0)Q\ (> 0) is put in a narrow channel with a cross-section AA (Fig.1(a)). We ignore the interaction between particles, effects of the wall, the fluid, the ions therein, and gravity.

By applying a uniform electric field EE in the xx-direction, particles are transported and their concentration n(x)n(x) (particle number per unit volume) becomes non-uniform (Fig.1(b)). When EE is removed, this non-uniformity gradually disappears. This is due to Brownian motion of particles. If n(x)n(x) is not uniform, the numbers of right-going and left-going particles may differ (Fig.1(c)). This generates a particle flux JD(x)J_D(x), the mean number of particles flowing at xx along the xx-axis per unit cross-sectional area and unit time. This flux is known to satisfy

JD(x)=Ddndx(x),(4)J_D(x) = -D\,\frac{dn}{dx}(x), \qquad (4)

where DD is called the diffusion coefficient.

Now let's assume, for simplicity, that half of the particles have velocity +v+v and the other half have velocity v-v. Let N+(x0)N_+(x_0) be the number of particles with velocity +v+v that cross x0x_0 from left to right per unit cross-sectional area and unit time. For particles with velocity +v+v to cross x0x_0 in the time interval δ\delta, they should be in the shaded region of Fig.1(c). Since δ\delta is small, we have n(x)n(x0)+(xx0)dndx(x0)n(x) \simeq n(x_0) + (x - x_0)\frac{dn}{dx}(x_0) in this region.

We define N(x0)N_-(x_0) as the counterpart of N+(x0)N_+(x_0) for the velocity v-v. With this, we have JD(x0)=N+(x0)N(x0)J_D(x_0) = \langle N_+(x_0) - N_-(x_0) \rangle. According to Eq.(3), we have v2=C\langle v^2 \rangle = C.

Now we discuss the effect of osmotic pressure Π\Pi. It is given by Π=nNART=nkT\Pi = \frac{n}{N_A}RT = nkT with the Avogadro constant NAN_A, the gas constant RR, temperature TT, and the Boltzmann constant k=RNAk = \frac{R}{N_A}. Let us consider the non-uniform concentration formed under the electric field EE (Fig.1(b)). Since n(x)n(x) depends on xx, so does Π(x)\Pi(x). Then the forces due to Π(x)\Pi(x) and Π(x+Δx)\Pi(x+\Delta x) must be balanced with the total force from the field EE acting on the particles (Fig.2). Here we consider small Δx\Delta x, so that n(x)n(x) can be regarded as constant over this range, while n(x+Δx)n(x)Δxdndx(x)n(x+\Delta x) - n(x) \simeq \Delta x\,\frac{dn}{dx}(x).

Let us discuss the balance of the flux now. Besides the flux JD(x)J_D(x) due to the Brownian motion, there is also a flux due to the electric field, JQ(x)J_Q(x). It is given by

JQ(x)=n(x)u,(5)J_Q(x) = n(x)u, \qquad (5)

where uu is the terminal velocity of particles driven by the field.

Part D. Mean square displacement (2.4 points)

Suppose we observed the Brownian motion of an isolated, spherical colloidal particle with radius a=5.0 μma = 5.0\ \mu\mathrm{m} in water. Figure 3 shows the histogram of displacements Δx\Delta x measured in the xx-direction at every interval Δt=60 s\Delta t = 60\ \mathrm{s}. The friction coefficient is given by γ=6πaη\gamma = 6\pi a\eta with water viscosity η=8.9×104 Pas\eta = 8.9 \times 10^{-4}\ \mathrm{Pa \cdot s} and the temperature was T=25 CT = 25\ {}^{\circ}\mathrm{C}.

Δx (μm)\Delta x\ (\mu\mathrm{m})10-109-98-87-76-65-54-4
NcountN_{\mathrm{count}}000241115
Δx (μm)\Delta x\ (\mu\mathrm{m})3-32-21-100112233
NcountN_{\mathrm{count}}35445868595026
Δx (μm)\Delta x\ (\mu\mathrm{m})4455667788991010
NcountN_{\mathrm{count}}15543010

Now we extend the model in Part B to describe the motion of a particle with charge QQ under an electric field EE. The particle velocity v(t)v(t) considered in Eq.(2) should be replaced by v(t)=u+vn (tn1<ttn)v(t) = u + v_n\ (t_{n-1} < t \leq t_n) with vnv_n satisfying Eq. (3) and uu being the terminal velocity considered in Eq.(5).

Next, we consider swimming microbes (Fig.4(a)), in one dimension for simplicity (Fig.4(b)). These are spherical particles with radius aa. They swim at velocity either +u0+u_0 or u0-u_0, the sign chosen randomly at every time interval δ0\delta_0 without correlation. The observed motion is a combination of displacements due to swimming and those due to the Brownian motion of a spherical particle.

Part E. Water purification (1.5 points)

Here we discuss the purification of water including colloid-like soil particles, by adding electrolytes to coagulate them. Particles interact through van der Waals force and electrostatic force, the latter including effects of both surface charges and the surrounding layer of oppositely charged ions (such ions and their layer are called counter-ions and the electric double layer, respectively; see Fig.6(a)). As a result, the interaction potential for particle distance dd (Fig.6(b)) is given by

U(d)=Ad+Bϵ(kT)2q2ed/λ,(6)U(d) = -\frac{A}{d} + \frac{B\epsilon(kT)^2}{q^2}e^{-d/\lambda}, \qquad (6)

where AA and BB are positive constants, ϵ\epsilon is the dielectric constant of water, and λ\lambda is the thickness of the electric double layer. Assuming that charges of ions are ±q\pm q, we have

λ=ϵkT2NAq2e,(7)\lambda = \sqrt{\frac{\epsilon kT}{2N_A q^2 e}}, \qquad (7)

where ee is the molar concentration of ion.

Three panels: (a) a 3D rectangular channel along the x-axis containing spherical colloidal particles, with cross-section A marked at the right end; (b) a flat channel with an arrow E pointing in the +x direction and a concentration gradient shaded from dark (left) to light (right); (c) a channel with circles representing particles, arrows showing exchanges across a shaded vertical strip of width vδ around x₀, with an arrow vδ above the strip.
Fig.1: Setting for Part C.
A rectangle representing a segment of the channel between positions x and x+Δx (dashed lines). An arrow E points to the right above the box. Arrow labeled Π(x)A points right into the segment at x, arrow labeled Π(x+Δx)A points left into the segment at x+Δx, and a right-pointing arrow between them represents the electric force.
Fig.2: Force balance.
Histogram of displacement Δx (μm) from −10 to 10 in 1 μm bins; vertical axis N_count from 0 to 70. The distribution is bell-shaped, peaked at Δx = 0 with about 68 counts.
Fig.3: Histogram of displacements.
(a) Five green circular microbes with flagella connected by blue line segments showing successive swimming runs in different directions, ending in a blue arrow. (b) A one-dimensional x-axis with a microbe and horizontal double-headed arrows of various lengths representing alternating +u₀ and −u₀ swims.
Fig.4: (a) Motion of microbes. (b) Its one-dimensional version.
Log-log plot of ⟨Δx(t)²⟩ (μm²) versus t (sec) from 10⁻³ to 10² s. An orange curve rises from about 10⁻⁴ μm² to about 10³ μm², with three dashed guide lines of different slopes: shallow at small t, steep (slope 2) at intermediate t, and intermediate slope at large t.
Fig.5: Mean square displacement of the microbes.
(a) Two grey spheres each carrying plus signs on their surfaces, surrounded by circled minus signs (counter-ions) forming an electric double layer. (b) Two grey spheres separated by a distance d indicated with a double arrow.
Fig.6: (a) Surface charges of colloidal particles and counter-ions. (b) Definition of the distance d.

A.1 Consider that a water molecule collides with the particle at t=t0t = t_0, giving impulse I0I_0, and F(t)=0F(t) = 0 afterward. If v(t)=0v(t) = 0 before the collision, v(t)=v0e(tt0)/τv(t) = v_0 e^{-(t-t_0)/\tau} for t>t0t > t_0. Determine v0v_0 and τ\tau, using I0I_0 and necessary parameters in Eq.(1). [0,8 т.]

A.2 Actually, water molecules collide with the particle one after another. Suppose the iith collision gives the impulse IiI_i at time tit_i and determine v(t)v(t) on condition that t>0t > 0 and v(0)=0v(0) = 0. Also give the inequality specifying the range of tit_i that needs to be considered for a given tt. In the answer sheet, it is not necessary to specify this range in the expression for v(t)v(t). [0,8 т.]

B.1 Determine Δx(t)\langle \Delta x(t) \rangle and Δx(t)2\langle \Delta x(t)^2 \rangle using CC, δ\delta, and tt. [1 т.]

B.2 The quantity Δx(t)2\langle \Delta x(t)^2 \rangle is called the mean square displacement (MSD). It is a characteristic observable of the Brownian motion, which corresponds to the limiting case δ0\delta \to 0. From this, we can show CδαC \propto \delta^{\alpha} and Δx(t)2tβ\langle \Delta x(t)^2 \rangle \propto t^{\beta}. Determine the values of α\alpha and β\beta. [0,8 т.]

C.1 Express N+(x0)N_+(x_0) using necessary quantities from vv, δ\delta, n(x0)n(x_0), and dndx(x0)\frac{dn}{dx}(x_0). [0,5 т.]

C.2 Determine JD(x0)J_D(x_0) using necessary quantities from CC, δ\delta, n(x0)n(x_0), and dndx(x0)\frac{dn}{dx}(x_0). Using this and Eq.(4), express DD in terms of CC and δ\delta, and Δx(t)2\langle \Delta x(t)^2 \rangle in terms of DD and tt. [0,7 т.]

C.3 Express dndx(x)\frac{dn}{dx}(x) using n(x)n(x), TT, QQ, EE, and kk. [0,5 т.]

C.4 To determine uu, we use Eq.(1) with Fext(t)=QEF_{\mathrm{ext}}(t) = QE. Since v(t)v(t) is fluctuating, we consider v(t)\langle v(t) \rangle. Assuming v(0)=0\langle v(0) \rangle = 0 and using F(t)=0\langle F(t) \rangle = 0, evaluate v(t)\langle v(t) \rangle and obtain u=limtv(t)u = \lim_{t \to \infty} \langle v(t) \rangle. [0,5 т.]

C.5 The flux balance reads JD(x)+JQ(x)=0J_D(x) + J_Q(x) = 0. Express the diffusion coefficient DD in terms of kk, γ\gamma, and TT. [0,5 т.]

D.1 Estimate the value of NAN_A without using the fact that it is the Avogadro constant, up to two significant digits from the data in Fig.3. The gas constant is R=8.31 J/KmolR = 8.31\ \mathrm{J/K \cdot mol}. Do not use the value of the Boltzmann constant kk given in General Instructions. As for the Avogadro constant, you might obtain a value different from that in General Instructions. [1 т.]

D.2 Express the MSD Δx(t)2\langle \Delta x(t)^2 \rangle in terms of uu, DD, and tt. Obtain approximate power laws for small tt and large tt, as well as the characteristic time tt_* where this change occurs. Draw a rough graph of MSD in a log-log plot, indicating the approximate location of t.t_*. [0,8 т.]

D.3 Figure 5 displays the MSD Δx(t)2\langle \Delta x(t)^2 \rangle of those microbes, showing different power laws for small, large, and intermediate tt, as indicated by dashed lines. Obtain the power law for each time range and express it using necessary quantities from DD, u0u_0, δ0\delta_0, and tt. [0,6 т.]

E.1 Addition of sodium chloride (NaCl) to the suspension causes colloidal particles to coagulate. Determine the lowest concentration cc of NaCl necessary for coagulation. It is sufficient to consider two particles without thermal fluctuations, i.e., F(t)=0F(t) = 0 in Eq.(1), and assume that the terminal velocity for the given potential force is reached instantaneously. [1,5 т.]

Решение

Покажи официалното решение


Оригинал в Архива: IPhO_2023_Q1.pdf · официални решения: IPhO_2023_S1.pdf