前沿拓展：

win10裝lammps

黑屏的原因和處理方法：此問題通常是由顯卡驅(qū)動(dòng)程序（雙顯卡）的兼容性問題引起的?？梢試L試關(guān)閉電源并多次重新啟動(dòng)計(jì)算機(jī)?？梢栽?/span>短時(shí)間內(nèi)重復(fù)三次后進(jìn)入winre（Windows恢復(fù)環(huán)境），但有些用戶可能需要直接斷開電源。

一、啟動(dòng)后會(huì)出法解決問題，可以嘗試選擇返回到windows疑難解答中的上一個(gè)系統(tǒng)還原點(diǎn)，以便系統(tǒng)可以返回到上一個(gè)兼容性問題。

跟據(jù)過度態(tài)理論，體系中某一group原子的擴(kuò)散系數(shù)與均方位移具有以下公式（1）和（2）的關(guān)系，公式（1）中r(t)為原子位移，d為維度，t為測(cè)量原子擴(kuò)散的時(shí)間。公式（3）為擴(kuò)散系數(shù)與溫度T的關(guān)系，其中Ea為擴(kuò)散能壘，kB為波爾滋蔓常數(shù)。因此只需要計(jì)算一定溫度范圍內(nèi)擴(kuò)散系數(shù)與溫度的關(guān)系，就可以通過公式（3）擬合擴(kuò)散能壘。

本文筆者將利用LAMMPS以H在Si@SiO2異質(zhì)結(jié)處擴(kuò)散為例，示范通過過渡態(tài)求300K，H擴(kuò)散系數(shù)的過程。同時(shí)筆者會(huì)給出in文件，讀者可經(jīng)過簡(jiǎn)單的修改，求出一定溫度范圍的擴(kuò)散系數(shù)與溫度的數(shù)值關(guān)系從而擬合擴(kuò)散能壘。

in文件腳本如下：

### Basic setup ###variable t equal 300units metaldimension 3boundary p p patom_style atomicnewton on### Structure ###read_data pos1.datneighbor 1 bincreate_atoms 3 single 5 5 22### Potentials ###pair_style hybrid tersoff lj/cut 10.0pair_coeff * * tersoff SiO.tersoff Si O pair_coeff 1 3 lj/cut 0.0057673 3.1988pair_coeff 2 3 lj/cut 0.0022281 2.8446pair_coeff 3 3 lj/cut 0.001908 2.5711mass 1 28.085mass 2 16mass 3 1.008### Group ###group H type 3### Minimization ###thermo 100thermo_style custom step time vol atoms temp ke pe lx pressfix 1 all box/relax iso 0.0 vmax 0.01min_style cgminimize 1e-25 1e-25 100000 100000### NPT1 ###reset_timestep 0thermo 1000thermo_style custom step dt time atoms temp ke pe vol press lx ly lz velocity all create $t 92137453 dist gaussianfix NPT all npt temp $t $t 1 iso 0.0 0.0 1.0 drag 0.2timestep 0.001run 10000reset_timestep 0compute Msd H msd fix s all vector 1 c_Msd[4]variable MSD equal sum(f_s)/1000 variable stp equal step*dt### NPT2 ###fix msd1 all print 1000 “${stp} ${MSD}” file msd_$t.txt screen nothermo 1000thermo_style custom step temp ke pe vol press lx v_MSDtimestep 0.001run 500000unfix NPT

運(yùn)行完畢后通過對(duì)計(jì)算所得MSD與時(shí)間得關(guān)系進(jìn)msd_300.txt文件拖入Oringlab行擬合，可得300 K下H在Si@SiO2異質(zhì)結(jié)處擴(kuò)散系數(shù)為8.295 ?2/ps=8.295*108 m2/s。

————————————————————–

PS：由于時(shí)間倉(cāng)促，可能有錯(cuò)誤之處，歡迎各位同行交流學(xué)習(xí)，同時(shí)若有疑問、計(jì)算技術(shù)討論、合作可以聯(lián)系筆者郵箱：1924311399@qq.com謝謝！

**推薦：計(jì)算運(yùn)維鳥

拓展知識(shí)：

win10裝lammps

踏進(jìn)材料科學(xué)的大門的我，需要先安裝一些模擬軟件：在win10系統(tǒng)中，安裝虛擬機(jī)virtualbox。在虛擬機(jī)中運(yùn)行CentOS系統(tǒng)，在CentOS系統(tǒng)中安裝fftw，mpich，Lammps。

在win10中在virtualbox中安裝CentOS的過程—— https://www.jianshu.com/p/e1791a8abb2d
linux指令學(xué)習(xí)—— https://www.runoob.com/linux/linux-file-content-manage.html
yum指令是什么—— https://wangchujiang.com/linux-command/c/yum.html
centos7 vim/vi的使用方法—— https://blog.csdn.net/capecape/article/details/78503497
centos下vi 命令用法—— https://www.cnblogs.com/winfosy/p/9478303.html
fftw是什么—— https://baike.baidu.com/item/FFTW/9974883?fr=aladdin
mpich是什么—— https://baike.baidu.com/item/MPICH/7488372?fr=aladdin
centos 下make命令的使用—— https://blog.csdn.net/programercch/article/details/56480168
Virtualbox主機(jī)和虛擬機(jī)之間文件夾共享及雙向拷貝—— https://blog.csdn.net/zhaipengfei1231/article/details/78471038?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-1.channel_param&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-1.channel_param
VirtualBox中如何實(shí)現(xiàn)主機(jī)和虛擬機(jī)之間的文件拖放（Drag-and-Drop）—— https://www.crifan.com/summary_virtualbox_how_to_drap_and_drop_file_between_host_and_guest_vm/
Linux wget 命令用法詳解—— https://www.jianshu.com/p/59bb131bc2ab
mpich安裝（前面的準(zhǔn)備工作千萬(wàn)要先做好?。?https://blog.csdn.net/sinat_30967935/article/details/82988659
CentOS 7 安裝gcc、g++、gdb教程—— https://blog.csdn.net/alan_gaohaodong/article/details/79656177
在CentOS下安裝并行l(wèi)ammps教程—— https://tieba.baidu.com/p/5883557632?red_tag=3335137741
CentOS普通用戶獲得管理員權(quán)限—— https://blog.csdn.net/sy_bz/article/details/52399049

在這里衷心地感謝所有善良熱心無(wú)私分享各種解決方案的作者！各位都是人間天使！