• Bulletin of the Korean Chemical Society
• /
• 제23권3호
• /
• pp.441-446
• /
• 2002

We have presented the results of thermodynamic, structural and dynamic properties of liquid benzene, toluene, and p-xylene in canonical (NVT) ensemble at 293.15 K by molecular dynamics (MD) simulations. The molecular model adopted for these molecules is a combination of the rigid body treatment for the benzene ring and an atomistically detailed model for the methyl hydrogen atoms. The calculated pressures are too low in the NVT ensemble MD simulations. The various thermodynamic properties reflect that the intermolecular interactions become stronger as the number of methyl group attached into the benzene ring increases. The pronounced nearest neighbor peak in the center of mass g(r) of liquid benzene at 293.15 K, provides the interpretation that nearest neighbors tend to be perpendicular. Two self-diffusion coefficients of liquid benzene at 293.15 K calculated from MSD and VAC function are in excellent agreement with the experimental measures. The self-diffusion coefficients of liquid toluene also agree well with the experimental ones for toluene in benzene and for toluene in cyclohexane.

• Carbon letters
• /
• 제11권4호
• /
• pp.316-324
• /
• 2010

The effects of geometrical parameters on mechanical properties of graphite-vinylester nanocomposites and their constituents (matrix, reinforcement and interface) are studied using molecular dynamics (MD) simulations. Young's modulii of 1.3 TPa and 1.16 TPa are obtained for graphene layer and for graphite layers respectively. Interfacial shear strength resulting from the molecular dynamic (MD) simulations for graphene-vinylester is found to be 256 MPa compared to 126 MPa for graphitevinylester. MD simulations prove that exfoliation improves mechanical properties of graphite nanoplatelet vinylester nanocomposites. Also, the effects of bromination on the mechanical properties of vinylester and interfacial strength of the graphene.brominated vinylester nanocomposites are investigated. MD simulation revealed that, although there is minimal effect of bromination on mechanical properties of pure vinylester, bromination tends to enhance interfacial shear strength between graphite-brominated vinylester/graphene-brominated vinylester in a considerable magnitude.

• International Journal of Automotive Technology
• /
• 제5권3호
• /
• pp.215-219
• /
• 2004

A vehicle structure needs to be more precisely analyzed because of complexities and varieties. Structural fatigue which is generated by fluctuations of stresses during the service life of a mechanical system is the primary concern in the structural design for safety. A fatigue life is difficult to obtain in structural components during the service life of mechanical systems since the fluctuating stress contributes to fatigue. This study introduces new procedures to measure the lethargy coefficient and to predict the fatigue life of a mechanical structure by using molecular dynamic simulation. A lethargy coefficient is the total defect-estimating coefficient, which was obtained by using the results of a simple tensile test in this study. With this lethargy coefficient, fatigue life was estimated. The proposed method will be useful in predicting the fatigue life of a structurally-modified vehicle design. The effectiveness of the proposed method using lethargy coefficient measurement to predict the fatigue life of a structure was examined by applying this method to predict the fatigue life of SS41 steel, used extensively as material of vehicle structures. Two types of specimen such as pre-cracked plate and simple plate is discussed. equation of fatigue life using the lethargy coefficient and failure time, both obtained from a simple tensile test, will be useful in engineering. This measurement and prediction technology will be extended for use in analysis of any geometric shapes of modified automotive structures.

• 대한금속재료학회지
• /
• 제47권1호
• /
• pp.32-37
• /
• 2009

Molecular dynamic simulations on the structural evolution of the Si(001) vicinal surfaces, which are tilted with respect to [100] and [110] directions were performed by using the empirical Tersoff potential. Tersoff potential was implemented at LAMMPS code and confirmed to describe the properties of Si. When the steps are generated along [100] direction, symmetric dimer rows formed with respect to the step edges. On the other hand, when the steps are generated along [110] direction, alternating dimer rows form with respect to the step edges. The configurational differences between the two vicinal surfaces were discussed in terms of the surface diffusion and the possibility of preventing step bunching for the (001) vicinal surface tilted along [100] direction was suggested.

• 한국진공학회지
• /
• 제15권4호
• /
• pp.410-420
• /
• 2006

Epitaxial 성장에서 screening 및 steering 효과 등과 같은 증착과정 중 나타나는 dynamic effects를 kinetic Monte Carlo 시뮬레이션으로 고찰하였다. 증착원자와 토대 원자와의 상호 작용을 엄밀하게 고려하기 위해 이 시뮬레이션 프로그램에 molecular dynamics 시뮬레이션을 결합시켰다. 기울어진 각도로 증착 시킬 경우 표면의 형상은 1) 증착 각도가 기울어짐에 따라 거칠기가 증가한다는 것, 2) 비대칭적인 언덕의 형성된다는 것, 그리고 3) 언덕의 면방향에 따라 비대칭적인 기울기를 갖는다는 세가지 특징을 보았다. 증착 각도나 온도 의존성에 대한 시뮬레이션 결과는 기존의 실험 결과와 잘 일치하는 것을 확인하였다. 기울어진 각도로 증착했을 때 나타나는 이러한 결과는 steering과 screening 효과에 따른 초기 증착 밀도의 불균일에 기인함을 알 수 있었고, 여기서 steering 효과가 screening 효과 보다 주요한 역할을 하는 것을 보았다. 시뮬레이션 계산에서 확인된 새로운 결과는 기울기 선택 (slope selection)이 이루어졌어도 언덕의 각 면이 단일한 기울기로 형성되어 있지 않고 여러 종류의 facet가 섞인 형태이며, 따라서 기울기 선택이 바로 facet 선택을 의미하지는 않는다는 것이다.

