← กลับไปยังบทความทั้งหมด

สัมประสิทธิ์วิเรียล (Virial Coefficient) การแก้ไขกฎแก๊สอุดมคติ

สรุปใจความสำคัญ

  • สัมประสิทธิ์วิเรียลใช้เพื่อแก้ไขกฎแก๊สอุดมคติเพื่อให้สามารถอธิบายพฤติกรรมของแก๊สจริงได้แม่นยำขึ้น
  • สัมประสิทธิ์วิเรียลตัวที่สอง (B2) ขึ้นอยู่กับปฏิสัมพันธ์ระหว่างอนุภาคเพียงสองตัว (Pair Interaction)
  • ในขีดจำกัดคลาสสิก สัมประสิทธิ์วิเรียลจะขึ้นอยู่กับศักย์ปฏิสัมพันธ์ระหว่างอนุภาคเท่านั้น โดยไม่ขึ้นกับพลังงานจลน์

ในทางกลศาสตร์สถิติ สัมประสิทธิ์วิเรียล (Virial Coefficients)

B_i
สัมประสิทธิ์วิเรียล B_i
ปรากฏเป็นสัมประสิทธิ์ในการขยายตัวของวิเรียล (Virial Expansion) ของความดันในระบบที่มีอนุภาคจำนวนมากในรูปของความหนาแน่น ซึ่งเป็นการให้ค่าการแก้ไขอย่างเป็นระบบต่อกฎแก๊สอุดมคติ

สัมประสิทธิ์เหล่านี้เป็นลักษณะเฉพาะของศักย์ปฏิสัมพันธ์ (Interaction Potential) ระหว่างอนุภาค และโดยทั่วไปจะขึ้นอยู่กับอุณหภูมิ โดยที่สัมประสิทธิ์วิเรียลตัวที่สอง

B_2
สัมประสิทธิ์วิเรียลตัวที่สอง B_2
จะขึ้นอยู่กับปฏิสัมพันธ์แบบคู่ (Pair Interaction) ระหว่างอนุภาคเท่านั้น ในขณะที่สัมประสิทธิ์ตัวที่สาม
B_3
สัมประสิทธิ์วิเรียลตัวที่สาม B_3
จะขึ้นอยู่กับปฏิสัมพันธ์แบบ 2 และ 3 อนุภาค (Non-additive 3-body interactions) และเป็นเช่นนี้ต่อไปสำหรับสัมประสิทธิ์ลำดับที่สูงขึ้น

การอนุพัทธ์ (Derivation)

ขั้นตอนแรกในการหาค่าสัมประสิทธิ์วิเรียลคือการขยายแบบคลัสเตอร์ (Cluster Expansion) ของฟังก์ชันแบ่งส่วนแกรนด์แคนอนิคัล (Grand Canonical Partition Function)

Grand Canonical Partition Function
สมการฟังก์ชันแบ่งส่วนแกรนด์แคนอนิคัล

โดยที่:

  • p
    p
    ความดัน (Pressure)
    คือความดัน
  • V
    V
    ปริมาตร (Volume)
    คือปริมาตรของภาชนะที่บรรจุอนุภาค
  • kB
    k_B
    ค่าคงที่โบลต์ซมันน์
    คือค่าคงที่โบลต์ซมันน์
  • T
    T
    อุณหภูมิ (Temperature)
    คืออุณหภูมิสัมบูรณ์
  • λ
    fugacity
    ฟูกาซิตี้ (Fugacity)
    คือฟูกาซิตี้ โดยที่ μ
    mu
    ศักย์เคมี (Chemical Potential)
    คือศักย์เคมี
  • Qn
    Q_n
    ฟังก์ชันแบ่งส่วนแคนอนิคัล
    คือฟังก์ชันแบ่งส่วนแคนอนิคัลของระบบย่อยที่มีอนุภาค n ตัว
    n
    จำนวนอนุภาค n
Q_n expression
นิยามของ Q_n ในรูปของ Trace

โดยที่ H(1, 2, ..., n)

H
แฮมิลโทเนียน
คือแฮมิลโทเนียน (ตัวดำเนินการพลังงาน) ของระบบย่อยที่มีอนุภาค n ตัว ซึ่งเป็นผลรวมของพลังงานจลน์ของอนุภาคและพลังงานศักย์รวม (พลังงานปฏิสัมพันธ์) ของอนุภาค n ตัว ซึ่งรวมถึงปฏิสัมพันธ์แบบคู่และอาจรวมถึงปฏิสัมพันธ์แบบ 3 ตัวหรือมากกว่านั้น

