|
局所上昇法 (きょくしょじょうしょうほう、) は、計算化学や計算物理学において主に分子シミュレーション(分子動力学法 (MD) およびモンテカルロ法 (MC) など)に用いられる手法である。1994年に分子動力学法において配座空間上における探索を加速するために Huber, Torda, van Gunsteren により開発され〔、分子動力学シミュレーション用ソフトウェア 96 以降で利用可能である。この手法は、配座フラッディング法〔と共に、分子シミュレーションに履歴依存性を導入した初めての例である。Engkvist–Karlström法〔、適応バイアス力法〔、ワン・ランダウ法、メタダイナミクス法、適応バイアス分子動力学法〔、適応反応座標力法〔、局所上昇アンブレラサンプリング法〔などの多くの新規手法が同様の原理に基いて開発された。この手法の基礎原理は、シミュレーション中に履歴に依存するポテンシャルエネルギー項を導入し、既にサンプリング済みの配座の再実現を阻害することで新しい配座の発見確率を上げるというものである。この手法はタブーサーチ法の連続版と見做すことができる。 == アルゴリズム == === 基礎ステップ === このアルゴリズムの基礎となるステップは、現在の分子配座を不利にするような小さな斥力的ポテンシャルエネルギーを加え、別の配座が見付かる確率を上げるというものである。このためには、関心のある配座変数を規定する自由度のサブセット を選ぶ必要がある。これは、通常は配座に関連する二面角とするが、原理的にはデカルト座標 の微分可能関数であればどんなものでも構わない。 このアルゴリズムでは、実際のポテンシャルエネルギーにバイアスエネルギーを印加し、以下のように変形させる。 局所上昇バイアス はシミュレーション時間 に依存し、シミュレーション開始時にはゼロとする () 。そして、以下の小さな斥力的関数を徐々に足し上げていく。 ここで、 はスケーリング定数、 は を満たす多次元斥力的関数である。 結果として得られるバイアスポテンシャルは、次のような総和で表わされる。 足し上げる斥力関数の数を減らすため、格子点を中心とする関数のみを用いる手法が広く用いられている。元々は、 として多変数ガウス関数が用いられていた。しかし、ガウス関数は全空間で非零であり、かつ格子点上のガウス関数を足し上げることでできてしまう人為的なエネルギー地形を避けるため、有限次の多項式関数を用いることがより良いと考えられるようになっている〔〔。 抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)』 ■ウィキペディアで「局所上昇法」の詳細全文を読む スポンサード リンク
|