ในกลศาสตร์เชิงสถิติโมเดลแบบน้ำแข็งหรือโมเดลหกจุดยอดเป็นกลุ่มของโมเดลจุดยอดสำหรับโครงผลึกที่มีพันธะไฮโดรเจน โมเดลแรกดังกล่าวได้รับการแนะนำโดยLinus Paulingในปี 1935 เพื่ออธิบายเอนโทรปีที่เหลืออยู่ของน้ำแข็ง^{[ 1 ]}มีการเสนอรูปแบบต่างๆ เป็นโมเดลของผลึกเฟอร์โรอิเล็กทริก^{[ 2 ]}และแอนติเฟอร์โรอิเล็กทริก^{[ 3 ] บางชนิด}

ในปี พ.ศ. 2510 Elliott H. Liebได้ค้นพบวิธีแก้ปัญหาที่แน่นอนสำหรับแบบจำลองน้ำแข็งสองมิติที่เรียกว่า "น้ำแข็งสี่เหลี่ยม" ^{[ 4 ]}วิธีแก้ปัญหาที่แน่นอนในสามมิติเป็นที่ทราบกันเฉพาะสำหรับสถานะ "แช่แข็ง" พิเศษเท่านั้น^{[ 5 ]}

แบบจำลองประเภทน้ำแข็งเป็นแบบจำลองแลตทิซที่กำหนดบนแลตทิซที่มีเลขการประสานงาน 4 นั่นคือ จุดยอดแต่ละจุดของแลตทิซเชื่อมต่อด้วยขอบกับ "เพื่อนบ้านที่ใกล้ที่สุด" สี่จุด สถานะของแบบจำลองประกอบด้วยลูกศรบนขอบแต่ละด้านของแลตทิซ โดยที่จำนวนลูกศรที่ชี้เข้าด้านในที่จุดยอดแต่ละจุดมี 2 ข้อจำกัดนี้เกี่ยวกับการกำหนดค่าลูกศรเรียกว่ากฎน้ำแข็งใน แง่ของ ทฤษฎีกราฟสถานะต่างๆ คือการวางแนวแบบออยเลอร์ ของกราฟแบบไม่มีทิศทาง ปกติ 4 พื้นฐานฟังก์ชันพาร์ติชันยังนับจำนวนการไหล 3 ที่ไม่มีที่ใดเป็นศูนย์ด้วย^{[ 6 ]}

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

ณ จุดยอดใดๆ จะมีการจัดเรียงลูกศรได้หกแบบที่สอดคล้องกับกฎไอซ์ (ซึ่งเป็นที่มาของชื่อ "แบบจำลองหกจุดยอด") การจัดเรียงที่ถูกต้องสำหรับโครงตาข่ายสี่เหลี่ยม (สองมิติ) มีดังต่อไปนี้:

พลังงานของสถานะหนึ่งๆ นั้นเข้าใจกันว่าเป็นฟังก์ชันของการจัดเรียงตัวที่แต่ละจุดยอด สำหรับโครงตาข่ายสี่เหลี่ยมจัตุรัส จะถือว่าพลังงานรวมทั้งหมด${\displaystyle E}$ได้รับจาก

${\displaystyle E=n_{1}\epsilon _{1}+n_{2}\epsilon _{2}+\ldots +n_{6}\epsilon _{6},}$

สำหรับค่าคงที่บางค่า${\displaystyle \epsilon _{1},\ldots ,\epsilon _{6}}$, ที่ไหน${\displaystyle n_{i}}$ในที่นี้ หมายถึงจำนวนจุดยอดที่มี${\displaystyle i}$การกำหนดค่าจากรูปด้านบน ค่า${\displaystyle \epsilon _{i}}$คือพลังงานที่เกี่ยวข้องกับหมายเลขการกำหนดค่าจุดยอด${\displaystyle i}$.

มีเป้าหมายเพื่อคำนวณฟังก์ชันการแบ่งส่วน${\displaystyle Z}$ของแบบจำลองประเภทน้ำแข็ง ซึ่งกำหนดโดยสูตร

${\displaystyle Z=\sum \exp(-E/k_{\rm {B}}T),}$

