• 한국세라믹학회지
• /
• 제53권3호
• /
• pp.295-300
• /
• 2016

To investigate the thermal conductivity and the structural properties of naturally cooled excimer-laser annealed Si, molecular-dynamics (MD) simulations have been performed. The thermal conductivity of crystalline Si (c-Si) was measured by direct method at 1000 K. Steady-state heat flow was measured using a stationary temperature profile; significant deviations from Fourier's law were not observed. Reliable processes for measuring the thermal conductivity of c-Si were presented. A natural cooling process to admit heat flow from molten Si (l-Si) to c-Si was performed using an MD cell with a size of $48.9{\times}48.9{\times}97.8{\AA}^3$. During the cooling process, the temperature of the bottom $10{\AA}$ of the MD cell was controlled at 300 K. The results suggest that the natural cooling system described the static structural property of amorphous Si (a-Si) well.

• Macromolecular Research
• /
• 제12권2호
• /
• pp.178-188
• /
• 2004

We performed molecular dynamics (MD) simulations at 300 K on a series of poly(benzyl ether) (PBE) dendrimers having a different core functionalities. We used the rotational isomeric state Metropolis Monte Carlo (RMMC) method to construct the initial configuration in a periodic boundary cell (PBC) before the MD simulations were undertaken. To elucidate the effects that the structural features have on the chain dimension, the overall internal structure, and the morphology, we monitored the radii of gyration, R$\sub$g/ and the conformational changes during the simulations. The PBE dendrimers in a glassy state adopted less-extended structures when compared with the conformations obtained from the RMMC calculations. We found that R$\sub$g/ of the PBE dendrimer depends on the molecular weight, M, according to the relation, R$\sub$g/∼M$\^$0.22/. The radial distributions of the dendrimers were developed identically in the PBC, irrespective of the core functionality. A gradual decrease in radial density resulted from the fact that the terminal branch ends are distributed all over the molecule, except for the core region.

• Interaction and multiscale mechanics
• /
• 제1권2호
• /
• pp.303-315
• /
• 2008

In this paper, thermodynamical properties of crystalline silicon under strain are calculated using classical molecular dynamics (MD) simulations based on the Tersoff interatomic potential. The Helmholtz free energy of the silicon crystal under strain is calculated by using the ensemble method developed by Frenkel and Ladd (1984). To account for quantum corrections under strain in the classical MD simulations, we propose an approach where the quantum corrections to the internal energy and the Helmholtz free energy are obtained by using the corresponding energy deviation between the classical and quantum harmonic oscillators. We calculate the variation of thermodynamic properties with temperature and strain and compare them with results obtained by using the quasi-harmonic model in the reciprocal space.

• Bulletin of the Korean Chemical Society
• /
• 제22권8호
• /
• pp.847-849
• /
• 2001

The OSS2(Ojame-Shavitt-Singer 2)[L. Ojame et al., J. Chem. Phys. 109, 5547 (1998)] model as a dissociable water model is examined in order to study the dynamics of H+ in water. MD simulations for 216 water system, 215 water + H+ ion system, and 215 water + OH- ion system using the OSS2 model at 298.15 K with the use of Ewald summation are carried out. The calculated O-H radial distribution functions for these systems are essentially the same and are in very good agreement with that obtained by Ojame.

• 한국광물학회지
• /
• 제26권2호
• /
• pp.119-127
• /
• 2013

핵폐기물을 고화시키는 재료로 사용하는 붕규산염(borosilicate) 유리의 용해는 지층 처분장에 처리된 고준위 방사성 폐기물의 생태계 유출을 결정할 수 있는 중요한 화학반응이다. 습식 실험에서 유리의 용해속도(dissolution rate)는 유리 화학조성에 의해 크게 좌우되는 것이 관찰된다. 유리의 bulk 구조를 규명한 분광분석 실험에 의하면 유리의 화학조성과 분자수준(molecular-level) 구조(예: $SiO_4$ 사면체의 연결구조와 B 원소의 배위구조) 사이의 상관관계가 존재한다. 따라서 화학조성에 따른 유리 용해도의 차이는 조성에 따른 bulk 내부구조의 변화로 이해되어 왔다. 그런데 유리 표면은 수용액과 계면을 이루면서 용해 과정에서 가장 직접적으로 반응하는 부분이기 때문에, 화학조성에 따른 표면구조 변화에 대한 지식 또한 필요하다. 본 논문에서는 분자 동역학(molecular dynamics, MD) 시뮬레이션을 사용하여 4가지의 다른 화학조성을 가지는 소듐붕규산염 유리($xNa_2O{\cdot}B_2O_3{\cdot}ySiO_2$ 화학조성)에 대하여 bulk 구조와 실험으로 얻기 어려운 표면(surface) 구조를 연구하였다. MD 시뮬레이션은 유리 표면의 화학조성과 분자수준 구조가 bulk의 것과 매우 상이한 결과를 보여준다. 본 연구의 MD 시뮬레이션 결과는 화학조성에 따른 유리 용해도(특히 초기 용해과정)는 bulk 구조의 변화보다 유리 표면구조의 변화에 의해 크게 좌우될 수 있다는 표면구조에 대한 이해의 중요성을 역설한다.

• 한국광물학회지
• /
• 제29권4호
• /
• pp.209-220
• /
• 2016

