 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (2678)
- >考研 (327)
- >導(dǎo)師招生 (231)
- >蟲(chóng)友互識(shí) (103)
- >文獻(xiàn)求助 (47)
- >考博 (38)
- >論文投稿 (37)
- >碩博家園 (33)
- >休閑灌水 (23)
- >博后之家 (22)
- >基金申請(qǐng) (16)
- >找工作 (16)
- >文獻(xiàn)求助完結(jié)子版 (14)
- >教師之家 (12)
- >招聘信息布告欄 (11)
- >論文道賀祈福 (11)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 1906 | 回復(fù): 6 | ||
zhangguangping木蟲(chóng) (著名寫手)
|
[求助]
【求助】DDscat軟件中關(guān)于獲得target的lab frame坐標(biāo)系中坐標(biāo)的問(wèn)題
|
|
首先（這一塊跟問(wèn)題沒(méi)什么關(guān)系，可以不看），非常感謝@donkeypku的資源帖(http://www.gaoyang168.com/bbs/viewthread.php?tid=3038851&fpage=2&target=blank)，介紹了DDA方法的軟件之一DDscat。我本是準(zhǔn)備要研究納米粒子的光學(xué)性質(zhì)的。起初選擇了Meep，先不說(shuō)Meep安裝起來(lái)多么費(fèi)勁，光Meep的輸入文件我學(xué)習(xí)了2個(gè)周還是一頭霧水。Meep的輸出文件準(zhǔn)備要借助于別的語(yǔ)言來(lái)寫。我的天啊，感覺(jué)被一層層的牽引到學(xué)習(xí)語(yǔ)言上去了。當(dāng)然我不否Meep是一個(gè)上手慢，但是掌握之后比較好用的軟件。在論壇上搜索Meep的資源之后，在一個(gè)帖子(http://www.gaoyang168.com/bbs/viewthread.php?tid=3044900&fpage=1)中看到@donkeypku提到DDA,因此我開(kāi)始將注意力轉(zhuǎn)移到DDA上。經(jīng)過(guò)3天的學(xué)習(xí)，基本掌握了DDscat的計(jì)算，并且還計(jì)算了一些自己想算的東西，結(jié)果還是比較滿意的。 現(xiàn)在得到的結(jié)果在ParaView中看都是TF坐標(biāo)系下面，我想把計(jì)算的結(jié)果轉(zhuǎn)換到LF坐標(biāo)系下面，得到數(shù)據(jù)與具體LF坐標(biāo)(x,y,z)_LF的對(duì)應(yīng)。 在DDscat中，有Lab Frame （LF）和Target Frame (TF)兩個(gè)坐標(biāo)系，我現(xiàn)在對(duì)這兩個(gè)坐標(biāo)系的表示有些疑問(wèn)。 按照DDscat 7.3.0 User Guide p34 §18中提到的，在TF坐標(biāo)系中，target的每一個(gè)dipole的坐標(biāo)都是表示為用d（dipole的等效邊長(zhǎng)，dipole的體積V=d*d*d）做單位，有點(diǎn)類似于分?jǐn)?shù)坐標(biāo)的意味（不過(guò)不同）。如果在shape.dat文件中，某一個(gè)dipole的坐標(biāo)是(i,j,k)，那么它在TF坐標(biāo)系中的物理坐標(biāo)（真實(shí)坐標(biāo)）按照 p34 §18中敘述的應(yīng)該是(x,y,z)_TF=[(i,j,k)+(x0,y0,z0)]*d，其中，(x0,y0,z0)是(i,j,k)=(0,0,0)在TF中的位置(單位是d)。要想得到(i,j,k)在LF的物理坐標(biāo)，需要按照p35 §19.1中給出的TF坐標(biāo)跟LF轉(zhuǎn)換關(guān)系來(lái)轉(zhuǎn)換。 但是按照p39 §21.1中對(duì)shape.dat的解釋，看到在對(duì)target中的每一個(gè)diople的描寫的時(shí)候，diople之間的間距可以不一定是d，并且是x,y,z三個(gè)方向的間距可以不同。如下： 1.000000 0.000000 0.000000 = A_1 vector 0.000000 1.000000 0.000000 = A_2 vector 1.000000 1.000000 1.000000 = lattice spacings (d_x,d_y,d_z)/d 0.00000 0.00000 0.00000 = lattice offset x0(1-3) = (x_TF,y_TF,z_TF)/d for dipole 0 0 0 其中上面第三行給出dipole點(diǎn)陣列在x,y,z方向的間隔，還是用d做單位。關(guān)于這個(gè)地方的解釋，手冊(cè)上沒(méi)有多說(shuō)。因此想搞清楚第三行的含義，去看代碼。reashp.f90中寫到DX(3)=dx/d, dy/d, dz/d [where d=(dx*dy*dz)^{1/3}]。因此，我認(rèn)為這兒的DX(1), DX(2),DX(3)給出的是dipole在TF坐標(biāo)下面x,y,z方向的間距(單位d)，這個(gè)也正如shape.dat文件描述的一樣：lattice spacings. 如果上面的理解是對(duì)的，那么 p34 §18中求某一個(gè)dipole在TF坐標(biāo)系中的坐標(biāo)的公式是錯(cuò)誤的，應(yīng)該是修改為：(x,y,z)_TF=[(i,j,k)+(x0,y0,z0)]*d x_TF=(i*DX(1)+x0)*d y_TF=(j*DX(2)+y0)*d z_TF=(k*DX(3)+z0)*d 也就是說(shuō)既然dipole在xyz方向的間距不在是相同的d，而分別是dx/d,dy/d,dz/d那么應(yīng)該計(jì)算dipole的坐標(biāo)的時(shí)候應(yīng)該考慮上這個(gè)不同的間距。 請(qǐng)高手指點(diǎn)！ @donkeypku @mahong123 @martinbible @bin2boy @tianwaiflyke @mp3gs @lkm8992 @79084787 @Ra冉 @jiangcong @zyy_323 @dpzxwswswws @rosefly88 @frank_zhan @xinyzy @思坦 @xiegangmai |
回復(fù)此樓» 猜你喜歡
化工學(xué)院（學(xué)院正式調(diào)劑群+聯(lián)系方式）
已經(jīng)有0人回復(fù)
-
化工學(xué)院（學(xué)院官方群）-院士團(tuán)隊(duì)招生-學(xué)生自由度高-培養(yǎng)模式成熟高效
已經(jīng)有0人回復(fù)
-
金屬材料論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有279人回復(fù)
-
黑龍江大學(xué)分析化學(xué)招調(diào)劑
已經(jīng)有0人回復(fù)
-
武漢高校國(guó)家級(jí)人才團(tuán)隊(duì)招收材料，化學(xué)，紡織等專業(yè)碩士調(diào)劑生
已經(jīng)有26人回復(fù)
-
紡大化工學(xué)院（官方群）-院士團(tuán)隊(duì)招生-學(xué)生自由度高發(fā)展好-培養(yǎng)模式成熟
已經(jīng)有0人回復(fù)
-
溫州大學(xué)化材學(xué)院王娟課題組招生
已經(jīng)有5人回復(fù)
-
學(xué)院官方群-院士團(tuán)隊(duì)招生-學(xué)生自由度高-培養(yǎng)模式成熟已有10余名學(xué)生赴985海外深造
已經(jīng)有0人回復(fù)
-
武漢紡織大學(xué)化工學(xué)院官方群-院士團(tuán)隊(duì)招生-學(xué)生自由度高-10余名學(xué)生赴985海外深造
已經(jīng)有0人回復(fù)
-
南昌航空大學(xué)蔣華麟教授課題組招收化學(xué)、環(huán)境、碳達(dá)峰碳中和及相關(guān)專業(yè)碩士信息
已經(jīng)有0人回復(fù)
-
湖北師范大學(xué)復(fù)試調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】離散偶極近似（ddscat7.3）
已經(jīng)有12人回復(fù)
- 【求助】ddscat的有關(guān)內(nèi)容
已經(jīng)有2人回復(fù)
zhangguangping
木蟲(chóng) (著名寫手)
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲(chóng)號(hào): 529624
- 注冊(cè): 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
2樓2014-11-13 12:03:46
zhangguangping
木蟲(chóng) (著名寫手)
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲(chóng)號(hào): 529624
- 注冊(cè): 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
3樓2014-11-18 16:51:17
_moving_
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 13
- 在線: 3.8小時(shí)
- 蟲(chóng)號(hào): 2761144
- 注冊(cè): 2013-10-29
- 專業(yè): 光學(xué)
4樓2014-11-19 16:59:16
zhangguangping
木蟲(chóng) (著名寫手)
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲(chóng)號(hào): 529624
- 注冊(cè): 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2014-11-20 06:20:53
_moving_
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 13
- 在線: 3.8小時(shí)
- 蟲(chóng)號(hào): 2761144
- 注冊(cè): 2013-10-29
- 專業(yè): 光學(xué)
6樓2014-11-21 16:13:33
_moving_
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 13
- 在線: 3.8小時(shí)
- 蟲(chóng)號(hào): 2761144
- 注冊(cè): 2013-10-29
- 專業(yè): 光學(xué)
7樓2014-11-21 17:18:00
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhangguangping 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +11 | minht 2026-03-29 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +7 | 披星河 2026-03-30 | 7/350 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +15 | kikiki7 2026-03-30 | 15/750 |
|
|
|
[考研] 一志愿中海洋材料357 +3 | 麥恩莉. 2026-03-30 | 3/150 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +5 | 朵朵- 2026-03-26 | 7/350 |
|
|
|
[考研] 309求調(diào)劑 +15 | 誰(shuí)不是少年 2026-03-29 | 15/750 |
|
|
|
[考研] 085601材料工程找調(diào)劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
|
[考研] 283求調(diào)劑（080500） +14 | A child 2026-03-27 | 14/700 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|