分子動力学プログラムMXDORTO用可視化ソフトウェアMDVISの開発 [Published online J. Comput. Chem., Jpn.,, 14, 147-151, by J-STAGE]

[Published online Journal of Computer Chemistry, Japan14, 147-151, by J-STAGE]
<Title:> 分子動力学プログラムMXDORTO用可視化ソフトウェアMDVISの開発
<Author(s):> 大石 駿介, 大川 政志
<Corresponding author E-Mill:> mokawa(at)numazu-ct.ac.jp

<Abstract:> 酸化物のシミュレーションに有用である分子動力学(MD)シミュレーションプログラムMXDORTO用の結果を可視化するソフトワエア Visualizer for Molecular Dynamics Simulation (MDVIS)を開発した. MDVISは,メイン,データ,オペレーション,チャートの4つのフォームで構成されている.メインフォームには3次元コンピュータグラフィックス (3DCG)によりMDシミュレーションで記録されている原子座標を連続的に表示する. オペレーションフォームの設定により特定の場所のみを表示したり,周期境界条件を考慮した拡張ができる.さらに選択した原子の軌跡を表示させる機能や選択 した2原子間の原子間距離を計算する機能を備えた.3DCGの画面に同期している温度,圧力,格子定数,密度,内部エネルギーなどのデータを数値 で表示するデータフォームとシミュレーション時間全般での変化を表示させるチャートフォームを実装した.
<Keywords:> Molecular dynamics simulation, Visualization software
<URL:> https://www.jstage.jst.go.jp/article/jccj/14/4/14_2015-0010/_article/-char/ja/

Frenkel Defects and Interstitial Atoms in Periclase (MgO) at High Temperature by Molecular Dynamics Simulation [Published online J. Comput. Chem., Jpn.,, 14, 117-123, by J-STAGE]

[Published online Journal of Computer Chemistry, Japan14, 117-123, by J-STAGE]
<Title:> Frenkel Defects and Interstitial Atoms in Periclase (MgO) at High Temperature by Molecular Dynamics Simulation
<Author(s):> Yasuhiro UEDA, Akira MIYAKE
<Corresponding author E-Mill:> miya(at)kueps.kyoto-u.ac.jp

<Abstract:> An analysis program was newly developed to explore the vacancies and interstitial atoms and to determine atomic migration, i.e., atomic diffusion. We applied this program to the results calculated by the molecular dynamics (MD) simulation of periclase (MgO) in which Schottky defects (vacancies) were not initially introduced. Generation, migration and extinction of Frenkel defects and interstitial atoms for only magnesium ion were first observed at high temperature in this MD system and they strongly corresponded to the change of the mean square distance (MSD) of magnesium ion in MD system. On the other hand, we could not observe Frenkel defects and interstitial atoms for oxygen ion and MSD value of oxygen ion had almost constant value. Generation, migration and extinction of Frenkel defects and interstitial atoms cannot be ignored for the diffusion process at high temperature.
<Keywords:> Keyword Periclase (MgO), Molecular dynamics simulation, Frenkel defect, Interstitial atom, Diffusion
<URL:> https://www.jstage.jst.go.jp/article/jccj/14/4/14_2015-0014/_article/-char/ja/

Molecular Dynamics Simulation of the Behavior of Beryllium Diffusion in Corundum [Published online J. Comput. Chem., Jpn.,, 14, 111-116, by J-STAGE]

[Published online Journal of Computer Chemistry, Japan14, 111-116, by J-STAGE]
<Title:> Molecular Dynamics Simulation of the Behavior of Beryllium Diffusion in Corundum
<Author(s):> Jun KAWANO, Muneyuki OKUYAMA, Takeshi MIYATA
<Corresponding author E-Mill:> j-kawano(at)mail.sci.hokudai.ac.jp

