(57)【特許請求の範囲】
【請求項１】
コンピュータによって実行される電気伝導特性予測方法であって、該コンピュータが、
バンド計算によって、記憶部に記憶された原子空孔欠陥を有するグラフェンシートの構造モデルにおける逆格子空間中の任意のｋ点において伝導方向のバンド構造を取得し、
所定のエネルギー範囲にある各バンドの伝導方向のｋ点間の傾きの最大値を取得し、
原子空孔欠陥のない理想的なグラフェンシートの前記所定のエネルギー範囲にあるバンドの傾きで該最大値を割ることによって透過係数を算出し、
前記伝導方向に垂直な方向のｋ点毎に前記バンド構造と前記最大値とを取得して算出した前記透過係数を合計して該垂直な方向のｋ点の総数で規格化し、
規格化された前記透過係数を所定の印加電圧に対応するエネルギー範囲で積分することによって電流情報を取得して前記記憶部に出力する
ことを特徴とする電気伝導特性予測方法。
【請求項２】
所定の原子空孔欠陥を有するグラフェンシートについて、二次元周期で欠陥が並んだ前記構造モデルを作成して、前記記憶部に記憶する
ことを特徴とする請求項１記載の電気伝導特性予測方法。
【請求項３】
前記構造モデルの伝導方向の周期サイズは３４ｎｍ以上である
ことを特徴とする請求項２記載の電気伝導特性予測方法。
【請求項４】
前記伝導方向に垂直な方向のｋ点数を、前記伝導方向の逆格子ベクトル上にｋ点が約０．０４Ｂｏｈｒ−１の間隔で並ぶように定められている
ことを特徴とする請求項１乃至３のいずれか一項記載の電気伝導特性予測方法。
【請求項５】
バンド計算によって、記憶部に記憶された原子空孔欠陥を有するグラフェンシートの構造モデルにおける逆格子空間中の任意のｋ点において伝導方向のバンド構造を取得し、
所定のエネルギー範囲にある各バンドの伝導方向のｋ点間の傾きの最大値を取得し、
原子空孔欠陥のない理想的なグラフェンシートの前記所定のエネルギー範囲にあるバンドの傾きで該最大値を割ることによって透過係数を算出し、
前記伝導方向に垂直な方向のｋ点毎に前記バンド構造と前記最大値とを取得して算出した前記透過係数を合計して該垂直な方向のｋ点の総数で規格化し、
規格化された前記透過係数を所定の印加電圧に対応するエネルギー範囲で積分することによって電流情報を取得して前記記憶部に出力する
処理をコンピュータに実行させるプログラム。
【請求項６】
前記所定のエネルギー範囲は、フェルミ準位を含み、理想的なグラフェンシートの前記透過係数が一定となるエネルギー範囲であることを特徴とする請求項１乃至４のいずれか一項記載の電気伝導特性予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、グラフェンシートを使った電子材料の電気伝導特性を予測する技術に関する。
【背景技術】
【０００２】
グラフェンは、例えば、黒鉛結晶中において炭素の六角形格子を構成するｓｐ２結合をした炭素原子よりなる原子層であるが、散乱の効果を抑制できれば室温でも２０００００ｃｍ2Ｖ-1ｓ-1を超える非常に大きな電子移動度を達成可能であることから、グラフェンのシートを使って超高速電子装置を作製する研究がなされている。
【０００３】
グラフェンシートはその製造過程において原子空孔欠陥が導入される可能性がある。また、非特許文献1のように意図的に原子空孔欠陥を導入する場合もある。そのような欠陥を含むグラフェンシートの電気伝導特性を知ることは電子装置へ応用する上で非常に重要である。
【０００４】
電子材料の原子レベルでの電気伝導特性を調べる方法としては、非平衡グリーン関数法に基づいた量子伝導計算が従来用いられる。この方法は、（１）まず半無限のリード電極で特性を調べたい対象の系を挟み込んだ構造のモデルを用意する。（２）次に、印加電圧に対応するポテンシャル勾配をそのリード電極間に発生させ、セルフコンシステントにポテンシャル電子状態を収束させる。（３）それにより得られたグリーン関数を用い、ランダウア公式に従って、逆格子空間の各ｋ点における伝導方向の透過係数を求める。（４）これをすべてのｋ点について求めて合計し、電流値を計算する。印加電圧を変えて以上の（２）−（４）を繰り返すことで電流電圧特性（電気伝導特性）を得ることができる。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】J. Bai et al., "Graphen nanomesh", Nature Nanotechnology, March 2010, vol.5, p.190
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上述した非平衡グリーン関数法に基づいた量子伝導計算の方法では、特に（２）のセルフコンシステントにポテンシャルと電子状態を収束させるプロセス、及びグリーン関数を作るプロセスに高コストの計算量を要する。
【０００７】
例えば、原子空孔欠陥を含むグラフェンシートの電気伝導特性を調べるためには、少なくとも１００原子程度から成る構造モデルを用意する必要があるが、その場合に任意の印加電圧下での電流値を十分な精度で計算するためには、例えば約１００ＭＰＩ並列の計算機を使って約２時間を要する。
【０００８】
よって、本発明の目的は、様々な欠陥パターンによる電気伝導特性の違いを探索したい場合や、更に多くの原子数を必要とする構造モデルを用いたい場合には、上述したような計算方法より低コストで電気伝導特性を予測することである。
【課題を解決するための手段】
【０００９】
開示の技術は、コンピュータによって実行される電気伝導特性予測方法であって、該コンピュータが、バンド計算によって、憶部に記憶された原子空孔欠陥を有するグラフェンシートの構造モデルにおける逆格子空間中の任意のｋ点において伝導方向のバンド構造を取得し、所定のエネルギー範囲にある各バンドの伝導方向のｋ点間の傾きの最大値を取得し、原子空孔欠陥のない理想的なグラフェンシートの前記所定のエネルギー範囲にあるバンドの傾きで該最大値を割ることによって透過係数を算出し、前記伝導方向に垂直な方向のｋ点毎に前記バンド構造と前記最大値とを取得して算出した前記透過係数を合計して該垂直な方向のｋ点の総数で規格化し、規格化された前記透過係数を所定の印加電圧に対応するエネルギー範囲で積分することによって電流情報を取得して前記記憶部に出力する。
【００１０】
また、上記課題を解決するための手段として、コンピュータに電気伝導特性を予測させるためのプログラム、そのプログラムを記録した記録媒体、及び、電気伝導特性予測装置とすることもできる。
【発明の効果】
【００１１】
開示の技術では、グラフェンシート中に原子空孔欠陥を含む種々の構造に関して、グラフェンシートの電流電圧特性を関連技術より低コストの計算で電気伝導特性を予測することができる。
【図面の簡単な説明】
【００１２】
【図１】コンピュータ装置のハードウェア構成を示す図である。
【図２】電気伝導特性を予測するための機能構成例を示す図である。
【図３】電気伝導特性予測処理における動作例を示す図である。
【図４】手順（２）での処理を説明するためのフローチャート図である。
【図５】手順（３）及び手順（４）での処理を説明するためのフローチャート図である。
【図６】手順（５）及び手順（６）での処理を説明するためのフローチャート図である。
【図７】第１の実施例において調査したい欠陥を含むグラフェンシートを用いた電子装置の模式図である。
【図８】第１の実施例における欠陥が二次元周期で並んだ構造を示す平面図である。
【図９】グラフェンの伝導方向のバンド構造を説明するための図である。
【図１０】透過係数の比較を示す図である。
【図１１】電流電圧特性の計算結果の比較例を示す図である。
【図１２】第１の実施例において、二次元周期構造の伝導方向周期ｎＡに対する電流電圧特性の依存性を示す図である。
【図１３】第１の実施例において、伝導方向に垂直な方向のｋ点数ｎｋＺに対する電流電圧特性の依存性を示す図である。
【図１４】第２の実施例に係る、調査したい欠陥を含むグラフェンシートを用いた電子装置の模式図である。
【図１５】第２の実施例における電流電圧特性の計算結果を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明の実施の形態を図面に基づいて説明する。本実施の形態に係る、原子空孔欠陥を含むグラフェンシートの電気伝導特性予測の方法について、以下に箇条書きで順を追って説明する。
【００１４】
（１）調べたい（解析対象の）原子空孔欠陥を含むグラフェンシートについて、二次元周期で欠陥が並んだ構造モデルを作成する。
【００１５】
（２）バンド計算を行い、逆格子空間中の任意のk点において伝導方向のバンド構造を求める。
【００１６】
（３）フェルミ準位ＥＦ前後０．８ｅＶ、即ち、
【００１７】
【数１】
0005928169.tif 000002
程度のエネルギー範囲にある各バンドについて傾きｄＥ／ｄｋの最大値を求め、ｄｂ（Ｅ）とする。
【００１８】
（４）バンド計算により予め求めておいた原子空孔欠陥のない理想的なグラフェンシートのフェルミ準位ＥＦ前後０．８ｅＶ、即ち、
【００１９】
【数２】
0005928169.tif 000003
のバンドの傾きｄａと手順（３）で求めたｄｂ（Ｅ）を用いて、透過係数Ｔｋ（Ｅ）をＴｋ（Ｅ）＝ｄｂ（Ｅ）／ｄａの式に従い算出する。
【００２０】
ここで、発明者は、バンド計算によって得られる理想的なグラフェンシートのバンドがＥＦ付近で直線的であり傾きが一定であることと、関連技術の方法で計算した理想的なグラフェンシートの透過係数がＥＦ付近で１となることに着目している。
【００２１】
伝導方向に垂直な方向のｋ点全てについて（２）から（４）の手順を行う。
【００２２】
（５）伝導方向に垂直な方向のｋ点全てについて（２）から（４）の手順を行った後、得られた透過係数Ｔｋ（Ｅ）を合計してｋ点数ｎｋＺで規格化する。
【００２３】
【数３】
0005928169.tif 000004
（６）Ｔｔｏｔ（Ｅ）を任意の印加電圧に対応するエネルギー範囲で積分して電流値を計算し、電流電圧特性を得る。
【００２４】
上述した手順（２）から手順（４）により、任意の構造の欠陥を含むグラフェンシートの電流電圧特性を計算することが可能となる。手順（１）から手順（４）で透過係数を求める方法が新規技術である。求めた透過係数を用いて、手順（５）、手順（６）で電流電圧特性を計算すること、また、手順（１）、手順（２）を透過係数の計算に適用することは知られていない。
【００２５】
上記手順（１）から手順（６）による原子空孔欠陥を含むグラフェンシートの電気伝導特性予測処理は、図１に示されるようなハードウェア構成を有するコンピュータ装置１００によって実行される。
【００２６】
図１は、コンピュータ装置のハードウェア構成を示す図である。図１において、コンピュータ装置１００は、コンピュータによって制御される端末であって、ＣＰＵ（Central Processing Unit）１１と、主記憶装置１２と、補助記憶装置１３と、入力装置１４と、表示装置１５と、出力装置１６と、通信Ｉ／Ｆ（インターフェース）１７と、ドライブ１８とを有し、バスＢに接続される。
【００２７】
ＣＰＵ１１は、主記憶装置１２に格納されたプログラムに従ってコンピュータ装置１００を制御する。主記憶装置１２には、ＲＡＭ（Random Access Memory）等が用いられ、ＣＰＵ１１にて実行されるプログラム、ＣＰＵ１１での処理に必要なデータ、ＣＰＵ１１での処理にて得られたデータ等を格納する。また、主記憶装置１２の一部の領域が、ＣＰＵ１１での処理に利用されるワークエリアとして割り付けられている。
【００２８】
補助記憶装置１３には、ハードディスクドライブが用いられ、各種処理を実行するためのプログラム等のデータを格納する。補助記憶装置１３に格納されているプログラムの一部が主記憶装置１２にロードされ、ＣＰＵ１１に実行されることによって、各種処理が実現される。記憶部１３０は、主記憶装置１２及び／又は補助記憶装置１３を有する。
【００２９】
入力装置１４は、マウス、キーボード等を有し、ユーザがコンピュータ装置１００による処理に必要な各種情報を入力するために用いられる。表示装置１５は、ＣＰＵ１１の制御のもとに必要な各種情報を表示する。出力装置１６は、プリンタ等を有し、ユーザからの指示に応じて各種情報を出力するために用いられる。通信Ｉ／Ｆ１７は、例えばインターネット、ＬＡＮ（Local Area Network）等に接続し、外部装置との間の通信制御をするための装置である。
【００３０】
コンピュータ装置１００によって行われる処理を実現するプログラムは、例えば、ＣＤ−ＲＯＭ（Compact Disc Read-Only Memory）等の記憶媒体１９によってコンピュータ装置１００に提供される。即ち、プログラムが保存された記憶媒体１９がドライブ１８にセットされると、ドライブ１８が記憶媒体１９からプログラムを読み出し、その読み出されたプログラムがバスＢを介して補助記憶装置１３にインストールされる。そして、プログラムが起動されると、補助記憶装置１３にインストールされたプログラムに従ってＣＰＵ１１がその処理を開始する。尚、プログラムを格納する媒体としてＣＤ−ＲＯＭに限定するものではなく、コンピュータが読み取り可能な媒体であればよい。コンピュータ読取可能な記憶媒体として、ＣＤ−ＲＯＭの他に、ＤＶＤディスク、ＵＳＢメモリ等の可搬型記録媒体、フラッシュメモリ等の半導体メモリであっても良い。
【００３１】
また、コンピュータ装置１００によって行われる処理を実現するプログラムが、通信Ｉ／Ｆ１７を介して外部装置から提供されてもよい。或いは、外部装置へ該プログラムを提供し、後述される各処理は外部装置で実現されるように構成してもよい。通信Ｉ／Ｆ１７による通信は無線又は有線に限定されるものではない。
【００３２】
電気伝導特性を予測するための機能構成について、図２及び図３で説明する。図２は、電気伝導特性を予測するための機能構成例を示す図である。図３は、電気伝導特性予測処理における動作例を示す図である。
【００３３】
図２において、コンピュータ装置１００は、入力パラメータ取得部４１と、構造モデル作成部４２と、バンド構造算出部４３と、透過係数算出部４４と、電流電圧特性算出部４５とを有する。入力パラメータ取得部４１と、構造モデル作成部４２と、バンド構造算出部４３と、透過係数算出部４４と、電流電圧特性算出部４５とは、ＣＰＵ１１が対応するプログラムを実行することによる処理によって実現される。
【００３４】
また、コンピュータ装置１００の記憶部１３０には、入力パラメータ５１と、構造モデル５２と、規定値５３と、出力データ５４とが格納される。
【００３５】
入力パラメータ取得部４１は、電気伝導特性予測処理を行う際に必要となる入力パラメータ５１を、ユーザから取得する処理部である。入力パラメータ取得部４１は、入力パラメータ５１を取得するために、入力パラメータ５１を設定可能なユーザインタフェースを表示装置１５に表示する。ユーザから取得した入力パラメータ５１は、記憶部１３０に格納される。
【００３６】
構造モデル作成部４２は、記憶部１３０に格納された入力パラメータ５１に基づいて、構造モデル５２を作成する処理部である。入力パラメータ５１に基づいて、ユーザが調べたい（解析対象の）原子空孔欠陥を含むグラフェンシートについて、二次元周期で欠陥が並んだ構造モデル５２が作成され、記憶部１３０に格納される。構造モデル作成部４２での処理は、上記手順（１）に相当する。
【００３７】
バンド構造算出部４３は、構造モデル５２を用いて、逆格子空間中の任意のk点において伝導方向のバンド構造を算出する処理部である。バンド構造算出部４３は、上記手順（２）に相当する。バンド構造算出部４３によって、原子空孔欠陥のない理想的なグラフェンシートのフェルミ準位ＥＦ前後０．８ｅＶのバンドの傾きｄａが算出される。
【００３８】
透過係数算出部４４は、透過係数Ｔｋ（Ｅ）を算出する処理部であり、上記手順（３）及び（４）に相当する。透過係数算出部４４は、フェルミ準位ＥＦ前後０．８ｅＶ程度のエネルギー範囲にある各バンドについて傾きｄＥ／ｄｋの最大値ｄｂ（Ｅ）を求める（手順（３））。そして、透過係数算出部４４は、バンド構造算出部４３によって算出されたバンドの傾きｄａと、手順（３）で求めた傾きｄＥ／ｄｋの最大値ｄｂ（Ｅ）とを用いて、透過係数Ｔｋ（Ｅ）をＴｋ（Ｅ）＝ｄｂ（Ｅ）／ｄａの式に従い算出する（手順（４））。
【００３９】
図３に示すように、伝導方向に垂直な方向のｋ点全てについて、ｋ点毎に、バンド構造算出部４３と透過係数算出部４４とによって（２）から（４）の手順が行われる。得られた透過係数Ｔｋ（Ｅ）は合計され、ｋ点数ｎｋＺで規格化される。
【００４０】
電流電圧特性算出部４５は、規格化された透過係数Ｔｔｏｔ（Ｅ）を、入力パラメータ５１で指定される印加電圧に対応するエネルギー範囲で積分して電流値を計算し、電流電圧特性を得る処理部である。電流電圧特性算出部４５は、上記手順（６）に相当する。電流電圧特性算出部４５による算出結果（電流密度Ｃｄｅｎ）を含む出力データ５４が記憶部１３０に出力され格納される。
【００４１】
入力パラメータ５１には、原子座標、アームチェア方向のユニットセルサイズｎＡ（＝８０）、ジグザグ方向のユニットセルサイズｎＺ、伝導方向（アームチェア方向）のｋ点数ｎｋＡ（＝１０）、伝導方向に垂直な方向（ジグザグ方向）のｋ点数ｎｋＺ（〜ｎＺ／３５）、電流電圧特性予測に用いるバンドの本数ｎｂ（〜ｎＡ＊ｎＺ／１０）、調べたい印加電圧の刻みｄＶｂ（＝０．１Ｖ）、及び、調べたい印加電圧の最大値Ｖｂ＿ｍａｘ（＝１．５Ｖ）等が含まれる。ここで、（  ）内の値は推奨値であるが、ユーザによって変更可能である。
【００４２】
入力パラメータ５１のうち、原子座標、ユニットセルサイズｎＡ及びｎＺ、ｋ点数ｎｋＡ及びｎｋＺ等が、構造モデル作成部４２によって構造モデル５２が作成される際に、参照される。
【００４３】
入力パラメータ５１のうち、電流電圧特性予測に用いるバンドの本数ｎｂ、調べたい印加電圧の刻みｄＶｂ、及び、調べたい印加電圧の最大値Ｖｂ＿ｍａｘ等が、電流電圧特性算出部４５によって参照される。
【００４４】
また、規定値５３には、グラフェンのＣ−Ｃボンド長
【００４５】
【数４】
0005928169.tif 000005
、強束縛近似計算で用いるホッピング積分値「２．５５ｅＶ」、欠陥のないグラフェンのフェミル準位付近のバンドの傾きｄa＝１０ｅＶ  Ｂｏｈｒ等が含まれる。規定値５３は、透過係数算出部４４による透過係数Ｔｋ（Ｅ）を算出する際に参照される。
【００４６】
出力データ５４は、電流電圧特性算出部４５による算出結果を示し、電流密度Ｃｄｅｎ等の値を含む。
【００４７】
以下に、手順（２）から（６）に係る処理ついて、図４から図６で説明する。図４は、手順（２）での処理を説明するためのフローチャート図である。図５は、手順（３）及び手順（４）での処理を説明するためのフローチャート図である。図６は、手順（５）及び手順（６）での処理を説明するためのフローチャート図である。
【００４８】
図４及び図５において、伝導方向に垂直な方向（ジグザグ方向）の各ｋ点毎（ｉｋｚ番目毎）に、ＣＰＵ１１は、バンド構造算出部４３を起動して、構造モデル５２におけるジグザグ方向のｋ点での手順（２）の処理を行い、更に、透過係数算出部４４を起動して、手順（３）及び手順（４）での処理を行う。図６において、ＣＰＵ１１は、電流電圧特性算出部４５を起動して、手順（５）及び手順（６）での処理を行う。
【００４９】
先ず、図４において、バンド構造算出部４３による処理が開始され、バンド構造算出部４３はｊにゼロを設定して初期化する（ステップＳ１１）。そして、バンド構造算出部４３は、強束縛近似計算によって、ｍ番目のバンドのｉｋＡ番目の逆格子点（０≦ｋ≦πの範囲、全部でｎｋＡ個）夫々について固有エネルギーｅｉｇを計算する（ステップＳ１２）。
【００５０】
バンド構造算出部４３は、計算した固有エネルギーｅｉｇがゼロより大であるか否かを判断する（ステップＳ１３）。固有エネルギーｅｉｇがゼロより大であると判断した場合、ｊを１インクリメントして（ステップＳ１４）、ｅ（ｊ，ｉｋＡ）に固有エネルギーｅｉｇを設定する（ステップＳ１５）。そして、バンド構造算出部４３は、ステップＳ１６へと進む。
【００５１】
一方、ステップＳ１３において、固有エネルギーｅｉｇがゼロ以下であると判断した場合、バンド構造算出部４３は、最後のｍ番目のバンドであるか否かを判断する（ステップＳ１６）。最後のｍ番目のバンドでない場合、バンド構造算出部４３は、ｍを１インクリメントして、ステップＳ１２へと戻り、上記同様の処理を行う。
【００５２】
一方、ステップＳ１６において、最後のｍ番目のバンドであると判断した場合、バンド構造算出部４３による手順（２）での処理は終了し、手順（３）及び手順（４）による処理へと進む（図５）。
【００５３】
図５中、ステップＳ２３が手順（３）に相当し、ステップＳ２８が手順（４）に相当する。図５において、透過係数算出部４４による処理が開示され、透過係数算出部４４はバンド数ｉｂ及びｋ点数ｉｋＡに１を設定する（ステップＳ２１）。
【００５４】
透過係数算出部４４は、傾きの最大値ｄｂ＿ｍａｘ（ｉｂ）及びエネルギー値ｅｎｅ（ｉｋｚ，ｉｂ）をゼロに初期化する（ステップＳ２２）。
【００５５】
そして、透過係数算出部４４は、現在のｋ点（ｉｋＡ）と次のｋ点（ｉｋＡ＋１）との距離ｄｋＡと、傾きｄｂとを下記式により計算する（ステップＳ２３）。
【００５６】
【数５】
0005928169.tif 000006
【００５７】
【数６】
0005928169.tif 000007
傾きｄｂは、次のｋ点（ｉｋＡ＋１）のエネルギー値ｅ（ｉｂ，ｉｋＡ＋１）と現在のｋ点（ｉｋＡ）のエネルギー値ｅ（ｉｂ，ｉｋＡ）との差を距離ｄｋＡで割ることにより求められる。
【００５８】
次に、透過係数算出部４４は、傾きｄｂが最大値ｄｂ＿ｍａｘ（ｉｂ）より大であるか否かを判断する（ステップＳ２４）。大である場合、透過係数算出部４４は、最大値ｄｂ＿ｍａｘ（ｉｂ）に傾きｄｂを設定し、ｅｎｅ（ｉｋｚ，ｉｂ）にｅ（ｉｂ，ｉｋＡ）を設定して（ステップＳ２５）、ステップＳ２６へと進む。
【００５９】
【数７】
0005928169.tif 000008
【００６０】
【数８】
0005928169.tif 000009
【００６１】
一方、ステップＳ２４において、傾きｄｂが最大値ｄｂ＿ｍａｘ（ｉｂ）以下であると判断した場合、透過係数算出部４４は、現在のｋ点（ｉｋＡ）がｋ点数ｎｋＡから１減算した値に等しいか否かを判断する（ステップＳ２６）。等しくない場合、透過係数算出部４４は、ｉｋＡを１インクリメントして（ステップＳ２７）、ステップＳ２３へ戻り、上記同様の処理を繰り返す。
【００６２】
一方、等しい場合、透過係数算出部４４は、伝導方向に垂直な方向のｋ点（ｉｋＺ）における透過係数ｔｒを計算する（ステップＳ２８）。透過係数ｔｒは、最大値ｄｂ＿ｍａｘ（ｉｂ）をグラフェンのフェミル準位付近のバンドの傾きｄaで割ることによって得られる。
【００６３】
【数９】
0005928169.tif 000010
【００６４】
そして、透過係数算出部４４は、バンド番号ｉｂがバンド総数ｎｂと等しいか否かを判断する（ステップＳ２９）。バンド番号ｉｂがバンド総数ｎｂと等しくない場合、バンド総数ｎｂに達していないため、透過係数算出部４４は、バンド番号ｉｂを１インクリメントして（ステップＳ３０）、ステップＳ２２へと戻り、上記同様の処理を繰り返す。
【００６５】
一方、バンド番号ｉｂがバンド総数ｎｂと等しい場合、伝導方向に垂直なジグザグ方向のｋ点番号ｉｋＺでの透過係数算出部４４による処理は、ＣＰＵ１１によって一旦終了する。ＣＰＵ１１は、伝導方向に垂直なジグザグ方向のｋ点番号ｉｋＺがジグザグ方向のｋ点の総数ｎｋＺ以下であるか否かを判断する（ステップＳ３１）。ジグザグ方向のｋ点番号ｉｋＺが総数ｎｋＺ以下である場合、ＣＰＵ１１は、ステップＳ１１（図４）へ戻り、上述したように、バンド構造算出部４３を起動して、次のジグザグ方向のｋ点での手順（２）の処理を行わせ、更に、透過係数算出部４４を起動して、手順（３）及び手順（４）での処理を行わせる。
【００６６】
一方、ジグザグ方向のｋ点番号ｉｋＺが総数ｎｋＺを超えた場合、ＣＰＵ１１は、電流電圧特性算出部４５を起動して、ステップＳ４１（図６）へと進む。
【００６７】
図６において、電流電圧特性算出部４５による処理が開始され、電流電圧特性算出部４５は、ｉＶｂに１を設定し（ステップＳ４１）、電流値Ｃｕｒｒにゼロを設定して初期化する（ステップＳ４２）。
【００６８】
そして、電流電圧特性算出部４５は、エネルギー値ｅｎｅ（ｉｋｚ，ｉｂ）がｄＶｂ＊ｉＶｂ／２以下であるか否かを判断する（ステップＳ４３）。ｄＶｂは、記憶部１３０に格納されている入力パラメータ５１から取得する。ｉＶｂは、印加電圧の刻み回数である。
【００６９】
エネルギー値ｅｎｅ（ｉｋｚ，ｉｂ）がｄＶｂ＊ｉＶｂ／２より大きい場合、電流電圧特性算出部４５は、電流値Ｃｕｒｒに、ｉｋｚでの透過係数とエネルギー値とに基づくバンド間のエネルギー値を付加して（ステップＳ４４）、ステップＳ４５へと進む。一方、ｄＶｂ＊ｉＶｂ／２以下の場合、電流電圧特性算出部４５は、ステップＳ４５へと進む。
【００７０】
電流電圧特性算出部４５は、バンド番号ｉｂがバンド総数から１減算した値と等しいか否かを判断する（ステップＳ４５）。等しくない場合、バンド番号ｉｂを１インクリメントして（ステップＳ４６）、ステップＳ４３へと戻り、上記同様の処理を繰り返す。
【００７１】
一方、等しい場合、電流電圧特性算出部４５は、ジグザグ方向のｋ点番号ｉｋｚがジグザグ方向のｋ点総数ｎｋｚに等しいか否かを判断する（ステップＳ４７）。等しくない場合、電流電圧特性算出部４５は、ジグザグ方向のｋ点番号ｉｋｚを１インクリメントして（ステップＳ４８）、ステップＳ４３へと戻り、上記同様の処理を繰り返す。
【００７２】
一方、ジグザグ方向のｋ点総数ｎｋｚに等しい場合、電流電圧特性算出部４５は、電流密度Ｃｄｅｎを算出する（ステップＳ４９）。
【００７３】
【数１０】
0005928169.tif 000011
電流密度Ｃｄｅｎは、上記数１０によって求められる。
【００７４】
電流電圧特性算出部４５は、印加電圧Ｖｂ（＝ｄＶｂ＊ｉＶｂ）を計算し、計算した印加電圧と、ステップＳ４９で算出した電流密度Ｃｄｅｎとを含む出力データを記憶部１３０に格納する（ステップＳ５０）。
【００７５】
そして、電流電圧特性算出部４５は、印加電圧Ｖｂが、記憶部１３０に格納されている入力パラメータ５１の調べたい印加電圧の最大値Ｖｂ＿ｍａｘより小さいか否かを判断する（ステップＳ５１）。最大値Ｖｂ＿ｍａｘより小さい場合、電流電圧特性算出部４５は、印加電圧の刻み回数ｉＶｂを１インクリメントして（ステップＳ５２）、ステップＳ４２へ戻り、上述同様の処理を繰り返す。一方、印加電圧Ｖｂが最大値Ｖｂ＿ｍａｘ以上を示す場合、ＣＰＵ１１は、電流電圧特性算出部４５による処理は終了し、この電気伝導特性予測処理全体を終了する。
【００７６】
上述した電気伝導特性予測処理による実施例を以下に説明する。
【００７７】
［第１の実施例］
図７は、第１の実施例において調査したい欠陥を含むグラフェンシートを用いた電子装置の模式図である。第１の実施例では、グラフェンシート１０Ｇ中に欠陥１０Ａが存在し（構造ｂ）、その欠陥１０Ａを図７の左右から挟み込む形で設けられたグラフェン電極１０Ｓとグラフェン電極１０Ｄとの間に任意の電圧が印加された場合の両電極間に流れる電流を、上述した手順（１）から（６）に従って計算する。
【００７８】
伝導方向は、図７の左右方向で、グラフェンシート１０Ｇのアームチェア方向である。なお、図７、及び、以降の図では、簡単のために欠陥１０Ａの周りのＣ原子は未結合手を終端せずに描いているが、第１の実施例中では未結合手は全て終端されているものとして計算している。
【００７９】
手順（１）において、構造モデル作成部４２によって作成された、二次元周期で欠陥１０Ａが並んだ構造モデルを図８に示す。手順（２）以降は、図８中に示したユニットセル１１について計算が行われる。後述の通り、このユニットセル１１の伝導方向、即ち、アームチェア方向の周期サイズｎＡは精度を決定するパラメータであるが、ここではｎＡ＝８０（長さ約３４ｎｍ）を用いる。図８のユニットセル１１はｎＡ＝５である。
【００８０】
手順（２）のバンド計算の手法としては、例えば、上述したように強束縛近似法を用いる。第一原理法等のその他のバンド計算手法を用いてもよいが、計算コスト面で強束縛近似法を用いるのが最も有利である。強束縛近似計算で用いる基底は炭素原子当たり１つのｐｚ軌道で、相互作用は最近接原子間のみを考慮する。第一原理計算から求めた理想的なグラフェンのバンド構造を再現するように、ホッピング積分の値は２．５５ｅＶを採用する。
【００８１】
図９は、グラフェンの伝導方向のバンド構造を説明するための図である。ユニットセル１１と同じセルサイズで欠陥１０Ａのない理想的なグラフェン（構造ａ）と、欠陥１０Ａを含むグラフェン（構造ｂ）のユニットセル１１とについて求めた、Ｋ（０，２／３π，０）点からＫＸ（π，２／３π，０）点までの伝導方向のバンド構造をそれぞれ図９（Ａ）及び図９（Ｂ）に示す。ここで用いている強束縛近似計算によって得られるグラフェンのバンド構造はＥＦを軸に対称である。
【００８２】
図９（Ａ）に示したように、理想的なグラフェン（構造ａ）は０＜Ｅ−ＥＦ＜０．８のエネルギー範囲では傾きｄＥ／ｄｋが一定であり、その値をｄａとする。第１の実施例の場合、ｄａ＝１０ｅＶ  Ｂｏｈｒである。
【００８３】
一方、図９（Ｂ）に示した、欠陥１０Ａを含むグラフェンは、ＥＦ付近で傾きがゼロに近いフラットなバンドを持ち、エネルギーが上がるに連れて傾きｄａに近いバンドを持つようになり、傾きの最大値ｄｂを示すようになる。
【００８４】
図９（Ｂ）で得られた各バンドにおいて、傾きが最大となるエネルギー値Ｅ−ＥＦとその時の傾きｄｂ（Ｅ−ＥＦ）を求める。これが手順（３）である。
【００８５】
次に、手順（４）として、透過係数Ｔｋ（Ｅ−ＥＦ）をＴｋ（Ｅ−ＥＦ）＝ｄｂ（Ｅ−ＥＦ）／ｄａの式に従って算出する。
【００８６】
図１０は、透過係数の比較を示す図である。図１０中、縦軸に透過係数Ｔｋ（Ｅ−ＥＦ）を示し、横軸にエネルギーＥ−ＥＦ（ｅＶ）を示す。図１０において、理想的なグラフェン（（構造ａ）Ｉｄｅａｌ）、本第１の実施例（（構造ｂ）Ｎｅｗ）、及び、本第１の実施例と同じ構造について関連技術を用いた場合（（構造ｂ）Ｒｅｌａｔｅｄ  Ａｒｔ）の夫々の透過係数Ｔｋ（Ｅ−ＥＦ）の計算結果が示されている。理想的なグラフェン（（構造ａ）Ｉｄｅａｌ）は、図９（Ａ）で示されるバンドの傾きｄａを示している。
【００８７】
図１０から、本実施の形態に係る計算によって得られる透過係数は、関連技術である非平衡グリーン関数法に基づいた量子伝導計算による結果をよく再現することが分かる。第１の実施例において、図１０にＩｄｅａｌとして示したように理想的なグラフェン（構造Ａ）が０＜Ｅ−ＥＦ＜０．８のエネルギー範囲でＴｋ（Ｅ−ＥＦ）＝１であることを利用している。
【００８８】
（構造ｂ）関連技術では、図１０の透過係数の算出結果を得るのにＡＭＤ  Ｏｐｔｅｒｏｎ  ２．４ＧＨｚ  ４コア×２の構成のコンピュータ装置を使って１０４ＭＰＩ並列で１２０分を要しているが、第１の実施例では、ＡＭＤ  Ｏｐｔｅｒｏｎ  ３．０ＧＨｚ  ２コア×２の構成のマシン（以降マシンＡと呼ぶことにする）を使って１ＭＰＩ並列で８分しかかかっていない。単純に見積もって、計算コストを（３．０＊１＊８）／（２．４＊１０４＊１２０）＝１／１２４８にできたと言える。
【００８９】
最後に、手順（５）及び手順（６）において、以上の手順から得た透過係数を使って、電流電圧特性が算出される。手順（５）で必要となる伝導方向に垂直な方向のｋ点数ｎｋＺは後述の通り精度を決定するパラメータであるが、ここではｎｋＺ＝７を用いる。手順（６）で、例えば印加電圧Ｖｂ＝１Ｖのときの電流密度を求めたい場合は、
【００９０】
【数１１】
0005928169.tif 000012
を
【００９１】
【数１２】
0005928169.tif 000013
のエネルギー範囲で積分し、その積分値をｕｐスピンとｄｏｗｎスピンの両方を考慮して２倍する。
【００９２】
第１の実施例ではＥＦを軸に対称な透過係数であるから、０＜Ｅ−ＥＦ＜０．５でＴｔｏｔ（Ｅ−ＥＦ）の積分を行い、その積分値を４倍した値をＶｂ＝１Ｖのときの電流密度とする。
【００９３】
このようにして得られた電流電圧特性を図１１に示す。図１１は、電流電圧特性の計算結果の比較例を示す図である。図１１中、縦軸に電流密度（μＡ／ｎｍ）を示し、横軸にバイアス電圧Ｖｂ（Ｖ）を示す。
【００９４】
図１１から、０＜Ｖｂ＜１Ｖの範囲においては、本第１の実施例（（構造ｂ）Ｎｅｗ）は関連技術（（構造ｂ）Ｒｅｌａｔｅｄ  Ａｒｔ）と誤差最大１割程度の精度で電流電圧特性を予測できていることが分かる。透過係数が求まってから後のプロセスに要する計算コストは、本第１の実施例と関連技術でほとんど差がない。すなわち、本第１の実施例は関連技術の約１／１０００の計算コストで図１１の結果を得ることができる。
【００９５】
以上のように、第１の実施例を用いることで、欠陥を含むグラフェンシートの電流電圧特性を予測する際に、計算精度をほとんど損なうことなく、大幅な計算コストの削減が可能である。
【００９６】
ここで、第１の実施例において精度を決定するパラメータｎＡ及びｎｋＺについて説明する。
【００９７】
先ず、手順（１）で用意する図８中のユニットセル１１の伝導方向、すなわちアームチェア方向の周期サイズｎＡを上から順に１６０、８０、４０と変化させたときの電流電圧特性を図１２に比較した。図１２は、第１の実施例において、二次元周期構造の伝導方向周期ｎＡに対する電流電圧特性の依存性を示す図である。図１２中、縦軸に電流密度（μＡ／ｎｍ）を示し、横軸にバイアス電圧Ｖｂ（Ｖ）を示す。
【００９８】
図１２より、ｎＡが大きいほど、図１０におけるエネルギーの刻みを細かくすることができ、精度が良くなる。ただし、本方法では計算コストのほとんどはバンド計算において固有値の対角化に要する時間であるため、ｎＡを大きくすると、（ｎＡ）３でスケーリングして計算コストが高くなる。
【００９９】
例えば、マシンＡで１ＭＰＩ並列の場合、ｎＡ＝１６０、８０、４０に対して、ｋ点１つあたりの透過係数を求める手順（４）終了までの計算時間はそれぞれ６０分、８分、１分である。ｎＡが２倍になると、約２３＝８倍の計算コストになっていることが分かる。第１の実施例では精度と計算コストのバランスからｎＡ＝８０を選択した。
【０１００】
次に、手順（５）で必要となる伝導方向に垂直な方向のｋ点数ｎｋＺを上から順に１３、７、５、３と変化させたときの電流電圧特性の比較を図１３に示す。図１３は、第１の実施例において、伝導方向に垂直な方向のｋ点数ｎｋＺに対する電流電圧特性の依存性を示す図である。図１３中、縦軸に電流密度（μＡ／ｎｍ）を示し、横軸にバイアス電圧Ｖｂ（Ｖ）を示す。
【０１０１】
図１３より、ｎｋＺが大きいほど精度が高いことが分かる。計算コストのスケーリングはｎｋＺに等倍である。透過係数を求める手順（４）までの計算はｋ点毎に独立して行えるため、簡単に並列処理できる。第１の実施例では精度と計算コストのバランスからｎｋＺ＝７を選択した。この場合、伝導方向の逆格子ベクトル上にｋ点が約０．０４Ｂｏｈｒ−１の間隔で並ぶことになる。
【０１０２】
［第２の実施例］
図１４は、第２の実施例に係る、調査したい欠陥を含むグラフェンシートを用いた電子装置の模式図である。図１４（Ａ）に構造ｃのグラフェンシート２０Ｇと、図１４（Ｂ）に構造ｄのグラフェンシート３０Ｇを示す。
【０１０３】
図１４（Ａ）に示すグラフェンシート２０Ｇは、欠陥２０Ａが存在する構造ｃを有し、欠陥２０Ａを左右から挟み込むようにグラフェン電極２０Ｓ及び２０Ｄが設けられる。図１４（Ｂ）においても同様に、グラフェンシート３０Ｇは欠陥３０Ａが存在する構造ｄを有し、欠陥３０Ａを左右から挟み込むようにグラフェン電極３０Ｓ及び３０Ｄが設けられる。
【０１０４】
第２の実施例では、図１４（Ａ）のグラフェンシート２０Ｇを挟むグラフェン電極２０Ｓと２０Ｄとの間、図１４（Ｂ）のグラフェンシート３０Ｇを挟むグラフェン電極３０Ｓと３０Ｄとの間に任意の電圧が印加された場合の両電極間に流れる電流を上述した手順（１）から手順（６）に従って計算する。
【０１０５】
伝導方向は、図１４（Ａ）及び図１４（Ｂ）の左右方向で、グラフェンシート２０Ｇ及び３０Ｇのアームチェア方向である。なお、図１４（Ａ）及び図１４（Ｂ）では、簡単のために、欠陥２０Ａ及び欠陥３０Ａの周りのＣ原子は、未結合手を終端せずに描いているが、第２の実施例中では未結合手は全て終端されているものとして計算している。
【０１０６】
精度を決定するパラメータとして、ｎＡ＝８０、ｎｋ＝７を用いて、第２の実施例と同様の方法で求めた電流電圧特性を図１５に示す。図１５は、第２の実施例における電流電圧特性の計算結果を示す図である。図１５中、縦軸に電流密度（μＡ／ｎｍ）を示し、横軸にバイアス電圧Ｖｂ（Ｖ）を示す。
【０１０７】
第１の実施例同様、第２の実施例での計算結果（Ｎｅｗ）は関連技術（Ｒｅｌａｔｅｄ  Ａｒｔ）による計算結果をよく再現していることが分かる。計算コストも第１の実施例とほぼ同じで、手順（４）終了までの計算時間は例えばマシンＡの１ＭＰＩ並列でどちらも８分である。
【０１０８】
上述したように、本実施の形態により、グラフェンシート中に原子空孔欠陥を含む構造について、関連技術による非平衡グリーン関数法に基づいた量子伝導計算を用いる方法と同程度の精度を保ちつつ計算コストを削減することができる。また、様々な欠陥パターンによる電気伝導特性の違いの探索や、数百原子以上から成る構造モデルを要する欠陥を含むグラフェンシートの電気伝導特性の調査も可能となる。
【０１０９】
本発明は、具体的に開示された実施例に限定されるものではなく、特許請求の範囲から逸脱することなく、種々の変形や変更が可能である。
【０１１０】
以上の第１から第２の実施例を含む実施形態に関し、更に以下の付記を開示する。
（付記１）
コンピュータによって実行される電気伝導特性予測方法であって、該コンピュータが、
バンド計算によって、前記記憶部に記憶された原子空孔欠陥を有するグラフェンシートの構造モデルにおける逆格子空間中の任意のｋ点において伝導方向のバンド構造を取得し、
所定のエネルギー範囲にある各バンドについて傾きの最大値を取得し、
原子空孔欠陥のない理想的なグラフェンシートの前記所定のエネルギー範囲にあるバンドの傾きで該最大値を割ることによって透過係数を算出し、
前記伝導方向に垂直な方向のｋ点毎に前記バンド構造と前記最大値とを取得して算出した前記透過係数を合計してｋ点数で規格化し、
所定の印加電圧に対応するエネルギー範囲で積分することによって電流情報を取得して前記記憶部に出力する
ことを特徴とする電気伝導特性予測方法。
（付記２）
所定の原子空孔欠陥を有するグラフェンシートについて、二次元周期で欠陥が並んだ前記構造モデルを作成して、前記記憶部に記憶する
ことを特徴とする付記１記載の電気伝導特性予測方法。
（付記３）
前記構造モデルの伝導方向の周期サイズは３４ｎｍ以上である
ことを特徴とする付記２記載の電気伝導特性予測方法。
（付記４）
前記伝導方向に垂直な方向のｋ点数を、前記伝導方向の逆格子ベクトル上にｋ点が約０．０４Ｂｏｈｒ−１の間隔で並ぶように定められている
ことを特徴とする付記１乃至３のいずれか一項記載の電気伝導特性予測方法。
（付記５）
バンド計算によって、記憶部に記憶された原子空孔欠陥を有するグラフェンシートの構造モデルにおける逆格子空間中の任意のｋ点において伝導方向のバンド構造を取得し、
所定のエネルギー範囲にある各バンドについて傾きの最大値を取得し、
原子空孔欠陥のない理想的なグラフェンシートの前記所定のエネルギー範囲にあるバンドの傾きで該最大値を割ることによって透過係数を算出し、
前記伝導方向に垂直な方向のｋ点毎に前記バンド構造と前記最大値とを取得して算出した前記透過係数を合計してｋ点数で規格化し、
所定の印加電圧に対応するエネルギー範囲で積分することによって電流情報を取得して前記記憶部に出力する
処理をコンピュータに実行させるプログラム。
（付記６）
原子空孔欠陥を有するグラフェンシートの構造モデルを記憶した記憶部と、
バンド計算によって、前記記憶部に記憶された前記構造モデルにおける逆格子空間中の任意のｋ点において伝導方向のバンド構造を取得するバンド構造算出部と、
所定のエネルギー範囲にある各バンドについて傾きの最大値を取得し、原子空孔欠陥のない理想的なグラフェンシートの前記所定のエネルギー範囲にあるバンドの傾きで該最大値を割ることによって透過係数を算出する透過係数算出部と、
前記伝導方向に垂直な方向のｋ点毎に、前記バンド構造算出部と透過係数算出部とによる処理を行って得られた透過係数を合計してｋ点数で規格化し、所定の印加電圧に対応するエネルギー範囲で積分することによって電流情報を取得して前記記憶部に出力する電流電圧特性算出部と
を有することを特徴とする電気伝導特性予測装置。
【符号の説明】
【０１１１】
１１          ＣＰＵ
１２          主記憶装置
１３          補助記憶装置
１４          入力装置
１５          表示装置
１６          出力装置
１７          通信Ｉ／Ｆ
１８          ドライブ
１９          記憶媒体
４１          入力パラメータ取得部
４２          構造モデル作成部
４３          バンド構造算出部
４４          透過係数算出部
４５          電流電圧特性算出部
５１          入力パラメータ
５２          構造モデル
５３          規定値
５４          出力データ
１００        コンピュータ装置
１３０        記憶部
【図１】
0005928169.tif 000014
【図２】
0005928169.tif 000015
【図３】
0005928169.tif 000016
【図４】
0005928169.tif 000017
【図５】
0005928169.tif 000018
【図６】
0005928169.tif 000019
【図７】
0005928169.tif 000020
【図８】
0005928169.tif 000021
【図９】
0005928169.tif 000022
【図１０】
0005928169.tif 000023
【図１１】
0005928169.tif 000024
【図１２】
0005928169.tif 000025
【図１３】
0005928169.tif 000026
【図１４】
0005928169.tif 000027
【図１５】
0005928169.tif 000028