RISTニュース No．
５８（２０１５）

セシウム高選択性吸着材ゼオライトに対する
数値シミュレション研究の進展
第一原理計算及びモンテカルロ法による吸着性能評価
Rec
entSt
at
uso
fNume
r
i
c
alSi
mul
at
i
o
nSt
udi
e
sf
o
rZe
o
l
i
t
e
sas
Hi
ghl
ySe
l
e
c
t
i
veCe
s
i
um Ads
o
r
be
nt
sbyFi
r
s
t
Pr
i
nc
i
pl
e
s
Cal
c
ul
at
i
o
nandMo
nt
eCal
r
oMe
t
ho
d
日本原子力研究開発機構 システム計算科学センター
中村 博樹、奥村 雅彦、町田 昌彦
東京電力福島第一原子力発電所（１F）事故以降、放射性セシウムの優れた吸着材として、ゼオライトやフェロシアン化物等が注目され多くの研究が行われてきた。それらは実際、１Fでの汚染水処理や環境水のセシウム浄化等に用いられ、その優れた性能は広く認知されている。しかしながら、その高い性能の微視的メカニズムは十分に研究されておらず、更なる性能向上を目指す際に必要な材料学的指針は未だ十分ではない。また、実験により様々な吸着材候
補物質の性能が評価されてきたが、実験条件によって得られるデータは異なる等、物質本来の性能評価は容易でなく、計算科学による精度の高い性能評価が求められている。本稿では最新の計算科学的手法を用いることで、上記の課題に取り組んできた著者らの主な研究開発の成果を示し、今後の方向性についても議論する。
§１． はじめに日々発生する汚染水の成分に注目すると、

２
０１１年３月１１日に発生した東北地方太平洋

セシウム（Cs），ストロンチウム（Sr），トリ

沖地震に起因する東京電力福島第一原子力発

チウム（３H）が主要な放射性元素であるが、電所（１Ｆ）事故では、原子炉内にあった核Csについては、これまでに設置したCsの分

燃料や周囲の構造材料が冷却材喪失によって離除去装置によって、その濃度を十分に低減

破損・融解し、圧力容器下部やその下へと崩することが可能である。一方、Sr，３Hやその落したものと見られている。崩落後もデブリ

他の微量の放射性元素の分離除去はCs
程容と呼ばれる燃料を含む部分については、熱が易ではないため、タンクに大量の汚染水とし発生し続けるため、冷却水循環施設が設置さてやむなく貯蔵されてきた。しかし、最近でれ、事故後４年余りが経つ現在でも、炉内のは、３Hを除く核種を分離除去可能とする多核温度を抑制すべく稼働し続けている。しか種 除 去 設 備（Advanc
e
dLi
qui
dPr
o
c
es
s
i
ng

し、注入された冷却水は炉内を巡る間に高濃System， ＡＬＰＳ）により汚染水の処理が進度に汚染される他、事故後に注水された海水

められつつある。

や絶え間なく浸出し原子炉建屋に溜まる地下

汚染水の処理に当たっては、対象放射性元水等の影響を受け、様々な成分を有する汚染素に対する高性能の分離除去を可能とする吸水となる。その上、地下水は日々絶えること

着材等の存在が必須であり、その材料の開発

なく浸出するため、汚染水の量は増え続け、及び選定においては、種々の試験によって多その累積量は膨大な量となり、その処理は１Ｆ

くのデータが取得されてきた。しかし、こう

廃炉の工程に対して大きな足枷となってきた。

した一連の開発過程には多大な時間やコスト−33−

RISTニュース No．
５８（２０１５）がかかる他、優れた材料を見逃している可能

着分離した後、選択性の極めて高い材料を用性もある。例えば、試験等で得られるデータ

いて、Cs
濃度を十分に低減するといった手順

は、条件を変えると吸着性が異なってしまうが取られる等の工夫を施す必要がある。こう場合もあり、物質本来の吸着性能を反映してした材料特性に依存した吸着材の選択に際しいないこともある。そこで、開発過程の短縮ては、その性能を詳細に知る必要がある一やコストの低減を図り、且つ科学的見地から方、材料科学としてみた場合、その性能の起材料の系統的評価を効率的に実施するために源を知り、その特徴的な性能の更なる向上のは、計算科学による研究開発手段が有効と著ために、根本的な吸着形態等の基礎知見や材者らは判断し、その手法確立を目指して研究料が本来有する種々の物性データを取得する

開発を進めてきた。計算科学は、現在、高精ことが重要であることが分かる。今現在、汚度な第一原理計算により吸着機構の詳細を明

染水処理については、種々の経験や物性情報

らかにできる一方、分子動力学法やモンテカが蓄積され、多核種除去施設等の稼働によるルロ法に代表される原子・分子の古典力学に効率的な汚染水処理へと結実してきたが、そ基づく手法も進展し、精度は落ちるが様々なの処理は、今後も長く続く廃炉工程において吸着材の吸着性能を推定する上で極めて有望

欠かすことができず、その基礎となる関連材