ฟังก์ชันแบ่งส่วนแกรนด์แคนอนิคัล Ξ

Xi
สัญลักษณ์ Ξ
สามารถขยายได้ในรูปของผลรวมของผลกระทบจากคลัสเตอร์หนึ่งตัว, สองตัว, ฯลฯ การขยายตัวของวิเรียลจะได้มาจากการสังเกตว่า ln Ξ
ln Xi
ลอการิทึมของ Ξ
มีค่าเท่ากับ pV/(kBT)
pV/kBT
ความสัมพันธ์ระหว่างความดันและฟังก์ชันแบ่งส่วน

จากกระบวนการนี้ เราสามารถอนุพัทธ์ค่าสัมประสิทธิ์วิเรียลได้ดังนี้:

B_2 formula
สูตรคำนวณสัมประสิทธิ์วิเรียลตัวที่สอง B_2
B_3 formula
สูตรคำนวณสัมประสิทธิ์วิเรียลตัวที่สาม B_3

นี่คือนิยามทางสถิติควอนตัมที่รวมพลังงานจลน์ไว้ด้วย โปรดสังเกตว่าฟังก์ชันแบ่งส่วนของอนุภาคตัวเดียว Q1

Q_1
ฟังก์ชันแบ่งส่วน Q_1
ประกอบด้วยพจน์พลังงานจลน์เท่านั้น ในขีดจำกัดคลาสสิก ℏ = 0
hbar=0
ขีดจำกัดคลาสสิก
ตัวดำเนินการพลังงานจลน์จะสลับที่กับตัวดำเนินการพลังงานศักย์ และพลังงานจลน์ในตัวเศษและตัวส่วนจะหักล้างกันเอง ทำให้ Trace (tr) กลายเป็นอินทิกรัลเหนือปริภูมิโครงสร้าง (Configuration Space) ส่งผลให้สัมประสิทธิ์วิเรียลแบบคลาสสิกขึ้นอยู่กับปฏิสัมพันธ์ระหว่างอนุภาคเท่านั้น และคำนวณได้จากอินทิกรัลเหนือพิกัดของอนุภาค

การอนุพัทธ์สัมประสิทธิ์วิเรียลที่สูงกว่า B3

B_3 repeat
สัมประสิทธิ์วิเรียลตัวที่สาม
จะกลายเป็นปัญหาทางคอมบิเนทอริกที่ซับซ้อนอย่างรวดเร็ว หากใช้การประมาณแบบคลาสสิกและละทิ้งปฏิสัมพันธ์แบบไม่เพิ่มพูน (Non-additive interactions) หากมีอยู่ คอมบิเนทอริกนี้สามารถจัดการได้ด้วยแผนภาพ (Graphical method) ตามที่ Joseph E. Mayer และ Maria Goeppert-Mayer ได้นำเสนอเป็นครั้งแรก

คำถามที่พบบ่อย

สัมประสิทธิ์วิเรียลคืออะไร?

คือค่าสัมประสิทธิ์ที่ใช้ในการขยายตัวของความดันของระบบอนุภาคในรูปของความหนาแน่น เพื่อแก้ไขกฎแก๊สอุดมคติให้เข้ากับพฤติกรรมของแก๊สจริง

B2 และ B3 แตกต่างกันอย่างไร?

B2 ขึ้นอยู่กับปฏิสัมพันธ์ระหว่างอนุภาคเพียง 2 ตัว ในขณะที่ B3 ขึ้นอยู่กับปฏิสัมพันธ์ระหว่างอนุภาค 2 และ 3 ตัวพร้อมกัน

ทำไมต้องใช้การขยายตัวของวิเรียลในทางฟิสิกส์?

เพราะในความเป็นจริง อนุภาคของแก๊สมีขนาดและมีแรงดึงดูดหรือแรงผลักระหว่างกัน ซึ่งกฎแก๊สอุดมคติไม่สามารถอธิบายได้ สัมประสิทธิ์วิเรียลจึงถูกนำมาใช้เพื่อให้คำนวณความดันได้ถูกต้องตามความเป็นจริงมากขึ้น