特開 2024-152042 希釈剤選択支援装置、希釈剤選択支援方法、およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2024152042

(43)【公開日】2024-10-25

(54)【発明の名称】希釈剤選択支援装置、希釈剤選択支援方法、およびプログラム

(51)【国際特許分類】

   B01D 11/04 20060101AFI20241018BHJP

   C22B 3/02 20060101ALI20241018BHJP

   C22B 3/26 20060101ALI20241018BHJP

【ＦＩ】

B01D11/04 B

C22B3/02

C22B3/26

【審査請求】未請求

【請求項の数】7

【出願形態】ＯＬ

(21)【出願番号】P 2023065952

(22)【出願日】2023-04-13

(71)【出願人】

【識別番号】000183303

【氏名又は名称】住友金属鉱山株式会社

(74)【代理人】

【識別番号】100107766

【弁理士】

【氏名又は名称】伊東 忠重

(74)【代理人】

【識別番号】100070150

【弁理士】

【氏名又は名称】伊東 忠彦

(72)【発明者】

【氏名】西原 泰孝

【テーマコード（参考）】

4D056

4K001

【Ｆターム（参考）】

4D056AB03

4D056AC01

4D056AC02

4D056AC03

4D056CA01

4D056CA39

4D056CA40

4D056DA10

4K001AA07

4K001AA16

4K001DB26

(57)【要約】      （修正有）

【課題】溶媒抽出法における希釈剤の選択を支援可能な希釈剤選択支援装置を提供する。
【解決手段】希釈剤選択支援装置は、抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得する分子情報取得部と、前記複数の希釈剤分子のそれぞれを使用した場合における有機相と水相とを含む立体モデルを生成する立体モデル生成部と、前記立体モデルにおいて分子動力学計算を行い、前記有機相と前記水相との界面における圧力テンソル、および抽出対象の金属イオンと前記金属イオンの周囲に存在する前記抽出剤分子の座標を算出する分子動力学計算部と、前記複数の希釈剤分子のそれぞれについて算出された前記圧力テンソルに基づいて、前記界面における界面張力を算出する界面張力算出部と、前記複数の希釈剤分子のそれぞれについて算出された前記抽出剤分子の前記座標に基づいて、前記金属イオンを中心とする動径分布関数を算出する動径分布関数算出部と、を備える。
【選択図】図１

【特許請求の範囲】

【請求項1】

溶媒抽出法における希釈剤の選択を支援するための希釈剤選択支援装置であって、
抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得する分子情報取得部と、
前記複数の希釈剤分子のそれぞれを使用した場合における有機相と水相とを含む立体モデルを生成する立体モデル生成部と、
前記立体モデルにおいて分子動力学計算を行い、前記有機相と前記水相との界面における圧力テンソル、および抽出対象の金属イオンと前記金属イオンの周囲に存在する前記抽出剤分子の座標を算出する分子動力学計算部と、
前記複数の希釈剤分子のそれぞれについて算出された前記圧力テンソルに基づいて、前記界面における界面張力を算出する界面張力算出部と、
前記複数の希釈剤分子のそれぞれについて算出された前記抽出剤分子の前記座標に基づいて、前記金属イオンを中心とする動径分布関数を算出する動径分布関数算出部と、を備える、希釈剤選択支援装置。

【請求項2】

算出された前記界面張力、および前記動径分布関数に基づいて希釈剤を選択する希釈剤選択部をさらに備える、請求項１に記載の希釈剤選択支援装置。

【請求項3】

前記立体モデル生成部は、構造最適化計算によって前記希釈剤分子がランダムに配置された前記立体モデルを生成する、請求項１または請求項２に記載の希釈剤選択支援装置。

【請求項4】

前記分子動力学計算部は、昇温計算、ＮＶＴ計算およびＮＰＴ計算の少なくともいずれかを行って、前記圧力テンソルを算出する、請求項１または請求項２に記載の希釈剤選択支援装置。

