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.

קרינה תרמית וקיבול החום של מוצקים

החידה: שתי קטסטרופות, הנחה אחת

כוונו תנור מטבח ל-500 K500\ \mathrm{K} והביטו פנימה. חשוך בו. הדפנות שלו קורנות, ואפשר להרגיש זאת על הפנים, אבל שום דבר אינו נראה: הקרינה היא אינפרה-אדומה, והשיא שלה סביב 6 μm6\ \mu\mathrm{m}.

עכשיו הפעילו את הפיזיקה של מודולים 04 עד 12, ביושר. פנים התנור הוא קופסה מלאה בגלים אלקטרומגנטיים. כל גל עומד שנכנס בקופסה הוא מתנד, וחלוקת האנרגיה השווה מעניקה לכל מתנד kBT\kB T. אין גל קצר ביותר שנכנס, ולכן אין מתנד אחרון: מספר האופנים ליחידת תדירות גדל כריבוע התדירות, והאנרגיה גדלה איתו. ב-500 K500\ \mathrm{K} התורה הקלאסית מכניסה לתחום האור הנראה לבדו יותר אנרגיה ממה שהתנור כולו מחזיק בפועל — פי תשעה-עשר אלף — ועוד יותר לאולטרה-סגול, ללא גבול. כל תנור חם היה צריך להיות מנורה אולטרה-סגולה בעלת הספק אינסופי.

לדפנות הפלדה של התנור יש בעיה שנייה. מוצק של NN אטומים הוא 3N3N מתנדים, אחד לכל כיוון תנודה של כל אטום. כל אחד מהם, שוב, מחזיק kBT\kB T, ולכן קיבול החום הוא 3NkB3N\kB — ערך דולונג–פטי, בערך 25 J mol−1 K−125\ \mathrm{J\,mol^{-1}\,K^{-1}} — בכל טמפרטורה. מדדו אותו. כסף בטמפרטורת החדר קרוב לכך. לכסף ב-4 K4\ \mathrm{K} יש קיבול חום קטן כמעט פי אלף ותשע מאות. קיבול החום של יהלום הוא רבע מערך דולונג–פטי בטמפרטורת החדר.

מודול 12 כבר פגש את התשובה עבור מתנד יחיד, וקרא לה קפיאה. המודול הזה לוקח את אותו רעיון לכל אופן של השדה האלקטרומגנטי ולכל אופן של גביש, ומוצא שהוא מציל את שניהם. לאורך הדרך הוא פוגש חידה קטנה יותר שמכשילה כמעט את כולם: לשאלה “באיזה אורך גל השמש מאירה בעוצמה הרבה ביותר?” יש שתי תשובות נכונות, ואיזו מהן מקבלים תלוי באופן שבו משרטטים את הגרף.

לנבא לפני שמחשבים

התחייבו לתשובה לכל אחת מן השאלות האלה לפני שתמשיכו לקרוא. כתבו אותן.

  1. תנור מחומם מ-300 K300\ \mathrm{K} ל-600 K600\ \mathrm{K}. האם ההספק הכולל שהדפנות שלו מקרינות גדל בגורם 2, 4, 8 או 16?

  2. בשרטוט ליחידת אורך גל, הספקטרום של השמש מגיע לשיא סביב 500 nm500\ \mathrm{nm}. בשרטוט ליחידת תדירות, האם השיא שלו נמצא בתדירות של אור באורך גל 500 nm500\ \mathrm{nm}, c/500 nm=600 THzc/500\ \mathrm{nm} = 600\ \mathrm{THz}?

  3. יהלום מקורר מ-300 K300\ \mathrm{K} ל-30 K30\ \mathrm{K}. האם קיבול החום שלו יורד בגורם של בערך 10, או בכמה סדרי גודל? שרטטו את C(T)C(T).

  4. הטמפרטורה של חלל מוכפלת. האם מספר הפוטונים שבתוכו מוכפל?

לחקור את המודל

במעבדה יש שני כלי חקירה. הראשון משרטט את ספקטרום הגוף השחור בטמפרטורה שאתם בוחרים, ליחידת תדירות וליחידת אורך גל זה לצד זה, עם הניבוי הקלאסי לפי דרישה, ומסמן היכן כל עקומה מגיעה לשיא — והיכן כל שיא נוחת כשמעבירים אותו לגרף השני. השני משרטט את קיבול החום של מודלי איינשטיין ודבאי מול נתונים מדודים של כסף מ-1 K1\ \mathrm{K} עד 300 K300\ \mathrm{K}, עם מחוון לפרמטר היחיד של כל מודל. אחר כך המעבדה מריצה את הניסויים שהמודול הזה בנוי עליהם: הנוסחה הקוונטית נבנית מחדש מהסכום של מודול 12, האנרגיה הכוללת מחושבת באינטגרציה נומרית, שני השיאים מאותרים בחיפוש ממצה, שני המודלים מותאמים לכסף, וקרינת הרקע הקוסמית עוברת התאמה.

ספקטרום הגוף השחור בזמן שהטמפרטורה מטפסת מ-3000 K3000\ \mathrm{K} ל-8000 K8000\ \mathrm{K}, עם עצירה ב-5772 K5772\ \mathrm{K} של השמש. משמאל: אנרגיה ליחידת נפח ליחידת תדירות, uνu_\nu, כנגד תדירות (כחול). מימין: אנרגיה ליחידת נפח ליחידת אורך גל, uλu_\lambda, כנגד אורך גל (אדום). התחום הצבוע הוא אור נראה, מ-380 עד 750 nm750\ \mathrm{nm}. כל עקומה משורטטת מנורמלת לגובה השיא שלה, כדי שאפשר יהיה להשוות צורות ומיקומים; האנרגיה שמתחת לעקומות האמיתיות גדלה פי חמישים לאורך הסריקה, כמו T4T^4. הקו הרציף בכל לוח מסמן את השיא שלו; הקו המקווקו מסמן את השיא של הלוח השני, מועבר דרך λ=c/ν\lambda = c/\nu. הקו המקווקו לעולם אינו נוחת על הרציף: בכל טמפרטורה עקומת התדירות מגיעה לשיא בפוטון שאורך הגל שלו גדול פי 1.76 מאורך הגל שבו עקומת אורך הגל מגיעה לשיא — היחס νmax⁡λmax⁡/c=0.568\nu_{\max}\lambda_{\max}/c = 0.568 המודפס מעל הלוחות. שתי העקומות מתארות את אותה קרינה, והאנרגיה שמתחתן זהה.

הריצו את המעבדה עכשיו, ואז חזרו. ארבעה דברים כדאי לעשות לפני שתמשיכו לקרוא:

לגזור את התוצאה

מערכת וגבול. שתי מערכות עם מבנה אחד. חלל בנפח VV שדפנותיו מוחזקות בטמפרטורה TT, המכיל קרינה אלקטרומגנטית; וגביש של NN אטומים בטמפרטורה TT. שתיהן אוספים של אופנים הרמוניים בלתי תלויים, שכל אחד מהם במגע עם אמבט חום — הדפנות עבור הקרינה, שאר הגביש עבור כל אחת מהתנודות שלו.

משתנים בלתי תלויים. הטמפרטורה TT והנפח VV. כל השאר — האנרגיה, הספקטרום, לחץ הקרינה, מספר הפוטונים — הוא פונקציה של שני אלה.

מוסכמת הסימן. מוסכמת הקורס, dU=δQ+δWon\mathrm{d}U = \dbar Q + \dbar \Won, תקפה. היא מופעלת פעם אחת, בשלב 5: דחיסת קרינה דורשת עבודה, δWon=−P dV>0\dbar \Won = -P\,\mathrm{d}V > 0 כאשר dV<0\mathrm{d}V < 0, מפני שלקרינה יש לחץ.

יחידות. SI לאורך כל הדרך, עם kB\kB, hh ו-cc כתובים במפורש; ℏ=h/2π\hbar = h/2\pi. לאופן בתדירות ν\nu (בהרץ) יש תדירות זוויתית ω=2πν\omega = 2\pi\nu וקוונט hν=ℏωh\nu = \hbar\omega. ספקטרומים הם אנרגיות ליחידת נפח ליחידת תדירות, ביחידות J s m−3\mathrm{J\,s\,m^{-3}}, או ליחידת אורך גל, ביחידות J m−4\mathrm{J\,m^{-4}}.

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