โดยผลรวมนั้นคำนวณจากทุกสถานะของแบบจำลอง${\displaystyle E}$คือพลังงานของรัฐ${\displaystyle k_{\rm {B}}}$คือค่าคงที่ของโบลต์ซมันน์และ${\displaystyle T}$คืออุณหภูมิของระบบ

โดยทั่วไปแล้ว เรามักสนใจขีดจำกัดทางเทอร์โมไดนามิกส์ซึ่งจำนวนนั้น${\displaystyle N}$เมื่อจำนวนจุดยอดเข้าใกล้ค่าอนันต์ ในกรณีนั้น จะต้องประเมินพลังงานอิสระต่อจุดยอด แทน${\displaystyle f}$ในขีดจำกัดเมื่อ${\displaystyle N\to \infty }$, ที่ไหน${\displaystyle f}$ได้รับจาก

${\displaystyle f=-k_{\rm {B}}TN^{-1}\log Z.}$

ในทำนองเดียวกัน เราจะประเมินฟังก์ชันการแบ่งส่วนต่อจุดยอดแต่ละจุด${\displaystyle W}$ในขีดจำกัดทางเทอร์โมไดนามิกส์ ซึ่ง

${\displaystyle W=Z^{1/N}.}$

ค่านิยม${\displaystyle f}$และ${\displaystyle W}$มีความสัมพันธ์กันโดย

${\displaystyle f=-k_{\rm {B}}T\log W.}$

ผลึกจริงหลายชนิดที่มีพันธะไฮโดรเจนเป็นไปตามแบบจำลองน้ำแข็ง รวมถึงน้ำแข็ง^{[ 1 ]}และโพแทสเซียมไดไฮโดรเจนฟอสเฟตKH PO^{[ 2 ]} (KDP) อันที่จริง ผลึกดังกล่าวเป็นแรงบันดาลใจให้เกิดการศึกษาแบบจำลองประเภทน้ำแข็ง

ในน้ำแข็ง อะตอมออกซิเจนแต่ละอะตอมจะเชื่อมต่อด้วยพันธะกับไฮโดรเจนสี่อะตอม และแต่ละพันธะจะมีอะตอมไฮโดรเจนหนึ่งอะตอมอยู่ระหว่างออกซิเจนปลายสุด ไฮโดรเจนจะอยู่ในตำแหน่งสมมาตรสองตำแหน่ง ซึ่งไม่มีตำแหน่งใดอยู่ตรงกลางของพันธะ Pauling โต้แย้ง^{[ 1 ]}ว่าการจัดเรียงตัวของอะตอมไฮโดรเจนที่อนุญาตนั้นทำให้มีไฮโดรเจนสองอะตอมอยู่ใกล้กับออกซิเจนแต่ละอะตอมเสมอ ทำให้สภาพแวดล้อมในบริเวณนั้นเลียนแบบโมเลกุลของน้ำH Oดังนั้น ถ้าเราถือว่าอะตอมออกซิเจนเป็นจุดยอดของโครงผลึก และพันธะไฮโดรเจนเป็นขอบของโครงผลึก และถ้าเราลากลูกศรบนพันธะที่ชี้ไปยังด้านที่อะตอมไฮโดรเจนอยู่ น้ำแข็งก็จะสอดคล้องกับแบบจำลองน้ำแข็ง เหตุผลที่คล้ายกันนี้สามารถนำมาใช้เพื่อแสดงว่า KDP ก็สอดคล้องกับแบบจำลองน้ำแข็งเช่นกัน

ในช่วงไม่กี่ปีที่ผ่านมา มีการสำรวจแบบจำลองประเภทน้ำแข็งเพื่อใช้อธิบายสปินน้ำแข็ง ไพโรคลอร์ ^{[ 7 ]}และระบบสปินน้ำแข็งเทียม^{[ 8 ] [ 9 ]}ซึ่งความขัดแย้งทางเรขาคณิตในการปฏิสัมพันธ์ระหว่างโมเมนต์แม่เหล็ก แบบสองสถานะ ("สปิน") นำไปสู่การเลือกการกำหนดค่าสปินแบบ "กฎน้ำแข็ง" เมื่อไม่นานมานี้ การเปรียบเทียบดังกล่าวได้รับการขยายออกไปเพื่อสำรวจสถานการณ์ที่ระบบสปินน้ำแข็งอาจได้รับการอธิบายอย่างแม่นยำโดยแบบจำลอง Rys F ^{[ 10 ] [ 11 ] [ 12 ] [ 13 ]}

