Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

אוסף בעיות — קרינה תרמית וקיבול החום של מוצקים

אוסף בעיות: קרינה תרמית וקיבול החום של מוצקים

בעיות בסגנון בחינה. פתרו אותן עם עט לפני שתיגעו במחשב; בעיות 4 ו-5 נועדו להסתיים במעבדה, ובעיות 6 ו-7 הן אתגרים. הפתרונות ומחווני הבדיקה נמצאים עם חומרי המרצה ובמכוון אינם באתר הזה.

קבועים: h=6.626×10−34 J sh = 6.626 \times 10^{-34}\ \mathrm{J\,s}, kB=1.381×10−23 J K−1\kB = 1.381 \times 10^{-23}\ \mathrm{J\,K^{-1}}, c=2.998×108 m s−1c = 2.998 \times 10^{8}\ \mathrm{m\,s^{-1}}, NA=6.022×1023 mol−1N_A = 6.022 \times 10^{23}\ \mathrm{mol^{-1}}, σ=5.670×10−8 W m−2 K−4\sigma = 5.670 \times 10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}}, וכן a=4σ/c=7.566×10−16 J m−3 K−4a = 4\sigma/c = 7.566 \times 10^{-16}\ \mathrm{J\,m^{-3}\,K^{-4}}.

בעיה 1 — שני התנאים של וין

(א) החל מחוק פלאנק uν=(8πh/c3) ν3/(ehν/kBT−1)u_\nu = (8\pi h/c^3)\,\nu^3/(e^{h\nu/\kB T} - 1), הראו שהמקסימום של uνu_\nu נמצא ב-x=hν/(kBT)x = h\nu/(\kB T) המקיים x=3(1−e−x)x = 3(1 - e^{-x}).

(ב) גזרו את uλu_\lambda מתוך uνu_\nu בעזרת uλ ∣dλ∣=uν ∣dν∣u_\lambda\,|\mathrm{d}\lambda| = u_\nu\,|\mathrm{d}\nu|, והראו שהמקסימום שלה נמצא ב-y=hc/(λkBT)y = hc/(\lambda \kB T) המקיים y=5(1−e−y)y = 5(1 - e^{-y}).

(ג) פתרו את שתי המשוואות עד ארבע ספרות משמעותיות באיטרציה x←3(1−e−x)x \leftarrow 3(1 - e^{-x}) החל מ-x=3x = 3 (ובאופן דומה עבור yy). מדוע האיטרציה הזו מתכנסת, ומדוע כל כך מהר?

(ד) הראו ש-νmax⁡λmax⁡/c\nu_{\max}\lambda_{\max}/c אינו תלוי ב-TT וחשבו את ערכו. סטודנט אומר: “אחד משני החוקים של וין חייב להיות שגוי, כי לפוטון שבשיא יכול להיות רק אורך גל אחד.” ענו לסטודנט בשלושה משפטים.

(ה) ספקטרום משורטט ליחידת ln⁡ν\ln\nu, כלומר כ-ν uν\nu\,u_\nu. מצאו את התנאי למקסימום שלו ואת ה-λT\lambda T המתאים. הראו שאותה עקומה מתקבלת משרטוט של λ uλ\lambda\,u_\lambda, והסבירו מדוע הבחירה הזו נותנת שיא יחיד לשני המשתנים.

בעיה 2 — גז הפוטונים

(א) בעזרת התוצאה של מודול 04, שלפיה גז איזוטרופי מפעיל לחץ P=13n⟨p v⟩P = \tfrac13 n \langle p\,v \rangle, הראו שלקרינה מתקיים P=u/3P = u/3, ושלגז חד-אטומי לא-יחסותי מתקיים P=23uP = \tfrac23 u. איזו תכונה של הפוטון יוצרת את ההבדל?

(ב) בנפח קבוע, δQ=dU=T dS\dbar Q = \mathrm{d}U = T\,\mathrm{d}S. החל מ-u=aT4u = aT^4, גזרו את S=43aVT3S = \tfrac43 aVT^3. היכן בגזירה השתמשתם בחוק השלישי, והאם התוצאה עקבית איתו?

(ג) חלל מלא קרינה מתפשט לאט ובאופן אדיאבטי עד לנפח כפול מנפחו. באיזה גורם יורדת הטמפרטורה שלו, ובאיזה גורם יורדת האנרגיה שלו? לאן הלכה האנרגיה? השתמשו במוסכמת הסימנים של הקורס, וציינו את הסימן של העבודה.