שלב 1: אופן אחד, נוסחה אחת

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

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

Z=∑n=0∞e−nhν/(kBT)=11−e−hν/(kBT),Z = \sum_{n=0}^{\infty} e^{-n h\nu/(\kB T)} = \frac{1}{1 - e^{-h\nu/(\kB T)}} ,

והאנרגיה הממוצעת שלו היא מינוס הנגזרת של ln⁡Z\ln Z לפי β=1/(kBT)\beta = 1/(\kB T).

שני הגבולות חשובים.

הגבול השני הזה הוא התרופה, וכדאי לומר מה הוא אינו. אופן שעבר קפיאה אינו דומם. במצב היסוד שלו מתנד שומר על תנועת נקודת האפס שלו; מה שהוא אינו יכול לעשות הוא לקבל אנרגיה בכמויות קטנות מ-hνh\nu. קפיאה היא מחיר שאי אפשר לשלם, לא דממה.

שלב 2: ספירת האופנים של חלל

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

k=πL (nx,ny,nz),nx,ny,nz=1,2,3,…\mathbf{k} = \frac{\pi}{L}\,(n_x, n_y, n_z) , \qquad n_x, n_y, n_z = 1, 2, 3, \ldots

והתדירות שלו היא ν=c ∣k∣/(2π)=c ∣n∣/(2L)\nu = c\,|\mathbf{k}|/(2\pi) = c\,|\mathbf{n}|/(2L). כל k\mathbf{k} מותר נושא שני קיטובים בלתי תלויים. ספירת אופנים היא ספירת הנקודות של סריג של מספרים שלמים. האופנים שהתדירות שלהם מתחת ל-ν\nu הם נקודות הסריג שבתוך כדור ברדיוס ∣n∣=2Lν/c|\mathbf{n}| = 2L\nu/c, המוגבלות לאוקטנט שבו שלושת המספרים השלמים חיוביים. כשהרדיוס גדול יש נקודת סריג אחת ליחידת נפח, והספירה היא נפח האוקטנט:

N(ν)=2×18×4π3(2Lνc)3=8πV3c3 ν3.N(\nu) = 2 \times \frac{1}{8} \times \frac{4\pi}{3}\left(\frac{2L\nu}{c}\right)^3 = \frac{8\pi V}{3c^3}\,\nu^3 .

הקירוב שב"לספור את הנפח במקום את הנקודות" ניתן לבדיקה, והמעבדה בודקת אותו: הספירה מחמיצה את נקודות הסריג שעל שלוש הפאות השטוחות של האוקטנט — גלים שאחד האינדקסים שלהם שווה לאפס, שקיימים אבל נושאים קיטוב אחד בלבד — והשגיאה היא יחס של שטח פנים לנפח, בגודל יחסי של בערך (9/4)/∣n∣(9/4)/|\mathbf{n}|. לחלל של 1 cm1\ \mathrm{cm} בטמפרטורת החדר, לאופנים החשובים יש ∣n∣|\mathbf{n}| מסדר גודל של אלף.

עכשיו תנו לכל אופן את ה-kBT\kB T הקלאסי. האנרגיה ליחידת נפח ליחידת תדירות היא

uνRJ=g(ν)V kBT=8πν2c3 kBT,u_\nu^{\mathrm{RJ}} = \frac{g(\nu)}{V}\,\kB T = \frac{8\pi \nu^2}{c^3}\,\kB T ,

חוק ריילי–ג’ינס. הוא מתאר היטב חללים ממשיים בתדירות נמוכה. אבל הוא גדל כמו ν2\nu^2 ללא גבול, והאנרגיה הכוללת היא האינטגרל שלו על פני כל התדירויות.

הספקטרום הקלאסי (אדום, מקווקו) מול הקוונטי (כחול) בזמן שהטמפרטורה עולה מ-1000 K1000\ \mathrm{K} ל-6000 K6000\ \mathrm{K}, עם תחום האור הנראה צבוע; שתי העקומות מנורמלות כך שלשיא של העקומה הקוונטית יש אותו גובה בכל פריים. השתיים מסכימות בתדירות נמוכה, שם קוונט זול וכל אופן מחזיק קרוב ל-kBT\kB T. מעל השיא העקומה הקוונטית מתהפכת ויורדת, מפני שאופנים שהקוונט שלהם עולה על כמה kBT\kB T עוברים קפיאה; העקומה הקלאסית אינה מתהפכת, והשטח שמתחתיה אינסופי. התדירות שבה דרכיהן נפרדות עולה ביחס ישר לטמפרטורה.

שלב 3: חוק פלאנק, בשני משתנים

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

בתדירות נמוכה הוא הופך לריילי–ג’ינס, שם התורה הקלאסית צדקה. בתדירות גבוהה הוא יורד כמו ν3e−hν/(kBT)\nu^3 e^{-h\nu/(\kB T)}, שם התורה הקלאסית טעתה. פלאנק מצא את הנוסחה ב-1900 בהתאמה לספקטרומים מדודים, ואז מצא את הקוונטיזציה שמייצרת אותה; זוהי התוצאה העתיקה ביותר של הפיזיקה הקוונטית.

ספקטרומים משורטטים גם כנגד אורך גל. המעבר מאחד לשני נראה כמו הצבה, והוא אינו כזה, מפני ש-uνu_\nu אינה פונקציה של הקרינה לבדה. היא צפיפות: האנרגיה בפס צר היא uν dνu_\nu\,\mathrm{d}\nu, ורק המכפלה הזו פיזיקלית.

בהחלה על חוק פלאנק,

uλ(λ,T)=8πhcλ5 1ehc/(λkBT)−1.u_\lambda(\lambda, T) = \frac{8\pi h c}{\lambda^5}\,\frac{1}{e^{hc/(\lambda \kB T)} - 1} .

עכשיו מצאו את שני השיאים. ליחידת תדירות, כתבו x=hν/(kBT)x = h\nu/(\kB T); אז uν∝x3/(ex−1)u_\nu \propto x^3/(e^x - 1), והשוואת הנגזרת שלה לאפס נותנת

x=3(1−e−x),x=2.821,x = 3\left(1 - e^{-x}\right) , \qquad x = 2.821 ,

ולכן

νmax⁡=2.821 kBTh.\nu_{\max} = 2.821\,\frac{\kB T}{h} .

ליחידת אורך גל, כתבו y=hc/(λkBT)y = hc/(\lambda \kB T); אז uλ∝y5/(ey−1)u_\lambda \propto y^5/(e^y - 1) — שתי חזקות נוספות, מהיעקוביאן — ואותו חישוב נותן

y=5(1−e−y),y=4.965,y = 5\left(1 - e^{-y}\right) , \qquad y = 4.965 ,

ולכן

λmax⁡T=hc4.965 kB=2.898×10−3 m K.\lambda_{\max} T = \frac{hc}{4.965\,\kB} = 2.898 \times 10^{-3}\ \mathrm{m\,K} .

עבור השמש, ב-5772 K5772\ \mathrm{K}: λmax⁡=502 nm\lambda_{\max} = 502\ \mathrm{nm}, אור ירוק. אבל νmax⁡=339 THz\nu_{\max} = 339\ \mathrm{THz}, התדירות של אור באורך גל 883 nm883\ \mathrm{nm} — באינפרה-אדום, בלתי נראה. שתי הטענות נכונות. אף אחת מהן אינה “המקום שבו האור של השמש באמת מגיע לשיא”, מפני שאין מקום כזה: השיא שייך לגרף, לא לאור. היעקוביאן c/λ2c/\lambda^2 גדול ביותר באורכי גל קצרים, שם ננומטר אחד משתרע על פני טרה-הרצים רבים, ולכן ליחידת אורך גל הקצה הכחול של הספקטרום מקבל משקל רב יותר והשיא זז לעברו. בחירה שלישית — ליחידת ln⁡ν\ln \nu, יחסי תדירות שווים, שהיא גם ליחידת ln⁡λ\ln \lambda — מגיעה לשיא בין השתיים האחרות, ב-x=3.921x = 3.921, 636 nm636\ \mathrm{nm} עבור השמש.

שלב 4: חוק סטפן–בולצמן

