• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to project page

2017 Fiscal Year Research-status Report

一般化カットオフ法と合成系密度力学による多様な分子動力学シミュレーションの実現

Research Project

Project/Area Number 17K05143
Research InstitutionOsaka University

Principal Investigator

福田 育夫  大阪大学, たんぱく質研究所, 招へい准教授 (40643185)

Co-Investigator(Kenkyū-buntansha) 森次 圭  横浜市立大学, 生命医科学研究科, 特任准教授 (80599506)
Project Period (FY) 2017-04-01 – 2020-03-31
Keywords分子動力学 / 相互作用計算 / 力学系 / 運動方程式 / サンプリング / 数値計算 / 数値積分
Outline of Annual Research Achievements

分子動力学 (Molecular Dynamics; MD) シミュレーションは,生体高分子系,機能材料系等さまざまなマクロな系の分子レベルでの解析を可能にし,我々の生活に大きく影響を与え得る技術の一つになってきている.しかし,現状のMD計算法においては,シミュレーション条件について,技術的問題に起因した様々な制約が存在する.これを解消するために,具体的に,従来のカットオフ(相互作用の打ち切り)の考え方を一般化する方法,および対象とする物理系に他の任意の環境系を容易かつ統計力学的厳密にカップルさせる安定な数値解法を開発した.これにより,多様なシミュレーションを精密かつ安定に実現する方法の確立に資する.
一般化カットオフ法は,系の総電荷が零でなければならない,或いは系が周期的でなければならない,といった従来のさまざまな制約を解消するために開発する方法である.本年度は,基本方法である1次元積分表現の検討を行なった.また,この方法の方向性を吟味するために,関連手法であり,我々が既に開発したZero-multipole summation (ZM) 法についての詳細検証を行った.
物理系に他の任意の環境系をカップルさせる理論的可能性に関しては,既に我々が発明した合成系密度力学 (Double Density Dynamics; DDD) 法により基礎付けられている.但し,十分安定なシミュレーションを実行するための技術的課題があり,それを段階的に解消していく事が必要である.本年度は,サンプリングをより強力にして安定なシミュレーションを実現する技術を開発した.

Current Status of Research Progress
Current Status of Research Progress

2: Research has progressed on the whole more than it was originally planned.

Reason

一般化カットオフ法について,1次元積分表現の検討を行なった.これにより,予定していた連続近似を用いた相互作用関数を用いる方向性が正しい事が分かった.一方,その数学的表式を吟味していく中で,遠方の寄与を繰り込む新しい方向性にも気づいたので,その可能性についての検討も行った.
ZM法の吟味としては、境界条件に由来する誤差の定量的見積もりを行った.本来非周期的である系に,周期性を仮定した従来のエバルト法を用いて周期境界条件を課すと一定のアーティファクトが生じること,及びZM法では,そのような仮定が不要な分,誤差が小さくなることを,水中のペプチド系にて,種々の量の自由エネルギーを測ることで明らかにした.また,精度とともに,ZM法の計算速度性能の詳細な調査を行った.
DDD法では,対象となる物理系を取り囲む環境の温度を動的に変えつつ,各温度の実現確率を任意に定めることができ,この実現確率と,各温度が実現した際の物理系変数の実現する条件付き確率との積として,全系―物理系と温度(環境)系―の実現確率を表す.DDD法は全系のこの任意の確率密度に比例する不変測度を持つような常微分方程式で定義される運動方程式であり,これを解く事で物理状態の強化されたサンプリングを行う.本年度は,サンプリングをより強力にして安定なシミュレーションを実現する技術を開発した.具体的には,まず,近似分配関数の情報を逆温度の設定確率分布関数に取り入れる技術を開発した.次に,収束因子の導入を行った.これは物理量期待値の有限性を保証するために運動方程式に入れる因子であり,温度分布の両端での振る舞いを安定化させる効果を持つ.これらの方法をテストモデル系に適用した数値検証を経て,得られた逆温度定確率分布関数の精度が実際に大きく改善されることを,水中のペプチド系にて検証した.加えて,任意のパラメタに対応可能な安定な数値積分法を構築した.

Strategy for Future Research Activity

一般化カットオフ法に関して、当初想定していたものとは別の一般化の道筋にもたどり着いたので、その可能性について引き続き探っていく。それとともに,予定していた3次元積分表現についての理論的考察を進める.
本年度手開発したDDD法の安定化法は,当初予定していた方式(の一つ)と必ずしも一致しないが,その再考察及び多数のシミュレーション結果から推測された方法であり,有力だと考えている.なお,予定していたもう一方の方式は,上記の方式に矛盾するものではなく,寧ろ,さらなる安定化と一般的をもたらす方法と考えらえるため,詳細な検討を行いたい.また,DDD法において揺らがせる温度を,物理系に付随する任意パラメタにまで一般化する仕組みを検討する.

Causes of Carryover

(理由) 研究発表・打ち合わせのための旅費を優先させたため.
(使用計画)予定していたクライアントマシンを購入する.

  • Research Products

    (7 results)

All 2018 2017

All Journal Article (2 results) (of which Peer Reviewed: 2 results) Presentation (4 results) (of which Int'l Joint Research: 1 results,  Invited: 2 results) Book (1 results)

  • [Journal Article] Enhanced Sampling of Molecular Dynamics Simulations of a Polyalanine Octapeptide: Effects of the Periodic Boundary Conditions on Peptide Conformation2018

    • Author(s)
      Kasahara Kota、Sakuraba Shun、Fukuda Ikuo
    • Journal Title

      The Journal of Physical Chemistry B

      Volume: 122 Pages: 2495~2503

    • DOI

      DOI: 10.1021/acs.jpcb.7b10830

    • Peer Reviewed
  • [Journal Article] Performance Evaluation of the Zero-Multipole Summation Method in Modern Molecular Dynamics2018

    • Author(s)
      Shun Sakuraba、Ikuo Fukuda
    • Journal Title

      Journal of Computational Chemistry

      Volume: 印刷中 Pages: 印刷中

    • Peer Reviewed
  • [Presentation] 合成系密度力学の開発2018

    • Author(s)
      福田育夫,森次圭
    • Organizer
      日本物理学会
  • [Presentation] Non-Ewald method for accurately and efficiently calculating electrostatic interactions in molecular simulations2017

    • Author(s)
      Ikuo Fukuda, Narutoshi Kamiya, Kota Kasahara, Han Wang, Shun Sakuraba, Haruki Nakamura
    • Organizer
      “Conformational Ensembles from Experimental Data and Computer Simulations” meeting sponsored by the Biophysical Society
    • Int'l Joint Research
  • [Presentation] Non-Ewald methods for simply and precisely calculating electrostatic interactions in molecular simulations2017

    • Author(s)
      Ikuo Fukuda
    • Organizer
      JCUP VIII
    • Invited
  • [Presentation] Zero-multipole summation method: non-Ewald method for simply and precisely calculating electrostatic interactions in molecular simulations2017

    • Author(s)
      Ikuo Fukuda
    • Organizer
      Seminar at Science for Life Laboratory, Stockholm University
    • Invited
  • [Book] 長距離力をいかに精度よく効率的に見積もるか―零多重極子和法による静電相互作用計算―2017

    • Author(s)
      福田育夫
    • Total Pages
      7
    • Publisher
      日本物理学会誌

URL: 

Published: 2018-12-17  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi