【実施例】
【0018】
以下に本発明の実施例を示し、本発明をより具体的に説明する。なお、これらは説明のための単なる例示であって、本発明はこれらによって何ら制限されるものではない。
【0019】
本実施例で用いた計算の条件を以下にまとめる。
【0020】
(分子シミュレーション解析用のソフトウェア)
本計算には、第一原理電子状態計算ソフトウェアであるAccelrys社製Dmol
3 version 6.1を用いた。
【0021】
(触媒表面モデル)
触媒表面モデルとして、ルチル型酸化チタン(110)面のスラブモデル(本発明では「酸化チタンスラブモデル」と言うことがある)を用いた。酸化チタンスラブモデルは、(4unit x 2unit)の大きさの面を持ち、深さ方向に4層とし、3次元周期境界条件を課した。下層の2層の原子をルチル型酸化チタンの結晶位置に固定し計算を行った。例えば
図1を用いて説明すると、図の中での上面が触媒表面であり、酸素分子は触媒表面のチタン原子に吸着し、吸着酸素分子となる。
【0022】
(計算条件)
spin polarized density functional theoryを基にしており、汎関数はGGA−RPBEを用いた。各原子に対して、Effective Core Potentialsを与え、計算基底関数はDNPを用い、K点のサンプリングはΓ点のみで行った。
【0023】
(構造決定工程)
本実施例で利用した酸化チタンスラブモデル構造および物理量は以下の構造決定工程と呼ぶ4つの工程それぞれより得た。
【0024】
構造決定工程(1)
酸素分子が吸着していない、酸化チタンスラブモデルで、チタン元素、あるいはかつ酸素元素の一部が他の元素に置換されている構造を初期構造とし、構造最適化を実行し、最適化された構造を得た。
【0025】
構造決定工程(2)
構造決定工程(1)で得られた構造の触媒表面に、酸素分子を近づけた構造を初期構造とし、構造最適化を実行し、最適化された構造を得た。前記最適化された構造の全エネルギーの値を所得した。さらに吸着酸素分子のマリケン電荷を、触媒表面に安定に吸着した酸素分子が持つマイナス電荷として所得した。
【0026】
構造決定工程(3)
構造決定工程(2)で得られた構造の吸着酸素分子の酸素間距離を伸ばした構造を初期構造とし、構造最適化を実行し、酸素原子間距離とTi‐O間の距離を比較したときに、酸素原子間距離の方が長い最適化された構造を得た。
【0027】
構造決定工程(4)
構造決定工程(2)および(3)で得られたそれぞれの最適化された構造を酸素分子解離反応の始状態および終状態とし、LST/QST法を実行し、構造決定工程(2)で得られた最適化された構造と構造決定工程(3)で得られた最適化された構造を結ぶ遷移状態の構造を得た。該遷移状態の構造の全エネルギーの値を所得した。該全エネルギーの値から構造決定工程(2)で得られた全エネルギーの値を差し引き、それを酸素結合解離の活性化障壁とした。
【0028】
(活性化障壁)
構造決定工程(4)で得られる活性化障壁は酸素解離の活性化障壁である。酸素解離の活性化障壁が低い酸素還元触媒が、すなわち、触媒能の高い酸素還元触媒である。
【0029】
実施例1:
酸素元素の一部がF、ClあるいはBr元素に置換されており、かつその酸素元素の置換割合が3.125%である酸化チタンスラブモデルにおいて、前記構造決定工程(1)から(2)を実行し吸着した酸素分子のもつマイナス電荷を得た。
【0030】
比較例1:
実施例1と同じ酸化チタンスラブモデルで、前記構造決定工程(1)から(4)を実行し酸素結合解離の活性化障壁を得た。
【0031】
実施例2:
チタン元素の一部がSb、BiあるいはOs元素に置換されており、かつ、そのチタン元素の置換割合が6.25%である酸化チタンスラブモデルにおいて、前記構造決定工程(1)から(2)を実行し吸着した酸素分子のもつマイナス電荷を得た。
【0032】
比較例2:
実施例2と同じ酸化チタンスラブモデルで、前記構造決定工程(1)から(4)を実行し酸素結合解離の活性化障壁を得た。
【0033】
実施例3:
チタン元素の一部がSb元素に置換されており、かつそのチタン元素の置換割合が1.56%である酸化チタンスラブモデルにおいて、前記構造決定工程(1)から(2)を実行し吸着した酸素分子のもつマイナス電荷を得た。
【0034】
比較例3:
実施例3と同じ酸化チタンスラブモデルで、前記構造決定工程(1)から(4)を実行し酸素結合解離の活性化障壁を得た。
【0035】
実施例4:
元素の置換がされていない酸化チタンスラブモデルにおいて、前記構造決定工程(1)から(2)を実行し吸着した酸素分子のもつマイナス電荷を得た。
【0036】
比較例4:
実施例4と同じ酸化チタンスラブモデルで、前記構造決定工程(1)から(4)を実行し酸素結合解離の活性化障壁を得た。
【0037】
各実施例及び各比較例で用いた酸化チタンスラブモデルの構造を
図1から
図4として表1に示した。
図1の酸化チタンスラブモデルでは、F、ClあるいはBr元素が触媒表面から3および4層目に4個存在し、
図2の酸化チタンスラブモデルではSb、BiあるいはOs元素が触媒表面から3および4層目に4個存在する。
図3の酸化チタンスラブモデルでは、Sb元素が触媒表面から4層目に1個存在する。
図4の酸化チタンスラブモデルでは、元素が置換されていない。
【0038】
また、各実施例及び各比較例で得られた、吸着した酸素分子がもつマイナス電荷および酸素結合解離の活性化障壁を表1にまとめた。なお、前記吸着酸素分子が持つマイナス電荷を得る時間は、前記活性化障壁の計算時間に比べ、いずれも約1/4であった。
【0039】
【表1】
【0040】
表1に示されるように、吸着した酸素分子が持つマイナス電荷と酸素結合解離の活性化障壁とを比較すると、前記マイナス電荷が大きくなると、前記活性化障壁が減少していくことが確認できる。ゆえに、活性化障壁エネルギーを計算せずとも、吸着酸素分子が持つマイナス電荷を評価することで、酸化還元触媒の触媒能を評価することができる。
【0041】
なお、活性化障壁の計算は、構造決定工程(3)で示した通り、遷移状態の構造を決定する必要があり、この計算時間は構造決定工程(2)などで得られる最適化された構造を決定するための計算時間と比較し非常に多くの時間が必要である。それに比べ、本発明の評価法ではより短い時間で酸素還元触媒の活性能を評価することができる。