- 追加された行はこの色です。
- 削除された行はこの色です。
前田 大聖（まえだ たいせい）~
仙台高専~
2018年8月6日〜%%9月14日%%11月2日
#contents
*困っていること [#bfd36b9b]
*日報 [#h5fade90]
**11月2日 [#z572e777]
やったこと
-直方体の厚さを140[nm]として，直方体の幅と長さの変化させたシミュレーションの解析
--それぞれの直方体の長さの変化を30[nm]，幅の変化を30[nm]としたものと，
長さの変化を30[nm]，幅の変化なしのものと，
長さの変化を30[nm]，幅の変化を10[nm]とした3つのパターンで比較を行った
--長さの変化が30[nm]，幅の変化が30[nm]のものが最も円偏光の偏りが大きかった
---ただし，楕円率のヒストグラムの山が高くなるだけでピークは楕円率が0.5程度の位置であった
-インターンシップのまとめの作成
**11月1日 [#h28c54f2]
やったこと
-直方体の厚さを140[nm]とし，3つの直方体の長辺を270[nm],300[nm],330[nm]としたときのシミュレーションの解析
--3つの直方体の長辺が300[nm]のものに比べて若干円偏光の偏りが増えたように見える
--それぞれの直方体の電荷分布をプロットしたところ，波源に一番近い直方体とほかの直方体とで電荷が回る向きが変わっていた
--それぞれの直方体の表面のエネルギーをプロットしたところ，
波源に一番近い直方体の表面よりほかの直方体の表面の方がエネルギーが高かった
明日やること
-長辺の長さをさらに変えた状態でシミュレーションの解析
**10月31日 [#n5d5d960]
やったこと
-直方体を厚くとった状態のシミュレーションの解析
--直方体の厚さを60[nm],100[nm],140[nm]とで比較すると100[nm]と140[nm]は60[nm]より
明らかに円偏光の偏りが増えていた
--100[nm]と140[nm]では140[nm]の方がわずかに円偏光の偏りが増えていた
-明日やること
--直方体の厚さを140[nm]にとり3つの直方体を長辺を少しずつ変化させたシミュレーションの解析
**10月30日 [#l46fa3e6]
やったこと
-間隔を保ったまま直方体を薄くした状態の解析
--厚さを薄くする前のものほど円偏光の偏りがなかった
---中心距離の間隔や直方体同士の間隔より厚さを薄くしたことで円偏光の偏りが減った
--光源と反対方向の直方体から離れた面では直線偏光しかでていなかった
明日やること
-直方体を厚くした状態のシミュレーションの解析
**10月25日 [#sb1420f0]
やったこと
-直方体の厚さを薄くした八木アンテナをもした形状のシミュレーションの解析
--直方体を薄くする前と薄くした後で比較すると薄くする前の方が円偏光が偏っていた
---直方体を薄くしたが直方体どうしの中心距離が変わっていないため，
直方体同士の間隔が厚さを薄くする前に比べて広がったからかもしれない
-直方体を薄くして，直方体同士の間隔を薄くする前と同じ間隔でシミュレーション
来週やること
-間隔を保ったまま直方体を薄くした状態の解析
**10月24日 [#y138a267]
やったこと
-八木アンテナを模した形状でシミュレーションの解析
--エネルギー×楕円率を時間経過でプロットとともにヒストグラムも表示させるプログラムの開発
---今後このプログラムを用いてシミュレーションの評価をしていきたい
-直方体の厚さを今までより薄くしてたシミュレーションの実行
明日やること
-直方体の厚さを薄くしたシミュレーションの解析
**10月23日 [#o28def57]
やったこと
-八木アンテナを模した形状でシミュレーションの解析
--meepから直接エネルギーを出力させそのデータをもとに反射波のみのエネルギーをプロットしようと試みた
---今まではmeepが出力するx方向やy方向の電場から計算をしてエネルギーを算出していたがプロットに時間がかかりすぎていたので，meepから直接エネルギーを出力させることでプロットするプログラムの処理を減らすため
--反射波のみのエネルギーを取り出すのは難しかったが，金の直方体があるときとないときのエネルギーの差をプロットすること反射が及ぼすエネルギーをプロットすることができた
--丸岡さんのおかげで電場からエネルギーを算出しながらプロットするプログラムも5分程度実行ができるようになった
---フレームが切り替わっても前のプロットが消去できていなかったため，フレームが進むごとにプロットがたまって処理が重くなっていた
明日やること
-楕円率×エネルギーのプロットとヒストグラムの作成
**10月22日 [#p24a97fa]
やったこと
-八木アンテナを模した形状でシミュレーション
--3つ並べる直方体の長さを1:3の比を保ちつつ長さを変えた
--電場ベクトルには3つとも同じ直方体を使ったものとの差は見えなかった
--電荷分布は画像からは違いがわからなかったが誘起される電荷の最大値が
同じものを3つ並べたものより大きくなっていた
--楕円率は画像から大きな差はあまり見られなかったが
楕円率を横軸にとったヒストグラムを見るとヒストグラムのピークや楕円率の平均値が違っていた
明日やること
-楕円率×エネルギーのヒストグラムの作成
-解析から八木アンテナの変化させるパラメーターの模索
**10月19日 [#qf6fc162]
やったこと
-昨日の八木アンテナを模した形状でシミュレーションの解析の続き
--電場ベクトルのプロット
---直方体の周りにできる渦が表面では直方体の傾けた向きに伝播するのに対して
表面から離した面では同心円状に伝播していた
---電場ベクトルをベクトルのプロットをしたことで昨日の楕円率を測った時間の区間はまだ定常状態になっていない可能性がでた
--楕円率のプロットのし直し
---昨日のプロットした時間の区間より10周期後で楕円率をし直した
---昨日の楕円率のマップとは違う画像であったが，特徴は同じであった
--電荷密度のプロット
---電荷は直方体が1つのものと同じ向きで回転していた
---誘起される電荷量が増加していた
--エネルギーのプロット
---エネルギーの最大値に合わせてプロットしたため，直方体の端っこの強い部分しか色として現れなかった
---エネルギーのプロットのプログラムは7時間くらいかかっているので改善が必要
来週やること
-八木アンテナを模した形状のシミュレーション
--直方体は1対3の比率を保ったままパラメーターの変更
-楕円率を数値としてどれくらいの幅を持っているか確かめる
-エネルギーのプロットを短い時間で終わるように改善する
**10月18日 [#jc81ddf1]
やったこと
-以下のような八木アンテナを模した形状でシミュレーションの解析
--直方体を1つだけのときほど円偏光の偏りは見られなかったが
表面と中心から200[nm]離した時の偏りの差が大きく開いていた
---この差を維持したまま全体を左円偏光に偏らせたい
--表面と表面から離した面での楕円率のプロット
---表面から離れるほど面全体が左円偏光
http://flex.phys.tohoku.ac.jp/~maeda/fig/yagi-Rec-eps-xy-plane.png
http://flex.phys.tohoku.ac.jp/~maeda/fig/yagi-Rec-eps-xz-plane.png~
(濃い黒い部分が金を配置した位置)~
(左の図・・・横軸：x軸，縦軸：y軸)~
(右の図・・・横軸；x軸，縦軸；z軸)~
明日やること
-今日のシミュレーションの解析の続き
--表面や表面から離した面の電場ベクトルのプロット
--電荷分布のプロット
--エネルギーのプロット
**10月17日 [#gcbac14a]
やったこと
-以下の形状のシミュレーションの解析
--直方体を1つだけ配置した時より左円偏光の偏りが減っていた
--楕円率は中心から400[nm]離したxy平面では左円偏光と直線偏光しか出ていなかった
---エネルギーのプロットが済んでないので弱い左円偏光が出ている可能性もある
http://flex.phys.tohoku.ac.jp/~maeda/fig/3Rec-eps.png~
(横軸:x軸，縦軸:y軸，濃い黒い部分が金を配置した位置)~
(金の直方体とx軸のなす角60°)
明日やること
-大きさを変えずに直方体をz軸方向に並べてシミュレーション
-エネルギーのプロット
**10月16日 [#rd910c12]
やったこと
-円偏光が強く現れそうな形のアイデア出し
-出したアイデアからシミュレーションの実行
--解析は明日やる
明日やること
-今日行ったシミュレーションの解析
**10月15日 [#vb1fc1e6]
やったこと
-波源をPML中に含んだ状態で直方体の表面と表面から離した面の解析
--10月11日のシミュレーションの設定を変えたことで表面の振る舞いが変化した原因は
波源の大きさを変えたことに来るものだと分かった
--表面から離れるほどに円偏光の割合が多くなることは変わらなかった
--PMLの厚さを変えたことによる円偏光の割合はほとんど変わらなかった
--600[nm]の波長の入射光を用いて方位角が30°付近で最も強いピークが現れ，
直方体の中心から400[nm]離した面では60%の割合で左円偏光が現れていることが確認できた
明日やること
-左円偏光だけもしくは右円偏光だけを強く発生させる構造の模索
**10月12日 [#ha5254e8]
やったこと
-直方体を配置しないシミュレーションの実行と解析
--直線偏光を入れてもわずかに円偏光の割合があった（0.3%程度）
--波源から離れるほど円偏光の割合が増加した
---PMLで吸収できなかった反射からの足し合わせによるもの？
-光源をPML中に入れるのと入れない場合でのシミュレーションの実行と解析
--PML中に入れない場合では斜めの干渉縞ができていた
来週やること
-表面から離した平面でのエネルギーのプロット
**10月11日 [#u794c37f]
やったこと
-直方体の表面から離した平面での楕円率×エネルギーをプロット
--シミュレーション設定の昨日をプロットしたデータと変えた
（PMLの厚さと光源の大きさを変更）ことにより
表面での振る舞いが変わった
--表面から離れるほど円偏光の割合が多くなった
明日やること
-何も置かない状態でPMLからの反射波がどの程度あるのかを確認
-直方体を置かないシミュレーションでの楕円率が0になるのかを確認
-直方体を置かないシミュレーションでの楕円率×エネルギーのプロット
-表面から離した平面でのエネルギーのプロット
**10月10日 [#oddbecd1]
やったこと
-直方体の表面の楕円率×エネルギーをプロット
--S字のときと同様に入射光（直線偏光）の方位角依存性があった
--どの方位角でもsin関数のような軌跡を描いており，0°，90°，180°では
楕円率×エネルギーが0になった
--500[nm]の波長で45°の直線偏光の時が最も左円偏光を発し，
135°の直線偏光の時が最も右円偏光を発していることが分かった
-直方体の表面から離したxy平面で観測するシミュレーションの実行中
-xz平面でのエネルギーのプロット
--作図のプログラムが時間がかかるので明日結果を確かめる
明日やること
-直方体の表面から離した平面での楕円率×エネルギーのプロット
-xz平面でのエネルギーのプロットの確認とプログラム・図の評価
**10月9日 [#xafacdb6]
やったこと
-S字と逆S字について楕円率×エネルギーをプロット
--入射光（直線偏光）の方位角依存性があった
--S字と逆Sのグラフは点対称であった
-長方形についても楕円率×エネルギーをプロット
--間違った計算をしているシミュレーションがあったため，シミュレーションのやり直し中
明日やること
-長方形について楕円率×エネルギーをプロット
-xy平面だけでなく，xz平面やyz平面で切り取り，円偏光がz軸方向に伝播するか確かめる
**10月5日 [#s4ad14e3]
やったこと
-反射波と入射波の位相差を求めるプログラム作成
--長軸同士の時間のインデックス差から位相差をプロットするが，
それぞれ二つの長軸ベクトルがあるため，4通りの組み合わせあるが，
それらをまったく考慮せずにプログラムを組んでいたため，
位相差が周期ごとに変化していた
-楕円率×エネルギーから左円偏光と右円偏光がどれくらい出ているか調べた
--入射波の直線偏光の方位角の依存性があった
来週やること
-長方形に対しても円偏光の差を調べる
-反射波と入射波の位相差を求めるプログラム完成
**10月4日 [#wa9f8561]
やったこと
-反射波と入射波の位相差を求めるプログラム作成
--永遠に実行するというバグはなくなった
-S字のシミュレーションで左円偏光と右円偏光がどれくらい出ているかを調べた
--左円偏光，右円偏光が出ている場所をカウントしていたが，
電場強度も考慮するべきだった
明日やること
-強度×楕円率から，左円偏光と右円偏光がどれくらい出ているか調べる
-反射波と入射波の位相差を求めるプログラムの作成した図に対して考察
-全部のベクトルを足したときにそのベクトルの和が入射した直線偏光になるのか調べる
**10月3日 [#tcbac060]
やったこと
-反射波と入射波の位相差を求めるプログラムの作成
--完成していない
-昨日の発表の振り返り
明日やること
-反射波と入射波の位相差を求めるプログラムの完成
-S字の左円偏光と右円偏光が同じ分だけ出ているか確かめる
**10月1日 [#ee270619]
やったこと
-明日の発表に向けての発表練習
--スライドの作りに注意するポイントや発表する際に気を付けることなどアドバイスをもらった
---アウトラインは詳しく書く
---図に対して過不足なく説明の言葉を入れる
---結果に対しては口頭の説明だけでなく，スライドにもメッセージを残す
明日までにやること
-指摘された箇所を直しと発表練習
**9月28日 [#y51eb06b]
やったこと
-光の吸収率のプログラムミス発見及び修正
--吸収率は"1-(透過率+反射率)"で計算し，反射率の計算が常に0になっていたため
吸収率と思って求めていたのは"1-透過率"であった
-正方形に対して300[nm]から800[nm]まで変化させたときの吸収率波長依存性をプロット
--500[nm]のところだけ4%の吸収率が得られ，ほかの波長では吸収率が負になった
-S字構造に対して電荷分布のプロット
--S字の様々な場所で+の電荷と-の電荷が回転しているように見えた
来週やること
-印加した電場と金に発生する電場の位相差を楕円率をマッピングしたような形でプロット
-波長変化による電荷分布を見る
**9月27日 [#xe0905f6]
やったこと
-正方形に対して800[nm]の光を当ててシミュレーション
--光の吸収率は約77%であった
--昨日の600[nm]のシミュレーションではpml層を300[nm]分とっていたが、
今日のシミュレーションでは800[nm]を用いたので半波長の400[nm]分とった
---波長による比較を行いたいので波長以外の条件を同じにするべきだった
-現在、300[nm]から700[nm]まで波長を100[nm]ごとに変化させてシミュレーション中
--PMLはすべて400[nm]分とった
明日やること
-波長を300[nm]から700[nm]まで変化させたときのシミュレーションの結果の考察
**9月26日 [#h293ec4f]
やったこと
-正方形に対してのシミュレーション
--吸収がとれるように正方形の構造の前後に通った光のエネルギーためておく領域(flux)をセットしてシミュレーションをした
--300[nm]の正方形に対して600[nm]を中心とするガウス波源を入れて計算したところ、
約83%吸収されていた
明日やること
-fluxを離してシミュレーション
-波長を変えて吸収の違いを見る
-より正確な吸収率のプロット
--吸収率はジュール熱に変換されるので電場の2乗を時間経過とともに足し合わせることで積分したことになり各点での吸収率をプロットできる
**9月25日 [#nf84f163]
やったこと
-60[nm]、6000[nm]の波長を用いたシミュレーション
--60[nm]の波長では長方形の構造がダイポールアンテナのような振る舞いをしていることが確認された
-系のメッシュサイズを大きくても使う波長の10分の1程度にすることを学んだ
明日やること
-電場の吸収強度の図を作成できるようなシミュレーションを行い、吸収強度をプロットする
**9月21日 [#k5125145]
やったこと
-波長を変えて長方形状のシミュレーション
--昨日は600[nm]の波長を用いたので、500[nm]と700[nm]でシミュレーションを行った
--波長を長くすると全体的に現れる渦のようなものが大きくなった
来週やること
-波長を10倍や10分の1など極端に変えた状態でのシミュレーション
-正方形でのシミュレーション
**9月20日 [#yc7f9e0e]
やったこと
-長方形状シミュレーション
--pmlを厚くとるような系で再度シミュレーションを行った
--電場ベクトルと楕円率のプロットを行った
--長方形状の外側に円偏光が現れているのが確認された
明日やること
-波長を変えて長方形状のシミュレーションと解析
**9月19日 [#vd00b41b]
やったこと
-tube61のメモリ確認
--メモリの組み合わせを変えてtube61が起動するか確かめた
---様々なパターンで確認した結果どのメモリも認識され、すべて差した状態でもtube61が起動できた
---前回起動できなかったのはメモリがきちんと差さっていなかったためだと思われる
-長方形状のシミュレーション
--上記のtube61の作業に時間がかかりまだ作図はしていない
--pmlの厚さを入射光の半波長程度にするとよいとのことだったので、そちらのシミュレーションしたい
明日やること
-今日のシミュレーションの作図
-pmlの厚さを半波長程度にした場合のシミュレーション
**9月18日 [#ad01d96e]
やったこと
-S字の幅を太くしたものの楕円率のプロット
--S字の太さを変えても楕円率に大きな違いはなかった
-FLEXに接続
--WindowsUpdateをしたらVcXsrvでも接続できるようになった
-細くしたS字のシミュレーション
--作図はまだ
明日やること
-細くしたS字の結果の作図
-下の図のような2つの棒状のものでのそれぞれシミュレーション~
http://flex.phys.tohoku.ac.jp/~maeda/fig/Picture1.png
**9月14日 [#q9f28541]
やったこと
-S字の幅を太くしたものでエネルギーの時間経過に対するプロット
--太くする前のものと比べてエネルギーが小さくなっていた
-FLEXに接続
--VcXsrvからではなくWindows上のUbuntuからsshにより接続
--WindowsUpdateが来ていたので更新後にVcXsrvから再度ログインできるか確かめる
来週やること
-太くしたS字で楕円率のマッピング
-細くしたS字でシミュレーション
**9月13日 [#cf1a98c8]
やったこと
-学会の見学と手伝い
-修正したプログラムのソースコードのチェック
--VcXsrvからFLEXへ接続できなったため実際に動かしてチェックができなかった
(エクスプローラーなどではFLEX上のファイルを閲覧、コピーなどは可能な状態）
---設定ファイルを再度作り直すも接続ができない
明日やること
-再びVcXsrvでFLEXに接続できるようにする
-エネルギーの時間経過をプロットするプログラムの完成
**9月12日 [#sbafb24e]
やったこと
-昨日のエネルギープロットのプログラムの修正
--同じようなバグで困っている人の[[記事:https://stackoverflow.com/questions/39472017/how-to-animate-the-colorbar-in-matplotlib]]を発見し，問題解決中
-学会の見学と手伝い
--MEEPでシミュレーションしている古田先生のポスターを見学しに行った
明日やること
-エネルギーをプロットするプログラムの改善
**9月11日 [#j100fb0a]
やったこと
-S字の幅を太くした系でシミュレーション
--電場を時間経過出力する際にバグが発生したのでそれを解決中
---カラーバーが時間経過毎に出力され途中からカラーだけの画像になってしまっている
-学会の見学と手伝い
明日やること
-上記のバグ修正
-シミュレーションの考察
**9月10日 [#ef81d0a1]
やったこと
-S字と逆S字の直上での電場ベクトル、楕円率、エネルギーの時間平均をプロットした
--どちらの場合でもS字以外の部分で円偏光が見られた
--エネルギーは時間平均をプロットするより、時間経過をプロットした方がよいと
アドバイスをもらった
-電場ベクトルのアニメーションを今まではpng画像をプログラムで作り、Photoscapeというソフトウェアでアニメーションgif化していたのを、プログラムから直接アニメーションgifを作成することができた
-アニメーションgifから動画への変換もPythonで実現できた
--[[gifから動画の作成方法:http://flex.phys.tohoku.ac.jp/japanese/pukiwiki/index.php?matplotlib#k4236f9e]]
明日やること
-エネルギーの時間経過のプロット
-S字の形を変化させた場合のシミュレーション（直上とメッシュ1つ分だけ上）
--Sの幅を変える
**9月7日 [#v2eb92bf]
やったこと
-昨日のエラーの改善
--本来であればS字の構造があるものとないものと２つでシミュレーションを行い、
2つのシミュレーションの差から電場ベクトル、楕円率、電場強度を求めるのだが、
2つのシミュレーションが全く同じシミュレーションであったため差をとって0になったり、0除算をしていた
--そのためシミュレーションのやり直しをした
-丸岡さんとFDTD法の本を使って輪読をした
--FDTD法が扱うマクスウェル方程式では磁流を定義したマクスウェル方程式を使っていることを学んだ
-FNTG55の準備の手伝い
来週やること
-今日再び行ったシミュレーションの解析
-丸岡さんとの輪読の準備
-様々な状態でのシミュレーションと解析
**9月6日 [#g3f9d565]
やったこと
-電場ベクトル、楕円率、電場強度をプロットするプログラムはすべて別々だったので
1つのプログラムに統合
-S字構造の直上を見るシミュレーション(照射する直線偏光の角度を変化をさせつつ)
--プロットするプログラムの実行中にエラーが発生し、画像の作成ができなった
---以下のエラーメッセージが楕円率の算出する部分で表示された
RuntimeWarning: invalid value encountered in double_scalars
明日やること
-楕円率を算出する部分で想定外の値が入らないかのチェックをする
-昨日の「明日やること」のシミュレーションをして置き、データを解析する準備をしておく
**9月5日 [#j391b9a0]
やったこと
-小さい系で照射する直線偏光の角度を変化ながらシミュレーション
--S字構造の物質をグラフェンから金に変更した
--S字構造を領域の中心付近に移動
--S字構造の大きさが全体で700nmとし、照射する光の波長を400nmとした
--PMLを厚くとった
-シミュレーションの結果
--グラフェンのときは違った結果が観測された
--楕円率が0.5以下しか現れず、円偏光が発生していなかった
---グラフェンでのシミュレーションを行ったときに比べて、
S字構造が波長に対して相対的大きくなったからかもしれない
--電場の強度をプロットした
---S字の内部より、輪っかのようになっている部分の方が強度が大きかった
明日やること
-波長を変えてシミュレーションしてみる
-S字の直上の電場ベクトルを観察して回転しているかを見る
-楕円率のカラーバーをどんな画像に対しても-1から1の範囲でプロットする
-逆S字でのも同様のシミュレーションをしてみる
**9月4日 [#bbe0a4cd]
やったこと
-シミュレーションの時間を短くするために系を小さくした状態でも同様の結果が得られるかの検証
--同様の結果が得られなかった
---S字構造をPMLから離すべきだった
---S字構造はグラフェンである必要がなかった
---S字構造に比べて波長が短いもの光を照射するべきだった
---resolution(メッシュの細かさ)を下げたためS字構造からメッシュ一つ分離れた場所が遠くなった影響が出たのかもしれない
明日やること
-今日の反省を生かして小さい系でのシミュレーションをし直す
--S字構造をPMLから離す
--S字構造の物質を金や銀に変える
--S字構造の大きさより小さい波長の光（可視光領域）を使う
**9月3日 [#o20c26bf]
やったこと
-左円偏光と右円偏光を分ける2つのプログラムの差の考察
--y成分が正の長軸ベクトルと短軸ベクトルの時間のインデックスの差が
すべて1/4周期でなく、3/4周期のペアもあり、それを考慮せずにプログラムを組んでいたためだった
--3/4周期のペアも考慮することで、2つのプログラムが出力する画像が一致した
-0°から90°まで9°刻みで直線偏光を回転させたときのシミュレーション結果を図として作成した
--直線偏光が回転するとそれに伴い左円偏光と右円偏光が発生する位置が動いていることが確認できた
明日やること
-入手する直線偏向を90°から360°までの回転させてシミュレーション
-電場強度をマッピングするプログラムの開発
**8月31日 [#e3cbb761]
やったこと
-直線偏向を0°から90°まで9°刻みで回転させた場合のS字のシミュレーション
--シミュレーションの実行ファイルを同時にTubeに投げるとTubeが落ちてしまうため
プログラムを2つに分け、Tube60とTube61に実行させた
--実行途中でTube61が落ちたりしたためシミュレーションの結果を図として作成することはできなかった
-左円偏光と右円偏光を分ける2つのプログラムの差の考察
--一昨日書いたプログラムにはミスがあったのでその部分の修正
---修正しても2つのプログラムが出力する画像は同じにならなかった
来週やること
-直線偏向を0°から90°まで回転させたS字のシミュレーションの結果の考察
-左円偏光と右円偏光を分ける2つのプログラムが出力する画像の差がソースコードのどこから由来するものかを調べる
--2つのプログラムが判定している短軸と長軸のベクトルの絶対値をとる
--楕円を描いているか知るために、電場ベクトルの軌跡のマッピングする
--短軸と長軸長軸ベクトルの時間に対するインデックスの差が1/4周期になっているかの確認
--軸を45°傾けても同じ結果になるか確かめてみる
**8月30日 [#z9cdaf0c]
やったこと
-y軸反転させたS字構造の左円偏光と右円偏光の観察
--入射波をx軸に振動するパターンとy軸に振動するパターンの2つ行った
-短軸ベクトルと1/4周期前のベクトルで外積をとることで左円偏光と右円偏光の区別をつけるプログラムの作成
--このプログラムで作成された画像と昨日作成した画像が違う画像だったので画像の出力どのプログラムが正しいのか検証する必要がある
明日やること
-S字構造に角度を変えた直線偏向を当てた時の円偏光の観察
**8月29日 [#g79d7301]
やったこと
-楕円率と外積を用いて左円偏光と右円偏光の区別をするプログラムの作成
--プログラムをややこしく作ったため何をするプログラムかわかりにくいものとなった
--短軸ベクトルと1/4周期前のベクトルで外積をとることで左円偏光と右円偏光の区別をつけるようにする
明日やること
-y軸反転させたS字構造の左円偏光と右円偏光の観察
-左円偏光と右円偏光の区別するプログラムの簡易化(1/4周期前のベクトルを用いる)
**8月28日 [#u7bb7891]
やったこと
-楕円率を求めるプログラムの作成
-楕円率に応じて直線偏光と円偏光を区別をつけることができた
-楕円率のマッピングに誘電率のマッピングを重ねることができた
明日やること
-楕円率を求める際の長軸ベクトルと短軸のベクトルの外積をとることで
左円偏光と右円偏光の区別をするプログラムの作成
**8月27日 [#wa3ae2d9]
やったこと
-MEEPのプログラムの(output-epsilon)が出力するデータが複素誘電率なのか複素誘電率の実部だけなのかを検証した
--複素誘電率の実部だけを出力していた
-S字から出る円偏光の右円偏光だけを抽出するプログラムの作成(不完全)
--直線偏光と思われる部分まで右円偏光として抽出されていた
---直線偏光と思われる部分でも完全な直線偏光ではないとため、少しでも偏りが生じると円偏光として扱われるようなプログラムを組んでいたため
---楕円率を求め、楕円率から円偏光と直線偏光を分ける
明日やること
-楕円率を求めるプログラムの作成し、楕円率に応じて直線偏光と円偏光を分ける
--直線偏光と円偏光に分けることができてから外積を用いて右円偏光と左円偏光に分ける
**8月24日 [#a13a1173]
やったこと
-誘電率をマッピングするプログラムの作成(不完全)
--下の図のようなS字ではなく不完全なS字が出力される~
http://flex.phys.tohoku.ac.jp/~maeda/fig/not-perfect-S-eps.png
--複素誘電率を引っ張ってくる必要がある
来週やること
-S字のシミュレーションが岡本先生の岡本先生の論文と同じような結果になること確かめる
--電場ベクトルが回転する向き(時計回りか反時計回りか)でプロットするベクトルの色を変える
-S以外のアルファベットなどに対してのシミュレーションを行い，円偏光が観測されるかどうか調べる
**8月23日 [#kecc9ac6]
やったこと
-S字構造のグラフェンからの反射波の電場マッピングするプログラムの作成
--右円偏光と左円偏光がそれぞれ確認できた
-S字構造のグラフェンの観察面をグラフェンさらに離した状態のシミュレーションの実行
--図の作成はまだ至っていない
明日やること
-電場のマップではどこにS字構造があるか正確に把握できないため
電場のマップに誘電率のマップを差し込むプログラムの作成
-今日実行したシミュレーションの結果の図の作成
**8月22日 [#i9b3ddca]
やったこと
-S字構造のグラフェンからメッシュ1つ分だけ離れた位置での電場の観察シミュレーション
--グラフェンからの反射波のみを観察したいのでグラフェンがある状態とない状態でシミュレーション
---グラフェンがある状態でx方向とy方向にそれぞれ正に強く出る部分と負に強く出る部分が観測できた
---グラフェンがある状態とない状態で差をとり、グラフェンからの反射波のみを取り出すことができた
明日やること
-以下のどちらかをやる予定
--観測地点をグラフェンからさらに離してみる
--S字構造に円偏光の光を入れてみる
困っていること
-PukiWikiに画像をアップロードしようしたが、管理者を知らないためアップロードできなかった
**8月21日 [#ic418e87]
やったこと
-S字構造のグラフェンと光源の間での電場の観測シミュレーション
--観測は各成分電場をカラーマップに落とし込んで観測したのだが、
端の部分の電場が大きすぎて過ぎて実際に見たかったグラフェンがある場所で電場状況を見ることができなかった
---カラーマップに落とす前に全体の電場の値を見るべきだった
---端の部分を無視して再度電場の観察をする
明日やること
-S字構造のグラフェンからメッシュ1つ分だけ離れた位置で電場の観察シミュレーション
**8月20日 [#wd007b5e]
やったこと
-誘電体のマルチレイヤーのシミュレーションの実行
--グラフェンの光吸収率が50%になった
-条件を変えた状態でのシミュレーションの実行
--t61の電源が落ちていてローカルのPCを使ったため何回もシミュレーションする時間がなく、プログラムのミスに気付くのに時間がかかった
--現在、出力結果から画像に電場ベクトルを挿入するプログラムを実行中
-HDF5ファイルの理解が進んだ
--PythonでHDF5ファイルを扱う際のh5pyライブラリの使い方を学んだ
明日やること
-S字のグラフェンから少し離れたところでのグラフェンからの反射した電場の観測
-アニメーションgifからmp4への変換の検討
--PowerPointに張り付けたアニメーションgifを再生するとPCによっては固まることがあるため
**8月17日 [#od59ed4c]
やったこと
-S字のシミュレーションの実行
--電場ベクトルが回転する様子が確認できた
--条件を変えた場合のシミュレーションを実行中
分かったこと
-プログラムの理解
-1次元の計算なら(set! dimensions 1)とすると計算スピードが向上することが分かった
**8月16日 [#u58f2b52]
やったこと
-昨日、帰る前に動かしていたプログラムが途中で停止していたため、再度実行
--初めに１００万回計算させ、桁が10桁落ちるまで反復計算をさせるプログラム
--プログラムの半分で約7.5時間かかった
--１００００回計算させるものでは、30分弱で計算が終わった
-計算中の時間でMEEPでPythonが使えるようにした
--理由
---Gaucheではソースコードが括弧だらけで可読性が悪い
---MEEPはPythonもサポートしているし、Gaucheに比べてソースコードが見やすい
--方法
---[[PukiWikiへ:http://flex.phys.tohoku.ac.jp/japanese/pukiwiki/index.php?%C5%C5%BC%A7%B3%A6%A5%B7%A5%DF%A5%E5%A5%EC%A1%BC%A5%B7%A5%E7%A5%F3#fc543fa9]]
あしたやること
-S字のプログラムの実行
**8月15日 [#i9f3a00c]
やったこと
-グラフェンの前後に誘電体を挟むことで
グラフェンの光吸収率を求めるプログラムの実行
--収束まで多くの反復計算が必要だったためまだ収束した結果は得られてない
--帰る直前に十分な反復計算をさせるプログラムを実行させる
**8月9日 [#p92f2ab1]
やったこと
-丸岡さんがもらったグラフェンの吸収率を求めるプログラムの実行と考察
--単位を合わせるためにシート伝導率にグラフェンの厚さを考慮する必要があるが
グラフェンの厚さはプログラム上の厚さでよいことが分かった
--2.36%という近しい値を得ることができた
-MEEPの理解が進んだ
--プログラムにコメントが多く入っていたため
明日やること
-グラフェンの間にいくつか誘電体を挟んだ状態でグラフェンの吸収率を求める
**8月8日 [#cb5bae23]
やったこと
-丸岡さんからもらったグラフェンの吸収率を求めるプログラムを参考にして、
グラフェンの吸収率を求めるプログラム作成（未完成）
本当にやるべきだったこと
-丸岡さんからもらったプログラムの実行と結果の考察
学んだこと
-人から言われたことを100%理解する、理解できなかったら理解できるまで質問する
-Yes,NoクエスチョンではYesかNoで答える
--×「〜はできましたか」ー「途中までやってたんですけど、〜のせいでまだ途中です」
--〇「〜はできましたか」ー「できませんでした。理由は〜です。」
-処理に使う変数や定数はプログラムの初めにすべて定義するよりも
処理の直前付近で定義した方が可読性がよくなる
--ただし、プログラムの核となるような変数はプログラムの初めて定義しておく
-プログラムには十分にコメントを入れる
-よいプログラムの条件
--正しい
--わかりやすい
--使いやすい
--柔軟性、汎用性が高い
--速い
--小さい
**8月7日 [#xed8c74a]
やったこと
-昨日(8月6日)の簡単なシミュレーションの結果の考察
-MEEPの理解(少しだけ)
-Gaucheの理解(少しだけ)
-flexへの接続完了
--グループポリシーの設定の変更
--詳しくは[[こちら:https://support.microsoft.com/de-ch/help/4046019/guest-access-smb2-disabled-by-default-in-windows-10-server-2016]]
明日の課題
-グラフェンの吸収率2.3%の再現シミュレーション
-出力画像にタイムステップの表示
**8月6日 [#c14efaf1]
やったこと
-Skypeのインストール
-MEEPをwindowsに入れるためにubuntuのインストール
-MEEPのインストール
-ubuntuにemacsとgaucheのインストール
-Xmingのインストール、Xmingでflexにログイン
-MEEPを使った簡単なシミュレーションを実行した
困っていること
-windowsとflex間でファイル転送ができない
--windows上にflexのショートカットを作成しようとするもflexが見つからないとエラーがでる
![[PukiWiki]](image/c60.png "[PukiWiki]")