な手段が存在するため、計算科学が上記課題

料の物性研究を更に進展させる必要がある。

の解決に大きく貢献できる可能性が高い。本また、こうした科学的知見の蓄積は、１Ｆ事稿では、現時点で利用可能な上記の先端的手故処理そして廃炉に向けた汚染水処理に係る法により、どれほど吸着材開発に貢献可能

喫緊の問題だけでなく、核燃料廃棄物等の再か、或いは新たな材料開発が可能かについ処理やその保管等とも関連し、今後、日本がて、著者らの研究成果を示し議論する。

直面せざるを得ない更に大きな問題において

１Ｆ事故後の汚染水処理に当たり、Csの分も、必須な研究開発課題であり、その研究の離除去はその放射性レベルの高さから、最も意義は極めて高いと言える。

優先度が高く、当初よりゼオライト、チタン

以上、研究開発の現状や意義について述べ

酸シリケート、フェロシアン化合物等が活用てきたが、実際にどのような指針を立て研究され、それらの高い分離除去性能が報告されを進める必要があるかについて考えよう。こてきた［１−４］。これらの材料の特徴をまとれまでの汚染水処理やその他の処理においてめると、イオン交換によりCsを選択的に吸着分かった重要な事実は、まず、如何に優れたするが、その選択性が高い上、吸着速度も速性能を持つ材料を開発できるかである一方、く、水溶液中からの分離除去性能が極めて高

それを長期に渡り安定にプラントとして稼働い。しかし、廃棄物処理という観点からは、させ、その後、保管管理できるかであること

その後に残る廃棄物の安定性等も鑑みて、総が分かる。材料が持つ性能発現の微視的メカ合的見地から材料を選択する必要があり、一ニズムや結晶構造が持つ機械的特性等の基礎

長一短であることが分かる。事実、選択性、物性を取得することは、最近、計算機の発達速さと言う観点では、フェロシアン化合物＞

により進歩してきた第一原理計算が極めて有チタン酸シリケート＞ゼオライトの順となっ効であるが、その計算例と成果については本ているが［１］、取扱いの容易さや廃棄物とし稿の前半（§
２）にて示す。但し、第一原理計てみた場合の安定性においては、ゼオライト

算で扱える原子・分子数は限られており、基が優れているため、高濃度の汚染水に対して本的に絶対零度での電子構造を詳細に計算すは、まずゼオライトにて大半のセシウムを吸る手段であることに注意する必要がある。つ−34−

RISTニュース No．
５８（２０１５）

まり、吸着材を工学的に評価する際に重要な題を議論する。§
４では、これらの成果を基

データとなるイオン交換等温線、固液分配係に、放射性物質の分離除去という課題におい数、吸脱着速度定数等の材料が示す熱力学デて、計算科学の先導による材料開発という計ータを直接評価することは困難であり、これ

算材料科学者の夢がどこまで実現可能か、そらの量を直接、計算科学にて計算できることして、今後どのような研究開発が必要かにつが次のステップとなる。もし、計算科学手法

いて議論する。

により、これらの量を既知の材料にて再現可能とするだけでなく、未知の材料においても§２．第一原理計算とゼオライト（モルデナ

予言可能となれば、計算科学は実験を代替すイト）の高いセシウム吸着選択性メカることができると言えよう。本稿の後半（§
３）

ニズム

では、この取組みの一つとして、モンテカル

１Ｆ事故直後の福島原発において発生したロ法を用いたイオン交換等温線の計算例を報

重要な課題は、大量に発生した汚染水から半告する。

減期が３０年程度のCs
−１
３
７やSr
−９
０を如何に効

本稿の構成は以下の通りである。§
２にて、率的に取り除くかであった。従来、燃料の再第一原理計算によるゼオライトを題材とし処理工程では、上記の核分裂生成物の除去のて、セシウムに対する高選択性のメカニズム

ため、ゼオライトを吸着剤として使うことがについて報告した後、今後の材料開発につい提案されていたが、スリーマイル島での原発ての方向性や計算科学を用いる上での課題等

事故や今回の福島原発でも、実際に汚染水処について議論する。§
３では、上で記したよ理プラントにゼオライトが投入され、その吸うに、材料の工学的評価において指標となる着剤としての性能は十分に認知されたと言え物性データを直接、原子・分子レベルのシる。その一方、より高性能な吸着剤探索のた

ミュレーションから得る試みとしてのモンテ

め、１F事故後、原子力学会の有志によって、カルロ法についての成果を報告し、今後の課

種々のゼオライトや他候補材料に対し、様々

図１

ゼオライトの一種「モルデナイト」の結晶構造（左）とCsを吸着する８員環の構造（右）。
黄色球はケイ素、青球はアルミニウム、赤球は酸素。

−35−

RISTニュース No．
５８（２０１５）な吸着剤の試験結果が報告されている［５］。

細孔に対して、セシウムとナトリウムの選択

本§
２では、実際に、Csに対して高い選択性性を調べるための高精度第一原理計算によるを有することで知られるゼオライトの一種で解析結果を報告する［８］。具体的な計算方法

