 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2332 | 回復: 4 | |||
牛排小黑嬌新蟲 (初入文壇)
|
[求助]
勢函數(shù)解析問題
|
|
我想做Mo-Ag的合金擴散模擬，但是并沒有現(xiàn)成的勢函數(shù)文件，我查閱文獻查到了有人擬合過這個合金的勢函數(shù)。 勢函數(shù)形式是FS勢格式和EAM勢差不多的，我看lammps手冊這個勢函數(shù)可以通過一個DYNAMO的程序解析出來，不知道您是不是也做過這個勢函數(shù)的解析文件，我想讓您幫我看一下我的程序是不 是有問題，我之前弄出來的模擬結(jié)果只能運行一行。希望尚老師您能抽時間幫我看下，萬分感謝！ #include <stdio.h> #include <math.h> #include <stdlib.h> #include <unistd.h> #include <string> #include <stdarg.h> #include <string.h> #include "vector" using namespace std; #define EOK 0 #define ERROR -1 double phi(double r,double cc,double c0,double c1,double c2,double c3,double c4) { double V_r; if(r<=cc) { V_r= (r-cc)*(r-cc)*(c0+c1*r+c2*r*r+c3*pow(r,3)+c4*pow(r,4)); } else V_r=0.0; return (V_r); } double rho(double r,double d,double beta) { double rho_r; if(r<=d) rho_r= (r-d)*(r-d) + beta*beta*pow((r-d),4); else rho_r=0.0; return (rho_r); } double Frho(double rho,double AA) { double F_rho; F_rho = -1.0*AA*pow(rho,0.5); return (F_rho); } int main (void) { int Nr = 10000; double rmax =9.8; 這邊取值的標準是？ double dr = rmax/(double)Nr; int Nrho = 10000; double rhomax = 45.5; double drho = rhomax/(double)Nrho; int i; string LAMMPSFilename="AgMo.eam.fs"; FILE *LAMMPSFile = fopen ((char*)LAMMPSFilename.data(), "w"  ;if (!LAMMPSFile) exit (ERROR); // Header for setfl format fprintf (LAMMPSFile, \ "#-> LAMMPS Potential File in DYNAMO 86 setfl Format <-#\n"\ "# Title :Embedded-Atom Method by Arthur F. Woter 1994\n"\ "# Implemented by Nize(2015) nize19890627@gmail.com\n"\ "2 Ag Mo\n"\ "%d %20.20f %d %20.20f %20.20f\n",Nrho,drho,Nr,dr,rmax); fprintf (LAMMPSFile,"47 107.8682 4.0853 FCC\n" ;// Embedding function and density function for Ag for(i=0;i<Nrho;i++) fprintf (LAMMPSFile,"%20.20e\n",Frho((double)i*drho,0.325514)); // Density function for Ag for (i = 0; i < Nr; i++) fprintf (LAMMPSFile, "%20.20e\n",rho((double)i*dr,4.41,-1.293394)); fprintf (LAMMPSFile,"42 95.94 3.1472 BCC\n" ;// Embedding function and density function for Mo for(i=0;i<Nrho;i++) fprintf (LAMMPSFile,"%20.20e\n",Frho((double)i*drho,1.848648)); // Density function for Mo for (i = 0; i < Nr; i++) fprintf (LAMMPSFile, "%20.20e\n",rho((double)i*dr,4.14,0.0)); // Pair potential for Ag-Ag for (i = 0; i < Nr; i++) fprintf (LAMMPSFile, "%20.20e\n", (double)i*dr*phi((double)i*dr,4.76,10.681200,-12.045170,5.203072,-1.013304,0.0742308)); // Pair potential for Ag-Mo for (i = 0; i < Nr; i++) fprintf (LAMMPSFile, "%20.20e\n", (double)i*dr*phi((double)i*dr,4.50,44.406810,-45.490260,15.605110,-1.793704,0.0)); // Pair potential for Mo-Mo for (i = 0; i < Nr; i++) fprintf (LAMMPSFile, "%20.20e\n", (double)i*dr*phi((double)i*dr,3.2572,47.980660,-34.099240,5.832293,0.101749,0.0203934)); fclose (LAMMPSFile); // system("pause" ;return (EOK); 知道的朋友幫我看看。急求 謝謝啦 QQ圖片20150626105650.png QQ圖片20150626105951.png QQ圖片20150702114556.jpg |
回復此樓» 收錄本帖的淘帖專輯推薦
 分子動力學模擬 |
» 猜你喜歡
336求調(diào)劑
已經(jīng)有3人回復
-
306求調(diào)劑
已經(jīng)有9人回復
-
工科0856求調(diào)劑
已經(jīng)有4人回復
-
工科材料085601 279求調(diào)劑
已經(jīng)有10人回復
-
081700 調(diào)劑 267分
已經(jīng)有4人回復
-
276求調(diào)劑。有半年電池和半年高分子實習經(jīng)歷
已經(jīng)有9人回復
-
263求調(diào)劑
已經(jīng)有9人回復
-
求老師收我
已經(jīng)有3人回復
-
招08考數(shù)學
已經(jīng)有12人回復
-
接收2026碩士調(diào)劑(學碩+專碩)
已經(jīng)有4人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 存在這樣的函數(shù)嗎
已經(jīng)有11人回復
- 有沒有人知道塑性功怎么計算
已經(jīng)有10人回復
- 關于函數(shù)共軛的問題，求助
已經(jīng)有6人回復
- 急求解答，關于力場的電荷參數(shù)�。�！
已經(jīng)有4人回復
- 尋找擬合函數(shù)（由0增長，趨近于常值，并且0處的導數(shù)值為0等諸多條件）
已經(jīng)有7人回復
- 分子中電荷變化種種
已經(jīng)有82人回復
- 至少500金幣求解答一個問題-磷氧鍵（主要是P=O雙鍵，還有P-O單鍵）鍵能多少？
已經(jīng)有12人回復
- 請教密度泛函理論中屏蔽庫倫勢的問題
已經(jīng)有14人回復
- Gulp 擬合勢函數(shù) 初值設置的問題
已經(jīng)有12人回復
- 求概率密度函數(shù)的解析式
已經(jīng)有9人回復
- 【求助】第一性原理與密度泛函理論的關系是什么
已經(jīng)有15人回復
- 【求助】RDF分析-lammps
已經(jīng)有12人回復
- 【交流】歸一化結(jié)構(gòu)因子ESEL的上限和下限值怎么設置？
已經(jīng)有3人回復
|
本帖內(nèi)容被屏蔽 |
2樓2015-09-14 18:51:44
3樓2015-09-24 17:08:17
weihui333
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 134.1
- 帖子: 17
- 在線: 30.8小時
- 蟲號: 2731250
- 注冊: 2013-10-17
- 性別: GG
- 專業(yè): 化學反應工程
4樓2015-09-26 08:01:55
5樓2016-09-11 14:02:20
相關版塊跳轉(zhuǎn)
我要訂閱樓主 牛排小黑嬌 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調(diào)劑 +4 | niko- 2026-03-22 | 4/200 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +3 | jiajunX 2026-03-22 | 3/150 |
|
|
|
[考研] 環(huán)境學碩288求調(diào)劑 +6 | 皮皮皮123456 2026-03-22 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南大，0703化學，分數(shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿深大，0703化學，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研] 生物學調(diào)劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|