<Abstract:> Molecular dynamics (MD) simulations are highly useful for analyzing atomic behavior during diffusion, especially in systems that are difficult to investigate experimentally. The focus of the present study was the diffusion behavior of Be in corundum, which was analyzed by MD calculations. First, we derived new potential parameter sets for O, Al, and Be. This parameter set was verified to well reproduce the structures and properties of corundum, bromellite, and chrysoberyl. Based on MD simulations of corundum containing Be as interstitial atoms, where the simulations were performed using the newly derived potential parameters, the diffusion coefficient was estimated to be approximately 10-7 cm2/s at around 2100 K. This is consistent with previously published experimental results, which confirms the validity of the MD simulation. The present calculations also reveal the detailed atomic movement, where Be atoms jump between Al sites and/or interstitial sites, and that the activation energy of this process is approximately 1.1 × 102 kJ/mol.
<Keywords:> Diffusion, Molecular dynamics simulation, Corundum, Beryllium, Gemology
<URL:> https://www.jstage.jst.go.jp/article/jccj/14/4/14_2015-0016/_article/-char/ja/

Interatomic Potential Model for Molecular Dynamics Simulation of Lithium Borate Melts/Glasses [Published online J. Comput. Chem., Jpn.,, 14, 139-146, by J-STAGE]

[Published online Journal of Computer Chemistry, Japan14, 139-146, by J-STAGE]
<Title:> Interatomic Potential Model for Molecular Dynamics Simulation of Lithium Borate Melts/Glasses
<Author(s):> Naoya SAWAGUCHI, Kakeru YAMAGUCHI, Makoto SASAKI, Katsuyuki KAWAMURA
<Corresponding author E-Mill:> nasawa(at)mmm.muroran-it.ac.jp

<Abstract:> An improved interatomic potential model was proposed for molecular dynamics simulations of lithium borate melt/glass systems. Charge of ion was reconsidered and a new composition dependent ionic charge model was suggested. A new three-body potential model controlling B-O-B angles was also proposed. The three-body term functioned to avoid square network ring consisted of B-O bonds, without preventing the change of boron coordination number between three and four. The edge-shared tetrahedra of four-coordination boron observed in the previous simulation were cleared by applying this three-body potential model.
<Keywords:> Molecular dynamics, Lithium borate, Glass, Three-body potential, Partial ionic charge model
<URL:> https://www.jstage.jst.go.jp/article/jccj/14/4/14_2015-0017/_article/-char/ja/

分子動力学法によるSiO2組成フォージャサイトの熱膨張変化の研究 [Published online J. Comput. Chem., Jpn.,, 14, 105-110, by J-STAGE]

[Published online Journal of Computer Chemistry, Japan14, 105-110, by J-STAGE]
<Title:> 分子動力学法によるSiO2組成フォージャサイトの熱膨張変化の研究
<Author(s):> 大川 政志, 堤 涼, 大石 駿介, 松本 泰誠, 山口 力
<Corresponding author E-Mill:> mokawa(at)numazu-ct.ac.jp

<Abstract:> マイクロ孔を有するゼオライトは,相転移や負の熱膨張など興味深い熱的挙動を示すことが知られている.様々な温度での三次元細孔の挙動を知ることは触媒反 応の点から重要である.そこで本研究では立方晶のフォージャサイト(FAU)の構造を有するSiO2組成ゼオライトを分子 動力学法により100 Kから1100 Kの範囲で計算を行い,格子定数および細孔に相当するスーパーケージのリング構造の変化を検討した.この温度領域において立方晶の構造は維持され格子定数 aは温度上昇とともに減少し,負の熱膨張を示した.二体相関関数からは第一近接のSiとSiの原子間距離は2種類ありそれらの温度に対する挙動が 異なり高温での負の熱膨張はSi-Si距離が減少することにより引き起こされていることが分かった.スーパーケージの12員環のリング径の平均値 は温度に対してほとんど変化しなかったが,温度上昇に伴い振幅は大きくなり600 K以上では10%以上変化することが分かった.
<Keywords:> molecular dynamics simulation, zeolite, faujasite, negative thermal expansion, ring structure
<URL:> https://www.jstage.jst.go.jp/article/jccj/14/4/14_2015-0020/_article/-char/ja/

Molecular Dynamics Simulation of Layerd Double Hydroxide (LDH) Mineral: Hydrotalcite and Water Systems [Published online J. Comput. Chem., Jpn.,, 14, 131-138, by J-STAGE]