あるモルデナイトに着目し、その高選択性のを示そう。まず、図１に示してあるように、

メカニズムを明らかにする。

モルデナイト中の細孔の一つである８員環に

一般にゼオライトによるセシウムやストロ注目する。この細孔は３Å程度の半径を有しンチウムの吸着は、イオン交換によって起こており、セシウムと酸素のイオン半径（１．
６９ることが知られている［６］。ゼオライトといÅ と１．
４０Å）を考慮すると、セシウムの吸着

う物質は、基本的にはケイ素原子と酸素原子に最適と考えられる。この細孔の内部に対からなる骨格を有し、この骨格中にいくつもし、９×９×９のメッシュを作成し、各メッの細孔が存在する（図１参照）。この骨格を形

シュ点上にてセシウムまたはナトリウムをお成するケイ素原子（＋４価）の一部が、アルいたときのエネルギーを第一原理計算で評価

ミニウム原子（＋３価）に置き換わることにした。この際、アルミニウムとケイ素の比をよって骨格が負に帯電する（図１、２参照）。

考慮するため、８員環を構成するケイ素原子

この際、電気的中性を保つために、ナトリウ（０〜４個）をアルミニウム原子（０〜４個）

ム（＋１価）などの陽イオンが骨格内の細孔に置換した場合の計算も行なった。アルミニに留まるが、このゼオライト中の陽イオンがウムが１つ増えた場合、陽イオンも１つ増やセシウムイオンを含む水溶液に触れると、ナ

さないと結晶全体の電気的な中性が壊れてしトリウムイオンとのイオン交換によりセシウ

まうため、骨格内の電子数を調整することでム吸着が起こる。

これに対応した。この方法は仮想結晶近似とゼオライトはその骨格の構造だけで１００以呼ばれる手法である。尚、第一原理計算は、上の種類があり［７］、その種類によって吸着精度の高さや優れた計算性能を有することでしやすい陽イオンが変わる。例えば、セシウ

定評の高いVASPコード［９］を用いて行なっ

ムを選択的に吸着するゼオライトとしては、た。モルデナイト中のセシウム吸着サイトとモルデナイトやチャバサイトが知られている

しては、８員環以外にも複数あり得るが、最が、それらに共通する陽イオンの選択性を制もセシウム吸着に適したものとして８員環に御する要因として次の事が知られている。

のみ着目した。現実の吸着性能等を評価する（１）陽イオンを吸着する細孔の大きさ、（２）場合は、すべての吸着サイトでの交換エネル

アルミニウムとケイ素の比。要因（１）に関ギーを考慮する必要があるが、主要な吸着メ

しては細孔の大きさと陽イオンの大きさが関カニズムのみを明確にすることを目的とした係していることが容易に想像できる一方、要ため、複数のサイトの存在については触れな

因（２）が選択性に影響を与える機構は明確い。一方、次の§３では、これらの効果を含めになっていない。ここではこれらの要因が、た、総合的な評価を可能とするモンテカルロ

どのように陽イオン選択性と関係しているか

法を用いた結果を報告する。

を第一原理計算により明らかにし、よりセシ

第一原理計算により得られた結果を図２にウム選択性の高いゼオライトの開発に向け、示す。モルデナイトの８員環内にセシウム及微視的レベルの理解から得られる材料開発のびナトリウムイオンを配置したときのエネル

方向性を議論する。以下では、セシウム高選

ギーの分布を図２の左及び右側に各々示し択性を有するゼオライトの一種「モルデナイ

た。これらの図では、エネルギー最小値でのト」の最も重要な細孔の一つに注目し、このエネルギーを基準として、濃淡でそのエネル−36−

RISTニュース No．
５８（２０１５）

図２ モルデナイト８員環におけるCsとNaのエネルギー等高線
マップ。８員環を構成するSi
原子を置換したAl
原子の数ごとに表している。NaとCsに対するポテンシャルの概略図も示す（下図）。球の色と元素の関係は図１参照。

−37−

RISTニュース No．
５８（２０１５）

ギー値（濃：小、淡：大）を表わしている。

ポテンシャルが深くなり安定化されることを実際の計算は３次元的に行なわれているが、意味しており、Csの場合その条件が満たされ

ここではエネルギーが最小値をとる点を含む

る一方、Naの場合は、全体から受ける相互作

２次元的な面のみを示した（図２の８員環の用ポテンシャルが分散されてしまい、力の重

穴を塞ぐような２次元面上でナトリウムイオ畳効果が半減してしまう。これが、Naの方

ンもセシウムイオンも安定化することに着が、イオン半径が小さく電荷密度が大きいに目）。図２から分かるように、セシウムは８も関わらず、クーロン相互作用によるエネル員環のほぼ中央の１点にエネルギーの極小域

ギー利得を得られず、Csと交換した方が安定が位置するのに対し、ナトリウムの場合は、となり、イオン交換が進行する主たる理由でアルミニウム原子数に依存するが、極小域がある。上記のように、今回考慮した８員環で分裂している場合が多い。このように極小域

