COSMO法のソースを表示

[[File:4-Nitrobenzoic-acid-elpot-3D-vdW.png|thumb|COSMO法により計算した4-ニトロ安息香酸の電荷密度分布]]
'''COSMO法'''(COSMOほう、'''CO'''nductor-like '''S'''creening '''MO'''delの略)<ref name=":0" /><ref name=":1" />とは、[[分子]]と[[溶媒]]の[[静電相互作用]]を決定するための計算法である。[[計算化学]]において[[溶媒効果]]をモデル化するためによく用いられる。

COSMO法は溶媒を[[誘電率]]''&epsilon;''を持つ連続体として扱うため、連続体溶媒モデルの一種として分類される。他の連続体溶媒モデルと同様に、COSMO法では溶媒を分子キャビティーの外側から溶質を取り囲む誘電連続体として近似する。ほとんどの場合において、分子キャビティーは[[ファンデルワールス半径]]のおよそ1.2倍の半径を持つ原子を中心とした球の組み合わせにより構成される。実際の計算においては、分子キャビティーの表面はセグメント（例：六角形・五角形・三角形）により近似される。

他の連続体溶媒モデルとは異なり、COSMO法では溶質の極性に起因する連続体の分極電荷を{{訳語疑問点範囲|date=2018-09-24|scaled-conductor近似}}により計算する。もし仮に溶媒が理想的な導体であるならば、キャビティー表面の[[電位|電位差]]は0になるはずである。量子化学計算などから溶質分子の[[電荷|電荷分布]]が分かっている場合、表面セグメント上の電荷''q''*を決定することができる。 有限の誘電率を持つ溶媒に対しては、この電荷は次式の係数''&fnof;''(''&epsilon;'')によってより小さくなるようにスケールされる。

:<math>q = f(\varepsilon) q^*.</math>

係数''&fnof;''(''&epsilon;'')は次式で近似される。

:<math>f(\varepsilon)=\frac{\varepsilon-1}{\varepsilon+x},</math>

ここで''x''の値は中性分子で0.5、イオンでは0.0になるように選ばれる。詳細は原著論文を参照のこと<ref name=":1">{{Cite book|title=From Quantum Chemistry to Fluid Phase Thermodynamics and Drug Design|last=Klamt|first=Andreas|publisher=Elsevier|year=2005|isbn=9780444519948|location=Boston, MA, USA|pages=}}</ref>。

このようにして決定された溶媒の電荷''q''と既知である溶質分子の電荷分布から、溶媒と溶質分子間の相互作用エネルギーが計算できる。

COSMO法は電荷密度が計算可能な全ての理論化学的手法（例えば[[半経験的分子軌道法]]・[[ハートリー-フォック法]]・[[密度汎関数法]]など）において利用可能である<ref name=":0">{{cite journal |last=A.|first=Klamt|last2=G.|first2=Schüürmann|date=1993|title=COSMO: a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient|url=http://pubs.rsc.org/en/content/articlepdf/1993/p2/p29930000799|journal=J. Chem. Soc.|publisher=Perkin Trans.2|volume=2|issue=5|pages=799–805|doi=10.1039/P29930000799|access-date=25 October 2015}}</ref>。

== 異なるバージョンとその実装==
COSMO法は[[Amsterdam Density Functional|ADF]]、 [[GAMESS]]<ref>{{cite web|url=https://www.msg.chem.iastate.edu/gamess/GAMESS_Manual/intro.pdf|format=pdf|title=GAMESS User's Guide|publisher=Iowa State University|accessdate=2019-02-11}}</ref>、[[GAUSSIAN]]、[[MOPAC]]、{{仮リンク|NWChem|en|NWChem}}、[[TURBOMOLE]]、{{仮リンク|Q-Chem|en|Q-Chem}}といった数多くの[[量子化学および固体物理計算ソフトの一覧|量子化学ソフトウェア]]に実装されている。[[分極連続体モデル]](PCM)のCOSMO法版もまた開発されている。

実装により、分子キャビティーの構成方法、計算に用いられる半径の値、キャビティー表面を表現するセグメントの形状、スケーリング関数''&fnof;''(''&epsilon;'')中の''x''の値などが異なることがある{{要出典|date=2019-02}}。

== 他手法との比較 ==
電荷密度の{{ill2|多重極展開|en|multipole expansion}}による手法の適用範囲が小分子や亜球形・楕円形の分子に制限されている一方で、COSMO法には大きな分子や不揃いな形をした分子にも適用できるというメリットがある。

誘電率について明確な境界条件を課す分極連続体モデル(PCM)とは対照的にCOSMO法では近似的なスケーリング関数f(''&epsilon;'')を用いる。スケーリング関数が近似であるにもかかわらず、COSMO法は関連する誤差を減少させ、いわゆる{{訳語疑問点範囲|date=2018-09-25|outlying charge}}をより正確に表現することができることが判明している。 COSMO法と積分方程式表式化PCM(IEFPCM)を比較したところ<ref name=":2">{{Cite journal|last=Klamt|first=A.|last2=Moya|first2=C.|last3=Palomar|first3=J.|date=2015|title=A Comprehensive Comparison of the IEFPCM and SS(V)PE Continuum Solvation Methods with the COSMO Approach|doi=10.1021/acs.jctc.5b00601|journal=Journal of Chemical Theory and Computation|volume=11 |issue=9|pages=4220–4225|via=}}</ref>実験データのばらつきと比べて計算手法の違いに起因する計算値の違いは小さいことが明らかにされている。複数の連続体溶媒モデルの間の詳細な違いよりも、溶媒を連続体として扱っているがために無視される[[水素結合]]や再配向といった効果に起因する誤差の方が実験データの再現により大きな影響を及ぼす。

== 脚注 ==
{{reflist}}

== 関連項目 ==
* [[分極連続体モデル]]
* {{ill2|COSMO-RS法|en|COSMO-RS}}

[[Category:計算化学]]