צפיפות האנרגיה הכוללת אינה תלויה בבחירת המשתנה. באינטגרציה של חוק פלאנק על פני התדירות עם x=hν/(kBT)x = h\nu/(\kB T),

u=∫0∞uν dν=8π(kBT)4h3c3∫0∞x3ex−1 dx.u = \int_0^\infty u_\nu\,\mathrm{d}\nu = \frac{8\pi (\kB T)^4}{h^3 c^3} \int_0^\infty \frac{x^3}{e^x - 1}\,\mathrm{d}x .

האינטגרל שנותר הוא מספר טהור. הוא שווה ל-π4/15=6.494\pi^4/15 = 6.494 (הוא 6 ζ(4)6\,\zeta(4), כאשר ζ\zeta היא פונקציית זטא של רימן); המודול הזה מצטט את הערך, והמעבדה בודקת אותו נומרית.

הגורם c/4c/4 בא מספירה של התורה הקינטית של מה שחוצה חור: מחצית הפוטונים הסמוכים אליו נעים לעבר הדופן, ורכיב המהירות הממוצע שלהם בניצב לה הוא c/2c/2, ולכן השטף הוא רבע מהצפיפות כפול המהירות. סטפן מצא את החזקה הרביעית ב-1879 מקריאת מדידות של אחרים, כחוק אמפירי; בולצמן גזר אותה מהתרמודינמיקה ב-1884, אבל לא יכול היה לקבל את aa, שדורש את hh. החזקה הרביעית פירושה שהטמפרטורה היא מנוף רב-עוצמה: הכפלתה מכפילה את ההספק המוקרן פי שש-עשרה.

מכשירים המכוונים אל מקור מודדים את הרדיאנס הספקטרלי, ההספק ליחידת שטח, ליחידת תדירות, ליחידת זווית מרחבית, Bν=(c/4π) uνB_\nu = (c/4\pi)\,u_\nu, ביחידות W m−2 Hz−1 sr−1\mathrm{W\,m^{-2}\,Hz^{-1}\,sr^{-1}}. יש לו אותה צורה כמו ל-uνu_\nu.

מספר הפוטונים מתקבל באותה דרך, עם חזקה אחת פחות של xx:

n=∫0∞g(ν)V ⟨n⟩ dν=8π(kBThc)3∫0∞x2ex−1 dx=16π ζ(3)(kBThc)3,n = \int_0^\infty \frac{g(\nu)}{V}\,\langle n \rangle\,\mathrm{d}\nu = 8\pi\left(\frac{\kB T}{hc}\right)^3 \int_0^\infty \frac{x^2}{e^x - 1}\,\mathrm{d}x = 16\pi\,\zeta(3)\left(\frac{\kB T}{hc}\right)^3 ,

עם ζ(3)=1.202\zeta(3) = 1.202: בערך 2.03×107 T32.03 \times 10^{7}\,T^3 פוטונים למטר מעוקב, כאשר TT בקלווין. מספר הפוטונים גדל כמו T3T^3. פוטונים אינם נשמרים; חימום חלל יוצר אותם.

שלב 5: גז הפוטונים

קרינה בחלל היא גז — של פוטונים — ויש לה לחץ, אנטרופיה ופוטנציאל כימי. כל אחד מהם דורש שורה.

לחץ. מודול 04 מצא את הלחץ של גז כשטף תנע: לגז איזוטרופי, P=13 n ⟨p v⟩P = \tfrac13\,n\,\langle p\,v \rangle, שליש מהצפיפות של המכפלה של תנע ומהירות. לחלקיק לא-יחסותי pv=mv2p v = m v^2, פעמיים האנרגיה הקינטית, ו-P=23uP = \tfrac23 u. לפוטון יש v=cv = c ותנע p=E/cp = E/c, ולכן pv=Ep v = E, ומכאן

P=u3=aT43.P = \frac{u}{3} = \frac{a T^4}{3} .

הוא קטן בטמפרטורות רגילות — 2 μPa2\ \mu\mathrm{Pa} ב-300 K300\ \mathrm{K} — והוא מבצע עבודה: דחיסת החלל דורשת δWon=−P dV\dbar \Won = -P\,\mathrm{d}V, חיובית כשהנפח מתכווץ. זהו לחץ הקרינה, הלחץ שמחזיק את הכוכבים המסיביים ביותר.

אנטרופיה. בנפח קבוע לא נעשית עבודה, ולכן δQ=dU=T dS\dbar Q = \mathrm{d}U = T\,\mathrm{d}S (מודול 09). קיבול החום הוא CV=dU/dT=4aVT3C_V = \mathrm{d}U/\mathrm{d}T = 4aVT^3, ובאינטגרציה של CV/TC_V/T מהאפס המוחלט כלפי מעלה,

S=∫0T4aVT′3T′ dT′=43 aVT3.S = \int_0^T \frac{4aVT'^3}{T'}\,\mathrm{d}T' = \frac43\,a V T^3 .

בהתפשטות אדיאבטית איטית SS נשאר קבוע, ולכן VT3VT^3 קבוע ו-T∝V−1/3T \propto V^{-1/3}: קרינה מתקררת ביחס הפוך לגודל המכל שלה. זו הסיבה שהקרינה שנותרה מהיקום המוקדם נמצאת עכשיו ב-2.7 K2.7\ \mathrm{K}.

פוטנציאל כימי.

שלב 6: מוצק איינשטיין, נמדד סוף סוף

מודול 01 נתן לשני מוצקי איינשטיין להחליף קוונטים; מודול 09 בנה את האנטרופיה של מוצק איינשטיין; מודול 12 מצא את הטמפרטורה שלו מתוך פונקציית החלוקה של המתנד. המודל תמיד היה אותו מודל: NN אטומים, שכל אחד מהם מתנודד בשלושה כיוונים, כ-3N3N מתנדים בלתי תלויים החולקים תדירות אחת. כאן הוא סוף סוף נותן ניבוי שקלורימטר יכול לבדוק.

עם x=θE/Tx = \theta_E/T, שלב 1 נותן לכל מתנד ℏωE/(ex−1)\hbar\omega_E/(e^x - 1) מעל מצב היסוד שלו, ולכן

U=3NℏωEeθE/T−1,U = \frac{3N\hbar\omega_E}{e^{\theta_E/T} - 1} ,

ובגזירה לפי TT:

זהו קיבול החום של המתנד ממודול 12 כפול 3N3N; המעבדה בודקת ששתיהן אותה פונקציה עד הספרה האחרונה. הגבולות שלו:

פיתוח הטמפרטורה הגבוהה גם פורע חוב. מודול 01 השתמש בקיבול חום קבוע, kB\kB אחד למתנד, והזהיר שזהו קירוב. השגיאה ידועה עכשיו: שבר (θE/T)2/12(\theta_E/T)^2/12 כאשר TT גדול, 8%8\% ב-T=θET = \theta_E, 28%28\% ב-T=θE/2T = \theta_E/2. התמונה של CC קבוע בטוחה מעל בערך θE\theta_E, ונכשלת, בפער מעריכי, הרבה מתחתיה.

וגבול הטמפרטורה הנמוכה נכשל במבחן אחר. קיבולי חום מדודים של גבישים מבדדים אכן מתאפסים בטמפרטורה נמוכה, אבל כמו T3T^3, לא מעריכית. מעריכי יורד מהר בהרבה מכל חזקה: מודל איינשטיין מנבא קיבול חום קטן בהרבה מדי בקור. האשמה היא בתדירות היחידה שלו. גביש ממשי גם מתנודד לאט — גלי קול באורך גל ארוך — ואופן בתדירות נמוכה מספיק לעולם אינו עובר קפיאה.

שלב 7: התיקון של דבאי

התנודות של גביש הן גלים. באורכי גל ארוכים הן קול, הנע במהירות csc_s, עם שלושה קיטובים — שניים רוחביים, אחד אורכי — במקום השניים של האור. הקוונטים שלהן נקראים פונונים, וכמו פוטונים מספרם אינו נשמר. המודל של דבאי סופר אותם בדיוק כפי ששלב 2 ספר פוטונים, עם csc_s במקום cc ושלושה קיטובים במקום שניים. בתדירות זוויתית,

