• Journal of the KSME
• /
• v.44 no.3
• /
• pp.55-63
• /
• 2004

이 글에서는 마이크로와 나노스케일의 해석에 사용하는 수치모사 방법인 직접모사 몬테 카를로 (Direct Simulation Monte Carlo : DSMC)방법과 분자동역학(Molecular Dynamics: MD)과이 관계에 대하여 설명한다.

• Proceeding of EDISON Challenge
• /
• 2014.03a
• /
• pp.347-358
• /
• 2014

본 논문에서는 양자 역학적 분자 동역학(Quantum Mechanical/Molecular Mechanical-Molecular Dynamics, QM/MM-MD)을 통해 수용액에 녹아 있는 Potassium Thiocyanate의 dynamics를 연구했다. Umbrella sampling technique을 활용하여 association/dissociation에 해당하는 Free energy surface를 구했다. 두 개의 Free energy minimum이 녹아 있는 두 이온의 center of mass 사이의 거리가 $4{\AA}$일 때와 $5{\sim}6{\AA}$ 부근일 때 나타났으며 $4{\AA}$일 때 더 안정 했다. 본 논문에서는 $4{\AA}$일 때를 Contact Ion Pair(CIP) $6{\AA}$일 때를 Dissociation Ion Pair(DlP)라고 칭했다. 이 minimum들이 무엇인 지를 밝혀 내기 위해 추가 연구를 수행하였다. Free energy 상에서 가장 안정 할 때(CIP) solute인 Potassium thiocyanate의 구조를 살펴 봤더니 Potassium ion은 Thiocyanate ion의 Sulfur보다 Nitrogen side를 선호하였다. 그 원인을 알아보기 위해 salvation shell의 구조를 Radial distribution function을 통해 살펴 봤더니 물 분자가 Nitrogen보다 Sulfur와 더 강한 상호작용을 하고 있었다. 그로 인해 Potassium ion이 Nitrogen을 선호한단 결과가 나온 것이다. 한편, 두 번째 minimum은 물 분자가 Potassium 이온과 Thiocyanate 이온 사이에 flexible하게 bridging을 하는 구조였다. 또한 단순 양자 계산을 통해서도 비슷한 구조를 얻을 수 있었다. 그러나 QM 계산은 0K에서 수행하는 것이기 때문에 엔트로피 효과가 없는 계산이지만 본 연구는 온도 300K로 실제 용매와 가깝게 수행함으로써 고정되어 있는 구조가 아니라 엔트로피와 엔탈피가 균형적으로 존재하는 실제 용액 속에서의 구조를 처음으로 보여주는 것이다.

• Bulletin of the Korean Chemical Society
• /
• v.23 no.3
• /
• pp.447-453
• /
• 2002

In this paper we have presented the results of thermodynamic, structural, and dynamic properties of model systems for liquid benzene, toluene and p-xylene in an isobaric-isothermal (NpT) ensemble at 283.15, 303.15, 323.15, and 343.15 K using molecular dynamics (MD) simulation. This work is initiated to compensate for our previous canonical (NVT) ensemble MD simulations [Bull. Kor. Chem. Soc. 2001, 23, 441] for the same systems in which the calculated pressures were too low. The calculated pressures in the NpT ensemble MD simulations are close to 1 atm and the volume of each system increases with increasing temperature. The first and second peaks in the center of mass g(r) diminish gradually and the minima increase as usual for the three liquids as the temperature increases. The three peaks of the site-site gC-C(r) at 283.15 K support the perpendicular structure of nearest neighbors in liquid benzene. Two self-diffusion coefficients of liquid benzene via the Einstein equation and via the Green-Kubo relation are in excellent agreement with the experimental measures. The self-diffusion coefficients of liquid toluene and p-xylene are in accord with the trend that the self-diffusion coefficient decreases with increasing number of methyl group. The friction constants calculated from the force auto-correlation (FAC) function with the assumption that the fast random force correlation ends at time which the FAC has the first negative value give a correct qualitative trends: decrease with increase of temperature and increase with the number of methyl group. The friction constants calculated from the FAC's are always less than those obtained from the friction-diffusion relation which reflects that the random FAC decays slower than the total FAC as described by Kubo [Rep. Prog. Phys. 1966, 29, 255].

• Korean Journal of Mineralogy and Petrology
• /
• v.33 no.3
• /
• pp.143-152
• /
• 2020

Birnessite is a layered manganese oxide mineral with ~7 Å of d-spacing. Because of its high cation exchange capacity, birnessite greatly impacts the chemical compositions of ground water and fluids in sediment pores. Understanding the cation exchange mechanisms requires atomistic investigations of the crystal structures and coordination environments of hydrated cations in the interlayer. In this study, we conducted classical molecular dynamics (MD) simulations, an atomistic simulation method of computational mineralogy, for triclinic Na-birnessite and K-birnessite whose chemical formula are from previous experiments. We report our MD simulation results of the crystal structures, coordination environments of Na⁺ and K⁺, and the polytypes of birnessite and compare them with available experimental results. The simulation results well reproduced experimental lattice parameters and provided atomic level information for the interlayer cation and water molecule sites that are difficult to distinguish in X-ray experiments. We also report that the polytype of the Mn octahedral sheets is identical between Na- and K-birnessite, but the cation positions differ from each other, demonstrating a correlation between the coordination environment of the interlayer cations and the crystal lattice parameters. This study shows that MD simulations are very promising in elucidating ion exchange reactions of birnessite.

• Journal of the Mineralogical Society of Korea
• /
• v.30 no.4
• /
• pp.161-172
• /
• 2017