(ד) אמתו את יחס אוילר U=TS−PV+μNU = TS - PV + \mu N עבור גז הפוטונים עם μ=0\mu = 0. אחר כך הסבירו בשני משפטים מדוע μ=0\mu = 0 נובע מכך שהדפנות בולעות ופולטות פוטונים, ולא מתכונה כלשהי של האור עצמו.

(ה) חשבו את הצפיפות המספרית של הפוטונים ב-300 K300\ \mathrm{K} מתוך n=16πζ(3)(kBT/hc)3n = 16\pi\zeta(3)(\kB T/hc)^3, עם ζ(3)=1.202\zeta(3) = 1.202, ואת האנרגיה הממוצעת לפוטון ביחידות של kBT\kB T. השוו את קיבול החום ליחידת נפח של הקרינה בחדר לזה של האוויר, ≈880 J m−3 K−1\approx 880\ \mathrm{J\,m^{-3}\,K^{-1}}.

בעיה 3 — איינשטיין ודבאי, בפיתוח לטור

(א) החל מ-U=3NℏωE/(eθE/T−1)U = 3N\hbar\omega_E/(e^{\theta_E/T} - 1), גזרו את קיבול החום של איינשטיין C=3NkB x2ex/(ex−1)2C = 3N\kB\,x^2 e^x/(e^x - 1)^2 עם x=θE/Tx = \theta_E/T.

(ב) פתחו אותו עבור xx קטן כדי להראות ש-C=3NkB(1−x2/12+O(x4))C = 3N\kB\left(1 - x^2/12 + O(x^4)\right). מה הניח מודול 01 לגבי קיבול החום של המוצק הזה, ומה גודל השגיאה של ההנחה הזו ב-T=2θET = 2\theta_E?

(ג) הראו שעבור xx גדול מתקיים C≈3NkB x2e−xC \approx 3N\kB\,x^2 e^{-x}. על פני איזה תחום טמפרטורות יורד CC בגורם 1000 אם θE=300 K\theta_E = 300\ \mathrm{K}, החל מ-T=50 KT = 50\ \mathrm{K}? השוו לחוק T3T^3 על פני אותו גורם.

(ד) במודל דבאי, C/(3NkB)=⟨x2ex/(ex−1)2⟩C/(3N\kB) = \langle x^2 e^x/(e^x - 1)^2 \rangle: פונקציית איינשטיין ממוצעת על פני האופנים במשקל 3ω2/ωD33\omega^2/\omega_D^3 בתחום 0<ω<ωD0 < \omega < \omega_D. השתמשו בסעיף (ב) כדי להראות שבטמפרטורה גבוהה C=3NkB(1−xD2/20+… )C = 3N\kB\left(1 - x_D^2/20 + \dots\right), כאשר xD=θD/Tx_D = \theta_D/T. מכאן מצאו את היחס θE/θD\theta_E/\theta_D שבו שני המודלים מסכימים עד הסדר הזה.

(ה) הראו שמודל דבאי נותן C∝T3C \propto T^3 בטמפרטורה נמוכה, בעזרת ∫0∞x4ex/(ex−1)2 dx=4π4/15\int_0^\infty x^4 e^x/(e^x - 1)^2\,\mathrm{d}x = 4\pi^4/15, ומצאו את המקדם. הסבירו פיזיקלית מדוע שום בחירה של θE\theta_E אינה מאפשרת למודל איינשטיין לעשות אותו דבר.

בעיה 4 — הכסף, כפי שנמדד

המעבדה טוענת את קיבול החום של כסף מ-1 עד 300 K300\ \mathrm{K}, כפי שהוא מופיע בטבלאות של NIST (סמית’ ופיקט, 1995), ביחידות J kg−1 K−1\mathrm{J\,kg^{-1}\,K^{-1}}; המסה המולרית של כסף היא 0.1078682 kg mol−10.1078682\ \mathrm{kg\,mol^{-1}}.

(א) שרטטו את C/TC/T כנגד T2T^2 עבור T≤6 KT \le 6\ \mathrm{K} והתאימו קו ישר. פרשו את נקודת החיתוך ואת השיפוע, וחלצו את המקדם האלקטרוני γ\gamma ואת טמפרטורת דבאי θD\theta_D, כל אחד עם אי-ודאות מההתאמה.

