這裡使用tip4p這種水計算其靜止介電係數與頻率依賴的介電常數,相關輸入檔案使用version 2018.1的gromacs自帶的tip4p.gro和tip4p.itp檔案。
在linux終端中輸入
gmx solvate -cs tip4p.gro -o conf.gro -box 4.11
4.54.9
-p topol.top
這裡這個盒子的尺寸是正好有3000個tip4p水分子
如果是要重複文獻中的實驗,則需要根據文獻修改以下過程中的.mdp檔案引數。特別注意約束方式和fs,還有就是模擬時間,以及最小鍵程截斷,這些都會影響最終的結果。詳細過程可以參考我的另一篇文章
gmx grompp -f minim.mdp -c conf.gro -p topol.top -o em.tpr
gmx mdrun -v -deffnm em
gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -o nvt.tpr
gmx mdrun -v -deffnm nvt
gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -p topol.top -o npt.tpr
gmx mdrun -v -deffnm npt
gmx dipoles -corr total -c -f npt.trr -s npt.tpr
gmx dielectric -f dipcorr.xvg -eps0 53.2129
-ffn aexp
這裡的eps0是偶極距計算中的靜止介電常數
xmgrace epsw.xvg
用grace顯示 mysql怎麼計算冪 在MySQL中計算數字的冪?
要計算數字的冪,請使用power 函式。讓我們首先建立乙個表 create table demotable amount int 使用插入命令在表中插入一些記錄 insert into demotable values 64 insert into demotable values 144 inse...
在Qt中計算文字的寬度與高度
在寫一些需要精確布局或者排版的軟體時候,我們需要知道文字的寬度和高度,保證我們在適當的座標中使用drawtext。例如,我在做乙個簡單的掃雷軟體時,希望把每個格仔附近的雷數給顯示在小格仔裡,如下圖所示,就需要這樣的功能。不同的字型,字寬字高都是不一樣的。所以計算文字寬度的是乙個和字型相關的乙個類qf...
在儲存過程中計算考核總分
alter procedure dbo zhkh sp checkcompute kcode char 10 種類 kjyear int,會計年度 period int 會計月度 asdeclare type int 型別 0 部門,1 員工 declare mcode char 10 模組編碼 d...