 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1211 | 回復(fù): 2 | |||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看 | |||
he66yuan88金蟲 (正式寫手)
|
[交流]
【求助】一個(gè)最最初級的問題：晶格中的原子坐標(biāo)（篇幅有點(diǎn)長，請見諒，請耐心） 已有2人參與
|
||
|
一個(gè)最最初級的問題：晶格中的原子坐標(biāo) 數(shù)據(jù)取自于LAMMPS中的例子：GayBerne ellipsoids in LJ background fluid 路徑是：lammps-17Feb11/examples/ellipse/log.ellipse.gayberne.10Sep10.linux.1 此例是一個(gè)GayBerne ellipsoids勢，橢圓體原子，npt系綜問題。 計(jì)算之后，log文件的前幾行為： units lj atom_style ellipsoid dimension 2 lattice sq 0.02 ###### Lattice spacing in x,y,z = 7.07107 7.07107 7.07107 ##### region box block 0 20 0 20 -0.5 0.5 ##### create_box 2 box Created orthogonal box = (0 0 -3.53553) to (141.421 141.421 3.53553) 1 by 1 by 1 processor grid create_atoms 1 box Created 400 atoms 其中用#####標(biāo)識的三行描述了晶格的構(gòu)造 lattice sq 0.02 ###### 給出了2d的正方形晶格，basis atom在正方形的左下角，約化的數(shù)密度rho*=0.02 region box block 0 20 0 20 -0.5 0.5 ##### 給出了晶格（box）的原點(diǎn)(0,0,0),及box的尺寸：20，20，1 Lattice spacing in x,y,z = 7.07107 7.07107 7.07107 #####是計(jì)算之后的晶格 那么，我的問題是：定義的和計(jì)算后得到的晶格似乎不合？我的理由是： 首先，對原子總數(shù)：假如一個(gè)20*20的正方形每0.02的約化距離一個(gè)原子的話，應(yīng)該是100*100=10000個(gè)原子，這樣與400不合 或者按照約畫數(shù)密度的定義：rho*=N/L^3=0.02，又因?yàn)閦方向厚度為1，所以求出的原子總數(shù)為8，也與400不合 如果，按照一個(gè)單位sigma一個(gè)原子的話，對于20*20的正方形的話，那就剛好是400個(gè)原子，但問題是這樣就與約化數(shù)密度rho*=0.02不合 其次，計(jì)算之后x,y,z = 7.07107 7.07107 7.07107，怎么看，我都看不出來，這個(gè)晶格矢量與b0x的尺寸20,20,1有什么關(guān)系 以上便是關(guān)于原子坐標(biāo)最最初級的問題。 還有一個(gè)問題也是關(guān)于原子坐標(biāo)的： 根據(jù)自己所學(xué)力學(xué)的觀點(diǎn)，只要構(gòu)型確定（即運(yùn)動學(xué)描述），加上力場確定（即動力學(xué)描述），一切相關(guān)的量即確定。 那么原子偏離原始晶格的尺寸以和外部邊界平衡，假如這種偏離造成晶胞變化，比如從正交晶系變?yōu)樾本�，那么這種偏離在npt/nph系綜 下為多大才是被允許的，怎樣描述這樣變化呢？這是否就是產(chǎn)生了相變？ 以上都是一些很基本的問題，由于自己非物理專業(yè)的，請高手們解答。謝謝 |
回復(fù)此樓» 猜你喜歡
299求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有6人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有8人回復(fù)
-
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有7人回復(fù)
-
289求調(diào)劑
已經(jīng)有13人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有9人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
調(diào)劑求院校招收
已經(jīng)有5人回復(fù)
-
調(diào)劑310
已經(jīng)有8人回復(fù)
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號: 659094
- 注冊: 2008-11-21
- 專業(yè): 金屬結(jié)構(gòu)材料
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
ghcacj(金幣+10, 模擬EPI+1): 精彩答疑，支持 2011-03-07 15:46:45
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
ghcacj(金幣+10, 模擬EPI+1): 精彩答疑，支持 2011-03-07 15:46:45
|
你大概沒有理解 lattice和region 這兩個(gè)命令的意思。 或者也可能沒有很好地理解約化密度的概念，或者沒有仔細(xì)閱讀說明書。 我們看看說明書，Lattice的解釋，有這么一段： Lattices of style sq and sq2 are 2d lattices that define a square unit cell with edge length = 1.0. This means a1 = 1 0 0 and a2 = 0 1 0. A sq lattice has 1 basis atom at the lower-left corner of the square. A sq2 lattice has 2 basis atoms, one at the corner and one at the center of the square. 。。。。 上一段是說， Lattice sq這個(gè)命令，生成二維的方形晶格元胞，其元胞邊長為1.0。對于sq 晶格，它的元胞里有一個(gè)原子，位于方形晶胞的左下角。 再往下看，有這一段，解釋后邊的scale這個(gè)數(shù)字： For unit style lj, the scale argument is the Lennard-Jones reduced density, typically written as rho*. LAMMPS converts this value into the multiplicative factor via the formula "factor^dim = rho/rho*", where rho = N/V with V = the volume of the lattice unit cell and N = the number of basis atoms in the unit cell (described below), and dim = 2 or 3 for the dimensionality of the simulation. Effectively, this means that if LJ particles of size sigma = 1.0 are used in the simulation, the lattice of particles will be at the desired reduced density. 看清紅字表示的公式。你現(xiàn)在是二維的晶格，dim=2，rho*=0.02是你輸入的。rho=N/V, N是元胞中的basis原子數(shù)，上面說了，是一個(gè)，N=1，V是元胞體積，上面也說了，生成的是邊長為1.0的方形元胞，所以有V=1.0，一算就知道 了，標(biāo)定因子factor=sqrt(50)=7.07107。 region box block 0 20 0 20 -0.5 0.5 這個(gè)意思是生成20*20的格子，你現(xiàn)在一個(gè)元胞里有一個(gè)原子，20*20的晶格里面多少原子？400個(gè)。 至于后邊關(guān)于相變的疑問，我還沒有能夠理解清楚你的問題。不大好回答。 [ Last edited by 老虎大王 on 2011-3-7 at 14:10 ] |
2樓2011-03-07 14:08:11
he66yuan88
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 549.8
- 散金: 20
- 紅花: 1
- 帖子: 458
- 在線: 100.4小時(shí)
- 蟲號: 1183270
- 注冊: 2011-01-04
- 性別: GG
- 專業(yè): 固體力學(xué)
3樓2011-03-07 14:41:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 he66yuan88 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +9 | 辰沐5211314 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +4 | Shern—- 2026-03-28 | 4/200 |
|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +10 | Mr. Z 2026-03-25 | 10/500 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|