점토광물은 대기부터 지하수에 이르는 크리티컬존(critical zone) 영역에서 금속과 탄소 순환을 결정짓는 역할을 한다. 계산광물학 연구방법 중에 하나인 분자동역학(molecular dynamics) 시뮬레이션은 지구물질을 원자단위로 계산하기 때문에, 점토광물의 물리화학적 현상들에 대해 원자수준의 자세한 정보를 제공할 수 있다. 이번 연구에서는 clayFF 힘 장(force field)을 사용한 분자동역학 시뮬레이션을 이팔면체 점토광물인 깁사이트(gibbsite, $Al(OH)_3$), 카올리나이트(kaolinite, $Al_2Si_2O_5(OH)_4$), 파이로필라이트(pyrophyllite, $Al_2Si_4O_{10}(OH)_2$)에 적용하여 300 K, 1기압조건에서 각 광물이 가지는 격자상수와 원자간 거리를 계산하고 실험결과와 비교하였다. 더불어 수산기의 방향성 및 수소결합의 양상 그리고 파워스펙트럼(power spectrum)을 추가적으로 계산하였다. 계산 결과, 격자상수는 기존의 실험결과와 약 0.1~3.7% 미만의 오차율을 보였으며, 원자간 거리는 실험결과와 약 5% 미만의 차이를 가졌다. 깁사이트나 카올리나이트의 팔면체층 표면에 존재하는 수산기가 가지는 신축진동파수(stretching vibrational wavenumber)는 실험값 보다 약 $200-300cm^{-1}$ 높게 계산되었지만, 팔면체층 표면에 존재하는 수산기들의 상대적 크기의 경향은 기존 실험 결과와 일치하였다. 팔면체층 표면의 수산기가 (001)면과 이루는 각도도 기존 실험결과와 상당부분 일치한 반면에 내부 수산기의 경우는 다소 차이를 보였다. ClayFF를 사용한 분자동역학 시뮬레이션 연구 방법은 이팔면체 점토광물 표면 내(층간) 금속이온 흡착에 대한 수산기의 역할을 규명하는데 유용한 연구방법이 될 수 있음을 시사한다.

• Bulletin of the Korean Chemical Society
• /
• 제18권5호
• /
• pp.501-509
• /
• 1997

In recent papers[Bull. Kor. Chem. Soc. 1996, 17, 735; ibid 1997, 18, 478] we reported results of molecular dynamics (MD) simulations for the thermodynamic, structural, and dynamic properties of liquid normal alkanes, from n-butane to n-heptadecane, using three different models. Two of the three classes of models are collapsed atomic models while the third class is an atomistically detailed model. In the present paper we present results of MD simulations for the corresponding properties of liquid branched-chain alkanes using the same models. The thermodynamic property reflects that the intermolecular interactions become weaker as the shape of the molecule tends to approach that of a sphere and the surface area decreases with branching. Not like observed in the straight-chain alkanes, the structural properties of model Ⅲ from the site-site radial distribution function, the distribution functions of the average end-to-end distance and the root-mean-squared radii of gyration are not much different from those of models Ⅰ and Ⅱ. The branching effect on the self diffusion of liquid alkanes is well predicted from our MD simulation results but not on the viscosity and thermal conductivity.

• Bulletin of the Korean Chemical Society
• /
• 제19권10호
• /
• pp.1068-1072
• /
• 1998

We present results of molecular dynamic (MD) simulations for the segmental motion of liquid n-butane as the base case for a consistent study for conformational transition from one rotational isomeric state to another in long chains of liquid n-alkanes. The behavior of the hazard plots for n-butane obtained from our MD simulations are compared with that for n-butane of Brownian dynamics study. The MD results for the conformational transition of n-butane by a Poisson process form the total first passage times are different from those from the separate t-g and g-t first passage times. This poor agreement is probably due to the failure of the detailed balance between the fractions of trans and gauche. The enhancement of the transitions t-g and g-t at short time regions are also discussed.

• Coupled systems mechanics
• /
• 제6권1호
• /
• pp.41-54
• /
• 2017

Shape-memory alloys (SMA) have interesting behaviors and important mechanical properties due to the solid-solid phase transformation. These phenomena are dominated by the evolution of microstructures. In recent years, the microstructures in SMAs have been studied extensively and modeled using molecular dynamics (MD) simulations. However, it remains difficult to identify the crystal variants in the simulation results, which consist of large numbers of atoms. In the present work, a method is developed to identify the austenite phase and the monoclinic martensite crystal variants in MD results. The transformation matrix of each lattice is calculated to determine the corresponding crystal variant. Evolution of the volume fraction of the crystal variants and the microstructure in Ni-Ti SMAs under thermal and mechanical boundary conditions are examined. The method is validated by comparing MD-simulated interface normals with theoretical solutions. In addition, the results show that, in certain cases, the interatomic potential used in the current study leads to inconsistent monoclinic lattices compared with crystallographic theory. Thus, a specific modification is applied and the applicability of the potential is discussed.

• Bulletin of the Korean Chemical Society
• /
• 제34권3호
• /
• pp.753-758
• /
• 2013

Molecular dynamics (MD) simulations have been used to study the diffusion behavior of small gas molecules ($CO_2$) in polyethylene terephthalate (PET)/polylactide (PLA) blends. The Flory-Huggins interaction parameters (${\chi}$) determined from the cohesive energy densities are smaller than the critical value of Flory-Huggins interaction parameters (${\chi}_{critical}$), and that indicates the good compatibility of PET/PLA blends. The diffusion coefficients of $CO_2$ are determined via MD simulations at 298 K. That the order of diffusion coefficients is correlated with the availably fractional free volume (FFV) of $CO_2$ in the PET/PLA blends means that the FFV plays a vital role in the diffusion behavior of $CO_2$ molecules in PET/PLA blends. The slopes of the log (MSD) as a function of log (t) are close to unity over the entire composition range of PET/PLA blends, which confirmes the feasibility of MD approach reaches the normal diffusion regime of $CO_2$ in PET/PLA blends.