(ב) החסירו את γT\gamma T מכל הטבלה. מדדו את השיפוע בגרף לוג-לוג של מה שנותר בין 1 ל-5 K5\ \mathrm{K}, ושל הנתונים הגולמיים בין 1 ל-3 K3\ \mathrm{K}. הסבירו את ההבדל.

(ג) התאימו את מודל דבאי ואת מודל איינשטיין, פרמטר אחד לכל אחד, לקיבול החום של הסריג בין 4 ל-100 K100\ \mathrm{K}, במזעור סכום הריבועים של השאריות היחסיות. דווחו על כל טמפרטורה מותאמת ועל שורש הממוצע של ריבועי השאריות. מדוע שאריות יחסיות ולא מוחלטות?

(ד) לכל נקודת נתונים, פתרו עבור ה-θD\theta_D שנוסחת דבאי צריכה כדי לשחזר אותה, ושרטטו את ה-θD(T)\theta_D(T) ה"רץ" הזה מ-2 עד 120 K120\ \mathrm{K}. מה הסחיפה מלמדת אתכם על מודל דבאי? מדוע ההליך נכשל מעל כ-150 K150\ \mathrm{K}?

(ה) מעל כ-260 K260\ \mathrm{K} ה-CPC_P הנמדד עולה על 3R3R. ציינו שתי סיבות והעריכו איזו מהן גדולה יותר ב-300 K300\ \mathrm{K}, בהינתן שלכסף α=5.7×10−5 K−1\alpha = 5.7 \times 10^{-5}\ \mathrm{K^{-1}} (מקדם ההתפשטות הנפחית), κT=1.0×10−11 Pa−1\kappa_T = 1.0 \times 10^{-11}\ \mathrm{Pa^{-1}} ונפח מולרי 1.03×10−5 m3 mol−11.03 \times 10^{-5}\ \mathrm{m^3\,mol^{-1}}.

בעיה 5 — קרינת הרקע הקוסמית

המעבדה טוענת את ספקטרום המונופול של COBE/FIRAS: 43 תדירויות, הנתונות כמספרי גל ν~\tilde\nu ביחידות cm−1\mathrm{cm^{-1}} (כך ש-ν=100 c ν~\nu = 100\,c\,\tilde\nu), עם העוצמה הסגולית ביחידות MJy sr−1\mathrm{MJy\,sr^{-1}} (1 MJy sr−1=10−20 W m−2 Hz−1 sr−11\ \mathrm{MJy\,sr^{-1}} = 10^{-20}\ \mathrm{W\,m^{-2}\,Hz^{-1}\,sr^{-1}}) ואת אי-הוודאות שלה ביחידות kJy sr−1\mathrm{kJy\,sr^{-1}}, יחידות קטנות פי אלף.

(א) התאימו את הטמפרטורה של ספקטרום הרדיאנס של פלאנק Bν=(c/4π) uνB_\nu = (c/4\pi)\,u_\nu במזעור של χ2\chi^2. דווחו על TT, על אי-הוודאות הסטטיסטית שלו ועל χ2\chi^2 לדרגת חופש. מה בודק המספר האחרון שהראשון אינו בודק?

(ב) קראו את השיא של הנתונים עצמם, בתדירות וביחידות cm−1\mathrm{cm^{-1}}. המירו את הספקטרום לצפיפות ליחידת אורך גל ומצאו היכן היא מגיעה לשיא. הסבירו את הגורם בין שני אורכי הגל.

(ג) קובץ הנתונים מציין שעמודת הספקטרום שלו נבנתה כעקומת פלאנק של 2.725 K2.725\ \mathrm{K} ועוד השאריות הנמדדות. מה משתמע מכך לגבי מה שההתאמה שלכם בסעיף (א) יכולה ואינה יכולה לבסס? מה מגביל את הטמפרטורה המוחלטת שפורסמה?

(ד) חשבו את הצפיפות המספרית של פוטוני הרקע ואת צפיפות האנרגיה שלהם, והשוו את האחרונה לצפיפות האנרגיה של אור הכוכבים בגלקסיה שלנו, כ-10−13 J m−310^{-13}\ \mathrm{J\,m^{-3}}.

(ה) הטמפרטורה של קרינת הרקע יורדת ביחס הפוך לגודל היקום. בעזרת בעיה 2(ג), הסבירו מדוע הספקטרום נשאר ספקטרום פלאנק מדויק בזמן שהוא מתקרר, ולא רק מוסט.

