• Current Applied Physics
• /
• v.18 no.11
• /
• pp.1313-1319
• /
• 2018

The coil-to-globule behavior of poly{${\gamma}$-2-[2-(2methoxyethoxy)ethoxy]ethoxy-3-caprolactone} (PMEEECL) as a ${\gamma}$-substituted poly (${\varepsilon}$-caprolactone) was investigated via atomistic molecular dynamics (MD) simulation. For this purpose, radius of gyration, end-to-end distance and radial distribution function of the chain in the presence of water were calculated. Consequently, the lower critical solution temperature (LCST) of PMEEECL chain at which the coil-to-globule transition takes place, was determined in each calculated parameter curve. The simulation results indicated that the LCST of PMEEECL was occurred at close to 320 K, which is in a good agreement with previous experimental results. Additionally, the appearance of sudden change in both Flory-Huggins interaction parameter (${\chi}$) and interaction energy between the PMEEECL chain and water molecules at about 320 K confirmed the calculated LCST result. The radial distribution function (RDF) results showed that the affinity of the PMEEECL side chain to water molecules is lower than its backbone.

• 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.

• Bulletin of the Korean Chemical Society
• /
• v.33 no.3
• /
• pp.901-905
• /
• 2012

In this study, we investigated the pure geometrical effect of porous materials in gas adsorption using the grand canonical Monte Carlo simulations of primitive gas-pore models with various pore geometries such as planar, cylindrical, and random pore geometries. Although the model does not possess atomistic level details of porous materials, our simulation results provided many insightful information in the effect of pore geometry on the adsorption behavior of gas molecules. First, the surface curvature of porous materials plays a significant role in the amount of adsorbed gas molecules: the concave surface such as in cylindrical pores induces more attraction between gas molecules and pore, which results in the enhanced gas adsorption. On the contrary, the convex surface of random pores gives the opposite effect. Second, this geometrical effect shows a nonmonotonic dependence on the gas-pore interaction strength and length. Third, as the external gas pressure is increased, the change in the gas adsorption due to pore geometry is reduced. Finally, the pore geometry also affects the collision dynamics of gas molecules. Since our model is based on primitive description of fluid molecules, our conclusion can be applied to any fluidic systems including reactant-electrode systems.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.15 no.4
• /
• pp.162-167
• /
• 2005

As the properties of porous materials during the drying process relate to the atomistic defects of heterogeneous materials such as dislocation, grain, grain boundary, pore, etc., the knowledge of nano-scale analysis is needed in order to accurately analyze the drying process for porous materials. In this study, the atomic behavior of porous materials Is statically predicted by using the molecular dynamics simulation and the nano-scale material properties are computed. The elastic modulus, thermal expansion coefficient, and volumetric heat capacity numerically found from the molecular dynamics simulation are compared with those of experiment and theory and proved the accuracy.

• Composites Research
• /
• v.24 no.3
• /
• pp.17-24
• /
• 2011

The EAM simulation of nanoindentation was performed to investigate misfit dislocation slip in the Ag/Cu. The film layer, whose thickness in the range of 2-5nm, was indented by a spherical indenter with the N$\'{o}$se-Hoover thermostat condition. The simulation shows that the indentation position relative to misfit dislocation (MFD) has the effect on the dislocation, glide up or cross slip, for Ag film layer thickness less than 4 nm. Elastic energy variation during MFDs slip was revealed to be a key factor for the softening of Ag/Cu. The critical film layer thickness was evaluated for each case of Ag/Cu according to the spline extrapolation technique.

• Interaction and multiscale mechanics
• /
• v.6 no.3
• /
• pp.301-316
• /
• 2013

The influence of indenter geometry on nanoindentation was studied using a static molecular dynamics simulation. Dislocation nucleation, dislocation locks, and dislocation movements during nanoindentation into Al (001) were studied. Spherical, rectangular, and Berkovich indenters were modeled to study the material behaviors and dislocation activities induced by their different shapes. We found that the elastic responses for the three cases agreed well with those predicted from elastic contact theory. Complicated stress fields were generated by the rectangular and Berkovich indenters, leading to a few uncommon nucleation and dislocation processes. The calculated mean critical resolved shear stresses for the Berkovich and rectangular indenters were lower than the theoretical strength. In the Berkovich indenter case, an amorphous region was observed directly below the indenter tip. In the rectangular indenter case, we observed that some dislocation loops nucleated on the plane. Furthermore, a prismatic loop originating from inside the material glided upward to create a mesa on the indenting surface. We observed an unusual softening phenomenon in the rectangular indenter case and proposed that heterogeneously nucleating dislocations are responsible for this.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.34-37
• /
• 2009

본 논문에서는 원자적 계산(atomistic calculation)을 이용한 나노박막의 평형상태(self-equilibrium state)에 대한 해석기법을 제시한다. 두께가 얇은 나노박막은 표면 응력(surface stress)에 의한 영향으로 원자간 거리가 벌크상태의 거리보다 작아진다. 두께가 얇은 나노박막에서의 원자 사이의 거리는 표면 응력과 탄성계수들의 표현식으로 계산이 가능하며, 본 논문에서는 {100}, {111}, {110} 표면을 가지는 나노박막의 평형상태의 해석을 위한 해석적 방법을 제시한다. 원자 사이의 거리를 계산하기 위해서는 보다 정확한 표면 응력의 계산방법이 필요하다. 본 연구에서는 나노박막의 평형상태에 대한 해석을 위해 surface relaxation model을 제시하고, 이 모델을 이용하여 표면응력(surface stress)과 표면강성계수(surface stiffness tensor)와 같은 surface parameter의 계산을 수행한다. 본 논문에서 제시된 surface relaxation model을 검증하기 위하여 분자동역학 전산모사(molecular dynamics simulation)의 수치 결과를 제시하고, 본 연구에서 계산한 equilibrium strain과 비교 검증한다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.4 s.175
• /
• pp.830-838
• /
• 2000

The paper addresses an application of molecular dynamics technique for fracture mechanics. Molecular dynamics simulation is an atomistic approach, while typical numerical methods such as finite element methods are macroscopic. Using the potential functions, which express the energy of a molecular system, a virtual specimen with molecules is set up and the trajectory of every molecule can be calculated by Newton's equation of motion. Several three-dimensional models with various types of cracks are considered. The stress intensity factors, the sizes of plastic zone as well as the dislocation emission are sought to be compared with the analytical solutions, which result in good agreement.

• Proceedings of the IEEK Conference
• /
• 2003.07b
• /
• pp.731-734
• /
• 2003

An atomistic process modeling, Kinetic Monte Carlo simulation, has the advantage of being both conceptually simple and extremely powerful. Instead of diffusion equations, it is based on the definitions of the interactions between individual atoms and defects. Those interactions can be derived either directly from molecular dynamics, first principles calculations, or from experiment. In this paper, as a simple illustration of the kinetic Monte Carlo we simulate defects (self-interstitials and vacancies) diffusion after ion implantation in Si crystalline.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1083-1088
• /
• 2003

Molecular dynamics simulations on the deformation behavior of single-walled carbon nanotube are performed. Formation energies of CNT's by interatomic potentials are computed and compared with ab initio results. Bending and axial compression are applied under lattice statics and NVT ensemble conditions. Specifically, we focus on the mechanism of kink formation in bending. The simulation results are comprehensively explained in the framework of atomistic energetics. The effects of temperature and chirality on the deformation of carbon nanotube are also studied.