はセシウムの方がナトリウムより安定して吸が分裂するかどうかが吸着における安定性に着されることが微視的に確かめられた。もち大きな影響を与えると考えられる。図２の下ろん、より小さい細孔ではイオン半径の小さ部にある概念図にあるように、８員環を構成なナトリウムの方がより安定に吸着することする酸素原子からのポテンシャルの極小点がが想像できる。こうして、今回用いた８員環一つに重なるほうが、極小点が分裂した場合の半径と同じ３Åの半径を持つ細孔がセシウに比べて、より吸着エネルギーが低くなる。

ム吸着に適していることが分かった。

これは８員環全体から受ける相互作用ポテン

次にナトリウムとセシウムのイオン交換エ

シャルの極小点が一点に集中するため、よりネルギーを実際に評価した。イオン交換エネ

図３ モルデナイトの８員環におけるNaとCsの交換エネルギーとNa−O、Cs−Oの結合距離。
Al原子の置換数の関数として表している。

−38−

RISTニュース No．
５８（２０１５）ルギー EEXはいる）。また、アルミニウム原子数が１と３の場合は、セシウム選択性が幾分悪くなっているが、この交換エネルギーの変化にはナトで表わされる。ここでEMOR （陽イオン）はリウムと酸素原子の距離との相関が見られ最安定点に陽イオンを含んだ場合のモルデナ

る。つまり、セシウム選択性が悪くなるの

イトのエネルギーであり、ESは陽イオンがは、（１）
アルミニウム原子が増えた場合と、固体として存在する場合のエネルギーであ（２）陽イオンと酸素原子の距離が近づいた場る。この交換エネルギーが低い方がセシウム合である。これらは以下のような議論から理に対する選択性が高く、逆に交換エネルギー解できる。まず、ケイ素がアルミニウムに置が正になるとナトリウムを選択的に吸着する

換されることによってゼオライトの骨格の負ようになる。但し、ここではエネルギーの基の電荷が増大する。これにより酸素原子の実準として固体のエネルギー ESを用いたが、効負電荷も増大する。酸素原子と陽イオンの

本来は水和イオンのエネルギーを基準とした間には静電気力による引力が働いているが、方が現実的である。但し、水和エネルギーを

酸素の負電荷の増大に伴い、この引力もまた採用したとしても、上記の定性的議論は変わ大きくなる。その一方、ナトリウムとセシウらない。図３に交換エネルギーをアルミニウ

ムの場合、イオン半径の小さいナトリウムのム原子の数の関数として示した。全体的な傾方がより酸素に近づくことができる。それ向としてはアルミニウム数の増加に伴い、セ故、静電気力の増大の影響を強く受けるのもシウムの選択性が劣化している様子が分かる

ナトリウムの方であり（図４参照）、結果とし（交換エネルギーがAl
原子数と共に増大しててアルミニウム原子数が増えることによって

図４

Si
原子をAl
原子に置換したときの近接O原子がNa,Csに及ぼすポテンシャルの変化を示した概略図。Oに隣接するSiをAlに置換したことによってOの電荷が変化し、NaとCsに対するOのポテンシャルが細線（赤）から太線（青）に変化する。Naの方がこの変化の影響を大きく受
ける。

−39−

RISTニュース No．
５８（２０１５）得られるエネルギー利得はセシウムよりナト

た、条件（３）は本研究にて得られた新たな提リウムの方が大きくなる。これが、アルミニ案である。条件（２）
だけでは、アルミニウム

ウムが増えた場合にセシウム選択性が比較的原子が単純に減った方がよいということにな

悪化する原因である。つまり、相対的にナトるが、その場合、吸着できる陽イオンの数もリウムの細孔内の安定性が増すことでセシウ減ることになり、吸着剤としての能力を下げ

ムの選択性が悪くなるのである。その一方、ることにも繋がる。その一方、条件（３）を満陽イオンと酸素の距離が近いときにセシウム

たせば、アルミニウムの数を減らす必要もな

選択性が悪くなる原因も同様の説明で理解でく、高いセシウム選択性を保持しつつ、吸着きる。図４においてアルミニウム原子数が１

容量を更に大きくできる可能性がある。今後と３の時に、陽イオンと酸素の距離が縮まっは実際のゼオライトの結晶構造に対して、上たのは、アルミニウム原子が８員環の片側に記の条件を満たす物質の探索を行ない、計算偏ったため、負電荷もアルミニウムリッチな科学によるセシウム高選択性ゼオライトの開側に偏って、よりゼオライト骨格に陽イオン発に繋げて行きたい。また、これらの結果かが近づいたためである。このようなアルミニ

ら、他の放射性元素イオンの選択的吸着につウムの偏りがなければ、陽イオン・酸素間距

いても同様の議論から推察可能であることが離を縮小せずとも、安定点に陽イオンを位置分かる。

させることが可能となり、その方が、セシウ§３．モンテカルロ法と種々のゼオライトのム選択性が高まるはずである。アルミニウム