Clay minerals play a major role in the geochemical cycles of metals in the Critical Zone, the Earth surface-layer ranging from the groundwater bottom to the tree tops. Atomistic scale research of the very fine particles can help understand the fundamental mechanisms of the important geochemical processes and possibly apply to development of hybrid nanomaterials. Molecular dynamics (MD) simulations can provide atomistic level insights into the crystal structures of clay minerals and the chemical reactivity. Classical MD simulations use a force field which is a parameter set of interatomic pair potentials. The ClayFF force field has been widely used in the MD simulations of dioctahedral clay minerals as the force field was developed mainly based on dioctahedral phyllosilicates. The ClayFF is often used also for trioctahedral mineral simulations, but disagreement exits in selection of the interatomic potential parameters, particularly for Mg atom-types of the octahedral sheet. In this study, MD simulations were performed for trioctahedral clay minerals such as brucite, lizardite, and talc, to test how the two different Mg atom types (i.e., 'mgo' or 'mgh') affect the simulation results. The structural parameters such as lattice parameters and interatomic distances were relatively insensitive to the choice of the parameter, but the vibrational power spectra of hydroxyls were more sensitive to the choice of the parameter particularly for lizardite.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2011.02a
• /
• pp.21-21
• /
• 2011

A series of Quasi-Elastic Neutron Scattering (QENS) experiments helps us to understand the single-particle (hydrogen atom) dynamics of a globular protein and its hydration water and strong coupling between them. We also performed Molecular Dynamics (MD) simulations on a realistic model of the hydrated hen-egg Lysozyme powder having two proteins in the periodic box. We found the existence of a Fragile-to-Strong dynamic Crossover (FSC) phenomenon in hydration water around a protein occurring at TL=$225{\pm}5K$ by analyzing Intermediate Scattering Function (ISF). On lowering of the temperature toward FSC, the structure of hydration water makes a transition from predominantly the High Density Liquid (HDL) form, a more fluid state, to predominantly the Low Density Liquid (LDL) form, a less fluid state, derived from the existence of a liquid?liquid critical point at an elevated pressure. We showed experimentally and confirmed theoretically that this sudden switch in the mobility of the hydration water around a protein triggers the dynamic transition (so-called glass transition) of the protein, at a temperature TD=220 K. Mean Square Displacement (MSD) is the important factor to show that the FSC is the key to the strong coupling between a protein and its hydration water by suggesting TL${\fallingdotseq}$TD. MD simulations with TIP4P force field for water were performed to understand hydration level dependency of the FSC temperature. We added water molecules to increase hydration level of the protein hydration water, from 0.30, 0.45, 0.60 and 1.00 (1.00 is the bulk water). These confirm the existence of the FSC and the hydration level dependence of the FSC temperature: FSC temperature is decreased upon increasing hydration level. We compared the hydration water around Lysozyme, B-DNA and RNA. Similarity among those suggests that the FSC and this coupling be universal for globular proteins, biopolymers.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.1
• /
• pp.34-41
• /
• 2005

An experimental study of molecular motions on a liquid-vapor interface is limited due to micro-scale characteristics of a system with an angstrom or a nanometer size Therefore, in recent, many studies for micro-scale systems have been conducted by a computer simulation because it is free from experimental limitations. In this study, through the molecular dynamic (MD) method. molecular behavior was clarified on a liquid-vapor interface and a criterion to distinguish between liquid and vapor was suggested by a potential energy and the number of neighboring molecules. At an interface. the potential energy of a molecule was increased but the number of neighboring molecules was decreased when the molecule moved into a vapor region from a liquid region, and vice versa.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2002.10b
• /
• pp.59-60
• /
• 2002

In this digest, we briefly review our current molecular dynamics (MD) simulations that utilize both the reactive empirical bond order potential (REBO) and the adaptive intermolecular REBO (AIREBO) potential energy functions. The AIREBO potential includes intermolecular interactions, so that self·assembled monolayers, and liquids, can be modeled. We have examined the mechanical and tribological properties of model self assembled monolayers and amorphous carbon films. Self-assembled monolayers are modeled by covalently bonding hydrocarbon chains to diamond substrates. Because the REBO potentials can model chemical reactions, specific compression and sliding induced chemical reactions were identified.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.332-336
• /
• 2007

Temperature is an essential process variable in nanoimprint lithography(NIL) where the temperature varies between room temperature and above the glass transition temperature. To simulate NIL process, we employ both the Nose-Poincare method for temperature controlled molecular dynamics(MD) and force field for polymer material i.e. polymethyl methacrylate(PMMA), which is most widely selected as NIL resist. Nose-Poincare method, which convinces the conservation of Hamiltonian structure and time-reversal symmetry, overcomes the drawbacks inherent in the conventional methods such as Nose thermostat and Nose-Hoover thermostat. Thus, this method exhibits enhanced numerical stability even when the temperature fluctuation is large. To describe PMMA, we adopt the force field which account for bond stretch, bending, torsion, inversion, partial charge, and van der Waals energy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.8
• /
• pp.559-564
• /
• 2009

The variation in the shape of water droplet moving on atomistically smooth solid surface in the presence of a constant body force is simulated using molecular dynamics simulation. We investigated how the advancing and receding contact angle of the moving water droplet changes on a solid surface having various characteristic energies. From the MD simulation results, we obtained the density profile defined as the number of water molecules at a given position. Then, assuming the water droplet periphery to be a circle, we calculated the contact angles by using a nonlinear fitting of the half-density contour line. The present simulation clearly shows the different profile of the advancing and receding contact angle for these three different interaction potential between the water droplet and the solid surface.