(57)【特許請求の範囲】
【請求項１】
少なくとも反応Ａとその競争反応Ｘが生じる反応系で反応Ａを優先的に起こさせる触媒を探索する触媒の探索方法であって、
反応Ａの反応種が前記触媒に吸着した状態のギブスエネルギーＧ1Aと、反応Ａの生成種が前記触媒に吸着した状態のギブスエネルギーＧ2Aの差ΔＧA＝Ｇ2A−Ｇ1Aを算出し、
反応Ｘの反応種が前記触媒に吸着した状態のギブスエネルギーＧ1Xと、反応Ｘの生成種が前記触媒に吸着した状態のギブスエネルギーＧ2Xの差ΔＧX＝Ｇ2X−Ｇ1Xを算出し、
前記ΔＧXおよびΔＧAの値に基づき、複数の化合物を前記触媒の候補として選定し、前記複数の化合物のΔＧX及びΔＧAをそれぞれ算出し、前記各化合物の−（ΔＧA−ΔＧX）の値を算出し、当該値が大きい化合物を触媒として選定することを特徴とする触媒の探索方法。
【請求項２】
請求項１に記載された触媒の探索方法において、
前記反応種がそれぞれ前記触媒に吸着した時の各吸着エネルギーの和ΣＥreactと、前記生成種がそれぞれ前記触媒に吸着したときの各吸着エネルギーの和ΣＥproductと、前記反応種および生成種のエネルギー差ΔＥgasとを算出し、前記ギブスエネルギーの差ΔＧをΣＥreact−ΣＥpro−ΔＥgasとして算出することを特徴とする触媒の探索方法。
【請求項３】
請求項１ないし２のいずれかに記載された触媒の探索方法において、
反応種と生成種のエネルギー差、または前記反応種または生成種の触媒への吸着エネルギーを第一原理計算を利用して求めることを特徴とする触媒の探索方法。
【請求項４】
請求項１ないし３のいずれかに記載された触媒の探索方法において、
前記反応ＡがＣＯを還元剤としたＮＯｘ還元反応であり、前記反応Ｘが酸素によるＣＯの酸化反応であることを特徴とする触媒の探索方法。
【請求項５】
計算機よりなり、少なくとも反応Ａとその競争反応Ｘが生じる反応系で反応Ａを優先的に起こさせる触媒を探索する触媒の探索装置であって、
反応Ａの反応種が前記触媒に吸着した状態のギブスエネルギーＧ1Aと、反応Ａの生成種が前記触媒に吸着した状態のギブスエネルギーＧ2Aとの差ΔＧA＝Ｇ2A−Ｇ1Aを算出し、反応Ｘの反応種が前記触媒に吸着した状態のギブスエネルギーＧ1Xと、反応Ｘの生成種が前記触媒に吸着した状態のギブスエネルギーＧ2Xとの差ΔＧX＝Ｇ2X−Ｇ1Xを算出するΔＧ算出手段と、
前記ΔＧXおよびΔＧAの値に基づき、複数の化合物を前記触媒の候補として選定し、前記複数の化合物のΔＧX及びΔＧAをそれぞれ算出し、前記各化合物の−（ΔＧA−ΔＧX）の値を算出し、当該値が大きい化合物を触媒として触媒を選定する選択手段とを有することを特徴とする触媒の探索装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、シミュレーションにより特定の反応を選択的に起こさせる触媒を探索する方法に関する。
【背景技術】
【０００２】
従来の触媒開発法としては、研究者が経験により触媒元素を選定し、試験により評価するという流れを繰り返し、試行錯誤の上、触媒材料を探索してきていた。しかし、その際にかかる労力や材料コストは膨大である。そこで、より効率的に触媒を探索できる方法が求められてきている。例えば特開２００５−２７０８５２号公報（特許文献１）には、一の触媒の電子の分布状態として電子のエネルギーをパラメータ化し、代替となる触媒の探索を行うことについて開示している。
【０００３】
【特許文献１】特開２００５−２７０８５２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
探索の手法の可能性の一つが計算化学的手法である。しかしながら、計算化学的な手法は、材料の物性値予測等に利用されてきているが、反応の予測問題に関しては、未だ実用段階に至っていない。また、特に二種類以上の競争反応が起こりうる系では、触媒が高選択率であることが求められる。
【０００５】
二つ以上の競争反応が起こりうる系において、目的の反応のみが生じやすい材料を触媒とするが、その探索方法は経験的であり、現在まで確立された方法は存在しない。また、このような競争反応系において、試行錯誤により望む反応の選択性を向上させる高選択率の触媒を探索することは非常に労力を要する。そこで本発明の課題は、二種類以上の競争反応が進行する可能性がある系（競争反応系）において、特定の反応を選択的に起こさせ、もしくは望む反応の選択性を向上させる触媒を探索する手法を提供することにある。
【課題を解決するための手段】
【０００６】
上記課題を解決する本発明の特徴は、二つ以上の競争反応が起こる系で、求める反応を他の競争反応よりも優先的に起こさせる触媒を、求める反応及び競争反応のそれぞれの反応種／生成種が触媒に吸着した状態のギブスエネルギーを利用することにより、シミュレーション計算し探索する方法にある。本発明では、二つ以上の競争反応が起きる場合、本来、複雑な素過程を含む競争反応について、各反応の進行しやすさを、触媒表面に反応物が吸着した際のギブスエネルギーと、触媒表面に生成物が吸着した際のギブスエネルギーの差としてあらわせることに着目した。これらを比較した結果に基づき触媒材料となる可能性が高い物質を選定することができる。
【０００７】
本発明の触媒の探索方法は、求める反応Ａと、競争反応Ｘとの競争反応が起こる系において、反応Ａを優先的に生じさせる触媒を計算機によるシミュレーションで探索する方法であって、反応Ａの反応種が触媒に吸着した状態のギブスエネルギーＧ1Aと、反応Ａの生成種が触媒に吸着した状態のギブスエネルギーＧ2Aと、競争反応Ｘの反応種が触媒に吸着した状態のギブスエネルギーＧ1Xと、競争反応Ｘの生成種が触媒に吸着した状態のギブスエネルギーＧ2Xとより、Ｇ2AとＧ1Aとの差ΔＧA（ΔＧA＝Ｇ2A−Ｇ1A）、およびＧ2XとＧ1Xとの差ΔＧX（ΔＧX＝Ｇ2X−Ｇ1X）を求める工程と、ΔＧAとΔＧXとを比較し、ΔＧAがΔＧXよりも大きい材料を選定する工程とを有する。
【０００８】
また、上記のギブスエネルギーＧはそれぞれ算出してもよいが、反応にかかわる分子の吸着エネルギーと脱離エネルギーを利用すると、前述のギブスエネルギー差を比較的簡易に概算できることを見出した。
【０００９】
従って本発明は、反応に用いられる各反応種が触媒の表面に吸着した際に得られる各吸着エネルギーを算出し、前記各反応種の吸着エネルギーを合計したΣＥreactを算出し、反応で生成する各生成種が触媒の表面に吸着した際に得られる各吸着エネルギーを算出し、前記各生成種の吸着エネルギーを合計したΣＥproductを算出し、反応種と生成種のエネルギー差（反応エネルギー）ΔＥgasを算出しΣＥreactよりΣＥproとΔＥgasとを減算して前記ΔＧを求める（ΔＧ＝ΣＥreact−ΣＥpro−ΔＥgas）ことが好ましい。なお、前記吸着エネルギーや反応エネルギーは、例えば第一原理計算を利用して求めることができる。
【００１０】
上記手法によって、ΔＧAがΔＧXよりも大きい触媒材料を選定することにより、反応Ａが進行しやすい触媒の探求が容易になる。また、上記シミュレーションは計算機で行われるため、事前にシステムを構築し、計算機と組み合わせた触媒探索装置としておくことで容易に触媒の探求が可能である。
【発明の効果】
【００１１】
本発明によれば、高い反応の選択率を有する触媒を提供できる。また、選択率の高い触媒を開発する際にかかるコスト（労力，時間，実験材料費用など）の低減ができる。
【発明を実施するための最良の形態】
【００１２】
次式（１）で表される反応と、次式（２）で表される二つの競争反応において、式（１）における反応のみを選択的に生じやすい触媒候補を探索する場合について述べる。
【００１３】
Ａ＋Ｂ → Ｃ＋Ｄ                                                ・・・（１）
Ａ＋Ｅ → Ｆ＋Ｇ                                                ・・・（２）
一方、式（１）（２）における反応が触媒の存在下で生じる場合について式（３）と（４）に示す。例えば、式（１）では、反応物ＡとＢが触媒表面に吸着し、触媒表面上で反応して生成物ＣとＤとなる。反応式（３）と（４）では、触媒表面に反応種が吸着した状態を記号（ａ）で表している。反応式（１）の反応物ＡとＢが触媒表面に吸着した状態をＡ（ａ）＋Ｂ（ａ）とし、生成物であるＣとＤが表面に吸着した状態をＣ（ａ）＋Ｄ（ａ）で表す。また、反応式（２）の反応物ＡとＥが表面に吸着した状態をＡ（ａ）＋Ｅ（ａ）、生成物であるＦとＧが表面に吸着した状態をＦ（ａ）＋Ｇ（ａ）で表す。
【００１４】
Ａ（ａ）＋Ｂ（ａ） → Ｃ（ａ）＋Ｄ（ａ）                        ・・・（３）
Ａ（ａ）＋Ｅ（ａ） → Ｆ（ａ）＋Ｇ（ａ）                        ・・・（４）
式（３）と式（４）それぞれの触媒表面反応の熱力学的安定性を示すため、反応の前後のギブスエネルギー差ΔＧ3，ΔＧ4を検討する。式（３）（４）の各状態のギブスエネルギーをＧ（Ａ，Ｂ），Ｇ（Ｃ，Ｄ），Ｇ（Ａ，Ｅ），Ｇ（Ｆ，Ｇ）で表す。各反応における触媒表面に反応物が吸着した際のギブスエネルギーと触媒表面に生成物が吸着した際のギブスエネルギーの差は式（５）と式（６）にて計算される。
【００１５】
ΔＧ3＝Ｇ（Ｃ，Ｄ）−Ｇ（Ａ，Ｂ）                               ・・・（５）
ΔＧ4＝Ｇ（Ｆ，Ｇ）−Ｇ（Ａ，Ｅ）                               ・・・（６）
このギブスエネルギー差ΔＧの値が小さい程、反応は進行しやすい。従って、二つの競争反応（式（３）及び式（４））において、求める反応である式（３）に対応するΔＧ（式（５））が小さく、阻害反応（競争反応）のΔＧが大きい程、触媒材料として優れていると判断する。従って、反応ギブスエネルギーを概算し、触媒を探索できる。
【００１６】
図１は、触媒が存在するときの反応物Ａと反応物Ｂから生成物Ｃと生成物Ｄが生じる反応のエネルギーダイアグラムを示す図である。この反応物Ａと反応物Ｂが触媒表面に吸着した状態と生成物Ｃと生成物Ｄが吸着した状態のエネルギー差でΔＧを見積もり、触媒の反応進行しやすさとした。
【００１７】
なお、各反応における吸着状態のギブスエネルギーは、第一原理計算による吸着状態のエネルギーとして計算できる。各反応における反応物と生成物の吸着エネルギーを計算し、表面上の反応物と表面上の生成物のギブスエネルギー差を概算することで、各競争反応において大小関係を比較し、触媒材料の選定を行う。
【００１８】
例えば反応式（３）の場合は、図１から分かるように、式（７）によりΔＧ3を算出する。
【００１９】
ΔＧ3＝Ｅads(Ｃ)＋Ｅads(Ｄ)−（Ｅads(Ａ)−Ｅads(Ｂ)）−ΔＥ1    ・・・（７）
ΔＥ1は反応式（１）の反応エネルギーである。また、Ｅads（Ａ），Ｅads（Ｂ），Ｅads（Ｃ），Ｅads（Ｄ）は、それぞれＡ，Ｂ，Ｃ，Ｄの触媒表面への吸着エネルギーである。
【実施例】
【００２０】
以下実施例では、窒素酸化物を窒素ガスに還元する触媒の探索について説明する。本実施例は、前記反応ＡがＣＯを還元剤としたＮＯｘ還元反応であり、前記競争反応Ｘが酸素によるＣＯの酸化反応である場合を想定している。
【００２１】
ボイラや自動車等による窒素酸化物は大気汚染物質のひとつであり、規制が強化されてきている。現在、ボイラから排出される窒素酸化物は、アンモニアを還元剤とした脱硝触媒法により窒素に還元されている。
【００２２】
しかし、アンモニアの調達，貯蓄等に必要な運転コストを抑えるために、同じ排ガス中に存在する一酸化炭素を還元剤とした脱硝触媒の開発がされている。酸素存在下での脱硝触媒による窒素酸化物の還元反応（式（８））では、窒素酸化物の還元のほか、酸素により還元剤である一酸化炭素の酸化反応（式（９））も起こってしまう。
【００２３】
ＮＯ（ａ）＋ＣＯ（ａ） → １／２Ｎ2（ａ）＋ＣＯ2（ａ）          ・・・（８）
ＣＯ（ａ）＋１／２Ｏ2（ａ） → ＣＯ2（ａ）                      ・・・（９）
酸素存在下で窒素酸化物を一酸化炭素によって還元するためには、一酸化炭素と酸素の反応よりも、一酸化炭素による窒素酸化物の還元反応が優先的に起きるような触媒を使用する必要がある。このような酸素に影響されない触媒が検討されているものの、実用可能な触媒はまだ見つかっていない。このように二種類以上の競争反応が起こりうる系は、触媒が高選択率であることが求められる。
【００２４】
そこで、種々の金属を触媒材料として想定し、選択率の高い触媒を探索した。
【００２５】
まず、触媒表面における一酸化炭素と窒素酸化物の反応（式（８））の反応ギブスエネルギーΔＧNOと、触媒表面における一酸化炭素と酸素の反応（式（９））の反応ギブスエネルギーΔＧO2を算出した。簡易に反応ギブスエネルギーΔＧNO，ΔＧO2を見積もるために、以下の手段を利用した。まず式（１０），式（１１）に示すガスの反応を想定した。（ｇ）はガス状態であることを示す。
【００２６】
ＮＯ（ｇ）＋ＣＯ（ｇ） → １／２Ｎ2（ｇ）＋ＣＯ2（ｇ）        ・・・（１０）
ＣＯ（ｇ）＋１／２Ｏ2（ｇ） → ＣＯ2（ｇ）                    ・・・（１１）
次に、式（１０），式（１１）に示すガスの反応エネルギーをそれぞれΔＥN2(g)，ΔＥCO2(g)とし、これをバンド計算により求めた。さらに、各ガスの触媒表面への吸着エネルギーＥadsをバンド計算により求めた。ＮＯの触媒表面への吸着エネルギーをＥads（ＮＯ）、ＣＯの吸着エネルギーをＥads（ＣＯ）、Ｎ2の吸着エネルギーをＥads（Ｎ2）、ＣＯ2の吸着エネルギーをＥads（ＣＯ2）、Ｏ2の吸着エネルギーをＥads（Ｏ2）と示す。
【００２７】
上記の反応エネルギー，吸着エネルギーを用いて、反応ギブスエネルギーΔＧNO（式（１２））ΔＧO2（式（１３））を算出した。
【００２８】
ΔＧNO＝１／２Ｅads（Ｎ2）＋Ｅads（ＣＯ2）−（Ｅads（ＮＯ）＋Ｅads（ＣＯ））
−ΔＥN2(g)                                           ・・・（１２）
ΔＧO2＝Ｅads（ＣＯ2）−（Ｅads（ＣＯ）＋１／２Ｅads（Ｏ2））−ΔＥCO2(g)
・・・（１３）
ΔＧNO，ΔＧO2を比較し、反応ギブスエネルギーは値が低い方が熱力学的に有利に進行するので、ΔＧO2＞ΔＧNOが成り立てば、その触媒は有効であると判断する。特に、−（ΔＧNO−ΔＧO2）が正で大きい程、触媒性能が高いことが予想される。
【００２９】
この方法を利用して、白金，金，イリジウム，ロジウム、及びパラジウムについて、ＮＯ，ＣＯ，Ｎ2，ＣＯ2、及びＯ2の吸着エネルギーをバンド計算法により計算した。バンド計算には密度汎関数法の局所密度近似を利用した。結果を図２に記載する。
【００３０】
これらの結果と先に計算したΔＥN2(g)，ΔＥCO2(g)から、−（ΔＧNO−ΔＧO2）を見積もった。図３に結果を記載する。各触媒表面における、ＮＯとＣＯの反応ギブスエネルギーを概算した値、ＣＯとＯ2の反応ギブスエネルギーを概算した値、それらの差にマイナスをかけたもので二つの競争反応の進行しやすさを予測する値を示している。
【００３１】
続いて、実際のＮＯの浄化率を測定し、計算結果の検証を行った。アルミナに、白金，金，イリジウム，ロジウム、及びパラジウムをそれぞれ０.１５mol％担持し、ＮＯ浄化触媒とした。各触媒について、ＮＯｘ浄化率を評価した。入り口ガス組成は、ＮＯが１００ppm、ＣＯが２０００ppm、Ｏ2が６％であり、残りはＮ2とした。ＳＶは１５０００（ｈ-1）とした。反応温度は３５０℃とした。ＮＯｘの浄化率試験結果を図４に示す。実際のＮＯｘの浄化率の試験結果と、先ほど計算から見積もった各触媒の−（ΔＧNO−ΔＧO2）とを比較した結果を図５に示す。−（ΔＧNO−ΔＧO2）がゼロより大きいものが高い性能を示しており、計算結果と実測値の間には相関があった。図５より、−（ΔＧNO−ΔＧO2）がＮＯの選択還元反応の起きやすさ（ＮＯ還元率）を表す値であることが明らかである。
【００３２】
上述のように、触媒材料の選定前に、複数ある競争反応の表面における反応ギブスエネルギーを第一原理計算により求め、比較することで、複数ある触媒材料化合物の候補より、優位な結果を得られる可能性が高いものを選定することができる。
【図面の簡単な説明】
【００３３】
【図１】触媒が存在するときの反応物Ａと反応物Ｂから生成物Ｃと生成物Ｄが生じる反応のエネルギーダイアグラムを示す図である。
【図２】各種触媒，各種反応化合物のそれぞれの組み合わせの吸着エネルギーの計算結果を示す図である。
【図３】各種触媒表面の反応ギブスエネルギーと、二つの競争反応の進行しやすさを比較する図である。
【図４】各触媒の３５０℃におけるＮＯ還元率の測定結果を示す。
【図５】−（ΔＧNO−ΔＧO2）の値と、各触媒のＮＯ還元率の関係を表すグラフである。
【図１】
0004861950.tif 000002
【図２】
0004861950.tif 000003
【図３】
0004861950.tif 000004
【図４】
0004861950.tif 000005
【図５】
0004861950.tif 000006