吸着等温線

数が２と４の場合が、正に偏りのない状況に対応（図３）しており、セシウム選択性を高前§２において、ゼオライトの吸着サイト

めるにはアルミニウムを、吸着する細孔の周となる細孔に着目し、吸着における選択性の囲に均等に配置することが必要となる。このメカニズムを第一原理計算を用いて明らかにような詳細な制御が可能かどうかは今の所、したが、§３ではモンテカルロ法を採用する

十分な知見を著者らは有していないが、ドーことで、実際に測定され、工学的評価指標と

プ原子位置のランダム性が十分に保障されるなるイオン交換等温線を如何にしてシミュレことが必要条件であり、合成に十分な時間をーションで再現し、既存の材料の再評価と新かける必要性を示していると言えよう。以しい材料開発に繋げていけるかについて議論

上、ここまでの議論から、セシウムに対するしよう。

吸着の選択性を高めるため、ゼオライトに対

ゼオライトは一般にイオン交換を起こしてし次の３つの条件を提案することができる。

吸着分離を実現するが、§２にて議論したよ
うに、吸着選択性に際し、二つの重要な因子（１）半径３Å程度の細孔を多く有すること。

があることが分かる。一つ目は、ミクロなス（２）アルミニウムとケイ素の比が比較的小さ

ケールでの孔構造であり、二つ目は、Si/Al比である。SiがAlに置換されると、それだけ電いこと。
（３）アルミニウム原子がセシウムを吸着する

荷がマイナスとなるため、それを補償するた

細孔の周囲に均等に分布すること。

めにアルカリ元素が孔内に留まることになる。つまり、Si/Al
比が小さいと保持するア

条件（１）
（２）
はこれまでにも経験的に知られルカリ元素の原子数は多くなる一方、Si/Alていたが、今回の第一原理計算による解析比が大きいとその逆となる。§
２で対象としで、そのメカニズムはクリアになった。またモルデナイトは、Si/Al
比が比較的大きい−40−

RISTニュース No．
５８（２０１５）

図５

ゼオライトのフレーム構造。
（a）ゼオライトY、（b）ゼオライトA、（c）モルデナイト。赤、緑の球はそれぞれ酸素原子、ケイ素原子を表す。ゼオライトYとゼオライトAにはフレーム構
造の模式図（黄色の図）を併記した。

一方、ゼオライトＹやゼオライトＡでは小さ上げる。これまでに、モンテカルロ法は様々いことから、経験的にSi/Al比が大きい方がな目的で多用されてきたが、ゼオライトのイCsの吸着及び選択性が良いという知見がオン交換等温線を計算により初めて求めたのあった。しかし、こうした知見はあくまで経

は、Je
f
f
r
o
ye
tal
．
［１
，
０］である。彼らはゼオ

験的なものであり、本§
３では、この量を変化

ライトＹを対象とし、それについては実験とさせることで、実際に吸着性能の指標である

良く一致した結果を報告している。彼らの採イオン交換等温線が実際、どのように変化す用したモンテカルロ法は、基本的には、一般

るかについて議論したい。

的なモンテカルロ法と共通であるが、準グラ

まず、本§３で対象とする複数のゼオライ

ンドカノニカルモンテカルロ法と呼ばれるもトの結晶構造を示す（図５）
［７］
。ゼオライ

のであり、数種類の競合イオンの存在下で効トＹ，モルデナイト、ゼオライトＡ、そして率良くイオン交換等温線を計算することがで

ZK−４（これは、ゼオライトＡのSi
/
Al
比を変

きる。この手法では、スワップ移動と呼ばれ化させたものである）を計算対象として取りる手法を通常のモンテカルロ法に加えてい−41−

RISTニュース No．
５８（２０１５）る。この方法は、複数種存在する陽イオンにことが検証できた。こうして、著者らはモル対し、その総数は変えずに、適宜、２種類のデナイトを始めとする他のゼオライトにも、イオン種の位置を交換する。そして、交換に

上記手法を適用し、特にCsのイオン交換等温より生じたゼオライト内での陽イオンの新配線に着目し、ゼオライト間の差やSi/Al比を置については、陽イオンの種毎の化学ポテン

変えた場合のイオン交換等温線の変化を調べ

シャルを基に吸着（交換）前後のエネルギーることとした［１２］。

を計算し、モンテカルロ法により棄却したり

モンテカルロ法を用いた計算結果を示そ採用したりする。更にゼオライト以外の水分う。結果を示す前に、まず計算するイオン交子の運動（平行移動と回転）や陽イオンの運

換等温線の意味について簡単に説明しよう。

動を考慮する。以上が用いた手法の概要だ

図６は、典型的なイオン交換等温線の例を示

が、Jef
f
r
oyetal．
，はゼオライトＹにおいてしている。横軸は対象としている液体中に含Rb、Cs、KとNaのイオン交換におけるイオ

まれる全陽イオンの中での対象とする陽イオ

