|
昨年度までの研究で修得した量子渦の交流電場応答についての数値計算手法を発展させた。この数値計算では、Keldysh形式の準古典グリーン関数を超伝導の秩序変数、不純物自己エネルギー、及び、電気化学ポテンシャルと自己無撞着に解く必要がある。しかし、秩序変数と電気化学ポテンシャルを求めるためには、準粒子励起のエネルギーについて超伝導ギャップよりも十分に大きいエネルギーカットオフまでエネルギー積分を実行する必要があり、従来の計算方法では高エネルギーまでの緻密なKeldyshグリーン関数を要する膨大な数値計算が必要であった。本年度はこのような数値計算の困難を解決するために計算方法の改善を行った。新たに開発した数値計算方法では、Keldyshグリーン関数を松原形式のグリーン関数への解析接続が可能な遅延グリーン関数、先進グリーン関数と解析接続できない異常関数に分ける。遅延グリーン関数と先進グリーン関数は松原グリーン関数に解析接続できるので、エネルギー積分を実行する代りに有限個の松原周波数についての和をとればよい。異常関数についてはエネルギー積分を実行する必要があるが、低温では交流電場の周波数以下のエネルギーを持つ準粒子にしか異常関数は寄与しないので、エネルギー積分の積分範囲を小さくできる。このようにして、秩序変数と電気化学ポテンシャルを松原グリーン関数の有限和とKeldyshグリーン関数の異常関数部分についての交流周波数までのエネルギー積分から計算することで、数値計算コストを抑えることができる。この数値計算方法は国内会議において発表しており、関連の深い研究者と情報交換を行っている。
|
