الفيزياء الجامعية — السنة الثانية · Bachelor Year 2
24انتشار الجسيمات
ألقِ بلورة صباغ في ماء ساكن وراقب: تتكوّن حولها غيمة ملونة، وتكبر، وتلين حروفها، وتنتشر — في دقيقة على مليمتر، وفي ساعة على سنتيمتر، وفي أسبوع على الكأس كلها. ولا شيء يدفع الصباغ؛ بل يحمله التزاحم الدائم للجزيئات، الذي يرسل كل جسيم في مسير عشوائي ويحمله، في المتوسط، من حيث الكثرة إلى حيث القلّة. وهذا هو الانتشار، وهو أبطأ النقل وأعمّه: فهو يغذّي كل خلية بالأكسجين، ويشوب كل ترانزستور بالبور، ويقسّي الفولاذ بالكربون، ويحمل عطرًا عبر غرفة — وإن كان، كما سنرى، لا في الزمن الذي قد يُظنّ. وهذا الفصل يعطي الانتشار قانونه (فيك)، ومعادلته (من موازنة جسيمية)، وحلوله المميّزة وسلالمه الزمنية (L∼Dt)، وأصله المجهري في المسير العشوائي — وهو يفسّر أيضًا لمَ تعرف معادلة الانتشار، خلافًا لكل معادلات الميكانيك، جهةَ الزمن. وسيعيد الفصل التالي استعمال ذلك كله من أجل الحرارة.
حبر يُطلَق في ماء ساكن: تتشوّش الغيمة الحادّة وتنتشر كلما انتشرت جزيئاتها — على المليمترات في دقيقة، وعلى السنتيمترات في ساعة.
24.1 كثافة الجسيمات والتدفق وقانون فيك
تعريف 24.1(الكثافة وتيار الجسيمات)
من أجل نوع من الجسيمات (جزيئات، أو شوارد، أو ذرات في جسم صلب) تكون الكثافة العدديةn(M,t) هي عدد الجسيمات في وحدة الحجم حول M (بالوحدة m−3؛ والتركيز المولي هو c=n/NA). وكثافة تيار الجسيماتjN هي الشعاع الذي يجعل عدد الجسيمات العابرة عنصر سطح موجَّهًا dS في dt مساويًا jN⋅dSdt (بالوحدة m−2s−1): وتدفق الجسيمات عبر سطح S هو ΦN=∬SjN⋅dS، أي عدد الجسيمات في الثانية عبر S. ومن أجل جسيمات يحملها مائع يتحرك بالسرعة v، يكون jN=nv (وهو الحمل)؛ والانتشار هو النقل الذي يبقى في مائع ساكن.
مبرهنة 24.2(قانون فيك)
في وسط ساكن، حيث الكثافة غير منتظمة، يظهر تيار جسيمات، متناسب مع متدرج الكثافة ومعاكس له:
jN=−Dgradn,
حيث معامل الانتشارD>0 (بالوحدة m2/s) يتعلق بالنوع المنتشر، وبالوسط وبالحرارة. والجسيمات تنزل مع انحدار متدرج الكثافة، من المناطق الأكثف إلى الأرقّ، بمعدل يحدده D.
برهان. ظاهراتي (أي قانون تجربة، كقانون أوم): خطي في المتدرج من أجل المتدرجات الصغيرة، ومتماثل المناحي في وسط متماثل المناحي، وذو الإشارة التي تفرضها التجربة. ونموذج المسير العشوائي في القسم 24.4 يستنتجه، ويستنتج D معه، من الحركة الجزيئية. ∎
مثال 24.3(رتب مقادير D)
الغازات: D∼1×10−5m2/s (بخار الماء في الهواء 2.5×10−5m2/s، وجزيء عطر 5×10−6m2/s). والسوائل: D∼1×10−9m2/s (الأكسجين في الماء 2×10−9m2/s، والسكر 5×10−10m2/s، وبروتين 1×10−10m2/s). والأجسام الصلبة: ضئيل و يتصاعد بشدة مع الحرارة، D=D0e−Ea/kBT (البور في السيليسيوم: 1.5×10−17m2/s عند 1100∘C، وصغير صغرًا لا يُقاس في حرارة الغرفة — ولهذا يحفظ الترانزستور، متى صُنع، مقطع تشويبه عقودًا). أي أربع رتب من المقدار من الغاز إلى السائل، وثماني أو أكثر من السائل إلى الصلب.
24.2 الموازنة الجسيمية ومعادلة الانتشار
مبرهنة 24.4(الموازنة الجسيمية المحلية)
إذا لم تُحدَث الجسيمات ولم تُتلَف،
∂t∂n+divjN=0;
ومع منبع يحدث σ جسيمًا في وحدة الحجم والزمن (تفاعل كيميائي، أو امتصاص)، يكون ∂tn+divjN=σ. وفي بُعد واحد (حيث تتعلق الكثافة والتيار بالمتغير x وحده): ∂tn+∂xjN=0.
برهان. خذ الشريحة بين x و x+dx، ومقطعها S: فهي تحوي nSdx جسيمًا؛ وفي dt يدخل jN(x)Sdt عبر الوجه الأيسر ويخرج jN(x+dx)Sdt عبر الأيمن، ومن ثَمّ ∂tnSdxdt=−(jN(x+dx)−jN(x))Sdt=−∂xjNdxSdt. وفي ثلاثة أبعاد يعطي العدّ نفسه على صندوق صغير التباعد (أي التدفق الخارج من سطح مغلق في وحدة الحجم، الفصل 11)، أو مباشرةً: من أجل أي حجم ثابت V، d/dt∭Vndτ=−∬∂VjN⋅dS=−∭VdivjNdτ بمبرهنة أوستروغرادسكي. ∎
الموازنة في بُعد واحد: الشريحة بين x و x+dx تكسب ما يدخل من اليسار وتفقد ما يخرج من اليمين؛ والفرق هو تغيّر n في داخلها.
مبرهنة 24.5(معادلة الانتشار)
بجمع قانون فيك والموازنة، من أجل D منتظم،
∂t∂n=DΔn+σ,وفي بُعد واحد∂t∂n=D∂x2∂2n+σ.
ومعادلة الانتشار هذه خطية (فالحلول تتراكب)، من المرتبة الأولى في الزمن والثانية في الفضاء — وليست معادلة أمواج: فليس لها سرعة انتشار، بل علاقة مميّزة بين الطول والزمن،
L∼Dt,t∼DL2:
فالانتشار على ضعف المسافة يستغرق أربعة أضعاف الزمن.
برهان.∂tn=−div(−Dgradn)+σ=Ddivgradn+σ=DΔn+σ. وأما السلالم: فقارن n/t بالمقدار Dn/L2. ∎
ملاحظة 24.6(الانتشار غير عكوس)
بدّل t بالمقدار −t في معادلة الأمواج، ∂t2u=c2∂x2u، فلا تتغير: ففيلم موجة يُدار بالمقلوب موجة ممكنة. وافعل ذلك في معادلة الانتشار فتنقلب إشارة الطرف الأيسر: فإن n(x,−t)ليس حلًّا. والغيمة التي تنتشر طبيعية؛ أما التي تتجمع تلقائيًّا في بلورة فلا — فللانتشار سهم زمن، وهو سهم المبدأ الثاني (والفصل 25 يحسب الإنتروبيا التي يحدثها). وتبيّن المعادلة أيضًا أن الانتشار يملّس: فحيث تكون n محليًّا عظمى (∂x2n<0) تتناقص، وحيث تكون صغرى تتزايد.
مثال 24.7(كم يستغرق ذلك؟)
t∼L2/D. والأكسجين عبر خلية، L=10µm في الماء: 10−10/2×10−9=0.05s — فالانتشار يغذّي الخلية بيسر؛ وعبر 1mm من نسيج: 500s — أي بطيء جدًّا، ولهذا لا يكون شيء حيّ أثخن من كسر من المليمتر بلا أوعية دموية. والسكر عبر فنجان شاي غير مقلّب، L=5cm: 2.5×10−3/5×10−10=5×106s، أي شهران — فقلّبه. وعطر عبر غرفة بالانتشار وحده، 5m في الهواء: 25/10−5=3×106s، أي شهر؛ وأنت تشمّه في دقيقة لأن الهواء يتحرك: والحمل ينقل، والانتشار لا يفعل إلا إتمام المليمترات الأخيرة.
24.3 الحلول المميّزة
قضية 24.8(النظام المستقر: الغشاء)
بين خزّانين محفوظين عند n1 و n2، عبر غشاء سماكته e ومساحته S، تكون الكثافة المستقرة خطية،
فللغشاء مقاومة انتشاريةe/DS، تمامًا كما للسلك ℓ/γS — والمقاومات على التسلسل تتجمع، وغشاء رقيق صغير D قد يغلب مع ذلك.
برهان. مستقر، وفي بُعد واحد، وبلا منبع: ∂x2n=0، فتكون n تآلفية؛ و jN=−Ddn/dx منتظم. ∎
قضية 24.9(النظام المستقر: الهندسة الكروية)
حول كرة نصف قطرها a يُحفَظ سطحها عند n0، في وسط تؤول فيه n→0 بعيدًا، تكون الكثافة المستقرة
n(r)=n0ra,ΦN=4πDan0:
فالتيار الكلي المطلَق (أو الممتصّ، بقلب الإشارات) ينمو مع نصف قطر الكرة لا مع مساحتها — فهندسة المصبّ الانتشاري هي هندسة سعة كهربائية ساكنة.
برهان.Δn=r21drd(r2dn/dr)=0 يعطي n=A+B/r؛ والشروط الحدية تثبّت A=0 و B=n0a؛ ثم ΦN=−4πr2Ddn/dr=4πDn0a عند كل r (فلا تراكم في النظام المستقر). ∎
المقاطع المستقرة. على اليسار: عبر غشاء تكون الكثافة خطية والتيار منتظمًا. وعلى اليمين: حول كرة تتناقص الكثافة وفق 1/r؛ والتيار الكلي 4πDan0 متناسب مع نصف القطر.
قضية 24.10(الغاوسي المنتشر)
N جسيمًا في وحدة المساحة تُطلَق عند t=0 على المستوي x=0 من وسط لانهائي تنتشر وفق
n(x,t)=4πDtNexp(−4Dtx2):
أي غاوسي انحرافه المعياري σ(t)=2Dt، ومساحته ثابتة N (فالجسيمات محفوظة) وقمته متناقصة ∝1/t. وفي ثلاثة أبعاد، تنتشر N جسيمًا تُطلَق عند نقطة وفق n=N(4πDt)−3/2exp(−r2/4Dt)، مع ⟨r2⟩=6Dt. ونبضة جسيمات تُطلَق عند سطح نصف فضاء (لا تستطيع مغادرته) تعطي ضعف العبارة أعلاه من أجل x>0.
برهان. عوّض: مع u=x2/4Dt، يكون ∂tn=n(−1/2t+u/t) و D∂x2n=Dn(−1/2Dt+x2/4D2t2)=n(−1/2t+u/t)؛ فهما متساويان. والتنظيم ∫ndx=N ينتج من ∫e−x2/4Dtdx=4πDt، و ∫x2ndx/N=2Dt. وأما الوحدانية (بالنبضة الابتدائية المعطاة) فمقبولة. ∎
الحلّ الغاوسي عند ثلاثة أزمنة: العرض ينمو وفق t، والقمة تهبط وفق 1/t، والمساحة (أي عدد الجسيمات) تبقى نفسها.
قضية 24.11(تركيز سطحي ثابت)
نصف فضاء x>0 فارغ ابتداءً، ويُحفَظ سطحه عند الكثافة n0 منذ t=0 (كغاز في تماس مع جسم صلب يحلّه)، يمتلئ وفق
n(x,t)=n0erfc(2Dtx),erfc(u)=π2∫u∞e−s2ds,
وهو تابع الخطأ المتمّم، الذي يهبط من 1 عند u=0 إلى 0.16 عند u=1 و 0.005 عند u=2: فعمق النفوذ هو ∼2Dt من جديد، والعدد الكلي للجسيمات الممتصّة في وحدة المساحة هو 2n0Dt/π.
برهان. ابحث عن n=f(u) مع u=x/2Dt: فتصير المعادلة f′′+2uf′=0، ومن ثَمّ f′∝e−u2 وتكون f تابع خطأ؛ والشرطان f(0)=n0 و f(∞)=0 ينتقيان erfc. والعدد الممتصّ هو ∫0∞ndx=n02Dt∫0∞erfc(u)du=2n0Dt/π. ∎
مثال 24.12(تقسية الفولاذ)
قطعة فولاذية تُحفَظ عند 900∘C في غاز غني بالكربون تمتصّ الكربون عند سطحها؛ ومع D=5×10−12m2/s للكربون في الحديد الساخن، تعطي أربع ساعات 2Dt=0.5mm: أي قشرة قاسية بنصف مليمتر على لبّ متين — وهي تقسية السطح في المسنّنات والمحامل، مؤقَّتة بمقطع تابع الخطأ المتمّم.
24.4 الصورة المجهرية: المسير العشوائي
قضية 24.13(المسير العشوائي ومعامل الانتشار)
جسيم يخطو خطوة طولها ℓ في جهة عشوائية كل τ (كجزيء بين تصادمين) له، بعد N=t/τ خطوة، متوسط مربع انتقال
D=2τℓ2(بُعد واحد),D=6τℓ2=6ℓv∗∼3ℓv∗(ثلاثة أبعاد، مع v∗=ℓ/τ؛ والنظرية الحركية تعطي 31).
ومسافة الماشي تنمو وفق t لا وفق t: فليذهب ضعف البعد يحتاج إلى أربعة أضعاف عدد الخطوات.
برهان.x=∑iϵiℓ مع ϵi=±1 مستقلة (في بُعد واحد): ⟨x2⟩=∑i,j⟨ϵiϵj⟩ℓ2=Nℓ2، إذ تتوسط الحدود المتقاطعة إلى الصفر. وتوزيع x من أجل N الكبيرة يؤول إلى غاوسي (وهي مبرهنة النهاية المركزية، وبرهانها من درس الاحتمالات) — ولهذا يكون القانون العياني معادلة الانتشار؛ وأما قانون فيك نفسه فينتج من عدّ الماشين العابرين مستويًا من الجهتين: jN≈−21ℓv∗∂xn (ومتوسط أدقّ يعطي 1/3 في ثلاثة أبعاد). ∎
مسير عشوائي من 400 خطوة طولها ℓ: لم يهم الماشي إلا نحو 20ℓ عن بدايته — وهو قانون N الذي يجعل الانتشار بطيئًا على المسافات الطويلة وسريعًا على القصيرة.
مثال 24.14(من الجزيئات إلى D)
الهواء في حرارة الغرفة: المسار الحر الوسطي ℓ≈70nm، والسرعة الوسطى v∗≈500m/s: D≈ℓv∗/3=1.2×10−5m2/s، كما يُقاس. ولينتشر جزيء 1m يحتاج إلى N=(L/ℓ)2=2×1014 تصادم، أي Nℓ/v∗≈3×104s — أي ثماني ساعات لرحلة كان يقطعها في 2ms لو طار مستقيمًا. وفي سائل، تحقق كرة نصف قطرها a تتقاذفها الجزيئات علاقة ستوكس–أينشتاين D=kBT/6πηa (أينشتاين 1905): فمن أجل السكر، a≈0.4nm في الماء، D=5×10−10m2/s؛ ومن أجل حبيبة 1µm، 4×10−13m2/s، أي انتقال ميكرومتر في الثانية — وهي الحركة البراونية التي قاسها بيران ليعدّ عدد أفوغادرو. وفي جسم صلب لا تقفز الذرة إلى موقع مجاور إلا حين يزوّدها تقلّب حراري بطاقة التنشيط: D=D0e−Ea/kBT، ويتضاعف كل بضع عشرات من الكلفنات قرب 1000∘C.
طريقة 24.15(تقديرات الانتشار)
(1) عيّن D (الغاز 10−5، والسائل 10−9، والصلب بأرينيوس). (2) والزمن أو المسافة: t∼L2/D، و L∼Dt (مع المضبوط 2Dt لعرض الغاوسي، و 2Dt لعمق تابع الخطأ المتمّم). (3) والمستقر: خطي في شريحة (Rd=e/DS)، و 1/r حول كرة (Φ=4πDaΔn)، و lnr حول أسطوانة. (4) والعابر: غاوسي من أجل نبضة، وتابع خطأ متمّم من أجل سطح محفوظ؛ وتراكَب. (5) واسأل ألا يغلب الحمل. (6) والتحقق المجهري: D∼ℓv∗/3.
24.5 تمارين
تمرين 24.1★
أوجد أزمنة الانتشار L2/D: من أجل مكعب سكر في قعر فنجان 5cm (D=5×10−10m2/s)؛ وعطر عبر غرفة 5m (D=5×10−6m2/s)؛ والأكسجين عبر 1mm من نسيج وعبر خلية 10µm (D=2×10−9m2/s). فأيّ هذه يفعله الانتشار فعلًا؟
حل
حل التمرين 24.1.
L2/D: السكر 5×106s (أي شهران)؛ والعطر مثله؛ والنسيج 500s؛ والخلية 0.05s. والانتشار يفعل الأخيرة فعلًا (ويفعل ما قبلها بالكاد)؛ والحمل يفعل الباقي.
تمرين 24.2★
غشاء سماكته 1µm ومساحته 1cm2، و D=1×10−11m2/s في داخله، يفصل محلولًا عند 1mol/m3 عن ماء نقي. أوجد كثافة التيار، والتدفق المولي وعدد الجزيئات في الثانية؛ والمقاومة الانتشارية؛ ومع غشاءين كهذا على التسلسل.
حل
حل التمرين 24.2.
j=DΔn/e=1×10−5mol/m2/s؛ ومضروبًا في ×1×10−4m2: 1×10−9mol/s، أي 6×1014 جزيء في الثانية؛ و Rd=e/DS=1×109s/m3؛ ومع اثنين على التسلسل: النصف.
تمرين 24.3★
طبقة رقيقة من صباغ (N=1×1020m−2) تُطلَق في الماء، و D=1×10−9m2/s. أوجد العرض 2Dt وكثافة القمة بعد 1s، و 1h، و 1d. ومتى تهبط القمة إلى 1% من قيمتها عند 1s؟
حل
حل التمرين 24.3.
2Dt: 45µm، 2.7mm، 13mm؛ والقمة N/4πDt: 9×1023m−3، 1.5×1022m−3، 3×1021m−3؛ و 1% بعد زمن أطول 104 مرة: 1×104s.
تمرين 24.4★
جزيئات الهواء: ℓ=70nm، و v∗=500m/s. قدّر D؛ وعدد التصادمات وزمن الانتشار على 1m؛ وقارن ذلك بالطيران المستقيم. وأعد ذلك لجزيء في الماء (ℓ≈0.1nm، و v∗≈500m/s): أيكون التقدير صحيحًا، ولمَ ليس بالضبط؟
حل
حل التمرين 24.4.
D≈1.2×10−5m2/s؛ و N=(L/ℓ)2=2×1014؛ و Nℓ/v∗=3×104s في وجه 2ms. وأما الماء: فالتقدير يعطي 2×10−8m2/s، أي عشرة أضعاف القيم المقيسة — ففي السائل يخشخش الجزيء في قفص جيرانه وتكون خطواته المتتالية متعاكسة الارتباط؛ فالخطوة الفعّالة أقصر.
تمرين 24.5★★
الأكسجين في نسيج. شريحة نسيج سماكتها 2a تستهلك الأكسجين بالمعدل المنتظم q (في وحدة الحجم)؛ ووجهاها محفوظان عند n0. (أ) اكتب المعادلة المستقرة مع المصبّ وحلّها. (ب) وشرط بلوغ الأكسجين المركز. (ج) ومن أجل D=2×10−9m2/s، و n0=0.2mol/m3، و q=0.01mol/m3/s: أوجد أعظم a. (د) واستنتج ما يخصّ تباعد الشعيرات.
حل
حل التمرين 24.5.
(أ) Dn′′=q: n=n0−q(a2−x2)/2D. (ب) و n(0)≥0: a≤2Dn0/q. (ج) و 0.28mm. (د) فلا تستطيع خلية أن تبعد أكثر من بضع مئات من الميكرومترات عن شعيرة (وعمليًّا 50–100µm في النسيج الفعّال).
تمرين 24.6★★
حبيبة تنحلّ. كرة سكر نصف قطرها a=1mm في ماء ساكن تحفظ سطحها عند كثافة الإشباع ns=5800mol/m3؛ والماء نقي بعيدًا؛ و D=5×10−10m2/s. (أ) أوجد المقطع المستقر والتيار الكلي. (ب) وعدد المولات في الحبيبة (والكثافة 1590kg/m3، والكتلة المولية 342g/mol) وزمن الانحلال (خذ التيار ثابتًا). (ج) ولمَ يكون الزمن الحقيقي أطول، ولمَ يفيد التقليب هذه الفائدة؟ (د) وبيّن أن المقطع المستقر من أجل أسطوانة طويلة جدًّا لوغاريتمي وأن المسألة لا حلّ لها في وسط لانهائي.
حل
حل التمرين 24.6.
(أ) n=nsa/r، و Φ=4πDans=3.6×10−8mol/s. (ب) و 1.95×10−5mol؛ و 540s. (ج) وقشرة 1/r تحتاج إلى ∼a2/D≈2000s لتتكوّن والحبيبة تتقلص؛ والتقليب يستبدل بالقشرة التي سماكتها ∼a طبقةً حدية δ≪a ويضرب التدفق في a/δ. (د) و (rn′)′=0: n=A+Blnr، وهو لا يمكن أن ينعدم في اللانهاية: فلا حالة مستقرة — إذ تظل غيمة الأسطوانة تكبر (لوغاريتميًّا).
تمرين 24.7★★
الحفظ واللاعكوسية. (أ) من المعادلة، بيّن أن ∫ndx ثابت وأن d⟨x2⟩/dt=2D (كامِل بالتجزئة، و n→0 في اللانهاية). (ب) وبيّن أن n(x,−t) لا تحقق المعادلة. (ج) وماذا يصير الحلّ الغاوسي من أجل t<0؟ (د) وبيّن أن ∫n2dx لا يمكن إلا أن يتناقص: فالانتشار يسطّح.
حل
حل التمرين 24.7.
(أ) d/dt∫n=D[n′]−∞∞=0؛ و d/dt∫x2n=D∫x2n′′=2D∫n=2DN. (ب) و ∂t[n(x,−t)]=−D∂x2n: فتنقلب الإشارة. (ج) ومن أجل t<0 يكون “العرض” 2Dt سالبًا: فلا وجود لحالة كهذه — إذ لا يمكن ردّ النبضة عن انتشارها. (د) و d/dt∫n2=2D∫nn′′=−2D∫n′2≤0.
تمرين 24.8★★
الكربنة. الكربون في الحديد الساخن، D=5×10−12m2/s عند 900∘C؛ والسطح محفوظ عند n0. (أ) تحقق من أن n0erfc(x/2Dt) يحلّ المعادلة. (ب) وأوجد العمق الذي تكون فيه n=0.1n0 بعد 4h (و erfc(1.16)=0.1). (ج) والزمن اللازم لضعف العمق. (د) وعند 950∘C يكون D أكبر بمقدار 2.5 مرة: أوجد الزمن الموفَّر.
حل
حل التمرين 24.8.
(أ) القضية 24.11. (ب) و x=2×1.16Dt=0.62mm. (ج) و 4×: 16h. (د) و 4/2.5=1.6h.
تمرين 24.9★★
ستوكس–أينشتاين.D=kBT/6πηa، وللماء η=1×10−3Pas، وعند 300K. (أ) أوجد D من أجل a=1µm و الانتقال التربيعي الوسطي في 1s و 1min. (ب) و D من أجل بروتين، a=3nm. (ج) وللسكر D=5×10−10m2/s: أوجد نصف القطر الفعّال. (د) وكيف يعطي قياس ⟨x2⟩ لحبيبة تحت المجهر عددَ أفوغادرو؟
حل
حل التمرين 24.9.
(أ) D=2.2×10−13m2/s؛ و 2Dt: 0.66µm، 5µm. (ب) و 7×10−11m2/s. (ج) و a=kBT/6πηD=0.44nm. (د) و ⟨x2⟩=2Dt=RTt/3πηaNA: فكل مقدار سوى NA مقيس.
تمرين 24.10★★★
تأخّر الغشاء. غشاء سماكته e، فارغ ابتداءً، يُعرَّض على أحد وجهيه للكثافة n1 عند t=0، والوجه الآخر محفوظ عند 0. (أ) ما التدفق المستقر النهائي؟ (ب) وبرهن أن التدفق على الوجه البعيد يرتفع خلال زمن ∼e2/D (والتأخّر المضبوط e2/6D). (ج) ولصاقة دواء مع e=20µm، و D=1×10−13m2/s عبر الطبقة الخارجية للجلد: أوجد زمن التأخّر. (د) ولمَ يفيد التأخّر في قياس D ويفيد التدفق المستقر في قياس D× الانحلالية؟
حل
حل التمرين 24.10.
(أ) Dn1S/e. (ب) وتحتاج الجسيمات إلى ∼e2/D لتعبر. (ج) و e2/6D=670s، أي إحدى عشرة دقيقة. (د) والتأخّر يعطي D وحده؛ والتدفق المستقر يعطي D× (الانحلالية): أي قياسان، ومجهولان.
تمرين 24.11★★★
الماصّ التامّ. كرة نصف قطرها a تمتصّ كل جسيم يلامسها، في وسط عند n∞ بعيدًا. (أ) أوجد المقطع المستقر والتيار الكلي الملتقَط. (ب) وجرثومة، a=1µm، في محلول سكر عند n∞=6×1020m−3، و D=5×10−10m2/s: أوجد الجزيئات الملتقَطة في الثانية. (ج) ويغطّي سطحها Nr مستقبلًا ماصًّا صغيرًا نصف قطره s، والباقي عاكس: وكل مستقبل وحده يلتقط ≈4Dsn∞ (كقرص)؛ بيّن أن الخلية كلها تلتقط قرابة العظمى متى صار Nrs≫a، أي بكسر ضئيل من سطحها مغطًّى. (د) وعلّق: لمَ تستطيع الخلايا تحمّل آلاف المستقبلات المختلفة.
حل
حل التمرين 24.11.
(أ) n=n∞(1−a/r)، و Φ=4πDan∞. (ب) و 3.8×106 في الثانية. (ج) وتعمل المستقبلات عمل مواصلات على التفرّع ثم على التسلسل مع القشرة الكروية: Φ≈4πDan∞⋅Nrs/(Nrs+πa)؛ ونصف العظمى عند Nrs=πa: فمع s=1nm، يكون Nr≈3000، مغطّيًا Nrs2/4a2≈10−3 من السطح. (د) فكل نوع من المستقبلات يحتاج إلى مساحة مهمَلة ليكون الالتقاط شبه أعظمي: فتستطيع الخلية مراقبة آلاف المواد معًا.
تمرين 24.12★★★
حلّها على شبكة. قسّم الفضاء إلى خلايا حجمها Δx والزمن إلى خطوات Δt؛ واكتب nik+1=nik+α(ni+1k−2nik+ni−1k) مع α=DΔt/Δx2. (أ) برّر ذلك من المعادلة. (ب) وفسّره مسيرًا عشوائيًّا حين α=1/2. (ج) وبيّن أنه من أجل α>1/2 تنمو كثافة تتناوب +,−,+,− من خلية إلى خلية: أي عدم استقرار. (د) ومن أجل مقطع البور في المسألة أدناه (D=1.5×10−17m2/s، و Δx=10nm)، أوجد أكبر Δt مستقر وعدد الخطوات في ساعة.
حل
حل التمرين 24.12.
(أ) فرق أمامي في t، وفرق ثانٍ مركزي في x. (ب) ومع α=1/2: nik+1=(ni−1k+ni+1k)/2 — أي إن كل جسيم يقفز يمينًا أو يسارًا باحتمال 1/2. (ج) ومن أجل ni=(−1)i، nk+1=(1−4α)nk: و ∣1−4α∣>1 حين α>1/2. (د) و Δt≤Δx2/2D=3.3s؛ ونحو 1100 خطوة.
24.6 مسألة: رقاقة مشوبة وغرفة معطَّرة
مسألة 24.1
مسألة نهاية الأسبوع — الانتشار في جسم صلب وفي غاز
الجزء الأول — دفع البور في السيليسيوم. ينتشر البور في السيليسيوم وفق D=D0e−Ea/kBT، مع D0=7.6×10−5m2/s، و Ea=3.46eV؛ و kB=8.62×10−5eV/K. وقد رُسّبت جرعة Q=1×1018m−2 من البور في طبقة رقيقة جدًّا على سطح رقاقة تشويبها الخلفي nB=1×1021m−3؛ ثم تُسخَّن (“تُدفَع إلى الداخل”) عند 1100∘C ساعةً واحدة. والسطح يعكس البور (فلا هروب).
احسب D عند 1100∘C وعند 1000∘C. وبأي عامل يتغير على هذه 100K؟
ولمَ يكون المقطع بعد الدفع n(x,t)=(Q/πDt)exp(−x2/4Dt) لا الغاوسي في القضية 24.10؟ وتحقق من أن تكامله على x>0 هو Q.
والتركيز السطحي بعد ساعة.
وعمق الوصلةxj هو حيث n=nB: احسبه.
وكيف يتغير xj إذا دام الدفع أربع ساعات؟ (وانتبه: فالتركيز السطحي يتغير أيضًا.)
والتغير النسبي في D من أجل خطأ 10K في الحرارة؛ وضبط الحرارة اللازم ليكون xj بدقة 1%.
وعند حرارة الغرفة، D: قدّر الزمن اللازم لينتقل المقطع بتباعد ذري واحد، واستنتج.
الجزء الثاني — الترسيب المسبق. وقد أُدخلت الجرعة نفسها عند 950∘C من غاز يحفظ السطح عند حدّ الانحلالية n0=2×1026m−3 مدة 30min.
أوجد D عند 950∘C و Dt من أجل 30min.
والمقطع في نهاية الترسيب المسبق؛ والعمق الذي تكون عنده n=nB (و erfc(3.2)≈6×10−6).
والجرعة المدخَلة؛ وقارنها بجرعة Q في الجزء الأول.
ولمَ يُرسَّب مسبقًا عند حرارة أدنى ويُدفَع عند حرارة أعلى؟
وبرّر أنه يمكن، في أثناء الدفع، معاملة الطبقة المرسَّبة مسبقًا كأنها رقيقة لانهائيًّا.
الجزء الثالث — عطر في غرفة. قطرة عطر (1mg، وكتلتها المولية 150g/mol) تتبخر دفعة واحدة في ركن غرفة ساكنة عند 300K و 1bar. وخذ قطر تصادم جزيء العطر مع الهواء d=0.5nmوالكثافة العددية للهواء na=P/kBT.
وفسّر في فقرة كيف تنتج الحركة العكوسة للجزيئات معادلةَ الانتشار غير العكوسة.
ولخّص: سلّم الطول، وسلّم الزمن، والمقدار الوحيد الذي يميّز غازًا وسائلًا وجسمًا صلبًا ساخنًا.
حل
حل المسألة 24.1.
1.kBT=0.118eV، و Ea/kBT=29.2: D=1.6×10−17m2/s؛ وعند 1000∘C: 1.6×10−18m2/s — أي بعامل 10.
2. يُعالَج السطح العاكس بصورة النبضة المرآوية، فتتضاعف السعة على x>0؛ و ∫0∞=(Q/πDt)⋅214πDt=Q.
3.Dt=5.8×10−14m2؛ و πDt=0.43µm؛ و ns=2.3×1024m−3.
4.e−x2/4Dt=4.3×10−4: xj2=4Dt×7.75، و xj=1.3µm.
5. ينتصف ns، ويهبط اللوغاريتم إلى 7.05، ويتضاعف 4Dt أربع مرات: xj=2.6µm — أي أقلّ من الضعف.
6.δD/D=(Ea/kBT)δT/T=21%؛ و xj∝Dt تقريبًا، فيحتاج 1% على xj إلى 2% على D: أي 1K.
7.Ea/kBT=134: D≈10−62m2/s؛ وتباعد ذري واحد (0.25nm) في a2/D∼1043s: أي مجمَّد.
8.Ea/kBT=32.8: D=4.3×10−19m2/s؛ و Dt=28nm.
9.n0erfc(x/2Dt)؛ و nB/n0=5×10−6: x≈6.4Dt=0.18µm.
10.2n0Dt/π=6.3×1018m−2: أي ستة أضعاف Q — فالجرعة تنمو وفق t ويحددها الزمن.
11. الحرارة المنخفضة: D صغير، والجرعة يضبطها الزمن، والعمق ضحل؛ والحرارة العالية: تُدفَع الجرعة الثابتة عميقًا وسريعًا.
12.0.18µm في وجه 1.3µm: أي رقيقة.
13.4×1018 جزيء؛ و na=2.4×1025m−3.
14.ℓ=37nm؛ و v∗=206m/s.
15.D≈2.5×10−6m2/s.
16.L2/D=107s، أي أشهر؛ ومع الانسياق: 50s.
17.n=8N(4πDt)−3/2e−r2/4Dt؛ وهو أعظم عند t∗=r2/6D=6.7×104s، أي ثماني عشرة ساعة.
18.nmax=8N(4πDt∗)−3/2e−3/2≈2×1018m−3، أي فوق العتبة بكثير: فيُشَمّ، في النهاية.
19.r2=4Dtln106=12m2: أي 3.5m.
20.v∗t/ℓ≈5×1014.
21.⟨x⟩=0، و ⟨x2⟩=Nℓ2=ℓ2t/τ؛ و D=ℓ2/2τ.
22.(N/2N)/2N: 1/2، 3/8، 5/16 — أي متناقص (وفق 1/πN/2): فالماشي يبتعد وفق N.
23. تتجمع الغيمة في نقطة — وهو ما لا يُرى قط؛ وهو ينتهك معادلة الانتشار، لا قوانين الميكانيك.
24. كل تصادم عكوس، لكن الحالة الابتدائية (أي الجسيمات كلها معًا) استثنائية: فكل تاريخ مجهري تقريبًا منها ينتشر، ويحتاج العكس إلى تآمر السرعات كلها. ومتوسط المسير العشوائي لا يبقي إلا ما هو نمطي ويطرح تلك المعلومة؛ فاللاعكوسية إحصائية.
25.L∼Dt، و t∼L2/D؛ و D∼ℓv∗/3، أي الخطوة بين التصادمات: 10−5، 10−9، 10−17m2/s.