[Published online Journal of Computer Chemistry, Japan14, 131-138, by J-STAGE]
<Title:> Molecular Dynamics Simulation of Layerd Double Hydroxide (LDH) Mineral: Hydrotalcite and Water Systems
<Author(s):> Katsuyuki KAWAMURA
<Corresponding author E-Mill:> kats(at)okayama-u.ac.jp

<Abstract:> Some clay minerals are used and expected as barrier materials for engineering and in nature. Hydrotalcite, Mg6Al2[(OH)16|CO3]4H2O and a kind of LDH (layered double hydroxides), is one of the most effective candidates for the anion adsorbents and the barrier. In this study, the behavior of hydrotalcite was investigated by means of the molecular dynamics method. Cl and I-hydrotalcite – water systems were simulated for various mineral/water ratios. The structure and dynamic properties are predicted. Water at the surface of hydrotalcite shows the electric double layer composed of Stern layer of one H2O molecular layer thickness and large self-diffusion coefficient of H2O and diffusion layer of 2.5 nm thickness at the interface.
<Keywords:> Keyword Clay mineral, Molecular dynamics simulation, Physical property, Structure, Mineral-water, Interface
<URL:> https://www.jstage.jst.go.jp/article/jccj/14/4/14_2015-0012/_article/-char/ja/

分子動力学プログラムMXDORTO用可視化ソフトウェアMDVISの開発 [Published online in advanced , by J-STAGE]

[Advanced Published online Journal of Computer Chemistry, Japan, by J-STAGE]
<Title:> 分子動力学プログラムMXDORTO用可視化ソフトウェアMDVISの開発
<Author(s):> 大石 駿介, 大川 政志
<Corresponding author E-Mill:> mokawa(at)numazu-ct.ac.jp

<Abstract:> 酸化物のシミュレーションに有用である分子動力学(MD)シミュレーションプログラムMXDORTO用の結果を可視化するソフトワエア Visualizer for Molecular Dynamics Simulation (MDVIS)を開発した. MDVISは,メイン,データ,オペレーション,チャートの4つのフォームで構成されている.メインフォームには3次元コンピュータグラフィックス (3DCG)によりMDシミュレーションで記録されている原子座標を連続的に表示する. オペレーションフォームの設定により特定の場所のみを表示したり,周期境界条件を考慮した拡張ができる.さらに選択した原子の軌跡を表示させる機能や選択 した2原子間の原子間距離を計算する機能を備えた.3DCGの画面に同期している温度,圧力,格子定数,密度,内部エネルギーなどのデータを数値 で表示するデータフォームとシミュレーション時間全般での変化を表示させるチャートフォームを実装した.
<Keywords:> Molecular dynamics simulation, Visualization software
<URL:> https://www.jstage.jst.go.jp/article/jccj/advpub/0/advpub_2015-0010/_article/-char/ja/

分子動力学法によるSiO2組成フォージャサイトの熱膨張変化の研究 [Published online in advanced , by J-STAGE]

[Advanced Published online Journal of Computer Chemistry, Japan, by J-STAGE]
<Title:> 分子動力学法によるSiO2組成フォージャサイトの熱膨張変化の研究
<Author(s):> 大川 政志, 堤 涼, 大石 駿介, 松本 泰誠, 山口 力
<Corresponding author E-Mill:> mokawa(at)numazu-ct.ac.jp

<Abstract:> マイクロ孔を有するゼオライトは,相転移や負の熱膨張など興味深い熱的挙動を示すことが知られている.様々な温度での三次元細孔の挙動を知ることは触媒反 応の点から重要である.そこで本研究では立方晶のフォージャサイト(FAU)の構造を有するSiO2組成ゼオライトを分子 動力学法により100 Kから1100 Kの範囲で計算を行い,格子定数および細孔に相当するスーパーケージのリング構造の変化を検討した.この温度領域において立方晶の構造は維持され格子定数 aは温度上昇とともに減少し,負の熱膨張を示した.二体相関関数からは第一近接のSiとSiの原子間距離は2種類ありそれらの温度に対する挙動が 異なり高温での負の熱膨張はSi-Si距離が減少することにより引き起こされていることが分かった.スーパーケージの12員環のリング径の平均値 は温度に対してほとんど変化しなかったが,温度上昇に伴い振幅は大きくなり600 K以上では10%以上変化することが分かった.
<Keywords:> molecular dynamics simulation, zeolite, faujasite, negative thermal expansion, ring structure
<URL:> https://www.jstage.jst.go.jp/article/jccj/advpub/0/advpub_2015-0020/_article/-char/ja/