บนโครงตาข่ายสี่เหลี่ยม พลังงาน${\displaystyle \epsilon _{1},\ldots ,\epsilon _{6}}$ค่าพลังงานจุดยอดที่เกี่ยวข้องกับการกำหนดค่าจุดยอด 1-6 จะกำหนดความน่าจะเป็นสัมพัทธ์ของสถานะต่างๆ และด้วยเหตุนี้จึงสามารถส่งผลต่อพฤติกรรมระดับมหภาคของระบบได้ ต่อไปนี้คือตัวเลือกทั่วไปสำหรับพลังงานจุดยอดเหล่านี้

กฎน้ำแข็งได้รับการแนะนำโดย Linus Pauling ในปี พ.ศ. 2478 เพื่ออธิบายเอนโทรปีที่เหลืออยู่ของน้ำแข็งที่วัดโดยWilliam F. Giauqueและ JW Stout ^{[ 14 ]}เอนโทรปีที่เหลืออยู่${\displaystyle S}$ปริมาณน้ำแข็งคำนวณได้จากสูตร

${\displaystyle S=k_{\rm {B}}\log Z=k_{\rm {B}}\,N\,\log W,}$

ที่ไหน${\displaystyle k_{\rm {B}}}$คือค่าคงที่ของโบลต์ซมันน์${\displaystyle N}$คือจำนวนอะตอมออกซิเจนในก้อนน้ำแข็ง ซึ่งโดยทั่วไปถือว่ามีค่ามาก ( ขีดจำกัดทางเทอร์โมไดนามิก ) และ${\displaystyle Z=W^{N}}$คือจำนวนการจัดเรียงตัวของอะตอมไฮโดรเจนตามกฎน้ำแข็งของพอลลิง หากไม่มีกฎน้ำแข็ง เราจะมี${\displaystyle W=4}$เนื่องจากจำนวนอะตอมไฮโดรเจนคือ${\displaystyle 2N}$และไฮโดรเจนแต่ละอะตอมมีตำแหน่งที่เป็นไปได้สองตำแหน่ง พอลลิงประเมินว่ากฎน้ำแข็งจะลดสิ่งนี้ลงเหลือ${\displaystyle W=1.5}$ซึ่งเป็นตัวเลขที่สอดคล้องกับค่าที่วัดได้จาก Giauque-Stout เป็นอย่างดี${\displaystyle S}$อาจกล่าวได้ว่าการคำนวณของพอลลิงนั้น${\displaystyle S}$การใช้ กลศาสตร์เชิงสถิติกับน้ำแข็งเป็นหนึ่งในวิธีการประยุกต์ใช้ กลศาสตร์เชิงสถิติ กับสารจริงที่ง่ายที่สุดแต่แม่นยำที่สุดเท่าที่เคยมีมา คำถามที่ยังคงอยู่คือ เมื่อพิจารณาจากแบบจำลองแล้ว การคำนวณของพอลลิงเกี่ยวกับน้ำแข็งนั้นถูกต้องหรือไม่${\displaystyle W}$ซึ่งเป็นการประมาณค่าอย่างคร่าวๆ จะต้องได้รับการยืนยันด้วยการคำนวณอย่างเข้มงวด นี่กลายเป็นปัญหาสำคัญใน วิชา คณิตศาสตร์เชิงการจัดเรียง

แบบจำลองทั้งสามมิติและสองมิติได้รับการคำนวณเชิงตัวเลขโดย John F. Nagle ในปี พ.ศ. 2509 ^{[ 15 ]}ซึ่งพบว่า${\displaystyle W=1.50685\pm 0.00015}$ในสามมิติและ${\displaystyle W=1.540\pm 0.001}$ในสองมิติ ทั้งสองค่าใกล้เคียงกับการคำนวณคร่าวๆ ของพอลลิงอย่างน่าทึ่ง ซึ่งก็คือ 1.5

