Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
志賀 基之; Tuckerman, M. E.*
Journal of Physical Chemistry Letters, 9(21), p.6207 - 6214, 2018/11
パーセンタイル：100(Chemistry, Physical)反応経路の予測は、理論・計算化学の最も重要な課題の一つである。本論文では、最急降下法と最緩上昇法の組み合わせを用いて、化学反応の自由エネルギーランドマーク、すなわち最小および鞍点を自動探索するための新手法を提案する。この手法を実証するため、ベンゾシクロブテンの開環反応や水溶液中のSN2反応への応用例を紹介する。
鈴木 健人*; 宮崎 貴暉*; 高柳 敏幸*; 志賀 基之
Physical Chemistry Chemical Physics, 20(41), p.26489 - 26499, 2018/11
パーセンタイル：100(Chemistry, Physical)ヘリウムクラスターの直接光イオン化に続く短時間過程について、DIMモデルに基づく経路積分分子動力学(PIMD)およびリングポリマー分子動力学(RPMD)シミュレーションで調べた。PIMDシミュレーションによって実験で得られた非対称で幅広いイオン化スペクトルが再現され、その原因はHe原子のエネルギー準位の不均一性にあることがわかった。また、RPMDシミュレーションから、高励起状態にあるイオン化ヘリウムクラスターは、非断熱的電荷移動を通じて高速電子状態緩和した後、ゆっくりと構造緩和することがわかった。
志賀 基之
Reference Module in Chemistry, Molecular Sciences and Chemical Engineering, p.1 - 22, 2018/11
原子核の量子効果は水素などの軽元素を含む分子系において重要である。例えば、水素結合系における大きな零点振動およびプロトントンネリングは核量子効果に由来する。本レビュー論文では、虚時間の経路積分理論に基づいて核量子効果を考慮した経路積分分子動力学法について説明する。この方法では、複雑分子系の量子揺らぎにおいて数値的厳密解が得られる。また、その拡張であるセントロイドおよびリングポリマー分子動力学法は近似的量子動力学法である。第一原理電子状態計算と組み合わせることで、多電子・原子核系全体の量子シミュレーションが可能になる。
小座間 瑛記*; 安達 サディア*; 高柳 敏幸*; 志賀 基之
Chemistry; A European Journal, 24(48), p.12716 - 12721, 2018/08
パーセンタイル：100(Chemistry, Multidisciplinary)量子論的な経路積分分子動力学シミュレーションによって、ヘリウムクラスター中の三価アクチニウムカチオン(Ac
He
)の分子構造を計算した。安定なAc-He複合体が作られる低温(1-3K)において、ヘリウム原子の原子核量子効果は重要な役割を果たす。計算の結果、第一溶媒和圏を構成するヘリウム原子はD
対称をもって配置され、その配位数は18個にも及ぶことがわかった。このAc-He
複合体は、クラスターサイズ(
)が大きくなるにつれて強固になり、安定化する。なお、本研究は相対論的第一原理計算によりAc-He相互作用の正確な記述に基づくものである。
町田 昌彦; 加藤 孝一郎*; 志賀 基之
Journal of Chemical Physics, 148(10), p.102324_1 - 102324_11, 2018/03
被引用回数：5 パーセンタイル：8.9(Chemistry, Physical)水,重水,トリチウム水の液体状態における構造を第一原理経路積分分子動力学法を用いて研究した。用いた計算手法は、電子と原子核のすべてを量子力学的に取り扱うものであり、現時点で最も正確な計算手法の一つである。シミュレーションした結果は、実験とおおよそ一致しており、核量子効果により、同位体間の構造上の違いが明らかにされた。特に特筆すべき結果は、水分子内のOH共有結合中のHが隣接水分子との水素結合上の中点まで移動できる一方、重水素,トリチウムの順に中点まで行ける確率が減少するという傾向を得たことである。しかしながら、実験との完全な一致は未だ得られておらず、密度汎関数法の改良が必須であることが示唆された。また、最近の研究成果である密度汎関数に対するvan der Waals補正は液体状態における同位体効果を研究する上で必須であることが分かった。
君塚 肇*; 尾方 成信*; 志賀 基之
Physical Review B, 97(1), p.014102_1 - 014102_11, 2018/01
被引用回数：1 パーセンタイル：55.95(Materials Science, Multidisciplinary)パラジウム(Pd)中の水素(H)の高い拡散性のメカニズムを理解するため、有限歪み下における電子と核の両方の量子性を考慮した第一原理シミュレーションに基づいて、面心立方晶PdにおけるHの格子間拡散プロセスを研究した。その結果、八面体サイトの核量子効果のために、温度の低下とともに水素移動の活性化障壁が大幅に増加することが明らかになった。しかしながら、拡大した格子歪みのもとでは、四面体サイトがより安定化し、核量子効果はあまり顕著にならなかった。これは、歪みが増加するにつれて、拡散機構が量子的拡散から古典的拡散へ徐々に変化することを意味する。
関 悠佑*; 高柳 敏幸*; 志賀 基之
Physical Chemistry Chemical Physics, 19(21), p.13798 - 13806, 2017/06
被引用回数：5 パーセンタイル：33.94(Chemistry, Physical)500個のヘリウム原子からなる低温クラスターに埋め込まれたAg原子の光励起ダイナミクスを理解するため、リングポリマー分子動力学シミュレーションを行った。本研究では、
Pスピン軌道の電子波動関数の時間依存シュレーディンガー方程式とRPMD方程式を結合したエーレンフェスト平均場近似を用いた。このシミュレーションから、Ag原子は光励起後100psの時間でヘリウムクラスターから大部分が放出されることが分かった。放出したAg原子の平均速度は60
70m/sと見積もられ、実験結果と定性的に一致した。
樋山 みやび*; 志賀 基之; 古賀 信明*; 杉野 修*; 秋山 英文*; 野口 良史*
Physical Chemistry Chemical Physics, 19(15), p.10028 - 10035, 2017/04
被引用回数：2 パーセンタイル：68.71(Chemistry, Physical)本研究では、水溶液中のオキシルシフェリンアニオン異性体における特徴的な水和構造が電子吸収スペクトルに与える影響を検討した。第一原理分子動力学シミュレーションから得られた水和構造のカノニカルアンサンブルを用いて、室温でのケト-, エノール-およびエノラート型水性オキシルシフェリンアニオンの電子吸収スペクトルをQM/MM法で計算したところ、スペクトルシフトおよび広がりに関する実験結果をよく再現した。これは、オキシルシフェリンアニオンの分子振動とその水和構造の変動を考慮していない、静的な量子化学計算とフランク-コンドン近似に基づく従来法とは対照的な結果である。
Ruiz-Barragan, S.*; Ribas Ari
o, J.*; 志賀 基之
Physical Chemistry Chemical Physics, 18(47), p.32438 - 32447, 2016/12
被引用回数：3 パーセンタイル：71.73(Chemistry, Physical)反応場として熱水の使用は、グリーンケミストリー分野において非常に有望な技術である。これを活用するにあたり、熱水における反応機構を解明することが重要である。本研究では、熱水における2,5-ヘキサンジオールの脱水の反応機構について、ストリング法, メタダイナミクス法, 分子動力学法という三つの異なる第一原理シミュレーションを用いて研究を行った。その結果、この反応はプロトン化、結合交替と脱プロトン化を含む反応全体において、安定な中間体を形成することなく、連続的に進行することがわかった。この際、周囲の水の作る水素結合ネットワークは、反応の開始時と終了時における効率的なプロトンリレーを促す上で重要な役割を持っていることが明らかになった。この反応は36kcal/molのエネルギー障壁のあるSN2経路を取ることがわかり、実験で観測されている高い立体選択性と反応速度に一致する結果を得た。
野口 良史*; 樋山 みやび*; 志賀 基之; 杉野 修*; 秋山 英文*
Journal of Physical Chemistry B, 120(34), p.8776 - 8783, 2016/09
被引用回数：3 パーセンタイル：80.81(Chemistry, Physical)室温の水溶液中におけるオキシルシフェリンアニオンの安定性(ケト,エノール、およびエノラート異性体)を調べるため、ナノ秒スケールの大規模な第一原理分子動力学シミュレーションを行った。その結果、従来の量子化学計算とは対照的に、ケト型よりもエノール型の方がわずかに安定で、最近の幾つかの実験的結論に合致した。このシミュレーションから、ケト型が特に疎水的でと水との界面にキャビティが形成することや、疎水性が水分子数とともに減少する様子が明らかになった。このことから、安定な異性体は水和クラスターの大きさに依存するものと予想される。
箕島 裕介*; 関 悠佑*; 高柳 敏幸*; 志賀 基之
Chemical Physics, 472, p.1 - 8, 2016/06
被引用回数：1 パーセンタイル：91.53(Chemistry, Physical)放射線によって生じた電子がDNAに付着して生体損傷のもととなる際に分子構造がどのように変化するのか、基礎的問題として興味をもたれている。本研究では、グアニン-シトシン塩基対の電子付着過程について、経験的原子価結合モデルに基づく半古典リングポリマー分子動力学法を用いて理論的に解析した。その結果、塩基対の負イオンは、双極子束縛状態から原子価束縛状態へ短時間で変化し、その後塩基間のプロトン移動を引き起こすことがわかった。この過程における同位体置換効果や温度効果について詳しく調べた。プロトン移動の動的過程は、原子核の量子性を伴うものであることがわかった。
太田 幸宏; Ruiz-Barragan, S.*; 町田 昌彦; 志賀 基之
Chemical Physics Letters, 648, p.119 - 123, 2016/03
パーセンタイル：100(Chemistry, Physical)自動化学反応経路探索のソフトウエアGRRMと半経験的分子軌道計算法ソフトウエアMOPACのインターフェースツールの開発に成功した。このツールを使ってセルロースに対するイオンの吸着についてのテストを実施した所、計算効率は密度汎関数理論を用いたときより遥かに高速となり、複雑かつ巨大な分子に対する構造や反応経路の探索においては、極めて有効なツールとなることが分かった。
simulationsRuiz-Barragan, S.*; 石村 和也*; 志賀 基之
Chemical Physics Letters, 646, p.130 - 135, 2016/02
被引用回数：4 パーセンタイル：62.17(Chemistry, Physical)第一原理シミュレーションを階層的に並列化した新しい計算コードPIMD-SMASHを開発した。これによって、数十分子程度の水クラスター系について、経路積分分子動力学シミュレーションを効率的に実施できることを示した。また、開発されたコードを利用して、第一原理ストリング法によりシクロペンタジエンとブテノールのディールス・アルダー反応を調べた。その結果、水の水素結合によって反応のエネルギー障壁が下がることが見出され、この反応が水溶液中で加速するメカニズムを説明した。
本田 知大*; 箕島 裕介*; 横井 悠輝*; 高柳 敏幸*; 志賀 基之
Chemical Physics Letters, 625, p.174 - 178, 2015/04
被引用回数：3 パーセンタイル：77.78(Chemistry, Physical)密度汎関数理論および分子動力学法により、グアニン-シトシン(G-C)塩基対への電子付着の動力学的ふるまいについて研究を行った。G-Cの負イオンのポテンシャルエネルギーは、長距離補正密度汎関数計算の情報をもとに経験的原子価結合法でモデル化した。半古典的リングポリマー分子動力学法の結果、初期において双極子束縛状態にある負イオンは、0.1ピコ秒程度で原子価束縛状態へと容易に移行した後、10ピコ秒程度でプロトン移動が起こることがわかった。一方、古典分子動力学法の結果でも同じ傾向が見られたが、反応にかかる時間はずっと遅かった。このことから、低エネルギー電子によるDNA損傷過程のダイナミクスにおいて、原子核の量子効果が重要な役割を果たしていると考えられる。
Vela, S.*; Deumal, M.*; 志賀 基之; Novoa, J.*; Ribas-Arino, J.*
Chemical Science, 6(4), p.2371 - 2381, 2015/04
被引用回数：15 パーセンタイル：32.6(Chemistry, Multidisciplinary)これまで、分子磁性体における磁気的性質はX線結晶解析などで得られる一つの分子構造で決まるものと考えられてきた。しかし、第一原理計算に基づく計算化学的手法によりこれを検証したところ、一点の構造における磁性と、露に温度ゆらぎを取り入れた磁性の計算は、必ずしも一致しないことがわかった。特に、結晶中におけるラジカルが大振幅振動を起こすような、静的な描像に限界がある。そのような分子磁性体である1,3,5,トリチア-2,4,6トリアザペンタレニル(TTTA)系を例としてこれを実証する。
Vela, S.*; Mota, F.*; Deumal, M.*; 水津 理恵*; 珠玖 良昭*; 水野 麻人*; 阿波賀 邦夫*; 志賀 基之; Novoa, J.*; Ribas-Arino, J.*
Nature Communications (Internet), 5, p.4411_1 - 4411_9, 2014/07
近年、同じ熱力学的条件下で二つの安定状態を有する特性のもつ新しい有機ラジカル材料が注目を集めている。なかでも、中性のラジカルであるTTTAは室温下でヒステリシス性のスピン相転移を起こすジチアゾリル系の化合物として知られている。その機構を解明するため、第一原理分子動力学シミュレーションを行ったところ、重層構造をもったTTTA分子の単位がペアで交換する動力学が見られた。これが振動エントロピーの変化に大きな役割を果たし、相転移の引き金になっていることが明らかになった。
志賀 基之; Masia, M.*
Journal of Chemical Physics, 139(14), p.144103_1 - 144103_14, 2013/10
被引用回数：6 パーセンタイル：67.33(Chemistry, Physical)交換対称性理論に基づくマルチスケールシミュレーションの境界条件を利用して、一連のアルカリイオン水溶液のハイブリッドQM/MMシミュレーションを行った。本発表では、QM/MMシミュレーションにおける統計サンプリングの問題を、多重時間スケールアルゴリズムによって克服する方法を述べ、本方法を第一原理QM/MM経路積分分子動力学法に拡張することで、注目する部分系についてすべて量子的に扱った計算を実行した結果について言及する。特に、イオン-水間の動径分布関数やイオン周囲の水の配位角度など、水和構造を再現するのに、電子相関効果を適切に考慮した計算が必要であるという知見を得たことを指摘する。
藤崎 弘士*; 志賀 基之; 森次 圭*; 木寺 詔紀*
Journal of Chemical Physics, 139(5), p.054117_1 - 054117_9, 2013/08
被引用回数：7 パーセンタイル：62.07(Chemistry, Physical)分子シミュレーションでは、化学反応やタンパク質の折り畳みなど、構造変化がゆっくりで稀にしか起きない現象の扱い方が議論されている。これまで様々なアプローチが提案されてきたが、本論文では、シミュレーションを効率よく実行する方法の一つを提案する。これは、マルチスケール拡張サンプリング法を一般化し、粗視化された自由度と連動させたLangevin方程式を解くことで、Onsager-Machlup作用に基づいた反応経路計算の効率化を図るものである。この手法を朝倉-大沢相互作用のポリマー模型に応用し、計算効率については遷移経路サンプリング法など従来の方法と比較した結果を報告する。
志賀 基之; Masia, M.*
Journal of Chemical Physics, 139(4), p.044120_1 - 044120_8, 2013/07
被引用回数：14 パーセンタイル：38.82(Chemistry, Physical)本論文では、拡散系に対応した新しい多階層分子シミュレーション法を提案する。この手法では、溶液のように、階層間移動のある拡散粒子系において、その移動を妨げる緩やかな拘束力をかけて、粒子の混合を抑制する。このとき、粒子の交換対称性に則った拘束力を与えることで、その影響が系の熱力学的性質に全く及ばないようにすることができることを示す。簡単なモデルでこれを理論的に証明するとともに、応用として、セシウムイオン水溶液の多階層分子シミュレーションを、古典力場モデルの範囲で試験的に行うことで提案手法の検証を行った。
and Ag clusters小泉 亮人*; 立川 仁典*; 志賀 基之
Chemical Physics, 419, p.44 - 49, 2013/06
被引用回数：1 パーセンタイル：95.9(Chemistry, Physical)水和金属イオンに対して原子核の量子効果及び温度効果を明らかにするため、第一原理分子動力学法,経路積分分子動力学法及びリングポリマー分子動力学法を用いてAg(H
O)
(=1-4)とCu(HO)の分子構造と赤外スペクトルを調べた。その結果、量子効果及び温度効果による分子振動の非調和性は分子系や振動モードに依存していることが見いだされた。低振動モードは、複雑なモードカップリングの存在により、振動スペクトルは水和分子が多くなるにつれ幅広いバンドとなる。これに対してHOH変角振動モードのピークは、非調和性により系統的に赤方にシフトする。OH伸縮振動モードでも同様な赤方シフトが見られたが、シフトの大きさは系により依存する。したがって、上記のようなスペクトルの非調和性は、調和振動の一様なスケーリング則で表すことはできないことが結論づけられる。