【請求項5】

前記複数の希釈剤分子のそれぞれについての希釈剤を示す名称と、算出された前記界面張力を示す情報と、前記動径分布関数を示す情報とを関連付けて表示する表示部をさらに備える、請求項１または請求項２に記載の希釈剤選択支援装置。

【請求項6】

溶媒抽出法における希釈剤の選択を支援するためのコンピュータが実行する希釈剤選択支援方法であって、
抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得するステップと、
前記複数の希釈剤分子のそれぞれを使用した場合における有機相と水相とを含む立体モデルを生成するステップと、
前記立体モデルにおいて分子動力学計算を行って、前記有機相と前記水相の界面における圧力テンソル、および抽出対象の金属イオンと前記金属イオンの周囲に存在する前記抽出剤分子の座標を算出するステップと、
前記複数の希釈剤分子のそれぞれについて算出された前記圧力テンソルに基づいて、前記界面における界面張力を算出するステップと、
前記複数の希釈剤分子のそれぞれについて算出された前記座標に基づいて、前記金属イオンを中心とする動径分布関数を算出するステップと、を備える、希釈剤選択支援方法。

【請求項7】

溶媒抽出法における希釈剤の選択を支援するためのコンピュータに、
抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得するステップと、
前記複数の希釈剤分子のそれぞれを使用した場合における有機相と水相とを含む立体モデルを生成するステップと、
前記立体モデルにおいて分子動力学計算を行って、前記有機相と前記水相の界面における圧力テンソル、および抽出対象の金属イオンと前記金属イオンの周囲に存在する前記抽出剤分子の座標を算出するステップと、
前記複数の希釈剤分子のそれぞれについて算出された前記圧力テンソルに基づいて、前記界面における界面張力を算出するステップと、
前記複数の希釈剤分子のそれぞれについて算出された前記座標に基づいて、前記金属イオンを中心とする動径分布関数を算出するステップと、を実行させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、希釈剤選択支援装置、希釈剤選択支援方法、およびプログラムに関する。

【背景技術】

【0002】

水溶液に含まれる金属イオンを選択的に抽出する方法の一つに溶媒抽出法（ＳＸ：Ｓｏｌｖｅｎｔ Ｅｘｔｒａｃｔｉｏｎ）がある。溶媒抽出法では互いに混じり合わない２相が用いられるが、有機相（有機溶媒）と水相（水溶液）の組み合わせが広く用いられている。水相中の金属イオンは、有機相に添加された抽出剤と２相界面付近で結合し、有機相中へ移動する。そのため、界面の安定性（相分離性）は金属イオンの抽出平衡や抽出速度に大きな影響を与える。

【0003】

相分離性は有機溶媒（希釈剤とも呼ばれる）や抽出剤の組み合わせで変化し、希釈剤の比重や界面張力、粘度が大きく影響する。また、系の温度とも関連する。特に、界面張力が大きいときは相分離性がよい。

【0004】

そこで、界面張力に基づいて、溶媒抽出に用いる希釈剤を選択する方法が研究されている。例えば、特許文献１には、相分離した状態の超薄膜を測定することからなる遠心液膜法によって界面張力を測定する装置が開示されている。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特許第３７５１４３７号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

しかし、上述した従来の技術を用いて、実験的に界面張力を測定して希釈剤を選択することは、候補となる希釈剤が多くなると時間とコストがかかってしまうという問題がある。また、相分離性以外の要素も加味して、より最適な希釈剤の選択が求められるようになっている。

【0007】

本発明は、上記事情に鑑みてなされたものであり、溶媒抽出法における希釈剤の選択を支援可能な希釈剤選択支援装置を提供することを目的としている。

【課題を解決するための手段】

【0008】

上記目的を達成するため、本発明の一態様に係る希釈剤選択支援装置は、
溶媒抽出法における希釈剤の選択を支援するための希釈剤選択支援装置であって、
抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得する分子情報取得部と、
前記複数の希釈剤分子のそれぞれを使用した場合における有機相と水相とを含む立体モデルを生成する立体モデル生成部と、
前記立体モデルにおいて分子動力学計算を行い、前記有機相と前記水相との界面における圧力テンソル、および抽出対象の金属イオンと前記金属イオンの周囲に存在する前記抽出剤分子の座標を算出する分子動力学計算部と、
前記複数の希釈剤分子のそれぞれについて算出された前記圧力テンソルに基づいて、前記界面における界面張力を算出する界面張力算出部と、
前記複数の希釈剤分子のそれぞれについて算出された前記抽出剤分子の前記座標に基づいて、前記金属イオンを中心とする動径分布関数を算出する動径分布関数算出部と、を備える。

【発明の効果】

【0009】

溶媒抽出法における希釈剤の選択を支援可能な希釈剤選択支援装置を提供することができる。

【図面の簡単な説明】

【0010】

【図1】希釈剤選択支援装置の機能構成の一例を示す図である。

【図2】希釈剤選択支援処理の流れの一例を示すフローチャートである。

【図3】希釈剤選択支援装置のハードウェア構成の一例を示す図である。

【図4】立体モデルの一例を示す図である。

【図5】抽出剤の分子モデルの一例を示す図である。

【図6】界面張力算出結果画面の一例を示す図である。

【図7】界面張力算出結果画面の一例を示す図である。

【発明を実施するための形態】

【0011】

以下に、図面を参照して、本発明の実施の形態（本実施の形態）について説明する。なお、本発明はこれらの例示に限定されるものではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内での全ての変更が含まれることが意図される。
［希釈剤選択支援装置］
（希釈剤選択支援装置の機能構成）
図１は、希釈剤選択支援装置の機能構成の一例を示す図である。

【0012】

本実施の形態に係る希釈剤選択支援装置１０は、溶媒抽出法における希釈剤の選択を支援するための装置である。具体的には、希釈剤選択支援装置１０は、分子情報取得部１１と、立体モデル生成部１２と、分子動力学計算部１３と、界面張力算出部１４と、動径分布関数算出部１５と、を備える。希釈剤選択支援装置１０は、必要に応じて、表示部１６や、希釈剤選択部１７を備えることもできる。

【0013】

分子情報取得部１１は、抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得する。分子情報は、分子の化学構造を示す情報であって、分子を構成する原子の数と種類、原子同士の関係を示すデータ等を含む。原子同士の関係を示すデータは、原子の大きさ、結合する原子と結合の角度等を含むデータである。

【0014】

立体モデル生成部１２は、複数の希釈剤分子のそれぞれを使用した場合における有機相と水相を持つ系の立体モデルを生成する。立体モデルとは、各分子の配置を示すデータである。有機相には希釈剤分子と抽出剤分子が含まれ、水相には水分子（必要に応じてカウンターイオン）が含まれる。具体的には、立体モデル生成部１２は、候補となる複数の希釈剤分子から選択した、各希釈剤分子について、分子情報取得部１１が取得した分子情報に基づいて、有機相と水相を持つ系における各分子の配置を示すデータを、立体モデルとして生成する。ここで、立体モデル生成部１２は、各分子が空間内に分散してランダムに配置される立体モデルを生成する。

【0015】

分子動力学計算部１３は、候補となる希釈剤分子である複数の希釈剤分子ごとに生成された立体モデルにおいて分子動力学計算を行う。分子動力学計算部１３は、有機相と水相との界面における圧力テンソル、および抽出対象の金属イオンと、該金属イオンの周囲に存在する抽出剤分子の座標を算出する。

【0016】

界面張力算出部１４は、候補となる複数の希釈剤分子のそれぞれについて算出された圧力テンソルに基づいて、２相の界面における界面張力を算出する。

【0017】