g(ω) dω=3V2π2cs3 ω2 dω.g(\omega)\,\mathrm{d}\omega = \frac{3V}{2\pi^2 c_s^3}\,\omega^2\,\mathrm{d}\omega .

לגביש של NN אטומים יש רק 3N3N אופנים — גל קצר מפעמיים המרווח הבין-אטומי אינו תנועה חדשה — ולכן דבאי קוטע את הספקטרום בתדירות ωD\omega_D שבה הספירה מגיעה ל-3N3N:

∫0ωDg(ω) dω=3N⟹ωD3=6π2 NV cs3,θD=ℏωDkB.\int_0^{\omega_D} g(\omega)\,\mathrm{d}\omega = 3N \qquad\Longrightarrow\qquad \omega_D^3 = 6\pi^2\,\frac{N}{V}\,c_s^3 , \qquad \theta_D = \frac{\hbar\omega_D}{\kB} .

בסכימת שלב 1 על פני הספקטרום הזה ובגזירה, עם x=ℏω/(kBT)x = \hbar\omega/(\kB T):

ל-T3T^3 יש הסבר ישיר. ב-TT נמוך רק אופנים עם ℏω≲kBT\hbar\omega \lesssim \kB T פעילים; יש בערך (kBT/ℏcs)3V(\kB T/\hbar c_s)^3 V כאלה, כמו לפוטונים; כל אחד מחזיק בערך kBT\kB T. לכן U∝T4U \propto T^4 ו-C∝T3C \propto T^3 — חוק סטפן–בולצמן של הקול.

מודל דבאי (כחול) ומודל איינשטיין (אדום) מול קיבול החום המדוד של כסף (שחור, אחרי הסרת ה-γT\gamma T של האלקטרונים), חלקי 3R3R. טמפרטורת דבאי נסרקת מ-120 עד 320 K320\ \mathrm{K} ומתייצבת על 218 K218\ \mathrm{K}; טמפרטורת איינשטיין קשורה אליה דרך θE=3/5 θD\theta_E = \sqrt{3/5}\,\theta_D, כך ששני המודלים מסכימים בחום. בצירים לינאריים שתי העקומות נבדלות רק סמוך לבסיס העלייה, וכל אחת מהן נראית כתיאור סביר. אחר כך הצירים עוברים לסקאלה לוגריתמית. עכשיו המדידות שמתחת לכ-10 K10\ \mathrm{K} נופלות על קו ישר בשיפוע 3 עד 1 K1\ \mathrm{K}, דבאי עוקב אחריהן, ואיינשטיין נופל מעריכית מתחת לכ-30 K30\ \mathrm{K} — נמוך מדי פי 550 בערך ב-10 K10\ \mathrm{K}, ויוצא מתחתית הגרף זמן קצר אחר כך. גרף לינארי מסתיר את ההבדל שמכריע בין המודלים.

שני תיקונים עומדים בין המודלים האלה לבין טבלת מדידות.

לאמת חישובית

כל מספר שלהלן בא מ-thermolab.radiation. שום דבר במודול הזה אינו אקראי; המספרים דטרמיניסטיים. אי-הוודאויות המצוטטות להתאמות של הכסף באות מהפיזור של הנתונים סביב עקומת המודל; זו של קרינת הרקע באה משגיאות המדידה שפורסמו. המעבדה חוזרת על כל ניסוי.

1. האכלוס, נבנה מחדש מהסכום של מודול 12. לאופן עם קוונט 4×10−21 J4 \times 10^{-21}\ \mathrm{J}, הצורה הסגורה ⟨n⟩=1/(ex−1)\langle n \rangle = 1/(e^{x} - 1) וסכום בולצמן מפורש על פני 4000 רמות, בעזרת הקוד של מודול 12, מסכימים עד 10-12 ב-30, 300 ו-3000 K3000\ \mathrm{K}. הנוסחה של שלב 1 היא פונקציית החלוקה של מודול 12, בלי שום תוספת.

2. ספירת גלים עומדים אחד אחד. במניית השלשות של מספרים שלמים בקובייה, אופן אחר אופן, והשוואה לספירה החלקה 8πVν3/(3c3)8\pi V\nu^3/(3c^3): הספירה המפורשת קטנה ממנה ב-11.2%11.2\% כאשר ∣n∣=20|\mathbf{n}| = 20, ב-2.8%2.8\% ב-80 וב-0.70%0.70\% ב-320. החוסר קטן בחצי בכל פעם שהרדיוס מוכפל — אופני הפאות החסרים, בערך (9/4)/∣n∣(9/4)/|\mathbf{n}| — והנוסחה החלקה היא הגבול שלו.

3. חוק סטפן–בולצמן באינטגרציה נומרית.

4. אנרגיה אחת, שתי צפיפויות, שני שיאים. אינטגרציה של uλu_\lambda על פני אורך הגל — משתנה אחר, על רשת לוגריתמית, שאין לה שום דבר משותף עם הרשת של פריט 3 — משחזרת את aT4aT^4 עד 10-16, ולכן מסכימה עם האינטגרל על פני התדירות של פריט 3 עד השגיאה של האינטגרל ההוא עצמו, 6×10−116 \times 10^{-11}, ב-2.725, 300 ו-5772 K5772\ \mathrm{K}. חיפוש ממצה על פני 200 001200\,001 נקודות מוצא את המקסימום של uνu_\nu ב-2.8214 kBT/h2.8214\,\kB T/h ואת זה של uλu_\lambda ב-hc/(4.9651 kBT)hc/(4.9651\,\kB T); המכפלה שלהם היא 0.5683 c0.5683\,c בכל טמפרטורה. האנרגיה מסכימה והשיאים אינם מסכימים, וזה בדיוק מה שצפיפות עושה.

5. חם יותר אינו בהיר יותר. הכפלת הטמפרטורה מ-300 ל-600 K600\ \mathrm{K} מכפילה את האנרגיה הכוללת פי 16.000. אבל בתדירות קבועה היחס בין שני הספקטרומים אינו מספר אחד: הוא 2 בקצה התדירות הנמוכה, שם שתי העקומות מצייתות לריילי–ג’ינס, 120 ב-60 THz60\ \mathrm{THz}, והוא גדל ללא גבול בזנב הרחוק, שם לגוף החם יותר יש פוטונים שהגוף הקר אינו יכול להרשות לעצמו. שום שינוי קנה מידה אנכי אינו הופך עקומה אחת לשנייה. השיא זז, מ-9.66 ל-4.83 μm4.83\ \mu\mathrm{m}.

6. שני המוצקים, זה מול זה. שני קיבולי החום מגיעים ל-3NkB3N\kB בטמפרטורה גבוהה, איינשטיין כמו 1−x2/121 - x^2/12 ודבאי כמו 1−x2/201 - x^2/20, וכל אחד מאושר עד הסדר הרביעי. בקור, השיפוע בצירים לוגריתמיים של עקומת דבאי מתחת ל-θD/50\theta_D/50 הוא 3 בתוך 10-15, והעקומה שווה שם לחוק ה-T3T^3 עד 10-11. בחלון שבין θ/40\theta/40 ל-θ/10\theta/10 השיפוע של עקומת איינשטיין הוא בערך 19. המודלים נבדלים זה מזה בשיפוע.

7. מול המציאות: כסף. הנתונים הם קיבולי החום המומלצים של כסף מ-1 עד 300 K300\ \mathrm{K} שבטבלאות המכון הלאומי לתקנים וטכנולוגיה של ארצות הברית (Smith and Fickett, 1995), מומרים לערכים למול.

שני לוחות. משמאל, גרף בצירים לוגריתמיים של קיבול החום של כסף כנגד טמפרטורה מ-1 K עד 300 K: נקודות נתונים שחורות עולות בתלילות ומתייצבות סמוך לקו אופקי מקווקו ב-3R. עקומת דבאי כחולה עוקבת אחרי הנקודות מבערך 4 K ומעלה; עקומת איינשטיין אדומה נופלת הרבה מתחתיהן בטמפרטורה נמוכה. מתחתיו, רצועה של שאריות מראה את הנקודות הכחולות בתוך כ-10 אחוזים ואת הנקודות האדומות צוללות למינוס 100 אחוז בטמפרטורה נמוכה. מימין, הספקטרום של קרינת הרקע הקוסמית כנגד תדירות: נקודות שחורות על עקומת פלאנק כחולה שמגיעה לשיא סביב 5.3 סנטימטר בחזקת מינוס אחת, עם רצועת שאריות מתחת שמפוזרת סביב אפס.

