توابع انرژی پتانسیل دقیقی برای بخارات فلزات قلیایی در دسترس نیست. به دلیل محدوده وسیع دمایی مایع در فلزات قلیایی مذاب، هیچ تابع انرژی پتانسیل تحلیلی، که قادر باشد خواص فلزات قلیایی مایع را پیشبینی کند شناخته نشده است. در این کار توابع انرژی پتانسیل برای سزیم با استفاده از روش تکرار شونده بولتزمن محاسبه شده است. در این روش ساختار مایع مطابق با ساختار معلوم تجربی، در یک سیکل تکرار شونده تنظیم می شود. در انجام شبیه سازی دینامیک مولکولی با یک پتانسیل حدسی اولیه، ساختار محاسبه شده و سپس با ساختار تجربی مقایسه می شود. سپس تابع انرژی پتانسیل نیز از همین طریق تنظیم می شود، تا بتواند ساختارهای قابل مقایسه با نوع تجربی متناظر بدست آورد. در این مطالعه منحنی عددی انرژی پتانسیل برای سزیم را در محدوده دماییk 1873-573 محاسبه کرده ایم. به نظر می رسد که منحنی های انرژی پتانسیل ارائه شده، قادراند ساختارهای قابل قیاس با نمونه های تجربی را پیشبینی کنند.