 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2335 | 回復(fù): 4 | |||
牛排小黑嬌新蟲 (初入文壇)
|
[求助]
勢(shì)函數(shù)解析問題
|
|
我想做Mo-Ag的合金擴(kuò)散模擬，但是并沒有現(xiàn)成的勢(shì)函數(shù)文件，我查閱文獻(xiàn)查到了有人擬合過這個(gè)合金的勢(shì)函數(shù)。 勢(shì)函數(shù)形式是FS勢(shì)格式和EAM勢(shì)差不多的，我看lammps手冊(cè)這個(gè)勢(shì)函數(shù)可以通過一個(gè)DYNAMO的程序解析出來，不知道您是不是也做過這個(gè)勢(shì)函數(shù)的解析文件，我想讓您幫我看一下我的程序是不 是有問題，我之前弄出來的模擬結(jié)果只能運(yùn)行一行。希望尚老師您能抽時(shí)間幫我看下，萬分感謝！ #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; 這邊取值的標(biāo)準(zhǔn)是？ 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 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 分子動(dòng)力學(xué)模擬 |
» 猜你喜歡
0854 考研調(diào)劑 招生了！AI 方向
已經(jīng)有14人回復(fù)
-
328求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有5人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有14人回復(fù)
-
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
0854電子信息求調(diào)劑 324
已經(jīng)有3人回復(fù)
-
303求調(diào)劑
已經(jīng)有3人回復(fù)
-
336求調(diào)劑
已經(jīng)有3人回復(fù)
-
306求調(diào)劑
已經(jīng)有9人回復(fù)
-
工科0856求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 存在這樣的函數(shù)嗎
已經(jīng)有11人回復(fù)
- 有沒有人知道塑性功怎么計(jì)算
已經(jīng)有10人回復(fù)
- 關(guān)于函數(shù)共軛的問題，求助
已經(jīng)有6人回復(fù)
- 急求解答，關(guān)于力場(chǎng)的電荷參數(shù)！��！
已經(jīng)有4人回復(fù)
- 尋找擬合函數(shù)（由0增長(zhǎng)，趨近于常值，并且0處的導(dǎo)數(shù)值為0等諸多條件）
已經(jīng)有7人回復(fù)
- 分子中電荷變化種種
已經(jīng)有82人回復(fù)
- 至少500金幣求解答一個(gè)問題-磷氧鍵（主要是P=O雙鍵，還有P-O單鍵）鍵能多少？
已經(jīng)有12人回復(fù)
- 請(qǐng)教密度泛函理論中屏蔽庫倫勢(shì)的問題
已經(jīng)有14人回復(fù)
- Gulp 擬合勢(shì)函數(shù) 初值設(shè)置的問題
已經(jīng)有12人回復(fù)
- 求概率密度函數(shù)的解析式
已經(jīng)有9人回復(fù)
- 【求助】第一性原理與密度泛函理論的關(guān)系是什么
已經(jīng)有15人回復(fù)
- 【求助】RDF分析-lammps
已經(jīng)有12人回復(fù)
- 【交流】歸一化結(jié)構(gòu)因子ESEL的上限和下限值怎么設(shè)置？
已經(jīng)有3人回復(fù)
|
本帖內(nèi)容被屏蔽 |
2樓2015-09-14 18:51:44
3樓2015-09-24 17:08:17
weihui333
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 134.1
- 帖子: 17
- 在線: 30.8小時(shí)
- 蟲號(hào): 2731250
- 注冊(cè): 2013-10-17
- 性別: GG
- 專業(yè): 化學(xué)反應(yīng)工程
4樓2015-09-26 08:01:55
涼yeah
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 34.4
- 帖子: 97
- 在線: 87.8小時(shí)
- 蟲號(hào): 4065768
- 注冊(cè): 2015-09-11
- 專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
5樓2016-09-11 14:02:20
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 牛排小黑嬌 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +8 | 材料學(xué)257求調(diào)劑 2026-03-23 | 9/450 |
|
|---|---|---|---|---|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 293求調(diào)劑 +12 | zjl的號(hào) 2026-03-16 | 17/850 |
|
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|