משמאל: קיבול החום המולרי של כסף CPC_P (נקודות שחורות), אחרי הפחתת ה-γT\gamma T של האלקטרונים, מול עקומת דבאי הטובה ביותר (כחול, θD=218 K\theta_D = 218\ \mathrm{K}) ועקומת איינשטיין הטובה ביותר (אדום, θE=123 K\theta_E = 123\ \mathrm{K}), שתיהן מותאמות בין 4 ל-100 K100\ \mathrm{K}; הקו המקווקו הוא 3R3R. הרצועה שמתחת מראה את השארית היחסית של כל מודל. מימין: קרינת הרקע הקוסמית כפי שנמדדה במכשיר FIRAS של COBE (נקודות שחורות; פסי השגיאה שלהן, הקטנים מהנקודות, אינם משורטטים) ועקומת פלאנק המותאמת (כחול), עם השאריות מתחת, מוגדלות, ביחידות של אי-הוודאות של כל נקודה.

8. מול המציאות: קרינת הרקע הקוסמית. מכשיר FIRAS של COBE מדד את הקרינה הממלאת את היקום ב-43 תדירויות בין 68 ל-640 GHz640\ \mathrm{GHz}. התאמה בת פרמטר אחד של חוק פלאנק, ברדיאנס, נותנת T=2.725 015±0.000 008 KT = 2.725\,015 \pm 0.000\,008\ \mathrm{K}, עם χ2=45.1\chi^2 = 45.1 עבור 42 דרגות חופש: עקומת פלאנק היא הצורה הנכונה, בתוך הדיוק של המדידה. כל הנקודות מלבד שתי הרועשות ביותר נמצאות בתוך 60 kJy sr−160\ \mathrm{kJy\,sr^{-1}} מהעקומה, 1.5 חלקים מתוך 104 מהשיא שלה. שתי הסתייגויות שייכות למספר הזה. אי-הוודאות המצוטטת היא סטטיסטית בלבד, וה-χ2\chi^2 מתייחס לשגיאות שפורסמו כבלתי תלויות, מה שהן אינן לגמרי; הכיול של FIRAS מגביל את הטמפרטורה המוחלטת ל-2.72548±0.00057 K2.72548 \pm 0.00057\ \mathrm{K}. והארכיון בונה את עמודת הספקטרום שלו כעקומת פלאנק של 2.725 K2.725\ \mathrm{K} ועוד השאריות המדודות, ולכן התאמה משחזרת את הטמפרטורה שהשאריות נלקחו ביחס אליה: מה שמערך הנתונים הזה בודק הוא הצורה, והצורה היא חוק פלאנק. FIRAS משרטט את הספקטרום שלו ליחידת תדירות, שם הוא מגיע לשיא ב-160 GHz160\ \mathrm{GHz}, אורך גל של 1.87 mm1.87\ \mathrm{mm}; ליחידת אורך גל אותה קרינה מגיעה לשיא ב-1.06 mm1.06\ \mathrm{mm}.

9. יהלום ועופרת, בטמפרטורת החדר. במדידה ב-298 K298\ \mathrm{K}, קיבול החום של יהלום הוא 6.12 J mol−1 K−16.12\ \mathrm{J\,mol^{-1}\,K^{-1}} (DeSorbo, 1953; Victor, 1962) — 25%25\% מדולונג–פטי — ושל עופרת 26.84 (NIST-JANAF), 8%8\% מעליו, כשהעודף הוא CP−CVC_P - C_V ואי-הרמוניות. האטומים הקלים והקשרים הנוקשים של היהלום נותנים לו טמפרטורת דבאי סביב 2230 K2230\ \mathrm{K} מנתוני טמפרטורה נמוכה; זו של עופרת היא בערך 105 K105\ \mathrm{K}. עם הערכים האלה מודל דבאי מציב את היהלום על 16%16\% מ-3R3R בטמפרטורת החדר, נמוך מדי: ה-θD\theta_D שמשחזר את הערך בטמפרטורת החדר הוא 1860 K1860\ \mathrm{K}. גם טמפרטורת דבאי של יהלום נודדת עם TT.

ניסויי ההפרכה. לכל תפיסה שגויה שהמודול הזה מתעמת איתה יש מדידה לעיל שהיא אינה שורדת. “שיא הספקטרום הוא תכונה של הקרינה”: פריט 4, שבו שני השיאים נבדלים בגורם 0.568 בכל טמפרטורה בעוד האנרגיות מסכימות עד 6×10−116 \times 10^{-11}, ופריט 8, שבו קרינת הרקע הקוסמית מגיעה לשיא ב-1.87 או ב-1.06 mm1.06\ \mathrm{mm} בהתאם לציר. “גופים חמים יותר הם גרסאות בהירות יותר של אותו ספקטרום”: פריט 5, שבו הכפלת TT מכפילה את האנרגיה פי שש-עשרה בעוד היחס בין שני הספקטרומים רץ מ-2 ומעלה ללא גבול. “האטומים של מוצק קר מפסיקים לזוז, ולכן הוא אינו אוגר חום”: פריטים 6 ו-7 — קיבול החום מתאפס כשהקוונטים נעשים יקרים מכדי לשלם עליהם, בקצב שנקבע לפי מספר האופנים הזולים שיש לגביש, בעוד תנועת נקודת האפס ששלב 1 הניח בצד נמשכת באפס המוחלט.

להעביר את הרעיון

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

במבט קדימה:

דוגמאות פתורות

נסו כל אחת לפני שתפתחו את הפתרון שלה.

דוגמה פתורה 1 — שני השיאים של השמש. פני השמש קורנים כגוף שחור ב-5772 K5772\ \mathrm{K}. מצאו את אורך הגל שבו uλu_\lambda מגיע לשיא, ואת התדירות שבה uνu_\nu מגיע לשיא. לאיזה אורך גל של אור מתאימה התדירות הזו?

פתרון

λmax⁡=2.898×10−3/5772=502 nm\lambda_{\max} = 2.898 \times 10^{-3}/5772 = 502\ \mathrm{nm}. νmax⁡=2.821 kBT/h=2.821×1.381×10−23×5772/6.626×10−34=3.39×1014 Hz\nu_{\max} = 2.821\,\kB T/h = 2.821 \times 1.381 \times 10^{-23} \times 5772/6.626 \times 10^{-34} = 3.39 \times 10^{14}\ \mathrm{Hz}, התדירות של אור באורך גל c/νmax⁡=883 nmc/\nu_{\max} = 883\ \mathrm{nm}, באינפרה-אדום הקרוב. היחס 502/883=0.568502/883 = 0.568, כפי שהוא חייב להיות בכל טמפרטורה.

דוגמה פתורה 2 — תנור, מוכפל. הדפנות של תנור מועלות מ-300 ל-600 K600\ \mathrm{K}. באיזה גורם גדל ההספק שהן מקרינות, והיכן מגיע לשיא הספקטרום ליחידת אורך גל לפני ואחרי?

פתרון

ההספק משתנה כמו T4T^4: 24=162^4 = 16. השיא זז מ-2.898×10−3/300=9.66 μm2.898 \times 10^{-3}/300 = 9.66\ \mu\mathrm{m} ל-4.83 μm4.83\ \mu\mathrm{m}. ספקטרום שרק נעשה בהיר יותר היה שומר על השיא שלו וגדל בגורם 2, יחס הטמפרטורות — וזה מה שהוא עושה בקצה האינפרה-אדום הרחוק, ובשום מקום אחר.

דוגמה פתורה 3 — הפוטונים בחדר. כמה פוטונים יש בסנטימטר מעוקב של חדר ב-300 K300\ \mathrm{K}, ובסנטימטר מעוקב של החלל הבין-גלקטי, המלא בקרינת הרקע של 2.725 K2.725\ \mathrm{K}? מהי האנרגיה הממוצעת של פוטון רקע?

פתרון

