馬pi精

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 32.4
帖子: 69
在線: 17.9小時(shí)
蟲號(hào): 1417222

[交流] Ising Model 蒙特卡洛模擬比熱和磁化率曲線

各位寫過Ising Model模擬算法的大手們~大家畫出來的比熱和磁化率曲線和理論曲線相符嗎?
在下用C寫了一個(gè)算法, 用的是Metropolice算法, 但是畫出來的比熱和磁化率的曲線在溫度0~1的時(shí)候有一個(gè)小峰, 如圖所示. 各位請(qǐng)忽略峰位不太對(duì), 這個(gè)圖不是方格子的, 方格子的圖一時(shí)沒照出來, 但是問題是一樣的...不同曲線是加了不同的外磁場.
苦思許久沒有找出程序的錯(cuò)誤, 竊以為是能量和磁化強(qiáng)度的漲落在低溫時(shí)都很小, 相近的溫度漲落值是很接近的, 從比熱的計(jì)算公式
Cv~(-^2)/T^2可見, 分子很接近, 分母越小, Cv越大....T-T 磁化率也有相似的情況, 但是在低溫時(shí)系統(tǒng)不應(yīng)該有很大的比熱(代表漲落)T^T
不知道大家遇到這個(gè)情況呢, 如果遇到了大家是怎么解決的? 還是只能是因?yàn)槲业某绦虺鲥e(cuò)了?
謝謝大家了~~

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

文集

» 猜你喜歡

276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷已經(jīng)有8人回復(fù)
接收2026碩士調(diào)劑(學(xué)碩+專碩) 已經(jīng)有3人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有3人回復(fù)
081700 調(diào)劑 267分已經(jīng)有3人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有3人回復(fù)
336化工調(diào)劑已經(jīng)有4人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有10人回復(fù)
263求調(diào)劑已經(jīng)有6人回復(fù)
298求調(diào)劑已經(jīng)有8人回復(fù)
考研化學(xué)308分求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

2273測試鈦合金的循環(huán)伏安曲線沒有峰已經(jīng)有7人回復(fù)
關(guān)于Ising模型的蒙特卡洛模擬已經(jīng)有6人回復(fù)
隨行標(biāo)準(zhǔn)曲線和現(xiàn)配標(biāo)準(zhǔn)曲線意思是一個(gè)樣的么�。。。� 已經(jīng)有8人回復(fù)
想擬合一條簡單曲線，請(qǐng)教一下大家~ 已經(jīng)有6人回復(fù)
菌種生長曲線的差異比較，用何種統(tǒng)計(jì)學(xué)方法？謝謝！已經(jīng)有3人回復(fù)
求助��！求一程序，用matlab程序做，用蒙特卡洛方法模擬已經(jīng)有10人回復(fù)
怎么用Origin軟件做三相圖的擬合曲線？已經(jīng)有5人回復(fù)
鋰離子電池充放電曲線的困惑已經(jīng)有5人回復(fù)
【討論】CV曲線贗電容相關(guān) 已經(jīng)有23人回復(fù)
【求助】如何使用gromacs中的constraint 已經(jīng)有8人回復(fù)
【求助】水蒸氣的比熱容是多少已經(jīng)有6人回復(fù)
【求助】Ising model 二維正方格子的磁化強(qiáng)度計(jì)算已經(jīng)有20人回復(fù)
【討論】關(guān)于不同z值處的水能形成氫鍵數(shù)的腳本討論已經(jīng)有4人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2012-08-07 00:35:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hrandy

金蟲 (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 947.3
帖子: 342
在線: 186.7小時(shí)
蟲號(hào): 577250

★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
華麗的飄過: 金幣+2, 3q 2012-10-15 01:07:49

你采用初態(tài)為高溫?zé)o序結(jié)構(gòu)，然后從高溫降到低溫試試。如果以零溫?zé)o序?yàn)槌鯌B(tài)，計(jì)算開始不久就會(huì)形成磁疇，這種疇一旦形成，低溫下很難改變疇內(nèi)自旋方向。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

馬pi精

2樓2012-08-07 16:54:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

馬pi精

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 32.4
帖子: 69
在線: 17.9小時(shí)
蟲號(hào): 1417222

送鮮花一朵

引用回帖:

2樓: Originally posted by hrandy at 2012-08-07 16:54:23
你采用初態(tài)為高溫?zé)o序結(jié)構(gòu)，然后從高溫降到低溫試試。如果以零溫?zé)o序?yàn)槌鯌B(tài)，計(jì)算開始不久就會(huì)形成磁疇，這種疇一旦形成，低溫下很難改變疇內(nèi)自旋方向。

謝謝回復(fù)~這個(gè)圖已經(jīng)是按照你這么說的做了呀...而且相鄰兩個(gè)溫度相差很小, 這才讓我糾結(jié)...

贊一下

回復(fù)此樓

4樓2012-08-07 22:45:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mywtfmp3

銀蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 242.7
帖子: 71
在線: 30.1小時(shí)
蟲號(hào): 1021073

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

關(guān)注一下樓主，我是新手

贊一下

回復(fù)此樓

6樓2012-08-09 00:54:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

馬pi精

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 32.4
帖子: 69
在線: 17.9小時(shí)
蟲號(hào): 1417222

自己頂一下...跪求各位大手指點(diǎn)啊T-T

贊一下

回復(fù)此樓

7樓2012-08-17 23:06:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhuanhang

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 2
在線: 6.6小時(shí)
蟲號(hào): 1619389

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

0溫時(shí)，應(yīng)該不會(huì)有相變，以我的理解，你應(yīng)該貼code

贊一下

回復(fù)此樓

8樓2012-10-14 17:09:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

蒙特卡羅

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 1207
帖子: 669
在線: 148.4小時(shí)
蟲號(hào): 3034532

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

內(nèi)容已刪除

回復(fù)此樓

9樓2014-09-19 16:05:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhongwei2284

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 14.8
帖子: 80
在線: 81.7小時(shí)
蟲號(hào): 1820687

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

這個(gè)有可能只是一種計(jì)算中出現(xiàn)的邊界效應(yīng)。而且你的這個(gè)曲線的臨界值出現(xiàn)在4左右了，實(shí)際上的不是2.3左右嗎？用的不一樣的公式？方法？

贊一下

回復(fù)此樓

10樓2014-09-22 22:02:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

鄭克寧3樓

2012-08-07 17:18 回復(fù)

華麗的飄過: 請(qǐng)勿純表情回復(fù)，多少說點(diǎn)什么吧 2012-10-15 01:07:32

strive1231235樓

2012-08-08 09:15 回復(fù)

華麗的飄過: 請(qǐng)勿純表情回復(fù) 2012-10-15 01:07:16

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	4/200	2026-03-23 09:59 by bingchuan
[考研] 291求調(diào)劑 +5	孅華 2026-03-22	5/250	2026-03-23 09:20 by haoshis
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4	幸運(yùn)的醬醬 2026-03-22	5/250	2026-03-22 20:15 by edmund7
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 化學(xué)調(diào)劑 +5	yzysaa 2026-03-21	5/250	2026-03-21 22:12 by peike
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學(xué)員8dgXkO
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[基金申請(qǐng)] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客