動径分布関数算出部１５は、候補となる複数の希釈剤分子のそれぞれについて算出された、抽出対象の金属イオンの周囲に存在する抽出剤分子の座標に基づいて、金属イオンを中心とする動径分布関数を算出する。

【0018】

表示部１６は、界面張力算出部１４や、動径分布関数算出部１５での算出結果を表示できる。表示部１６は、例えば候補である複数の希釈剤分子のそれぞれについての希釈剤を示す名称と、それぞれの希釈剤分子を使用した立体モデルを用いて算出された界面張力を示す情報と、動径分布関数を示す情報と、を関連付けて表示できる。これによって、選択を支援する情報をユーザに提示することができる。

【0019】

表示部１６は、例えばディスプレイ等の画面に上記情報を表示できる。なお、表示部１６は、他の装置、例えばユーザの操作する端末等とし、該端末等に上記情報を表示しても良い。

【0020】

希釈剤選択部１７は、算出された界面張力、および動径分布関数に基づいて、希釈剤を選択できる。希釈剤選択支援装置１０が希釈剤選択部１７を有することで、界面張力算出部１４や、動径分布関数算出部１５での算出結果に基づいて、希釈剤を自動的に選択できる。特に希釈剤選択支援装置１０が希釈剤選択部１７を有することで、例えば希釈剤の選択におけるユーザ間のばらつきを抑え、均一な基準による選択も可能となる。

【0021】

界面張力が大きいと、相分離性に優れる。

【0022】

また、金属イオンと負電荷をもつ原子の動径分布関数では、金属イオンに結合している原子の位置を示す第一ピーク位置が金属イオンに近い方が、すなわち数値が小さい方が、有機相中での金属錯体の安定性が高くなる。第二ピーク位置等で表される金属イオンに結合している原子以外の原子の位置が、金属イオンから離れている方が、すなわち数値が大きい方が、有機相中での金属錯体の安定性が高くなる。

【0023】

このため、希釈剤選択部１７は、複数の希釈剤分子から、例えば界面張力が大きく、かつ動径分布関数の第一ピーク位置の値が小さくなる希釈剤分子を選択することができる。また、必要に応じて、動径分布関数の第二ピーク位置の値についても加味し、第二ピーク位置の値が大きくなる希釈剤分子を選択することもできる。

【0024】

なお、界面張力や、動径分布関数に関して、予め基準値を設定しておき、希釈剤選択部１７は例えば該基準値に近い値を取る希釈剤や、該基準値から一定の範囲内にある値を取る希釈剤を選択しても良い。

【0025】

なお、希釈剤選択部１７は、候補となる複数の希釈剤分子から最適な希釈剤分子を１つ選択しても良く、複数の好適な希釈剤分子を選択しても良い。さらに、表示部１６は、選択された希釈剤分子や、該希釈剤分子に対応する希釈剤を示す情報を表示してもよい。

【0026】

（希釈剤選択支援装置の動作）
次に、希釈剤選択支援装置１０の動作について、図面を参照して説明する。図２には、希釈剤選択支援処理の流れの一例を示すフローチャート２０を示している。

【0027】

ユーザの操作を受けて、分子情報取得部１１は、抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得する（ステップＳ１１）。

【0028】

次に、立体モデル生成部１２は、候補となる複数の希釈剤分子のうちの一つを選択して、選択された希釈剤分子の分子情報に基づいて、有機相と水相を持つ系の立体モデルを生成する（ステップＳ１２）。ここで、立体モデル生成部１２は、分子力学（ＭＭ：Ｍｏｌｅｃｕｌａｒ Ｍｅｃｈａｎｉｃｓ）計算や分子動力学（ＭＤ：Ｍｏｌｅｃｕｌａｒ Ｄｙｎａｍｉｃｓ）計算、量子化学（ＱＭ：Ｑｕａｎｔｕｍ Ｍｅｃｈａｎｉｃｓ）計算等による構造最適化計算を行ってもよい。立体モデル生成部１２は、上記構造最適化計算によって、希釈剤分子がランダムに配置された立体モデルを生成することもできる。これによって、立体モデルの作成において、すべての分子の配置の指定をする必要がないため、簡易な操作によって立体モデルの生成を実現することができる。