n=2.03×107 T3 m−3n = 2.03 \times 10^{7}\,T^3\ \mathrm{m^{-3}} נותן 5.5×1014 m−35.5 \times 10^{14}\ \mathrm{m^{-3}} ב-300 K300\ \mathrm{K}, 5.5×1085.5 \times 10^{8} לסנטימטר מעוקב, ו-4.1×108 m−34.1 \times 10^{8}\ \mathrm{m^{-3}} — 411 לסנטימטר מעוקב — ב-2.725 K2.725\ \mathrm{K}. האנרגיה הממוצעת לפוטון היא u/n=(π4/30ζ(3)) kBT=2.70 kBTu/n = (\pi^4/30\zeta(3))\,\kB T = 2.70\,\kB T, 6.3×10−4 eV6.3 \times 10^{-4}\ \mathrm{eV} עבור קרינת הרקע.

דוגמה פתורה 4 — קיבול החום הקבוע של מודול 01. מודול 01 התייחס למוצק איינשטיין כאילו יש לו C=kBC = \kB למתנד. איזה חלק מזה נותן מודל איינשטיין ב-T=θE/2T = \theta_E/2, ב-T=θET = \theta_E וב-T=5θET = 5\theta_E?

פתרון

C/(3NkB)=x2ex/(ex−1)2C/(3N\kB) = x^2 e^x/(e^x - 1)^2. ב-x=2x = 2: 4×7.389/6.3892=0.7244 \times 7.389/6.389^2 = 0.724. ב-x=1x = 1: 0.921. ב-x=0.2x = 0.2: 0.9967, לעומת הטור 1−0.04/12=0.99671 - 0.04/12 = 0.9967. הקירוב של מודול 01 טוב עד 8%8\% ב-T=θET = \theta_E ומתדרדר במהירות מתחת לזה.

דוגמה פתורה 5 — כסף ב-4 K. מתוך θD=226 K\theta_D = 226\ \mathrm{K} ו-γ=0.641 mJ mol−1 K−2\gamma = 0.641\ \mathrm{mJ\,mol^{-1}\,K^{-2}}, נבאו את קיבול החום המולרי של כסף ב-4 K4\ \mathrm{K}, והשוו לערך שבטבלה, 0.0134 J mol−1 K−10.0134\ \mathrm{J\,mol^{-1}\,K^{-1}}.

פתרון

סריג: 1944×(4/226)3=1944×5.54×10−6=0.0108 J mol−1 K−11944 \times (4/226)^3 = 1944 \times 5.54 \times 10^{-6} = 0.0108\ \mathrm{J\,mol^{-1}\,K^{-1}}. אלקטרונים: 0.641×10−3×4=0.00260.641 \times 10^{-3} \times 4 = 0.0026. סך הכול 0.0133, לעומת 0.0134 שנמדד. בלי האלקטרונים הניבוי היה נמוך ב-20%20\%; בטמפרטורה הזו הם חמישית מהתשובה.

דוגמה פתורה 6 — ההספק של השמש. הרדיוס של השמש הוא 6.957×108 m6.957 \times 10^8\ \mathrm{m} וטמפרטורת פני השטח שלה 5772 K5772\ \mathrm{K}. העריכו את הלומינוסיות שלה.

פתרון

L=4πR2σT4=4π(6.957×108)2×5.670×10−8×57724=3.83×1026 WL = 4\pi R^2 \sigma T^4 = 4\pi (6.957 \times 10^8)^2 \times 5.670 \times 10^{-8} \times 5772^4 = 3.83 \times 10^{26}\ \mathrm{W}. הערך שאומץ על ידי האיגוד האסטרונומי הבין-לאומי הוא 3.828×1026 W3.828 \times 10^{26}\ \mathrm{W} — לא צירוף מקרים: “הטמפרטורה האפקטיבית” של השמש, 5772 K5772\ \mathrm{K}, מוגדרת כטמפרטורת הגוף השחור שמקרין את ההספק המדוד שלה.

