 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 1907 | 回復: 6 | ||
zhangguangping木蟲 (著名寫手)
|
[求助]
【求助】DDscat軟件中關于獲得target的lab frame坐標系中坐標的問題
|
|
首先（這一塊跟問題沒什么關系，可以不看），非常感謝@donkeypku的資源帖(http://www.gaoyang168.com/bbs/viewthread.php?tid=3038851&fpage=2&target=blank)，介紹了DDA方法的軟件之一DDscat。我本是準備要研究納米粒子的光學性質(zhì)的。起初選擇了Meep，先不說Meep安裝起來多么費勁，光Meep的輸入文件我學習了2個周還是一頭霧水。Meep的輸出文件準備要借助于別的語言來寫。我的天啊，感覺被一層層的牽引到學習語言上去了。當然我不否Meep是一個上手慢，但是掌握之后比較好用的軟件。在論壇上搜索Meep的資源之后，在一個帖子(http://www.gaoyang168.com/bbs/viewthread.php?tid=3044900&fpage=1)中看到@donkeypku提到DDA,因此我開始將注意力轉(zhuǎn)移到DDA上。經(jīng)過3天的學習，基本掌握了DDscat的計算，并且還計算了一些自己想算的東西，結(jié)果還是比較滿意的。 現(xiàn)在得到的結(jié)果在ParaView中看都是TF坐標系下面，我想把計算的結(jié)果轉(zhuǎn)換到LF坐標系下面，得到數(shù)據(jù)與具體LF坐標(x,y,z)_LF的對應。 在DDscat中，有Lab Frame （LF）和Target Frame (TF)兩個坐標系，我現(xiàn)在對這兩個坐標系的表示有些疑問。 按照DDscat 7.3.0 User Guide p34 §18中提到的，在TF坐標系中，target的每一個dipole的坐標都是表示為用d（dipole的等效邊長，dipole的體積V=d*d*d）做單位，有點類似于分數(shù)坐標的意味（不過不同）。如果在shape.dat文件中，某一個dipole的坐標是(i,j,k)，那么它在TF坐標系中的物理坐標（真實坐標）按照 p34 §18中敘述的應該是(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的物理坐標，需要按照p35 §19.1中給出的TF坐標跟LF轉(zhuǎn)換關系來轉(zhuǎn)換。 但是按照p39 §21.1中對shape.dat的解釋，看到在對target中的每一個diople的描寫的時候，diople之間的間距可以不一定是d，并且是x,y,z三個方向的間距可以不同。如下： 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點陣列在x,y,z方向的間隔，還是用d做單位。關于這個地方的解釋，手冊上沒有多說。因此想搞清楚第三行的含義，去看代碼。reashp.f90中寫到DX(3)=dx/d, dy/d, dz/d [where d=(dx*dy*dz)^{1/3}]。因此，我認為這兒的DX(1), DX(2),DX(3)給出的是dipole在TF坐標下面x,y,z方向的間距(單位d)，這個也正如shape.dat文件描述的一樣：lattice spacings. 如果上面的理解是對的，那么 p34 §18中求某一個dipole在TF坐標系中的坐標的公式是錯誤的，應該是修改為：(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 也就是說既然dipole在xyz方向的間距不在是相同的d，而分別是dx/d,dy/d,dz/d那么應該計算dipole的坐標的時候應該考慮上這個不同的間距。 請高手指點！ @donkeypku @mahong123 @martinbible @bin2boy @tianwaiflyke @mp3gs @lkm8992 @79084787 @Ra冉 @jiangcong @zyy_323 @dpzxwswswws @rosefly88 @frank_zhan @xinyzy @思坦 @xiegangmai |
回復此樓» 猜你喜歡
化工學院（學院正式調(diào)劑群+聯(lián)系方式）
已經(jīng)有0人回復
-
化工學院（學院官方群）-院士團隊招生-學生自由度高-培養(yǎng)模式成熟高效
已經(jīng)有0人回復
-
無機化學論文潤色/翻譯怎么收費?
已經(jīng)有57人回復
-
黑龍江大學分析化學招調(diào)劑
已經(jīng)有0人回復
-
武漢高校國家級人才團隊招收材料，化學，紡織等專業(yè)碩士調(diào)劑生
已經(jīng)有26人回復
-
紡大化工學院（官方群）-院士團隊招生-學生自由度高發(fā)展好-培養(yǎng)模式成熟
已經(jīng)有0人回復
-
溫州大學化材學院王娟課題組招生
已經(jīng)有5人回復
-
學院官方群-院士團隊招生-學生自由度高-培養(yǎng)模式成熟已有10余名學生赴985海外深造
已經(jīng)有0人回復
-
武漢紡織大學化工學院官方群-院士團隊招生-學生自由度高-10余名學生赴985海外深造
已經(jīng)有0人回復
-
南昌航空大學蔣華麟教授課題組招收化學、環(huán)境、碳達峰碳中和及相關專業(yè)碩士信息
已經(jīng)有0人回復
-
湖北師范大學復試調(diào)劑
已經(jīng)有0人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 【求助】離散偶極近似（ddscat7.3）
已經(jīng)有12人回復
- 【求助】ddscat的有關內(nèi)容
已經(jīng)有2人回復
zhangguangping
木蟲 (著名寫手)
- 應助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
2樓2014-11-13 12:03:46
zhangguangping
木蟲 (著名寫手)
- 應助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
3樓2014-11-18 16:51:17
4樓2014-11-19 16:59:16
zhangguangping
木蟲 (著名寫手)
- 應助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2014-11-20 06:20:53
6樓2014-11-21 16:13:33
7樓2014-11-21 17:18:00
相關版塊跳轉(zhuǎn)
我要訂閱樓主 zhangguangping 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817化工學碩調(diào)劑 +6 | 努力上岸中！ 2026-03-31 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 266分，求材料相關專業(yè)調(diào)劑 +8 | 哇呼哼呼哼 2026-03-30 | 10/500 |
|
|
[考研] 一志愿南開大學0710生物學359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學 材料與化工（085600），已過六級。 +9 | XX小鄧 2026-03-29 | 9/450 |
|
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研] 環(huán)境科學與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6 | 來財aa 2026-03-25 | 6/300 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 本科新能源科學與工程，一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 考研一志愿蘇州大學初始315（英一）求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|