【0029】

次に、分子動力学計算部１３は、分子動力学計算によって、界面における圧力テンソル、抽出剤分子の座標を算出する（ステップＳ１３）。

【0030】

圧力テンソルは、希釈剤分子が存在する空間にあらかじめ規定されたＸＹＺ軸によってテンソル（Ｐｘｘ，Ｐｙｙ，Ｐｚｚ）と表される。

【0031】

分子動力学計算部１３は、例えば昇温計算、ＮＶＴ計算およびＮＰＴ計算の少なくともいずれかを行って、圧力テンソルを算出できる。これによって、大量の候補に対する計算を試みることが可能となり、コストと手間を軽減できる。

【0032】

具体的には例えば、分子動力学計算部１３は、生成された立体モデルにおいて、室温まで徐々に温度を上昇させた状態をシミュレートする昇温計算を実行できる。次に、分子動力学計算部１３は、原子数、体積、温度がそれぞれ規定の値に保持されるようにシミュレートするＮＶＴ計算を行うことができる。続いて、分子動力学計算部１３は、原子数、圧力、温度がそれぞれ規定の値に保持されるようにシミュレートするＮＰＴ計算を行い、ＮＰＴ計算における最後の０．１ｎｓのデータに基づいて界面における圧力テンソルを算出できる。

【0033】

分子動力学計算部１３は、抽出剤分子の座標、具体的には、抽出対象の金属イオンと、該金属イオンの周囲に存在する抽出剤分子の座標を算出できる。分子動力学計算を行うことで、分子動力学計算に供した立体モデルに含まれる各原子、分子の座標の経時変化が得られる。このため、分子動力学計算の結果から、抽出剤分子の座標を求められる。

【0034】

次に、界面張力算出部１４は、圧力テンソルに基づいて、界面張力を算出する（ステップＳ１４）。具体的には、界面張力算出部１４は、下記の式（１）によって界面張力σを算出する。

【0035】

σ＝（Ｌ／２）×｛Ｐｚｚ－（１／２）×（Ｐｘｘ＋Ｐｙｙ）｝・・・（１）
ここで、（Ｐｘｘ，Ｐｙｙ，Ｐｚｚ）は圧力テンソル、Ｌは界面の法線方向のセルサイズである。

【0036】

また、動径分布関数算出部１５は、抽出剤分子の座標に基づいて、動径分布関数（ＲＤＦ）を算出する（ステップＳ１５）。

【0037】

動径分布関数ｇ（ｒ）はある粒子からの距離に他の粒子が存在する確率を表す。

【0038】

ある粒子からの距離ｒと距離ｒ＋ｄｒとの間の球殻内に存在する粒子の数をｎ（ｒ）としたとき、球殻内の粒子密度はｎ（ｒ）／（４πｒ^２ｄｒ）と書ける。このため、該粒子密度を平均密度ρで割ったものが動径分布関数ｇ（ｒ）となり、下記の式（２）により表される。

【0039】

ｇ（ｒ）＝＜ｎ（ｒ）＞／（４πｒ^２ｄｒρ）・・・（２）
式（２）中の＜Ａ＞は、Ａの時間平均、粒子平均を表わす。

【0040】

次に、希釈剤選択支援装置１０は、候補となる複数の希釈剤分子について、全ての希釈剤を計算したか否かを判定する（ステップＳ１６）。希釈剤選択支援装置１０は、候補となる複数の希釈剤分子のうち、計算していない希釈剤があると判定すると（ステップＳ１６：ＮＯ）、ステップＳ１２に戻り、候補となる複数の希釈剤分子から、まだ計算していない希釈剤分子を１つ選択して、処理を実行する。

【0041】