粗視化分子動力学法におけるダブルネットワークゲルのモデリング手法の開発と破壊プロセスへの応用 [Published online J. Comput. Chem., Jpn.,, 14, 94-95, by J-STAGE]

[Published online Journal of Computer Chemistry, Japan14, 94-95, by J-STAGE]
<Title:> 粗視化分子動力学法におけるダブルネットワークゲルのモデリング手法の開発と破壊プロセスへの応用
<Author(s):> 斎藤 圭祐, 樋口 祐次, 尾澤 伸樹, 久保 百司
<Corresponding author E-Mill:> momoji(at)imr.tohoku.ac.jp

<Abstract:> We examine the mechanism of fracture processes in a double-network (DN) gel model by using a coarse-grained molecular dynamics simulation. Initially, we develop a modeling method for DN gel containing both slightly and highly cross-linked networks, and then stretch the DN gel model. During stretching, the highly cross-linked network begins to dissociate at a strain of 1.0, increasing the stress. At strains from 4.0 to 5.0, the slightly and highly cross-linked networks simultaneously dissociate and the stress decreases. Then, the dissociation of the highly cross-linked network stops and only the slightly cross-linked network dissociates at a strain of 12.0, while the stress remains almost the same. We reveal that characteristics of each type of network gradually appear in the DN gel. Next, we change the polymer chain length to reveal its influence on the mechanical properties of the gel. An increase in the length of the slightly cross-linked network chains improves the strength of the DN gel, whereas that of the highly cross-linked network chains does not affect its strength. An increase in the slightly cross-linked network chain length increases the number of entanglements, leading to the increase in strength.
<Keywords:> Polymer, Double-network gel, Fracture process, Mechanical strength, Coarse-grained molecular dynamics
<URL:> https://www.jstage.jst.go.jp/article/jccj/14/3/14_2015-0053/_article/-char/ja/

The Nature of Si-O-Si Bonding via Molecular Orbital Calculations [Published online in advanced , by J-STAGE]

[Advanced Published online Journal of Computer Chemistry, Japan, by J-STAGE]
<Title:> The Nature of Si-O-Si Bonding via Molecular Orbital Calculations
<Author(s):> Fumiya NORITAKE, Katsuyuki KAWAMURA
<Corresponding author E-Mill:> noritake(at)eqchem.s.u-tokyo.ac.jp

<Abstract:> The nature of Si-O bonding and Si-O-Si bridging is discussed using molecular orbital calculations. We found the equilibrium geometries for two pyrosilisic acid molecules (C2V and 60° torsion) using Moller-Plesset perturbation theory and 6-311G (d,p) split valence basis set. The bent configuration of the Si-O-Si angle in equilibrium geometries can be explained by the balance of Coulombic repulsion between SiO4 tetrahedra and the energy of lone pair orbitals belonging to bridging oxygen atom without concerning the contribution of d-p π-bonding from the results of natural bonding orbital analysis. The energy surfaces of two pyrosilisic acid molecules with varying Si-O length to the bridging oxygen and Si-O-Si angle were calculated and we found the relationship between Si-O length to the bridging and Si-O-Si bridging angle. The Si-O bonding strengthens with increasing Si-O-Si angle because of stabilization in energy of Si-O bonding orbital with decreasing the hybridization index λ in spλ orbital of bridging oxygen and increase of Coulombic interaction between Si and bridging oxygen atom.
<Keywords:> Keyword Siloxanes, Silicates, Molecular Orbital Calculation, Natural Bonding Orbital Analysis
<URL:> https://www.jstage.jst.go.jp/article/jccj/advpub/0/advpub_2015-0009/_article/-char/ja/