av影院在线观看-久久久久久国产精品视频-亚洲综合第一页-日日夜夜天天-快射视频网-国产精品国产-欧美三级午夜理伦三级中视频-国产日韩未满十八禁止观看-亚洲一区国产-亚洲成av人片-久草综合视频-天天射影院-91精品一区二区-国产又粗又硬-黄色91免费-亚洲精品无码永久在线观看-中文字幕人妻高清乱码

Document

當(dāng)前位置：文庫百科 ? 文章詳情

分子動力學(xué)(LAMMPS)——實(shí)例講解：NEMD方法計算晶格熱導(dǎo)率

來源：科學(xué)10分鐘時間：2021-05-18 17:58:03 瀏覽：13780次

熱導(dǎo)率又稱導(dǎo)熱系數(shù)，反映物質(zhì)的熱傳導(dǎo)能力，其定義為單位溫度梯度（在1m長度內(nèi)溫度降低1K）在單位時間內(nèi)經(jīng)單位導(dǎo)熱面所傳遞的熱量。材料的熱導(dǎo)率是由電子熱導(dǎo)率和晶格熱導(dǎo)率共同貢獻(xiàn)的。在溫度較高的情況下，材料的熱導(dǎo)率主要由電子熱導(dǎo)率主導(dǎo)；在溫度較低的情況，主要由晶格熱導(dǎo)率主導(dǎo)。

實(shí)驗中測材料熱導(dǎo)率是在材料兩端施加熱流，此時可以測得沿?zé)崃鞣较虻臋M截面上的熱流通量J以及溫度梯度?T ，而J和?T 具有公式（1）的線性關(guān)系，比例系數(shù)K 就是材料的熱導(dǎo)率。

J = -k ?T （1）

本文我們將使用由美國能源部桑迪亞國家實(shí)驗室設(shè)計開發(fā)的分子動學(xué)軟件包LAMMPS[1]復(fù)現(xiàn)文獻(xiàn)[2]中用分子動力學(xué)計算金屬鎢晶格熱導(dǎo)率的方法。文獻(xiàn)[2]采用的是非平衡分子動力學(xué)(NEMD) [3-6]方法計算晶格熱導(dǎo)率。計算過程中的模擬盒子示意圖如圖1所示，當(dāng)模擬盒子在某一溫度下平衡后，將模擬盒子中間作為熱端，施加熱流?ε；模擬盒子兩端作為冷端對熱流進(jìn)行抽取，在抽取的過程當(dāng)中體系中原子的動量是從熱端到冷端不斷交換的，此時體系處于一個非平衡的狀態(tài)。因此該過程也叫非平衡分子動力學(xué)(NEMD)過程。

在施加熱流后，將模擬盒子沿?zé)崃鞣较蚍殖稍S多層，每個特定的timesteps (如10000步)將每一層的溫度進(jìn)行統(tǒng)計，當(dāng)傳熱過程穩(wěn)定后 (這里是指統(tǒng)計的溫度保持基本穩(wěn)定，一般在加熱流數(shù)萬步到十萬步后)，對特定timesteps間隔統(tǒng)計的分層溫度做平均并進(jìn)行線性擬合得到如圖2所示的溫度梯度曲線（黑色）。得到溫度梯度后根據(jù)公式（2）得到晶格熱導(dǎo)率k。公式（2）中A為模擬盒子的橫截面積，M為熱交換的頻率，t為單位時間。最后，結(jié)合公式（1）和公式（2）得出金屬鎢的晶格熱導(dǎo)率。筆者取的?ε/t =1eV/ps，A=16?²（筆者這里設(shè)置的有些不合理但對計算結(jié)果影響不大），M =10。筆者計算得出W在300K時延<100>方向的熱導(dǎo)率是16.7 W/KM，和文獻(xiàn)的計算結(jié)果17.28034 W/KM基本一致[2]。

圖1. NEMD計算熱導(dǎo)率的三維模型示意圖

圖2. 模擬盒子分層溫度統(tǒng)計曲線(紅色),溫度梯度擬合曲線(黑色)

在計算過程當(dāng)中，LAMMPS的in文件中的傳熱過程設(shè)置腳本段如下：

unfix 1

fix 2 all nve

fix hot all heat 10 ${power} region hot

fix cold all heat 10 -${power} region cold

其中，unfix 1命令取消的是前面作特定溫度下平衡的系綜。其余三段為熱流為傳熱過程的設(shè)置。

