|
2004/11/19 (金) SPRシミュレーター
相変わらず化学系らしからぬ毎日を送っています。これも何かの運命?というわけで、表面プラズモン共鳴の、Kretchmann configurationにおける3層のFresnel equationを使って、角度に対する光吸収をシミュレーションしてみることにしました。 とりあえず使い慣れたExcelVBAでやってみようと、式や変数定義といったプログラムを100行ほどで書き終えて、いざ計算実行!パラメータとしては、入射光の角振動数、光の速度、ガラスプリズム、金薄膜、センシング層(液体)それぞれの誘電率が入ります。・・・ところが、とある問題にぶつかりました。パラメータによっては計算を実行するとエラーが出て停止してしまうことがあり、とりあえず変数をみながら原因を探っていきますと・・・なんと負の数の平方根を求めようとしておりました。そう、その式は複素演算をしなければならないのです。VBAでは複素数の計算式や変数は扱えないので、やろうと思うと配列か何かでさらに自分で式を組み立てていく必要があり、時間×能力<やるべき仕事になってしまうので、そこはExcelの複素計算機能のお力を借りることにしました。 頭の痛くなるような式をいくつか書き、計算実行・・・あれ、うまくいかない・・・あ、ここ間違ってた・・・こんどこそ・・・・う〜ん・・・・・・結局まだできてません・・・
|