ในปี พ.ศ. 2510 Lieb พบวิธีแก้ปัญหาที่แม่นยำของแบบจำลองน้ำแข็งสองมิติสามแบบ ได้แก่ แบบจำลองน้ำแข็ง^{[ 4 ]}และ Rys${\displaystyle F}$แบบจำลอง^{[ 16 ]}และแบบจำลอง KDP ^{[ 17 ]}วิธีแก้ปัญหาสำหรับแบบจำลองน้ำแข็งให้ค่าที่แน่นอนของ${\displaystyle W}$ในสองมิติเช่น

${\displaystyle W_{2D}=\left({\frac {4}{3}}\right)^{3/2}=1.5396007....}$

ซึ่งรู้จักกันในชื่อ ค่าคง ที่กำลังสองของน้ำแข็งของ Lieb

ต่อมาในปี พ.ศ. 2510 ที. บิล ซัทเธอร์แลนด์ได้ขยายวิธีแก้ปัญหาของ Lieb สำหรับแบบจำลองน้ำแข็งสามแบบเฉพาะเจาะจงไปสู่วิธีแก้ปัญหาที่แม่นยำทั่วไปสำหรับแบบจำลองน้ำแข็งแบบตารางสี่เหลี่ยมที่สอดคล้องกับสมมติฐานสนามศูนย์^{[ 18 ]}

ต่อมาในปี พ.ศ. 2510 CP Yang ^{[ 19 ]}ได้ขยายวิธีแก้ปัญหาของ Sutherland ไปสู่วิธีแก้ปัญหาที่แม่นยำสำหรับแบบจำลองน้ำแข็งแบบตารางสี่เหลี่ยมในสนามไฟฟ้าแนวนอน

ในปี พ.ศ. 2512 จอห์น เนเกิล ได้หาคำตอบที่แน่นอนสำหรับแบบจำลอง KDP เวอร์ชันสามมิติ สำหรับช่วงอุณหภูมิที่เฉพาะเจาะจง^{[ 5 ]}สำหรับอุณหภูมิดังกล่าว แบบจำลองจะ "หยุดนิ่ง" ในแง่ที่ว่า (ในขีดจำกัดทางเทอร์โมไดนามิก) พลังงานต่อจุดยอดและเอนโทรปีต่อจุดยอดเป็นศูนย์ทั้งคู่ นี่เป็นคำตอบที่แน่นอนเพียงคำตอบเดียวที่ทราบสำหรับแบบจำลองประเภทน้ำแข็งสามมิติ

แบบจำลองแปดจุดยอดซึ่งได้รับการแก้ไขอย่างแม่นยำแล้วนั้น เป็นการวางนัยทั่วไปของแบบจำลองหกจุดยอด (แบบตารางสี่เหลี่ยม) โดยในการกู้คืนแบบจำลองหกจุดยอดจากแบบจำลองแปดจุดยอด ให้ตั้งค่าพลังงานสำหรับการกำหนดค่าจุดยอด 7 และ 8 เป็นอนันต์ แบบจำลองหกจุดยอดได้รับการแก้ไขในบางกรณีที่แบบจำลองแปดจุดยอดไม่สามารถแก้ไขได้ ตัวอย่างเช่น วิธีแก้ปัญหาของ Nagle สำหรับแบบจำลอง KDP สามมิติ^{[ 5 ]}และวิธีแก้ปัญหาของ Yang สำหรับแบบจำลองหกจุดยอดในสนามแนวนอน^{[ 19 ]}