希釈剤選択支援装置１０が、候補となる複数の希釈剤分子について、全ての希釈剤を計算したと判定すると（ステップＳ１６：ＹＥＳ）、該希釈剤選択支援装置１０はフローを終えることもできるが、以下のステップＳ１７以降を実施することもできる。

【0042】

具体的には例えば表示部１６が、界面張力、動径分布関数の算出結果を表示することができる（ステップＳ１７）。

【0043】

また、希釈剤選択部１７が希釈剤分子を選択して、表示部１６が選択された希釈剤分子に対応した希釈剤を示す情報を表示することもできる（ステップＳ１８）。

【0044】

希釈剤選択部１７は、１つの希釈剤を選択してもよいし、複数の希釈剤を選択してもよく、推奨する希釈剤の順位を基準に基づいて決定してもよい。

【0045】

なお、希釈剤選択支援装置１０は、希釈剤選択部１７を備えていなくてもよく、ステップＳ１７およびステップＳ１８の処理を実行しなくてもよい。その場合、ユーザは、表示された界面張力や動径分布関数の算出結果に基づいて、希釈剤を選択すればよい。

【0046】

次に、希釈剤選択支援装置１０のハードウェア構成について説明する。図３は、希釈剤選択支援装置のハードウェア構成の一例を示す図である。

【0047】

希釈剤選択支援装置１０は、コンピュータによって構成され、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１０１、主記憶装置１０２、補助記憶装置１０３、入力装置１０４、表示装置１０５、通信インターフェース装置１０６、ドライブ装置１０７を備える。これらの各装置は、バスで接続されている。

【0048】

ＣＰＵ１０１は、希釈剤選択支援装置１０の動作を制御する主制御部であり、主記憶装置１０２に格納されたプログラムを読みだして実行することで、後述する各種の機能を実現する。

【0049】

主記憶装置１０２は、希釈剤選択支援装置１０の起動時に補助記憶装置１０３からプログラムを読み出して格納する。補助記憶装置１０３は、インストールされたプログラムを格納すると共に、後述する各種機能に必要なファイル、データ等を格納する。

【0050】

入力装置１０４は、各種の情報の入力を行うための装置であり、例えばキーボードやポインティングデバイス等により実現される。表示装置１０５は、各種の情報の表示を行うためものであり、例えばディスプレイ等により実現される。通信インターフェース装置１０６は、ＬＡＮカード等を含み、他の装置等との接続の為に用いられる。
［希釈剤選択支援方法］
本実施形態の希釈剤選択支援方法について説明する。本実施形態の希釈剤選択支援方法は、例えば既述の希釈剤選択支援装置を用いて実施できる。このため、既に説明した事項は説明を省略する。

【0051】

本実施形態の希釈剤選択支援方法は、溶媒抽出法における希釈剤の選択を支援するためのコンピュータが実行する方法であり、図２を用いて説明したように、以下の各ステップを有することができる。

【0052】

抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得するステップ。

【0053】

複数の希釈剤分子のそれぞれを使用した場合における有機相と水相とを含む立体モデルを生成するステップ。

【0054】

立体モデルにおいて分子動力学計算を行って、有機相と水相の界面における圧力テンソル、および抽出対象の金属イオンと金属イオンの周囲に存在する前記抽出剤分子の座標を算出するステップ。

【0055】

複数の希釈剤分子のそれぞれについて算出された圧力テンソルに基づいて、界面における界面張力を算出するステップ。

【0056】

複数の希釈剤分子のそれぞれについて算出された座標に基づいて、金属イオンを中心とする動径分布関数を算出するステップ。

【0057】

本実施形態の希釈剤選択支援方法は、必要に応じて、図２に示したステップＳ１７や、ステップＳ１８を有することもできる。

【0058】

各ステップについては、希釈剤選択支援装置において説明したため、説明を省略する。
［プログラム］
本実施形態のプログラムについて説明する。

【0059】