ン交換等温線を再現することに成功していンの割合（陽イオン全体、例えばNaイオンとる。著者らは彼らの文献［１０］に従い、提案

Cs
イオンの和に対してのCs
イオンの割合）でされたモンテカルロ手法を独立に組み込んであり、縦軸は対象とするイオン（本稿ではCs

コードを作成し並列化も行った。検証のたイオン）がゼオライトに含まれる割合である。

め、彼らの手法が採用されている商用コード

等温線とは、常に一定温度の下、溶液中での

Gi
bbs
［１
１］を用いて、著者らの作成したコ

Cs
イオンの割合が増すと、その一部はゼオラードと文献で対象としていたゼオライトＹに

イトに吸着するが、その際、Cs
イオンの割合適用したところ、両者共に文献で得られたイが増す程、吸着量も増加するため、右上がり

オン交換等温線を再現できたことから、著者の線を描くことになる。その際、ゼオライト

らのオリジナルコードも正常に動作しているがCs
イオンに対し選択性が高く、少しのCs

図６ イオン交換等温曲線の概略図。縦軸がゼオライト中の陽イオンXの割合を表し、横軸は溶液中での割合を表す。赤線のようになるとき、吸着材はXへの高選択性を示したことになる。

−42−

RISTニュース No．
５８（２０１５）

イオンでもゼオライトに吸着されてしまう場

一方、ゼオライトAやZK-４の立ち上がりは

合は、赤線のような線が得られる。そして、

悪く、緑線と似たパターンを示しており、選Cs
イオンの割合が増えると、その多くは吸着択性が高くないことが見て取れる。また、ゼされるが、割合が十分大きくなると、今度はオライトYは、モルデナイト程の選択性はな

ゼオライトの方が飽和してしまい、右肩上がいが、ある程度の選択吸着力を有しているこりの線が飽和するような線が得られる。そのとが分かる。これらの結果において、実線や一方、選択性が低い場合は、Cs
イオンの割合点線らは、シミュレーションにより得られたが十分な値となって初めて吸着し始めるた

イオン交換等温線である一方、ドット（上向め、緑線のような線が得られるのである。

き△（灰色）と下向き▽（青色））は実験によモンテカルロ法により得られた各ゼオライ

り得られたデータであり、シミュレーション

トが示すイオン交換等温線（NaとXのイオンが実験結果（種々のゼオライトに対して）を交換）を図７に示す。図中、上に並べた図はかなりの精度で再現できていることが分か

各ゼオライトのCs
に対するイオン交換等温る。実験結果は、データの分散が大きいが、線であり、下に並べた図は、カリウム（K）

異なるゼオライトに対する傾向の差異は、に対するイオン交換等温線である（CsもKもはっきりと見て取ることができる一方、シ

競 合 イ オ ン はNa）。選 択 性 が 高 く １F事 故ミュレーションも、その傾向の違いをはっき

後、様々な場面にて実際に用いられているモりと再現していることが分かる。シミュレー

ルデナイトの場合、Cs
イオンの割合が僅かでションによってこのように十分な再現力がああっても吸着するため、先ほど説明した赤線

ることが分かったのは、最近の研究の成果でのパターンを示していることが分かる。その

あり、Je
f
f
r
o
yらによる先進的な手法開発と著

図７

各ゼオライトのNaに対するXイオンのイオン交換等温曲線。
（a）ゼオライトY（X=Cs）、（b））、（d）モルデナイト（X=K）、（e）ゼオライトゼオライトY（X=K）、（c）モルデナイト（X=−Cs）、（h）ZK−４（X=K）のイオン交A（X=Cs）、（f）ゼオライトA（X=K）、（g）ZK−４（X＝Cs（g）と（h）で は、Si/Al
比を５換 に 対 応。各 図 中、線 が 計 算 結 果、点 が 実 験 結 果［１３−１６］。
/
３，３，７と変えている。

−43−

RISTニュース No．
５８（２０１５）

者らによる系統的計算により初めて確認されの違いとして現れることは、ミクロレベルのた。次に詳細な計算結果の分析になるが、図物性の差異が材料機能を左右するという関係

７（g）にて、ゼオライトZK−４のSi/
Al
比を性を如実に示しており、ミクロレベルでの制変化した場合の吸着等温線の変化を示した。

御が機能材料では重要であることが分かる。

Si/Al
比が小さい場合は、吸着性能は低いが、Si/Al比を上昇させると、その性能はある程§４．計算科学による吸着材開発の将来度向上するが飽和してしまうことが見て取れ§
２，３において、第一原理計算によるゼオラる。この原因を説明するに当たり、まず、何

イトの陽イオンの吸着選択性メカニズムとモ

故Si
/
Al
比を上昇させると性能が向上するか

ンテカルロ法を用いて得られるイオン交換等に つ い て 説 明 し よ う。Si/Al比 が 小 さ い 場

温線の計算結果（実験結果の再現と予測）に合、Alが含まれる割合が大きく、ゼオライト

ついて述べてきた。まず、Csイオンに対する