• 폴리머
• /
• 제35권6호
• /
• pp.513-519
• /
• 2011

침투 가능한 구형 모델 유체의 충돌 특성을 고찰하고자 분자 동력학 방법을 이용한 전산 모사를 수행하였다. 이로부터 다양한 범위의 입자 충전 분율 ${\phi}$ 및 척력적 에너지 ${\varepsilon}^*$ 조건에 대한 충돌 빈도수, 평균 자유 행로, 강체형 반사 충돌각 및 연체형 침투 충돌각의 분포도, 유효 충전 분율 등을 계산하였다. 낮은 척력적 에너지 조건에서는 연체형 충돌이 주된 특성이나, 반면 높은 척력적 에너지 조건에서는 강체형 충돌이 주된 요인이었다. 매우 흥미롭게도, 전체 충돌 빈도수에 대한 연체형 충돌비(또는, 강체형 충돌비)는 정준 앙상블의 몬테카를로 전산 모사에서 받아들임 확률(또는, 되돌림 확률)로 표시되는 볼쯔만 인자와 직접적으로 관계되었다. 이와 같은 거동 입자들의 동적 충돌 특성들은 ${\varepsilon}^*{\geqq}3.0$ 및 ${\phi}{\geqq}0.7$ 범위의 높은 척력적 에너지 및 높은 충전 분율 조건에서 제한적이었으며, 이는 침투성 구형 모델에서 나타나는 클러스터 형성 구조를 함축하고 있다.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2003년도 추계학술대회논문집
• /
• pp.339-346
• /
• 2003

Molecular dynamic simulations of nano indentation on single-crystal silicon (100) surface were performed using diamond indentor. Silicon substrate and diamond indentor were modeled diamond structure with Tersoff potential model. Phase transformation of silicon, incipient plastic deformation, change of incident temperature distribution are investigated through the change of potential energy distribution, displacement-load diagram, the change of kinetic energy distribution and displacements of silicon atoms. Phase transformation is highly localized and consists of a high-density region surrounding the tip. Axial load linearly increased according to the indenting depth. Number of atoms with high kinetic energy increased at the interface between substrate and indentor tip.

• 설비공학논문집
• /
• 제19권12호
• /
• pp.883-888
• /
• 2007

As a starting point of investigating what molecular dynamic simulations can reveal about the nature of atomic level of heating and cooling process, argon described by the LJ potential is considered. Stepwise heating and cooling of constant rates are simulated in the NPT (constant number, pressure and temperature) ensemble. Hysteresis is found due to the superheating and supercooling. Drastic change of volume and energy is involved with phase change, but the melting point can not be obtained by simply observing the changes of these quantities. Since liquid and solid phases can co-exist at the same temperature, Gibbs free energy should be calculated to find the temperature where the Gibbs free energy of liquid is equal to that of the solid since the equilibrium state is the state of minimum Gibbs free energy. The obtained melting temperature, $T^*=0.685$, is close to that of the experiment with only 2% error.

• Carbon letters
• /
• 제16권3호
• /
• pp.183-191
• /
• 2015

Gas transport through graphene-derived membranes has gained much interest recently due to its promising potential in filtration and separation applications. In this work, we explore Kr-85 gas radionuclide sequestration from natural air in nanoporous graphene oxide membranes in which different sizes and geometries of pores were modeled on the graphene oxide sheet. This was done using atomistic simulations considering mean-squared displacement, diffusion coefficient, number of crossed species of gases through nanoporous graphene oxide, and flow through interlayer galleries. The results showed that the gas features have the densest adsorbed zone in nanoporous graphene oxide, compared with a graphene membrane, and that graphene oxide was more favorable than graphene for Kr separation. The aim of this paper is to show that for the well-defined pore size called P-7, it is possible to separate Kr-85 from a gas mixture containing Kr-85, O₂ and N₂. The results would benefit the oil industry among others.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2010년도 제39회 하계학술대회 초록집
• /
• pp.93-93
• /
• 2010

In this talk, I will introduce a reactive force field (ReaxFF) molecular dynamics (MD) simulation. In contrast to common MD simulations with empirical FFs, we can predict chemical reactions (bond breaking and formation) in large scale systems with the ReaxFF simulation where all of the ReaxFF parameters are from quantum mechanical calculations such as density functional theory to provide high accuracy. Accordingly, the ReaxFF simulation provides both accuracy of quantum mechanical calculations and description of large scale systems of atomistic simulations at the same time. Here, I will first discuss a theory in the ReaxFF including the differences from other empirical FFs, and then show several applications for studying chemical reactions of SiHx radicals on Si surfaces, which is an important issue in Si process.