בעיה 6 — מאזן האנרגיה של כוכב לכת (אתגר)

כדור הארץ מקבל S=1361 W m−2S = 1361\ \mathrm{W\,m^{-2}} מהשמש ומחזיר חלק A=0.30A = 0.30 ממנו.

(א) הסבירו את הגורם 4 ב-σTe4=(1−A)S/4\sigma T_e^4 = (1 - A)S/4, וחשבו את TeT_e.

(ב) הוסיפו שכבת אטמוספרה אחת, השקופה לאור השמש ובולעת את כל הקרינה התרמית. כתבו את מאזן האנרגיה של השכבה ושל פני השטח, והראו ש-Ts=21/4TeT_s = 2^{1/4}T_e. הכלילו ל-nn שכבות.

(ג) כעת תנו לשכבה לבלוע רק חלק ε\varepsilon מהקרינה התרמית העוברת דרכה (ולפי חוק קירכהוף, לפלוט εσTa4\varepsilon\sigma T_a^4 מכל אחד מצדדיה). הראו ש-Ts=Te [2/(2−ε)]1/4T_s = T_e\,[2/(2 - \varepsilon)]^{1/4}, ומצאו את ה-ε\varepsilon שנותן את ה-288 K288\ \mathrm{K} הנצפים.

(ד) בעזרת חוק ההסחה של וין, הסבירו מדוע סביר להתייחס לשכבה כשקופה לאור השמש אך בולעת את הקרינה של כוכב הלכת עצמו, וציינו דרך אחת שבה אטמוספרה אמיתית אינה שכבה יחידה בטמפרטורה יחידה.

בעיה 7 — טמפרטורת דבאי ממהירות הקול (אתגר)

מודל דבאי נותן את θD\theta_D במונחי הצפיפות המספרית n=N/Vn = N/V ומהירות קול ממוצעת, ωD3=6π2n cs3\omega_D^3 = 6\pi^2 n\,c_s^3, θD=ℏωD/kB\theta_D = \hbar\omega_D/\kB.

(א) הראו שאם המהירות האורכית vLv_L ושתי המהירויות הרוחביות vTv_T שונות זו מזו, צפיפות האופנים מחייבת את הממוצע 3/cs3=1/vL3+2/vT33/c_s^3 = 1/v_L^3 + 2/v_T^3.

(ב) ליהלום צפיפות 3515 kg m−33515\ \mathrm{kg\,m^{-3}} ומסה מולרית 12.011 g mol−112.011\ \mathrm{g\,mol^{-1}}; קחו vL=18.0 km s−1v_L = 18.0\ \mathrm{km\,s^{-1}} ו-vT=12.4 km s−1v_T = 12.4\ \mathrm{km\,s^{-1}}. חשבו את θD\theta_D והשוו לערך ≈2230 K\approx 2230\ \mathrm{K} המתקבל מקיבול החום בטמפרטורה נמוכה.

(ג) חזרו על החישוב עבור כסף: צפיפות 10 490 kg m−310\,490\ \mathrm{kg\,m^{-3}}, מסה מולרית 107.87 g mol−1107.87\ \mathrm{g\,mol^{-1}}, vL=3.65 km s−1v_L = 3.65\ \mathrm{km\,s^{-1}}, vT=1.66 km s−1v_T = 1.66\ \mathrm{km\,s^{-1}}. השוו ל-θD\theta_D שמצאתם בבעיה 4(א).

(ד) החישוב בסעיפים (ב) ו-(ג) השתמש רק בתכונות אלסטיות — סרגל, מאזניים ושעון עצר — ובכל זאת ניבא קיבול חום. הסבירו אילו אופנים מובטח שהוא מתאר נכון, ומדוע ההתאמה טובה ביותר עבור ה-θD\theta_D של הטמפרטורה הנמוכה.

(ה) בעזרת בעיה 3(ד), העריכו את טמפרטורת איינשטיין של יהלום שתתאים למודל דבאי בטמפרטורה גבוהה. השוו אותה ל-θE\theta_E שמשחזר את קיבול החום הנמדד של יהלום בטמפרטורת החדר, 6.12 J mol−1 K−16.12\ \mathrm{J\,mol^{-1}\,K^{-1}}, כ-1310 K1310\ \mathrm{K}. מה אומר הפער על מודל בעל תדירות אחת?