骨格は負の電荷を十分持つため、陽イオンを高い選択性を有するモルデナイトのメカニズ

たくさん保持している。このような場合、陽ムについては、孔内サイズとイオンのサイズ

イオンをトラップするミクロの孔の周辺にてが適合することにより、Cs
イオンが感じるエ

Alが高い密度で分布するため、孔とイオンのネルギー安定点（極小点）が一点に集約する

サイズが一致することによる整合効果に対ことで、より深くなり安定化することが分し、負電荷と陽電荷との相互作用による単純

かった。この効果は、細孔がトラップされたなエネルギー利得が優勢となり、よりNaのCs
イオンを取り囲み、周囲からの相互作用が方が相互作用による利得を得てしまうためで重畳することで、より深いエネルギー極小点ある。つまり、細孔サイズによる選択性が働を作るというメカニズムである。これは、単

きつつも、電荷間の直接の作用によりそれがなる点電荷間の１対１の１次元的相互作用でマスクされてしまう。一方、Si
/
Al
比を高く

なく、高次元で多体問題として現れる効果でするとゼオライトの骨格が有する負電荷の量ある。このような絶妙な効果を調べるためには減り、より骨格による孔のサイズによる適は、第一原理計算が最も信頼性が高く、電荷合効果が顕著となり、イオンの選択性が卓越の分極や電荷の移動などの効果も含めて、細すると考えられる。しかし、Si
/
Al
比が増加

孔構造とイオン間の相互作用を詳細に調べるしても、ZK−４の孔は２種類存在し、一方はことができる。その一方、第一原理計算が有Naイオン、もう一方はCsイオンに適合して効に機能するのは、ゼオライトの一部の構造いるため、Csの選択性はいずれ飽和してしまだけを取り出し、その一部に集中した場合のう。図７（ｇ）の結果はそれを如実に示していみであり、計算の精度をある程度犠牲にしてることが分かる。尚、これらの説明の一部はも、競合するイオンを含む水溶液に浸された

基本的には§
２の第一原理計算結果の説明に

ゼオライトという現実の状況を模擬し、現実

おいて、より詳細に説明したが、モンテカルに得られる評価指標であるイオン交換等温線

ロ法でも同じ傾向が得られていることが分かを再現できるシミュレーション技術も極めてる。更にモンテカルロ法を用いることで、複有用であることが分かる。理想を言えば、第数タイプの細孔が存在する場合の吸着性能に一原理計算により得られる精度の高さをそのついても、その性能評価をトータルな性能とまま維持し、モンテカルロ法等のマクロな物

して得ることができることを指摘する。この性評価計算を実行するのが望ましいが、第一ようなミクロの差異が、イオン交換等温線の原理計算の枠組みにて系全体を計算するにような工学的評価指標においても重要な性質

は、膨大な計算量が必要であり、原子・分子−44−

RISTニュース No．
５８（２０１５）間の力を与えるポテンシャルの表式は簡単なれるイオン交換等温線等の熱力学レベルのデモデルであることが望ましい。従って、現時ータをモンテカルロ法にて再現することに成点では、メカニズム評価については第一原理

功し、異なるゼオライトが示す性質の違いを計算、実験との直接の対応や材料の工学的評再現した他、これまで経験的な知見に留まっ価については、古典的な原子・分子間力を基

ていたSi
/
Al
比に依存する吸着性能の変化をにしたモンテカルロ法等のレベルで計算する

微視的な要因と結びつけて説明することがでことが妥協点であり、両者の利点を有する手きた。これらの成果を基に、今後は他の材料法の一つであるQM/
MM法等の適用も今後の（例えばフェロシアン化物やチタンシリケー

課題と考えられる。但し、モンテカルロ法

ト、そして、グラフェンや有機化合物、生体は、稀に起こるイオン交換を即時に実行した分子等）についても吸着メカニズムや機能評

後、状態の採用や棄却を判定する等、高いポ価を進め、種々の知見を蓄積することで、今テンシャルを稀に超えて実現する状態も評価

後の材料開発に貢献したいと考えている。しできるといったメリットがあり、時間発展をかしながら、計算科学は日々進歩しており、追跡する際にぶつかる時間の壁をどのように本稿で示していない計算手法や考え方により超えるかという課題もあることが分かる。つ

抜本的な進歩があるかもしれない。また、計まり、現在の計算力の限界下にて工学的評価

算機のリソース環境においても「京」等スーが可能な計算科学を実施するに当たっては、

パーコンピュータの進展によって、これまで重要な過程を見定める先験的知見が必須であ不可能とされてきた計算が可能となり、よりり、専門家による知見の蓄積を文献情報から

現実に即した計算結果が示され、実際の材料

見出し、的確にモンテカルロ法等のシミュレ開発に貢献する成果も多々報告されているとーションに組み入れることが必須であるこという現実がある。今後は、第一原理計算によが分かる。このような研究はスーパーコン

る微視的メカニズムの探求、工学的評価に資

ピュータ「京」やその後に予定されているエするための計算手法の改良、そして、超並列