計算結(jié)束后，最關(guān)鍵的是對模擬盒子每一層原子溫度的統(tǒng)計，為此我們寫了一個python小腳本，可以直接運(yùn)行python3 thermal_cond.py profile.txt 200 25 45 > result.txt輕松處理，python的用法，及函數(shù)的定義，腳本thermal_cond.py中都有注釋。處理后得到文本文件result.txt，該文件中的信息是模擬盒子每一層溫度與層中心到模擬盒子左端距離的關(guān)系。直接將result.txt拖到作圖軟件中，或用gnuplot等程序作圖即可獲得火山峰狀的溫度統(tǒng)計曲線。最后截取其中的直線段擬合即可得溫度梯度?T，本文的?T=-0.6K/?。

最后，介紹幾個可以使結(jié)果更優(yōu)的方法：

（1）將模擬盒子三個方向邊長，建成W單包三個方向邊長的倍數(shù)，這樣模擬盒子具有完美晶體的周期性。

（2）將采集溫度的模擬時間盡量增大，這樣能夠采集更多的樣本，統(tǒng)計的溫度就會更加準(zhǔn)確。

該算例以及溫度梯度處理腳本可以在網(wǎng)盤下載：

鏈接：https://pan.baidu.com/s/1zkpULJ2T12NFkgnQOt5KOA

提取碼：CSGO

參考文獻(xiàn)：

[1] Plimpton S 1995 J. Comput. Phys. 117 1

[2] Baoqin Fu,Wensheng Lai,Yue Yuan, et al. 2012 J. Nucl. Mater. 427 1-3

[3] P. K. Schelling, S.R. Phillpot, P. Keblinski. 2002 Phys. Rev. B 65 144306.

[4] A. Maiti, G.D. Mahan, S.T. Pantelides. 1997 Solid State Commun. 102 517.

[5] A. Maeda. 1995 Phys. Rev. E 52 34.

[6] A.R. Abramson, C.L. Tien, A. Majumdar. 2002 J. Heat Transfer 124 963.

文章精選

共價鍵vs離子鍵？ELF、DORI兩種方法玩轉(zhuǎn)化學(xué)鍵分析
【技能】Win10系統(tǒng)Bash中快速安裝和使用LAMMPS的方法
直播回放《Lammps分子動力學(xué)課程》第1-5節(jié)，讓你零基礎(chǔ)快速入門（附PPT課件）
物理VS化學(xué)吸附？RDG，讓你的文章與眾不同
基于Quantum Espresso的第一性原理計算建模教程（四）：模型的建立

你的點(diǎn)贊和在看，我都認(rèn)真當(dāng)成喜歡

上一篇：?根正苗紅，發(fā)文量大，影響因子高，這本據(jù)說只有大牛才可以灌水的期刊，十分不錯

下一篇：大牛指路丨崔屹、張強(qiáng)、樓雄文和麥立強(qiáng)最新成果速遞

相關(guān)文章

基礎(chǔ)理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結(jié)構(gòu)）

2020-05-03

手把手教你用ChemDraw 畫化學(xué)結(jié)構(gòu)式：基礎(chǔ)篇

2021-06-19

晶體結(jié)構(gòu)可視化軟件 VESTA使用教程（下篇）

2021-01-22

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（上）

2019-10-25

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（下）

2019-10-25

XRD的基本原理與應(yīng)用

2020-11-03

舉報有獎

TEL: 191-3608-6524

如：在網(wǎng)絡(luò)上惡意使用“測試狗”等相關(guān)關(guān)鍵詞誤導(dǎo)用戶點(diǎn)擊、惡意盜用測試狗商標(biāo)、冒稱官方工作人員等情形，請您向我們舉報，經(jīng)查實(shí)后，我們將給予您獎勵。

舉報內(nèi)容：

200

上傳附件：

文件格式不正確，請重新上傳文件格式不正確，請重新上傳文件格式不正確，請重新上傳

文件格式：jpg、jpeg、png、gif、tif、doc、docx、ppt、pptx、xls、xlsx、pdf、zip、rar

聯(lián)系方式

姓名

電話

提交意見

意見反饋

Suggestions

您可以在此留下您寶貴的意見，您的意見或問題反饋將會成為我們不斷改進(jìn)的動力。

意見類型

測試服務(wù)

網(wǎng)站功能

財務(wù)報賬

其他類型

意見內(nèi)容

200

聯(lián)系方式

姓名

電話

提交意見

+

你好，很高興為您服務(wù)！

發(fā)送

項目推薦

舉報有獎

意見反饋