בחנו את הבנתכם

  1. לאופן של השדה האלקטרומגנטי יש קוונט h nu = 3 k_B T. מהי האנרגיה הממוצעת שלו מעל מצב היסוד, ביחידות של k_B T?

    A. 3/(e^3 - 1) = 0.157 -- האופן נמצא רוב הזמן במצב היסוד שלו, כי קוונט שלם עולה פי שלושה מהאנרגיה התרמית האופיינית. B. 1 -- כל אופן נושא k_B T ללא תלות בתדירות שלו. C. 3 -- האופן מחזיק בדיוק קוונט אחד. D. 0 -- אי אפשר לעולם לעורר קוונט הגדול מ-k_B T.

    בדקו את התשובה

    התשובה הנכונה: A — 3/(e^3 - 1) = 0.157 -- האופן נמצא רוב הזמן במצב היסוד שלו, כי קוונט שלם עולה פי שלושה מהאנרגיה התרמית האופיינית.

    • A. נכון. = h nu/(e^(h nu/(k_B T)) - 1) = 3 k_B T/(e^3 - 1) = 0.157 k_B T. האופן אינו קפוא לגמרי -- בערך פעם אחת מעשרים הוא מחזיק קוונט -- אבל הוא מחזיק הרבה פחות מה-k_B T הקלאסי, ואופנים בתדירות גבוהה יותר מחזיקים פחות באופן אקספוננציאלי.

    • B. זו חלוקת האנרגיה השווה, ההנחה הקלאסית שהמודול הזה מפיל. היא הגבול של התוצאה הקוונטית כאשר h nu << k_B T; ב-h nu = 3 k_B T היא מעריכה ביתר את האנרגיה פי שישה, וכשמחילים אותה על כל אופן היא נותנת את הקטסטרופה האולטרה-סגולה.

    • C. מספר הקוונטים אינו קבוע; הוא עובר פלוקטואציות, והממוצע שלו הוא = 1/(e^3 - 1) = 0.052. בטמפרטורה הזו אופן מחזיק קוונט אחד רק בכ-5% מהזמן, ואפס קוונטים ברוב הזמן.

    • D. קוונטים גדולים מדוכאים באופן אקספוננציאלי, לא אסורים. גורם בולצמן לקוונט אחד הוא e^-3 = 0.05, קטן אבל לא אפס; הקפיאה היא דעיכה אקספוננציאלית חלקה, לא קטיעה חדה.

  2. הפיזיקה הקלאסית מנבאת שהקרינה בחלל חם מחזיקה אנרגיה אינסופית. מה, במדויק, מתבדר?

    A. מספר האופנים ליחידת תדירות גדל כמו nu^2 ללא גבול, וחלוקת האנרגיה השווה נותנת לכל אחד מהם k_B T, ולכן האינטגרל של 8 pi nu^2 k_B T / c^3 על פני כל התדירויות אינסופי. B. פוטונים נעים במהירות האור, ולכן האנרגיה שלהם אינסופית. C. הדפנות בולעות ופולטות מחדש ללא גבול, ולכן האנרגיה מצטברת לעד. D. האנרגיה של כל אופן גדלה עם התדירות שלו, ולכן האופנים הגבוהים ביותר נושאים אנרגיה אינסופית.

    בדקו את התשובה

    התשובה הנכונה: A — מספר האופנים ליחידת תדירות גדל כמו nu^2 ללא גבול, וחלוקת האנרגיה השווה נותנת לכל אחד מהם k_B T, ולכן האינטגרל של 8 pi nu^2 k_B T / c^3 על פני כל התדירויות אינסופי.

    • A. נכון. ספירת האופנים בחלל נותנת g(nu) = 8 pi V nu^2 / c^3: אין גל עומד קצר ביותר, ולכן אין אופן אחרון. כל אופן הוא מתנד שמחזיק k_B T באופן קלאסי. התיקון הקוונטי שומר על מספר האופנים ומשנה את האנרגיה לאופן, שיורדת אקספוננציאלית ברגע ש-h nu עולה על כמה k_B T.

    • B. האנרגיה של פוטון היא h nu, סופית בכל תדירות. ההתבדרות היא בעיית ספירה: אינסוף אופנים, שכל אחד מהם מקבל אותה אנרגיה קלאסית סופית.

    • C. בשיווי משקל הדפנות פולטות בדיוק כמה שהן בולעות; שום דבר אינו מצטבר. הקטסטרופה הקלאסית היא טענה על מצב שיווי המשקל עצמו, לא על תהליך שיוצא משליטה.

    • D. באופן קלאסי כל אופן נושא אותו k_B T, ללא תלות בתדירות שלו. מה שגדל ללא גבול הוא מספר האופנים, לא האנרגיה של כל אחד מהם.

  3. כשהוא משורטט ליחידת אורך גל, הספקטרום של השמש (גוף שחור בסביבות 5772 K) מגיע לשיא בכ-500 nm. כשהוא משורטט ליחידת תדירות, היכן נמצא השיא שלו?

    A. בכ-340 THz, התדירות של אור באורך גל 880 nm -- פוטון אחר, כי u_nu ו-u_lambda הן צפיפויות הקשורות זו לזו דרך היעקוביאן c / lambda^2, שמזיז את המקסימום. B. ב-c / 500 nm = 600 THz -- אותו ספקטרום, אותו שיא. C. ב-1200 THz, כי התדירות נמצאת ביחס הפוך לאורך הגל, ולכן שיא באורך גל קצר הוא שיא בתדירות גבוהה. D. אין שיא ליחידת תדירות; u_nu גדלה ללא גבול.

    בדקו את התשובה

    התשובה הנכונה: A — בכ-340 THz, התדירות של אור באורך גל 880 nm -- פוטון אחר, כי u_nu ו-u_lambda הן צפיפויות הקשורות זו לזו דרך היעקוביאן c / lambda^2, שמזיז את המקסימום.

    • A. נכון. nu_max = 2.821 k_B T/h = 339 THz, ו-c/nu_max = 883 nm. שני השיאים נבדלים זה מזה בגורם nu_max lambda_max / c = 0.568 בכל טמפרטורה. האנרגיה בכל פס היא אותה אנרגיה בשני התיאורים; מיקומו של “השיא” הוא תכונה של הגרף, לא של האור.

    • B. זו התפיסה השגויה שהשיא הוא תכונה של הקרינה לבדה. u_nu ו-u_lambda אינן אותה פונקציה עם תווית אחרת: u_lambda = u_nu c / lambda^2, והגורם c / lambda^2 משתנה מנקודה לנקודה, ולכן המקסימומים נופלים על פוטונים שונים. הכלי החוקר מציג את שני הסמנים; הסמן הממופה לעולם אינו נוחת על השני.

    • C. ההיפוך נכון להמרת אורך הגל של פוטון אחד לתדירות שלו (500 nm הם 600 THz, לא 1200), אבל השאלה היא היכן הצפיפות ליחידת תדירות היא הגדולה ביותר, וזה מחייב את היעקוביאן, לא רק את ההיפוך.

    • D. זו העקומה הקלאסית של ריילי–ג’ינס. ל-u_nu של פלאנק יש מקסימום ב-x = h nu / (k_B T) = 2.821, ומעבר לו היא יורדת אקספוננציאלית.

  4. טמפרטורת פני השטח של סיריוס A היא כ-9940 K. אם מתייחסים אליו כאל גוף שחור, באיזה אורך גל, בננומטרים, מגיע הספקטרום שלו ליחידת אורך גל לשיא? (השתמשו ב-lambda_max T = 2.898e-3 m K.)

    בדקו את התשובה

    התשובה הנכונה: 292 (± 2)

    lambda_max = 2.898e-3 / 9940 = 2.92e-7 m = 292 nm, באולטרה-סגול. ליחידת תדירות השיא הוא nu_max = 2.821 k_B T/h = 5.84e14 Hz, התדירות של אור באורך גל 513 nm; היחס nu_max lambda_max / c = 0.568, כמו בכל טמפרטורה. אם ענית 513, חישבת את שיא התדירות והמרת אותו לאורך גל.

  5. הדפנות של תנור מחוממות מ-300 K ל-600 K. באיזה גורם גדל ההספק הכולל שהן מקרינות?

    בדקו את התשובה

    התשובה הנכונה: 16 (± 0.1)

    ההספק גדל כמו sigma T^4, ולכן הגורם הוא 2^4 = 16. אם ענית 2, הנחת שגוף חם יותר מקרין את אותו ספקטרום, רק בהיר יותר ביחס ישר לטמפרטורה שלו: זה נכון רק בקצה הרחוק של התדירויות הנמוכות. הספקטרום גם מרוויח פוטונים בתדירויות שהגוף הקר יותר לא יכול היה להרשות לעצמו, והשיא שלו נע מ-9.66 ל-4.83 מיקרומטרים. אם ענית 8, ספרת פוטונים (T^3) במקום אנרגיה.

  6. מדוע הפוטנציאל הכימי של גז הפוטונים בחלל הוא אפס?

    A. הדפנות בולעות ופולטות פוטונים בחופשיות, ולכן מספר הפוטונים אינו נשמר; הוא מתכוונן כך שהאנרגיה החופשית תהיה מינימלית ב-T וב-V קבועים, שם (dF/dN)_{T,V} = mu = 0. B. מספר הפוטונים נשמר, ולכן mu נקבע על ידו ובמקרה יוצא אפס. C. פוטונים הם חסרי מסה, ופוטנציאל כימי הוא אנרגיה של מסה. D. פוטונים אינם מקיימים אינטראקציה זה עם זה, ו-mu מודד אינטראקציות.

    בדקו את התשובה

    התשובה הנכונה: A — הדפנות בולעות ופולטות פוטונים בחופשיות, ולכן מספר הפוטונים אינו נשמר; הוא מתכוונן כך שהאנרגיה החופשית תהיה מינימלית ב-T וב-V קבועים, שם (dF/dN)_{T,V} = mu = 0.

    • A. נכון. מספר שנשמר הוא אילוץ, והפוטנציאל הכימי שלו הוא המחיר של הוספת עוד אחד; מספר שאינו נשמר הוא משתנה שמתייצב במינימום של F. תוצאה שאפשר לבדוק: מספר הפוטונים אינו קבוע אלא גדל כמו T^3 כשמחממים את החלל, ו-U = TS - PV מתקיים בדיוק.

    • B. מספר הפוטונים אינו נשמר בחלל שיש לו דפנות: חימום החלל יוצר פוטונים (n גדל כמו T^3), וקירור משמיד אותם. אילו היה נשמר, mu היה מקבל את הערך שקובע את המספר -- בדרך כלל לא אפס -- כמו לאטומים של מודול 17.

    • C. הפוטנציאל הכימי הוא המחיר, באנרגיה חופשית, של הוספת חלקיק ב-T וב-V קבועים, לא אנרגיית מנוחה. היעדר המסה לבדו אינו קובע mu = 0; אי-השימור קובע. אור בתיבה עם מראות מושלמות ובלי שום דבר שיבלע אותו היה משמר את מספר הפוטונים.

    • D. לגז של חלקיקים נשמרים ללא אינטראקציה, כמו הגז האידיאלי של מודול 13, יש mu = k_B T ln(n / n_Q) תקין לחלוטין ושונה מאפס. הסיבה היא אי-השימור, לא היעדר האינטראקציות.

  7. במודל איינשטיין, מהו קיבול החום ב-T = theta_E / 2, ביחידות של ערך דולונג–פטי 3 N k_B?

    בדקו את התשובה

    התשובה הנכונה: 0.724 (± 0.005)

    עם x = theta_E / T = 2, C/(3 N k_B) = x^2 e^x/(e^x - 1)^2 = 4 e^2/(e^2 - 1)^2 = 0.724. ב-T = 5 theta_E (x = 0.2) הוא 0.9967, בהתאם לטור 1 - x^2/12: דולונג–פטי הוא הגבול של טמפרטורה גבוהה. אם ענית 0.541, השתמשת ב-x^2 e^-x, הקירוב של טמפרטורה נמוכה, מחוץ לתחום התקפות שלו.

  8. קיבול החום של גביש מבדד, שנמדד מתחת ל-10 K, יורד בשיפוע לוג-לוג 3. איזה מודל שורד, ומדוע C יורד בכלל?

    A. דבאי -- לגביש יש אופני קול בתדירות נמוכה כרצוננו, שהקוונטים שלהם נשארים זולים מספיק; ככל ש-T יורד, פחות אופנים (בערך V (k_B T / hbar c_s)^3 מהם) יכולים לקבל אנרגיה, ולכן C יורד כמו T^3. B. איינשטיין -- שני המודלים מתאפסים ב-T = 0, ולכן כל אחד מהם מתאים. C. אף אחד מהם -- האטומים קפואים ללא תנועה בטמפרטורה נמוכה, ולכן אי אפשר לאגור אנרגיה ו-C צונח לאפס בבת אחת. D. איינשטיין, בתנאי ש-theta_E נבחר קטן מספיק.

    בדקו את התשובה

    התשובה הנכונה: A — דבאי -- לגביש יש אופני קול בתדירות נמוכה כרצוננו, שהקוונטים שלהם נשארים זולים מספיק; ככל ש-T יורד, פחות אופנים (בערך V (k_B T / hbar c_s)^3 מהם) יכולים לקבל אנרגיה, ולכן C יורד כמו T^3.

    • A. נכון. הספקטרום של דבאי מכיל תמיד אופנים עם hbar omega < k_B T, וספירתם נותנת T^3, כמו שספירת פוטונים נותנת את ה-T^3 של מספר הפוטונים. C יורד מפני שהקוונטים נעשים יקרים מדי, אופן אחר אופן; האטומים שומרים על תנועת נקודת האפס שלהם.

    • B. שניהם מתאפסים, אבל ה-C של איינשטיין יורד כמו x^2 e^-x, אקספוננציאלית, הרבה יותר מהר מכל חזקה: שיפוע הלוג-לוג שלו אינו קבוע, והוא כ-19 בין theta/40 ל-theta/10. השיפוע מבחין בין המודלים, והנתונים בוחרים בדבאי.

    • C. זו התפיסה השגויה שאטומים קרים מפסיקים לנוע. כל מתנד שומר על אנרגיית נקודת האפס שלו hbar omega/2 ב-T = 0, וקיבול החום יורד בהדרגה, כמו T^3, לא בבת אחת. C קטן מפני שהאמבט רק לעתים רחוקות מספק קוונט שלם, לא מפני שהתנועה נפסקה.

    • D. theta_E קטן יותר דוחה את הקריסה אבל אינו משנה את צורתה: מתחת לכ-theta_E / 5 עקומת איינשטיין עדיין יורדת אקספוננציאלית. שום תדירות יחידה אינה נותנת חוק חזקה.

  9. הטמפרטורה של חלל סגור מוכפלת. באיזה גורם משתנה מספר הפוטונים שבתוכו?

    בדקו את התשובה

    התשובה הנכונה: 8 (± 0.05)

    n = 16 pi zeta(3) (k_B T / h c)^3, ולכן הכפלת T כופלת את מספר הפוטונים ב-2^3 = 8. הפוטונים נוצרים על ידי הדפנות ואינם נשמרים. האנרגיה גדלה פי 16 (T^4), ולכן האנרגיה הממוצעת לפוטון, 2.70 k_B T, מוכפלת. אם ענית 2, התייחסת לפוטונים כמו למולקולות של גז בכלי אטום, שאת מספרן חימום אינו יכול לשנות.

  10. בארבעה או חמישה משפטים: ציינו את ההנחה הקלאסית היחידה שעומדת מאחורי הקטסטרופה האולטרה-סגולה וגם מאחורי הכישלון של דולונג–פטי בטמפרטורה נמוכה, הסבירו כיצד שינוי יחיד מציל את שניהם, והסבירו מדוע קיבול החום של גביש יורד לעבר אפס אף שהאטומים שלו לעולם אינם מפסיקים לנוע.