クサスケールの計算機を用いたとしても、明計算を駆使した計算力不足の打破といった側らかに必要であり、今後も、大規模計算によ面を有機的に活かし、上記の取り組みを続る直接計算手法による研究開発と共に、物質

け、微力ながら計算科学が福島復興の一助と材料に対するなるよう、そして、日本が直面している原子

センスを磨く研究も進めていく必要性がある。

力が抱える大きな課題の解決に力を尽くして行ければと考えている。

§５．まとめ謝辞

本稿では放射性物質の分離除去に関し、吸着材（ゼオライト）に対する計算科学による本研究に際し、村松健、吉田善幸、大貫敏

研究成果を記し、今後の研究の方向性につい

彦、河村雄行、山岸功、大西康夫各氏より多ても議論した。第一原理計算によるゼオライ

くのご教示をいただきました。ここに感謝致

トの吸着メカニズム解明においては、細孔全します。計算は全て原子力機構・システム計

体で陽イオンを取り囲むことで得られる相互

算科学センターのBX９００を用いて行いまし

作用の重畳効果が重要であることが明らかにた。谷正之センター長を始めとする全てのス

され、計算科学の重要性を示すことができた

タッフの支援に感謝致します。また、計算にと著者らは考えている。また、吸着材の工学

際し、（株）菱化システムの藤原敦志、森健人

的価値を評価するに当たり、実験にて測定さ両氏より技術的支援をいただきました。両氏−45−

RISTニュース No．
５８（２０１５）に感謝致します。最後に本研究の一部は、原Mac
hida: J. Phys
. Soc
. Jpn.８２,子力機構・福島技術本部による支援の下、行

０
２
３
８
０
１(
２
０
１
３)
.
［９］G.KresseandJ.Fur
t
hmül
l
e
r
:Phys
.

われました。特に、廃棄土壌減容プロジェク

Re
v.B５
４,
１
１
１
６
９(
１
９
９
６)
.

トの支援に対して、リーダーの矢板毅氏を始［１０］M. Jeffroy, A. Bout
in and A. H.

めとするそのメンバー全員に感謝致します。

Fuchs,J.Phys.Chem.B,
１
１
５,
１
５
０
５
９１
５
０
６
６(
２
０
１
１)
.

参考文献
［１］三村均，日本イオン交換学会誌 ２５，４５［１１］Me
de
AGi
bbs l
i
c
e
ns
e
d f
r
o
m: I
FP
Ene
r
gi
e
s no
uve
l
l
e
s and Labo
r
at
o
i
r
e

（２０１
４）
．
［２］三村 均，佐藤 修彰，桐島 陽，日本イオ

deChi
mi
e
Phys
i
que
,Uni
ve
r
s
i
t
e de

ン交換学会誌 ２２，９
６（２０
１１）．

de
A ve
r
s
i
o
n
Par
i
s
Sud, CNRS; Me［３］三村 均，山岸 功，日本イオン交換学会

２.
１
３ ２
０
１
３ Mat
e
r
i
al
s De
s
i
gn, I
nc
.
,
Ange
lFi
r
e
,NM,USA.

誌 ２３，６（２０１２）．

［１２］H. Nakamura, M. Okumura, M.

［４］三村 均，山岸 功，日本イオン交換学会

Machida,RSC Adv.４,５
２
７
５
７５
２
７
６
１

誌２
３，２９（２０
１２）．

(
２
０
１
４)
.

［５］日 本 原 子 力 学 会 バ ッ ク エ ン ド 部 会：
ht
t
p:
/
/
www.
nuc
e
ae
s
j
.
o
r
g/

［１
３］H.S.Sherry,J.Phys.Chem.
,１
９６
６,［６］ S.
M.Aue
r
bac
h,K.
A.Car
r
ado
,P.
K.
Dut
t
a(
Eds
.
)
: Handbo
o
ko
fZe
o
l
i
t
e

７
０,
１
１
５
８１
１
６
８.
［１
４］F.Wo
l
f
,H.F¨
ur
t
i
g and H.Kno
l
l
,

Sc
i
enc
e and Te
c
hno
l
o
gy (
Mar
c
e
l
De
kker
,Ne
w Yo
r
k,
２
０
０
３)
.

Chem.Te
c
hn.
,
１
９
７
１,
２
３,
２
７
３.
［１
５］L.L.Ames,Can.Mi
ne
r
al
.
,１
９
６
５,８,
３
２
５３
３
３.

［７］Ch.Baer
l
o
c
he
r
,L.B.Mc
Cus
ke
rand
s
o
n, At
l
as o
f Ze
o
l
i
t
e
D. H. Ol

［１６］R.M.Barrer,L.V.C.ReesandD.J.

Fr
ame
wo
r
kType
s
,El
s
e
vi
e
r
,Ams
t
e
r
-

War
d,Pr
o
c
.R.So
c
.Lo
nd.A,１
９
６
３,

dam,
２００７.

２
７
３,
１
８
０１
９
７.

［８］H.Nakamura,M.Okumur
a,andM.

−46−