本実施形態のプログラムは、溶媒抽出法における希釈剤の選択を支援するためのコンピュータに、以下の各ステップを実行させるためのプログラムである。本実施形態のプログラムは、コンピュータを希釈剤選択支援装置で、図１を用いて説明した各部として機能させることができる。

【0060】

抽出剤分子および候補となる複数の希釈剤分子の分子情報を取得するステップ。

【0061】

複数の希釈剤分子のそれぞれを使用した場合における有機相と水相とを含む立体モデルを生成するステップ。

【0062】

立体モデルにおいて分子動力学計算を行って、有機相と水相の界面における圧力テンソル、および抽出対象の金属イオンと前記金属イオンの周囲に存在する前記抽出剤分子の座標を算出するステップ。

【0063】

複数の希釈剤分子のそれぞれについて算出された圧力テンソルに基づいて、前記界面における界面張力を算出するステップ。

【0064】

複数の希釈剤分子のそれぞれについて算出された座標に基づいて、金属イオンを中心とする動径分布関数を算出するステップ。

【0065】

本実施形態のプログラムは、必要に応じて、図２に示したステップＳ１７や、ステップＳ１８をコンピュータに実行させることもできる。

【0066】

各ステップについては、希釈剤選択支援装置において説明したため、説明を省略する。

【0067】

本実施形態に係るプログラムは、希釈剤選択支援装置１０を制御する各種プログラムの少なくとも一部とすることができる。プログラムは、例えば記憶媒体１０８（図３を参照）の配布やネットワークからのダウンロード等によって提供できる。プログラムを記録した記憶媒体１０８は、ＣＤ－ＲＯＭ、フレキシブルディスク、光磁気ディスク等の様に情報を光学的、電気的或いは磁気的に記録する記憶媒体、ＲＯＭ、フラッシュメモリ等の様に情報を電気的に記録する半導体メモリ等、様々なタイプの記憶媒体を用いることができる。

【0068】

また、プログラムは、プログラムを記録した記憶媒体１０８がドライブ装置１０７にセットされると、記憶媒体１０８からドライブ装置１０７を介して補助記憶装置１０３にインストールされる。ネットワークからダウンロードされたプログラムは、通信インターフェース装置１０６を介して補助記憶装置１０３にインストールできる。

【実施例0069】

以下に、実施例によって本発明をさらに詳細に説明するが、本発明は、これらの実施例によってなんら限定されるものではない。

【0070】

本実施の形態に係る希釈剤選択支援装置１０の実施例について、説明する。

【0071】

抽出剤分子をＬＩＸ６３、候補となる複数の希釈剤分子を（ａ）ヘキサン（Ｃ_６Ｈ_１４）、（ｂ）シクロヘキサン（Ｃ_６Ｈ_１２）、（ｃ）ベンゼン（Ｃ_６Ｈ_６）として、実施した結果を以下に示す。

【0072】

図２に示した希釈剤選択支援処理のフローチャート２０に沿って実施した。

【0073】

まず、ステップＳ１１において、分子情報取得部１１が、候補となる複数の希釈剤分子、および抽出剤分子の分子情報を取得した。

【0074】

次いで、ステップＳ１２において、立体モデル生成部１２が、有機相と水相を持つ系の立体モデルを生成した。

【0075】

図４は、立体モデルの一例を示す図である。立体モデル生成部１２によって生成される立体モデルは、例えば、有機相４１と水相４２を含むボックス４０がＸＹ方向に複数個並ぶように配置されたものとしてもよい。また、ボックス４０のＸ方向の長さＬｘ、およびＹ方向の長さＬｙは、例えばそれぞれ５０オングストロームであり、ボックス４０のＺ方向の長さＬｚは、例えば１００オングストロームである。

【0076】

立体モデル生成部１２は、例えば、候補となる希釈剤分子を約６００分子含むボックス中に約２５分子の抽出剤分子（ＬＩＸ６３）が配置された有機相４１と、水約４３００分子を含む水相４２と、の２相を、ボックス４０内に生成する。