בדקו את התשובה

שני הניבויים נשענים על חלוקת האנרגיה השווה: כל אופן -- כל גל עומד בחלל, כל תנודה של גביש -- מחזיק בממוצע k_B T, תהא התדירות שלו אשר תהא. עם אינסוף אופני חלל זה נותן אנרגיה אינסופית; עם 3N אופני גביש זה נותן C = 3 N k_B בכל הטמפרטורות. השינוי היחיד הוא קוונטיזציה: אופן בתדירות nu מחזיק אנרגיה רק בקוונטים של h nu, ולכן האנרגיה הממוצעת שלו היא h nu/(e^(h nu/(k_B T)) - 1), שהיא k_B T עבור קוונטים זולים וקטנה אקספוננציאלית עבור יקרים. אופני חלל בתדירות גבוהה, וכשקר גם אופני הגביש, אינם יכולים להרשות לעצמם קוונט ונושרים מהחשבון. האטומים שומרים על תנועת נקודת האפס שלהם ב-T = 0; מה שנעלם הוא היכולת שלהם לקבל כמויות קטנות של אנרגיה, וזה בדיוק מה שקיבול חום מודד. תשובה חזקה מוסיפה שהאופנים בתדירות נמוכה של דבאי משאירים כמה אופנים בהישג יד, ונותנים T^3 במקום אקספוננט.

בעיות בסגנון בחינה למודול הזה — אלה ששווה לפתור עם עט — מרוכזות באוסף הבעיות.

להסביר במילים שלכם

ענו בכמה משפטים על כל אחת.

  1. מדוע “השיא של ספקטרום הגוף השחור” אינו גודל מוגדר היטב עד שנוקבים במשתנה ספקטרלי? מה כן מוגדר היטב?

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

  3. ציינו את ההנחה הקלאסית האחת שמאחורי הקטסטרופה האולטרה-סגולה ומאחורי הכישלון של דולונג–פטי בטמפרטורה נמוכה, והסבירו כיצד שינוי יחיד בה מציל את שניהם.

  4. במשפט אחד: מדוע הפוטנציאל הכימי של פוטונים הוא אפס?

מתקדם: מאזן האנרגיה של כוכב לכת

כוכב הלכת החשוף. כדור הארץ מיירט אור שמש בחתך שלו, πR2\pi R^2, וקורן מכל פני השטח שלו, 4πR24\pi R^2. הקבוע הסולרי, שטף אור השמש במרחק של כדור הארץ, הוא S=1361 W m−2S = 1361\ \mathrm{W\,m^{-2}}, וכוכב הלכת מחזיר חלק A≈0.3A \approx 0.3 ממנו, האלבדו שלו. אם פני השטח קורנים כגוף שחור בטמפרטורה אחת TeT_e, באיזון בין מה שנבלע למה שנפלט,

σTe4=(1−A) S4=238 W m−2,Te=255 K.\sigma T_e^4 = \frac{(1 - A)\,S}{4} = 238\ \mathrm{W\,m^{-2}} , \qquad T_e = 255\ \mathrm{K} .

כלומר −19 ∘C-19\,{}^\circ\mathrm{C}. הטמפרטורה הממוצעת המדודה של פני השטח היא בערך 288 K288\ \mathrm{K}, 15 ∘C15\,{}^\circ\mathrm{C}.

שכבה בולעת אחת. הוסיפו אטמוספרה, באידיאליזציה של שכבה יחידה שמעבירה אור שמש ובולעת את כל הקרינה התרמית. חוק ההסחה של וין הוא מה שהופך את האידיאליזציה לסבירה: אור השמש מגיע לשיא ליחידת אורך גל סביב 0.5 μm0.5\ \mu\mathrm{m}, כוכב לכת של 255 K255\ \mathrm{K} סביב 11 μm11\ \mu\mathrm{m}, וגזים כמו אדי מים ופחמן דו-חמצני בולעים חזק בתחום השני וכמעט שאינם בולעים בראשון. השכבה, בטמפרטורה TaT_a, קורנת גם כלפי מעלה וגם כלפי מטה. בשיווי משקל החלק העליון של השכבה חייב לפלוט לחלל את מה שכוכב הלכת בולע, והשכבה חייבת לפלוט כמה שהיא מקבלת מפני השטח:

σTa4=(1−A) S4,σTs4=2 σTa4,\sigma T_a^4 = \frac{(1 - A)\,S}{4} , \qquad \sigma T_s^4 = 2\,\sigma T_a^4 ,

ולכן Ta=TeT_a = T_e ופני השטח נמצאים ב-Ts=21/4 Te=303 KT_s = 2^{1/4}\,T_e = 303\ \mathrm{K}: חמים יותר, מפני שהם מקבלים את הקרינה של השכבה כלפי מטה בנוסף לזו של השמש. עם nn שכבות כאלה, Ts=(n+1)1/4 TeT_s = (n + 1)^{1/4}\,T_e.

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