แบบจำลองน้ำแข็งนี้ให้ 'ตัวอย่างคัดค้าน' ที่สำคัญในกลศาสตร์เชิงสถิติ: พลังงานอิสระจำนวนมากในขีดจำกัดทางเทอร์โมไดนามิก ขึ้นอยู่กับเงื่อนไขขอบเขต^{[ 20 ]}แบบจำลองนี้ได้รับการแก้ไขเชิงวิเคราะห์สำหรับเงื่อนไขขอบเขตแบบคาบ แบบแอนตี้คาบ แบบเฟอร์โรแมกเนติก และเงื่อนไขขอบเขตผนังโดเมน แบบจำลองหกจุดยอดที่มีเงื่อนไขขอบเขตผนังโดเมนบนโครงตาข่ายสี่เหลี่ยมมีความสำคัญเฉพาะในคณิตศาสตร์เชิงการจัดเรียง ช่วยในการแจงนับเมทริกซ์เครื่องหมายสลับในกรณีนี้ ฟังก์ชันพาร์ติชันสามารถแสดงเป็นดีเทอร์มิแนนต์ของเมทริกซ์ (ซึ่งมีมิติเท่ากับขนาดของโครงตาข่าย) แต่ในกรณีอื่น ๆ การแจงนับของ${\displaystyle W}$มันไม่ได้ออกมาในรูปแบบที่เรียบง่ายและปิดสนิทเช่นนั้น

เห็นได้ชัดว่าใหญ่ที่สุด${\displaystyle W}$กำหนดโดย เงื่อนไขขอบเขต อิสระ (ไม่มีข้อจำกัดใดๆ เกี่ยวกับการกำหนดค่าบนขอบเขต) แต่เหมือนกัน${\displaystyle W}$เกิดขึ้นในขีดจำกัดทางเทอร์โมไดนามิกสำหรับเงื่อนไขขอบเขตแบบคาบ^{[ 21 ]}ดังที่ใช้ในตอนแรกเพื่อหาอนุพันธ์${\displaystyle W_{2D}}$.

จำนวนสถานะของแบบจำลองประเภทน้ำแข็งบนขอบภายในของการรวมกันของสี่เหลี่ยมจัตุรัสที่เชื่อมต่อกันอย่างง่ายแบบจำกัดบนโครงข่ายนั้น เท่ากับหนึ่งในสามของจำนวนวิธีในการระบายสีสี่เหลี่ยมจัตุรัส 3 สี โดยที่ไม่มีสี่เหลี่ยมจัตุรัสที่อยู่ติดกันสองรูปใดมีสีเดียวกัน ความสัมพันธ์ระหว่างสถานะนี้เป็นผลงานของแอนดรูว์ เลนาร์ด และแสดงได้ดังนี้ หากสี่เหลี่ยมจัตุรัสมีสีi = 0, 1 หรือ 2 ลูกศรบนขอบไปยังสี่เหลี่ยมจัตุรัสที่อยู่ติดกันจะชี้ไปทางซ้ายหรือขวา (ตามผู้สังเกตในสี่เหลี่ยมจัตุรัส) ขึ้นอยู่กับว่าสีในสี่เหลี่ยมจัตุรัสที่อยู่ติดกันคือi + 1 หรือi − 1 mod 3 มี 3 วิธีที่เป็นไปได้ในการระบายสีสี่เหลี่ยมจัตุรัสเริ่มต้นที่กำหนดไว้ และเมื่อเลือกสีเริ่มต้นแล้ว จะได้ความสัมพันธ์แบบ 1:1 ระหว่างการระบายสีและการจัดเรียงลูกศรที่ตรงตามเงื่อนไขประเภทน้ำแข็ง

• แบบจำลองแปดจุดยอด

• Lieb, EH; Wu, FY (1972), "แบบจำลองเฟอร์โรอิเล็กทริกสองมิติ", ใน C. Domb; MS Green (บรรณาธิการ), การเปลี่ยนเฟสและปรากฏการณ์วิกฤต , เล่ม 1, นิวยอร์ก: Academic Press, หน้า331–490 
• Baxter, Rodney J. (1982), แบบจำลองที่แก้ได้อย่างแม่นยำในกลศาสตร์เชิงสถิติ (PDF) , ลอนดอน: Academic Press Inc. [Harcourt Brace Jovanovich Publishers], ISBN 978-0-12-083180-7MR 0690578เก็บถาวรจากต้นฉบับ(PDF)เมื่อ 2021-04-14 เรียกดูเมื่อ2012-08-12