【0077】

図５は、抽出剤の分子モデル５０の一例を示す図である。図５に示す分子モデル５０は、抽出剤分子ＬＩＸ６３の分子モデルとして、分子情報取得部１１が取得した分子情報（各原子の配置や結合の関係を示す情報）を図示したものである。

【0078】

ステップＳ１３において、分子動力学計算部１３は、立体モデル生成部１２が生成した立体モデルにおいて分子動力学計算を行い、圧力テンソル、抽出剤分子の座標を算出した。

【0079】

そして、ステップＳ１４において、界面張力算出部１４は、得られた圧力テンソルに基づいて、各立体モデルでの界面張力を算出した。また、ステップＳ１５において、動径分布関数算出部１５は、得られた抽出剤分子の座標に基づいて動径分布関数を算出した。

【0080】

候補となる全ての希釈剤分子について、上記各ステップを実施した後、ステップＳ１７において、表示部１６が界面張力、動径分布関数の算出結果を表示した。

【0081】

図６は抽出対象の金属イオンがＣｏ^２＋の場合の、図７は抽出対象の金属イオンがＭｎ^２＋の場合の、界面張力、動径分布関数についての算出結果画面の一例を示す図である。図６に示した算出結果画面６０、図７に示した算出結果画面７０は、希釈剤選択支援処理のステップＳ１７において表示部１６が表示する画面の一例である。算出結果画面６０、算出結果画面７０には、各希釈剤の名称と、各希釈剤について算出された界面張力を示す情報と、動径分布関数を示す情報とが関連付けられて表示されている。

【0082】

また、表示部１６は、希釈剤選択支援処理のステップＳ１８において、希釈剤選択部１７によって選択された希釈剤を示す情報を表示することもできる。

【0083】

例えば、算出結果画面６０によれば、動径分布関数について、第一ピーク位置と第二ピーク位置とは、希釈剤によらず、大きな違いが見られないことから、界面張力が最も大きい、ヘキサンを希釈剤として選択できる。表示部１６は、算出結果画面６０において選択されたヘキサンを、他と異なる文字色で表示してもよい。

【0084】

例えば、算出結果画面７０によれば、動径分布関数について、第一ピーク位置は、希釈剤によらず、大きな違いが見られないことから、界面張力、および動径分布関数の第二ピーク位置の値が最も大きい、ヘキサンを希釈剤として選択できる。表示部１６は、算出結果画面７０において選択されたヘキサンを、他と異なる文字色で表示してもよい。

【0085】

以上に説明した本実施の形態に係る希釈剤選択支援装置１０によれば、分子動力学計算によって界面張力や、動径分布関数を算出し、相分離性、有機相中での金属錯体の安定性の指標とすることで、希釈剤の効率的な選択を支援することができる。

【0086】

以上、本実施の形態に基づき本発明の説明を行ってきたが、上記実施の形態に示した要件に本発明が限定されるものではない。これらの点に関しては、本発明の主旨をそこなわない範囲で変更することができ、その応用形態に応じて適切に定めることができる。

【符号の説明】

【0087】

１０ 希釈剤選択支援装置
１１ 分子情報取得部
１２ 立体モデル生成部
１３ 分子動力学計算部
１４ 界面張力算出部
１５ 動径分布関数算出部
１６ 表示部
１７ 希釈剤選択部
１０１ ＣＰＵ
１０２ 主記憶装置
１０３ 補助記憶装置
１０４ 入力装置
１０５ 表示装置
１０６ 通信インターフェース装置
１０７ ドライブ装置
１０８ 記憶媒体
Ｓ１１～Ｓ１８ ステップ
４０ ボックス
４１ 有機相
４２ 水相
Ｌｘ Ｘ方向の長さ
Ｌｙ Ｙ方向の長さ
Ｌｚ Ｚ方向の長さ
５０ 分子モデル
６０、７０ 算出結果画面

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】